Source author record

Matthew L. N. Ashby

Matthew L. N. Ashby 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

43works
4topics
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

43 published item(s)

preprint2022arXiv

Rapid build-up of the stellar content in the protocluster core SPT2349$-$56 at $z\,{=}\,4.3$

The protocluster SPT2349$-$56 at $z\,{=}\,4.3$ contains one of the most actively star-forming cores known, yet constraints on the total stellar mass of this system are highly uncertain. We have therefore carried out deep optical and infrared observations of this system, probing rest-frame ultraviolet to infrared wavelengths. Using the positions of the spectroscopically-confirmed protocluster members, we identify counterparts and perform detailed source deblending, allowing us to fit spectral energy distributions in order to estimate stellar masses. We show that the galaxies in SPT2349$-$56 have stellar masses proportional to their high star-formation rates, consistent with other protocluster galaxies and field submillimetre galaxies (SMGs) around redshift 4. The galaxies in SPT2349$-$56 have on average lower molecular gas-to-stellar mass fractions and depletion timescales than field SMGs, although with considerable scatter. We construct the stellar-mass function for SPT2349$-$56 and compare it to the stellar-mass function of $z\,{=}\,1$ galaxy clusters, finding consistent shapes between the two. We measure rest-frame galaxy ultraviolet half-light radii from our HST-F160W imaging, finding that on average the galaxies in our sample are similar in size to typical star-forming galaxies at these redshifts. However, the brightest HST-detected galaxy in our sample, found near the luminosity-weighted centre of the protocluster core, remains unresolved at this wavelength. Hydrodynamical simulations predict that the core galaxies will quickly merge into a brightest cluster galaxy, thus our observations provide a direct view of the early formation mechanisms of this class of object.

preprint2022arXiv

The galaxy starburst/main-sequence bimodality over five decades in stellar mass at z ~ 3-6.5

We study the relation between stellar mass (M*) and star formation rate (SFR) for star-forming galaxies over approximately five decades in stellar mass (5.5 <~ log10(M*/Msun) <~ 10.5) at z ~ 3-6.5. This unprecedented coverage has been possible thanks to the joint analysis of blank non-lensed fields (COSMOS/SMUVS) and cluster lensing fields (Hubble Frontier Fields) which allow us to reach very low stellar masses. Previous works have revealed the existence of a clear bimodality in the SFR-M* plane with a star-formation Main Sequence and a starburst cloud at z ~ 4-5. Here we show that this bimodality extends to all star-forming galaxies and is valid in the whole redshift range z ~ 3-6.5. We find that starbursts constitute at least 20% of all star-forming galaxies with M* >~ 10^9 Msun at these redshifts and reach a peak of 40% at z=4-5. More importantly, 60% to 90% of the total SFR budget at these redshifts is contained in starburst galaxies, indicating that the starburst mode of star-formation is dominant at high redshifts. Almost all the low stellar-mass starbursts with log10(M*/Msun) <~ 8.5 have ages comparable to the typical timescales of a starburst event, suggesting that these galaxies are being caught in the process of formation. Interestingly, galaxy formation models fail to predict the starburst/main-sequence bimodality and starbursts overall, suggesting that the starburst phenomenon may be driven by physical processes occurring at smaller scales than those probed by these models.

preprint2020arXiv

A proto-pseudobulge in ESO 320-G030 fed by a massive molecular inflow driven by a nuclear bar

Galaxies with nuclear bars are believed to efficiently drive gas inward, generating a nuclear starburst and possibly an active galactic nucleus (AGN). We confirm this scenario for the isolated, double-barred, luminous infrared galaxy ESO 320-G030 based on an analysis of Herschel and ALMA spectroscopic observations. Herschel/PACS and SPIRE observations of ESO 320-G030 show absorption/emission in 18 lines of H2O, which we combine with the ALMA H2O 423-330 448 GHz line (Eupper~400 K) and continuum images to study the nuclear region. Radiative transfer models indicate that 3 nuclear components are required to account for the H2O and continuum data. An envelope, with R~130-150 pc, T_dust~50 K, and N_H2~2x10^{23} cm^{-2}, surrounds a nuclear disk with R~40 pc and tau_100um~1.5-3 (N_H2~2x10^{24} cm^{-2}) and an extremely compact (R~12 pc), warm (~100 K), and buried (tau_100um>5, N_H2>~5x10^{24} cm^{-2}) core component. The three nuclear components account for 70% of the galaxy L_IR (SFR~16-18 Msun yr^{-1}). The nucleus is fed by a molecular inflow observed in CO 2-1 with ALMA, which is associated with the nuclear bar. With decreasing radius (r=450-225 pc), the mass inflow rate increases up to ~20 Msun yr^{-1}, which is similar to the nuclear SFR, indicating that the starburst is sustained by the inflow. At lower r, the inflow is best probed by the far-infrared OH ground-state doublets, with an estimated inflow rate of ~30 Msun yr^{-1}. The short timescale of ~20 Myr for nuclear gas replenishment indicates quick secular evolution, and indicates that we are witnessing an intermediate stage (<100 Myr) proto-pseudobulge fed by a massive inflow that is driven by a strong nuclear bar. We also apply the H2O model to the Herschel far-infrared spectroscopic observations of H2^{18}O, OH, $^{18}OH, OH+, H2O^+, H3O^+, NH, NH2, NH3, CH, CH^+, ^{13}CH^+, HF, SH, and C3, and estimate their abundances.

preprint2016arXiv

Metallicity and Age of the Stellar Stream Around the Disk Galaxy NGC 5907

