Draft version April 5, 2019 Typeset using LATEX twocolumn style in AASTeX62

The Polycyclic Aromatic Hydrocarbon Mass Fraction on a 10 pc scale in the Magellanic Clouds

Jer´ emy´ Chastenet,1 Karin Sandstrom,1 I-Da Chiang (江宜達), 1 Adam K. Leroy,2 Dyas Utomo,2 Caroline Bot,3 Karl D. Gordon,4, 5 Bruce T. Draine,6 Yasuo Fukui,7, 8 Toshikazu Onishi,9 and Kisetsu Tsuge8

1Center for Astrophysics and Space Sciences, Department of Physics, University of California, San Diego 9500 Gilman Drive, La Jolla, CA 92093, USA 2Department of Astronomy, The Ohio State University, 4055 McPherson Laboratory, 140 West 18th Ave, Columbus, OH 43210, USA 3Observatoire astronomique de Strasbourg, Universit´ede Strasbourg, CNRS, UMR 7550, 11 rue de l'Universit´e,F-67000 Strasbourg, France 4Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD, 21218, USA 5Sterrenkundig Observatorium, Universiteit Gent, Gent, Belgium 6Princeton University Observatory, Peyton Hall, Princeton, NJ 08544-1001, USA 7Institute for Advanced Research, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8601, Japan 8Department of Physics, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8601, Japan 9Department of Physical Science, Graduate School of Science, Osaka Prefecture University, 1-1 Gakuen-cho, Naka-ku, Sakai, Osaka 599-8531, Japan

ABSTRACT We present maps of the dust properties in the Small and Large Magellanic Clouds (SMC, LMC) from fitting Spitzer and Herschel observations with the Draine & Li(2007) dust model. We derive the abundance of the small carbonaceous grain (or polycyclic aromatic hydrocarbon; PAH) component. The global PAH fraction (qPAH, the fraction of the dust mass in the form of PAHs) is smaller in the +0.3 +1.4 SMC (1.0−0.3 %) than in the LMC (3.3−1.3 %). We measure the PAH fraction in different gas phases (H II regions, ionized gas outside of H II regions, molecular gas, and diffuse neutral gas). H II regions appear as distinctive holes in the spatial distribution of the PAH fraction. In both , the PAH fraction in the diffuse neutral medium is higher than in the ionized gas, but similar to the molecular gas. Even at equal radiation field intensity, the PAH fraction is lower in the ionized gas than in the diffuse neutral gas. We investigate the PAH life-cycle as a function of metallicity between the two galaxies. The PAH fraction in the diffuse neutral medium of the LMC is similar to that of the Milky Way (∼ 4.6 %), while it is significantly lower in the SMC. Plausible explanations for the higher PAH fraction in the diffuse neutral medium of the LMC compared to the SMC include: a more effective PAH production by fragmentation of large grains at higher metallicity, and/or the growth of PAHs in molecular gas.

Keywords: ISM: abundances – (ISM:) dust – (galaxies:) Magellanic Clouds

1. INTRODUCTION through photoelectric heating, which is the main mech- arXiv:1904.02705v1 [astro-ph.GA] 4 Apr 2019 Dust grains have a major impact on the energy bal- anism that heats the neutral gas (Wolfire et al. 1995). ance and chemistry of the interstellar medium (ISM), The effectiveness with which dust grains play these key and therefore are critical to the evolution of a . roles in the ISM depends on their intrinsic properties, e.g. size, charge, and chemical composition. To under- Acting as a favored surface for H2 formation, they are a key agent in the chemical balance of the ISM (Le Page stand the effects of dust grains on the ISM, we need to et al. 2009; Le Bourlot et al. 2012; Bron et al. 2014). understand their properties. Dust grains are also an efficient heat source for the ISM In this regard, the smallest grains are of particular interest. The small carbonaceous grain component is thought to be in the form of polycylic aromatic hydro- Corresponding author: J´er´emy Chastenet carbons (PAHs; Leger & Puget 1984; Allamandola et al. [email protected] 1985, 1989). PAHs play an important role in the photo- 2 Chastenet, Sandstrom et al. electric heating of the ISM, the efficiency of which de- Draine et al.(2007) measured qPAH in the SINGS pends on their UV absorption cross-section and grain galaxy sample and found a wide range of qPAH, from charge (Bakes & Tielens 1994; Weingartner & Draine ∼ 0.5 % in dwarf galaxies, up to almost ∼ 5 % in spi- 2001b). ral galaxies. Their results suggested a dependence be- PAHs are also widely considered to be responsible for tween the PAH fraction and the metallicity of a galaxy, the mid-IR (MIR) emission features. Their emission where qPAH drops at lower metallicity. Some stud- dominates the MIR through broad emission bands at ies have suggested a threshold in metallicity, around 3.3, 6.2, 7.7, 11.3, and 17 µm. These prominent features 12+log(O/H) ∼ 8.0 − 8.2, at which the PAH abundance trace the vibrational modes of the C–C, C–C–C, and varies drastically (Draine et al. 2007). On the other C–H bonds in PAHs, and can be used to probe the ion- hand, R´emy-Ruyer et al.(2015) find a power-law re- ization and size distribution of the PAH population (for lationship between metallicity and qPAH, rather than a review, see Tielens 2008). Because the mid-IR bands a step-function. There are several hypotheses for the from PAHs are preferentially excited by higher energy dearth of PAHs in low-metallicity environments. PAHs photons, they are often considered as a tracer of could be exposed to more intense and/or harder far-UV formation (Peeters et al. 2004). The intensity of some radiation fields due to the overall decrease in dust shield- of the emission bands even allows for detection at high ing, and suffer from a more efficient selective photo- redshift, and can be used to determine the star forma- destruction (e.g. Madden et al. 2006). PAHs could also tion earlier in the history of the Universe (Sajina et al. form in the dense ISM (Zhukovska et al. 2016), and this 2009; Siana et al. 2009; Shipley et al. 2016). At inter- process could be less effective at lower metallicities. Or mediate redshifts, they are found to trace molecular gas the low PAH abundance could be the sign of a lower effi- (Cortzen et al. 2019). PAHs are also candidate carriers ciency of forming PAH-like dust, due to particular stellar of the “2175 A˚ bump”, seen in extinction (Mathis 1994; evolution at low metallicity (Galliano et al. 2008). Steglich et al. 2010). This intriguing feature shows vari- In order to study both the metallicity trend, and dis- ations both in width and intensity between lines of sight, entangle the various local influences on the PAH fraction between galaxies and within the same galaxy (Gordon in a galaxy, we focus on two of the closest galaxies, the et al. 2003). To understand ISM thermal balance, mid- Magellanic Clouds (MCs). The Small Magellanic Cloud IR emission and UV absorption in galaxies, we must (SMC) lies at about 62 kpc (Graczyk et al. 2014), and understand the life cycle of PAHs. has a metallicity of ∼ 1/5 Z (Russell & Dopita 1992). Studies have found evidence of changes in PAH prop- The is closer, at about 50 kpc erties as a function of galaxy properties, particularly a (Walker 2012), and has a higher metallicity, of ∼ 1/2 Z PAH deficiency at low metallicity (Engelbracht et al. (Russell & Dopita 1992). The MCs therefore represent 2005; Madden et al. 2006; Draine et al. 2007; Galliano attractive targets for detailed studies of dust proper- et al. 2008; Sandstrom et al. 2010; Paradis et al. 2011; ties (e.g. Leroy et al. 2007; Bernard et al. 2008; Paradis R´emy-Ruyer et al. 2015). The intensity of the MIR fea- et al. 2009; Bot et al. 2010; Israel et al. 2010; Sandstrom tures decreases in these environments relative to the to- et al. 2010; Galliano et al. 2011; Paradis et al. 2011; tal IR emission. This suggests that there is a change Sandstrom et al. 2012; Gordon et al. 2014; Roman-Duval in the dust composition, with a lower abundance of et al. 2014; Tchernyshyov et al. 2015). Their respective the grains responsible for the MIR emission, compared metallicities bracket a threshold at which PAH proper- to the larger grains, emitting primarily at far-infrared ties are thought to vary significantly (Draine et al. 2007). (FIR) wavelengths. In the Draine & Li(2007) dust Here, we present a study of the PAH fraction across model, the PAH fraction is defined as the fraction of both MCs at 10 pc resolution, with the same dust grain the total dust mass in grains with less than 103 car- model, in order to compare to the MW and other nearby bon atoms, and is hereafter labeled qPAH. The Galactic galaxies. Thanks to the MCs proximity, we can resolve 1 diffuse neutral medium qPAH lies around 4.6 % (Li & ISM structures, such as H II regions, in the FIR. This Draine 2001; Weingartner & Draine 2001a). allows us to provide detailed maps of 10 pc scale PAH abundance, and investigate its variation as a function of the dominant ISM phase. Throughout this paper, we 1 It is difficult to estimate a systematic uncertainty on qPAH. It depends on the physics of PAHs (e.g., UV, optical and IR cross assume a constant metallicity across each galaxy. How- sections, broad continuum PAH emission, which show a large scat- ever, recent studies by Fukui et al.(2017) and Tsuge ter in their theoretical values), and modeling assumption (e.g., et al.(2019) suggested that the H I ridge of the LMC, starlight spectral shape). south of the star forming complex 30 Doradus, is mainly SMC-stripped gas from an colliding event. PAHs in the MCs 3

