Search

We present the properties of a massive, large, dusty, metal-rich, star-forming galaxy at z_spec=6.73. GOODSN-100182 was observed with JWST/NIRSpec as part of the AURORA survey, and is also covered by public multi-wavelength HST and JWST imaging. While the large mass of GOODSN-100182 (~10^10 M_sun) was indicated prior to JWST, NIRCam rest-optical imaging now reveals the presence of an extended disk (r_eff~1.5 kpc). In addition, the NIRSpec R~1000 spectrum of GOODSN-100182 includes the detection of a large suite of rest-optical nebular emission lines ranging in wavelength from [OII]3727 up to [NII]6583. The ratios of Balmer lines suggest significant dust attenuation (E(B-V)_gas=0.40+0.10/-0.09), consistent with the red rest-UV slope inferred for GOODSN-100182 (beta=-0.50+/-0.09). The star-formation rate based on dust-corrected H-alpha emission is log(SFR(H-alpha)/ M_sun/yr)=2.02+0.13/-0.14, well above the z~7 star-forming main sequence in terms of specific SFR. Strikingly, the ratio of [NII]6583/H-alpha emission suggests almost solar metallicity, as does the ratio ([OIII]5007/H-beta)/([NII]6583/H-alpha) and the detection of the faint [FeII]4360 emission feature, whereas the [OIII]5007/[OII]3727 ratio suggests roughly 50% solar metallicity. Overall, the excitation and ionization properties of GOODSN-100182 more closely resemble those of typical star-forming galaxies at z~2-3 rather than z~7. Based on public spectroscopy of the GOODS-N field, we find that GOODSN-100182 resides within a significant galaxy overdensity, and is accompanied by a spectroscopically-confirmed neighbor galaxy. GOODSN-100182 demonstrates the existence of mature, chemically-enriched galaxies within the first billion years of cosmic time, whose properties must be explained by galaxy formation models., Comment: 16 pages, 13 figures, submitted to ApJ

The H{\alpha} nebular emission line is an optimal tracer for recent star formation in galaxies. With the advent of JWST, this line has recently become observable at z>3 for the first time. We present a catalog of 1013 H{\alpha} emitters at 3.73 obtained based purely on spectroscopic data, robustly tracing galaxy star formation rates (SFRs) beyond the peak of the cosmic star formation history. We compare our results with theoretical predictions from three different simulations and find good agreement at z~4-6. The UV LFs of this spectroscopically-confirmed sample are in good agreement with pre-JWST measurements obtained with photometrically-selected objects. Finally, we derive SFR functions and integrate these to compute the evolution of the cosmic star-formation rate densities across z~4-6, finding values in good agreement with recent UV estimates from Lyman-break galaxies, which imply a continuous decrease in SFR density by a factor of 3x over z~4 to z~6. Our work shows the power of NIRCam grism observations to efficiently provide new tests for early galaxy formation models based on emission line statistics., Comment: 17 pages, 14 figures

We use JWST/NIRSpec observations from the Assembly of Ultradeep Rest-optical Observations Revealing Astrophysics (AURORA) survey to constrain the shape of the nebular attenuation curve of a star-forming galaxy at z=4.41, GOODSN-17940. We utilize 11 unblended HI recombination lines to derive the attenuation curve spanning optical to near-infrared wavelengths (3751-9550 \r{A}). We then leverage a high-S/N spectroscopic detection of the rest-frame ultraviolet continuum in combination with rest-UV photometric measurements to constrain the shape of the curve at ultraviolet wavelengths. While this UV constraint is predominantly based on stellar emission, the large measured equivalent widths of H$\alpha$ and H$\beta$ indicate that GOODSN-17940 is dominated by an extremely young stellar population <10 Myr in age such that the UV stellar continuum experiences the same attenuation as the nebular emission. The resulting combined nebular attenuation curve spans 1400-9550 \r{A} and has a shape that deviates significantly from commonly assumed dust curves in high-redshift studies. Relative to the Milky Way, SMC, and Calzetti curves, the new curve has a steeper slope at long wavelengths ($\lambda>5000$ \r{A}) while displaying a similar slope across blue-optical wavelengths ($\lambda=3750-5000$ \r{A}). In the ultraviolet, the new curve is shallower than the SMC and Calzetti curves and displays no significant 2175 \r{A} bump. This work demonstrates that the most commonly assumed dust curves are not appropriate for all high-redshift galaxies. These results highlight the ability to derive nebular attenuation curves for individual high-redshift sources with deep JWST/NIRSpec spectroscopy, thereby improving the accuracy of physical properties inferred from nebular emission lines., Comment: 21 pages, 6 figures, 1 table, submitted to ApJ

We present results on the emission-line properties of z=1.4-7.5 star-forming galaxies in the Assembly of Ultradeep Rest-optical Observations Revealing Astrophysics (AURORA) Cycle 1 JWST/NIRSpec program. Based on its depth, continuous wavelength coverage from 1--5 microns, and medium spectral resolution (R~1000), AURORA includes detections of a large suite of nebular emission lines spanning a broad range in rest wavelength. We investigate the locations of AURORA galaxies in multiple different emission-line diagrams, including traditional "BPT" diagrams of [OIII]/Hbeta vs. [NII]/Halpha, [SII]/Halpha, and [OI]/Halpha, and the "ionization-metallicity" diagram of [OIII]/[OII] (O32) vs. ([OIII]+[OII])/Hbeta (R23). We also consider a bluer rest-frame "ionization-metallicity" diagram introduced recently to characterize z>10 galaxies: [NeIII]/[OII] vs. ([NeIII]+[OII])/Hdelta; as well as longer-wavelength diagnostic diagrams extending into the rest-frame near-IR: [OIII]/Hbeta vs. [SIII]/[SII] (S32); and HeI/Pagamma and [SIII]/Pagamma vs. [FeII]/Pabeta. With a significant boost in signal-to-noise and large, representative samples of individual galaxy detections, the AURORA emission-line diagrams presented here definitively confirm a physical picture in which chemically-young, alpha-enhanced, massive stars photoionize the ISM in distant galaxies with a harder ionizing spectrum at fixed nebular metallicity than in their z~0 counterparts. We also uncover previously unseen evolution prior to z~2 in the [OIII]/Hbeta vs. [NII]/Halpha diagram, which motivates deep NIRSpec observations at even higher redshift. Finally, we present the first statistical sample of rest-frame near-IR emission-line diagnostics in star-forming galaxies at high redshift. In order to truly interpret rest-frame near-IR line ratios including [FeII], we must obtain better constraints on dust depletion in the high-redshift ISM., Comment: 17 pages, 8 figures, submitted to ApJ

Combining the public JWST/NIRCam imaging programs CEERS, PRIMER and JADES, spanning a total area of $\sim500\,{\rm arcmin}^2$, we obtain a sample of $>$30,000 galaxies at $z_{\rm phot}\sim4-9$ that allows us to perform a complete, rest-optical selected census of the galaxy population at $z>3$. Comparing the stellar mass $M_*$ and the UV-slope $\beta$ distributions between JWST- and HST-selected samples, we generally find very good agreement and no significant biases. Nevertheless, JWST enables us to probe a new population of UV-red galaxies that was missing from previous HST-based Lyman Break Galaxy (LBG) samples. We measure galaxy stellar mass functions (SMFs) at $z\sim4-9$ down to limiting masses of $10^{7.5}-10^{8.5}\,{\rm M_\odot}$, finding steep low mass slopes over the entire redshift range, reaching values of $\alpha\approx-2$ at $z\gtrsim6$. At the high-mass end, UV-red galaxies dominate at least out to $z\sim6$. The implied redshift evolution of the SMF suggests a rapid build-up of massive dust-obscured or quiescent galaxies from $z\sim6$ to $z\sim4$ as well as an enhanced efficiency of star formation towards earlier times ($z\gtrsim6$). Finally, we show that the galaxy mass density grows by a factor $\sim20\times$ from $z\sim9$ to $z\sim4$. Our results emphasize the importance of rest-frame optically-selected samples in inferring accurate distributions of physical properties and studying the mass build-up of galaxies in the first 1.5 Gyr of cosmic history., Comment: 31 pages, 17 figures, published in MNRAS

