Researcher profile

A. G. G. M. Tielens

A. G. G. M. Tielens contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

astro-ph.GA astro-ph.SR astro-ph.CO astro-ph astro-ph.EP physics.chem-ph astro-ph.IM

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author

59works

0followers

7topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

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.

preprint2022arXiv

Dust Formation in Astrophysical Environments: The Importance of Kinetics

Astronomical observations and analysis of stardust isolated from meteorites have revealed a highly diverse interstellar and circumstellar grain inventory, including a wide range of amorphous materials and crystalline compounds (silicates and carbon). This diversity reflects the wide range of stellar sources injecting solids into the interstellar medium each with its own physical characteristics such as density, temperature and elemental composition and highlights the importance of kinetics rather than thermodynamics in the formation of these compounds. Based upon the extensive literature on soot formation in terrestrial settings, detailed kinetic pathways have been identified for the formation of carbon dust in C-rich stellar ejecta. These have been incorporated in astronomical models for these environments. In recent years, the chemical routes in the nucleation of oxides and silicates have been the focus of much astronomical research. These aspects of stardust formation will be reviewed and lessons for dust formation in planetary atmospheres will be drawn with the emphasis on the influence of kinetics on the characteristics and structure of dust in these environments.

preprint2022arXiv

FEEDBACK from the NGC7538 HII region

Context: How do expanding HII regions interact with their environmental cloud? This is one of the central questions driving the SOFIA legacy program FEEDBACK. Here, we present a case study toward the prototypical H{\sc ii} region NGC7538. Methods: With SOFIA we mapped an area of ~210'^2 around NGC7538 in the [CII] line at 1.9THz. Complementary observed atomic carbon [CI] and high-J CO(8-7) data as well as archival NIR/FIR, cm continuum, CO(3-2) and HI data are folded into the analysis. Results: While the overall [CII] morphology follows the general ionized gas, the channel maps show multiple bubble-like structures with sizes on the order of ~80-100" (~1.0-1.28pc). While at least one of them may be an individual feedback bubble driven by the main exciting sources of the region, the other bubble-morphologies may also be due to the intrinsically porous structure of the HII region. An analysis of the expansion velocities around 10km s^{-1} indicates that thermal expansion is not sufficient but that wind-driving from the central O-stars is required. The most blue-shifted [CII] component has barely any molecular or atomic counterparts. At the interface to the molecular cloud, we find a typical photon-dominated region (PDR) with a bar-shape. Ionized, atomic and molecular carbon show a layered structure in this PDR. The carbon in the PDR is dominated by its ionized form with atomic and molecular masses of ~0.45+-0.1M_{\odot} and ~1.2+-0.1M_{\odot}, respectively, compared to the ionized carbon in the range of 3.6-9.7M_{\odot}. Conclusions: The NGC7538 HII region exhibits a diverse set of sub-structures that interact with each other as well as with the adjacent cloud. Compared to other recent [CII] observations of HII regions (e.g., Orion Veil, RCW120, RCW49), bubble-shape morphologies revealed in [CII] emission, indicative of expanding shells, are recurring structures of PDRs.

preprint2022arXiv

Filamentary structures of ionized gas in Cygnus X

Ionized gas probes the influence of massive stars on their environment. The Cygnus X region (d~1.5 kpc) is one of the most massive star forming complexes in our Galaxy, in which the Cyg OB2 association (age of 3-5 Myr and stellar mass $2 \times 10^{4}$ M$_{\odot}$) has a dominant influence. We observe the Cygnus X region at 148 MHz using the Low Frequency Array (LOFAR) and take into account short-spacing information during image deconvolution. Together with data from the Canadian Galactic Plane Survey, we investigate the morphology, distribution, and physical conditions of low-density ionized gas in a $4^{\circ} \times 4^{\circ}$ (100 pc $\times$ 100 pc) region at a resolution of 2' (0.9 pc). The Galactic radio emission in the region analyzed is almost entirely thermal (free-free) at 148 MHz, with emission measures of $10^3 < EM~{\rm[pc~cm^{-6}]} < 10^6$. As filamentary structure is a prominent feature of the emission, we use DisPerSE and FilChap to identify filamentary ridges and characterize their radial ($EM$) profiles. The distribution of radial profiles has a characteristic width of 4.3 pc and a power-law distribution ($β= -1.8 \pm 0.1$) in peak $EM$ down to our completeness limit of 4200 pc cm$^{-6}$. The electron densities of the filamentary structure range from $10 < n_e~{\rm[cm^{-3}]} < 400$ with a median value of 35 cm$^{-3}$, remarkably similar to [N II] surveys of ionized gas. Cyg OB2 may ionize at most two-thirds of the total ionized gas and the ionized gas in filaments. More than half of the filamentary structures are likely photoevaporating surfaces flowing into a surrounding diffuse (~5 cm$^{-3}$) medium. However, this is likely not the case for all ionized gas ridges. A characteristic width in the distribution of ionized gas points to the stellar winds of Cyg OB2 creating a fraction of the ionized filaments through swept-up ionized gas or dissipated turbulence.

preprint2022arXiv

Polycyclic Aromatic Hydrocarbon emission model in photodissociation regions I: Application to the 3.3, 6.2, and 11.2 $μ$m bands

We present a charge distribution based model that computes the infrared spectrum of polycyclic aromatic hydrocarbon (PAH) molecules using recent measurements or quantum chemical calculations of specific PAHs. The model is applied to a sample of well-studied photodissociation regions (PDRs) with well-determined physical conditions (the radiation field strength, $G_{0}$, electron density $n_{e}$, and the gas temperature, $T_{\rm gas}$). Specifically, we modelled the emission of five PAHs ranging in size from 18 to 96 carbon atoms, over a range of physical conditions characterized by the ionization parameter $γ= G_{0}\times T_{\rm gas}^{1/2}/n_{e}$. The anions emerge as the dominant charge carriers in low $γ$ ($< 2\times 10^{2}$) environments, neutrals in the intermediate $γ$ ($10^{3} - 10^{4}$) environments, and cations in the high $γ$ ($ > 10^{5}$) environments. Furthermore, the PAH anions and cations exhibit similar spectral characteristics. The similarity in the cationic and anionic spectra translates into the interpretation of the 6.2/(11.0+11.2) band ratio, with high values of this ratio associated with large contributions from either cations or anions. The model's predicted values of 6.2/(11.0+11.2) and 3.3/6.2 compared well to the observations in the PDRs NGC 7023, NGC 2023, the horsehead nebula, the Orion bar, and the diffuse ISM, demonstrating that changes in the charge state can account for the variations in the observed PAH emission. We also reassess the diagnostic potential of the 6.2/(11.0+11.2) vs 3.3/(11.0+11.2) ratios and show that without any prior knowledge about $γ$, the 3.3/(11.0+11.2) can predict the PAH size, but the 6.2/(11.0+11.2) cannot predict the $γ$ of the astrophysical environment.

preprint2022arXiv

Self-absorption in [CII], $^{12}$CO, and HI in RCW120. Building up a geometrical and physical model of the region

Revealing the 3D dynamics of HII regions and their associated molecular clouds is important for understanding the longstanding problem as to how stellar feedback affects the density structure and kinematics of the interstellar medium. We employed observations of the HII region RCW 120 in [CII], observed within the SOFIA legacy program FEEDBACK, and the $^{12}$CO and $^{13}$CO (3$\to$2) lines, obtained with APEX. In addition we used HI data from the Southern Galactic Plane Survey. Two radiative transfer models were used to fit the observed data. A line profile analysis with the 1D non-LTE radiative transfer code SimLine proves that the CO emission cannot stem from a spherically symmetric molecular cloud configuration. With a two-layer multicomponent model, we then quantified the amount of warm background and cold foreground gas. There is a deficit of CO emission along the line-of-sight toward the center of the HII region which indicates that the HII region is associated with a flattened molecular cloud. Self-absorption in the CO line may hide signatures of infalling and expanding molecular gas. The [CII] emission arises from an expanding [CII] bubble and from the PDRs. A significant part of [CII] emission is absorbed in a cool (~60-100 K), low-density (<500 cm$^{-3}$) atomic foreground layer with a thickness of a few parsec. We propose that the RCW 120 HII region formed in a flattened molecular cloud and is now bursting out of its parental cloud. The compressed surrounding molecular layer formed a torus around the spherically expanding HII bubble. This scenario can possibly be generalized for other HII bubbles and would explain the observed "flat" structure of molecular clouds associated with HII bubbles. We suggest that the [CII] absorption observed in many star-forming regions is at least partly caused by low-density, cool, HI-envelopes surrounding the molecular clouds.

preprint2022arXiv

The Infrared Database of Extragalactic Observables from Spitzer. -- II. The Database & The Diagnostic Power of Crystalline Silicate Features in Galaxy Spectra

We present the Infrared Database of Extragalactic Observables from Spitzer (IDEOS), a homogeneous, publicly available, database of 77 fitted mid-infrared observables in the 5.4-36um range, comprising measurements for 3335 galaxies observed in the low-resolution staring mode of the Infrared Spectrometer onboard the Spitzer Space Telescope. Among the included observables are PAH fluxes and their equivalent widths, the strength of the 9.8um silicate feature, emission line fluxes, solid-state features, rest frame continuum fluxes, synthetic photometry, and a mid-infrared spectral classification. The IDEOS spectra were selected from the Cornell Atlas of Spitzer-IRS Sources. To our surprise we have detected at a >95% confidence level crystalline silicates in the spectra of 786 IDEOS galaxies. The detections range from single band detections to detections of all fitted crystalline bands (16, 19, 23, 28 and 33um). We find the strength of the crystalline silicate bands to correlate with the amorphous silicate strength, and the change from an emission to an absorption feature to occur at higher obscuration as the wavelength of the crystalline silicate band is longer. These observed characteristics are consistent with an origin for the amorphous and crystalline silicate features in a centrally heated dust geometry, either an edge-on disk or a cocoon. We find the 23 and 33um crystalline silicate bands to be well-suited to classify the obscuration level of galactic nuclei, even in the presence of strong circumnuclear star formation. Based on our detection statistics, we conclude that crystalline silicates are a common component of the interstellar medium of galactic nuclei.

preprint2022arXiv

The SOFIA FEEDBACK Legacy Survey: Dynamics and mass ejection in the bipolar HII region RCW 36

