This paper is one in a series that explores the importance of radius as a second parameter in galaxy evolution. The topic investigated here is the relationship between star formation rate (SFR) and galaxy radius ($R_{\rm e}$) for main-sequence star-forming galaxies. The key observational result is that, over a wide range of stellar mass and redshift in both CANDELS and SDSS, there is little trend between SFR and $R_{\rm e}$ at fixed stellar mass. The Kennicutt-Schmidt law, or any similar density-related star formation law, then implies that smaller galaxies must have lower gas fractions than larger galaxies (at fixed $M_{\ast}$), and this is supported by observations of local star-forming galaxies. We investigate the implication by adopting the equilibrium "bathtub" model: the ISM gas mass is assumed to be constant over time and the net star formation rate is the difference between the accretion rate of gas onto the galaxy from the halo and the outflow rate due to winds. To match the observed null correlation between SFR and radius, the bathtub model requires that smaller galaxies at fixed mass have weaker galactic winds. Our hypothesis is that galaxies are a 2-dimensional family whose properties are set mainly by halo mass and concentration. Galaxy radius and accretion rate plausibly both depend on halo concentration, which predicts how wind strength should vary with $R_{\rm e}$ and SFR.

We present a WFC3 F160W ($H$-band) selected catalog in the CANDELS/GOODS-N field containing photometry from the ultraviolet (UV) to the far-infrared (IR), photometric redshifts and stellar parameters derived from the analysis of the multi-wavelength data. The catalog contains 35,445 sources over the 171 arcmin$^{2}$ of the CANDELS F160W mosaic. The 5$\sigma$ detection limits (within an aperture of radius 0\farcs17) of the mosaic range between $H=27.8$, 28.2 and 28.7 in the wide, intermediate and deep regions, that span approximately 50\%, 15\% and 35\% of the total area. The multi-wavelength photometry includes broad-band data from UV (U band from KPNO and LBC), optical (HST/ACS F435W, F606W, F775W, F814W, and F850LP), near-to-mid IR (HST/WFC3 F105W, F125W, F140W and F160W, Subaru/MOIRCS Ks, CFHT/Megacam K, and \spitzer/IRAC 3.6, 4.5, 5.8, 8.0 $\mu$m) and far IR (\spitzer/MIPS 24$\mu$m, HERSCHEL/PACS 100 and 160$\mu$m, SPIRE 250, 350 and 500$\mu$m) observations. In addition, the catalog also includes optical medium-band data (R$\sim50$) in 25 consecutive bands, $\lambda=500$ to 950~nm, from the SHARDS survey and WFC3 IR spectroscopic observations with the G102 and G141 grisms (R$\sim210$ and 130). The use of higher spectral resolution data to estimate photometric redshifts provides very high, and nearly uniform, precision from $z=0-2.5$. The comparison to 1,485 good quality spectroscopic redshifts up to $z\sim3$ yields $\Delta z$/(1+$z_{\rm spec}$)$=$0.0032 and an outlier fraction of $\eta=$4.3\%. In addition to the multi-band photometry, we release added-value catalogs with emission line fluxes, stellar masses, dust attenuations, UV- and IR- based star formation rates and rest-frame colors.

Stellar masses in future James Webb Space Telescope (JWST) deep blank-fields will be mainly derived fitting the spectral energy distribution with theoretical galaxy templates. We investigate the uncertainties and biases of the stellar masses derived by using the LePhare code for spectral energy distribution fitting and the Yggdrasil theoretical templates. We consider a sample of mock galaxies at z$=$7-10 with mock JWST observations with S/N$_{F150W}\geqslant$ 10. Our goal is to provide a list of statistical stellar mass corrections to include on the stellar mass derivation for different output galaxy properties and JWST filter combinations to correct for template degeneracies. Median statistical stellar mass corrections vary from -0.83 dex to 0.87 dex, while 25$\%$ (75$\%$) quartiles range from -0.83 dex (-0.67 dex) to 0.51 dex (0.88 dex), depending on filter combinations and galaxy models. The most challenging cases are galaxies with nebular emission lines, especially the ones that are wrongly identified as galaxies without, relative dust-free galaxies and galaxies with small metallicities (i.e. Z$=1/50Z_{\odot}$). The stellar mass estimation of galaxies correctly identified without emission lines is generally fine, except at z$=$10 when considering only the 8 NIRCam bands, which make the MIRI bands very valuable. We have tested our stellar mass corrections using the public JAGUAR galaxy catalogue, deriving that the average discrepancy in the recovered stellar mass distribution decreases by 20-50$\%$ at z$>$7 after the correction. We found that without the stellar-mass corrections the number of low-mass galaxies (M$^{*}<10^{7}M_{\odot}$) is overestimated, which can potentially lead to systematic errors in the calculation of the galaxy stellar mass function faint-end slope at high z.

Galaxy mergers are expected to have a significant role in the mass assembly of galaxies in the early Universe, but there are very few observational constraints on the merger history of galaxies at $z>2$. We present the first study of galaxy major mergers (mass ratios $>$ 1:4) in mass-selected samples out to $z\approx6$. Using all five fields of the HST/CANDELS survey and a probabilistic pair count methodology that incorporates the full photometric redshift posteriors and corrections for stellar mass completeness, we measure galaxy pair-counts for projected separations between 5 and 30 kpc in stellar mass selected samples at $9.7 < \log_{10}(\rm{M}_{*}/\rm{M}_{\odot}) < 10.3$ and $\log_{10}(\rm{M}_{*}/\rm{M}_{\odot}) > 10.3$. We find that the major merger pair fraction rises with redshift to $z\approx6$ proportional to $(1+z)^{m}$, with $m = 0.8\pm0.2$ ($m = 1.8\pm0.2$) for $\log_{10}(\rm{M}_{*} / \rm{M}_{\odot}) > 10.3$ ($9.7 < \log_{10}(\rm{M}_{*}/\rm{M}_{\odot}) < 10.3$). Investigating the pair fraction as a function of mass ratio between 1:20 and 1:1, we find no evidence for a strong evolution in the relative numbers of minor to major mergers out to $z<3$. Using evolving merger timescales we find that the merger rate per galaxy ($\mathcal{R}$) rises rapidly from $0.07\pm 0.01$ Gyr$^{-1}$ at $z < 1$ to $7.6\pm 2.7$ Gyr$^{-1}$ at $z = 6$ for galaxies at $\log_{10}(\rm{M}_{*}/\rm{M}_{\odot}) > 10.3$. The corresponding co-moving major merger rate density remains roughly constant during this time, with rates of $\Gamma \approx 10^{-4}$ Gyr$^{-1}$ Mpc$^{-3}$. Based on the observed merger rates per galaxy, we infer specific mass accretion rates from major mergers that are comparable to the specific star-formation rates for the same mass galaxies at $z>3$ - observational evidence that mergers are as important a mechanism for building up mass at high redshift as in-situ star-formation.

Studying the formation and evolution of galaxies at the earliest cosmic times, and their role in reionization, requires the deepest imaging possible. Ultra-deep surveys like the HUDF and HFF have pushed to mag \mAB$\,\sim\,$30, revealing galaxies at the faint end of the LF to $z$$\,\sim\,$9$\,-\,$11 and constraining their role in reionization. However, a key limitation of these fields is their size, only a few arcminutes (less than a Mpc at these redshifts), too small to probe large-scale environments or clustering properties of these galaxies, crucial for advancing our understanding of reionization. Achieving HUDF-quality depth over areas $\sim$100 times larger becomes possible with a mission like the Wide Field Infrared Survey Telescope (WFIRST), a 2.4-m telescope with similar optical properties to HST, with a field of view of $\sim$1000 arcmin$^2$, $\sim$100$\times$ the area of the HST/ACS HUDF. This whitepaper motivates an Ultra-Deep Field survey with WFIRST, covering $\sim$100$\,-\,$300$\times$ the area of the HUDF, or up to $\sim$1 deg$^2$, to \mAB$\,\sim\,$30, potentially revealing thousands of galaxies and AGN at the faint end of the LF, at or beyond $z$\,$\sim$\,9$\,-\,$10 in the epoch of reionization, and tracing their LSS environments, dramatically increasing the discovery potential at these redshifts. (Note: This paper is a somewhat expanded version of one that was submitted as input to the Astro2020 Decadal Survey, with this version including an Appendix (which exceeded the Astro2020 page limits), describing how the science drivers for a WFIRST Ultra Deep Field might map into a notional observing program, including the filters used and exposure times needed to achieve these depths.)

We constrain the stellar population properties of a sample of 52 massive galaxies, with stellar mass log Ms>10.5, over the redshift range 0.5<z<2 by use of observer-frame optical and near-infrared slitless spectra from HST's ACS and WFC3 grisms. The deep exposures (~100 ks) allow us to target individual spectra of massive galaxies to F160W=22.5AB. Our spectral fitting approach uses a set of six base models adapted to the redshift and spectral resolution of each observation, and fits the weights of the base models, including potential dust attenuation, via an MCMC method. Our sample comprises a mixed distribution of quiescent (19) and star-forming galaxies (33). We quantify the width of the age distribution (Dt) that is found to dominate the variance of the retrieved parameters according to Principal Component Analysis. The population parameters follow the expected trend towards older ages with increasing mass, and Dt appears to weakly anti-correlate with stellar mass, suggesting a more efficient star formation at the massive end. As expected, the redshift dependence of the relative stellar age (measured in units of the age of the Universe at the source) in the quiescent sample rejects the hypothesis of a single burst (aka monolithic collapse). Radial colour gradients within each galaxy are also explored, finding a wider scatter in the star-forming subsample, but no conclusive trend with respect to the population parameters.

Recent near-IR power-spectra and panchromatic Extragalactic Background Light measurements provide upper limits on the near-IR surface brightness (SB>31 mag/arcsec^2) that may come from Pop III stars and accretion disks around resulting stellar-mass black holes (BHs) in the epoch of First Light (z=7-17). Physical parameters for zero metallicity Pop III stars at z>7 can be estimated from MESA stellar evolution models through helium-depletion, and for BH accretion disks from quasar microlensing results and multicolor accretion models. Second-generation stars can form at higher multiplicity, so that BH accretion disks may be fed by Roche-lobe overflow from lower-mass companions in their AGB stage. The near-IR SB constraints can be used to calculate the number of caustic transits behind lensing clusters that JWST and the 25~39 m ground-based telescopes may detect for both Pop III stars and stellar mass BH accretion disks. Because Pop III stars and stellar mass BH accretion disks have sizes of a few x 10^-11 arcsec at z>7, typical caustic magnifications can be mu=10^4~10^5, with rise times of hours and decline times of < 1 year for cluster transverse velocities of v<1000 km/s. Microlensing by intracluster medium objects can modify transit magnifications, and lengthen visibility times. Depending on BH masses, accretion-disk radii and feeding efficiencies, stellar-mass BH accretion-disk caustic transits could outnumber those from Pop III stars. To observe Pop III caustic transits directly may require monitoring 3~30 lensing clusters to AB< 29 mag over a decade or more. Such a program must be started with JWST in Cycle 1, and -- depending on the role of microlensing in the Intra Cluster Light -- should be continued for decades with the GMT and TMT, where JWST and the ground-based telescopes each will play a unique and strongly complementary role.

Galaxies grow their supermassive black holes in concert with their stars, although the relationship between these major galactic components is poorly understood. Observations of the cosmic growth of stars and black holes in galaxies suffer from disjoint samples and the strong effects of dust attenuation. The thermal infrared holds incredible potential for simultaneously measuring both the star formation and black hole accretion rates in large samples of galaxies covering a wide range of physical conditions. Spitzer demonstrated this potential at low redshift, and by observing some of the most luminous galaxies at z~2. JWST will apply these methods to normal galaxies at these epochs, but will not be able to generate large spectroscopic samples or access the thermal infrared at high-redshift. An order of magnitude gap in our wavelength coverage will persist between JWST and ALMA. A large, cold infrared telescope can fill this gap to determine when (in cosmic time), and where (within the cosmic web), stars and black holes co-evolve, by measuring these processes simultaneously in statistically complete and unbiased samples of galaxies to z>8. A next-generation radio interferometer will have the resolution and sensitivity to measure star-formation and nuclear accretion in even the dustiest galaxies. Together, the thermal infrared and radio can uniquely determine how stars and supermassive blackholes co-evolve in galaxies over cosmic time.

We derive direct measurement gas-phase metallicities of $7.4 < 12 + \log(O/H) < 8.4$ for 14 low-mass Emission Line Galaxies (ELGs) at $0.3 < z < 0.8$ identified in the Faint Infrared Grism Survey (FIGS). We use deep slitless G102 grism spectroscopy of the Hubble Ultra Deep Field (HUDF), dispersing light from all objects in the field at wavelengths between 0.85 and 1.15 microns. We run an automatic search routine on these spectra to robustly identify 71 emission line sources, using archival data from VLT/MUSE to measure additional lines and confirm redshifts. We identify 14 objects with $0.3 < z < 0.8$ with measurable O[III]$\lambda$4363 \AA\ emission lines in matching VLT/MUSE spectra. For these galaxies, we derive direct electron-temperature gas-phase metallicities with a range of $7.4 < 12 + \log(O/H) < 8.4$. With matching stellar masses in the range of $10^{7.9} M_{\odot} < M_{\star} < 10^{10.4} M_{\odot}$, we construct a mass-metallicity (MZ) relation and find that the relation is offset to lower metallicities compared to metallicities derived from alternative methods (e.g.,$R_{23}$, O3N2, N2O2) and continuum selected samples. Using star formation rates (SFR) derived from the $H\alpha$ emission line, we calculate our galaxies' position on the Fundamental Metallicity Relation (FMR), where we also find an offset toward lower metallicities. This demonstrates that this emission-line-selected sample probes objects of low stellar masses but even lower metallicities than many comparable surveys. We detect a trend suggesting galaxies with higher Specific Star Formation (SSFR) are more likely to have lower metallicity. This could be due to cold accretion of metal-poor gas that drives star formation, or could be because outflows of metal-rich stellar winds and SNe ejecta are more common in galaxies with higher SSFR.

We present the results from the application of a two-dimensional emission line detection method, EMission-line two-Dimensional (EM2D), to the near-infrared G102 grism observations obtained with the Wide-Field Camera 3 (WFC3) as part of the Cycle 22 {\em Hubble Space Telescope} Treasury Program: the Faint Infrared Grism Survey (FIGS). Using the EM2D method, we have assembled a catalog of emission line galaxies (ELGs) with resolved star formation from each of the four FIGS fields. Not only can one better assess the global properties of ELGs, but the EM2D method allows for the analysis and an improved study of the individual emission-line region {\it within} each galaxy. This paper includes a description of the methodology, advantages, and the first results of the EM2D method applied to ELGs in FIGS. The advantage of 2D emission line measurements includes significant improvement of galaxy redshift measurements, approaching the level of accuracy seen in high-spectral-resolution data, but with greater efficiency; and the ability to identify and measure the properties of multiple sites of star-formation and over scales of $\sim$ 1 kpc within individual galaxies out to z $\sim$ 4. The EM2D method also significantly improves the reliability of high-redshift ($z\sim7$) Lyman-$\alpha$ detections. Coupled with the wide field of view and high efficiency of space-based grism observations, EM2D provides a noteworthy improvement on the physical parameters that can be extracted from grism observations.

This is the first in a series of papers examining the demographics of star-forming galaxies at $0.2<z<2.5$ in CANDELS. We study 9,100 galaxies from GOODS-S and UDS having published values of redshifts, masses, star-formation rates (SFRs), and dust attenuation ($A_V$) derived from UV-optical SED fitting. In agreement with previous works, we find that the $UVJ$ colors of a galaxy are closely correlated with its specific star-formation rate (SSFR) and $A_V$. We define rotated $UVJ$ coordinate axes, termed $S_\mathrm{SED}$ and $C_\mathrm{SED}$, that are parallel and perpendicular to the star-forming sequence and derive a quantitative calibration that predicts SSFR from $C_\mathrm{SED}$ with an accuracy of ~0.2 dex. SFRs from UV-optical fitting and from UV+IR values based on Spitzer/MIPS 24 $\mu\mathrm{m}$ agree well overall, but systematic differences of order 0.2 dex exist at high and low redshifts. A novel plotting scheme conveys the evolution of multiple galaxy properties simultaneously, and dust growth, as well as star-formation decline and quenching, exhibit "mass-accelerated evolution" ("downsizing"). A population of transition galaxies below the star-forming main sequence is identified. These objects are located between star-forming and quiescent galaxies in $UVJ$ space and have lower $A_V$ and smaller radii than galaxies on the main sequence. Their properties are consistent with their being in transit between the two regions. The relative numbers of quenched, transition, and star-forming galaxies are given as a function of mass and redshift.

We improve the accuracy of photometric redshifts by including low-resolution spectral data from the G102 grism on the Hubble Space Telescope, which assists in redshift determination by further constraining the shape of the broadband Spectral Energy Disribution (SED) and identifying spectral features. The photometry used in the redshift fits includes near-IR photometry from FIGS+CANDELS, as well as optical data from ground-based surveys and HST ACS, and mid-IR data from Spitzer. We calculated the redshifts through the comparison of measured photometry with template galaxy models, using the EAZY photometric redshift code. For objects with F105W $< 26.5$ AB mag with a redshift range of $0 < z < 6$, we find a typical error of $\Delta z = 0.03 * (1+z)$ for the purely photometric redshifts; with the addition of FIGS spectra, these become $\Delta z = 0.02 * (1+z)$, an improvement of 50\%. Addition of grism data also reduces the outlier rate from 8\% to 7\% across all fields. With the more-accurate spectrophotometric redshifts (SPZs), we searched the FIGS fields for galaxy overdensities. We identified 24 overdensities across the 4 fields. The strongest overdensity, matching a spectroscopically identified cluster at $z=0.85$, has 28 potential member galaxies, of which 8 have previous spectroscopic confirmation, and features a corresponding X-ray signal. Another corresponding to a cluster at $z=1.84$ has 22 members, 18 of which are spectroscopically confirmed. Additionally, we find 4 overdensities that are detected at an equal or higher significance in at least one metric to the two confirmed clusters.

The rate of major galaxy-galaxy merging is theoretically predicted to steadily increase with redshift during the peak epoch of massive galaxy development ($1{\leq}z{\leq}3$). We use close-pair statistics to objectively study the incidence of massive galaxies (stellar $M_{1}{\geq}2{\times}10^{10}M_{\odot}$) hosting major companions ($1{\leq}M_{1}/M_{2}{\leq}4$; i.e., $<$4:1) at six epochs spanning $0{<}z{<}3$. We select companions from a nearly complete, mass-limited ($\geq5{\times}10^{9}M_{\odot}$) sample of 23,696 galaxies in the five CANDELS fields and the SDSS. Using $5-50$ kpc projected separation and close redshift proximity criteria, we find that the major companion fraction $f_{\mathrm{mc}}(z)$ based on stellar mass-ratio (MR) selection increases from 6% ($z{\sim}0$) to 16% ($z{\sim}0.8$), then turns over at $z{\sim}1$ and decreases to 7% ($z{\sim}3$). Instead, if we use a major F160W flux ratio (FR) selection, we find that $f_{\mathrm{mc}}(z)$ increases steadily until $z=3$ owing to increasing contamination from minor (MR$>$4:1) companions at $z>1$. We show that these evolutionary trends are statistically robust to changes in companion proximity. We find disagreements between published results are resolved when selection criteria are closely matched. If we compute merger rates using constant fraction-to-rate conversion factors ($C_{\mathrm{merg,pair}}{=}0.6$ and $T_{\mathrm{obs,pair}}{=}0.65\mathrm{Gyr}$), we find that MR rates disagree with theoretical predictions at $z{>}1.5$. Instead, if we use an evolving $T_{\mathrm{obs,pair}}(z){\propto}(1+z)^{-2}$ from Snyder et al., our MR-based rates agree with theory at $0{<}z{<}3$. Our analysis underscores the need for detailed calibration of $C_{\mathrm{merg,pair}}$ and $T_{\mathrm{obs,pair}}$ as a function of redshift, mass and companion selection criteria to better constrain the empirical major merger history.

We present the results of an unbiased search for Ly{\alpha} emission from continuum-selected 6 < z < 8 galaxies. Our dataset consists of 160 orbits of G102 slitless grism spectroscopy obtained with the Hubble Space Telescope (HST) Wide Field Camera 3 as part of the Faint Infrared Grism Survey (FIGS; PI: Malhotra), which obtains deep slitless spectra of all sources in four fields, and was designed to minimize contamination in observations of previously-identified high-redshift galaxy candidates. The FIGS data can potentially spectroscopically confirm the redshifts of galaxies, and as Ly{\alpha} emission is resonantly scattered by neutral gas, FIGS can also constrain the ionization state of the intergalactic medium (IGM) during the epoch of reionization. These data have sufficient depth to detect Ly{\alpha} emission in this epoch, as Tilvi et al. (2016) have published the FIGS detection of previously known (Finkelstein et al. 2013) Ly{\alpha} emission at z = 7.51. The FIGS data use five separate roll-angles of HST to mitigate the contamination by nearby galaxies. We created a method that accounts for and removes the contamination from surrounding galaxies, and also removes any dispersed continuum light from each individual spectrum (Pirzkal et al. 2017). We searched for significant (> 4{\sigma}) emission lines using two different automated detection methods, free of any visual inspection biases. Applying these methods on photometrically-selected high-redshift candidates between 6 < z < 8 we find two emission lines, one previously published by Tilvi et al. (2016), and a new line at 1.028{\mu}m. We identify this lines as Ly{\alpha} at z = 7.452 +/- 0.003. This newly spectroscopically confirmed galaxy has the highest Ly{\alpha} rest-frame equivalent width (EW Ly{\alpha}) yet published at z > 7 (140.3 +/- 19.0{\AA}).

While major mergers have long been proposed as a driver of both AGN activity and the M-sigma relation, studies of moderate to high redshift Seyfert-luminosity AGN hosts have found little evidence for enhanced rates of interactions. However, both theory and observation suggest that while these AGN may be fueled by stochastic accretion and secular processes, high-luminosity, high-redshift, and heavily obscured AGN are the AGN most likely to be merger-driven. To better sample this population of AGN, we turn to infrared selection in the CANDELS/COSMOS field. Compared to their lower-luminosity and less obscured X-ray-only counterparts, IR-only AGN (luminous, heavily obscured AGN) are more likely to be classified as either irregular (50$^{+12}_{-12}$% vs. 9$^{+5}_{-2}$%) or asymmetric (69$^{+9}_{-13}$% vs. 17$^{+6}_{-4}$%) and are less likely to have a spheroidal component (31$^{+13}_{-9}$% vs. 77$^{+4}_{-6}$%). Furthermore, IR-only AGN are also significantly more likely than X-ray-only AGN (75$^{+8}_{-13}$% vs. 31$^{+6}_{-6}$%) to be classified either as interacting or merging in a way that significantly disturbs the host galaxy or disturbed though not clearly interacting or merging, which potentially represents the late stages of a major merger. This suggests that while major mergers may not contribute significantly to the fueling of Seyfert luminosity AGN, interactions appear to play a more dominant role in the triggering and fueling of high-luminosity heavily obscured AGN.

The analysis of galaxies on the star formation rate - stellar mass (SFR-M*) plane is a powerful diagnostic for galaxy evolution at different cosmic times. We consider a sample of 24463 galaxies from the CANDELS/GOODS-S survey to conduct a detailed analysis of the SFR-M* relation at redshifts 0.5$\leqslant z<$3 over more than three dex in stellar mass. To obtain SFR estimates, we utilise mid- and far-IR photometry when available, and rest-UV fluxes for all the other galaxies. We perform our analysis in different redshift bins, with two different methods: 1) a linear regression fitting of all star-forming galaxies, defined as those with specific star formation rates $\rm log_{10}(sSFR/yr^{-1}) > -9.8$, similarly to what is typically done in the literature; 2) a multi-Gaussian decomposition to identify the galaxy main sequence (MS), the starburst sequence and the quenched galaxy cloud. We find that the MS slope becomes flatter when higher stellar mass cuts are adopted, and that the apparent slope change observed at high masses depends on the SFR estimation method. In addition, the multi-Gaussian decomposition reveals the presence of a starburst population which increases towards low stellar masses and high redshifts. We find that starbursts make up ~5% of all galaxies at z=0.5-1.0, while they account for ~16% of galaxies at 2$<z<$3 with log$_{10}(M^{*})=$8.25-11.25. We conclude that the dissection of the SFR-M* in multiple components over a wide range of stellar masses is necessary to understand the importance of the different modes of star formation through cosmic time.

We examine the fraction of massive ($M_{*}>10^{10} M_{\odot}$), compact star-forming galaxies (cSFGs) that host an active galactic nucleus (AGN) at $z\sim2$. These cSFGs are likely the direct progenitors of the compact quiescent galaxies observed at this epoch, which are the first population of passive galaxies to appear in large numbers in the early Universe. We identify cSFGs that host an AGN using a combination of Hubble WFC3 imaging and Chandra X-ray observations in four fields: the Chandra Deep Fields, the Extended Groth Strip, and the UKIDSS Ultra Deep Survey field. We find that $39.2^{+3.9}_{-3.6}$\% (65/166) of cSFGs at $1.4<z<3.0$ host an X-ray detected AGN. This fraction is 3.2 times higher than the incidence of AGN in extended star-forming galaxies with similar masses at these redshifts. This difference is significant at the $6.2\sigma$ level. Our results are consistent with models in which cSFGs are formed through a dissipative contraction that triggers a compact starburst and concurrent growth of the central black hole. We also discuss our findings in the context of cosmological galaxy evolution simulations that require feedback energy to rapidly quench cSFGs. We show that the AGN fraction peaks precisely where energy injection is needed to reproduce the decline in the number density of cSFGs with redshift. Our results suggest that the first abundant population of massive, quenched galaxies emerged directly following a phase of elevated supermassive black hole growth and further hints at a possible connection between AGN and the rapid quenching of star formation in these galaxies.

We present an analysis of 15 Type Ia supernovae (SNe Ia) at redshift z > 1 (9 at 1.5 < z < 2.3) recently discovered in the CANDELS and CLASH Multi-Cycle Treasury programs using WFC3 on the Hubble Space Telescope. We combine these SNe Ia with a new compilation of 1050 SNe Ia, jointly calibrated and corrected for simulated survey biases to produce accurate distance measurements. We present unbiased constraints on the expansion rate at six redshifts in the range 0.07 < z < 1.5 based only on this combined SN Ia sample. The added leverage of our new sample at z > 1.5 leads to a factor of ~3 improvement in the determination of the expansion rate at z = 1.5, reducing its uncertainty to ~20%, a measurement of H(z=1.5)/H0=2.67 (+0.83,-0.52). We then demonstrate that these six measurements alone provide a nearly identical characterization of dark energy as the full SN sample, making them an efficient compression of the SN Ia data. The new sample of SNe Ia at z > 1 usefully distinguishes between alternative cosmological models and unmodeled evolution of the SN Ia distance indicators, placing empirical limits on the latter. Finally, employing a realistic simulation of a potential WFIRST SN survey observing strategy, we forecast optimistic future constraints on the expansion rate from SNe Ia.

We study galactic star-formation activity as a function of environment and stellar mass over 0.5<z<2.0 using the FourStar Galaxy Evolution (ZFOURGE) survey. We estimate the galaxy environment using a Bayesian-motivated measure of the distance to the third nearest neighbor for galaxies to the stellar mass completeness of our survey, $\log(M/M_\odot)>9 (9.5)$ at z=1.3 (2.0). This method, when applied to a mock catalog with the photometric-redshift precision ($\sigma_z / (1+z) \lesssim 0.02$), recovers galaxies in low- and high-density environments accurately. We quantify the environmental quenching efficiency, and show that at z> 0.5 it depends on galaxy stellar mass, demonstrating that the effects of quenching related to (stellar) mass and environment are not separable. In high-density environments, the mass and environmental quenching efficiencies are comparable for massive galaxies ($\log (M/M_\odot)\gtrsim$ 10.5) at all redshifts. For lower mass galaxies ($\log (M/M)_\odot) \lesssim$ 10), the environmental quenching efficiency is very low at $z\gtrsim$ 1.5, but increases rapidly with decreasing redshift. Environmental quenching can account for nearly all quiescent lower mass galaxies ($\log(M/M_\odot) \sim$ 9-10), which appear primarily at $z\lesssim$ 1.0. The morphologies of lower mass quiescent galaxies are inconsistent with those expected of recently quenched star-forming galaxies. Some environmental process must transform the morphologies on similar timescales as the environmental quenching itself. The evolution of the environmental quenching favors models that combine gas starvation (as galaxies become satellites) with gas exhaustion through star-formation and outflows ("overconsumption"), and additional processes such as galaxy interactions, tidal stripping and disk fading to account for the morphological differences between the quiescent and star-forming galaxy populations.

