Source author record

Romain Teyssier

Romain Teyssier appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.CO astro-ph.GA astro-ph.IM astro-ph astro-ph.SR astro-ph.HE physics.comp-ph hep-ph physics.flu-dyn Artificial Intelligence Data Structures and Algorithms Machine Learning physics.plasm-ph

Catalog footprint

What is connected

97works

13topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

A Subgrid Turbulent Mean Field Dynamo Model for Cosmological Galaxy Formation Simulations

Magnetic fields have been included in cosmological simulations of galaxy formation only recently, due to advances in numerical techniques and a better understanding of the galaxy formation physics. In this paper, we develop a new subgrid model for the turbulent dynamo that takes place in the supersonic interstellar medium in star-forming galaxies. It is based on a mean-field approach that computes the turbulent kinetic energy at unresolved scales (the so-called subgrid scales) and modifies the induction equation to account for the corresponding $α$ dynamo. Our subgrid model depends on one free parameter, the quenching parameter, that controls the saturation of the subgrid dynamo. Thanks to this mean-field approach, we can now model the fast amplification of the magnetic field inside turbulent star-forming galaxies, without relying on artificially strong initial fields or without using prohibitively expensive high-resolution simulations. We show that the evolution of the magnetic field in our zoom-in Milky Way-like galaxy is consistent with a simple picture, in which the field is in equipartition with the turbulent kinetic energy inside the star-forming disc, with a field strength around 10 $μ$G at low redshift, while at the same time strong galactic outflows fill the halo with a slightly weaker magnetic field, whose strength (10 nG) is consistent will the ideal MHD dilution factor. Our results are in good agreement with recent theoretical and numerical predictions. We also compare our simulation with Faraday depth observations at both low and high redshift, seeing overall good agreement with some caveats. Our model naturally predicts stronger magnetic fields at high redshift (around 100 $μ$G in the galaxy and 1 $μ$G in the halo), but also stronger depolarisation effects due to stronger turbulence at early time.

preprint2022arXiv

Cosmological simulations of the same spiral galaxy: the impact of baryonic physics

The interplay of star formation and supernova (SN) feedback in galaxy formation is a key element for understanding galaxy evolution. Since these processes occur at small scales, it is necessary to have sub-grid models that recover their evolution and environmental effects at the scales reached by cosmological simulations. We simulate the same spiral galaxy inhabiting a Milky Way (MW) size halo in a cosmological environment changing the sub-grid models for SN feedback and star formation. We test combinations of the Schmidt law and a multi-freefall based star formation with delayed cooling feedback or mechanical feedback. We reach a resolution of 35 pc in a zoom-in box of 36 Mpc. For this, we use the code RAMSES with the implementation of gas turbulence in time and trace the local hydrodynamical features of the star-forming gas. Finally, we compare the galaxies at redshift 0 with global and interstellar medium observations in the MW and local spiral galaxies. The simulations show successful comparisons with observations. Nevertheless, diverse galactic morphologies are obtained from different numerical implementations. We highlight the importance of detailed modelling of the star formation and feedback processes, especially when increasing the resolution of simulations. Future improvements could alleviate the degeneracies exhibited in our simulated galaxies under different sub-grid models.

preprint2022arXiv

Cosmology with one galaxy?

Galaxies can be characterized by many internal properties such as stellar mass, gas metallicity, and star-formation rate. We quantify the amount of cosmological and astrophysical information that the internal properties of individual galaxies and their host dark matter halos contain. We train neural networks using hundreds of thousands of galaxies from 2,000 state-of-the-art hydrodynamic simulations with different cosmologies and astrophysical models of the CAMELS project to perform likelihood-free inference on the value of the cosmological and astrophysical parameters. We find that knowing the internal properties of a single galaxy allow our models to infer the value of $Ω_{\rm m}$, at fixed $Ω_{\rm b}$, with a $\sim10\%$ precision, while no constraint can be placed on $σ_8$. Our results hold for any type of galaxy, central or satellite, massive or dwarf, at all considered redshifts, $z\leq3$, and they incorporate uncertainties in astrophysics as modeled in CAMELS. However, our models are not robust to changes in subgrid physics due to the large intrinsic differences the two considered models imprint on galaxy properties. We find that the stellar mass, stellar metallicity, and maximum circular velocity are among the most important galaxy properties to determine the value of $Ω_{\rm m}$. We believe that our results can be explained taking into account that changes in the value of $Ω_{\rm m}$, or potentially $Ω_{\rm b}/Ω_{\rm m}$, affect the dark matter content of galaxies. That effect leaves a distinct signature in galaxy properties to the one induced by galactic processes. Our results suggest that the low-dimensional manifold hosting galaxy properties provides a tight direct link between cosmology and astrophysics.

preprint2022arXiv

Internal dark matter structure of the most massive galaxy clusters since redshift 1