We present [CII] 158 $μ$m and [OI] 63 $μ$m observations of the bipolar HII region RCW 36 in the Vela C molecular cloud, obtained within the SOFIA legacy project FEEDBACK, which is complemented with APEX $^{12/13}$CO(3-2) and Chandra X-ray (0.5-7 keV) data. This shows that the molecular ring, forming the waist of the bipolar nebula, expands with a velocity of 1 - 1.9 km s$^{-1}$. We also observe an increased linewidth in the ring indicating that turbulence is driven by energy injection from the stellar feedback. The bipolar cavity hosts blue-shifted expanding [CII] shells at 5.2$\pm$0.5$\pm$0.5 km s$^{-1}$ (statistical and systematic uncertainty) which indicates that expansion out of the dense gas happens non-uniformly and that the observed bipolar phase might be relatively short ($\sim$0.2 Myr). The X-ray observations show diffuse emission that traces a hot plasma, created by stellar winds, in and around RCW 36. At least 50 \% of the stellar wind energy is missing in RCW 36. This is likely due to leakage which is clearing even larger cavities around the bipolar RCW 36 region. Lastly, the cavities host high-velocity wings in [CII] which indicates relatively high mass ejection rates ($\sim$5$\times$10$^{-4}$ M$_{\odot}$ yr$^{-1}$). This could be driven by stellar winds and/or radiation pressure, but remains difficult to constrain. This local mass ejection, which can remove all mass within 1 pc of RCW 36 in 1-2 Myr, and the large-scale clearing of ambient gas in the Vela C cloud indicates that stellar feedback plays a significant role in suppressing the star formation efficiency (SFE).

preprint2021arXiv

[CII] 158$μ$m emission from Orion A. II. Photodissociation region physics

The [CII] 158$μ$m fine-structure line is the dominant cooling line of moderate-density photodissociation regions (PDRs) illuminated by moderately bright far-ultraviolet (FUV) radiation fields. We aim to understand the origin of [CII] emission and its relation to other tracers of gas and dust in PDRs. One focus is a study of the heating efficiency of interstellar gas as traced by the [CII] line to test models of the photoelectric heating of neutral gas by polycyclic aromatic hydrocarbon (PAH) molecules and very small grains. We make use of a one-square-degree map of velocity-resolved [CII] line emission toward the Orion Nebula complex, and split this out into the individual spatial components, the expanding Veil Shell, the surface of OMC4, and the PDRs associated with the compact HII region of M43 and the reflection nebula NGC 1977. We employed Herschel far-infrared photometric images to determine dust properties. Moreover, we compared with Spitzer mid-infrared photometry to trace hot dust and large molecules, and velocity-resolved IRAM 30m CO(2-1) observations of the molecular gas. The [CII] intensity is tightly correlated with PAH emission in the IRAC 8$μ$m band and far-infrared emission from warm dust. The correlation between [CII] and CO(2-1) is very different in the four subregions and is very sensitive to the detailed geometry. Constant-density PDR models are able to reproduce the observed [CII], CO(2-1), and integrated far-infrared (FIR) intensities. We observe strong variations in the photoelectric heating efficiency in the Veil Shell behind the Orion Bar and these variations are seemingly not related to the spectral properties of the PAHs. The [CII] emission from the Orion Nebula complex stems mainly from moderately illuminated PDR surfaces. Future observations with the James Webb Space Telescope can shine light on the PAH properties that may be linked to these variations.

preprint2021arXiv

IR photofragmentation of the Phenyl Cation: Spectroscopy and Fragmentation Pathways

We present the gas-phase infrared spectra of the phenyl cation, phenylium, in its perprotio C$_6$H$_5^+$ and perdeutero C$_6$D$_5^+$ forms, in the 260-1925 cm$^{-1}$ (5.2-38 $μ$m) spectral range, and investigate the observed photofragmentation. The spectral and fragmentation data were obtained using Infrared Multiple Photon Dissociation (IRMPD) spectroscopy within a Fourier Transform Ion Cyclotron Resonance Mass Spectrometer (FTICR MS) located inside the cavity of the free electron laser FELICE (Free Electron Laser for Intra-Cavity Experiments). The $^1$A$_1$ singlet nature of the phenylium ion is ascertained by comparison of the observed IR spectrum with DFT calculations, using both harmonic and anharmonic frequency calculations. To investigate the observed loss of predominantly [2C,nH] (n=2-4) fragments, we explored the potential energy surface (PES) to unravel possible isomerization and fragmentation reaction pathways. The lowest energy pathways toward fragmentation include direct H elimination, and a combination of facile ring-opening mechanisms ($\leq2.4$ eV), followed by elimination of H or CCH$_2$. Energetically, all H-loss channels found are more easily accessible than CCH$_2$-loss. Calculations of the vibrational density of states for the various intermediates show that at high internal energies, ring opening is the thermodynamically the most advantageous, eliminating direct H-loss as a competing process. The observed loss of primarily [2C,2H] can be explained through entropy calculations that show favored loss of [2C,2H] at higher internal energies.

preprint2020arXiv

Breakdown products of gaseous polycyclic aromatic hydrocarbons investigated with infrared ion spectroscopy

We report on a common fragment ion formed during the electron-ionization-induced fragmentation of three different three-ring polycyclic aromatic hydrocarbons (PAHs), fluorene (C$_{13}$H$_{10}$), 9,10-dihydrophenanthrene (C$_{14}$H$_{12}$), and 9,10-dihydroanthracene (C$_{14}$H$_{12}$). The infrared spectra of the mass-isolated product ions with $m/z=165$ were obtained in a Fourier transform ion cyclotron resonance mass spectrometer whose cell was placed inside the optical cavity of an infrared free-electron laser, thus providing the high photon fluence required for efficient infrared multiple-photon dissociation. The infrared spectra of the $m/z=165$ species generated from the three different precursors were found to be similar, suggesting the formation of a single C$_{13}$H$_{9}^+$ isomer. Theoretical calculations using density functional theory (DFT) revealed the fragment's identity as the closed-shell fluorenyl cation. Decomposition pathways from each parent precursor to the fluorenyl ion are proposed on the basis of DFT calculations. The identification of a single fragmentation product from three different PAHs supports the notion of the existence of common decomposition pathways of PAHs in general and can aid in understanding the fragmentation chemistry of astronomical PAH species.

preprint2020arXiv

Expanding bubbles in Orion A: [CII] observations of M42, M43, and NGC 1977

The Orion Molecular Cloud is the nearest massive-star forming region. Massive stars have profound effects on their environment due to their strong radiation fields and stellar winds. Velocity-resolved observations of the [CII] $158\,μ\mathrm{m}$ fine-structure line allow us to study the kinematics of UV-illuminated gas. Here, we present a square-degree-sized map of [CII] emission from the Orion Nebula complex obtained by the upGREAT instrument onboard SOFIA, covering the entire Orion Nebula (M42) plus M43 and the nebulae NGC 1973, 1975, and 1977. We compare the stellar characteristics of these three regions with the kinematics of the expanding bubbles surrounding them. The bubble blown by the O7V star $θ^1$ Ori C in the Orion Nebula expands rapidly, at $13\,\mathrm{km\,s^{-1}}$. Simple analytical models reproduce the characteristics of the hot interior gas and the neutral shell of this wind-blown bubble and give us an estimate of the expansion time of $0.2\,\mathrm{Myr}$. M43 with the B0.5V star NU Ori also exhibits an expanding bubble structure, with an expansion velocity of $6\,\mathrm{km\,s^{-1}}$. Comparison with analytical models for the pressure-driven expansion of H\,{\sc ii} regions gives an age estimate of $0.02\,\mathrm{Myr}$. The bubble surrounding NGC 1973, 1975, and 1977 with the central B1V star 42 Orionis expands at $1.5\,\mathrm{km\,s^{-1}}$, likely due to the over-pressurized ionized gas as in the case of M43. We derive an age of $0.4\,\mathrm{Myr}$ for this structure. We conclude that the bubble of the Orion Nebula is driven by the mechanical energy input by the strong stellar wind from $θ^1$ Ori C, while the bubbles associated with M43 and NGC 1977 are caused by the thermal expansion of the gas ionized by their central later-type massive stars.

preprint2020arXiv

FEEDBACK: a SOFIA Legacy Program to Study Stellar Feedback in Regions of Massive Star Formation

FEEDBACK is a SOFIA legacy program dedicated to study the interaction of massive stars with their environment. It performs a survey of 11 galactic high mass star forming regions in the 158 $μ$m (1.9 THz) line of CII and the 63 $μ$m (4.7 THz) line of OI. We employ the 14 pixel LFA and 7 pixel HFA upGREAT instrument to spectrally resolve (0.24 MHz) these FIR structure lines. With an observing time of 96h, we will cover $\sim$6700 arcmin$^2$ at 14.1$''$ angular resolution for the CII line and 6.3$''$ for the OI line. The observations started in spring 2019 (Cycle 7). Our aim is to understand the dynamics in regions dominated by different feedback processes from massive stars such as stellar winds, thermal expansion, and radiation pressure, and to quantify the mechanical energy injection and radiative heating efficiency. The CII line provides the kinematics of the gas and is one of the dominant cooling lines of gas for low to moderate densities and UV fields. The OI line traces warm and high-density gas, excited in photodissociations regions with a strong UV field or by shocks. The source sample spans a broad range in stellar characteristics from single OB stars, to small groups of O stars, to rich young stellar clusters, to ministarburst complexes. It contains well-known targets such as Aquila, the Cygnus X region, M16, M17, NGC7538, NGC6334, Vela, and W43 as well as a selection of HII region bubbles, namely RCW49, RCW79, and RCW120. These CII maps, together with the less explored OI 63 $μ$m line, provide an outstanding database for the community. They will be made publically available and will trigger further studies and follow-up observations.

preprint2020arXiv

Molecular globules in the Veil bubble of Orion. IRAM 30m 12CO, 13CO, and C18O 2-1 expanded maps of Orion A

