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 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.

Mid-infrared spectroscopy of lensed galaxies at 1 < z < 3: The nature of sources near the MIPS confusion limit

We present Spitzer IRS mid-infrared spectra for 15 gravitationally lensed, 24 μm-selected galaxies, and combine the results with four additional very faint galaxies with IRS spectra in the literature. The median intrinsic 24 μm flux density of the sample is 130 μJy, enabling a systematic survey of the spectral properties of the very faint 24 μm sources that dominate the number counts of Spitzer cosmological surveys. Six of the 19 galaxy spectra (32%) show the strong mid-IR continuua expected of AGNs; X-ray detections confirm the presence of AGNs in three of these cases, and reveal AGNs in two other galaxies. These results suggest that nuclear accretion may contribute more flux to faint 24 μm-selected samples than previously assumed. Almost all the spectra show some aromatic (PAH) emission features; the measured aromatic flux ratios do not show evolution from z = 0. In particular, the high signal-to-noise mid-IR spectrum of SMM J163554.2+661225 agrees remarkably well with low-redshift, lower luminosity templates. We compare the rest-frame 8 μm and total infrared luminosities of star-forming galaxies, and find that the behavior of this ratio with total IR luminosity has evolved modestly from z = 2 to z = 0. Since the high aromatic-to-continuum flux ratios in these galaxies rule out a dominant contribution by AGNs, this finding implies systematic evolution in the structure and/or metallicity of infrared sources with redshift. It also has implications for the estimates of star-forming rates inferred from 24 μm measurements, in the sense that at z ~ 2, a given observed frame 24 μm luminosity corresponds to a lower bolometric luminosity than would be inferred from low-redshift templates of similar luminosity at the corresponding rest wavelength.

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.

Are long gamma-ray bursts biased tracers of star formation? Clues from the host galaxies of the Swift/BAT6 complete sample of bright LGRBs II. Star formation rates and metallicities at z < 1

Aims. Long gamma-ray bursts (LGRBs) are associated with the deaths of massive stars and might therefore be a potentially powerful tool for tracing cosmic star formation. However, especially at low redshifts (z< 1.5) LGRBs seem to prefer particular types of environment. Our aim is to study the host galaxies of a complete sample of bright LGRBs to investigate the effect of the environment on GRB formation. Methods. We studied host galaxy spectra of the Swift/BAT6 complete sample of 14 z< 1 bright LGRBs. We used the detected nebular emission lines to measure the dust extinction, star formation rate (SFR), and nebular metallicity (Z) of the hosts and supplemented the data set with previously measured stellar masses M_*. The distributions of the obtained properties and their interrelations (e.g. mass-metallicity and SFR-M_* relations) are compared to samples of field star-forming galaxies. Results. We find that LGRB hosts at z< 1 have on average lower SFRs than if they were direct star formation tracers. By directly comparing metallicity distributions of LGRB hosts and star-forming galaxies, we find a good match between the two populations up to 12 +log (O/H)~8.4−8.5, after which the paucity of metal-rich LGRB hosts becomes apparent. The LGRB host galaxies of our complete sample are consistent with the mass-metallicity relation at similar mean redshift and stellar masses. The cutoff against high metallicities (and high masses) can explain the low SFR values of LGRB hosts. We find a hint of an increased incidence of starburst galaxies in the Swift/BAT6 z< 1 sample with respect to that of a field star-forming population. Given that the SFRs are low on average, the latter is ascribed to low stellar masses. Nevertheless, the limits on the completeness and metallicity availability of current surveys, coupled with the limited number of LGRB host galaxies, prevents us from investigating more quantitatively whether the starburst incidence is such as expected after taking into account the high-metallicity aversion of LGRB host galaxies.

Aperture effects on the oxygen abundance determinations from Califa data

