Mid-infrared molecular hydrogen (H$_2$) emission is a powerful cooling agent in galaxy mergers and in radio galaxies; it is a potential key tracer of gas evolution and energy dissipation associated with mergers, star formation, and accretion onto supermassive black holes. We detect mid-IR H$_2$ line emission in at least one rotational transition in 91\% of the 214 Luminous Infrared Galaxies (LIRGs) observed with Spitzer as part of the Great Observatories All-sky LIRG Survey (GOALS). We use H$_2$ excitation diagrams to estimate the range of masses and temperatures of warm molecular gas in these galaxies. We find that LIRGs in which the IR emission originates mostly from the Active Galactic Nuclei (AGN) have about 100K higher H$_2$ mass-averaged excitation temperatures than LIRGs in which the IR emission originates mostly from star formation. Between 10 and 15\% of LIRGs have H$_2$ emission lines that are sufficiently broad to be resolved or partially resolved by the high resolution modules of Spitzer's Infrared Spectrograph (IRS). Those sources tend to be mergers and contain AGN. This suggests that a significant fraction of the H$_2$ line emission is powered by AGN activity through X-rays, cosmic rays, and turbulence. We find a statistically significant correlation between the kinetic energy in the H$_2$ gas and the H$_2$ to IR luminosity ratio. The sources with the largest warm gas kinetic energies are mergers. We speculate that mergers increase the production of bulk in-flows leading to observable broad H$_2$ profiles and possibly denser environments.

The true importance of the warm, AGN-driven outflows for the evolution of galaxies remains uncertain. Measurements of the radial extents of the outflows are key for quantifying their masses and kinetic powers, and also establishing whether the AGN outflows are galaxy-wide. Therefore, as part of a larger project to investigate the significance of warm, AGN-driven outflows in the most rapidly evolving galaxies in the local universe, here we present deep Hubble Space Telescope ( HST) narrow-band [OIII]$\lambda$5007 observations of a complete sample of 8 nearby ULIRGs with optical AGN nuclei. Combined with the complementary information provided by our ground-based spectroscopy, the HST images show that the warm gas outflows are relatively compact for most of the objects in the sample: in three objects the outflow regions are barely resolved at the resolution of HST ($0.065 < R_{[OIII]} < 0.12$ kpc); in a further four cases the outflows are spatially resolved but with flux weighted mean radii in the range $0.65 < R_{[OIII]} < 1.2$ kpc; and in only one object (Mrk273) is there clear evidence for a more extended outflow, with a maximum extent of $R_{[OIII]}\sim5$ kpc. Overall, our observations show little evidence for the galaxy-wide outflows predicted by some models of AGN feedback.

We report molecular gas observations of IRAS 20100-4156 and IRAS 03158+4227, two local ultraluminous infrared galaxies (ULIRGs) hosting some of the fastest and most massive molecular outflows known. Using ALMA and PdBI observations, we spatially resolve the CO(1-0) emission from the outflowing molecular gas in both and find maximum outflow velocities of $ v_{\rm max} \sim 1600$ and $\sim 1700$ km/s for IRAS 20100-4156 and IRAS 03158+4227, respectively. We find total gas mass outflow rates of $\dot M_{\rm OF} \sim 670$ and $\sim 350$ Msun/yr, respectively, corresponding to molecular gas depletion timescales $\tau^{\rm dep}_{\rm OF} \sim 11$ and $\sim 16$ Myr. This is nearly 3 times shorter than the depletion timescales implied by star formation, $\tau^{\rm dep}_{\rm SFR} \sim 33$ and $\sim 46$ Myr, respectively. To determine the outflow driving mechanism, we compare the starburst ($L_{*}$) and AGN ($L_{\rm AGN}$) luminosities to the outflowing energy and momentum fluxes, using mid-infrared spectral decomposition to discern $L_{\rm AGN}$. Comparison to other molecular outflows in ULIRGs reveals that outflow properties correlate similarly with $L_{*}$ and $L_{\rm IR}$ as with $L_{\rm AGN}$, indicating that AGN luminosity alone may not be a good tracer of feedback strength and that a combination of AGN and starburst activity may be driving the most powerful molecular outflows. We also detect the OH 1.667 GHz maser line from both sources and demonstrate its utility in detecting molecular outflows.

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

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.

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.

The galaxy NGC 4418 harbours a compact ($<20$ pc) core with a very high bolometric luminosity ($\sim10^{11}$L$_\odot$). As most of the galaxy's energy output comes from this small region, it is of interest to determine what fuels this intense activity. An interaction with VV 655 has been proposed, where gas aquired by NGC 4418 could trigger intense star formation and/or black hole accretion in the centre. We aim to constrain the interaction hypothesis by studying neutral hydrogen structures around the two galaxies. We present observations at 1.4 GHz with the Very Large Array of radio continuum as well as emission and absorption from atomic hydrogen. Gaussian distributions are fitted to observed HI emission and absorption spectra. An atomic HI bridge is seen in emission, connecting NGC 4418 to VV 655. While NGC 4418 is bright in continuum emission and seen in HI absorption, VV 655 is barely detected in the continuum but show bright HI emission (M$_\mathrm{HI}\sim10^9$ M$_\odot$). We estimate SFRs from 1.4 GHz of 3.2 M$_\odot$ yr$^{-1}$ and 0.13 M$_\odot$ yr$^{-1}$ for NGC 4418 and VV 655 respectively. Systemic HI velocities of 2202$\pm$20 km s$^{-1}$ (emission) and 2105.4$\pm$10 km s$^{-1}$ (absorption) are measured for VV 655 and NGC 4418 respectively. Redshifted HI absorption is seen towards NGC 4418, suggesting gas infall. Blueshifted HI-emission is seen north-west of NGC 4418, which we interpret as a continuation of the outflow previously discussed by Sakamoto et al. (2013). The morphology and velocity structure seen in HI is consistent with an interaction scenario, where gas was transferred from VV 655 to NGC 4418, and may fuel the activity in the centre.

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.

We present the detection of four far-infrared fine-structure oxygen lines, as well as strong upper limits for the CO(2-1) and [N II] 205 um lines, in 3C 368, a well-studied radio-loud galaxy at z = 1.131. These new oxygen lines, taken in conjunction with previously observed neon and carbon fine-structure lines, suggest a powerful active galactic nucleus (AGN), accompanied by vigorous and extended star formation. A starburst dominated by O8 stars, with an age of ~6.5 Myr, provides a good fit to the fine-structure line data. This estimated age of the starburst makes it nearly concurrent with the latest episode of AGN activity, suggesting a link between the growth of the supermassive black hole and stellar population in this source. We do not detect the CO(2-1) line, down to a level twelve times lower than the expected value for star forming galaxies. This lack of CO line emission is consistent with recent star formation activity if the star-forming molecular gas has low metallicity, is highly fractionated (such that CO is photodissociated through much of the clouds), or is chemically very young (such that CO has not yet had time to form). It is also possible, though we argue unlikely, that the ensemble of fine structure lines are emitted from the region heated by the AGN.

We present the Herschel SPIRE Fourier Transform Spectroscopy (FTS) atlas for a complete flux limited sample of local Ultra-Luminous Infra-Red Galaxies as part of the HERschel ULIRG Survey (HERUS). The data reduction is described in detail and was optimized for faint FTS sources with particular care being taken with the subtraction of the background which dominates the continuum shape of the spectra. Special treatment in the data reduction has been given to any observation suffering from artefacts in the data caused by anomalous instrumental effects to improve the final spectra. Complete spectra are shown covering $200 - 671\mu$m with photometry in the SPIRE bands at 250$\mu$m, 350$\mu$m and 500$\mu$m. The spectra include near complete CO ladders for over half of our sample, as well as fine structure lines from [CI] 370 $\mu$m, [CI] 609 $\mu$m, and [NII] 205 $\mu$m. We also detect H$_{2}$O lines in several objects. We construct CO Spectral Line Energy Distributions (SLEDs) for the sample, and compare their slopes with the far-infrared colours and luminosities. We show that the CO SLEDs of ULIRGs can be broadly grouped into three classes based on their excitation. We find that the mid-J (5$<$J$<$8) lines are better correlated with the total far-infrared luminosity, suggesting that the warm gas component is closely linked to recent star-formation. The higher J transitions do not linearly correlate with the far-infrared luminosity, consistent with them originating in hotter, denser gas unconnected to the current star-formation. {\bf We conclude that in most cases more than one temperature components are required to model the CO SLEDs.}

