We present the Spitzer IRAC/MUSYC Public Legacy Survey in the Extended CDF-South (SIMPLE), which consists of deep IRAC observations covering the ~1,600 arcmin^2 area surrounding GOODS-S. The limiting magnitudes of the SIMPLE IRAC mosaics typically are 23.8, 23.6, 21.9, and 21.7, at 3.6 um, 4.5 um, 5.8 um, and 8.0 um, respectively (5-sigma total point source magnitudes in AB). The SIMPLE IRAC images are combined with the 10'x15' GOODS IRAC mosaics in the center. We give detailed descriptions of the observations, data reduction, and properties of the final images, as well as the detection and photometry methods used to build a catalog. Using published optical and near-infrared data from the Multiwavelength Survey by Yale-Chile (MUSYC), we construct an IRAC-selected catalog, containing photometry in UBVRIz'JHK, [3.6 um], [4.5 um], [5.8 um], and [8.0 um]. The catalog contains 43,782 sources with S/N > 5 at 3.6 um, 19,993 of which have 13-band photometry. We compare this catalog to the publicly available MUSYC and FIREWORKS catalogs and discuss the differences. Using a high signal-to-noise sub-sample of 3,391 sources with ([3.6] + [4.5])/2 < 21.2, we investigate the star formation rate history of massive galaxies out to z ~ 1.8. We find that at z ~ 1.8 at least 30% +/-7% of the most massive galaxies (Mstar > 10^11 Msol) are passively evolving, in agreement with earlier results from surveys covering less area.

We present deep optical 18-medium-band photometry from the Subaru telescope over the ~30' x 30' Extended Chandra Deep Field-South (ECDF-S), as part of the Multiwavelength Survey by Yale-Chile (MUSYC). This field has a wealth of ground- and space-based ancillary data, and contains the GOODS-South field and the Hubble Ultra Deep Field. We combine the Subaru imaging with existing UBVRIzJHK and Spitzer IRAC images to create a uniform catalog. Detecting sources in the MUSYC BVR image we find ~40,000 galaxies with R_AB<25.3, the median 5 sigma limit of the 18 medium bands. Photometric redshifts are determined using the EAZY code and compared to ~2000 spectroscopic redshifts in this field. The medium band filters provide very accurate redshifts for the (bright) subset of galaxies with spectroscopic redshifts, particularly at 0.1 < z < 1.2 and at z > 3.5. For 0.1 < z < 1.2, we find a 1 sigma scatter in \Delta z/(1+z) of 0.007, similar to results obtained with a similar filter set in the COSMOS field. As a demonstration of the data quality, we show that the red sequence and blue cloud can be cleanly identified in rest-frame color-magnitude diagrams at 0.1 < z < 1.2. We find that ~20% of the red-sequence-galaxies show evidence of dust-emission at longer rest-frame wavelengths. The reduced images, photometric catalog, and photometric redshifts are provided through the public MUSYC website.

We study the properties of a sample of 211 heavily-obscured Active Galactic Nucleus (AGN) candidates in the Extended Chandra Deep Field-South selecting objects with f_24/f_R>1000 and R-K>4.5. Of these, 18 were detected in X-rays and found to be obscured AGN with neutral hydrogen column densities of ~10^23 cm^-2. In the X-ray undetected sample, the following evidence suggests a large fraction of heavily-obscured (Compton Thick) AGN: (i) The stacked X-ray signal of the sample is strong, with an observed ratio of soft to hard X-ray counts consistent with a population of ~90% heavily obscured AGN combined with 10% star-forming galaxies. (ii) The X-ray to mid-IR ratios for these sources are significantly larger than that of star-forming galaxies and ~2 orders of magnitude smaller than for the general AGN population, suggesting column densities of N_H>5x10^24 cm^-2. (iii) The Spitzer near- and mid-IR colors of these sources are consistent with those of the X-ray-detected sample if the effects of dust self-absorption are considered. Spectral fitting to the rest-frame UV/optical light (dominated by the host galaxy) returns stellar masses of ~10^11 M_sun and <E(B-V)> =0.5, and reveals evidence for a significant young stellar population, indicating that these sources are experiencing considerable star-formation. This sample of heavily-obscured AGN candidates implies a space density at z~2 of ~10^-5 Mpc^-3, finding a strong evolution in the number of L_X>10^44 erg/s sources from z=1.5 to 2.5, possibly consistent with a short-lived heavily-obscured phase before an unobscured quasar is visible.

