Great Observatories All-sky LIRG Survey Research Articles

THE LOCAL [C ii] 158 μm EMISSION LINE LUMINOSITY FUNCTION

ABSTRACT 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 ∼107–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.

Morphological classification of local luminous infrared galaxies

We present an analysis of the morphological classification of 89 luminous infrared galaxies (LIRGs) from the Great Observatories All-sky LIRG Survey (GOALS) sample using non-parametric coefficients and compare their morphology as a function of wavelength. We rely on images obtained in the optical (B- and I-band) as well as in the infrared (H-band and 5.8$\mu$m). Our classification is based on the calculation of $Gini$ and the second order of light ($M_{20}$) non-parametric coefficients which we explore as a function of stellar mass ($M_\star$), infrared luminosity ($L_{IR}$) and star formation rate (SFR). We investigate the relation between $M_{20}$, the specific SFR (sSFR) and the dust temperature ($T_{dust}$) in our galaxy sample. We find that $M_{20}$ is a better morphological tracer than $Gini$, as it allows to distinguish systems formed by double systems from isolated and post-merger LIRGs. The multi-wavelength analysis allows us to identify a region in the $Gini$-$M_{20}$ parameter space where ongoing mergers reside, regardless of the band used to calculate the coefficients. In particular when measured in the H-band, this region can be used to identify ongoing mergers, with a minimal contamination from LIRGs in other stages. We also find that while the sSFR is positively correlated with $M_{20}$ when measured in the mid-infrared, i.e. star-bursting galaxies show more compact emission, it is anti-correlated with the B-band based $M_{20}$. We interpret this as the spatial decoupling between obscured and un-obscured star formation, whereby the ultraviolet/optical size of a LIRGs experience an intense dust enshrouded central starburst is larger than in the one in the mid-infrared since the contrast between the nuclear to the extended disk emission is smaller in the mid-infrared. This has important implications for high redshift surveys of dusty sources. [abridged]

GALAXY MERGERS DRIVE SHOCKS: AN INTEGRAL FIELD STUDY OF GOALS GALAXIES

We present an integral field spectroscopic study of radiative shocks in 27 nearby ultraluminous and luminous infrared galaxies (U/LIRGs) from the Great Observatory All-sky LIRG Survey, a subset of the Revised Bright Galaxy Sample. Our analysis of the resolved spectroscopic data from the Wide Field Spectrograph focuses on determining the detailed properties of the emission-line gas, including a careful treatment of multicomponent emission-line profiles. The resulting information obtained from the spectral fits is used to map the kinematics of the gas, sources of ionizing radiation, and feedback present in each system. The resulting properties are tracked as a function of merger stage. Using emission-line flux ratios and velocity dispersions, we find evidence for widespread, extended shock excitation in many local U/LIRGs. These low-velocity shocks become an increasingly important component of the optical emission lines as a merger progresses. We find that shocks may account for as much as half of the Hα luminosity in the latest-stage mergers in our sample. We discuss some possible implications of our result and consider the presence of active galactic nuclei and their effects on the spectra in our sample.

EXCITATION MECHANISMS FOR HCN(1–0) AND HCO+(1–0) IN GALAXIES FROM THE GREAT OBSERVATORIES ALL-SKY LIRG SURVEY*