We present deep new GTC/OSIRIS narrow-band images and optical WHT/ISIS long-slit spectroscopy of the merging system Mrk273 that show a spectacular extended halo of warm ionised gas out to a radius of $\sim45$ kpc from the system nucleus. Outside of the immediate nuclear regions (r > 6 kpc), there is no evidence for kinematic disturbance in the ionised gas: in the extended regions covered by our spectroscopic slits the emission lines are relatively narrow (FWHM $\lesssim$ 350 km$\rm s^{-1}$) and velocity shifts small (|$\Delta$V| $\lesssim{} $250 km$\rm s^{-1}$). This is despite the presence of powerful near-nuclear outflows (FWHM > 1000 km$\rm s^{-1}$; |$\Delta$V| > 400 km$\rm s^{-1}$; r < 6 kpc). Diagnostic ratio plots are fully consistent with Seyfert 2 photo-ionisation to the NE of the nuclear region, however to the SW the plots are more consistent with low-velocity radiative shock models. The kinematics of the ionised gas, combined with the fact that the main structures are aligned with low-surface-brightness tidal continuum features, are consistent with the idea that the ionised halo represents tidal debris left over from a possible triple-merger event, rather than a reservoir of outflowing gas.

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 ({\sigma}(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.

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.

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.

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.

We have recently detected the [CII] 157.7 micron line in eight star forming galaxies at redshifts 1 to 2 using the redshift(z) Early Universe Spectrometer (ZEUS). Our sample targets star formation dominant sources detected in PAH emission. This represents a significant addition to [CII] observations during the epoch of peak star formation. We have augmented this survey with observations of the [OI] 63 micron line and far infrared photometry from the PACS and SPIRE Herschel instruments as well as Spitzer IRS spectra from the literature showing PAH features. Our sources exhibit above average gas heating efficiency, many with both [OI]/FIR and [CII]/FIR ~1% or more. The relatively strong [CII] emission is consistent with our sources being dominated by star formation powered PDRs, extending to kpc scales. We suggest that the star formation mode in these systems follows a Schmidt-Kennicutt law similar to local systems, but at a much higher rate due to molecular gas surface densities 10 to 100 times that of local star forming systems. The source of the high molecular gas surface densities may be the infall of neutral gas from the cosmic web. In addition to the high [CII]/FIR values, we also find high [CII]/PAH ratios and, in at least one source, a cool dust temperature. This source, SWIRE 4-5, bears a resemblance in these diagnostics to shocked regions of Stephan's Quintet, suggesting that another mode of [CII] excitation in addition to normal photoelectric heating may be contributing to the observed [CII] line.

It remains uncertain which continuum and emission line diagnostics best indicate the bolometric powers of active galactic nuclei (AGNs), especially given the attenuation caused by the circumnuclear material and the possible contamination by components related to star formation. Here we use mid-IR spectra along with multiwavelength data to investigate the merit of various diagnostics of AGN radiative power, including the mid-IR [Ne III] lambda25.89 mum and [O IV] lambda25.89 mum fine-structure lines, the optical [O III] lambda5007 forbidden line, and mid-IR 24 mum, 5 GHz radio, and X-ray continuum emission, for complete samples of 46 2Jy radio galaxies (0.05 < z < 0.7) and 17 3CRR FRII radio galaxies (z < 0.1). We find that the mid-IR [O IV] line is the most reliable indicator of AGN power for powerful radio-loud AGNs. By assuming that the [O IV] is emitted isotropically, and comparing the [O III] and 24 mum luminosities of the broad- and narrow-line AGNs in our samples at fixed [O IV] luminosity, we show that the [O III] and 24 mum emission are both mildly attenuated in the narrow-line compared to the broad-line objects by a factor of #2. However, despite this attenuation, the [O III] and 24 mum luminosities are better AGN power indicators for our sample than either the 5 GHz radio or the X-ray continuum luminosities. We also detect the mid-IR 9.7 mum silicate feature in the spectra of many objects but not ubiquitously: at least 40% of the sample shows no clear evidence for these features. We conclude that, for the majority of powerful radio galaxies, the mid-IR lines are powered by AGN photoionization.

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.

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.

A summary is presented for 130 galaxies observed with the Herschel PACS instrument to measure fluxes for the [CII] 158 um emission line. Sources cover a wide range of active galactic nucleus to starburst classifications, as derived from polycyclic aromatic hydrocarbon (PAH) strength measured with the Spitzer Infrared Spectrograph. Redshifts from [CII] and line to continuum strengths (equivalent width of [CII]) are given for the full sample, which includes 18 new [CII] flux measures. Calibration of L([CII)]) as a star formation rate (SFR) indicator is determined by comparing [CII] luminosities with mid-infrared [NeII] and [NeIII] emission line luminosities; this gives the same result as determining SFR using bolometric luminosities of reradiating dust from starbursts: log SFR = log L([CII)]) - 7.0, for SFR in solar masses per year and L([CII]) in solar luminosities. We conclude that L([CII]) can be used to measure SFR in any source to a precision of ~ 50%, even if total source luminosities are dominated by an AGN component. The line to continuum ratio at 158 um, EW([CII]), is not significantly greater for starbursts (median EW([CII]) = 1.0 um) compared to composites and AGN (median EW([CII]) = 0.7 um), showing that the far infrared continuum at 158 um scales with [CII] regardless of classification. This indicates that the continuum at 158 um also arises primarily from the starburst component within any source, giving log SFR = log vLv(158 um) - 42.8 for SFR in solar masses per year and vLv(158 um) in erg per sec.

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.

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/}

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.

We present the data and our analysis of MIR fine-structure emission lines detected in Spitzer/IRS high-res spectra of 202 local LIRGs observed as part of the GOALS project. We detect emission lines of [SIV], [NeII], [NeV], [NeIII], [SIII]18.7, [OIV], [FeII], [SIII]33.5, and [SiII]. Over 75% of our galaxies are classified as starburst (SB) sources in the MIR. We compare ratios of the emission line fluxes to stellar photo- and shock-ionization models to constrain the gas properties in the SB nuclei. Comparing the [SIV]/[NeII] and [NeIII]/[NeII] ratios to the Starburst99-Mappings III models with an instantaneous burst history, the line ratios suggest that the SB in our LIRGs have ages of 1-4.5Myr, metallicities of 1-2Z_sun, and ionization parameters of 2-8e7cm/s. Based on the [SIII]/[SIII] ratios, the electron density in LIRG nuclei has a median electron density of ~300cm-3 for sources above the low density limit. We also find that strong shocks are likely present in 10 SB sources. A significant fraction of the GOALS sources have resolved neon lines and 5 show velocity differences of >200km/s in [NeIII] or [NeV] relative to [NeII]. Furthermore, 6 SB and 5 AGN LIRGs show a trend of increasing line width with ionization potential, suggesting the possibility of a compact energy source and stratified ISM in their nuclei. We confirm a strong correlation between the [NeII]+[NeIII] emission, as well as [SIII]33.5, with both the IR luminosity and the 24um warm dust emission measured from the spectra. Finally, we find no correlation between the hardness of the radiation field or the line width and the ratio of the total IR to 8um emission (IR8). This may be because the IR luminosity and the MIR fine-structure lines are sensitive to different timescales over the SB, or that IR8 is more sensitive to the geometry of the warm dust region than the radiation field producing the HII region emission.

