Researcher profile

Ilse de Looze

Ilse de Looze contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
17works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2026arXiv

REBELS-MOSFIRE: Weak CIII] Emission is Typical Among Extremely UV-bright, Massive Galaxies at $z\sim7$

We present Keck/MOSFIRE H-band spectroscopic measurements covering the [CIII]1907, CIII]1909 doublet for a sample of 8 z~7 spectroscopically-confirmed star-forming galaxies drawn from the Reionization Era Bright Emission Line Survey (REBELS). This REBELS-MOSFIRE sample is notable for its bright median UV luminosity (Muv=-22.5 AB) and large median stellar mass (log(Mstar/Msun)=9.2). Although three sources show tentative evidence of a CIII] detection, we obtain no confident detections for any of the 8 REBELS-MOSFIRE sources. The median [CIII]1907+CIII]1909 3-sigma upper limit in equivalent width (EW) for the REBELS-MOSFIRE sample is 6.5 AA, and a stack of their H-band MOSFIRE spectra yields a non-detection with an associated 3-sigma upper limit of 2.6 AA. These upper limits fall significantly below the CIII] EW measured in a composite spectrum of representative z~7 star-forming galaxies, as well as those measured for notable early star-forming galaxies such as GN-z11, GHZ2, GS-z12, and RXCJ2248-ID. The lack of strong CIII] emission can be understood within the context of the stellar populations of the REBELS galaxies, as well as the ionization conditions and gas-phase metallicity implied by rest-frame optical spectroscopic properties ([OIII]+Hb EWs, and [OIII]5007/[OII]3727 and [NeIII]3869/[OII]3727 line ratios). The REBELS-MOSFIRE sample represents the higher-mass, higher-metallicity, lower-excitation tail of the z~7 galaxy population, whose ionizing properties must be fully characterized to constrain the role of star-forming galaxies during cosmic reionization.

preprint2022arXiv

B-fields in Star-Forming Region Observations (BISTRO): Magnetic Fields in the Filamentary Structures of Serpens Main

We present 850 $μ$m polarimetric observations toward the Serpens Main molecular cloud obtained using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields In STar-forming Region Observations (BISTRO) survey. These observations probe the magnetic field morphology of the Serpens Main molecular cloud on about 6000 au scales, which consists of cores and six filaments with different physical properties such as density and star formation activity. Using the histogram of relative orientation (HRO) technique, we find that magnetic fields are parallel to filaments in less dense filamentary structures where $N_{H_2} < 0.93\times 10^{22}$ cm$^{-2}$ (magnetic fields perpendicular to density gradients), while being perpendicular to filaments (magnetic fields parallel to density gradients) in dense filamentary structures with star formation activity. Moreover, applying the HRO technique to denser core regions, we find that magnetic field orientations change to become perpendicular to density gradients again at $N_{H_2} \approx 4.6 \times 10^{22}$ cm$^{-2}$. This can be interpreted as a signature of core formation. At $N_{H_2} \approx 16 \times 10^{22}$ cm$^{-2}$ magnetic fields change back to being parallel to density gradients once again, which can be understood to be due to magnetic fields being dragged in by infalling material. In addition, we estimate the magnetic field strengths of the filaments ($B_{POS} = 60-300~μ$G)) using the Davis-Chandrasekhar-Fermi method and discuss whether the filaments are gravitationally unstable based on magnetic field and turbulence energy densities.

preprint2022arXiv

Dust masses for a large sample of core-collapse supernovae from optical emission line asymmetries: dust formation on 30-year timescales

