Breaking the curve with candels: a Bayesian approach to reveal the non-universality of the dust-attenuation law at high redshift

Dust attenuation affects nearly all observational aspects of galaxy evolution, yet very little is known about the form of the dust-attenuation law in the distant universe. Here, we model the spectral energy distributions of galaxies at z ~ 1.5–3 from CANDELS with rest-frame UV to near-IR imaging under different assumptions about the dust law, and compare the amount of inferred attenuated light with the observed infrared (IR) luminosities. Some individual galaxies show strong Bayesian evidence in preference of one dust law over another, and this preference agrees with their observed location on the plane of infrared excess (IRX, L_TIR/L_UV) and UV slope (β). We generalize the shape of the dust law with an empirical model, A_ λ,σ =E(B-V)k_ λ (λ / λ v)^ σ where k_λ is the dust law of Calzetti et al., and show that there exists a correlation between the color excess E(B-V) and tilt δ with δ =(0.62±0.05)log(E(B-V))+(0.26±0.02). Galaxies with high color excess have a shallower, starburst-like law, and those with low color excess have a steeper, SMC-like law. Surprisingly, the galaxies in our sample show no correlation between the shape of the dust law and stellar mass, star formation rate, or β. The change in the dust law with color excess is consistent with a model where attenuation is caused by scattering, a mixed star–dust geometry, and/or trends with stellar population age, metallicity, and dust grain size. This rest-frame UV-to-near-IR method shows potential to constrain the dust law at even higher redshifts (z>3).

Far-infrared constraints on the contamination by dust-obscured galaxies of high-z dropout searches

The spectral energy distributions (SED) of dusty galaxies at intermediate redshift may look similar to very high-redshift galaxies in the optical/near infrared (NIR) domain. This can lead to the contamination of high-redshift galaxy searches based on broad-band optical/NIR photometry by lower redshift dusty galaxies because both kind of galaxies cannot be distinguished. The contamination rate could be as high as 50%. This work shows how the far-infrared (FIR) domain can help to recognize likely low-z interlopers in an optical/NIR search for high-z galaxies. We analyze the FIR SEDs of two galaxies that are proposed to be very high-redshift (z > 7) dropout candidates based on deep Hawk-I/VLT observations. The FIR SEDs are sampled with PACS/Herschel at 100 and 160 μm, with SPIRE/Herschel at 250, 350 and 500 μm and with LABOCA/APEX at 870 μm. We find that redshifts > 7 would imply extreme FIR SEDs (with dust temperatures >100 K and FIR luminosities >10^13 L_⊙). At z ~ 2, instead, the SEDs of both sources would be compatible with those of typical ultra luminous infrared galaxies or submillimeter galaxies. Considering all available data for these sources from visible to FIR we re-estimate the redshifts and find z ~ 1.6–2.5. Owing to the strong spectral breaks observed in these galaxies, standard templates from the literature fail to reproduce the visible-to-near-IR part of the SEDs even when additional extinction is included. These sources strongly resemble dust-obscured galaxies selected in Spitzer observations with extreme visible-to-FIR colors, and the galaxy GN10 at z = 4. Galaxies with similar SEDs could contaminate other high-redshift surveys.

Star formation in a stellar mass-selected sample of galaxies to z=3 from the GOODS-NICMOS Survey