We explore observational and theoretical constraints on how galaxies might transition between the "star-forming main sequence" (SFMS) and varying "degrees of quiescence" out to $z=3$. Our analysis is focused on galaxies with stellar mass $M_*>10^{10}M_{\odot}$, and is enabled by GAMA and CANDELS observations, a semi-analytic model (SAM) of galaxy formation, and a cosmological hydrodynamical "zoom in" simulation with momentum-driven AGN feedback. In both the observations and the SAM, transition galaxies tend to have intermediate S\'ersic indices, half-light radii, and surface stellar mass densities compared to star-forming and quiescent galaxies out to $z=3$. We place an observational upper limit on the average population transition timescale as a function of redshift, finding that the average high-redshift galaxy is on a "fast track" for quenching whereas the average low-redshift galaxy is on a "slow track" for quenching. We qualitatively identify four physical origin scenarios for transition galaxies in the SAM: oscillations on the SFMS, slow quenching, fast quenching, and rejuvenation. Quenching timescales in both the SAM and the hydrodynamical simulation are not fast enough to reproduce the quiescent population that we observe at $z\sim3$. In the SAM, we do not find a clear-cut morphological dependence of quenching timescales, but we do predict that the mean stellar ages, cold gas fractions, SMBH masses, and halo masses of transition galaxies tend to be intermediate relative to those of star-forming and quiescent galaxies at $z<3$.

The bulk of the stellar growth over cosmic time is dominated by IR luminous galaxies at cosmic noon (z=1-2), many of which harbor a hidden active galactic nucleus (AGN). We use state of the art infrared color diagnostics, combining Spitzer and Herschel observations, to separate dust-obscured AGN from dusty star forming galaxies (SFGs) in the CANDELS and COSMOS surveys. We calculate 24 micron counts of SFGs, AGN/star forming "Composites", and AGN. AGN and Composites dominate the counts above 0.8 mJy at 24 micron, and Composites form at least 25% of an IR sample even to faint detection limits. We develop methods to use the Mid-Infrared Instrument (MIRI) on JWST to identify dust-obscured AGN and Composite galaxies from z~1-2. With the sensitivity and spacing of MIRI filters, we will detect >4 times as many AGN hosts than with Spitzer/IRAC criteria. Any star formation rates based on the 7.7 micron PAH feature (likely to be applied to MIRI photometry) must be corrected for the contribution of the AGN, or the SFR will be overestimated by ~35% for cases where the AGN provides half the IR luminosity and ~50% when the AGN accounts for 90% of the luminosity. Finally, we demonstrate that our MIRI color technique can select AGN with an Eddington ratio of $\lambda_{\rm Edd}\sim0.01$ and will identify AGN hosts with a higher sSFR than X-ray techniques alone. JWST/MIRI will enable critical steps forward in identifying and understanding dust-obscured AGN and the link to their host galaxies.

The Faint Infrared Grism Survey (FIGS) is a deep Hubble Space Telescope (HST) WFC3/IR (Wide Field Camera 3 Infrared) slitless spectroscopic survey of four deep fields. Two fields are located in the Great Observatories Origins Deep Survey-North (GOODS-N) area and two fields are located in the Great Observatories Origins Deep Survey-South (GOODS-S) area. One of the southern fields selected is the Hubble Ultra Deep Field. Each of these four fields were observed using the WFC3/G102 grism (0.8$\mu m$-1.15$\mu m$ continuous coverage) with a total exposure time of 40 orbits (~ 100 kilo-seconds) per field. This reaches a 3 sigma continuum depth of ~26 AB magnitudes and probes emission lines to $\approx 10^{-17}\ erg\ s^{-1} \ cm^{-2}$. This paper details the four FIGS fields and the overall observational strategy of the project. A detailed description of the Simulation Based Extraction (SBE) method used to extract and combine over 10000 spectra of over 2000 distinct sources brighter than m_F105W=26.5 mag is provided. High fidelity simulations of the observations is shown to significantly improve the background subtraction process, the spectral contamination estimates, and the final flux calibration. This allows for the combination of multiple spectra to produce a final high quality, deep, 1D-spectra for each object in the survey.

Imaging with the James Webb Space Telescope (JWST) will allow for observing the bulk of distant galaxies at the epoch of reionisation. The recovery of their properties, such as age, color excess E(B-V), specific star formation rate (sSFR) and stellar mass, will mostly rely on spectral energy distribution fitting, based on the data provided by JWST's two imager cameras, namely the Near Infrared Camera (NIRCam) and the Mid Infrared Imager (MIRI). In this work we analyze the effect of choosing different combinations of NIRCam and MIRI broad-band filters, from 0.6 {\mu}m to 7.7 {\mu}m, on the recovery of these galaxy properties. We performed our tests on a sample of 1542 simulated galaxies, with known input properties, at z=7-10. We found that, with only 8 NIRCam broad-bands, we can recover the galaxy age within 0.1 Gyr and the color excess within 0.06 mag for 70% of the galaxies. Besides, the stellar masses and sSFR are recovered within 0.2 and 0.3 dex, respectively, at z=7-9. Instead, at z=10, no NIRCam band traces purely the {\lambda}> 4000 {\AA} regime and the percentage of outliers in stellar mass (sSFR) increases by > 20% (> 90%), in comparison to z=9. The MIRI F560W and F770W bands are crucial to improve the stellar mass and the sSFR estimation at z=10. When nebular emission lines are present, deriving correct galaxy properties is challenging, at any redshift and with any band combination. In particular, the stellar mass is systematically overestimated in up to 0.3 dex on average with NIRCam data alone and including MIRI observations improves only marginally the estimation.

We derive relations between the effective radii $R_{\rm{eff}}$ of galaxies and the virial radii $R_{200c}$ of their dark matter halos over the redshift range $0 < z < 3$. For galaxies, we use the measured sizes from deep images taken with \emph{Hubble Space Telescope} for the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey; for halos, we use the inferred sizes from abundance matching to cosmological dark matter simulations via a stellar mass--halo mass (SMHM) relation. For this purpose, we derive a new SMHM relation based on the same selection criteria and other assumptions as for our sample of galaxies with size measurements. As a check on the robustness of our results, we also derive $R_{\rm{eff}}$--$R_{200c}$ relations for three independent SMHM relations from the literature. We find that galaxy $R_{\rm{eff}}$ is proportional on average to halo $R_{200c}$, confirming and extending to high redshifts the $z=0$ results of Kravtsov. Late-type galaxies (with low S\'ersic index and high specific star formation rate [sSFR]) follow a linear $R_{\rm{eff}}$--$R_{200c}$ relation, with effective radii at $0.5 < z < 3$ close to those predicted by simple models of disk formation; at $z < 0.5$, the sizes of late-type galaxies appear to be slightly below this prediction. Early-type galaxies (with high S\'ersic index and low sSFR) follow a roughly parallel $R_{\rm{eff}}$--$R_{200c}$ relation, $\sim$ 0.2--0.3 dex below the one for late-type galaxies. Our observational results, reinforced by recent hydrodynamical simulations, indicate that galaxies grow quasi-homologously with their dark matter halos.

We present a 0.4-8$\mu$m multi-wavelength photometric catalog in the Extended Groth Strip (EGS) field. This catalog is built on the Hubble Space Telescope (HST) WFC3 and ACS data from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS), and it incorporates the existing HST data from the All-wavelength Extended Groth strip International Survey (AEGIS) and the 3D-HST program. The catalog is based on detections in the F160W band reaching a depth of F160W=26.62 AB (90% completeness, point-sources). It includes the photometry for 41457 objects over an area of $\approx 206$ arcmin$^2$ in the following bands: HST ACS F606W and F814W; HST WFC3 F125W, F140W and F160W; CFHT/Megacam $u^*$, $g'$, $r'$, $i'$ and $z'$; CFHT/WIRCAM $J$, $H$ and $K_\mathrm{S}$; Mayall/NEWFIRM $J1$, $J2$, $J3$, $H1$, $H2$, $K$; Spitzer IRAC $3.6\mu$m, $4.5\mu$m, $5.8\mu$m and $8.0\mu$m. We are also releasing value-added catalogs that provide robust photometric redshifts and stellar mass measurements. The catalogs are publicly available through the CANDELS repository.

The reionization of the Universe is one of the most important topics of present day astrophysical research. The most plausible candidates for the reionization process are star-forming galaxies, which according to the predictions of the majority of the theoretical and semi-analytical models should dominate the HI ionizing background at z~3. We aim at measuring the Lyman continuum escape fraction, which is one of the key parameters to compute the contribution of star-forming galaxies to the UV background. We have used ultra-deep U-band imaging (U=30.2mag at 1sigma) by LBC/LBT in the CANDELS/GOODS-North field, as well as deep imaging in COSMOS and EGS fields, in order to estimate the Lyman continuum escape fraction of 69 star-forming galaxies with secure spectroscopic redshifts at 3.27<z<3.40 to faint magnitude limits (L=0.2L*, or equivalently M1500~-19). We have measured through stacks a stringent upper limit (<1.7% at 1sigma) for the relative escape fraction of HI ionizing photons from bright galaxies (L>L*), while for the faint population (L=0.2L*) the limit to the escape fraction is ~10%. We have computed the contribution of star-forming galaxies to the observed UV background at z~3 and we have found that it is not enough to keep the Universe ionized at these redshifts, unless their escape fraction increases significantly (>10%) at low luminosities (M1500>-19). We compare our results on the Lyman continuum escape fraction of high-z galaxies with recent estimates in the literature and discuss future prospects to shed light on the end of the Dark Ages. In the future, strong gravitational lensing will be fundamental to measure the Lyman continuum escape fraction down to faint magnitudes (M1500~-16) which are inaccessible with the present instrumentation on blank fields.

We present the rest-1500\AA\ UV luminosity functions (LF) for star-forming galaxies during the cosmic \textit{high noon} -- the peak of cosmic star formation rate at $1.5<z<3$. We use deep NUV imaging data obtained as part of the \textit{Hubble} Ultra-Violet Ultra Deep Field (UVUDF) program, along with existing deep optical and NIR coverage on the HUDF. We select F225W, F275W and F336W dropout samples using the Lyman break technique, along with samples in the corresponding redshift ranges selected using photometric redshifts and measure the rest-frame UV LF at $z\sim1.7,2.2,3.0$ respectively, using the modified maximum likelihood estimator. We perform simulations to quantify the survey and sample incompleteness for the UVUDF samples to correct the effective volume calculations for the LF. We select galaxies down to $M_{UV}=-15.9,-16.3,-16.8$ and fit a faint-end slope of $\alpha=-1.20^{+0.10}_{-0.13}, -1.32^{+0.10}_{-0.14}, -1.39^{+0.08}_{-0.12}$ at $1.4<z<1.9$, $1.8<z<2.6$, and $2.4<z<3.6$, respectively. We compare the star formation properties of $z\sim2$ galaxies from these UV observations with results from H\alpha\ and UV$+$IR observations. We find a lack of high SFR sources in the UV LF compared to the H\alpha\ and UV$+$IR, likely due to dusty SFGs not being properly accounted for by the generic $IRX-\beta$ relation used to correct for dust. We compute a volume-averaged UV-to-H\alpha\ ratio by \textit{abundance matching} the rest-frame UV LF and H\alpha\ LF. We find an increasing UV-to-H\alpha\ ratio towards low mass galaxies ($M_\star \lesssim 5\times10^9$ M$_\odot$). We conclude that this could be due to a larger contribution from starbursting galaxies compared to the high-mass end.

We present an investigation of clumpy galaxies in the Hubble Ultra Deep Field at 0.5 $\leq z \leq$ 1.5 in the rest-frame far-ultraviolet (FUV) using HST WFC3 broadband imaging in F225W, F275W, and F336W. An analysis of 1,404 galaxies yields 209 galaxies that host 403 kpc-scale clumps. These host galaxies appear to be typical star-forming galaxies, with an average of 2 clumps per galaxy and reaching a maximum of 8 clumps. We measure the photometry of the clumps, and determine the mass, age, and star formation rates (SFR) utilizing the SED-fitting code FAST. We find that clumps make an average contribution of 19% to the total rest-frame FUV flux of their host galaxy. Individually, clumps contribute a median of 5% to the host galaxy SFR and an average of $\sim$4% to the host galaxy mass, with total clump contributions to the host galaxy stellar mass ranging widely from less than 1% up to 93%. Clumps in the outskirts of galaxies are typically younger, with higher star formation rates, than clumps in the inner regions. The results are consistent with clump migration theories in which clumps form through violent gravitational instabilities in gas-rich turbulent disks, eventually migrate toward the center of the galaxies, and coalesce into the bulge.

We present a multi-wavelength photometric catalog in the COSMOS field as part of the observations by the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). The catalog is based on Hubble Space Telescope Wide Field Camera 3 (HST/WFC3) and Advanced Camera for Surveys (ACS) observations of the COSMOS field (centered at RA: $10^h00^m28^s$, Dec:$+02^{\circ}12^{\prime}21^{\prime\prime}$). The final catalog has 38671 sources with photometric data in forty two bands from UV to the infrared ($\rm \sim 0.3-8\,\mu m$). This includes broad-band photometry from the HST, CFHT, Subaru, VISTA and Spitzer Space Telescope in the visible, near infrared and infrared bands along with intermediate and narrow-band photometry from Subaru and medium band data from Mayall NEWFIRM. Source detection was conducted in the WFC3 F160W band (at $\rm 1.6\,\mu m$) and photometry is generated using the Template FITting algorithm. We further present a catalog of the physical properties of sources as identified in the HST F160W band and measured from the multi-band photometry by fitting the observed spectral energy distributions of sources against templates.

Galaxies at high redshifts provide a valuable tool to study cosmic dawn, and therefore it is crucial to reliably identify these galaxies. Here, we present an unambiguous and first simultaneous detection of both the Lyman-alpha emission and the Lyman break from a z = 7.512+/- 0.004 galaxy, observed in the Faint Infrared Grism Survey (FIGS). These spectra, taken with G102 grism on Hubble Space Telescope (HST), show a significant emission line detection (6 sigma) in multiple observational position angles (PA), with total integrated Ly{\alpha} line flux of 1.06+/- 0.12 e10-17erg s-1cm-2. The line flux is nearly a factor of four higher than the previous MOSFIRE spectroscopic observations of faint Ly{\alpha} emission at {\lambda} = 1.0347{\mu}m, yielding z = 7.5078+/- 0.0004. This is consistent with other recent observations implying that ground-based near-infrared spectroscopy underestimates total emission line fluxes, and if confirmed, can have strong implications for reionization studies that are based on ground-based Lyman-{\alpha} measurements. A 4-{\sigma} detection of the NV line in one PA also suggests a weak Active Galactic Nucleus (AGN), potentially making this source the highest-redshift AGN yet found. Thus, this observation from the Hubble Space Telescope clearly demonstrates the sensitivity of the FIGS survey, and the capability of grism spectroscopy to study the epoch of reionization.

The determination of galaxy redshifts in James Webb Space Telescope (JWST)'s blank-field surveys will mostly rely on photometric estimates, based on the data provided by JWST's Near-Infrared Camera (NIRCam) at 0.6-5.0 {\mu}m and Mid Infrared Instrument (MIRI) at {\lambda}>5.0 {\mu}m. In this work we analyse the impact of choosing different combinations of NIRCam and MIRI broad-band filters (F070W to F770W), as well as having ancillary data at {\lambda}<0.6 {\mu}m, on the derived photometric redshifts (zphot) of a total of 5921 real and simulated galaxies, with known input redshifts z=0-10. We found that observations at {\lambda}<0.6 {\mu}m are necessary to control the contamination of high-z samples by low-z interlopers. Adding MIRI (F560W and F770W) photometry to the NIRCam data mitigates the absence of ancillary observations at {\lambda}<0.6 {\mu}m and improves the redshift estimation. At z=7-10, accurate zphot can be obtained with the NIRCam broad bands alone when S/N>=10, but the zphot quality significantly degrades at S/N<=5. Adding MIRI photometry with one magnitude brighter depth than the NIRCam depth allows for a redshift recovery of 83-99%, depending on SED type, and its effect is particularly noteworthy for galaxies with nebular emission. The vast majority of NIRCam galaxies with [F150W]=29 AB mag at z=7-10 will be detected with MIRI at [F560W, F770W]<28 mag if these sources are at least mildly evolved or have spectra with emission lines boosting the mid-infrared fluxes.

We present quantified visual morphologies of approximately 48,000 galaxies observed in three Hubble Space Telescope legacy fields by the Cosmic And Near-infrared Deep Extragalactic Legacy Survey (CANDELS) and classified by participants in the Galaxy Zoo project. 90% of galaxies have z < 3 and are observed in rest-frame optical wavelengths by CANDELS. Each galaxy received an average of 40 independent classifications, which we combine into detailed morphological information on galaxy features such as clumpiness, bar instabilities, spiral structure, and merger and tidal signatures. We apply a consensus-based classifier weighting method that preserves classifier independence while effectively down-weighting significantly outlying classifications. After analysing the effect of varying image depth on reported classifications, we also provide depth-corrected classifications which both preserve the information in the deepest observations and also enable the use of classifications at comparable depths across the full survey. Comparing the Galaxy Zoo classifications to previous classifications of the same galaxies shows very good agreement; for some applications the high number of independent classifications provided by Galaxy Zoo provides an advantage in selecting galaxies with a particular morphological profile, while in others the combination of Galaxy Zoo with other classifications is a more promising approach than using any one method alone. We combine the Galaxy Zoo classifications of "smooth" galaxies with parametric morphologies to select a sample of featureless disks at 1 < z < 3, which may represent a dynamically warmer progenitor population to the settled disk galaxies seen at later epochs.

Although there has been much progress in understanding how galaxies evolve, we still do not understand how and when they stop forming stars and become quiescent. We address this by applying our galaxy spectral energy distribution models, which incorporate physically motivated star formation histories (SFHs) from cosmological simulations, to a sample of quiescent galaxies at $0.2<z<2.1$. A total of 845 quiescent galaxies with multi-band photometry spanning rest-frame ultraviolet through near-infrared wavelengths are selected from the CANDELS dataset. We compute median SFHs of these galaxies in bins of stellar mass and redshift. At all redshifts and stellar masses, the median SFHs rise, reach a peak, and then decline to reach quiescence. At high redshift, we find that the rise and decline are fast, as expected because the Universe is young. At low redshift, the duration of these phases depends strongly on stellar mass. Low-mass galaxies ($\log(M_{\ast}/M_{\odot})\sim9.5$) grow on average slowly, take a long time to reach their peak of star formation ($\gtrsim 4$ Gyr), and the declining phase is fast ($\lesssim 2$ Gyr). Conversely, high-mass galaxies ($\log(M_{\ast}/M_{\odot})\sim11$) grow on average fast ($\lesssim 2$ Gyr), and, after reaching their peak, decrease the star formation slowly ($\gtrsim 3$ Gyr). These findings are consistent with galaxy stellar mass being a driving factor in determining how evolved galaxies are, with high-mass galaxies being the most evolved at any time (i.e., downsizing). The different durations we observe in the declining phases also suggest that low- and high-mass galaxies experience different quenching mechanisms that operate on different timescales.

[abridged] We quantify the evolution of the stellar mass functions of star-forming and quiescent galaxies as a function of morphology from $z\sim 3$ to the present. Our sample consists of ~50,000 galaxies in the CANDELS fields ($\sim880$ $arcmin^2$), which we divide into four main morphological types, i.e. pure bulge dominated systems, pure spiral disk dominated, intermediate 2-component bulge+disk systems and irregular disturbed galaxies. Our main results are: Star-formation: At $z\sim 2$, 80\% of the stellar mass density of star-forming galaxies is in irregular systems. However, by $z\sim 0.5$, irregular objects only dominate at stellar masses below $10^9M\odot$. A majority of the star-forming irregulars present at $z\sim 2$ undergo a gradual transformation from disturbed to normal spiral disk morphologies by $z\sim 1$ without significant interruption to their star-formation. Rejuvenation after a quenching event does not seem to be common except perhaps for the most massive objects. Quenching: We confirm that galaxies reaching a stellar mass of $M_*\sim10^{10.8}M_\odot$ ($M^*$) tend to quench. Also, quenching implies the presence of a bulge: the abundance of massive red disks is negligible at all redshifts over 2~dex in stellar mass. However the dominant quenching mechanism evolves. At $z>2$, the SMF of quiescent galaxies above $M^*$ is dominated by compact spheroids. Quenching at this early epoch destroys the disk and produces a compact remnant unless the star-forming progenitors at even higher redshifts are significantly more dense. At $1<z<2$, the majority of newly quenched galaxies are disks with a significant central bulge. This suggests that mass-quenching at this epoch starts from the inner parts and preserves the disk. At $z<1$, the high mass end of the passive SMF is globally in place and the evolution mostly happens at stellar masses below $10^{10}M_\odot$.

The Frontier Fields are a director's discretionary time campaign with HST and the Spitzer Space Telescope to see deeper into the universe than ever before. The Frontier Fields combine the power of HST and Spitzer with the natural gravitational telescopes of massive high-magnification clusters of galaxies to produce the deepest observations of clusters and their lensed galaxies ever obtained. Six clusters - Abell 2744, MACSJ0416.1-2403, MACSJ0717.5+3745, MACSJ1149.5+2223, Abell S1063, and Abell 370 - were selected based on their lensing strength, sky darkness, Galactic extinction, parallel field suitability, accessibility to ground-based facilities, HST, Spitzer and JWST observability, and pre-existing ancillary data. These clusters have been targeted by the HST ACS/WFC and WFC3/IR with coordinated parallels of adjacent blank fields for over 840 HST orbits. The Spitzer Space Telescope has dedicated > 1000 hours of director's discretionary time to obtain IRAC 3.6 and 4.5 micron imaging to ~26.5, 26.0 ABmag 5-sigma point-source depths in the six cluster and six parallel Frontier Fields. The Frontier Field parallel fields are the second-deepest observations thus far by HST with ~29th ABmag 5-sigma point source depths in seven optical - near-infrared bandpasses. Galaxies behind the Frontier Field cluster lenses experience typical magnification factors of a few, with small regions near the critical curves magnified by factors 10-100. Therefore, the Frontier Field cluster HST images achieve intrinsic depths of ~30-33 magnitudes over very small volumes. Early studies of the Frontier Fields have probed galaxies fainter than any seen before during the epoch of reionization 6 < z < 10, mapped out the cluster dark matter to unprecedented resolution, and followed lensed transient events.

Current observational evidence suggests that the star formation rate (SFR) efficiency of neutral atomic hydrogen gas measured in Damped Ly-alpha Systems (DLAs) at z~3 is more than 10 times lower than predicted by the Kennicutt-Schmidt (KS) relation. To understand the origin of this deficit, and to investigate possible evolution with redshift and galaxy properties, we measure the SFR efficiency of atomic gas at z~1, z~2, and z~3 around star-forming galaxies. We use new robust photometric redshifts in the Hubble Ultra Deep Field to create galaxy stacks in these three redshift bins, and measure the SFR efficiency by combining DLA absorber statistics with the observed rest-frame UV emission in the galaxies' outskirts. We find that the SFR efficiency of HI gas at z>1 is ~1-3% of that predicted by the KS relation. Contrary to simulations and models that predict a reduced SFR efficiency with decreasing metallicity and thus with increasing redshift, we find no significant evolution in the SFR efficiency with redshift. Our analysis instead suggests that the reduced SFR efficiency is driven by the low molecular content of this atomic-dominated phase, with metallicity playing a secondary effect in regulating the conversion between atomic and molecular gas. This interpretation is supported by the similarity between the observed SFR efficiency and that observed in local atomic-dominated gas, such as in the outskirts of local spiral galaxies and local dwarf galaxies.

We present estimates of intrinsic scatter in the Star Formation Rate (SFR) - Stellar Mass (M*) correlation in the redshift range 0.5 < z < 3.0 and in the mass range 10^7 < M* < 10^11 Msun. We utilize photometry in the Hubble Ultradeep Field (HUDF12), Ultraviolet Ultra Deep Field (UVUDF) campaigns and CANDELS/GOODS-S. We estimate SFR, M* from broadband Spectral Energy Distributions (SEDs) and the best available redshifts. The maximum depth of the HUDF photometry (F160W 29.9 AB, 5 sigma depth) probes the SFR-M* correlation down to M* ~ 10 ^7 Msun, a factor of 10-100X lower in M* than previous studies, and comparable to dwarf galaxies in the local universe. We find the slope of the SFR-M* relationship to be near unity at all redshifts and the normalization to decrease with cosmic time. We find a moderate increase in intrinsic scatter with cosmic time from 0.2 to 0.4 dex across the epoch of peak cosmic star formation. None of our redshift bins show a statistically significant increase in intrinsic scatter at low mass. However, it remains possible that intrinsic scatter increases at low mass on timescales shorter than ~ 100 Myr. Our results are consistent with a picture of gradual and self-similar assembly of galaxies across more than three orders of magnitude in stellar mass from as low as 10^7 Msun.

Important but rare and subtle processes driving galaxy morphology and star-formation may be missed by traditional spiral, elliptical, irregular or S\'ersic bulge/disk classifications. To overcome this limitation, we use a principal component analysis of non-parametric morphological indicators (concentration, asymmetry, Gini coefficient, $M_{20}$, multi-mode, intensity and deviation) measured at rest-frame $B$-band (corresponding to HST/WFC3 F125W at 1.4 $< z <$ 2) to trace the natural distribution of massive ($>10^{10} M_{\odot}$) galaxy morphologies. Principal component analysis (PCA) quantifies the correlations between these morphological indicators and determines the relative importance of each. The first three principal components (PCs) capture $\sim$75 per cent of the variance inherent to our sample. We interpret the first principal component (PC) as bulge strength, the second PC as dominated by concentration and the third PC as dominated by asymmetry. Both PC1 and PC2 correlate with the visual appearance of a central bulge and predict galaxy quiescence. PC1 is a better predictor of quenching than stellar mass, as as good as other structural indicators (S\'ersic-n or compactness). We divide the PCA results into groups using an agglomerative hierarchical clustering method. Unlike S\'ersic, this classification scheme separates compact galaxies from larger, smooth proto-elliptical systems, and star-forming disk-dominated clumpy galaxies from star-forming bulge-dominated asymmetric galaxies. Distinguishing between these galaxy structural types in a quantitative manner is an important step towards understanding the connections between morphology, galaxy assembly and star-formation.

We present the results of a study investigating the sizes and morphologies of redshift 4 < z < 8 galaxies in the CANDELS GOODS-S, HUDF and HUDF parallel fields. Based on non-parametric measurements and incorporating a careful treatment of measurement biases, we quantify the typical size of galaxies at each redshift as the peak of the log-normal size distribution, rather than the arithmetic mean size. Parameterizing the evolution of galaxy half-light radius as $r_{50} \propto (1+z)^n$, we find $n = -0.20 \pm 0.26$ at bright UV-luminosities ($0.3L_{*(z=3)} < L < L_*$) and $n = -0.47 \pm 0.62$ at faint luminosities ($0.12L_* < L < 0.3L_*$). Furthermore, simulations based on artificially redshifting our z~4 galaxy sample show that we cannot reject the null hypothesis of no size evolution. We show that this result is caused by a combination of the size-dependent completeness of high-redshift galaxy samples and the underestimation of the sizes of the largest galaxies at a given epoch. To explore the evolution of galaxy morphology we first compare asymmetry measurements to those from a large sample of simulated single S\'ersic profiles, in order to robustly categorise galaxies as either `smooth' or `disturbed'. Comparing the disturbed fraction amongst bright ($M_{UV} \leq -20$) galaxies at each redshift to that obtained by artificially redshifting our z~4 galaxy sample, while carefully matching the size and UV-luminosity distributions, we find no clear evidence for evolution in galaxy morphology over the redshift interval 4 < z < 8. Therefore, based on our results, a bright ($M_{UV} \leq -20$) galaxy at z~6 is no more likely to be measured as `disturbed' than a comparable galaxy at z~4, given the current observational constraints.

We study the relationships between galaxy environments and galaxy properties related to disk (re)growth, considering two highly complete samples that are approximately baryonic mass limited into the high-mass dwarf galaxy regime, the Environmental COntext (ECO) catalog (data release herein) and the B-semester region of the REsolved Spectroscopy Of a Local VolumE (RESOLVE) survey. We quantify galaxy environments using both group identification and smoothed galaxy density field methods. We use by-eye and quantitative morphological classifications plus atomic gas content measurements and estimates. We find that blue early-type (E/S0) galaxies, gas-dominated galaxies, and UV-bright disk host galaxies all become distinctly more common below group halo mass ~10^11.5 Msun, implying that this low group halo mass regime may be a preferred regime for significant disk growth activity. We also find that blue early-type and blue late-type galaxies inhabit environments of similar group halo mass at fixed baryonic mass, consistent with a scenario in which blue early types can regrow late-type disks. In fact, we find that the only significant difference in the typical group halo mass inhabited by different galaxy classes is for satellite galaxies with different colors, where at fixed baryonic mass red early and late types have higher typical group halo masses than blue early and late types. More generally, we argue that the traditional morphology-environment relation (i.e., that denser environments tend to have more early types) can be largely attributed to the morphology-galaxy mass relation for centrals and the color-environment relation for satellites.

We have undertaken an ambitious program to visually classify all galaxies in the five CANDELS fields down to H<24.5 involving the dedicated efforts of 65 individual classifiers. Once completed, we expect to have detailed morphological classifications for over 50,000 galaxies up to z<4 over all the fields. Here, we present our detailed visual classification scheme, which was designed to cover a wide range of CANDELS science goals. This scheme includes the basic Hubble sequence types, but also includes a detailed look at mergers and interactions, the clumpiness of galaxies, $k$-corrections, and a variety of other structural properties. In this paper, we focus on the first field to be completed -- GOODS-S. The wide area coverage spanning the full field includes 7634 galaxies that have been classified by at least three different people. In the deep area of the field, 2534 galaxies have been classified by at least five different people at three different depths. With this paper, we release to the public all of the visual classifications in GOODS-S along with the GUI that we developed to classify galaxies. We find that the level of agreement among classifiers is good and depends on both the galaxy magnitude and the galaxy type, with disks showing the highest level of agreement and irregulars the lowest. A comparison of our classifications with the Sersic index and rest-frame colors shows a clear separation between disk and spheroid populations. Finally, we explore morphological k-corrections between the V-band and H-band observations and find that a small fraction (84 galaxies in total) are classified as being very different between these two bands. These galaxies typically have very clumpy and extended morphology or are very faint in the V-band.