Stellar streams have become central to studies of the interaction histories of nearby galaxies. To characterize the most prominent parts of the stellar stream around the well-known nearby (d = 17 Mpc) edge-on disk galaxy NGC 5907, we have obtained and analyzed new, deep gri Subaru/Suprime-Cam and 3.6 micron Spitzer/Infrared Array Camera (IRAC) observations. Combining the near-infrared 3.6 micron data with visible-light images allows us to use a long wavelength baseline to estimate the metallicity and age of the stellar population along a ~60 kpc long segment of the stream. We have fitted the stellar spectral energy distribution (SED) with a single-burst stellar population synthesis model and we use it to distinguish between the proposed satellite accretion and minor/major merger formation models of the stellar stream around this galaxy. We conclude that a massive minor merger (stellar mass ratio of at least 1:8) can best account for the metallicity of -0.3 inferred along the brightest parts of the stream.

preprint2016arXiv

Quasar Variability in the Mid-Infrared

The Decadal IRAC Bootes Survey is a mid-IR variability survey of the ~9 sq. deg. of the NDWFS Bootes Field and extends the time baseline of its predecessor, the Spitzer Deep, Wide-Field Survey (SDWFS), from 4 to 10 years. The Spitzer Space Telescope visited the field five times between 2004 and 2014 at 3.6 and 4.5 microns. We provide the difference image analysis photometry for a half a million mostly extragalactic sources. In mid-IR color-color plane, sources with quasar colors constitute the largest variability class (75%), 16% of the variable objects have stellar colors and the remaining 9% have the colors of galaxies. Adding the fifth epoch doubles the number of variable active galactic nuclei (AGNs) for the same false positive rates as in SDWFS, or increases the number of sources by 20% while decreasing the false positive rates by factors of 2-3 for the same variability amplitude. We quantify the ensemble mid-IR variability of ~1500 spectroscopically confirmed AGNs using single power-law structure functions (SFs), which we find to be steeper (index $γ=0.45$) than in the optical ($γ=0.3$), leading to much lower amplitudes at short time-lags. This provides evidence for large emission regions, smoothing out any fast UV/optical variations, as the origin of infrared quasar variability. The mid-IR AGN SF slope $γ$ seems to be uncorrelated with both the luminosity and rest-frame wavelength, while the amplitude shows an anti-correlation with the luminosity and a correlation with the rest-frame wavelength.

preprint2016arXiv

The Evolution of the Galaxy Stellar Mass Function at z= 4-8: A Steepening Low-mass-end Slope with Increasing Redshift

We present galaxy stellar mass functions (GSMFs) at $z=$ 4-8 from a rest-frame ultraviolet (UV) selected sample of $\sim$4500 galaxies, found via photometric redshifts over an area of $\sim$280 arcmin$^2$ in the CANDELS/GOODS fields and the Hubble Ultra Deep Field. The deepest Spitzer/IRAC data yet-to-date and the relatively large volume allow us to place a better constraint at both the low- and high-mass ends of the GSMFs compared to previous space-based studies from pre-CANDELS observations. Supplemented by a stacking analysis, we find a linear correlation between the rest-frame UV absolute magnitude at 1500 Å ($M_{\rm UV}$) and logarithmic stellar mass ($\log M_*$) that holds for galaxies with $\log(M_*/M_{\odot}) \lesssim 10$. We use simulations to validate our method of measuring the slope of the $\log M_*$-$M_{\rm UV}$ relation, finding that the bias is minimized with a hybrid technique combining photometry of individual bright galaxies with stacked photometry for faint galaxies. The resultant measured slopes do not significantly evolve over $z=$ 4-8, while the normalization of the trend exhibits a weak evolution toward lower masses at higher redshift. We combine the $\log M_*$-$M_{\rm UV}$ distribution with observed rest-frame UV luminosity functions at each redshift to derive the GSMFs, finding that the low-mass-end slope becomes steeper with increasing redshift from $α=-1.55^{+0.08}_{-0.07}$ at $z=4$ to $α=-2.25^{+0.72}_{-0.35}$ at $z=8$. The inferred stellar mass density, when integrated over $M_*=10^8$-$10^{13} M_{\odot}$, increases by a factor of $10^{+30}_{-2}$ between $z=7$ and $z=4$ and is in good agreement with the time integral of the cosmic star formation rate density.

preprint2016arXiv

The extended disc and halo of the Andromeda galaxy observed with Spitzer-IRAC

We present the first results from an extended survey of the Andromeda galaxy (M31) using 41.1 h of observations by Spitzer-IRAC at 3.6 and 4.5$\,μ$m. This survey extends previous observations to the outer disc and halo, covering total lengths of 4$.\!\!^\circ$4 and 6$.\!\!^\circ$6 along the minor and major axes, respectively. We have produced surface brightness profiles by combining the integrated light from background-corrected maps with stellar counts from a new catalogue of point sources. Using auxiliary catalogues we have carried out a statistical analysis in colour-magnitude space to discriminate M31 objects from foreground Milky Way stars and background galaxies. The catalogue includes 426,529 sources, of which 66 per cent have been assigned probability values to identify M31 objects with magnitude depths of [3.6]$\,=\,$19.0$\,\pm\,$0.2, [4.5]$\,=\,$18.7$\,\pm\,$0.2. We discuss applications of our data for constraining the stellar mass and characterising point sources in the outer radii.

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

An Increasing Stellar Baryon Fraction in Bright Galaxies at High Redshift

