Radio continuum properties of luminous infrared galaxies. Identifying the presence of an AGN in the radio

Luminous infrared galaxies are systems enshrouded in dust, which absorbs most of their optical/UV emission and re-radiates it in the mid- and far-infrared. Radio observations are largely unaffected by dust obscuration, enabling us to study the central regions of LIRGs in an unbiased manner. The main goal of this project is to examine how the radio properties of local LIRGs relate to their infrared spectral characteristics. Here we present an analysis of the radio continuum properties of a subset of the Great Observatories All-sky LIRG Survey (GOALS), which consists of 202 nearby systems (z<0.088). Our radio sample consists of 35 systems, or 46 individual galaxies, that were observed at both 1.49 and 8.44 GHz with the VLA with a resolution of about 1 arcsec (FWHM). The aim of the project is to use the radio imagery to probe the central kpc of these LIRGs in search of active galactic nuclei. We used the archival data at 1.49 and 8.44 GHz to create radio-spectral-index maps using the standard relation between flux density Sv and frequency v, S~v^-a, where a is the radio spectral index. By studying the spatial variations in a, we classified the objects as radio-AGN, radio-SB, and AGN/SB (a mixture). We identified the presence of an active nucleus using the radio morphology, deviations from the radio/infrared correlation, and spatially resolved spectral index maps, and then correlated this to the usual mid-infrared ([NeV]/[NeII] and [OIV]/[NeII] line ratios and EQW of the 6.2 um PAH feature) and optical (BPT diagram) AGN diagnostics. We find that 21 out of the 46 objects in our sample are radio-AGN, 9 are classified as starbursts (SB), and 16 are AGN/SB. After comparing to other AGN diagnostics we find 3 objects out of the 46 that are identified as AGN based on the radio analysis, but are not classified as such based on the mid-infrared and optical AGN diagnostics presented in this study.Comment: 33 pages, 7 figures, 5 tables, to appear in A&

A 33 GHz Survey of Local Major Mergers: Estimating the Sizes of the Energetically Dominant Regions from High Resolution Measurements of the Radio Continuum

We present Very Large Array observations of the 33 GHz radio continuum emission from 22 local ultraluminous and luminous infrared (IR) galaxies (U/LIRGs). These observations have spatial (angular) resolutions of 30--720 pc (0.07"-0.67") in a part of the spectrum that is likely to be optically thin. This allows us to estimate the size of the energetically dominant regions. We find half-light radii from 30 pc to 1.7 kpc. The 33 GHz flux density correlates well with the IR emission, and we take these sizes as indicative of the size of the region that produces most of the energy. Combining our 33 GHz sizes with unresolved measurements, we estimate the IR luminosity and star formation rate per area, and the molecular gas surface and volume densities. These quantities span a wide range (4 dex) and include some of the highest values measured for any galaxy (e.g., $\mathrm{\Sigma_{SFR}^{33GHz} \leq 10^{4.1} M_{\odot} yr^{-1} kpc^{-2}}$). At least $13$ sources appear Compton thick ($\mathrm{N_{H}^{33GHz} \geq 10^{24} cm^{-2}}$). Consistent with previous work, contrasting these data with observations of normal disk galaxies suggests a nonlinear and likely multi-valued relation between SFR and molecular gas surface density, though this result depends on the adopted CO-to-H$_{2}$ conversion factor and the assumption that our 33 GHz sizes apply to the gas. 11 sources appear to exceed the luminosity surface density predicted for starbursts supported by radiation pressure and supernovae feedback, however we note the need for more detailed observations of the inner disk structure. U/LIRGs with higher surface brightness exhibit stronger [{\sc Cii}] 158$\mu$m deficits, consistent with the suggestion that high energy densities drive this phenomenon.Comment: 32 pages, 11 figures, 7 tables. Accepted for publication in Ap

Star-forming Clumps in Local Luminous Infrared Galaxies