We present the results from a stellar population modeling analysis of a sample of 162 z=4.5, and 14 z=5.7 Lyman alpha emitting galaxies (LAEs) in the Bootes field, using deep Spitzer/IRAC data at 3.6 and 4.5 um from the Spitzer Lyman Alpha Survey, along with Hubble Space Telescope NICMOS and WFC3 imaging at 1.1 and 1.6 um for a subset of the LAEs. This represents one of the largest samples of high-redshift LAEs imaged with Spitzer IRAC. We find that 30/162 (19%) of the z=4.5 LAEs and 9/14 (64%) of the z=5.7 LAEs are detected at >3-sigma in at least one IRAC band. Individual z=4.5 IRAC-detected LAEs have a large range of stellar mass, from 5x10^8 to 10^11 Msol. One-third of the IRAC-detected LAEs have older stellar population ages of 100 Myr - 1 Gyr, while the remainder have ages < 100 Myr. A stacking analysis of IRAC-undetected LAEs shows this population to be primarily low mass (8 -- 20 x 10^8 Msol) and young (64 - 570 Myr). We find a correlation between stellar mass and the dust-corrected ultraviolet-based star-formation rate (SFR) similar to that at lower redshifts, in that higher mass galaxies exhibit higher SFRs. However, the z=4.5 LAE correlation is elevated 4-5 times in SFR compared to continuum-selected galaxies at similar redshifts. The exception is the most massive LAEs which have SFRs similar to galaxies at lower redshifts suggesting that they may represent a different population of galaxies than the traditional lower-mass LAEs, perhaps with a different mechanism promoting Lyman alpha photon escape.

We derive the total cold gas, atomic hydrogen, and molecular gas masses of approximately 24 000 galaxies covering four decades in stellar mass at redshifts 0.5 < z < 3.0, taken from the CANDELS survey. Our inferences are based on the inversion of a molecular hydrogen based star formation law, coupled with a prescription to separate atomic and molecular gas. We find that: 1) there is an increasing trend between the inferred cold gas (HI and H2), HI, and H2 mass and the stellar mass of galaxies down to stellar masses of 10^8 Msun already in place at z = 3; 2) the molecular fractions of cold gas increase with increasing stellar mass and look-back time; 3) there is hardly any evolution in the mean HI content of galaxies at fixed stellar mass; 4) the cold gas fraction and relative amount of molecular hydrogen in galaxies decrease at a relatively constant rate with time, independent of stellar mass; 5) there is a large population of low-stellar mass galaxies dominated by atomic gas. These galaxies are very gas rich, but only a minor fraction of their gas is molecular; 6) the ratio between star-formation rate (SFR) and inferred total cold gas mass (HI + H2) of galaxies (i.e., star-formation efficiency; SFE) increases with star-formation at fixed stellar masses. Due to its simplicity, the presented approach is valuable to assess the impact of selection biases on small samples of directly-observed gas masses and to extend scaling relations down to stellar mass ranges and redshifts that are currently difficult to probe with direct measurements of gas content.

The spatial fluctuations of the extragalactic background light trace the total emission from all stars and galaxies in the Universe. A multi-wavelength study can be used to measure the integrated emission from first galaxies during reionization when the Universe was about 500 million years old. Here we report arcminute-scale spatial fluctuations in one of the deepest sky surveys with the Hubble Space Telescope in five wavebands between 0.6 and 1.6 $\mu$m. We model-fit the angular power spectra of intensity fluctuation measurements to find the ultraviolet luminosity density of galaxies at $z$ > 8 to be $\log \rho_{\rm UV} = 27.4^{+0.2}_{-1.2}$ erg s$^{-1}$ Hz$^{-1}$ Mpc$^{-3}$ $(1\sigma)$. This level of integrated light emission allows for a significant surface density of fainter primeval galaxies that are below the point source detection level in current surveys.

We present a robust measurement and analysis of the rest-frame ultraviolet (UV) luminosity function at z=4-8. We use deep Hubble Space Telescope imaging over the CANDELS/GOODS fields, the Hubble Ultra Deep Field and the Year 1 Hubble Frontier Field deep parallel observations. These surveys provides an effective volume of 0.6-1.2 x 10^6 Mpc^3 over this epoch, allowing us to perform a robust search for faint (M_UV=-18) and bright (M_UV < -21) galaxies. We select candidate galaxies using a well-tested photometric redshift technique with careful screening of contaminants, finding a sample of 7446 galaxies at 3.5<z<8.5, with >1000 galaxies at z~6-8. We measure the luminosity function using a Markov Chain Monte Carlo analysis to measure robust uncertainties. At the faint end our results agree with previous studies, yet we find a higher abundance of UV-bright galaxies at z>6, with M* ~ -21 at z>5, different than that inferred based on previous trends at lower redshift. At z=8, a single power-law provides an equally good fit to the UV luminosity function, while at z=6 and 7, an exponential cutoff at the bright-end is moderately preferred. We compare to semi-analytical models, and find that the lack of evolution in M* is consistent with models where the impact of dust attenuation on the bright-end of the luminosity function decreases at higher redshift. We measure the evolution of the cosmic star-formation rate density, correcting for dust attenuation, and find that it declines as (1+z)^(-4.3 +/- 0.5) at z>4, consistent with observations at z>9. Our observations are consistent with a reionization history that starts at z>10, completes at z>6, and reaches a midpoint (x_HII = 0.5) at 6.7<z<9.4. Finally, our observations predict that the abundance of bright z=9 galaxies is likely higher than previous constraints, though consistent with recent estimates of bright z~10 galaxies. [abridged]

The Spitzer-Cosmic Assembly Deep Near-Infrared Extragalactic Legacy Survey (S-CANDELS; PI G. Fazio) is a Cycle 8 Exploration Program designed to detect galaxies at very high redshifts (z > 5). To mitigate the effects of cosmic variance and also to take advantage of deep coextensive coverage in multiple bands by the Hubble Space Telescope Multi-Cycle Treasury Program CANDELS, S-CANDELS was carried out within five widely separated extragalactic fields: the UKIDSS Ultra-Deep Survey, the Extended Chandra Deep Field South, COSMOS, the HST Deep Field North, and the Extended Groth Strip. S-CANDELS builds upon the existing coverage of these fields from the Spitzer Extended Deep Survey (SEDS) by increasing the integration time from 12 hours to a total of 50 hours but within a smaller area, 0.16 square degrees. The additional depth significantly increases the survey completeness at faint magnitudes. This paper describes the S-CANDELS survey design, processing, and publicly-available data products. We present IRAC dual-band 3.6+4.5 micron catalogs reaching to a depth of 26.5 AB mag. Deep IRAC counts for the roughly 135,000 galaxies detected by S-CANDELS are consistent with models based on known galaxy populations. The increase in depth beyond earlier Spitzer/IRAC surveys does not reveal a significant additional contribution from discrete sources to the diffuse Cosmic Infrared Background (CIB). Thus it remains true that only roughly half of the estimated CIB flux from COBE/DIRBE is resolved.

In this paper we perform a comprehensive study of the main sources of random and systematic errors in stellar mass measurement for galaxies using their Spectral Energy Distributions (SEDs). We use mock galaxy catalogs with simulated multi-waveband photometry (from U-band to mid-infrared) and known redshift, stellar mass, age and extinction for individual galaxies. Given different parameters affecting stellar mass measurement (photometric S/N ratios, SED fitting errors, systematic effects, the inherent degeneracies and correlated errors), we formulated different simulated galaxy catalogs to quantify these effects individually. We studied the sensitivity of stellar mass estimates to the codes/methods used, population synthesis models, star formation histories, nebular emission line contributions, photometric uncertainties, extinction and age. For each simulated galaxy, the difference between the input stellar masses and those estimated using different simulation catalogs, $\Delta\log(M)$, was calculated and used to identify the most fundamental parameters affecting stellar masses. We measured different components of the error budget, with the results listed as follows: (1). no significant bias was found among different codes/methods, with all having comparable scatter; (2). A source of error is found to be due to photometric uncertainties and low resolution in age and extinction grids; (3). The median of stellar masses among different methods provides a stable measure of the mass associated with any given galaxy; (4). The deviations in stellar mass strongly correlate with those in age, with a weaker correlation with extinction; (5). the scatter in the stellar masses due to free parameters are quantified, with the sensitivity of the stellar mass to both the population synthesis codes and inclusion of nebular emission lines studied.

We present photometry and derived redshifts from up to eleven bandpasses for 9927 galaxies in the Hubble Ultra Deep field (UDF), covering an observed wavelength range from the near-ultraviolet (NUV) to the near-infrared (NIR) with Hubble Space Telescope observations. Our Wide Field Camera 3 (WFC3)/UV F225W, F275W, and F336W image mosaics from the ultra-violet UDF (UVUDF) imaging campaign are newly calibrated to correct for charge transfer inefficiency, and use new dark calibrations to minimize background gradients and pattern noise. Our NIR WFC3/IR image mosaics combine the imaging from the UDF09 and UDF12 campaigns with CANDELS data to provide NIR coverage for the entire UDF field of view. We use aperture-matched point-spread function corrected photometry to measure photometric redshifts in the UDF, sampling both the Lyman break and Balmer break of galaxies at z~0.8-3.4, and one of the breaks over the rest of the redshift range. Our comparison of these results with a compilation of robust spectroscopic redshifts shows an improvement in the galaxy photometric redshifts by a factor of two in scatter and a factor three in outlier fraction over previous UDF catalogs. The inclusion of the new NUV data is responsible for a factor of two decrease in the outlier fraction compared to redshifts determined from only the optical and NIR data, and improves the scatter at z<0.5 and at z>2. The panchromatic coverage of the UDF from the NUV through the NIR yields robust photometric redshifts of the UDF, with the lowest outlier fraction available.

We present the discovery of two galaxy overdensities in the Hubble Space Telescope UDF: a proto-cluster, HUDFJ0332.4-2746.6 at $z = 1.84 \pm 0.01$, and a group, HUDFJ0332.5-2747.3 at $z =1.90 \pm 0.01$. Assuming viralization, the velocity dispersion of HUDFJ0332.4-2746.6 implies a mass of $M_{200}= (2.2 \pm 1.8) \times 10^{14} M_{\odot}$, consistent with the lack of extended X-ray emission. Neither overdensity shows evidence of a red sequence. About $50\%$ of their members show interactions and/or disturbed morphologies, which are signatures of merger remnants or disk instability. Most of their ETGs have blue colors and show recent star-formation. These observations reveal for the first time large fractions of spectroscopically confirmed star-forming blue ETGs in proto-clusters at $z\approx 2$. These star-forming ETGs are most likely among the progenitors of the quiescent population in clusters at more recent epochs. Their mass-size relation is consistent with that of passive ETGs in clusters at $z\sim0.7-1.5$. If these galaxies are the progenitors of cluster ETGs at these lower redshifts, their size would evolve according to a similar mass-size relation. It is noteworthy that quiescent ETGs in clusters at $z=1.8-2$ also do not show any significant size evolution over this redshift range, contrary to field ETGs. The ETG fraction is $\lesssim 50\%$, compared to the typical quiescent ETG fraction of $\approx 80\%$ in cluster cores at $z< 1$. The fraction, masses, and colors of the newly discovered ETGs imply that other cluster ETGs will be formed/accreted at a later time.

We have conducted a detailed object-by-object study of a mass-complete (M*>10^11 M_sun) sample of 56 galaxies at 1.4 < z < 2 in the GOODS-South field, showing that an accurate de-blending in MIPS/24um images is essential to properly assign to each galaxy its own star formation rate (SFR), whereas an automatic procedure often fails. This applies especially to galaxies with SFRs below the Main Sequence (MS) value, which may be in their quenching phase. After that, the sample splits evenly between galaxies forming stars within a factor of 4 of the MS rate (~45%), and sub-MS galaxies with SFRs ~10-1000 times smaller (~55%). We did not find a well defined class of intermediate, transient objects below the MS, suggesting that the conversion of a massive MS galaxy into a quenched remnant may take a relatively short time (<1 Gyr), though a larger sample should be analyzed in the same way to set precise limits on the quenching timescale. X-ray detected AGNs represent a ~30% fraction of the sample, and are found among both star-forming and quenched galaxies. The morphological analysis revealed that ~50% of our massive objects are bulge-dominated, and almost all MS galaxies with a relevant bulge component host an AGN. We also found sub-MS SFRs in many bulge-dominated systems, providing support to the notion that bulge growth, AGN activity and quenching of star formation are closely related to each other.

We use data from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey to study how the spatial variation in the stellar populations of galaxies relate to the formation of galaxies at $1.5 < z < 3.5$. We use the Internal Color Dispersion (ICD), measured between the rest-frame UV and optical bands, which is sensitive to age (and dust attenuation) variations in stellar populations. The ICD shows a relation with the stellar masses and morphologies of the galaxies. Galaxies with the largest variation in their stellar populations as evidenced by high ICD have disk-dominated morphologies (with S\'{e}rsic indexes $< 2$) and stellar masses between $10 < \mathrm{Log~M/ M_\odot}< 11$. There is a marked decrease in the ICD as the stellar mass and/or the S\'ersic index increases. By studying the relations between the ICD and other galaxy properties including sizes, total colors, star-formation rate, and dust attenuation, we conclude that the largest variations in stellar populations occur in galaxies where the light from newly, high star-forming clumps contrasts older stellar disk populations. This phase reaches a peak for galaxies only with a specific stellar mass range, $10 < \mathrm{Log~M/ M_\odot} < 11$, and prior to the formation of a substantial bulge/spheroid. In contrast, galaxies at higher or lower stellar masses, and/or higher S\'{e}rsic index ($n > 2$) show reduced ICD values, implying a greater homogeneity of their stellar populations. This indicates that if a galaxy is to have both a quiescent bulge along with a star forming disk, typical of Hubble Sequence galaxies, this is most common for stellar masses $10 < \mathrm{Log~M/M_\odot} < 11$ and when the bulge component remains relatively small ($n<2$).

We present the public release of the stellar mass catalogs for the GOODS-S and UDS fields obtained using some of the deepest near-IR images available, achieved as part of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) project. We combine the effort from ten different teams, who computed the stellar masses using the same photometry and the same redshifts. Each team adopted their preferred fitting code, assumptions, priors, and parameter grid. The combination of results using the same underlying stellar isochrones reduces the systematics associated with the fitting code and other choices. Thanks to the availability of different estimates, we can test the effect of some specific parameters and assumptions on the stellar mass estimate. The choice of the stellar isochrone library turns out to have the largest effect on the galaxy stellar mass estimates, resulting in the largest distributions around the median value (with a semi interquartile range larger than 0.1 dex). On the other hand, for most galaxies, the stellar mass estimates are relatively insensitive to the different parameterizations of the star formation history. The inclusion of nebular emission in the model spectra does not have a significant impact for the majority of galaxies (less than a factor of 2 for ~80% of the sample). Nevertheless, the stellar mass for the subsample of young galaxies (age < 100 Myr), especially in particular redshift ranges (e.g., 2.2 < z < 2.4, 3.2 < z < 3.6, and 5.5 < z < 6.5), can be seriously overestimated (by up to a factor of 10 for < 20 Myr sources) if nebular contribution is ignored.

In order to derive the AGN contribution to the cosmological ionizing emissivity we have selected faint AGN candidates at $z>4$ in the CANDELS GOODS-South field which is one of the deepest fields with extensive multiwavelength coverage from Chandra, HST, Spitzer and various groundbased telescopes. We have adopted a relatively novel criterion. As a first step high redshift galaxies are selected in the NIR $H$ band down to very faint levels ($H\leq27$) using reliable photometric redshifts. This corresponds at $z>4$ to a selection criterion based on the galaxy rest-frame UV flux. AGN candidates are then picked up from this parent sample if they show X-ray fluxes above a threshold of $F_X\sim 1.5\times 10^{-17}$ cgs (0.5-2 keV). We have found 22 AGN candidates at $z>4$ and we have derived the first estimate of the UV luminosity function in the redshift interval $4<z<6.5$ and absolute magnitude interval $-22.5\lesssim M_{1450} \lesssim -18.5$ typical of local Seyfert galaxies. The faint end of the derived luminosity function is about two/four magnitudes fainter at $z\sim 4-6$ than that derived from previous UV surveys. We have then estimated ionizing emissivities and hydrogen photoionization rates in the same redshift interval under reasonable assumptions and after discussion of possible caveats, the most important being the large uncertainties involved in the estimate of photometric redshift for sources with featureless, almost power-law SEDs and/or low average escape fraction of ionizing photons from the AGN host galaxies. We argue that, under reasonable evaluations of possible biases, the probed AGN population can produce at $z=4-6.5$ photoionization rates consistent with that required to keep highly ionized the intergalactic medium observed in the Lyman-$\alpha$ forest of high redshift QSO spectra, providing an important contribution to the cosmic reionization.

We combine HST/WFC3 imaging and G141 grism observations from the CANDELS and 3D-HST surveys to produce a catalog of grism spectroscopic redshifts for galaxies in the CANDELS/GOODS-South field. The WFC3/G141 grism spectra cover a wavelength range of 1.1<lambda<1.7 microns with a resolving power of R~130 for point sources, thus providing rest-frame optical spectra for galaxies out to z~3.5. The catalog is selected in the H-band (F160W) and includes both galaxies with and without previously published spectroscopic redshifts. Grism spectra are extracted for all H-band detected galaxies with H<24 and a CANDELS photometric redshift z_phot > 0.6. The resulting spectra are visually inspected to identify emission lines and redshifts are determined using cross-correlation with empirical spectral templates. To establish the accuracy of our redshifts, we compare our results against high-quality spectroscopic redshifts from the literature. Using a sample of 411 control galaxies, this analysis yields a precision of sigma_NMAD=0.0028 for the grism-derived redshifts, which is consistent with the accuracy reported by the 3D-HST team. Our final catalog covers an area of 153 square arcmin and contains 1019 redshifts for galaxies in GOODS-S. Roughly 60% (608/1019) of these redshifts are for galaxies with no previously published spectroscopic redshift. These new redshifts span a range of 0.677 < z < 3.456 and have a median redshift of z=1.282. The catalog contains a total of 234 new redshifts for galaxies at z>1.5. In addition, we present 20 galaxy pair candidates identified for the first time using the grism redshifts in our catalog, including four new galaxy pairs at z~2, nearly doubling the number of such pairs previously identified.

Quenched galaxies at z>2 are nearly all very compact relative to z~0, suggesting a physical connection between high stellar density and efficient, rapid cessation of star-formation. We present restframe UV spectra of Lyman-break galaxies (LBGs) at z~3 selected to be candidate progenitors of quenched galaxies at z~2 based on their compact restframe optical sizes and high surface density of star-formation. We compare their UV properties to those of more extended LBGs of similar mass and star formation rate (non-candidates). We find that candidate progenitors have faster ISM gas velocities and higher equivalent widths of interstellar absorption lines, implying larger velocity spread among absorbing clouds. Candidates deviate from the relationship between equivalent widths of Lyman-alpha and interstellar absorption lines in that their Lyman-alpha emission remains strong despite high interstellar absorption, possibly indicating that the neutral HI fraction is patchy such that Lyman-alpha photons can escape. We detect stronger CIV P-Cygni features (emission and absorption) and HeII emission in candidates, indicative of larger populations of metal rich Wolf-Rayet stars compared to non-candidates. The faster bulk motions, broader spread of gas velocity, and Lyman-alpha properties of candidates are consistent with their ISM being subject to more energetic feedback than non-candidates. Together with their larger metallicity (implying more evolved star-formation activity) this leads us to propose, if speculatively, that they are likely to quench sooner than non-candidates, supporting the validity of selection criteria used to identify them as progenitors of z~2 passive galaxies. We propose that massive, compact galaxies undergo more rapid growth of stellar mass content, perhaps because the gas accretion mechanisms are different, and quench sooner than normally-sized LBGs at these early epochs.

The galaxy stellar mass function (GSMF) at high-z provides key information on star-formation history and mass assembly in the young Universe. We aimed to use the unique combination of deep optical/NIR/MIR imaging provided by HST, Spitzer and the VLT in the CANDELS-UDS, GOODS-South, and HUDF fields to determine the GSMF over the redshift range 3.5<z<7.5. We utilised the HST WFC3/IR NIR imaging from CANDELS and HUDF09, reaching H~27-28.5 over a total area of 369 arcmin2, in combination with associated deep HST ACS optical data, deep Spitzer IRAC imaging from the SEDS programme, and deep Y and K-band VLT Hawk-I images from the HUGS programme, to select a galaxy sample with high-quality photometric redshifts. These have been calibrated with more than 150 spectroscopic redshifts in the range 3.5<z<7.5, resulting in an overall precision of sigma_z/(1+z)~0.037. We have determined the low-mass end of the high-z GSMF with unprecedented precision, reaching down to masses as low as M*~10^9 Msun at z=4 and ~6x10^9 Msun at z=7. We find that the GSMF at 3.5<z<7.5 depends only slightly on the recipes adopted to measure the stellar masses, namely the photo-z, the SFHs, the nebular contribution or the presence of AGN on the parent sample. The low-mass end of the GSMF is steeper than has been found at lower redshifts, but appears to be unchanged over the redshift range probed here. Our results are very different from previous GSMF estimates based on converting UV galaxy luminosity functions into mass functions via tight M/L relations. Integrating our evolving GSMF over mass, we find that the growth of stellar mass density is barely consistent with the time-integral of the SFR density over cosmic time at z>4. These results confirm the unique synergy of the CANDELS+HUDF, HUGS, and SEDS surveys for the discovery and study of moderate/low-mass galaxies at high redshifts.

We combine photometry from the UDS, and CANDELS UDS and CANDELS GOODS-S surveys to construct the galaxy stellar mass function probing both the low and high mass end accurately in the redshift range 0.3<z<3. The advantages of using a homogeneous concatenation of these datasets include meaningful measures of environment in the UDS, due to its large area (0.88 deg^2), and the high resolution deep imaging in CANDELS (H_160 > 26.0), affording us robust measures of structural parameters. We construct stellar mass functions for the entire sample as parameterised by the Schechter function, and find that there is a decline in the values of phi and of alpha with higher redshifts, and a nearly constant M* up to z~3. We divide the galaxy stellar mass function by colour, structure, and environment and explore the links between environmental over-density, morphology, and the quenching of star formation. We find that a double Schechter function describes galaxies with high Sersic index (n>2.5), similar to galaxies which are red or passive. The low-mass end of the n>2.5 stellar mass function is dominated by blue galaxies, whereas the high-mass end is dominated by red galaxies. This hints that possible links between morphological evolution and star formation quenching are only present in high-mass galaxies. This is turn suggests that there are strong mass dependent quenching mechanisms. In addition, we find that the number density of high mass systems is elevated in dense environments, suggesting that an environmental process is building up massive galaxies quicker in over densities than in lower densities.

We present photometric redshifts and associated probability distributions for all detected sources in the Extended Chandra Deep Field South (ECDFS). The work makes use of the most up-to-date data from the Cosmic Assembly Near-IR Deep Legacy Survey (CANDELS) and the Taiwan ECDFS Near-Infrared Survey (TENIS) in addition to other data. We also revisit multi-wavelength counterparts for published X-ray sources from the 4Ms-CDFS and 250ks-ECDFS surveys, finding reliable counterparts for 1207 out of 1259 sources ($\sim 96\%$). Data used for photometric redshifts include intermediate-band photometry deblended using the TFIT method, which is used for the first time in this work. Photometric redshifts for X-ray source counterparts are based on a new library of AGN/galaxy hybrid templates appropriate for the faint X-ray population in the CDFS. Photometric redshift accuracy for normal galaxies is 0.010 and for X-ray sources is 0.014, and outlier fractions are $4\%$ and $5.4\%$ respectively. The results within the CANDELS coverage area are even better as demonstrated both by spectroscopic comparison and by galaxy-pair statistics. Intermediate-band photometry, even if shallow, is valuable when combined with deep broad-band photometry. For best accuracy, templates must include emission lines.

We present the results of a new, ultra-deep, near-infrared imaging survey executed with the Hawk-I imager at the ESO VLT, of which we make all the data public. This survey, named HUGS (Hawk-I UDS and GOODS Survey), provides deep, high-quality imaging in the K and Y bands over the CANDELS UDS and GOODS-South fields. We describe here the survey strategy, the data reduction process, and the data quality. HUGS delivers the deepest and highest quality K-band images ever collected over areas of cosmological interest, and ideally complements the CANDELS data set in terms of image quality and depth. The seeing is exceptional and homogeneous, confined to the range 0.38"-0.43". In the deepest region of the GOODS-S field, (which includes most of the HUDF) the K-band exposure time exceeds 80 hours of integration, yielding a 1-sigma magnitude limit of ~28.0 mag/sqarcsec. In the UDS field the survey matches the shallower depth of the CANDELS images reaching a 1-sigma limit per sq.arcsec of ~27.3mag in the K band and ~28.3mag in the Y-band, We show that the HUGS observations are well matched to the depth of the CANDELS WFC3/IR data, since the majority of even the faintest galaxies detected in the CANDELS H-band images are also detected in HUGS. We present the K-band galaxy number counts produced by combining the HUGS data from the two fields. We show that the slope of the number counts depends sensitively on the assumed distribution of galaxy sizes, with potential impact on the estimated extra-galactic background light (abridged).

We study the relationship between the structure and star-formation rate (SFR) of X-ray selected low and moderate luminosity active galactic nuclei (AGNs) in the two Chandra Deep Fields, using Hubble Space Telescope imaging from the Cosmic Assembly Near Infrared Extragalactic Legacy Survey (CANDELS) and deep far-infrared maps from the PEP+GOODS-Herschel survey. We derive detailed distributions of structural parameters and FIR luminosities from carefully constructed control samples of galaxies, which we then compare to those of the AGNs. At z~1, AGNs show slightly diskier light profiles than massive inactive (non-AGN) galaxies, as well as modestly higher levels of gross galaxy disturbance (as measured by visual signatures of interactions and clumpy structure). In contrast, at z~2, AGNs show similar levels of galaxy disturbance as inactive galaxies, but display a red central light enhancement, which may arise due to a more pronounced bulge in AGN hosts or due to extinguished nuclear light. We undertake a number of tests of these alternatives, but our results do not strongly favour one interpretation over the other. The mean SFR and its distribution among AGNs and inactive galaxies are similar at z>1.5. At z<1, however, clear and significant enhancements are seen in the SFRs of AGNs with bulge-dominated light profiles. These trends suggest an evolution in the relation between nuclear activity and host properties with redshift, towards a minor role for mergers and interactions at z>1.5.

The formation of bars in disk galaxies is a tracer of the dynamical maturity of the population. Previous studies have found that the incidence of bars in disks decreases from the local Universe to z ~ 1, and by z > 1 simulations predict that bar features in dynamically mature disks should be extremely rare. Here we report the discovery of strong barred structures in massive disk galaxies at z ~ 1.5 in deep rest-frame optical images from CANDELS. From within a sample of 876 disk galaxies identified by visual classification in Galaxy Zoo, we identify 123 barred galaxies. Selecting a sub-sample within the same region of the evolving galaxy luminosity function (brighter than L*), we find that the bar fraction across the redshift range 0.5< z < 2 (f_bar = 10.7 +6.3 -3.5% after correcting for incompleteness) does not significantly evolve. We discuss the implications of this discovery in the context of existing simulations and our current understanding of the way disk galaxies have evolved over the last 11 billion years.

We measure new estimates for the galaxy stellar mass function and star formation rates for samples of galaxies at $z \sim 4,~5,~6~\&~7$ using data in the CANDELS GOODS South field. The deep near-infrared observations allow us to construct the stellar mass function at $z \geq 6$ directly for the first time. We estimate stellar masses for our sample by fitting the observed spectral energy distributions with synthetic stellar populations, including nebular line and continuum emission. The observed UV luminosity functions for the samples are consistent with previous observations, however we find that the observed $M_{UV}$ - M$_{*}$ relation has a shallow slope more consistent with a constant mass to light ratio and a normalisation which evolves with redshift. Our stellar mass functions have steep low-mass slopes ($\alpha \approx -1.9$), steeper than previously observed at these redshifts and closer to that of the UV luminosity function. Integrating our new mass functions, we find the observed stellar mass density evolves from $\log_{10} \rho_{*} = 6.64^{+0.58}_{-0.89}$ at $z \sim 7$ to $7.36\pm0.06$ $\text{M}_{\odot} \text{Mpc}^{-3}$ at $z \sim 4$. Finally, combining the measured UV continuum slopes ($\beta$) with their rest-frame UV luminosities, we calculate dust corrected star-formation rates (SFR) for our sample. We find the specific star-formation rate for a fixed stellar mass increases with redshift whilst the global SFR density falls rapidly over this period. Our new SFR density estimates are higher than previously observed at this redshift.

