Dense Photodissociation Regions Research Articles

First detection of CF+ in the Large Magellanic Cloud

Context. CF+ has been established as a valuable diagnostic tool for investigating photodissociation regions (PDRs) and fluorine abundances in the Milky Way. However, its role in extragalactic environments remains largely uncharted. Aims. Our objective is to explore the significance of CF+ in the Large Magellanic Cloud (LMC) and assess its utility as a probe for examining C+ and fluorine abundances in external galaxies. Methods. We performed pointed CF+ observations toward an active star-forming region, N113 in the LMC, using the Atacama Pathfinder EXperiment 12 m submillimeter telescope. Results. We report the first discovery of CF+ in the LMC through the successful detection of the CF+ (2→1) and (3→2) lines. The excitation models indicate that CF+ emission originates from dense PDRs characterized by an H2 number density of (0.5–7.9) × 104 cm−3 in N113. Our observations provide the first constraint on the fluorine abundance in molecular clouds in the LMC, ≲1.7 × 10−9. This value is about an order of magnitude lower than those previously measured toward red giants in the LMC, indicative of fluorine deficiency in the molecular gas. The estimated column density ratio between C+ and CF+ appears to be lower than the anticipated equilibrium ratio derived from the fluorine abundance in red giants. Both phenomena can be explained by the deficiency of CF+ caused by the freeze-out of its primary chemical precursor, HF, onto dust grains. Conclusions. The deficiency of CF+ within molecular clouds suggests that the measurements presented in this work serve exclusively as conservative estimates, establishing lower bounds for both the fluorine abundance and C+ column densities in external galaxies.

PDRs4All

Context. The James Webb Space Telescope (JWST) has captured the most detailed and sharpest infrared (IR) images ever taken of the inner region of the Orion Nebula, the nearest massive star formation region, and a prototypical highly irradiated dense photo-dissociation region (PDR). Aims. We investigate the fundamental interaction of far-ultraviolet (FUV) photons with molecular clouds. The transitions across the ionization front (IF), dissociation front (DF), and the molecular cloud are studied at high-angular resolution. These transitions are relevant to understanding the effects of radiative feedback from massive stars and the dominant physical and chemical processes that lead to the IR emission that JWST will detect in many Galactic and extragalactic environments. Methods. We utilized NIRCam and MIRI to obtain sub-arcsecond images over ~150″ and 42″ in key gas phase lines (e.g., Pa α, Br α, [FeII] 1.64 µm, H2 1−0 S(1) 2.12 µm, 0–0 S(9) 4.69 µm), aromatic and aliphatic infrared bands (aromatic infrared bands at 3.3–3.4 µm, 7.7, and 11.3 µm), dust emission, and scattered light. Their emission are powerful tracers of the IF and DF, FUV radiation field and density distribution. Using NIRSpec observations the fractional contributions of lines, AIBs, and continuum emission to our NIRCam images were estimated. A very good agreement is found for the distribution and intensity of lines and AIBs between the NIRCam and NIRSpec observations. Results. Due to the proximity of the Orion Nebula and the unprecedented angular resolution of JWST, these data reveal that the molecular cloud borders are hyper structured at small angular scales of ~0.1–1″ (~0.0002–0.002 pc or ~40–400 au at 414 pc). A diverse set of features are observed such as ridges, waves, globules and photoevaporated protoplanetary disks. At the PDR atomic to molecular transition, several bright features are detected that are associated with the highly irradiated surroundings of the dense molecular condensations and embedded young star. Toward the Orion Bar PDR, a highly sculpted interface is detected with sharp edges and density increases near the IF and DF. This was predicted by previous modeling studies, but the fronts were unresolved in most tracers. The spatial distribution of the AIBs reveals that the PDR edge is steep and is followed by an extensive warm atomic layer up to the DF with multiple ridges. A complex, structured, and folded H0/H2 DF surface was traced by the H2 lines. This dataset was used to revisit the commonly adopted 2D PDR structure of the Orion Bar as our observations show that a 3D “terraced” geometry is required to explain the JWST observations. JWST provides us with a complete view of the PDR, all the way from the PDR edge to the substructured dense region, and this allowed us to determine, in detail, where the emission of the atomic and molecular lines, aromatic bands, and dust originate. Conclusions. This study offers an unprecedented dataset to benchmark and transform PDR physico-chemical and dynamical models for the JWST era. A fundamental step forward in our understanding of the interaction of FUV photons with molecular clouds and the role of FUV irradiation along the star formation sequence is provided.

