CANDELS: The progenitors of compact quiescent galaxies at z ~ 2

We combine high-resolution Hubble Space Telescope/WFC3 images with multi-wavelength photometry to track the evolution of structure and activity of massive (M_*> 10^10 M_☉) galaxies at redshifts z = 1.4-3 in two fields of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey. We detect compact, star-forming galaxies (cSFGs) whose number densities, masses, sizes, and star formation rates (SFRs) qualify them as likely progenitors of compact, quiescent, massive galaxies (cQGs) at z = 1.5-3. At z≲2, cSFGs present SFR = 100-200 M_☉ yr^–1, yet their specific star formation rates (sSFR ~ 10^–9 yr^–1) are typically half that of other massive SFGs at the same epoch, and host X-ray luminous active galactic nuclei (AGNs) 30 times (~30%) more frequently. These properties suggest that cSFGs are formed by gas-rich processes (mergers or disk-instabilities) that induce a compact starburst and feed an AGN, which, in turn, quench the star formation on dynamical timescales (few 10^8 yr). The cSFGs are continuously being formed at z = 2-3 and fade to cQGs down to z ~ 1.5. After this epoch, cSFGs are rare, thereby truncating the formation of new cQGs. Meanwhile, down to z = 1, existing cQGs continue to enlarge to match local QGs in size, while less-gas-rich mergers and other secular mechanisms shepherd (larger) SFGs as later arrivals to the red sequence. In summary, we propose two evolutionary tracks of QG formation: an early (z≲2), formation path of rapidly quenched cSFGs fading into cQGs that later enlarge within the quiescent phase, and a late-arrival (z≳2) path in which larger SFGs form extended QGs without passing through a compact state.

Testing diagnostics of nuclear activity and star formation in galaxies at z > 1

We present some of the first science data with the new Keck/MOSFIRE instrument to test the effectiveness of different AGN/SF diagnostics at z ~ 1.5. MOSFIRE spectra were obtained in three H-band multi-slit masks in the GOODS-S field, resulting in 2 hr exposures of 36 emission-line galaxies. We compare X-ray data with the traditional emission-line ratio diagnostics and the alternative mass-excitation and color-excitation diagrams, combining new MOSFIRE infrared data with previous HST/WFC3 infrared spectra (from the 3D-HST survey) and multiwavelength photometry. We demonstrate that a high [O III]/Hβ ratio is insufficient as an active galactic nucleus (AGN) indicator at z > 1. For the four X-ray-detected galaxies, the classic diagnostics ([O III]/Hβ versus [N II]/Hα and [S II]/Hα) remain consistent with X-ray AGN/SF classification. The X-ray data also suggest that "composite" galaxies (with intermediate AGN/SF classification) host bona fide AGNs. Nearly ~2/3 of the z ~ 1.5 emission-line galaxies have nuclear activity detected by either X-rays or the classic diagnostics. Compared to the X-ray and line ratio classifications, the mass-excitation method remains effective at z > 1, but we show that the color-excitation method requires a new calibration to successfully identify AGNs at these redshifts.

The relation between cool cluster cores and Herschel-detected star formation in brightest cluster galaxies

We present far-infrared (FIR) analysis of 68 brightest cluster galaxies (BCGs) at 0.08 < z < 1.0. Deriving total infrared luminosities directly from Spitzer and Herschel photometry spanning the peak of the dust component (24-500 μm), we calculate the obscured star formation rate (SFR). 22^+6.2 _–5.3% of the BCGs are detected in the far-infrared, with SFR = 1-150 M ☉ yr^–1. The infrared luminosity is highly correlated with cluster X-ray gas cooling times for cool-core clusters (gas cooling time <1 Gyr), strongly suggesting that the star formation in these BCGs is influenced by the cluster-scale cooling process. The occurrence of the molecular gas tracing Hα emission is also correlated with obscured star formation. For all but the most luminous BCGs (L_TIR > 2 × 10^11 L_☉), only a small (≤0.4 mag) reddening correction is required for SFR(Hα) to agree with SFR_FIR. The relatively low Hα extinction (dust obscuration), compared to values reported for the general star-forming population, lends further weight to an alternate (external) origin for the cold gas. Finally, we use a stacking analysis of non-cool-core clusters to show that the majority of the fuel for star formation in the FIR-bright BCGs is unlikely to originate from normal stellar mass loss.

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.