We present Herschel observations of six fine-structure lines in 25 Ultraluminous Infrared Galaxies at z<0.27. The lines, [O III]52, [N III]57, [O I]63, [N II]122, [O I]145, and [C II]158, are mostly single gaussians with widths <600 km s-1 and luminosities of 10^7 - 10^9 Solar. There are deficits in the [O I]63/L_IR, [N II]/L_IR, [O I]145/L_IR, and [C II]/L_IR ratios compared to lower luminosity systems. The majority of the line deficits are consistent with dustier H II regions, but part of the [C II] deficit may arise from an additional mechanism, plausibly charged dust grains. This is consistent with some of the [C II] originating from PDRs or the ISM. We derive relations between far-IR line luminosities and both IR luminosity and star formation rate. We find that [N II] and both [O I] lines are good tracers of IR luminosity and star formation rate. In contrast, [C II] is a poor tracer of IR luminosity and star formation rate, and does not improve as a tracer of either quantity if the [C II] deficit is accounted for. The continuum luminosity densities also correlate with IR luminosity and star formation rate. We derive ranges for the gas density and ultraviolet radiation intensity of 10^1 < n < 10^2.5 and 10^2.2 < G_0 < 10^3.6, respectively. These ranges depend on optical type, the importance of star formation, and merger stage. We do not find relationships between far-IR line properties and several other parameters; AGN activity, merger stage, mid-IR excitation, and SMBH mass. We conclude that these far-IR lines arise from gas heated by starlight, and that they are not strongly influenced by AGN activity.

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.

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.

We present mid-infrared spectra and photometry of thirteen redshift 0.4<z<1 dust-reddened quasars obtained with Spitzer IRS and MIPS. We compare properties derived from their infrared spectral energy distributions (intrinsic AGN luminosity and far-infrared luminosity from star formation) to the host luminosities and morphologies from HST imaging, and black hole masses estimated from optical and/or near-infrared spectroscopy. Our results are broadly consistent with models in which most dust reddened quasars are an intermediate phase between a merger-driven starburst triggering a completely obscured AGN, and a normal, unreddened quasar. We find that many of our objects have high accretion rates, close to the Eddington limit. These objects tend to fall below the black hole mass -- bulge luminosity relation as defined by local galaxies, whereas most of our low accretion rate objects are slightly above the local relation, as typical for normal quasars at these redshifts. Our observations are therefore most readily interpreted in a scenario in which galaxy stellar mass growth occurs first by about a factor of three in each merger/starburst event, followed sometime later by black hole growth by a similar amount. We do not, however, see any direct evidence for quasar feedback affecting star formation in our objects, for example in the form of a relationship between accretion rate and star formation. Five of our objects, however, do show evidence for outflows in the OIII 5007 Angstrom emission line profile, suggesting that the quasar activity is driving thermal winds in at least some members of our sample.

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.

Results are presented for [CII] 158 micron line fluxes observed with the Herschel PACS instrument in 112 sources with both starburst and AGN classifications, of which 102 sources have confident detections. Results are compared with mid-infrared spectra from the Spitzer Infrared Spectrometer and with L(IR) from IRAS fluxes; AGN/starburst classifications are determined from equivalent width of the 6.2 micron PAH feature. It is found that the [CII] line flux correlates closely with the flux of the 11.3 micron PAH feature independent of AGN/starburst classification, log [f([CII] 158 micron)/f(11.3 micron PAH)] = -0.22 +- 0.25. It is concluded that [CII] line flux measures the photodissociation region associated with starbursts in the same fashion as the PAH feature. A calibration of star formation rate for the starburst component in any source having [CII] is derived comparing [CII] luminosity L([CII]) to L(IR) with the result that log SFR = log L([CII)]) - 7.08 +- 0.3, for SFR in solar masses per year and L([CII]) in solar luminosities. The decreasing ratio of L([CII]) to L(IR) in more luminous sources (the "[CII] deficit") is shown to be a consequence of the dominant contribution to L(IR) arising from a luminous AGN component because the sources with largest L(IR) and smallest L([CII])/L(IR) are AGN.

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.

We present deep Spitzer/IRS spectra for complete samples of 46 2Jy radio galaxies (0.05<z<0.7) and 19 3CRR FRII radio galaxies (z<0.1), and use the detection of polycyclic aromatic hydrocarbon (PAH) features to examine the incidence of contemporaneous star formation and radio-loud AGN activity. Our analysis reveals PAH features in only a minority (30%) of the objects with good IRS spectra. Using the wealth of complementary data available for the 2Jy and 3CRR samples we make detailed comparisons between a range of star formation diagnostics: optical continuum spectroscopy, mid- to far-IR (MFIR) color, far-IR excess and PAH detection. There is good agreement between the various diagnostic techniques: most candidates identified to have star formation activity on the basis of PAH detection are also identified using at least two of the other techniques. We find that only 35% of the combined 2Jy and 3CRR sample show evidence for recent star formation activity (RSFA) at optical and/or MFIR wavelengths. This result argues strongly against the idea of a close link between starburst and powerful radio-loud AGN activity, reinforcing the view that, although a large fraction of powerful radio galaxies may be triggered in galaxy interactions, only a minority are triggered at the peaks of star formation activity in major, gas-rich mergers. However, we find that compact radio sources (D < 15 kpc) show a significantly higher incidence of RSFA (>75%) than their more extended counterparts (=15 -- 25%). We discuss this result in the context of a possible bias towards the selection of compact radio sources triggered in gas-rich environments.

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

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.

F00183-7111 is one of the most extreme Ultra-Luminous Infrared Galaxies known. Here we present a VLBI image which shows that F00183-7111 is powered by a combination of a radio-loud Active Galactic Nucleus surrounded by vigorous starburst activity. Although already radio-loud, the quasar jets are only 1.7 kpc long, boring through the dense gas and starburst activity that confine them. We appear to be witnessing this remarkable source in the brief transition period between merging starburst and radio-loud "quasar-mode" accretion.

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.

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.

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 ({\tau}_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.

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 {\mu}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.

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

We present the first principal component analysis (PCA) applied to a sample of 119 Spitzer Infrared Spectrograph (IRS) spectra of local ultraluminous infrared galaxies (ULIRGs) at z<0.35. The purpose of this study is to objectively and uniquely characterise the local ULIRG population using all information contained in the observed spectra. We have derived the first three principal components (PCs) from the covariance matrix of our dataset which account for over 90% of the variance. The first PC is characterised by dust temperatures and the geometry of the mix of source and dust. The second PC is a pure star formation component. The third PC represents an anti-correlation between star formation activity and a rising AGN. Using the first three PCs, we are able to accurately reconstruct most of the spectra in our sample. Our work shows that there are several factors that are important in characterising the ULIRG population, dust temperature, geometry, star formation intensity, AGN contribution, etc. We also make comparison between PCA and other diagnostics such as ratio of the 6.2 microns PAH emission feature to the 9.7 micron silicate absorption depth and other observables such as optical spectral type.

The "Extreme starbursts in the local universe" workshop was held at the Insituto de Astrofisica de Andalucia in Granada, Spain on 21-25 June 2010. Bearing in mind the advent of a new generation of facilities such as JWST, Herschel, ALMA, eVLA and eMerlin, the aim of the workshop was to bring together observers and theorists to review the latest results. The purpose of the workshop was to address the following issues: what are the main modes of triggering extreme starbursts in the local Universe? How efficiently are stars formed in extreme starbursts? What are the star formation histories of local starburst galaxies? How well do the theoretical simulations model the observations? What can we learn about starbursts in the distant Universe through studies of their local counterparts? How important is the role of extreme starbursts in the hierarchical assembly of galaxies? How are extreme starbursts related to the triggering of AGN in the nuclei of galaxies? Overall, 41 talks and 4 posters with their corresponding 10 minutes short talks were presented during the workshop. In addition, the workshop was designed with emphasis on discussions, and therefore, there were 6 discussion sessions of up to one hour during the workshop. Here is presented a summary of the purposes of the workshop as well as a compilation of the abstracts corresponding to each of the presentations. The summary and conclusions of the workshop along with a description of the future prospects by Sylvain Veilleux can be found in the last section of this document. A photo of the assistants is included.

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).

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.

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.

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$\mu$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$\mu$m/[SIII]33.48$\mu$m as a function of [NeIII]15.56$\mu$m/[NeII]12.81$\mu$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$\alpha$) vs. log([OIII]/H$\beta$) 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$\mu$m/[NeII]12.81$\mu$m vs. [SIV]10.51$\mu$m/[SIII]18.71$\mu$m is not as efficient in separating AGNs from BCDs and pure starbursts. (abridged)

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.

