Using $0.2^{\prime \prime}$ ($\sim3$ pc) ALMA images of vibrationally excited HC$_3$N emission (HC$_3$N$^*$) we reveal the presence of $8$ unresolved Super Hot Cores (SHCs) in the inner $160$ pc of NGC\,253. Our LTE and non-LTE modelling of the HC$_3$N$^*$ emission indicate that SHCs have dust temperatures of $200-375$ K, relatively high H$_2$ densities of $1-6\times 10^{6}$ cm$^{-3}$ and high IR luminosities of $0.1-1\times 10^8$ L$_\odot$. As expected from their short lived phase ($\sim 10^4$ yr), all SHCs are associated with young Super Star Clusters (SSCs). We use the ratio of luminosities form the SHCs (protostar phase) and from the free-free emission (ZAMS star phase), to establish the evolutionary stage of the SSCs. The youngest SSCs, with the larges ratios, have ages of a few $10^4$ yr (proto-SSCs) and the more evolved SSCs are likely between $10^5$ and $10^6$ yr (ZAMS-SSCs). The different evolutionary stages of the SSCs are also supported by the radiative feedback from the UV radiation as traced by the HNCO/CS ratio, with this ratio being systematically higher in the young proto-SSCs than in the older ZAMS-SSCs. We also estimate the SFR and the SFE of the SSCs. The trend found in the estimated SFE ($\sim40\%$ for proto-SSCs and $>85\%$ for ZAMS-SSCs) and in the gas mass reservoir available for star formation, one order of magnitude higher for proto-SSCs, suggests that star formation is still going on in proto-SSCs. We also find that the most evolved SSCs are located, in projection, closer to the center of the galaxy than the younger proto-SSCs, indicating an inside-out SSC formation scenario.

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.

Our current knowledge of star formation and accretion luminosity at high-redshift (z>3-4), as well as the possible connections between them, relies mostly on observations in the rest-frame ultraviolet (UV), which are strongly affected by dust obscuration. Due to the lack of sensitivity of past and current infrared (IR) instrumentation, so far it has not been possible to get a glimpse into the early phases of the dust-obscured Universe. Among the next generation of IR observatories, SPICA, observing in the 12-350 micron range, will be the only facility that can enable us to make the required leap forward in understanding the obscured star-formation rate and black-hole accretion rate densities (SFRD and BHARD, respectively) with respect to what Spitzer and Herschel achieved in the mid- and far-IR at z<3. In particular, SPICA will have the unique ability to trace the evolution of the obscured SFRD and BHARD over cosmic time, from the peak of their activity back to the reionisation epoch (i.e., 3<z<6-7), where its predecessors had severe limitations. Here we discuss the potential of both deep and shallow photometric surveys performed with the SPICA mid-IR instrument (SMI), enabled by the very low level of impact of dust obscuration in a band centred at 34 micron. These unique unbiased photometric surveys that SPICA will perform will be followed up by observations both with the SPICA spectrometers and with other facilities at shorter and longer wavelengths, with the aim to fully characterise the evolution of AGNs and star-forming galaxies after re-ionisation.

IR spectroscopy in the range 12-230 micron with the SPace IR telescope for Cosmology and Astrophysics (SPICA) will reveal the physical processes that govern the formation and evolution of galaxies and black holes through cosmic time, bridging the gap between the James Webb Space Telescope (JWST) and the new generation of Extremely Large Telescopes (ELTs) at shorter wavelengths and the Atacama Large Millimeter Array (ALMA) at longer wavelengths. SPICA, with its 2.5-m telescope actively-cooled to below 8K, will obtain the first spectroscopic determination, in the mid-IR rest-frame, of both the star-formation rate and black hole accretion rate histories of galaxies, reaching lookback times of 12 Gyr, for large statistically significant samples. Densities, temperatures, radiation fields and gas-phase metallicities will be measured in dust-obscured galaxies and active galactic nuclei (AGN), sampling a large range in mass and luminosity, from faint local dwarf galaxies to luminous quasars in the distant Universe. AGN and starburst feedback and feeding mechanisms in distant galaxies will be uncovered through detailed measurements of molecular and atomic line profiles. SPICA's large-area deep spectrophotometric surveys will provide mid-IR spectra and continuum fluxes for unbiased samples of tens of thousands of galaxies, out to redshifts of z~6. Furthermore, SPICA spectroscopy will uncover the most luminous galaxies in the first few hundred million years of the Universe, through their characteristic dust and molecular hydrogen features.

