Source author record

Bahram Mobasher

Bahram Mobasher appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

114works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

114 published item(s)

preprint2025arXiv

COSMOS2025: The COSMOS-Web galaxy catalog of photometry, morphology, redshifts, and physical parameters from JWST, HST, and ground-based imaging

We present COSMOS2025, the COSMOS-Web catalog of photometry, morphology, photometric redshifts and physical parameters for more than 700,000 galaxies in the Cosmic Evolution Survey (COSMOS) field. This catalog is based on our \textit{James Webb Space Telescope} 255\,h COSMOS-Web program, which provides deep near-infrared imaging in four NIRCam (F115W, F150W, F277W, F444W) and one MIRI (F770W) filter over the central $\sim 0.54 {\, \rm deg}^2$ ($\sim 0.2 {\, \rm deg}^2$ for MIRI) in COSMOS. These data are combined with ground- and space-based data to derive photometric measurements of NIRCam-detected sources using both fixed-aperture photometry (on the space-based bands) and a profile-fitting technique on all 37 bands spanning 0.3-8 micron. We provide morphology for all sources from complementary techniques including profile fitting and machine-learning classification. We derive photometric redshifts, physical parameters and non-parametric star formation histories from spectral energy distribution (SED) fitting. The catalog has been extensively validated against previous COSMOS catalogs and other surveys. Photometric redshift accuracy measured using spectroscopically confirmed galaxies out to $z\sim9$ reaches $σ_{\rm MAD} = 0.012$ at $m_{\rm F444W}<28$ and remains at $σ_{\rm MAD} \lesssim 0.03$ as a function of magnitude, color, and galaxy type. This represents a factor of $\sim 2$ improvement at 26 AB mag compared to COSMOS2020. The catalog is approximately 80\% complete at $\log(M_{\star}/{\rm M}_{\odot}) \sim 9$ at $z \sim 10$ and at $\log(M_{\star}/{\rm M}_{\odot}) \sim 7$ at $z \sim 0.2$, representing a gain of 1\,dex compared to COSMOS2020. COSMOS2025 represents the definitive COSMOS-Web catalog. It is provided with complete documentation, together with redshift probability distributions, and it is ready for scientific exploitation today.

preprint2022arXiv

CO Emission, Molecular Gas, and Metallicity in Main-Sequence Star-Forming Galaxies at $z\sim2.3$

We present observations of CO(3-2) in 13 main-sequence $z=2.0-2.5$ star-forming galaxies at $\log(M_*/M_{\odot})=10.2-10.6$ that span a wide range in metallicity (O/H) based on rest-optical spectroscopy. We find that CO(3-2)/SFR decreases with decreasing metallicity, implying that the CO luminosity per unit gas mass is lower in low-metallicity galaxies at $z\sim2$. We constrain the CO-to-H$_2$ conversion factor ($α_{\text{CO}}$) and find that $α_{\text{CO}}$ inversely correlates with metallicity at $z\sim2$. We derive molecular gas masses ($M_{\text{mol}}$) and characterize the relations among $M_*$, SFR, $M_{\text{mol}}$, and metallicity. At $z\sim2$, $M_{\text{mol}}$ increases and molecular gas fraction ($M_{\text{mol}}$/$M_*$) decrease with increasing $M_*$, with a significant secondary dependence on SFR. Galaxies at $z\sim2$ lie on a near-linear molecular KS law that is well-described by a constant depletion time of 700 Myr. We find that the scatter about the mean SFR-$M_*$, O/H-$M_*$, and $M_{\text{mol}}$-$M_*$ relations is correlated such that, at fixed $M_*$, $z\sim2$ galaxies with larger $M_{\text{mol}}$ have higher SFR and lower O/H. We thus confirm the existence of a fundamental metallicity relation at $z\sim2$ where O/H is inversely correlated with both SFR and $M_{\text{mol}}$ at fixed $M_*$. These results suggest that the scatter of the $z\sim2$ star-forming main sequence, mass-metallicity relation, and $M_{\text{mol}}$-$M_*$ relation are primarily driven by stochastic variations in gas inflow rates. We place constraints on the mass loading of galactic outflows and perform a metal budget analysis, finding that massive $z\sim2$ star-forming galaxies retain only 30% of metals produced, implying that a large mass of metals resides in the circumgalactic medium.

preprint2022arXiv

H$α$ Morphologies of Star Clusters in 16 LEGUS Galaxies: Constraints on HII region evolution timescales

The analysis of star cluster ages in tandem with the morphology of their HII regions can provide insight into the processes that clear a cluster's natal gas, as well as the accuracy of cluster ages and dust reddening derived from Spectral Energy Distribution (SED) fitting. We classify 3757 star clusters in 16 nearby galaxies according to their H$α$ morphology (concentrated, partially exposed, no emission), using Hubble Space Telescope (HST) imaging from the Legacy ExtraGalactic Ultraviolet Survey (LEGUS). We find: 1) The mean SED ages of clusters with concentrated (1-2 Myr) and partially exposed HII region morphologies (2-3 Myr) indicate a relatively early onset of gas clearing and a short (1-2 Myr) clearing timescale. 2) The reddening of clusters can be overestimated due to the presence of red supergiants, which is a result of stochastic sampling of the IMF in low mass clusters. 3) The age-reddening degeneracy impacts the results of the SED fitting - out of 1408 clusters with M$_*$ $\geq$ 5000 M$_{\odot}$, we find that at least 46 (3%) have SED ages which appear significantly underestimated or overestimated based on H$α$ and their environment, while the total percentage of poor age estimates is expected to be several times larger. 4) Lastly, we examine the dependence of the morphological classifications on spatial resolution. At HST resolution, our conclusions are robust to the distance range spanned by the sample (3-10 Mpc). However, analysis of ground-based H$α$ images shows that compact and partially exposed morphologies frequently cannot be distinguished from each other.

preprint2022arXiv

Implications of a Temperature Dependent IMF II: An Updated View of the Star-Forming Main Sequence

The stellar initial mass function (IMF) is predicted to depend upon the temperature of gas in star-forming molecular clouds. The introduction of an additional parameter, $T_{IMF}$ , into photometric template fitting, allows galaxies to be fit with a range of IMFs. Three surprising new features appear: (1) most star-forming galaxies are best fit with a bottom-lighter IMF than the Milky Way; (2) most star-forming galaxies at fixed redshift are fit with a very similar IMF; and (3) the most massive star-forming galaxies at fixed redshift instead exhibit a less bottom-light IMF, similar to that measured in quiescent galaxies. Additionally, since stellar masses and star formation rates both depend on the IMF, these results slightly modify the resulting relationship, while yielding similar qualitative characteristics to previous studies.

preprint2022arXiv

Implications of a Temperature Dependent IMF III: Mass Growth and Quiescence

The stellar initial mass function (IMF) is predicted to depend upon the temperature of gas in star-forming molecular clouds. The introduction of an additional parameter, $T_{IMF}$ , into photometric template fitting, suggest most galaxies obey an IMF top-heavier than the Galactic IMF. The implications of these revised fits on mass functions, quiescence and turnoff are discussed. At all redshifts the highest mass galaxies become quiescent first with the turnoff mass decreasing towards the present. The synchronous turnoff mass across galaxies suggests quiescence is driven by universal mechanisms rather than by stochastic or environmental processes.

preprint2022arXiv

IR SED and Dust Masses of Sub-solar Metallicity Galaxies at z~2.3

We present results from ALMA 1.2mm continuum observations of a sample of 27 star-forming galaxies at z=2.1-2.5 from the MOSFIRE Deep Evolution Field (MOSDEF) survey. These galaxies have gas-phase metallicity and star-formation rate measurements from Hb, [OIII], Ha, and [NII]. Using stacks of Spitzer, Herschel, and ALMA photometry (rest-frame ~ 8-400$μ$m), we examine the IR SED of high-redshift subsolar metallicity (~0.5 $Z_{\odot}$) LIRGs. We find that the data agree well with an average SED template of higher luminosity local low-metallicity dwarf galaxies (reduced $χ^2$ of 1.8). When compared with the commonly used templates for solar-metallicity local galaxies or high-redshift LIRGs and ULIRGs, even in the most favorable case (with reduced $χ^2$ of 2.8), the templates are rejected at >98% confidence level. The broader and hotter IR SED of both the local dwarfs and high-redshift subsolar metallicity galaxies may result from different grain properties, a clumpy dust geometry, or a harder/more intense ionizing radiation field that heats the dust to higher temperatures. The obscured SFR indicated by the FIR emission of the subsolar metallicity galaxies is only ~ 60% of the total SFR, which is considerably lower than that of the local LIRGs with ~ 96-97% obscured fractions. Due to the evolving IR SED shape, the local LIRG templates fit to mid-IR data can overestimate the Rayleigh-Jeans tail measurements at z~2 by a factor of 2-20, and these templates underestimate IR luminosities if fit to the observed ALMA fluxes by >0.4dex. At a given stellar mass or metallicity, dust masses at z~2.3 are an order of magnitude higher than those at z~0. Given the predicted molecular gas mass fractions, the observed z~2.3 dust-to-stellar mass ratios suggest lower dust-to-molecular gas masses than in local galaxies at the same metallicity.

preprint2022arXiv

Kinematics and Metallicity of Red Giant Branch Stars in the Northeast Shelf of M31

We obtained Keck/DEIMOS spectra of 556 individual red giant branch stars in 4 spectroscopic fields spanning $13-31$ projected kpc along the Northeast (NE) shelf of M31. We present the first detection of a complete wedge pattern in the space of projected M31-centric radial distance versus line-of-sight velocity for this feature, which includes the returning stream component of the shelf. This wedge pattern agrees with expectations of a tidal shell formed in a radial merger and provides strong evidence in favor of predictions of Giant Stellar Stream (GSS) formation models in which the NE shelf originates from the second orbital wrap of the tidal debris. The observed concentric wedge patterns of the NE, West (W), and Southeast (SE) shelves corroborate this interpretation independently of the models. We do not detect a kinematical signature in the NE shelf region corresponding to an intact progenitor core, favoring GSS formation models in which the progenitor is completely disrupted. The shelf's photometric metallicity distribution implies that it is dominated by tidal material, as opposed to the phase-mixed stellar halo or the disk. The metallicity distribution ([Fe/H]$_{\rm phot}$ = $-0.42$ $\pm$ $0.01$) also matches the GSS, and consequently the W and SE shelves, further supporting a direct physical association between the tidal features.

preprint2022arXiv

Reconciling the Results of the z~2 MOSDEF and KBSS-MOSFIRE Surveys

The combination of the MOSDEF and KBSS-MOSFIRE surveys represents the largest joint investment of Keck/MOSFIRE time to date, with ~3000 galaxies at 1.4<=z<=3.8, roughly half of which are at z~2. MOSDEF is photometric- and spectroscopic-redshift selected with a rest-optical magnitude limit, while KBSS-MOSFIRE is primarily selected based on rest-UV colors and a rest-UV magnitude limit. Analyzing both surveys in a uniform manner with consistent spectral-energy-distribution (SED) models, we find that the MOSDEF z~2 targeted sample has a higher median M_* and redder rest U-V color than the KBSS-MOSFIRE z~2 targeted sample, and a smaller median SED-based SFR and sSFR (SFR(SED) and sSFR(SED)). Specifically, MOSDEF targeted a larger population of red galaxies with U-V and V-J >=1.25, while KBSS-MOSFIRE contains more young galaxies with intense star formation. Despite these differences in the z~2 targeted samples, the subsets of the surveys with multiple emission lines detected and analyzed in previously published work are much more similar. All median host-galaxy properties with the exception of stellar population age -- i.e., M_*, SFR(SED), sSFR(SED), A_V, and UVJ colors -- agree within the uncertainties. Additionally, when uniform emission-line fitting and stellar Balmer absorption correction techniques are applied, there is no significant offset between the two samples in the [OIII]$λ$5008/H$β$ vs. [NII]$λ$6585/H$α$ diagnostic diagram, in contrast to previously-reported discrepancies. We can now combine the MOSDEF and KBSS-MOSFIRE surveys to form the largest z~2 sample with moderate-resolution rest-optical spectra and construct the fundamental scaling relations of star-forming galaxies during this important epoch.

preprint2022arXiv

The MOSDEF Survey: Implications of the Lack of Evolution in the Dust Attenuation-Mass Relation to z~2

We investigate the relationship between dust attenuation and stellar mass ($M_*$) in star-forming galaxies over cosmic time. For this analysis, we compare measurements from the MOSFIRE Deep Evolution Field (MOSDEF) survey at $z\sim2.3$ and the Sloan Digital Sky Survey (SDSS) at $z\sim0$, augmenting the latter optical dataset with both UV Galaxy Evolution Explorer (GALEX) and mid-infrared Wide-field Infrared Survey Explorer (WISE) photometry from the GALEX-SDSS-WISE Catalog. We quantify dust attenuation using both spectroscopic measurements of H$α$ and H$β$ emission lines, and photometric measurements of the rest-UV stellar continuum. The H$α$/H$β$ ratio is used to determine the magnitude of attenuation at the wavelength of H$α$, $A_{{\rm H}α}$. Rest-UV colors and spectral-energy-distribution fitting are used to estimate $A_{1600}$, the magnitude of attenuation at a rest wavelength of 1600Å. As in previous work, we find a lack of significant evolution in the relation between dust attenuation and $M_*$ over the redshift range $z\sim0$ to $z\sim2.3$. Folding in the latest estimates of the evolution of $M_{\rm dust}$, $({M_{\rm dust}}/{M_{\rm gas}})$, and gas surface density at fixed $M_*$, we find that the expected $M_{\rm dust}$ and dust mass surface density are both significantly higher at $z\sim2.3$ than at $z\sim0$. These differences appear at odds with the lack of evolution in dust attenuation. To explain the striking constancy in attenuation vs. $M_*$, it is essential to determine the relationship between metallicity and $({M_{\rm dust}}/{M_{\rm gas}})$, the dust mass absorption coefficient, and dust geometry, and the evolution of these relations and quantities from $z\sim0$ to $z\sim2.3$.

preprint2022arXiv

The MOSDEF-LRIS Survey: Connection between Galactic-scale Outflows and the Properties of $z$~2 star-forming Galaxies

We investigate the conditions that facilitate galactic-scale outflows using a sample of 155 typical star-forming galaxies at $z$~2 drawn from the MOSFIRE Deep Evolution Field (MOSDEF) survey. The sample includes deep rest-frame UV spectroscopy from the Keck Low-Resolution Imaging Spectrometer (LRIS), which provides spectral coverage of several low-ionisation interstellar (LIS) metal absorption lines and Ly$α$ emission. Outflow velocities are calculated from the centroids of the LIS absorption and/or Ly$α$ emission, as well as the highest-velocity component of the outflow from the blue wings of the LIS absorption lines. Outflow velocities are found to be marginally correlated or independent of galaxy properties, such as star-formation rate (SFR) and star-formation rate surface density ($Σ_{\rm SFR}$). Outflow velocity scales with SFR as a power-law with index 0.24, which suggests that the outflows may be primarily driven by mechanical energy generated by supernovae explosions, as opposed to radiation pressure acting on dusty material. On the other hand, outflow velocity and $Σ_{\rm SFR}$ are not significantly correlated, which may be due to the limited dynamic range of $Σ_{\rm SFR}$ probed by our sample. The relationship between outflow velocity and $Σ_{\rm SFR}$ normalised by stellar mass ($Σ_{\rm sSFR}$), as a proxy for gravitational potential, suggests that strong outflows (e.g., > 200 km s$^{-1}$) appear ubiquitous above a threshold of log($Σ_{\rm sSFR}/\rm{yr}^{-1}\ \rm{kpc}^{-2}$) ~ -11.3, and that above this threshold, outflow velocity uncouples from $Σ_{\rm sSFR}$. These results highlight the need for higher resolution spectroscopic data and spatially resolved imaging to test the driving mechanisms of outflows predicted by theory.

preprint2021arXiv

The Effects of Stellar Population and Gas Covering Fraction on the Emergent Lyman Alpha Emission of High-Redshift Galaxies

We perform joint modeling of the composite rest-frame far-UV (FUV) and optical spectra of redshift 1.85<z<3.49 star-forming galaxies to deduce key properties of the massive stars, ionized ISM, and neutral ISM, with the aim of investigating the principal factors affecting the production and escape of Ly-alpha (Lya) photons. Our sample consists of 136 galaxies with deep Keck/LRIS and MOSFIRE spectra covering, respectively, Ly-beta through CIII] 1907, 1909; and [OII], [NeIII], H-beta, [OIII], H-alpha, [NII], and [SII]. Spectral and photoionization modeling indicate that the galaxies are uniformly consistent with stellar population synthesis models that include the effects of stellar binarity. Over the dynamic range of our sample, there is little variation in stellar and nebular abundance with Lya equivalent width, W(Lya), and only a marginal anti-correlation between age and W(Lya). The inferred range of ionizing spectral shapes is insufficient to solely account for the variation in W(Lya). Rather, the covering fraction of optically-thick HI appears to be the principal factor modulating the escape of Lya, with most of the Lya photons in down-the-barrel observations of galaxies escaping through low-column-density or ionized channels in the ISM. Our analysis shows that a high star-formation-rate surface density, Sigma_SFR, particularly when coupled with a low galaxy potential (i.e., low stellar mass), can aid in reducing the covering fraction and ease the escape of Lya photons. We conclude with a discussion of the implications of our results for the escape of ionizing radiation at high redshift.

preprint2021arXiv

The MOSDEF Survey: A Comprehensive Analysis of the Rest-optical Emission-line Properties of $z\sim 2.3$ Star-forming Galaxies

We analyze the rest-optical emission-line spectra of $z\sim2.3$ star-forming galaxies in the complete MOSFIRE Deep Evolution Field (MOSDEF) survey. In investigating the origin of the well-known offset between the sequences of high-redshift and local galaxies in the [O III]5008/H$β$ vs. [N II]6585/H$α$ ("[N II] BPT") diagram, we define two populations of $z\sim2.3$ MOSDEF galaxies. These include the "high" population that is offset towards higher [O III]5008/H$β$ and/or [N II]6585/H$α$ with respect to the local SDSS sequence and the "low" population that overlaps the SDSS sequence. These two groups are also segregated within the [O III]5008/H$β$ vs. [S II]6718,6733/H$α$ and the [O III]4960,5008/[O II]3727,3730 (O$_{32}$) vs. ([O III]4960,5008+[O II]3727,3730)/H$β$ (R$_{23}$) diagram, which suggests qualitatively that star-forming regions in the more offset galaxies are characterized by harder ionizing spectra at fixed nebular oxygen abundance. We also investigate many galaxy properties of the split sample and find that the "high" sample is on average smaller in size and less massive, but has higher specific star-formation rate and star-formation-rate surface density values and is slightly younger compared to the "low" population. From Cloudy+BPASS photoionization models, we estimate that the "high" population has a lower stellar metallicity (i.e., harder ionizing spectrum) but slightly higher nebular metallicity and higher ionization parameter compared to the "low" population. While the "high" population is more $α$-enhanced (i.e., higher $α$/Fe) than the "low" population, both samples are significantly more $α$-enhanced compared to local star-forming galaxies with similar rest-optical line ratios. These differences must be accounted for in all high-redshift star-forming galaxies -- not only those "offset" from local excitation sequences.