Modelling the red-blue asymmetries seen in the broad emission lines of core-collapse supernovae (CCSNe) is a powerful technique to quantify total dust mass formed in the ejecta at late times ($>5$ years after outburst) when ejecta dust temperatures become too low to be detected by mid-IR instruments. Following our success in using the Monte Carlo radiative transfer code DAMOCLES to measure the dust mass evolution in SN~1987A and other CCSNe, we present the most comprehensive sample of dust mass measurements yet made with DAMOCLES, for CCSNe aged between four and sixty years after outburst. Our sample comprises of multi-epoch late-time optical spectra taken with the Gemini GMOS and VLT X-Shooter spectrographs, supplemented by archival spectra. For the fourteen CCSNe that we have modelled, we confirm a dust mass growth with time that can be fit by a sigmoid curve which is found to saturate beyond an age of $\sim30$ years, at a mass of 0.23$^{+0.17}_{-0.12}$ M$_\odot$. An expanded sample including dust masses found in the literature for a further eleven CCSNe and six CCSN remnants, the dust mass at saturation is found to be 0.42$^{+0.09}_{-0.05}$~M$_\odot$. Uncertainty limits for our dust masses were determined from a Bayesian analysis using the affine invariant Markov Chain Monte Carlo ensemble sampler emcee with DAMOCLES. The best-fitting line profile models for our sample all required grain radii between 0.1 and 0.5 $μ$m. Our results are consistent with CCSNe forming enough dust in their ejecta to significantly contribute to the dust budget of the Universe.

preprint2022arXiv

The ALMA REBELS Survey: Specific Star-Formation Rates in the Reionization Era

We present specific star-formation rates for 40 UV-bright galaxies at $z\sim7-8$ observed as part of the Reionization Era Bright Emission Line Survey (REBELS) ALMA large program. The sSFRs are derived using improved measures of SFR and stellar masses, made possible by measurements of far-infrared (FIR) continuum emission and [CII]-based spectroscopic redshifts. For each source in the sample, we derive stellar masses from SED fitting and total SFRs from calibrations of the UV and FIR emission. The median sSFR is $18_{-5}^{+7}$ Gyr$^{-1}$, significantly larger than literature measurements lacking constraints in the FIR. The increase in sSFR reflects the larger obscured SFRs we derive from the dust continuum relative to that implied by the UV+optical SED. We suggest that such differences may reflect spatial variations in dust across these luminous galaxies, with the component dominating the FIR distinct from that dominating the UV. We demonstrate that the inferred stellar masses (and hence sSFRs) are strongly-dependent on the assumed star formation history (SFH) in reionization-era galaxies. When large sSFR galaxies are modeled with non-parametric SFHs, the derived stellar masses can increase by an order of magnitude relative to constant star formation models, owing to the presence of a significant old stellar population that is outshined by the recent burst. The [CII] line widths in the largest sSFR systems are often very broad, suggesting dynamical masses that are easily able to accommodate the dominant old stellar population suggested by non-parametric models. Regardless of these systematic uncertainties in the derived parameters, we find that the sSFR increases rapidly toward higher redshifts for massive galaxies ($9.6<\log(\rm M_*/M_{\odot})<9.8$), with a power law that goes as $(1+z)^{1.7\pm0.3}$, broadly consistent with expectations from the evolving baryon accretion rates.

preprint2022arXiv

The First Maps of $κ_{d}$ -- the Dust Mass Absorption Coefficient -- in Nearby Galaxies, with DustPedia

The dust mass absorption coefficient, $κ_{d}$, is the conversion function used to infer physical dust masses from observations of dust emission. However, it is notoriously poorly constrained, and it is highly uncertain how it varies, either between or within galaxies. Here we present the results of a proof-of concept study, using the DustPedia data for two nearby face-on spiral galaxies M74 (NGC 628) and M83 (NGC 5236), to create the first ever maps of $κ_{d}$ in galaxies. We determine $κ_{d}$ using an empirical method that exploits the fact that the dust-to-metals ratio of the interstellar medium is constrained by direct measurements of the depletion of gas-phase metals. We apply this method pixel-by-pixel within M74 and M83, to create maps of $κ_{d}$. We also demonstrate a novel method of producing metallicity maps for galaxies with irregularly-sampled measurements, using the machine learning technique of Gaussian process regression. We find strong evidence for significant variation in $κ_{d}$. We find values of $κ_{d}$ at 500 $μ$m spanning the range 0.11-0.25 ${\rm m^{2}\,kg^{-1}}$ in M74, and 0.15-0.80 ${\rm m^{2}\,kg^{-1}}$ in M83. Surprisingly, we find that $κ_{d}$ shows a distinct inverse correlation with the local density of the interstellar medium. This inverse correlation is the opposite of what is predicted by standard dust models. However, we find this relationship to be robust against a large range of changes to our method - only the adoption of unphysical or highly unusual assumptions would be able to suppress it.

preprint2021arXiv