In Tombesi et al. (2015), we reported the first direct evidence for a quasar accretion disk wind driving a massive molecular outflow. The target was F11119+3257, an ultraluminous infrared galaxy (ULIRG) with unambiguous type-1 quasar optical broad emission lines. The energetics of the accretion disk wind and molecular outflow were found to be consistent with the predictions of quasar feedback models where the molecular outflow is driven by a hot energy-conserving bubble inflated by the inner quasar accretion disk wind. However, this conclusion was uncertain because the energetics were estimated from the optically thick OH 119 um transition profile observed with Herschel. Here, we independently confirm the presence of the molecular outflow in F11119+3257, based on the detection of broad wings in the CO(1-0) profile derived from ALMA observations. The broad CO(1-0) line emission appears to be spatially extended on a scale of at least ~7 kpc from the center. Mass outflow rate, momentum flux, and mechanical power of (80-200) R_7^{-1} M_sun/yr, (1.5-3.0) R_7^{-1} L_AGN/c, and (0.15-0.40)% R_7^{-1} L_AGN are inferred from these data, assuming a CO-to-H_2 conversion factor appropriate for a ULIRG (R_7 is the radius of the outflow normalized to 7 kpc and L_AGN is the AGN luminosity). These rates are time-averaged over a flow time scale of 7x10^6 yrs. They are similar to the OH-based rates time-averaged over a flow time scale of 4x10^5 yrs, but about a factor 4 smaller than the local ("instantaneous"; <10^5 yrs) OH-based estimates cited in Tombesi et al. The implications of these new results are discussed in the context of time-variable quasar-mode feedback and galaxy evolution. The need for an energy-conserving bubble to explain the molecular outflow is also re-examined.

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 the results from a systematic search for galactic-scale, molecular (OH 119 $\mu$m) outflows in a sample of 52 Local Volume ($d < 50$ Mpc) Burst Alert Telescope detected active galactic nuclei (BAT AGN) with \emph{Herschel}-PACS. We combine the results from our analysis of the BAT AGN with the published \emph{Herschel}/PACS data of 43 nearby ($z<0.3$) galaxy mergers, mostly ultraluminous infrared galaxies (ULIRGs) and QSOs. The objects in our sample of BAT AGN have, on average, $\sim 10-100$ times lower AGN luminosities, star formation rates (SFRs), and stellar masses than those of the ULIRG and QSO sample. OH 119 $\mu$m is detected in 42 of our BAT AGN targets. Evidence for molecular outflows (i.e. OH absorption profiles with median velocities more blueshifted than $-$50 km s$^{-1}$ and/or blueshifted wings with 84-percentile velocities less than $-$300 km s$^{-1}$) is seen in only four BAT AGN (NGC~7479 is the most convincing case). Evidence for molecular inflows (i.e. OH absorption profiles with median velocities more redshifted than 50 km s$^{-1}$) is seen in seven objects, although an inverted P-Cygni profile is detected unambiguously in only one object (Circinus). Our data show that both the starburst and AGN contribute to driving OH outflows, but the fastest OH winds require AGN with quasar-like luminosities. We also confirm that the total absorption strength of OH 119 $\mu$m is a good proxy for dust optical depth as it correlates strongly with the 9.7 $\mu$m silicate absorption feature, a measure of obscuration originating in both the nuclear torus and host galaxy disk.