We estimate the UV continuum slope, beta, for 923 galaxies in the range 1 < z < 8 in the Hubble Ultradeep Field (HUDF). These data include 460 galaxies at 1 < z < 2 down to an absolute magnitude M_{UV} = -14 (~0.006 L*_{z=1}; 0.02 L*_{z=0}), comparable to dwarf galaxies in the local universe. We combine deep HST/UVIS photometry in F225W, F275W, F336W wavebands (UVUDF) with recent data from HST/WFC3/IR (HUDF12). Galaxies in the range 1 < z < 2 are significantly bluer than local dwarf galaxies. We find their mean (median) values <beta> = -1.382 (-1.830) +/- 0.002 (random) +/- 0.1 (systematic). We find comparable scatter in beta (standard deviation = 0.43) to local dwarf galaxies and 30% larger scatter than z > 2 galaxies. We study the trends of beta with redshift and absolute magnitude for binned sub-samples and find a modest color-magnitude relation, dbeta/dM = -0.11 +/- 0.01 and no evolution in dbeta/dM with redshift. A modest increase in dust reddening with redshift and luminosity, Delta E(B-V) ~ 0.1, and a comparable increase in the dispersion of dust reddening at z < 2, appears likely to explain the observed trends. At z > 2, we find trends that are consistent with previous works; combining our data with the literature in the range 1 < z < 8, we find a color evolution with redshift, dbeta/dz = -0.09 +/-0.01 for low luminosity (0.05 L*_{z=3}), and dbeta/dz = -0.06 +/-0.01 for medium luminosity (0.25 L*_{z=3}) galaxies.

We use CANDELS imaging, 3D-HST spectroscopy, and Chandra X-ray data to investigate if active galactic nuclei (AGNs) are preferentially fueled by violent disk instabilities funneling gas into galaxy centers at 1.3<z<2.4. We select galaxies undergoing gravitational instabilities using the number of clumps and degree of patchiness as proxies. The CANDELS visual classification system is used to identify 44 clumpy disk galaxies, along with mass-matched comparison samples of smooth and intermediate morphology galaxies. We note that, despite being being mass-matched and having similar star formation rates, the smoother galaxies tend to be smaller disks with more prominent bulges compared to the clumpy galaxies. The lack of smooth extended disks is probably a general feature of the z~2 galaxy population, and means we cannot directly compare with the clumpy and smooth extended disks observed at lower redshift. We find that z~2 clumpy galaxies have slightly enhanced AGN fractions selected by integrated line ratios (in the mass-excitation method), but the spatially resolved line ratios indicate this is likely due to extended phenomena rather than nuclear AGNs. Meanwhile the X-ray data show that clumpy, smooth, and intermediate galaxies have nearly indistinguishable AGN fractions derived from both individual detections and stacked non-detections. The data demonstrate that AGN fueling modes at z~1.85 - whether violent disk instabilities or secular processes - are as efficient in smooth galaxies as they are in clumpy galaxies.

Spectroscopic + photometric redshifts, stellar mass estimates, and rest-frame colors from the 3D-HST survey are combined with structural parameter measurements from CANDELS imaging to determine the galaxy size-mass distribution over the redshift range 0<z<3. Separating early- and late-type galaxies on the basis of star-formation activity, we confirm that early-type galaxies are on average smaller than late-type galaxies at all redshifts, and find a significantly different rate of average size evolution at fixed galaxy mass, with fast evolution for the early-type population, R_eff ~ (1+z)^-1.48, and moderate evolution for the late-type population, R_eff ~ (1+z)^-0.75. The large sample size and dynamic range in both galaxy mass and redshift, in combination with the high fidelity of our measurements due to the extensive use of spectroscopic data, not only fortify previous results, but also enable us to probe beyond simple average galaxy size measurements. At all redshifts the slope of the size-mass relation is shallow, R_eff ~ M_star^0.22, for late-type galaxies with stellar mass >3x10^9 M_sol, and steep, R_eff M_star^0.75, for early-type galaxies with stellar mass >2x10^10 M_sol. The intrinsic scatter is <~0.2 dex for all galaxy types and redshifts. For late-type galaxies, the logarithmic size distribution is not symmetric, but skewed toward small sizes: at all redshifts and masses a tail of small late-type galaxies exists that overlaps in size with the early-type galaxy population. The number density of massive (~10^11 M_sol), compact (R_eff < 2 kpc) early-type galaxies increases from z=3 to z=1.5-2 and then strongly decreases at later cosmic times.

We have constructed a mass-selected sample of Mstar>10^11Msolar galaxies at 1<z<3 in the CANDELS UDS and COSMOS fields and have decomposed these systems into their separate bulge and disk components according to their H(160)-band morphologies. By extending this analysis to multiple bands we have been able to conduct individual bulge and disk component SED fitting which has provided us with stellar-mass and star-formation rate estimates for the separate bulge and disk components. These have been combined with size measurements to explore the evolution of these massive high-redshift galaxies. By utilising the new decomposed stellar-mass estimates, we confirm that the bulge components display a stronger size evolution than the disks. This can be seen from both the fraction of bulge components which lie below the local relation and the median sizes of the bulge components, where the bulges are a median factor of 2.93+/-0.32 times smaller than similarly massive local galaxies at 1<z<2 and 3.41+/-0.58 smaller at 2<z<3; for the disks the corresponding factors are 1.65+/-0.14 and 1.99+/-0.25. Moreover, by splitting our sample into the passive and star-forming bulge and disk sub-populations and examining their sizes as a fraction of their present-day counter-parts, we find that the star-forming and passive bulges are equally compact, star-forming disks are larger, while the passive disks have intermediate sizes. This trend is not evident when classifying galaxy morphology on the basis of single-Sersic fits and adopting the overall star-formation rates. Finally, by evolving the star-formation histories of the passive disks back to the redshifts when the passive disks were last active, we show that the passive and star-forming disks have consistent sizes at the relevant epoch. These trends need to be reproduced by any mechanisms which attempt to explain the morphological evolution of galaxies.

We present the results of a new and improved study of the morphological and spectral evolution of massive galaxies over the redshift range 1<z<3. Our analysis is based on a bulge-disk decomposition of 396 galaxies with Mstar>10^11 Msolar from the CANDELS WFC3/IR imaging within the COSMOS and UKIDSS UDS survey fields. We find that, by modelling the H(160) image of each galaxy with a combination of a de Vaucouleurs bulge (Sersic index n=4) and an exponential disk (n=1), we can then lock all derived morphological parameters for the bulge and disk components, and successfully reproduce the shorter-wavelength J(125), i(814), v(606) HST images simply by floating the magnitudes of the two components. This then yields sub-divided 4-band HST photometry for the bulge and disk components which, with no additional priors, is well described by spectrophotometric models of galaxy evolution. Armed with this information we are able to properly determine the masses and star-formation rates for the bulge and disk components, and find that: i) from z=3 to z=1 the galaxies move from disk-dominated to increasingly bulge-dominated, but very few galaxies are pure bulges/ellipticals by z=1; ii) while most passive galaxies are bulge-dominated, and most star-forming galaxies disk-dominated, 18+/-5% of passive galaxies are disk-dominated, and 11+/-3% of star-forming galaxies are bulge-dominated, a result which needs to be explained by any model purporting to connect star-formation quenching with morphological transformations; iii) there exists a small but significant population of pure passive disks, which are generally flatter than their star-forming counterparts (whose axial ratio distribution peaks at b/a~0.7); iv) flatter/larger disks re-emerge at the highest star-formation rates, consistent with recent studies of sub-mm galaxies, and with the concept of a maximum surface-density for star-formation activity.

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) was a multi-cycle treasury program on the Hubble Space Telescope (HST) that surveyed a total area of ~0.25 deg^2 with ~900 HST orbits spread across 5 fields over 3 years. Within these survey images we discovered 65 supernovae (SN) of all types, out to z~2.5. We classify ~24 of these as Type Ia SN (SN Ia) based on host-galaxy redshifts and SN photometry (supplemented by grism spectroscopy of 6 SN). Here we present a measurement of the volumetric SN Ia rate as a function of redshift, reaching for the first time beyond z=2 and putting new constraints on SN Ia progenitor models. Our highest redshift bin includes detections of SN that exploded when the universe was only ~3 Gyr old and near the peak of the cosmic star-formation history. This gives the CANDELS high-redshift sample unique leverage for evaluating the fraction of SN Ia that explode promptly after formation (<500 Myr). Combining the CANDELS rates with all available SN Ia rate measurements in the literature we find that this prompt SN Ia fraction is fP=0.53 +0.09 -0.10 (stat) +0.10 -0.26 (sys), consistent with a delay time distribution that follows a simple t^{-1} power law for all times t>40 Myr. However, a mild tension is apparent between ground-based low-z surveys and space-based high-z surveys. In both CANDELS and the sister HST program CLASH, we find a low rate of SN Ia at z>1. This could be a hint that prompt progenitors are in fact relatively rare, accounting for only ~20% of all SN Ia explosions -- though further analysis and larger samples will be needed to examine that suggestion.

Exploiting the deep high-resolution imaging of all 5 CANDELS fields, and accurate redshift information provided by 3D-HST, we investigate the relation between structure and stellar populations for a mass-selected sample of 6764 galaxies above 10^10 Msun, spanning the redshift range 0.5 < z < 2.5. For the first time, we fit 2-dimensional models comprising a single Sersic fit and two-component (i.e., bulge + disk) decompositions not only to the H-band light distributions, but also to the stellar mass maps reconstructed from resolved stellar population modeling. We confirm that the increased bulge prominence among quiescent galaxies, as reported previously based on rest-optical observations, remains in place when considering the distributions of stellar mass. Moreover, we observe an increase of the typical Sersic index and bulge-to-total ratio (with median B/T reaching 40-50%) among star-forming galaxies above 10^11 Msun. Given that quenching for these most massive systems is likely to be imminent, our findings suggest that significant bulge growth precedes a departure from the star-forming main sequence. We demonstrate that the bulge mass (and ideally knowledge of the bulge and total mass) is a more reliable predictor of the star-forming versus quiescent state of a galaxy than the total stellar mass. The same trends are predicted by the state-of-the-art semi-analytic model by Somerville et al. In the latter, bulges and black holes grow hand in hand through merging and/or disk instabilities, and AGN-feedback shuts off star formation. Further observations will be required to pin down star formation quenching mechanisms, but our results imply they must be internal to the galaxies and closely associated with bulge growth.

We present the rest-frame UV wavelength dependence of the Petrosian-like half-light radius ($r_{50}$), and the concentration parameter for a sample of 198 star-forming galaxies at 0.5 < z < 1.5. We find a ~5% decrease in $r_{50}$ from 1500 \AA\ to 3000 \AA, with half-light radii at 3000 \AA\ ranging from 0.6 kpc to 6 kpc. We also find a decrease in concentration of ~0.07 (1.9 < $C_{3000}$ < 3.9). The lack of a strong relationship between $r_{50}$ and wavelength is consistent with a model in which clumpy star formation is distributed over length scales comparable to the galaxy's rest-frame optical light. While the wavelength dependence of $r_{50}$ is independent of size at all redshifts, concentration decreases more sharply in the far-UV (~1500 \AA) for large galaxies at z ~ 1. This decrease in concentration is caused by a flattening of the inner ~20% of the light profile in disk-like galaxies, indicating that the central regions have different UV colors than the rest of the galaxy. We interpret this as a bulge component with older stellar populations and/or more dust. The size-dependent decrease in concentration is less dramatic at z ~ 2, suggesting that bulges are less dusty, younger, and/or less massive than the rest of the galaxy at higher redshifts.

We present the results of a study investigating the rest-frame ultra-violet (UV) spectral slopes of redshift z~5 Lyman-break galaxies (LBGs). By combining deep Hubble Space Telescope imaging of the CANDELS and HUDF fields with ground-based imaging from the UKIDSS Ultra Deep Survey (UDS), we have produced a large sample of z~5 LBGs spanning an unprecedented factor of >100 in UV luminosity. Based on this sample we find a clear colour-magnitude relation (CMR) at z~5, such that the rest-frame UV slopes (beta) of brighter galaxies are notably redder than their fainter counterparts. We determine that the z~5 CMR is well described by a linear relationship of the form: d beta = (-0.12 +/- 0.02) d Muv, with no clear evidence for a change in CMR slope at faint magnitudes (i.e. Muv > -18.9). Using the results of detailed simulations we are able, for the first time, to infer the intrinsic (i.e. free from noise) variation of galaxy colours around the CMR at z~5. We find significant (12 sigma) evidence for intrinsic colour variation in the sample as a whole. Our results also demonstrate that the width of the intrinsic UV slope distribution of z~5 galaxies increases from Delta(beta)=0.1 at Muv=-18 to Delta(beta)=0.4 at Muv=-21. We suggest that the increasing width of the intrinsic galaxy colour distribution and the CMR itself are both plausibly explained by a luminosity independent lower limit of beta=-2.1, combined with an increase in the fraction of red galaxies in brighter UV-luminosity bins.

We derive physical properties of 10 submillimeter galaxies located in the CANDELS coverage of the GOODS-S field. The galaxies were first identified as submillimeter sources with the LABOCA bolometer and subsequently targeted for 870um continuum observation with ALMA. The high angular resolution of the ALMA imaging allows secure counterparts to be identified in the CANDELS multiband dataset. The CANDELS data provide deep photometric data from UV through near-infrared wavelengths. Using synthetic spectral energy distributions, we derive photometric redshifts, stellar masses, extinction, ages, and the star formation history. The redshift range is z=1.65-4.76, with two of the galaxies located at z>4. Two SMG counterparts have stellar masses 2-3 orders of magnitude lower than the rest. The remaining SMG counterparts have stellar masses around 1x10^11 Msun. The stellar population in the SMGs is typically older than the expected duration of the submillimeter phase, suggesting that the star formation history of submillimeter galaxies is more complex than a single burst. Non-parametric morphology indices suggest that the SMG counterparts are among the most asymmetric systems compared with galaxies of the same stellar mass and redshift. The HST images shows that 3 of the SMGs are associated with on-going mergers. The remaining counterparts are isolated. Estimating the dust and molecular gas mass from the submm fluxes, and comparing with our stellar masses shows that the molecular gas mass fraction of SMGs is ~28% and that the final stellar mass is likely to be (1-2)x10^11 Msun.

The processes that trigger Active Galactic Nuclei (AGN) remain poorly understood. While lower luminosity AGN may be triggered by minor disturbances to the host galaxy, stronger disturbances are likely required to trigger luminous AGN. Major wet mergers of galaxies are ideal environments for AGN triggering since they provide large gas supplies and galaxy scale torques. There is however little observational evidence for a strong connection between AGN and major mergers. We analyse the morphological properties of AGN host galaxies as a function of AGN and host galaxy luminosity and compare them to a carefully matched sample of control galaxies. AGN are X-ray selected in the redshift range 0.5 < z < 0.8 and have luminosities 41 < log(L_X [erg/s]) < 44.5. 'Fake AGN' are simulated in the control galaxies by adding point sources with the magnitude of the matched AGN. We find that AGN host and control galaxies have comparable assymetries, Sersic indices and ellipticities at restframe ~950nm. AGN host galaxies show neither higher average asymmetries nor higher fractions of very disturbed objects. There is no increase in the prevalence of merger signatures with AGN luminosity. At 95% confidence we find that major mergers are responsible for <6% of all AGN in our sample as well as <40% of the highest luminosity AGN log(L_X [erg/s]) ~ 43.5). Major mergers therefore either play only a very minor role in the triggering of AGN in the luminosity range studied or time delays are too long for merger features to remain visible.

Spectroscopic observations from the Large Binocular Telescope and the Very Large Telescope reveal kinematically narrow lines (~50 km/s) for a sample of 14 Extreme Emission Line Galaxies (EELGs) at redshifts 1.4 < z < 2.3. These measurements imply that the total dynamical masses of these systems are low (< 3x10^9 M_sun). Their large [O III] 5007 equivalent widths (500-1100 Angstroms) and faint blue continuum emission imply young ages of 10-100 Myr and stellar masses of 10^8-10^9 M_sun, confirming the presence of a violent starburst. The dynamical masses represent the first such determinations for low-mass galaxies at z > 1. The stellar mass formed in this vigorous starburst phase represents a large fraction of the total (dynamical) mass, without a significantly massive underlying population of older stars. The occurrence of such intense events in shallow potentials strongly suggests that supernova-driven winds must be of critical importance in the subsequent evolution of these systems.

Out of several dozen z > 7 candidate galaxies observed spectroscopically, only five have been confirmed via Lyman-alpha emission, at z=7.008, 7.045, 7.109, 7.213 and 7.215. The small fraction of confirmed galaxies may indicate that the neutral fraction in the intergalactic medium (IGM) rises quickly at z > 6.5, as Lyman-alpha is resonantly scattered by neutral gas. However, the small samples and limited depth of previous observations makes these conclusions tentative. Here we report the results of a deep near-infrared spectroscopic survey of 43 z > 6.5 galaxies. We detect only a single galaxy, confirming that some process is making Lyman-alpha difficult to detect. The detected emission line at 1.0343 um is likely to be Lyman-alpha emission, placing this galaxy at a redshift z = 7.51, an epoch 700 million years after the Big Bang. This galaxy's colors are consistent with significant metal content, implying that galaxies become enriched rapidly. We measure a surprisingly high star formation rate of 330 Msol/yr, more than a factor of 100 greater than seen in the Milky Way. Such a galaxy is unexpected in a survey of our size, suggesting that the early universe may harbor more intense sites of star-formation than expected.

We analyze the resolved stellar populations of 473 massive star-forming galaxies at 0.7 < z < 1.5, with multi-wavelength broad-band imaging from CANDELS and Halpha surface brightness profiles at the same kiloparsec resolution from 3D-HST. Together, this unique data set sheds light on how the assembled stellar mass is distributed within galaxies, and where new stars are being formed. We find the Halpha morphologies to resemble more closely those observed in the ACS I band than in the WFC3 H band, especially for the larger systems. We next derive a novel prescription for Halpha dust corrections, which accounts for extra extinction towards HII regions. The prescription leads to consistent SFR estimates and reproduces the observed relation between the Halpha/UV luminosity ratio and visual extinction, both on a pixel-by-pixel and on a galaxy-integrated level. We find the surface density of star formation to correlate with the surface density of assembled stellar mass for spatially resolved regions within galaxies, akin to the so-called 'main sequence of star formation' established on a galaxy-integrated level. Deviations from this relation towards lower equivalent widths are found in the inner regions of galaxies. Clumps and spiral features, on the other hand, are associated with enhanced Halpha equivalent widths, bluer colors, and higher specific star formation rates compared to the underlying disk. Their Halpha/UV luminosity ratio is lower than that of the underlying disk, suggesting the ACS clump selection preferentially picks up those regions of elevated star formation activity that are the least obscured by dust. Our analysis emphasizes that monochromatic studies of galaxy structure can be severely limited by mass-to-light ratio variations due to dust and spatially inhomogeneous star formation histories.

We present an overview of a 90-orbit Hubble Space Telescope treasury program to obtain near ultraviolet imaging of the Hubble Ultra Deep Field using the Wide Field Camera 3 UVIS detector with the F225W, F275W, and F336W filters. This survey is designed to: (i) Investigate the episode of peak star formation activity in galaxies at 1<z<2.5; (ii) Probe the evolution of massive galaxies by resolving sub-galactic units (clumps); (iii) Examine the escape fraction of ionizing radiation from galaxies at z~2-3; (iv) Greatly improve the reliability of photometric redshift estimates; and (v) Measure the star formation rate efficiency of neutral atomic-dominated hydrogen gas at z~1-3. In this overview paper, we describe the survey details and data reduction challenges, including both the necessity of specialized calibrations and the effects of charge transfer inefficiency. We provide a stark demonstration of the effects of charge transfer inefficiency on resultant data products, which when uncorrected, result in uncertain photometry, elongation of morphology in the readout direction, and loss of faint sources far from the readout. We agree with the STScI recommendation that future UVIS observations that require very sensitive measurements use the instrument's capability to add background light through a "post-flash". Preliminary results on number counts of UV-selected galaxies and morphology of galaxies at z~1 are presented. We find that the number density of UV dropouts at redshifts 1.7, 2.1, and 2.7 is largely consistent with the number predicted by published luminosity functions. We also confirm that the image mosaics have sufficient sensitivity and resolution to support the analysis of the evolution of star-forming clumps, reaching 28-29th magnitude depth at 5 sigma in a 0.2 arcsecond radius aperture depending on filter and observing epoch.

We use GOODS and CANDELS images to identify progenitors of massive (log M > 10 Msun) compact "early-type" galaxies (ETGs) at z~1.6. Since merging and accretion increase the size of the stellar component of galaxies, if the progenitors are among known star-forming galaxies, these must be compact themselves. We select candidate progenitors among compact Lyman-break galaxies at z~3 based on their mass, SFR and central stellar density and find that these account for a large fraction of, and possibly all, compact ETGs at z~1.6. We find that the average far-UV SED of the candidates is redder than that of the non-candidates, but the optical and mid-IR SED are the same, implying that the redder UV of the candidates is inconsistent with larger dust obscuration, and consistent with more evolved (aging) star-formation. This is in line with other evidence that compactness is a sensitive predictor of passivity among high-redshift massive galaxies. We also find that the light distribution of both the compact ETGs and their candidate progenitors does not show any extended "halos" surrounding the compact "core", both in individual images and in stacks. We argue that this is generally inconsistent with the morphology of merger remnants, even if gas-rich, as predicted by N-body simulations. This suggests that the compact ETGs formed via highly dissipative, mostly gaseous accretion of units whose stellar components are very small and undetected in the HST images, with their stellar mass assembling in-situ, and that they have not experienced any major merging until the epoch of observations at z~1.6.

Using spectroscopy from the Large Binocular Telescope and imaging from the Hubble Space Telescope we discovered the first strong galaxy lens at z(lens)>1. The lens has a secure photometric redshift of z=1.53+/-0.09 and the source is spectroscopically confirmed at z=3.417. The Einstein radius (0.35"; 3.0 kpc) encloses 7.6 x 10^10 Msol, with an upper limit on the dark matter fraction of 60%. The highly magnified (40x) source galaxy has a very small stellar mass (~10^8 Msol) and shows an extremely strong [OIII]_5007A emission line (EW_0 ~ 1000A) bolstering the evidence that intense starbursts among very low-mass galaxies are common at high redshift.

We present results from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) photometric redshift methods investigation. In this investigation, the results from eleven participants, each using a different combination of photometric redshift code, template spectral energy distributions (SEDs) and priors, are used to examine the properties of photometric redshifts applied to deep fields with broad-band multi-wavelength coverage. The photometry used includes U-band through mid-infrared filters and was derived using the TFIT method. Comparing the results, we find that there is no particular code or set of template SEDs that results in significantly better photometric redshifts compared to others. However, we find codes producing the lowest scatter and outlier fraction utilize a training sample to optimize photometric redshifts by adding zero-point offsets, template adjusting or adding extra smoothing errors. These results therefore stress the importance of the training procedure. We find a strong dependence of the photometric redshift accuracy on the signal-to-noise ratio of the photometry. On the other hand, we find a weak dependence of the photometric redshift scatter with redshift and galaxy color. We find that most photometric redshift codes quote redshift errors (e.g., 68% confidence intervals) that are too small compared to that expected from the spectroscopic control sample. We find that all codes show a statistically significant bias in the photometric redshifts. However, the bias is in all cases smaller than the scatter, the latter therefore dominates the errors. Finally, we find that combining results from multiple codes significantly decreases the photometric redshift scatter and outlier fraction. We discuss different ways of combining data to produce accurate photometric redshifts and error estimates.

We present a UV-to-mid infrared multi-wavelength catalog in the CANDELS/GOODS-S field, combining the newly obtained CANDELS HST/WFC3 F105W, F125W, and F160W data with existing public data. The catalog is based on source detection in the WFC3 F160W band. The F160W mosaic includes the data from CANDELS deep and wide observations as well as previous ERS and HUDF09 programs. The mosaic reaches a 5$\sigma$ limiting depth (within an aperture of radius 0.17 arcsec) of 27.4, 28.2, and 29.7 AB for CANDELS wide, deep, and HUDF regions, respectively. The catalog contains 34930 sources with the representative 50% completeness reaching 25.9, 26.6, and 28.1 AB in the F160W band for the three regions. In addition to WFC3 bands, the catalog also includes data from UV (U-band from both CTIO/MOSAIC and VLT/VIMOS), optical (HST/ACS F435W, F606W, F775W, F814W, and F850LP), and infrared (HST/WFC3 F098M, VLT/ISAAC Ks, VLT/HAWK-I Ks, and Spitzer/IRAC 3.6, 4.5, 5.8, 8.0 $\mu$m) observations. The catalog is validated via stellar colors, comparison with other published catalogs, zeropoint offsets determined from the best-fit templates of the spectral energy distribution of spectroscopically observed objects, and the accuracy of photometric redshifts. The catalog is able to detect unreddened star-forming (passive) galaxies with stellar mass of 10^{10}M_\odot at a 50% completeness level to z$\sim$3.4 (2.8), 4.6 (3.2), and 7.0 (4.2) in the three regions. As an example of application, the catalog is used to select both star-forming and passive galaxies at z$\sim$2--4 via the Balmer break. It is also used to study the color--magnitude diagram of galaxies at 0<z<4.

Projected axis ratio measurements of 880 early-type galaxies at redshifts 1<z<2.5 selected from CANDELS are used to reconstruct and model their intrinsic shapes. The sample is selected on the basis of multiple rest-frame colors to reflect low star-formation activity. We demonstrate that these galaxies as an ensemble are dust-poor and transparent and therefore likely have smooth light profiles, similar to visually classified early-type galaxies. Similar to their present-day counterparts, the z>1 early-type galaxies show a variety of intrinsic shapes; even at a fixed mass, the projected axis ratio distributions cannot be explained by the random projection of a set of galaxies with very similar intrinsic shapes. However, a two-population model for the intrinsic shapes, consisting of a triaxial, fairly round population, combined with a flat (c/a~0.3) oblate population, adequately describes the projected axis ratio distributions of both present-day and z>1 early-type galaxies. We find that the proportion of oblate versus triaxial galaxies depends both on the galaxies' stellar mass, and - at a given mass - on redshift. For present-day and z<1 early-type galaxies the oblate fraction strongly depends on galaxy mass. At z>1 this trend is much weaker over the mass range explored here (10^10<M*/M_sun<10^11), because the oblate fraction among massive (M*~10^11 M_sun) was much higher in the past: 0.59+-0.10 at z>1, compared to 0.20+-0.02 at z~0.1. In contrast, the oblate fraction among low-mass early-type galaxies (log(M*/M_sun)<10.5) increased toward the present, from 0.38+-0.11 at z>1 to 0.72+-0.06 at z=0. [Abridged]

We present evidence for a strong relationship between galaxy size and environment for the quiescent population in the redshift range 1 < z < 2. Environments were measured using projected galaxy overdensities on a scale of 400 kpc, as determined from ~ 96,000 K-band selected galaxies from the UKIDSS Ultra Deep Survey (UDS). Sizes were determined from ground-based K-band imaging, calibrated using space-based CANDELS HST observations in the centre of the UDS field, with photometric redshifts and stellar masses derived from 11-band photometric fitting. From the resulting size-mass relation, we confirm that quiescent galaxies at a given stellar mass were typically ~ 50 % smaller at z ~ 1.4 compared to the present day. At a given epoch, however, we find that passive galaxies in denser environments are on average significantly larger at a given stellar mass. The most massive quiescent galaxies (M_stellar > 2 x 10^11 M_sun) at z > 1 are typically 50 % larger in the highest density environments compared to those in the lowest density environments. Using Monte Carlo simulations, we reject the null hypothesis that the size-mass relation is independent of environment at a significance > 4.8 sigma for the redshift range 1 < z < 2. In contrast, the evidence for a relationship between size and environment is much weaker for star-forming galaxies.

We discuss the state of the assembly of the Hubble Sequence in the mix of bright galaxies at redshift 1.4< z \le 2.5 with a large sample of 1,671 galaxies down to H_{AB}~26, selected from the HST/ACS and WFC3 images of the GOODS--South field obtained as part of the GOODS and CANDELS observations. We investigate the relationship between the star formation properties and morphology using various parametric diagnostics, such as the Sersic light profile, Gini (G), M_{20}, Concentration (C), Asymmetry (A) and multiplicity parameters. Our sample clearly separates into massive, red and passive galaxies versus less massive, blue and star forming ones, and this dichotomy correlates very well with the galaxies' morphological properties. Star--forming galaxies show a broad variety of morphological features, including clumpy structures and bulges mixed with faint low surface brightness features, generally characterized by disky-type light profiles. Passively evolving galaxies, on the other hand, very often have compact light distribution and morphology typical of today's spheroidal systems. We also find that artificially redshifted local galaxies have a similar distribution with z~2galaxies in a G-M_{20} plane. Visual inspection between the rest-frame optical and UV images show that there is a generally weak morphological k-correction for galaxies at z~2, but the comparison with non-parametric measures show that galaxies in the rest-frame UV are somewhat clumpier than rest-frame optical. Similar general trends are observed in the local universe among massive galaxies, suggesting that the backbone of the Hubble sequence was already in place at z~2.

