Photodissociation Regions Research Articles

A Q-band Line Survey toward Orion KL Using the Tianma Radio Telescope

We have conducted a line survey toward Orion KL using the Q-band receiver of the Tianma 65 m radio telescope (TMRT), covering 34.8–50 GHz with a velocity resolution between 0.79 and 0.55 km s−1, respectively. The observations reach a sensitivity of the level of 1–8 mK, proving that the TMRT is sensitive for conducting deep-line surveys. In total, 597 Gaussian features are extracted. Among them, 177 radio recombination lines (RRLs) are identified, including 126, 40, and 11 RRLs of hydrogen, helium, and carbon, with a maximum Δn of 16, 7, and 3, respectively. The carbon RRLs are confirmed to originate from photodissociation regions with a V LSR ∼ 9 km s−1. In addition, 371 molecular transitions of 53 molecular species are identified. Twenty-one molecular species of this survey were not firmly detected in the Q band by Rizzo et al., including species such as H2CS, HCOOH, C2H5OH, CO, H2CCO, CH3CHO, CH2OCH2, HCN υ 2 = 1, and CH3OCHO υ t = 1. In particular, the vibrationally excited states of ethyl cyanide (C2H5CN υ13/υ21) are for the first time firmly detected in the Q band. NH3 (15,15) and (16,16) are identified, and they are so far the highest transitions of the NH3 inversion lines detected toward Orion KL. All of the identified lines can be reproduced by a radiative transfer model.

Contribution of polycyclic aromatic hydrocarbon ionization to neutral gas heating in galaxies: model versus observations

The ionization of polycyclic aromatic hydocarbons (PAHs) by ultraviolet (UV) photons from massive stars is expected to account for a large fraction of the heating of neutral gas in galaxies. Evaluation of this proposal, however, has been limited by our ability to directly compare observational diagnostics to the results of a molecular model describing PAH ionization. The objective of this article is to take advantage of the most recent values of molecular parameters derived from laboratory experiments and quantum chemical calculations on PAHs and provide a detailed comparison between modeled values and observational diagnostics for the PAH charge state and the heating efficiency for PAHs. Despite the use of a simple analytical model, we obtain a good agreement between model results and observational diagnostics over a wide range of radiation fields and physical conditions, in environments such as star-forming regions, galaxies, and protoplanetary disks. In addition, we found that the modeled photoelectric heating rates by PAHs are close to the observed cooling rates given by the gas emission. These results show that PAH ionization is the main source of neutral gas heating in these environments. The results of our photoelectric heating model by PAHs can thus be used to assess the contribution of UV radiative heating in galaxies (vs. e.g. shock). We conclude on the importance of implementing the physics of PAH ionization in astrophysical codes, which are developed, for example, for the evaporating surfaces of protoplanetary disks, the diffuse interstellar medium, and the photodissociation regions associated with star-forming regions in the local and distant universe. We provide the empirical formulas and Python code to calculate the heating rates and heating efficiencies for PAHs.

Efficiency of the top-down polycyclic aromatic hydrocarbon-to-fullerene conversion in ultraviolet irradiated environments

ABSTRACT Polycyclic aromatic hydrocarbons (PAHs) and fullerenes play a major role in the physics and chemistry of the interstellar medium (ISM). Based on a number of recent experimental and theoretical investigations we developed a model in which PAHs are subject to photo-dissociation (carbon and hydrogen loss) and hydrogenation. We take into account that dehydrogenated PAHs may fold into closed structures – fullerenes. Fullerenes, in their turn, can be also hydrogenated, becoming fulleranes, and photo-dissociated, losing carbon and hydrogen atoms. The carbon loss leads to shrinking of fullerene cages to smaller ones. We calculate the abundance of PAHs and fullerenes of different sizes and hydrogenation level depending on external conditions: the gas temperature, intensity of radiation field, number density of hydrogen atoms, carbon atoms, and electrons. We highlight the conditions, which are favourable for fullerene formation from PAHs, and we conclude that this mechanism works not only in H-poor environment but also at modest values of hydrogen density up to 104 cm−3. We found that fulleranes can be formed in the ISM, although the fraction of carbon atoms locked in them can be maximum around 10−9. We applied our model to two photo-dissociation regions, Orion Bar and NGC 7023. We compare our estimates of the fullerene abundance and synthetic band intensities in these objects with the observations and conclude that our model gives good results for the closest surroundings of ionizing stars. We also demonstrate that additional fullerene formation channels should operate along with ultraviolet (UV)-induced formation to explain abundance of fullerenes far from UV sources.