We present a study of the star-forming properties of a stellar mass-selected sample of galaxies in the GOODS (Great Observatories Origins Deep Survey) NICMOS Survey (GNS), based on deep Hubble Space Telescope (HST) imaging of the GOODS North and South fields. Using a stellar mass-selected sample, combined with HST/ACS and Spitzer data to measure both ultraviolet (UV) and infrared-derived star formation rates (SFRs), we investigate the star forming properties of a complete sample of ∼1300 galaxies down to log M_*= 9.5 at redshifts 1.5 < z < 3. Eight per cent of the sample is made up of massive galaxies with M_*≥ 10^11 M_⊙. We derive optical colours, dust extinctions and UV and infrared SFR to determine how the SFR changes as a function of both stellar mass and time. Our results show that SFR increases at higher stellar mass such that massive galaxies nearly double their stellar mass from star formation alone over the redshift range studied, but the average value of SFR for a given stellar mass remains constant over this ∼2 Gyr period. Furthermore, we find no strong evolution in the SFR for our sample as a function of mass over our redshift range of interest; in particular we do not find a decline in the SFR among massive galaxies, as is seen at z < 1. The most massive galaxies in our sample (log M_*≥ 11) have high average SFRs with values SFR_UV, corr= 103 ± 75 M_⊙ yr^−1, and yet exhibit red rest-frame (U−B) colours at all redshifts. We conclude that the majority of these red high-redshift massive galaxies are red due to dust extinction. We find that A_2800 increases with stellar mass, and show that between 45 and 85 per cent of massive galaxies harbour dusty star formation. These results show that even just a few Gyr after the first galaxies appear, there are strong relations between the global physical properties of galaxies, driven by stellar mass or another underlying feature of galaxies strongly related to the stellar mass.

Improving the identification of high-z Herschel sources with position priors and optical/NIR and FIR/mm photometric redshifts

We present preliminary results about the detection of high redshift (U)LIRGs in the Bullet cluster field by the PACS and SPIRE instruments within the Herschel Lensing Survey (HLS) Program. We describe in detail a photometric procedure designed to recover robust fluxes and deblend faint Herschel sources near the confusion noise. The method is based on the use of the positions of Spitzer/MIPS 24 μm sources as priors. Our catalogs are able to reliably (5σ) recover galaxies with fluxes above 6 and 10 mJy in the PACS 100 and 160 μm channels, respectively, and 12 to 18 mJy in the SPIRE bands. We also obtain spectral energy distributions covering the optical through the far-infrared/millimeter spectral ranges of all the Herschel detected sources, and analyze them to obtain independent estimations of the photometric redshift based on either stellar population or dust emission models. We exemplify the potential of the combined use of Spitzer position priors plus independent optical and IR photometric redshifts to robustly assign optical/NIR counterparts to the sources detected by Herschel and other (sub-)mm instruments.

The far-infrared/submillimeter properties of galaxies located behind the Bullet cluster

The Herschel Lensing Survey (HLS) takes advantage of gravitational lensing by massive galaxy clusters to sample a population of high-redshift galaxies which are too faint to be detected above the confusion limit of current far-infrared/submillimeter telescopes. Measurements from 100-500 μm bracket the peaks of the far-infrared spectral energy distributions of these galaxies, characterizing their infrared luminosities and star formation rates. We introduce initial results from our science demonstration phase observations, directed toward the Bullet cluster (1E0657-56). By combining our observations with LABOCA 870 μm and AzTEC 1.1 mm data we fully constrain the spectral energy distributions of 19 MIPS 24 μm-selected galaxies which are located behind the cluster. We find that their colors are best fit using templates based on local galaxies with systematically lower infrared luminosities. This suggests that our sources are not like local ultra-luminous infrared galaxies in which vigorous star formation is contained in a compact highly dust-obscured region. Instead, they appear to be scaled up versions of lower luminosity local galaxies with star formation occurring on larger physical scales.

Local luminous infrared galaxies: Spatially resolved mid-infrared observations with Spitzer/IRS