Recent observations have shown that the characteristic luminosity of the rest-frame ultraviolet (UV) luminosity function does not significantly evolve at 4 < z < 7 and is approximately M*_UV ~ -21. We investigate this apparent non-evolution by examining a sample of 178 bright, M_UV < -21 galaxies at z=4 to 7, analyzing their stellar populations and host halo masses. Including deep Spitzer/IRAC imaging to constrain the rest-frame optical light, we find that M*_UV galaxies at z=4-7 have similar stellar masses of log(M/Msol)=9.6-9.9 and are thus relatively massive for these high redshifts. However, bright galaxies at z=4-7 are less massive and have younger inferred ages than similarly bright galaxies at z=2-3, even though the two populations have similar star formation rates and levels of dust attenuation. We match the abundances of these bright z=4-7 galaxies to halo mass functions from the Bolshoi Lambda-CDM simulation to estimate the halo masses. We find that the typical halo masses in ~M*_UV galaxies decrease from log(M_h/Msol)=11.9 at z=4 to log(M_h/Msol)=11.4 at z=7. Thus, although we are studying galaxies at a similar mass across multiple redshifts, these galaxies live in lower mass halos at higher redshift. The stellar baryon fraction in units of the cosmic mean Omega_b/Omega_m rises from 5.1% at z=4 to 11.7% at z=7; this evolution is significant at the ~3-sigma level. This rise does not agree with simple expectations of how galaxies grow, and implies that some effect, perhaps a diminishing efficiency of feedback, is allowing a higher fraction of available baryons to be converted into stars at high redshifts.

preprint2015arXiv

Merger Signatures in the Dynamics of Star-forming Gas

Spatially resolved kinematics have been used to determine the dynamical status of star-forming galaxies with ambiguous morphologies, and constrain the importance of galaxy interactions during the assembly of galaxies. However, measuring the importance of interactions or galaxy merger rates requires knowledge of the systematics in kinematic diagnostics and the visible time with merger indicators. We analyze the dynamics of star-forming gas in a set of binary merger hydrodynamic simulations with stellar mass ratios of 1:1 and 1:4. We find that the evolution of kinematic asymmetries traced by star-forming gas mirrors morphological asymmetries derived from mock optical images, in which both merger indicators show the largest deviation from isolated disks during strong interaction phases. Based on a series of simulations with various initial disk orientations, orbital parameters, gas fractions, and mass ratios, we find that the merger signatures are visible for ~0.2-0.4 Gyr with kinematic merger indicators but can be approximately twice as long for equal-mass mergers of massive gas-rich disk galaxies designed to be analogs of z~2-3 submillimeter galaxies. Merger signatures are most apparent after the second passage and before the black holes coalescence, but in some cases they persist up to several hundred Myr after coalescence. About 20-60% of the simulated galaxies are not identified as mergers during the strong interaction phase, implying that galaxies undergoing violent merging process do not necessarily exhibit highly asymmetric kinematics in their star-forming gas. The lack of identifiable merger signatures in this population can lead to an underestimation of merger abundances in star-forming galaxies, and including them in samples of star-forming disks may bias the measurements of disk properties such as intrinsic velocity dispersion.

preprint2015arXiv

Star-forming blue ETGs in two newly discovered galaxy overdensities in the HUDF at z=1.84 and 1.9: unveiling the progenitors of passive ETGs in cluster cores

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

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

Variations of the ISM Compactness Across the Main Sequence of Star-Forming Galaxies: Observations and Simulations

(abridged) The majority of star-forming galaxies follow a simple empirical correlation in the star formation rate (SFR) versus stellar mass ($M_*$) plane, usually referred to as the star formation Main Sequence (MS). Here we combine a set of hydro-dynamical simulations of interacting galactic disks with state-of-the-art radiative transfer codes to analyze how the evolution of mergers is reflected upon the properties of the MS. We present \textsc{Chiburst}, a Markov Chain Monte Carlo (MCMC) Spectral Energy Distribution (SED) code that fits the multi-wavelength, broad-band photometry of galaxies and derives stellar masses, star formation rates, and geometrical properties of the dust distribution. We apply this tool to the SEDs of simulated mergers and compare the derived results with the reference output from the simulations. Our results indicate that changes in the SEDs of mergers as they approach coalescence and depart from the MS are related to an evolution of dust geometry in scales larger than a few hundred parsecs. This is reflected in a correlation between the specific star formation rate (sSFR), and the compactness parameter $\mathcal{C}$, that parametrizes this geometry and hence the evolution of dust temperature ($T_{\rm{dust}}$) with time. As mergers approach coalescence, they depart from the MS and increase their compactness, which implies that moderate outliers of the MS are consistent with late-type mergers. By further applying our method to real observations of Luminous Infrared Galaxies (LIRGs), we show that the merger scenario is unable to explain these extreme outliers of the MS. Only by significantly increasing the gas fraction in the simulations are we able to reproduce the SEDs of LIRGs.

preprint2014arXiv

3C 220.3: a radio galaxy lensing a submillimeter galaxy

Herschel Space Observatory photometry and extensive multiwavelength followup have revealed that the powerful radio galaxy 3C 220.3 at z=0.685 acts as a gravitational lens for a background submillimeter galaxy (SMG) at z=2.221. At an observed wavelength of 1mm, the SMG is lensed into three distinct images. In the observed near infrared, these images are connected by an arc of 1.8" radius forming an Einstein half-ring centered near the radio galaxy. In visible light, only the arc is apparent. 3C 220.3 is the only known instance of strong galaxy-scale lensing by a powerful radio galaxy not located in a galaxy cluster and therefore it offers the potential to probe the dark matter content of the radio galaxy host. Lens modeling rejects a single lens, but two lenses centered on the radio galaxy host A and a companion B, separated by 1.5", provide a fit consistent with all data and reveal faint candidates for the predicted fourth and fifth images. The model does not require an extended common dark matter halo, consistent with the absence of extended bright X-ray emission on our Chandra image. The projected dark matter fractions within the Einstein radii of A (1.02") and B (0.61") are about 0.4 +/- 0.3 and 0.55 +/- 0.3. The mass to i-band light ratios of A and B, M/L ~ 8 +/- 4 Msun/Lsun, appear comparable to those of radio-quiet lensing galaxies at the same redshift in the CASTLES, LSD, and SL2S samples. The lensed SMG is extremely bright with observed f(250um) = 440mJy owing to a magnification factor mu~10. The SMG spectrum shows luminous, narrow CIV 154.9nm emission, revealing that the SMG houses a hidden quasar in addition to a violent starburst. Multicolor image reconstruction of the SMG indicates a bipolar morphology of the emitted ultraviolet (UV) light suggestive of cones through which UV light escapes a dust-enshrouded nucleus.

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