Our paper is laid out as follows: we first describe the remove as much of the contaminating point sources as data used in this study, in Section2. The dust grain possible, we perform a fit with the un-masked images. model used to fit the dust IR emission is detailed in Sec- Using these fit results, we mask where there is an tion3, before giving the results of the fit, and showing evident bias in the parameter maps due to a point source the variations of the PAH fraction with environment, in in the image. We mask all the selected sources by re- Section4. Finally, Section5 is dedicated to discussion placing the value within a small radius from the source and interpretation of the PAH abundance variations in with the averaged value of the local diffuse emission. the MCs. Prior to fitting we perform several additional steps of image processing. We correct the IRAC and SPIRE images with extended source factors, as suggested by 2. DATA the IRAC Instrument Handbook3 and the KINGFISH 2.1. Infrared User Guide4. Then, all observations are convolved to the We combine observations from the Spitzer Space Tele- SPIRE 500 resolution (∼ 36”) using the Aniano et al. scope (Werner et al. 2004) and the Herschel Space Ob- (2011) convolution kernels. Even though the clouds are servatory (Pilbratt et al. 2010) to cover the mid- through at relatively high Galactic latitude, the observations still far-IR spectral energy distribution (SED). In the mid- suffer from a non-negligible emission from the Milky- infrared (MIR), we use observations from the Spitzer Way cirrus. We follow Gordon et al.(2014) and Chas- Legacy program SAGE (Surveying the Agents of the tenet et al.(2017) to remove this foreground emission. Galaxy), which observed both the LMC (SAGE-LMC; We convert the foreground H I MW column density map Meixner et al. 2006) and the SMC (SAGE-SMC; Gordon to a dust column density using coefficients 1.47 × 10−3, et al. 2011). The final images produced by SAGE-SMC 9.32 × 10−4, 1.34 × 10−2, 4.28 × 10−2, 2.53 × 10−2, include deeper measurements in the main star form- 2.63 × 10−1, 1.36, 1.07, 1.85, 1.20, 0.62, 0.25 MJy/sr ing regions of the SMC, from the S3MC program (Bo- (1020 H I/cm2)−1 from 3.6 to 500 µm respectively, and latto et al. 2007). In the far-infrared (FIR) and sub- subtract the foreground cirrus. We then perform a back- millimeter range, the Herschel Key Project HERITAGE ground subtraction to get rid of residual emission from (the Herschel Inventory of the Agents of Galaxy Evolu- the cosmic infrared background, zodiacal light and mo- tion; Meixner et al. 2013, 2015) surveyed both clouds. saicing offsets. Background regions are selected by eye: This leads to a total of 11 photometric bands included we visually identify portions of the images where we can in our study, at 3.6, 4.5, 5.8, 8.0 (Spitzer, IRAC; Fazio avoid contamination from dust emission from the tar- et al. 2004), 24, 70 (Spitzer, MIPS; Rieke et al. 2004), get galaxies. In the LMC, these are chosen to be at 100, 1602 (Herschel, PACS; Poglitsch et al. 2010), 250, the edges of the images. In the SMC, we avoid the 350, and 500 µm (Herschel, SPIRE; Griffin et al. 2010). SMC Bar and wing to select the background pixels. The goal of this study is to model dust emission from We fit and subtract a tilted plane, which removes the the MIR to the sub-millimeter range. However, there gradient across the background. All images were then is a significant contribution from stars in the Spitzer projected with the final pixel grid sampling the point IRAC bands. For the average stellar populations in a spread function with approximately independent pixels, pixel, we can account for stellar emission with a simple that is lpixel ∼ 42”, which corresponds to a pixel size of assumption of a 5,000 K blackbody (Draine et al. 2007, ∼ 12 pc in the SMC and ∼ 10 pc in the LMC. After the and see next Section), but due to the proximity of the final projection, a background covariance matrix Cbkg Magellanic Clouds, occasionally a resolution element is was constructed from the background pixels, to quantify dominated by a very bright source that is not well mod- the correlations between noise in different photometric eled by a blackbody. In particular, this can be young bands (Gordon et al. 2014). stellar objects or evolved stars (Woods et al. 2011; Jones et al. 2017). To avoid contamination by these sources, 2.2. Additional data we mask out the bright point sources that show up in In this study, we are interested in possible correlations the short wavelength bands. We perform our own mask- between the fitted dust properties and other components ing of the brightest sources in these bands, in two steps. of the ISM. We use the SHASSA survey (the Southern A first set of sources is simply chosen by looking at the images, and selecting the brightest stars. Second, to 3 https://irsa.ipac.caltech.edu/data/SPITZER/docs/irac/ iracinstrumenthandbook/29/ 2 We use the PACS160 instead of MIPS160 band, for the sake 4 http://irsa.ipac.caltech.edu/data/Herschel/KINGFISH/ of resolution. docs/KINGFISH DR3.pdf 4 Chastenet, Sandstrom et al.

H-Alpha Sky Survey Atlas; Gaustad et al. 2001) to in- Table 1. Fitting parameters vestigate the spatial distribution of the ionized gas from Hα emission5. We use the smoothed maps of the LMC Parameter Unit Range Step 0 a (field 013), and SMC (field 510), at 4 angular resolution. Umin [0.1, 50] Uneven spacing The maps were projected onto the final pixel grid of our b log(qPAH) % [-1.0, 0.88] 0.0725 data set, thus oversampling the point spread function of the Hα emission data. log10(γ) - [-3.3, 0] 0.1 12 2 We use CO(J = 1 − 0) maps from the NANTEN log10(Σd)M /pc [-2.0, 0.7] 0.15

(Fukui et al. 1991) survey of the SMC (Mizuno et al. log10(Ω∗) - [-2.0, 2.7] 0.15 2001) and the LMC (Fukui et al. 2008) to trace the spa- a tial distribution of the molecular gas (∼ 30 resolution). Umin ∈ {0.1, 0.12, 0.15, 0.17, 0.2, 0.25, 0.3, 0.35, 0.4, 0.5, 0.6, 0.7, 0.8, 1.0, 1.2, 1.5, 1.7, 2.0, 2.5, 3.0, 3.5, 4.0, In Section 5.2, we will use a 3-σ detection threshold 5.0, 6.0, 7.0, 8.0, 10.0, 12.0, 15.0, 17.0, 20.0, 25.0, 30.0, in CO integrated intensity to define the “molecular gas 35.0, 40.0, 50.0}. phase”. We determine this threshold by computing the b standard deviation σ in a region where we do not find We also include qPAH = 0. any detections by eye. We will consider the molecular gas phase as every pixel above 3-σ. In the SMC, we find a 3-σ value of 0.3 K km s−1, and 0.75 K km s−1 in the SMC −2 −1 −1 LMC. Assuming αCO = 76 M pc (K km s ) , and α = 2. In the end, we have 5 free parameters: the LMC −2 −1 −1 and αCO = 10 M pc (K km s ) (Jameson et al. minimum radiation field Umin, γ the fraction of the dust 2016), this corresponds to a molecular gas surface den- mass heated by the power-law distribution of radiation −2 −2 sity of 22.8 M pc in the SMC, and 7.5 M pc in fields, the PAH fraction qPAH, the dust surface density the LMC. Σd, and the scaling parameter of stellar surface bright- ness, Ω∗, which adjusts a 5,000 K blackbody to match 3. SPECTRAL ENERGY DISTRIBUTION FITTING the observed starlight continuum in the shortest wave- METHODOLOGY length bands. See Table1 for the boundaries of each We use the Draine & Li(2007) model to fit the dust parameter (Range), and the sampling (Step). The pa- emission in our data set, from 3.6 to 500 µm (with up- rameter values ensure a sampling fine enough to resolve dates similar as those from Aniano et al. 2012). The the 1-D likelihood functions and were determined after model has a size distribution of grains divided into a several iterations of the fit. carbonaceous and a silicate component. The fraction of Using the fitted parameters, we calculate, in each pixel the dust grain mass made up of PAHs with less than j, U and fPDR, as described by Aniano et al.(2012). U 3 measures the dust-mass-weighted average radiation field 10 carbon atoms is given by the qPAH parameter. Like Draine et al.(2007) and Aniano et al.(2012) we fit the intensity, generally given by: SED using the MW RV = 3.1 grain model. This model has been found to provide good fits in low metallicity Uj = (1 − γj) Umin,j + γj  conditions by Draine et al.(2007) and Sandstrom et al.   2−αj 2−αj !  αj − 1 Umax − Umin,j  , if α 6= 1, α 6= 2; (2010) and allows us to include qPAH as a fit parameter,  1−αj 1−αj  αj − 2 U − U which is critical to our study. × max min,j In each pixel j, the dust is heated by a range of radia-  ln(Umax/Uminj) Umin,j , if α = 2. tion field intensities described by the parameter U, a di-  1 − Umin,j/Umax mensionless factor scaling the Mathis et al.(1983) 10 kpc (1) Milky Way interstellar radiation field. In each pixel, a fraction (1−γ) of the dust grain mass is heated by a radi- fPDR is the fraction of the dust luminosity produced by 2 ation field of intensity Umin. The remaining fraction γ is regions where U ≥ UPDR = 10 : heated by a power-law distribution of Umin ≤ U ≤ Umax (see equations 8, 9 and 10 of Aniano et al. 2012), with a fPDR = LPDR/Ldust (2) 7 power-law exponent α. In this study, we fix Umax = 10 , with