We present the first detections of the 3um water ice and 3.4um amorphous hydrocarbon (HAC) absorption features in z~2 ULIRGs. These are based on deep rest-frame 2-8um Spitzer IRS spectra of 11 sources selected for their appreciable silicate absorption. The HAC-to-silicate ratio for our z~2 sources is typically higher by a factor of 2-5 than that observed in the Milky Way. This HAC `excess' suggests compact nuclei with steep temperature gradients as opposed to predominantly host obscuration. Beside the above molecular absorption features, we detect the 3.3um PAH emission feature in one of our sources with three more individual spectra showing evidence for it. Stacking analysis suggests that water ice, hydrocarbons, and PAH are likely present in the bulk of this sample even when not individually detected. The most unexpected result of our study is the lack of clear detections of the 4.67um CO gas absorption feature. Only three of the sources show tentative signs of this feature and at significantly lower levels than has been observed in local ULIRGs. Overall, we find that the closest local analogs to our sources, in terms of 3-4um color, HAC-to-silicate and ice-to-silicate ratios, as well as low PAH equivalent widths are sources dominated by deeply obscured nuclei. Such sources form only a small fraction of ULIRGs locally and are commonly believed to be dominated by buried AGN. Our sample suggests that, in absolute number, such buried AGN are at least an order of magnitude more common at z~2 than today. The presence of PAH suggests that significant levels of star-formation are present even if the obscured AGN typically dominate the power budget.

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.

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.

We apply methods from Bayesian inferencing and graph theory to a dataset of 102 mid-infrared spectra, and archival data from the optical to the millimeter, to construct an evolutionary paradigm for z<0.4 infrared-luminous galaxies (ULIRGs). We propose that the ULIRG lifecycle consists of three phases. The first phase lasts from the initial encounter until approximately coalescence. It is characterized by homogeneous mid-IR spectral shapes, and IR emission mainly from star formation, with a contribution from an AGN in some cases. At the end of this phase, a ULIRG enters one of two evolutionary paths depending on the dynamics of the merger, the available quantities of gas, and the masses of the black holes in the progenitors. On one branch, the contributions from the starburst and the AGN to the total IR luminosity decline and increase respectively. The IR spectral shapes are heterogeneous, likely due to feedback from AGN-driven winds. Some objects go through a brief QSO phase at the end. On the other branch, the decline of the starburst relative to the AGN is less pronounced, and few or no objects go through a QSO phase. We show that the 11.2 micron PAH feature is a remarkably good diagnostic of evolutionary phase, and identify six ULIRGs that may be archetypes of key stages in this lifecycle.

We present preliminary results from the highest available signal-to-noise rest-frame 2-8um spectra of z~2 ULIRGs. Our 10 targets are selected for their deep silicate absorption features based on previous shallower IRS spectra. The goal of this follow-up program is: 1) allow for a more accurate analysis of inner/hot dust continuum, 2) detecting the 3.3um and 6.2um PAH features, and 3) detecting molecular absorption features such as due to water ice and hydrocarbons (HACs). We find that the 3.4um HAC absorption feature is observed in four sources, while the 3.05um water ice feature is observed in three of the sources. The HAC detectability is higher and ice detectability lower than expected from local ULIRGs, but consistent with a more AGN-dominated sample such as this one. Where ice is detected, the ice-to-silicate ratio is somewhat lower than many local ULIRGs implying on average thinner ice mantles. One source shows the, to our knowledge, highest redshift reported detection of the 3.3um PAH feature (along with a previously detected 6.2um feature). The strength of the 3.3um feature is as expected for a starburst-dominated ULIRG.

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.

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.

We present mid-infrared spectra of thirty two high redshift ultraluminous infrared galaxies, selected via the stellar photospheric feature at rest-frame 1.6um, and an observed-frame 24um flux of >500muJy. Nearly all the sample reside in a redshift range of <z>=1.71+/-0.15, and have rest-frame 1-1000um luminosities of 10^12.9 - 10^13.8 Lsun. Most of the spectra exhibit prominent polycyclic aromatic hydrocarbon emission features, and weak silicate absorption, consistent with a starburst origin for the IR emission. Our selection method appears to be a straightforward and efficient way of finding distant, IR-luminous, star-forming galaxies in narrow redshift ranges. There is however evidence that the mid-IR spectra of our sample differ systematically from those of local ULIRGs; our sample have comparable PAH equivalent widths but weaker apparent silicate absorption, and (possibly) enhanced PAH 6.2um/7.7um and 6.2um/11.2um flux ratios. Furthermore, the composite mid-IR spectrum of our sample is almost identical to that of local starbursts with IR luminosities of 10^10-10^11 Lsun rather than that of local ULIRGs. These differences are consistent with a reduced dust column, which can plausibly be obtained via some combination of (1) star formation that is extended over spatial scales of 1-4Kpc, and (2) star formation in unusually gas-rich regions.

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.

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.

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.

(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.

We present mid-IR spectral decomposition of a sample of 48 Spitzer-selected ULIRGs spanning z~1-3 and likely L_IR~10^12-10^13Lsun. Our study aims at quantifying the star-formation and AGN processes in these sources which recent results suggest have evolved strongly between the observed epoch and today. To do this, we study the mid-IR contribution of PAH emission, continuum, and extinction. About 3/4 of our sample are continuum- (i.e. AGN) dominated sources, but ~60% of these show PAH emission, suggesting the presence of star-formation activity. These sources have redder mid-IR colors than typical optically-selected quasars. About 25% of our sample have strong PAH emission, but none are likely to be pure starbursts as reflected in their relatively high 5um hot dust continua. However, their steep 30um-to-14um slopes suggest that star-formation might dominate the total infrared luminosity. Six of our z~2 sources have EW6.2>~0.3um and L_14um>~10^12Lsun (implying L_IR>~10^13Lsun). At these luminosities, such high EW6.2 ULIRGs do not exist in the local Universe. We find a median optical depth at 9.7um of <tau_9.7>=1.4. This is consistent with local IRAS-selected ULIRGs, but differs from early results on SCUBA-selected z~2 ULIRGs. Similar to local ULIRGs about 25% of our sample show extreme obscuration (tau_9.7>~3) suggesting buried nuclei. In general, we find that our sources are similar to local ULIRGs, but are an order of magnitude more luminous. It is not clear whether our z~2 ULIRGs are simply scaled-up versions of local ULIRGs, or subject to fundamentally different physical processes.

The study of the dominant population of high redshift IR-luminous galaxies (10^11 - 10^12 Lsun at 1<z<3), requires observation of sources at the ~0.1 mJy level in the mid-IR. We present the deepest spectra taken to date with the Infrared Spectrograph (IRS) on the Spitzer Space Telescope. We targeted two faint (f24~0.15 mJy) sources in the Southern GOODS field at z=1.09 and z=2.69. Spectra of the lower redshift target were taken in the observed-frame 8--21 micron range, while the spectrum of the higher redshift target covered 21--37 microns. We also present the spectra of two secondary sources within the slit. We detect strong PAH emission in all four targets, and compare the spectra to those of local galaxies observed by the IRS. The z=1.09 source appears to be a typical, star-formation dominated IR-luminous galaxy, while the z=2.69 source is a composite source with strong star formation and a prominent AGN. The IRAC colors of this source show no evidence of rest-frame near-infrared stellar photospheric emission. We demonstrate that an AGN which contributes only a small (~10%) fraction of the bolometric luminosity can produce enough hot dust emission to overwhelm the near-IR photospheric emission from stars. Such sources would be excluded from photometric surveys which rely on the near-IR bump to identify starbursts, leading to an underestimate of the star formation rate density.

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.

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.

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.

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.

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.

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.

We analyze extensive data on NGC1377 from the Spitzer Infrared Nearby Galaxies Survey (SINGS). Within the category of nascent starbursts, that we previously selected by their infrared to radio continuum ratios in large excess of the average and their hot dust, NGC1377 has the largest infrared excess yet measured. Optical imaging reveals a morphological distortion suggestive of a recent accretion event. Infrared spectroscopy reveals a compact and opaque source dominated by a hot, self-absorbed continuum (tau ~ 20 in the 10 micron silicate band). We provide physical evidence against non-stellar activity being the heating source. HII regions are detected through the single [NeII] line, probing <1% of the ionizing radiation. Not only is the optical depth very high, but >85% of ionizing photons are suppressed by dust. The only other detected emission features are molecular hydrogen lines, arguably excited mainly by shocks, besides photodissociation regions, and weak aromatic bands. The new observations support our interpretation in terms of an extremely young starburst (<1 Myr). More generally, galaxies deficient in radio synchrotron are likely observed within a few Myr of the onset of a starburst and after a long quiescence, prior to the replenishment of the ISM with cosmic rays. The similar infrared-radio properties of NGC1377 and some infrared-luminous galaxies suggest that NGC1377 constitutes an archetype to better understand starburst evolution. Although rare locally because observed in a brief evolutionary stage, nascent starbursts may represent a non-negligible fraction of merger-induced starbursts dominating deep infrared counts. Since they differ dramatically from usual starburst templates, they bear important consequences for the interpretation of deep surveys.

We present the first MIR spectrum of the z=2.2856 ultraluminous, infrared galaxy FSC 10214+4724, obtained with the Infrared Spectrograph onboard the Spitzer Space Telescope. The spectrum spans a rest wavelength range of 2.3-11.5 microns, covering a number of key diagnostic emission and absorption features. The most prominent feature in the IRS spectrum is the silicate emission at rest-frame 10 microns. We also detect an unresolved emission line at a rest wavelength of 7.65 microns which we identify with [NeVI], and a slightly resolved feature at 5.6 microns identified as a blend of [Mg VII] and [Mg V]. There are no strong PAH emission features in the FSC 10214+4724 spectrum. We place a limit of 0.1 micron on the equivalent width of 6.2 micron PAH emission but see no evidence of a corresponding 7.7 micron feature. Semi-empirical fits to the spectral energy distribution suggest about 45% of the bolometric luminosity arises from cold 50 K dust, half arises from warm (190 K) dust, and the remainder, 5%, originates from hot (640 K) dust. The hot dust is required to fit the blue end of the steep MIR spectrum. The combination of a red continuum, strong silicate emission, little or no PAH emission, and no silicate absorption, makes FSC 10214+4724 unlike most other ULIRGs or AGN observed thus far with IRS. These apparently contradictory properties may be explained by an AGN which is highly magnified by the lens, masking a (dominant) overlying starburst with unusually weak PAH emission.

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.

We report on our first results from a mid-infrared spectroscopic study of ISM features in a sample of deeply obscured ULIRG nuclei using the InfraRed Spectrograph (IRS) on the Spitzer Space Telescope. The spectra are extremely rich and complex, revealing absorption features of both amorphous and crystalline silicates, aliphatic hydrocarbons, water ice and gas phase bands of hot CO and warm C_2H_2, HCN and CO_2. PAH emission bands were found to be generally weak and in some cases absent. The features are probing a dense and warm environment in which crystalline silicates and water ice are able to survive but volatile ices, commonly detected in Galactic dense molecular clouds, cannot. If powered largely by star formation, the stellar density and conditions of the gas and dust have to be extreme not to give rise to the commonly detected emission features associated with starburst.

We present mid-infrared spectra of the nearby, Ultraluminous Infrared Galaxy NGC 6240 taken with the IRS on Spitzer. The spectrum of NGC 6240 is dominated by strong fine-structure lines, rotational H2 lines, and PAH emission features. The H2 line fluxes suggest molecular gas at a variety of temperatures. A simple two-temperature fit to the S(0) through S(7) lines implies a mass of 6.7E6 solar masses at T=957K and 1.6E9 solar masses at T=164K, or about 15% of the total molecular gas mass in this system. Noteably, we have detected the [NeV] 14.3 micron emission line, with a flux of 5E-14 erg/cm^2/s, providing the first direct detection of the buried AGN in the mid-infrared. Modelling of the total SED from near to far-infrared wavelengths requires the presence of a hot dust (T=700K) component, which we also associate with the buried AGN. The small [NeV]/[NeII] and [NeV]/IR flux ratios, the relative fraction of hot dust emission, as well as the large 6.2 micron PAH EQW, are all consistent with an apparent AGN contribution of only 3-5% to the bolometric luminosity. However, correcting the measured [NeV] flux by the extinction implied by the silicate optical depth and our SED fitting, suggests an intrinsic fractional AGN contribution to the bolometric luminosity of about 20-24% in NGC 6240, which lies within the range implied by fits to the hard X-ray spectrum.

We present a {\em Spitzer}-based mid-infrared study of a large sample of Blue Compact Dwarf galaxies (BCD) using the Infrared Spectrograph (IRS), including the first mid-IR spectrum of IZw18, the archetype for the BCD class and among the most metal poor galaxies known. We show the spectra of Polycyclic Aromatic Hydrocarbon (PAH) emission in low-metallicity environment. We find that the equivalent widths (EW) of PAHs at 6.2, 7.7, 8.6 and 11.2 $\mu$m are generally weaker in BCDs than in typical starburst galaxies and that the fine structure line ratio, [NeIII]/[NeII], has a weak anti-correlation with the PAH EW. A much stronger anti-correlation is shown between the PAH EW and the product of the [NeIII]/[NeII] ratio and the UV luminosity density divided by the metallicity. We conclude that PAH EW in metal-poor high-excitation environments is determined by a combination of PAH formation and destruction effects.

Silicates are an important component of interstellar dust and the structure of these grains -- amorphous versus crystalline -- is sensitive to the local physical conditions. We have studied the infrared spectra of a sample of ultra-luminous infrared galaxies. Here, we report the discovery of weak, narrow absorption features at 11, 16, 19, 23, and 28 microns, characteristic of crystalline silicates, superimposed on the broad absorption bands at 10 and 18 microns due to amorphous silicates in a subset of this sample. These features betray the presence of forsterite (Mg_2SiO_4), the magnesium-rich end member of the olivines. Previously, crystalline silicates have only been observed in circumstellar environments. The derived fraction of forsterite to amorphous silicates is typically 0.1 in these ULIRGs. This is much larger than the upper limit for this ratio in the interstellar medium of the Milky Way, 0.01. These results suggest that the timescale for injection of crystalline silicates into the ISM is short in a merger-driven starburst environment (e.g., as compared to the total time to dissipate the gas), pointing towards massive stars as a prominent source of crystalline silicates. Furthermore, amorphization due to cosmic rays, which is thought to be of prime importance for the local ISM, lags in vigorous starburst environments.

We report the spectroscopic detection of silicate emission at 10 and 18 microns in five PG quasars, the first detection of these two features in galaxies outside the Local Group. This finding is consistent with the unification model for Active Galactic Nuclei (AGNs), which predicts that an AGN torus seen pole-on should show a silicate emission feature in the mid-infrared. The strengths of the detected silicate emission features range from 0.12 to 1.25 times the continuum at 10 microns and from 0.20 to 0.79 times the continuum at 18 microns. The silicate grain temperatures inferred from the ratio of 18-to-10 micron silicate features under the assumption of optically thin emission range from 140 to 220K.

Prominent Fe K line emission is detected at around 6.7 keV in the XMM-Newton spectrum of the ultraluminous infrared galaxy Arp220. The continuum emission in the 2.5-10 keV band is flat and a few other spectral features are suggested. The large EW of the Fe K line poses a problem with interpreting the hard X-ray emission as integrated X-ray binary emission. A thermal emission spectrum with a temperature of kT~7 keV modified by some absorption, which would be expected from an internally shocked hot bubble in a starburst region, can describe the spectrum. An ensemble of radio supernovae in a dense environment, as suggested from VLBI imaging, could be another possibility, if frequent occurence of such powerful supernovae is sustained. However, the apparent lack of X-ray binary emission does not match the high supernova rate required by both interpretations. Highly photoionized, low-density gas illuminated by a hidden Compton-thick AGN is a possible alternative, which can be tested by presence/absence of radiative recombination continua in better quality data expected from a forthcoming observation.

In order to constrain the properties of dense and warm gas around active galactic nuclei, we have searched Infrared Space Observatory spectra of local active galactic nuclei for the signature of the 4.7micron fundamental ro-vibrational band of carbon monoxide. Low resolution spectra of 31 AGN put upper limits on the presence of wide absorption bands corresponding to absorption by large columns of warm and dense gas against the nuclear dust continuum. High resolution (R~2500) spectra of NGC 1068 detect no significant absorption or emission in individual lines, to a 3sigma limit of 7% of the continuum. The limits set on CO absorption in local AGN are much lower than the recent Spitzer Space Telescope detection of strong CO absorption by dense and warm gas in the obscured ultraluminous infrared galaxy IRAS F00183-7111, despite evidence for dense material on parsec scales near an AGN in both types of objects. This suggests that such deep absorptions are not intimately related to the obscuring `torus' material invoked in local AGN, but rather are a signature of the peculiar conditions in the circumnuclear region of highly obscured infrared galaxies like IRAS F00183-7111. They may reflect fully covered rather than torus geometries.

