Outflowing OH$^+$ in Markarian 231: the ionization rate of the molecular gas

The oxygen-bearing molecular ions OH+, H2O+, and H3O+ are key species that probe the ionization rate of (partially) molecular gas that is ionized by X-rays and cosmic rays permeating the interstellar medium. We report Herschel far-infrared and submillimeter spectroscopic observations of OH+ in Mrk 231, showing both ground-state P-Cygni profiles, and excited line profiles with blueshifted absorption wings extending up to ~1000 km s^{-1}. In addition, OH+ probes an excited component peaking at central velocities, likely arising from the torus probed by the OH centimeter-wave megamaser. Four lines of H2O+ are also detected at systemic velocities, but H3O+ is undetected. Based on our earlier OH studies, we estimate an abundance ratio of OH/OH+~5-10 for the outflowing components and ~20 for the torus, and an OH+ abundance relative to H nuclei of ~>10^{-7}. We also find high OH+/H2O+ and OH+/H3O+ ratios, both are ~>4 in the torus and ~>10-20 in the outflowing gas components. Chemical models indicate that these high OH+ abundances relative to OH, H2O+, and H3O+ are characteristic of gas with a high ionization rate per unit density, zeta/n_H~(1-5)x10^{-17} cm^3 s^{-1} and ~(1-2)x10^{-16} cm^3 s^{-1} for the above components, respectively, and an ionization rate of zeta~(0.5-2)x10^{-12} s^{-1}. X-rays appear to be unable to explain the inferred ionization rate, and thus we suggest that low-energy (10-400 MeV) cosmic-rays are primarily responsible for the ionization with dot{M}_{CR}~0.01 M_{sun} yr^{-1} and dot{E}_{CR}~10^{44} erg s^{-1}, the latter corresponding to 1% of the AGN luminosity and similar to the energetics of the molecular outflow. We suggest that cosmic-rays accelerated in the forward shock associated with the molecular outflow are responsible for the ionization, as they diffuse through the outflowing molecular phase downstream.

Volatile carbon locking and release in protoplanetary disks. A study of TW Hya and HD 100546

The composition of planetary solids and gases is largely rooted in the processing of volatile elements in protoplanetary disks. To shed light on the key processes, we carry out a comparative analysis of the gas-phase carbon abundance in two systems with a similar age and disk mass, but different central stars: HD 100546 and TW Hya. We combine our recent detections of C$^{0}$ in these disks with observations of other carbon reservoirs (CO, C$^{+}$, C$_{2}$H) and gas mass and warm gas tracers (HD, O$^{0}$), as well as spatially resolved ALMA observations and the spectral energy distribution. The disks are modelled with the DALI 2D physical-chemical code. Stellar abundances for HD 100546 are derived from archival spectra. Upper limits on HD emission from HD 100546 place an upper limit on the total disk mass of $leq0.1,M_{odot}$. The gas-phase carbon abundance in the atmosphere of this warm Herbig disk is at most moderately depleted compared to the interstellar medium, with [C]/[H]$_{rm gas}=(0.1-1.5)times 10^{-4}$. HD 100546 itself is a $lambda,$Bo{o}tis star, with solar abundances of C and O but a strong depletion of rock-forming elements. In the gas of the T Tauri disk TW Hya, both C and O are strongly underabundant, with [C]/[H]$_{rm gas}=(0.2-5.0)times 10^{-6}$ and C/O${>}1$. We discuss evidence that the gas-phase C and O abundances are high in the warm inner regions of both disks. Our analytical model, including vertical mixing and a grain size distribution, reproduces the observed [C]/[H]$_{rm gas}$ in the outer disk of TW Hya and allows to make predictions for other systems.

WISH VI. Constraints on UV and X-ray irradiation from a survey of hydrides in low- to high-mass YSOs