We investigate the evolution of the dark matter density profiles of the most massive galaxy clusters in the Universe. Using a `zoom-in' procedure on a large suite of cosmological simulations of total comoving volume of $3\,(h^{-1}\,\rm Gpc)^3$, we study the 25 most massive clusters in four redshift slices from $z\sim 1$ to the present. The minimum mass is $M_{500} > 5.5 \times 10^{14}$ M$_{\odot}$ at $z=1$. Each system has more than two million particles within $r_{500}$. Once scaled to the critical density at each redshift, the dark matter profiles within $r_{500}$ are strikingly similar from $z\sim1$ to the present day, exhibiting a low dispersion of 0.15 dex, and showing little evolution with redshift in the radial logarithmic slope and scatter. They have the running power law shape typical of the NFW-type profiles, and their inner structure, resolved to $3.8\,h^{-1}$ comoving kpc at $z=1$, shows no signs of converging to an asymptotic slope. Our results suggest that this type of profile is already in place at $z>1$ in the highest-mass haloes in the Universe, and that it remains exceptionally robust to merging activity.

preprint2022arXiv

On the origin of surprisingly cold gas discs in galaxies at high redshift

We address the puzzling observational indications for very "cold" galactic discs at redshifts $z \gtrsim 3$, an epoch when discs are expected to be highly perturbed. Using a high-resolution cosmological zoom-in simulation, we identify such a cold disc at $z\sim 3.5$, with a rotation velocity to velocity dispersion ratio of $v_ϕ/σ_r \simeq 5$ for the total gas. It forms as a result of a period of intense accretion of co-planar, co-rotating gas via cold cosmic-web streams. This thin disc survives for $\sim 5$ orbital periods, after which it is disrupted by mergers and counter-rotating streams, longer but consistent with our estimate that a galaxy of this mass ($M_\star\sim10^{10}\mathrm{M_\odot}$) typically survives merger-driven spin flips for $\sim 2-3$ orbital periods. We find that $v_ϕ/σ_r$ is highly sensitive to the tracer used to perform the kinematic analysis. While it is $v_ϕ/σ_r \simeq 3.5$ for atomic HI gas, it is $v_ϕ/σ_r \simeq 8$ for molecular CO and H$_2$. This reflects the confinement of molecular gas to cold, dense clouds that reside near the disc mid-plane, while the atomic gas is spread into a turbulent and more extended thicker disc. The proposed mechanisms is a theoretical proposal that has not been validated yet with proper statistical measurements and it remains unclear whether it occurs frequently enough to explain the multiple discoveries of cold gas disks in high-z galaxies.

preprint2022arXiv

Robust field-level inference with dark matter halos

We train graph neural networks on halo catalogues from Gadget N-body simulations to perform field-level likelihood-free inference of cosmological parameters. The catalogues contain $\lesssim$5,000 halos with masses $\gtrsim 10^{10}~h^{-1}M_\odot$ in a periodic volume of $(25~h^{-1}{\rm Mpc})^3$; every halo in the catalogue is characterized by several properties such as position, mass, velocity, concentration, and maximum circular velocity. Our models, built to be permutationally, translationally, and rotationally invariant, do not impose a minimum scale on which to extract information and are able to infer the values of $Ω_{\rm m}$ and $σ_8$ with a mean relative error of $\sim6\%$, when using positions plus velocities and positions plus masses, respectively. More importantly, we find that our models are very robust: they can infer the value of $Ω_{\rm m}$ and $σ_8$ when tested using halo catalogues from thousands of N-body simulations run with five different N-body codes: Abacus, CUBEP$^3$M, Enzo, PKDGrav3, and Ramses. Surprisingly, the model trained to infer $Ω_{\rm m}$ also works when tested on thousands of state-of-the-art CAMELS hydrodynamic simulations run with four different codes and subgrid physics implementations. Using halo properties such as concentration and maximum circular velocity allow our models to extract more information, at the expense of breaking the robustness of the models. This may happen because the different N-body codes are not converged on the relevant scales corresponding to these parameters.

preprint2022arXiv

The driving mode of shock-driven turbulence

Turbulence in the interstellar medium (ISM) is crucial in the process of star formation. Shocks produced by supernova explosions, jets, radiation from massive stars, or galactic spiral-arm dynamics are amongst the most common drivers of turbulence in the ISM. However, it is not fully understood how shocks drive turbulence, in particular whether shock driving is a more solenoidal(rotational, divergence-free) or a more compressive (potential, curl-free) mode of driving turbulence. The mode of turbulence driving has profound consequences for star formation, with compressive driving producing three times larger density dispersion, and an order of magnitude higher star formation rate than solenoidal driving. Here, we use hydrodynamical simulations of a shock inducing turbulent motions in a structured, multi-phase medium. This is done in the context of a laser-induced shock, propagating into a foam material, in preparation for an experiment to be performed at the National Ignition Facility (NIF). Specifically, we analyse the density and velocity distributions in the shocked turbulent medium, and measure the turbulence driving parameter $b=(σ^{2 Γ}_{ρ/\langle ρ\rangle}-1)^{1/2} (1-σ_{ρ\langle ρ\rangle}^{-2})^{-1/2}\mathcal{M}^{-1}Γ^{-1/2}$ with the density dispersion $σ_{ρ/ \langle ρ\rangle}$, the turbulent Mach number $\mathcal{M}$, and the polytropic exponent $Γ$. Purely solenoidal and purely compressive driving correspond to $b \sim 1/3$ and $b \sim 1$, respectively. Using simulations in which a shock is driven into a multi-phase medium with structures of different sizes and $Γ< 1$, we find $b \sim 1$ for all cases, showing that shock-driven turbulence is consistent with strongly compressive driving.

preprint2022arXiv

The Nature of High ${\rm [OIII]}_{\rm 88μm}$/${\rm [CII]}_{\rm 158μm}$ Galaxies in the Epoch of Reionization: Low Carbon Abundance and a Top-Heavy IMF?

ALMA observations of $z>6$ galaxies have revealed abnormally high [OIII]$_{\rm 88μm}$/[CII]$_{\rm 158μm}$ ratios and [CII]$_{\rm 158μm}$ deficits compared to local galaxies. Numerous solutions have been proposed including differences in C and O abundance ratios, observational bias, and differences in ISM properties. In order to elucidate the underlying physics that drives this high-redshift phenomenon, we employ SPHINX$^{20}$, a state-of-the-art, cosmological radiation-hydrodynamics simulation, that resolves detailed ISM properties of thousands of galaxies in the epoch of reionization. We find that the observed $z>6$ [OIII]-SFR and [CII]-SFR relations can only be reproduced when the C/O abundance ratio is $\sim8\times$ lower than Solar and the total metal production is $\sim4\times$ higher than that of a Kroupa IMF. This implies that high-redshift galaxies are potentially primarily enriched by low-metallicity core-collapse supernovae with a more top-heavy IMF. As AGB stars and type-Ia supernova begin to contribute to the galaxy metallicity, both the [CII]-SFR and [CII] luminosity functions are predicted to converge to observed values at $z\sim4.5$. While we demonstrate that ionisation parameter, LyC escape fraction, ISM gas density, and CMB attenuation all drive galaxies towards higher [OIII]/[CII], observed values at $z>6$ can only be reproduced with substantially lower C/O abundances compared to Solar. The combination of [CII] and [OIII] can be used to predict the values of ionisation parameter, ISM gas density, and LyC escape fraction and we provide estimates of these quantities for nine observed $z>6$ galaxies. Finally, we demonstrate that [OI]$_{\rm 63μm}$ can be used as a replacement for [CII] in high-redshift galaxies where [CII] is unobserved and argue that more observation time should be used to target [OI] at $z>6$. (Abridged)

preprint2022arXiv

The short ionizing photon mean free path at z=6 in Cosmic Dawn III, a new fully-coupled radiation-hydrodynamical simulation of the Epoch of Reionization

Recent determinations of the mean free path of ionising photons (mfp) in the intergalactic medium (IGM) at $\rm z=6$ are lower than many theoretical predictions. To gain insight into this issue, we investigate the evolution of the mfp in our new massive fully coupled radiation hydrodynamics cosmological simulation of reionization: Cosmic Dawn III (CoDaIII). CoDaIII's scale ($\rm 94^3 \, cMpc^3$) and resolution ($\rm 8192^3$ grid) make it particularly suitable to study the evolution of the IGM during Reionization. The simulation was performed with RAMSES-CUDATON on Summit, and used 131072 processors coupled to 24576 GPUs, making it the largest Reionization simulation, and largest RAMSES simulation ever performed. A superior agreement with global constraints on Reionization is obtained in CoDaIII over CoDaII especially for the evolution of the neutral hydrogen fraction and the cosmic photo-ionization rate, thanks to an improved calibration, later end of reionization ($\rm z=5.6$), and higher spatial resolution. Analyzing the mfp, we find that CoDaIII reproduces the most recent observations very well, from $\rm z=6$ to $\rm z=4.6$. We show that the distribution of the mfp in CoDaIII is bimodal, with short (neutral) and long (ionized) mfp modes, respectively, due to the patchiness of reionization and the co-existence of neutral versus ionized regions during Reionization. The neutral mode peaks at sub-kpc to kpc scales of mfp, while the ionized mode peak evolves from $\rm 0.1 Mpc/h$ at $\rm z=7$ to $\sim 10$ Mpc/h at $\rm z=5.2$. Computing the mfp as the average of the ionized mode provides the best match to the recent observational determinations. The distribution reduces to a single neutral (ionized) mode at $\rm z>13$ ($\rm z<5$).

preprint2022arXiv

Toward fully compressible numerical simulations of stellar magneto-convection with the RAMSES code

Numerical simulations of magneto-convection have greatly expanded our understanding of stellar interiors and stellar magnetism. Recently, fully compressible hydrodynamical simulations of full-star models have demonstrated the feasibility of studying the excitation and propagation of pressure and internal gravity waves in stellar interiors, which would allow for a direct comparison with asteroseismological measurements. However, the impact of magnetic fields on such waves has not been taken into account yet in three-dimensional simulations. We conduct a proof of concept for the realization of three-dimensional, fully compressible, magneto-hydrodynamical numerical simulations of stellar interiors with the RAMSES code. We adapted the RAMSES code to deal with highly subsonic turbulence, typical of stellar convection, by implementing a well-balanced scheme in the numerical solver. We then ran and analyzed three-dimensional hydrodynamical and magneto-hydrodynamical simulations with different resolutions of a plane-parallel convective envelope on a Cartesian grid. Both hydrodynamical and magneto-hydrodynamical simulations develop a quasi-steady, turbulent convection layer from random density perturbations introduced over the initial profiles. The convective flows are characterized by small-amplitude fluctuations around the hydrodynamical equilibrium of the stellar interior, which is preserved over the whole simulation time. Using our compressible well-balanced scheme, we were able to model flows with Mach numbers as low as $\mathcal{M} \sim 10^{-3}$, but even lower Mach number flows are possible in principle. In the magneto-hydrodynamical runs, we observe an exponential growth of magnetic energy consistent with the action of a small-scale dynamo. (Abridged)

preprint2022arXiv

Turbulence Generation by Shock Interaction with a Highly Non-Uniform Medium

An initially planar shock wave propagating into a medium of non-uniform density will be perturbed, leading to the generation of post-shock velocity perturbations. Using numerical simulations we study this phenomenon in the case of highly-non-uniform density (order-unity normalized variance, $σ_ρ/\overlineρ \sim 1$) and strong shocks (shock Mach numbers $\overline{M}_s \gtrsim 10$). This leads to a highly disrupted shock and a turbulent post-shock flow. We simulate this interaction for a range of shock drives and initial density configurations meant to mimic those which might be presently achieved in experiments. Theoretical considerations lead to scaling relations, which are found to reasonably predict the post-shock turbulence properties. The turbulent velocity dispersion and turbulent Mach number are found to depend on the pre-shock density dispersion and shock speed in a manner consistent with the linear Richtymer-Meshkov instability prediction. We also show a dependence of the turbulence generation on the scale of density perturbations. The post-shock pressure and density, which can be substantially reduced relative to the unperturbed case, are found to be reasonably predicted by a simplified analysis that treats the extended shock transition region as a single normal shock.

preprint2021arXiv

Forming early-type galaxies without AGN feedback: a combination of merger-driven outflows and inefficient star formation

Regulating the available gas mass inside galaxies proceeds through a delicate balance between inflows and outflows, but also through the internal depletion of gas due to star formation. At the same time, stellar feedback is the internal engine that powers the strong outflows. Since star formation and stellar feedback are both small scale phenomena, we need a realistic and predictive subgrid model for both. We describe the implementation of supernova momentum feedback and star formation based on the turbulence of the gas in the RAMSES code. For star formation, we adopt the so-called multi-freefall model. The resulting star formation efficiencies can be significantly smaller or bigger than the traditionally chosen value of $1\%$. We apply these new numerical models to a prototype cosmological simulation of a massive halo that features a major merger which results in the formation of an early-type galaxy without using AGN feedback. We find that the feedback model provides the first order mechanism for regulating the stellar and baryonic content in our simulated galaxy. At high redshift, the merger event pushes gas to large densities and large turbulent velocity dispersions, such that efficiencies come close to $10\%$, resulting in large SFR. We find small molecular gas depletion time during the starburst, in perfect agreement with observations. Furthermore, at late times, the galaxy becomes quiescent with efficiencies significantly smaller than $1\%$, resulting in small SFR and long molecular gas depletion time.

preprint2021arXiv

Parameter inference with non-linear galaxy clustering: accounting for theoretical uncertainties

We implement EuclidEmulator (version 1), an emulator for the non-linear correction of the matter power spectrum, into the MCMC forecasting code MontePython. We compare the performance of Halofit, HMCode, and EuclidEmulator1, both at the level of power spectrum prediction and at the level of posterior probability distributions of the cosmological parameters, for different cosmological models and different galaxy power spectrum wave number cut-offs. We confirm that the choice of the power spectrum predictor has a non-negligible effect on the computed sensitivities when doing cosmological parameter forecasting, even for a conservative wave number cut-off of $0.2\,h\,{\rm Mpc}^{-1}$. We find that EuclidEmulator1 is on average up to $17\%$ more sensitive to the cosmological parameters than the other two codes, with the most significant improvements being for the Hubble parameter of up to $42\%$ and the equation of state of dark energy of up to $26\%$, depending on the case. In addition, we point out that the choice of the power spectrum predictor contributes to the risk of computing a significantly biased mean cosmology when doing parameter estimations. For the four tested scenarios we find biases, averaged over the cosmological parameters, of between 0.5 and 2$σ$ (from below $1σ$ up to $6σ$ for individual parameters). This paper provides a proof of concept that this risk can be mitigated by taking a well-tailored theoretical uncertainty into account as this allows to reduce the bias by a factor of 2 to 5, depending on the case under consideration, while keeping posterior credibility contours small: the standard deviations are amplified by a factor of $\leq1.4$ in all cases.

preprint2020arXiv

Baryonic effects for weak lensing. Part I. Power spectrum and covariance matrix

Baryonic feedback effects lead to a suppression of the weak lensing angular power spectrum on small scales. The poorly constrained shape and amplitude of this suppression is an important source of uncertainties for upcoming cosmological weak lensing surveys such as Euclid or LSST. In this first paper in a series of two, we use simulations to build a Euclid-like tomographic mock data-set for the cosmic shear power spectrum and the corresponding covariance matrix, which are both corrected for baryonic effects following the baryonification method of Schneider et al. (2019). In addition, we develop an emulator to obtain fast predictions of the baryonic power suppression, allowing us to perform a likelihood inference analysis for a standard $Λ$CDM cosmology with both cosmological and astrophysical parameters. Our main findings are the following: (i) ignoring baryonic effects leads to a greater than 5$σ$ bias on the cosmological parameters $Ω_m$ and $σ_8$; (ii) restricting the analysis to the largest scales, that are mostly unaffected by baryons, makes the bias disappear, but results in a blow-up of the $Ω_m$-$σ_8$ contour area by more than a factor of 10; (iii) ignoring baryonic effects on the covariance matrix does not significantly affect cosmological parameter estimates; (iv) while the baryonic suppression is mildly cosmology dependent, this effect does not noticeably modify the posterior contours. Overall, we conclude that including baryonic uncertainties in terms of nuisance parameters results in unbiased and surprisingly tight constraints on cosmology.

preprint2020arXiv

Baryonic effects for weak lensing. Part II. Combination with X-ray data and extended cosmologies

An accurate modelling of baryonic feedback effects is required to exploit the full potential of future weak-lensing surveys such as Euclid or LSST. In this second paper in a series of two, we combine Euclid-like mock data of the cosmic shear power spectrum with an eROSITA X-ray mock of the cluster gas fraction to run a combined likelihood analysis including both cosmological and baryonic parameters. Following the first paper of this series, the baryonic effects (based on the baryonic correction model of Schneider et al. 2019) are included in both the tomographic power spectrum and the covariance matrix. However, this time we assume the more realistic case of a $Λ$CDM cosmology with massive neutrinos, and we consider several extensions of the currently favoured cosmological model. For the standard $Λ$CDM case, we show that including X-ray data reduces the uncertainties on the sum of the neutrino mass by $\sim30$ percent, while there is only a mild improvement on other parameters such as $Ω_m$ and $σ_8$. As extensions of $Λ$CDM, we consider the cases of a dynamical dark energy model (wCDM), a $f(R)$ gravity model (fRCDM), and a mixed dark matter model ($Λ$MDM) with both a cold and a warm/hot dark matter component. We find that combining weak lensing with X-ray data only leads to a mild improvement of the constraints on the additional parameters of wCDM, while the improvement is more substantial for both fRCDM and $Λ$MDM. Ignoring baryonic effects in the analysis pipeline leads to significant false-detections of either phantom dark energy or a light subdominant dark matter component. Overall we conclude that for all cosmologies considered, a general parametrisation of baryonic effects is both necessary and sufficient to obtain tight constraints on cosmological parameters.

preprint2020arXiv

Galactic ionising photon budget during the Epoch of Reionisation in the Cosmic Dawn II simulation

Cosmic Dawn ("CoDa") II yields the first statistically-meaningful determination of the relative contribution to reionization by galaxies of different halo mass, from a fully-coupled radiation-hydrodynamics simulation of the epoch of reionization large enough ($\sim$ 100 Mpc) to model global reionization while resolving the formation of all galactic halos above $\sim 10^8 M_\odot$. Cell transmission inside high-mass haloes is bi-modal -- ionized cells are transparent, while neutral cells absorb the photons their stars produce - and the halo escape fraction $f_{esc}$ reflects the balance of star formation rate ("SFR") between these modes. The latter is increasingly prevalent at higher halo mass, driving down $f_{esc}$ (we provide analytical fits to our results), whereas halo escape luminosity, proportional to $f_{esc} \times$SFR, increases with mass. Haloes with dark matter masses within $6.10^{8} M_\odot < M_h < 3.10^{10} M_\odot$ produce $\sim 80$% of the escaping photons at z=7, when the Universe is 50% ionized, making them the main drivers of cosmic reionization. Less massive haloes, though more numerous, have low SFRs and contribute less than 10% of the photon budget then, despite their high $f_{esc}$. High mass haloes are too few and too opaque, contributing $<10$% despite their high SFRs. The dominant mass range is lower (higher) at higher (lower) redshift, as mass function and reionization advance together (e.g. at z$=8.5$, x$_{\rm HI}=0.9$, $M_h < 5.10^9 M_\odot$ haloes contributed $\sim$80%). Galaxies with UV magnitudes $M_{AB1600}$ between $-12$ and $-19$ dominated reionization between z$=6$ and 8.

preprint2020arXiv

Rapid Filamentary Accretion as the Origin of Extended Thin Discs

Galactic outflows driven by stellar feedback are crucial for explaining the inefficiency of star formation in galaxies. Although strong feedback can promote the formation of galactic discs by limiting star formation at early times and removing low angular momentum gas, it is not understood how the same feedback can result in diverse objects such as elliptical galaxies or razor thin spiral galaxies. We investigate this problem using cosmological zoom-in simulations of two galaxies forming within $10^{12}~\mathrm{M_\odot}$ halos with almost identical mass accretion histories and halo spin parameters. However, the two resulting galaxies end up with very different bulge-to-disc ratios at $z = 0$. At $z>1.5$, the two galaxies feature a surface density of star formation $Σ_{\rm SFR}\simeq 10~\mathrm{M_\odot}~{\rm yr}^{-1}~{\rm kpc}^{-2}$, leading to strong outflows. After the last starburst episode, both galaxies feature a dramatic gaseous disc growth from 1~kpc to 5~kpc during 1~Gyr, a decisive event we dub "the Grand Twirl". After this event, the evolutionary tracks diverge strongly, with one galaxy ending up as a bulge-dominated galaxy, whereas the other ends up as a disc-dominated galaxy. The origins of this dichotomy are the angular momentum of the accreted gas, and whether it adds constructively to the initial disc angular momentum. The build-up of this extended disc leads to a rapid lowering of $Σ_{\rm SFR}$ by over two orders of magnitude with $Σ_{\rm SFR} \lesssim 0.1~\mathrm{M_\odot}~{\rm yr}^{-1}~{\rm kpc}^{-2}$, in remarkable agreement with what is derived from Milky Way stellar populations. As a consequence, supernovae explosions are spread out and cannot launch galactic outflows anymore, allowing for the persistence of a thin, gently star forming, extended disc.

preprint2020arXiv

The AGORA high-resolution galaxy simulations comparison project: Public data release

As part of the AGORA High-resolution Galaxy Simulations Comparison Project (Kim et al. 2014, 2016) we have generated a suite of isolated Milky Way-mass galaxy simulations using 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical galaxy formation community. In these simulations we adopted identical galactic disk initial conditions, and common physics models (e.g., radiative cooling and ultraviolet background by a standardized package). Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production were carefully constrained. Here we release the simulation data to be freely used by the community. In this release we include the disk snapshots at 0 and 500Myr of evolution per each code as used in Kim et al. (2016), from simulations with and without star formation and feedback. We encourage any member of the numerical galaxy formation community to make use of these resources for their research - for example, compare their own simulations with the AGORA galaxies, with the common analysis yt scripts used to obtain the plots shown in our papers, also available in this release.

preprint2019arXiv

On the origin of the peak of the stellar initial mass function: exploring the tidal screening theory

Classical theories for the stellar initial mass function (IMF) predict a peak mass which scales with the properties of the molecular cloud. In this work, we explore a new theory proposed by Lee & Hennebelle (2018). The idea is that the tidal field around first Larson cores prevents the formation of other collapsing clumps within a certain radius. The protostar can then freely accrete the gas within this radius. This leads to a peak mass of roughly $10 \, M_{\mathrm{1LC}}$, independent of the parent cloud properties. Using simple analytical arguments, we derive a collapse condition for clumps located close to a protostar. We then study the tidal field and the corresponding collapse condition using a series of numerical simulations. We find that the tidal field around protostars is indeed strong enough to prevent clumps from collapsing unless they have high enough densities. For each newly formed protostar, we determine the region in which tidal screening is dominant. We call this the tidal bubble. The mass within this bubble is our estimate for the final mass of the star. Using this formalism, we are able to construct a very good prediction for the final IMF in our simulations. Not only do we correctly predict the peak, but we are also able to reproduced the high and low mass end of the IMF. We conclude that tidal forces are important in determining the final mass of a star and might be the dominant effect in setting the peak mass of the IMF.

preprint2016arXiv

A new method to quantify the effects of baryons on the matter power spectrum

Future large-scale galaxy surveys have the potential to become leading probes for cosmology provided the influence of baryons on the total mass distribution is understood well enough. As hydrodynamical simulations strongly depend on details in the feedback implementations, no unique and robust predictions for baryonic effects currently exist. In this paper we propose a baryonic correction model that modifies the density field of dark-matter-only $N$-body simulations to mimic the effects of baryons from any underlying adopted feedback recipe. The model assumes haloes to consist of 4 components: 1- hot gas in hydrostatical equilibrium, 2- ejected gas from feedback processes, 3- central galaxy stars, and 4- adiabatically relaxed dark matter, which all modify the initial dark-matter-only density profiles. These altered profiles allow to define a displacement field for particles in $N$-body simulations and to modify the total density field accordingly. The main advantage of the baryonic correction model is to connect the total matter density field to the observable distribution of gas and stars in haloes, making it possible to parametrise baryonic effects on the matter power spectrum. We show that the most crucial quantities are the mass fraction of ejected gas and its corresponding ejection radius. The former controls how strongly baryons suppress the power spectrum, while the latter provides a measure of the scale where baryonic effects become important. A comparison with X-ray and Sunyaev-Zel'dovich cluster observations suggests that baryons suppress wave modes above $k\sim0.5$ h/Mpc with a maximum suppression of 10-25 percent around $k\sim 2$ h/Mpc. More detailed observations of the gas in the outskirts of groups and clusters are required to decrease the large uncertainties of these numbers.

preprint2016arXiv

Cosmic Dawn (CoDa): the First Radiation-Hydrodynamics Simulation of Reionization and Galaxy Formation in the Local Universe

Cosmic reionization by starlight from early galaxies affected their evolution, thereby impacting reionization, itself. Star formation suppression, for example, may explain the observed underabundance of Local Group dwarfs relative to N-body predictions for Cold Dark Matter. Reionization modelling requires simulating volumes large enough [~(100Mpc)^3] to sample reionization "patchiness", while resolving millions of galaxy sources above ~10^8 Msun , combining gravitational and gas dynamics with radiative transfer. Modelling the Local Group requires initial cosmological density fluctuations pre-selected to form the well-known structures of the local universe today. Cosmic Dawn ("CoDa") is the first such fully-coupled, radiation-hydrodynamics simulation of reionization of the local universe. Our new hybrid CPU-GPU code, RAMSES-CUDATON, performs hundreds of radiative transfer and ionization rate-solver timesteps on the GPUs for each hydro-gravity timestep on the CPUs. CoDa simulated (91Mpc)^3 with 4096^3 particles and cells, to redshift 4.23, on ORNL supercomputer Titan, utilizing 8192 cores and 8192 GPUs. Global reionization ended slightly later than observed. However, a simple temporal rescaling which brings the evolution of ionized fraction into agreement with observations also reconciles ionizing flux density, cosmic star formation history, CMB electron scattering optical depth and galaxy UV luminosity function with their observed values. Photoionization heating suppressed the star formation of haloes below ~2 x 10^9 Msun , decreasing the abun- dance of faint galaxies around MAB_1600 = [-10,-12]. For most of reionization, star formation was dominated by haloes between 10^10 - 10^11 Msun , so low-mass halo suppression was not reflected by a distinct feature in the global star formation history. (Abridged)

preprint2016arXiv

Matter power spectrum and the challenge of percent accuracy

Future galaxy surveys require one percent precision in the theoretical knowledge of the power spectrum over a large range including very nonlinear scales. While this level of accuracy is easily obtained in the linear regime with perturbation theory, it represents a serious challenge for small scales where numerical simulations are required. In this paper we quantify the precision of present-day $N$-body methods, identifying main potential error sources from the set-up of initial conditions to the measurement of the final power spectrum. We directly compare three widely used $N$-body codes, Ramses, Pkdgrav3, and Gadget3 which represent three main discretisation techniques: the particle-mesh method, the tree method, and a hybrid combination of the two. For standard run parameters, the codes agree to within one percent at $k\leq1$ $h\,\rm Mpc^{-1}$ and to within three percent at $k\leq10$ $h\,\rm Mpc^{-1}$. We also consider the bispectrum and show that the reduced bispectra agree at the sub-percent level for $k\leq 2$ $h\,\rm Mpc^{-1}$. In a second step, we quantify potential errors due to initial conditions, box size, and resolution using an extended suite of simulations performed with our fastest code Pkdgrav3. We demonstrate that the simulation box size should not be smaller than $L=0.5$ $h^{-1}\rm Gpc$ to avoid systematic finite-volume effects (while much larger boxes are required to beat down the statistical sample variance). Furthermore, a maximum particle mass of $M_{\rm p}=10^{9}$ $h^{-1}\rm M_{\odot}$ is required to conservatively obtain one percent precision of the matter power spectrum. As a consequence, numerical simulations covering large survey volumes of upcoming missions such as DES, LSST, and Euclid will need more than a trillion particles to reproduce clustering properties at the targeted accuracy.

preprint2016arXiv

nIFTY galaxy cluster simulations III: The Similarity & Diversity of Galaxies & Subhaloes

We examine subhaloes and galaxies residing in a simulated LCDM galaxy cluster ($M^{\rm crit}_{200}=1.1\times10^{15}M_\odot/h$) produced by hydrodynamical codes ranging from classic Smooth Particle Hydrodynamics (SPH), newer SPH codes, adaptive and moving mesh codes. These codes use subgrid models to capture galaxy formation physics. We compare how well these codes reproduce the same subhaloes/galaxies in gravity only, non-radiative hydrodynamics and full feedback physics runs by looking at the overall subhalo/galaxy distribution and on an individual objects basis. We find the subhalo population is reproduced to within $\lesssim10\%$ for both dark matter only and non-radiative runs, with individual objects showing code-to-code scatter of $\lesssim0.1$ dex, although the gas in non-radiative simulations shows significant scatter. Including feedback physics significantly increases the diversity. Subhalo mass and $V_{max}$ distributions vary by $\approx20\%$. The galaxy populations also show striking code-to-code variations. Although the Tully-Fisher relation is similar in almost all codes, the number of galaxies with $10^{9}M_\odot/h\lesssim M_*\lesssim 10^{12}M_\odot/h$ can differ by a factor of 4. Individual galaxies show code-to-code scatter of $\sim0.5$ dex in stellar mass. Moreover, strong systematic differences exist, with some codes producing galaxies $70\%$ smaller than others. The diversity partially arises from the inclusion/absence of AGN feedback. Our results combined with our companion papers demonstrate that subgrid physics is not just subject to fine-tuning, but the complexity of building galaxies in all environments remains a challenge. We argue even basic galaxy properties, such as the stellar mass to halo mass, should be treated with errors bars of $\sim0.2-0.4$ dex.

preprint2016arXiv

nIFTy Galaxy Cluster simulations IV: Quantifying the Influence of Baryons on Halo Properties

Building on the initial results of the nIFTy simulated galaxy cluster comparison, we compare and contrast the impact of baryonic physics with a single massive galaxy cluster, run with 11 state-of-the-art codes, spanning adaptive mesh, moving mesh, classic and modern SPH approaches. For each code represented we have a dark matter only (DM) and non-radiative (NR) version of the cluster, as well as a full physics (FP) version for a subset of the codes. We compare both radial mass and kinematic profiles, as well as global measures of the cluster (e.g. concentration, spin, shape), in the NR and FP runs with that in the DM runs. Our analysis reveals good consistency (<= 20%) between global properties of the cluster predicted by different codes when integrated quantities are measured within the virial radius R200. However, we see larger differences for quantities within R2500, especially in the FP runs. The radial profiles reveal a diversity, especially in the cluster centre, between the NR runs, which can be understood straightforwardly from the division of codes into classic SPH and non-classic SPH (including the modern SPH, adaptive and moving mesh codes); and between the FP runs, which can also be understood broadly from the division of codes into those that include AGN feedback and those that do not. The variation with respect to the median is much larger in the FP runs with different baryonic physics prescriptions than in the NR runs with different hydrodynamics solvers.

preprint2016arXiv

nIFTy galaxy cluster simulations V: Investigation of the Cluster Infall Region

We examine the properties of the galaxies and dark matter haloes residing in the cluster infall region surrounding the simulated $Λ$CDM galaxy cluster studied by Elahi et al. (2016) at z=0. The $1.1\times10^{15}h^{-1}\text{M}_{\odot}$ galaxy cluster has been simulated with eight different hydrodynamical codes containing a variety of hydrodynamic solvers and subgrid schemes. All models completed a dark-matter only, non-radiative and full-physics run from the same initial conditions. The simulations contain dark matter and gas with mass resolution $m_{\text{DM}}=9.01\times 10^8h^{-1}\text{M}_{\odot}$ and $m_{\text{gas}}=1.9\times 10^8h^{-1}\text{M}_{\odot}$ respectively. We find that the synthetic cluster is surrounded by clear filamentary structures that contain ~60% of haloes in the infall region with mass ~$10^{12.5} - 10^{14} h^{-1}\text{M}_{\odot}$, including 2-3 group-sized haloes ($> 10^{13}h^{-1}\text{M}_{\odot}$). However, we find that only ~10% of objects in the infall region are subhaloes residing in haloes, which may suggest that there is not much ongoing preprocessing occurring in the infall region at z=0. By examining the baryonic content contained within the haloes, we also show that the code-to-code scatter in stellar fraction across all halo masses is typically ~2 orders of magnitude between the two most extreme cases, and this is predominantly due to the differences in subgrid schemes and calibration procedures that each model uses. Models that do not include AGN feedback typically produce too high stellar fractions compared to observations by at least ~1 order of magnitude.

preprint2016arXiv

PKDGRAV3: Beyond Trillion Particle Cosmological Simulations for the Next Era of Galaxy Surveys

We report on the successful completion of a 2 trillion particle cosmological simulation to z=0 run on the Piz Daint supercomputer (CSCS, Switzerland), using 4000+ GPU nodes for a little less than 80h of wall-clock time or 350,000 node hours. Using multiple benchmarks and performance measurements on the US Oak Ridge National Laboratory Titan supercomputer, we demonstrate that our code PKDGRAV3, delivers, to our knowledge, the fastest time-to-solution for large-scale cosmological N-body simulations. This was made possible by using the Fast Multipole Method in conjunction with individual and adaptive particle time steps, both deployed efficiently (and for the first time) on supercomputers with GPU-accelerated nodes. The very low memory footprint of PKDGRAV3 allowed us to run the first ever benchmark with 8 trillion particles on Titan, and to achieve perfect scaling up to 18000 nodes and a peak performance of 10 Pflops.

preprint2016arXiv

RHAPSODY-G simulations II - Baryonic growth and metal enrichment in massive galaxy clusters

We study the evolution of the stellar component and the metallicity of both the intracluster medium and of stars in massive ($M_{\rm vir}\approx 6\times 10^{14}$ M$_{\odot}/h$) simulated galaxy clusters from the Rhapsody-G suite in detail and compare them to observational results. The simulations were performed with the AMR code RAMSES and include the effect of AGN feedback at the sub-grid level. AGN feedback is required to produce realistic galaxy and cluster properties and plays a role in mixing material in the central regions and regulating star formation in the central galaxy. In both our low and high resolution runs with fiducial stellar yields, we find that stellar and ICM metallicities are a factor of two lower than in observations. We find that cool core clusters exhibit steeper metallicity gradients than non-cool core clusters, in qualitative agreement with observations. We verify that the ICM metallicities measured in the simulation can be explained by a simple "regulator" model in which the metallicity is set by a balance of stellar yield and gas accretion. It is plausible that a combination of higher resolution and higher metal yield in AMR simulation would allow the metallicity of simulated clusters to match observed values; however this hypothesis needs to be tested with future simulations. Comparison to recent literature highlights that results concerning the metallicity of clusters and cluster galaxies might depend sensitively on the scheme chosen to solve the hydrodynamics.

preprint2015arXiv

Black hole evolution: I. Supernova-regulated black hole growth

The growth of a supermassive black hole (BH) is determined by how much gas the host galaxy is able to feed it, which in turn is controlled by the cosmic environment, through galaxy mergers and accretion of cosmic flows that time how galaxies obtain their gas, but also by internal processes in the galaxy, such as star formation and feedback from stars and the BH itself. In this paper, we study the growth of a 10^12 Msun halo at z=2, which is the progenitor of al group of galaxies at z=0, and of its central BH by means of a high-resolution zoomed cosmological simulation, the Seth simulation. We study the evolution of the BH driven by the accretion of cold gas in the galaxy, and explore the efficiency of the feedback from supernovae (SNe). For a relatively inefficient energy input from SNe, the BH grows at the Eddington rate from early times, and reaches self-regulation once it is massive enough. We find that at early cosmic times z>3.5, efficient feedback from SNe forbids the formation of a settled disc as well as the accumulation of dense cold gas in the vicinity of the BH and starves the central compact object. As the galaxy and its halo accumulate mass, they become able to confine the nuclear inflows provided by major mergers and the BH grows at a sustained near-to-Eddington accretion rate. We argue that this mechanism should be ubiquitous amongst low-mass galaxies, corresponding to galaxies with a stellar mass below <10^9 Msun in our simulations.

preprint2015arXiv

Distribution of streaming rates into high-redshift galaxies

We study the accretion along streams from the cosmic web into high-redshift massive galaxies using three sets of AMR hydro-cosmological simulations. We find that the streams keep a roughly constant accretion rate as they penetrate into the halo centre. The mean accretion rate follows the mass and redshift dependence predicted for haloes by the EPS approximation, dM / dt is proportional to Mvir^{1.25} (1 + z)^{2.5}. The distribution of the accretion rates can well be described by a sum of two Gaussians, the primary corresponding to "smooth inflow" and the secondary to "mergers". The same functional form was already found for the distributions of specific star formation rates in observations. The mass fraction in the smooth component is 60 - 90 %, insensitive to redshift or halo mass. The simulations with strong feedback show clear signs of reaccretion due to recycling of galactic winds. The mean accretion rate for the mergers is a factor 2 - 3 larger than that of the smooth component. The standard deviation of the accretion rate is 0.2 - 0.3 dex, showing no trend with mass or redshift. For the smooth component it is 0.12 - 0.24 dex.

preprint2015arXiv

Galaxies that Shine: radiation-hydrodynamical simulations of disk galaxies

Radiation feedback is typically implemented using subgrid recipes in hydrodynamical simulations of galaxies. Very little work has so far been performed using radiation-hydrodynamics (RHD), and there is no consensus on the importance of radiation feedback in galaxy evolution. We present RHD simulations of isolated galaxy disks of different masses with a resolution of 18 pc. Besides accounting for supernova feedback, our simulations are the first galaxy-scale simulations to include RHD treatments of photo-ionisation heating and radiation pressure, from both direct optical/UV radiation and multi-scattered, re-processed infrared (IR) radiation. Photo-heating smooths and thickens the disks and suppresses star formation about as much as the inclusion of ("thermal dump") supernova feedback does. These effects decrease with galaxy mass and are mainly due to the prevention of the formation of dense clouds, as opposed to their destruction. Radiation pressure, whether from direct or IR radiation, has little effect, but for the IR radiation we show that its impact is limited by our inability to resolve the high optical depths for which multi-scattering becomes important. While artificially boosting the IR optical depths does reduce the star formation, it does so by smoothing the gas rather than by generating stronger outflows. We conclude that although higher-resolution simulations, and potentially also different supernova implementations, are needed for confirmation, our findings suggest that radiation feedback is more gentle and less effective than is often assumed in subgrid prescriptions.

preprint2015arXiv

Galaxy Evolution: Modeling the Role of Non-thermal Pressure in the Interstellar medium

Galaxy evolution depends strongly on the physics of the interstellar medium (ISM). Motivated by the need to incorporate the properties of the ISM in cosmological simulations we construct a simple method to include the contribution of non-thermal components in the calculation of pressure of interstellar gas. In our method we treat three non-thermal components - turbulence, magnetic fields and cosmic rays - and effectively parametrize their amplitude. We assume that the three components settle into a quasi-steady-state that is governed by the star formation rate, and calibrate their magnitude and density dependence by the observed Radio-FIR correlation, relating synchrotron radiation to star formation rates of galaxies. We implement our model in single cell numerical simulation of a parcel of gas with constant pressure boundary conditions and demonstrate its effect and potential. Then, the non-thermal pressure model is incorporated into RAMSES and hydrodynamic simulations of isolated galaxies with and without the non-thermal pressure model are presented and studied. Specifically, we demonstrate that the inclusion of realistic non-thermal pressure reduces the star formation rate by an order of magnitude and increases the gas depletion time by as much. We conclude that the non-thermal pressure can prolong the star formation epoch and achieve consistency with observations without invoking artificially strong stellar feedback.

preprint2015arXiv

Milking the spherical cow: on aspherical dynamics in spherical coordinates

Galaxies and the dark matter halos that host them are not spherically symmetric, yet spherical symmetry is a helpful simplifying approximation for idealised calculations and analysis of observational data. The assumption leads to an exact conservation of angular momentum for every particle, making the dynamics unrealistic. But how much does that inaccuracy matter in practice for analyses of stellar distribution functions, collisionless relaxation, or dark matter core-creation? We provide a general answer to this question for a wide class of aspherical systems; specifically, we consider distribution functions that are "maximally stable", i.e. that do not evolve at first order when external potentials (which arise from baryons, large scale tidal fields or infalling substructure) are applied. We show that a spherically-symmetric analysis of such systems gives rise to the false conclusion that the density of particles in phase space is ergodic (a function of energy alone). Using this idea we are able to demonstrate that: (a) observational analyses that falsely assume spherical symmetry are made more accurate by imposing a strong prior preference for near-isotropic velocity dispersions in the centre of spheroids; (b) numerical simulations that use an idealised spherically-symmetric setup can yield misleading results and should be avoided where possible; and (c) triaxial dark matter halos (formed in collisionless cosmological simulations) nearly attain our maximally-stable limit, but their evolution freezes out before reaching it.

preprint2015arXiv

nIFTy galaxy cluster simulations II: radiative models

We have simulated the formation of a massive galaxy cluster (M$_{200}^{\rm crit}$ = 1.1$\times$10$^{15}h^{-1}M_{\odot}$) in a $Λ$CDM universe using 10 different codes (RAMSES, 2 incarnations of AREPO and 7 of GADGET), modeling hydrodynamics with full radiative subgrid physics. These codes include Smoothed-Particle Hydrodynamics (SPH), spanning traditional and advanced SPH schemes, adaptive mesh and moving mesh codes. Our goal is to study the consistency between simulated clusters modeled with different radiative physical implementations - such as cooling, star formation and AGN feedback. We compare images of the cluster at $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. We find that, with respect to non-radiative simulations, dark matter is more centrally concentrated, the extent not simply depending on the presence/absence of AGN feedback. The scatter in global quantities is substantially higher than for non-radiative runs. Intriguingly, adding radiative physics seems to have washed away the marked code-based differences present in the entropy profile seen for non-radiative simulations in Sembolini et al. (2015): radiative physics + classic SPH can produce entropy cores. Furthermore, the inclusion/absence of AGN feedback is not the dividing line -as in the case of describing the stellar content- for whether a code produces an unrealistic temperature inversion and a falling central entropy profile. However, AGN feedback does strongly affect the overall stellar distribution, limiting the effect of overcooling and reducing sensibly the stellar fraction.

preprint2015arXiv

PHEW: a parallel segmentation algorithm for three-dimensional AMR datasets - application to structure detection in self-gravitating flows

We introduce PHEW (Parallel HiErarchical Watershed), a new segmentation algorithm to detect structures in astrophysical fluid simulations, and its implementation into the adaptive mesh refinement (AMR) code ramses. PHEW works on the density field defined on the adaptive mesh, and can thus be used on the gas density or the dark matter density after a projection of the particles onto the grid. The algorithm is based on a "watershed" segmentation of the computational volume into dense regions, followed by a merging of the segmented patches based on the saddle point topology of the density field. PHEW is capable of automatically detecting connected regions above the adopted density threshold, as well as the entire set of substructures within. Our algorithm is fully parallel and uses the MPI library. We describe in great detail the parallel algorithm and perform a scaling experiment which proves the capability of phew to run efficiently on massively parallel systems.

preprint2015arXiv

Rhapsody-G simulations: galaxy clusters as baryonic closed boxes and the covariance between hot gas and galaxies

Within a sufficiently large cosmic volume, conservation of baryons implies a simple `closed box' view in which the sum of the baryonic components must equal a constant fraction of the total enclosed mass. We present evidence from Rhapsody-G hydrodynamic simulations of massive galaxy clusters that the closed-box expectation may hold to a surprising degree within the interior, non-linear regions of haloes. At a fixed halo mass, we find a significant anti-correlation between hot gas mass fraction and galaxy mass fraction (cold gas + stars), with a rank correlation coefficient of -0.69 within $R_{500c}$. Because of this anti-correlation, the total baryon mass serves as a low-scatter proxy for total cluster mass. The fractional scatter of total baryon fraction scales approximately as $0.02 (Δ_c/100)^{0.6}$, while the scatter of either gas mass or stellar mass is larger in magnitude and declines more slowly with increasing radius. We discuss potential observational tests using cluster samples selected by optical and hot gas properties; the simulations suggest that joint selection on stellar and hot gas has potential to achieve 5% scatter in total halo mass.