We present observations of the starburst galaxy NGC 7714 with the Infrared Spectrograph IRS on board the Spitzer Space Telescope. The spectra yield a wealth of ionic and molecular features that allow a detailed characterization of its properties. NGC 7714 has an HII region-like spectrum with strong PAH emission features. We find no evidence for an obscured active galactic nucleus, and with [NeIII]/[NeII]~0.73, NGC7714 lies near the upper end of normal-metallicity starburst galaxies. With very little slicate absorption and a temperature of the hottest dust component of 340K, NGC 7714 is the perfect template for a young, unobscured starburst

We report the detection of strong absorption and weak emission features in the 4--27 micron Spitzer-IRS spectrum of the distant ultraluminous infrared galaxy (ULIRG) IRAS F00183--7111 (z=0.327). The absorption features of CO2 and CO gas, water ice, hydrocarbons and silicates are indicative of a strongly obscured (A[9.6]>=5.4; A[V]>=90) and complex line of sight through both hot diffuse ISM and shielded cold molecular clouds towards the nuclear power source. From the profile of the 4.67 micron CO fundamental vibration mode we deduce that the absorbing gas is dense (n~10^6 cm^-3) and warm (720 K) and has a CO column density of ~10^19.5 cm^-2, equivalent to N[H]~10^23.5 cm^-2. The high temperature and density, as well as the small infered size (<0.03pc), locates this absorbing gas close to the power source of this region. Weak emission features of molecular hydrogen, PAHs and Ne+, likely associated with star formation, are detected against the 9.7 micron silicate feature, indicating an origin away from the absorbing region. Based on the 11.2 micron PAH flux, we estimate the star formation component to be responsible for up to 30% of the IR luminosity of the system. While our mid-infrared spectrum shows no tell-tale signs of AGN activity, the similarities to the mid-infrared spectra of deeply obscured sources (e.g. NGC4418) and AGN hot dust (e.g. NGC1068), as well as evidence from other wavelength regions, suggest that the power source hiding behind the optically thick dust screen may well be a buried AGN.