Determining star formation rates for infrared galaxies

We show that measures of star formation rates (SFRs) for infrared galaxies using either single-band 24 μm or extinction-corrected Paα luminosities are consistent in the total infrared luminosity = L(TIR) ~ 10^10 L_☉ range. MIPS 24 μm photometry can yield SFRs accurately from this luminosity upward: SFR(M_☉ yr^–1) = 7.8 × 10^–10 L(24 μm, L_☉) from L(TIR) = 5× 10^9 L_☉ to 10^11 L_☉ and SFR = 7.8 × 10^–10 L(24 μm, L_☉)(7.76 × 10^–11 L(24))^0.048 for higher L(TIR). For galaxies with L(TIR) ≥ 10^10 L_☉, these new expressions should provide SFRs to within 0.2 dex. For L(TIR) ≥ 10^11 L_☉, we find that the SFR of infrared galaxies is significantly underestimated using extinction-corrected Paα (and presumably using any other optical or near-infrared recombination lines). As a part of this work, we constructed spectral energy distribution templates for eleven luminous and ultraluminous purely star forming infrared galaxies and over the spectral range 0.4 μm to 30 cm. We use these templates and the SINGS data to construct average templates from 5 μm to 30 cm for infrared galaxies with L(TIR) = 5× 10^9 to 10^13 L_☉. All of these templates are made available online.

The star formation and extinction coevolution of UV-selected galaxies over 0.05 < z < 1.2

We use a new stacking technique to obtain mean mid-IR and far-IR to far-UV flux ratios over the rest-frame near-UV, near-IR color-magnitude diagram. We employ COMBO-17 redshifts and COMBO-17 optical, GALEX far- and near-UV, and Spitzer IRAC and MIPS mid-IR photometry. This technique permits us to probe the infrared excess (IRX), the ratio of far-IR to far-UV luminosity, and the specific star formation rate (SSFR) and their coevolution over 2 orders of magnitude of stellar mass and over redshift 0.1 < z < 1.2. We find that the SSFR and the characteristic mass (Script M_0) above which the SSFR drops increase with redshift (downsizing). At any given epoch, the IRX is an increasing function of mass up to Script M_0. Above this mass the IRX falls, suggesting gas exhaustion. In a given mass bin below Script M_0, the IRX increases with time in a fashion consistent with enrichment. We interpret these trends using a simple model with a Schmidt-Kennicutt law and extinction that tracks gas density and enrichment. We find that the average IRX and SSFR follow a galaxy age parameter ξ, which is determined mainly by the galaxy mass and time since formation. We conclude that blue-sequence galaxies have properties which show simple, systematic trends with mass and time such as the steady buildup of heavy elements in the interstellar media of evolving galaxies and the exhaustion of gas in galaxies that are evolving off the blue sequence. The IRX represents a tool for selecting galaxies at various stages of evolution.

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.

Sub-kiloparsec alma imaging of compact star-forming galaxies at z ~ 2.5: revealing the formation of dense galactic cores in the progenitors of compact quiescent galaxies

