In this paper we study the filamentary substructure of 3.3 $\mu$m PAH emission from JWST/NIRCam observations in the base of the M82 star-burst driven wind. We identify plume-like substructure within the PAH emission with widths of $\sim$50 pc. Several of the plumes extend to the edge of the field-of-view, and thus are at least 200-300 pc in length. In this region of the outflow, the vast majority ($\sim$70\%) of PAH emission is associated with the plumes. We show that those structures contain smaller scale "clouds" with widths that are $\sim$5-15 pc, and they are morphologically similar to the results of "cloud-crushing" simulations. We estimate the cloud-crushing time-scales of $\sim$0.5-3 Myr, depending on assumptions. We show this time scale is consistent with a picture in which these observed PAH clouds survived break-out from the disk rather than being destroyed by the hot wind. The PAH emission in both the midplane and the outflow is shown to tightly correlate with that of Pa$\alpha$ emission (from HST/NICMOS data), at the scale of both plumes and clouds, though the ratio of PAH-to-Pa$\alpha$ increases at further distances from the midplane. Finally, we show that the outflow PAH emission is suppressed in regions of the M82 wind that are bright in X-ray emission. Overall, our results are broadly consistent with a picture in which cold gas in galactic outflows is launched via hierarchically structured plumes, and those small scale clouds are more likely to survive the wind environment when collected into the larger plume structure.

We introduce a prescription for estimating the flux of the 7.7 micron and 11.3 micron polycyclic aromatic hydrocarbon (PAH) features from broadband JWST/MIRI images. Probing PAH flux with MIRI imaging data has advantages in field of view, spatial resolution, and sensitivity compared with MIRI spectral maps, but comparisons with spectra are needed to calibrate these flux estimations over a wide variety of environments. For 267 MIRI/MRS spectra from independent regions in the four luminous infrared galaxies (LIRGs) in the Great Observatories All-sky LIRG Survey (GOALS) early release science program, we derive synthetic filter photometry and directly compare estimated PAH fluxes to those measured from detailed spectral fits. We find that for probing PAH 7.7 micron, the best combination of filters is F560W, F770W, and either F1500W or F2100W, and the best for PAH 11.3 micron is F560W, F1000W, F1130W, and F1500W. The prescription with these combinations yields predicted flux densities that typically agree with values from spectral decomposition within 7 PAH 7.7 and 11.3 micron, respectively.

We measured the spectral properties of a sample of 20 galaxies at z similar to 0.35 selected for having surprisingly red JWST/NIRCAM F200W-F444W colors. Of these, 19 galaxies were observed with JWST/NIRSpec in the PRISM configuration, while the remaining galaxy was observed with the high-resolution gratings. Of the 20 galaxies in our sample, 17 exhibit strong 3.3 mu m polycyclic aromatic hydrocarbon (PAH) emission (equivalent width (EW) (PAH(3.3)) >= 0.03 mu m). In these galaxies, the strength of the color excess does not depend on environment and correlates with EW(PAH(3.3)). Nonetheless, the presence of the PAH(3.3) alone cannot fully explain the color excess, as an EW of similar to 0.1 mu m is able to increase the color of galaxies by only 0.13 mag. A contribution from a hot dust component is required to explain the excess. Both the PAH(3.3) EW and flux correlate with the H alpha EW and flux, suggesting that they are produced by the same mechanism. Five of the galaxies of our sample showing PAH(3.3) would be classified as passive based on broadband rest frame colors ((B-V) and/or UVJ diagrams) and are hence "faux passive". Of these, three galaxies have a significantly lower EW(PAH(3.3)) given their color and also have low EW(H alpha), and we tentatively conclude that this behavior is due to the presence of an active galactic nucleus. The three galaxies with no PAH(3.3) in emission have passive spectra, as do the eight galaxies in our sample with normal F200W-F444W colors. We therefore conclude that the PAH(3.3) feature is linked to dust-enshrouded star formation. The dust-corrected star formation rate (SFR) from PAH(3.3) is a factor of 3.5 higher than the SFR obtained from H alpha, suggesting that these galaxies are characterized by significant amounts of dust.