We present our early results of the mapping of the nucleus of the starburst galaxy NGC 253 and its immediate surroundings using the Infrared Spectrograph onboard the Spitzer Space Telescope. The map is centered on the nucleus of the galaxy and spans the inner 800 x 688 pc2. The spatial distribution of the \neiii line at 15.55 microns and the polycyclic aromatic hydrocarbons feature at 11.3 microns peaks at the center while the purely rotationnal transition of molecular hydrogen at 17.03 microns is strong over several slit positions. We perform a brief investigation of the implications of these measurement on the properties of the star formation in this region using theories developped to explain the deficiency of massive stars in starbursts.

We present spectra taken with the Infrared Spectrograph on Spitzer covering the 5-38micron region of three Ultraluminous Infrared Galaxies (ULIRGs): Mrk 1014 (z=0.163), and Mrk 463 (z=0.051), and UGC 5101 (z=0.039). The continua of UGC 5101 and Mrk 463 show strong silicate absorption suggesting significant optical depths to the nuclei at 10microns. UGC 5101 also shows the clear presence of water ice in absorption. PAH emission features are seen in both Mrk 1014 and UGC 5101, including the 16.4micron line in UGC 5101. The fine structure lines are consistent with dominant AGN power sources in both Mrk 1014 and Mrk 463. In UGC 5101 we detect the [NeV] 14.3micron emission line providing the first direct evidence for a buried AGN in the mid-infrared. The detection of the 9.66micron and 17.03micron H$_{2}$ emission lines in both UGC 5101 and Mrk 463 suggest that the warm molecular gas accounts for 22% and 48% of the total molecular gas masses in these galaxies.

SBS0335-052 is a blue compact dwarf galaxy (BCD) with one of the lowest known metallicities, Z$\sim$Z$_{\sun}$/41, making it a local example of how primordial starburst galaxies and their precursors might appear. A spectrum obtained with the Infrared Spectrograph (IRS) on the Spitzer Space Telescope clearly shows silicate absorption features, emission lines of [SIV] and [NeIII], and puts strong upper limits on the PAH emission features. The observed low resolution spectrum (R~90) extends from 5.3 to 35microns and peaks at ~28microns. The spectrum is compared to IRS observations of the prototypical starburst nucleus NGC7714. SBS0335-052 is quite unlike normal starburst galaxies, which show strong PAH bands, low ionization emission lines, and a continuum peak near 80microns. The continuum difference for $\lambda >30 \mu$m implies a substantial reduction in the mass of cold dust. If the spectrum of this very low metallicity galaxy is representative of star forming galaxies at higher redshifts, it may be difficult to distinguish them from AGNs which also show relatively featureless flat spectra in the mid-IR.

IR emission bands at 3.3, 6.2, 7.7, 8.6 and 11.3 um are generally attributed to IR fluorescence from (mainly) FUV pumped PAHs. As such, they trace the FUV stellar flux and are a measure of star formation. We examined the IR spectral characteristics of Galactic star forming regions, normal and starburst galaxies, AGNs and ULIRGs. The goal is to analyze if PAH bands are a good qualitative and/or quantitative tracer of star formation and hence the application of PAH bands as a diagnostic in order to identify the dominant processes contributing to the IR emission from Seyfert's and ULIRGs. We develop a MIR/FIR diagnostic and compare it to known diagnostics, with these also applied to the Galactic sample. This diagnostic is based on the FIR normalized 6.2 um PAH flux and the FIR normalized 6.2 um continuum flux. The Galactic sources form a sequence spanning a range of 3 orders of magnitude, from embedded compact HII regions to exposed PDRs and the (D)ISM. The variation in the 6.2 um PAH/continuum ratio is relative small. Normal and starburst galaxies ressemble exposed PDRs. While Seyfert-2's coincide with the starburst trend, Seyfert-1's are displaced by at least a factor 10 in 6.2 um continuum flux. ULIRGs show a diverse spectral appearance (AGN hot dust continuum, starburst-like or strong dust obscuration in the nucleus). ULIRGs also seems to have more prominent FIR emission than either starburst galaxies or AGNs. We discuss the observed variation in the Galactic sample in view of the evolutionary state and the PAH/dust abundance and the use of PAHs as quantitative tracers of star formation activity. We find that PAHs may be better suited as a tracer of B stars, which dominate the Galactic stellar energy budget, than as a tracer of massive star formation (O stars).

We use mid-infrared spectral decomposition to separate the 6micron mid-infrared AGN continuum from the host emission in the ISO low resolution spectra of 71 active galaxies and compare the results to observed and intrinsic 2-10keV hard X-ray fluxes from the literature. We find a correlation between mid-infrared luminosity and absorption corrected hard X-ray luminosity, but the scatter is about an order of magnitude, significantly larger than previously found with smaller statistics. Main contributors to this scatter are likely variations in the geometry of absorbing dust, and AGN variability in combination with non-simultaneous observations. There is no significant difference between Type 1 and Type 2 objects in the average ratio of mid-infrared and hard X-ray emission, a result which is not consistent with the most simple version of a unified scheme in which an optically and geometrically thick torus dominates the mid-infrared AGN continuum. Most probably, significant non-torus contributions to the AGN mid-IR continuum are masking the expected difference between the two types of AGN.

