Photodissociation Region Models Research Articles

Constraining cloud parameters using high density gas tracers in galaxies

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.

Modelling mechanical heating in star-forming galaxies: CO and13CO Line ratios as sensitive probes

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}$.

RELATIONS WITH CO ROTATIONAL LADDERS OF GALAXIES ACROSS THE HERSCHEL SPIRE ARCHIVE

ABSTRACT We present a catalog of all CO (J = 4−3 through J = 13−12), [C i], and [N ii] lines available from extragalactic spectra from the Herschel SPIRE Fourier Transform Spectrometer (FTS) archive combined with observations of the low-J CO lines from the literature and from the Arizona Radio Observatory. This work examines the relationships between L FIR, , and L CO/L CO,1−0. We also present a new method for estimating probability distribution functions from marginal signal-to-noise ratio Herschel FTS spectra, which takes into account the instrumental “ringing” and the resulting highly correlated nature of the spectra. The slopes of log(L FIR) versus log( ) are linear for all mid- to high-J CO lines and slightly sublinear if restricted to (ultra)luminous infrared galaxies ((U)LIRGs). The mid- to high-J CO luminosity relative to CO J = 1−0 increases with increasing L FIR, indicating higher excitement of the molecular gas, although these ratios do not exceed ∼180. For a given bin in L FIR, the luminosities relative to CO J = 1−0 remain relatively flat from J = 6−5 through J = 13−12, across three orders of magnitude of L FIR. A single component theoretical photodissociation region (PDR) model cannot match these flat SLED shapes, although combinations of PDR models with mechanical heating added qualitatively match the shapes, indicating the need for further comprehensive modeling of the excitation processes of warm molecular gas in nearby galaxies.

MID-J CO SHOCK TRACING OBSERVATIONS OF INFRARED DARK CLOUDS. III. SLED FITTING

ABSTRACT Giant molecular clouds contain supersonic turbulence that can locally heat small fractions of gas to over 100 K. We run shock models for low-velocity, C-type shocks propagating into gas with densities between 103 and 105 cm−3 and find that CO lines are the most important cooling lines. Comparison to photodissociation region (PDR) models indicates that mid-J CO lines (J = 8 7 and higher) should be dominated by emission from shocked gas. In Papers I and II we presented CO J = 3 2, 8 7, and 9 8 observations toward four primarily quiescent clumps within infrared dark clouds. Here we fit PDR models to the combined spectral line energy distributions and show that the PDR models that best fit the low-J CO emission underpredict the mid-J CO emission by orders of magnitude, strongly hinting at a hot gas component within these clumps. The low-J CO data clearly show that the integrated intensities of both the CO J = 8 7 and 9 8 lines are anomalously high, such that the line ratio can be used to characterize the hot gas component. Shock models are reasonably consistent with the observed mid-J CO emission, with models with densities near cm−3 providing the best agreement. Where this mid-J CO is detected, the mean volume filling factor of the hot gas is 0.1%. Much of the observed mid-J CO emission, however, is also associated with known protostars and may be due to protostellar feedback.

PDR Emission from the Arched-Filaments and Nearby Positions

AbstractWe investigate the physical conditions of the gas, atomic and molecular, in the filaments in the context of Photo-Dissociation Regions (PDRs) using the KOSMA-PDR mode of clumpy clouds. We also compare the [CII] vs. [NII] integrated intensity predictions in Abel et al. 2005 for HII regions and adjacent PDRs in the Galactic disk, and check for their applicability under the extreme physical conditions present in the GC. Our preliminary results show that observed integrated intensities are well reproduced by the PDR model. The gas is exposed to a relatively low Far-UV field between 102 – 103 Draine fields. The total volume hydrogen density is well constrained between 104 – 105 cm−3. The hydrogen ionization rate due to cosmic-rays varies between 10−15 and 4× 10−15 s−1, with the highest value ~ 10−14 s−1 found towards G0.07+0.04. Our results show that the line-of-sight contribution to the total distance of the filaments to the Arches Cluster is not negligible. The spatial distribution of the [CII]/[NII] ratio shows that the integrated intensity ratios are fairly homogeneously distributed for values below 10 in energy units. Calculations including variation on the [C/N] abundance ratio show that tight constraints on this ratio are needed to reproduce the observations.