We present HST narrowband near-infrared imaging of Paα and Paβ emission of 48 local luminous infrared galaxies (LIRGs) from the Great Observatories All-Sky LIRG Survey. These data allow us to measure the properties of 810 spatially resolved star-forming regions (59 nuclei and 751 extranuclear clumps) and directly compare their properties to those found in both local and high-redshift star-forming galaxies. We find that in LIRGs the star-forming clumps have radii ranging from ~90 to 900 pc and star formation rates (SFRs) of ~1 × 10⁻³ to 10 M⊙ yr⁻¹, with median values for extranuclear clumps of 170 pc and 0.03 M⊙ yr⁻¹. The detected star-forming clumps are young, with a median stellar age of 8.7 Myr, and have a median stellar mass of 5 × 10⁵ M ⊙. The SFRs span the range of those found in normal local star-forming galaxies to those found in high-redshift star-forming galaxies at z = 1–3. The luminosity function of the LIRG clumps has a flatter slope than found in lower-luminosity, star-forming galaxies, indicating a relative excess of luminous star-forming clumps. In order to predict the possible range of star-forming histories and gas fractions, we compare the star-forming clumps to those measured in the MassiveFIRE high-resolution cosmological simulation. The star-forming clumps in MassiveFIRE cover the same range of SFRs and sizes found in the local LIRGs and have total gas fractions that extend from 10% to 90%. If local LIRGs are similar to these simulated galaxies, we expect that future observations with ALMA will find a large range of gas fractions, and corresponding star formation efficiencies, among the star-forming clumps in LIRGs

Extended [CII] Emission in Local Luminous Infrared Galaxies

We present Herschel/PACS observations of extended [CII]157.7{\mu}m line emission detected on ~ 1 - 10 kpc scales in 60 local luminous infrared galaxies (LIRGs) from the Great Observatories All-sky LIRG Survey (GOALS). We find that most of the extra-nuclear emission show [CII]/FIR ratios >~ 4 x 10^-3, larger than the mean ratio seen in the nuclei, and similar to those found in the extended disks of normal star-forming galaxies and the diffuse inter-stellar medium (ISM) of our Galaxy. The [CII] "deficits" found in the most luminous local LIRGs are therefore restricted to their nuclei. There is a trend for LIRGs with warmer nuclei to show larger differences between their nuclear and extra-nuclear [CII]/FIR ratios. We find an anti-correlation between [CII]/FIR and the luminosity surface density, {\Sigma}_IR, for the extended emission in the spatially-resolved galaxies. However, there is an offset between this trend and that found for the LIRG nuclei. We use this offset to derive a beam filling-factor for the star-forming regions within the LIRG disks of ~ 6 % relative to their nuclei. We confront the observed trend to photo-dissociation region (PDR) models and find that the slope of the correlation is much shallower than the model predictions. Finally, we compare the correlation found between [CII]/FIR and {\Sigma}_IR with measurements of high-redshift starbursting IR-luminous galaxies.Comment: Accepted for publication in ApJ Letters. 5 pages, 2 figure

A Herschel/PACS Far Infrared Line Emission Survey of Local Luminous Infrared Galaxies

We present an analysis of [OI]63, [OIII]88, [NII]122 and [CII]158 far-infrared (FIR) fine-structure line observations obtained with Herschel/PACS, for ~240 local luminous infrared galaxies (LIRGs) in the Great Observatories All-sky LIRG Survey (GOALS). We find pronounced declines -deficits- of line-to-FIR-continuum emission for [NII]122, [OI]63 and [CII]158 as a function of FIR color and infrared luminosity surface density, $\Sigma_{\rm IR}$. The median electron density of the ionized gas in LIRGs, based on the [NII]122/[NII]205 ratio, is $n_{\rm e}$ = 41 cm$^{-3}$. We find that the dispersion in the [CII]158 deficit of LIRGs is attributed to a varying fractional contribution of photo-dissociation-regions (PDRs) to the observed [CII]158 emission, f([CII]PDR) = [CII]PDR/[CII], which increases from ~60% to ~95% in the warmest LIRGs. The [OI]63/[CII]158PDR ratio is tightly correlated with the PDR gas kinetic temperature in sources where [OI]63 is not optically-thick or self-absorbed. For each galaxy, we derive the average PDR hydrogen density, $n_{\rm H}$, and intensity of the interstellar radiation field, in units of G$_0$, and find G$_0$/$n_{\rm H}$ ratios ~0.1-50 cm$^3$, with ULIRGs populating the upper end of the distribution. There is a relation between G$_0$/$n_{\rm H}$ and $\Sigma_{\rm IR}$, showing a critical break at $\Sigma_{\rm IR}^{\star}$ ~ 5 x 10$^{10}$ Lsun/kpc$^2$. Below $\Sigma_{\rm IR}^{\star}$, G$_0$/$n_{\rm H}$ remains constant, ~0.32 cm$^3$, and variations in $\Sigma_{\rm IR}$ are driven by the number density of star-forming regions within a galaxy, with no change in their PDR properties. Above $\Sigma_{\rm IR}^{\star}$, G$_0$/$n_{\rm H}$ increases rapidly with $\Sigma_{\rm IR}$, signaling a departure from the typical PDR conditions found in normal star-forming galaxies towards more intense/harder radiation fields and compact geometries typical of starbursting sources.Comment: 23 pages, 12 figures. Submitted to ApJ, including 2nd round of referee comments. Data tables can be found at: http://goals.ipac.caltech.edu/ and will also appear on the electronic version of the journa

