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-IR Spitzer IRS spectra from 5-38um of 244 LIRG nuclei. The detailed fits and high quality spectra allow for characterization of the individual PAH features, warm molecular hydrogen emission, and optical depths for silicate dust grains and water ices. We find that starbursting LIRGs, which make up the majority of GOALS, are very consistent in their MIR properties (i.e. tau_9.7um, tau_ice, neon line and PAH feature ratios). However, as their PAH EQW decreases, usually an indicator of an increasingly dominant AGN, LIRGs cover a larger spread in these MIR parameters. The contribution from PAHs to the total L(IR) in LIRGs varies from 2-29% and LIRGs prior to their first encounter show higher L(PAH)/L(IR) ratios on average. We observe a correlation between the strength of the starburst (IR8) and the PAH fraction at 8um but not with the 7.7 to 11.3 PAH ratio, suggesting the fractional 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 obscured sources (s_9.7um < -1.24). Ice absorption features are observed in 11% (56%) of GOALS LIRGs (ULIRGs). 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%) show excess H2 emission that is inconsistent with PDR models and may be excited by shocks or AGN-induced outflows.Comment: 21 pages with 20 figures plus 2 table

Spitzer Quasar and ULIRG Evolution Study (QUEST). IV. Comparison of 1-Jy Ultraluminous Infrared Galaxies with Palomar-Green Quasars

We report the results from a comprehensive study of 74 ultraluminous infrared galaxies (ULIRGs) and 34 Palomar-Green (PG) quasars within z ~ 0.3$ observed with the Spitzer Infrared Spectrograph (IRS). The contribution of nuclear activity to the bolometric luminosity in these systems is quantified using six independent methods that span a range in wavelength and give consistent results within ~ +/-10-15% on average. The average derived AGN contribution in ULIRGs is ~35-40%, ranging from ~15-35% among "cool" (f_25/f_60 =< 0.2) optically classified HII-like and LINER ULIRGs to ~50 and ~75% among warm Seyfert 2 and Seyfert 1 ULIRGs, respectively. This number exceeds ~80% in PG QSOs. ULIRGs fall in one of three distinct AGN classes: (1) objects with small extinctions and large PAH equivalent widths are highly starburst-dominated; (2) systems with large extinctions and modest PAH equivalent widths have larger AGN contributions, but still tend to be starburst-dominated; and (3) ULIRGs with both small extinctions and small PAH equivalent widths host AGN that are at least as powerful as the starbursts. The AGN contributions in class 2 ULIRGs are more uncertain than in the other objects, and we cannot formally rule out the possibility that these objects represent a physically distinct type of ULIRGs. A morphological trend is seen along the sequence (1)-(2)-(3), in general agreement with the standard ULIRG - QSO evolution scenario and suggestive of a broad peak in extinction during the intermediate stages of merger evolution. However, the scatter in this sequence, implies that black hole accretion, in addition to depending on the merger phase, also has a strong chaotic/random component, as in local AGN. (abridged)Comment: 61 pages, 39 figures, 16 tables, accepted for publication in ApJS, June 2009 issue. Unabbreviated version can be found at http://www.astro.umd.edu/~veilleux/pubs/quest4.pd

Mid-Infrared Spectral Diagnostics of Luminous Infrared Galaxies

We present a statistical analysis of the mid-infrared (MIR) spectra of 248 luminous infrared (IR) galaxies (LIRGs) which comprise the Great Observatories All-sky LIRG Survey (GOALS) observed with the Infrared Spectrograph (IRS) on-board the Spitzer Space Telescope. The GOALS sample enables a direct measurement of the relative contributions of star-formation and active galactic nuclei (AGN) to the total IR emission from a large sample of local LIRGs. The AGN contribution to the MIR emission (f-AGN) is estimated by employing several diagnostics based on the properties of the [NeV], [OIV] and [NeII] fine structure gas emission lines, the 6.2 microns PAH and the shape of the MIR continuum. We find that 18% of all LIRGs contain an AGN and that in 10% of all sources the AGN contributes more than 50% of the total IR luminosity. Summing up the total IR luminosity contributed by AGN in all our sources suggests that AGN supply ~12% of the total energy emitted by LIRGs. The average spectrum of sources with an AGN looks similar to the average spectrum of sources without an AGN, but it has lower PAH emission and a flatter MIR continuum. AGN dominated LIRGs have higher IR luminosities, warmer MIR colors and are found in interacting systems more often than pure starbursts LIRGs. However we find no linear correlations between these properties and f-AGN. We used the IRAC colors of LIRGs to confirm that finding AGN on the basis of their MIR colors may miss ~40% of AGN dominated (U)LIRGsComment: accepted for publication in ApJ, 34 pages, 12 figure

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

