Whenever gas is blown out of a galaxy, chances are that it contains some cosmic dust. This dust is an important part of the metals budget for the circumgalactic and intergalactic media, and traces the outflow and stripping history of the galaxy. The dust is also interesting in its own right, as dust plays an essential role in many astrophysical processes. We have only begun to learn about circumgalactic dust, and in particular we do not know how (and when) it escapes its host galaxy. Here we describe the prospects for measuring the dust mass and properties around many individual galaxies, which will form the basis for understanding how the dust got there.

We present the results from an analysis of deep Herschel Space Observatory observations of six nearby dwarf galaxies known to host galactic-scale winds. The superior far-infrared sensitivity and angular resolution of Herschel have allowed detection of cold circumgalactic dust features beyond the stellar components of the host galaxies traced by Spitzer 4.5 $\mu$m images. Comparisons of these cold dust features with ancillary data reveal an imperfect spatial correlation with the ionized gas and warm dust wind components. We find that typically $\sim$10-20\% of the total dust mass in these galaxies resides outside of their stellar disks, but this fraction reaches $\sim$60\% in the case of NGC 1569. This galaxy also has the largest metallicity (O/H) deficit in our sample for its stellar mass. Overall, the small number of objects in our sample precludes drawing strong conclusions on the origin of the circumgalactic dust. We detect no statistically significant trends with star formation properties of the host galaxies, as might be expected if the dust were lifted above the disk by energy inputs from on-going star formation activity. Although a case for dust entrained in a galactic wind is seen in NGC 1569, in all cases, we cannot rule out the possibility that some of the circumgalactic dust might be associated instead with gas accreted or removed from the disk by recent galaxy interaction events, or that it is part of the outer gas-rich portion of the disk that lies below the sensitivity limit of the Spitzer 4.5 $\mu$m data.

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 combine the Herschel Space Observatory PACS and SPIRE photometry with archival WISE photometry to construct the spectral energy distributions (SED) for over 300 local ($z < 0.05$), ultra-hard X-ray (14 - 195 keV) selected active galactic nuclei (AGN) from the Swift Burst Alert Telescope (BAT) 58 month catalogue. Using a simple analytical model that combines an exponentially cut-off powerlaw with a single temperature modified blackbody, we decompose the SEDs into a host-galaxy and AGN component. We calculate dust masses, dust temperatures, and star-formation rates (SFR) for our entire sample and compare them to a stellar mass-matched sample of local non-AGN galaxies. We find AGN host galaxies have systematically higher dust masses, dust temperatures, and SFRs due to the higher prevalence of late-type galaxies to host an AGN, in agreement with previous studies of the Swift/BAT AGN. We provide a scaling to convert X-ray luminosities into 8 - 1000 $\mu$m AGN luminosities, as well as determine the best mid-to-far IR colors for identifying AGN dominated galaxies in the IR regime. We find that for nearly 30 per cent of our sample, the 70 $\mu$m emission contains a significant contribution from the AGN ($> 0.5$), especially at higher luminosities ($L_{14-195\,\rm{keV}} > 10^{42.5}$ ergs s$^{-1}$). Finally, we measure the local SFR-AGN luminosity relationship, finding a slope of 0.18, large scatter (0.37 dex), and no evidence for an upturn at high AGN luminosity. We conclude with a discussion on the implications of our results within the context of galaxy evolution with and without AGN feedback.

Context. The LIRG Zw 049.057 contains a compact obscured nucleus where a considerable amount of the galaxy's luminosity is generated. This nucleus contains a dusty environment that is rich in molecular gas. One approach to probing this kind of environment and to revealing what is hidden behind the dust is to study the rotational lines of molecules that couple well with the IR radiation emitted by the dust. Methods. We observed Zw 049.057 with PACS and SPIRE onboard the Herschel Space Observatory in rotational lines of H2O, H218O, OH, 18OH, and [O I]. We modeled the unresolved core of the galaxy using a spherically symmetric radiative transfer code. Results. We present the full SPIRE FTS spectrum of Zw 049.057, along with relevant spectral scans in the PACS range. We find that a minimum of two different components (nuclear and extended) are required in order to account for the rich molecular line spectrum. The nuclear component has a radius of 10-30 pc, a very high infrared surface brightness (1e14 Lsun/kpc2), warm dust (Td > 100 K), and a very large H2 column density (NH2 = 1e24-1e25 cm-2). The modeling also indicates high nuclear H2O (5e-6) and OH (4e-6) abundances relative to H2 as well as a low 16O/18O-ratio of 50-100. We also find a prominent infall signature in the [O I] line. We tentatively detect a 500 km/s outflow in the H2O 313->202 line. Conclusions. The high surface brightness of the core indicates the presence of either a buried active galactic nucleus or a very dense nuclear starburst.nThe H2O abundance is comparable to that of other compact (ultra-)luminous infrared galaxies such as NGC 4418 and Arp 220 - and also to hot cores in the Milky Way. The enhancement of 18O is a possible indicator that the nucleus of Zw 049.057 is in a similar evolutionary stage as the nuclei of Arp 220 - and more advanced than NGC 4418.