Nano-grain depletion in photon-dominated regions

Context. Carbonaceous nano-grains play a fundamental role in the physico-chemistry of the interstellar medium (ISM) and especially of photon-dominated regions (PDRs). Their properties vary with the local physical conditions and affect the local chemistry and dynamics. Aims. We aim to highlight the evolution of carbonaceous nano-grains in three different PDRs and propose a scenario of dust evolution as a response to the physical conditions. Methods. We used Spitzer/IRAC (3.6, 4.5, 5.8, and 8 µm) and Spitzer/MIPS (24 µm) together with Herschel/PACS (70 µm) to map dust emission in IC63 and the Orion Bar. To assess the dust properties, we modelled the dust emission in these regions using the radiative transfer code SOC together with the THEMIS dust model. Results. Regardless of the PDR, we find that nano-grains are depleted and that their minimum size is larger than in the diffuse ISM (DISM), which suggests that the mechanisms that lead nano-grains to be photo-destroyed are very efficient below a given critical size limit. The evolution of the nano-grain dust-to-gas mass ratio with both G0 and the effective temperature of the illuminating star indicates a competition between the nano-grain formation through the fragmentation of larger grains and nano-grain photo-destruction. We modelled dust collisions driven by radiative pressure with a classical 1D approach to show that this is a viable scenario for explaining nano-grain formation through fragmentation and, thus, the variations observed in nano-grain dust-to-gas mass ratios from one PDR to another. Conclusions. We find a broad variation in the nano-grain dust properties from one PDR to another, along with a general trend of nano-grain depletion in these regions. We propose a viable scenario of nano-grain formation through fragmentation of large grains due to radiative pressure-induced collisions.

Disentangling emission from star-forming regions in the Magellanic Clouds: Linking [O iii]λ88 µm and 24 µm

Context. The [O iii]λ88 µm line is observed in many galaxies including our neighboring Magellanic Clouds and is a well-known tracer of H ii regions, while the 24 µm continuum emission has often been used to trace warm dust in the ionized phases of galaxies. The association of both the [O iii]λ88 µm line and 24 µm in galaxies to star formation motivates this study to determine their observational relation. Aims. This study explores the link between the [O iii]λ88 µm and 24 µm continuum in star-forming regions in the Magellanic Clouds. We also explore the local conditions driving the relation between those tracers. Methods. We compared observations with 1D Cloudy models consisting of an H ii region plus a photodissociation region (PDR) component, varying the stellar age, the initial density (at the illuminated edge of the cloud), and the ionization parameter. We introduced a new parameter, cPDR, to quantify the proportion of emission arising from PDRs and that with an origin in H ii regions along each line of sight. We used the ratio ([C ii]+[O i])/[O iii] as a proxy for the ratio of PDR versus H ii region emission, and compared it to the [O iii]/24 µm ratio. The use of [O iii]/24 µm and [O iii]/70 µm together allowed us to constrain the models most efficiently. Results. We find a correlation over at least 3 orders of magnitude in [O iii]λ88 µm and 24 µm continuum. This correlation is seen for spatially resolved maps of the Magellanic Cloud regions as well as unresolved galaxy-wide low metallicity galaxies of the Dwarf Galaxy Survey. We also find that most of the regions have low proportions of PDRs along the lines of sight (<12%), while a limited area of some of the mapped regions can reach 30–50%. For most lines of sight within the star-forming regions we have studied in the Magellanic Clouds, H ii regions are the dominant phase. Conclusions. We propose the use of the correlation between the [O iii]λ88 µm and 24 µm continuum as a new predictive tool to estimate, for example, the [O iii]λ88 µm when the 24 µm continuum is available or inversely. This can be especially useful to prepare for Atacama Large Milimeter Array (ALMA) observations of [O iii]λ88 µm in high-z galaxies. The simple and novel method we developed may also provides a way to disentangle different phases along the line of sight, when other 3D information is not available.