preprint2014arXiv

A fast, robust, and simple implicit method for adaptive time-stepping on adaptive mesh-refinement grids

Implicit solvers present strong limitations when used on supercomputing facilities and in particular for adaptive mesh-refinement codes. We present a new method for implicit adaptive time-stepping on adaptive mesh refinement-grids. We implement it in the radiation hydrodynamics solver we designed for the RAMSES code for astrophysical purposes and, more particularly, for protostellar collapse. We briefly recall the radiation hydrodynamics equations and the adaptive time-stepping methodology used for hydrodynamical solvers. We then introduce the different types of boundary conditions (Dirichlet, Neumann, and Robin) that are used at the interface between levels and present our implementation of the new method in the RAMSES code. The method is tested against classical diffusion and radiation hydrodynamics tests, after which we present an application for protostellar collapse. We show that using Dirichlet boundary conditions at level interfaces is a good compromise between robustness and accuracy and that it can be used in structure formation calculations. The gain in computational time over our former unique time step method ranges from factors of 5 to 50 depending on the level of adaptive time-stepping and on the problem. We successfully compare the old and new methods for protostellar collapse calculations that involve highly nonlinear physics. We have developed a simple but robust method for adaptive time-stepping of implicit scheme on adaptive mesh-refinement grids. It can be applied to a wide variety of physical problems that involve diffusion processes.

