Observing gas and dust in simulations of star formation with Monte Carlo radiation transport on Voronoi meshes

Abstract

Ionising feedback from massive stars dramatically affects the interstellar medium local to star forming regions. Numerical simulations are now starting to include enough complexity to produce morphologies and gas properties that are not too dissimilar from observations. The comparison between the density fields produced by hydrodynamical simulations and observations at given wavelengths relies however on photoionisation/chemistry and radiative transfer calculations. We present here an implementation of Monte Carlo radiation transport through a Voronoi tessellation in the photoionisation and dust radiative transfer code MOCASSIN. We show for the first time a synthetic spectrum and synthetic emission line maps of an hydrodynamical simulation of a molecular cloud affected by massive stellar feedback. We show that the approach on which previous work is based, which remapped hydrodynamical density fields onto Cartesian grids before performing radiative transfer/photoionisation calculations, results in significant errors in the temperature and ionisation structure of the region. Furthermore, we describe the mathematical process of tracing photon energy packets through a Voronoi tessellation, including optimisations, treating problematic cases and boundary conditions. We perform various benchmarks using both the original version of MOCASSIN and the modified version using the Voronoi tessellation. We show that for uniform grids, or equivalently a cubic lattice of cell generating points, the new Voronoi version gives the same results as the original Cartesian-grid version of MOCASSIN for all benchmarks. For non-uniform initial conditions, such as using snapshots from Smoothed Particle Hydrodynamics simulations, we show that the Voronoi version performs better than the Cartesian grid version, resulting in much better resolution in dense regions.

12 pages, 4 figures, 3 tables. Accepted for publication by MNRAS

References in corpus (7)

• Before the first supernova: combined effects of HII regions and winds on molecular clouds
• X-ray irradiated protoplanetary disk atmospheres I: Predicted emission line spectrum and photoevaporation
• The timing and location of dust formation in the remnant of SN 1987A
• X-ray enabled MOCASSIN: a 3D code for photoionized media
• Dust yields in clumpy SN shells: SN 1987A revisited
• Theoretical calculations of the HI, HeI and HeII free-bound continuou emission spectra
• The effects of spatially distributed ionisation sources on the temperature structure of HII region

Cited by in corpus (14)

• Monte Carlo Radiative Transfer
• Synthetic observations of star formation and the interstellar medium
• Modelling massive-star feedback with Monte Carlo radiation hydrodynamics: photoionization and radiation pressure in a turbulent cloud
• Arepo-MCRT: Monte Carlo Radiation Hydrodynamics on a Moving Mesh
• Insights from Synthetic Star-forming Regions: I. Reliable Mock Observations from SPH Simulations
• Fast and accurate Voronoi density gridding from Lagrangian hydrodynamics data
• The Voronoi tessellation method in astronomy
• SPAMCART: a code for smoothed particle Monte Carlo radiative transfer
• CMacIonize 2.0: a novel task-based approach to Monte Carlo radiation transfer
• Monte Carlo radiative transfer with explicit absorption to simulate absorption, scattering, and stimulated emission
• Optical depth in polarised Monte Carlo radiative transfer
• Modelling of ionising feedback with Smoothed Particle Hydrodynamics and Monte Carlo Radiative Transfer on a Voronoi grid
• Tetrahedral grids in Monte Carlo radiative transfer
• 3-D structures of planetary nebulae