M. Spaans

Diagnostics of irradiated dense gas in galaxy nuclei - II. A grid of XDR and PDR models

Aims. The nuclei of active galaxies harbor massive young stars, an accreting central black hole, or both. In order to determine the physical conditions that pertain to molecular gas close to t he sources of radiation, numerical models are constructed. Methods. These models iteratively determine the thermal and chemical balance of molecular gas that is exposed to X-rays (1-100 keV) and far-ultraviolet radiation (6-13.6 eV), as a functi on of depth. Results. We present a grid of XDR and PDR models that span ranges in density (10 2 10 6.5 cm 3 ), irradiation (10 0.5 10 5 G0 and FX = 1.6 × 10 2 160 erg cm 2 s 1 ) and column density (3 × 10 21 1 × 10 25 cm 2 ). Predictions are made for the most important atomic fine-structure lines, e.g., [CII], [OI], [ CI], [SiII], and for molecular species like HCO + , HCN, HNC, CS and SiO up to J = 4, CO and 13 CO up to J = 16, and column densities for CN, CH, CH + , HCO, HOC + , NO and N2H + . We find that surface temperatures are higher (lower) in PDRs compared to XDRs for densities > 10 4 ( 1) for XDRs (PDRs) if the density exceeds 10 5 cm 3 and if the column density is larger than 10 23 cm 2 . For columns less than 10 22.5 cm 2 the XDR HCN/HCO + 1-0 ratio becomes larger than one, although the individual HCN 1-0 and HCO + 1-0 line intensities are weaker. For modest densities, n = 10 4 10 5 cm 3 , and strong radiation fields ( > 100 erg s 1 cm 2 ), HCN/HCO + ratios can become larger in XDRs than PDRs as well. Also, the HCN/CO 1-0 ratio is typically smaller in XDRs, and the HCN emission in XDRs is boosted with respect to CO only for high (column) density gas, with columns in excess of 10 23 cm 2 and densities larger than 10 4 cm 3 . Furthermore, CO is typically warmer in XDRs than in PDRs, for the same total energy input. This leads to higher CO J=N+1-N/CO 1-0, N � 1, line ratios in XDRs. In particular, lines with N � 10, like CO(16-15) and CO(10-9) observable with HIFI/Herschel, discriminate very well between XDRs and PDRs. This is crucial since the XDR/AGN contribution will typically be of a much smaller (possibly beam diluted) angular scale and a 10-25% PDR contribution can already suppress XDR distinguishing features involving HCN/HCO+ and HNC/HCN. For possible future observations, column density ratios indicate that CH, CH + , NO, HOC + and HCO

Molecular line emission in NGC 1068 imaged with ALMA - I. An AGN-driven outflow in the dense molecular gas

Aims. We investigate the fueling and the feedback of star formation and nuclear activity in NGC 1068, a nearby (D = 14 Mpc) Seyfert 2 barred galaxy, by analyzing the distribution and kinematics of the molecular gas in the disk. We aim to understand if and how gas accretion can self-regulate.Methods. We have used the Atacama Large Millimeter Array (ALMA) to map the emission of a set of dense molecular gas (n(H2) 1056 cm3) tracers (CO(3-2), CO(6-5), HCN(4-3), HCO+(4-3), and CS(7-6)) and their underlying continuum emission in the central r ∼ 2 kpc of NGC 1068 with spatial resolutions ∼0:3000:500 (∼20-35 pc for the assumed distance of D = 14 Mpc). Results. The sensitivity and spatial resolution of ALMA give an unprecedented detailed view of the distribution and kinematics of the dense molecular gas (n(H2) ≈ 1056cm3) in NGC 1068. Molecular line and dust continuum emissions are detected from a r ∼ 200 pc off-centered circumnuclear disk (CND), from the 2.6 kpc-diameter bar region, and from the r ∼ 1:3 kpc starburst (SB) ring. Most of the emission in HCO+, HCN, and CS stems from the CND. Molecular line ratios show dramatic order-of-magnitude changes inside the CND that are correlated with the UV/X-ray illumination by the active galactic nucleus (AGN), betraying ongoing feedback. We used the dust continuum fluxes measured by ALMA together with NIR/MIR data to constrain the properties of the putative torus using CLUMPY models and found a torus radius of 20+6 10 pc. The Fourier decomposition of the gas velocity field indicates that rotation is perturbed by an inward radial flow in the SB ring and the bar region. However, the gas kinematics from r ∼ 50 pc out to r ∼ 400 pc reveal a massive (Mmol ∼ 2:7+0:9 1:2 × 107 M) outflow in all molecular tracers. The tight correlation between the ionized gas outflow, the radio jet, and the occurrence of outward motions in the disk suggests that the outflow is AGN driven. Conclusions. The molecular outflow is likely launched when the ionization cone of the narrow line region sweeps the nuclear disk. The outflow rate estimated in the CND, dM=dt ∼ 63+21 37 M yr1, is an order of magnitude higher than the star formation rate at these radii, confirming that the outflow is AGN driven. The power of the AGN is able to account for the estimated momentum and kinetic luminosity of the outflow. The CND mass load rate of the CND outflow implies a very short gas depletion timescale of ≤1 Myr. The CND gas reservoir is likely replenished on longer timescales by efficient gas inflow from the outer disk.

