Showing 1–46 of 46 results for author: Nersesian, A

The Size-Mass relation at Rest-Frame $1.5μ$m from JWST/NIRCam in the COSMOS-WEB and PRIMER-COSMOS fields

Authors: Marco Martorano, Arjen van der Wel, Maarten Baes, Eric F. Bell, Gabriel Brammer, Marijn Franx, Angelos Nersesian

Abstract: We present the galaxy stellar mass - size relation in the rest-frame near-IR ($1.5~μ{\text{m}}$) and its evolution with redshift up to $z=2.5$. Sérsic 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 e… ▽ More We present the galaxy stellar mass - size relation in the rest-frame near-IR ($1.5~μ{\text{m}}$) and its evolution with redshift up to $z=2.5$. Sérsic 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.5<z<2.5$. The general pace of size evolution is the same in the near-IR as previously demonstrated in the optical, with slower evolution ($R_{\text{e}} \propto (1+z)^{-0.7}$) for $L^*$~star-forming galaxies and faster evolution ($R_{\text{e}} \propto (1+z)^{-1.3}$) for $L^*$~quiescent galaxies. Massive ($M_\star>10^{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. △ Less

Submitted 25 June, 2024; originally announced June 2024.

Comments: Accepted for publication

Euclid: Early Release Observations -- Overview of the Perseus cluster and analysis of its luminosity and stellar mass functions

Authors: J. -C. Cuillandre, M. Bolzonella, A. Boselli, F. R. Marleau, M. Mondelin, J. G. Sorce, C. Stone, F. Buitrago, Michele Cantiello, K. George, N. A. Hatch, L. Quilley, F. Mannucci, T. Saifollahi, R. Sánchez-Janssen, F. Tarsitano, C. Tortora, X. Xu, H. Bouy, S. Gwyn, M. Kluge, A. Lançon, R. Laureijs, M. Schirmer, Abdurro'uf , et al. (177 additional authors not shown)

Abstract: The Euclid ERO programme targeted the Perseus cluster of galaxies, gathering deep data in the central region of the cluster over 0.7 square degree, corresponding to approximately 0.25 r_200. The data set reaches a point-source depth of IE=28.0 (YE, JE, HE = 25.3) AB magnitudes at 5 sigma with a 0.16" and 0.48" FWHM, and a surface brightness limit of 30.1 (29.2) mag per square arcsec. The exception… ▽ More The Euclid ERO programme targeted the Perseus cluster of galaxies, gathering deep data in the central region of the cluster over 0.7 square degree, corresponding to approximately 0.25 r_200. The data set reaches a point-source depth of IE=28.0 (YE, JE, HE = 25.3) AB magnitudes at 5 sigma with a 0.16" and 0.48" FWHM, and a surface brightness limit of 30.1 (29.2) mag per square arcsec. The exceptional depth and spatial resolution of this wide-field multi-band data enable the simultaneous detection and characterisation of both bright and low surface brightness galaxies, along with their globular cluster systems, from the optical to the NIR. This study advances beyond previous analyses of the cluster and enables a range of scientific investigations summarised here. We derive the luminosity and stellar mass functions (LF and SMF) of the Perseus cluster in the Euclid IE band, thanks to supplementary u,g,r,i,z and Halpha data from the CFHT. We adopt a catalogue of 1100 dwarf galaxies, detailed in the corresponding ERO paper. We identify all other sources in the Euclid images and obtain accurate photometric measurements using AutoProf or AstroPhot for 138 bright cluster galaxies, and SourceExtractor for half a million compact sources. Cluster membership for the bright sample is determined by calculating photometric redshifts with Phosphoros. Our LF and SMF are the deepest recorded for the Perseus cluster, highlighting the groundbreaking capabilities of the Euclid telescope. Both the LF and SMF fit a Schechter plus Gaussian model. The LF features a dip at M(IE)=-19 and a faint-end slope of alpha_S = -1.2 to -1.3. The SMF displays a low-mass-end slope of alpha_S = -1.2 to -1.35. These observed slopes are flatter than those predicted for dark matter halos in cosmological simulations, offering significant insights for models of galaxy formation and evolution. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Submitted to A&A, 44 pages, 35 figures, Part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

Euclid. I. Overview of the Euclid mission

Authors: Euclid Collaboration, Y. Mellier, Abdurro'uf, J. A. Acevedo Barroso, A. Achúcarro, J. Adamek, R. Adam, G. E. Addison, N. Aghanim, M. Aguena, V. Ajani, Y. Akrami, A. Al-Bahlawan, A. Alavi, I. S. Albuquerque, G. Alestas, G. Alguero, A. Allaoui, S. W. Allen, V. Allevato, A. V. Alonso-Tetilla, B. Altieri, A. Alvarez-Candal, A. Amara, L. Amendola , et al. (1086 additional authors not shown)

Abstract: The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14… ▽ More The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue`Euclid on Sky'

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

Observational signatures of the dust size evolution in isolated galaxy simulations

Authors: Kosei Matsumoto, Hiroyuki Hirashita, Kentaro Nagamine, Stefan van der Giessen, Leonard E. C. Romano, Monica Relaño, Ilse De Looze, Maarten Baes, Angelos Nersesian, Peter Camps, Kuan-chou Hou, Yuri Oku

Abstract: We aim to provide observational signatures of the dust size evolution in the ISM, in particular exploring indicators of polycyclic aromatic hydrocarbon (PAH) mass fraction ($q_{PAH}$) defined as the mass fraction of PAHs relative to total dust grains. Additionally, we validate our dust evolution model by comparing the observational signatures from our simulations to those from observations. We mod… ▽ More We aim to provide observational signatures of the dust size evolution in the ISM, in particular exploring indicators of polycyclic aromatic hydrocarbon (PAH) mass fraction ($q_{PAH}$) defined as the mass fraction of PAHs relative to total dust grains. Additionally, we validate our dust evolution model by comparing the observational signatures from our simulations to those from observations. We model the evolution of grain size distribution of Milky Way-like and NGC 628-like galaxies representing star-forming galaxies with a hydrodynamic simulation code, GADGET4-OSAKA, which considers dust production and interstellar processing. Furthermore, we perform 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 correlates with $q_{PAH}$ and can be used as an indicator of PAH mass fraction. However, this ratio is influenced by the radiation field. As another indicator, we suggest the 8 um-to-total infrared intensity ratio ($νI_ν(8 μm)/I$(TIR)), which tightly correlates with $q_{PAH}$. Furthermore, we explore the spatial evolution of $q_{PAH}$ in the simulated Milky Way-like galaxy using $νI_ν(8 μm)/I$(TIR). We find that the spatially resolved $q_{PAH}$ increases with metallicity at lower metallicity (Z<0.2 Zsun) due to the interplay between accretion and shattering while it decreases with metallicity at higher metallicity (Z>0.2 Zsun) due to coagulation. Finally, we compare the above indicators in the NGC 628-like simulation with those observed in NGC 628. Consequently, our simulation underestimates the PAH mass fraction throughout the entire galaxy. This is probably because PAH is too efficiently lost by coagulation in the interstellar medium in our model, which suggests that the inhibition of coagulation of the PAHs is key to enhancing PAH formation. △ Less

Submitted 4 February, 2024; originally announced February 2024.

Comments: 16 pages, 12 figures, submit to A&A

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

Less is less: photometry alone cannot predict the observed spectral indices of $z\sim1$ galaxies from the LEGA-C spectroscopic survey

Authors: Angelos Nersesian, Arjen van der Wel, Anna Gallazzi, Joel Leja, Rachel Bezanson, Eric F. Bell, Francesco D'Eugenio, Anna de Graaff, Yasha Kaushal, Marco Martorano, Michael Maseda, Stefano Zibetti

Abstract: 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<z<1$. Broad-band photometry of the same galaxies, dra… ▽ More 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<z<1$. Broad-band photometry of the same galaxies, drawn from the recent COSMOS2020 catalog, is used to predict the optical spectra with the spectral energy distribution (SED) fitting code Prospector and the MILES stellar library. The observed and predicted spectra are compared in terms of two age and metallicity-sensitive absorption features (H$δ_\mathrm{A}$ and Fe4383). The global bimodality of star-forming and quiescent galaxies in photometric space is recovered with the model spectra. But the presence of a systematic offset in the Fe4383 line strength and the weak correlation between the observed and modeled line strength imply that accurate age or metallicity determinations cannot be inferred from photometry alone. For now we caution that photometry-based estimates of stellar population properties are determined mostly by the modeling approach and not the physical properties of galaxies, even when using the highest-quality photometric datasets and state-of-the-art fitting techniques. When exploring a new physical parameter space (i.e. redshift or galaxy mass) high-quality spectroscopy is always needed to inform the analysis of photometry. △ Less

Submitted 27 October, 2023; originally announced October 2023.

Comments: 13 pages, 8 figures, accepted 26 October 2023

IAS/CEA Evolution of Dust in Nearby Galaxies (ICED): the spatially-resolved dust properties of NGC4254

Authors: L. Pantoni, R. Adam, P. Ade, H. Ajeddig, P. André, E. Artis, H. Aussel, M. Baes, A. Beelen, A. Benoît, S. Berta, L. Bing, O. Bourrion, M. Calvo, A. Catalano, M. De Petris, F. -X. Désert, S. Doyle, E. F. C. Driessen, G. Ejlali, F. Galliano, A. Gomez, J. Goupy, A. P. Jones, C. Hanser , et al. (35 additional authors not shown)

Abstract: We present the first preliminary results of the project \textit{ICED}, focusing on the face-on galaxy NGC4254. We use the millimetre maps observed with NIKA2 at IRAM-30m, as part of the IMEGIN Guaranteed Time Large Program, and of a wide collection of ancillary data (multi-wavelength photometry and gas phase spectral lines) that are publicly available. We derive the global and local properties of… ▽ More We present the first preliminary results of the project \textit{ICED}, focusing on the face-on galaxy NGC4254. We use the millimetre maps observed with NIKA2 at IRAM-30m, as part of the IMEGIN Guaranteed Time Large Program, and of a wide collection of ancillary data (multi-wavelength photometry and gas phase spectral lines) that are publicly available. We derive the global and local properties of interstellar dust grains through infrared-to-radio spectral energy distribution fitting, using the hierarchical Bayesian code HerBIE, which includes the grain properties of the state-of-the-art dust model, THEMIS. Our method allows us to get the following dust parameters: dust mass, average interstellar radiation field, and fraction of small grains. Also, it is effective in retrieving the intrinsic correlations between dust parameters and interstellar medium properties. We find an evident anti-correlation between the interstellar radiation field and the fraction of small grains in the centre of NGC4254, meaning that, at strong radiation field intensities, very small amorphous carbon grains are efficiently destroyed by the ultra-violet photons coming from newly formed stars, through photo-desorption and sublimation. We observe a flattening of the anti-correlation at larger radial distances, which may be driven by the steep metallicity gradient measured in NGC4254. △ Less

Submitted 10 October, 2023; originally announced October 2023.

Comments: to appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

Journal ref: EPJ Web of Conferences 293 (2024) 00038

Exploring the interstellar medium of NGC 891 at millimeter wavelengths using the NIKA2 camera

Authors: S. Katsioli, R. Adam, P. Ade, H. Ajeddig, P. André, E. Artis, H. Aussel, M. Baes, A. Beelen, A. Benoît, S. Berta, L. Bing, O. Bourrion, M. Calvo, A. Catalano, C. J. R. Clark, I. De Looze, M. De Petris, F. -X. Désert, S. Doyle, E. F. C. Driessen, G. Ejlali, M. Galametz, F. Galliano, A. Gomez , et al. (39 additional authors not shown)

Abstract: In the framework of the IMEGIN Large Program, we used the NIKA2 camera on the IRAM 30-m telescope to observe the edge-on galaxy NGC 891 at 1.15 mm and 2 mm and at a FWHM of 11.1" and 17.6", respectively. Multiwavelength data enriched with the new NIKA2 observations fitted by the HerBIE SED code (coupled with the THEMIS dust model) were used to constrain the physical properties of the ISM. Emission… ▽ More In the framework of the IMEGIN Large Program, we used the NIKA2 camera on the IRAM 30-m telescope to observe the edge-on galaxy NGC 891 at 1.15 mm and 2 mm and at a FWHM of 11.1" and 17.6", respectively. Multiwavelength data enriched with the new NIKA2 observations fitted by the HerBIE SED code (coupled with the THEMIS dust model) were used to constrain the physical properties of the ISM. Emission originating from the diffuse dust disk is detected at all wavelengths from mid-IR to mm, while mid-IR observations reveal warm dust emission from compact HII regions. Indications of mm excess emission have also been found in the outer parts of the galactic disk. Furthermore, our SED fitting analysis constrained the mass fraction of the small (< 15 Angstrom) dust grains. We found that small grains constitute 9.5% of the total dust mass in the galactic plane, but this fraction increases up to ~ 20% at large distances (|z| > 3 kpc) from the galactic plane. △ Less

Submitted 6 October, 2023; originally announced October 2023.

Comments: To appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

Journal ref: EPJ Web of conferences 293 (2024) 00026

Constraining Millimeter Dust Emission in Nearby Galaxies with NIKA2: the case of NGC2146 and NGC2976

Authors: G. Ejlali, R. Adam, P. Ade, H. Ajeddig, P. André, E. Artis, H. Aussel, M. Baes, A. Beelen, Benoît, S. Berta, L. Bing, O. Bourrion, M. Calvo, A. Catalano, M. De Petris, F. -X. Désert, S. Doyle, E. F. C. Driessen, F. Galliano, A. Gomez, J. Goupy, A. P. Jones, C. Hanser, A. Hughes , et al. (35 additional authors not shown)

Abstract: This study presents the first millimeter continuum mapping observations of two nearby galaxies, the starburst spiral galaxy NGC2146 and the dwarf galaxy NGC2976, at 1.15 mm and 2 mm using the NIKA2 camera on the IRAM 30m telescope, as part of the Guaranteed Time Large Project IMEGIN. These observations provide robust resolved information about the physical properties of dust in nearby galaxies by… ▽ More This study presents the first millimeter continuum mapping observations of two nearby galaxies, the starburst spiral galaxy NGC2146 and the dwarf galaxy NGC2976, at 1.15 mm and 2 mm using the NIKA2 camera on the IRAM 30m telescope, as part of the Guaranteed Time Large Project IMEGIN. These observations provide robust resolved information about the physical properties of dust in nearby galaxies by constraining their FIR-radio SED in the millimeter domain. After subtracting the contribution from the CO line emission, the SEDs are modeled spatially using a Bayesian approach. Maps of dust mass surface density, temperature, emissivity index, and thermal radio component of the galaxies are presented, allowing for a study of the relations between the dust properties and star formation activity (using observations at 24$μ$m as a tracer). We report that dust temperature is correlated with star formation rate in both galaxies. The effect of star formation activity on dust temperature is stronger in NGC2976, an indication of the thinner interstellar medium of dwarf galaxies. Moreover, an anti-correlation trend is reported between the dust emissivity index and temperature in both galaxies. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: To appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

Journal ref: EPJ Web of conferences 293 (2024) 00016

TODDLERS: A new UV-mm emission library for star-forming regions. 1. Integration with SKIRT and public release

Authors: Anand Utsav Kapoor, Maarten Baes, Arjen van der Wel, Andrea Gebek, Peter Camps, Angelos Nersesian, Sharon E. Meidt, Aaron Smith, Sebastien Vicens, Francesco D'Eugenio, Marco Martorano, Daniela Barrientos, Nina Sanches Sartorio

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

Submitted 20 October, 2023; v1 submitted 30 September, 2023; originally announced October 2023.

Comments: Accepted for publication in MNRAS on 27/09/23

The stratification of ISM properties in the edge-on galaxy NGC 891 revealed by NIKA2

Authors: S. Katsioli, E. M. Xilouris, C. Kramer, R. Adam, P. Ade, H. Ajeddig, P. André, E. Artis, H. Aussel, M. Baes, A. Beelen, A. Benoît, S. Berta, L. Bing, O. Bourrion, M. Calvo, A. Catalano, C. J. R. Clark, I. De Looze, M. De Petris, F. -X. Désert, S. Doyle, E. F. C. Driessen, G. Ejlali, M. Galametz , et al. (38 additional authors not shown)

Abstract: As the millimeter wavelength range remains a largely unexplored spectral region for galaxies, the IMEGIN large program aims to map the millimeter continuum emission of 22 nearby galaxies at 1.15 and 2 mm. Using the high-resolution maps produced by the NIKA2 camera, we explore the existence of very cold dust and take possible contamination by free-free and synchrotron emission into account. We stud… ▽ More As the millimeter wavelength range remains a largely unexplored spectral region for galaxies, the IMEGIN large program aims to map the millimeter continuum emission of 22 nearby galaxies at 1.15 and 2 mm. Using the high-resolution maps produced by the NIKA2 camera, we explore the existence of very cold dust and take possible contamination by free-free and synchrotron emission into account. We study the IR-to-radio emission coming from different regions along the galactic plane and at large vertical distances. New observations of NGC 891, using the NIKA2 camera on the IRAM 30m telescope, along with a suite of observations at other wavelengths were used to perform a multiwavelength study of the spectral energy distribution in the interstellar medium in this galaxy. This analysis was performed globally and locally, using the advanced hierarchical Bayesian fitting code, HerBIE, coupled with the THEMIS dust model. Our dust modeling is able to reproduce the near-IR to millimeter emission of NGC 891, with the exception of an excess at a level of 25% obtained by the NIKA2 observations in the outermost parts of the disk. The radio continuum and thermal dust emission are distributed differently in the disk and galaxy halo. Different dusty environments are also revealed by a multiwavelength investigation of the emission features. Our detailed decomposition at millimeter and centimeter wavelengths shows that emission at 1 mm is purely originated by dust. Radio components become progressively important with increasing wavelengths. Finally, we find that emission arising from small dust grains accounts for ~ 9.5% of the total dust mass, reaching up to 20% at large galactic latitudes. Shock waves in the outflows that shatter the dust grains might explain this higher fraction of small grains in the halo. △ Less

Submitted 15 September, 2023; originally announced September 2023.

Comments: 15 pages, 10 figures, Accepted for publication in A&A

Journal ref: A&A 679, A7 (2023)

Rest-Frame Near-Infrared Radial Light Profiles up to z=3 from JWST/NIRCam: Wavelength Dependence of the Sérsic Index

Authors: Marco Martorano, Arjen van der Wel, Eric F. Bell, Marijn Franx, Katherine E. Whitaker, Angelos Nersesian, Sedona H. Price, Maarten Baes, Katherine A. Suess, Erica J. Nelson, Tim B. Miller, Rachel Bezanson, Gabriel Brammer

Abstract: We examine the wavelength dependence of radial light profiles based on Sérsic 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 wavel… ▽ More We examine the wavelength dependence of radial light profiles based on Sérsic 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$). △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: 15 pages, 7 figures, Submitted to ApJ. Comments welcome

Journal ref: ApJ 957 46 (2023)

A census of star formation histories of massive galaxies at 0.6 < z < 1 from spectro-photometric modeling using Bagpipes and Prospector

Authors: Yasha Kaushal, Angelos Nersesian, Rachel Bezanson, Arjen van der Wel, Joel Leja, Adam Carnall, Stefano Zibetti, Gourav Khullar, Marijn Franx, Adam Muzzin, Anna De Graaff, Camilla Pacifici, Katherine E. Whitaker, Eric F. Bell, Marco Martorano

Abstract: 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.… ▽ More 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 ($τ_{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. △ Less

Submitted 9 November, 2023; v1 submitted 7 July, 2023; originally announced July 2023.

Comments: 28 pages, 12 figures, 4 tables, Accepted to ApJ

The Gas-Phase Mass--Metallicity Relation for Massive Galaxies at $z\sim0.7$ with the LEGA-C Survey

Authors: Zach J. Lewis, Brett H. Andrews, Rachel Bezanson, Michael Maseda, Eric F. Bell, Romeel Davé, Francesco D'Eugenio, Marijn Franx, Anna Gallazzi, Anna de Graaff, Yasha Kaushal, Angelos Nersesian, Jeffrey A. Newman, Arjen van der Wel, Po-Feng Wu

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

Submitted 24 April, 2023; originally announced April 2023.

Comments: 10 pages, 4 figures, 1 table

Evolution in the orbital structure of quiescent galaxies from MAGPI, LEGA-C and SAMI surveys: direct evidence for merger-driven growth over the last 7 Gy

Authors: Francesco D'Eugenio, Arjen van der Wel, Joanna M. Piotrowska, Rachel Bezanson, Edward N. Taylor, Jesse van de Sande, William M. Baker, Eric F. Bell, Sabine Bellstedt, Joss Bland-Hawthorn, Asa F. L. Bluck, Sarah Brough, Julia J. Bryant, Matthew Colless, Luca Cortese, Scott M. Croom, Caro Derkenne, Pieter van Dokkum, Deanne Fisher, Caroline Foster, Anna Gallazzi, Anna de Graaff, Brent Groves, Josha van Houdt, Claudia del P. Lagos , et al. (15 additional authors not shown)

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

Submitted 9 March, 2023; originally announced March 2023.

Comments: 19 pages, 9 figures Accepted for publication in MNRAS

Different higher-order kinematics between star-forming and quiescent galaxies based on the SAMI, MAGPI and LEGA-C surveys

Authors: Francesco D'Eugenio, Arjen van der Wel, Caro Derkenne, Josha van Houdt, Rachel Bezanson, Edward N. Taylor, Jesse van de Sande, William M. Baker, Eric F. Bell, Joss Bland-Hawthorn, Asa F. L. Bluck, Sarah Brough, Julia J. Bryant, Matthew Colless, Luca Cortese, Scott M. Croom, Pieter van Dokkum, Deanne Fisher, Caroline Foster, Amelia Fraser-McKelvie, Anna Gallazzi, Anna de Graaff, Brent Groves, Claudia del P. Lagos, Tobias J. Looser , et al. (16 additional authors not shown)

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

Submitted 7 March, 2023; originally announced March 2023.

Comments: 26 pages, 15 figures Accepted for publication in MNRAS

Nonparametric galaxy morphology from stellar and nebular emission with the CALIFA sample

Authors: Angelos Nersesian, Stefano Zibetti, Francesco D'Eugenio, Maarten Baes

Abstract: 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… ▽ More 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]$λ$5007, H$α$, and [NII]$λ$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 Å. △ Less

Submitted 3 March, 2023; originally announced March 2023.

Journal ref: A&A 673, A63 (2023)

A simple spectroscopic technique to identify rejuvenating galaxies

Authors: Junyu Zhang, Yijia Li, Joel Leja, Katherine E. Whitaker, Angelos Nersesian, Rachel Bezanson, Arjen van der Wel

Abstract: 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 formati… ▽ More 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$β$, H$γ$, and/or H$δ$ absorption EWs $\lesssim 3$Å 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. △ Less

Submitted 13 July, 2023; v1 submitted 18 November, 2022; originally announced November 2022.

Comments: Published in ApJ, 16 pages, 9 figures

The resolved scaling relations in DustPedia: Zooming in on the local Universe

Authors: Viviana Casasola, Simone Bianchi, Laura Magrini, Aleksandr V. Mosenkov, Francesco Salvestrini, Maarten Baes, Francesco Calura, Letizia P. Cassara', Christopher J. R. Clark, Edvige Corbelli, Jacopo Fritz, Frederic Galliano, Elisabetta Liuzzo, Suzanne Madden, Angelos Nersesian, Francesca Pozzi, Sambit Roychowdhury, Ivano Baronchelli, Matteo Bonato, Carlotta Gruppioni, Lara Pantoni

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

Submitted 28 October, 2022; originally announced October 2022.

Comments: 42 pages, 11 figures and 5 tables in the main text, 2 figures and 1 table in Appendix. Accepted for publication in A&A

Journal ref: A&A 668, A130 (2022)

Star-forming early- and quiescent late-type galaxies in the local Universe

Authors: E. -D. Paspaliaris, E. M. Xilouris, A. Nersesian, S. Bianchi, I. Georgantopoulos, V. A. Masoura, G. E. Magdis, M. Plionis

Abstract: The general consensus is that LTGs undergo intense star-formation activity, while ETGs are mostly inactive. We question this general rule and investigate the existence of star-forming ETGs and quiescent LTGs in the local Universe. By computing the physical properties of 2,209 such galaxies in the GAMA survey being morphologically classified and using information on their structural properties as w… ▽ More The general consensus is that LTGs undergo intense star-formation activity, while ETGs are mostly inactive. We question this general rule and investigate the existence of star-forming ETGs and quiescent LTGs in the local Universe. By computing the physical properties of 2,209 such galaxies in the GAMA survey being morphologically classified and using information on their structural properties as well as the density of their local environment, we seek for understanding the differences from their 'typical' counterparts. We separate galaxies into subsets based on their dominant ionising process, making use of criteria based on the WH$_α$ width and the [NII/H$_α$] ratio. Taking advantage of the SED fitting code CIGALE we derive galaxy properties, such as the $M_\text{star}$, $M_\text{dust}$, and SFR and also estimate the unattenuated and the dust-absorbed stellar emission, for both the young and old stellar populations. Ongoing star-formation activity is found in 47% of ETGs and 8% of LTGs are quiescent. The star-forming E galaxies, together with the LBSs, constitute a population that follows very well the SFMS of spiral galaxies. The fraction of the luminosity originating from young stars in the star-forming ETGs is quite substantial ($\sim$ 25%) and similar to that of the star-forming LTGs. Investigating possible differences between star-forming and quiescent galaxies we find that the intrinsic shape of the SED of the star-forming galaxies is, on average, very similar for all morphological types. Concerning their structural parameters, quiescent galaxies tend to show larger values of the Sérsic index and larger $R_\text{eff}$ (compared to star-forming galaxies). Finally, we find that star-forming galaxies preferably reside in lower-density environments compared to the quiescent ones, which exhibit a higher percentage of sources being members of groups. △ Less

Submitted 27 September, 2022; originally announced September 2022.

Journal ref: A&A 669, A11 (2023)

Dust emissivity in resolved spiral galaxies

Authors: S. Bianchi, V. Casasola, E. Corbelli, F. Galliano, L. Magrini, A. Nersesian, F. Salvestrini, M. Baes, L. P. Cassara, C. J. R. Clark, I. De Looze, A. P. Jones, S. C. Madden, A. Mosenkov, N. Ysard

Abstract: Context: The far-infrared (FIR) and sub-millimeter (submm) emissivity of the Milky Way (MW) cirrus is an important benchmark for dust grain models. Dust masses in other galaxies are generally derived from the FIR/submm using the emission properties of these MW-calibrated models. Aims: We seek to derive the FIR/submm emissivity in nine nearby spiral galaxies to check its compatibility with MW cirru… ▽ More Context: The far-infrared (FIR) and sub-millimeter (submm) emissivity of the Milky Way (MW) cirrus is an important benchmark for dust grain models. Dust masses in other galaxies are generally derived from the FIR/submm using the emission properties of these MW-calibrated models. Aims: We seek to derive the FIR/submm emissivity in nine nearby spiral galaxies to check its compatibility with MW cirrus measurements. Methods: We obtained values of the emissivity at 70 to 500 um, using maps of dust emission from the Herschel satellite and of gas surface density from the THINGS and HERACLES surveys on a scale generally corresponding to 440 pc. We studied the variation of the emissivity with the surface brightness ratio I(250um)/I(500um), a proxy for the intensity of the interstellar radiation field heating the dust. Results: We find that the average value of the emissivity agrees with MW estimates for pixels sharing the same color as the cirrus, namely, for I(250um)/I(500um) = 4.5. For I(250um)/I(500um) > 5, the measured emissivity is instead up to a factor ~2 lower than predicted from MW dust models heated by stronger radiation fields. Regions with higher I(250um)/I(500um) are preferentially closer to the galactic center and have a higher overall (stellar+gas) surface density and molecular fraction. The results do not depend strongly on the adopted CO-to-molecular conversion factor and do not appear to be affected by the mixing of heating conditions. Conclusions: Our results confirm the validity of MW dust models at low density, but are at odds with predictions for grain evolution in higher density environments. If the lower-than-expected emissivity at high I(250um)/I(500um) is the result of intrinsic variations in the dust properties, it would imply an underestimation of the dust mass surface density of up to a factor ~2 when using current dust models. △ Less

Submitted 27 June, 2022; originally announced June 2022.

Comments: 16 pages, 11 figures, A&A accepted

Journal ref: A&A 664, A187 (2022)

Dust Emission in Galaxies at Millimeter Wavelengths: Cooling of star forming regions in NGC6946

Authors: G. Ejlali, R. Adam, P. Ade, H. Ajeddig, P. André, E. Artis, H. Ausse, A. Beelen, A. Benoît, S. Berta, L. Bing, O. Bourrion, M. Calvo, A. Catalano, I. de Looze, M. De Petris, F. -X. Désert, S. Doyle, E. F. C. Driessen, M. Galametz, F. Galliano, A. Gomez, J. Goupy, A. P. Jones, A. Hughes , et al. (32 additional authors not shown)

Abstract: Interstellar dust plays an important role in the formation of molecular gas and the heating and cooling of the interstellar medium. The spatial distribution of the mm-wavelength dust emission from galaxies is largely unexplored. The NIKA2 Guaranteed Time Project IMEGIN (Interpreting the Millimeter Emission of Galaxies with IRAM and NIKA2) has recently mapped the mm emission in the grand design spi… ▽ More Interstellar dust plays an important role in the formation of molecular gas and the heating and cooling of the interstellar medium. The spatial distribution of the mm-wavelength dust emission from galaxies is largely unexplored. The NIKA2 Guaranteed Time Project IMEGIN (Interpreting the Millimeter Emission of Galaxies with IRAM and NIKA2) has recently mapped the mm emission in the grand design spiral galaxy NGC6946. By subtracting the contributions from the free-free, synchrotron, and CO line emission, we map the distribution of the pure dust emission at 1:15mm and 2mm. Separating the arm/interarm regions, we find a dominant 2mm emission from interarms indicating the significant role of the general interstellar radiation field in heating the cold dust. Finally, we present maps of the dust mass, temperature, and emissivity index using the Bayesian MCMC modeling of the spectral energy distribution in NGC6946. △ Less

Submitted 6 November, 2021; originally announced November 2021.

Comments: To appear in the Proceedings of the International Conference entitled "mm Universe @ NIKA2", Rome (Italy), June 2021, EPJ Web of conferences

Exploring the millimetre emission in nearby galaxies: analysis of the edge-on galaxy NGC 891

Authors: S. Katsioli, R. Adam, P. Ade, H. Ajeddig, P. André, E. Artis, H. Aussel, A. Beelen, A. Benoît, S. Berta, L. Bing, O. Bourrion, M. Calvo, A. Catalano, I. De Looze, M. De Petris, F. -X. Désert, S. Doyle, E. F. C. Driessen, G. Ejlali, M. Galametz, F. Galliano, A. Gomez, J. Goupy, A. P. Jones , et al. (32 additional authors not shown)

Abstract: New observations of the edge-on galaxy NGC 891, at 1.15 and 2 mm obtained with the IRAM 30-m telescope and the NIKA2 camera, within the framework of the IMEGIN (Interpreting the Millimetre Emission of Galaxies with IRAM and NIKA2) Large Program, are presented in this work. By using multiwavelength maps (from the mid-IR to the cm wavelengths) we perform SED fitting in order to extract the physical… ▽ More New observations of the edge-on galaxy NGC 891, at 1.15 and 2 mm obtained with the IRAM 30-m telescope and the NIKA2 camera, within the framework of the IMEGIN (Interpreting the Millimetre Emission of Galaxies with IRAM and NIKA2) Large Program, are presented in this work. By using multiwavelength maps (from the mid-IR to the cm wavelengths) we perform SED fitting in order to extract the physical properties of the galaxy on both global and local ($\sim$kpc) scales. For the interpretation of the observations we make use of a state-of-the-art SED fitting code, HerBIE (HiERarchical Bayesian Inference for dust Emission). The observations indicate a galaxy morphology, at mm wavelengths, similar to that of the cold dust emission traced by sub-mm observations and to that of the molecular gas. The contribution of the radio emission at the NIKA2 bands is very small (negligible at 1.15 mm and $\sim10\%$ at 2 mm) while it dominates the total energy budget at longer wavelengths (beyond 5 mm). On local scales, the distribution of the free-free emission resembles that of the dust thermal emission while the distribution of the synchrotron emission shows a deficiency along the major axis of the disc of the galaxy. △ Less

Submitted 2 November, 2021; originally announced November 2021.

Comments: To appear in the Proceedings of the International Conference entitled "mm Universe @ NIKA2", Rome (Italy), June 2021, EPJ Web of conferences

Modelling the cold dust in nearby spiral galaxies with radiative transfer

Authors: Angelos Nersesian, Maarten Baes, Suzanne C. Madden

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

Submitted 2 November, 2021; originally announced November 2021.

Comments: To appear in the Proceedings of the International Conference entitled "mm Universe @ NIKA2", Rome (Italy), June 2021, EPJ Web of conferences

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.

The physical properties of local (U)LIRGs: a comparison with nearby early- and late-type galaxies

Authors: E. -D. Paspaliaris, E. M. Xilouris, A. Nersesian, V. A. Masoura, M. Plionis, I. Georgantopoulos, S. Bianchi, S. Katsioli, G. Mountrichas

Abstract: In order to pinpoint the place of the (U)LIRGs in the local Universe we examine the properties of a sample of 67 such systems and compare them with those of 268 ETGs and 542 LTGs from the DustPedia database. We make use of multi-wavelength photometric data and the CIGALE SED fitting code to extract their physical parameters. The median SEDs as well as the values of the derived parameters were comp… ▽ More In order to pinpoint the place of the (U)LIRGs in the local Universe we examine the properties of a sample of 67 such systems and compare them with those of 268 ETGs and 542 LTGs from the DustPedia database. We make use of multi-wavelength photometric data and the CIGALE SED fitting code to extract their physical parameters. The median SEDs as well as the values of the derived parameters were compared to those of the local ETGs and LTGs. In addition to that, (U)LIRGs were divided into seven classes, according to the merging stage of each system, and variations in the derived parameters were investigated. (U)LIRGs occupy the `high-end' on the dust and stellar mass, and SFR in the local Universe with median values of 5.2$\times10^7~M_{\odot}$, 6.3$\times10^{10}~M_{\odot}$ and 52$~M_{\odot}$yr$^{-1}$, respectively. The PDR-dust emission in (U)LIRGs is 11.7% of the total dust luminosity, significantly higher than ETGs (1.6%) and the LTGs (5.2%). The median value of the dust temperature in (U)LIRGs is 32 K, which is higher compared to both the ETGs (28 K) and the LTGs (22 K). Small differences, in the derived parameters, are seen for the seven merging classes of our sample of (U)LIRGs with the most evident one being on the star-formation rate, where in systems in late merging stages the median SFR reaches up to 99 M$_{\odot}$ yr$^{-1}$ compared to 26 M$_{\odot}$ yr$^{-1}$ for the isolated ones. In contrast to the local normal galaxies where old stars dominate the stellar emission, the young stars in (U)LIRGs contribute with 64% of their luminosity to the total stellar luminosity. The fraction of the dust-absorbed stellar luminosity is extremely high in (U)LIRGs (78%) compared to 7% and 25% in ETGs and ETGs, respectively. The fraction of the stellar luminosity used to heat up the dust grains is very high in (U)LIRGs, while 74% of the dust emission comes from the young stars. △ Less

Submitted 7 February, 2021; originally announced February 2021.

Journal ref: A&A 649, A137 (2021)

A Nearby Galaxy Perspective on Dust Evolution. Scaling relations and constraints on the dust build-up in galaxies with the DustPedia and DGS samples

Authors: Frédéric Galliano, Angelos Nersesian, Simone Bianchi, Ilse De Looze, Sambit Roychowdhury, Maarten Baes, Viviana Casasola, Letizia, P. Cassará, Wouter Dobbels, Jacopo Fritz, Maud Galametz, Anthony P. Jones, Suzanne C. Madden, Aleksandr Mosenkov, Emmanuel M. Xilouris, Nathalie Ysard

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

Submitted 2 January, 2021; originally announced January 2021.

Comments: 64 pages, 36 figures, accepted for publication in Astronomy & Astrophysics

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)

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 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)

