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

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

Gas Accretion and Giant Lyman-alpha Nebulae

Several decades of observations and discoveries have shown that high-redshift AGN and massive galaxies are often surrounded by giant Lyman-alpha nebulae extending in some cases up to 500 kpc in size. In this review, I discuss the properties of the such nebulae discovered at z>2 and their connection with gas flows in and around the galaxies and their halos. In particular, I show how current observations are used to constrain the physical properties and origin of the emitting gas in terms of the Lyman-alpha photon production processes and kinematical signatures. These studies suggest that recombination radiation is the most viable scenario to explain the observed Lyman-alpha luminosities and Surface Brightness for the large majority of the nebulae and imply that a significant amount of dense, ionized and cold clumps should be present within and around the halos of massive galaxies. Spectroscopic studies suggest that, among the giant Lyman-alpha nebulae, the one associated with radio-loud AGN should have kinematics dominated by strong, ionized outflows within at least the inner 30-50 kpc. Radio-quiet nebulae instead present more quiescent kinematics compatible with stationary situation and, in some cases, suggestive of rotating structures. However, definitive evidences for accretion onto galaxies of the gas associated with the giant Lyman-alpha emission are not unambiguously detected yet. Deep surveys currently ongoing using other bright, non-resonant lines such as Hydrogen H-alpha and HeII1640 will be crucial to search for clearer signatures of cosmological gas accretion onto galaxies and AGN.Comment: Invited review to appear in Gas Accretion onto Galaxies, Astrophysics and Space Science Library, eds. A. J. Fox & R. Dave', to be published by Springe

Gravitational Instabilities in a proto-solar like disc I.: Dynamics and Chemistry

To date, most simulations of the chemistry in protoplanetary discs have used 1 + 1D or 2D axisymmetric α-disc models to determine chemical compositions within young systems. This assumption is inappropriate for non-axisymmetric, gravitationally unstable discs, which may be a significant stage in early protoplanetary disc evolution. Using 3D radiative hydrodynamics, we have modelled the physical and chemical evolution of a 0.17 M⊙ self-gravitating disc over a period of 2000 yr. The 0.8 M⊙ central protostar is likely to evolve into a solar-like star, and hence this Class 0 or early Class I young stellar object may be analogous to our early Solar system. Shocks driven by gravitational instabilities enhance the desorption rates, which dominate the changes in gas-phase fractional abundances for most species. We find that at the end of the simulation, a number of species distinctly trace the spiral structure of our relatively low-mass disc, particularly CN. We compare our simulation to that of a more massive disc, and conclude that mass differences between gravitationally unstable discs may not have a strong impact on the chemical composition. We find that over the duration of our simulation, successive shock heating has a permanent effect on the abundances of HNO, CN and NH3, which may have significant implications for both simulations and observations. We also find that HCO+ may be a useful tracer of disc mass. We conclude that gravitational instabilities induced in lower mass discs can significantly, and permanently, affect the chemical evolution, and that observations with high-resolution instruments such as Atacama Large Millimeter/submillimeter Array (ALMA) offer a promising means of characterizing gravitational instabilities in protosolar discs

Observational Diagnostics of Gas Flows: Insights from Cosmological Simulations

Galactic accretion interacts in complex ways with gaseous halos, including galactic winds. As a result, observational diagnostics typically probe a range of intertwined physical phenomena. Because of this complexity, cosmological hydrodynamic simulations have played a key role in developing observational diagnostics of galactic accretion. In this chapter, we review the status of different observational diagnostics of circumgalactic gas flows, in both absorption (galaxy pair and down-the-barrel observations in neutral hydrogen and metals; kinematic and azimuthal angle diagnostics; the cosmological column density distribution; and metallicity) and emission (Lya; UV metal lines; and diffuse X-rays). We conclude that there is no simple and robust way to identify galactic accretion in individual measurements. Rather, progress in testing galactic accretion models is likely to come from systematic, statistical comparisons of simulation predictions with observations. We discuss specific areas where progress is likely to be particularly fruitful over the next few years.Comment: Invited review to appear in Gas Accretion onto Galaxies, Astrophysics and Space Science Library, eds. A. J. Fox & R. Dave, to be published by Springer. Typos correcte

STARBENCH: the D-type expansion of an HII region

STARBENCH is a project focused on benchmarking and validating different star formation and stellar feedback codes. In this first STARBENCH paper we perform a comparison study of the D-type expansion of an H II region. The aim of this work is to understand the differences observed between the 12 participating numerical codes against the various analytical expressions examining the D-type phase of H II region expansion. To do this, we propose two well-defined tests which are tackled by 1D and 3D grid- and smoothed particle hydrodynamics-based codes. The first test examines the ‘early phase’ D-type scenario during which the mechanical pressure driving the expansion is significantly larger than the thermal pressure of the neutral medium. The second test examines the ‘late phase’ D-type scenario during which the system relaxes to pressure equilibrium with the external medium. Although they are mutually in excellent agreement, all 12 participating codes follow a modified expansion law that deviates significantly from the classical Spitzer solution in both scenarios. We present a semi-empirical formula combining the two different solutions appropriate to both early and late phases that agrees with high-resolution simulations to ≲ 2 per cent. This formula provides a much better benchmark solution for code validation than the Spitzer solution. The present comparison has validated the participating codes and through this project we provide a data set for calibrating the treatment of ionizing radiation hydrodynamics codes

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’