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.

We study in this paper three variants of the high-order Discontinuous Galerkin (DG) method with Runge-Kutta (RK) time integration for the induction equation, analysing their ability to preserve the divergence free constraint of the magnetic field. To quantify divergence errors, we use a norm based on both a surface term, measuring global divergence errors, and a volume term, measuring local divergence errors. This leads us to design a new, arbitrary high-order numerical scheme for the induction equation in multiple space dimensions, based on a modification of the Spectral Difference (SD) method [1] with ADER time integration [2]. It appears as a natural extension of the Constrained Transport (CT) method. We show that it preserves $\nabla\cdot\vec{B}=0$ exactly by construction, both in a local and a global sense. We compare our new method to the 3 RKDG variants and show that the magnetic energy evolution and the solution maps of our new SD-ADER scheme are qualitatively similar to the RKDG variant with divergence cleaning, but without the need for an additional equation and an extra variable to control the divergence errors. [1] Liu Y., Vinokur M., Wang Z.J. (2006) Discontinuous Spectral Difference Method for Conservation Laws on Unstructured Grids. In: Groth C., Zingg D.W. (eds) Computational Fluid Dynamics 2004. Springer, Berlin, Heidelberg [2] Dumbser M., Castro M., Par\'es C., Toro E.F (2009) ADER schemes on unstructured meshes for nonconservative hyperbolic systems: Applications to geophysical flows. In: Computers & Fluids, Volume 38, Issue 9