Investigation of Dual Active Nuclei, Outflows, Shock-Heated Gas, and Young Star Clusters in Markarian 266

Results of observations with the Spitzer, Hubble, GALEX, Chandra, and XMM-Newton space telescopes are presented for the Luminous Infrared Galaxy (LIRG) merger Mrk 266. The SW (Seyfert 2) and NE (LINER) nuclei reside in galaxies with Hubble types SBb (pec) and S0/a (pec), respectively. Both galaxies have L > L*, and they are inferred to each contain a ~2.5x10^8 M_sun black hole. Mrk 266 SW is likely the primary source of a bright Fe K-alpha line detected from the system, consistent with the reflection-dominated X-ray spectrum of a heavily obscured AGN. Optical knots embedded in an arc with aligned radio continuum radiation, combined with luminous H_2 line emission, provide evidence for a radiative bow shock in an AGN-driven outflow surrounding the NE nucleus. Soft X-ray emission modeled as shock-heated plasma is co-spatial with radio continuum emission between the galaxies. Mid-IR diagnostics indicate roughly equal contributions of AGN and starburst radiation powering the bolometric luminosity. Approximately 120 star clusters have been detected, with most having estimated ages < 50 Myr. Detection of 24 micron emission extending ~34 arcsec (20 kpc) north of the galaxies is interpreted as ~2x10^7 M_sun of dust entrained in an outflowing superwind; at optical wavelengths this region is resolved into a fragmented morphology indicative of Rayleigh-Taylor instabilities in an expanding shell of ionized gas. Mrk 266 demonstrates that the dust "blow-out" phase can begin in a LIRG well before the galaxies fully coalesce during a subsequent ULIRG phase, and rapid gas consumption in luminous dual AGNs (kpc scale separations) early in the merger process may explain the paucity of detected binary QSOs (sub-pc scale orbits) in large surveys. An evolutionary sequence is proposed representing a progression from dual to binary AGNs, accompanied by an increase in observed L_x/L_ir ratios by 10^4 or more.Comment: Accepted for publication in the Astronomical Journal, 40 pages, 31 figures, 15 table

HERUS: the far-IR/submm spectral energy distributions of local ULIRGs and photometric atlas

We present the Herschel-SPIRE photometric atlas for a complete flux limited sample of 43 local ultraluminous infrared galaxies (ULIRGs), selected at 60 μm by IRAS, as part of the HERschel ULIRG Survey (HERUS). Photometry observations were obtained using the SPIRE instrument at 250, 350, and 500 μm. We describe these observations, present the results, and combine the new observations with data from IRAS to examine the far-infrared spectral energy distributions (SEDs) of these sources. We fit the observed SEDs of HERUS objects with a simple parametrized modified blackbody model, where temperature and emissivity β are free parameters. We compare the fitted values to those of non-ULIRG local galaxies, and find, in agreement with earlier results, that HERUS ULIRGs have warmer dust (median temperature T = 37.9 ± 4.7 K compared to 21.3 ± 3.4 K) but a similar β distribution (median β = 1.7 compared to 1.8) to the Herschel reference sample (HRS, Cortese et al. 2014) galaxies. Dust masses are found to be in the range of 107.5–109 M⊙, significantly higher than that of HRS sources. We compare our results for local ULIRGs with higher redshift samples selected at 250 and 850 μm. These latter sources generally have cooler dust and/or redder 100-to-250  μm colours than our 60 μm-selected ULIRGs. We show that this difference may in part be the result of the sources being selected at different wavelengths rather than being a simple indication of rapid evolution in the properties of the population

Mid-infrared spectroscopy of infrared-luminous galaxies at z~0.5-3