GN-z11 is the highest redshift galaxy spectroscopically confirmed with the Hubble Space Telescope (HST). Previous measurements of the effective radius of GN-z11 utilized galfit, which is not optimized to measure structural parameters for such a faint, distant object. Using a new software program called forcepho on HST data for the first time, we derive a size from images in the F160W band obtained both from the complete CANDELS survey and additional midcycle observations in order to contribute to the knowledge base on the size evolution, size-luminosity, and size-mass relation of early galaxies. We find a half-light radius mean of 0''.036 \(\pm\) 0''.006 corresponding to a physical size of 0.15 \(\pm\) 0.025 kpc. This size, smaller than the point spread function, is dramatically smaller than previous estimates with shallower HST data using galfit but consistent with recent measurements using forcepho on new JWST data arXiv:2302.07234. Such a small size, combined with the JWST/NIRSpec spectroscopic observations arXiv:2305.12492, suggests that GN-z11's high luminosity is dominated by an AGN., Comment: 4 pages, 1 figure

We present the H{\alpha} luminosity function (LF) derived from a large sample of Lyman break galaxies at z {\sim} 4.5 over the GOODS-South and North fields. This study makes use of the new, full-depth Spitzer/IRAC [3.6] and [4.5] imaging from the GOODS Re-ionization Era wide-Area Treasury from the Spitzer program. The H{\alpha} flux is derived from the offset between the continuum flux estimated from the best-fit spectral energy distribution, and the observed photometry in IRAC [3.6]. From these measurements, we build the H{\alpha} LF and study its evolution providing the best constraints of this property at high redshift, where spectroscopy of H{\alpha} is not yet available. Schechter parameterizations of the H{\alpha} LF show a decreasing evolution of {\Phi^\star} with redshift, increasing evolution in L{^\star}, and no significant evolution in the faint-end slope at high z. We find that star formation rates (SFRs) derived from H{\alpha} are higher than those derived from the rest-frame UV for low SFR galaxies but the opposite happens for the highest SFRs. This can be explained by lower mass galaxies (also lower SFR) having, on average, rising star formation histories (SFHs), while at the highest masses the SFHs may be declining. The SFR function is steeper, and because of the excess SFR(H{\alpha}) compared to SFR(UV) at low SFRs, the SFR density estimated from H{\alpha} is higher than the previous estimates based on UV luminosities., Comment: Astrophysical Journal, Accepted, 17 pages, 16 figures

We present the first constraints on the prevalence of z>10 galaxies in the Hubble Ultra Deep Field (HUDF) leveraging new NIRCam observations from JEMS (JWST Extragalactic Medium-band Survey). These NIRCam observations probe redward of 1.6$\mu$m, beyond the wavelength limit of HST, allowing us to search for galaxies to z>10. These observations indicate that the highest redshift candidate identified in the HUDF09 data with HST, UDFj-39546284, has a redshift of z>11.5, as had been suggested in analyses of the HUDF12/XDF data. This has now been confirmed with JWST NIRSpec. This source is thus the most distant galaxy discovered by HST in its >30 years of operation. Additionally, we identify nine other z~8-13 candidate galaxies over the HUDF, two of which are new discoveries that appear to lie at z~11-12. We use these results to characterize the evolution of the UV luminosity function (LF) from z~15 to z~8.7. While our LF results at z~8.7 and z~10.5 are consistent with previous findings over the HUDF, our new LF estimates at z~12.6 are higher than other results in the literature, potentially pointing to a milder evolution in the UV luminosity density from z~12.6. We emphasize that our LF results are uncertain given the small number of z~12.6 sources and limited volume probed. The new NIRCam data also indicate that the faint z~8-13 galaxies in the HUDF/XDF show blue UV-continuum slopes beta~-2.7, high specific star formation rates ~24.5 Gyr**-1, and high EW (~1300A) [OIII]+H$\beta$ emission, with two z~8.5 sources showing [OIII]+H$\beta$ EWs of ~2300 A., Comment: 21 pages, 14 figures, 7 tables, resubmitted to MNRAS, after replying to the referee report

$JWST$'s first glimpse of the $z>10$ Universe has yielded a surprising abundance of luminous galaxy candidates. Here we present the most extreme of these systems: CEERS-1749. Based on $0.6-5\mu$m photometry, this strikingly luminous ($\approx$26 mag) galaxy appears to lie at $z\approx17$. This would make it an $M_{\rm{UV}}\approx-22$, $M_{\rm{\star}}\approx5\times10^{9}M_{\rm{\odot}}$ system that formed a mere $\sim220$ Myrs after the Big Bang. The implied number density of this galaxy and its analogues challenges virtually every early galaxy evolution model that assumes $\Lambda$CDM cosmology. However, there is strong environmental evidence supporting a secondary redshift solution of $z\approx5$: all three of the galaxy's nearest neighbors at $<2.5$" have photometric redshifts of $z\approx5$. Further, we show that CEERS-1749 may lie in a $z\approx5$ protocluster that is $\gtrsim5\times$ overdense compared to the field. Intense line emission at $z\approx5$ from a quiescent galaxy harboring ionized gas, or from a dusty starburst, may provide satisfactory explanations for CEERS-1749's photometry. The emission lines at $z\approx5$ conspire to boost the $>2\mu$m photometry, producing an apparent blue slope as well as a strong break in the SED. Such a perfectly disguised contaminant is possible only in a narrow redshift window ($\Delta z\lesssim0.1$), implying that the permitted volume for such interlopers may not be a major concern for $z>10$ searches, particularly when medium-bands are deployed. If CEERS-1749 is confirmed to lie at $z\approx5$, it will be the highest-redshift quiescent galaxy, or one of the lowest mass dusty galaxies of the early Universe detected to-date. Both redshift solutions of this intriguing galaxy hold the potential to challenge existing models of early galaxy evolution, making spectroscopic follow-up of this source critical., Comment: Submitted to ApJL. Fig. 1 shows the "superposed" state, Fig. 2 presents the potential z~5 protocluster, Fig. 3 shows some atypical z~5 SEDs, Fig. 4 studies the morphology, Fig. 5 summarizes the extraordinary situation at z~17, and Fig. 6 argues for why such interlopers may not be a major issue for z>10 searches. Comments warmly welcomed

We present the deepest constraints yet on the median rest-UV+optical SED of $z\sim10$ galaxies, prior to JWST science operations. We constructed stacks based on four robust $J_{125}$-dropouts, previously identified across the GOODS fields. We used archival HST/WFC3 data and the full depth Spitzer/IRAC mosaics from the GREATS program, the deepest coverage at $\sim3-5\mu$m to date. The most remarkable feature of the SED is a blue IRAC $[3.6]-[4.5]=-0.18\pm0.25$ mag color. We also find a nearly flat $H_{160}-[3.6]=0.07\pm0.22$ mag color, corresponding to a UV slope $\beta= -1.92\pm0.25$. This is consistent with previous studies, and indicative of minimal dust absorption. The observed blue IRAC color and SED fitting suggest that $z\sim10$ galaxies have very young (few $\times10$ Myr) stellar populations, with $80\%$ of stars being formed in the last $\lesssim 160$ Myr ($2\sigma$). While an exciting result, the uncertainties on the SED are too large to allow us to place strong constraints on the presence of a nebular continuum in $z\sim10$ galaxies (as might be suggested by the blue $[3.6]-[4.5] < 0$ mag color). The resulting sSFR is consistent with the specific accretion rate of dark matter halos, indicative of a star-formation efficiency showing quite limited evolution at such early epochs., Comment: Submitted to ApJ. 9 pages, 3 figures. Comments welcome