We use a recently introduced statistic called {\em Integrated Bispectrum} (IB) to probe the gravity-induced non-Gaussianity at the level of the bispectrum from weak lensing convergence or $\kappa$ maps. We generalize the concept of the IB to spherical coordinates. This result is next connected to the response function approach. Finally, we use the Euclid Flagship simulations to compute the IB as a function of redshift and wave number. We also outline how the IB can be computed using a variety of analytical approaches including the ones based on Effective Field Theory (EFT), {\em Halo models} and models based on the {\em Separate Universe approach} in projection or two-dimension (2D). Comparing these results against simulations we find that the existing theoretical models tend to over-predict the numerical value of the IB. We emphasize the role of the finite volume effect in the numerical estimation of the IB. We introduced the concept of squeezed and collapsed tripsectrum for 2D $\kappa$ maps. We derive the IB for many parameterized theories of modified gravity including the Horndeskii and beyond-Horndeskii theories specifically for the non-degenerate scenarios that are also known as the Gleyzes-Langlois-Piazza-Venizzi or GPLV theories. In addition, the cosmological models with clustering quintessence and models involving massive neutrinos are also derived.

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.

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 $\Lambda$CDM cosmology with massive neutrinos, and we consider several extensions of the currently favoured cosmological model. For the standard $\Lambda$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 $\Omega_m$ and $\sigma_8$. As extensions of $\Lambda$CDM, we consider the cases of a dynamical dark energy model (wCDM), a $f(R)$ gravity model (fRCDM), and a mixed dark matter model ($\Lambda$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 $\Lambda$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.

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 $\Lambda$CDM cosmology with both cosmological and astrophysical parameters. Our main findings are the following: (i) ignoring baryonic effects leads to a greater than 5$\sigma$ bias on the cosmological parameters $\Omega_m$ and $\sigma_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 $\Omega_m$-$\sigma_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.

Galactic outflows driven by stellar feedback are crucial for explaining the inefficiency of galaxy formation. 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. In this work, 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. At $z>1.5$, the two galaxies feature a surface density of star formation $\Sigma_{\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 epoch we dub "the Grand Twirl". After this episode, the evolutionary tracks diverge, 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 $\Sigma_{\rm SFR}$ by over two orders of magnitude with $\Sigma_{\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 in the solar neighbourhood. 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.

Stars in globular clusters (GCs) lose a non negligible amount of mass during their post-main sequence evolution. This material is then expected to build up a substantial intra-cluster medium (ICM) within the GC. However, the observed gas content in GCs is a couple of orders of magnitude below these expectations. Here we follow the evolution of this stellar wind material through hydrodynamical simulations to attempt to reconcile theoretical predictions with observations. We test different mechanisms proposed in the literature to clear out the gas such as ram-pressure stripping by the motion of the GC in the Galactic halo medium and ionisation by UV sources. We use the code ramses to run 3D hydrodynamical simulations to study for the first time the ICM evolution within discretised multi-mass GC models including stellar winds and full radiative transfer. We find that the inclusion of both ram-pressure and ionisation is mandatory to explain why only a very low amount of ionised gas is observed in the core of GCs. The same mechanisms operating in ancient GCs that clear the gas could also be efficient at younger ages, meaning that young GCs would not be able to retain gas and form multiple generations of stars as assumed in many models to explain "multiple populations". However, this rapid clearing of gas is consistent with observations of young massive clusters.

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.

Dark matter dominates the properties of large cosmological structures such as galaxy clusters, and the mass profiles of the dark matter have been measured for these equilibrated structures for years using X-rays, lensing or galaxy velocities. A new method has been proposed, which should allow us to estimate a dynamical property of the dark matter, namely the velocity anisotropy. For the gas a similar velocity anisotropy is zero due to frequent collisions, however, the collisionless nature of dark matter allows it to be non-trivial. Numerical simulations have for years found non-zero and radially varying dark matter velocity anisotropies. Here we employ the method proposed by Hansen and Pifaretti (2007), and developed by Host et al. (2009) to estimate the dark matter velocity anisotropy in the bright galaxy cluster Perseus, to near 5 times the radii previously obtained. We find the dark matter velocity anisotropy to be consistent with the results of numerical simulations, however, still with large error-bars. At half the virial radius we find the velocity anisotropy to be non-zero at 1.7 standard deviations, lending support to the collisionless nature of dark matter.

We introduce the "Engineering Dwarfs at Galaxy formation's Edge" (EDGE) project to study the cosmological formation and evolution of the smallest galaxies in the Universe. In this first paper, we explore the effects of resolution and sub-grid physics on a single low mass halo ($M_{\rm halo}=10^{9}~M_\odot$), simulated to redshift $z=0$ at a mass and spatial resolution of $\sim 20~M_\odot$ and $\sim 3$ pc. We consider different star formation prescriptions, supernova feedback strengths and on-the-fly radiative transfer (RT). We show that RT changes the mode of galactic self-regulation at this halo mass, suppressing star formation by causing the interstellar and circumgalactic gas to remain predominantly warm ($\sim 10^4$ K) even before cosmic reionisation. By contrast, without RT, star formation regulation occurs only through starbursts and their associated vigorous galactic outflows. In spite of this difference, the entire simulation suite (with the exception of models without any feedback) matches observed dwarf galaxy sizes, velocity dispersions, $V$-band magnitudes and dynamical mass-to-light-ratios. This is because such structural scaling relations are predominantly set by the host dark matter halo, with the remaining model-to-model variation being smaller than the observational scatter. We find that only the stellar mass-metallicity relation differentiates the galaxy formation models. Explosive feedback ejects more metals from the dwarf, leading to a lower metallicity at a fixed stellar mass. We conclude that the stellar mass-metallicity relation of the very smallest galaxies provides a unique constraint on galaxy formation physics.

Cosmic Dawn II (CoDa II) is a new, fully-coupled radiation-hydrodynamics simulation of cosmic reionization and galaxy formation and their mutual impact, to redshift $z < 6$. With $4096^3$ particles and cells in a 94 Mpc box, it is large enough to model global reionization and its feedback on galaxy formation while resolving all haloes above $10^8$ M$_{\odot}$. Using the same hybrid CPU-GPU code RAMSES-CUDATON as CoDa I in Ocvirk et al. (2016), CoDa II modified and re-calibrated the subgrid star-formation algorithm, making reionization end earlier, at $z \gtrsim 6$, thereby better matching the observations of intergalactic Lyman-alpha opacity from quasar spectra and electron-scattering optical depth from cosmic microwave background fluctuations. CoDa II predicts a UV continuum luminosity function in good agreement with observations of high-z galaxies, especially at $z = 6$. As in CoDa I, reionization feedback suppresses star formation in haloes below $\sim 2 \times 10^9$ M$_{\odot}$, though suppression here is less severe, a possible consequence of modifying the star-formation algorithm. Suppression is environment-dependent, occurring earlier (later) in overdense (underdense) regions, in response to their local reionization times. Using a constrained realization of $\Lambda$CDM constructed from galaxy survey data to reproduce the large-scale structure and major objects of the present-day Local Universe, CoDa II serves to model both global and local reionization. In CoDa II, the Milky Way and M31 appear as individual islands of reionization, i.e. they were not reionized by the progenitor of the Virgo cluster, nor by nearby groups, nor by each other.

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.

The cold molecular gas in contemporary galaxies is structured in discrete cloud complexes. These giant molecular clouds (GMCs), with $10^4$-$10^7$ solar masses and radii of 5-100 parsecs, are the seeds of star formation. Highlighting the molecular gas structure at such small scales in distant galaxies is observationally challenging. Only a handful of molecular clouds were reported in two extreme submillimetre galaxies at high redshift. Here we search for GMCs in a typical Milky Way progenitor at z = 1.036. Using the Atacama Large Millimeter/submillimeter Array (ALMA), we mapped the CO(4-3) emission of this gravitationally lensed galaxy at high resolution, reading down to 30 parsecs, which is comparable to the resolution of CO observations of nearby galaxies. We identify 17 molecular clouds, characterized by masses, surface densities and supersonic turbulence all of which are 10-100 times higher than present-day analogues. These properties question the universality of GMCs and suggest that GMCs inherit their properties from ambient interstellar medium. The measured cloud gas masses are similar to the masses of stellar clumps seen in the galaxy in comparable numbers. This corroborates the formation of molecular clouds by fragmentation of distant turbulent galactic gas disks, which then turn into stellar clumps ubiquitously observed in galaxies at cosmic noon.

We simulate the flux emitted from galaxy halos in order to quantify the brightness of the circumgalactic medium (CGM). We use dedicated zoom-in cosmological simulations with the hydrodynamical Adaptive Mesh Refinement code RAMSES, which are evolved down to z=0 and reach a maximum spatial resolution of 380 $h^{-1}$pc and a gas mass resolution up to 1.8$\times 10^{5} h^{-1} \rm{M}_{\odot}$ in the densest regions. We compute the expected emission from the gas in the CGM using CLOUDY emissivity models for different lines (e.g. Ly$\alpha$, CIV, OVI, CVI, OVIII) considering UV background fluorescence, gravitational cooling and continuum emission. In the case of Ly$\alpha$ we additionally consider the scattering of continuum photons. We compare our predictions to current observations and find them to be in good agreement at any redshift after adjusting the Ly$\alpha$ escape fraction. We combine our mock observations with instrument models for FIREBall-2 (UV balloon spectrograph) and HARMONI (visible and NIR IFU on the ELT) to predict CGM observations with either instrument and optimise target selections and observing strategies. Our results show that Ly$\alpha$ emission from the CGM at a redshift of 0.7 will be observable with FIREBall-2 for bright galaxies (NUV$\sim$18 mag), while metal lines like OVI and CIV will remain challenging to detect. HARMONI is found to be well suited to study the CGM at different redshifts with various tracers.

We review the numerical techniques for ideal and non-ideal magneto-hydrodynamics (MHD) used in the context of star formation simulations. We outline the specific challenges offered by modeling star forming environments, which are dominated by supersonic and super-Alfvenic turbulence in a radiative, self-gravitating fluid. These conditions are rather unique in physics and engineering and pose particularly severe restrictions on the robustness and accuracy of numerical codes. One striking aspect is the formation of collapsing fluid elements leading to the formation of singularities that represent point-like objects, namely the proto-stars. Although a few studies have attempted to resolve the formation of the first and second Larson cores, resolution limitations force us to use sink particle techniques, with sub-grid models to compute the accretion rates of mass, momentum and energy, as well as their ejection rate due to radiation and jets from the proto-stars. We discuss the most popular discretisation techniques used in the community, namely smoothed particle hydrodynamics, finite difference and finite volume methods, stressing the importance to maintain a divergence-free magnetic field. We discuss how to estimate the truncation error of a given numerical scheme, and its importance in setting the magnitude of the numerical diffusion. This can have a strong impact on the outcome of these MHD simulations, where both viscosity and resistivity are implemented at the grid scale. We then present various numerical techniques to model non-ideal MHD effects, such as Ohmic and ambipolar diffusion, as well as the Hall effect. These important physical ingredients are posing strong challenges in term of resolution and time stepping. For the latter, several strategies are discussed to overcome the limitations due to prohibitively small time steps (abridged).

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 this 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 star burst, 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.

Feedback processes from baryons are expected to strongly affect weak-lensing observables of current and future cosmological surveys. In this paper we present a new parametrisation of halo profiles based on gas, stellar, and dark matter density components. This parametrisation is used to modify outputs of gravity-only $N$-body simulations (following the prescription of Schneider and Teyssier [1]) in order to mimic baryonic effects on the matter density field. The resulting baryonic correction model relies on a few well motivated physical parameters and is able to reproduce the redshift zero clustering signal of hydrodynamical simulations at two percent accuracy below $k\sim10$ h/Mpc. A detailed study of the baryon suppression effects on the matter power spectrum and the weak lensing shear correlation reveals that the signal is dominated by two parameters describing the slope of the gas profile in haloes and the maximum radius of gas ejection. We show that these parameters can be constrained with the observed gas fraction of galaxy groups and clusters from X-ray data. Based on these observations we predict a beyond percent effect on the power spectrum above $k=0.2-1.0$ h/Mpc with a maximum suppression of 15-25 percent around $k\sim 10$ h/Mpc. As a result, the weak lensing angular shear power spectrum is suppressed by 15-25 percent at scales beyond $\ell\sim 100-600$ and the shear correlations $\xi_{+}$ and $\xi_{-}$ are affected at the 10-25 percent level below 5 and 50 arc-minutes, respectively. The relatively large uncertainties of these predictions are a result of the poorly known hydrostatic mass bias of current X-ray observations as well as the generic difficulty to observe the low density gas outside of haloes.

We present a new power spectrum emulator named EuclidEmulator that estimates the nonlinear correction to the linear dark matter power spectrum. It is based on a spectral decomposition method called polynomial chaos expansion. All steps in the construction of the emulator have been tested and optimized: the large high-resolution N-body simulations carried out with PKDGRAV3 were validated using a simulation from the Euclid Flagship campaign and demonstrated to have converged up to wavenumbers $k\approx 5\,h\,{\rm Mpc}^{-1}$ for redshifts $z\leq 5$. The emulator is constructed using the uncertainty quantification software UQLab and it has been optimized first by creating mock emulators based on Takahashi's HALOFIT. We show that it is possible to successfully predict the performance of the final emulator in this way prior to performing any N-body simulations. We provide a C-code to calculate the nonlinear correction at a relative accuracy of $\sim0.3\%$ with respect to N-body simulations within 50 ms. The absolute accuracy of the final nonlinear power spectrum is comparable to one obtained with N-body simulations, i.e. $\sim 1\%$ for $k\lesssim 1\,h\,{\rm Mpc}^{-1}$ and $z\lesssim 3.5$. This enables efficient forward modeling in the nonlinear regime allowing for maximum likelihood estimation of cosmological parameters. EuclidEmulator has been compared to HALOFIT and CosmicEmu, an alternative emulator based on the Mira-Titan Universe, and shown to be more accurate than these other approaches. This work paves a new way for optimal construction of future emulators that also consider other cosmological observables, use higher resolution input simulations and investigate higher dimensional cosmological parameter spaces.

The implementation of ACACIA (\textbf{a}nother \textbf{c}atchy \textbf{a}cronym for \textbf{c}lump-based mergertree \textbf{i}dentification in r\textbf{a}mses), an algorithm to generate dark matter halo merger trees on the fly, into the AMR code RAMSES is presented. The algorithm is fully parallel using MPI, works on the fly, and tracks dark matter substructures individually through particle IDs. Once a clump merges into another, it is still being tracked through the last identified most tightly bound particle of that clump, allowing to check at later snapshots whether the identified merging event truly was one, or whether a misidentification by the density field gradient based clump finding algorithm might have occurred, as well as to track orphan galaxies. The influence on the merger trees of various definitions of substructure and the maximal number of particles tracked per clump have been tested. Furthermore, with the known formation history of dark matter structure, galaxies can be introduced in a simulation containing only dark matter particles in a forward modelling approach through use of a parametrised stellar-mass-to-halo-mass (SMHM) relation. The obtained predicted stellar mass functions of central galaxies from $z\sim 0 - 8$ and correlation functions at $z\sim0$ are compared to observations. Considering that the mock galaxy catalogues were obtained using simulations without particularly high resolutions ($512^3$ particles), the results show satisfactory agreement with observational data. We have however demonstrated that our new merger tree algorithm and the generation of the corresponding mock galaxy catalogue can be performed on the fly, using the RAMSES built-in clump finder.

We introduce non-equilibrium molecular hydrogen chemistry into the radiation hydrodynamics code Ramses-RT. This is an adaptive mesh refinement grid code with radiation hydrodynamics that couples the thermal chemistry of hydrogen and helium to moment-based radiative transfer with the Eddington tensor closure model. The H2 physics that we include are formation on dust grains, gas phase formation, formation by three-body collisions, collisional destruction, photodissociation, photoionization, cosmic ray ionization, and self-shielding. In particular, we implement the first model for H2 self-shielding that is tied locally to moment-based radiative transfer by enhancing photodestruction. This self-shielding from Lyman-Werner line overlap is critical to H2 formation and gas cooling. We can now track the non-equilibrium evolution of molecular, atomic, and ionized hydrogen species with their corresponding dissociating and ionizing photon groups. Over a series of tests we show that our model works well compared to specialized photodissociation region codes. We successfully reproduce the transition depth between molecular and atomic hydrogen, molecular cooling of the gas, and a realistic Stromgren sphere embedded in a molecular medium. In this paper we focus on test cases to demonstrate the validity of our model on small scales. Our ultimate goal is to implement this in large-scale galactic simulations.

We present a multi-scale analysis of molecular hydrogen in a Milky Way-like simulated galaxy. Our census covers the gas content of the entire disc, to radial profiles and the Kennicutt-Schmidt relation, to a study of its molecular clouds, and finally down to a cell-by-cell analysis of the gas phases. Where observations are available we find agreement. A significant fraction of the H2 gas is in low-density regions mixed with atomic hydrogen and would therefore be difficult to observe. We use the molecular addition to ramses-rt, an adaptive mesh refinement grid code with the hydrodynamics coupled to moment-based radiative transfer. Three resolutions of the same galaxy detail the effects it has on H2 formation, with grid cells sized 97, 24, and 6.1 pc. Only the highest resolution yields gas densities high enough to host significant H2 fractions, and resolution is therefore key to simulating H2. Apart our pieces of galactic analysis are disparate, but assembled they provide a cohesive portrait of H2 in the interstellar medium. H2 chemistry on the atomic scale is sufficient to generate its dynamics throughout an entire galaxy.

Equilibrium or stationary solutions usually proceed through the exact balance between hyperbolic transport terms and source terms. Such equilibrium solutions are affected by truncation errors that prevent any classical numerical scheme from capturing the evolution of small amplitude waves of physical significance. In order to overcome this problem, we compare two commonly adopted strategies: going to very high order and reduce drastically the truncation errors on the equilibrium solution, or design a specific scheme that preserves by construction the equilibrium exactly, the so-called well-balanced approach. We present a modern numerical implementation of these two strategies and compare them in details, using hydrostatic but also dynamical equilibrium solutions of several simple test cases. Finally, we apply our methodology to the simulation of a protoplanetary disc in centrifugal equilibrium around its star and model its interaction with an embedded planet, illustrating in a realistic application the strength of both methods.

The interaction of gas-rich galaxies with the intra-cluster medium (ICM) of galaxy clusters has a remarkable impact on their evolution, mainly due to the gas loss associated with this process. In this work, we use an idealised, high-resolution simulation of a Virgo-like cluster, run with RAMSES and with dynamics reproducing that of a zoom cosmological simulation, to investigate the interaction of infalling galaxies with the ICM. We find that the tails of ram pressure stripped galaxies give rise to a population of up to more than a hundred clumps of molecular gas lurking in the cluster. The number count of those clumps varies a lot over time -- they are preferably generated when a large galaxy crosses the cluster (M$_{200c} > 10^{12}$ M$_\odot$), and their lifetime ($\lesssim 300$ Myr) is small compared to the age of the cluster. We compute the intracluster luminosity associated with the star formation which takes place within those clumps, finding that the stars formed in all of the galaxy tails combined amount to an irrelevant contribution to the intracluster light. Surprisingly, we also find in our simulation that the ICM gas significantly changes the composition of the gaseous disks of the galaxies: after crossing the cluster once, typically 20% of the cold gas still in those disks comes from the ICM.

Photoheating associated with reionization suppressed star formation in low-mass galaxies. Reionization was inhomogeneous, however, affecting different regions at different times. To establish the causal connection between reionization and suppression, we must take this local variation into account. We analyze the results of CoDa (`Cosmic Dawn') I, the first fully-coupled radiation-hydrodynamical simulation of reionization and galaxy formation in the Local Universe, in a volume large enough to model reionization globally but with enough resolving power to follow all atomic-cooling galactic halos in that volume. For every halo identified at a given time, we find the redshift at which the surrounding IGM reionized, along with its instantaneous star formation rate (`SFR') and baryonic gas-to-dark matter ratio ($M_\text{gas}/M_\text{DM}$). The average SFR per halo with $M < 10^9 \text{ M}_\odot$ was steady in regions not yet reionized, but declined sharply following local reionization. For $M > 10^{10} \text{ M}_\odot$, this SFR continued through local reionization, increasing with time, instead. For $10^9 < M < 10^{10} \text{ M}_\odot$, the SFR generally increased modestly through reionization, followed by a modest decline. In general, halo SFRs were higher for regions that reionized earlier. A similar pattern was found for $M_\text{gas}/M_\text{DM}$, which declined sharply following local reionization for $M < 10^9 \text{ M}_\odot$. Local reionization time correlates with local matter overdensity, which determines the local rates of structure formation and ionizing photon consumption. The earliest patches to develop structure and reionize ultimately produced more stars than they needed to finish and maintain their own reionization, exporting their `surplus' starlight to help reionize regions that developed structure later.

We present the SPHINX suite of cosmological adaptive mesh refinement simulations, the first radiation-hydrodynamical simulations to simultaneously capture large-scale reionization and the escape of ionizing radiation from thousands of resolved galaxies. Our $5$ and $10$ co-moving Mpc volumes resolve haloes down to the atomic cooling limit and model the inter-stellar medium with better than $\approx10$ pc resolution. The project has numerous goals in improving our understanding of reionization and making predictions for future observations. In this first paper we study how the inclusion of binary stars in computing stellar luminosities impacts reionization, compared to a model that includes only single stars. Owing to the suppression of galaxy growth via strong feedback, our galaxies are in good agreement with observational estimates of the galaxy luminosity function. We find that binaries have a significant impact on the timing of reionization: with binaries, our boxes are $99.9$ percent ionized by volume at $z\approx 7$, while without them our volumes fail to reionize by $z=6$. These results are robust to changes in volume size, resolution, and feedback efficiency. The escape of ionizing radiation from individual galaxies varies strongly and frequently. On average, binaries lead to escape fractions of $\approx 7-10$ percent, about $3.5$ times higher than with single stars only. The higher escape fraction is a result of a shallower decline in ionizing luminosity with age, and is the primary reason for earlier reionization, although the higher integrated luminosity with binaries also plays a sub-dominant role.

Arguably the main challenge of galactic magnetism studies is to explain how the interstellar medium of galaxies reaches energetic equipartition despite the extremely weak cosmic primordial magnetic fields that are originally predicted to thread the inter-galactic medium. Previous numerical studies of isolated galaxies suggest that a fast dynamo amplification might suffice to bridge the gap spanning many orders of magnitude in strength between the weak early Universe magnetic fields and the ones observed in high redshift galaxies. To better understand their evolution in the cosmological context of hierarchical galaxy growth, we probe the amplification process undergone by the cosmic magnetic field within a spiral galaxy to unprecedented accuracy by means of a suite of constrained transport magnetohydrodynamical adaptive mesh refinement cosmological zoom simulations with different stellar feedback prescriptions. A galactic turbulent dynamo is found to be naturally excited in this cosmological environment, being responsible for most of the amplification of the magnetic energy. Indeed, we find that the magnetic energy spectra of simulated galaxies display telltale inverse cascades. Overall, the amplification process can be divided in three main phases, which are related to different physical mechanisms driving galaxy evolution: an initial collapse phase, an accretion-driven phase, and a feedback-driven phase. While different feedback models affect the magnetic field amplification differently, all tested models prove to be subdominant at early epochs, before the feedback-driven phase is reached. Thus the three-phase evolution paradigm is found to be quite robust vis-a-vis feedback prescriptions.

We present a two-dimensional Cartesian code based on high order discontinuous Galerkin methods, implemented to run in parallel over multiple GPUs. A simple planet-disc setup is used to compare the behaviour of our code against the behaviour found using the FARGO3D code with a polar mesh. We make use of the time dependence of the torque exerted by the disc on the planet as a mean to quantify the numerical viscosity of the code. We find that the numerical viscosity of the Keplerian flow can be as low as a few $10^{-8}r^2\Omega$, $r$ and $\Omega$ being respectively the local orbital radius and frequency, for fifth order schemes and resolution of $\sim 10^{-2}r$. Although for a single disc problem a solution of low numerical viscosity can be obtained at lower computational cost with FARGO3D (which is nearly an order of magnitude faster than a fifth order method), discontinuous Galerkin methods appear promising to obtain solutions of low numerical viscosity in more complex situations where the flow cannot be captured on a polar or spherical mesh concentric with the disc.

Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly-formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low density region, and between more diffusive and less diffusive schemes in the high density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.

Today's galaxies experienced cosmic reionization at different times in different locations. For the first time, reionization ($50\%$ ionized) redshifts, $z_R$, at the location of their progenitors are derived from new, fully-coupled radiation-hydrodynamics simulation of galaxy formation and reionization at $z > 6$, matched to N-body simulation to z = 0. Constrained initial conditions were chosen to form the well-known structures of the local universe, including the Local Group and Virgo, in a (91 Mpc)$^3$ volume large enough to model both global and local reionization. Reionization simulation CoDa I-AMR, by CPU-GPU code EMMA, used (2048)$^3$ particles and (2048)$^3$ initial cells, adaptively-refined, while N-body simulation CoDa I-DM2048, by Gadget2, used (2048)$^3$ particles, to find reionization times for all galaxies at z = 0 with masses $M(z=0)\ge 10^8 M_\odot$. Galaxies with $M(z=0) \gtrsim 10^{11} M_\odot$ reionized earlier than the universe as a whole, by up to $\sim$ 500 Myrs, with significant scatter. For Milky-Way-like galaxies, $z_R$ ranged from 8 to 15. Galaxies with $M(z=0) \lesssim 10^{11} M_\odot$ typically reionized as late or later than globally-averaged $50\%$ reionization at $\langle z_R\rangle =7.8$, in neighborhoods where reionization was completed by external radiation. The spread of reionization times within galaxies was sometimes as large as the galaxy-to-galaxy scatter. The Milky Way and M31 reionized earlier than global reionization but later than typical for their mass, neither dominated by external radiation. Their most massive progenitors at $z>6$ had $z_R$ = 9.8 (MW) and 11 (M31), while their total masses had $z_R$ = 8.2 (both).

While feedback from Active Galactic Nuclei (AGN) is an important heating source in the centre of galaxy clusters, it is still unclear how the feedback energy is injected into the intracluster medium (ICM) and what role different numerical approaches play. Here, we compare four hydrodynamical schemes in idealized simulations of a rising bubble inflated by AGN feedback in a hot stratified ICM: (traditional) smoothed particle hydrodynamics (TSPH), a pressure flavour of SPH (PSPH), a meshless finite mass (MFM) scheme, as well as an Eulerian code with adaptive mesh refinement. In the absence of magnetic fields, the bubble is Kelvin-Helmholtz unstable on short enough time scales to dissolve it fully in the ICM, which is captured by MFM and RAMSES simulations, while in the TSPH simulation the bubble survives. When the ICM is turbulent, mixing of the bubble with the ICM is accelerated. This occurs if the numerical scheme can capture the instabilities well. The differences in the evolution of the bubble has a surprisingly small influence on the thermal structure of the ICM. However, in the simulations with MFM and RAMSES the bubble disruption leads to turbulent stirring of the ICM which is suppressed in SPH. In the latter the thermal energy remains trapped in the bubble and is transported to large radii. We discuss if the choice of hydrodynamical schemes can lead to systematic differences in the outcomes of cosmological simulations.

We have recently improved our model of active galactic nucleus (AGN) by attaching the supermassive black hole (SMBH) to a massive nuclear star cluster (NSC). Here we study the effects of this new model in massive, gas-rich galaxies with several simulations of different feedback recipes with the hydrodynamics code RAMSES. These simulations are compared to a reference simulation without any feedback, in which the cooling halo gas is quickly consumed in a burst of star formation. In the presence of strong supernovae (SN) feedback, we observe the formation of a galactic fountain that regulates star formation over a longer period, but without halting it. If only AGN feedback is considered, as soon as the SMBH reaches a critical mass, strong outflows of hot gas are launched and prevent the cooling halo gas from reaching the disk, thus efficiently halting star formation, leading to the so-called "quenching". If both feedback mechanisms act in tandem, we observe a non-linear coupling, in the sense that the dense gas in the supernovae-powered galactic fountain is propelled by the hot outflow powered by the AGN at much larger radii than without AGN. We argue that these particular outflows are able to unbind dense gas from the galactic halo, thanks to the combined effect of SN and AGN feedback. We speculate that this mechanism occurs at the end of the fast growing phase of SMBH, and is at the origin of the dense molecular outflows observed in many massive high-redshift galaxies.

Cosmic far-infrared background (CFIRB) is a powerful probe of the history of star formation rate (SFR) and the connection between baryons and dark matter across cosmic time. In this work, we explore to which extent the CFIRB anisotropies can be reproduced by a simple physical framework for galaxy evolution, the gas regulator (bathtub) model. This model is based on continuity equations for gas, stars, and metals, taking into account cosmic gas accretion, star formation, and gas ejection. We model the large-scale galaxy bias and small-scale shot noise self-consistently, and we constrain our model using the CFIRB power spectra measured by Planck. Because of the simplicity of the physical model, the goodness of fit is limited. We compare our model predictions with the observed correlation between CFIRB and gravitational lensing, bolometric infrared luminosity functions, and submillimetre source counts. The strong clustering of CFIRB indicates a large galaxy bias, which corresponds to haloes of mass $10^{12.5} M_\odot$ at z=2, higher than the mass associated with the peak of the star formation efficiency. We also find that the far-infrared luminosities of haloes above $10^{12} M_\odot$ are higher than the expectation from the SFR observed in ultraviolet and optical surveys.

Understanding the origin of strong galactic outflows and the suppression of star formation in dwarf galaxies is a key problem in galaxy formation. Using a set of radiation-hydrodynamic simulations of an isolated dwarf galaxy embedded in a $10^{10}\,M_\odot$ halo, we show that the momentum transferred from resonantly scattered Lyman-$\alpha$ (Lya) photons is an important source of stellar feedback which can shape the evolution of galaxies. We find that Lya feedback suppresses star formation by a factor of two in metal-poor galaxies by regulating the dynamics of star-forming clouds before the onset of supernova explosions (SNe). This is possible because each Lya photon resonantly scatters and imparts 10-300 times greater momentum than in the single scattering limit. Consequently, the number of star clusters predicted in the simulations is reduced by a factor of $\sim 5$, compared to the model without the early feedback. More importantly, we find that galactic outflows become weaker in the presence of strong Lya radiation feedback, as star formation and associated SNe become less bursty. We also examine a model in which radiation field is arbitrarily enhanced by a factor of up to 10, and reach the same conclusion. The typical mass loading factors in our metal-poor dwarf system are estimated to be $\sim5-10$ near the mid plane, while it is reduced to $\sim1$ at larger radii. Finally, we find that the escape of ionizing radiation and hence the reionization history of the Universe is unlikely to be strongly affected by Lya feedback.

Most stars in the Galaxy are believed to be formed within star clusters from collapsing molecular clouds. However, the complete process of star formation, from the parent cloud to a gas-free star cluster, is still poorly understood. We perform radiation-hydrodynamical simulations of the collapse of a turbulent molecular cloud using the RAMSES-RT code. Stars are modelled using sink particles, from which we self-consistently follow the propagation of the ionising radiation. We study how different feedback models affect the gas expulsion from the cloud and how they shape the final properties of the emerging star cluster. We find that the star formation efficiency is lower for stronger feedback models. Feedback also changes the high mass end of the stellar mass function. Stronger feedback also allows the establishment of a lower density star cluster, which can maintain a virial or sub-virial state. In the absence of feedback, the star formation efficiency is very high, as well as the final stellar density. As a result, high energy close encounters make the cluster evaporate quickly. Other indicators, such as mass segregation, statistics of multiple systems and escaping stars confirm this picture. Observations of young star clusters are in best agreement with our strong feedback simulation.

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.

We consider the dynamics in and near galaxy clusters. Gas, dark matter and galaxies are presently falling into the clusters between approximately 1 and 5 virial radii. At very large distances, beyond 10 virial radii, all matter is following the Hubble flow, and inside the virial radius the matter particles have on average zero radial velocity. The cosmological parameters are imprinted on the infall profile of the gas, however, no method exists, which allows a measurement of it. We consider the results of two cosmological simulations (using the numerical codes RAMSES and Gadget) and find that the gas and dark matter radial velocities are very similar. We derive the relevant dynamical equations, in particular the generalized hydrostatic equilibrium equation, including both the expansion of the Universe and the cosmological background. This generalized gas equation is the main new contribution of this paper. We combine these generalized equations with the results of the numerical simulations to estimate the contribution to the measured cluster masses from the radial velocity: inside the virial radius it is negligible, and inside two virial radii the effect is below 40%, in agreement the earlier analyses for DM. We point out how the infall velocity in principle may be observable, by measuring the gas properties to distance of about two virial radii, however, this is practically not possible today.

Magnetic fields are widely observed in the Universe in virtually all astrophysical objects, from individual stars to entire galaxies, even in the intergalactic medium, but their specific generation has long been debated. Due to the development of more realistic models of galaxy formation, viable scenarios are emerging to explain cosmic magnetism, thanks to both deeper observations and more efficient and accurate computer simulations. We present here a new cosmological high-resolution zoom-in magnetohydrodynamic (MHD) simulation, using the adaptive mesh refinement (AMR) technique, of a dwarf galaxy with an initially weak and uniform magnetic seed field that is amplified by a small-scale dynamo driven by supernova-induced turbulence. As first structures form from the gravitational collapse of small density fluctuations, the frozen-in magnetic field separates from the cosmic expansion and grows through compression. In a second step, star formation sets in and establishes a strong galactic fountain, self-regulated by supernova explosions. Inside the galaxy, the interstellar medium becomes highly turbulent, dominated by strong supersonic shocks, as demonstrated by the spectral analysis of the gas kinetic energy. In this turbulent environment, the magnetic field is quickly amplified via a small-scale dynamo process and is finally carried out into the circumgalactic medium by a galactic wind. This realistic cosmological simulation explains how initially weak magnetic seed fields can be amplified quickly in early, feedback-dominated galaxies, and predicts, as a consequence of the small scale dynamo process, that high-redshift magnetic fields are likely to be dominated by their small scale components.

Star formation from the interstellar medium of galactic disks is a basic process controlling the evolution of galaxies. Understanding the star formation rate in a local patch of a disk with a given gas mass is thus an important challenge for theoretical models. Here we simulate a kiloparsec region of a disk, following the evolution of self-gravitating molecular clouds down to subparsec scales, as they form stars that then inject feedback energy by dissociating and ionizing UV photons and supernova explosions. We assess the relative importance of each feedback mechanism. We find that $\rm H_2$-dissociating feedback results in the largest absolute reduction in star formation compared to the run with no feedback. Subsequently adding photoionization feedback produces a more modest reduction. Our fiducial models that combine all three feedback mechanisms yield, without fine-tuning, star formation rates that are in excellent agreement with observations, with $\rm H_2$-dissociating photons playing a crucial role. Models that only include supernova feedback---a common method in galaxy evolution simulations---settle to similar star formation rates, but with very different temperature and chemical states of the gas, and with very different spatial distributions of young stars.

Magnetic fields in galaxies are believed to be the result of dynamo amplification of initially weak seed fields, reaching equipartition strength inside the interstellar medium. The small-scale dynamo appears to be a viable mechanism to explain observations of strong magnetic fields in present-day and high-redshift galaxies, considering the extreme weakness of seed fields predicted by battery mechanisms or primordial fields. Performing high-resolution adaptive mesh magneto-hydrodynamic simulations of a small mass, isolated cooling halo with an initial magnetic seed field strength well below equipartition, we follow the small-scale dynamo amplification from supernova-induced turbulence up to saturation of the field. We find that saturation occurs when the average magnetic pressure reaches only 3 % to 5 % of the turbulent pressure. The magnetic energy growth transitions from exponential to linear, and finally comes to halt. The saturation level increases slightly with grid resolution. These results are in good agreement with theoretical predictions for magnetic Prandtl numbers of order $\mathrm{Pr_M} \sim 1$ and turbulent Mach numbers of order $\mathrm{M} \sim 10$. When we suppress supernova feedback after our simulation has reached saturation, we find that turbulence decays and that the gas falls back onto a thin disk with the magnetic field in local equipartition. We propose a scenario in which galactic magnetic fields are amplified from weak seed fields in the early stages of the Universe to sub-equipartition fields, owing to the turbulent environment of feedback-dominated galaxies at high redshift, and are evolved further in a later stage up to equipartition, as galaxies transformed into more quiescent, large spiral disks.

The origin and evolution of magnetic fields in the Universe is still an open question. Their observations in galaxies suggest strong magnetic fields already at high redshift as well as at present time. However, neither primordial magnetic fields nor battery processes can account for such high field strengths, which implies the presence of a dynamo process with rapid growth rates in high-redshift galaxies and subsequent maintenance against decay. We investigate the particular role played by feedback mechanisms in creating strong fluid turbulence, allowing for a magnetic dynamo to emerge. Performing magnetohydrodynamic simulations of isolated cooling gas halos, we compare the magnetic field evolution for various initial field topologies and various stellar feedback mechanisms. We find that feedback can indeed drive strong gas turbulence and dynamo action. We see typical properties of Kolmogorov turbulence with a $k^{-5/3}$ kinetic energy spectrum, as well as a small-scale dynamo, with a $k^{3/2}$ magnetic energy spectrum predicted by Kazantsev dynamo theory. We also investigate simulations with a final quiescent phase. As turbulence decreases, the galactic fountain settles into a thin, rotationally supported disk. The magnetic field develops a large-scale, well-ordered structure with even symmetry, which is in good agreement with magnetic field observations of nearby spirals. Our findings suggest that weak initial seed fields were first amplified by a small-scale dynamo during a violent, feedback-dominated early phase in the galaxy formation history, followed by a more quiescent evolution, where the fields have slowly decayed or were maintained via large-scale dynamo action.

We introduce a new model for the formation and evolution of supermassive black holes (SMBHs) in the RAMSES code using sink particles, improving over previous work the treatment of gas accretion and dynamical evolution. This new model is tested against a suite of high-resolution simulations of an isolated, gas-rich, cooling halo. We study the effect of various feedback models on the SMBH growth and its dynamics within the galaxy. In runs without any feedback, the SMBH is trapped within a massive bulge and is therefore able to grow quickly, but only if the seed mass is chosen larger than the minimum Jeans mass resolved by the simulation. We demonstrate that, in the absence of supernovae (SN) feedback, the maximum SMBH mass is reached when Active Galactic Nucleus (AGN) heating balances gas cooling in the nuclear region. When our efficient SN feedback is included, it completely prevents bulge formation, so that massive gas clumps can perturb the SMBH orbit, and reduce the accretion rate significantly. To overcome this issue, we propose an observationally motivated model for the joint evolution of the SMBH and a parent nuclear star cluster (NSC), which allows the SMBH to remain in the nuclear region, grow fast and resist external perturbations. In this scenario, however, SN feedback controls the gas supply and the maximum SMBH mass now depends on the balance between AGN heating and gravity. We conclude that SMBH/NSC co-evolution is crucial for the growth of SMBH in high-z galaxies, the progenitors of massive elliptical today.

The effect of baryons on the matter power spectrum is likely to have an observable effect for future galaxy surveys, like Euclid or LSST. As a first step towards a fully predictive theory, we investigate the effect of non-radiative hydrodynamics on the structure of galaxy groups sized halos, which contribute the most to the weak lensing power spectrum. We perform high resolution (more than one million particles per halo and one kilo-parsec resolution) non-radiative hydrodynamical zoom-in simulations of a sample of 16 halos, comparing the profiles to popular analytical models. We find that the total mass profile is well fitted by a Navarro, Frenk & White model, with parameters slightly modified from the dark matter only simulation. We also find that the Komatsu & Seljak hydrostatic solution provides a good fit to the gas profiles, with however significant deviations, arising from strong turbulent mixing in the core and from non-thermal, turbulent pressure support in the outskirts. The turbulent energy follows a shallow, rising linear profile with radius, and correlates with the halo formation time. Using only three main structural halo parameters as variables (total mass, concentration parameter and central gas density), we can predict with an accuracy better than 20% the individual gas density and temperature profiles. For the average total mass profile, which is relevant for power spectrum calculations, we even reach an accuracy of 1%. The robustness of these predictions has been tested against resolution effects, different types of initial conditions and hydrodynamical schemes.

We compare 5 sub-grid models for supernova (SN) feedback in adaptive mesh refinement (AMR) simulations of isolated dwarf and L-star disk galaxies with 20-40 pc resolution. The models are thermal dump, stochastic thermal, 'mechanical' (injecting energy or momentum depending on the resolution), kinetic, and delayed cooling feedback. We focus on the ability of each model to suppress star formation and generate outflows. Our highest-resolution runs marginally resolve the adiabatic phase of the feedback events, which correspond to 40 SN explosions, and the first three models yield nearly identical results, possibly indicating that kinetic and delayed cooling feedback converge to wrong results. At lower resolution all models differ, with thermal dump feedback becoming inefficient. Thermal dump, stochastic, and mechanical feedback generate multiphase outflows with mass loading factors $\beta \ll 1$, which is much lower than observed. For the case of stochastic feedback we compare to published SPH simulations, and find much lower outflow rates. Kinetic feedback yields fast, hot outflows with $\beta\sim 1$, but only if the wind is in effect hydrodynamically decoupled from the disk by using a large bubble radius. Delayed cooling generates cold, dense and slow winds with $\beta> 1$, but large amounts of gas occupy regions of temperature-density space with short cooling times. We conclude that either our resolution is too low to warrant physically motivated models for SN feedback, that feedback mechanisms other than SNe are important, or that other aspects of galaxy evolution, such as star formation, require better treatment.

Galaxy mergers are believed to trigger strong starbursts. This is well assessed by observations in the local Universe. However the efficiency of this mechanism has poorly been tested so far for high redshift, actively star forming, galaxies. We present a suite of pc-resolution hydrodynamical numerical simulations to compare the star formation process along a merging sequence of high and low z galaxies, by varying the gas mass fraction between the two models. We show that, for the same orbit, high-redshift gas-rich mergers are less efficient than low-redshift ones at producing starbursts: the star formation rate excess induced by the merger and its duration are both around 10 times lower than in the low gas fraction case. The mechanisms that account for the star formation triggering at low redshift - the increased compressive turbulence, gas fragmentation, and central gas inflows - are only mildly, if not at all, enhanced for high gas fraction galaxy encounters. Furthermore, we show that the strong stellar feedback from the initially high star formation rate in high redshift galaxies does not prevent an increase of the star formation during the merger. Our results are consistent with the observed increase of the number of major mergers with increasing redshift being faster than the respective increase in the number of starburst galaxies.

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.

We examine the properties of the galaxies and dark matter haloes residing in the cluster infall region surrounding the simulated $\Lambda$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.

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)

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.

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.

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.

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.

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.

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 $\Lambda$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.

We have implemented non-ideal Magneto-Hydrodynamics (MHD) effects in the Adaptive Mesh Refinement (AMR) code RAMSES, namely ambipolar diffusion and Ohmic dissipation, as additional source terms in the ideal MHD equations. We describe in details how we have discretized these terms using the adaptive Cartesian mesh, and how the time step is diminished with respect to the ideal case, in order to perform a stable time integration. We have performed a large suite of test runs, featuring the Barenblatt diffusion test, the Ohmic diffusion test, the C-shock test and the Alfven wave test. For the latter, we have performed a careful truncation error analysis to estimate the magnitude of the numerical diffusion induced by our Godunov scheme, allowing us to estimate the spatial resolution that is required to address non-ideal MHD effects reliably. We show that our scheme is second-order accurate, and is therefore ideally suited to study non-ideal MHD effects in the context of star formation and molecular cloud dynamics.

The strong time-dependence of the dynamics of galactic bars yields a complex and rapidly evolving distribution of dense gas and star forming regions. Although bars mainly host regions void of any star formation activity, their extremities can gather the physical conditions for the formation of molecular complexes and mini-starbursts. Using a sub-parsec resolution hydrodynamical simulation of a Milky Way-like galaxy, we probe these conditions to explore how and where bar (hydro-)dynamics favours the formation or destruction of molecular clouds and stars. The interplay between the kpc-scale dynamics (gas flows, shear) and the parsec-scale (turbulence) is key to this problem. We find a strong dichotomy between the leading and trailing sides of the bar, in term of cloud fragmentation and in the age distribution of the young stars. After orbiting along the bar edge, these young structures slow down at the extremities of the bar, where orbital crowding increases the probability of cloud-cloud collision. We find that such events increase the Mach number of the cloud, leading to an enhanced star formation efficiency and finally the formation of massive stellar associations, in a fashion similar to galaxy-galaxy interactions. We highlight the role of bar dynamics in decoupling young stars from the clouds in which they form, and discuss the implications on the injection of feedback into the interstellar medium, in particular in the context of galaxy formation.

We present the Rhapsody-G suite of cosmological hydrodynamic AMR zoom simulations of ten massive galaxy clusters at the $M_{\rm vir}\sim10^{15}\,{\rm M}_\odot$ scale. These simulations include cooling and sub-resolution models for star formation and stellar and supermassive black hole feedback. The sample is selected to capture the whole gamut of assembly histories that produce clusters of similar final mass. We present an overview of the successes and shortcomings of such simulations in reproducing both the stellar properties of galaxies as well as properties of the hot plasma in clusters. In our simulations, a long-lived cool-core/non-cool core dichotomy arises naturally, and the emergence of non-cool cores is related to low angular momentum major mergers. Nevertheless, the cool-core clusters exhibit a low central entropy compared to observations, which cannot be alleviated by thermal AGN feedback. For cluster scaling relations we find that the simulations match well the $M_{500}-Y_{500}$ scaling of Planck SZ clusters but deviate somewhat from the observed X-ray luminosity and temperature scaling relations in the sense of being slightly too bright and too cool at fixed mass, respectively. Stars are produced at an efficiency consistent with abundance matching constraints and central galaxies have star formation rates consistent with recent observations. While our simulations thus match various key properties remarkably well, we conclude that the shortcomings strongly suggest an important role for non-thermal processes (through feedback or otherwise) or thermal conduction in shaping the intra-cluster medium.

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.

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 (\Delta_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.

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.

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.

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.

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.

We describe and test an updated version of radiation-hydrodynamics (RHD) in the RAMSES code, that includes three new features: i) radiation pressure on gas, ii) accurate treatment of radiation diffusion in an unresolved optically thick medium, and iii) relativistic corrections that account for Doppler effects and work done by the radiation to first order in v/c. We validate the implementation in a series of tests, which include a morphological assessment of the M1 closure for the Eddington tensor in an astronomically relevant setting, dust absorption in a optically semi-thick medium, direct pressure on gas from ionising radiation, convergence of our radiation diffusion scheme towards resolved optical depths, correct diffusion of a radiation flash and a constant luminosity radiation, and finally, an experiment from Davis et al. of the competition between gravity and radiation pressure in a dusty atmosphere, and the formation of radiative Rayleigh-Taylor instabilities. With the new features, RAMSES-RT can be used for state-of-the-art simulations of radiation feedback from first principles, on galactic and cosmological scales, including not only direct radiation pressure from ionising photons, but also indirect pressure via dust from multi-scattered IR photons reprocessed from higher-energy radiation, both in the optically thin and thick limits.

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.

We model the intensity of emission lines from the CO molecule, based on hydrodynamic simulations of spirals, mergers, and high-redshift galaxies with very high resolutions (3pc and 10^3 Msun) and detailed models for the phase-space structure of the interstellar gas including shock heating, stellar feedback processes and galactic winds. The simulations are analyzed with a Large Velocity Gradient (LVG) model to compute the local emission in various molecular lines in each resolution element, radiation transfer and opacity effects, and the intensity emerging from galaxies, to generate synthetic spectra for various transitions of the CO molecule. This model reproduces the known properties of CO spectra and CO-to-H2 conversion factors in nearby spirals and starbursting major mergers. The high excitation of CO lines in mergers is dominated by an excess of high-density gas, and the high turbulent velocities and compression that create this dense gas excess result in broad linewidths and low CO intensity-to-H2 mass ratios. When applied to high-redshift gas-rich disks galaxies, the same model predicts that their CO-to-H2 conversion factor is almost as high as in nearby spirals, and much higher than in starbursting mergers. High-redshift disk galaxies contain giant star-forming clumps that host a high-excitation component associated to gas warmed by the spatially-concentrated stellar feedback sources, although CO(1-0) to CO(3-2) emission is overall dominated by low-excitation gas around the densest clumps. These results overall highlight a strong dependence of CO excitation and the CO-to-H2 conversion factor on galaxy type, even at similar star formation rates or densities. The underlying processes are driven by the interstellar medium structure and turbulence and its response to stellar feedback, which depend on global galaxy structure and in turn impact the CO emission properties.

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.

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.

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$\sigma$. 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.

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.

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.

Models of modified gravity offer promising alternatives to the concordance $\Lambda$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 $\Lambda$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.

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.

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.

Metal poor globular clusters (MPGCs) are a unique probe of the early universe, in particular the reionization era. Systems of globular clusters in galaxy clusters are particularly interesting as it is in the progenitors of galaxy clusters that the earliest reionizing sources first formed. Although the exact physical origin of globular clusters is still debated, it is generally admitted that globular clusters form in early, rare dark matter peaks (Moore et al. 2006; Boley et al. 2009). We provide a fully numerical analysis of the Virgo cluster globular cluster system by identifying the present day globular cluster system with exactly such early, rare dark matter peaks. A popular hypothesis is that that the observed truncation of blue metal poor globular cluster formation is due to reionization (Spitler et al. 2012; Boley et al. 2009; Brodie & Strader 2006); adopting this view, constraining the formation epoch of MPGCs provides a complementary constraint on the epoch of reionization. By analyzing both the line of sight velocity dispersion and the surface density distribution of the present day distribution we are able to constrain the redshift and mass of the dark matter peaks. We find and quantify a dependence on the chosen line of sight of these quantities, whose strength varies with redshift, and coupled with star formation efficiency arguments find a best fitting formation mass and redshift of $\simeq 5 \times 10^8 \rm{M}_\odot$ and $z\simeq 9$. We predict $\simeq 300$ intracluster MPGCs in the Virgo cluster. Our results confirm the techniques pioneered by Moore et al. (2006) when applied to the the Virgo cluster and extend and refine the analytic results of Spitler et al. (2012) numerically.

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 $\sigma_8$ and $\Omega_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.

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.

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.

We introduce the AGORA project, a comprehensive numerical study of well-resolved galaxies within the LCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of 8 galaxies with halo masses M_vir ~= 1e10, 1e11, 1e12, and 1e13 Msun at z=0 and two different ("violent" and "quiescent") assembly histories. The numerical techniques and implementations used in this project include the smoothed particle hydrodynamics codes GADGET and GASOLINE, and the adaptive mesh refinement codes ART, ENZO, and RAMSES. The codes will share common initial conditions and common astrophysics packages including UV background, metal-dependent radiative cooling, metal and energy yields of supernovae, and stellar initial mass function. These are described in detail in the present paper. Subgrid star formation and feedback prescriptions will be tuned to provide a realistic interstellar and circumgalactic medium using a non-cosmological disk galaxy simulation. Cosmological runs will be systematically compared with each other using a common analysis toolkit, and validated against observations to verify that the solutions are robust - i.e., that the astrophysical assumptions are responsible for any success, rather than artifacts of particular implementations. The goals of the AGORA project are, broadly speaking, to raise the realism and predictive power of galaxy simulations and the understanding of the feedback processes that regulate galaxy "metabolism." The proof-of-concept dark matter-only test of the formation of a galactic halo with a z=0 mass of M_vir ~= 1.7e11 Msun by 9 different versions of the participating codes is also presented to validate the infrastructure of the project.

We aim at addressing the questions related to galaxy mass assembly through major and minor wet merging processes in the redshift range 1<z<2. A consequent fraction of Milky Way like galaxies are thought to have undergone an unstable clumpy phase at this early stage. Using the adaptive mesh refinement code RAMSES, with a recent physically-motivated implementation of stellar feedback, we build the Merging and Isolated high-Redshift Adaptive mesh refinement Galaxies (MIRAGE) sample. It is composed of 20 mergers and 3 isolated idealized disks simulations with global physical properties in accordance with the 1<z<2 mass complete sample MASSIV. The numerical hydrodynamical resolution reaches 7 parsecs in the smallest Eulerian cells. Our simulations include: star formation, metal line cooling, metallicity advection, and a recent implementation of stellar feedback which encompasses OB-type stars radiative pressure, photo-ionization heating, and supernovae. The initial conditions are set to match the z~2 observations, thanks to a new public code DICE. The numerical resolution allows us to follow the formation and evolution of giant clumps formed in-situ from Jeans instabilities triggered by high initial gas fraction. The star formation history of isolated disks shows stochastic star formation rate, which proceeds from the complex behavior of the giant clumps. Our minor and major gas-rich merger simulations do not trigger starbursts, suggesting a saturation of the star formation in a turbulent and clumpy interstellar medium fed by substantial accretion from the circum-galactic medium. Our simulations are close to the normal regime of the disk-like star formation on a Schmidt-Kennicutt diagram. The mass-size relation and its rate of evolution matches observations, suggesting that the inside-out growth mechanisms of the stellar disk do not necessarily require to be achieved through a cold accretion.

The organization of the interstellar medium in disk galaxies obeys the large scale dynamics set by kpc-size structures. Improving our knowledge of how the dense, molecular gas is distributed in a disk is an important step in our understanding of star formation at galactic scale. Using a recently published simulation of a Milky Way-like galaxy, we explore the formation and dynamical organization of the star forming gas in a proto-typical disk. Along spiral arms, we report the formation of regularly spaced clouds, called beads on a string and spurs. The former form through gravitational instabilities while the later originate from Kelvin-Helmholtz instabilities. We propose that the co-existence of both structures in the same galaxy can be explained by a different role of the disk dynamics, depending on the location within the disk. In particular, we highlight the impact of the pitch angle of the spiral arm in the development of either type of structure.

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.

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.

We present a self-consistent hydrodynamical simulation of a Milky Way-like galaxy, at the resolution of 0.05 pc. The model includes star formation and a new implementation of stellar feedback through photo-ionization, radiative pressure and supernovae. The simulation resolves the structure of the interstellar medium at subparsec resolution for a few cloud lifetimes, and at 0.05 pc for about a cloud crossing time. Turbulence cascade and gravitation from the kpc scales are de facto included in smaller structures like molecular clouds. We show that the formation of a bar influences the dynamics of the central ~100 pc by creating resonances. At larger radii, the spiral arms host the formation of regularly spaced clouds: beads on a string and spurs. These instabilities pump turbulent energy into the gas, generally in the supersonic regime. Because of asymmetric drift, the supernovae explode outside of their gaseous nursery, which diminishes the effect of feedback on the structure of clouds. The evolution of clouds is thus mostly due to fragmentation and gas consumption, regulated mainly by supersonic turbulence. The transition from turbulence supported to self-gravitating gas is detected in the gas density probability distribution function at ~2000 cm^-3. The power spectrum density suggests that gravitation governs the hierarchical organisation of structures from the galactic scale down to a few parsecs.

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.

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.

In a recent study (Martizzi et al. 2012), we used cosmological simulations to show that active galactic nuclei (AGN) feedback on the gas distribution in clusters of galaxies can be important in determining the spatial distribution of stars and dark matter in the central regions of these systems. The hierarchical assembly of dark matter, baryons and black holes obscures the physical mechanism behind the restructuring process. Here we use idealized simulations to follow the response of a massive dark matter halo as we feed the central black hole with a controlled supply of cold gas. This removes most of the complexity taking place in the cosmological simulations that may have biased our previous study. We confirm our previous results: gas heated and expelled from the central regions of the halo by AGN feedback can return after cooling; repeated cycles generate gravitational potential fluctuations responsible for irreversible modifications of the dark matter mass profile. The main result is the expulsion of large amounts of baryons and dark matter from the central regions of the halo. According to the work presented here, outflow induced fluctuations represent the only mechanism able to efficiently create dark matter cores in clusters of galaxies.

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)

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.