Cosmic ray-driven magnetohydrodynamic waves in magnetized self-gravitating dusty molecular clouds

ABSTRACT The impact of galactic cosmic rays (CRs) in terms of CR pressure and parallel CR diffusion has been investigated on the low-frequency magnetohydrodynamic (MHD) waves and linear gravitational instability in the typical dusty plasma environment of molecular clouds (MCs). The dusty fluid model is formulated by combining the equations of the magnetized electrons/ions and dust particles, including the CR effects. The interactions between CR fluid and gravitating magnetized dusty plasma have been studied with the help of modified dispersion properties of the MHD waves and instabilities using the hydrodynamic fluid–fluid (CR–plasma) approach. CR diffusion affects the coupling of CR pressure-driven mode with dust-Alfvén MHD mode and causes damping in the MHD waves. It persists in its effect along the direction of the magnetic field and is diminished across the magnetic field. The phase-speed diagram shows that for super-Alfvénic wave, the slow mode becomes the intermediate Alfvén mode. The fundamental Jeans instability criterion remains unaffected due to CR effects, but in the absence of CR diffusion, the effects of dust-acoustic speed and CR pressure-driven wave speed are observed in the instability criterion. It is found that CR pressure stabilizes while CR diffusion destabilizes the growth rates of Jeans instability and significantly affects the gravitational collapse of dusty MCs. The charged dust grains play a dominant role in the sub-Alfvénic and super-Alfvénic MHD waves and the collapse of MCs, triggering gravitational instability. The consequences have been discussed to understand the gravitational instability in the dense photodissociation regions of dusty MCs.

Globules and pillars in Cygnus X

IRAS 20319+3958 in Cygnus X South is a rare example of a free-floating globule (mass ~240M⊙, length ~1.5 pc) with an internal H IIregion created by the stellar feedback of embedded intermediate-mass stars, in particular, one Herbig Be star. In Schneider et al. 2012, (A&A, 542, L18) and Djupvik et al. 2017, (A&A, 599, A37), we proposed that the emission of the far-infrared (FIR) lines of [C II] at 158 μm and [O I] at 145 μm in the globule head are mostly due to an internal photodissociation region (PDR). Here, we present aHerschel/HIFI [C II] 158 μm map of the whole globule and a large set of other FIR lines (mid-to high-JCO lines observed withHerschel/PACS and SPIRE, the [O I] 63 μm line and the12CO 16→15 line observed with upGREAT on SOFIA), covering the globule head and partly a position in the tail. The [C II] map revealed that the whole globule is probably rotating. Highly collimated, high-velocity [C II] emission is detected close to the Herbig Be star. We performed a PDR analysis using the KOSMA-τPDR code for one position in the head and one in the tail. The observed FIR lines in the head can be reproduced with a two-component model: an extended, non-clumpy outer PDR shell and a clumpy, dense, and thin inner PDR layer, representing the interface between the H IIregion cavity and the external PDR. The modelled internal UV field of ~2500G°is similar to what we obtained from theHerschelFIR fluxes, but lower than what we estimated from the census of the embedded stars. External illumination from the ~30 pc distant Cyg OB2 cluster, producing an UV field of ~150–600G°as an upper limit, is responsible for most of the [C II] emission. For the tail, we modelled the emission with a non-clumpy component, exposed to a UV-field of around 140G°.

The magnetic field in the dense photodissociation region of DR 21

ABSTRACT Measuring interstellar magnetic fields is extremely important for understanding their role in different evolutionary stages of interstellar clouds and star formation. However, detecting the weak field is observationally challenging. We present measurements of the Zeeman effect in the 1665 and 1667 MHz (18 cm) lines of the hydroxyl radical (OH) lines towards the dense photodissociation region (PDR) associated with the compact H ii region DR 21 (Main). From the OH 18 cm absorption, observed with the Karl G. Jansky Very Large Array, we find that the line-of-sight magnetic field in this region is ∼0.13 mG. The same transitions in maser emission towards the neighbouring DR 21(OH) and W 75S-FR1 regions also exhibit the Zeeman splitting. Along with the OH data, we use [C ii] 158 μm line and hydrogen radio recombination line data to constrain the physical conditions and the kinematics of the region. We find the OH column density to be ∼3.6 × 1016(Tex/25 K) cm−2, and that the 1665 and 1667 MHz absorption lines are originating from the gas where OH and C+ are co-existing in the PDR. Under reasonable assumptions, we find the measured magnetic field strength for the PDR to be lower than the value expected from the commonly discussed density–magnetic field relation while the field strength values estimated from the maser emission are roughly consistent with the same. Finally, we compare the magnetic field energy density with the overall energetics of DR 21’s PDR and find that, in its current evolutionary stage, the magnetic field is not dynamically important.