Hydrides are simple compounds containing one or a few hydrogen atoms bonded to a heavier atom. They are fundamental precursor molecules in cosmic chemistry and many hydride ions have become observable in high quality for the first time thanks to the Herschel Space Observatory. Ionized hydrides, such as CH+ and OH+, and also HCO+ that affect the chemistry of molecules such as water, provide complementary information on irradiation by far UV (FUV) or X-rays and gas temperature. The targeted lines of CH+, OH+, H2O+, C+ and CH are detected mostly in blue-shifted absorption. H3O+ and SH+ are detected in emission and only toward some high-mass objects. The observed line parameters and correlations suggest two different origins, related to gas entrained by the outflows and to the circumstellar envelope. The column density ratios of CH+/OH+ are estimated from chemical slab models, assuming that the H2 density is given by the specific density model of each object at the beam radius. For the low-mass YSOs the observed ratio can be reproduced for an FUV flux of 2-400 times the ISRF at the location of the molecules. In two high-mass objects, the UV flux is 20-200 times the ISRF derived from absorption lines, and 300-600 ISRF using emission lines. If the FUV flux required for low-mass objects originates at the central protostar, a substantial FUV luminosity, up to 1.5 L_sun, is required. There is no molecular evidence for X-ray induced chemistry in the low-mass objects on the observed scales of a few 1000 AU. For high-mass regions, the FUV flux required to produce the observed molecular ratios is smaller than the unattenuated flux expected from the central object(s) at the Herschel beam radius. This is consistent with an FUV flux reduced by circumstellar extinction or by bloating of the protostar.

Observations and modelling of CO and [CI] in disks. First detections of [CI] and constraints on the carbon abundance

The gas-solid budget of carbon in protoplanetary disks is related to the composition of the cores and atmospheres of the planets forming in them. The key gas-phase carbon carriers CO, C$^{0}$ and C$^{+}$ can now be observed in disks. The gas-phase carbon abundance in disks has not yet been well characterized, we aim to obtain new constraints on the [C]/[H] ratio in a sample of disks, and to get an overview of the strength of [CI] and warm CO emission. We carried out a survey of the CO$,6$--$5$ and [CI]$,1$--$0$ and $2$--$1$ lines towards $37$ disks with APEX, and supplemented it with [CII] data from the literature. The data are interpreted using a grid of models produced with the DALI code. We also investigate how well the gas-phase carbon abundance can be determined in light of parameter uncertainties. The CO$,6$--$5$ line is detected in $13$ out of $33$ sources, the [CI]$,1$--$0$ in $6$ out of $12$, and the [CI]$,2$--$1$ in $1$ out of $33$. With deep integrations, the first unambiguous detections of [CI]~$1$--$0$ in disks are obtained, in TW~Hya and HD~100546. Gas-phase carbon abundance reductions of a factor $5$--$10$ or more can be identified robustly based on CO and [CI] detections. The atomic carbon detection in TW~Hya confirms a factor $100$ reduction of [C]/[H]$_{rm gas}$ in that disk, while the data are consistent with an ISM-like carbon abundance for HD~100546. In addition, BP~Tau, T~Cha, HD~139614, HD~141569, and HD~100453 are either carbon-depleted or gas-poor disks. The low [CI]~$2$--$1$ detection rates in the survey mostly reflect insufficient sensitivity to detect T~Tauri disks. The Herbig~Ae/Be disks with CO and [CII] upper limits below the models are debris disk like systems. A roughly order of magnitude increase in sensitivity compared to our survey is required to obtain useful constraints on the gas-phase [C]/[H] ratio in most of the targeted systems.

On the asymmetry of the OH ro-vibrational lines in HD 100546

We present multi-epoch high-spectral resolution observations with VLT/CRIRES of the OH doublet $^2Pi_{3/2}$ P4.5 (1+,1-) (2.934 $mu$m) towards the protoplanetary disk around HD 100546. The OH doublet is detected at all epochs and is spectrally resolved while nearby H$_2$O lines remains undetected. The OH line velocity profile is different in the three datasets: in the first epoch (April 2012, PA=26$^{circ}$) the OH lines are symmetric and line broadening is consistent with the gas being in Keplerian rotation around the star. No OH emission is detected within a radius of 8-11 au from the star: the line emitting region is similar in size and extent to that of the CO ro-vibrational lines. In the other two epochs (March 2013 and April 2014, PA=90$^{circ}$ and 10$^{circ}$, respectively) the OH lines appear asymmetric and fainter compared to April 2012. We investigate the origin of these line asymmetries which were taken by previous authors as evidence for tidal interaction between an (unseen) massive planet and the disk. We show that the observed asymmetries can be fully explained by a misalignment of the slit of order 0farcs04-0farcs20 with respect to the stellar position. The disk is spatially resolved and the slit misalignment is likely caused by the extended dust emission which is brighter than the stellar photosphere at near-infrared wavelengths which is the wavelength used for the pointing. This can cause the photo-center of HD 100546 to be mis-aligned with the stellar position at near-infrared wavelengths.