We study the global efficiency of star formation in high resolution hydrodynamical simulations of gas discs embedded in isolated early-type and spiral galaxies. Despite using a universal local law to form stars in the simulations, we find that the early-type galaxies are offset from the spirals on the large-scale Kennicutt relation, and form stars 2 to 5 times less efficiently. This offset is in agreement with previous results on morphological quenching: gas discs are more stable against star formation when embedded in early-type galaxies due to the lower disc self-gravity and increased shear. As a result, these gas discs do not fragment into dense clumps and do not reach as high densities as in the spiral galaxies. Even if some molecular gas is present, the fraction of very dense gas (above 10^4 cm-3) is significantly reduced, which explains the overall lower star formation efficiency. We also analyse a sample of local early-type and spiral galaxies, measuring their CO and HI surface densities and their star formation rates as determined by their non-stellar 8um emission. As predicted by the simulations, we find that the early-type galaxies are offset from the Kennicutt relation compared to the spirals, with a twice lower efficiency. Finally, we validate our approach by performing a direct comparison between models and observations. We run a simulation designed to mimic the stellar and gaseous properties of NGC524, a lenticular galaxy, and find a gas disc structure and global star formation rate in good agreement with the observations. Morphological quenching thus seems to be a robust mechanism, and is also consistent with other observations of a reduced star formation efficiency in early-type galaxies in the COLD GASS survey. This lower efficiency of star formation is not enough to explain the formation of the whole Red Sequence, but can contribute to the reddening of some galaxies.

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.

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.

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, [{\alpha}/Fe]-[Fe/H] and Type Ia/II supernovae rates.