Compression Disturbances in the Atomic Zone of Photodissociation Regions of Interstellar Gas

Compression Disturbances in the Atomic Zone of Photodissociation Regions of Interstellar Gas

From the Circumnuclear Disk in the Galactic Center to thick, obscuring tori of AGNs

The high accretion rates needed to fuel the central black hole in a galaxy can be achieved via viscous torques in thick disks and rings, which can be resolved by millimeter interferometry within the inner ∼20 pc of the active galaxy NGC 1068 at comparable scales and sensitivity to single dish observations of the Circumnuclear Disk (CND) in the Galactic Center. To interpret observations of these regions and determine the physical properties of their gas distribution, we present a modeling effort that includes the following: (i) simple dynamical simulations involving partially inelastic collisions between disk gas clouds; (ii) an analytical model of a turbulent clumpy gas disk calibrated by the dynamical model and observations; (iii) local turbulent and cosmic ray gas heating and cooling via H2O, H2, and CO emission; and (iv) determination of the molecular abundances. We also consider photodissociation regions (PDRs) where gas is directly illuminated by the central engine. We compare the resulting model datacubes of the CO, HCN, HCO+, and CS brightness temperatures to available observations. In both cases the kinematics can be explained by one or two clouds colliding with a preexisting ring, in a prograde sense for the CND and retrograde for NGC 1068. And, with only dense disk clouds, the line fluxes can be reproduced to within a factor of about two. To avoid self-absorption of the intercloud medium, turbulent heating at the largest scales, comparable to the disk height, has to be decreased by a factor of 50–200. Our models indicate that turbulent mechanical energy input is the dominant gas-heating mechanism within the thick gas disks. Turbulence is maintained by the gain of potential energy via radial gas accretion, which is itself enhanced by the collision of the infalling cloud. In NGC 1068, we cannot exclude that intercloud gas significantly contributes to the molecular line emission. In this object, while the bulk of the X-ray radiation of the active galactic nucleus is absorbed in a layer of Compton-thick gas inside the dust sublimation radius, the optical and UV radiation may enhance the molecular line emission from photodissociation regions by ∼50% at the inner edge of the gas ring. Infrared pumping may also increase the HCN(3−2) line flux throughout the gas ring by about a factor of two. Our models support the scenario of infalling gas clouds onto preexisting gas rings in galactic centers, and it is viable and consistent with available observations of the CND in the Galactic Center and the dense gas distribution within the inner 20 pc of NGC 1068.

Spinning Nanoparticles Impacted by C-shock: Implications for Radio-millimeter Emission from Star-forming Regions

We investigate the impact of anomalous microwave emission (AME) on the radio-millimeter spectral energy distribution for three typical interstellar medium (ISM) conditions surrounding star-forming regions—cold neutral medium, warm neutral medium, and photodissociation regions—by comparing the emissivities of three major contributors: free–free, thermal dust emission, and AME. In particular, for spinning nanoparticles (i.e., potential carriers of AME), we consider a known grain destruction mechanism due to a centrifugal force from spin-up processes caused by collisions between dust grains and supersonic neutral streams in a magnetized shock (C-shock). We demonstrate that, if the ISM in a magnetic field is impacted by a C-shock developed by a supernova explosion in the early phase of massive star formation (≲10 Myr), AME can be significantly or almost entirely suppressed relative to free–free and thermal dust continuum emission if the grain tensile strength is small enough. This study may shed light on explaining the rare observations of AME from extragalactic star-forming regions preferentially observed from massive star clusters and suggest a scenario of “the rise and fall of AME” in accordance with the temporal evolution of star-forming regions.