We present spatially resolved Atacama Large Millimeter/submillimeter Array (ALMA) 870 μm dust continuum maps of six massive, compact, dusty star-forming galaxies at z ~ 2.5. These galaxies are selected for their small rest-frame optical sizes (r_e,F160W ~ 1.6 kpc) and high stellar mass densities that suggest that they are direct progenitors of compact quiescent galaxies at z ~ 2. The deep observations yield high far-infrared (FIR) luminosities of L_IR = 10^12.3-12.8 L_⨀ and star formation rates (SFRs) of SFR = 200–700 M_⊙ yr^−1, consistent with those of typical star-forming "main sequence" galaxies. The high spatial resolution (FWHM ~ 0 12–0 18) ALMA and Hubble Space Telescope photometry are combined to construct deconvolved, mean radial profiles of their stellar mass and (UV+IR) SFR. We find that the dusty, nuclear IR–SFR overwhelmingly dominates the bolometric SFR up to r ~ 5 kpc, by a factor of over 100× from the unobscured UV–SFR. Furthermore, the effective radius of the mean SFR profile (r_e,SFR ~ 1 kpc) is ~30% smaller than that of the stellar mass profile. The implied structural evolution, if such nuclear starburst last for the estimated gas depletion time of Δt = ±100 Myr, is a 4×increase of the stellar mass density within the central 1 kpc and a 1.6× decrease of the half-mass–radius. This structural evolution fully supports dissipation-driven, formation scenarios in which strong nuclear starbursts transform larger, star-forming progenitors into compact quiescent galaxies.

Space density distribution of galaxies in the absolute magnitude - rotation velocity plane: a volume-complete Tully-Fisher relation from CALIFA stellar kinematics

We measured the distribution in absolute magnitude - circular velocity space for a well-defined sample of 199 rotating galaxies of the Calar Alto Legacy Integral Field Area Survey (CALIFA) using their stellar kinematics. Our aim in this analysis is to avoid subjective selection criteria and to take volume and large-scale structure factors into account. Using stellar velocity fields instead of gas emission line kinematics allows including rapidly rotating early-type galaxies. Our initial sample contains 277 galaxies with available stellar velocity fields and growth curve r-band photometry. After rejecting 51 velocity fields that could not be modelled because of the low number of bins, foreground contamination, or significant interaction, we performed Markov chain Monte Carlo modelling of the velocity fields, from which we obtained the rotation curve and kinematic parameters and their realistic uncertainties. We performed an extinction correction and calculated the circular velocity v_circ accounting for the pressure support of a given galaxy. The resulting galaxy distribution on the M-r - v(circ) plane was then modelled as a mixture of two distinct populations, allowing robust and reproducible rejection of outliers, a significant fraction of which are slow rotators. The selection effects are understood well enough that we were able to correct for the incompleteness of the sample. The 199 galaxies were weighted by volume and large-scale structure factors, which enabled us to fit a volume-corrected Tully-Fisher relation (TFR). More importantly, we also provide the volume-corrected distribution of galaxies in the M_r - v_circ plane, which can be compared with cosmological simulations. The joint distribution of the luminosity and circular velocity space densities, representative over the range of -20 > M_r > -22 mag, can place more stringent constraints on the galaxy formation and evolution scenarios than linear TFR fit parameters or the luminosity function alone.

Stellar populations in bulges of spiral galaxies

We present a few results of a work in progress tackling the radial spectroscopic properties of bulges of spiral galaxies.

Limits to dark matter annihilation cross-section from a combined analysis of MAGIC and Fermi-LAT observations of dwarf satellite galaxies

We present the first joint analysis of gamma-ray data from the MAGIC Cherenkov telescopes and the Fermi Large Area Telescope (LAT) to search for gamma-ray signals from dark matter annihilation in dwarf satellite galaxies. We combine 158 hours of Segue 1 observations with MAGIC with 6-year observations of 15 dwarf satellite galaxies by the Fermi-LAT. We obtain limits on the annihilation cross-section for dark matter particle masses between 10 GeV and 100 TeV – the widest mass range ever explored by a single gamma-ray analysis. These limits improve on previously published Fermi-LAT and MAGIC results by up to a factor of two at certain masses. Our new inclusive analysis approach is completely generic and can be used to perform a global, sensitivity-optimized dark matter search by combining data from present and future gamma-ray and neutrino detectors.

The star formation history of CALIFA galaxies: Radial structures

