sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Meijerink, R.

Normalized to: Meijerink, R.

76 article(s) in total. 381 co-authors, from 1 to 44 common article(s). Median position in authors list is 4,0.

[1]  [pdf] - 1622341
X-ray radiative transfer in protoplanetary disks - The role of dust and X-ray background fields
Comments: 16 pages, 12 figures, accepted for publication in A&A
Submitted: 2017-11-20
The X-ray luminosities of T Tauri stars are about two to four orders of magnitude higher than the luminosity of the contemporary Sun. As these stars are born in clusters, their disks are not only irradiated by their parent star but also by an X-ray background field produced by the cluster members. We aim to quantify the impact of X-ray background fields produced by young embedded clusters on the chemical structure of disks. Further, we want to investigate the importance of the dust for X-ray radiative transfer in disks. We present a new X-ray radiative transfer module for the radiation thermo-chemical disk code ProDiMo, which includes X-ray scattering and absorption by both the gas and dust component. For the X-ray radiative transfer, we consider irradiation by the star and by X-ray background fields. To study the impact of X-rays on the chemical structure of disks we use the well-established disk ionization tracers N2H+ and HCO+. For evolved dust populations, X-ray opacities are mostly dominated by the gas; only for photon energies $E\gtrsim5-10\,$keV, dust opacities become relevant. Consequently, the local disk X-ray radiation field is only affected in dense regions close to the disk midplane. X-ray background fields can dominate the local X-ray disk ionization rate for disk radii $r\gtrsim20\,$au. However, the N2H+ and HCO+ column densities are only significantly affected in case of low cosmic-ray ionization rates, or if the background flux is at least a factor of ten higher than the flux level expected for clusters typical for the solar vicinity. Observable signatures of X-ray background fields in low-mass star-formation regions, like Taurus, are only expected for cluster members experiencing a strong X-ray background field. For the majority of the cluster members, the X-ray background field has only little impact on the disk chemical structure.
[2]  [pdf] - 1598093
Molecular gas in three z~7 quasar host galaxies
Comments: 12 pages, 6 figures. Accepted for publication in ApJ
Submitted: 2017-07-17
We present ALMA band 3 observations of the CO(6-5), CO(7-6), and [CI] 369micron emission lines in three of the highest redshift quasar host galaxies at 6.6<z<6.9. These measurements constitute the highest-redshift CO detections to date. The target quasars have previously been detected in [CII] 158micron emission and the underlying far-infrared (FIR) dust continuum. We detect (spatially unresolved, at a resolution of >2", or >14kpc) CO emission in all three quasar hosts. In two sources, we detect the continuum emission around 400micron (rest-frame), and in one source we detect [CI] at low significance. We derive molecular gas reservoirs of (1-3)x10^10 M_sun in the quasar hosts, i.e. approximately only 10 times the mass of their central supermassive black holes. The extrapolated [CII]-to-CO(1-0) luminosity ratio is 2500-4200, consistent with measurements in galaxies at lower redshift. The detection of the [CI] line in one quasar host galaxy and the limit on the [CI] emission in the other two hosts enables a first characterization of the physical properties of the interstellar medium in z~7 quasar hosts. In the sources, the derived global CO/[CII]/[CI] line ratios are consistent with expectations from photodissociation regions (PDR), but not X-ray dominated regions (XDR). This suggest that quantities derived from the molecular gas and dust emission are related to ongoing star-formation activity in the quasar hosts, providing further evidence that the quasar hosts studied here harbor intense starbursts in addition to their active nucleus.
[3]  [pdf] - 1538971
The outflow of gas from the Centaurus A circumnuclear disk: atomic spectral line maps from Herschel-PACS and APEX
Comments: 14 pages, 9 figures. Accepted by Astronomy & Astrophysics
Submitted: 2016-11-17
The physical state of the gas in the central 500 pc of NGC~5128 (the radio galaxy Centaurus A - Cen A), was investigated using the far-infrared fine-structure lines of carbon, oxygen, and nitrogen, as well as the CO(4-3) molecular line. The circumnuclear disk (CND) is traced by emission from dust and the neutral gas ([CI] and CO). A gas outflow with a line-of-sight velocity of 60 km/s is evident in both species. The center of the CND is bright in [OI], [OIII], and [CII]; [OI]63mu emission dominates that of [CII] even though it is absorbed with optical depths of 1.0-1.5. The outflow is well-traced by the [NII] and [NIII] lines and also seen in the [CII] and [OIII] lines that peak in the center. Ionized gas densities are moderate in the CND and low everywhere else. Neutral gas densities range from 4000 per cm3 (outflow, extended thin disk ETD) to 20 000 per cm3 (CND). The CND radiation field is weak compared to the ETD starburst field. The outflow has a much stronger radiation field. The total mass of all the CND gas is 9 x 10^(7) M(o) and the mass of the outflowing gas is only 15%-30% of that. The outflow most likely originates from the shock-dominated CND cavity surrounding the central black hole. With a factor of three uncertainty, the mass outflow rate is about 2 M(o)/yr, a thousand times higher than the accretion rate of the black hole. Without replenishment, the CND will be depleted in 15-120 million years. However, the outflow velocity is well below the escape velocity.
[4]  [pdf] - 1532341
The mid-IR water and silicate relation in protoplanetary disks
Submitted: 2016-10-31
Mid-IR water lines from protoplanetary disks around T Tauri stars have a detection rate of 50\%. Models have identified multiple physical properties of disks such as dust-to-gas mass ratio, dust size power law distribution, disk gas mass, disk inner radius, and disk scale height as potential explanation for the current detection rate. We search for a connection between mid-IR water line fluxes and the strength of the 10~$\mu$m silicate feature. We analyse observed water line fluxes from three blends and compute the 10~$\mu$m silicate feature strength from Spitzer spectra. We use a series of published models, exploring disk dust and gas properties, and the effects of different stars. The models also show that the increasing stellar luminosity enhance simultaneously the strength of this dust feature and the water lines fluxes. No correlation is found between the observed mid-IR water lines and the 10~$\mu$m silicate. Our sample shows the same difference in the peak strength between amorphous and crystalline silicates that was noted in earlier studies, but our models do not support this intrinsic difference in silicate peak strength. Individual properties of our model series are not able to reproduce the most extreme observations, suggesting that more complex dust properties are required. A parametrized settling prescription is able to boost the peak strength by a factor 2 for the standard model. Water line fluxes are unrelated to the composition of the dust. The pronounced regular trends seen in the model results are washed out in the data due to the larger diversity in stellar and disk properties compared to our model series. The disks with with weaker mid-IR water line fluxes are depleted in gas or enhanced in dust in the inner 10~au. In the case of gas depleted disks, settling produces very strong 10~$\mu$m silicate features, with strong peak strength.
[5]  [pdf] - 1531656
HCO+ Detection of Dust-Depleted Gas in the Inner Hole of the LkCa 15 Pre-Transitional Disk
Comments: ApJ accepted, 46 pages (including the appendix), 18 figures
Submitted: 2016-09-19
LkCa 15 is an extensively studied star in the Taurus region known for its pre-transitional disk with a large inner cavity in dust continuum and normal gas accretion rate. The most popular hypothesis to explain the LkCa 15 data invokes one or more planets to carve out the inner cavity, while gas continues to flow across the gap from the outer disk onto the central star. We present spatially unresolved HCO+ J=4-3 observations of the LkCa 15 disk from the JCMT and model the data with the ProDiMo code. We find that: (1) HCO+ line-wings are clearly detected, certifying the presence of gas in the cavity within <50 AU of the star. (2) Reproducing the observed line-wing flux requires both a significant suppression of cavity dust (by a factor >10^4 compared to the ISM) and a substantial increase in the gas scale-height within the cavity (H_0/R_0 ~ 0.6). An ISM dust-to-gas ratio (d:g=10^-2) yields too little line-wing flux regardless of the scale-height or cavity gas geometry, while a smaller scale-height also under predicts the flux even with a reduced d:g. (3) The cavity gas mass is consistent with the surface density profile of the outer disk extended inwards to the sublimation radius (corresponding to mass M_d ~ 0.03 M_sun), and masses lower by a factor >10 appear to be ruled out.
[6]  [pdf] - 1521344
Modelling Mechanical Heating in Star-Forming Galaxies: CO and 13CO Line Ratios as Sensitive Probes
Comments: 17 pages, 12 figures
Submitted: 2016-09-08
We apply photo-dissociation region (PDR) molecular line emission models, that have varying degrees of enhanced mechanical heating rates, to the gaseous component of simulations of star-forming galaxies taken from the literature. Snapshots of these simulations are used to produce line emission maps for the rotational transitions of the CO molecule and its 13CO isotope up to J = 4-3. We consider two galaxy models: a small disk galaxy of solar metallicity and a lighter dwarf galaxy with 0.2 \zsun metallicity. Elevated excitation temperatures for CO(1 - 0) correlate positively with mechanical feedback, that is enhanced towards the central region of both model galaxies. The emission maps of these model galaxies are used to compute line ratios of CO and 13CO transitions. These line ratios are used as diagnostics where we attempt to match them These line ratios are used as diagnostics where we attempt to match them to mechanically heated single component (i.e. uniform density, Far-UV flux, visual extinction and velocity gradient) equilibrium PDR models. We find that PDRs ignoring mechanical feedback in the heating budget over-estimate the gas density by a factor of 100 and the far-UV flux by factors of ~10 - 1000. In contrast, PDRs that take mechanical feedback into account are able to fit all the line ratios for the central < 2 kpc of the fiducial disk galaxy quite well. The mean mechanical heating rate per H atom that we recover from the line ratio fits of this region varies between $10^{-27}$ -- $10^{-26}$~erg s$^{-1}$.
[7]  [pdf] - 1521345
Constraining cloud parameters using high density gas tracers in galaxies
Comments: 13 pages, 11 figures
Submitted: 2016-09-08
Far-infrared molecular emission is an important tool used to understand the excitation mechanisms of the gas in the inter-stellar medium of star-forming galaxies. In the present work, we model the emission from rotational transitions with critical densities n >~ 10^4 cm-3. We include 4-3 < J <= 15-14 transitions of CO and 13CO, in addition to J <= 7-6 transitions of HCN, HNC, and HCO+ on galactic scales. We do this by re-sampling high density gas in a hydrodynamic model of a gas-rich disk galaxy, assuming that the density field of the interstellar medium of the model galaxy follows the probability density function (PDF) inferred from the resolved low density scales. We find that in a narrow gas density PDF, with a mean density of ~10 cm-3 and a dispersion \sigma = 2.1 in the log of the density, most of the emission of molecular lines, emanates from the 10-1000 cm-3 part of the PDF. We construct synthetic emission maps for the central 2 kpc of the galaxy and fit the line ratios of CO and 13CO up to J = 15-14, as well as HCN, HNC, and HCO+ up to J = 7-6, using one photo-dissociation region (PDR) model. We attribute the goodness of the one component fits for our model galaxy to the fact that the distribution of the luminosity, as a function of density, is peaked at gas densities between 10 and 1000 cm-3. We explore the impact of different log-normal density PDFs on the distribution of the line-luminosity as a function of density, and we show that it is necessary to have a broad dispersion, corresponding to Mach numbers >~ 30 in order to obtain significant emission from n > 10^4 cm-3 gas. Such Mach numbers are expected in star-forming galaxies, LIRGS, and ULIRGS. By fitting line ratios of HCN(1-0), HNC(1-0), and HCO+(1-0) for a sample of LIRGS and ULIRGS using mechanically heated PDRs, we constrain the Mach number of these galaxies to 29 < M < 77.
[8]  [pdf] - 1483503
Multi-phase Nature of a Radiation-Driven Fountain with Nuclear Starburst in a Low-mass Active Galactic Nucleus
Comments: 9 pages, 5 figures. Accepted for ApJL. A movie file for Fig.5b can be downloaded from
Submitted: 2016-08-24
The structures and dynamics of molecular, atomic, and ionized gases are studied around a low-luminosity active galactic nucleus (AGN) with a small ($2\times 10^6 M_\odot$) black hole using 3D radiation hydrodynamic simulations. We studied, for the first time, the non-equilibrium chemistry for the X-ray dominated region in the "radiation-driven fountain" (Wada 2012) with supernova feedback. A double hollow cone structure is naturally formed without postulating a thick "torus" around a central source. The cone is occupied with an inhomogeneous, diffuse ionized gas and surrounded by a geometrically thick ($h/r \gtrsim 1$) atomic gas. Dense molecular gases are distributed near the equatorial plane, and energy feedback from supernovae enhances their scale height. Molecular hydrogen exists in a hot phase ( > 1000 K) as well as in a cold ( < 100 K), dense ( >$10^3$ cm$^{-3}$) phase. The velocity dispersion of H$_2$ in the vertical direction is comparable to the rotational velocity, which is consistent with near infrared observations of nearby Seyfert galaxies. Using 3D radiation transfer calculations for the dust emission, we find polar emission in the mid-infrared band (12$\mu m$), which is associated with bipolar outflows, as suggested in recent interferometric observations of nearby AGNs. If the viewing angle for the nucleus is larger than 75 deg, the spectral energy distribution (~ 2 -- 60 $\mu m$) of this model is consistent with that of the Circinus galaxy. The multi-phase interstellar medium observed in optical/infrared and X-ray observations is also discussed.
[9]  [pdf] - 1531028
The SCUBA-2 Cosmology Legacy Survey: 850um maps, catalogues and number counts
Comments: Submitted to MNRAS, comments welcome. Catalogue and maps at
Submitted: 2016-07-13
We present a catalogue of nearly 3,000 submillimetre sources detected at 850um over ~5 square degrees surveyed as part of the James Clerk Maxwell Telescope (JCMT) SCUBA-2 Cosmology Legacy Survey (S2CLS). This is the largest survey of its kind at 850um, probing a meaningful cosmic volume at the peak of star formation activity and increasing the sample size of submillimetre galaxies selected at 850um by an order of magnitude. We describe the wide 850um survey component of S2CLS, which covers the key extragalactic survey fields: UKIDSS-UDS, COSMOS, Akari-NEP, Extended Groth Strip, Lockman Hole North, SSA22 and GOODS-North. The average 1-sigma depth of S2CLS is 1.2 mJy/beam, approaching the SCUBA-2 850um confusion limit, which we determine to be ~0.8 mJy/beam. We measure the single dish 850um number counts to unprecedented accuracy, reducing the Poisson errors on the differential counts to approximately 4% at S_850~3mJy. With several independent fields, we investigate field-to-field variance, finding that the number counts on 0.5-1 degree scales are generally within 50% of the S2CLS mean for S_850>3mJy, with scatter consistent with the Poisson and estimated cosmic variance uncertainties, although there is a marginal (2-sigma) density enhancement in the GOODS-North field. The observed number counts are in reasonable agreement with recent phenomenological and semi-analytic models. Finally, the large solid angle of S2CLS allows us to measure the bright-end counts: at S_850>10mJy there are approximately ten sources per square degree, and we detect the distinctive up-turn in the number counts indicative of the detection of local sources of 850um emission and strongly lensed high-redshift galaxies. Here we describe the data collection and reduction procedures and present calibrated maps and a catalogue of sources; these are made publicly available.
[10]  [pdf] - 1426367
Surface chemistry in photodissociation regions
Comments: 18 pages and 17 figures
Submitted: 2016-04-22
The presence of dust can strongly affect the chemical composition of the interstellar medium. We model the chemistry in photodissociation regions (PDRs) using both gas-phase and dust-phase chemical reactions. Our aim is to determine the chemical compositions of the interstellar medium (gas/dust/ice) in regions with distinct (molecular) gas densities that are exposed to radiation fields with different intensities. We have significantly improved the Meijerink PDR code by including 3050 new gas-phase chemical reactions and also by implementing surface chemistry. In particular, we have included 117 chemical reactions occurring on grain surfaces covering different processes, such as adsorption, thermal desorption, chemical desorption, two-body reactions, photo processes, and cosmic-ray processes on dust grains. We obtain abundances for different gas and solid species as a function of visual extinction, depending on the density and radiation field. We also analyse the rates of the formation of CO2 and H2O ices in different environments. In addition, we study how chemistry is affected by the presence/absence of ice mantles (bare dust or icy dust) and the impact of considering different desorption probabilities. The type of substrate (bare dust or icy dust) and the probability of desorption can significantly alter the chemistry occurring on grain surfaces, leading to differences of several orders of magnitude in the abundances of gas-phase species, such as CO, H2CO, and CH3OH. The type of substrate, together with the density and intensity of the radiation field, also determine the threshold extinction to form ices of CO2 and H2O. We also conclude that H2CO and CH3OH are mainly released into the gas phase of low, far-ultraviolet illuminated PDRs through chemical desorption upon two-body surface reactions, rather than through photodesorption.
[11]  [pdf] - 1351565
Consistent dust and gas models for protoplanetary disks. I. Disk shape, dust settling, opacities, and PAHs
Comments: accepted by A&A, 36 pages, 30 figures, 5 tables
Submitted: 2015-11-11
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. We propose new standard dust opacities for disk models, we present a simplified treatment of PAHs sufficient to reproduce the PAH emission features, and we suggest using a simple treatment of dust settling. We roughly adjust parameters to obtain a model that predicts typical Class II T Tauri star continuum and line observations. 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 continuum and line observables, in particular on the SED, mm-slope, continuum visibilities, and emission lines including [OI] 63um, high-J CO lines, (sub-)mm CO isotopologue lines, and CO fundamental ro-vibrational lines. We find that evolved dust properties (large grains) often needed to fit the SED, have important consequences for disk chemistry and heating/cooling balance, leading to stronger 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 shield the gas from stellar UV radiation because of their strong absorption and negligible scattering opacities. 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 some key properties of the disks, such as disk shape and mass, opacities, and the dust/gas ratio, by simultaneously fitting continuum and line observations.
[12]  [pdf] - 1331062
Mid-IR spectra of Pre-Main Sequence Herbig stars: an explanation for the non-detections of water lines
Submitted: 2015-10-27
The mid-IR detection rate of water lines in disks around Herbig stars disks is about 5\%, while it is around 50\% for disks around TTauri stars. The reason for this is still unclear. In this study, we want to find an explanation for the different detection rates between low mass and high mass pre-main-sequence stars (PMSs) in the mid-IR regime. We run disk models with stellar parameters adjusted to spectral types B9 through M2, using the radiation thermo-chemical disk modeling code ProDiMo. We produce convolved spectra at the resolution of Spitzer IRS, JWST MIRI and VLT VISIR spectrographs. We apply random noise derived from typical Spitzer spectra for a direct comparison with observations. The strength of the mid-IR water lines correlates directly with the luminosity of the central star. We explored a small parameter space around a standard disk model, considering dust-to-gas mass ratio, disk gas mass, mixing coefficient for dust settling, flaring index, dust maximum size and size power law distribution index. The models show that it is possible to suppress the water emission, however, current observations are not sensitive enough to detect mid-IR lines in disks for most of the explored parameters. The presence of noise in the spectra, combined with the high continuum flux (noise level is proportional to the continuum flux), is the most likely explanation for the non detections towards Herbig stars. Mid-IR spectra with resolution higher than 20000 are needed to investigate water in protoplanetary disks. Intrinsic differences in disk structure, e.g. inner gaps, gas-to-dust ratio, dust size and distribution, and inner disk scale height, between Herbig and TTauri star disks are able to explain a lower water detection rate in disks around Herbig stars.
[13]  [pdf] - 1296355
Understanding the water emission in the mid- and far-IR from protoplanetary disks around T~Tauri stars
Submitted: 2015-10-06
We investigate which properties of protoplanetary disks around T Tauri stars affect the physics and chemistry in the regions where mid- and far-IR water lines originate and their respective line fluxes. We search for diagnostics for future observations. With the code ProDiMo, we build a series of models exploring a large parameter space, computing rotational and rovibrational transitions of water in nonlocal thermodynamic equilibrium (non-LTE). We select a sample of transitions in the mid- IR regime and the fundamental ortho and para water transitions in the far-IR. We investigate the chemistry and the local physical conditions in the line emitting regions. We calculate Spitzer spectra for each model and compare far-IR and mid-IR lines. In addition, we use mid-IR colors to tie the water line predictions to the dust continuum. Parameters affecting the water line fluxes in disks by more than a factor of three are : the disk gas mass, the dust-to-gas mass ratio, the dust maximum grain size, ISM(InterStellarMedium) UV radiation field, the mixing parameter of Dubrulle settling, the disk flaring parameter, and the dust size distribution. The first four parameters affect the mid-IR lines much more than the far-IR lines. A key driver behind water spectroscopy is the dust opacity, which sets the location of the water line emitting region. We identify three types of parameters. Parameters, such as dust-to-gas ratio, ISM radiation field, and dust size distribution, affect the mid-IR lines more, while the far-IR transitions are more affected by the flaring index. The gas mass greatly affects lines in both regimes. Higher spectral resolution and line sensitivities, like from the James Webb Space Telescope, are needed to detect a statistically relevant sample of individual water lines to distinguish further between these types of parameters.
[14]  [pdf] - 1245790
The SCUBA-2 Cosmology Legacy Survey: ALMA resolves the bright-end of the sub-millimeter number counts
Comments: 14 pages, 6 figures. ApJ in press
Submitted: 2015-05-19
We present high-resolution 870-um ALMA continuum maps of 30 bright sub-millimeter sources in the UKIDSS UDS field. These sources are selected from deep, 1-square degrees 850-um maps from the SCUBA--2 Cosmology Legacy Survey, and are representative of the brightest sources in the field (median SCUBA2 flux S_850=8.7+/-0.4 mJy). We detect 52 sub-millimeter galaxies (SMGs) at >4-sigma significance in our 30 ALMA maps. In 61+/-17% 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 L_IR>10^12 Lo. The brightest SMG contributes on average 80+/-4% of the single-dish flux density, and in the ALMA maps containing >=2 SMGs the secondary SMG contributes 25+/-3% of the integrated ALMA flux. We construct source counts and show that multiplicity boosts the apparent single-dish cumulative counts by 20% at S_870>7.5mJy, and by 60% at S_870>12mJy. 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.6mJy. The break corresponds to a luminosity of ~6x10^12Lsol or a star-formation rate of ~1000Mo/yr. For the typical sizes of these SMGs, which are resolved in our ALMA data with r=1.2+/-0.1kpc, this yields a limiting SFR density of ~100Msol/yr/kpc2. Finally, the number density of S_870>2mJy 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 must be physically associated.
[15]  [pdf] - 1224166
Probing stellar accretion with mid-infrared hydrogen lines
Submitted: 2015-01-25
In this paper we investigate the origin of the mid-infrared (IR) hydrogen recombination lines for a sample of 114 disks in different evolutionary stages (full, transitional and debris disks) collected from the {\it Spitzer} archive. We focus on the two brighter {H~{\sc i}} lines observed in the {\it Spitzer} spectra, the {H~{\sc i}}(7-6) at 12.37$\mu$m and the {H~{\sc i}}(9-7) at 11.32$\mu$m. We detect the {H~{\sc i}}(7-6) line in 46 objects, and the {H~{\sc i}}(9-7) in 11. We compare these lines with the other most common gas line detected in {\it Spitzer} spectra, the {[Ne~{\sc iii}]} at 12.81$\mu$m. We argue that it is unlikely that the {H~{\sc i}} emission originates from the photoevaporating upper surface layers of the disk, as has been found for the {[Ne~{\sc iii}]} lines toward low-accreting stars. Using the {H~{\sc i}}(9-7)/{H~{\sc i}}(7-6) line ratios we find these gas lines are likely probing gas with hydrogen column densities of 10$^{10}$-10$^{11}$~cm$^{-3}$. The subsample of objects surrounded by full and transitional disks show a positive correlation between the accretion luminosity and the {H~{\sc i}} line luminosity. These two results suggest that the observed mid-IR {H~{\sc i}} lines trace gas accreting onto the star in the same way as other hydrogen recombination lines at shorter wavelengths. A pure chromospheric origin of these lines can be excluded for the vast majority of full and transitional disks.We report for the first time the detection of the {H~{\sc i}}(7-6) line in eight young (< 20~Myr) debris disks. A pure chromospheric origin cannot be ruled out in these objects. If the {H~{\sc i}}(7-6) line traces accretion in these older systems, as in the case of full and transitional disks, the strength of the emission implies accretion rates lower than 10$^{-10}$M$_{\odot}$/yr. We discuss some advantages of extending accretion indicators to longer wavelengths.
[16]  [pdf] - 1223995
The Herschel Comprehensive (U)LIRG Emission Survey (HerCULES): CO Ladders, fine structure lines, and neutral gas cooling
Submitted: 2015-01-13
(Ultra) Luminous Infrared Galaxies ((U)LIRGs) are objects characterized by their extreme infrared (8-1000 $\mu$m) luminosities ($L_{LIRG}>10^{11} $L$_\odot$ and $L_{ULIRG}>10^{12}$ L$_\odot$). The Herschel Comprehensive ULIRG Emission Survey (HerCULES; PI van der Werf) presents a representative flux-limited sample of 29 (U)LIRGs that spans the full luminosity range of these objects (10$^{11}\leq L_\odot \geq10^{13}$). With the \emph{Herschel Space Observatory}, we observe [CII] 157 $\mu$m, [OI] 63 $\mu$m, and [OI] 145 $\mu$m line emission with PACS, CO J=4-3 through J=13-12, [CI] 370 $\mu$m, and [CI] 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 [OI] 63 $\mu$m emission line is self absorbed. Comparing the CO excitation to the IRAS 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 [CII], [SiII], [OI], and [CI] lines in the objects with the highest far IR fluxes, but do not observe this for CO $4\leq J_{upp}\leq13$. 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 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.
[17]  [pdf] - 1451827
High-lying OH absorption, [C II] deficits, and extreme $L_{\mathrm{FIR}}/M_{\mathrm{H2}}$ ratios in galaxies
Comments: 11 pages, 7 figures, accepted for publication in The Astrophysical Journal
Submitted: 2014-12-15
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.
[18]  [pdf] - 1223235
The SCUBA-2 Cosmology Legacy Survey: ALMA resolves the rest-frame far-infrared emission of sub-millimeter galaxies
Comments: 16 pages, 7 figures, accepted for publication in ApJ
Submitted: 2014-11-18
We present high-resolution (0.3'') ALMA 870um imaging of 52 sub-millimeter galaxies (SMGs) in the Ultra Deep Survey (UDS) field and investigate the size and morphology of the sub-millimeter (sub-mm) emission on 2-10kpc scales. We derive a median intrinsic angular size of FWHM=0.30$\pm$0.04'' for the 23 SMGs in the sample detected at a signal-to-noise ratio (SNR) >10. Using the photometric redshifts of the SMGs we show that this corresponds to a median physical half-light diameter of 2.4$\pm$0.2kpc. A stacking analysis of the SMGs detected at an SNR <10 shows they have sizes consistent with the 870um-bright SMGs in the sample. We compare our results to the sizes of SMGs derived from other multi-wavelength studies, and show that the rest-frame ~250um sizes of SMGs are consistent with studies of resolved 12CO (J=3-2 to 7-6) emission lines, but that sizes derived from 1.4GHz imaging appear to be approximately two times larger on average, which we attribute to cosmic ray diffusion. The rest-frame optical sizes of SMGs are around four times larger than the sub-millimeter sizes, indicating that the star formation in these galaxies is compact relative to the pre-existing stellar distribution. The size of the starburst region in SMGs is consistent with the majority of the star formation occurring in a central region, a few kpc in extent, with a median star formation rate surface density of 90$\pm$30Msol/yr/kpc$^2$, which may suggest that we are witnessing an intense period of bulge growth in these galaxies.
[19]  [pdf] - 1451813
ALMA Observations of Warm Dense Gas in NGC 1614 --- Breaking of Star Formation Law in the Central kpc
Comments: 13 pages, 10 figures; accepted by ApJ
Submitted: 2014-11-04
We present ALMA Cycle-0 observations of the CO (6-5) line emission and of the 435um dust continuum emission in the central kpc of NGC 1614, a local luminous infrared galaxy (LIRG) at a distance of 67.8 Mpc (1 arcsec = 329 pc). The CO emission is well resolved by the ALMA beam (0".26 x 0".20) into a circum-nuclear ring, with an integrated flux of f_{CO(6-5)} = 898 (+-153) Jy km/s, which is 63(+-12)% of the total CO(6-5) flux measured by Herschel. The molecular ring, located between 100pc < r < 350pc from the nucleus, looks clumpy and includes seven unresolved (or marginally resolved) knots with median velocity dispersion of 40 km/s. These knots are associated with strong star formation regions with \Sigma_{SFR} 100 M_\sun/yr/kpc^{2} and \Sigma_{Gas} 1.0E4 M_\sun/pc^{2}. The non-detections of the nucleus in both the CO (6-5) line emission and the 435um continuum rule out, with relatively high confidence, a Compton-thick AGN in NGC 1614. Comparisons with radio continuum emission show a strong deviation from an expected local correlation between \Sigma_{Gas} and \Sigma_{SFR}, indicating a breakdown of the Kennicutt-Schmidt law on the linear scale of 100 pc.
[20]  [pdf] - 903674
Songlines from Direct Collapse Seed Black Holes: Effects of X-rays on Black Hole Growth and Stellar Populations
Comments: 8 pages, 6 figures. Accepted for publication in ApJ, with minor changes to submitted version
Submitted: 2014-09-01, last modified: 2014-10-13
In the last decade, the growth of supermassive black holes (SMBHs) has been intricately linked to galaxy formation and evolution and is a key ingredient in the assembly of galaxies. To investigate the origin of SMBHs, we perform cosmological simulations that target the direct collapse black hole (DCBH) seed formation scenario in the presence of two different strong Lyman-Werner (LW) background fields. These simulations include the X-ray irradiation from a central massive black hole (MBH), $\rm{H}_2$ self-shielding and stellar feedback from metal-free and metal-enriched stars. We find in both simulations that local X-ray feedback induces metal-free star formation $\sim 0.5$ Myr after the MBH forms. The MBH accretion rate reaches a maximum of $10^{-3}$ $M_{\odot}$ yr$^{-1}$ in both simulations. However, the duty cycle differs which is derived to be $6\%$ and $50\%$ for high and low LW cases, respectively. The MBH in the high LW case grows only $\sim 6\%$ in 100 Myr compared to $16\%$ in the low LW case. We find that the maximum accretion rate is determined by the local gas thermodynamics whereas the duty cycle is determined by the large scale gas dynamics and gas reservoir. We conclude that radiative feedback from the central MBH plays an important role in star formation in the nuclear regions and stifling initial MBH growth, relative to the typical Eddington rate argument, and that initial MBH growth might be affected by the local LW radiation field.
[21]  [pdf] - 1215713
Star formation relations and CO SLEDs across the J-ladder and redshift
Comments: (18 pages, 6 figures. Accepted for publication in ApJ)
Submitted: 2014-07-16
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.
[22]  [pdf] - 863072
Molecular gas heating in Arp 299
Submitted: 2014-07-08
Understanding the heating and cooling mechanisms in nearby (Ultra) luminous infrared galaxies can give us insight into the driving mechanisms in their more distant counterparts. Molecular emission lines play a crucial role in cooling excited gas, and recently, with Herschel Space Observatory we have been able to observe the rich molecular spectrum. CO is the most abundant and one of the brightest molecules in the Herschel wavelength range. CO transitions are observed with Herschel, and together, these lines trace the excitation of CO. We study Arp 299, a colliding galaxy group, with one component harboring an AGN and two more undergoing intense star formation. For Arp 299 A, we present PACS spectrometer observations of high-J CO lines up to J=20-19 and JCMT observations of $^{13}$CO and HCN to discern between UV heating and alternative heating mechanisms. There is an immediately noticeable difference in the spectra of Arp 299 A and Arp 299 B+C, with source A having brighter high-J CO transitions. This is reflected in their respective spectral energy line distributions. We find that photon-dominated regions (PDRs) are unlikely to heat all the gas since a very extreme PDR is necessary to fit the high-J CO lines. In addition, this extreme PDR does not fit the HCN observations, and the dust spectral energy distribution shows that there is not enough hot dust to match the amount expected from such an extreme PDR. Therefore, we determine that the high-J CO and HCN transitions are heated by an additional mechanism, namely cosmic ray heating, mechanical heating, or X-ray heating. We find that mechanical heating, in combination with UV heating, is the only mechanism that fits all molecular transitions. We also constrain the molecular gas mass of Arp 299 A to 3e9 Msun and find that we need 4% of the total heating to be mechanical heating, with the rest UV heating.
[23]  [pdf] - 1209572
Extended [CII] Emission in Local Luminous Infrared Galaxies
Comments: Accepted for publication in ApJ Letters. 5 pages, 2 figures
Submitted: 2014-05-15
We present Herschel/PACS observations of extended [CII]157.7{\mu}m line emission detected on ~ 1 - 10 kpc scales in 60 local luminous infrared galaxies (LIRGs) from the Great Observatories All-sky LIRG Survey (GOALS). We find that most of the extra-nuclear emission show [CII]/FIR ratios >~ 4 x 10^-3, larger than the mean ratio seen in the nuclei, and similar to those found in the extended disks of normal star-forming galaxies and the diffuse inter-stellar medium (ISM) of our Galaxy. The [CII] "deficits" found in the most luminous local LIRGs are therefore restricted to their nuclei. There is a trend for LIRGs with warmer nuclei to show larger differences between their nuclear and extra-nuclear [CII]/FIR ratios. We find an anti-correlation between [CII]/FIR and the luminosity surface density, {\Sigma}_IR, for the extended emission in the spatially-resolved galaxies. However, there is an offset between this trend and that found for the LIRG nuclei. We use this offset to derive a beam filling-factor for the star-forming regions within the LIRG disks of ~ 6 % relative to their nuclei. We confront the observed trend to photo-dissociation region (PDR) models and find that the slope of the correlation is much shallower than the model predictions. Finally, we compare the correlation found between [CII]/FIR and {\Sigma}_IR with measurements of high-redshift starbursting IR-luminous galaxies.
[24]  [pdf] - 1208820
ALMA Observations of Warm Molecular Gas and Cold Dust in NGC 34
Comments: 10 pages, 13 figures, accepted for publication in Astrophysical Journal
Submitted: 2014-04-04
We present ALMA Cycle-0 observations of the CO (6-5) line emission (rest-frame frequency = 691.473 GHz) and of the 435$\mu m$ dust continuum emission in the nuclear region of NGC 34, a local luminous infrared galaxy (LIRG) at a distance of 84 Mpc (1" = 407 pc) which contains a Seyfert 2 active galactic nucleus (AGN) and a nuclear starburst. The CO emission is well resolved by the ALMA beam ($\rm 0."26\times 0."23$), with an integrated flux of $\rm f_{CO~(6-5)} = 1004\; (\pm 151) \; Jy\; km\; s^{-1}$. Both the morphology and kinematics of the CO (6-5) emission are rather regular, consistent with a compact rotating disk with a size of 200 pc. A significant emission feature is detected on the red-shifted wing of the line profile at the frequency of the $\rm H^{13}CN\; (8-7)$ line, with an integrated flux of $\rm 17.7 \pm 2.1 (random) \pm 2.7 (sysmatic)\; Jy\;km\; s^{-1}$. However, it cannot be ruled out that the feature is due to an outflow of warm dense gas with a mean velocity of $\rm 400\; km\; s^{-1}$. The continuum is resolved into an elongated configuration, and the observed flux corresponds to a dust mass of $\rm M_{dust} = 10^{6.97\pm 0.13}\; M_{sun}$. An unresolved central core ($\rm radius \simeq 50\; pc$) contributes $28\%$ of the continuum flux and $19\%$ of the CO (6-5) flux, consistent with insignificant contributions of the AGN to both emissions. Both the CO (6-5) and continuum spatial distributions suggest a very high gas column density ($\rm >= 10^4\; M_{sun}\; pc^{-2}$) in the nuclear region at $\rm radius <= 100\; pc$.
[25]  [pdf] - 1208637
Diagnostics of the molecular component of PDRs with mechanical heating. II: line intensities and ratios
Comments: 38 pages, to appear in A&A
Submitted: 2014-03-27
CO observations in active galactic nuclei and star-bursts reveal high kinetic temperatures. Those environments are thought to be very turbulent due to dynamic phenomena such as outflows and high supernova rates. We investigate the effect of mechanical heating (MH) on atomic fine-structure and molecular lines, and their ratios. We use those ratios as a diagnostic to constrain the amount of MH in an object and also study its significance on estimating the H2 mass. Equilibrium PDRs models were used to compute the thermal and chemical balance for the clouds. The equilibria were solved for numerically using the optimized version of the Leiden PDR-XDR code. Large velocity gradient calculations were done as post-processing on the output of the PDR models using RADEX. High-J CO line ratios are very sensitive to MH. Emission becomes at least one order of magnitude brighter in clouds with n~10^5~cm^-3 and a star formation rate of 1 Solar Mass per year (corresponding to a MH rate of 2 * 10^-19 erg cm^-3 s^-1). Emission of low-J CO lines is not as sensitive to MH, but they do become brighter in response to MH. Generally, for all of the lines we considered, MH increases excitation temperatures and decreases the optical depth at the line centre. Hence line ratios are also affected, strongly in some cases. Ratios involving HCN are a good diagnostic for MH, such as HCN(1-0)/CO(1-0) and HCN(1-0)/HCO^+(1-0). Both ratios increase by a factor 3 or more for a MH equivalent to > 5 percent of the surface heating, as opposed to pure PDRs. The first major conclusion is that low-J to high-J intensity ratios will yield a good estimate of the MH rate (as opposed to only low-J ratios). The second one is that the MH rate should be taken into account when determining A_V or equivalently N_H, and consequently the cloud mass. Ignoring MH will also lead to large errors in density and radiation field estimates.
[26]  [pdf] - 841573
Disk evolution, element abundances and cloud properties of young gas giant planets
Comments: 32 pages, accepted for publication in the OA journal 'Life', Special issue 'Planet Formation and the Rise of Life'. Life (ISSN 2075-1729)
Submitted: 2014-03-18
We discuss the chemical pre-conditions for planet formation, in terms of gas and ice abundances in a protoplanetary disk, as function of time and position, and the resulting chemical composition and cloud properties in the atmosphere when young gas giant planets form. Large deviations between the abundances of the host star and its gas giants seem likely to occur if the planet formation follows the core-accretion scenario. These deviations stem from the separate evolution of gas and dust in the disk, where the dust forms the planet cores, followed by the final run-away accretion of the left-over gas. ProDiMo protoplanetary disk models are used to predict the chemical evolution of gas and ice in the midplane. We find that cosmic rays play a crucial role in slowly un-blocking the CO, where the liberated oxygen forms water, which then freezes out quickly. Therefore, the C/O ratio in the gas phase is found to gradually increase with time, in a region bracketed by the water and CO ice-lines. In this regions, C/O is found to approach unity after about 5 Myrs, scaling with the cosmic ray ionisation rate. We then explore how the atmospheric chemistry and cloud properties in young gas giants are affected when the non-solar C/O ratios predicted by the disk models are assumed. The DRIFT cloud formation model is applied to study the formation of atmospheric clouds under the influence of varying primordial element abundances and its feedback onto the local gas. We demonstrate that element depletion by cloud formation plays a crucial role in converting an oxygen-rich atmosphere gas into carbon-rich gas when non-solar, primordial element abundances are considered as suggested by disk model.
[27]  [pdf] - 833348
[OI] disk emission in the Taurus star forming region
Comments: 9 pages, accepted for publication in A&A
Submitted: 2014-02-11
The structure of protoplanetary disks is thought to be linked to the temperature and chemistry of their dust and gas. Whether the disk is flat or flaring depends on the amount of radiation that it absorbs at a given radius, and on the efficiency with which this is converted into thermal energy. The understanding of these heating and cooling processes is crucial to provide a reliable disk structure for the interpretation of dust continuum emission and gas line fluxes. Especially in the upper layers of the disk, where gas and dust are thermally decoupled, the infrared line emission is strictly related to the gas heating/cooling processes. We aim to study the thermal properties of the disk in the oxygen line emission region, and to investigate the relative importance of X-ray (1-120 Angstrom) and far-UV radiation (FUV, 912-2070 Angstrom) for the heating balance there. We use [OI] 63 micron line fluxes observed in a sample of protoplanetary disks of the Taurus/Auriga star forming region and compare it to the model predictions presented in our previous work. The data were obtained with the PACS instrument on board the Herschel Space Observatory as part of the Herschel Open Time Key Program GASPS (GAS in Protoplanetary diskS), published in Howard et al. (2013). Our theoretical grid of disk models can reproduce the [OI] absolute fluxes and predict a correlation between [OI] and the sum Lx+Lfuv. The data show no correlation between the [OI] line flux and the X-ray luminosity, the FUV luminosity or their sum. The data show that the FUV or X-ray radiation has no notable impact on the region where the [OI] line is formed. This is in contrast with what is predicted from our models. Possible explanations are that the disks in Taurus are less flaring than the hydrostatic models predict, and/or that other disk structure aspects that were left unchanged in our models are important. ..abridged..
[28]  [pdf] - 1203188
The molecular circumnuclear disk (CND) in Centaurus A: A multi-transition CO and [CI] survey with Herschel, APEX, JCMT, and SEST
Comments: 14 pages, 5 figures; A&A in press
Submitted: 2014-02-05
We present new CO and C^o line measurements of the compact circumnuclear disk in the center of NGC 128 (Centaurus~A) obtained with the Herschel Space Observatory, as well as SEST, JCMT, and APEX. The Cen A center CO ladder is quite different from those of either star-burst galaxies or AGNs. In addition, the relative intensity of the central Cen A [CI] emission lines is much greate than that in any other galaxy. The CO surface brightness of the compact circumnuclear disk (CND) is significantly higher than that of the much more extended thin disk (ETD) in the same line of sight. Our LVG and PDR/XDR models suggest that much of the CND gas is relatively cool (25 - 80 K) and not very dense (~ 300 cm^{-3}) if the heating is by UV photons, although there is some gas in both the CND and the ETD with a much higher density of ~30 000 cm^{-3}. Finally, there is also high-excitation, high-density phase in the CND (but not in the ETD), either in the form of an extreme PDR but more likely in the form of an XDR. The total gas mass of the Cen A CND is 8.4 x 10^{7} M(sun), uncertain by a factor of two. The CO-H2 conversion factor is 4 x 10^{20} K km/s, also within a factor of two.
[29]  [pdf] - 791951
Probing the gaseous disk of T Tau N with CN 5-4 lines
Comments: 6 pages, 3 figures, accepted for publication on ApJ Letters
Submitted: 2014-01-30
We present spectrally resolved Herschel/HIFI observations of the young multiple system T Tau in atomic and molecular lines. While CO, H2O, [C II], and SO lines trace the envelope and the outflowing gas up to velocities of 33 km/s with respect to systemic, the CN 5-4 hyperfine structure lines at 566.7, 566.9 GHz show a narrow double-peaked profile centered at systemic velocity, consistent with an origin in the outer region of the compact disk of T Tau N. Disk modeling of the T Tau N disk with the thermo-chemical code ProDiMo produces CN line fluxes and profiles consistent with the observed ones and constrain the size of the gaseous disk (R_out = 110 (+10, -20) AU) and its inclination (i = 25 \pm 5 degree). The model indicates that the CN lines originate in a disk upper layer at 40-110 AU from the star, which is irradiated by the stellar UV field and heated up to temperatures of 50-700 K. With respect to previously observed CN 2-1 millimeter lines, the CN 5-4 lines appear to be less affected by envelope emission, due to their larger critical density and excitation temperature. Hence, high-J CN lines are a unique confusion-free tracer of embedded disks, such as the disk of T Tau N.
[30]  [pdf] - 1202860
Radiative and mechanical feedback into the molecular gas of NGC 253
Submitted: 2014-01-20
Starburst galaxies are undergoing intense periods of star formation. Understanding the heating and cooling mechanisms in these galaxies can give us insight to the driving mechanisms that fuel the starburst. Molecular emission lines play a crucial role in the cooling of the excited gas. With SPIRE on the Herschel Space Observatory we have observed the rich molecular spectrum towards the central region of NGC 253. CO transitions from J=4-3 to 13-12 are observed and together with low-J line fluxes from ground based observations, these lines trace the excitation of CO. By studying the CO excitation ladder and comparing the intensities to models, we investigate whether the gas is excited by UV radiation, X-rays, cosmic rays, or turbulent heating. Comparing the $^{12}$CO and $^{13}$CO observations to large velocity gradient models and PDR models we find three main ISM phases. We estimate the density, temperature,and masses of these ISM phases. By adding $^{13}$CO, HCN, and HNC line intensities, we are able to constrain these degeneracies and determine the heating sources. The first ISM phase responsible for the low-J CO lines is excited by PDRs, but the second and third phases, responsible for the mid to high-J CO transitions, require an additional heating source. We find three possible combinations of models that can reproduce our observed molecular emission. Although we cannot determine which of these are preferable, we can conclude that mechanical heating is necessary to reproduce the observed molecular emission and cosmic ray heating is a negligible heating source. We then estimate the mass of each ISM phase; $6\times 10^7$ M$_\odot$ for phase 1 (low-J CO lines), $3\times 10^7$ M$_\odot$ for phase 2 (mid-J CO lines), and $9\times 10^6$ M$_\odot$ for phase 3 (high-J CO lines) for a total system mass of $1\times10^{8}$ M$_\odot$.
[31]  [pdf] - 1179370
Gas lines from the 5-Myr old optically thin disk around HD141569A. Herschel observations and modeling
Comments: accepted for publication in A&A
Submitted: 2013-09-19
At the distance of 99-116 pc, HD141569A is one of the nearest HerbigAe stars that is surrounded by a tenuous disk, probably in transition between a massive primordial disk and a debris disk. We observed the fine-structure lines of OI at 63 and 145 micron and the CII line at 157 micron with the PACS instrument onboard the Herschel Space Telescope as part of the open-time large programme GASPS. We complemented the atomic line observations with archival Spitzer spectroscopic and photometric continuum data, a ground-based VLT-VISIR image at 8.6 micron, and 12CO fundamental ro-vibrational and pure rotational J=3-2 observations. We simultaneously modeled the continuum emission and the line fluxes with the Monte Carlo radiative transfer code MCFOST and the thermo-chemical code ProDiMo to derive the disk gas- and dust properties assuming no dust settling. The models suggest that the oxygen lines are emitted from the inner disk around HD141569A, whereas the [CII] line emission is more extended. The CO submillimeter flux is emitted mostly by the outer disk. Simultaneous modeling of the photometric and line data using a realistic disk structure suggests a dust mass derived from grains with a radius smaller than 1 mm of 2.1E-7 MSun and from grains with a radius of up to 1 cm of 4.9E-6 MSun. We constrained the polycyclic aromatic hydrocarbons (PAH) mass to be between 2E-11 and 1..4E-10 MSun assuming circumcircumcoronene (C150H30) as the representative PAH. The associated PAH abundance relative to hydrogen is lower than those found in the interstellar medium (3E-7) by two to three orders of magnitude. The disk around HD141569A is less massive in gas (2.5 to 4.9E-4 MSun or 67 to 164 MEarth) and has a flat opening angle (<10%). [abridged]
[32]  [pdf] - 1173528
The SCUBA-2 Cosmology Legacy Survey: demographics of the 450 \mu m-population
Comments: 20 pages, 12 figures, accepted for publication in MNRAS
Submitted: 2013-08-20
We investigate the multi-wavelength properties of a sample of 450-\mu m selected sources from the SCUBA-2 Cosmology Legacy Survey (S2CLS). A total of 69 sources were identified above 4\sigma\ in deep SCUBA-2 450-\mu m observations overlapping the UDS and COSMOS fields and covering 210 sq. arcmin to a typical depth of \sigma 450=1.5 mJy. Reliable cross identification are found for 58 sources (84 per cent) in Spitzer and Hubble Space Telescope WFC3/IR data. The photometric redshift distribution (dN/dz) of 450\mu m-selected sources is presented, showing a broad peak in the redshift range 1<z<3, and a median of z=1.4. Combining the SCUBA-2 photometry with Herschel SPIRE data from HerMES, the submm spectral energy distribution (SED) is examined via the use of modified blackbody fits, yielding aggregate values for the IR luminosity, dust temperature and emissivity of <LIR>=10^12 +/- 0.8 L_sol, <T_D>=42 +/- 11 K and <\beta_D>=1.6 +/- 0.5, respectively. The relationship between these SED parameters and the physical properties of galaxies is investigated, revealing correlations between T_D and LIR and between \beta_D and both stellar mass and effective radius. The connection between star formation rate and stellar mass is explored, with 24 per cent of 450 \mu m sources found to be ``star-bursts'', i.e. displaying anomalously high specific SFRs. However, both the number density and observed properties of these ``star-burst'' galaxies are found consistent with the population of normal star-forming galaxies.
[33]  [pdf] - 1173277
Uncertainties in water chemistry in disks: An application to TW Hya
Comments: accepted for publication in A&A
Submitted: 2013-08-08
Context. This paper discusses the sensitivity of water lines to chemical processes and radiative transfer for the protoplanetary disk around TW Hya. The study focuses on the Herschel spectral range in the context of new line detections with the PACS instrument from the Gas in Protoplanetary Systems project (GASPS). Aims. The paper presents an overview of the chemistry in the main water reservoirs in the disk around TW Hya. It discusses the limitations in the interpretation of observed water line fluxes. Methods. ... (abbreviated) Results. We report new line detections of p-H2O (3_22-2_11) at 89.99 micron and CO J=18-17 at 144.78 micron for the disk around TW Hya. Disk modeling shows that the far-IR fine structure lines ([OI], [CII]) and molecular submm lines are very robust to uncertainties in the chemistry, while the water line fluxes can change by factors of a few. The water lines are optically thick, sub-thermally excited and can couple to the background continuum radiation field. The low-excitation water lines are also sensitive to uncertainties in the collision rates, e.g. with neutral hydrogen. The gas temperature plays an important role for the [OI] fine structure line fluxes, the water line fluxes originating from the inner disk as well as the high excitation CO, CH+ and OH lines. Conclusions. Due to their sensitivity on chemical input data and radiative transfer, water lines have to be used cautiously for understanding details of the disk structure. Water lines covering a wide range of excitation energies provide access to the various gas phase water reservoirs (inside and outside the snow line) in protoplanetary disks and thus provide important information on where gas-phase water is potentially located. Experimental and/or theoretical collision rates for H2O with atomic hydrogen are needed to diminish uncertainties from water line radiative transfer.
[34]  [pdf] - 1172966
Using [CI] to probe the Interstellar Medium in z~2.5 Sub-Millimeter Galaxies
Comments: 21 pages, 12 figures, accepted for publication in MNRAS
Submitted: 2013-07-24
We present new [CI](1-0) and 12CO(4-3) Plateau de Bure Interferometer observations of five Sub-Millimeter Galaxies (SMGs) and combine these with all available [CI](1-0) literature detections in SMGs to probe the gas distribution within a sample of 14 systems. We explore the [CI](1-0) properties of the SMG population, particularly investigating the ratio of the [CI](1-0) luminosity to various 12CO transition and far-infrared luminosities. We find that the SMGs with new observations extend the spread of L[CI](1-0)/LFIR to much higher values than found before, with our complete sample providing a good representation of the diverse z>2 SMG population. We compare the line ratios to the outputs of photodissociation region (PDR) models to constrain the average density (n) and radiation field (in terms of the local field value, G0) of the SMGs. We find the SMGs are most comparable to local ULIRGs in G0 and n, however a significant tail of 5 of the 14 SMGs are likely best compared to less compact, local starburst galaxies, providing new evidence that many SMGs have extended star formation distributions and are therefore not simply scaled up versions of local ULIRGs. We explore the limitations of using simple PDR models to understand [CI], which may be concomitant with the bulk H2 mass rather than PDR-distributed. We therefore also assess [CI] as a tracer of H2, finding that for our sample SMGs, the H2 masses derived from [CI] are often consistent with those determined from low excitation 12CO. We conclude that [CI] observations provide a useful tool to probe the bulk gas and gas processes occurring within merging SMGs, however more detailed, resolved observations are required to fully exploit [CI] as a diagnostic.
[35]  [pdf] - 1172611
Explaining the [CII]158um Deficit in Luminous Infrared Galaxies - First Results from a Herschel/PACS Study of the GOALS Sample
Comments: 15 pages, 8 figures. Accepted for publication in ApJ. Table 1 will be available in the electronic version of the paper and also soon at
Submitted: 2013-07-09
We present the first results of a survey of the [CII]158um emission line in 241 luminous infrared galaxies (LIRGs) comprising the Great Observatories All-sky Survey (GOALS) sample, obtained with the PACS instrument on board Herschel. The [CII] luminosities of the LIRGs in GOALS range from ~10^7 to 2x10^9 Lsun. We find that LIRGs show a tight correlation of [CII]/FIR with far-IR flux density ratios, with a strong negative trend spanning from ~10^-2 to 10^-4, as the average temperature of dust increases. We find correlations between the [CII]/FIR ratio and the strength of the 9.7um silicate absorption feature as well as with the luminosity surface density of the mid-IR emitting region (Sigma_MIR), suggesting that warmer, more compact starbursts have substantially smaller [CII]/FIR ratios. Pure star-forming (SF) LIRGs have a mean [CII]/FIR ~ 4x10^-3, while galaxies with low 6.2um PAH equivalent widths (EWs), indicative of the presence of active galactic nuclei (AGN), span the full range in [CII]/FIR. However, we show that even when only pure SF galaxies are considered, the [CII]/FIR ratio drops by an order of magnitude, from 10^-2 to 10^-3, with Sigma_MIR and Sigma_IR, implying that the [CII] luminosity is not a good indicator of the star formation rate (SFR) for most LIRGs, for it does not scale linearly with the warm dust emission. Moreover, even in LIRGs in which we detect an AGN in the mid-IR, the majority (2/3) of galaxies show [CII]/FIR >= 10^-3 typical of high 6.2um PAH EW sources, suggesting that most AGNs do not contribute significantly to the far-IR emission. We provide an empirical relation between the [CII]/FIR and the specific SFR (SSFR) for SF LIRGs. Finally, we present predictions for the starburst size based on the observed [CII] and far-IR luminosities which should be useful for comparing with results from future surveys of high-redshift galaxies with ALMA and CCAT.
[36]  [pdf] - 1171959
Tracing the Physical Conditions in Active Galactic Nuclei with Time-Dependent Chemistry
Comments: Accepted for publication at the Journal of Physical Chemistry A "Oka Festschrift: Celebrating 45 Years of Astrochemistry"
Submitted: 2013-06-11
We present an extension of the code ProDiMo that allows for a modeling of processes pertinent to active galactic nuclei and to an ambient chemistry that is time dependent. We present a proof-of-concept and focus on a few astrophysically relevant species, e.g., H+, H2+ and H3+; C+ and N+; C and O; CO and H2O; OH+, H2O+ and H3O+; HCN and HCO+. We find that the freeze-out of water is strongly suppressed and that this affects the bulk of the oxygen and carbon chemistry occurring in AGN. The commonly used AGN tracer HCN/HCO+ is strongly time-dependent, with ratios that vary over orders of magnitude for times longer than 10^4 years. Through ALMA observations this ratio can be used to probe how the narrow-line region evolves under large fluctuations in the SMBH accretion rate. Strong evolutionary trends, on time scales of 10^4-10^8 years, are also found in species such as H3O+, CO, and H2O. These reflect, respectively, time dependent effects in the ionization balance, the transient nature of the production of molecular gas, and the freeze-out/sublimation of water.
[37]  [pdf] - 1165797
The Response of Metal Rich Gas to X-Ray Irradiation from a Massive Black Hole at High Redshift: Proof of Concept
Comments: 27 pages, 11 figures. Resubmitted to ApJ
Submitted: 2013-04-07
Observational studies show that there is a strong link between the formation and evolution of galaxies and the growth of supermassive black holes (SMBH) at their centers. However, the underlying physics of this observed relation is poorly understood. In order to study the effects of X-ray radiation on the surroundings of the black hole, we implement X-ray Dominated Region (XDR) physics into Enzo and use the radiation transport module Moray to calculate the radiative transfer for a polychromatic spectrum. In this work, we investigate the effects of X-ray irradiation, produced by a central massive black hole (MBH) with a mass of M = 5x10^4 M_(solar), on ambient gas with solar and zero metallicity. We find that in the solar metallicity case, due to high opacity of the metals, the energy deposition rate in the central region (< 20 pc) is high and hence the temperatures in the center are on the order of 10^(5-7) K. Moreover, due to the cooling ability and high intrinsic opacity of solar metallicity gas, column densities of 10^(24) cm^(-2) are reached at a radius of 20 pc from the MBH. These column densities are about 3 orders of magnitudes higher than in the zero metallicity case. Furthermore, in the zero metallicity case an X-ray induced H II region is formed already after 5.8 Myr. This causes a significant outflow of gas (~8x10^6 M_(solar) from the central region, with the gas reaching outflow velocities up to ~100 km s^(-1). At later times, ~23 Myr after we insert the MBH, we find that the solar metallicity case also develops an X-ray induced H II region, but delayed by ~17 Myr.
[38]  [pdf] - 1159548
Water vapor in the protoplanetary disk of DG Tau
Comments: 6 pages, 2 figures, 3 tables, letter accepted for publication on ApJ
Submitted: 2013-02-06
Water is key in the evolution of protoplanetary disks and the formation of comets and icy/water planets. While high excitation water lines originating in the hot inner disk have been detected in several T Tauri stars (TTSs), water vapor from the outer disk, where most of water ice reservoir is stored, was only reported in the closeby TTS TW Hya. We present spectrally resolved Herschel/HIFI observations of the young TTS DG Tau in the ortho- and para- water ground-state transitions at 557, 1113 GHz. The lines show a narrow double-peaked profile, consistent with an origin in the outer disk, and are ~19-26 times brighter than in TW Hya. In contrast, CO and [C II] lines are dominated by emission from the envelope/outflow, which makes H2O lines a unique tracer of the disk of DG Tau. Disk modeling with the thermo-chemical code ProDiMo indicates that the strong UV field, due to the young age and strong accretion of DG Tau, irradiates a disk upper layer at 10-90 AU from the star, heating it up to temperatures of 600 K and producing the observed bright water lines. The models suggest a disk mass of 0.015-0.1 Msun, consistent with the estimated minimum mass of the solar nebula before planet formation, and a water reservoir of ~1e2-1e3 Earth oceans in vapour, and ~100 times larger in the form of ice. Hence, this detection supports the scenario of ocean delivery on terrestrial planets by impact of icy bodies forming in the outer disk.
[39]  [pdf] - 1158120
Evidence for CO shock excitation in NGC 6240 from Herschel SPIRE spectroscopy
Comments: Accepted by Astrophysical Journal Letters
Submitted: 2012-11-28
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.
[40]  [pdf] - 1158122
The SCUBA-2 Cosmology Legacy Survey: blank-field number counts of 450um-selected galaxies and their contribution to the cosmic infrared background
Comments: 9 pages, 5 figures, submitted to MNRAS
Submitted: 2012-11-28
The first deep blank-field 450um map (1-sigma~1.3mJy) from the SCUBA-2 Cosmology Legacy Survey (S2CLS), conducted with the James Clerk Maxwell Telescope (JCMT) is presented. Our map covers 140 arcmin^2 of the COSMOS field, in the footprint of the HST CANDELS area. Using 60 submillimetre galaxies (SMGs) detected at >3.75-sigma, we evaluate the number counts of 450um-selected galaxies with flux densities S_450>5mJy. The 8-arcsec JCMT beam and high sensitivity of SCUBA-2 now make it possible to directly resolve a larger fraction of the cosmic infrared background (CIB, peaking at ~200um) into the individual galaxies responsible for its emission than has previously been possible at this wavelength. At S_450>5mJy we resolve (7.4[+/-]0.7)x10^-2 MJy/sr of the CIB at 450um (equivalent to 16[+/-]7% of the absolute brightness measured by COBE at this wavelength) into point sources. A further ~40% of the CIB can be recovered through a statistical stack of 24um emitters in this field, indicating that the majority (~60%) of the CIB at 450um is emitted by galaxies with S_450>2mJy. The average redshift of 450um emitters identified with an optical/near-infrared counterpart is estimated to be <z>=1.3, implying that the galaxies in the sample are in the ultraluminous class (L_IR~1.1x10^12 L_sun). If the galaxies contributing to the statistical stack lie at similar redshifts, then the majority of the CIB at 450um is emitted by galaxies in the LIRG class with L_IR>3.6x10^11 L_sun.
[41]  [pdf] - 1151095
FUV and X-ray irradiated protoplanetary disks: a grid of models II - Gas diagnostic line emission
Comments: 15 pages, 11 figures. Accepted by A&A
Submitted: 2012-09-04
Most of the mass in protoplanetary disks is in the form of gas. The study of the gas and its diagnostics is of fundamental importance in order to achieve a detailed description of the thermal and chemical structure of the disk. The radiation from the central star (from optical to X-ray wavelengths) and viscous accretion are the main source of energy and dominates the disk physics and chemistry in its early stages. This is the environment in which the first phases of planet formation will proceed. We investigate how stellar and disk parameters impact the fine-structure cooling lines [NeII], [ArII], [OI], [CII] and H2O rotational lines in the disk. These lines are potentially powerful diagnostics of the disk structure and their modelling permits a thorough interpretation of the observations carried out with instrumental facilities such as Spitzer and Herschel. Following Aresu et al. (2011), we computed a grid of 240 disk models, in which the X-ray luminosity, UV-excess luminosity, minimum dust grain size, dust size distribution power law and surface density distribution power law, are systematically varied. We solve self-consistently for the disk vertical hydrostatic structure in every model and apply detailed line radiative transfer to calculate line fluxes and profiles for a series of well known mid- and far-infrared cooling lines. The [OI] 63 micron line flux increases with increasing FUV luminosity when Lx < 1e30 erg/s, and with increasing X-ray luminosity when LX > 1e30 erg/s. [CII] 157 micron is mainly driven by FUV luminosity via C+ production, X-rays affect the line flux to a lesser extent. [NeII] 12.8 micron correlates with X-rays; the line profile emitted from the disk atmosphere shows a double-peaked component, caused by emission in the static disk atmosphere, next to a high velocity double-peaked component, caused by emission in the very inner rim. (abridged)
[42]  [pdf] - 1150878
FUV and X-ray irradiated protoplanetary disks: a grid of models I. The disk structure
Comments: 47 pages, accepted by Astronomy and Astrophysics, a high resolution version of the paper is located at
Submitted: 2012-08-24
Context. Planets are thought to eventually form from the mostly gaseous (~99% of the mass) disks around young stars. The density structure and chemical composition of protoplanetary disks are affected by the incident radiation field at optical, FUV, and X-ray wavelengths, as well as by the dust properties. Aims. The effect of FUV and X-rays on the disk structure and the gas chemical composition are investigated. This work forms the basis of a second paper, which discusses the impact on diagnostic lines of, e.g., C+, O, H2O, and Ne+ observed with facilities such as Spitzer and Herschel. Methods. A grid of 240 models is computed in which the X-ray and FUV luminosity, minimum grain size, dust size distribution, and surface density distribution are varied in a systematic way. The hydrostatic structure and the thermo-chemical structure are calculated using ProDiMo. Results. The abundance structure of neutral oxygen is stable to changes in the X-ray and FUV luminosity, and the emission lines will thus be useful tracers of the disk mass and temperature. The C+ abundance distribution is sensitive to both X-rays and FUV. The radial column density profile shows two peaks, one at the inner rim and a second one at a radius r=5-10 AU. Ne+ and other heavy elements have a very strong response to X-rays, and the column density in the inner disk increases by two orders of magnitude from the lowest (LX = 1e29 erg/s) to the highest considered X-ray flux (LX = 1e32 erg/s). FUV confines the Ne+ ionized region to areas closer to the star at low X-ray luminosities (LX = 1e29 erg/s). H2O abundances are enhanced by X-rays due to higher temperatures in the inner disk and higher ionization fractions in the outer disk. The line fluxes and profiles are affected by the effects on these species, thus providing diagnostic value in the study of FUV and X-ray irradiated disks around T Tauri stars. (abridged)
[43]  [pdf] - 1117486
Diagnostics of the Molecular Component of PDRs with Mechanical Heating
Comments: 26 pages, 14 figures in the text and 13 figures as supplementary material. Accepted for publication in A&A
Submitted: 2012-03-22
Context. Multitransition CO observations of galaxy centers have revealed that significant fractions of the dense circumnuclear gas have high kinetic temperatures, which are hard to explain by pure photon excitation, but may be caused by dissipation of turbulent energy. Aims. We aim to determine to what extent mechanical heating should be taken into account while modelling PDRs. To this end, the effect of dissipated turbulence on the thermal and chemical properties of PDRs is explored. Methods. Clouds are modelled as 1D semi-infinite slabs whose thermal and chemical equilibrium is solved for using the Leiden PDR-XDR code. Results. In a steady-state treatment, mechanical heating seems to play an important role in determining the kinetic temperature of the gas in molecular clouds. Particularly in high-energy environments such as starburst galaxies and galaxy centers, model gas temperatures are underestimated by at least a factor of two if mechanical heating is ignored. The models also show that CO, HCN and H2 O column densities increase as a function of mechanical heating. The HNC/HCN integrated column density ratio shows a decrease by a factor of at least two in high density regions with n \sim 105 cm-3, whereas that of HCN/HCO+ shows a strong dependence on mechanical heating for this same density range, with boosts of up to three orders of magnitude. Conclusions. The effects of mechanical heating cannot be ignored in studies of the molecular gas excitation whenever the ratio of the star formation rate to the gas density is close to, or exceeds, 7 \times 10-6 M yr-1 cm4.5 . If mechanical heating is not included, predicted column densities are underestimated, sometimes even by a few orders of magnitude. As a lower bound to its importance, we determined that it has non-negligible effects already when mechanical heating is as little as 1% of the UV heating in a PDR.
[44]  [pdf] - 1091563
Gas modelling in the disc of HD 163296
Comments: 20 pages, 15 figures. Accepted for publication in A&A
Submitted: 2011-11-10
We present detailed model fits to observations of the disc around the Herbig Ae star HD 163296. This well-studied object has an age of ~ 4 Myr, with evidence of a circumstellar disc extending out to ~ 540AU. We use the radiation thermo-chemical disc code ProDiMo to model the gas and dust in the circumstellar disc of HD 163296, and attempt to determine the disc properties by fitting to observational line and continuum data. These include new Herschel/PACS observations obtained as part of the open-time key program GASPS (Gas in Protoplanetary Systems), consisting of a detection of the [OI]63mic line and upper limits for several other far infrared lines. We complement this with continuum data and ground-based observations of the 12CO 3-2, 2-1 and 13CO J=1-0 line transitions, as well as the H2 S(1) transition. We explore the effects of stellar ultraviolet variability and dust settling on the line emission, and on the derived disc properties. Our fitting efforts lead to derived gas/dust ratios in the range 9-100, depending on the assumptions made. We note that the line fluxes are sensitive in general to the degree of dust settling in the disc, with an increase in line flux for settled models. This is most pronounced in lines which are formed in the warm gas in the inner disc, but the low excitation molecular lines are also affected. This has serious implications for attempts to derive the disc gas mass from line observations. We derive fractional PAH abundances between 0.007 and 0.04 relative to ISM levels. Using a stellar and UV excess input spectrum based on a detailed analysis of observations, we find that the all observations are consistent with the previously assumed disc geometry.
[45]  [pdf] - 1091525
Enhanced H2O formation through dust grain chemistry in X-ray exposed environments
Comments: 13 pages, 8 figures, accepted by A&A
Submitted: 2011-11-09
The ULIRG Mrk 231 exhibits very strong water rotational lines between \lambda = 200-670\mu m, comparable to the strength of the CO rotational lines. High redshift quasars also show similar CO and H2O line properties, while starburst galaxies, such as M82, lack these very strong H2O lines in the same wavelength range, but do show strong CO lines. We explore the possibility of enhancing the gas phase H2O abundance in X-ray exposed environments, using bare interstellar carbonaceous dust grains as a catalyst. Cloud-cloud collisions cause C and J shocks, and strip the grains of their ice layers. The internal UV field created by X-rays from the accreting black hole does not allow to reform the ice. We determine formation rates of both OH and H2O on dust grains, having temperature T_dust=10-60 K, using both Monte Carlo as well as rate equation method simulations. The acquired formation rates are added to our X-ray chemistry code, that allows us to calculate the thermal and chemical structure of the interstellar medium near an active galactic nucleus. We derive analytic expressions for the formation of OH and H2O on bare dust grains as a catalyst. Oxygen atoms arriving on the dust are released into the gas phase under the form of OH and H2O. The efficiencies of this conversion due to the chemistry occurring on dust are of order 30 percent for oxygen converted into OH and 60 percent for oxygen converted into H_2O between T_dust=15-40 K. At higher temperatures, the efficiencies rapidly decline. When the gas is mostly atomic, molecule formation on dust is dominant over the gas-phase route, which is then quenched by the low H2 abundance. Here, it is possible to enhance the warm (T> 200 K) water abundance by an order of magnitude in X-ray exposed environments. This helps to explain the observed bright water lines in nearby and high-redshift ULIRGs and Quasars.
[46]  [pdf] - 1077512
Water vapor emission reveals a highly obscured, star forming nuclear region in the QSO host galaxy APM08279+5255 at z=3.9
Comments: accepted for publication in Astrophysical Journal Letters; typos fixed
Submitted: 2011-06-23, last modified: 2011-09-27
We present the detection of four rotational emission lines of water vapor, from energy levels Eu/k= 101 - 454 K, in the gravitationally lensed z=3.9 QSO host galaxy APM08279+5255. While the lowest H2O lines are collisionally excited in clumps of warm, dense gas (density of hydrogen nuclei n_H=(3.1 +/- 1.2) x 10^6 cm^-3, gas temperature T_g ~ 105 +/- 21 K), we find that the excitation of the higher lines is dominated by the intense local infrared radiation field. Since only collisionally excited emission contributes to gas cooling, we conclude that H2O is not a significant coolant of the warm molecular gas. Our excitation model requires the radiatively excited gas to be located in an extended region of high 100 micron opacity (tau_100 = 0.9 +/- 0.2). Locally, such extended infrared-opaque regions are found only in the nuclei of ultraluminous infrared galaxies. We propose a model where the infrared-opaque circumnuclear cloud, which is penetrated by the X-ray radiation field of the QSO nucleus, contains clumps of massive star formation where the H2O emission originates. The radiation pressure from the intense local infrared radiation field exceeds the thermal gas pressure by about an order of magnitude, suggesting close to Eddington-limited star formation in these clumps.
[47]  [pdf] - 1076700
X-ray Ionization of Heavy Elements Applied to Protoplanetary Disks
Comments: Accepted for publication in ApJ
Submitted: 2011-05-17
The consequences of the Auger effect on the population of heavy element ions are analyzed for the case of relatively cool gas irradiated by keV X-rays, with intended applications to the accretion disks of young stellar ob jects. Highly charged ions are rapidly reduced to the doubly-charged state in neutral gas, so the aim here is to derive the production rates for these singly- and doubly-charged ions and to specify their transformation by recombination, charge transfer, and molecular reactions. The theory is illustrated by calculations of the abundances of eleven of the most cosmically abundant heavy elements in a model of a typical T Tauri star disk. The physical properties of the gas are determined with an X-ray irradiated thermal-chemical model, which shows that the disk atmosphere consists of a hot atmosphere overlaying the mainly cool body of the disk. There is a warm transition layer where hydrogen, carbon and oxygen are transformed from atomic to molecular species and the ionization drops by several orders of magnitude. The doubly charged ions are then largely confined to the hot outer layers of the disk. The nature of the dominant ions below the transition depends sensitively on the poorly constrained abundances of the heavy elements. Observational consequences and connections with the active layers of the magneto-rotational instability are briefly discussed.
[48]  [pdf] - 1042860
H3O+ Line Emission from Starbursts and AGNs
Comments: 8 pages, 2 figures, accepted for publication in Astronomy and Astrophysics journal
Submitted: 2011-01-04
We use the H3O+ molecule to investigate the impact of starburst and AGN activity on the chemistry of the molecular interstellar medium. Using the JCMT, we have observed the 3+_2 - 2-_2 364 GHz line of p-H3O+ towards the centers of seven active galaxies. We have detected p-H3O+ towards IC342, NGC253, NGC1068, NGC4418, and NGC6240. Upper limits were obtained for IRAS15250 and Arp299. We find large H3O+ abundances (N(H3O+)/N(H2)>10^{-8}) in all detected galaxies apart from in IC342 where it is about one order of magnitude lower. We note, however, that uncertainties in N(H3O+) may be significant due to lack of definite information on source size and excitation. We furthermore compare the derived N(H3O+) with N(HCO+) and find that the H3O+ to HCO+ column density ratio is large in NGC1068 (24), moderate in NGC4418 and NGC253 (4-5), slightly less than unity in NGC6240 (0.7) and lowest in IC342 (0.2-0.6). We compare our results with models of X-ray and photon dominated regions (XDRs and PDRs). For IC342 we find that a starburst PDR chemistry can explain the observed H3O+ abundance. For the other galaxies, the large H3O+ columns are generally consistent with XDR models. In particular for NGC1068 the elevated N(H3O+)/N(HCO+) ratio suggests a low column density XDR. For NGC4418 however, large HC3N abundances are inconsistent with the XDR interpretation. An alternative possibility is that H3O+ forms through H2O evaporating off dust grains and reacting with HCO+ in warm, dense gas. This scenario could also potentially fit the results for NGC253. Further studies of the excitation and distribution of H3O+ - as well as Herschel observations of water abundances - will help to further constrain the models.
[49]  [pdf] - 1042190
X-ray impact on the protoplanetary disks around T Tauri stars
Comments: 6 pages, accepted for publication in Astronomy and Astrophysics (Interstellar and circumstellar matter section, date of acceptance is 24/11/2010)
Submitted: 2010-11-25, last modified: 2010-11-29
Context: T Tauri stars have X-ray luminosities ranging from L_X = 10^28-10^32 erg/s. These luminosities are similar to UV luminosities (L_UV 10^30-10^31 erg/s) and therefore X-rays are expected to affect the physics and chemistry of the upper layers of their surrounding protoplanetary disks. Aim: The effects and importance of X-rays on the chemical and hydrostatic structure of protoplanetary disks are investigated, species tracing X-ray irradiation (for L_X >= 10^29 erg/s) are identified and predictions for [OI], [CII] and [NII] fine structure line fluxes are provided. Methods: We have implemented X-ray physics and chemistry into the chemo-physical disk code ProDiMo. We include Coulomb heating and H2 ionization as heating processes and primary and secondary ionization due to X-rays in the chemistry. Results: X-rays heat up the gas causing it to expand in the optically thin surface layers. Neutral molecular species are not much affected in their abundance and spatial distribution, but charged species such as N+, OH+, H2O+ and H3O+ show enhanced abundances in the disk surface. Conclusions: Coulomb heating by X-rays changes the vertical structure of the disk, yielding temperatures of ~ 8000 K out to distances of 50 AU. The chemical structure is altered by the high electron abundance in the gas in the disk surface, causing an efficient ion-molecule chemistry. The products of this, OH+, H2O+ and H3O+, are of great interest for observations of low-mass young stellar objects with the Herschel Space Observatory. [OI] (at 63 and 145 mic) and [CII] (at 158 mic) fine structure emission are only affected for L_X > 10^30 erg/s.
[50]  [pdf] - 1041065
Star Formation in Extreme Environments: The Effects of Cosmic Rays and Mechanical Heating
Comments: 13 pages, 8 figures, A&A accepted
Submitted: 2010-10-05
Context: Molecular data of extreme environments, such as Arp 220, but also NGC 253, show evidence for extremely high cosmic ray (CR) rates (10^3-10^4 * Milky Way) and mechanical heating from supernova driven turbulence. Aims: The consequences of high CR rates and mechanical heating on the chemistry in clouds are explored. 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 to 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. HCN and H2O abundances are boosted for very high mechanical heating rates, while ionized species are relatively unaffected. OH+ and H2O+ are enhanced for very high CR rates zeta > 5 * 10^-14 s^-1. A combination of OH+, OH, H2O+, H2O, and H3O+ trace the CR rates, and are able to distinguish between enhanced cosmic rays and X-rays.
[51]  [pdf] - 1033326
Excitation of the molecular gas in the nuclear region of M82
Comments: Accepted for publication in A&A Letters
Submitted: 2010-06-25, last modified: 2010-07-05
We present high resolution HIFI spectroscopy of the nucleus of the archetypical starburst galaxy M82. Six 12CO lines, 2 13CO lines and 4 fine-structure lines are detected. Besides showing the effects of the overall velocity structure of the nuclear region, the line profiles also indicate the presence of multiple components with different optical depths, temperatures and densities in the observing beam. The data have been interpreted using a grid of PDR models. It is found that the majority of the molecular gas is in low density (n=10^3.5 cm^-3) clouds, with column densities of N_H=10^21.5 cm^-2 and a relatively low UV radiation field (GO = 10^2). The remaining gas is predominantly found in clouds with higher densities (n=10^5 cm^-3) and radiation fields (GO = 10^2.75), but somewhat lower column densities (N_H=10^21.2 cm^-2). The highest J CO lines are dominated by a small (1% relative surface filling) component, with an even higher density (n=10^6 cm^-3) and UV field (GO = 10^3.25). These results show the strength of multi-component modeling for the interpretation of the integrated properties of galaxies.
[52]  [pdf] - 1033229
A Spitzer survey of mid-infrared molecular emission from protoplanetary disks I: Detection rates
Comments: 19 pages, accepted for publication in ApJ
Submitted: 2010-06-21
We present a Spitzer InfraRed Spectrometer search for 10-36 micron molecular emission from a large sample of protoplanetary disks, including lines from H2O, OH, C2H2, HCN and CO2. This paper describes the sample and data processing and derives the detection rate of mid-infrared molecular emission as a function of stellar mass. The sample covers a range of spectral type from early M to A, and is supplemented by archival spectra of disks around A and B stars. It is drawn from a variety of nearby star forming regions, including Ophiuchus, Lupus and Chamaeleon. In total, we identify 22 T Tauri stars with strong mid-infrared H2O emission. Integrated water line luminosities, where water vapor is detected, range from 5x10^-4 to 9x10^-3 Lsun, likely making water the dominant line coolant of inner disk surfaces in classical T Tauri stars. None of the 5 transitional disks in the sample show detectable gaseous molecular emission with Spitzer upper limits at the 1% level in terms of line-to-continuum ratios (apart from H2). We find a strong dependence on detection rate with spectral type; no disks around our sample of 25 A and B stars were found to exhibit water emission, down to 1-2% line-to-continuum ratios, in the mid-infrared, while almost 2/3 of the disks around K stars show sufficiently intense water emission to be detected by Spitzer. Some Herbig Ae/Be stars show tentative H2O/OH emission features beyond 20 micron at the 1-2 level, however, and one of them shows CO2 in emission. We argue that the observed differences between T Tauri disks and Herbig Ae/Be disks is due to a difference in excitation and/or chemistry depending on spectral type and suggest that photochemistry may be playing an important role in the observable characteristics of mid-infrared molecular line emission from protoplanetary disks.
[53]  [pdf] - 1032623
Herschel observations of water vapour in Markarian 231
Comments: 5 pages, 3 figures
Submitted: 2010-05-20
The Ultra Luminous InfraRed Galaxy Mrk 231 reveals up to seven rotational lines of water (H2O) in emission, including a very high-lying (E_{upper}=640 K) line detected at a 4sigma level, within the Herschel/SPIRE wavelength range, whereas PACS observations show one H2O line at 78 microns in absorption, as found for other H2O lines previously detected by ISO. The absorption/emission dichotomy is caused by the pumping of the rotational levels by far-infrared radiation emitted by dust, and subsequent relaxation through lines at longer wavelengths, which allows us to estimate both the column density of H2O and the general characteristics of the underlying far-infrared continuum source. Radiative transfer models including excitation through both absorption of far-infrared radiation emitted by dust and collisions are used to calculate the equilibrium level populations of H2O and the corresponding line fluxes. The highest-lying H2O lines detected in emission, with levels at 300-640 K above the ground state, indicate that the source of far-infrared radiation responsible for the pumping is compact (radius=110-180 pc) and warm (T_{dust}=85-95 K), accounting for at least 45% of the bolometric luminosity. The high column density, N(H2O)~5x10^{17} cm^{-2}, found in this nuclear component, is most probably the consequence of shocks/cosmic rays, an XDR chemistry, and/or an "undepleted chemistry" where grain mantles are evaporated. A more extended region, presumably the inner region of the 1-kpc disk observed in other molecular species, could contribute to the flux observed in low-lying H2O lines through dense hot cores, and/or shocks. The H2O 78 micron line observed with PACS shows hints of a blue-shifted wing seen in absorption, possibly indicating the occurrence of H2O in the prominent outflow detected in OH (Fischer et al., this volume).
[54]  [pdf] - 1032599
First results of the Herschel Key Program 'Dust, Ice and Gas in Time': Dust and Gas Spectroscopy of HD 100546
Submitted: 2010-05-19
We present far-infrared spectroscopic observations, taken with the Photodetector Array Camera and Spectrometer (PACS) on the Herschel Space Observatory, of the protoplanetary disk around the pre-main-sequence star HD 100546. These observations are the first within the DIGIT Herschel key program, which aims to follow the evolution of dust, ice, and gas from young stellar objects still embedded in their parental molecular cloud core, through the final pre-main-sequence phases when the circumstellar disks are dissipated. Our aim is to improve the constraints on temperature and chemical composition of the crystalline olivines in the disk of HD 100546 and to give an inventory of the gas lines present in its far-infrared spectrum. The 69 \mu\m feature is analyzed in terms of position and shape to derive the dust temperature and composition. Furthermore, we detected 32 emission lines from five gaseous species and measured their line fluxes. The 69 \mu\m emission comes either from dust grains with ~70 K at radii larger than 50 AU, as suggested by blackbody fitting, or it arises from ~200 K dust at ~13 AU, close to the midplane, as supported by radiative transfer models. We also conclude that the forsterite crystals have few defects and contain at most a few percent iron by mass. Forbidden line emission from [CII] at 157 \mu\m and [OI] at 63 and 145 \mu\m, most likely due to photodissociation by stellar photons, is detected. Furthermore, five H2O and several OH lines are detected. We also found high-J rotational transition lines of CO, with rotational temperatures of ~300 K for the transitions up to J=22-21 and T~800 K for higher transitions.
[55]  [pdf] - 1032510
Black hole accretion and star formation as drivers of gas excitation and chemistry in Mrk231
Comments: 5 pages, 2 figures, accepted for publication in Astronomy & Astrophysics Special Issue on Herschel first results
Submitted: 2010-05-17
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.
[56]  [pdf] - 1032479
Dust, Ice and Gas in Time (DIGIT) Herschel program first results: A full PACS-SED scan of the gas line emission in protostar DK Cha
Comments: 4 Page letter, To appear in A&A special issue on Herschel
Submitted: 2010-05-14
DK Cha is an intermediate-mass star in transition from an embedded configuration to a star plus disk stage. We aim to study the composition and energetics of the circumstellar material during this pivotal stage. Using the Range Scan mode of PACS on the Herschel Space Observatory, we obtained a spectrum of DK Cha from 55 to 210 micron as part of the DIGIT Key Program. Almost 50 molecular and atomic lines were detected, many more than the 7 lines detected in ISO-LWS. Nearly the entire ladder of CO from J=14-13 to 38-37 (E_u/k = 4080 K), water from levels as excited as E_u/k = 843 K, and OH lines up to E_u/k = 290 K were detected. The continuum emission in our PACS SED scan matches the flux expected from a model consisting of a star, a surrounding disk of 0.03 Solar mass, and an envelope of a similar mass, supporting the suggestion that the object is emerging from its main accretion stage. Molecular, atomic, and ionic emission lines in the far-infrared reveal the outflow's influence on the envelope. The inferred hot gas can be photon-heated, but some emission could be due to C-shocks in the walls of the outflow cavity.
[57]  [pdf] - 901629
Radiative transfer models of mid-infrared H2O lines in the Planet-forming Region of Circumstellar Disks
Comments: 12 pages, accepted for publication in ApJ
Submitted: 2009-09-04
The study of warm molecular gas in the inner regions of protoplanetary disks is of key importance for the study of planet formation and especially for the transport of H2O and organic molecules to the surfaces of rocky planets/satellites. Recent Spitzer observations have shown that the mid-infrared spectra of protoplanetary disks are covered in emission lines due to water and other molecules. Here, we present a non-LTE 2D radiative transfer model of water lines in the 10-36 mum range that can be used to constrain the abundance structure of water vapor, given an observed spectrum, and show that an assumption of local thermodynamic equilibrium (LTE) does not accurately estimate the physical conditions of the water vapor emission zones. By applying the model to published Spitzer spectra we find that: 1) most water lines are subthermally excited, 2) the gas-to-dust ratio must be one to two orders of magnitude higher than the canonical interstellar medium ratio of 100-200, and 3) the gas temperature must be higher than the dust temperature, and 4) the water vapor abundance in the disk surface must be truncated beyond ~ 1 AU. A low efficiency of water formation below ~ 300 K may naturally result in a lower water abundance beyond a certain radius. However, we find that chemistry, may not be sufficient to produce an abundance drop of many orders of magnitude and speculate that the depletion may also be caused by vertical turbulent diffusion of water vapor from the superheated surface to regions below the snow line, where the water can freeze out and be transported to the midplane as part of the general dust settling. Such a vertical cold finger effect is likely to be efficient due to the lack of a replenishment mechanism of large, water-ice coated dust grains to the disk surface.
[58]  [pdf] - 1017153
A new raytracer for modeling AU-scale imaging of lines from protoplanetary disks
Comments: 14 pages, accepted for publication in ApJ
Submitted: 2009-08-20
The material that formed the present-day Solar System originated in feeding zones in the inner Solar Nebula located at distances within ~20 AU from the Sun, known as the planet-forming zone. Meteoritic and cometary material contain abundant evidence for the presence of a rich and active chemistry in the planet-forming zone during the gas-rich phase of Solar System formation. It is a natural conjecture that analogs can be found amoung the zoo of protoplanetary disks around nearby young stars. The study of the chemistry and dynamics of planet formation requires: 1) tracers of dense gas at 100-1000 K and 2) imaging capabilities of such tracers with 5-100 (0.5-20 AU) milli-arcsec resolution, corresponding to the planet-forming zone at the distance of the closest star-forming regions. Recognizing that the rich infrared (2-200 micron) molecular spectrum recently discovered to be common in protoplanetary disks represents such a tracer, we present a new general raytracing code, RADLite, that is optimized for producing infrared line spectra and images from axisymmetric structures. RADLite can consistently deal with a wide range of velocity gradients, such as those typical for the inner regions of protoplanetary disks. The code is intended as a backend for chemical and excitation codes, and can rapidly produce spectra of thousands of lines for grids of models for comparison with observations. Such radiative transfer tools will be crucial for constraining both the structure and chemistry of planet-forming regions, including data from current infrared imaging spectrometers and extending to the Atacama Large Millimeter Array and the next generation of Extremely Large Telescopes, the James Webb Space Telescope and beyond.
[59]  [pdf] - 24702
Formation of Water in the Warm Atmospheres of Protoplanetary Disks
Submitted: 2009-05-27
The gas-phase chemistry of water in protoplanetary disks is analyzed with a model based on X-ray heating and ionization of the disk atmosphere. Several uncertain processes appear to play critical roles in generating the column densities of warm water that are detected from disks at infrared wavelengths. The dominant factors are the reactions that form molecular hydrogen, including formation on warm grains, and the ionization and heating of the atmosphere. All of these can work together to produce a region of high water abundances in the molecular transition layer of the inner disk atmosphere, where atoms are transformed into molecules, the temperature drops from thousands to hundreds of Kelvins, and the ionization begins to be dominated by the heavy elements. Grain formation of molecular hydrogen and mechanical heating of the atmosphere can play important roles in this region and directly affect the amount of warm water in protoplanetary disk atmospheres. Thus it may be possible to account for the existing measurements of water emission from Tauri disks without invoking transport of water from cooler to warmer regions. The hydroxyl radical OH is under-abundant in this model of disk atmospheres and requires consideration of additional production and excitation processes.
[60]  [pdf] - 23402
Training the Next Generation of Astronomers
Comments: 9 pp + cover page; position paper submitted to the Astro2010 Decadal Survey Subcommittee on the State of the Profession; small typographical corrections made; submitted version at
Submitted: 2009-04-16
While both society and astronomy have evolved greatly over the past fifty years, the academic institutions and incentives that shape our field have remained largely stagnant. As a result, the astronomical community is faced with several major challenges, including: (1) the training that we provide does not align with the skills that future astronomers will need, (2) the postdoctoral phase is becoming increasingly demanding and demoralizing, and (3) our jobs are increasingly unfriendly to families with children. Solving these problems will require conscious engineering of our profession. Fortunately, this Decadal Review offers the opportunity to revise outmoded practices to be more effective and equitable. The highest priority of the Subcommittee on the State of the Profession should be to recommend specific, funded activities that will ensure the field meets the challenges we describe.
[61]  [pdf] - 23239
High Resolution Spectroscopy of [NeII] Emission from AA Tau and GM Aur
Comments: Accepted for publication in ApJ
Submitted: 2009-04-09
We present high resolution (R=80,000) spectroscopy of [NeII] emission from two young stars, GM Aur and AA Tau, which have moderate to high inclinations. The emission from both sources appears centered near the stellar velocity and is broader than the [NeII] emission measured previously for the face-on disk system TW Hya. These properties are consistent with a disk origin for the [NeII] emission we detect, with disk rotation (rather than photoevaporation or turbulence in a hot disk atmosphere) playing the dominant role in the origin of the line width. In the non-face-on systems, the [NeII] emission is narrower than the CO fundamental emission from the same sources. If the widths of both diagnostics are dominated by Keplerian rotation, this suggests that the [NeII] emission arises from larger disk radii on average than does the CO emission. The equivalent width of the [NeII] emission we detect is less than that of the spectrally unresolved [NeII] feature in the Spitzer spectra of the same sources. Variability in the [NeII] emission or the mid-infrared continuum, a spatially extended [NeII] component, or a very (spectrally) broad [NeII] component might account for the difference in the equivalent widths.
[62]  [pdf] - 17292
Rotational Line Emission from Water in Protoplanetary Disks
Comments: 5 pages, 4 figures. Accepted by ApJL
Submitted: 2008-10-09
Circumstellar disks provide the material reservoir for the growth of young stars and for planet formation. We combine a high-level radiative transfer program with a thermal-chemical model of a typical T Tauri star disk to investigate the diagnostic potential of the far-infrared lines of water for probing disk structure. We discuss the observability of pure rotational H2O lines with the Herschel Space Observatory, specifically the residual gas where water is mainly frozen out. We find that measuring both the line profile of the ground 110-101 ortho-H2O transition and the ratio of this line to the 312-303 and 221-212 line can provide information on the gas phase water between 5-100 AU, but not on the snow line which is expected to occur at smaller radii.
[63]  [pdf] - 14408
Mechanical feedback in the molecular ISM of luminous IR galaxies
Comments: 4 pages, 3 figures, published as Letter to the Editor in A&A (see
Submitted: 2008-07-10, last modified: 2008-08-20
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 observedsystems 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.
[64]  [pdf] - 10997
On the Detection of High Redshift Black Holes with ALMA through CO and H2 Emission
Comments: 12 pages, 3 figures, ApJL accepted; minor typographical changes
Submitted: 2008-03-17, last modified: 2008-03-31
Many present-day galaxies are known to harbor supermassive, >10^6 Mo, black holes. These central black holes must have grown through accretion from less massive seeds in the early universe. The molecules CO and H2 can be used to trace this young population of accreting massive black holes through the X-ray irradiation of ambient gas. The X-rays drive a low-metallicity ion-molecule chemistry that leads to the formation and excitation of CO and H2 in 100<T<1,000 K gas. H2 traces very low metallicity gas, ~10^-3 solar or less, while some pollution by metals, ~10^-2 solar or more, must have taken place to form CO. Strong CO J>15 and H2 S(0) & S(1) emission is found that allows one to constrain ambient conditions. Comparable line strengths cannot be produced by FUV or cosmic ray irradiation. Weak, but perhaps detectable, H3+ (2,2) - (1,1) emission is found and discussed. The models predict that black hole masses larger than 10^5 Mo can be detected with ALMA, over a redshift range of 5-20, provided that the black holes radiate close to Eddington.
[65]  [pdf] - 8815
Extragalactic H3O+: Some Consequences
Comments: Five pages, one figure; contribution to proceedings of conference "Far-infrared observations of the interstellar medium", December 2007, Bad Honnef
Submitted: 2008-01-09
We discuss some implications of our recent detection of extragalactic H3O+: the location of the gas in M82, the origin of energetic radiation in M82, and the possible feedback effects of star formation on the cosmic ray flux in galaxies.
[66]  [pdf] - 7627
Atomic Diagnostics of X-ray Irradiated Protoplanetary Disks
Comments: Accepted by ApJ
Submitted: 2007-12-01
We study atomic line diagnostics of the inner regions of protoplanetary disks with our model of X-ray irradiated disk atmospheres which was previously used to predict observable levels of the NeII and NeIII fine-structure transitions at 12.81 and 15.55mum. We extend the X-ray ionization theory to sulfur and calculate the fraction of sulfur in S, S+, S2+ and sulfur molecules. For the D'Alessio generic T Tauri star disk, we find that the SI fine-structure line at 25.55mum is below the detection level of the Spitzer Infrared Spectrometer (IRS), in large part due to X-ray ionization of atomic S at the top of the atmosphere and to its incorporation into molecules close to the mid-plane. We predict that observable fluxes of the SII 6718/6732AA forbidden transitions are produced in the upper atmosphere at somewhat shallower depths and smaller radii than the neon fine-structure lines. This and other forbidden line transitions, such as the OI 6300/6363AA and the CI 9826/9852AA lines, serve as complementary diagnostics of X-ray irradiated disk atmospheres. We have also analyzed the potential role of the low-excitation fine-structure lines of CI, CII, and OI, which should be observable by SOFIA and Herschel.
[67]  [pdf] - 7007
Detection of extragalactic H3O+
Comments: Accepted by A&A Letters; 4-5 pages depending on paper format; two b/w figures
Submitted: 2007-11-14
The H3O+ molecule probes the oxygen chemistry and the ionization rate of dense circumnuclear gas in galaxies. In particular, recent H3O+ observations show variations in the cosmic-ray ionization rate by factors of $>$10 within our Galaxy. Using the JCMT, we have observed the 364 GHz line of p-H3O+ in the centers of M82 and Arp 220. In Arp 220, the line profile suggests that the emission originates in the Western nucleus. In M82, both the eastern molecular peak and the circumnuclear region contribute to the emission. The derived column densities, abundances, and H3O+ / H2O ratios indicate ionization rates similar to or even exceeding that in the Galactic Center. Model calculations of the chemistry of irradiated molecular gas indicate a likely origin of this high ionization rate in the extended, evolved starburst of M82. In contrast, irradiation by X-rays from the AGN disk is the most likely model for Arp 220.
[68]  [pdf] - 5508
The Dense Gas in the Central Kiloparsec of NGC 6946
Comments: ApJ accepted. 11 pages, 13 figures
Submitted: 2007-09-30
We present observations of the HCN and HCO+ J=1-0 transitions in the center of the nearby spiral galaxy NGC 6946 made with the BIMA and CARMA interferometers. Using the BIMA SONG CO map, we investigate the change in the I_HCN/I_CO and I_ HCO/I_CO integrated intensity ratios as a function of radius in the central kiloparsec of the galaxy, and find that they are strongly concentrated at the center. We use the 2MASS K_S band image to find the stellar surface density, and then construct a map of the hydrostatic midplane pressure. We apply a PDR model to the observed I_HCN/I_HCO+ integrated intensity ratio to calculate the number density of molecular hydrogen in the dense gas tracer emitting region, and find that it is roughly constant at 10^5 cm^-3 across our map. We explore two hypotheses for the distribution of the dense gas. If the HCN and HCO+ emission comes from self-gravitating density peaks inside of a less dense gas distribution, there is a linear proportionality between the internal velocity dispersion of the dense gas and the size of the density peak. Alternatively, the HCN and HCO+ emission could come from dense gas homogeneously distributed throughout the center and bound by ambient pressure, similar to what is observed toward the center of the Milky Way. We find both of these hypotheses to be plausible. We fit the relationships between I_HCN, I_HCO+, and I_CO. Correlations between the hydrostatic midplane pressure and I_HCN and I_HCO+ are demonstrated, and power law fits are provided. We confirm the validity of a relation found by Blitz & Rosolowsky (2006) between pressure and the molecular to atomic gas ratio in the high hydrostatic midplane pressure regime (10^6-10^8 cm^-3 K).
[69]  [pdf] - 5205
The irradiated ISM of ULIRGs
Comments: 5 pages, 6 figures, to appear in: IAU Symposium 242 Astrophysical Masers and their Environments
Submitted: 2007-09-21
The nuclei of ULIRGs harbor massive young stars, an accreting central black hole, or both. Results are presented for molecular gas that is exposed to X-rays (1-100 keV, XDRs) and far-ultraviolet radiation (6-13.6 eV, PDRs). Attention is paid to species like HCO+, HCN, HNC, OH, H2O and CO. Line ratios of HCN/HCO+ and HNC/HCN discriminate between PDRs and XDRs. Very high J (>10) CO lines, observable with HIFI/Herschel, discriminate very well between XDRs and PDRs. In XDRs, it is easy to produce large abundances of warm (T>100 K) H2O and OH. In PDRs, only OH is produced similarly well.
[70]  [pdf] - 1779
CO+ in M 82: A Consequence of Irradiation by X-rays
Comments: added discussion on more recent X-ray observations
Submitted: 2007-05-31, last modified: 2007-06-02
Based on its strong CO+ emission it is argued that the M 82 star-burst galaxy is exposed to a combination of FUV and X-ray radiation. The latter is likely to be the result of the star-burst superwind, which leads to diffuse thermal emission at ~0.7 keV, and a compact hard, 2-10 keV, source (but not an AGN). Although a photon-dominated region (FUV) component is clearly present in the nucleus of M 82, and capable of forming CO+, only X-ray irradiated gas of density 10^3-10^5 cm^-3 can reproduce the large, ~(1-4)x10^13 cm^-2, columns of CO+ that are observed toward the proto-typical star-burst M 82. The total X-ray luminosity produced by M 82 is weak, ~10^41 erg s^-1, but this is sufficient to drive the formation of CO+.
[71]  [pdf] - 338144
A PDR-Code Comparison Study
Comments: 28 pages, 16 figures, see all benchmark data under:
Submitted: 2007-02-08
We present a comparison between independent computer codes, modeling the physics and chemistry of interstellar photon dominated regions (PDRs). Our goal was to understand the mutual differences in the PDR codes and their effects on the physical and chemical structure of the model clouds, and to converge the output of different codes to a common solution. A number of benchmark models have been created, covering low and high gas densities and far ultraviolet intensities. The benchmark models were computed in two ways: one set assuming constant temperatures, thus testing the consistency of the chemical network and photo-processes, and a second set determining the temperature selfconsistently. We investigated the impact of PDR geometry and agreed on the comparison of results from spherical and plane-parallel PDR models. We identified a number of key processes governing the chemical network which have been treated differently in the various codes, and defined a proper common treatment. We established a comprehensive set of reference models for ongoing and future PDR model bench-marking and were able to increase the agreement in model predictions for all benchmark models significantly.
[72]  [pdf] - 85784
Diagnostics of irradiated dense gas in galaxy nuclei. II. A grid of XDR and PDR models
Comments: 34 pages, 47 figures, accepted by A&A, a version with all figures is available at
Submitted: 2006-10-12, last modified: 2006-10-13
\abridged 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 the sources of radiation, numerical models are constructed. These models determine the thermal and chemical balance of molecular gas that is exposed to X-rays and far-ultraviolet radiation, as a function of depth. We present a grid of XDR and PDR models that span ranges in density, irradiation and column density. We find that the fine-structure line ratios of e.g. [SiII] 35mum/[CII] 158 mum are larger in XDRs than in PDRs, for a given density, column and irradiation strength. We find that the line ratios HCN/HCO+ and HNC/HCN, as well as the column density ratio CN/HCN, discriminate between PDRs and XDRs. The HCN/HCO+ 1-0 ratio is <1 (>1) for XDRs (PDRs) if the density exceeds 10^5 cm^-3 and if the column density is >10^23 cm^-2. For columns <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. 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 columns >10^{23} cm^{-2} and densities <10^4 cm^-3. 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. Lines with N>=10, like CO(16-15) and CO(10-9) observable with HIFI/Herschel, discriminate very well between XDRs and PDRs. Column density ratios indicate that CH, CH+, NO, HOC+ and HCO are good PDR/XDR discriminators.
[73]  [pdf] - 316513
Irradiated ISM: Discriminating between Cosmic Rays and X-rays
Comments: accepted by Astrophysical Journal Letters
Submitted: 2006-09-07
The ISM of active galaxy centers is exposed to a combination of cosmic ray, FUV and X-ray radiation. We apply PDR models to this ISM with both `normal' and highly elevated (5\times 10^{-15}s^-1) cosmic-ray rates and compare the results to those obtained for XDRs. Our existing PDR-XDR code is used to construct models over a 10^3-10^5 cm^-3 density range and for 0.16-160 erg s^-1 cm^-2 impingent fluxes. We obtain larger high J (J>10) CO ratios in PDRs when we use the highly elevated cosmic ray rate, but these are always exceeded by the corresponding XDR ratios. The [CI] 609 mum/13CO(2-1) line ratio is boosted by a factor of a few in PDRs with n~10^3 cm^-3 exposed to a high cosmic ray rate. At higher densities ratios become identical irrespective of cosmic ray flux, while XDRs always show elevated [CI] 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 cosmic ray irradiation conditions, and XDRs. Hence, the HIFI instrument on Herschel, which can detect these CO lines, will be crucial in the study of active galaxies.
[74]  [pdf] - 71245
Diagnostics of irradiated gas in galaxy nuclei. I: A Far-ultraviolet and X-ray dominated region code
Comments: 21 pages, accepted by A&A
Submitted: 2005-02-22
We present a far-ultraviolet (PDR) and an X-ray dominated region (XDR) code. We include and discuss thermal and chemical processes that pertain to irradiated gas. An elaborate chemical network is used and a careful treatment of PAHs and H2 formation, destruction and excitation is included. For both codes we calculate four depth-dependent models for different densities and radiation fields, relevant to conditions in starburst galaxies and active galactic nuclei. A detailed comparison between PDR and XDR physics is made for total gas column densities between ~10^20 and ~10^25 cm^-2. We show cumulative line intensities for a number of fine-structure lines (e.g., [CII], [OI], [CI], [SiII], [FeII]), as well as cumulative column densities and column density ratios for a number of species (e.g., CO/H2, CO/C, HCO+/HCN, HNC/HCN). The comparison between the results for the PDRs and XDRs shows that column density ratios are almost constant up to N_H=10^22 cm^-2 for XDRs, unlike those in PDRs. For example, CO/C in PDRs changes over four orders of magnitude from the edge to N_H=10^22 cm^-2. The CO/C and CO/H2 ratios are lower in XDRs at low column densities and rise at N_H > 10^23 cm^-2. At most column densities N_H > 10^21.5 cm^-2, HNC/HCN ratios are lower in XDRs too, but they show a more moderate increase at higher N_H.
[75]  [pdf] - 67592
A submillimeter exponential disk in M~51: evidence for an extended cold dust disk
Comments: accepted by A&A, 8 pages
Submitted: 2004-09-21
A 850 micron map of the interacting spiral galaxy M51 shows well-defined spiral arms, closely resembling the structures seen in CO and HI emission. However, most of the 850 micron emission originates in an underlying exponential disk, a component that has not been observed before in a face-on galaxy at these wavelengths. The scale-length of this disk is 5.45 kpc, which is somewhat larger than the scale-length of the stellar disk, but somewhat smaller than that of atomic hydrogen. Its profile can not be explained solely by a radial disk temperature gradient but requires the underlying dust to have an exponential distribution as well. This reinforces the view that the submm emission from spiral galaxy disks traces total hydrogen column density, i.e.the sum of H2 and HI. A canonical gas-to-dust ratio of 100+/-26 is obtained for kappa(850)=1.2 g**-1 cm**2, where kappa(850) is the dust opacity at 850 micron.
[76]  [pdf] - 363033
The post-AGB evolution of AGB mass loss variations
Comments: 12 pages, 13 figures, accepted by A&A
Submitted: 2003-05-20
We present new numerical hydrodynamical modelling of the evolution of Asymptotic Giant Branch (AGB) mass loss fluctuations during the post-AGB/Planetary Nebula phase. These models show that after ionization, the observable effects of the mass loss fluctuations disappear in a few thousand years, consistent with the fact that only few PNe have been found to be surrounded by `rings'. We derive the observational characteristics of these rings, and compare them to reported observations of the rings around NGC 6543, finding a good match of emission properties and line shapes. We predict small variations in the observable electron temperatures.