JWST observations of the Horsehead photon-dominated region I. First results from multi-band near- and mid-infrared imaging

Abstract: The JWST has captured the sharpest IR images ever taken of the Horsehead nebula, a prototypical moderately irradiated PDR that is fully representative of most of the UV-illuminated molecular gas in the Milky Way and star-forming galaxies. We investigate the impact of FUV radiation of a molecular cloud and constrain the structure of the edge of the PDR and its illumination conditions. We used NIRCam and MIRI to obtain 17 broadband and 6 narrowband maps from 0.7 to 28 $μ$m. We mapped the dust emission, scattered light, and several gas phase lines. We also used HST-WFC3 maps at 1.1 and 1. 6 $μ$m, along with HST-STIS spectroscopic observations of the H$α$ line. We probed the structure of the edge of the Horsehead and resolved its spatial complexity. We detected a network of faint striated features extending perpendicularly to the PDR front into the H\,II region in filters sensitive to nano-grain emission and light scattered by larger grains. This may indeed figure as the first detection of the entrainment of dust particles in the evaporative flow. The map of the 1-0 S(1) line of H$_2$ presents sharp sub-structures on scales as small as 1.5 arcsec. The ionization and dissociation fronts appear at distances 1-2 arcsec behind the edge of the PDR and seem to spatially coincide, indicating a thickness of the neutral atomic layer below 100 au. All broadband maps present strong color variations which can be explained by dust attenuation. Deviations of the emissions in the H$α$, Pa$α,$ and Br$α$ lines also indicate dust attenuation. With a very simple model, we derive the main features of the extinction curve. A small excess of extinction at 3 $μ$m may be attributed to icy H$_2$O mantles onto grains. In all lines of sight crossing the inner regions of the Horsehead, it appears that dust attenuation is non-negligible over the entire spectral range.

• The paper presents high-resolution JWST imaging that reveals co-spatial ionization and dissociation fronts in the Horsehead PDR with separations under 0.25″.
• It uses data from 17 NIRCam and 9 MIRI filters along with HST observations to detail dust entrainment, attenuation profiles, and molecular tracers at the PDR interface.
• The findings constrain PDR models by demonstrating the impact of oblique FUV illumination and photoevaporative flows on dust dynamics and radiative transfer.

JWST Multi-band Infrared Imaging of the Horsehead PDR: Structural and Attenuation Properties

Introduction

The Horsehead Nebula serves as a canonical photon-dominated region (PDR), characterized by moderate far-ultraviolet (FUV) irradiation and representative of the ultraviolet-illuminated molecular gas in the Milky Way and extragalactic star-forming systems. JWST, with its unprecedented near- and mid-infrared spatial resolution and spectral coverage, enables a transformative analysis of PDR boundaries down to scales of 0.1–1" ($\sim$40–400 AU at 400 pc). This study presents a comprehensive analysis of JWST-NIRCam and MIRI imaging, coupled with HST-WFC3 and STIS observations, to probe the detailed dust and gas structure, extinction processes, and energetic feedback at the Horsehead molecular cloud edge (2404.15816).

Observational Methodology

The study utilizes 17 NIRCam broadband and 6 narrowband filters (0.7–5 $\mu$m), 9 MIRI broadband filters (5–28 $\mu$m), and selected HST imaging/spectroscopy for H$\alpha$ calibration. Mosaics were constructed perpendicular to the PDR front, resolving ionized and molecular tracers (e.g., H$\alpha$, Pa$\alpha$, Br$\alpha$, H$_2$ 1–0 S(1)) and a suite of aromatic/aliphatic dust features and continua. Data reduction involved custom pipeline applications for optimal inter-pixel capacitance correction, precise astrometric alignment (GAIA DR3), and suppression of $1/f$ detector noise. Background correction and continuum subtraction protocols were critical for accurate profile extraction of both emission lines and dust spectral features.

Physical Structure and Stratification of the PDR

Sub-arcsecond Spatial Morphology and Entrained Dust

JWST resolves complex, striated morphology perpendicular to the PDR front, visible in both NIRCam dust-sensitive bands (e.g., F335M, F770W) and HST scattered light. The morphology and photometric ratios indicate detection of dust—possibly nano-grains and PAH-sized carriers—entrained in the photoevaporative outflow. This constitutes compelling evidence for the dynamical feedback exerted by FUV-driven flows on dust population transport within a photoionized interface.