We present a study of new 7.7-11.3 μm data obtained with the James Webb Space Telescope Mid-InfraRed Instrument in the starburst galaxy M 82. In particular, we focus on the dependency of the integrated CO(1-0) line intensity on the MIRI-F770W and MIRI-F1130W filter intensities to investigate the correlation between H_2 content and the 7.7 and 11.3 μm features from polycyclic aromatic hydrocarbons (PAH) in M 82's outflows. To perform our analysis, we identify CO clouds using the archival ^12CO(J=1-0) NOEMA moment 0 map within 2 kpc from the center of M 82, with sizes ranging between ∼21 and 270 pc; then, we compute the CO-to-PAH relations for the 306 validated CO clouds. On average, the power-law slopes for the two relations in M 82 are lower than what is seen in local main-sequence spirals. In addition, there is a moderate correlation between I_ 7.7μ m /I_ 11.3μ m for most of the CO cloud groups analyzed in this work. Our results suggest that the extreme conditions in M 82 translate into CO not tracing the full budget of molecular gas in smaller clouds, perhaps as a consequence of photoionization and/or emission suppression of CO molecules due to hard radiation fields from the central starburst.

Context.Mid-infrared emission features are important probes of the properties of ionized gas and hot or warm molecular gas, which are difficult to probe at other wavelengths. The Orion Bar photodissociation region (PDR) is a bright, nearby, and frequently studied target containing large amounts of gas under these conditions. Under the “PDRs4All” Early Release Science Program for JWST, a part of the Orion Bar was observed with MIRI integral field unit (IFU) spectroscopy, and these high-sensitivity IR spectroscopic images of very high angular resolution (0.2″) provide a rich observational inventory of the mid-infrared (MIR) emission lines, while resolving the HIIregion, the ionization front, and multiple dissociation fronts.Aims.We list, identify, and measure the most prominent gas emission lines in the Orion Bar using the new MIRI IFU data. An initial analysis summarizes the physical conditions of the gas and demonstrates the potential of these new data and future IFU observations with JWST.Methods.The MIRI IFU mosaic spatially resolves the substructure of the PDR, its footprint cutting perpendicularly across the ionization front and three dissociation fronts. We performed an up-to-date data reduction, and extracted five spectra that represent the ionized, atomic, and molecular gas layers. We identified the observed lines through a comparison with theoretical line lists derived from atomic data and simulated PDR models. The identified species and transitions are summarized in the main table of this work, with measurements of the line intensities and central wavelengths.Results.We identified around 100 lines and report an additional 18 lines that remain unidentified. The majority consists of HIrecombination lines arising from the ionized gas layer bordering the PDR. The HIline ratios are well matched by emissivity coefficients from H recombination theory, but deviate by up to 10% because of contamination by HeIlines. We report the observed emission lines of various ionization stages of Ne, P, S, Cl, Ar, Fe, and Ni. We show how the NeIII/NeII, SIV/SIII, and ArIII/ArIIratios trace the conditions in the ionized layer bordering the PDR, while FeIII/FeIIand NiIII/NiIIexhibit a different behavior, as there are significant contributions to FeIIand NiIIfrom the neutral PDR gas. We observe the pure-rotational H2lines in the vibrational ground state from 0–0S(1) to 0–0S(8), and in the first vibrationally excited state from 1–1S(5) to 1–1 S(9). We derive H2excitation diagrams, and for the three observed dissociation fronts, the rotational excitation can be approximated with one thermal (~700 K) component representative of an average gas temperature, and one nonthermal component (~2700 K) probing the effect of UV pumping. We compare these results to an existing model of the Orion Bar PDR, and find that the predicted excitation matches the data qualitatively, while adjustments to the parameters of the PDR model are required to reproduce the intensity of the 0–0 S (6) to S (8) lines.