We measure, for the first time, the median equivalent width (EW) of H$\alpha$+[N II] in star-forming galaxies at $z\sim8$. Our estimate leverages the unique photometric depth of the Spitzer/IRAC $5.8\mu$m-band mosaics (probing $\approx 5500 - 7100$ A at $z\sim8$) of the GOODS Reionization Era Wide Area Treasury from Spitzer (GREATS) program. We median stacked the stamps of $102$ Lyman-break galaxies in the $3.6, 4.5, 5.8$ and $8.0\mu$m bands, after carefully removing potential contamination from neighbouring sources. We infer an extreme rest-frame EW$_0$(H$\alpha$+[N II])$=2328^{+1326}_{-1127}$ A from the measured red $[3.6]-[5.8]=0.82\pm0.27$ mag, consistent with young ($\lesssim10^7$ yr) average stellar population ages at $z\sim8$. This implies an ionizing photon production efficiency of $\log(\xi_{\mathrm{ion},0}/\mathrm{erg\ Hz}^{-1})=25.97^{+0.18}_{-0.28}$. Such a high value for photo production, similar to the highest values found at $z\lesssim4$, indicates that only modest escape fractions $f_\mathrm{esc}\lesssim0.3$ (at $2\sigma$) are sufficient for galaxies brighter than $M_\mathrm{UV}<-18$ mag to re-ionize the neutral Hydrogen at $z\sim8$. This requirement is relaxed even more to $f_\mathrm{esc}\le 0.1$ when considering galaxies brighter than $M_\mathrm{UV}\approx -13$ mag, consistent with recent luminosity functions and as typically assumed in studies addressing re-ionization. These exceptional results clearly indicate that galaxies can be the dominant source of reionizing photons, and provide us with an exciting glimpse into what we might soon learn about the early universe, and particularly about the Reionization Epoch, from forthcoming JWST/MIRI and NIRCam programs., Comment: Submitted to ApJ. 12 pages, 5 figures (main) + appendix. Comments welcome

We present the deepest Spitzer/IRAC $3.6$, $4.5$, $5.8$ and $8.0\mu$m wide-area mosaics yet over the GOODS-N and GOODS-S fields as part of the GOODS Re-ionization Era wide-Area Treasury from Spitzer (GREATS) project. We reduced and mosaicked in a self-consistent way observations taken by the 11 different Spitzer/IRAC programs over the two GOODS fields from 12 years of Spitzer cryogenic and warm mission data. The cumulative depth in the $3.6\mu$m and $4.5\mu$m bands amounts to $\sim 4260$ hr, $\sim 1220$ hr of which are new very deep observations from the GREATS program itself. In the deepest area, the full-depth mosaics reach $\gtrsim200$ hr over an area of $\sim100$ arcmin$^2$, corresponding to a sensitivity of $\sim29$ AB magnitude at $3.6\mu$m ($1\sigma$ for point sources). Archival cryogenic $5.8\mu$m and $8.0\mu$m band data (a cumulative 976 hr) are also included in the release. The mosaics are projected onto the tangential plane of CANDELS/GOODS at a $0.3''$ pixel$^{-1}$ scale. This paper describes the methodology enabling, and the characteristics of, the public release of the mosaic science images, the corresponding coverage maps in the four IRAC bands, and the empirical Point-Spread Functions (PSFs). These PSFs enable mitigation of the source blending effects by taking into account the complex position-dependent variation in the IRAC images. The GREATS data products are in the Infrared Science Archive (IRSA). We also release the deblended $3.6$-to-$8.0\mu$m photometry for $9192$ Lyman-Break galaxies at $z\sim3.5-10$. GREATS will be the deepest mid-infrared imaging until JWST and, as such, constitutes a major resource for characterizing early galaxy assembly., Comment: Accepted for publication on ApJS. Mosaics, PSFs and catalogs are publicly available at https://home.strw.leidenuniv.nl/~stefanon/GREATS/

The highly neutral inter-galactic medium (IGM) during the Epoch of Reionization (EoR) is expected to suppress Ly$\alpha$ emission with damping-wing absorption, causing nearly no Ly$\alpha$ detection from star-forming galaxies at $z{\sim}8$. However, spectroscopic observations of the 4 brightest galaxies (${\rm H}_{160}{\sim}25$ mag) at these redshifts do reveal prominent Ly$\alpha$ line, suggesting locally ionised IGM. In this paper, we explore the Ly$\alpha$ IGM transmission and environment of bright galaxies during the EoR using the Meraxes semi-analytic model. We find brighter galaxies to be less affected by damping-wing absorption as they are effective at ionizing surrounding neutral hydrogen. Specifically, the brightest sources (${\rm H}_{160}{\lesssim}25.5$ mag) lie in the largest ionized regions in our simulation, and have low attenuation of their Ly$\alpha$ from the IGM (optical depth ${<}1$). Fainter galaxies (25.5 mag${<}{\rm H}_{160}{<}27.5$ mag) have transmission that depends on UV luminosity, leading to a lower incidence of Ly$\alpha$ detection at fainter magnitudes. This luminosity-dependent attenuation explains why Ly$\alpha$ has only been observed in the brightest galaxies at $z{\sim}8$. Follow-up observations have revealed counterparts in the vicinity of these confirmed $z{\sim}8$ Ly$\alpha$ emitters. The environments of our modelled analogues agree with these observations in the number of nearby galaxies, which is a good indicator of whether Ly$\alpha$ can be detected among fainter galaxies. At the current observational limit, galaxies with ${\ge}2$--5 neighbours within $2'{\times}2'$ are ${\sim}2$--3 times more likely to show Ly$\alpha$ emission. JWST will discover an order of magnitude more neighbours, revealing ${\gtrsim}50$ galaxies in the largest ionizing bubbles and facilitating direct study of reionization morphology., Comment: 10 pages, 3 figures, submitted to MNRAS, comments are welcome

We make use of ALMA continuum observations of $15$ luminous Lyman-break galaxies at $z$$\sim$$7$$-$$8$ to probe their dust-obscured star-formation. These observations are sensitive enough to probe to obscured SFRs of $20$ $M_{\odot}$$/$$yr$ ($3\sigma$). Six of the targeted galaxies show significant ($\geq$$3$$\sigma$) dust continuum detections, more than doubling the number of known dust-detected galaxies at $z$$>$$6.5$. Their IR luminosities range from $2.7$$\times$$10^{11}$ $L_{\odot}$ to $1.1$$\times$$10^{12}$ $L_{\odot}$, equivalent to obscured SFRs of $20$ to $105$ $M_{\odot}$$/$$yr$. We use our results to quantify the correlation of the infrared excess IRX on the UV-continuum slope $\beta_{UV}$ and stellar mass. Our results are most consistent with an SMC attenuation curve for intrinsic $UV$-slopes $\beta_{UV,intr}$ of $-2.63$ and most consistent with an attenuation curve in-between SMC and Calzetti for $\beta_{UV,intr}$ slopes of $-2.23$, assuming a dust temperature $T_d$ of $50$ K. Our fiducial IRX-stellar mass results at $z$$\sim$$7$$-$$8$ are consistent with marginal evolution from $z$$\sim$$0$. We then show how both results depend on $T_d$. For our six dust-detected sources, we estimate their dust masses and find that they are consistent with dust production from SNe if the dust destruction is low ($<$$90$%). Finally we determine the contribution of dust-obscured star formation to the star formation rate density for $UV$ luminous ($<$$-$$21.5$ mag: $\gtrsim$$1.7$$L_{UV} ^*$) $z$$\sim$$7$$-$$8$ galaxies, finding that the total SFR density at $z$$\sim$$7$ and $z$$\sim$$8$ from bright galaxies is $0.18_{-0.10}^{+0.08}$ dex and $0.20_{-0.09}^{+0.05}$ dex higher, respectively, i.e. $\sim$$\frac{1}{3}$ of the star formation in $\gtrsim$$1.7$$L_{UV} ^*$ galaxies at $z$$\sim$$7$$-$$8$ is obscured by dust., Comment: 21 pages, 12 figures, 2 tables. Submitted to ApJ