Luminous Infrared (IR) Galaxies (LIRGs, L_IR=10^11-10 L_⨀) are an important cosmological class of galaxies as they are the main contributors to the co-moving star formation rate density of the universe at z=1. In this paper we present a guaranteed time observation (GTO) Spitzer InfraRed Spectrograph (IRS) program aimed to obtain spectral mapping of a sample of 14 local d<76Mpc LIRGs. The data cubes map, at least, the central 20arcsec X 20arcsec to 30 arcsec X 30 arcsec regions of the galaxies, and use all four IRS modules covering the full 5-38 μ m spectral range. The final goal of this project is to characterize fully the mid-IR properties of local LIRGs as a first step to understanding their more distant counterparts. In this paper we present the first results of this GTO program. The IRS spectral mapping data allow us to build spectral maps of the bright mid-IR emission lines (e.g., [Ne II] 12.81 μ m, [Ne III]15.56 μ m, [S III] 18.71 μ m, H_2 at 17 μ m), continuum, the 6.2 and 11.3 μ m polycyclic aromatic hydrocarbon (PAH) features, and the 9.7 μ m silicate feature, as well as to extract 1D spectra for regions of interest in each galaxy. The IRS data are used to obtain spatially resolved measurements of the extinction using the 9.7 μ m silicate feature, and to trace star forming regions using the neon lines and the PAH features. We also investigate a number of active galactic nuclei (AGN) indicators, including the presence of high excitation emission lines and a strong dust continuum emission at around 6 9.7 μ m . We finally use the integrated Spitzer/IRS spectra as templates of local LIRGs. We discuss several possible uses for these templates, including the calibration of the star formation rate of IR-bright galaxies at high redshift. We also predict the intensities of the brightest mid-IR emission lines for LIRGs as a function of redshift, and compare them with the expected sensitivities of future space IR missions.

The host galaxies and black holes of typical z~0.5-1.4 AGNs

We study the stellar and star formation properties of the host galaxies of 58 X-ray-selected AGNs in the GOODS portion of the Chandra Deep Field South (CDF-S) region at z ~ 0.5-1.4. The AGNs are selected such that their rest-frame UV to near-infrared spectral energy distributions (SEDs) are dominated by stellar emission; i.e., they show a prominent 1.6 μm bump, thus minimizing the AGN emission "contamination." This AGN population comprises approximately 50% of the X-ray-selected AGNs at these redshifts. We find that AGNs reside in the most massive galaxies at the redshifts probed here. Their characteristic stellar masses (M_* ~ 7.8 × 10^10 and M_* ~ 1.2 × 10^11 M_☉ at median redshifts of 0.67 and 1.07, respectively) appear to be representative of the X-ray-selected AGN population at these redshifts and are intermediate between those of local type 2 AGNs and high-redshift (z ~ 2) AGNs. The inferred black hole masses (M_BH ~ 2 × 10^8 M_☉) of typical AGNs are similar to those of optically identified quasars at similar redshifts. Since the AGNs in our sample are much less luminous (L_2–10 keV < 10^44 erg s^−1) than quasars, typical AGNs have low Eddington ratios (η ~ 0.01-0.001). This suggests that, at least at intermediate redshifts, the cosmic AGN "downsizing" is due to both a decrease in the characteristic stellar mass of typical host galaxies and less efficient accretion. Finally, there is no strong evidence in AGN host galaxies for either highly suppressed star formation (expected if AGNs played a role in quenching star formation) or elevated star formation when compared to mass-selected (i.e., IRAC-selected) galaxies of similar stellar masses and redshifts.

The SCUBA HAlf degree extragalactic survey - III. Identification of radio and mid-infrared counterparts to submillimetre galaxies