Strong winds and ultraviolet (UV) radiation from O-type stars disrupt and ionize their molecular core birthplaces, sweeping up material into parsec-size shells. Owing to dissociation by starlight, the thinnest shells are expected to host low molecular abundances and therefore little star formation. Here, we expand previous maps taken with the IRAM 30m telescope and present square-degree 12CO and 13CO (J=2-1) maps of the wind-driven "Veil bubble'' that surrounds the Trapezium cluster and its natal Orion molecular core (OMC). Although widespread and extended CO emission is largely absent from the Veil, we show that several CO "globules'' exist and are embedded in the [CII]158um-bright shell that confines the bubble. This includes the first detection of quiescent CO at negative LSR velocities in Orion. Given the harsh UV irradiation conditions in this translucent material, the detection of CO globules is surprising. These globules are small (R=7,100 AU), not massive (M=0.3M_Sun), and are moderately dense: n_ H=4x10^4 cm^-3 (median values). They are confined by the external pressure of the shell, P_ext/k~10^7 cm^-3 K, and are likely magnetically supported. They are either transient objects formed by instabilities or have detached from pre-existing molecular structures, sculpted by the passing shock associated with the expanding shell and by UV radiation from the Trapezium. Some represent the first stages in the formation of small pillars, others of isolated small globules. Although their masses do not suggest they will form stars, one globule matches the position of a known YSO. The lack of extended CO in the "Veil shell'' demonstrates that feedback from massive stars expels, agitates, and reprocesses most of the disrupted molecular cloud gas, thereby limiting the star-formation rate in the region. The presence of globules is a result of this feedback.

preprint2020arXiv

The H$_2$O Spectrum of the Massive Protostar AFGL 2136 IRS 1 from 2 to 13 $μ$m at High Resolution: Probing the Circumstellar Disk

We have observed the massive protostar AFGL 2136 IRS 1 in multiple wavelength windows in the near-to-mid-infrared at high ($\sim3$ km s$^{-1}$) spectral resolution using VLT+CRIRES, SOFIA+EXES, and Gemini North+TEXES. There is an abundance of H$_2$O absorption lines from the $ν_1$ and $ν_3$ vibrational bands at 2.7 $μ$m, from the $ν_2$ vibrational band at 6.1 $μ$m, and from pure rotational transitions near 10-13 $μ$m. Analysis of state-specific column densities derived from the resolved absorption features reveals that an isothermal absorbing slab model is incapable of explaining the relative depths of different absorption features. In particular, the strongest absorption features are much weaker than expected, indicating optical depth effects resulting from the absorbing gas being well-mixed with the warm dust that serves as the "background" continuum source at all observed wavelengths. The velocity at which the strongest H$_2$O absorption occurs coincides with the velocity centroid along the minor axis of the compact disk in Keplerian rotation recently observed in H$_2$O emission with ALMA. We postulate that the warm regions of this dust disk dominate the continuum emission at near-to-mid infrared wavelengths, and that H$_2$O and several other molecules observed in absorption are probing this disk. Absorption line profiles are not symmetric, possibly indicating that the warm dust in the disk that produces the infrared continuum has a non-uniform distribution similar to the substructure observed in 1.3 mm continuum emission.

preprint2019arXiv

Searching for the largest bound atoms in space

(abridged) Radio recombination lines (RRLs) at frequencies $ν$ < 250 MHz trace the cold, diffuse phase of the ISM. Next generation low frequency interferometers, such as LOFAR, MWA and the future SKA, with unprecedented sensitivity, resolution, and large fractional bandwidths, are enabling the exploration of the extragalactic RRL universe. We observed the radio quasar 3C 190 (z~1.2) with the LOFAR HBA. In reducing this data for spectroscopic analysis, we have placed special emphasis on bandpass calibration. We devised cross-correlation techniques to significantly identify the presence of RRLs in a low frequency spectrum. We demonstrate the utility of this method by applying it to existing low-frequency spectra of Cassiopeia A and M 82, and to the new observations of 3C 190. RRLs have been detected in the foreground of 3C 190 at z = 1.12355 (assuming a carbon origin), owing to the first detection of RRLs outside of the local universe (first reported in Emig et al. 2019). Towards the Galactic supernova remnant Cas A, we uncover three new detections: (1) C$ε$-transitions ($Δ$n = 5) for the first time at low radio frequencies, (2) H$α$-transitions at 64 MHz with a FWHM of 3.1 km/s, the most narrow and one of the lowest frequency detections of hydrogen to date, and (3) C$α$ at v$_{LSR}$ = 0 km/s in the frequency range 55-78 MHz for the first time. Additionally we recover C$α$, C$β$, C$γ$, and C$δ$ from the -47 km/s and -38 km/s components. In the nearby starburst galaxy, M 82, we do not find a significant feature. Our current searches for RRLs in LOFAR observations are limited to narrow (< 100 km/s) features, owing to the relatively small number of channels available for continuum estimation. Future strategies making use of larger contiguous frequency coverage would aid calibration to deeper sensitivities and broader features.

preprint2015arXiv

A bimodal dust grain distribution in the IC 434 HII region

Recent studies of dust in the interstellar medium have challenged the capabilities and validity of current dust models, indicating that the properties of dust evolve as it transits between different phases of the interstellar medium. We conduct a multi-wavelength study of the dust emission from the ionized gas of the IC 434 emission nebula, and combine this with modeling, from large scales that provide insight into the history of the IC 434/L1630 region, to small scales that allow us to infer quantitative properties of the dust content inside the H II region. The dust enters the H II region through momentum transfer with a champagne flow of ionized gas, set up by a chance encounter between the L1630 molecular cloud and the star cluster of $σ$ Ori. We observe two clearly separated dust populations inside the ionized gas, that show different observational properties, as well as contrasting optical properties. Population A is colder ($\sim$ 25 K) than predicted by widely-used dust models, its temperature is insensitive to an increase of the impinging radiation field, is momentum-coupled to the gas, and efficiently absorbs radiation pressure to form a dust wave at 1.0 pc ahead of $σ$ Ori AB. Population B is characterized by a constant [20/30] flux ratio throughout the HII region, heats up to $\sim$ 75 K close to the star, and is less efficient in absorbing radiation pressure, forming a dust wave at 0.1 pc from the star. We conclude that the dust inside IC 434 is bimodal. The characteristics of population A are remarkable and can not be explained by current dust models. Population B are grains that match the classical description of spherical, compact dust. Our results confirm recent work that stress the importance of variations in the dust properties between different regions of the interstellar medium.

preprint2015arXiv

The Physics of the Cold Neutral Medium: Low-frequency Carbon Radio Recombination Lines with the Square Kilometre Array

The Square Kilometre Array (SKA) will transform our understanding of the role of the cold, atomic gas in galaxy evolution. The interstellar medium (ISM) is the repository of stellar ejecta and the birthsite of new stars and, hence, a key factor in the evolution of galaxies over cosmic time. Cold, diffuse, atomic clouds are a key component of the ISM, but so far this phase has been difficult to study, because its main tracer, the HI 21 cm line, does not constrain the basic physical information of the gas (e.g., temperature, density) well. The SKA opens up the opportunity to study this component of the ISM through a complementary tracer in the form of low-frequency (<350 MHz) carbon radio recombination lines (CRRL). These CRRLs provide a sensitive probe of the physical conditions in cold, diffuse clouds. The superb sensitivity, large field of view, frequency resolution and coverage of the SKA allows for efficient surveys of the sky, that will revolutionize the field of low-frequency recombination line studies. By observing these lines with the SKA we will be able determine the thermal balance, chemical enrichment, and ionization rate of the cold, atomic medium from degree-scales down to scales corresponding to individual clouds and filaments in our Galaxy, the Magellanic Clouds and beyond. Furthermore, being sensitive only to the cold, atomic gas, observations of low-frequency CRRLs with the SKA will aid in disentangling the warm and cold constituents of the HI 21 cm emission.

preprint2014arXiv

Depletion of chlorine into HCl ice in a protostellar core

The freezeout of gas-phase species onto cold dust grains can drastically alter the chemistry and the heating-cooling balance of protostellar material. In contrast to well-known species such as carbon monoxide (CO), the freezeout of various carriers of elements with abundances $<10^{-5}$ has not yet been well studied. Our aim here is to study the depletion of chlorine in the protostellar core, OMC-2 FIR 4. We observed transitions of HCl and H2Cl+ towards OMC-2 FIR 4 using the Herschel Space Observatory and Caltech Submillimeter Observatory facilities. Our analysis makes use of state of the art chlorine gas-grain chemical models and newly calculated HCl-H$_{2}$ hyperfine collisional excitation rate coefficients. A narrow emission component in the HCl lines traces the extended envelope, and a broad one traces a more compact central region. The gas-phase HCl abundance in FIR 4 is 9e-11, a factor of only 0.001 that of volatile elemental chlorine. The H2Cl+ lines are detected in absorption and trace a tenuous foreground cloud, where we find no depletion of volatile chlorine. Gas-phase HCl is the tip of the chlorine iceberg in protostellar cores. Using a gas-grain chemical model, we show that the hydrogenation of atomic chlorine on grain surfaces in the dark cloud stage sequesters at least 90% of the volatile chlorine into HCl ice, where it remains in the protostellar stage. About 10% of chlorine is in gaseous atomic form. Gas-phase HCl is a minor, but diagnostically key reservoir, with an abundance of <1e-10 in most of the protostellar core. We find the 35Cl/37Cl ratio in OMC-2 FIR 4 to be 3.2\pm0.1, consistent with the solar system value.

preprint2014arXiv

Discovery of Carbon Radio Recombination Lines in M82

Carbon radio recombination lines (RRLs) at low frequencies (<=500 MHz) trace the cold, diffuse phase of the interstellar medium, which is otherwise difficult to observe. We present the detection of carbon RRLs in absorption in M82 with LOFAR in the frequency range of 48-64 MHz. This is the first extragalactic detection of RRLs from a species other than hydrogen, and below 1 GHz. Since the carbon RRLs are not detected individually, we cross-correlated the observed spectrum with a template spectrum of carbon RRLs to determine a radial velocity of 219 +- 9 km/s . Using this radial velocity, we stack 22 carbon-α transitions from quantum levels n = 468-508 to achieve an 8.5 sigma detection. The absorption line profile exhibits a narrow feature with peak optical depth of 0.003 and FWHM of 31 km/s. Closer inspection suggests that the narrow feature is superimposed on a broad, shallow component. The total line profile appears to be correlated with the 21 cm H I line profile reconstructed from H I absorption in the direction of supernova remnants in the nucleus. The narrow width and centroid velocity of the feature suggests that it is associated with the nuclear starburst region. It is therefore likely that the carbon RRLs are associated with cold atomic gas in the direction of the nucleus of M82.

