Source author record

F. F. S. van der Tak

F. F. S. van der Tak 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

53works
5topics
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

53 published item(s)

preprint2022arXiv

Breaking Orion's Veil bubble with fossil outflows

The role of feedback in the self-regulation of star formation is a fundamental question in astrophysics. The Orion Nebula is the nearest site of ongoing and recent massive star formation. It is a unique laboratory for the study of stellar feedback. Recent SOFIA [CII] 158 $μ$m observations revealed an expanding bubble, the Veil shell, being powered by stellar winds and ionization feedback. We have identified a protrusion-like substructure in the Northwest portion of the Orion Veil Shell that may indicate additional feedback mechanisms that are highly directional. Our goal is to investigate the origin of the protrusion by quantifying its possible driving mechanisms. We use the [CII] 158 $μ$m map of the Orion Nebula obtained with the upGREAT instrument onboard SOFIA. The spectral and spatial resolution of the observations are 0.3 km/s and 16 arcsec, respectively. We consider three possible origins for this protrusion: Fossil outflow cavities created by jets/outflows during the protostellar accretion phase, pre-existing clumpiness in the OMC-1 core, and the stellar wind during the main sequence phase. Based on the energetics and the morphology, we conclude that the northwestern part of the pre-existing cloud was locally perturbed by outflows ejected from massive protostars in the Trapezium cluster. This suggests that the protrusion of the Veil is the result of mechanical rather than radiative feedback. Furthermore, we argue that the location of the protrusion is a suitable place to break the Orion Veil owing to the photo-ablation from the walls of the protrusion. We conclude that the outflows of massive protostars can influence the morphology of the future \hii\,region and even cause breakages in the ionization front. Specifically, the interaction of stellar winds of main-sequence stars with the molecular core pre-processed by the protostellar jet is important.

preprint2022arXiv

Dents in the Veil: Protostellar feedback in Orion

Interest in stellar feedback has recently increased because new studies suggest that radiative and mechanical feedback from young massive stars regulate the physical and chemical composition of the interstellar medium (ISM) significantly. Recent SOFIA [CII] 158 micron observations of the Orion Veil revealed that the expanding bubble is powered by stellar winds and influenced by previously active molecular outflows of ionizing massive stars. We aim to investigate the mechanical feedback on the whole Veil shell by searching for jets/outflows interacting with the Veil shell and determining the origin/driving mechanisms of these collisions. In the light of these findings, as well as the momenta of the dents and their dynamical timescales, we propose that the dents are created by the interaction of collimated jets/outflows from protostars with luminosities ranging from 10$^3$ to 10$^4$ $L_\odot$ indicating B-type stars in the Orion star-forming cloud with the surrounding Veil shell. However, it is challenging to pinpoint the driving stars as they may have moved from the original ejection points of the jets/outflows. We conclude that the dynamics of the expanding Veil shell is influenced not just by the O-type stars in the Trapezium cluster, but also by less massive stars, especially B-type, in the Orion Nebula. Mechanical feedback from protostars with a range of masses appears to play an important role in determining the morphology of [HII] regions and injecting turbulence into the medium.

preprint2020arXiv

A search for radio jets from massive young stellar objects. Association of radio jets with H2O and CH3OH masers

Recent theoretical and observational studies debate the similarities between the formation process of high-mass (>8 Msun) and low-mass stars. The formation of low-mass star formation is directly associated with the presence of disks and jets. According to this scenario, radio jets are expected to be common in high-mass star-forming regions. We aim to increase the number of known radio jets in high-mass star forming regions by searching for radio jet candidates at radio continuum wavelengths. We have used the Karl G. Jansky Very Large Array (VLA) to observe 18 high-mass star-forming regions in the C band (6 cm, ~1.0 arcsec resolution) and K band (1.3 cm, ~0.3 arcsec resolution). We have searched for radio jet candidates by studying the association of radio continuum sources with shock activity signposts. We have identified 7 as the most probable radio jets. The radio luminosity of the radio jet candidates is correlated with the bolometric luminosity and the outflow momentum rate. About 7-36% of the radio jet candidates are associated with non-thermal emission. The radio jet candidates associated with 6.7 GHz CH3OH maser emission are preferentially thermal winds and jets, while a considerable fraction of radio jet candidates associated with H2O masers show non-thermal emission, likely due to strong shocks. Our sample of 18 regions is divided in 8 less evolved, infrared-dark regions and 10 more evolved, infrared-bright regions. We have found that ~71% of the identified radio jet candidates are located in the more evolved regions. Similarly, 25% of the less evolved regions harbor one of the most probable radio jets, while up to 50% of the more evolved regions contain one of these radio jet candidates. This suggests that the detection of radio jets in high-mass star forming regions is larger in slightly more evolved regions.

preprint2020arXiv

Diagnosing the interstellar medium of galaxies with far-infrared emission lines I. The [C II] 158 microns line at z~0

Atomic fine structure lines have been detected in the local Universe and at high redshifts over the past decades. The [C II] emission line at 158 $μ$m is an important observable as it provides constraints on the interstellar medium (ISM) cooling processes. We develop a physically motivated framework to simulate the production of far-infrared line emission from galaxies in a cosmological context. This first paper sets out our methodology and describes its first application, simulating the [C II] 158 $μ$m line emission in the local Universe. We combine the output from EAGLE cosmological hydrodynamical simulations with a multi-phase model of the ISM. Gas particles are divided into three phases: dense molecular gas, neutral atomic gas and diffuse ionised gas (DIG). We estimate the [C II] line emission from the three phases using a set of Cloudy cooling tables. Our results agree with previous findings regarding the contribution of these three ISM phases to the [C II] emission. Our model shows good agreement with the observed ${\rm L_{[C II]}}$-star formation rate (SFR) relation in the local Universe within 0.4 dex scatter. The fractional contribution to the [C II] line from different ISM phases depends on the total SFR and metallicity. The neutral gas phase dominates the [C II] emission in galaxies with $\rm{SFR}\sim0.01$-$1\,\rm{M_{\odot}}\,\rm{yr^{-1}}$, but the ionised phase dominates at lower SFRs. Galaxies above solar metallicity exhibit lower ${\rm L_{[C II]}}$/SFR ratios for the neutral phase. In comparison, the ${\rm L_{[C II]}}$/SFR ratio in the DIG is stable when metallicity varies. We suggest that the reduced size of the neutral clouds, caused by increased SFRs, is the likely cause for the ${\rm L_{[C II]}}$ deficit at high infrared luminosities, although EAGLE simulations do not reach these luminosities at $z=0$.

preprint2020arXiv