preprint2021arXiv

The MOSDEF Survey: Environmental dependence of the gas-phase metallicity of galaxies at $1.4 \leq z \leq 2.6$

Using the near-IR spectroscopy of the MOSFIRE Deep Evolution Field (MOSDEF) survey, we investigate the role of local environment in the gas-phase metallicity of galaxies. The local environment measurements are derived from accurate and uniformly calculated photometric redshifts with well-calibrated probability distributions. Based on rest-frame optical emission lines, [NII]$\lambda6584$ and H$α$, we measure gas-phase oxygen abundance of 167 galaxies at $1.37\leq z\leq1.7$ and 303 galaxies at $2.09\leq z\leq2.61$, located in diverse environments. We find that at $z\sim1.5$, the average metallicity of galaxies in overdensities with $M_*\sim10^{9.8}M_\odot, 10^{10.2}M_\odot$ and $10^{10.8}M_\odot$ is higher relative to their field counterparts by $0.094\pm0.051$, $0.068\pm0.028$ and $0.052\pm0.043$ dex, respectively. However, this metallicity enhancement does not exist at higher redshift, $z\sim2.3$, where, compared to the field galaxies, we find $0.056\pm0.043$, $0.056\pm0.028$ and $0.096\pm 0.034$ dex lower metallicity for galaxies in overdense environments with $M_*\sim10^{9.8}M_\odot, 10^{10.2}M_\odot$ and $10^{10.7}M_\odot$, respectively. Our results suggest that, at $1.37\leq z\leq2.61$, the variation of mass-metallicity relation with local environment is small ($<0.1$dex), and reverses at $z\sim2$. Our results support the hypothesis that, at the early stages of cluster formation, owing to efficient gas cooling, galaxies residing in overdensities host a higher fraction of pristine gas with prominent primordial gas accretion, which lowers their gas-phase metallicity compared to their coeval field galaxies. However, as the Universe evolves to lower redshifts ($z\lesssim2$), shock-heated gas in overdensities cannot cool down efficiently, and galaxies become metal-rich rapidly due to the suppression of pristine gas inflow and re-accretion of metal-enriched outflows in overdensities.

preprint2020arXiv

A comprehensive study of H$α$ emitters at $z \sim$ 0.62 in the DAWN survey: the need for deep and wide regions

We present new estimates of the luminosity function (LF) and star formation rate density (SFRD) for an H$α$ selected sample at $z\sim0.62$ from the Deep And Wide Narrow-band (DAWN) survey. Our results are based on a new H$α$ sample in the extended COSMOS region (compared to Coughlin et al. 2018) with the inclusion of flanking fields, resulting in a total area coverage of $\sim$1.5 deg$^2$. A total of 241 H$α$ emitters were selected based on robust selection criteria using spectro-photometric redshifts and broadband color-color classification. We explore the effect of different dust correction prescriptions by calculating the LF and SFRD using a constant dust extinction correction, A{$_{\textrm{H}α}=1$} mag, a luminosity-dependent correction, and a stellar-mass dependent correction. The resulting H$α$ LFs are well fitted using Schechter functions with best-fit parameters: L$^*=10^{42.24}$ erg s$^{-1}$, $ϕ^*=10^{-2.85}$ Mpc$^{-3}$, $α= -1.62$ for constant dust correction, L$^*=10^{42.31}$ erg s$^{-1}$, $ϕ^*=10^{-2.8}$ Mpc$^{-3}$, $α=-1.39$ for luminosity-dependent dust correction, and L$^*=10^{42.36}$ erg s$^{-1}$, $ϕ^*=10^{-2.91}$ Mpc$^{-3}$, $α= -1.48$, for stellar mass-dependent dust correction. The deep and wide nature of the DAWN survey effectively samples H$α$ emitters over a wide range of luminosities, thereby providing better constraints on both the faint and bright end of the LF. Also, the SFRD estimates $ρ_{\textrm{SFR}}=10^{-1.39}$ M$_{\odot}$yr$^{-1}$Mpc$^{-3}$ (constant dust correction), $ρ_{\textrm{SFR}}=10^{-1.47}$ M$_{\odot}$yr$^{-1}$Mpc$^{-3}$ (luminosity-dependent dust correction), and $ρ_{\textrm{SFR}}=10^{-1.49}$ M$_{\odot}$yr$^{-1}$Mpc$^{-3}$ (stellar mass-dependent dust correction) are in good agreement with the evolution of SFRD across redshifts ($0 < z < 2$) seen from previous H$α$ surveys.

preprint2020arXiv

Bridging between the integrated and resolved main sequence of star formation

The position of galaxies on the stellar mass, star formation rate plane with respect to the star-forming main sequence at each redshift is a convenient way to infer where the galaxy is in its evolution compared to the rest of the population. We use Hubble Space Telescope high resolution images in the GOODS-S field from the CANDELS survey and fit multi wavelength lights in resolution elements of galaxies with stellar population synthesis models. We then construct resolved kpc-scale stellar mass, star formation rate surface density curves for galaxies at z~1. Fitting these resolved main sequence curves with Schechter functions, we parameterize and explain the multi-wavelength structure of galaxies with three variables: phi*, alpha, and M*. For quenched galaxies below the main sequence, we find an average high mass slope (alpha) of the resolved main sequence curves to be ~ -0.4. The scatter of this slope is higher among the lower mass star forming galaxies and those above the main sequence compared to quenched galaxies, due to lack of an evolved bulge. Our findings agree well with an inside-out quenching of star-formation. We find that the knee of the Schechter fits (M*) for galaxies below the main sequence occurs at lower stellar mass surface densities compared to star forming galaxies, which hints at how far quenching has proceeded outwards.

preprint2020arXiv

Dependence of the IRX-$β$ dust attenuation relation on metallicity and environment

We use a sample of star-forming field and protocluster galaxies at z=2.0-2.5 with Keck/MOSFIRE K-band spectra, a wealth of rest-frame UV photometry, and Spitzer/MIPS and Herschel/PACS observations, to dissect the relation between the ratio of IR to UV luminosity (IRX) versus UV slope ($β$) as a function of gas-phase metallicity (12+log(O/H)~8.2-8.7). We find no significant dependence of the IRX-$β$ trend on environment. However, we find that at a given $β$, IRX is highly correlated with metallicity, and less correlated with mass, age, and sSFR. We conclude that, of the physical properties tested here, metallicity is the primary physical cause of the IRX-$β$ scatter, and the IRX correlation with mass is presumably due to the mass dependence on metallicity. Our results indicate that the UV attenuation curve steepens with decreasing metallicity, and spans the full range of slope possibilities from a shallow Calzetti-type curve for galaxies with the highest metallicity in our sample (12+log(O/H)~8.6) to a steep SMC-like curve for those with 12+log(O/H)~8.3. Using a Calzetti (SMC) curve for the low (high) metallicity galaxies can lead to up to a factor of 3 overestimation (underestimation) of the UV attenuation and obscured SFR. We speculate that this change is due to different properties of dust grains present in the ISM of low- and high-metallicity galaxies.

preprint2020arXiv

Effects of Stellar Feedback on Stellar and Gas Kinematics of Star-Forming Galaxies at 0.6<z<1.0

Recent zoom-in cosmological simulations have shown that stellar feedback can flatten the inner density profile of the dark matter halo in low-mass galaxies. A correlation between the stellar/gas velocity dispersion ($σ_{star}$, $σ_{gas}$) and the specific star formation rate (sSFR) is predicted as an observational test of the role of stellar feedback in re-shaping the dark matter density profile. In this work we test the validity of this prediction by studying a sample of star-forming galaxies at $0.6<z<1.0$ from the LEGA-C survey, which provides high signal-to-noise measurements of stellar and gas kinematics. We find that a weak but significant correlation between $σ_{star}$ (and $σ_{gas}$) and sSFR indeed exists for galaxies in the lowest mass bin (M$_\ast\sim10^{10}\,$M$_\odot$). This correlation, albeit with a $\sim$35% scatter, holds for different tracers of star formation, and becomes stronger with redshift. This result generally agrees with the picture that at higher redshifts star formation rate was generally higher, and galaxies at M$_\ast\lesssim10^{10}\,$M$_\odot$ have not yet settled into a disk. As a consequence, they have shallower gravitational potentials more easily perturbed by stellar feedback. The observed correlation between $σ_{star}$ (and $σ_{gas}$) and sSFR supports the scenario predicted by cosmological simulations, in which feedback-driven outflows cause fluctuations in the gravitation potential that flatten the density profiles of low-mass galaxies.

preprint2020arXiv

Resolved Multi-element Stellar Chemical Abundances in the Brightest Quiescent Galaxy at z $\sim$ 2

Measuring the chemical composition of galaxies is crucial to our understanding of galaxy formation and evolution models. However, such measurements are extremely challenging for quiescent galaxies at high redshifts, which have faint stellar continua and compact sizes, making it difficult to detect absorption lines and nearly impossible to spatially resolve them. Gravitational lensing offers the opportunity to study these galaxies with detailed spectroscopy that can be spatially resolved. In this work, we analyze deep spectra of MRG-M0138, a lensed quiescent galaxy at z = 1.98 which is the brightest of its kind, with an H-band magnitude of 17.1. Taking advantage of full spectral fitting, we measure $[{\rm Mg/Fe}]=0.51\pm0.05$, $[\rm{Fe/H}]=0.26\pm0.04$, and, for the first time, the stellar abundances of 6 other elements in this galaxy. We further constrained, also for the first time in a $z\sim2$ galaxy, radial gradients in stellar age, [Fe/H], and [Mg/Fe]. We detect no gradient in age or [Mg/Fe] and a slightly negative gradient in [Fe/H], which has a slope comparable to that seen in local early-type galaxies. Our measurements show that not only is MRG-M0138 very Mg-enhanced compared to the centers of local massive early-type galaxies, it is also very iron rich. These dissimilar abundances suggest that even the inner regions of massive galaxies have experienced significant mixing of stars in mergers, in contrast to a purely inside-out growth model. The abundance pattern observed in MRG-M0138 challenges simple galactic chemical evolution models that vary only the star formation timescale and shows the need for more elaborate models.

preprint2020arXiv

Searching for z > 6.5 Analogs Near the Peak of Cosmic Star Formation

Strong [OIII]$λλ$4959,5007+H$β$ emission appears to be typical in star-forming galaxies at z > 6.5. As likely contributors to cosmic reionization, these galaxies and the physical conditions within them are of great interest. At z > 6.5, where Ly$α$ is greatly attenuated by the intergalactic medium, rest-UV metal emission lines provide an alternative measure of redshift and also constraints on the physical properties of star-forming regions and massive stars. We present the first statistical sample of rest-UV line measurements in z $\sim$ 2 galaxies selected as analogs of those in the reionization era based on [OIII]$λλ$4959,5007 EW or rest-frame U-B color. Our sample is drawn from the 3D-HST Survey and spans the redshift range 1.36 $\leqslant$ z $\leqslant$ 2.49. We find that the median Ly$α$ and CIII]$λλ$1907,1909 EWs of our sample are significantly greater than those of z $\sim$ 2 UV-continuum-selected star-forming galaxies. Measurements from both individual and composite spectra indicate a monotonic, positive correlation between CIII] and [OIII], while a lack of trend is observed between Ly$α$ and [OIII] at [OIII] EW < 1000$\unicode{x212B}$. At higher [OIII] EW, extreme Ly$α$ emission starts to emerge. Using stacked spectra, we find that Ly$α$ and CIII] are significantly enhanced in galaxies with lower metallicity. Two objects in our sample appear comparable to z > 6.5 galaxies with exceptionally strong rest-UV metal line emission. These objects have significant CIV$λλ$1548,1550, HeII$λ$1640, and OIII]$λλ$1661,1665 emission in addition to intense Ly$α$ or CIII]. Detailed characterization of these lower-redshift analogs provides unique insights into the physical conditions in z > 6.5 star-forming regions, motivating future observations of reionization-era analogs at lower redshifts.

preprint2020arXiv

Selection of massive evolved galaxies at $3 \leq z \leq 4.5$ in the CANDELS fields

Using the CANDELS photometric catalogs for the HST/ACS and WFC3, we identified massive evolved galaxies at $3 < z < 4.5$, employing three different selection methods. We find the comoving number density of these objects to be $\sim 2 \times 10^{-5}$ and $8 \times 10^{-6}Mpc^{-3}$ after correction for completeness for two redshift bins centered at $z=3.4, 4.7$. We quantify a measure of how much confidence we should have for each candidate galaxy from different selections and what are the conservative error estimates propagated into our selection. Then we compare the evolution of the corresponding number densities and their stellar mass density with numerical simulations, semi-analytical models, and previous observational estimates, which shows slight tension at higher redshifts as the models tend to underestimate the number and mass densities. By estimating the average halo masses of the candidates ($M_h \approx 4.2, 1.9, 1.3 \times 10^{12} M_\odot$ for redshift bins centered at $z=3.4, 4.1, 4.7$), we find them to be consistent with halos that were efficient in turning baryons to stars and were relatively immune to the feedback effects and on the verge of transition into hot-mode accretion. This can suggest the relative cosmological starvation of the cold gas followed by an overconsumption phase in which the galaxy consumes the available cold gas rapidly as one of the possible drivers for the quenching of the massive evolved population at high redshift.

preprint2020arXiv

Similar Scaling Relations for the Gas Content of Galaxies across Environments to z ~ 3.5

We study the effects of the local environment on the molecular gas content of a large sample of log($M_{*}$/$M_{\odot}$) $\gtrsim$ 10 star-forming and starburst galaxies with specific star-formation rates (sSFRs) on and above the main sequence (MS) to $z$ $\sim$ 3.5. ALMA observations of the dust continuum in the COSMOS field are used to estimate molecular gas masses at $z$ $\approx$ 0.5-3.5. We also use a local universe sample from the ALFALFA HI survey after converting it into molecular masses. The molecular mass ($M_{ISM}$) scaling relation shows a dependence on $z$, $M_{*}$, and sSFR relative to the MS, but no dependence on environmental overdensity $Δ$ ($M_{ISM}$ $\propto$ $Δ^{0.03}$). Similarly, gas mass fraction (f$_{gas}$) and depletion timescale ($τ$) show no environmental dependence to $z$ $\sim$ 3.5. At $\langle z\rangle$ $\sim$ 1.8, the average $\langle M_{ISM}\rangle$,$\langle$f$_{gas}\rangle$, and $\langle τ\rangle$ in densest regions is (1.6$\pm$0.2)$\times$10$^{11}$ $M_{\odot}$, 55$\pm$2%, and 0.8$\pm$0.1 Gyr, respectively, similar to those in the lowest density bin. Independent of the environment, f$_{gas}$ decreases and $τ$ increases with increasing cosmic time. Cosmic molecular mass density ($ρ$) in the lowest density bins peaks at $z$ $\sim$ 1-2, and this peak happens at $z$ $<$ 1 in densest bins. This differential evolution of $ρ$ across environments is likely due to the growth of the large-scale structure with cosmic time. Our results suggest that the molecular gas content and the subsequent star-formation activity of log($M_{*}$/$M_{\odot}$) $\gtrsim$ 10 star-forming and starburst galaxies is primarily driven by internal processes, and not by their local environment since $z$ $\sim$ 3.5.

preprint2020arXiv

Spectroscopic Confirmation of a Coma Cluster Progenitor at z ~ 2.2

We report the spectroscopic confirmation of a new protocluster in the COSMOS field at $z$ $\sim$ 2.2, COSMOS Cluster 2.2 (CC2.2), originally identified as an overdensity of narrowband selected H$α$ emitting candidates. With only two masks of Keck/MOSFIRE near-IR spectroscopy in both $H$ ($\sim$ 1.47-1.81 $μ$m) and $K$ ($\sim$ 1.92-2.40 $μ$m) bands ($\sim$ 1.5 hour each), we confirm 35 unique protocluster members with at least two emission lines detected with S/N $>$ 3. Combined with 12 extra members from the zCOSMOS-deep spectroscopic survey (47 in total), we estimate a mean redshift and a line-of-sight velocity dispersion of $z_{mean}$=2.23224 $\pm$ 0.00101 and $σ_{los}$=645 $\pm$ 69 km s$^{-1}$ for this protocluster, respectively. Assuming virialization and spherical symmetry for the system, we estimate a total mass of $M_{vir}$ $\sim$ $(1-2) \times$10$^{14}$ $M_{\odot}$ for the structure. We evaluate a number density enhancement of $δ_{g}$ $\sim$ 7 for this system and we argue that the structure is likely not fully virialized at $z$ $\sim$ 2.2. However, in a spherical collapse model, $δ_{g}$ is expected to grow to a linear matter enhancement of $\sim$ 1.9 by $z$=0, exceeding the collapse threshold of 1.69, and leading to a fully collapsed and virialized Coma-type structure with a total mass of $M_{dyn}$($z$=0) $\sim$ 9.2$\times$10$^{14}$ $M_{\odot}$ by now. This observationally efficient confirmation suggests that large narrowband emission-line galaxy surveys, when combined with ancillary photometric data, can be used to effectively trace the large-scale structure and protoclusters at a time when they are mostly dominated by star-forming galaxies.

preprint2020arXiv

The BUFFALO HST Survey

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

preprint2020arXiv

The CatWISE Preliminary Catalog: Motions from ${\it WISE}$ and ${\it NEOWISE}$ Data

