-
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.
-
Reduction of dust radial drift by turbulence in protoplanetary disks
Authors:
Fabiola Antonietta Gerosa,
Jérémie Bec,
Héloïse Méheut,
Anand Utsav Kapoor
Abstract:
Dust particles in protoplanetary disks, lacking support from pressure, rotate at velocities exceeding those of the surrounding gas. Consequently, they experience a head-wind from the gas that drives them toward the central star. Radial drift occurs on timescales much shorter than those inferred from disk observations or those required for dust to aggregate and form planets. Additionally, turbulenc…
▽ More
Dust particles in protoplanetary disks, lacking support from pressure, rotate at velocities exceeding those of the surrounding gas. Consequently, they experience a head-wind from the gas that drives them toward the central star. Radial drift occurs on timescales much shorter than those inferred from disk observations or those required for dust to aggregate and form planets. Additionally, turbulence is often assumed to amplify the radial drift of dust in planet-forming disks when modeled through an effective viscous transport. However, the local interactions between turbulent eddies and particles are known to be significantly more intricate than in a viscous fluid. Our objective is to elucidate and characterize the dynamic effects of Keplerian turbulence on the mean radial and azimuthal velocities of dust particles. We employ 2D shearing-box incompressible simulations of the gas, which is maintained in a developed turbulent state while rotating at a sub-Keplerian speed. Dust is modeled as Lagrangian particles set at a Keplerian velocity, therefore experiencing a radial force toward the star through drag. Turbulent eddies are found to reduce the radial drift, while simultaneously enhancing the azimuthal velocities of small particles. This dynamic behavior arises from the modification of dust trajectories due to turbulent eddies.
△ Less
Submitted 17 April, 2024;
originally announced April 2024.
-
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.
-
Spatially resolved mock observations of stellar kinematics: full radiative transfer treatment of simulated galaxies
Authors:
Daniela Barrientos Acevedo,
Arjen van der Wel,
Maarten Baes,
Robert J. J Grand,
Anand Utsav Kapoor,
Peter Camps,
Anna de Graaff,
Caroline M. S. Straatman,
Rachel Bezanson
Abstract:
We present a framework to build realistic mock spectroscopic observations for state-of-the-art hydrodynamical simulations, using high spectral resolution stellar population models and full radiative transfer treatment with SKIRT. As a first application we generate stellar continuum mock observations for the Auriga cosmological zoom simulations emulating integral-field observations from the SAMI Ga…
▽ More
We present a framework to build realistic mock spectroscopic observations for state-of-the-art hydrodynamical simulations, using high spectral resolution stellar population models and full radiative transfer treatment with SKIRT. As a first application we generate stellar continuum mock observations for the Auriga cosmological zoom simulations emulating integral-field observations from the SAMI Galaxy survey. We perform spectral fitting on our synthetic cubes and compute the resulting rotation velocity ($V_{\rm{rot}}$) and velocity dispersion within 1$R_{\text{e}}$ ($σ_{\text{e}}$) for a sub-set of the Auriga sample. We find that the kinematics produced by Auriga are in good agreement with the observations from the SAMI Galaxy survey after taking into account the effects of dust and the systematics produced by the observation limitations. We also explore the effects of seeing convolution, inclination, and attenuation on the line-of-sight velocity distribution. For highly inclined galaxies, these effects can lead to an artificial decrease in the measured $V/σ$ by nearly a factor two (after inclination correction). We also demonstrate the utility of our method for high-redshift galaxies by emulating spatially resolved continuum spectra from the LEGA-C survey and, looking forward, E-ELT HARMONI. Our framework represents a crucial link between the ground truth for stellar populations and kinematics in simulations and the observed stellar continuum observations at low and high redshift.
△ Less
Submitted 13 June, 2023;
originally announced June 2023.
-
The atomic-to-molecular hydrogen transition in the TNG50 simulation: Using realistic UV fields to create spatially resolved HI maps
Authors:
Andrea Gebek,
Maarten Baes,
Benedikt Diemer,
W. J. G. de Blok,
Dylan Nelson,
Anand Utsav Kapoor,
Peter Camps,
Omphile Rabyang,
Lerothodi Leeuw
Abstract:
Cold gas in galaxies provides a crucial test to evaluate the realism of cosmological hydrodynamical simulations. To extract the atomic and molecular hydrogen properties of the simulated galaxy population, postprocessing methods taking the local UV field into account are required. We improve upon previous studies by calculating realistic UV fields with the dust radiative transfer code SKIRT to mode…
▽ More
Cold gas in galaxies provides a crucial test to evaluate the realism of cosmological hydrodynamical simulations. To extract the atomic and molecular hydrogen properties of the simulated galaxy population, postprocessing methods taking the local UV field into account are required. We improve upon previous studies by calculating realistic UV fields with the dust radiative transfer code SKIRT to model the atomic-to-molecular transition in TNG50, the highest-resolution run of the IllustrisTNG suite. Comparing integrated quantities such as the HI mass function, we study to what detail the UV field needs to be modelled in order to calculate realistic cold gas properties. We then evaluate new, spatially resolved comparisons for cold gas in galaxies by exploring synthetic maps of atomic hydrogen at redshift zero and compare them to 21-cm observations of local galaxies from the WHISP survey. In terms of non-parametric morphologies, we find that TNG50 HI maps are less concentrated than their WHISP counterparts (median $ΔC\approx0.3$), due in part to central HI deficits related to the ejective character of supermassive black hole feedback in TNG. In terms of the HI column density distribution function, we find discrepancies between WHISP and IllustrisTNG that depend on the total HI abundance in these datasets as well as the postprocessing method. To fully exploit the synergy between cosmological simulations and upcoming deep HI/H2 data, we advocate the use of accurate methods to estimate the UV radiation field and to generate mock maps.
△ Less
Submitted 13 March, 2023;
originally announced March 2023.
-
UV to submillimetre luminosity functions of TNG50 galaxies
Authors:
Ana Trcka,
Maarten Baes,
Peter Camps,
Anand Utsav Kapoor,
Dylan Nelson,
Annalisa Pillepich,
Daniela Barrientos,
Lars Hernquist,
Federico Marinacci,
Mark Vogelsberger
Abstract:
We apply the radiative transfer (RT) code SKIRT on a sample of ~14000 low-redshift (z<= 0.1) galaxies extracted from the TNG50 simulation to enable an apples-to-apples comparison with observations. The RT procedure is calibrated via comparison of a subsample of TNG50 galaxies with the DustPedia observational sample: we compare several luminosity and colour scaling relations and spectral energy dis…
▽ More
We apply the radiative transfer (RT) code SKIRT on a sample of ~14000 low-redshift (z<= 0.1) galaxies extracted from the TNG50 simulation to enable an apples-to-apples comparison with observations. The RT procedure is calibrated via comparison of a subsample of TNG50 galaxies with the DustPedia observational sample: we compare several luminosity and colour scaling relations and spectral energy distributions in different specific SFR bins. We consistently derive galaxy luminosity functions for the TNG50 simulation in 14 broadband filters from UV to submillimetre wavelengths and investigate the effects of the aperture, orientation, radiative transfer recipe, and numerical resolution. We find that, while our TNG50+RT fiducial model agrees well with the observed luminosity functions at the knee (+/- 0.04 dex typical agreement), the TNG50+RT luminosity functions evaluated within 5R_1/2 are generally higher than observed at both the faint and bright ends, by 0.004 (total IR)-0.27 (UKIDSS H) dex and 0.12 (SPIRE250)-0.8 (GALEX FUV) dex, respectively. A change in the aperture does affect the bright end of the luminosity function, easily by up to 1 dex depending on the choice. However, we also find that the galaxy luminosity functions of a worse-resolution run of TNG50 (TNG50-2, with 8 times worse mass resolution than TNG50, similar to TNG100) are in better quantitative agreement with observational constraints. Finally, we publicly release the photometry for the TNG50 sample in 53 broadbands from FUV to submillimetre, in three orientations and four apertures, as well as galaxy spectral energy distributions.
△ Less
Submitted 12 August, 2022;
originally announced August 2022.
-
High-resolution synthetic UV-submm images for Milky Way-mass simulated galaxies from the ARTEMIS project
Authors:
Peter Camps,
Anand Utsav Kapoor,
Ana Trčka,
Andreea S. Font,
Ian G. McCarthy,
James Trayford,
Maarten Baes
Abstract:
We present redshift-zero synthetic dust-aware observations for the 45 Milky Way-mass simulated galaxies of the ARTEMIS project, calculated with the SKIRT radiative transfer code. The post-processing procedure includes components for star-forming regions, stellar sources, and diffuse dust. We produce and publicly release realistic high-resolution images for 50 commonly-used broadband filters from u…
▽ More
We present redshift-zero synthetic dust-aware observations for the 45 Milky Way-mass simulated galaxies of the ARTEMIS project, calculated with the SKIRT radiative transfer code. The post-processing procedure includes components for star-forming regions, stellar sources, and diffuse dust. We produce and publicly release realistic high-resolution images for 50 commonly-used broadband filters from ultraviolet to sub-millimetre wavelengths and for 18 different viewing angles. We compare the simulated ARTEMIS galaxies to observed galaxies in the DustPedia database with similar stellar mass and star formation rate, and to synthetic observations of the simulated galaxies of the Auriga project produced in previous work using a similar post-processing technique. In all cases, global galaxy properties are derived using SED fitting. We find that, similar to Auriga, the post-processed ARTEMIS galaxies generally reproduce the observed scaling relations for global fluxes and physical properties, although dust extinction at FUV/UV wavelengths is underestimated and representative dust temperatures are lower than observed. At a resolved scale, we compare multi-wavelength non-parametric morphological properties of selected disc galaxies across the data sets. We find that the ARTEMIS galaxies largely reproduce the observed morphological trends as a function of wavelength, although they appear to be more clumpy and less symmetrical than observed. We note that the ARTEMIS and Auriga galaxies occupy adjacent regions in the specific star formation versus stellar mass plane, so that the synthetic observation data sets supplement each other.
△ Less
Submitted 11 March, 2022;
originally announced March 2022.
-
A new analytical scattering phase function for interstellar dust
Authors:
Maarten Baes,
Peter Camps,
Anand Utsav Kapoor
Abstract:
Context: Properly modelling scattering by interstellar dust grains requires a good characterisation of the scattering phase function. The Henyey-Greenstein phase function has become the standard for describing anisotropic scattering by dust grains, but it is a poor representation of the real scattering phase function outside the optical range. Aims: We investigate alternatives for the Henyey-Green…
▽ More
Context: Properly modelling scattering by interstellar dust grains requires a good characterisation of the scattering phase function. The Henyey-Greenstein phase function has become the standard for describing anisotropic scattering by dust grains, but it is a poor representation of the real scattering phase function outside the optical range. Aims: We investigate alternatives for the Henyey-Greenstein phase function that would allow the scattering properties of dust grains to be described. Our goal is to find a balance between realism and complexity: the scattering phase function should be flexible enough to provide an accurate fit to the scattering properties of dust grains over a wide wavelength range, and it should be simple enough to be easy to handle, especially in the context of radiative transfer calculations. Methods: We fit various analytical phase functions to the scattering phase function corresponding to the BARE-GR-S model, one of the most popular and commonly adopted models for interstellar dust. We weigh the accuracy of the fit against the number of free parameters in the analytical phase functions. Results: We confirm that the Henyey-Greenstein phase functions poorly describe scattering by dust grains, particularly at ultraviolet (UV) wavelengths, with relative differences of up to 50%. The Draine phase function alleviates this problem at near-infrared (NIR) wavelengths, but not in the UV. The two-term Reynolds-McCormick phase function, recently advocated in the context of light scattering in nanoscale materials and aquatic media, describes the BARE-GR-S data very well, but its five free parameters are degenerate. We propose a simpler phase function, the two-term ultraspherical-2 (TTU2) phase function, that also provides an excellent fit to the BARE-GR-S phase function over the entire UV-NIR wavelength range. (Abridged)
△ Less
Submitted 3 February, 2022;
originally announced February 2022.
-
High-resolution synthetic UV-submm images for simulated Milky Way-type galaxies from the Auriga project
Authors:
Anand Utsav Kapoor,
Peter Camps,
Maarten Baes,
Ana Trcka,
Robert J. J. Grand,
Arjen van der Wel,
Luca Cortese,
Ilse De Looze,
Daniela Barrientos
Abstract:
We present redshift-zero synthetic observational data considering dust attenuation and dust emission for the thirty galaxies of the Auriga project, calculated with the SKIRT radiative transfer code. The post-processing procedure includes components for star-forming regions, stellar sources, and diffuse dust taking into account stochastic heating of dust grains. This allows us to obtain realistic h…
▽ More
We present redshift-zero synthetic observational data considering dust attenuation and dust emission for the thirty galaxies of the Auriga project, calculated with the SKIRT radiative transfer code. The post-processing procedure includes components for star-forming regions, stellar sources, and diffuse dust taking into account stochastic heating of dust grains. This allows us to obtain realistic high-resolution broadband images and fluxes from ultraviolet to sub-millimeter wavelengths. For the diffuse dust component, we consider two mechanisms for assigning dust to gas cells in the simulation. In one case, only the densest or the coldest gas cells are allowed to have dust, while in the other case this condition is relaxed to allow a larger number of dust-containing cells. The latter approach yields galaxies with a larger radial dust extent and an enhanced dust presence in the inter-spiral regions. At a global scale, we compare Auriga galaxies with observations by deriving dust scaling relations using SED fitting. At a resolved scale, we make a multi-wavelength morphological comparison with nine well-resolved spiral galaxies from the DustPedia observational database. We find that for both dust assignment methods, although the Auriga galaxies show a good overall agreement with observational dust properties, they exhibit a slightly higher specific dust mass. The multi-wavelength morphological analysis reveals a good agreement between the Auriga and the observed galaxies in the optical wavelengths. In the mid and far-infrared wavelengths, Auriga galaxies appear smaller and more centrally concentrated in comparison to their observed counterparts. We publicly release the multi-observer images and fluxes in 50 commonly used broadband filters.
△ Less
Submitted 14 July, 2021;
originally announced July 2021.
-
Polarised emission from aligned dust grains in nearby galaxies: predictions from the Auriga simulations
Authors:
Bert Vandenbroucke,
Maarten Baes,
Peter Camps,
Anand Utsav Kapoor,
Daniela Barrientos,
Jean-Philippe Bernard
Abstract:
(Aims) In this work, we predict the far-infrared polarisation signal emitted by non-spherical dust grains in nearby galaxies. We determine the angular resolution and sensitivity required to study the magnetic field configuration in these galaxies. (Methods) We post-process a set of Milky Way like galaxies from the Auriga project, assuming a dust mix consisting of spheroidal dust grains that are pa…
▽ More
(Aims) In this work, we predict the far-infrared polarisation signal emitted by non-spherical dust grains in nearby galaxies. We determine the angular resolution and sensitivity required to study the magnetic field configuration in these galaxies. (Methods) We post-process a set of Milky Way like galaxies from the Auriga project, assuming a dust mix consisting of spheroidal dust grains that are partially aligned with the model magnetic field. We constrain our dust model using Planck 353 GHz observations of the Milky Way. This model is then extrapolated to shorter wavelengths that cover the peak of interstellar dust emission and to observations of arbitrarily oriented nearby Milky Way like galaxies. (Results) Assuming an intrinsic linear polarisation fraction that does not vary significantly with wavelength for wavelengths longer than 50 micron, we predict a linear polarisation fraction with a maximum of $10-15\%$ and a median value of $\approx{}7\%$ for face-on galaxies and $\approx{}3\%$ for edge-on galaxies. The polarisation fraction anti-correlates with the line of sight density and with the angular dispersion function which expresses the large scale order of the magnetic field perpendicular to the line of sight. The maximum linear polarisation fraction agrees well with the intrinsic properties of the dust model. The true magnetic field orientation can be traced along low density lines of sight when it is coherent along the line of sight. These results also hold for nearby galaxies, where a coherent magnetic field structure is recovered over a range of different broad bands. (Conclusions) Polarised emission from non-spherical dust grains accurately traces the large scale structure of the galactic magnetic field in Milky Way like galaxies, with expected maximum linear polarisation fractions of $10-15\%$. To resolve this maximum, a spatial resolution of at least 1 kpc is required.
△ Less
Submitted 5 July, 2021;
originally announced July 2021.
-
Effects of Spatial Discretization in Lyman-alpha Line Radiation Transfer Simulations
Authors:
Peter Camps,
Christoph Behrens,
Maarten Baes,
Anand Utsav Kapoor,
Robert Grand
Abstract:
We describe the addition of Lyman-alpha resonant line transfer to our dust continuum radiation transfer code SKIRT, verifying our implementation with published results for spherical problems and using some self-designed three-dimensional setups. We specifically test spatial discretization through various grid types, including hierarchical octree grids and unstructured Voronoi tessellations. We the…
▽ More
We describe the addition of Lyman-alpha resonant line transfer to our dust continuum radiation transfer code SKIRT, verifying our implementation with published results for spherical problems and using some self-designed three-dimensional setups. We specifically test spatial discretization through various grid types, including hierarchical octree grids and unstructured Voronoi tessellations. We then use a radiation transfer post-processing model for one of the spiral galaxies produced by the Auriga cosmological zoom simulations to investigate the effect of spatial discretization on the synthetic observations. We find that the calculated Lyman-alpha line profiles exhibit an extraordinarily strong dependence on the type and resolution of the spatial grid, rendering the results untrustworthy at best. We attribute this effect to the large gradients in the hydrogen density distribution over small distances, which remain significantly under-resolved in the input model. We therefore argue that further research is needed to determine the required spatial resolution of a hydrodynamical simulation snapshot to enable meaningful Lyman-alpha line transfer post-processing.
△ Less
Submitted 1 June, 2021;
originally announced June 2021.