Enrichment of the HR 8799 planets by minor bodies and dust

In the Solar System, minor bodies and dust deliver various materials to planetary surfaces. Several exoplanetary systems are known to host inner and outer belts, analogues of the main asteroid belt and the Kuiper belt. We study the possibility that exominor bodies and exodust deliver volatiles and refractories to the exoplanets in the system HR8799 by performing N-body simulations. The model consists of the host star, four giant planets (HR8799 e, d, c, and b), 650000 test particles representing the inner belt, and 1450000 test particles representing the outer belt. Moreover we modelled dust populations that originate from both belts. Within a million years, the two belts evolve towards the expected dynamical structure (also derived in other works), where mean-motion resonances with the planets carve the analogues of Kirkwood gaps. We find that, after this point, the planets suffer impacts by objects from the inner and outer belt at rates that are essentially constant with time, while dust populations do not contribute significantly to the delivery process. We convert the impact rates to volatile and refractory delivery rates using our best estimates of the total mass contained in the belts and their volatile and refractory content. Over their lifetime, the four giant planets receive between $10^{-4}$ and $10^{-3}M_\bigoplus$ of material from both belts. The total amount of delivered volatiles and refractories, ${5\times10^{-3}\textrm{M}_\bigoplus}$, is small compared to the total mass of the planets, $11\times10^{3}\textrm{M}_\bigoplus$. However, if the planets were formed to be volatile-rich, their exogenous enrichment in refractory material may well be significant and observable, for example with JWST-MIRI. If terrestrial planets exist within the snow line of the system, volatile delivery would be an important astrobiological mechanism and may be observable as atmospheric trace gases.

preprint2020arXiv

The subsurface habitability of small, icy exomoons

Assuming our Solar System as typical, exomoons may outnumber exoplanets. If their habitability fraction is similar, they would thus constitute the largest portion of habitable real estate in the Universe. Icy moons in our Solar System, such as Europa and Enceladus, have already been shown to possess liquid water, a prerequisite for life on Earth. We intend to investigate under what circumstances small, icy moons may sustain subsurface oceans and thus be "subsurface habitable". We pay specific attention to tidal heating. We made use of a phenomenological approach to tidal heating. We computed the orbit averaged flux from both stellar and planetary (both thermal and reflected stellar) illumination. We then calculated subsurface temperatures depending on illumination and thermal conduction to the surface through the ice shell and an insulating layer of regolith. We adopted a conduction only model, ignoring volcanism and ice shell convection as an outlet for internal heat. In doing so, we determined at which depth, if any, ice melts and a subsurface ocean forms. We find an analytical expression between the moon's physical and orbital characteristics and the melting depth. Since this expression directly relates icy moon observables to the melting depth, it allows us to swiftly put an upper limit on the melting depth for any given moon. We reproduce the existence of Enceladus' subsurface ocean; we also find that the two largest moons of Uranus (Titania & Oberon) could well sustain them. Our model predicts that Rhea does not have liquid water. Habitable exomoon environments may be found across an exoplanetary system, largely irrespective of the distance to the host star. Small, icy subsurface habitable moons may exist anywhere beyond the snow line. This may, in future observations, expand the search area for extraterrestrial habitable environments beyond the circumstellar habitable zone.

preprint2019arXiv

ALMA and ROSINA detections of phosphorus-bearing molecules: the interstellar thread between star-forming regions and comets

To understand how Phosphorus-bearing molecules are formed in star-forming regions, we have analysed ALMA observations of PN and PO towards the massive star-forming region AFGL 5142, combined with a new analysis of the data of the comet 67P/Churyumov-Gerasimenko taken with the ROSINA instrument onboard Rosetta. The ALMA maps show that the emission of PN and PO arises from several spots associated with low-velocity gas with narrow linewidths in the cavity walls of a bipolar outflow. PO is more abundant than PN in most of the spots, with the PO/PN ratio increasing as a function of the distance to the protostar. Our data favor a formation scenario in which shocks sputter phosphorus from the surface of dust grains, and gas-phase photochemistry induced by UV photons from the protostar allows efficient formation of the two species in the cavity walls. Our analysis of the ROSINA data has revealed that PO is the main carrier of P in the comet, with PO/PN>10. Since comets may have delivered a significant amount of prebiotic material to the early Earth, this finding suggests that PO could contribute significantly to the phosphorus reservoir during the dawn of our planet. There is evidence that PO was already in the cometary ices prior to the birth of the Sun, so the chemical budget of the comet might be inherited from the natal environment of the Solar System, which is thought to be a stellar cluster including also massive stars.

preprint2016arXiv

The evolutionary status of dense cores in the NGC 1333 IRAS 4 region

Protostellar evolution, following the formation of the protostar is becoming reasonably well characterized, but the evolution from a prestellar core to a protostar is not well known, although the first hydrostatic core (FHSC) must be a pivotal step. NGC 1333 IRAS 4C is a potentially very young object, that we directly compare with the nearby Class 0 IRAS 4A and IRAS 4B. Observational constraints are provided by spectral imaging from the JCMT Spectral Legacy Survey (330-373 GHz) and continuum and line observations from CARMA. We present integrated intensity and velocity maps of several species, including CO, H2CO and CH3OH. The velocity of an observed outflow, the degree of CO depletion, the deuteration of DCO+/HCO+ and gas kinetic temperatures are observational signatures that we present. We report differences between the three sources in four aspects: a) the kinetic temperature is much lower towards IRAS 4C, b) the line profiles of the detected species show strong outflow activity towards IRAS 4A and IRAS 4B but not towards IRAS 4C, c) the HCN/HNC is < 1 towards IRAS 4C, which confirms the cold nature of the source, d) the degree of CO depletion and the deuteration are the lowest towards the warmest of the sources, IRAS 4B. IRAS 4C seems to be in a different evolutionary state than the IRAS 4A and IRAS 4B sources. We can probably exclude the FHSC stage due to the relatively low Lsmm/Lbol (~6\%) and we investigate the earliest accretion phase of Class 0 stage and the transition between Class 0 to Class I. Our results do not show a consistent scenario for either case, with the major issue being the absence of outflow activity and the cold nature of IRAS 4C. The number of FHSC candidates in Perseus is ~10 times higher than current models predict, which suggests that the lifespan of these objects is greater than 1000 years, possibly due to an accretion rate lower than 4x10^-5 msun/yr.

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

Infrared dark clouds on the far side of the Galaxy