Star formation in extreme environments: The effects of cosmic rays and mechanical heating

Context. The molecular interstellar medium in extreme environments, such as Arp 220, but also NGC 253 appears to have extremely high cosmic ray (CR) rates (10 3 −10 4 × Milky Way) and substantial mechanical heating from supernova driven turbulence. Aims. We explore the consequences of high CR rates and mechanical heating on the chemistry of the clouds. Methods. PDR model predictions are made for low, n = 10 3 , and high, n = 10 5.5 cm −3 , density clouds using well-tested chemistry and radiation transfer codes. Column densities of relevant species are discussed, and special attention is given to water-related species. Fluxes are shown for fine-structure lines of O, C + ,C , and N + , and molecular lines of CO, HCN, HNC, and HCO + . A comparison is made with an X-ray dominated region model. Results. Fine-structure lines of [CII], [CI], and [OI] are remarkably similar for different mechanical heating and CR rates, when already exposed to large amounts of UV. Both HCN and H2O abundances are boosted for very high mechanical heating rates, while ionized species are relatively unaffected. Both OH + and H2O + are enhanced for very high CR rates ζ ≥ 5 × 10 −14 s −1 . A combination of OH + , OH, H2O + ,H 2O, and H3O + traces the CR rates, and is able to distinguish between enhanced cosmic rays and X-rays.

Evidence for CO Shock Excitation in NGC 6240 from Herschel SPIRE Spectroscopy

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 [C I] and [N II] 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 H-2 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 x 10(4) cm(-3). We find that the solution best fitting the H-2 lines is degenerate. The shock velocities and number densities range between v(s) = 17-47 km s(-1) and n(H) = 10(7)-5x10(4) cm(-3), respectively. The H-2 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 (less than or similar to 1%) of the interstellar medium is exposed to the high-velocity shocks. This implies that the gas is rapidly losing its highly turbulent motions. We argue that a high CO line-to-continuum ratio is a key diagnostic for the presence of shocks.

Consistent dust and gas models for protoplanetary disks I. Disk shape, dust settling, opacities, and PAHs

We propose a set of standard assumptions for the modelling of Class II and III protoplanetary disks, which includes detailed continuum radiative transfer, thermo-chemical modelling of gas and ice, and line radiative transfer from optical to cm wavelengths. The first paper of this series focuses on the assumptions about the shape of the disk, the dust opacities, dust settling, and polycyclic aromatic hydrocarbons (PAHs). In particular, we propose new standard dust opacities for disk models, we present a simplified treatment of PAHs in radiative equilibrium which is sufficient to reproduce the PAH emission features, and we suggest using a simple yet physically justified treatment of dust settling. We roughly adjust parameters to obtain a model that predicts continuum and line observations that resemble typical multi-wavelength continuum and line observations of Class II T Tauri stars. We systematically study the impact of each model parameter (disk mass, disk extension and shape, dust settling, dust size and opacity, gas/dust ratio, etc.) on all mainstream continuum and line observables, in particular on the SED, mm-slope, continuum visibilities, and emission lines including [OI] 63 μm, high-J CO lines, (sub-)mm CO isotopologue lines, and CO fundamental ro-vibrational lines. We find that evolved dust properties, i.e. large grains, often needed to fit the SED, have important consequences for disk chemistry and heating/cooling balance, leading to stronger near- to far-IR emission lines in general. Strong dust settling and missing disk flaring have similar effects on continuum observations, but opposite effects on far-IR gas emission lines. PAH molecules can efficiently shield the gas from stellar UV radiation because of their strong absorption and negligible scattering opacities in comparison to evolved dust. The observable millimetre-slope of the SED can become significantly more gentle in the case of cold disk midplanes, which we find regularly in our T Tauri models. We propose to use line observations of robust chemical tracers of the gas, such as O, CO, and H2, as additional constraints to determine a number of key properties of the disks, such as disk shape and mass, opacities, and the dust/gas ratio, by simultaneously fitting continuum and line observations.