We present new JWST-NIRSpec IFS data for the luminous infrared galaxy NGC7469: a nearby (70.6Mpc) active galaxy with a Sy 1.5 nucleus that drives a highly ionized gas outflow and a prominent nuclear star-forming ring. Using the superb sensitivity and high spatial resolution of the JWST instrument NIRSpec-IFS, we investigate the role of the Seyfert nucleus in the excitation and dynamics of the circumnuclear gas. Our analysis focuses on the [Fe ii], H2, and hydrogen recombination lines that trace the radiation/shocked-excited molecular and ionized ISM around the AGN. We investigate the gas excitation through H2/Br{\gamma} and [Fe ii]/Pa\b{eta} emission line ratios and find that photoionization by the AGN dominates within the central 300 pc of the galaxy and together with a small region show ing signatures of shock-heated gas; these shock-heated regions are likely associated with a compact radio jet. In addition, the velocity field and velocity dispersion maps reveal complex gas kinematics. Rotation is the dominant feature, but we also identify non-circular motions consistent with gas inflows as traced by the velocity residuals and the spiral pattern in the Pa{\alpha} velocity dispersion map. The inflow is consistent with the mass outflow rate and two orders of magnitude higher than the AGN accretion rate. The compact nuclear radio jet has enough power to drive the highly ionized outflow. This scenario suggests that the inflow and outflow are in a self-regulating feeding-feedback process, with a contribution from the radio jet helping to drive the outflow.

Mid-infrared emission features probe the properties of ionized gas, and hot or warm molecular gas. The Orion Bar is a frequently studied photodissociation region (PDR) containing large amounts of gas under these conditions, and was observed with the MIRI IFU aboard JWST as part of the "PDRs4All" program. The resulting IR spectroscopic images of high angular resolution (0.2") reveal a rich observational inventory of mid-IR emission lines, and spatially resolve the substructure of the PDR, with a mosaic cutting perpendicularly across the ionization front and three dissociation fronts. We extracted five spectra that represent the ionized, atomic, and molecular gas layers, and measured the most prominent gas emission lines. An initial analysis summarizes the physical conditions of the gas and the potential of these data. We identified around 100 lines, report an additional 18 lines that remain unidentified, and measured the line intensities and central wavelengths. The H I recombination lines originating from the ionized gas layer bordering the PDR, have intensity ratios that are well matched by emissivity coefficients from H recombination theory, but deviate up to 10% due contamination by He I lines. We report the observed emission lines of various ionization stages of Ne, P, S, Cl, Ar, Fe, and Ni, and show how certain line ratios vary between the five regions. We observe the pure-rotational H$_2$ lines in the vibrational ground state from 0-0 S(1) to 0-0 S(8), and in the first vibrationally excited state from 1-1 S(5) to 1-1 S(9). We derive H$_2$ excitation diagrams, and approximate the excitation with one thermal (~700 K) component representative of an average gas temperature, and one non-thermal component (~2700 K) probing the effect of UV pumping. We compare these results to an existing model for the Orion Bar PDR and highlight the differences with the observations.

We use deep Spitzer mid-infrared spectroscopic maps of radial strips across three nearby galaxies with well-studied metallicity gradients (M101, NGC 628, and NGC 2403) to explore the physical origins of the observed deficit of polycyclic aromatic hydrocarbons (PAHs) at subsolar metallicity (i.e., the PAH-metallicity relation or PZR). These maps allow us to trace the evolution of all PAH features from 5-18 mu m as metallicity decreases continuously from solar (Z circle dot) to 0.2 Z circle dot. The total PAH-to-dust luminosity ratio remains relatively constant until reaching a threshold of similar to 2/3 Z circle dot, below which it declines smoothly but rapidly. The PZR has been attributed to preferential destruction of the smallest grains in the hard radiation environments found at low metallicity. In this scenario, a decrease in emission from the shortest-wavelength PAH features is expected. In contrast, we find a steep decline in long-wavelength power below Z circle dot, especially in the 17 mu m feature, with the shorter-wavelength PAH bands carrying an increasingly large fraction of power at low metallicity. We use newly developed grain models to reproduce the observed PZR trends, including these variations in fractional PAH feature strengths. The model that best reproduces the data employs an evolving grain size distribution that shifts to smaller sizes as metallicity declines. We interpret this as a result of inhibited grain growth at low metallicity, suggesting continuous replenishment in the interstellar medium is the dominant process shaping the PAH grain population in galaxies.

The JWST has captured the most detailed and sharpest infrared 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). We investigate the fundamental interaction of far-ultraviolet 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. 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. A complex, structured, and folded DF surface was traced by the H2 lines. This dataset was used to revisit the commonly adopted 2D PDR structure of the Orion Bar. 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.

