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 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 ionizing photons, but also indirect pressure via dust from multiscattered IR photons reprocessed from higher-energy radiation, both in the optically thin and thick limits.
We present a new implementation of radiation hydrodynamics (RHD) in the adaptive mesh refinement (AMR) code ramses. The multigroup radiative transfer (RT) is performed on the AMR grid with a ...first-order Godunov method using the M1 closure for the Eddington tensor, and is coupled to the hydrodynamics via non-equilibrium thermochemistry of hydrogen and helium. This moment-based approach has the great advantage that the computational cost is independent of the number of radiative sources - it can even deal with continuous regions of emission such as bound-free emission from gas. As it is built directly into ramses, the RT takes natural advantage of the refinement and parallelization strategies already in place. Since we use an explicit advection solver for the radiative transport, the time-step is restricted by the speed of light - a severe limitation that can be alleviated using the so-called reduced speed of light approximation. We propose a rigorous framework to assess the validity of this approximation in various conditions encountered in cosmology and galaxy formation. We finally perform with our newly developed code a complete suite of RHD tests, comparing our results to other RHD codes. The tests demonstrate that our code performs very well and is ideally suited for exploring the effect of radiation on current scenarios of structure and galaxy formation.
The physics of diffuse gas accretion and the properties of the cold and hot modes of accretion on to proto-galaxies between z= 2 and 5.4 is investigated using the large cosmological simulation ...performed with the ramses code on the MareNostrum supercomputing facility. Galactic winds, chemical enrichment, ultraviolet background heating and radiative cooling are taken into account in this very high resolution simulation. Using accretion-weighted temperature histograms, we have performed two different measurements of the thermal state of the gas accreted towards the central galaxy. The first measurement, performed using accretion-weighted histograms on a spherical surface of radius 0.2Rvir centred on the densest gas structure near the halo centre of mass, is a good indicator of the presence of an accretion shock in the vicinity of the galactic disc. We define the hot shock mass, Mshock, as the typical halo mass separating cold dominated from hot dominated accretion in the vicinity of the galaxy. The second measurement is performed by radially averaging histograms between 0.2Rvir and Rvir, in order to detect radially extended structures such as gas filaments: this is a good proxy for detecting cold streams feeding the central galaxy. We define Mstream as the transition mass separating cold dominated from hot dominated accretion in the outer halo, marking the disappearance of these cold streams. We find a hot shock transition mass of Mshock= 1011.6M⊙ (dark matter), with no significant evolution with redshift. Conversely, we find that Mstream increases sharply with z. Our measurements are in agreement with the analytical predictions of Birnboim & Dekel and Dekel & Birnboim, if we correct their model by assuming low metallicity(≤10−3Z⊙)for the filaments, correspondingly to our measurements. Metal enrichment of the intergalactic medium is therefore a key ingredient in determining the transition mass from cold to hot dominated diffuse gas accretion. We find that the diffuse cold gas supply at the inner halo stops at z= 2 for objects with stellar masses of about 1011.1M⊙, which is close to the quenching mass determined observationally by Bundy et al. However, its evolution with z is not well constrained, making it difficult to rule out or confirm the need for an additional feedback process such as active galactic nuclei.
Context.The hierarchical model of galaxy formation, despite its many successes, still overpredicts the baryons fraction locked in galaxies as a condensed phase. The efficiency of supernovae feedback, ...proposed a long time ago as a possible solution for this so-called “overcooling” problem, is still under debate, mainly because modelling supernovae explosions within a turbulent interstellar medium, while capturing realistic large scale flows around the galaxy is a very demanding task. Aims. Our goal is to study the effect of supernovae feedback on a disc galaxy, taking into account the impact of infalling gas on both the star formation history and the corresponding outflow structure, the apparition of a supernovae-driven wind being highly sensitive to the halo mass, the galaxy spin and the star formation efficiency. Methods. We model our galaxies as cooling and collapsing NFW spheres. The dark matter component is modelled as a static external potential, while the baryon component is described by the Euler equations using the AMR code RAMSES. Metal-dependent cooling and supernovae-heating are also implemented using state-of-the-art recipes coming from cosmological simulations. We allow for three parameters to vary: the halo circular velocity, the spin parameter and the star formation efficiency. Results. We found that the ram pressure of infalling material is the key factor limiting the apparition of galactic winds. We obtain a very low feedback efficiency, with supernovae to wind energy conversion factor around one percent, so that only low circular velocity galaxies give rise to strong winds. For massive galaxies, we obtain a galactic fountain, for which we discuss the observational properties. Conclusions. We conclude that for quiescent isolated galaxies, galactic winds appear only in very low mass systems. Although this can quite efficiently enrich the IGM with metals, they do not carry away enough cold material to solve the overcooling problem.
A large-scale hydrodynamical cosmological simulation, Horizon-AGN, is used to investigate the alignment between the spin of galaxies and the cosmic filaments above redshift 1.2. The analysis of more ...than 150 000 galaxies per time step in the redshift range 1.2 < z < 1.8 with morphological diversity shows that the spin of low-mass blue galaxies is preferentially aligned with their neighbouring filaments, while high-mass red galaxies tend to have a perpendicular spin. The reorientation of the spin of massive galaxies is provided by galaxy mergers, which are significant in their mass build-up. We find that the stellar mass transition from alignment to misalignment happens around 3 × 1010 M⊙. Galaxies form in the vorticity-rich neighbourhood of filaments, and migrate towards the nodes of the cosmic web as they convert their orbital angular momentum into spin. The signature of this process can be traced to the properties of galaxies, as measured relative to the cosmic web. We argue that a strong source of feedback such as active galactic nuclei is mandatory to quench in situ star formation in massive galaxies and promote various morphologies. It allows mergers to play their key role by reducing post-merger gas inflows and, therefore, keeping spins misaligned with cosmic filaments.
Massive galaxies in the young Universe, ten billion years ago, formed stars at surprising intensities. Although this is commonly attributed to violent mergers, the properties of many of these ...galaxies are incompatible with such events, showing gas-rich, clumpy, extended rotating disks not dominated by spheroids. Cosmological simulations and clustering theory are used to explore how these galaxies acquired their gas. Here we report that they are 'stream-fed galaxies', formed from steady, narrow, cold gas streams that penetrate the shock-heated media of massive dark matter haloes. A comparison with the observed abundance of star-forming galaxies implies that most of the input gas must rapidly convert to stars. One-third of the stream mass is in gas clumps leading to mergers of mass ratio greater than 1:10, and the rest is in smoother flows. With a merger duty cycle of 0.1, three-quarters of the galaxies forming stars at a given rate are fed by smooth streams. The rarer, submillimetre galaxies that form stars even more intensely are largely merger-induced starbursts. Unlike destructive mergers, the streams are likely to keep the rotating disk configuration intact, although turbulent and broken into giant star-forming clumps that merge into a central spheroid. This stream-driven scenario for the formation of discs and spheroids is an alternative to the merger picture.
We model the intensity of emission lines from the CO molecule, based on hydrodynamic simulations of spirals, mergers, and high-redshift galaxies with very high resolutions (3 pc and 103 M⊙) and ...detailed models for the phase-space structure of the interstellar gas including shock heating, stellar feedback processes, and galactic winds. The simulations are analyzed with a large velocity gradient (LVG) model to compute the local emission in various molecular lines in each resolution element, radiation transfer, opacity effect, and the intensity emerging from galaxies to generate synthetic spectra for various transitions of the CO molecule. This model reproduces the known properties of CO spectra and CO-to-H2 conversion factors in nearby spirals and starbursting major mergers. The high excitation of CO lines in mergers is dominated by an excess of high-density gas, and the high turbulent velocities and compression that create this dense gas excess result in broad linewidths and low CO intensity-to-H2 mass ratios. When applied to high-redshift gas-rich disks galaxies, the same model predicts that their CO-to-H2 conversion factor is almost as high as in nearby spirals, and much higher than in starbursting mergers. High-redshift disk galaxies contain giant star-forming clumps that host a high-excitation component associated to gas warmed by the spatially concentrated stellar feedback sources, although CO(1−0) to CO(3−2) emission is dominated overall by low-excitation gas around the densest clumps. These results generally highlight a strong dependence of CO excitation and the CO-to-H2 conversion factor on galaxy type, even at similar star formation rates or densities. The underlying processes are driven by the interstellar medium structure and turbulence and its response to stellar feedback, which depend on global galaxy structure and in turn affect the CO emission properties.
Context. Stars, and more particularly massive stars, have a drastic impact on galaxy evolution. Yet the conditions in which they form and collapse are still not fully understood. Aims. In particular, ...the influence of the magnetic field on the collapse of massive clumps is relatively unexplored, it is therefore of great relevance in the context of the formation of massive stars to investigate its impact. Methods. We perform high resolution, MHD simulations of the collapse of one hundred solar masses, turbulent and magnetized clouds, with the adaptive mesh refinement code RAMSES. We compute various quantities such as mass distribution, magnetic field, and angular momentum within the collapsing core and study the episodic outflows and the fragmentation that occurs during the collapse. Results. The magnetic field has a drastic impact on the cloud evolution. We find that magnetic braking is able to substantially reduce the angular momentum in the inner part of the collapsing cloud. Fast and episodic outflows are being launched with typical velocities of the order of 1−3 km s-1, although the highest velocities can be as high as 20−40 km s-1. The fragmentation in several objects is reduced in substantially magnetized clouds with respect to hydrodynamical ones by a factor of the order of 1.5−2. Conclusions. We conclude that magnetic fields have a significant impact on the evolution of massive clumps. In combination with radiation, magnetic fields largely determine the outcome of massive core collapse. We stress that numerical convergence of MHD collapse is a challenging issue. In particular, numerical diffusion appears to be important at high density and therefore could possibly lead to an overestimation of the number of fragments.
ABSTRACT
Supernova ejecta and stellar winds are believed to produce interstellar dust grains with relatively large sizes. Smaller grains can be produced via the shattering of large grains that have ...been stochastically accelerated. To understand this stochastic acceleration, we have implemented novel magnetohydrodynamic (MHD) particle-in-cell methods into the astrophysical fluid code ramses. We treat dust grains as a set of massive ‘superparticles’ that experience aerodynamic drag and Lorentz force. We subject our code to a range of numerical tests designed to validate our method in different physical conditions, as well as to illustrate possible mechanisms by which grains can be accelerated. As a final test as well as a foundation for future work, we present the results of decaying dusty MHD turbulence simulations with grain parameters chosen to resemble 1–2 μm grains in typical cold neutral medium conditions. We find that in these conditions, these grains can be effectively accelerated to well beyond their shattering velocities. This is true for both electrically charged and neutral grains. While the peak of the gas-grain relative drift velocity distribution is higher for neutral grains, the drift velocity distribution of charged grains exhibits an extended exponential tail out to much greater velocities. Even so, the shapes of the distributions are such that the extra gas-grain coupling provided by the Lorentz force offers grains relative protection from shattering. We also discuss the connection between our simulations and the relatively pristine ∼μm-sized pre-solar grains that do not appear to have undergone significant wear in their lifetimes.