<figure><img src="image_1.png" /></figure>

HAL
open science

# Far-reaching dust distribution in galaxy discs

M. w. l. Smith, S. a. Eales, I. de Looze, M. Baes, G. j. Bendo, S. Bianchi, M. Boquien, A. Boselli, V. Buat, L. Ciesla, et al.

<figure><img src="image_2.png" /></figure>

## ► To cite this version:

M. w. l. Smith, S. a. Eales, I. de Looze, M. Baes, G. j. Bendo, et al. Far-reaching dust distribution in galaxy discs. Monthly Notices of the Royal Astronomical Society: Letters, 2016, 462, pp.331-344. 10.1093/mnras/stw1611. hal-01433405

HAL Id: hal-01433405

https://hal.science/hal-01433405

Submitted on 9 Sep 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_3.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_4.png" /></figure>

# Far-reaching dust distribution in galaxy discs

Matthew W. L. Smith,¹* Stephen A. Eales,¹* Ilse De Looze,²,³ Maarten Baes,²

George J. Bendo,⁴ Simone Bianchi,⁵ Médéric Boquien,³,⁶ Alessandro Boselli,⁷

Veronique Buat,⁷ Laure Ciesla,⁸,⁹ Marcel Clemens,¹⁰ David L. Clements,¹¹

Asantha R. Cooray,¹² Luca Cortese,¹³ Jonathan I. Davies,¹ Jacopo Fritz,²,¹⁴

Haley L. Gomez,¹ Thomas M. Hughes,²,¹⁵ Oskar Ł. Karczewski,¹⁶ Nanyao Lu,¹⁷

Seb J. Oliver,¹⁶ Aurélie Remy-Ruyer,¹⁸ Luigi Spinoglio¹⁹ and Sebastien Viaene²

*Affiliations are listed at the end of the paper*

Accepted 2016 July 4. Received 2016 July 4; in original form 2015 December 6

## ABSTRACT

In most studies of dust in galaxies, dust is only detected from its emission to approximately the optical radius of the galaxy. By combining the signal of 110 spiral galaxies observed as part of the Herschel Reference Survey, we are able to improve our sensitivity by an order of magnitude over that for a single object. Here we report the direct detection of dust from its emission that extends out to at least twice the optical radius. We find that the distribution of dust is consistent with an exponential at all radii with a gradient of ~-1.7 dex $R_{25}^{-1}$. Our dust temperature declines linearly from ~25 K in the centre to 15 K at $R_{25}$ from where it remains constant out to ~2.0 $R_{25}$. The surface density of dust declines with radius at a similar rate to the surface density of stars but more slowly than the surface density of the star-formation rate. Studies based on dust extinction and reddening of high-redshift quasars have concluded that there are substantial amounts of dust in intergalactic space. By combining our results with the number counts and angular correlation function from the SDSS, we show that with Milky Way-type dust we can explain the reddening of the quasars by the dust within galactic discs alone. Given the uncertainties in the properties of any intergalactic dust, we cannot rule out its existence, but our results show that statistical investigations of the dust in galactic haloes that use the reddening of high-redshift objects must take account of the dust in galactic discs.

**Key words:** galaxies: ISM – galaxies: spiral – submillimetre: ISM.

## 1 INTRODUCTION

Recent studies of the extinction of quasars (Ménard et al. 2010a) and the reddening of galaxies (Peek, Ménard & Corrales 2015) imply that approximately 50 per cent of the interstellar dust in the Universe lies outside galactic discs. However, we still know remarkably little about how far out the dust in galactic discs extends.

Spiral galaxies have huge discs of atomic hydrogen that extend much further than the standard optical radius (Bigiel & Blitz 2012, $R_{25}$, the radius where the optical brightness corresponds to a B-band brightness of 25 mag arcsec$^{-2}$), but the atomic hydrogen is unprocessed material that may have fallen on the galaxies from outside. Studies into the extent of molecular gas using the carbon monoxide (CO) line, even when combining the results from a very

large number of galaxies, have only managed to trace the molecular gas out to $1.1 R_{25}$ (Schruba et al. 2011).

There have been many attempts to detect dust around nearby galaxies by using both its absorption and reflection properties in the optical, and its emission properties in the far-infrared. Observations of edge-on galaxies have been used to look for dust above the plane of the galaxy. Hodges-Kluck & Bregman (2014) found that UV emission extends between 5 and 20 kpc from the mid-plane of the galaxies in their sample, and was consistent with the emission being light from the disc scattered by dust above the disc. They found that the halo dust has an exponential form with scale-heights ~2.5–5 kpc, much larger than typical scaleheights of stellar thick discs (0.3–1.5 kpc; Hodges-Kluck & Bregman 2014). Bocchio et al. (2016) combined *Herschel* and *Spitzer* data for the edge-on galaxy NGC 891, finding extended emission above the disc with a scaleheight of 1.44 kpc.

The extent of dust discs has been investigated by looking for the effects of dust extinction of background sources (e.g. quasars,

<footnote>* E-mail: matthew.smith@astro.cf.ac.uk (MWLS); Steve.Eales@astro.cf.ac.uk (SAE)</footnote>
background galaxies) through the outskirts of a nearby galaxy. Hol-
werda et al. (2009) have shown in a pair of occulting galaxies that
dust extinction can be reliably detected out to $1.5 R_{25}$. Attempts to
detect dust in the outer disc from its emission were undertaken with
the early infrared space telescopes, *IRAS* and *ISO*. However, they
were limited by sensitivity, wavelength coverage and resolution.
Nelson, Zaritsky & Cutri (1998) obtained a $2\sigma$ detection at $1.5 R_{25}$
using *IRAS* by a stacking analysis, and Alton et al. (1998) used *ISO*
observations of a sample of eight nearby galaxies to show that the
dust has a larger exponential scalelength than the optical emission.
Recently *Spitzer* and *Herschel* data have been used to create radial
profiles of the dust for the nearby galaxies in the SINGS/KINGFISH
sample (Muñoz-Mateos et al. 2009; Hunt et al. 2015). The KING-
FISH radial profiles presented in Hunt et al. (2015) typically extend
between 1.0 and $1.5 R_{25}$ at 250 µm, across a variety of morpholog-
ical types.

In this paper, we aim to go beyond the typical detection radius
of dust for an individual galaxy, by combining the signal from a
large sample of galaxies. Since the dust in the outskirts of galaxies
is likely to be very cool ($\le 15$ K) with low surface brightness,
the *Herschel Space Observatory* (Pilbratt et al. 2010) is the only
telescope with enough sensitivity to have the potential to directly
detect this emission. We have selected galaxies from the Herschel
Reference Survey (HRS; Boselli et al. 2010), the largest targeted
survey of nearby galaxies with *Herschel*. These galaxies all have
either spiral or irregular morphology and lie face-on, giving us the
best possible physical resolution. In Section 2 we present the data
used for this analysis. Section 3 describes our method for creating
the averaged radial profiles of the dust and other galaxy components.
Section 4 presents our results, and the conclusions are presented in
Section 5.

## 2 THE SAMPLE

The HRS is a volume-limited sample of 322 galaxies with distances
between 15 and 25 Mpc (Boselli et al. 2010). To ensure a reasonable
detection rate with Herschel, the HRS also has a 2MASS K-band
magnitude limit of $K_{\text{Stot}} \le 12$ for late-type galaxies. Observations
with the SPIRE instrument (Griffin et al. 2010) at 250, 350 and
500 µm and with the PACS instrument (Poglitsch et al. 2010) at
100 and 160 µm exist for all the galaxies. 83 objects fall within
the Virgo cluster and as such have been observed by the Herschel
Virgo Cluster Survey (HeViCS; Davies et al. 2010, 2012). For all
late-type galaxies the HRS observations were designed to cover an
area of at least 1.5 times the optical diameter with a minimum size
of 4 arcmin $\times$ 4 arcmin. In practice, the SPIRE images generally
cover a much larger area around each galaxy, allowing us to trace
the dust to larger radii. Each observation of a late-type galaxy has
an instrumental $1\sigma$ sensitivity of 5.64, 5.65, and 6.60 mJy beam$^{-1}$
(0.57, 0.32, and 0.18 MJy sr$^{-1}$) for the 250, 350 and 500 µm bands,
respectively. The instrumental noise is comparable to the confusion
noise of 5.8, 6.3, and 6.8 mJy beam$^{-1}$ (Nguyen et al. 2010) for
the 250, 350 and 500 µm bands, respectively. Full details about the
survey can be found in Boselli et al. (2010).

