Showing 1–23 of 23 results for author: Trčka, A

The many colors of the TNG100 simulation

Authors: Andrea Gebek, Ana Trčka, Maarten Baes, Marco Martorano, Annalisa Pillepich, Anand Utsav Kapoor, Angelos Nersesian, Arjen van der Wel

Abstract: We apply the 3D dust radiative transfer code SKIRT to the low-redshift ($z\leq0.1$) galaxy population in the TNG100 cosmological simulation, the fiducial run of the IllustrisTNG project. We compute global fluxes and spectral energy distributions (SEDs) from the far-ultraviolet to the sub-millimeter for $\approx\,$60 000 galaxies. Our post-processing methodology follows the study of Trčka et al. (2… ▽ More We apply the 3D dust radiative transfer code SKIRT to the low-redshift ($z\leq0.1$) galaxy population in the TNG100 cosmological simulation, the fiducial run of the IllustrisTNG project. We compute global fluxes and spectral energy distributions (SEDs) from the far-ultraviolet to the sub-millimeter for $\approx\,$60 000 galaxies. Our post-processing methodology follows the study of Trčka et al. (2022) of the higher-resolution TNG50 simulation. We verify that TNG100 reproduces observational luminosity functions at low redshifts to excellent precision, unlike TNG50. Additionally, we test the realism of our TNG100 plus SKIRT fluxes by comparing various flux and color relations to data from the GAMA survey. TNG100 broadly reproduces the observed distributions, but we predict ultraviolet colors that are too blue by $\approx\,$0.4 mag, possibly related to the extinction in the star-forming regions subgrid model not being selective enough. Furthermore, we find that the simulated galaxies exhibit mid-infrared fluxes elevated by up to $\approx\,$0.5 mag that we attribute to overly effective stochastic heating of the diffuse dust. All synthetic broadband fluxes and SEDs are made publicly available in three orientations and four apertures, and can readily be used to study TNG100 galaxies in a mock observational fashion. △ Less

Submitted 30 May, 2024; v1 submitted 8 May, 2024; originally announced May 2024.

Comments: Main text 14 pages, 5 figures, accepted to MNRAS. Generated data products are publicly available at www.tng-project.org/gebek24

The TNG50-SKIRT Atlas: wavelength dependence of the effective radius

Authors: Maarten Baes, Aleksandr Mosenkov, Raymond Kelly, Abdurro'uf, Nick Andreadis, Sena Bokona Tulu, Peter Camps, Abdissa Tassama Emana, Jacopo Fritz, Andrea Gebek, Inja Kovacic, Antonio La Marca, Marco Martorano, Angelos Nersesian, Vicente Rodriguez-Gomez, Crescenzo Tortora, Ana Trcka, Bert Vander Meulen, Arjen van der Wel, Lingyu Wang

Abstract: Galaxy sizes correlate with many other important properties of galaxies, and the cosmic evolution of galaxy sizes is an important observational diagnostic for constraining galaxy evolution models. The effective radius is probably the most widely used indicator of galaxy size. We used the TNG50-SKIRT Atlas to investigate the wavelength dependence of the effective radius of galaxies at optical and n… ▽ More Galaxy sizes correlate with many other important properties of galaxies, and the cosmic evolution of galaxy sizes is an important observational diagnostic for constraining galaxy evolution models. The effective radius is probably the most widely used indicator of galaxy size. We used the TNG50-SKIRT Atlas to investigate the wavelength dependence of the effective radius of galaxies at optical and near-infrared (NIR) wavelengths. We find that, on average, the effective radius in every band exceeds the stellar mass effective radius, and that this excess systematically decreases with increasing wavelength. The optical g-band (NIR Ks-band) effective radius is on average 58% (13%) larger than the stellar mass effective radius. Effective radii measured from dust-obscured images are systematically larger than those measured from dust-free images, although the effect is limited (8.7% in the g-band, 2.1% in the Ks-band). We find that stellar population gradients are the dominant factor (about 80%) in driving the wavelength dependence of the effective radius, and that differential dust attenuation is a secondary factor (20%). Comparing our results to recent observational data, we find offsets in the absolute values of the median effective radii, up to 50% for the population of blue galaxies. We find better agreement in the slope of the wavelength dependence of the effective radius, with red galaxies having a slightly steeper slope than green-blue galaxies. Comparing our effective radii with those of galaxies from the Siena Galaxy Atlas in separate bins in z-band absolute magnitude and g-z colour, we find excellent agreement for the reddest galaxies, but again significant offsets for the blue populations: up to 70% for galaxies around Mz=-21.5. This difference in median effective radius for the bluer galaxies is most probably due to (abridged...). △ Less

Submitted 8 January, 2024; originally announced January 2024.

Comments: 11 pages, 10 figures, accepted for publication in A&A

The TNG50-SKIRT Atlas: post-processing methodology and first data release

Authors: Maarten Baes, Andrea Gebek, Ana Trcka, Peter Camps, Arjen van der Wel, Abdurro'uf, Nick Andreadis, Sena Bokona Tulu, Abdissa Tassama Emana, Jacopo Fritz, Raymond Kelly, Inja Kovacic, Antonio La Marca, Marco Martorano, Aleksandr Mosenkov, Angelos Nersesian, Vicente Rodriguez-Gomez, Crescenzo Tortora, Bert Vander Meulen, Lingyu Wang

Abstract: Galaxy morphology is a powerful diagnostic to assess the realism of cosmological hydrodynamical simulations. Determining the morphology of simulated galaxies requires the generation of synthetic images through 3D radiative transfer post-processing that properly accounts for different stellar populations and interstellar dust attenuation. We use the SKIRT code to generate the TNG50-SKIRT Atlas, a s… ▽ More Galaxy morphology is a powerful diagnostic to assess the realism of cosmological hydrodynamical simulations. Determining the morphology of simulated galaxies requires the generation of synthetic images through 3D radiative transfer post-processing that properly accounts for different stellar populations and interstellar dust attenuation. We use the SKIRT code to generate the TNG50-SKIRT Atlas, a synthetic UV to near-infrared broadband image atlas for a complete stellar-mass selected sample of 1154 galaxies extracted from the TNG50 cosmological simulation at $z=0$. The images have a high spatial resolution (100 pc) and a wide field of view (160 kpc). In addition to the dust-obscured images, we also release dust-free images and physical parameter property maps with matching characteristics. As a sanity check and preview application we discuss the UVJ diagram of the galaxy sample. We investigate the effect of dust attenuation on the UVJ diagram and find that it affects both the star-forming and the quiescent galaxy populations. The quiescent galaxy region is polluted by younger and star-forming highly inclined galaxies, while dust attenuation induces a separation in inclination of the star-forming galaxy population, with low-inclination galaxies remaining at the blue side of the diagram and high-inclination galaxies systematically moving towards the red side. This image atlas can be used for a variety of other applications, including galaxy morphology studies and the investigation of local scaling relations. We publicly release the images and parameter maps, and we invite the community to use them. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Comments: 13 pages, 7 figures, accepted for publication in A&A

UV to submillimetre luminosity functions of TNG50 galaxies

Authors: Ana Trcka, Maarten Baes, Peter Camps, Anand Utsav Kapoor, Dylan Nelson, Annalisa Pillepich, Daniela Barrientos, Lars Hernquist, Federico Marinacci, Mark Vogelsberger

Abstract: We apply the radiative transfer (RT) code SKIRT on a sample of ~14000 low-redshift (z<= 0.1) galaxies extracted from the TNG50 simulation to enable an apples-to-apples comparison with observations. The RT procedure is calibrated via comparison of a subsample of TNG50 galaxies with the DustPedia observational sample: we compare several luminosity and colour scaling relations and spectral energy dis… ▽ More We apply the radiative transfer (RT) code SKIRT on a sample of ~14000 low-redshift (z<= 0.1) galaxies extracted from the TNG50 simulation to enable an apples-to-apples comparison with observations. The RT procedure is calibrated via comparison of a subsample of TNG50 galaxies with the DustPedia observational sample: we compare several luminosity and colour scaling relations and spectral energy distributions in different specific SFR bins. We consistently derive galaxy luminosity functions for the TNG50 simulation in 14 broadband filters from UV to submillimetre wavelengths and investigate the effects of the aperture, orientation, radiative transfer recipe, and numerical resolution. We find that, while our TNG50+RT fiducial model agrees well with the observed luminosity functions at the knee (+/- 0.04 dex typical agreement), the TNG50+RT luminosity functions evaluated within 5R_1/2 are generally higher than observed at both the faint and bright ends, by 0.004 (total IR)-0.27 (UKIDSS H) dex and 0.12 (SPIRE250)-0.8 (GALEX FUV) dex, respectively. A change in the aperture does affect the bright end of the luminosity function, easily by up to 1 dex depending on the choice. However, we also find that the galaxy luminosity functions of a worse-resolution run of TNG50 (TNG50-2, with 8 times worse mass resolution than TNG50, similar to TNG100) are in better quantitative agreement with observational constraints. Finally, we publicly release the photometry for the TNG50 sample in 53 broadbands from FUV to submillimetre, in three orientations and four apertures, as well as galaxy spectral energy distributions. △ Less

Submitted 12 August, 2022; originally announced August 2022.

Comments: The full galaxy SEDs and the broadband fluxes at z <= 0.1 are made available at https://www.tng- project.org/trcka22 The luminosity function plots are presented on http://luminosity-functions.herokuapp.com/

High-resolution synthetic UV-submm images for Milky Way-mass simulated galaxies from the ARTEMIS project

Authors: Peter Camps, Anand Utsav Kapoor, Ana Trčka, Andreea S. Font, Ian G. McCarthy, James Trayford, Maarten Baes

Abstract: We present redshift-zero synthetic dust-aware observations for the 45 Milky Way-mass simulated galaxies of the ARTEMIS project, calculated with the SKIRT radiative transfer code. The post-processing procedure includes components for star-forming regions, stellar sources, and diffuse dust. We produce and publicly release realistic high-resolution images for 50 commonly-used broadband filters from u… ▽ More We present redshift-zero synthetic dust-aware observations for the 45 Milky Way-mass simulated galaxies of the ARTEMIS project, calculated with the SKIRT radiative transfer code. The post-processing procedure includes components for star-forming regions, stellar sources, and diffuse dust. We produce and publicly release realistic high-resolution images for 50 commonly-used broadband filters from ultraviolet to sub-millimetre wavelengths and for 18 different viewing angles. We compare the simulated ARTEMIS galaxies to observed galaxies in the DustPedia database with similar stellar mass and star formation rate, and to synthetic observations of the simulated galaxies of the Auriga project produced in previous work using a similar post-processing technique. In all cases, global galaxy properties are derived using SED fitting. We find that, similar to Auriga, the post-processed ARTEMIS galaxies generally reproduce the observed scaling relations for global fluxes and physical properties, although dust extinction at FUV/UV wavelengths is underestimated and representative dust temperatures are lower than observed. At a resolved scale, we compare multi-wavelength non-parametric morphological properties of selected disc galaxies across the data sets. We find that the ARTEMIS galaxies largely reproduce the observed morphological trends as a function of wavelength, although they appear to be more clumpy and less symmetrical than observed. We note that the ARTEMIS and Auriga galaxies occupy adjacent regions in the specific star formation versus stellar mass plane, so that the synthetic observation data sets supplement each other. △ Less

Submitted 11 March, 2022; originally announced March 2022.

Comments: Accepted for publication by MNRAS

High-resolution synthetic UV-submm images for simulated Milky Way-type galaxies from the Auriga project

Authors: Anand Utsav Kapoor, Peter Camps, Maarten Baes, Ana Trcka, Robert J. J. Grand, Arjen van der Wel, Luca Cortese, Ilse De Looze, Daniela Barrientos

Abstract: We present redshift-zero synthetic observational data considering dust attenuation and dust emission for the thirty galaxies of the Auriga project, calculated with the SKIRT radiative transfer code. The post-processing procedure includes components for star-forming regions, stellar sources, and diffuse dust taking into account stochastic heating of dust grains. This allows us to obtain realistic h… ▽ More We present redshift-zero synthetic observational data considering dust attenuation and dust emission for the thirty galaxies of the Auriga project, calculated with the SKIRT radiative transfer code. The post-processing procedure includes components for star-forming regions, stellar sources, and diffuse dust taking into account stochastic heating of dust grains. This allows us to obtain realistic high-resolution broadband images and fluxes from ultraviolet to sub-millimeter wavelengths. For the diffuse dust component, we consider two mechanisms for assigning dust to gas cells in the simulation. In one case, only the densest or the coldest gas cells are allowed to have dust, while in the other case this condition is relaxed to allow a larger number of dust-containing cells. The latter approach yields galaxies with a larger radial dust extent and an enhanced dust presence in the inter-spiral regions. At a global scale, we compare Auriga galaxies with observations by deriving dust scaling relations using SED fitting. At a resolved scale, we make a multi-wavelength morphological comparison with nine well-resolved spiral galaxies from the DustPedia observational database. We find that for both dust assignment methods, although the Auriga galaxies show a good overall agreement with observational dust properties, they exhibit a slightly higher specific dust mass. The multi-wavelength morphological analysis reveals a good agreement between the Auriga and the observed galaxies in the optical wavelengths. In the mid and far-infrared wavelengths, Auriga galaxies appear smaller and more centrally concentrated in comparison to their observed counterparts. We publicly release the multi-observer images and fluxes in 50 commonly used broadband filters. △ Less

Submitted 14 July, 2021; originally announced July 2021.

Comments: 19 pages, 20 figures, accepted for publication in MNRAS

Probing the spectral shape of dust emission with the DustPedia galaxy sample

Authors: Angelos Nersesian, Wouter Dobbels, Manolis E. Xilouris, Maarten Baes, Simone Bianchi, Viviana Casasola, Christopher J. R. Clark, Ilse De Looze, Frédéric Galliano, Suzanne C. Madden, Aleksandr V. Mosenkov, Evangelos-D. Paspaliaris, Ana Trčka

Abstract: The objective of this paper is to understand the variance of the far-infrared (FIR) spectral energy distribution (SED) of the DustPedia galaxies, and its link with the stellar and dust properties. An interesting aspect of the dust emission is the inferred FIR colours which could inform us about the dust content of galaxies, and how it varies with the physical conditions within galaxies. However, t… ▽ More The objective of this paper is to understand the variance of the far-infrared (FIR) spectral energy distribution (SED) of the DustPedia galaxies, and its link with the stellar and dust properties. An interesting aspect of the dust emission is the inferred FIR colours which could inform us about the dust content of galaxies, and how it varies with the physical conditions within galaxies. However, the inherent complexity of dust grains as well as the variety of physical properties depending on dust, hinder our ability to utilise their maximum potential. We use principal component analysis (PCA) to explore new hidden correlations with many relevant physical properties such as the dust luminosity, dust temperature, dust mass, bolometric luminosity, star-formation rate (SFR), stellar mass, specific SFR, dust-to-stellar mass ratio, the fraction of absorbed stellar luminosity by dust (f_abs), and metallicity. We find that 95% of the variance in our sample can be described by two principal components (PCs). The first component controls the wavelength of the peak of the SED, while the second characterises the width. The physical quantities that correlate better with the coefficients of the first two PCs, and thus control the shape of the FIR SED are: the dust temperature, the dust luminosity, the SFR, and f_abs. Finally, we find a weak tendency for low-metallicity galaxies to have warmer and broader SEDs, while on the other hand high-metallicity galaxies have FIR SEDs that are colder and narrower. △ Less

Submitted 9 July, 2021; originally announced July 2021.

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

Authors: Angelos Nersesian, Sebastien Viaene, Ilse De Looze, Maarten Baes, Emmanuel M. Xilouris, Matthew W. L. Smith, Simone Bianchi, Viviana Casasola, Letizia P. Cassara, Christopher J. R. Clark, Wouter Dobbels, Jacopo Fritz, Frederic Galliano, Suzanne C. Madden, Aleksandr V. Mosenkov, Ana Trcka

Abstract: 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 underlined the role of the various stellar populations heating the diffuse dust. Beyond the expected heating through ongoing star formation, both older stellar popul… ▽ More 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 underlined the role of the various stellar populations heating the diffuse dust. Beyond the expected heating through ongoing star formation, both older stellar population (> 8Gyr) and even AGN can contribute energy to the infrared emission of diffuse dust. Here, we examine how the radiation of an external heating source, like the less massive galaxy NGC5195, in the M51 interacting system, could affect the heating of the diffuse dust of its parent galaxy, NGC5194, and vice versa. To quantify the exchange of energy between the two galaxies we use SKIRT, a state-of-the-art Monte Carlo radiative transfer code. In the interest of modelling, the assumed centre-to-centre distance separation between the two galaxies is 10kpc. Our model reproduces the global SED of the system, and it closely matches the observed images. 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 NGC5194 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 NGC5195. Locally, we find that the regions of NGC5194 closer to NGC5195 are significantly affected by the radiation field of the latter, with the absorbed energy fraction rising up to 38%. The contribution of NGC5195 remains under the percentage level in the outskirts of the disc of NGC5194. This is the first time that the heating of the diffuse dust by a companion galaxy is quantified in a nearby interacting system. △ Less

Submitted 18 September, 2020; v1 submitted 15 September, 2020; originally announced September 2020.

Journal ref: A&A 643, A90 (2020)

Nonparametric galaxy morphology from UV to submm wavelengths

Authors: Maarten Baes, Angelos Nersesian, Viviana Casasola, Simone Bianchi, Letizia P. Cassarà, Christopher J. R. Clark, Ilse De Looze, Wouter Dobbels, Jacopo Fritz, Maud Galametz, Frédéric Galliano, Suzanne C. Madden, Aleksandr V. Mosenkov, Sébastien Viaene, Ana Trčka, Emmanuel M. Xilouris

Abstract: 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 c… ▽ More 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. △ Less

Submitted 17 July, 2020; v1 submitted 4 July, 2020; originally announced July 2020.

Comments: 9 pages, 5 figures, accepted for publication in A&A

Journal ref: A&A 641, A119 (2020)

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

Authors: S. Viaene, A. Nersesian, J. Fritz, S. Verstocken, M. Baes, S. Bianchi, V. Casasola, L. Cassarà, C. Clark, J. Davies, I. De Looze, P. De Vis, W. Dobbels, M. Galametz, F. Galliano, A. Jones, S. Madden, A. Mosenkov, A. Trcka, E. M. Xilouris, N. Ysard

Abstract: 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 have 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 spira… ▽ More 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 have 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. In this particular study, we focus on NGC 1068, which contains an active galactic nucleus (AGN). The effect of diffuse dust heating by AGN (beyond the torus) was so far only investigated for quasars. This additional dust heating source further contaminates the broadband fluxes on which classic galaxy modelling tools rely 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 submm. We find a strong wavelength dependency of AGN contamination to the broadband fluxes. It peaks in the MIR, drops in the FIR, but 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 100 pc, 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 can be stronger in galaxies with a more luminous AGN. This could significantly impact the derived star formation rates and ISM masses for such systems. △ Less

Submitted 4 May, 2020; originally announced May 2020.

Comments: 18 pages, 11 figures, accepted for publication in A&A

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

Authors: Angelos Nersesian, Sam Verstocken, Sebastien Viaene, Maarten Baes, Emmanuel M. Xilouris, Simone Bianchi, Viviana Casasola, Christopher J. R. Clark, Jonathan I. Davies, Ilse De Looze, Pieter De Vis, Wouter Dobbels, Jacopo Fritz, Maud Galametz, Frederic Galliano, Anthony P. Jones, Suzanne C. Madden, Aleksandr V. Mosenkov, Ana Trcka, Nathalie Ysard

Abstract: 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… ▽ More 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. △ Less

Submitted 7 April, 2020; originally announced April 2020.

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

Authors: Sam Verstocken, Angelos Nersesian, Maarten Baes, Sébastien Viaene, Simone Bianchi, Viviana Casasola, Christopher J. R. Clark, Jonathan I. Davies, Ilse De Looze, Pieter De Vis, Wouter Dobbels, FrédÉric Galliano, Anthony P. Jones, Suzanne C. Madden, Aleksandr V. Mosenkov, Ana Trčka, Emmanuel M. Xilouris

Abstract: 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 stel… ▽ More 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) △ Less

