Searching for molecular outflows in Hyper-Luminous Infrared Galaxies

We present constraints on the molecular outflows in a sample of five Hyper-Luminous Infrared Galaxies using Herschel observations of the OH doublet at 119 {\mu}m. We have detected the OH doublet in three cases: one purely in emission and two purely in absorption. The observed emission profile has a significant blueshifted wing suggesting the possibility of tracing an outflow. Out of the two absorption profiles, one seems to be consistent with the systemic velocity while the other clearly indicates the presence of a molecular outflow whose maximum velocity is about ~1500 km/s. Our analysis shows that this system is in general agreement with previous results on Ultra-luminous Infrared Galaxies and QSOs, whose outflow velocities do not seem to correlate with stellar masses or starburst luminosities (star formation rates). Instead the galaxy outflow likely arises from an embedded AGN.Comment: Accepted for publication in MNRAS. 13 pages, 11 figures, 4 table

Chemically Distinct Nuclei and Outflowing Shocked Molecular Gas in Arp 220

We present the results of interferometric spectral line observations of Arp 220 at 3.5mm and 1.2mm from the Plateau de Bure Interferometer (PdBI), imaging the two nuclear disks in H$^{13}$CN$(1 - 0)$ and $(3 - 2)$, H$^{13}$CO$^+(1 - 0)$ and $(3 - 2)$, and HN$^{13}$C$(3 - 2)$ as well as SiO$(2 - 1)$ and $(6 - 5)$, HC$^{15}$N$(3 - 2)$, and SO$(6_6 - 5_5)$. The gas traced by SiO$(6 - 5)$ has a complex and extended kinematic signature including a prominent P Cygni profile, almost identical to previous observations of HCO$^+(3 - 2)$. Spatial offsets $0.1''$ north and south of the continuum centre in the emission and absorption of the SiO$(6 - 5)$ P Cygni profile in the western nucleus (WN) imply a bipolar outflow, delineating the northern and southern edges of its disk and suggesting a disk radius of $\sim40$ pc, consistent with that found by ALMA observations of Arp 220. We address the blending of SiO$(6 - 5)$ and H$^{13}$CO$^+(3 - 2)$ by considering two limiting cases with regards to the H$^{13}$CO$^+$ emission throughout our analysis. Large velocity gradient (LVG) modelling is used to constrain the physical conditions of the gas and to infer abundance ratios in the two nuclei. Our most conservative lower limit on the [H$^{13}$CN]/[H$^{13}$CO$^+$] abundance ratio is 11 in the WN, cf. 0.10 in the eastern nucleus (EN). Comparing these ratios to the literature we argue on chemical grounds for an energetically significant AGN in the WN driving either X-ray or shock chemistry, and a dominant starburst in the EN.Comment: 28 pages, 17 figures, accepted to Ap

Spectroscopic FIR mapping of the disk and galactic wind of M82 with Herschel-PACS

[Abridged] We present maps of the main cooling lines of the neutral atomic gas ([OI] at 63 and 145 micron and [CII] at 158 micron) and in the [OIII] 88 micron line of the starburst galaxy M82, carried out with the PACS spectrometer on board the Herschel satellite. By applying PDR modeling we derive maps of the main ISM physical parameters, including the [CII] optical depth, at unprecedented spatial resolution (~300 pc). We can clearly kinematically separate the disk from the outflow in all lines. The [CII] and [OI] distributions are consistent with PDR emission both in the disk and in the outflow. Surprisingly, in the outflow, the atomic and the ionized gas traced by the [OIII] line both have a deprojected velocity of ~75 km/s, very similar to the average velocity of the outflowing cold molecular gas (~ 100 km/s) and several times smaller than the outflowing material detected in Halpha (~ 600 km/s). This suggests that the cold molecular and neutral atomic gas and the ionized gas traced by the [OIII] 88 micron line are dynamically coupled to each other but decoupled from the Halpha emitting gas. We propose a scenario where cold clouds from the disk are entrained into the outflow by the winds where they likely evaporate, surviving as small, fairly dense cloudlets (n_H\sim 500-1000 cm^-3, G_0\sim 500- 1000, T_gas\sim300 K). We show that the UV photons provided by the starburst are sufficient to excite the PDR shells around the molecular cores. The mass of the neutral atomic gas in the outflow is \gtrsim 5-12x 10^7 M_sun to be compared with that of the molecular gas (3.3 x 10^8 M_sun) and of the Halpha emitting gas (5.8 x 10^6 M_sun). The mass loading factor, (dM/dt)/SFR, of the molecular plus neutral atomic gas in the outflow is ~ 2. Energy and momentum driven outflow models can explain the data equally well, if all the outflowing gas components are taken into account.Comment: 26 pages, 23 figures, 4 Tables, Accepted for publication in Astronomy & Astrophysic