## 3 METHODOLOGY

### 3.1 Radial profiles and averaging

To obtain the highest sensitivity to extended dust we average the
emission from a large sample of HRS galaxies to create an average

surface-brightness profile. We included galaxies in the sample if
they met the following criteria:

(i) The galaxy morphology must be classified as a spiral or irreg-
ular. The classifications are the same as used in other HRS studies
(Cortese et al. 2012a). Ellipticals and S0s were not included as they
are known to contain much less dust and have much smaller dust
discs than spirals and irregulars (Smith et al. 2012a).

(ii) The galaxy cannot have an inclination greater than 60°. The inclinations were calculated from the major and minor diameters measured by the HRS team (Cortese et al. 2012b) using images from the Sloan Digital Sky Survey (SDSS).

(iii) The galaxy must have been detected by Herschel with emis-
sion that is clearly extended (Ciesla et al. 2012).

(iv) There must not be another extended Herschel source that
could affect the measured radial profile (i.e. the sources must be
separated by a distance of at least 6.0 R<sub>25</sub>).

(v) The Herschel image must not have significant contamination
from diffuse emission from Galactic interstellar dust.

(vi) The image of each galaxy must extend to a distance of at least 2.5 R<sub>25</sub> from the centre of the galaxy to make it possible to remove any residual background emission that is not from the galaxy itself.

Once all the criteria had been applied, there were 117 galaxies,
providing an order of magnitude increase in sensitivity over what is
possible for a single object.

We obtained a surface-brightness profile for each galaxy by cal-
culating the distance to each pixel from the centre of the galaxy
and then expressing this distance as a fraction of the optical radius
($R_{25}$). We calculated this distance using the coordinates of the pixel,
the optical radius, the inclination (estimated from the eccentricity),
the distance to the galaxy (Cortese et al. 2012a) and the optical
centre. A line of fixed radius would appear as an ellipse on a SPIRE
map. We extended the profile out until a radial annulus reaches the
edge of the image. Before calculating the average profile for the
whole sample, we removed a constant offset from each galaxy’s
profile using the mean value of the surface brightness between 2.2
$\le R/R_{25} < 2.7$ (a couple of profiles do not extend as far as $2.7 R_{25}$,
but contain enough pixels for an accurate background subtraction).
This choice of background subtraction means by construction we
are only sensitive to dust within $2.2 R_{25}$. We then placed all the pix-
els from all images into radial bins and calculated the mean value
in each bin and its error ($\sigma/\sqrt{N}$, where N is the number of
data points and $\sigma$ is the standard deviation). To ensure we account for
all uncertainties (e.g. background variations), we also ran a simula-
tion by carrying out an identical analysis using the same number of
positions and the same radial bins as in our real analysis at random
positions on the HeViCS map (but avoiding Virgo cluster objects).
We then performed the same stacking analysis, and re-ran the pro-
cess a 1000 times. By measuring the standard deviations in each
bin, we get a robust estimate for the minimum error in each radial
bin. We then use this uncertainty if it is larger than the error above.

The results are shown in Fig. 1. The red-dashed lines in Fig. 1
show the $3\sigma$ detection level found from the background simulations
described in the previous paragraph. The median optical radius for
the sample is 1.25 arcmin ($R_{25}$ ranges from 0.33 to 4.83 arcmin)
and the median distance is 17.0 Mpc. Therefore, the median value
of $R_{25}$ corresponds to a physical radius of 6.2 kpc.

### 3.2 Reliability

To check that our average profile is representative of the sample of
galaxies and not just dominated by a small subset of objects, we
<figure><img src="image_5.png" /><figcaption>**Figure 1.** The average emission from dust plotted against the radial distance for the late-type galaxies in the Herschel Reference Survey. Before averaging, the radial distance has been normalized for each galaxy by the optical radius ($R_{25}$) of the galaxy. The blue data points show the average radial profile for each wavelength. The wavelength and number of galaxies that has been used in each plot are shown in the top-right corner. The green data points represent the profile of a point source processed using the same method and scaled by the median value of $R_{25}$ for the sample. The dashed red line shows the $3\sigma$ sensitivity limit found from simulations using the background regions of the Herschel Virgo Cluster Survey map.</figcaption></figure>
performed a Monte Carlo test. We randomly split the sample into four subsets and measured a profile for each subset as described above. For each profile we measured the radius where the surface brightness falls below 0.002 mJy arcsec⁻² at 250 µm. We repeated this process 10 000 times and found no subset where the profile is significantly steeper than the one we show in Fig. 1, showing our result is statistically robust.

An important question is whether at large radii, where the surface brightness is very low, we are effectively measuring the point spread function (PSF) of the observations. To test this, we used SPIRE calibration observations of Neptune. We found the average surface-brightness profile of the PSF using the same procedure as for the galaxies, scaling it to the same physical units as the galaxies using the median $R_{25}$ of the sample. The Neptune profile has been normalized to the second data point of the averaged galaxy profile and is shown in green in Fig. 1. In the figure we see that in all bands the average galaxy profile is significantly higher at all radii than the expected flux from a point source. The shallow slope of the galaxy profile compared to the PSF precludes the possibility that dust detected at large *R* is from the Airy rings of the brighter inner parts of the galaxy. Using the median $R_{25}$ to scale the Neptune profile probably overestimates the effect of the PSF as the galaxies with larger angular sizes contribute more pixels to the average profile.

As a final check that we are not simply detecting at large radii the wings of the PSF, we repeated the analysis for the 45 galaxies with $R_{25} > 1.5$ arcmin (the largest is 4.83 arcmin). We find a much larger difference between the average surface-brightness profile and the wings of the PSF (Fig. 2). Fig. 3 shows the morphological-type and stellar mass distributions for the entire HRS, our initial sample of 117 galaxies, and our sample of galaxies with $R_{25} > 1.5$ arcmin. Our initial sample of 117 galaxies has a similar distribution to the entire HRS for both morphological type (excluding early-types) and stellar mass, but our sample with $R_{25} > 1.5$ arcmin is biased towards higher stellar masses. As mentioned in Section 2, 83 of our galaxies fall within the Virgo cluster, and so the galaxies may be deficient in dust because they are deficient in H I (Cortese et al. 2010). However, we saw little difference in the results when we split the sample into classes based on measurements of the H I deficiency. To minimize the effect of the PSF, we used this sample for the rest of the analysis described in this paper. The median optical radius for this sample is 2.1 arcmin and the median distance is 17.0 Mpc. Therefore, the median value of $R_{25}$ corresponds to a physical radius of 10.3 kpc.

### 3.3 Model fitting procedure

To fit analytic models to the measured average surface-brightness profiles in Fig. 2, we need to take into account the PSF, which is large at submillimetre wavelengths. We fitted three different models to the 250–500 µm emission:

(i) A single exponential. This has two free parameters: a scale-length and the normalization of the function.

(ii) Double exponential (sum of two exponentials). This has four free parameters: the two normalizations and scalelength of each function.

(iii) Broken exponential (an exponential where at a given radius the scalelength changes). This has four free parameters: the overall normalization, the scalelength of each exponential, and the radius dividing the regions fitted by the two functions.

The fitting procedure was similar to a standard profile fitting procedure with a chi-squared minimization except that we performed a full 2D-convolution by projecting our model on to an image,

convolving the model image with the PSF, creating a model radial profile and then comparing it to the observed radial profile. To decide what is the best model, we compared the reduced chi-squared value of the three models. The fits were performed over the region in which there is at least a $3\sigma$ detection at 250 µm, except we did not fit the model to the central two data points (equivalent to central 30 arcsec) to ensure there is no effect from the nuclei. We found that the broken-exponential model provided a significantly better fit at 250 µm but at 350 and 500 µm there was no statistical reason to prefer the more complex models over the single exponential (Fig. 2). Virtually all iterations of the Monte Carlo test described in Section 3.2 show evidence for the shallower slope at large radii suggesting we are not being influenced by a small subset of galaxies.