AlFoCS + F3D II: unexpectedly low gas-to-dust ratios in the Fornax galaxy cluster

We combine observations from ALMA, ATCA, MUSE, andHerschel to study gas-to-dust ratios in 15 Fornax cluster galaxies detected in the FIR/sub-mm by Herschel and observed by ALMA as part of the ALMA Fornax Cluster Survey (AlFoCS). The sample spans a stellar mass range of 8.3 $\leq$ log (M$_*$ / M$_\odot$) $\leq$ 11.16, and a variety of morphological types. We use gas-phase metallicities derived from MUSE observations (from the Fornax3D survey) to study these ratios as a function of metallicity, and to study dust-to-metal ratios, in a sub-sample of nine galaxies. We find that gas-to-dust ratios in Fornax galaxies are systematically lower than those in field galaxies at fixed stellar mass/metallicity. This implies that a relatively large fraction of the metals in these Fornax systems is locked up in dust, which is possibly due to altered chemical evolution as a result of the dense environment. The low ratios are not only driven by HI deficiencies, but H$_2$-to-dust ratios are also significantly decreased. This is different in the Virgo cluster, where low gas-to-dust ratios inside the virial radius are driven by low HI-to-dust ratios, while H$_2$-to-dust ratios are increased. Resolved observations of NGC1436 show a radial increase in H$_2$-to-dust ratio, and show that low ratios are present throughout the disc. We propose various explanations for the low H$_2$-to-dust ratios in the Fornax cluster, including the more efficient stripping of H$_2$ compared to dust, more efficient enrichment of dust in the star formation process, and altered ISM physics in the cluster environment.

preprint2021arXiv

The Case for Thermalization as a Contributor to the [CII] Deficit

The [CII] deficit, which describes the observed decrease in the ratio of [CII] 158 micron emission to continuum infrared emission in galaxies with high star formation surface densities, places a significant challenge to the interpretation of [CII] detections from across the observable universe. In an attempt to further decode the cause of the [CII] deficit, the [CII] and dust continuum emission from 18 Local Volume galaxies has been split based on conditions within the interstellar medium where it originated. This is completed using the Key Insights in Nearby Galaxies: a Far-Infrared Survey with Herschel (KINGFISH) and Beyond the Peak (BtP) surveys and the wide-range of wavelength information, from UV to far-infrared emission lines, available for a selection of star-forming regions within these samples. By comparing these subdivided [CII] emissions to isolated infrared emission and other properties, we find that the thermalization (collisional de-excitation) of the [CII] line in HII regions plays a significant role in the deficit observed in our sample.

preprint2020arXiv

A Complete Catalogue of Dusty Supernova Remnants

We search for far-infrared (FIR) counterparts of known supernova remnants (SNRs) in the Galactic plane (360 degrees in longitude and b = +/- 1 deg ) at 70 - 500 micron with Herschel. We detect dust signatures in 39 SNRs out of 190, made up of 13 core-collapse supernovae (CCSNe), including 4 Pulsar Wind Nebulae (PWNe), and 2 Type Ia SNe. A further 24 FIR detected SNRs have unknown types. We confirm the FIR detection of ejecta dust within G350.1-0.3, adding to the known sample of ~10 SNRs containing ejecta dust. We discover dust features at the location of a radio core at the centre of G351.2+0.1, indicating FIR emission coincident with a possible Crab-like compact object, with dust temperature and mass of Td = 45.8 K and Md = 0.18 solar mass, similar to the PWN G54.1+0.3. We show that the detection rate is higher among young SNRs. We produce dust temperature maps of 11 SNRs and mass maps of those with distance estimates, finding dust at temperatures 15 < Td < 40 K. If the dust is heated by shock interactions the shocked gas must be relatively cool and/or have a low density to explain the observed low grain temperatures.

preprint2020arXiv

AlFoCS + Fornax3D: resolved star formation in the Fornax cluster with ALMA and MUSE

