Similar Scaling Relations for the Gas Content of Galaxies Across Environments to z ~ 3.5

We study the effects of the local environment on the molecular gas content of a large sample of log(M*/M⊙) ≳ 10 star-forming and starburst galaxies with specific star formation rates (sSFRs) on and above the main sequence (MS) to z ~ 3.5. ALMA observations of the dust continuum in the COSMOS field are used to estimate molecular gas masses at z ≈ 0.5–3.5. We also use a local universe sample from the ALFALFA H I survey after converting it into molecular masses. The molecular mass (M_(ISM)) scaling relation shows a dependence on z, M *, and sSFR relative to the MS, but no dependence on environmental overdensity Δ(M_(ISM) ∝ Δ^(0.03)). Similarly, gas mass fraction (f_(gas)) and depletion timescale (τ) show no environmental dependence to z ~ 3.5. At〈z〉~ 1.8, the average〈M_(ISM) 〉,〈f_(gas) 〉, and〈τ〉in densest regions is (1.6 ± 0.2) × 10^(11) M⊙, 55 ± 2%, and 0.8 ± 0.1 Gyr, respectively, similar to those in the lowest density bin. Independent of the environment, f_(gas) decreases and τincreases with increasing cosmic time. Cosmic molecular mass density (ρ) in the lowest density bins peaks at z ~ 1–2, and this peak happens at z < 1 in densest bins. This differential evolution of ρ across environments is likely due to the growth of the large-scale structure with cosmic time. Our results suggest that the molecular gas content and the subsequent star formation activity of log(M*/M⊙) ≳ 10 star-forming and starburst galaxies is primarily driven by internal processes, and not by their local environment since z ~ 3.5

The Local [C ii] 158 μm Emission Line Luminosity Function

We present, for the first time, the local [C ii] 158 μm emission line luminosity function measured using a sample of more than 500 galaxies from the Revised Bright Galaxy Sample. [C ii] luminosities are measured from the Herschel PACS observations of the Luminous Infrared Galaxies (LIRGs) in the Great Observatories All-sky LIRG Survey and estimated for the rest of the sample based on the far-infrared (far-IR) luminosity and color. The sample covers 91.3% of the sky and is complete at S_(60 μm) > 5.24 Jy. We calculate the completeness as a function of [C ii] line luminosity and distance, based on the far-IR color and flux densities. The [C ii] luminosity function is constrained in the range ~10^(7–9) L_⊙ from both the 1/V_(max) and a maximum likelihood methods. The shape of our derived [C ii] emission line luminosity function agrees well with the IR luminosity function. For the CO(1-0) and [C ii] luminosity functions to agree, we propose a varying ratio of [C ii]/CO(1-0) as a function of CO luminosity, with larger ratios for fainter CO luminosities. Limited [C ii] high-redshift observations as well as estimates based on the IR and UV luminosity functions are suggestive of an evolution in the [C ii] luminosity function similar to the evolution trend of the cosmic star formation rate density. Deep surveys using the Atacama Large Millimeter Array with full capability will be able to confirm this prediction

Resolving the AGN and host emission in the mid-infrared using a model-independent spectral decomposition

We present results on the spectral decomposition of 118 Spitzer Infrared Spectrograph (IRS) spectra from local active galactic nuclei (AGN) using a large set of Spitzer/IRS spectra as templates. The templates are themselves IRS spectra from extreme cases where a single physical component (stellar, interstellar, or AGN) completely dominates the integrated mid-infrared emission. We show that a linear combination of one template for each physical component reproduces the observed IRS spectra of AGN hosts with unprecedented fidelity for a template fitting method, with no need to model extinction separately. We use full probability distribution functions to estimate expectation values and uncertainties for observables, and find that the decomposition results are robust against degeneracies. Furthermore, we compare the AGN spectra derived from the spectral decomposition with sub-arcsecond resolution nuclear photometry and spectroscopy from ground-based observations. We find that the AGN component derived from the decomposition closely matches the nuclear spectrum, with a 1-sigma dispersion of 0.12 dex in luminosity and typical uncertainties of ~0.19 in the spectral index and ~0.1 in the silicate strength. We conclude that the emission from the host galaxy can be reliably removed from the IRS spectra of AGN. This allows for unbiased studies of the AGN emission in intermediate and high redshift galaxies -currently inaccesible to ground-based observations- with archival Spitzer/IRS data and in the future with the Mid-InfraRed Instrument of the James Webb Space Telescope. The decomposition code and templates are available at http://www.denebola.org/ahc/deblendIRS.Comment: 16 pages, 15 figures, 2 tables, accepted for publication in Ap

Understanding the 8 micron vs. Pa-alpha relationship on sub-arcsecond scales in Luminous Infrared Galaxies