A self-consistent model for the evolution of the gas produced in the debris disc of β Pictoris

This paper presents a self-consistent model for the evolution of gas produced in the debris disc of β Pictoris. Our model proposes that atomic carbon and oxygen are created from the photodissociation of CO, which is itself released from volatile-rich bodies in the debris disc due to grain–grain collisions or photodesorption. While the CO lasts less than one orbit, the atomic gas evolves by viscous spreading resulting in an accretion disc inside the parent belt and a decretion disc outside. The temperature, ionization fraction and population levels of carbon and oxygen are followed with the photodissociation region model cloudy, which is coupled to a dynamical viscous α model. We present new gas observations of β Pic, of C i observed with Atacama Pathfinder EXperiment and O i observed with Herschel, and show that these along with published C ii and CO observations can all be explained with this new model. Our model requires a viscosity α > 0.1, similar to that found in sufficiently ionized discs of other astronomical objects; we propose that the magnetorotational instability is at play in this highly ionized and dilute medium. This new model can be tested from its predictions for high-resolution ALMA observations of C i. We also constrain the water content of the planetesimals in β Pic. The scenario proposed here might be at play in all debris discs and this model could be used more generally on all discs with C, O or CO detections.

Efficient ortho-para conversion of H2on interstellar grain surfaces

Context: Fast surface conversion between ortho- and para-H2 has been observed in laboratory studies, and this mechanism has been proposed to play a role in the control of the ortho-para ratio in the interstellar medium. Observations of rotational lines of H2 in Photo-Dissociation Regions (PDRs) have indeed found significantly lower ortho-para ratios than expected at equilibrium. The mechanisms controlling the balance of the ortho-para ratio in the interstellar medium thus remain incompletely understood, while this ratio can affect the thermodynamical properties of the gas (equation of state, cooling function). Aims: We aim to build an accurate model of ortho-para conversion on dust surfaces based on the most recent experimental and theoretical results, and to validate it by comparison to observations of H2 rotational lines in PDRs. Methods: We propose a statistical model of ortho-para conversion on dust grains with fluctuating dust temperatures, based on a master equation approach. This computation is then coupled to full PDR models and compared to PDR observations. Results: We show that the observations of rotational H2 lines indicate a high conversion efficiency on dust grains, and that this high efficiency can be accounted for if taking dust temperature fluctuations into account with our statistical model of surface conversion. Simpler models neglecting the dust temperature fluctuations do not reach the high efficiency deduced from the observations. Moreover, this high efficiency induced by dust temperature fluctuations is quite insensitive to the values of microphysical parameters of the model. Conclusions: Ortho-para conversion on grains is thus an efficient mechanism in most astrophysical conditions that can play a significant role in controlling the ortho-para ratio.

HerschelPACS and SPIRE spectroscopy of the photodissociation regions associated with S 106 and IRAS 23133+6050

