Source author record

V. Ossenkopf

V. Ossenkopf appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

49works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

49 published item(s)

preprint2016arXiv

Understanding star formation in molecular clouds III. Probability distribution functions of molecular lines in Cygnus X

Column density (N) PDFs serve as a powerful tool to characterize the physical processes that influence the structure of molecular clouds. Star-forming clouds can best be characterized by lognormal PDFs for the lower N range and a power-law tail for higher N, commonly attributed to turbulence and self-gravity and/or pressure, respectively. We report here on PDFs obtained from observations of 12CO, 13CO, C18O, CS, and N2H+ in the Cygnus X North region and compare to a PDF derived from dust observations with the Herschel satellite. The PDF of 12CO is lognormal for Av~1-30, but is cut for higher Av due to optical depth effects. The PDFs of C18O and 13CO are mostly lognormal up for Av~1-15, followed by excess up to Av~40. Above that value, all CO PDFs drop, most likely due to depletion. The high density tracers CS and N2H+ exhibit only a power law distribution between Av~15 and 400, respectively. The PDF from dust is lognormal for Av~3-15 and has a power-law tail up to Av~500. Absolute values for the molecular line column densities are, however, rather uncertain due to abundance and excitation temperature variations. Taken the dust PDF face value, we 'calibrate' the molecular line PDF of CS to the one of the dust and determined an abundance [CS]/[H2] of 10^-9. The slopes of the power-law tails of the CS, N2H+, and dust PDFs are consistent with free-fall collapse of filaments and clumps. A quasi static configuration of filaments and clumps can possibly also account for the observed N-PDFs, as long as they have a sufficiently condensed density structure and external ram pressure by gas accretion is provided. The somehow flatter slopes of N2H+ and CS can reflect an abundance change and/or subthermal excitation at low column densities.

preprint2015arXiv

C$_2$H observations toward the Orion Bar

C$_2$H is one of the first radicals to be detected in the interstellar medium. Its higher rotational transitions have recently become available with the Herschel Space Observatory. We aim to constrain the physical parameters of the C$_2$H emitting gas toward the Orion Bar. We analyse the C$_2$H line intensities measured toward the Orion Bar CO$^+$ Peak and Herschel/HIFI maps of C$_2$H, CH, and HCO$^+$, and a NANTEN map of [CI]. We interpret the observed C$_2$H emission using radiative transfer and PDR models. Five rotational transitions of C$_2$H have been detected in the HIFI frequency range toward the CO$^+$ peak. A single component rotational diagram gives a rotation temperature of ~64 K and a beam-averaged C$_2$H column density of 4$\times$10$^{13}$ cm$^{-2}$. The measured transitions cannot be explained by any single parameter model. According to a non-LTE model, most of the C$_2$H column density produces the lower-$N$ C$_2$H transitions and traces a warm ($T_{\rm{kin}}$ ~ 100-150 K) and dense ($n$(H$_2$)~10$^5$-10$^6$ cm$^{-3}$) gas. A small fraction of the C$_2$H column density is required to reproduce the intensity of the highest-$N$ transitions ($N$=9-8 and N=10-9) originating from a high density ($n$(H$_2$)~5$\times$10$^6$ cm$^{-3}$) hot ($T_{\rm{kin}}$ ~ 400 K) gas. The total beam-averaged C$_2$H column density in the model is 10$^{14}$ cm$^{-2}$. Both the non-LTE radiative transfer model and a simple PDR model representing the Orion Bar with a plane-parallel slab of gas and dust suggest, that C$_2$H cannot be described by a single pressure component, unlike the reactive ion CH$^+$, which was previously analysed toward the Orion Bar CO$^+$ peak. The physical parameters traced by the higher rotational transitions ($N$=6-5,...,10-9) of C$_2$H may be consistent with the edges of dense clumps exposed to UV radiation near the ionization front of the Orion Bar.

preprint2015arXiv

Detection of a large fraction of atomic gas not associated with star-forming material in M17 SW

We probe the column densities and masses traced by the ionized and neutral atomic carbon with spectrally resolved maps, and compare them to the diffuse and dense molecular gas traced by [C I] and low-$J$ CO lines toward the star-forming region M17SW. We mapped a 4.1pc x 4.7pc region in the [C I] 609 m$μ$ line using the APEX telescope, as well as the CO isotopologues with the IRAM 30m telescope. We analyze the data based on velocity channel maps that are 1 km/s wide. We correlate their spatial distribution with that of the [C II] map obtained with SOFIA/GREAT. Optically thin approximations were used to estimate the column densities of [C I] and [C II] in each velocity channel. The spatial distribution of the [C I] and all CO isotopologues emission was found to be associated with that of [C II] in about 20%-80% of the mapped region, with the high correlation found in the central (15-23 km/s ) velocity channels. The excitation temperature of [C I] ranges between 40 K and 100 K in the inner molecular region of M17 SW. Column densities in 1 km/s channels between ~10$^{15}$ and ~10$^{17}$ cm$^{-2}$ were found for [C I]. Just ~20% of the velocity range (~40 km/s) that the [C II] line spans is associated with the star-forming material traced by [C I] and CO. The total gas mass estimated from the [C II] emission gives a lower limit of ~4.4x10$^3$ $M_{\odot}$. At least 64% of this mass is not associated with the star-forming material in M17SW. We also found that about 36%, 17%, and 47% of the [C II] emission is associated with the HII, HI, and H_2 regimes, respectively. Comparisons with the H41$α$ line shows an ionization region mixed with the neutral and part of the molecular gas, in agreement with the clumped structure and dynamical processes at play in M17SW. These results are also relevant to extra-galactic studies in which [C II] is often used as a tracer of star-forming material.

preprint2015arXiv

Detection of two power-law tails in the probability distribution functions of massive GMCs

We report the novel detection of complex high-column density tails in the probability distribution functions (PDFs) for three high-mass star-forming regions (CepOB3, MonR2, NGC6334), obtained from dust emission observed with Herschel. The low column density range can be fit with a lognormal distribution. A first power-law tail starts above an extinction (Av) of ~6-14. It has a slope of alpha=1.3-2 for the rho~r^-alpha profile for an equivalent density distribution (spherical or cylindrical geometry), and is thus consistent with free-fall gravitational collapse. Above Av~40, 60, and 140, we detect an excess that can be fitted by a flatter power law tail with alpha>2. It correlates with the central regions of the cloud (ridges/hubs) of size ~1 pc and densities above 10^4 cm^-3. This excess may be caused by physical processes that slow down collapse and reduce the flow of mass towards higher densities. Possible are: 1. rotation, which introduces an angular momentum barrier, 2. increasing optical depth and weaker cooling, 3. magnetic fields, 4. geometrical effects, and 5. protostellar feedback. The excess/second power-law tail is closely linked to high-mass star-formation though it does not imply a universal column density threshold for the formation of (high-mass) stars.

preprint2015arXiv

Disentangling the excitation conditions of the dense gas in M17 SW

We probe the chemical and energetic conditions in dense gas created by radiative feedback through observations of multiple CO, HCN and HCO$^+$ transitions toward the dense core of M17 SW. We used the dual band receiver GREAT on board the SOFIA airborne telescope to obtain maps of the $J=16-15$, $J=12-11$, and $J=11-10$ transitions of $^{12}$CO. We compare these maps with corresponding APEX and IRAM 30m telescope data for low- and mid-$J$ CO, HCN and HCO$^+$ emission lines, including maps of the HCN $J=8-7$ and HCO$^+$ $J=9-8$ transitions. The excitation conditions of $^{12}$CO, HCO$^+$ and HCN are estimated with a two-phase non-LTE radiative transfer model of the line spectral energy distributions (LSEDs) at four selected positions. The energy balance at these positions is also studied. We obtained extensive LSEDs for the CO, HCN and HCO$^+$ molecules toward M17 SW. The LSED shape, particularly the high-$J$ tail of the CO lines observed with SOFIA/GREAT, is distinctive for the underlying excitation conditions. The critical magnetic field criterion implies that the cold cloudlets at two positions are partially controlled by processes that create and dissipate internal motions. Supersonic but sub-Alfvénic velocities in the cold component at most selected positions indicates that internal motions are likely MHD waves. Magnetic pressure dominates thermal pressure in both gas components at all selected positions, assuming random orientation of the magnetic field. The magnetic pressure of a constant magnetic field throughout all the gas phases can support the total internal pressure of the cold components, but it cannot support the internal pressure of the warm components. If the magnetic field scales as $B \propto n^{2/3}$, then the evolution of the cold cloudlets at two selected positions, and the warm cloudlets at all selected positions, will be determined by ambipolar diffusion.

