266 research outputs found
A scheme for radiation pressure and photon diffusion with the M1 closure in RAMSES-RT
We describe and test an updated version of radiation-hydrodynamics (RHD) in
the RAMSES code, that includes three new features: i) radiation pressure on
gas, ii) accurate treatment of radiation diffusion in an unresolved optically
thick medium, and iii) relativistic corrections that account for Doppler
effects and work done by the radiation to first order in v/c. We validate the
implementation in a series of tests, which include a morphological assessment
of the M1 closure for the Eddington tensor in an astronomically relevant
setting, dust absorption in a optically semi-thick medium, direct pressure on
gas from ionising radiation, convergence of our radiation diffusion scheme
towards resolved optical depths, correct diffusion of a radiation flash and a
constant luminosity radiation, and finally, an experiment from Davis et al. of
the competition between gravity and radiation pressure in a dusty atmosphere,
and the formation of radiative Rayleigh-Taylor instabilities. With the new
features, RAMSES-RT can be used for state-of-the-art simulations of radiation
feedback from first principles, on galactic and cosmological scales, including
not only direct radiation pressure from ionising photons, but also indirect
pressure via dust from multi-scattered IR photons reprocessed from
higher-energy radiation, both in the optically thin and thick limits.Comment: 25 pages, 13 figures, accepted for publication in MNRAS. Revised to
match published versio
A scheme for radiation pressure and photon diffusion with the M1 closure in ramses-rt
We describe and test an updated version of radiation-hydrodynamics in the ramses code, that includes three new features: (i) radiation pressure on gas, (ii) accurate treatment of radiation diffusion in an unresolved optically thick medium, and (iii) relativistic corrections that account for Doppler effects and work done by the radiation to first order in v/c. We validate the implementation in a series of tests, which include a morphological assessment of the M1 closure for the Eddington tensor in an astronomically relevant setting, dust absorption in an optically semithick medium, direct pressure on gas from ionizing radiation, convergence of our radiation diffusion scheme towards resolved optical depths, correct diffusion of a radiation flash and a constant luminosity radiation, and finally, an experiment from Davis etal. of the competition between gravity and radiation pressure in a dusty atmosphere, and the formation of radiative Rayleigh-Taylor instabilities. With the new features, ramses-rt can be used for state-of-the-art simulations of radiation feedback from first principles, on galactic and cosmological scales, including not only direct radiation pressure from ionizing photons, but also indirect pressure via dust from multiscattered IR photons reprocessed from higher-energy radiation, both in the optically thin and thick limit
The energy and dynamics of trapped radiative feedback with stellar winds
In this paper, we explore the significant, non-linear impact that stellar winds have on HâII regions. We perform a parameter study using three-dimensional radiative magnetohydrodynamic simulations of wind and ultraviolet radiation feedback from a 35 Mâ star formed self-consistently in a turbulent, self-gravitating cloud, similar to the Orion Nebula (M42) and its main ionizing source Ξ1 Ori C. Stellar winds suppress early radiative feedback by trapping ionizing radiation in the shell around the wind bubble. Rapid breakouts of warm photoionized gas (âchampagne flowsâ) still occur if the star forms close to the edge of the cloud. The impact of wind bubbles can be enhanced if we detect and remove numerical overcooling caused by shocks crossing grid cells. However, the majority of the energy in the wind bubble is still lost to turbulent mixing between the wind bubble and the gas around it. These results begin to converge if the spatial resolution at the wind bubble interface is increased by refining the grid on pressure gradients. Wind bubbles form a thin chimney close to the star, which then expands outwards as an extended plume once the wind bubble breaks out of the dense core the star formed in, allowing them to expand faster than a spherical wind bubble. We also find wind bubbles mixing completely with the photoionized gas when the HâII region breaks out of the cloud as a champagne flow, a process we term âhot champagneâ
Kiloparsec-scale Simulations of Star Formation in Disk Galaxies. IV. Regulation of Galactic Star Formation Rates by Stellar Feedback
Star formation from the interstellar medium of galactic disks is a basic process controlling the evolution of galaxies. Understanding the star formation rate in a local patch of a disk with a given gas mass is thus an important challenge for theoretical models. Here we simulate a kiloparsec region of a disk, following the evolution of self-gravitating molecular clouds down to subparsec scales, as they form stars that then inject feedback energy by dissociating and ionizing UV photons and supernova explosions. We assess the relative importance of each feedback mechanism. We find that H2-dissociating feedback results in the largest absolute reduction in star formation compared to the run with no feedback. Subsequently adding photoionization feedback produces a more modest reduction. Our fiducial models that combine all three feedback mechanisms yield, without fine-tuning, star formation rates that are in excellent agreement with observations, with H2-dissociating photons playing a crucial role. Models that only include supernova feedbackâa common method in galaxy evolution simulationsâsettle to similar star formation rates, but with very different temperature and chemical states of the gas, and with very different spatial distributions of young stars
Understanding the escape of LyC and Lyα photons from turbulent clouds
Understanding the escape of Lyman continuum (LyC) and Lyman alpha (Lya)
photons from molecular clouds is one of the keys to constraining the
reionization history of the Universe and the evolution of galaxies at high
redshift. Using a set of radiation-hydrodynamic simulations with adaptive mesh
refinement, we investigate how photons propagate and escape from turbulent
clouds with different masses, star formation efficiencies (SFEs), and
metallicities, as well as with different models of stellar spectra and
supernova feedback. We find that the escape fractions in both LyC and Lya are
generally increasing with time if the cloud is efficiently dispersed by
radiation and supernova feedback. When the total SFE is low (1% of the cloud
mass), 0.1-5% of LyC photons leave the metal-poor cloud, whereas the fractions
increase to 20-70% in clouds with a 10% SFE. LyC photons escape more
efficiently if gas metallicity is lower, if the upper mass limit in the stellar
initial mass function is higher, if binary interactions are allowed in the
evolution of stars, or if additional strong radiation pressure, such as Lya
pressure, is present. As a result, the number of escaping LyC photons can
easily vary by a factor of on cloud scales. The escape fractions of Lya
photons are systemically higher (60-80%) than those of LyC photons despite
large optical depths at line centre (). Scattering of Lya
photons is already significant on cloud scales, leading to double-peaked
profiles with peak separations of during
the initial stage of the cloud evolution, while it becomes narrower than
in the LyC bright phase. Comparisons
with observations of low-redshift galaxies suggest that Lya photons require
further interactions with neutral hydrogen to reproduce their velocity offset
for a given LyC escape fraction
Probing cosmic dawn with emission lines: predicting infrared and nebular line emission for ALMA and JWST
Infrared and nebular lines provide some of our best probes of the physics
regulating the properties of the interstellar medium (ISM) at high-redshift.
However, interpreting the physical conditions of high-redshift galaxies
directly from emission lines remains complicated due to inhomogeneities in
temperature, density, metallicity, ionisation parameter, and spectral hardness.
We present a new suite of cosmological, radiation-hydrodynamics simulations,
each centred on a massive Lyman-break galaxy that resolves such properties in
an inhomogeneous ISM. Many of the simulated systems exhibit transient but well
defined gaseous disks that appear as velocity gradients in [CII]~158.6m
emission. Spatial and spectral offsets between [CII]~158.6m and
[OIII]~88.33m are common, but not ubiquitous, as each line probes a
different phase of the ISM. These systems fall on the local [CII]-SFR relation,
consistent with newer observations that question previously observed
[CII]~158.6m deficits. Our galaxies are consistent with the nebular line
properties of observed galaxies and reproduce offsets on the BPT and
mass-excitation diagrams compared to local galaxies due to higher star
formation rate (SFR), excitation, and specific-SFR, as well as harder spectra
from young, metal-poor binaries. We predict that local calibrations between
H and [OII]~3727 luminosity and galaxy SFR apply up to , as
do the local relations between certain strong line diagnostics (R23 and
[OIII]~5007/H) and galaxy metallicity. Our new simulations are well
suited to interpret the observations of line emission from current (ALMA and
HST) and upcoming facilities (JWST and ngVLA)
New Methods for Identifying Lyman Continuum Leakers and Reionization-Epoch Analogues
Identifying low-redshift galaxies that emit Lyman continuum radiation (LyC leakers) is one of the primary, indirect methods of studying galaxy formation in the epoch of reionization. However, not only has it proved challenging to identify such systems, it also remains uncertain whether the low-redshift LyC leakers are truly âanaloguesâ of the sources that reionized the Universe. Here, we use high-resolution cosmological radiation hydrodynamics simulations to examine whether simulated galaxies in the epoch of reionization share similar emission line properties to observed LyC leakers at z ⌠3 and z ⌠0. We find that the simulated galaxies with high LyC escape fractions (fesc) often exhibit high O32 and populate the same regions of the R23âO32 plane as z ⌠3 LyC leakers. However, we show that viewing angle, metallicity, and ionization parameter can all impact where a galaxy resides on the O32âfesc plane. Based on emission line diagnostics and how they correlate with fesc, lower metallicity LyC leakers at z ⌠3 appear to be good analogues of reionization-era galaxies. In contrast, local [SâII]-deficient galaxies do not overlap with the simulated high-redshift LyC leakers on the SâII BaldwinâPhillipsâTerlevich (BPT) diagram; however, this diagnostic may still be useful for identifying leakers. We use our simulated galaxies to develop multiple new diagnostics to identify LyC leakers using infrared and nebular emission lines. We show that our model using only [CâII]158 ÎŒm and [OâIII]88 ÎŒm can identify potential leakers from non-leakers from the local Dwarf Galaxy Survey. Finally, we apply this diagnostic to known high-redshift galaxies and find that MACS 1149_JD1 at z = 9.1 is the most likely galaxy to be actively contributing to the reionization of the Universe
The SPHINX cosmological simulations of the first billion years: The impact of binary stars on reionization
We present the SPHINX suite of cosmological adaptive mesh refinement
simulations, the first radiation-hydrodynamical simulations to simultaneously
capture large-scale reionization and the escape of ionizing radiation from
thousands of resolved galaxies. Our and co-moving Mpc volumes resolve
haloes down to the atomic cooling limit and model the inter-stellar medium with
better than pc resolution. The project has numerous goals in
improving our understanding of reionization and making predictions for future
observations. In this first paper we study how the inclusion of binary stars in
computing stellar luminosities impacts reionization, compared to a model that
includes only single stars. Owing to the suppression of galaxy growth via
strong feedback, our galaxies are in good agreement with observational
estimates of the galaxy luminosity function. We find that binaries have a
significant impact on the timing of reionization: with binaries, our boxes are
percent ionized by volume at , while without them our
volumes fail to reionize by . These results are robust to changes in
volume size, resolution, and feedback efficiency. The escape of ionizing
radiation from individual galaxies varies strongly and frequently. On average,
binaries lead to escape fractions of percent, about times
higher than with single stars only. The higher escape fraction is a result of a
shallower decline in ionizing luminosity with age, and is the primary reason
for earlier reionization, although the higher integrated luminosity with
binaries also plays a sub-dominant role
- âŠ