Deconstructing the Galaxy Stellar Mass Function with UKIDSS and CANDELS: the Impact of Colour, Structure and Environment

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

preprint2014arXiv

High-precision Photometric Redshifts from Spitzer/IRAC: Extreme [3.6]-[4.5] Colors Identify Galaxies in the Redshift Range z~6.6-6.9

One of the most challenging aspects of studying galaxies in the z>~7 universe is the infrequent confirmation of their redshifts through spectroscopy, a phenomenon thought to occur from the increasing opacity of the intergalactic medium to Lya photons at z>6.5. The resulting redshift uncertainties inhibit the efficient search for [C II] in z~7 galaxies with sub-mm instruments such as ALMA, given their limited scan speed for faint lines. One means by which to improve the precision of the inferred redshifts is to exploit the potential impact of strong nebular emission lines on the colors of z~4-8 galaxies as observed by Spitzer/IRAC. At z~6.8, galaxies exhibit IRAC colors as blue as [3.6]-[4.5] ~-1, likely due to the contribution of [O III]+Hb to the 3.6 mum flux combined with the absence of line contamination in the 4.5 mum band. In this paper we explore the use of extremely blue [3.6]-[4.5] colors to identify galaxies in the narrow redshift window z~6.6-6.9. When combined with an I-dropout criterion, we demonstrate that we can plausibly select a relatively clean sample of z~6.8 galaxies. Through a systematic application of this selection technique to our catalogs from all five CANDELS fields, we identify 20 probable z~6.6-6.9 galaxies. We estimate that our criteria select the ~50% strongest line emitters at z~6.8 and from the IRAC colors we estimate a typical [O III]+Hb rest-frame equivalent width of 1085A for this sample. The small redshift uncertainties on our sample make it particularly well suited for follow-up studies with facilities such as ALMA.

preprint2014arXiv

PdBI Cold Dust Imaging of Two Extremely Red H-[4.5]>4 Galaxies Discovered with SEDS and CANDELS

We report Plateau de Bure Interferometer (PdBI) 1.1 mm continuum imaging towards two extremely red H-[4.5]>4 (AB) galaxies at z>3, which we have previously discovered making use of Spitzer SEDS and Hubble Space Telescope CANDELS ultra-deep images of the Ultra Deep Survey field. One of our objects is detected on the PdBI map with a 4.3 sigma significance, corresponding to Snu(1.1mm)=(0.78 +/- 0.18) mJy. By combining this detection with the Spitzer 8 and 24 micron photometry for this source, and SCUBA2 flux density upper limits, we infer that this galaxy is a composite active galactic nucleus/star-forming system. The infrared (IR)-derived star formation rate is SFR~(200 +/- 100) Msun/yr, which implies that this galaxy is a higher-redshift analogue of the ordinary ultra-luminous infrared galaxies (ULIRGs) more commonly found at z~2-3. In the field of the other target, we find a tentative 3.1 sigma detection on the PdBI 1.1 mm map, but 3.7 arcsec away of our target position, so it likely corresponds to a different object. In spite of the lower significance, the PdBI detection is supported by a close SCUBA2 3.3 sigma detection. No counterpart is found on either the deep SEDS or CANDELS maps, so, if real, the PdBI source could be similar in nature to the sub-millimetre source GN10. We conclude that the analysis of ultra-deep near- and mid-IR images offers an efficient, alternative route to discover new sites of powerful star formation activity at high redshifts.

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

Simulated Galaxy Interactions as Probes of Merger Spectral Energy Distributions

We present the first systematic comparison of ultraviolet-millimeter spectral energy distributions (SEDs) of observed and simulated interacting galaxies. Our sample is drawn from the Spitzer Interacting Galaxy Survey, and probes a range of galaxy interaction parameters. We use 31 galaxies in 14 systems which have been observed with Herschel, Spitzer, GALEX, and 2MASS. We create a suite of GADGET-3 hydrodynamic simulations of isolated and interacting galaxies with stellar masses comparable to those in our sample of interacting galaxies. Photometry for the simulated systems is then calculated with the SUNRISE radiative transfer code for comparison with the observed systems. For most of the observed systems, one or more of the simulated SEDs match reasonably well. The best matches recover the infrared luminosity and the star formation rate of the observed systems, and the more massive systems preferentially match SEDs from simulations of more massive galaxies. The most morphologically distorted systems in our sample are best matched to simulated SEDs close to coalescence, while less evolved systems match well with SEDs over a wide range of interaction stages, suggesting that an SED alone is insufficient to identify interaction stage except during the most active phases in strongly interacting systems. This result is supported by our finding that the SEDs calculated for simulated systems vary little over the interaction sequence.

preprint2014arXiv

The HETDEX Pilot Survey V: The Physical Origin of Lyman-alpha Emitters Probed by Near-infrared Spectroscopy