preprint2015arXiv

Temperatures of dust and gas in S~140

In dense parts of interstellar clouds (> 10^5 cm^-3), dust & gas are expected to be in thermal equilibrium, being coupled via collisions. However, previous studies have shown that the temperatures of the dust & gas may remain decoupled even at higher densities. We study in detail the temperatures of dust & gas in the photon-dominated region S 140, especially around the deeply embedded infrared sources IRS 1-3 and at the ionization front. We derive the dust temperature and column density by combining Herschel PACS continuum observations with SOFIA observations at 37 $μ$m and SCUBA at 450 $μ$m. We model these observations using greybody fits and the DUSTY radiative transfer code. For the gas part we use RADEX to model the CO 1-0, CO 2-1, 13CO 1-0 and C18O 1-0 emission lines mapped with the IRAM-30m over a 4' field. Around IRS 1-3, we use HIFI observations of single-points and cuts in CO 9-8, 13CO 10-9 and C18O 9-8 to constrain the amount of warm gas, using the best fitting dust model derived with DUSTY as input to the non-local radiative transfer model RATRAN. We find that the gas temperature around the infrared sources varies between 35 and 55K and that the gas is systematically warmer than the dust by ~5-15K despite the high gas density. In addition we observe an increase of the gas temperature from 30-35K in the surrounding up to 40-45K towards the ionization front, most likely due to the UV radiation from the external star. Furthermore, detailed models of the temperature structure close to IRS 1 show that the gas is warmer and/or denser than what we model. Finally, modelling of the dust emission from the sub-mm peak SMM 1 constrains its luminosity to a few ~10^2 Lo. We conclude that the gas heating in the S 140 region is very efficient even at high densities, most likely due to the deep UV penetration from the embedded sources in a clumpy medium and/or oblique shocks.

preprint2015arXiv

Understanding star formation in molecular clouds I. Effects of line-of-sight contamination on the column density structure

Column-density maps of molecular clouds are one of the most important observables in the context of molecular cloud- and star-formation (SF) studies. With the Herschel satellite it is now possible to determine the column density from dust emission. We use observations and simulations to demonstrate how LOS contamination affects the column density probability distribution function (PDF). We apply a first-order approximation (removing a constant level) to the molecular clouds of Auriga, Maddalena, Carina and NGC3603. In perfect agreement with the simulations, we find that the PDFs become broader, the peak shifts to lower column densities, and the power-law tail of the PDF flattens after correction. All PDFs have a lognormal part for low column densities with a peak at Av~2, a deviation point (DP) from the lognormal at Av(DP)~4-5, and a power-law tail for higher column densities. Assuming a density distribution rho~r^-alpha, the slopes of the power-law tails correspond to alpha(PDF)=1.8, 1.75, and 2.5 for Auriga, Carina, and NGC3603 (alpha~1.5-2 is consistent gravitational collapse). We find that low-mass and high-mass SF clouds display differences in the overall column density structure. Massive clouds assemble more gas in smaller cloud volumes than low-mass SF ones. However, for both cloud types, the transition of the PDF from lognormal shape into power-law tail is found at the same column density (at Av~4-5 mag). Low-mass and high-mass SF clouds then have the same low column density distribution, most likely dominated by supersonic turbulence. At higher column densities, collapse and external pressure can form the power-law tail. The relative importance of the two processes can vary between clouds and thus lead to the observed differences in PDF and column density structure.

preprint2015arXiv

Understanding star formation in molecular clouds II. Signatures of gravitational collapse of IRDCs

We analyse column density and temperature maps derived from Herschel dust continuum observations of a sample of massive infrared dark clouds (G11.11-0.12, G18.82-0.28, G28.37+0.07, G28.53-0.25). We disentangle the velocity structure of the clouds using 13CO 1-0 and 12CO 3-2 data, showing that these IRDCs are the densest regions in massive giant molecular clouds and not isolated features. The probability distribution function (PDF) of column densities for all clouds have a power-law distribution over all (high) column densities, regardless of the evolutionary stage of the cloud: G11.11-0.12, G18.82-0.28, and G28.37+0.07 contain (proto)-stars, while G28.53-0.25 shows no signs of star formation. This is in contrast to the purely log-normal PDFs reported for near/mid-IR extinction maps. We only find a log-normal distribution for lower column densities, if we perform PDFs of the column density maps of the whole GMC in which the IRDCs are embedded. By comparing the PDF slope and the radial column density profile, we attribute the power law to the effect of large-scale gravitational collapse and to local free-fall collapse of pre- and protostellar cores. Independent from the PDF analysis, we find infall signatures in the spectral profiles of 12CO for G28.37+0.07 and G11.11-0.12, supporting the scenario of gravitational collapse. IRDCs are the densest regions within GMCs, which may be the progenitors of massive stars or clusters. At least some of the IRDCs are probably the same features as ridges (high column density regions with N>1e23 cm-2 over small areas), which were defined for nearby IR-bright GMCs. Because IRDCs are only confined to the densest (gravity dominated) cloud regions, the PDF constructed from this kind of a clipped image does not represent the (turbulence dominated) low column density regime of the cloud.

preprint2015arXiv

Velocity resolved [CII] spectroscopy of the center and the BCLMP302 region of M33 (HerM33es)

We aim to understand the contribution of the ionized, atomic and molecular phases of the ISM to the [CII] emission from clouds near the dynamical center and the BCLMP302 HII region in the north of the nearby galaxy M33 at a spatial resolution of 50pc. We combine high resolution [CII] spectra taken with the HIFI spectrometer onboard the Herschel satellite with [CII] Herschel-PACS maps and ground-based observations of CO(2-1) and HI. All data are at a common spatial resolution of 50pc. Typically, the [CII] lines have widths intermediate between the narrower CO(2-1) and broader HI line profiles. We decomposed the [CII] spectra in terms of contribution from molecular and atomic gas detected in CO(2-1) and HI, respectively. We find that the relative contribution of molecular and atomic gas traced by CO(2-1) and HI varies depends mostly on the local physical conditions and geometry. We estimate that 11-60% and 5-34% of the [CII] intensities in the center and in BCLMP302, respectively, arise at velocities showing no CO(2-1) or HI emission and could arise in CO-dark molecular gas. The deduced strong variation in the [CII] emission not associated with CO and HI cannot be explained in terms of differences in A_v, far-ultraviolet radiation field, and metallicity between the two studied regions. Hence the relative amounts of diffuse (CO-dark) and dense molecular gas possibly vary on spatial scales smaller than 50pc. Based on the emission measure observed at radio wavelengths we estimate the contribution of ionized gas at a few positions to lie between 10-25%. The correlations between the intensities of tracers corresponding to the same velocity range as [CII], differ from the correlation derived from PACS data. The results in this paper emphasize the need for velocity-resolved observations to discern the contribution of different components of the ISM to [CII] emission. (abridged)

preprint2015arXiv

Velocity-resolved [CII] emission and [CII]/FIR Mapping along Orion with Herschel

We present the first 7.5'x11.5' velocity-resolved map of the [CII]158um line toward the Orion molecular cloud-1 (OMC-1) taken with the Herschel/HIFI instrument. In combination with far-infrared (FIR) photometric images and velocity-resolved maps of the H41alpha hydrogen recombination and CO J=2-1 lines, this data set provides an unprecedented view of the intricate small-scale kinematics of the ionized/PDR/molecular gas interfaces and of the radiative feedback from massive stars. The main contribution to the [CII] luminosity (~85%) is from the extended, FUV-illuminated face of the cloud G_0>500, n_H>5x10^3 cm^-3) and from dense PDRs (G_0~10^4, n_H~10^5 cm^-3) at the interface between OMC-1 and the HII region surrounding the Trapezium cluster. Around 15% of the [CII] emission arises from a different gas component without CO counterpart. The [CII] excitation, PDR gas turbulence, line opacity (from [13CII]) and role of the geometry of the illuminating stars with respect to the cloud are investigated. We construct maps of the [CII]/FIR and FIR/M_Gas ratios and show that [CII]/FIR decreases from the extended cloud component (10^-2-10^-3) to the more opaque star-forming cores (10^-3-10^-4). The lowest values are reminiscent of the "[CII] deficit" seen in local ultra-luminous IR galaxies hosting vigorous star formation. Spatial correlation analysis shows that the decreasing [CII]/FIR ratio correlates better with the column density of dust through the molecular cloud than with FIR/M_Gas. We conclude that the [CII] emitting column relative to the total dust column along each line of sight is responsible for the observed [CII]/FIR variations through the cloud.