preprint2014arXiv

A Systematic Look at the Effects of Radiative Feedback on Disc Galaxy Formation

Galaxy formation models and simulations rely on various feedback mechanisms to reproduce the observed baryonic scaling relations and galaxy morphologies. Although dwarf galaxy and giant elliptical properties can be explained using feedback from supernova and active galactic nuclei, Milky Way-sized galaxies still represent a challenge to current theories of galaxy formation. In this paper, we explore the possible role of feedback from stellar radiation in regulating the main properties of disk galaxies such as our own Milky Way. We have performed a suite of cosmological simulations of the same $\sim10^{12} {\rm M}_{\odot}$ halo selected based on its rather typical mass accretion history. We have implemented radiative feedback from young stars using a crude model of radiative transfer for ultraviolet (UV) and infrared (IR) radiation. However, the model is realistic enough such that the dust opacity plays a direct role in regulating the efficiency of our feedback mechanism. We have explored various models for the dust opacity, assuming different constant dust temperatures, as well as a varying dust temperature model. We find that while strong radiative feedback appears as a viable mechanism to regulate the stellar mass fraction in massive galaxies, it also prevents the formation of discs with reasonable morphologies. In models with strong stellar radiation feedback, stellar discs are systematically too thick while the gas disc morphology is completely destroyed due to the efficient mixing between the feedback-affected gas and its surroundings. At the resolution of our simulation suite, we find it impossible to preserve spiral disc morphology while at the same time expelling enough baryons to satisfy the abundance matching constraints.

preprint2014arXiv

Baryonic and dark matter distribution in cosmological simulations of spiral galaxies

We study three cosmological hydrodynamical simulations of Milky Way(MW)-sized halos including a comparison with the dark matter(DM)-only counterparts. We find one of our simulated galaxies with interesting MW-like features. Thanks to a consistently tuned star formation rate and supernovae feedback we obtain an extended disk and a flat rotation curve with a satisfying circular velocity and a reasonable DM density in the solar neighbourhood. Mimicking observational methods, we re-derive the stellar mass and obtain stellar-to-halo mass ratios reduced by more than 50\%. We show the interaction between the baryons and the dark matter which is first contracted by star formation and then cored by feedback processes. Indeed, we report an unprecedentedly observed effect in the DM density profile consisting of a central core combined with an adiabatic contraction at larger galactic radii. The cores obtained are typically $\sim$ 5 kpc large. Moreover, this also impacts the DM density at the solar radius. In our simulation resembling most to the MW, the density is raised from 0.23 GeV/cm$^3$ in the dark matter only run to 0.36 GeV/cm$^3$ (spherical shell) or 0.54 GeV/cm$^3$ (circular ring) in the hydrodynamical run. Studying the subhalos, the dark matter within luminous satellites is also affected by baryonic processes and exihibits cored profiles whereas dark satellites are cuspy. We find a shift in mass compared to DM-only simulations and obtain, for halos in the lower MW mass range, a distribution of luminous satellites comparable to the MW spheroidal dwarf galaxies.

preprint2014arXiv

Baryonic effects on weak-lensing two-point statistics and its cosmological implications

We develop an extension of \textit{the Halo Model} that describes analytically the corrections to the matter power spectrum due to the physics of baryons. We extend these corrections to the weak-lensing shear angular power spectrum. Within each halo, our baryonic model accounts for: 1) a central galaxy, the major stellar component whose properties are derived from abundance matching techniques; 2) a hot plasma in hydrostatic equilibrium and 3) an adiabatically-contracted dark matter component. This analytic approach allows us to compare our model to the dark-matter-only case. Our basic assumptions are tested against the hydrodynamical simulations of Martizzi et. al. (2014), with which a remarkable agreement is found. Our baryonic model has only one free parameter, $M_{\rm crit}$, the critical halo mass that marks the transition between feedback-dominated halos, mostly devoid of gas, and gas rich halos, in which AGN feedback effects become weaker. We explore the entire cosmological parameter space, using the angular power spectrum in three redshift bins as the observable, assuming a Euclid-like survey. We derive the corresponding constraints on the cosmological parameters, as well as the possible bias introduced by neglecting the effects of baryonic physics. We find that, up to $\ell_{max}$=4000, baryonic physics plays very little role in the cosmological parameters estimation. However, if one goes up to $\ell_{max}$=8000, the marginalized errors on the cosmological parameters can be significantly reduced, but neglecting baryonic physics can lead to bias in the recovered cosmological parameters up to 10$σ$. These biases are removed if one takes into account the main baryonic parameter, $M_{\rm crit}$, which can also be determined up to 1-2\%, along with the other cosmological parameters.

preprint2014arXiv

Brightest Cluster Galaxies in Cosmological Simulations with Adaptive Mesh Refinement: Successes and Failures

A large sample of cosmological hydrodynamical zoom-in simulations with Adaptive Mesh Refinement (AMR) is analysed to study the properties of simulated Brightest Cluster Galaxies (BCGs). Following the formation and evolution of BCGs requires modeling an entire galaxy cluster, because the BCG properties are largely influenced by the state of the gas in the cluster and by interactions and mergers with satellites. BCG evolution is also deeply influenced by the presence of gas heating sources such as Active Galactic Nuclei (AGNs) that prevent catastrophic cooling of large amounts of gas. We show that AGN feedback is one of the most important mechanisms in shaping the properties of BCGs at low redshift by analysing our statistical sample of simulations with and without AGN feedback. When AGN feedback is included BCG masses, sizes, star formation rates and kinematic properties are closer to those of the observed systems. Some small discrepancies are observed only for the most massive BCGs and in the fraction of star-forming BCGs, effects that might be due to physical processes that are not included in our model.

preprint2014arXiv

Chameleon $f(R)$ gravity on the Virgo cluster scale

Models of modified gravity offer promising alternatives to the concordance $Λ$CDM cosmology to explain the late-time acceleration of the universe. A popular such model is $f(R)$ gravity, in which the Ricci scalar in the Einstein-Hilbert action is replaced by a general function of it. We study the $f(R)$ model of Hu & Sawicki (2007), which recovers standard General Relativity in high density regimes, while reproducing the desired late-time acceleration at cosmological scales. We run a suite of high resolution zoom simulations using the ECOSMOG code to examine the effect of $f(R)$ gravity on the properties of a halo that is analogous to the Virgo cluster. We show that the velocity dispersion profiles can potentially discriminate between $f(R)$ models and $Λ$CDM, and provide complementary analysis of lensing signal profiles to explore the possibility to further distinguish the different $f(R)$ models. Our results confirm the techniques explored by Cabre et al. (2012) to quantify the effect of environment in the behavior of $f(R)$ gravity, and we extend them to study halo satellites at various redshifts. We find that the modified gravity effects in our models are most observable at low redshifts, and that effects are generally stronger for satellites far from the center of the main halo. We show that the screening properties of halo satellites trace very well that of dark matter particles, which means that low-resolution simulations in which subhalos are not very well resolved can in principle be used to study satellite properties. We discuss observables, particularly for halo satellites, that can potentially be used to constrain the observational viability of $f(R)$ gravity.

preprint2014arXiv

Dancing in the dark: galactic properties trace spin swings along the cosmic web

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.10^10 M_sun. 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.

preprint2014arXiv

Numerical cosmology on the GPU with Enzo and Ramses

A number of scientific numerical codes can currently exploit GPUs with remarkable performance. In astrophysics, Enzo and Ramses are prime examples of such applications. The two codes have been ported to GPUs adopting different strategies and programming models, Enzo adopting CUDA and Ramses using OpenACC. We describe here the different solutions used for the GPU implementation of both cases. Performance benchmarks will be presented for Ramses. The results of the usage of the more mature GPU version of Enzo, adopted for a scientific project within the CHRONOS programme, will be summarised.

preprint2014arXiv

The biasing of baryons on the cluster mass function and cosmological parameter estimation

We study the effect of baryonic processes on the halo mass function in the galaxy cluster mass range using a catalogue of 153 high resolution cosmological hydrodynamical simulations performed with the AMR code ramses. We use the results of our simulations within a simple analytical model to gauge the effects of baryon physics on the halo mass function. Neglect of AGN feedback leads to a significant boost in the cluster mass function similar to that reported by other authors. However, including AGN feedback not only gives rise to systems that are similar to observed galaxy clusters, but they also reverse the global baryonic effects on the clusters. The resulting mass function is closer to the unmodified dark matter halo mass function but still contains a mass dependent bias at the 5-10% level. These effects bias measurements of the cosmological parameters, such as $σ_8$ and $Ω_m$. For current cluster surveys baryonic effects are within the noise for current survey volumes, but forthcoming and planned large SZ, X-ray and multi-wavelength surveys will be biased at the percent level by these processes. The predictions for the halo mass function including baryonic effects need to be carefully studied with larger and improved simulations. However, simulations of full cosmological boxes with the resolution we achieve and including AGN feedback are still computationally challenging.

preprint2014arXiv

The role of turbulence in star formation laws and thresholds

The Schmidt-Kennicutt relation links the surface densities of gas to the star formation rate in galaxies. The physical origin of this relation, and in particular its break, i.e. the transition between an inefficient regime at low gas surface densities and a main regime at higher densities, remains debated. Here, we study the physical origin of the star formation relations and breaks in several low-redshift galaxies, from dwarf irregulars to massive spirals. We use numerical simulations representative of the Milky Way, the Large and the Small Magellanic Clouds with parsec up to subparsec resolution, and which reproduce the observed star formation relations and the relative variations of the star formation thresholds. We analyze the role of interstellar turbulence, gas cooling, and geometry in drawing these relations, at 100 pc scale. We suggest in particular that the existence of a break in the Schmidt- Kennicutt relation could be linked to the transition from subsonic to supersonic turbulence and is independent of self-shielding effects. This transition being connected to the gas thermal properties and thus to the metallicity, the break is shifted toward high surface densities in metal-poor galaxies, as observed in dwarf galaxies. Our results suggest that together with the collapse of clouds under self-gravity, turbulence (injected at galactic scale) can induce the compression of gas and regulate star formation.

preprint2014arXiv

Towards a more realistic sink particle algorithm for the RAMSES code

We present a new sink particle algorithm developed for the Adaptive Mesh Refinement code RAMSES. Our main addition is the use of a clump finder to identify density peaks and their associated regions (the peak patches). This allows us to unambiguously define a discrete set of dense molecular cores as potential sites for sink particle formation. Furthermore, we develop a new scheme to decide if the gas in which a sink could potentially form, is indeed gravitationally bound and rapidly collapsing. This is achieved using a general integral form of the virial theorem, where we use the curvature in the gravitational potential to correctly account for the background potential. We detail all the necessary steps to follow the evolution of sink particles in turbulent molecular cloud simulations, such as sink production, their trajectory integration, sink merging and finally the gas accretion rate onto an existing sink. We compare our new recipe for sink formation to other popular implementations. Statistical properties such as the sink mass function, the average sink mass and the sink multiplicity function are used to evaluate the impact that our new scheme has on accurately predicting fundamental quantities such as the stellar initial mass function or the stellar multiplicity function.

preprint2013arXiv

Beyond the nuclear starburst? Clustered star formation in major mergers

Recent simulation work has successfully captured the formation of the star clusters that have been observed in merging galaxies. These studies, however, tend to focus on studying extreme starbursts, such as the Antennae galaxies. We aim to establish whether there is something special occurring in these extreme systems or whether the mechanism for cluster formation is present in all mergers to a greater or lesser degree. We undertake a general study of merger-induced star formation in a sample of 5 pc resolution adaptive mesh refinement simulations of low redshift equal-mass mergers with randomly-chosen orbital parameters. We find that there is an enhanced mass fraction of very dense gas that appears as the gas density probability density function evolves during the merger. This finding has implications for the interpretation of some observations; a larger mass fraction of dense gas could account for the enhanced HCN/CO ratios seen in ULIRGs and predicts that alpha_CO is lower in mergers, as for a given mass of H_2, CO emission will increase in a denser environment. We also find that as the star formation rate increases, there is a correlated peak in the velocity dispersion of the gas, which we attribute to increasing turbulence driven by the interaction itself. Star formation tends to be clumpy: in some cases there is extended clumpy star formation, but even when star formation is concentrated within the inner kpc (i.e. what may be considered a nuclear starburst) it still often has a clumpy, rather than a smooth, distribution. We find no strong evidence for a clear bimodality in the Kennicutt-Schmidt relation for the average mergers simulated here. Instead, they are typically somewhat offset above the predicted quiescent relation during their starbursts.

preprint2013arXiv

RAMSES-RT: Radiation hydrodynamics in the cosmological context

We present a new implementation of radiation hydrodynamics (RHD) in the adaptive mesh refinement (AMR) code RAMSES. The multi-group 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 large 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.

preprint2013arXiv

Simulating gamma-ray binaries with a relativistic extension of RAMSES

Gamma-ray binaries are composed of a massive star and a rotation-powered pulsar with a highly relativistic wind. The collision between the winds from both objects creates a shock structure where particles are accelerated, resulting in the observed high energy emission. We study the impact of special relativity on the structure and stability of the colliding wind region and highlight the differences with colliding winds from massive stars. We focus on evolution with increasing values of the Lorentz factor of the pulsar wind, keeping in mind that current simulations are unable to reach the expected values of the pulsar wind Lorentz factors by orders of magnitude. We use high resolution numerical simulations with a relativistic extension to the hydrodynamics code RAMSES we have developed. Using 2D simulations, we focus on the region close to the binary, neglecting orbital motion. We use different values of the Lorentz factor of the pulsar wind, up to 16. We find analytic scaling relations between stellar wind collisions and gamma-ray binaries. They provide the position of the contact discontinuity. The position of the shocks strongly depends on the Lorentz factor, the relativistic wind is more collimated than expected based on non-relativistic simulations. Beyond a certain distance, the shocked flow is accelerated to its initial velocity and follows adiabatic expansion. We provide guidance for extrapolation towards more realistic values of the Lorentz factor of the pulsar wind. We extended the adaptive mesh refinement code RAMSES to relativistic hydrodynamics. This code is suited for the study of astrophysical objects such as pulsar wind nebulae, gamma-ray bursts or relativistic jets and will be part of the next public release of RAMSES. Using this code we performed simulations of gamma-ray binaries, highlighting the limits and possibilities of current hydrodynamic models of such systems.