Dust emissivity and absorption cross section in DustPedia late-type galaxies

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

Abstract: Aims: We compare the far-infrared to sub-millimetre dust emission properties measured in high Galactic latitude cirrus with those determined in a sample of 204 late-type DustPedia galaxies. The aim is to verify if it is appropriate to use Milky Way dust properties to derive dust masses in external galaxies. Methods: We used Herschel observations and atomic and molecular gas masses to estimate the… ▽ More Aims: We compare the far-infrared to sub-millimetre dust emission properties measured in high Galactic latitude cirrus with those determined in a sample of 204 late-type DustPedia galaxies. The aim is to verify if it is appropriate to use Milky Way dust properties to derive dust masses in external galaxies. Methods: We used Herschel observations and atomic and molecular gas masses to estimate the disc-averaged dust emissivity at 250 micrometres, and from this, the absorption cross section per H atom and per dust mass. The emissivity requires one assumption, which is the CO-to-H_2 conversion factor, and the dust temperature is additionally required for the absorption cross section per H atom; yet another constraint on the dust-to-hydrogen ratio D/H, depending on metallicity, is required for the absorption cross section dust mass. Results: We find epsilon(250) = 0.82 +/- 0.07 MJy sr^-1 (1E20 H cm^-2)^-1 for galaxies with 4 < F(250)/F(500) < 5. This depends only weakly on the adopted CO-to-H_2 conversion factor. The value is almost the same as that for the Milky Way at the same colour ratio. Instead, for F(250)/F(500) > 6, epsilon(250) is lower than predicted by its dependence on the heating conditions. The reduction suggests a variation in dust emission properties for spirals of earlier type, higher metallicity, and with a higher fraction of molecular gas. When the standard emission properties of Galactic cirrus are used for these galaxies, their dust masses might be underestimated by up to a factor of two. Values for the absorption cross sections at the Milky Way metallicity are also close to those of the cirrus. Mild trends of the absorption cross sections with metallicity are found, although the results depend on the assumptions made. △ Less