Determining an accurate position for a submillimetre (submm) galaxy (SMG) is the crucial step that enables us to move from the basic properties of an SMG sample - source counts and 2D clustering - to an assessment of their detailed, multiwavelength properties, their contribution to the history of cosmic star formation and their links with present-day galaxy populations. In this paper, we identify robust radio and/or infrared (IR) counterparts, and hence accurate positions, for over two-thirds of the SCUBA HAlf-Degree Extragalactic Survey (SHADES) Source Catalogue, presenting optical, 24-μm and radio images of each SMG. Observed trends in identification rate have given no strong rationale for pruning the sample. Uncertainties in submm position are found to be consistent with theoretical expectations, with no evidence for significant additional sources of error. Employing the submm/radio redshift indicator, via a parametrization appropriate for radio-identified SMGs with spectroscopic redshifts, yields a median redshift of 2.8 for the radio-identified subset of SHADES, somewhat higher than the median spectroscopic redshift. We present a diagnostic colour-colour plot, exploiting Spitzer photometry, in which we identify regions commensurate with SMGs at very high redshift. Finally, we find that significantly more SMGs have multiple robust counterparts than would be expected by chance, indicative of physical associations. These multiple systems are most common amongst the brightest SMGs and are typically separated by 2-6 arcsec, similar to 15-20/sin i kpc at z~ 2, consistent with early bursts seen in merger simulations.

Emission-line galaxy surveys with GTC: The Hα Universe

Emission line galaxies are the most easily detected and studied objects in the high redshift Universe. They are being used to trace the evolution of critical observables of the Universe such as Star Formation Rate densities, starburst properties and abundances. Most of the research is being done using [OII]3727 and UV lines, but Hα is still one the best tracers for Star Formation Rate and physical properties of current star-forming galaxies. As a complementary contribution to studies of galaxy evolution, our team has focused into a long-term project to study the population of Hα-selected star-forming galaxies of the Universe at different redshifts. In 1995 we first determined the local Hα luminosity function, and from it the Star Formation Rate density (SFRd) of the local Universe. We then, using narrow-band imaging in the optical, extended this measurement to z ≃ 0.24 and z ≃ 0.4. Working in the near-infrared, GTC will be a very powerful tool to study the evolution of the Hα emission-line galaxy populations at different redshifts. We will both quantify the SFRd evolution and characterize the star-forming galaxy populations by directly comparing the same observables at all redshifts up to z ≃ 2.5.

Episodic star formation in a group of LAEs at z=5.07

We are undertaking a search for high-redshift low-luminosity Lyman Alpha sources in the SHARDS (Survey for High-z Absorption Red and Dead Sources) survey. Among the pre-selected Lyman Alpha sources two candidates were spotted, located 3.19 arcsec apart, and tentatively at the same redshift. Here, we report on the spectroscopic confirmation with Gran Telescopio Canarias of the Lyman Alpha emission from this pair of galaxies at a confirmed spectroscopic redshifts of z=5.07. Furthermore, one of the sources is interacting/merging with another close companion that looks distorted. Based on the analysis of the spectroscopy and additional photometric data, we infer that most of the stellar mass of these objects was assembled in a burst of star formation 100 Myr ago. A more recent burst (2 Myr old) is necessary to account for the measured Lyman Alpha flux. We claim that these two galaxies are good examples of Lyman Alpha sources undergoing episodic star formation. Besides, these sources very likely constitute a group of interacting Lyman Alpha emitters (LAEs).

[C II] and 12^CO(1-0) emission maps in HLSJ091828.6+514223: a strongly lensed interacting system AT z = 5.24

We present Submillimeter Array [C II] 158 μm and Karl G. Jansky Very Large Array 12^CO(1-0) line emission maps for the bright, lensed, submillimeter source at z = 5.2430 behind A 773: HLSJ091828.6+514223 (HLS0918). We combine these measurements with previously reported line profiles, including multiple 12^CO rotational transitions, [C I], water, and [N II], providing some of the best constraints on the properties of the interstellar medium in a galaxy at z > 5. HLS0918 has a total far-infrared (FIR) luminosity L_FIR(8–1000 μm) = (1.6 ± 0.1) × 10^14 L_☉ μ^–1, where the total magnification μ_total = 8.9 ± 1.9, via a new lens model from the [C II] and continuum maps. Despite a HyLIRG luminosity, the FIR continuum shape resembles that of a local LIRG. We simultaneously fit all of the observed spectral line profiles, finding four components that correspond cleanly to discrete spatial structures identified in the maps. The two most redshifted spectral components occupy the nucleus of a massive galaxy, with a source-plane separation <1 kpc. The reddest dominates the continuum map (demagnified L_FIR, component = (1.1 ± 0.2) × 10^13 L_☉) and excites strong water emission in both nuclear components via a powerful FIR radiation field from the intense star formation. A third star-forming component is most likely a region of a merging companion (ΔV ~ 500 km s^–1) exhibiting generally similar gas properties. The bluest component originates from a spatially distinct region and photodissociation region analysis suggests that it is lower density, cooler, and forming stars less vigorously than the other components. Strikingly, it has very strong [N II] emission, which may suggest an ionized, molecular outflow. This comprehensive view of gas properties and morphology in HLS0918 previews the science possible for a large sample of high-redshift galaxies once ALMA attains full sensitivity.