ATLASGAL-selected massive clumps in the inner Galaxy

Aims. We study ten molecular transitions obtained from an unbiased 3 mm molecular line survey using the IRAM 30 m telescope toward 409 compact dust clumps identified by the APEX Telescope Large Area Survey of the Galaxy (ATLASGAL) to better understand the photodissociation regions (PDRs) associated with these clumps. The main goal of this study is to investigate whether the abundances of the selected molecules show any variations resulting from the PDR chemistry in different clump environments. Methods. We selected HCO, HOC+, C2H, c-C3H2, CN, H13CN, HC15N, and HN13C as PDR tracers, and H13CO+ and C18O as dense gas tracers. By using estimated optical depths of C2H and H13CN and assuming optically thin emission for other molecular transitions, we derived the column densities of those molecules and their abundances. To assess the influence of the presence and strength of ultra-violet radiation, we compare abundances of three groups of the clumps: HII regions, infrared bright non-HII regions, and infrared dark non-HII regions. Results. We detected C18O, H13CO+, C2H, c-C3H2, CN, and HN13C toward most of the observed dust clumps (detection rate >94%), and H13CN is also detected with a detection rate of 75%. On the other hand, HCO and HC15N show detection rates of 32 and 39%, respectively, toward the clumps, which are mostly associated with HII region sources: detection rates of HCO and HC15N toward the HII regions are 66 and 79%. We find that the abundances of HCO, CN, C2H, and c-C3H2 decrease as the H2 column density increases, indicating high visual extinction, while those of high-density tracers (i.e., H13CO+ and HC15N) are constant. In addition, N(HCO)/N(H13CO+) ratios significantly decrease as H2 column density increases, and, in particular, 82 clumps have X(HCO) ≳ 10−10 and N(HCO)/N(H13CO+) ≳ 1, which are indications of far-ultraviolet (FUV) chemistry. This suggests the observed HCO abundances are likely associated with FUV radiation illuminating the PDRs. We also find that high N(c-C3H2)/N(C2H) ratios found for HII regions that have high HCO abundances (≳10−10) are associated with more evolved clumps with high Lbol/Mclump. This trend might be associated with grain-surface processes, which determine the initial abundances of these molecules, and time-dependent effects in the clumps corresponding to the envelopes around dense PDRs and HII regions. In addition, some fraction of the measured abundances of the small hydrocarbons of the HII sources may be the result of the photodissociation of PAH molecules.

The PDR structure and kinematics around the compact H ii regions S235 A and S235 C with [C ii], [13C ii], [O i], and HCO+ line profiles

ABSTRACT The aim of this work is to study structure and gas kinematics in the photodissociation regions (PDRs) around the compact H ii regions S235 A and S235 C. We observe the [C ii], [13C ii], and [O i] line emission, using SOFIA/upGREAT, and complement them by data of HCO+ and CO. We use the [13C ii] line to measure the optical depth of the [C ii] emission, and find that the [C ii] line profiles are influenced by self-absorption, while the [13C ii] line remains unaffected by these effects. Hence, for dense PDRs, [13C ii] emission is a better tracer of gas kinematics. The optical depth of the [C ii] line is up to 10 in S235 A. We find an expanding motion of the [C ii]-emitting layer of the PDRs into the front molecular layer in both regions. Comparison of the gas and dust columns shows that gas components visible neither in the [C ii] nor in low-J CO lines may contribute to the total column across S235 A. We test whether the observed properties of the PDRs match the predictions of spherical models of expanding H ii region + PDR + molecular cloud. Integrated intensities of the [13C ii], [C ii], and [O i] lines are well represented by the model, but the models do not reproduce the double-peaked [C ii] line profiles due to an insufficient column density of C+. The model predicts that the [O i] line could be a more reliable tracer of gas kinematics, but the foreground self-absorbing material does not allow using it in the considered regions.

