Modelling of ionising feedback with Smoothed Particle Hydrodynamics and Monte Carlo Radiative Transfer on a Voronoi grid
arXiv:2109.13953 · doi:10.1093/mnras/stab2178
Abstract
The ionising feedback of young massive stars is well known to influence the dynamics of the birth environment and hence plays an important role in regulating the star formation process in molecular clouds. For this reason, modern hydrodynamics codes adopt a variety of techniques accounting for these radiative effects. A key problem hampering these efforts is that the hydrodynamics are often solved using smoothed particle hydrodynamics (SPH), whereas radiative transfer is typically solved on a grid. Here we present a radiation-hydrodynamics (RHD) scheme combining the SPH code Phantom and the Monte Carlo Radiative Transfer (MCRT) code CMacIonize, using the particle distribution to construct a Voronoi grid on which the MCRT is performed. We demonstrate that the scheme successfully reproduces the well-studied problem of D-type H II region expansion in a uniform density medium. Furthermore, we use this simulation setup to study the robustness of the RHD code with varying choice of grid structure, density mapping method, and mass and temporal resolution. To test the scheme under more realistic conditions, we apply it to a simulated star-forming cloud reminiscing those in the Central Molecular Zone of our galaxy, in order to estimate the amount of ionised material that a single source could create. We find that a stellar population of several is needed to noticeably ionise the cloud. Based on our results, we formulate a set of recommendations to guide the numerical setup of future and more complex simulations of star forming clouds.
arXiv comment22 pages, 21 figures, 4 tables; accepted by MNRAS (31 July 2021)
References in corpus (52)
- E pur si muove: Galiliean-invariant cosmological hydrodynamical simulations on a moving mesh
- SPLASH: An interactive visualisation tool for Smoothed Particle Hydrodynamics simulations
- Phantom: A smoothed particle hydrodynamics and magnetohydrodynamics code for astrophysics
- The Effects of Stellar Rotation. II. A Comprehensive Set of Starburst99 Models
- Modelling discontinuities and Kelvin-Helmholtz instabilities in SPH
- Massive stars in massive clusters II: Disruption of bound clusters by photoionization
- Fast and inefficient star formation due to short-lived molecular clouds and rapid feedback
- Dense gas in the Galactic central molecular zone is warm and heated by turbulence
- Dispersal of molecular clouds by ionising radiation
- Destruction of protoplanetary disks in the Orion Nebula Cluster
- Photoionising feedback in star cluster formation
- Enzo+Moray: Radiation Hydrodynamics Adaptive Mesh Refinement Simulations with Adaptive Ray Tracing
- Dynamical Expansion of Ionization and Dissociation Front around a Massive Star. II. On the Generality of Triggered Star Formation
- The thermal state of molecular clouds in the Galactic Center: evidence for non-photon-driven heating
- G0.253+0.016: a molecular cloud progenitor of an Arches-like cluster
- Adaptive Ray Tracing for Radiative Transfer around Point Sources
- Comparing molecular gas across cosmic time-scales: the Milky Way as both a typical spiral galaxy and a high-redshift galaxy analogue
- Radiation driven implosion and triggered star formation
- Dynamical HII Region Evolution in Turbulent Molecular Clouds
- Ionizing feedback from massive stars in massive clusters: Fake bubbles and untriggered star formation
- First hydrodynamics simulations of radiation forces and photoionization feedback in massive star formation
- A cluster in the making: ALMA reveals the initial conditions for high-mass cluster formation
- TRAPHIC - Radiative Transfer for Smoothed Particle Hydrodynamics Simulations
- Radiation Driven Implosion of Molecular Cloud Cores
- The modelling of feedback in star formation simulations
- StarBench: The D-type expansion of an HII region
- The Survey of Water and Ammonia in the Galactic Center (SWAG): Molecular Cloud Evolution in the Central Molecular Zone
- The dynamical evolution of molecular clouds near the Galactic Centre - II. Spatial structure and kinematics of simulated clouds
- Ionizing Radiation in Smoothed Particle Hydrodynamics
- Which feedback mechanisms dominate in the high-pressure environment of the Central Molecular Zone?
- Tracing the Conversion of Gas into Stars in Young Massive Cluster Progenitors
- SILCC-ZOOM: The early impact of ionizing radiation on forming molecular clouds
- Radiation-hydrodynamical simulations of massive star formation using Monte Carlo radiative transfer: I. Algorithms and numerical methods
- Lyman-alpha radiation hydrodynamics of galactic winds before cosmic reionization
- Photoevaporation of Jeans-unstable molecular clumps
- The dynamical evolution of molecular clouds near the Galactic Centre -- III. Tidally--induced star formation in protocluster clouds
- Star formation in 'the Brick': ALMA reveals an active proto-cluster in the Galactic centre cloud G0.253+0.016
- The fragmentation of expanding shells II: Thickness matters
- Massive star feedback in clusters: variation of the FUV interstellar radiation field in time and space
- The classical D-type expansion of spherical H II regions
- The moving mesh code Shadowfax
- Observing gas and dust in simulations of star formation with Monte Carlo radiation transport on Voronoi meshes
- Insights from Synthetic Star-forming Regions: I. Reliable Mock Observations from SPH Simulations
- TREVR: A general radiative transfer algorithm
- Fast and accurate Voronoi density gridding from Lagrangian hydrodynamics data
- Massive star formation via torus accretion: the effect of photoionization feedback
- SPAMCART: a code for smoothed particle Monte Carlo radiative transfer
- Pre-supernova feedback mechanisms drive the destruction of molecular clouds in nearby star-forming disc galaxies
- The fragmentation of expanding shells III: Oligarchic accretion and the mass spectrum of fragments
- Photoionization feedback in turbulent molecular clouds
- Native Synthetic Imaging of Smoothed Particle Hydrodynamics density fields using gridless Monte Carlo Radiative Transfer
- SWIFT: Maintaining weak-scalability with a dynamic range of in time-step size to harness extreme adaptivity