BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Wojciech A. Hellwing

Wojciech A. Hellwing contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.COastro-ph.GAhep-phgr-qcastro-ph.IMhep-thphysics.hist-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
32works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

32 published item(s)

preprint2022arXiv

COWS: A filament finder for Hessian cosmic web identifiers

The large scale galaxy and matter distribution is often described by means of the cosmic web made up of voids, sheets, filaments and knots. Many different recipes exist for identifying this cosmic web. Here we focus on a sub-class of cosmic web identifiers, based on the analysis of the Hessian matrix, and proposed a method, called COsmic Web Skeleton (COWS), of separating a set of filaments cells into an ensemble of individual discreet filaments. Specifically, a thinning algorithm is applied to velocity shear tensor based cosmic web (V-web) to identify the spine of the filaments. This results in a set of filaments with well defined end-point and length. It is confirmed that these sit at local density ridges and align with the appropriate direction defined by the underlying velocity field. The radial density profile of these curved cylindrical filaments, as well as the distribution of their lengths is also examined. The robustness of all results is checked against resolution and the V-web threshold. The code for the COWS method has been made publicly available.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Universality of the halo mass function in modified gravity cosmologies

We study the halo mass function (HMF) in modified gravity (MG) models using a set of large $N$-body simulations -- the ELEPHANT suite. We consider two popular beyond-general relativity scenarios: the Hu-Sawicki chameleon $f(R)$ model and the normal branch of the Dvali-Gabadadze-Porrati braneworld (nDGP). We show that in MG, analytic formulation based on the Press-Schechter framework offers a grossly inaccurate description of the HMF. We find, however, that once the HMF is expressed in terms of the dimensionless multiplicity function, it approximately assumes a redshift-independent universal character for all the models. Exploiting this property, we propose universal fits for the MG HMF in terms of their fractional departures from the $\mathrm{Λ\text{CDM}}$ case. We find two enclosed formulas, one for $f(R)$ and another for nDGP, that provide a reliable description of the HMF over the mass range covered by the simulations. These are accurate to a few percent with respect to the $N$-body data. We test the extrapolation potential of our fits against separate simulations with a different cosmological background and mass resolution, and find very good accuracy, within $\sim 10\%$. A particularly interesting finding from our analysis is a Gaussian-like shape of the HMF deviation that seems to appear universally across the whole $f(R)$ family, peaking at a mass variance scale characteristic for each $f(R)$ variant. We attribute this behavior to the specific physics of the environmentally-dependent chameleon screening models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The skewness of z=0.5 redshift-space galaxy distribution in Modified Gravity

We study the reduced skewness, $S_{3,g}\equiv\barξ_{3,g}/\barξ_{2,g}^2$ of galaxy distribution at $z=0.5$ in two families of modfied gravity models: the Hu-Sawicki $f(R)$-gravity and normal-branch of Dvali-Gabadadze-Porrati (nDGP) models. We use a set of mock galaxy catalogues specifally designed to match CMASS spectroscopic galaxy sample. For the first time we investigate the third reduced moment of such galaxy distributions both in the redshift space. Our analysis confirms that the signal previously indicated only for dark matter halo catalogues persists also in realistic mock galaxy samples. This result offers a possibility to extract a potential modified gravity signal in $S_3$ from spectroscopic galaxy data without a need for a very precise and self-consistent RSD models constructed for each and every modified gravity scenario separately. We show that the relative deviations from $Λ$CDM~ $S_{3,g}$ of various modified gravity models can vary from $7$ down to $\sim 2-3\%$ effects. Albeit, the effect looks small, we show that for considered models it can foster a $2-3σ$ falsification. Finally we argue that galaxy sample of a significantly higher number density should provide even stronger constraints by limiting shot-noise effects affecting the $S_{3,g}$ estimates at small comoving separations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Constraints on the identity of the dark matter from strong gravitational lenses

The cold dark matter (CDM) cosmological model unambigously predicts that a large number of haloes should survive as subhaloes when they are accreted into a larger halo. The CDM model would be ruled out if such substructures were shown not to exist. By contrast, if the dark matter consists of warm particles (WDM), then below a threshold mass that depends on the particle mass far fewer substructures would be present. Finding subhaloes below a certain mass would then rule out warm particle masses below some value. Strong gravitational lensing provides a clean method to measure the subhalo mass function through distortions in the structure of Einstein rings and giant arcs.Using mock lensing observations constructed from high-resolution N-body simulations, we show that measurements of approximately 100 strong lens systems with a detection limit of $10^7 h^{-1} M_{\odot}$ would clearly distinguish CDM from WDM in the case where this consists of 7 keV sterile neutrinos such as those that might be responsible for the 3.5 keV X-ray emission line recently detected in galaxies and clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Copernicus Complexio: a high-resolution view of the small-scale Universe