The Buried Starburst in the Interacting Galaxy II Zw 096 as Revealed by the Spitzer Space Telescope

An analysis of data from the Spitzer Space Telescope, Hubble Space Telescope, Chandra X-ray Observatory, and AKARI Infrared Astronomy Satellite is presented for the z = 0.036 merging galaxy system II Zw 096 (CGCG 448-020). Because II Zw 096 has an infrared luminosity of log(L_(IR)/L_☉) = 11.94, it is classified as a Luminous Infrared Galaxy (LIRG), and was observed as part of the Great Observatories All-sky LIRG Survey (GOALS). The Spitzer data suggest that 80% of the total infrared luminosity comes from an extremely compact, red source not associated with the nuclei of the merging galaxies. The Spitzer mid-infrared spectra indicate no high-ionization lines from a buried active galactic nucleus in this source. The strong detection of the 3.3 μm and 6.2 μm polycyclic aromatic hydrocarbon emission features in the AKARI and Spitzer spectra also implies that the energy source of II Zw 096 is a starburst. Based on Spitzer infrared imaging and AKARI near-infrared spectroscopy, the star formation rate is estimated to be 120 M_☉ yr^(-1) and >45 M_☉ yr^(-1), respectively. Finally, the high-resolution B-, I-, and H-band images show many star clusters in the interacting system. The colors of these clusters suggest at least two populations—one with an age of 1-5 Myr and one with an age of 20-500 Myr, reddened by 0-2 mag of visual extinction. The masses of these clusters span a range between 10^6 and 10^8 M_☉. This starburst source is reminiscent of the extranuclear starburst seen in NGC 4038/9 (the Antennae Galaxies) and Arp 299 but approximately an order of magnitude more luminous than the Antennae. The source is remarkable in that the off-nuclear infrared luminosity dominates the entire system

Mid-Infrared Properties of Nearby Luminous Infrared Galaxies I: Spitzer IRS Spectra for the GOALS Sample

The Great Observatories All-Sky LIRG Survey (GOALS) is a multiwavelength study of luminous infrared galaxies (LIRGs) in the local universe. Here we present low resolution Spitzer spectra covering 5-38um and provide a basic analysis of the mid-IR spectral properties for nearby LIRGs. In a companion paper, we discuss detailed fits to the spectra. The GOALS sample of 244 nuclei in 180 luminous and 22 ultraluminous IR galaxies represents a complete subset of the IRAS RBGS and covers a range of merger stages, morphologies and spectral types. The majority (>60%) of GOALS LIRGs have high 6.2um PAH equivalent widths (EQW > 0.4um) and low levels of silicate absorption (s_9.7um >-1.0). There is a general trend among the U/LIRGs for silicate depth and MIR slope to increase with LIR. U/LIRGs in the late stages of a merger also have on average steeper MIR slopes and higher levels of dust obscuration. Together these trends suggest that as gas & dust is funneled towards the center of a coalescing merger, the nuclei become more compact and obscured. The sources that depart from these correlations have very low PAH EQW (EQW < 0.1um) consistent with their MIR emission being dominated by an AGN. The most heavily dust obscured sources are the most compact in their MIR emission, suggesting that the obscuring (cool) dust is associated with the outer regions of the starburst. As the merger progresses a marked decline is seen for the fraction of high EQW (star formation dominated) sources while the fraction of composite sources increases but the fraction of AGN-dominated sources remains low. When compared to the MIR spectra of submillimeter galaxies (SMGs) at z~2, the average GOALS LIRG is more absorbed at 9.7um and has more PAH emission. However, when the AGN contributions to both the local LIRGs and the high-z SMGs are removed, the average local starbursting LIRG closely resembles the starbursting SMGs.Comment: Accepted for publication in ApJ

