The nature of the UV halo around the spiral galaxy NGC 3628

Thanks to deep UV observations with GALEX and Swift, diffuse UV haloes have recently been discovered around galaxies. Based on UV-optical colours, it has been advocated that the UV haloes around spiral galaxies are due to UV radiation emitted from the disc and scattered off dust grains at high latitudes. Detailed UV radiative transfer models that take into account scattering and absorption can explain the morphology of the UV haloes, and they require the presence of an additional thick dust disc next the to traditional thin disc for half of the galaxies in their sample. We test whether such an additional thick dust disc agrees with the observed infrared emission in NGC 3628, an edge-on galaxy with a clear signature of a thick dust disc. We extend the far-ultraviolet radiative transfer models to full-scale panchromatic models. Our model, which contains no fine-tuning, can almost perfectly reproduce the observed spectral energy distribution from UV to mm wavelengths. These results corroborate the interpretation of the extended UV emission in NGC 3628 as scattering off dust grains, and hence of the presence of a substantial amount of diffuse extra-planar dust. A significant caveat, however, is the geometrical simplicity and non-uniqueness of our model: other models with a different geometrical setting could lead to a similar spectral energy distribution. More detailed radiative transfer simulations that compare the model results to images from UV to submm wavelengths are a way to break this degeneracy, as are UV polarisation measurements.Comment: 6 pages, 2 figures, accepted for publication in Astronomy & Astrophysic

First Resolved Dust Continuum Measurements of Individual Giant Molecular Clouds in the Andromeda Galaxy

© 2020 The American Astronomical Society.In our local Galactic neighborhood, molecular clouds are best studied using a combination of dust measurements, to determine robust masses, sizes, and internal structures of the clouds, and molecular-line observations to determine cloud kinematics and chemistry. We present here the first results of a program designed to extend such studies to nearby galaxies beyond the Magellanic Clouds. Utilizing the wideband upgrade of the Submillimeter Array (SMA) at 230 GHz, we have obtained the first continuum detections of the thermal dust emission on sub-GMC scales (∼15 pc) within the Andromeda galaxy (M31). These include the first resolved continuum detections of dust emission from individual giant molecular clouds (GMCs) beyond the Magellanic Clouds. Utilizing a powerful capability of the SMA, we simultaneously recorded CO(2-1) emission with identical (u, v) coverage, astrometry, and calibration, enabling the first measurements of the CO conversion factor, α CO(2-1), toward individual GMCs across an external galaxy. Our direct measurement yields an average CO-to-dust mass conversion factor of α' CO-dust = 0.042 ± 0.018 M o (K km s -1 pc 2) -1 for the J = 2-1 transition. This value does not appear to vary with galactocentric radius. Assuming a constant gas-to-dust ratio of 136, the resulting α CO = 5.7 ± 2.4 M o (K km s -1 pc 2) -1 for the 2-1 transition is in excellent agreement with that of GMCs in the Milky Way, given the uncertainties. Finally, using the same analysis techniques, we compare our results with observations of the local Orion molecular clouds, placed at the distance of M31 and simulated to appear as they would if observed by the SMA.Peer reviewedFinal Published versio

An extremely low gas-to-dust ratio in the dust-lane lenticular galaxy NGC 5485

Date of Acceptance: 21/07/2014Peer reviewe

High-resolution radiative transfer modelling of M33