Carbon radio recombination lines from gigahertz to megahertz frequencies towards Orion A

Context. The combined use of carbon radio recombination lines (CRRLs) and the 158μm-[CII] line is a powerful tool for the study of the energetics and physical conditions (e.g., temperature and density) of photodissociation regions (PDRs). However, there are few observational studies that exploit this synergy. Aims. Here we explore the relation between CRRLs and the 158μm-[CII] line in light of new observations and models. Methods. We present new and existing observations of CRRLs in the frequency range 0.15–230 GHz with ALMA, VLA, the GBT, Effelsberg 100m, and LOFAR towards Orion A (M 42). We complement these observations with SOFIA observations of the 158μm-[CII] line. We studied two PDRs: the foreground atomic gas, known as the Veil, and the dense PDR between the HII region and the background molecular cloud. Results. In the Veil we are able to determine the gas temperature and electron density, which we use to measure the ionization parameter and the photoelectric heating efficiency. In the dense PDR, we are able to identify a layered PDR structure at the surface of the molecular cloud to the south of the Trapezium cluster. There we find that the radio lines trace the colder portion of the ionized carbon layer, the C+/C/CO interface. By modeling the emission of the 158μm-[CII] line and CRRLs as arising from a PDR we derive a thermal pressure >5 × 107 K cm−3 and a radiation field G0 ≈ 105 close to the Trapezium. Conclusions. This work provides additional observational support for the use of CRRLs and the 158μm-[CII] line as complementary tools to study dense and diffuse PDRs, and highlights the usefulness of CRRLs as probes of the C+/C/CO interface.

Formation of interstellar SH+from vibrationally excited H2: Quantum study of S++ H2⇄ SH++ H reaction and inelastic collision

The rate constants for the formation, destruction, and collisional excitation of SH+are calculated from quantum mechanical approaches using two new SH+2potential energy surfaces (PESs) of4A″ and2A″ electronic symmetry. The PESs were developed to describe all adiabatic states correlating to the SH+(3Σ−) + H(2S) channel. The formation of SH+through the S++ H2reaction is endothermic by ≈9860 K, and requires at least two vibrational quanta on the H2molecule to yield significant reactivity. Quasi-classical calculations of the total formation rate constant for H2(v = 2) are in very good agreement with the quantum results above 100 K. Further quasi-classical calculations are then performed forv = 3, 4, and 5 to cover all vibrationally excited H2levels significantly populated in dense photodissociation regions (PDR). The new calculated formation and destruction rate constants are two to six times larger than the previous ones and have been introduced in the Meudon PDR code to simulate the physical and illuminating conditions in the Orion bar prototypical PDR. New astrochemical models based on the new molecular data produce four times larger SH+column densities, in agreement with those inferred from recent ALMA observations of the Orion bar.

CNO Emission of an Unlensed Submillimeter Galaxy at z = 4.3

Abstract We present the results from Atacama Large Millimeter/submillimeter Array observations of [N ii] 205 μm, [C ii] 158 μm, and [O iii] 88 μm lines in an unlensed submillimeter galaxy at z = 4.3, COSMOS-AzTEC-1, hosting a compact starburst core with an effective radius of ∼1 kpc. The [C ii] and [N ii] emission are spatially resolved in 0.″3-resolution (1 kpc in radius). The kinematic properties of the [N ii] emission are consistent with those of the CO(4–3) and [C ii] emission, suggesting that the ionized gas feels the same gravitational potential as the associated molecular gas and photodissociation regions (PDRs). On the other hand, the spatial extent is different among the lines and dust continuum: the [C ii] emitting gas is the most extended and the dust is the most compact, leading to a difference of the physical conditions in the interstellar medium. We derive the incident far-ultraviolet flux and the hydrogen gas density through PDR modeling by properly subtracting the contribution of ionized gas to the total [C ii] emission. The observed [C ii] emission is likely produced by dense PDRs with cm−3 and G 0 = 103.5–3.75 in the central 1 kpc region and cm−3 and G 0 = 103.25–3.5 in the central 3 kpc region. We have also successfully measured the line ratio of [O iii]/[N ii] in the central 3 kpc region of COSMOS-AzTEC-1 at z = 4.3, which is the highest redshift where both nitrogen and oxygen lines are detected. Under the most likely physical conditions, the measured luminosity ratio of L [O iii]/L [N ii] = 6.4 ± 2.2 indicates a near solar metallicity with Z gas = 0.7–1.0 Z ⊙, suggesting a chemically evolved system at z = 4.3.