CatWISE is a program to catalog sources selected from combined ${\it WISE}$ and ${\it NEOWISE}$ all-sky survey data at 3.4 and 4.6 $μ$m (W1 and W2). The CatWISE Preliminary Catalog consists of 900,849,014 sources measured in data collected from 2010 to 2016. This dataset represents four times as many exposures and spans over ten times as large a time baseline as that used for the AllWISE Catalog. CatWISE adapts AllWISE software to measure the sources in coadded images created from six-month subsets of these data, each representing one coverage of the inertial sky, or epoch. The catalog includes the measured motion of sources in 8 epochs over the 6.5 year span of the data. From comparison to ${\it Spitzer}$, the SNR=5 limits in magnitudes in the Vega system are W1=17.67 and W2=16.47, compared to W1=16.96 and W2=16.02 for AllWISE. From comparison to ${\it Gaia}$, CatWISE positions have typical accuracies of 50 mas for stars at W1=10 mag and 275 mas for stars at W1=15.5 mag. Proper motions have typical accuracies of 10 mas yr$^{-1}$ and 30 mas yr$^{-1}$ for stars with these brightnesses, an order of magnitude better than from AllWISE. The catalog is available in the WISE/NEOWISE Enhanced and Contributed Products area of the NASA/IPAC Infrared Science Archive.

preprint2020arXiv

The MOSDEF Survey: Calibrating the relationship between H$α$ star-formation rate and radio continuum luminosity at $1.4 < z < 2.6$

The observed empirical relation between the star-formation rates (SFR) of low-redshift galaxies and their radio continuum luminosity offers a potential means of measuring SFR in high redshift galaxies that is unaffected by dust obscuration. In this study, we make the first test for redshift evolution in the SFR-radio continuum relation at high redshift using dust-corrected H$α$ SFR. Our sample consists of 178 galaxies from the MOSFIRE Deep Evolution Field (MOSDEF) Survey at $1.4 < z < 2.6$ with rest-frame optical spectroscopy and deep 1.5 GHz radio continuum observations from the Karl G. Jansky Very Large Array (VLA) GOODS North field. Using a stacking analysis we compare the observed radio continuum luminosities with those predicted from the dust-corrected H$α$ SFR assuming a range of $z\sim0$ relations. We find no evidence for a systematic evolution with redshift, when stacking the radio continuum as a function of dust-corrected H$α$ SFR and when stacking both optical spectroscopy and radio continuum as a function of stellar mass. We conclude that locally calibrated relations between SFR and radio continuum luminosity remain valid out to $z\sim 2$.

preprint2020arXiv

The MOSDEF Survey: Kinematic and Structural Evolution of Star-Forming Galaxies at $1.4\leq z\leq 3.8$

We present ionized gas kinematics for 681 galaxies at $z\sim 1.4-3.8$ from the MOSFIRE Deep Evolution Field survey, measured using models which account for random galaxy-slit misalignments together with structural parameters derived from CANDELS Hubble Space Telescope (HST) imaging. Kinematics and sizes are used to derive dynamical masses. Baryonic masses are estimated from stellar masses and inferred gas masses from dust-corrected star formation rates (SFRs) and the Kennicutt-Schmidt relation. We measure resolved rotation for 105 galaxies. For the remaining 576 galaxies we use models based on HST imaging structural parameters together with integrated velocity dispersions and baryonic masses to statistically constrain the median ratio of intrinsic ordered to disordered motion, $V/σ_{V,0}$. We find that $V/σ_{V,0}$ increases with increasing stellar mass and decreasing specific SFR (sSFR). These trends may reflect marginal disk stability, where systems with higher gas fractions have thicker disks. For galaxies with detected rotation we assess trends between their kinematics and mass, sSFR, and baryon surface density ($Σ_{\mathrm{bar},e}$). Intrinsic dispersion correlates most with $Σ_{\mathrm{bar},e}$ and velocity correlates most with mass. By comparing dynamical and baryonic masses, we find that galaxies at $z\sim 1.4-3.8$ are baryon dominated within their effective radii ($R_E$), with Mdyn/Mbaryon increasing over time. The inferred baryon fractions within $R_E$, $f_{\mathrm{bar}}$, decrease over time, even at fixed mass, size, or surface density. At fixed redshift, $f_{\mathrm{bar}}$ does not appear to vary with stellar mass but increases with decreasing $R_E$ and increasing $Σ_{\mathrm{bar},e}$. For galaxies at $z\geq2$, the median inferred baryon fractions generally exceed 100%. We discuss possible explanations and future avenues to resolve this tension.

preprint2020arXiv

The MOSDEF Survey: the Variation of the Dust Attenuation Curve with Metallicity

We derive the UV-optical stellar dust attenuation curve of galaxies at z=1.4-2.6 as a function of gas-phase metallicity. We use a sample of 218 star-forming galaxies, excluding those with very young or heavily obscured star formation, from the MOSFIRE Deep Evolution Field (MOSDEF) survey with H$α$, H$β$, and [NII]$λ6585$ spectroscopic measurements. We constrain the shape of the attenuation curve by comparing the average flux densities of galaxies sorted into bins of dust obscuration using Balmer decrements, i.e., H$α$-to-H$β$ luminosities. The average attenuation curve for the high-metallicity sample (12+log(O/H)>8.5, corresponding to $M_*\gtrsim10^{10.4}\,M_{\odot}$) has a shallow slope, identical to that of the Calzetti local starburst curve, and a significant UV 2175A extinction bump that is $\sim 0.5\times$ the strength of the Milky Way bump. On the other hand, the average attenuation curve of the low-metallicity sample (12+log(O/H) $\sim 8.2-8.5$) has a steeper slope similar to that of the SMC curve, only consistent with the Calzetti slope at the $3σ$ level. The UV bump is not detected in the low-metallicity curve, indicating the relative lack of the small dust grains causing the bump at low metallicities. Furthermore, we find that on average the nebular reddening (E(B-V)) is a factor of 2 times larger than that of the stellar continuum for galaxies with low metallicities, while the nebular and stellar reddening are similar for galaxies with higher metallicities. The latter is likely due to a high surface density of dusty clouds embedding the star forming regions but also reddening the continuum in the high-metallicity galaxies.

preprint2020arXiv

The MOSDEF-LRIS Survey: The Interplay Between Massive Stars and Ionized Gas in High-Redshift Star-Forming Galaxies

We present a joint analysis of rest-UV and rest-optical spectra obtained using Keck/LRIS and Keck/MOSFIRE for a sample of 62 star-forming galaxies at $z\sim2.3$. We divide our sample into 2 bins based on their location in the [OIII]/Hb vs. [NII]/Ha BPT diagram, and perform the first differential study of the rest-UV properties of massive ionizing stars as a function of rest-optical emission-line ratios. Fitting BPASS stellar population synthesis models, including nebular continuum emission, to our rest-UV spectra, we find that high-redshift galaxies offset towards higher [OIII]/Hb and [NII]/Ha have younger ages (log(Age/yr)=$7.20^{+0.57}_{-0.20}$) and lower stellar metallicities ($Z_*=0.0010^{+0.0011}_{-0.0003}$) resulting in a harder ionizing spectrum, compared to the galaxies in our sample that lie on the local BPT star-forming sequence (log(Age/yr)=$8.57^{+0.88}_{-0.84}$, $Z_*=0.0019\pm0.0006$). Additionally, we find that the offset galaxies have an ionization parameter of log(U)=$-3.04^{+0.06}_{-0.11}$ and nebular metallicity of 12+log(O/H)=$8.40^{+0.06}_{-0.07}$, and the non-offset galaxies have log(U)=$-3.11\pm0.08$ and 12+log(O/H)=$8.30^{+0.05}_{-0.06}$. The stellar and nebular metallicities derived for our sample imply that the galaxies offset from the local BPT relation are more alpha-enhanced ($7.28^{+2.52}_{-2.82}$ O/Fe$_\odot$) compared to those consistent with the local sequence ($3.04^{+0.95}_{-0.54}$ O/Fe$_\odot$). However, even galaxies that are entirely consistent with the local nebular excitation sequence appear to be alpha-enhanced -- in contrast with typical local systems. Such differences must be considered when estimating gas-phase O/H at high redshift based on strong emission-line ratios. Specifically, a similarity in the location of high-redshift and local galaxies in the BPT diagram may not be indicative of a similarity in their physical properties.

preprint2020arXiv

The SFR-radius connection: data and implications for wind strength and halo concentration

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.

preprint2019arXiv

Large Scale Structures in the CANDELS Fields: The Role of the Environment in Star Formation Activity

We present a robust method, weighted von Mises kernel density estimation, along with boundary correction to reconstruct the underlying number density field of galaxies. We apply this method to galaxies brighter than $\rm HST/F160w\le 26$ AB mag at the redshift range of $0.4\leq z \leq 5$ in the five CANDELS fields (GOODS-N, GOODS-S, EGS, UDS, and COSMOS). We then use these measurements to explore the environmental dependence of the star formation activity of galaxies. We find strong evidence of environmental quenching for massive galaxies ($\rm M \gtrsim 10^{11} \rm {M}_\odot$) out to $z\sim 3.5$ such that an over-dense environment hosts $\gtrsim 20\%$ more massive quiescent galaxies compared to an under-dense region. We also find that environmental quenching efficiency grows with stellar mass and reaches $\sim 60\%$ for massive galaxies at $z\sim 0.5$. The environmental quenching is also more efficient in comparison to the stellar mass quenching for low mass galaxies ($\rm M \lesssim 10^{10} \rm {M}_\odot$) at low and intermediate redshifts ($z\lesssim 1.2$). Our findings concur thoroughly with the "over-consumption" quenching model where the termination of cool gas accretion (cosmological starvation) happens in an over-dense environment and the galaxy starts to consume its remaining gas reservoir in depletion time. The depletion time depends on the stellar mass and could explain the evolution of environmental quenching efficiency with the stellar mass.

preprint2019arXiv

Spatially resolved properties of galaxies from CANDELS+MUSE: Radial extinction profile and insights on quenching

Studying the internal processes of individual galaxies at kilo-parsec scales is crucial in enhancing our understanding of galaxy formation and evolution processes. In this work, we investigate the distribution of star formation rate (SFR), specific SFR (sSFR), and dust attenuation across individual galaxies for a sample of 32 galaxies selected from the MUSE-Wide Survey at 0.1 $< \textit{z} <$ 0.42 with a dynamic range in stellar masses between $10^{7.7}$ and $10^{10.3}$ M$_{\odot}$. We take advantage of the high spatial resolution of the MUSE integral-field spectrograph and measure reliable spatially resolved H$α$ and H$β$ emission line maps for individual galaxies. We also derive resolved stellar mass, SFR and dust maps using pixel-by-pixel SED fitting on high resolution multi-band HST/ACS and HST/WFC3 data from the CANDELS survey. By combining these, we analyze the radial profile of various physical parameters across these galaxies. We observe a radial dependence in both stellar and nebular color excess profiles peaking at the inner regions of galaxies. We also find the color excess profiles to most strongly correlate with the integrated sSFRs of galaxies. The median sSFR$_{\mathrm{H}α}$ radial profiles of galaxies in our sample show a 0.8 dex increase from the central regions outward. This increase compared to the almost flat median radial profile of sSFR$_{\mathrm{SED}}$, which traces longer timescales of star formation, is in favor of the inside-out quenching of star formation. We bring further evidence for this quenching scenario from the locus of different subregions of galaxies on the SFR-M$_{*}$ and sSFR-M$_{*}$ relations.

preprint2019arXiv

The MOSDEF survey: a census of AGN-driven ionized outflows at $z = 1.4-3.8$

Using data from the MOSFIRE Deep Evolution Field (MOSDEF) survey, we present a census of AGN-driven ionized outflows in a sample of 159 AGNs at $1.4 \le z \le 3.8$. The sample spans AGN bolometric luminosities of $10^{44-47} \mathrm{~erg~s}^{-1}$ and includes both quiescent and star-forming galaxies extending across three orders of magnitude in stellar mass. We identify and characterize outflows from the \hbeta, [OIII], \halpha ~and [NII] emission line spectra. We detect outflows in $17\%$ of the AGNs, seven times more often than in a mass-matched sample of inactive galaxies in MOSDEF. The outflows are fast and galaxy-wide, with velocities of $\sim 400-3500 ~\mathrm{km~s}^{-1}$ and spatial extents of $0.3-11.0$ kpc. The incidence of outflows among AGNs is independent of the stellar mass of the host galaxy, with outflows detected in both star-forming and quiescent galaxies. This suggests that outflows exist across different phases in galaxy evolution. We investigate relations between outflow kinematic, spatial, and energetic properties and both AGN and host galaxy properties. Our results show that AGN-driven outflows are widespread in galaxies along the star-forming main sequence. The mass-loading factors of the outflows are typically $0.1-1$ and increase with AGN luminosity, capable of exceeding unity at $L_\mathrm{AGN} \gtrsim 10^{46.3} \mathrm{~erg~s}^{-1}$. In these more luminous sources the ionized outflow alone is likely sufficient to regulate star formation, and when combined with outflowing neutral and molecular gas may be able to quench star formation in their host galaxies.

preprint2019arXiv

The MOSDEF Survey: The Metallicity Dependence of X-ray Binary Populations at $z\sim2$

Population synthesis models predict that high-mass X-ray binary (HMXB) populations produced in low metallicity environments should be more X-ray luminous, a trend supported by studies of nearby galaxies. This trend may be responsible for the observed increase of the X-ray luminosity ($L_{\mathrm{X}}$) per star formation rate (SFR) with redshift due to the decrease of metallicity ($Z$) at fixed stellar mass as a function of redshift. To test this hypothesis, we use a sample of 79 $z\sim2$ star-forming galaxies with oxygen abundance measurements from the MOSDEF survey, which obtained rest-frame optical spectra for $\sim1500$ galaxies in the CANDELS fields at $1.37<z<3.80$. Using Chandra data from the AEGIS-X Deep, Deep Field North, and Deep Field South surveys, we stack the X-ray data at the galaxy locations in bins of redshift and $Z$ because the galaxies are too faint to be individually detected. In agreement with previous studies, the average $L_{\mathrm{X}}$/SFR of our $z\sim2$ galaxy sample is enhanced by $\approx0.4-0.8$ dex relative to local HMXB $L_{\mathrm{X}}$-SFR scaling relations. Splitting our sample by $Z$, we find that $L_{\mathrm{X}}$/SFR and $Z$ are anti-correlated with 97% confidence. This observed $Z$ dependence for HMXB-dominated galaxies is consistent both with the local $L_{\mathrm{X}}$-SFR-$Z$ relation and a subset of population synthesis models. Although the statistical significance of the observed trends is weak due to the low X-ray statistics, these results constitute the first direct evidence connecting the redshift evolution of $L_{\mathrm{X}}$/SFR and the $Z$ dependence of HMXBs.

preprint2016arXiv

A massive, quiescent, population II galaxy at a redshift of 2.1

Unlike spiral galaxies such as the Milky Way, the majority of the stars in massive elliptical galaxies were formed in a short period early in the history of the Universe. The duration of this formation period can be measured using the ratio of magnesium to iron abundance ([Mg/Fe]), which reflects the relative enrichment by core-collapse and type Ia supernovae. For local galaxies, [Mg/Fe] probes the combined formation history of all stars currently in the galaxy, including younger and metal-poor stars that were added during late-time mergers. Therefore, to directly constrain the initial star-formation period, we must study galaxies at earlier epochs. The most distant galaxy for which [Mg/Fe] had previously been measured is at z~1.4, with [Mg/Fe]=0.45(+0.05,-0.19). A slightly earlier epoch (z~1.6) was probed by stacking the spectra of 24 massive quiescent galaxies, yielding an average [Mg/Fe] of 0.31+/-0.12. However, the relatively low S/N of the data and the use of index analysis techniques for both studies resulted in measurement errors that are too large to allow us to form strong conclusions. Deeper spectra at even earlier epochs in combination with analysis techniques based on full spectral fitting are required to precisely measure the abundance pattern shortly after the major star-forming phase (z>2). Here we report a measurement of [Mg/Fe] for a massive quiescent galaxy at z=2.1. With [Mg/Fe]=0.59+/-0.11, this galaxy is the most Mg-enhanced massive galaxy found so far, having twice the Mg enhancement of similar-mass galaxies today. The abundance pattern of the galaxy is consistent with enrichment exclusively by core-collapse supernovae and with a star-formation timescale of 0.1-0.5 Gyr - characteristics that are similar to population II stars in the Milky Way. With an average past SFR of 600-3000 Msol/yr, this galaxy was among the most vigorous star-forming galaxies in the Universe.

preprint2016arXiv

Effects of Local Environment and Stellar Mass on Galaxy Quenching out to z~3

We study the effects of local environment and stellar mass on galaxy properties using a mass complete sample of quiescent and star-forming systems in the COSMOS field at $z\lesssim$ 3. We show that at $z\lesssim$ 1, the median star-formation rate (SFR) and specific SFR (sSFR) of all galaxies depend on environment, but they become independent of environment at $z\gtrsim$ 1. However, we find that only for \textit{star-forming} galaxies, the median SFR and sSFR are similar in different environments, regardless of redshift and stellar mass. We find that the quiescent fraction depends on environment at $z\lesssim$ 1, and on stellar mass out to $z\sim$ 3. We show that at $z\lesssim$ 1, galaxies become quiescent faster in denser environments and that the overall environmental quenching efficiency increases with cosmic time. Environmental and mass quenching processes depend on each other. At $z\lesssim$ 1, denser environments more efficiently quench galaxies with higher masses (log($M/M_{\odot}$)$\gtrsim$ 10.7), possibly due to a higher merger rate of massive galaxies in denser environments, and that mass quenching is more efficient in denser regions. We show that the overall mass quenching efficiency ($ε_{mass}$) for more massive galaxies (log($M/M_{\odot}$)$\gtrsim$ 10.2) rises with cosmic time until $z\sim$ 1 and flattens out since then. However, for less massive galaxies, the rise in $ε_{mass}$ continues to the present time. Our results suggest that environmental quenching is only relevant at $z\lesssim$ 1, likely a fast process, whereas mass quenching is the dominant mechanism at $z\gtrsim$ 1, with a possible stellar feedback physics.

preprint2016arXiv

Evolution of Intrinsic Scatter in the SFR-Stellar Mass Correlation at 0.5<z<3

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.

preprint2016arXiv

The MOSDEF Survey: Detection of [OIII]$λ$4363 and the direct-method oxygen abundance of a star-forming galaxy at z=3.08

