Chemical stratification in the Orion Bar: JCMT Spectral Legacy Survey observations

Photon-dominated regions (PDRs) are expected to show a layered structure in molecular abundances and emerging line emission, which is sensitive to the physical structure of the region as well as the UV radiation illuminating it. We aim to study this layering in the Orion Bar, a prototypical nearby PDR with a favorable edge-on geometry. We present new maps of 2 by 2 arcminute fields at 14-23 arcsecond resolution toward the Orion Bar in the SO 8_8-9_9, H2CO 5_(1,5)-4_(1,4), 13CO 3-2, C2H 4_(9/2)-3_(7/2) and 4_(7/2)-3_(5/2), C18O 2-1 and HCN 3-2 transitions. The data reveal a clear chemical stratification pattern. The C2H emission peaks close to the ionization front, followed by H2CO and SO, while C18O, HCN and 13CO peak deeper into the cloud. A simple PDR model reproduces the observed stratification, although the SO emission is predicted to peak much deeper into the cloud than observed while H2CO is predicted to peak closer to the ionization front than observed. In addition, the predicted SO abundance is higher than observed while the H2CO abundance is lower than observed. The discrepancies between the models and observations indicate that more sophisticated models, including production of H2CO through grain surface chemistry, are needed to quantitatively match the observations of this region.

💡 Research Summary

The Orion Bar is one of the nearest and best‑studied photon‑dominated regions (PDRs), and its nearly edge‑on orientation makes it an ideal laboratory for testing the predicted layered chemistry that arises when a strong ultraviolet (UV) field impinges on a molecular cloud. In this work the authors present new, high‑resolution (14–23 arcsec) maps of a 2′ × 2′ field covering six key transitions: SO 8₈–9₉, H₂CO 5₁,₅–4₁,₄, ¹³CO 3–2, C₂H 4₉⁄₂–3₇⁄₂ and 4₇⁄₂–3₅⁄₂, C¹⁸O 2–1, and HCN 3–2. The data were obtained as part of the JCMT Spectral Legacy Survey (SLS) and processed to produce integrated intensity maps that can be directly compared with PDR model predictions.

The observational results reveal a clear chemical stratification. The C₂H emission peaks right at the ionization front (IF), indicating that this radical is formed in the very surface layer that is directly exposed to the UV radiation, most likely through photodissociation of larger hydrocarbons or PAH fragments. Moving a few arcseconds deeper into the cloud, the emission from H₂CO and SO reaches its maximum. H₂CO is traditionally thought to arise both from gas‑phase reactions (e.g., CH₃ + O → H₂CO + H) and from the hydrogenation of CO on dust‑grain mantles followed by UV‑induced desorption. The observed offset suggests that grain‑surface formation and subsequent photodesorption are important contributors in the Orion Bar. SO, a sulfur‑bearing molecule, peaks at a similar depth; its formation is efficient at temperatures above ~100 K where atomic S reacts with O or OH, consistent with a moderately warm, partially shielded layer.

Deeper still, the emission from the CO isotopologues (C¹⁸O, ¹³CO) and HCN reaches its maximum. These species require higher visual extinctions (A_V ≈ 5–10 mag) to survive the UV field and are therefore tracers of the denser, more shielded interior of the bar. HCN, in particular, is known to be enhanced in high‑density gas (n > 10⁵ cm⁻³) and moderate temperatures (~50 K), matching the conditions inferred for this region.

To interpret the observations, the authors employed a simple, one‑dimensional steady‑state PDR model (similar to the Meudon or KOSMA‑τ codes) with a density profile appropriate for the Orion Bar and an incident UV field of χ ≈ 10⁴ Draine units. The model reproduces the overall ordering of the peaks (C₂H → H₂CO/SO → CO isotopologues/HCN) but fails in two critical aspects. First, the model predicts that SO should peak much farther inside the cloud than observed, while H₂CO should peak closer to the IF. Second, the absolute abundances are mismatched: the modeled SO column density is 2–3 times higher than measured, whereas the modeled H₂CO column density is roughly half of the observed value.

These discrepancies point to missing physics in the simple model. The most plausible explanation is the neglect of grain‑surface chemistry and UV‑induced desorption processes. H₂CO is known to be efficiently produced on icy mantles via successive hydrogenation of CO (CO → HCO → H₂CO → CH₃OH). In a strong UV field, a fraction of this ice‑bound H₂CO can be released into the gas phase, shifting its emission peak deeper than pure gas‑phase chemistry would predict. Conversely, sulfur chemistry is notoriously uncertain; the model likely overestimates the efficiency of S → SO conversion in the warm interior and underestimates the destruction pathways (e.g., SO + H → S + OH) that become important near the IF.

The authors conclude that more sophisticated modeling is required. A multi‑dimensional (2‑D or 3‑D) treatment that includes realistic density and temperature gradients, coupled with a full grain‑surface network (including CO hydrogenation, sulfur depletion, and photodesorption yields), would be needed to simultaneously reproduce the spatial offsets and absolute abundances of all observed species. Future observations with higher angular resolution (≤5″) from facilities such as ALMA, as well as complementary infrared data from JWST to directly probe ice features, will provide the necessary constraints to refine these models.

In summary, this study delivers the first high‑resolution, multi‑species map of the Orion Bar that clearly demonstrates chemical stratification predicted for PDRs. While a basic PDR model captures the general layering, the mismatches in SO and H₂CO highlight the importance of grain‑surface processes and the need for more complex, physically realistic simulations to fully understand the chemistry of irradiated molecular clouds.