We present the results from a VLT/SINFONI and Keck/NIRSPEC near-infrared spectroscopic survey of 16 Lyman-alpha emitters (LAEs) at $z$ = 2.1 - 2.5 in the COSMOS and GOODS-N fields discovered from the HETDEX Pilot Survey. We detect rest-frame optical nebular lines (H$α$ and/or [OIII]$λ$5007) for 10 of the LAEs and measure physical properties, including the star formation rate (SFR), gas-phase metallicity, gas-mass fraction, and Ly$α$ velocity offset. We find that LAEs may lie below the mass-metallicity relation for continuum-selected star-forming galaxies at the same redshift. The LAEs all show velocity shifts of Ly$α$ relative to the systemic redshift ranging between +85 and +296 km s$^{-1}$ with a mean of +180 km s$^{-1}$. This value is smaller than measured for continuum-selected star-forming galaxies at similar redshifts. The Ly$α$ velocity offsets show a moderate correlation with the measured star formation rate (2.5$σ$), but no significant correlations are seen with the SFR surface density, specific SFR, stellar mass, or dynamical mass ($\lesssim$ 1.5$σ$). Exploring the role of dust, kinematics of the interstellar medium (ISM), and geometry on the escape of Ly$α$ photons, we find no signature of selective quenching of resonantly scattered Ly$α$ photons. However, we also find no evidence that a clumpy ISM is enhancing the Ly$α$ equivalent width. Our results suggest that the low metallicity in LAEs may be responsible for yielding an environment with a low neutral hydrogen column density as well as less dust, easing the escape of Ly$α$ photons over that in continuum-selected star-forming galaxies.

preprint2014arXiv

The mass evolution of the first galaxies: stellar mass functions and star formation rates at $4 < z < 7$ in the CANDELS GOODS-South field

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

preprint2014arXiv

The Spitzer South Pole Telescope Deep Field Survey: Linking galaxies and halos at z=1.5

We present an analysis of the clustering of high-redshift galaxies in the recently completed 94 deg$^2$ Spitzer-SPT Deep Field survey. Applying flux and color cuts to the mid-infrared photometry efficiently selects galaxies at $z\sim1.5$ in the stellar mass range $10^{10}-10^{11}M_\odot$, making this sample the largest used so far to study such a distant population. We measure the angular correlation function in different flux-limited samples at scales $>6^{\prime \prime}$ (corresponding to physical distances $>0.05$ Mpc) and thereby map the one- and two-halo contributions to the clustering. We fit halo occupation distributions and determine how the central galaxy's stellar mass and satellite occupation depend on the halo mass. We measure a prominent peak in the stellar-to-halo mass ratio at a halo mass of $\log(M_{\rm halo} / M_\odot) = 12.44\pm0.08$, 4.5 times higher than the $z=0$ value. This supports the idea of an evolving mass threshold above which star formation is quenched. We estimate the large-scale bias in the range $b_g=2-4$ and the satellite fraction to be $f_\mathrm{sat}\sim0.2$, showing a clear evolution compared to $z=0$. We also find that, above a given stellar mass limit, the fraction of galaxies that are in similar mass pairs is higher at $z=1.5$ than at $z=0$. In addition, we measure that this fraction mildly increases with the stellar mass limit at $z=1.5$, which is the opposite of the behavior seen at low-redshift.

preprint2014arXiv

The total infrared luminosity may significantly overestimate the star formation rate of recently quenched galaxies

The total infrared (IR) luminosity is very useful for estimating the star formation rate (SFR) of galaxies, but converting the IR luminosity into an SFR relies on assumptions that do not hold for all galaxies. We test the effectiveness of the IR luminosity as an SFR indicator by applying it to synthetic spectral energy distributions generated from three-dimensional hydrodynamical simulations of isolated disc galaxies and galaxy mergers. In general, the SFR inferred from the IR luminosity agrees well with the true instantaneous SFR of the simulated galaxies. However, for the major mergers in which a strong starburst is induced, the SFR inferred from the IR luminosity can overestimate the instantaneous SFR during the post-starburst phase by greater than two orders of magnitude. Even though the instantaneous SFR decreases rapidly after the starburst, the stars that were formed in the starburst remain dust-obscured and thus produce significant IR luminosity. Consequently, use of the IR luminosity as an SFR indicator may cause one to conclude that post-starburst galaxies are still star-forming, whereas in reality, star formation was recently quenched.

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

Global Star Formation Rates and Dust Emission Over the Galaxy Interaction Sequence

We measured and modeled the spectral energy distributions (SEDs) in 28 bands from the ultraviolet to the far-infrared (FIR) for 31 interacting galaxies in 14 systems. The sample is drawn from the Spitzer Interacting Galaxy Survey, which probes a range of galaxy interaction parameters at multiple wavelengths with an emphasis on the infrared bands. The subset presented in this paper consists of all galaxies for which FIR Herschel SPIRE observations are publicly available. Our SEDs combine the Herschel photometry with multi-wavelength data from Spitzer, GALEX, Swift UVOT, and 2MASS. While the shapes of the SEDs are broadly similar across our sample, strongly interacting galaxies typically have more mid-infrared emission relative to their near-infrared and FIR emission than weakly or moderately interacting galaxies. We modeled the full SEDs to derive host galaxy star formation rates (SFR), specific star formation rates (sSFR), stellar masses, dust temperatures, dust luminosities, and dust masses. We find increases in the dust luminosity and mass, SFR, and cold (15-25 K) dust temperature as the interaction progresses from moderately to strongly interacting and between non- interacting and strongly interacting galaxies. We also find increases in the SFR between weakly and strongly interacting galaxies. In contrast, the sSFR remains unchanged across all the interaction stages. The SFR derived from the SED modeling agrees well with rates estimated by proportionality relations that depend on infrared emission. (abridged)

preprint2013arXiv

To Stack or Not to Stack: Spectral Energy Distribution Properties of Lya-Emitting Galaxies at z=2.1