The KOSMA-τ PDR model

Numerical models of photodissociation regions (PDRs) are an essential tool to quantitatively understand observations of massive star forming regions through simulations. Few mature PDR models are available and the Cologne KOSMA-τ PDR model is the only sophisticated model that uses a spherical cloud geometry, thereby allowing us to simulate clumpy PDRs. We present the current status of the code as a reference for modelers and for observers that plan to apply KOSMA-τ to interpret their data. For the numerical solution of the chemical problem, we present a superior Newton-Raphson stepping algorithm and discuss strategies to numerically stabilize the problem and speed up the iterations. The chemistry in KOSMA-τ is upgraded to include the full surface chemistry in an up-to-date formulation and we discuss a novel computation of branching ratios in chemical desorption reactions. The high dust temperature in PDRs leads to a selective freeze-out of oxygen-bearing ice species due to their higher condensation temperatures and we study changes in the ice mantle structures depending on the PDR parameters, in particular the impinging ultraviolet field. Selective freeze-out can produce enhanced C abundances and higher gas temperatures, resulting in a fine-structure line emission of atomic carbon [C i] enhanced by up to 50% if surface reactions are considered. We show how recent Atacama Large Millimeter Array (ALMA) observations of HCO+ emission in the Orion Bar with high spatial resolution on the scale of individual clumps can be interpreted in the context of nonstationary, clumpy PDR ensembles. Additionally, we introduce WL-PDR, a simple plane-parallel PDR model written in Mathematica to act as a numerical testing environment of PDR modeling aspects.

The role of highly vibrationally excited H2 initiating the nitrogen chemistry

The formation of hydrides by gas-phase reactions between H2 and a heavy element atom is a very selective process. Reactions with ground-state neutral carbon, oxygen, nitrogen, and sulfur atoms are very endoergic and have high energy barriers because the H2 molecule has to be fragmented before a hydride bond is formed. In cold interstellar clouds, these barriers exclude the formation of CH, OH, NH, and SH radicals through hydrogen abstraction reactions. Here we study a very energetically unfavorable process, the reaction of N(4S) atoms with H2 molecules. We calculated the reaction rate coefficient for H2 in different vibrational levels, using quantum methods for v = 0−7 and quasi-classical methods up to v =12; for comparison purposes, we also calculated the rate coefficients of the analogous reaction S (3P)+ H2(v) → SH + H. Owing to the high energy barrier, these rate coefficients increase with v and also with the gas temperature. We implemented the new rates in the Meudon photodissociation region (PDR) code and studied their effect on models with different ultraviolet (UV) illumination conditions. In strongly UV-irradiated dense gas (Orion Bar conditions), the presence of H2 in highly vibrationally excited levels (v ≥ 7) enhances the NH abundance by two orders of magnitude (at the PDR surface) compared to models that use the thermal rate coefficient for reaction N(4S) + H2 → NH + H. The increase in NH column density, N(NH), across the PDR is a factor of ~25. We investigate the excitation and detectability of submillimeter NH rotational emission lines. Being a hydride, NH excitation is very subthermal (Trot ≪ Tk) even in warm and dense gas. We explore existing Herschel/HIFI observations of the Orion Bar and Horsehead PDRs. We report a 3σ emission feature at the ~974 GHz frequency of the NH NJ = 12 − 01 line toward the Bar. The emission level implies N(NH) ≃ 1013 cm−2, which is consistent with PDR models using the new rate coefficients for reactions between N and UV-pumped H2. This formation route dominates over hydrogenation reactions involving the less abundant N+ ion. JWST observations will quantify the amount and reactivity of UV-pumped H2 in many interstellar and circumstellar environments.