We compare the results of thirteen cosmological gasdynamical codes used to simulate the formation of a galaxy in the LCDM structure formation paradigm. The various runs differ in their hydrodynamical treatment (SPH, moving-mesh and AMR) but share the same initial conditions and adopt their latest published model of cooling, star formation and feedback. Despite the common halo assembly history, we find large code-to-code variations in the stellar mass, size, morphology and gas content of the galaxy at z=0, due mainly to the different implementations of feedback. Compared with observation, most codes tend to produce an overly massive galaxy, smaller and less gas-rich than typical spirals, with a massive bulge and a declining rotation curve. A stellar disk is discernible in most simulations, though its prominence varies widely from code to code. There is a well-defined trend between the effects of feedback and the severity of the disagreement with observation. Models that are more effective at limiting the baryonic mass of the galaxy come closer to matching observed galaxy scaling laws, but often to the detriment of the disk component. Our conclusions hold at two different numerical resolutions. Some differences can also be traced to the numerical techniques: more gas seems able to cool and become available for star formation in grid-based codes than in SPH. However, this effect is small compared to the variations induced by different feedback prescriptions. We conclude that state-of-the-art simulations cannot yet uniquely predict the properties of the baryonic component of a galaxy, even when the assembly history of its host halo is fully specified. Developing feedback algorithms that can effectively regulate the mass of a galaxy without hindering the formation of high-angular momentum stellar disks remains a challenge.

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.

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.

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.

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 \beta^{1/2} compared to the Bondi value, where \beta 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.

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.