We present a focused study of radially-resolved varying PAH emission in the low-luminosity AGN-host NGC 4138 using deep Spitzer/IRS spectral maps. Using new model PAH spectra, we investigate whether these variations could be associated with changes to the PAH grain size distribution due to photodestruction by the AGN. Separately, we model the effects of the varying radiation field within NGC 4138, and we use this model to predict the corresponding changes in the PAH emission spectrum. We find that PAH band ratios are strongly variable with radius in this galaxy with short-to-long wavelength band ratios peaking in the starburst ring. The changing mix of starlight appears to have a considerable effect on the trends in these band ratios, and our radiation model predicts the shapes of these trends. However, the amplitude of observed variation is  2.5 times larger than predicted for some ratios. A cutoff of small grains in the PAH size distribution, as has been suggested for AGN, together with changes in PAH ionization fraction could explain the behavior of the shorter bands, but this model fails to reproduce longer band behaviors. Additionally, we find that short-to-long wavelength PAH band ratios increase slightly within  270pc of the center, suggesting that the AGN may directly influence PAH emission there.

(Abridged) We investigate the impact of radiative feedback from massive stars on their natal cloud and focus on the transition from the HII region to the atomic PDR (crossing the ionisation front (IF)), and the subsequent transition to the molecular PDR (crossing the dissociation front (DF)). We use high-resolution near-IR integral field spectroscopic data from NIRSpec on JWST to observe the Orion Bar PDR as part of the PDRs4All JWST Early Release Science Program. The NIRSpec data reveal a forest of lines including, but not limited to, HeI, HI, and CI recombination lines, ionic lines, OI and NI fluorescence lines, Aromatic Infrared Bands (AIBs including aromatic CH, aliphatic CH, and their CD counterparts), CO2 ice, pure rotational and ro-vibrational lines from H2, and ro-vibrational lines HD, CO, and CH+, most of them detected for the first time towards a PDR. Their spatial distribution resolves the H and He ionisation structure in the Huygens region, gives insight into the geometry of the Bar, and confirms the large-scale stratification of PDRs. We observe numerous smaller scale structures whose typical size decreases with distance from Ori C and IR lines from CI, if solely arising from radiative recombination and cascade, reveal very high gas temperatures consistent with the hot irradiated surface of small-scale dense clumps deep inside the PDR. The H2 lines reveal multiple, prominent filaments which exhibit different characteristics. This leaves the impression of a "terraced" transition from the predominantly atomic surface region to the CO-rich molecular zone deeper in. This study showcases the discovery space created by JWST to further our understanding of the impact radiation from young stars has on their natal molecular cloud and proto-planetary disk, which touches on star- and planet formation as well as galaxy evolution.

Aromatic infrared bands (AIBs) are a set of broad emission bands at 3.3, 6.2, 7.7, 8.6, 11.2, and 12.7 $\mu$m, seen in the infrared spectra of most galaxies. With JWST, the 3.3 $\mu$m AIB can in principle be detected up to a redshift of $\sim$ 7. Relating the evolution of the 3.3 $\mu$m AIB to local physical properties of the ISM is thus of paramount importance. By applying a dedicated machine learning algorithm to JWST NIRSpec observations of the Orion Bar photodissociation region obtained as part of the PDRs4All Early Release Science (ERS) program, we extracted two template spectra capturing the evolution of the AIB-related emission in the 3.2-3.6 $\mu$m range, which includes the AIB at 3.3 $\mu$m and its main satellite band at 3.4 $\mu$m. In the Orion Bar, we analyze the spatial distribution of the templates and their relationship with the fluorescent emission of H$_2$ in the near infrared. We find that one template ("AIB$_{\rm Irrad}$") traces regions of neutral atomic gas with strong far-UV fields, while the other template ("AIB$_{\rm Shielded}$") corresponds to shielded regions with lower FUV fields and a higher molecular gas fraction. We then show that these two templates can be used to fit the NIRSpec AIB-related spectra of nearby galaxies. The relative weight of the two templates (AIB$_{\rm Irrad/Shielded}$) is a tracer of the radiative feedback from massive stars on the ISM. We derive an estimate of AIB$_{\rm Irrad/Shielded}$ in a $z$ = 4.22 lensed galaxy, and find that it has a lower value than for local galaxies. This pilot study illustrates how a detailed analysis of AIB emission in nearby regions can be used to probe the physical conditions of the extragalactic ISM.