We present a detailed analysis of deep far-infrared observations of the nearby edge-on star-forming galaxy NGC 4631 obtained with the Herschel Space Observatory. Our PACS images at 70 and 160 um show a rich complex of filaments and chimney-like features that extends up to a projected distance of 6 kpc above the plane of the galaxy. The PACS features often match extraplanar Halpha, radio-continuum, and soft X-ray features observed in this galaxy, pointing to a tight disk-halo connection regulated by star formation. On the other hand, the morphology of the colder dust component detected on larger scale in the SPIRE 250, 350, and 500 um data matches the extraplanar H~I streams previously reported in NGC 4631 and suggests a tidal origin. The PACS 70/160 ratios are elevated in the central ~3.0 kpc region above the nucleus of this galaxy (the "superbubble"). A pixel-by-pixel analysis shows that dust in this region has a higher temperature and/or an emissivity with a steeper spectral index (beta > 2) than the dust in the disk, possibly the result of the harsher environment in the superbubble. Star formation in the disk seems energetically insufficient to lift the material out of the disk, unless it was more active in the past or the dust-to-gas ratio in the superbubble region is higher than the Galactic value. Some of the dust in the halo may also have been tidally stripped from nearby companions or lifted from the disk by galaxy interactions.

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.

Far-Infrared (FIR) photometry from the the Photodetector Array Camera and Spectrometer (PACS) on the Herschel Space Observatory is presented for 313 nearby, hard X-ray selected galaxies from the 58-month Swift Burst Alert Telescope (BAT) Active Galactic catalog. The present data do not distinguish between the FIR luminosity distributions at 70 and 160um for Seyfert 1 and Seyfert 2 galaxies. This result suggests that if the FIR emission is from the nuclear obscuring material surrounding the accretion disk, then it emits isotropically, independent of orientation. Alternatively, a significant fraction of the 70 and 160um could be from star formation, independent of AGN type. Using a non-parametric test for partial correlation with censored data, we find a statistically significant correlation between the AGN intrinsic power (in the 14-195 keV band ) and the FIR emission at 70 and 160um for Seyfert 1 galaxies. We find no correlation between the 14-195 keV and FIR luminosities in Seyfert 2 galaxies. The observed correlations suggest two possible scenarios: (i) if we assume that the FIR luminosity is a good tracer of star formation, then there is a connection between star formation and the AGN at sub-kiloparsec scales, or (ii) dust heated by the AGN has a statistically significant contribution to the FIR emission. Using a Spearman rank-order analysis, the 14-195 keV luminosities for the Seyfert 1 and 2 galaxies are weakly statistically correlated with the F70/F160 ratios.