Photoinduced polycyclic aromatic hydrocarbon dehydrogenation

The physical and chemical conditions in photodissociation regions (PDRs) are largely determined by the influence of far ultraviolet radiation. Far-UV photons can efficiently dissociate molecular hydrogen, a process that must be balanced at the HI/H2 interface of the PDR. Given that reactions involving hydrogen atoms in the gas phase are highly inefficient under interstellar conditions, H2 formation models mostly rely on catalytic reactions on the surface of dust grains. Additionally, molecular hydrogen formation in polycyclic aromatic hydrocarbons (PAHs) through the Eley–Rideal mechanism has been considered as well, although it has been found to have low efficiency in PDR fronts. In a previous work, we have described the possibility of efficient H2 release from medium to large sized PAHs upon photodissociation, with the exact branching between H-/H2-loss reactions being molecule dependent. Here, we investigate the astrophysical relevance of this process, by using a model for the photofragmentation of PAHs under interstellar conditions. We focus on three PAHs cations (coronene, ovalene, and circumcoronene), which represent three possibilities in the branching of atomic and molecular hydrogen losses. We find that, for ovalene (H2-loss dominated) the rate coefficient for H2 formation reaches values of the same order as H2 formation in dust grains. This result suggests that this hitherto disregarded mechanism can account, at least partly, for the high level of molecular hydrogen formation in dense PDRs.

The ALMA view of UV-irradiated cloud edges: unexpected structures and processes

AbstractFar-UV photons (FUV, E < 13.6 eV) from hot massive stars regulate, or at least influence, the heating, ionization, and chemistry of most of the neutral interstellar medium (H i and H2 clouds). Investigating the interaction between FUV radiation and interstellar matter (molecules, atoms and grains) thus plays an important role in astrochemistry.The Orion Bar, an interface region between the Orion A molecular cloud and the H ii region around the Trapezium cluster, is a textbook example of a strongly illuminated dense PDR (photodissociation region). The Bar is illuminated by a FUV field of a few 104 times the mean interstellar radiation field. Because of its proximity and nearly edge-on orientation, it provides a very good template to investigate the chemical content, structure, and dynamics of a strongly irradiated molecular cloud edge. We have used ALMA to mosaic a small field of the Bar where the critical transition from atomic to molecular gas takes place. These observations provide an unprecedented sharp view of this transition layer (≲ 1″ resolution or ≲ 414 AU). The resulting images (so far in the rotational emission of CO, HCO+, H13CO+, SO+, SO, and reactive ions SH+ and HOC+) show the small-scale structure in gas density and temperature, and the steep abundance gradients. The images reveal a pattern of high-density substructures, photo-ablative gas flows and instabilities at the edge of the molecular cloud. These first ALMA images thus show a more complex morphology than the classical clump/interclump static model of a PDR.In order to quantify the chemical content in strongly FUV-irradiated gas, we have also used the IRAM-30 m telescope to carry out a complete line-survey of the illuminated edge of the Bar in the millimeter domain. Our observations reveal the presence of complex organic molecules (and precursors) that were not expected in such a harsh environment. In particular, we have reported the first detection of the unstable cis conformer of formic acid (HCOOH) in the ISM. The energy barrier to internal rotation (the conversion from trans to cis) is approximately 4827 cm−1 (≈7000 K). Hence, this detection is surprising. The low inferred trans-to-cis abundance ratio of 2.8±1.0 supports a photoswitching mechanism: a given conformer absorbs a FUV stellar photon that radiatively excites the molecule to electronic states above the interconversion barrier. Subsequent fluorescent decay leaves the molecule in a different conformer form. This mechanism, which we have specifically studied with ab initio quantum calculations, was not considered so far in astrochemistry although it can affect the structure of a variety of molecules in PDRs.