5 The maps are in units of dR. ln(Umax/UPDR) 6 −2 −1 −1 LPDR,j =P0,j(qPAH) Md,j γj , 1 R = 10 /4π photons cm s sr . U −1 − U−1 1 R = 5.661 × 10−18 erg s−1 cm−2 arcsec−2 min,j max (3) Ldust,j =P0,j(qPAH) Md,j Uj. PAHs in the MCs 5

Here P (q ) is the power radiated per unit dust mass, 0 PAH 102 SMC when U = 1, Md is the dust mass. LMC Aside from their pixel-to-pixel distribution, we mea- sure galaxy-average values of the fitted parameters, 1

] 10

weighted according to the dust mass or luminosity dis- r s tribution. Following Aniano et al.(2012), we define the / y dust-mass averaged value of the mass fraction of PAHs, J M [ 100

PN I qPAH,j Md,j hq i = j , (4) PAH PN j Md,j and similarly, the dust-mass averaged starlight intensity, 10 1 Fit Data PN U j Md,j 1 2 hUi = j . (5) 10 10 PN m j Md,j Wavelength [ ] The average value of fPDR, i.e. dust luminosity Figure 1. Examples of the fitting results in two pixels in weighted, is expressed as the diffuse ISM of the SMC (blue) and the LMC (orange). PN The circles mark the data SEDs with their 3-σ errors, and LPDR,j hf i = j . (6) the solid lines show the best fit model, in each case. PDR PN j Ldust,j grid, weighted by the likelihood function. The realiza- We assume that P (q ) variations are small enough 0 PAH tions render noise properties more accurately than find- pixel-to-pixel to be negligible (Aniano et al. 2012). ing only the maximum of the likelihood function. To th 3.1. Fitting and uncertainties measure values like hqPAHi and their associated 16 and th In Figure1 we show two examples of the data in the 84 percentiles, we prefer the realization method over SMC and the LMC fit by the best model (i.e. maximum the expectation value (Gordon et al. 2014; Chiang et al. likelihood), in two pixels of the diffuse ISM. The error 2018; Utomo et al. 2019), since it draws samples without bars at short wavelengths show that the errors are small marginalizing the likelihood over all but one parameter dimension. Since we compute weighted-values with both enough to strongly constrain qPAH. Residuals show that the fits are good in both galaxies at short wavelengths, qPAH and Σd, it is better to use the full-dimension likeli- especially 8 µm. At long wavelengths, the fractional hood distribution to avoid losing information. We build residuals are mostly negative, indicating that the model a large number of realization maps for each parameter generally overestimates the data (residuals peak at less simultaneously (both the fitted parameters and the cal- than 10 %). culated parameters like U and fPDR), and use them to The fitting is done with the DustBFF tool (Gordon calculate the mean in each pixel. This is what is pre- et al. 2014), which determines the n-dimensional a pos- sented in Figures2 and3. teriori likelihood distribution of the parameters. It is a From the realizations of a single pixel or group of pix- Bayesian fitting tool that uses flat priors on the param- els, we can determine its statistical error by measuring eter distributions, and a covariance matrix, built from the standard deviation from a large number of realiza- background pixels and instrument errors, to propagate tions. When combining many pixels together we can uncertainties. We use the same calibration errors on propagate these uncertainties to calculate the error on Spitzer and Herschel instruments as those in Gordon the mean values. To represent the intrinsic scatter of a parameter within a specific region of the galaxy (e.g. et al.(2014) and Chastenet et al.(2017). The diagonal th th elements of the covariance matrix include both uncor- Figures7 or8), we use the 16 and 84 percentiles of related (or statistical) and correlated errors, and the the dust mass (or dust luminosity for fPDR) weighted non-diagonal elements measure only correlated errors distribution of said parameter. For all averaged param- between the photometric bands. We refer the reader eters, e.g. hqPAHi, quoted in the rest of the paper, the to Gordon et al.(2014) for further details on the covari- statistical uncertainties are generally very small, due to the large numbers of pixels begin averaged together. We ance matrix and the DustBFF fitting. th th We use the 5-dimensional (for 5 fitting parameters) instead list the ± the 16 and 84 percentiles of the likelihood function to build realizations of the parame- dust mass- or luminosity-weighted distributions. ter maps: a realization is a sample of the 5-dimensional 6 Chastenet, Sandstrom et al.

4. RESULTS colder dust, affecting the PAH fraction. In addition, the 3 4.1. Fitting parameter maps extent of the S MC Spitzer maps used in that paper did not allow as accurate a MW foreground removal as Figures2 and3 show the results from fitting the pixels enabled by the expanded SAGE-SMC coverage. This of the SMC and the LMC maps, respectively, as well may have resulted in an oversubtraction of actual SMC as the computed values of U and fPDR, as described emission in the mid-IR bands, decreasing the Sandstrom in Section3. The overall distribution of Σ d shifts et al.(2010) qPAH value. towards higher values in the LMC than in the SMC, as In the LMC, Paradis et al.(2009) found an enhanced expected because of the larger dust-to-gas ratio of the PAH fraction with respect to the large grain abundance former (Gordon et al. 2014). A significant difference is in the stellar bar. Our results do not show a increased also seen in the qPAH parameter, which is clearly higher qPAH in this region, as other regions of the LMC show in the LMC than in the SMC (a larger version of the the same PAH fraction as in the optical bar. This differ- qPAH map is shown in Figure4). We find dust-mass- ence could be the result of including the Herschel bands, weighted average for each galaxy of and thereby more accurately measuring qPAH. It may also be the result of using a different dust model (e.g. hqSMCi = 1.0+0.3 %, and PAH −0.3 the Desert et al.(1990) model in the paper by Paradis LMC +1.4 hqPAH i = 3.3−1.3 %. et al. 2009), or a different stellar continuum model. We find values of the global PDR fraction (i.e. the Both of these values fall within the range seen by Draine fraction of the dust luminosity produced by regions et al.(2007) in the SINGS sample at the relevant metal- where U ≥ 102): licities. By studying a sample of low-metallicity galax- SMC +0.01 ies, R´emy-Ruyer et al.(2015) found a power-law relation hfPDR i = 3.29−0.02 %, and between the metallicity and the PAH fraction (see their hf LMCi = 9.17+0.04 %. Equation 5). If we apply this relation to the MCs metal- PDR −0.04 licities (SMC: 12+log(O/H) ∼ 8.1; LMC: 12+log(O/H) These values are below the average values found in NGC ∼ 8.3; Russell & Dopita 1992), we find: 628 and NGC 6946, two nearby, resolved galaxies, by Aniano et al.(2012, ∼ 11 %). However, it is not surpris- RR15 0.69 % ≤ qPAH (SMC) ≤ 3.47 %, and ing as the spatial scales in their study are coarser than RR15 for the MCs: due to limited resolution, the luminosity- 1.25 % ≤ qPAH (LMC) ≤ 6.31 %. weighted fPDR is biased towards high-luminosity values. Our average PAH fractions found in this work fall well The resulting weighted-average is then more sensitive to within these ranges as well. these high-luminosity regions, due to the blending of sig- In both galaxies over most of the area, Umin ∼ U. nal. This is one of the key results of the recent study Thus we only show the fitting parameter Umin in Fig- by Utomo et al.(2019). We do find high fPDR values in ures2 and3. We calculate the dust-mass-weighted av- star forming regions, up to ∼ 60 % in 30 Dor (LMC), +0.4 erage of U and find hUi = 1.2−0.4 in the SMC and and ∼ 50 % in N66 (SMC). +1.6 hUi = 1.6−1.0 in the LMC. Utomo et al.(2019) recently The spatial distributions of Ω∗ trace the regions with studied the distributions of mass and temperature in high stellar density (Zaritsky et al. 2002, 2004). It shows four nearby galaxies, including the SMC and the LMC, the old stellar spheroid population of the SMC, as well using a single-temperature modified blackbody model. as the optical bar in the LMC. They find the distribution of dust mass as a function of 4.2. PAHs in the Magellanic Clouds radiation field intensity peaks at values of Upeak = 1.1 in the SMC and Upeak = 1.8 in the LMC, which are In Figure4, we show the maps of the qPAH parameter consistent with our results. in the LMC (top) and the SMC (bottom). The contours Sandstrom et al.(2010) found a dust-mass-weighted are Hα emission from the SHASSA survey (Gaustad PAH fraction of ∼ 0.6 % in the SMC (compared to the et al. 2001) at level of ∼ 1.5 × 10−16 erg/s/cm2/arcsec2 1.0 % in this study). One possible explanation for our (∼ 300 dR; solid line; the reason for choosing this value lower fraction is the limited coverage of the dust emission is discussed later in this Section). The labeled circles are SED in Sandstrom et al.(2010) compared to this paper: H II regions as identified in Lopez et al.(2014); the radii here, the longest wavelength is 500 µm while it is 160 µm correspond to those given in their Table 1, and defined in Sandstrom et al.(2010). With the addition of the as the limit where they enclose 90 % of the Hα emission Herschel bands, we are able to better constrain the total of the source. There are other known H II regions in dust mass and temperature, particularly in regions with the Magellanic Clouds; however, identification of H II PAHs in the MCs 7

log(Umin) log(qPAH [%]) log( ) 0.70 0.25 0.20 0.65 1.10 0.8 0.4 0.0 0.4 0.8 3.3 2.8 2.3 1.8 1.3

1 kpc 1 kpc 1 kpc

2 log( d [M pc ]) log( * ) fPDR 2.00 1.40 0.75 0.10 0.50 1.60 0.65 0.30 1.25 2.20 0.00 0.05 0.10 0.15 0.20

1 kpc 1 kpc 1 kpc

Figure 2. Results of the fitting in the SMC for the minimum radiation field, Umin, the PAH mass fraction, qPAH (see Figure4, and Section 4.2), the weight of dust mass heated by a power-law combination of radiation field, γ, the total dust surface density, 2 Σd, the scaling of stellar surface brightness, Ω∗, and the derived fraction of dust luminosity where U > 10 , fPDR. regions and their boundaries is not straightforward, and The implications of these observations for the PAH life- we use only the ones in Lopez et al.(2014) for the sake of cycle will be discussed further in Section 5.2. We define homogeneity. From these maps we draw several conclu- the four phases as follows: sions which we discuss in the following (sub)sections: 1) • Ionized gas toward H II regions. This is sim- qPAH varies dramatically between the SMC and LMC, 2) ply defined by the pixels falling within the radii qPAH shows variations within each galaxy and 3) a pri- of H II regions from Lopez et al.(2014). In mary driver for the variation in qPAH within each galaxy the LMC, there is a clear drop in the dust-mass- appears to be correlated with the presence of ionized gas HII weighted PAH fraction, and hqPAHi reaches only as traced by Hα. +1.1 ∼ 1.8−1.3 %, i.e. slightly more than half of To investigate the variations of qPAH within and be- the galaxy average. In the SMC, the dust-mass- tween the galaxies, we make four separations with re- +0.3 weighted PAH fraction reaches ∼ 0.8−0.5 %. We spect to the ISM gas phase. In Table2, we report hqPAHi note that the harder radiation field in and near in these four phases. There is no overlap between these H II regions, which is not taken into account in values as we use only pixels that do not fall in two gas our fitting, would lead us to overestimate q in phase definitions. In short, we observe the following: PAH these regions. The values found here are therefore conservative.

