Search

We present spatially resolved, simple stellar population equivalent ages, stellar metallicities, and abundance ratios for 456 massive ($10.3\lesssim\log(\mathrm{M}_*/\mathrm{M}_\odot)\lesssim11.8$) quiescent galaxies at $0.6\lesssim z\lesssim1.0$ from the Large Early Galaxy Astrophysics Census, derived using full-spectrum models. Typically, we find flat age and [Mg/Fe] gradients, and negative [Fe/H] gradients, implying iron-rich cores. We also estimate intrinsic [Fe/H] gradients via forward modelling. We examine the observed gradients in three age bins. Younger quiescent galaxies typically have negative [Fe/H] gradients and positive age gradients, possibly indicating a recent central starburst. Additionally, this finding suggests that photometrically measured flat colour gradients in young quiescent galaxies are the result of the positive age and negative metallicity gradients cancelling each other. For older quiescent galaxies, the age gradients become flat and [Fe/H] gradients weaken, though remain negative. Thus, negative colour gradients at older ages are likely driven by metallicity gradients. The diminishing age gradient may result from the starburst fading. Furthermore, the persistence of the [Fe/H] gradients may suggest that the outskirts are simultaneously built up by mergers with lower metallicity satellites. On the other hand, the gradients could be inherited from the star-forming phase, in which case mergers may not be needed to explain our findings. This work illustrates the need for resolved spectroscopy, instead of just photometry, to measure stellar population gradients. Extending these measurements to higher redshift is imperative for understanding how stellar populations in quiescent galaxies are assembled over cosmic time., Comment: Accepted for publication in MNRAS; minor typesetting corrections after copyediting

We present the galaxy stellar mass - size relation in the rest-frame near-IR ($1.5~\mu{\text{m}}$) and its evolution with redshift up to $z=2.5$. S\'ersic profiles are measured for $\sim$ $26\,000$ galaxies with stellar masses $M_\star > 10^9~{\text{M}}_\odot$ from JWST/NIRCam F277W and F444W imaging provided by the COSMOS-WEB and PRIMER surveys, using coordinates, redshifts, colors and stellar mass estimates from the COSMOS2020 catalog. The new rest-frame near-IR effective radii are generally smaller than previously measured rest-frame optical sizes, on average by 0.14~dex, with no significant dependence on redshift. For quiescent galaxies this size offset does not depend on stellar mass, but for star-forming galaxies the offset increases from -0.1~dex at $M_\star = 10^{9.5}~{\text{M}}_\odot$ to -0.25~dex at $M_\star > 10^{11}~{\text{M}}_\odot$. That is, we find that the near-IR stellar mass - size relation for star-forming galaxies is flatter in the rest-frame near-IR than in the rest-frame optical at all redshifts $0.510^{11}~{\text{M}}_\odot$) star-forming galaxies evolve in size almost as fast as quiescent galaxies. Low-mass ($M_\star<10^{10}~{\text{M}}_\odot$)~quiescent galaxies evolve as slow as star-forming galaxies. Our main conclusion is that the size evolution narrative as it has emerged over the past two decades does not radically change when accessing with JWST the rest-frame near-IR, a better proxy of the underlying stellar mass distribution., Comment: Accepted for publication. V2: includes a new table (table 1) of the data produced and available in machine-readable format on the journal version of the article

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\v{c}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., Comment: Main text 14 pages, 5 figures, accepted to MNRAS. Generated data products are publicly available at www.tng-project.org/gebek24

We aim to provide observational signatures of the dust size evolution in the ISM. In particular, we explore indicators of the polycyclic aromatic hydrocarbon (PAH) mass fraction ($q_{PAH}$), defined as the mass fraction of PAHs relative to total dust grains. In addition, we validate our dust evolution model by comparing the observational signatures from our simulations to observations. We used the hydrodynamic simulation code, GADGET4-OSAKA to model the dust properties of Milky Way-like and NGC 628-like galaxies representing star-forming galaxies. This code incorporates the evolution of grain size distributions driven by dust production and interstellar processing. Furthermore, we performed post-processing dust radiative transfer with SKIRT based on the simulations to predict the observational properties. We find that the intensity ratio between 8 um and 24 um is correlated with $q_{PAH}$ and can be used as an indicator of PAH mass fraction. However, this ratio is influenced by the radiation field. We suggest the 8 um-to-total infrared intensity ratio ($\nu I_\nu(8 \mu m)/I$(TIR)) as another indicator, since it is tightly correlated with $q_{PAH}$. Furthermore, we explored the spatially resolved $q_{PAH}$ in the simulated Milky Way-like galaxy using $\nu I_\nu(8 \mu m)/I$(TIR). We find that the spatially resolved $q_{PAH}$ increases with metallicity at metallicity at Z<0.2 Zsun due to the interplay between accretion and shattering while it decreases at Z>0.2 Zsun because of coagulation. Finally, we compared the above indicators in the NGC 628-like simulation with those observed in NGC 628 by recent observations. Consequently, we find that our simulation underestimates the PAH mass fraction throughout the entire galaxy by a factor of $\sim 8$ on average. This could be due to the efficient loss of PAHs by coagulation in our model., Comment: 16 pages, 12 figures, accepted by A&A

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...)., Comment: 11 pages, 10 figures, accepted for publication in A&A

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., Comment: 13 pages, 7 figures, accepted for publication in A&A