preprint2015arXiv

Wavelet-based cross-correlation analysis of structure scaling in turbulent clouds

We propose a statistical tool to compare the scaling behaviour of turbulence in pairs of molecular cloud maps. Using artificial maps with well defined spatial properties, we calibrate the method and test its limitations to ultimately apply it to a set of observed maps. We develop the wavelet-based weighted cross-correlation (WWCC) method to study the relative contribution of structures of different sizes and their degree of correlation in two maps as a function of spatial scale, and the mutual displacement of structures in the molecular cloud maps. We test the WWCC for circular structures having a single prominent scale and fractal structures showing a self-similar behavior without prominent scales. Observational noise and a finite map size limit the scales where the cross-correlation coefficients and displacement vectors can be reliably measured. For fractal maps containing many structures on all scales, the limitation from the observational noise is negligible for signal-to-noise ratios >5. (abbrev). Application of the WWCC to the observed line maps of the giant molecular cloud G333 allows to add specific scale information to the results obtained earlier using the principle component analysis. It confirms the chemical and excitation similarity of $^{13}$CO and C$^{18}$O on all scales, but shows a deviation of HCN at scales of up to 7' (~7 pc). This can be interpreted as a chemical transition scale. The largest structures also show a systematic offset along the filament, probably due to a large-scale density gradient. The WWCC can compare correlated structures in different maps of molecular clouds identifying scales that represent structural changes such as chemical and phase transitions and prominent or enhanced dimensions.

preprint2014arXiv

Carbon in different phases ([CII], [CI], and CO) in infrared dark clouds: Cloud formation signatures and carbon gas fractions

Context: How do molecular clouds form out of the atomic phase? And what are the relative fractions of carbon in the ionized, atomic and molecular phase? These are questions at the heart of cloud and star formation. Methods: Using multiple observatories from Herschel and SOFIA to APEX and the IRAM 30m telescope, we mapped the ionized, atomic and molecular carbon ([CII]@1900GHz, [CI]@492GHz and C18O(2-1)@220GHz) at high spatial resolution (12"-25") in four young massive infrared dark clouds (IRDCs). Results: The three carbon phases were successfully mapped in all four regions, only in one source the [CII] line remained a non-detection. Both the molecular and atomic phases trace the dense structures well, with [CI] also tracing material at lower column densities. [CII] exhibits diverse morphologies in our sample, from compact to diffuse structures probing the cloud environment. In at least two out of the four regions, we find kinematic signatures strongly indicating that the dense gas filaments have formed out of a dynamically active and turbulent atomic/molecular cloud, potentially from converging gas flows. The atomic-to-molecular carbon gas mass ratios are low between 7% and 12% with the lowest values found toward the most quiescent region. In the three regions where [CII] is detected, its mass is always higher by a factor of a few than that of the atomic carbon. The ionized carbon emission depends as well on the radiation field, however, we also find strong [CII] emission in a region without significant external sources, indicating that other processes, e.g., energetic gas flows can contribute to the [CII] excitation as well.

preprint2014arXiv

Deuteration around the ultracompact HII region Mon R2

The massive star-forming region Mon R2 hosts the closest ultra-compact HII region that can be spatially resolved with current single-dish telescopes. We used the IRAM-30m telescope to carry out an unbiased spectral survey toward two important positions (namely IF and MP2), in order to studying the chemistry of deuterated molecules toward Mon R2. We found a rich chemistry of deuterated species at both positions, with detections of C2D, DCN, DNC, DCO+, D2CO, HDCO, NH2D, and N2D+ and their corresponding hydrogenated species and isotopologs. Our high spectral resolution observations allowed us to resolve three velocity components: the component at 10 km/s is detected at both positions and seems associated with the layer most exposed to the UV radiation from IRS 1; the component at 12 km/s is found toward the IF position and seems related to the molecular gas; finally, a component at 8.5 km/s is only detected toward the MP2 position, most likely related to a low-UV irradiated PDR. We derived the column density of all the species, and determined the deuterium fractions (Dfrac). The values of Dfrac are around 0.01 for all the observed species, except for HCO+ and N2H+ which have values 10 times lower. The values found in Mon R2 are well explained with pseudo-time-dependent gas-phase model in which deuteration occurs mainly via ion-molecule reactions with H2D+, CH2D+ and C2HD+. Finally, the [H13CN]/[HN13C] ratio is very high (~11) for the 10 km/s component, which also agree with our model predictions for an age of ~0.01-0.1 Myr. The deuterium chemistry is a good tool for studying star-forming regions. The low-mass star-forming regions seem well characterized with Dfrac(N2H+) or Dfrac(HCO+), but it is required a complete chemical modeling to date massive star-forming regions, because the higher gas temperature together with the rapid evolution of massive protostars.

preprint2014arXiv

Herschel Survey of Galactic OH+, H2O+, and H3O+: Probing the Molecular Hydrogen Fraction and Cosmic-Ray Ionization Rate

In diffuse interstellar clouds the chemistry that leads to the formation of the oxygen bearing ions OH+, H2O+, and H3O+ begins with the ionization of atomic hydrogen by cosmic rays, and continues through subsequent hydrogen abstraction reactions involving H2. Given these reaction pathways, the observed abundances of these molecules are useful in constraining both the total cosmic-ray ionization rate of atomic hydrogen (zeta_H) and molecular hydrogen fraction, f(H2). We present observations targeting transitions of OH+, H2O+, and H3O+ made with the Herschel Space Observatory along 20 Galactic sight lines toward bright submillimeter continuum sources. Both OH+ and H2O+ are detected in absorption in multiple velocity components along every sight line, but H3O+ is only detected along 7 sight lines. From the molecular abundances we compute f(H2) in multiple distinct components along each line of sight, and find a Gaussian distribution with mean and standard deviation 0.042+-0.018. This confirms previous findings that OH+ and H2O+ primarily reside in gas with low H2 fractions. We also infer zeta_H throughout our sample, and find a log-normal distribution with mean log(zeta_H)=-15.75, (zeta_H=1.78x10^-16 s^-1), and standard deviation 0.29 for gas within the Galactic disk, but outside of the Galactic center. This is in good agreement with the mean and distribution of cosmic-ray ionization rates previously inferred from H3+ observations. Ionization rates in the Galactic center tend to be 10--100 times larger than found in the Galactic disk, also in accord with prior studies.

preprint2013arXiv

Kinematics of the ionized-to-neutral interfaces in Monoceros R2

Context. Monoceros R2 (Mon R2), at a distance of 830 pc, is the only ultra-compact H ii region (UC H ii) where its associated photon-dominated region (PDR) can be resolved with the Herschel Space Observatory. Aims. Our aim is to investigate observationally the kinematical patterns in the interface regions (i.e., the transition from atomic to molecular gas) associated with Mon R2. Methods. We used the HIFI instrument onboard Herschel to observe the line profiles of the reactive ions CH+, OH+ and H2O+ toward different positions in Mon R2. We derive the column density of these molecules and compare them with gas-phase chemistry models. Results. The reactive ion CH+ is detected both in emission (at central and red-shifted velocities) and in absorption (at blue-shifted velocities). OH+ is detected in absorption at both blue- and red-shifted velocities, with similar column densities. H2O+ is not detected at any of the positions, down to a rms of 40 mK toward the molecular peak. At this position, we find that the OH+ absorption originates in a mainly atomic medium, and therefore is associated with the most exposed layers of the PDR. These results are consistent with the predictions from photo-chemical models. The line profiles are consistent with the atomic gas being entrained in the ionized gas flow along the walls of the cavity of the H ii region. Based on this evidence, we are able to propose a new geometrical model for this region. Conclusions. The kinematical patterns of the OH+ and CH+ absorption indicate the existence of a layer of mainly atomic gas for which we have derived, for the first time, some physical parameters and its dynamics.

preprint2013arXiv

Large scale IRAM 30m CO-observations in the giant molecular cloud complex W43

