A radiative transfer scheme for cosmological reionization based on a local Eddington tensor

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.Comment: 14 pages, 13 figures, submitted to MNRA

Free-floating molecular clumps and gas mixing: hydrodynamic aftermaths of the intracluster−-interstellar medium interaction

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.Comment: 9 pages, 6 figures. Accepted for publication in MNRA

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&

Chameleon f(R) gravity on the Virgo cluster scale

Models of modified gravity offer promising alternatives to the concordance Λ cold dark matter (Λ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, which recovers standard general relativity in high-density regimes, while reproducing the desired late time acceleration at cosmological scales. We run a suite of high-resolution zoom simulations using the ecosmog code to examine the effect of f(R) gravity on the properties of a halo that is analogous to the Virgo cluster. We show that the velocity dispersion profiles can potentially discriminate between f(R) models and ΛCDM, and provide complementary analysis of lensing signal profiles to explore the possibility to further distinguish the different f(R) models. Our results confirm the techniques explored by Cabré etal. to quantify the effect of environment in the behaviour 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 centre 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 subhaloes 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) gravit

Evolution of the mass, size, and star formation rate in high-redshift merging galaxies MIRAGE - A new sample of simulations with detailed stellar feedback

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.Comment: 18 pages, 12 figures. Accepted in A&

Optical properties of BiTeBr and BiTeCl

We present a comparative study of the optical properties - reflectance, transmission and optical conductivity - and Raman spectra of two layered bismuth-tellurohalides BiTeBr and BiTeCl at 300 K and 5 K, for light polarized in the a-b planes. Despite different space groups, the optical properties of the two compounds are very similar. Both materials are doped semiconductors, with the absorption edge above the optical gap which is lower in BiTeBr (0.62 eV) than in BiTeCl (0.77 eV). The same Rashba splitting is observed in the two materials. A non-Drude free carrier contribution in the optical conductivity, as well as three Raman and two infrared phonon modes, are observed in each compound. There is a dramatic difference in the highest infrared phonon intensity for the two compounds, and a difference in the doping levels. Aspects of the strong electron-phonon interaction are identified. Several interband transitions are assigned, among them the low-lying absorption $\beta$ which has the same value 0.25 eV in both compounds, and is caused by the Rashba spin splitting of the conduction band. An additional weak transition is found in BiTeCl, caused by the lower crystal symmetry.Comment: Accepted 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