Context: Infrared dark clouds are the coldest and densest portions of giant molecular clouds. The most massive ones represent some of the most likely birthplaces for the next generation of massive stars in the Milky Way. Because a strong mid-IR background is needed to make them appear in absorption, they are usually assumed to be nearby. Aims: We use THz absorption spectroscopy to solve the distance ambiguity associated with kinematic distances for the IR-dark clouds in the TOP100 ATLASGAL sample, a flux-limited selection of massive clumps in different evolutionary phases of star formation. Methods: The para-H2O ground state transition at 1113.343 GHz, observed with Herschel/HIFI, was used to investigate the occurrence of foreground absorption along the line of sight directly towards infrared-dark clouds. Additional consistency checks were performed using MALT90 and HiGAL archival data and targeted Mopra and APEX spectroscopic observations. Results: We report the first discovery of five IRDCs in the TOP100 lying conclusively at the far kinematic distance, showing that the mere presence of low-contrast mid-IR absorption is not sufficient to unequivocally resolve the near/far ambiguity in favour of the former. All IRDCs are massive and actively forming high-mass stars; four of them also show infall signatures. Conclusions: We give a first estimate of the fraction of dark sources at the far distance (~11% in the TOP100) and describe their appearance and properties. The assumption that all dark clouds lie at the near distance may lead, in some cases, to underestimating masses, sizes, and luminosities, possibly causing clouds to be missed that will form very massive stars and clusters.

preprint2015arXiv

Physical and chemical differentiation of the luminous star-forming region W49A - Results from the JCMT Spectral Legacy Survey

The massive and luminous star-forming region W49A is a well known Galactic candidate to probe the physical conditions and chemistry similar to those expected in external starburst galaxies. We aim to probe the physical and chemical structure of W49A on a spatial scale of ~0.8 pc based on the JCMT Spectral Legacy Survey, which covers the frequency range between 330 and 373 GHz. The wide 2x2 arcminutes field and the high spectral resolution of the HARP instrument on JCMT provides information on the spatial structure and kinematics of the cloud. For species where multiple transitions are available, we estimate excitation temperatures and column densities. We detected 255 transitions corresponding to 60 species in the 330-373 GHz range at the center position of W49A. Excitation conditions can be probed for 16 molecules. The chemical composition suggests the importance of shock-, PDR-, and hot core chemistry. Many molecular lines show a significant spatial extent across the maps including high density tracers (e.g. HCN, HNC, CS, HCO+) and tracers of UV-irradiation (e.g. CN and C2H). Large variations are seen between the sub-regions with mostly blue-shifted emission toward the Eastern tail, mostly red-shifted emission toward the Northern clump, and emission peaking around the expected source velocity toward the South-west clump. A comparison of column density ratios of characteristic species observed toward W49A to Galactic PDRs suggests that while the chemistry toward the W49A center is driven by a combination of UV-irradiation and shocks, UV-irradiation dominates for the Northern Clump, Eastern tail, and South-west clump regions. A comparison to a starburst galaxy and an AGN suggests similar C2H, CN, and H2CO abundances (with respect to the dense gas tracer 34CS) between the ~0.8 pc scale probed for W49A and the >1 kpc regions in external galaxies with global star-formation.

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

Water in star-forming regions with Herschel (WISH) V. The physical conditions in low-mass protostellar outflows revealed by multi-transition water observations

Context: Outflows are an important part of the star formation process as both the result of ongoing active accretion and one of the main sources of mechanical feedback on small scales. Water is the ideal tracer of these effects because it is present in high abundance in various parts of the protostar. Method: We present \textit{Herschel} HIFI spectra of multiple water-transitions towards 29 nearby Class 0/I protostars as part of the WISH Survey. These are decomposed into different Gaussian components, with each related to one of three parts of the protostellar system; quiescent envelope, cavity shock and spot shocks in the jet and at the base of the outflow. We then constrain the excitation conditions present in the two outflow-related components. Results: Water emission is optically thick but effectively thin, with line ratios that do not vary with velocity, in contrast to CO. The physical conditions of the cavity and spot shocks are similar, with post-shock H$_{2}$ densities of order 10$^{5}-$10$^{8}$\,cm$^{-3}$ and H$_{2}$O column densities of order 10$^{16}-$10$^{18}$\,cm$^{-2}$. H$_{2}$O emission originates in compact emitting regions: for the spot shocks these correspond to point sources with radii of order 10-200\,AU, while for the cavity shocks these come from a thin layer along the outflow cavity wall with thickness of order 1-30\,AU. Conclusions: Water emission at the source position traces two distinct kinematic components in the outflow; J shocks at the base of the outflow or in the jet, and C shocks in a thin layer in the cavity wall. Class I sources have similar excitation conditions to Class 0 sources, but generally smaller line-widths and emitting region sizes. We suggest that it is the velocity of the wind driving the outflow, rather than the decrease in envelope density or mass, that is the cause of the decrease in H$_{2}$O intensity between Class 0 and I.

preprint2014arXiv

A necklace of dense cores in the high-mass star forming region G35.20-0.74N: ALMA observations

The present study aims at characterizing the massive star forming region G35.20N, which is found associated with at least one massive outflow and contains multiple dense cores, one of them recently found associated with a Keplerian rotating disk. We used ALMA to observe the G35.20N region in the continuum and line emission at 350 GHz. The observed frequency range covers tracers of dense gas (e.g. H13CO+, C17O), molecular outflows (e.g. SiO), and hot cores (e.g. CH3CN, CH3OH). The ALMA 870 um continuum emission map reveals an elongated dust structure (0.15 pc long and 0.013 pc wide) perpendicular to the large-scale molecular outflow detected in the region, and fragmented into a number of cores with masses 1-10 Msun and sizes 1600 AU. The cores appear regularly spaced with a separation of 0.023 pc. The emission of dense gas tracers such as H13CO+ or C17O is extended and coincident with the dust elongated structure. The three strongest dust cores show emission of complex organic molecules characteristic of hot cores, with temperatures around 200 K, and relative abundances 0.2-2x10^(-8) for CH3CN and 0.6-5x10^(-6) for CH3OH. The two cores with highest mass (cores A and B) show coherent velocity fields, with gradients almost aligned with the dust elongated structure. Those velocity gradients are consistent with Keplerian disks rotating about central masses of 4-18 Msun. Perpendicular to the velocity gradients we have identified a large-scale precessing jet/outflow associated with core B, and hints of an east-west jet/outflow associated with core A. The elongated dust structure in G35.20N is fragmented into a number of dense cores that may form massive stars. Based on the velocity field of the dense gas, the orientation of the magnetic field, and the regularly spaced fragmentation, we interpret this elongated structure as the densest part of a 1D filament fragmenting and forming massive stars.

preprint2014arXiv

A non-equilibrium ortho-to-para ratio of water in the Orion PDR