Photodissociation regions (PDRs) contain a large fraction of all of the interstellar matter in galaxies. Classical examples include the boundaries between ionized regions and molecular clouds in regions of massive star formation, marking the point where all of the photons energetic enough to ionize hydrogen have been absorbed. In this paper we determine the physical properties of the PDRs associated with the star forming regions IRAS 23133+6050 and S 106 and present them in the context of other Galactic PDRs associated with massive star forming regions. We employ Herschel PACS and SPIRE spectroscopic observations to construct a full 55-650 {\mu}m spectrum of each object from which we measure the PDR cooling lines, other fine- structure lines, CO lines and the total far-infrared flux. These measurements are then compared to standard PDR models. Subsequently detailed numerical PDR models are compared to these predictions, yielding additional insights into the dominant thermal processes in the PDRs and their structures. We find that the PDRs of each object are very similar, and can be characterized by a two-phase PDR model with a very dense, highly UV irradiated phase (n $\sim$ 10^6 cm^(-3), G$_0$ $\sim$ 10^5) interspersed within a lower density, weaker radiation field phase (n $\sim$ 10^4 cm^(-3), G$_0$ $\sim$ 10^4). We employed two different numerical models to investigate the data, firstly we used RADEX models to fit the peak of the $^{12}$CO ladder, which in conjunction with the properties derived yielded a temperature of around 300 K. Subsequent numerical modeling with a full PDR model revealed that the dense phase has a filling factor of around 0.6 in both objects. The shape of the $^{12}$CO ladder was consistent with these components with heating dominated by grain photoelectric heating. An extra excitation component for the highest J lines (J > 20) is required for S 106.

Mid-JCO shock tracing observations of infrared dark clouds. I.

Infrared dark clouds (IRDCs) are dense, molecular structures in the interstellar medium that can harbour sites of high-mass star formation. IRDCs contain supersonic turbulence, which is expected to generate shocks that locally heat pockets of gas within the clouds. We present observations of the CO J = 8-7, 9-8, and 10-9 transitions, taken with the Herschel Space Observatory, towards four dense, starless clumps within IRDCs (C1 in G028.37+00.07, F1 and F2 in G034.43+0007, and G2 in G034.77-0.55). We detect the CO J = 8-7 and 9-8 transitions towards three of the clumps (C1, F1, and F2) at intensity levels greater than expected from photodissociation region (PDR) models. The average ratio of the 8-7 to 9-8 lines is also found to be between 1.6 and 2.6 in the three clumps with detections, significantly smaller than expected from PDR models. These low line ratios and large line intensities strongly suggest that the C1, F1, and F2 clumps contain a hot gas component not accounted for by standard PDR models. Such a hot gas component could be generated by turbulence dissipating in low velocity shocks.

Erratum: Mid-J CO observations of Perseus B1-East 5: evidence for turbulent dissipation via low-velocity shocks

Giant molecular clouds contain supersonic turbulence and magnetohydrodynamic simulations predict that this turbulence should decay rapidly. Such turbulent dissipation has the potential to create a warm (T ~100 K) gas component within a molecular cloud. We present observations of the CO J = 5-4 and 6-5 transitions, taken with the Herschel Space Observatory, towards the Perseus B1-East 5 region. We combine these new observations with archival measurements of lower rotational transitions and fit photodissociation region models to the data. We show that Perseus B1-E5 has an anomalously large CO J = 6-5 integrated intensity, consistent with a warm gas component existing within the region. This excess emission is consistent with predictions for shock heating due to the dissipation of turbulence in low velocity shocks with the shocks having a volume filling factor of 0.15 per cent. We find that B1-E has a turbulent energy dissipation rate of 3.5 x 10$^{32}$ erg / s and a dissipation time-scale that is only a factor of 3 larger than the flow crossing time-scale.

Mid-J CO observations of Perseus B1-East 5: evidence for turbulent dissipation via low-velocity shocks

Giant molecular clouds contain supersonic turbulence and magnetohydrodynamic simulations predict that this turbulence should decay rapidly. Such turbulent dissipation has the potential to create a warm (T ~100 K) gas component within a molecular cloud. We present observations of the CO J = 5-4 and 6-5 transitions, taken with the Herschel Space Observatory, towards the Perseus B1-East 5 region. We combine these new observations with archival measurements of lower rotational transitions and fit photodissociation region models to the data. We show that Perseus B1-E5 has an anomalously large CO J = 6-5 integrated intensity, consistent with a warm gas component existing within the region. This excess emission is consistent with predictions for shock heating due to the dissipation of turbulence in low velocity shocks with the shocks having a volume filling factor of 0.15 per cent. We find that B1-E has a turbulent energy dissipation rate of 3.5 x 10$^{32}$ erg / s and a dissipation time-scale that is only a factor of 3 larger than the flow crossing time-scale.