We present ISO PHOT-S spectra of four galaxies known or suspected to host a central AGN, selected from the initial Iras/Rosat sample of Boller et al. (1992). Two of them had no obvious Seyfert features in their previous optical spectra: IRAS 14201+2956, and IRAS 21582+1018. The latter was bright enough to also allow SWS observations around selected neon lines, to establish its excitation. While both PHOT-S spectra are characteristic of starburst-dominated galaxies, the neon line ratios in IRAS 21582+1018 indicate the presence of a hard excitation source. New, high-resolution, optical spectra show only a weak, broad component around Halpha, classifying now these two objects as Sey 1.9 galaxies. The two other galaxies observed are the NLS1 galaxies Mrk 359 and Mrk 1388. Their ISO spectra however do not reveal the typical, strong PAH features found in the starburst galaxies and are more like those of standard Seyferts. These results show therefore that, although IR observations were expected to be able to always reveal the presence of an active nucleus by piercing through the central obscuration, the result may be ambiguous: the broad band IR energy distribution can still be dominated by starburts located in a circumnuclear region, and the AGN appear only in specific observations (high-excitation lines in the IR, or optical spectra with better quality than classification spectra). The obscuration needs however to be patchy rather than complete, to explain the detection of the high-excitation lines or broad Balmer wings. Only high-energy observations can then establish the strength of the central AGN and the amount of extinction with certainty.

We have re-analyzed the 6-12 micron ISO spectrum of the ultra-luminous infrared galaxy Arp220 with the conclusion that it is not consistent with that of a scaled up version of a typical starburst. Instead, both template fitting with spectra of the galaxies NGC4418 and M83 and with dust models suggest that it is best represented by combinations of a typical starburst component, exhibiting PAH emission features, and a heavily absorbed dust continuum which contributes ~40% of the 6-12 micron flux and likely dominates the luminosity. Of particular significance relative to previous studies of Arp220 is the fact that the emission feature at 7.7 micron comprises both PAH emission and a broader component resulting from ice and silicate absorption against a heavily absorbed continuum. Extinction to the PAH emitting source, however, appears to be relatively low. We tentatively associate the PAH emitting and heavily dust/ice absorbed components with the diffuse emission region and the two compact nuclei respectively identified by Soifer et al. (2002) in their higher spatial resolution 10 micron study. Both the similarity of the absorbed continuum with that of the embedded Galactic protostars and results of the dust models imply that the embedded source(s) in Arp220 could be powered by, albeit extremely dense, starburst activity. Due to the high extinction, it is not possible with the available data to exclude that AGN(s) also contribute some or all of the observed luminosity. In this case, however, the upper limit measured for its hard X-ray emission would require Arp220 to be the most highly obscured AGN known.

We investigate the mid-infrared radiation of galaxies that are powered by a starburst or by an AGN. For this end, we compare the spectra obtained at different spatial scales in a sample of infrared bright galaxies. ISO observations which include emission of the nucleus as well as most of the host galaxy are compared with TIMMI2 spectra of the nuclear region. We find that ISO spectra are generally dominated by strong PAH bands. However, this is no longer true when inspecting the mid-infrared emission of the pure nucleus. Here PAH emission is detected in starbursts whereas it is significantly reduced or completely absent in AGNs. A physical explanation of these new observational results is presented by examining the temperature fluctuation of a PAH after interaction with a photon. It turns out that the hardness of the radiation field is a key parameter for quantifying the photo-destruction of small grains. Our theoretical study predicts PAH evaporation in soft X-ray environments. Radiative transfer calculations of clumpy starbursts and AGN corroborate the observational fact that PAH emission is connected to starburst activity whereas PAHs are destroyed near an AGN. The radiative transfer models predict for starbursts a much larger mid-infrared size than for AGN. This is confirmed by our TIMMI2 acquisition images: We find that the mid-infrared emission of Seyferts is dominated by a compact core while most of the starbursts are spatially resolved.

We present Infrared Space Observatory mid- and far-infrared spectroscopy of the merging galaxy NGC 6240, an object presenting many aspects of importance for the role of star formation and AGN activity in [ultra]luminous infrared galaxies. The mid-infrared spectrum shows starburst indicators in the form of low excitation fine-structure line emission and aromatic `PAH' features. A strong high excitation [OIV] line is observed which most likely originates in the Narrow Line Region of an optically obscured AGN. NGC 6240 shows extremely powerful emission in the pure rotational lines of molecular hydrogen. We argue that this emission is mainly due to shocks in its turbulent central gas component and its starburst superwind. The total shock cooling in infrared emission lines accounts for ~0.6% of the bolometric luminosity, mainly through rotational H_2 emission and the [OI] 63micron line. We analyse several ways of estimating the luminosities of the starburst and the AGN in NGC 6240 and suggest that the contributions to its bolometric luminosity are most likely in the range 50-75% starburst and 25-50% AGN.

The composition of ice grains provides an important tool for the study of the molecular environment of star forming regions. Using ISAAC at the VLT to obtain spectra around 4.65 microns we have detected for the first time `XCN' and CO ice in an extragalactic environment: the nuclear region of the nearby dusty starburst/AGN galaxy NGC4945. The profile of the solid CO band reveals the importance of thermal processing of the ice while the prominence of the XCN band attests to the importance of energetic processing of the ice by FUV radiation and/or energetic particles. In analogy to the processing of ices by embedded protostars in our Galaxy, we attribute the processing of the ices in the center of NGC4945 to ongoing massive star formation. Our M-band spectrum also shows strong HI Pfund-beta and H2 0-0 S(9) line emission and gas phase CO absorption lines. The HI, H2, PAH, gas phase CO and the ices seem to be embedded in a rotating molecular disk which is undergoing vigorous star formation. Recently, strong OCN- absorption has been detected in the spectrum of the Galactic center star GC:IRS19. The most likely environment for the OCN- absorption is the strongly UV-exposed GC molecular ring. The presence of processed ice in the center of NGC4945 and our Galactic center leads us to believe that processed ice may be a common characteristic of dense molecular material in star forming galactic nuclei.

Mid infrared spectra provide a powerful probe of the conditions in dusty galactic nuclei. They variously contain emission features associated with star forming regions and absorptions by circumnuclear silicate dust plus ices in cold molecular cloud material. Here we report the detection of 6-8um water ice absorption in 18 galaxies observed by ISO. While the mid-IR spectra of some of these galaxies show a strong resemblance to the heavily absorbed spectrum of NGC 4418, other galaxies in this sample also show weak to strong PAH emission. The 18 ice galaxies are part of a sample of 103 galaxies with good S/N mid-IR ISO spectra. Based on our sample we find that ice is present in most of the ULIRGs, whereas it is weak or absent in the large majority of Seyferts and starburst galaxies. This result is consistent with the presence of larger quantities of molecular material in ULIRGs as opposed to other galaxy types. Like NGC 4418, several of our ice galaxy spectra show a maximum near 8um that is not or only partly due to PAH emission. While this affects only a small part of the galaxy population studied by ISO, it stresses the need for high S/N data and refined diagnostic methods, to properly discriminate spectra dominated by PAH emission and spectra dominated by heavy obscuration. The spectral variation from PAH emission to absorbed continuum emission near 8um shows strong similarities with Galactic star forming clouds. This leads us to believe that our classification of ice galaxy spectra might reflect an evolutionary sequence from strongly obscured beginnings of star formation (and AGN activity) to a less enshrouded stage of advanced star formation (and AGN activity), as the PAHs get stronger and the broad 8um feature weakens.

We present low resolution mid infrared spectra of 16 ultraluminous infrared galaxies (ULIRGs) obtained with the CVF spectroscopy mode of ISOCAM on board the Infrared Space Observatory ISO. Our sample completes previous ISO spectroscopy of ultra- and hyperluminous infrared galaxies towards higher luminosities. The combined samples cover an infrared luminosity range of \~10^{12 -13.1} Lo. For about half of the high luminosity ULIRGs studied here, strong aromatic emission bands suggest starburst dominance. Other spectra are dominated by a strong AGN-related continuum with weak superposed emission features of uncertain nature. An improved method to quantitatively characterize the relative contribution of star formation and AGN activity to the mid-infrared emission of ULIRGs is presented. As dominant source of the bolometric luminosity, starbursts prevail at the lower end and AGNs at the higher end of this range. The transition between mostly starburst and mostly AGN powered occurs at ~10^{12.4}$ to $10^{12.5} Lo, and individual luminous starbursts are found up to ~10^{12.65} Lo.