We present a full analysis of the Probing Evolution And Reionization Spectroscopically (PEARS) slitess grism spectroscopic data obtained with the Advanced Camera for Surveys on HST. PEARS covers fields within both the Great Observatories Origins Deep Survey (GOODS) North and South fields, making it ideal as a random survey of galaxies, as well as the availability of a wide variety of ancillary observations to support the spectroscopic results. Using the PEARS data we are able to identify star forming galaxies within the redshift volume 0< z<1.5. Star forming regions in the PEARS survey are pinpointed independently of the host galaxy. This method allows us to detect the presence of multiple emission line regions (ELRs) within a single galaxy. 1162 Ha, [OIII] and/or [OII] emission lines have been identified in the PEARS sample of ~906 galaxies down to a limiting flux of ~1e-18 erg/s/cm^2. The ELRs have also been compared to the properties of the host galaxy, including morphology, luminosity, and mass. From this analysis we find three key results: 1) The computed line luminosities show evidence of a flattening in the luminosity function with increasing redshift; 2) The star forming systems show evidence of disturbed morphologies, with star formation occurring predominantly within one effective (half-light) radius. However, the morphologies show no correlation with host stellar mass; and 3) The number density of star forming galaxies with M_* > 1e9} M_sun decreases by an order of magnitude at z<0.5 relative to the number at 0.5<z<0.9 in support of the argument for galaxy downsizing.

We present global structural parameter measurements of 109,533 unique, H_F160W-selected objects from the CANDELS multi-cycle treasury program. Sersic model fits for these objects are produced with GALFIT in all available near-infrared filters (H_F160W, J_F125W and, for a subset, Y_F105W). The parameters of the best-fitting Sersic models (total magnitude, half-light radius, Sersic index, axis ratio, and position angle) are made public, along with newly constructed point spread functions for each field and filter. Random uncertainties in the measured parameters are estimated for each individual object based on a comparison between multiple, independent measurements of the same set of objects. To quantify systematic uncertainties we create a mosaic with simulated galaxy images with a realistic distribution of input parameters and then process and analyze the mosaic in an identical manner as the real data. We find that accurate and precise measurements -- to 10% or better -- of all structural parameters can typically be obtained for galaxies with H_F160W < 23, with comparable fidelity for basic size and shape measurements for galaxies to H_F160W ~ 24.5.

In this paper we present a detailed study of the structures and morphologies of a sample of 1188 massive galaxies with Mstar>10^10Msun between redshifts z=1-3 within the Ultra Deep Survey (UDS) region of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) field. Using this sample we determine how galaxy structure and morphology evolve with time. We visually classify our sample into disks, ellipticals and peculiar systems and correct for redshift effects on our classifications through simulations. We find evolution in the fractions of galaxies at a given visual classification as a function of redshift. The peculiar population is dominant at z>2 with a substantial spheroid population, and a negligible disk population. We compute the transition redshift, ztrans, where the combined fraction of spheroids and disks is equal to that of peculiar galaxies, as ztrans=1.86+/-0.62 for galaxies in our stellar mass range. We find that this changes as a function of stellar mass, with Hubble-type systems becoming dominant at higher redshifts for higher mass galaxies (ztrans=2.22+/-0.82), than for the lower mass galaxies (ztrans=1.73+/-0.57). Higher mass galaxies become morphologically settled before their lower mass counterparts, a form of morphological downsizing. We furthermore compare our visual classifications with Sersic index, the concentration, asymmetry and clumpiness (CAS) parameters, star formation rate and rest frame U-B colour. We find links between the colour of a galaxy, its star formation rate and how extended or peculiar it appears. Finally, we discuss the negligible z>2 disk fraction based on visual morphologies and speculate that this is an effect of forming disks appearing peculiar through processes such as violent disk instabilities or mergers. We conclude that to properly define high redshift morphology a new and more exact classification scheme is needed.

We present observations of a luminous galaxy at redshift z=6.573 --- the end of the reioinization epoch --- which has been spectroscopically confirmed twice. The first spectroscopic confirmation comes from slitless HST ACS grism spectra from the PEARS survey (Probing Evolution And Reionization Spectroscopically), which show a dramatic continuum break in the spectrum at restframe 1216 A wavelength. The second confirmation is done with Keck + DEIMOS. The continuum is not clearly detected with ground-based spectra, but high wavelength resolution enables the Lyman alpha emission line profile to be determined. We compare the line profile to composite line profiles at redshift z=4.5. The Lyman alpha line profile shows no signature of a damping wing attenuation, confirming that the intergalactic gas is ionized at redshift z=6.57. Spectra of Lyman breaks at yet higher redshifts will be possible using comparably deep observations with IR-sensitive grisms, even at redshifts where Lyman alpha is too attenuated by the neutral IGM to be detectable using traditional spectroscopy from the ground.

We have exploited the HST CANDELS WFC3/IR imaging to study the properties of (sub-)mm galaxies in GOODS-South. After using the deep radio and Spitzer imaging to identify galaxy counterparts for the (sub-)mm sources, we have used the new CANDELS data in two ways. First, we have derived improved photometric redshifts and stellar masses, confirming that the (sub-)mm galaxies are massive (<M*>=2.2x10^11 M_solar) galaxies at z=1-3. Second, we have exploited the depth and resolution of the WFC3/IR imaging to determine the sizes and morphologies of the galaxies at rest-frame optical wavelengths, fitting two-dimensional axi-symmetric Sersic models. Crucially, the WFC3/IR H-band imaging enables modelling of the mass-dominant galaxy, rather than the blue high-surface brightness features which often dominate optical (rest-frame UV) images of (sub-)mm galaxies, and can confuse visual morphological classification. As a result of this analysis we find that >95% of the rest-frame optical light in almost all of the (sub-)mm galaxies is well-described by either a single exponential disk, or a multiple-component system in which the dominant constituent is disk-like. We demonstrate that this conclusion is consistent with the results of high-quality ground-based K-band imaging, and explain why. The massive disk galaxies which host luminous (sub-)mm emission are reasonably extended (r_e=4 kpc), consistent with the sizes of other massive star-forming disks at z~2. In many cases we find evidence of blue clumps within the sources, with the mass-dominant disk becoming more significant at longer wavelengths. Finally, only a minority of the sources show evidence for a major galaxy-galaxy interaction. Taken together, these results support the view that most (sub-)mm galaxies at z~2 are simply the most extreme examples of normal star-forming galaxies at that era.

We present the multiwavelength - ultraviolet to mid-infrared - catalog of the UKIDSS Ultra-Deep Survey (UDS) field observed as part of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). Based on publicly available data, the catalog includes: the CANDELS data from the Hubble Space Telescope (near-infrared WFC3 F125W and F160W data and visible ACS F606W and F814W data), u-band data from CFHT/Megacam, B, V, Rc, i' and z' band data from Subaru/Suprime-Cam, Y and Ks band data from VLT/HAWK-I, J, H and K bands data from UKIDSS (Data Release 8), and Spitzer/IRAC data (3.6, 4.5 from SEDS, 5.8 and 8.0um from SpUDS). The present catalog is F160W-selected and contains 35932 sources over an area of 201.7 square arcmin and includes radio and X-ray detected sources and spectroscopic redshifts available for 210 sources.

We have made a serendipitous discovery of a massive cD galaxy at z=1.096 in a candidate rich cluster in the HUDF area of GOODS-South. This brightest cluster galaxy is the most distant cD galaxy confirmed to date. Ultra-deep HST/WFC3 images reveal an extended envelope starting from ~10 kpc and reaching ~70 kpc in radius along the semi-major axis. The spectral energy distributions indicate that both its inner component and outer envelope are composed of an old, passively-evolving stellar population. The cD galaxy lies on the same mass-size relation as the bulk of quiescent galaxies at similar redshifts. The cD galaxy has a higher stellar mass surface density but a similar velocity dispersion to those of more-massive, nearby cDs. If the cD galaxy is one of the progenitors of today's more massive cDs, its size and stellar mass have had to increase on average by factors of $3.4\pm1.1$ and $3.3\pm1.3$ over the past ~8 Gyrs, respectively. Such increases in size and stellar mass without being accompanied by significant increases in velocity dispersion are consistent with evolutionary scenarios driven by both major and minor dry mergers. If such cD envelopes originate from dry mergers, our discovery of even one example proves that some BCGs entered the dry merger phase at epochs earlier than z=1. Our data match theoretical models which predict that the continuance of dry mergers at z<1 can result in structures similar to those of massive cD galaxies seen today. Moreover, our discovery is a surprise given that the extreme depth of the HUDF is essential to reveal such an extended cD envelope at z>1 and, yet, the HUDF covers only a minuscule region of sky. Adding that cDs are rare, Our serendipitous discovery hints that such cDs may be more common than expected. [Abridged]

We present the discovery of a Type Ia supernova (SN) at redshift $z = 1.914$ from the CANDELS multi-cycle treasury program on the \textit{Hubble Space Telescope (HST)}. This SN was discovered in the infrared using the Wide-Field Camera 3, and it is the highest-redshift Type Ia SN yet observed. We classify this object as a SN\,Ia by comparing its light curve and spectrum with those of a large sample of Type Ia and core-collapse supernovae (SNe). Its apparent magnitude is consistent with that expected from the $\Lambda$CDM concordance cosmology. We discuss the use of spectral evidence for classification of $z > 1.5$ SNe\,Ia using {\it HST} grism simulations, finding that spectral data alone can frequently rule out SNe\,II, but distinguishing between SNe\,Ia and SNe\,Ib/c can require prohibitively long exposures. In such cases, a quantitative analysis of the light curve may be necessary for classification. Our photometric and spectroscopic classification methods can aid the determination of SN rates and cosmological parameters from the full high-redshift CANDELS SN sample.

Star-forming galaxies at redshifts z>6 are likely responsible for the reionization of the universe, and it is important to study the nature of these galaxies. We present three candidates for z~7 Lyman-break galaxies (LBGs) from a 155 arcmin^2 area in the CANDELS/COSMOS field imaged by the deep FourStar Galaxy Evolution (zFourGE) survey. The FourStar medium-band filters provide the equivalent of R~10 spectroscopy, which cleanly distinguishes between z~7 LBGs and brown dwarf stars. The distinction between stars and galaxies based on an object's angular size can become unreliable even when using HST imaging; there exists at least one very compact z~7 candidate (FWHM~0.5-1 kpc) that is indistinguishable from a point source. The medium-band filters provide narrower redshift distributions compared with broad-band-derived redshifts. The UV luminosity function derived using the three z~7 candidates is consistent with previous studies, suggesting an evolution at the bright end (MUV -21.6 mag) from z~7 to z~5. Fitting the galaxies' spectral energy distributions, we predict Lyman-alpha equivalent widths for the two brightest LBGs, and find that the presence of a Lyman-alpha line affects the medium-band flux thereby changing the constraints on stellar masses and UV spectral slopes. This illustrates the limitations of deriving LBG properties using only broad-band photometry. The derived specific star-formation rates for the bright LBGs are ~13 per Gyr, slightly higher than the lower-luminosity LBGs, implying that the star-formation rate increases with stellar mass for these galaxies.

We study the evolution of the number density, as a function of the size, of passive early-type galaxies with a wide range of stellar masses 10^10<M*/Msun<10^11.5) from z~3 to z~1, exploiting the unique dataset available in the GOODS-South field, including the recently obtained WFC3 images as a part of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). In particular, we select a sample of 107 massive (M*>10^10 M_sun), passive (SSFR<10^-2 Gyr^-1) and morphologically spheroidal galaxies at 1.2<z<3, taking advantage of the panchromatic dataset available for GOODS, including VLT, CFHT, Spitzer, Chandra and HST ACS+WFC3 data. We find that at 1<z<3 the passively evolving early-type galaxies are the reddest and most massive objects in the Universe, and we prove that a correlation between mass, morphology, color and star-formation activity is already in place at that epoch. We measure a significant evolution in the mass-size relation of passive early-type galaxies (ETGs) from z~3 to z~1, with galaxies growing on average by a factor of 2 in size in a 3 Gyr timescale only. We witness also an increase in the number density of passive ETGs of 50 times over the same time interval. We find that the first ETGs to form at z>2 are all compact or ultra-compact, while normal sized ETGs (meaning ETGs with sizes comparable to those of local counterparts of the same mass) are the most common ETGs only at z<1. The increase of the average size of ETGs at 0<z<1 is primarily driven by the appearance of new large ETGs rather than by the size increase of individual galaxies.

We use multi-band spatially resolved photometry from the Cosmic Assembly Near-IR Deep Legacy Survey (CANDELS) in the 4 MSec Chandra Deep Field-South (CDF-S) to explore the nuclear and extended colors, color gradients and stellar populations of X-ray selected AGN host galaxies out to z=3. Based on a study of their central light, we develop X-ray based criteria to exclude objects with strong AGN contamination. We use stellar masses from the FIREWORKS database to understand and account for stellar mass selection effects, and carefully study, for the first time, the resolved host galaxy properties of AGNs at z~2 in their rest-frame optical light without substantial nuclear contamination. AGN hosts span a sizable range of stellar masses, colors and color gradients at these redshifts. Their colors, color gradients and stellar population properties are very similar to inactive galaxies of the same stellar mass. At z~1, we find a slightly narrower range in host colors compared to inactive galaxies, as well as hints of more recent star-formation. These differences are weaker or non-existent among AGN hosts at z~2. We discuss the importance of AGN driven feedback in the quenching of galaxies at z>~1 and speculate on possible evolution in the relationship between black hole accretion and the host galaxy towards high redshifts.

We have used high-resolution, HST WFC3/IR, near-infrared imaging to conduct a detailed bulge-disk decomposition of the morphologies of ~200 of the most massive (M_star > 10^11 M_solar) galaxies at 1<z<3 in the CANDELS-UDS field. We find that, while such massive galaxies at low redshift are generally bulge-dominated, at redshifts 1<z<2 they are predominantly mixed bulge+disk systems, and by z>2 they are mostly disk-dominated. Interestingly, we find that while most of the quiescent galaxies are bulge-dominated, a significant fraction (25-40%) of the most quiescent galaxies, have disk-dominated morphologies. Thus, our results suggest that the physical mechanisms which quench star-formation activity are not simply connected to those responsible for the morphological transformation of massive galaxies.

One key goal of the Hubble Space Telescope Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey is to track galaxy evolution back to z ~ 8. Its two-tiered "wide and deep" strategy bridges significant gaps in existing near-infrared surveys. Here we report on z ~ 8 galaxy candidates selected as F105W-band dropouts in one of its deep fields, which covers 50.1 square arcmin to 4 ks depth in each of three near-infrared bands in the Great Observatories Origins Deep Survey southern field. Two of our candidates have J<26.2 mag, and are > 1 mag brighter than any previously known F105W-dropouts. We derive constraints on the bright-end of the rest-frame ultraviolet luminosity function of galaxies at z ~ 8, and show that the number density of such very bright objects is higher than expected from the previous Schechter luminosity function estimates at this redshift. Another two candidates are securely detected in Spitzer Infrared Array Camera images, which are the first such individual detections at z ~ 8. Their derived stellar masses are on the order of a few x 10^9 M_sun, from which we obtain the first measurement of the high-mass end of the galaxy stellar mass function at z ~ 8. The high number density of very luminous and very massive galaxies at z ~ 8, if real, could imply a large stellar-to-halo mass ratio and an efficient conversion of baryons to stars at such an early time.

The exploration of the relation between galaxy sizes and other physical parameters has provided important clues for understanding galaxy formation. We use the CANDELS Deep+Wide surveys in the GOODS-South, UDS and EGS fields, complemented by data from the HUDF09 program, to address the relation between size and luminosity at z\sim7. We select 153 z-band drop-out galaxies in six different fields characterized by a wide combination of depth and areal coverage, ideally suited to sample without biases the observed size-magnitude plane. Detailed simulations allow us to derive the completeness as a function of size and magnitude and to quantify measurements errors/biases. We find that the half light radius distribution function of z\sim7 galaxies fainter than J=26.6 is peaked at <0.1 arcsec (or equivalently 0.5 kpc proper), while at brighter magnitudes high-z galaxies are typically larger than ~0.15 arcsec. We also find a well defined size-luminosity relation, Rh\simL^1/2. We compute the Luminosity Function in the HUDF and P12HUDF fields, finding large spatial variation on the number density of faint galaxies. Adopting the size distribution and the size-luminosity relation found for faint galaxies at z=7, we derive a mean slope of -1.7\pm0.1 for the luminosity function of LBGs at this redshift. Using this LF, we find that the amount of ionizing photons cannot keep the Universe re-ionized if the IGM is clumpy (C_HII>3) and the Lyman continuum escape fraction of high-z LBGs is relatively low (f_esc<0.3). With future CANDELS data, we can put severe limits to the role of galaxies in the reionization of the Universe.

Using deep 100-160 micron observations in GOODS-S from the GOODS-H survey, combined with HST/WFC3 NIR imaging from CANDELS, we present the first morphological analysis of a complete, FIR selected sample of 52 ULIRGs at z~2. We also make use of a comparison sample of galaxies without Herschel detections but with the same z and magnitude distribution. Our visual classifications of these two samples indicate that the fraction of objects with disk and spheroid morphologies is roughly the same but that there are significantly more mergers, interactions, and irregular galaxies among the ULIRGs. The combination of disk and irregular/interacting morphologies suggests that early stage interactions and minor mergers could play an important role in ULIRGs at z~2. We compare these fractions with those of a z~1 sample across a wide luminosity range and find that the fraction of disks decreases systematically with L_IR while the fraction of mergers and interactions increases, as has been observed locally. At comparable luminosities, the fraction of ULIRGs with various morphological classifications is similar at z~2 and z~1. We investigate the position of the ULIRGs, along with 70 LIRGs, on the specific star formation rate versus redshift plane, and find 52 systems to be starbursts (lie more than a factor of 3 above the main sequence relation). The morphologies of starbursts are dominated by interacting and merging systems (50%). If irregular disks are included as potential minor mergers, then we find that up to 73% of starbursts are involved in a merger or interaction at some level. Although the final coalescence of a major merger may not be required for the high luminosities of ULIRGs at z~2 as is the case locally, the large fraction of interactions at all stages and potential minor mergers suggest that the high star formation rates of ULIRGs are still largely externally triggered at z~2.

We study the evolution of galaxy rest-frame ultraviolet (UV) colors in the epoch 4 < z < 8. We use new wide-field near-infrared data in GOODS-S from the CANDELS, HUDF09 and ERS programs to select galaxies via photometric redshift measurements. Our sample consists of 2812 candidate galaxies at z > 3.5, including 113 at z = 7 to 8. We fit the observed spectral energy distribution to a suite of synthetic stellar population models, and measure the value of the UV spectral slope (beta) from the best-fit model spectrum. The median value of beta evolves significantly from -1.82 (+0.00,-0.04) at z = 4, to -2.37 (+0.26,-0.06) at z = 7. Additionally, we find that faint galaxies at z = 7 have beta = -2.68 (+0.39,-0.24) (~ -2.4 after correcting for observational bias); this is redder than previous claims in the literature, and does not require "exotic" stellar populations to explain their colors. This evolution can be explained by an increase in dust extinction, with the timescale consistent with low-mass AGB stars forming the bulk of the dust. We find no significant (< 2-sigma) correlation between beta and M_UV when measuring M_UV at a consistent rest-frame wavelength of 1500 A. This is particularly true at bright magnitudes, though our results do show evidence for a weak correlation at faint magnitudes when galaxies in the HUDF are considered separately, hinting that dynamic range in sample luminosities may play a role. We do find a strong correlation between beta and the stellar mass at all redshifts, in that more massive galaxies exhibit redder colors. The most massive galaxies in our sample have red colors at each redshift, implying that dust can build up quickly in massive galaxies, and that feedback is likely removing dust from low-mass galaxies at z > 7. Thus the stellar-mass - metallicity relation, previously observed up to z ~ 3, may extend out to z = 7 - 8.

We report the results of a study exploring the stellar populations of 13 luminous (L>L*), spectroscopically confirmed, galaxies in the redshift interval 5.5<z<6.5, all with WFC3/IR and IRAC imaging from the HST/CANDELS and Spitzer/SEDS surveys. Based on fitting the observed photometry with SED templates covering a wide range of different star-formation histories, and a self consistent treatment of Lyman-alpha emission, we find that the derived stellar masses lie within the range 10^9 Msun < M*< 10^10 Msun and are robust to within a factor of two. In contrast, we confirm previous reports that the ages of the stellar populations are poorly constrained. Although the best-fitting models for three objects have ages >= 300 Myr, the degeneracies introduced by dust extinction mean that only two of these objects actually require a >300 Myr old stellar population to reproduce the observed photometry. Moreover, when considering only smoothly-varying star-formation histories, we observe a clear tension between the data and models such that a galaxy SED template with an old age is often chosen in order to try and fit objects with blue UV-slopes but red UV-to-optical colours. To break this tension we explore SED fitting with two-component models (burst plus on-going star-formation) and allow for nebular emission. On average, the inclusion of nebular emission leads to lower stellar-mass estimates (median offset 0.18 dex), moderately higher specific star-formation rates, and allows for a wider range of plausible stellar ages. However, based on our SED modelling, we find no strong evidence for extremely young ages in our sample (<50 Myr). Finally, considering all of the different star-formation histories explored, we find that the median best-fitting ages are of the order 200-300 Myr and that the objects with the tightest constraints indicate ages in the range 50-200 Myr (Abridged).

Using a sample of BzK-selected galaxies at z~2 identified from the CFHT/WIRCAM near-infrared survey of GOODS-North, we discuss the relation between star formation rate (SFR), specific star formation rate (SSFR), and stellar mass (M_{*}), and the clustering of galaxies as a function of these parameters. For star-forming galaxies (sBzKs), the UV-based SFR, corrected for extinction, scales with the stellar mass as SFR ~ M_{*}^{alpha} with alpha = 0.74+/-0.20 down to M_{*} ~ 10^{9} M_{solar}, indicating a weak dependence on the stellar mass of the SSFR. We also measure the angular correlation function and hence infer the correlation length for sBzK galaxies as a function of M_{*}, SFR, and SSFR, as well as K-band apparent magnitude. We show that passive galaxies (pBzKs) are more strongly clustered than sBzK galaxies at a given stellar mass, mirroring the color-density relation seen at lower redshifts. We also find that the correlation length of sBzK galaxies ranges from 4 to 20 h^{-1}Mpc, being a strong function of M_{K}, M_{*}, and SFR. On the other hand, the clustering dependence on SSFR changes abruptly at 2x10^{-9} yr^{-1}, which is the typical value for "main sequence" star-forming galaxies at z~2. We show that the correlation length reaches a minimum at this characteristic value, and is larger for galaxies with both smaller and larger SSFRs; a dichotomy that is only marginally implied from the predictions of the semi-analytical models. Our results suggest that there are two types of environmental effects at work at z~2. Stronger clustering for relatively quiescent galaxies implies that the environment has started to play a role in quenching star formation. At the same time, stronger clustering for galaxies with elevated SSFRs ("starbursts") might be attributed to an increased efficiency for galaxy interactions and mergers in dense environments.

We combine high-resolution HST/WFC3 images with multi-wavelength photometry to track the evolution of structure and activity of massive (log(M*) > 10) galaxies at redshifts z = 1.4 - 3 in two fields of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). We detect compact, star-forming galaxies (cSFGs) whose number densities, masses, sizes, and star formation rates qualify them as likely progenitors of compact, quiescent, massive galaxies (cQGs) at z = 1.5 - 3. At z > 2 most cSFGs have specific star-formation rates (sSFR = 10^-9 yr^-1) half that of typical, massive SFGs at the same epoch, and host X-ray luminous AGN 30 times (~30%) more frequently. These properties suggest that cSFGs are formed by gas-rich processes (mergers or disk-instabilities) that induce a compact starburst and feed an AGN, which, in turn, quench the star formation on dynamical timescales (few 10^8 yr). The cSFGs are continuously being formed at z = 2 - 3 and fade to cQGs by z = 1.5. After this epoch, cSFGs are rare, thereby truncating the formation of new cQGs. Meanwhile, down to z = 1, existing cQGs continue to enlarge to match local QGs in size, while less-gas-rich mergers and other secular mechanisms shepherd (larger) SFGs as later arrivals to the red sequence. In summary, we propose two evolutionary scenarios of QG formation: an early (z > 2), fast-formation path of rapidly-quenched cSFGs that evolve into cQGs that later enlarge within the quiescent phase, and a slow, late-arrival (z < 2) path for SFGs to form QGs without passing through a compact state.

We have used deep, HST, near-IR imaging to study the morphological properties of the most massive galaxies at high z, modelling the WFC3/IR H-band images of the ~200 galaxies in the CANDELS-UDS field with 1 < z_phot < 3, and stellar masses M_star > 10^11 M_sun. We have used both single-Sersic and bulge+disk models, have investigated the errors/biases introduced by uncertainties in the background and the PSF, and have obtained formally-acceptable model fits to >90% of the galaxies. Our results indicate that these massive galaxies at 1 < z < 3 lie both on and below the local size-mass relation, with a median R_e~2.6 kpc, a factor of ~2.3 smaller than comparably-massive local galaxies. Moreover, we find that bulge-dominated objects in particular show evidence for a growing bimodality in the size-mass relation with increasing z, and by z > 2 the compact bulges display effective radii a factor ~4 smaller than local ellipticals of comparable mass. These trends appear to extend to the bulge components of disk-dominated galaxies, and vice versa. We also find that, while such massive galaxies at low z are bulge-dominated, at 1 < z < 2 they are predominantly mixed bulge+disk systems, and by z > 2 they are mostly disk-dominated. The majority of the disk-dominated galaxies are actively forming stars, but this is also true for many of the bulge-dominated systems. Interestingly, however, while most of the quiescent galaxies are bulge-dominated, we find that a significant fraction (25-40%) of the most quiescent galaxies have disk-dominated morphologies. Thus, while our results show that the massive galaxy population is undergoing dramatic changes at this crucial epoch, they also suggest that the physical mechanisms which quench star-formation activity are not simply connected to those responsible for the morphological transformation of massive galaxies into present-day giant ellipticals.

We present measurements of the specific ultraviolet luminosity density from a sample of 483 galaxies at 6<z<8. These galaxies were selected from new deep near-infrared HST imaging from the CANDELS, HUDF09 and ERS programs. In contrast to the majority of previous analyses, which assume that the distribution of galaxy ultraviolet (UV) luminosities follows a Schechter distribution, and that the distribution continues to luminosities far below our observable limit, we investigate the contribution to reionization from galaxies which we can observe, free from these assumptions. We find that the observable population of galaxies can sustain a fully reionized IGM at z=6, if the average ionizing photon escape fraction (f_esc) is ~30%. A number of previous studies have measured UV luminosity densities at these redshifts that vary by 5X, with many concluding that galaxies could not complete reionization by z=6 unless a large population of galaxies fainter than the detection limit were invoked, or extremely high values of f_esc were present. The observed UV luminosity density from our observed galaxy samples at z=7-8 is not sufficient to maintain a fully reionized IGM unless f_esc>50%. Combining our observations with constraints on the emission rate of ionizing photons from Ly-alpha forest observations at z=6, we can constrain f_esc<34% (2-sigma) if the observed galaxies are the only contributors to reionization, or <13% (2-sigma) if the luminosity function extends to M_UV = -13. These escape fractions are sufficient to complete reionization by z=6. These constraints imply that the volume ionized fraction of the IGM becomes less than unity at z>7, consistent with a number of complementary reionization probes. If faint galaxies dominate reionization, future JWST observations will probe deep enough to see them, providing an indirect constraint on the ionizing photon escape fraction [abridged].