We use the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS) GOODS-S multi-wavelength catalog to identify counterparts for 20 Lya Emitting (LAE) galaxies at z=2.1. We build several types of stacked Spectral Energy Distributions (SEDs) of these objects. We combine photometry to form average and median flux-stacked SEDs, and postage stamp images to form average and median image-stacked SEDs. We also introduce scaled flux stacks that eliminate the influence of variation in overall brightness. We use the SED fitting code SpeedyMC to constrain the physical properties of individual objects and stacks. Our LAEs at z = 2.1 have stellar masses ranging from 2x10^7 Msun - 8x10^9 Msun (median = 3x10^8 Msun), ages ranging from 4 Myr to 500 Myr (median =100 Myr), and E(B-V) between 0.02 and 0.24 (median = 0.12). We do not observe strong correlations between Lya equivalent width (EW) and stellar mass, age, or E(B-V). The Lya radiative transfer (q) factors of our sample are predominantly close to one and do not correlate strongly with EW or E(B-V), implying that Lya radiative transfer prevents Lya photons from resonantly scattering in dusty regions. The SED parameters of the flux stacks match the average and median values of the individual objects, with the flux-scaled median SED performing best with reduced uncertainties. Median image-stacked SEDs provide a poor representation of the median individual object, and none of the stacking methods captures the large dispersion of LAE properties.

preprint2013arXiv

Tracing the evolution of Active Galactic Nuclei host galaxies over the last 9 Gyrs of Cosmic time

We present the results of a combined galaxy population analysis for the host galaxies of active galactic nuclei (AGN) identified at 0 < z < 1.4 within the SDSS, Bootes and DEEP2 surveys. We identified AGN in a uniform and unbiased manner at X-ray, infrared and radio wavelengths. Supermassive black holes undergoing radiatively-efficient accretion (detected as X-ray and/or infrared AGN) appear to be hosted in a separate and distinct galaxy population than AGN undergoing powerful mechanically dominated accretion (radio AGN). Consistent with some previous studies, radiatively efficient AGN appear to be preferentially hosted in modest star-forming galaxies, with little dependence on AGN or galaxy luminosity. AGN exhibiting radio-emitting jets due to mechanically-dominated accretion are almost exclusively observed in massive, passive galaxies. Crucially, we now provide strong evidence that the observed host-galaxy trends are independent of redshift. In particular, these different accretion-mode AGN have remained as separate galaxy populations throughout the last 9 Gyr. Furthermore, it appears that galaxies hosting AGN have evolved along the same path as galaxies that are not hosting AGN with little evidence for distinctly separate evolution.

preprint2012arXiv

A measurement of the intrahalo light fraction with near-infrared background anisotropies

Unresolved near-infrared background anisotropies are expected to have contributions from the earliest galaxies during reionization and faint, dwarf galaxies at intermediate redshifts. Previous measurements were unable to conclusively pinpoint the dominant origin because they did not sample spatial scales that were sufficiently large to distinguish between these two possibilities. Here we report a measurement of the anisotropy power spectrum from sub-arcminute to one degree angular scales and find the clustering amplitude to be larger than the model predictions involving the two existing explanations. As the shot-noise level of the power spectrum is consistent with that expected from faint galaxies, a new source population on the sky is not necessary to explain the observations. A physical mechanism that increases the clustering amplitude, however, is needed. Motivated by recent results related to the extended stellar light profile in dark matter halos, we consider the possibility that the fluctuations originate from diffuse intrahalo stars of all galaxies. We find that the measured power spectrum can be explained by an intrahalo light fraction of 0.07 to 0.2 % relative to the total luminosity in dark matter halos of masses log(M/M_Sun) ~ 9 to 12 at redshifts of ~ 1 to 4.

preprint2012arXiv

The Clustering of Extremely Red Objects

We measure the clustering of Extremely Red Objects (EROs) in ~8 deg^2 of the NOAO Deep Wide Field Survey Boötes field in order to establish robust links between ERO z~1.2 and local galaxy z<0.1 populations. Three different color selection criteria from the literature are analyzed to assess the consequences of using different criteria for selecting EROs. Specifically, our samples are (R-K_s)>5.0 (28,724 galaxies), (I-K_s)>4.0 (22,451 galaxies) and (I-[3.6])>5.0 (64,370 galaxies). Magnitude-limited samples show the correlation length (r_0) to increase for more luminous EROs, implying a correlation with stellar mass. We can separate star-forming and passive ERO populations using the (K_s-[24]) and ([3.6]-[24]) colors to K_s=18.4 and [3.6]=17.5, respectively. Star-forming and passive EROs in magnitude limited samples have different clustering properties and host dark halo masses, and cannot be simply understood as a single population. Based on the clustering, we find that bright passive EROs are the likely progenitors of >4L^* elliptical galaxies. Bright EROs with ongoing star formation were found to occupy denser environments than star-forming galaxies in the local Universe, making these the likely progenitors of >L^* local ellipticals. This suggests that the progenitors of massive >4L^* local ellipticals had stopped forming stars by z>1.2, but that the progenitors of less massive ellipticals (down to L^*) can still show significant star formation at this epoch.

preprint2012arXiv

The Faint End of the Cluster Galaxy Luminosity Function at High Redshift

We measure the faint end slope of the galaxy luminosity function (LF) for cluster galaxies at 1<z<1.5 using Spitzer IRAC data. We investigate whether this slope, alpha, differs from that of the field LF at these redshifts, and with the cluster LF at low redshifts. The latter is of particular interest as low-luminosity galaxies are expected to undergo significant evolution. We use seven high-redshift spectroscopically confirmed galaxy clusters drawn from the IRAC Shallow Cluster Survey to measure the cluster galaxy LF down to depths of M* + 3 (3.6 microns) and M* + 2.5 (4.5 microns). The summed LF at our median cluster redshift (z=1.35) is well fit by a Schechter distribution with alpha[3.6] = -0.97 +/- 0.14 and alpha[4.5] = -0.91 +/- 0.28, consistent with a flat faint end slope and is in agreement with measurements of the field LF in similar bands at these redshifts. A comparison to alpha in low-redshift clusters finds no statistically significant evidence of evolution. Combined with past studies which show that M* is passively evolving out to z~1.3, this means that the shape of the cluster LF is largely in place by z~1.3. This suggests that the processes that govern the build up of the mass of low-mass cluster galaxies have no net effect on the faint end slope of the cluster LF at z<1.3.