We report the detection of absorption features in the 6-8 micron region superimposed on a featureless mid-infrared continuum in NGC 4418. For several of these features this is the first detection in an external galaxy. We compare the absorption spectrum of NGC 4418 to that of embedded massive protostars and the Galactic centre, and attribute the absorption features to ice grains and to hydrogenated amorphous carbon grains. From the depth of the ice features, the powerful central source responsible for the mid-infrared emission must be deeply enshrouded. Since this emission is warm and originates in a compact region, an AGN must be hiding in the nucleus of NGC 4418.

We have observed the central region of the nearby starburst galaxy NGC4945 with the mid-infrared spectrometers SWS and PHT-S aboard ISO. We do not find any evidence for the existence of the powerful AGN, inferred from hard X-ray observations. The upper limits on our AGN tracers [NeV] 14.32um & 24.32um and [NeVI] 7.65um imply an A(V)>160 towards the NLR, assuming the NLR to be of equal strength as in the Circinus galaxy. Other possibilities are discussed. The starburst excitation indicators [NeIII]15.56um/[NeII]12.81um and L(bol)/L(lyc) suggest that the starburst in the central region is at least 5x10^6 yrs old, and that it accounts for at least half of the nuclear bolometric luminosity. The starburst might well power the entire bolometric luminosity, but the available constraints are also consistent with an up to 50% contribution of the embedded AGN. With PHT-S, at a resolution of approximately 90, we detect strong absorption features of water ice, and, for the first time in an external galaxy, of CO_2 and CO. The same PHT-S spectrum also reveals strong emission from the family of PAH features. Finally, we have observed and detected several pure rotational and ro-vibrational H_2 lines, two of which, the (0-0) S(0) & S(1) lines, allow us to determine the excitation temperature (160K) and warm H_2 mass (2.4x10^7 Msun). The low excitation temperature shows Orion-like shocks not to be representative for the entire emission of the central region of the galaxy and fairly normal PDRs to be perhaps more typical.

We present new ISO-SWS data for a sample of 27 starburst galaxies, and with these data examine the issues of formation and evolution of the most massive stars in starburst galaxies. Using starburst models which incorporate time evolution, new stellar atmosphere models for massive stars, and a starburst model geometry derived from observations of the prototypical starburst M82, we model the integrated mid-infrared line ratio [NeIII](15.6 microns)/[NeII](12.8 microns). This line ratio is sensitive to the hardness of the stellar energy distribution and therefore to the most massive stars present. We conclude from our models, with consideration of recent determinations of the stellar census in local, high-mass star forming regions, that the [NeIII]/[NeII] ratios we measure are consistent with the formation of massive (~50-100 solar mass) stars in most starbursts. In this framework, the low nebular excitation inferred from the measured line ratios can be attributed to aging effects. By including estimates of the ratio of infrared-to-Lyman continuum luminosity for the galaxies in our sample, we further find that most starbursts are relatively short-lived (1-10 million years), only a few O-star lifetimes. We discuss a possible cause of such short events: the effectiveness of stellar winds and supernovae in destroying the starburst environment.

We present a low resolution mid-infrared spectroscopic survey of an unbiased sample of 62 Ultraluminous Infrared Galaxies (ULIRGs) using ISOPHOT-S on board ISO. For comparison we also present ISOPHOT-S spectra for 23 AGN and 15 starburst and normal galaxies. The line-to-continuum ratio of the 7.7 $\mu$m PAH emission feature is used as a discriminator between starburst and AGN activity in ULIRGs. We find that the majority of ULIRGs are predominantly powered by starbursts. The ratio L$_{PAH}/$L$_{IR}$ for starburst dominated ULIRGs is similar to the one found for template starbursts.The shapes of the PAH features are sometimes unusual.Extinction has a noticeable effect on the PAH spectra of ULIRG starbursts. We have obtained high resolution near-infrared imaging for the Southern ISOPHOT-S ULIRGs to investigate their evolution stage. The majority (68%) of the ULIRGs imaged are double systems and all show distorted morphologies. Using separations for our work and from the literature, we have examined whether ULIRGs that are advanced mergers are more AGN-like. We have found no such evidence contrary to what is postulated by the classical evolutionary scenario. No correlation is found between the stage of merger in ULIRGs and their infrared luminosity.We also find that the total mass of interstellar gas, as estimated from the CO (1-0) luminosity, does not decrease with decreasing merger separation. When both an AGN and a starburst occur concurrently in ULIRGs, we find that the starburst dominates the luminosity output. We propose that the available gas reservoir and the individual structure of the interacting galaxies plays a major role in the evolution of the system.

We report the first results of a low resolution mid-infrared spectroscopic survey of an unbiased, far-infrared selected sample of 60 ultraluminous infrared galaxies, using ISOPHOT-S on board ISO. We use the ratio of the 7.7um `PAH' emission feature to the local continuum as a discriminator between starburst and AGN activity. About 80% of all the ULIRGs are found to be predominantly powered by star formation but the fraction of AGN powered objects increases with luminosity. Observed ratios of the PAH features in ULIRGs differ slightly from those in lower luminosity starbursts. This can be plausibly explained by the higher extinction and/or different physical conditions in the interstellar medium of ULIRGs. The PAH feature-to-continuum ratio is anticorrelated with the ratio of feature-free 5.9um continuum to the IRAS 60um continuum, confirming suggestions that strong mid-IR continuum is a prime AGN signature. The location of starburst-dominated ULIRGs in such a diagram is consistent with previous ISO-SWS spectroscopy which implies significant extinction even in the mid-infrared. We have searched for indications that ULIRGs which are advanced mergers might be more AGN-like, as postulated by the classical evolutionary scenario. No such trend has been found amongst those objects for which near infrared images are available to assess their likely merger status.

We have obtained ISOPHOT-S low resolution mid-infrared spectra of a sample of 60 Ultraluminous Infrared Galaxies (ULIRGs). We use the strength of the `PAH' mid-infrared features as a discriminator between starburst and AGN activity, and to probe for evolutionary effects. We focus on the fact that observed ratios of PAH features in ULIRGs differ slightly from those in lower luminosity starbursts. We suggest that such PAH ratio changes relate to the conditions in the interstellar medium in these galaxies, and in particular to extinction.

We report the detection of faint emission in the high-excitation [OIV] 25.90um line in a number of starburst galaxies, from observations obtained with the Short Wavelength Spectrometer (SWS) on board ISO. Further observations of M82 spatially resolve the [OIV] emitting region. Detection of this line in starbursts is surprising since it is not produced in measurable quantities in HII regions around hot main-sequence stars, the dominant energy source of starburst galaxies. We discuss various models for the formation of this line. [OIV] that is spatially resolved by ISO cannot originate in a weak AGN and must be due to very hot stars or ionizing shocks related to the starburst activity. For low-excitation starbursts like M82, shocks are the most plausible source of [OIV] emission.

We present an ISO SWS and ISOPHOT-S, mid-infrared spectroscopic survey of 15 ultra-luminous IRAS galaxies. We combine the survey results with a detailed case study, based on near-IR and mm imaging spectroscopy, of one of the sample galaxies (UGC 5101). We compare the near- and mid-IR characteristics of these ultra-luminous galaxies to ISO and literature data of thirty starburst and active galactic nuclei (AGN), template galaxies. We find that 1) 70-80% of the ultra-luminous IRAS galaxies in our sample are predominantly powered by recently formed massive stars. 20-30% are powered by a central AGN. These conclusions are based on a new infrared 'diagnostic diagram' involving the ratio of high to low excitation mid-IR emission lines on the one hand, and on the strength of the 7.7um PAH feature on the other hand. 2) at least half of the sources probably have simultaneously an active nucleus and starburst activity in a 1-2 kpc diameter circum-nuclear disk/ring. 3) the mid-infrared emitting regions are highly obscured. After correction for these extinctions, we estimate that the star forming regions in ULIRGs have ages between 10^7 and 10^8 years, similar to but somewhat larger than those found in lower luminosity starburst galaxies. 4) in the sample we have studied there is no obvious trend for the AGN component to dominate in the most compact, and thus most advanced mergers. Instead, at any given time during the merger evolution, the time dependent compression of the circum-nuclear interstellar gas, the accretion rate onto the central black hole and the associated radiation efficiency may determine whether star formation or AGN activity dominates the luminosity of the system.