Ionization and Dissociation Fronts: Co-spatiality and Thin Atomic Layer

Brightness profiles in H and H$_2$ lines localize the ionization (IF) and dissociation fronts (DF) to within 1–2" of the molecular cloud edge, with an upper limit on the IF–DF separation below 0.25" (100 AU). The observed spatial coincidence of IF and DF is consistent with models predicting merging of ionization and dissociation fronts under moderate FUV fields and high gas pressure but implies negligible atomic layer thickness, a regime distinct from classical, high-excitation PDRs (e.g., Orion Bar). Substructure in H$_2$ filaments on scales down to 1.5" is resolved, confirming the heterogeneity in local gas density and illumination field at the PDR boundary.

Impact of Oblique Illumination and Attenuation

Profile analyses and color metric gradients indicate that the exciting star ($\sigma$-Orionis) imposes an oblique illumination geometry, illuminating the Horsehead from behind. This configuration amplifies differential dust attenuation effects, responsible for both deviations in recombination line flux scaling (relative to case B) and significant color gradients across the PDR. The attenuation is non-negligible over the entire JWST spectral window and scales with projected column density through the PDR.

Quantitative Attenuation and Extinction Curve Properties

Extraction of Extinction Law Features

Using a thin-screen attenuation model, the spectra of optical depth differences extracted along several spatial cuts reveal:

• A decreasing extinction from 0.7 to $\sim$5 $\mu$m,
• Pronounced excesses at 9.7 $\mu$m and 18 $\mu$m, associated with the Si-O vibrational features,
• A detectable 3 $\mu$m excess in the F300M filter, ascribed to H$_2$O ice mantles on grains in shielded regions.

The derived attenuation profiles are congruent with independent measurements via other tracers (mm continuum, ro-vibrational H$_2$ transitions), substantiating the reliability of the simple screen model at the PDR interface.

Implications for Dust Evolution and Radiative Transfer

The smooth variation of color ratios and the match to empirical extinction curves at most wavelengths provides strong evidence that dust attenuation dominates the observed color gradients, with secondary impacts from variations in nano-grain (PAH) abundance and potentially multiple scattering contributions. The ultra-high angular resolution permits, for the first time, the direct mapping of attenuation profiles across a resolved PDR and highlights the necessity of incorporating attenuation effects even for IR-bright, nearby star-forming environments.

Theoretical and Practical Implications

Feedback, ISM Evolution, and Model Constraints

The results provide stringent empirical constraints on PDR modeling—specifically, the need to account for dynamically thin atomic layers, the role of moderate illumination geometry, and fine-scale variations in dust and gas properties. The observed entrainment of dust in the photoevaporative flow reinforces the paradigm that stellar feedback redistributes, and potentially removes, small dust and PAH populations from star-forming interfaces, with implications for ISM enrichment and star formation regulation. Attenuation properties derived herein will inform the interpretation of spatially unresolved PDRs in external galaxies and the deduction of their intrinsic properties from integrated IR SEDs.

Future Prospects

The imaging presented forms the foundation for follow-up JWST IFU spectroscopic studies (NIRSpec, MIRI-MRS), enabling decomposition of physical conditions (density, temperature, pressure), characterization of dust evolution pathways, and direct spectral disentanglement of emission, absorption, and scattering components. Extension of this analysis to other PDR archetypes will clarify the universality of the observed phenomena and refine theoretical models of irradiated ISM regions in a range of galactic environments.

Conclusion

This JWST program demonstrates the capacity of multi-band near- and mid-IR imaging to resolve and quantify the complex interplay of ionization, dissociation, dust processing, and attenuation at a canonical PDR boundary. The data support a scenario of simultaneous IF/DF alignment, pervasive dust attenuation, and morphological evidence for dust entrainment in photoevaporative flows. Theoretical modeling, synergistic IFU spectroscopy, and comparative studies across PDR environments are warranted as the next steps to unravel the lifecycle of dust and gas under UV feedback in the molecular ISM (2404.15816).