We combine data from ALMA and MUSE to study the resolved (~300 pc scale) star formation relation (star formation rate vs. molecular gas surface density) in cluster galaxies. Our sample consists of 9 Fornax cluster galaxies, including spirals, ellipticals, and dwarfs, covering a stellar mass range of ~10^8.8 - 10^11 M_Sun. CO(1-0) and extinction corrected Halpha were used as tracers for the molecular gas mass and star formation rate, respectively. We compare our results with Kennicutt (1998) and Bigiel et al. (2008). Furthermore, we create depletion time maps to reveal small-scale variations in individual galaxies. We explore these further in FCC290, using the &#39;uncertainty principle for star formation&#39; (Kruijssen & Longmore, 2014a) to estimate molecular cloud lifetimes, which we find to be short (<10 Myr) in this galaxy. Galaxy-averaged depletion times are compared with other parameters such as stellar mass and cluster-centric distance. We find that the star formation relation in the Fornax cluster is close to those from Kennicutt (1998) and Bigiel et al. (2008}), but overlaps mostly with the shortest depletion times predicted by Bigiel et al. (2008). This slight decrease in depletion time is mostly driven by dwarf galaxies with disturbed molecular gas reservoirs close to the virial radius. In FCC90, a dwarf galaxy with a molecular gas tail, we find that depletion times are a factor >~10 higher in its tail than in its stellar body.

preprint2020arXiv

ALMA Observations of the Molecular Clouds in NGC 625