EVIDENCE FROM SOFIA IMAGING OF POLYCYCLIC AROMATIC HYDROCARBON FORMATION ALONG A RECENT OUTFLOW IN NGC 7027

ABSTRACT We report spatially resolved (FWHM ∼ 3.″8–4.″6) mid-IR imaging observations of the planetary nebula (PN) NGC 7027 taken with the 2.5 m telescope on board the Stratospheric Observatory for Infrared Astronomy (SOFIA). Images of NGC 7027 were acquired at 6.3, 6.6, 11.1, 19.7, 24.2, 33.6, and 37.1 using the Faint Object Infrared Camera for the SOFIA Telescope (FORCAST). The observations reveal emission from polycyclic aromatic hydrocarbons (PAHs) and warm dust ( K) from the illuminated inner edge of the molecular envelope surrounding the ionized gas and central star. The DustEM code was used to fit the spectral energy distribution of fluxes obtained by FORCAST and the archival infrared spectrum of NGC 7027 acquired by the Short Wavelength Spectrometer (SWS) on the Infrared Space Observatory (ISO). Best-fit dust models provide a total dust mass of , where carbonaceous large (a = 1.5 μm) and very small ( ) grains, and PAHs ( ) compose 96.5, 2.2, and 1.3% of the dust by mass, respectively. The 37 μm optical depth map shows minima in the dust column density at regions in the envelope that are coincident with a previously identified collimated outflow from the central star. The optical depth minima are also spatially coincident with enhancements in the 6.2 μm PAH feature, which is derived from the 6.3 and 6.6 μm maps. We interpret the spatial anti-correlation of the dust optical depth and PAH 6.2 μm feature strength and their alignment with the outflow from the central star as evidence of dust processing and rapid PAH formation via grain–grain collisions in the post-shock environment of the dense ( ) photo-dissociation region and molecular envelope.

Further ALMA observations and detailed modeling of the Red Rectangle.

We aim to study the rotating and expanding gas in the Red Rectangle, which is a well known object that recently left the asymptotic giant branch (AGB) phase. We analyze the properties of both components and the relation between them. Rotating disks have been very elusive in post-AGB nebulae, in which gas is almost always found to be in expansion. We present new high-quality ALMA observations of C17O J=6-5 and H13CN J=4-3 line emission and results from a new reduction of already published 13CO J=3-2 data. A detailed model fitting of all the molecular line data, including previous maps and single-dish observations of lines of CO, CII, and CI, was performed using a sophisticated code that includes an accurate nonlocal treatment of radiative transfer in 2D. These observations (of low- and high-opacity lines requiring various degrees of excitation) and the corresponding modeling allowed us to deepen the analysis of the nebular properties. We also stress the uncertainties, particularly in the determination of the boundaries of the CO-rich gas and some properties of the outflow. We confirm the presence of a rotating equatorial disk and an outflow, which is mainly formed of gas leaving the disk. The mass of the disk is ~ 0.01 M⊙, and that of the CO-rich outflow is around ten times smaller. High temperatures of ≳ 100 K are derived for most components. From comparison of the mass values, we roughly estimate the lifetime of the rotating disk, which is found to be of about 10000 yr. Taking data of a few other post-AGB composite nebulae into account, we find that the lifetimes of disks around post-AGB stars typically range between 5000 and more than 20000 yr. The angular momentum of the disk is found to be high, ~ 9 M⊙ AU km s-1, which is comparable to that of the stellar system at present. Our observations of H13CN show a particularly wide velocity dispersion and indicate that this molecule is only abundant in the inner Keplerian disk, at ≲ 60 AU from the stellar system. We suggest that HCN is formed in a dense photodissociation region (PDR) due to the UV excess known to be produced by the stellar system, following chemical mechanisms that are well established for interstellar medium PDRs and disks orbiting young stars. We further suggest that this UV excess could lead to an efficient formation and excitation of PAHs and other C-bearing macromolecules, whose emission is very intense in the optical counterpart.

AN IN-DEPTH VIEW OF THE MID-INFRARED PROPERTIES OF POINT SOURCES AND THE DIFFUSE ISM IN THE SMC GIANT H II REGION, N66