We present new stellar mass functions at $z\sim6$, $z\sim7$, $z\sim8$, $z\sim9$ and, for the first time, $z\sim10$, constructed from $\sim800$ Lyman-Break galaxies previously identified over the XDF/UDF, parallels and the five CANDELS fields. Our study is distinctive due to (1) the much deeper ($\sim200$ hour) wide-area Spitzer/IRAC imaging at $3.6\mu$m and $4.5\mu$m from the GOODS Re-ionization Era wide Area Treasury from Spitzer (GREATS) program and (2) consideration of $z\sim6-10$ sources over a $3\times$ larger area than previous HST+Spitzer studies. The Spitzer/IRAC data enable $\ge2\sigma$ rest-frame optical detections for an unprecedented $50\%$ of galaxies down to a stellar mass limit of $\sim10^{8}\mathcal{M}_\odot$ across all redshifts. Schechter fits to our volume densities suggest a combined evolution in characteristic mass $\mathcal{M}^*$ and normalization factor $\phi^*$ between $z\sim6$ and $z\sim8$. The stellar mass density (SMD) increases by $\sim1000\times$ in the $\sim500$ Myr between $z\sim10$ and $z\sim6$, with indications of a steeper evolution between $z\sim10$ and $z\sim8$, similar to the previously-reported trend of the star-formation rate density. Strikingly, abundance matching to the Bolshoi-Planck simulation indicates halo mass densities evolving at approximately the same rate as the SMD between $z\sim10$ and $z\sim4$. Our results show that the stellar-to-halo mass ratios, a proxy for the star-formation efficiency, do not change significantly over the huge stellar mass build-up occurred from $z\sim10$ to $z\sim6$, indicating that the assembly of stellar mass closely mirrors the build-up in halo mass in the first $\sim1$ Gyr of cosmic history. JWST is poised to extend these results into the "first galaxy" epoch at $z\gtrsim10$., Comment: Submitted to ApJ, comments welcome

Deep rest-optical observations are required to accurately constrain the stellar populations of $z\sim8$ galaxies. Due to significant limitations in the availability of such data for statistically complete samples, observational results have been limited to modest numbers of bright or lensed sources. To revolutionize the present characterization of $z\sim8$ galaxies, we exploit the ultradeep ($\sim27$ mag, $3\sigma$) Spitzer/IRAC $3.6\mu$m and $4.5\mu$m data, probing the rest-frame optical at $z\sim8$, over $\sim200$ arcmin$^2$ of the GOODS fields from the recently completed GOODS Re-ionization Era wide-Area Treasury from Spitzer (GREATS) program, combined with observations in the CANDELS UDS and COSMOS fields. We stacked $\gtrsim100$ $z\sim8$ Lyman-Break galaxies in four bins of UV luminosity ($M_\mathrm{UV}\sim -20.7$ to $-18.4$) and study their $H_\mathrm{160}-[3.6]$ and $[3.6]-[4.5]$ colors. We find young ages ($\lesssim100$ Myr) for the three faintest stacks, inferred from their blue $H_\mathrm{160}-[3.6]\sim 0$ mag colors, consistent with a negative Balmer break. Meanwhile, the redder $H_\mathrm{160}-[3.6]$ color seen in the brightest stack is suggestive of slightly older ages. We explored the existence of a correlation between the UV luminosity and age, and find either no trend or fainter galaxies being younger. The stacked SEDs also exhibit very red $[3.6]-[4.5]\sim0.5$ mag colors, indicative of intense [OIII]+H$\beta$ nebular emission and SFR. The correspondingly high specific star-formation rates, sSFR$\gtrsim10$Gyr$^{-1}$, are consistent with recent determinations at similar redshifts and higher luminosities, and support the co-evolution between the sSFR and the specific halo mass accretion rate., Comment: Submitted to ApJ. Comments welcome. See Fig. 2 & 3 for the stacks, Fig. 4 & 6 for the colors/age and Fig. 8 for the sSFR evolution