We test whether we can predict optical spectra from deep-field photometry of distant galaxies. Our goal is to perform a comparison in data space, highlighting the differences between predicted and observed spectra. The Large Early Galaxy Astrophysics Census (LEGA-C) provides high-quality optical spectra of thousands of galaxies at redshift $0.6

We present and publicly release a new star-forming regions emission library TODDLERS (Time evolution of Observables including Dust Diagnostics and Line Emission from Regions containing young Stars) for the publicly available radiative transfer code SKIRT. The library generation involves the spherical evolution of a homogeneous gas cloud around a young stellar cluster that accounts for stellar feedback processes including stellar winds, supernovae, and radiation pressure, as well as the gravitational forces on the gas. The semi-analytical evolution model is coupled with the photoionization code Cloudy to calculate time-dependent UV-mm spectral energy distributions (SEDs) from star-forming regions of varying metallicity, star-formation efficiency, birth-cloud density, and mass. The calculated SEDs include the stellar, nebular, and dust continuum emission along with a wide range of emission lines originating from H ii, photo-dissociation, and molecular gas regimes tabulated at high resolution. The SEDs incorporated in SKIRT are generated by calculating a stellar-mass normalized luminosity, which assumes that each emission source is composed of a power-law population of star-forming clouds. When compared to the previous treatment of star-forming regions in SKIRT, TODDLERS shows a better agreement with low-redshift observational data in the IR wavelength range while offering a more comprehensive line-emission support. This paves the way for a variety of applications using simulated galaxies at low and high redshift., Comment: Accepted for publication in MNRAS on 27/09/23

We examine the wavelength dependence of radial light profiles based on S\'ersic index $n$ measurements of 1067 galaxies with M$_*\geq$ 10$^{9.5}$M$_\odot$ and in the redshift range $0.5 < z < 3$. The sample and rest-frame optical light profiles are drawn from CANDELS$+$3D-HST; rest-frame near-infrared light profiles are inferred from CEERS JWST/NIRCam imaging. $n$ shows only weak dependence on wavelength, regardless of redshift, galaxy mass and type: on average, star-forming galaxies have $n = 1-1.5$ and quiescent galaxies have $n = 3-4$ in the rest-frame optical and near-infrared. The strong correlation at all wavelengths between $n$ and star-formation activity implies a physical connection between the radial stellar mass profile and star-formation activity. The main caveat is that the current sample is too small to discern trends for the most massive galaxies (M$_* > 10^{11}M_\odot$)., Comment: 15 pages, 7 figures, Submitted to ApJ. Comments welcome

We present individual star-formation histories of $\sim3000$ massive galaxies (log($\mathrm{M_*/M_{\odot}}$) > 10.5) from the Large Early Galaxy Astrophysics Census (LEGA-C) spectroscopic survey at a lookback time of $\sim$7 billion years and quantify the population trends leveraging 20hr-deep integrated spectra of these $\sim$ 1800 star-forming and $\sim$ 1200 quiescent galaxies at 0.6 < $z$ < 1.0. Essentially all galaxies at this epoch contain stars of age < 3 Gyr, in contrast with older massive galaxies today, facilitating better recovery of previous generations of star formation at cosmic noon and earlier. We conduct spectro-photometric analysis using parametric and non-parametric Bayesian SPS modeling tools - Bagpipes and Prospector to constrain the median star-formation histories of this mass-complete sample and characterize population trends. A consistent picture arises for the late-time stellar mass growth when quantified as $t_{50}$ and $t_{90}$, corresponding to the age of the universe when galaxies formed 50\% and 90\% of their total stellar mass, although the two sets of models disagree at the earliest formation times (e.g. $t_{10}$). Our results reveal trends in both stellar mass and stellar velocity dispersion as in the local universe - low-mass galaxies with shallower potential wells grow their stellar masses later in cosmic history compared to high-mass galaxies. Unlike local quiescent galaxies, the median duration of late-time star-formation ($\tau_{SF,late}$ = $t_{90}$ - $t_{50}$) does not consistently depend on the stellar mass. This census sets a benchmark for future deep spectro-photometric studies of the more distant universe., Comment: 28 pages, 12 figures, 4 tables, Accepted to ApJ