We investigate the star formation history of the universe using FIREWORKS, a multiwavelength survey of the CDFS. We study the evolution of the specific star formation rate (sSFR) with redshift in different mass bins from z = 0 to z ~ 3. We find that the sSFR increases with redshift for all masses. The logarithmic increase of the sSFR with redshift is nearly independent of mass, but this cannot yet be verified at the lowest-mass bins at z > 0.8, due to incompleteness. We convert the sSFRs to a dimensionless growth rate to facilitate a comparison with a semi-analytic galaxy formation model that was implemented on the Millennium Simulation. The model predicts that the growth rates and sSFRs increase similarly with redshift for all masses, consistent with the observations. However, we find that for all masses, the inferred observed growth rates increase more rapidly with redshift than the model predictions. We discuss several possible causes for this discrepancy, ranging from field-to-field variance, conversions to SFR, and shape of the IMF. We find that none of these can solve the discrepancy completely. We conclude that the models need to be adapted to produce the steep increase in growth rate between redshift z=0 and z=1.

We present a new K-selected, optical-to-near-infrared photometric catalog of the Extended Chandra Deep Field South (ECDFS), making it publicly available to the astronomical community. The dataset is founded on publicly available imaging, supplemented by original zJK imaging data obtained as part of the MUltiwavelength Survey by Yale-Chile (MUSYC). The final photometric catalog consists of photometry derived from nine band U-K imaging covering the full 0.5x0.5 sq. deg. of the ECDFS, plus H band data for approximately 80% of the field. The 5sigma flux limit for point-sources is K = 22.0 (AB). This is also the nominal completeness and reliability limit of the catalog: the empirical completeness for 21.75 < K < 22.00 is 85+%. We have verified the quality of the catalog through both internal consistency checks, and comparisons to other existing and publicly available catalogs. As well as the photometric catalog, we also present catalogs of photometric redshifts and restframe photometry derived from the ten band photometry. We have collected robust spectroscopic redshift determinations from published sources for 1966 galaxies in the catalog. Based on these sources, we have achieved a (1sigma) photometric redshift accuracy of Dz/(1+z) = 0.036, with an outlier fraction of 7.8%. Most of these outliers are X-ray sources. Finally, we describe and release a utility for interpolating restframe photometry from observed SEDs, dubbed InterRest. Particularly in concert with the wealth of already publicly available data in the ECDFS, this new MUSYC catalog provides an excellent resource for studying the changing properties of the massive galaxy population at z < 2. (Abridged)

We present the first results of our optical spectroscopy program aimed to provide redshifts and identifications for the X-ray sources in the Extended Chandra Deep Field South. A total of 339 sources were targeted using the IMACS spectrograph at the Magellan telescopes and the VIMOS spectrograph at the VLT. We measured redshifts for 186 X-ray sources, including archival data and a literature search. We find that the AGN host galaxies have on average redder rest-frame optical colors than non-active galaxies, and that they live mostly in the "green valley". The dependence of the fraction of AGN that are obscured on both luminosity and redshift is confirmed at high significance and the observed AGN space density is compared with the expectations from existing luminosity functions. These AGN show a significant difference in the mid-IR to X-ray flux ratio for obscured and unobscured AGN, which can be explained by the effects of dust self-absorption on the former. This difference is larger for lower luminosity sources, which is consistent with the dust opening angle depending on AGN luminosity.

We study the evolution of the star formation rate (SFR) of mid-infrared (IR) selected galaxies in the extended Chandra Deep Field South (E-CDFS). We use a combination of U-K GaBoDS and MUSYC data, deep IRAC observations from SIMPLE, and deep MIPS data from FIDEL. This unique multi-wavelength data set allows us to investigate the SFR history of massive galaxies out to redshift z ~ 1.8. We determine star formation rates using both the rest-frame ultraviolet luminosity from young, hot stars and the total IR luminosity of obscured star formation obtained from the MIPS 24 um flux. We find that at all redshifts the galaxies with higher masses have substantially lower specific star formation rates than lower mass galaxies. The average specific star formation rates increase with redshift, and the rate of incline is similar for all galaxies (roughly (1+z)^{n}, n = 5.0 +/- 0.4). It does not seem to be a strong function of galaxy mass. Using a subsample of galaxies with masses M_*> 10^11 M_sun, we measured the fraction of galaxies whose star formation is quenched. We consider a galaxy to be in quiescent mode when its specific star formation rate does not exceed 1/(3 x t_H), where t_H is the Hubble time. The fraction of quiescent galaxies defined as such decreases with redshift out to z ~ 1.8. We find that, at that redshift, 19 +/-9 % of the M_* > 10^11 M_sun sources would be considered quiescent according to our criterion.