We present for the first time predictions for UV line emission of intergalactic and circumgalactic gas from Adaptive Mesh Resolution (AMR) Large Scale Structure (LSS) simulations at redshifts 0.3<z<1.2, with specific emphasis on observability with current and near-future UV instrumentation. In three transitions of interest (Lya, OVI and CIV) there is a clear bimodality in the type of objects : the overwhelming majority of flux stems from discrete, compact sources, while a much larger volume fraction is filled by more tenuous gas. We characterise both object types with regard to number density, physical size and shape, brightness, luminosity, velocity structure, mass, temperature, ionisation state, and metal content. Degrading AMR grids to characteristic resolutions of available (such as FIREBall) or foreseeable instrumentation, allows to assess which inferences can be drawn from currently possible observations, and set foundations to prepare observing strategies for future missions. In general, the faint emission of the IGM and filamentary structure remains beyond capabilities of instruments with only short duration exposure potential (stratospheric balloons), even for optimistic assumptions for Lya, while the yet fainter metal line transitions for these structures will remain challenging for long duration exposures (space-based telescopes), mostly due to low metallicity pushing them more than three orders of magnitudes in brightness below Lya radiation. For the circum-galactic medium (CGM) the situation is more promising, and it is foreseeable that in the near future we will not only just dectect such sources, but the combination of all three lines in addition to velocity information will yield valuable insight into the physical processes at hand, illuminating important mechanisms during the formation of galaxies and their backreaction onto the IGM from whence they formed. (abrigded)

Euclid is a space-based survey mission from the European Space Agency designed to understand the origin of the Universe's accelerating expansion. It will use cosmological probes to investigate the nature of dark energy, dark matter and gravity by tracking their observational signatures on the geometry of the universe and on the cosmic history of structure formation. The mission is optimised for two independent primary cosmological probes: Weak gravitational Lensing (WL) and Baryonic Acoustic Oscillations (BAO). The Euclid payload consists of a 1.2 m Korsch telescope designed to provide a large field of view. It carries two instruments with a common field-of-view of ~0.54 deg2: the visual imager (VIS) and the near infrared instrument (NISP) which contains a slitless spectrometer and a three bands photometer. The Euclid wide survey will cover 15,000 deg2 of the extragalactic sky and is complemented by two 20 deg2 deep fields. For WL, Euclid measures the shapes of 30-40 resolved galaxies per arcmin2 in one broad visible R+I+Z band (550-920 nm). The photometric redshifts for these galaxies reach a precision of dz/(1+z) < 0.05. They are derived from three additional Euclid NIR bands (Y, J, H in the range 0.92-2.0 micron), complemented by ground based photometry in visible bands derived from public data or through engaged collaborations. The BAO are determined from a spectroscopic survey with a redshift accuracy dz/(1+z) =0.001. The slitless spectrometer, with spectral resolution ~250, predominantly detects Ha emission line galaxies. Euclid is a Medium Class mission of the ESA Cosmic Vision 2015-2025 programme, with a foreseen launch date in 2019. This report (also known as the Euclid Red Book) describes the outcome of the Phase A study.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Stars and more particularly massive stars, have a drastic impact on galaxy evolution. Yet the conditions in which they form and collapse are still not fully understood. In particular, the influence of the magnetic field on the collapse of massive clumps is relatively unexplored, it is thus of great relevance in the context of the formation of massive stars to investigate its impact. We perform high resolution, MHD simulations of the collapse of hundred solar masses, turbulent and magnetized clouds, using the adaptive mesh refinement code RAMSES. We compute various quantities such as mass distribution, magnetic field and angular momentum within the collapsing core and study the episodic outflows and the fragmentation that occurs during the collapse. The magnetic field has a drastic impact on the cloud evolution. We find that magnetic braking is able to substantially reduce the angular momentum in the inner part of the collapsing cloud. Fast and episodic outflows are being launched with typical velocities of the order of 3-5 km s$^{-1}$ although the highest velocities can be as high as 30-40 km s$^{-1}$. The fragmentation in several objects, is reduced in substantially magnetized clouds with respect to hydrodynamical ones by a factor of the order of 1.5-2. We conclude that magnetic fields have a significant impact on the evolution of massive clumps. In combination with radiation, magnetic fields largely determine the outcome of massive core collapse. We stress that numerical convergence of MHD collapse is a challenging issue. In particular, numerical diffusion appears to be important at high density therefore possibly leading to an over-estimation of the number of fragments.

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.

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.

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.

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.

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.

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.

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.

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}, $\rho_0$. This quantity is found to be always larger than the average density in a spherical shell of same radius $\bar{\rho}_0$, which is the quantity inferred from dynamical measurements in the Galaxy, and to vary in the range $\rho_0/\bar{\rho}_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.

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.

We study the properties of ISM substructure and turbulence in hydrodynamic (AMR) galaxy simulations with resolutions up to 0.8 pc and 5x10^3 Msun. We analyse the power spectrum of the density distribution, and various components of the velocity field. We show that the disk thickness is about the average Jeans scale length, and is mainly regulated by gravitational instabilities. From this scale of energy injection, a turbulence cascade towards small-scale is observed, with almost isotropic small-scale motions. On scales larger than the disk thickness, density waves are observed, but there is also a full range of substructures with chaotic and strongly non-isotropic gas velocity dispersions. The power spectrum of vorticity in an LMC-sized model suggests that an inverse cascade of turbulence might be present, although energy input over a wide range of scales in the coupled gaseous+stellar fluid could also explain this quasi-2D regime on scales larger than the disk scale height. Similar regimes of gas turbulence are also found in massive high-redshift disks with high gas fractions. Disk properties and ISM turbulence appear to be mainly regulated by gravitational processes, both on large scales and inside dense clouds. Star formation feedback is however essential to maintain the ISM in a steady state by balancing a systematic gas dissipation into dense and small clumps. Our galaxy simulations employ a thermal model based on a barotropic Equation of State (EoS) aimed at modelling the equilibrium of gas between various heating and cooling processes. Denser gas is typically colder in this approach, which is shown to correctly reproduce the density structures of a star-forming, turbulent, unstable and cloudy ISM down to scales of a few parsecs.

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.

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.

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.

We use high-resolution cosmological hydrodynamical AMR simulations to predict the characteristics of La emission from the cold gas streams that fed galaxies in massive haloes at high redshift. The La luminosity in our simulations is powered by the release of gravitational energy as gas flows from the intergalactic medium into the halo potential wells. The UV background contributes only <20% to the gas heating. The La emissivity is due primarily to electron-impact excitation cooling radiation in gas ~2x10^4K. We calculate the La emissivities assuming collisional ionisation equilibrium (CIE) at all gas temperatures. The simulated streams are self-shielded against the UV background, so photoionisation and recombination contribute negligibly to the La line formation. We produce theoretical maps of the La surface brightnesses, assuming that ~85% of the La photons are directly observable. We find that typical haloes of mass Mv~10^12-13 Msun at z~3 emit as La blobs (LABs) with luminosities 10^43-44 erg/s. Most of the La comes from the extended narrow, partly clumpy, inflowing, cold streams that feed the growing galaxies. The predicted LAB morphology is therefore irregular, with dense clumps and elongated extensions. The linewidth is expected to range from 10^2 to more than 10^3 km/s with a large variance. The typical La surface brightness profile is proportional to r^-1.2 where r is the distance from the halo centre. Our simulated LABs are similar in luminosity, morphology and extent to the observed LABs, with distinct kinematic features. The predicted La luminosity function is consistent with observations, and the predicted areas and linewidths roughly recover the observed scaling relations. This mechanism for producing LABs appears inevitable in many high-z galaxies. Some of the LABs may thus be regarded as direct detections of the cold streams that drove galaxy evolution at high z.

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.

We extract at redshift z=0 a Milky Way sized object including gas, stars and dark matter (DM) from a recent, high-resolution cosmological N-body simulation with baryons. Its resolution is sufficient to witness the formation of a rotating disk and bulge at the center of the halo potential. The phase-space structure of the central galactic halo reveals the presence of a dark disk component, that is co-rotating with the stellar disk. At the Earth's location, it contributes to around 25% of the total DM local density, whose value is rho_DM ~ 0.37 GeV/cm^3. The velocity distributions also show strong deviations from pure Gaussian and Maxwellian distributions, with a sharper drop of the high velocity tail. We give a detailed study of the impact of these features on the predictions for DM signals in direct detection experiments. In particular, the question of whether the modulation signal observed by DAMA is or is not excluded by limits set by other experiments (CDMS, XENON and CRESST...) is re-analyzed and compared to the case of a standard Maxwellian halo, in both the elastic and the inelastic scattering scenarios. We find that the compatibility between DAMA and the other experiments is improved. In the elastic scenario, the DAMA modulation signal is slightly enhanced in the so-called channeling region, as a result of several effects. For the inelastic scenario, the improvement of the fit is mainly attributable to the departure from a Maxwellian distribution at high velocity.

We study the dissipative effects of baryon physics on cosmic statistics at small scales using a cosmological simulation of a (50 Mpc/h)^3 volume of the universe. The MareNostrum simulation was performed using the AMR code RAMSES, and includes most of the physical ingredients which are part of the current theory of galaxy formation, such as metal-dependent cooling and UV heating, subgrid modelling of the ISM, star formation and supernova feedback. We re-ran the same initial conditions for a dark matter only universe, as a reference point for baryon-free cosmic statistics. In this paper, we present the measured small-scale amplification of sigma^2 and S_3 due to baryonic physics and their interpretation in the framework of the halo model. As shown in recent studies, the effect of baryons on the matter power spectrum can be accounted for at scales k <~ 10 h/Mpc by modifying the halo concentration parameter. We propose to extend this result by using a halo profile which is a linear combination of a NFW profile for the dark matter and an exponential disk profile mimicking the baryonic component at the heart of the halo. This halo profile form is physically motivated, and depends on two parameters, the mass fraction f_d of baryons in the disk, and the ratio lambda_d of the disk's characteristic scale to the halo's virial radius. We find this composite profile to reproduce both the small-scale variance and skewness boosts measured in the simulation up to k ~ 10^2 h/Mpc for physically meaningful values of the parameters f_d and lambda_d. Although simulations like the one presented here usually suffer from various problems when compared to observations, our modified halo model could be used as a fitting model to improve the determination of cosmological parameters from weak lensing convergence spectra and skewness measurements.

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.

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.

The correlation between the large-scale distribution of galaxies and their spectroscopic properties at z=1.5 is investigated using the Horizon MareNostrum cosmological run. We have extracted a large sample of 10^5 galaxies from this large hydrodynamical simulation featuring standard galaxy formation physics. Spectral synthesis is applied to these single stellar populations to generate spectra and colours for all galaxies. We use the skeleton as a tracer of the cosmic web and study how our galaxy catalogue depends on the distance to the skeleton. We show that galaxies closer to the skeleton tend to be redder, but that the effect is mostly due to the proximity of large haloes at the nodes of the skeleton, rather than the filaments themselves. This effects translate into a bimodality in the colour distribution of our sample. The origin of this bimodality is investigated and seems to follow from the ram pressure stripping of satellite galaxies within the more massive clusters of the simulation. The virtual catalogues (spectroscopical properties of the MareNostrum galaxies at various redshifts) are available online at http://www.iap.fr/users/pichon/MareNostrum/catalogues

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.

It is established that both radiative transfer and magnetic field have a strong impact on the collapse and the fragmentation of prestellar dense cores, but no consistent calculation exists yet at such scales. We present original AMR calculations including magnetic field (in the ideal MHD limit) and radiative transfer, within the Flux Limited Diffusion approximation, of the collapse of a 1 solar mass dense core. We compare the results with calculations performed with a barotropic EOS. We show that radiative transfer has an important impact on the collapse and the fragmentation, through the cooling or heating of the gas, and is complementary of the magnetic field. A larger field yields a stronger magnetic braking, increasing the accretion rate and thus the effect of the radiative feedback. Even for a strongly magnetized core, where the dynamics of the collapse is dominated by the magnetic field, radiative transfer is crucial to determine the temperature and optical depth distributions, two potentially accessible observational diagnostics. A barotropic EOS cannot account for realistic fragmentation. The diffusivity of the numerical scheme, however, is found to strongly affect the output of the collapse, leading eventually to spurious fragmentation. Both radiative transfer and magnetic field must be included in numerical calculations of star formation to obtain realistic collapse configurations and observable signatures. Nevertheless, the numerical resolution and the robustness of the solver are of prime importance to obtain reliable results. When using an accurate solver, the fragmentation is found to always remain inhibited by the magnetic field, at least in the ideal MHD limit, even when radiative transfer is included.