The massive end of the gas-phase mass--metallicity relation (MZR) is a sensitive probe of active galactic nuclei (AGN) feedback that is a crucial but highly uncertain component of galaxy evolution models. In this paper, we extend the $z\sim0.7$ MZR by $\sim$0.5 dex up to log$(M_\star/\textrm{M}_\odot)\sim11.1$. We use extremely deep VLT VIMOS spectra from the Large Early Galaxy Astrophysics Census (LEGA-C) survey to measure metallicities for 145 galaxies. The LEGA-C MZR matches the normalization of the $z\sim0.8$ DEEP2 MZR where they overlap, so we combine the two to create an MZR spanning from 9.3 to 11.1 log$(M_\star/\textrm{M}_\odot)$. The LEGA-C+DEEP2 MZR at $z\sim0.7$ is offset to slightly lower metallicities (0.05-0.13 dex) than the $z\sim0$ MZR, but it otherwise mirrors the established power law rise at low/intermediate stellar masses and asymptotic flattening at high stellar masses. We compare the LEGA-C+DEEP2 MZR to the MZR from two cosmological simulations (IllustrisTNG and SIMBA), which predict qualitatively different metallicity trends for high-mass galaxies. This comparison highlights that our extended MZR provides a crucial observational constraint for galaxy evolution models in a mass regime where the MZR is very sensitive to choices about the implementation of AGN feedback., Comment: 10 pages, 4 figures, 1 table

We present the first study of spatially integrated higher-order stellar kinematics over cosmic time. We use deep rest-frame optical spectroscopy of quiescent galaxies at redshifts z=0.05, 0.3 and 0.8 from the SAMI, MAGPI and LEGA-C surveys to measure the excess kurtosis $h_4$ of the stellar velocity distribution, the latter parametrised as a Gauss-Hermite series. Conservatively using a redshift-independent cut in stellar mass ($M_\star = 10^{11}\,{\rm M}_\odot$), and matching the stellar-mass distributions of our samples, we find 7 $\sigma$ evidence of $h_4$ increasing with cosmic time, from a median value of 0.019$\pm$0.002 at z=0.8 to 0.059$\pm$0.004 at z=0.06. Alternatively, we use a physically motivated sample selection, based on the mass distribution of the progenitors of local quiescent galaxies as inferred from numerical simulations; in this case, we find 10 $\sigma$ evidence. This evolution suggests that, over the last 7 Gyr, there has been a gradual decrease in the rotation-to-dispersion ratio and an increase in the radial anisotropy of the stellar velocity distribution, qualitatively consistent with accretion of gas-poor satellites. These findings demonstrate that massive galaxies continue to accrete mass and increase their dispersion support after becoming quiescent., Comment: 19 pages, 9 figures Accepted for publication in MNRAS

We present the first statistical study of spatially integrated non-Gaussian stellar kinematics spanning 7 Gyr in cosmic time. We use deep, rest-frame optical spectroscopy of massive galaxies (stellar mass $M_\star > 10^{10.5} {\rm M}_\odot$) at redshifts z = 0.05, 0.3 and 0.8 from the SAMI, MAGPI and LEGA-C surveys, to measure the excess kurtosis $h_4$ of the stellar velocity distribution, the latter parametrised as a Gauss-Hermite series. We find that at all redshifts where we have large enough samples, $h_4$ anti-correlates with the ratio between rotation and dispersion, highlighting the physical connection between these two kinematic observables. In addition, and independently from the anti-correlation with rotation-to-dispersion ratio, we also find a correlation between $h_4$ and $M_\star$, potentially connected to the assembly history of galaxies. In contrast, after controlling for mass, we find no evidence of independent correlation between $h_4$ and aperture velocity dispersion or galaxy size. These results hold for both star-forming and quiescent galaxies. For quiescent galaxies, $h_4$ also correlates with projected shape, even after controlling for the rotation-to-dispersion ratio. At any given redshift, star-forming galaxies have lower $h_4$ compared to quiescent galaxies, highlighting the link between kinematic structure and star-forming activity., Comment: 26 pages, 15 figures Accepted for publication in MNRAS