We present high resolution (0."4) IRAM PdBI and ALMA mm and submm observations of the (ultra) luminous infrared galaxies ((U)LIRGs) IRAS17208-0014, Arp220, IC860 and Zw049.057 that reveal intense line emission from vibrationally excited ($\nu_2$=1) J=3-2 and 4-3 HCN. The emission is emerging from buried, compact (r<17-70 pc) nuclei that have very high implied mid-infrared surface brightness $>$$5\times 10^{13}$ L$_{\odot}$ kpc$^{-2}$. These nuclei are likely powered by accreting supermassive black holes (SMBHs) and/or hot (>200 K) extreme starbursts. Vibrational, $\nu_2$=1, lines of HCN are excited by intense 14 micron mid-infrared emission and are excellent probes of the dynamics, masses, and physical conditions of (U)LIRG nuclei when H$_2$ column densities exceed $10^{24}$ cm$^{-2}$. It is clear that these lines open up a new interesting avenue to gain access to the most obscured AGNs and starbursts. Vibrationally excited HCN acts as a proxy for the absorbed mid-infrared emission from the embedded nuclei, which allows for reconstruction of the intrinsic, hotter dust SED. In contrast, we show strong evidence that the ground vibrational state ($\nu$=0), J=3-2 and 4-3 rotational lines of HCN and HCO$^+$ fail to probe the highly enshrouded, compact nuclear regions owing to strong self- and continuum absorption. The HCN and HCO$^+$ line profiles are double-peaked because of the absorption and show evidence of non-circular motions - possibly in the form of in- or outflows. Detections of vibrationally excited HCN in external galaxies are so far limited to ULIRGs and early-type spiral LIRGs, and we discuss possible causes for this. We tentatively suggest that the peak of vibrationally excited HCN emission is connected to a rapid stage of nuclear growth, before the phase of strong feedback.

Powerful winds driven by active galactic nuclei (AGN) are often invoked to play a fundamental role in the evolution of both supermassive black holes (SMBHs) and their host galaxies, quenching star formation and explaining the tight SMBH-galaxy relations. Recent observations of large-scale molecular outflows in ultra-luminous infrared galaxies (ULIRGs) have provided the evidence to support these studies, as they directly trace the gas out of which stars form. Theoretical models suggest an origin of these outflows as energy-conserving flows driven by fast AGN accretion disk winds. Previous claims of a connection between large-scale molecular outflows and AGN activity in ULIRGs were incomplete because they were lacking the detection of the putative inner wind. Conversely, studies of powerful AGN accretion disk winds to date have focused only on X-ray observations of local Seyferts and a few higher redshift quasars. Here we show the clear detection of a powerful AGN accretion disk wind with a mildly relativistic velocity of 0.25c in the X-ray spectrum of IRAS F11119+3257, a nearby (z = 0.189) optically classified type 1 ULIRG hosting a powerful molecular outflow. The AGN is responsible for ~80% of the emission, with a quasar-like luminosity of L_AGN = 1.5x10^46 erg/s. The energetics of these winds are consistent with the energy-conserving mechanism, which is the basis of the quasar mode feedback in AGN lacking powerful radio jets.

We report the detection of far-infrared (FIR) CO rotational emission from nearby active galactic nuclei (AGN) and starburst galaxies, as well as several merging systems and Ultra-Luminous Infrared Galaxies (ULIRGs). Using Herschel-PACS, we have detected transitions in the J$_{upp}$ = 14 - 20 range ($\lambda \sim$ 130 - 185 $\mu$m, $\nu \sim$ 1612 - 2300 GHz) with upper limits on (and in two cases, detections of) CO line fluxes up to J$_{upp}$ = 30. The PACS CO data obtained here provide the first well-sampled FIR extragalactic CO SLEDs for this range, and will be an essential reference for future high redshift studies. We find a large range in the overall SLED shape, even amongst galaxies of similar type, demonstrating the uncertainties in relying solely on high-J CO diagnostics to characterize the excitation source of a galaxy. Combining our data with low-J line intensities taken from the literature, we present a CO ratio-ratio diagram and discuss its potential diagnostic value in distinguishing excitation sources and physical properties of the molecular gas. The position of a galaxy on such a diagram is less a signature of its excitation mechanism, than an indicator of the presence (or absence) of warm, dense molecular gas. We then quantitatively analyze the CO emission from a subset of the detected sources with Large Velocity Gradient (LVG) radiative transfer models to fit the CO SLEDs. Using both single-component and two-component LVG models to fit the kinetic temperature, velocity gradient, number density and column density of the gas, we derive the molecular gas mass and the corresponding CO-to-H$_2$ conversion factor, $\alpha_{CO}$, for each respective source. For the ULIRGs we find $\alpha$ values in the canonical range 0.4 - 5 M$_\odot$/(K kms$^{-1}$pc$^2$), while for the other objects, $\alpha$ varies between 0.2 and 14.} Finally, we compare our best-fit LVG model ..