ABSTRACT We present new Institut de Radioastronomie Millimétrique (IRAM) 30 m spectroscopic observations of the ∼88 GHz band, including emission from the CCH ( N = 1 → 0 ) ?> multiplet, HCN ( J = 1 → 0 ) ?> , HCO + ( J = 1 → 0 ) ?> , and HNC ( J = 1 → 0 ) ?> , for a sample of 58 local luminous and ultraluminous infrared galaxies from the Great Observatories All-sky LIRG Survey (GOALS). By combining our new IRAM data with literature data and Spitzer/IRS spectroscopy, we study the correspondence between these putative tracers of dense gas and the relative contribution of active galactic nuclei (AGNs) and star formation to the mid-infrared luminosity of each system. We find the HCN (1–0) emission to be enhanced in AGN-dominated systems ( ⟨ ?> L HCN ( 1 – 0 ) ′ ?> / L HCO + ( 1 – 0 ) ′ ?> ⟩ = 1.84 ?> ), compared to composite and starburst-dominated systems ( ⟨ ?> L HCN ( 1 – 0 ) ′ ?> / L HCO + ( 1 – 0 ) ′ ?> ⟩ = 1.14 , ?> and 0.88, respectively). However, some composite and starburst systems have L HCN ( 1 – 0 ) ′ ?> / L HCO + ( 1 – 0 ) ′ ?> ratios comparable to those of AGNs, indicating that enhanced HCN emission is not uniquely associated with energetically dominant AGNs. After removing AGN-dominated systems from the sample, we find a linear relationship (within the uncertainties) between log10( L HCN ( 1 – 0 ) ′ ?> ) and log10(L IR), consistent with most previous findings. L HCN ( 1 – 0 ) ′ ?> /L IR, typically interpreted as the dense-gas depletion time, appears to have no systematic trend with L IR for our sample of luminous and ultraluminous infrared galaxies, and has significant scatter. The galaxy-integrated L HCN ( 1 – 0 ) ′ ?> and L HCO + ( 1 – 0 ) ′ ?> emission do not appear to have a simple interpretation in terms of the AGN dominance or the star formation rate, and are likely determined by multiple processes, including density and radiative effects.

Radio continuum properties of luminous infrared galaxies

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.

THE FUV TO NEAR-IR MORPHOLOGIES OF LUMINOUS INFRARED GALAXIES IN THE GOALS SAMPLE

We compare the morphologies of a sample of 20 LIRGs from the Great Observatories All-sky LIRG Survey (GOALS) in the FUV, B, I and H bands, using the Gini (G) and M20 parameters to quantitatively estimate the distribution and concentration of flux as a function of wavelength. HST images provide an average spatial resolution of ~80 pc. While our LIRGs can be reliably classified as mergers across the entire range of wavelengths studied here, there is a clear shift toward more negative M20 (more bulge-dominated) and a less significant decrease in G values at longer wavelengths. We find no correlation between the derived FUV G-M20 parameters and the global measures of the IR to FUV flux ratio, IRX. Given the fine resolution in our HST data, this suggests either that the UV morphology and IRX are correlated on very small scales, or that the regions emitting the bulk of the IR emission emit almost no FUV light. We use our multi-wavelength data to simulate how merging LIRGs would appear from z~0.5-3 in deep optical and near-infrared images such as the HUDF, and use these simulations to measure the G-M20 at these redshifts. Our simulations indicate a noticeable decrease in G, which flattens at z >= 2 by as much as 40%, resulting in mis-classifying our LIRGs as disk-like, even in the rest-frame FUV. The higher redshift values of M20 for the GOALS sources do not appear to change more than about 10% from the values at z~0. The change in G-M20 is caused by the surface brightness dimming of extended tidal features and asymmetries, and also the decreased spatial resolution which reduced the number of individual clumps identified. This effect, seen as early as z~0.5, could easily lead to an underestimate of the number of merging galaxies at high-redshift in the rest-frame FUV.

MID-INFRARED PROPERTIES OF LUMINOUS INFRARED GALAXIES. II. PROBING THE DUST AND GAS PHYSICS OF THE GOALS SAMPLE

The Great Observatories All-sky LIRG Survey (GOALS) is a comprehensive, multiwavelength study of luminous infrared galaxies (LIRGs) in the local universe. Here, we present the results of a multi-component, spectral decomposition analysis of the low-resolution mid-infrared (MIR) Spitzer Infrared Spectrograph spectra from 5–38 μm of 244 LIRG nuclei. The detailed fits and high-quality spectra allow for characterization of the individual polycyclic aromatic hydrocarbon (PAH) features, warm molecular hydrogen emission, and optical depths for both silicate dust grains and water ices. We find that starbursting LIRGs, which make up the majority of the GOALS sample, are very consistent in their MIR properties (i.e., τ9.7 μm, τice, neon line ratios, and PAH feature ratios). However, as their EQW6.2 μm decreases, usually an indicator of an increasingly dominant active galactic nucleus (AGN), LIRGs cover a larger spread in these MIR parameters. The contribution from PAH emission to the total IR luminosity (L(PAH)/L(IR)) in LIRGs varies from 2%–29% and LIRGs prior to their first encounter show significantly higher L(PAH)/L(IR) ratios on average. We observe a correlation between the strength of the starburst (represented by IR8 = LIR/L8 μm) and the PAH fraction at 8 μm but no obvious link between IR8 and the 7.7 to 11.3 PAH ratio, suggesting that the fractional photodissociation region (PDR) emission, and not the overall grain properties, is associated with the rise in IR8 for galaxies off the starburst main sequence. We detect crystalline silicate features in ∼6% of the sample but only in the most obscure sources (s9.7 μm < −1.24). Ice absorption features are observed in ∼11% (56%) of GOALS LIRGs (ULIRGs) in sources with a range of silicate depths. Most GOALS LIRGs have L(H2)/L(PAH) ratios elevated above those observed for normal star-forming galaxies and exhibit a trend for increasing L(H2)/L(PAH) ratio with increasing L(H2). While star formation appears to be the dominant process responsible for exciting the H2 in most of the GOALS galaxies, a subset of LIRGs (∼10%) shows excess H2 emission that is inconsistent with PDR models and may be excited by shocks or AGN-induced outflows.