A deep Herschel/PACS observation of CO(40-39) in NGC 1068: a search for the molecular torus

Emission from high-J CO lines in galaxies has long been proposed as a tracer of X-ray dominated regions (XDRs) produced by AGN. Of particular interest is the question of whether the obscuring torus, which is required by AGN unification models, can be observed via high-J CO cooling lines. Here we report on the analysis of a deep Herschel-PACS observation of an extremely high J CO transition (40-39) in the Seyfert 2 galaxy NGC 1068. The line was not detected, with a derived 3$\sigma$ upper limit of $2 \times 10^{-17}\,\text{W}\,\text{m}^{-2}$. We apply an XDR model in order to investigate whether the upper limit constrains the properties of a molecular torus in NGC 1068. The XDR model predicts the CO Spectral Line Energy Distributions for various gas densities and illuminating X-ray fluxes. In our model, the CO(40-39) upper limit is matched by gas with densities $\sim 10^{6}-10^{7}\,\text{cm}^{-3}$, located at $1.6-5\,\text{pc}$ from the AGN, with column densities of at least $10^{25}\,\text{cm}^{-2}$. At such high column densities, however, dust absorbs most of the CO(40-39) line emission at $\lambda = 65.69\, \mu$m. Therefore, even if NGC 1068 has a molecular torus which radiates in the CO(40-39) line, the dust can attenuate the line emission to below the PACS detection limit. The upper limit is thus consistent with the existence of a molecular torus in NGC 1068. In general, we expect that the CO(40-39) is observable in only a few AGN nuclei (if at all), because of the required high gas column density, and absorption by dust.Comment: 22 pages, accepted for publication in Ap

Excited OH+, H2O+, and H3O+ in NGC 4418 and Arp 220

We report on Herschel/PACS observations of absorption lines of OH+, H2O+ and H3O+ in NGC 4418 and Arp 220. Excited lines of OH+ and H2O+ with E_lower of at least 285 and \sim200 K, respectively, are detected in both sources, indicating radiative pumping and location in the high radiation density environment of the nuclear regions. Abundance ratios OH+/H2O+ of 1-2.5 are estimated in the nuclei of both sources. The inferred OH+ column and abundance relative to H nuclei are (0.5-1)x10^{16} cm-2 and \sim2x10^{-8}, respectively. Additionally, in Arp 220, an extended low excitation component around the nuclear region is found to have OH+/H2O+\sim5-10. H3O+ is detected in both sources with N(H3O+)\sim(0.5-2)x10^{16} cm-2, and in Arp 220 the pure inversion, metastable lines indicate a high rotational temperature of ~500 K, indicative of formation pumping and/or hot gas. Simple chemical models favor an ionization sequence dominated by H+ - O+ - OH+ - H2O+ - H3O+, and we also argue that the H+ production is most likely dominated by X-ray/cosmic ray ionization. The full set of observations and models leads us to propose that the molecular ions arise in a relatively low density (\gtrsim10^4 cm-3) interclump medium, in which case the ionization rate per H nucleus (including secondary ionizations) is zeta>10^{-13} s-1, a lower limit that is severalx10^2 times the highest rate estimates for Galactic regions. In Arp 220, our lower limit for zeta is compatible with estimates for the cosmic ray energy density inferred previously from the supernova rate and synchrotron radio emission, and also with the expected ionization rate produced by X-rays. In NGC 4418, we argue that X-ray ionization due to an AGN is responsible for the molecular ion production.Comment: 24 pages, 13 figures. Accepted for publication in Astronomy & Astrophysic