We introduce Copernicus Complexio (COCO), a high-resolution cosmological N-body simulation of structure formation in the $Λ{\rm CDM}{}$ model. COCO follows an approximately spherical region of radius $\sim 17.4h^{-1}\,{\rm Mpc}$ embedded in a much larger periodic cube that is followed at lower resolution. The high resolution volume has a particle mass of $1.135\times10^5h^{-1}{\rm M}_{\odot}$ (60 times higher than the Millennium-II simulation). COCO gives the dark matter halo mass function over eight orders of magnitude in halo mass; it forms $\sim 60$ haloes of galactic size, each resolved with about 10 million particles. We confirm the power-law character of the subhalo mass function, $\bar{N}(>μ)\proptoμ^{-s}$, down to a reduced subhalo mass $M_{sub}/M_{200}\equivμ=10^{-6}$, with a best-fit power-law index, $s=0.94$, for hosts of mass $\langle M_{200}\rangle=10^{12}h^{-1}{\rm M}_{\odot}$. The concentration-mass relation of COCO haloes deviates from a single power law for masses $M_{200}<\textrm{a few}\times 10^{8}h^{-1}{\rm M}_{\odot}$, where it flattens, in agreement with results by Sanchez-Conde et al. The host mass invariance of the reduced maximum circular velocity function of subhaloes, $ν\equiv V_{max}/V_{200}$, hinted at in previous simulations, is clearly demonstrated over five orders of magnitude in host mass. Similarly, we find that the average, normalised radial distribution of subhaloes is approximately universal (i.e. independent of subhalo mass), as previously suggested by the Aquarius simulations of individual haloes. Finally, we find that at fixed physical subhalo size, subhaloes in lower mass hosts typically have lower central densities than those in higher mass hosts.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The effect of baryons on redshift space distortions and cosmic density and velocity fields in the EAGLE simulation

We use the EAGLE galaxy formation simulation to study the effects of baryons on the power spectrum of the total matter and dark matter distributions and on the velocity fields of dark matter and galaxies. On scales $k{\stackrel{>}{{}_\sim}} 4{h\,{\rm Mpc}^{-1}}$ the effect of baryons on the amplitude of the total matter power spectrum is greater than $1\%$. The back-reaction of baryons affects the density field of the dark matter at the level of $\sim3\%$ on scales of $1\leq k/({h\,{\rm Mpc}^{-1}})\leq 5$. The dark matter velocity divergence power spectrum at $k{\stackrel{<}{{}_\sim}}0.5{h\,{\rm Mpc}^{-1}}$ is changed by less than $1\%$. The 2D redshift-space power spectrum is affected at the level of $\sim6\%$ at $|\vec{k}|{\stackrel{>}{{}_\sim}} 1{h\,{\rm Mpc}^{-1}}$ (for $μ>0.5$), but for $|\vec{k}|\leq 0.4{h\,{\rm Mpc}^{-1}}$ it differs by less than $1\%$. We report vanishingly small baryonic velocity bias for haloes: the peculiar velocities of haloes with $M_{200}>3\times10^{11}{{\rm M}_{\odot}}$ (hosting galaxies with $M_{*}>10^9{{\rm M}_{\odot}}$) are affected at the level of at most $1~$km/s, which is negligible for $1\%$-precision cosmology. We caution that since EAGLE overestimates cluster gas fractions it may also underestimate the impact of baryons, particularly for the total matter power spectrum. Nevertheless, our findings suggest that for theoretical modelling of redshift space distortions and galaxy velocity-based statistics, baryons and their back-reaction can be safely ignored at the current level of observational accuracy. However, we confirm that the modelling of the total matter power spectrum in weak lensing studies needs to include realistic galaxy formation physics in order to achieve the accuracy required in the precision cosmology era.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Mass-Concentration-Redshift Relation of Cold and Warm Dark Matter Halos

