Initial Conditions for Large Cosmological Simulations

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.Comment: 12 pages, 11 figures, submitted to ApJ

Infall near clusters of galaxies: comparing gas and dark matter velocity profiles

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.Comment: 7 pages, 3 figures, to appear in MNRA

Collapse, outflows and fragmentation of massive, turbulent and magnetized prestellar barotropic cores

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.Comment: accepted for publication in A&

The effect of baryons on the variance and the skewness of the mass distribution in the Universe at small scales

We study the dissipative effects of baryon physics on cosmic statistics at small scales using a cosmological simulation of a (50 Mpc h−1)3 volume of universe. The MareNostrum simulation was performed using the adaptive mesh refinement (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 interstellar medium, star formation and supernova feedback. We reran 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 σ2 and S3 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−1 by modifying the halo concentration parameter. We propose to extend this result by using a composite halo profile, which is a linear combination of a Navarro, Frenk and White profile for the dark matter component and an exponential disc 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 disc and the ratio λd of the disc'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∼ 102 h Mpc−1 for physically meaningful values of the parameters f d and λ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 measurement

Antimatter cosmic rays from dark matter annihilation: First results from an N-body experiment

[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.Comment: 19 pages, 9 figures. Backgrounds included and additional comments and figures on the positron fraction. Accepted for publication in PR

Systematic uncertainties in the determination of the local dark matter density

A precise determination of the local dark matter density and an accurate control over the corresponding uncertainties are of paramount importance for Dark Matter (DM) searches. Using very recent high-resolution numerical simulations of a Milky Way like object, we study the systematic uncertainties that affect the determination of the local dark matter density based on dynamical measurements in the Galaxy. In particular, extracting from the simulation with baryons the orientation of the Galactic stellar disk with respect to the DM distribution, we study the DM density for an observer located at $\sim$8 kpc from the Galactic center {\it on the stellar disk}, $\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.Comment: 6 pages, 3 figures, matches published versio

The distant galaxy cluster CL0016+16: X-ray analysis up to R200R_{200}

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.Comment: 11 pages, 17 figures, Accepted for publication in A&

Optical spectroscopy and the nature of the insulating state of rare-earth nickelates

Using a combination of spectroscopic ellipsometry and DC transport measurements, we determine the temperature dependence of the optical conductivity of NdNiO$_3$ and SmNiO$_{3}$ films. The optical spectra show the appearance of a characteristic two-peak structure in the near-infrared when the material passes from the metal to the insulator phase. Dynamical mean-field theory calculations confirm this two-peak structure, and allow to identify these spectral changes and the associated changes in the electronic structure. We demonstrate that the insulating phase in these compounds and the associated characteristic two-peak structure are due to the combined effect of bond-disproportionation and Mott physics associated with half of the disproportionated sites. We also provide insights into the structure of excited states above the gap.Comment: 12 pages, 13 figure

Star Formation Near Photodissociation Regions: Detection of a Peculiar Protostar Near Ced 201

We present the detection and characterization of a peculiar low-mass protostar (IRAS 22129+7000) located ~0.4 pc from Ced 201 Photodissociation Region (PDR) and ~0.2 pc from the HH450 jet. The cold circumstellar envelope surrounding the object has been mapped through its 1.2 mm dust continuum emission with IRAM-30m/MAMBO. The deeply embedded protostar is clearly detected with Spitzer/MIPS (70 um), IRS (20-35 um) and IRAC (4.5, 5.8, and 8 um) but also in the K_s band (2.15 um). Given the large "near- and mid-IR excess" in its spectral energy distribution, but large submillimeter-to-bolometric luminosity ratio (~2%), IRAS 22129+7000 must be a transition Class 0/I source and/or a multiple stellar system. Targeted observations of several molecular lines from CO, 13CO, C18O, HCO+ and DCO+ have been obtained. The presence of a collimated molecular outflow mapped with the CSO telescope in the CO J=3-2 line suggests that the protostar/disk system is still accreting material from its natal envelope. Indeed, optically thick line profiles from high density tracers such as HCO+ J=1-0 show a red-shifted-absorption asymmetry reminiscent of inward motions. We construct a preliminary physical model of the circumstellar envelope (including radial density and temperature gradients, velocity field and turbulence) that reproduces the observed line profiles and estimates the ionization fraction. The presence of both mechanical and (non-ionizing) FUV-radiative input makes the region an interesting case to study triggered star formation