We have studied the radial structure of the stellar mass surface density (μ∗) and stellar population age as a function of the total stellar mass and morphology for a sample of 107 galaxies from the CALIFA survey. We applied the fossil record method based on spectral synthesis techniques to recover the star formation history (SFH), resolved in space and time, in spheroidal and disk dominated galaxies with masses from 10^9 to 10^12 M_⊙. We derived the half-mass radius, and we found that galaxies are on average 15% more compact in mass than in light. The ratio of half-mass radius to half-light radius (HLR) shows a dual dependence with galaxy stellar mass; it decreases with increasing mass for disk galaxies, but is almost constant in spheroidal galaxies. In terms of integrated versus spatially resolved properties, we find that the galaxy-averaged stellar population age, stellar extinction, and μ_∗ are well represented by their values at 1 HLR. Negative radial gradients of the stellar population ages are present in most of the galaxies, supporting an inside-out formation. The larger inner (≤1 HLR) age gradients occur in the most massive (10^11 M_⊙) disk galaxies that have the most prominent bulges; shallower age gradients are obtained in spheroids of similar mass. Disk and spheroidal galaxies show negative μ∗ gradients that steepen with stellar mass. In spheroidal galaxies, μ∗ saturates at a critical value (~7 × 10^2 M_⊙/pc^2 at 1 HLR) that is independent of the galaxy mass. Thus, all the massive spheroidal galaxies have similar local μ_∗ at the same distance (in HLR units) from the nucleus. The SFH of the regions beyond 1 HLR are well correlated with their local μ_∗, and follow the same relation as the galaxy-averaged age and μ_∗; this suggests that local stellar mass surface density preserves the SFH of disks. The SFH of bulges are, however, more fundamentally related to the total stellar mass, since the radial structure of the stellar age changes with galaxy mass even though all the spheroid dominated galaxies have similar radial structure in μ_∗. Thus, galaxy mass is a more fundamental property in spheroidal systems, while the local stellar mass surface density is more important in disks.

The nature of LINER galaxies: Ubiquitous hot old stars and rare accreting black holes

Context. Galaxies, which often contain ionised gas, sometimes also exhibit a so-called low-ionisation nuclear emission line region (LINER). For 30 years, this was attributed to a central mass-accreting supermassive black hole (more commonly known as active galactic nucleus, AGN) of low luminosity, making LINER galaxies the largest AGN sub-population, which dominate in numbers over higher luminosity Seyfert galaxies and quasars. This, however, poses a serious problem. While the inferred energy balance is plausible, many LINERs clearly do not contain any other independent signatures of an AGN. Aims. Using integral field spectroscopic data from the CALIFA survey, we compare the observed radial surface brightness profiles with what is expected from illumination by an AGN. Methods. Essential for this analysis is a proper extraction of emission lines, especially weak lines, such as Balmer H beta lines, which are superposed on an absorption trough. To accomplish this, we use the GANDALF code, which simultaneously fits the underlying stellar continuum and emission lines. Results. For 48 galaxies with LINER-like emission, we show that the radial emission-line surface brightness profiles are inconsistent with ionisation by a central point-source and hence cannot be due to an AGN alone. Conclusions. The most probable explanation for the excess LINER-like emission is ionisation by evolved stars during the short but very hot and energetic phase known as post-AGB. This leads us to an entirely new interpretation. Post-AGB stars are ubiquitous and their ionising effect should be potentially observable in every galaxy with the gas present and with stars older than ~1 Gyr unless a stronger radiation field from young hot stars or an AGN outshines them. This means that galaxies with LINER-like emission are not a class defined by a property but rather by the absence of a property. It also explains why LINER emission is observed mostly in massive galaxies with old stars and little star formation.

The ionized gas in the CALIFA early-type galaxies I. Mapping two representative cases: NGC 6762 and NGC 5966