Far-infrared molecular lines from Low- to High-Mass Star Forming Regions observed with Herschel

(Abridged) We study the response of the gas to energetic processes associated with high-mass star formation and compare it with studies on low- and intermediate-mass young stellar objects (YSOs) using the same methods. The far-IR line emission and absorption of CO, H$_2$O, OH, and [OI] reveals the excitation and the relative contribution of different species to the gas cooling budget. Herschel-PACS spectra covering 55-190 um are analyzed for ten high-mass star forming regions of various luminosities and evolutionary stages at spatial scales of ~10^4 AU. Radiative transfer models are used to determine the contribution of the envelope to the far-IR CO emission. The close environments of high-mass YSOs show strong far-IR emission from molecules, atoms, and ions. Water is detected in all 10 objects even up to high excitation lines. CO lines from J=14-13 up to typically 29-28 show a single temperature component, Trot~300 K. Typical H$_2$O temperatures are Trot~250 K, while OH has Trot~80 K. Far-IR line cooling is dominated by CO (~75 %) and to a smaller extent by OI (~20 %), which increases for the most evolved sources. H$_2$O is less important as a coolant for high-mass sources because many lines are in absorption. Emission from the envelope is responsible for ~45-85 % of the total CO luminosity in high-mass sources compared with only ~10 % for low-mass YSOs. The highest-J lines originate most likely from shocks, based on the strong correlation of CO and H$_2$O with physical parameters of the sources from low- to high-masses. Excitation of warm CO is very similar for all mass regimes, whereas H$_2$O temperatures are ~100 K higher for high-mass sources than the low-mass YSOs. Molecular cooling is ~4 times more important than cooling by [OI]. The total far-IR line luminosity is about 10$^{-3}$ and 10$^{-5}$ times lower than the dust luminosity for the low- and high-mass YSOs.

Probing the radial temperature structure of protoplanetary disks with Herschel/HIFI

Herschel/HIFI spectroscopic observations of CO J=10-9, CO J=16-15 and [CII] towards HD 100546 are presented. The objective is to resolve the velocity profile of the lines to address the emitting region of the transitions and directly probe the distribution of warm gas in the disk. The spectra reveal double-peaked CO line profiles centered on the systemic velocity, consistent with a disk origin. The J=16-15 line profile is broader than that of the J=10-9 line, which in turn is broader than those of lower J transitions (6-5, 3-2, observed with APEX), thus showing a clear temperature gradient of the gas with radius. A power-law flat disk model is used to fit the CO line profiles and the CO rotational ladder simultaneously, yielding a temperature of T_0=1100 pm 350 K (at r_0 = 13 AU) and an index of q=0.85 pm 0.1 for the temperature radial gradient. This indicates that the gas has a steeper radial temperature gradient than the dust (mean q_{dust} ~ 0.5), providing further proof of the thermal decoupling of gas and dust at the disk heights where the CO lines form. The [CII] line profile shows a strong single-peaked profile red-shifted by 0.5 km s-1 compared to the systemic velocity. We conclude that the bulk of the [CII] emission has a non-disk origin (e.g., remnant envelope or diffuse cloud).

DIGIT survey of far-infrared lines from protoplanetary disks I