This work explores in detail the relation between the 8 micron and the Pa-alpha emissions for 122 HII regions identified in a sample of 10 low-z LIRGs with nearly constant metallicity (12 + log (O/H) ~ 8.8). We use Gemini/T-ReCS high-spatial resolution (<~ 0.4" ~ 120 pc for the average distance of 60 Mpc of our sample) mid-infrared imaging (at 8.7 micron or 10.3 micron) together with HST/NICMOS continuum and Pa-alpha images. The LIRG HII regions extend the L_8micron vs. L_Pa-alpha relation found for HII knots in the high-metallicity SINGS galaxies by about two orders of magnitude to higher luminosities. Since the metallicity of the LIRG sample is nearly constant, we can rule out this effect as a cause for the scatter seen in the relationship. In turn, it is attributed to two effects: age and PAH features. The L_8micron/L_Pa-alpha ratio, which varies by a factor of ten for the LIRG HII regions, is reproduced by a model with instantaneous star formation and ages ranging from ~ 4 to 7.5 Myr. The remaining dispersion around the model predictions for a given age is probably due to differential contributions of the PAH features (the 8.6 micron, in our case) to the 8 micron emission from galaxy to galaxy.Comment: 16 pages, 9 figures, accepted for publication in ApJ; paper with full-resolution figures can be found at: http://damir.iem.csic.es/extragalactic

Measuring the Heating and Cooling of the Interstellar Medium at High Redshift: PAH and [C II] Observations of the Same Star-forming Galaxies at z ∼ 2

Star formation depends critically on cooling mechanisms in the interstellar medium (ISM); however, thermal properties of gas in galaxies at the peak epoch of star formation (z ~ 2) remain poorly understood. A limiting factor in understanding the multiphase ISM is the lack of multiple tracers detected in the same galaxies, such as Polycyclic Aromatic Hydrocarbon (PAH) emission, a tracer of a critical photoelectric heating mechanism in interstellar gas, and [C ii] 158 μm fine-structure emission, a principal coolant. We present ALMA Band 9 observations targeting [C ii] in six z ~ 2 star-forming galaxies with strong Spitzer IRS detections of PAH emission. All six galaxies are detected in dust continuum and marginally resolved. We compare the properties of PAH and [C ii] emission, and constrain their relationship as a function of total infrared luminosity (L_(IR)) and IR surface density. [C ii] emission is detected in one galaxy at high signal-to-noise (34σ), and we place a secure upper limit on a second source. The rest of our sample are not detected in [C ii] likely due to redshift uncertainties and narrow ALMA bandpass windows. Our results are consistent with the deficit in [C ii]/L_(IR) and PAH/L_(IR) observed in the literature. However, the ratio of [C ii] to PAH emission at z ~ 2 is possibly much lower than what is observed in nearby dusty star-forming galaxies. This could be the result of enhanced cooling via [O i] at high-z, hotter gas and dust temperatures, and/or a reduction in the photoelectric efficiency, in which the coupling between interstellar radiation and gas heating is diminished

The Role of Star Formation and an AGN in Dust Heating of z = 0.3–2.8 Galaxies. I. Evolution with Redshift and Luminosity

We characterize infrared spectral energy distributions of 343 (ultra)luminous infrared galaxies from z = 0.3–2.8. We diagnose the presence of an active galactic nucleus (AGN) by decomposing individual Spitzer mid-IR spectroscopy into emission from star formation and an AGN-powered continuum; we classify sources as star-forming galaxies (SFGs), AGNs, or composites. Composites comprise 30% of our sample and are prevalent at faint and bright S_(24), making them an important source of IR AGN emission. We combine spectroscopy with multiwavelength photometry, including Herschel imaging, to create three libraries of publicly available templates (2–1000 μm). We fit the far-IR emission using a two-temperature modified blackbody to measure cold and warm dust temperatures (T_c and T_w). We find that T_c does not depend on mid-IR classification, while T_w shows a notable increase as the AGN grows more luminous. We measure a quadratic relationship between mid-IR AGN emission and total AGN contribution to LIR. AGNs, composites, and SFGs separate in S_8/S_(3.6) and S_(250)/S_(24), providing a useful diagnostic for estimating relative amounts of these sources. We estimate that >40% of IR-selected samples host an AGN, even at faint selection thresholds (S_(24) > 100 μJy). Our decomposition technique and color diagnostics are relevant given upcoming observations with the James Webb Space Telescope

Measuring the Heating and Cooling of the Interstellar Medium at High Redshift: PAH and [C II] Observations of the Same Star-forming Galaxies at z ∼ 2

Star formation depends critically on cooling mechanisms in the interstellar medium (ISM); however, thermal properties of gas in galaxies at the peak epoch of star formation (z ~ 2) remain poorly understood. A limiting factor in understanding the multiphase ISM is the lack of multiple tracers detected in the same galaxies, such as Polycyclic Aromatic Hydrocarbon (PAH) emission, a tracer of a critical photoelectric heating mechanism in interstellar gas, and [C ii] 158 μm fine-structure emission, a principal coolant. We present ALMA Band 9 observations targeting [C ii] in six z ~ 2 star-forming galaxies with strong Spitzer IRS detections of PAH emission. All six galaxies are detected in dust continuum and marginally resolved. We compare the properties of PAH and [C ii] emission, and constrain their relationship as a function of total infrared luminosity (L_(IR)) and IR surface density. [C ii] emission is detected in one galaxy at high signal-to-noise (34σ), and we place a secure upper limit on a second source. The rest of our sample are not detected in [C ii] likely due to redshift uncertainties and narrow ALMA bandpass windows. Our results are consistent with the deficit in [C ii]/L_(IR) and PAH/L_(IR) observed in the literature. However, the ratio of [C ii] to PAH emission at z ~ 2 is possibly much lower than what is observed in nearby dusty star-forming galaxies. This could be the result of enhanced cooling via [O i] at high-z, hotter gas and dust temperatures, and/or a reduction in the photoelectric efficiency, in which the coupling between interstellar radiation and gas heating is diminished

Star-formation histories of local luminous infrared galaxies

We present the analysis of the integrated spectral energy distribution (SED) from the ultraviolet (UV) to the far-infrared and H$\alpha$ of a sample of 29 local systems and individual galaxies with infrared (IR) luminosities between 10^11 Lsun and 10^11.8 Lsun. We have combined new narrow-band H$\alpha$+[NII] and broad-band g, r optical imaging taken with the Nordic Optical Telescope (NOT), with archival GALEX, 2MASS, Spitzer, and Herschel data. The SEDs (photometry and integrated H$\alpha$ flux) have been fitted with a modified version of the MAGPHYS code using stellar population synthesis models for the UV-near-IR range and thermal emission models for the IR emission taking into account the energy balance between the absorbed and re-emitted radiation. From the SED fits we derive the star-formation histories (SFH) of these galaxies. For nearly half of them the star-formation rate appears to be approximately constant during the last few Gyrs. In the other half, the current star-formation rate seems to be enhanced by a factor of 3-20 with respect to that occured ~1 Gyr ago. Objects with constant SFH tend to be more massive than starbursts and they are compatible with the expected properties of a main-sequence (M-S) galaxy. Likewise, the derived SFHs show that all our objects were M-S galaxies ~1 Gyr ago with stellar masses between 10^10.1 and 10^11.5 Msun. We also derived from our fits the average extinction (A_v=0.6-3 mag) and the polycyclic aromatic hydrocarbons (PAH) luminosity to L(IR) ratio (0.03-0.16). We combined the A_v with the total IR and H$\alpha$ luminosities into a diagram which can be used to identify objects with rapidly changing (increasing or decreasing) SFR during the last 100 Myr.Comment: 16 pages + online material, accepted for publication in A&

Low frequency radio continuum imaging and SED modeling of 11 LIRGs: radio-only and FUV to radio bands

We present the detailed analysis of 11 local luminous infrared galaxies (LIRGs) from ultraviolet through far-infrared to radio ($\sim$70 MHz to $\sim$15 GHz) bands. We derive the astrophysical properties through spectral energy distribution (SED) modeling using the Code Investigating GALaxy Emission (CIGALE) and UltraNest codes. The radio SEDs include our new observations at 325 and 610 MHz from the GMRT and the measurements from public archives. Our main results are (1) radio SEDs show turnovers and bends, (2) the synchrotron spectral index of the fitted radio spectra ranges between $-$0.5 and $-$1.7, and (3) the infrared luminosity, dust mass, dust temperature, stellar mass, star-formation rates (SFRs) and AGN fraction obtained from CIGALE falls in the range exhibited by galaxies of the same class. The ratio of 60$\mu$m infrared and 1.4 GHz radio luminosity, the 1.4 GHz thermal fraction, and emission measure range between 2.1 and 2.9, 0.1% and 10%, 0.02 and 269.5$\times$10$^{6}$ cm$^{-6}$ pc, respectively. We conclude that the turnovers seen in the radio SEDs are due to free-free absorption; this is supported by the low AGN fraction derived from the CIGALE analysis. The decomposed 1.4 GHz thermal and nonthermal radio luminosities allowed us to compute the star formation rate (SFR) using scaling relations. A positive correlation is observed between the SFR$_{IR}$ obtained 10 Myr ago (compared to 100 Myr ago) and 1.4 GHz radio (total and nonthermal) because similar synchrotron lifetimes are expected for typical magnetic field strengths observed in these galaxies ($\approx$50$\mu$G).Comment: ApJ accepted. Comments are welcom