EXTENDED [C II] 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.

WARM MOLECULAR GAS IN LUMINOUS INFRARED GALAXIES

We present our initial results on the CO rotational spectral line energy distribution (SLED) of the $J$ to $J$$-$1 transitions from $J=4$ up to $13$ from Herschel SPIRE spectroscopic observations of 65 luminous infrared galaxies (LIRGs) in the Great Observatories All-Sky LIRG Survey (GOALS). The observed SLEDs change on average from one peaking at $J \le 4$ to a broad distribution peaking around $J \sim\,$6$-$7 as the IRAS 60-to-100 um color, $C(60/100)$, increases. However, the ratios of a CO line luminosity to the total infrared luminosity, $L_{\rm IR}$, show the smallest variation for $J$ around 6 or 7. This suggests that, for most LIRGs, ongoing star formation (SF) is also responsible for a warm gas component that emits CO lines primarily in the mid-$J$ regime ($5 \lesssim J \lesssim 10$). As a result, the logarithmic ratios of the CO line luminosity summed over CO (5$-$4), (6$-$5), (7$-$6), (8$-$7) and (10$-$9) transitions to $L_{\rm IR}$, $\log R_{\rm midCO}$, remain largely independent of $C(60/100)$, and show a mean value of $-4.13$ ($\equiv \log R^{\rm SF}_{\rm midCO}$) and a sample standard deviation of only 0.10 for the SF-dominated galaxies. Including additional galaxies from the literature, we show, albeit with small number of cases, the possibility that galaxies, which bear powerful interstellar shocks unrelated to the current SF, and galaxies, in which an energetic active galactic nucleus contributes significantly to the bolometric luminosity, have their $R_{\rm midCO}$ higher and lower than $R^{\rm SF}_{\rm midCO}$, respectively.

COMPOSITE SPECTRA IN MERGING U/LIRGs CAUSED BY SHOCKS

We present a key result from our optical integral field spectroscopic survey of 27 nearby ultra luminous and luminous infrared galaxies (U/LIRGs) from the Great Observatory All-sky LIRG Survey. Using spatially resolved multi-component emission line fitting to trace the emission line ratios and velocity dispersion of the ionized gas, we quantify for the first time the widespread shock ionization in gas-rich merging U/LIRGs. Our results show a fractional contribution to the total observed \Ha flux from radiative shocks increasing from a few percent during early merger stages to upwards of 60% of the observed optical emission line flux in late stage mergers. We compare our resolved spectroscopy to nuclear spectra and find that 3/4 of the galaxies in our sample that would be classified as "composite" based on optical spectroscopy are primarily characterized by a combination of star formation and merger-driven shocks. Our results have important implications for the interpretation of "composite" rest-frame optical spectra of U/LIRGs as starburst+AGN, as the shock emission combined with star formation can mimic "composite" optical spectra in the absence of any contribution from an AGN.

AGN and Starbursts in Dusty Galaxy Mergers: Insights from the Great Observatories All-sky LIRG Survey