Submitted 7 April, 2020; originally announced April 2020.

Comments: 23 pages, 16 figures, accepted for publication in A&A

Journal ref: A&A 637, A24 (2020)

Infrared luminosity functions and dust mass functions in the EAGLE simulation

Authors: Maarten Baes, Ana Trčka, Peter Camps, James Trayford, Antonios Katsianis, Lucia Marchetti, Tom Theuns, Mattia Vaccari, Bert Vandenbroucke

Abstract: We present infrared luminosity functions and dust mass functions for the EAGLE cosmological simulation, based on synthetic multi-wavelength observations generated with the SKIRT radiative transfer code. In the local Universe, we reproduce the observed infrared luminosity and dust mass functions very well. Some minor discrepancies are encountered, mainly in the high luminosity regime, where the EAG… ▽ More We present infrared luminosity functions and dust mass functions for the EAGLE cosmological simulation, based on synthetic multi-wavelength observations generated with the SKIRT radiative transfer code. In the local Universe, we reproduce the observed infrared luminosity and dust mass functions very well. Some minor discrepancies are encountered, mainly in the high luminosity regime, where the EAGLE-SKIRT luminosity functions mildly but systematically underestimate the observed ones. The agreement between the EAGLE-SKIRT infrared luminosity functions and the observed ones gradually worsens with increasing lookback time. Fitting modified Schechter functions to the EAGLE-SKIRT luminosity and dust mass functions at different redshifts up to $z=1$, we find that the evolution is compatible with pure luminosity/mass evolution. The evolution is relatively mild: within this redshift range, we find an evolution of $L_{\star,250}\propto(1+z)^{1.68}$, $L_{\star,\text{TIR}}\propto(1+z)^{2.51}$ and $M_{\star,\text{dust}}\propto(1+z)^{0.83}$ for the characteristic luminosity/mass. For the luminosity/mass density we find $\varepsilon_{250}\propto(1+z)^{1.62}$, $\varepsilon_{\text{TIR}}\propto(1+z)^{2.35}$ and $ρ_{\text{dust}}\propto(1+z)^{0.80}$, respectively. The mild evolution of the dust mass density is in relatively good agreement with observations, but the slow evolution of the infrared luminosity underestimates the observed luminosity evolution significantly. We argue that these differences can be attributed to increasing limitations in the radiative transfer treatment due to increasingly poorer resolution, combined with a slower than observed evolution of the SFR density in the EAGLE simulation and the lack of AGN emission in our EAGLE-SKIRT post-processing recipe. △ Less

