Showing 1–50 of 61 results for author: Camps, P

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

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

Self-consistent dust and non-LTE line radiative transfer with SKIRT

Authors: Kosei Matsumoto, Peter Camps, Maarten Baes, Frederik De Ceuster, Keiichi Wada, Takao Nakagawa, Kentaro Nagamine

Abstract: We introduce Monte Carlo-based non-LTE line radiative transfer calculations in the 3D dust radiative transfer code SKIRT, which was originally set up as a dust radiative transfer code. By doing so, we develop a generic and powerful 3D radiative transfer code that can self-consistently generate spectra with molecular and atomic lines against the underlying continuum. We test the accuracy of the non… ▽ More We introduce Monte Carlo-based non-LTE line radiative transfer calculations in the 3D dust radiative transfer code SKIRT, which was originally set up as a dust radiative transfer code. By doing so, we develop a generic and powerful 3D radiative transfer code that can self-consistently generate spectra with molecular and atomic lines against the underlying continuum. We test the accuracy of the non-LTE line radiative transfer module in the extended SKIRT code using standard benchmarks. We find excellent agreement between the SKIRT results, the published benchmark results, and results obtained using the ray-tracing non-LTE line radiative transfer code MAGRITTE, which validates our implementation. We apply the extended SKIRT code on a 3D hydrodynamic simulation of a dusty AGN torus model and generate multi-wavelength images with CO rotational-line spectra against the underlying dust continuum. We find that the low-J CO emission traces the geometrically thick molecular torus, whereas the higher-J CO lines originate from the gas with high kinetic temperature located in the innermost regions of the torus. Comparing the calculations with and without dust radiative transfer, we find that higher-J CO lines are slightly attenuated by the surrounding cold dust when seen edge-on. This shows that atomic and molecular lines can experience attenuation, an effect that is particularly important for transitions at mid- and near-infrared wavelengths. Therefore, our self-consistent dust and non-LTE line radiative transfer calculations can help interpret the observational data from Herschel, ALMA, and JWST. △ Less

Submitted 5 September, 2023; originally announced September 2023.

Comments: 12 pages, 9 figures, Accepted for publication in A&A