The authors thank the anonymous reviewer, whose comments have certainly improved the manuscript. The authors also thank Peter Camps for valuable technical advice, along with all of the participants at the recent SKIRT meeting for comments and discussions. I.D.L. gratefully acknowledges the supports of the Research Foundation – Flanders (FWO). M.W.L.S acknowledges funding from the UK Science and Technology Facilities Council consolidated grant ST/K000926/1. M.R. acknowledges support by the research projects AYA2014-53506-P and AYA2017-84897P from the Spanish Ministerio de Econom´ıa y Competitividad. This research made use of MONTAGE (http://montage.ipac.caltech.edu/), which is funded by the National Science Foundation under Grant Number ACI-1440620, and was previously funded by the National Aeronautics and Space Administration’s Earth Science Technology Office, Computation Technologies Project, under Cooperative Agreement Number NCC5-626 between NASA and the California Institute of Technology. This research has made use of Astropy, a community-developed core PYTHON package for Astronomy (http: //www.astropy.org/; Astropy Collaboration et al. 2013, 2018). This research has made use of NumPy (http://www.numpy.org/; van der Walt, Colbert & Varoquaux 2011), SciPy (http://www.scipy.org/), andMatPlotLib (http://matplotlib.org/; Hunter 2007). This research made use of APLpy, an open-source plotting package for PYTHON (https://aplpy.github.io/; Robitaille & Bressert 2012).In this work, we characterize the contributions from both ongoing star formation and the ambient radiation field in Local Group galaxy M33, as well as estimate the scale of the local dust-energy balance (i.e. the scale at which the dust is re-emitting starlight generated in that same region) in this galaxy through high-resolution radiative transfer (RT) modelling, with defined stellar and dust geometries. We have characterized the spectral energy distribution (SED) of M33 from UV to sub-mm wavelengths, at a spatial scale of 100 pc. We constructed input maps of the various stellar and dust geometries for use in the RTmodelling. By modifying our dust mix (fewer very small carbon grains and a lower silicate-to-carbon ratio as compared to the Milky Way), we can much better fit the sub-mm dust continuum. Using this new dust composition, we find that we are able to well reproduce the observed SED of M33 using our adopted model. In terms of stellar attenuation by dust, we find a reasonably strong, broad UV bump, as well as significant systematic differences in the amount of dust attenuation when compared to standard SED modelling.We also find discrepancies in the residuals of the spiral arms versus the diffuse interstellar medium (ISM), indicating a difference in properties between these two regimes. The dust emission is dominated by heating due to the young stellar populations at all wavelengths (∼80 per cent at 10 μm to ∼50 per cent at 1 mm). We find that the local dust-energy balance is restored at spatial scales greater than around 1.5 kpc.Funding from the UK Science and Technology Facilities Council consolidated grant ST/K000926/1. M.R. acknowledges support by the research projects AYA2014-53506-P and AYA2017-84897P from the Spanish Ministerio de Economía y Competitividad. This research made use of MONTAGE (http://montage.ipac.caltech.edu/), which is funded by the National Science Foundation under Grant Number ACI-1440620, and was previously funded by the National Aeronautics and Space Administration’s Earth Science Technology Office, Computation Technologies Project, under Cooperative Agreement Number NCC5-626 between NASA and the California Institute of Technolog

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 M 74 (NGC 628) and M 83 (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 M 74 and M 83, 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 m^{2 kg^{-1}} in M 74, and 0.15-0.80 m^{2 kg^{-1}} in M 83. 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

High-resolution, 3D radiative transfer modelling : V. A detailed model of the M 51 interacting pair

Context. Investigating the dust heating mechanisms in galaxies provides a deeper understanding of how the internal energy balance drives their evolution. Over the last decade radiative transfer simulations based on the Monte Carlo method have emphasised the role of the various stellar populations heating the diffuse dust. Beyond the expected heating through ongoing star formation, older stellar populations (>= 8 Gyr) and even active galactic nuclei can both contribute energy to the infrared emission of diffuse dust.Aims. In this particular study we examine how the radiation of an external heating source, such as the less massive galaxy NGC 5195 in the M 51 interacting system, could affect the heating of the diffuse dust of its parent galaxy NGC 5194, and vice versa. Our goal is to quantify the exchange of energy between the two galaxies by mapping the 3D distribution of their radiation field.Methods. We used SKIRT, a state-of-the-art 3D Monte Carlo radiative transfer code, to construct the 3D model of the radiation field of M 51, following the methodology defined in the DustPedia framework. In the interest of modelling, the assumed centre-to-centre distance separation between the two galaxies is similar to 10 kpc.Results. Our model is able to reproduce the global spectral energy distribution of the system, and it matches the resolved optical and infrared images fairly well. In total, 40.7% of the intrinsic stellar radiation of the combined system is absorbed by dust. Furthermore, we quantify the contribution of the various dust heating sources in the system, and find that the young stellar population of NGC 5194 is the predominant dust-heating agent, with a global heating fraction of 71.2%. Another 23% is provided by the older stellar population of the same galaxy, while the remaining 5.8% has its origin in NGC 5195. Locally, we find that the regions of NGC 5194 closer to NGC 5195 are significantly affected by the radiation field of the latter, with the absorbed energy fraction rising up to 38%. The contribution of NGC 5195 remains under the percentage level in the outskirts of the disc of NGC 5194. This is the first time that the heating of the diffuse dust by a companion galaxy is quantified in a nearby interacting system

Herschel observations of edge-on spirals (Heroes) IV. Dust energy balance problem

We present results of the detailed dust energy balance study for the seven large edge-on galaxies in the HEROES sample using three-dimensional (3D) radiative transfer (RT) modelling. Based on available optical and near-infrared (NIR) observations of the HEROES galaxies, we derive the 3D distribution of stars and dust in these galaxies. For the sake of uniformity, we apply the same technique to retrieve galaxy properties for the entire sample: we use a stellar model consisting of a Sérsic bulge and three double-exponential discs (a superthin disc for a young stellar population and thin and thick discs for old populations). For the dust component, we adopt a double-exponential disc with the new THEMIS dust-grain model. We fit oligochromatic RT models to the optical and NIR images with the fitting algorithm fitskirt and run panchromatic simulations with the skirt code at wavelengths ranging from ultraviolet to submillimeter. We confirm the previously stated dust energy balance problem in galaxies: for the HEROES galaxies, the dust emission derived from our RT calculations underestimates the real observations by a factor 1.5-4 for all galaxies except NGC 973 and NGC 5907 (apparently, the latter galaxy has a more complex geometry than we used). The comparison between our RT simulations and the observations at mid-infrared-submillimetre wavelengths shows that most of our galaxies exhibit complex dust morphologies (possible spiral arms, star-forming regions, more extended dust structure in the radial and vertical directions). We suggest that, in agreement with results from the literature, the large- and small-scale structure is the most probable explanation for the dust energy balance problem

High-resolution, 3D radiative transfer modelling : IV. AGN-powered dust heating in NGC 1068

The star formation rate and the mass of interstellar medium (ISM) have a high predictive power for the future evolution of a galaxy. Nevertheless, deriving such properties is not straightforward. Dust emission, an important diagnostic of star formation and ISM mass throughout the Universe, can be powered by sources unrelated to ongoing star formation. In the framework of the DustPedia project we set out to disentangle the radiation of the ongoing star formation from that of the older stellar populations. This is done through detailed 3D radiative transfer simulations of face-on spiral galaxies. We take special care in modelling the morphological features present for each source of radiation. In this particular study, we focus on NGC 1068, which in addition contains an active galactic nucleus (AGN). The effect of diffuse dust heating by an AGN (beyond the torus) has so far only been investigated for quasars. This additional dust heating source further contaminates the broadband fluxes that are used by classic galaxy modelling tools to derive physical properties. We aim to fit a realistic model to the observations of NGC 1068 and quantify the contribution of the several dust-heating sources. Our model is able to reproduce the global spectral energy distribution of the galaxy. It matches the resolved optical and infrared images fairly well, but deviates in the UV and the submillimetre (submm). This is partly due to beam smearing effects, but also because the input dust distribution is not sufficiently peaked in the centre. We find that AGN contamination of the broadband fluxes has a strong dependency on wavelength. It peaks in the mid-infrared, drops in the far-infrared, and then rises again at submm wavelengths. We quantify the contribution of the dust-heating sources in each 3D dust cell and find a median value of 83% for the star formation component. The AGN contribution is measurable at the percentage level in the disc, but quickly increases in the inner few hundred parsecs, peaking above 90%. This is the first time the phenomenon of an AGN heating the diffuse dust beyond its torus is quantified in a nearby star-forming galaxy. NGC 1068 only contains a weak AGN, meaning this effect could be stronger in galaxies with a more luminous AGN. This could significantly impact the derived star formation rates and ISM masses for such systems