This paper aims to provide aperture corrections for emission lines in a sample of spiral galaxies from the Calar Alto Legacy Integral Field Area Survey (CALIFA) database. In particular, we explore the behavior of the log([O III] λ5007/Hβ)/([N II] λ6583/Hα) (O3N2) and log[N II] lambda 6583/Hα (N2) flux ratios since they are closely connected to different empirical calibrations of the oxygen abundances in star-forming galaxies. We compute the median growth curves of Hα, Hα/Hβ, O3N2, and N-2 up to 2.5R(50) and 1.5 disk R-eff. These distances cover most of the optical spatial extent of the CALIFA galaxies. The growth curves simulate the effect of observing galaxies through apertures of varying radii. We split these growth curves by morphological types and stellar masses to check if there is any dependence on these properties. The median growth curve of the Hα flux shows a monotonous increase with radius with no strong dependence on galaxy inclination, morphological type, and stellar mass. The median growth curve of the Hα/HβH ratio monotonically decreases from the center toward larger radii, showing for small apertures a maximum value of ≈10% larger than the integrated one. It does not show any dependence on inclination, morphological type, and stellar mass. The median growth curve of N-2 shows a similar behavior, decreasing from the center toward larger radii. No strong dependence is seen on the inclination, morphological type, and stellar mass. Finally, the median growth curve of O3N2 increases monotonically with radius, and it does not show dependence on the inclination. However, at small radii it shows systematically higher values for galaxies of earlier morphological types and for high stellar mass galaxies. Applying our aperture corrections to a sample of galaxies from the SDSS survey at 0.02 ≤ z ≤ 0.3 shows that the average difference between fiber-based and aperture-corrected oxygen abundances, for different galaxy stellar mass and redshift ranges, reaches typically to ≈11%, depending on the abundance calibration used. This average difference is found to be systematically biased, though still within the typical uncertainties of oxygen abundances derived from empirical calibrations. Caution must be exercised when using observations of galaxies for small radii (e.g., below 0.5 R_eff) given the high dispersion shown around the median growth curves. Thus, the application of these median aperture corrections to derive abundances for individual galaxies is not recommended when their fluxes come from radii much smaller than either R_50 or R_eff.

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.

Mass assembly and morphological transformations since z ∼ 3 from CANDELS

We quantify the evolution of the stellar mass functions (SMFs) of star-forming and quiescent galaxies as a function of morphology from z ∼ 3 to the present. Our sample consists of ∼50 000 galaxies in the CANDELS fields (∼880 arcmin^2), which we divide into four main morphological types, i.e. pure bulge-dominated systems, pure spiral disc-dominated, intermediate two-component bulge+disc systems and irregular disturbed galaxies. At z ∼ 2, 80 per cent of the stellar mass density of star-forming galaxies is in irregular systems. However, by z ∼ 0.5, irregular objects only dominate at stellar masses below 10^9 M_⊙. A majority of the star-forming irregulars present at z ∼ 2 undergo a gradual transformation from disturbed to normal spiral disc morphologies by z ∼ 1 without significant interruption to their star formation. Rejuvenation after a quenching event does not seem to be common except perhaps for the most massive objects, because the fraction of bulge-dominated star-forming galaxies with M^*/M_⊙ > 10^10.7 reaches 40 per cent at z < 1. Quenching implies the presence of a bulge: the abundance of massive red discs is negligible at all redshifts over 2 dex in stellar mass. However, the dominant quenching mechanism evolves. At z > 2, the SMF of quiescent galaxies above M^* is dominated by compact spheroids. Quenching at this early epoch destroys the disc and produces a compact remnant unless the star-forming progenitors at even higher redshifts are significantly more dense. At 1 < z < 2, the majority of newly quenched galaxies are discs with a significant central bulge. This suggests that mass quenching at this epoch starts from the inner parts and preserves the disc. At z < 1, the high-mass end of the passive SMF is globally in place and the evolution mostly happens at stellar masses below 10^10 M_⊙. These low-mass galaxies are compact, bulge-dominated systems, which were environmentally quenched: destruction of the disc through ram-pressure stripping is the likely process.

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.

Mass-metallicity relation explored with CALIFA I. Is there a dependence on the star-formation rate?

We studied the global and local ℳ-Z relation based on the first data available from the CALIFA survey (150 galaxies). This survey provides integral field spectroscopy of the complete optical extent of each galaxy (up to 2−3 effective radii), with a resolution high enough to separate individual H II regions and/or aggregations. About 3000 individual H II regions have been detected. The spectra cover the wavelength range between [OII]3727 and [SII]6731, with a sufficient signal-to-noise ratio to derive the oxygen abundance and star-formation rate associated with each region. In addition, we computed the integrated and spatially resolved stellar masses (and surface densities) based on SDSS photometric data. We explore the relations between the stellar mass, oxygen abundance and star-formation rate using this dataset. We derive a tight relation between the integrated stellar mass and the gas-phase abundance, with a dispersion lower than the one already reported in the literature (σ_Δlog (O/H) = 0.07 dex). Indeed, this dispersion is only slightly higher than the typical error derived for our oxygen abundances. However, we found no secondary relation with the star-formation rate other than the one induced by the primary relation of this quantity with the stellar mass. The analysis for our sample of ~3000 individual H II   regions confirms (i) a local mass-metallicity relation and (ii) the lack of a secondary relation with the star-formation rate. The same analysis was performed with similar results for the specific star-formation rate. Our results agree with the scenario in which gas recycling in galaxies, both locally and globally, is much faster than other typical timescales, such like that of gas accretion by inflow and/or metal loss due to outflows. In essence, late-type/disk-dominated galaxies seem to be in a quasi-steady situation, with a behavior similar to the one expected from an instantaneous recycling/closed-box model.

DUst around NEarby Stars. The survey observational results