We present the highest resolution (1&#34;) $^{12}$CO(1-0) observations of molecular gas in the dwarf starburst galaxy NGC 625, to date, obtained with ALMA. Molecular gas is distributed in discrete clouds within an area of $0.4$ kpc$^2$ and does not have well-ordered large-scale motions. We measure a total molecular mass in NGC 625 of $5.3\times 10^6 M_\odot$, assuming a Milky Way CO-to-H$_2$ conversion factor. We use the CPROPS package to identify molecular clouds and measure their properties. The 19 resolved CO clouds have a median radius of 20 pc, a median line width 2.5 km s$^{-1}$, and a median surface density of 169 $M_\odot$pc$^{-2}$. Larson scaling relations suggest that molecular clouds in NGC 625 are mostly in virial equilibrium. Comparison of our high-resolution CO observations with a star formation rate map, inferred from ancillary optical observations observations, suggests that about 40% of the molecular clouds coincide with the brightest HII regions. These bright HII regions have a range of molecular gas depletion timescales, all within a factor of $\sim 3$ of the global depletion time in NGC 625 of 106-134 Myr. The highest surface density molecular clouds towards the southwest of the galaxy, in a region we call the Butterfly, do not show strong star formation activity and suggest a depletion time scale longer than 5 Gyr.

preprint2020arXiv

High-resolution, 3D radiative transfer modelling III. The DustPedia barred galaxies

Context: Dust in late-type galaxies in the local Universe is responsible for absorbing approximately one third of the energy emitted by stars. It is often assumed that dust heating is mainly attributable to the absorption of UV and optical photons emitted by the youngest (<= 100 Myr) stars. Consequently, thermal re-emission by dust at FIR wavelengths is often linked to the star-formation activity of a galaxy. However, several studies argue that the contribution to dust heating by much older stars might be more significant. Advances in radiation transfer (RT) simulations finally allow us to actually quantify the heating mechanisms of diffuse dust by the stellar radiation field. Aims: As one of the main goals in the DustPedia project, we have constructed detailed 3D stellar and dust RT models for nearby galaxies. We analyse the contribution of the different stellar populations to the dust heating in four face-on barred galaxies: NGC1365, M83, M95, and M100. We aim to quantify the fraction directly related to young stars, both globally and on local scales, and to assess the influence of the bar on the heating fraction. Results: We derive global attenuation laws for each galaxy and confirm that galaxies of high sSFR have shallower attenuation curves and weaker UV bumps. On average, 36.5% of the bolometric luminosity is absorbed by dust. We report a clear effect of the bar structure on the radial profiles of the dust-heating fraction by the young stars, and the dust temperature. We find that the young stars are the main contributors to the dust heating, donating, on average ~59% of their luminosity to this purpose throughout the galaxy. This dust-heating fraction drops to ~53% in the bar region and ~38% in the bulge region where the old stars are the dominant contributors to the dust heating. We also find a strong link between the heating fraction by the young stars and the sSFR.

preprint2020arXiv

High-resolution, 3D radiative transfer modelling. II. The early-type spiral galaxy M81

Interstellar dust absorbs stellar light very efficiently and thus shapes the energetic output of galaxies. Studying the impact of different stellar populations on the dust heating remains hard because it requires decoupling the relative geometry of stars and dust, and involves complex processes as scattering and non-local dust heating. We aim to constrain the relative distribution of dust and stellar populations in the spiral galaxy M81 and create a realistic model of the radiation field that describes the observations. Investigating the dust-starlight interaction on local scales, we want to quantify the contribution of young and old stellar populations to the dust heating. We aim to standardise the setup and model selection of such inverse radiative transfer simulations so this can be used for comparable modelling of other nearby galaxies. We present a semi-automated radiative transfer modelling pipeline that implements the necessary steps such as the geometric model construction and the normalisation of the components through an optimisation routine. We use the Monte Carlo radiative transfer code SKIRT to calculate a self-consistent, panchromatic model of the interstellar radiation field. By looking at different stellar populations independently, we can quantify to what extent different stellar age populations contribute to the dust heating. Our method takes into account the effects of non-local heating. We obtain a realistic 3D radiative transfer model of the face-on galaxy M81. We find that only 50.2\% of the dust heating can be attributed to young stellar populations. We confirm a tight correlation between the specific star formation rate and the heating fraction by young stellar populations, both in sky projection and in 3D, also found for radiative transfer models of M31 and M51. We conclude that... (abridged)

preprint2020arXiv

Nonparametric galaxy morphology from UV to submm wavelengths

We present the first nonparametric morphological analysis of a set of spiral galaxies from UV to submm wavelengths. Our study is based on high-quality multi-wavelength imaging for nine well-resolved spiral galaxies from the DustPedia database, combined with nonparametric morphology indicators calculated in a consistent way using the {\tt{StatMorph}} package. We measure the half-light radius, the concentration index, the asymmetry index, the smoothness index, the Gini coefficient and the $M_{20}$ indicator in various wavebands from UV to submm wavelengths, as well as in stellar mass, dust mass and star formation rate maps. We find that the interstellar dust in galaxies is distributed in a more extended, less centrally concentrated, more asymmetric, and more clumpy way than the stars. This is particularly evident when comparing morphological indicators based on the stellar mass and dust mass maps. This should serve as a warning sign against treating the dust in galaxies as a simple smooth component. We argue that the nonparametric galaxy morphology of galaxies from UV to submm wavelengths is an interesting test for cosmological hydrodynamics simulations.

preprint2020arXiv

Reproducing the Universe: a comparison between the EAGLE simulations and the nearby DustPedia galaxy sample

We compare the spectral energy distributions (SEDs) and inferred physical properties for simulated and observed galaxies at low redshift. We exploit UV-submillimetre mock fluxes of ~7000 z=0 galaxies from the EAGLE suite of cosmological simulations, derived using the radiative transfer code SKIRT. We compare these to ~800 observed galaxies in the UV-submillimetre range, from the DustPedia sample of nearby galaxies. To derive global properties, we apply the SED fitting code CIGALE consistently to both data sets, using the same set of ~80 million models. The results of this comparison reveal overall agreement between the simulations and observations, both in the SEDs and in the derived physical properties, with a number of discrepancies. The optical and far-infrared regimes, and the scaling relations based upon the global emission, diffuse dust and stellar mass, show high levels of agreement. However, the mid-infrared fluxes of the EAGLE galaxies are overestimated while the far-UV domain is not attenuated enough, compared to the observations. We attribute these discrepancies to a combination of galaxy population differences between the samples, and limitations in the subgrid treatment of star-forming regions in the EAGLE-SKIRT post-processing recipe. Our findings show the importance of detailed radiative transfer calculations and consistent comparison, and provide suggestions for improved numerical models.

preprint2019arXiv

The ISM scaling relations in DustPedia late-type galaxies: A benchmark study for the Local Universe

The purpose of this work is the characterization of the main scaling relations between all the ISM components (dust, atomic/molecular/total gas), gas-phase metallicity, and other galaxy properties, such as Mstar and galaxy morphology, for late-type galaxies in the Local Universe. This study is performed by extracting late-type galaxies from the entire DustPedia sample and by exploiting the large and homogeneous dataset available thanks to the DustPedia project. The sample consists of 436 galaxies with morphological stage from T = 1 to 10, Mstar from 6 x 10^7 to 3 x 10^11 Msun, SFR from 6 x 10^(-4) to 60 Msun/yr, and 12 + log(O/H) from 8 to 9.5. The scaling relations involving the molecular gas are studied by assuming both a constant and a metallicity-dependent XCO. The analysis has been performed by means of the survival analysis technique. We confirm that the dust mass correlates very well with the total gas mass, and find -- for the first time -- that the dust mass correlates better with the atomic gas mass than the molecular one. We characterize important mass ratios such as gas fraction, molecular-to-atomic gas mass ratio, dust-to-total gas mass ratio (DGR), and dust-to-stellar mass ratio. Only the assumption of a metallicity-dependent XCO reproduces the expected decrease of the DGR with increasing morphological stage and decreasing gas-phase metallicity, with a slope of about 1. DGR, gas-phase metallicity, and the dust-to-stellar mass ratio are, for our galaxy sample, directly linked to the galaxy morphology. The molecular-to-atomic gas mass ratio and the DGR show a positive correlation for low molecular gas fractions, but for molecular gas rich galaxies this trend breaks down. This trend has never been found previously, to our knowledge. It provides new constraints for theoretical models of galaxy evolution and a reference for high-redshift studies.

preprint2019arXiv

Using [CII] 158 micron Emission From Isolated ISM Phases as a Star-Formation Rate Indicator

The brightest observed emission line in many star-forming galaxies is the [CII] 158 micron line, making it detectable up to z~7. In order to better understand and quantify the [CII] emission as a tracer of star-formation, the theoretical ratio between the [NII] 205 micron emission and the [CII] 158 micron emission has been employed to empirically determine the fraction of [CII] emission that originates from the ionized and neutral phases of the ISM. Sub-kiloparsec measurements of the [CII] 158 micron and [NII] 205 micron line in nearby galaxies have recently become available as part of the Key Insights in Nearby Galaxies: a Far Infrared Survey with Herschel (KINGFISH) and Beyond the Peak (BtP) programs. With the information from these two far-infrared lines along with the multi-wavelength suite of KINGFISH data, a calibration of the [CII] emission line as a star formation rate indicator and a better understanding of the [CII] deficit are pursued. [CII] emission is also compared to PAH emission in these regions to compare photoelectric heating from PAH molecules to cooling by [CII] in the neutral and ionized phases of the ISM. We find that the [CII] emission originating in the neutral phase of the ISM does not exhibit a deficit with respect to the infrared luminosity and is therefore preferred over the [CII] emission originating in the ionized phase of the ISM as a star formation rate indicator for the normal star-forming galaxies included in this sample.

preprint2018arXiv

A Catalogue of Galactic Supernova Remnants in the far-infrared: revealing ejecta dust in pulsar wind nebulae

We search for far-infrared (FIR) counterparts of known supernova remnants (SNRs) in the Galactic plane (10 deg <| l |< 60 deg) at 70-500 micron using the Herschel Infrared Galactic Plane Survey (Hi-GAL). Of 71 sources studied, we find that 29 (41 per cent) SNRs have a clear FIR detection of dust emission associated with the SNR. Dust from 8 of these is in the central region, and 4 indicate pulsar wind nebulae (PWNe) heated ejecta dust. A further 23 have dust emission in the outer shell structures which is potentially related to swept up material. Many Galactic SNe have dust signatures but we are biased towards detecting ejecta dust in young remnants and those with a heating source (shock or PWN). We estimate the dust temperature and mass contained within three PWNe, G11.2-0.3, G21.5-0.9, and G29.7-0.3 using modified blackbody fits. To more rigorously analyse the dust properties at various temperatures and dust emissivity index beta, we use point process mapping (PPMAP). We find significant quantities of cool dust (at 20-40 K) with dust masses of Md = 0.34 +/- 0.14 solar mass, Md = 0.29 +/- 0.08 solar mass, and Md = 0.51 +/- 0.13 solar mass for G11.2-0.3, G21.5-0.9, and G29.7-0.3 respectively. We derive the dust emissivity index for the PWN ejecta dust in G21.5-0.3 to be beta = 1.4 +/- 0.5 compared to dust in the surrounding medium where beta = 1.8 +/- 0.1.