MID-INFRARED PROPERTIES OF LUMINOUS INFRARED GALAXIES. II. PROBING THE DUST AND GAS PHYSICS OF THE GOALS SAMPLE

The Great Observatories All-sky LIRG Survey (GOALS) is a comprehensive, multiwavelength study of luminous infrared galaxies (LIRGs) in the local universe. Here, we present the results of a multi-component, spectral decomposition analysis of the low-resolution mid-infrared (MIR) Spitzer Infrared Spectrograph spectra from 5–38 μm of 244 LIRG nuclei. The detailed fits and high-quality spectra allow for characterization of the individual polycyclic aromatic hydrocarbon (PAH) features, warm molecular hydrogen emission, and optical depths for both silicate dust grains and water ices. We find that starbursting LIRGs, which make up the majority of the GOALS sample, are very consistent in their MIR properties (i.e., τ9.7 μm, τice, neon line ratios, and PAH feature ratios). However, as their EQW6.2 μm decreases, usually an indicator of an increasingly dominant active galactic nucleus (AGN), LIRGs cover a larger spread in these MIR parameters. The contribution from PAH emission to the total IR luminosity (L(PAH)/L(IR)) in LIRGs varies from 2%–29% and LIRGs prior to their first encounter show significantly higher L(PAH)/L(IR) ratios on average. We observe a correlation between the strength of the starburst (represented by IR8 = LIR/L8 μm) and the PAH fraction at 8 μm but no obvious link between IR8 and the 7.7 to 11.3 PAH ratio, suggesting that the fractional photodissociation region (PDR) emission, and not the overall grain properties, is associated with the rise in IR8 for galaxies off the starburst main sequence. We detect crystalline silicate features in ∼6% of the sample but only in the most obscure sources (s9.7 μm < −1.24). Ice absorption features are observed in ∼11% (56%) of GOALS LIRGs (ULIRGs) in sources with a range of silicate depths. Most GOALS LIRGs have L(H2)/L(PAH) ratios elevated above those observed for normal star-forming galaxies and exhibit a trend for increasing L(H2)/L(PAH) ratio with increasing L(H2). While star formation appears to be the dominant process responsible for exciting the H2 in most of the GOALS galaxies, a subset of LIRGs (∼10%) shows excess H2 emission that is inconsistent with PDR models and may be excited by shocks or AGN-induced outflows.

Gas properties in the disc of NGC 891 from Herschel far-infrared spectroscopy

AbstractWe investigate the physical properties of the interstellar gas in the nearby edge-on spiral galaxy NGC 891, using Herschel PACS/SPIRE observations of the most important far-infrared (FIR) cooling lines – [Cii] 158 μm, [Nii] 122, 205 μm, [Oi] 63, 145 μm, and [Oiii] 88 μm – obtained as part of the Very Nearby Galaxy Survey (P. I.: C. D. Wilson). We compare our observations to the predictions of a photo dissociation region (PDR) model to determine the gas density, n, and the strength of the incident FUV radiation field, G0, on a pixel-by-pixel basis. The majority of PDRs in NGC 891's disc exhibit properties similar to the physical conditions found in the spiral arm and inter-arm regions of the face-on M51 galaxy. We estimate a stronger FUV field in the far north-eastern side than compared to the rest of the disc.

Warm molecular gas temperature distribution in six local infrared bright Seyfert galaxies