We present a nonparametric morphology analysis of the stellar continuum and nebular emission lines for a sample of local galaxies. We explore the dependence of the various morphological parameters on wavelength and morphological type. Our goal is to quantify the difference in morphology between the stellar and nebular components. We derive the nonparametric morphological indicators of 364 galaxies from the CALIFA Survey. To calculate those indicators, we apply the StatMorph package on the high-quality integral field spectroscopic data cubes, as well as to the most prominent nebular emission-line maps, namely [OIII]$\lambda$5007, H$\alpha$, and [NII]$\lambda$6583. We show that the physical size of galaxies, M$_{20}$ index, and concentration have a strong gradient from blue to red optical wavelengths. We find that the light distribution of the nebular emission is less concentrated than the stellar continuum. A comparison between the nonparametric indicators and the galaxy physical properties revealed a very strong correlation of the concentration with the specific star-formation rate and morphological type. Furthermore, we explore how the galaxy inclination affects our results. We find that edge-on galaxies show a more rapid change in physical size and concentration with increasing wavelength due to the increase in optical free path. We conclude that the apparent morphology of galaxies originates from the pure stellar distribution, but the morphology of the ISM presents differences with respect to the morphology of the stellar component. Our analysis also highlights the importance of dust attenuation and galaxy inclination in the measurement of nonparametric morphological indicators, especially in the the wavelength range 4000-5000 \r{A}.

Rejuvenating galaxies are unusual galaxies that fully quench and then subsequently experience a "rejuvenation" event to become star-forming once more. Rejuvenation rates vary substantially in models of galaxy formation: 10%-70% of massive galaxies are expected to experience rejuvenation by z = 0. Measuring the rate of rejuvenation is therefore important for calibrating the strength of star formation feedback mechanisms. However, these observations are challenging because rejuvenating systems blend in with normal star-forming galaxies in broadband photometry. In this paper, we use the galaxy spectral energy distribution (SED)-fitting code Prospector to search for observational markers that distinguish normal star-forming galaxies from rejuvenating galaxies. We find that rejuvenating galaxies have smaller Balmer absorption line equivalent widths (EWs) than star-forming galaxies. This is analogous to the well-known "K + A" or post-starburst galaxies, which have strong Balmer absorption due to A-stars dominating the light: in this case, rejuvenating systems have a lack of A-stars, instead resembling "O - A" systems. We find star-forming galaxies that have H$\beta$, H$\gamma$, and/or H$\delta$ absorption EWs $\lesssim 3${\AA} corresponds to a highly pure selection of rejuvenating systems. Interestingly, while this technique is highly effective at identifying mild rejuvenation, "strongly" rejuvenating systems remain nearly indistinguishable from star-forming galaxies due to the well-known stellar outshining effect. We conclude that measuring Balmer absorption line EWs in star-forming galaxy populations is an efficient method to identify rejuvenating populations, and discuss several techniques to either remove or resolve the nebular emission which typically lies on top of these absorption lines., Comment: Published in ApJ, 16 pages, 9 figures

We perform a homogeneous analysis of an unprecedented set of spatially resolved scaling relations (SRs) between ISM components and other properties in the range of scales 0.3-3.4 kpc. We also study some ratios: dust-to-stellar, dust-to-gas, and dust-to-metal. We use a sample of 18 large, spiral, face-on DustPedia galaxies. All the SRs are moderate/strong correlations except the dust-HI SR that does not exist or is weak for most galaxies. The SRs do not have a universal form but each galaxy is characterized by distinct correlations, affected by local processes and galaxy peculiarities. The SRs hold starting from 0.3 kpc, and if a breaking down scale exists it is < 0.3 kpc. By evaluating all galaxies at 3.4 kpc, differences due to peculiarities of individual galaxies are cancelled out and the corresponding SRs are consistent with those of whole galaxies. By comparing subgalactic and global scales, the most striking result emerges from the SRs involving ISM components: the dust-total gas SR is a good correlation at all scales, while the dust-H2 and dust-HI SRs are good correlations at subkpc/kpc and total scales, respectively. For the other explored SRs, there is a good agreement between small and global scales and this may support the picture where the main physical processes regulating the properties and evolution of galaxies occur locally. Our results are consistent with the hypothesis of self-regulation of the SF process. The analysis of subgalactic ratios shows that they are consistent with those derived for whole galaxies, from low to high z, supporting the idea that also these ratios could be set by local processes. Our results highlight the heterogeneity of galaxy properties and the importance of resolved studies on local galaxies in the context of galaxy evolution. They also provide observational constraints to theoretical models and updated references for high-z studies., Comment: 42 pages, 11 figures and 5 tables in the main text, 2 figures and 1 table in Appendix. Accepted for publication in A&A