We present James Webb Space Telescope (JWST) integral-field spectrograph observations of NGC 3256, a local infrared-luminous late-stage merging system with two nuclei roughly 1 kpc apart, both of which have evidence of cold molecular outflows. Using JWST/NIRSpec and Mid-Infrared Instrument data sets, we investigate this morphologically complex system on spatial scales of <100 pc, where we focus on the warm molecular H-2 gas surrounding the nuclei. We detect collimated outflowing warm H-2 gas originating from the southern nucleus, though we do not find significant outflowing H-2 gas surrounding the northern nucleus. We measure maximum intrinsic outflow velocities of similar to 1000 km s(-1), which extend out to a distance of 0.7 kpc. Based on H-2 S(7)/S(1) ratios, we find a larger fraction of warmer gas near the southern nucleus, which decreases with increasing distance from the nucleus, signifying the southern nucleus as a primary source of H-2 heating. The gas mass of the warm H-2 outflow component is estimated to be M (warm ,out) = (1.4 +/- 0.2) x 10(6) M (circle dot), as much as 6% of the cold H2 mass estimated using Atacama Large Millimeter/submillimeter Array CO data. The outflow timescale is about 7 x 10(5) yr, resulting in a mass outflow rate M-warm,M-out=2.0 +/- 0.8M(circle dot)yr(-1) and kinetic power P (warm,out) similar to 4 x 10(41) erg s(-1). Lastly, regions within our 3.'' 0 x 3.'' 0 NIRSpec data, where the outflowing gas reside show high [Fe II]/Pa beta and H-2/Br gamma line ratios, indicate enhanced mechanical heating caused by the outflows. The fluxes and ratios of polycyclic aromatic hydrocarbons (PAHs) in these regions are not significantly different compared to those elsewhere in the disk, suggesting the outflows may not significantly alter the PAH ionization state or grain size.

The enormous increase in mid-IR sensitivity and spatial and spectral resolution provided by the JWST spectrographs enables, for the first time, detailed extragalactic studies of molecular vibrational bands. This opens an entirely new window for the study of the molecular interstellar medium in luminous infrared galaxies (LIRGs). We present a detailed analysis of rovibrational bands of gas-phase CO, H _2 O, C _2 H _2 , and HCN toward the heavily obscured eastern nucleus of the LIRG VV 114, as observed by NIRSpec and the medium resolution spectrograph on the Mid-InfraRed Instrument (MIRI MRS). Spectra extracted from apertures of 130 pc in radius show a clear dichotomy between the obscured active galactic nucleus (AGN) and two intense starburst regions. We detect the 2.3 μ m CO bandheads, characteristic of cool stellar atmospheres, in the star-forming regions, but not toward the AGN. Surprisingly, at 4.7 μ m, we find highly excited CO ( T _ex ≈ 700–800 K out to at least rotational level J = 27) toward the star-forming regions, but only cooler gas ( T _ex ≈ 200 K) toward the AGN. We conclude that only mid-infrared pumping through the rovibrational lines can account for the equilibrium conditions found for CO and H _2 O in the deeply embedded starbursts. Here, the CO bands probe regions with an intense local radiation field inside dusty young massive star clusters or near the most massive young stars. The lack of high-excitation molecular gas toward the AGN is attributed to geometric dilution of the intense radiation from the bright point source. An overview of the relevant excitation and radiative transfer physics is provided in an appendix.