(abridged)The focus of this work is to study mid-infrared point sources and the diffuse interstellar medium (ISM) in the low-metallicity (~0.2 solar) giant HII region N66 using the Spitzer Space Telescope's Infrared Spectrograph. We study 14 targeted infrared point sources as well as spectra of the diffuse ISM that is representative of both the photodissociation regions (PDRs) and the HII regions. Among the point source spectra, we spectroscopically confirm that the brightest mid-infrared point source is a massive embedded young stellar object, we detect silicates in emission associated with two young stellar clusters, and we observe spectral features of a known B[e] supergiant that are more commonly associated with Herbig Be stars. In the diffuse ISM, we provide additional evidence that the very small grain population is being photodestroyed in the hard radiation field. The 11.3 um PAH complex emission exhibits an unexplained centroid shift in both the point source and ISM spectra that should be investigated at higher signal-to-noise and resolution. Unlike studies of other regions, the 6.2 um and 7.7 um band fluxes are decoupled; the data points cover a large range of I7.7/I11.3 PAH ratio values within a narrow band of I6.2/I11.3 ratio values. Furthermore, there is a spread in PAH ionization, being more neutral in the dense PDR where the radiation field is relatively soft, but ionized in the diffuse ISM/PDR. By contrast, the PAH size distribution appears to be independent of local ionization state. Important to unresolved studies of extragalactic low-metallicity star-forming regions, we find that emission from the infrared-bright point sources accounts for only 20-35% of the PAH emission from the entire region. These results make a comparative dataset to other star-forming regions with similarly hard and strong radiation fields.

Carbon fractionation in photo-dissociation regions

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

A complete model of CH+rotational excitation including radiative and chemical pumping processes