The ortho-to-para ratio (OPR) of H$_2$O is thought to be sensitive to the temperature of water formation. The OPR of H$_2$O is thus useful to study the formation mechanism of water. We investigate the OPR of water in the Orion PDR (Photon-dominated region), at the Orion Bar and Orion S positions, using data from {\it Herschel}/HIFI. We detect the ground-state lines of ortho- and para-H$_2$$^{18}$O in the Orion Bar and Orion S and we estimate the column densities using LTE and non-LTE methods. Based on our calculations, the ortho-to-para ratio (OPR) in the Orion Bar is 0.1 $-$ 0.5, which is unexpectedly low given the gas temperature of $\sim$ 85 K, and also lower than the values measured for other interstellar clouds and protoplanetary disks. Toward Orion S, our OPR estimate is below 2. This low OPR at 2 positions in the Orion PDR is inconsistent with gas phase formation and with thermal evaporation from dust grains, but it may be explained by photodesorption.

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.

preprint2014arXiv

Observations of water with Herschel/HIFI toward the high-mass protostar AFGL 2591

Water is a sensitive tracer of physical conditions in star-forming regions because of its large abundance variations between hot and cold regions. We use spectrally resolved observations of rotational lines of H$_2$O and its isotopologs with Herschel/HIFI to constrain the physical conditions of the water emitting region toward the high-mass protostar AFGL2591. We use analytical estimates and rotation diagrams to estimate T$_{ex}$ and column densities of H$_2$O of the envelope, the outflow, and a foreground cloud. Furthermore, we use the non-LTE radiative transfer code to estimate the temperature and volume density of the H$_2$O emitting gas. Assuming LTE, we estimate an T$_{ex}$ of 42 K and a column density of 2$\times$10$^{14}$ cm$^{-2}$ for the envelope and 45 K and 4$\times$10$^{13}$ cm$^{-2}$ for the outflow, in beams of 4" and 30", respectively. Non-LTE models indicate a kinetic temperature of 60-230 K and a volume density of 7$\times$10$^6$-10$^8$ cm$^{-3}$ for the envelope, and a kinetic temperature of 70-90 K and a gas density of 10$^7$-10$^8$ cm$^{-3}$ for the outflow. The o/p ratio of the foreground absorption is 1.9$\pm$0.4, suggesting a low temperature. In contrast, the o/p ratio seen in absorption by the outflow is 3.5$\pm$1.0, as expected for warm gas. The water abundance in the envelope is 10$^{-9}$, similar to the low values found for other high- and low-mass protostars, suggesting that this abundance is constant during the embedded phase of high-mass star formation. The water abundance in the outflow is 10$^{-10}$, which is 10$\times$ lower than in the envelope and in the outflows of high- and low-mass protostars. Since beam size effects can only increase this estimate by a factor of 2, we suggest that the water in the outflow is affected by dissociating UV radiation due to the low extinction in the outflow lobe.

preprint2014arXiv

The HIFI spectral survey of AFGL 2591 (CHESS). II. Summary of the survey

This paper presents the richness of submillimeter spectral features in the high-mass star forming region AFGL 2591. As part of the CHESS (Chemical Herschel Survey of Star Forming Regions) Key Programme, AFGL 2591 was observed by the Herschel/HIFI instrument. The spectral survey covered a frequency range from 480 up to 1240 GHz as well as single lines from 1267 to 1901 GHz (i.e. CO, HCl, NH3, OH and [CII]). Rotational and population diagram methods were used to calculate column densities, excitation temperatures and the emission extents of the observed molecules associated with AFGL 2591. The analysis was supplemented with several lines from ground-based JCMT spectra. From the HIFI spectral survey analysis a total of 32 species were identified (including isotopologues). In spite of the fact that lines are mostly quite week, 268 emission and 16 absorption lines were found (excluding blends). Molecular column densities range from 6e11 to 1e19 cm-2 and excitation temperatures range from 19 to 175 K. One can distinguish cold (e.g. HCN, H2S, NH3 with temperatures below 70 K) and warm species (e.g. CH3OH, SO2) in the protostellar envelope.

preprint2013arXiv

A candidate circumbinary Keplerian disk in G35.20-0.74 N: A study with ALMA

We report on ALMA observations of continuum and molecular line emission with 0.4" resolution towards the high-mass star forming region G35.20-0.74 N. Two dense cores are detected in typical hot-core tracers, such as CH3CN, which reveal velocity gradients. In one of these cores, the velocity field can be fitted with an almost edge-on Keplerian disk rotating about a central mass of 18 Msun. This finding is consistent with the results of a recent study of the CO first overtone bandhead emission at 2.3mum towards G35.20-0.74 N. The disk radius and mass are >2500 au and 3 Msun. To reconcile the observed bolometric luminosity (3x10^4 Lsun) with the estimated stellar mass of 18 Msun, we propose that the latter is the total mass of a binary system.

preprint2013arXiv

Dense molecular cocoons in the massive protocluster W3 IRS5: a test case for models of massive star formation

[Context] Two competing models describe the formation of massive stars in objects like the Orion Trapezium. In the turbulent core accretion model, the resulting stellar masses are directly related to the mass distribution of the cloud condensations. In the competitive accretion model, the gravitational potential of the protocluster captures gas from the surrounding cloud for which the individual cluster members compete. [Aims] With high resolution submillimeter observations of the structure, kinematics, and chemistry of the proto-Trapezium cluster W3 IRS5, we aim to determine which mode of star formation dominates. [Methods] We present 354 GHz Submillimeter Array observations at resolutions of 1"-3" (1800-5400 AU) of W3 IRS5. ...... [Results] The observations show five emission peaks (SMM1-5). SMM1 and SMM2 contain massive embedded stars (~20 Msun); SMM3-5 are starless or contain low-mass stars (<8 Msun). The inferred densities are high, >= 10^7 cm^-3, but the core masses are small, 0.2-0.6 Msun. The detected molecular emission reveals four different chemical zones. ...... [Conclusions] The proto-Trapezium cluster W3 IRS5 is an ideal test case to discriminate between models of massive star formation. Either the massive stars accrete locally from their local cores; in this case the small core masses imply that W3 IRS5 is at the very end stages (1000 yr) of infall and accretion, or the stars are accreting from the global collapse of a massive, cluster forming core. We find that the observed masses, densities and line widths observed toward W3 IRS 5 and the surrounding cluster forming core are consistent with the competitive accretion of gas at rates of Macc~10^-4 Msun yr^-1 by the massive young forming stars. ......

preprint2013arXiv