[abridged] We present far-infrared spectroscopic observations of PMS stars taken with Herschel/PACS as part of the DIGIT key project. The sample includes 22 Herbig AeBe and 8 T Tauri sources. Multiple atomic fine structure and molecular lines are detected at the source position: [OI], [CII], CO, OH, H_2O, CH^+. The most common feature is the [OI] 63micron line detected in almost all of the sources followed by OH. In contrast with CO, OH is detected toward both Herbig AeBe groups (flared and non-flared sources). An isothermal LTE slab model fit to the OH lines indicates column densities of 10^13 < N_OH < 10^16 cm^-2, emitting radii 15 < r < 100 AU and excitation temperatures 100 < T_ex < 400 K. The OH emission thus comes from a warm layer in the disk at intermediate stellar distances. Warm H_2O emission is detected through multiple lines toward the T Tauri systems AS 205, DG Tau, S CrA and RNO 90 and three Herbig AeBe systems HD 104237, HD 142527, HD 163296 (through line stacking). Overall, Herbig AeBe sources have higher OH/H_2O abundance ratios across the disk than do T Tauri disks, from near- to far-infrared wavelengths. Far-infrared CH^+ emission is detected toward HD 100546 and HD 97048. The slab model suggests moderate excitation (T_ex ~ 100 K) and compact (r ~ 60 AU) emission in the case of HD 100546. The [CII] emission is spatially extended in all sources where the line is detected. This suggests that not all [CII] emission is associated with the disk and that there is a substantial contribution from diffuse material around the young stars. The flux ratios of the atomic fine structure lines are consistent with a disk origin for the oxygen lines for most of the sources.

OH far-infrared emission from low- and intermediate-mass protostars surveyed with Herschel-PACS

OH is a key species in the water chemistry of star-forming regions, because its presence is tightly related to the formation and destruction of water. This paper presents OH observations from 23 low- and intermediate-mass young stellar objects obtained with the PACS integral field spectrometer on-board Herschel in the context of the Water In Star-forming Regions with Herschel (WISH) key program. Most low-mass sources have compact OH emission (< 5000 AU scale), whereas the OH lines in most intermediate-mass sources are extended over the whole PACS detector field-of-view (> 20000 AU). The strength of the OH emission is correlated with various source properties such as the bolometric luminosity and the envelope mass, but also with the OI and H2O emission. Rotational diagrams for sources with many OH lines show that the level populations of OH can be approximated by a Boltzmann distribution with an excitation temperature at around 70 K. Radiative transfer models of spherically symmetric envelopes cannot reproduce the OH emission fluxes nor their broad line widths, strongly suggesting an outflow origin. Slab excitation models indicate that the observed excitation temperature can either be reached if the OH molecules are exposed to a strong far-infrared continuum radiation field or if the gas temperature and density are sufficiently high. Using realistic source parameters and radiation fields, it is shown for the case of Ser SMM1 that radiative pumping plays an important role in transitions arising from upper level energies higher than 300 K. The compact emission in the low-mass sources and the required presence of a strong radiation field and/or a high density to excite the OH molecules points towards an origin in shocks in the inner envelope close to the protostar.

Warm H2O and OH in the disk around the Herbig star HD 163296

We present observations of far-infrared (50-200 micron) OH and H2O emission of the disk around the Herbig Ae star HD 163296 obtained with Herschel/PACS in the context of the DIGIT key program. In addition to strong [OI] emission, a number of OH doublets and a few weak highly excited lines of H2O are detected. The presence of warm H2O in this Herbig disk is confirmed by a line stacking analysis, enabled by the full PACS spectral scan, and by lines seen in Spitzer data. The line fluxes are analyzed using an LTE slab model including line opacity. The water column density is 10^14 - 10^15 cm^-2, and the excitation temperature is 200-300 K implying warm gas with a density n > 10^5 cm^-3. For OH we find a column density of 10^14 - 2x10^15 cm^-2 and T_ex ~ 300-500 K. For both species we find an emitting region of r ~ 15-20 AU from the star. We argue that the molecular emission arises from the protoplanetary disk rather than from an outflow. This far-infrared detection of both H2O and OH contrasts with near- and mid-infrared observations, which have generally found a lack of water in the inner disk around Herbig AeBe stars due to strong photodissociation of water. Given the similarity in column density and emitting region, OH and H2O emission seems to arise from an upper layer of the disk atmosphere of HD 163296, probing a new reservoir of water. The slightly lower temperature of H2O compared to OH suggests a vertical stratification of the molecular gas with OH located higher and water deeper in the disk, consistent with thermo-chemical models.