preprint2014arXiv

Dust composition and mass-loss return from the luminous blue variable R71 in the LMC

We present an analysis of mid-and far-infrared (IR) spectrum and spectral energy distribution (SED) of the LBV R71 in the LMC.This work aims to understand the overall contribution of high-mass LBVs to the total dust-mass budget of the interstellar medium (ISM) of the LMC and compare this with the contribution from low-mass asymptotic giant branch (AGB) stars. As a case study, we analyze the SED of R71. We compiled all the available photometric and spectroscopic observational fluxes from various telescopes for a wide wavelength range (0.36 -- 250\,$μ$m). We determined the dust composition from the spectroscopic data, and derived the ejected dust mass, dust mass-loss rate, and other dust shell properties by modeling the SED of R71. We noted nine spectral features in the dust shell of R71 by analyzing Spitzer spectroscopic data. Among these, we identified three new crystalline silicate features. We computed our model spectrum by using 3D radiative transfer code MCMax. Our model calculation shows that dust is dominated by amorphous silicates, with some crystalline silicates, metallic iron, and a very tiny amount of polycyclic aromatic hydrocarbon (PAH) molecules. The presence of both silicates and PAHs indicates that the dust has a mixed chemistry. We derived a dust mass of 0.01 M$_\odot$, from which we arrive at a total ejected mass of $\approx$ 5 M$_\odot$. This implies a time-averaged dust mass-loss rate of 2.5$\times$10$^{-6}$ M$_\odot$\,yr$^{-1}$ with an explosion about 4000 years ago. We assume that the other five confirmed dusty LBVs in the LMC loose mass at a similar rate, and estimate the total contribution to the mass budget of the LMC to be $\approx$ 10$^{-5}$ M$_\odot$\,yr$^{-1}$, which is comparable to the contribution by all the AGB stars in the LMC. Based on our analysis on R71, we speculate that LBVs as a class may be an important dust source in the ISM of the LMC.

preprint2014arXiv

Exact Bound-Bound Gaunt Factor Values for Quantum Levels up to n=2000

Comparison of observations of radio recombination lines in the interstellar medium with theoretical models can be used to constrain electron temperature and density of the gas. An important component of the models is spontaneous transition rates between bound levels. Calculating these rates relies on accurate bound-bound oscillator strengths, which can be cast in terms of the Gaunt factor. The Gaunt factor contains terminating hypergeometric functions that cannot be calculated with sufficient accuracy for high quantum levels ($n \gtrsim 50$) by standard machine-precision methods. Methods to overcome the accuracy problem have been developed, which include asymptotic expansions and recursion relations. These methods, used in astrophysical models to calculate oscillator strengths, can introduce errors, sometimes up to as much as $\sim 8$ per cent. Detections of radio recombination lines with the new Low Frequency Array (LOFAR) has prompted an examination of theoretical models of the interstellar medium. We revisit the calculation of the Gaunt factor, employing modern arbitrary-precision computational methods to tabulate the Gaunt factor for transitions up to quantum level $n=2000$, sufficient to model low frequency Carbon radio recombination lines. The calculations provide a relative error of $\sim3\times 10^{-4}$ when compared to more detailed calculations including relativistic corrections. Our values for the Gaunt factor are provided for download in a tabular format to be used for a wide range of applications.

preprint2014arXiv

Herschel Planetary Nebula Survey (HerPlaNS) - First Detection of OH+ in Planetary Nebulae

We report the first detections of OH$^+$ emission in planetary nebulae (PNe). As part of an imaging and spectroscopy survey of 11 PNe in the far-IR using the PACS and SPIRE instruments aboard the Herschel Space Observatory, we performed a line survey in these PNe over the entire spectral range between 51 and 672$μ$m to look for new detections. OH$^+$ rotational emission lines at 152.99, 290.20, 308.48, and 329.77$μ$m were detected in the spectra of three planetary nebulae: NGC 6445, NGC 6720, and NGC 6781. Excitation temperatures and column densities derived from these lines are in the range of 27 to 47 K and 2$\times$10$^{10}$ to 4 $\times$10$^{11}$ cm$^{-2}$, respectively. In PNe, the OH+ rotational line emission appears to be produced in the photodissociation region (PDR) in these objects. The emission of OH+ is observed only in PNe with hot central stars (T$_{eff}$ > 100000 K), suggesting that high-energy photons may play a role in the OH+ formation and its line excitation in these objects, as it seems to be the case for ultraluminous galaxies.

preprint2014arXiv

Location and sizes of forsterite grains in protoplanetary disks: interpretation from the Herschel DIGIT programme

The spectra of protoplanetary disks contain mid- and far- infrared emission features produced by forsterite dust grains. The spectral features contain information about the forsterite temperature, chemical composition and grain size. We aim to characterize how the 23 and 69 micron features can be used to constrain the physical locations of forsterite in disks. We check for consistency between two independent forsterite temperature measurements: the 23/69 feature strength ratio and the shape of the 69 micron band. We performed radiative transfer modeling to study the effect of disk properties to the forsterite spectral features. Temperature-dependent forsterite opacities were considered in self-consistent models to compute forsterite emission from protoplanetary disks. Modelling grids are presented to study the effects of grain size, disk gaps, radial mixing and optical depth to the forsterite features. Independent temperature estimates derived from the 23/69 feature strength ratio and the 69 micron band shape are most inconsistent for HD141569 and Oph IRS 48. A case study of the disk of HD141569 shows two solutions to fit the forsterite spectrum. A model with T ~ 40 K, iron-rich (~0-1 % Fe) and 1 micron forsterite grains, and a model with warmer (T ~ 100 K), iron-free, and larger (10 micron) grains. We find that for disks with low upper limits of the 69 micron feature (most notably in flat, self-shadowed disks), the forsterite must be hot, and thus close to the star. We find no correlation between disk gaps and the presence or absence of forsterite features. We argue that the 69 micron feature of the evolved transitional disks HD141569 and Oph IRS 48 is most likely a tracer of larger (i.e. ~10 micron) forsterite grains.

preprint2014arXiv

Polycyclic Aromatic Hydrocarbons with armchair edges and the 12.7 μm band

In this Letter we report the results of density functional theory calculations on medium-sized neutral Polycyclic Aromatic Hydrocarbon (PAH) molecules with armchair edges. These PAH molecules possess strong C-H stretching and bending modes around 3 μm and in the fingerprint region (10-15 μm), and also strong ring deformation modes around 12.7 μm. Perusal of the entries in the NASA Ames PAHs Database shows that ring deformation modes of PAHs are common - although generally weak. We then propose that armchair PAHs with NC >65 are responsible for the 12.7 μm Aromatic Infrared Band in HII regions and discuss astrophysical implications in the context of the PAH life-cycle.

preprint2014arXiv

Pore evolution in interstellar ice analogues: simulating the effects of temperature increase

Context. The level of porosity of interstellar ices - largely comprised of amorphous solid water (ASW) - contains clues on the trapping capacity of other volatile species and determines the surface accessibility that is needed for solid state reactions to take place. Aims. Our goal is to simulate the growth of amorphous water ice at low temperature (10 K) and to characterize the evolution of the porosity (and the specific surface area) as a function of temperature (from 10 to 120 K). Methods. Kinetic Monte Carlo simulations are used to mimic the formation and the thermal evolution of pores in amorphous water ice. We follow the accretion of gas-phase water molecules as well as their migration on surfaces with different grid sizes, both at the top growing layer and within the bulk. Results. We show that the porosity characteristics change substantially in water ice as the temperature increases. The total surface of the pores decreases strongly while the total volume decreases only slightly for higher temperatures. This will decrease the overall reaction efficiency, but in parallel, small pores connect and merge, allowing trapped molecules to meet and react within the pores network, providing a pathway to increase the reaction efficiency. We introduce pore coalescence as a new solid state process that may boost the solid state formation of new molecules in space and has not been considered so far.

preprint2014arXiv

Radiation-pressure-driven dust waves inside bursting interstellar bubbles

Massive stars drive the evolution of the interstellar medium through their radiative and mechanical energy input. After their birth, they form bubbles of hot gas surrounded by a dense shell. Traditionally, the formation of bubbles is explained through the input of a powerful stellar wind, even though direct evidence supporting this scenario is lacking. Here we explore the possibility that interstellar bubbles seen by the Spitzer- and Herschel space telescopes, blown by stars with log(L/L_sun) < 5.2, form and expand due to the thermal pressure accompanying ionization of the surrounding gas. We show that density gradients in the natal cloud or a puncture in the swept up shell lead to an ionized gas flow through the bubble into the general interstellar medium, which is traced by a dust wave near the star, demonstrating the importance of radiation pressure during this phase. Dust waves provide a natural explanation for the presence of dust inside H II bubbles, offer a novel method to study dust in H II regions and provide direct evidence that bubbles are relieving their pressure into the ISM through a champagne flow, acting as a probe of the radiative interaction of a massive star with its surroundings. We create a parameter space connecting the ambient density, the ionizing source luminosity, and the position of the dust wave, while using the well-studied H II bubbles RCW 120 and RCW 82 as benchmarks of our model. Finally, we briefly examine the implications of our study for the environments of super star clusters formed in UltraLuminous InfraRed Galaxies (ULIRGs), merging galaxies and the early Universe, which occur in very luminous and dense environments and where radiation pressure is expected to dominate the dynamical evolution.

preprint2014arXiv

The Mass-Loss Return from Evolved Stars to the Large Magellanic Cloud II: Dust Properties for Oxygen-Rich Asymptotic Giant Branch Stars

