Context. Midinfrared 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 m. The most sensitive highestresolution infrared spectral imaging data ever taken of the prototypical PDR the Orion Bar have been captured by JWST. These highquality data allow for an unprecedentedly detailed view of AIBs. BR  Aims We provide an inventory of the AIBs found in the Orion Bar along with midIR template spectra from five distinct regions in the Bar the molecular PDR i.e. the three HSUB2SUB dissociation fronts the atomic PDR and the H II region. BR  Methods We used JWST NIRSpec IFU and MIRI MRS observations of the Orion Bar from the JWST Early Release Science Program PDRs4All ID 1288. We extracted five template spectra to represent the morphology and environment of the Orion Bar PDR. We investigated and characterised the AIBs in these template spectra. We describe the variations among them here. BR  Results 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 subcomponents. The Orion Bar spectra are dominated by the wellknown AIBs at 3.3 6.2 7.7 8.6 11.2 and 12.7 m with welldefined profiles. In addition the spectra display a wealth of weaker features and subcomponents. The widths of many AIBs show clear and systematic variations being narrowest in the atomic PDR template but showing a clear broadening in the H II region template while the broadest bands are found in the three dissociation front templates. In addition the relative strengths of AIB subcomponents vary among the template spectra as well. All AIB profiles are characteristic of class A sources as designated by Peeters 2022 AampA 390 1089 except for the 11.2 m AIB profile deep in the molecular zone which belongs to class BSUB11.2SUB. Furthermore the observations show that the subcomponents that contribute to the 5.75 7.7 and 11.2 m AIBs become much weaker in the PDR surface layers. We attribute this to the presence of small more labile carriers in the deeper PDR layers that are photolysed away in the harsh radiation field near the surface. The 3.311.2 AIB intensity ratio decreases by about 40 between the dissociation fronts and the H II region indicating a shift in the polycyclic aromatic hydrocarbon PAH size distribution to larger PAHs in the PDR surface layers also likely due to the effects of photochemistry. The observed broadening of the bands in the molecular PDR is consistent with an enhanced importance of smaller PAHs since smaller PAHs attain a higher internal excitation energy at a fixed photon energy. BR  Conclusions Spectralimaging observations of the Orion Bar using JWST yield key insights into the photochemical evolution of PAHs such as the evolution responsible for the shift of 11.2 m AIB emission from class BSUB11.2SUB in the molecular PDR to class ASUB11.2SUB 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 socalled grandPAHs. P The 5 template spectra are available at the CDS via anonymous ftp to A hrefhttpscdsarc.cds.unistra.frcdsarc.cds.unistra.frA ftp130.79.128.5 or via A hrefhttpscdsarc.cds.unistra.frvizbincatJAA685A75httpscdsarc.cds.unistra.frvizbincatJAA685A75A P Tim Lee sadly passed away on Nov. 3 2022.
['astrochemistry', 'infrared: ISM', 'ISM: molecules', 'ISM: individual objects: Orion Bar', 'photon-dominated region (PDR)', 'techniques: spectroscopic', 'Astrophysics - Astrophysics of Galaxies', 'Astrophysics - Solar and Stellar Astrophysics']
PDRs4All. IV. An embarrassment of riches Aromatic infrared bands in the Orion Bar
{'Aromatic infrared bands in the Orion Bar': "Ryan Chown 1 , 2 , Ameek Sidhu 1 , 2 , Els Peeters 1 , 2 , 3 , Alexander G. G. M. Tielens 4 , 5 , Jan Cami 1 , 2 , 3 , Olivier Berné 6 , Emilie Habart 7 , Felipe Alarcón 8 , Amélie Canin 6 , Ilane Schroetter 6 , Boris Trahin 7 , Dries Van De Putte 9 , Alain Abergel 7 , Edwin A. Bergin 8 , Jeronimo Bernard-Salas 10 , 11 , Christiaan Boersma 12 , Emeric Bron 13 , Sara Cuadrado 14 , Emmanuel Dartois 15 , Daniel Dicken 16 , Meriem El-Yajouri 7 , Asunción Fuente 17 , Javier R. Goicoechea 14 , Karl D. Gordon 9 , 18 , Lina Issa 6 , Christine Joblin 4 , Olga Kannavou 7 , Baria Khan 1 , Ozan Lacinbala 19 , David Languignon 13 , Romane Le Gal 20 , 21 , Alexandros Maragkoudakis 12 , Raphael Meshaka 7 , 13 , Yoko Okada 22 , Takashi Onaka 23 , 24 , Sofia Pasquini 1 , Marc W. Pound 5 , Massimo Robberto 9 , 25 , Markus Röllig 26 , 27 , Bethany Schefter 1 , 2 , Thiébaut Schirmer 7 , 28 , Sílvia Vicente 29 , Mark G. Wolfire 5 , Marion Zannese 7 , Isabel Aleman 30 , Louis Allamandola 12 , 31 , Rebecca Auchettl 32 , Giuseppe Antonio Baratta 33 , Salma Bejaoui 12 , Partha P. Bera 12 , 29 , John H. Black 28 , Francois Boulanger 34 , Jordy Bouwman 35 , 36 , 37 , Bernhard Brandl 4 , 38 , Philippe Brechignac 15 , Sandra Brünken 39 , Mridusmita Buragohain 40 , Andrew Burkhardt 41 , Alessandra Candian 42 , Stéphanie Cazaux 43 , Jose Cernicharo 14 , Marin Chabot 44 , Shubhadip Chakraborty 45 , 46 , Jason Champion 6 , Sean W.J. Colgan 12 , Ilsa R. Cooke 47 , Audrey Coutens 6 , Nick L.J. Cox 10 , 11 , Karine Demyk 6 , Jennifer Donovan Meyer 48 , Sacha Foschino 6 , Pedro García-Lario 49 , Lisseth Gavilan 12 , Maryvonne Gerin 50 , Carl A. Gottlieb 51 , Pierre Guillard 52 , 53 , Antoine Gusdorf 34 , 50 , Patrick Hartigan 54 , Jinhua He 55 , 56 , 57 , Eric Herbst 58 , Liv Hornekaer 59 , Cornelia Jäger 60 , Eduardo Janot-Pacheco 61 , Michael Kaufman 62 , Francisca Kemper 63 , 64 , 65 , Sarah Kendrew 66 , Maria S. Kirsanova 67 , Pamela Klaassen 16 , Sun Kwok 68 , Álvaro Labiano 69 , Thomas S.-Y. Lai 70 , Timothy J. Lee 12 ⋆ , Bertrand Lefloch 20 , Franck Le Petit 13 , Aigen Li 71 , Hendrik Linz 72 , Cameron J. Mackie 73 , 74 , Suzanne C. Madden 75 , Joëlle Mascetti 76 , Brett A. McGuire 48 , 77 , Pablo Merino 78 , Elisabetta R. Micelotta 79 , Karl Misselt 80 , Jon A. Morse 81 , Giacomo Mulas 33 , 6 , Naslim Neelamkodan 82 , Ryou Ohsawa 83 , Alain Omont 84 , Roberta Paladini 85 , Maria Elisabetta Palumbo 33 , Amit Pathak 86 , Yvonne J. Pendleton 87 , Annemieke Petrignani 88 , Thomas Pino 15 , Elena Puga 66 , Naseem Rangwala 12 , Mathias Rapacioli 89 , Alessandra Ricca 12 , 3 , Julia Roman-Duval 9 , Joseph Roser 3 , 12 , Evelyne Roue ff 13 , Gaël Rouillé 90 , Farid Salama 12 , Dinalva A. Sales 91 , Karin Sandstrom 92 , Peter Sarre 93 , Ella Sciamma-O'Brien 12 , Kris Sellgren 94 , Sachindev S. Shenoy 95 , David Teyssier 69 , Richard D. Thomas 96 , Aditya Togi 72 , Laurent Verstraete 7 , Adolf N. Witt 97 , Alwyn Wootten 48 , Henning Zettergren 96 , Yong Zhang 98 , Ziwei E. Zhang 99 , and Junfeng Zhen 100  \n(A ffi liations can be found after the references)  \nReceived 15 April 2023; accepted 1 August 2023", 'ABSTRACT': "Context. 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 µ m. The most sensitive, highest-resolution infrared spectral imaging data ever taken of the prototypical PDR, the Orion Bar, have been captured by JWST . These high-quality data allow for an unprecedentedly detailed view of AIBs.  \nAims. 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 (i.e. the three H2 dissociation fronts), the atomic PDR, and the H ii region.  \nMethods. We used JWST NIRSpec IFU and MIRI MRS observations of the Orion Bar from the JWST Early Release Science Program, PDRs4All (ID: 1288). We extracted five template spectra to represent the morphology and environment of the Orion Bar PDR. We investigated and characterised the AIBs in these template spectra. We describe the variations among them here.  \nResults. 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. The Orion Bar spectra are dominated by the well-known AIBs at 3.3, 6.2, 7.7, 8.6, 11.2, and 12.7 µ m with well-defined profiles. In addition, the spectra display a wealth of weaker features and subcomponents. The widths of many AIBs show clear and systematic variations, being narrowest in the atomic PDR template, but showing a clear broadening in the H ii region template while the broadest bands are found in the three dissociation front templates. In addition, the relative strengths of AIB (sub-)components vary among the template spectra as well. All AIB profiles are characteristic of class A sources as designated by Peeters et al. (2002a), except for the 11.2 µ m AIB profile deep in the molecular zone, which belongs to class B11 . 2. Furthermore, the observations show that the sub-components that contribute to the 5.75, 7.7, and 11.2 µ m AIBs become much weaker in the PDR surface layers. We attribute this to the presence of small, more labile carriers in the deeper PDR layers that are photolysed away in the harsh radiation field near the surface. The 3.3 / 11.2 AIB intensity ratio decreases by about 40% between the dissociation fronts and the H ii region, indicating a shift in the polycyclic aromatic hydrocarbon (PAH) size distribution to larger PAHs in the PDR surface layers, also likely due to the e ff ects of photochemistry. The  \nobserved broadening of the bands in the molecular PDR is consistent with an enhanced importance of smaller PAHs since smaller PAHs attain a higher internal excitation energy at a fixed photon energy.  \nConclusions. Spectral-imaging observations of the Orion Bar using JWST yield key insights into the photochemical evolution of PAHs, such as the evolution responsible for the shift of 11.2 µ m AIB emission from class B11 . 2 in the molecular PDR to class A11 . 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'.  \nKey words. astrochemistry - infrared: ISM - ISM: molecules - ISM: individual objects: Orion Bar - ISM: photon-dominated region (PDR) techniques: spectroscopic", '1. Introduction': "A major component of the infrared (IR) emission near starforming regions in the Universe consists of a set of broad emission features at 3.3, 6.2, 7.7, 8.6, 11.2, and 12.7 µ m (e.g. Tielens 2008, and references therein). These mid-IR emission features, referred to as aromatic infrared bands (AIBs), are generally attributed to vibrational emission from polycyclic aromatic hydrocarbons and related species upon absorption of interstellar far-ultraviolet (FUV; 6-13.6 eV) photons (Léger & Puget 1984; Allamandola et al. 1985). The AIB spectrum is very rich and consists of the main bands listed above and a plethora of weaker emission features. Moreover, many AIBs are in fact blends of strong and weak bands (e.g. Peeters et al. 2004a). The AIB emission is known to vary from source to source and spatially within extended sources in terms of the profile and relative intensities of the features (e.g. Joblin et al. 1996; Hony et al. 2001; Berné et al. 2007; Sandstrom et al. 2010; Boersma et al. 2012; Candian et al. 2012; Stock & Peeters 2017; Peeters et al. 2017). These remarkably widespread emission features have been described in many diverse astronomical sources, including protoplanetary disks (e.g. Meeus et al. 2001; Vicente et al. 2013), H ii regions (e.g. Bregman 1989; Peeters et al. 2002b), reflection nebulae (e.g. Peeters et al. 2002a; Werner et al. 2004), planetary nebulae (e.g. Gillett et al. 1973; Bregman 1989; Beintema et al. 1996), the interstellar medium (ISM) of galaxies ranging from the Milky Way (Boulanger et al. 1996), the Magellanic Clouds (e.g. Vermeij et al. 2002; Sandstrom et al. 2010), starburst galaxies, luminous and ultra-luminous IR galaxies, and high-redshift galaxies (e.g. Genzel et al. 1998; Lutz et al. 1998; Peeters et al. 2004b; Yan et al. 2005; Galliano et al. 2008), as well as in the harsh environments of galactic nuclei (e.g. Smith et al. 2007; Esquej et al. 2014; Jensen et al. 2017).  \nA useful observational proxy for studying AIBs is the spectroscopic classification scheme devised by Peeters et al. (2002a) which classifies each individual AIBs based on their profile shapes and precise peak positions (classes A, B, and C). While the AIBs observed in a given source generally belong to the same class, this is not always the case. In particular, the classes in the 6 to 9 µ m region do not always correspond to those of the 3.3 and 11.2 AIBs (van Diedenhoven et al. 2004). Class A sources are the most common - they exhibit the 'classical' AIBs, with a 6.2 µ mAIBthat peaks between 6.19 and 6.23 µ m, a 7.7 µ mcomplex in which the 7.6 µ m sub-peak is stronger than the 7.8 µ m sub-peak, and the 8.6 µ m feature peaks at 8.6 µ m. Class B sources can be slightly redshifted compared to class A, while at the same time the 7.7 µ mcomplex peaks between 7.8 and 8 µ m. Class C sources show a very broad emission band peaking near 8.2 µ m, and typically do not exhibit the 6.2 or 7.7 µ mAIBs.  \nThese three classes were found to show a strong correlation with the type of object considered. The most common AIB spec-  \ntrum, class A, is identified in the spectra of photodissociation regions (PDRs), H ii regions, reflection nebulae, the ISM, and galaxies. The most widely used template for class A sources has been the spectrum of the Orion Bar (Peeters et al. 2002a; van Diedenhoven et al. 2004). Class B sources are isolated Herbig Ae / Be stars and a few evolved stars; in fact, evolved star spectra can belong to either of the classes. Class C sources include post-AGB and Herbig Ae / Be stars, as well as a few T-Tauri disks (Peeters et al. 2002a; Bouwman et al. 2008; Shannon & Boersma 2019). More recent work has developed analogous classification schemes for other AIBs and has included a new class D (e.g. van Diedenhoven et al. 2004; Sloan et al. 2014; Matsuura et al. 2014).  \nObserved variations in AIBs reflect changes in the molecular properties of the species responsible for the AIB emission (charge, size, and molecular structure; e.g. Joblin et al. 1996; Berné et al. 2007; Pilleri et al. 2012; Boersma et al. 2013; Candian & Sarre 2015; Peeters et al. 2017; Robertson 1986; Dartois et al. 2004; Pino et al. 2008; Godard et al. 2011; Jones et al. 2013), which are set by the local physical conditions (including FUV radiation field strength, G 0, gas temperature, and density, n (H); e.g. Bakes et al. 2001; Galliano et al. 2008; Pilleri et al. 2012, 2015; Stock et al. 2016; Schirmer et al. 2020; Schirmer et al. 2022; Sidhu et al. 2022; Knight et al. 2022b; Murga et al. 2022). The observed variability in AIB emission thus implies that the population responsible for their emission is not static, but undergoes photochemical evolution.  \nObservations using space-based IR observatories - in particular the Short-Wavelength Spectrometer (SWS; de Graauw et al. 1996) on board the Infrared Space Observatory ( ISO ; Kessler et al. 1996) and the Infrared Spectrograph (IRS; Houck et al. 2004) on board the Spitzer Space Telescope (Werner et al. 2004) - have revealed the richness of AIBs (for a review see e.g. Peeters et al. 2004a; Tielens 2008). However, obtaining a full understanding of the photochemical evolution underlying AIBs has been limited by insu ffi cient spatial and spectral resolution (e.g. Spitzer-IRS ) or by limited sensitivity and spatial resolution (e.g. ISO / SWS) of these IR facilities.  \nJWST is set to unravel the observed complexity of AIBs, as it o ff ers access to the full wavelength range of importance for AIB studies at medium spectral resolution and at unprecedented spatial resolution and sensitivity. JWST is able to resolve, for the first time, where and how the photochemical evolution of polycyclic aromatic related species, the carriers of AIBs, occurs while providing a detailed view of the resulting AIB spectral signatures. The PDRs4All JWST Early Release Science Program observed the prototypical highly irradiated PDR, the Orion Bar (Berné et al. 2022; Habart et al. 2023; Peeters et al. 2023). The Orion Bar PDR has a G0 which varies with position from about 1 × 10 4 to 4 × 10 4 Habings (e.g. Marconi et al. 1998; Peeters et al. 2023) and it has a gas density which varies from of a few 10 4 cm -3 in the atomic PDR to ∼ 10 6 cm -3 in the molecular region (e.g. Parmar et al. 1991; Tauber et al. 1994; Young Owl  \net al. 2000; Bernard-Salas et al. 2012a; Goicoechea et al. 2016; Joblin et al. 2018; Habart et al. 2023). Given the proximity of Orion (414 pc; Menten et al. 2007), the PDRs4All dataset takes full advantage of JWST 's spatial resolution to showcase the AIB emission in unprecedented detail.  \nIn this paper, we present five MIRI-MRS template spectra representing key regions of the Orion Bar PDR. Combined with corresponding JWST NIRSpec-IFU template spectra (Peeters et al. 2023), we present an updated inventory and characterization of the AIB emission in this important reference source. We describe the observations, data reduction, and the determination of the underlying continuum in our template spectra in Sect. 2. In Sect. 3, we give a detailed account of the observed AIB bands and sub-components along with their vibrational assignments. We compare our findings with previous works and discuss the AIB profiles and the AIB variability in the Orion Bar in Sect. 4. Finally, we summarize our results and narrate a picture of the origins and evolution of the AIB emission in Sect. 5.", '2.1. MIRI-MRS observations and data reduction': "On 30 January 2023, JWST observed the Orion Bar PDR with the Mid-Infrared Instrument (MIRI) in medium resolution spectroscopy (MRS) mode (Wells et al. 2015; Argyriou et al. 2023) as part of the PDRs4All Early Release Science program (Berné et al. 2022). We obtained a 1 × 9 pointing mosaic in all four MRS channels (channels 1, 2, 3, and 4), and all three sub-bands within each channel (short, medium, and long). We applied a 4-point dither optimised for extended sources and use the FASTR1 readout pattern adapted for bright sources. We integrated for 521.7 s using 47 groups per integration and 4 integrations. The resulting datacube thus spans the full MRS wavelength range (4.90 to 27.90 µ m) with a spectral resolution ranging from R ∼ 3700 in channel 1 to ∼ 1700 in channel 4 and a spatial resolution of 0.207 '' at short wavelengths to 0.803 '' at long wavelengths, corresponding to 86 and 332 AU, respectively at the distance of the Orion Nebula.  \nThe mosaic was positioned to overlap the PDRs4All JWST Near Infrared Spectrograph (NIRSpec) IFU (Böker et al. 2022) observations of the Orion Bar (Peeters et al. 2023, Fig. 1). Given the di ff erent fields of view of the MRS channels ( ∼ 3 '' in channel 1 to ∼ 7 '' in channel 4), the spatial footprint with full wavelength coverage is limited by the field of view (FOV) of channel 1. The footprint shown in Fig. 1 represents the area with full MRS wavelength coverage, noting that the MRS data in channels 2-4 exist beyond the area shown, but we choose to use only the sub-set of data with full wavelength coverage. The NIRSpec IFU and MIRI MRS datasets combined provide a perpendicular cross-section of the Orion Bar from the H ii region to the molecular zone at very high spatial and spectral resolution from 0.97 to 27.9 µ m.  \nWe reduced the MIRI-MRS data using version 1.9.5.dev10 + g04688a77 of the JWST pipeline 1 , and JWST Calibration Reference Data System 2 (CRDS) context 1041. We ran the JWST pipeline with default parameters except the following. The master background subtraction, outlier detection, fringe- and residual-fringe correction steps were all turned on. Cubes were built using the drizzle algorithm. The pipeline combined all pointings for each sub-band, resulting in 12 cubes  \n(4 channels of 3 sub-bands each) covering the entire field of view.  \nWe stitched the 12 sub-band cubes into a single cube by reprojecting all of the cubes onto a common spatial grid using channel 1 short as a reference. We then scaled the spectra to match in flux where they overlap using channel 2 long as the reference. This stitching algorithm is part of the 'Haute Couture' algorithm described in Canin et al., (in preparation).  \nWhile the pipeline and reference files produce high-quality data products, a few artefacts still remain in the data. The most important artifacts for our analysis are fringes and flux calibration that are not yet finalised (see Appendix B). Neither of these artefacts have a strong impact on our results, as discussed in Appendix B, although we do limit our analysis to wavelengths ≤ 15 µ mdue to the presence of artefacts beyond this range.", '2.2. Extracting template spectra from key regions': 'We extracted MIRI template spectra using the same extraction apertures as Peeters et al. (2023) 3 . These apertures are selected to represent the key physical zones of the Orion Bar PDR: the H ii region, the atomic PDR, and the three bright H i / H2 dissociation fronts (DF 1, DF 2, and DF 3) corresponding to three molecular hydrogen (H2) filaments that were identified in the NIRSpec FOV (Fig. 1). We emphasize that the remaining areas in the MRS spectral map will be analysed at a later time. We note that the AIB emission detected in the H ii region template originates from the background PDR. Combined with the NIRSpec templates of Peeters et al. (2023), these spectra capture all of the AIB emission in each of the five regions. In this paper, we focus on the inventory and characterization of the AIBs found in these template spectra. We refer to Peeters et al. (2023) and Van De Putte et al. (2023) for the inventories of the gas lines from atoms and small molecules extracted from NIRSpec and MIRI MRS data, respectively. For a detailed description of the Orion Bar PDR morphology as seen by JWST we refer to Habart et al. (2023) and to Peeters et al. (2023).', '2.3. Measuring the underlying continuum': 'The AIB emission is perched on top of the continuum emission from stochastically heated very small grains (e.g. Smith et al. 2007). Di ff erent spectral decomposition methods have deduced additional emission components referred to as: 1) emission from evaporating very small grains (based on the blind signal separation method; Berné et al. 2007; Pilleri et al. 2012; Foschino et al. 2019) and 2) plateau emission due to large PAHs, PAH clusters and nanoparticles (Bregman et al. 1989; Roche et al. 1989; Peeters et al. 2012; Boersma et al. 2014; Sloan et al. 2014; Peeters et al. 2017).  \nIn order to identify and characterize AIBs, we subtracted estimates of the continuum emission in each template spectrum. We computed a linear continuum for NIRSpec and a spline continuum anchored at selected wavelengths for MIRI data. Furthermore, we adopted the same anchor points for all five templates. While our measurements of the full-width at half-maximum (FWHM) of AIBs (see Table 1) do depend on the selected continuum, we note that our main goal - to catalog AIBs and their sub-components qualitatively - does not require highly precise estimates of the continuum.  \nFig. 1. PDRs4All MIRI MRS and NIRSpec footprints (dashed and solid white boundaries, respectively), and spectral extraction apertures (black boxes in the right panel) on top of a composite NIRCam image of the Orion Bar (data from Habart et al. 2023). DF 1, DF 2, and DF 3 are H2 dissociation fronts as designated in Habart et al. (2023). Red, green, and blue are encoded as F335M (AIB), F470N-F480M (H2 emission), and F187N (Paschen α ), respectively.  \n<!-- image -->  \nThe FWHM 4 of each AIB complex was measured by normalizing the continuum-subtracted template spectrum by the peak intensity of the AIB and then calculating the FWHM of the entire AIB complex, that is, without taking into consideration blends, components, and / or sub-components that make up the AIB complex. The measured FWHM strongly depends on the estimated continuum emission; however, this does not impact qualitative trends in FWHM from template to template. The integrated flux of each AIB was computed from the continuumsubtracted spectra. We refer to Peeters et al. (2023) for details on how the 3.3 and 3.4 µ mAIB fluxes were measured.', '3. AIB characteristics and assignments': "The superb quality of the Orion Bar observations combined with the increased spectral resolution compared to prior IR space observations reveals an ever-better characterization of the AIBs in terms of sub-components, multiple components making up a 'single' band, and the precise shapes of the band profiles (see Fig. 2 for an overview and Figs. 3 and 4 for zoom-ins of selected AIBs in the five template spectra).  \nWe o ff er a detailed description of the spectral characteristics of the AIB emission as seen by JWST as well as current vibrational assignments in Sects. 3.1 to 3.6. The detailed AIB inventory is listed in Table A.1. We note that we consider all spectrally resolved emission features to be candidate AIBs. To assess whether a candidate AIB is real or an artefact, we compare the template spectra in the location of the candidate AIB to the spectrum of the calibration standard star 10 Lac (see Appendix B for details). Due to the very high signal-to-noise ratio  \n(S / N), the spectra reveal an abundance of weak features, either as standalone features, or as shoulders of other bands. Occasionally these shoulders are only visible as a change in the slope along the wing of a stronger AIB and, in such cases, we estimated the central wavelength of the weak AIB visually based on the AIB profile of the main component.  \nHereafter, all mentions of nominal AIBs, given in boldface in col. 1 of Table A.1, namely 3.3, 3.4, 5.25, 5.75, 6.2, 7.7, 8.6, 11.2, 12.0, 12.7, and 13.5 µ m, do not indicate the precise peak positions of these AIBs. The precise peak positions of these nominal AIBs are reported in terms of wavelength in col. 3 of Table A.1. We note that we converted the positions in wavelength to wavenumber by rounding to the nearest integer in units of cm -1 and so the precision of the reported wavenumbers does not reflect the instrumental precision of the peak position of the AIBs.", '3.1. The 3.2-3.5 µ m (3125-2860 cm -1 ) range': "The 3 µ m spectral region is dominated by the 3.3 and 3.4 µ m AIBs that peak at 3.29 and 3.4 µ m, respectively. While some studies have found that the peak of the 3.3 µ mAIB shifts toward longer wavelengths (van Diedenhoven et al. 2004), our measurements of the peak position of this band in each template are consistent with the nominal value of 3.29 µ m. The band profiles, however, show some slight di ff erences among the templates. The width of the 3.29 µ mbandvaries slightly (see Table 1 and Fig. 5). The templates in order of increasing 3.29 µ m band width are: atomic PDR, DF 1, H ii region, DF 2, and DF 3. While the small increase in width on the red side may be attributed to underlying broad plateau emission (see Peeters et al. 2023), the blue wing broadens by ∼ 5 cm -1 on a total band width of ∼ 37 . 5 cm -1 (Table 1 and Peeters et al. 2023).  \nFig. 2. AIB spectrum as seen by JWST using the Orion Bar atomic PDR template spectrum (Sect. 2.2) as an example. Red shaded regions indicate emission from AIBs while blue curves indicate the underlying continuum. Figure is adapted from Peeters et al. (2004a).  \n<!-- image -->  \nThe 3.29 µ m band is characteristic for the CH stretching mode in PAHs. The peak position is somewhat dependent on molecular structure, for example the number of adjacent hydrogens on a ring and steric hindrance between opposing hydrogens in so-called 'bay' regions. Molecular symmetry has a more important e ff ect as it controls the number of allowed transitions and the range over which IR activity is present (Maltseva et al. 2015, 2016). For a given PAH, the initial excitation energy has a very minor influence on the peak position (Mackie et al. 2022). Earlier works have also demonstrated this minor influence (Joblin et al. 1995; Pech et al. 2002). Additionally, Tokunaga & Bernstein (2021) found that the peak position and width of the 3.3 µ m feature must be fitted simultaneously as both depend on the carrier. There is also a weak dependence of peak position on the charge state, but since the CH stretch is very weak in cations (Allamandola et al. 1999; Peeters et al. 2002a), this is of no consequence. These modes are very much influenced by resonance e ff ects with combination bands involving CC modes and CHin-plane bending modes (Mackie et al. 2015, 2016). Overall, in the emission spectra of highly excited species, the di ff erences in peak position mentioned here will be too subtle compared to the impact of molecular symmetry when attempting to identify the carrier(s). The observed narrow width of the 3.3 µ mAIB implies then emission by very symmetric PAHs (Pech et al. 2002; Ricca et al. 2012; Mackie et al. 2022).  \nThe very weak shoulder on the blue side, namely, at ≃ 3 . 246 µ m, has the same strength relative to the main feature in all template spectra, suggesting it is part of the same emission complex. Its peak wavelength may point toward the stretching mode of aromatic CH groups in bay regions or, alternatively, the e ff ect of resonant interaction in a specific species (van Diedenhoven et al. 2004; Candian et al. 2012; Mackie et al. 2015) or aromatic  \nCH in polyaromatic carbon clusters (Dubosq et al. 2023). In a recent analysis of the 3.3 µ m AIB in the Red Rectangle, Tokunaga et al. (2022) found di ff erences in the spectra of this source compared to earlier analyses (e.g. Tokunaga et al. 1991; Candian et al. 2012) due to the treatment of Pfund emission lines from the standard star. Candian et al. (2012) fit the 3.3 µ m AIB in each spaxel of their IFU cube with two components and analysed spatial variations in the integrated intensities of these components. While issues with the standard star spectrum would a ff ect all spectra in the cube (Tokunaga et al. 2022), spatial variations in integrated intensities should not be a ff ected.  \nThe AIB spectra reveal a plethora of bands longward of the 3.29 µ m feature between ≃ 3 . 4 and ≃ 3 . 6 µ m (Table A.1; Peeters et al. 2023; Sloan et al. 1997). As the relative strengths of these sub-components show variations from source to source and within sources (e.g. Joblin et al. 1996; Pilleri et al. 2015; Peeters et al. 2023), they are generally ascribed to di ff erent emitting groups on PAHs. Here, we note that the emission profile of the 3.4 µ m band varies between the five templates, broadening to longer wavelength (Peeters et al. 2023), indicating the presence of multiple components in the main 3.4 µ m band at 3.395, 3.403, and 3.424 µ m. The other bands do not show such profile variations. Bands in this wavelength range are due to the CH stretching mode in aliphatic groups and assignments to methyl (CH3) groups attached to PAHs and to superhydrogenated PAHs have been proposed (Joblin et al. 1996; Bernstein et al. 1996; Maltseva et al. 2018; Buragohain et al. 2020; Pla et al. 2020). As for the aromatic CH stretching mode, the peak position is sensitive to resonances with combination bands of CC modes and CH in-plane bending modes (Mackie et al. 2018). Typically, methylated PAHs show a prominent band around 3.4 µ m, but its peak position falls within a wide range, ≃ 0 . 17 µ m(Maltseva  \n3.7  \nFig. 3. Zoom-ins on the template spectra at wavelength regions centered on the 3.3 µ m AIB (Peeters et al. 2023, top), the 6.2 µ m AIB (second from top), the 7.7 µ mAIB(second from bottom), and the 11.2 µ mAIB(bottom). Each spectrum (on an F ν scale) is normalised by the peak surface brightness of the indicated AIB on the y-axes in each panel. The vertical tick marks indicate the positions of identified (blue) or tentative (black) AIBs and components (see Table A.1 and main text). A post-pipeline correction for residual artifacts was performed for Ch2-long (10.02-11.70 µ m), Ch3-medium (13.34-15.57 µ m), and Ch3-long (15.41-17.98 µ m). For further details, see Appendix B. Red dashed vertical ticks indicate the wavelengths where we switch from using data from one MRS sub-band to another. Continued in Fig. 4. Article number, page 6 of 25  \n<!-- image -->  \nFig. 4. Continued from Fig. 3. From left to right, top to bottom: Zoom-ins on the template spectra normalised by the peak flux of the 5.25, 5.75 and 5.878, 6.2, 8.6, 11.2 and 12.0, 12.7, 13.5, and 14.2 µ mAIBs (indicated in the y-axis label of each panel). The panels show wavelength ranges that are also shown in Fig. 3, except for the panels that are centered on the 13.5 and 14.2 µ mAIBs (small panels in the lower right). These figures illustrate the overall similarity and subtle di ff erences in AIB profiles from region to region.  \n<!-- image -->  \net al. 2018; Buragohain et al. 2020). For hydrogenated PAHs, the main activity is at slightly longer wavelengths, ≃ 3 . 5 µ mwithin a somewhat narrower range ( ≃ 0 . 05 µ m). As the extra hydrogens in superhydrogenated PAHs are relatively weakly bound (1.41.8 eV; Bauschlicher & Ricca 2014), astronomical models imply that superhydrogenated PAHs quickly lose all these sp 3 hydrogens in strongly irradiated PDRs (e.g. when G 0 / n (H) > 0 . 03; Andrews et al. 2016).  \nFor further analysis of the AIB emission in this region, including the many weaker features listed in Table A.1, we refer to Peeters et al. (2023).", '3.2. The 5-6 µ m (1600-2000 cm -1 ) range': 'In the 5-6 µ m region, previous observations have revealed two moderately weak AIB features at approximately 5.25 and 5.75 µ m (Table A.1; Allamandola et al. 1989a; Boersma et al. 2009b). The 5.25 µ m band (Fig. 4) consists of a broad  \nFig. 5. Variation from region to region in the ratios of AIB integrated fluxes (top) and the widths of AIBs (bottom; data from Table 1). All values are normalised to those from the atomic PDR. The ordering of the regions on the x-axis is arbitrary, but was chosen so that the data roughly show an increasing trend in FWHM from left to right (Sect. 4.3).  \n<!-- image -->  \nblue shoulder centered at ∼ 5 . 18 µ m and extending to about 5.205 µ m, followed by a sharp blue rise to a peak at 5.236 µ m and a strong red wing extending to about 5.38 µ m. A detailed inspection of the profiles reveals a very weak feature at ∼ 5 . 30 µ m superposed on the red wing. Comparing the five template spectra, we conclude that the 5.25 µ m feature broadens, in particular on the red side, with the narrowest feature seen in the atomic PDR, and then increasing in width in the H ii region, DF 1, DF 2, and DF 3 (Table 1 and Fig. 5). Besides the broadening, the observed profiles are very similar for the five templates, implying that the main feature consists of a single band.  \nInspection of the template spectra (Fig. 4) reveals that the 5.75 µ m band is a blend of three bands at 5.642, 5.699, and 5.755 µ m (e.g. comparing the atomic PDR and DF 3 spectra in Fig. 4). The MIRI MRS spectra clearly exhibit a new, symmetric feature at 5.878 µ m. We also report a tentative detection of two very weak features at 5.435 and 5.535 µ m.  \nThe spectra of PAHs show weak combination bands in this wavelength range generated by modes of the same type, for example, out-of-plane (OOP) modes (Boersma et al. 2009b,a).  \nTable 1. FWHMof AIBs and AIB complexes measured from five template spectra extracted from NIRSpec-IFU (first row) and MIRI-MRS (rows 2-6) observations of the Orion Bar.  \nNotes. The FWHM of the 8.6 µ m AIB cannot be determined, for the peak intensity of the plateau between the 7.7 µ m and 8.6 µ m AIBs is greater than half the peak intensity of the 8.6 µ m AIB (see the third panel of Fig. 3).  \nColumns: (1) central wavelength of AIB; (2) central wavenumber of AIB; (3)-(7) FWHM (cm -1 ) of AIB complex from indicated template spectrum; (8) peak-to-peak variation in FWHM (cm -1 ) observed in the templates (see Fig. 5).  \nCombination bands involving in-plane modes occur at shorter wavelengths (3 . 8 -4 . 4 µ m) and are typically an order of magnitude weaker (Mackie et al. 2015, 2016). Combination bands involving the OOP bending modes typically result in a spectrum with two relatively simple AIBs near 5.25 and 5.75 µ m. For small PAHs, the ratio of the intrinsic strength of these bands to the OOP modes increases linearly with PAH size (Lemmens et al. 2019). This ratio increases further for the larger PAHs studied in Lemmens et al. (2021). However, whether the correlation continues linearly is yet to be confirmed.', '3.3. The 6.2 µ m (1610 cm -1 ) AIB': 'The interstellar 6.2 µ m band is one of the main AIBs. The profile peaks at 6.212 µ m (1610 cm -1 ) and has a steep blue rise, a pronounced red wing, and a blue broad shoulder centered at ∼ 6 . 07 µ m. Comparing the five template spectra, we conclude that the feature broadens toward the blue side at the same time that the red wing becomes more pronounced (by about 8.2 cm -1 on a total width of ≃ 33 . 6 cm -1 ; Table 1). Pending confirmation, the peak position possibly varies (its value ranges from 6.2115 to 6.2161 µ m). There is a distinct weaker feature at 6.024 µ m (1660 cm -1 ) superposed on the blue shoulder. This symmetric feature has a constant width and varies in intensity independently of the main feature (see Fig. 4). This suggests that the 6.024 µ m band is an independent component. We note that the observed strength variations of the 6.024 µ m band do not a ff ect the conclusion on the broadening of the blue side of the 6.2 µ m band. There is a very weak feature perched on the red wing at 6.395 µ m (1564 cm -1 ) in the template of the atomic region. It may be obscured by the stronger red wing in the other template spectra. A very subtle change in slope of the red wing may also be present near 6.5 µ min some templates (e.g. the atomic PDR).  \nPure aromatic CC stretching modes fall between 6.1 and 6.5 µ m. In the 6-9 µ m wavelength range, the number of bands and their precise positions will depend on charge, molecular structure, size, and heterogeneity. In particular, their intrinsic strengths are very sensitive to the charge state of the species, increasing by about a factor of 10 for cations (Allamandola et al. 1999; Peeters et al. 2002a; Bauschlicher et al. 2008). Given the  \nobserved strength of the 6.2 µ m band relative to the CH stretching and OOP bending modes that dominate the neutral spectra, this AIB is attributed to PAH cations (Allamandola et al. 1999). In the past, the peak position was somewhat of an enigma. In early comparisons with harmonic calculations, this band arose at too red a wavelength in PAH cations. This problem was compounded by the adoption of a redshift of 15 cm -1 to account for anharmonic e ff ects during the emission cascade (for a summary, see Bauschlicher et al. 2009). However, model studies have revealed that anharmonicity introduces a red wing on the profile but does not lead to an appreciable redshift of the peak (Mackie et al. 2022). Recent experimental and quantum chemical studies of neutral, symmetric PAHs have shown that the mismatch between the experimental and interstellar 6.2 µ m band positions is less severe than thought (Lemmens et al. 2021). Furthermore, quantum chemical studies on PAH cations have employed the cc-pVTZ basis set that better accounts for treatment of polarization in PAHs. With this basis set used in density functional theory calculations, the calculated peak position of the aromatic CC stretching mode in cations is in much better agreement with the observations (Ricca et al. 2021), but this still needs to be confirmed by experimental studies on PAH cations. The discrepancy noted in earlier studies between the peak position of the 6.2 µ mAIB and the aromatic CC stretch in PAHs has prompted a number of suggestions. Specifically, incorporation of heteroatoms such as N into the ring backbone or coordination of atoms such as Si, the presence of aliphatic structures, protonated PAHs, and / or (pentagonal) defects will induce blue shifts in the peak position of this mode (Hudgins et al. 2005; Pino et al. 2008; Joalland et al. 2009; Carpentier et al. 2012; Galué 2014; Tsuge et al. 2018; Wenzel et al. 2022; Rap et al. 2022). Further studies are warranted to assess whether these suggestions are still relevant.  \nThe observed 6.024 µ m band is at too short a wavelength to be an aromatic CC stretching vibration. Rather, this position is characteristic of the C = Ostretch in conjugated carbonyl groups; that is, as quinones or attached to aromatic rings (Allamandola et al. 1989b; Sarre 2019). This band has not been the focus in many quantum chemical studies. We also note that the very weak feature at 6.395 µ m is likely another aromatic CC stretching mode.', '3.4. The 7.7 µ m (1300 cm -1 ) AIB complex': 'It has been well established that the 7.7 µ m AIB is a blend of several features (Bregman 1989; Cohen et al. 1989; Peeters et al. 2002a). The JWST spectra reveal that the main component at 7.626 µ m is accompanied by moderately strong bands at 7.8 and 7.85 µ m. The 7.8 µ m component appears narrower in DF 2 and DF 3, peaking near 7.743 µ m, although this may arise due to di ff erences in the red wing of the 7.626 µ m component or this may reflect the lack of a di ff erent component at 7.775 µ m present in the atomic PDR, the H ii region and DF 1. In any case, given the observed variations between the templates, these bands are independent components. The 7.7 µ mAIBcomplex as a whole broadens significantly from the atomic to the molecular region (by about 18.8 cm -1 , Table 1). In addition to these moderately strong components, there are also weak features at 7.24 and 7.43 µ mand between the 7.7 and 8.6 complexes at 8.223 and 8.330 µ m. Very weak features are also present at shorter wavelengths (6.638, 6.711, 6.850, 6.943, 7.05, and 7.10 µ m).  \nBands in this wavelength range are due to modes with a mixed character of CC stretching and CH in-plane bending vibrations. As mentioned in Sect. 3.3, the strength of these modes  \nis very dependent on the charge state of the species and the interstellar 7.7 µ mAIB is generally attributed to PAH cations. The spectra of very symmetric PAHs become more complex with increasing size and the main band(s) in the 7 . 5 -8 . 0 µ mrange shift systematically with size toward longer wavelength from about 7.6 to about 7.8 µ m or even larger (Bauschlicher et al. 2008, 2009; Ricca et al. 2012). These quantum chemical calculations point toward compact PAHs in the size range 24 to 100 C atoms as the carriers and probably more toward the smaller size for the main 7.626 µ m band and slightly larger for the two moderate components. Detailed spectral decompositions of ISO-SWS and Spitzer-IRS observations agree with these conclusions (Joblin et al. 2008; Shannon & Boersma 2019). The very weak features at 7.24 and 7.43 µ m are likely also CC stretching modes. The CH deformation modes of aliphatic groups also occur around 6.8 and 7.2 µ m, but these modes are weaker compared to the CH stretching modes of aliphatic groups around 3.4 µ m(Wexler 1967; Yang et al. 2016; Dartois et al. 2005) and, given the weakness of the 3.4 µ mAIB, we deem that identification unlikely for the weak features at 7.24 and 7.43 µ m detected in all templates. The very weak features at 6.850 and 6.943 µ mare only present in DF 2 and DF 3. As both these templates also show the strongest 3.4 µ m emission, these bands may arise from CH deformation modes of aliphatic groups (Wexler 1967; Arnoult et al. 2000).', '3.5. The 8.6 µ m (1160 cm -1 ) AIB': "This AIB peaks at 8.60 µ m (1163 cm -1 ). The apparent shift toward shorter wavelengths in the DF 3 spectrum as well as the apparent broadening of the band are likely caused by the change in the underlying 'continuum' due to the 7.7 µ m AIB and / or plateau emission and / or very small grain emission. The change in slope in the blue wing at 8.46 and 8.54 µ m suggests the presence of more than one component in this AIB. However, these components seem to be very weak compared to the main band. There is a similar change in slope at 8.74 µ m and potentially at 8.89 µ m in all template spectra and this likely has a similar origin. Since these features are very weak, we label them as tentative.  \nThe 8.6 µ m AIB is due to CH in-plane bending modes in PAHs, but this mode has a large CC stretching admixture. The intensity of this band increases significantly and it shifts to longer wavelength, producing the very prominent band that appears near 8.5 µ min the spectra of large ( N C ∼ 100, N C being the number of C atoms in a PAH molecule) compact PAHs (Bauschlicher et al. 2008). For even larger, compact PAHs, this band starts to dominate the spectra in the 7 -9 µ m range and these species are excluded as carriers of the typical AIB emission (Ricca et al. 2012). In large polyaromatic and aliphatic systems, the geometrical distortions of the C-C backbone and defects, partly related to the hydrogen content, shift the position of this band (Carpentier et al. 2012; Dartois et al. 2020). The weaker features on the blue side of the main band may be due to somewhat smaller and / or less symmetric PAHs while the longer wavelength feature may be due to a minor amount of somewhat larger symmetric, compact PAHs.", '3.6. The 10-20 µ m (500-1000 cm -1 ) range': "This wavelength range is dominated by the strong AIB at 11.2 µ m, a moderately strong AIB at 12.7 µ m and a plethora of weaker AIBs at 10.95, 11.005, 12.0, 13.5, 13.95, 14.21, and 16.43 µ m. The 11.2 µ m AIB clearly displays two components  \nat 11.207 and 11.25 µ m, along with a tentative component at 11.275 µ m. The AIB peaks at the first component (11.207 µ m) in the atomic PDR, the H ii region, and DF 1, while it peaks at the second component (11.25 µ m) in DFs 2 and 3. These two components may shift to longer wavelengths in DFs 2 and 3, however, such shifts are still yet to be confirmed. The relative strengths of these two components vary across the five templates, indicating they are independent components. The combined 11.2 µ m profile is asymmetric with a steep blue rise and a red wing. The AIB broadens significantly (by ∼ 5 . 9 cm -1 on a total width of ∼ 12 . 2 cm -1 ; Table 1) through the atomic PDR, DF 1, the H ii region, DF 2, and DF 3 in increasing order. This broadening is driven by changes in the red wing though similar but very small changes in the steepness of the blue side are present. An additional weaker component may be present on the red wing at 11.275 µ m. In addition, similar to the 5.25 µ m AIB, the 11.2 µ m AIB displays a broad blue, slow-rising, shoulder from ∼ 10 . 4 µ m up to the start of the steep blue wing. A wellknown weaker AIB is present at 11.005 µ m and is superposed on this blue shoulder.  \nThe 12.0 µ m band peaks at 11.955 µ m and may have a second component at 12.125 µ m. The template spectra furthermore display elevated emission between the red wing of the 11.2 µ m band and the 12.2 µ m band (see e.g. DF 3) suggestive of more complex AIB emission than expected based on the presence of these two bands. However, due to an artefact at 12.2 µ m(see Appendix B), confirmation of the second component, the 12.0 µ m profile, and this elevated emission between the 11.2 and 12.0 µ m bands requires further improvements to the calibration.  \nThe 12.7 µ mband is very complex displaying a terraced blue wing and a steep red decline. It peaks at 12.779 µ mexcept in the atomic PDR where it peaks at a second component at 12.729 µ m. The strengths of these components vary independently from each other. Three additional terraces are located near 12.38, 12.52, and 12.625 µ m and a red shoulder near 12.98 µ m suggests the presence of an additional component. Given the complexity of the 12.7 µ m band, spatial-spectral JWST maps are required to understand its spectral decomposition into its numerous components. The entire 12.7 µ mcomplex significantly broadens largely on the blue side but also on the red side. We report a broadening (by ∼ 14 . 2 cm -1 on a total width of ∼ 21 . 9 cm -1 ; Table 1 and Fig. 5) through the atomic PDR, DF 1, the H ii region, DF 2, and DF 3 in increasing order.  \nThe 13.5 µ m band peaks at 13.55 µ m and may be accompanied by two additional components at 13.50 and 13.62 µ m. The 13.5 µ m band seems to broaden as well. We note that several artefacts exist just longwards of this band (Appendix B), hampering its analysis. Hence, future improvements to the calibration and additional observations on a wider range of sources will have to confirm this broadening. These artefacts also limit the detection and analysis of bands in the 14-15 µ m range. We detected a band at 14.21 µ mand potentially at 13.95 µ m, although the latter is just to the red of the artefact at 13.92 µ m.  \nBands in the 11-14 µ m range are attributed to CH OOP bending modes. The peak position and pattern of these bands is very characteristic for the molecular edge structure of the PAH; that is, the number of adjacent H's 5 . The bands making up the 11.2 µ m AIB can be ascribed to neutral species with solo H's (Hony et al. 2001; Bauschlicher et al. 2008). The cationic solo  \nH OOP band falls at slightly shorter wavelength than the corresponding solo H OOP band of neutral PAHs and the 11.0 µ m AIB has been attributed to cations (Hudgins & Allamandola 1999; Hony et al. 2001; Rosenberg et al. 2011). The 12.7 AIB complex is due to either duo H's in neutral PAHs or trio H's in cations. For species with both solo and duo H's, coupling of the duo with the solo CH OOP modes splits the former into two bands. The sub-components in the 12.7 µ mAIB may reflect this coupling and / or it may be caused by contributions of more than one species with duo's.  \nThe weak 12.0 µ m AIB can be attributed to OOP modes of duo H's, while the 13.5 µ m AIBs likely have an origin in OOPmodes of quartet H's in pendant aromatic rings (Hony et al. 2001; Bauschlicher et al. 2008). The weak bands near 14.2 µ m could be due to OOP modes of quintet H's. Alternatively, for larger PAHs, CCC skeletal modes are present in this wavelength range (Ricca et al. 2012).  \nWealso detected a band at 16.43 µ m. Other weaker bands are present in this region but due to calibration issues (Appendix B), we refrained from characterizing them.", '4.1. Comparison to previous observations': "Overall, in terms of spectral inventory, the observed AIB emission in the 3 µ m range is consistent with prior high-quality ground-based observations of the Orion Bar (e.g. Sloan et al. 1997). Likewise, the main characteristics of the AIB emission are also detected in prior observations of the Orion Bar carried out with ISO-SWS (Verstraete et al. 2001; Peeters et al. 2002a; van Diedenhoven et al. 2004) and Spitzer-IRS in short-low mode (Knight et al. 2022a). Furthermore, in retrospect, many (weaker) bands and sub-components of the AIB emission seen by JWST may also be recognised in the ISO-SWS observation of the Orion Bar, but they were too weak and too close to the S / N limit to be reported in previous works. However, as these JWST data have an unparalleled combination of extremely high S / N, spectral resolution and, most importantly, superb spatial resolution, these spectral imaging data reveal already known bands and subcomponents in unprecedented detail allowing for a much improved characterization of the AIB emission. In addition, these spectral imaging data reveal previously unreported components (blends) and sub-components of the AIB emission (indicated in Table A.1 and discussed in Sect. 3).  \nThe AIBs at 5.75, 7.7, 8.6, 11.2, and 12.7 µ m have complex sub-components. Boersma et al. (2009b) noted that the 5.75 AIB has an unusual profile, resembling a blended double-peaked feature. The JWST template spectra indicate that the band is composed of three components with variable strengths. New components are also seen in the 12.7 µ m band. Shannon et al. (2016) reported that the 12.7 µ m band shifts to longer wavelengths at larger distances from the illuminating star in reflection nebulae. This is consistent with the behavior of this band in the Orion Bar reported here where it reflects the relative intensities of the two components at 12.729 and 12.779 µ m. These authors also reported a change in the blue wing. The JWST data now characterizes the components (i.e. terraces) in the blue wing and their relative intensities.  \nWhile the sub-components of the 8.6 µ m AIB have not, to our knowledge, been reported in the literature, several studies detail sub-components in the 7.7, 11.2, and 12.7 µ m AIBs. The 7.7 µ mAIB complex is composed of two main sub-components at ∼ 7.626 and ∼ 7.8 µ m (Cohen et al. 1989; Bregman 1989;  \nPeeters et al. 2002a). The 7.7 µ m AIB complex is distinguished into four classes (A, B, C, and D) primarily based on its peak position (Peeters et al. 2002a; Sloan et al. 2014; Matsuura et al. 2014). Spectral-spatial imaging has revealed that the (class A) 7.7 µ m profile varies within extended ISM-type sources and depends on the local physical conditions: the 7.8 µ m component gains in prominence relative to the 7.626 µ m component and is accompanied by increased emission 'between' the 7.7 µ m and 8.6 µ m AIBs in regions with less harsh radiation fields (Bregman & Temi 2005; Berné et al. 2007; Pilleri et al. 2012; Boersma et al. 2014; Peeters et al. 2017; Stock & Peeters 2017; Foschino et al. 2019; Knight et al. 2022b). Our findings using the JWST Orion Bar templates (Fig. 3) are consistent with these past results. Pilleri et al. (2012) attributed this to an increased contribution of evaporating very small grains (eVSGs). In addition, the JWST data reveal that the 7.8 µ m component is composed of three components whose relative contribution varies.  \nLikewise, the 11.2 µ m AIB has been classified into class A11 . 2, B11 . 2, and A(B)11 . 2. Class A11 . 2 peaks in the 11.2011.24 µ mrange and displays a less pronounced red wing relative to the peak intensity (corresponding to a FWHM of ∼ 0.17 µ m), class B11 . 2 peaks at ∼ 11.25 µ m and shows a more pronounced red wing (FWHM of ∼ 0.20 µ m), and class A(B)11 . 2 is a mix with a peak position as that of class A11 . 2 and prominence of its red wing as that of class B11 . 2 (resulting in a FWHM of ∼ 0.21 µ m; van Diedenhoven et al. 2004). As for the 7.7 AIB complex, ISM-type sources display a class A11 . 2 profile. Recent spectralimaging data however indicated that the 11.2 µ m profile shifts to slightly longer wavelengths accompanied with a stronger red wing relative to the peak intensity in two (out of 17) positions of the Orion Veil (Boersma et al. 2012) and in two reflection nebulae (Boersma et al. 2013; Shannon et al. 2016). These authors classified these profile variations as a shift from class A11 . 2 to class A(B)11 . 2, which, in the case of the two reflection nebulae, occurred when moving away from the illuminating star. Boersma et al. (2014) linked the change in the 11.2 µ mprofile to a change in the 7.7 µ mAIB complex (probed by the 11.2 / 11.3 and 7.6 / 7.8 intensity ratios, respectively). A change in peak position along with a broadening of the profile is consistent with the JWST templates of the Orion Bar. As discussed in Sect. 3, the change in the peak position of the 11.2 µ mAIBreflects the relative importance of two components at 11.207 and 11.25 µ m that are now clearly discerned in the JWST data. Furthermore, thanks to the increased spectral and spatial resolution, we conclude that DF 3 belongs to class B11 . 2 (Fig. 6). Hence, the Orion Bar exhibits class A11 . 2 profiles near the surface of the PDR which evolved from class B11 . 2 profiles deeper in the molecular zone.  \nThe ISO-SWS observations of the Orion Bar (taken in a 14 '' × 20 '' aperture) resemble the atomic PDR template, even when centered on DF 3 6 . This resemblance is due to the fact that the AIB emission is significantly stronger in the atomic PDR compared to the molecular PDR (Habart et al. 2023; Peeters et al. 2023) and it dominates the emission within the large ISO / SWS aperture. Hence, the JWST spectrum of the atomic PDR in the Orion Bar (Fig. 2) serves as the updated, high-resolution, more detailed template spectrum for class A AIB emission. The DF 2 and DF 3 templates, which probe regions deep in the molecular PDR, no longer exhibit class A11 . 2 profiles while the 3.3, 6.2, 7.7, and 8.6 µ m AIBs still clearly belong to class A. A similar situation, where individual targets are found to belong to two classes, has been reported for two targets: the planetary nebula Hb 5 and  \nFig. 6. Comparison of the 11.2 µ m profile in the five template spectra with a class A 11.2 µ mprofile represented by the ISO-SWS spectrum of the Orion Bar H2S1 (van Diedenhoven et al. 2004, top panel) and a class B 11.2 µ m profile represented by the ISO-SWS spectrum of HD 44179 (van Diedenhoven et al. 2004, bottom panel).  \n<!-- image -->  \nthe Circinus galaxy (van Diedenhoven et al. 2004). These authors furthermore found that the other two galaxies in their sample display class A profiles for the 3.3, 6.2, 7.7, and 8.6 µ mAIBs, while displaying a class A(B)11 . 2 AIB profile. This suggests that, out of the main AIBs, the 11.2 µ m AIB is the cleanest indicator of the shift from class B to class A.", '4.2. AIB profiles': "Broadly speaking, the prominent AIBs can be separated into three groups: 1) bands with a steep blue rise and a pronounced red wing. The 5.25, 6.2, and 11.2 µ m AIBs are clear examples. The profiles also often show a shoulder on the blue side, which is considerably weaker than the red wing; 2) bands that are clear blends of multiple components. This group includes the 3.4, 5.75, 7.7, and 12.7 µ m AIBs. They typically comprise three or more sub-components.; and 3) bands that seem to be symmetric, often resembling Gaussian profiles. This group includes the 3.3, 5.878, and 8.6 µ m AIBs, as well as the 6.024 and 11.005 µ mAIBs. These divisions are not entirely strict. The  \ngroup 1 AIBs typically have very weak features perched on the red and / or blue side of the profile. Likewise, it is conceivable that each of the blended components in the group 2 AIBs might have an intrinsic profile with relatively sharp blue rise and a more gradual red wing that is obfuscated by blending. For example, while the 3.403 µ mAIB is often blended with a feature at 3.424 µ m, this component is very weak and the 3.403 µ m profile resembles that of group 1. We note that while the 11.2 µ m AIB is also a blend of three components, the character of the profile is dominated by the presence of a steep blue rise and a pronounced red wing rather than the presence of the sub-components. Therefore, we list the 11.2 µ mAIB in group 1.  \nProfiles with a steep blue rise and a pronounced red wing are characteristic for the e ff ects of anharmonicity (Pech et al. 2002; Mackie et al. 2022). Detailed models have been developed that follow the emission cascade for a highly excited, single PAH and that include the e ff ects of anharmonic interactions based on quantum chemical calculations (Mackie et al. 2022). These models do not contain free parameters besides the size of the emitting species (i.e. the average excitation level after absorption of a FUV photon) and the resulting profiles agree qualitatively well with the observations of group 1 AIBs (Mackie et al. 2022). The results show that the wavelength extent of the red wing depends on the details of the anharmonic coupling coe ffi cients with other modes. The strength of the wing relative to the peak emission is sensitive to the excitation level of the emitting species after absorption of the UV photon (i.e. the initial average energy per mode) and the cascade process (e.g. how fast the energy is 'leaking' away through radiative cooling). The steepness of the blue rise is controlled by rotational broadening. Analysis of the profile of the 11.2 µ mAIBobserved by ISO / SWSsuggests emission by a modestly sized PAH ( N C ∼ 30; Mackie et al. 2022) but that conclusion has to be reassessed given the presence of more than one component in this AIB in the JWST spectra of the Orion Bar.  \nNot all bands will show equally prominent anharmonic profiles. In particular, the far-infrared modes in small PAHs are very harmonic in nature (Lemmens et al. 2020) and their profiles would not develop red wings. Likewise, the aromatic CH stretching modes are very susceptible to resonances with combination bands (Maltseva et al. 2015, 2016; Mackie et al. 2015, 2016) and this interaction dominates their profiles (Mackie et al. 2022).", '4.3. AIB variability in Orion': "Spatial-spectral maps carry much promise to untangle the complexity of the AIBs and possibly link observed variations to the presence of specific carriers. The first forays into this field were based on Spitzer spectral maps. Analysis of the spatial behaviour of individual AIBs and AIB components revealed their interdependence as well as new components (e.g. Boersma et al. 2009b; Peeters et al. 2012; Boersma et al. 2013, 2014; Shannon et al. 2016; Peeters et al. 2017). Analyses based on blind signal separation methods have uncovered several distinct components and spectral details (Berné et al. 2007; Joblin et al. 2008; Pilleri et al. 2012; Foschino et al. 2019), but the increased spectral and spatial resolution, as well as the higher sensitivity of JWST , can be expected to take this to a new level. Indeed, some of the spectral details uncovered by previous spatial-spectral studies are now directly detected in the presented JWST data of the Orion Bar. Applications of blind signal separation techniques, as well as spectral fitting using PAHFIT (Smith et al. 2007) and the Python PAH Database (Matthew J. Shannon & Christiaan Boersma 2018) on the full spectral map of the Orion Bar may be a promising ground  \nfor additional detections and potential identifications. Here, we address spatial-spectral variations on inspection of the five template spectra. More detailed analyses are deferred to future studies.  \nWhile the spectra are rich in components (Table A.1), there is little diversity between the template spectra. All templates show evidence for each sub-components, except possibly for a few very weak bands whose presence may be easily lost in the profiles of nearby strong bands. The most obvious variations are the increased prominence of the sub-components in the 7.7 µ mAIB at 7.743 and 7.85 µ mand in the 11.2 µ mAIB at 11.25 µ min the DF 2 and DF 3 spectra, and the variation in the relative strength of the sub-components of the weak 5.75 µ mband. Similarly, the width of many AIBs (3.3, 5.25, 6.2, 7.7, 11.2, and 12.7) broadens significantly (Table 1). This broadening is also systematic, with the FWHM being smallest in the atomic PDR, increasing subsequently in DF 1, followed by the H ii region, then DF 2, and finally DF 3 (Table 1 and Fig. 5). The only exception to this systematic trend is the H ii region having a FWHM smaller than DF 1 (but larger than the atomic PDR) for the 5.25 and 7.7 µ m AIBs. As pointed out in Sect. 2.2, we note that the AIB emission in the H ii region originates from the background PDR.  \nThere is also some region-to-region variation in the relative strengths of the main AIBs (Fig. 5). The largest variations are seen in the 3.4 µ m AIB to 3.3 µ m AIB ratio ( ∼ 100% greater in DF 3 than in the atomic PDR) and in the 3.3 µ mAIB to 11.2 µ m AIB ratio ( ∼ 40% greater in DF 3 than in the atomic PDR). Both of these ratios are largest in DF 2 and DF 3, and decrease in the atomic region. Variations in the strength of the CC modes and in-plane CH bending modes (6.2, 7.7, and 8.6 µ m) relative to the 11.2 µ m OOP modes are more modest (at the 10-20% level). Larger variations occur in the relative strength of the moderate bands, namely, the 5.25, 5.878, 6.024, and 11.955 µ m bands which are much more pronounced in DF 2 and DF 3 than in the atomic zone. As noted in Sect. 4.1, the ISO-SWS data of the Orion Bar resembles the atomic PDR template. Hence, the range in spectral variability within class A AIBs is well represented by the five templates for all AIBs except the 11.2 µ m AIB. For the 11.2 µ mAIB, the presented data not only showcase the class A AIB variability but also the shifts from class A to class A(B) and then to class B. It is expected that future JWST observations probing a large range of physical conditions and environments further extend the spectral variability in the AIB emission.  \nThe anharmonic profile of bands due to smaller PAHs will tend to have a less steep blue rise due to the increase in the rotational broadening as well as a slight increase in the width and a more pronounced red wing due the higher internal excitation for the same photon energy (Mackie et al. 2022; Tielens 2021). Hence, the presence of somewhat smaller PAHs may be at the origin (of some) of the overall profile variations in the 3.4, 5.25, 6.2, and 11.2 µ m AIBs in the JWST template spectra. We note that the spectrum of the DF 1 template resembles much more the atomic region template than the DF 2 and DF 3 dissociation front templates (Figs. 3 and 4, as well as Table 1). This is likely due to the terraced-field-like structure of the molecular PDR resulting in a strongly enhanced line-of-sight visual extinction through the foreground atomic PDR toward DF 1 compared to DF 2 and DF 3 (Habart et al. 2023; Peeters et al. 2023). Hence, a large contribution from the atomic region in the foreground is contributing to the emission toward DF 1. The overall similarity of the spectra suggests that the PAH family is very robust but has a small amount of additional species in the DF 2 and DF 3 zones that is not present in the surface layers of the PDR. In the Orion Bar, the PDR material is advected from the molecular zone to  \nthe ionization front at about 1 km s -1 (Pabst et al. 2019) over ≃ 20 , 000 yr. In that period, a PAH will have absorbed some 10 8 UV photons and, yet, apparently the e ff ect on the composition of the interstellar PAH family is only minor as it only results in a change in the prominence of the sub-components in the 7.7 and 11.2 µ mAIBs. This likely reflects that for moderate-to-large PAHs ( N C ≳ 30), photofragmentation is a minor channel compared to IR emission and, moreover, when fragmentation occurs, the H-loss channel dominates over C loss (Allain et al. 1996a,b; Zhen et al. 2014b,a; Wenzel et al. 2020) and then rapidly followed by rehydrogenation with abundant atomic H (Montillaud et al. 2013; Andrews et al. 2016). We note that the UV field increases by about two orders of magnitude between the H2 dissociation front and the PDR surface but the atomic H abundance increases by a similar factor as H2 is increasingly photolysed near the surface. Hence, the ratio of the local FUV field to the atomic hydrogen density, G 0 / n (H), which controls the photoprocessing (Andrews et al. 2016), does not vary much among the five template regions. We thus suggest that the additional species in the deeper layers of the Orion Bar causing the increased prominence of sub-components in the 5.75 (at 5.755 µ m), 7.7, and 11.2 µ m AIBs, are aromatic species and / or functional groups and / or pendant rings that are more susceptible to photolysis.  \nPhotolytic processing of the PAH family with position in the Orion Bar may also leave its imprint on the PAH size distribution and this will a ff ect the relative strength of the AIBs. The 3.3 / 11.2 AIB ratio has long been used as an indicator of the size of the emitting species (Allamandola et al. 1989b; Pech et al. 2002; Ricca et al. 2012; Mori et al. 2012; Croiset et al. 2016; Maragkoudakis et al. 2020; Knight et al. 2021, 2022a) as this ratio is controlled by the 'excitation temperature' of the emitting species and, hence, for a fixed FUV photon energy, by the size. For the five JWST template regions, the 3.3 / 11.2 AIB ratio is observed to vary by about 40%, being largest in DF 3 and decreasing toward DF 2, followed by DF 1, the atomic PDR and the H ii region. The variation in this ratio is slightly less than what is observed in the reflection nebulae NGC 7023 and in the larger Orion region (e.g. the Orion Bar and the Veil region beyond the Orion Bar) and corresponds to an increase in the typical size of the emitting species by about 40% toward the surface (Croiset et al. 2016; Knight et al. 2021, 2022a; Murga et al. 2022). Hence, we link the decreased prominence of the sub-component in the 5.75 µ m (at 5.755 µ m), 7.7 µ m (at 7.743 and 7.775 µ m), and 11.2 µ m AIBs (at 11.275 µ m) as well as the variation in the 3.3 / 11.2 AIB ratio to the e ff ects of photolysis as material is advected from the deeper layers of the PDR to the surface.  \nVariations in the 6.2 / 11.2 ratio (Fig. 5) are generally attributed to variations in the ionised fraction of PAH (Peeters et al. 2002a; Galliano et al. 2008; Stock et al. 2016; Boersma et al. 2018). This ratio is only 12% stronger in DF 3 than the surface of the PDR. The limited variations in the 6.2 / 11.2 AIB ratio is at odds with those measured by Spitzer and ISO . Specifically, this ratio is observed to increase by about 50% across the Orion Bar when approaching the Trapezium cluster (Knight et al. 2022a). Moreover, Galliano et al. (2008) measured an increase in this ratio by almost a factor of 2 over (a much wider swath of) the Orion Bar. The PAH ionization balance is controlled by the ionization parameter, G 0 T 1 / 2 / ne with G 0, T , and ne the intensity of the FUV field, the gas temperature, and the electron density. The high spatial resolution of JWST allows for a clear separation of the emission at the H2 dissociation fronts and the PDR surface. The limited variation in the 6.2 / 11.2 AIB ratio is somewhat surprising because the PAH ionizing photon flux di ff ers by about a factor of 40 between the dissociation fronts (located at  \nAV = 2 mag) and the PDR surface; the gas temperature will also increase (slightly) toward the surface, while the electron abundance remains constant over this region (Tielens & Hollenbach 1985). It is also possible that PAH cations contribute an appreciable amount to the 11.2 µ m band (Shannon et al. 2016; Boersma et al. 2018). Further modelling will be important to fully understand the complexity of the Orion Bar (A. Sidhu et al., in prep.).  \nThe 3.4 / 3.3 AIB ratio is observed to decrease by about 100% from DF 3 to the atomic PDR (Fig. 5). Such variations have also been seen in other nebulae and attributed to photofragmentation processes (Joblin et al. 1996). The 3.4 µ m AIB is attributed to aliphatic CH modes in the form of a minor amount of H bonded to sp 3 C atoms either in the form of methyl groups or as superhydrogenated PAHs (Schutte et al. 1993; Bernstein et al. 1996; Joblin et al. 1996). The abundance of superhydrogenated PAHs is expected to be very small throughout the Orion Bar as such extra H's are readily lost in the strong FUV radiation field (Andrews et al. 2016). Methyl groups are also more easily photolysed than aromatic H's (energy barriers are 3.69, 4.00, and 4.47 eV for CH2-H, -CH3 and aromatic H loss, respectively; Tielens 2021). Recent experiments report that, for cations, this methyl group photolysis can lead to quite stable tropylium formation (loss of H followed by isomerization to a seven-membered ring; Jochims et al. 1999; Zhen et al. 2016; Wenzel et al. 2022). However, further investigation is required to firmly establish the importance of this fragmentation route for conditions present in the interstellar medium. If borne out, the reaction of the tropylium cation with atomic H has a calculated barrier of 3.2 kcal mol -1 (1600 K; Bullins et al. 2009) and, hence, under warm, dense Hrich conditions, the methyl group could be reformed. Hence, in a 'suitable' PDR, the species may cycle back and forth between a methyl functional group and the tropylium structure until eventually -CH3 loss occurs. In any case, it can be expected that the stronger UV field nearer to the surface will reduce the number of CH methyl groups compared to the number of CH aromatic bonds. Further experimental and quantum chemical studies will have to address the competition between the various channels involved in the chemistry of methylated PAHs in PDRs.  \nThe increased importance of PAH photolysis near the surface of the Orion Bar PDR is in line with the GrandPAH hypothesis that only the most resilient species can sustain the harsh conditions of strong FUV radiation fields. Thus, a limited number of compact, large PAHs will dominate the interstellar PAH family in these conditions (Andrews et al. 2015; Tielens 2013). If the conditions are right, large PAHs may even be stripped of all their H's and isomerize to the fullerene, C60 (Boersma et al. 2012; Berné et al. 2015; Zhen et al. 2014a). We also realize that the presence of somewhat smaller PAHs and the increased importance of aliphatic functional groups deep in the PDR may reflect the importance of ion-molecule and / or radical chemistry during the preceding dark cloud core phase modifying and / or forming PAHs in a bottom-up scenario akin to that proposed for the formation of benzonitrile, indene, and cyanonaphthalene (McGuire et al. 2018; Cernicharo et al. 2021; McGuire et al. 2021).", '5. Conclusions': 'The superb sensitivity and spectral resolution of JWST have revealed an ever-better characterization of the AIBs in the Orion Bar in terms of sub-components, multiple components making up a \'single\' band, and band profiles. In addition, the unprecedented spatial resolution of the spectral imaging data showcases the interdependence of the numerous AIB components.  \nWe extracted five template spectra in apertures positioned on the H ii region, the atomic PDR, and the three dissociation fronts DFs 1, 2, and 3. The spectra display a wealth of detail and many weak features have now been firmly identified, their peak positions quantified, and their profiles established. At the same time, the spectra are really very simple. There are a limited number of strong bands with well defined peak positions and red-shaded profiles characteristic for anharmonic interactions. And there is little diversity between the templates. A modest variation is observed in the relative intensities of the main AIBs and of the sub-components of an AIB as well as a systematic broadening of the FWHM of many AIBs (smallest in the atomic PDR and largest in DF 3). Consequently, these templates demonstrate the spectral variations in the class A AIB emission as well as the shift from class B11 . 2 (DF 3) to class A11 . 2 (atomic PDR). The comparison of the template spectra with the ISO-SWS spectrum of the Orion Bar underscores that the spectrum of the atomic region is the "poster child" for the class A spectrum (Fig. 2; Peeters et al. 2002a; van Diedenhoven et al. 2004). This comparison also demonstrates that in a large aperture, PDRs such as Orion are expected to show class A spectra. Conversely, PDRs with more gentle physical conditions (e.g. in the DF 3) are expected to display a slightly modified class A AIB spectrum (except for the 11.2 µ m AIB), showcasing broader AIBs and an increased prominence of minor sub-components with respect to the exemplar class A spectrum. In the case of the 11.2 µ m AIB, more gentle physical conditions broaden the AIB and increase the prominence of minor sub-components seen in class A, resulting in a class B 11.2 µ m AIB profile. Hence, the templates suggest a shift from class B11 . 2 (DF 3) to class A11 . 2 (atomic PDR). Further modelling of the PDR physics and chemistry may help to pinpoint the physical and chemical processes that drive these spatial-spectral variations in the Orion Bar. Furthermore, we expect that similar studies of a variety of sources with JWST will provide deeper insight in the origin of the A, B, C, and D classes identified by using measurements obtained with ISO and Spitzer (Peeters et al. 2002a; van Diedenhoven et al. 2004; Sloan et al. 2014; Matsuura et al. 2014). In any case, the spatial-spectral variations in the Orion Bar provide a framework in which AIB spectra of extragalactic regions of massive star formation can be analysed in terms of the physical conditions in their PDRs.  \nAn analysis of the Spitzer spectra of a variety of objects revealed that the mid-IR spectra at the brightest spots in PDRs show remarkably similar AIBs and this has been taken to imply that in the harsh conditions of these positions, the PAH family is dominated by a few species that can withstand the harsh conditions in PDRs (Andrews et al. 2015; Tielens 2013). For now, the limited diversity in the AIB characteristics of the Orion Bar templates points in the same direction. Indeed, a broad distribution of PAHs would result in much more sub-structure and variation behavior. In addition, the disappearance of the (weak) 11.25 µ m component of the 11.2 µ m AIB in the PDR surface layers implies that photochemistry is important: only the most robust species survive in the harsh conditions at the surface of the PDR. We note that while the cosmic AIB emission can be classified in four classes (A, B, C, and D), interstellar AIB emission invariably belongs to class A. The Orion Bar spectrum is the "poster child" of the class A spectrum - out of all the AIBs found in the Orion Bar template spectra, only the 11.2 µ m AIB shows some indication of a class B contribution. Only very distinctly di ff erent classes of objects with unique histories display classes B, C, and D AIB emission. This too implies that the interstellar PAH family consists of a small set of very robust species.  \nMoreover, we conclude that the profiles of the 5.25, 6.2, and 11.2 µ m AIBs are controlled by anharmonicity, rather than by blending of a number of bands, while variations in the widths of these bands in the di ff erent template spectra are related to variations in the excitation of the emitting species in those positions. As a corollary, this implies that these bands are likely dominated by emission of a single carrier, further supporting the GrandPAH hypothesis (Mackie et al. 2022). However, it remains to be seen whether this spectral similarity still holds when a much larger range of objects is investigated at the higher spectral resolution of JWST .  \nThe much higher spatial resolution of JWST provides further insight in the processes that might be relevant for the composition of the PAH family. Specifically, as argued in Sect. 4, the decreased prominence of the minor features in the 5.75, 7.7, and 11.2 µ m AIBs in the atomic region indicates the loss of a subpopulation of the PAH family (Berné & Tielens 2012; Montillaud et al. 2013; Berné et al. 2015; Andrews et al. 2016) or the loss of very small grains (Pilleri et al. 2012, 2015) in this region. Likewise, the decrease in the methyl group coverage, as evidenced in the variation of the 3.4 / 3.3 AIB ratio, indicates the loss of more loosely bound functional groups in the surface layers or their conversion to aromatic moieties. Models suggest that photolysis of PAHs is controlled by the strength of the UV field over the atomic hydrogen density (Montillaud et al. 2013; Andrews et al. 2016; Berné et al. 2015) and over much of the Orion Bar PDR, G 0 / n (H) ≃ 1 (Tielens & Hollenbach 1985; BernardSalas et al. 2012b), implying that the PAHs that are lost from the PAH family are small ( N C ≲ 50).  \nHence, the picture that emerges from the analysis of the template spectra is the increased importance of FUV processing in PDR surface layers, resulting in a \'weeding out\' of the weakest links of the PAH family. These less resistant species are possibly formed from small hydrocarbons during the preceding dark cloud phase of the region in a bottom-up chemical scenario (Cuadrado et al. 2015; McGuire et al. 2021), feeding the gas with small hydrogen rich photolytically produced species (Alata et al. 2015). The UV processing of the PAH family will start with the loss of the smallest PAHs but as the PDR material is advected to the surface, larger and larger PAHs will become susceptible to photoprocessing. Similar scenarios have been proposed recently with data not obtained with JWST , albeit with slightly di ff erent numbers (e.g. N C ≲ 50 in Murga et al. 2022). In favourable conditions, the processing of very large PAHs ( N C ≳ 60) may lead to the formation of C60 (Berné & Tielens 2012; Berné et al. 2015). While the present data around 18.9 µ mare still marred by instrumental artefacts, future searches for the signature of fullerenes in the MIRI spectra of these surface layers may reveal whether such a scenario plays a role under the conditions of the Orion Bar.  \nJWST is poised to obtain high-quality spectral imaging observations of a large sample of environments probing the full range of physical conditions that are relevant for AIB emission. These observations promise to capture the complexity of AIB emission with the unprecedented detail that is needed in order to advance our understanding of the photochemical evolution of large carbonaceous molecules.  \nAcknowledgements. Wethank the referee, Alan Tokunaga, for constructive comments on the manuscript. We are very grateful to the JWST Help Desk for their support with pipeline and calibration issues that were encountered while writing this paper. This work is based on observations made with the NASA / ESA / CSA James Webb Space Telescope. The data were obtained from the Mikulski Archive for Space Telescopes at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-03127 for JWST. These observations are associated with program #1288. Support for program #1288 was provided by NASA  \nthrough a grant from the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-03127. EP and JC acknowledge support from the University of Western Ontario, the Institute for Earth and Space Exploration, the Canadian Space Agency (CSA; 22JWGO1- 16), and the Natural Sciences and Engineering Research Council of Canada. Studies of interstellar PAHs at Leiden Observatory (AT) are supported by a Spinoza premie from the Dutch Science Agency, NWO. CB is grateful for an appointment at NASA Ames Research Center through the San José State University Research Foundation (80NSSC22M0107). Part of this work was supported by the Programme National \'Physique et Chimie du Milieu Interstellaire\' (PCMI) of CNRS / INSU with INC / INP co-funded by CEA and CNES. JRG and SC thank the Spanish MCINN for funding support under grant PID2019-106110GB-I00. TO is supported by JSPS Bilateral Program, Grant Number 120219939. Work by YO and MR is carried out within the Collaborative Research Centre 956, sub-project C1, funded by the Deutsche Forschungsgemeinschaft (DFG) - project ID 184018867. This work was also partly supported by the Spanish program Unidad de Excelencia María de Maeztu CEX2020001058-M, financed by MCIN / AEI / 10.13039 / 501100011033. NN is funded by the United Arab Emirates University (UAEU) through UAEU Program for Advanced Research (UPAR) grant G00003479. AP (Amit Pathak) would like to acknowledge financial support from Department of Science and Technology SERB via Core Research Grant (DST-CRG) grant (SERB-CRG / 2021 / 000907), Institutes of Eminence (IoE) incentive grant, BHU (incentive / 2021- 22 / 32439), Banaras Hindu University, Varanasi and thanks the Inter-University Centre for Astronomy and Astrophysics, Pune for associateship. This work is sponsored (in part) by the CAS, through a grant to the CAS South America Center for Astronomy (CASSACA) in Santiago, Chile. AR gratefully acknowledges support from the directed Work Package at NASA Ames titled: \'Laboratory Astrophysics -The NASA Ames PAH IR Spectroscopic Database\'. MB acknowledges DST for the DST INSPIRE Faculty fellowship. HZ acknowledges support from the Swedish Research Council (contract No 2020-03437). P.M. acknowledges grants EUR2021-122006, TED2021-129416A-I00 and PID2021-125309OA-I00 funded by MCIN / AEI / 10.13039 / 501100011033 and European Union NextGenerationEU / PRTR. MSK is funded by RSCF, grant number 21-12-00373.', 'References': "Alata, I., Jallat, A., Gavilan, L., et al. 2015, A&A, 584, A123 Allain, T., Leach, S., & Sedlmayr, E. 1996a, A&A, 305, 602 Allain, T., Leach, S., & Sedlmayr, E. 1996b, A&A, 305, 616 Allamandola, L. J., Bregman, J. D., Sandford, S. A., et al. 1989a, ApJ, 345, L59 Allamandola, L. J., Hudgins, D. M., & Sandford, S. A. 1999, ApJ, 511, L115 Allamandola, L. J., Tielens, A. G. G. M., & Barker, J. R. 1985, ApJ, 290, L25 Allamandola, L. J., Tielens, A. G. G. M., & Barker, J. R. 1989b, ApJS, 71, 733 Andrews, H., Boersma, C., Werner, M. W., et al. 2015, ApJ, 807, 99 Andrews, H., Candian, A., & Tielens, A. G. G. M. 2016, A&A, 595, A23 Argyriou, I., Glasse, A., Law, D. R., et al. 2023, arXiv e-prints,  \n- arXiv:2303.13469 Argyriou, I., Wells, M., Glasse, A., et al. 2020, A&A, 641, A150 Arnoult, K. M., Wdowiak, T. J., & Beegle, L. W. 2000, ApJ, 535, 815 Bakes, E. L. O., Tielens, A. G. G. M., & Bauschlicher, C. W. 2001, ApJ, 556, 501 Bauschlicher, C. W., Peeters, E., & Allamandola, L. J. 2009, ApJ, 697, 311 Bauschlicher, C. W. J. & Ricca, A. 2014, Theoretical Chemistry Accounts, 133, 1454 Bauschlicher, Charles W., J., Peeters, E., & Allamandola, L. J. 2008, ApJ, 678, 316 Beintema, D. A., van den Ancker, M. E., Molster, F. J., et al. 1996, A&A, 315, L369 Bernard-Salas, J., Habart, E., Arab, H., et al. 2012a, A&A, 538, A37 Bernard-Salas, J., Habart, E., Arab, H., et al. 2012b, A&A, 538, A37 Berné, O., Habart, É., Peeters, E., et al. 2022, PASP, 134, 054301 Berné, O., Joblin, C., Deville, Y., et al. 2007, A&A, 469, 575 Berné, O., Montillaud, J., & Joblin, C. 2015, A&A, 577, A133 Berné, O. & Tielens, A. G. G. M. 2012, Proceedings of the National Academy of Science, 109, 401 Bernstein, M. P., Sandford, S. A., & Allamandola, L. J. 1996, ApJ, 472, L127 Boersma, C., Bregman, J., & Allamandola, L. J. 2014, ApJ, 795, 110 Boersma, C., Bregman, J., & Allamandola, L. J. 2018, ApJ, 858, 67 Boersma, C., Bregman, J. D., & Allamandola, L. J. 2013, ApJ, 769, 117 Boersma, C., Mattioda, A. L., Bauschlicher, C. W., et al. 2009a, ApJ, 694, 704 Boersma, C., Mattioda, A. L., Bauschlicher, C. W., J., et al. 2009b, ApJ, 690, 1208  \nBoersma, C., Rubin, R. H., & Allamandola, L. J. 2012, ApJ, 753, 168 Bohlin, R. C., Gordon, K. D., & Tremblay, P. E. 2014, PASP, 126, 711 Böker, T., Arribas, S., Lützgendorf, N., et al. 2022, A&A, 661, A82 Boulanger, F., Reach, W. T., Abergel, A., et al. 1996, A&A, 315, L325  \n- Bouwman, J., Henning, T., Hillenbrand, L. A., et al. 2008, ApJ, 683, 479\n- Bregman, J. 1989, in IAU Symp. 135: Interstellar Dust, Vol. 135, 109\n- Bregman, J. & Temi, P. 2005, ApJ, 621, 831\n- Bregman, J. D., Allamandola, L. J., Witteborn, F. C., Tielens, A. G. G. M., & Geballe, T. R. 1989, ApJ, 344, 791\n- Bullins, K. W., Huang, T. T. S., & Kirkby, S. J. 2009, International Journal of Quantum Chemistry, 109, 1322  \nBuragohain, M., Pathak, A., Sakon, I., & Onaka, T. 2020, ApJ, 892, 11  \nCami, J., Bernard-Salas, J., Peeters, E., & Malek, S. E. 2010, Science, 329, 1180  \nCandian, A., Kerr, T. H., Song, I. O., McCombie, J., & Sarre, P. J. 2012, MNRAS,  \n426, 389  \nCandian, A. & Sarre, P. J. 2015, MNRAS, 448, 2960 Carpentier, Y., Féraud, G., Dartois, E., et al. 2012, A&A, 548, A40 Cernicharo, J., Agúndez, M., Cabezas, C., et al. 2021, A&A, 649, L15 Cohen, M., Tielens, A. G. G. M., Bregman, J., et al. 1989, ApJ, 341, 246 Croiset, B. A., Candian, A., Berné, O., & Tielens, A. G. G. M. 2016, A&A, 590,  \n- A26\n- Cuadrado, S., Goicoechea, J. R., Pilleri, P., et al. 2015, A&A, 575, A82 Dartois, E., Charon, E., Engrand, C., Pino, T., & Sandt, C. 2020, A&A, 637, A82 Dartois, E., Muñoz Caro, G. M., Debo ffl e, D., & d'Hendecourt, L. 2004, A&A, 423, L33\n- Dartois, E., Muñoz Caro, G. M., Debo ffl e, D., Montagnac, G., & D'Hendecourt, L. 2005, A&A, 432, 895  \nde Graauw, T., Haser, L. N., Beintema, D. A., et al. 1996, A&A, 315, L49  \n- Dubosq, C., Pla, P., Dartois, E., & Simon, A. 2023, A&A, 670, A175\n- Esquej, P., Alonso-Herrero, A., González-Martín, O., et al. 2014, ApJ, 780, 86  \nFoschino, S., Berné, O., & Joblin, C. 2019, A&A, 632, A84  \n- Galliano, F., Madden, S. C., Tielens, A. G. G. M., Peeters, E., & Jones, A. P. 2008, ApJ, 679, 310  \nGalué, H. Á. 2014, Chem. Sci., 5, 2667 Genzel, R., Lutz, D., Sturm, E., et al. 1998, ApJ, 498, 579 Gillett, F. C., Forrest, W. J., & Merrill, K. M. 1973, ApJ, 183, 87 Godard, M., Féraud, G., Chabot, M., et al. 2011, A&A, 529, A146 Goicoechea, J. R., Pety, J., Cuadrado, S., et al. 2016, Nature, 537, 207 Gordon, K. D., Bohlin, R., Sloan, G. C., et al. 2022, AJ, 163, 267 Habart, É., Peeters, E., Berné, O., et al. 2023, A&A, 000, 000 Hony, S., Van Kerckhoven, C., Peeters, E., et al. 2001, A&A, 370, 1030 Houck, J. R., Roellig, T. L., van Cleve, J., et al. 2004, ApJS, 154, 18 Hudgins, D. M. & Allamandola, L. J. 1999, ApJ, 516, L41  \n- Hudgins, D. M., Bauschlicher, Charles W., J., & Allamandola, L. J. 2005, ApJ, 632, 316  \nJensen, J. J., Hönig, S. F., Rakshit, S., et al. 2017, MNRAS, 470, 3071  \n- Joalland, B., Simon, A., Marsden, C. J., & Joblin, C. 2009, A&A, 494, 969\n- Joblin, C., Boissel, P., Leger, A., D'Hendecourt, L., & Defourneau, D. 1995, A&A, 299, 835\n- Joblin, C., Bron, E., Pinto, C., et al. 2018, A&A, 615, A129\n- Joblin, C., Szczerba, R., Berné, O., & Szyszka, C. 2008, A&A, 490, 189 Joblin, C., Tielens, A. G. G. M., Allamandola, L. J., & Geballe, T. R. 1996, ApJ, 458, 610\n- Jochims, H. W., Baumgärtel, H., & Leach, S. 1999, ApJ, 512, 500\n- Jones, A. P., Fanciullo, L., Köhler, M., et al. 2013, A&A, 558, A62\n- Kessler, M. F., Steinz, J. A., Anderegg, M. E., et al. 1996, A&A, 315, L27\n- Knight, C., Peeters, E., Stock, D. J., Vacca, W. D., & Tielens, A. G. G. M. 2021, ApJ, 918, 8\n- Knight, C., Peeters, E., Tielens, A. G. G. M., & Vacca, W. D. 2022a, MNRAS, 509, 3523\n- Knight, C., Peeters, E., Wolfire, M., & Stock, D. J. 2022b, MNRAS, 510, 4888 Léger, A. & Puget, J. L. 1984, A&A, 137, L5\n- Lemmens, A. K., Rap, D. B., Thunnissen, J. M. M., et al. 2020, The Journal of Physical Chemistry Letters, 11, 8997, pMID: 33035060\n- Lemmens, A. K., Rap, D. B., Thunnissen, J. M. M., et al. 2019, A&A, 628, A130 Lemmens, A. K., Rijs, A. M., & Buma, W. J. 2021, ApJ, 923, 238\n- Lutz, D., Spoon, H. W. W., Rigopoulou, D., Moorwood, A. F. M., & Genzel, R. 1998, ApJ, 505, L103\n- Mackie, C. J., Candian, A., Huang, X., et al. 2015, J. Chem. Phys., 143, 224314 Mackie, C. J., Candian, A., Huang, X., et al. 2018, Physical Chemistry Chemical Physics (Incorporating Faraday Transactions), 20, 1189\n- Mackie, C. J., Candian, A., Huang, X., et al. 2016, Journal of Chemical Physics, 145, 084313\n- Mackie, C. J., Candian, A., Lee, T. J., & Tielens, A. G. G. M. 2022, The Journal of Physical Chemistry A, 126, 3198, pMID: 35544706\n- Maltseva, E., Mackie, C. J., Candian, A., et al. 2018, Astron. Astrophys., 610, A65\n- Maltseva, E., Petrignani, A., Candian, A., et al. 2015, The Astrophysical Journal, 814, 23\n- Maltseva, E., Petrignani, A., Candian, A., et al. 2016, The Astrophysical Journal, 831, 58\n- Maragkoudakis, A., Peeters, E., & Ricca, A. 2020, MNRAS, 494, 642\n- Marconi, A., Testi, L., Natta, A., & Walmsley, C. M. 1998, A&A, 330, 696  \n- Matsuura, M., Bernard-Salas, J., Lloyd Evans, T., et al. 2014, MNRAS, 439, 1472\n- Matthew J. Shannon & Christiaan Boersma. 2018, in Proceedings of the 17th Python in Science Conference, ed. Fatih Akici, David Lippa, Dillon Niederhut, & M. Pacer, 99 - 105  \nMcGuire, B. A., Burkhardt, A. M., Kalenskii, S., et al. 2018, Science, 359, 202 McGuire, B. A., Loomis, R. A., Burkhardt, A. M., et al. 2021, Science, 371, 1265 Meeus, G., Waters, L. B. F. M., Bouwman, J., et al. 2001, A&A, 365, 476  \n- Menten, K. M., Reid, M. J., Forbrich, J., & Brunthaler, A. 2007, A&A, 474, 515 Montillaud, J., Joblin, C., & Toublanc, D. 2013, A&A, 552, A15\n- Mori, T. I., Sakon, I., Onaka, T., et al. 2012, ApJ, 744, 68\n- Murga, M. S., Kirsanova, M. S., Wiebe, D. S., & Boley, P. A. 2022, MNRAS, 509, 800  \nPabst, C., Higgins, R., Goicoechea, J. R., et al. 2019, Nature, 565, 618  \nParmar, P. S., Lacy, J. H., & Achtermann, J. M. 1991, ApJ, 372, L25 Pech, C., Joblin, C., & Boissel, P. 2002, A&A, 388, 639  \nPeeters, E., Allamandola, L. J., Hudgins, D. M., Hony, S., & Tielens, A. G. G. M.  \n2004a, in Astronomical Society of the Pacific Conference Series, Vol. 309,  \nAstrophysics of Dust, ed. A. N. Witt, G. C. Clayton, & B. T. Draine, 141  \n- Peeters, E., Bauschlicher, Charles W., J., Allamandola, L. J., et al. 2017, ApJ, 836, 198\n- Peeters, E., Berné, O., Habart, É., et al. 2023, A&A, 000, 000\n- Peeters, E., Hony, S., Van Kerckhoven, C., et al. 2002a, A&A, 390, 1089  \nPeeters, E., Martín-Hernández, N. L., Damour, F., et al. 2002b, A&A, 381, 571  \n- Peeters, E., Spoon, H. W. W., & Tielens, A. G. G. M. 2004b, ApJ, 613, 986\n- Peeters, E., Tielens, A. G. G. M., Allamandola, L. J., & Wolfire, M. G. 2012, ApJ, 747, 44  \nPilleri, P., Joblin, C., Boulanger, F., & Onaka, T. 2015, A&A, 577, A16  \nPilleri, P., Montillaud, J., Berné, O., & Joblin, C. 2012, A&A, 542, A69  \nPino, T., Dartois, E., Cao, A. T., et al. 2008, A&A, 490, 665  \n- Pla, P., Wang, Y., Martín, F., & Alcamí, M. 2020, ApJ, 899, 18\n- Rap, D. B., Schrauwen, J. G. M., Marimuthu, A. N., Redlich, B., & Brünken, S. 2022, Nature Astronomy, 6, 1059\n- Ricca, A., Bauschlicher, Charles W., J., Boersma, C., Tielens, A. G. G. M., & Allamandola, L. J. 2012, ApJ, 754, 75  \nRicca, A., Boersma, C., & Peeters, E. 2021, ApJ, 923, 202  \nRobertson, J. 1986, Advances in Physics, 35, 317  \nRoche, P. F., Aitken, D. K., & Smith, C. H. 1989, MNRAS, 236, 485  \n- Rosenberg, M. J. F., Berné, O., Boersma, C., Allamandola, L. J., & Tielens, A. G. G. M. 2011, A&A, 532, A128\n- Sandstrom, K. M., Bolatto, A. D., Draine, B. T., Bot, C., & Stanimirovi'c, S. 2010, ApJ, 715, 701  \nSarre, P. J. 2019, MNRAS, 490, L17  \n- Schirmer, T., Abergel, A., Verstraete, L., et al. 2020, Astronomy & Astrophysics, 639, A144\n- Schirmer, T., Ysard, N., Habart, E., et al. 2022, A&A, 666, A49\n- Schutte, W. A., Tielens, A. G. G. M., & Allamandola, L. J. 1993, ApJ, 415, 397 Sellgren, K., Werner, M. W., Ingalls, J. G., et al. 2010, ApJ, 722, L54\n- Shannon, M. J. & Boersma, C. 2019, ApJ, 871, 124\n- Shannon, M. J., Stock, D. J., & Peeters, E. 2016, ApJ, 824, 111\n- Sidhu, A., Tielens, A. G. G. M., Peeters, E., & Cami, J. 2022, MNRAS, 514, 342 Sloan, G. C., Bregman, J. D., Geballe, T. R., Allamandola, L. J., & Woodward, E. 1997, ApJ, 474, 735\n- Sloan, G. C., Lagadec, E., Zijlstra, A. A., et al. 2014, ApJ, 791, 28\n- Smith, J. D. T., Draine, B. T., Dale, D. A., et al. 2007, ApJ, 656, 770\n- Stock, D. J., Choi, W. D. Y., Moya, L. G. V., et al. 2016, ApJ, 819, 65\n- Stock, D. J. & Peeters, E. 2017, ApJ, 837, 129\n- Tauber, J. A., Tielens, A. G. G. M., Meixner, M., & Goldsmith, P. F. 1994, ApJ, 422, 136\n- Tielens, A. G. G. M. 2008, ARA&A, 46, 289\n- Tielens, A. G. G. M. 2013, Reviews of Modern Physics, 85, 1021\n- Tielens, A. G. G. M. 2021, Molecular Astrophysics (Cambridge: Cambridge University Press)\n- Tielens, A. G. G. M. & Hollenbach, D. 1985, ApJ, 291, 722\n- Tokunaga, A. T. & Bernstein, L. S. 2021, ApJ, 916, 52\n- Tokunaga, A. T., Bernstein, L. S., & Connelley, M. 2022, ApJ, 939, 86\n- Tokunaga, A. T., Sellgren, K., Smith, R. G., et al. 1991, ApJ, 380, 452\n- Tsuge, M., Tseng, C.-Y., & Lee, Y.-P. 2018, Phys. Chem. Chem. Phys., 20, 5344 Van De Putte, D., Trahin, B., Habart, É., Peeters, E., & Berné, O. 2023, A&A, 000, 000\n- van Diedenhoven, B., Peeters, E., Van Kerckhoven, C., et al. 2004, ApJ, 611, 928 Vermeij, R., Peeters, E., Tielens, A. G. G. M., & van der Hulst, J. M. 2002, A&A, 382, 1042  \nVerstraete, L., Pech, C., Moutou, C., et al. 2001, A&A, 372, 981  \n- Vicente, S., Berné, O., Tielens, A. G. G. M., et al. 2013, ApJ, 765, L38\n- Wells, M., Pel, J. W., Glasse, A., et al. 2015, PASP, 127, 646\n- Wenzel, G., Joblin, C., Giuliani, A., et al. 2020, A&A, 641, A98\n- Wenzel, G., Simon, A., Banhatti, S., et al. 2022, Journal of Molecular Spectroscopy, 385, 111620\n- Werner, M. W., Roellig, T. L., Low, F. J., et al. 2004, ApJS, 154, 1  \nArticle number, page 16 of 25  \nWexler, A. S. 1967, Applied Spectroscopy Reviews, 1, 29 Yan, L., Chary, R., Armus, L., et al. 2005, ApJ, 628, 604 Yang, X. J., Glaser, R., Li, A., & Zhong, J. X. 2016, MNRAS, 462, 1551 Young Owl, R. C., Meixner, M. M., Wolfire, M., Tielens, A. G. G. M., & Tauber, J. 2000, ApJ, 540, 886 Zhen, J., Castellanos, P., Linnartz, H., & Tielens, A. 2016, Molecular Astrophysics, 5, 1 Zhen, J., Castellanos, P., Paardekooper, D. M., Linnartz, H., & Tielens, A. G. G. M. 2014a, ApJ, 797, L30 Zhen, J., Paardekooper, D. M., Candian, A., Linnartz, H., & Tielens, A. G. G. M.  \n- 2014b, Chemical Physics Letters, 592, 211\n- 1 Department of Physics & Astronomy, The University of Western Ontario, London ON N6A 3K7, Canada; e-mail: rchown3@uwo.ca\n- 2 Institute for Earth and Space Exploration, The University of Western Ontario, London ON N6A 3K7, Canada\n- 3 Carl Sagan Center, SETI Institute, 339 Bernardo Avenue, Suite 200, Mountain View, CA 94043, USA\n- 4 Leiden Observatory, Leiden University, P.O. Box 9513, 2300 RA Leiden, The Netherlands\n- 5 Astronomy Department, University of Maryland, College Park, MD20742, USA\n- 6 Institut de Recherche en Astrophysique et Planétologie, Université Toulouse III - Paul Sabatier, CNRS, CNES, 9 Av. du colonel Roche, 31028 Toulouse Cedex 04, France\n- 7 Institut d'Astrophysique Spatiale, Université Paris-Saclay, CNRS, Bˆ a timent 121, 91405 Orsay Cedex, France\n- 8 Department of Astronomy, University of Michigan, 1085 South University Avenue, Ann Arbor, MI 48109, USA\n- 9 Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD, 21218, USA\n- 10 ACRI-ST, Centre d'Etudes et de Recherche de Grasse (CERGA), 10 Av. Nicolas Copernic, F-06130 Grasse, France\n- 11 INCLASS Common Laboratory., 10 Av. Nicolas Copernic, 06130 Grasse, France\n- 12 NASA Ames Research Center, MS 245-6, Mo ff ett Field, CA 94035-1000, USA\n- 13 LERMA, Observatoire de Paris, PSL Research University, CNRS, Sorbonne Universités, F-92190 Meudon, France\n- 14 Instituto de Física Fundamental (CSIC), Calle Serrano 121-123, 28006, Madrid, Spain\n- 15 Institut des Sciences Moléculaires d'Orsay, CNRS, Université Paris-Saclay, Bˆ a timent 520, 91405 Orsay Cedex, France\n- 16 UKAstronomy Technology Centre, Royal Observatory Edinburgh, Blackford Hill, Edinburgh EH9 3HJ, UK\n- 17 Observatorio Astronómico Nacional (OAN,IGN), Alfonso XII, 3, E-28014 Madrid, Spain  \n18  \nSterrenkundig Observatorium, Universiteit Gent, Gent, Belgium  \n- 19 KU Leuven Quantum Solid State Physics (QSP), Celestijnenlaan 200d - box 2414, 3001 Leuven, Belgium\n- 20 Institut de Planétologie et d'Astrophysique de Grenoble (IPAG), Université Grenoble Alpes, CNRS, F-38000 Grenoble, France\n- 21 Institut de Radioastronomie Millimétrique (IRAM), 300 Rue de la Piscine, F-38406 Saint-Martin d'Hères, France\n- 22 I. Physikalisches Institut der Universität zu Köln, Zülpicher Straße 77, 50937 Köln, Germany\n- 23 Department of Astronomy, Graduate School of Science, The University of Tokyo, 7-3-1 Bunkyo-ku, Tokyo 113-0033, Japan\n- 24 Department of Physics, Faculty of Science and Engineering, Meisei University, 2-1-1 Hodokubo, Hino, Tokyo 191-8506, Japan\n- 25 Johns Hopkins University, 3400 N. Charles Street, Baltimore, MD, 21218, USA\n- 26 Physikalischer Verein - Gesellschaft für Bildung und Wissenschaft, Robert-Mayer-Straße 2, 60325 Frankfurt am Main, Germany\n- 27 Institut für Angewandte Physik. Goethe-Universität Frankfurt, Max-von-Laue-Str. 1, 60438 Frankfurt am Main, Germany\n- 28 Department of Space, Earth and Environment, Chalmers University of Technology, Onsala Space Observatory, SE-439 92 Onsala, Sweden  \n- 29 Instituto de Astrofísica e Ciências do Espaço, Tapada da Ajuda, Edifício Leste, 2 · Piso, P-1349-018 Lisboa, Portugal\n- 30 Instituto de Física e Química, Universidade Federal de Itajubá, Av. BPS 1303, Pinheirinho, 37500-903, Itajubá, MG, Brazil\n- 31 Bay Area Environmental Research Institute, Mo ff ett Field, CA 94035, USA\n- 32 Australian Synchrotron, Australian Nuclear Science and Technology Organisation (ANSTO), Victoria, Australia\n- 33 INAF - Osservatorio Astrofisico di Catania, Via Santa Sofia 78, 95123 Catania, Italy\n- 34 Laboratoire de Physique de l'École Normale Supérieure, ENS, Université PSL, CNRS, Sorbonne Université, Université de Paris, 75005, Paris, France\n- 35 Laboratory for Atmospheric and Space Physics, University of Colorado, Boulder, CO 80303, USA\n- 36 Department of Chemistry, University of Colorado, Boulder, CO 80309, USA\n- 37 Institute for Modeling Plasma, Atmospheres, and Cosmic Dust (IMPACT), University of Colorado, Boulder, CO 80303, USA\n- 38 Faculty of Aerospace Engineering, Delft University of Technology, Kluyverweg 1, 2629 HS Delft, The Netherlands\n- 39 FELIX Laboratory, Institute for Molecules and Materials, Radboud University, Toernooiveld 7, 6525ED Nijmegen, The Netherlands\n- 40 School of Physics, University of Hyderabad, Hyderabad, Telangana 500046, India\n- 41 Department of Earth, Environment, and Physics Worcester State University, Worcester, MA, USA\n- 42 Anton Pannekoek Institute for Astronomy (API), University of Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands\n- 43 Delft University of Technology, Delft, The Netherlands\n- 44 Laboratoire de Physique des deux infinis Irène Joliot-Curie, Université Paris-Saclay, CNRS / IN2P3, Bˆ a timent 104, 91405 Orsay Cedex, France\n- 45 Department of Chemistry, GITAM school of Science, GITAM Deemed to be University, Bangalore, India\n- 46 Institut de Physique de Rennes, UMR CNRS 6251, Université de Rennes 1, Campus de Beaulieu, 35042 Rennes Cedex, France\n- 47 Department of Chemistry, The University of British Columbia, Vancouver, British Columbia, Canada\n- 48 National Radio Astronomy Observatory (NRAO), 520 Edgemont Road, Charlottesville, VA 22903, USA\n- 49 European Space Astronomy Centre (ESAC / ESA), Villanueva de la Cañada, E-28692 Madrid, Spain\n- 50 Observatoire de Paris, PSL University, Sorbonne Université, LERMA, 75014, Paris, France\n- 51 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge MA 02138, USA\n- 52 Sorbonne Université, CNRS, UMR 7095, Institut d'Astrophysique de Paris, 98bis bd Arago, 75014 Paris, France\n- 53 Institut Universitaire de France, Ministère de l'Enseignement Supérieur et de la Recherche, 1 rue Descartes, 75231 Paris Cedex 05, France\n- 54 Department of Physics and Astronomy, Rice University, Houston TX, 77005-1892, USA\n- 55 Yunnan Observatories, Chinese Academy of Sciences, 396 Yangfangwang, Guandu District, Kunming, 650216, China\n- 56 Chinese Academy of Sciences South America Center for Astronomy, National Astronomical Observatories, CAS, Beijing 100101, China\n- 57 Departamento de Astronomía, Universidad de Chile, Casilla 36-D, Santiago, Chile\n- 58 Departments of Chemistry and Astronomy, University of Virginia, Charlottesville, Virginia 22904, USA\n- 59 InterCat and Dept. Physics and Astron., Aarhus University, Ny Munkegade 120, 8000 Aarhus C, Denmark\n- 60 Laboratory Astrophysics Group of the Max Planck Institute for Astronomy at the Friedrich Schiller University Jena, Institute of Solid State Physics, Helmholtzweg 3, 07743 Jena, Germany\n- 61 Instituto de Astronomia, Geofísica e Ciências Atmosféricas, Universidade de São Paulo, 05509-090 São Paulo, SP, Brazil\n- 62 Department of Physics and Astronomy, San José State University, San Jose, CA 95192, USA\n- 63 Institut de Ciencies de l'Espai (ICE, CSIC), Can Magrans, s / n, E08193 Bellaterra, Barcelona, Spain\n- 64 ICREA, Pg. Lluìs Companys 23, E-08010 Barcelona, Spain\n- 65 Institut d'Estudis Espacials de Catalunya (IEEC), E-08034 Barcelona, Spain\n- 66 European Space Agency, Space Telescope Science Institute, 3700 San Martin Drive, Baltimore MD 21218, USA\n- 67 Institute of Astronomy, Russian Academy of Sciences, 119017, Pyatnitskaya str., 48 , Moscow, Russia\n- 68 Department of Earth, Ocean, & Atmospheric Sciences, University of British Columbia, Canada V6T 1Z4\n- 69 Telespazio UK for ESA, ESAC, E-28692 Villanueva de la Cañada, Madrid, Spain\n- 70 IPAC, California Institute of Technology, Pasadena, CA, USA\n- 71 Department of Physics and Astronomy, University of Missouri, Columbia, MO 65211, USA\n- 72 Max Planck Institute for Astronomy, Königstuhl 17, 69117 Heidelberg, Germany\n- 73 Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California, USA\n- 74 Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California - Berkeley, Berkeley, California, USA\n- 75 AIM, CEA, CNRS, Université Paris-Saclay, Université Paris Diderot, Sorbonne Paris Cité, 91191 Gif-sur-Yvette, France\n- 76 Institut des Sciences Moléculaires, CNRS, Université de Bordeaux, 33405 Talence, France\n- 77 Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139, USA\n- 78 Instituto de Ciencia de Materiales de Madrid (CSIC), Sor Juana Ines de la Cruz 3, E28049, Madrid, Spain\n- 79 Department of Physics, PO Box 64, 00014 University of Helsinki, Finland\n- 80 Steward Observatory, University of Arizona, Tucson, AZ 857210065, USA\n- 81 BoldlyGo Institute, 31 W 34TH ST FL 7 STE 7159, New York, NY 10001\n- 82 Department of Physics, College of Science, United Arab Emirates University (UAEU), Al-Ain, 15551, UAE\n- 83 National Astronomical Observatory of Japan, National Institutes of Natural Science, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan\n- 84 Sorbonne Université, CNRS, UMR 7095, Institut d'Astrophysique de Paris, 98bis bd Arago, 75014 Paris, France\n- 85 California Institute of Technology, IPAC, 770, S. Wilson Ave., Pasadena, CA 91125, USA\n- 86 Department of Physics, Institute of Science, Banaras Hindu University, Varanasi 221005, India\n- 87 University of Central Florida, Orlando, FL 32765\n- 88 Van't Ho ff Institute for Molecular Sciences, University of Amsterdam, PO Box 94157, 1090 GD, Amsterdam, The Netherlands\n- 89 Laboratoire de Chimie et Physique Quantiques LCPQ / FERMI, UMR5626, Université de Toulouse (UPS) and CNRS, Toulouse, France\n- 90 Laboratory Astrophysics Group of the Max Planck Institute for Astronomy at the Friedrich Schiller University Jena, Institute of Solid State Physics, Helmholtzweg 3, 07743 Jena, Germany\n- 91 Instituto de Matemática, Estatística e Física, Universidade Federal do Rio Grande, 96201-900, Rio Grande, RS, Brazil\n- 92 Center for Astrophysics and Space Sciences, Department of Physics, University of California, San Diego, 9500 Gilman Drive, La Jolla, CA 92093, USA\n- 93 School of Chemistry, The University of Nottingham, University Park, Nottingham, NG7 2RD, United Kingdom\n- 94 Astronomy Department, Ohio State University, Columbus, OH 43210 USA\n- 95 Space Science Institute, 4765 Walnut St., R203, Boulder, CO 80301  \n- 96 Department of Physics, Stockholm University, SE-10691 Stockholm, Sweden\n- 97 Department of Physics, Texas State University, San Marcos, TX 78666 USA\n- 98 School of Physics and Astronomy, Sun Yat-sen University, 2 Da Xue Road, Tangjia, Zhuhai 519000, Guangdong Province, China\n- 99 Star and Planet Formation Laboratory, RIKEN Cluster for Pioneering Research, Hirosawa 2-1, Wako, Saitama 351-0198, Japan\n- 100 Institute of Deep Space Sciences, Deep Space Exploration Laboratory, Hefei 230026, China", 'Appendix A: Inventory of AIBs': 'Table A.1 shows the inventory of AIBs identified in this study.  \nTable A.1. Catalogue of AIB features detected in the five template spectra of the Orion Bar. The method used to select these features is described in Sect. 3. The precise peak positions of nominal AIBs - whose wavelengths appear in boldface in this table - are indicated in Col. 3. We note that we converted the positions in wavelength to wavenumber by rounding to the nearest integer in units of cm -1 , and so the precision of the reported wavenumbers do not reflect the precision of the peak position of the AIBs.  \nRyan Chown et al.: PDRs4All IV: the Aromatic Infrared Bands in the Orion Bar  \nTable A.1. continued.', 'A & A proofs: manuscript no. 46662corr': "Notes. Columns: (1) Wavelength of AIB; nominal AIB names are listed in boldface; (2) Wavenumber of AIB; (3) Wavelength of component; (4) Wavenumber of component; (1)-(4) Tentative detection indicated by '?'; (5) New detection; (6) Peak intensity relative to the 3.3, 7.7 or 11.2 peak intensity for bands in the 3-4, 5-10, and 10-15 µ m range, respectively. The categories are vs ( > 65%), s ( > 35%), ms ( > 15%), mw ( > 5%), w ( > 2%), and vw ( < 2%), where v = very, m = moderate, s = strong, and w = weak. The relative peak intensity of AIBs strongly depend on the employed continuum and spectral decomposition; (7) symmetric profile (symm); asymmetric profile (asymm); blended with another band / component at x µ m (blend (x)); change of slope (slope change); (8) vibrational assignment. ( a ) From Peeters et al. (2023). ( b ) Estimated assuming the component is superimposed on other AIB component(s) such as, e.g., a blue or red wing or shoulder, i.e. we assume a local continuum. ( c ) We only list the strongest band in the 15-20 µ mregion as it su ff ers from calibration issues (Appendix B). ( d ) Influenced by artefact(s) (Appendix B).", 'Appendix B: The impact of residual artefacts in the data': "At the time of writing, there are two types of artifacts in the data that sometimes have appearances similar to AIBs. Firstly, the constructive and destructive interference of layers in the MIRI detector arrays leads to periodic variations in the observed signal versus wavelength called 'fringes' (Argyriou et al. 2020). The latest versions of the JWST pipeline and calibration reference files greatly reduce, but do not completely remove all fringes. These residual fringes can sometimes be di ffi cult to distinguish from AIB components.  \nSecondly, at the moment, the absolute flux calibration of JWST (Gordon et al. 2022) is based on a single A star (private communication with the JWST Help Desk). Intrinsic spectral di ff erences between the star that was used for calibration and the true spectrum of a given source can lead to artifacts in the observed spectrum of that source. Some artifacts are more complicated couplings between instrumental e ff ects and incomplete absolute flux calibration.  \nNon-exhaustive list of wavelengths of artifacts:  \n- 1. Excess emission: 14.3, 14.35 µ m\n- 2. Excess absorption: ∼ 12.2, 13.71-13.77, 13.77-13.8, 13.8513.95 µ m\n- 3. Fringes: 10-12 µ m with an amplitude of about 5% (a ff ects red tail of 11.2)  \nWhile future versions of the pipeline and reference files will fix fringes and calibration issues, we have developed a temporary fix for our analysis. We use MIRI-MRS observations of 10 Lac, which is an O9V star from the JWST CALSPEC 7 program (Bohlin et al. 2014). These data were reduced with the same versions of the JWST pipeline and CRDS files that were used when reducing the Orion Bar data. Then we extracted the MIRI spectrum of 10 Lac using apertures whose radii increase linearly with wavelength. A comparison of the extracted spectrum and the CALSPECmodelspectrum of 10 Lac revealed o ff sets in addition to fringes and absorption / emission artifacts. To deal with these o ff sets which vary from subband to subband, we performed linear regression on the extracted flux versus the model flux within a given subband, and then corrected the observed spectrum using the best-fit parameters for that subband. We then compute the ratio of the corrected, observed 10 Lac flux, to the model flux  \nC ν ≡ F 10 Lac , obs . ν F 10 Lac , model , (B.1)  \nν  \nwhere F 10 Lac , obs . ν is the (corrected) observed 10 Lac spectrum, and F 10 Lac , model ν is the model spectrum. Artifacts in the observed 10 Lac spectrum show up as deviations in C away from 1.0.  \nν  \nWe use C ν to correct the Orion template spectra F obs . ,  \nν  \nF corr . ν [MJy sr -1 ] = F obs . ν / GLYPH<0> A chan ., band( C ν -1) + 1 GLYPH<1> , (B.2)  \nwhere A chan ., band is a parameter to adjust the amplitude of the correction for a given channel and subband. We apply corrections to Channel 2 LONG, Channel 3 MEDIUM, and Channel 3 LONG, using A chan ., band = 0 . 5 , 0 . 25 , 1 . 0, respectively (these factors were chosen manually). For all other wavelength ranges, no correction was applied either because no fringes were visible, or  \nthe correction introduced noise, making it di ffi cult to tell if the fringes were actually corrected or simply buried in noise.  \nExamples of uncorrected template spectra are shown in Fig. B.1. Fringes and absorption / emission artifacts in 10 Lac (the green curve) are detected at the same wavelengths in the Orion data. In Fig. B.2 we show the H ii region template spectrum before and after applying the correction described above. The correction e ffi ciently removes fringes (the first three panels), excess absorption (fourth panel), and excess emission artifacts (bottom panel) which would have otherwise been incorrectly classified as AIB emission. We also conclude that the terraces in the 12.7 µ mAIBare real. We note however that some residual calibration issues remain, in particular at the longer wavelengths. We therefore restricted our analysis to wavelengths shorter than 15 µ m with the exception of the addition of the moderately strong 16.4 µ m AIB. We note that an artefact is also present at 18.9 µ m. This wavelength corresponds to the strongest band of C60 emission (Cami et al. 2010; Sellgren et al. 2010). Unfortunately, due to the artefact, we are currently unable to confirm (or refute) the detection of C60 in these templates. In addition, the 12.2 µ m absorption feature is due to a spectral leak (estimated to be ∼ 3% of the 6.1 µ m signal; JWST Help Desk). This leak also a ff ected the calibration star resulting in an incorrect photometric correction in the current pipeline. Due to the di ff erent spectral shape of Orion with respect to the calibration star, this resulted in the observed absorption feature and our applied correction is unable to correct this. This artefact influences the 12.0 and 12.7 µ mAIBs.", 'Appendix C: Continuum determination': 'Here we show the continuum curves and the anchor points that were used to estimate them. Figure C.1 shows the MIRI MRS template spectra and continua, while Figure C.2 shows the same but for NIRSpec data.  \nFor the 3 µ mregion (NIRSpec templates), we fit a linear continuum to the wavelength ranges [3.05, 3.07], [3.667, 3.689], and [3.7068, 3.720], chosen to avoid emission lines. For the MIR (MIRI templates), we compute a cubic spline interpolation (tension of 1.0, IDL) anchored at selected wavelengths to estimate the continuum emission. We selected anchor points that avoid all (weak) AIBs in all five templates. This resulted in anchor points at 5.011, 5.042, 5.073, 5.487, 5.948, 6.576, 6.785, 7.013, 9.131, 9.440, 10.191, 10.344, 14.642, 14.913, 15.178, 15.461, 19.349, and 19.632 µ m. The resulting continuum is shown in Fig. C.1. Our choice to use the same anchor points for all templates leads to slightly-overestimated continuum at 5.5 µ m in the DF 3 and at 7.1 µ m in the H ii region. Subtracting slightly-overestimated continuum from the templates will a ff ect the resulting AIB emission in some cases.  \nFig. B.1. Zoom-ins on continuum-subtracted template spectra (not corrected for calibration issues, i.e. F obs . ν in Eq. B.2), compared with the MIRI-MRS spectrum of the absolute calibration star 10 Lac divided by the model for this star (green). The positions of identified AIB features are shown with vertical lines. The bracketed regions indicate excess absorption or emission artifacts due to incomplete flux calibration (Appendix B).  \n<!-- image -->  \nFig. B.2. F corr . ν (black) and F obs . ν (blue) for the H ii region template spectrum, showing improvements from the correction (Eq. B.2). The positions of identified AIB features are shown with vertical lines. The bracketed regions indicate excess absorption or emission artifacts due to incomplete flux calibration (Appendix B).  \n<!-- image -->  \nFig. C.1. The MIRI-MRS template spectra with their respective continua. The anchor points for each continuum are represented by filled circles. See Sect. 2.3 and Appendix C for details.  \n<!-- image -->  \nFig. C.2. The NIRSpec template spectra (grating g395h, filter f290lp) from Peeters et al. (2023) along with their respective continua, zoomed in on the wavelength range that is relevant to the 3.3 µ m AIB. See Sect. 2.3 and Appendix C for details.  \n<!-- image -->'}