Showing 1–50 of 77 results for author: Spoon, H W W

Spectroscopic Constraints on the Mid-Infrared Attenuation Curve: I -- Attenuation Model using PAH Emission

Authors: Thomas S. -Y. Lai, J. D. T. Smith, Els Peeters, Henrik W. W. Spoon, Shunsuke Baba, Masatoshi Imanishi, Takao Nakagawa

Abstract: We introduce a novel model to spectroscopically constrain the mid-infrared (MIR) extinction/attenuation curve from 3--17 um, using Polycyclic Aromatic Hydrocarbon (PAH) emission drawn from an AKARI-Spitzer extragalactic cross-archival dataset. Currently proposed MIR extinction curves vary significantly in their slopes toward the near-infrared, and the variation of the strengths and shapes of the 9… ▽ More We introduce a novel model to spectroscopically constrain the mid-infrared (MIR) extinction/attenuation curve from 3--17 um, using Polycyclic Aromatic Hydrocarbon (PAH) emission drawn from an AKARI-Spitzer extragalactic cross-archival dataset. Currently proposed MIR extinction curves vary significantly in their slopes toward the near-infrared, and the variation of the strengths and shapes of the 9.7 um and 18 um silicate absorption features make MIR spectral modeling and interpretation challenging, particularly for heavily obscured galaxies. By adopting the basic premise that PAH bands have relatively consistent intrinsic ratios within dusty starbursting galaxies, we can, for the first time, empirically determine the overall shape of the MIR attenuation curve by measuring the differential attenuation at specific PAH wavelengths. Our attenuation model shows PAH emission in most (U)LIRGs is unambiguously subjected to attenuation, and we find strong evidence that PAH bands undergo differential attenuation as obscuration increases. Compared to pre-existing results, the MIR attenuation curve derived from the model favors relatively gray continuum absorption from 3-8 $μ$m and silicate features with intermediate strength at 9.7 um but with stronger than typical 18 um opacity. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: 15 pages, 5 figures. Paper re-submitted to ApJ after addressing the first round of comments by the referee

A Detailed Look at the Most Obscured Galactic Nuclei in the Mid-Infrared

Authors: F. R. Donnan, D. Rigopoulou, I. García-Bernete, M. Pereira-Santaella, A. Alonso-Herrero, P. F. Roche, S. Aalto, A. Hernán-Caballero, H. W. W. Spoon

Abstract: Context. Compact Obscured Nuclei (CONs) are an extreme phase of galaxy evolution where rapid supermassive black hole growth and$/$or compact star-forming activity is completely obscured by gas and dust. Aims. We investigate the properties of CONs in the mid-infrared and explore techniques aimed at identifying these objects such as through the equivalent width (EW) ratios of their Polycyclic Aromat… ▽ More Context. Compact Obscured Nuclei (CONs) are an extreme phase of galaxy evolution where rapid supermassive black hole growth and$/$or compact star-forming activity is completely obscured by gas and dust. Aims. We investigate the properties of CONs in the mid-infrared and explore techniques aimed at identifying these objects such as through the equivalent width (EW) ratios of their Polycyclic Aromatic Hydrocarbon (PAH) features. Methods. We model Spitzer spectra by decomposing the continua into nuclear and star-forming components from which we then measure the nuclear optical depth, $τ_N$, of the $9.8 μ$m silicate absorption feature. We also use Spitzer spectral maps to investigate how PAH EW ratios vary with aperture size for objects hosting CONs. Results. We find that the nuclear optical depth, $τ_N$, strongly correlates with the HCN-vib emission line in the millimetre for CONs with a Pearson correlation coefficient of 0.91. We find the PAH EW ratios technique to be effective at selecting CONs and robust against highly inclined galaxies where strong dust lanes may mimic a CON like spectrum by producing a high $τ_N$. Our analysis of the Spitzer spectral maps showed that the efficacy of the PAH EW ratios to isolate CONs is reduced when there is a strong star-forming component from the host galaxy. In addition, we find that the use of the inferred nuclear optical depth is a reliable method to identify CONs in $36^{+8}_{-7}\%$ of ULIRGs and $17^{+3}_{-3}\%$ of LIRGs, consistent with previous work. Conclusions. We confirm mid-IR spectra to be a powerful diagnostic of CONs where the increased sensitivity of JWST will allow identification of CONs at cosmic noon revealing this extreme but hidden phase of galaxy evolution. △ Less

Submitted 17 November, 2022; originally announced November 2022.

Comments: Accepted for publication in A&A. 23 pages, 15 figures, 3 tables

Journal ref: A&A 669, A87 (2023)

The Obscured Nucleus and Shocked Environment of VV 114E Revealed by JWST/MIRI Spectroscopy

Authors: F. R. Donnan, I. García-Bernete, D. Rigopoulou, M. Pereira-Santaella, A. Alonso-Herrero, P. F. Roche, A. Hernán-Caballero, H. W. W. Spoon

Abstract: Compact Obscured Nuclei (CONs) potentially hide extreme supermassive black hole (SMBH) growth behind large column densities of gas/dust. We present a spectroscopic analysis of the heavily obscured nucleus and the surrounding environment of the eastern region of the nearby ($z = 0.02007$) interacting galaxy VV 114 with the JWST Mid-InfraRed Instrument (MIRI). We model the spectrum from 4.9 - 28… ▽ More Compact Obscured Nuclei (CONs) potentially hide extreme supermassive black hole (SMBH) growth behind large column densities of gas/dust. We present a spectroscopic analysis of the heavily obscured nucleus and the surrounding environment of the eastern region of the nearby ($z = 0.02007$) interacting galaxy VV 114 with the JWST Mid-InfraRed Instrument (MIRI). We model the spectrum from 4.9 - 28 $μ$m to extract Polycyclic Aromatic Hydrocarbon (PAH) emission and the underlying obscured continuum. We find that the NE nucleus (A) is highly obscured where the low PAH equivalent width (EW) ratio, EW(12.7)/EW(11.3), reveals a dust enshrouded continuum source. This is confirmed by decomposing the continuum into nuclear and star-forming where the nuclear component is found to be typical of CONs. The 11.3/6.2 PAH flux ratio is consistent with originating in star-forming regions rather than typical AGN. The second nucleus (B) is much less obscured, with PAH flux ratios also typical of star-forming regions. We do not detect any high ionisation lines such as [Ne V] or [Ne VI] which suggests that if an AGN is present it must be highly obscured. Additionally, we detect a shock front south of the secondary nucleus (B) in the [Fe II] (5.34 $μ$m) line and in warm molecular hydrogen. The 6.2 PAH emission does not spatially coincide with the low-J transitions of H$_2$ but rather appears strong at the shock front which may suggest destruction of the ionised PAHs in the post-shock gas behind the shock front. △ Less

Submitted 15 December, 2022; v1 submitted 10 October, 2022; originally announced October 2022.

Comments: 15 pages, 15 figures. Accepted for publication in MNRAS

A new look at local ultraluminous infrared galaxies: the atlas and radiative transfer models of their complex physics

Authors: A. Efstathiou, D. Farrah, J. Afonso, D. L. Clements, E. González-Alfonso, M. Lacy, S. Oliver, V. Papadopoulou Lesta, C. Pearson, D. Rigopoulou, M. Rowan-Robinson, H. W. W. Spoon, A. Verma, L. Wang