Star Formation Histories of the LEGUS Dwarf Galaxies (I): recent History of NGC1705, NGC4449 and Holmberg II

We use HST observations from the Legacy Extragalactic UV Survey to reconstruct the recent star formation histories (SFHs) of three actively star-forming dwarf galaxies, NGC4449, Holmberg II and NGC1705, from their UV color-magnitude diagrams (CMDs). We apply a CMD fitting technique using two independent sets of stellar isochrones, PARSEC-COLIBRI and MIST, to assess the uncertainties related to stellar evolution modelling. Irrespective of the adopted stellar models, all the three dwarfs are found to have had almost constant star formation rates (SFRs) in the last 100-200 Myr, with modest enhancements (a factor of $\sim$2) above the 100 Myr-averaged-SFR. Significant differences among the three dwarfs are found in the overall SFR, the timing of the most recent peak and the SFR$/$area. The Initial Mass Function (IMF) of NGC1705 and Holmberg II is consistent with a Salpeter slope down to $\approx$ 5 M$_{\odot}$, whereas it is slightly flatter, s$=-2.0$, in NGC4449. The SFHs derived with the two different sets of stellar models are consistent with each other, except for some quantitative details, attributable to their input assumptions. They also share the drawback that all synthetic diagrams predict a clear separation in color between upper main sequence and helium burning stars, which is not apparent in the data. Since differential reddening, significant in NGC4449, or unresolved binaries don't appear to be sufficient to fill the gap, we suggest this calls for a revision of both sets of stellar evolutionary tracks.Comment: 22 pages, 17 figures, accepted for publication on Ap

The Spatial Extent of (U)LIRGs in the mid-Infrared I: The Continuum Emission

We present an analysis of the extended mid-infrared (MIR) emission of the Great Observatories All-Sky LIRG Survey (GOALS) sample based on 5-15um low resolution spectra obtained with the IRS on Spitzer. We calculate the fraction of extended emission as a function of wavelength for the galaxies in the sample, FEE_lambda. We can identify 3 general types of FEE_lambda: one where it is constant, one where features due to emission lines and PAHs appear more extended than the continuum, and a third which is characteristic of sources with deep silicate absorption at 9.7um. More than 30% of the galaxies have a median FEE_lambda larger than 0.5 implying that at least half of their MIR emission is extended. Luminous Infrared Galaxies (LIRGs) display a wide range of FEE in their warm dust continuum (0<=FEE_13.2um<=0.85). The large values of FEE_13.2um that we find in many LIRGs suggest that their extended MIR continuum emission originates in scales up to 10kpc. The mean size of the LIRG cores at 13.2um is 2.6kpc. However, once the LIR of the systems reaches the threshold of ~10^11.8Lsun, all sources become clearly more compact, with FEE_13.2um<=0.2, and their cores are unresolved. Our estimated upper limit for the core size of ULIRGs is less than 1.5kpc. The analysis indicates that the compactness of systems with LIR>~10^11.25Lsun strongly increases in those classified as mergers in their final stage of interaction. The FEE_13.2um is also related to the contribution of an active galactic nucleus (AGN) to the MIR. Galaxies which are more AGN-dominated are less extended, independently of their LIR. We finally find that the extent of the MIR continuum emission is correlated with the far-IR IRAS log(f_60um/f_100um) color. This enables us to place a lower limit to the area in a galaxy from where the cold dust emission may originate, a prediction which can be tested soon with the Herschel Space Telescope.Comment: 18 pages, 8 figures, accepted for publication in Ap