Before deriving dust temperatures, we first convolved all images to the same resolution as the SPIRE 500 µm images which have the lowest angular resolution of all our images (using the method described in Aniano et al. 2011). We then calculated empirical average radial profiles in the SPIRE bands for the 45 HRS galaxies with $R_{25} > 1.5$ arcmin, using the same method as above. Dust temperatures tend to increase towards the centres of galaxies (Smith et al. 2010, 2012b; Foyle et al. 2012; Mentuch Cooper et al. 2012), which means that data at shorter wavelengths than the SPIRE wavelengths are necessary for accurate temperature measurements. We determined average surface-brightness profiles at 100 and 160 µm from the images of the HRS galaxies made with the PACS camera (Poglitsch et al. 2010; Cortese et al. 2014) on Herschel. However, when estimating temperatures we only used PACS points at $R < 0.5R_{25}$, since background variations in the PACS images appear to dominate at larger radii. While only using three data points to find a temperature is not ideal, at $R > 0.5R_{25}$ the dust is colder and so it becomes easier to estimate a dust temperature from the three SPIRE measurements alone than at smaller radii.

We estimated dust temperatures by fitting a modified blackbody to the five or three flux measurements for each radial bin, taking full account of the correlated flux errors between the different Herschel bands (Smith et al. 2012b). As we were limited to only three data points for $R > 0.5R_{25}$ we had to assume a constant dust-emissivity index ($\beta$). We chose to use $\beta = 2$ which is the typical value found in resolved studies of nearby galaxies (Eales et al. 2012; Smith et al. 2012b; Hughes et al. 2014; Kirkpatrick et al. 2014). The uncertainties on the dust temperature were calculated by a Monte Carlo approach; in each radial bin every data point was re-sampled using its value and uncertainty (taken correct account of the part of the calibration uncertainty that is correlated between bands) and the modified blackbody fit performed again. From the distribution of temperature values obtained from the 1000 runs of the Monte Carlo simulation, we estimated the uncertainty in the dust temperature (typically ~2 K). Our temperature profile is shown in Fig. 4. We fitted a simple model to the data points, where the dust temperature decreases linearly with radius until it reaches a minimum. We then derived the dust surface density by combining our best-fitting model for the dust emission and our model of the dust temperature. We assume a value of 0.192 m² kg⁻¹ at 350 µm (Draine 2003) for the dust absorption coefficient ($\kappa_v$).

This value is chosen to match the values used in other HRS and HeViCS papers, and is based on Milky Way-type dust. For this work we have assumed that the value of $\kappa$ is constant out to a radius of 2.0 $R_{25}$, as metallicity gradients in galaxies tend to be flat (Moustakas et al. 2010; Sánchez et al. 2014) and that we are measuring dust in the discs of galaxies. However, as very little is known about dust in galaxy outskirts it is possible that the value of $\kappa$ could have a significant systematic change with radius. The absolute
<figure><img src="image_6.png" /><figcaption>**Figure 2.** The best-fitting models to the measured average surface-brightness profiles of the sample of large galaxies (45 objects). The curved red dashed line shows the $3\sigma$ limit (Section 3.1) with the blue points showing the measurements that are above the limit. As in Fig. 1, the green data points represent the profile of a point source, processed using the same method and scaled by the median value of $R_{25}$ for the sample of large galaxies. The three models are described in Section 3.3 and are only fit to the blue data points. The colour key for the three models is shown in the legend. Only one model is shown at 350 and 500 µm because the other more complex models did not provide a significantly better fit to the data (see description in Section 3.3). The dashed and the solid lines show the models before and after they have been convolved with the PSF.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 3.** The distributions of morphology and stellar mass for the entire HRS (blue), our initial sample of 117 late-type galaxies (green) and the sample of 45 galaxies with $R_{25} > 1.5$ arcmin (red). The morphological classification is the Goldmine (Gavazzi et al. 2003) classification system, where types 3–10 represent Sa to Sdm galaxies, 11–14 are irregular/peculiar galaxies and 14–17 are BCD’s. Morphological types −3 to 2 represent early-type galaxies which are not used in this study; the grey-dashed line shows the boundary between early-type and late-type galaxies. The stellar masses are taken from Cortese et al. (2012b). Note that 10 galaxies of the HRS do not have stellar mass determinations, with 7 and 1 of these in the initial and large sample, respectively, and so are not included in the plot.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 4.** The relationship between dust temperature and galactocentric radius from fitting a modified blackbody to the spectral energy distribution of the dust emission (blue points). The dashed line shows a model in which the dust temperature declines linearly with radius until it reaches a minimum. The green line shows the prediction of a model of how the dust in a disc is heated by starlight (Section 4.1).</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>**Figure 5.** The surface density of dust plotted against radius, obtained by combining our fitted dust emission models with the temperature model shown in Fig. 4. The solid lines show our best estimates of the dust surface density for our best-fitting dust emission models at 250 and 500 µm, with the shaded regions showing the $1\sigma$ uncertainties. The dashed cyan and yellow lines are the best-fitting exponential models to the empirical dust surface-density distributions. The grey vertical dashed line shows the radius below which the data points are not fit to ensure there is no effect from galactic nuclei. The green data point is the one overlapping data point from the quasar reddening study of Ménard et al. (2010a), and the green dotted line is their estimate of the relationship between dust surface density and radius, which varies with radius as $r^{-0.8}$ (green line). If Ménard et al. (2010a) does represent the true distribution of dust our disagreement in surface densities may be the result of the range of radii that contribute to the Ménard et al. (2010a) data point or our background subtraction.</figcaption></figure>

value of $\kappa_v$ is also an important consideration when we compare our results with optical reddening studies and so we discuss the assumptions and uncertainties of $\kappa_v$ in Section 4.4.

Fig. 5 shows the results. The two solid lines show the surface density of dust we derive when starting from the best-fitting models to the 250 and 500 µm surface-brightness profiles. The dashed lines show fits of single-exponential models to these empirical dust column-density distributions.

## 3.4 Radial variation of other galaxy phases

This work requires comparing our dust profile to the other constituents which make up galaxies. For the molecular gas and total gas we used the results available in the literature (Schruba et al. 2011; Bigiel & Blitz 2012), derived from the HERACLES (Leroy et al. 2009) and THINGS (Walter et al. 2008) samples. This could introduce a potential bias due to the different selection criteria used for HERACLES and THINGS relative to the HRS, but equivalent data are not available for the HRS sample.

To derive the radial distributions of stellar surface density and star-formation rate (SFR) surface density, we used two alternative data sets. *Spitzer* and *GALEX* data exist for many of the galaxies (see Table 1), and we used this data to determine empirical average radial distributions of stellar surface density and SFR surface

density to compare with the radial distributions of the other galaxy components. However, since 24 µm *Spitzer* data does not exist for the whole sample and a few of the 3.6 µm images are smaller than ideal, we also estimated profiles using WISE data, which does exist for the entire sample. In the case of the WISE data, we first convolved all the images to the same resolution as the 250 µm image to make the comparison with the dust surface density as fair as possible. We did not do this when using the *Spitzer* data as the higher resolution profile provides an independent check of the model fitting applied to the *Herschel* data. Table 1 summarizes the data sets. We obtained very similar results whichever data set we started with.

To find the radial variation of the surface density of the stars, we assumed the 3.6 µm emission is proportional to the column density of the stars. We used *Spitzer* 3.6 µm observations which are available for 44 of the 45 in our sample of large galaxies and were obtained as part of the S\*G survey (Sheth et al. 2010). Our second data set was the WISE 3.4 µm data which exists for all 45 galaxies in the sample. In the case of this data set, we first convolved the images to the same resolution as the 250 µm images before estimating the average radial distribution of stellar surface density.