AbstractThe Great Observatories All-sky LIRG Survey (GOALS) is combining imaging and spectroscopic data from the Herschel, Spitzer, Hubble, GALEX, Chandra, and XMM-Newton space telescopes augmented with extensive ground-based observations in a multiwavelength study of approximately 180 Luminous Infrared Galaxies (LIRGs) and 20 Ultraluminous Infrared Galaxies (ULIRGs) that comprise a statistically complete subset of the 60μm-selected IRAS Revised Bright Galaxy Sample. The objects span the full range of galaxy environments (giant isolated spirals, wide and close pairs, minor and major mergers, merger remnants) and nuclear activity types (Seyfert 1, Seyfert 2, LINER, starburst/HII), with proportions that depend strongly on the total infrared luminosity. I will review the science motivations and present highlights of recent results selected from over 25 peer-reviewed journal articles published recently by the GOALS Team. Statistical investigations include detection of high-ionization Fe K emission indicative of deeply embedded AGN, comparison of UV and far-IR properties, investigations of the fraction of extended emission as a function of wavelength derived from mid-IR spectroscopy, mid-IR spectral diagnostics and spectral energy distributions revealing the relative contributions of AGN and starbursts to powering the bolometric luminosity, and quantitative structure analyses that delineate the evolution of stellar bars and nuclear stellar cusps during the merger process. Multiwavelength dissections of individual systems have unveiled large populations of young star clusters and heavily obscured AGN in early-stage (II Zw 96), intermediate-stage (Mrk 266, Mrk 273), and late-stage (NGC 2623, IC 883) mergers. A recently published study that matches numerical simulations to the observed morphology and gas kinematics in mergers has placed four systems on a timeline spanning 175-260 million years after their first passages, and modeling of additional (U)LIRGs is underway. A very brief description of ongoing work with Herschel and ALMA will be given. The talk will conclude with a discussion of the demographics of dual AGN (kpc-scale orbits) in the GOALS sample, including the difficulty of their detection and confirmation, a proposed sequence representing a progression from dual AGN to binary AGNs (sub-pc scale orbits), and implications for the scarcity of confirmed binary QSOs in recent large surveys. Grant support from NASA is gratefully acknowledged.

The build-up of nuclear stellar cusps in extreme starburst galaxies and major mergers

Nuclear stellar cusps are defined as central excess light component in the stellar light profiles of galaxies and are suggested to be stellar relics of intense compact starbursts in the central ∼100–500 pc region of gas-rich major mergers. Here, we probe the build-up of nuclear cusps during the actual starburst phase for a complete sample of luminous infrared galaxy (LIRG) systems (85 LIRGs, with 11.4 < log [LIR/L⊙] < 12.5) in the Great Observatories All-sky LIRG Survey sample. Cusp properties are derived via 2D fitting of the nuclear stellar light imaged in the near-infrared (NIR) by the Hubble Space Telescope and have been combined with mid-infrared (IR) diagnostics for active galactic nucleus (AGN)/starburst characterization. We find that nuclear stellar cusps are resolved in 76 per cent of LIRGs (merger and non-interacting galaxies). The cusp strength and luminosity increase with far-IR luminosity (excluding AGN) and merger stage, confirming theoretical models that starburst activity is associated with the build-up of nuclear stellar cusps. Evidence for ultracompact nuclear starbursts is found in ∼13 per cent of LIRGs, which have a strong unresolved central NIR light component but no significant contribution of an AGN. The nuclear NIR surface density (measured within 1 kpc radius) increases by a factor of ∼5 towards late merger stages. A careful comparison to local early-type galaxies with comparable masses reveals (a) that local (U)LIRGs have a significantly larger cusp fraction and (b) that the majority of the cusp LIRGs have host galaxy luminosities (H band) similar to core ellipticals which are roughly one order in magnitude larger than those for cusp ellipticals.

HUBBLE SPACE TELESCOPEACS IMAGING OF THE GOALS SAMPLE: QUANTITATIVE STRUCTURAL PROPERTIES OF NEARBY LUMINOUS INFRARED GALAXIES WITHLIR&gt; 1011.4L☉