NIR spectroscopic survey of protostellar jets in the star-forming region IC 1396N

Context. The bright-rimmed cloud IC 1396N, associated with an intermediate-mass star-forming region, hosts a number of CO, molecular hydrogen, and Herbig-Haro (HHs) outflows powered by a set of millimetre compact sources. Aims. The aim of this work is to characterise the kinematics and physical conditions of the H2 emission features spread throughout the IC 1396N region. The features appear as chains of knots with a jet-like morphology and trace different H2 outflows. We also obtain further information about (and an identification of) the driving sources. Methods. Low-resolution, long-slit near-infrared spectra were acquired with the NICS camera at the TNG telescope, using grisms KB (R ~ 1200), HK, and JH (R ~ 500). Several slit pointings and position angles were used throughout the IC 1396N region in order to sample a number of the H2 knots that were previously detected in deep H2 2.12 μm images. Results. The knots exhibit rich ro–vibrational spectra of H2, consistent with shock-excited excitation, from which radial velocities and relevant physical conditions of the IC 1396N H2 outflows were derived. These also allowed estimating extinction ranges towards several features. [Fe ii] emission was only detected towards a few knots that also display unusually high H2 1–0 S(3)/S(1) flux ratios. The obtained radial velocities confirm that most of the outflows are close to the plane of the sky. Nearby knots in the same chain often display different radial velocities, both blue–shifted and red–shifted, which we interpret as due to ubiquitous jet precession in the driving sources or the development of oblique shocks. One of the chains (strand A, i.e. knots A1 to A15) appears as a set of features trailing a leading bow-shock structure consistent with the results of 3D magneto-hydrodynamical models. The sides of the leading bow shock (A15) exhibit different radial velocities. We discuss possible explanations. Our data cannot confirm whether strands A and B have both originated in the intermediate mass young stellar object [BGE2002] BIMA 2 because a simple model of a precessing jet cannot account for their locations. Conclusions. Near-infrared spectroscopy has confirmed that most of the H2 ro-vibrational emission in IC 1396N is shock-excited rather than uv-excited in photon-dominated regions. It has shown a complex kinematical structure in most strands of emitting knots as well.

Fragmented atomic shell around S187 H ii region and its interaction with molecular and ionized gas

ABSTRACT The environment of S187, a nearby H ii region (1.4 ± 0.3 kpc), is analyzed. A surrounding shell has been studied in the H i line, molecular lines, and also in infrared and radio continua. We report the first evidence of a clumpy H i environment in its photodissociation region. A background radio galaxy enables the estimation of the properties of cold atomic gas. The estimated atomic mass fraction of the shell is ∼260 M⊙, the median spin temperature is ∼50 K, the shell size is ∼4 pc with typical wall width around 0.2 pc. The atomic shell consists of ∼100 fragments. The fragment sizes correlate with mass with a power-law index of 2.39–2.50. The S187 shell has a complex kinematical structure, including the expanding quasi-spherical layer, molecular envelope, an atomic sub-bubble inside the shell and two dense cores (S187 SE and S187 NE) at different stages of evolution. The atomic sub-bubble inside the shell is young, contains a Class II young stellar object and OH maser in the centre and the associated YSOs in the walls of the bubble. S187 SE and S187 NE have similar masses (∼1200 and ∼900 M⊙, respectively). S187 SE is embedded into the atomic shell and has a number of associated objects, including high-mass protostars, outflows, maser sources, and other indicators of ongoing star formation. No YSOs inside S187 NE were detected, but indications of compression and heating by the H ii region exist.

Methanol masers in NGC 253 with ALCHEMI