The Herschel Lensing Survey (HLS): Overview

The Herschel Lensing Survey (HLS) will conduct deep PACS and SPIRE imaging of ~40 massive clusters of galaxies. The strong gravitational lensing power of these clusters will enable us to penetrate through the confusion noise, which sets the ultimate limit on our ability to probe the Universe with Herschel. Here we present an overview of our survey and a summary of the major results from our science demonstration phase (SDP) observations of the Bullet cluster (z = 0.297). The SDP data are rich and allow us to study not only the background high-redshift galaxies (e. g., strongly lensed and distorted galaxies at z = 2.8 and 3.2) but also the properties of cluster-member galaxies. Our preliminary analysis shows a great diversity of far-infrared/submillimeter spectral energy distributions (SEDs), indicating that we have much to learn with Herschel about the properties of galaxy SEDs. We have also detected the Sunyaev-Zel'dovich (SZ) effect increment with the SPIRE data. The success of this SDP program demonstrates the great potential of the Herschel Lensing Survey to produce exciting results in a variety of science areas.

A bright z=5.2 lensed submillimeter galaxy in the field of Abell 773 HLSJ091828.6+514223

During our Herschel Lensing Survey (HLS) of massive galaxy clusters, we have discovered an exceptionally bright source behind the z = 0.22 cluster Abell 773, which appears to be a strongly lensed submillimeter galaxy (SMG) at z = 5.2429. This source is unusual compared to most other lensed sources discovered by Herschel so far, because of its higher submm flux (∼200 mJy at 500 μm) and its high redshift. The dominant lens is a foreground z = 0.63 galaxy, not the cluster itself. The source has a far-infrared (FIR) luminosity of L_FIR = 1.1 × 10^14/μ L_⨀, where μ is the magnification factor, likely ∼11. We report here the redshift identification through CO lines with the IRAM-30 m, and the analysis of the gas excitation, based on CO(7–6), CO(6–5), CO(5–4) detected at IRAM and the CO(2–1) at the EVLA. All lines decompose into a wide and strong red component, and a narrower and weaker blue component, 540 km s^−1 apart. Assuming the ultraluminous galaxy (ULIRG) CO-to-H_2 conversion ratio, the H_2 mass is 5.8×10^11/μ M_⨀, of which one third is in a cool component. From the CI(^3P_2−^3 P_1) line we derive a C_I/H_2 number abundance of 6 × 10^−5 similar to that in other ULIRGs. The H_2O_p(2, 0, 2−1, 1, 1) line is strong only in the red velocity component, with an intensity ratio I(H_2O)/I(CO) ∼ 0.5, suggesting a strong local FIR radiation field, possibly from an active nucleus (AGN) component. We detect the [NII]205 μm line for the first time at high-z. It shows comparable blue and red components, with a strikingly broad blue one, suggesting strong ionized gas flows.

Investigating evidence for different black hole accretion modes since redshift z ∼ 1

