Unsupervised feature-learning for galaxy SEDs with denoising autoencoders

With the increasing number of deep multi-wavelength galaxy surveys, the spectral energy distribution (SED) of galaxies has become an invaluable tool for studying the formation of their structures and their evolution. In this context, standard analysis relies on simple spectro-photometric selection criteria based on a few SED colors. If this fully supervised classification already yielded clear achievements, it is not optimal to extract relevant information from the data. In this article, we propose to employ very recent advances in machine learning, and more precisely in feature learning, to derive a data-driven diagram. We show that the proposed approach based on denoising autoencoders recovers the bi-modality in the galaxy population in an unsupervised manner, without using any prior knowledge on galaxy SED classification. This technique has been compared to principal component analysis (PCA) and to standard color/color representations. In addition, preliminary results illustrate that this enables the capturing of extra physically meaningful information, such as redshift dependence, galaxy mass evolution and variation over the specific star formation rate. PCA also results in an unsupervised representation with physical properties, such as mass and sSFR, although this representation separates out. less other characteristics (bimodality, redshift evolution) than denoising autoencoders.Comment: 11 pages and 15 figures. To be published in A&

Missing GRB host galaxies in deep mid-infrared observations: implications on the use of GRBs as star formation tracers

We report on the first mid-infrared observations of 16 GRB host galaxies performed with the Spitzer Space Telescope, and investigate the presence of evolved stellar populations and dust-enshrouded star-forming activity associated with GRBs. Only a very small fraction of our sample is detected by Spitzer, which is not consistent with recent works suggesting the presence of a GRB host population dominated by massive and strongly-starbursting galaxies (SFR >~ 100M[sun]yr^–1). Should the GRB hosts be representative of star-forming galaxies at high redshift, models of galaxy evolution indicate that >~ 50% of GRB hosts would be easily detected at the depth of our mid-infrared observations. Unless our sample suffers from a strong observational bias which remains to be understood, we infer in this context that the GRBs identified with the current techniques can not be directly used as unbiased probes of the global and integrated star formation history of the Universe

The galaxy density environment of gamma-ray burst host galaxies

We analyze cross-correlation functions between Gamma-Ray Burst (GRB) hosts and surrounding galaxies. We have used data obtained with the Very Large Telescope at Cerro Paranal (Chile), as well as public Hubble Space Telescope data. Our results indicate that Gamma-Ray Burst host galaxies do not reside in high galaxy density environments. Moreover, the host-galaxy cross-correlations show a relatively low amplitude. Our results are in agreement with the cross-correlation function between star-forming galaxies and surrounding objects in the HDF-N.Comment: 6 figures, accepted for publication in Ap

Infrared Spectral Energy Distributions of z~0.7 Star-Forming Galaxies

We analyze the infrared (IR) spectral energy distributions (SEDs) for 10micron < lambda(rest) < 100micron for ~600 galaxies at z~0.7 in the extended Chandra Deep Field South by stacking their Spitzer 24, 70 and 160micron images. We place interesting constraints on the average IR SED shape in two bins: the brightest 25% of z~0.7 galaxies detected at 24micron, and the remaining 75% of individually-detected galaxies. Galaxies without individual detections at 24micron were not well-detected at 70micron and 160micron even through stacking. We find that the average IR SEDs of z~0.7 star-forming galaxies fall within the diversity of z~0 templates. While dust obscuration Lir/Luv seems to be only a function of star formation rate (SFR; ~ Lir+Luv), not of redshift, the dust temperature of star-forming galaxies (with SFR ~ 10 solar mass per year) at a given IR luminosity was lower at z~0.7 than today. We suggest an interpretation of this phenomenology in terms of dust geometry: intensely star-forming galaxies at z~0 are typically interacting, and host dense centrally-concentrated bursts of star formation and warm dust temperatures. At z~0.7, the bulk of intensely star-forming galaxies are relatively undisturbed spirals and irregulars, and we postulate that they have large amounts of widespread lower-density star formation, yielding lower dust temperatures for a given IR luminosity. We recommend what IR SEDs are most suitable for modeling intermediate redshift galaxies with different SFRs.Comment: 12 pages, 11 figures, 2 tables, accepted for publication in Ap

A Far-infrared Characterization of 24 μm Selected Galaxies at 0 < z < 2.5 using Stacking at 70 μm and 160 μm in the COSMOS Field

We present a study of the average properties of luminous infrared galaxies detected directly at 24 μm in the COSMOS field using a median stacking analysis at 70 μm and 160 μm. Over 35,000 sources spanning 0 ≤ z ≤ 3 and 0.06 mJy ≤ S_(24) ≤ 3.0 mJy are stacked, divided into bins of both photometric redshift and 24 μm flux. We find no correlation of S_(70)/S_(24) flux density ratio with S_(24), but find that galaxies with higher S_(24) have a lower S_(160)/S_(24) flux density ratio. These observed ratios suggest that 24 μm selected galaxies have warmer spectral energy distributions (SEDs) at higher mid-IR fluxes, and therefore have a possible higher fraction of active galactic nuclei. Comparisons of the average S_(70)/S_(24) and S_(160)/S_(24) colors with various empirical templates and theoretical models show that the galaxies detected at 24 μm are consistent with "normal" star-forming galaxies and warm mid-IR galaxies such as Mrk 231, but inconsistent with heavily obscured galaxies such as Arp 220. We perform a χ^2 analysis to determine best-fit galactic model SEDs and total IR luminosities for each of our bins. We compare our results to previous methods of estimating L IR and find that previous methods show considerable agreement over the full redshift range, except for the brightest S_(24) sources, where they overpredict the bolometric IR luminosity at high redshift, most likely due to their warmer dust SED. We present a table that can be used as a more accurate and robust method for estimating bolometric infrared luminosity from 24 μm flux densities

