No Arabic abstract
Galaxies selected at 170um by the ISO FIRBACK survey represent the brightest ~10% of the Cosmic Infrared Background. Examining their nature in detail is therefore crucial for constraining models of galaxy evolution. Here we combine Spitzer archival data with previous near-IR, far-IR, and sub-mm observations of a representative sample of 22 FIRBACK galaxies spanning three orders of magnitude in infrared luminosity. We fit a flexible, multi-component, empirical SED model of star-forming galaxies designed to model the entire ~1-1000um wavelength range. The fits are performed with a Markov Chain Monte Carlo (MCMC) approach, allowing for meaningful uncertainties to be derived. This approach also highlights degeneracies such as between Td and beta, which we discuss in detail. From these fits and standard relations we derive: L_IR, L_PAH, SFR, tau_V, M_star, M_dust, Td, and beta. We look at a variety of correlations between these and combinations thereof in order to examine the physical nature of these galaxies. Our conclusions are supplemented by morphological examination of the sources, and comparison with local samples. We find the bulk of our sample to be consistent with fairly standard size and mass disk galaxies with somewhat enhanced star-formation relative to local spirals, but likely not bona fide starbursts. A few higher-z LIGs and ULIGs are also present, but contrary to expectation, they are weak mid-IR emitters and overall are consistent with star-formation over an extended cold region rather than concentrated in the nuclear regions. We discuss the implications of this study for understanding populations detected at other wavelengths, such as the bright 850um SCUBA sources or the faint Spitzer 24um sources.
We compare observed far infra-red/sub-millimetre (FIR/sub-mm) galaxy spectral energy distributions (SEDs) of massive galaxies ($M_{star}gtrsim10^{10}$ $h^{-1}$M$_{odot}$) derived through a stacking analysis with predictions from a new model of galaxy formation. The FIR SEDs of the model galaxies are calculated using a self-consistent model for the absorption and re-emission of radiation by interstellar dust based on radiative transfer calculations and global energy balance arguments. Galaxies are selected based on their position on the specific star formation rate (sSFR) - stellar mass ($M_{star}$) plane. We identify a main sequence of star-forming galaxies in the model, i.e. a well defined relationship between sSFR and $M_star$, up to redshift $zsim6$. The scatter of this relationship evolves such that it is generally larger at higher stellar masses and higher redshifts. There is remarkable agreement between the predicted and observed average SEDs across a broad range of redshifts ($0.5lesssim zlesssim4$) for galaxies on the main sequence. However, the agreement is less good for starburst galaxies at $zgtrsim2$, selected here to have elevated sSFRs$>10times$ the main sequence value. We find that the predicted average SEDs are robust to changing the parameters of our dust model within physically plausible values. We also show that the dust temperature evolution of main sequence galaxies in the model is driven by star formation on the main sequence being more burst-dominated at higher redshifts.
The dominant source of electromagnetic energy in the Universe today (over ultraviolet, optical and near-infrared wavelengths) is starlight. However, quantifying the amount of starlight produced has proven difficult due to interstellar dust grains which attenuate some unknown fraction of the light. Combining a recently calibrated galactic dust model with observations of 10,000 nearby galaxies we find that (integrated over all galaxy types and orientations) only (11 +/- 2)% of the 0.1 micron photons escape their host galaxies; this value rises linearly (with log(lambda)) to (87 +/- 3)% at 2.1 micron. We deduce that the energy output from stars in the nearby Universe is (1.6+/-0.2) x 10^{35} W Mpc^{-3} of which (0.9+/-0.1) x 10^{35} W Mpc^{-3} escapes directly into the inter-galactic medium. Some further ramifications of dust attenuation are discussed, and equations that correct individual galaxy flux measurements for its effect are provided.
We present the first large, unbiased sample of Lyman Break Galaxies (LBGs) at z ~ 1. Far ultraviolet-dropout (1530 A) galaxies in the Chandra Deep Field South have been selected using GALEX data. This first large sample in the z ~ 1 universe provides us with a high quality reference sample of LBGs. We analyzed the sample from the UV to the IR using GALEX, SPITZER, ESO and HST data. The morphology (obtained from GOODS data) of 75 % of our LBGs is consistent with a disk. The vast majority of LBGs with an IR detection are also Luminous Infrared Galaxies (LIRGs). As a class, the galaxies not detected at 24 microns are an order of magnitude fainter relative to the UV compared with those detected individually, suggesting that there may be two types of behavior within the sample. For the IR-bright galaxies, there is an apparent upper limit for the UV dust attenuation and this upper limit is anti-correlated with the observed UV luminosity. Previous estimates of dust attenuations based on the ultraviolet slope are compared to new ones based on the FIR/UV ratio (for LBGs detected at 24 microns), which is usually a more reliable estimator. Depending on the calibration we use to estimate the total IR luminosity, beta-based attenuations A_{FUV} are larger by 0.2 to 0.6 mag. than the ones estimated from FIR/UV ratio. Finally, for IR-bright LBGs, median estimated beta-based SFRs are 2-3 times larger than the total SFRs estimated as SFR_{TOT} = SFR_{UV} + SFR_{IR} while IR-based SFRs provide values below SFR_{TOT} by 15 - 20 %. We use a stacking method to statistically constrain the 24 microns flux of LBGs non individually detected. The results suggest that these LBGs do not contain large amounts of dust.
We present the results of IRS low-resolution spectroscopy of 51 Seyfert galaxies, part of a large Spitzer observing program to determine the mid-to-far infrared spectral energy distributions of a well-defined sample of 87 nearby, 12 micron-selected Seyferts. We find that the spectra clearly divide into groups based on their continuum shapes and spectral features. The infrared spectral types appear to be related to the Seyfert types. Some features are clearly related to a starburst contribution to the IR spectrum, while the observed power-law continuum shapes, attributed to the AGN, may be dust or non-thermal emission. Principal component analysis results suggest that the relative contribution of starburst emission is the dominant cause of variance in the spectra. We find that the Sy 2s show on average stronger starburst contributions than the Sy 1s.
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 sigma). 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 rho(L, z) = (1 + z)^q rho(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.