Context. Methanol masers of Class I (collisionally pumped) and Class II (radiatively pumped) have been studied in great detail in our Galaxy in a variety of astrophysical environments such as shocks and star-forming regions and are they are helpful to analyze the properties of the dense interstellar medium. However, the study of methanol masers in external galaxies is still in its infancy. Aims. Our main goal is to search for methanol masers in the central molecular zone (CMZ; inner 500 pc) of the nearby starburst galaxy NGC 253. Methods. Covering a frequency range between 84 and 373 GHz (λ = 3.6–0.8 mm) at high angular (1.″6 ∼ 27 pc) and spectral (∼8–9 km s−1) resolution with ALCHEMI (ALMA Comprehensive High-resolution Extragalactic Molecular Inventory), we have probed different regions across the CMZ of NGC 253. In order to look for methanol maser candidates, we employed the rotation diagram method and a set of radiative transfer models. Results. We detect for the first time masers above 84 GHz in NGC 253, covering an ample portion of the J−1 → (J − 1)0 − E line series (at 84, 132, 229, and 278 GHz) and the J0 → (J − 1)1 − A series (at 95, 146, and 198 GHz). This confirms the presence of the Class I maser line at 84 GHz, which was already reported, but now being detected in more than one location. For the J−1 → (J− 1)0 − E line series, we observe a lack of Class I maser candidates in the central star-forming disk. Conclusions. The physical conditions for maser excitation in the J−1 → (J − 1)0 − E line series can be weak shocks and cloud-cloud collisions as suggested by shock tracers (SiO and HNCO) in bi-symmetric shock regions located in the outskirts of the CMZ. On the other hand, the presence of photodissociation regions due to a high star-formation rate would be needed to explain the lack of Class I masers in the very central regions.

FAST Search for Circumstellar Atomic Hydrogen. I. The Young Planetary Nebula IC 4997

Using the Five-hundred-meter Aperture Spherical radio Telescope in Guizhou, China, we detect the 21 cm neutral atomic hydrogen absorption in the young planetary nebula IC 4997. The absorption arises from a shell that is also associated with Na i D lines. The H i shell has a mass of 1.46 × 10−2 M ⊙ and a dynamic age of 990 yr. The column density of H i is estimated to be 7.1 × 1020 cm−2, which can be well explained in terms of a photodissociation region around the ionized nebula, limited by the self-shielding of H2. We find that the atomic-to-ionized hydrogen ratio is 0.6, suggesting that H i substantially contributes to the overall nebular mass.

Photodissociation and X-Ray-Dominated Regions

The radiation from stars and active galactic nuclei (AGNs) creates photodissociation regions (PDRs) and X-ray-dominated regions (XDRs), where the chemistry or heating is dominated by far-ultraviolet (FUV) radiation or X-ray radiation, respectively. PDRs include a wide range of environments, from the diffuse interstellar medium (ISM) to dense star-forming regions. XDRs are found in the center of galaxies hosting AGNs, in protostellar disks, and in the vicinity of X-ray binaries. In this review, we describe the dominant thermal, chemical, and radiation transfer processes in PDRs and XDRs, as well as give a brief description of models and their use for analyzing observations. We then present recent results from Milky Way, nearby extragalactic, and high-redshift observations. Several important results include the following: ▪ Velocity-resolved PDR lines reveal the kinematics of the neutral atomic gas and provide constraints on the stellar feedback process. Their interpretation is, however, in dispute, as observations suggest a prominent role for stellar winds, whereas they are much less important in theoretical models. ▪ A significant fraction of molecular mass resides in CO-dark gas especially in low-metallicity and/or highly irradiated environments. ▪ The CO ladder and [Ci]/[Cii] ratios can determine if FUV or X rays dominate the ISM heating of extragalactic sources. ▪ With Atacama Large Millimeter/submillimeter Array, PDR and XDR tracers are now routinely detected on galactic scales over cosmic time. This makes it possible to link the star-formation history of the Universe to the evolution of the physical and chemical properties of the gas.

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

ABSTRACT 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, G0, electron density ne, and the gas temperature, Tgas). 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 $\gamma = G_{0}\times T_{\rm gas}^{1/2}/n_{e}$. The anions emerge as the dominant charge carriers in low γ (<2 × 102) environments, neutrals in the intermediate γ (103 – 104) environments, and cations in the high γ (>105) 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 interstellar medium, 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) versus 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.