The SCUBA-2 Cosmology Legacy Survey: ALMA Resolves the Bright-end of the Sub-millimeter Number Counts

We present high-resolution 870 μm Atacama Large Millimeter/sub-millimeter Array (ALMA) continuum maps of 30 bright sub-millimeter sources in the UKIDSS UDS field. These sources are selected from deep, 1 degree2 850 μm maps from the SCUBA-2 Cosmology Legacy Survey, and are representative of the brightest sources in the field (median SSCUBA-2= 8.7 ± 0.4 mJy). We detect 52 sub-millimeter galaxies (SMGs) at >4σ significance in our 30 ALMA maps. In 61-15+19% of the ALMA maps the single-dish source comprises a blend of ≥2 SMGs, where the secondary SMGs are Ultra-luminous Infrared Galaxies (ULIRGs) with LIR ≳ 1012 L⊙. The brightest SMG contributes on average 80-2+6% of the single-dish flux density, and in the ALMA maps containing ≥2 SMGs the secondary SMG contributes 25-5+1% of the integrated ALMA flux. We construct source counts and show that multiplicity boosts the apparent single-dish cumulative counts by 20% at S870 > 7.5 mJy, and by 60% at S870 > 12 mJy. We combine our sample with previous ALMA studies of fainter SMGs and show that the counts are well-described by a double power law with a break at 8.5 ± 0.6 mJy. The break corresponds to a luminosity of ∼6 × 1012 L⊙ or a star formation rate (SFR) of ∼103 M⊙ yr-1. For the typical sizes of these SMGs, which are resolved in our ALMA data with Re = 1.2 ± 0.1 kpc, this yields a limiting SFR density of ∼100 M⊙ yr-1 kpc-2 Finally, the number density of S870 ≳ 2 mJy SMGs is 80 ± 30 times higher than that derived from blank-field counts. An over-abundance of faint SMGs is inconsistent with line-of-sight projections dominating multiplicity in the brightest SMGs, and indicates that a significant proportion of these high-redshift ULIRGs are likely to be physically associated.

Irradiated ISM: Discriminating between Cosmic Rays and X-Rays

The interstellar medium (ISM) at the centers of active galaxies is exposed to a combination of cosmic-ray, far-ultraviolet (FUV), and X-ray radiation. We apply photodissociation region (PDR) models to this ISM with both normal and highly elevated (5 × 10-15 s-1) cosmic-ray (CR) rates and compare the results to those obtained for X-ray dissociation regions (XDRs). Our existing PDR-XDR code is used to construct models over a 103-105 cm-3 density range and for 0.16-160 ergs s-1 cm-2 impingent fluxes. We obtain larger high-J (J > 10) CO ratios in PDRs when we use the highly elevated CR rate, but these are always exceeded by the corresponding XDR ratios. The [C I] 609 μm/13CO (2-1) line ratio is boosted by a factor of a few in PDRs with n ~ 103 cm-3 exposed to a high CR rate. At higher densities, ratios become identical irrespective of CR flux, while XDRs always show elevated [C I] emission per CO column. The HCN/CO and HCN/HCO+ line ratios, combined with high-J CO emission lines, are good diagnostics to distinguish between PDRs, under either low or high CR irradiation conditions, and XDRs. Hence, the Heterodyne Instrument for the Far Infrared (HIFI) on the Herschel Space Observatory, which can detect these CO lines, will be crucial in the study of active galaxies.

Mechanical feedback in the molecular ISM of luminous IR galaxies