We model multi-wavelength broadband UBVIJHKs and Spitzer IRAC and MIPS photometry and IRS spectra from the SAGE and SAGE-Spec observing programs of two oxygen-rich asymptotic giant branch (O-rich AGB) stars in the Large Magellanic Cloud (LMC) using radiative transfer models of dust shells around stars. We chose a star from each of the bright and faint O-rich AGB populations found by earlier studies of the SAGE sample in order to derive a baseline set of dust properties to be used in the construction of an extensive grid of radiative transfer models of the O-rich AGB stars found in the SAGE surveys. From the bright O-rich AGB population we chose HV 5715, and from the faint O-rich AGB population we chose SSTISAGE1C J052206.92-715017.6 (SSTSAGE052206). We found the complex indices of refraction of oxygen-deficient silicates from Ossenkopf et al (1992) and a "KMH"-like grain size distribution with gamma of -3.5, a_min of 0.01 microns, and a_0 of 0.1 microns to be reasonable dust properties for these models. There is a slight indication that the dust around the faint O-rich AGB may be more silica-rich than that around the bright O-rich AGB. Simple models of gas emission suggest a relatively extended gas envelope for the faint O-rich AGB star modeled, consistent with the relatively large dust shell inner radius for the same model. The models compute the dust mass-loss rates for the two stars to be 2.0x10^-9 M_Sun/yr and 2.3x10^-9 M_Sun/yr, respectively. When a dust-to-gas mass ratio of 0.002 is assumed for SSTSAGE052206 and HV 5715, the dust mass-loss rates imply total mass-loss rates of 1.0x10^-6 M_Sun/yr and 1.2x10^-6 M_Sun/yr, respectively. These properties of the dust shells and stars, as inferred from our models of the two stars, are found to be consistent with properties observed or assumed by detailed studies of other O-rich AGB stars in the LMC and elsewhere.

preprint2012arXiv

The enigmatic nature of the circumstellar envelope and bow shock surrounding Betelgeuse as revealed by Herschel. I. Evidence of clumps, multiple arcs, and a linear bar-like structure

Context. The interaction between stellar winds and the interstellar medium (ISM) can create complex bow shocks. The photometers on board the Herschel Space Observatory are ideally suited to studying the morphologies of these bow shocks. Aims. We aim to study the circumstellar environment and wind-ISM interaction of the nearest red supergiant, Betelgeuse. Methods. Herschel PACS images at 70, 100, and 160 micron and SPIRE images at 250, 350, and 500 micron were obtained by scanning the region around Betelgeuse. These data were complemented with ultraviolet GALEX data, near-infrared WISE data, and radio 21 cm GALFA-HI data. The observational properties of the bow shock structure were deduced from the data and compared with hydrodynamical simulations. Results. The infrared Herschel images of the environment around Betelgeuse are spectacular, showing the occurrence of multiple arcs at 6-7 arcmin from the central target and the presence of a linear bar at 9 arcmin. Remarkably, no large-scale instabilities are seen in the outer arcs and linear bar. The dust temperature in the outer arcs varies between 40 and 140 K, with the linear bar having the same colour temperature as the arcs. The inner envelope shows clear evidence of a non-homogeneous clumpy structure (beyond 15 arcsec), probably related to the giant convection cells of the outer atmosphere. The non-homogeneous distribution of the material even persists until the collision with the ISM. A strong variation in brightness of the inner clumps at a radius of 2 arcmin suggests a drastic change in mean gas and dust density some 32 000 yr ago. Using hydrodynamical simulations, we try to explain the observed morphology of the bow shock around Betelgeuse. Conclusions: [abbreviated]

preprint2011arXiv

Fusion between ground- and space-based mid-infrared observations

Mid-infrared astronomy (operating at wavelengths ranging from 2 to 25 $μ$m) has progressed significantly in the last decades, thanks to the improvement of detector techniques and the growing diameter of telescopes. Space observatories benefit from the absence of atmospheric absorption, allowing to reach the very high sensitivities needed to perform 3D hyperspectral observations, but telescopes are limited in diameter ($< 1$ meter) and therefore provide observations at low angular resolution (typically a few seconds of arc). On the other hand, ground-based facilities suffer from strong atmospheric absorption but use large telescopes (above 8m diameter) to perform sub-arcsecond angular resolution imaging through selected windows in the mid-infrared range. In this Paper, we present a method based on Lee and Seung's Non-negative Matrix Factorization (NMF) to merge data from space and ground based mid-infrared (mid-IR) telescopes in order to combine the best sensitivity, spectral coverage and angular resolution. We prove the efficiency of this technique when applied to real mid-IR astronomical data. We suggest that this method can be applied to any combination of low and high spatial resolution positive hyperspectral datasets, as long as the spectral variety of the data allows decomposition into components using NMF.

preprint2011arXiv

Polycyclic aromatic hydrocarbon processing in a hot gas

Context: PAHs are thought to be a ubiquitous and important dust component of the interstellar medium. However, the effects of their immersion in a hot (post-shock) gas have never before been fully investigated. Aims: We study the effects of energetic ion and electron collisions on PAHs in the hot post-shock gas behind interstellar shock waves. Methods: We calculate the ion-PAH and electron-PAH nuclear and electronic interactions, above the carbon atom loss threshold, in H II regions and in the hot post-shock gas, for temperatures ranging from 10^3 to 10^8 K. Results: PAH destruction is dominated by He collisions at low temperatures (T < 3x10^4 K), and by electron collisions at higher temperatures. Smaller PAHs are destroyed faster for T < 10^6 K, but the destruction rates are roughly the same for all PAHs at higher temperatures. The PAH lifetime in a tenuous hot gas (n_H ~ 0.01 cm^-3, T ~ 10^7 K), typical of the coronal gas in galactic outflows, is found to be about thousand years, orders of magnitude shorter than the typical lifetime of such objects. Conclusions: In a hot gas, PAHs are principally destroyed by electron collisions and not by the absorption of X-ray photons from the hot gas. The resulting erosion of PAHs occurs via C_2 loss from the periphery of the molecule, thus preserving the aromatic structure. The observation of PAH emission from a million degree, or more, gas is only possible if the emitting PAHs are ablated from dense, entrained clumps that have not yet been exposed to the full effect of the hot gas.

preprint2011arXiv

The Spitzer Spectroscopic Survey of the Small Magellanic Cloud (S4MC): Probing the Physical State of Polycyclic Aromatic Hydrocarbons in a Low-Metallicity Environment

We present results of mid-infrared spectroscopic mapping observations of six star-forming regions in the Small Magellanic Cloud from the Spitzer Spectroscopic Survey of the SMC (S4MC). We detect the mid-IR emission from polycyclic aromatic hydrocarbons (PAHs) in all of the mapped regions, greatly increasing the range of environments where PAHs have been spectroscopically detected in the SMC. We investigate the variations of the mid-IR bands in each region and compare our results to studies of the PAH bands in the SINGS sample and in a sample of low-metallicity starburst galaxies. PAH emission in the SMC is characterized by low ratios of the 6-9 micron features relative to the 11.3 micron feature and weak 8.6 and 17.0 micron features. Interpreting these band ratios in the light of laboratory and theoretical studies, we find that PAHs in the SMC tend to be smaller and less ionized than those in higher metallicity galaxies. Based on studies of PAH destruction, we argue that a size distribution shifted towards smaller PAHs cannot be the result of processing in the interstellar medium, but instead reflects differences in the formation of PAHs at low metallicity. Finally, we discuss the implications of our observations for our understanding of the PAH life-cycle in low-metallicity galaxies---namely that the observed deficit of PAHs may be a consequence of PAHs forming with smaller average sizes and therefore being more susceptible to destruction under typical interstellar medium conditions.

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

Cold dust clumps in dynamically hot gas

We present clumps of dust emission from Herschel observations of the Large Magellanic Cloud (LMC) and their physical and statistical properties. We catalog cloud features seen in the dust emission from Herschel observations of the LMC, the Magellanic type irregular galaxy closest to the Milky Way, and compare these features with HI catalogs from the ATCA+Parkes HI survey. Using an automated cloud-finding algorithm, we identify clouds and clumps of dust emission and examine the cumulative mass distribution of the detected dust clouds. The mass of cold dust is determined from physical parameters that we derive by performing spectral energy distribution fits to 250, 350, and 500 micronm emission from SPIRE observations using DUSTY and GRASIL radiative transfer calculation with dust grain size distributions for graphite/silicate in low-metallicity extragalactic environments. The dust cloud mass spectrum follows a power law distribution with an exponent of gamma=-1.8 for clumps larger than 400 solar mass and is similar to the HI mass distribution. This is expected from the theory of ISM structure in the vicinity of star formation.

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

Determining dust temperatures and masses in the Herschel era: the importance of observations longward of 200 micron

The properties of the dust grains (e.g., temperature and mass) can be derived from fitting far-IR SEDs (>100 micron). Only with SPIRE on Herschel has it been possible to get high spatial resolution at 200 to 500 micron that is beyond the peak (~160 micron) of dust emission in most galaxies. We investigate the differences in the fitted dust temperatures and masses determined using only <200 micron data and then also including >200 micron data (new SPIRE observations) to determine how important having >200 micron data is for deriving these dust properties. We fit the 100 to 350 micron observations of the Large Magellanic Cloud (LMC) point-by-point with a model that consists of a single temperature and fixed emissivity law. The data used are existing observations at 100 and 160 micron (from IRAS and Spitzer) and new SPIRE observations of 1/4 of the LMC observed for the HERITAGE Key Project as part of the Herschel Science Demonstration phase. The dust temperatures and masses computed using only 100 and 160 micron data can differ by up to 10% and 36%, respectively, from those that also include the SPIRE 250 & 350 micron data. We find that an emissivity law proportional to lambda^-1.5 minimizes the 100-350 micron fractional residuals. We find that the emission at 500 micron is ~10% higher than expected from extrapolating the fits made at shorter wavelengths. We find the fractional 500 micron excess is weakly anti-correlated with MIPS 24 micron flux and the total gas surface density. This argues against a flux calibration error as the origin of the 500 micron excess. Our results do not allow us to distinguish between a systematic variation in the wavelength dependent emissivity law or a population of very cold dust only detectable at lambda > 500 micron for the origin of the 500 micron excess.

preprint2010arXiv

Excitation of the molecular gas in the nuclear region of M82