Journal ref: A&A 678, A175 (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.

Comments: 17 pages, 11 figures. Accepted for publication in MNRAS

X-ray radiative transfer in full 3D with SKIRT

Authors: Bert Vander Meulen, Peter Camps, Marko Stalevski, Maarten Baes

Abstract: Models of active galactic nuclei (AGN) suggest that their circumnuclear media are complex with clumps and filaments, while recent observations hint towards polar extended structures of gas and dust, as opposed to the classical torus paradigm. The X-ray band forms an interesting window to study these circumnuclear media in great detail. In this work, we extended the radiative transfer code SKIRT wi… ▽ More Models of active galactic nuclei (AGN) suggest that their circumnuclear media are complex with clumps and filaments, while recent observations hint towards polar extended structures of gas and dust, as opposed to the classical torus paradigm. The X-ray band forms an interesting window to study these circumnuclear media in great detail. In this work, we extended the radiative transfer code SKIRT with the X-ray processes that govern the broadband X-ray spectra of obscured AGN, to study the structure of AGN circumnuclear media in full 3D, based on their reflected X-ray emission. We extended the SKIRT code with Compton scattering on free electrons, photo-absorption and fluorescence by cold atomic gas, scattering on bound electrons, and extinction by dust. This includes a novel treatment of extreme-forward scattering by dust, and a detailed description of anomalous Rayleigh scattering. To verify our X-ray implementation, we performed the first dedicated benchmark of X-ray torus models, comparing five X-ray radiative transfer codes. Finally, we illustrated the 3D nature of the code by producing synthetic X-ray images and spectra of clumpy torus models. SKIRT forms a powerful new tool to model AGN circumnuclear media in full 3D from X-ray to millimetre wavelengths, and is now publicly available. In the X-ray regime, we find an excellent agreement with the simulation results of the MYTorus and RefleX codes, which validates our X-ray implementation. We find some discrepancies with other codes, which motivates the need for a robust framework that can handle non-linear 3D radiative transfer effects. The new X-ray functionalities of the SKIRT code allow for uncomplicated access to a broad suite of 3D X-ray models for AGN that can easily be tested and modified. This will be particularly useful with the advent of X-ray microcalorimeter observations. △ Less

Submitted 20 April, 2023; originally announced April 2023.

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

Journal ref: A&A 674, A123 (2023)

Simulated observations of star formation regions: infrared evolution of globally collapsing clouds

Authors: Jesús M. Jáquez-Domínguez, Roberto Galván-Madrid, Jacopo Fritz, Manuel Zamora-Avilés, Peter Camps, Gustavo Bruzual, Maarten Baes, Yuxin Lin, Enrique Vázquez-Semadeni

Abstract: The direct comparison between hydrodynamical simulations and observations is needed to improve the physics included in the former and test biases in the latter. Post-processing radiative transfer and synthetic observations are now the standard way to do this. We report on the first application of the \texttt{SKIRT} radiative transfer code to simulations of a star-forming cloud. The synthetic obser… ▽ More The direct comparison between hydrodynamical simulations and observations is needed to improve the physics included in the former and test biases in the latter. Post-processing radiative transfer and synthetic observations are now the standard way to do this. We report on the first application of the \texttt{SKIRT} radiative transfer code to simulations of a star-forming cloud. The synthetic observations are then analyzed following traditional observational workflows. We find that in the early stages of the simulation, stellar radiation is inefficient in heating dust to the temperatures observed in Galactic clouds, thus the addition of an interstellar radiation field is necessary. The spectral energy distribution of the cloud settles rather quickly after $\sim3$ Myr of evolution from the onset of star formation, but its morphology continues to evolve for $\sim8$ Myr due to the expansion of \textsc{Hii} regions and the respective creation of cavities, filaments, and ridges. Modeling synthetic \textit{Herschel} fluxes with 1- or 2-component modified black bodies underestimates total dust masses by a factor of $\sim2$. Spatially-resolved fitting recovers up to about $70\%$ of the intrinsic value. This ``missing mass'' is located in a very cold dust component with temperatures below $10$ K, which does not contribute appreciably to the far-infrared flux. This effect could bias real observations if such dust exists in large amounts. Finally, we tested observational calibrations of the SFR based on infrared fluxes and concluded that they are in agreement when compared to the intrinsic SFR of the simulation averaged over $\sim100$ Myr. △ Less

Submitted 10 April, 2023; originally announced April 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.

Comments: 21 pages, 10 figures (main text). Accepted to MNRAS

Effects of dust sources on dust attenuation properties in IllustrisTNG galaxies at $z\sim 7$

Authors: Yuan-Ming Hsu, Hiroyuki Hirashita, Yen-Hsing Lin, Peter Camps, Maarten Baes

Abstract: Dust emission from high-redshift galaxies gives us a clue to the origin and evolution of dust in the early Universe. Previous studies have shown that different sources of dust (stellar dust production and dust growth in dense clouds) predict different ultraviolet (UV) extinction curves for galaxies at $z\sim 7$ but that the observed attenuation curves depend strongly on the geometry of dust and st… ▽ More Dust emission from high-redshift galaxies gives us a clue to the origin and evolution of dust in the early Universe. Previous studies have shown that different sources of dust (stellar dust production and dust growth in dense clouds) predict different ultraviolet (UV) extinction curves for galaxies at $z\sim 7$ but that the observed attenuation curves depend strongly on the geometry of dust and star distributions. Thus, we perform radiative transfer calculations under the dust-stars geometries computed by a cosmological hydrodynamic simulation (IllustrisTNG). This serves to investigate the dust attenuation curves predicted from `realistic' geometries. We choose objects with stellar mass and star formation rate appropriate for Lyman break galaxies at $z\sim 7$. We find that the attenuation curves are very different from the original extinction curves in most of the galaxies. This makes it difficult to constrain the dominant dust sources from the observed attenuation curves. We further include infrared dust emission in the analysis and plot the infrared excess (IRX)-UV spectral slope ($β$) diagram. We find that different sources of dust cause different IRX-$β$ relations for the simulated galaxies. In particular, if dust growth is the main source of dust, a variation of dust-to-metal ratio causes a more extended sequence with smaller IRX in the IRX-$β$ diagram. Thus, the comprehensive analysis of the abundances of dust and metals, the UV slope, and the dust emission could provide a clue to the dominant dust sources in the Universe. △ Less

Submitted 9 December, 2022; originally announced December 2022.

Comments: 11 pages, 6 figures, accepted for publication in MNRAS

UV to submillimetre luminosity functions of TNG50 galaxies

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

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

Submitted 12 August, 2022; originally announced August 2022.

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

Monte Carlo radiative transfer with explicit absorption to simulate absorption, scattering, and stimulated emission

Authors: Maarten Baes, Peter Camps, Kosei Matsumoto

Abstract: Context: The Monte Carlo method is probably the most widely used approach to solve the radiative transfer problem, especially in a general 3D geometry. The physical processes of emission, absorption, and scattering are easily incorporated in the Monte Carlo framework. Net stimulated emission, or absorption with a negative cross section, does not fit this method, however. Aims: We explore alteratio… ▽ More Context: The Monte Carlo method is probably the most widely used approach to solve the radiative transfer problem, especially in a general 3D geometry. The physical processes of emission, absorption, and scattering are easily incorporated in the Monte Carlo framework. Net stimulated emission, or absorption with a negative cross section, does not fit this method, however. Aims: We explore alterations to the standard photon packet life cycle in Monte Carlo radiative transfer that allow the treatment of net stimulated emission without loss of generality or efficiency. Methods: We present the explicit absorption technique that allows net stimulated emission to be handled efficiently. It uses the scattering rather than the extinction optical depth along a photon packet's path to randomly select the next interaction location, and offers a separate, deterministic treatment of absorption. We implemented the technique in a special-purpose Monte Carlo code for a two-stream 1D radiative transfer problem and in the fully featured 3D code SKIRT, and we studied its overall performance using quantitative statistical tests. Results: Our special-purpose code is capable of recovering the analytical solutions to the two-stream problem in all regimes, including the one of strong net stimulated emission. The implementation in SKIRT is straightforward, as the explicit absorption technique easily combines with the variance reduction and acceleration techniques already incorporated. In general, explicit absorption tends to improve the efficiency of the Monte Carlo routine in the regime of net absorption. Conclusions: Explicit absorption allows the treatment of net stimulated emission in Monte Carlo radiative transfer, it interfaces smoothly with... (abridged) △ Less

Submitted 12 August, 2022; originally announced August 2022.

Comments: 14 pages, accepted for publication in A&A

Journal ref: A&A 666, A101 (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.

Comments: Accepted for publication by MNRAS

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.

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

Journal ref: A&A 659, A149 (2022)

Geometry effects on dust attenuation curves with different grain sources at high redshift

Authors: Yen-Hsing Lin, Hiroyuki Hirashita, Peter Camps, Maarten Baes

Abstract: Dust has been detected in high-redshift ($z>5$) galaxies but its origin is still being debated. Dust production in high-redshift galaxies could be dominated by stellar production or by accretion (dust growth) in the interstellar medium. Previous studies have shown that these two dust sources predict different grain size distributions, which lead to significantly different extinction curves. In thi… ▽ More Dust has been detected in high-redshift ($z>5$) galaxies but its origin is still being debated. Dust production in high-redshift galaxies could be dominated by stellar production or by accretion (dust growth) in the interstellar medium. Previous studies have shown that these two dust sources predict different grain size distributions, which lead to significantly different extinction curves. In this paper, we investigate how the difference in the extinction curves affects the dust attenuation properties of galaxies by performing radiative transfer calculations. To examine the major effects of the dust--stars distribution geometry, we adopt two representative cases in spherical symmetry: the well-mixed geometry (stars and dust are homogeneously mixed) and the two-layer geometry (young stars are more concentrated in the centre). In both cases, we confirm that the attenuation curve can be drastically steepened by scattering and by different optical depths between young and old stellar populations, and can be flattened by the existence of unobscured stellar populations. We can reproduce similar attenuation curves even with very different extinction curves. Thus, we conclude that it is difficult to distinguish the dust sources only with attenuation curves. However, if we include information on dust emission and plot the IRX (infrared excess)--$β$ (ultraviolet spectral slope) relation, different dust sources predict different positions in the IRX--$β$ diagram. A larger $β$ is preferred under a similar IRX if dust growth is the dominant dust source. △ Less

Submitted 7 September, 2021; originally announced September 2021.

Comments: 11 pages, 4 figures, accepted for publication in MNRAS

Journal ref: MNRAS, 507, 2755-2765 (2021)

The structure of the Milky Way based on unWISE 3.4$μ$m integrated photometry

Authors: Aleksandr V. Mosenkov, Sergey S. Savchenko, Anton A. Smirnov, Peter Camps

Abstract: We present a detailed analysis of the Galaxy structure using an unWISE wide-field image at $3.4μ$m. We perform a 3D photometric decomposition of the Milky Way taking into account i) the projection of the Galaxy on the celestial sphere and ii) that the observer is located within the Galaxy at the solar radius. We consider a large set of photometric models starting with a pure disc model and ending… ▽ More We present a detailed analysis of the Galaxy structure using an unWISE wide-field image at $3.4μ$m. We perform a 3D photometric decomposition of the Milky Way taking into account i) the projection of the Galaxy on the celestial sphere and ii) that the observer is located within the Galaxy at the solar radius. We consider a large set of photometric models starting with a pure disc model and ending with a complex model which consists of thin and thick discs plus a boxy-peanut-shaped bulge. In our final model, we incorporate many observed features of the Milky Way, such as the disc flaring and warping, several over-densities in the plane, and the dust extinction. The model of the bulge with the corresponding X-shape structure is obtained from N-body simulations of a Milky Way-like galaxy. This allows us to retrieve the parameters of the aforementioned stellar components, estimate their contribution to the total Galaxy luminosity, and constrain the position angle of the bar. The mass of the thick disc in our models is estimated to be 0.4-1.3 of that for the thin disc. The results of our decomposition can be directly compared to those obtained for external galaxies via multicomponent photometric decomposition. △ Less

Submitted 23 August, 2021; originally announced August 2021.

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

First Light And Reionisation Epoch Simulations (FLARES) III: The properties of massive dusty galaxies at cosmic dawn

Authors: Aswin P. Vijayan, Stephen M. Wilkins, Christopher C. Lovell, Peter A. Thomas, Peter Camps, Maarten Baes, James Trayford, Jussi Kuusisto, William J. Roper

Abstract: Using the First Light And Reionisation Epoch Simulations (\textsc{Flares}) we explore the dust driven properties of massive high-redshift galaxies at $z\in[5,10]$. By post-processing the galaxy sample using the radiative transfer code \textsc{skirt} we obtain the full spectral energy distribution. We explore the resultant luminosity functions, IRX-$β$ relations as well as the luminosity-weighted d… ▽ More Using the First Light And Reionisation Epoch Simulations (\textsc{Flares}) we explore the dust driven properties of massive high-redshift galaxies at $z\in[5,10]$. By post-processing the galaxy sample using the radiative transfer code \textsc{skirt} we obtain the full spectral energy distribution. We explore the resultant luminosity functions, IRX-$β$ relations as well as the luminosity-weighted dust temperatures in the Epoch of Reionisation (EoR). We find that most of our results are in agreement with the current set of observations, but under-predict the number densities of bright IR galaxies, which are extremely biased towards the most overdense regions. We see that the \textsc{Flares} IRX-$β$ relation (for $5\le z\le8$) predominantly follows the local starburst relation. The IRX shows an increase with stellar mass, plateauing at the high-mass end ($\sim10^{10}$M$_{\odot}$) and shows no evolution in the median normalisation with redshift. We also look at the dependence of the peak dust temperature ($T_{\mathrm{peak}}$) on various galaxy properties including the stellar mass, IR luminosity and sSFR, finding the correlation to be strongest with sSFR. The luminosity-weighted dust temperatures increase towards higher redshifts, with the slope of the $T_{\mathrm{peak}}$ - redshift relation showing a higher slope than the lower redshift relations obtained from previous observational and theoretical works. The results from \textsc{Flares}, which is able to provide a better statistical sample of high-redshift galaxies compared to other simulations, provides a distinct vantage point for the high-redshift Universe. △ Less

Submitted 10 March, 2022; v1 submitted 2 August, 2021; originally announced August 2021.

Comments: 19 pages, 16 figures, Published in MNRAS. Data link: https://flaresimulations.github.io/data.html

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

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

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

Submitted 14 July, 2021; originally announced July 2021.

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

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.

Comments: 22 pages, 19 figures, accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 653, A34 (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.

Comments: Accepted for publication in The Astrophysical Journal

SpheCow: flexible dynamical models for galaxies and dark matter haloes

Authors: Maarten Baes, Peter Camps, Bert Vandenbroucke

Abstract: Simple but flexible dynamical models are useful for many purposes, including serving as the starting point for more complex models or numerical simulations of galaxies, clusters, or dark matter haloes. We present SpheCow, a new light-weight and flexible code that allows one to easily explore the structure and dynamics of any spherical model. Assuming an isotropic or Osipkov-Merritt anisotropic orb… ▽ More Simple but flexible dynamical models are useful for many purposes, including serving as the starting point for more complex models or numerical simulations of galaxies, clusters, or dark matter haloes. We present SpheCow, a new light-weight and flexible code that allows one to easily explore the structure and dynamics of any spherical model. Assuming an isotropic or Osipkov-Merritt anisotropic orbital structure, the code can automatically calculate the dynamical properties of any model with either an analytical density profile or an analytical surface density profile as starting point. We have extensively validated SpheCow using a combination of comparisons to analytical and high-precision numerical calculations, as well as the calculation of inverse formulae. SpheCow contains readily usable implementations for many standard models, including the Plummer, Hernquist, NFW, Einasto, Sérsic, and Nuker models. The code is publicly available as a set of C++ routines and as a Python module, and it is designed to be easily extendable, in the sense that new models can be added in a straightforward way. We demonstrate this by adding two new families of models in which either the density slope or the surface density slope is described by an algebraic sigmoid function. We advocate the use of the SpheCow code to investigate the full dynamical structure for models for which the distribution function cannot be expressed analytically and to explore a much wider range of models than is possible using analytical models alone. △ Less

Submitted 1 June, 2021; v1 submitted 25 May, 2021; originally announced May 2021.

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

Journal ref: A&A 652, A36 (2021)

The dynamical structure of broken power-law and double power-law models for dark matter haloes

Authors: Maarten Baes, Peter Camps

Abstract: Galaxy kinematics and gravitational lensing are two complementary ways to constrain the distribution of dark matter on galaxy scales. The typical dark matter density profiles adopted in dynamical studies cannot easily be adopted in lensing studies. Ideally, a mass model should be used that has the global characteristics of realistic dark matter distributions, and that allows for an analytical calc… ▽ More Galaxy kinematics and gravitational lensing are two complementary ways to constrain the distribution of dark matter on galaxy scales. The typical dark matter density profiles adopted in dynamical studies cannot easily be adopted in lensing studies. Ideally, a mass model should be used that has the global characteristics of realistic dark matter distributions, and that allows for an analytical calculation of the magnifications and deflection angles. A simple model with these properties, the broken-power-law (BPL) model, has very recently been introduced. We examine the dynamical structure of the family of BPL models. We derive simple closed expressions for basic dynamical properties, and study the distribution function under the assumption of velocity isotropy. We find that none of the BPL models with realistic parameters has an isotropic distribution function that is positive over the entire phase space, implying that the BPL models cannot be supported by an isotropic velocity distribution, or models with a more radially anisotropic orbital structure. This result limits the attractiveness of the BPL family as a tool for lensing studies to some degree. More generally, we find that not all members of the general family of double power-law or Zhao models, often used to model dark matter haloes, can be supported by an isotropic or radially anisotropic distribution function. In other words, the distribution function may become negative even for spherically symmetric models with a well-behaved density profile. △ Less

Submitted 2 March, 2021; originally announced March 2021.

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

CMacIonize 2.0: a novel task-based approach to Monte Carlo radiation transfer

Authors: Bert Vandenbroucke, Peter Camps

Abstract: (Context) Monte Carlo radiative transfer (MCRT) is a widely used technique to model the interaction between radiation and a medium, and plays an important role in astrophysical modelling and when comparing those models with observations. (Aims) In this work, we present a novel approach to MCRT that addresses the challenging memory access patterns of traditional MCRT algorithms, which hinder optima… ▽ More (Context) Monte Carlo radiative transfer (MCRT) is a widely used technique to model the interaction between radiation and a medium, and plays an important role in astrophysical modelling and when comparing those models with observations. (Aims) In this work, we present a novel approach to MCRT that addresses the challenging memory access patterns of traditional MCRT algorithms, which hinder optimal performance of MCRT simulations on modern hardware with a complex memory architecture. (Methods) We reformulate the MCRT photon packet life cycle as a task-based algorithm, whereby the computation is broken down into small tasks that are executed concurrently. Photon packets are stored in intermediate buffers, and tasks propagate photon packets through small parts of the computational domain, moving them from one buffer to another in the process. (Results) Using the implementation of the new algorithm in the photoionization MCRT code CMacIonize 2.0, we show that the decomposition of the MCRT grid into small parts leads to a significant performance gain during the photon packet propagation phase, which constitutes the bulk of an MCRT algorithm, as a result of better usage of memory caches. Our new algorithm is a factor 2 to 4 faster than an equivalent traditional algorithm and shows good strong scaling up to 30 threads. We briefly discuss how our new algorithm could be adjusted or extended to other astrophysical MCRT applications. (Conclusions) We show that optimising the memory access patterns of a memory-bound algorithm such as MCRT can yield significant performance gains. △ Less

Submitted 26 June, 2020; originally announced June 2020.

Comments: 18 pages, 21 figures, accepted for publication in Astronomy & Astrophysics. The source code of CMacIonize is publicly available from https://github.com/bwvdnbro/CMacIonize

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

CosTuuM: polarized thermal dust emission by magnetically oriented spheroidal grains

Authors: Bert Vandenbroucke, Maarten Baes, Peter Camps

Abstract: We present the new open source C++-based Python library CosTuuM that can be used to generate infrared absorption and emission coefficients for arbitrary mixtures of spheroidal dust grains that are (partially) aligned with a magnetic field. We outline the algorithms underlying the software, demonstrate the accuracy of our results using benchmarks from literature, and use our tool to investigate som… ▽ More We present the new open source C++-based Python library CosTuuM that can be used to generate infrared absorption and emission coefficients for arbitrary mixtures of spheroidal dust grains that are (partially) aligned with a magnetic field. We outline the algorithms underlying the software, demonstrate the accuracy of our results using benchmarks from literature, and use our tool to investigate some commonly used approximative recipes. We find that the linear polarization fraction for a partially aligned dust grain mixture can be accurately represented by an appropriate linear combination of perfectly aligned grains and grains that are randomly oriented, but that the commonly used picket fence alignment breaks down for short wavelengths. We also find that for a fixed dust grain size, the absorption coefficients and linear polarization fraction for a realistic mixture of grains with various shapes cannot both be accurately represented by a single representative grain with a fixed shape, but that instead an average over an appropriate shape distribution should be used. Insufficient knowledge of an appropriate shape distribution is the main obstacle in obtaining accurate optical properties. CosTuuM is available as a standalone Python library and can be used to generate optical properties to be used in radiative transfer applications. △ Less

Submitted 16 June, 2020; originally announced June 2020.

Comments: 25 pages, 21 figures, accepted for publication in The Astronomical Journal, CosTuuM is available from https://github.com/SKIRT/CosTuuM

Infrared luminosity functions and dust mass functions in the EAGLE simulation

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

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

Submitted 6 April, 2020; originally announced April 2020.

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

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

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

Abstract: We compare the spectral energy distributions (SEDs) and inferred physical properties for simulated and observed galaxies at low redshift. We exploit UV-submillimetre mock fluxes of ~7000 z=0 galaxies from the EAGLE suite of cosmological simulations, derived using the radiative transfer code SKIRT. We compare these to ~800 observed galaxies in the UV-submillimetre range, from the DustPedia sample o… ▽ More We compare the spectral energy distributions (SEDs) and inferred physical properties for simulated and observed galaxies at low redshift. We exploit UV-submillimetre mock fluxes of ~7000 z=0 galaxies from the EAGLE suite of cosmological simulations, derived using the radiative transfer code SKIRT. We compare these to ~800 observed galaxies in the UV-submillimetre range, from the DustPedia sample of nearby galaxies. To derive global properties, we apply the SED fitting code CIGALE consistently to both data sets, using the same set of ~80 million models. The results of this comparison reveal overall agreement between the simulations and observations, both in the SEDs and in the derived physical properties, with a number of discrepancies. The optical and far-infrared regimes, and the scaling relations based upon the global emission, diffuse dust and stellar mass, show high levels of agreement. However, the mid-infrared fluxes of the EAGLE galaxies are overestimated while the far-UV domain is not attenuated enough, compared to the observations. We attribute these discrepancies to a combination of galaxy population differences between the samples, and limitations in the subgrid treatment of star-forming regions in the EAGLE-SKIRT post-processing recipe. Our findings show the importance of detailed radiative transfer calculations and consistent comparison, and provide suggestions for improved numerical models. △ Less

Submitted 27 March, 2020; originally announced March 2020.

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

SKIRT 9: redesigning an advanced dust radiative transfer code to allow kinematics, line transfer and polarization by aligned dust grains

Authors: Peter Camps, Maarten Baes

Abstract: The open source SKIRT Monte Carlo radiative transfer code has been used for more than 15 years to model the interaction between radiation and dust in various astrophysical systems. In this work, we present version 9 of the code, which has been substantially redesigned to support long-term objectives. We invite interested readers to participate in the development, testing and application of new fea… ▽ More The open source SKIRT Monte Carlo radiative transfer code has been used for more than 15 years to model the interaction between radiation and dust in various astrophysical systems. In this work, we present version 9 of the code, which has been substantially redesigned to support long-term objectives. We invite interested readers to participate in the development, testing and application of new features such as including gas media types in addition to dust, performing line transfer in addition to continuum radiation transfer, and modeling polarization by non-spherical dust grains aligned by magnetic fields. We describe the major challenges involved in preparing the code for these and other extensions, as well as their resolution, including a completely new treatment of wavelengths to support kinematics. SKIRT 9 properly runs over 400 handcrafted functional tests and successfully performs all relevant benchmarks. The source code and all documentation is publicly available for use and ready for further collaborative development. △ Less

Submitted 2 March, 2020; originally announced March 2020.

Comments: 24 pages, 14 figures, accepted for publication in Astronomy and Computing

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

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

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

Submitted 16 January, 2020; originally announced January 2020.

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

Optical depth in polarised Monte Carlo radiative transfer

Authors: Maarten Baes, Christian Peest, Peter Camps, Ralf Siebenmorgen

Abstract: Context: The Monte Carlo method is the most widely used method to solve radiative transfer problems in astronomy, especially in a fully general 3D geometry. A crucial concept in any Monte Carlo radiative transfer code is the random generation of the next interaction location. In polarised Monte Carlo radiative transfer with aligned non-spherical grains, the nature of dichroism complicates the conc… ▽ More Context: The Monte Carlo method is the most widely used method to solve radiative transfer problems in astronomy, especially in a fully general 3D geometry. A crucial concept in any Monte Carlo radiative transfer code is the random generation of the next interaction location. In polarised Monte Carlo radiative transfer with aligned non-spherical grains, the nature of dichroism complicates the concept of optical depth. Aims: We investigate in detail the relation between optical depth and the optical properties and density of the attenuating medium in polarised Monte Carlo radiative transfer codes that take into account dichroic extinction. Methods: Based on solutions for the radiative transfer equation, we discuss the optical depth scale in polarised radiative transfer with spheroidal grains. We compare the dichroic optical depth to the extinction and total optical depth scale. Results: In a dichroic medium, the optical depth is not equal to the usual extinction optical depth, nor to the total optical depth. For representative values of the optical properties of dust grains, the dichroic optical depth can differ from the extinction or total optical depth by several ten percent. A closed expression for the dichroic optical depth cannot be given, but it can be derived efficiently through an algorithm that is based on the analytical result corresponding to elongated grains with a uniform grain alignment. Conclusions: Optical depth is more complex in dichroic media than in systems without dichroic attenuation, and this complexity needs to be considered when generating random free path lengths in Monte Carlo radiative transfer simulations. There is no benefit in using approximations instead of the dichroic optical depth. △ Less

Submitted 30 July, 2019; originally announced July 2019.

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

Journal ref: A&A 630, A61 (2019)

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

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

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

Submitted 6 May, 2019; originally announced May 2019.

Comments: 23 pages, 10 figures, Accepted to ApJ

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

The nature of sub-millimeter and highly star-forming galaxies in the EAGLE simulation

Authors: Stuart McAlpine, Ian Smail, Richard G. Bower, Mark A. Swinbank, James W. Trayford, Tom Theuns, Maarten Baes, Peter Camps, Robert A. Crain, Joop Schaye

Abstract: We exploit EAGLE, a cosmological hydrodynamical simulation, to reproduce the selection of the observed sub-millimeter (submm) galaxy population by selecting the model galaxies at $z \geq 1$ with mock submm fluxes $S_{850} \geq 1$ mJy. There is a reasonable agreement between the galaxies within this sample and the properties of the observed submm population, such as their star formation rates (SFRs… ▽ More We exploit EAGLE, a cosmological hydrodynamical simulation, to reproduce the selection of the observed sub-millimeter (submm) galaxy population by selecting the model galaxies at $z \geq 1$ with mock submm fluxes $S_{850} \geq 1$ mJy. There is a reasonable agreement between the galaxies within this sample and the properties of the observed submm population, such as their star formation rates (SFRs) at $z<3$, redshift distribution and many integrated galaxy properties. We find that the bulk of the $S_{850} \geq 1$ mJy model population is at $z = 2.5$, and that they are massive galaxies ($M_* \sim 10^{11}$ Msol) with high dust masses ($M_{\mathrm{dust}} \sim 10^{8}$ Msol), gas fractions ($f_{\mathrm{gas}} \approx 50$%) and SFRs ($\dot M_* \approx 100$ Msol/yr). They have major and minor merger fractions similar to the general population, suggesting that mergers are not the primary driver of the model submm galaxies. Instead, the $S_{850} \geq 1$ mJy model galaxies yield high SFRs primarily because they maintain a significant gas reservoir as a result of hosting an undermassive black hole. In addition, we find that not all highly star-forming EAGLE galaxies have submm fluxes $S_{850} > 1$ mJy. Thus, we investigate the nature of $z \geq 1$ highly star-forming Submm-Faint galaxies (i.e., $\dot M_* \geq 80$ Msol/yr but $S_{850}< 1$ mJy). We find they are similar to the model submm galaxies; being gas rich and hosting undermassive black holes, however they are typically lower mass ($M_* \sim 10^{10}$ Msol) and are at higher redshifts ($z>4$). These typically higher-$z$ galaxies show stronger evidence for having been triggered by major mergers, and critically, they are likely missed by current submm surveys due to their higher dust temperatures. This suggests a potentially even larger contribution to the SFR density at $z > 3$ from dust-obscured systems than implied by current observations. △ Less

Submitted 16 January, 2019; originally announced January 2019.

Comments: 15 pages, 9 figures, submitted to MNRAS

The Failure of Monte Carlo Radiative Transfer at Medium to High Optical Depths

Authors: Peter Camps, Maarten Baes

Abstract: Computer simulations of photon transport through an absorbing and/or scattering medium form an important research tool in astrophysics. Nearly all software codes performing such simulations for three-dimensional geometries employ the Monte Carlo radiative transfer method, including various forms of biasing to accelerate the calculations. Because of the probabilistic nature of the Monte Carlo techn… ▽ More Computer simulations of photon transport through an absorbing and/or scattering medium form an important research tool in astrophysics. Nearly all software codes performing such simulations for three-dimensional geometries employ the Monte Carlo radiative transfer method, including various forms of biasing to accelerate the calculations. Because of the probabilistic nature of the Monte Carlo technique, the outputs are inherently noisy, but it is often assumed that the average values provide the physically correct result. We show that this assumption is not always justified. Specifically, we study the intensity of radiation penetrating an infinite, uniform slab of material that absorbs and scatters the radiation with equal probability. The basic Monte Carlo radiative transfer method, without any biasing mechanisms, starts to break down for transverse optical depths above ~20 because so few of the simulated photon packets reach the other side of the slab. When including biasing techniques such as absorption/scattering splitting and path length stretching, the simulated photon packets do reach the other side of the slab but the biased weights do not necessarily add up to the correct solution. While the noise levels seem to be acceptable, the average values sometimes severely underestimate the correct solution. Detecting these anomalies requires the judicious application of statistical tests, similar to those used in the field of nuclear particle transport, possibly in combination with convergence tests employing consecutively larger numbers of photon packets. In any case, for transverse optical depths above ~75 the Monte Carlo methods used in our study fail to solve the one-dimensional slab problem, implying the need for approximations such as a modified random walk. △ Less

Submitted 24 May, 2018; originally announced May 2018.

Comments: Accepted for publication in the ApJ; 13 pages, 6 figures

HERschel Observations of Edge-on Spirals (HEROES). IV. Dust energy balance problem

Authors: Aleksandr V. Mosenkov, Flor Allaert, Maarten Baes, Simone Bianchi, Peter Camps, Christopher J. R. Clark, Marjorie Decleir, Gert De Geyter, Ilse De Looze, Jacopo Fritz, Gianfranco Gentile, Benne W. Holwerda, Thomas M. Hughes, Fraser Lewis, Matthew W. L. Smith, Joris Verstappen, Sam Verstocken, Sébastien Viaene

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

Submitted 24 April, 2018; originally announced April 2018.

Comments: 53 pages, 31 figures, Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 616, A120 (2018)

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

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

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

Submitted 15 December, 2017; originally announced December 2017.

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

SKIRT: hybrid parallelization of radiative transfer simulations

Authors: Sam Verstocken, Dries Van De Putte, Peter Camps, Maarten Baes

Abstract: We describe the design, implementation and performance of the new hybrid parallelization scheme in our Monte Carlo radiative transfer code SKIRT, which has been used extensively for modeling the continuum radiation of dusty astrophysical systems including late-type galaxies and dusty tori. The hybrid scheme combines distributed memory parallelization, using the standard Message Passing Interface (… ▽ More We describe the design, implementation and performance of the new hybrid parallelization scheme in our Monte Carlo radiative transfer code SKIRT, which has been used extensively for modeling the continuum radiation of dusty astrophysical systems including late-type galaxies and dusty tori. The hybrid scheme combines distributed memory parallelization, using the standard Message Passing Interface (MPI) to communicate between processes, and shared memory parallelization, providing multiple execution threads within each process to avoid duplication of data structures. The synchronization between multiple threads is accomplished through atomic operations without high-level locking (also called lock-free programming). This improves the scaling behavior of the code and substantially simplifies the implementation of the hybrid scheme. The result is an extremely flexible solution that adjusts to the number of available nodes, processors and memory, and consequently performs well on a wide variety of computing architectures. △ Less

Submitted 12 May, 2017; originally announced May 2017.

Comments: 21 pages, 20 figures

Optical colours and spectral indices of $z=0.1$ EAGLE galaxies with 3D dust radiative transfer code SKIRT

Authors: James W. Trayford, Peter Camps, Tom Theuns, Maarten Baes, Richard G. Bower, Robert A. Crain, Madusha L. P. Gunawardhana, Matthieu Schaller, Joop Schaye, Carlos S. Frenk

Abstract: We present mock optical images, broad-band and H$α$ fluxes, and D4000 spectral indices for 30,145 galaxies from the EAGLE hydrodynamical simulation at redshift $z=0.1$, modelling dust with the SKIRT Monte Carlo radiative transfer code. The modelling includes a subgrid prescription for dusty star-forming regions, with both the subgrid obscuration of these regions and the fraction of metals in diffu… ▽ More We present mock optical images, broad-band and H$α$ fluxes, and D4000 spectral indices for 30,145 galaxies from the EAGLE hydrodynamical simulation at redshift $z=0.1$, modelling dust with the SKIRT Monte Carlo radiative transfer code. The modelling includes a subgrid prescription for dusty star-forming regions, with both the subgrid obscuration of these regions and the fraction of metals in diffuse interstellar dust calibrated against far-infrared fluxes of local galaxies. The predicted optical colours as a function of stellar mass agree well with observation, with the SKIRT model showing marked improvement over a simple dust screen model. The orientation dependence of attenuation is weaker than observed because EAGLE galaxies are generally puffier than real galaxies, due to the pressure floor imposed on the interstellar medium. The mock H$α$ luminosity function agrees reasonably well with the data, and we quantify the extent to which dust obscuration affects observed H$α$ fluxes. The distribution of D4000 break values is bimodal, as observed. In the simulation, 20$\%$ of galaxies deemed `passive' for the SKIRT model, i.e. exhibiting D4000 $> 1.8$, are classified `active' when ISM dust attenuation is not included. The fraction of galaxies with stellar mass greater than $10^{10}$ M$_\odot$ that are deemed passive is slightly smaller than observed, which is due to low levels of residual star formation in these simulated galaxies. Colour images, fluxes and spectra of EAGLE galaxies are to be made available through the public EAGLE database. △ Less

Submitted 5 May, 2017; originally announced May 2017.

Comments: Accepted for publication in MNRAS, 33 pages, 21 figures, 3 Tables

TRUST I: A 3D externally illuminated slab benchmark for dust radiative transfer

Authors: K. D. Gordon, M. Baes, S. Bianchi, P. Camps, M. Juvela, R. Kuiper, T. Lunttila, K. A. Misselt, G. Natale, T. Robitaille, J. Steinacker

Abstract: The radiative transport of photons through arbitrary three-dimensional (3D) structures of dust is a challenging problem due to the anisotropic scattering of dust grains and strong coupling between different spatial regions. The radiative transfer problem in 3D is solved using Monte Carlo or Ray Tracing techniques as no full analytic solution exists for the true 3D structures. We provide the first… ▽ More The radiative transport of photons through arbitrary three-dimensional (3D) structures of dust is a challenging problem due to the anisotropic scattering of dust grains and strong coupling between different spatial regions. The radiative transfer problem in 3D is solved using Monte Carlo or Ray Tracing techniques as no full analytic solution exists for the true 3D structures. We provide the first 3D dust radiative transfer benchmark composed of a slab of dust with uniform density externally illuminated by a star. This simple 3D benchmark is explicitly formulated to provide tests of the different components of the radiative transfer problem including dust absorption, scattering, and emission. This benchmark includes models with a range of dust optical depths fully probing cases that are optically thin at all wavelengths to optically thick at most wavelengths. This benchmark includes solutions for the full dust emission including single photon (stochastic) heating as well as two simplifying approximations: One where all grains are considered in equilibrium with the radiation field and one where the emission is from a single effective grain with size-distribution-averaged properties. A total of six Monte Carlo codes and one Ray Tracing code provide solutions to this benchmark. Comparison of the results revealed that the global SEDs are consistent on average to a few percent for all but the scattered stellar flux at very high optical depths. The image results are consistent within 10%, again except for the stellar scattered flux at very high optical depths. The lack of agreement between different codes of the scattered flux at high optical depths is quantified for the first time. We provide the first 3D dust radiative transfer benchmark and validate the accuracy of this benchmark through comparisons between multiple independent codes and detailed convergence tests. △ Less

Submitted 21 April, 2017; originally announced April 2017.

Comments: 21 pages, 18 figures, Astronomy & Astrophysics, in press

Journal ref: A&A 603, A114 (2017)

Analytical expressions and numerical evaluation of the luminosity distance in a flat cosmology

Authors: Maarten Baes, Peter Camps, Dries Van De Putte

Abstract: Accurate and efficient methods to evaluate cosmological distances are an important tool in modern precision cosmology. In a flat $Λ$CDM cosmology, the luminosity distance can be expressed in terms of elliptic integrals. We derive an alternative and simple expression for the luminosity distance in a flat $Λ$CDM based on hypergeometric functions. Using a timing experiment we compare the computation… ▽ More Accurate and efficient methods to evaluate cosmological distances are an important tool in modern precision cosmology. In a flat $Λ$CDM cosmology, the luminosity distance can be expressed in terms of elliptic integrals. We derive an alternative and simple expression for the luminosity distance in a flat $Λ$CDM based on hypergeometric functions. Using a timing experiment we compare the computation time for the numerical evaluation of the various exact formulae, as well as for two approximate fitting formulae available in the literature. We find that our novel expression is the most efficient exact expression in the redshift range $z\gtrsim1$. Ideally, it can be combined with the expression based on Carlson's elliptic integrals in the range $z\lesssim1$ for high precision cosmology distance calculations over the entire redshift range. On the other hand, for practical work where relative errors of about 0.1% are acceptable, the analytical approximation proposed by Adachi & Kasai (2012) is a suitable alternative. △ Less

Submitted 28 February, 2017; originally announced February 2017.

Comments: 4 pages, 1 figure, accepted for publication in MNRAS

Polarization in Monte Carlo radiative transfer and dust scattering polarization signatures of spiral galaxies

Authors: Christian Peest, Peter Camps, Marko Stalevski, Maarten Baes, Ralf Siebenmorgen

Abstract: Polarization is an important tool to further the understanding of interstellar dust and the sources behind it. In this paper we describe our implementation of polarization that is due to scattering of light by spherical grains and electrons in the dust Monte Carlo radiative transfer code SKIRT. In contrast to the implementations of other Monte Carlo radiative transfer codes, ours uses co-moving re… ▽ More Polarization is an important tool to further the understanding of interstellar dust and the sources behind it. In this paper we describe our implementation of polarization that is due to scattering of light by spherical grains and electrons in the dust Monte Carlo radiative transfer code SKIRT. In contrast to the implementations of other Monte Carlo radiative transfer codes, ours uses co-moving reference frames that rely solely on the scattering processes. It fully supports the peel-off mechanism that is crucial for the efficient calculation of images in 3D Monte Carlo codes. We develop reproducible test cases that push the limits of our code. The results of our program are validated by comparison with analytically calculated solutions. Additionally, we compare results of our code to previously published results. We apply our method to models of dusty spiral galaxies at near-infrared and optical wavelengths. We calculate polarization degree maps and show them to contain signatures that trace characteristics of the dust arms independent of the inclination or rotation of the galaxy. △ Less

Submitted 23 February, 2017; originally announced February 2017.

Comments: 15 pages, 10 figures, accepted by Astronomy and Astrophysics

Journal ref: A&A 601, A92 (2017)

The Herschel Exploitation of Local Galaxy Andromeda (HELGA) VII: A SKIRT radiative transfer model and insights on dust heating

Authors: S. Viaene, M. Baes, A. Tamm, E. Tempel, G. Bendo, J. A. D. L. Blommaert, M. Boquien, A. Boselli, P. Camps, A. Cooray, I. De Looze, P. De Vis, J. A. Fernandez-Ontiveros, J. Fritz, M. Galametz, G. Gentile, S. Madden, M. W. L. Smith, L. Spinoglio, S. Verstocken

Abstract: The radiation of stars heats dust grains in the diffuse interstellar medium and in star-forming regions in galaxies. Modelling this interaction provides information on dust in galaxies, a vital ingredient for their evolution. It is not straightforward to identify the stellar populations heating the dust, and to link attenuation to emission on a sub-galactic scale. Radiative transfer models are abl… ▽ More The radiation of stars heats dust grains in the diffuse interstellar medium and in star-forming regions in galaxies. Modelling this interaction provides information on dust in galaxies, a vital ingredient for their evolution. It is not straightforward to identify the stellar populations heating the dust, and to link attenuation to emission on a sub-galactic scale. Radiative transfer models are able to simulate this dust-starlight interaction in a realistic, three-dimensional setting. We investigate the dust heating mechanisms on a local and global galactic scale, using the Andromeda galaxy (M31) as our laboratory. We perform a series of panchromatic radiative transfer simulations of Andromeda with our code SKIRT. The high inclination angle of M31 complicates the 3D modelling and causes projection effects. However, the observed morphology and flux density are reproduced fairly well from UV to sub-millimeter wavelengths. Our model reveals a realistic attenuation curve, compatible with previous, observational estimates. We find that the dust in M31 is mainly (91 % of the absorbed luminosity) heated by the evolved stellar populations. The bright bulge produces a strong radiation field and induces non-local heating up to the main star-forming ring at 10 kpc. The relative contribution of unevolved stellar populations to the dust heating varies strongly with wavelength and with galactocentric distance.The dust heating fraction of unevolved stellar populations correlates strongly with NUV-r colour and specific star formation rate. These two related parameters are promising probes for the dust heating sources at a local scale. △ Less

Submitted 27 September, 2016; originally announced September 2016.

Comments: 16 pages + appendix, 13 figures, accepted for publication in A&A

Journal ref: A&A 599, A64 (2017)

Radiative transfer in disc galaxies $-$ V. The accuracy of the KB approximation

Authors: Dukhang Lee, Maarten Baes, Kwang-Il Seon, Peter Camps, Sam Verstocken, Wonyong Han

Abstract: We investigate the accuracy of an approximate radiative transfer technique that was first proposed by Kylafis & Bahcall (hereafter the KB approximation) and has been popular in modelling dusty late-type galaxies. We compare realistic galaxy models calculated with the KB approximation with those of a three-dimensional Monte Carlo radiative transfer code SKIRT. The SKIRT code fully takes into accoun… ▽ More We investigate the accuracy of an approximate radiative transfer technique that was first proposed by Kylafis & Bahcall (hereafter the KB approximation) and has been popular in modelling dusty late-type galaxies. We compare realistic galaxy models calculated with the KB approximation with those of a three-dimensional Monte Carlo radiative transfer code SKIRT. The SKIRT code fully takes into account of the contribution of multiple scattering whereas the KB approximation calculates only single scattered intensity and multiple scattering components are approximated. We find that the KB approximation gives fairly accurate results if optically thin, face-on galaxies are considered. However, for highly inclined ($i \gtrsim 85^{\circ}$) and/or optically thick (central face-on optical depth $\gtrsim1$) galaxy models, the approximation can give rise to substantial errors, sometimes, up to $\gtrsim 40\%$. Moreover, it is also found that the KB approximation is not always physical, sometimes producing infinite intensities at lines of sight with high optical depth in edge-on galaxy models. There is no "simple recipe" to correct the errors of the KB approximation that is universally applicable to any galaxy models. Therefore, it is recommended that the full radiative transfer calculation be used, even though it's slower than the KB approximation. △ Less

Submitted 6 September, 2016; originally announced September 2016.

Comments: 10 pages, 6 figures, accepted for publication in MNRAS

Far-infrared and dust properties of present-day galaxies in the EAGLE simulations

Authors: Peter Camps, James W. Trayford, Maarten Baes, Tom Theuns, Matthieu Schaller, Joop Schaye

Abstract: The EAGLE cosmological simulations reproduce the observed galaxy stellar mass function and many galaxy properties. In this work, we study the dust-related properties of present-day EAGLE galaxies through mock observations in the far-infrared and submm wavelength ranges obtained with the 3D dust radiative transfer code SKIRT. To prepare an EAGLE galaxy for radiative transfer processing, we derive a… ▽ More The EAGLE cosmological simulations reproduce the observed galaxy stellar mass function and many galaxy properties. In this work, we study the dust-related properties of present-day EAGLE galaxies through mock observations in the far-infrared and submm wavelength ranges obtained with the 3D dust radiative transfer code SKIRT. To prepare an EAGLE galaxy for radiative transfer processing, we derive a diffuse dust distribution from the gas particles and we re-sample the star-forming gas particles and the youngest star particles into star-forming regions that are assigned dedicated emission templates. We select a set of redshift-zero EAGLE galaxies that matches the K-band luminosity distribution of the galaxies in the Herschel Reference Survey (HRS), a volume-limited sample of about 300 normal galaxies in the Local Universe. We find overall agreement of the EAGLE dust scaling relations with those observed in the HRS, such as the dust-to-stellar mass ratio versus stellar mass and versus NUV-r colour relations. A discrepancy in the f_250/f_350 versus f_350/f_500 submm colour-colour relation implies that part of the simulated dust is insufficiently heated, likely because of limitations in our sub-grid model for star-forming regions. We also investigate the effect of adjusting the metal-to-dust ratio and the covering factor of the photodissociation regions surrounding the star-forming cores. We are able to constrain the important dust-related parameters in our method, informing the calculation of dust attenuation for EAGLE galaxies in the UV and optical domain. △ Less

Submitted 15 July, 2016; originally announced July 2016.

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

Pinwheels in the sky, with dust: 3D modeling of the Wolf-Rayet 98a environment

Authors: Tom Hendrix, Rony Keppens, Allard Jan van Marle, Peter Camps, Maarten Baes, Zakaria Meliani

Abstract: The Wolf-Rayet 98a (WR 98a) system is a prime target for interferometric surveys, since its identification as a "rotating pinwheel nebulae", where infrared images display a spiral dust lane revolving with a 1.4 year periodicity. WR 98a hosts a WC9+OB star, and the presence of dust is puzzling given the extreme luminosities of Wolf-Rayet stars. We present 3D hydrodynamic models for WR 98a, where du… ▽ More The Wolf-Rayet 98a (WR 98a) system is a prime target for interferometric surveys, since its identification as a "rotating pinwheel nebulae", where infrared images display a spiral dust lane revolving with a 1.4 year periodicity. WR 98a hosts a WC9+OB star, and the presence of dust is puzzling given the extreme luminosities of Wolf-Rayet stars. We present 3D hydrodynamic models for WR 98a, where dust creation and redistribution are self-consistently incorporated. Our grid-adaptive simulations resolve details in the wind collision region at scales below one percent of the orbital separation (~4 AU), while simulating up to 1300 AU. We cover several orbital periods under conditions where the gas component alone behaves adiabatic, or is subject to effective radiative cooling. In the adiabatic case, mixing between stellar winds is effective in a well-defined spiral pattern, where optimal conditions for dust creation are met. When radiative cooling is incorporated, the interaction gets dominated by thermal instabilities along the wind collision region, and dust concentrates in clumps and filaments in a volume-filling fashion, so WR 98a must obey close to adiabatic evolutions to demonstrate the rotating pinwheel structure. We mimic Keck, ALMA or future E-ELT observations and confront photometric long-term monitoring. We predict an asymmetry in the dust distribution between leading and trailing edge of the spiral, show that ALMA and E-ELT would be able to detect fine-structure in the spiral indicative of Kelvin-Helmholtz development, and confirm the variation in photometry due to the orientation. Historic Keck images are reproduced, but their resolution is insufficient to detect the details we predict. △ Less

Submitted 30 May, 2016; originally announced May 2016.

Comments: Accepted for publication in mnras

HERschel Observations of Edge-on Spirals (HEROES). III. Dust energy balance study of IC 2531

Authors: Aleksandr V. Mosenkov, Flor Allaert, Maarten Baes, Simone Bianchi, Peter Camps, Gert De Geyter, Ilse De Looze, Jacopo Fritz, Gianfranco Gentile, Thomas M. Hughes, Fraser Lewis, Joris Verstappen, Sam Verstocken, Sébastien Viaene

Abstract: We investigate the dust energy balance for the edge-on galaxy IC 2531, one of the seven galaxies in the HEROES sample. We perform a state-of-the-art radiative transfer modelling based, for the first time, on a set of optical and near-infrared galaxy images. We show that taking into account near-infrared imaging in the modelling significantly improves the constraints on the retrieved parameters of… ▽ More We investigate the dust energy balance for the edge-on galaxy IC 2531, one of the seven galaxies in the HEROES sample. We perform a state-of-the-art radiative transfer modelling based, for the first time, on a set of optical and near-infrared galaxy images. We show that taking into account near-infrared imaging in the modelling significantly improves the constraints on the retrieved parameters of the dust content. We confirm the result from previous studies that including a young stellar population in the modelling is important for explaining the observed stellar energy distribution. However, the discrepancy between the observed and modelled thermal emission at far-infrared wavelengths, the so-called dust energy balance problem, is still present: the model underestimates the observed fluxes by a factor of about two. We compare two different dust models, and find that dust parameters and thus the spectral energy distribution in the infrared domain are sensitive to the adopted dust model. In general, the THEMIS model reproduces the observed emission in the infrared wavelength domain better than the popular Zubko et al. BARE-GR-S model. Our study of IC 2531 is a pilot case for detailed and uniform radiative transfer modelling of the entire HEROES sample, which will shed more light on the strength and origins of the dust energy balance problem. △ Less

Submitted 20 May, 2016; originally announced May 2016.

Comments: 15 pages, 8 figures, Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 592, A71 (2016)

Composite biasing in Monte Carlo radiative transfer

Authors: Maarten Baes, Karl D. Gordon, Tuomas Lunttila, Simone Bianchi, Peter Camps, Mika Juvela, Rolf Kuiper

Abstract: Biasing or importance sampling is a powerful technique in Monte Carlo radiative transfer, and can be applied in different forms to increase the accuracy and efficiency of simulations. One of the drawbacks of the use of biasing is the potential introduction of large weight factors. We discuss a general strategy, composite biasing, to suppress the appearance of large weight factors. We use this comp… ▽ More Biasing or importance sampling is a powerful technique in Monte Carlo radiative transfer, and can be applied in different forms to increase the accuracy and efficiency of simulations. One of the drawbacks of the use of biasing is the potential introduction of large weight factors. We discuss a general strategy, composite biasing, to suppress the appearance of large weight factors. We use this composite biasing approach for two different problems faced by current state-of-the-art Monte Carlo radiative transfer codes: the generation of photon packages from multiple components, and the penetration of radiation through high optical depth barriers. In both cases, the implementation of the relevant algorithms is trivial and does not interfere with any other optimisation techniques. Through simple test models, we demonstrate the general applicability, accuracy and efficiency of the composite biasing approach. In particular, for the penetration of high optical depths, the gain in efficiency is spectacular for the specific problems that we consider: in simulations with composite path length stretching, high accuracy results are obtained even for simulations with modest numbers of photon packages, while simulations without biasing cannot reach convergence, even with a huge number of photon packages. △ Less

Submitted 25 March, 2016; originally announced March 2016.

Comments: 12 pages, accepted for publication in A&A

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

The EAGLE simulations of galaxy formation: public release of halo and galaxy catalogues

Authors: Stuart McAlpine, John C. Helly, Matthieu Schaller, James W. Trayford, Yan Qu, Michelle Furlong, Richard G. Bower, Robert A. Crain, Joop Schaye, Tom Theuns, Claudio Dalla Vecchia, Carlos S. Frenk, Ian G. McCarthy, Adrian Jenkins, Yetli Rosas-Guevara, Simon D. M. White, Maarten Baes, Peter Camps, Gerard Lemson

Abstract: We present the public data release of halo and galaxy catalogues extracted from the EAGLE suite of cosmological hydrodynamical simulations of galaxy formation. These simulations were performed with an enhanced version of the GADGET code that includes a modified hydrodynamics solver, time-step limiter and subgrid treatments of baryonic physics, such as stellar mass loss, element-by-element radiativ… ▽ More We present the public data release of halo and galaxy catalogues extracted from the EAGLE suite of cosmological hydrodynamical simulations of galaxy formation. These simulations were performed with an enhanced version of the GADGET code that includes a modified hydrodynamics solver, time-step limiter and subgrid treatments of baryonic physics, such as stellar mass loss, element-by-element radiative cooling, star formation and feedback from star formation and black hole accretion. The simulation suite includes runs performed in volumes ranging from 25 to 100 comoving megaparsecs per side, with numerical resolution chosen to marginally resolve the Jeans mass of the gas at the star formation threshold. The free parameters of the subgrid models for feedback are calibrated to the redshift z=0 galaxy stellar mass function, galaxy sizes and black hole mass - stellar mass relation. The simulations have been shown to match a wide range of observations for present-day and higher-redshift galaxies. The raw particle data have been used to link galaxies across redshifts by creating merger trees. The indexing of the tree produces a simple way to connect a galaxy at one redshift to its progenitors at higher redshift and to identify its descendants at lower redshift. In this paper we present a relational database which we are making available for general use. A large number of properties of haloes and galaxies and their merger trees are stored in the database, including stellar masses, star formation rates, metallicities, photometric measurements and mock gri images. Complex queries can be created to explore the evolution of more than 10^5 galaxies, examples of which are provided in appendix. (abridged) △ Less

Submitted 25 August, 2016; v1 submitted 5 October, 2015; originally announced October 2015.

Comments: 12 pages, 3 figures + appendices. The simulation data is available at http://www.eaglesim.org/database.php (comments welcome). Paper accepted to AstronomyComputing. Note there is a typo in the published version, Section 4.1, replace 3/2 with 2/3 for velocity dispersion

Journal ref: Astronomy and Computing (2016), pp. 72-89

Benchmarking the Calculation of Stochastic Heating and Emissivity of Dust Grains in the Context of Radiative Transfer Simulations

Authors: Peter Camps, Karl Misselt, Simone Bianchi, Tuomas Lunttila, Christophe Pinte, Giovanni Natale, Mika Juvela, Joerg Fischera, Michael P. Fitzgerald, Karl Gordon, Maarten Baes, Juergen Steinacker

Abstract: We define an appropriate problem for benchmarking dust emissivity calculations in the context of radiative transfer (RT) simulations, specifically including the emission from stochastically heated dust grains. Our aim is to provide a self-contained guide for implementors of such functionality, and to offer insights in the effects of the various approximations and heuristics implemented by the part… ▽ More We define an appropriate problem for benchmarking dust emissivity calculations in the context of radiative transfer (RT) simulations, specifically including the emission from stochastically heated dust grains. Our aim is to provide a self-contained guide for implementors of such functionality, and to offer insights in the effects of the various approximations and heuristics implemented by the participating codes to accelerate the calculations. The benchmark problem definition includes the optical and calorimetric material properties, and the grain size distributions, for a typical astronomical dust mixture with silicate, graphite and PAH components; a series of analytically defined radiation fields to which the dust population is to be exposed; and instructions for the desired output. We process this problem using six RT codes participating in this benchmark effort, and compare the results to a reference solution computed with the publicly available dust emission code DustEM. The participating codes implement different heuristics to keep the calculation time at an acceptable level. We study the effects of these mechanisms on the calculated solutions, and report on the level of (dis)agreement between the participating codes. For all but the most extreme input fields, we find agreement within 10% across the important wavelength range from 3 to 1000 micron. We conclude that the relevant modules in RT codes can and do produce fairly consistent results for the emissivity spectra of stochastically heated dust grains. △ Less

Submitted 17 June, 2015; originally announced June 2015.

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

Journal ref: A&A 580, A87 (2015)

SKIRT: the design of a suite of input models for Monte Carlo radiative transfer simulations

Authors: Maarten Baes, Peter Camps

Abstract: The Monte Carlo method is the most popular technique to perform radiative transfer simulations in a general 3D geometry. The algorithms behind and acceleration techniques for Monte Carlo radiative transfer are discussed extensively in the literature, and many different Monte Carlo codes are publicly available. On the contrary, the design of a suite of components that can be used for the distributi… ▽ More The Monte Carlo method is the most popular technique to perform radiative transfer simulations in a general 3D geometry. The algorithms behind and acceleration techniques for Monte Carlo radiative transfer are discussed extensively in the literature, and many different Monte Carlo codes are publicly available. On the contrary, the design of a suite of components that can be used for the distribution of sources and sinks in radiative transfer codes has received very little attention. The availability of such models, with different degrees of complexity, has many benefits. For example, they can serve as toy models to test new physical ingredients, or as parameterised models for inverse radiative transfer fitting. For 3D Monte Carlo codes, this requires algorithms to efficiently generate random positions from 3D density distributions. We describe the design of a flexible suite of components for the Monte Carlo radiative transfer code SKIRT. The design is based on a combination of basic building blocks (which can be either analytical toy models or numerical models defined on grids or a set of particles) and the extensive use of decorators that combine and alter these building blocks to more complex structures. For a number of decorators, e.g. those that add spiral structure or clumpiness, we provide a detailed description of the algorithms that can be used to generate random positions. Advantages of this decorator-based design include code transparency, the avoidance of code duplication, and an increase in code maintainability. Moreover, since decorators can be chained without problems, very complex models can easily be constructed out of simple building blocks. Finally, based on a number of test simulations, we demonstrate that our design using customised random position generators is superior to a simpler design based on a generic black-box random position generator. △ Less

Submitted 28 May, 2015; originally announced May 2015.

Comments: 15 pages, 4 figures, accepted for publication in Astronomy and Computing

Dust energy balance study of two edge-on spiral galaxies in the Herschel-ATLAS survey

Authors: Gert De Geyter, Maarten Baes, Ilse De Looze, George J. Bendo, Nathan Bourne, Peter Camps, Asantha Cooray, Gianfranco De Zotti, Loretta Dunne, Simon Dye, Steve A. Eales, Jacopo Fritz, Cristina Furlanetto, Gianfranco Gentile, Thomas M. Hughes, Rob J. Ivison, Steve J. Maddox, Michał J. Michałowski, Matthew W. L. Smith, Elisabetta Valiante, Sébastien Viaene

Abstract: Interstellar dust in galaxies can be traced either through its extinction effects on the star light, or through its thermal emission at infrared wavelengths. Recent radiative transfer studies of several nearby edge-on galaxies have found an apparent inconsistency in the dust energy balance: the radiative transfer models that successfully explain the optical extinction underestimate the observed fl… ▽ More Interstellar dust in galaxies can be traced either through its extinction effects on the star light, or through its thermal emission at infrared wavelengths. Recent radiative transfer studies of several nearby edge-on galaxies have found an apparent inconsistency in the dust energy balance: the radiative transfer models that successfully explain the optical extinction underestimate the observed fluxes by an average factor of three. We investigate the dust energy balance for IC4225 and NGC5166, two edge-on spiral galaxies observed by the Herschel Space Observatory in the frame of the H-ATLAS survey. We start from models which were constrained from optical data and extend them to construct the entire spectral energy distribution of our galaxies. These predicted values are subsequently compared to the observed far-infrared fluxes. We find that including a young stellar population in the modelling is necessary as it plays a non-negligible part in the heating of the dust grains. While the modelling approach for both galaxies is nearly identical, we find two very different results. As is often seen in other edge-on spiral galaxies, the far-infrared emission of our radiative transfer model of IC4225 underestimates the observed fluxes by a factor of about three. For NGC5166 on the other hand, we find that both the predicted spectral energy distribution as well as the simulated images match the observations particularly well. We explore possible reasons for this difference and conclude that it is unlikely that one single mechanism is the cause of the dust energy balance problem in spiral galaxies. We discuss the different approaches that can be considered in order to get a conclusive answer on the origin this discrepancy. △ Less

Submitted 18 May, 2015; originally announced May 2015.

Comments: 12 pages, 4 figures, accepted for publication in MNRAS