preprint2013arXiv

The long lives of giant clumps and the birth of outflows in gas-rich galaxies at high redshift

Star-forming disk galaxies at high redshift are often subject to violent disk instability, characterized by giant clumps whose fate is yet to be understood. The main question is whether the clumps disrupt within their dynamical timescale (<50Myr), like molecular clouds in today's galaxies, or whether they survive stellar feedback for more than a disk orbital time (~300Myr) in which case they can migrate inward and help building the central bulge. We present 3.5-7pc resolution AMR simulations of high-redshift disks including photo-ionization, radiation pressure, and supernovae feedback (Renaud et al. 2013, and Perret et al., this astro-ph issue). Our modeling of radiation pressure determines the mass loading and initial velocity of winds from basic physical principles. We find that the giant clumps produce steady outflow rates comparable to and sometimes somewhat larger than their star formation rate, with velocities largely sufficient to escape galaxy. The clumps also lose mass, especially old stars, by tidal stripping, and the stellar populations contained in the clumps hence remain relatively young (<=200Myr), as observed. The clumps survive gaseous outflows and stellar loss, because they are wandering in gas-rich turbulent disks from which they can re-accrete gas at high rates compensating for outflows and tidal stripping, overall keeping realistic and self-regulated gaseous and stellar masses. Our simulations produce gaseous outflows with velocities, densities and mass loading consistent with observations, and at the same time suggest that the giant clumps survive for hundreds of Myr and complete their migration to the center of high-redshift galaxies, without rapid dispersion and reformation of clumps.

preprint2013arXiv

Volume Rendering of AMR Simulations

High-resolution simulations often rely on the Adaptive Mesh Resolution (AMR) technique to optimize memory consumption versus attainable precision. While this technique allows for dramatic improvements in terms of computing performance, the analysis and visualization of its data outputs remain challenging. The lack of effective volume renderers for the octree-based AMR used by the RAMSES simulation program has led to the development of the solutions presented in this paper. Two custom algorithms are discussed, based on the splatting and the ray-casting techniques. Their usage is illustrated in the context of the visualization of a high-resolution, 6000-processor simulation of a Milky Way-like galaxy. Performance obtained in terms of memory management and parallelism speedup are presented.

preprint2012arXiv

A diversity of progenitors and histories for isolated spiral galaxies

We analyze a suite of 33 cosmological simulations of the evolution of Milky Way-mass galaxies in low-density environments. Our sample spans a broad range of Hubble types at z=0, from nearly bulgeless disks to bulge-dominated galaxies. Despite the fact that a large fraction of the bulge is typically in place by z=1, we find no significant correlation between the morphology at z=1 and at z=0. The z=1 progenitors of disk galaxies span a range of morphologies, including smooth disks, unstable disks, interacting galaxies and bulge-dominated systems. By z=0.5, spiral arms and bars are largely in place and the progenitor morphology is correlated with the final morphology. We next focus on late-type galaxies with a bulge-to-total ratio B/T<0.3 at z=0. These show a correlation between B/T at z=0 and the mass ratio of the largest merger at z<2, as well as with the gas accretion rate at z>1. We find that the galaxies with the lowest B/T tend to have a quiet baryon input history, with no major mergers at z<2, and with a low and constant gas accretion rate that keeps a stable angular-momentum direction. More violent merger or gas accretion histories lead to galaxies with more prominent bulges. Most disk galaxies have a bulge Sersic index n<2. The galaxies with the highest bulge Sersic index tend to have histories of intense gas accretion and disk instability rather than active mergers.

preprint2012arXiv

Co-planar streams, pancakes, and angular-momentum exchange in high-z disc galaxies

We study the feeding of massive galaxies at high redshift through streams from the cosmic web using the Mare Nostrum hydro-cosmological simulation. Our statistical sample consists of 350 dark-matter haloes of ~10^12 Msun at z = 2.5. We find that ~70% of the influx into the virial radius Rv is in narrow streams covering 10% of the virial shell. On average 64% of the stream influx is in one stream, and 95% is in three dominant streams. The streams that feed a massive halo tend to lie in a plane that extends from half to a few Rv, hereafter "the stream plane" (SP). The streams are typically embedded in a thin sheet of low-entropy gas, a Zel'dovich pancake, which carries ~20% of the influx into Rv. The filaments-in-a-plane configuration about the massive haloes at the nodes of the cosmic web differs from the large- scale structure of the web where the filaments mark the intersections of slanted sheets. The stream plane is only weakly aligned with the angular momentum (AM) near Rv, consistent with the fact that typically 80% of the AM is carried by one dominant stream. The galactic disc plane shows a weak tendency to be perpendicular to the large-scale SP, consistent with tidal-torque theory. Most interesting, the direction of the disc AM is only weakly correlated with the AM direction at Rv. This indicates a significant AM exchange at the interphase between streams and disc in the greater environment of the disc inside an "AM sphere" of radius ~0.3Rv . The required large torques are expected based on the perturbed morphology and kinematics within this interaction sphere. This AM exchange may or may not require a major modification of the standard disc modeling based on AM conservation, depending on the extent to which the amplitude of the disc AM is affected, which is yet to be studied.

preprint2012arXiv

Cusp-core transformations in dwarf galaxies: observational predictions

The presence of a dark matter core in the central kiloparsec of many dwarf galaxies has been a long standing problem in galaxy formation theories based on the standard cold dark matter paradigm. Recent cosmological simulations, based on Smooth Particle Hydrodynamics and rather strong feedback recipes have shown that it was indeed possible to form extended dark matter cores using baryonic processes related to a more realistic treatment of the interstellar medium. Using adaptive mesh refinement, together with a new, stronger supernovae feedback scheme that we have recently implemented in the RAMSES code, we show that it is also possible to form a prominent dark matter core within the well-controlled framework of an isolated, initially cuspy, 10 billion solar masses dark matter halo. Although our numerical experiment is idealized, it allows a clean and unambiguous identification of the dark matter core formation process. Our dark matter inner profile is well fitted by a pseudo-isothermal profile with a core radius of 800 pc. The core formation mechanism is consistent with the one proposed recently by Pontzen & Governato. We highlight two key observational predictions of all simulations that find cusp-core transformations: (i) a bursty star formation history with peak to trough ratio of 5 to 10 and a duty cycle comparable to the local dynamical time; and (ii) a stellar distribution that is hot with v/sigma=1. We compare the observational properties of our model galaxy with recent measurements of the isolated dwarf WLM. We show that the spatial and kinematical distribution of stars and HI gas are in striking agreement with observations, supporting the fundamental role played by stellar feedback in shaping both the stellar and dark matter distribution.

preprint2012arXiv

Hydrodynamics of galaxy mergers with supermassive black holes: is there a last parsec problem ?

We study the formation of a supermassive black hole (SMBH) binary and the shrinking of the separation of the two holes to sub-pc scales starting from a realistic major merger between two gas-rich spiral galaxies with mass comparable to our Milky Way. The simulations, carried out with the Adaptive Mesh Refinement (AMR) code RAMSES, are capable of resolving separations as small as 0.1 pc. The collision of the two galaxies produces a gravo-turbulent rotating nuclear disk with mass (10^9 Msun) and size (60 pc) in excellent agreement with previous SPH simulations with particle splitting that used a similar setup (Mayer et al. 2007) but were limited to separations of a few parsecs. The AMR results confirm that the two black holes sink rapidly as a result of dynamical friction onto the gaseous background, reaching a separation of 1 pc in less than 10^7 yr. We show that the dynamical friction wake is well resolved by our model and we find good agreement with analytical predictions of the drag force as a function of the Mach number. Below 1 pc, black hole pairing slows down significantly, as the relative velocity between the sinking SMBH becomes highly subsonic and the mass contained within their orbit falls below the mass of the binary itself, rendering dynamical friction ineffective. In this final stage, the black holes have not opened a gap as the gaseous background is highly pressurized in the center. Non-axisymmetric gas torques do not arise to restart sinking in absence of efficient dynamical friction, at variance with previous calculations using idealized equilibrium nuclear disk models. (abridged)

preprint2012arXiv

Indirect dark matter searches: towards a consistent top-bottom approach for studying the gamma-ray signals and associated backgrounds

While dark matter (DM) is the key ingredient for a successful theory of structure formation, its microscopic nature remains elusive. Indirect detection may provide a powerful test for some strongly motivated DM particle models. Nevertheless, astrophysical backgrounds are usually expected with amplitudes and spectral features similar to the chased signals. On galactic scales, these backgrounds arise from interactions of cosmic rays (CRs) with the interstellar gas, both being difficult to infer and model in detail from observations. Moreover, the associated predictions unavoidably come with theoretical errors, which are known to be significant. We show that a trustworthy guide for such challenging searches can be obtained by exploiting the full information contained in cosmological simulations of galaxies, which now include baryonic gas dynamics and star formation. We further insert CR production and transport from the identified supernova events and fully calculate the CR distribution in a simulated galaxy. We focus on diffuse gamma-rays, and self-consistently calculate both the astrophysical galactic emission and the dark matter signal. We notably show that adiabatic contraction does not necessarily induce large signal-to-noise ratios in galactic centers, and could anyway be traced from the astrophysical background itself. We finally discuss how all this may be used as a generic diagnostic tool for galaxy formation.

preprint2012arXiv

Radiatively Efficient Magnetized Bondi Accretion

We have carried out a numerical study of the effect of large scale magnetic fields on the rate of accretion from a uniform, isothermal gas onto a resistive, stationary point mass. Only mass, not magnetic flux, accretes onto the point mass. The simulations for this study avoid complications arising from boundary conditions by keeping the boundaries far from the accreting object. Our simulations leverage adaptive refinement methodology to attain high spatial fidelity close to the accreting object. Our results are particularly relevant to the problem of star formation from a magnetized molecular cloud in which thermal energy is radiated away on time scales much shorter than the dynamical time scale. Contrary to the adiabatic case, our simulations show convergence toward a finite accretion rate in the limit in which the radius of the accreting object vanishes, regardless of magnetic field strength. For very weak magnetic fields, the accretion rate first approaches the Bondi value and then drops by a factor ~ 2 as magnetic flux builds up near the point mass. For strong magnetic fields, the steady-state accretion rate is reduced by a factor ~ 0.2 β^{1/2} compared to the Bondi value, where βis the ratio of the gas pressure to the magnetic pressure. We give a simple expression for the accretion rate as a function of the magnetic field strength. Approximate analytic results are given in the Appendixes for both time-dependent accretion in the limit of weak magnetic fields and steady-state accretion for the case of strong magnetic fields.

preprint2012arXiv

RAMSES-CH: A New Chemodynamical Code for Cosmological Simulations

We present a new chemodynamical code - Ramses-CH - for use in simulating the self-consistent evolution of chemical and hydrodynamical properties of galaxies within a fully cosmological framework. We build upon the adaptive mesh refinement code Ramses, which includes a treatment of self-gravity, hydrodynamics, star formation, radiative cooling, and supernovae feedback, to trace the dominant isotopes of C, N, O, Ne, Mg, Si, and Fe. We include the contribution of Type Ia and II supernovae, in addition to low- and intermediate-mass asymptotic giant branch stars, relaxing the instantaneous recycling approximation. The new chemical evolution modules are highly flexible and portable, lending themselves to ready exploration of variations in the underpining stellar and nuclear physics. We apply Ramses-CH to the cosmological simulation of a typical L* galaxy, demonstrating the successful recovery of the basic empirical constraints regarding, [α/Fe]-[Fe/H] and Type Ia/II supernovae rates.

preprint2012arXiv

The effects of baryon physics, black holes and AGN feedback on the mass distribution in clusters of galaxies

The spatial distribution of matter in clusters of galaxies is mainly determined by the dominant dark matter component, however, physical processes involving baryonic matter are able to modify it significantly. We analyse a set of 500 pc resolution cosmological simulations of a cluster of galaxies with mass comparable to Virgo, performed with the AMR code RAMSES. We compare the mass density profiles of the dark, stellar and gaseous matter components of the cluster that result from different assumptions for the subgrid baryonic physics and galaxy formation processes. First, the prediction of a gravity only N-body simulation is compared to that of a hydrodynamical simulation with standard galaxy formation recipes, then all results are compared to a hydrodynamical simulation which includes thermal AGN feedback from Super Massive Black Holes (SMBH). We find the usual effects of overcooling and adiabatic contraction in the run with standard galaxy formation physics, but very different results are found when implementing SMBHs and AGN feedback. Star formation is strongly quenched, producing lower stellar densities throughout the cluster, and much less cold gas is available for star formation at low redshifts. At redshift z = 0 we find a flat density core of radius 10 kpc in both of the dark and stellar matter density profiles. We specu- late on the possible formation mechanisms able to produce such cores and we conclude that they can be produced through the coupling of different processes: (I) dynamical friction from the decay of black hole orbits during galaxy mergers; (II) AGN driven gas outflows producing fluctuations of the gravitational potential causing the removal of collisionless matter from the central region of the cluster; (III) adiabatic expansion in response to the slow expulsion of gas from the central region of the cluster during the quiescent mode of AGN activity.

preprint2012arXiv

The formation of the brightest cluster galaxies in cosmological simulations: the case for AGN feedback

We use 500 pc resolution cosmological simulations of a Virgo-like galaxy cluster to study the properties of the brightest cluster galaxy (BCG) that forms at the center of the halo. We compared two simulations; one incorporating only supernovae feedback and a second that also includes prescriptions for black hole growth and the resulting AGN feedback from gas accretion. As previous work has shown, with supernovae feedback alone we are unable to reproduce any of the observed properties of massive cluster ellipticals. The resulting BCG is rotating quickly, has a high Sersic index, a strong mass excess in the center and a total central density profile falling more steeply than isothermal. Furthermore, it is far too efficient at converting most of the available baryons into stars which is strongly constrained by abundance matching. With a treatment of black hole dynamics and AGN feedback the BCG properties are in good agreement with data: they rotate slowly, have a cored surface density profile, a flat or rising velocity dispersion profile and a low stellar mass fraction. The AGN provides a new mechanism to create cores in luminous elliptical galaxies; the core expands due to the combined effects of heating from dynamical friction of sinking massive black holes and AGN feedback that ejects gaseous material from the central regions.

preprint2011arXiv

A simple multigrid scheme for solving the Poisson equation with arbitrary domain boundaries

We present a new multigrid scheme for solving the Poisson equation with Dirichlet boundary conditions on a Cartesian grid with irregular domain boundaries. This scheme was developed in the context of the Adaptive Mesh Refinement (AMR) schemes based on a graded-octree data structure. The Poisson equation is solved on a level-by-level basis, using a "one-way interface" scheme in which boundary conditions are interpolated from the previous coarser level solution. Such a scheme is particularly well suited for self-gravitating astrophysical flows requiring an adaptive time stepping strategy. By constructing a multigrid hierarchy covering the active cells of each AMR level, we have designed a memory-efficient algorithm that can benefit fully from the multigrid acceleration. We present a simple method for capturing the boundary conditions across the multigrid hierarchy, based on a second-order accurate reconstruction of the boundaries of the multigrid levels. In case of very complex boundaries, small scale features become smaller than the discretization cell size of coarse multigrid levels and convergence problems arise. We propose a simple solution to address these issues. Using our scheme, the convergence rate usually depends on the grid size for complex grids, but good linear convergence is maintained. The proposed method was successfully implemented on distributed memory architectures in the RAMSES code, for which we present and discuss convergence and accuracy properties as well as timing performances.