Deuterated water in the solar-type protostars NGC 1333 IRAS 4A and IRAS 4B

Aims. The aim of this paper is to study deuterated water in the solar-type protostars NGC1333 IRAS4A and IRAS4B, to compare their HDO abundance distribution with other star-forming regions, and to constrain their HDO/H2O ratios. Methods. Using the Herschel/HIFI instrument as well as ground-based telescopes, we observed several HDO lines covering a large excitation range (Eup/k=22-168 K) towards these protostars and an outflow position. Non-LTE radiative transfer codes were then used to determine the HDO abundance profiles in these sources. Results. The HDO fundamental line profiles show a very broad component, tracing the molecular outflows, in addition to a narrower emission component and a narrow absorbing component. In the protostellar envelope of NGC1333 IRAS4A, the HDO inner (T>100 K) and outer (T<100 K) abundances with respect to H2 are estimated at 7.5x10^{-9} and 1.2x10^{-11}, respectively, whereas, in NGC1333 IRAS4B, they are 1.0x10^{-8} and 1.2x10^{-10}, respectively. Similarly to the low-mass protostar IRAS16293-2422, an absorbing outer layer with an enhanced abundance of deuterated water is required to reproduce the absorbing components seen in the fundamental lines at 465 and 894 GHz in both sources. This water-rich layer is probably extended enough to encompass the two sources as well as parts of the outflows. In the outflows emanating from NGC1333 IRAS4A, the HDO column density is estimated at about (2-4)x10^{13} cm^{-2}, leading to an abundance of about (0.7-1.9)x10^{-9}. An HDO/H2O ratio between 7x10^{-4} and 9x10^{-2} is derived in the outflows. In the warm inner regions of these two sources, we estimate the HDO/H2O ratios at about 1x10^{-4}-4x10^{-3}. This ratio seems higher (a few %) in the cold envelope of IRAS4A, whose possible origin is discussed in relation to formation processes of HDO and H2O.

preprint2013arXiv

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

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

preprint2013arXiv

Herschel-HIFI observations of high-J CO and isotopologues in star-forming regions: from low- to high-mass

Context; Our understanding of the star formation process has traditionally been confined to certain mass or luminosity boundaries because most studies focus only on low-, intermediate- or high-mass star-forming regions. As part of the "Water In Star-forming regions with Herschel" (WISH) key program, water and other important molecules, such as CO and OH, have been observed in 51 embedded young stellar objects (YSOs). The studied sample covers a range of luminosities from <1 to >10^5 L_sol. Aims; We analyse the CO line emission towards a large sample of protostars in terms of both line intensities and profiles. Methods; Herschel-HIFI spectra of the 12CO 10-9, 13CO 10-9 and C18O 5-4, 9-8 and 10-9 lines are analysed for a sample of 51 YSOs. In addition, JCMT spectra of 12CO 3-2 and C18O 3-2 extend this analysis to cooler gas components. Results; All observed CO and isotopologue spectra show a strong linear correlation between the logarithms of the line and bolometric luminosities across six orders of magnitude on both axes. This suggests that the high-J CO lines primarily trace the amount of dense gas associated with YSOs. This relation can be extended to larger (extragalactic) scales. The majority of the detected 12CO line profiles can be decomposed into a broad and a narrow Gaussian component, while the C18O spectra are mainly fitted with a single Gaussian. A broadening of the line profile is also observed from pre-stellar cores to embedded protostars, which is due mostly to non-thermal motions (turbulence/infall). The widths of the broad 12CO 3-2 and 10-9 velocity components correlate with those of the narrow C18O 9-8 profiles, suggesting that the entrained outflowing gas and envelope motions are related independent of the mass of the protostar. These results indicate that physical processes in protostellar envelopes have similar characteristics across the studied luminosity range.

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 HIFI spectral survey of AFGL2591 (CHESS). I. Highly excited linear rotor molecules in the high-mass protostellar envelope

We aim to reveal the gas energetics in the circumstellar environment of the prototypical high-mass protostellar object AFGL2591 using space-based far-infrared observations of linear rotor molecules. Rotational spectral line signatures of CO, HCO+, CS, HCN and HNC from a 490-1240 GHz survey with Herschel/HIFI, complemented by ground-based JCMT and IRAM 30m spectra, cover transitions with E(up)/k between 5 and ~300 K (750K for 12C16O, using selected frequency settings up to 1850 GHz). The resolved spectral line profiles are used to separate and study various kinematic components. The line profiles show two emission components, the widest and bluest of which is attributed to an approaching outflow and the other to the envelope. We find evidence for progressively more redshifted and wider line profiles from the envelope gas with increasing energy level, qualitatively explained by residual outflow contribution picked up in the systematically decreasing beam size. Integrated line intensities for each species decrease as E(up)/k increases from <50 to 700K. We constrain the following: n(H2)~10^5-10^6 cm^-3 and T~60-200K for the outflow gas; T=9-17K and N(H2)~3x10^21 cm^-2 for a known foreground absorption cloud; N(H2)<10^19 cm^-2 for a second foreground component. Our spherical envelope radiative transfer model systematically underproduces observed line emission at E(up)/k > 150 K for all species. This indicates that warm gas should be added to the model and that the model's geometry should provide low optical depth pathways for line emission from this warm gas to escape, for example in the form of UV heated outflow cavity walls viewed at a favorable inclination angle. Physical and chemical conditions derived for the outflow gas are similar to those in the protostellar envelope, possibly indicating that the modest velocity (<10 km/s) outflow component consists of recently swept-up gas.

preprint2012arXiv

Extended warm and dense gas towards W49A: starburst conditions in our Galaxy?

The star formation rates in starburst galaxies are orders of magnitude higher than in local star-forming regions, and the origin of this difference is not well understood. We use sub-mm spectral line maps to characterize the physical conditions of the molecular gas in the luminous Galactic star-forming region W49A and compare them with the conditions in starburst galaxies. We probe the temperature and density structure of W49A using H_2CO and HCN line ratios over a 2'x2' (6.6x6.6 pc) field with an angular resolution of 15" (~0.8 pc) provided by the JCMT Spectral Legacy Survey. We analyze the rotation diagrams of lines with multiple transitions with corrections for optical depth and beam dilution, and estimate excitation temperatures and column densities. Comparing the observed line intensity ratios with non-LTE radiative transfer models, our results reveal an extended region (about 1'x1', equivalent to ~3x3 pc at the distance of W49A) of warm (> 100 K) and dense (>10^5 cm^-3) molecular gas, with a mass of 2x10^4 - 2x10^5 M_Sun (by applying abundances derived for other regions of massive star-formation). These temperatures and densities in W49A are comparable to those found in clouds near the center of the Milky Way and in starburst galaxies. The highly excited gas is likely to be heated via shocks from the stellar winds of embedded, O-type stars or alternatively due to UV irradiation, or possibly a combination of these two processes. Cosmic rays, X-ray irradiation and gas-grain collisional heating are less likely to be the source of the heating in the case of W49A.

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