PDRs4All: A JWST Early Release Science Program on Radiative Feedback from Massive Stars

Massive stars disrupt their natal molecular cloud material through radiative and mechanical feedback processes. These processes have profound effects on the evolution of interstellar matter in our Galaxy and throughout the universe, from the era of vigorous star formation at redshifts of 1–3 to the present day. The dominant feedback processes can be probed by observations of the Photo-Dissociation Regions (PDRs) where the far-ultraviolet photons of massive stars create warm regions of gas and dust in the neutral atomic and molecular gas. PDR emission provides a unique tool to study in detail the physical and chemical processes that are relevant for most of the mass in inter- and circumstellar media including diffuse clouds, proto-planetary disks, and molecular cloud surfaces, globules, planetary nebulae, and star-forming regions. PDR emission dominates the infrared (IR) spectra of star-forming galaxies. Most of the Galactic and extragalactic observations obtained with the James Webb Space Telescope (JWST) will therefore arise in PDR emission. In this paper we present an Early Release Science program using the MIRI, NIRSpec, and NIRCam instruments dedicated to the observations of an emblematic and nearby PDR: the Orion Bar. These early JWST observations will provide template data sets designed to identify key PDR characteristics in JWST observations. These data will serve to benchmark PDR models and extend them into the JWST era. We also present the Science-Enabling products that we will provide to the community. These template data sets and Science-Enabling products will guide the preparation of future proposals on star-forming regions in our Galaxy and beyond and will facilitate data analysis and interpretation of forthcoming JWST observations.

Breaking Orion's Veil with fossil outflows

Context. 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 [C II] 158 μm observations have revealed an expanding bubble, the Veil shell, that is powered by stellar winds and ionization feedback. Aims. We have identified a protrusion-like substructure in the northwestern portion of the Orion Veil shell that may indicate additional, highly directional feedback mechanisms. Our goal is to investigate the origin of the protrusion by quantifying its possible driving mechanisms. Methods. We use the [C II] 158 μm map of the Orion Nebula obtained with the upGREAT instrument on board SOFIA. The spectral and spatial resolution of the observations are 0.3 km s−1 and 16", respectively. The velocity-resolved [C II] observations allow us to construct position-velocity (pv) diagrams to measure the morphology and the expansion velocity of the protrusion. For the morphology, we also use new observations of 12CO and 13CO J = 2-1 (to trace molecular gas), Spitzer 8 μm observations (to trace the far-ultraviolet-illuminated surfaces of photodissociation regions), and Hα observations (to trace ionized gas). For the kinematics, we perform a line-profile analysis of [C II], 13CO, and 12CO at 12 positions, covering the entire protrusion. To quantify the stellar feedback, we estimate the mass of the protrusion by fitting the dust thermal emission. We compare the kinetic energy with the stellar wind of θ1 Ori C and the momentum of the outflows of massive protostars to investigate the driving mechanism of the protrusion. Results. The pv diagrams reveal two half-shells expanding at velocities of +6 km s−1 and +12 km s−1. We find that the protrusion has a diameter of ~1.3 pc with a ~45 M⊙ shell expanding at +12 km s−1 at the northwestern rim of the Veil. The thickness of the expanding shell is ~0.1 pc. Using the mass in the limb-brightened shell and the maximum expansion velocity, we calculate the kinetic energy and the momentum of the protrusion to be ~7 × 1046 erg and 540 M⊙ km s−1, respectively. We consider three possible origins for this protrusion: fossil outflow cavities created by jets and outflows during the protostellar accretion phase, preexisting "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 preexisting 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 Orion's Veil owing to the photo-ablation from the walls of the protrusion. Conclusions. We conclude that the outflows of massive protostars can influence the morphology of the future H II region and even cause breakages in the ionization front. Specifically, the interaction of stellar winds of main-sequence stars with the molecular core preprocessed by the protostellar jet is important.