preprint2011arXiv

AGN feedback using AMR cosmological simulations

Feedback processes are thought to solve some of the long-standing issues of the numerical modelling of galaxy formation: over-cooling, low angular momentum, massive blue galaxies, extra-galactic enrichment, etc. The accretion of gas onto super-massive black holes in the centre of massive galaxies can release tremendous amounts of energy to the surrounding medium. We show, with cosmological Adaptive Mesh Refinement simulations, how the growth of black holes is regulated by the feedback from Active Galactic Nuclei using a new dual jet/heating mechanism. We discuss how this large amount of feedback is able to modify the cold baryon content of galaxies, and perturb the properties of the hot plasma in their vicinity.

preprint2011arXiv

Black Hole growth and AGN obscuration by instability-driven inflows in high-redshift disk galaxies fed by cold streams

Disk galaxies at high redshift have been predicted to maintain high gas surface densities due to continuous feeding by intense cold streams leading to violent gravitational instability, transient features and giant clumps. Gravitational torques between the perturbations drive angular momentum out and mass in, and the inflow provides the energy for keeping strong turbulence. We use analytic estimates of the inflow for a self-regulated unstable disk at a Toomre stability parameter Q~1, and isolated galaxy simulations capable of resolving the nuclear inflow down to the central parsec. We predict an average inflow rate ~10 Msun/yr through the disk of a 10^11 Msun galaxy, with conditions representative of z~2 stream-fed disks. The inflow rate scales with disk mass and (1+z)^{3/2}. It includes clump migration and inflow of the smoother component, valid even if clumps disrupt. This inflow grows the bulge, while only a fraction ~ 10^-3 of it needs to accrete onto a central black hole (BH), in order to obey the observed BH-bulge relation. A galaxy of 10^11 Msun at z~2 is expected to host a BH of ~10^8 Msun, accreting on average with moderate sub-Eddington luminosity L_X ~ 10^42-43 erg/s, accompanied by brighter episodes when dense clumps coalesce. We note that in rare massive galaxies at z~6, the same process may feed 10^9 Msun BH at the Eddington rate. High central gas column densities can severely obscure AGN in high-redshift disks, possibly hindering their detection in deep X-ray surveys.

preprint2011arXiv

Comparing Numerical Methods for Isothermal Magnetized Supersonic Turbulence

We employ simulations of supersonic super-Alfvenic turbulence decay as a benchmark test problem to assess and compare the performance of nine astrophysical MHD methods actively used to model star formation. The set of nine codes includes: ENZO, FLASH, KT-MHD, LL-MHD, PLUTO, PPML, RAMSES, STAGGER, and ZEUS. We present a comprehensive set of statistical measures designed to quantify the effects of numerical dissipation in these MHD solvers. We compare power spectra for basic fields to determine the effective spectral bandwidth of the methods and rank them based on their relative effective Reynolds numbers. We also compare numerical dissipation for solenoidal and dilatational velocity components to check for possible impacts of the numerics on small-scale density statistics. Finally, we discuss convergence of various characteristics for the turbulence decay test and impacts of various components of numerical schemes on the accuracy of solutions. We show that the best performing codes employ a consistently high order of accuracy for spatial reconstruction of the evolved fields, transverse gradient interpolation, conservation law update step, and Lorentz force computation. The best results are achieved with divergence-free evolution of the magnetic field using the constrained transport method, and using little to no explicit artificial viscosity. Codes which fall short in one or more of these areas are still useful, but they must compensate higher numerical dissipation with higher numerical resolution. This paper is the largest, most comprehensive MHD code comparison on an application-like test problem to date. We hope this work will help developers improve their numerical algorithms while helping users to make informed choices in picking optimal applications for their specific astrophysical problems.

preprint2011arXiv

Disc Heating: Comparing the Milky Way with Cosmological Simulations

We present the analysis of a suite of simulations run with different particle-and grid-based cosmological hydrodynamical codes and compare them with observational data of the Milky Way. This is the first study to make comparisons of properties of galaxies simulated with particle and grid-based codes. Our analysis indicates that there is broad agreement between these different modelling techniques. We study the velocity dispersion - age relation for disc stars at z=0 and find that four of the simulations are more consistent with observations by Holmberg et al. (2008) in which the stellar disc appears to undergo continual/secular heating. Two other simulations are in better agreement with the Quillen & Garnett (2001) observations that suggest a "saturation" in the heating profile for young stars in the disc. None of the simulations have thin discs as old as that of the Milky Way. We also analyse the kinematics of disc stars at the time of their birth for different epochs in the galaxies' evolution and find that in some simulations old stars are born cold within the disc and are subsequently heated, while other simulations possess old stellar populations which are born relatively hot. The models which are in better agreement with observations of the Milky Way's stellar disc undergo significantly lower minor-merger/assembly activity after the last major merger - i.e. once the disc has formed. All of the simulations are significantly "hotter" than the Milky Way disc; on top of the effects of mergers, we find a "floor" in the dispersion that is related to the underlying treatment of the heating and cooling of the interstellar medium, and the low density threshold which such codes use for star formation. This finding has important implications for all studies of disc heating that use hydrodynamical codes.

preprint2011arXiv

ECOSMOG: An Efficient Code for Simulating Modified Gravity

We introduce a new code, ECOSMOG, to run N-body simulations for a wide class of modified gravity and dynamical dark energy theories. These theories generally have one or more new dynamical degrees of freedom, the dynamics of which are governed by their (usually rather nonlinear) equations of motion. Solving these non-linear equations has been a great challenge in cosmology. Our code is based on the RAMSES code, which solves the Poisson equation on adaptively refined meshes to gain high resolutions in the high-density regions. We have added a solver for the extra degree(s) of freedom and performed numerous tests for the f(R) gravity model as an example to show its reliability. We find that much higher efficiency could be achieved compared with other existing mesh/grid-based codes thanks to two new features of the present code: (1) the efficient parallelisation and (2) the usage of the multigrid relaxation to solve the extra equation(s) on both the regular domain grid and refinements, giving much faster convergence even under much more stringent convergence criteria. This code is designed for performing high-accuracy, high-resolution and large-volume cosmological simulations for modified gravity and general dark energy theories, which can be utilised to test gravity and the dark energy hypothesis using the upcoming and future deep and high-resolution galaxy surveys.

preprint2011arXiv

Grid of Lya radiation transfer models for the interpretation of distant galaxies

Lya is a key diagnostic for numerous observations of distant star-forming galaxies. It's interpretation requires, however, detailed radiation transfer models. We provide an extensive grid of 3D radiation transfer models simulating the Lya and UV continuum radiation transfer in the interstellar medium of star-forming galaxies. We have improved our Monte Carlo MCLya code, and have used it to compute a grid of 6240 radiation transfer models for homogeneous spherical shells containing HI and dust surrounding a central source. The simulations cover a wide range of parameter space. We present the detailed predictions from our models including in particular the Lya escape fraction fesc, the continuum attenuation, and detailed Lya line profiles. The Lya escape fraction is shown to depend strongly on dust content, but also on other parameters (HI column density and radial velocity). The predicted line profiles show a great diversity of morphologies ranging from broad absorption lines to emission lines with complex features. The results from our simulations are distributed in electronic format. Our models should be of use for the interpretation of observations from distant galaxies, for other simulations, and should also serve as an important base for comparison for future, more refined, radiation transfer models.

preprint2011arXiv

How AGN feedback and metal cooling shape cluster entropy profiles

Observed clusters of galaxies essentially come in two flavors: non cool core clusters characterized by an isothermal temperature profile and a central entropy floor, and cool-core clusters where temperature and entropy in the central region are increasing with radius. Using cosmological resimulations of a galaxy cluster, we study the evolution of its intracluster medium (ICM) gas properties, and through them we assess the effect of different (sub-grid) modelling of the physical processes at play, namely gas cooling, star formation, feedback from supernovae and active galactic nuclei (AGN). More specifically we show that AGN feedback plays a major role in the pre-heating of the proto-cluster as it prevents a high concentration of mass from collecting in the center of the future galaxy cluster at early times. However, AGN activity during the cluster's later evolution is also required to regulate the mass flow into its core and prevent runaway star formation in the central galaxy. Whereas the energy deposited by supernovae alone is insufficient to prevent an overcooling catastrophe, supernovae are responsible for spreading a large amount of metals at high redshift, enhancing the cooling efficiency of the ICM gas. As the AGN energy release depends on the accretion rate of gas onto its central black hole engine, the AGN responds to this supernova enhanced gas accretion by injecting more energy into the surrounding gas, and as a result increases the amount of early pre-heating. We demonstrate that the interaction between an AGN jet and the ICM gas that regulates the growth of the AGN's BH, can naturally produce cool core clusters if we neglect metals. However, as soon as metals are allowed to contribute to the radiative cooling, only the non cool core solution is produced.

preprint2011arXiv

Hydrodynamics of high-redshift galaxy collisions: From gas-rich disks to dispersion-dominated mergers and compact spheroids

Disk galaxies at high redshift (z~2) are characterized by high fractions of cold gas, strong turbulence, and giant star-forming clumps. Major mergers of disk galaxies at high redshift should then generally involve such turbulent clumpy disks. Merger simulations, however, model the ISM as a stable, homogeneous, and thermally pressurized medium. We present the first merger simulations with high fractions of cold, turbulent, and clumpy gas. We discuss the major new features of these models compared to models where the gas is artificially stabilized and warmed. Gas turbulence, which is already strong in high-redshift disks, is further enhanced in mergers. Some phases are dispersion-dominated, with most of the gas kinetic energy in the form of velocity dispersion and very chaotic velocity fields, unlike merger models using a thermally stabilized gas. These mergers can reach very high star formation rates, and have multi-component gas spectra consistent with SubMillimeter Galaxies. Major mergers with high fractions of cold turbulent gas are also characterized by highly dissipative gas collapse to the center of mass, with the stellar component following in a global contraction. The final galaxies are early-type with relatively small radii and high Sersic indices, like high-redshift compact spheroids. The mass fraction in a disk component that survives or re-forms after a merger is severely reduced compared to models with stabilized gas, and the formation of a massive disk component would require significant accretion of external baryons afterwards. Mergers thus appear to destroy extended disks even when the gas fraction is high, and this lends further support to smooth infall as the main formation mechanism for massive disk galaxies.

preprint2011arXiv

Radiation hydrodynamics with Adaptive Mesh Refinement and application to prestellar core collapse. I Methods

Radiative transfer has a strong impact on the collapse and the fragmentation of prestellar dense cores. We present the radiation-hydrodynamics solver we designed for the RAMSES code. The method is designed for astrophysical purposes, and in particular for protostellar collapse. We present the solver, using the co-moving frame to evaluate the radiative quantities. We use the popular flux limited diffusion approximation, under the grey approximation (one group of photon). The solver is based on the second-order Godunov scheme of RAMSES for its hyperbolic part, and on an implicit scheme for the radiation diffusion and the coupling between radiation and matter. We report in details our methodology to integrate the RHD solver into RAMSES. We test successfully the method against several conventional tests. For validation in 3D, we perform calculations of the collapse of an isolated 1 M_sun prestellar dense core, without rotation. We compare successfully the results with previous studies using different models for radiation and hydrodynamics. We have developed a full radiation hydrodynamics solver in the RAMSES code, that handles adaptive mesh refinement grids. The method is a combination of an explicit scheme and an implicit scheme, accurate to the second-order in space. Our method is well suited for star formation purposes. Results of multidimensional dense core collapse calculations with rotation are presented in a companion paper.

preprint2011arXiv

Self-regulated growth of supermassive black holes by a dual jet/heating AGN feedback mechanism: methods, tests and implications for cosmological simulations

We develop a new sub-grid model for the growth of supermassive Black Holes (BHs) and their associated Active Galactic Nuclei (AGN) feedback in hydrodynamical cosmological simulations. Assuming that BHs are created in the early stages of galaxy formation, they grow by mergers and accretion of gas at a Eddington-limited Bondi accretion rate. However this growth is regulated by AGN feedback which we model using two different modes: a quasar-heating mode when accretion rates onto the BHs are comparable to the Eddington rate, and a radio-jet mode at lower accretion rates. In other words, our feedback model deposits energy as a succession of thermal bursts and jet outflows depending on the properties of the gas surrounding the BHs. We assess the plausibility of such a model by comparing our results to observational measurements of the coevolution of BHs and their host galaxy properties, and check their robustness with respect to numerical resolution. We show that AGN feedback must be a crucial physical ingredient for the formation of massive galaxies as it appears to be the only physical mechanism able to efficiently prevent the accumulation of and/or expel cold gas out of halos/galaxies and significantly suppress star formation. Our model predicts that the relationship between BHs and their host galaxy mass evolves as a function of redshift, because of the vigorous accretion of cold material in the early Universe that drives Eddington-limited accretion onto BHs. Quasar activity is also enhanced at high redshift. However, as structures grow in mass and lose their cold material through star formation and efficient BH feedback ejection, the AGN activity in the low-redshift Universe becomes more and more dominated by the radio mode, which powers jets through the hot circum-galactic medium.

preprint2011arXiv

The impact of ISM turbulence, clustered star formation and feedback on galaxy mass assembly through cold flows and mergers

Two of the dominant channels for galaxy mass assembly are cold flows (cold gas supplied via the filaments of the cosmic web) and mergers. How these processes combine in a cosmological setting, at both low and high redshift, to produce the whole zoo of galaxies we observe is largely unknown. Indeed there is still much to understand about the detailed physics of each process in isolation. While these formation channels have been studied using hydrodynamical simulations, here we study their impact on gas properties and star formation (SF) with some of the first simulations that capture the multiphase, cloudy nature of the interstellar medium (ISM), by virtue of their high spatial resolution (and corresponding low temperature threshold). In this regime, we examine the competition between cold flows and a supernovae(SNe)-driven outflow in a very high-redshift galaxy (z {\approx} 9) and study the evolution of equal-mass galaxy mergers at low and high redshift, focusing on the induced SF. We find that SNe-driven outflows cannot reduce the cold accretion at z {\approx} 9 and that SF is actually enhanced due to the ensuing metal enrichment. We demonstrate how several recent observational results on galaxy populations (e.g. enhanced HCN/CO ratios in ULIRGs, a separate Kennicutt Schmidt (KS) sequence for starbursts and the population of compact early type galaxies (ETGs) at high redshift) can be explained with mechanisms captured in galaxy merger simulations, provided that the multiphase nature of the ISM is resolved.

preprint2010arXiv

Cell Count Moments in the Halo Model

We study cell count moments up to fifth order of the distributions of haloes, of halo substructures as a proxy for galaxies, and of mass in the context of the halo model and compare theoretical predictions to the results of numerical simulations. On scales larger than the size of the largest cluster, we present a simple point cluster model in which results depend only on cluster-cluster correlations and on the distribution of the number of objects within a cluster, or cluster occupancy. The point cluster model leads to expressions for moments of galaxy counts in which the volume-averaged moments on large scales approach those of the halo distribution and on smaller scales exhibit hierarchical clustering with amplitudes $S_k$ determined by moments of the occupancy distribution. In this limit, the halo model predictions are purely combinatoric, and have no dependence on halo profile, concentration parameter, or potential asphericity. The full halo model introduces only two additional effects: on large scales, haloes of different mass have different clustering strengths, introducing relative bias parameters; and on the smallest scales, halo structure is resolved and details of the halo profile become important, introducing shape-dependent form factors. Because of differences between discrete and continuous statistics, the hierarchical amplitudes for galaxies and for mass behave differently on small scales even if galaxy number is exactly proportional to mass, a difference that is not necessarily well described in terms of bias.

preprint2010arXiv

Galactic star formation in parsec-scale resolution simulations

The interstellar medium (ISM) in galaxies is multiphase and cloudy, with stars forming in the very dense, cold gas found in Giant Molecular Clouds (GMCs). Simulating the evolution of an entire galaxy, however, is a computational problem which covers many orders of magnitude, so many simulations cannot reach densities high enough or temperatures low enough to resolve this multiphase nature. Therefore, the formation of GMCs is not captured and the resulting gas distribution is smooth, contrary to observations. We investigate how star formation (SF) proceeds in simulated galaxies when we obtain parsec-scale resolution and more successfully capture the multiphase ISM. Both major mergers and the accretion of cold gas via filaments are dominant contributors to a galaxy's total stellar budget and we examine SF at high resolution in both of these contexts.