We have carried out extensive and detailed photoionization modeling to successfully constrain the locations of different emission-line galaxies in optical and mid-infrared diagnostic diagrams. Our model grids cover a wide range in parameter space for the active galaxy continuum and starburst galaxies with different stellar population laws and metallicities. We compare the predicted AGN and star-formation mid-infrared line ratios [Ne III]15.56mm/[Ne II]12.81mm and [O IV]25.89mm/[Ne III]15.56mm to the observed values, and find that the best fit for the AGN is via a two-zone approximation. This two-zone approximation is a combination of a matter-bounded component, where [Ne III] and [O IV] are emitted efficiently, and a radiation-bounded component that maximizes [Ne II] emission. We overlay the predictions from this two-zone approximation onto the optical [O III]l5007/Hbeta and [N II]l6583/Halpha diagnostic diagram derived from the Sloan Digital Sky Survey, to find that the high-density and low-ionization radiation-bounded component in our two-zone AGN approximation model provides a good lower limit for [N II] emission. This establishes a new theoretical demarcation line for the minimum AGN contribution in this diagram. This new classification results by a factor of ~1.4 in a higher AGN population than predictions derived from previous divisions of star-forming galaxies. Similarly, we define a maximum AGN contribution in the [O III]/Hbeta and [N II]/Halpha diagram by using a two-zone approximation within a parameter range typical of the narrow-line region.

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 Effective Eddington Limit for dusty gas surrounding AGN is lower than the canonical Eddington limit for hydrogen gas. Previous results from the Swift/BAT 9-month catalogue suggested that in the overwhelming majority of local AGN, the dusty absorbing gas is below this Effective Eddington limit, implying that radiation pressure is insufficient to blow away the absorbing clouds. We present an analysis of three objects from that sample which were found to be close to the Effective Eddington limit (NGC454, 2MASX J03565655-4041453 and XSS J05054-2348), using newly obtained XMM-Newton data. We use the X-ray data to better constrain the absorbing column density, and supplement them with XMM optical monitor (OM) data, infrared Spitzer and Herschel data where available to construct a broad-band spectral energy distribution to estimate refined bolometric luminosities and Eddington ratios for these three objects. The new XMM-Newton observations show all three objects moving away from the region expected for short-lived absorption in the N_H-\lambda_{Edd} plane into the `long-lived absorption' region. We find our conclusions robust to different methods for estimating the bolometric luminosity and Eddington ratio. Interestingly, 2MASX J03565655-4041453 and XSS J05054-2348 now exhibit complex X-ray spectra, at variance with previous analyses of their Swift/XRT data. We find evidence for absorption variability in NGC 454 and 2MASX J03565655-4041453, perhaps implying that although the radiation pressure from the central engine is insufficient to cause clearly detectable outflows, it may cause absorption variations over longer timescales. However, more robust black hole mass estimates would improve the accuracy of the Eddington ratio estimates for these objects.

The blue compact dwarf galaxy I Zw 18 is one of the most metal poor systems known in the local Universe (12 + log(O/H) $=$ 7.17). In this work we study I Zw 18 using data from {\it Spitzer}, {\it Herschel Space Telescope} and IRAM Plateau de Bure Interferometer. Our data set includes the most sensitive maps of I Zw 18, to date, in both, the far infrared and the CO $J=1\rightarrow0$ transition. We use dust emission models to derive a dust mass upper limit of only M$_{dust}\leq1.1\times10^4$ M$_{\odot}$ ($3\sigma$ limit). This upper limit is driven by the non-detection at 160 $\mu$m, and it is a factor of 4-10 times smaller than previous estimates (depending upon the model used). We also estimate an upper limit to the total dust-to-gas mass ratio of M$_{Dust}$/M$_{gas}\leq5.0\times10^{-5}$. If a linear correlation between the dust-to-gas mass ratio and metallicity (measure as O/H) were to hold, we would expect a ratio of 3.9$\times10^{-4}$. We also show that the infrared SED is similar to that of starbursting systems.

The shape of the spectral energy distribution of active galaxies in the EUV--soft X-ray band (13.6 eV to 1 keV) is uncertain because obscuration by dust and gas can hamper our view of the continuum. To investigate the shape of the spectral energy distribution in this energy band, we have generated a set of photoionization models which reproduce the small dispersion found in correlations between high-ionization mid-infrared emission lines in a sample of hard X-ray selected AGN. Our calculations show that a broken power-law continuum model is sufficient to reproduce the [Ne V]14.32 mm/[NeIII], [Ne V]24.32mm/[O IV]25.89mm and [O IV] 25.89mm/[Ne III] ratios, and does not require the addition of a "big bump" EUV model component. We constrain the EUV--soft X-ray slope, alpha_i, to be between 1.5 -- 2.0 and derive a best fit of alpha_i ~ 1.9 for Seyfert 1 galaxies, consistent with previous studies of intermediate redshift quasars. If we assume a blue bump model, most sources in our sample have derived temperatures between T_{BB}=10^{5.18} K to 10^{5.7} K, suggesting that the peak of this component spans a large range of energies extending from ~ lambda 600A to lambda 1900A. In this case, the best fitting peak energy that matches the mid-infrared line ratios of Seyfert 1 galaxies occurs between ~ lambda 700--1000A. Despite the fact that our results do not rule out the presence of an EUV bump, we conclude that our power-law model produces enough photons with energies > 4 Ry to generate the observed amount of mid-infrared emission in our sample of BAT AGN.

We compare mid-infrared emission-line properties, from high-resolution Spitzer spectra of a hard X-ray (14 -- 195 keV) selected sample of nearby (z < 0.05) AGN detected by the Burst Alert Telescope (BAT) aboard Swift. The luminosity distribution for the mid-infrared emission-lines, [O IV] 25.89 micron, [Ne II] 12.81 micron, [Ne III] 15.56 micron and [Ne V] 14.32/24.32 micron, and hard X-ray continuum show no differences between Seyfert 1 and Seyfert 2 populations, however six newly discovered BAT AGNs are under-luminous in [O IV], most likely the result of dust extinction in the host galaxy. The overall tightness of the mid-infrared correlations and BAT fluxes and luminosities suggests that the emission lines primarily arise in gas ionized by the AGN. We also compare the mid-infrared emission-lines in the BAT AGNs with those from published studies of ULIRGs, PG QSOs, star-forming galaxies and LINERs. We find that the BAT AGN sample fall into a distinctive region when comparing the [Ne III]/[Ne II] and the [O IV]/[Ne III] ratios. These line ratios are lower in sources that have been previously classified in the mid-infrared/optical as AGN than those found for the BAT AGN, suggesting that, in our X-ray selected sample, the AGN represents the main contribution to the observed line emission. These ratios represent a new emission line diagnostic for distinguishing between AGN and star forming galaxies.

We have studied the relationship between the nuclear (high-resolution) radio emission, at 8.4 GHz (3.6 cm) and 1.4 GHz (20 cm), the [O IV] 25.89um, [Ne III] 15.56um and [Ne II] 12.81um emission lines and the black hole mass accretion rate for a sample of Seyfert galaxies. In order to characterize the radio contribution for the Seyfert nuclei we used the 8.4GHz/[O IV] ratio, assuming that [O IV] scales with the luminosity of the AGN. From this we find that Seyfert 1's (i.e., Seyfert 1.0's, 1.2's, and 1.5's) and Seyfert 2's (i.e., Seyfert 1.8's, 1.9's, and 2.0's) have similar radio contributions, relative to the AGN. On the other hand, sources in which the [Ne II] emission is dominated either by the AGN or star formation have statistically different radio contributions, with star formation dominated sources more "radio loud", by a factor of ~2.8 on average, than AGN dominated sources. We show that star formation dominated sources with relatively larger radio contribution have smaller mass accretion rates. Overall, we suggest that 8.4GHz/[O IV], or alternatively, 1.4GHz/[O IV] ratios, can be used to characterize the radio contribution, relative to the AGN, without the limitation of previous methods that rely on optical observables.

We have studied the relationship between the high- and low-ionization [O IV] lambda 25.89 micron, [Ne III] lambda 15.56 micron and [Ne II] lambda 12.81 micron emission lines with the aim of constraining the active galactic nuclei (AGN) and star formation contributions for a sample of 103 Seyfert galaxies. We used the [O IV] and [Ne II] emission as tracers for the AGN power and star formation to investigate the ionization state of the emission-line gas. We find that Seyfert 2 galaxies have, on average, lower [O IV]/[Ne II] ratios than those of Seyfert 1 galaxies. This result suggests two possible scenarios: 1) Seyfert 2 galaxies have intrinsically weaker AGN, or 2) Seyfert 2 galaxies have relatively higher star formation rates than Seyfert 1 galaxies. We estimate the fraction of [Ne II] directly associated with the AGN and find that Seyfert 2 galaxies have a larger contribution from star formation, by a factor of ~1.5 on average, than what is found in Seyfert 1 galaxies. Using the stellar component of [Ne II] as a tracer of the current star formation we found similar star formation rates in Seyfert 1 and Seyfert 2 galaxies. We examined the mid- and far-infrared continua and find that [Ne II] is well correlated with the continuum luminosity at 60 micron and that both [Ne III] and [O IV] are better correlated with the 25 micron luminosities than with the continuum at longer wavelengths, suggesting that the mid-infrared continuum luminosity is dominated by the AGN, while the far-infrared luminosity is dominated by star formation. Overall, these results test the unified model of AGN, and suggest that the differences between Seyfert galaxies cannot be solely due to viewing angle dependence.

We have studied the relationship between the [O IV] lambda 25.89 micron emission line luminosities, obtained from Spitzer spectra, the X-ray continua in the 2-10 keV band, primarily from ASCA, and the 14-195 keV band obtained with the SWIFT/Burst Alert Telescope (BAT), for a sample of nearby (z < 0.08) Seyfert galaxies. For comparison, we have examined the relationship between the [O III] 5007, the 2-10 keV and the 14-195 keV luminosities for the same set of objects. We find that both the [O IV] and [O III] luminosities are well-correlated with the BAT luminosities. On the other hand, the [O III] luminosities are better-correlated with 2-10 keV luminosities than are those of [O IV]. When comparing [O IV] and [O III] luminosities for the different types of galaxies, we find that the Seyfert 2's have significantly lower [O III] to [O IV] ratios than the Seyfert 1's. We suggest that this is due to more reddening of the narrow line region (NLR) of the Seyfert 2's. Assuming Galactic dust to gas ratios, the average amount of extra reddening corresponds to a hydrogen column density of ~ few times 10^21 cm^-2, which is a small fraction of the X-ray absorbing columns in the Seyfert 2's. The combined effects of reddening and the X-ray absorption are the probable reason why the [O III] versus 2-10 keV correlation is better than the [O IV] versus 2-10 keV, since the [O IV] emission line is much less affected by extinction. Overall, we find the [O IV] to be an accurate and truly isotropic indicator of the power of the AGN. This suggests that it can be useful in deconvolving the contribution of the AGN and starburst to the spectrum of Compton-thick and/or X-ray weak sources.