The largest galaxies acquire their mass early on, when the Universe is still youthful. Cold streams violently feed these young galaxies a vast amount of fresh gas, resulting in very efficient star formation. Using a well resolved hydrodynamical simulation of galaxy formation, we demonstrate that these mammoth galaxies are already in place a couple of billion years after the Big Bang. Contrary to local starforming galaxies, where dust re-emits a large part of the stellar ultraviolet (UV) light at infrared and sub-millimetre wavelengths, our self-consistent modelling of dust extinction predicts that a substantial fraction of UV photons should escape from primordial galaxies. Such a model allows us to compute reliably the number of high redshift objects as a function of luminosity, and yields galaxies whose UV luminosities closely match those measured in the deepest observational surveys available. This agreement is remarkably good considering our admittedly still simple modelling of the interstellar medium (ISM) physics. The luminosity functions (LF) of virtual UV luminous galaxies coincide with the existing data over the whole redshift range from 4 to 7, provided cosmological parameters are set to their currently favoured values. Despite their considerable emission at short wavelengths, we anticipate that the counterparts of the brightest UV galaxies will be detected by future sub-millimetre facilities like ALMA

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.

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.

Within the context of upcoming full-sky lensing surveys, the edge-preserving non- linear algorithm Aski is presented. Using the framework of Maximum A Posteriori inversion, it aims at recovering the full-sky convergence map from surveys with masks. It proceeds in two steps: CCD images of crowded galactic fields are deblurred using automated edge-preserving deconvolution; once the reduced shear is estimated, the convergence map is also inverted via an edge- preserving method. For the deblurring, it is found that when the observed field is crowded, this gain can be quite significant for realistic ground-based surveys when both positivity and edge-preserving penalties are imposed during the iterative deconvolution. For the convergence inversion, the quality of the reconstruction is investigated on noisy maps derived from the horizon N-body simulation, with and without Galactic cuts, and quantified using one-point statistics, power spectra, cluster counts, peak patches and the skeleton. It is found that the reconstruction is able to interpolate and extrapolate within the Galactic cuts/non-uniform noise; its sharpness-preserving penalization avoids strong biasing near the clusters of the map; it reconstructs well the shape of the PDF as traced by its skewness and kurtosis; the geometry and topology of the reconstructed map is close to the initial map as traced by the peak patch distribution and the skeleton's differential length; the two-points statistics of the recovered map is consistent with the corresponding smoothed version of the initial map; the distribution of point sources is also consistent with the corresponding smoothing, with a significant improvement when edge preserving prior is applied. The contamination of B-modes when realistic Galactic cuts are present is also investigated. Leakage mainly occurs on large scales.

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.

Weak gravitational lensing provides a unique method to map directly the dark matter in the Universe. The majority of lensing analyses uses the two-point statistics of the cosmic shear field to constrain the cosmological model yielding degeneracies, such as that between sigma_8 and Omega_M respectively the r.m.s. of the mass fluctuations at a scale of 8 Mpc/h and the matter density parameter both at z = 0. However, the two-point statistics only measure the Gaussian properties of the field and the weak lensing field is non-Gaussian. It has been shown that the estimation of non-Gaussian statistics on weak lensing data can improve the constraints on cosmological parameters. In this paper, we systematically compare a wide range of non-Gaussian estimators in order to determine which one provides tighter constraints on the cosmological parameters. These statistical methods include skewness, kurtosis and the Higher Criticism test in several sparse representations such as wavelet and curvelet; as well as the bispectrum, peak counting and a new introduced statistic called Wavelet Peak Counting (WPC). Comparisons based on sparse representations show that the wavelet transform is the most sensitive to non-Gaussian cosmological structures. It appears also that the best statistic for non-Gaussian characterization in weak lensing mass maps is the WPC. Finally, we show that the sigma_8 -Omega_m degeneracy could be even better broken if the WPC estimation is performed on weak lensing mass maps filtered by the wavelet method, MRLens.

With increasingly large data sets, weak lensing measurements are able to measure cosmological parameters with ever greater precision. However this increased accuracy also places greater demands on the statistical tools used to extract the available information. To date, the majority of lensing analyses use the two point-statistics of the cosmic shear field. These can either be studied directly using the two-point correlation function, or in Fourier space, using the power spectrum. But analyzing weak lensing data inevitably involves the masking out of regions or example to remove bright stars from the field. Masking out the stars is common practice but the gaps in the data need proper handling. In this paper, we show how an inpainting technique allows us to properly fill in these gaps with only $N \log N$ operations, leading to a new image from which we can compute straight forwardly and with a very good accuracy both the pow er spectrum and the bispectrum. We propose then a new method to compute the bispectrum with a polar FFT algorithm, which has the main advantage of avoiding any interpolation in the Fourier domain. Finally we propose a new method for dark matter mass map reconstruction from shear observations which integrates this new inpainting concept. A range of examples based on 3D N-body simulations illustrates the results.

We have performed a 70 billion dark-matter particles N-body simulation in a 2 $h^{-1}$ Gpc periodic box, using the concordance, cosmological model as favored by the latest WMAP3 results. We have computed a full-sky convergence map with a resolution of $\Delta \theta \simeq 0.74$ arcmin$^{2}$, spanning 4 orders of magnitude in angular dynamical range. Using various high-order statistics on a realistic cut sky, we have characterized the transition from the linear to the nonlinear regime at $\ell \simeq 1000$ and shown that realistic galactic masking affects high-order moments only below $\ell < 200$. Each domain (Gaussian and non-Gaussian) spans 2 decades in angular scale. This map is therefore an ideal tool for testing map-making algorithms on the sphere. As a first step in addressing the full map reconstruction problem, we have benchmarked in this paper two denoising methods: 1) Wiener filtering applied to the Spherical Harmonics decomposition of the map and 2) a new method, called MRLens, based on the modification of the Maximum Entropy Method on a Wavelet decomposition. While the latter is optimal on large spatial scales, where the signal is Gaussian, MRLens outperforms the Wiener method on small spatial scales, where the signal is highly non-Gaussian. The simulated full-sky convergence map is freely available to the community to help the development of new map-making algorithms dedicated to the next generation of weak-lensing surveys.

The massive galaxies in the young universe, ten billion years ago, formed stars at surprising intensities. Although this is commonly attributed to violent mergers, the properties of many of these galaxies are incompatible with such events, showing gas-rich, clumpy, extended rotating disks not dominated by spheroids (Genzel et al. 2006, 2008). Cosmological simulations and clustering theory are used to explore how these galaxies acquired their gas. Here we report that they are stream-fed galaxies, formed from steady, narrow, cold gas streams that penetrate the shock-heated media of massive dark matter haloes (Dekel & Birnboim 2006; Keres et al. 2005). A comparison with the observed abundance of star-forming galaxies implies that most of the input gas must rapidly convert to stars. One-third of the stream mass is in gas clumps leading to mergers of mass ratio greater than 1:10, and the rest is in smoother flows. With a merger duy cycle of 0.1, three-quarters of the galaxies forming stars at a given rate are fed by smooth streams. The rarer, submillimetre galaxies that form stars even more intensely are largely merger-induced starbursts. Unlike destructive mergers, the streams are likely to keep the rotating disk configuration intact, although turbulent and broken into giant star-forming clumps that merge into a central spheroid (Noguchi 1999; Genzel et al. 2008, Elmegreen, Bournaud & Elmegreen 2008, Dekel, Sari & Ceverino 2009). This stream-driven scenario for the formation of disks and spheroids is an alternative to the merger picture.

Observations of turbulent velocity dispersions in the HI component of galactic disks show a characteristic floor in galaxies with low star formation rates and within individual galaxies the dispersion profiles decline with radius. We carry out several high resolution adaptive mesh simulations of gaseous disks embedded within dark matter haloes to explore the roles of cooling, star-formation, feedback, shearing motions and baryon fraction in driving turbulent motions. In all simulations the disk slowly cools until gravitational and thermal instabilities give rise to a multi-phase medium in which a large population of dense self-gravitating cold clouds are embedded within a warm gaseous phase that forms through shock heating. The diffuse gas is highly turbulent and is an outcome of large scale driving of global non-axisymmetric modes as well as cloud-cloud tidal interactions and merging. At low star-formation rates these processes alone can explain the observed HI velocity dispersion profiles and the characteristic value of ~10 km/s observed within a wide range of disk galaxies. Supernovae feedback creates a significant hot gaseous phase and is an important driver of turbulence in galaxies with a star-formation rate per unit area >10^-3 M_sun/yr/kpc^2.

[Abridged]. We take advantage of the galaxy-like 3D dark matter map extracted from the HORIZON Project results to calculate the positron and antiproton fluxes from dark matter annihilation, in a model-independent approach as well as for dark matter particle benchmarks relevant at the LHC scale (from supersymmetric and extra-dimensional theories). Such a study is dedicated to a better estimate of the theoretical uncertainties affecting predictions, while the PAMELA and GLAST satellites are currently taking data which will soon provide better observational constraints. We discuss the predictions of the antiproton and positron fluxes, and of the positron fraction as well, as compared to the current data. We finally discuss the limits of the Nbody framework in describing the dark matter halo of our Galaxy.

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.

(abridged) The physics of diffuse gas accretion and the properties of the cold and hot modes of accretion onto proto-galaxies between z=2 and z=5.4 is investigated using the large cosmological simulation performed with the RAMSES code on the MareNostrum supercomputing facility. Galactic winds, chemical enrichment, UV background heating and radiative cooling are taken into account in this very high resolution simulation. Using accretion--weighted temperature histograms, we have perfomed two different measurements of the thermal state of the gas accreted towards the central galaxy. The first measurement, performed using accretion--weighted histograms on a spherical surface of radius 0.2 Rvir centred on the densest gas structure near the halo centre of mass, is a good indicator of the presence of an accretion shock in the vicinity of the galactic disc. We define the hot shock mass, Mshock, as the typical halo mass separating cold dominated from hot dominated accretion in the vicinity of the galaxy. The second measurement is performed by radially averaging histograms between 0.2 Rvir and Rvir, in order to detect radially extended structures such as gas filaments: this is a good proxy for detecting cold streams feeding the central galaxy. We define Mstream as the transition mass separating cold dominated from hot dominated accretion in the outer halo, marking the disappearance of these cold streams. We find a hot shock transition mass of Mshock=10^{11.6} Msun (dark matter), with no significant evolution with redshift. Conversely, we find that Mstream increases sharply with z. Our measurements are in agreement with the analytical predictions of Birnboim & Dekel (2003) and Dekel & Birnboim (2006), if we correct their model by assuming low metallicity (< 10^{-3} Zsun) for the filaments, correspondingly to our measurements.

In this paper, we estimate the gamma-ray and neutrino fluxes coming from dark matter annihilation in a Milky Way framework provided by a recent N-BODY HORIZON simulation. We first study the characteristics of the simulation and highlight the mass distribution within the galactic halo. The general dark matter density has a typical $r^{-3}$ power law for large radii, but the inner behaviour is poorly constrained below the resolution of the simulation ($\sim 200$ pc). We identify clumps and subclumps and analyze their distribution, as well as their internal structure. Inside the clumps, the power law is rather universal, $r^{-2.5}$ in the outer part with again strong uncertainties for smaller radii, especially for light clumps. We show a full-sky map of the astrophysical contribution to the gamma-ray or neutrino fluxes in this N-body framework. Using quite model independent and general assumptions for the high energy physics part, we evaluate the possible absolute fluxes and show some benchmark regions for the experiments GLAST, EGRET, and a km3 size extension of ANTARES like the KM3NeT project. While individual clumps seem to be beyond detection reach, the galactic center region is promising and GLAST could be sensitive to the geometry and the structure of its dark matter distribution. The detection by a km3 version of ANTARES is, however, more challenging due to a higher energy threshold. We also point out that the lack of resolution leaves the inner structure of subhalos poorly constrained. Using the same clump spectrum and mass fraction, a clump luminosity boost of order ten can be achieved with a steeper profile in the inner part of the sub-halos.

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.

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).

We present a detailed study of the X-ray properties of the galaxy cluster Abell 548b (z=0.04), using XMM-Newton data, and discuss the connection between the thermal properties and the presence of two extended relic radio sources located at the cluster periphery. We wish to analyze the dynamical state of the cluster and confirm the presence of a major merger. We will discuss the merger effects on the extended nonthermal emission. From the analysis of the temperature distribution and of the surface brightness profiles, we find evidence of a shock in the northern cluster region, just before the location of the two extended peripheral relics. From the optical analysis, we find that the cluster galaxies show a large $\sigma_V$, and two components are needed to fit their velocity distribution. Observational results were compared with a cluster simulation. The maps of gas temperature and density distribution from the simulation agree with the observational data in the case of a cluster merger nearly perpendicular to the plane of the sky and in the state after the maximum core collapse. We conclude that we are observing a galaxy cluster in a major merger phase, just after the maximum core collapse. The mass ratio is about 1:2, and the merger collision is nearly perpendicular to the plane of the sky. A shock is present in the northern cluster region, and it is very likely responsible for the electron reacceleration and the magnetic field amplification that will give cause the cluster relics. The relative position of the shock and the relics is strongly affected by projection effects.

This technical paper describes a software package that was designed to produce initial conditions for large cosmological simulations in the context of the Horizon collaboration. These tools generalize E. Bertschinger's Grafic1 software to distributed parallel architectures and offer a flexible alternative to the Grafic2 software for ``zoom'' initial conditions, at the price of large cumulated cpu and memory usage. The codes have been validated up to resolutions of 4096^3 and were used to generate the initial conditions of large hydrodynamical and dark matter simulations. They also provide means to generate constrained realisations for the purpose of generating initial conditions compatible with, e.g. the local group, or the SDSS catalog.

Various observations of magnetic fields in the Intra-Cluster Medium (ICM), most of the time restricted to cluster cores, point towards field strength of the order of a few microG (synchrotron radiation from radio relics and radio halos, inverse Compton radiation in X-rays and Faraday rotation measure of polarised background sources). Both the origin and the spatial structure of galaxy clusters magnetic fields are still under debate. In particular, the radial profile of the magnetic field, from the core of clusters to their outskirts, is of great importance for cosmic rays propagation within the Cosmic Web. In this letter, we highlight the importance of cooling processes in amplifying the magnetic field in the core of galaxy clusters up to one order of magnitude above the typical amplification obtained for a pure adiabatic evolution. We have performed a "zoom'' cosmological simulation of a 3 keV cluster, including dark matter and gas dynamics, atomic cooling, UV heating and star formation using the newly developed MHD solver in the AMR code RAMSES. Magnetic field amplification proceeds mainly through gravitational contraction. Shearing motions due to turbulence provide additional amplification in the outskirts of the cluster, while magnetic reconnection during mergers causes magnetic field dissipation in the core. Cooling processes have a strong impact on the magnetic field structure in the cluster. First, due to the sharp rise of the gas density in the centre, gravitational amplification is significantly amplified, when compared to the non--radiative run. Second, due to cooling processes, shearing motions are much stronger in the core than in the adiabatic case, leading to additional field amplification and no significant magnetic reconnection.

