PDRs4All. V. Modelling the dust evolution across the illuminated edge of the Orion Bar

We study the emission of dust grains within the Orion Bar - a well-known, highly far-UV (FUV)-irradiated PDR. The Orion Bar because of its edge-on geometry provides an exceptional benchmark for characterizing dust evolution and the associated driving processes under varying physical conditions. Our goal is to constrain the local properties of dust by comparing its emission to models. Taking advantage of the recent JWST PDRs4All data, we follow the dust emission as traced by JWST NIRCam (at 3.35 and 4.8 micron) and MIRI (at 7.7, 11.3, 15.0, and 25.5 micron), along with NIRSpec and MRS spectroscopic observations. First, we constrain the minimum size and hydrogen content of carbon nano-grains from a comparison between the observed dust emission spectra and the predictions of the THEMIS dust model coupled to the numerical code DustEM. Using this dust model, we then perform 3D radiative transfer simulations of dust emission with the SOC code and compare to data obtained along well chosen profiles across the Orion Bar. The JWST data allows us, for the first time, to spatially resolve the steep variation of dust emission at the illuminated edge of the Orion Bar PDR. By considering a dust model with carbonaceous nano-grains and submicronic coated silicate grains, we derive unprecedented constraints on the properties of across the Orion Bar. To explain the observed emission profiles with our simulations, we find that the nano-grains must be strongly depleted with an abundance (relative to the gas) 15 times less than in the diffuse ISM. The NIRSpec and MRS spectroscopic observations reveal variations in the hydrogenation of the carbon nano-grains. The lowest hydrogenation levels are found in the vicinity of the illuminating stars suggesting photo-processing while more hydrogenated nano-grains are found in the cold and dense molecular region, potentially indicative of larger grains.


Introduction
Interstellar dust is an essential component of interstellar matter (ISM) and plays a key role in the formation of stars and protoplanetary disks.Ubiquitous in the ISM with a dust-to-gas mass ratio of a percent, dust is formed of tiny solid grains with sizes between a few 0.1 nm and a few 100 nm.Analysis of ISM observations has long shown that dust grains control the transfer of radiation, heat the gas via the photoelectric effect and form efficiently important molecules like H 2 on their surface.Dust grains also carry a significant part of the ISM charge and follow evolution.This evolution may change the abundances of the dust populations but also the structure (porous, amorphous) and composition of the grains.These latter changes are however more difficult to derive from observations because they require detailed and spectrally extended dust properties for models of dust analogues.Conversely, abundance and size distribution changes of grains are directly reflected in the dust emission.In this context, former work suggested that the main impact of dust evolution is to change the relative abundance of the small (1 nm in size) to large (100 nm in size) grains (Arab et al. 2012;Schirmer et al. 2020).In addition, the observed large variation of carbon gas-phase abundance (see references in Compiègne et al. (2011)) suggests that this element is efficiently cycled in and out of grains by accretion and destruction processes (Jones 2014).The gas-phase abundance of carbon is thus enhanced by energetic processes (shocks or UV light) which break up dust into smaller carbonaceous grains.This active carbon cycle and extinction observations (Parvathi et al. 2012) suggest that the small dust population is in fact carbonaceous nano-grains.Proposed to explain the aromatic infrared bands (AIBs) observed between 3 and 20 µm, these nano-grains must be, at least partially, aromatic and hydrogenated.Plausible interstellar nano-grains range from compact, symmetric species such as Polycyclic Aromatic Hydrocarbons (PAHs) or fullerenes (Leger & Puget 1984;Allamandola et al. 1985;Sellgren et al. 2010) to more disordered ones such as amorphous hydrogenated carbon or a-C:H (Jones et al. 2013) .It must be noted here that models predict that nanograins (i) dominate the gas photoelectric heating rate (Bakes & Tielens 1994;Habart et al. 2001;Berné et al. 2022a), (ii) carry most of the surface for H 2 formation (Bron et al. 2014;Jones & Habart 2015), and (iii) are the carriers of the UV dust extinction (Weingartner & Draine 2001).Describing the physical state of the gas that will form stars thus requires knowledge of the properties and abundance variations of nano-grains.
Ever since the advent of UV and IR data in astronomy, dust models have been designed to explain the main observables (emission and extinction) without evolution.Recent data collected during the ISO, Akari, Spitzer, Planck, and Herschel missions have brought a first look at dust evolution however limited by the low angular resolution (e.g.Compiègne et al. 2008;Arab et al. 2012;Mori et al. 2012Mori et al. , 2014;;Pilleri et al. 2015;Schirmer et al. 2022).The advent of JWST data offers new opportunities to study and model dust evolution.In our Galaxy and in external galaxies, interstellar matter close to massive stars are privileged regions to target because of the prevailing high excitation conditions (strong UV radiation field and high gas density) which trigger dust evolution and also provide enough emission for detailed observations.The evolution of such regions and of more quiescent ones exposed to the standard interstellar radiation field (Mathis et al. 1983) is due to the radiative heating, this is why they are called photon-dominated regions (PDRs).Bright PDRs such as the Orion Bar are active regions where the star formation process is best studied (Hollenbach & Tielens 1999).Within the framework of an Early Release Science (ERS) observation program, the JWST space IR facility has observed gas and dust emissions towards the Orion Bar -the brightest PDR in the sky -at high angular resolution in imaging and spectroscopy (Berné et al. 2022b;Habart et al. 2023;Peeters et al. 2023), providing an ideal data set to study the evolution of dust and its nature.Using these data, we analyse the dust emission across the Orion Bar PDR and compare it to radiative transfer simulations in order to derive the effect of dust evolution.
The structure of the paper is as follows.In Sect.2, we provide an overview of previous studies of the Orion Bar.We recall the JWST spectroscopic and photometric observations in Sect.3. In Sect.4, we detail the THEMIS dust model and the radiative transfer code used to simulate dust emission.
We compare the modelled dust emission to observations and discuss the resulting constraints on dust properties in Sects.5 and 6.Finally, in Sect.7, we discuss our results and conclude in Sect.8.