We use a suite of cosmological simulations to study the mass-concentration-redshift relation, $c({\rm M},z)$, of dark matter halos. Our simulations include standard $Λ$-cold dark matter (CDM) models, and additional runs with truncated power spectra, consistent with a thermal warm dark matter (WDM) scenario. We find that the mass profiles of CDM and WDM halos are self-similar and well approximated by the Einasto profile. The $c({\rm M},z)$ relation of CDM halos is monotonic: concentrations decrease with increasing virial mass at fixed redshift, and decrease with increasing redshift at fixed mass. The mass accretion histories (MAHs) of CDM halos are also scale-free, and can be used to infer concentrations directly. These results do not apply to WDM halos: their MAHs are not scale-free because of the characteristic scale imposed by the power-spectrum suppression. Further, the WDM $c({\rm M},z)$ relation is non-monotonic: concentrations peak at a mass scale dictated by the truncation scale, and decrease at higher and lower masses. We show that the assembly history of a halo can still be used to infer its concentration, provided that the total mass of its progenitors is considered (the &#34;collapsed mass history&#34;; CMH), rather than just that of its main ancestor. This exploits the scale-free nature of CMHs to derive a simple scaling that reproduces the mass-concentration-redshift relation of both CDM and WDM halos over a vast range of halo masses and redshifts. Our model therefore provides a robust account of the mass, redshift, cosmology and power spectrum dependence of dark matter halo concentrations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Exploring the liminality: properties of haloes and subhaloes in borderline $f(R)$ gravity

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 consistent with current data. We find that the halo mass function in this model has a maximum 20% enhancement compared with the $Λ$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 $Λ$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 $Λ$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 $Λ$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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Modified Gravity N-body Code Comparison Project

Self-consistent ${\it N}$-body simulations of modified gravity models are a key ingredient to obtain rigorous constraints on deviations from General Relativity using large-scale structure observations. This paper provides the first detailed comparison of the results of different ${\it N}$-body codes for the $f(R)$, DGP, and Symmetron models, starting from the same initial conditions. We find that the fractional deviation of the matter power spectrum from $Λ$CDM agrees to better than $1\%$ up to $k \sim 5-10~h/{\rm Mpc}$ between the different codes. These codes are thus able to meet the stringent accuracy requirements of upcoming observational surveys. All codes are also in good agreement in their results for the velocity divergence power spectrum, halo abundances and halo profiles. We also test the quasi-static limit, which is employed in most modified gravity ${\it N}$-body codes, for the Symmetron model for which the most significant non-static effects among the models considered are expected. We conclude that this limit is a very good approximation for all of the observables considered here.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Planes of satellite galaxies: when exceptions are the rule

The detection of planar structures within the satellite systems of both the Milky Way (MW) and Andromeda (M31) has been reported as being in stark contradiction to the predictions of the standard cosmological model ($Λ$CDM). Given the ambiguity in defining a planar configuration, it is unclear how to interpret the low incidence of the MW and M31 planes in $Λ$CDM. We investigate the prevalence of satellite planes around galactic mass haloes identified in high resolution cosmological simulations. We find that planar structures are very common, and that ~10% of $Λ$CDM haloes have even more prominent planes than those present in the Local Group. While ubiquitous, the planes of satellite galaxies show a large diversity in their properties. This precludes using one or two systems as small scale probes of cosmology, since a large sample of satellite systems is needed to obtain a good measure of the object-to-object variation. This very diversity has been misinterpreted as a discrepancy between the satellite planes observed in the Local Group and $Λ$CDM predictions. In fact, ~10% of $Λ$CDM galactic haloes have planes of satellites that are as infrequent as the MW and M31 planes. The look-elsewhere effect plays an important role in assessing the detection significance of satellite planes and accounting for it leads to overestimating the significance level by a factor of 30 and 100 for the MW and M31 systems, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The COpernicus COmplexio: Statistical Properties of Warm Dark Matter Haloes

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A clear and measurable signature of modified gravity in the galaxy velocity field

The velocity field of dark matter and galaxies reflects the continued action of gravity throughout cosmic history. We show that the low-order moments of the pairwise velocity distribution, $v_{12}$, are a powerful diagnostic of the laws of gravity on cosmological scales. In particular, the projected line-of-sight galaxy pairwise velocity dispersion, $σ_{12}(r)$, is very sensitive to the presence of modified gravity. Using a set of high-resolution N-body simulations we compute the pairwise velocity distribution and its projected line-of-sight dispersion for a class of modified gravity theories: the chameleon \fR gravity and Galileon gravity (cubic and quartic). The velocities of dark matter halos with a wide range of masses would exhibit deviations from General Relativity at the $(5-10)σ$ level. We examine strategies for detecting these deviations in galaxy redshift and peculiar velocity surveys. If detected, this signature would be a &#34;smoking gun&#34; for modified gravity.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

DM haloes in the fifth-force cosmology

We investigate how long-range scalar interactions affect the properties of dark matter haloes. For doing so we employ the ReBEL model which implements an additional interaction between dark matter particles. On the phenomenological level this is equivalent to a modification of gravity. We analyse the differences between five ReBEL models and $Λ$CDM using a series of high resolution cosmological simulations. Emphasis is placed on investigating how halo properties change in the presence of a fifth force. We report that the density profile of ReBEL haloes is well described by the NFW profile but with mean concentrations from $5\%$ to a few times higher than the standard $Λ$CDM value. We also find a slight increase of the halo spin for haloes more massive than $5\times10^{11}\M_{\odot}$, reflecting a higher rotational support of those haloes due to scalar forces. In addition, the dark matter haloes in our models are more spherical than their counterparts in $Λ$CDM. The ReBEL haloes are also more virialised, with a large difference from $Λ$CDM for strong fifth forces and a much smaller change for weak scalar interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dynamics of pairwise motions in the Cosmic Web

We present results of analysis of the dark matter (DM) pairwise velocity statistics in different Cosmic Web environments. We use the DM velocity and density field from the Millennium 2 simulation together with the NEXUS+ algorithm to segment the simulation volume into voxels uniquely identifying one of the four possible environments: nodes, filaments, walls or cosmic voids. We show that the PDFs of the mean infall velocities $v_{12}$ as well as its spatial dependence together with the perpendicular and parallel velocity dispersions bear a significant signal of the large-scale structure environment in which DM particle pairs are embedded. The pairwise flows are notably colder and have smaller mean magnitude in wall and voids, when compared to much denser environments of filaments and nodes. We discuss on our results, indicating that they are consistent with a simple theoretical predictions for pairwise motions as induced by gravitational instability mechanism. Our results indicate that the Cosmic Web elements are coherent dynamical entities rather than just temporal geometrical associations. In addition it should be possible to observationally test various Cosmic Web finding algorithms by segmenting available peculiar velocity data and studying resulting pairwise velocity statistics

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Halo model and halo properties in Galileon gravity cosmologies

We investigate the performance of semi-analytical modelling of large-scale structure in Galileon gravity cosmologies using results from N-body simulations. We focus on the Cubic and Quartic Galileon models that provide a reasonable fit to CMB, SNIa and BAO data. We demonstrate that the Sheth-Tormen mass function and linear halo bias can be calibrated to provide a very good fit to our simulation results. We also find that the halo concentration-mass relation is well fitted by a power law. The nonlinear matter power spectrum computed in the halo model approach is found to be inaccurate in the mildly nonlinear regime, but captures reasonably well the effects of the Vainshtein screening mechanism on small scales. In the Cubic model, the screening mechanism hides essentially all of the effects of the fifth force inside haloes. In the case of the Quartic model, the screening mechanism leaves behind residual modifications to gravity, which make the effective gravitational strength time-varying and smaller than the standard value. Compared to normal gravity, this causes a deficiency of massive haloes and leads to a weaker matter clustering on small scales. For both models, we show that there are realistic halo occupation distributions of Luminous Red Galaxies that can match both the observed large-scale clustering amplitude and the number density of these galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Milky Way mass constraints from the Galactic satellite gap

We use the distribution of maximum circular velocities, $V_{max}$, of satellites in the Milky Way (MW) to constrain the virial mass, $M_{200}$, of the Galactic halo under an assumed prior of a $Λ$CDM universe. This is done by analysing the subhalo populations of a large sample of halos found in the Millennium II cosmological simulation. The observation that the MW has at most three subhalos with $V_{max}\ge30 km/s$ requires a halo mass $M_{200}\le1.4\times10^{12} M_\odot$, while the existence of the Magellanic Clouds (assumed to have $V_{max}\ge60 km/s$) requires $M_{200}\ge1.0\times10^{12} M_\odot$. The first of these conditions is necessary to avoid the &#34;too-big-to-fail&#34; problem highlighted by Boylan-Kolchin et al., while the second stems from the observation that massive satellites like the Magellanic Clouds are rare. When combining both requirements, we find that the MW halo mass must lie in the range $0.25 \le M_{200}/(10^{12} M_\odot) \le 1.4$ at $90\%$ confidence. The gap in the abundance of Galactic satellites between $30 km/s\le V_{max} \le 60 km/s$ places our galaxy in the tail of the expected satellite distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nonlinear structure formation in Nonlocal Gravity

We study the nonlinear growth of structure in nonlocal gravity models with the aid of N-body simulation and the spherical collapse and halo models. We focus on a model in which the inverse-squared of the d&#39;Alembertian operator acts on the Ricci scalar in the action. For fixed cosmological parameters, this model differs from $Λ{\rm CDM}$ by having a lower late-time expansion rate and an enhanced and time-dependent gravitational strength ($\sim 6\%$ larger today). Compared to $Λ{\rm CDM}$ today, in the nonlocal model, massive haloes are slightly more abundant (by $\sim 10\%$ at $M \sim 10^{14} M_{\odot}/h$) and concentrated ($\approx 8\%$ enhancement over a range of mass scales), but their linear bias remains almost unchanged. We find that the Sheth-Tormen formalism describes the mass function and halo bias very well, with little need for recalibration of free parameters. The fitting of the halo concentrations is however essential to ensure the good performance of the halo model on small scales. For $k \gtrsim 1 h/{\rm Mpc}$, the amplitude of the nonlinear matter and velocity divergence power spectra exhibits a modest enhancement of $\sim 12\%$ to $15\%$, compared to $Λ{\rm CDM}$ today. This suggests that this model might only be distinguishable from $Λ{\rm CDM}$ by future observational missions. We point out that the absence of a screening mechanism may lead to tensions with Solar System tests due to local time variations of the gravitational strength, although this is subject to assumptions about the local time evolution of background averaged quantities.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nonlinear structure formation in the Cubic Galileon gravity model

We model the linear and nonlinear growth of large scale structure in the Cubic Galileon gravity model, by running a suite of N-body cosmological simulations using the {\tt ECOSMOG} code. Our simulations include the Vainshtein screening effect, which reconciles the Cubic Galileon model with local tests of gravity. In the linear regime, the amplitude of the matter power spectrum increases by $\sim 25%$ with respect to the standard $Λ$CDM model today. The modified expansion rate accounts for $\sim 20%$ of this enhancement, while the fifth force is responsible for only $\sim 5%$. This is because the effective unscreened gravitational strength deviates from standard gravity only at late times, even though it can be twice as large today. In the nonlinear regime ($k \gtrsim 0.1 h\rm{Mpc}^{-1}$), the fifth force leads to only a modest increase ($\lesssim 8%$) in the clustering power on all scales due to the very efficient operation of the Vainshtein mechanism. Such a strong effect is typically not seen in other models with the same screening mechanism. The screening also results in the fifth force increasing the number density of halos by less than 10%, on all mass scales. Our results show that the screening does not ruin the validity of linear theory on large scales which anticipates very strong constraints from galaxy clustering data. We also show that, whilst the model gives an excellent match to CMB data on small angular scales ($l \gtrsim 50$), the predicted integrated Sachs-Wolf effect is in tension with Planck/WMAP results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Subhalo statistics of galactic halos: beyond the resolution limit

We study the substructure population of Milky Way (MW)-mass halos in the $Λ$CDM cosmology using a novel procedure to extrapolate subhalo number statistics beyond the resolution limit of N-body simulations. The technique recovers the mean and the variance of the subhalo abundance, but not its spatial distribution. It extends the dynamic range over which precise statistical predictions can be made by the equivalent of performing a simulation with 50 times higher resolution, at no additional computational cost. We apply this technique to MW-mass halos, but it can easily be applied to halos of any mass. We find up to $20\%$ more substructures in MW-mass halos than found in previous studies. Our analysis lowers the mass of the MW halo required to accommodate the observation that the MW has only three satellites with a maximum circular velocity $V_{max}\ge30 km/s$ in the $Λ$CDM cosmology. The probability of having a subhalo population similar to that in the MW is $20\%$ for a virial mass, $M_{200}=1\times10^{12} M_\odot$ and practically zero for halos more massive than $M_{200}=2\times10^{12} M_\odot$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Testing the quasi-static approximation in $f(R)$ gravity simulations

Numerical simulations in modified gravity have commonly been performed under the quasi-static approximation -- that is, by neglecting the effect of time derivatives in the equation of motion of the scalar field that governs the fifth force in a given modified gravity theory. To test the validity of this approximation, we analyse the case of $f(R)$ gravity beyond this quasi-static limit, by considering effects, if any, these terms have in the matter and velocity divergence cosmic fields. To this end, we use the adaptive mesh refinement code ECOSMOG to study three variants ($|f_{R}|= 10^{-4}[$F4$], 10^{-5}[$F5$]$ and $10^{-6}[$F6$]$) of the Hu-Sawicki $f(R)$ gravity model, each of which refers to a different magnitude for the scalar field that generates the fifth force. We find that for F4 and F5, which show stronger deviations from standard gravity, a low-resolution simulation is enough to conclude that time derivatives make a negligible contribution to the matter distribution. The F6 model shows a larger deviation from the quasi-static approximation, but one that diminishes when re-simulated at higher resolution. We therefore come to the conclusion that the quasi-static approximation is valid for the most practical applications in $f(R)$ cosmologies.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Alpha, Betti and the Megaparsec Universe: on the Topology of the Cosmic Web

We study the topology of the Megaparsec Cosmic Web in terms of the scale-dependent Betti numbers, which formalize the topological information content of the cosmic mass distribution. While the Betti numbers do not fully quantify topology, they extend the information beyond conventional cosmological studies of topology in terms of genus and Euler characteristic. The richer information content of Betti numbers goes along the availability of fast algorithms to compute them. For continuous density fields, we determine the scale-dependence of Betti numbers by invoking the cosmologically familiar filtration of sublevel or superlevel sets defined by density thresholds. For the discrete galaxy distribution, however, the analysis is based on the alpha shapes of the particles. These simplicial complexes constitute an ordered sequence of nested subsets of the Delaunay tessellation, a filtration defined by the scale parameter, $α$. As they are homotopy equivalent to the sublevel sets of the distance field, they are an excellent tool for assessing the topological structure of a discrete point distribution. In order to develop an intuitive understanding for the behavior of Betti numbers as a function of $α$, and their relation to the morphological patterns in the Cosmic Web, we first study them within the context of simple heuristic Voronoi clustering models. Subsequently, we address the topology of structures emerging in the standard LCDM scenario and in cosmological scenarios with alternative dark energy content. The evolution and scale-dependence of the Betti numbers is shown to reflect the hierarchical evolution of the Cosmic Web and yields a promising measure of cosmological parameters. We also discuss the expected Betti numbers as a function of the density threshold for superlevel sets of a Gaussian random field.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Betti numbers of Gaussian fields

We present the relation between the genus in cosmology and the Betti numbers for excursion sets of three- and two-dimensional smooth Gaussian random fields, and numerically investigate the Betti numbers as a function of threshold level. Betti numbers are topological invariants of figures that can be used to distinguish topological spaces. In the case of the excursion sets of a three-dimensional field there are three possibly non-zero Betti numbers; $β_0$ is the number of connected regions, $β_1$ is the number of circular holes, and $β_2$ is the number of three-dimensional voids. Their sum with alternating signs is the genus of the surface of excursion regions. It is found that each Betti number has a dominant contribution to the genus in a specific threshold range. $β_0$ dominates the high-threshold part of the genus curve measuring the abundance of high density regions (clusters). $β_1$ dominates the genus near the median thresholds which measures the topology of negatively curved iso-density surfaces, and $β_2$ corresponds to the low-threshold part measuring the void abundance. We average the Betti number curves (the Betti numbers as a function of the threshold level) over many realizations of Gaussian fields and find that both the amplitude and shape of the Betti number curves depend on the slope of the power spectrum $n$ in such a way that their shape becomes broader and their amplitude drops less steeply than the genus as $n$ decreases. This behaviour contrasts with the fact that the shape of the genus curve is fixed for all Gaussian fields regardless of the power spectrum. Even though the Gaussian Betti number curves should be calculated for each given power spectrum, we propose to use the Betti numbers for better specification of the topology of large scale structures in the universe.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Hierarchical clustering and the BAO signature

In this contribution we present the preliminary results regarding the non-linear BAO signal in higher-order statistics of the cosmic density field. We use ensembles of N-body simulations to show that the non-linear evolution changes the amplitudes of the BAO signal, but has a negligible effect on the scale of the BAO feature. The latter observation accompanied by the fact that the BAO feature amplitude roughly doubles as one moves to higher orders, suggests that the higher-order correlation amplitudes can be used as probe of the BAO signal.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Hierarchical clustering in chameleon $f(R)$ gravity

We use a suite of high resolution state-of-the-art N-body Dark Matter simulations of chameleon $f(R)$ gravity to study the higher order volume averaged correlation functions $\bar{ξ_n}$ together with the hierarchical $n$-th order correlation amplitudes $S_n=\barξ_n/\barξ_2^{n-1}$ and density distribution functions (PDF). We show that under the non-linear modifications of gravity the hierarchical scaling of the reduced cumulants is preserved. This is however characterised by significant changes of both the $\bar{ξ_n}$ and $S_n$&#39;s values and their scale dependence with respect to General Relativity gravity (GR). In addition, we measure a significant increase of the non linear $σ_8$ parameter reaching $14, 5$ and 0.5% in excess of the GR value for the three flavours of our $f(R)$ models. We further note that the values of the reduced cumulants up to order $n=9$ are significantly increased in $f(R)$ gravity for all our models at small scales $R\simlt 30\hmpc$. In contrast the values of the hierarchical amplitudes, $S_n$&#39;s, are smaller in $f(R)$ indicating that the modified gravity density distribution functions are deviating from the GR case. Furthermore we find that the redshift evolution of relative deviations of the $f(R)$ hierarchical correlation amplitudes is fastest at high and moderate redshifts $1\leq z \leq4$. The growth of these deviations significantly slows down in the low redshift universe. We also compute the PDFs and show that for scales below $\sim 20\hmpc$ they are significantly shifted in $f(R)$ gravity towards the low densities. Finally we discuss the implications of our theoretical predictions for measurements of the hierarchical clustering in galaxy redshifts surveys, including the important problems of the galaxy biasing and redshifts space distortions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Professor &#34;Nonlinear&#34;: Obituary and memoir for Roman Juszkiewicz (1952-2012)

On 28th and 29th of January 2013 we held an international meeting in Zielona Góra (Poland) honouring the first anniversary of premature passing away^(1) of Professor Roman Juszkiewicz. We have celebrated an opening of a new seminar room at the University of Zielona Góra commemorated to the memory of Roman Juszkiewicz and we have shared our anecdotes and memories of this great scientist and friend. Here we want to present a limited and short memoir and obituary for Roman Juszkiewicz. 1: Roman Juszkiewicz passed away on 28th of January 2012

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Simulating the quartic Galileon gravity model on adaptively refined meshes

We develop a numerical algorithm to solve the high-order nonlinear derivative-coupling equation associated with the quartic Galileon model, and implement it in a modified version of the RAMSES N-body code to study the effect of the Galileon field on the large-scale matter clustering. The algorithm is tested for several matter field configurations with different symmetries, and works very well. This enables us to perform the first simulations for a quartic Galileon model which provides a good fit to the cosmic microwave background (CMB) anisotropy, supernovae and baryonic acoustic oscillations (BAO) data. Our result shows that the Vainshtein mechanism in this model is very efficient in suppressing the spatial variations of the scalar field. However, the time variation of the effective Newtonian constant caused by the curvature coupling of the Galileon field cannot be suppressed by the Vainshtein mechanism. This leads to a significant weakening of the strength of gravity in high-density regions at late times, and therefore a weaker matter clustering on small scales. We also find that without the Vainshtein mechanism the model would have behaved in a completely different way, which shows the crucial role played by nonlinearities in modified gravity theories and the importance of performing self-consistent N-body simulations for these theories.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Skewness as a probe of Baryon Acoustic Oscillations

In this study we show that the skewness S_3 of the cosmic density field contains a significant and potentially detectable and clean imprint of Baryonic Acoustic Oscillations. Although the BAO signal in the skewness has a lower amplitude than second order measures like the two-point correlation function and power spectrum, it has the advantage of a considerably lower sensitivity to systematic influences. Because it lacks a direct dependence on bias if this concerns simple linear bias, skewness will be considerably less beset by uncertainties due to galaxy bias. Also, it has a weaker sensitivity to redshift distortion effects. We use perturbation theory to evaluate the magnitude of the effect on the volume-average skewness, for various cosmological models. One important finding of our analysis is that the skewness BAO signal occurs at smaller scales than that in second order statistics. For an LCDM spectrum with WMAP7 normalization, the BAO feature has a maximum wiggle amplitude of ~3% and appears at a scale of ~82Mpc/h. We conclude that the detection of BAO wiggles in future extensive galaxy surveys via the skewness of the observed galaxy distribution may provide us with a useful, and potentially advantageous, measure of the nature of Dark Energy.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The nonlinear matter and velocity power spectra in f(R) gravity

We study the matter and velocity divergence power spectra in a f(R) gravity theory and their time evolution measured from several large-volume N-body simulations with varying box sizes and resolution. We find that accurate prediction of the matter power spectrum in f(R) gravity places stronger requirements on the simulation than is the case with LCDM, because of the nonlinear nature of the fifth force. Linear perturbation theory is shown to be a poor approximation for the f(R) models, except when the chameleon effect is very weak. We show that the relative differences from the fiducial LCDM model are much more pronounced in the nonlinear tail of the velocity divergence power spectrum than in the matter power spectrum, which suggests that future surveys which target the collection of peculiar velocity data will open new opportunities to constrain modified gravity theories. A close investigation of the time evolution of the power spectra shows that there is a pattern in the evolution history, which can be explained by the properties of the chameleon-type fifth force in f(R) gravity. Varying the model parameter |f_R0|, which quantifies the strength of the departure from standard gravity, mainly varies the epoch marking the onset of the fifth force, as a result of which the different f(R) models are in different stages of the same evolutionary path at any given time

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Probing Dark Energy with Alpha Shapes and Betti Numbers

We introduce a new descriptor of the weblike pattern in the distribution of galaxies and matter: the scale dependent Betti numbers which formalize the topological information content of the cosmic mass distribution. While the Betti numbers do not fully quantify topology, they extend the information beyond conventional cosmological studies of topology in terms of genus and Euler characteristic used in earlier analyses of cosmological models. The richer information content of Betti numbers goes along with the availability of fast algorithms to compute them. When measured as a function of scale they provide a &#34;Betti signature&#34; for a point distribution that is a sensitive yet robust discriminator of structure. The signature is highly effective in revealing differences in structure arising in different cosmological models, and is exploited towards distinguishing between different dark energy models and may likewise be used to trace primordial non-Gaussianities. In this study we demonstrate the potential of Betti numbers by studying their behaviour in simulations of cosmologies differing in the nature of their dark energy.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Boosting hierarchical structure formation with scalar-interacting dark matter

We investigate the effect of long-range scalar interactions in dark matter (DM) models of cosmic structure formation with a particular focus on the formation times of haloes. Utilising $N$-body simulations with $512^3$ DM particles we show that in our models dark matter haloes form substantially earlier: tracing objects up to redshift $z\sim6$ we find that the formation time, as characterised by the redshift $z_{1/2}$ at which the halo has assembled half of its final mass, is gradually shifted from $z_{1/2}\approx 1.83$ in the fiducial LCDM model to $z_{1/2}\approx 2.54$ in the most extreme self-interaction model. This is accompanied by a shift of the redshift that marks the transition between merger and steady accretion epochs from $z_{*}\approx 4.32$ in the \lcdm\ halos to $z_{*}\approx 6.39$ in our strongest interaction model. In other words, the scalar-interacting model employed in this work produces more structures at high redshifts, prolonging at the same time the steady accretion phases. These effects taken together can help the LCDM model to account for a high redshift reionisation as indicated by the WMAP data and can alleviate issues related to the survival of the thin-disk dominated galaxies at low redshifts.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Galactic halos in cosmology with long-range scalar DM interaction

Based on a set of cosmological N-body simulations we analyze properties of the dark matter haloes (DM) in a galaxy mass range ($10^{11} - 10^{13} h^{-1}M_{\odot}$) in modified $\lcdm$ cosmology with additional dynamically screened scalar interactions in DM sector. Our simulations show that scalar interactions support picture of the Island Universe. Rapid structure formation processes are shifted into higher redshifts resulting in a much smaller accretion and merging rates for galactic haloes at low redshifts. Finally, we present how this &#34;fifth&#34; force affects halo properties, like density profile, triaxiality, ellipticities and the spin parameter.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Hierarchy of N-point functions in the LCDM and ReBEL cosmologies

In this work we investigate higher order statistics for the $\lcdm$ and ReBEL scalar-interacting dark matter models by analyzing $180\hmpc$ dark matter N-body simulation ensembles. The N-point correlation functions and the related hierarchical amplitudes, such as skewness and kurtosis, are computed using the Count-In-Cells method. Our studies demonstrate that the hierarchical amplitudes $S_n$ of the scalar-interacting dark matter model significantly deviate from the values in the $\lcdm$ cosmology on scales comparable and smaller then the screening length $r_s$ of a given scalar-interacting model. The corresponding additional forces that enhance the total attractive force exerted on dark matter particles at galaxy scales lowers the values of the hierarchical amplitudes $S_n$. We conclude that hypothetical additional exotic interactions in the dark matter sector should leave detectable markers in the higher-order correlation statistics of the density field. We focussed in detail on the redshift evolution of the dark matter field&#39;s skewness and kurtosis. From this investigation we find that the deviations from the canonical $\lcdm$ model introduced by the presence of the ``fifth&#39;&#39; force attain a maximum value at redshifts $0.5<z<2$. We therefore conclude that moderate redshift data are better suited for setting observational constraints on the investigated ReBEL models.

0 citations0 reviewsNo code linkedOpen access