We present the Super Eight galaxies - a set of very luminous, high-redshift ($7.1

We present 16 new ultrabright $H_{AB}\lesssim25$ galaxy candidates at z~8 identified over the COSMOS/UltraVISTA field. The new search takes advantage of the deepest-available ground-based optical and near-infrared observations, including the DR3 release of UltraVISTA and full-depth Spitzer/IRAC observations from the SMUVS and SPLASH programs. Candidates are selected using Lyman-break criteria, combined with strict optical non-detection and SED-fitting criteria, minimizing contamination by low-redshift galaxies and low-mass stars. HST/WFC3 coverage from the DASH program reveals that one source evident in our ground-based near-IR data has significant substructure and may actually correspond to 3 separate z~8 objects, resulting in a sample of 18 galaxies, 10 of which seem to be fairly robust (with a >97% probability of being at z>7). The UV-continuum slope $\beta$ for the bright z~8 sample is $\beta=-2.2\pm0.6$, bluer but still consistent with that of similarly bright galaxies at z~6 ($\beta=-1.55\pm0.17$) and z~7 ($\beta=-1.75\pm0.18$). Their typical stellar masses are 10$^{9.1^{+0.5}_{-0.4}}M_{\odot}$, with the SFRs of $32^{+44}_{-32}M_{\odot}$/year, specific SFR of $4^{+8}_{-4}$ Gyr$^{-1}$, stellar ages of $\sim22^{+69}_{-22}$\,Myr, and low dust content A$_V=0.15^{+0.30}_{-0.15}$ mag. Using this sample we constrain the bright end of the z~8 UV luminosity function (LF). When combined with recent empty field LF estimates at z~8-9, the resulting z~8 LF can be equally well represented by either a Schechter or a double power-law (DPL) form. Assuming a Schechter parameterization, the best-fit characteristic magnitude is $M^*= -20.95^{+0.30}_{-0.35}$ mag with a very steep faint end slope $\alpha=-2.15^{+0.20}_{-0.19}$. These new candidates include amongst the brightest yet found at these redshifts, 0.5-1.0 mag brighter than found over CANDELS, providing excellent targets for follow-up studies., Comment: This version reflects the accepted manuscript

High-resolution multi-wavelength photometry is crucial to explore the spatial distribution of star formation in galaxies and understand how these evolve. To this aim, in this paper we exploit the deep, multi-wavelength Hubble Space Telescope (HST) data available in the central parts of the GOODS fields and study the distribution of star formation activity and mass in galaxies located at different positions with respect to the Main Sequence (MS) of star-forming galaxies. Our sample consists of galaxies with stellar mass $\geq 10^{9.5} M_{\odot}$ in the redshift range 0.2 $ \leq z \leq 1.2$. Exploiting 10-band photometry from the UV to the near-infrared at HST resolution, we derive spatially resolved maps of galaxies properties, such as stellar mass and star formation rate and specific star formation rate, with a resolution of $\sim 0.16$ arcsec. We find that the star formation activity is centrally enhanced in galaxies above the MS and centrally suppressed below the MS, with quiescent galaxies (1 dex below the MS) characterised by the highest suppression. The sSFR in the outer region does not show systematic trends of enhancement or suppression above or below the MS. The distribution of mass in MS galaxies indicates that bulges are growing when galaxies are still on the MS relation. Galaxies below the MS are more bulge-dominated with respect to MS counterparts at fixed stellar mass, while galaxies in the upper envelope are more extended and have S\'ersic indexes that are always smaller than or comparable to their MS counterparts. The suppression of star formation activity in the central region of galaxies below the MS hints at \textit{inside-out} quenching, as star formation is still ongoing in the outer regions., Comment: 17 pages, 12 figures, accepted for publication on A&A

We investigate the dependence of galaxy clustering at $z \sim 4 - 7$ on UV-luminosity and stellar mass. Our sample consists of $\sim$ 10,000 Lyman-break galaxies (LBGs) in the XDF and CANDELS fields. As part of our analysis, the $M_\star - M_{\rm UV}$ relation is estimated for the sample, which is found to have a nearly linear slope of $d\log_{10} M_\star / d M_{\rm UV} \sim 0.44$. We subsequently measure the angular correlation function and bias in different stellar mass and luminosity bins. We focus on comparing the clustering dependence on these two properties. While UV-luminosity is only related to recent starbursts of a galaxy, stellar mass reflects the integrated build-up of the whole star formation history, which should make it more tightly correlated with halo mass. Hence, the clustering segregation with stellar mass is expected to be larger than with luminosity. However, our measurements suggest that the segregation with luminosity is larger with $\simeq 90\%$ confidence (neglecting contributions from systematic errors). We compare this unexpected result with predictions from the \textsc{Meraxes} semi-analytic galaxy formation model. Interestingly, the model reproduces the observed angular correlation functions, and also suggests stronger clustering segregation with luminosity. The comparison between our observations and the model provides evidence of multiple halo occupation in the small scale clustering., Comment: 10 pages, 6 figures, 2 tables, accepted for publication in MNRAS

We report on the discovery of three especially bright candidate $z_{phot} \gtrsim 8$ galaxies. Five sources were targeted for follow-up with HST/WFC3, selected from a larger sample of 16 bright ($24.8 \lesssim H\lesssim25.5$~mag) candidate $z\gtrsim 8$ LBGs identified over the 1.6 degrees$^2$ of the COSMOS/UltraVISTA field. These were identified as Y and J dropouts by leveraging the deep (Y-to-$K_{S} \sim 25.3-24.8$~mag, $5\sigma$) NIR data from the UltraVISTA DR3 release, deep ground based optical imaging from the CFHTLS and Subaru Suprime Cam programs and Spitzer/IRAC mosaics combining observations from the SMUVS and SPLASH programs. Through the refined spectral energy distributions, which now also include new HyperSuprime Cam g, r, i, z and Y band data, we confirm that 3/5 galaxies have robust $z_{phot}\sim8.0-8.7$, consistent with the initial selection. The remaining 2/5 galaxies have a nominal $z_{phot}\sim2$. However, if we use the HST data alone, these objects have increased probability of being at $z\sim9$. Furthermore, we measure mean UV continuum slopes $\beta=-1.91\pm0.26$ for the three $z\sim8-9$ galaxies, marginally bluer than similarly luminous $z\sim4-6$ in CANDELS but consistent with previous measurements of similarly luminous galaxies at $z\sim7$. The circularized effective radius for our brightest source is $0.9\pm0.2$ kpc, similar to previous measurements for a bright $z\sim11$ galaxy and bright $z\sim7$ galaxies. Finally, enlarging our sample to include the six brightest $z\sim8$ LBGs identified over UltraVISTA (i.e., including three other sources from Labbe et al. 2017, in prep.) we estimate for the first time the volume density of galaxies at the extreme bright ($M_{UV}\sim-22$~mag) end of the $z\sim8$ UV LF. Despite this exceptional result, the still large statistical uncertainties do not allow us to discriminate between a Schechter and a double power-law form., Comment: Resubmitting after addressing the referee's comments. Appendix B includes the additional 3 sources from Labbe+2017, used for the LF estimate

The earliest galaxies are expected to emerge in the first billion years of the Universe during the Epoch of Reionization. However, both the spectroscopic confirmation of photometrically-selected galaxies at this epoch and the characterization of their early dynamical state has been hindered by the lack of bright, accessible lines to probe the velocity structure of their interstellar medium. We present the first ALMA spectroscopic confirmation of such sources at z > 6 using the far-infrared [C II]{\lambda}157.74{\mu}m emission line, and, for the first time, measurement of the velocity structure, for two galaxies at z = 6.8540+/-0.0003 and z = 6.8076+/-0.0002. Remarkably, the [C II] line luminosity from these galaxies is higher than previously found in `normal' star-forming galaxies at z > 6.5. This suggests that we are sampling a part of the galaxy population different from the galaxies found through detection of the Ly{\alpha} line. The luminous and extended [C II] detections reveal clear velocity gradients that, if interpreted as rotation, would suggest these galaxies have turbulent, yet rotation-dominated disks, with similar stellar-to-dynamical mass fractions as observed for H{\alpha} emitting galaxies 2 Gyr later at cosmic noon. Our novel approach for confirming galaxies during Reionization paves the way for larger studies of distant galaxies with spectroscopic redshifts from ALMA. Particularly important, this opens up opportunities for high angular-resolution [C II] dynamics in galaxies less than one billion years after the Big Bang., Comment: 14 pages, 5 figures, 1 table, published in Nature (submitted version)

We use a newly assembled large sample of 3,545 star-forming galaxies with secure spectroscopic, grism, and photometric redshifts at z=1.5-2.5 to constrain the relationship between UV slope (beta) and dust attenuation (L(IR)/L(UV)=IRX). Our sample benefits from the combination of deep Hubble WFC3/UVIS photometry from the Hubble Deep UV (HDUV) Legacy survey and existing photometric data compiled in the 3D-HST survey, and extends the range of UV luminosity and beta probed in previous UV-selected samples. IRX is measured using stacks of deep Herschel/PACS 100 and 160 micron data, and the results are compared with predictions of the IRX-beta relation for different assumptions of the stellar population model and obscuration curve. We find that z=1.5-2.5 galaxies have an IRX-beta relation that is consistent with the predictions for an SMC extinction curve if we invoke sub-solar metallicity models that are currently favored for high-redshift galaxies, while the commonly assumed starburst attenuation curve over-predicts the IRX at a given beta by a factor of ~3. The IRX of high-mass (M*>10^9.75 Msun) galaxies is a factor of >4 larger than that of low-mass galaxies, lending support for the use of stellar mass as a proxy for attenuation. The commonly observed trend of fainter galaxies having bluer beta may simply reflect bluer intrinsic UV slopes for such galaxies, rather than lower obscurations. The IRX-beta for young/low-mass galaxies implies a dust curve that is steeper than the SMC, suggesting a lower attenuation at a given beta relative to older/more massive galaxies. The lower attenuations and higher ionizing photon output implied by low metallicity stellar population models point to Lyman continuum production efficiencies, xi_ion, that may be elevated by a factor of ~2 relative to the canonical value for L* galaxies, aiding in their ability to keep the universe ionized at z~2. [Abridged], Comment: 22 pages, 3 tables, 11 figures, submitted to the Astrophysical Journal

Using Hubble data, including new grism spectra, Oesch et al. recently identified GN-z11, an $M_\textrm{UV}$=-21.1 galaxy at $z$=11.1 (just 400Myr after the big bang). With an estimated stellar mass of $\sim$10$^9$M$_{\odot}$, this galaxy is surprisingly bright and massive, raising questions as to how such an extreme object could form so early in the Universe. Using \Meraxes{}, a semi-analytic galaxy-formation model developed as part of the Dark-ages Reionization And Galaxy-formation Observables from Numerical Simulations (DRAGONS) programme, we investigate the potential formation mechanisms and eventual fate of GN-z11. The volume of our simulation is comparable to that of the discovery observations and possesses two analogue galaxies of similar luminosity to this remarkably bright system. Existing in the two most massive subhaloes at $z$=11.1 ($M_\textrm{vir}$=1.4$\times 10^{11}$M$_{\odot}$ and 6.7$\times 10^{10}$M$_{\odot}$), our model analogues show excellent agreement with all available observationally derived properties of GN-z11. Although they are relatively rare outliers from the full galaxy population at high-$z$, they are no longer the most massive or brightest systems by $z$=5. Furthermore, we find that both objects possess relatively smooth, but extremely rapid mass growth histories with consistently high star formation rates and UV luminosities at $z{>}11$, indicating that their brightness is not a transient, merger-driven feature. Our model results suggest that future wide-field surveys with the \textit{James Webb Space Telescope} may be able to detect the progenitors of GN-z11 analogues out to $z{\sim}$13--14, pushing the frontiers of galaxy-formation observations to the early phases of cosmic reionization and providing a valuable glimpse of the first galaxies to reionize the Universe on large scales., Comment: Updated to match published version. Minor correction to figure 2, with no resulting changes to findings. 8 pages, 6 figures, published in MNRAS

In this paper we study a key phase in the formation of massive galaxies: the transition of star forming galaxies into massive (M_stars~10^11 Msun), compact (r_e~1 kpc) quiescent galaxies, which takes place from z~3 to z~1.5. We use HST grism redshifts and extensive photometry in all five 3D-HST/CANDELS fields, more than doubling the area used previously for such studies, and combine these data with Keck MOSFIRE and NIRSPEC spectroscopy. We first confirm that a population of massive, compact, star forming galaxies exists at z~2, using K-band spectroscopy of 25 of these objects at 2.0

We present direct observational evidence for star formation quenching in galaxy groups in the redshift range 0

We present deep WFC3 grism observations of the candidate z~12 galaxy UDFj-39546284 in the HST Ultra Deep Field (UDF), by combining spectroscopic data from the 3D-HST and CANDELS surveys. The total exposure time is 40.5 ks and the spectrum covers 1.10 < {\lambda} < 1.65{\mu}m. We search for faint emission lines by cross-correlating the 2D G141 spectrum with the observed H160 morphology, a technique that is unique to slitless spectroscopy at HST resolution. We find a 2.7{\sigma} detection of an emission line at 1.599 {\mu}m---just redward of the JH140 filter---with flux 3.5 +/- 1.3x10-18 erg/s/cm2. Assuming the line is real, it contributes 110 +/- 40% of the observed H160 flux and has an observed equivalent width > 7300 {\AA}. If the line is confirmed, it could be Ly-{\alpha} at z = 12.12. However, a more plausible interpretation, given current results, could be a lower redshift feature such as [O III]{\lambda}4959,5007 at z = 2.19. We find two other 3D-HST [O III] emitters within 1000 km/s of that redshift in the GOODS-South field. Additional support for this interpretation comes from the discovery of a bright "[O III] blob" with a secure G141 grism redshift of z = 1.605. This object has a strikingly large observed equivalent width of nearly 9000 {\AA} that results in similar "dropout" colors as UDFj-39546284., Comment: Accepted for publication in the Astrophysical Journal Letters

We present stellar mass surface density profiles of a mass-selected sample of 177 galaxies at 0.5 < z < 2.5, obtained using very deep HST optical and near-infrared data over the GOODS-South field, including recent CANDELS data. Accurate stellar mass surface density profiles have been measured for the first time for a complete sample of high-redshift galaxies more massive than 10^10.7 M_sun. The key advantage of this study compared to previous work is that the surface brightness profiles are deconvolved for PSF smoothing, allowing accurate measurements of the structure of the galaxies. The surface brightness profiles account for contributions from complex galaxy structures such as rings and faint outer disks. Mass profiles are derived using radial rest-frame u-g color profiles and a well-established empirical relation between these colors and the stellar mass-to-light ratio. We derive stellar half-mass radii from the mass profiles, and find that these are on average ~25% smaller than rest-frame g band half-light radii. This average size difference of 25% is the same at all redshifts, and does not correlate with stellar mass, specific star formation rate, effective surface density, Sersic index, or galaxy size. Although on average the difference between half-mass size and half-light size is modest, for approximately 10% of massive galaxies this difference is more than a factor two. These extreme galaxies are mostly extended, disk-like systems with large central bulges. These results are robust, but could be impacted if the central dust extinction becomes high. ALMA observations can be used to explore this possibility. These results provide added support for galaxy growth scenarios wherein massive galaxies at these epochs grow by accretion onto their outer regions., Comment: 11 pages, 8 figures, 3 tables, accepted for publication in ApJ

For the first time, we study the evolution of the stellar mass-size relation for star-forming galaxies from z ~ 4 to z ~ 7 from Hubble-WFC3/IR camera observations of the HUDF and Early Release Science (ERS) field. The sizes are measured by determining the best fit model to galaxy images in the rest-frame 2100 \AA \ with the stellar masses estimated from SED fitting to rest-frame optical (from Spitzer/IRAC) and UV fluxes. We show that the stellar mass-size relation of Lyman-break galaxies (LBGs) persists, at least to z ~ 5, and the median size of LBGs at a given stellar mass increases towards lower redshifts. For galaxies with stellar masses of 9.5

We confirm the detection of 3 groups in the Lynx supercluster, at z~1.3, and give their redshifts and masses. We study the properties of the group galaxies as compared to the central clusters, RXJ0849+4452 and RXJ0848+4453, selecting 89 galaxies in the clusters and 74 galaxies in the groups. We morphologically classify galaxies by visual inspection, noting that our early-type galaxy (ETG) sample would have been contaminated at the 30% -40% level by simple automated classification methods (e.g. based on Sersic index). In luminosity selected samples, both clusters and groups show high fractions of Sa galaxies. The ETG fractions never rise above ~50% in the clusters, which is low compared to the fractions observed in clusters at z~1. However, ETG plus Sa fractions are similar to those observed for ETGs in clusters at z~1. Bulge-dominated galaxies visually classified as Sas might also be ETGs with tidal features or merger remnants. They are mainly red and passive, and span a large range in luminosity. Their star formation seems to have been quenched before experiencing a morphological transformation. Because their fraction is smaller at lower redshifts, they might be the spiral population that evolves into ETGs. For mass-selected samples, the ETG fraction show no significant evolution with respect to local clusters, suggesting that morphological transformations occur at lower masses and densities. The ETG mass-size relation shows evolution towards smaller sizes at higher redshift in both clusters and groups, while the late-type mass-size relation matches that observed locally. The group ETG red sequence shows lower zero points and larger scatters than in clusters, both expected to be an indication of a younger galaxy population. The estimated age difference is small when compared to the difference in age at different galaxy masses., Comment: ApJ, submitted - referee report answered - iterating with the referee

We combine recent estimates of dust extinction at z~4-7 with UV luminosity function (LF) determinations to derive star formation rate (SFR) functions at z~4, 5, 6 and 7. SFR functions provide a more physical description of galaxy build-up at high redshift and allow for direct comparisons to SFRs at lower redshifts determined by a variety of techniques. Our SFR functions are derived from well-established z~4-7 UV LFs, UV-continuum slope trends with redshift and luminosity, and IRX-beta relations. They are well-described by Schechter relations. We extend the comparison baseline for SFR functions to z~2 by considering recent determinations of the H{\alpha} and mid-IR luminosity functions. The low-end slopes of the SFR functions are flatter than for the UV LFs, \Delta\alpha\sim+0.13, and show no clear evolution with cosmic time (z~0-7). In addition, we find that the characteristic value SFR* from the Schechter fit to SFR function exhibits consistent, and substantial, linear growth as a function of redshift from ~5 M_sun/yr at z~8, 650 Myr after the Big Bang, to ~100 M_sun/yr at z~2, ~2.5 Gyr later. Recent results at z~10, close to the onset of galaxy formation, are consistent with this trend. The uniformity of this evolution is even greater than seen in the UV LF over the redshift range z~2-8, providing validation for our dust corrections. These results provide strong evidence that galaxies build up uniformly over the first 3 Gyr of cosmic time., Comment: Added an appendix, 1 figure and 3 tables: 9 pages, 5 figures, 4 tables, ApJ, in press

We present first results from the 3D-HST program, a near-IR spectroscopic survey performed with the Wide Field Camera 3 on the Hubble Space Telescope. We have used 3D-HST spectra to measure redshifts and Halpha equivalent widths for a stellar mass-limited sample of 34 galaxies at 110^11 M(sun) in the COSMOS, GOODS, and AEGIS fields. We find that a substantial fraction of massive galaxies at this epoch are forming stars at a high rate: the fraction of galaxies with Halpha equivalent widths >10 A is 59%, compared to 10% among SDSS galaxies of similar masses at z=0.1. Galaxies with weak Halpha emission show absorption lines typical of 2-4 Gyr old stellar populations. The structural parameters of the galaxies, derived from the associated WFC3 F140W imaging data, correlate with the presence of Halpha: quiescent galaxies are compact with high Sersic index and high inferred velocity dispersion, whereas star-forming galaxies are typically large two-armed spiral galaxies, with low Sersic index. Some of these star forming galaxies might be progenitors of the most massive S0 and Sa galaxies. Our results challenge the idea that galaxies at fixed mass form a homogeneous population with small scatter in their properties. Instead we find that massive galaxies form a highly diverse population at z>1, in marked contrast to the local Universe., Comment: Accepted for publication in ApJ Letters

We utilize for the first time HST ACS imaging to examine the structural properties of galaxies in the rest-frame U-V versus V-J diagram (i.e., the UVJ diagram) using a sample at 0.610.25). The use of the UVJ diagram as a tool to distinguish quiescent galaxies from star forming galaxies (SFGs) is becoming more common due to its ability to separate red quiescent galaxies from reddened SFGs. Quiescent galaxies occupy a small and distinct region of UVJ color space and we find most of them to have concentrated profiles with high Sersic indices (n>2.5) and smooth structure characteristic of early-type systems. SFGs populate a broad, but well-defined sequence of UVJ colors and are comprised of objects with a mix of Sersic indices. Interestingly, most UVJ-selected SFGs with high Sersic indices also display structure due to dust and star formation typical of the n<2.5 SFGs and late-type systems. Finally, we find that the position of a SFG on the sequence of UVJ colors is determined to a large degree by the mass of the galaxy and its inclination. Systems that are closer to edge-on generally display redder colors and lower [OII]3727 luminosity per unit mass as a consequence of the reddening due to dust within the disks. We conclude that the two main features seen in UVJ color space correspond closely to the traditional morphological classes of early and late-type galaxies., Comment: 7 pages, 4 figures, submitted to ApJ Letters

We use ultradeep HST WFC3/IR imaging of the HUDF to investigate the rest-frame optical morphologies of a mass-selected sample of galaxies at z~2. We find a large variety of galaxy morphologies, ranging from large, blue, disk-like galaxies to compact, red, early-type galaxies. We derive rest-frame u-g color profiles for these galaxies and show that most z~2 galaxies in our sample have negative color gradients such that their cores are red. Although these color gradients may partly be caused by radial variations in dust content, they point to the existence of older stellar populations in the centers of z~2 galaxies. This result is consistent with an "inside-out" scenario of galaxy growth. We find that the median color gradient is fairly constant with redshift: [Delta(u-g_{rest})/Delta(log r)]_{median} = -0.46, -0.44 and -0.49 for z~2, z~1 and z=0, respectively. Using structural parameters derived from surface brightness profiles we confirm that at z~2 galaxy morphology correlates well with specific star formation rate. At the same mass, star forming galaxies have larger effective radii, bluer rest-frame u-g colors and lower Sersic indices than quiescent galaxies. These correlations are very similar to those at lower redshift, suggesting that the relations that give rise to the Hubble sequence at z=0 are already in place for massive galaxies at this early epoch., Comment: Accepted for publication in ApJ Letters, 7 pages, 4 figures

We present deep near-infrared (NIR) medium-bandwidth photometry over the wavelength range 1-1.8 microns in the All-wavelength Extended Groth strip International Survey (AEGIS) and Cosmic Evolution Survey (COSMOS) fields. The observations were carried out as part of the NEWFIRM Medium-Band Survey (NMBS), an NOAO survey program on the Mayall 4m telescope on Kitt Peak using the NOAO Extremely Wide-Field Infrared Imager (NEWFIRM). In this paper, we describe the full details of the observations, data reduction and photometry for the survey. We also present a public K-selected photometric catalog, along with accurate photometric redshifts. The redshifts are computed with 37 (20) filters in the COSMOS (AEGIS) fields, combining the NIR medium-bandwidth data with existing ultraviolet (UV; Galaxy Evolution Explorer), visible and NIR (Canada-France-Hawaii Telescope and Subaru) and mid-IR (Spitzer/IRAC) imaging. We find excellent agreement with publicly available spectroscopic redshifts, with sigma_z/(1+z)~1-2% for ~4000 galaxies at z=0-3. The NMBS catalogs contain ~13,000 galaxies at z>1.5 with accurate photometric redshifts and rest-frame colors. Due to the increased spectral resolution obtained with the five NIR medium-band filters, the median 68% confidence intervals of the photometric redshifts of both quiescent and star-forming galaxies are a factor of ~2 times smaller when comparing catalogs with medium-band NIR photometry to NIR broadband photometry. We show evidence for a clear bimodal color distribution between quiescent and star-forming galaxies that persists to z~3, a higher redshift than has been probed so far., Comment: All NMBS data products and a high resolution version of paper are available for download at http://www.astro.yale.edu/nmbs; Accepted for publication in ApJ; 24 pages, 21 figures, 4 tables

We study the star formation rates (SFRs) of galaxies as a function of local galaxy density at 0.61.8x10^{10} Msun to conduct our main analysis. With three different SFR indicators, (1) Spitzer MIPS 24-micron imaging, (2) SED fitting, and (3) [OII]3727 emission, we find the median specific SFR (SSFR) and SFR to decline from the low-density field to the cores of groups and a rich cluster. For the SED and [OII] based SFRs, the decline in SSFR is roughly an order of magnitude while for the MIPS based SFRs, the decline is a factor of ~4. We find approximately the same magnitude of decline in SSFR even after removing the sample of galaxies near the cluster. Galaxies in groups and a cluster at these redshifts therefore have lower star formation (SF) activity than galaxies in the field, as is the case at z~0. We investigated whether the decline in SFR with increasing density is caused by a change in the proportion of quiescent and star forming galaxies (SFGs) or by a decline in the SFRs of SFGs. Using the rest-frame U-V and V-J colors to distinguish quiescent galaxies from SFGs we find the fraction of quiescent galaxies increases from ~32% to 79% from low to high density. In addition, we find the SSFRs of SFGs, selected based on U-V and V-J colors, to decline with increasing density by factors of ~5-6 for the SED and [OII] based SFRs. The MIPS based SSFRs for SFGs decline with a shallower slope. The order of magnitude decline in the SSFR-density relation at 0.6

We have measured velocity dispersions for a sample of 36 galaxies with J < 21.2 or Mr < -20.6 mag in MS1054-03, a massive cluster of galaxies at z = 0.83. Our data are of uniformly high quality down to our selection limit, our 16-hour exposures typically yielding errors of only \delta(dispersion)~10% for L* and fainter galaxies. By combining our measurements with data from the literature, we have 53 cluster galaxies with measured dispersions, and HST/ACS-derived sizes, colors and surface brightnesses. This sample is complete for the typical L* galaxy at z~1, unlike most previous z~1 cluster samples which are complete only for the massive cluster members (>1e11 M_sun). We find no evidence for a change in the tilt of the fundamental plane (FP). Nor do we find evidence for evolution in the slope of the color-dispersion relation and M/L_B-dispersion relations; measuring evolution at a fixed dispersion should minimize the impact of size evolution found in other work. The M/L_B at fixed dispersion evolves by \Delta log10 M/L_B=-0.50 +/- 0.03 between z=0.83 and z=0.02 or d(log10 M/L_B)=-0.60 +/- 0.04 dz, and we find \Delta (U-V)_z=-0.24 +/- 0.02 mag at fixed dispersion in the rest-frame, matching the expected evolution in M/L_B within 2.25 standard deviations. The implied formation redshift from both the color and M/L_B evolution is z*=2.0 +/- 0.2 +/- 0.3 (sys), during the epoch in which the cosmic star-formation activity peaked, with the systematic uncertainty showing the dependence of z* on the assumptions we make about the stellar populations. The lack of evolution in either the tilt of the FP or in the M/L- and color-dispersion relations imply that the formation epoch depends weakly on mass, ranging from z*=2.3 +1.3 -0.3 at 300 km/s to z*=1.7 +0.3 -0.2 at 160 km/s and implies that the IMF similarly varies slowly with galaxy mass., Comment: revised; typos corrected, references updated, and other cosmetic changes to meet ApJ format ApJ accepted, 22 pages in emulate ApJ format, 8 color figures, 1 b/w figure

[Abridged] We use the NEWFIRM Medium-Band Survey (NMBS) to characterize the properties of a mass-complete sample of 14 galaxies at 3.02.5x10^11 Msun, and to derive more accurate measurements of the high-mass end of the stellar mass function (SMF) of galaxies at z=3.5, with significantly reduced contributions from photometric redshift errors and cosmic variance to the total error budget of the SMF. The typical very massive galaxy at z=3.5 is red and faint in the observer's optical, with median r=26.1, and rest-frame U-V=1.6. About 60% of the sample have optical colors satisfying either the U- or the B-dropout color criteria, although ~50% of these galaxies have r>25.5. About 30% of the sample has SFRs from SED modeling consistent with zero. However, >80% of the sample is detected at 24 micron, with total infrared luminosities in the range (0.5-4.0)x10^13 Lsun. This implies the presence of either dust-enshrouded starburst activity (with SFRs of 600-4300 Msun/yr) and/or highly-obscured active galactic nuclei (AGN). The contribution of galaxies with M_star>2.5x10^11 Msun to the total stellar mass budget at z=3.5 is ~8%. We find an evolution by a factor of 2-7 and 3-22 from z~5 and z~6, respectively, to z=3.5. The previously found disagreement at the high-mass end between observed and model-predicted SMFs is now significant at the 3sigma level. However, systematic uncertainties dominate the total error budget, with errors up to a factor of ~8 in the densities, bringing the observed SMF in marginal agreement with the predicted SMF. Additional systematic uncertainties on the high-mass end could be introduced by either 1) the intense star-formation and/or the very common AGN activities as inferred from the MIPS 24 micron detections, and/or 2) contamination by a significant population of massive, old, and dusty galaxies at z~2.6., Comment: 20 pages, 11 figures. Accepted in ApJ. Minor changes to colors of figures to match accepted version

We present very deep Wide Field Camera 3 (WFC3) photometry of a massive, compact galaxy located in the Hubble Ultra Deep Field. This quiescent galaxy has a spectroscopic redshift z=1.91 and has been identified as an extremely compact galaxy by Daddi et al. 2005. We use new H-F160W imaging data obtained with Hubble Space Telescope/WFC3 to measure the deconvolved surface brightness profile to H = 28 mag arcsec**-2. We find that the surface brightness profile is well approximated by an n=3.7 Sersic profile. Our deconvolved profile is constructed by a new technique which corrects the best-fit Sersic profile with the residual of the fit to the observed image. This allows for galaxy profiles which deviate from a Sersic profile. We determine the effective radius of this galaxy: r_e=0.42 +- 0.14 kpc in the observed H-F160W-band. We show that this result is robust to deviations from the Sersic model used in the fit. We test the sensitivity of our analysis to faint "wings" in the profile using simulated galaxy images consisting of a bright compact component and a faint extended component. We find that due to the combination of the WFC3 imaging depth and our method's sensitivity to extended faint emission we can accurately trace the intrinsic surface brightness profile, and that we can therefore confidently rule out the existence of a faint extended envelope around the observed galaxy down to our surface brightness limit. These results confirm that the galaxy lies a factor of 10 off from the local mass-size relation., Comment: 6 pages, 4 figures. Accepted for publication in ApJ Letters.

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 fields

We study the growth of massive galaxies from z=2 to the present using data from the NEWFIRM Medium Band Survey. The sample is selected at a constant number density of n=2x10^-4 Mpc^-3, so that galaxies at different epochs can be compared in a meaningful way. We show that the stellar mass of galaxies at this number density has increased by a factor of ~2 since z=2, following the relation log(M)=11.45-0.15z. In order to determine at what physical radii this mass growth occurred we construct very deep stacked rest-frame R-band images at redshifts z=0.6, 1.1, 1.6, and 2.0. These image stacks of typically 70-80 galaxies enable us to characterize the stellar distribution to surface brightness limits of ~28.5 mag/arcsec^2. We find that massive galaxies gradually built up their outer regions over the past 10 Gyr. The mass within a radius of r=5 kpc is nearly constant with redshift whereas the mass at 5-75 kpc has increased by a factor of ~4 since z=2. Parameterizing the surface brightness profiles we find that the effective radius and Sersic n parameter evolve as r_e~(1+z)^-1.3 and n~(1+z)^-1.0 respectively. The data demonstrate that massive galaxies have grown mostly inside-out, assembling their extended stellar halos around compact, dense cores with possibly exponential radial density distributions. Comparing the observed mass evolution to the average star formation rates of the galaxies we find that the growth is likely dominated by mergers, as in-situ star formation can only account for ~20% of the mass build-up from z=2 to z=0. The main uncertainties in this study are possible redshift-dependent systematic errors in the total stellar masses and the conversion from light-weighted to mass-weighted radial profiles., Comment: Accepted for publication in the Astrophysical Journal. 26 pages, 13 figures in main text

We present optical spectroscopic follow-up of a sample of Distant Red Galaxies (DRGs) with K < 22.5 (Vega), selected by J-K > 2.3, in the Hubble Deep Field South, the MS 1054-03 field, and the Chandra Deep Field South. Spectroscopic redshifts were obtained for 15 DRGs. Only 2 out of 15 DRGs are located at z < 2, suggesting a high efficiency to select high-redshift sources. From other spectroscopic surveys in the CDFS targeting intermediate to high redshift populations selected with different criteria, we find spectroscopic redshifts for a further 30 DRGs. We use the sample of spectroscopically confirmed DRGs to establish the high quality (scatter in \Delta z/(1+z) of ~ 0.05) of their photometric redshifts in the considered deep fields, as derived with EAZY (Brammer et al. 2008). Combining the spectroscopic and photometric redshifts, we find that 74% of DRGs with K < 22.5 lie at z > 2. The combined spectroscopic and photometric sample is used to analyze the distinct intrinsic and observed properties of DRGs at z < 2 and z > 2. In our photometric sample to K < 22.5, low-redshift DRGs are brighter in K than high-redshift DRGs by 0.7 mag, and more extincted by 1.2 mag in Av. Our analysis shows that the DRG criterion selects galaxies with different properties at different redshifts. Such biases can be largely avoided by selecting galaxies based on their rest-frame properties, which requires very good multi-band photometry and high quality photometric redshifts., Comment: Accepted for publication in the Astrophysical Journal, 13 pages, 8 figures, 5 tables

We examine the star formation rates (SFRs) of galaxies in a redshift slice encompassing the z=0.834 cluster RX J0152.7-1357. We used a low-dispersion prism in the Inamori Magellan Areal Camera and Spectrograph (IMACS) to identify galaxies with z<23.3 AB mag in diverse environments around the cluster out to projected distances of ~8 Mpc from the cluster center. We utilize a mass-limited sample (M>2x10^{10} M_sun) of 330 galaxies that were imaged by Spitzer MIPS at 24 micron to derive SFRs and study the dependence of specific SFR (SSFR) on stellar mass and environment. We find that the SFR and SSFR show a strong decrease with increasing local density, similar to the relation at z~0. Our result contrasts with other work at z~1 that find the SFR-density trend to reverse for luminosity-limited samples. These other results appear to be driven by star-formation in lower mass systems (M~10^{10} M_sun). Our results imply that the processes that shut down star-formation are present in groups and other dense regions in the field. Our data also suggest that the lower SFRs of galaxies in higher density environments may reflect a change in the ratio of star-forming to non-star-forming galaxies, rather than a change in SFRs. As a consequence, the SFRs of star-forming galaxies, in environments ranging from small groups to clusters, appear to be similar and largely unaffected by the local processes that truncate star-formation at z~0.8., Comment: 5 pages, 3 figures, accepted for publication in ApJL

We present Hubble Space Telescope NIC2 morphologies of a spectroscopic sample of massive galaxies at z~2.3, by extending our sample of 9 compact quiescent galaxies (r_e~0.9 kpc) with 10 massive emission-line galaxies. The emission-line galaxies are classified by the nature of their ionized emission; there are six star-forming galaxies and four galaxies hosting an active galactic nucleus (AGN). The star-forming galaxies are the largest among the emission-line galaxies, with a median size of re = 2.8 kpc. The three galaxies with the highest star formation rates (> 100M \odot/yr) have irregular and clumpy morphologies. The AGN host galaxies are more similar to the compact quiescent galaxies in terms of their structures (re~1.1 kpc) and spectral energy distributions. The total sample clearly separates into two classes in a color-mass diagram: the large star-forming galaxies that form the blue cloud, and the compact quiescent galaxies on the red sequence. However, it is unclear how or even if the two classes are evolutionary related. Three out of six massive star-forming galaxies have dense cores and thus may passively evolve into compact galaxies due to fading of outer star-forming regions. For these galaxies a reverse scenario, in which compact galaxies grow inside-out by star formation is also plausible. We do caution though that the sample is small. Nonetheless, it is evident that a Hubble sequence of massive galaxies with strongly correlated galaxy properties is already in place at z > 2., Comment: Minor changes to match published version. Animated versions of some figures can be found at https://www.cfa.harvard.edu/~mkriek/Presentation_material.html. The reduced NIC2 images and corresponding PSFs of the full galaxy sample can be downloaded at https://www.cfa.harvard.edu/~mkriek/gnirs/nicmos/