We present high resolution HIFI spectroscopy of the nucleus of the archetypical starburst galaxy M82. Six 12CO lines, 2 13CO lines and 4 fine-structure lines are detected. Besides showing the effects of the overall velocity structure of the nuclear region, the line profiles also indicate the presence of multiple components with different optical depths, temperatures and densities in the observing beam. The data have been interpreted using a grid of PDR models. It is found that the majority of the molecular gas is in low density (n=10^3.5 cm^-3) clouds, with column densities of N_H=10^21.5 cm^-2 and a relatively low UV radiation field (GO = 10^2). The remaining gas is predominantly found in clouds with higher densities (n=10^5 cm^-3) and radiation fields (GO = 10^2.75), but somewhat lower column densities (N_H=10^21.2 cm^-2). The highest J CO lines are dominated by a small (1% relative surface filling) component, with an even higher density (n=10^6 cm^-3) and UV field (GO = 10^3.25). These results show the strength of multi-component modeling for the interpretation of the integrated properties of galaxies.

preprint2010arXiv

HD 95881: A gas rich to gas poor transition disk?

Context. Based on the far infrared excess the Herbig class of stars is divided into a group with flaring circumstellar disks (group I) and a group with flat circumstellar disks (group II). Dust sedimentation is generally proposed as an evolution mechanism to transform flaring disks into flat disks. Theory predicts that during this process the disks preserve their gas content, however observations of group II Herbig Ae stars demonstrate a lack of gas. Aims. We map the spatial distribution of the gas and dust around the group II Herbig Ae star HD 95881. Methods. We analyze optical photometry, Q-band imaging, infrared spectroscopy, and K and N-band interferometric spectroscopy. We use a Monte Carlo radiative transfer code to create a model for the density and temperature structure which quite accurately reproduces all the observables. Results. We derive a consistent picture in which the disk consists of a thick puffed up inner rim and an outer region which has a flaring gas surface and is relatively void of 'visible' dust grains. Conclusions. HD 95881 is in a transition phase from a gas rich flaring disk to a gas poor self-shadowed disk.

preprint2010arXiv

Herschel Inventory of The Agents of Galaxy Evolution (HERITAGE): the Large Magellanic Cloud dust

The HERschel Inventory of The Agents of Galaxy Evolution (HERITAGE) of the Magellanic Clouds will use dust emission to investigate the life cycle of matter in both the Large and Small Magellanic Clouds (LMC and SMC). Using the Herschel Space Observatory's PACS and SPIRE photometry cameras, we imaged a 2x8 square degree strip through the LMC, at a position angle of ~22.5 degrees as part of the science demonstration phase of the Herschel mission. We present the data in all 5 Herschel bands: PACS 100 and 160 μm and SPIRE 250, 350 and 500 μm. We present two dust models that both adequately fit the spectral energy distribution for the entire strip and both reveal that the SPIRE 500 μm emission is in excess of the models by 6 to 17%. The SPIRE emission follows the distribution of the dust mass, which is derived from the model. The PAH-to-dust mass (f_PAH) image of the strip reveals a possible enhancement in the LMC bar in agreement with previous work. We compare the gas mass distribution derived from the HI 21 cm and CO J=1-0 line emission maps to the dust mass map from the models and derive gas-to-dust mass ratios (GDRs). The dust model, which uses the standard graphite and silicate optical properties for Galactic dust, has a very low GDR = 65(+15,-18) making it an unrealistic dust model for the LMC. Our second dust model, which uses amorphous carbon instead of graphite, has a flatter emissivity index in the submillimeter and results in a GDR = 287(+25,-42) that is more consistent with a GDR inferred from extinction.

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

Polycyclic aromatic hydrocarbon processing in interstellar shocks

Context: PAHs appear to be an ubiquitous interstellar dust component but the effects of shocks waves upon them have never been fully investigated. Aims: To study the effects of energetic (~0.01-1 keV) ion (H, He and C) and electron collisions on PAHs in interstellar shock waves.Methods: We calculate the ion-PAH and electron-PAH nuclear and electronic interactions, above the threshold for carbon atom loss from a PAH, in 50-200 km/s shock waves in the warm intercloud medium. Results: Interstellar PAHs (Nc = 50) do not survive in shocks with velocities greater than 100 km/s and larger PAHs (Nc = 200) are destroyed for shocks with velocities greater/equal to 125 km/s. For shocks in the ~75 - 100 km/s range, where destruction is not complete, the PAH structure is likely to be severely denatured by the loss of an important fraction (20-40%) of the carbon atoms. We derive typical PAH lifetimes of the order of a few x10^8 yr for the Galaxy. These results are robust and independent of the uncertainties in some key parameters that have yet to be well-determined experimentally. Conclusions: The observation of PAH emission in shock regions implies that that emission either arises outside the shocked region or that those regions entrain denser clumps that, unless they are completely ablated and eroded in the shocked gas, allow dust and PAHs to survive in extreme environments.

preprint2010arXiv

SPITZER-IRS spectral fitting of discs around binary post-AGB stars - Corrigendum

Recently, we have discovered an error in our Monte-Carlo spectral fitting routine, more specifically where the errors on the fluxes were rescaled to get a reduced chi2 of 1. The rescaled errors were too big, resulting in too wide a range of good fits in our 100 step Monte-Carlo routine. This problem affects Figs. 7-9 and Tables A.1, A.2 in Gielen et al. (2008), Table 3 in Gielen et al. (2009a), and Table 4 in Gielen et al. (2009b). We corrected for this error and present the new values and errors in the tables below. The new values and errors nearly all fall within the old error range. Our best chi2 values and overall former scientific results are not affected. With these new errors some possible new trends in the dust parameters might be observed. These will be discussed in an upcoming paper where we extend the sample presented in Gielen et al. (2008) with newly obtained SPITZER-IRS data.

preprint2010arXiv

Spitzer's view on aromatic and aliphatic hydrocarbon emission in Herbig Ae stars

The chemistry of astronomical hydrocarbons, responsible for the well-known infrared emission features detected in a wide variety of targets, remains enigmatic. Here we focus on the group of young intermediate-mass Herbig Ae stars. We have analyzed the aliphatic and polycyclic aromatic hydrocarbon (PAH) emission features in the infrared spectra of a sample of 53 Herbig Ae stars, obtained with the Infrared Spectrograph aboard the Spitzer Space Telescope. We confirm that the PAH-to-stellar luminosity ratio is higher in targets with a flared dust disk. However, a few sources with a flattened dust disk still show relatively strong PAH emission. Since PAH molecules trace the gas disk, this indicates that gas disks may still be flared, while the dust disk has settled due to grain growth. There are indications that the strength of the 11.3 um feature also depends on dust disk structure, with flattened disks being less bright in this feature. We confirm that the CC bond features at 6.2 and 7.8 um shift to redder wavelengths with decreasing stellar effective temperature. Moreover, we show that this redshift is accompanied by a relative increase of aliphatic CH emission and a decrease of the aromatic 8.6 um CH feature strength. Cool stars in our sample are surrounded by hydrocarbons with a high aliphatic/aromatic CH ratio and a low aromatic CH/CC ratio, and vice versa for the hot stars. We conclude that, while the overall hydrocarbon emission strength depends on the dust disk's geometry, the relative differences seen in the IR emission features in disks around Herbig Ae stars are mainly due to chemical differences of the hydrocarbon molecules induced by the stellar UV field. Strong UV flux reduces the aliphatic component and emphasizes the spectral signature of the aromatic molecules in the IR spectra.

preprint2010arXiv

The SAGE-Spec Spitzer Legacy program: The life-cycle of dust and gas in the Large Magellanic Cloud

The SAGE-Spec Spitzer Legacy program is a spectroscopic follow-up to the SAGE-LMC photometric survey of the Large Magellanic Cloud carried out with the Spitzer Space Telescope. We present an overview of SAGE-Spec and some of its first results. The SAGE-Spec program aims to study the life cycle of gas and dust in the Large Magellanic Cloud, and to provide information essential to the classification of the point sources observed in the earlier SAGE-LMC photometric survey. We acquired 224.6 hours of observations using the InfraRed Spectrograph and the SED mode of the Multiband Imaging Photometer for Spitzer. The SAGE-Spec data, along with archival Spitzer spectroscopy of objects in the Large Magellanic Cloud, are reduced and delivered to the community. We discuss the observing strategy, the specific data reduction pipelines applied and the dissemination of data products to the scientific community. Initial science results include the first detection of an extragalactic "21 um" feature towards an evolved star and elucidation of the nature of disks around RV Tauri stars in the Large Magellanic Cloud. Towards some young stars, ice features are observed in absorption. We also serendipitously observed a background quasar, at a redshift of z~0.14, which appears to be host-less.

preprint2009arXiv

A Spitzer Space Telescope far-infrared spectral atlas of compact sources in the Magellanic Clouds. I. The Large Magellanic Cloud

[abridged] We present 52-93 micron spectra obtained with Spitzer in the MIPS-SED mode, of a representative sample of luminous compact far-IR sources in the LMC. These include carbon stars, OH/IR AGB stars, post-AGB objects and PNe, RCrB-type star HV2671, OH/IR red supergiants WOHG064 and IRAS05280-6910, B[e] stars IRAS04530-6916, R66 and R126, Wolf-Rayet star Brey3a, Luminous Blue Variable R71, supernova remnant N49, a large number of young stellar objects, compact HII regions and molecular cores, and a background galaxy (z~0.175). We use the spectra to constrain the presence and temperature of cold dust and the excitation conditions and shocks within the neutral and ionized gas, in the circumstellar environments and interfaces with the surrounding ISM. Evolved stars, including LBV R71, lack cold dust except in some cases where we argue that this is swept-up ISM. This leads to an estimate of the duration of the prolific dust-producing phase ("superwind") of several thousand years for both RSGs and massive AGB stars, with a similar fractional mass loss experienced despite the different masses. We tentatively detect line emission from neutral oxygen in the extreme RSG WOHG064, with implications for the wind driving. In N49, the shock between the supernova ejecta and ISM is revealed by its strong [OI] 63-micron emission and possibly water vapour; we estimate that 0.2 Msun of ISM dust was swept up. Some of the compact HII regions display pronounced [OIII] 88-micron emission. The efficiency of photo-electric heating in the interfaces of ionized gas and molecular clouds is estimated at 0.1-0.3%. We confirm earlier indications of a low nitrogen content in the LMC. Evidence for solid state emission features is found in both young and evolved object; some of the YSOs are found to contain crystalline water ice.

preprint2009arXiv

Chemical depletion in the Large Magellanic Cloud: RV Tauri stars and the photospheric feedback from their dusty discs