Submitted 6 April, 2020; originally announced April 2020.

Comments: 13 pages, 6 pages, accepted for publication in MNRAS

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

Authors: Ana Trčka, Maarten Baes, Peter Camps, Sharon E. Meidt, James Trayford, Simone Bianchi, Viviana Casasola, Letizia P. Cassarà, Ilse De Looze, Pieter De Vis, Wouter Dobbels, Jacopo Fritz, Maud Galametz, Frédéric Galliano, Antonios Katsianis, Suzanne C. Madden, Aleksandr V. Mosenkov, Angelos Nersesian, Sébastien Viaene, Emmanuel M. Xilouris

Abstract: 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 o… ▽ More 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. △ Less

Submitted 27 March, 2020; originally announced March 2020.

Comments: 17 pages, 14 figures, accepted for publication in MNRAS

The high redshift SFR-M* relation is sensitive to the employed star formation rate and stellar mass indicators: Towards addressing the tension between observations and simulations

Authors: A. Katsianis, V Gonzalez, D. Barrientos, X. Yang, C. D. P. Lagos, J. Schaye, P. Camps, A. Trčka, M. Baes, M. Stalevski, G. A. Blanc, T. Theuns

Abstract: There is a severe tension between the observed star formation rate (SFR) - stellar mass (${\rm M}_{\star}$) relations reported by different authors at $z = 1-4$. In addition, the observations have not been successfully reproduced by state-of-the-art cosmological simulations which tend to predict a factor of 2-4 smaller SFRs at a fixed ${\rm M}_{\star}$. We examine the evolution of the SFR… ▽ More There is a severe tension between the observed star formation rate (SFR) - stellar mass (${\rm M}_{\star}$) relations reported by different authors at $z = 1-4$. In addition, the observations have not been successfully reproduced by state-of-the-art cosmological simulations which tend to predict a factor of 2-4 smaller SFRs at a fixed ${\rm M}_{\star}$. We examine the evolution of the SFR$-{\rm M}_{\star}$ relation of $z = 1-4 $ galaxies using the SKIRT simulated spectral energy distributions of galaxies sampled from the EAGLE simulations. We derive SFRs and stellar masses by mimicking different observational techniques. We find that the tension between observed and simulated SFR$-{\rm M}_{\star}$ relations is largely alleviated if similar methods are used to infer the galaxy properties. We find that relations relying on infrared wavelengths (e.g. 24 ${\rm μm}$, MIPS - 24, 70 and 160 ${\rm μm}$ or SPIRE - 250, 350, 500 ${\rm μm}$) have SFRs that exceed the intrinsic relation by 0.5 dex. Relations that rely on the spectral energy distribution fitting technique underpredict the SFRs at a fixed stellar mass by -0.5 dex at $z \sim 4$ but overpredict the measurements by 0.3 dex at $z \sim 1$. Relations relying on dust-corrected rest-frame UV luminosities, are flatter since they overpredict/underpredict SFRs for low/high star forming objects and yield deviations from the intrinsic relation from 0.10 dex to -0.13 dex at $z \sim 4$. We suggest that the severe tension between different observational studies can be broadly explained by the fact that different groups employ different techniques to infer their SFRs. △ Less