Complex hybridization patterns in European pond turtles (Emys orbicularis) in the Pyrenean Region

Hybrid zones are natural laboratories allowing insights in genetic processes like lineage diversification, speciation and introgression. Using large sampling, 15 microsatellite loci and a mitochondrial marker, we examined the Pyrenean contact zone of three pond turtle taxa (Emys orbicularis orbicularis, E. o. galloitalica, E. o. occidentalis). The Pyrenees are a biogeographically important region separating many lineages endemic to the Iberian Peninsula from their Western European counterparts. We discovered limited admixture, reflecting a complex biogeographic scenario. Simulations using Approximate Bayesian Computing supported that E. o. orbicularis invaded the Iberian Peninsula in the Holocene, circumventing the Pyrenees along the Mediterranean coast, and hybridized in the northern peninsula with the local coastal subspecies galloitalica, and to a lesser extent, with occidentalis. While E. o. occidentalis, and in particular E. o. orbicularis, expanded their ranges considerably during Holocene warming, E. o. galloitalica remained largely confined to its former Iberian refuge. Admixture among the three taxa is surprisingly low, and a future taxonomic investigation that includes the unstudied subspecies of E. orbicularis from North Africa, Eastern Europe and Asia has to determine whether their current status properly reflects their evolutionary divergence or whether certain taxa should be regarded as full species

Herschel PACS Spectroscopic Diagnostics of Local ULIRGs: Conditions and Kinematics in Mrk 231

In this first paper on the results of our Herschel PACS survey of local Ultraluminous Infrared Galaxies (ULIRGs), as part of our SHINING survey of local galaxies, we present far-infrared spectroscopy of Mrk 231, the most luminous of the local ULIRGs, and a type 1 broad absorption line AGN. For the first time in a ULIRG, all observed far-infrared fine-structure lines in the PACS range were detected and all were found to be deficient relative to the far infrared luminosity by 1 - 2 orders of magnitude compared with lower luminosity galaxies. The deficits are similar to those for the mid-infrared lines, with the most deficient lines showing high ionization potentials. Aged starbursts may account for part of the deficits, but partial covering of the highest excitation AGN powered regions may explain the remaining line deficits. A massive molecular outflow, discovered in OH and 18OH, showing outflow velocities out to at least 1400 km/sec, is a unique signature of the clearing out of the molecular disk that formed by dissipative collapse during the merger. The outflow is characterized by extremely high ratios of 18O / 16O suggestive of interstellar medium processing by advanced starbursts.Comment: Accepted for publication in the Astronomy and Astrophysics Herschel Special Issue, 5 pages, 4 figure

Molecular tracers of PDR-dominated galaxies

Photon-dominated regions (PDRs) are powerful molecular line emitters in external galaxies. They are expected in galaxies with high rates of massive star formation due to either starburst (SB) events or starburst coupled with active galactic nuclei (AGN) events. We have explored the PDR chemistry for a range of physical conditions representing a variety of galaxy types. Our main result is a demonstration of the sensitivity of the chemistry to changes in the physical conditions. We adopt crude estimates of relevant physical parameters for several galaxy types and use our models to predict suitable molecular tracers of those conditions. The set of recommended molecular tracers differs from that which we recommended for use in galaxies with embedded massive stars. Thus, molecular observations can in principle be used to distinguish between excitation by starburst and by SB+AGN in distant galaxies. Our recommendations are intended to be useful in preparing Herschel and ALMA proposals to identify sources of excitation in galaxies.Comment: 18 pages, 6 figures, Accepted in Ap

