Spatial field reconstruction with INLA: Application to simulated galaxies

Abstract

Aims. Monte Carlo Radiative Transfer (MCRT) simulations are a powerful tool for understanding the role of dust in astrophysical systems and its influence on observations. However, due to the strong coupling of the radiation field and medium across the whole computational domain, the problem is non-local and non-linear and such simulations are computationally expensive in case of realistic 3D inhomogeneous dust distributions. We explore a novel technique for post-processing MCRT output to reduce the total computational run time by enhancing the output of computationally less expensive simulations of lower-quality. Methods. We combine principal component analysis (PCA) and non-negative matrix factorization (NMF) as dimensionality reduction techniques together with Gaussian Markov random fields and the Integrated nested Laplace approximation (INLA), an approximate method for Bayesian inference, to detect and reconstruct the non-random spatial structure in the images of lower signal-to-noise or with missing data. Results. We test our methodology using synthetic observations of a galaxy from the SKIRT Auriga project - a suite of high resolution magneto-hydrodynamic Milky Way-sized galaxies simulated in cosmological environment by 'zoom-in' technique. With this approach, we are able to reproduce high photon number reference images $\sim 5$ times faster with median residuals below $\sim 20\%$.

To be published in Numerical methods and codes of Astronomy and Astrophysics

References in corpus (19)

• The EAGLE project: Simulating the evolution and assembly of galaxies and their environments
• Introducing the Illustris Project: Simulating the coevolution of dark and visible matter in the Universe
• The EAGLE simulations of galaxy formation: calibration of subgrid physics and model variations
• Properties of galaxies reproduced by a hydrodynamic simulation
• Simulating galaxy formation with black hole driven thermal and kinetic feedback
• Introducing the Illustris Project: the evolution of galaxy populations across cosmic time
• SKIRT: an Advanced Dust Radiative Transfer Code with a User-Friendly Architecture
• Modelling the Pan-Spectral Energy Distribution of Starburst Galaxies: IV The Controlling Parameters of the Starburst SED
• Monte Carlo Radiative Transfer
• Dissecting the active galactic nucleus in Circinus -- I. Peculiar mid-IR morphology explained by a dusty hollow cone
• Dissecting the active galactic nucleus in Circinus -- II. A thin dusty disc and a polar outflow on parsec scales
• Unsupervised star, galaxy, qso classification: Application of HDBSCAN
• Galaxy Evolution in all Five CANDELS Fields and IllustrisTNG: Morphological, Structural, and the Major Merger Evolution to $z\sim 3$
• The two formation pathways of S0 galaxies
• High-resolution synthetic UV-submm images for simulated Milky Way-type galaxies from the Auriga project
• High-resolution synthetic UV-submm images for Milky Way-mass simulated galaxies from the ARTEMIS project
• Early galaxy growth: mergers or gravitational instability?
• The dust attenuation law in $z\sim 6$ quasars
• Unmixing methods based on nonnegativity and weakly mixed pixels for astronomical hyperspectral datasets

Cited by in corpus (2)

• $R_V$ from multi-waveband galaxy polarimetry in supernovae vicinity
• EmulART: Emulating Radiative Transfer -- A pilot study on autoencoder based dimensionality reduction for radiative transfer models