We present measurements of the electron-temperature based oxygen abundance for a highly star-forming galaxy at z=3.08, COSMOS-1908. This is the highest redshift at which [OIII]$λ$4363 has been detected, and the first time that this line has been measured at z>2. We estimate an oxygen abundance of 12+log(O/H)$=8.00^{+0.13}_{-0.14}$. This galaxy is a low-mass ($10^{9.3}$ M$_{\odot}$), highly star-forming ($\sim50$ M$_{\odot}$ yr$^{-1}$) system that hosts a young stellar population ($\sim160$ Myr). We investigate the physical conditions of the ionized gas in COSMOS-1908 and find that this galaxy has a high ionization parameter, little nebular reddening ($E(B-V)_{\rm gas}<0.14$), and a high electron density ($n_e\sim500$ cm$^{-3}$). We compare the ratios of strong oxygen, neon, and hydrogen lines to the direct-method oxygen abundance for COSMOS-1908 and additional star-forming galaxies at z=0-1.8 with [OIII]$λ$4363 measurements, and show that galaxies at z$\sim$1-3 follow the same strong-line correlations as galaxies in the local universe. This agreement suggests that the relationship between ionization parameter and O/H is similar for z$\sim$0 and high-redshift galaxies. These results imply that metallicity calibrations based on lines of oxygen, neon, and hydrogen do not strongly evolve with redshift and can reliably estimate abundances out to z$\sim$3, paving the way for robust measurements of the evolution of the mass-metallicity relation to high redshift.

preprint2016arXiv

The MOSDEF Survey: Dynamical and Baryonic Masses and Kinematic Structures of Star-Forming Galaxies at $1.4 \leq z \leq 2.6$

We present H$α$ gas kinematics for 178 star-forming galaxies at z~2 from the MOSFIRE Deep Evolution Field survey. We have developed models to interpret the kinematic measurements from fixed-angle multi-object spectroscopy, using structural parameters derived from CANDELS HST/F160W imaging. For 35 galaxies we measure resolved rotation with a median $(V/σ_{V,0})_{R_E}=2.1$. We derive dynamical masses from the kinematics and sizes and compare them to baryonic masses, with gas masses estimated from dust-corrected H$α$ star formation rates (SFRs) and the Kennicutt-Schmidt relation. When assuming that galaxies with and without observed rotation have the same median $(V/σ_{V,0})_{R_E}$, we find good agreement between the dynamical and baryonic masses, with a scatter of $σ_{RMS}=0.34$dex and a median offset of $Δ\log_{10}M=0.04$dex. This comparison implies a low dark matter fraction (8% within an effective radius) for a Chabrier initial mass function (IMF), and disfavors a Salpeter IMF. Moreover, the requirement that $M_{dyn}/M_{baryon}$ should be independent of inclination yields a median value of $(V/σ_{V,0})_{R_E}=2.1$ for galaxies without observed rotation. If instead we treat the galaxies without detected rotation as early-type galaxies, the masses are also in reasonable agreement ($Δ\log_{10}M=-0.07$dex, $σ_{RMS}=0.37$dex). The inclusion of gas masses is critical in this comparison; if gas masses are excluded there is an increasing trend of $M_{dyn}/M_{*}$ with higher specific SFR (SSFR). Furthermore, we find indications that $V/σ$ decreases with increasing H$α$ SSFR for our full sample, which may reflect disk settling. We also study the Tully-Fisher relation and find that at fixed stellar mass $S_{0.5}=(0.5V_{2.2}^2+σ_{V,0}^2)^{1/2}$ was higher at earlier times. At fixed baryonic mass, we observe the opposite trend. [abridged]

preprint2016arXiv

The MOSDEF Survey: The strong agreement between H-alpha and UV-to-FIR star formation rates for z~2 star-forming galaxies

We present the first direct comparison between Balmer line and panchromatic SED-based SFRs for z~2 galaxies. For this comparison we used 17 star-forming galaxies selected from the MOSFIRE Deep Evolution Field (MOSDEF) survey, with $3σ$ detections for H$α$ and at least two IR bands (Spitzer/MIPS 24$μ$m and Herschel/PACS 100 and 160$μ$m, and in some cases Herschel/SPIRE 250, 350, and 500$μ$m). The galaxies have total IR (8-1000$μ$m) luminosities of $\sim10^{11.4}-10^{12.4}\,\textrm{L}_\odot$ and star-formation rates (SFRs) of $\sim30-250\,\textrm{M}_\odot\,\mathrm{yr^{-1}}$. We fit the UV-to-far-IR SEDs with flexible stellar population synthesis (FSPS) models - which include both stellar and dust emission - and compare the inferred SFRs with the SFR(H$α$,H$β$) values corrected for dust attenuation using Balmer decrements. The two SFRs agree with a scatter of 0.17 dex. Our results imply that the Balmer decrement accurately predicts the obscuration of the nebular lines and can be used to robustly calculate SFRs for star-forming galaxies at z~2 with SFRs up to $\sim200\,\textrm{M}_\odot\,\mathrm{yr^{-1}}$. We also use our data to assess SFR indicators based on modeling the UV-to-mid-IR SEDs or by adding SFR(UV) and SFR(IR), for which the latter is based on the mid-IR only or on the full IR SED. All these SFRs show a poorer agreement with SFR(H$α$,H$β$) and in some cases large systematic biases are observed. Finally, we show that the SFR and dust attenuation derived from the UV-to-near-IR SED alone are unbiased when assuming a delayed exponentially declining star-formation history.

preprint2016arXiv

The Nature of H$β+$[O{\sc iii}] and [O{\sc ii}] emitters to $z \sim 5$ with HiZELS: stellar mass functions and the evolution of EWs

We investigate the properties of $\sim7000$ narrow-band selected galaxies with strong H$β+$[OIII] and [OII] nebular emission lines from the High-$z$ Emission Line Survey (HiZELS) between $z \sim 0.8 - 5.0$. Our sample covers a wide range in stellar mass ($M\sim10^{7.5 - 12.0}$ M$_\odot$), rest-frame equivalent widths (EW$_\mathrm{rest}\sim 10 - 10^5$ Å), and line luminosities ($L\sim10^{40.5 - 43.2}$ erg s$^{-1}$). We measure the H$β+$[OIII] and [OII]-selected stellar mass functions out to $z \sim 3.5$ and find that both $M_\star$ and $ϕ_\star$ increases with cosmic time, which may be due to the [OIII] selection including an increasing fraction of AGN at lower redshifts. The [OII]-selected stellar mass functions show a constant $M_\star\sim10^{11.6}$ M$_\odot$ and a strong, increasing evolution with cosmic time in $ϕ_\star$ in line with H$α$ studies. We also investigate the EW$_\mathrm{rest}$ evolution as a function of redshift with a fixed mass range (10$^{9.5 - 10.0}$ M$_\odot$) and find an increasing trend best represented by $(1+z)^{3.81\pm0.14}$ and $(1+z)^{2.72\pm0.19}$ up to $z\sim 2$ and $\sim 3$ for H$β$+[OIII] and [OII] emitters, respectively. This is the first time that the EW evolution has been directly measured for H$β+$[OIII] and [OII] emitters up to these redshifts. There is evidence for a slower evolution for $z > 2$ in the H$β+$[OIII] EW and a decreasing trend for $z > 3$ in the [OII] EW evolution, which would imply low [OII] EW at the highest redshifts and higher [OIII]/[OII] line ratios. This suggests that the ionization parameter at higher redshift may be significantly higher than the local Universe. Our results set the stage for future near-IR space-based spectroscopic surveys to test our extrapolated predictions and also produce $z > 5$ measurements to constrain the high-$z$ end of the EW$_\mathrm{rest}$ and [OIII]/[OII] evolution.

preprint2015arXiv

A Comparative Study of Density Field Estimation for Galaxies: New Insights into the Evolution of Galaxies with Environment in COSMOS out to z~3

It is well-known that galaxy environment has a fundamental effect in shaping its properties. We study the environmental effects on galaxy evolution, with an emphasis on the environment defined as the local number density of galaxies. The density field is estimated with different estimators (weighted adaptive kernel smoothing, 10$^{th}$ and 5$^{th}$ nearest neighbors, Voronoi and Delaunay tessellation) for a K$_{s}<$24 sample of $\sim$190,000 galaxies in the COSMOS field at 0.1$<$z$<$3.1. The performance of each estimator is evaluated with extensive simulations. We show that overall, there is a good agreement between the estimated density fields using different methods over $\sim$2 dex in overdensity values. However, our simulations show that adaptive kernel and Voronoi tessellation outperform other methods. Using the Voronoi tessellation method, we assign surface densities to a mass complete sample of quiescent and star-forming galaxies out to z$\sim$3. We show that at a fixed stellar mass, the median color of quiescent galaxies does not depend on their host environment out to z$\sim$3. We find that the number and stellar mass density of massive ($>$10$^{11}$M$_{\odot}$) star-forming galaxies have not significantly changed since z$\sim$3, regardless of their environment. However, for massive quiescent systems at lower redshifts (z$\lesssim$1.3), we find a significant evolution in the number and stellar mass densities in denser environments compared to lower density regions. Our results suggest that the relation between stellar mass and local density is more fundamental than the color-density relation and that environment plays a significant role in quenching star formation activity in galaxies at z$\lesssim$1.

preprint2015arXiv

A Correlation Between Lyα Spectral Line Profile and Rest-Frame UV Morphology

We explore the relationship between the spectral shape of the Lyα emission and the UV morphology of the host galaxy using a sample of 304 Lyα-emitting BV i-dropouts at 3 < z < 7 in the GOODS and COSMOS fields. Using our extensive reservoir of high-quality Keck DEIMOS spectra combined with HST WFC3 data, we measure the Lyα line asymmetries for individual galaxies and compare them to axial ratios measured from observed J- and H-band (restframe UV) images. We find that the Lyα skewness exhibits a large scatter at small elongation (a/b < 2), and this scatter decreases as axial ratio increases. Comparison of this trend to radiative transfer models and various results from literature suggests that these high-redshift Lyα emitters are not likely to be intrinsically round and symmetric disks, but they probably host galactic outflows traced by Lyα emitting clouds. The ionizing sources are centrally located, with the optical depth a good indicator of the absorption and scattering events on the escape path of Lyα photons from the source. Our results find no evidence for evolution in Lyα asymmetry or axial ratio with look-back time.

preprint2015arXiv

A Critical Assessment of Stellar Mass Measurement Methods

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, $Δ\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.

preprint2015arXiv

A WFC3 Grism Emission Line Redshift Catalog in the GOODS-South Field

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.

preprint2015arXiv

CANDELS Visual Classifications: Scheme, Data Release, and First Results

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.

preprint2015arXiv

CF-HiZELS, a 10 deg$^2$ emission-line survey with spectroscopic follow-up: Hα, [OIII]+Hβ and [OII] luminosity functions at z=0.8, 1.4 and 2.2

We present results from the largest contiguous narrow-band survey in the near-infrared. We have used WIRCam/CFHT and the lowOH2 filter (1.187$\pm$0.005 $μ$m) to survey ~10 deg$^2$ of contiguous extragalactic sky in the SA22 field. A total of ~6000 candidate emission-line galaxies are found. We use deep ugrizJK data to obtain robust photometric redshifts. We combine our data with the High-redshift Emission Line Survey (HiZELS), explore spectroscopic surveys (VVDS, VIPERS) and obtain our own spectroscopic follow-up with KMOS, FMOS and MOSFIRE to derive large samples of high-redshift emission-line selected galaxies: 3471 H$α$ emitters at z=0.8, 1343 [OIII]+H$β$ emitters at z=1.4 and 572 [OII] emitters at z=2.2. We probe co-moving volumes of >10$^6$ Mpc$^3$ and find significant over-densities, including an 8.5$σ$ (spectroscopically confirmed) over-density of H$α$ emitters at z=0.81. We derive H$α$, [OIII]+H$β$ and [OII] luminosity functions at z=0.8,1.4,2.2, respectively, and present implications for future surveys such as Euclid. Our uniquely large volumes/areas allow us to sub-divide the samples in thousands of randomised combinations of areas and provide a robust empirical measurement of sample/cosmic variance. We show that surveys for star-forming/emission-line galaxies at a depth similar to ours can only overcome cosmic-variance (errors <10%) if they are based on volumes >5x10$^{5}$ Mpc$^{3}$; errors on $L^*$ and $ϕ^*$ due to sample (cosmic) variance on surveys probing ~10$^4$ Mpc$^{3}$ and ~10$^5$ Mpc$^{3}$ are typically very high: ~300% and ~40-60%, respectively.

preprint2015arXiv

Evidence for PopIII-like stellar populations in the most luminous Lyman-$α$ emitters at the epoch of re-ionisation: spectroscopic confirmation

Faint Lyman-$α$ (Ly$α$) emitters become increasingly rarer towards the re-ionisation epoch (z~6-7). However, observations from a very large (~5deg$^2$) Ly$α$ survey at z=6.6 (Matthee et al. 2015) show that this is not the case for the most luminous emitters. Here we present follow-up observations of the two most luminous z~6.6 Ly$α$ candidates in the COSMOS field: `MASOSA' and `CR7'. We used X-SHOOTER, SINFONI and FORS2 (VLT), and DEIMOS (Keck), to confirm both candidates beyond any doubt. We find redshifts of z=6.541 and z=6.604 for MASOSA and CR7, respectively. MASOSA has a strong detection in Ly$α$ with a line width of $386\pm30$ km/s (FWHM) and with high EW$_0$ (>200 Å), but it is undetected in the continuum. CR7, with an observed Ly$α$ luminosity of $10^{43.93\pm0.05}$erg/s is the most luminous Ly$α$ emitter ever found at z>6. CR7 reveals a narrow Ly$α$ line with $266\pm15$ km/s FWHM, being detected in the NIR (rest-frame UV, with $β=-2.3\pm0.1$) with an excess in $J$, and also strongly detected in IRAC/Spitzer. We detect a narrow HeII1640$Å$ emission line ($6σ$) which explains the excess seen in the $J$ band photometry (EW$_0$~80 Å). We find no other emission lines from the UV to the NIR in our X-SHOOTER spectra (HeII/OIII]1663A>3 and HeII/CIII]1908A>2.5). We find that CR7 is best explained by a combination of a PopIII-like population which dominates the rest-frame UV and the nebular emission, and a more normal stellar population which dominates the mass. HST/WFC3 observations show that the light is indeed spatially separated between a very blue component, coincident with Ly$α$ and HeII emission, and two red components (~5 kpc away), which dominate the mass. Our findings are consistent with theoretical predictions of a PopIII wave, with PopIII star formation migrating away from the original sites of star formation.

preprint2015arXiv

Evolution of the H$β$+[OIII] and [OII] luminosity functions and the [OII] star-formation history of the Universe up to $z$ ~ 5 from HiZELS

We investigate the evolution of the H$β$+[OIII] and [OII] luminosity functions from $z \sim 0.8$ to $\sim5$ in four redshift slices per emission line using data from the High-{\it z} Emission Line Survey (HiZELS). This is the first time that the H$β$+[OIII] and [OII] luminosity functions have been studied at these redshifts in a self-consistent analysis. This is also the largest sample of [OII] and H$β$+[OIII] emitters (3475 and 3298 emitters, respectively) in this redshift range, with large co-moving volumes $\sim 1 \times 10^6$ Mpc$^{-3}$ in two independent volumes (COSMOS and UDS), greatly reducing the effects of cosmic variance. The emitters were selected by a combination of photometric redshift and color-color selections, as well as spectroscopic follow-up, including recent spectroscopic observations using DEIMOS and MOSFIRE on the Keck Telescopes and FMOS on Subaru. We find a strong increase in $L_\star$ and a decrease in $ϕ_\star$ for both H$β$+[OIII] and [OII] emitters. We derive the [OII] star-formation history of the Universe since $z\sim5$ and find that the cosmic SFRD rises from $z \sim 5$ to $\sim 3$ and then drops towards $z \sim 0$. We also find that our star-formation history is able to reproduce the evolution of the stellar mass density up to $z\sim 5$ based only on a single tracer of star-formation. When comparing the H$β$+[OIII] SFRDs to the [OII] and H$α$ SFRD measurements in the literature, we find that there is a remarkable agreement, suggesting that the H$β$+[OIII] sample is dominated by star-forming galaxies at high-$z$ rather than AGNs.

preprint2015arXiv

Identification of the brightest Lyα emitters at z=6.6: implications for the evolution of the luminosity function in the re-ionisation era

Using wide field narrow-band surveys, we provide a new measurement of the $z=6.6$ Lyman-$α$ Emitter (LAE) luminosity function (LF), which constraints the bright end for the first time. We use a combination of archival narrow-band NB921 data in UDS and new NB921 measurements in SA22 and COSMOS/UltraVISTA, all observed with the Subaru telescope, with a total area of $\sim 5$ deg$^2$. We exclude lower redshift interlopers by using broad-band optical and near-infrared photometry and also exclude three supernovae with data split over multiple epochs. Combining the UDS and COSMOS samples we find no evolution of the bright end of the Ly$α$ LF between $z=5.7$ and $6.6$, which is supported by spectroscopic follow-up, and conclude that sources with \emph{Himiko}-like luminosity are not as rare as previously thought, with number densities of $\sim 1.5\times10^{-5}$ Mpc$^{-3}$. Combined with our wide-field SA22 measurements, our results indicate a non-Schechter-like bright end of the LF at $z=6.6$ and a different evolution of \emph{observed} faint and bright LAEs. This differential evolution is also seen in the spectroscopic follow-up of UV selected galaxies and is now also confirmed for Ly$α$ emitters, and we argue that it may be an effect of re-ionisation. Using a toy-model, we show that such differential evolution of the LF is expected, since brighter sources are able to ionise their surroundings earlier, such that Ly$α$ photons are able to escape. Our targets are excellent candidates for detailed follow-up studies and provide the possibility to give a unique view on the earliest stages in the formation of galaxies and re-ionisation process.

preprint2015arXiv

Mapping the Galaxy Color-Redshift Relation: Optimal Photometric Redshift Calibration Strategies for Cosmology Surveys

Calibrating the photometric redshifts of >10^9 galaxies for upcoming weak lensing cosmology experiments is a major challenge for the astrophysics community. The path to obtaining the required spectroscopic redshifts for training and calibration is daunting, given the anticipated depths of the surveys and the difficulty in obtaining secure redshifts for some faint galaxy populations. Here we present an analysis of the problem based on the self-organizing map, a method of mapping the distribution of data in a high-dimensional space and projecting it onto a lower-dimensional representation. We apply this method to existing photometric data from the COSMOS survey selected to approximate the anticipated Euclid weak lensing sample, enabling us to robustly map the empirical distribution of galaxies in the multidimensional color space defined by the expected Euclid filters. Mapping this multicolor distribution lets us determine where - in galaxy color space - redshifts from current spectroscopic surveys exist and where they are systematically missing. Crucially, the method lets us determine whether a spectroscopic training sample is representative of the full photometric space occupied by the galaxies in a survey. We explore optimal sampling techniques and estimate the additional spectroscopy needed to map out the color-redshift relation, finding that sampling the galaxy distribution in color space in a systematic way can efficiently meet the calibration requirements. While the analysis presented here focuses on the Euclid survey, similar analysis can be applied to other surveys facing the same calibration challenge, such as DES, LSST, and WFIRST.