We present a focused study of radially resolved varying PAH emission in the low-luminosity active galactic nucleus (AGN)-host NGC 4138 using deep Spitzer/infrared spectrograph spectral maps. Using new model PAH spectra, we investigate whether these variations could be associated with changes to the PAH grain size distribution due to photodestruction by the AGN. Separately, we model the effects of the varying radiation field within NGC 4138, and we use this model to predict the corresponding changes in the PAH emission spectrum. We find that PAH band ratios are strongly variable with radius in this galaxy with short-to-long wavelength band ratios peaking in the starburst ring. The changing mix of starlight appears to have a considerable effect on the trends in these band ratios, and our radiation model predicts the shapes of these trends. However, the amplitude of observed variation is ∼2.5 × larger than predicted for some ratios. A cutoff of small grains in the PAH size distribution, as has been suggested for AGN, together with changes in PAH ionization fraction could explain the behavior of the shorter bands, but this model fails to reproduce longer band behaviors. Additionally, we find that short-to-long wavelength PAH band ratios increase slightly within ∼270 pc of the center, suggesting that the AGN may directly influence PAH emission there.

We introduce a novel model to spectroscopically constrain the mid-infrared (MIR) extinction/attenuation curve from 3--17 um, using Polycyclic Aromatic Hydrocarbon (PAH) emission drawn from an AKARI-Spitzer extragalactic cross-archival dataset. Currently proposed MIR extinction curves vary significantly in their slopes toward the near-infrared, and the variation of the strengths and shapes of the 9.7 um and 18 um silicate absorption features make MIR spectral modeling and interpretation challenging, particularly for heavily obscured galaxies. By adopting the basic premise that PAH bands have relatively consistent intrinsic ratios within dusty starbursting galaxies, we can, for the first time, empirically determine the overall shape of the MIR attenuation curve by measuring the differential attenuation at specific PAH wavelengths. Our attenuation model shows PAH emission in most (U)LIRGs is unambiguously subjected to attenuation, and we find strong evidence that PAH bands undergo differential attenuation as obscuration increases. Compared to pre-existing results, the MIR attenuation curve derived from the model favors relatively gray continuum absorption from 3-8 $\mu$m and silicate features with intermediate strength at 9.7 um but with stronger than typical 18 um opacity.

We present new observations of the central 1 kpc of the M82 starburst obtained with the James Webb Space Telescope near-infrared camera instrument at a resolution θ ∼ 0.″05–0.″1 (∼1–2 pc). The data comprises images in three mostly continuum filters (F140M, F250M, and F360M), and filters that contain [Fe ii ] (F164N), H _2 v = 1 → 0 (F212N), and the 3.3 μ m polycyclic aromatic hydrocarbon (PAH) feature (F335M). We find prominent plumes of PAH emission extending outward from the central starburst region, together with a network of complex filamentary substructures and edge-brightened bubble-like features. The structure of the PAH emission closely resembles that of the ionized gas, as revealed in Paschen α and free–free radio emission. We discuss the origin of the structure, and suggest the PAHs are embedded in a combination of neutral, molecular, and photoionized gas.

(Abridged) Mid-infrared observations of photodissociation regions (PDRs) are dominated by strong emission features called aromatic infrared bands (AIBs). The most prominent AIBs are found at 3.3, 6.2, 7.7, 8.6, and 11.2 $\mu$m. The most sensitive, highest-resolution infrared spectral imaging data ever taken of the prototypical PDR, the Orion Bar, have been captured by JWST. We provide an inventory of the AIBs found in the Orion Bar, along with mid-IR template spectra from five distinct regions in the Bar: the molecular PDR, the atomic PDR, and the HII region. We use JWST NIRSpec IFU and MIRI MRS observations of the Orion Bar from the JWST Early Release Science Program, PDRs4All (ID: 1288). We extract five template spectra to represent the morphology and environment of the Orion Bar PDR. The superb sensitivity and the spectral and spatial resolution of these JWST observations reveal many details of the AIB emission and enable an improved characterization of their detailed profile shapes and sub-components. While the spectra are dominated by the well-known AIBs at 3.3, 6.2, 7.7, 8.6, 11.2, and 12.7 $\mu$m, a wealth of weaker features and sub-components are present. We report trends in the widths and relative strengths of AIBs across the five template spectra. These trends yield valuable insight into the photochemical evolution of PAHs, such as the evolution responsible for the shift of 11.2 $\mu$m AIB emission from class B$_{11.2}$ in the molecular PDR to class A$_{11.2}$ in the PDR surface layers. This photochemical evolution is driven by the increased importance of FUV processing in the PDR surface layers, resulting in a "weeding out" of the weakest links of the PAH family in these layers. For now, these JWST observations are consistent with a model in which the underlying PAH family is composed of a few species: the so-called 'grandPAHs'.