Submitted 27 September, 2019; originally announced September 2019.

Comments: 13 pages, 6 figures, 2 tables, A&A accepted

Journal ref: A&A 631, A102 (2019)

DustPedia - the relationships between stars, gas and dust for galaxies residing in different environments

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

Abstract: We use a sub-set of the DustPedia galaxy sample (461 galaxies) to investigate the effect the environment has had on galaxies. We consider Virgo cluster and field samples and also assign a density contrast parameter to each galaxy, as defined by the local density of SDSS galaxies. We consider their chemical evolution (using M_{Dust}/M_{Baryon} and M_{Gas}/M_{Baryon}), their specific star formation… ▽ More We use a sub-set of the DustPedia galaxy sample (461 galaxies) to investigate the effect the environment has had on galaxies. We consider Virgo cluster and field samples and also assign a density contrast parameter to each galaxy, as defined by the local density of SDSS galaxies. We consider their chemical evolution (using M_{Dust}/M_{Baryon} and M_{Gas}/M_{Baryon}), their specific star formation rate (SFR/M_{Stars}), star formation efficiency (SFR/M_{Gas}), stars-to-dust mass ratio (M_{Stars}/M_{Dust}), gas-to-dust mass ratio (M_{Gas}/M_{Dust}) and the relationship between star formation rate per unit mass of dust and dust temperature (SFR/M_{Dust} and T_{Dust}). Late type galaxies (later than Sc) in all of the environments can be modelled using simple closed box chemical evolution and a simple star formation history (SFR(t) \propto t\exp{-t/τ}). For earlier type galaxies the physical mechanisms that give rise to their properties are clearly much more varied and require a more complicated model (mergers, gas in or outflow). However, we find little or no difference in the properties of galaxies of the same morphological type within the cluster, field or with different density contrasts. It appears that it is morphology, how and whenever this is laid down, and consistent internal physical processes that primarily determine the derived properties of galaxies in the DustPedia sample and not processes related to differences in the local environment. △ Less