preprint2015arXiv

Nebular and Stellar Dust Extinction Across the Disk of Emission-Line Galaxies on Small (KPC) Scales

We investigate resolved kpc-scale stellar and nebular dust distribution in eight star-forming galaxies at z~0.4 in the GOODS fields. Constructing the observed Spectral Energy Distributions (SEDs) per pixel, based on seven bands photometric data from HST/ACS and WFC3, we performed pixel-by-pixel SED fits to population synthesis models and estimated small-scale distribution of stellar dust extinction. We use Halpha / Hbeta nebular emission line ratios from Keck/DEIMOS high resolution spectra at each spatial resolution element to measure the amount of attenuation faced by ionized gas at different radii from the center of galaxies. We find a good agreement between the integrated and median of resolved color excess measurements in our galaxies. The ratio of integrated nebular to stellar dust extinction is always greater than unity, but does not show any trend with stellar mass or star formation rate. We find that inclination plays an important role in the variation of the nebular to stellar excess ratio. The stellar color excess profiles are found to have higher values at the center compared to outer parts of the disk. However, for lower mass galaxies, a similar trend is not found for the nebular color excess. We find that the nebular color excess increases with stellar mass surface density. This explains the absence of radial trend in the nebular color excess in lower mass galaxies which lack a large radial variation of stellar mass surface density. Using standard conversions of star formation rate surface density to gas mass surface density, and the relation between dust mass surface density and color excess, we find no significant variation in the dust to gas ratio in regions with high gas mass surface densities, over the scales probed in this study.

preprint2015arXiv

Probing the Physical Properties of z=4.5 Lyman Alpha Emitters with Spitzer

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.

preprint2015arXiv

Spectroscopic Study of Star-forming Galaxies in Filaments and the Field at $z\sim$0.5: Evidence for Environmental Dependence of Electron Density

We study the physical properties of a spectroscopic sample of 28 star-forming galaxies in a large filamentary structure in the COSMOS field at $z\sim$0.53, with spectroscopic data taken with the Keck/DEIMOS spectrograph, and compare them with a control sample of 30 field galaxies. We spectroscopically confirm the presence of a large galaxy filament ($\sim$ 8 Mpc), along which five confirmed X-ray groups exist. We show that within the uncertainties, the ionization parameter, equivalent width (EW), EW versus specific star-formation rate (sSFR) relation, EW versus stellar mass relation, line-of-sight velocity dispersion, dynamical mass, and stellar-to-dynamical mass ratio are similar for filament and field star-forming galaxies. However, we show that on average, filament star-forming galaxies are more metal-enriched ($\sim$ 0.1$-$0.15 dex), possibly due to the inflow of the already enriched intrafilamentary gas into filament galaxies. Moreover, we show that electron densities are significantly lower (a factor of $\sim$17) in filament star-forming systems compared to those in the field, possibly because of a longer star-formation timescale for filament star-forming galaxies. Our results highlight the potential pre-processing role of galaxy filaments and intermediate-density environments on the evolution of galaxies, which has been highly underestimated.

preprint2015arXiv

The Evolution of the Galaxy Rest-Frame Ultraviolet Luminosity Function Over the First Two Billion Years

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]

preprint2015arXiv

The MOSDEF Survey: Dissecting the star-formation rate vs. stellar mass relation using H$α$ and H$β$ emission lines at z ~ 2

We present results on the star-formation rate (SFR) versus stellar mass ($M_*$) relation (i.e., the "main sequence") among star-forming galaxies at 1.37$\leq$$z$$\leq$2.61 using the MOSFIRE Deep Evolution Field (MOSDEF) survey. Based on a sample of 261 galaxies with H$α$ and H$β$ spectroscopy, we have estimated robust dust-corrected instantaneous SFRs over a large range in $M_*$ ($\sim10^{9.5}-10^{11.5}M_\odot$). We find a correlation between log(SFR(H$α$)) and log($M_*$) with a slope of 0.65$\pm$0.08 (0.58$\pm$0.10) at 1.4<z<2.6 (2.1<z<2.6). We find that different assumptions for the dust correction, such as using the color-excess of the stellar continuum to correct the nebular lines, sample selection biases against red star-forming galaxies, and not accounting for Balmer absorption can yield steeper slopes of the log(SFR)-log($M_*$) relation. Our sample is immune from these biases as it is rest-frame optically selected, H$α$ and H$β$ are corrected for Balmer absorption, and the H$α$ luminosity is dust-corrected using the nebular color-excess computed from the Balmer decrement. The scatter of the log(SFR(H$α$))-log($M_*$) relation, after accounting for the measurement uncertainties, is 0.31 dex at 2.1<z<2.6, which is 0.05 dex larger than the scatter in log(SFR(UV))-log($M_*$). Based on comparisons to a simulated SFR-$M_*$ relation with some intrinsic scatter, we argue that in the absence of direct measurements of galaxy-to-galaxy variations in the attenuation/extinction curves and the IMF, one cannot use the difference in the scatter of the SFR(H$α$)- and SFR(UV)-$M_*$ relations to constrain the stochasticity of star formation in high-redshift galaxies.

preprint2015arXiv

The MOSDEF Survey: Electron Density and Ionization Parameter at $z\sim2.3$

Using observations from the MOSFIRE Deep Evolution Field (MOSDEF) survey, we investigate the physical conditions of star-forming regions in $z\sim2.3$ galaxies, specifically the electron density and ionization state. From measurements of the [O II]$λλ$3726,3729 and [S II]$λλ$6716,6731 doublets, we find a median electron density of $\sim250$ cm$^{-3}$ at $z\sim2.3$, an increase of an order of magnitude compared to measurements of galaxies at $z\sim0$. While $z\sim2.3$ galaxies are offset towards significantly higher O$_{32}$ values relative to local galaxies at fixed stellar mass, we find that the high-redshift sample follows a similar distribution to the low-metallicity tail of the local distribution in the O$_{32}$ vs. R$_{23}$ and O3N2 diagrams. Based on these results, we propose that $z\sim2.3$ star-forming galaxies have the same ionization parameter as local galaxies at fixed metallicity. In combination with simple photoionization models, the position of local and $z\sim2.3$ galaxies in excitation diagrams suggests that there is no significant change in the hardness of the ionizing spectrum at fixed metallicity from $z\sim0$ to $z\sim2.3$. We find that $z\sim2.3$ galaxies show no offset compared to low-metallicity local galaxies in emission line ratio diagrams involving only lines of hydrogen, oxygen, and sulfur, but show a systematic offset in diagrams involving [N II]$λ$6584. We conclude that the offset of $z\sim2.3$ galaxies from the local star-forming sequence in the [N II] BPT diagram is primarily driven by elevated N/O at fixed O/H compared to local galaxies. These results suggest that the local gas-phase and stellar metallicity sets the ionization state of star-forming regions at $z\sim0$ and $z\sim2$.

preprint2015arXiv

The MOSDEF Survey: Excitation Properties of $z\sim 2.3$ Star-forming Galaxies

We present results on the excitation properties of z~2.3 galaxies using early observations from the MOSFIRE Deep Evolution Field (MOSDEF) Survey. With its coverage of the full suite of strong rest-frame optical emission lines, MOSDEF provides an unprecedented view of the rest-frame optical spectra of a representative sample of distant star-forming galaxies. We investigate the locations of z~2.3 MOSDEF galaxies in multiple emission-line diagnostic diagrams. These include the [OIII]/Hb vs. [NII]/Ha and [OIII]/Hb vs. [SII]/Ha "BPT" diagrams, as well as the O_32 vs. R_23 excitation diagram. We recover the well-known offset in the star-forming sequence of high-redshift galaxies in the [OIII]/Hb vs. [NII]/Ha BPT diagram relative to SDSS star-forming galaxies. However, the shift for our rest-frame optically selected sample is less significant than for rest-frame-UV selected and emission-line selected galaxies at z~2. Furthermore, we find that the offset is mass-dependent, only appearing within the low-mass half of the z~2.3 MOSDEF sample, where galaxies are shifted towards higher [NII]/Ha at fixed [OIII]/Hb. Within the [OIII]/Hb vs. [SII]/Ha and O_32 vs. R_23 diagrams, we find that z~2.3 galaxies are distributed like local ones, and therefore attribute the shift in the [OIII]/Hb vs. [NII]/Ha BPT diagram to elevated N/O abundance ratios among lower-mass (M_*<10^10 M_sun) high-redshift galaxies. The variation in N/O ratios calls into question the use at high redshift of oxygen abundance indicators based on nitrogen lines, but the apparent invariance with redshift of the excitation sequence in the O_32 vs. R_23 diagram paves the way for using the combination of O_32 and R_23 as an unbiased metallicity indicator over a wide range in redshift. This indicator will allow for an accurate characterization of the shape and normalization of the mass-metallicity relationship over more than 10 Gyr.

preprint2015arXiv

The MOSDEF Survey: Measurements of Balmer Decrements and the Dust Attenuation Curve at Redshifts z~1.4-2.6

We present results on the dust attenuation curve of z~2 galaxies using early observations from the MOSFIRE Deep Evolution Field (MOSDEF) survey. Our sample consists of 224 star-forming galaxies with nebular spectroscopic redshifts in the range z= 1.36-2.59 and high S/N measurements of, or upper limits on, the H-alpha and H-beta emission lines obtained with Keck/MOSFIRE. We construct composite SEDs of galaxies in bins of specific SFR and Balmer optical depth in order to directly constrain the dust attenuation curve from the UV through near-IR for typical star-forming galaxies at high redshift. Our results imply an attenuation curve that is very similar to the SMC extinction curve at wavelengths redward of 2500 Angstroms. At shorter wavelengths, the shape of the curve is identical to that of the Calzetti relation, but with a lower normalization (R_V). Hence, the new attenuation curve results in SFRs that are ~20% lower, and log stellar masses that are 0.16 dex lower, than those obtained with the Calzetti attenuation curve. Moreover, we find that the difference in the reddening---and the total attenuation---of the ionized gas and stellar continuum correlates strongly with SFR, such that for dust-corrected SFRs larger than 20 Msun/yr assuming a Chabrier IMF, the nebular emission lines suffer an increasing degree of obscuration relative to the continuum. A simple model that can account for these trends is one in which the UV through optical stellar continuum is dominated by a population of less reddened stars, while the nebular line and bolometric luminosities become increasingly dominated by dustier stellar populations for galaxies with large SFRs, as a result of the increased dust enrichment that accompanies such galaxies. Consequently, UV- and SED-based SFRs may underestimate the total SFR at even modest levels of ~20 Msun/yr. [Abridged]

preprint2015arXiv

The MOSDEF Survey: Optical AGN Diagnostics at z~2.3

We present results from the MOSFIRE Deep Evolution Field (MOSDEF) survey on rest-frame optical AGN identification and completeness at z~2.3. With our sample of 50 galaxies and 10 X-ray and IR-selected AGN with measured H-beta, [OIII], H-alpha, and [NII] emission lines, we investigate the location of AGN in the BPT, MEx (mass-excitation), and CEx (color-excitation) diagrams. We find that the BPT diagram works well to identify AGN at z~2.3 and that the z~0 AGN/star-forming galaxy classifications do not need to shift substantially at z~2.3 to robustly separate these populations. However, the MEx diagram fails to identify all of the AGN identified in the BPT diagram, and the CEx diagram is substantially contaminated at high redshift. We further show that AGN samples selected using the BPT diagram have selection biases in terms of both host stellar mass and stellar population, in that AGN in low mass and/or high specific star formation rate galaxies are difficult to identify using the BPT diagram. These selection biases become increasingly severe at high redshift, such that optically-selected AGN samples at high redshift will necessarily be incomplete. We also find that the gas in the narrow-line region appears to be more enriched than gas in the host galaxy for at least some MOSDEF AGN. However, AGN at z~2 are generally less enriched than local AGN with the same host stellar mass.

preprint2015arXiv

The MOSFIRE Deep Evolution Field (MOSDEF) Survey: Rest-Frame Optical Spectroscopy for ~1500 H-Selected Galaxies at 1.37 < z < 3.8

In this paper we present the MOSFIRE Deep Evolution Field (MOSDEF) survey. The MOSDEF survey aims to obtain moderate-resolution (R=3000-3650) rest-frame optical spectra (~3700-7000 Angstrom) for ~1500 galaxies at 1.37<z<3.80 in three well-studied CANDELS fields: AEGIS, COSMOS, and GOODS-N. Targets are selected in three redshift intervals: 1.37<z<1.70, 2.09<z<2.61, and 2.95<z<3.80, down to fixed H_AB (F160W) magnitudes of 24.0, 24.5 and 25.0, respectively, using the photometric and spectroscopic catalogs from the 3D-HST survey. We target both strong nebular emission lines (e.g., [OII], Hbeta, [OIII], 5008, Halpha, [NII], and [SII]) and stellar continuum and absorption features (e.g., Balmer lines, Ca-II H and K, Mgb, 4000 Angstrom break). Here we present an overview of our survey, the observational strategy, the data reduction and analysis, and the sample characteristics based on spectra obtained during the first 24 nights. To date, we have completed 21 masks, obtaining spectra for 591 galaxies. For ~80% of the targets we derive a robust redshift from either emission or absorption lines. In addition, we confirm 55 additional galaxies, which were serendipitously detected. The MOSDEF galaxy sample includes unobscured star-forming, dusty star-forming, and quiescent galaxies and spans a wide range in stellar mass (~10^9-10^11.5 Msol) and star formation rate (~10^0-10^3 Msol/yr). The spectroscopically confirmed sample is roughly representative of an H-band limited galaxy sample at these redshifts. With its large sample size, broad diversity in galaxy properties, and wealth of available ancillary data, MOSDEF will transform our understanding of the stellar, gaseous, metal, dust, and black hole content of galaxies during the time when the universe was most active.

preprint2015arXiv

The Rate of Core Collapse Supernovae to Redshift 2.5 From The CANDELS and CLASH Supernova Surveys

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) and Cluster Lensing And Supernova survey with Hubble (CLASH) multi-cycle treasury programs with the Hubble Space Telescope (HST) have provided new opportunities to probe the rate of core-collapse supernovae (CCSNe) at high redshift, now extending to $z\approx2.5$. Here we use a sample of approximately 44 CCSNe to determine volumetric rates, $R_{CC}$, in six redshift bins in the range $0.1<z<2.5$. Together with rates from our previous HST program, and rates from the literature, we trace a more complete history of $R_{CC}(z)$, with $R_{CC}=0.72\pm0.06$ yr$^{-1}$ Mpc$^{-3}$ 10$^{-4}$ $h_{70}^{3}$ at $z<0.08$, and increasing to $3.7^{+3.1}_{-1.6}$ yr$^{-1}$ Mpc$^{-3}$ 10$^{-4}$ $h_{70}^{3}$ to $z\approx2.0$. The statistical precision in each bin is several factors better than than the systematic error, with significant contributions from host extinction, and average peak absolute magnitudes of the assumed luminosity functions for CCSN types. Assuming negligible time delays from stellar formation to explosion, we find these composite CCSN rates to be in excellent agreement with cosmic star formation rate density (SFRs) derived largely from dust-corrected rest-frame UV emission, with a scaling factor of $k=0.0091\pm0.0017\,M^{-1}_{\odot}$, and inconsistent (to $>95\%$ confidence) with SFRs from IR luminous galaxies, or with SFR models that include simple evolution in the initial mass function over time. This scaling factor is expected if the fraction of the IMF contributing to CCSN progenitors is in the 8 to 50 $M_{\odot}$ range. It is not supportive, however, of an upper mass limit for progenitors at $<20\,M_{\odot}$.

preprint2015arXiv

The Relation Between SFR and Stellar Mass for Galaxies at 3.5 $\le z\le$ 6.5 in CANDELS

Distant star-forming galaxies show a correlation between their star formation rates (SFR) and stellar masses, and this has deep implications for galaxy formation. Here, we present a study on the evolution of the slope and scatter of the SFR-stellar mass relation for galaxies at $3.5\leq z\leq 6.5$ using multi-wavelength photometry in GOODS-S from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) and Spitzer Extended Deep Survey. We describe an updated, Bayesian spectral-energy distribution fitting method that incorporates effects of nebular line emission, star formation histories that are constant or rising with time, and different dust attenuation prescriptions (starburst and Small Magellanic Cloud). From $z$=6.5 to $z$=3.5 star-forming galaxies in CANDELS follow a nearly unevolving correlation between stellar mass and SFR that follows SFR $\sim$ $M_\star^a$ with $a = 0.54 \pm 0.16$ at $z\sim 6$ and $0.70 \pm 0.21$ at $z\sim 4$. This evolution requires a star formation history that increases with decreasing redshift (on average, the SFRs of individual galaxies rise with time). The observed scatter in the SFR-stellar mass relation is tight, $σ(\log \mathrm{SFR}/\mathrm{M}_\odot$ yr$^{-1})< 0.3\ - $ 0.4 dex, for galaxies with $\log M_\star/\mathrm{M}_\odot > 9$ dex. Assuming that the SFR is tied to the net gas inflow rate (SFR $\sim$ $\dot{M}_\mathrm{gas}$), then the scatter in the gas inflow rate is also smaller than 0.3$-$0.4 dex for star-forming galaxies in these stellar mass and redshift ranges, at least when averaged over the timescale of star formation. We further show that the implied star formation history of objects selected on the basis of their co-moving number densities is consistent with the evolution in the SFR-stellar mass relation.

preprint2015arXiv

Two Type Ia Supernovae at Redshift ~2 : Improved Classification and Redshift Determination with Medium-band Infrared Imaging