The development of parallel supercomputers allows today the detailed study of the collapse and the fragmentation of prestellar cores with increasingly accurate numerical simulations. Thanks to the advances in sub-millimeter observations, a wide range of observed initial conditions enable us to study the different modes of low-mass star formation. The challenge for the simulations is to reproduce the observational results. Two main numerical methods, namely AMR and SPH, are widely used to simulate the collapse and the fragmentation of prestellar cores. We compare thoroughly these two methods within their standard framework. We use the AMR code RAMSES and the SPH code DRAGON. Our physical model is as simple as possible and consists of an isothermal sphere rotating around the z-axis. We first study the conservation of angular momentum as a function of the resolution. Then, we explore a wide range of simulation parameters to study the fragmentation of prestellar cores. There seems to be a convergence between the two methods, provided resolution in each case is sufficient. Resolution criteria adapted to our physical cases, in terms of resolution per Jeans mass, for an accurate description of the formation of protostellar cores are deduced from the present study. This convergence is encouraging for future work in simulations of low-mass star formation, providing the aforementioned criteria are fulfilled. Higher resolution figures can be downloaded at http://www-dapnia.cea.fr/Projets/COAST/paper_amrvssph.pdf

To study the mass distribution of galaxy clusters up to their Virial radius, CL0016+16 seems to be a good candidate,since it is a bright massive cluster, previously considered as being dynamically relaxed. Using XMM-Newton observations of CL0016+16, we performed a careful X-ray background analysis, and we detected convincingly its X-ray emission up to $R_{200}$. We then studied its dynamical state with a detailed 2D temperature and surface brightness analysis of the inner part of the cluster. Using the assumption of both spherical symmetry and hydrostatic equilibrium (HE) we can determine the main cluster parameters: total mass, temperature profile, surface brightness profile and $\beta$-parameter. We also build a temperature map which clearly exhibits departure from spherical symmetry in the centre. To estimate the influence of these perturbations onto our total mass estimate, we also compute the total mass in the framework of the HE approach, but this time with various temperature profiles obtained in different directions. These various total mass estimates are consistent with each other. The temperature perturbations are clear signatures of ongoing merger activity. We also find significant residuals after subtracting the emissivity map by a 2D $\beta$-model fit. We conclude that, although CL0016+16 shows clear signs of merger activity and departure from spherical symmetry in the centre, its X-ray emissivity can be detected up to $R_{200}$ and the corresponding mass $M_{200}$ can be computed directly. It is therefore a good candidate to study cosmological scaling laws as predicted by the theory.

Abridged. A large fraction of stars are found in binary systems. It is therefore important for our understanding of the star formation process, to investigate the fragmentation of dense molecular cores. We study the influence of the magnetic field, ideally coupled to the gas, on the fragmentation in multiple systems of collapsing cores. We present high resolution numerical simulations performed with the RAMSES MHD code starting with a uniform sphere in solid body rotation and a uniform magnetic field parallel to the rotation axis. We pay particular attention to the strength of the magnetic field and interpret the results using the analysis presented in a companion paper. The results depend much on the amplitude, $A$, of the perturbations seeded initially. For a low amplitude, $A=0.1$, we find that for values of the mass-to-flux over critical mass-to-flux ratio, $\mu$, as high as $\mu = 20$, the centrifugally supported disk which fragments in the hydrodynamical case, is stabilized and remains axisymmetric. Detailed investigations reveals that this is due to the rapid growth of the toroidal magnetic field induced by the differential motions within the disk. For values of $\mu$ smaller $\simeq 5$, corresponding to larger magnetic intensities, there is no centrifugally supported disk because of magnetic braking. When the amplitude of the perturbation is equal to $A=0.5$, each initial peak develops independently and the core fragments for a large range of $\mu$. Only for values of $\mu$ close to 1 is the magnetic field able to prevent the fragmentation. Since a large fraction of stars are binaries, the results of low magnetic intensities preventing the fragmentation in case of weak perturbations, is problematic. We discuss three possible mechanisms...

The development of parallel supercomputers allows today the detailed study of the collapse and the fragmentation of prestellar cores with increasingly accurate numerical simulations. Thanks to the advances in sub-millimeter observations, a wide range of observed initial conditions enable us to study the different modes of low-mass star formation. The challenge for the simulations is to reproduce the observational results. Two main numerical methods, namely AMR and SPH, are widely used to simulate the collapse and the fragmentation of prestellar cores. We compare here thoroughly these two methods with numerical resolution requirements deduced from previous studies. Our physical model is as simple as possible, and consists of an isothermal sphere rotating around the z-axis. We first study the conservation of angular momentum as a function of the resolution. Then, we explore a wide range of simulation parameters to study the fragmentation of prestellar cores. There seems to be a convergence between the two methods, provided resolution in each case is sufficient. Resolution criteria adapted to our physical cases, in terms of resolution per Jeans mass, for an accurate description of the formation of protostellar cores are deduced from the present study. This convergence is encouraging for future work in simulations of low-mass star formation, providing the aforementioned criteria are fulfilled.

A radiative transfer scheme is presented, based on a moment description of the equation of radiative transfer and the so-called ``M1 closure model'' for the Eddington tensor. This model features a strictly hyperbolic transport step for radiation: it has been implemented using standard Godunov--like techniques in a new code called ATON. Coupled to simple models of ionization chemistry and photo-heating, ATON is able to reproduce the results of other schemes on a various set of standard tests such as the expansion of a HII region, the shielding of the radiation by dense clumps and cosmological ionization by multiple sources. Being simple yet robust, such a scheme is intended to be naturally and easily included in grid--based cosmological fluid solvers.

We studied the effect of supernovae feedback on a disk galaxy, taking into account the impact of infalling gas on both the star formation history and the corresponding outflow structure, the apparition of a supernovae-driven wind being highly sensitive to the halo mass, the galaxy spin and the star formation efficiency. We model our galaxies as cooling and collapsing NFW spheres. The dark matter component is modelled as a static external potential, while the baryon component is described by the Euler equations using the AMR code RAMSES. Metal-dependent cooling and supernovae-heating are also implemented using state-of-the-art recipes coming from cosmological simulations. We allow for 3 parameters to vary: the halo circular velocity, the spin parameter and the star formation efficiency. We found that the ram pressure of infalling material is the key factor limiting the apparition of galactic winds. We obtain a very low feedback efficiency, with supernovae to wind energy conversion factor around one percent, so that only low cicrular velocity galaxies give rise to strong winds. For massive galaxies, we obtain a galatic fountain, for which we discuss the observational properties. We conclude that for quiescent isolated galaxies, galactic winds appear only in very low mass systems. Although that can quite efficiently enrich the IGM with metals, they don't carry away enough cold material to solve the overcooling problem.

We use three-dimensional high-resolution adaptive-mesh-refinement simulations to investigate if mechanical feedback from active galactic nucleus jets can halt a massive cooling flow in a galaxy cluster and give rise to a self-regulated accretion cycle. We start with a 3 x 10^9MSun black hole at the centre of a spherical halo with the mass of the Virgo cluster. Initially, all the baryons are in a hot intracluster medium in hydrostatic equilibrium within the dark matter's gravitational potential. The black hole accretes the surrounding gas at the Bondi rate and a fraction of the accretion power is returned into the intracluster medium mechanically through the production of jets. The accretion, initially slow (~0.0002MSun/yr), becomes catastrophic, as the gas cools and condenses in the dark matter's potential. Therefore, it cannot prevent the cooling catastrophe at the centre of the cluster. However, after this rapid phase, where the accretion rate reaches a peak of ~0.2MSun/yr, the cavities inflated by the jets become highly turbulent. The turbulent mixing of the shock-heated gas with the rest of the intracluster medium puts a quick end to this short-lived rapid-growth phase. After dropping by almost two orders of magnitudes, the black hole accretion rate stabilises at ~0.006MSun/yr, without significant variations for several billions of years, indicating that a self-regulated steady-state has been reached. This accretion rate corresponds to a negligible increase of the black hole mass over the age of the Universe, but is sufficient to create a quasi-equilibrium state in the cluster core.

One way to access the aggregated power of a collection of heterogeneous machines is to use a grid middleware, such as DIET, GridSolve or NINF. It addresses the problem of monitoring the resources, of handling the submissions of jobs and as an example the inherent transfer of input and output data, in place of the user. In this paper we present how to run cosmological simulations using the RAMSES application along with the DIET middleware. We will describe how to write the corresponding DIET client and server. The remainder of the paper is organized as follows: Section 2 presents the DIET middleware. Section 3 describes the RAMSES cosmological software and simulations, and how to interface it with DIET. We show how to write a client and a server in Section 4. Finally, Section 5 presents the experiments realized on Grid'5000, the French Research Grid, and we conclude in Section 6.

Observations give evidences of the presence of metals in the intergalactic medium (IGM). The stars responsible for transforming hydrogen and helium into more complex atoms do not form outside the galaxies in the standard scenario of galaxy formation. Supernovae-driven winds and their associated feedback was proposed as a possible solution to explain such enrichment of the IGM. It turned out that a proper modelling of supernovae explosions within a turbulent interstellar medium (ISM) is a difficult task. Recent advances have been obtained using a multiphase approach to solve for the thermal state of the ISM, plus some additional recipes to account for the kinetic effect of supernovae on the galactic gas. We briefly describe here our implementation of supernovae feedback within the RAMSES code, and apply it to the formation and evolution of isolated galaxies of various masses and angular momenta. We have explored under what conditions a galactic wind can develop, if one considers only a quiescent mode of star formation. We have also characterized the distribution and evolution of metallicity in the gas outflow spreading in the IGM.

We have carried out a hydrodynamical code comparison study of interacting multiphase fluids. The two commonly used techniques of grid and smoothed particle hydrodynamics (SPH) show striking differences in their ability to model processes that are fundamentally important across many areas of astrophysics. Whilst Eulerian grid based methods are able to resolve and treat important dynamical instabilities, such as Kelvin-Helmholtz or Rayleigh-Taylor, these processes are poorly or not at all resolved by existing SPH techniques. We show that the reason for this is that SPH, at least in its standard implementation, introduces spurious pressure forces on particles in regions where there are steep density gradients. This results in a boundary gap of the size of the SPH smoothing kernel over which information is not transferred.

In this paper, we present a new method to perform numerical simulations of astrophysical MHD flows using the Adaptive Mesh Refinement framework and Constrained Transport. The algorithm is based on a previous work in which the MUSCL--Hancock scheme was used to evolve the induction equation. In this paper, we detail the extension of this scheme to the full MHD equations and discuss its properties. Through a series of test problems, we illustrate the performances of this new code using two different MHD Riemann solvers (Lax-Friedrich and Roe) and the need of the Adaptive Mesh Refinement capabilities in some cases. Finally, we show its versatility by applying it to two completely different astrophysical situations well studied in the past years: the growth of the magnetorotational instability in the shearing box and the collapse of magnetized cloud cores. We have implemented this new Godunov scheme to solve the ideal MHD equations in the AMR code RAMSES. It results in a powerful tool that can be applied to a great variety of astrophysical problems, ranging from galaxies formation in the early universe to high resolution studies of molecular cloud collapse in our galaxy.

We present the transverse flux correlation function of the Lyman-alpha forest in quasar absorption spectra at z~2.1 from VLT-FORS and VLT-UVES observations of a total of 32 pairs of quasars; 26 pairs with separations in the range 0.6<theta<4 arcmin and 6 pairs with 4<theta<10 arcmin. Correlation is detected at the 3sigma level up to separations of the order of ~4 arcmin (or ~4.4 h^{-1} Mpc comoving at z=2.1 for Omega_{m}=0.3 and Omega_{Lambda}=0.7). We have, furthermore, measured the longitudinal correlation function at a somewhat higher mean redshift (z=2.39) from 20 lines of sight observed with high spectral resolution and high signal-to-noise ratio with VLT-UVES. We compare the observed transverse and longitudinal correlation functions to that obtained from numerical simulations and illustrate the effect of spectral resolution, thermal broadening and peculiar motions. The shape and correlation length of the correlation functions are in good agreement with those expected from absorption by the filamentary and sheet-like structures in the photoionized warm intergalactic medium predicted in CDM-like models for structures formation. Using a sample of 139 C_{iv} systems detected along the lines of sight toward the pairs of quasars we also investigate the transverse correlation of metals on the same scales. The observed transverse correlation function of intervening C_{iv} absorption systems is consistent with that of a randomly distributed population of absorbers. This is likely due to the small number of pairs with separation less than 2 arcmin. We detect, however, a significant overdensity of systems in the sightlines towards the quartet Q0103-294A&B, Q0102-2931 and Q0102-293 which extends over the redshift range 1.5<z<2.2 and an angular scale larger than 10 arcmin.

The bulk of the extragalactic background between 10 keV and 10 GeV is likely to be explained by the emission of Seyfert galaxies, type Ia supernovae, and blazars. However, as revealed by the INTEGRAL satellite, the bulge of our galaxy is an intense source of a 511 keV gamma-ray line, indicating the production of a large number of positrons that annihilate. The origin of the latter is debated, and they could be produced, in particular, by the (S- or P-wave) annihilation of light Dark Matter particles into e+e-. In any case, the cumulated effect of similar sources at all redshifts could lead to a new background of hard X-ray and soft gamma-ray photons. On the basis of the hierarchical model of galaxy formation, we compute analytically the SNIa contribution to the background, and add it to Seyfert and blazars emission models. We find that any extra contribution to this unresolved background at 511 keV should be lower than about 4 keV/cm^2/s/sr. We also estimate analytically the extragalactic background due to Dark Matter annihilation, increasing the accuracy of the earlier computations. Indeed, we take into account the large positron escape fraction from low mass dark matter halos, unable to confine a dense and magnetized interstellar medium. Our new background estimate turns out to be one order of magnitude lower, so that the hypothesis of a light Dark Matter candidate remains compatible with the observed extragalactic background for a wider range of particle masses and cross-sections.