Context.Extraordinary observations of the Orion Bar by JWST have shown, for the first time, the incredible richness of polycyclic aromatic hydrocarbon (PAH) emission bands and their variation on very small scales. These variations are the result of photochemical evolution of the PAH carrier.Aims.We aim to probe the photochemical evolution of PAHs across the key zones of the ideal photodissociation region (PDR) that is the Orion Bar using unsupervised machine learning.Methods.We used JWST NIRSpec IFU and MIRI MRS observations of the Orion Bar from the JWST Early Release Science programme PDRs4All (ID: 1288). We levered bisectingk-means clustering to generate highly detailed spatial maps of the spectral variability in the 3.2–3.6, 5.95–6.6, 7.25–8.95, and 10.9–11.63 μm wavelength regions. We analysed and subsequently described the variations in the cluster profiles and connected them to the conditions of the physical locations from which they arise. We interpreted the origin of the observed variations with respect to the following key zones: the HIIregion, the atomic PDR zone, and the layers of the molecular PDR zone stratified by the first, second, and third dissociation fronts (DF 1, DF 2, and DF 3, respectively).Results.Observed PAH emission exhibits spectral variation that is highly dependent on the spatial position in the PDR. We find the 8.6 μm band to behave differently than all other bands, which vary systematically with one another. Notably, we find a uniform variation in the 3.4–3.6 μm bands and 3.4/3.3 intensity ratio. We attribute the carrier of the 3.4–3.6 μm bands to a single side group attached to very similarly sized PAHs. Further, cluster profiles reveal a transition between characteristic profile classes of the 11.2 μm feature from the atomic to the molecular PDR zones. We find the carriers of each of the profile classes to be independent, and reason the latter to be PAH clusters existing solely deep in the molecular PDR. Clustering also reveals a connection between the 11 .2 and 6.2 μm bands and that clusters generated from variation in the 10.9–11.63 μm region can be used to recover those in the 5.95–6.6 μm region.Conclusions.Clustering is a powerful and comprehensive tool for characterising PAH spectral variability on both spatial and spectral scales. For individual bands as well as global spectral behaviours, we find ultraviolet processing to be the most important driver of the evolution of PAHs and their spectral signatures in the Orion Bar PDR.

The early -type star gamma Cas illuminates the reflection nebulae IC 59 and IC 63, creating two photodissociation regions (PDRs). Uncertainties about the distances to the nebulae and the resulting uncertainty about the density of the radiation fields incident on their surfaces have hampered the study of these PDRs during the past three decades. We employed far -ultra violet (UV) - optical nebula - star colour differences of dust -scattered light to infer the locations of the nebulae relative to the plane of the sky containing gamma Cas, finding IC 63 to be positioned behind the star and IC 59 in front of the star. To obtain the linear distances of the nebulae relative to gamma Cas, we fit far -infrared archi v al Herschel flux data for IC 59 and IC 63 with modified blackbody curves and relate the resulting dust temperatures with the luminosity of gamma Cas, yielding approximate distances of 4.15 pc for IC 59 and 2.3 pc for IC 63. With these distances, using updated far -UV flux data in the 6-13.6 eV range for gamma Cas with two recent determinations of the interstellar extinction for gamma Cas, we estimate that the far -UV radiation density at the surface of IC 63 takes on values of G(0) = 58 or G(0) = 38 with respecti ve v alues for E(B - V ) for gamma Cas of 0.08 and 0.04 mag. This is a substantial reduction from the range 150 <= G(0 )<= 650 used for IC 63 during the past three decades. The corresponding, e ven lo wer ne w v alues for IC 59 are G(0) = 18 and G(0) = 12.