We present two supernovae (SNe) discovered with the Hubble Space Telescope (HST) in the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS), an HST multi-cycle treasury program. We classify both objects as Type Ia SNe and find redshifts of z = 1.80+-0.02 and 2.26 +0.02 -0.10, the latter of which is the highest redshift Type Ia SN yet seen. Using light curve fitting we determine luminosity distances and find that both objects are consistent with a standard Lambda-CDM cosmological model. These SNe were observed using the HST Wide Field Camera 3 infrared detector (WFC3-IR), with imaging in both wide- and medium-band filters. We demonstrate that the classification and redshift estimates are significantly improved by the inclusion of single-epoch medium-band observations. This medium-band imaging approximates a very low resolution spectrum (lambda/delta lambda ~ 100) which can isolate broad spectral absorption features that differentiate Type Ia SNe from their most common core collapse cousins. This medium-band method is also insensitive to dust extinction and (unlike grism spectroscopy) it is not affected by contamination from the SN host galaxy or other nearby sources. As such, it can provide a more efficient - though less precise - alternative to IR spectroscopy for high-z SNe.

preprint2014arXiv

A Study of Massive and Evolved Galaxies at High Redshift

We use data taken as part of HST/WFC3 observations of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) to identify massive and evolved galaxies at 3<z<4.5. This is performed using the strength of the Balmer break feature at rest-frame 3648A, which is a diagnostic of the age of the stellar population in galaxies. Using WFC3 H-band selected catalog for the CANDELS GOODS-S field and deep multi-waveband photometry from optical (HST) to mid-infrared (Spitzer) wavelengths, we identify a population of old and evolved post-starburst galaxies based on the strength of their Balmer breaks (Balmer Break Galaxies- BBGs). The galaxies are also selected to be bright in rest-frame near-IR wavelengths and hence, massive. We identify a total of 16 BBGs. Fitting the spectral energy distribution (SED) of the BBGs show that the candidate galaxies have average estimated ages of ~800 Myr and average stellar masses of ~5x10^10 M_sun, consistent with being old and massive systems. Two of our BBG candidates are also identified by the criteria that is sensitive to star forming galaxies (LBG selection). We find a number density of ~3.2x10^-5 Mpc^-3 for the BBGs corresponding to a mass density of ~2.0x10^6 M_sun/Mpc^3 in the redshift range covering the survey. Given the old age and the passive evolution, it is argued that some of these objects formed the bulk of their mass only a few hundred million years after the Big Bang.

preprint2014arXiv

CANDELS/GOODS-S, CDFS, ECDFS: Photometric Redshifts For Normal and for X-Ray-Detected Galaxies

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.

preprint2014arXiv

Cosmic Web and Star Formation Activity in Galaxies at z~1

We investigate the role of the delineated cosmic web/filaments on the star formation activity by exploring a sample of 425 narrow-band selected Hα emitters, as well as 2846 color-color selected underlying star-forming galaxies for a large scale structure (LSS) at z=0.84 in the COSMOS field from the HiZELS survey. Using the scale-independent Multi-scale Morphology Filter (MMF) algorithm, we are able to quantitatively describe the density field and disentangle it into its major components: fields, filaments and clusters. We show that the observed median star formation rate (SFR), stellar mass, specific star formation rate (sSFR), the mean SFR-Mass relation and its scatter for both Hα emitters and underlying star-forming galaxies do not strongly depend on different classes of environment, in agreement with previous studies. However, the fraction of Hα emitters varies with environment and is enhanced in filamentary structures at z~1. We propose mild galaxy-galaxy interactions as the possible physical agent for the elevation of the fraction of Hα star-forming galaxies in filaments. Our results show that filaments are the likely physical environments which are often classed as the "intermediate" densities, and that the cosmic web likely plays a major role in galaxy formation and evolution which has so far been poorly investigated.

preprint2014arXiv

Kpc-scale Properties of Emission-line Galaxies

We perform a detailed study of the resolved properties of emission-line galaxies at kpc-scale to investigate how small-scale and global properties of galaxies are related. 119 galaxies with high-resolution Keck/DEIMOS spectra are selected to cover a wide range in morphologies over the redshift range 0.2<z<1.3. Using the HST/ACS and HST/WFC3 imaging data taken as a part of the CANDELS project, for each galaxy we perform SED fitting per resolution element, producing resolved rest-frame U-V color, stellar mass, star formation rate, age and extinction maps. We develop a technique to identify blue and red "regions" within individual galaxies, using their rest-frame color maps. As expected, for any given galaxy, the red regions are found to have higher stellar mass surface densities and older ages compared to the blue regions. Furthermore, we quantify the spatial distribution of red and blue regions with respect to both redshift and stellar mass, finding that the stronger concentration of red regions toward the centers of galaxies is not a significant function of either redshift or stellar mass. We find that the "main sequence" of star forming galaxies exists among both red and blue regions inside galaxies, with the median of blue regions forming a tighter relation with a slope of 1.1+/-0.1 and a scatter of ~0.2 dex compared to red regions with a slope of 1.3+/-0.1 and a scatter of ~0.6 dex. The blue regions show higher specific Star Formation Rates (sSFR) than their red counterparts with the sSFR decreasing since z~1, driver primarily by the stellar mass surface densities rather than the SFRs at a giver resolution element.

preprint2014arXiv

Physical Properties of Emission-Line Galaxies at z ~ 2 from Near-Infrared Spectroscopy with Magellan FIRE

We present results from near-infrared spectroscopy of 26 emission-line galaxies at z ~ 2 obtained with the FIRE spectrometer on the Magellan Baade telescope. The sample was selected from the WISP survey, which uses the near-infrared grism of the Hubble Space Telescope Wide Field Camera 3 to detect emission-line galaxies over 0.3 < z < 2.3. Our FIRE follow-up spectroscopy (R~5000) over 1.0-2.5 micron permits detailed measurements of physical properties of the z~2 emission-line galaxies. Dust-corrected star formation rates for the sample range from ~5-100 M_sun yr-1. We derive a median metallicity for the sample of ~0.45 Z_sun, and the estimated stellar masses range from ~10^8.5 - 10^9.5 M_sun. The average ionization parameters measured for the sample are typically much higher than what is found for local star-forming galaxies. We derive composite spectra from the FIRE sample, from which we infer typical nebular electron densities of ~100-400 cm^-3. Based on the location of the galaxies and composite spectra on BPT diagrams, we do not find evidence for significant AGN activity in the sample. Most of the galaxies as well as the composites are offset in the BPT diagram toward higher [O III]/H-beta at a given [N II]/H-alpha, in agreement with other observations of z > 1 star-forming galaxies, but composite spectra derived from the sample do not show an appreciable offset from the local star-forming sequence on the [O III]/H-beta versus [S II]/H-alpha diagram. We infer a high nitrogen-to-oxygen abundance ratio from the composite spectrum, which may contribute to the offset of the high-redshift galaxies from the local star-forming sequence in the [O III]/H-beta versus [N II]/H-alpha diagram. We speculate that the elevated nitrogen abundance could result from substantial numbers of Wolf-Rayet stars in starbursting galaxies at z~2. (Abridged)

preprint2014arXiv

Properties of Submillimeter Galaxies in the CANDELS GOODS-S Field

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.

preprint2014arXiv

Rapid Decline of Lyman-alpha Emission Toward the Reionization Era

The observed deficit of strongly Lyman-alpha emitting galaxies at z>6.5 is attributed to either increasing neutral hydrogen in the intergalactic medium (IGM) and/or to the evolving galaxy properties. To investigate this, we have performed very deep near-IR spectroscopy of z>7 galaxies using MOSFIRE on the Keck-I Telescope. We measure the Lyman-alpha fraction at z~8 (combined photometric redshift peak at z=7.7) using two methods. First, we derived NLyα/Ntot directly using extensive simulations to correct for incompleteness. Second, we used a Bayesian formalism (introduced by Treu et al. 2012) that compares the z>7 galaxy spectra to models of the Lyman-alpha equivalent width (WLyα) distribution at z~6. We explored two simple evolutionary scenarios: smooth evolution where Lyman-alpha is attenuated in all galaxies by a constant factor (perhaps owing to processes from galaxy evolution or a slowly increasing IGM opacity), and patchy evolution where Lyman-alpha is blocked in some fraction of galaxies (perhaps due to the IGM being opaque along only some fraction of sightlines). The Bayesian formalism places stronger constraints compared with the direct method. Combining our data with that in the literature we find that at z~8 the Lyman-alpha fraction has dropped by a factor >3(84% confidence interval) using both the smooth and patchy scenarios compared to the z~6 values. Furthermore, we find a tentative evidence that the data favor the patchy scenario over smooth (with "positive" Bayesian evidence), extending trends observed at z~7 to higher redshift. If this decrease is a result of reionization as predicted by theory, then our data imply the volume averaged neutral hydrogen fraction in the IGM to be >0.3 suggesting that the reionization of the universe is in progress at z~8.

preprint2014arXiv

Sub-Megaparsec Individual Photometric Redshift Estimation from Cosmic Web Constraints

We present a method, PhotoWeb, for estimating photometric redshifts of individual galaxies, and their equivalent distance, with megaparsec and even sub-megaparsec accuracy using the Cosmic Web as a constraint over photo-z estimates. PhotoWeb redshift errors for individual galaxies are of the order of delta_z = 0.0007, compared to errors of delta_z = 0.02 for current photo-z techniques. The mean redshift error is of the order of 0.00005-0.0004 compared to mean errors in the range delta_z =z 0.001-0.01 for the best available photo-z estimates in the literature. Current photo-z techniques based on the spectral energy distribution of galaxies and projected clustering produce redshift estimates with large errors due to the poor constraining power the galaxy's spectral energy distribution and projected clustering can provide. The Cosmic Web, on the other hand, provides the strongest constraints on the position of galaxies. The network of walls, filaments and voids occupy ~%10 of the volume of the Universe, yet they contain ~%95 of galaxies. The cosmic web, being a cellular system with well-defined boundaries, sets a restricted set of intermittent positions a galaxy can occupy along a given line-of-sight. Using the information in the density field computed from spectroscopic redshifts we can narrow the possible locations of a given galaxy along the line of sight from a single broad probability distribution (from photo-z) to one or a few narrow peaks. Our first results improve previous photo-z errors by more than one order of magnitude allowing sub-megaparsec errors in some cases. Such accurate estimates for tens of millions of galaxies will allow unprecedented galaxy-LSS studies. In this work we apply our technique to the SDSS photo-z galaxy sample and discuss its performance and future improvements.

preprint2014arXiv

The HST/ACS Coma Cluster Survey - VII. Structure and Assembly of Massive Galaxies in the Center of the Coma Cluster

We constrain the assembly history of galaxies in the projected central 0.5 Mpc of the Coma cluster by performing structural decomposition on 69 massive (M_star >= 10^9 M_sun) galaxies using high-resolution F814W images from the HST Treasury Survey of Coma. Each galaxy is modeled with up to three Sersic components having a free Sersic index n. After excluding the two cDs in the projected central 0.5 Mpc of Coma, 57% of the galactic stellar mass in the projected central 0.5 Mpc of Coma resides in classical bulges/ellipticals while 43% resides in cold disk-dominated structures. Most of the stellar mass in Coma may have been assembled through major (and possibly minor) mergers. Hubble types are assigned based on the decompositions, and we find a strong morphology-density relation; the ratio of (E+S0):spirals is (91.0%):9.0%. In agreement with earlier work, the size of outer disks in Coma S0s/spirals is smaller compared with lower-density environments captured with SDSS (Data Release 2). Among similar-mass clusters from a hierarchical semi-analytic model, no single cluster can simultaneously match all the global properties of the Coma cluster. The model strongly overpredicts the mass of cold gas and underpredicts the mean fraction of stellar mass locked in hot components over a wide range of galaxy masses. We suggest that these disagreements with the model result from missing cluster physics (e.g., ram-pressure stripping), and certain bulge assembly modes (e.g., mergers of clumps). Overall, our study of Coma underscores that galaxy evolution is not solely a function of stellar mass, but also of environment.

preprint2014arXiv

The interstellar medium and feedback in the progenitors of the compact passive galaxies at z~2

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.

preprint2014arXiv

The MOSDEF Survey: Mass, Metallicity, and Star-formation Rate at z~2.3

We present results on the z~2.3 mass-metallicity relation (MZR) using early observations from the MOSFIRE Deep Evolution Field (MOSDEF) survey. We use an initial sample of 87 star-forming galaxies with spectroscopic coverage of Hβ, [OIII]λ5007, Hα, and [NII]λ6584 rest-frame optical emission lines, and estimate the gas-phase oxygen abundance based on the N2 and O3N2 strong-line indicators. We find a positive correlation between stellar mass and metallicity among individual z~2.3 galaxies using both the N2 and O3N2 indicators. We also measure the emission-line ratios and corresponding oxygen abundances for composite spectra in bins of stellar mass. Among composite spectra, we find a monotonic increase in metallicity with increasing stellar mass, offset ~0.15-0.3 dex below the local MZR. When the sample is divided at the median star-formation rate (SFR), we do not observe significant SFR dependence of the z~2.3 MZR among either individual galaxies or composite spectra. We furthermore find that z~2.3 galaxies have metallicities ~0.1 dex lower at a given stellar mass and SFR than is observed locally. This offset suggests that high-redshift galaxies do not fall on the local "fundamental metallicity relation" among stellar mass, metallicity, and SFR, and may provide evidence of a phase of galaxy growth in which the gas reservoir is built up due to inflow rates that are higher than star-formation and outflow rates. However, robust conclusions regarding the gas-phase oxygen abundances of high-redshift galaxies await a systematic reappraisal of the application of locally calibrated metallicity indicators at high redshift.

preprint2014arXiv

Type Ia Supernova Rate Measurements to Redshift 2.5 from CANDELS : Searching for Prompt Explosions in the Early Universe

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.

preprint2013arXiv

A Critical Assessment of Photometric Redshift Methods: A CANDELS Investigation

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.

preprint2013arXiv

CANDELS Multi-wavelength Catalogs: Source Detection and Photometry in the GOODS-South Field

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$σ$ 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 $μ$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.

preprint2013arXiv

Hot-Dust (690K) Luminosity Density and its Evolution in the last 7.5Gyr

[Abridged] We study the contribution of hot-dust to the luminosity density of galaxies and its evolution with cosmic time. Using the Spitzer-IRAC data in the COSMOS field, we estimate the contribution from hot-dust at rest-frame 4.2um (from ~0 < z < ~0.2 up to ~0.5 < z < ~0.9). This wavelength corresponds to black-body temperature of ~690K. The contribution due to stellar emission is estimated from the rest-frame 1.6um luminosity (assumed to result from stellar emission alone) and subtracted from the mid-infrared luminosity of galaxies to measure hot-dust emission. In order to attempt the study of the 3.3um-PAH feature, we use the rest-frame 4.2um to infer the hot-dust flux at 3.3um. This study is performed for different spectral types of galaxies: early-type, late-type, starburst, and IR-selected AGN. We find that: (a) the decrease of the hot-dust luminosity density since ~0.5 < z < ~1 is steeper (by at least ~0.5dex) compared to that of the cold-dust, giving support to the scenario where galaxy obscuration increases with redshift, as already proposed in the literature; (b) hot-dust and PAH emission evolution seems to be correlated with stellar mass, where rest-frame 1.6um luminous non-AGN galaxies (i.e., massive systems) show a stronger decrement (with decreasing redshift) in hot-dust and PAH emission than the less luminous (less massive) non-AGN galaxies; (c) despite comprising < ~3% of the total sample, AGN contribute as much as a third to the hot-dust luminosity density at z < 1 and clearly dominate the bright-end of the total hot-dust Luminosity Density Function at ~0.5 < z < ~0.9; (d) the average dust-to-total luminosity ratio increases with redshift, while PAH-to-total luminosity ratio remains fairly constant; (e) at M_1.6 > -25, the dust-to-total and PAH-to-total luminosity ratios increase with decreasing luminosity, but deeper data is required to confirm this result.

preprint2013arXiv

The dependency of AGN infrared color-selection on source luminosity and obscuration: An observational perspective in CDFS and COSMOS

AIMS. This work addresses the AGN IR-selection dependency on intrinsic source luminosity and obscuration, in order to identify and charaterise biases which could affect conclusions in studies. METHODS. We study IR-selected AGN in the Chandra Deep Field South (CDFS) survey and in the Cosmological Survey (COSMOS). The AGN sample is divided into low and high X-ray luminosity classes, and into unobscured (type-1) and obscured (type-2) classes by means of X-ray and optical spectroscopy data. Specifically in the X-rays regime, we adopt the intrinsic luminosity taking into account the estimated column density (N_H). We also take the chance to highlight important differences resulting from adopting different methods to assess AGN obscuration. RESULTS. In agreement with previous studies, we also find that AGN IR-selection efficiency shows a decrease with decreasing source AGN X-ray luminosity. For the intermediate-luminosity AGN population (43.3<~log(L_X[erg/s])<~44$), the efficiency also worsens with increasing obscuration. The same sample also shows an evolution with cosmic time of the obscured fraction at the highest X-ray luminosities, independently of adopted type-1/type-2 classification method. CONCLUSIONS. We confirm that AGN IR-selection is genuinely biased toward unobscured AGN, but only at intermediate luminosities. At the highest luminosities, where AGN IR-selection is more efficient, there is no obscuration bias. We show that type-1 AGN are intrinsically more luminous than type-2 AGN only at z<~1.6, thus resulting in more type-1 AGN being selected the shallower the IR survey is. [ABRIDGED: Conclusion on survey strategy. See manuscript.] Finally, when the James Webb Space Telescope comes online, the broad-band filters 2.0um, 4.4um, 7.7um, and 18um will be key to disentangle AGN from non-AGN dominated SEDs at depths where spectroscopy becomes impractical.

preprint2013arXiv

The Discovery of the Most Distant Known Type Ia Supernova at Redshift 1.914

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 $Λ$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.

preprint2013arXiv

The stellar mass function of star-forming galaxies and the mass-dependent SFR function since z=2.23 from HiZELS

