skip to main content

Title: Simulating dust grain-radiation coupling on a moving mesh
ABSTRACT We present a model for the interaction between dust and radiation fields in the radiation hydrodynamic code arepo-rt, which solves the moment-based radiative transfer equations on an unstructured moving mesh. Dust is directly treated using live simulation particles, each of which represent a population of grains that are coupled to hydrodynamic motion through a drag force. We introduce methods to calculate radiation pressure on and photon absorption by dust grains. By including a direct treatment of dust, we are able to calculate dust opacities and update radiation fields self-consistently based on the local dust distribution. This hybrid scheme coupling dust particles to an unstructured mesh for radiation is validated using several test problems with known analytic solutions, including dust driven via spherically symmetric flux from a constant luminosity source and photon absorption from radiation incident on a thin layer of dust. Our methods are compatible with the multifrequency scheme in arepo-rt, which treats UV, optical photons as single scattered and IR photons as multi scattered. At IR wavelengths, we model heating of and thermal emission from dust. Dust and gas are not assumed to be in local thermodynamic equilibrium but transfer energy through collisional exchange. We estimate dust temperatures more » by balancing these dust-radiation and dust-gas energy exchange rates. This framework for coupling dust and radiation can be applied in future radiation hydrodynamic simulations of galaxy formation. « less
; ; ; ; ;
Award ID(s):
1814259 1814053 1909933 1909831
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
1344 to 1354
Sponsoring Org:
National Science Foundation
More Like this

    Hydrogen emission lines can provide extensive information about star-forming galaxies in both the local and high-redshift Universe. We present a detailed Lyman continuum (LyC), Lyman-α (Lyα), and Balmer line (Hα and Hβ) radiative transfer study of a high-resolution isolated Milky Way simulation using the state-of-the-art Arepo-RT radiation hydrodynamics code with the SMUGGLE galaxy formation model. The realistic framework includes stellar feedback, non-equilibrium thermochemistry accounting for molecular hydrogen, and dust grain evolution in the interstellar medium (ISM). We extend our publicly available Cosmic Lyα Transfer (COLT) code with photoionization equilibrium Monte Carlo radiative transfer and various methodology improvements for self-consistent end-to-end (non-)resonant line predictions. Accurate LyC reprocessing to recombination emission requires modelling pre-absorption by dust ($f_\text{abs} \approx 27.5\,\rm{per\,\,cent}$), helium ionization ($f_\text{He} \approx 8.7\,\rm{per\,\,cent}$), and anisotropic escape fractions ($f_\text{esc} \approx 7.9\,\rm{per\,\,cent}$), as these reduce the available budget for hydrogen line emission ($f_\text{H} \approx 55.9\,\rm{per\,\,cent}$). We investigate the role of the multiphase dusty ISM, disc geometry, gas kinematics, and star formation activity in governing the physics of emission and escape, focusing on the time variability, gas-phase structure, and spatial spectral, and viewing angle dependence of the emergent photons. Isolated disc simulations are well-suited for comprehensive observational comparisons with local Hα surveys, butmore »would require a proper cosmological circumgalactic medium (CGM) environment as well as less dust absorption and rotational broadening to serve as analogs for high-redshift Lyα emitting galaxies. Future applications of our framework to next-generation cosmological simulations of galaxy formation including radiation-hydrodynamics that resolve ≲10 pc multiphase ISM and ≲1 kpc CGM structures will provide crucial insights and predictions for current and upcoming Lyα observations.

    « less
  2. Abstract

    We present a formulation and numerical algorithm to extend the scheme for gray radiation magnetohydrodynamics (MHD) developed by Jiang to include the frequency dependence via the multigroup approach. The entire frequency space can be divided into an arbitrary number of groups in the lab frame, and we follow the time-dependent evolution of frequency-integrated specific intensities along discrete rays inside each group. Spatial transport of photons is done in the lab frame while all the coupling terms are solved in the fluid rest frame. Lorentz transformation is used to connect different frames. The radiation transport equation is solved fully implicitly in time while the MHD equations are evolved explicitly so that time step is not limited by the speed of light. A finite volume approach is used for transport in both spatial and frequency spaces to conserve the radiation energy density and momentum. The algorithm includes photon absorption, electron scattering, as well as Compton scattering, which is calculated by solving the Kompaneets equation. The algorithm is accurate for a wide range of optical depth conditions and can handle both radiation-pressure- and gas-pressure-dominated flows. It works for both Cartesian and curvilinear coordinate systems with adaptive mesh refinement. We provide a varietymore »of test problems including a radiating sphere, shadow test, absorption of a moving gas, Bondi-type flows, as well as a collection of test problems for thermal and bulk Compton scattering. We also discuss examples where frequency dependence can make a big difference compared with the gray approach.

    « less
  3. ABSTRACT Feedback from accreting supermassive black holes (BHs), active galactic nuclei (AGNs), is now a cornerstone of galaxy formation models. In this work, we present radiation-hydrodynamic simulations of radiative AGN feedback using the novel arepo-rt code. A central BH emits radiation at a constant luminosity and drives an outflow via radiation pressure on dust grains. Utilizing an isolated Navarro–Frenk–White (NFW) halo we validate our set-up in the single- and multiscattering regimes, with the simulated shock front propagation in excellent agreement with the expected analytic result. For a spherically symmetric NFW halo, an examination of the simulated outflow properties with radiation collimation demonstrates a decreasing mass outflow rate and momentum flux, but increasing kinetic power and outflow velocity with decreasing opening angle. We then explore the impact of a central disc galaxy and the assumed dust model on the outflow properties. The contraction of the halo during the galaxy’s formation and modelling the production of dust grains result in a factor 100 increase in the halo’s optical depth. Radiation then couples momentum more efficiently to the gas, driving a stronger shock and producing a mass-loaded $\sim \!10^{3}\, \mathrm{M}_{\odot }\, \mathrm{yr}^{-1}$ outflow with a velocity of $\sim \!2000\, \mathrm{km}\, \mathrm{s}^{-1}$. However, themore »inclusion of dust destruction mechanisms, like thermal sputtering, leads to the rapid destruction of dust grains within the outflow, reducing its properties below the initial NFW halo. We conclude that radiative AGN feedback can drive outflows, but a thorough numerical and physical treatment is required to assess its true impact.« less
  4. Abstract We present self-consistent radiation hydrodynamic simulations of hydrogen reionization performed with arepo-rt complemented by a state-of-the-art galaxy formation model. We examine how photoheating feedback, due to reionization, shapes the galaxies properties. Our fiducial model completes reionization by z ≈ 6 and matches observations of the Ly α forest, the cosmic microwave background electron scattering optical depth, the high-redshift ultraviolet (UV) luminosity function, and stellar mass function. Contrary to previous works, photoheating suppresses star formation rates by more than $50{{\ \rm per\ cent}}$ only in haloes less massive than ∼108.4 M⊙ (∼108.8 M⊙) at z = 6 (z = 5), suggesting inefficient photoheating feedback from photons within galaxies. The use of a uniform UV background that heats up the gas at z ≈ 10.7 generates an earlier onset of suppression of star formation compared to our fiducial model. This discrepancy can be mitigated by adopting a UV background model with a more realistic reionization history. In the absence of stellar feedback, photoheating alone is only able to quench haloes less massive than ∼109 M⊙ at z ≳ 5, implying that photoheating feedback is sub-dominant in regulating star formation. In addition, stellar feedback, implemented as a non-local galactic wind scheme in the simulations, weakens the strength of photoheating feedback by reducing the amountmore »of stellar sources. Most importantly, photoheating does not leave observable imprints in the UV luminosity function, stellar mass function, or the cosmic star formation rate density. The feasibility of using these observables to detect imprints of reionization therefore requires further investigation.« less

    We introduce the thesan project, a suite of large volume ($L_\mathrm{box} = 95.5 \, \mathrm{cMpc}$) radiation-magnetohydrodynamic simulations that simultaneously model the large-scale statistical properties of the intergalactic medium during reionization and the resolved characteristics of the galaxies responsible for it. The flagship simulation has dark matter and baryonic mass resolutions of $3.1 \times 10^6\, {\rm M_\odot }$ and $5.8 \times 10^5\, {\rm M_\odot }$, respectively. The gravitational forces are softened on scales of 2.2 ckpc with the smallest cell sizes reaching 10 pc at z = 5.5, enabling predictions down to the atomic cooling limit. The simulations use an efficient radiation hydrodynamics solver (arepo-rt) that precisely captures the interaction between ionizing photons and gas, coupled to well-tested galaxy formation (IllustrisTNG) and dust models to accurately predict the properties of galaxies. Through a complementary set of medium resolution simulations we investigate the changes to reionization introduced by different assumptions for ionizing escape fractions, varying dark matter models, and numerical convergence. The fiducial simulation and model variations are calibrated to produce realistic reionization histories that match the observed evolution of the global neutral hydrogen fraction and electron scattering optical depth to reionization. They also match a wealth of high-redshift observationally inferred data, including themore »stellar-to-halo-mass relation, galaxy stellar mass function, star formation rate density, and the mass–metallicity relation, despite the galaxy formation model being mainly calibrated at z = 0. We demonstrate that different reionization models give rise to varied bubble size distributions that imprint unique signatures on the 21 cm emission, especially on the slope of the power spectrum at large spatial scales, enabling current and upcoming 21 cm experiments to accurately characterize the sources that dominate the ionizing photon budget.

    « less