Submitted 16 January, 2020; originally announced January 2020.

Comments: 15 pages, 6 Figures, 4 tables, accepted for publication in MNRAS, Figures 1, 2, 3, 4 summarize results

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

Authors: Viviana Casasola, Simone Bianchi, Pieter De Vis, Laura Magrini, Edvige Corbelli, Christopher J. R. Clark, Jacopo Fritz, Angelos Nersesian, Sebastien Viaene, Maarten Baes, Letizia P. Cassara', Jon Davies, Ilse De Looze, Wouter Dobbels, Maud Galametz, Frederic Galliano, Anthony P. Jones, Suzanne C. Madden, Aleksandr V. Mosenkov, Ana Trcka, Emmanuel Xilouris

Abstract: 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 a… ▽ More 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. △ Less

Submitted 19 December, 2019; v1 submitted 20 November, 2019; originally announced November 2019.

Comments: Accepted for publication in A&A, 21 pages, 14 Figures, 7 Tables, 2 Appendices (with 3 Tables). Published version, unchanged results with respect to the previously posted version

Journal ref: A&A 633, A100 (2020)

Predicting the global far-infrared SED of galaxies via machine learning techniques

Authors: W. Dobbels, M. Baes, S. Viaene, S. Bianchi, J. I. Davies, V. Casasola, C. J. R. Clark, J. Fritz, M. Galametz, F. Galliano, A. Mosenkov, A. Nersesian, A. Trčka