We explore a large uniformly selected sample of H$α$ selected star-forming galaxies (SFGs) at z=0.40,0.84,1.47,2.23 to unveil the evolution of the star formation rate (SFR) function and the stellar mass function. We find strong evolution in the SFR function, with the typical SFR of SFGs declining exponentially with time in the last 11Gyrs as SFR$^*$(T[Gyr])=10$^{4.23/T+0.37}$ M$_{\odot}$yr$^{-1}$, but with no evolution in the faint-end slope, $α\approx$-1.6. The stellar mass function of SFGs, however, reveals little evolution: $α \approx$-1.4, M$^*$$\approx$10$^{11.2\pm0.2}$ M$_{\odot}$ and just a slight increase of $\approx$2.3x in $Φ^*$ from z=2.23 to z=0.4. The stellar mass density within SFGs has been roughly constant since z=2.23 at $\approx$10$^{7.65\pm0.08}$ M$_{\odot}$Mpc$^{-3}$, comprising $\approx$100% of the stellar mass density in all galaxies at z=2.23, and declining to $\approx$20% by z=0.4, driven by the rise of the passive population. We find that SFGs with M$\approx$10$^{10.0\pm0.2}$ M$_{\odot}$ contribute most to the SFR density ($ρ_{\rm SFR}$) per dlog$_{10}$M, and that there is no significant evolution in the fractional contribution from SFGs of different masses to $ρ_{\rm SFR}$ or $ρ_{\rm SFR}$(dlog$_{10}$M)$^{-1}$ since z=2.23. Instead, we show that the decline of SFR$^*$ and of $ρ_{\rm SFR}$ are primarily driven by an exponential decline in SFRs at all masses. Our results have important implications not only on how SFGs need to be quenched across cosmic time, but also on the driver(s) of the exponential decline in SFR$^*$ from $\approx$66 M$_{\odot}$yr$^{-1}$ to $\approx$5 M$_{\odot}$yr$^{-1}$ since z=2.23.

preprint2012arXiv

A new infrared color criterion for the selection of 0<z<7 AGN: application to deep fields and implications for JWST surveys

[Abridged] It is widely accepted that observations at mid-infrared (mid-IR) wavelengths enable the selection of galaxies with nuclear activity, which may not be revealed even in the deepest X-ray surveys. In this work new near- and mid-IR color diagnostics are explored, aiming for improved efficiency - better completeness and less contamination - in selecting AGN out to very high redshifts. We restrict our study to the James Webb Space Telescope wavelength range (0.6-27um). The criteria are created based on the predictions by state-of-the-art galaxy and AGN templates covering a wide variety of galaxy properties, and tested against control samples with deep multi-wavelength coverage (ranging from the X-rays to radio frequencies). We show that the colors Ks-[4.5], [4.5]-[8.0], and [8.0]-[24] are ideal as AGN/non-AGN diagnostics at, respectively, z<~1, 1<~z<~2.5, and z>~2.5-3. However, when the source redshift is unknown, these colors should be combined. We thus develop an improved IR criterion (using Ks and IRAC bands, KI) as a new alternative at z<~2.5. KI does not show improved completeness (50-60% overall) in comparison to commonly used IRAC-based AGN criteria, but is less affected by non-AGN contamination (revealing a >50-90% level of successful AGN selection). We also propose KIM (using Ks, IRAC, and MIPS-24um bands, KIM), which aims to select AGN hosts from local distances to as far back as the end of reionization (0<z<~7) with reduced non-AGN contamination. However, the necessary testing-constraints and the small control-sample sizes prevent the confirmation of its improved efficiency at z<~2.5. Overall, KIM shows a ~30-40% completeness and a >70-90% level of successful AGN selection. KI and KIM are built to be reliable against a ~10-20% error in flux, are based on existing filters, and are suitable for immediate use.

preprint2012arXiv

A Type Ia Supernova at Redshift 1.55 in Hubble Space Telescope Infrared Observations from CANDELS

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.

preprint2012arXiv

CANDELS: The Evolution of Galaxy Rest-Frame Ultraviolet Colors from z = 8 to 4

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.

preprint2012arXiv

Evolution of the Quasar Luminosity Function Over 3 < z < 5 in the COSMOS Survey Field

We investigate the high-redshift quasar luminosity function (QLF) down to an apparent magnitude of I(AB) = 25 in the Cosmic Evolution Survey (COSMOS). Careful analysis of the extensive COSMOS photometry and imaging data allows us to identify and remove stellar and low-redshift contaminants, enabling a selection that is nearly complete for type-1 quasars at the redshifts of interest. We find 155 likely quasars at z > 3.1, 39 of which have prior spectroscopic confirmation. We present our sample in detail and use these confirmed and likely quasars to compute the rest-frame UV QLF in the redshift bins 3.1 < z < 3.5 and 3.5 < z < 5. The space density of faint quasars decreases by roughly a factor of four from z \sim 3.2 to z \sim 4, with faint-end slopes of β \sim -1.7 at both redshifts. The decline in space density of faint optical quasars at z > 3 is similar to what has been found for more luminous optical and X-ray quasars. We compare the rest-frame UV luminosity functions found here with the X-ray luminosity function at z > 3, and find that they evolve similarly between z \sim 3.2 and z \sim 4; however, the different normalizations imply that roughly 75% of X-ray bright active galactic nuclei (AGN) at z \sim 3 - 4 are optically obscured. This fraction is higher than found at lower redshift and may imply that the obscured, type-2 fraction continues to increase with redshift at least to z \sim 4. Finally, the implications of the results derived here for the contribution of quasars to cosmic reionization are discussed.

preprint2012arXiv

Lyα Emission from High Redshift Sources in COSMOS

We investigate spectroscopically measured Lyα equivalent widths and escape fractions of 244 sources of which 95 are Lyman Break Galaxies (LBGs) and 106 Lyman Alpha Emitters (LAEs) at z~4.2, z~4.8, and z~5.6 selected from intermediate and narrow band observations. The sources were selected from the Cosmic Evolution Survey (COSMOS), and observed with the DEIMOS spectrograph. We find that the distribution of equivalent widths shows no evolution with redshift for both the LBG selected sources and the intermediate/narrowband LAEs. We also find that the Lyα escape fraction of intermediate and narrow band LAEs is on average higher and has a larger variation than the escape fraction of LBG selected sources. The escape fraction does not show a dependence with redshift. Similar to what has been found for LAEs at low redshifts, the sources with the highest extinctions show the lowest escape fractions. The range of escape fractions increases with decreasing extinction. This is evidence that the dust extinction is the most important factor affecting the escape of Lyα photons, but at low extinctions other factors such as HI covering fraction and gas kinematics can be just as effective at inhibiting the escape of Lyα photons.

preprint2012arXiv

Near-Infrared Survey of the GOODS-North Field: Search for Luminous Galaxy Candidates at z=>6.5

We present near-infrared (NIR; J & Ks) survey of the Great Observatories Origins Deep Survey-North (GOODS-N) field. The publicly available imaging data were obtained using the MOIRCS instrument on the 8.2m Subaru and the WIRCam instrument on the 3.6m Canada-France-Hawaii Telescope (CFHT). These observations fulfill a serious wavelength gap in the GOODS-N data - i.e., lack of deep NIR observations. We combine the Subaru/MOIRCS and CFHT/WIRCam archival data to generate deep J and Ks band images, covering the full GOODS-N field (~169 sq. arcmin) to an AB magnitude limit of ~25 mag (3sigma). We applied z'-band dropout color selection criteria, using the NIR data generated here. We have identified two possible Lyman Break Galaxy (LBG) candidates at z\gtrsim6.5 with J\lesssim24.5. The first candidate is a likely LBG at z\sim6.5 based on a weak spectral feature tentatively identified as Lyalpha line in the deep Keck/DEIMOS spectrum, while the second candidate is a possible LBG at z\sim7 based on its photometric redshift. These z'-dropout objects, if confirmed, are among the brightest such candidates found so far. At z\gtrsim6.5, their star formation rate is estimated as 100-200 solar mass per year. If they continue to form stars at this rate, they assemble a stellar mass of ~5x10^10 solar mass after about 400 million years, becoming the progenitors of massive galaxies observed at z\sim5. We study the implication of the z'-band dropout candidates discovered here, in constraining the bright-end of the luminosity function and understanding the nature of high redshift galaxies.

preprint2012arXiv

SDSS 0956+5128: A Broad-line Quasar with Extreme Velocity Offsets

We report on the discovery of a Type 1 quasar, SDSS 0956+5128, with a surprising combination of extreme velocity offsets. SDSS 0956+5128 is a broad-lined quasar exhibiting emission lines at three substantially different redshifts: a systemic redshift of z ~ 0.714 based on narrow emission lines, a broad MgII emission line centered 1200 km/s bluer than the systemic velocity, at z ~ 0.707, and broad Hαand Hβemission lines centered at z ~ 0.690. The Balmer line peaks are 4100 km/s bluer than the systemic redshift. There are no previously known objects with such an extreme difference between broad MgII and broad Balmer emission. The two most promising explanations are either an extreme disk emitter or a high-velocity black hole recoil. However, neither explanation appears able to explain all of the observed features of SDSS 0956+5128, so the object may provide a challenge to our general understanding of quasar physics.

preprint2012arXiv

The Dwarf Starburst Host Galaxy of a Type Ia SN at z = 1.55 from CANDELS

We present VLT/X-shooter observations of a high redshift, type Ia supernova host galaxy, discovered with HST/WFC3 as part of the CANDELS Supernova project. The galaxy exhibits strong emission lines of Lyα, [O II], Hβ, [O III], and Hα at z = 1.54992(+0.00008-0.00004). From the emission-line fluxes and SED fitting of broad-band photometry we rule out AGN activity and characterize the host galaxy as a young, low mass, metal poor, starburst galaxy with low intrinsic extinction and high Lyα escape fraction. The host galaxy stands out in terms of the star formation, stellar mass, and metallicity compared to its lower redshift counterparts, mainly because of its high specific star-formation rate. If valid for a larger sample of high-redshift SN Ia host galaxies, such changes in the host galaxy properties with redshift are of interest because of the potential impact on the use of SN Ia as standard candles in cosmology.

preprint2012arXiv

The Extended Hubble Space Telescope Supernova Survey: The Rate of Core Collapse Supernovae to z~1

We use a sample of 45 core collapse supernovae detected with the Advanced Camera for Surveys on-board the Hubble Space Telescope to derive the core collapse supernova rate in the redshift range 0.1<z<1.3. In redshift bins centered on <z>=0.39, <z>=0.73, and <z>=1.11, we find rates 3.00 {+1.28}{-0.94}{+1.04}{-0.57}, 7.39 {+1.86}{-1.52}{+3.20}{-1.60}, and 9.57 {+3.76}{-2.80}{+4.96}{-2.80}, respectively, given in units yr^{-1} Mpc^{-3} 10^{-4} h_{70}^3. The rates have been corrected for host galaxy extinction, including supernovae missed in highly dust enshrouded environments in infrared bright galaxies. The first errors represent statistical while the second are the estimated systematic errors. We perform a detailed discussion of possible sources of systematic errors and note that these start to dominate over statistical errors at z>0.5, emphasizing the need to better control the systematic effects. For example, a better understanding of the amount of dust extinction in the host galaxies and knowledge of the supernova luminosity function, in particular the fraction of faint M > -15 supernovae, is needed to better constrain the rates. When comparing our results with the core collapse supernova rate based on the star formation rate, we find a good agreement, consistent with the supernova rate following the star formation rate, as expected.

preprint2012arXiv

The HST/ACS Coma Cluster Survey. VIII. Barred Disk Galaxies in the Core of the Coma Cluster