Physical and Chemical Properties of the Molecular Gas Associated with the Mid-infrared Bubble S156

Abstract Using archival infrared data of GLIMPSE, Hi-GAL, and molecular line data of SEDIGISM, MSGPCOS, and MALT90, we investigate the physical and chemical properties of the molecular gas associated with the mid-infrared bubble S156. By the method of spectral energy distribution, we made H2 column density and dust temperature maps of this region. We find two clouds with masses of 5.4 ± 1.1 × 104 M ⊙ and 2.2 ± 0.5 × 104 M ⊙, respectively. In both of the two clouds, the 13CO (2–1/1–0) and 13CO (2–1)/N2H+ (1–0) intensity ratios are enhanced on the boundary. Cloud A has a clear dust temperature gradient decreasing from the boundary to the outside region. Our analysis indicates cloud A is mainly influenced by the feedback from S156, while cloud B is affected both by S156 and the G305 complex. We also find the 13CO and C2H emissions tend to be brighter in photon dominated regions (PDRs), while N2H+ is brighter in the regions of cold gas that is far away. Furthermore, we make the abundance maps of C2H and N2H+ of cloud A. We find the abundance of C2H is enhanced in the region facing ionizing stars and it decreases steadily moving away from them. On the other hand, the abundance of N2H+ increases from the ionized boundary to the cold gas outward. These phenomena indicate C2H prefers to be produced in hot gas such as PDRs, while N2H+ could be destroyed by it. Our study also suggests the abundance ratio of C2H to N2H+ could trace PDRs in the late stages of massive star formation.

FEEDBACK from the NGC 7538 H II region

Context. The interaction of expanding H II regions with their environmental clouds is one of the central questions driving the Stratospheric Observatory for Infrared Astronomy (SOFIA) legacy program FEEDBACK. Aims. We want to understand the interaction of the prototypical NGC 7538 H II region with the neighboring molecular cloud hosting several active star-forming regions. Methods. Using the SOFIA, we mapped an area of ~210′2 (~125 pc2) around NGC 7538 in the velocity-resolved ionized carbon fine-structure line [CII] at 1.9 THz (158 μm). Complementary observed atomic carbon [CI] at 492 GHz and high-J CO(8–7) data, as well as archival near- and far-infrared, cm continuum, CO(3–2), and HI data are folded into the analysis. Results. The ionized carbon [CII] data reveal rich morphological and kinematic structures. While the overall morphology follows the general ionized gas that is also visible in the radio continuum emission, the channel maps show multiple bubble-like structures with sizes on the order of ~80–100″ (~1.0–1.28 pc). While at least one of them may be an individual feedback bubble driven by the main exciting sources of the NGC 7538 H II region (the O3 and O9 stars IRS6 and IRS5), the other bubble-like morphologies may also be due to the intrinsically porous structure of the H II region. An analysis of the expansion velocities around 10 km s−1 indicates that thermal expansion is not sufficient but that wind-driving from the central O-stars is required. The region exhibits a general velocity gradient across, but we also identify several individual velocity components. 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 C+, atomic C0 and molecular carbon CO show a layered structure in this PDR. The carbon in the PDR is dominated by its ionized C+ form with atomic C0 and molecular CO masses of ~0.45 ± 0.1 M⊙ and ~1.2 ± 0.1 M⊙, respectively, compared to the ionized carbon C+ in the range of 3.6−9.7 M⊙. This bar-shaped PDR exhibits a velocity-gradient across, indicating motions along the line of sight toward the observer. Conclusions. Even if it is shown to be dominated by two nearby exciting sources (IRS6 and IRS5), the NGC 7538 H II region exhibits a diverse set of substructures that interact with each other as well as with the adjacent cloud. Compared to other recent [CII] observations of H II regions (e.g., Orion Veil, RCW120, RCW49), the bubble-shape morphologies revealed in [CII] emission that are indicative of expanding shells are recurring structures of PDRs.