Archeops is a balloon--borne experiment widely inspired by the Planck satellite and by its High Frequency Instrument (HFI). It is mainly dedicated to measure the Cosmic Microwave Background (CMB) temperature anisotropies at high angular resolution (about 12 arcminutes) over a large fraction of the sky (around 30 %) in the millimetre and submillimetre range at 143, 217, 353 and 545 GHz. Further, the Archeops 353 GHz channel consists of three pairs of polarized sensitive bolometers designed to detect the polarized diffuse emission of Galactic dust. We present in this paper the update of the instrumental setup as well as the inflight performance for the last Archeops flight campaign in February 2002 from Kiruna (Sweden). We also describe the processing and analysis of the Archeops time ordered data for that campaign which lead to the measurement of the CMB anisotropies power spectrum in the multipole range l=10-700 (Benoit et al. 2003a, Tristram et al. 2005) and to the first measurement of the dust polarized emission at large angular scales and its polarized. We present maps of 30 % of the sky of the Galactic emission, including the Galactic plane, in the four Archeops channels at 143, 217, 353 and 545 GHz and maps of the CMB anisotropies at 143 and 217 GHz. These are the firstever available sub--degree resolution maps in the millimetre and submillimetre range of the large angular-scales Galactic dust diffuse emission and CMB temperature anisotropies respectively.

We propose to extend the well-known MUSCL-Hancock scheme for Euler equations to the induction equation modeling the magnetic field evolution in kinematic dynamo problems. The scheme is based on an integral form of the underlying conservation law which, in our formulation, results in a ``finite-surface'' scheme for the induction equation. This naturally leads to the well-known ``constrained transport'' method, with additional continuity requirement on the magnetic field representation. The second ingredient in the MUSCL scheme is the predictor step that ensures second order accuracy both in space and time. We explore specific constraints that the mathematical properties of the induction equations place on this predictor step, showing that three possible variants can be considered. We show that the most aggressive formulations (referred to as C-MUSCL and U-MUSCL) reach the same level of accuracy as the other one (referred to as Runge-Kutta), at a lower computational cost. More interestingly, these two schemes are compatible with the Adaptive Mesh Refinement (AMR) framework. It has been implemented in the AMR code RAMSES. It offers a novel and efficient implementation of a second order scheme for the induction equation. We have tested it by solving two kinematic dynamo problems in the low diffusion limit. The construction of this scheme for the induction equation constitutes a step towards solving the full MHD set of equations using an extension of our current methodology.

We derive the transverse flux correlation function in the Lyman-alpha forest at z~2.1 from VLT-FORS observations of a total of 32 pairs of quasars. The shape and correlation length of the transverse correlation function are in good agreement with the paradigm of intergalactic medium predicted in CDM-like models for structures formation. Using a sample of 139 Civ systems detected along the lines of sight toward the pairs of quasars we investigate the transverse correlation of metals on the same scales. We find that the correlation function is consistent with that of a randomly distributed population of Civ systems. However, we detect an important overdensity of systems in front of a quartet.

We have compared the cosmic microwave background (CMB) temperature anisotropy maps made from one-year time ordered data (TOD) streams that simulated observations of the originally planned 100 GHz Planck Low Frequency Instrument (LFI). The maps were made with three different codes. Two of these, ROMA and MapCUMBA, were implementations of maximum-likelihood (ML) map-making, whereas the third was an implementation of the destriping algorithm. The purpose of this paper is to compare these two methods, ML and destriping, in terms of the maps they produce and the angular power spectrum estimates derived from these maps. The difference in the maps produced by the two ML codes was found to be negligible. As expected, ML was found to produce maps with lower residual noise than destriping. In addition to residual noise, the maps also contain an error which is due to the effect of subpixel structure in the signal on the map-making method. This error is larger for ML than for destriping. If this error is not corrected a bias will be introduced in the power spectrum estimates. This study is related to Planck activities.

Using a series of high-resolution N-body hydrodynamical numerical simulations, we investigate several scenarios for the evolution of the baryon budget in galactic halos. We derive individual halo star formation history (SFH), as well as the global star formation rate in the universe. We develop a simple analytical model that allows us to compute surprisingly accurate predictions, when compared to our simulations, but also to other simulations presented in Springel & Hernquist (2003). The model depends on two main parameters: the star formation time scale t* and the wind efficiency eta_w. We also compute, for halos of a given mass, the baryon fraction in each of the following phases: cold disc gas, hot halo gas and stars. Here again, our analytical model predictions are in good agreement with simulation results, if one correctly takes into account finite resolution effect. We compare predictions of our analytical model to several observational constraints, and conclude that a very narrow range of the model parameters is allowed. The important role played by galactic winds is outlined, as well as a possible `superwind' scenario in groups and clusters. The `anti-hierarchical' behavior of observed SFH is well reproduced by our best model with t*=3Gyr and eta_w=1.5. We obtain in this case a present-day cosmic baryon budget of Omega*= 0.004, Omega_cold=0.0004, Omega_hot=0.01 and Omega_back=0.02 (diffuse background).

We present the first determination of the Galactic polarized emission at 353 GHz by Archeops. The data were taken during the Arctic night of February 7, 2002 after the balloon--borne instrument was launched by CNES from the Swedish Esrange base near Kiruna. In addition to the 143 GHz and 217 GHz frequency bands dedicated to CMB studies, Archeops had one 545 GHz and six 353 GHz bolometers mounted in three polarization sensitive pairs that were used for Galactic foreground studies. We present maps of the I, Q, U Stokes parameters over 17% of the sky and with a 13 arcmin resolution at 353 GHz (850 microns). They show a significant Galactic large scale polarized emission coherent on the longitude ranges [100, 120] and [180, 200] deg. with a degree of polarization at the level of 4-5%, in agreement with expectations from starlight polarization measurements. Some regions in the Galactic plane (Gem OB1, Cassiopeia) show an even stronger degree of polarization in the range 10-20%. Those findings provide strong evidence for a powerful grain alignment mechanism throughout the interstellar medium and a coherent magnetic field coplanar to the Galactic plane. This magnetic field pervades even some dense clouds. Extrapolated to high Galactic latitude, these results indicate that interstellar dust polarized emission is the major foreground for PLANCK-HFI CMB polarization measurement.

We analyze the cosmological constraints that Archeops places on adiabatic cold dark matter models with passive power-law initial fluctuations. Because its angular power spectrum has small bins in l and large l coverage down to COBE scales, Archeops provides a precise determination of the first acoustic peak in terms of position at multipole l_peak=220 +- 6, height and width. An analysis of Archeops data in combination with other CMB datasets constrains the baryon content of the Universe, Omega(b)h^2 = 0.022 (+0.003,-0.004), compatible with Big-Bang nucleosynthesis and with a similar accuracy. Using cosmological priors obtainedfrom recent non-CMB data leads to yet tighter constraints on the total density, e.g. Omega(tot)=1.00 (+0.03,-0.02) using the HST determination of the Hubble constant. An excellent absolute calibration consistency is found between Archeops and other CMB experiments, as well as with the previously quoted best fit model.The spectral index n is measured to be 1.04 (+0.10,-0.12) when the optical depth to reionization, tau, is allowed to vary as a free parameter, and 0.96 (+0.03,-0.04) when tau is fixed to zero, both in good agreement with inflation.

We present a determination by the Archeops experiment of the angular power spectrum of the cosmic microwave background anisotropy in 16 bins over the multipole range l=15-350. Archeops was conceived as a precursor of the Planck HFI instrument by using the same optical design and the same technology for the detectors and their cooling. Archeops is a balloon-borne instrument consisting of a 1.5 m aperture diameter telescope and an array of 21 photometers maintained at ~100 mK that are operating in 4 frequency bands centered at 143, 217, 353 and 545 GHz. The data were taken during the Arctic night of February 7, 2002 after the instrument was launched by CNES from Esrange base (Sweden). The entire data cover ~ 30% of the sky.This first analysis was obtained with a small subset of the dataset using the most sensitive photometer in each CMB band (143 and 217 GHz) and 12.6% of the sky at galactic latitudes above 30 degrees where the foreground contamination is measured to be negligible. The large sky coverage and medium resolution (better than 15 arcminutes) provide for the first time a high signal-to-noise ratio determination of the power spectrum over angular scales that include both the first acoustic peak and scales probed by COBE/DMR. With a binning of Delta(l)=7 to 25 the error bars are dominated by sample variance for l below 200. A companion paper details the cosmological implications.

We show the XMM-Newton guaranted time observations of 2 clusters of galaxies in a different stage of merger.

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.

We report here the first full sky component separation and CMB power spectrum estimation using a Wiener filtering technique on simulated data from the upcoming MAP experiment, set to launch in early 2001. The simulations included contributions from the three dominant astrophysical components expected in the five MAP spectral bands, namely CMB radiation, Galactic dust, and synchrotron emission. We assumed a simple homogeneous and isotropic white noise model and performed our analysis up to a spherical harmonic multipole lmax=512 on the fraction of the sky defined by b>20 degrees. We find that the reconstruction errors are reasonably well fitted by a Gaussian with an rms of 24 $\mu $K, but with significant deviations in the tails. Our results further support the predictions on the resulting CMB power spectrum of a previous estimate by Bouchet and Gispert (1999), which entailed a number of assumptions this work removes.

We present a self consistent method to perfom a joint analysis of Sunyaev-Zel'dovich and weak gravitational lensing observation of galaxy clusters. The spatial distribution of the cluster main constituents is described by a perturbative approach. Assuming the hydrostatic equilibrium and the equation of state, we are able to deduce, from observations, maps of projected gas density and gas temperature. The method then naturally entails a X-ray emissivity prediction which can be compared to observed X-ray emissivity maps. When tested on simulated clusters (noise free), this prediction turns out to be in very good agreement with the simulated surface brightness. The simulated and predicted surface brightness images have a correlation coefficient higher than 0.9 and the total flux differ by 0.9 % or 9 % in the two simulated clusters we studied. The method should be easily used on real data in order to provide a physical description of the cluster physics and of its constituents. The tests performed show that we can recover the amount and the spatial distributions of both the baryonic and non-baryonic material with an accuracy better than 10 %. So, in principle, in it might indeed help to alleviate some well known bias affecting, eg baryon fraction measurements.

The data analysis of current Cosmic Microwave Background (CMB) experiments like BOOMERanG or MAXIMA poses severe challenges which already stretch the limits of current (super-) computer capabilities, if brute force methods are used. In this paper we present a practical solution to the optimal map making problem which can be used directly for next generation CMB experiments like ARCHEOPS and TopHat, and can probably be extended relatively easily to the full PLANCK case. This solution is based on an iterative multi-grid Jacobi algorithm which is both fast and memory sparing. Indeed, if there are N_tod data points along the one dimensional timeline to analyse, the number of operations is O(N_tod ln N_tod) and the memory requirement is O(N_tod). Timing and accuracy issues have been analysed on simulated ARCHEOPS and TopHat data, and we discuss as well the issue of the joint evaluation of the signal and noise statistical properties.

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.

We investigate the influence of non equilibrium thermodynamics on cosmological structure formation. In this paper, we consider the collapse of planar perturbations usually called "Zel'dovich pancakes". We have developed for that purpose a new two fluids (gas and dark matter) hydrodynamical code, with three different thermodynamical species: electrons, ions and neutral particles (T_e\ne T_i \ne T_n). We describe in details the complex structure of accretion shock waves. We include several relevant processes for a low density, high temperature, collisional plasma such as non-equilibrium chemical reactions, cooling, shock heating, thermal energy equipartition between electrons, ions and neutral particles and electronic conduction. We find two different regions in the pancake structure: a thermal precursor ahead of the compression front and an equipartition wave after the compression front where electrons and ions temperatures differ significantly. This complex structure may have two interesting consequences: pre-heating of unshocked regions in the vicinity of massive X-ray clusters and ions and electrons temperatures differences in the outer regions of X-rays clusters.

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.

We investigate the role of radiative cooling within the core of large X-ray clusters using multi-phase hydrodynamics. We developed for that purpose a spherically symmetric hydrodynamical code, coupled to a "fluid model" that describes accurately the dark matter component. Cooling is included using a self-consistent multi-phase approach, leading to "cooled gas" mass deposition throughout the flow. We simulate the collapse and the subsequent evolution of a Coma-like X-ray cluster, avoiding the well-known "cooling catastrophe". The total mass profile of our simulated cluster is very similar to the "universal" profile proposed by Navarro, Frenk & White (1995). More interestingly, we also obtain a quasi-isothermal temperature profile, which is a direct consequence of multi-phase cooling within such a potential well.

We simulate the formation of a large X-ray cluster using a fully 3D hydrodynamical code coupled to a Particle-Mesh scheme which models the dark matter component. We focus on a possible decoupling between electrons and ions temperatures. We then solve the energy transfer equations between electrons, ions and neutrals without assuming thermal equilibrium between the three gases (T_e <> T_i <> T_n). We solve self-consistently the chemical equations for an hydrogen/helium primordial plasma without assuming ionization-recombination equilibrium. We find that the electron temperature differs from the true dynamical temperature by 20% at the Virial radius of our simulated cluster. This could lead marginally to an underestimate of the total mass in the outer regions of large X-ray clusters.

We investigate the formation of spherical cosmological structures following both dark matter and gas components. We focus on the dynamical aspect of the collapse assuming an adiabatic, $\gamma = 5/3$, fully ionized primordial plasma. We use for that purpose a fully Lagrangian hydrodynamical code designed to describe highly compressible flows in spherical geometry. We investigate also a "fluid approach" to describe the mean physical quantities of the dark matter flow. We test its validity for a wide range of initial density contrast. We show that an homogeneous isentropic core forms in the gas distribution, surrounded by a self-similar hydrostatic halo, with much higher entropy generated by shock dissipation. We derive analytical expressions for the size, density and temperature of the core, as well as for the surrounding halo. We show that, unless very efficient heating processes occur in the intergalactic medium, we are unable to reproduce within adiabatic models the typical core sizes in X-ray clusters. We also show that, for dynamical reasons only, the gas distribution is naturally antibiased relative to the total mass distribution, without invoking any reheating processes. This could explain why the gas fraction increases with radius in very large X-ray clusters. As a preparation for the next study devoted to the thermodynamical aspect of the collapse, we investigate the initial entropy level required to solve the core problem in X-ray clusters.

We reconsider the collapse of cosmic structures in an Einstein-de Sitter Universe, using the self similar initial conditions of Fillmore & Goldreich (1984). We first derive a new approximation to describe the dark matter dynamics in spherical geometry, that we refer to the "fluid approach". This method enables us to recover the self-similarity solutions of Fillmore & Goldreich for dark matter. We derive also new self-similarity solutions for the gas. We thus compare directly gas and dark matter dynamics, focusing on the differences due to their different dimensionalities in velocity space. This work may have interesting consequences for gas and dark matter distributions in large galaxy clusters, allowing to explain why the total mass profile is always steeper than the X-ray gas profile. We discuss also the shape of the dark matter density profile found in N-body simulations in terms of a change of dimensionality in the dark matter velocity space. The stable clustering hypothesis has been finally considered in the light of this analytical approach.

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 $\delta_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.