Chandra data in the COSMOS, AEGIS-XD and 4 Ms Chandra Deep Field South are combined with multiwavelength photometry available in those fields to determine the rest-frame U − V versus V − J colours of X-ray AGN hosts in the redshift intervals 0.1 < z < 0.6 (mean z¯=0.40) and 0.6 < z < 1.2 (mean z¯=0.85). This combination of colours provides an effective and least model-dependent means of separating quiescent from star-forming, including dust reddened, galaxies. Morphological information emphasizes differences between AGN populations split by their U − V versus V − J colours. AGN in quiescent galaxies consist almost exclusively of bulges, while star-forming hosts are equally split between early- and late-type hosts. The position of AGN hosts on the U − V versusV − J diagram is then used to set limits on the accretion density of the Universe associated with evolved and star-forming systems independent of dust induced biases. It is found that most of the black hole growth at z≈ 0.40 and 0.85 is associated with star-forming hosts. Nevertheless, a non-negligible fraction of the X-ray luminosity density, about 15–20 per cent, at both z¯=0.40 and 0.85, is taking place in galaxies in the quiescent region of the U − V versus V − J diagram. For the low-redshift sub-sample, 0.1 < z < 0.6, we also find tentative evidence, significant at the 2σ level, that AGN split by their U − V and V − J colours have different Eddington ratio distributions. AGN in blue star-forming hosts dominate at relatively high Eddington ratios. In contrast, AGN in red quiescent hosts become increasingly important as a fraction of the total population towards low Eddington ratios. At higher redshift, z > 0.6, such differences are significant at the 2σ level only for sources with Eddington ratios ≳ 10^− 3. These findings are consistent with scenarios in which diverse accretion modes are responsible for the build-up of supermassive black holes at the centres of galaxies. We compare these results with the predictions of theGALFORM semi-analytic model for the cosmological evolution of AGN and galaxies. This model postulates two black hole fuelling modes, the first is linked to star formation events and the second takes place in passive galaxies. GALFORM predicts that a substantial fraction of the black hole growth at z < 1 is associated with quiescent galaxies, in apparent conflict with the observations. Relaxing the strong assumption of the model that passive AGN hosts have zero star formation rate could bring those predictions in better agreement with the data.

The XMM deep survey in the CDF-S. VI. Obscured AGN selected as infrared power-law galaxies

Context. Accretion onto supermassive black holes is believed to occur mostly in obscured active galactic nuclei (AGN). Such objects are proving rather elusive in surveys of distant galaxies, including those at X-ray energies. Aims. Our main goal is to determine whether the revised IRAC criteria of Donley et al. (2012, ApJ, 748, 142; objects with an infrared (IR) power-law spectral shape), are effective at selecting X-ray type-2 AGN (i.e., absorbed N_H > 10^22 cm^-2). Methods. We present the results from the X-ray spectral analysis of 147 AGN selected by cross-correlating the highest spectral quality ultra-deep XMM-Newton and the Spitzer/IRAC catalogues in the Chandra Deep Field South. Consequently it is biased towards sources with high S/N X-ray spectra. In order to measure the amount of intrinsic absorption in these sources, we adopt a simple X-ray spectral model that includes a power-law modified by intrinsic absorption at the redshift of each source and a possible soft X-ray component. Results. We find 21/147 sources to be heavily absorbed but the uncertainties in their obscuring column densities do not allow us to confirm their Compton-Thick nature without resorting to additional criteria. Although IR power-law galaxies are less numerous in our sample than IR non-power-law galaxies (60 versus 87 respectively), we find that the fraction of absorbed (N_H^intr > 10^22 cm^-2) AGN is significantly higher (at about 3 sigma level) for IR-power-law sources (similar to 2/3) than for those sources that do not meet this IR selection criteria (~1/2). This behaviour is particularly notable at low luminosities, but it appears to be present, although with a marginal significance, at all luminosities. Conclusions. We therefore conclude that the IR power-law method is efficient in finding X-ray-absorbed sources. We would then expect that the long-sought dominant population of absorbed AGN is abundant among IR power-law spectral shape sources not detected in X-rays.