Submitted 30 April, 2019; originally announced April 2019.

Comments: Accepted for publication in Astronomy and Astrophysics

Revealing the dust attenuation properties on resolved scales in NGC628 with SWIFT UVOT data

Authors: Marjorie Decleir, Ilse De Looze, Médéric Boquien, Maarten Baes, Sam Verstocken, Daniela Calzetti, Laure Ciesla, Jacopo Fritz, Rob Kennicutt, Angelos Nersesian, Mathew Page

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

Submitted 15 March, 2019; originally announced March 2019.

Comments: 26 pages, 16 figures, 4 tables, accepted for publication in MNRAS

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)

A systematic metallicity study of DustPedia galaxies reveals evolution in the dust-to-metal ratios

Authors: P. De Vis, A. Jones, S. Viaene, V. Casasola, C. J. R. Clark, M. Baes, S. Bianchi, L. P. Cassara, J. I. Davies, I. De Looze, M. Galametz, F. Galliano, S. Lianou, S. Madden, A. Manilla-Robles, A. V. Mosenkov, A. Nersesian, S. Roychowdhury, E. M. Xilouris, N. Ysard

Abstract: Observations of evolution in the dust-to-metal ratio allow us to constrain the dominant dust processing mechanisms. In this work, we present a study of the dust-to-metal and dust-to-gas ratios in a sub-sample of ~500 DustPedia galaxies. Using literature and MUSE emission line fluxes, we derived gas-phase metallicities (oxygen abundances) for over 10000 individual regions and determine characterist… ▽ More Observations of evolution in the dust-to-metal ratio allow us to constrain the dominant dust processing mechanisms. In this work, we present a study of the dust-to-metal and dust-to-gas ratios in a sub-sample of ~500 DustPedia galaxies. Using literature and MUSE emission line fluxes, we derived gas-phase metallicities (oxygen abundances) for over 10000 individual regions and determine characteristic metallicities for each galaxy. We study how the relative dust, gas, and metal contents of galaxies evolve by using metallicity and gas fraction as proxies for evolutionary state. The global oxygen abundance and nitrogen-to-oxygen ratio are found to increase monotonically as galaxies evolve. Additionally, unevolved galaxies (gas fraction > 60%, metallicity 12 + log(O/H) < 8.2) have dust-to-metal ratios that are about a factor of 2.1 lower (a factor of six lower for galaxies with gas fraction > 80%) than the typical dust-to-metal ratio (Md/MZ ~ 0.214) for more evolved sources. However, for high gas fractions, the scatter is larger due to larger observational uncertainties as well as a potential dependence of the dust grain growth timescale and supernova dust yield on local conditions and star formation histories. We find chemical evolution models with a strong contribution from dust grain growth describe these observations reasonably well. The dust-to-metal ratio is also found to be lower for low stellar masses and high specific star formation rates (with the exception of some sources undergoing a starburst). Finally, the metallicity gradient correlates weakly with the HI-to-stellar mass ratio, the effective radius and the dust-to-stellar mass ratio, but not with stellar mass. △ Less

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