Kinematics of the inner thousand AU region around the young massive star AFGL 2591-VLA3: a massive disk candidate?

[Context] Recent detections of disks around young high-mass stars support the idea of massive star formation through accretion rather than coalescence, but the detailed kinematics in the equatorial region of the disk candidates is not well known, which limits our understanding of the accretion process. [Aims] This paper explores the kinematics of the gas around a young massive star with millimeter-wave interferometry to improve our understanding of the formation of massive stars though accretion. [Methods] We use Plateau de Bure interferometric images to probe the environment of the nearby (~1 kpc) and luminous (~20000 Lsun) high-mass (10-16 Msun) young star AFGL 2591-VLA3 in continuum and in lines of HDO, H218O and SO2 in the 115 and 230 GHz bands. Radiative transfer calculations are employed to investigate the kinematics of the source. [Results] At ~0.5" (500 AU) resolution, the line images clearly resolve the velocity field of the central compact source (diameter of ~ 800 AU) and show linear velocity gradients in the northeast-southwest direction. Judging from the disk-outflow geometry, the observed velocity gradient results from rotation and radial expansion in the equatorial region of VLA3. Radiative transfer calculations suggest that the velocity field is consistent with sub-Keplerian rotation plus Hubble-law like expansion. The line profiles of the observed molecules suggest a layered structure, with HDO emission arising from the disk mid-plane, H218O from the warm mid-layer, and SO2 from the upper disk. [Conclusions] We propose AFGL 2591-VLA3 as a new massive disk candidate, with peculiar kinematics. The rotation of this disk is sub-Keplerian, probably due to magnetic braking, while the stellar wind may be responsible for the expansion of the disk. The expansion motion [...]

preprint2012arXiv

Multi-line detection of O2 toward rho Oph A

Models of pure gas-phase chemistry in well-shielded regions of molecular clouds predict relatively high levels of molecular oxygen, O2, and water, H2O. Contrary to expectation, the space missions SWAS and Odin found only very small amounts of water vapour and essentially no O2 in the dense star-forming interstellar medium. Only toward rho Oph A did Odin detect a weak line of O2 at 119 GHz in a beam size of 10 arcmin. A larger telescope aperture such as that of the Herschel Space Observatory is required to resolve the O2 emission and to pinpoint its origin. We use the Heterodyne Instrument for the Far Infrared aboard Herschel to obtain high resolution O2 spectra toward selected positions in rho Oph A. These data are analysed using standard techniques for O2 excitation and compared to recent PDR-like chemical cloud models. The 487.2GHz line was clearly detected toward all three observed positions in rho Oph A. In addition, an oversampled map of the 773.8GHz transition revealed the detection of the line in only half of the observed area. Based on their ratios, the temperature of the O2 emitting gas appears to vary quite substantially, with warm gas (> 50 K) adjacent to a much colder region, where temperatures are below 30 K. The exploited models predict O2 column densities to be sensitive to the prevailing dust temperatures, but rather insensitive to the temperatures of the gas. In agreement with these model, the observationally determined O2 column densities seem not to depend strongly on the derived gas temperatures, but fall into the range N(O2) = (3 to >6)e15/cm^2. Beam averaged O2 abundances are about 5e-8 relative to H2. Combining the HIFI data with earlier Odin observations yields a source size at 119 GHz of about 4 - 5 arcmin, encompassing the entire rho Oph A core.

preprint2012arXiv

The Abundance, Ortho/Para Ratio, and Deuteration of Water in the High-Mass Star Forming Region NGC 6334 I

We present Herschel/HIFI observations of 30 transitions of water isotopologues toward the high-mass star forming region NGC 6334 I. The line profiles of H_2^{16}O, H_2^{17}O, H_2^{18}O, and HDO show a complex pattern of emission and absorption components associated with the embedded hot cores, a lower-density envelope, two outflow components, and several foreground clouds, some associated with the NGC 6334 complex, others seen in projection against the strong continuum background of the source. Our analysis reveals an H2O ortho/para ratio of 3 +/- 0.5 in the foreground clouds, as well as the outflow. The water abundance varies from ~10^{-8} in the foreground clouds and the outer envelope to ~10^{-6} in the hot core. The hot core abundance is two orders of magnitude below the chemical model predictions for dense, warm gas, but within the range of values found in other Herschel/HIFI studies of hot cores and hot corinos. This may be related to the relatively low gas and dust temperature (~100 K), or time dependent effects, resulting in a significant fraction of water molecules still locked up in dust grain mantles. The HDO/H_2O ratio in NGC 6334 I, ~2 10^{-4}, is also relatively low, but within the range found in other high-mass star forming regions.

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

First hyperfine resolved far-infrared OH spectrum from a star-forming region