A Survey of HC3N in Extragalactic Sources - Is HC3N a Tracer of Activity in ULIRGs?

Context. HC3N is a molecule that is mainly associated with Galactic star-forming regions, but it has also been detected in extragalactic environments. Aims. To present the first extragalactic survey of HC3N, when combining earlier data from the literature with six new single-dish detections, and to compare HC3N with other molecular tracers (HCN, HNC), as well as other properties (silicate absorption strength, IR flux density ratios, C II flux, and megamaser activity). Methods. We present mm IRAM 30 m, OSO 20 m, and SEST observations of HC3N rotational lines (mainly the J = 10-9 transition) and of the J = 1-0 transitions of HCN and HNC. Our combined HC3N data account for 13 galaxies (excluding the upper limits reported for the non-detections), while we have HCN and HNC data for more than 20 galaxies. Results. A preliminary definition "HC3N-luminous galaxy" is made based upon the HC3N/HCN ratio. Most (~80 %) HC3N-luminous galaxies seem to be deeply obscured galaxies and (U)LIRGs. A majority (~60 % or more) of the HC3N-luminous galaxies in the sample present OH mega- or strong kilomaser activity. A possible explanation is that both HC3N and OH megamasers need warm dust for their excitation. Alternatively, the dust that excites the OH megamaser offers protection against UV destruction of HC3N. A high silicate absorption strength is also found in several of the HC3N-luminous objects, which may help the HC3N to survive. Finally, we find that a high HC3N/HCN ratio is related to a high dust temperature and a low C II flux.Comment: 17 pages, 12 figures, to be published in Astronomy & Astrophysic

Herschel-PACS Observations of Far-IR CO Line Emission in NGC 1068: Highly Excited Molecular Gas in the Circumnuclear Disk

We report the detection of far-IR CO rotational emission from the prototypical Seyfert 2 galaxy NGC 1068. Using Herschel-PACS, we have detected 11 transitions in the J_upper=14-30 (E_upper/k_B = 580-2565 K) range, all of which are consistent with arising from within the central 10" (700 pc). The detected transitions are modeled as arising from 2 different components: a moderate excitation (ME) component close to the galaxy systemic velocity, and a high excitation (HE) component that is blueshifted by ~80 km s^{-1}. We employ a large velocity gradient (LVG) model and derive n_H2~10^{5.6} cm^{-3}, T_kin~170 K, and M_H2~10^{6.7} M_sun for the ME component, and n_H2~10^{6.4} cm^{-3}, T_kin~570 K, and M_H2~10^{5.6} M_sun for the HE component, although for both components the uncertainties in the density and mass are plus/minus (0.6-0.9) dex. We compare the CO line profiles with those of other molecular tracers observed at higher spatial and spectral resolution, and find that the ME transitions are consistent with these lines arising in the ~200 pc diameter ring of material traced by H_2 1-0 S(1) observations. The blueshift of the HE lines may also be consistent with the bluest regions of this H_2 ring, but a better kinematic match is found with a clump of infalling gas ~40 pc north of the AGN. We discuss the prospects of placing the HE component near the AGN, and conclude that while the moderate thermal pressure precludes an association with the ~1 pc radius H_2O maser disk, the HE component could potentially be located only a few parsecs more distant from the AGN, and might then provide the N_H~10^{25} cm^{-2} column obscuring the nuclear hard X-rays. Finally, we also report sensitive upper limits extending up to J_upper=50, which place constraints on a previous model prediction for the CO emission from the X-ray obscuring torus. [Abridged]Comment: 20 pages, 10 figures. Accepted for publication in Ap