preprint2010arXiv

Jet-regulated cooling catastrophe

We present the first implementation of Active Galactic Nuclei (AGN) feedback in the form of momentum driven jets in an Adaptive Mesh Refinement (AMR) cosmological resimulation of a galaxy cluster. The jets are powered by gas accretion onto Super Massive Black Holes (SMBHs) which also grow by mergers. Throughout its formation, the cluster experiences different dynamical states: both a morphologically perturbed epoch at early times and a relaxed state at late times allowing us to study the different modes of BH growth and associated AGN jet feedback. BHs accrete gas efficiently at high redshift (z>2), significantly pre-heating proto-cluster halos. Gas-rich mergers at high redshift also fuel strong, episodic jet activity, which transports gas from the proto-cluster core to its outer regions. At later times, while the cluster relaxes, the supply of cold gas onto the BHs is reduced leading to lower jet activity. Although the cluster is still heated by this activity as sound waves propagate from the core to the virial radius, the jets inefficiently redistribute gas outwards and a small cooling flow develops, along with low-pressure cavities similar to those detected in X-ray observations. Overall, our jet implementation of AGN feedback quenches star formation quite efficiently, reducing the stellar content of the central cluster galaxy by a factor 3 compared to the no AGN case. It also dramatically alters the shape of the gas density profile, bringing it in close agreement with the beta model favoured by observations, producing quite an isothermal galaxy cluster for gigayears in the process. However, it still falls short in matching the lower than Universal baryon fractions which seem to be commonplace in observed galaxy clusters.

preprint2010arXiv

Magnetised winds in dwarf galaxies

The origin and the magnitude of the inter-galactic magnetic field is of primordial importance in the global picture of magnetic field evolution, as it is considered to be the missing link between galactic magnetic fields and cluster magnetic fields on much larger scales. We are testing whether dwarf galaxies are good candidates to explain the enrichment of the IGM: after their discs form and trigger galactic dynamos, supernova feedback will launch strong winds, expelling magnetic field lines in the IGM. We have performed MHD simulations of an isolated dwarf galaxy, forming self-consistently inside a cooling halo. Using the RAMSES code, we have for the first time simulated the formation of a magnetised supernova-driven galactic outflow. This simulation is an important step towards a more realistic modelling using fully cosmological simulations. Our simulations reproduce well the observed properties of magnetic fields in spiral galaxies. The formation and the evolution of our simulated disc lead to a strong magnetic field amplification: the magnetic field in the final wind bubble is one order of magnitude larger than the initial value. The magnetic field in the disc, essentially toroidal, is growing linearly with time as a consequence of differential rotation. We discuss the consequence of this simple mechanism on the cosmic evolution of the magnetic field: we propose a new scenario for the evolution of the magnetic field, with dwarf galaxies playing a key role in amplifying and ejecting magnetic energy in the IGM, resulting in what we call a "Cosmic Dynamo" that could contribute to the rather high field strengths observed in galaxies and clusters today.

preprint2010arXiv

Mass Distribution in Galaxy Clusters: the Role of AGN Feedback

We use 1 kpc resolution cosmological AMR simulations of a Virgo-like galaxy cluster to investigate the effect of feedback from supermassive black holes (SMBH) on the mass distribution of dark matter, gas and stars. We compared three different models: (i) a standard galaxy formation model featuring gas cooling, star formation and supernovae feedback, (ii) a "quenching" model for which star formation is artificially suppressed in massive halos and finally (iii) the recently proposed AGN feedback model of Booth & Schaye (2009). Without AGN feedback (even in the quenching case), our simulated cluster suffers from a strong overcooling problem, with a stellar mass fraction significantly above observed values in M87. The baryon distribution is highly concentrated, resulting in a strong adiabatic contraction (AC) of dark matter. With AGN feedback, on the contrary, the stellar mass in the bright central galaxy (BCG) lies below observational estimates and the overcooling problem disappears. The stellar mass of the BCG is seen to increase with increasing mass resolution, suggesting that our stellar masses converges to the correct value from below. The gas and total mass distributions are in striking agreement with observations. We also find a slight deficit (~10%) of baryons at the virial radius, due to the effect of AGN-driven shock waves pushing gas to Mpc scales and beyond. This baryon deficit results in a slight adiabatic expansion of the dark matter distribution, that can be explained quantitatively by AC theory.

preprint2010arXiv

Radiative, magnetic and numerical feedbacks on small-scale fragmentation

Radiative feedback and magnetic field are understood to have a strong impact on the protostellar collapse. We present high resolution numerical calculations of the collapse of a 1 solar mass dense core in solid body rotation, including both radiative transfer and magnetic field. Using typical parameters for low-mass cores, we study thoroughly the effect of radiative transfer and magnetic field on the first core formation and fragmentation. We show that including the two aforementioned physical processes does not correspond to the simple picture of adding them separately. The interplay between the two is extremely strong, via the magnetic braking and the radiation from the accretion shock.

preprint2010arXiv

Reionization simulations powered by GPUs I: the structure of the Ultraviolet radiation field

We present a set of cosmological simulations with radiative transfer in order to model the reionization history of the Universe. Galaxy formation and the associated star formation are followed self-consistently with gas and dark matter dynamics using the RAMSES code, while radiative transfer is performed as a post-processing step using a moment-based method with M1 closure relation in the ATON code. The latter has been ported to a multiple Graphics Processing Units (GPU) architecture using CUDA + MPI, resulting in an overall acceleration (x80) that allows us to tackle radiative transfer problems at resolution of 1024^3 + 2 levels of refinement for the hydro adaptive grid and 1024^3 for the RT cartesian grid. We observe a good convergence between our different resolution runs as long as the effects of finite resolution on the star formation history are properly taken into account. We also show that the neutral fraction depends on the total mass density, in a way close to the predictions of photoionization equilibrium, as long as the effect of self-shielding is included in the background radiation model. However we still fail at reproducing the z=6 constraints on the H neutral fraction and the intensity of the UV background. In order to account for unresolved density fluctuations, we added a simple clumping factor model. Using our most spatially resolved simulation (12.5 Mpc/h-1024^3) to calibrate our subgrid model, we have resimulated our largest box (100 Mpc/h 1024^3), successfully reproducing the observed level of H neutral fraction at z=6. We don't reproduce the photoionization rate inferred from the same observations. We argue that this discrepancy could be explained by the fact that the average radiation intensity and the average neutral fraction depends on different regions of the gas density distribution, so that one quantity cannot be simply deduced from the other.

preprint2010arXiv

Simulation of the growth of the 3D Rayleigh-Taylor instability in Supernova Remnants using an expanding reference frame

Context: The Rayleigh-Taylor instabilities generated by the deceleration of a supernova remnant during the ejecta-dominated phase are known to produce finger-like structures in the matter distribution which modify the geometry of the remnant. The morphology of supernova remnants is also expected to be modified when efficient particle acceleration occurs at their shocks. Aims: The impact of the Rayleigh-Taylor instabilities from the ejecta-dominated to the Sedov-Taylor phase is investigated over one octant of the supernova remnant. We also study the effect of efficient particle acceleration at the forward shock on the growth of the Rayleigh-Taylor instabilities. Methods: We modified the Adaptive Mesh Refinement code RAMSES to study with hydrodynamic numerical simulations the evolution of supernova remnants in the framework of an expanding reference frame. The adiabatic index of a relativistic gas between the forward shock and the contact discontinuity mimics the presence of accelerated particles. Results: The great advantage of the super-comoving coordinate system adopted here is that it minimizes numerical diffusion at the contact discontinuity, since it is stationary with respect to the grid. We propose an accurate expression for the growth of the Rayleigh-Taylor structures that connects smoothly the early growth to the asymptotic self-similar behaviour. Conclusions: The development of the Rayleigh-Taylor structures is affected, although not drastically, if the blast wave is dominated by cosmic rays. The amount of ejecta that makes it into the shocked interstellar medium is smaller in the latter case. If acceleration occurs at both shocks the extent of the Rayleigh-Taylor structures is similar but the reverse shock is strongly perturbed.

preprint2010arXiv

Star formation in galaxy mergers: ISM turbulence, dense gas excess, and scaling relations for disks and starbusts

Galaxy interactions and mergers play a significant, but still debated and poorly understood role in the star formation history of galaxies. Numerical and theoretical models cannot yet explain the main properties of merger-induced starbursts, including their intensity and their spatial extent. Usually, the mechanism invoked in merger-induced starbursts is a global inflow of gas towards the central kpc, resulting in a nuclear starburst. We show here, using high-resolution AMR simulations and comparing to observations of the gas component in mergers, that the triggering of starbursts also results from increased ISM turbulence and velocity dispersions in interacting systems. This forms cold gas that are denser and more massive than in quiescent disk galaxies. The fraction of dense cold gas largely increases, modifying the global density distribution of these systems, and efficient star formation results. Because the starbursting activity is not just from a global compacting of the gas to higher average surface densities, but also from higher turbulence and fragmentation into massive and dense clouds, merging systems can enter a different regime of star formation compared to quiescent disk galaxies. This is in quantitative agreement with recent observations suggesting that disk galaxies and starbursting systems are not the low-activity end and high-activity end of a single regime, but actually follow different scaling relations for their star formation.

preprint2010arXiv

Systematic uncertainties in the determination of the local dark matter density

A precise determination of the local dark matter density and an accurate control over the corresponding uncertainties are of paramount importance for Dark Matter (DM) searches. Using very recent high-resolution numerical simulations of a Milky Way like object, we study the systematic uncertainties that affect the determination of the local dark matter density based on dynamical measurements in the Galaxy. In particular, extracting from the simulation with baryons the orientation of the Galactic stellar disk with respect to the DM distribution, we study the DM density for an observer located at $\sim$8 kpc from the Galactic center {\it on the stellar disk}, $ρ_0$. This quantity is found to be always larger than the average density in a spherical shell of same radius $\barρ_0$, which is the quantity inferred from dynamical measurements in the Galaxy, and to vary in the range $ρ_0/\barρ_0=1.01-1.41$. This suggests that the actual dark matter density in the solar neighbourhood is on average 21\% larger than the value inferred from most dynamical measurements, and that the associated systematic errors are larger than the statistical errors recently discussed in the literature.

preprint2010arXiv

The circum-galactic gas around cosmologically simulated disks

We analyze the physical properties and infall rates of the circum-galactic gas around disks obtained in multi-resolved, cosmological, AMR simulations. At intermediate and low redshifts, disks are embedded into an extended, hot, tenuous corona that contributes largely in fueling the disk with non-enriched gas whereas the accretion of enriched gas from tidal streams occurs throughout episodic events. We derive an infall rate close to the disk of the same value as the one of the star formation rate in the disk and its temporal evolution as a function of galacto-centric radius nicely shows that the growth of galactic disks proceeds according to an inside-out formation scenario.

preprint2010arXiv

The driving mechanism of starbursts in galaxy mergers

We present hydrodynamic simulations of a major merger of disk galaxies, and study the ISM dynamics and star formation properties. High spatial and mass resolutions of 12pc and 4x10^4 M_sol allow to resolve cold and turbulent gas clouds embedded in a warmer diffuse phase. We compare to lower resolution models, where the multiphase ISM is not resolved and is modeled as a relatively homogeneous and stable medium. While merger-driven bursts of star formation are generally attributed to large-scale gas inflows towards the nuclear regions, we show that once a realistic ISM is resolved, the dominant process is actually gas fragmentation into massive and dense clouds and rapid star formation therein. As a consequence, star formation is more efficient by a factor of up to 10 and is also somewhat more extended, while the gas density probability distribution function (PDF) rapidly evolves towards very high densities. We thus propose that the actual mechanism of starburst triggering in galaxy collisions can only be captured at high spatial resolution and when the cooling of gas is modeled down to less than 10^3 K. Not only does our model reproduce the properties of the Antennae system, but it also explains the ``starburst mode'' revealed recently in high-redshift mergers compared to quiescent disks.

preprint2010arXiv

The formation of disc galaxies in a LCDM universe

We study the formation of disc galaxies in a fully cosmological framework using adaptive mesh refinement simulations. We perform an extensive parameter study of the main subgrid processes that control how gas is converted into stars and the coupled effect of supernovae feedback. We argue that previous attempts to form disc galaxies have been unsuccessful because of the universal adoption of strong feedback combined with high star formation efficiencies. Unless extreme amounts of energy are injected into the interstellar medium during supernovae events, these star formation parameters result in bulge dominated S0/Sa galaxies as star formation is too efficient at z~3. We show that a low efficiency of star-formation more closely models the subparsec physical processes, especially at high redshift. We highlight the successful formation of extended disc galaxies with scale lengths r_d=4-5 kpc, flat rotation curves and bulge to disc ratios of B/D~1/4. Not only do we resolve the formation of a Milky Way-like spiral galaxy, we also observe the secular evolution of the disc as it forms a pseudo-bulge. The disc properties agree well with observations and are compatible with the photometric and baryonic Tully-Fisher relations, the Kennicutt-Schmidt relation and the observed angular momentum content of spiral galaxies. We conclude that underlying small-scale star formation physics plays a larger role than previously considered in simulations of galaxy formation.

preprint2010arXiv

The Large-Scale Orientations of Disk Galaxies

We use a 380 h-1 pc resolution hydrodynamic AMR simulation of a cosmic filament to investigate the orientations of a sample of ~100 well-resolved galactic disks spanning two orders of magnitude in both stellar and halo mass. We find: (i) At z=0, there is an almost perfect alignment at a median angle of 18 deg, in the inner dark matter halo regions where the disks reside, between the spin vector of the gaseous and stellar galactic disks and that of their inner host haloes. The alignment between galaxy spin and spin of the entire host halo is however significantly weaker, ranging from a median of ~46 deg at z=1 to ~50 deg at z=0. (ii) The most massive galaxy disks have spins preferentially aligned so as to point along their host filaments. (iii) The spin of disks in lower-mass haloes shows, at redshifts above z~0.5 and in regions of low environmental density, a clear signature of alignment with the intermediate principal axis of the large-scale tidal field. This behavior is consistent with predictions of linear tidal torque theory. This alignment decreases with increasing environmental density, and vanishes in the highest density regions. Non-linear effects in the high density environments are plausibly responsible for establishing this density-alignment correlation. We expect that our numerical results provide important insights for both understanding intrinsic alignment in weak lensing from the astrophysical perspective and formation and evolution processes of galactic disks in a cosmological context.

preprint2009arXiv

3D simulations of supernova remnants evolution including non-linear particle acceleration

If a sizeable fraction of the energy of supernova remnant shocks is channeled into energetic particles (commonly identified with Galactic cosmic rays), then the morphological evolution of the remnants must be distinctly modified. Evidence of such modifications has been recently obtained with the Chandra and XMM-Newton X-ray satellites. To investigate these effects, we coupled a semi-analytical kinetic model of shock acceleration with a 3D hydrodynamic code (by means of an effective adiabatic index). This enables us to study the time-dependent compression of the region between the forward and reverse shocks due to the back reaction of accelerated particles, concomitantly with the development of the Rayleigh-Taylor hydrodynamic instability at the contact discontinuity. Density profiles depend critically on the injection level eta of particles: for eta up to about 10^-4 modifications are weak and progressive, for eta of the order of 10^-3 modifications are strong and immediate. Nevertheless, the extension of the Rayleigh-Taylor unstable region does not depend on the injection rate. A first comparison of our simulations with observations of Tycho's remnant strengthens the case for efficient acceleration of protons at the forward shock.

preprint2009arXiv

Disk formation and the origin of clumpy galaxies at high redshift