Abstract: Dust plays an important role in shaping a galaxy's spectral energy distribution (SED). It absorbs ultraviolet (UV) to near-infrared (NIR) radiation and re-emits this energy in the far-infrared (FIR). The FIR is essential to understand dust in galaxies. However, deep FIR observations require a space mission, none of which are still active today. We aim to infer the FIR emission across six Herschel… ▽ More Dust plays an important role in shaping a galaxy's spectral energy distribution (SED). It absorbs ultraviolet (UV) to near-infrared (NIR) radiation and re-emits this energy in the far-infrared (FIR). The FIR is essential to understand dust in galaxies. However, deep FIR observations require a space mission, none of which are still active today. We aim to infer the FIR emission across six Herschel bands, along with dust luminosity, mass, and effective temperature, based on the available UV to mid-infrared (MIR) observations. We also want to estimate the uncertainties of these predictions, compare our method to energy balance SED fitting, and determine possible limitations of the model. We propose a machine learning framework to predict the FIR fluxes from 14 UV-MIR broadband fluxes. We used a low redshift sample by combining DustPedia and H-ATLAS, and extracted Bayesian flux posteriors through SED fitting. We trained shallow neural networks to predict the far-infrared fluxes, uncertainties, and dust properties. We evaluated them on a test set using a root mean square error (RMSE) in log-space. Our results (RMSE = 0.19 dex) significantly outperform UV-MIR energy balance SED fitting (RMSE = 0.38 dex), and are inherently unbiased. We can identify when the predictions are off, for example when the input has large uncertainties on WISE 22, or when the input does not resemble the training set. The galaxies for which we have UV-FIR observations can be used as a blueprint for galaxies that lack FIR data. This results in a 'virtual FIR telescope', which can be applied to large optical-MIR galaxy samples. This helps bridge the gap until the next FIR mission. △ Less