Abstract: We present the ultraviolet to submillimetre spectral energy distributions (SEDs) of the HERschel Ultra Luminous Infrared Galaxy Survey (HERUS) sample of 42 local ultraluminous infrared galaxies (ULIRGs) and fit them with a Markov chain Monte Carlo (MCMC) code using the CYprus models for Galaxies and their NUclear Spectra (CYGNUS) radiative transfer models for starbursts, active galactic nucleus (A… ▽ More We present the ultraviolet to submillimetre spectral energy distributions (SEDs) of the HERschel Ultra Luminous Infrared Galaxy Survey (HERUS) sample of 42 local ultraluminous infrared galaxies (ULIRGs) and fit them with a Markov chain Monte Carlo (MCMC) code using the CYprus models for Galaxies and their NUclear Spectra (CYGNUS) radiative transfer models for starbursts, active galactic nucleus (AGN) tori and host galaxy. The Spitzer IRS spectroscopy data are included in the fitting. Our bayesian SED fitting method takes comparable time to popular energy balance methods but it is more physically motivated and versatile. All HERUS galaxies harbor high rates of star formation but we also find bolometrically significant AGN in all of the galaxies of the sample. We estimate the correction of the luminosities of the AGN in the ULIRGs due to the anisotropic emission of the torus and find that it could be up to a factor of $\sim10$ for nearly edge-on tori. We present a comparison of our results with the smooth torus model of Fritz et al. and the two-phase models of Siebenmorgen et al. and SKIRTOR. We find that the CYGNUS AGN torus models fit significantly better the SEDs of our sample compared to all other models. We find no evidence that strong AGN appear either at the beginning or end of a starburst episode or that starbursts and AGN affect each other. IRAS 01003-2238 and Mrk 1014 show evidence for dual AGN in their SED fits suggesting a minimum dual AGN fraction in the sample of 5%. △ Less

Submitted 29 April, 2022; originally announced May 2022.

Comments: Published in Monthly Notices of the Royal Astronomical Society. The CYGNUS models used in this paper are available at https://arc.euc.ac.cy/cygnus/

Journal ref: 2022MNRAS.512.5183E

A technique to select the most obscured galaxy nuclei

Authors: I. García-Bernete, D. Rigopoulou, S. Aalto, H. W. W. Spoon, A. Hernán-Caballero, A. Efstathiou, P. F. Roche, S. König

Abstract: Compact obscured nuclei (CONs) are mainly found in local U/LIRGs. In the local Universe, these sources are generally selected through the detection of the HCN-vib (3-2) emission line at submillimetre wavelengths. In this work, we present a diagnostic method to select deeply buried nuclei based on mid-infrared (mid-IR) polycyclic aromatic hydrocarbons (PAHs) and continuum ratios. Using Spitzer/IRS… ▽ More Compact obscured nuclei (CONs) are mainly found in local U/LIRGs. In the local Universe, these sources are generally selected through the detection of the HCN-vib (3-2) emission line at submillimetre wavelengths. In this work, we present a diagnostic method to select deeply buried nuclei based on mid-infrared (mid-IR) polycyclic aromatic hydrocarbons (PAHs) and continuum ratios. Using Spitzer/IRS spectra of a representative sample of local ULIRGs (z<0.27), we examine their PAH and underlying continuum emission ratios. For deeply embedded sources, we find that the 9.7 micron silicate absorption band has a particularly pronounced effect on the 11.3 micron PAH feature. The low flux level in the nuclear silicate absorption band enhances the 11.3 micron PAH feature contrast (high PAH equivalent width) compared to that of the other PAH features. The technique has been extended to include the use of the continuum ratios. However, the latter are affected both by the extinction coming from the host galaxy as well as the nuclear region, whereas the foreground extinction is cancelled out when using the PAH equivalent width ratios. We apply our method to the HERUS and GOALS samples and classify as CON candidates 14 ULIRGs and 10 LIRGs, corresponding to 30% of ULIRGs and 7% of LIRGs from these samples. We find that the observed continuum ratios of CON-dominated sources can be explained by assuming torus models with a tapered disk geometry and a smooth dust distribution. This suggests that the nuclear dusty structure of CONs has an extremely high dust coverage. We also demonstrate that the use of mid-IR color-color diagrams is an effective way to select CON-dominated sources at different redshifts. In particular, the combination of filters of the JWST/MIRI will enable the selection of CONs out to z~1.5. This will allow extending the selection of CONs to high redshifts where U/LIRGs are more numerous. △ Less

Submitted 28 March, 2022; originally announced March 2022.

Comments: 13 pages, 11 figures, 1 table, accepted for publication in A&A

Journal ref: A&A 663, A46 (2022)

The Infrared Database of Extragalactic Observables from Spitzer. -- II. The Database & The Diagnostic Power of Crystalline Silicate Features in Galaxy Spectra

Authors: H. W. W. Spoon, A. Hernán Caballero, D. Rupke, L. B. F. M. Waters, V. Lebouteiller, A. G. G. M. Tielens, T. Loredo, Y. Su, V. Viola

Abstract: We present the Infrared Database of Extragalactic Observables from Spitzer (IDEOS), a homogeneous, publicly available, database of 77 fitted mid-infrared observables in the 5.4-36um range, comprising measurements for 3335 galaxies observed in the low-resolution staring mode of the Infrared Spectrometer onboard the Spitzer Space Telescope. Among the included observables are PAH fluxes and their equ… ▽ More We present the Infrared Database of Extragalactic Observables from Spitzer (IDEOS), a homogeneous, publicly available, database of 77 fitted mid-infrared observables in the 5.4-36um range, comprising measurements for 3335 galaxies observed in the low-resolution staring mode of the Infrared Spectrometer onboard the Spitzer Space Telescope. Among the included observables are PAH fluxes and their equivalent widths, the strength of the 9.8um silicate feature, emission line fluxes, solid-state features, rest frame continuum fluxes, synthetic photometry, and a mid-infrared spectral classification. The IDEOS spectra were selected from the Cornell Atlas of Spitzer-IRS Sources. To our surprise we have detected at a >95% confidence level crystalline silicates in the spectra of 786 IDEOS galaxies. The detections range from single band detections to detections of all fitted crystalline bands (16, 19, 23, 28 and 33um). We find the strength of the crystalline silicate bands to correlate with the amorphous silicate strength, and the change from an emission to an absorption feature to occur at higher obscuration as the wavelength of the crystalline silicate band is longer. These observed characteristics are consistent with an origin for the amorphous and crystalline silicate features in a centrally heated dust geometry, either an edge-on disk or a cocoon. We find the 23 and 33um crystalline silicate bands to be well-suited to classify the obscuration level of galactic nuclei, even in the presence of strong circumnuclear star formation. Based on our detection statistics, we conclude that crystalline silicates are a common component of the interstellar medium of galactic nuclei. △ Less

Submitted 6 March, 2022; originally announced March 2022.

Comments: 46 pages, 34 figures, accepted for publication in ApJS. Associated spectra and data table can be downloaded from the IDEOS portal at http://ideos.astro.cornell.edu

All the PAHs: an AKARI-Spitzer Cross Archival Spectroscopic Survey of Aromatic Emission in Galaxies

Authors: Thomas S. -Y. Lai, J. D. T. Smith, Shunsuke Baba, Henrik W. W. Spoon, Masatoshi Imanishi

Abstract: We present a large sample of 2.5-38 $μm$ galaxy spectra drawn from a cross-archival comparison in the AKARI-Spitzer Extragalactic Spectral Survey (ASESS), and investigate a subset of 113 star-forming galaxies with prominent polycyclic aromatic hydrocarbon (PAH) emission spanning a wide range of star formation properties. With AKARI's extended 2.5-5 $μm$ wavelength coverage, we self-consistently mo… ▽ More We present a large sample of 2.5-38 $μm$ galaxy spectra drawn from a cross-archival comparison in the AKARI-Spitzer Extragalactic Spectral Survey (ASESS), and investigate a subset of 113 star-forming galaxies with prominent polycyclic aromatic hydrocarbon (PAH) emission spanning a wide range of star formation properties. With AKARI's extended 2.5-5 $μm$ wavelength coverage, we self-consistently model for the first time all PAH emission bands using a modified version of PAHFIT. We find $L_\mathrm{PAH\,3.3}/L_\mathrm{IR}$ $\sim$ 0.1% and the 3.3 $μm$ PAH feature contributes $\sim$1.5-3% to the total PAH power -- somewhat less than earlier dust models have assumed. We establish a calibration between 3.3 $μm$ PAH emission and star formation rate, but also find regimes where it loses reliability, including at high luminosity and low metallicity. The 3.4 $μm$ aliphatic emission and a broad plateau feature centered at 3.47 $μm$ are also modeled. As the shortest wavelength PAH feature, 3.3 $μm$ is susceptible to attenuation, leading to a factor of $\sim$ 3 differences in the inferred star formation rate at high obscuration with different assumed attenuation geometries. Surprisingly, $L_\mathrm{PAH\, 3.3} / L_\mathrm{Σ\, PAH}$ shows no sign of decline at high luminosities, and the low metallicity dwarf galaxy II Zw 40 exhibits an unusually strong 3.3 $μm$ band; both results suggest either that the smallest PAHs are better able to survive under intense radiation fields than presumed, or that PAH emission is shifted to shorter wavelengths in intense and high energy radiation environments. A photometric surrogate for 3.3 $μm$ PAH luminosity using JWST/NIRCam is provided and found to be highly reliable at low redshift. △ Less

Submitted 10 October, 2020; originally announced October 2020.

Comments: 29 pages, 18 figures, accepted in ApJ

Extinction in the 11.2 micron PAH band and the low L_11.2/L_IR in ULIRGs

Authors: Antonio Hernan-Caballero, Henrik W. W. Spoon, Almudena Alonso-Herrero, Evanthia Hatziminaoglou, Georgios E. Magdis, Pablo G. Perez-Gonzalez, Miguel Pereira-Santaella, Santiago Arribas, Isabella Cortzen, Alvaro Labiano, Javier Piqueras, Dimitra Rigopoulou

Abstract: We present a method for recovering the intrinsic (extinction-corrected) luminosity of the 11.2 micron PAH band in galaxy spectra. Using 105 high S/N Spitzer/IRS spectra of star-forming galaxies, we show that the equivalent width ratio of the 12.7 and 11.2 micron PAH bands is independent on the optical depth, with small dispersion of ~5% indicative of a nearly constant intrinsic flux ratio R_int =… ▽ More We present a method for recovering the intrinsic (extinction-corrected) luminosity of the 11.2 micron PAH band in galaxy spectra. Using 105 high S/N Spitzer/IRS spectra of star-forming galaxies, we show that the equivalent width ratio of the 12.7 and 11.2 micron PAH bands is independent on the optical depth, with small dispersion of ~5% indicative of a nearly constant intrinsic flux ratio R_int = (f_12.7/f_11.2)_int = 0.377 +/- 0.020. Conversely, the observed flux ratio, R_obs = (f_12.7/f_11.2)_obs strongly correlates with the silicate strength (S_sil) confirming that differences in R_obs reflect variation in the optical depth. The relation between R_obs and S_sil reproduces predictions for the Galactic Centre extinction law but disagrees with other laws. We calibrate the total extinction affecting the 11.2 micron PAH from R_obs, which we apply to another sample of 215 galaxies with accurate measurements of the total infrared luminosity (L_IR) to investigate the impact of extinction on L_11.2/L_IR. Correlation between L_11.2/L_IR and R_obs independently on L_IR suggests that increased extinction explains the well known decrease in the average L_11.2/L_IR at high L_IR. The extinction-corrected L_11.2 is proportional to L_IR in the range L_IR/L_sun = 10^9--10^13. These results consolidate L_11.2 as a robust tracer of star formation in galaxies. △ Less

Submitted 1 September, 2020; originally announced September 2020.

Comments: 14 pages, 9 figures, 2 tables. Accepted for publication in MNRAS

Feedback and feeding in the context of galaxy evolution with SPICA: direct characterization of molecular outflows and inflows

Authors: E. González-Alfonso, L. Armus, F. J. Carrera, V. Charmandaris, A. Efstathiou, E. Egami, J. A. Fernández-Ontiveros, J. Fischer, G. L. Granato, C. Gruppioni, E. Hatziminaoglou, M. Imanishi, N. Isobe, H. Kaneda, D. Koziel-Wierzbowska, M. A. Malkan, J. Martin-Pintado, S. Mateos, H. Matsuhara, G. Miniutti, T. Nakagawa, F. Pozzi, F. Rico-Villas, G. Rodighiero, P. Roelfsema , et al. (6 additional authors not shown)

Abstract: A far-infrared observatory such as the {\it SPace Infrared telescope for Cosmology and Astrophysics} ({\it SPICA}), with its unprecedented spectroscopic sensitivity, would unveil the role of feedback in galaxy evolution during the last $\sim10$ Gyr of the Universe ($z=1.5-2$), through the use of far- and mid-infrared molecular and ionic fine structure lines that trace outflowing and infalling gas.… ▽ More A far-infrared observatory such as the {\it SPace Infrared telescope for Cosmology and Astrophysics} ({\it SPICA}), with its unprecedented spectroscopic sensitivity, would unveil the role of feedback in galaxy evolution during the last $\sim10$ Gyr of the Universe ($z=1.5-2$), through the use of far- and mid-infrared molecular and ionic fine structure lines that trace outflowing and infalling gas. Outflowing gas is identified in the far-infrared through P-Cygni line shapes and absorption blueshifted wings in molecular lines with high dipolar moments, and through emission line wings of fine-structure lines of ionized gas. We quantify the detectability of galaxy-scale massive molecular and ionized outflows as a function of redshift in AGN-dominated, starburst-dominated, and main-sequence galaxies, explore the detectability of metal-rich inflows in the local Universe, and describe the most significant synergies with other current and future observatories that will measure feedback in galaxies via complementary tracers at other wavelengths. △ Less

Submitted 6 October, 2017; originally announced October 2017.

Comments: This paper belongs to the SPICA Special Issue on PASA. Accepted for publication in PASA

The active nucleus of the ULIRG IRAS F00183-7111 viewed by NuSTAR

Authors: K. Iwasawa, H. W. W. Spoon, A. Comastri, R. Gilli, G. Lanzuisi, E. Piconcelli, C. Vignali, M. Brusa, S. Puccetti

Abstract: We present an X-ray study of the ultra-luminous infrared galaxy IRAS F00183-7111 (z=0.327), using data obtained from NuSTAR, Chandra X-ray Observatory, Suzaku and XMM-Newton. The Chandra imaging shows that a point-like X-ray source is located at the nucleus of the galaxy at energies above 2 keV. However, the point source resolves into diffuse emission at lower energies, extending to the east, wher… ▽ More We present an X-ray study of the ultra-luminous infrared galaxy IRAS F00183-7111 (z=0.327), using data obtained from NuSTAR, Chandra X-ray Observatory, Suzaku and XMM-Newton. The Chandra imaging shows that a point-like X-ray source is located at the nucleus of the galaxy at energies above 2 keV. However, the point source resolves into diffuse emission at lower energies, extending to the east, where the extranuclear [O III] emission, presumably induced by a galactic-scale outflow, is present. The nuclear source is detected by NuSTAR up to the rest-frame 30 keV. The strong, high-ionization Fe K line, first seen by XMM-Newton, and subsequently by Suzaku and Chandra, is not detected in the NuSTAR data. The line flux appears to have been declining continuously between 2003 and 2016, while the continuum emission remained stable to within 30%. The X-ray continuum below 10 keV is characterised by a hard spectrum caused by cold absorption of nH ~1e23 cm-2, compatible to that of the silicate absorption at 9.7 micron, and a broad absorption feature around 8 keV which we attribute to a high-ionization Fe K absorption edge. The latter is best described by a blueshifted, high-ionization (log xi ~3) absorber. No extra hard component, which would arise from a Compton-thick source, is seen in the NuSTAR data. While a pure reflection scenario (with a totally hidden central source) is viable, direct emission from the central source of L(2-10 keV) = 2e44 erg/s, behind layers of cold and hot absorbing gas may be an alternative explanation. In this case, the relative X-ray quietness (Lx/L_AGN ~6e-3), the high-ionization Fe line, strong outflows inferred from various observations, and other similarities to the well-studied ULIRG/QSO Mrk 231 point that the central source in this ULIRG might be accreting close to the Eddington limit. △ Less

Submitted 6 September, 2017; originally announced September 2017.

Comments: 9 pages, 7 figures, Accepted for publication in A&A

Journal ref: A&A 606, A117 (2017)

Molecular outflows in local ULIRGs: energetics from multi-transition OH analysis

Authors: E. González-Alfonso, J. Fischer, H. W. W. Spoon, K. P. Stewart, M. L. N. Ashby, S. Veilleux, H. A. Smith, E. Sturm, D. Farrah, N. Falstad, M. Meléndez, J. Graciá-Carpio, A. W. Janssen, V. Lebouteiller

Abstract: We report on the energetics of molecular outflows in 14 local Ultraluminous Infrared Galaxies (ULIRGs) that show unambiguous outflow signatures (P-Cygni profiles or high-velocity absorption wings) in the far-infrared lines of OH measured with the Herschel/PACS spectrometer. Detection of both ground-state (at 119 and 79 um) and one or more radiatively-excited (at 65 and 84 um) lines allows us to mo… ▽ More We report on the energetics of molecular outflows in 14 local Ultraluminous Infrared Galaxies (ULIRGs) that show unambiguous outflow signatures (P-Cygni profiles or high-velocity absorption wings) in the far-infrared lines of OH measured with the Herschel/PACS spectrometer. Detection of both ground-state (at 119 and 79 um) and one or more radiatively-excited (at 65 and 84 um) lines allows us to model the nuclear gas (<~300 pc) as well as the more extended components using spherically symmetric radiative transfer models. The highest molecular outflow velocities are found in buried sources, in which slower but massive expansion of the nuclear gas is also observed. With the exception of a few outliers, the outflows have momentum fluxes of (2-5)xL_IR/c and mechanical luminosities of (0.1-0.3)% of L_IR. The moderate momentum boosts in these sources (<~3) suggest that the outflows are mostly momentum-driven by the combined effects of AGN and nuclear starbursts, as a result of radiation pressure, winds, and supernovae remnants. In some sources (~20%), however, powerful (10^{10.5-11} Lsun) AGN feedback and (partially) energy-conserving phases are required, with momentum boosts in the range 3-20. These outflows appear to be stochastic strong-AGN feedback events that occur throughout the merging process. In a few sources, the outflow activity in the innermost regions has subsided in the last ~1 Myr. While OH traces the molecular outflows at sub-kpc scales, comparison of the masses traced by OH with those previously inferred from tracers of more extended outflowing gas suggests that most mass is loaded (with loading factors of Mdot/SFR=1-10) from the central galactic cores (a few x 100 pc). Outflow depletion timescales are <10^8 yr, shorter than the gas consumption timescales by factors of 1.1-15, and are anti-correlated with the AGN luminosity. △ Less

Submitted 11 February, 2017; v1 submitted 24 December, 2016; originally announced December 2016.

Journal ref: The Astrophysical Journal, 2017, 836, 11

The Infrared Database of Extragalactic Observables from Spitzer I: the redshift catalog

Authors: Antonio Hernán-Caballero, Henrik W. W. Spoon, Vianney Lebouteiller, David S. N. Rupke, Donald P. Barry

Abstract: This is the first of a series of papers on the Infrared Database of Extragalactic Observables from Spitzer (IDEOS). In this work we describe the identification of optical counterparts of the infrared sources detected in Spitzer Infrared Spectrograph (IRS) observations, and the acquisition and validation of redshifts. The IDEOS sample includes all the spectra from the Cornell Atlas of Spitzer/IRS S… ▽ More This is the first of a series of papers on the Infrared Database of Extragalactic Observables from Spitzer (IDEOS). In this work we describe the identification of optical counterparts of the infrared sources detected in Spitzer Infrared Spectrograph (IRS) observations, and the acquisition and validation of redshifts. The IDEOS sample includes all the spectra from the Cornell Atlas of Spitzer/IRS Sources (CASSIS) of galaxies beyond the Local Group. Optical counterparts were identified from correlation of the extraction coordinates with the NASA Extragalactic Database (NED). To confirm the optical association and validate NED redshifts, we measure redshifts with unprecedented accuracy on the IRS spectra (σ(dz/(1+z))=0.0011) by using an improved version of the maximum combined pseudo-likelihood method (MCPL). We perform a multi-stage verification of redshifts that considers alternate NED redshifts, the MCPL redshift, and visual inspection of the IRS spectrum. The statistics is as follows: the IDEOS sample contains 3361 galaxies at redshift 0<z<6.42 (mean: 0.48, median: 0.14). We confirm the default NED redshift for 2429 sources and identify 124 with incorrect NED redshifts. We obtain IRS-based redshifts for 568 IDEOS sources without optical spectroscopic redshifts, including 228 with no previous redshift measurements. We provide the entire IDEOS redshift catalog in machine-readable formats. The catalog condenses our compilation and verification effort, and includes our final evaluation on the most likely redshift for each source, its origin, and reliability estimates. △ Less

Submitted 23 November, 2015; originally announced November 2015.

Comments: 11 pages, 6 figures, 1 table. Accepted for publication in MNRAS. Full redshift table in machine-readable format available at http://ideos.astro.cornell.edu/redshifts.html

CASSIS: The Cornell Atlas of Spitzer/Infrared Spectrograph Sources. II. High-resolution observations

Authors: V. Lebouteiller, D. J. Barry, C. Goes, G. C. Sloan, H. W. W. Spoon, D. W. Weedman, J. Bernard-Salas, J. R. Houck

Abstract: The Infrared Spectrograph (IRS) on board the Spitzer Space Telescope observed about 15,000 objects during the cryogenic mission lifetime. Observations provided low-resolution (R~60-127) spectra over ~5-38um and high-resolution (R~600) spectra over ~10-37um. The Cornell Atlas of Spitzer/IRS Sources (CASSIS) was created to provide publishable quality spectra to the community. Low-resolution spectra… ▽ More The Infrared Spectrograph (IRS) on board the Spitzer Space Telescope observed about 15,000 objects during the cryogenic mission lifetime. Observations provided low-resolution (R~60-127) spectra over ~5-38um and high-resolution (R~600) spectra over ~10-37um. The Cornell Atlas of Spitzer/IRS Sources (CASSIS) was created to provide publishable quality spectra to the community. Low-resolution spectra have been available in CASSIS since 2011, and we present here the addition of the high-resolution spectra. The high-resolution observations represent approximately one third of all staring observations performed with the IRS instrument. While low-resolution observations are adapted to faint objects and/or broad spectral features (e.g., dust continuum, molecular bands), high-resolution observations allow more accurate measurements of narrow features (e.g., ionic emission lines) as well as a better sampling of the spectral profile of various features. Given the narrow aperture of the two high-resolution modules, cosmic ray hits and spurious features usually plague the spectra. Our pipeline is designed to minimize these effects through various improvements. A super sampled point-spread function was created in order to enable the optimal extraction in addition to the full aperture extraction. The pipeline selects the best extraction method based on the spatial extent of the object. For unresolved sources, the optimal extraction provides a significant improvement in signal-to-noise ratio over a full aperture extraction. We have developed several techniques for optimal extraction, including a differential method that eliminates low-level rogue pixels (even when no dedicated background observation was performed). The updated CASSIS repository now includes all the spectra ever taken by the IRS, with the exception of mapping observations. △ Less

Submitted 25 June, 2015; originally announced June 2015.

Comments: Accepted on May 11th 2015 in ApJS

Journal ref: Lebouteiller et al., 2015, ApJS 218, 21

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

Authors: Antonio Hernán-Caballero, Almudena Alonso-Herrero, Evanthia Hatziminaoglou, Henrik W. W. Spoon, Cristina Ramos Almeida, Tanio Díaz Santos, Sebastian F. Hönig, Omaira González-Martín, Pilar Esquej

Abstract: 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 linea… ▽ More 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. △ Less

Submitted 20 February, 2015; originally announced February 2015.

Comments: 16 pages, 15 figures, 2 tables, accepted for publication in ApJ

High-lying OH absorption, [C II] deficits, and extreme $L_{\mathrm{FIR}}/M_{\mathrm{H2}}$ ratios in galaxies

Authors: E. González-Alfonso, J. Fischer, E. Sturm, J. Graciá-Carpio, S. Veilleux, M. Meléndez, D. Lutz, A. Poglitsch, S. Aalto, N. Falstad, H. W. W. Spoon, D. Farrah, A. Blasco, C. Henkel, A. Contursi, A. Verma, M. Spaans, H. A. Smith, M. L. N. Ashby, S. Hailey-Dunsheath, S. García-Burillo, J. Martín-Pintado, P. van der Werf, R. Meijerink, R. Genzel

Abstract: Herschel/PACS observations of 29 local (Ultra-)Luminous Infrared Galaxies, including both starburst and AGN-dominated sources as diagnosed in the mid-infrared/optical, show that the equivalent width of the absorbing OH 65 um Pi_{3/2} J=9/2-7/2 line (W_{eq}(OH65)) with lower level energy E_{low}~300 K, is anticorrelated with the [C ii]158 um line to far-infrared luminosity ratio, and correlated wit… ▽ More Herschel/PACS observations of 29 local (Ultra-)Luminous Infrared Galaxies, including both starburst and AGN-dominated sources as diagnosed in the mid-infrared/optical, show that the equivalent width of the absorbing OH 65 um Pi_{3/2} J=9/2-7/2 line (W_{eq}(OH65)) with lower level energy E_{low}~300 K, is anticorrelated with the [C ii]158 um line to far-infrared luminosity ratio, and correlated with the far-infrared luminosity per unit gas mass and with the 60-to-100 um far-infrared color. While all sources are in the active L_{IR}/M_{H2}>50 Lsun/Msun mode as derived from previous CO line studies, the OH65 absorption shows a bimodal distribution with a discontinuity at L_{FIR}/M_{H2}~100 Lsun/Msun. In the most buried sources, OH65 probes material partially responsible for the silicate 9.7 um absorption. Combined with observations of the OH 71 um Pi_{1/2} J=7/2-5/2 doublet (E_{low}~415 K), radiative transfer models characterized by the equivalent dust temperature, Tdust, and the continuum optical depth at 100 um, tau_{100}, indicate that strong [C ii]158 um deficits are associated with far-IR thick (tau_{100}>~0.7, N_{H}>~10^{24} cm^{-2}), warm (Tdust>~60 K) structures where the OH 65 um absorption is produced, most likely in circumnuclear disks/tori/cocoons. With their high L_{FIR}/M_{H2} ratios and columns, the presence of these structures is expected to give rise to strong [C ii] deficits. W_{eq}(OH65) probes the fraction of infrared luminosity arising from these compact/warm environments, which is >~30-50% in sources with high W_{eq}({OH65}). Sources with high W_{eq}({OH65}) have surface densities of both L_{IR} and M_{H2} higher than inferred from the half-light (CO or UV/optical) radius, tracing coherent structures that represent the most buried/active stage of (circum)nuclear starburst-AGN co-evolution. △ Less

Submitted 15 December, 2014; originally announced December 2014.

Comments: 11 pages, 7 figures, accepted for publication in The Astrophysical Journal

Extended Warm Gas in the ULIRG Mrk273: Galactic Outflows and Tidal Debris

Authors: Javier Rodriguez Zaurin, Clive N. Tadhunter, David S. N. Rupke, Sylvain Veilleux, Henrik W. W. Spoon, Marco Chiaberge, Cristina Ramos Almeida, Dan Batcheldor, William B. Sparks

Abstract: We present new HST ACS medium- and narrow-band images and long-slit, optical (4000 - 7200A) spectra obtained using the Isaac Newton Telescope (INT) on La Palma, of the merging system Mrk273. The HST observations sample the [OIII]4959,5007 emission from the galaxy and the nearby continuum. The images show that the morphologies of the extended continuum and the ionised gas emission from the galaxy a… ▽ More We present new HST ACS medium- and narrow-band images and long-slit, optical (4000 - 7200A) spectra obtained using the Isaac Newton Telescope (INT) on La Palma, of the merging system Mrk273. The HST observations sample the [OIII]4959,5007 emission from the galaxy and the nearby continuum. The images show that the morphologies of the extended continuum and the ionised gas emission from the galaxy are decoupled, extending almost perpendicular to each other. In particular, we detect for the first time a spectacular structure of ionised gas in the form of filaments extending ~23 kpc to the east of the nuclear region. The quiescent ionised gas kinematics at these locations suggests that these filaments are tidal debris left over from a secondary merger event that are illuminated by an AGN in the nuclear regions. The images also reveal a complex morphology in the nuclear region of the galaxy for both the continuum and the [OIII] emission. Kinematic disturbance, in the form of broad (FWHM > 500 km s-1) and/or strongly shifted (abs(\DeltaV) >150 km s-1) emission line components, is found at almost all locations within a radius of ~4 kpc to the east and west of the northern nucleus. We fit the profiles of all the emission lines of different ionisation with a kinematic model using up to 3 Gaussian components. From these fits we derive diagnostic line ratios that are used to investigate the ionisation mechanisms at the different locations in the galaxy. We show that, in general, the line ratios are consistent with photoionization by an AGN as the main ionisation mechanism. Finally, the highest surface brightness [OIII] emission is found in a compact region that is coincident with the so-called SE nuclear component. The compactness, kinematics and emission line ratios of this component suggest that it is a separate nucleus with its own AGN. △ Less

Submitted 23 July, 2014; originally announced July 2014.

Comments: Accepted for Publication in A&A

Mid-Infrared Properties of Luminous Infrared Galaxies II: Probing the Dust and Gas Physics of the GOALS Sample

Authors: Sabrina Stierwalt, Lee Armus, Vassilis Charmandaris, Tanio Diaz-Santos, Jason Marshall, Aaron Evans, Sebastian Haan, Justin Howell, Kazushi Iwasawa, Dongchan Kim, Eric J. Murphy, Jeff A. Rich, Henrik W. W. Spoon, Hanae Inami, Andreea Petric, Vivian U

Abstract: 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 indiv… ▽ More 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. △ Less

Submitted 15 June, 2014; originally announced June 2014.

Comments: 21 pages with 20 figures plus 2 tables

The Mrk 231 molecular outflow as seen in OH

Authors: E. González-Alfonso, J. Fischer, J. Graciá-Carpio, N. Falstad, E. Sturm, M. Meléndez, H. W. W. Spoon, A. Verma, R. I. Davies, D. Lutz, S. Aalto, E. Polisensky, A. Poglitsch, S. Veilleux, A. Contursi

Abstract: We report on the Herschel/PACS observations of OH in Mrk 231, with detections in 9 doublets observed within the PACS range, and present radiative transfer models for the outflowing OH. Signatures of outflowing gas are found in up to 6 OH doublets with different excitation requirements. At least two outflowing components are identified, one with OH radiatively excited, and the other with low excita… ▽ More We report on the Herschel/PACS observations of OH in Mrk 231, with detections in 9 doublets observed within the PACS range, and present radiative transfer models for the outflowing OH. Signatures of outflowing gas are found in up to 6 OH doublets with different excitation requirements. At least two outflowing components are identified, one with OH radiatively excited, and the other with low excitation, presumably spatially extended. Particularly prominent, the blue wing of the absorption detected in the in-ladder 2Pi_{3/2} J=9/2-7/2 OH doublet at 65 um, with E_lower=290 K, indicates that the excited outflowing gas is generated in a compact and warm (circum)nuclear region. Because the excited, outflowing OH gas in Mrk 231 is associated with the warm, far-IR continuum source, it is likely more compact (diameter of 200-300 pc) than that probed by CO and HCN. Nevertheless, its mass-outflow rate per unit of solid angle as inferred from OH is similar to that previously derived from CO, >~70x(2.5x10^{-6}/X_{OH}) Msun yr^{-1} sr^{-1}, where X_{OH} is the OH abundance relative to H nuclei. In spherical symmetry, this would correspond to >~850x(2.5x10^{-6}/X_{OH}) Msun yr^{-1}, though significant collimation is inferred from the line profiles. The momentum flux of the excited component attains ~15 L_{AGN}/c, with an OH column density of (1.5-3)x10^{17} cm^-2 and a mechanical luminosity of ~10^{11} Lsun. The detection of very excited OH peaking at central velocities indicates the presence of a nuclear reservoir of gas rich in OH, plausibly the 130-pc scale circumnuclear torus previously detected in OH megamaser emission, that may be feeding the outflow. An exceptional ^{18}OH enhancement, with OH/^{18}OH<~30 at both central and blueshifted velocities, is likely the result of interstellar-medium processing by recent starburst/SNe activity. △ Less

Submitted 31 October, 2013; v1 submitted 11 October, 2013; originally announced October 2013.

Comments: 20 pages, 14 figures, accepted for publication in Astronomy & Astrophysics

Learning the Fundamental MIR Spectral Components of Galaxies with Non-Negative Matrix Factorisation

Authors: P. D. Hurley, S. Oliver, D. Farrah, V. Lebouteiller, H. W. W. Spoon

Abstract: The mid-infrared (MIR) spectra observed with the \textit{Spitzer} Infrared Spectrograph (IRS) provide a valuable dataset for untangling the physical processes and conditions within galaxies. This paper presents the first attempt to blindly learn fundamental spectral components of MIR galaxy spectra, using non-negative matrix factorisation (NMF). NMF is a recently developed multivariate technique… ▽ More The mid-infrared (MIR) spectra observed with the \textit{Spitzer} Infrared Spectrograph (IRS) provide a valuable dataset for untangling the physical processes and conditions within galaxies. This paper presents the first attempt to blindly learn fundamental spectral components of MIR galaxy spectra, using non-negative matrix factorisation (NMF). NMF is a recently developed multivariate technique shown to be successful in blind source separation problems. Unlike the more popular multivariate analysis technique, principal component analysis, NMF imposes the condition that weights and spectral components are non-negative. This more closely resembles the physical process of emission in the mid-infrared, resulting in physically intuitive components. By applying NMF to galaxy spectra in the Cornell Atlas of Spitzer/IRS sources (CASSIS), we find similar components amongst different NMF sets. These similar components include two for AGN emission and one for star formation. [... ABBREVIATED...] We show an NMF set with seven components can reconstruct the general spectral shape of a wide variety of objects, though struggle to fit the varying strength of emission lines. We also show that the seven components can be used to separate out different types of objects. We model this separation with Gaussian Mixtures modelling and use the result to provide a classification tool. We also show the NMF components can be used to separate out the emission from AGN and star formation regions and define a new star formation/AGN diagnostic which is consistent with all mid-infrared diagnostics already in use but has the advantage that it can be applied to mid-infrared spectra with low signal to noise or with limited spectral range. The 7 NMF components and code for classification are made public on arxiv and are available at: \url{https://github.com/pdh21/NMF_software/} △ Less

Submitted 2 October, 2013; originally announced October 2013.

Comments: 23 pages, 19 figures, accepted for publication in MNRAS

Herschel observations and a model for IRAS 08572+3915: a candidate for the most luminous infrared galaxy in the local (z < 0.2) Universe

Authors: A. Efstathiou, C. Pearson, D. Farrah, D. Rigopoulou, J. Gracia-Carpio, A. Verma, H. W. W. Spoon, J. Afonso, J. Bernard-Salas, D. L. Clements, A. Cooray, D. Cormier, M. Etxaluze, J. Fischer, E. Gonzalez-Alfonso, P. Hurley, V. Lebouteiller, S. J. Oliver, M. Rowan-Robinson, E. Sturm

Abstract: We present Herschel photometry and spectroscopy, carried out as part of the Herschel ULIRG survey (HERUS), and a model for the infrared to submillimetre emission of the ultraluminous infrared galaxy IRAS 08572+3915. This source shows one of the deepest known silicate absorption features and no polycyclic aromatic hydrocarbon (PAH) emission. The model suggests that this object is powered by an acti… ▽ More We present Herschel photometry and spectroscopy, carried out as part of the Herschel ULIRG survey (HERUS), and a model for the infrared to submillimetre emission of the ultraluminous infrared galaxy IRAS 08572+3915. This source shows one of the deepest known silicate absorption features and no polycyclic aromatic hydrocarbon (PAH) emission. The model suggests that this object is powered by an active galactic nucleus (AGN) with a fairly smooth torus viewed almost edge-on and a very young starburst. According to our model the AGN contributes about 90% of the total luminosity of 1.1 x 10^13 Lo, which is about a factor of five higher than previous estimates. The large correction of the luminosity is due to the anisotropy of the emission of the best fit torus. Similar corrections may be necessary for other local and high-z analogs. This correction implies that IRAS 08572+3915 at a redshift of 0.05835 may be the nearest hyperluminous infrared galaxy and probably the most luminous infrared galaxy in the local (z < 0.2) Universe. IRAS 08572+3915 shows a low ratio of [CII] to IR luminosity (log L_[CII]/L_{IR} < -3.8) and a [OI]63um to [CII]158um line ratio of about 1 that supports the model presented in this letter. △ Less

Submitted 26 September, 2013; originally announced September 2013.

Comments: MNRAS Letters accepted

Diagnostics of AGN-driven Molecular Outflows in ULIRGs from Herschel-PACS Observations of OH at 119um

Authors: H. W. W. Spoon, D. Farrah, V. Lebouteiller, E. Gonzalez-Alfonso, J. Bernard-Salas, T. Urrutia, D. Rigopoulou, M. S. Westmoquette, H. A. Smith, J. Afonso, C. Pearson, D. Cormier, A. Efstathiou, C. Borys, A. Verma, M. Etxaluze, D. L. Clements

Abstract: We report on our observations of the 79 and 119um doublet transitions of OH for 24 local (z<0.262) ULIRGs observed with Herschel-PACS as part of the Herschel ULIRG Survey (HERUS). Some OH119 profiles display a clear P-Cygni shape and therefore imply outflowing OH gas, other profiles are predominantly in absorption or are completely in emission. We find that the relative strength of the OH emission… ▽ More We report on our observations of the 79 and 119um doublet transitions of OH for 24 local (z<0.262) ULIRGs observed with Herschel-PACS as part of the Herschel ULIRG Survey (HERUS). Some OH119 profiles display a clear P-Cygni shape and therefore imply outflowing OH gas, other profiles are predominantly in absorption or are completely in emission. We find that the relative strength of the OH emission component decreases as the silicate absorption increases. This locates the OH outflows inside the obscured nuclei. The maximum outflow velocities for our sources range from less than 100 to 2000 km/s, with 15/24 (10/24) sources showing OH absorption at velocities exceeding 700 km/s (1000 km/s). Three sources show maximum OH outflow velocities exceeding that of Mrk231. Since outflow velocities above 500-700 km/s are thought to require an active galactic nucleus (AGN) to drive them, about 2/3 of our ULIRG sample may host AGN-driven molecular outflows. This finding is supported by the correlation we find between the maximum OH outflow velocity and the IR-derived bolometric AGN luminosity. No such correlation is found with the IR-derived star formation rate. The highest outflow velocities are found among sources which are still deeply embedded. We speculate that the molecular outflows in these sources may be in an early phase of disrupting the nuclear dust veil before these sources evolve into less obscured AGN. Four of our sources show high-velocity wings in their [C II] fine-structure line profiles implying neutral gas outflow masses of at least 2-4.5 x 10^8 Msun. △ Less

Submitted 23 July, 2013; originally announced July 2013.

Comments: 23 pages, 12 figures, 4 tables

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

Authors: S. Stierwalt, L. Armus, J. A. Surace, H. Inami, A. O. Petric, T. Diaz-Santos, S. Haan, V. Charmandaris, J. Howell, D. C. Kim, J. Marshall, J. M. Mazzarella, H. W. W. Spoon, S. Veilleux, A. Evans, D. B. Sanders, P. Appleton, G. Bothun, C. R. Bridge, B. Chan, D. Frayer, K. Iwasawa, L. J. Kewley, S. Lord, B. F. Madore , et al. (8 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 18 February, 2013; originally announced February 2013.

Comments: Accepted for publication in ApJS

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

Authors: J. M. Mazzarella, K. Iwasawa, T. Vavilkin, L. Armus, D. -C. Kim, G. Bothun, A. S. Evans, H. W. W. Spoon, S. Haan, J. H. Howell, S. Lord, J. A. Marshall, C. M. Ishida, C. K. Xu, A. Petric, D. B. Sanders, J. A. Surace, P. Appleton, B. H. P. Chan, D. T. Frayer, H. Inami, E. Ye. Khachikian, B. F. Madore, G. C. Privon, E. Sturm , et al. (2 additional authors not shown)

Abstract: 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 likel… ▽ More 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. △ Less

Submitted 15 August, 2012; originally announced August 2012.

Comments: Accepted for publication in the Astronomical Journal, 40 pages, 31 figures, 15 tables

Multi-wavelength GOALS Observations of Star Formation and Active Galactic Nucleus Activity in the Luminous Infrared Galaxy IC 883

Authors: F. Modica, T. Vavilkin, A. S. Evans, D. C. Kim, J. M. Mazzarella, K. Iwasawa, A. Petric, J. H. Howell, J. A. Surace, L. Armus, H. W. W. Spoon, D. B. Sanders, J. E. Barnes

Abstract: 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… ▽ More 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μ$m) are both indicative of a relatively weak AGN. (abridged) △ Less

Submitted 14 November, 2011; originally announced November 2011.

Comments: 16 pages, 12 Figures. Accepted for publication in AJ

Herschel/PACS spectroscopy of NGC 4418 and Arp 220: H2O, H2^{18}O, OH, ^{18}OH, O I, HCN and NH3

Authors: E. González-Alfonso, J. Fischer, J. Graciá-Carpio, E. Sturm, S. Hailey-Dunsheath, D. Lutz, A. Poglitsch, A. Contursi, H. Feuchtgruber, S. Veilleux, H. W. W. Spoon, A. Verma, N. Christopher, R. Davies, A. Sternberg, R. Genzel, L. Tacconi

Abstract: Herschel/PACS spectroscopy of the luminous infrared galaxies NGC4418 and Arp220 reveals high excitation in H2O, OH, HCN, and NH3. In NGC4418, absorption lines were detected with E_low>800 K (H2O), 600 K (OH), 1075 K (HCN), and 600 K (NH3), while in Arp220 the excitation is somewhat lower. While outflow signatures in moderate excitation lines are seen in Arp220 as reported in previous studies, in N… ▽ More Herschel/PACS spectroscopy of the luminous infrared galaxies NGC4418 and Arp220 reveals high excitation in H2O, OH, HCN, and NH3. In NGC4418, absorption lines were detected with E_low>800 K (H2O), 600 K (OH), 1075 K (HCN), and 600 K (NH3), while in Arp220 the excitation is somewhat lower. While outflow signatures in moderate excitation lines are seen in Arp220 as reported in previous studies, in NGC4418 the lines tracing its outer regions are redshifted relative to the nucleus, suggesting an inflow with Mdot<~12 Msun yr^{-1}. Both galaxies have warm (Tdust>~100 K) nuclear continuum components, together with a more extended component that is much more prominent and massive in Arp220. A chemical dichotomy is found in both sources: on the one hand, the nuclear regions have high H2O abundances, ~10^{-5}, and high HCN/H2O and HCN/NH3 column density ratios of 0.1-0.4 and 2-5, respectively, indicating a chemistry typical of evolved hot cores where grain mantle evaporation has occurred. On the other hand, the high OH abundance, with OH/H2O ratios of ~0.5, indicates the effects of X-rays and/or cosmic rays. The nuclear media have surface brightnesses >~10^{13} Lsun/kpc^2 and are estimated to be thick (N_H>~10^{25} cm^{-2}). While NGC4418 shows weak absorption in H2^{18}O and ^{18}OH, with a ^{16}O-to-^{18}O ratio of >~250-500, the strong absorption of the rare isotopologues in Arp220 indicates ^{16}O-to-^{18}O of 70-130. Further away from the nuclear regions, the H2O abundance decreases to <~10^{-7} and the OH/H2O ratio is 2.5-10. Despite the different scales of NGC4418, Arp220, and Mrk231, preliminary evidence is found for an evolutionary sequence from infall, hot-core like chemistry, and solar oxygen isotope ratio to high velocity outflow, disruption of the hot core chemistry and cumulative high mass stellar processing of 18O. △ Less

Submitted 12 February, 2012; v1 submitted 6 September, 2011; originally announced September 2011.

Comments: 31 pages, 20 figures. Accepted for publication in Astronomy and Astrophysics

CASSIS: The Cornell Atlas of Spitzer/IRS Sources

Authors: V. Lebouteiller, D. J. Barry, H. W. W. Spoon, J. Bernard-Salas, G. C. Sloan, J. R. Houck, D. W. Weedman

Abstract: We present the spectral atlas of sources observed in low resolution with the Infrared Spectrograph on board the Spitzer Space Telescope. More than 11,000 distinct sources were extracted using a dedicated algorithm based on the SMART software with an optimal extraction (AdOpt package). These correspond to all 13,000 low resolution observations of fixed objects (both single source and cluster observ… ▽ More We present the spectral atlas of sources observed in low resolution with the Infrared Spectrograph on board the Spitzer Space Telescope. More than 11,000 distinct sources were extracted using a dedicated algorithm based on the SMART software with an optimal extraction (AdOpt package). These correspond to all 13,000 low resolution observations of fixed objects (both single source and cluster observations). The pipeline includes image cleaning, individual exposure combination, and background subtraction. A particular attention is given to bad pixel and outlier rejection at the image and spectra levels. Most sources are spatially unresolved so that optimal extraction reaches the highest possible signal-to-noise ratio. For all sources, an alternative extraction is also provided that accounts for all of the source flux within the aperture. CASSIS provides publishable quality spectra through an online database together with several important diagnostics, such as the source spatial extent and a quantitative measure of detection level. Ancillary data such as available spectroscopic redshifts are also provided. The database interface will eventually provide various ways to interact with the spectra, such as on-the-fly measurements of spectral features or comparisons among spectra. △ Less

Submitted 16 August, 2011; originally announced August 2011.

Comments: Accepted for publication in ApJ Supplement Series

The deeply obscured AGN of NGC4945 I. Spitzer-IRS maps of [Ne V], [Ne II], H2 0-0 S(1), S(2), and other tracers

Authors: J. P. Pérez-Beaupuits, H. W. W. Spoon, M. Spaans, J. D. Smith

Abstract: The nearly edge-on galaxy NGC4945 is one of the closest galaxies where an AGN and starburst coexist, and is one of the brightest sources at 100 keV. Near and mid-infrared spectroscopy have shown very strong obscuration of its central region, rivaled only in strength by some of the most deeply obscured ULIRGs. We aim to determine the spatial distribution of ISM features in the central 426x426 pc^2… ▽ More The nearly edge-on galaxy NGC4945 is one of the closest galaxies where an AGN and starburst coexist, and is one of the brightest sources at 100 keV. Near and mid-infrared spectroscopy have shown very strong obscuration of its central region, rivaled only in strength by some of the most deeply obscured ULIRGs. We aim to determine the spatial distribution of ISM features in the central 426x426 pc^2 of NGC4945. We map the central region of NGC4945 in three Spitzer-IRS modules (SH, SL and LL). We produce maps of the flux distribution of the starburst tracers [Ne II], [Ne III], [S III] and [S IV] at 12.81, 15.56, 18.71 and 10.51 mum, respectively, and a map of the AGN narrow-line region tracer [Ne V] at 14.32 mum. We also mapped the S(1), S(2) and S(3) pure rotational lines of H2, which trace the distribution of warm molecular hydrogen. We obtained an extinction map (A_V) based on the apparent strength of the 9.7 mum silicate absorption feature. Our A_V map traces the contours of the starburst ring but the highest extinction (A_V(9.85 mum)~60) is found at the nucleus. Within the uncertainty of the astrometry all emission lines are found to peak on the nucleus, except for the warm molecular hydrogen emission which shows a maximum 60-100 pc NW of the nucleus. We favour a scenario in which the lower H2 0-0 S(1) and S(2) rotational lines originate mainly from an unobscured extra-nuclear component associated with the super-wind cone observed in the HST NICMOS map of the H2 1-0 S(1) vibrational line. For the [Ne V] emission we infer an attenuation of a factor 12-160 (A_V=55-112) based on a comparison of the ratio of our [Ne V] flux and the absorption-corrected 14-195 keV Swift-BAT flux to the average [Ne V]/BAT ratio for Seyfert 1 nuclei. The high attenuation indicates that [Ne V] and [O IV] cannot be used as extinction-free tracers of AGN power in galaxies with deeply buried nuclei. △ Less

Submitted 1 July, 2011; originally announced July 2011.

Comments: 19 pages, 16 figures, accepted for publication in A&A

Mid-infrared properties of OH megamaser host galaxies. II: Analysis and modeling of the maser environment

Authors: Kyle W. Willett, Jeremy Darling, Henrik W. W. Spoon, Vassilis Charmandaris, Lee Armus

Abstract: We present a comparison of Spitzer IRS data for 51 OH megamaser (OHM) hosts and 15 non-masing ULIRGs. 10-25% of OHMs show evidence for the presence of an AGN, significantly lower than the estimated AGN fraction from previous optical and radio studies. Non-masing ULIRGs have a higher AGN fraction (50-95%) than OHMs, although some galaxies in both samples show evidence of co-existing starbursts and… ▽ More We present a comparison of Spitzer IRS data for 51 OH megamaser (OHM) hosts and 15 non-masing ULIRGs. 10-25% of OHMs show evidence for the presence of an AGN, significantly lower than the estimated AGN fraction from previous optical and radio studies. Non-masing ULIRGs have a higher AGN fraction (50-95%) than OHMs, although some galaxies in both samples show evidence of co-existing starbursts and AGN. Radiative transfer models of the dust environment reveal that non-masing galaxies tend to have clumpy dust geometries commonly associated with AGN, while OHMs have deeper absorption consistent with a smooth, thick dust shell. Statistical analyses show that the major differences between masing and non-masing ULIRGs in the mid-IR relate to the optical depth and dust temperature, which we measure using the 9.7 um silicate depth and 30-20 um spectral slope from the IRS data. Dust temperatures of 40-80 K derived from the IRS data are consistent with predictions of OH pumping models and with a minimum T_dust required for maser production. The best-fit dust opacities (τ_V ~ 100 - 400), however, are nearly an order of magnitude larger than those predicted for OH inversion, and suggest that modifications to the model may be required. These diagnostics offer the first detailed test of an OHM pumping model based only on the properties of its host galaxy and provide important restrictions on the physical conditions relevant to OHM production. △ Less

Submitted 25 January, 2011; originally announced January 2011.

Comments: 15 pages, 12 figures; accepted to ApJ

Mid-infrared properties of OH megamaser host galaxies. I: Spitzer IRS low- and high-resolution spectroscopy

Authors: Kyle W. Willett, Jeremy Darling, Henrik W. W. Spoon, Vassilis Charmandaris, Lee Armus

Abstract: We present mid-infrared spectra and photometry from the Infrared Spectrograph on the Spitzer Space Telescope for 51 OH megamasers (OHMs), along with 15 galaxies confirmed to have no megamaser emission above L_OH = 10^2.3 L_sun. The majority of galaxies display moderate-to-deep 9.7 um amorphous silicate absorption, with OHM galaxies showing stronger average absorption and steeper 20-30 um continuum… ▽ More We present mid-infrared spectra and photometry from the Infrared Spectrograph on the Spitzer Space Telescope for 51 OH megamasers (OHMs), along with 15 galaxies confirmed to have no megamaser emission above L_OH = 10^2.3 L_sun. The majority of galaxies display moderate-to-deep 9.7 um amorphous silicate absorption, with OHM galaxies showing stronger average absorption and steeper 20-30 um continuum emission than non-masing galaxies. Emission from multiple polycyclic aromatic hydrocarbons (PAHs), especially at 6.2, 7.7, and 11.3 um, is detected in almost all systems. Fine-structure atomic emission (including [Ne II], [Ne III], [S III], and [S IV]) and multiple H2 rotational transitions are observed in more than 90% of the sample. A subset of galaxies show emission from rarer atomic lines, such as [Ne V], [O IV], and [Fe II]. 50% of the OHMs show absorption from water ice and hydrogenated amorphous carbon grains (HACs), while absorption features from CO2, HCN, C2H2, and crystalline silicates are also seen in several OHMs. Column densities of OH derived from 34.6 um OH absorption are similar to those derived from 1667 MHz OH absorption in non-masing galaxies, indicating that the abundance of masing molecules is similar for both samples. This data paper presents full mid-infrared spectra for each galaxy, along with measurements of line fluxes and equivalent widths, absorption feature depths, and spectral indices. △ Less

Submitted 25 January, 2011; originally announced January 2011.

Comments: 28 pages, 10 figures; accepted to ApJS. Ancillary data includes full IRS spectra of the complete sample

Molecules as tracers of galaxy evolution: an EMIR survey. I. Presentation of the data and first results

Authors: F. Costagliola, S. Aalto, M. I. Rodriguez, S. Muller, H. W. W. Spoon, S. Martín, M. A. Peréz-Torres, A. Alberdi, J. E. Lindberg, F. Batejat, E. Jütte, P. van der Werf, F. Lahuis

Abstract: We investigate the molecular gas properties of a sample of 23 galaxies in order to find and test chemical signatures of galaxy evolution and to compare them to IR evolutionary tracers. Observation at 3 mm wavelengths were obtained with the EMIR broadband receiver, mounted on the IRAM 30 m telescope on Pico Veleta, Spain. We compare the emission of the main molecular species with existing models of… ▽ More We investigate the molecular gas properties of a sample of 23 galaxies in order to find and test chemical signatures of galaxy evolution and to compare them to IR evolutionary tracers. Observation at 3 mm wavelengths were obtained with the EMIR broadband receiver, mounted on the IRAM 30 m telescope on Pico Veleta, Spain. We compare the emission of the main molecular species with existing models of chemical evolution by means of line intensity ratios diagrams and principal component analysis. We detect molecular emission in 19 galaxies in two 8 GHz-wide bands centred at 88 and 112 GHz. The main detected transitions are the J=1-0 lines of CO, 13CO, HCN, HNC, HCO+, CN, and C2H. We also detect HC3N J=10-9 in the galaxies IRAS 17208, IC 860, NGC 4418, NGC 7771, and NGC 1068. The only HC3N detections are in objects with HCO+/HCN<1 and warm IRAS colours. Galaxies with the highest HC3N/HCN ratios have warm IRAS colours (60/100 μm>0.8). The brightest HC3N emission is found in IC 860, where we also detect the molecule in its vibrationally excited state.We find low HNC/HCN line ratios (<0.5), that cannot be explained by existing PDR or XDR chemical models. Bright HC3N emission in HCO+-faint objects may imply that these are not dominated by X-ray chemistry. Thus the HCN/HCO+ line ratio is not, by itself, a reliable tracer of XDRs. Bright HC3N and faint HCO+ could be signatures of embedded starformation, instead of AGN activity. △ Less

Submitted 11 January, 2011; originally announced January 2011.

Ultra-deep Spitzer Mid-Infrared Spectroscopy of LIRGs and ULIRGs at z ~ 1-2

Authors: D. Fadda, L. Yan, G. Lagache, A. Sajina, D. Lutz, S. Wuyts, D. T. Frayer, D. Marcillac, E. Le Floc'h, K. Caputi, H. W. W. Spoon, S. Veilleux, A. Blain, G. Helou

Abstract: We present ultra-deep mid-IR spectra of 48 infrared-luminous galaxies in the GOODS-South field obtained with the InfraRed Spectrograph (IRS) on the Spitzer Space Telescope. These galaxies are selected among faint infrared sources (0.14 - 0.5 mJy at 24 um) in two redshift bins (0.76-1.05 and 1.75-2.4) to sample the major contributors to the cosmic infrared background at the most active epochs. We e… ▽ More We present ultra-deep mid-IR spectra of 48 infrared-luminous galaxies in the GOODS-South field obtained with the InfraRed Spectrograph (IRS) on the Spitzer Space Telescope. These galaxies are selected among faint infrared sources (0.14 - 0.5 mJy at 24 um) in two redshift bins (0.76-1.05 and 1.75-2.4) to sample the major contributors to the cosmic infrared background at the most active epochs. We estimate redshifts for 92% of the sample using PAH and Si absorption features. Only few of these galaxies (5% at z~1 and 12% at z~2) have their total infrared luminosity dominated by emission from AGN. The averaged mid-IR spectra of the z~1 LIRGs and of the z~2 ULIRGs are very similar to the averaged spectrum of local starbursts and HII-like ULIRGs, respectively. We find that 6.2um PAH equivalent widths reach a plateau of ~1 um for L(24 mu) < 1E11 L(sun). At higher luminosities, EW (6.2 mu) anti-correlates with L(24 um). Intriguingly, high-z ULIRGs and SMG lie above the local EW (6.2 um) - L(24 um) relationship suggesting that, at a given luminosity, high-z ULIRGs have AGN contributions to their dust emission lower than those of local counterparts. A quantitative analysis of their morphology shows that most of the luminous IR galaxies have morphologies similar to those of IR-quiet galaxies at the same redshift. All z~2 ULIRGs of our sample are IR-excess BzK galaxies and most of them have L(FIR)/L(1600A) ratios higher than those of starburst galaxies at a given UV slope. The ``IR excess'' (Daddi et al. 2007) is mostly due to strong 7.7 um PAH emission and under-estimation of UV dust extinction. On the basis of the AGN-powered L (6 um) continuum measured directly from the mid-IR spectra, we estimate an average intrinsic X-ray AGN luminosity of L(2-10 keV) = (0.1 +/- 0.6) 1E43 erg/s, a value substantially lower than the prediction by Daddi et al. (2007). △ Less

Submitted 14 June, 2010; originally announced June 2010.

Comments: 45 pages, 36 figures, 2 tables, accepted for publication on ApJ

The extraordinary mid-infrared spectral properties of FeLoBAL Quasars

Authors: D. Farrah, T. Urrutia, M. Lacy, V. Lebouteiller, H. W. W. Spoon, J. Bernard-Salas, N. Connolly, J. Afonso, B. Connolly, J. Houck

Abstract: We present mid-infrared spectra of six FeLoBAL QSOs at 1<z<1.8, taken with the Spitzer space telescope. The spectra span a range of shapes, from hot dust dominated AGN with silicate emission at 9.7 microns, to moderately obscured starbursts with strong Polycyclic Aromatic Hydrocarbon (PAH) emission. The spectrum of one object, SDSS 1214-0001, shows the most prominent PAHs yet seen in any QSO at an… ▽ More We present mid-infrared spectra of six FeLoBAL QSOs at 1<z<1.8, taken with the Spitzer space telescope. The spectra span a range of shapes, from hot dust dominated AGN with silicate emission at 9.7 microns, to moderately obscured starbursts with strong Polycyclic Aromatic Hydrocarbon (PAH) emission. The spectrum of one object, SDSS 1214-0001, shows the most prominent PAHs yet seen in any QSO at any redshift, implying that the starburst dominates the mid-IR emission with an associated star formation rate of order 2700 solar masses per year. With the caveats that our sample is small and not robustly selected, we combine our mid-IR spectral diagnostics with previous observations to propose that FeLoBAL QSOs are at least largely comprised of systems in which (a) a merger driven starburst is ending, (b) a luminous AGN is in the last stages of burning through its surrounding dust, and (c) which we may be viewing over a restricted line of sight range. △ Less

Submitted 19 May, 2010; originally announced May 2010.

Comments: ApJ, accepted

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

Authors: Hanae Inami, L. Armus, J. A. Surace, J. M. Mazzarella, A. S. Evans, D. B. Sanders, J. H. Howell, A. Petric, T. Vavilkin, K. Iwasawa, S. Haan, E. J. Murphy, S. Stierwalt, P. N. Appleton, J. E. Barnes, G. Bothun, C. R. Bridge, B. Chan, V. Charmandaris, D. T. Frayer, L. J. Kewley, D. C. Kim, S. Lord, B. F. Madore, J. A. Marshall , et al. (9 additional authors not shown)

Abstract: 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_sun) = 11.94, it is classified as a Luminous Infrared Galaxy (LIRG), and was observed as part of the Great Observatories… ▽ More 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_sun) = 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 micron and 6.2 micron PAH 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_sun/yr and > 45 M_sun/yr, 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 magnitudes of visual extinction. The masses of these clusters span a range between 10^6-10^8 M_sun. This starburst source is reminiscent of the extra-nuclear 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 enitre system. △ Less

Submitted 20 April, 2010; originally announced April 2010.

Comments: 46 pages, 10 figures, accepted for publication in AJ

Probing the Excitation of Extreme Starbursts: High Resolution Mid-IR Spectroscopy of Blue Compact Dwarfs

Authors: Lei Hao, Yanling Wu, V. Charmandaris, H. W. W. Spoon, J. Bernard-Salas, D. Devost, V. Lebouteiller, J. R. Houck

Abstract: We present an analysis of the mid-infrared emission lines for a sample of 12 low metallicity Blue Compact Dwarf (BCD) galaxies based on high resolution observations obtained with Infrared Spectrograph on board the {\rm Spitzer} Space Telescope. We compare our sample with a local sample of typical starburst galaxies and active galactic nuclei (AGNs), to study the ionization field of starbursts ov… ▽ More We present an analysis of the mid-infrared emission lines for a sample of 12 low metallicity Blue Compact Dwarf (BCD) galaxies based on high resolution observations obtained with Infrared Spectrograph on board the {\rm Spitzer} Space Telescope. We compare our sample with a local sample of typical starburst galaxies and active galactic nuclei (AGNs), to study the ionization field of starbursts over a broad range of physical parameters and examine its difference from the one produced by AGN. The high-ionization line [OIV]25.89$μ$m is detected in most of the BCDs, starbursts, and AGNs in our sample. We propose a diagnostic diagram of the line ratios [OIV]25.89$μ$m/[SIII]33.48$μ$m as a function of [NeIII]15.56$μ$m/[NeII]12.81$μ$m which can be useful in identifying the principal excitation mechanism in a galaxy. Galaxies in this diagram split naturally into two branches. Classic AGNs as well as starburst galaxies with an AGN component populate the upper branch, with stronger AGNs displaying higher [NeIII]/[NeII] ratios. BCDs and pure starbursts are located in the lower branch. We find that overall the placement of galaxies on this diagram correlates well with their corresponding locations in the log([NII]/H$α$) vs. log([OIII]/H$β$) diagnostic diagram, which has been widely used in the optical. The two diagrams provide consistent classifications of the excitation mechanism in a galaxy. On the other hand, the diagram of [NeIII]15.56$μ$m/[NeII]12.81$μ$m vs. [SIV]10.51$μ$m/[SIII]18.71$μ$m is not as efficient in separating AGNs from BCDs and pure starbursts. (abridged) △ Less

Submitted 27 August, 2009; originally announced August 2009.

Comments: 15 pages, including 1 table and 14 figures. Accepted by ApJ

A Spitzer high resolution mid-infrared spectral atlas of starburst galaxies

Authors: J. Bernard-Salas, H. W. W. Spoon, V. Charmandaris, V. Lebouteiller, D. Farrah, D. Devost, B. R. Brandl, Yanling Wu, L. Armus, L. Hao, G. C. Sloan, D. Weedman, J. R. Houck

Abstract: We present an atlas of Spitzer/IRS high resolution (R~600) 10-37um spectra for 24 well known starburst galaxies. The spectra are dominated by fine-structure lines, molecular hydrogen lines, and emission bands of polycyclic aromatic hydrocarbons. Six out of the eight objects with a known AGN component show emission of the high excitation [NeV] line. This line is also seen in one other object (NGC… ▽ More We present an atlas of Spitzer/IRS high resolution (R~600) 10-37um spectra for 24 well known starburst galaxies. The spectra are dominated by fine-structure lines, molecular hydrogen lines, and emission bands of polycyclic aromatic hydrocarbons. Six out of the eight objects with a known AGN component show emission of the high excitation [NeV] line. This line is also seen in one other object (NGC4194) with, a priori, no known AGN component. In addition to strong polycyclic aromatic hydrocarbon emission features in this wavelength range (11.3, 12.7, 16.4um), the spectra reveal other weak hydrocarbon features at 10.6, 13.5, 14.2um, and a previously unreported emission feature at 10.75um. An unidentified absorption feature at 13.7um is detected in many of the starbursts. We use the fine-structure lines to derive the abundance of neon and sulfur for 14 objects where the HI 7-6 line is detected. We further use the molecular hydrogen lines to sample the properties of the warm molecular gas. Several basic diagrams characterizing the properties of the sample are also shown. We have combined the spectra of all the pure starburst objects to create a high S/N template, which is available to the community. △ Less

Submitted 19 August, 2009; originally announced August 2009.

Comments: 25 pages (emulate apj), 6 tables, 14 figures, Accepted for publication in ApJS

Discovery of strongly blue shifted mid-infrared [NeIII] and [NeV] emission in ULIRGs

Authors: H. W. W. Spoon, J. Holt

Abstract: We report the discovery of blue shifted (delta(V) > 200 km/s) mid-infrared [NeIII] and/or [NeV] emission in 25 out of 82 ULIRGs (30% of our sample). The incidence of blue shifted [NeV] emission is even higher (59%) among the sources with a [NeV] detection -- the tell-tale signature of an active galactic nucleus (AGN). Sixteen ULIRGs in our sample, eleven of which are optically classified as AGN,… ▽ More We report the discovery of blue shifted (delta(V) > 200 km/s) mid-infrared [NeIII] and/or [NeV] emission in 25 out of 82 ULIRGs (30% of our sample). The incidence of blue shifted [NeV] emission is even higher (59%) among the sources with a [NeV] detection -- the tell-tale signature of an active galactic nucleus (AGN). Sixteen ULIRGs in our sample, eleven of which are optically classified as AGN, have [NeIII] blue shifts above 200 km/s. A comparison of the line profiles of their 12.81um [NeII], 15.56um [NeIII] and 14.32um [NeV] lines reveals the ionization of the blue shifted gas to increase with blue shift, implying decelerating outflows in a stratified medium, photo-ionized by the AGN. The strong correlation of the line width of the [NeIII] line with the radio luminosity indicates that interaction of expanding radio jets with the dense ISM surrounding the AGN may explain the observed neon line kinematics for the strongest radio sources in this sample. △ Less

Submitted 24 July, 2009; originally announced July 2009.

Comments: Accepted for publication by ApJ Letters. 15 pages, 4 figures

GOALS: The Great Observatories All-Sky LIRG Survey

Authors: L. Armus, J. M. Mazzarella, A. S. Evans, J. A. Surace, D. B. Sanders, K. Iwasawa, D. T. Frayer, J. H. Howell, B. Chan, A. O. Petric, T. Vavilkin, D. C. Kim, S. Haan, H. Inami, E. J. Murphy, P. N. Appleton, J. E. Barnes, G. Bothun, C. R. Bridge, V. Charmandaris, J. B. Jensen, L. J. Kewley, S. Lord, B. F. Madore, J. A. Marshall , et al. (9 additional authors not shown)

Abstract: The Great Observatories All-sky LIRG Survey (GOALS) combines data from NASA's Spitzer, Chandra, Hubble and GALEX observatories, together with ground-based data into a comprehensive imaging and spectroscopic survey of over 200 low redshift Luminous Infrared Galaxies (LIRGs). The LIRGs are a complete subset of the IRAS Revised Bright Galaxy Sample (RBGS). The LIRGs targeted in GOALS span the full… ▽ More The Great Observatories All-sky LIRG Survey (GOALS) combines data from NASA's Spitzer, Chandra, Hubble and GALEX observatories, together with ground-based data into a comprehensive imaging and spectroscopic survey of over 200 low redshift Luminous Infrared Galaxies (LIRGs). The LIRGs are a complete subset of the IRAS Revised Bright Galaxy Sample (RBGS). The LIRGs targeted in GOALS span the full range of nuclear spectral types defined via traditional optical line-ratio diagrams as well as interaction stages. They provide an unbiased picture of the processes responsible for enhanced infrared emission in galaxies in the local Universe. As an example of the analytic power of the multi-wavelength GOALS dataset, we present data for the interacting system VV 340 (IRAS F14547+2449). Between 80-95% of the total far-infrared emission (or about 5E11 solar luminosities) originates in VV 340 North. While the IRAC colors of VV 340 North and South are consistent with star-forming galaxies, both the Spitzer IRS and Chandra ACIS data indicate the presence of a buried AGN in VV 340 North. The GALEX far and near-UV fluxes imply a extremely large infrared "excess" (IRX) for the system (IR/FUV = 81) which is well above the correlation seen in starburst galaxies. Most of this excess is driven by VV 340 N, which alone has an IR excess of nearly 400. The VV 340 system seems to be comprised of two very different galaxies - an infrared luminous edge-on galaxy (VV 340 North) that dominates the long-wavelength emission from the system and which hosts a buried AGN, and a face-on starburst (VV 340 South) that dominates the short-wavelength emission. △ Less

Submitted 3 June, 2009; v1 submitted 28 April, 2009; originally announced April 2009.

Comments: 17 pages, 2 tables, 7 postscript figures. Accepted for publication in PASP. Updated manuscript includes complete source table (Table 1)

High-velocity neon line emission from the ULIRG IRAS F00183-7111: revealing the optically obscured base of a nuclear outflow

Authors: H. W. W. Spoon, L. Armus, J. A. Marshall, J. Bernard-Salas, D. Farrah, V. Charmandaris, B. R. Kent

Abstract: We report the first mid-IR detection of highly disturbed ionized gas in the ultraluminous infrared galaxy IRAS F00183-7111. The gas, traced by the 12.81um [NeII] and 15.56um [NeIII] lines, spans a velocity range of-3500 to 3000 km/s with respect to systemic velocity. Optical and near-IR spectroscopic studies show no evidence for similarly high velocity gas components in forbidden lines at shorte… ▽ More We report the first mid-IR detection of highly disturbed ionized gas in the ultraluminous infrared galaxy IRAS F00183-7111. The gas, traced by the 12.81um [NeII] and 15.56um [NeIII] lines, spans a velocity range of-3500 to 3000 km/s with respect to systemic velocity. Optical and near-IR spectroscopic studies show no evidence for similarly high velocity gas components in forbidden lines at shorter wavelengths. We interpret this as the result of strong extinction (Av=10-50) on the high-velocity gas, which identifies the base of the outflow traced in 5007A [OIII] as a plausible origin. Unusual excitation conditions are implied by a comparison of the mid-infrared low-excitation neon line emission and the PAH emission for a sample of 56 ULIRGs. For IRAS F00183, the neon/PAH ratio is 8 times higher than the average ratio. Similar mid-infrared kinematic and excitation characteristics are found for only 2 other ULIRGs in our sample: IRAS 12127-1412NE and IRAS 13451+1232. Both sources have an elevated neon/PAH ratio and exhibit pronounced blue wings in their 15.56um [NeIII] line profiles. IRAS 13451 even shows a strongly blue shifted and broad 14.32um [NeV] line. While for IRAS 13451 the observed [NeIII]/[NeII] and [NeV]/[NeII] line ratios indicate exposure of the blue shifted gas to direct radiation from the AGN, for IRAS F00183 and 12127 the neon line ratios are consistent with an origin in fast shocks in a high-density environment, and with an evolutionary scenario in which strongly blue shifted [Ne II] and [Ne III] emission trace the (partial) disruption of the obscuring medium around a buried AGN. The detection of strongly blue shifted [Ne V] emission in IRAS 13451 would then indicate this process to be much further advanced in this ULIRG than in IRAS F00183 and 12127, where this line is undetected. △ Less

Submitted 10 November, 2008; originally announced November 2008.

Comments: 13 pages, 10 figures, accepted for publication in the Astrophysical Journal

Journal ref: Astrophys.J.693:1223-1235,2009

Silicates in Ultra-Luminous Infrared Galaxies

Authors: M. M. Sirocky, N. A. Levenson, M. Elitzur, H. W. W. Spoon, L. Armus

Abstract: We analyze the mid-infrared (MIR) spectra of ultraluminous infrared galaxies (ULIRGs) observed with the Spitzer Space Telescope's Infrared Spectrograph. Dust emission dominates the MIR spectra of ULIRGs, and the reprocessed radiation that emerges is independent of the underlying heating spectrum. Instead, the resulting emission depends sensitively on the geometric distribution of the dust, which… ▽ More We analyze the mid-infrared (MIR) spectra of ultraluminous infrared galaxies (ULIRGs) observed with the Spitzer Space Telescope's Infrared Spectrograph. Dust emission dominates the MIR spectra of ULIRGs, and the reprocessed radiation that emerges is independent of the underlying heating spectrum. Instead, the resulting emission depends sensitively on the geometric distribution of the dust, which we diagnose with comparisons of numerical simulations of radiative transfer. Quantifying the silicate emission and absorption features that appear near 10 and 18um requires a reliable determination of the continuum, and we demonstrate that including a measurement of the continuum at intermediate wavelength (between the features) produces accurate results at all optical depths. With high-quality spectra, we successfully use the silicate features to constrain the dust chemistry. The observations of the ULIRGs and local sightlines require dust that has a relatively high 18/10um absorption ratio of the silicate features (around 0.5). Specifically, the cold dust of Ossenkopf et al. (1992) is consistent with the observations, while other dust models are not. We use the silicate feature strengths to identify two families of ULIRGs, in which the dust distributions are fundamentally different. Optical spectral classifications are related to these families. In ULIRGs that harbor an active galactic nucleus, the spectrally broad lines are detected only when the nuclear surroundings are clumpy. In contrast, the sources of lower ionization optical spectra are deeply embedded in smooth distributions of optically thick dust. △ Less

Submitted 30 January, 2008; originally announced January 2008.

Comments: 38 pages, 11 figures; to appear in ApJ v679 (May 20)

PAH Emission from Ultraluminous Infrared Galaxies

Authors: V. Desai, L. Armus, H. W. W. Spoon, V. Charmandaris, J. Bernard-Salas, B. R. Brandl, D. Farrah, B. T. Soifer, H. I. Teplitz, P. M. Ogle, D. Devost, S. J. U. Higdon, J. A. Marshall, J. R. Houck

Abstract: We explore the relationships between the Polycyclic Aromatic Hydrocarbon (PAH) feature strengths, mid-infrared continuum luminosities, far-infrared spectral slopes, optical spectroscopic classifications, and silicate optical depths within a sample of 107 ULIRGs observed with the Infrared Spectrograph on the Spitzer Space Telescope. The detected 6.2 micron PAH equivalent widths (EQWs) in the samp… ▽ More We explore the relationships between the Polycyclic Aromatic Hydrocarbon (PAH) feature strengths, mid-infrared continuum luminosities, far-infrared spectral slopes, optical spectroscopic classifications, and silicate optical depths within a sample of 107 ULIRGs observed with the Infrared Spectrograph on the Spitzer Space Telescope. The detected 6.2 micron PAH equivalent widths (EQWs) in the sample span more than two orders of magnitude (0.006-0.8 micron), and ULIRGs with HII-like optical spectra or steep far-infrared spectral slopes (S_{25} / S_{60} < 0.2) typically have 6.2 micron PAH EQWs that are half that of lower-luminosity starbursts. A significant fraction (~40-60%) of HII-like, LINER-like, and cold ULIRGs have very weak PAH EQWs. Many of these ULIRGs also have large (tau_{9.7} > 2.3) silicate optical depths. The far-infrared spectral slope is strongly correlated with PAH EQW, but not with silicate optical depth. In addition, the PAH EQW decreases with increasing rest-frame 24 micron luminosity. We argue that this trend results primarily from dilution of the PAH EQW by continuum emission from dust heated by a compact central source, probably an AGN. High luminosity, high-redshift sources studied with Spitzer appear to have a much larger range in PAH EQW than seen in local ULIRGs, which is consistent with extremely luminous starburst systems being absent at low redshift, but present at early epochs. △ Less

Submitted 27 July, 2007; originally announced July 2007.

Comments: 15 pages, 9 Figures; Accepted for publication in ApJ

Decomposing Dusty Galaxies. I. Multi-Component Spectral Energy Distribution Fitting

Authors: J. A. Marshall, T. L. Herter, L. Armus, V. Charmandaris, H. W. W. Spoon, J. Bernard-Salas, J. R. Houck

Abstract: We present a new multi-component spectral energy distribution (SED) decomposition method and use it to analyze the ultraviolet to millimeter wavelength SEDs of a sample of dusty infrared-luminous galaxies. SEDs are constructed from spectroscopic and photometric data obtained with the Spitzer Space Telescope, in conjunction with photometry from the literature. Each SED is decomposed into emission… ▽ More We present a new multi-component spectral energy distribution (SED) decomposition method and use it to analyze the ultraviolet to millimeter wavelength SEDs of a sample of dusty infrared-luminous galaxies. SEDs are constructed from spectroscopic and photometric data obtained with the Spitzer Space Telescope, in conjunction with photometry from the literature. Each SED is decomposed into emission from populations of stars, an AGN accretion disk, PAHs, atomic and molecular lines, and distributions of graphite and silicate grains. Decompositions of the SEDs of the template starburst galaxies NGC7714 and NGC2623 and the template AGNs PG0804+761 and Mrk463 provide baseline properties to aid in quantifying the strength of star-formation and accretion in the composite systems NGC6240 and Mrk1014. We find that obscured radiation from stars is capable of powering the total dust emission from NGC6240, although we cannot rule out a contribution from a deeply embedded AGN visible only in X-rays. The decomposition of Mrk1014 is consistent with ~65% of its power emerging from an AGN and ~35% from star-formation. We suggest that many of the variations in our template starburst SEDs may be explained in terms of the different mean optical depths through the clouds of dust surrounding the young stars within each galaxy. Prompted by the divergent far-IR properties of our template AGNs, we suggest that variations in the relative orientation of their AGN accretion disks with respect to the disks of the galaxies hosting them may result in different amounts of AGN-heated cold dust emission emerging from their host galaxies. We estimate that 30-50% of the far-IR and PAH emission from Mrk1014 may originate from such AGN-heated material in its host galaxy disk. △ Less

Submitted 19 July, 2007; originally announced July 2007.

Comments: 27 pages, 12 figures. Accepted for publication in the ApJ

High resolution mid-infrared spectroscopy of ultraluminous infrared galaxies

Authors: D. Farrah, J. Bernard-Salas, H. W. W. Spoon, B. T. Soifer, L. Armus, B. Brandl, V. Charmandaris, V. Desai, S. Higdon, D. Devost, J. Houck

Abstract: (Abridged) We present R~600, 10-37um spectra of 53 ULIRGs at z<0.32, taken using the IRS on board Spitzer. All of the spectra show fine structure emission lines of Ne, O, S, Si and Ar, as well as molecular Hydrogen lines. Some ULIRGs also show emission lines of Cl, Fe, P, and atomic Hydrogen, and/or absorption features from C_2H_2, HCN, and OH. We employ diagnostics based on the fine-structure l… ▽ More (Abridged) We present R~600, 10-37um spectra of 53 ULIRGs at z<0.32, taken using the IRS on board Spitzer. All of the spectra show fine structure emission lines of Ne, O, S, Si and Ar, as well as molecular Hydrogen lines. Some ULIRGs also show emission lines of Cl, Fe, P, and atomic Hydrogen, and/or absorption features from C_2H_2, HCN, and OH. We employ diagnostics based on the fine-structure lines, as well as the EWs and luminosities of PAH features and the strength of the 9.7um silicate absorption feature (S_sil), to explore the power source behind the infrared emission in ULIRGs. We show that the IR emission from the majority of ULIRGs is powered mostly by star formation, with only ~20% of ULIRGs hosting an AGN with a comparable or greater IR luminosity than the starburst. The detection of the 14.32um [NeV] line in just under half the sample however implies that an AGN contributes significantly to the mid-IR flux in ~42% of ULIRGs. The emission line ratios, luminosities and PAH EWs are consistent with the starbursts and AGN in ULIRGs being more extincted, and for the starbursts more compact, versions of those in lower luminosity systems. The excitations and electron densities in the NLRs of ULIRGs appear comparable to those of lower luminosity starbursts, though there is evidence that the NLR gas in ULIRGs is more dense. We show that the combined luminosity of the 12.81um [NeII] and 15.56um [NeIII] lines correlates with both IR luminosity and the luminosity of the 6.2 micron and 11.2 micron PAH features in ULIRGs, and use this to derive a calibration between PAH luminosity and star formation rate. Finally, we show that ULIRGs with 0.8 < S_sil < 2.4 are likely to be powered mainly by star formation, but that ULIRGs with S_sil < 0.8, and possibly those with S_sil > 2.4, contain an IR-luminous AGN. △ Less

Submitted 4 June, 2007; originally announced June 2007.

Comments: 62 pages in preprint format, 4 tables, 23 figures. ApJ accepted

Journal ref: Astrophys.J.667:149-169,2007

Infrared Molecular Starburst Fingerprints in Deeply Obscured (U)LIRG Nuclei

Authors: F. Lahuis, H. W. W. Spoon, A. G. G. M. Tielens, S. D. Doty, L. Armus, V. Charmandaris, J. R. Houck, P. Stäuber, E. F. van Dishoeck

Abstract: High resolution spectra of the Spitzer Space Telescope show vibration-rotation absorption bands of gaseous C2H2, HCN, and CO2 molecules toward a sample of deeply obscured (U)LIRG nuclei. The observed bands reveal the presence of dense (n>~ 10^7 cm^-3), warm (T_ex = 200-700 K) molecular gas with high column densities of these molecules ranging from a few 10^15 - 10^17 cm^-2. Abundances relative t… ▽ More High resolution spectra of the Spitzer Space Telescope show vibration-rotation absorption bands of gaseous C2H2, HCN, and CO2 molecules toward a sample of deeply obscured (U)LIRG nuclei. The observed bands reveal the presence of dense (n>~ 10^7 cm^-3), warm (T_ex = 200-700 K) molecular gas with high column densities of these molecules ranging from a few 10^15 - 10^17 cm^-2. Abundances relative to H2, inferred from the silicate optical depth, range from ~10^-7 to 10^-6 and show no correlation with temperature. Theoretical studies show that the high abundances of both C2H2 and HCN exclude a X-ray dominated region (XDR) associated with the toroid surrounding an AGN as the origin of this dense warm molecular gas. Galactic massive protostars in the so-called Hot Core phase have similar physical characteristics with comparable high abundances of C2H2, HCN, and CO2 in the hot phase. However, the abundances of C2H2 and HCN and the C2H2/CO2 and HCN/CO2 ratios are much higher toward the (U)LIRGs in the cooler (T_ex <= 400 K) phase. We suggest that the warm dense molecular gas revealed by the mid-IR absorption lines is associated with a phase of deeply embedded star formation where the extreme pressures and densities of the nuclear starburst environment have inhibited the expansion of HII regions and the global disruption of the star forming molecular cloud cores, and `trapped' the star formation process in an `extended' Hot Core phase. △ Less

Submitted 28 June, 2007; v1 submitted 28 December, 2006; originally announced December 2006.

Comments: 22 pages, 4 figures, published in ApJ, 659, 296 v2. updated hydrogen column densities

Journal ref: Astrophys.J.659:296-304,2007

The distribution of silicate strength in Spitzer spectra of AGNs and ULIRGs

Authors: Lei Hao, D. W. Weedman, H. W. W. Spoon, J. A. Marshall, N. A. Levenson, M. Elitzur, J. R. Houck

Abstract: A sample of 196 AGNs and ULIRGs observed by the Infrared Spectrograph (IRS) on Spitzer is analyzed to study the distribution of the strength of the 9.7 micron silicate feature. Average spectra are derived for quasars, Seyfert 1 and Seyfert 2 AGNs, and ULIRGs. We find that quasars are characterized by silicate features in emission and Seyfert 1s equally by emission or weak absorption. Seyfert 2s… ▽ More A sample of 196 AGNs and ULIRGs observed by the Infrared Spectrograph (IRS) on Spitzer is analyzed to study the distribution of the strength of the 9.7 micron silicate feature. Average spectra are derived for quasars, Seyfert 1 and Seyfert 2 AGNs, and ULIRGs. We find that quasars are characterized by silicate features in emission and Seyfert 1s equally by emission or weak absorption. Seyfert 2s are dominated by weak silicate absorption, and ULIRGs are characterized by strong silicate absorption (mean apparent optical depth about 1.5). Luminosity distributions show that luminosities at rest frame 5.5 micron are similar for the most luminous quasars and ULIRGs and are almost 10^5 times more luminous than the least luminous AGN in the sample. The distributions of spectral characteristics and luminosities are compared to those of optically faint infrared sources at z~2 being discovered by the IRS, which are also characterized by strong silicate absorption. It is found that local ULIRGs are a similar population, although they have lower luminosities and somewhat stronger absorption compared to the high redshift sources. △ Less

Submitted 18 December, 2006; originally announced December 2006.

Comments: Accepted for publication on ApJL

Journal ref: Astrophys.J.655:L77-L80,2007

Mid-IR Galaxy Classification Based on Silicate Obscuration and PAH Equivalent Width

Authors: H. W. W. Spoon, J. A. Marshall, J. R. Houck, M. Elitzur, L. Hao, L. Armus, B. R. Brandl, V. Charmandaris

Abstract: We present a new diagnostic diagram for mid-infrared spectra of infrared galaxies based on the equivalent width of the 6.2 micron PAH emission feature and the strength of the 9.7 micron silicate feature. Based on the position in this diagram we classify galaxies into 9 classes ranging from continuum-dominated AGN hot dust spectra and PAH-dominated starburst spectra to absorption-dominated spectr… ▽ More We present a new diagnostic diagram for mid-infrared spectra of infrared galaxies based on the equivalent width of the 6.2 micron PAH emission feature and the strength of the 9.7 micron silicate feature. Based on the position in this diagram we classify galaxies into 9 classes ranging from continuum-dominated AGN hot dust spectra and PAH-dominated starburst spectra to absorption-dominated spectra of deeply obscured galactic nuclei. We find that galaxies are systematically distributed along two distinct branches: one of AGN and starburst-dominated spectra and one of deeply obscured nuclei and starburst-dominated spectra. The separation into two branches likely reflects a fundamental difference in the dust geometry in the two sets of sources: clumpy versus non-clumpy obscuration. Spectra of ULIRGs are found along the full length of both branches, reflecting the diverse nature of the ULIRG family. △ Less

Submitted 29 November, 2006; originally announced November 2006.

Comments: 5 pages, 3 figures. Accepted for publication in ApJ Letters

Journal ref: Astrophys.J.654:L49-L52,2006

Deep Mid-Infrared Silicate Absorption as a Diagnostic of Obscuring Geometry Toward Galactic Nuclei

Authors: N. A. Levenson, M. M. Sirocky, L. Hao, H. W. W. Spoon, J. A. Marshall, M. Elitzur, J. R. Houck

Abstract: The silicate cross section peak near 10um produces emission and absorption features in the spectra of dusty galactic nuclei observed with the Spitzer Space Telescope. Especially in ultraluminous infrared galaxies, the observed absorption feature can be extremely deep, as IRAS 08572+3915 illustrates. A foreground screen of obscuration cannot reproduce this observed feature, even at large optical… ▽ More The silicate cross section peak near 10um produces emission and absorption features in the spectra of dusty galactic nuclei observed with the Spitzer Space Telescope. Especially in ultraluminous infrared galaxies, the observed absorption feature can be extremely deep, as IRAS 08572+3915 illustrates. A foreground screen of obscuration cannot reproduce this observed feature, even at large optical depth. Instead, the deep absorption requires a nuclear source to be deeply embedded in a smooth distribution of material that is both geometrically and optically thick. In contrast, a clumpy medium can produce only shallow absorption or emission, which are characteristic of optically-identified active galactic nuclei. In general, the geometry of the dusty region and the total optical depth, rather than the grain composition or heating spectrum, determine the silicate feature's observable properties. The apparent optical depth calculated from the ratio of line to continuum emission generally fails to accurately measure the true optical depth. The obscuring geometry, not the nature of the embedded source, also determines the far-IR spectral shape. △ Less

Submitted 14 November, 2006; originally announced November 2006.

Comments: To appear in ApJL

Journal ref: Astrophys.J.654:L45-L48,2007

Observations of Ultraluminous Infrared Galaxies with the Infrared Spectrograph on the Spitzer Space Telescope II: The IRAS Bright Galaxy Sample

Authors: L. Armus, V. Charmandaris, J. Bernard-Salas, H. W. W. Spoon, J. A. Marshall, S. J. U Higdon, V. Desai, H. I. Teplitz, L. Hao, D. Devost, B. R. Brandl, Y. Wu, G. C. Sloan, B. T. Soifer, J. R. Houck, T. L. Herter

Abstract: We present spectra taken with the Infrared Spectrograph on Spitzer covering the 5-38 micron region of the ten Ultraluminous Infrared Galaxies (ULIRGs) found in the IRAS Bright Galaxy Sample. Among the BGS ULIRGs, we find a factor of 50 spread in the rest-frame mid to far-infrared spectral slope. The 9.7 micron silicate optical depths range from less than 0.4 more than 4.2, implying line of sight… ▽ More We present spectra taken with the Infrared Spectrograph on Spitzer covering the 5-38 micron region of the ten Ultraluminous Infrared Galaxies (ULIRGs) found in the IRAS Bright Galaxy Sample. Among the BGS ULIRGs, we find a factor of 50 spread in the rest-frame mid to far-infrared spectral slope. The 9.7 micron silicate optical depths range from less than 0.4 more than 4.2, implying line of sight extinctions of A(V) ~ 8 - 78 mag. There is evidence for water ice and hydrocarbon absorption and C2H2 and HCN absorption features in four and possibly six of the 10 BGS ULIRGs, indicating shielded molecular clouds and a warm, dense ISM. We have detected [NeV] emission in three of the ten BGS ULIRGs, at flux levels of 5-18E-14 erg/cm^2/sec and [NeV] 14.3/[NeII] 12.8 line flux ratios of 0.12-0.85. The remaining BGS ULIRGs have limits on their [NeV]/[NeII] line flux ratios which range from less than 0.15 to less than 0.01. Among the BGS ULIRGs, the AGN fractions implied by either the [NeV]/[NeII] or [OIV]/[NeII] line flux ratios (or their upper limits) are significantly lower than implied by the MIR slope or EQW of the 6.2 micron PAH feature. Fitting the SEDs, we see evidence for hot (T > 300K) dust in five of the BGS ULIRGs, with the fraction of hot dust to total dust luminosity ranging from ~1-23%, before correcting for extinction. When integrated over the IRAC-8, IRS blue peakup, and MIPS-24 filter bandpasses, the IRS spectra imply very blue colors for some ULIRGs at z ~ 1.3. This is most extreme for sources with significant amounts of warm dust and deep silicate absorption. △ Less

Submitted 6 October, 2006; originally announced October 2006.

Comments: accepted for publication in the Astrophysical Journal

Journal ref: Astrophys.J.656:148-167,2007

The Mid-IR Properties of Starburst Galaxies from Spitzer-IRS Spectroscopy

Authors: B. R. Brandl, J. Bernard-Salas, H. W. W. Spoon, D. Devost, G. C. Sloan, S. Guilles, Y. Wu, J. R. Houck, L. Armus, D. W. Weedman, V. Charmandaris, P. N. Appleton, B. T. Soifer, L. Hao, J. A. Marshall, S. J. Higdon, T. L. Herter

Abstract: We present 5-38um mid-infrared spectra at a spectral resolution of R~65-130 of a large sample of 22 starburst nuclei taken with the Infrared Spectrograph IRS on board the Spitzer Space Telescope. The spectra show a vast range in starburst SEDs. The silicate absorption ranges from essentially no absorption to heavily obscured systems with an optical depth of tau(9.8um)~5. The spectral slopes can… ▽ More We present 5-38um mid-infrared spectra at a spectral resolution of R~65-130 of a large sample of 22 starburst nuclei taken with the Infrared Spectrograph IRS on board the Spitzer Space Telescope. The spectra show a vast range in starburst SEDs. The silicate absorption ranges from essentially no absorption to heavily obscured systems with an optical depth of tau(9.8um)~5. The spectral slopes can be used to discriminate between starburst and AGN powered sources. The monochromatic continuum fluxes at 15um and 30um enable a remarkably accurate estimate of the total infrared luminosity of the starburst. We find that the PAH equivalent width is independent of the total starburst luminosity L_IR as both continuum and PAH feature scale proportionally. However, the luminosity of the 6.2um feature scales with L_IR and can be used to approximate the total infrared luminosity of the starburst. Although our starburst sample covers about a factor of ten difference in the [NeIII]/[NeII] ratio, we found no systematic correlation between the radiation field hardness and the PAH equivalent width or the 7.7um/11.3um PAH ratio. These results are based on spatially integrated diagnostics over an entire starburst region, and local variations may be ``averaged out''. It is presumably due to this effect that unresolved starburst nuclei with significantly different global properties appear spectrally as rather similar members of one class of objects. △ Less

Submitted 1 September, 2006; originally announced September 2006.

Comments: 22 pages, accepted for publication in ApJ, a high-resolution version is available from http://www.strw.leidenuniv.nl/~brandl/IRS_starbursts.pdf

Journal ref: Astrophys.J.653:1129-1144,2006; Erratum-ibid.665:884-,2007

A Spitzer Infrared Spectrograph Survey of Warm Molecular Hydrogen in Ultra-luminous Infrared Galaxies

Authors: S. J. U. Higdon, L. Armus, J. L. Higdon, B. T. Soifer, H. W. W. Spoon

Abstract: We have conducted a survey of Ultra-luminous Infrared Galaxies (ULIRGs) with the Infrared Spectrograph on the Spitzer Space Telescope, obtaining spectra from 5.0-38.5um for 77 sources with 0.02<z <0.93. Observations of the pure rotational H2 lines S(3) 9.67um, S(2) 12.28um, and S(1) 17.04um are used to derive the temperature and mass of the warm molecular gas. We detect H2 in 77% of the sample,… ▽ More We have conducted a survey of Ultra-luminous Infrared Galaxies (ULIRGs) with the Infrared Spectrograph on the Spitzer Space Telescope, obtaining spectra from 5.0-38.5um for 77 sources with 0.02<z <0.93. Observations of the pure rotational H2 lines S(3) 9.67um, S(2) 12.28um, and S(1) 17.04um are used to derive the temperature and mass of the warm molecular gas. We detect H2 in 77% of the sample, and all ULIRGs with F(60um)>2Jy. The average warm molecular gas mass is ~2x10^8solar-masses. High extinction, inferred from the 9.7um silicate absorption depth, is not observed along the line of site to the molecular gas. The derived H2 mass does not depend on F(25um)/F(60um), which has been used to infer either starburst or AGN dominance. Similarly, the molecular mass does not scale with the 25 or 60um luminosities. In general, the H2 emission is consistent with an origin in photo-dissociation regions associated with star formation. We detect the S(0) 28.22um emission line in a few ULIRGs. Including this line in the model fits tends to lower the temperature by ~50-100K, resulting in a significant increase in the gas mass. The presence of a cooler component cannot be ruled out in the remainder of our sample, for which we do not detect the S(0) line. The measured S(7) 5.51um line fluxes in six ULIRGs implies ~3x10^6 solar-masses of hot (~1400K) H2. The warm gas mass is typically less than 1% of the cold gas mass derived from CO observations. △ Less

Submitted 14 June, 2006; v1 submitted 15 May, 2006; originally announced May 2006.

Comments: Accepted ApJ 01 September 2006, v648n1 issue. 14 pages 12 figures IRAS 06361-6217 the f25/f60 ratio is 0.10 not 1.00

IRS spectra of two ultraluminous infrared galaxies at z=1.3

Authors: V. Desai, L. Armus, B. T. Soifer, D. W. Weedman, S. Higdon, C. Bian, C. Borys, H. W. W. Spoon, V. Charmandaris, K. Brand, M. J. I. Brown, A. Dey, J. Higdon, J. Houck, B. T. Jannuzi, E. Le Floc'h, M. L. N. Ashby, H. A. Smith

Abstract: We present low-resolution (64 < R < 124) mid-infrared (8--38 micron) Spitzer/IRS spectra of two z~1.3 ultraluminous infrared galaxies (LFIR~10^13) discovered in a Spitzer/MIPS survey of the Bootes field of the NOAO Deep Wide-Field Survey (NDWFS). MIPS J142824.0+352619 is a bright 160 micron source with a large infrared-to-optical flux density ratio and a possible lensing amplification of <~10. T… ▽ More We present low-resolution (64 < R < 124) mid-infrared (8--38 micron) Spitzer/IRS spectra of two z~1.3 ultraluminous infrared galaxies (LFIR~10^13) discovered in a Spitzer/MIPS survey of the Bootes field of the NOAO Deep Wide-Field Survey (NDWFS). MIPS J142824.0+352619 is a bright 160 micron source with a large infrared-to-optical flux density ratio and a possible lensing amplification of <~10. The 6.2, 7.7, 11.3, and 12.8 micron PAH emission bands in its IRS spectrum indicate a redshift of z~1.3. The large equivalent width of the 6.2 micron PAH feature indicates that at least 50% of the mid-infrared energy is generated in a starburst, an interpretation that is supported by a large [NeII]/[NeIII] ratio and a low upper limit on the X-ray luminosity. SST24 J142827.19+354127.71 has the brightest 24 micron flux (10.55 mJy) among optically faint (R > 20) galaxies in the NDWFS. Its mid-infrared spectrum lacks emission features, but the broad 9.7 micron silicate absorption band places this source at z~1.3. Given this redshift, SST24 J142827.19+354127.71 has among the largest rest-frame 5 micron luminosities known. The similarity of its SED to those of known AGN-dominated ULIRGs and its lack of either PAH features or large amounts of cool dust indicate that the powerful mid-infrared emission is dominated by an active nucleus rather than a starburst. Our results illustrate the power of the IRS in identifying massive galaxies in the ``redshift desert'' and in discerning their power sources. Because they are bright, MIPS J142824.0+352619 (pending future observations to constrain its lensing amplification) and SST24 J142827.19+354127.71 are useful z>1 templates of a high luminosity starburst and AGN, respectively. △ Less

Submitted 13 December, 2005; v1 submitted 8 December, 2005; originally announced December 2005.

Comments: 9 pages, 4 figures; accepted for publication in ApJ; several small clarifications have been made

Journal ref: Astrophys.J.641:133-139,2006