Comments: 24 pages, 15 figures, published in A&A, bibref: 2019A&A...623A...5D, This Version 2 has increased accuracy in Table B1 in appendix B

Journal ref: A&A 623, A5 (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)

Long-term optical variability of high-mass X-ray binaries. II. Spectroscopy

Authors: P. Reig, A. Nersesian, A. Zezas, L. Gkouvelis, M. J. Coe

Abstract: We present the results of our monitoring program to study the long-term variability of the Halpha line in high-mass X-ray binaries. We have carried out the most complete optical spectroscopic study of the global properties of high-mass X-ray binaries so far with the analysis of more than 1100 spectra of 20 sources. Our aim is to characterise the optical variability timescales and study the interac… ▽ More We present the results of our monitoring program to study the long-term variability of the Halpha line in high-mass X-ray binaries. We have carried out the most complete optical spectroscopic study of the global properties of high-mass X-ray binaries so far with the analysis of more than 1100 spectra of 20 sources. Our aim is to characterise the optical variability timescales and study the interaction between the neutron star and the accreting material. Our results can be summarised as follows: i) we find that Be/X-ray binaries with narrow orbits are more variable than systems with long orbital periods, ii) we show that a Keplerian distribution of the gas particles provides a good description of the disks in Be/X-ray binaries, as it does in classical Be stars, iii) a decrease in the Halpha equivalent width is generally observed after major X-ray outbursts, iv) we confirm that the Halpha equivalent width correlates with disk radius, v) while systems with supergiant companions display, multi-structured profiles, most of the Be/X-ray binaries show at some epoch double-peak asymmetric profiles, indicating that density inhomogeneities is a common property in the disk of Be/X-ray binaries, vi) the profile variability (V/R ratio) timescales are shorter and the Halpha equivalent width are smaller in Be/X-ray binaries than in isolated Be stars, and vii) we provide new evidence that the disk in Be/X-ray binaries is on average denser than in classical Be stars. △ Less

Submitted 25 April, 2016; v1 submitted 28 March, 2016; originally announced March 2016.

Comments: new version after language editting and proofs. Reference updated. Accepted for publication in A&A

Journal ref: A&A 590, A122 (2016)