We aim to give a full description of the distribution and location of dense molecular clouds in the giant molecular cloud complex W43. It has previously been identified as one of the most massive star-forming regions in our Galaxy. To trace the moderately dense molecular clouds in the W43 region, we initiated an IRAM 30m large program, named W43-HERO, covering a large dynamic range of scales (from 0.3 to 140 pc). We obtained on-the-fly-maps in 13CO (2-1) and C18O (2-1) with a high spectral resolution of 0.1 km/s and a spatial resolution of 12". These maps cover an area of ~1.5 square degrees and include the two main clouds of W43, as well as the lower density gas surrounding them. A comparison with Galactic models and previous distance calculations confirms the location of W43 near the tangential point of the Scutum arm at a distance from the Sun of approximately 6 kpc. The resulting intensity cubes of the observed region are separated into sub-cubes, centered on single clouds which are then analyzed in detail. The optical depth, excitation temperature, and H2 column density maps are derived out of the 13CO and C18O data. These results are then compared with those derived from Herschel dust maps. The mass of a typical cloud is several 10^4 solar masses while the total mass in the dense molecular gas (>100 cm^-3) in W43 is found to be about 1.9e6 solar masses. Probability distribution functions obtained from column density maps derived from molecular line data and Herschel imaging show a log-normal distribution for low column densities and a power-law tail for high densities. A flatter slope for the molecular line data PDF may imply that those selectively show the gravitationally collapsing gas.

preprint2013arXiv

Spatial distribution of small hydrocarbons in the neighborhood of the Ultra Compact HII region Monoceros R2

We study the chemistry of small hydrocarbons in the photon-dominated regions (PDRs) associated with the ultra-compact HII region Mon R2. Our goal is to determine the variations of the abundance of small hydrocarbons in a high-UV irradiated PDR and investigate their chemistry. We present an observational study of CH, CCH and c-C$_3$H$_2$ in Mon R2 combining data obtained with the IRAM 30m telescope and Herschel. We determine the column densities of these species, and compare their spatial distributions with that of polycyclic aromatic hydrocarbon (PAH). We compare the observational results with different chemical models to explore the relative importance of gas-phase, grain-surface and time-dependent chemistry in these environments. The emission of the small hydrocarbons show different patterns. The CCH emission is extended while CH and c-C$_3$H$_2$ are concentrated towards the more illuminated layers of the PDR. The ratio of the column densities of c-C$_3$H$_2$ and CCH shows spatial variations up to a factor of a few, increasing from $N(c-C$_3$H$_2$)/N(CCH)\approx0.004$ in the envelope to a maximum of $\sim0.015-0.029$ towards the 8$μ$m emission peak. Comparing these results with other galactic PDRs, we find that the abundance of CCH is quite constant over a wide range of G$_0$, whereas the abundance of c-C$_3$H$_2$ is higher in low-UV PDRs. In Mon R2, the gas-phase steady-state chemistry can account relatively well for the abundances of CH and CCH in the most exposed layers of the PDR, but falls short by a factor of 10 to reproduce c-C$_3$H$_2$. In the molecular envelope, time-dependent effects and grain surface chemistry play a dominant role in determining the hydrocarbons abundances. Our study shows that CCH and c-C$_3$H$_2$ present a complex chemistry in which UV photons, grain-surface chemistry and time dependent effects contribute to determine their abundances.

preprint2013arXiv

The chemistry of ions in the Orion Bar I. - CH+, SH+, and CF+: The effect of high electron density and vibrationally excited H2 in a warm PDR surface

The abundances of interstellar CH+ and SH+ are not well understood as their most likely formation channels are highly endothermic. Using data from Herschel, we study the formation of CH+ and SH+ in a typical high UV-illumination photon-dominated region (PDR), the Orion Bar. Herschel/HIFI provides velocity-resolved data of CH+ 1-0 and 2-1 and three hyperfine transitions of SH+. Herschel/PACS provides information on the excitation and spatial distribution of CH+ (up to J=6-5). The widths of the CH+ 2-1 and 1-0 transitions are of ~5 km/s, significantly broader than the typical width of dense gas tracers in the Orion Bar (2-3 km/s) and are comparable to the width of tracers of the interclump medium such as C+ and HF. The detected SH+ transitions are narrower compared to CH+ and have line widths of 3 km/s, indicating that SH+ emission mainly originates in denser condensations. Non-LTE radiative transfer models show that electron collisions affect the excitation of CH+ and SH+, and that reactive collisions need to be taken into account to calculate the excitation of CH+. Comparison to PDR models shows that CH+ and SH+ are tracers of the warm surface region (AV<1.5) of the PDR with temperatures between 500-1000 K. We have also detected the 5-4 transition of CF+ (FWHM=1.9 km/s) with an intensity that is consistent with previous observations of lower-J CF+ transitions toward the Orion Bar. A comparison to PDR models indicate that the internal vibrational energy of H2 can explain the formation of CH+ for typical physical conditions in the Orion Bar near the ionization front. H2 vibrational excitation is the most likely explanation of SH+ formation as well. The abundance ratios of CH+ and SH+ trace the destruction paths of these ions, and through that, indirectly, the ratios of H, H2 and electron abundances as a function of depth into the cloud.

preprint2013arXiv

The Effects of Radiative Transfer on the PDFs of Molecular MHD Turbulence

We study the effects of radiative transfer on the Probability Distribution Functions (PDFs) of simulations of magnetohydrodynamic turbulence in the widely studied $^{13}$CO 2-1 transition. We find that the integrated intensity maps generally follow a log-normal distribution, with the cases that have $τ\approx 1$ best matching the PDF of the column density. We fit a 2D variance-sonic Mach number relationship to our logarithmic PDFs of the form $σ_{ln(Σ/Σ_0)}^2=A\times ln(1+b^2{\cal M}_s^2)$ and find that, for parameter $b=1/3$, parameter $A$ depends on the radiative transfer environment. We also explore the variance, skewness, and kurtosis of the linear PDFs finding that higher moments reflect both higher sonic Mach number and lower optical depth. Finally, we apply the Tsallis incremental PDF function and find that the fit parameters depend on both Mach numbers, but also are sensitive to the radiative transfer parameter space, with the $τ\approx 1$ case best fitting the incremental PDF of the true column density. We conclude that, for PDFs of low optical depth cases, part of the gas is always sub-thermally excited so that the spread of the line intensities exceeds the spread of the underlying column densities and hence the PDFs do not reflect the true column density. Similarly, PDFs of optically thick cases are dominated by the velocity dispersion and therefore do not represent the true column density PDF. Thus, in the case of molecules like carbon monoxide, the dynamic range of intensities, structures observed and consequently, the observable PDFs, are less determined by turbulence and more-often determined by radiative transfer effects.

preprint2013arXiv

The Turbulence Power Spectrum in Optically Thick Interstellar Clouds

The Fourier power spectrum is one of the most widely used statistical tools to analyze the nature of magnetohydrodynamic turbulence in the interstellar medium. Lazarian & Pogosyan (2004) predicted that the spectral slope should saturate to -3 for an optically thick medium and many observations exist in support of their prediction. However, there have not been any numerical studies to-date testing these results. We analyze the spatial power spectrum of MHD simulations with a wide range of sonic and Alfvénic Mach numbers, which include radiative transfer effects of the $^{13}$CO transition. We confirm numerically the predictions of Lazarian & Pogosyan (2004) that the spectral slope of line intensity maps of an optically thick medium saturates to -3. Furthermore, for very optically thin supersonic CO gas, where the density or CO abundance values are too low to excite emission in all but the densest shock compressed gas, we find that the spectral slope is shallower than expected from the column density. Finally, we find that mixed optically thin/thick CO gas, which has average optical depths on order of unity, shows mixed behavior: for super-Alfvénic turbulence, the integrated intensity power spectral slopes generally follow the same trend with sonic Mach number as the true column density power spectrum slopes. However, for sub-Alfvénic turbulence the spectral slopes are steeper with values near -3 which are similar to the very optically thick regime.

preprint2012arXiv

Full SED fitting with the KOSMA-τ PDR code - I. Dust modelling

We revised the treatment of interstellar dust in the KOSMA-τ PDR model code to achieve a consistent description of the dust-related physics in the code. The detailed knowledge of the dust properties is then used to compute the dust continuum emission together with the line emission of chemical species. We coupled the KOSMA-τ PDR code with the MCDRT (multi component dust radiative transfer) code to solve the frequency-dependent radiative transfer equations and the thermal balance equation in a dusty clump under the assumption of spherical symmetry, assuming thermal equilibrium in calculating the dust temperatures, neglecting non-equilibrium effects. We updated the calculation of the photoelectric heating and extended the parametrization range for the photoelectric heating toward high densities and UV fields. We revised the computation of the H2 formation on grain surfaces to include the Eley-Rideal effect, thus allowing for high-temperature H2 formation. We demonstrate how the different optical properties, temperatures, and heating and cooling capabilities of the grains influence the physical and chemical structure of a model cloud. The most influential modification is the treatment of H2 formation on grain surfaces that allows for chemisorption. This increases the total H2 formation significantly and the connected H2 formation heating provides a profound heating contribution in the outer layers of the model clumps. The contribution of PAH surfaces to the photoelectric heating and H2 formation provides a boost to the temperature of outer cloud layers, which is clearly traced by high-J CO lines. Increasing the fraction of small grains in the dust size distribution results in hotter gas in the outer cloud layers caused by more efficient heating and cooler cloud centers, which is in turn caused by the more efficient FUV extinction.