As part of the ongoing CALIFA survey, we have conducted a thorough bidimensional analysis of the ionized gas in two E/S0 galaxies, NGC 6762 and NGC 5966, aiming to shed light on the nature of their warm ionized ISM. Specifically, we present optical (3745–7300 Å) integral field spectroscopy obtained with the PMAS/PPAK integral field spectrophotometer. Its wide field-of-view (1′ × 1′) covers the entire optical extent of each galaxy down to faint continuum surface brightnesses. To recover the nebular lines, we modeled and subtracted the underlying stellar continuum from the observed spectra using the STARLIGHT spectral synthesis code. The pure emission-line spectra were used to investigate the gas properties and determine the possible sources of ionization. We show the advantages of IFU data in interpreting the complex nature of the ionized gas in NGC 6762 and NGC 5966. In NGC 6762, the ionized gas and stellar emission display similar morphologies, while the emission line morphology is elongated in NGC 5966, spanning ~6 kpc, and is oriented roughly orthogonal to the major axis of the stellar continuum ellipsoid. Whereas gas and stars are kinematically aligned in NGC 6762, the gas is kinematically decoupled from the stars in NGC 5966. A decoupled rotating disk or an “ionization cone” are two possible interpretations of the elongated ionized gas structure in NGC 5966. The latter would be the first “ionization cone” of such a dimension detected within a weak emission-line galaxy. Both galaxies have weak emission-lines relative to the continuum[EW(Hα) ≲ 3 Å] and have very low excitation, log([OIII]λ5007/Hβ) ≲ 0.5. Based on optical diagnostic ratios ([OIII]λ5007/Hβ, [NII]λ6584/Hα, [SII]λ6717, 6731/Hα, [OI]λ6300/Hα), both objects contain a LINER nucleus and an extended LINER-like gas emission. The emission line ratios do not vary significantly with radius or aperture, which indicates that the nebular properties are spatially homogeneous. The gas emission in NGC 6762 can be best explained by photoionization by pAGB stars without the need of invoking any other excitation mechanism. In the case of NGC 5966, the presence of a nuclear ionizing source seems to be required to shape the elongated gas emission feature in the “ionization cone” scenario, although ionization by pAGB stars cannot be ruled out. Further study of this object is needed to clarify the nature of its elongated gas structure.

Pathways to quiescence: SHARDS view on the star formation histories of massive quiescent galaxies at 1.0 < z < 1.5

We present star formation histories (SFHs) for a sample of 104 massive (stellar mass M > 10^10 M_⊙) quiescent galaxies (MQGs) at z = 1.0–1.5 from the analysis of spectrophotometric data from the Survey for High-z Absorption Red and Dead Sources (SHARDS) and HST/WFC3 G102 and G141 surveys of the GOODS-North field, jointly with broad-band observations from ultraviolet (UV) to far-infrared (far-IR). The sample is constructed on the basis of rest-frame UVJ colours and specific star formation rates (sSFRs = SFR/Mass). The spectral energy distributions (SEDs) of each galaxy are compared to models assuming a delayed exponentially declining SFH. A Monte Carlo algorithm characterizes the degeneracies, which we are able to break taking advantage of the SHARDS data resolution, by measuring indices such as MgUV and D4000. The population of MQGs shows a duality in their properties. The sample is dominated (85 per cent) by galaxies with young mass-weighted ages, t_M t_M < 2 Gyr, short star formation time-scales, 〈τ〉 ∼ 60–200 Myr, and masses log(M/M_⊙) ∼ 10.5. There is an older population (15 per cent) with t_M t_M = 2–4 Gyr, longer star formation time-scales, 〈τ〉∼ 400 Myr, and larger masses, log(M/M_⊙) ∼ 10.7. The SFHs of our MQGs are consistent with the slope and the location of the main sequence of star-forming galaxies at z > 1.0, when our galaxies were 0.5–1.0 Gyr old. According to these SFHs, all the MQGs experienced a luminous infrared galaxy phase that lasts for ∼500 Myr, and half of them an ultraluminous infrared galaxy phase for ∼100 Myr. We find that the MQG population is almost assembled at z ∼ 1, and continues evolving passively with few additions to the population.