We obtained high resolution (0."25 to 0."90) observations of HCN and HCO+ J=3-2 of the ultraluminous QSO galaxy Mrk231 with the IRAM Plateau de Bure Interferometer. We find luminous HCN and HCO+ 3-2 emission in the main disk and we detect compact (r<90 pc) vibrationally excited HCN 3-2, v2=1f emission centred on the nucleus. The velocity field of the vibrationally excited HCN is strongly inclined (PA=155 deg.) compared to the east-west rotation of the main disk. The nuclear molecular mass is estimated to 8e8 Msun with an average N(H2)of 1.2e24 cm-2. Prominent, spatially extended (>350 pc) line wings are found for HCN 3-2 with velocities +-750 km/s. Line ratios indicate that the emission is emerging in dense gas n=1e4 - 5e5 cm-3 of elevated HCN abundance X(HCN)=1e-8 to 1e-6. High X(HCN) also allows for the emission to originate in gas of more moderate density. We tentatively detect nuclear emission from the reactive ion HOC+ with HCO+/HOC+=10-20. The HCN v2=1f line emission is consistent with the notion of a hot, dusty, warped inner disk of Mrk231 where the v2=1f line is excited by bright mid-IR 14 micron continuum. We estimate the vibrational temperature T_vib to 200-400 K. We propose that 50% of the main HCN emission may have its excitation affected by the radiation field through IR pumping of the vibrational ground state. The HCN emission in the line wings, however, is more extended and thus likely not strongly affected by IR pumping. Our results reveal that dense clouds survive (and/or are formed) in the AGN outflow on scales of at least several hundred pc before evaporating or collapsing. The elevated HCN abundance in the outflow is consistent with warm chemistry possibly related to shocks and/or X-ray irradiated gas. An upper limit to the mass and momentum flux is 4e8 Msun and 12L_AGN/c, respectively, and we discuss possible driving mechanisms for the dense outflow.