OH is an important molecule in the H2O chemistry and the cooling budget of star-forming regions. The goal of the Herschel key program `Water in Star-forming regions with Herschel' (WISH) is to study H2O and related species during protostellar evolution. Our aim in this letter is to assess the origin of the OH emission from star-forming regions and constrain the properties of the emitting gas. High-resolution observations of the OH 2Pi1/2 J = 3/2-1/2 triplet at 1837.8 GHz (163.1 micron) towards the high-mass star-forming region W3 IRS 5 with the Heterodyne Instrument for the Far-Infrared (HIFI) on Herschel reveal the first hyperfine velocity-resolved OH far-infrared spectrum of a star-forming region. The line profile of the OH emission shows two components: a narrow component (FWHM approx. 4-5 km/s) with partially resolved hyperfine structure resides on top of a broad (FWHM approx. 30 km/s) component. The narrow emission agrees well with results from radiative transfer calculations of a spherical envelope model for W3 IRS 5 with a constant OH abundance of approx. 8e-9. Comparison with H2O yields OH/H2O abundance ratios of around 1e-3 for T > 100 K and around unity for T < 100K, consistent with the current picture of the dense cloud chemistry with freeze-out and photodesorption. The broad component is attributed to outflow emission. An abundance ratio of OH/H2O > 0.028 in the outflow is derived from comparison with results of water line modeling. This ratio can be explained by a fast J-type shock or a slower UV-irradiated C-type shock.

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.

preprint2011arXiv

The JCMT Spectral Legacy Survey: physical structure of the molecular envelope of the high-mass protostar AFGL2591

The understanding of the formation process of massive stars (>8 Msun) is limited, due to theoretical complications and observational challenges. We investigate the physical structure of the large-scale (~10^4-10^5 AU) molecular envelope of the high-mass protostar AFGL2591 using spectral imaging in the 330-373 GHz regime from the JCMT Spectral Legacy Survey. Out of ~160 spectral features, this paper uses the 35 that are spatially resolved. The observed spatial distributions of a selection of six species are compared with radiative transfer models based on a static spherically symmetric structure, a dynamic spherical structure, and a static flattened structure. The maps of CO and its isotopic variations exhibit elongated geometries on scales of ~100", and smaller scale substructure is found in maps of N2H+, o-H2CO, CS, SO2, CCH, and methanol lines. A velocity gradient is apparent in maps of all molecular lines presented here, except SO, SO2, and H2CO. We find two emission peaks in warm (Eup~200K) methanol separated by 12", indicative of a secondary heating source in the envelope. The spherical models are able to explain the distribution of emission for the optically thin H13CO+ and C34S, but not for the optically thick HCN, HCO+, and CS, nor for the optically thin C17O. The introduction of velocity structure mitigates the optical depth effects, but does not fully explain the observations, especially in the spectral dimension. A static flattened envelope viewed at a small inclination angle does slightly better. We conclude that a geometry of the envelope other than an isotropic static sphere is needed to circumvent line optical depth effects. We propose that this could be achieved in envelope models with an outflow cavity and/or inhomogeneous structure at scales smaller than ~10^4 AU. The picture of inhomogeneity is supported by observed substructure in at least six species.

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

First detection of ND in the solar-mass protostar IRAS16293-2422

In the past decade, much progress has been made in characterising the processes leading to the enhanced deuterium fractionation observed in the ISM and in particular in the cold, dense parts of star forming regions such as protostellar envelopes. Very high molecular D/H ratios have been found for saturated molecules and ions. However, little is known about the deuterium fractionation in radicals, even though simple radicals often represent an intermediate stage in the formation of more complex, saturated molecules. The imidogen radical NH is such an intermediate species for the ammonia synthesis in the gas phase. Herschel/HIFI represents a unique opportunity to study the deuteration and formation mechanisms of such species, which are not observable from the ground. We searched here for the deuterated radical ND in order to determine the deuterium fractionation of imidogen and constrain the deuteration mechanism of this species. We observed the solar-mass Class 0 protostar IRAS16293-2422 with the heterodyne instrument HIFI as part of the Herschel key programme CHESS (Chemical HErschel Surveys of Star forming regions). The deuterated form of the imidogen radical ND was detected and securely identified with 2 hyperfine component groups of its fundamental transition in absorption against the continuum background emitted from the nascent protostar. The 3 groups of hyperfine components of its hydrogenated counterpart NH were also detected in absorption. We derive a very high deuterium fractionation with an [ND]/[NH] ratio of between 30 and 100%. The deuterium fractionation of imidogen is of the same order of magnitude as that in other molecules, which suggests that an efficient deuterium fractionation mechanism is at play. We discuss two possible formation pathways for ND, by means of either the reaction of N+ with HD, or deuteron/proton exchange with NH.

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 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

Herschel/HIFI observations of spectrally resolved methylidyne signatures toward the high-mass star-forming core NGC6334I

In contrast to extensively studied dense star-forming cores, little is known about diffuse gas surrounding star-forming regions. We study molecular gas in the high-mass star-forming region NGC6334I, which contains diffuse, quiescent components that are inconspicuous in widely used molecular tracers such as CO. We present Herschel/HIFI observations of CH toward NGC6334I observed as part of the CHESS key program. HIFI resolves the hyperfine components of its J=3/2-1/2 transition, observed in both emission and absorption. The CH emission appears close to the systemic velocity of NGC6334I, while its measured linewidth of 3 km/s is smaller than previously observed in dense gas tracers such as NH3 and SiO. The CH abundance in the hot core is 7 10^-11, two to three orders of magnitude lower than in diffuse clouds. While other studies find distinct outflows in, e.g., CO and H2O toward NGC6334I, we do not detect outflow signatures in CH. To explain the absorption signatures, at least two absorbing components are needed at -3.0 and +6.5 km/s with N(CH)=7 10^13 and 3 10^13 cm^-2. Two additional absorbing clouds are found at +8.0 and 0.0 km/s, both with N(CH)=2 10^13 cm^-2. Turbulent linewidths for the four absorption components vary between 1.5 and 5.0 km/s in FWHM. We constrain physical properties of our CH clouds by matching our CH absorbers with other absorption signatures. In the hot core, molecules such as H2O and CO trace gas that is heated and dynamically influenced by outflow activity, whereas CH traces more quiescent material. The four CH absorbers have column densities and turbulent properties consistent with diffuse clouds: two are located near NGC6334, and two are unrelated foreground clouds. Local density and dynamical effects influence the chemical composition of physical components of NGC6334, causing some components to be seen in CH but not in other tracers, and vice versa.

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

Magnetic field structure in a high-mass outflow/disk system

To characterize the magnetic field structure of the outflow and core region within a prototypical high-mass star-forming region, we analyzed polarized CO(3-2) - for the first time observed with the Submillimeter Array -- as well as 880mum submm continuum emission from the high-mass outflow/disk system IRAS18089-1732. Both emission features with polarization degrees at a few percent level indicate that the magnetic field structure is largely aligned with the outflow/jet orientation from the small core scales to the larger outflow scales. Although quantitative estimates are crude, the analysis indicates that turbulent energy dominates over magnetic energy. The data also suggest a magnetic field strength increase from the lower-density envelope to the higher-density core.

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 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

Tracing early evolutionary stages of high-mass star formation with molecular lines

Despite its major role in the evolution of the interstellar medium, the formation of high-mass stars (M > 10 Msol) is still poorly understood. Two types of massive star cluster precursors, the so-called Massive Dense Cores (MDCs), have been observed, which differ in their mid-infrared brightness. The origin of this difference is not established and could be the result of evolution, density, geometry differences, or a combination of these. We compare several molecular tracers of physical conditions (hot cores, shocks) observed in a sample of mid-IR weak emitting MDCs with previous results obtained in a sample of exclusively mid-IR bright MDCs. The aim is to understand the differences between these two types of object. We present single-dish observations of HDO, H2O-18, SO2 and CH3OH lines at lambda = 1.3 - 3.5 mm. We study line profiles and estimate abundances of these molecules, and use a partial correlation method to search for trends in the results. The detection rates of thermal emission lines are found to be very similar between mid-IR quiet and bright objects. The abundances of H2O, HDO (1E-13 to 1E-9 in the cold outer envelopes), SO2 and CH3OH differ from source to source but independently of their mid-IR flux. In contrast, the methanol class I maser emission, a tracer of outflow shocks, is found to be strongly anti-correlated with the 12 micron source brightnesses. The enhancement of the methanol maser emission in mid-IR quiet MDCs may indicate a more embedded nature. Since total masses are similar between the two samples, we suggest that the matter distribution is spherical around mid-IR quiet sources but flattened around mid-IR bright ones. In contrast, water emission is associated with objects containing a hot molecular core, irrespective of their mid-IR brightness. These results indicate that the mid-IR brightness of MDCs is an indicator of their evolutionary stage.

preprint2010arXiv

W49A: A starburst triggered by expanding shells

W49A is a giant molecular cloud which harbors some of the most luminous embedded clusters in the Galaxy. However, the explanation for this starburst-like phenomenon is still under debate. Methods. We investigated large-scale Spitzer mid-infrared images together with a Galatic Ring Survey 13CO J = 1-0 image, complemented with higher resolution (~ 11 arcsec) 13CO J = 2-1 and C18O J = 2-1 images over a ~ 15 x 13 pc^2 field obtained with the IRAM 30m telescope. Two expanding shells have been identified in the mid-infrared images, and confirmed in the position-velocity diagrams made from the 13CO J = 2-1 and C18O J = 2-1 data. The mass of the averaged expanding shell, which has an inner radius of ~ 3.3 pc and a thickness of ~ 0.41 pc, is about 1.9 x 10^4 M*. The total kinetic energy of the expanding shells is estimated to be ~ 10^49 erg which is probably provided by a few massive stars, whose radiation pressure and/or strong stellar winds drive the shells. The expanding shells are likely to have a common origin close to the two ultracompact Hii regions (source O and source N), and their expansion speed is estimated to be ~ 5 km/s, resulting in an age of ~ 3-7 x 10^5 years. In addition, on larger (~ 35 x 50 pc^2) scales, remnants of two gas ejections have been identified in the 13CO J = 1 - 0 data. Both ejections seem to have the same center as the expanding shells with a total energy of a few times 10^50 erg. The main driving mechanism for the gas ejections is unclear, but likely related to the mechanism which triggers the starburst in W49A.

preprint2010arXiv

Water abundances in high-mass protostellar envelopes: Herschel observations with HIFI

We derive the dense core structure and the water abundance in four massive star-forming regions which may help understand the earliest stages of massive star formation. We present Herschel-HIFI observations of the para-H2O 1_11-0_00 and 2_02-1_11 and the para-H2-18O 1_11-0_00 transitions. The envelope contribution to the line profiles is separated from contributions by outflows and foreground clouds. The envelope contribution is modelled using Monte-Carlo radiative transfer codes for dust and molecular lines (MC3D and RATRAN), with the water abundance and the turbulent velocity width as free parameters. While the outflows are mostly seen in emission in high-J lines, envelopes are seen in absorption in ground-state lines, which are almost saturated. The derived water abundances range from 5E-10 to 4E-8 in the outer envelopes. We detect cold clouds surrounding the protostar envelope, thanks to the very high quality of the Herschel-HIFI data and the unique ability of water to probe them. Several foreground clouds are also detected along the line of sight. The low H2O abundances in massive dense cores are in accordance with the expectation that high densities and low temperatures lead to freeze-out of water on dust grains. The spread in abundance values is not clearly linked to physical properties of the sources.

preprint2009arXiv

Probing X-ray irradiation in the nucleus of NGC 1068 with observations of high-J lines of dense gas tracers

With the incorporation of high-J molecular lines, we aim to constrain the physical conditions of the dense gas in the central region of the Seyfert 2 galaxy NGC 1068 and to determine signatures of the AGN or the starburst contribution. We used the James Clerk Maxwell Telescope to observe the J=4-3 transition of HCN, HNC, and HCO+, as well as the CN N_J=2_{5/2}-1_{3/2} and N_J=3_{5/2}-2_{5/2}, in NGC 1068. We estimate the excitation conditions of HCN, HNC, and CN, based on the line intensity ratios and radiative transfer models. We find that the bulk emission of HCN, HNC, CN, and the high-J HCO+ emerge from dense gas n(H_2)>=10^5 cm^-3). However, the low-J HCO+ lines (dominating the HCO+ column density) trace less dense (n(H_2)<10^5 cm^-3) and colder (T_K<=20 K) gas, whereas the high-J HCO+ emerges from warmer (>30 K) gas than the other molecules. The HCO+ J=4-3 line intensity, compared with the lower transition lines and with the HCN J=4-3 line, support the influence of a local XDR environment. The estimated N(CN)/N(HCN)~1-4 column density ratios are indicative of an XDR/AGN environment with a possible contribution of grain-surface chemistry induced by X-rays or shocks.

preprint2009arXiv

The ISO Long Wavelength Spectrometer line spectrum of VY Canis Majoris and other oxygen-rich evolved stars

The far-infrared spectra of circumstellar envelopes around various oxygen-rich stars were observed using the ISO Long Wavelength Spectrometer (LWS). We have examined high signal-to-noise ISO LWS observations of the luminous supergiant star, VY CMa, with the aim of identifying all of the spectral lines. By paying particular attention to water lines, we aim to separate the lines due to other species, in particular, to prepare for forthcoming observations that will cover the same spectral range using Herschel PACS and at higher spectral resolution using Herschel HIFI and SOFIA. We have developed a fitting method to account for blended water lines using a simple weighting scheme to distribute the flux. We have applied this approach to several other stars which we compare with VY CMa. We present line fluxes for the unblended H2O and CO lines, and present detections of several possible nu_2=1 vibrationally excited water lines. We also identify blended lines of OH, one unblended and several blended lines of NH3, and one possible detection of H3O+.

preprint2004arXiv

An atomic and molecular database for analysis of submillimetre line observations

Atomic and molecular data for the transitions of a number of astrophysically interesting species are summarized, including energy levels, statistical weights, Einstein A-coefficients and collisional rate coefficients. Available collisional data from quantum chemical calculations and experiments are extrapolated to higher energies. These data, which are made publically available through the WWW at http://www.strw.leidenuniv.nl/~moldata, are essential input for non-LTE line radiative transfer programs. An online version of a computer program for performing statistical equilibrium calculations is also made available as part of the database. Comparisons of calculated emission lines using different sets of collisional rate coefficients are presented. This database should form an important tool in analyzing observations from current and future (sub)millimetre and infrared telescopes.