We use HST/WFC3 imaging from the CANDELS Multicycle Treasury Survey, in conjunction with the Sloan Digital Sky Survey, to explore the evolution of galactic structure for galaxies with stellar masses >3e10M_sun from z=2.2 to the present epoch, a time span of 10Gyr. We explore the relationship between rest-frame optical color, stellar mass, star formation activity and galaxy structure. We confirm the dramatic increase from z=2.2 to the present day in the number density of non-star-forming galaxies above 3e10M_sun reported by others. We further find that the vast majority of these quiescent systems have concentrated light profiles, as parametrized by the Sersic index, and the population of concentrated galaxies grows similarly rapidly. We examine the joint distribution of star formation activity, Sersic index, stellar mass, inferred velocity dispersion, and stellar surface density. Quiescence correlates poorly with stellar mass at all z<2.2. Quiescence correlates well with Sersic index at all redshifts. Quiescence correlates well with `velocity dispersion' and stellar surface density at z>1.3, and somewhat less well at lower redshifts. Yet, there is significant scatter between quiescence and galaxy structure: while the vast majority of quiescent galaxies have prominent bulges, many of them have significant disks, and a number of bulge-dominated galaxies have significant star formation. Noting the rarity of quiescent galaxies without prominent bulges, we argue that a prominent bulge (and perhaps, by association, a supermassive black hole) is an important condition for quenching star formation on galactic scales over the last 10Gyr, in qualitative agreement with the AGN feedback paradigm.

We perform a detailed analysis of the resolved colors and stellar populations of a complete sample of 323 star-forming galaxies at 0.5 < z < 1.5, and 326 star-forming galaxies at 1.5 < z < 2.5 in the ERS and CANDELS-Deep region of GOODS-South. Galaxies were selected to be more massive than 10^10 Msun and have specific star formation rates above 1/t_H. We model the 7-band optical ACS + near-IR WFC3 spectral energy distributions of individual bins of pixels, accounting simultaneously for the galaxy-integrated photometric constraints available over a longer wavelength range. We analyze variations in rest-frame color, stellar surface mass density, age, and extinction as a function of galactocentric radius and local surface brightness/density, and measure structural parameters on luminosity and stellar mass maps. We find evidence for redder colors, older stellar ages, and increased dust extinction in the nuclei of galaxies. Big star-forming clumps seen in star formation tracers are less prominent or even invisible on the inferred stellar mass distributions. Off-center clumps contribute up to ~20% to the integrated SFR, but only 7% or less to the integrated mass of all massive star-forming galaxies at z ~ 1 and z ~ 2, with the fractional contributions being a decreasing function of wavelength used to select the clumps. The stellar mass profiles tend to have smaller sizes and M20 coefficients, and higher concentration and Gini coefficients than the light distribution. Our results are consistent with an inside-out disk growth scenario with brief (100 - 200 Myr) episodic local enhancements in star formation superposed on the underlying disk. Alternatively, the young ages of off-center clumps may signal inward clump migration, provided this happens efficiently on the order of an orbital timescale.

We present a study on Spectral Energy Distributions, Morphologies, and star formation for an IRAC-selected extremely red object sample in the GOODS Chandra Deep Field-South. This work was enabled by new HST/WFC3 near-IR imaging from the CANDELS survey as well as the deepest available X-ray data from Chandra 4 Ms observations. This sample consists of 133 objects with the 3.6um limiting magnitude of [3.6] = 21.5, and is approximately complete for galaxies with M >10^{11}M_sun at 1.5 < z < 2.5. We classify this sample into two types, quiescent and star-forming galaxies, in the observed infrared color-color ([3.6]-[24] vs K-[3.6]) diagram. The further morphological study of this sample show a consistent result with the observed color classification. The classified quiescent galaxies are bulge dominated and star-forming galaxies in the sample have disk or irregular morphologies. Our observed infrared color classification is also consistent with the rest-frame color (U-V vs V-J) classification. We also found that quiescent and star-forming galaxies are well separated in the nonparametric morphology parameter (Gini vs M_{20}) diagram measuring their concentration and clumpiness: quiescent galaxies have Gini coefficient higher than 0.58 and star forming galaxies have Gini coefficient lower that 0.58. We argue that the star formation quenching process must lead to or be accompanied by the increasing galaxy concentration. One prominent morphological feature of this sample is that disks are commonly seen in this massive galaxy sample at 1.5 < z < 2.5: 30% of quiescent galaxies and 70% of star forming galaxies with M >10^{11}M_sun have disks in their rest-frame optical morphologies. The prevalence of these extended, relatively undisturbed disks challenges the merging scenario as the main mode of massive galaxy formation.

We investigate how star-forming galaxies typically assemble their masses at high redshift. Using the deep multi-wavelength coverage of the GOODS dataset, we measure stellar mass of a large sample of star-forming galaxies at z~4 and 5, and make a robust determination of stellar mass function (SMF). We report a broad correlation between stellar mass and UV luminosity, such that more UV-luminous galaxies are, on average, more massive. However, the correlation has a substantial intrinsic scatter evidenced by a non-negligible number of UV-faint but massive galaxies. Furthermore, the low-mass end of the SMF does not rise as steeply as the UV luminosity function (alpha_UVLF} -(1.7-1.8), alpha_SMF -(1.3-1.4)) of the same galaxies. In a smooth formation scenario where star formation (SF) is sustained at the observed rates for a long time, these galaxies would have accumulated more mass (by a factor of ~3) than observed and therefore the SMF would mirror more closely that of the UVLF. The relatively shallow slope of the SMF is due to the fact that many of the UV-selected galaxies are not massive enough, and therefore are too faint in their rest-frame optical bands, to be detected in the current observations. Our results favor a episodic formation history in which SF in low-mass galaxies vary significantly over time, a scenario favored by galaxy clustering. Our findings for the UV-faint galaxies are in contrast with those found for more UV-luminous galaxies, which exhibit tighter SFR-M_star correlations. The discrepancy may suggest that galaxies at different luminosities may have different evolutionary paths.

A new set of color selection criteria (VJL) analogous with the BzK method is designed to select both star-forming galaxies (SFGs) and passively-evolving galaxies (PEGs) at 2.3<z<3.5 by using rest-frame UV--optical (V-J vs. J-L) colors. The criteria are thoroughly tested with theoretical stellar population synthesis models and real galaxies with spectroscopic redshifts to evaluate their efficiency and contamination. We apply the well-tested VJL criteria to the HST/WFC3 Early Release Science field and study the physical properties of selected galaxies. The redshift distribution of selected SFGs peaks at z~2.7, slightly lower than that of Lyman Break Galaxies at z~3. Comparing the observed mid-infrared fluxes of selected galaxies with the prediction of pure stellar emission, we find that our VJL method is effective at selecting massive dusty SFGs that are missed by the Lyman Break Technique. About half of the star formation in massive (M_{star}>10^{10}M_{Sun}) galaxies at 2.3<z<3.5 is contributed by dusty (extinction E(B-V)>0.4) SFGs, which however, only account for ~20% of the number density of massive SFGs. We also use the mid-infrared fluxes to clean our PEG sample, and find that galaxy size can be used as a secondary criterion to effectively eliminate the contamination of dusty SFGs. The redshift distribution of the cleaned PEG sample peaks at z~2.5. We find 6 PEG candidates at z>3 and discuss possible methods to distinguish them from dusty contamination. We conclude that at least part of our candidates are real PEGs at z~3, implying that this type of galaxies began to form their stars at z>5. We measure the integrated stellar mass density of PEGs at z~2.5 and set constraints on it at z>3. We find that the integrated stellar mass density grows by at least about factor of 10 in 1 Gyr at 3<z<5 and by another factor of 10 in next 3.5 Gyr (1<z<3).

We have analysed a sample of 25 extremely red H-[4.5]>4 galaxies, selected using 4.5 micron data from the Spitzer SEDS survey and deep H-band data from the Hubble Space Telescope (HST) CANDELS survey, over ~180 square arcmin of the UKIDSS Ultra Deep Survey (UDS) field. Our aim is to investigate the nature of this rare population of mid-infrared (mid-IR) sources that display such extreme near-to-mid-IR colours. Using up to 17-band photometry (U through 8.0 microns), we have studied in detail their spectral energy distributions, including possible degeneracies in the photometric redshift/internal extinction (zphot-Av) plane. Our sample appears to include sources of very different nature. Between 45% and 75% of them are dust-obscured, massive galaxies at 3<zphot<5. All of the 24 micron-detected sources in our sample are in this category. Two of these have S(24 micron)>300 microJy, which at 3<zphot<5 suggests that they probably host a dust-obscured active galactic nucleus (AGN). Our sample also contains four highly obscured (Av>5) sources at zphot<1. Finally, we analyse in detail two zphot~6 galaxy candidates, and discuss their plausibility and implications. Overall, our red galaxy sample contains the tip of the iceberg of a larger population of z>3 galaxies to be discovered with the future James Webb Space Telescope.

The spectral energy distribution (SED) of a galaxy contains information on the galaxy's physical properties, and multi-wavelength observations are needed in order to measure these properties via SED fitting. In planning these surveys, optimization of the resources is essential. The Fisher Matrix formalism can be used to quickly determine the best possible experimental setup to achieve the desired constraints on the SED fitting parameters. However, because it relies on the assumption of a Gaussian likelihood function, it is in general less accurate than other slower techniques that reconstruct the probability distribution function (PDF) from the direct comparison between models and data. We compare the uncertainties on SED fitting parameters predicted by the Fisher Matrix to the ones obtained using the more thorough PDF fitting techniques. We use both simulated spectra and real data, and consider a large variety of target galaxies differing in redshift, mass, age, star formation history, dust content, and wavelength coverage. We find that the uncertainties reported by the two methods agree within a factor of two in the vast majority (~ 90%) of cases. If the age determination is uncertain, the top-hat prior in age used in PDF fitting to prevent each galaxy from being older than the Universe needs to be incorporated in the Fisher Matrix, at least approximately, before the two methods can be properly compared. We conclude that the Fisher Matrix is a useful tool for astronomical survey design.

We present Hubble Space Telescope Wide Field Camera 3 slitless grism spectroscopy of 28 emission-line galaxies at z~2, in the GOODS-S region of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). The high sensitivity of these grism observations, with 1-sigma detections of emission lines to f > 2.5x10^{-18} erg/s/cm^2, means that the galaxies in the sample are typically ~7 times less massive (median M_* = 10^{9.5} M_sun) than previously studied z~2 emission-line galaxies. Despite their lower mass, the galaxies have OIII/Hb ratios which are very similar to previously studied z~2 galaxies and much higher than the typical emission-line ratios of local galaxies. The WFC3 grism allows for unique studies of spatial gradients in emission lines, and we stack the two-dimensional spectra of the galaxies for this purpose. In the stacked data the OIII emission line is more spatially concentrated than the Hb emission line with 98.1 confidence. We additionally stack the X-ray data (all sources are individually undetected), and find that the average L(OIII)/L(0.5-10 keV) ratio is intermediate between typical z~0 obscured active galaxies and star-forming galaxies. Together the compactness of the stacked OIII spatial profile and the stacked X-ray data suggest that at least some of these low-mass, low-metallicity galaxies harbor weak active galactic nuclei.

We use the spatially-resolved, multi-band photometry in the GOODS South field acquired by the CANDELS project to constrain the nature of candidate Lyman continuum (LyC) emitters at redshift z~3.7 identified using ultra-deep imaging below the Lyman limit (1-sigma limit of ~30 AB in a 2" diameter aperture). In 18 candidates, out of a sample of 19 with flux detected at >3-sigma level, the light centroid of the candidate LyC emission is offset from that of the LBG by up to 1.5". We fit the SED of the LyC candidates to spectral population synthesis models to measure photometric redshifts and the stellar population parameters. We also discuss the differences in the UV colors between the LBG and the LyC candidates, and how to estimate the escape fraction of ionizing radiation (f_esc) in cases, like in most of our galaxies, where the LyC emission is spatially offset from the host galaxy. In all but one case we conclude that the candidate LyC emission is most likely due to lower redshift interlopers. Based on these findings, we argue that the majority of similar measurements reported in the literature need further investigation before it can be firmly concluded that LyC emission is detected. Our only surviving LyC candidate is a LBG at z=3.795, which shows the bluest (B-V) color among LBGs at similar redshift, a stellar mass of M~2 x 10^9 Msun, weak interstellar absorption lines and a flat UV spectral slope with no Lya in emission. We estimate its f_esc to be in the range 25%-100%, depending on the dust and intergalactic attenuation.

We report the discovery of a Type Ia supernova (SNIa) at redshift z=1.55 with the infrared detector of the Wide Field Camera 3 (WFC3-IR) on the Hubble Space Telescope (HST). This object was discovered in CANDELS imaging data of the Hubble Ultra Deep Field, and followed as part of the CANDELS+CLASH Supernova project, comprising the SN search components from those two HST multi-cycle treasury programs. This is the highest redshift SNIa with direct spectroscopic evidence for classification. It is also the first SN Ia at z>1 found and followed in the infrared, providing a full light curve in rest-frame optical bands. The classification and redshift are securely defined from a combination of multi-band and multi-epoch photometry of the SN, ground-based spectroscopy of the host galaxy, and WFC3-IR grism spectroscopy of both the SN and host. This object is the first of a projected sample at z>1.5 that will be discovered by the CANDELS and CLASH programs. The full CANDELS+CLASH SN Ia sample will enable unique tests for evolutionary effects that could arise due to differences in SN Ia progenitor systems as a function of redshift. This high-z sample will also allow measurement of the SN Ia rate out to z~2, providing a complementary constraint on SN Ia progenitor models.

PDS 144 is a pair of Herbig Ae stars that are separated by 5.35" on the sky. It has previously been shown to have an A2Ve Herbig Ae star viewed at 83\circ inclination as its northern member and an A5Ve Herbig Ae star as its southern member. Direct imagery revealed a disk occulting PDS 144 N - the first edge-on disk observed around a Herbig Ae star. The lack of an obvious disk in direct imagery suggested PDS 144 S might be viewed face-on or not physically associated with PDS 144 N. Multi-epoch Hubble Space Telescope imagery of PDS 144 with a 5 year baseline demonstrates PDS 144 N & S are comoving and have a common proper motion with TYC 6782-878-1. TYC 6782-878-1 has previously been identified as a member of Upper Sco sub-association A at d = 145 \pm 2 pc with an age of 5-10 Myr. Ground-based imagery reveals jets and a string of Herbig-Haro knots extending 13' (possibly further) which are aligned to within 7\circ \pm 6\circ on the sky. By combining proper motion data and the absence of a dark mid-plane with radial velocity data, we measure the inclination of PDS 144 S to be i = 73\circ \pm 7\circ. The radial velocity of the jets from PDS 144 N & S indicates they, and therefore their disks, are misaligned by 25\circ \pm 9\circ. This degree of misalignment is similar to that seen in T Tauri wide binaries.

This paper describes the Hubble Space Telescope imaging data products and data reduction procedures for the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). This survey is designed to document the evolution of galaxies and black holes at $z\sim1.5-8$, and to study Type Ia SNe beyond $z>1.5$. Five premier multi-wavelength sky regions are selected, each with extensive multiwavelength observations. The primary CANDELS data consist of imaging obtained in the Wide Field Camera 3 / infrared channel (WFC3/IR) and UVIS channel, along with the Advanced Camera for Surveys (ACS). The CANDELS/Deep survey covers \sim125 square arcminutes within GOODS-N and GOODS-S, while the remainder consists of the CANDELS/Wide survey, achieving a total of \sim800 square arcminutes across GOODS and three additional fields (EGS, COSMOS, and UDS). We summarize the observational aspects of the survey as motivated by the scientific goals and present a detailed description of the data reduction procedures and products from the survey. Our data reduction methods utilize the most up to date calibration files and image combination procedures. We have paid special attention to correcting a range of instrumental effects, including CTE degradation for ACS, removal of electronic bias-striping present in ACS data after SM4, and persistence effects and other artifacts in WFC3/IR. For each field, we release mosaics for individual epochs and eventual mosaics containing data from all epochs combined, to facilitate photometric variability studies and the deepest possible photometry. A more detailed overview of the science goals and observational design of the survey are presented in a companion paper.

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) is designed to document the first third of galactic evolution, over the approximate redshift (z) range 8--1.5. It will image >250,000 distant galaxies using three separate cameras on the Hubble Space Telescope, from the mid-ultraviolet to the near-infrared, and will find and measure Type Ia supernovae at z>1.5 to test their accuracy as standardizable candles for cosmology. Five premier multi-wavelength sky regions are selected, each with extensive ancillary data. The use of five widely separated fields mitigates cosmic variance and yields statistically robust and complete samples of galaxies down to a stellar mass of 10^9 M_\odot to z \approx 2, reaching the knee of the ultraviolet luminosity function (UVLF) of galaxies to z \approx 8. The survey covers approximately 800 arcmin^2 and is divided into two parts. The CANDELS/Deep survey (5\sigma\ point-source limit H=27.7 mag) covers \sim 125 arcmin^2 within GOODS-N and GOODS-S. The CANDELS/Wide survey includes GOODS and three additional fields (EGS, COSMOS, and UDS) and covers the full area to a 5\sigma\ point-source limit of H \gtrsim 27.0 mag. Together with the Hubble Ultra Deep Fields, the strategy creates a three-tiered "wedding cake" approach that has proven efficient for extragalactic surveys. Data from the survey are nonproprietary and are useful for a wide variety of science investigations. In this paper, we describe the basic motivations for the survey, the CANDELS team science goals and the resulting observational requirements, the field selection and geometry, and the observing design. The Hubble data processing and products are described in a companion paper.

We present results from a search for gravitationally lensed galaxies present in the Hubble Space Telescope (HST)/Wide Field Camera-3 (WFC3) images of the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). We present one bona fide lens system in UDS and two compact lens candidates in the GOODS-S field. The lensing system in UDS involves two background galaxies, one at z=1.847 lensed to an arc and a counterimage, and the second at a photometric redshift of z=2.32^{+0.10}_{-0.06} lensed to a double image. We reconstruct the lensed sources in the source plane and find in each of the two cases the sources can be separated to a pair of galaxies. The sources responsible for the arc are compact with effective radii of 0.3 to 0.4 kpc in WFC3 J_{125}-band and a total stellar mass and a star-formation rate of 2.1_{-0.4}^{+2.4} times 10^7 M_sun and 2.3_{-1.7}^{+ 0.6} M_sun yr^{-1}, respectively.The abnormally high H_{160}-band flux of this source is likely due to OIII emission lines with a rest-frame equivalent width about 700 Angstroms for OIII 5007 Angstroms. The sources responsible for the double image have corresponding values of about 0.4 to 0.5 kpc, 1.4_{-0.8}^{+1.9} times 10^9 M_sun, and 8.7_{-7.0}^{+11.1} M_sun yr^{-1}. Once completed CANDELS is expected to contain about 15 lensing systems and will allow statistical studies on both lensing mass profiles and z ~ 2 lensed galaxies.

We present the recent merger history of massive galaxies in a spectroscopically-confirmed proto-cluster at z=1.62. Using HST WFC3 near-infrared imaging from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS), we select cluster galaxies and z ~ 1.6 field galaxies with M_star >= 3 x 10^10 M_sun, and determine the frequency of double nuclei or close companions with projected separations less than 20 kpc co-moving and stellar mass ratios between 1:1 and roughly 10:1. We find that four out of five spectroscopically-confirmed massive proto-cluster galaxies have double nuclei, and 42 +13/-25 % of all M_star >= 3 x 10^10 M_sun cluster candidates are either in close pair systems or have double nuclei. In contrast, only 4.5 +/- 2.6% of the field galaxies are in close pair/double nuclei systems. The implied merger rate per massive galaxy in the proto-cluster is 3-10 times higher than the merger rate of massive field galaxies at z ~ 1.6, depending upon the assumed mass ratios. Close pairs in the cluster have minor merger stellar mass ratios (M_primary:M_satellite ~ 6:1), while the field pairs are typically major mergers with stellar mass ratios between 1:1 and 4:1. At least half of the cluster mergers are dissipationless, as indicated by their red colors and low 24 micron fluxes. Two of the double-nucleated cluster members have X-ray detected AGN with L_x > 10^43 erg/s, and are strong candidates for dual or offset super-massive black holes. We conclude that the massive z = 1.62 proto-cluster galaxies are undergoing accelerated assembly relative to the field population, and discuss the implications for galaxy evolution in proto-cluster environments.

We identify an abundant population of extreme emission line galaxies (EELGs) at redshift z~1.7 in the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS) imaging from Hubble Space Telescope/Wide Field Camera 3 (HST/WFC3). 69 EELG candidates are selected by the large contribution of exceptionally bright emission lines to their near-infrared broad-band magnitudes. Supported by spectroscopic confirmation of strong [OIII] emission lines -- with rest-frame equivalent widths ~1000\AA -- in the four candidates that have HST/WFC3 grism observations, we conclude that these objects are galaxies with 10^8 Msol in stellar mass, undergoing an enormous starburst phase with M_*/(dM_*/dt) of only ~15 Myr. These bursts may cause outflows that are strong enough to produce cored dark matter profiles in low-mass galaxies. The individual star formation rates and the co-moving number density (3.7x10^-4 Mpc^-3) can produce in ~4 Gyr much of the stellar mass density that is presently contained in 10^8-10^9 Msol dwarf galaxies. Therefore, our observations provide a strong indication that many or even most of the stars in present-day dwarf galaxies formed in strong, short-lived bursts, mostly at z>1.

Using HST/WFC3 imaging taken as part of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS), we examine the role that major galaxy mergers play in triggering active galactic nuclei (AGN) activity at z~2. Our sample consists of 72 moderate-luminosity (Lx ~ 1E42-1E44 erg/s) AGN at 1.5<z<2.5 that are selected using the 4 Msec Chandra observations in the Chandra Deep Field South, the deepest X-ray observations to date. Employing visual classifications, we have analyzed the rest-frame optical morphologies of the AGN host galaxies and compared them to a mass-matched control sample of 216 non-active galaxies at the same redshift. We find that most of the AGN reside in disk galaxies (51.4%), while a smaller percentage are found in spheroids (27.8%). Roughly 16.7% of the AGN hosts have highly disturbed morphologies and appear to be involved in a major merger or interaction, while most of the hosts (55.6%) appear relatively relaxed and undisturbed. These fractions are statistically consistent with the fraction of control galaxies that show similar morphological disturbances. These results suggest that the hosts of moderate-luminosity AGN are no more likely to be involved in an ongoing merger or interaction relative to non-active galaxies of similar mass at z~2. The high disk fraction observed among the AGN hosts also appears to be at odds with predictions that merger-driven accretion should be the dominant AGN fueling mode at z~2, even at moderate X-ray luminosities. Although we cannot rule out that minor mergers are responsible for triggering these systems, the presence of a large population of relatively undisturbed disk-like hosts suggests that secular processes play a greater role in fueling AGN activity at z~2 than previously thought.

We report on the evolution of the number density and size of early-type galaxies from z~2 to z~0. We select a sample of 563 massive (M>10^{10} Msun), passively evolving (SSFR<10^{-2} Gyr^{-1}) and morphologically spheroidal galaxies at 0<z<2.5, using the panchromatic photometry and spectroscopic redshifts available in the GOODS fields. We combine ACS and WFC3 HST images to study the morphology of our galaxies in their optical rest-frame in the whole 0<z<2.5 range. We find that throughout the explored redshift range the passive galaxies selected with our criteria have weak morphological K-correction, with size being slightly smaller in the optical than in the UV rest-frame (by ~20 and ~10 at z>1.2 and z<1.2, respectively). We measure a significant evolution of the mass-size relation of early-type galaxies, with the fractional increment that is almost independent on the stellar mass. Early-type galaxies (ETGs) formed at z>1 appear to be preferentially small, and the evolution of the mass-size relation at z<1 is driven by both the continuous size growth of the compact galaxies and the appearance of new ETGs with large sizes. We also find that the number density of all passive early-type galaxies increases rapidly, by a factor of 5, from z~2 to z~1, and then more mildly by another factor of 1.5 from z~1 to z~0. We interpret these results as the evidence that the bulk of the ETGs are formed at 1<z<3 through a mechanism that leaves very compact remnants. At z<1 the compact ETGs grow gradually in size, becoming normal size galaxies, and at the same time new ETGs with normal-large sizes are formed.

We report the detection of color gradients in six massive (stellar mass > 10^{10} M_{sun}) and passively evolving (specific SFR < 10^{-11}/yr) galaxies at redshift 1.3<z<2.5 identified in the HUDF using HST ACS and WFC3/IR images. After matching different PSFs, we obtain color maps and multi-band optical/near-IR photometry (BVizYJH) in concentric annuli, from the smallest resolved radial (~1.7 kpc) up to several times the H-band effective radius. We find that the inner regions of these galaxies have redder rest-frame UV--optical colors than the outer parts. The slopes of the color gradients mildly depend on the overall dust obscuration and rest-frame (U-V) color, with more obscured or redder galaxies having steeper color gradients. The z~2 color gradients are also steeper than those of local early-types. The gradient of a single parameter (age, extinction or metallicity) cannot fully explain the observed color gradients. Fitting spatially resolved HST seven-band photometry to stellar population synthesis models, we find that, regardless of assumptions for metallicity gradient, the redder inner regions of the galaxies have slightly higher dust obscuration than the bluer outer regions, although the magnitude depends on the assumed extinction law. The derived age gradient depends on the assumptions for metallicity gradient. We discuss the implications of a number of assumptions for metallicity gradient on the formation and evolution of these galaxies. We find that the evolution of the mass--size relationship from z~2 to z~0 cannot be driven by in--situ extended star formation, implying that accretion or merger is mostly responsible for the evolution. The lack of a correlation between color gradient and stellar mass argues against the metallicity gradient predicted by the monolithic collapse, which would require significant major mergers to evolve into the one observed at z~0. (Abridged)

We compute black hole masses and bolometric luminosities for 87 obscured AGN in the redshift range 0.25 < z < 1.25, selected from the GOODS deep multi-wavelength survey fields via their X-ray emission. We fit the optical images and obtain morphological parameters for the host galaxy, separating the galaxy from its central point source, thereby obtaining a four-band optical SED for each active nucleus. We calculate bolometric luminosities for these AGN by reddening a normalized mean SED of GOODS broad-line AGN to match the observed central point-source SED of each obscured AGN. This estimate of Lbol has a smaller spread than simple bolometric corrections to the X-ray luminosity or direct integration of the observed multi-wavelength SED, suggesting it is a better measure. We estimate central black hole masses from the bulge luminosities. The black hole masses span a wide range, 7 x 10^6 M_sun to 6 x 10^9 M_sun; the median black hole mass is 5 x 10^8 M_sun. The majority of these AGN have L/L_Edd < 0.01, and we detect no significant evolution of the mean Eddington ratio to z = 1.25. This implies that the bulk of black hole growth in these obscured AGN must have occurred at z > 1 and that we are observing these AGN in a slow- or no-growth state.

We present spectroscopy of 76 emission-line galaxies (ELGs) in CDF-S taken with the LDSS3 spectrograph on Magellan Telescope. These galaxies are selected to have emission lines with ACS grism data in the Hubble Space Telescope Probing Evolution and Reionization Spectroscopically (PEARS) grism Survey. The ACS grism spectra cover the wavelength range 6000-9700 \AA\ and most PEARS grism redshifts are based on a single emission line + photometric redshifts from broad-band colors; the Magellan spectra cover a wavelength range from 4000 {\AA} to 9000 {\AA}, and provide a check on redshifts derived from PEARS data. We find an accuracy of $\sigma_z$ = 0.006 for the ACS grism redshifts with only one catastrophic outlier. We probe for AGN in our sample via several different methods. In total we find 7 AGNs and AGN candidates out of 76 galaxies. Two AGNs are identified from the X-ray full-band luminosity, $L_{X-ray,FB}>10^{43}$ erg$\;$s$^{-1}$, the line widths and the power-law continuum spectra. Two unobscured faint AGN candidates are identified from the X-ray full-band luminosity $L_{X-ray,FB}\sim10^{41}$ erg$\;$s$^{-1}$, the hardness ratio and the column density, and the emission-line and X-ray derived SFRs. Two candidates are classified based on the line ratio of [NII]\lambda6584/H$\alpha$ versus [OIII]$\lambda$5007/H$\beta$ (BPT diagram), which are between the empirical and theoretical demarcation curves, i.e, the transition region from star-forming galaxies to AGNs. One AGN candidate is identified from the high-ionization emission line HeII{\AA}4686.

We use the deepest and the most comprehensive photometric data currently available for GOODS-South galaxies to measure their photometric redshifts. The photometry includes VLT/VIMOS (U-band), HST/ACS (F435W, F606W, F775W, and F850LP bands), VLT/ISAAC (J-, H-, and Ks-bands), and four Spitzer/IRAC channels (3.6, 4.5, 5.8, and 8.0 micron). The catalog is selected in the z-band (F850LP) and photometry in each band is carried out using the recently completed TFIT algorithm, which performs PSF matched photometry uniformly across different instruments and filters, despite large variations in PSFs and pixel scales. Photometric redshifts are derived using the GOODZ code, which is based on the template fitting method using priors. The code also implements "training" of the template SED set, using available spectroscopic redshifts in order to minimize systematic differences between the templates and the SEDs of the observed galaxies. Our final catalog covers an area of 153 sq. arcmin and includes photometric redshifts for a total of 32,505 objects. The scatter between our estimated photometric and spectroscopic redshifts is sigma=0.040 with 3.7% outliers to the full z-band depth of our catalog, decreasing to sigma=0.039 and 2.1% outliers at a magnitude limit m(z)<24.5. This is consistent with the best results previously published for GOODS-S galaxies, however, the present catalog is the deepest yet available and provides photometric redshifts for significantly more objects to deeper flux limits and higher redshifts than earlier works. Furthermore, we show that the photometric redshifts estimated here for galaxies selected as dropouts are consistent with those expected based on the Lyman break technique.

We present the details and early results from a deep near-infrared survey utilising the NICMOS instrument on the Hubble Space Telescope centred around massive M_* > 10^11 M_0 galaxies at 1.7 < z < 2.9 found within the Great Observatories Origins Deep Survey (GOODS) fields. The GOODS NICMOS Survey (GNS) was designed to obtain deep F160W (H-band) imaging of 80 of these massive galaxies, as well as other colour selected objects such as Lyman-break drop-outs, BzK objects, Distant Red Galaxies, EROs, Spitzer Selected EROs, BX/BM galaxies, as well as sub-mm galaxies. We present in this paper details of the observations, our sample selection, as well as a description of features of the massive galaxies found within our survey fields. This includes: photometric redshifts, rest-frame colours, and stellar masses. We furthermore provide an analysis of the selection methods for finding massive galaxies at high redshifts, including colour selection, and how galaxy populations selected through different methods overlap. We find that a single colour selection method cannot locate all of the massive galaxies, with no one method finding more than 70 percent. We however find that the combination of these colour methods finds nearly all the massive galaxies, as selected by photometric redshifts with the exception of apparently rare blue massive galaxies. By investigating the rest-frame (U-B) vs. M_B diagram for these galaxies we furthermore show that there exists a bimodality in colour-magnitude space at z < 2, driven by stellar mass, such that the most massive galaxies are systematically red up to z~2.5, while lower mass galaxies tend to be blue. We also discuss the number densities for galaxies with stellar masses M_* > 10^11 M_0, whereby we find an increase of a factor of eight between z = 3 and z = 1.5, demonstrating that this is an epoch when massive galaxies establish most of their mass.

Variability is a property shared by practically all AGN. This makes variability selection a possible technique for identifying AGN. Given that variability selection makes no prior assumption about spectral properties, it is a powerful technique for detecting both low-luminosity AGN in which the host galaxy emission is dominating and AGN with unusual spectral properties. In this paper, we will discuss and test different statistical methods for the detection of variability in sparsely sampled data that allow full control over the false positive rates. We will apply these methods to the GOODS North and South fields and present a catalog of variable sources in the z band in both GOODS fields. Out of 11931 objects checked, we find 155 variable sources at a significance level of 99.9%, corresponding to about 1.3% of all objects. After rejection of stars and supernovae, 139 variability selected AGN remain. Their magnitudes reach down as faint as 25.5 mag in z. Spectroscopic redshifts are available for 22 of the variability selected AGN, ranging from 0.046 to 3.7. The absolute magnitudes in the rest-frame z-band range from ~ -18 to -24, reaching substantially fainter than the typical luminosities probed by traditional X-ray and spectroscopic AGN selection in these fields. Therefore, this is a powerful technique for future exploration of the evolution of the faint end of the AGN luminosity function up to high redshifts.

We present the highest redshift detections of resolved Lyman alpha emission, using Hubble Space Telescope/ACS F658N narrowband-imaging data taken in parallel with the Wide Field Camera 3 Early Release Science program in the GOODS CDF-S. We detect Lyman alpha emission from three spectroscopically confirmed z = 4.4 Lyman alpha emitting galaxies (LAEs), more than doubling the sample of LAEs with resolved Lyman alpha emission. Comparing the light distribution between the rest-frame ultraviolet continuum and narrowband images, we investigate the escape of Lyman alpha photons at high redshift. While our data do not support a positional offset between the Lyman alpha and rest-frame ultraviolet (UV) continuum emission, the half-light radii in two out of the three galaxies are significantly larger in Lyman alpha than in the rest-frame UV continuum. This result is confirmed when comparing object sizes in a stack of all objects in both bands. Additionally, the narrowband flux detected with HST is significantly less than observed in similar filters from the ground. These results together imply that the Lyman alpha emission is not strictly confined to its indigenous star-forming regions. Rather, the Lyman alpha emission is more extended, with the missing HST flux likely existing in a diffuse outer halo. This suggests that the radiative transfer of Lyman alpha photons in high-redshift LAEs is complicated, with the interstellar-medium geometry and/or outflows playing a significant role in galaxies at these redshifts.

We use the multi-epoch, mid-infrared Spitzer Deep, Wide-Field Survey to investigate the variability of 474,179 objects in 8.1 deg^2 of the NDWFS Bootes field. We perform a Difference Image Analysis of the four available epochs between 2004 and 2008, focusing on the deeper 3.6 and 4.5 micron bands. We find that 1.1% of the studied sample meet our standard selection criteria for being classed as a variable source. We require that the 3.6 and 4.5 micron light-curves are strongly correlated (r>0.8) and that their joint variance exceeds that for all sources with the same magnitude by 2 sigma. We then examine the mid-IR colors of the variable sources and match them with X-ray sources from the XBootes survey, radio catalogs, 24 micron-selected AGN candidates, and spectroscopically identified AGNs from the AGN and Galaxy Evolution Survey (AGES). Based on their mid-IR colors, most of the variable sources are AGNs (76%), with smaller contributions from stars (11%), galaxies (6%), and unclassified objects. Most of the stellar, galaxy and unclassified sources are false positives. For our standard selection criteria, 11-12% of the mid-IR counterparts to X-ray sources, 24 micron-selected AGN candidates and spectroscopically identified AGNs show variability. Mid-IR AGN variability can be well described by a single power-law structure function with a power-law index of 0.5 at both 3.6 and 4.5 microns, and an amplitude of 0.1 mag on rest-frame time scales of 2 years. The variability amplitude is higher for shorter rest-frame wavelengths and lower luminosities. (Abridged)

We present near-IR images, obtained with the Hubble Space Telescope (HST) and the WFC3/IR camera, of six passive and massive galaxies at redshift 1.3<z<2.4 (SSFR<10^{-2} Gyr^{-1}; stellar mass M~10^{11} M_{sun}), selected from the Great Observatories Origins Deep Survey (GOODS). These images, which have a spatial resolution of ~1.5 kpc, provide the deepest view of the optical rest-frame morphology of such systems to date. We find that the light profile of these galaxies is regular and well described by a Sersic model with index typical of today's spheroids. Their size, however, is generally much smaller than today's early types of similar stellar mass, with four out of six galaxies having r_e ~ 1 kpc or less, in quantitative agreement with previous similar measures made at rest-frame UV wavelengths. The images reach limiting surface brightness mu~26.5 mag arcsec^{-2} in the F160W bandpass; yet, there is no evidence of a faint halo in the galaxies of our sample, even in their stacked image. We also find that these galaxies have very weak "morphological k-correction" between the rest-frame UV (from the ACS z-band), and the rest--frame optical (WFC3 H-band): the Sersic index, physical size and overall morphology are independent or only mildly dependent on the wavelength, within the errors.

Using HST/ACS slitless grism spectra from the PEARS program, we study the stellar populations of morphologically selected early-type galaxies in the GOODS North and South fields. The sample - extracted from a visual classification of the (v2.0) HST/ACS images and restricted to redshifts z>0.4 - comprises 228 galaxies (F775W<24 ABmag) out to z~1.3 over 320 arcmin2, with a median redshift zM=0.75. This work significantly increases our previous sample from the GRAPES survey in the HUDF (18 galaxies over ~11 arcmin2; Pasquali et al. 2006b). The grism data allow us to separate the sample into `red' and `blue' spectra, with the latter comprising 15% of the total. Three different grids of models parameterising the star formation history are used to fit the low-resolution spectra. Over the redshift range of the sample - corresponding to a cosmic age between 5 and 10 Gyr - we find a strong correlation between stellar mass and average age, whereas the **spread** of ages (defined by the RMS of the distribution) is roughly ~1 Gyr and independent of stellar mass. The best-fit parameters suggest it is formation epoch and not formation timescale, that best correlates with mass in early-type galaxies. This result - along with the recently observed lack of evolution of the number density of massive galaxies - motivates the need for a channel of (massive) galaxy formation bypassing any phase in the blue cloud, as suggested by the simulations of Dekel et al. (2009).