• Non-H II region ionized gas. We distinguish the mol diff. neutr. out−HII HII qPAH ∼ qPAH > qPAH > qPAH (7) ionized gas inside and outside H II regions, by 8 Chastenet, Sandstrom et al.

log(Umin) log(qPAH [%]) log( ) 0.70 0.25 0.20 0.65 1.10 0.8 0.4 0.0 0.4 0.8 3.3 2.8 2.3 1.8 1.3

1 kpc 1 kpc 1 kpc

2 log( d [M pc ]) log( * ) fPDR 2.00 1.40 0.75 0.10 0.50 1.60 0.65 0.30 1.25 2.20 0.00 0.05 0.10 0.15 0.20

1 kpc 1 kpc 1 kpc

Figure 3. Results of the fitting in the LMC for Umin, qPAH, γ,Σd, and Ω∗, and the derived fPDR. selecting pixels whose Hα surface brightness is in the SMC, similar to the values in the diffuse −16 2 2 above IHα ∼ 1.5 × 10 erg/s/cm /arcsec but neutral medium. not in identified H II regions. Although this ion- ized gas is in a more diffuse phase than the gas • Diffuse neutral gas. This is defined with the in HII regions, we avoid identifying it as “diffuse regions that fall into none of the above cate- ionized gas (DIG)” due to the specific ways that gories: with Hα emission lower than ∼ 1.5 × DIG is defined in nearby galaxies (for reviews, see 10−16 erg/s/cm2/arcsec2, which also means out- Mathis 2000; Haffner et al. 2009, and references side of an H II region, and below the 3-σ CO de- therein). We discuss this further in Section 5.3. tection. We use the dust-mass-weighted PAH frac- out−HII +1.1 tion in the diffuse neutral medium as a reference We find hqPAH i = 2.9−1.2 % in the LMC and out−HII +0.3 ref +0.2 hq i = 0.9 % in the SMC. value for each cloud. We find hqPAHi = 1.1−0.3 % PAH −0.3 ref +0.6 in the SMC, and hqPAHi = 4.1−0.8 % in the LMC.

• Molecular gas. We use 12CO(J = 1 − 0) maps (Section2) to trace the molecular gas. 4.3. qPAH in H II Regions We define this phase with every pixel above The H II regions, as identified in Lopez et al.(2014), the 3-σ detection threshold, and IHα . 1.5 × appear as minima in the qPAH map, and suggest that −16 2 2 mol. 10 erg/s/cm /arcsec . We find hqPAHi ∼ PAHs are destroyed inside of H II regions. In Figure4, +1.3 mol. +0.1 4.3−0.9 % in the LMC, and hqPAHi ∼ 1.1−0.2 % we can see that most H II regions, marked by the black PAHs in the MCs 9

7

6 -66°00'00.0" N48 5 N55

4 -67°00'00.0" N4 N30 N59 3

-68°00'00.0" ] % [

2

N105 H A

N157 P Dec (J2000)

N144 q -69°00'00.0" N160 N79 1 N180 -70°00'00.0"

N191 -71°00'00.0" N206

1 kpc 0 45m00.00s 30m00.00s 15m00.00s 5h00m00.00s 4h45m00.00s RA (J2000)

3.0 1 kpc N71 2.5

-72°00'00.0" N80 N78 N76 2.0 N66 1.5 30'00.0" ] % [

1.0 H N50 A P

Dec (J2000) N63

-73°00'00.0" q N36

DEM74 0.5 N19 N84 30'00.0" N22 N90 N13 N51 N17 -74°00'00.0" 0.0 24m00.00s 12m00.00s 1h00m00.00s 0h48m00.00s RA (J2000)

SMC Figure 4. Maps of the qPAH parameter in the LMC (top) and the SMC (bottom). We find dust-mass-averaged value of hqPAH i = LMC −16 2 2 1.0 %, and hqPAH i = 3.3 %. The contours are Hα emission from the SHASSA survey at ∼ 1.5 × 10 erg/s/cm /arcsec . The labeled circles are H II regions identified by Lopez et al.(2014). See Section 4.2 for details. 10 Chastenet, Sandstrom et al.

Table 2. hqPAHi (in %) in different gas phases of each galaxy. LMC 4.0 N11 N30 SMC LMC 2.5 N44 N48 +0.3 +1.4

Global average 1.0−0.3 3.3−1.3 ] 1.6 N55