Aims: By studying the photospheric abundances of 4 RV Tauri stars in the LMC, we test whether the depletion pattern of refractory elements, seen in similar Galactic sources, is also common for extragalactic sources. Since this depletion process probably only occurs through interaction with a stable disc, we investigate the circumstellar environment of these sources. Methods: A detailed photospheric abundance study was performed using high-resolution UVES optical spectra. To study the circumstellar environment we use photometric data to construct the spectral energy distributions of the stars, and determine the geometry of the circumstellar environment, whereas low-resolution Spitzer-IRS infrared spectra are used to trace its mineralogy. Results: Our results show that, also in the LMC, the photospheres of RV Tauri stars are commonly affected by the depletion process, although it can differ significantly in strength from source to source. From our detailed disc modelling and mineralogy study, we find that this process, as in the Galaxy, appears closely related to the presence of a stable Keplerian disc. The newly studied extragalactic objects have similar observational characteristics as Galactic post-AGB binaries surrounded by a dusty disc, and are therefore also believed to be part of a binary system. One source shows a very small infrared excess, atypical for a disc source, but still has evidence for depletion. We speculate this could point to the presence of a very evolved disc, similar to debris discs seen around young stellar objects.

preprint2009arXiv

Dust coagulation and fragmentation in molecular clouds. I. How collisions between dust aggregates alter the dust size distribution

In dense molecular clouds collisions between dust grains alter the ISM-dust size distribution. We study this process by inserting the results from detailed numerical simulations of two colliding dust aggregates into a coagulation model that computes the dust size distribution with time. All collisional outcomes -- sticking, fragmentation (shattering, breakage, and erosion) -- are included and the effects on the internal structure of the aggregates are also tabulated. The dust aggregate evolution model is applied to an homogeneous and static cloud of temperature 10 K and gas densities between 10^3 and 10^7 cm^-3. The coagulation is followed locally on timescales of ~10^7 yr. We find that the growth can be divided into two stages: a growth dominated phase and a fragmentation dominated phase. Initially, the mass distribution is relatively narrow and shifts to larger sizes with time. At a certain point, dependent on the material properties of the grains as well as on the gas density, collision velocities will become sufficiently energetic to fragment particles, halting the growth and replenishing particles of lower mass. Eventually, a steady state is reached, where the mass distribution is characterized by a mass spectrum of approximately equal amount of mass per logarithmic size bin. The amount of growth that is achieved depends on the cloud's lifetime. If clouds exist on free-fall timescales the effects of coagulation on the dust size distribution are very minor. On the other hand, if clouds have long-term support mechanisms, the impact of coagulation is important, resulting in a significant decrease of the opacity on timescales longer than the initial collision timescale between big grains.

preprint2009arXiv

Ice chemistry in embedded young stellar objects in the Large Magellanic Cloud

We present spectroscopic observations of a sample of 15 embedded young stellar objects (YSOs) in the Large Magellanic Cloud (LMC). These observations were obtained with the Spitzer Infrared Spectrograph (IRS) as part of the SAGE-Spec Legacy program. We analyze the two prominent ice bands in the IRS spectral range: the bending mode of CO_2 ice at 15.2 micron and the ice band between 5 and 7 micron that includes contributions from the bending mode of water ice at 6 micron amongst other ice species. The 5-7 micron band is difficult to identify in our LMC sample due to the conspicuous presence of PAH emission superimposed onto the ice spectra. We identify water ice in the spectra of two sources; the spectrum of one of those sources also exhibits the 6.8 micron ice feature attributed to ammonium and methanol. We model the CO_2 band in detail, using the combination of laboratory ice profiles available in the literature. We find that a significant fraction (> 50%) of CO_2 ice is locked in a water-rich component, consistent with what is observed for Galactic sources. The majority of the sources in the LMC also require a pure-CO_2 contribution to the ice profile, evidence of thermal processing. There is a suggestion that CO_2 production might be enhanced in the LMC, but the size of the available sample precludes firmer conclusions. We place our results in the context of the star formation environment in the LMC.

preprint2009arXiv

Microscopic simulation of methanol and formaldehyde ice formation in cold dense cores

Methanol and its precursor formaldehyde are among the most studied organic molecules in the interstellar medium and are abundant in the gaseous and solid phases. We recently developed a model to simulate CO hydrogenation via H atoms on interstellar ice surfaces, the most important interstellar route to H2CO and CH3OH, under laboratory conditions. We extend this model to simulate the formation of both organic species under interstellar conditions, including freeze-out from the gas and hydrogenation on surfaces. Our aim is to compare calculated abundance ratios with observed values and with the results of prior models. Simulations under different conditions, including density and temperature, have been performed. We find that H2CO and CH3OH form efficiently in cold dense cores or the cold outer envelopes of young stellar objects. The grain mantle is found to have a layered structure with CH3OH on top. The species CO and H2CO are found to exist predominantly in the lower layers of ice mantles where they are not available for hydrogenation at late times.

preprint2009arXiv

Physical Conditions in the Ionized Gas of 30 Doradus

We present a mid-infrared spectroscopic data cube of the central part of 30 Doradus, observed with Spitzer's IRS and MIPS/SED mode. Aromatic dust emission features and emission lines from molecular and atomic hydrogen are detected but not particularly strong. The dominant spectral features are emission lines from moderately ionized species of argon, neon, and sulphur, which are used to determine the physical conditions in the ionized gas. The ionized gas excitation shows strong variations on parsec scales, some of which can plausibly be associated with individual hot stars. We fit the ionic line strengths with photoionization and shock models, and find that photoionization dominates in the region. The ionization parameter U traces the rim of the central bubble, as well as highlighting isolated sources of ionization, and at least one quiescent clump. The hardness of the ionizing radiation field T_rad reveals several "hot spots" that are either the result of individual very hot stars or trace the propagation of the diffuse ionizing field through the surrounding neutral cloud. Consistent with other measurements of giant molecular hydrogen regions, log(U) ranges between -3 and -0.75, and T_rad between 30000 and 85000K.

preprint2009arXiv

Polycyclic aromatic hydrocarbons with SPICA

Thanks to high sensitivity and angular resolution and broad spectral coverage, SPICA will offer a unique opportunity to better characterize the nature of polycyclic aromatic hydrocarbons (PAHs) and very small grains (VSGs), to better use them as probes of astrophysical environments. The angular resolution will enable to probe the chemical frontiers in the evolution process from VSGs to neutral PAHs, to ionized PAHs and to "Grand-PAHs" in photodissotiation regions and HII regions, as a function of G$_0$/n (UV radiation field / density). High sensitivity will favor the detection of the far-IR skeletal emission bands of PAHs, which provide specific fingerprints and could lead to the identification of individual PAHs. This overall characterization will allow to use PAH and VSG populations as tracers of physical conditions in spatially resolved protoplanetary disks and nearby galaxies (using mid-IR instruments), and in high redshift galaxies (using the far-IR instrument), thanks to the broad spectral coverage SPICA provides. Based on our previous experience with ISO and Spitzer we discuss how these goals can be reached.

preprint2009arXiv

The Long-Lived Disks in the Eta Chamaeleontis Cluster

We present IRS spectra and revised MIPS photometry for the 18 members of the Eta Chamaeleontis cluster. Aged 8 Myr, the Eta Cha cluster is one of the few nearby regions within the 5-10 Myr age range, during which the disk fraction decreases dramatically and giant planet formation must come to an end. For the 15 low-mass members, we measure a disk fraction ~50%, high for their 8 Myr age, and 4 of the 8 disks lack near-IR excesses, consistent with the empirical definition of "transition'' disks. Most of the disks are comparable to geometrically flat disks. The comparison with regions of different ages suggests that at least some of the "transition" disks may represent the normal type of disk around low-mass stars. Therefore, their flattened structure and inner holes may be related to other factors (initial masses of the disk and the star, environment, binarity), rather than to pure time evolution. We analyze the silicate dust in the disk atmosphere, finding moderate crystalline fractions (~10-30%) and typical grain sizes ~1-3 micron, without any characteristic trend in the composition. These results are common to other regions of different ages, suggesting that the initial grain processing occurs very early in the disk lifetime (<1 Myr). Large grain sizes in the disk atmosphere cannot be used as a proxy for age, but are likely related to higher disk turbulence. The dust mineralogy varies between the 8-12micron and the 20-30 micron features, suggesting high temperature dust processing and little radial mixing. Finally, the analysis of IR and optical data on the B9 star Eta Cha reveals that it is probably surrounded by a young debris disk with a large inner hole, instead of being a classical Be star.

preprint2009arXiv

Unusual Dust Emission from Planetary Nebulae in the Magellanic Clouds

We present a Spitzer Space Telescope spectroscopic study of a sample of 25 planetary nebulae in the Magellanic Clouds. The low-resolution modules are used to analyze the dust features present in the infrared spectra. This study complements a previous work by the same authors where the same sample was analyzed in terms of neon and sulfur abundances. Over half of the objects (14) show emission of polycyclic aromatic hydrocarbons, typical of carbon-rich dust environments. We compare the hydrocarbon emission in our objects to those of Galactic HII regions and planetary nebulae, and LMC/SMC HII regions. Amorphous silicates are seen in just two objects, enforcing the now well-known-fact that oxygen-rich dust is less common at low metallicities. Besides these common features, some planetary nebulae show very unusual dust. Nine objects show a strong silicon carbide feature at 11um and twelve of them show magnesium sulfide emission starting at 25um. The high percentage of spectra with silicon carbide in the Magellanic Clouds is not common. Two objects show a broad band which may be attributed to hydrogenated amorphous carbon and weak low-excitation atomic lines. It is likely that these nebulae are very young. The spectra of the remaining eight nebulae are dominated by the emission of fine-structure lines with a weak continuum due to thermal emission of dust, although in a few cases the S/N in the spectra is low, and weak dust features may not have been detected.

preprint2007arXiv

Infrared Molecular Starburst Fingerprints in Deeply Obscured (U)LIRG Nuclei