We present results of a search for emission-line galaxies in the Southern Fields of the Hubble Space Telescope PEARS (Probing Evolution And Reionization Spectroscopically) grism survey. The PEARS South Fields consist of five ACS pointings (including the Hubble Ultra Deep Field) with the G800L grism for a total of 120 orbits, revealing thousands of faint object spectra in the GOODS-South region of the sky. Emission-line galaxies (ELGs) are one subset of objects that are prevalent among the grism spectra. Using a 2-dimensional detection and extraction procedure, we find 320 emission lines orginating from 226 galaxy "knots'' within 192 individual galaxies. Line identification results in 118 new grism-spectroscopic redshifts for galaxies in the GOODS-South Field. We measure emission line fluxes using standard Gaussian fitting techniques. At the resolution of the grism data, the H-beta and [OIII] doublet are blended. However, by fitting two Gaussian components to the H-beta and [OIII] features, we find that many of the PEARS ELGs have high [OIII]/H-beta ratios compared to other galaxy samples of comparable luminosities. The star-formation rates (SFRs) of the ELGs are presented, as well as a sample of distinct giant star-forming regions at z~0.1-0.5 across individual galaxies. We find that the radial distances of these HII regions in general reside near the galaxies' optical continuum half-light radii, similar to those of giant HII regions in local galaxies.

We apply Bayesian statistics with prior probabilities of galaxy surface luminosity (SL) to improve photometric redshifts. We apply the method to a sample of 1266 galaxies with spectroscopic redshifts in the GOODS North and South fields at 0.1 < z < 2.0. We start with spectrophotometric redshifts (SPZs) based on Probing Evolution and Reionization Spectroscopically grism spectra, which cover a wavelength range of 6000-9000A, combined with (U)BViz(JHK) broadband photometry in the GOODS fields. The accuracy of SPZ redshifts is estimated to be \sigma (\Delta(z))=0.035 with an systematic offset of -0.026, where \Delta(z)=\Delta z / (1+z), for galaxies in redshift range of 0.5 < z < 1.25. The addition of the SL prior probability helps break the degeneracy of SPZ redshifts between low redshift 4000 A break galaxies and high-redshift Lyman break galaxies which are mostly catastrophic outliers. For the 1138 galaxies at z < 1.6, the fraction of galaxies with redshift deviation \Delta (z) > 0.2 is reduced from 15.0% to 10.4%, while the rms scatter of the fractional redshift error does not change much.

The Spitzer Deep, Wide-Field Survey (SDWFS) is a four-epoch infrared survey of ten square degrees in the Bootes field of the NOAO Deep Wide-Field Survey using the IRAC instrument on the Spitzer Space Telescope. SDWFS, a Cycle four Spitzer Legacy project, occupies a unique position in the area-depth survey space defined by other Spitzer surveys. The four epochs that make up SDWFS permit -- for the first time -- the selection of infrared-variable and high proper motion objects over a wide field on timescales of years. Because of its large survey volume, SDWFS is sensitive to galaxies out to z~3 with relatively little impact from cosmic variance for all but the richest systems. The SDWFS datasets will thus be especially useful for characterizing galaxy evolution beyond z~1.5. This paper explains the SDWFS observing strategy and data processing, presents the SDWFS mosaics and source catalogs, and discusses some early scientific findings. The publicly-released, full-depth catalogs contain 6.78, 5.23, 1.20, and 0.96 x 10e5 distinct sources detected to the average 5-sigma, 4" diameter, aperture-corrected limits of 19.77, 18.83, 16.50, and 15.82 Vega mag at 3.6, 4.5, 5.8, and 8.0 micron, respectively. The SDWFS number counts and color-color distribution are consistent with other, earlier Spitzer surveys. At the 6 min integration time of the SDWFS IRAC imaging, more than 50% of isolated FIRST radio sources and more than 80% of on-axis XBootes sources are detected out to 8.0 micron. Finally, we present the four highest proper motion IRAC-selected sources identified from the multi-epoch imaging, two of which are likely field brown dwarfs of mid-T spectral class.

We present stellar population modeling results for 10 newly discovered Lyman alpha emitting galaxies (LAEs), as well as four previously known LAEs at z ~ 4.5 in the Chandra Deep Field - South. We fit stellar population models to these objects in order to learn specifically if there exists more than one class of LAE. Past observational and theoretical evidence has shown that while many LAEs appear to be young, they may be much older, with Lyman alpha EWs enhanced due to resonant scattering of Lyman alpha photons in a clumpy interstellar medium (ISM). Our results show a large range of stellar population age (3 - 500 Myr), stellar mass (1.6 x 10^8 - 5.0 x 10^10 Msol) and dust extinction (A_1200 = 0.3 - 4.5 mag), broadly consistent with previous studies. With such a large number of individually analyzed objects, we have looked at the distribution of stellar population ages in LAEs for the first time, and we find a very interesting bimodality, in that our objects are either very young (< 15 Myr) or old (> 450 Myr). This bimodality may be caused by dust, and it could explain the Lyman alpha duty cycle which has been proposed in the literature. We find that eight of the young objects are best fit with a clumpy ISM. We find that dust geometry appears to play a large role in shaping the SEDs that we observe, and that it may be a major factor in the observed Lyman alpha equivalent width distribution in high redshift Lyman alpha galaxies, although other factors (i.e. outflows) may be in play. We conclude that 12 out of our 14 LAEs are dusty star-forming galaxies, with the other two LAEs being evolved galaxies.

We analyze a new gravitational lens, OAC-GL J1223-1239, serendipitously found in a deep I-band image of the Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS). The lens is a L_*, edge-on S0 galaxy at z=0.4656. The gravitational arc has a radius of 0.42 arcsec. We have determined the total mass and the dark matter (DM) fraction within the Einstein radius as a function of the lensed source redshift, which is presently unknown. For z ~ 1.3, which is in the middle of the redshift range plausible for the source according to some external constraints, we find the central velocity dispersion to be ~180 km/s. With this value, close to that obtained by means of the Faber-Jackson relation at the lens redshift, we compute a 30% DM fraction within the Einstein radius (given the uncertainty in the source redshift, the allowed range for the DM fraction is 25-35 % in our lensing model). When compared with the galaxies in the local Universe, the lensing galaxy, OAC-GL J1223-1239 seems to fall in the transition regime between massive DM dominated galaxies and lower-mass, DM deficient systems.

We present the faintest spectroscopically confirmed sample of redshift four and five Lyman break galaxies to date. The sample is based on slitless grism spectra of the Hubble Ultra Deep Field region from the GRAPES (Grism ACS Program for Extragalactic Science) and PEARS (Probing Evolution and Reionization Spectroscopically) projects, using the G800L grism on the HST Advanced Camera for Surveys. We report here confirmations of 39 galaxies, pre-selected as candidate Lyman break galaxies using photometric selection criteria. We compare a "traditional" V-dropout selection to a more liberal one (with V-i > 0.9), and find that the traditional criteria are about 64% complete and 81% reliable. We also study the Lyman alpha emission properties of our sample. We find that Lyman alpha emission is detected in about 1/4 of the sample, and that our broad-band color selected sample includes 55% of previously published line-selected Lyman alpha sources. Finally, we examine our stacked 2D spectra. We demonstrate that strong, spatially extended (arcsecond scale) Lyman alpha emission is not a generic property of these Lyman break galaxies, but that a modest extension of the Lyman alpha photosphere (compared to the starlight) may be present in those galaxies with prominent Lyman alpha emission.

We present spectroscopic observations of 26 galaxies of type E and S0, based on their blue morphologies, located in voids by the study of Grogin & Geller (1999). Measurements of redshift, velocity dispersion, and four Lick line indices, Mgb, Fe5270, Fe5335, and Hbeta with their errors are given for all of these galaxies, along with Hbeta, [OIII], Halpha, and [NII] emission line strengths for a subset of these objects. These sources are brighter than M star for low density regions and tend to be bluer than their counterpart early-type objects in high density regions. Using the models of Thomas et al. (2003) gives metal abundances and ages with a median alpha enhancement, [alpha/Fe] = +0.13, and median metals abundance, [Z/H] = +0.22, values comparable to those found for E and S0 galaxies in clusters, but with a wider spread in [Z/H] towards low values. If the emission line subsample is interpreted as younger, the proportion of young objects is higher than for early-types in higher density regions. There is a significant incidence of sources in the sample with emission lines in their spectra (46% with Hbeta and [OIII] and 69% with Halpha or [NII]) as well as shells and rings in their morphologies (19%). The diagnostic log[NII/Halpha], log[OIII/Hbeta] diagram places 10 of 12 emission line galaxies in or near the star forming and liner region and two among the Seyferts. The Halpha fluxes indicate star formation rates of 0.2 to 1.0 Msun per yr. The percentage of these early-type void galaxies undergoing star formation appears to be higher compared to their cluster counterparts and the range of ages wider.

We report results from a search for massive and evolved galaxies at z>5 in the GOODS southern field. Combining HST ACS, VLT ISAAC and Spitzer IRAC photometric data, we develop a color selection technique to identify candidates for being evolved galaxies at high redshifts. The color selection is primarily based on locating the Balmer-break using the K- and 3.6micron bands. Stellar population synthesis models are fitted to the SEDs of these galaxies to identify the final sample. We find 11 candidates with photometric redshifts in the range 4.9 < z < 6.5, dominated by an old stellar population, with ages 0.2-1.0 Gyr, and stellar masses in the range (0.5 - 5) 10^{11} Msun. The majority of the stars in these galaxies were formed at z > 9. One candidate has a spectroscopically confirmed redshift, in good agreement with our photometric redshift. The galaxies are very compact, with half-light radii in the observed K-band smaller than ~2 kpc. Seven of the 11 candidates are also detected at 24micron with the MIPS instrument on Spitzer. The 24micron emission could be interpreted as PAH emission from a dusty starburst at z~2-3, however, it is also consistent with the presence of an obscured AGN at z>5. We estimate the completeness of the Balmer break galaxy sample to be ~40%. The comoving number density of galaxies with a stellar mass >10^{11} Msun, at an average redshift z=5.2, is 3.9 10^{-5} Mpc^{-3} (no-MIPS sample: 1.4 10^{-5} Mpc^{-3}). The corresponding stellar mass density is 8 10^{6} Msun/Mpc^3 (no-MIPS sample: 6.2 10^6 Msun/Mpc^3).

Equivalent widths (EWs) observed in high-redshift Lyman alpha galaxies could be stronger than the EW intrinsic to the stellar population if dust is present residing in clumps in the inter-stellar medium (ISM). In this scenario, continuum photons could be extinguished while the Lyman alpha photons would be resonantly scattered by the clumps, eventually escaping the galaxy. We investigate this radiative transfer scenario with a new sample of six Lyman alpha galaxy candidates in the GOODS CDF-S, selected at z = 4.4 with ground-based narrow-band imaging obtained at CTIO. Grism spectra from the HST PEARS survey confirm that three objects are at z = 4.4, and that another object contains an active galactic nuclei (AGN). If we assume the other five (non-AGN) objects are at z = 4.4, they have rest-frame EWs from 47 -- 190 A. We present results of stellar population studies of these objects, constraining their rest-frame UV with HST and their rest-frame optical with Spitzer. Out of the four objects which we analyzed, three objects were best-fit to contain stellar populations with ages on the order of 1 Myr and stellar masses from 3 - 10 x 10^8 solar masses, with dust in the amount of A_1200 = 0.9 - 1.8 residing in a quasi-homogeneous distribution. However, one object (with a rest EW ~ 150 A) was best fit by an 800 Myr, 6.6 x 10^9 solar mass stellar population with a smaller amount of dust (A_1200 = 0.4) attenuating the continuum only. In this object, the EW was enhanced ~ 50% due to this dust. This suggests that large EW Lyman alpha galaxies are a diverse population. Preferential extinction of the continuum in a clumpy ISM deserves further investigation as a possible cause of the overabundance of large-EW objects that have been seen in narrow-band surveys in recent years.

Approximately 20-30% of 1.4<z<2.5 galaxies with K<22 (Vega) detected with Spitzer MIPS at 24um show excess mid-IR emission relative to that expected based on the rates of star formation measured from other multiwavelength data.These galaxies also display some near-IR excess in Spitzer IRAC data, with a spectral energy distribution peaking longward of 1.6um in the rest frame, indicating the presence of warm-dust emission usually absent in star forming galaxies. Stacking Chandra data for the mid-IR excess galaxies yields a significant hard X-ray detection at rest-frame energies >6.2 keV. The stacked X-ray spectrum rises steeply at >10 keV, suggesting that these sources host Compton-thick Active Galactic Nuclei (AGNs) with column densities N_H~10^{24} cm^-2 and an average, unobscured X-ray luminosity L_{2-8 keV}~(1-4)x10^43 erg/s. Their sky density(~3200 deg^-2) and space density (~2.6x10^-4 Mpc^-3) are twice those of X-ray detected AGNs at z~2, and much larger than those of previously-known Compton thick sources at similar redshifts. The mid-IR excess galaxies are part of the long sought-after population of distant heavily obscured AGNs predicted by synthesis models of the X-ray background. The fraction of mid-IR excess objects increases with galaxy mass, reaching ~50-60% for M~10^11 M_sun, an effect likely connected with downsizing in galaxy formation. The ratio of theinferred black hole growth rate from these Compton-thick sources to the global star formation rate at z=2 is similar to the mass ratio of black holes to stars inlocal spheroids, implying concurrent growth of both within the precursors oftoday's massive galaxies.

Cosmic strings were postulated by Kibble in 1976 and, from a theoretical point of view, their existence finds support in modern superstring theories, both in compactification models and in theories with extended additional dimensions. Their eventual discovery would lead to significant advances in both cosmology and fundamental physics. One of the most effective ways to detect cosmic strings is through their lensing signatures which appear to be significantly different from those introduced by standard lenses (id est, compact clumps of matter). In 2003, the discovery of the peculiar object CSL-1 (Sazhin et al.2003) raised the interest of the physics community since its morphology and spectral features strongly argued in favour of it being the first case of gravitational lensing by a cosmic string. In this paper we provide a detailed description of the expected observational effects of a cosmic string and show, by means of simulations, the lensing signatures produced on background galaxies. While high angular resolution images obtained with HST, revealed that CSL-1 is a pair of interacting ellipticals at redshift 0.46, it represents a useful lesson to plan future surveys.

We present an Ultraviolet (UV) selected sample of 268 objects in the two fields of the Great Observatories Origins Deep Survey (GOODS). We used the parallel observations taken with WFPC2 in the U--band (F300W) which covered 88% of the GOODS fields to identify sources and selected only objects with GOODS/ACS counterparts. Spectroscopic redshifts for 95 of these sources are available and we have used the multiwavelength GOODS data to estimate photometric redshifts for the others. Most of the objects are between 0.2<z<0.8. We used the spectral types obtained by the photometric redshift fitting to identify the starburst galaxies. We have also visually checked all objects and looked for tidal effects and nearby companions. We find that (i) 45% of the UV-selected galaxies are starbursts, (ii) nearly 75% of the starbursts have tidal tails or show some peculiarity typical of interactions or mergers, (iii) ~50% have companions within an area of 5x5 arcsec. The UV-selected sample has an average rest-frame M_B=--19.9 +- 0.1. The bluest objects in the sample (U-B < 0.2 and B-V < 0.1) are at 1.1<z<1.9 and have peculiar morphologies that resembles either tadpoles, chains, or double-clump galaxies. Starbursts with tadpole or clump morphology at z=0.8-1.3 have sizes comparable to LBGs and compact Ultraviolet-luminous galaxies (UVLGs).

On January 12 2006, the Hubble Space Telescope observed the peculiar double extragalactic object CSL-1, suspected to be the result of gravitational lensing by a cosmic string. The high resolution image shows that the object is actually a pair of interacting giant elliptical galaxies. In spite of the weird similarities of the energy and light distributions and of the radial velocities of the two components, CSL-1 is not the lensing of an elliptical galaxy by a cosmic string.

We present mid-infrared observations of AGN in the GOODS fields, performed with the Spitzer Space Telescope. These are the deepest infrared and X-ray fields to date and cover a total area of ~0.1 square degrees. AGN are selected on the basis of their hard (2-8 keV) X-ray emission. The median AGN infrared luminosity is at least 10 times larger than the median for normal galaxies with the same redshift distribution, suggesting that the infrared emission is dominated by the central nucleus. The X-ray to infrared luminosity ratios of GOODS AGN, most of which are at 0.5<z<1.5, are similar to the values obtained for AGN in the local Universe. The observed infrared flux distribution has an integral slope of ~1.5 and there are 1000 sources per square degree brighter than ~50 uJy at 3-6 microns. The counts approximately match the predictions of models based on AGN unification, in which the majority of AGN are obscured. This agreement confirms that the faintest X-ray sources, which are dominated by the host galaxy light in the optical, are obscured AGN. Using these Spitzer data, the AGN contribution to the extragalactic infrared background light is calculated by correlating the X-ray and infrared catalogues. This is likely to be a lower limit given that the most obscured AGN are missed in X-rays. We estimate the contribution of AGN missed in X-rays, using a population synthesis model, to be ~45% of the observed AGN contribution, making the AGN contribution to the infrared background at most ~2-10% in the 3-24 micron range, depending on wavelength, lower than most previous estimates. The AGN contribution to the infrared background remains roughly constant with source flux in the IRAC bands but decreases with decreasing flux in the MIPS 24 um band, where the galaxy population becomes more important.

We present results from a search for high-redshift J--band ``dropout'' galaxies in the portion of the GOODS southern field that is covered by extremely deep imaging from the Hubble Ultradeep Field (HUDF).Using observations at optical, near-infrared and mid-infrared wavelengths from the Hubble and Spitzer Space Telescopes and the ESO-VLT, we search for very massive galaxies at high redshifts and find one particularly remarkable candidate. Its spectral energy distribution is consistent with a galaxy at z ~ 6.5 and a stellar mass of 6x10e11 M(sun) (for a Salpeter IMF). We interpret a prominent photometric break between the near-infrared and Spitzer bandpasses as the 3646A Balmer discontinuity. The best-fitting models have low reddening and ages of several hundred Myr, placing the formation of the bulk of the stars at z > 9. Alternative models of dusty galaxies at z ~ 2.5 are possible but provide significantly poorer fits. The object is detected with Spitzer at 24 micron. This emission originats from an obscured active nucleus or star formation. We present optical and near-infrared spectroscopy which has, thus far, failed to detect any spectral features. This helps limit the solution in which the galaxy is a starburst or active galaxy at z ~ 2.5. If the high-redshift interpretation is correct, this object would be an example of a galaxy that formed by a process strongly resembling traditional models of monolithic collapse, in a way which a very large mass of stars formed within a remarkably short period of time, at very high redshift.

We use data from the first epoch of observations with the IRAC/Spitzer for the GOODS to detect and study a collection of LBGs at z ~ 6 to 5 in the HUDF, six of which have spectroscopic confirmation. At these redshifts, IRAC samples rest-frame optical light in the range 0.5 to 0.8 um, where the effects of dust extinction are smaller and the sensitivity to light from evolved stars is greater than at shorter, rest-frame UV observable from the ground or with the HST. As such, it provides useful constraints on the ages and masses of these galaxies' stellar populations. We find that the SEDs for many of these galaxies are best fitted by models of stellar populations with masses of a few x 1e10 M_sun, and with ages of a few hundred Myrs, values quite similar to those derived for typical LBGs at z ~ 3. When the universe was only 1 Gyr old, some galaxies had already formed a mass of stars approaching that of the present-day Milky Way, and that they started forming those stars at z > 7, and in some cases much earlier. We find that the lower limits to the space density for galaxies in this mass range are consistent with predictions from recent hydrodynamic simulations of structure formation in a LCDM universe. All objects in our samples are consistent with having solar metallicity, suggesting that they might have already been significantly polluted by metals. The values for dust reddening derived from the model fitting are low or zero, and we find that some of the galaxies have rest-frame UV colors that are even bluer than those predicted by the stellar population models to which we compare them. These colors might be attributed to the presence of very massive stars (> 100 M_sun), or by weaker intergalactic HI absorption than what is commonly assumed.

We investigate morphological structure parameters and local environments of distant moderate-luminosity active galactic nucleus (AGN) host galaxies in the overlap between the HST/ACS observations of the Great Observatories Origins Deep Survey (GOODS) and the two Chandra Deep Fields. We compute near-neighbor counts and BViz asymmetry (A) and concentration (C) indices for ~35,500 GOODS/ACS galaxies complete to z_850 ~ 26.6, including the resolved hosts of 322 X-ray-selected AGNs. Distributions of (1) z_850 asymmetry for 130 AGN hosts at z_850 < 23 and (2) near-neighbor counts for 173 AGN hosts at z_850 < 24 are both consistent with non-AGN control samples. This implies no close connection between recent galaxy mergers and moderate-luminosity AGN activity out to appreciable look-back times (z < 1.3), approaching the epoch of peak AGN activity in the universe. The distribution of z_850 concentrations for the AGN hosts is offset by +0.5 compared to the non-AGN, a 6.4-sigma discrepancy much larger than can be explained by the possible influence of unresolved emission from the AGN or a circumnuclear starburst. The local universe association between AGN and bulge-dominated galaxies thus persists to substantial look-back time. We discuss implications in the context of the low-redshift supermassive central black hole mass correlation with host galaxy properties, including concentration.

We constrain the X-ray emission properties of Lyman break galaxies (LBGs) at z ~ 3-6 using the 2 Ms Chandra Deep Field-North and 1 Ms Chandra Deep Field-South. Large samples of LBGs were discovered using HST as part of the Great Observatories Origins Deep Survey (GOODS). Deep optical and X-ray imaging over the GOODS fields have allowed us to place the most significant constraints on the X-ray properties of LBGs to date. Mean X-ray properties of 449, 1734, 629, and 247 LBGs with z ~ 3, 4, 5, and 6, respectively, were determined using stacking techniques. When stacked, we detect X-ray emission from LBGs at z ~ 3 (7 sigma) and from an optically bright subset (brightest 25%) of LBGs at z ~ 4 (3 sigma); the latter is the highest redshift detection yet for "normal" galaxies in the X-ray band. The effective exposure times for these stacked observations are 0.7 and 0.5 Gs, respectively. The derived average rest-frame 2.0-8.0 keV luminosities are 1.5 and 1.4 x 10^41 erg s^-1, respectively. X-ray emission from these LBGs is likely due to high mass X-ray binaries (HMXBs) and Type II supernovae; the corresponding star formation rates are 85-240 M_sol yr^-1. The X-ray to B-band mean luminosity ratio (Lx/LB) at z~3 is somewhat elevated with respect to that measured for starburst galaxies in the local Universe. When stacking full samples of LBGs at z ~ 4, 5, and 6 we do not obtain significant detections and derive rest-frame 2.0-8.0 keV luminosity upper limits of 0.9, 2.8, and 7.1 x 10^41 erg s^-1, respectively. These upper limits constrain any widespread AGN activity in these objects to be modest at best. Furthermore, we find that ~ 0.5% of our LBGs from z ~ 3-6 are detected individually in the X-ray band. Continued ...

The deep X-ray, optical, and far-infrared fields that constitute GOODS are sensitive to obscured AGN (N_H>10^{22} cm^{-2}) at the quasar epoch (z~2-3), as well as to unobscured AGN as distant as z~7. Luminous X-ray emission is a sign of accretion onto a supermassive black hole and thus reveals all but the most heavily obscured AGN. We combine X-ray luminosity functions with appropriate spectral energy distributions for AGN to model the X-ray, optical and far-infrared flux distributions of the X-ray sources in the GOODS fields. A simple model based on the unified paradigm for AGN, with ~3 times as many obscured AGN as unobscured, successfully reproduces the z-band flux distributions measured in the deep HST ACS observations on the GOODS North and South fields. This model is also consistent with the observed spectroscopic and photometric redshift distributions once selection effects are considered. The previously reported discrepancy between observed spectroscopic redshift distributions and the predictions of population synthesis models for the X-ray background can be explained by bias against the most heavily obscured AGN generated both by X-ray observations and the identification of sources via optical spectroscopy. We predict the AGN number counts for Spitzer MIPS 24 um and IRAC 3.6-8 um observations in the GOODS fields, which will verify whether most AGN in the early Universe are obscured in the optical. Such AGN should be very bright far-infrared sources and include some obscured AGN missed even by X-ray observations.