A {\it Hubble Space Telescope} ({\it HST}) / Advanced Camera for Surveys (ACS) study of the structural properties of 85 luminous and ultraluminous ($L_{\rm IR} > 10^{11.4}$ L$_\odot$) infrared galaxies (LIRGs and ULIRGs) in the Great Observatories All-sky LIRG Survey (GOALS) sample is presented. Two-dimensional GALFIT analysis has been performed on F814W "{\it I}-band" images to decompose each galaxy, as appropriate, into bulge, disk, central PSF and stellar bar components. The fraction of bulge-less disk systems is observed to be higher in LIRGs (35%) than in ULIRGs (20%), with the disk+bulge systems making up the dominant fraction of both LIRGs (55%) and ULIRGs (45%). Further, bulge+disk systems are the dominant late-stage merger galaxy type and are the dominant type for LIRGs and ULIRGs at almost every stage of galaxy-galaxy nuclear separation. The mean {\it I}-band host absolute magnitude of the GOALS galaxies is $-22.64\pm$0.62 mag (1.8$^{+1.4}_{-0.4}$ L$^*_I$), and the mean bulge absolute magnitude in GOALS galaxies is about 1.1 magnitude fainter than the mean host magnitude. Almost all ULIRGs have bulge magnitudes at the high end (-20.6 to - 23.5 mag) of the GOALS bulge magnitude range. Mass ratios in the GOALS binary systems are consistent with most of the galaxies being the result of major mergers, and an examination of the residual-to-host intensity ratios in GOALS binary systems suggests that smaller companions suffer more tidal distortion than the larger companions. We find approximately twice as many bars in GOALS disk+bulge systems (32.8%) than in pure-disk mergers (15.9%) but most of the disk+bulge systems that contain bars are disk-dominated with small bulges. The bar-to-host intensity ratio, bar half-light radius, and bar ellipticity in GOALS galaxies are similar to those found in nearby spiral galaxies.(abridged)

MID-INFRARED PROPERTIES OF NEARBY LUMINOUS INFRARED GALAXIES. I. SPITZER INFRARED SPECTROGRAPH 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.

SPECTRAL ENERGY DISTRIBUTIONS OF LOCAL LUMINOUS AND ULTRALUMINOUS INFRARED GALAXIES

Luminous and ultraluminous infrared galaxies ((U)LIRGs) are the most extreme star forming galaxies in the universe. The local (U)LIRGs provide a unique opportunity to study their multi-wavelength properties in detail for comparison to their more numerous counterparts at high redshifts. We present common large aperture photometry at radio through X-ray wavelengths, and spectral energy distributions (SEDs) for a sample of 53 nearby LIRGs and 11 ULIRGs spanning log (LIR/Lsun) = 11.14-12.57 from the flux-limited Great Observatories All-sky LIRG Survey (GOALS). The SEDs for all objects are similar in that they show a broad, thermal stellar peak and a dominant FIR thermal dust peak, where nuLnu(60um) / nuLnu(V) increases from ~2-30 with increasing LIR. When normalized at IRAS-60um, the largest range in the luminosity ratio, R(lambda)=log[nuLnu(lambda)/nuLnu(60um)] observed over the full sample is seen in the Hard X-rays (HX=2-10 keV). A small range is found in the Radio (1.4GHz), where the mean ratio is largest. Total infrared luminosities, LIR(8-1000um), dust temperatures, and dust masses were computed from fitting thermal dust emission modified blackbodies to the mid-infrared (MIR) through submillimeter SEDs. The new results reflect an overall ~0.02 dex lower luminosity than the original IRAS values. Total stellar masses were computed by fitting stellar population synthesis models to the observed near-infrared (NIR) through ultraviolet (UV) SEDs. Mean stellar masses are found to be log(M/Msun) = 10.79+/-0.40. Star formation rates have been determined from the infrared (SFR_IR~45Msun/yr) and from the monochromatic UV luminosities (SFR_UV~1.3Msun/yr), respectively. Multiwavelength AGN indicators have be used to select putative AGN: about 60% of the ULIRGs would have been classified as an AGN by at least one of the selection criteria.

Probing the Build-Up of Stellar Mass in the Center of IR Luminous Major Mergers with HST