We present results on low-resolution mid-infrared (MIR) spectra of 70 infrared-luminous galaxies obtained with the Infrared Spectrograph (IRS) onboard Spitzer. We selected sources from the European Large Area Infrared Survey (ELAIS) with S15 > 0.8 mJy and photometric or spectroscopic z > 1. About half of the sample are QSOs in the optical, while the remaining sources are galaxies, comprising both obscured AGN and starbursts. We classify the spectra using well-known infrared diagnostics, as well as a new one that we propose, into three types of source: those dominated by an unobscured AGN (QSOs), obscured AGN, and starburst-dominated sources. Starbursts concentrate at z ~ 0.6-1.0 favored by the shift of the 7.7-micron PAH band into the selection 15 micron band, while AGN spread over the 0.5 < z < 3.1 range. Star formation rates (SFR) are estimated for individual sources from the luminosity of the PAH features. An estimate of the average PAH luminosity in QSOs and obscured AGN is obtained from the composite spectrum of all sources with reliable redshifts. The estimated mean SFR in the QSOs is 50-100 Mo yr^-1, but the implied FIR luminosity is 3-10 times lower than that obtained from stacking analysis of the FIR photometry, suggesting destruction of the PAH carriers by energetic photons from the AGN. The SFR estimated in obscured AGN is 2-3 times higher than in QSOs of similar MIR luminosity. This discrepancy might not be due to luminosity effects or selection bias alone, but could instead indicate a connection between obscuration and star formation. However, the observed correlation between silicate absorption and the slope of the near- to mid-infrared spectrum is compatible with the obscuration of the AGN emission in these sources being produced in a dust torus.Comment: 32 pages, 24 figures, 15 tables, accepted for publication in MNRA

Neural cognitive control moderates the association between insular risk processing and risk-taking behaviors via perceived stress in adolescents

Adolescence is a critical period for the initiation of risk-taking behaviors. We examined the longitudinal interplay between neural correlates of risk processing and cognitive control in predicting risk-taking behaviors via stress. The sample consisted of 167 adolescents (53% males) who were assessed twice (MAgeTime1 = 14.13, MAgeTime2 = 15.05). Neural risk processing was operationalized as blood-oxygen-level-dependent (BOLD) responses in the anterior insula during a lottery choice task and neural cognitive control as BOLD responses during an inhibitory control task. Adolescents reported on perceived stress and risk-taking behaviors. Structural equation modeling analyses indicated that low insular risk processing predicted increases in perceived stress, while perceived stress did not predict changes in insular risk processing across one year. Moreover, significant moderation by neural cognitive control indicated that low insular risk processing predicted increases in risk-taking behaviors via increases in perceived stress among adolescents with poor neural cognitive control, but not among adolescents with good neural cognitive control. The results suggest that risk processing in the anterior insular cortex plays an important role in stress experience and risk-taking behaviors particularly for vulnerable adolescents with poor neural cognitive control

Optical identification of XMM sources in the CFHTLS

We present optical spectroscopic identifications of X-ray sources in ~3 square degrees of the XMM-Large Scale Structure survey (XMM-LSS), also covered by the Canada France Hawaii Telescope Legacy Survey (CFHTLS), obtained with the AAOmega instrument at the Anglo Australian Telescope. In a flux limited sample of 829 point like sources in the optical band with g' <~22 mag and the 0.5-2 keV flux > 1x10^{-15}erg/cm^2/s, we observed 695 objects and obtained reliable spectroscopic identification for 489 sources, ~59% of the overall sample. We therefore increase the number of identifications in this field by a factor close to five. Galactic stellar sources represent about 15% of the total (74/489). About 55% (267/489) are broad-line Active Galactic Nuclei (AGNs) spanning redshifts between 0.15 and 3.87 with a median value of 1.68. The optical-to-X-ray spectral index of the broad-line AGNs is 1.47, typical of optically-selected Type I quasars and is found to correlate with the rest frame X-ray and optical monochromatic luminosities at 2 keV and 2500 angstroms respectively. Consistent with previous studies, we find alpha_ox not to be correlated with z. In addition, 32 and 116 X-ray sources are, respectively absorption and emission-line galaxies at z<0.76. From a line ratio diagnostic diagram it is found that in about 50% of these emission line galaxies, the emission lines are powered significantly by the AGN. Thirty of the XMM sources are detected at one or more radio frequencies. In addition, 24 sources have ambiguous identification: in 8 cases, two XMM sources have a single optical source within 6 arcsecs of each of them, whereas, 2 and 14 XMM sources have, respectively, 3 and 2 possible optical sources within 6 arcsecs of each of them.Comment: 15 pages, 14 figures, 5 tables, accepted for publication in MNRA