Using early data from the Infrared Array Camera (IRAC) on the Spitzer Space Telescope, taken for the Great Observatories Origins Deep Survey (GOODS), we identify and study objects that are well-detected at 3.6um, but are very faint (and in some cases, invisible) in the Hubble Ultra Deep Field (HUDF) ACS and NICMOS images and in very deep VLT Ks-band imaging. We select a sample of 17 objects with f_nu(3.6um)/f_nu(z_{850})>20. The analysis of their spectral energy-distributions (SEDs) from 0.4 to 8.0um shows that the majority of these objects cannot be satisfactorily explained without a well-evolved stellar population. We find that most of them can be well fitted by a simple two-component model, where the primary component represents a massive, old population that dominates the strong IR emission, while the secondary component represents a low amplitude, on-going star-formation process that accounts for the weak optical fluxes. Their estimated photometric redshifts (z_p) range from 1.6 to 2.9, with the median at z_p=2.4. For the simple star formation histories considered here, their corresponding stellar masses range from 0.1 to 1.6e11 M_sun for a Chabrier Initial Mass Function (IMF). Their median rest-frame Ks-band absolute magnitude is -22.9 mag in AB system, or 1.5xL*(K) for present--day elliptical galaxies. In the scenario of pure luminosity evolution, such objects may be direct progenitors for at least 14 to 51% of the local population of early--type galaxies. Due to the small cosmic volume of the HUDF, however, this simple estimate could be affected by other effects such as cosmic variance and the strong clustering of massive galaxies. A full analysis of the entire GOODS area is now underway to assess such effects.

We present results from the Hubble Higher-z Supernova Search, the first space-based open field survey for supernovae (SNe). In cooperation with the Great Observatories Origins Deep Survey, we have used the Hubble Space Telescope with the Advanced Camera for Surveys to cover 300 square arcmin in the area of the Chandra Deep Field South and the Hubble Deep Field North on five separate search epochs (separated by 45 day intervals) to a limiting magnitude of z'=26. These deep observations have allowed us to discover 42 SNe in the redshift range 0.2 < z < 1.6. As these data span a large range in redshift, they are ideal for testing the validity of Type Ia supernova progenitor models with the distribution of expected ``delay times,'' from progenitor star formation to SN Ia explosion, and the SN rates these models predict. Through a Bayesian maximum likelihood test, we determine which delay-time models best reproduce the redshift distribution of SNe Ia discovered in this survey. We find that models that require a large fraction of ``prompt'' (less than 2 Gyr) SNe Ia poorly reproduce the observed redshift distribution and are rejected at 95% confidence. We find that Gaussian models best fit the observed data for mean delay times in the range of 3 to 4 Gyr.

The Spitzer Space Telescope has undertaken the deepest ever observations of the 24 micron sky in the ELAIS-N1 field as part of GOODS Science Verification observations. We present the shape of the 24 micron source counts in the flux range 20-1000 microJy, discuss the redshift distribution and nature of these sources with particular emphasis on their near-infrared counterparts.

We studied the X-ray variability of sources detected in the Chandra Deep Field South (Giacconi et al. 2002), nearly all of which are low to moderate z AGN (Tozzi et al. 2001). We find that 45% of the sources with >100 counts exhibit significant variability on timescales ranging from a day up to a year. The fraction of sources found to be variable increases with observed flux, suggesting that >90% of all AGNs possess intrinsic variability. We also find that the fraction of variable sources appears to decrease with increasing intrinsic absorption; a lack of variability in hard, absorbed AGNs could be due to an increased contribution of reflected X-rays to the total flux. We do not detect significant spectral variability in the majority (~70%) of our sources. In half of the remaining 30%, the hardness ratio is anti-correlated with flux, mimicking the high/soft-low/hard states of galactic sources. The X-ray variability appears anti-correlated with the luminosity of the sources, in agreement with previous studies. High redshift sources, however, have larger variability amplitudes than expected from extrapolations of their low-z counterparts, suggesting a possible evolution in the accretion rate and/or size of the X-ray emitting region. Finally, we discuss some effects that may produce the observed decrease in the fraction of variable sources from z=0.5 out to z=2.

We examine the evolution of the sizes and number densities of disk galaxies using the high resolution images obtained by the Great Observatories Origins Deep Survey (GOODS) with the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope (HST). The multiwavelength (B,V,i,z) images allow us to classify galaxies based on their rest-frame B-band morphologies out to redshift z = 1.25. In order to minimize the effect of selection biases, we confine our analysis to galaxies which occupy the region of magnitude-size plane where the survey is almost 90% complete at all redshifts. The observed size distribution is consistent with a log-normal distribution as seen for the disk galaxies in the local Universe and does not show any significant evolution over the redshift range 0.25 < z < 1.25. We find that the number densities of disk galaxies remains fairly constant over this redshift range, although a modest evolution by a factor of four may be possible within the 2-sigma uncertainties.

Using U'- through Ks-band imaging data in the GOODS-South field, we construct a large, complete sample of 275 ``extremely red objects'' (EROs; K_s<22.0, R-K_s>3.35; AB), all with deep HST/ACS imaging in B_435, V_606, i_775, and z_850, and well-calibrated photometric redshifts. Quantitative concentration and asymmetry measurements fail to separate EROs into distinct morphological classes. We therefore visually classify the morphologies of all EROs into four broad types: ``Early'' (elliptical-like), ``Late'' (disk galaxies), ``Irregular'' and ``Other'' (chain galaxies and low surface brightness galaxies), and calculate their relative fractions and comoving space densities. For a broad range of limiting magnitudes and color thresholds, the relative number of early-type EROs is approximately constant at 33-44%, and the comoving space densities of Early- and Late-type EROs are comparable. Mean rest-frame spectral energy distributions (SEDs) at wavelengths between 0.1 and 1.2 um are constructed for all EROs. The SEDs are extremely similar in their range of shapes, independent of morphological type. The implication is that any differences between the broad-band SEDs of Early-type EROs and the other types are relatively subtle, and there is no robust way of photometrically distinguishing between different morphological types with usual optical/near-infrared photometry.

We use extensive multi-wavelength photometric data from the Great Observatories Origins Deep Survey (GOODS) to estimate photometric redshifts for a sample of 434 galaxies with spectroscopic redshifts in the Chandra Deep Field South. Using the Bayesian method, which incorporates redshift/magnitude priors, we estimate photometric redshifts for galaxies in the range 18 < R (AB) < 25.5, giving an rms scatter of 0.11. The outlier fraction is < 10%, with the outlier-clipped rms being 0.047. We examine the accuracy of photometric redshifts for several, special sub--classes of objects. The results for extremely red objects are more accurate than those for the sample as a whole, with rms of 0.051 and very few outliers (3%). Photometric redshifts for active galaxies, identified from their X-ray emission, have a dispersion of 0.104, with 10% outlier fraction, similar to that for normal galaxies. Employing a redshift/magnitude prior in this process seems to be crucial in improving the agreement between photometric and spectroscopic redshifts.

We study the rest-frame morphology and structural properties of optically selected starburst galaxies at redshift z < 1, using multi-waveband (BViz) high resolution images taken by the Advanced Camera for Surveys (ACS) as part of the Great Observatories Origins Deep Survey (GOODS). We classify galaxies into starburst, early and late types by comparing their observed spectral energy distributions (SEDs) with local templates. We find that early-type systems have significantly higher rest-frame B -band concentration indices and AGN fraction (> 25%) than late-type spirals and optically-selected starbursts. These results are consistent with the scenario in which early-epoch (z > 1) gas-rich dissipative processes (e.g., major mergers) have played an important role in developing large central concentrations in early-type E/Sa galaxies, leading to concurrent growth of central black holes and bulge formation in some of these early merger events. The starbursts have, on average, larger asymmetries than our control sample of normal galaxies, suggesting that a significant fraction of the starburst activity is tidally triggered.

We present quantitative morphological analyses of 37 HST/WFPC2 counterparts of X-ray sources in the 1 Ms Chandra Deep Field-South (CDFS). We investigate: 1) 1-D surface brightness profiles via isophotal ellipse fitting; 2) 2-D, PSF- convolved, bulge+disk+nucleus profile-fitting; 3) asymmetry and concentration indices compared with all ~3000 sources in our three WFPC2 fields; and 4) near- neighbor analyses comparing local environments of X-ray sources versus the field control sample. Significant nuclear point-source optical components appear in roughly half of the resolved HST/WFPC2 counterparts, showing a narrow range of F_X/F_{opt,nuc} consistent with the several HST-unresolved X-ray sources (putative type-1 AGN) in our fields. We infer roughly half of the HST/WFPC2 counterparts host unobscured AGN, which suggests no steep decline in the type-1/type-2 ratio out to the redshifts z~0.5-1 typical of our sources. The concentration indices of the CDFS counterparts are clearly larger on average than those of the field distribution, at 5-sigma, suggesting that the strong correlation between central black hole mass and host galaxy properties (including concentration index) observed in nearby galaxies is already evident by z~0.5-1. By contrast, the asymmetry index distribution of the 21 resolved CDFS sources at I<23 is indistinguishable from the I<23 field. Moreover, the frequency of I<23 near neighbors around the CDFS counterparts is not significantly different from the field sample. These results, combined with previous similar findings for local samples, suggest that recent merger/ interaction history is not a good indicator of AGN activity over a substantial range of look-back time.

Understanding the physical formation of the Hubble sequence remains one of the most important unsolved astrophysical problems. Searches for proto-disks and proto-ellipticals can now be effectively done using deep wide-field Hubble Space Telescope images taken with the new Advanced Camera for Surveys. Through an analysis of concentrations (C), asymmetries (A) and clumpiness values (S) (CAS) of galaxies found in the GOODS Field South, we are able to identify objects possibly forming onto the Hubble sequence. Using this approach, we detect a sizeable population of star forming luminous diffuse objects and star forming luminous asymmetric objects between redshifts 0.5 < z < 2. These galaxies have extremely low light concentrations, or high asymmetries, with absolute magnitudes M_B < -19. The luminous diffuse objects are found in abundance between z = 1 - 2, with fewer objects at z > 2 and z < 1. The luminous asymmetric objects are found at a similar abundance, with a peak at z~1. We argue that these galaxies are a subset of modern disks and ellipticals in formation. The co-moving volume density of the luminous diffuse objects between z = 1 - 2 is similar to the local density of bright disk galaxies, with values ~5*10^{5} Gpc^{-3}. The SEDs of these objects are mostly consistent with starbursts, or star-forming normal galaxies, with average uncorrected for extinction star formation rates of ~4 M_0 yr^{-1}. These galaxies also host 35-40% of the star formation activity at 1 < z < 2. We briefly discuss the implications of these objects for understanding the origin of the Hubble sequence.

We present the first results of the pure parallel near-UV observations with WFPC2 in the Chandra Deep Field South. The main goal is to search for the population of objects that contribute to the rise in the volume-averaged star formation rate at intermediate redshifts (z<1). A total of 34 objects, were identified in the WFPC2 images and their counterparts found in the ACS GOODS images. We use template fitting of the spectral energy distributions to classify the objects and found that 17 are starbursts, 9 are late-type, 2 are early-type galaxies, and 6 are stars; two of the late-types and one of the starbursts are X-ray sources. The colors of the starbursts were reproduced by the stellar evolutionary synthesis code Starburst99 which shows a mixed population of bursts with ages < 1 Gyr. Analysis of the light concentration, asymmetry and clumpiness, shows that this sample is a mixed bag, containing dwarf ellipticals, early and late-spirals and peculiar objects which resemble mergers in progress.

We have identified a sample of 10 highly reliable off-nuclear (> 2 kpc) X-ray sources at z=0.03-0.25 in late-type host galaxies within the two GOODS fields (i.e., the two Chandra Deep Fields). The combination of the superb spatial resolution and great depth of the HST ACS and Chandra ACIS coverage in the GOODS fields is critical to the identification of these sources. We extend the study of this enigmatic population up to higher redshifts and larger look-back times than has been possible before, finding that the fraction of optically luminous galaxies exhibiting luminous off-nuclear sources at z~0.1 is larger than at the current time (~36% versus ~8%). The X-ray luminosities of the GOODS off-nuclear X-ray sources are comparable to or slightly greater than the integrated point-source X-ray luminosities of nearby galaxies, providing a possible constraint on the X-ray Luminosity Function at higher star-formation rates. X-ray variability demonstrates that several of these sources are likely single accreting black holes; for the majority, however, the number of X-ray counts is too low for X-ray variability studies.

We present discoveries of SNe Ia at z > 1 and the photometric diagnostic used to discriminate them from other types of SNe detected during the GOODS Hubble Space Telescope Treasury program with the Advanced Camera for Surveys (ACS). Photometric redshift measurements of the hosts combined with deep f606w, f775w, and f850lp imaging discriminates hydrogen-rich SNe II from SNe I at z > 1 by exploiting the ultraviolet (UV) deficit in the energy distributions of SNe Ia. This sorting allows rapid follow-up of space-based discoveries. Subsequent spectroscopy of 11 GOODS SNe Ia obtained from the ground and with the grism on ACS confirmed the reliability of our photometric screening. We present the highest-redshift spectrum of any supernova published to date, SN Ia 2002fw at z = 1.3 observed near maximum brightness with the ACS grism. The lack of UV flux for true SNe Ia provides an effective tool for our ongoing efforts to build a sample of SNe Ia at 1.2 < z < 1.8 which will extend the useful range of the magnitude-redshift relation of SNe Ia.

We describe a possible new class of X-ray sources that have robust detections in ultra-deep Chandra data, yet have no detections at all in our deep multi-band GOODS Hubble Space Telescope (HST) ACS images, which represent the highest quality optical imaging obtained to date on these fields. These extreme X-ray / Optical ratio sources ("EXO"s) have values of Fx/Fopt at least an order of magnitude above those generally found for other AGN, even those that are harbored by reddened hosts. We thus infer two possible scenarios: (1) if these sources lie at redshifts z < 6, then their hosts need to be exceedingly underluminous, or more reddened, compared with other known sources; (2) if these sources lie above z ~ 6-7, such that even their Lyman-alpha emission is redshifted out of the bandpass of our ACS z(850) filter, then their optical and X-ray fluxes can be accounted for in terms of relatively normal L* hosts and moderate-luminosity AGN.

In the Chandra Deep Field South 1Msec exposure we have found, at redshift 3.700 +- 0.005, the most distant Type 2 AGN ever detected. It is the source with the hardest X-ray spectrum with redshift z>3. The optical spectrum has no detected continuum emission to a 3sigma detection limit of ~3 10^{-19} ergs/s/cm^2/AA and shows narrow lines of Ly_alpha, CIV, NV, HeII, OVI, [OIII], and CIII]. Their FWHM line widths have a range of ~700-2300 km/s with an average of approximately ~1500 km/s. The emitting gas is metal rich (Z ~2.5-3 Z_solar). In the X-ray spectrum of 130 counts in the 0.5-7 keV band there is evidence for intrinsic absorption with N_H > 10^{24} cm^{-2}. An iron K_alpha line with rest frame energy and equivalent width of ~6.4 keV and ~1 keV, respectively, in agreement with the obscuration scenario, is detected at a 2sigma level. If confirmed by our forthcoming XMM observations this would be the highest redshift detection of FeK_alpha. Depending on the assumed cosmology and the X-ray transfer model, the 2-10 keV rest frame luminosity corrected for absorption is ~10^{45 +- 0.5} ergs/s, which makes our source a classic example of the long sought Type 2 QSOs. From standard population synthesis models, these sources are expected to account for a relevant fraction of the black-hole-powered QSO distribution at high redshift.

In this Paper we present the source catalog obtained from a 942 ks exposure of the Chandra Deep Field South (CDFS), using ACIS-I on the Chandra X-ray Observatory. Catalog generation proceeded simultaneously using two different methods; a method of our own design using a modified version of the SExtractor algorithm, and a wavelet transform technique developed specifically for Chandra observations. The detection threshold has been set in order to have less than 10 spurious sources, as assessed by extensive simulations. We subdivided the catalog into four sections. The primary list consists of objects common to the two detection methods. Two secondary lists contain sources which were detected by: 1) the SExtractor algorithm alone and 2) the wavelet technique alone. The fourth list consists of possible diffuse or extended sources. The flux limits at the aimpoint for the soft (0.5--2 keV) and hard (2--10 keV) bands are 5.5E-17 erg/s/cm^2 and 4.5E-16 erg/s/cm^2 respectively. The total number of sources is 346; out of them, 307 were detected in the 0.5--2 keV band, and 251 in the 2--10 keV band. We also present optical identifications for the catalogued sources. Our primary optical data is R band imaging to a depth of R~26.5 (Vega). We found that the R-band/Chandra offsets are small, ~1 arcsec. Coordinate cross-correlation finds 85% of the Chandra sources covered in R to have counterparts within the 3-sigma error box (>~1.5 arcsec depending on off-axis angle and signal-to-noise). The unidentified fraction of sources, approximately ~10--15 %, is close to the limit expected from the observed X-ray flux to R-band ratio distribution for the identified sample.

We present the main results from our 940 ksec observation of the Chandra Deep Field South (CDFS), using the source catalog described in an accompanying paper (Giacconi et al. 2001). We extend the measurement of source number counts to 5.5e-17 erg/cm^2/s in the soft 0.5-2 keV band and 4.5e-16 erg/cm^2/s in the hard 2-10 keV band. The hard band LogN-LogS shows a significant flattening (slope~=0.6) below ~1e-14 erg/cm^2/s, leaving at most 10-15% of the X-ray background (XRB) to be resolved, the main uncertainty lying in the measurement of the total flux of the XRB. On the other hand, the analysis in the very hard 5-10 keV band reveals a relatively steep LogN-LogS (slope ~=1.3) down to 1e-15 erg/cm^2/s. Together with the evidence of a progressive flattening of the average X-ray spectrum near the flux limit, this indicates that there is still a non negligible population of faint hard sources to be discovered at energies not well probed by Chandra, which possibly contribute to the 30 keV bump in the spectrum of the XRB. We use optical redshifts and identifications, obtained with the VLT, for one quarter of the sample to characterize the combined optical and X-ray properties of the CDFS sample. Different source types are well separated in a parameter space which includes X-ray luminosity, hardness ratio and R-K color. Type II objects, while redder on average than the field population, have colors which are consistent with being hosted by a range of galaxy types. Type II AGN are mostly found at z<~1, in contrast with predictions based on AGN population synthesis models, thus suggesting a revision of their evolutionary parameters.

We present HST/WFPC2 observations of a well-defined sample of 40 X-ray sources with X-ray fluxes above the detection threshold of the full 1 Msec Chandra Deep Field South (CDFS). The sensitivity and spatial resolution of our HST observations are sufficient to detect the optical counterparts of 37 of the X-ray sources, yielding information on their morphologies and environments. In this paper we extend the results obtained in our previous study on the 300 ks CDFS X-ray data (Schreier et al. 2001, Paper I). Specifically, we show that the optical counterparts to the X-ray sources are divided into two distinct populations: 1) an optically faint group with relatively blue colors, similar to the faint blue field galaxy population, and 2) an optically brighter group, including resolved galaxies with average colors significantly redder than the corresponding bright field galaxy population. The brighter objects comprise a wide range of types, including early and late type galaxies, starbursts, and AGN. By contrast, we show that the faint blue X-ray population are most consistent with being predominantly Type 2 AGN of low to moderate luminosity, located at higher redshifts (z ~ 1 - 2). This conclusion is supported by luminosity function models of the various classes of objects. Hence, the combination of deep X-ray data with the high spatial resolution of HST are for the first time allowing us to probe the faint end of the AGN luminosity function at cosmologically interesting redshifts.

We present results from 300 ks of X-ray observations of the Chandra Deep Field South. The field of the four combined exposures is now 0.1035 deg^2 and we reach a flux limit of 10^{-16} erg s^{-1} cm^{-2} in the 0.5-2 keV soft band and 10^{-15} erg s^{-1} cm^{-2} in the 2-10 keV hard band, thus a factor 2 fainter than the previous 120 ks exposure. The total catalogue is composed of 197 sources including 22 sources detected only in the hard band, 51 only in the soft band, and 124 detected in both bands. We have now the optical spectra for 86 optical counterparts. We compute the total contribution to the X-ray background in the 2-10 keV band, which now amounts to (1.45\pm 0.15)*10^{-11} erg cm^{-2} s^{-1} deg^{-2} (after the inclusion of the ASCA sources to account for the bright end) to a flux limit of 10^{-15} erg s^{-1} cm^{-2}. This corresponds to 60-90% of the unresolved hard X-ray background (XRB), given the uncertainties on its actual value.[ABRIDGED]

We present preliminary results from imaging three HST/WFPC2 fields in V and I within the Chandra Deep Field South (CDFS). HST's sensitivity and resolution are sufficient to reveal optical counterparts for 24 of the 26 CDFS X-ray sources detected in the 300 ksec X-ray catalog and to determine the morphologies of most of these. We find that the X-ray sources comprise two apparently distinct populations of optical candidates: one optically faint (I > 24) with V - I colors consistent with the I > 24 field population; the other significantly brighter (I < 22) with colors redder than the I < 22 field population. More than 2/3 of the X-ray source counterparts are resolved galaxies. The brighter sources are mostly AGN, based on their high X-ray luminosity. The optically resolved sources in the brighter population have a very narrow range of V - I color and appear to be a mix of both late and early type morphologies at low to moderate redshift. We show that the second population, with fainter optical counterparts, can be explained as higher redshift Type 2 AGN.

We assess the morphology-density relation, H-alpha linewidth-density relation, and the effects of nearby companions for 300 galaxies selected from low-density environments of the Center for Astrophysics Redshift Survey. Both the morphological mix and the EW(H-alpha) distribution of galaxies at modest underdensites are indistinguishable from our control sample at modest overdensities (early-type and absorption-line fractions both ~35%). At the lowest densities: the absorption-line fraction drops to 15%; more galaxies have active star formation; and there are more irregular and peculiar morphologies. Our magnitude-limited redshift survey of 250 fainter projected companions reveals that the void galaxies are no less likely to have a neighbor in redshift space. However, the typical velocity separation of close pairs drops sharply at densities below half the mean, from >200 km/s to ~100 km/s. The EW(H-alpha) distribution for galaxies lacking nearby companions varies little with density, while galaxies with nearby companions are much more likely to have strong current star formation (H-alpha linewidths >20 \AA) in the voids. These results and the luminosity- and color-density relations of this sample (astro-ph/9910073) suggest that the formation and evolution of field galaxies are only sensitive to large-scale underdensity below a threshold of roughly half the mean density. The differences in galaxy properties at the lower densities may be explained by: 1) a relative scarcity of the small-scale primordial density enhancements needed to form massive early-type/absorption- line galaxies; 2) more effective galaxy encounters due to the lower velocity dispersion on small scales; and 3) more gas and dust available for triggered star formation, presuming that luminous void galaxies formed more recently.

We study the optical properties of a large sample of galaxies in low-density regions of the nearby universe. We make a 5/h Mpc-smoothed map of the galaxy density throughout the Center for Astrophysics Redshift Survey (CfA2) to identify galaxies within three prominent nearby ``voids'' with diameter >30/h Mpc. We augment the CfA2 void galaxy sample with fainter galaxies found in the same regions from the deeper Century and 15R Redshift Surveys. We obtain B and R CCD images and high-S/N longslit spectra for the resulting sample of 149 void galaxies, and also for 131 galaxies on the periphery of these voids. Here we describe the photometry for the sample, including B isophotal magnitudes and B-R colors. For the 149 galaxies which lie in regions below the mean survey density, the luminosity functions in B and R are well-fit by Schechter functions with respective parameters (alpha_B = -0.5 +/- 0.3, B_* = -18.9 +/- 0.2) and (alpha_R = -0.9 +/- 0.3, R_* = -20.4 +/- 0.3). The B luminosity function (LF) is consistent with typical survey LFs (e.g., the Southern Sky Redshift Survey), and the R LF is consistent with the Century Survey. The B and R LFs of 131 galaxies in the ``void periphery'', regions between the mean density and twice the mean, have similar Schechter parameters. The CfA2 LF is inconsistent with both samples at 3.5-sigma. When we narrow our analysis to the 46 galaxies in regions below half the mean density, the LF is steeper: alpha = -1.4 +/- 0.5. The typical survey LFs are inconsistent with this subsample at 2-sigma. The B-R color distribution of galaxies in the lowest-density regions is also shifted significantly (3-sigma) blue-ward of the higher density samples. The most luminous red galaxies (R < -21) are absent from the lowest density regions at 2.5-sigma.

We examine models of the mass distribution for the first known case of gravitational lensing. Several new sets of constraints are used, based on recent observations. We remodel the VLBI observations of the radio jets in the two images of Q0957+561, showing that the previously derived positions and uncertainties were incorrect. We use as additional constraints the candidate lensed pairs of galaxies identified recently with the Hubble Space Telescope. We explore a wider range of lens models than before, and find that the Hubble constant is not tightly constrained once elliptical lens models are considered. We also discuss the systematic uncertainties caused by the cluster containing the lens galaxy. We conclude that additional observations of the candidate lensed galaxies as well as direct measurements of the cluster mass profile are needed before an accurate value of the Hubble constant can be derived from this lens.

We compile 1113 redshifts (648 new measurements, 465 from the literature) for Zwicky catalogue galaxies in the region (-3.5deg < \delta < 8.5deg, 11.5h < \alpha < 13.5h). We include redshifts for 114 component objects in 78 Zwicky catalogue multiplets. The redshift survey in this region is 99.5% complete to the Zwicky catalogue limit, m_Zw=15.7. It is 99.9% complete to m_Zw = 15.5, the CfA Redshift Survey (CfA2) magnitude limit. The survey region is adjacent to the northern portion of CfA2, overlaps the northernmost slice of the Las Campanas Redshift Survey, includes the southern extent of the Virgo Cluster, and is roughly centered on the QSO 3C 273. As in other portions of the Zwicky catalogue, bright and faint galaxies trace the same large-scale structure.

We study the galaxy number density (smoothed on a 5h^{-1} Mpc scale) around 18 low-redshift Lyman-alpha absorbers previously observed with HST. The absorbers lie in the foregrounds of Mrk 335, Mrk 421, Mrk 501, I Zw 1, and 3C 273, all within regions where there are now complete redshift surveys to m_{Zw}=15.5. We construct a smoothed galaxy number density field from the redshift survey data and determine the distribution of densities at the Lyman-alpha absorber locations. We also find the distribution of galaxy number density for a variety of test samples: all galaxy locations within the Center for Astrophysics Redshift Survey (CfA2), CfA2 galaxy locations along randomly selected lines of sight (LOS), and randomly chosen redshifts along random LOS. The Lyman-alpha absorbers are present in dense regions of the survey, but occur far more frequently in underdense regions than do typical luminous galaxies. The distribution of smoothed galaxy density around the Lyman-alpha absorbers is inconsistent at the 4-sigma level with the density distribution around survey galaxies. It is highly consistent with a density distribution at randomly chosen redshifts along random LOS. This supports earlier evidence that the nearby, low column density (log N_{HI} < 14) Lyman-alpha forest systems are spatially distributed at random; they are not well correlated with the local large-scale structure.

We present a simple mass model for the lensing galaxy in the gravitationally lensed quasar 0957+561. The model is a generalization of the singular isothermal sphere and includes a core radius, $r_c$, and a power-law index, $\eta$, defined such that mass increases as $r^\eta$ at large radius. We approximate the galaxy cluster surrounding the lensing galaxy with a quadratic potential described by its convergence $\kappa$ and shear $\gamma$. We fit the model to a recent high resolution VLBI map of the two images of 0957+561. We obtain a tight constraint on the radial index, $1.07 < \eta <1.18$, which means that the lens galaxy is nearly isothermal with increasing mass-to-light ratio out to $15 h^{-1}$ kpc. We also obtain an upper limit on the core radius, $r_c < 330 h^{-1}$ pc. We use the model to calculate the Hubble constant $H_0$ as a function of the time delay $\Delta\tau_{BA}$ between the two images: $H_0 = \left({ 82.5^{+5.9} _{-3.0} }\right) (1 - \kappa) \left({ \Delta\tau_{BA}/1.1 \,{\rm yr} }\right)^{-1}$ km/s/Mpc. Once $\Delta \tau_{BA}$ is measured, this will provide an upper bound on $H_0$ since $\kappa$ cannot be negative. In addition, the model degeneracy due to $\kappa$ can be eliminated if the one-dimensional velocity dispersion $\sigma$ of the lensing galaxy is measured. We then have $H_0 = \left({ 82.5^{+8.7} _{-5.6} }\right) (\sigma / 322\,{\rm km/s})^2 \left({ \Delta \tau_{BA} /1.1 \,{\rm yr} }\right)^{-1}$ km/s/Mpc. We investigate the effects of ellipticity in the lensing galaxy and clumpiness in the lensing cluster and find that these cause little change in our results.