We present the mid-infrared colors of X-ray-detected AGN and explore mid-infrared selection criteria. Using a statistical matching technique, the likelihood ratio, over 900 IRAC counterparts were identified with a new MUSYC X-ray source catalog that includes ~1000 published X-ray sources in the Chandra Deep Field-South and Extended Chandra Deep Field-South. Most X-ray-selected AGN have IRAC spectral shapes consistent with power-law slopes, f_{nu} ~ nu^{alpha}, and display a wide range of colors, -2 < alpha < 2. Although X-ray sources typically fit to redder (more negative alpha) power-laws than non-X-ray detected galaxies, more than 50% do have flat or blue (galaxy-like) spectral shapes in the observed 3-8 micron band. Only a quarter of the X-ray selected AGN detected at 24 micron are well fit by featureless red power laws in the observed 3.6-24 micron, likely the subset of our sample whose infrared spectra are dominated by emission from the central AGN region. Most IRAC color-selection criteria fail to identify the majority of X-ray-selected AGN, finding only the more luminous AGN, the majority of which have broad emission lines. In deep surveys, these color-selection criteria select 10-20% of the entire galaxy population and miss many moderate luminosity AGN.

We investigate the stellar populations of a sample of 162 Lyman-alpha emitting galaxies (LAEs) at z = 3.1 in the Extended Chandra Deep Field South, using deep Spitzer IRAC data available from the GOODS and SIMPLE surveys to derive reliable stellar population estimates. We divide the LAEs according to their rest-frame near-IR luminosities into IRAC-detected and IRAC-undetected samples. About 70% of the LAEs are undetected in 3.6 micron down to [3.6] = 25.2 AB. Stacking analysis reveals that the average stellar population of the IRAC-undetected sample has an age of ~ 200 Myr and a mass of ~ 3x10^8 solar masses, consistent with the expectation that LAEs are mostly young and low-mass galaxies. On the other hand, the IRAC-detected LAEs are on average significantly older and more massive, with an average age > 1 Gyr and mass ~ 10^10 solar masses. Comparing the IRAC colors and magnitudes of the LAEs to z ~ 3 Lyman break galaxies (LBGs) shows that the IRAC-detected LAEs lie at the faint blue end of the LBG color-magnitude distribution, suggesting that IRAC-detected LAEs may be the low mass extension of the LBG population. We also present tentative evidence for a small fraction (~ 5%) of obscured AGN within the LAE sample. Our results suggest that LAEs posses a wide range of ages and masses. Additionally, the presence of evolved stellar populations inside LAEs suggests that the Lyman-alpha luminous phase of galaxies may either be a long-lasting or recurring phenomenon.

We investigate with unprecedented accuracy the correlations between the dynamical mass-to-light ratio M/L and other global observables of E and S0 galaxies. We construct two-integral Jeans and three-integral Schwarzschild dynamical models for a sample of 25 E/S0 galaxies with SAURON integral-field stellar kinematics. We find a tight correlation of the form (M/L)=(3.80+/-0.14)*(sigma/200 km/s)^(0.84+/-0.07) between the dynamical M/L (in the I-band) and the luminosity-weighted second moment (sigma) of the line-of-sight velocity-distribution within Re. The observed rms scatter in M/L for our sample is 18%, while the inferred intrinsic scatter is ~13%. The (M/L)-sigma relation can be included in the remarkable series of tight correlations between sigma and other galaxy global observables. The comparison of the observed correlations with the predictions of the Fundamental Plane (FP), and with simple virial estimates, shows that the `tilt' of the FP of early-type galaxies, is due to a real M/L variation, while structural and orbital non-homology have a negligible effect. The virial mass is a reliable estimator of the mass in the central regions of galaxies. The best-fitting virial relation has the form (M/L)_vir=(5.0+/-0.1)*Re*sigma^2/(L*G). The comparison of the dynamical M/L with the (M/L)_pop inferred from the analysis of the stellar population, indicates a median dark matter fraction in early-type galaxies of ~30% of the total mass inside one Re. (Abridged)

We investigate the well-known correlations between the dynamical mass-to-light ratio M/L and other global observables of elliptical (E) and lenticular (S0) galaxies. We construct two-integral Jeans and three-integral Schwarzschild dynamical models for a sample of 25 E/S0 galaxies with SAURON integral-field stellar kinematics to about one effective (half-light) radius R_e. The comparison of the dynamical M/L with the (M/L)_pop inferred from the analysis of the stellar population, indicates that dark matter in early-type galaxies contributes ~30% of the total mass inside one R_e, in agreement with previous studies, with significant variations from galaxy to galaxy. Our results suggest a variation in M/L at constant (M/L)_pop, which seems to be linked to the galaxy dynamics. We speculate that fast rotating galaxies have lower dark matter fractions than the slow rotating and generally more massive ones.