High resolution spectra of the Spitzer Space Telescope show vibration-rotation absorption bands of gaseous C2H2, HCN, and CO2 molecules toward a sample of deeply obscured (U)LIRG nuclei. The observed bands reveal the presence of dense (n>~ 10^7 cm^-3), warm (T_ex = 200-700 K) molecular gas with high column densities of these molecules ranging from a few 10^15 - 10^17 cm^-2. Abundances relative to H2, inferred from the silicate optical depth, range from ~10^-7 to 10^-6 and show no correlation with temperature. Theoretical studies show that the high abundances of both C2H2 and HCN exclude a X-ray dominated region (XDR) associated with the toroid surrounding an AGN as the origin of this dense warm molecular gas. Galactic massive protostars in the so-called Hot Core phase have similar physical characteristics with comparable high abundances of C2H2, HCN, and CO2 in the hot phase. However, the abundances of C2H2 and HCN and the C2H2/CO2 and HCN/CO2 ratios are much higher toward the (U)LIRGs in the cooler (T_ex <= 400 K) phase. We suggest that the warm dense molecular gas revealed by the mid-IR absorption lines is associated with a phase of deeply embedded star formation where the extreme pressures and densities of the nuclear starburst environment have inhibited the expansion of HII regions and the global disruption of the star forming molecular cloud cores, and `trapped' the star formation process in an `extended' Hot Core phase.

preprint2006arXiv

Massive-Star Supernovae as Major Dust Factories

We present late-time optical and mid-infrared observations of the Type-II supernova 2003gd in NGC 628. Mid-infrared excesses consistent with cooling dust in the ejecta are observed 499-678 days after outburst, and are accompanied by increasing optical extinction and growing asymmetries in the emission-line profiles. Radiative-transfer models show that up to 0.02 solar masses of dust has formed within the ejecta, beginning as early as 250 days after outburst. These observations show that dust formation in supernova ejecta can be efficient and that massive-star supernovae can be major dust producers throughout the history of the Universe.

preprint2006arXiv

Spitzer Survey of the Large Magellanic Cloud, Surveying the Agents of a Galaxy's Evolution (SAGE) I: Overview and Initial Results

We are performing a uniform and unbiased, ~7x7 degrees imaging survey of the Large Magellanic Cloud (LMC), using the IRAC and MIPS instruments on board the Spitzer Space Telescope in order to survey the agents of a galaxy's evolution (SAGE), the interstellar medium (ISM) and stars in the LMC. The detection of diffuse ISM with column densities >1.2x10^21 H cm^-2 permits detailed studies of dust processes in the ISM. SAGE's point source sensitivity enables a complete census of newly formed stars with masses >3 solar masses that will determine the current star formation rate in the LMC. SAGE's detection of evolved stars with mass loss rates >1x10^-8 solar masses per year will quantify the rate at which evolved stars inject mass into the ISM of the LMC. The observing strategy includes two epochs in 2005, separated by three months, that both mitigate instrumental artifacts and constrain source variability. The SAGE data are non-proprietary. The data processing includes IRAC and MIPS pipelines and a database for mining the point source catalogs, which will be released to the community in support of Spitzer proposal cycles 4 and 5. We present initial results on the epoch 1 data with a special focus on the N79 and N83 region. The SAGE epoch 1 point source catalog has ~4 million sources. The point source counts are highest for the IRAC 3.6 microns band and decrease dramatically towards longer wavelengths consistent with the fact that stars dominate the point source catalogs and that the dusty objects, e.g. young stellar objects and dusty evolved stars that detected at the longer wavelengths, are rare in comparison. We outline a strategy for identifying foreground MW stars, that may comprise as much as 18% of the source list, and background galaxies, that may comprise ~12% of the source list.

preprint2002arXiv

The rich 6 to 9 um spectrum of interstellar PAHs

Here, we present 6-9 um spectra of a sample of RNe, HII regions, YSOs, evolved stars and galaxies that show strong unidentified infrared bands, obtained with the ISO-SWS. The IR emission features show pronounced variations. 1) The 6.2 um feature shifts and clearly shows profile variations. 2) The 7.7 um complex is comprised of at least two subpeaks peaking at 7.6 and one longwards of 7.7 um. In some cases the main peak can apparently shift up to 8 um. Two sources do not exhibit a 7.7 um complex but instead show a broad emission feature at 8.22 um. 3) The 8.6 um feature has a symmetric profile in all sources and some sources exhibit this band at slightly longer wavelengths. For the 6.2, 7.7 and 8.6 um features, the sources have been classified independently based on their profile and peak position. The classes derived for these features are directly linked with each other. Furthermore, the observed 6--9 um spectrum depends on the type of object. We attribute the observed 6.2 um profile and peak position to the combined effect of a PAH family and anharmonicity with pure PAHs representing the 6.3 um component and substituted/complexed PAHs representing the 6.2 um component. The 7.6 um component is well reproduced by both pure and substituted/complexed PAHs but the 7.8 um component remains an enigma. In addition, the exact identification of the 8.22 um feature remains unknown. The observed variations in the characteristics of the IR emission bands are linked to the local physical conditions. Possible formation and evolution processes that may influence the interstellar PAH class are highlighted.

A. G. G. M. Tielens

Quick read

Decide how to stay connected

How to connect with this researcher

Open a focused conversation when the fit is right

See the researcher in context

Building this graph slice

59 published item(s)

Breaking Orion&#39;s Veil bubble with fossil outflows

Dents in the Veil: Protostellar feedback in Orion

Dust Formation in Astrophysical Environments: The Importance of Kinetics

FEEDBACK from the NGC7538 HII region

Filamentary structures of ionized gas in Cygnus X

Polycyclic Aromatic Hydrocarbon emission model in photodissociation regions I: Application to the 3.3, 6.2, and 11.2 $μ$m bands

Self-absorption in [CII], $^{12}$CO, and HI in RCW120. Building up a geometrical and physical model of the region

The Infrared Database of Extragalactic Observables from Spitzer. -- II. The Database & The Diagnostic Power of Crystalline Silicate Features in Galaxy Spectra

The SOFIA FEEDBACK Legacy Survey: Dynamics and mass ejection in the bipolar HII region RCW 36

[CII] 158$μ$m emission from Orion A. II. Photodissociation region physics

IR photofragmentation of the Phenyl Cation: Spectroscopy and Fragmentation Pathways

Breakdown products of gaseous polycyclic aromatic hydrocarbons investigated with infrared ion spectroscopy

Expanding bubbles in Orion A: [CII] observations of M42, M43, and NGC 1977

FEEDBACK: a SOFIA Legacy Program to Study Stellar Feedback in Regions of Massive Star Formation

Molecular globules in the Veil bubble of Orion. IRAM 30m 12CO, 13CO, and C18O 2-1 expanded maps of Orion A

The H$_2$O Spectrum of the Massive Protostar AFGL 2136 IRS 1 from 2 to 13 $μ$m at High Resolution: Probing the Circumstellar Disk

Searching for the largest bound atoms in space

A bimodal dust grain distribution in the IC 434 HII region

The Physics of the Cold Neutral Medium: Low-frequency Carbon Radio Recombination Lines with the Square Kilometre Array

Depletion of chlorine into HCl ice in a protostellar core

Discovery of Carbon Radio Recombination Lines in M82

Dust composition and mass-loss return from the luminous blue variable R71 in the LMC

Exact Bound-Bound Gaunt Factor Values for Quantum Levels up to n=2000

Herschel Planetary Nebula Survey (HerPlaNS) - First Detection of OH+ in Planetary Nebulae

Location and sizes of forsterite grains in protoplanetary disks: interpretation from the Herschel DIGIT programme

Polycyclic Aromatic Hydrocarbons with armchair edges and the 12.7 μm band

Pore evolution in interstellar ice analogues: simulating the effects of temperature increase

Radiation-pressure-driven dust waves inside bursting interstellar bubbles

The Mass-Loss Return from Evolved Stars to the Large Magellanic Cloud II: Dust Properties for Oxygen-Rich Asymptotic Giant Branch Stars

The enigmatic nature of the circumstellar envelope and bow shock surrounding Betelgeuse as revealed by Herschel. I. Evidence of clumps, multiple arcs, and a linear bar-like structure

Fusion between ground- and space-based mid-infrared observations

Polycyclic aromatic hydrocarbon processing in a hot gas

The Spitzer Spectroscopic Survey of the Small Magellanic Cloud (S4MC): Probing the Physical State of Polycyclic Aromatic Hydrocarbons in a Low-Metallicity Environment

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

Cold dust clumps in dynamically hot gas

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

Determining dust temperatures and masses in the Herschel era: the importance of observations longward of 200 micron

Excitation of the molecular gas in the nuclear region of M82

HD 95881: A gas rich to gas poor transition disk?

Herschel Inventory of The Agents of Galaxy Evolution (HERITAGE): the Large Magellanic Cloud dust

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

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

Polycyclic aromatic hydrocarbon processing in interstellar shocks

SPITZER-IRS spectral fitting of discs around binary post-AGB stars - Corrigendum

Spitzer&#39;s view on aromatic and aliphatic hydrocarbon emission in Herbig Ae stars

The SAGE-Spec Spitzer Legacy program: The life-cycle of dust and gas in the Large Magellanic Cloud

A Spitzer Space Telescope far-infrared spectral atlas of compact sources in the Magellanic Clouds. I. The Large Magellanic Cloud

Chemical depletion in the Large Magellanic Cloud: RV Tauri stars and the photospheric feedback from their dusty discs

Dust coagulation and fragmentation in molecular clouds. I. How collisions between dust aggregates alter the dust size distribution

Ice chemistry in embedded young stellar objects in the Large Magellanic Cloud

Microscopic simulation of methanol and formaldehyde ice formation in cold dense cores

Physical Conditions in the Ionized Gas of 30 Doradus

Polycyclic aromatic hydrocarbons with SPICA

The Long-Lived Disks in the Eta Chamaeleontis Cluster

Unusual Dust Emission from Planetary Nebulae in the Magellanic Clouds

Infrared Molecular Starburst Fingerprints in Deeply Obscured (U)LIRG Nuclei

Massive-Star Supernovae as Major Dust Factories

Spitzer Survey of the Large Magellanic Cloud, Surveying the Agents of a Galaxy&#39;s Evolution (SAGE) I: Overview and Initial Results

The rich 6 to 9 um spectrum of interstellar PAHs

Breaking Orion's Veil bubble with fossil outflows

Spitzer's view on aromatic and aliphatic hydrocarbon emission in Herbig Ae stars

Spitzer Survey of the Large Magellanic Cloud, Surveying the Agents of a Galaxy's Evolution (SAGE) I: Overview and Initial Results