Aims. Molecular emission lines originating in the nuclei of luminous infra-red galaxies are used to determine the physical properties of the nuclear ISM in these systems. Methods. A large observational database of molecular emission lines is compared with model predictions that include heating by UV and X-ray radiation, mechanical heating, and the effects of cosmic rays. Results. The observed line ratios and model predictions imply a separation of the observed systems into three groups: XDRs, UV-dominated high-density (n >= 10(5) cm(-3)) PDRs, and lower-density ( n = 10(4.5) cm-3) PDRs that are dominated by mechanical feedback. Conclusions. The division of the two types of PDRs follows naturally from the evolution of the star formation cycle of these sources, which evolves from deeply embedded young stars, resulting in high-density ( n >= 10(5) cm(-3)) PDRs, to a stage where the gas density has decreased ( n = 10(4.5) cm(-3)) and mechanical feedback from supernova shocks dominates the heating budget.

Star Formation Relations and CO Spectral Line Energy Distributions across the J-ladder and Redshift

We present FIR [50-300 mu m]-CO luminosity relations (i.e., log L-FIR = alpha log L(CO) + beta) for the full CO rotational ladder from J = 1-0 up to J = 13-12 for a sample of 62 local (z 10(11) L-circle dot) using data from Herschel SPIRE-FTS and ground-based telescopes. We extend our sample to high redshifts (z > 1) by including 35 submillimeter selected dusty star forming galaxies from the literature with robust CO observations, and sufficiently well-sampled FIR/submillimeter spectral energy distributions (SEDs), so that accurate FIR luminosities can be determined. The addition of luminous starbursts at high redshifts enlarge the range of the FIR-CO luminosity relations toward the high-IR-luminosity end, while also significantly increasing the small amount of mid-J/high-J CO line data (J = 5-4 and higher) that was available prior to Herschel. This new data set (both in terms of IR luminosity and J-ladder) reveals linear FIR-CO luminosity relations (i.e., a similar or equal to 1) for J = 1-0 up to J = 5-4, with a nearly constant normalization (beta similar to 2). In the simplest physical scenario, this is expected from the (also) linear FIR-(molecular line) relations recently 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 the J = 13-12 transition, we find an increasingly sublinear slope and higher normalization constant with increasing J. We argue that these are caused by a warm (similar to 100 K) and dense (>10(4) 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, 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.

The Herschel Comprehensive (U)LIRG Emission Survey (HERCULES): CO Ladders, Fine Structure Lines, and Neutral Gas Cooling

(Ultra) luminous infrared galaxies ((U)LIRGs) are objects characterized by their extreme infrared (8-1000 mu m) luminosities (L-LIRG > 10(11) L-circle dot and L-ULIRG > 10(12) L-circle dot). The Herschel Comprehensive ULIRG Emission Survey (PI: van derWerf) presents a representative flux-limited sample of 29 (U)LIRGs that spans the full luminosity range of these objects (10(11)L(circle dot) <= L-IR <= 10(13)L(circle dot)). With the Herschel Space Observatory, we observe [CII] 157 mu m, [O I] 63 mu m, and [O I] 145 mu m line emission with Photodetector Array Camera and Spectrometer, CO J = 4-3 through J = 13-12, [C I] 370 mu m, and [C I] 609 mu m with SPIRE, and low-J CO transitions with ground-based telescopes. The CO ladders of the sample are separated into three classes based on their excitation level. In 13 of the galaxies, the [O I] 63 mu m emission line is self absorbed. Comparing the CO excitation to the InfraRed Astronomical Satellite 60/100 mu m ratio and to far infrared luminosity, we find that the CO excitation is more correlated to the far infrared colors. We present cooling budgets for the galaxies and find fine-structure line flux deficits in the [C II], [Si II], [O I], and [C I] lines in the objects with the highest far IR fluxes, but do not observe this for CO 4 <= J(upp) <= 13. In order to study the heating of the molecular gas, we present a combination of three diagnostic quantities to help determine the dominant heating source. Using the CO excitation, the CO J = 1-0 linewidth, and the active galactic nucleus (AGN) contribution, we conclude that galaxies with large CO linewidths always have high-excitation CO ladders, and often low AGN contributions, suggesting that mechanical heating is important.