CO J=2-1 line emission in cluster galaxies at z~1: fueling star formation in dense environments

We present observations of CO J=2-1 line emission in infrared-luminous cluster galaxies at z~1 using the IRAM Plateau de Bure Interferometer. Our two primary targets are optically faint, dust-obscured galaxies (DOGs) found to lie within 2 Mpc of the centers of two massive (>10^14 Msun) galaxy clusters. CO line emission is not detected in either DOG. We calculate 3-sigma upper limits to the CO J=2-1 line luminosities, L'_CO < 6.08x10^9 and < 6.63x10^9 K km/s pc^2. Assuming a CO-to-H_2 conversion factor derived for ultraluminous infrared galaxies in the local Universe, this translates to limits on the cold molecular gas mass of M_H_2 < 4.86x10^9 Msun and M_H_2 < 5.30x10^9 Msun. Both DOGs exhibit mid-infrared continuum emission that follows a power-law, suggesting that an AGN contributes to the dust heating. As such, estimates of the star formation efficiencies in these DOGs are uncertain. A third cluster member with an infrared luminosity, L_IR < 7.4x10^11 Lsun, is serendipitously detected in CO J=2-1 line emission in the field of one of the DOGs located roughly two virial radii away from the cluster center. The optical spectrum of this object suggests that it is likely an obscured AGN, and the measured CO line luminosity is L'_CO = (1.94 +/- 0.35)x10^10 K km/s pc^2, which leads to an estimated cold molecular gas mass M_H_2 = (1.55+/-0.28)x10^10 Msun. A significant reservoir of molecular gas in a z~1 galaxy located away from the cluster center demonstrates that the fuel can exist to drive an increase in star-formation and AGN activity at the outskirts of high-redshift clusters.Comment: 22 pages, 4 figures; accepted for publication in Ap

A MMT/Hectospec Redshift Survey of 24 Micron Sources in the Spitzer First Look Survey

We present a spectroscopic survey using the MMT/Hectospec fiber spectrograph of 24 micron sources selected with the Spitzer Space Telescope in the Spitzer First Look Survey. We report 1296 new redshifts for 24 micron sources, including 599 with f(24micron) > 1 mJy. Combined with 291 additional redshifts for sources from the Sloan Digital Sky Survey (SDSS), our observing program was highly efficient and is ~90% complete for i' 1 mJy, and is 35% complete for i' < 20.5 mag and 0.3 mJy < f(24micron) < 1 mJy. Our Hectospec survey includes 1078 and 168 objects spectroscopically classified as galaxies and QSOs, respectively. Combining the Hectospec and SDSS samples, we find 24 micron-selected galaxies to z < 0.98 and QSOs to z < 3.6, with mean redshifts of = 0.27 and =1.1. As part of this publication, we include the redshift catalogs and the reduced spectra; these are also available online (http://mips.as.arizona.edu/~papovich/fls) and through the NASA/IPAC Infrared Science Archive (http://irsa.ipac.caltech.edu).Comment: Accepted for publication in the Astronomical Journal, AASTEX format, 23 pages, 7 figures (some in color). This replacement is the accepted version, and includes minor changes from previous version. Data tables and spectra available at http://mips.as.arizona.edu/~papovich/fls or at http://irsa.ipac.caltech.ed

The FUV to Near-IR Morphologies of Luminous Infrared Galaxies in the GOALS Sample

We compare the morphologies of a sample of 20 LIRGs from the Great Observatories All-sky LIRG Survey (GOALS) in the FUV, B, I and H bands, using the Gini (G) and M20 parameters to quantitatively estimate the distribution and concentration of flux as a function of wavelength. HST images provide an average spatial resolution of ~80 pc. While our LIRGs can be reliably classified as mergers across the entire range of wavelengths studied here, there is a clear shift toward more negative M20 (more bulge-dominated) and a less significant decrease in G values at longer wavelengths. We find no correlation between the derived FUV G-M20 parameters and the global measures of the IR to FUV flux ratio, IRX. Given the fine resolution in our HST data, this suggests either that the UV morphology and IRX are correlated on very small scales, or that the regions emitting the bulk of the IR emission emit almost no FUV light. We use our multi-wavelength data to simulate how merging LIRGs would appear from z~0.5-3 in deep optical and near-infrared images such as the HUDF, and use these simulations to measure the G-M20 at these redshifts. Our simulations indicate a noticeable decrease in G, which flattens at z >= 2 by as much as 40%, resulting in mis-classifying our LIRGs as disk-like, even in the rest-frame FUV. The higher redshift values of M20 for the GOALS sources do not appear to change more than about 10% from the values at z~0. The change in G-M20 is caused by the surface brightness dimming of extended tidal features and asymmetries, and also the decreased spatial resolution which reduced the number of individual clumps identified. This effect, seen as early as z~0.5, could easily lead to an underestimate of the number of merging galaxies at high-redshift in the rest-frame FUV.Comment: Accepted for publication in the Astronomical Journal. The total page count is 15 pages with 13 figures and 1 Tabl