AbstractInteractions and mergers are important drivers of galaxy evolution, transform spiral galaxies into massive ellipticals, and fuel both powerful starbursts and massive nuclear black holes. In particular one galaxy population, namely Luminous Infrared Galaxies (LIRGs), are believed to be responsible for most of the star formation that happened in the history of the universe (see e.g. Le Floch et al. 2005, Caputi et al. 2007, Magnelli et al. 2009), and hence represent a critical phase in the evolution of galaxies where most of the galaxies mass is building up. During a merger process, violent relaxation acts on stars present in gas-rich progenitor disks, while the centers are structured by the relics of dissipational, compact starbursts, imprinting a central “extra light” component or “cusp” into the surface brightness profiles of merger remnants. Our HST NICMOS/WFC3 imaging program of the 88 most luminous LIRGs in the Great Observatories Allsky LIRG Survey (GOALS, see Armus et al. 2009) shows that the central luminosity surface density in nearby LIRGs increases significantly along the merger sequence, indicating that the gas inflow fuels a central starburst and subsequently builds a compact stellar cusp (Haan et al. 2011). A large fraction of all galaxies in our sample possess double or multiple nuclei (~63%). Half of these double nuclei are not visible in the HST B-band images due to dust obscuration, which implies strong limitations on the ability to detect the true nuclear structures of luminous infrared galaxies at high-redshift (z &gt;2) and may explain some of the apparent discrepancy of the LIRG population and merger ratio between local and high-redshift galaxies. We find that ULIRGs (log[LIR/L⊙] &gt; 12.0) have significantly smaller nuclear separations than LIRGs (log[LIR/L⊙] = 11.4 — 12.0) with a median value of 1.2 kpc and 6.7 kpc, respectively. In our sample, merger (regardless of whether LIRG or ULIRG) seem to be prevalent at two time scales (based on the projected nuclear separation and mass ratio of the nuclei): First, at a remaining merger time scale of 0.3&lt;[t-tmerg]&lt;1.3 Gyr (53% of mergers in our sample), and second, at [t-tmerg]~ 0 (26%), likely representing the first passage of interacting galaxies and the final nuclear coalescence, respectively, with a post-merger time (starburst phase after the nuclei merged) of roughly 300 Myrs.

MULTI-WAVELENGTH GOALS OBSERVATIONS OF STAR FORMATION AND ACTIVE GALACTIC NUCLEUS ACTIVITY IN THE LUMINOUS INFRARED GALAXY IC 883

New optical HST, Spitzer, GALEX, and Chandra observations of the single-nucleus, luminous infrared galaxy (LIRG) merger IC 883 are presented. The galaxy is a member of the Great Observatories All-sky LIRG Survey (GOALS), and is of particular interest for a detailed examination of a luminous late-stage merger due to the richness of the optically-visible star clusters and the extended nature of the nuclear X-ray, mid-IR, CO and radio emission. In the HST ACS images, the galaxy is shown to contain 156 optically visible star clusters distributed throughout the nuclear regions and tidal tails of the merger, with a majority of visible clusters residing in an arc ~ 3-7 kpc from the position of the mid-infrared core of the galaxy. The luminosity functions of the clusters have an alpha_F435W ~ -2.17+/-0.22 and alpha_F814W ~ -2.01+/-0.21. Further, the colors and absolute magnitudes of the majority of the clusters are consistent with instantaneous burst population synthesis model ages in the range of a few x10^7 - 10^8 yrs (for 10^5 M_sun clusters), but may be as low as few x10^6 yrs with extinction factored in. The X-ray and mid-IR spectroscopy are indicative of predominantly starburst-produced nuclear emission, and the star formation rate is ~ 80 M_sun / yr. The kinematics of the CO emission and the morphology of both the CO and radio emission are consistent with the nuclear starburst being situated in a highly inclined disk 2 kpc in diameter with an infrared surface brightness mu_IR ~ 2x10^11 L_sun kpc^-2, a factor of 10 less than that of the Orion star-forming region. Finally, the detection of the [Ne V] 14.32 um emission line is evidence that an AGN is present. The faintness of the line (i.e., [Ne V] / [Ne II] um ~ 0.01) and the small equivalent width of the 6.2 um PAH feature ($= 0.39\mu$m) are both indicative of a relatively weak AGN. (abridged)

THE SPATIAL EXTENT OF (U)LIRGS IN THE MID-INFRARED. II. FEATURE EMISSION