Aims. Excitation of far-infrared and submillimetric molecular lines may originate from nonreactive collisions, chemical formation, or far infrared, near-infrared, and optical fluorescences. As a template, we investigate the impact of each of these processes on the excitation of the methylidyne cation CH+ and on the intensities of its rotational transitions recently detected in emission in dense photodissociation regions (PDRs) and in planetary nebulae. Methods. We have developed a nonlocal thermodynamic equilibrium (non-LTE) excitation model that includes the entire energy structure of CH+, i.e. taking into account the pumping of its vibrational and bound and unbound electronic states by near-infrared and optical photons. The model includes the theoretical cross-sections of nonreactive collisions with H, H2, He, and e-, and a Boltzmann distribution is used to describe the probability of populating the excited levels of CH+ during its chemical formation by hydrogenation of C+. To confirm our results we also performed an extensive analytical study, which we use to predict the main excitation process of several diatomic molecules, namely HF, HCl, SiO, CS, and CO. Results. At densities nH = 10^4 cm-3, the excitation of the rotational levels of CH+ is dominated by the radiative pumping of its electronic, vibrational, and rotational states if the intensities of the radiation field at \sim 0.4, \sim 4, and \sim 300 \mum are stronger than 10^5, 10^8, and 10^4 times those of the local interstellar radiation field (ISRF). Below these values, the chemical pumping is the dominant source of excitation of the J > 1 levels, even at high kinetic temperatures (\sim 1000 K). The far-infrared emission lines of CH+ observed in the Orion Bar and the NGC 7027 PDRs are consistent with the predictions of our excitation model assuming an incident far-ultraviolet (FUV) radiation field of \sim 3 \times 10^4 (in Draine's unit) and densities of \sim 5 \times 10^4 and \sim 2 \times 10^5 cm-3. In the case of NGC 7027, the estimate of the density is 10 to 100 times lower than those deduced by traditional excitation codes. Applying our model to other X1\Sigma+ ground state diatomic molecules, we find that HF, and SiO and HCl are the species the most sensitive to the radiative pumping of their vibrational and bound electronic states. In both cases, the minimal near-infrared and optical/UV radiation field intensities required to modify their rotational level populations are \sim 10^3 times those of the local ISRF at densities nH = 10^4 cm-3. All these results point towards interstellar and circumstellar media with densities lower than previously established and cast doubts on the clumpiness of well-studied molecular clouds.

The nature of the interstellar medium of the starburst low-metallicity galaxy Haro 11: a multi-phase model of the infrared emission

(abridged) Our goal is to describe the multi-phase ISM of the IR bright low-metallicity galaxy Haro 11, dissecting the photoionised and photodissociated gas components. We present observations of the mid- and far-IR fine-structure cooling lines obtained with the Spitzer/IRS and Herschel/PACS spectrometers. We use the spectral synthesis code Cloudy to methodically model the ionised and neutral gas from which these lines originate. We find that the mid- and far-IR lines account for ~1% of the total IR luminosity L_TIR. Haro 11 is undergoing a phase of intense star formation, as traced by the brightest line [OIII] 88um, with L_[OIII]/L_TIR ~0.3%, and high ratios of [NeIII]/[NeII] and [SIV]/[SIII]. Due to their different origins, the observed lines require a multi-phase modeling comprising: a compact HII region, dense fragmented photodissociation regions (PDRs), a diffuse extended low-ionisation/neutral gas which has a volume filling factor of at least 90%, and porous warm dust in proximity to the stellar source. For a more realistic picture of the ISM of Haro 11 we would need to model the clumpy source and gas structures. We combine these 4 model components to explain the emission of 17 spectral lines, investigate the global energy balance of the galaxy through its spectral energy distribution, and establish a phase mass inventory. While the ionic emission lines of Haro 11 essentially originate from the dense HII region component, a diffuse low-ionisation gas is needed to explain the [NeII], [NII], and [CII] line intensities. The [OIII] 88um line intensity is not fully reproduced by our model, hinting towards the possible presence of yet another low-density high-ionisation medium. The [OI] emission is consistent with a dense PDR of low covering factor, and we find no evidence for an X-ray dominated component. The PDR component accounts for only 10% of the [CII] emission.

Evaporating very small grains as tracers of the UV radiation field in photo-dissociation regions

Context. In photo-dissociation regions (PDRs), Polycyclic Aromatic Hydrocarbons (PAHs) could be produced by evaporation of Very Small Grains (VSGs) by the impinging UV radiation field from a nearby star. Aims. We investigate quantitatively the transition zone between evaporating Very Small Grains (eVSGs) and PAHs in several PDRs. Methods. We study the relative contribution of PAHs and eVSGs to the mid-IR emission in a wide range of excitation conditions. We fit the observed mid-IR emission of PDRs by using a set of template band emission spectra of PAHs, eVSGs and gas lines. The fitting tool PAHTAT (PAH Toulouse Astronomical Templates) is made available to the community as an IDL routine. From the results of the fit, we derive the fraction of carbon f_eVSG locked in eVSGs and compare it to the intensity of the local UV radiation field. Results. We show a clear decrease of f_eVSG with increasing intensity of the local UV radiation field, which supports the scenario of photo-destruction of eVSGs. Conversely, this dependence can be used to quantify the intensity of the UV radiation field for different PDRs, including non resolved ones. Conclusions. PAHTAT can be used to trace the intensity of the local UV radiation field in regions where eVSGs evaporate, which correspond to relatively dense (nH = [100, 10^5 ] cm-3) and UV irradiated PDRs (G0 = [100, 5x10^4]) where H2 emits in rotational lines.

SPITZERMAPPING OF POLYCYCLIC AROMATIC HYDROCARBON AND H2FEATURES IN PHOTODISSOCIATION REGIONS

The mid-infrared (MIR) spectra of dense photodissociation regions (PDRs) are typically dominated by emission from polycyclic aromatic hydrocarbons (PAHs) and the lowest pure rotational states of molecular hydrogen (H2); two species which are probes of the physical properties of gas and dust in intense UV radiation fields. We utilize the high angular resolution of the Infrared Spectrograph on the Spitzer Space Telescope to construct spectral maps of the PAH and H2 features for three of the best studied PDRs in the galaxy, NGC 7023, NGC 2023 and IC 63. We present spatially resolved maps of the physical properties, including the H2 ortho-to-para ratio, temperature, and G_o/n_H. We also present evidence for PAH dehydrogenation, which may support theories of H2 formation on PAH surfaces, and a detection of preferential self-shielding of ortho-H2. All PDRs studied exhibit average temperatures of ~500 - 800K, warm H2 column densities of ~10^20 cm^-2, G_o/n_H ~ 0.1 - 0.8, and ortho-to-para ratios of ~ 1.8. We find that while the average of each of these properties is consistent with previous single value measurements of these PDRs, when available, the addition of spatial resolution yields a diversity of values with gas temperatures as high as 1500 K, column densities spanning ~ 2 orders of magnitude, and extreme ortho-to-para ratios of < 1 and > 3.