preprint2012arXiv

The Nature of Extremely Red H-[4.5]>4 Galaxies revealed with SEDS and CANDELS

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

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

The Canada-France Ecliptic Plane Survey - Full Data Release: The orbital structure of the Kuiper belt

We report the orbital distribution of the trans-neptunian objects (TNOs) discovered during the Canada-France Ecliptic Plane Survey, whose discovery phase ran from early 2003 until early 2007. The follow-up observations started just after the first discoveries and extended until late 2009. We obtained characterized observations of 321 sq.deg. of sky to depths in the range g ~ 23.5--24.4 AB mag. We provide a database of 169 TNOs with high-precision dynamical classification and known discovery efficiency. Using this database, we find that the classical belt is a complex region with sub-structures that go beyond the usual splitting of inner (interior to 3:2 mean-motion resonance [MMR]), outer (exterior to 2:1 MMR), and main (in between). The main classical belt (a=40--47 AU) needs to be modeled with at least three components: the `hot' component with a wide inclination distribution and two `cold' components (stirred and kernel) with much narrower inclination distributions. The hot component must have a significantly shallower absolute magnitude (Hg) distribution than the other two components. With 95% confidence, there are 8000+1800-1600 objects in the main belt with Hg <= 8.0, of which 50% are from the hot component, 40% from the stirred component and 10% from the kernel; the hot component's fraction drops rapidly with increasing Hg. Because of this, the apparent population fractions depend on the depth and ecliptic latitude of a trans-neptunian survey. The stirred and kernel components are limited to only a portion of the main belt, while we find that the hot component is consistent with a smooth extension throughout the inner, main and outer regions of the classical belt; the inner and outer belts are consistent with containing only hot-component objects. The Hg <= 8.0 TNO population estimates are 400 for the inner belt and 10,000 for the outer belt within a factor of two.

preprint2011arXiv

The SEDs of Interacting Galaxies

The evolution of galaxies is greatly influenced by their interactions. As part of a program to study interacting galaxies, we have measured and modeled the spectral energy distri- butions (SEDs) from the ultraviolet (UV) to the far-infrared (FIR). We describe the constraints imposed on star formation histories by these SEDs, and the variations therein seen across the interaction sequence, and we compare the results of different star formation rate prescriptions applied to the data. The sample itself is based on the Spitzer Interacting Galaxy Survey (SIGS) of 111 galaxies in 50 systems, a project designed to probe a range of galaxy interaction parameters in the infrared. Our SEDs combine the Spitzer results with multiwavelength data from other missions, in particular GALEX and Herschel. The subset presented here is the sample for which FIR Herschel observations are currently publicly available.

preprint2011arXiv

UV-to-FIR analysis of Spitzer/IRAC sources in the Extended Groth Strip II: Photometric redshifts, Stellar masses and Star formation rates

Based on the ultraviolet to far-infrared photometry already compiled and presented in a companion paper (Barro et al. 2011a, Paper I), we present a detailed SED analysis of nearly 80,000 IRAC 3.6+4.5 micron selected galaxies in the Extended Groth Strip. We estimate photometric redshifts, stellar masses, and star formation rates separately for each galaxy in this large sample. The catalog includes 76,936 sources with [3.6] < 23.75 (85% completeness level of the IRAC survey) over 0.48 square degrees. The typical photometric redshift accuracy is Delta z/(1+z)=0.034, with a catastrophic outlier fraction of just 2%. We quantify the systematics introduced by the use of different stellar population synthesis libraries and IMFs in the calculation of stellar masses. We find systematic offsets ranging from 0.1 to 0.4 dex, with a typical scatter of 0.3 dex. We also provide UV- and IR-based SFRs for all sample galaxies, based on several sets of dust emission templates and SFR indicators. We evaluate the systematic differences and goodness of the different SFR estimations using the deep FIDEL 70 micron data available in the EGS. Typical random uncertainties of the IR-bases SFRs are a factor of two, with non-negligible systematic effects at z$\gtrsim$1.5 observed when only MIPS 24 micron data is available. All data products (SEDs, postage stamps from imaging data, and different estimations of the photometric redshifts, stellar masses, and SFRs of each galaxy) described in this and the companion paper are publicly available, and they can be accessed through our the web-interface utility Rainbow-navigator

preprint2010arXiv

Constraining the Outburst Properties of the SMBH in Fornax A through X-ray, Infrared, and Radio Observations

Combined Spitzer, Chandra, XMM-Newton, and VLA observations of the giant radio galaxy NGC 1316 (Fornax A) show a radio jet and X-ray cavities from AGN outbursts likely triggered by a merger with a late-type galaxy at least 0.4 Gyr ago. We detect a weak nucleus with an SED typical of a low-luminosity AGN with a bolometric luminosity of 2.4x10^42 erg/s. The Spitzer IRAC and MIPS images show dust emission strongest in regions with little or no radio emission. The large infrared luminosity relative to the galaxy's K-band luminosity implies an external origin for the dust. The dust mass implies that the merger spiral galaxy had a stellar mass of 1-6x10^10 M_sun and a gas mass of 2-4x10^9 M_sun. Chandra images show a small ~15"(1.6 kpc) cavity coincident with the radio jet, while the XMM-Newton image shows two large X-ray cavities lying 320"(34.8 kpc) east and west of the nucleus, each approximately 230"(25 kpc) in radius. The radio lobes lie at radii of 14.3'(93 kpc) and 15.6'(101 kcp), more distant from the nucleus than the detected X-ray cavities. The relative morphology of the large scale 1.4 GHz and X-ray emission suggests they were products of two distinct outbursts, an earlier one creating the radio lobes and a later one producing the X-ray cavities. Alternatively, if a single outburst created both the X-ray cavities and the radio lobes, this would require that the morphology is not fully defined by the 1.4 GHz emission. For the more likely two outburst scenario, we use the buoyancy rise times to estimate an age for the more recent outburst of 0.1 Gyr and use the associated PV work done by the expanding plasma to create the X-ray cavities to estimate the outburst's energy of 10^58 ergs. The present size and location of the large radio lobes implies that the AGN outburst that created them happened ~0.4 Gyr ago and released ~5x10^58 ergs. (abridged)