(ABRIDGED) We use high resolution (~0.1") F814W ACS images from the HST ACS Treasury survey of the Coma cluster at z~0.02 to study bars in massive disk galaxies (S0s), and in dwarf galaxies in the Coma core. Our study helps constrain the evolution of bars and disks in dense environments and provides a comparison point for studies in lower density environments and at higher redshifts. (1) We characterize the fraction and properties of bars in a sample of 32 bright (M_V <= -18, M_* > 10^9.5 M_sun) S0 galaxies, which dominate the population of massive disk galaxies in the Coma core. Measuring the S0 bar fraction must be handled carefully, as the results depend on the method used: the bar fraction for bright S0s in the Coma core is 50%+/-11%, 65%+/-11%, and 60%+/-11% for three methods of bar detection: strict ellipse fit criteria, relaxed ellipse fit criteria, and visual classification. (2) We compare the S0 bar fraction across different environments (Coma core, A901/902, Virgo). We find that the bar fraction among bright S0 galaxies does not show a statistically significant variation across environments spanning two orders of magnitude in galaxy number density (n~300-10,000 gal/Mpc^3). We speculate that the S0 bar fraction is not significantly enhanced in rich clusters because S0s in rich clusters are less prone to bar instabilities as they are dynamically hot and gas poor due to ram pressure stripping and accelerated star formation. In addition, high-speed encounters in rich clusters may be less effective than slow, strong encounters in inducing bars. (3) We analyze a sample of 333 faint (M_V > -18) dwarf galaxies in the Coma core. Using unsharp-masking, we find only 13 galaxies with bar and/or spiral structure. The paucity of disk structures in Coma dwarfs suggests that either disks are not common in these galaxies, or that any disks present are too hot to develop instabilities.

preprint2011arXiv

A massive proto-cluster of galaxies at a redshift of z {\approx} 5.3

Massive clusters of galaxies have been found as early as 3.9 Billion years (z=1.62) after the Big Bang containing stars that formed at even earlier epochs. Cosmological simulations using the current cold dark matter paradigm predict these systems should descend from "proto-clusters" - early over-densities of massive galaxies that merge hierarchically to form a cluster. These proto-cluster regions themselves are built-up hierarchically and so are expected to contain extremely massive galaxies which can be observed as luminous quasars and starbursts. However, observational evidence for this scenario is sparse due to the fact that high-redshift proto-clusters are rare and difficult to observe. Here we report a proto-cluster region 1 billion years (z=5.3) after the Big Bang. This cluster of massive galaxies extends over >13 Mega-parsecs, contains a luminous quasar as well as a system rich in molecular gas. These massive galaxies place a lower limit of >4x10^11 solar masses of dark and luminous matter in this region consistent with that expected from cosmological simulations for the earliest galaxy clusters.

preprint2011arXiv

CANDELS: Constraining the AGN-Merger Connection with Host Morphologies at z~2

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.

preprint2011arXiv

CANDELS: The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey

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σ 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σ 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.

preprint2011arXiv

CANDELS: The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey - The Hubble Space Telescope Observations, Imaging Data Products and Mosaics

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.

preprint2011arXiv

Spectroscopic Confirmation of Three z-Dropout Galaxies at z = 6.844 - 7.213: Demographics of Lyman-Alpha Emission in z ~ 7 Galaxies

We present the results of our ultra-deep Keck/DEIMOS spectroscopy of z-dropout galaxies in the SDF and GOODS-N. For 3 out of 11 objects, we detect an emission line at ~ 1um with a signal-to-noise ratio of ~ 10. The lines show asymmetric profiles with high weighted skewness values, consistent with being Lya, yielding redshifts of z=7.213, 6.965, and 6.844. Specifically, we confirm the z=7.213 object in two independent DEIMOS runs with different spectroscopic configurations. The z=6.965 object is a known Lya emitter, IOK-1, for which our improved spectrum at a higher resolution yields a robust skewness measurement. The three z-dropouts have Lya fluxes of 3 x 10^-17 erg s^-1 cm^-2 and rest-frame equivalent widths EW_0^Lya = 33-43A. Based on the largest spectroscopic sample of 43 z-dropouts that is the combination of our and previous data, we find that the fraction of Lya-emitting galaxies (EW_0^Lya > 25A) is low at z ~ 7; 17 +- 10% and 24 +- 12% for bright (Muv ~= -21) and faint (Muv ~= -19.5) galaxies, respectively. The fractions of Lya-emitting galaxies drop from z ~ 6 to 7 and the amplitude of the drop is larger for faint galaxies than for bright galaxies. These two pieces of evidence would indicate that the neutral hydrogen fraction of the IGM increases from z ~ 6 to 7, and that the reionization proceeds from high- to low-density environments, as suggested by an inside-out reionization model.

preprint2011arXiv

The UDF05 Follow-up of the Hubble Ultra Deep Field. III. The Luminosity Function at z~6

In this paper, we present a derivation of the rest-frame 1400A luminosity function (LF) at redshift six from a new application of the maximum likelihood method by exploring the five deepest HST/ACS fields, i.e., the HUDF, two UDF05 fields, and two GOODS fields. We work on the latest improved data products, which makes our results more robust than those of previous studies. We use un-binned data and thereby make optimal use of the information contained in the dataset. We focus on the analysis to a magnitude limit where the completeness is larger than 50% to avoid possibly large errors in the faint end slope that are difficult to quantify. We also take into account scattering in and out of the dropout sample due to photometric errors by defining for each object a probability that it belongs to the dropout sample. We find the best fit Schechter parameters to the z~6 LF are: alpha = 1.87 +/- 0.14, M* = -20.25 +/- 0.23, and phi*=1.77^{+0.62}_{-0.49} * 10^{-3} Mpc^{-3}. Such a steep slope suggests that galaxies, especially the faint ones, are possibly the main sources of ionizing photons in the universe at redshift six. We also combine results from all studies at z~6 to reach an agreement in 95% confidence level that -20.45<M*<-20.05 and -1.90<alpha<-1.55. The luminosity density has been found not to evolve significantly between z~6 and z~5, but considerable evolution is detected from z~6 to z~3.

preprint2010arXiv

A Detailed Study of Photometric Redshifts for GOODS-South Galaxies

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.

preprint2010arXiv

A dozen new galaxies caught in the act: Gas stripping and extended emission line regions in the Coma cluster

We present images of extended H-alpha clouds associated with 14 member galaxies in the Coma cluster obtained from deep narrow band imaging observations with Suprime-Cam at the Subaru Telescope. The parent galaxies of the extended H-alpha clouds are distributed farther than 0.2 Mpc from the peak of X-ray emission of the cluster. Most of the galaxies have colors bluer than g-r approx 0.5 and they account for 57% of the blue (g-r<0.5) bright (r<17.8 mag) galaxies in the central region of the Coma cluster. They reside near the red- and blue-shifted edges of the Coma cluster's radial velocity distribution. These findings suggest that the most of the parent galaxies were recently captured by the Coma cluster potential and are now infalling toward the cluster center with their disk gas being stripped off and producing the observed H-alpha clouds.

preprint2010arXiv

A multi-wavelength approach to the properties of Extremely Red Galaxy populations: I - Contribution to the Star Formation Rate density and AGN content

We present a multi-wavelength analysis of the properties of Extremely Red Galaxy (ERG) populations, selected in the GOODS-South/Chandra Deep Field South field. By using all the photometric and spectroscopic information available on large deep samples of EROs (645 sources), IEROs (294 sources), and DRGs (350 sources), we derive redshift distributions, identify AGN powered and Star-formation powered galaxies, and, using the radio observations of this field, estimate robust (AGN- and dust-unbiased) Star Formation Rate Densities (SFRD) for these populations. We also investigate the properties of "pure" (galaxies that conform to only one of the three ERG criteria considered) and "combined" (galaxies that verify all three criteria) sub-populations. Overall, a large number of AGN are identified (up to ~30%, based on X-rays, and mid-infrared criteria), the majority of which are type-2 (obscured) objects. Among ERGs with no evidence for AGN activity, we identify sub-populations covering a wide range of average star-formation rates, from below 10 Mo/yr to as high as 200 Mo/yr. Applying a redshift separation (1<z<2 and 2<z<3) we find significant evolution (an increase of a factor of 2 or higher) of SFRD for EROs and DRGs, while none is observed for IEROs. The former populations can contribute more than 20\% to the global SFRD at 2<z<3. The emission from AGN activity is typically not strong in the ERG population, with AGN increasing the average radio luminosity of ERG sub-populations by, nominally, less than 20%. AGN are common, however, and, if no discrimination is attempted, this could significantly increase the SFRD estimate (by over 100% in some cases). Thus, and while the contribution of star forming processes to the radio luminosity in galaxies with AGN remains uncertain, a comprehensive identification of AGN in these populations is necessary to obtain meaningful results.

preprint2010arXiv

Deep GALEX Observations of the Coma Cluster: Source Catalog and Galaxy Counts

We present a source catalog from deep 26 ks GALEX observations of the Coma cluster in the far-UV (FUV; 1530 A) and near-UV (NUV; 2310 A) wavebands. The observed field is centered 0.9 deg (1.6 Mpc) south-west of the Coma core, and has full optical photometric coverage with SDSS. The catalog consists of 9700 galaxies with GALEX and SDSS photometry, including 242 spectroscopically-confirmed Coma member galaxies that range from giant spirals and elliptical galaxies to dwarf irregular and early-type galaxies. The full multi-wavelength catalog (cluster plus background galaxies) is ~80% complete to NUV=23 and FUV=23.5, and has a limiting depth at NUV=24.5 and FUV=25.0 which corresponds to a star formation rate of ~0.001 Msun/yr at the distance of Coma. Our deep GALEX observations required a two-fold approach to generating a source catalog: we used a Bayesian deblending algorithm to measure faint and compact sources (using SDSS coordinates as a position prior), and relied on the GALEX pipeline catalog for bright/extended objects. We performed simulations to assess the influence that systematic effects (e.g. object blends, source confusion, Eddington Bias) have on source detection and photometry when using both methods. The Bayesian deblending method roughly doubles the number of source detections and provides reliable photometry to a few magnitudes deeper than the GALEX pipeline catalog. This method is also free from source confusion over the UV magnitude range studied here; conversely, we estimate that the GALEX pipeline catalogs are confusion limited at magnitudes fainter than NUV~23 and FUV~24. We have measured the total UV galaxy counts using our catalog and report a ~50% excess of counts across FUV=22-23.5 and NUV=21.5-23 relative to previous GALEX measurements, which is not attributed to cluster member galaxies. Our galaxy counts are a better match to deeper UV counts measured with HST.

preprint2010arXiv

The Herschel-SPIRE Legacy Survey (HSLS): the scientific goals of a shallow and wide submillimeter imaging survey with SPIRE

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

preprint2010arXiv

The HST/ACS Coma Cluster Survey. II. Data Description and Source Catalogs

The Coma cluster was the target of a HST-ACS Treasury program designed for deep imaging in the F475W and F814W passbands. Although our survey was interrupted by the ACS instrument failure in 2007, the partially completed survey still covers ~50% of the core high-density region in Coma. Observations were performed for 25 fields that extend over a wide range of cluster-centric radii (~1.75 Mpc) with a total coverage area of 274 arcmin^2. The majority of the fields are located near the core region of Coma (19/25 pointings) with six additional fields in the south-west region of the cluster. In this paper we present reprocessed images and SExtractor source catalogs for our survey fields, including a detailed description of the methodology used for object detection and photometry, the subtraction of bright galaxies to measure faint underlying objects, and the use of simulations to assess the photometric accuracy and completeness of our catalogs. We also use simulations to perform aperture corrections for the SExtractor Kron magnitudes based only on the measured source flux and half-light radius. We have performed photometry for ~73,000 unique objects; one-half of our detections are brighter than the 10-sigma point-source detection limit at F814W=25.8 mag (AB). The slight majority of objects (60%) are unresolved or only marginally resolved by ACS. We estimate that Coma members are 5-10% of all source detections, which consist of a large population of unresolved objects (primarily GCs but also UCDs) and a wide variety of extended galaxies from a cD galaxy to dwarf LSB galaxies. The red sequence of Coma member galaxies has a constant slope and dispersion across 9 magnitudes (-21<M_F814W<-13). The initial data release for the HST-ACS Coma Treasury program was made available to the public in 2008 August. The images and catalogs described in this study relate to our second data release.

preprint2010arXiv

The Spitzer Survey of Stellar Structure in Galaxies (S^4G)

The Spitzer Survey of Stellar Structure in Galaxies S^4G is an Exploration Science Legacy Program approved for the Spitzer post-cryogenic mission. It is a volume-, magnitude-, and size-limited (d < 40 Mpc, |b| > 30 degrees, m_(Bcorr) < 15.5, D25>1') survey of 2,331 galaxies using IRAC at 3.6 and 4.5 microns. Each galaxy is observed for 240 s and mapped to > 1.5 x D25. The final mosaicked images have a typical 1 sigma rms noise level of 0.0072 and 0.0093 MJy / sr at 3.6 and 4.5 microns, respectively. Our azimuthally-averaged surface brightness profile typically traces isophotes at mu_3.6 (AB) (1 sigma) ~ 27 mag arcsec^-2, equivalent to a stellar mass surface density of ~ 1 Msun pc^-2. S^4G thus provides an unprecedented data set for the study of the distribution of mass and stellar structures in the local Universe. This paper introduces the survey, the data analysis pipeline and measurements for a first set of galaxies, observed in both the cryogenic and warm mission phase of Spitzer. For every galaxy we tabulate the galaxy diameter, position angle, axial ratio, inclination at mu_3.6 (AB) = 25.5 and 26.5 mag arcsec^-2 (equivalent to ~ mu_B (AB) =27.2 and 28.2 mag arcsec^-2, respectively). These measurements will form the initial S^4G catalog of galaxy properties. We also measure the total magnitude and the azimuthally-averaged radial profiles of ellipticity, position angle, surface brightness and color. Finally, we deconstruct each galaxy using GALFIT into its main constituent stellar components: the bulge/spheroid, disk, bar, and nuclear point source, where necessary. Together these data products will provide a comprehensive and definitive catalog of stellar structures, mass and properties of galaxies in the nearby Universe.

preprint2009arXiv

Large Area Survey for z=7 Galaxies in SDF and GOODS-N: Implications for Galaxy Formation and Cosmic Reionization

We present results of our large-area survey for z'-band dropout galaxies at z=7 in a 1568 arcmin^2 sky area covering the SDF and GOODS-N fields. Combining our ultra-deep Subaru/Suprime-Cam z'- and y-band (lambda_eff=1um) images with legacy data of Subaru and HST, we have identified 22 bright z-dropout galaxies down to y=26, one of which has a spectroscopic redshift of z=6.96 determined from Lya emission. The z=7 luminosity function (LF) yields the best-fit Schechter parameters of phi*=0.69 +2.62/-0.55 x10^(-3) Mpc^(-3), Muv*=-20.10 +/-0.76 mag, and alpha=-1.72 +/-0.65, and indicates a decrease from z=6 at a >95% confidence level. This decrease is beyond the cosmic variance in our two fields, which is estimated to be a factor of <~2. We have found that the cosmic star formation rate density drops from the peak at z=2-3 to z=7 roughly by a factor of ~10 but not larger than ~100. A comparison with the reionization models suggests either that the Universe could not be totally ionized by only galaxies at z=7, or more likely that properties of galaxies at z=7 are different from those at low redshifts having, e.g., a larger escape fraction (>~0.2), and/or a flatter IMF. Our SDF z-dropout galaxies appear to form 60-Mpc long filamentary structures, and the z=6.96 galaxy with Lya emission is located at the center of an overdense region consisting of four UV bright dropout candidates, which might suggest an existence of a well-developed ionized bubble at z=7.

preprint2009arXiv

Optical and near infra-red colours as a discriminant of the age and metallicity of stellar populations

We present a comprehensive analysis of the ability of current stellar population models to reproduce the optical (ugriz) and near infra-red (JHK) colours of a small sample of well-studied nearby elliptical and S0 galaxies. We find broad agreement between the ages and metallicities derived using different population models, although different models show different systematic deviations from the measured broad-band fluxes. Although it is possible to constrain Simple Stellar Population models to a well defined area in age-metallicity space, there is a clear degeneracy between these parameters even with such a full range of precise colours. The precision to which age and metallicity can be determined independently, using only broad band photometry with realistic errors, is Delta{[Fe/H]} ~ 0.18 and Delta{log(Age)} ~ 0.25. To constrain the populations and therefore the star formation history further it will be necessary to combine broad-band optical-IR photometry with either spectral line indices, or else photometry at wavelengths outside of this range.

preprint2007arXiv

Dark matter maps reveal cosmic scaffolding

Ordinary baryonic particles (such as protons and neutrons) account for only one-sixth of the total matter in the Universe. The remainder is a mysterious "dark matter" component, which does not interact via electromagnetism and thus neither emits nor reflects light. As dark matter cannot be seen directly using traditional observations, very little is currently known about its properties. It does interact via gravity, and is most effectively probed through gravitational lensing: the deflection of light from distant galaxies by the gravitational attraction of foreground mass concentrations. This is a purely geometrical effect that is free of astrophysical assumptions and sensitive to all matter -- whether baryonic or dark. Here we show high fidelity maps of the large-scale distribution of dark matter, resolved in both angle and depth. We find a loose network of filaments, growing over time, which intersect in massive structures at the locations of clusters of galaxies. Our results are consistent with predictions of gravitationally induced structure formation, in which the initial, smooth distribution of dark matter collapses into filaments then into clusters, forming a gravitational scaffold into which gas can accumulate, and stars can be built.

preprint2007arXiv

Evolution of the Bar Fraction in COSMOS: Quantifying the Assembly of the Hubble Sequence

We have analyzed the redshift-dependent fraction of galactic bars over 0.2<z<0.84 in 2,157 luminous face-on spiral galaxies from the COSMOS 2-square degree field. Our sample is an order of magnitude larger than that used in any previous investigation, and is based on substantially deeper imaging data than that available from earlier wide-area studies of high-redshift galaxy morphology. We find that the fraction of barred spirals declines rapidly with redshift. Whereas in the local Universe about 65% of luminous spiral galaxies contain bars (SB+SAB), at z ~0.84 this fraction drops to about 20%. Over this redshift range the fraction of strong (SB) bars drops from about 30% to under 10%. It is clear that when the Universe was half its present age, the census of galaxies on the Hubble sequence was fundamentally different from that of the present day. A major clue to understanding this phenomenon has also emerged from our analysis, which shows that the bar fraction in spiral galaxies is a strong function of stellar mass, integrated color and bulge prominence. The bar fraction in very massive, luminous spirals is about constant out to z ~ 0.84 whereas for the low mass, blue spirals it declines significantly with redshift beyond z=0.3. There is also a slight preference for bars in bulge dominated systems at high redshifts which may be an important clue towards the co-evolution of bars, bulges and black holes. Our results thus have important ramifications for the processes responsible for galactic downsizing, suggesting that massive galaxies matured early in a dynamical sense, and not just as a result of the regulation of their star formation rate.

preprint2007arXiv

Redshifts of Emission Line Objects in the Hubble Ultra Deep Field

We present redshifts for 115 emission line objects in the Hubble Ultra Deep Field (HUDF) identified through the GRism ACS Program for Extragalactic Science (GRAPES) project using the slitless grism spectroscopy mode of the ACS Camera on the Hubble Space Telescope (HST). The sample was selected by an emission line search on all extracted 1-dimensional GRAPES spectra. We identify the emission lines using line wavelength ratios where multiple lines are detected in the grism wavelength range (5800A < lambda < 9600A), and using photometric redshift information where multiple lines are unavailable. We then derive redshifts using the identified lines. Our redshifts are accurate to delta(z) = 0.009, based on both statistical uncertainty estimates and comparison with published ground-based spectra. Over 40% of our sample is fainter than typical magnitude limits for ground-based spectroscopy (with i_{AB}>25 mag). Such emission lines would likely remain undiscovered without our deep survey. The emission line objects fall into 3 categories: 1) Most are low to moderate redshift galaxies (0 < z < 2), including many actively star forming galaxies with strong HII regions; 2) 9 are high redshift (4 < z < 7) Lyman-alpha emitters; and 3) at least 3 are candidate AGNs.

preprint2007arXiv

The Far-Infrared Luminosity Function from GOODS-N: Constraining the Evolution of Infrared Galaxies for z \leq 1

We present the IR luminosity function derived from ultra-deep 70 micron imaging of the GOODS-North field. The 70 micron observations are longward of the PAH and silicate features which complicate work in the MIR. We derive far-infrared luminosities for the 143 sources with S_{70} > 2 mJy (S/N > 3 σ). The majority (81%) of the sources have spectroscopic redshifts, and photometric redshifts are calculated for the remainder. The IR luminosity function at four redshifts (z ~ 0.28, 0.48, 0.78, and 0.97) is derived and compared to the local one. There is considerable degeneracy between luminosity and density evolution. If the evolving luminosity function is described as ρ(L, z) = (1 + z)^q ρ(L/(1 + z)^p, 0), we find q = -2.19p + 6.09. In the case of pure luminosity evolution, we find a best fit of p = 2.78^{+0.34}_{-0.32}. This is consistent with the results from 24 micron and 1.4 GHz studies. Our results confirm the emerging picture of strong evolution in LIRGs and ULIRGs at 0.4 < z < 1.1, but we find no evidence of significant evolution in the sub-LIRG (L < 10^{11} L_{\odot}) population for z < 0.4.

preprint2006arXiv

Infrared and Optical Observations of GRB 030115 and its Extremely Red Host Galaxy: Implications for Dark Bursts

We present near-infrared (nIR) and optical observations of the afterglow of GRB 030115. Discovered in an infrared search at Kitt Peak 5 hours after the burst trigger, this afterglow is amongst the faintest observed in the R-band at an early epoch, and exhibits very red colors, with $R-K\approx 6$. The magnitude of the optical afterglow of GRB 030115 is fainter than many upper limits for other bursts, suggesting that without early nIR observations it would have been classified as a ``dark'' burst. Both the color and optical magnitude of the afterglow are likely due to dust extinction and indicate that at least some optical afterglows are very faint due to dust along the line of sight. Multicolor {\it Hubble Space Telescope} observations were also taken of the host galaxy and the surrounding field. Photometric redshifts imply that the host, and a substantial number of faint galaxies in the field are at $z \sim 2.5$. The overdensity of galaxies is sufficiently great that GRB 030115 may have occurred in a rich high-redshift cluster. The host galaxy shows extremely red colors (R-K=5) and is the first GRB host to be classified as an Extremely Red Object (ERO). Some of the galaxies surrounding the host also show very red colors, while the majority of the cluster are much bluer, indicating ongoing unobscured star formation. As it is thought that much of high redshift star formation occurs in highly obscured environments it may be that GRB 030115 represent a transition object, between the relatively unobscured afterglows seen to date and a population which are very heavily extinguished, even in the nIR.

preprint2003arXiv

A Photometric and Spectroscopic Study of Dwarf and Giant Galaxies in the Coma Cluster - IV. The Luminosity Function

A large spectroscopic survey is constructed of galaxies in the Coma cluster. The survey covers a wide area (1 deg$^2$) to deep magnitudes (R ~ 19.5), covering both the core (high density) and outskirts (intermediate to low density) of the cluster. The spectroscopic sample consists of a total of 1191 galaxies, of which, 760 galaxies are confirmed members of the Coma cluster. A statistical technique is developed to correct the spectroscopic sample for incompletness. The corrected sample is then used to construct R-band luminosity function (LF) spanning a range of 7 magnitudes (-23 < M_R < -16) both at the core and outskirts of the cluster. Dependence of the LF on local environment in Coma is explored. The LFs are found to be the same, within the errors, between the inner and outer regions and close to those from recent measurements for field galaxies. The total B-band LF for the Coma cluster, fitted to a Schecter form is also measured and shows a dip at M(B) = -18 mag., in agreement with previous studies. The implications of this feature are discussed. The LF is studied in B-R color intervals and shows a steep faint-end slope for red (B-R > 1.35) galaxies, both at the core and outskirts of the cluster. This population of low luminosity red galaxies has a higher surface density than the blue (B-R < 1.35) star-forming population, dominating the faint-end of the Coma cluster LF. It is found that relative number of high surface brightness galaxies is larger at the cluster core, implying destruction of low surface brightness galaxies in dense core environment.

preprint2003arXiv

Identification of Type Ia Supernovae at Redshift 1.3 and Beyond with the Advanced Camera for Surveys on HST

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.

preprint1999arXiv

Galaxy Clustering at z ~ 2 and Halo Radii

The amplitude of the angular two-point galaxy correlation function w(θ) for galaxies at z~2 is estimated for galaxies in the Hubble Deep Field by using a U < 27 complete sub-sample. (i) It is confirmed that the amplitude of the correlation can be corrected for the integral constraint without having to make assumptions about the shape of the correlation function and by avoiding the introduction of linear error terms. The estimate using this technique is w(5'') = 0.10 \pm 0.09. (ii) If the biases introduced in faint galaxy selection due to obscuration by large objects are not corrected for by masking areas around them, then the estimate would be w(5'') =0.16\pm 0.07. (iii) The effective (3-D) galaxy pair separation at 5'' and this redshift range is ~ 25-250 /h kpc, so the correction to the spatial correlation function ξ(r) due to exclusion of overlapping galaxy dark matter haloes should be considered. For clustering stable in proper units in an Ω=1,λ=0 universe, our w(5\arcs) estimate (a) implies a present-day correlation length of r_0 ~ 2.6^{+1.1}_{-1.7}/h Mpc if halo overlapping is ignored, but (b) for a present-day correlation length of r_0=5.5/h Mpc implies that a typical halo exclusion radius is r_halo=70^{+420}_{-30}/h kpc. (iv) The decreasing correlation period (DCP) of a high initial bias in the spatial correlation function is not detected at this redshift. For an Ω=1,λ=0 universe and (proper) stable clustering, possible detections of the DCP in other work would imply that ξat redshifts greater than z_t = 1.7\pm0.9 would be [(1+z)/(1+z_t)]^{2.1\pm3.6} times higher than at z_t, which is consistent with our lack of a detection at z ~ 2.