Submitted 14 October, 2019; originally announced October 2019.

Comments: 25 pages, 16 figures, accepted by A&A. Interactive figures + jupyter notebooks at https://wdobbels.github.io/FIREnet/

Journal ref: A&A 634, A57 (2020)

An evolving and mass dependent ${\rm σ_{\text{sSFR}}}$-${\rm M_{\star}}$ relation for galaxies

Authors: Antonios Katsianis, Xianzhong Zheng, Valentino Gonzalez, Guillermo Blanc, Claudia del P. Lagos, Luke J. M. Davies, Peter Camps, Ana Trčka, Maarten Baes, Joop Schaye, James W. Trayford, Tom Theuns, Marko Stalevski

Abstract: The scatter (${\rmσ_{\text{sSFR}}}$) of the specific star formation rates (sSFRs) of galaxies is a measure of the diversity in their star formation histories (SFHs) at a given mass. In this paper we employ the EAGLE simulations to study the dependence of the ${\rm σ_{\text{sSFR}}}$ of galaxies on stellar mass (${\rm M_{\star}}$) through the ${\rm σ_{\text{sSFR}}}$-${\rm M_{\star}}$ relation in… ▽ More The scatter (${\rmσ_{\text{sSFR}}}$) of the specific star formation rates (sSFRs) of galaxies is a measure of the diversity in their star formation histories (SFHs) at a given mass. In this paper we employ the EAGLE simulations to study the dependence of the ${\rm σ_{\text{sSFR}}}$ of galaxies on stellar mass (${\rm M_{\star}}$) through the ${\rm σ_{\text{sSFR}}}$-${\rm M_{\star}}$ relation in $ {\rm z \sim 0-4}$. We find that the relation evolves with time, with the dispersion depending on both stellar mass and redshift. The models point to an evolving U-shape form for the ${\rm σ_{\text{sSFR}}}$-${\rm M_{\star}}$ relation with the scatter being minimal at a characteristic mass $M^{\star}$ of ${\rm 10^{9.5}}$ ${\rm M_{\odot}}$ and increasing both at lower and higher masses. This implication is that the diversity of SFHs increases towards both at the low- and high-mass ends. We find that active galactic nuclei feedback is important for increasing the ${\rm σ_{\text{sSFR}}}$ for high mass objects. On the other hand, we suggest that SNe feedback increases the ${\rm σ_{\text{sSFR}}}$ of galaxies at the low-mass end. We also find that excluding galaxies that have experienced recent mergers does not significantly affect the ${\rm σ_{\text{sSFR}}}$-${\rm M_{\star}}$ relation. Furthermore, we employ the combination of the EAGLE simulations with the radiative transfer code SKIRT to evaluate the effect of SFR/stellar mass diagnostics in the ${\rm σ_{\text{sSFR}}}$-${\rm M_{\star}}$ relation and find that the ${\rm SFR/M_{\star}}$ methodologies (e.g. SED fitting, UV+IR, UV+IRX-$β$) widely used in the literature to obtain intrinsic properties of galaxies have a large effect on the derived shape and normalization of the ${\rm σ_{\text{sSFR}}}$-${\rm M_{\star}}$ relation. △ Less

Submitted 6 May, 2019; originally announced May 2019.

Comments: 23 pages, 10 figures, Accepted to ApJ

Old and young stellar populations in DustPedia galaxies and their role in dust heating

Authors: A. Nersesian, E. M. Xilouris, S. Bianchi, F. Galliano, A. P. Jones, M. Baes, V. Casasola, L. P. Cassara, C. J. R. Clark, J. I. Davies, M. Decleir, W. Dobbels, I. De Looze, P. De Vis, J. Fritz, M. Galametz, S. C. Madden, A. V. Mosenkov, A. Trcka, S. Verstocken, S. Viaene, S. Lianou