To find the radial variation of the SFR surface density, we first created a SFR map for each galaxy by combining an obscured star-formation tracer (either 22 or 24 µm emission) and an un-obscured star formation tracer, the ultraviolet emission, using
<table><caption>**Table 1.** Exponential scalefactors.</caption><thead><tr><th>Component</th><th>Data source</th><th>Number of objects</th><th>Matched processing</th><th>Gradient (dex R<sub>25</sub><sup>-1</sup>)</th></tr></thead><tbody><tr><td>Dust</td><td>(250 µm and T)</td><td>45</td><td>Y</td><td>-1.72 ± 0.01</td></tr><tr><td>Dust</td><td>(500 µm and T)</td><td>45</td><td>Y</td><td>-1.75 ± 0.02</td></tr><tr><td>Stellar mass</td><td>WISE, 3.4 µm</td><td>45</td><td>Y</td><td>-1.78 ± 0.07</td></tr><tr><td>Stellar mass</td><td>S<sup>4</sup>G, 3.6 µm</td><td>44</td><td>N</td><td>-1.90 ± 0.03</td></tr><tr><td>SFR</td><td>WISE 22 µm and NUV</td><td>44</td><td>Y</td><td>-2.69 ± 0.05</td></tr><tr><td>SFR</td><td>MIPS 24 µm and NUV</td><td>28</td><td>N</td><td>-2.58 ± 0.05</td></tr><tr><td>Total Gas<sup>a</sup></td><td>H<sub>1</sub> and H<sub>2</sub></td><td>33</td><td>—</td><td>-0.71</td></tr><tr><td>H<sub>2</sub><sup>a</sup></td><td>—</td><td>33</td><td>—</td><td>-2.17</td></tr><tr><td>Metallicity<sup>a</sup></td><td>—</td><td>306</td><td>—</td><td>-0.16 (with σ = 0.12)</td></tr></tbody></table>

Notes. the gradients from exponential fits to data of different components of a galaxy. The 'Matched Processing' column lists whether the data have been convolved to the 250 µm resolution (or already have similar resolution) and have had the same procedure applied as for the dust profiles, or if the model has been fit to high-resolution data. Components marked <sup>a</sup> are gradients taken from the literature (Schruba et al. 2011; Bigiel & Blitz 2012; Sánchez et al. 2014) for different samples of galaxies.

published recipes (Leroy et al. 2008). Ideally we would have used the far-UV (UV) data from GALEX, but many of the observations are too shallow for this analysis. Instead, we used GALEX near-UV (NUV) images collected for the HRS (Cortese et al. 2012b) that exist for all but one of our sample, and assumed a conversion factor of 0.657 to convert NUV fluxes to FUV fluxes (the conversion factor was calculated from the global UV fluxes for the sample of large galaxies). To estimate the obscured component, we used the *Spitzer* 24 µm data (Bendo, Galliano & Madden 2012), which exists for about half the sample, and WISE 22 µm data, which exists for the whole sample. In the case of the WISE data, we first convolved the WISE and NUV images to the same resolution as the 250 µm images before producing the SFR maps for each galaxy. We also tested the effect of applying a small correction, using the 3.6 µm emission, to account for the stellar contribution to both the infrared and ultraviolet data (Leroy et al. 2008). This correction was found to have a negligible effect on our conclusions. To summarize, we were able to create high-resolution SFR maps (resolution of ~5 arcsec) for 28 of our 45 galaxies and maps at similar resolution to the SPIRE 250 µm images for 44 out of 45 galaxies.

Before producing the average distributions of either stellar surface density or SFR surface density, we masked any bright objects on the images, which is especially important due to contamination from stars in these wavebands, and also removed a linear background from the images. In creating the distributions we used very similar methods to the one used to create the radial distributions of submillimetre emission. The smaller sizes of the *Spitzer* 3.6 µm and 24 µm images meant that our default background annulus was reduced to 2.0–2.5 $R_{25}$. However, for about half the galaxies the 3.6 µm images were still too small for the background to be determined from this annulus. For these galaxies, we instead estimated the background on each image from an annulus with $1.0 < R/R_{25} < 1.5$ and then added a corrective offset to the profile calculated in the same region from the average profile for the other galaxies.

# 4 RESULTS AND DISCUSSION

## 4.1 Dust radial profile

The surface-brightness profile produced from combining the radial profiles from the 117 images (see Section 3.1) obtained with the Herschel SPIRE camera is shown in Fig. 1. Emission from the dust is detected at $>3\sigma$ to $R = 1.975 R_{25}$ at 250 µm and at $>3\sigma$ to $R = 1.75 R_{25}$ at 350 and 500 µm (the $3\sigma$ sensitivity limit is shown

by the red-dashed line in the figure). While an increase to $R \sim 2.0 R_{25}$ may not sound a significant improvement over studies that traced dust to $R \sim 1.2 R_{25}$, this corresponds to over a factor of 10 increase in sensitivity.

The average surface-brightness profile for the sample of large galaxies and the best-fitting models to the profile are shown in Fig. 2. We find that at 250 µm the average surface-brightness profile is best-fitted with a broken exponential, but that at 350 and 500 µm a single exponential is sufficient (for details of the models see Section 3.3).

The continuum emission from dust depends both on the column density and temperature of the dust, and so to determine the column density, we determined the average line-of-sight temperature of the dust using the method described in Section 3.3. On the assumption that the dust is heated by starlight, both the light from the old stellar population and the hot young stars in star-formation regions, we have used our stellar surface density, SFR surface density and dust-mass profiles to predict the dust temperature. We used the STELLAR KINEMATICS INCLUDING RADIATIVE TRANSFER (SKIRT) code (Camps & Baes 2015), which models the radiative transfer of starlight through dust using a Monte Carlo based approach. The results of the temperature prediction are shown in Fig. 4.

The best-fitting temperature (Fig. 4) declines linearly from $\approx 24$ K in the centre to approximately 15 K at a radius of $\approx R_{25}$ where the temperature appears to reach a minimum. Our radiative transfer model of the dust in a galaxy predicts a similar decline in the dust temperature with radius but does not predict this plateau in the dust temperature (Fig. 4). We note that fairly high dust temperatures well outside discs have been seen before, with Bocchio et al. (2016) finding that the dust around the edge-on galaxy NGC 891 still has a temperature of $\sim 19$ K 4 kpc above the plane.

Kashiwagi & Suto (2015) estimated that the cosmic UV background (Xu et al. 2005) would produce a minimum dust temperature of $\sim 10$ K. However, given dust radiates with a $T^6$ dependence, the UV background would have to be higher by $\sim 1.5^6$ to explain our minimum temperature of 15 K. Another possible heating source is from collisional heating in a hot X-ray gas. Bocchio (2014) found dust that had been pulled out of NGC 4338 was heated to 20–30 K in the hot X-ray gas in the Virgo Cluster, and Yamada & Kitayama (2005) show that dust could reach $\sim 30$ K in the centres of clusters. However, although spiral galaxies are known to have haloes of hot gas (e.g. Anderson, Churazov & Bregman 2016) it is unclear whether this gas could heat the dust to $\sim 15$ K.

We used our temperature estimates and the models that fit the surface-brightness profiles at the SPIRE wavelengths to determine
<figure><img src="image_10.png" /><figcaption>**Figure 6.** Average surface density versus radius for all components of a galaxy. The data points show the estimated average surface density for the HRS galaxies in each radial bin for dust (red) star-formation rate (green) and stellar mass (blue). The stellar mass and SFR data points are radial profiles measured from high-resolution data sets (≤6.0 arcsec full-width half maximum) from *Spitzer* and *GALEX*. The dashed lines of the same colour show the best-fitting exponential models for images which have been convolved to the same resolution as our 250 µm images and have had the same procedure applied to produce the average profiles. In this case we have used WISE data instead of *Spitzer* data as observations exist for the entire sample and cover a larger area. In the case of the total gas (black) and molecular gas (magenta) we have shown the exponential distributions derived from samples of local galaxies in the literature (Section 3.4).</figcaption></figure>