preprint2012arXiv

Herschel / HIFI observations of CO, H2O and NH3 in Mon R2

Context. Mon R2 is the only ultracompact HII region (UCHII) where the associated photon-dominated region (PDR) can be resolved with Herschel. Due to its brightness and proximity, it is the best source to investigate the chemistry and physics of highly UV-irradiated PDRs. Aims. Our goal is to estimate the abundance of H2O and NH3 in this region and investigate their origin. Methods. We present new observations obtained with HIFI and the IRAM-30m telescope. Using a large velocity gradient approach, we model the line intensities and derive an average abundance of H2O and NH3 across the region. Finally, we model the line profiles with a non-local radiative transfer model and compare these results with the abundance predicted by the Meudon PDR code. Results. The variations of the line profiles and intensities indicate complex geometrical and kinematical patterns. The H2O lines present a strong absorption at the ambient velocity and emission in high velocity wings towards the HII region. The spatial distribution of the o-H2^18O line shows that the its emission arises in the PDR surrounding the HII region. By modeling the o-H2^18O emission we derive a mean abundance of o-H2O of ~10^-8 relative to H2. The ortho-H2O abundance is however larger, ~1x10^-7, in the high velocity wings. Possible explanations for this larger abundance include an expanding hot PDR and/or an outflow. Ammonia seems to be present only in the envelope with an average abundance of ~2x10^-9 relative to H2. Conclusions. The Meudon PDR code can account for the measured water abundance in the high velocity gas as long as we assume that it originates from a <1 mag hot expanding layer of the PDR, i.e. that the outflow has only a minor contribution to this emission. To explain the abundances in the rest of the cloud the molecular freeze out and grain surface chemistry would need to be included.

preprint2012arXiv

Herschel/HIFI observations of [CII] and [13CII] in PDRs

Recently, we introduced detailed isotopic chemistry into the KOSMA-tau model for photon-dominated regions (PDRs) to give theoretical predictions for the abundance of the carbon isotopologues as a function of PDR parameters. Combined with radiative transfer computations for specific geometries, we estimated the possible intensity ratio of the [CII]/[13CII] lines. Here, we compare these predictions with new observations. We performed Herschel/HIFI observations of the [CII] 158micron line in a number of PDRs. In all sources we observed at least two hyperfine components of the [13CII] transition allowing to determine the [CII]/[13CII] intensity ratio, after some revision of the intrinsic hyperfine ratios. Comparing the intensity ratios with the results from the updated KOSMA-tau model, we identify cases dominated by chemical fractionation and cases dominated by the optical depth of the main isotopic line. An observable enhancement of the [CII]/[13CII] intensity ratio due to chemical fractionation depends mostly on geometry and velocity structure, and less on the gas density and radiation field. In our observations the [CII]/[13CII] ratio for the integrated line intensity was always dominated by the optical depth of the main isotopic line. However, an enhanced intensity ratio is found for particular velocity components in a few sources: the red-shifted material in the ultracompact HII region Mon R2, the wings of the turbulent profile in the Orion Bar, and possibly a blue wing in NGC7023. The mapping of the [13CII] lines in the Orion Bar allows to derive a C+ column density map confirming the temperature stratification of the C+ layer, in agreement with the chemical stratification of the Bar. The C+ column densities for all sources show that at the position of the [CII] peak emission, a dominant fraction of the gas-phase carbon is in the form of C+.

preprint2012arXiv

Modeling and Analysis Generic Interface for eXternal numerical codes (MAGIX)

The modeling and analysis generic interface for external numerical codes (MAGIX) is a model optimizer developed under the framework of the coherent set of astrophysical tools for spectroscopy (CATS) project. The MAGIX package provides a framework of an easy interface between existing codes and an iterating engine that attempts to minimize deviations of the model results from available observational data, constraining the values of the model parameters and providing corresponding error estimates. Many models (and, in principle, not only astrophysical models) can be plugged into MAGIX to explore their parameter space and find the set of parameter values that best fits observational/experimental data. MAGIX complies with the data structures and reduction tools of ALMA (Atacama Large Millimeter Array), but can be used with other astronomical and with non-astronomical data.

preprint2012arXiv

Spectral line survey of the ultracompact HII region Mon R2

Ultracompact (UC) HII regions constitute one of the earliest phases in the formation of a massive star and are characterized by extreme physical conditions (Go>10^5 Habing field and n>10^6 cm^-3). The UC HII Mon R2 is the closest one and therefore an excellent target to study the chemistry in these complex regions. We carried out a 3mm and 1mm spectral survey using the IRAM 30-m telescope towards three positions that represent different physical environments in Mon R2: (i) the ionization front (IF) at (0",0"); two peaks in the molecular cloud (ii) MP1 at the offset (+15",-15") and (iii) MP2 at the farther offset (0",40"). In addition, we carried out extensive modeling to explain the chemical differences between the three observed regions. We detected more than thirty different species. We detected SO+ and C4H suggesting that UV radiation plays an important role in the molecular chemistry of this region. We detected the typical PDR molecules CN, HCN, HCO, C2H, and c-C3H2. While the IF and the MP1 have a chemistry similar to that found in high UV field and dense PDRs like the Orion Bar, the MP2 is more similar to lower UV/density PDRs like the Horsehead nebula. We also detected complex molecules that are not usually found in PDRs (CH3CN, H2CO, HC3N, CH3OH and CH3C2H). Sulfur compounds CS, HCS+, C2S, H2CS, SO and SO2 and the deuterated species DCN and C2D were also identified. [DCN]/[HCN]=0.03 and [C2D]/[C2H]=0.05, are among the highest in warm regions. Our results show that the high UV/dense PDRs present a different chemistry from that of the low UV case. Abundance ratios like [CO+]/[HCO+] or [HCO]/[HCO+] are good diagnostics to differentiate between them. In Mon R2 we have the two classes of PDRs, a high UV PDR towards the IF and the adjacent molecular bar and a low-UV PDR which extends towards the north-west following the border of the cloud.

preprint2012arXiv

The ionized and hot gas in M17 SW: SOFIA/GREAT THz observations of [C II] and 12CO J=13-12

With new THz maps that cover an area of ~3.3x2.1 pc^2 we probe the spatial distribution and association of the ionized, neutral and molecular gas components in the M17 SW nebula. We used the dual band receiver GREAT on board the SOFIA airborne telescope to obtain a 5'.7x3'.7 map of the 12CO J=13-12 transition and the [C II] 158 um fine-structure line in M17 SW and compare the spectroscopically resolved maps with corresponding ground-based data for low- and mid-J CO and [C I] emission. For the first time SOFIA/GREAT allow us to compare velocity-resolved [C II] emission maps with molecular tracers. We see a large part of the [C II] emission, both spatially and in velocity, that is completely non-associated with the other tracers of photon-dominated regions (PDR). Only particular narrow channel maps of the velocity-resolved [C II] spectra show a correlation between the different gas components, which is not seen at all in the integrated intensity maps. These show different morphology in all lines but give hardly any information on the origin of the emission. The [C II] 158 um emission extends for more than 2 pc into the M17 SW molecular cloud and its line profile covers a broader velocity range than the 12CO J=13-12 and [C I] emissions, which we interpret as several clumps and layers of ionized carbon gas within the telescope beam. The high-J CO emission emerges from a dense region between the ionized and neutral carbon emissions, indicating the presence of high-density clumps that allow the fast formation of hot CO in the irradiated complex structure of M17 SW. The [C II] observations in the southern PDR cannot be explained with stratified nor clumpy PDR models.

preprint2012arXiv

The structure of hot gas in Cepheus B

By observing radiation-affected gas in the Cepheus B molecular cloud we probe whether the sequential star formation in this source is triggered by the radiation from newly formed stars. We used the dual band receiver GREAT onboard SOFIA to map [C II] and CO 13--12 and 11--10 in Cep B and compared the spatial distribution and the spectral profiles with complementary ground-based data of low-$J$ transitions of CO isotopes, atomic carbon, and the radio continuum. The interaction of the radiation from the neighboring OB association creates a large photon-dominated region (PDR) at the surface of the molecular cloud traced through the photoevaporation of C^+. Bright internal PDRs of hot gas are created around the embedded young stars, where we detect evidence of the compression of material and local velocity changes; however, on the global scale we find no indications that the dense molecular material is dynamically affected.

