The recent detection of a 3.5 keV X-ray line from the centres of galaxies and clusters by Bulbul et al. (2014a) and Boyarsky et al. (2014a) has been interpreted as emission from the decay of 7 keV sterile neutrinos which could make up the (warm) dark
matter (WDM). As part of the COpernicus COmplexio (COCO) programme, we investigate the properties of dark matter haloes formed in a high-resolution cosmological $N$-body simulation from initial conditions similar to those expected in a universe in which the dark matter consists of 7 keV sterile neutrinos. This simulation and its cold dark matter (CDM) counterpart have $sim13.4$bn particles, each of mass $sim 10^5, h^{-1} M_odot$, providing detailed information about halo structure and evolution down to dwarf galaxy mass scales. Non-linear structure formation on small scales ($M_{200}, leq, 2 times 10^9,h^{-1},M_odot$) begins slightly later in COCO-Warm than in COCO-Cold. The halo mass function at the present day in the WDM model begins to drop below its CDM counterpart at a mass $sim 2 times 10^{9},h^{-1},M_odot$ and declines very rapidly towards lower masses so that there are five times fewer haloes of mass $M_{200}= 10^{8},h^{-1},M_odot$ in COCO-Warm than in COCO-Cold. Halo concentrations on dwarf galaxy scales are correspondingly smaller in COCO-Warm, and we provide a simple functional form that describes its evolution with redshift. The shapes of haloes are similar in the two cases, but the smallest haloes in COCO-Warm rotate slightly more slowly than their CDM counterparts.
Radiation damage to space-based Charge-Coupled Device (CCD) detectors creates defects which result in an increasing Charge Transfer Inefficiency (CTI) that causes spurious image trailing. Most of the trailing can be corrected during post-processing,
by modelling the charge trapping and moving electrons back to where they belong. However, such correction is not perfect -- and damage is continuing to accumulate in orbit. To aid future development, we quantify the limitations of current approaches, and determine where imperfect knowledge of model parameters most degrade measurements of photometry and morphology. As a concrete application, we simulate $1.5times10^{9}$ worst case galaxy and $1.5times10^{8}$ star images to test the performance of the Euclid visual instrument detectors. There are two separable challenges: If the model used to correct CTI is perfectly the same as that used to add CTI, $99.68$ % of spurious ellipticity is corrected in our setup. This is because readout noise is not subject to CTI, but gets over-corrected during correction. Second, if we assume the first issue to be solved, knowledge of the charge trap density within $Deltarho/rho!=!(0.0272pm0.0005)$ %, and the characteristic release time of the dominant species to be known within $Deltatau/tau!=!(0.0400pm0.0004)$ % will be required. This work presents the next level of definition of in-orbit CTI calibration procedures for Euclid.
The scale of Baryon Acoustic Oscillations (BAO) imprinted in the matter power spectrum provides an almost-perfect standard ruler: it only suffers sub-percent deviations from fixed comoving length due to non-linear effects. We study the BAO shift in t
he large Horndeski class of gravitational theories and compute its magnitude in momentum space using second-order perturbation theory and a peak-background split. The standard prediction is affected by the modified linear growth, as well as by non-linear gravitational effects that alter the mode-coupling kernel. For covariant Galileon models, we find a $14-45%$ enhancement of the BAO shift with respect to standard gravity and a distinct time evolution depending on the parameters. Despite the larger values, the shift remains well below the forecasted precision of next-generation galaxy surveys. Models that produce significant BAO shift would cause large redshift-space distortions or affect the bispectrum considerably. Our computation therefore validates the use of the BAO scale as a comoving standard ruler for tests of general dark energy models.
Using N-body simulations, we measure the power spectrum of the effective dark matter density field, which is defined through the modified Poisson equation in $f(R)$ cosmologies. We find that when compared to the conventional dark matter power spectru
m, the effective power spectrum deviates more significantly from the $Lambda$CDM model. For models with $f_{R0}=-10^{-4}$, the deviation can exceed 150% while the deviation of the conventional matter power spectrum is less than 50%. Even for models with $f_{R0}=-10^{-6}$, for which the conventional matter power spectrum is very close to the $Lambda$CDM prediction, the effective power spectrum shows sizeable deviations. Our results indicate that traditional analyses based on the dark matter density field may seriously underestimate the impact of $f(R)$ gravity on galaxy clustering. We therefore suggest the use of the effective density field in such studies. In addition, based on our findings, we also discuss several possible methods of making use of the differences between the conventional and effective dark matter power spectra in $f(R)$ gravity to discriminate the theory from the $Lambda$CDM model.
We show that Solar System tests can place very strong constraints on K-mouflage models of gravity, which are coupled scalar field models with nontrivial kinetic terms that screen the fifth force in regions of large gravitational acceleration. In part
icular, the bounds on the anomalous perihelion of the Moon imposes stringent restrictions on the K-mouflage Lagrangian density, which can be met when the contributions of higher-order operators in the static regime are sufficiently small. The bound on the rate of change of the gravitational strength in the Solar System constrains the coupling strength $beta$ to be smaller than $0.1$. These two bounds impose tighter constraints than the results from the Cassini satellite and Big Bang Nucleosynthesis. Despite the Solar System restrictions, we show that it is possible to construct viable models with interesting cosmological predictions. In particular, relative to $Lambda$-CDM, such models predict percent-level deviations for the clustering of matter and the number density of dark matter haloes. This makes these models predictive and testable by forthcoming observational missions.
We investigate the properties of dark matter haloes and subhaloes in an $f(R)$ gravity model with $|f_{R0}|=10^{-6}$, using a very high-resolution N-body simulation. The model is a borderline between being cosmologically interesting and yet still con
sistent with current data. We find that the halo mass function in this model has a maximum 20% enhancement compared with the $Lambda$CDM predictions between $z=1$ and $z=0$. Because of the chameleon mechanism which screens the deviation from standard gravity in dense environments, haloes more massive than $10^{13}h^{-1}M_odot$ in this $f(R)$ model have very similar properties to haloes of similar mass in $Lambda$CDM, while less massive haloes, such as that of the Milky Way, can have steeper inner density profiles and higher velocity dispersions due to their weaker screening. The halo concentration is remarkably enhanced for low-mass haloes in this model due to a deepening of the total gravitational potential. Contrary to the naive expectation, the halo formation time $z_f$ is later for low-mass haloes in this model, a consequence of these haloes growing faster than their counterparts in $Lambda$CDM at late times and the definition of $z_f$. Subhaloes, especially those less massive than $10^{11}h^{-1}M_odot$, are substantially more abundant in this $f(R)$ model for host haloes less massive than $10^{13}h^{-1}M_odot$. We discuss the implications of these results for the Milky Way satellite abundance problem. Although the overall halo and subhalo properties in this borderline $f(R)$ model are close to their $Lambda$CDM predictions, our results suggest that studies of the Local Group and astrophysical systems, aided by high-resolution simulations, can be valuable for further tests of it.
We introduce the idea of {it effective} dark matter halo catalog in $f(R)$ gravity, which is built using the {it effective} density field. Using a suite of high resolution N-body simulations, we find that the dynamical properties of halos, such as th
e distribution of density, velocity dispersion, specific angular momentum and spin, in the effective catalog of $f(R)$ gravity closely mimic those in the $Lambda$CDM model. Thus, when using effective halos, an $f(R)$ model can be viewed as a $Lambda$CDM model. This effective catalog therefore provides a convenient way for studying the baryonic physics, the galaxy halo occupation distribution and even semi-analytical galaxy formation in $f(R)$ cosmologies.
We present a comprehensive derivation of linear perturbation equations for different matter species, including photons, baryons, cold dark matter, scalar fields, massless and massive neutrinos, in the presence of a generic conformal coupling. Startin
g from the Lagrangians, we show how the conformal transformation affects the dynamics. In particular, we discuss how to incorporate consistently the scalar coupling in the equations of the Boltzmann hierarchy for massive neutrinos and the subsequent fluid approximations. We use the recently proposed K-mouflage model as an example to demonstrate the numerical implementation of our linear perturbation equations. K-mouflage is a new mechanism to suppress the fifth force between matter particles induced by the scalar coupling, but in the linear regime the fifth force is unsuppressed and can change the clustering of different matter species in different ways. We show how the CMB, lensing potential and matter power spectra are affected by the fifth force, and find ranges of K-mouflage parameters whose effects could be seen observationally. We also find that the scalar coupling can have the nontrivial effect of shifting the amplitude of the power spectra of the lensing potential and density fluctuations in opposite directions, although both probe the overall clustering of matter. This paper can serve as a reference for those who work on generic coupled scalar field cosmology, or those who are interested in the cosmological behaviour of the K-mouflage model.
We investigate void properties in $f(R)$ models using N-body simulations, focusing on their differences from General Relativity (GR) and their detectability. In the Hu-Sawicki $f(R)$ modified gravity (MG) models, the halo number density profiles of v
oids are not distinguishable from GR. In contrast, the same $f(R)$ voids are more empty of dark matter, and their profiles are steeper. This can in principle be observed by weak gravitational lensing of voids, for which the combination of a spectroscopic redshift and a lensing photometric redshift survey over the same sky is required. Neglecting the lensing shape noise, the $f(R)$ model parameter amplitudes $|f_{R0}|=10^{-5}$ and $10^{-4}$ may be distinguished from GR using the lensing tangential shear signal around voids by 4 and 8$sigma$ for a volume of 1~(Gpc/$h$)$^3$. The line-of-sight projection of large-scale structure is the main systematics that limits the significance of this signal for the near future wide angle and deep lensing surveys. For this reason, it is challenging to distinguish $|f_{R0}|=10^{-6}$ from GR. We expect that this can be overcome with larger volume. The halo void abundance being smaller and the steepening of dark matter void profiles in $f(R)$ models are unique features that can be combined to break the degeneracy between $|f_{R0}|$ and $sigma_8$.
We explore voids in dark matter and halo fields from simulations of $Lambda$CDM and Hu-Sawicki $f(R)$ models. In $f(R)$ gravity, dark matter void abundances are greater than that of general relativity (GR). However, when using haloes to identify void
s, the differences of void abundances become much smaller, but can still be told apart, in principle, at the 2, 6 and 14 $sigma$ level for the $f(R)$ model parameter amplitudes of $|f_{R0}|=10^{-6}$, $10^{-5}$ and $10^{-4}$. In contrast, the abundance of large voids found using haloes in $f(R)$ gravity is lower than in GR. The more efficient halo formation in underdense regions makes $f(R)$ voids less empty of haloes. This counter intuitive result suggests that voids are not necessarily emptier in $f(R)$ if one looks at galaxies in voids. Indeed, the halo number density profiles of voids are not distinguishable from GR. However, the same $f(R)$ voids are more empty of dark matter. This can in principle be observed by weak gravitational lensing of voids, for which the combination of a spec-$z$ and a photo-$z$ survey over the same sky is necessary. For a volume of 1~(Gpc/$h$)$^3$, neglecting the lensing shape noise, $|f_{R0}|=10^{-5}$ and $10^{-4}$ may be distinguished from GR using the lensing tangential shear signal around voids by 4 and 8$sigma$. The line-of-sight projection of large-scale structure is the main systematics that limits the significance of this signal, limiting the constraining power for $|f_{R0}|=10^{-6}$. The halo void abundance being smaller and the steepening of dark matter void profiles in $f(R)$ models are unique features that can be combined to break the degeneracy between $|f_{R0}|$ and $sigma_8$. The outflow of mass from void centers and velocity dispersions are greater in $f(R)$. Model differences in velocity profiles imply potential powerful constraints of the model in phase space and in redshift space.