We present FIR-CO luminosity relations ($\log L_{\rm FIR} = \alpha \log L'_{\rm CO} + \beta$) for the full CO rotational ladder from J=1-0 to J=13-12 for 62 local (z < 0.1) (Ultra) Luminous Infrared Galaxies (LIRGs) using data from Herschel SPIRE-FTS and ground-based telescopes. We extend our sample to high redshifts (z > 1) by including 35 (sub)-millimeter selected dusty star forming galaxies from the literature with robust CO observations. The addition of luminous starbursts at high redshifts enlarge the range of the FIR-CO luminosity relations towards the high-IR-luminosity end while also significantly increasing the small amount of mid-/high-J CO line data available prior to Herschel. This new data-set (both in terms of IR luminosity and J-ladder) reveals linear FIR-CO luminosity relations ($\alpha \sim 1$) for J=1-0 up to J=5-4, with a nearly constant normalisation ($\beta \sim 2$). This is expected from the (also) linear FIR-(molecular line) relations found for the dense gas tracer lines (HCN and CS), as long as the dense gas mass fraction does not vary strongly within our (merger/starburst)-dominated sample. However from J=6-5 and up to J=13-12 we find an increasingly sub-linear slope and higher normalization constant with increasing J. We argue that these are caused by a warm (~100K) and dense ($>10^4{\rm cm^{-3}}$) gas component whose thermal state is unlikely to be maintained by star formation powered far-UV radiation fields (and thus is no longer directly tied to the star formation rate). We suggest that mechanical heating (e.g., supernova driven turbulence and shocks), and not cosmic rays, is the more likely source of energy for this component. The global CO spectral line energy distributions (SLEDs), which remain highly excited from J=6-5 up to J=13-12, are found to be a generic feature of the (U)LIRGs in our sample, and further support the presence of this gas component.

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 report the results from a systematic search for molecular (OH-119 um) outflows with Herschel-PACS in a sample of 43 nearby (z < 0.3) galaxy mergers, mostly ultraluminous infrared galaxies (ULIRGs) and QSOs. We find that the character of the OH feature (strength of the absorption relative to the emission) correlates with that of the 9.7-um silicate feature, a measure of obscuration in ULIRGs. Unambiguous evidence for molecular outflows, based on the detection of OH absorption profiles with median velocities more blueshifted than -50 km/sec, is seen in 26 (70%) of the 37 OH-detected targets, suggesting a wide-angle (~145 degrees) outflow geometry. Conversely, unambiguous evidence for molecular inflows, based on the detection of OH absorption profiles with median velocities more redshifted than +50 km/sec, is seen in only 4 objects, suggesting a planar or filamentary geometry for the inflowing gas. Terminal outflow velocities of ~-1000 km/sec are measured in several objects, but median outflow velocities are typically ~-200 km s^{-1}. While the outflow velocities show no statistically significant dependence on the star formation rate, they are distinctly more blueshifted among systems with large AGN fractions and luminosities [log (L_AGN / L_sun) > 11.8 +/- 0.3]. The quasars in these systems play a dominant role in driving the molecular outflows. In contrast, the most AGN dominated systems, where OH is seen purely in emission, show relatively modest OH line widths, despite their large AGN luminosities, perhaps indicating that molecular outflows subside once the quasar has cleared a path through the obscuring material.

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.

We present Herschel SPIRE FTS spectroscopy of the nearby luminous infrared galaxy NGC 6240. In total 20 lines are detected, including CO J=4-3 through J=13-12, 6 H2O rotational lines, and [CI] and [NII] fine-structure lines. The CO to continuum luminosity ratio is 10 times higher in NGC 6240 than Mrk 231. Although the CO ladders of NGC 6240 and Mrk 231 are very similar, UV and/or X-ray irradiation are unlikely to be responsible for the excitation of the gas in NGC 6240. We applied both C and J shock models to the H2 v=1-0 S(1) and v=2-1 S(1) lines and the CO rotational ladder. The CO ladder is best reproduced by a model with shock velocity v_s=10 km s^-1 and a pre-shock density n_H=5 * 10^4 cm^-3. We find that the solution best fitting the H2 lines is degenerate: The shock velocities and number densities range between v_s = 17 - 47 km s^-1 and n_H=10^7 - 5 * 10^4 cm^-3, respectively. The H2 lines thus need a much more powerful shock than the CO lines. We deduce that most of the gas is currently moderately stirred up by slow (10 km s^-1) shocks while only a small fraction (< 1 percent) of the ISM is exposed to the high velocity shocks. This implies that the gas is rapidly loosing its highly turbulent motions. We argue that a high CO line-to-continuum ratio is a key diagnostic for the presence of shocks.

Star formation is arguably the most important physical process in the cosmos. It is a fundamental driver of galaxy evolution and the ultimate source of most of the energy emitted by galaxies. A correct interpretation of star formation rate (SFR) measures is therefore essential to our understanding of galaxy formation and evolution. Unfortunately, however, no single SFR estimator is universally available or even applicable in all circumstances: the numerous galaxies found in deep surveys are often too faint (or too distant) to yield significant detections with most standard SFR measures, and until now there have been no global, multi-band observations of nearby galaxies that span all the conditions under which star-formation is taking place. To address this need in a systematic way, we have undertaken a multi-band survey of all types of star-forming galaxies in the local Universe. This project, the Star Formation Reference Survey (SFRS), is based on a statistically valid sample of 369 nearby galaxies that span all existing combinations of dust temperature, SFR, and specific SFR. Furthermore, because the SFRS is blind with respect to AGN fraction and environment it serves as a means to assess the influence of these factors on SFR. Our panchromatic global flux measurements (including GALEX FUV+NUV, SDSS ugriz, 2MASS JHKs, Spitzer 3-8{\mu}m, and others) furnish uniform SFR measures and the context in which their reliability can be assessed. This paper describes the SFRS survey strategy, defines the sample, and presents the multi-band photometry collected to date.

We analyse the far infrared dust emission from the Galactic Centre region, including the Circumnuclear Disk (CND) and other structures, using Herschel PACS and SPIRE photometric observations. These Herschel data are complemented by unpublished observations by the Infrared Space Observatory Long Wavelength Spectrometer (ISO LWS), which used parallel mode scans to obtain photometric images of the region with a larger beam than Herschel but with a complementary wavelength coverage and more frequent sampling with ten detectors observing at ten different wavelengths in the range from 46 to 180 \mum, where the emission peaks. We also include data from the MSX at 21.3 \mum for completeness. We model the combined ISO LWS continuum plus Herschel PACS and SPIRE photometric data toward the central 2 pc in Sgr A*, a region that includes the CND. We find that the FIR spectral energy distribution is best represented by a continuum that is the sum of three greybody curves from dust at temperatures of 90, 44.5, and 23 K. We obtain temperature and molecular hydrogen column density maps of the region. We estimate the mass of the inner part of the CND to be ~5.0x10e+4 Msum, with luminosities: Lcavity ~2.2x10e+6 Lsun and Lcnd ~1.5x10e+6 Lsun in the central 2 pc radius around SgrA* . We find from the Herschel and ISO data that the cold component of the dust dominates the total dust mass, with a contribution of ~3.2x10E+4 Msun; this important cold material had escaped the notice of earlier studies that relied on shorter wavelength observations. The hotter component disagrees with some earlier estimates, but is consistent with measured gas temperatures and with models that imply shock heating or turbulent effects are at work. We find that the dust grain sizes apparently change widely across the region, perhaps in response to the temperature variations, and we map that distribution.

We present a full high resolution SPIRE FTS spectrum of the nearby ultraluminous infrared galaxy Mrk231. In total 25 lines are detected, including CO J=5-4 through J=13-12, 7 rotational lines of H2O, 3 of OH+ and one line each of H2O+, CH+, and HF. We find that the excitation of the CO rotational levels up to J=8 can be accounted for by UV radiation from star formation. However, the approximately flat luminosity distribution of the CO lines over the rotational ladder above J=8 requires the presence of a separate source of excitation for the highest CO lines. We explore X-ray heating by the accreting supermassive black hole in Mrk231 as a source of excitation for these lines, and find that it can reproduce the observed luminosities. We also consider a model with dense gas in a strong UV radiation field to produce the highest CO lines, but find that this model strongly overpredicts the hot dust mass in Mrk231. Our favoured model consists of a star forming disk of radius 560 pc, containing clumps of dense gas exposed to strong UV radiation, dominating the emission of CO lines up to J=8. X-rays from the accreting supermassive black hole in Mrk231 dominate the excitation and chemistry of the inner disk out to a radius of 160 pc, consistent with the X-ray power of the AGN in Mrk231. The extraordinary luminosity of the OH+ and H2O+ lines reveals the signature of X-ray driven excitation and chemistry in this region.

The first high resolution Spitzer IRS 9-37um spectra of 29 Seyfert galaxies (about one quarter) of the 12um Active Galaxy Sample are presented and discussed. The high resolution spectroscopy was obtained with corresponding off-source observations. This allows excellent background subtraction, so that the continuum levels and strengths of weak emission lines are accurately measured. The result is several new combinations of emission line ratios, line/continuum and continuum/continuum ratios that turn out to be effective diagnostics of the strength of the AGN component in the IR emission of these galaxies. The line ratios [NeV]/[NeII], [OIV]/[NeII], already known, but also [NeIII]/[NeII] and [NeV]/[SiII] can all be effectively used to measure the dominance of the AGN. We extend the analysis, already done using the 6.2um PAH emission feature, to the equivalent width of the 11.25um PAH feature, which also anti-correlates with the dominance of the AGN. We measure that the 11.25um PAH feature has a constant ratio with the H_2 S(1) irrespective of Seyfert type, approximately 10 to 1. Using the ratio of accurate flux measurements at about 19um with the two spectrometer channels, having aperture areas differing by a factor 4, we measured the source extendness and correlated it with the emission line and PAH feature equivalent widths. The extendness of the source gives another measure of the AGN dominance and correlates both with the EWs of [NeII] and PAH emission. Using the rotational transitions of H$_2$ we were able to estimate temperatures (200-300K) and masses (1-10 x 10^6 M_sun), or significant limits on them, for the warm molecular component in the galaxies observed.

The first complete far-infrared spectrum (43-197um) of the Seyfert 2 galaxy NGC 1068 as observed with the LWS onboard ISO is reported. In addition to the expected ionic fine structure emission lines, the OH rotational lines at 79, 119 and 163um were all detected in emission. The line intensities were modelled together with ISO SWS and optical and UV line intensities from the literature, considering two independent emission components: the AGN component and the starburst component in the circumnuclear ring of 3 kpc in size. Using the UV to mid-IR emission line spectrum to constrain the nuclear ionizing continuum, we have confirmed previous results: a canonical power-law ionizing spectrum is a poorer fit than one with a deep absorption trough, while the presence of a big blue bump is ruled out. Based on starburst synthesis models, we were able to fit the mid- and far-infrared fine-structure line emission, with a low ionization parameter (logU=-3.5) and low densities (n=100 cm^-3). Combining AGN and starburst components, we succeeded in modeling the overall UV to far-IR atomic spectrum of NGC1068, reproducing the line fluxes to within a factor 2.0 on average with a standard deviation of 1.3 and the overall continuum as the sum of the contribution of the thermal dust emission in the ionized and neutral components. TheOH 119 um emission indicates that the line is collisionally excited and arises in a warm and dense region. The OH emission has been modeled using spherically symmetric, non-local, non-LTE radiative transfer models. The models indicate that the bulk of the emission arises from the nuclear region, although some extended contribution from the starburst is not ruled out. The OH abundance in the nuclear region is expected to be $\sim10^{-5}$, characteristic of X-ray dominated regions.

The formation and evolution of young low mass stars are characterized by important processes of mass loss and accretion ocurring in the innermost regions of their placentary circumstellar disks. Because of the large obscuration of these disks at optical and infrared wavelengths in the early protostellar stages (Class 0 Sources), they were previously detected only at radio wavelengths using interferometric techniques. We have detected with the Infrared Space Observatory (ISO) the mid-infrared emission associated with the Class 0 protostar VLA1 in the HH1-2 region located in the Orion nebula. The emission arises in the three wavelength windows at 5.3, 6.6 and 7.5 micras where the absorption due to ices and silicates has a local minimum that exposes the central parts of the youngest protostellar systems to mid-infrared investigations. The mid-infrared emission arises from a central source with 4 AU diameter at an averaged temperature of 700 K, deeply embedded in a dense region with a visual extinction of Av=80-100mag.

First detections of thermal water vapor absorption lines have been made toward Orion IRc2 using the Short Wavelength Spectrometer (SWS) on board the Infrared Space Observatory (ISO). Grating spectra covering wavelengths 25-45 micron yield 19 pure rotational lines, originating from energy levels 200-750 K above ground. Fabry-Perot spectra of 5 transitions resolve the line profiles and reveal the water gas kinematics. The fact that all lines are seen in absorption is in striking contrast with data from the ISO Long Wavelength Spectrometer (LWS), where the water lines appear in emission. At least one line displays a P-Cygni type profile, which suggests that the water is located in an expanding shell centered on or near IRc2. The expansion velocity is 18 km per second, in agreement with the value inferred from water maser observations by Genzel et al. (1981). Because the continuum is intense and likely formed in or near the water-containing gas, the excitation of the observed transitions is dominated by radiative processes. A simple, generalised curve-of-growth method is presented and used to analyze the data. A mean excitation temperature of 72 K and a total water column density of 1.5e18 molecules per cm squared are inferred, each with an estimated maximum uncertainty of 20%. Combined with the molecular hydrogen column density derived from ISO observations of the pure rotational lines, and an assumed temperature of 200-350 K, the inferred water abundance is 2e-4 to 5e-4 in the warm shocked gas. This abundance is similar to that found recently by Harwit et al. (1998) toward Orion using data from the LWS, but higher than that found for most other shocked regions by, for example, Liseau et al. (1996).

We report on sensitive far-infrared observations of $^{12}$CO pure rotational transitions in the OMC-1 core of Orion. The lines were observed with the Long Wavelength Spectrometer (LWS) in the grating mode on board the Infrared Space Observatory (ISO), covering the 43-197 $\mu$m wavelength range. The transitions from $J_{up}=14$ up to $J_{up}=19$ have been identified across the whole OMC-1 core and lines up to $J_{up}= 43$ have been detected towards the central region, KL/IRc2. In addition, we have taken high-quality spectra in the Fabry-Perot mode of some of the CO lines. In KL/IRc2 the lines are satisfactorily accounted for by a three-temperature model describing the plateau and ridge emission. The fluxes detected in the high-$J$ transitions ($J_{up} > 34$) reveal the presence of a very hot and dense gas component ($T=1500-2500$ K; $\rm N(CO)$=2\times 10^{17}\cmmd$), probably originating from some of the embedded sources previously observed in the $\rm H_2$ near-infrared lines. At all other positions in the OMC-1 core, we estimate kinetic temperatures $\geq 80$ K and as high as 150 K at some positions around IRc2, from a simple Large-Velocity Gradient model.

We present observations of the 5$_{1,5}\to4_{2,2}$ transition of water vapor at 325.15 GHz taken with the CSO telescope towards Orion IRc2. The emission is more extended than that of other molecular species such as CH$_3$OH. However, it is much less extended than the emission of water vapor at 183.31 GHz reported by Cernicharo et al (1994). A comparison of the line intensities at 325.15 GHz and 183.31 GHz puts useful constraints on the density and temperature of the emitting regions and allows an estimate of H_2O abundance, x(H_2O), of $\simeq10^{-4}$ in the Plateau and $\simeq10^{-6}-10^{-5}$ in the Ridge.