preprint2011arXiv

Detection of HF emission from the Orion Bar

The clumpy density structure of photon-dominated regions is well established, but the physical properties of the clumps and of the surrounding interclump medium are only approximately known. The aim of this paper is to constrain the physical and chemical conditions in the Orion Bar, a prototypical nearby photon-dominated region. We present observations of the HF J=1-0 line, which appears in emission toward the Orion Bar, and compare the brightness of the line to non-LTE radiative transfer calculations. The large width of the HF line suggests an origin of the emission in the interclump gas, but collisional excitation by H2 in the interclump gas underpredicts the observed line intensity by factors of 3-5. In contrast, an origin of the line in the dense clumps requires a density of ~10^9 cm^-3, 10-100 times higher than previous estimates, which is unlikely. However, electron impact excitation reproduces our observations for T = 100 K and n(e) = 10 cm^-3, as expected for the interclump gas. We conclude that HF emission is a signpost of molecular gas with a high electron density. Similar conditions may apply to active galactic nuclei where HF also appears in emission.

preprint2011arXiv

Herschel observations of Extra-Ordinary Sources: Methanol as a probe of physical conditions in Orion KL

We have examined methanol emission from Orion KL with of the {\em Herschel}/HIFI instrument, and detected two methanol bands centered at 524 GHz and 1061 GHz. The 524 GHz methanol band (observed in HIFI band 1a) is dominated by the isolated $Δ$J$=$0, K$=-4\rightarrow$-3, v$_t$$=$0 Q branch, and includes 25 E-type and 2 A-type transitions. The 1061 GHz methanol band (observed in HIFI band 4b) is dominated by the $Δ$J$=$0, K$=7\rightarrow$6, v$_t$$=$0 Q branch transitions which are mostly blended. We have used the isolated E-type v$_t$$=$0 methanol transitions to explore the physical conditions in the molecular gas. With HIFI's high velocity resolution, the methanol emission contributed by different spatial components along the line of sight toward Orion KL (hot core, low velocity flow, and compact ridge) can be distinguished and studied separately. The isolated transitions detected in these bands cover a broad energy range (upper state energy ranging from 80 K to 900 K), which provides a unique probe of the thermal structure in each spatial component. The observations further show that the compact ridge is externally heated. These observations demonstrate the power of methanol lines as probes of the physical conditions in warm regions in close proximity to young stars.

preprint2010arXiv

CHESS, Chemical Herschel surveys of star forming regions:Peering into the protostellar shock L1157-B1