% N59 [

+0.3 +1.1 H II regions 0.8−0.5 1.8−1.3 1.0 N79 H N105 1 +0.3 +1.1 A P 0.63 Non-H II region ionized gas 0.9−0.3 2.9−1.2 N119 q N144 2 +0.2 +0.6 0.4 Diffuse neutral 1.1−0.3 4.1−0.8 N157 3 +0.1 +1.3 0.25 N160 Molecular gas 1.1−0.2 4.3−0.9 N180 1 0.16 N191 All pixels above the lower limit in Hα, excluding N206 H II regions. 0.1 0.1 1.0 10.0 2 L14 All pixels below the CO detection, the limit in r/rH II Hα, and outside of H II regions. 3All pixels with CO detection not overlapping SMC 1.0 DEM74 with pixels with N13 −16 2 2 IHα & 1.5 × 10 erg/s/cm /arcsec . N17 0.63 N19 ] N22 % [

N36

H 0.4 N50 A

P N51 labeled circles, are indeed low in PAHs, with respect to q N63 0.25 N66 the abundance in the other parts of each galaxy. N76 We are interested in understanding how H II regions N78 0.16 N80 affect the PAH fraction. In Figure5, we show the radial N84 profiles of hq i for each H II region in the SMC (bot- PAH 0.1 tom) and the LMC (top). We see that hq i drops to 0.1 1.0 10.0 PAH L14 very low values only in a handful of cases (notably, 0 % r/rH II in 30 Doradus). This is expected, as the diffuse neutral Figure 5. Radial profiles of hqPAHi in the LMC (top) and gas projected along the line-of-sight contaminates our the SMC (bottom), from the centers of each H II region, measurement of the qPAH inside the H II region. The normalized by the radius given by Lopez et al.(2014). The dust emission along the line-of-sight reflects the full col- horizontal grey lines mark the average for each galaxy. The umn through the ISM of the galaxy, not just the H II maximum radii are chosen to be slightly higher than the radii region, so we do not expect the observed qPAH to be 0 %. of the spheres of influence, in order to show its full extent. In the particular case of 30 Doradus, a result consistent with 0 % PAH fraction could mean that the actual H II tion of the ionizing stars within an H II region scales as region dominates the entire line of sight through the one would expect for the Str¨omgrensphere (Str¨omgren LMC. 1939). In this particular case, assuming a constant For each H II region, we measure the radius at which gas density the radius of the Str¨omgrensphere grows hqPAHi goes back to the global average (horizontal grey with the ionizing photon production rate (as does the line in Figure5). In this Figure, we see that hqPAHi Hα luminosity) with a power-law coefficient of 1/3. If in each H II region goes back to the galaxy average at PAHs are destroyed in ionized gas by photo-destruction different radii. We treat these radii as the “spheres of or sputtering, then one would expect the region with a influence” of each H II region on the surrounding PAH deficit of PAHs relative to the galaxy average to grow fraction. In Figure6, we plot these radii against the to- as the size of the Str¨omgrensphere. In Binder & Povich tal Hα luminosity from within that radius. We find that (2018), the authors found a relation between the pop- these radii correlate very well with Hα luminosity with ulation of stars inside star forming regions of the MW, a power-law coefficient between 0.35 and 0.4, in both and the PAH fraction. They show that a single O6 star galaxies. If the Hα surface brightness were constant, is less effective at destroying PAHs than a population growing the radius would lead to a dependence of LHα of stars that extends to O2/O3, and Wolf-Rayet stars. in r2. We see in Figure6 that is not the case. It rather Glatzle et al.(2019) showed that the growth of H II suggests that the sphere of influence on the PAH frac- regions is related to the dust content within, including PAHs in the MCs 11

± ± 104 SMC: log(r) = (0.35 0.04) log(LH ) + (-11.0 1.5) 3 SMC LMC: log(r) = (0.40 ± 0.04) log(LH ) + (-12.7 ± 1.6) 2 LMC H g v

3 A 1 /3 a 10 1 P x y q

/

H A P q Radius [''] 2 1/3 10 x y

1 0.1 10 35 36 37 38 39 40 18 17 16 15 14 10 10 10 10 10 10 2 2 log(IH [erg/s/cm /arcsec ]) LH [erg/s]

Figure 7. Correlation between hqPAHi and Hα surface Figure 6. Correlations between the total Hα luminosity brightness, for the LMC (orange circles) and the SMC (blue and the radius at which hqPAHi returns to the global average squares). The PAH fraction is normalized to the global av- from the center of each H II regions. The slope found in erage in each galaxy. Arrows indicate that the value falls both cases is close to that expected for the expansion of a below hq i = 0.1. The vertical dashed gray line marks the 1/3 PAH typical Str¨omgrensphere. We add y = x slope for compar- limit between the diffuse neutral and non-H II region ionized ison (grey dashed lines). (The empty symbols mark region gas. where the calculated radius exceeds five times the one given in Lopez et al.(2014), and are not used to fit the line. The light lines are results from bootstrapping for uncertainties on ionized gas. Given that the distribution of H I is fairly the coefficients.) flat across both galaxies (e.g. Stanimirovic et al. 2000), it is possible that the turn-over at 300 dR occurs at a PAH abundance, by impacting the ionization fronts. A point where ionized gas starts to make a large contri- more detailed study of H II regions would be interesting bution to the total surface density. The limit marked by the vertical dashed line is lower than the Hα sur- to possibly link the initial hqPAHi value at small radius in Figure5, the properties of the ionizing star(s) within face brightness in H II regions, implying that PAHs may and the expansion of H II regions. undergo destruction even in ionized gas outside of H II regions.6 4.4. qPAH and Hα A decrease of PAH fraction in diffuse ionized gas has Even outside of the H II regions, we find a clear re- been suggested by previous works in the Milky Way lation between increasing Hα emission and decreasing Dobler et al.(2009) and Dong & Draine(2011). In Dong hqPAHi. Figure7 shows hqPAHi binned as a function & Draine(2011), the authors studied the H α-to-free-free of Hα surface brightness for the SMC (blue squares) emission in the MW diffuse ionized gas, and found that and the LMC (orange circles). Both correlations tend this ratio corresponds to a lower temperature than they to reach hqPAHi ∼ 0 % at high IHα values, indepen- were able to produce with their model using the MW dently of the average value of each galaxy. This agrees diffuse ISM PAH abundance. With a model consisting with the scenario where PAHs are destroyed in H II of ionized gas, recombining gas in the process of cool- regions. In Figure7, we identify the value IHα ∼ ing, and cool neutral gas, they manage to reproduce the 1.5 × 10−16 erg/s/cm2/arcsec2 (or ∼ 300 dR) as the observations by lowering (by a factor of ∼ 3) the PAH Hα emission level at which hqPAHi starts to rapidly de- fraction in photo-ionized regions, with respect to that of crease in the LMC (dashed gray line). We use this value the global average in the ISM. Our results, constrained to plot the solid contours on Figure4. This marks a from observed SEDs, agree with a scenario where the limit where the PAH fraction drops more steeply with PAH fraction is lower in ionized gas. Hα surface brightness. We use this value as a sepa- ration between the diffuse neutral gas and ionized gas. 4.5. qPAH and the radiation field Figure7 shows that the PAH fraction decreases even at low values of Hα surface brightness, but that it drops 6 Assuming an electron, and proton, density ne ∼ np = rapidly only after it reaches the vertical dashed line −3 −16 2 2 0.5 cm , a gas temperature T = 5000 K, and a flat H I sur- (∼ 1.5 × 10 erg/s/cm /arcsec , or 300 dR). Before 20 −2 face density in the LMC ΣHI = 20 × 10 cm , our limit in Hα that value, the PAH fraction decreases only mildly. To surface brightness corresponds to a fraction of ionized gas ∼ 12 %, first order the Hα emission traces the surface density of or an ionized zone depth ∼ 150 pc. 12 Chastenet, Sandstrom et al.

A high intensity of the radiation field is often quoted U 0.251 0.398 0.631 1.0 1.58 2.51 3.98 6.31 as a cause for enhanced destruction of PAHs. Here, we 7.0 6.0 diffuse neutral non-H II diffuse ionized SMC can test that scenario by looking at hqPAHi as a func- 5.0 LMC tion of Umin, or U. In our work, we only adjust for the 4.0 radiation field intensity, and not the hardness. In the 3.0 top panel of Figure8, we show the variations of hqPAHi ] % [ in bins of U, in the diffuse neutral medium (filled sym- 2.0 H A bols), and the non-H II region ionized gas (above the P limit in Hα emission and outside of H II regions; empty q symbols). We choose to compare the diffuse neutral 1.0 medium and the non-H II region ionized gas because ionized-gas related destruction processes should not oc- .

cur in the former. In the diffuse neutral medium of the t u

e 1.0 n

SMC, hqPAHiseems to increase slightly with Ualthough . f H f i A d with increasing uncertainty. This uptick in hq i may P 0.8 PAH q

/ be the result of an overestimation due to the fact that I I 0.6 we do not take into account a change in radiation field H t H u A o hardness (see Section 5.3.1) as U increases. It is possi- P 0.4 q ble that in low-metallicity environments, changes in the -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 radiation field hardness are more important, hence the log(U) observed increase of hqPAHi in the SMC. In the diffuse neutral medium in the LMC however, we can notice a Figure 8. Top: Variations of hqPAHi as a function of U derived from the fit, for the SMC (blue), and the LMC (or- decrease of hq i at U ∼ 3. Given the scatter, there are PAH ange). Pixels that fall above the limit in Hα emission, and only minor variations of hqPAHi up to log(U) = 0.5 (U outside of H II regions, are marked with empty symbols. Pix- ∼ 3). At this point, hqPAHi decreases noticeably. Below els that fall in the diffuse neutral medium, i.e. below 300 dR this, the intensity of the radiation field does not seem to and below the CO detection threshold, are marked with filled affect the PAH fraction. As expected, in the ionized gas symbols. Bottom: Fraction of hqPAHi in the ionized medium outside of H II regions, there seems to be a decrease of with respect to hqPAHi in the neutral medium. hqPAHi with U even for U< 1. This is not surprising be- cause in defining this phase, we selected the pixels where et al.(2015) 7. Our values also sit well within the scatter hqPAHi decreases in Figure7. When we control for the of the Draine et al.(2007) results, where a steep drop of intensity of the radiation field, we can see that the gas qPAH occurs at a metallicity threshold of 12+log(O/H) phase has an impact of the PAH fraction. For an iden- ∼ 8.1. tical U, the PAH fraction does not decrease as quickly, The degree to which our results agree with either the whether it is in the diffuse neutral or the non-H II region threshold or power-law dependence, however, depends ionized gas. We also point out that, even at the lowest on the gas phase. In the MW, studies have found dif- U, hqPAHi in the ionized medium never reaches that of fuse neutral medium qPAH value of ∼ 4.6 % (Li & Draine the neutral medium. 2001; Weingartner & Draine 2001a). For a fair compar- The bottom panel of Figure8 shows the ratio of PAH ison, this should be compared to our selection of the out−HII diff.neut. fraction in different gas phases hqPAH i/hqPAH i diffuse neutral medium in the LMC and SMC listed in (same color code than the top panel). It is interesting to Table2. In this case, the power-law boundaries found notice the increasing offset between the PAH fraction in by R´emy-Ruyer et al.(2015) underestimate the PAH the ionized gas and the neutral medium as U increases fraction in the LMC. Indeed, the LMC’s diffuse neutral and that this offset is the same in both galaxies. We gas qPAH is very similar to the Milky Way’s, while the discuss these results in the context of PAH destruction SMC falls a factor of ∼ 4 lower. If we treat the diffuse in Section 5.2. neutral gas as a reference value, where most destruc- tion processes are not operating, our results strengthen 5. DISCUSSION the interpretation that there is a threshold in metallic-

5.1. Dependence of qPAH on metallicity SMC 7 The difference between the global averages hqPAH i and Their study deals with integrated SED fitting. Their qPAH hqLMCi is in good agreement with the power-law depen- parameter is then similar to a luminosity-mass weighted hqPAHi PAH in our study. dence of qPAH with metallicity found by R´emy-Ruyer PAHs in the MCs 13 ity where qPAH decreases significantly, and that is lies redistribution of carbon-rich dust mass from large grains between SMC and LMC metallicities. into small grains, and even PAH-like fragments, is signif- icant. In Slavin et al.(2015), the authors found that the 5.2. Insights into the life-cycle of PAHs from redistribution of dust mass in smaller grains is impor- comparison of hqPAHi in different phases tant in supernova remnants shocks up to ∼ 200 km s−1. We use the separation of the gas phases in the MCs A critical question, however, is whether the hot shocked and their associated PAH fraction to investigate the life- gas is able to quickly sputter and destroy any PAH frag- cycle of PAHs in the ISM. We first go over some of ments that are created by shattering. Theoretical calcu- the possible scenarios for PAH formation and how they lations by Micelotta et al.(2010b) suggest that PAHs are might affect hqPAHi in each phase. quickly destroyed by collisions with energetic particles Theoretical studies have found that the atmospheres at shock velocities higher than 100 km s−1. This work of post-AGB stars could be a favored site for the forma- leads to the conclusion that shattering of dust grains, tion of PAHs (Cherchneff et al. 1992; Sloan et al. 2007). in < 100 km s−1 shocks could be a significant source PAHs produced by AGB stars should be primarily input of PAH-like grains. Because the shattering process is into the diffuse ISM, since the old stellar distribution is directly linked to the available carbon-rich dust mass, not closely related to the current location of star-forming one would expect a higher abundance of small grain regions or dense gas (e.g. see Sandstrom et al. 2010, for fragments if there is a higher dust-to-gas ratio. Given a comparison between the SMC old stellar distribution the observed dust-to-gas ratio in the Magellanic Clouds and ISM). If PAHs are primarily formed by AGB stars (∼ 0.003, ∼ 0.0008 in the LMC and SMC, respectively; and no additional destruction or production mechanisms Leroy et al. 2011; Gordon et al. 2014; Roman-Duval et al. are operating in the diffuse neutral gas, one could argue 2014), this would be consistent with our finding of higher ref ref that the hqPAHi we have identified should reflect the ef- hqPAHi in the LMC, compared to that of the SMC. In ref ficiency of PAH formation by AGB stars. Boyer et al. the former, the higher hqPAHi, assuming destruction of (2011) and Boyer et al.(2012) studied the stellar pop- PAHs is not significant in the diffuse neutral medium, ulation of both MCs, and found that the carbon-rich could be the result of more efficient shattering of large AGB population have similar dust production rates in carbonaceous grains. the SMC and in the LMC, despite their different metal- Another hypothesis for PAHs formation is growth in licities. Although other evolved stars such as red super- the molecular phase of the ISM (e.g. Sandstrom et al. giants seem to be more effective at producing dust in 2010, 2012; Zhukovska et al. 2016). In this scenario, the the LMC, they do not play a major role in the overall lower hqPAHi in the SMC would be related to less efficient dust production by evolved stars (e.g. Boyer et al. 2011). growth processes in the molecular gas, possibly due to Therefore, without some additional destruction mecha- the lower abundance of metals to accrete onto existing mol. nism in the diffuse neutral phase, we cannot explain the grains to form PAHs. We do see that the hqPAHi is lower ref different hqPAHi between the SMC and the LMC solely in the SMC than in the LMC. Previous work by Sand- with more effective formation of PAHs in AGB atmo- strom et al.(2010) found enhanced qPAH in dense regions spheres in the LMC8. of the SMC compared to the diffuse gas phases, which Another scenario for PAH formation suggested by the- was interpreted as evidence for growth in the molecular oretical studies (e.g. Jones et al. 1996) is shattering of gas and destruction operating in the diffuse ISM. We large carbon grains leading to the production of smaller do find higher PAH fraction in the molecular gas phase carbon grain fragments. Jones et al.(1996) argued that of each galaxy with respect to the global averages, (see the shattering of carbon grains can happen at velocity Table2; LMC: 4 .3 %; SMC: 1.1 %), but the hqPAHi as low as 1 km s−1, and is the prevalent process affect- is similar in the molecular and diffuse neutral phases −1 ing dust grains at shock velocities > 100 km s . Their of each galaxy. The observation of the same hqPAHi in study focused on the grain-grain interaction in the so- the molecular and diffuse neutral gas is consistent either called warm intercloud medium, similar to our definition with formation of PAHs in the diffuse ISM and incorpo- of diffuse neutral medium. There, they found that the ration into molecular clouds or its inverse, assuming no destruction processes operate differentially in these two phases. In the same molecular gas, accretion and coagu- 8 If the PAH-like dust production is equal in both clouds, a higher production of non-PAH dust in the SMC would lead to a lation of small grains onto big grains (e.g. Stepnik et al. lower PAH fraction. Based on the same works previously quoted, 2003; K¨ohleret al. 2012) could also lead to a decrease in there is no indication than the SMC stellar population is more the observed PAH fraction. In our case, observational effective at producing non-PAH dust than in the LMC, and we evidence for that mechanism would be lower PAH frac- disfavor this possibility. 14 Chastenet, Sandstrom et al. tion in the molecular regions. Since we do not see such PAHs in a more intense radiation field could be more variation, we do not favor this possibility. easily destroyed because they are more fragile, due to Observational studies have shown that the PAH size their ionization state. Ionized PAHs are less stable, and distribution and properties are sensitive to the local ra- more prone to losing, for example, H atoms (Montillaud diation field and gas ionization. This is seen in the vicin- et al. 2013). Weingartner & Draine(2001c) showed that ity of Milky Way H II regions and photo-dissociation the charge of small carbonaceous√ grains varies depend- regions (e.g. Bern´eet al. 2007; Compi`egneet al. 2007; ing on a parameter G T/ne, where G is the radiation Arab et al. 2012; Peeters et al. 2017) as well as in nearby field intensity, T the temperature of the gas, and ne the galaxies (Gordon et al. 2008; Paradis et al. 2011; Rela˜no electron density. To test whether PAH ionization may et al. 2016). PAH destruction can be accomplished in lead to higher destruction rates, the different behavior several ways, mediated by interaction with energetic of the neutral and ionized gas phases at fixed U from photons, sputtering by particles in hot gas, or chemi- Figure8 is of interest. In a bin of U, we control for cal reactions. A long-standing hypothesis to explain the the radiation field intensity, and it is the same in both PAH deficit at low metallicity is that PAHs are more the diffuse neutral and the non-H II diffuse ionized gas. readily destroyed in such conditions (e.g. Madden et al. Weingartner & Draine(2001c) showed that the gas tem- 2006; Galliano et al. 2008). A possible scenario for en- perature does not significantly affect the ionization state hanced destruction of PAHs at low metallicity is photo- of PAHs. The only parameter remaining then is ne. If destruction by the radiation field, either due to increased the electron density were lower in the non-H II ionized intensity or hardness. There are a number of observa- gas, compared to the diffuse neutral gas, then PAHs tional studies that have found correlations that agree could be more ionized in this medium, facilitating their with such a scenario. For example, Madden et al.(2006) destruction. While this seems counter-intuitive, given and Gordon et al.(2008) found that the PAH features that the non-H II region ionized gas is defined by its ion- disappear as the radiation field hardness increases in ization state, the difference in overall density between H II regions (measured from mid-IR neon and sulfur line the ionized and neutral phases and the low, but non- ratios). Theoretical studies have also shown that PAHs negligible, fractional ionization in the neutral medium are subject to sputtering and fragmentation in ionized could lead to the situation where ne is lower for the ion- gas (Micelotta et al. 2010a; Bocchio et al. 2012) because ized gas. Li & Draine(2001) showed that PAHs larger of electronic- and/or atomic interactions. There, projec- than 7 A˚ are more ionized in the Milky Way’s warm tile particles can reach a high velocity with respect to neutral medium than in the warm ionized medium. In out−HII that of the grains, leading to catastrophic collisions. In that case, the lower hqPAH i compared to that of the general, the regions of the galaxy with the hardest ra- diffuse neutral medium could not be explained by an diation fields will also be those where ionized gas exists easier destruction of PAHs because of their ionization. and the overall intensity of the radiation field is higher. Based on the comparison between the gas phases in Therefore, to distinguish between the potential destruc- the two galaxies, we could explain the offset in hqPAHi tion mechanisms, we need to attempt to separate the between the LMC and SMC with: (i) fragmentation of effects of these quantities. We do so using Hα emission large grains in the diffuse medium leading to more PAHs as a tracer of ionized gas, and our determination of U in the LMC due to its higher dust content, (ii) forma- as a tracer of the intensity of the radiation field, as pre- tion of PAHs in the molecular ISM and their injection sented in Section 4.5 and Figure8. We do not have a in the diffuse gas, assuming there is no preferential de- direct tracer for the hardness of the radiation field cov- struction in either molecular or diffuse gas. The forma- ering the full extent of the galaxies. We note the similar tion of PAHs in AGB stars cannot explain the difference out−HII diff.neut. trends of the ratios hqPAH i/hqPAH i of the SMC in PAH fraction between the SMC and the LMC. The (blue) and the LMC (orange), in the bottom panel of destruction mechanisms (e.g. photo-destruction by the Figure8. They suggest that metallicity does not have radiation field), suggested to be more effective in lower an impact on the relative efficiencies of the destruction metallicity galaxies, do not differ significantly between processes between the diffuse neutral medium and the the SMC and the LMC in this work (see bottom panel ionized gas outside of H II regions (in the same bin of of Figure8). radiation field). We do not test, however, the impact of metallicity on the global amount of each gas phase 5.3. Caveats in a galaxy, which might have an impact on the overall 5.3.1. Impact of the modeled radiation field destruction of PAHs. The Draine & Li(2007) model uses the radiation field described by Mathis et al.(1983), for the Milky Way at PAHs in the MCs 15 the galactocentric distance DG = 10 kpc. This limits 6. CONCLUSIONS the possible variations allowed in our fitting. Specif- We fit the dust emission SED in the Small and Large ically, the hardness of the radiation field is constant. Magellanic Clouds, using photometry from Spitzer and However, it is expected that the relative proportion of Herschel (3.6 to 500 µm), with the dust emission model UV vs optical photons will not be the same in H II from Draine & Li(2007). We provide maps of each regions where the ionizing stars produce more UV pho- fitting parameter at a 42” pixel size, i.e. ∼ 10 pc in the tons. Some studies have made adjustments to the ra- LMC and ∼ 12 pc in the SMC. diation field to address this issue (e.g. Galliano et al. We especially focus on the spatial distribution of the 2005; Salgado et al. 2016). In M33, Rela˜noet al.(2016) PAH fraction qPAH, the fraction of dust mass in grains studied the dust content of H II regions, and specifically with less than 103 carbon atoms. We find a global changed the radiation field to one with more UV pho- SMC dust-mass-weighted PAH fraction hqPAH i = 1.0 % and tons. They showed that this leads to lower PAH frac- LMC hqPAH i = 3.3 %, both lower than the Milky Way diffuse tion, by a factor up to 3. Paradis et al.(2011) studied neutral medium value (4.6 %). the impact of the radiation field on dust emission fitting We measure the PAH fraction in different gas phases, in the LMC, by adding a 4 Myr stellar population to distinguished by 12CO(J = 1 − 0) and Hα emission (Ta- the Mathis et al.(1983) radiation field. As expected, ble2) or the lack thereof. We use the diffuse neutral they found that, in H II regions, using a harder radi- −16 2 2 medium (IHα . 1.5 × 10 erg/s/cm /arcsec and no ation field leads to a decrease in the PAH abundance molecular gas detection) as a reference value and dis- estimation. This is linked to the PAH being more sen- cuss the relative hqPAHi in each gas phase. We find that sitive to the UV-visible part of the incident radiation the PAH fraction in the LMC diffuse neutral medium field (Li & Draine 2002). Harder radiation, i.e. more (4.1 %) is similar to that of the MW diffuse neutral energetic photons, would enhance their MIR emission. ISM, while in the SMC it is substantially lower than in In our study, this would make the offset between the the MW (1.1 %). The galaxy-averaged hqPAHi are con- SMC and the LMC (low-metallicity stars produce more sistent with both a power-law dependence of qPAH with UV photons; e.g. Eldridge et al. 2008), and the differ- metallicity (R´emy-Ruyer et al. 2015), as well as the ex- ence between diffuse-to-ionized PAH fraction, even more istence of a threshold around 12 + log(O/H) ∼ 8.1 at dramatic (this was also shown by Paradis et al. 2011, in which the PAH abundance changes rapidly (e.g., Draine the LMC). et al. 2007). But hqPAHi in the diffuse neutral gas favors the latter hypothesis. 5.3.2. Diffuse ionized gas We find evidence that hqPAHi is systematically low in The definition of ionized gas in the MW (called “warm all identified H II regions (Figures4,7). Additionally, ionized gas”, WIM; e.g. Reynolds 1984; Madsen et al. we show that the sphere of influence of the H II regions 2006), and in extra-galactic studies (called “diffuse ion- on the PAH fraction grows as one would expect for a ized gas”, DIG; e.g. Zurita et al. 2000) is a delicate sub- Str¨omgrensphere (Figure6). ject (see also McKee & Begelman 1990; Mathis 2000; We investigate possible metallicity-dependent PAH Haffner et al. 2009). Given the difficulty defining “dif- formation and destruction scenarios to explain the fuse ionized gas”, we decided not to make assumptions higher PAH fraction in the neutral medium of the LMC on the exact definition of the ionized medium outside compared to the SMC: of H II regions. A finer description of the local condi- tions of the ionized gas would require measurements of - We find higher hqPAHi in molecular gas with re- the electron/proton density, incident radiation field, and spect to the galaxy global averages, but similar to gas temperature. those of the diffuse neutral medium. This is con- sistent with formation of PAHs in the molecular 5.3.3. Metallicity variations across the Magellanic Clouds gas and injection in the diffuse neutral medium, or vice-versa. In this study, we assume the metallicity to be constant across each galaxy. If there were metallicity variations, - In each galaxy, we find a clear trend of qPAH with there should also be variations in qPAH. If the metal- Hα emission on global scale (Figure7). We find licity variations are not correlated with ISM phase, we a limit in Hα luminosity (or ionized gas fraction) would only expect to see enhanced scatter. If they were at which the PAH abundance starts to decrease, correlated with ISM phase, the results would likely be and interpret this as a destruction of PAHs in the different. We do not see a good reason to assume any ionized gas. This limit is lower than the typical metallicity variations would be ISM phase correlated. Hα emission in an H II regions, suggesting that 16 Chastenet, Sandstrom et al.

the ionized medium, even outside of H II regions, us the opportunity to study the PAH fraction at low still affects the PAH fraction through destruction metallicity in the early Universe. processes. - A radiation field intensity about twice that of the solar neighborhood seems to be enough to affect We are very grateful to the referee for their care- the PAH fraction, in both the neutral and ion- ful reading and their comments which greatly improved ized gas phase (Figure8). However, the ionized the paper. JC wishes to thank Jean-Philippe Bernard, medium always shows a lower PAH fraction than Olivier Bern´e and D´eborah Paradis for fruitful dis- the neutral medium, even at equal intensity of the cussions. The work of JC, KS, IC, AKL, and DU radiation field. is supported by NASA ADAP grants NNX16AF48G - We find that formation of PAHs through the frag- and NNX17AF39G and National Science Foundation mentation of large grains is a plausible explanation grant No. 1615728. The work of AKL and DU is par- for the higher hqPAHi in the LMC compared to the tially supported by the National Science Foundation un- SMC diffuse neutral medium. der Grants No. 1615105, 1615109, and 1653300. This paper used the Southern H-Alpha Sky Survey Atlas Future work will investigate the variations of the PAH (SHASSA), which is supported by the National Science fraction with metallicity in resolved, nearby galaxies. Foundation. The launch of the James Webb Space Telescope will allow for the detection of PAHs at higher redshifts, and give Facility: Herschel, Spitzer

REFERENCES

Allamandola, L. J., Tielens, A. G. G. M., & Barker, J. R. Cherchneff, I., Barker, J. R., & Tielens, A. G. G. M. 1992, 1985, ApJL, 290, L25 ApJ, 401, 269 —. 1989, ApJS, 71, 733 Chiang, I.-D., Sandstrom, K. M., Chastenet, J., et al. 2018, Aniano, G., Draine, B. T., Gordon, K. D., & Sandstrom, K. ApJ, 865, 117 2011, PASP, 123, 1218 Compi`egne,M., Abergel, A., Verstraete, L., et al. 2007, Aniano, G., Draine, B. T., Calzetti, D., et al. 2012, ApJ, A&A, 471, 205 756, 138 Cortzen, I., Garrett, J., Magdis, G., et al. 2019, MNRAS, Arab, H., Abergel, A., Habart, E., et al. 2012, A&A, 541, 482, 1618 A19 Desert, F.-X., Boulanger, F., & Puget, J. L. 1990, A&A, Bakes, E. L. O., & Tielens, A. G. G. M. 1994, ApJ, 427, 822 237, 215 Bernard, J.-P., Reach, W. T., Paradis, D., et al. 2008, AJ, Dobler, G., Draine, B., & Finkbeiner, D. P. 2009, ApJ, 699, 136, 919 1374 Bern´e,O., Joblin, C., Deville, Y., et al. 2007, A&A, 469, Dong, R., & Draine, B. T. 2011, ApJ, 727, 35 575 Draine, B. T., & Li, A. 2007, ApJ, 657, 810 Binder, B. A., & Povich, M. S. 2018, ApJ, 864, 136 Draine, B. T., Dale, D. A., Bendo, G., et al. 2007, ApJ, Bocchio, M., Micelotta, E. R., Gautier, A.-L., & Jones, 663, 866 A. P. 2012, A&A, 545, A124 Eldridge, J. J., Izzard, R. G., & Tout, C. A. 2008, MNRAS, Bolatto, A. D., Simon, J. D., Stanimirovi´c,S., et al. 2007, 384, 1109 ApJ, 655, 212 Engelbracht, C. W., Gordon, K. D., Rieke, G. H., et al. Bot, C., Ysard, N., Paradis, D., et al. 2010, A&A, 523, A20 2005, ApJL, 628, L29 Boyer, M. L., Srinivasan, S., van Loon, J. T., et al. 2011, Fazio, G. G., Hora, J. L., Allen, L. E., et al. 2004, ApJS, AJ, 142, 103 154, 10 Boyer, M. L., Srinivasan, S., Riebel, D., et al. 2012, ApJ, Fukui, Y., Ogawa, H., Kawabata, K., Mizuno, A., & 748, 40 Sugitani, K. 1991, in IAU Symposium, Vol. 148, The Bron, E., Le Bourlot, J., & Le Petit, F. 2014, A&A, 569, Magellanic Clouds, ed. R. Haynes & D. Milne, 105–106 A100 Fukui, Y., Tsuge, K., Sano, H., et al. 2017, PASJ, 69, L5 Chastenet, J., Bot, C., Gordon, K. D., et al. 2017, A&A, Fukui, Y., Kawamura, A., Minamidani, T., et al. 2008, 601, A55 ApJS, 178, 56 PAHs in the MCs 17

Galliano, F., Dwek, E., & Chanial, P. 2008, ApJ, 672, 214 Mathis, J. S. 1994, ApJ, 422, 176 Galliano, F., Madden, S. C., Jones, A. P., Wilson, C. D., & —. 2000, ApJ, 544, 347 Bernard, J.-P. 2005, A&A, 434, 867 Mathis, J. S., Mezger, P. G., & Panagia, N. 1983, A&A, Galliano, F., Hony, S., Bernard, J.-P., et al. 2011, A&A, 128, 212 536, A88 McKee, C. F., & Begelman, M. C. 1990, ApJ, 358, 392 Gaustad, J. E., McCullough, P. R., Rosing, W., & Van Meixner, M., Gordon, K. D., Indebetouw, R., et al. 2006, Buren, D. 2001, PASP, 113, 1326 AJ, 132, 2268 Glatzle, M., Ciardi, B., & Graziani, L. 2019, MNRAS, 482, Meixner, M., Panuzzo, P., Roman-Duval, J., et al. 2013, 321 AJ, 146, 62 Gordon, K. D., Clayton, G. C., Misselt, K. A., Landolt, —. 2015, AJ, 149, 88 A. U., & Wolff, M. J. 2003, ApJ, 594, 279 Micelotta, E. R., Jones, A. P., & Tielens, A. G. G. M. Gordon, K. D., Engelbracht, C. W., Rieke, G. H., et al. 2010a, A&A, 510, A37 2008, ApJ, 682, 336 —. 2010b, A&A, 510, A36 Gordon, K. D., Meixner, M., Meade, M. R., et al. 2011, AJ, Mizuno, N., Rubio, M., Mizuno, A., et al. 2001, PASJ, 53, 142, 102 L45 Gordon, K. D., Roman-Duval, J., Bot, C., et al. 2014, ApJ, Montillaud, J., Joblin, C., & Toublanc, D. 2013, A&A, 552, 797, 85 A15 Graczyk, D., Pietrzy´nski,G., Thompson, I. B., et al. 2014, Paradis, D., Paladini, R., Noriega-Crespo, A., et al. 2011, ApJ, 780, 59 ApJ, 735, 6 Griffin, M. J., Abergel, A., Abreu, A., et al. 2010, A&A, Paradis, D., Reach, W. T., Bernard, J.-P., et al. 2009, AJ, 518, L3 138, 196 Haffner, L. M., Dettmar, R.-J., Beckman, J. E., et al. 2009, Peeters, E., Bauschlicher, Jr., C. W., Allamandola, L. J., Reviews of Modern Physics, 81, 969 et al. 2017, ApJ, 836, 198 Israel, F. P., Wall, W. F., Raban, D., et al. 2010, A&A, Peeters, E., Spoon, H. W. W., & Tielens, A. G. G. M. 2004, 519, A67 ApJ, 613, 986 Jameson, K. E., Bolatto, A. D., Leroy, A. K., et al. 2016, Pilbratt, G. L., Riedinger, J. R., Passvogel, T., et al. 2010, ApJ, 825, 12 A&A, 518, L1 Jones, A. P., Tielens, A. G. G. M., & Hollenbach, D. J. Poglitsch, A., Waelkens, C., Geis, N., et al. 2010, A&A, 1996, ApJ, 469, 740 518, L2 Jones, O. C., Woods, P. M., Kemper, F., et al. 2017, Rela˜no,M., Kennicutt, R., Lisenfeld, U., et al. 2016, A&A, MNRAS, 470, 3250 595, A43 K¨ohler,M., Stepnik, B., Jones, A. P., et al. 2012, A&A, R´emy-Ruyer, A., Madden, S. C., Galliano, F., et al. 2015, 548, A61 A&A, 582, A121 Le Bourlot, J., Le Petit, F., Pinto, C., Roueff, E., & Roy, Reynolds, R. J. 1984, ApJ, 282, 191 F. 2012, A&A, 541, A76 Rieke, G. H., Young, E. T., Engelbracht, C. W., et al. 2004, Le Page, V., Snow, T. P., & Bierbaum, V. M. 2009, ApJ, ApJS, 154, 25 704, 274 Roman-Duval, J., Gordon, K. D., Meixner, M., et al. 2014, Leger, A., & Puget, J. L. 1984, A&A, 137, L5 ApJ, 797, 86 Leroy, A., Bolatto, A., Stanimirovic, S., et al. 2007, ApJ, Russell, S. C., & Dopita, M. A. 1992, ApJ, 384, 508 658, 1027 Sajina, A., Spoon, H., Yan, L., et al. 2009, ApJ, 703, 270 Leroy, A. K., Bolatto, A., Gordon, K., et al. 2011, ApJ, Salgado, F., Bern´e,O., Adams, J. D., et al. 2016, ApJ, 830, 737, 12 118 Li, A., & Draine, B. T. 2001, ApJ, 554, 778 Sandstrom, K. M., Bolatto, A. D., Draine, B. T., Bot, C., —. 2002, ApJ, 572, 232 & Stanimirovi´c,S. 2010, ApJ, 715, 701 Lopez, L. A., Krumholz, M. R., Bolatto, A. D., et al. 2014, Sandstrom, K. M., Bolatto, A. D., Bot, C., et al. 2012, ApJ, 795, 121 ApJ, 744, 20 Madden, S. C., Galliano, F., Jones, A. P., & Sauvage, M. Shipley, H. V., Papovich, C., Rieke, G. H., Brown, M. J. I., 2006, A&A, 446, 877 & Moustakas, J. 2016, ApJ, 818, 60 Madsen, G. J., Reynolds, R. J., & Haffner, L. M. 2006, Siana, B., Smail, I., Swinbank, A. M., et al. 2009, ApJ, 698, ApJ, 652, 401 1273 18 Chastenet, Sandstrom et al.

Slavin, J. D., Dwek, E., & Jones, A. P. 2015, ApJ, 803, 7 Weingartner, J. C., & Draine, B. T. 2001a, ApJ, 548, 296 Sloan, G. C., Jura, M., Duley, W. W., et al. 2007, ApJ, —. 2001b, ApJS, 134, 263 664, 1144 —. 2001c, ApJS, 134, 263 Stanimirovic, S., Staveley-Smith, L., van der Hulst, J. M., Werner, M. W., Roellig, T. L., Low, F. J., et al. 2004, et al. 2000, MNRAS, 315, 791 ApJS, 154, 1 Steglich, M., J¨ager, C., Rouill´e,G., et al. 2010, ApJL, 712, Wolfire, M. G., Hollenbach, D., McKee, C. F., Tielens, L16 A. G. G. M., & Bakes, E. L. O. 1995, ApJ, 443, 152 Stepnik, B., Abergel, A., Bernard, J.-P., et al. 2003, A&A, Woods, P. M., Oliveira, J. M., Kemper, F., et al. 2011, 398, 551 MNRAS, 411, 1597 Str¨omgren,B. 1939, ApJ, 89, 526 Zaritsky, D., Harris, J., Thompson, I. B., & Grebel, E. K. Tchernyshyov, K., Meixner, M., Seale, J., et al. 2015, ApJ, 2004, AJ, 128, 1606 811, 78 Zaritsky, D., Harris, J., Thompson, I. B., Grebel, E. K., & Tielens, A. G. G. M. 2008, ARA&A, 46, 289 Massey, P. 2002, AJ, 123, 855 Tsuge, K., Sano, H., Tachihara, K., et al. 2019, ApJ, 871, 44 Zhukovska, S., Dobbs, C., Jenkins, E. B., & Klessen, R. S. Utomo, D., Chiang, I., Leroy, A. K., Sandstrom, K. M., & 2016, ApJ, 831, 147 Chastenet, J. 2019, arXiv e-prints, arXiv:1902.08629 Zurita, A., Rozas, M., & Beckman, J. E. 2000, A&A, 363, 9 Walker, A. R. 2012, Ap&SS, 341, 43