Context. Debris discs are a consequence of the planet formation process and constitute the fingerprints of planetesimal systems. Their solar system counterparts are the asteroid and Edgeworth-Kuiper belts. Aims. The DUNES survey aims at detecting extra-solar analogues to the Edgeworth-Kuiper belt around solar-type stars, putting in this way the solar system into context. The survey allows us to address some questions related to the prevalence and properties of planetesimal systems. Methods. We used Herschel/PACS to observe a sample of nearby FGK stars. Data at 100 and 160 μm were obtained, complemented in some cases with observations at 70 μm, and at 250, 350 and 500 μm using SPIRE. The observing strategy was to integrate as deep as possible at 100 μm to detect the stellar photosphere. Results. Debris discs have been detected at a fractional luminosity level down to several times that of the Edgeworth-Kuiper belt. The incidence rate of discs around the DUNES stars is increased from a rate of ~12.1% ± 5% before Herschel to ~20.2% ± 2%. A significant fraction (~52%) of the discs are resolved, which represents an enormous step ahead from the previously known resolved discs. Some stars are associated with faint far-IR excesses attributed to a new class of cold discs. Although it cannot be excluded that these excesses are produced by coincidental alignment of background galaxies, statistical arguments suggest that at least some of them are true debris discs. Some discs display peculiar SEDs with spectral indexes in the 70–160 μm range steeper than the Rayleigh-Jeans one. An analysis of the debris disc parameters suggests that a decrease might exist of the mean black body radius from the F-type to the K-type stars. In addition, a weak trend is suggested for a correlation of disc sizes and an anticorrelation of disc temperatures with the stellar age.

Pipe3D, a pipeline to analyze integral field spectroscopy DATA: II. Analysis sequence and CALIFA dataproducts

We present PIPE3D, an analysis pipeline based on the FIT3D fitting tool, developed to explore the properties of the stellar populations and ionized gas of integral field spectroscopy (IFS) data. PIPE3D was created to provide coherent, simple to distribute, and comparable dataproducts, independently of the origin of the data, focused on the data of the most recent IFU surveys (e.g., CALIFA, MaNGA, and SAMI), and the last generation IFS instruments (e.g., MUSE). In this article we describe the different steps involved in the analysis of the data, illustrating them by showing the dataproducts derived for NGC 2916, observed by CALIFA and P-MaNGA. As a practical example of the pipeline we present the complete set of dataproducts derived for the 200 datacubes that comprises the V500 setup of the CALIFA Data Release 2 (DR2), making them freely available through the network. Finally, we explore the hypothesis that the properties of the stellar populations and ionized gas of galaxies at the effective radius are representative of the overall average ones, finding that this is indeed the case.

Pipe3D, a pipeline to analyse integral Field Spectroscopy data: I. New fitting philosophy of FIT3D

We present an improved version of FIT3D, a fitting tool for the analysis of the spectroscopic properties of the stellar populations and the ionized gas derived from moderate resolution spectra of galaxies. This tool was developed to analyze integral field spectroscopy data and it is the basis of PIPE3D, a pipeline used in the analysis of CALIFA, MaNGA, and SAMI data. We describe the philosophy and each step of the fitting procedure. We present an extensive set of simulations in order to estimate the precision and accuracy of the derived parameters for the stellar populations and the ionized gas. We report on the results of those simulations. Finally, we compare the results of the analysis using FIT3D with those provided by other widely used packages, and we find that the parameters derived by FIT3D are fully compatible with those derived using these other tools.

Low-resolution spectroscopy and spectral energy distributions of selected sources towards σ Orionis

Aims. We study in detail nine sources in the direction of the young σ Orionis cluster, which is considered to be a unique site for studying stellar and substellar formation. The nine sources were selected because of their peculiar properties, such as extremely-red infrared colours or excessively strong Hα emission for their blue optical colours. Methods. We acquired high-quality, low-resolution spectroscopy (R ∼ 500) of the nine targets with ALFOSC at the Nordic Optical Telescope. We also re-analysed [24]-band photometry from MIPS/Spitzer and compiled the highest quality photometric dataset available at the ViJHK_s passbands and the four IRAC/Spitzer channels, for constructing accurate spectral energy distributions between 0.55 and 24 μm. Results. The nine targets were classified into: one Herbig Ae/Be star with a scattering edge-on disc; two G-type stars; one X-ray flaring, early-M, young star with chromospheric Hα emission; one very low-mass, accreting, young spectroscopic binary; two young objects at the brown-dwarf boundary with the characteristics of classical T Tauri stars; and two emission-line galaxies, one undergoing star formation, and another whose spectral energy distribution is dominated by an active galactic nucleus. We also discovered three infrared sources associated with overdensities in a cold cloud of the cluster centre. Conclusions. Low-resolution spectroscopy and spectral energy distributions are a vital tool for measuring the physical properties and evolution of young stars and candidates in the σ Orionis cluster.