preprint2010arXiv

Infrared Properties of a Complete Sample of Star-Forming Dwarf Galaxies

We present a study of a large, statistically complete sample of star-forming dwarf galaxies using mid-infrared observations from the {\it Spitzer Space Telescope}. The relationships between metallicity, star formation rate (SFR) and mid-infrared color in these systems show that the galaxies span a wide range of properties. However, the galaxies do show a deficit of 8.0 \um\ polycyclic aromatic hydrocarbon emission as is apparent from the median 8.0 \um\ luminosity which is only 0.004 \lstarf\ while the median $B$-band luminosity is 0.05 \lstarb. Despite many of the galaxies being 8.0 \um\ deficient, there is about a factor of 4 more extremely red galaxies in the [3.6] $-$ [8.0] color than for a sample of normal galaxies with similar optical colors. We show correlations between the [3.6] $-$ [8.0] color and luminosity, metallicity, and to a lesser extent SFRs that were not evident in the original, smaller sample studied previously. The luminosity--metallicity relation has a flatter slope for dwarf galaxies as has been indicated by previous work. We also show a relationship between the 8.0 \um\ luminosity and the metallicity of the galaxy which is not expected given the competing effects (stellar mass, stellar population age, and the hardness of the radiation field) that influence the 8.0 \um\ emission. This larger sample plus a well-defined selection function also allows us to compute the 8.0 \um\ luminosity function and compare it with the one for the local galaxy population. Our results show that below 10$^{9}$ $L$\solar, nearly all the 8.0 \um\ luminosity density of the local universe arises from dwarf galaxies that exhibit strong \ha\ emission -- i.e., 8.0 \um\ and \ha\ selection identify similar galaxy populations despite the deficit of 8.0 \um\ emission observed in these dwarfs.

preprint2010arXiv

Mid-Infrared Variability from the Spitzer Deep, Wide-Field Survey

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

preprint2010arXiv

Ultracool Field Brown Dwarf Candidates Selected at 4.5 microns

We have identified a sample of cool field brown dwarf candidates using IRAC data from the Spitzer Deep, Wide-Field Survey (SDWFS). The candidates were selected from 400,000 SDWFS sources with [4.5] <= 18.5 mag and required to have [3.6]-[4.5] >= 1.5 and [4.5] - [8.0] <= 2.0 on the Vega system. The first color requirement selects objects redder than all but a handful of presently known brown dwarfs with spectral classes later than T7, while the second eliminates 14 probable reddened AGN. Optical detection of 4 of the remaining 18 sources implies they are likely also AGN, leaving 14 brown dwarf candidates. For two of the brightest candidates (SDWFS J143524.44+335334.6 and SDWFS J143222.82+323746.5), the spectral energy distributions including near-infrared detections suggest a spectral class of ~ T8. The proper motion is < 0.25 "/yr, consistent with expectations for a luminosity inferred distance of >70 pc. The reddest brown dwarf candidate (SDWFS J143356.62+351849.2) has [3.6] - [4.5]=2.24 and H - [4.5] > 5.7, redder than any published brown dwarf in these colors, and may be the first example of the elusive Y-dwarf spectral class. Models from Burrows et al. (2003) predict larger numbers of cool brown dwarfs should be found for a Chabrier (2003) mass function. Suppressing the model [4.5] flux by a factor of two, as indicated by previous work, brings the Burrows models and observations into reasonable agreement. The recently launched Wide-field Infrared Survey Explorer (WISE) will probe a volume ~40x larger and should find hundreds of brown dwarfs cooler than T7.

preprint2010arXiv

Warm Molecular Hydrogen Emission in Normal Edge-On Galaxies NGC 4565 and NGC 5907

We have observed warm molecular hydrogen in two nearby edge-on disk galaxies, NGC 4565 and NGC 5907, using the Spitzer high-resolution infrared spectrograph. The 0-0 S(0) 28.2 micron and 0-0 S(1) 17.0 micron pure rotational lines were detected out to 10 kpc from the center of each galaxy on both sides of the major axis, and in NGC 4565 the S(0) line was detected at r = 15 kpc on one side. This location lies beyond a steep drop in the radio continuum emission from cosmic rays in the disk. Despite indications that star formation activity decreases with radius, the H2 excitation temperature and the ratio of the H2 line and the far-IR luminosity surface densities, Sigma_L(H2}/Sigma_L(TIR}, change very little as a function of radius, even into the diffuse outer region of the disk of NGC 4565. This suggests that the source of excitation of the H2 operates over a large range of radii, and is broadly independent of the strength and relative location of UV emission from young stars. Although excitation in photodissociation regions is the most common explanation for the widespread H2 emission, cosmic ray heating or shocks cannot be ruled out. The inferred mass surface densities of warm molecular hydrogen in both edge-on galaxies differ substantially, being 4(-60) M_solar/pc^2 and 3(-50) M_solar/pc^2 at r = 10 kpc for NGC 4565 and NGC 5907, respectively. The higher values represent very unlikely point-source upper limits. The point source case is not supported by the observed emission distribution in the spectral slits. These mass surface densities cannot support the observed rotation velocities in excess of 200 km/s. Therefore, warm molecular hydrogen cannot account for dark matter in these disk galaxies, contrary to what was implied by a previous ISO study of the nearby edge-on galaxy NGC 891.