Abstract: Within the framework of the DustPedia project we investigate the properties of cosmic dust and its interaction with the stellar radiation (originating from different stellar populations) for 814 galaxies in the nearby Universe, all observed by the Herschel Space Observatory. We take advantage of the widely used galaxy SED fitting code CIGALE, properly adapted to include the state-of-the-art dust m… ▽ More Within the framework of the DustPedia project we investigate the properties of cosmic dust and its interaction with the stellar radiation (originating from different stellar populations) for 814 galaxies in the nearby Universe, all observed by the Herschel Space Observatory. We take advantage of the widely used galaxy SED fitting code CIGALE, properly adapted to include the state-of-the-art dust model THEMIS. Using the DustPedia photometry we determine the physical properties of the galaxies, such as, the dust and stellar mass, the star-formation rate, the bolometric luminosity as well as the unattenuated and the absorbed by dust stellar light, for both the old (> 200 Myr) and young (<= 200 Myr) stellar populations. We show how the mass of stars, dust, and atomic gas, as well as the star-formation rate and the dust temperature vary between galaxies of different morphologies and provide recipes to estimate these parameters given their Hubble stage (T). We find a mild correlation between the mass fraction of the small a-C(:H) grains with the specific star-formation rate. On average, young stars are very efficient in heating the dust, with absorption fractions reaching as high as ~77% of the total, unattenuated luminosity of this population. On the other hand, the maximum absorption fraction of old stars is ~24%. Dust heating in early-type galaxies is mainly due to old stars, up to a level of ~90%. Young stars progressively contribute more for `typical' spiral galaxies and they become the dominant source of dust heating for Sm type and irregular galaxies, donating up to ~60% of their luminosity to this purpose. Finally, we find a strong correlation of the dust heating fraction by young stars with morphology and the specific star-formation rate. △ Less

Submitted 16 April, 2019; v1 submitted 14 March, 2019; originally announced March 2019.

Comments: 22 pages, 18 figures

Journal ref: A&A 624, A80 (2019)

The cosmic spectral energy distribution in the EAGLE simulation

Authors: Maarten Baes, Ana Trčka, Peter Camps, Angelos Nersesian, James Trayford, Tom Theuns, Wouter Dobbels

Abstract: The cosmic spectral energy distribution (CSED) is the total emissivity as a function of wavelength of galaxies in a given cosmic volume. We compare the observed CSED from the UV to the submm to that computed from the EAGLE cosmological hydrodynamical simulation, post-processed with stellar population synthesis models and including dust radiative transfer using the SKIRT code. The agreement with th… ▽ More The cosmic spectral energy distribution (CSED) is the total emissivity as a function of wavelength of galaxies in a given cosmic volume. We compare the observed CSED from the UV to the submm to that computed from the EAGLE cosmological hydrodynamical simulation, post-processed with stellar population synthesis models and including dust radiative transfer using the SKIRT code. The agreement with the data is better than 0.15 dex over the entire wavelength range at redshift $z=0$, except at UV wavelengths where the EAGLE model overestimates the observed CSED by up to a factor 2. Global properties of the CSED as inferred from CIGALE fits, such as the stellar mass density, mean star formation density, and mean dust-to-stellar-mass ratio, agree to within better than 20 per cent. At higher redshift, EAGLE increasingly underestimates the CSED at optical-NIR wavelengths with the FIR/submm emissivity underestimated by more than a factor of 5 by redshift $z=1$. We believe that these differences are due to a combination of incompleteness of the EAGLE-SKIRT database, the small simulation volume and the consequent lack of luminous galaxies, and our lack of knowledge on the evolution of the characteristics of the interstellar dust in galaxies. The impressive agreement between the simulated and observed CSED at lower $z$ confirms that the combination of EAGLE and SKIRT dust processing yields a fairly realistic representation of the local Universe. △ Less

Submitted 28 January, 2019; v1 submitted 25 January, 2019; originally announced January 2019.

Comments: 14 pages, 5 figures, accepted for publication in MNRAS

Dust emission profiles of DustPedia galaxies

Authors: A. V. Mosenkov, M. Baes, S. Bianchi, V. Casasola, L. P. Cassarà, C. J. R. Clark, J. Davies, I. De Looze, P. De Vis, J. Fritz, M. Galametz, F. Galliano, A. P. Jones, S. Lianou, S. C. Madden, A. Nersesian, M. W. L. Smith, A. Trčka, S. Verstocken, S. Viaene, M. Vika, E. Xilouris

Abstract: Most radiative transfer models assume that dust in spiral galaxies is distributed exponentially. In this paper our goal is to verify this assumption by analysing the two-dimensional large-scale distribution of dust in galaxies from the DustPedia sample. For this purpose, we make use of Herschel imaging in five bands, from 100 to 500μm, in which the cold dust constituent is primarily traced and mak… ▽ More Most radiative transfer models assume that dust in spiral galaxies is distributed exponentially. In this paper our goal is to verify this assumption by analysing the two-dimensional large-scale distribution of dust in galaxies from the DustPedia sample. For this purpose, we make use of Herschel imaging in five bands, from 100 to 500μm, in which the cold dust constituent is primarily traced and makes up the bulk of the dust mass in spiral galaxies. For a subsample of 320 disc galaxies, we successfully perform a simultaneous fitting with a single Sérsic model of the Herschel images in all five bands using the multiband modelling code GALFITM. We report that the Sérsic index $n$, which characterises the shape of the Sérsic profile, lies systematically below 1 in all Herschel bands and is almost constant with wavelength. The average value at 250μm is $0.67\pm0.37$ (187 galaxies are fitted with $n_{250}\leq0.75$, 87 galaxies have $0.75<n_{250}\leq1.25$, and 46 - with $n_{250}>1.25$). Most observed profiles exhibit a depletion in the inner region (at $r<0.3-0.4$ of the optical radius $r_{25}$ ) and are more or less exponential in the outer part. We also find breaks in the dust emission profiles at longer distances $(0.5-0.6)r_{25}$ which are associated with the breaks in the optical and near-infrared. We assume that the observed deficit of dust emission in the inner galaxy region is related to the depression in the radial profile of the HI surface density in the same region because the atomic gas reaches high enough surface densities there to be transformed into molecular gas. If a galaxy has a triggered star formation in the inner region (for example, because of a strong bar instability, which transfers the gas inwards to the centre, or a pseudobulge formation), no depletion or even an excess of dust emission in the centre is observed. △ Less

Submitted 21 November, 2018; originally announced November 2018.

Comments: 24 pages, 17 figures, Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 622, A132 (2019)

The fraction of bolometric luminosity absorbed by dust in DustPedia galaxies

Authors: S. Bianchi, P. De Vis, S. Viaene, A. Nersesian, A. V. Mosenkov, E. M. Xilouris, M. Baes, V. Casasola, L. P. Cassarà, C. J. R. Clark, J. I. Davies, I. De Looze, W. Dobbels, M. Galametz, F. Galliano, A. P. Jones, S. Lianou, S. C. Madden, A. Trčka

Abstract: We study the fraction of stellar radiation absorbed by dust, f_abs, in 814 galaxies of different morphological types. The targets constitute the vast majority (93%) of the DustPedia sample, including almost all large (optical diameter larger than 1'), nearby (v <= 3000 km/s) galaxies observed with the Herschel Space Observatory. For each object, we model the spectral energy distribution from the u… ▽ More We study the fraction of stellar radiation absorbed by dust, f_abs, in 814 galaxies of different morphological types. The targets constitute the vast majority (93%) of the DustPedia sample, including almost all large (optical diameter larger than 1'), nearby (v <= 3000 km/s) galaxies observed with the Herschel Space Observatory. For each object, we model the spectral energy distribution from the ultraviolet to the sub-millimetre using the dedicated, aperture-matched DustPedia photometry and the fitting code CIGALE. The value of f_abs is obtained from the total luminosity emitted by dust and from the bolometric luminosity, which are estimated by the fit. On average, 19% of the stellar radiation is absorbed by dust in DustPedia galaxies. The fraction rises to 25% if only late-type galaxies are considered. The dependence of f_abs on morphology, showing a peak for Sb-Sc galaxies, is weak; it reflects a stronger, yet broad, positive correlation with the bolometric luminosity, which is identified for late-type, disk-dominated, high-specific-star-formation rate, gas-rich objects. We find no variation of f_abs with inclination, at odds with radiative transfer models of edge-on galaxies. These results call for a self-consistent modelling of the evolution of the dust mass and geometry along the build-up of the stellar content. We also provide template spectral energy distributions in bins of morphology and luminosity and study the variation of f_abs with stellar mass and specific star formation rate. We confirm that the local Universe is missing the high f_abs}, luminous and actively star-forming objects necessary to explain the energy budget in observations of the extragalactic background light. △ Less

Submitted 2 October, 2018; originally announced October 2018.

Comments: 21 pages, 14 figures, A&A accepted. Template SEDs available at the DustPedia Archive: http://dustpedia.astro.noa.gr/

Journal ref: A&A 620, A112 (2018)

Data release of UV to submm broadband fluxes for simulated galaxies from the EAGLE project

Authors: Peter Camps, Ana Trčka, James Trayford, Maarten Baes, Tom Theuns, Robert A. Crain, Stuart McAlpine, Matthieu Schaller, Joop Schaye

Abstract: We present dust-attenuated and dust emission fluxes for sufficiently resolved galaxies in the EAGLE suite of cosmological hydrodynamical simulations, calculated with the SKIRT radiative transfer code. The post-processing procedure includes specific components for star formation regions, stellar sources, and diffuse dust, and takes into account stochastic heating of dust grains to obtain realistic… ▽ More We present dust-attenuated and dust emission fluxes for sufficiently resolved galaxies in the EAGLE suite of cosmological hydrodynamical simulations, calculated with the SKIRT radiative transfer code. The post-processing procedure includes specific components for star formation regions, stellar sources, and diffuse dust, and takes into account stochastic heating of dust grains to obtain realistic broad-band fluxes in the wavelength range from ultraviolet to sub-millimeter. The mock survey includes nearly half a million simulated galaxies with stellar masses above 10^8.5 solar masses across six EAGLE models. About two thirds of these galaxies, residing in 23 redshift bins up to z=6, have a sufficiently resolved metallic gas distribution to derive meaningful dust attenuation and emission, with the important caveat that the same dust properties were used at all redshifts. These newly released data complement the already publicly available information about the EAGLE galaxies, which includes intrinsic properties derived by aggregating the properties of the smoothed particles representing matter in the simulation. We further provide an open source framework of Python procedures for post-processing simulated galaxies with the radiative transfer code SKIRT. The framework allows any third party to calculate synthetic images, SEDs, and broadband fluxes for EAGLE galaxies, taking into account the effects of dust attenuation and emission. △ Less

Submitted 15 December, 2017; originally announced December 2017.

Comments: 20 pages, 10 figures, accepted for publication in ApJS