the radial variation in the column density of the dust (Fig. 5). We find that the dust column density has an exponential decline with radius (gradient ~-1.7 dex $R_{25}^{-1}$), no matter which model of the dust emission we start with. Throughout this paper the gradients ($\alpha$) are given in units of dex $R_{25}^{-1}$ (i.e. a value of -1 corresponds to a decrease by a factor of 10 over a radial distance of $R_{25}$, and can be described by:

$$ \Sigma = \Sigma_0 10^{\alpha r}, \qquad (1) $$

where $\Sigma$ is the surface density, $\Sigma_0$ is the surface density when $r=0$, and $r$ is the radius in angular units. If we extrapolate our exponential surface density model to an infinite radius, we estimate that the dust between $1.0 R_{25} < R \le \infty$ constitutes only 9.3 per cent of all the total dust in the disc

## 4.2 Comparison with other galaxy components

Fig. 6 shows all galaxy components plotted on one diagram and Table 1 gives the gradients of the exponential models that best fit the radial distributions of stellar surface density, SFR surface density, dust surface density, molecular gas and total gas. In the case of stellar surface density and SFR surface density, we obtained very similar values for the gradient whether we started with the high-resolution images or the images convolved to the same resolution as the 250 µm images. The blue and green points in Fig. 6 show the empirical average radial distributions of stellar surface density and SFR surface density estimated from the high-resolution data sets. The blue and green dashed lines show the fits of single exponential models to the empirical radial distributions of stellar surface density

and SFR surface density derived from the low-resolution data sets. In this case, the empirical distributions are not shown in Fig. 6.

As expected, the average surface density of all galaxy components falls with distance from the centre of the galaxy. As expected, given the extended distribution of atomic gas around galaxies (Leroy et al. 2008), the total column density of gas declines much more slowly than the other components. The gradients of the distribution of stellar mass and of dust are in good agreement, but the SFR falls more rapidly with radius than either of the other two.

Hunt et al. (2015) also find similar distributions of the dust and stars, finding the scalelengths of KINGFISH spirals to have a value of $R_{250}/R_{25} \approx 0.37$ at 250 µm versus $R_{3.6}/R_{25} \approx 0.35$ at 3.6 µm ($1.17 \text{ dex } R_{25}^{-1}$ and $1.24 \text{ dex } R_{25}^{-1}$, respectively, in our gradient units). Our results show that the similarity between the distributions of dust and stellar mass extends out to $R \sim 1.8 R_{25}$. Our dust gradient of $-1.75 \text{ dex } R_{25}^{-1}$ is steeper than the value of $1.17 \text{ dex } R_{25}^{-1}$ for the KINGFISH spiral galaxies (Hunt et al. 2015), although the KINGFISH dust distributions are distributions of dust 250 µm emission rather than dust surface density. Muñoz-Mateos et al. (2009) find the median gradient of the radial distribution of the surface density of dust for the entire SINGS survey of $1.49 \text{ dex } R_{25}^{-1}$, significantly closer to our value.

It is still uncertain whether dust is mostly formed in evolved stars (Matsuura et al. 2009; Boyer et al. 2012), supernovae (Gomez et al. 2012; Indebetouw et al. 2014) or formed by grain growth in the interstellar medium (ISM; Ossenkopf 1993; Hirashita & Kuo 2011). We note that if dust grains are formed in supernovae and then fairly quickly destroyed, we would expect the surface densities of the dust and the SFR to have the same dependence on radius, which is not found.
## 4.3 Dust, gas and metallicity

Metallicity gradients in galaxies tend to be small. The average metallicity gradient for a sample of 306 galaxies measured with CALIFA is $-0.16 \text{ dex } R_{25}^{-1}$ (scatter of $\sim 0.12 \text{ dex } R_{25}^{-1}$; Sánchez et al. 2014), from observations of H II regions out to $1.3 R_{25}$ (there is possible evidence of a flattening at larger radii). For the galaxies in the SINGS survey the measured radial metallicity gradients are steeper, $\sim -0.42$ to $-0.33 \text{ dex } R_{25}^{-1}$ (depending on the calibration used; Moustakas et al. 2010). On the assumption that these metallicity measurements apply to all phases of the interstellar medium, we combine the gradient in metallicity from the CALIFA measurement and the gradient in the column density of all gas (Table 1) and find that the gradient in the column density of metals is $\sim -0.9 \text{ dex } R_{25}^{-1}$, again significantly flatter than the dust gradient of $-1.7 \text{ dex } R_{25}^{-1}$. Similar radial declines in the dust-to-metal ratio have been seen in M31 (Mattsson et al. 2014) and in SINGS galaxies (Mattsson & Andersen 2012). In closed-box chemical evolution models this radial decline can be explained by dust grain growth in the ISM (Mattsson, Andersen & Munkhammar 2012).

## 4.4 Are galaxies surrounded by extended haloes of dust?

Based on the angular correlation between the reddening of $\sim 85 000$ high-redshift quasars and the positions of $\sim 2.4 \times 10^7$ galaxies detected in the SDSS (York et al. 2000), Ménard et al. (2010a, henceforth M10) concluded that galaxies are surrounded by very extended haloes of dust. They found that the mean dust surface density varies with angular distance, $\theta$, from the position of an SDSS galaxy as

$$ \Sigma_d(\theta) = 2.5 \times 10^{-3} \left( \frac{\theta}{0.1 \text{ arcmin}} \right)^{-0.8} h (M_{\odot} \text{ pc}^{-2}), \quad (2) $$

(taken from fig. 8 of M10), with this distribution extending to a physical distance of several Mpc ($h$ is $H_0/100$). Kashiwagi & Suto (2015) have suggested that the reddening measurements are actually due, not to extended haloes of dust around individual galaxies, to the dust within galactic discs and the fact that the positions of galaxies are correlated. Their evidence for this claim are the results from stacking measurements around the positions of the same SDSS galaxies, using the far-infrared image of the sky made with IRAS. They, too, detect extended dust around the SDSS galaxies, but this time in emission, and by comparing the reddening and emission measurements, they estimate the temperature of the dust is $\sim 20$ K, which they argue is what one would expect for dust in the galactic discs rather than dust well outside galactic discs (however, see Section 4.1).

The strong practical motive for trying to distinguish between these two possibilities is the need to understand the effect of dust on experiments to measure cosmological parameters, such as experiments that use the magnitudes of Type Ia Supernovae to measure $\Omega_\Lambda$ and the equation of state of dark energy (Ménard, Kilbinger & Scranton 2010b). If all galaxies are surrounded by such extended haloes of dust, all measurements of the magnitudes of high-redshift objects need to be corrected for the effects of dust extinction and reddening. But, if the conclusion of Kashiwagi & Suto (2015) is correct, and almost all dust is confined to galaxy discs, only a small percentage of lines of sight actually pass through galactic discs (we quantify this below). We now discuss what our results can contribute to resolving the controversy. We borrow the nomenclature of Kashiwagi & Suto (2015) in which the first model is referred to

as the circumgalactic dust model (CGD model) and the latter model as the interstellar dust model (ISD model).

Our result that the dust column density of galaxies declines more rapidly with radius out to $R \sim 2 R_{25}$ than given by equation (2) (Fig. 5) does not allow us to discriminate between the ISD and the CGD models. This is because equation (2) is only valid beyond the minimum radius for which M10 have a reddening measurement which is only slightly less than $2 R_{25}$ (Fig. 5), where we are not sensitive. Therefore it is possible that the surface density of dust declines rapidly with radius out to $R \sim 2 R_{25}$ but then follows the more gentle relationship given by equation (2).

We now adopt the alternative approach of combining the statistical properties of the SDSS galaxies – their number counts and clustering properties – with the average dust surface density radial profile for the HRS galaxies (Fig. 5) to test whether we can produce a relationship similar to equation (2) from the dust within galaxy discs.

M10 used SDSS galaxies with $i < 21$, whereas our study is ultimately based on the photometric B-band, since $R_{25}$ is the distance from the centre of a galaxy out to an isophote with a surface brightness in the B-band of $25 \text{ mag arcsec}^{-2}$. To use a sample of galaxies as similar as possible to that used by M10 and to connect the sample to the standard definition of optical radius we used the galaxy number counts derived from the SDSS r-band images and an average B-r colour of 1.65 (Yasuda et al. 2001). Since the mean r-i colour measured by Yasuda et al. (2001) is $\sim 0.4$, we used the number counts with $r < 21.5$.

Of course, not all galaxies detected in the SDSS contain much dust, and in particular early-type galaxies contain small amounts of dust (Smith et al. 2012a). Eales et al. (2015) have recently measured the percentage of the stellar mass in the nearby Universe that is in late-type galaxies is $\sim 49$ per cent, using a definition of ‘late-type’ that means a late-type galaxy is almost certainly dominated by a disc. Informed by this result, we multiply the SDSS number counts by 0.5 to produce an estimate of the number counts of SDSS galaxies which contain a disc, $N_{\text{disc}}$. We then make the assumption that each of these galaxies has the dust surface density radial profile shown in Fig. 5. Although these are large assumptions, they are not as important as they appear because of the larger uncertainty in the amount of evolution of dust in galaxies, which we discuss below.

On the assumption that the surface brightness of a galaxy depends on radius with the form shown in equation (1), the total flux density out to the standard optical radius is given by

$$ I_{\text{tot}} = \frac{2\pi I_{25}}{(\alpha \log_e(10))^2} \left( (\alpha R_{25} \log_e(10) - 1) + 10^{-\alpha R_{25}} \right), \quad (3) $$

in which $I_{25}$ is the surface brightness corresponding to $25 \text{ mag arcsec}^{-2}$. Because $\alpha$ is given in units of $R_{25}$ (Table 1), this equation makes it possible to calculate $\theta_{25}$, and thus the area subtended on the sky, for a galaxy of any apparent magnitude. Pohlen & Trujillo (2006) have fit exponential models to the SDSS images of 90 face-on late-type galaxies. The mean value of $\alpha$ derived from the values given by them for the inner scalelength in the g-band for each galaxy (their table 3) is $-1.64 \pm 0.08 \text{ dex } R_{25}^{-1}$. Given the average value of $\alpha$ we can use equation (3) to calculate the optical radius of a galaxy of any apparent magnitude.

A very simple thing we can then do is calculate the fraction of the sky covered by the discs. As we can only detect dust out to a radius of $2 \times R_{25}$, we assume that a disc only extends out to this radius. Using this assumption, the number counts for late-type galaxies, $N_{\text{disc}}$, and the average value of $\alpha$ from Pohlen & Trujillo (2006), we estimate that the percentage of the sky covered by discs
<figure><img src="image_11.png" /><figcaption>**Figure 7.** The relationship between mean dust surface density and angular distance from the centre of an SDSS galaxy. The dashed black line shows the relationship estimated by M10 from looking at the correlation between the reddening of high-redshift quasars and the positions of SDSS galaxies. The solid red line shows our prediction for this relationship considering only the dust in galactic discs, with the shaded region showing the predictions if we use values for α in our model between the upper and lower 1σ errors in the mean value of α (see the text for details). The blue line shows our prediction if we incorporate the evolution in the dust masses of galaxies measured by Bourne et al. (2012).</figcaption></figure>

is 0.44 ± 0.05 per cent, with the errors obtained using the error for the mean value of α given above. This low percentage confirms the importance of distinguishing between the CGD and ISD models for making the correct corrections for dust extinction and reddening when carrying out cosmological experiments, because in the ISD model a very small fraction of the sky is covered by dust.

For a galaxy in which the dust surface density varies with radius with the relationship given in equation (1) the mean column density out to twice the optical radius ($2 \times \theta_{25}$) is given by:

$$ \langle \Sigma_{\text{d,galaxy}} \rangle = \frac{\Sigma_0 (10^{2\alpha} R_{25} (2\alpha R_{25} \log_e(10) - 1) + 1)}{2 (\alpha R_{25} \log_e(10))^2}. \quad (4) $$

We use the values for $\Sigma_0$ and $\alpha$ for the mean dust surface-density profile shown in Fig. 5: $\Sigma_0 = 0.6 M_{\odot} pc^{-2}$ and $\alpha = -1.75 \text{ dex } R_{25}^{-1}$. The expected average dust column density profile around SDSS galaxies is then given by:

$$ \Sigma_{\text{d}}(\theta) = \int_{0}^{21} \langle \Sigma_{\text{d,galaxy}} \rangle 4\pi (\theta_{25})^2 w(\theta, m) N_{\text{disc}}(m) dm. \quad (5) $$

In this equation, $w(\theta, m)$ is the angular correlation function for SDSS galaxies. Connolly et al. (2002) have measured the angular correlation function for SDSS galaxies as a function of r magnitude. They found that at all magnitudes the angular correlation function is well represented by $w(\theta) = A\theta^{-\gamma}$, with $\gamma \simeq 0.7$ and $\log_{10}A$ a linear function of r-band magnitude. In calculating $w(\theta, m)$, we assume the same value of $\gamma$ and estimate A at each magnitude from a linear fit to the data shown in fig. 3 of Connolly et al. (2002).

Fig. 7 shows our estimate of $\Sigma_{\text{d}}(\theta)$ from this procedure. The red solid line in the figure shows our prediction with the mean value of $\alpha$ we calculated from the results in Pohlen & Trujillo (2006)

and the shaded region shows the predictions if we use values for $\alpha$ between the $\pm 1 \sigma$ errors on the mean. The black dashed line shows the relationship measured by M10 (equation 2). This relationship is very similar to the result from the reddening study (M10).

This is without considering the effect of galactic evolution, as the mean redshift of the SDSS galaxies is $\simeq 0.3$. Two results, both from the Herschel-ATLAS survey (Eales et al. 2010), have shown that there is substantial evolution in the amount of dust in galaxies by this redshift. Dunne et al. (2011) have shown that the dust masses of galaxies have increased by a factor of $\sim 4.6$ by this redshift. This result is based on the small fraction of galaxies with very high dust masses, but Bourne et al. (2012) have carried out a stacking analysis on a large sample of SDSS galaxies which is sensitive to galaxies with much lower dust masses, finding that the average dust mass of optically selected blue galaxies of the same stellar mass has increased by a factor of $\simeq 3$ by a redshift of 0.3. In Fig. 7 we have shown our estimate of $\Sigma_{\text{d}}(\theta)$ scaled by a factor of 3.0 to allow for the effect of galaxy evolution. This now exceeds the relationship found in the reddening study.

An important consideration is the absolute and relative uncertainties in the values of $\kappa_{\nu}$ used for both our analysis of dust emission and the reddening study of M10. Our value of $\kappa_{\nu} = 0.192 \text{ m}^2 \text{ kg}^{-1}$ at 350 µm (Draine 2003) assumes Milky Way like carbonaceous and silicate dust, and is in the typical range of values used of $\kappa_{\nu}$ in local galaxies studies (see Clark et al. 2016, and references therein). Our $\kappa_{\nu}$ value is higher by a factor of $\sim 2$ than the latest empirical value derived in this study (Clark et al. 2016); using this value would roughly double the amount of dust we measure in the discs of galaxies. As we are only measuring the dust emission in the discs of galaxies these values should be reasonable values to assume out
to 2.0 $R_{25}$, particularly as the metallicity gradients of galaxies tends to be small. A greater unknown is what the extinction coefficient would be for dust in the IGM; M10 assume an SMC-type dust due to the lack of a 0.2 µm bump in low-ionization absorbers in the haloes of ~L* galaxies (e.g. York et al. 2006) and the fact that extinction curves for high-redshift galaxies also do not show the 0.2 µm bump. If larger Milky Way-type dust were assumed, the halo dust surface density implied by the M10 results would approximately double. Of course, the dust in intergalactic space might be very unlike the dust within galactic discs. Sputtering is likely to be important (Yamada & Kitayama 2005) and possibly only large grains would survive being expelled from galactic discs (Croft et al. 2000). We emphasize, however, that in our model 90 per cent of the dust contributing to the average surface density shown in Fig. 7 is produced by dust within a radius less than $R_{25}$, where our knowledge of the properties of dust is much better.

Of course there are also large uncertainties in our model, including the division between early-type and late-types, the assumption that all late-type galaxies have the same distribution of dust, and the amount of galactic evolution. Because of these uncertainties and because of the huge uncertainty as to whether the dust in intergalactic space is like dust in galactic discs, we cannot rule out the CGD model. All we can say is that with some assumptions about the properties of dust (e.g. Milky Way-type dust grains), the results of M10 might be explained by the dust within individual galaxies rather than very extended haloes of dust around all galaxies.

Despite using results from most sensitive ever far-infrared space telescope, and then combining the results in a powerful statistical stacking analysis, we have still only managed to detect the continuum emission from dust out to a distance of twice $R_{25}$ from the galactic centre. Therefore the only practical way to investigate the properties of dust at greater distances from centres of galaxies is still through statistical investigations of the reddening of high-redshift objects. Over three decades, investigations of the dust associated with damped Ly$\alpha$ quasar absorption-line systems have produced many significant detections of this dust and have even produced detailed measurements of the properties of this dust (e.g. Fall, Pei & McMahon 1989; Khare et al. 2004; York et al. 2006; Wild, Hewett & Pettini 2007; Vladilo, Prochaska & Wolfe 2008). There have also been detections through quasar absorption-line of gas and metals in the haloes of intervening individual galaxies, with projected distances of up to 100 kpc (e.g. Steidel et al. 2002; Chen, Kennicutt & Rauch 2005), although some galaxies fail to show any absorption (Bechtold & Ellingson 1992). Investigations of the correlation of quasar absorption-line systems with galaxy redshift surveys imply that the absorption-line gas is often a long way from the centre of a galaxy (Churchill, Steidel & Kacprzak 2005 and references therein; Pérez-Rafols et al. 2015), which suggests that much of the dust detected via reddening of quasars is also outside galactic discs. The method used by M10 is an elegant way of both measuring the properties of the dust and also its distribution around galaxies. We have shown that there could be a significant contribution to the reddening measured by M10 from dust within galactic discs.

Fortunately, there is actually a simple way of adapting the M10 technique to allow for this contamination. Even if the CGD and ISD models predict very similar relationships for how the mean reddening depends on distance from an SDSS galaxy, in the ISD case this relationship is caused by the small percentage of lines of sight that pass through galactic discs, whereas in the CGD case it is caused by all lines of sight. Therefore, a simple way to distinguish between the two hypotheses would be to measure the variance of the quasar reddening measurements, which should be much higher

for the ISD model than the CGD model – a test that should be easy to carry out using the measurements described by M10.

## 5 CONCLUSIONS

In this paper we present the results of combining images of a large sample of nearby galaxies made from observations with the Herschel Space Observatory. By combining the images, we determine how the average surface density of dust depends on radius. We find the following results:

(i) We detect the emission of dust at $>3\sigma$ threshold out to a radius of 2.0 $R_{25}$ at 250 µm and out to a radius of 1.75 $R_{25}$ at 350 and 500 µm, by combining images of 117 galaxies from the HRS. When we restrict the sample to the 45 galaxies with $R_{25} > 1.5$ arcmin we can trace dust at 250 µm out to a radius of 1.75 $R_{25}$ at $> 3\sigma$.

(ii) We fit modified blackbody spectral-energy distributions to the SPIRE and PACS (100–500 µm) radial profiles. We find the dust temperature declines from ~25 K in the centre of the galaxy to ~15 K at 1.0 $R_{25}$ where it remains constant out to 1.75 $R_{25}$.

(iii) We find the radial dust-mass distribution in a galaxy is consistent with an exponential profile out to a radius of 1.8 $R_{25}$. The gradient of the surface-density profile is $-1.7 \text{ dex } R_{25}^{-1}$.

(iv) We created profiles of the stellar mass (traced by 3.4 or 3.6 µm) and SFR (UV and Infrared) using images from Spitzer, GALEX and WISE. We found the stellar gradient is in good agreement with the dust gradient with the dust gradient between 1 per cent and 10 per cent larger (in dex $R_{25}^{-1}$) depending on the data set used. The SFR profile is significantly steeper than that of the dust with a gradient of ~$-2.6 \text{ dex } R_{25}^{-1}$.

(v) Our results show a much steeper gradient in the surface density of dust compared to metals (based on typical metallicity gradients for nearby galaxies), similar to the results for individual galaxies (e.g. Mattsson & Andersen 2012). In closed-box chemical evolution models this is a signature that dust originates from grain growth in the ISM.

(vi) We use our results for the radial distribution of dust in galaxies to see if the large extended haloes of dust around galaxies as described by M10 could have a significant contamination caused by the clustering of galaxies and the dust emission within individual galaxies. By using the number counts and angular correlation function from SDSS, we show that, assuming Milky Way-type dust, the dust within galaxy discs might explain the reddening results. Because of the large uncertainties in our analysis and the even larger uncertainty of whether the dust in intergalactic space is like dust in galactic discs, we cannot rule out the existence of intergalactic dust. We propose a simple method of testing whether there is actually intergalactic dust.

Our results from a statistical stacking analysis with images from the Herschel Space Observatory represent the most sensitive study of the extended dust emission around galaxies that will be possible until a new submillimetre space mission. It is likely to be impossible to improve on our results with ground-based telescopes because of the fluctuating emission from the atmosphere, which makes it difficult to detect low surface-brightness emission from galaxies.

## ACKNOWLEDGEMENTS

We thank everyone involved with the Herschel Space Observatory.

SPIRE has been developed by a consortium of institutes led by Cardiff University (UK) and including: University of Lethbridge (Canada); NAOC (China); CEA, LAM (France); IFSI, University
of Padua (Italy); IAC (Spain); Stockholm Observatory (Sweden); Imperial College London, RAL, UCL-MSSL, UKATC, University of Sussex (UK); and Caltech, JPL, NHSC, University of Colorado (USA). This development has been supported by national funding agencies: CSA (Canada); NAOC (China); CEA, CNES, CNRS (France); ASI (Italy); MCINN (Spain); SNSB (Sweden); STFC, UKSA (UK); and NASA (USA).

HIPE is a joint development by the Herschel Science Ground Segment Consortium, consisting of ESA, the NASA Herschel Science Center and the HIFI, PACS, and SPIRE consortia.

This research made use of Astropy, a community-developed core `PYTHON` package for Astronomy (Astropy Collaboration, 2013).

HLG acknowledges support from the European Research Council (ERC) in the form of Consolidator Grant COSMICDUST (ERC-2014-CoG-647939).

## REFERENCES

* Alton P. B. et al., 1998, A&A, 335, 807

* Anderson M. E., Churazov E., Bregman J. N., 2016, MNRAS, 445, 227

* Aniano G., Draine B. T., Gordon K. D., Sandstrom K., 2011, PASP, 123, 1218

* Bechtold J., Ellingson E., 1992, ApJ, 396, 20

* Bendo G. J., Galliano F., Madden S. C., 2012, MNRAS, 423, 197

* Bigiel F., Blitz L., 2012, ApJ, 756, 183

* Bocchio M., 2014, PhD thesis, Institut d'Astrophysique Spatiale (IAS), UMR 8617, CNRS/Université Paris-Sud, 91405 Orsay, France

* Bocchio M., Bianchi S., Hunt L. K., Schneider R., 2016, A&A, 586, A8

* Boselli A. et al., 2010, PASP, 122, 261

* Bourne N. et al., 2012, MNRAS, 421, 3027

* Boyer M. L. et al., 2012, ApJ, 748, 40

* Camps P., Baes M., 2015, Astron. Comput., 9, 20

* Chen H.-W. Kennicutt R. C., Jr, Rauch M., 2005, ApJ, 620, 703

* Churchill C., Steidel C., Kacprzak G., 2005, in Braun R., ed., ASP Conf. Ser. Vol. 331, Extra-Planar Gas. Astron. Soc. Pac., San Francisco, p. 387

* Ciesla L. et al., 2012, A&A, 543, A161

* Clark C. J. R., Schofield S. P., Gomez H. L., Davies J. I., 2016, MNRAS, 459, 1646

* Connolly A. J. et al., 2002, ApJ, 579, 42

* Cortese L. et al., 2010, A&A, 518, L49

* Cortese L. et al., 2012a, A&A, 540, A52

* Cortese L. et al., 2012b, A&A, 544, A101

* Cortese L. et al., 2014, MNRAS, 440, 942

* Croft R. A. C., Davé R., Hernquist L., Katz N., 2000, ApJ, 534, L123

* Davies J. I. et al., 2010, A&A, 518, L48

* Davies J. I. et al., 2012, MNRAS, 419, 3505

* Draine B. T., 2003, ARA&A, 41, 241

* Dunne L. et al., 2011, MNRAS, 417, 1510

* Eales S. et al., 2010, PASP, 122, 499

* Eales S. et al., 2012, ApJ, 761, 168

* Eales S. et al., 2015, MNRAS, 452, 3489

* Fall S. M., Pei Y. C., McMahon R. G., 1989, ApJ, 341, L5

* Foyle K. et al., 2012, MNRAS, 421, 2917

* Gavazzi G., Boselli A., Donati A., Franzetti P., Scodeggio M., 2003, A&A, 400, 451

* Gomez H. L. et al., 2012, ApJ, 760, 96

* Griffin M. J. et al., 2010, A&A, 518, L3

* Hirashita H., Kuo T.-M., 2011, MNRAS, 416, 1340

* Hodges-Kluck E., Bregman J. N., 2014, ApJ, 789, 131

* Holwerda B. W., Keel W. C., Williams B., Dalcanton J. J., de Jong R. S., 2009, AJ, 137, 3000

* Hughes T. M. et al., 2014, A&A, 565, A4

* Hunt L. K. et al., 2015, A&A, 576, A33

* Indebetouw R. et al., 2014, ApJ, 782, L2

* Kashiwagi T., Suto Y., 2015, MNRAS, 451, 4162

* Khare P., Kulkarni V. P., Lauroesch J. T., York D. G., Crotts A. P. S., Nakamura O., 2004, ApJ, 616, 86

* Kirkpatrick A. et al., 2014, ApJ, 789, 130

* Leroy A. K., Walter F., Brinks E., Bigiel F., de Blok W. J. G., Madore B., Thornley M. D., 2008, AJ, 136, 2782

* Leroy A. K. et al., 2009, AJ, 137, 4670

* Matsuura M. et al., 2009, MNRAS, 396, 918

* Mattsson L., Andersen A. C., 2012, MNRAS, 423, 38

* Mattsson L., Andersen A. C., Munkhammar J. D., 2012, MNRAS, 423, 26

* Mattsson L. et al., 2014, MNRAS, 444, 797

* Ménard B., Scranton R., Fukugita M., Richards G., 2010a, MNRAS, 405, 1025 (M10)

* Ménard B., Kilbinger M., Scranton R., 2010b, MNRAS, 406, 1815

* Mentuch Cooper E. et al., 2012, ApJ, 755, 165

* Moustakas J., Kennicutt R. C., Jr, Tremonti C. A., Dale D. A., Smith J.-D. T., Calzetti D., 2010, ApJS, 190, 233

* Muñoz-Mateos J. C. et al., 2009, ApJ, 701, 1965

* Nelson A. E., Zaritsky D., Cutri R. M., 1998, AJ, 115, 2273

* Nguyen H. T. et al., 2010, A&A, 518, L5

* Ossenkopf V., 1993, A&A, 280, 617

* Peek J. E. G., Ménard B., Corrales L., 2015, ApJ, 813, 7

* Pérez-Rafols I., Miralda-Escudé J., Lundgren B., Ge J., Petitjean P., Schneider D. P., York D. G., Weaver B. A., 2015, MNRAS, 447, 2784

* Pilbratt G. L. et al., 2010, A&A, 518, L1

* Poglitsch A. et al., 2010, A&A, 518, L2

* Pohlen M., Trujillo I., 2006, A&A, 454, 759

* Sánchez S. F. et al., 2014, A&A, 563, A49

* Schruba A. et al., 2011, AJ, 142, 37

* Sheth K. et al., 2010, PASP, 122, 1397

* Smith M. W. L. et al., 2010, A&A, 518, L51

* Smith M. W. L. et al., 2012a, ApJ, 748, 123

* Smith M. W. L. et al., 2012b, ApJ, 756, 40

* Steidel C. C., Kollmeier J. A., Shapley A. E., Churchill C. W., Dickinson M., Pettini M., 2002, ApJ, 570, 526

* Vladilo G., Prochaska J. X., Wolfe A. M., 2008, A&A, 478, 701

* Walter F., Brinks E., de Blok W. J. G., Bigiel F., Kennicutt R. C., Jr, Thornley M. D., Leroy A., 2008, AJ, 136, 2563

* Wild V., Hewett P. C., Pettini M., 2007, MNRAS, 374, 292

* Xu C. K. et al., 2005, ApJ, 619, L11

* Yamada K., Kitayama T., 2005, PASJ, 57, 611

* Yasuda N. et al., 2001, AJ, 122, 1104

* York D. G. et al., 2000, AJ, 120, 1579

* York D. G. et al., 2006, MNRAS, 367, 945

<footnote n="1">¹School of Physics and Astronomy, Cardiff University, The Parade, Cardiff CF24 3AA, UK</footnote>

<footnote n="2">²Sterrenkundig Observatorium, Universiteit Gent, Krijgslaan 281 S9, B-9000 Gent, Belgium</footnote>

<footnote n="3">³Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK</footnote>

<footnote n="4">⁴Jodrell Bank Centre for Astrophysics, School of Physics and Astronomy, Alan Turing Building, The University of Manchester, Oxford Road, Manchester M13 9PL, UK</footnote>

<footnote n="5">⁵INAF-Osservatorio Astrofisico di Arcetri, Largo Enrico Fermi 5, I-50125 Firenze, Italy</footnote>

<footnote n="6">⁶Unidad de Astronomía, Fac. de Ciencias Básicas, Universidad de Antofagasta, Avda. U. de Antofagasta 02800, Antofagasta, Chile</footnote>

<footnote n="7">⁷Laboratoire d'Astrophysique de Marseille-LAM, Université d'Aix-Marseille & CNRS, UMR7326, 38 Rue F. Joliot-Curie, F-13388 Marseille Cedex 13, France</footnote>

<footnote n="8">⁸Department of Physics, University of Crete, Heraklion GR-71003, Greece</footnote>

<footnote n="9">⁹Institute for Astronomy, Astrophysics, Space Applications and Remote Sensing, National Observatory of Athens, GR-15236 Penteli, Greece</footnote>

<footnote n="10">¹⁰INAF-Osservatorio Astronomico di Padova, Vicolo dell'Osservatorio 5, I-35122 Padova, Italy</footnote>

<footnote n="11">¹¹Astrophysics Group, Imperial College, Blackett Laboratory, Prince Consort Road, London SW7 2AZ, UK</footnote>
This paper has been typeset from a *TEX/LATEX* file prepared by the author.

<footnote n="12">¹²Center for Cosmology and the Department of Physics and Astronomy,
University of California, Irvine, CA 92697, USA</footnote>

<footnote n="13">¹³International Centre for Radio Astronomy Research, University of Western Australia, 35 Stirling Highway, Crawley WA 6009, Australia</footnote>

<footnote n="14">¹⁴Centro de Radioastronomía y Astrofísica, UNAM, Campus Morelia, A.P. 3-72, C.P. 58089, Mexico</footnote>

<footnote n="15">¹⁵Instituto de Física y Astronomía, Universidad de Valparaíso, Avda. Gran Bretaña 1111, Valparaíso, Chile</footnote>

<footnote n="16">¹⁶Astronomy Centre, Department of Physics and Astronomy, University of Sussex, Brighton BN1 9QH, UK</footnote>

<footnote n="17">¹⁷Infrared Processing and Analysis Center, California Institute of Technology, MS 100-22, Pasadena, CA 91125, USA</footnote>

<footnote n="18">¹⁸CEA, Laboratoire AIM, Irfu/SAp, Orme des Merisiers, F-91191 Gif-sur-Yvette, France</footnote>

<footnote n="19">¹⁹Istituto di Fisica dello Spazio Interplanetario, INAF, Via del Fosso del Cavaliere 100, I-00133 Roma, Italy</footnote>