The outflow driven by the low-mass class 0 protostar L1157 is the prototype of the so-called chemically active outflows. The bright bowshock B1 in the southern outflow lobe is a privileged testbed of magneto-hydrodynamical (MHD) shock models, for which dynamical and chemical processes are strongly interdependent. We present the first results of the unbiased spectral survey of the L1157-B1 bowshock, obtained in the framework of the key program "Chemical Herschel Surveys of Star Forming Regions" (CHESS). The main aim is to trace the warm and chemically enriched gas and to infer the excitation conditions in the shock region. The CO 5-4 and H2O lines have been detected at high-spectral resolution in the unbiased spectral survey of the HIFI-Band 1b spectral window (555-636 GHz), presented by Codella et al. in this volume. Complementary ground-based observations in the submm window help establish the origin of the emission detected in the main-beam of HIFI, and the physical conditions in the shock.}{Both lines exhibit broad wings, which extend to velocities much higher than reported up to now. We find that the molecular emission arises from two regions with distinct physical conditions: an extended, warm (100K), dense (3e5 cm-3) component at low-velocity, which dominates the water line flux in Band~1; a secondary component in a small region of B1 (a few arcsec) associated with high-velocity, hot (> 400 K) gas of moderate density ((1.0-3.0)e4 cm-3), which appears to dominate the flux of the water line at 179mu observed with PACS. The water abundance is enhanced by two orders of magnitude between the low- and the high-velocity component, from 8e-7 up to 8e-5. The properties of the high-velocity component agree well with the predictions of steady-state C-shock models.

preprint2010arXiv

Detection of interstellar oxidaniumyl: abundant H2O+ towards the star-forming regions DR21, Sgr B2, and NGC6334

We identify a prominent absorption feature at 1115 GHz, detected in first HIFI spectra towards high-mass star-forming regions, and interpret its astrophysical origin. The characteristic hyperfine pattern of the H2O+ ground-state rotational transition, and the lack of other known low-energy transitions in this frequency range, identifies the feature as H2O+ absorption against the dust continuum background and allows us to derive the velocity profile of the absorbing gas. By comparing this velocity profile with velocity profiles of other tracers in the DR21 star-forming region, we constrain the frequency of the transition and the conditions for its formation. In DR21, the velocity distribution of H2O+ matches that of the [CII] line at 158μ\m and of OH cm-wave absorption, both stemming from the hot and dense clump surfaces facing the HII-region and dynamically affected by the blister outflow. Diffuse foreground gas dominates the absorption towards Sgr B2. The integrated intensity of the absorption line allows us to derive lower limits to the H2O+ column density of 7.2e12 cm^-2 in NGC 6334, 2.3e13 cm^-2 in DR21, and 1.1e15 cm^-2 in Sgr B2.

preprint2010arXiv

Detection of OH$^+$ and H$_2$O$^+$ towards Orion~KL

We report observations of the reactive molecular ions OH$^+$, H$_2$O$^+$, and H$_3$O$^+$ towards Orion KL with Herschel/HIFI. All three $N=1-0$ fine-structure transitions of OH$^+$ at 909, 971, and 1033GHz and both fine-structure components of the doublet {\it ortho}-H$_2$O$^+$ $1_{11}-0_{00}$ transition at 1115 and 1139GHz were detected; an upper limit was obtained for H$_3$O$^+$. OH$^+$ and H$_2$O$^+$ are observed purely in absorption, showing a narrow component at the source velocity of 9 kms$^{-1}$, and a broad blueshifted absorption similar to that reported recently for HF and {\it para}-H$_{2}^{18}$O, and attributed to the low velocity outflow of Orion KL. We estimate column densities of OH$^+$ and H$_2$O$^+$ for the 9 km s$^{-1}$ component of $9 \pm 3 \times 10^{12}$cm$^{-2}$ and $7 \pm 2 \times 10^{12}$cm$^{-2}$, and those in the outflow of $1.9 \pm 0.7 \times 10^{13}$cm$^{-2}$ and $1.0 \pm 0.3 \times 10^{13}$cm$^{-2}$. Upper limits of $2.4\times 10^{12}$cm$^{-2}$ and $8.7\times 10^{12}$cm$^{-2}$ were derived for the column densities of {\it ortho} and {\it para}-H$_3$O$^+$ from transitions near 985 and 1657GHz. The column densities of the three ions are up to an order of magnitude lower than those obtained from recent observations of W31C and W49N. The comparatively low column densities may be explained by a higher gas density despite the assumption of a very high ionization rate.

preprint2010arXiv

Gas morphology and energetics at the surface of PDRs: new insights with Herschel observations of NGC 7023

We investigate the physics and chemistry of the gas and dust in dense photon-dominated regions (PDRs), along with their dependence on the illuminating UV field. Using Herschel-HIFI observations, we study the gas energetics in NGC 7023 in relation to the morphology of this nebula. NGC 7023 is the prototype of a PDR illuminated by a B2V star and is one of the key targets of Herschel. Our approach consists in determining the energetics of the region by combining the information carried by the mid-IR spectrum (extinction by classical grains, emission from very small dust particles) with that of the main gas coolant lines. In this letter, we discuss more specifically the intensity and line profile of the 158 micron (1901 GHz) [CII] line measured by HIFI and provide information on the emitting gas. We show that both the [CII] emission and the mid-IR emission from polycyclic aromatic hydrocarbons (PAHs) arise from the regions located in the transition zone between atomic and molecular gas. Using the Meudon PDR code and a simple transfer model, we find good agreement between the calculated and observed [CII] intensities. HIFI observations of NGC 7023 provide the opportunity to constrain the energetics at the surface of PDRs. Future work will include analysis of the main coolant line [OI] and use of a new PDR model that includes PAH-related species.

preprint2010arXiv

Herschel observations of deuterated water towards Sgr B2(M)

Observations of HDO are an important complement for studies of water, because they give strong constraints on the formation processes -- grain surfaces versus energetic process in the gas phase, e.g. in shocks. The HIFI observations of multiple transitions of HDO in Sgr~B2(M) presented here allow the determination of the HDO abundance throughout the envelope, which has not been possible before with ground-based observations only. The abundance structure has been modeled with the spherical Monte Carlo radiative transfer code RATRAN, which also takes radiative pumping by continuum emission from dust into account. The modeling reveals that the abundance of HDO rises steeply with temperature from a low abundance ($2.5\times 10^{-11}$) in the outer envelope at temperatures below 100~K through a medium abundance ($1.5\times 10^{-9}$) in the inner envelope/outer core, at temperatures between 100 and 200~K, and finally a high abundance ($3.5\times 10^{-9}$) at temperatures above 200~K in the hot core.

preprint2010arXiv

Herschel observations of EXtra-Ordinary Sources (HEXOS): Observations of H2O and its isotopologues towards Orion KL

We report the detection of more than 48 velocity-resolved ground rotational state transitions of H2(16)O, H2(18)O, and H2(17)O - most for the first time - in both emission and absorption toward Orion KL using Herschel/HIFI. We show that a simple fit, constrained to match the known emission and absorption components along the line of sight, is in excellent agreement with the spectral profiles of all the water lines. Using the measured H2(18)O line fluxes, which are less affected by line opacity than their H2(16)O counterparts, and an escape probability method, the column densities of H2(18)O associated with each emission component are derived. We infer total water abundances of 7.4E-5, 1.0E-5, and 1.6E-5 for the plateau, hot core, and extended warm gas, respectively. In the case of the plateau, this value is consistent with previous measures of the Orion-KL water abundance as well as those of other molecular outflows. In the case of the hot core and extended warm gas, these values are somewhat higher than water abundances derived for other quiescent clouds, suggesting that these regions are likely experiencing enhanced water-ice sublimation from (and reduced freeze-out onto) grain surfaces due to the warmer dust in these sources.

preprint2010arXiv

Herschel observations of EXtra-Ordinary Sources: The Terahertz spectrum of Orion KL seen at high spectral resolution

We present the first high spectral resolution observations of Orion KL in the frequency ranges 1573.4 - 1702.8 GHz (band 6b) and 1788.4 - 1906.8 GHz (band 7b) obtained using the HIFI instrument on board the Herschel Space Observatory. We characterize the main emission lines found in the spectrum, which primarily arise from a range of components associated with Orion KL including the hot core, but also see widespread emission from components associated with molecular outflows traced by H2O, SO2, and OH. We find that the density of observed emission lines is significantly diminished in these bands compared to lower frequency Herschel/HIFI bands.

preprint2010arXiv

Herschel observations of ortho- and para-oxidaniumyl (H2O+) in spiral arm clouds toward Sgr B2(M)

H2O+ has been observed in its ortho- and para- states toward the massive star forming core Sgr B2(M), located close to the Galactic center. The observations show absorption in all spiral arm clouds between the Sun and Sgr B2. The average o/p ratio of H2O+ in most velocity intervals is 4.8, which corresponds to a nuclear spin temperature of 21 K. The relationship of this spin temperature to the formation temperature and current physical temperature of the gas hosting H2O+ is discussed, but no firm conclusion is reached. In the velocity interval 0 to 60 km/s, an ortho/para ratio of below unity is found, but if this is due to an artifact of contamination by other species or real is not clear.

preprint2010arXiv

Herschel-HIFI observations of high-J CO lines in the NGC 1333 low-mass star-forming region

Herschel-HIFI observations of high-J lines (up to J_u=10) of 12CO, 13CO and C18O are presented toward three deeply embedded low-mass protostars, NGC 1333 IRAS 2A, IRAS 4A, and IRAS 4B, obtained as part of the Water In Star-forming regions with Herschel (WISH) key program. The spectrally-resolved HIFI data are complemented by ground-based observations of lower-J CO and isotopologue lines. The 12CO 10-9 profiles are dominated by broad (FWHM 25-30 km s^-1) emission. Radiative transfer models are used to constrain the temperature of this shocked gas to 100-200 K. Several CO and 13CO line profiles also reveal a medium-broad component (FWHM 5-10 km s^-1), seen prominently in H2O lines. Column densities for both components are presented, providing a reference for determining abundances of other molecules in the same gas. The narrow C18O 9-8 lines probe the warmer part of the quiescent envelope. Their intensities require a jump in the CO abundance at an evaporation temperature around 25 K, thus providing new direct evidence for a CO ice evaporation zone around low-mass protostars.

preprint2010arXiv

Herschel-PACS spectroscopy of the intermediate mass protostar NGC7129 FIRS 2

Aims: We present preliminary results of the first Herschel spectroscopic observations of NGC7129 FIRS2, an intermediate mass star-forming region. We attempt to interpret the observations in the framework of an in-falling spherical envelope. Methods: The PACS instrument was used in line spectroscopy mode (R=1000-5000) with 15 spectral bands between 63 and 185 microns. This provided good detections of 26 spectral lines seen in emission, including lines of H2O, CO, OH, O I, and C II. Results: Most of the detected lines, particularly those of H2O and CO, are substantially stronger than predicted by the spherical envelope models, typically by several orders of magnitude. In this paper we focus on what can be learned from the detected CO emission lines. Conclusions: It is unlikely that the much stronger than expected line emission arises in the (spherical) envelope of the YSO. The region hot enough to produce such high excitation lines within such an envelope is too small to produce the amount of emission observed. Virtually all of this high excitation emission must arise in structures such as as along the walls of the outflow cavity with the emission produced by a combination of UV photon heating and/or non-dissociative shocks.

preprint2010arXiv

Herschel/HIFI discovery of interstellar chloronium (H$_2$Cl$^+$)

We report the first detection of chloronium, H$_2$Cl$^+$, in the interstellar medium, using the HIFI instrument aboard the \emph{Herschel} Space Observatory. The $2_{12}-1_{01}$ lines of ortho-H$_2^{35}$Cl$^+$ and ortho-H$_2^{37}$Cl$^+$ are detected in absorption towards NGC~6334I, and the $1_{11}-0_{00}$ transition of para-H$_2^{35}$Cl$^+$ is detected in absorption towards NGC~6334I and Sgr~B2(S). The H$_2$Cl$^+$ column densities are compared to those of the chemically-related species HCl. The derived HCl/H$_2$Cl$^+$ column density ratios, $\sim$1--10, are within the range predicted by models of diffuse and dense Photon Dominated Regions (PDRs). However, the observed H$_2$Cl$^+$ column densities, in excess of $10^{13}$~cm$^{-2}$, are significantly higher than the model predictions. Our observations demonstrate the outstanding spectroscopic capabilities of HIFI for detecting new interstellar molecules and providing key constraints for astrochemical models.

preprint2010arXiv

HIFI observations of warm gas in DR21: Shock versus radiative heating

The molecular gas in the DR21 massive star formation region is known to be affected by the strong UV field from the central star cluster and by a fast outflow creating a bright shock. The relative contribution of both heating mechanisms is the matter of a long debate. By better sampling the excitation ladder of various tracers we provide a quantitative distinction between the different heating mechanisms. HIFI observations of mid-J transitions of CO and HCO+ isotopes allow us to bridge the gap in excitation energies between observations from the ground, characterizing the cooler gas, and existing ISO LWS spectra, constraining the properties of the hot gas. Comparing the detailed line profiles allows to identify the physical structure of the different components. In spite of the known shock-excitation of H2 and the clearly visible strong outflow, we find that the emission of all lines up to > 2 THz can be explained by purely radiative heating of the material. However, the new Herschel/HIFI observations reveal two types of excitation conditions. We find hot and dense clumps close to the central cluster, probably dynamically affected by the outflow, and a more widespread distribution of cooler, but nevertheless dense, molecular clumps.

preprint2010arXiv

Photon Dominated Regions in NGC 3603

Aims: We aim at deriving the excitation conditions of the interstellar gas as well as the local FUV intensities in the molecular cloud surrounding NGC 3603 to get a coherent picture of how the gas is energized by the central stars. Methods: The NANTEN2-4m submillimeter antenna is used to map the [CI] 1-0, 2-1 and CO 4-3, 7-6 lines in a 2' x 2' region around the young OB cluster NGC 3603 YC. These data are combined with C18O 2-1 data, HIRES-processed IRAS 60 and 100 micron maps of the FIR continuum, and Spitzer/IRAC maps. Results: The NANTEN2 observations show the presence of two molecular clumps located south-east and south-west of the cluster and confirm the overall structure already found by previous CS and C18O observations. We find a slight position offset of the peak intensity of CO and [CI], and the atomic carbon appears to be further extended compared to the molecular material. We used the HIRES far-infrared dust data to derive a map of the FUV field heating the dust. We constrain the FUV field to values of χ= 3 - 6 \times 10^3 in units of the Draine field across the clouds. Approximately 0.2 to 0.3 % of the total FUV energy is re-emitted in the [CII] 158 μm cooling line observed by ISO. Applying LTE and escape probability calculations, we derive temperatures (TMM1 = 43 K, TMM2 = 47 K), column densities (N(MM1) = 0.9 \times 10^22 cm^-2, N(MM2) = 2.5 \times 10^22 cm^-2) and densities (n(MM1) = 3 \times 10^3 cm^-3, n(MM2) = 10^3 -10^4 cm^-3) for the two observed molecular clumps MM1 and MM2. Conclusions: The cluster is strongly interacting with the ambient molecular cloud, governing its structure and physical conditions. A stability analysis shows the existence of gravitationally collapsing gas clumps which should lead to star formation. Embedded IR sources have already been observed in the outskirts of the molecular cloud and seem to support our conclusions.

preprint2010arXiv

Reversal of infall in SgrB2(M) revealed by Herschel/HIFI observations of HCN lines at THz frequencies

To investigate the accretion and feedback processes in massive star formation, we analyze the shapes of emission lines from hot molecular cores, whose asymmetries trace infall and expansion motions. The high-mass star forming region SgrB2(M) was observed with Herschel/HIFI (HEXOS key project) in various lines of HCN and its isotopologues, complemented by APEX data. The observations are compared to spherically symmetric, centrally heated models with density power-law gradient and different velocity fields (infall or infall+expansion), using the radiative transfer code RATRAN. The HCN line profiles are asymmetric, with the emission peak shifting from blue to red with increasing J and decreasing line opacity (HCN to H$^{13}$CN). This is most evident in the HCN 12--11 line at 1062 GHz. These line shapes are reproduced by a model whose velocity field changes from infall in the outer part to expansion in the inner part. The qualitative reproduction of the HCN lines suggests that infall dominates in the colder, outer regions, but expansion dominates in the warmer, inner regions. We are thus witnessing the onset of feedback in massive star formation, starting to reverse the infall and finally disrupting the whole molecular cloud. To obtain our result, the THz lines uniquely covered by HIFI were critically important.

preprint2010arXiv

The distribution of water in the high-mass star-forming region NGC 6334I

We present observations of twelve rotational transitions of H2O-16, H2O-18, and H2O-17 toward the massive star-forming region NGC 6334 I, carried out with Herschel/HIFI as part of the guaranteed time key program Chemical HErschel Surveys of Star forming regions (CHESS). We analyze these observations to obtain insights into physical processes in this region. We identify three main gas components (hot core, cold foreground, and outflow) in NGC 6334 I and derive the physical conditions in these components. The hot core, identified by the emission in highly excited lines, shows a high excitation temperature of 200 K, whereas water in the foreground component is predominantly in the ortho- and para- ground states. The abundance of water varies between 4 10^-5 (outflow) and 10^-8 (cold foreground gas). This variation is most likely due to the freeze-out of water molecules onto dust grains. The H2O-18/H2O-17 abundance ratio is 3.2, which is consistent with the O-18/O-17 ratio determined from CO isotopologues. The ortho/para ratio in water appears to be relatively low 1.6(1) in the cold, quiescent gas, but close to the equilibrium value of three in the warmer outflow material (2.5(0.8)).

preprint2010arXiv

The link between molecular cloud structure and turbulence

We aim to better understand how the spatial structure of molecular clouds is governed by turbulence. For that, we study the large-scale spatial distribution of low density molecular gas and search for characteristic length scales. We employ a 35 square degrees 13CO 1-0 molecular line survey of Cygnus X and visual extinction (A_V) maps of 17 Galactic clouds to analyse the spatial structure using the Delta-variance method. This sample contains a large variety of different molecular cloud types with different star forming activity. The Delta-variance spectra obtained from the A_V maps show differences between low-mass star-forming (SF) clouds and massive giant molecular clouds (GMC) in terms of shape of the spectrum and its power-law exponent beta. Low-mass SF clouds have a double-peak structure with characteristic size scales around 1 pc (though with a large scatter around this value) and 4 pc. GMCs show no characteristic scale in the A_V-maps, which can partly be ascribed to a distance effect due to a larger line-of-sight (LOS) confusion. The Delta-variance for Cygnus, determined from the 13CO survey, shows characteristic scales at 4 pc and 40 pc, either reflecting the filament structure and large-scale turbulence forcing or - for the 4 pc scale - the scale below which the 13CO 1-0 line becomes optically thick. Though there are different processes that can introduce characteristic scales, i.e. geometry, decaying turbulence the transition scale from supersonic to subsonic turbulence (the sonic scale), line-of-sight effects and energy injection due to expanding supernova shells, outflows, HII-regions, and the relative contribution of these effects strongly varies from cloud to cloud, it is remarkable that the resulting turbulent structure of molecular clouds shows similar characteristics.

preprint2010arXiv

The origin of the [C II] emission in the S140 PDRs - new insights from HIFI

Using Herschel's HIFI instrument we have observed [C II] along a cut through S140 and high-J transitions of CO and HCO+ at two positions on the cut, corresponding to the externally irradiated ionization front and the embedded massive star forming core IRS1. The HIFI data were combined with available ground-based observations and modeled using the KOSMA-tau model for photon dominated regions. Here we derive the physical conditions in S140 and in particular the origin of [C II] emission around IRS1. We identify three distinct regions of [C II] emission from the cut, one close to the embedded source IRS1, one associated with the ionization front and one further into the cloud. The line emission can be understood in terms of a clumpy model of photon-dominated regions. At the position of IRS1, we identify at least two distinct components contributing to the [C II] emission, one of them a small, hot component, which can possibly be identified with the irradiated outflow walls. This is consistent with the fact that the [C II] peak at IRS1 coincides with shocked H2 emission at the edges of the outflow cavity. We note that previously available observations of IRS1 can be well reproduced by a single-component KOSMA-tau model. Thus it is HIFI's unprecedented spatial and spectral resolution, as well as its sensitivity which has allowed us to uncover an additional hot gas component in the S140 region.

preprint2010arXiv

Variations in H2O+/H2O ratios toward massive star-forming regions

Early results from the Herschel Space Observatory revealed the water cation H2O+ to be an abundant ingredient of the interstellar medium. Here we present new observations of the H2O and H2O+ lines at 1113.3 and 1115.2 GHz using the Herschel Space Observatory toward a sample of high-mass star-forming regions to observationally study the relation between H2O and H2O+ . Nine out of ten sources show absorption from H2O+ in a range of environments: the molecular clumps surrounding the forming and newly formed massive stars, bright high-velocity outflows associated with the massive protostars, and unrelated low-density clouds along the line of sight. Column densities per velocity component of H2 O+ are found in the range of 10^12 to a few 10^13 cm-2 . The highest N(H2O+) column densities are found in the outflows of the sources. The ratios of H2O+/H2O are determined in a range from 0.01 to a few and are found to differ strongly between the observed environments with much lower ratios in the massive (proto)cluster envelopes (0.01-0.1) than in outflows and diffuse clouds. Remarkably, even for source components detected in H2O in emission, H2O+ is still seen in absorption.

preprint2010arXiv

Water production in comet 81P/Wild 2 as determined by Herschel/HIFI

The high spectral resolution and sensitivity of Herschel/HIFI allows for the detection of multiple rotational water lines and accurate determinations of water production rates in comets. In this letter we present HIFI observations of the fundamental 110-101 (557 GHz) ortho and 111-000 (1113 GHz) para rotational transitions of water in comet 81P/Wild 2 acquired in February 2010. We mapped the extent of the water line emission with five point scans. Line profiles are computed using excitation models which include excitation by collisions with electrons and neutrals and solar infrared radiation. We derive a mean water production rate of $1.0 \times 10^{28}$ molecules s$^{-1}$ at a heliocentric distance of 1.61 AU about 20 days before perihelion, in agreement with production rates measured from the ground using observations of the 18-cm OH lines. Furthermore, we constrain the electron density profile and gas kinetic temperature, and estimate the coma expansion velocity by fitting the water line shapes.

preprint2007arXiv

A PDR-Code Comparison Study

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.