The Bar PDR
The Orion Bar, a strongly ultraviolet (UV) irradiated PDR, is an escarpment of the Orion molecular cloud (OMC), the closest site of ongoing massive star formation.The gas density in the ambient molecular cloud is estimated to be n H = 0.5-1.0 × 10 5 cm −3 (Tielens & Hollenbach 1985;Hogerheijde et al. 1995).The Bar is illuminated by the O7-type star θ 1 Ori C, the brightest and most massive member of the Trapezium cluster at the heart of the Orion Nebula (e.g.O'Dell 2001).The Trapezium cluster creates a blister H ii, region that is gradually eroding into its parent cloud.A large cavity has been carved out of the molecular gas and the concave inner structure slopes to form the Orion Bar (Wen & O'dell 1995;O'Dell 2001).
The edge of the Orion Bar PDR corresponds to the ionization front (IF), where the ionized gas (H ii) recombines to neutral gas (H i).Right after this edge, photons have energies below 13.6 eV and the incident far-UV (FUV) radiation field is G 0 = 2-7 × 10 4 times the Habing unit (e.g.Peeters et al. 2023).This atomic H i layer extends over 10-20 ′′ (or 0.02-0.04pc) until dissociating FUV photons are sufficiently attenuated to let hydrogen become molecular (e.g.Habart et al. 2023) across a dissociation H/H 2 front (DF).At this position the visual extinction to the IF is A V ≃ 0.5-2 mag.Beyond the DF, between A V = 2 and 4 mag, the C + /C/CO transition takes place (Tauber et al. 1995) and the PDR becomes molecular.
Recently the ALMA and JWST images at high angular resolution (∼0.1-1 ′′ ) have revealed small-scale (∼0.004 pc) overdense structures, sharp edges, and bright emission from embedded proplyds (object 203-506;Champion et al. 2017), possibly induced by the intense UV-radiation field (Gorti & Hollenbach 2002;Tremblin et al. 2012).All these small-scale structures are superimposed over a rather homogeneous layer of atomic H i gas that extends almost 0.01 pc from the PDR edge (Habart et al. 2023).
Previous studies of dust evolution in PDRs relied on data with a limited angular resolution (e.g.Compiègne et al. 2008;Arab et al. 2012;Mori et al. 2012Mori et al. , 2014;;Pilleri et al. 2015;Schirmer et al. 2022) that lead to significant broadening and confusion of the emitting structures.In particular, in their preliminary study of dust evolution in the Orion Bar based on Spitzer and Herschel data (resolution of a few arcseconds), Schirmer et al. (2022) were only able to place upper and lower limits on the smallest size and abundance of nano-grains, respectively.We show below that a detailed modelling of high angular resolution JWST data (∼ 100 times higher than that of former observations) allows us to obtain quantitative constraints on the size distribution of nano-grains.In addition, original constraints on the optical properties of these nano-grains are derived from the JWST MRS and NIRSpec spectroscopic data.This dust characterisation is important for understanding the impact of radiative feedback on the dynamics and chemistry of the ISM (e.g.Schirmer et al. 2021).
Article number, page 2 of 18 M. Elyajouri et al.: Modelling the dust evolution across the illuminated egde of the Orion Bar

Description of the JWST data used in this study
The photometric and spectroscopic data used in our study were obtained from the ERS project PDRs4All (#1288) 'Radiative Feedback from Massive Stars as Traced by Multiband Imaging and Spectroscopic Mosaics'1 (Berné et al. 2022b).

Observations
The MIRI and NIRCam maps are presented in Habart et al. (2023), the MRS speçctra in Chown et al. (2023), and the NIRSpec spectra in Peeters et al. (2023).An in-depth description of the reduction procedure for these data, including corrections for instrument signatures, background removal, calibration, rectification, and final combination, can be found in these three papers.The fractional contributions of line, aromatic infrared bands (AIBs), and continuum emission to our NIRCam and MIRI images are examined in detail in the Science Enabling Product 4 article of Chown et al., in prep.

Imaging data
We use the F335M and F480M NIRCam filters.The F335M filter is dominated by the 3.3 µm AIB, with some contribution of the aliphatic bands at 3.4 µm and the continuum.The F300M filter traces the emission continuum on the shorter wavelength side of the F335M filter, but it is affected by scattered light.We also use the F480M filter which is dominated by the continuum, the contribution of H i lines, and H 2 line being negligible (Chown et al., in prep.).We also use MIRI data in the F770W and F1130W filters, which include AIBs, and the F1500W and F2550W filters, predominantly capturing continuum emission.
To showcase the data collected, Figs. 2 and 3 present images with our selection of NIRCam (F335M and F480M) and MIRI filters (F770W, F1130W, F1500W, and F2550W).These filters trace the 3.3, 3.4, 7.7, and 11.3 µm bands alongside the continuum emission.The F335M, F480M, F770W, and F1130W maps are dominated by the emission of very small carbonaceous grains, while larger grains contribute to the F1500W and F2550W filters.The contribution of the different grain sizes as a function of wavelength and radiation field strength is well illustrated in Figs.8-10 of Compiègne et al. (2011) andin Figs. 3 and5 of Schirmer et al. (2020).
The emission profiles perpendicular to the IF and the Bar are shown in Figs. 2 and 3, allowing us to probe dust emission variations across the PDR with increasing distance from the excitation source.These profiles provide insights into the density increase inside the PDR which occurs as the FUV radiation field decreases deeper into the Bar (see Habart et al. 2023, and Fig. 4 adapted from their study).

Spectroscopic data
As shown in their Fig. 1, Peeters et al. (2023) extracted NIR-Spec spectra from five strategically positioned apertures: in front of the ionization front (IF), at the peak of the 3.3 µm band, and at the H i/H 2 dissociation fronts (DF1, DF2, and DF3).These locations within the Orion Bar PDR are crucial as they represent the ionized, atomic, and warm molecular regions.Chown et al. (2023) extracted MIRI template spectra using the same aperture positions as Peeters et al. (2023).By combining these MIRI spectra with the NIRSpec templates, we effectively cap-tured the entirety of the AIB emission across the five regions of interest.In this paper, we specifically focus on three of the five combined NIRSpec and MRS template spectra representing key areas: the H II region, the atomic zone, and the DF3.

Methods and tools
In the following, we first explain the methodology employed for performing radiative transfer calculations.We then discuss the dust model we are using, which includes evolution from the diffuse interstellar medium to dense molecular clouds.Finally, we outline the assumptions made to represent the cloud geometry, the radiation field, and the free parameters, and an approach to constrain these parameters.

Radiative transfer
We use the radiative transfer program Scattering with OpenCL2 (SOC; Juvela 2019) which is based on the Monte-Carlo method, the simulation of large numbers of photon packages that represent the actual radiation field.This enables self-consistent modelling of the emission, radiation transport, and extinction throughout the model volume.This method relies on probability to determine the distance a photon travels before interacting with a grain, the type of interaction (scattering or absorption), and the direction of the photon in case of scattering.The calculations are performed with the coupling of the Monte Carlo radiative transfer code SOC and the dust emission and extinction code DustEM3 (Compiègne et al. 2011).Dust properties are defined by the absorption and scattering efficiencies Q abs and Q sca , and the asymmetry parameter of the Henyey-Greenstein scattering phase function g, which are determined using the dust model defined in the next section.

Dust model description
The Heterogeneous dust Evolution Model for Interstellar Solids4 (THEMIS; Jones et al. 2017) is built upon the foundations of the laboratory-measured properties of physically reasonable interstellar dust analogue materials.These include the family of hydrogenated amorphous carbon materials, from H-poor (a-C) to H-rich (a-C:H), collectively a-C(:H), and also amorphous olivine-type and pyroxene-type silicates with iron and iron sulphide nano-inclusions, a-Sil (Jones et al. 2013;Köhler et al. 2014).
In the initial model for the diffuse ISM, small grains up to 20 nm in size are entirely made up of aromatic-rich, hydrogenpoor amorphous carbon (a-C), with a power-law size distribution and an exponential cut-off and optical properties that are coherent and continuous through the entire size range.These a-C grains are responsible of both AIBs, generally attributed to the astrophysical PAHs, and infrared continuum.
Larger grains have a core/mantle (CM) structure (a-C:H/a-C CM grains) with a-C mantles (depth ∼ 20 nm) enveloping a-C:H cores, and a-Sil/a-C grains with a-C mantles (depth ∼ 5-10 nm) (Jones et al. 2013;Köhler et al. 2014).The two populations of larger grains have log-normal size distributions peaking between 100 and 200 nm.THEMIS takes into account the changes in dust properties through the effects of UV photo-processing, accretion, and coagulation, as the gas density evolves (Köhler ).In the transition from the diffuse ISM to dense clouds or in environments with attenuated radiation or efficient rehydrogenation, the carbonaceous mantles can retain their hydrogen content, leading to the formation of grains with two mantles (core/mantle/mantle grains, CMM).This formation can occur due to the coagulation of small aromatic-rich carbon grains onto larger grains and by the accretion of C and H atoms from the gas phase to form an a-C:H mantle.Further, within dense clouds, these CMM grains can coagulate into aggregates (AMM).Fi-nally, the formation of ice mantles (I) on the aggregates (AMMI) can occur in the densest regions, shielded from energetic photons and gas molecules can freeze out onto the grain surfaces, possibly allowing chemistry to proceed (Köhler et al. 2015).
Another key feature of THEMIS is the recognition of the inherent variability in the optical properties of hydrocarbon grains.As these or semi-conductor materials darken upon exposure to UV light and through thermal annealing, there is a decrease in the band gap or optical gap energy, E g (Iida et al. 1985;Smith Article number, page 4 of 18 1984).This property, and the corresponding effects on the optical properties, is crucial for understanding the evolutionary histories of hydrocarbon grains in the ISM (Duley 1996;Jones 2009Jones , 2012a)).The band gap E g correlates linearly with the fractional atomic hydrogen content X H (≃ E g /4.3; Tamor & Wu 1990), and the band gap inversely depends upon the number of aromatic rings in the aromatic domains (Robertson & O'Reilly 1987).
A key element of the THEMIS framework is a self-consistent treatment of the evolution of the dust material properties (size distribution, chemical composition and structure) as they react to and adjust to the local radiation field intensity and hardness, and to the gas density and dynamics.For the carbonaceous nano-grains, all these processes, interactions will lead to changing of the C/H in the grains and also to changing of the dust structure for more or less aliphatic-or aromatic-rich.As detailed in Jones (2012a,b,c), this is fully captured by the band gap in the optEC model.
In the forthcoming sections, we will test two variations of the THEMIS dust model.The dust populations in each case demonstrate specific size distributions and optical properties, indicative of the ISM's dynamic cosmic dust evolution.They are: -THEMIS for the diffuse ISM (DISM), featuring two dust populations and three components consisting of a-C, a-C(:H)/a-C, and a-Sil/a-C.
Article number, page 5 of 18 -The a-C + Mix(1:50) model, which substitutes THEMIS aggregates with spherical grains composed of 2/3 a-Sil and 1/3 a-C.These spherical grains exhibit a 50% porosity level and a log-normal size distribution with a minimum size a min,a-C of 0.5 nm, a peak size a 0 (free parameter in the following) and a maximum size a max of 4.9 µm.This case extends the size distribution to larger grains, using the method proposed by

Geometry and density profile
We take a plane-parallel geometry to represent an edge-on PDR, where all photons enter perpendicular to the PDR's edge.The cloud model is discretized onto a Cartesian grid consisting of N X × N Y × N Z cubes with each cube having a size of 0.001 pc (equivalent to 0.5 ′′ ).This cell size is three times smaller than the smallest scales observed in the profiles (see Fig. 2).We set N Z = 100, which corresponds to the number of cubes along the PDR-star axis, and N Y = 10 (to avoid edge effects), corresponding to the number of cubes along the axis Y perpendicular to both the PDR-star axis Z and the line-of-sight axis X.The number of cubes in the PDR along the line-of-sight, N X , depends on the value of the length of the Bar along the line-of-sight, l PDR , given by the equation: l PDR = N X × cell size.This is a free parameter which will be constrained.We take the density profile used in previous gas and dust models of PDR edges (Habart et al. 2005;Arab et al. 2012;Schirmer et al. 2020;Schirmer et al. 2022): where z is the distance from the edge of the PDR.The power-law index γ governs the steepness of the density front and is equal to 2.5.The maximum value of the density n 0 is reached at the depth z 0 and is assumed to remain constant for z > z 0 .
Article number, page 6 of 18

Radiation field
The PDR is illuminated by the O7-type star θ 1 Ori C (Sota et al. 2011), the most massive member of the Trapezium young stellar cluster, and its spectrum can be modelled by a T = 38 000 K blackbody.We fix the incident radiation entering the PDR so that, at the edge of the PDR, G 0 = 2.6 × 10 4 in Habing units (in agreement with estimates from UV-pumped IR-fluorescent lines of O i by Marconi et al. 1998 andPeeters et al. 2023 which indicate G 0 = 2.2-7.1 × 10 4 ).The star was assumed to be at a distance of 0.24 pc from the PDR front, which corresponds to the projected distance between the PDR and the star (about 2 ′ northwest of the Bar, e.g.Habart et al. 2023).This radiation field is calculated on a grid of 757 frequencies regularly sampled, except for frequencies close to bands where we use higher resolution in order to sufficiently resolve these features (see Fig. 6 for a schematic illustration of the PDR).

PSF convolution and band integration of our dust emission model
In order to compare our dust emission models with observational data, it is crucial to perform a post-processing treatment of the model outputs, including integration over specific photometric bands and convolution with the Point Spread Functions (PSFs) of the JWST.The NIRCam and MIRI PSFs are taken from WebbPSF5 (Perrin et al. 2014).Following the convolution, we integrate our model outputs over the different photometric bands using the photon-to-electron conversion efficiency.

Methodology and free parameters
Our methodology involves a grid-based analysis, where we explore a range of free parameters to construct a set of models using THEMIS and SOC.These models allow us to examine the influence of various parameters on the resulting dust emission profiles and compare them with photometric data from NIRCam and MIRI, as well as with spectra from NIRSpec and MRS.
The free parameters in our model6 incorporate elements associated with both nano-grain and larger grain properties, as well as the density profile: 1. Band gap energy (E g ) -hydrogenation state: is completely determined by the electronic structure of a-C(:H) nano-grains and is primarily and strongly constrained by the spectral data.This parameter particularly affects the 3.3-to-3.4µm band ratio and as such is the initial focus of our grid exploration using NIRSpec spectra.2. Abundance (M a-C /M H ): The a-C dust mass to gas ratio, representing the relative abundance of the smaller dust grains (∼ sub 20 nm radius), which determines the global level of emission from the NIR to the MIR. 3. Minimum size (a min,a-C ): Denotes the smallest a-C nanograin size in the distribution.This parameter impacts both the spectral shape and emission continuum slope.4. Slope (α): of the a-C power-law size distribution.Like a min,a-C , this plays a significant role in determining the spectral shape and emission continuum slope.5. Size a 0,Mix(1:50) : Peak of the log-normal size distribution for the large grains (Mix 1:50).This mostly impacts the dust temperature and therefore where the SED peaks in the MIR to FIR wavelength range.6. Distance (z 0 ) and density (n 0 ) are the density profile parameters already defined in Eq.(1).They influence the spatial emission profiles (steepness and width) for all of the NIR-Cam and MIRI filter fitting.7. Length of the PDR (l PDR ) along the line of sight does not affect the shape of the dust spectrum and is considered as a multiplying factor on the dust spectrum.The intensity increases linearly with l PDR without altering the shape of the dust spectrum (see Schirmer et al. 2020).
The influence of variations in the three parameters M a-C /M H , a min,a-C , and α on both the dust size distribution and the associated spectra in the optically thin limit are shown in Schirmer et al. (2020).Our challenge here is to derive realistic self-consistent models compatible with the observations in an eight-dimensional space defined by E g , M a-C /M H , a min,a-C , α, a 0,Mix(1:50) , z 0 , n 0 , and l PDR .To tackle this problem, we adopt a strategy that involves fixing or constraining certain parameters whenever possible.Among these, E g is particularly well constrained by NIR spectral data and thus forms the starting point of our grid exploration using NIRSpec spectra (see Sect. 5).By evaluating the spectral influence of various E g values, we are able to refine our exploration domain by narrowing down the parameter space.A similar approach is adopted for the parameters a min,a-C and α, due to their significant impact on the spectra (band shape and continuum slope).With E g , a min,a-C and α established, we proceed with a more targeted search, varying z 0 and a 0 within reasonable ranges (see Sect. 6.1).This approach generates a grid of models, allowing us to reduce the range of parameter values and complexity of our exploration grid space (Sect.6.2).The resulting emission profiles are then compared with the observed photometric data, aiming to identify the parameter combinations compatible with the observed profiles.

Adjustments of NIRspec template spectra
We begin by conducting simulations using the DustEM code without radiative transfer, in order to produce a grid of THEMIS Article number, page 8 of 18 SEDs for G 0 = 2.6 × 10 4 .The grid covers variations in the minimum grain size (a min,a-C ) from 0.35 to 0.8 nm, in increments of 0.05 nm, the power-law slope (α) from −12 to −5 with a step size of 0.25, and the hydrogenation state (E g ) from 0.01 to 0.1 eV.The modelled SED grid is then compared to three NIRSpec spectral templates corresponding to different regions in front of and inside the Bar (Peeters et al. 2023).These three spectra, shown in Fig. 5, correspond to (1) the ionized region in front of the Bar (top), (2) the atomic region of the Bar (middle), and (3) the warm molecular region (bottom), respectively.
In comparing the spectra, we note variations in the 3.3 µm band (associated with aromatics) and the 3.4 µm band (associated with aliphatics), as well as differences in the band-tocontinuum ratio.We initially compare our model with the DISM THEMIS model (with a minimum grain size of 0.4 nm and a 0.1 eV band gap, corresponding to a hydrogenation of X H ∼ 0.02).Although the continuum fits relatively well for some regions, we find discrepancies in the 3.4 µm band and the overall continuum level for other regions.In the H ii region (top panel of Fig. 5), the continuum is underestimated by a factor of ∼2, while on the atomic Bar (middle panel of Fig. 5), the continuum fits relatively well, but the 3.4 µm band is significantly overestimated (factor of ∼2).Regarding the molecular filament DF3 (bottom panel of Fig. 5), the fit is reasonably good for the continuum and the bands 7 .
To address the discrepancies found for the H ii and atomic regions, we vary the minimum grain size and the degree of hydrogenation through E g , in order to adjust the spectra.Figure 5 illustrates the effect of each parameter on the modelled SED.As expected, increasing the minimum size of the smallest and therefore the hottest grains decreases the SED intensity at the shortest wavelengths, having an effect on both the bands and the continuum.The effect of decreasing the band gap is most visible in the 3.3/3.4ratio as it corresponds to removing hydrogen atoms from the grain structure and thus to a decrease in the 3.4 µm aliphatic band (see also Jones 2012c; Jones et al. 2013 for details).The 3.4 µm intensity is modelled by the integrated intensity and the best agreement is achieved with the following parameters: -For the H ii region, E g = 0.03 eV and a min = 0.55 nm.
-For the atomic region, E g = 0.03 eV and a min = 0.475 nm.
-For the warm molecular region (DF3), E g = 0.07 eV and a min = 0.425 nm.
In order to adjust the ratios of the 3.3 and 3.4 µm bands and match the observed continuum level, it was necessary to decrease the band gap and increase the minimum size a min in the transition from denser shielded regions to more irradiated regions.This indicates that as they get closer to the star, the hydrogen content in the grain structure decreases, consistent with previous photo-processing studies (see for instance Jones 2014, and references therein).Additionally, we found that in the irradiated regions, slightly larger grains are needed because the smaller ones are apparently more easily destroyed.
It is important to highlight that in the H ii region, the contribution of free-free emission is likely to be significant.This could account for the observed underestimation of the continuum level attributed to dust emission, as suggested by Haraguchi At the H ii region, the observed emission encompasses both the gas and dust from the ionized region as well as the PDR face-on from OMC-1.However, our results indicate an evolution in both hydrogenation and grain size from DF3 to the H ii region, suggesting that emission from the ionized medium is a contributing factor.This observation underlines the complexity of the H ii region's emission profile, where multiple sources and processes are at play.Now and hereafter in Sect.6, we will focus on the PDR edge, that is the atomic region.To best determine the range of possible E g values in this region model results for different band gap energies (E g = 0.01, 0.03, 0.05, and 0.07 eV) are compared to the observed data in the middle panel of Fig. 5.This comparison highlights that the E g = 0.03 eV model best reproduces the observations in the atomic PDR region, closely following the trend of the data points and accurately reproducing the ratio between the 3.3 and 3.4 µm bands.In contrast, the models with E g ≥ 0.05 eV overestimate the intensity in the 3.4 µm band, while the model with E g = 0.01 eV underestimates the intensity in the same band.Thus, we assume a fixed value of E g = 0.03 eV for the nano-grain band gap.Moreover, based on the comparison between our models and the NIRSpec spectra in the atomic region, we found that the parameter space for a min,a-C and α to be explored to fit the NIRCam and MIRI data is the following: 0.4 ⩽ a min,a-C ⩽ 0.6 nm, −6 ⩽ α ⩽ −4.

Radiative transfer modelling in the atomic region
As part of our comprehensive approach, we have constructed a model grid for NIRCam and MIRI photometric data in order to model the brightness profiles observed at the edge of the PDR, in particular, the rise in emission just after the IF and the decay due to extinction8 .This zone is important because it is where the density of the atomic gas increases sharply and the extinction (in the UV illuminating direction) perpendicular to the Bar begins.In this atomic zone, the dust significantly contributes 9 (at least 50% for a density n 0 ≤ 1 × 10 5 ) to the gas heating via the photoelectric effect and it is crucial to provide constraints on its properties.The PDR edge was unresolved in the IR until now.Thanks to the JWST, we are now able to spatially resolve this zone, and we can put constraints on the dust properties.
part due to material along the line of sight located in front of the PDR, in the foreground face-on surface layer (as seen in geometry on Fig. 5 of Habart et al. 2023).This emission may originate in the compressed region that is still illuminated directly by the ionizing stars. 9Based on the PDR Meudon code, the photoelectric effect along with the heating by H 2 cascades (vibrational de-excitation) will dominate the heating in the pure atomic region of the Orion Bar.
6.1.Constraining z 0 and a 0,Mix(1:50) One of the key parameters we aim to constrain in our analysis is the depth threshold (z 0 ) where the density profile reaches its maximum n 0 .The determination of this parameter is largely independent of the dust properties, and is primarily reliant on the width of the dust-modelled emission.A value of z 0 = 0.003 ± 0.0005 pc is necessary to fit the brightness profiles.We note that z 0 is much smaller than the value of 0.025 pc deduced from the modelling of Spitzer and Herschel observations due to their lower spatial resolution (Schirmer et al. 2022).
Another important parameter is the size of the pseudoaggregates, a 0,Mix(1:50) .The variation of a 0,Mix(1:50) of the lognormal size distribution of the large grains has a significant effect on the dust SED at all wavelengths.This is because the dust size distribution affects the opacity from the visible to the sub-Article number, page 10 of 18 millimetre, which in turn affects the amount of radiation that is absorbed and re-emitted by dust. Figure 7 shows the SED variations with a 0,Mix(1:50) .Increasing a 0,Mix(1:50) broadens the lognormal grain size distribution.This results in a higher abundance of large grains, which are more efficient at absorbing longwavelength radiation (IR) and a smaller population of small dust grains, which are more efficient at absorbing short-wavelength radiation (UV).As a result, increasing the value of a 0,Mix(1:50) will tend to shift the peak of the SED to longer wavelengths, as the dust will be cooler.Regarding our imaging data, this parameter has a significant influence on the ratio between the emission in the 15 and 25.5 µm bands.Through an exploration of our preliminary grid, we concluded that a mean grain size (a 0,Mix(1:50) ) of 100 nm most accurately produces the brightness profiles observed at 15 and 25.5 µm.

Remaining free parameters
With E g , a 0,Mix(1:50) , z 0 now established, the remaining parameters to determine are a min,a-C , α, M a-C /M H , n 0 , and l PDR .To explore these five parameters, we examined a more extensive model grid, spanning the following values: 1. M a-C /M H varies from 0.0025 × 10 −2 to 0.17 × 10 −2 on a 10point linear grid.2. a min,a-C varies from 0.4 to 0.6 nm on a 10-point linear grid.3. α varies from −4.0 to −6.0 in steps of 0.2. 4. n 0 varies from 4 × 10 4 to 1.5 × 10 5 cm −3 in steps of 10 4 .5. l PDR varies from 0.08 to 0.14 pc in steps of 0.01 pc.
The maximum value for M a-C /M H corresponds to the value found in DISM.The range of a min,a-C and α are discussed in Sect. 5. Concerning n 0 and l PDR , the range is taken to include the estimates from Habart et al. (2023) and Peeters et al. (2023).In any case, we have checked that for the five parameters, all values outside these ranges do not allow us to reproduce the data.
We are actually using an early stage of data reduction, which may have some calibration errors.In addition, the 3D model of the illuminated edge of the Bar is extremely simple, with a plane-parallel geometry observed edge-on, and illuminated by a star in the plane of the sky.Therefore a χ 2 analysis to derive the 'best' model is premature.Instead, we opted for a grid search across the five parameters to be constrained (M a-C /M H , a min,a-C , α, n 0 , l PDR ) in order to derive a set of parameters physically coherent and which provides a satisfactory agreement between a self-consistent model and the data (brightness profiles and spectra).
We conducted multiple grid tests, varying the five parameters.Initially focusing on the DISM population, it was not possible to adjust the ratio between the 15 µm and 25.5 µm bands.Consequently, we explored a second dust population: the a-C + Mix(1:50) model (see Sect. 4.2).The model best parameters were selected in order to minimize the quadratic difference between observed and modelled peak intensities across all six bands.Numerous models were tested (combinations of the 5 free parameters), and the only model that simultaneously reproduced all bands within the uncertainty range is the one presented in Sect.6.3.Other models, like the one in Fig A .1, either overestimated or underestimated at least one band.Since our analysis encompassed all parameters of the dust model, we therefore believe that the simple geometry of the model is its primary limitation.

Comparison with NIRCam and MIRI brightness profiles
We took into account an offset for the observed emission which is crucial for the adjustment process.Specifically, the value of the offset is taken in the H ii region, just in front of the illuminated edge for each profile (shown as a dashed red line, top panel of Fig. 8).This offset may have several physical origins (OMC-1 seen face-on or H ii material, see Sec. 5) and its modelling is beyond the scope of this paper.Additionally, we incorporated 20% error bars to account for the calibration accuracy.
The model from our grid which provides the best agreement with the six brightness profiles is presented in Fig. 8.The a-C abundance M a-C /M H (0.0109 × 10 −2 ) is found to be 15 times less than in the diffuse ISM (0.17 × 10 −2 ).The minimum grain size a min,a-C = 0.4 nm and the slope α = −5 are consistent with the values found in the diffuse ISM (Jones et al. 2013).The gas density n 0 is estimated to be about 9 × 10 4 H cm −3 at the nanograin emission peak.This is in agreement with estimates made by Habart et al. (2023) and previous estimates from atomic gas FIR lines (e.g.[C ii] and [O i] fine-structure lines, Bernard-Salas et al. 2012) and estimates from the Raman-scattered wings of hydrogen Hα lines in the Bar (Henney 2021).Interestingly, this density is also close to the one derived at the DFs from NIRSpec data (Peeters et al. 2023).The value of l PDR = 0.12 pc we find is consistent with the value reported by Peeters et al. (2023).
Notably, for the first time, we are able to simultaneously fit the profiles in all photometric bands (3.35 and 4.8 µm NIRCam filters and 7.7, 11.3, 15, and 25.5 µm MIRI filters) within the calibration accuracy.This achievement is significant because the bands at 3. 35, 4.8, 7.7, and 11.3 µm cover both aromatic and aliphatic features and also the continuum, whereas the 15 and 25.5 µm bands only cover the dust continuum.The successful reproduction of the band-to-continuum ratio and band shapes in the Orion Bar suggests that the THEMIS model also works well in highly irradiated regions, despite the fact that it was constructed to explain DISM dust.

Comparison with NIRSpec and MRS spectra
Using the derived parameters, the bottom panel of Fig. 8 shows the computed spectrum at the distance of 0.003 pc from the PDR edge which corresponds to the peak emission in all photometric bands.In Fig. 8, we overplot the NIRSpec and MIRI spectra obtained within the atomic region of the Bar.Although the modelled and observed spectra are not obtained at exactly the same position on the sky (there is a slight offset of a few arcseconds, see Fig. 2), both trace material in the atomic region of the Bar, allowing for a direct comparison.
The comparison between the modelled and observed spectra shows a good agreement both in the NIR and mid-IR part of the SED, covering the range from 3 to 25 µm.This global dust modelling approach is promising and highlights the effectiveness of our approach.However, we also notice some differences between the model predictions and the observed data: 1.The lack of emission features around 5-6 µm.Jones (2012a,b,c) assigned bands to the various C-C and C-H bonds according to the data available at that time and only bands with well-determined origins can be included within THEMIS.This means that an update of the positions and strengths of the various bands will need to be made in the light of the JWST observations once more spectra of regions with various G 0 /n H will become available.This update will have to be based on these new observational constraints and, particularly, on new laboratory data and Article number, page 11 of 18 band assignments that build upon the available laboratory data (Dartois et al. 2004;Pino et al. 2008;Carpentier et al. 2012).
2. The model tends to overestimate the continuum around 7 µm, between the 6.2 and 7.7 µm features.This is again related to the need to refine the band assignments based upon future laboratory data.
3. Although we can accurately reproduce the features at 3.3-3.4µm in our model, we encounter challenges in adjusting the continuum around 3.5 µm.To address this, we consider the possibility of slightly increasing the value of a min,a-C , which leads to an improved fit between 3 and 4 µm as shown in Fig. C.1.However, such an adjustment results in a slightly less accurate fitting of the F480M profile.
4. The observations present a plateau above our model between 9 and 11 µm.This could be related to the presence of nano-silicates (e.g.Hensley & Draine 2021;Chastenet et al. 2021;Zeegers et al. 2023) which are not included in THEMIS.For such small grains, stochastic heating would be efficient (Guiu & Bromley 2022), possibly leading to observable mid-IR emission in the PDR.Cesarsky et al. (2000) reported the observation of amorphous silicates in Orion.More precisely, they detected their broad emission band, centred at 9.7 µm, in the whole H ii region and around the star θ 2 Ori A. Knight et al. (2022) observed silicate emission in the H ii region only and found that this emission plummeted towards the IF of the Orion Bar.They proposed that the silicate grains were large and, while irradiation by Lyman-α photons in the H ii region would heat them sufficiently to cause significant emission at 9.7 µm, the PDR environment behind the IF would prevent such heating and emission.Our best model does not include isolated silicates but only silicates mixed with the population of carbonaceous grains in the form of large porous grains (the Mix(1:50) type grains presented in Sect.4.2).This has the effect of strongly attenuating the mid-IR vibrational bands of the silicates (e.g.Ysard et al. 2019).Then, above ∼13 µm, the shape of the model continuum does not match the observations either.This may be due to calibration issues (Channel 4 of MIRI/MRS) or differences in the equilibrium temperature of the large grains.This could also again be due to our simple modelling of aggregate grains based on spherical particles where materials are mixed according to effective medium theory (see the description of Mix(1:50) type grains in Ysard et al. 2019, and Sect. 4.2 for details).However, calculating new aggregate optical properties would require a huge computational effort to explore the expected parameter space and is beyond the scope of this paper.
In conclusion, it is too early to draw definitive conclusions regarding the observed differences between the model results and observations10 .

Nano-grain hydrogenation variation
Our study provides key insights into the hydrogenation variation within the different regions of the Orion Bar.The analysis of the NIRSpec spectra shows variations in nano-grain properties from dense to irradiated regions.In comparison with the predicted THEMIS spectra for the DISM, we see that the continuum is underestimated in the H ii region in front of the Bar, while the 3.4 µm band is significantly overestimated in the atomic region.However, the fit is reasonably good in the less excited molecular filament.
To address these discrepancies, we adjusted the grain properties in order to bring the continuum level within an acceptable range and adjust the band ratios, the grain band gap or hydrogenation level needs to be decreased.This suggests that the closer the region is to the star, the lower the hydrogen content in the grain structure 11 .Additionally, slightly larger grains (i.e. with higher minimum size) are found to be necessary in the irradiated regions, which is consistent with smaller grains being less resilient to harsh UV photons (see for instance Schirmer et al. 2022, and references therein).
The sensitivity of grain size and hydrogenation to the environment (Jones 2012a,b,c) highlights the importance of considering these factors when analysing dust and gas observations in irradiated media.As shown for instance by Schirmer et al. (2021), the exact properties of the smallest hydrocarbons indeed impact the chemistry and dynamics in PDRs and thus our understanding of the location of the various ionization fronts 12 .The gas heating through photoelectric effect, mostly efficient for the smallest grains, as well as the H 2 molecule formation rate are also strongly dependent on a min,a-C , α, M a-C /M H , and X H .By allowing us to compare the variations in the nano-grain properties for various G 0 /n H , the JWST observations will help us to improve our understanding of the stellar feedback influence on star formation, one of the main goals of the ERS program 1288.

Comparisons with previous works
Schirmer et al. ( 2022) have proposed an explanation for observations of dust from mid-to far-IR wavelengths in the Orion Bar.Their study was based on photometric data from Spitzer IRAC and MIPS in the mid-IR, and from Herschel in the far-IR.Using the same methodology (DustEM, SOC, and THEMIS grain model), our results are comparable albeit with some distinctions.In their study, the degeneracy between a min,a-C and α did not allow them to draw a conclusion on the nano-grain minimum size.While Schirmer et al. (2022) could only set an upper limit of a min,a-C ⩽ 0.8 nm, our results suggest a smaller size of a min,a-C = 0.4 nm.Additionally, our model indicates a less steep grain size distribution (α = −5) compared to theirs (α ⩽ −5.5).They also found significant depletion of nano-grains (10 −5 ⩽ M a-C /M H ⩽ 4 × 10 −5 ) compared to the gas and to their level in the diffuse ISM (0.17 × 10 −2 ).However, the abundance of nano-grains they found is four times lower than in our results M a-C /M H = 1.09 × 10 −4 .
We believe that the major difference explaining the discrepancies between our two models is the use of NIRSpec spectroscopic data.The spectroscopic data lead us to consider hydrogenation in the Orion Bar more than three times lower than in the diffuse medium, with the main constraint being the 3.3 to 3.4 µm band ratio.Insofar as Schirmer et al. (2022) kept the value used for the diffuse medium, E g = 0.1 eV (X H ∼ 0.02), their model greatly overestimated the band at 3.4 µm (see middle panel in Fig. 5).As the IRAC channel at 3.6 µm includes the two AIBs at 3.3 and 3.4 µm, the emission in this channel without any change in the size distribution is therefore also overestimated.In view of the free parameters of the Schirmer et al. (2022) study, the only way to decrease the total flux in the IRAC channel at 3.6 µm is to increase a min,a-C , the variation in the α parameter making it pos-sible to keep a certain level of continuum in the IRAC channel at 4.5 µm.
The photo-dissociation of CH bonds by UV photons is known to be efficient in the ISM both from lab experiments and theory (e.g.Welch & Judge 1972;Gruzdkov et al. 1994;Alata et al. 2014;Jones 2014).Variations in the 3.4 to 3.3 µm band ratios observed toward NGC 7023 also suggested variations in the hydrogen content of subnanometric carbonaceous particles (or PAHs Pilleri et al. 2015).According to Jones (2014), the photoprocessing timescale of small a-C:H particles would be only a few thousand years for the radiation field illuminating the Orion Bar.The use of JWST spectroscopic data allowed us to clearly demonstrate the efficiency of the hydrogen loss in such a PDR.It shows the value of the JWST and the PDRs4All program, in particular, for understanding stellar feedback in the interstellar medium.The scenario of equilibrium between the formation and destruction of the hydrocarbons responsible for emission in the AIBs will therefore have to be analysed in the light of these new results.However, this is beyond the scope of this paper as this will require a range of G 0 /n H to assess the efficiencies of all the physical and chemical processes involved.The necessary observations are currently being made for the GTO on the Horsehead (G 0 =100) and NGC 7023 (G 0 =1500) PDRs.

Consequences of a 15-fold decrease in abundance of nano-grains on the PDR
The nano-grain depletion found in this work is expected to have a profound influence on the physics and chemistry of the region, as discussed also in Schirmer et al. (2021) in the case of the Horsehead nebula.Indeed in a PDR, the nano-grains directly control several key processes namely, the photoelectric gas heating, the UV extinction (see Fig. It is important to emphasize that this extinction curve is attributed to dust primarily situated in the atomic region of the Orion Bar, with a visual extinction in the UV illuminating direction (perpendicular to the Bar) A V ≲ 0.7 mag.This distinction is crucial as it underlines that the observed extinction is not a result of dust within the denser regions of the bar.

Conclusion
We have investigated the dust emission by comparing model predictions with JWST/NIRCam and MIRI imaging, NIRSpec and MRS spectroscopic observations.We focussed on the illuminated edge of the Orion Bar, a highly irradiated dense PDR.
We used radiative transfer modelling to investigate the spatial distribution of the dust emission along a cut through the Orion Bar PDR and compute the emerging spectrum.We firstly examined the spectra in the ionized, atomic, and warm molecular regions and then the spatial variations in all used filters at the PDR edge.This scheme is reliable, efficient, and allows us to probe important variations of the dust properties (hydrogenation and abundance).The model incorporated essential parameters such as the density profile, incident radiation dilution factor, and the dust population consisting of THEMIS nano-grains transiently heated and larger pseudo-aggregates in thermal equilibrium with the radiation field.Our main results can be summarized as follows: In the atomic region, the grains are characterized by a band gap energy of 0.03 eV, which is significantly smaller compared to the typical 0.1 eV in the diffuse ISM.This lower band gap energy indicates reduced hydrogenation levels in the grains, consistent with the NIRSpec observations of the Orion Bar.To our knowledge, this is the first time that hydrogenation variations can be traced through a PDR.For the molecular region, it was necessary to employ larger grains than in the diffuse ISM.We modelled profiles corresponding to the 3. 3,4.8 NIRCam,and 7.7,11.3,15,and 25 µm MIRI filters, subsequently checking the consistency with the MRS spectrum.
Importantly, the investigation reveals that the nano-grains in the PDR edge are approximately 15 times less abundant than in the diffuse medium.This significant variation challenges the majority of existing PDR gas models that neglect such changes in nano-grain abundance.It highlights the importance of considering these variations in order to accurately estimate gas properties, including temperature, H 2 formation efficiency, and line intensities, which in turn will impact our understanding of the ongoing chemistry and dynamics in such a region.
Due to the unprecedented high angular resolution of the JWST, the spatial distribution of the dust emission reveal a very sharp illuminated edge (on scales of ∼1 ′′ or 0.002 pc) with a strong density gradient in the atomic region, just behind the IF.This strong density rise can be due to the sharp decrease in gas temperature at the IF, especially if the thermal pressures in the ionized and neutral regions are comparable.

Fig. 1 .
Fig. 1.Inner region of the Orion Nebula as seen by the JWST's NIRCam instrument with north up and east left.This RGB composite image adapted from Habart et al. (2023) is produced by combining three images in different filters: F335M (red), F470N (green), and F187N (blue) that trace emission from hydrocarbons (AIBs), dust and molecular gas (H 2 0-0 S(9) line), and ionized gas (Pa α line), respectively.The red box delineates regions of particular interest zoomed in and shown in Figs. 2 and 3.The perpendicular (NIRSpec) cut is shown in black (blue) dashed line.The white boxes indicate the apertures used to extract the three template spectra from Peeters et al. (2023).

Fig. 2 .
Fig. 2. Spatial distribution of the NIRCam (F335M and F480M) and MIRI (F770W) filters tracing the AIBs at 3.35, 4.8, and 7.7 µm and continuum dust emission.Left: Maps in the NIRCam (F335M and F480M) and MIRI (F770W) filters on a rotated Bar where the ionizing radiation is incident from the left.The blue inclined line shows the cut obtained from the NIRSpec field and the black line indicates the cut perpendicular to the Bar.A black circle shows the position of the IF for the perpendicular cut.Right: Surface brightness profiles in the NIRCam and MIRI filters as a function of distance from the IF for perpendicular and NIRSpec cuts, with the average position of the IF marked by black vertical dash-dotted lines.Units are in MJy sr −1 .The observed very high pixel values are attributed to the NIRSpec cut crossing two proplyds, which affects the distribution of pixel intensities in this region.
Ysard et al. (2019) as shown in the left panel of Fig. A.1.

Fig. 4 .
Fig. 4. Adapted from Habart et al. (2023).Schematic diagram illustrating the different components observed in the AIB surface brightness profile perpendicular to the Orion Bar.These components are shown for filter F335M (same as in Fig.2).In this study, we focus on modelling the first steep increase at the IF position, followed by a slower decrease indicated by the red line.The positions of DF1, DF2, and DF3 correspond to peaks in H 2 emission.

Fig. 5 .
Fig. 5. Comparison of observed surface brightness and DustEM models normalized to the 3.3 µm observed band (without radiative transfer) in the H ii (upper), atomic (middle), and DF3 (lower) regions.The black data points represent observations from NIRSpec, while the coloured lines show model predictions with varying band gap energies.In each figure, the DISM model is shown in blue dashed line and the best models for the H ii, atomic, and DF3 regions are shown in magenta, red, orange respectively.Article number, page 7 of 18

Fig. 7 .
Fig. 7. Dust spectra for different values of a 0,Mix(1:50) from 7.90 to 177.90 nm.The spectra are computed with radiative transfer using the SOC program and convoluted at the resolution of NIRSpec and MRS.

Fig. 8 .
Fig. 8.Comparison of the observed (JWST) and modelled dust emission in the atomic PDR region of the Orion Bar.Top: Comparison of the observed and modelled dust emission profiles in the atomic PDR region of the Orion Bar, using the best set of dust parameters in six photometric bands (3.3 and 4.8 µm; NIRCam filters, and 7.7, 11.3, 15, and 25.5 µm; MIRI filters).The observed dust emission is shown in black and the model fit in red, with uncertainties of 20% displayed in grey.The observed emission that we model (indicated by the first steep increase at the IF position, followed by a slower decrease) is shown in solid black line, while the other components of the observed emission are represented by dashed lines as explained in Fig. 4. The cut is perpendicular across the Orion bar, as shown in Fig. 1.Middle: Wavelength coverage and photon-to-electron conversion efficiency (PCE) of the NIRCam and MIRI imager filters used in our study.Throughput refers to photon-to-electron conversion efficiency.The central wavelengths of the filters are indicated with dashed vertical lines.The name of each filter is indicated in the top figure.The values of central wavelengths and bandwidths are presented in Habart et al. (2023).Bottom: Comparison of the dust emission spectral energy distribution (SED) calculated using our radiative transfer model at the AIB peak emission at a distance of 0.003 pc (represented by the red line), with template spectra for the atomic PDR in the Orion Bar from Peeters et al. (2023) (shown in dark blue for NIRSpec and black for MRS) in the wavelength range of 2.75 to 28 µm.
B.1), and the grain surface where H 2 forms.As a result the temperature profile and H 2 emission from a PDR edge are expected to be affected when nano-grains are depleted.Preliminary simulations with an updated version of the Meudon PDR code where the dust model of this paper has been fully integrated and coupled to the gas, suggests that the depletion of nano-grains is reflected in the intensity ratio of rovibrational to pure rotational H 2 emission line(Meshaka et al., in  prep.).The depletion in nano-grain abundance significantly influences the UV extinction curve, as clearly shown in Fig. B.1.