Cosmic dust grains are one of the fundamental ingredients of the interstellar medium (ISM). In spite of their small contribution to the total mass budget, dust grains play a significant role in the physical and chemical evolution of galaxies. Over the past decades, a plethora of multi-wavelength data, from UV to far-infrared, has increased substantially our knowledge on the dust properties of nearby galaxies. Nevertheless, one regime of the spectrum, the mm range, remains relatively unexplored. Thanks to the new, high-resolution data in the mm range observed with the NIKA2 instrument and our radiative transfer framework, we aim to firmly characterise the physical properties of the very cold dust (<15K), and to quantify the importance of different emission mechanisms in the mm. So far, we have developed a methodology to use dust radiative transfer modelling and applied it to a small group of face-on spiral galaxies. The combination of the new NIKA2 data with our radiative transfer techniques would provide the right conditions to generate an accurate model of the interplay between starlight and dust in a sizeable sample of spatially-resolved nearby galaxies., Comment: To appear in the Proceedings of the International Conference entitled "mm Universe @ NIKA2", Rome (Italy), June 2021, EPJ Web of conferences

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.

Methods. We have modelled a sample of ~800 nearby galaxies, spanning a wide range of metallicity, gas fraction, specific star formation rate and Hubble stage. We have derived the dust properties of each object from its spectral energy distribution. Through an additional level of analysis, we have inferred the timescales of dust condensation in core-collapse supernova ejecta, grain growth in cold clouds and dust destruction by shock waves. Throughout this paper, we have adopted a hierarchical Bayesian approach, resulting in a single large probability distribution of all the parameters of all the galaxies, to ensure the most rigorous interpretation of our data. Results. We confirm the drastic evolution with metallicity of the dust-to-metal mass ratio (by two orders of magnitude), found by previous studies. We show that dust production by core-collapse supernovae is efficient only at very low-metallicity, a single supernova producing on average less than ~0.03 Msun/SN of dust. Our data indicate that grain growth is the dominant formation mechanism at metallicity above ~1/5 solar, with a grain growth timescale shorter than ~50 Myr at solar metallicity. Shock destruction is relatively efficient, a single supernova clearing dust on average in at least ~1200 Msun/SN of gas. These results are robust when assuming different stellar initial mass functions. In addition, we show that early-type galaxies are outliers in several scaling relations. This feature could result from grain thermal sputtering in hot X-ray emitting gas, an hypothesis supported by a negative correlation between the dust-to-stellar mass ratio and the X-ray photon rate per grain. Finally, we confirm the well-known evolution of the aromatic-feature-emitting grain mass fraction as a function of metallicity and interstellar radiation field intensity. Our data indicate the relation with metallicity is significantly stronger., Comment: 64 pages, 36 figures, accepted for publication in Astronomy & Astrophysics

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.

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., Comment: 9 pages, 5 figures, accepted for publication in A&A

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), Comment: 23 pages, 16 figures, accepted for publication in A&A

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.

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., Comment: 17 pages, 14 figures, accepted for publication in MNRAS

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., Comment: 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

Understanding how dust attenuation laws vary between and within galaxies is a key question if we want to reliably measure the physical properties of galaxies at both global and local scales. To shed new light on this question, we present a detailed study of the slope and bump strength of the attenuation law in the nearby spiral galaxy NGC628 at the resolved spatial scale of 325 pc. To do so, we have modelled a broad multi-wavelength dataset from the ultraviolet (UV) to the infrared (IR) with the state-of-the-art SED fitting code CIGALE, including SWIFT UVOT data for which we have developed a new optimized reduction pipeline. We find that the median dust attenuation curve of NGC628 is fairly steep, but not as steep as the SMC curve, and has a sub-MW-type UV bump. We observe intriguing variations within the galaxy, with regions of high $A_V$ exhibiting a shallower attenuation curve. We argue that the flattening of the curve is due to a dominance of absorption-over-scattering events at higher $A_V$. No trend between the bump strength and the IRAC 8.0 micron emission was found. However, this does not necessarily rule out PAHs as the main contributors to the UV bump., Comment: 26 pages, 16 figures, 4 tables, accepted for publication in MNRAS

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., Comment: 14 pages, 5 figures, accepted for publication in MNRAS