We present results from the second part of our analysis of the extended mid-infrared (MIR) emission of the Great Observatories All-Sky LIRG Survey (GOALS) sample based on 5-14 micron low-resolution spectra obtained with the IRS on Spitzer. We calculate the fraction of extended emission as a function of wavelength for all galaxies in the sample, FEE_lambda, and spatially separate the MIR spectrum of galaxies into their nuclear and extended components. We find that the [NeII] emission line is as compact as the hot dust MIR continuum, while the polycyclic aromatic hydrocarbon (PAH) emission is more extended. The 6.2 and 7.7 micron PAH emission is more compact than that of the 11.3 micron PAH, which is consistent with the formers being enhanced in a more ionized medium. The presence of an AGN or a powerful nuclear starburst increases the compactness of the hot dust MIR continuum, but has a negligible effect on the spatial extent of the PAH emission on kpc-scales. Globally, the spectra of the extended emission component are homogeneous for all galaxies in GOALS. This suggests that the physical properties of star formation taking place at distances farther than 1.5 kpc from the nuclei of (U)LIRGs are very similar, resembling local star-forming galaxies with L_IR < 10^11 Lsun, as well as star formation-dominated ULIRGs at z~2. In contrast, the MIR spectra of the nuclear component of local (U)LIRGs are very diverse. This implies that the observed variety of their integrated MIR properties arise, on average, only from the processes that are taking place in their cores.

THE NUCLEAR STRUCTURE IN NEARBY LUMINOUS INFRARED GALAXIES:HUBBLE SPACE TELESCOPENICMOS IMAGING OF THE GOALS SAMPLE

We present results of Hubble Space Telescope NICMOS H-band imaging of 73 of most luminous (i.e., log[L_IR/L_0]>11.4) Infrared Galaxies (LIRGs) in the Great Observatories All-sky LIRG Survey (GOALS). This dataset combines multi-wavelength imaging and spectroscopic data from space (Spitzer, HST, GALEX, and Chandra) and ground-based telescopes. In this paper we use the high-resolution near-infrared data to recover nuclear structure that is obscured by dust at optical wavelengths and measure the evolution in this structure along the merger sequence. A large fraction of all galaxies in our sample possess double nuclei (~63%) or show evidence for triple nuclei (~6%). Half of these double nuclei are not visible in the HST B-band images due to dust obscuration. The majority of interacting LIRGs have remaining merger timescales of 0.3 to 1.3 Gyrs, based on the projected nuclear separations and the mass ratio of nuclei. We find that the bulge luminosity surface density increases significantly along the merger sequence (primarily due to a decrease of the bulge radius), while the bulge luminosity shows a small increase towards late merger stages. No significant increase of the bulge Sersic index is found. LIRGs that show no interaction features have on average a significantly larger bulge luminosity, suggesting that non merging LIRGs have larger bulge masses than merging LIRGs. This may be related to the flux limited nature of the sample and the fact that mergers can significantly boost the IR luminosity of otherwise low luminosity galaxies. We find that the projected nuclear separation is significantly smaller for ULIRGs (median value of 1.2 kpc) than for LIRGs (mean value of 6.7 kpc), suggesting that the LIRG phase appears earlier in mergers than the ULIRG phase.

Modeling IR spectral energy distributions: a pilot study of starburst parameters and silicate absorption curves for some GOALS galaxies

This paper describes a pilot study into the spectral energy distribution (SED) fitting and the derivation of physical parameters for 19 galaxies observed as part of the Great Observatories All-sky LIRG Survey (GOALS) survey as observed with the Spitzer Space Telescope. For this we have used the pan-spectral fitting tools developed in a series of papers by Dopita and his co-workers. We show that the standard Lee and Draine ‘astronomical silicate’ model cannot provide a good fit to the silicate absorption features as observed in the heavily dust-extinguished (A V∼50 mag.) starbursts. We have derived an empirical fit to the ‘starburst silicate’ absorption in these objects. This absorption curve is consistent with the silicate grains being systematically larger in starburst environments than in the local Galactic interstellar medium. We demonstrate the sensitivity of the SED fitting to each of the fitted parameters, and derive these parameters for those galaxies which do not have an embedded AGN. This technique is simple and provides reasonably robust and uniform parameters for the starburst, especially as far as the star formation rate, population of old stars, compactness of the starburst region and total foreground extinction are concerned. However, the chemical abundances and the optical extinction cannot be reliably determined by this analysis, and optical SEDs will also be required to provide a complete characterization of the starburst region and of the surrounding galaxy.