We simultaneously analyze the spectral line energy distributions (SLEDs) of CO and H2 of six local luminous infrared (IR) Seyfert galaxies. For the CO SLEDs, we used new Herschel/SPIRE FTS data (from J=4-3 to J=13-12) and ground-based observations for the lower-J CO transitions. The H2 SLEDs were constructed using archival mid-IR Spitzer/IRS and near-IR VLT/SINFONI data for the rotational and ro-vibrational H2 transitions, respectively. In total, the SLEDs contain 26 transitions with upper level energies between 5 and 15000 K. A single, constant density, model (n$_{H_2}$ ~ 10$^{4.5-6}$ cm$^{-3}$) with a broken power-law temperature distribution reproduces well both the CO and H2 SLEDs. The power-law indices are $\beta_1$ ~ 1-3 for warm molecular gas (20 K < T < 100 K) and $\beta_2$ ~ 4-5 for hot molecular gas (T > 100 K). We show that the steeper temperature distribution (higher $\beta$) for hot molecular gas can be explained by shocks and photodissociation region (PDR) models, however, the exact $\beta$ values are not reproduced by PDR or shock models alone and a combination of both is needed. We find that the three major mergers among our targets have shallower temperature distributions for warm molecular gas than the other three spiral galaxies. This can be explained by a higher relative contribution of shock excitation, with respect to PDR excitation, for the warm molecular gas in these mergers. For only one of the mergers, IRASF 05189-2524, the shallower H2 temperature distribution differs from that of the spiral galaxies. The presence of a bright active galactic nucleus in this source might explain the warmer molecular gas observed.

EXTENDED [C II] EMISSION IN LOCAL LUMINOUS INFRARED GALAXIES

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.

THE CO-TO-H2CONVERSION FACTOR ACROSS THE PERSEUS MOLECULAR CLOUD

We derive the CO-to-H2 conversion factor, X_CO = N(H2)/I_CO, across the Perseus molecular cloud on sub-parsec scales by combining the dust-based N(H2) data with the I_CO data from the COMPLETE Survey. We estimate an average X_CO ~ 3 x 10^19 cm^-2 K^-1 km^-1 s and find a factor of ~3 variations in X_CO between the five sub-regions in Perseus. Within the individual regions, X_CO varies by a factor of ~100, suggesting that X_CO strongly depends on local conditions in the interstellar medium. We find that X_CO sharply decreases at Av < 3 mag but gradually increases at Av > 3 mag, with the transition occurring at Av where I_CO becomes optically thick. We compare the N(HI), N(H2), I_CO, and X_CO distributions with two models of the formation of molecular gas, a one-dimensional photodissociation region (PDR) model and a three-dimensional magnetohydrodynamic (MHD) model tracking both the dynamical and chemical evolution of gas. The PDR model based on the steady state and equilibrium chemistry reproduces our data very well but requires a diffuse halo to match the observed N(HI) and I_CO distributions. The MHD model generally matches our data well, suggesting that time-dependent effects on H2 and CO formation are insignificant for an evolved molecular cloud like Perseus. However, we find interesting discrepancies, including a broader range of N(HI), likely underestimated I_CO, and a large scatter of I_CO at small Av. These discrepancies likely result from strong compressions/rarefactions and density fluctuations in the MHD model.

Using [C i] to probe the interstellar medium in z ∼ 2.5 sub-millimeter galaxies★

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.

Herschelobservations of the Sagittarius B2 cores: Hydrides, warm CO, and cold dust