Observations of high redshift galaxies have revealed a multitude of large clumpy rapidly star-forming galaxies. Their formation scenario and their link to present day spirals is still unknown. In this Letter we perform adaptive mesh refinement simulations of disk formation in a cosmological context that are unrivalled in terms of mass and spatial resolution. We find that the so called "chain-galaxies" and "clump-clusters" are a natural outcome of early epochs of enhanced gas accretion from cold dense streams as well as tidally and ram-pressured stripped material from minor mergers and satellites. Through interaction with the hot halo gas, this freshly accreted cold gas settles into a large disk-like system, not necessarily aligned to an older stellar component, that undergoes fragmentation and subsequent star formation, forming large clumps in the mass range 10^7-10^9 M_sun. Galaxy formation is a complex process at this important epoch when most of the central baryons are being acquired through a range of different mechanisms - we highlight that a rapid mass loading epoch is required to fuel the fragmentation taking place in the massive arms in the outskirts of extended disks, an accretion mode that occurs naturally in the hierarchical assembly process at early epochs.

preprint2009arXiv

Morphological quenching of star formation: making early-type galaxies red

We point out a natural mechanism for quenching of star formation in early-type galaxies. It automatically links the color of a galaxy with its morphology and does not require gas consumption, removal or termination of gas supply. Given that star formation takes place in gravitationally unstable gas disks, it can be quenched when a disk becomes stable against fragmentation to bound clumps. This can result from the growth of a stellar spheroid, for instance by mergers. We present the concept of morphological quenching (MQ) using standard disk instability analysis, and demonstrate its natural occurrence in a cosmological simulation using an efficient zoom-in technique. We show that the transition from a stellar disk to a spheroid can be sufficient to stabilize the gas disk, quench star formation, and turn an early-type galaxy red and dead while gas accretion continues. The turbulence necessary for disk stability can be stirred up by sheared perturbations within the disk in the absence of bound star-forming clumps. While gas stripping processes are limited to dense groups and clusters, and other quenching mechanisms like AGN feedback, virial shock heating and gravitational heating, are limited to halos more massive than 10^12 Mo, the MQ can explain the appearance of red ellipticals even in less massive halos and in the field. The dense gas disks observed in some of today's red ellipticals may be the relics of this mechanism, whereas red galaxies with quenched gas disks are expected to be more frequent at high redshift.

preprint2009arXiv

Resolving Globular Cluster Formation within a Cosmological Context

We place constraints on the formation redshifts for blue globular clusters (BGCs), independent of the details of hydrodynamics and population III star formation. The observed radial distribution of BGCs in the Milky Way Galaxy suggests that they formed in biased dark matter halos at high redshift. As a result, simulations of a ~1 Mpc box up to z~10 must resolve BGC formation in LCDM. We find that most halo stars could be produced from destroyed BGCs and other low-mass clusters that formed at high redshift. We present a proof-of-concept simulation that captures the formation of globular-like star clusters.

preprint2008arXiv

Cosmic velocity--gravity relation in redshift space

We propose a simple way to estimate the parameter beta = Omega_m^(0.6)/b from three-dimensional galaxy surveys. Our method consists in measuring the relation between the cosmological velocity and gravity fields, and thus requires peculiar velocity measurements. The relation is measured *directly in redshift space*, so there is no need to reconstruct the density field in real space. In linear theory, the radial components of the gravity and velocity fields in redshift space are expected to be tightly correlated, with a slope given, in the distant observer approximation, by g / v = (1 + 6 beta / 5 + 3 beta^2 / 7)^(1/2) / beta. We test extensively this relation using controlled numerical experiments based on a cosmological N-body simulation. To perform the measurements, we propose a new and rather simple adaptive interpolation scheme to estimate the velocity and the gravity field on a grid. One of the most striking results is that nonlinear effects, including `fingers of God', affect mainly the tails of the joint probability distribution function (PDF) of the velocity and gravity field: the 1--1.5 sigma region around the maximum of the PDF is *dominated by the linear theory regime*, both in real and redshift space. This is understood explicitly by using the spherical collapse model as a proxy of nonlinear dynamics. Applications of the method to real galaxy catalogs are discussed, including a preliminary investigation on homogeneous (volume limited) `galaxy' samples extracted from the simulation with simple prescriptions based on halo and sub-structure identification, to quantify the effects of the bias between the galaxy and the total matter distibution, and of shot noise (ABRIDGED).

preprint2008arXiv

Hydrodynamical Adaptive Mesh Refinement Simulations of Disk Galaxies

To date, fully cosmological hydrodynamic disk simulations to redshift zero have only been undertaken with particle-based codes, such as GADGET, Gasoline, or GCD+. In light of the (supposed) limitations of traditional implementations of smoothed particle hydrodynamics (SPH), or at the very least, their respective idiosyncrasies, it is important to explore complementary approaches to the SPH paradigm to galaxy formation. We present the first high-resolution cosmological disk simulations to redshift zero using an adaptive mesh refinement (AMR)-based hydrodynamical code, in this case, RAMSES. We analyse the temporal and spatial evolution of the simulated stellar disks' vertical heating, velocity ellipsoids, stellar populations, vertical and radial abundance gradients (gas and stars), assembly/infall histories, warps/lopsideness, disk edges/truncations (gas and stars), ISM physics implementations, and compare and contrast these properties with our sample of cosmological SPH disks, generated with GCD+. These preliminary results are the first in our long-term Galactic Archaeology Simulation program.

preprint2008arXiv

Numerical simulations of galaxy evolution in cosmological context

Large volume cosmological simulations succeed in reproducing the large-scale structure of the Universe. However, they lack resolution and may not take into account all relevant physical processes to test if the detail properties of galaxies can be explained by the CDM paradigm. On the other hand, galaxy-scale simulations could resolve this in a robust way but do not usually include a realistic cosmological context. To study galaxy evolution in cosmological context, we use a new method that consists in coupling cosmological simulations and galactic scale simulations. For this, we record merger and gas accretion histories from cosmological simulations and re-simulate at very high resolution the evolution of baryons and dark matter within the virial radius of a target galaxy. This allows us for example to better take into account gas evolution and associated star formation, to finely study the internal evolution of galaxies and their disks in a realistic cosmological context. We aim at obtaining a statistical view on galaxy evolution from z = 2 to 0, and we present here the first results of the study: we mainly stress the importance of taking into account gas accretion along filaments to understand galaxy evolution.

preprint2001arXiv

Cosmological Hydrodynamics with Adaptive Mesh Refinement: a new high resolution code called RAMSES

A new N-body and hydrodynamical code, called RAMSES, is presented. It has been designed to study structure formation in the universe with high spatial resolution. The code is based on Adaptive Mesh Refinement (AMR) technique, with a tree based data structure allowing recursive grid refinements on a cell-by-cell basis. The N-body solver is very similar to the one developed for the ART code (Kravtsov et al. 97), with minor differences in the exact implementation. The hydrodynamical solver is based on a second-order Godunov method, a modern shock-capturing scheme known to compute accurately the thermal history of the fluid component. The accuracy of the code is carefully estimated using various test cases, from pure gas dynamical tests to cosmological ones. The specific refinement strategy used in cosmological simulations is described, and potential spurious effects associated to shock waves propagation in the resulting AMR grid are discussed and found to be negligible. Results obtained in a large N-body and hydrodynamical simulation of structure formation in a low density LCDM universe are finally reported, with 256^3 particles and 4.1 10^7 cells in the AMR grid, reaching a formal resolution of 8192^3. A convergence analysis of different quantities, such as dark matter density power spectrum, gas pressure power spectrum and individual haloes temperature profiles, shows that numerical results are converging down to the actual resolution limit of the code, and are well reproduced by recent analytical predictions in the framework of the halo model.

preprint2000arXiv

Numerical and Analytical Predictions for the Large-Scale Sunyaev-Zel'dovich Effect

The hot gas embedded in the large-scale structures in the Universe produces secondary fluctuations in the Cosmic Microwave Background (CMB). Because it is proportional to the gas pressure integrated along the line of sight, this effect, the thermal Sunyaev-Zel'dovich (SZ) effect, provides a direct measure of large-scale structure and of cosmological parameters. We study the statistical properties of this effect using both hydrodynamical simulations and analytical predictions from an extended halo model. The Adaptive Mesh Refinement scheme, used in the newly developed code RAMSES, provides a dynamic range of 4 order of magnitudes, and thus allows us to significantly improve upon earlier calculations. After accounting for the finite mass resolution and box size of the simulation, we find that the halo model agrees well with the simulations. We discuss and quantify the uncertainty in both methods, and thus derive an accurate prediction for the SZ power spectrum in the 10^2<l<10^5 range of multipole. We show how this combined analytical and numerical approach is essential for accuracy, and useful for the understanding of the physical processes and scales which contribute to large-scale SZ anisotropies.

preprint1998arXiv

Non-linear dynamics and mass function of cosmic structures. II Numerical results

We compare the mass functions obtained analytically, in the framework of an extended Press & Schechter (PS) formalism, in a previous paper to the numerical mass functions obtained in N -- body simulations, using different algorithms to define objects in the density field. After discussing the properties of the algorithms, we show that the mass function obtained using the friend -- of -- friend algorithm reproduces best the scaling behaviors predicted in the extended PS formalism. Following this statistical analysis, we show that it is possible in the framework of our extended PS formalism to reproduce the mass function but also, and for the first time, the initial statistical properties of structures and their collapse time. This allow to present a ``coherent'' picture of structure formation which can account for the initial, final and dynamical properties of structures.

preprint1997arXiv

Non-linear Dynamics and Mass Function of Cosmic Structures: I Analytical Results

We investigate some modifications to the Press & Schechter (1974) (PS) prescription resulting from shear and tidal effects. These modifications rely on more realistic treatments of the collapse process than the standard approach based on the spherical model. First, we show that the mass function resulting from a new approximate Lagrangian dynamic (Audit & Alimi 96), contains more objects at high mass, than the classical PS mass function and is well fitted by a PS-like function with a threshold density of $δ_c \simeq 1.4$. However, such a Lagrangian description can underestimate the epoch of structure formation since it defines it as the collapse of the first principal axis. We therefore suggest some analytical prescriptions, for computing the collapse time along the second and third principal axes, and we deduce the corresponding mass functions. The collapse along the third axis is delayed by the shear and the number of objects of high mass then decreases. Finally, we show that the shear also strongly affects the formation of low-mass halos. This dynamical effect implies a modification of the low-mass slope of the mass function and allows the reproduction of the observed luminosity function of field galaxies. In a companion paper, we present results of numerical simulations which complete this work.

Romain Teyssier

What is connected

Connect this record

See the researcher in context

Building this map preview

97 published item(s)

A Subgrid Turbulent Mean Field Dynamo Model for Cosmological Galaxy Formation Simulations

Cosmological simulations of the same spiral galaxy: the impact of baryonic physics

Cosmology with one galaxy?

Internal dark matter structure of the most massive galaxy clusters since redshift 1

On the origin of surprisingly cold gas discs in galaxies at high redshift

Robust field-level inference with dark matter halos

The driving mode of shock-driven turbulence

The Nature of High ${\rm [OIII]}_{\rm 88μm}$/${\rm [CII]}_{\rm 158μm}$ Galaxies in the Epoch of Reionization: Low Carbon Abundance and a Top-Heavy IMF?

The short ionizing photon mean free path at z=6 in Cosmic Dawn III, a new fully-coupled radiation-hydrodynamical simulation of the Epoch of Reionization

Toward fully compressible numerical simulations of stellar magneto-convection with the RAMSES code

Turbulence Generation by Shock Interaction with a Highly Non-Uniform Medium

Forming early-type galaxies without AGN feedback: a combination of merger-driven outflows and inefficient star formation

Parameter inference with non-linear galaxy clustering: accounting for theoretical uncertainties

Baryonic effects for weak lensing. Part I. Power spectrum and covariance matrix

Baryonic effects for weak lensing. Part II. Combination with X-ray data and extended cosmologies

Galactic ionising photon budget during the Epoch of Reionisation in the Cosmic Dawn II simulation

Rapid Filamentary Accretion as the Origin of Extended Thin Discs

The AGORA high-resolution galaxy simulations comparison project: Public data release

On the origin of the peak of the stellar initial mass function: exploring the tidal screening theory

A new method to quantify the effects of baryons on the matter power spectrum

Cosmic Dawn (CoDa): the First Radiation-Hydrodynamics Simulation of Reionization and Galaxy Formation in the Local Universe

Matter power spectrum and the challenge of percent accuracy

nIFTY galaxy cluster simulations III: The Similarity & Diversity of Galaxies & Subhaloes

nIFTy Galaxy Cluster simulations IV: Quantifying the Influence of Baryons on Halo Properties

nIFTy galaxy cluster simulations V: Investigation of the Cluster Infall Region

PKDGRAV3: Beyond Trillion Particle Cosmological Simulations for the Next Era of Galaxy Surveys

RHAPSODY-G simulations II - Baryonic growth and metal enrichment in massive galaxy clusters

Black hole evolution: I. Supernova-regulated black hole growth

Distribution of streaming rates into high-redshift galaxies

Galaxies that Shine: radiation-hydrodynamical simulations of disk galaxies

Galaxy Evolution: Modeling the Role of Non-thermal Pressure in the Interstellar medium

Milking the spherical cow: on aspherical dynamics in spherical coordinates

nIFTy galaxy cluster simulations II: radiative models

PHEW: a parallel segmentation algorithm for three-dimensional AMR datasets - application to structure detection in self-gravitating flows

Rhapsody-G simulations: galaxy clusters as baryonic closed boxes and the covariance between hot gas and galaxies

A fast, robust, and simple implicit method for adaptive time-stepping on adaptive mesh-refinement grids

A Systematic Look at the Effects of Radiative Feedback on Disc Galaxy Formation

Baryonic and dark matter distribution in cosmological simulations of spiral galaxies

Baryonic effects on weak-lensing two-point statistics and its cosmological implications

Brightest Cluster Galaxies in Cosmological Simulations with Adaptive Mesh Refinement: Successes and Failures

Chameleon $f(R)$ gravity on the Virgo cluster scale

Dancing in the dark: galactic properties trace spin swings along the cosmic web

Numerical cosmology on the GPU with Enzo and Ramses

The biasing of baryons on the cluster mass function and cosmological parameter estimation

The role of turbulence in star formation laws and thresholds

Towards a more realistic sink particle algorithm for the RAMSES code

Beyond the nuclear starburst? Clustered star formation in major mergers

RAMSES-RT: Radiation hydrodynamics in the cosmological context

Simulating gamma-ray binaries with a relativistic extension of RAMSES

The long lives of giant clumps and the birth of outflows in gas-rich galaxies at high redshift

Volume Rendering of AMR Simulations

A diversity of progenitors and histories for isolated spiral galaxies

Co-planar streams, pancakes, and angular-momentum exchange in high-z disc galaxies

Cusp-core transformations in dwarf galaxies: observational predictions

Hydrodynamics of galaxy mergers with supermassive black holes: is there a last parsec problem ?

Indirect dark matter searches: towards a consistent top-bottom approach for studying the gamma-ray signals and associated backgrounds

Radiatively Efficient Magnetized Bondi Accretion

RAMSES-CH: A New Chemodynamical Code for Cosmological Simulations

The effects of baryon physics, black holes and AGN feedback on the mass distribution in clusters of galaxies

The formation of the brightest cluster galaxies in cosmological simulations: the case for AGN feedback

A simple multigrid scheme for solving the Poisson equation with arbitrary domain boundaries

AGN feedback using AMR cosmological simulations

Black Hole growth and AGN obscuration by instability-driven inflows in high-redshift disk galaxies fed by cold streams

Comparing Numerical Methods for Isothermal Magnetized Supersonic Turbulence

Disc Heating: Comparing the Milky Way with Cosmological Simulations

ECOSMOG: An Efficient Code for Simulating Modified Gravity

Grid of Lya radiation transfer models for the interpretation of distant galaxies

How AGN feedback and metal cooling shape cluster entropy profiles

Hydrodynamics of high-redshift galaxy collisions: From gas-rich disks to dispersion-dominated mergers and compact spheroids

Radiation hydrodynamics with Adaptive Mesh Refinement and application to prestellar core collapse. I Methods

Self-regulated growth of supermassive black holes by a dual jet/heating AGN feedback mechanism: methods, tests and implications for cosmological simulations

The impact of ISM turbulence, clustered star formation and feedback on galaxy mass assembly through cold flows and mergers

Cell Count Moments in the Halo Model

Galactic star formation in parsec-scale resolution simulations