Sagittarius B2 (Sgr B2) is one of the most massive and luminous star-forming regions in the Galaxy and shows chemical and physical conditions similar to those in distant extragalactic starbursts. We present large-scale far-IR/submm photometric images and spectroscopic maps taken with the PACS and SPIRE instruments onboard Herschel. The spectra towards the Sgr B2 star-forming cores, B2(M) and B2(N), are characterized by strong CO line emission, emission lines from high-density tracers (HCN, HCO+, and H2S), [N II] 205 um emission from ionized gas, and absorption lines from hydride molecules (OH+, H2O+, H2O, CH+, CH, SH+, HF, NH, NH2, and NH3). The rotational population diagrams of CO suggest the presence of two gas temperature components: an extended warm component, which is associated with the extended envelope, and a hotter component, which is seen towards the B2(M) and B2(N) cores. As observed in other Galactic Center clouds, the gas temperatures are significantly higher than the dust temperatures inferred from photometric images. We determined far-IR and total dust masses in the cores. Non-local thermodynamic equilibrium models of the CO excitation were used to constrain the averaged gas density in the cores. A uniform luminosity ratio is measured along the extended envelope, suggesting that the same mechanism dominates the heating of the molecular gas at large scales. The detection of high-density molecular tracers and of strong [N II] 205 um line emission towards the cores suggests that their morphology must be clumpy to allow UV radiation to escape from the inner HII regions. Together with shocks, the strong UV radiation field is likely responsible for the heating of the hot CO component. At larger scales, photodissociation regions models can explain both the observed CO line ratios and the uniform L(CO)/LFIR luminosity ratios.

Carbon fractionation in photo-dissociation regions

We upgraded the chemical network from the UMIST Database for Astrochemistry 2006 to include isotopes such as ^{13}C and ^{18}O. This includes all corresponding isotopologues, their chemical reactions and the properly scaled reaction rate coefficients. We study the fractionation behavior of astrochemically relevant species over a wide range of model parameters, relevant for modelling of photo-dissociation regions (PDRs). We separately analyze the fractionation of the local abundances, fractionation of the total column densities, and fractionation visible in the emission line ratios. We find that strong C^+ fractionation is possible in cool C^+ gas. Optical thickness as well as excitation effects produce intensity ratios between 40 and 400. The fractionation of CO in PDRs is significantly different from the diffuse interstellar medium. PDR model results never show a fractionation ratio of the CO column density larger than the elemental ratio. Isotope-selective photo-dissociation is always dominated by the isotope-selective chemistry in dense PDR gas. The fractionation of C, CH, CH^+, and HCO^+ is studied in detail, showing that the fractionation of C, CH and CH^+ is dominated by the fractionation of their parental species. The light hydrides chemically derive from C^+, and, consequently, their fractionation state is coupled to that of C^+. The fractionation of C is a mixed case depending on whether formation from CO or HCO^+ dominates. Ratios of the emission lines of [C II], [C I], ^{13}CO, and H^{13}CO^+ provide individual diagnostics to the fractionation status of C^+, C, and CO.

A STUDY OF HEATING AND COOLING OF THE ISM IN NGC 1097 WITHHERSCHEL-PACS ANDSPITZER-IRS

NGC 1097 is a nearby Seyfert 1 galaxy with a bright circumnuclear starburst ring, a strong large-scale bar and an active nucleus. We present a detailed study of the spatial variation of the far infrared (FIR) [CII]158um and [OI]63um lines and mid-infrared H2 emission lines as tracers of gas cooling, and of the polycyclic aromatic hydrocarbon (PAH) bands as tracers of the photoelectric heating, using Herschel-PACS, and Spitzer-IRS infrared spectral maps. We focus on the nucleus and the ring, and two star forming regions (Enuc N and Enuc S). We estimated a photoelectric gas heating efficiency ([CII]158um+[OI]63um)/PAH in the ring about 50% lower than in Enuc N and S. The average 11.3/7.7um PAH ratio is also lower in the ring, which may suggest a larger fraction of ionized PAHs, but no clear correlation with [CII]158{\mu}m/PAH(5.5 - 14um) is found. PAHs in the ring are responsible for a factor of two more [CII]158um and [OI]63um emission per unit mass than PAHs in the Enuc S. SED modeling indicates that at most 25% of the FIR power in the ring and Enuc S can come from high intensity photodissociation regions (PDRs), in which case G0 ~ 10^2.3 and nH ~ 10^3.5 cm^-3 in the ring. For these values of G0 and nH PDR models cannot reproduce the observed H2 emission. Much of the the H2 emission in the starburst ring could come from warm regions in the diffuse ISM that are heated by turbulent dissipation or shocks.