BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Raul E. Angulo

Raul E. Angulo appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.COastro-ph.GAastro-ph.IMApplicationsastro-ph.SR

Catalog footprint

What is connected

41works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

41 published item(s)

preprint2025arXiv

J-PLUS: Turning Off the Bright Stars

Photometric surveys require precise point spread function (PSF) characterization, as it varies across filters and is crucial for accurate photometry and low surface brightness (LSB) studies. However, the small PSF size provided by default pipelines suits only barely resolved objects, making it difficult to analyze regions near bright stars (rendering those regions unusable). These components are then combined to generate a final PSF for each exposure and filter, spanning 15 mag arcsec-2 in surface brightness and 4 arcmin in radius in the broad bands. In narrow-band filters, the J-PLUS PSF exhibits two rings, whereas in broad-band filters, only one ring is observed. Additionally, the position of the ring shifts with filter wavelength: as the filters become redder, the ring radius increases. We find that there is no significant variation in the extended PSF observed as a function of time (within 2.5h) or position in the field of view. The radial profile of NGC 4212 (which is close to a star) is also studied before/after PSF-subtraction. We developed a novel method to determine the central coordinates of saturated stars, and classify stars without using Gaia magnitudes. Additionally, mirror reflections are automatically detected and masked. Furthermore, in combining different stars and various components of the PSF, we avoided the use of a fixed radius by introducing a new method that does not depend on radial measurements. Accurate characterization of the extended PSF and its subtraction improves sky subtraction, increases the effective area of the survey by about 10%, and enables the study of extended large LSB features in wide area surveys like J-PLUS. Our pipeline is published as free software (GNU GPLv3) an can be customized to other surveys such as J-PAS, where its impact will be even greater due to its depth. This paper is fully reproducible and produced from Commit 4860c70.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Accelerating Large-Scale-Structure data analyses by emulating Boltzmann solvers and Lagrangian Perturbation Theory

The linear matter power spectrum is an essential ingredient in all theoretical models for interpreting large-scale-structure observables. Although Boltzmann codes such as CLASS or CAMB are very efficient at computing the linear spectrum, the analysis of data usually requires $10^4$-$10^6$ evaluations, which means this task can be the most computationally expensive aspect of data analysis. Here, we address this problem by building a neural network emulator that provides the linear theory (total and cold) matter power spectrum in about one millisecond with 0.2% (0.5%) accuracy over redshifts $z \le 3$ ($z \le 9$), and scales $10^{-4} \le k \, [h {\rm Mpc^{-1}}] < 50$. We train this emulator with more than 200,000 measurements, spanning a broad cosmological parameter space that includes massive neutrinos and dynamical dark energy. We show that the parameter range and accuracy of our emulator is enough to get unbiased cosmological constraints in the analysis of a Euclid-like weak lensing survey. Complementing this emulator, we train 15 other emulators for the cross-spectra of various linear fields in Eulerian space, as predicted by 2nd-order Lagrangian Perturbation theory, which can be used to accelerate perturbative bias descriptions of galaxy clustering. Our emulators are specially designed to be used in combination with emulators for the nonlinear matter power spectrum and for baryonic effects, all of which are publicly available at http://www.dipc.org/bacco.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

J-PLUS: A catalogue of globular cluster candidates around the M81/M82/NGC3077 triplet of galaxies

Globular clusters (GCs) are proxies of the formation assemblies of their host galaxies. However, few studies exist targeting GC systems of spiral galaxies up to several effective radii. Through 12-band Javalambre Photometric Local Universe Survey (J-PLUS) imaging, we study the point sources around the M81/M82/NGC3077 triplet in search of new GC candidates. We develop a tailored classification scheme to search for GC candidates based on their similarity to known GCs via a principal components analysis (PCA) projection. Our method accounts for missing data and photometric errors. We report 642 new GC candidates in a region of 3.5 deg$^2$ around the triplet, ranked according to their Gaia astrometric proper motions when available. We find tantalising evidence for an overdensity of GC candidate sources forming a bridge connecting M81 and M82. Finally, the spatial distribution of the GC candidates $(g-i)$ colours is consistent with halo/intra-cluster GCs, i.e. it gets bluer as they get further from the closest galaxy in the field. We further employ a regression-tree based model to estimate the metallicity distribution of the GC candidates based on their J-PLUS bands. The metallicity distribution of the sample candidates is broad and displays a bump towards the metal-rich end. Our list increases the population of GC candidates around the triplet by 3-fold, stresses the usefulness of multi-band surveys in finding these objects, and provides a testbed for further studies analysing their spatial distribution around nearby (spirals) galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Priors on Lagrangian bias parameters from galaxy formation modelling

We study the relations among the parameters of the hybrid Lagrangian bias expansion model, fitting biased auto and cross power spectra up to $k_{\rm max} = 0.7 \, h \, \mathrm{Mpc}^{-1}$. We consider $\sim 8000$ halo and galaxy samples, with different halo masses, redshifts, galaxy number densities, and varying the parameters of the galaxy formation model. Galaxy samples are obtained through state-of-the-art extended subhalo abundance matching techniques and include both stellar-mass and star-formation-rate selected galaxies. All of these synthetic galaxies samples are publicly available at https://bacco.dipc.org/galpk.html. We find that the hybrid Lagrangian bias model provides accurate fits to all of our halo and galaxy samples. The coevolution relations between galaxy bias parameters, although roughly compatible with those obtained for haloes, show systematic shifts and larger scatter. We explore possible sources of this difference in terms of dependence on halo occupation and assembly bias of each sample. The bias parameter relations displayed in this work can be used as a prior for future Bayesian analyses employing the hybrid Lagrangian bias expansion model.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

J-PLUS: Stellar Parameters, C, N, Mg, Ca and [α/Fe] Abundances for Two Million Stars from DR1

Context. The Javalambre Photometric Local Universe Survey (J-PLUS) has obtained precise photometry in twelve specially designed filters for large numbers of Galactic stars. Deriving their precise stellar atmospheric parameters and individual elemental abundances is crucial for studies of Galactic structure, and the assembly history and chemical evolution of our Galaxy. Aims. Our goal is to estimate not only stellar parameters (effective temperature, Teff, surface gravity, log g, and metallicity, [Fe/H]), but also [α/Fe] and four elemental abundances ([C/Fe], [N/Fe], [Mg/Fe], and [Ca/Fe]) using data from J-PLUS DR1. Methods. By combining recalibrated photometric data from J-PLUS DR1, Gaia DR2, and spectroscopic labels from LAMOST, we design and train a set of cost-sensitive neural networks, the CSNet, to learn the non-linear mapping from stellar colors to their labels. Results. We have achieved precisions of δTeff {\sim}55K, δlogg{\sim}0.15dex, and δ[Fe/H]{\sim}0.07dex, respectively, over a wide range of temperature, surface gravity, and metallicity. The uncertainties of the abundance estimates for [α/Fe] and the four individual elements are in the range 0.04-0.08 dex. We compare our parameter and abundance estimates with those from other spectroscopic catalogs such as APOGEE and GALAH, and find an overall good agreement. Conclusions. Our results demonstrate the potential of well-designed, high-quality photometric data for determinations of stellar parameters as well as individual elemental abundances. Applying the method to J-PLUS DR1, we have obtained the aforementioned parameters for about two million stars, providing an outstanding data set for chemo-dynamic analyses of the Milky Way. The catalog of the estimated parameters is publicly accessible.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Large-scale dark matter simulations

We review the field of collisionless numerical simulations for the large-scale structure of the Universe. We start by providing the main set of equations solved by these simulations and their connection with General Relativity. We then recap the relevant numerical approaches: discretization of the phase-space distribution (focusing on N-body but including alternatives, e.g., Lagrangian submanifold and Schrödinger-Poisson) and the respective techniques for their time evolution and force calculation (Direct summation, mesh techniques, and hierarchical tree methods). We pay attention to the creation of initial conditions and the connection with Lagrangian Perturbation Theory. We then discuss the possible alternatives in terms of the micro-physical properties of dark matter (e.g., neutralinos, warm dark matter, QCD axions, Bose-Einstein condensates, and primordial black holes), and extensions to account for multiple fluids (baryons and neutrinos), primordial non-Gaussianity and modified gravity. We continue by discussing challenges involved in achieving highly accurate predictions. A key aspect of cosmological simulations is the connection to cosmological observables, we discuss various techniques in this regard: structure finding, galaxy formation and baryonic modelling, the creation of emulators and light-cones, and the role of machine learning. We finalise with a recount of state-of-the-art large-scale simulations and conclude with an outlook for the next decade.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Modelling galaxy clustering in redshift space with a Lagrangian bias formalism and $N$-body simulations

Improving the theoretical description of galaxy clustering on small scales is an important challenge in cosmology, as it can considerably increase the scientific return of forthcoming galaxy surveys -- e.g. tightening the bounds on neutrino masses and deviations from general relativity. In this paper, we propose and test a new model for the clustering of galaxies that is able to accurately describe redshift-space distortions even down to small scales. This model corresponds to a second-order perturbative Lagrangian bias expansion which is advected to Eulerian space employing a displacement field extracted from $N$-body simulations. Eulerian coordinates are then transformed into redshift space by directly employing simulated velocity fields augmented with nuisance parameters capturing various possible satellite fractions and intra-halo small-scale velocities. We quantify the accuracy of our approach against samples of physically-motivated mock galaxies selected according to either Stellar Mass (SM) or Star Formation Rate (SFR) at multiple abundances and at $z=0$ and $1$. We find our model describes the monopole, quadrupole, and hexadecapole of the galaxy-power spectra down to scales of $k\approx 0.6 [h/$Mpc] within the accuracy of our simulations. This approach could pave the way to significantly increase the amount of cosmological information to be extracted from future galaxy surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Modelling the large scale structure of the Universe as a function of cosmology and baryonic physics

We present and test a framework that models the three-dimensional distribution of mass in the Universe as a function of cosmological and astrophysical parameters. Our approach combines two different techniques: a rescaling algorithm that modifies the cosmology of gravity-only N-body simulations, and a baryonification algorithm which mimics the effects of astrophysical processes induced by baryons, such as star formation and AGN feedback. We show how this approach can accurately reproduce the effects of baryons on the matter power spectrum of various state-of-the-art hydro-dynamical simulations (EAGLE, Illustris, Illustris-TNG, Horizon-AGN, and OWLS,Cosmo-OWLS and BAHAMAS), to percent level from very large down to small, highly nonlinear scales, k= 5 h/Mpc, and from z=0 up to z=2. We highlight that, thanks to the heavy optimisation of the algorithms, we can obtain these predictions for arbitrary baryonic models and cosmology (including massive neutrinos and dynamical dark energy models) with an almost negligible CPU cost. Therefore, this approach is efficient enough for cosmological data analyses. With these tools in hand we explore the degeneracies between cosmological and astrophysical parameters in the nonlinear mass power spectrum. Our findings suggest that after marginalising over baryonic physics, cosmological constraints inferred from weak gravitational lensing should be moderately degraded.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Simulating the Complexity of the Dark Matter Sheet I: Numerical Algorithms

At early times dark matter has a thermal velocity dispersion of unknown amplitude which, for warm dark matter models, can influence the formation of nonlinear structure on observable scales. We propose a new scheme to simulate cosmologies with a small-scale suppression of perturbations that combines two previous methods in a way that avoids the numerical artefacts which have so far prevented either from producing fully reliable results. At low densities and throughout most of the cosmological volume, we represent the dark matter phase-sheet directly using high-accuracy interpolation, thereby avoiding the artificial fragmentation which afflicts particle-based methods in this regime. Such phase-sheet methods are, however, unable to follow the rapidly increasing complexity of the denser regions of dark matter haloes, so for these we switch to an N-body scheme which uses the geodesic deviation equation to track phase-sheet properties local to each particle. In addition, we present a novel high-resolution force calculation scheme based on an oct-tree of cubic force resolution elements which is well suited to approximate the force-field of our combined sheet+particle distribution. Our hybrid simulation scheme enables the first reliable simulations of the internal structure of low-mass haloes in a warm dark matter cosmology.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The manifestation of secondary halo bias on the galaxy population from IllustrisTNG300

We use the improved IllustrisTNG300 magneto-hydrodynamical cosmological simulation to revisit the effect that secondary halo bias has on the clustering of the central galaxy population. With a side length of 205 $h^{-1}$Mpc and significant improvements on the sub-grid model with respect to the previous Illustris boxes, IllustrisTNG300 allows us to explore the dependencies of galaxy clustering over a large cosmological volume and wide halo-mass range. We show, at high statistical significance, that the halo assembly bias signal (i.e., the secondary dependence of halo bias on halo formation redshift) manifests itself on the clustering of the central galaxy population when this is split by stellar mass, colour, specific star formation rate, and surface density. A significant detection is also obtained for galaxy size: at fixed halo mass, larger central galaxies are more tightly clustered than smaller central galaxies in haloes of mass M$_{\rm vir} \lesssim 10^{12.5}$ $h^{-1}$M$_{\odot}$. This effect, however, seems to be uncorrelated with halo formation time, unlike the rest of the secondary dependencies analysed. We also explore the transmission of the halo spin bias signal, i.e., the secondary dependence of halo bias on halo spin. Although galaxy spin retains little information about the total spin of the halo, the correlation is enough to produce a significant galaxy spin bias signal. We discuss possible ways to probe the spin bias effects with observations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Characterising Strong Lensing Galaxy Clusters using the Millennium-XXL and MOKA simulations

In this paper we investigate the strong lensing statistics in galaxy clusters. We extract dark matter haloes from the Millennium-XXL simulation, compute their Einstein radius distribution, and find a very good agreement with Monte Carlo predictions produced with the MOKA code. The distribution of the Einstein radii is well described by a log-normal distribution, with a considerable fraction of the largest systems boosted by different projection effects. We discuss the importance of substructures and triaxiality in shaping the size of the critical lines for cluster size haloes. We then model and interpret the different deviations, accounting for the presence of a Bright Central Galaxy (BCG) and two different stellar mass density profiles. We present scaling relations between weak lensing quantities and the size of the Einstein radii. Finally we discuss how sensible is the distribution of the Einstein radii on the cosmological parameters Ω_M-σ_8 finding that cosmologies with higher Ω_M and σ_8 possess a large sample of strong lensing clusters. The Einstein radius distribution may help distinguish Planck13 and WMAP7 cosmology at 3σ.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Cosmological $N$-body simulations with suppressed variance

We present and test a method that dramatically reduces variance arising from the sparse sampling of wavemodes in cosmological simulations. The method uses two simulations which are fixed (the initial Fourier mode amplitudes are fixed to the ensemble average power spectrum) and paired (with initial modes exactly out of phase). We measure the power spectrum, monopole and quadrupole redshift-space correlation functions, halo mass function and reduced bispectrum at $z=1$. By these measures, predictions from a fixed pair can be as precise on non-linear scales as an average over 50 traditional simulations. The fixing procedure introduces a non-Gaussian correction to the initial conditions; we give an analytic argument showing why the simulations are still able to predict the mean properties of the Gaussian ensemble. We anticipate that the method will drive down the computational time requirements for accurate large-scale explorations of galaxy bias and clustering statistics, enabling more precise comparisons with theoretical models, and facilitating the use of numerical simulations in cosmological data interpretation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Subhalo abundance matching and assembly bias in the EAGLE simulation

Subhalo abundance matching (SHAM) is a widely-used method to connect galaxies with dark matter structures in numerical simulations. SHAM predictions agree remarkably well with observations, yet they still lack strong theoretical support. We examine the performance, implementation, and assumptions of SHAM using the EAGLE project simulations. We find that $V_{\rm relax}$, the highest value of the circular velocity attained by a subhalo while it satisfies a relaxation criterion, is the subhalo property that correlates most strongly with galaxy stellar mass ($M_{\rm star}$). Using this parameter in SHAM, we retrieve the real-space clustering of EAGLE to within our statistical uncertainties on scales greater than $2$ Mpc for galaxies with $8.77<\log_{10}(M_{\rm star}[M_\odot])<10.77$. Conversely, clustering is overestimated by $30\%$ on scales below $2$ Mpc for galaxies with $8.77<\log_{10}(M_{\rm star}[M_\odot])<9.77$ because SHAM slightly overpredicts the fraction of satellites in massive haloes compared to EAGLE. The agreement is even better in redshift-space, where the clustering is recovered to within our statistical uncertainties for all masses and separations. Additionally, we analyse the dependence of galaxy clustering on properties other than halo mass, i.e. the assembly bias. We demonstrate assembly bias alters the clustering in EAGLE by $20\%$ and that $V_{\rm relax}$ captures its effect to within $15\%$. We trace small differences in the clustering to the failure of SHAM as typically implemented, i.e. the $M_{\rm star}$ assigned to a subhalo does not depend on i) its host halo mass, ii) whether it is a central or a satellite. In EAGLE we find that these assumptions are not completely satisfied.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The galaxy correlation function as a constraint on galaxy formation physics

We introduce methods which allow observed galaxy clustering to be used together with observed luminosity or stellar mass functions to constrain the physics of galaxy formation. We show how the projected two-point correlation function of galaxies in a large semi-analytic simulation can be estimated to better than ~10% using only a very small subsample of the subhalo merger trees. This allows measured correlations to be used as constraints in a Monte Carlo Markov Chain exploration of the astrophysical and cosmological parameter space. An important part of our scheme is an analytic profile which captures the simulated satellite distribution extremely well out to several halo virial radii. This is essential to reproduce the correlation properties of the full simulation at intermediate separations. As a first application, we use low-redshift clustering and abundance measurements to constrain a recent version of the Munich semi-analytic model. The preferred values of most parameters are consistent with those found previously, with significantly improved constraints and somewhat shifted "best" values for parameters that primarily affect spatial distributions. Our methods allow multi-epoch data on galaxy clustering and abundance to be used as joint constraints on galaxy formation. This may lead to significant constraints on cosmological parameters even after marginalising over galaxy formation physics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An adaptively refined phase-space element method for cosmological simulations and collisionless dynamics

N-body simulations are essential for understanding the formation and evolution of structure in the Universe. However, the discrete nature of these simulations affects their accuracy when modelling collisionless systems. We introduce a new approach to simulate the gravitational evolution of cold collisionless fluids by solving the Vlasov-Poisson equations in terms of adaptively refineable "Lagrangian phase space elements". These geometrical elements are piecewise smooth maps between Lagrangian space and Eulerian phase space and approximate the continuum structure of the distribution function. They allow for dynamical adaptive splitting to accurately follow the evolution even in regions of very strong mixing. We discuss in detail various one-, two- and three-dimensional test problems to demonstrate the performance of our method. Its advantages compared to N-body algorithms are: i) explicit tracking of the fine-grained distribution function, ii) natural representation of caustics, iii) intrinsically smooth gravitational potential fields, thus iv) eliminating the need for any type of ad-hoc force softening. We show the potential of our method by simulating structure formation in a warm dark matter scenario. We discuss how spurious collisionality and large-scale discreteness noise of N-body methods are both strongly suppressed, which eliminates the artificial fragmentation of filaments. Therefore, we argue that our new approach improves on the N-body method when simulating self-gravitating cold and collisionless fluids, and is the first method that allows to explicitly follow the fine-grained evolution in six-dimensional phase space.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An exploration of galaxy-galaxy lensing and galaxy clustering in the Millennium-XXL simulation

The combination of galaxy-galaxy lensing and galaxy clustering data has the potential to simultaneously constrain both the cosmological and galaxy formation models. In this paper we perform a comprehensive exploration of these signals and their covariances through a combination of analytic and numerical approaches. First, we derive analytic expressions for the projected galaxy correlation function and stacked tangential shear profile and their respective covariances, which include Gaussian and discreteness noise terms. Secondly, we measure these quantities from mock galaxy catalogues obtained from the Millennium-XXL simulation and semi-analytic models of galaxy formation. We find that on large scales (R>10 Mpc/h), the galaxy bias is roughly linear and deterministic. On smaller scales (R<5 Mpc/h) the bias is a complicated function of scale and luminosity, determined by the different spatial distribution and abundance of satellite galaxies present when different magnitude cuts are applied, as well as by the mass dependence of the host haloes on magnitude. Our theoretical model for the covariances provides a reasonably good description of the measured ones on small and large scales. However, on intermediate scales (1<R<10 Mpc/h), the predicted errors are ~2-3 times smaller, suggesting that the inclusion of higher-order, non-Gaussian terms in the covariance will be required for further improvements. Importantly, both our theoretical and numerical methods show that the galaxy-galaxy lensing and clustering signals have a non-zero cross-covariance matrix with significant bin-to-bin correlations. Future surveys aiming to combine these probes must take this into account in order to obtain unbiased and realistic constraints.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Bayesian redshift-space distortions correction from galaxy redshift surveys

We present a Bayesian reconstruction method which maps a galaxy distribution from redshift-space to real-space inferring the distances of the individual galaxies. The method is based on sampling density fields assuming a lognormal prior with a likelihood given by the negative binomial distribution function modelling stochastic bias. We assume a deterministic bias given by a power law relating the dark matter density field to the expected halo or galaxy field. Coherent redshift-space distortions are corrected in a Gibbs-sampling procedure by moving the galaxies from redshift-space to real-space according to the peculiar motions derived from the recovered density field using linear theory with the option to include tidal field corrections from second order Lagrangian perturbation theory. The virialised distortions are corrected by sampling candidate real-space positions (being in the neighbourhood of the observations along the line of sight), which are compatible with the bulk flow corrected redshift-space position adding a random dispersion term in high density collapsed regions. The latter are defined according to the eigenvalues of the Hessian. This approach presents an alternative method to estimate the distances to galaxies using the three dimensional spatial information, and assuming isotropy. Hence the number of applications is very broad. In this work we show the potential of this method to constrain the growth rate up to $k$ ~ 0.3 $h$ Mpc$^{-1}$. Furthermore it could be useful to correct for photo-metric redshift errors, and to obtain improved BAO reconstructions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Forecasts on neutrino mass constraints from the redshift-space two-point correlation function

We provide constraints on the accuracy with which the neutrino mass fraction, $f_ν$, can be estimated when exploiting measurements of redshift-space distortions, describing in particular how the error on neutrino mass depends on three fundamental parameters of a characteristic galaxy redshift survey: density, halo bias and volume. In doing this, we make use of a series of dark matter halo catalogues extracted from the BASICC simulation. The mock data are analysed via a Markov Chain Monte Carlo likelihood analysis. We find a fitting function that well describes the dependence of the error on bias, density and volume, showing a decrease in the error as the bias and volume increase, and a decrease with density down to an almost constant value for high density values. This fitting formula allows us to produce forecasts on the precision achievable with future surveys on measurements of the neutrino mass fraction. For example, a Euclid-like spectroscopic survey should be able to measure the neutrino mass fraction with an accuracy of $δf_ν \approx 6.7\times10^{-4}$, using redshift-space clustering once all the other cosmological parameters are kept fixed to the $Λ$CDM case.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Properties of Cosmic Velocity Fields

Understanding the velocity field is very important for modern cosmology: it gives insights to structure formation in general, and also its properties are crucial ingredients in modelling redshift-space distortions and in interpreting measurements of the kinetic Sunyaev-Zeldovich effect. Unfortunately, characterising the velocity field in cosmological N-body simulations is inherently complicated by two facts: i) The velocity field becomes manifestly multi-valued after shell-crossing and has discontinuities at caustics. This is due to the collisionless nature of dark matter. ii) N-body simulations sample the velocity field only at a set of discrete locations, with poor resolution in low-density regions. In this paper, we discuss how the associated problems can be circumvented by using a phase-space interpolation technique. This method provides extremely accurate estimates of the cosmic velocity fields and its derivatives, which can be properly defined without the need of the arbitrary "coarse-graining" procedure commonly used. We explore in detail the configuration-space properties of the cosmic velocity field on very large scales and in the highly nonlinear regime. In particular, we characterise the divergence and curl of the velocity field, present their one-point statistics, analyse the Fourier-space properties and provide fitting formulae for the velocity divergence bias relative to the non-linear matter power spectrum. We furthermore contrast some of the interesting differences in the velocity fields of warm and cold dark matter models. We anticipate that the high-precision measurements carried out here will help to understand in detail the dynamics of dark matter and the structures it forms.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The universality of the virial halo mass function and models for non-universality of other halo definitions

The abundance of galaxy clusters can constrain both the geometry and growth of structure in our Universe. However, this probe could be significantly complicated by recent claims of nonuniversality -- non-trivial dependences with respect to the cosmological model and redshift. In this work we analyse the dependance of the mass function on the way haloes are identified and establish if this can cause departures from universality. In order to explore this dependance, we use a set of different N-body cosmological simulations (Le SBARBINE simulations), with the latest cosmological parameters from the Planck collaboration; this first suite of simulations is followed by a lower resolution set, carried out with different cosmological parameters. We identify dark matter haloes using a Spherical Overdensity algorithm with varying overdensity thresholds (virial, 2000rho_c, 1000rho_c, 500rho_c, 200rho_c and 200rho_b) at all redshifts. We notice that, when expressed in term of the rescaled variable nu, the mass functionfor virial haloes is a nearly universal as a function of redshift and cosmology, while this is clearly not the case for the other overdensities we considered. We provide fitting functions for the halo mass function parameters as a function of overdensity, that allow to predict, to within a few percent accuracy, the halo mass function for a wide range of halo definitions, redshifts and cosmological models. We then show how the departures from universality associated with other halo definitions can be derived by combining the universality of the virial definition with the expected shape of the density profile of halos.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Characterizing the Best Cosmic Telescopes with the Millennium Simulations

Certain configurations of massive structures projected along the line of sight maximize the number of detections of gravitationally lensed $z\sim10$ galaxies. We characterize such lines of sight with the étendue $σ_μ$, the area in the source plane magnified over some threshold $μ$. We use the Millennium I and Millennium XXL cosmological simulations to determine the frequency of high $σ_μ$ beams on the sky, their properties, and efficient selection criteria. We define the best beams as having $σ_{μ>3} >2000$ arcsec$^2$, for a $z\sim10$ source plane, and predict $477 \pm 21$ such beams on the sky. The total mass in the beam and $σ_{μ>3}$ are strongly correlated. After controlling for total mass, we find a significant residual correlation between $σ_{μ>3}$ and the number of cluster-scale halos ($>10^{14} M_\odot h^{-1}$) in the beam. Beams with $σ_{μ>3} >2000$ arcsec$^2$, which should be best at lensing $z\sim10$ galaxies, are ten times more likely to contain multiple cluster-scale halos than a single cluster-scale halo. Beams containing an Abell 1689-like massive cluster halo often have additional structures along the line of sight, including at least one additional cluster-scale ($M_{200}>10^{14}M_\odot h^{-1}$) halo 28% of the time. Selecting beams with multiple, massive structures will lead to enhanced detection of the most distant and intrinsically faint galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Cosmological constraints from the CFHTLenS shear measurements using a new, accurate and flexible way of predicting nonlinear mass clustering

We explore the cosmological constraints from cosmic shear using a new way of modelling the non-linear matter correlation functions. The new formalism extends the method of Angulo & White (2010), which manipulates outputs of $N$-body simulations to represent the three-dimensional nonlinear mass distribution in different cosmological scenarios. We show that predictions from our approach for shear two-point correlations at $1$ to $300$ arcmin separations are accurate at the $\sim10$\% level, even for extreme changes in cosmology. For moderate changes, with target cosmologies similar to that preferred by analyses of recent Planck data, the accuracy is close to $\sim5$\%. We combine this approach with a MonteCarlo Markov Chain sampler to explore constraints on a $Λ$CDM model from the shear correlation functions measured in the Canada-France Hawaii Telescope Lensing Survey (CFHTLenS). We obtain constraints on the parameter combination $σ_8 (Ω_m/0.27)^{0.6} = 0.801 \pm 0.028$. Combined with results from CMB data, we obtain marginalised constraints on $σ_8 = 0.81 \pm 0.01$ and $Ω_m = 0.29 \pm 0.01$. These results are fully compatible with previous analyses, which supports the validity of our approach. We discuss the advantages of our method and the potential it offers, including a path to incorporate in detail the effects of baryons, among others effects, in future high-precision cosmological analyses.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Matched filter optimization of kSZ measurements with a reconstructed cosmological flow field

We develop and test a new statistical method to measure the kinematic Sunyaev-Zel'dovich (kSZ) effect. A sample of independently detected clusters is combined with the cosmic flow field predicted from a galaxy redshift survey in order to derive a matched filter that optimally weights the kSZ signal for the sample as a whole given the noise involved in the problem. We apply this formalism to realistic mock microwave skies based on cosmological $N$-body simulations, and demonstrate its robustness and performance. In particular, we carefully assess the various sources of uncertainty, cosmic microwave background primary fluctuations, instrumental noise, uncertainties in the determination of the velocity field, and effects introduced by miscentring of clusters and by uncertainties of the mass-observable relation (normalization and scatter). We show that available data (\plk\ maps and the MaxBCG catalogue) should deliver a $7.7σ$ detection of the kSZ. A similar cluster catalogue with broader sky coverage should increase the detection significance to $\sim 13σ$. We point out that such measurements could be binned in order to study the properties of the cosmic gas and velocity fields, or combined into a single measurement to constrain cosmological parameters or deviations of the law of gravity from General Relativity.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Extending the halo mass resolution of $N$-body simulations

We present a scheme to extend the halo mass resolution of N-body simulations of the hierarchical clustering of dark matter. The method uses the density field of the simulation to predict the number of sub-resolution dark matter haloes expected in different regions. The technique requires as input the abundance of haloes of a given mass and their average clustering, as expressed through the linear and higher order bias factors. These quantities can be computed analytically or, more accurately, derived from a higher resolution simulation as done here. Our method can recover the abundance and clustering in real- and redshift-space of haloes with mass below $\sim 7.5 \times 10^{13}h^{-1}M_{\odot}$ at $z=0$ to better than 10%. We demonstrate the technique by applying it to an ensemble of 50 low resolution, large-volume $N$-body simulations to compute the correlation function and covariance matrix of luminous red galaxies (LRGs). The limited resolution of the original simulations results in them resolving just two thirds of the LRG population. We extend the resolution of the simulations by a factor of 30 in halo mass in order to recover all LRGs. With existing simulations it is possible to generate a halo catalogue equivalent to that which would be obtained from a $N$-body simulation using more than 20 trillion particles; a direct simulation of this size is likely to remain unachievable for many years. Using our method it is now feasible to build the large numbers of high-resolution large volume mock galaxy catalogues required to compute the covariance matrices necessary to analyse upcoming galaxy surveys designed to probe dark energy.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Galaxy formation in WMAP1 and WMAP7 cosmologies

Using the technique of Angulo & White (2010) we scale the Millennium and Millennium-II simulations of structure growth in a LCDM universe from the cosmological parameters with which they were carried out (based on first-year results from the Wilkinson Microwave Anisotropy Probe, WMAP1) to parameters consistent with the seven-year WMAP data (WMAP7). We implement semi-analytic galaxy formation modelling on both simulations in both cosmologies to investigate how the formation, evolution and clustering of galaxies are predicted to vary with cosmological parameters. The increased matter density Omega_m and decreased linear fluctuation amplitude sigma8 in WMAP7 have compensating effects, so that the abundance and clustering of dark halos are predicted to be very similar to those in WMAP1 for z <= 3. As a result, local galaxy properties can be reproduced equally well in the two cosmologies by slightly altering galaxy formation parameters. The evolution of the galaxy populations is then also similar. In WMAP7, structure forms slightly later. This shifts the peak in cosmic star formation rate to lower redshift, resulting in slightly bluer galaxies at z=0. Nevertheless, the model still predicts more passive low-mass galaxies than are observed. For rp< 1Mpc, the z=0 clustering of low-mass galaxies is weaker for WMAP7 than for WMAP1 and closer to that observed, but the two cosmologies give very similar results for more massive galaxies and on large scales. At z>1 galaxies are predicted to be more strongly clustered for WMAP7. Differences in galaxy properties, including, clustering, in these two cosmologies are rather small up to redshift 3. Given that there are still considerable residual uncertainties in galaxy formation models, it is very difficult to distinguish WMAP1 from WMAP7 through observations of galaxy properties or their evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Galaxy formation on the largest scales: The impact of astrophysics on the BAO peak

We investigate the effects of galaxy formation on the baryonic acoustic oscillations (BAO) peak by applying semi-analytic modelling techniques to the Millennium-XXL, a $3 \times 10^{11}$ particle N-body simulation of similar volume to the future EUCLID survey. Our approach explicitly incorporates the effects of tidal fields and stochasticity on halo formation, as well as the presence of velocity bias, spatially correlated merger histories, and the connection of all these with the observable and physical properties of galaxies. We measure significant deviations in the shape of the BAO peak from the expectations of a linear bias model built on top of the nonlinear dark matter distribution. We find that the galaxy correlation function shows an excess close to the maximum of the BAO peak ($r\sim110 Mpc/h$) and a deficit at $r\sim90 Mpc/h$. Depending on the redshift, selection criteria and number density of the galaxy samples, these bias distortions can be up to 5% in amplitude. They are, however, largely absorbed by marginalization over nuisance parameters in current analytical modelling of the BAO peak in configuration space, in particular into the parameter that controls the broadening due to nonlinear evolution. As a result, the galaxy formation effects detected here are unlikely to bias the high-precision measurements planned by the upcoming generation of wide-field galaxy surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

How closely do baryons follow dark matter on large scales?

We investigate the large-scale clustering and gravitational interaction of baryons and dark matter (DM) over cosmic time using a set of collisionless N-body simulations. Both components, baryons and DM, are evolved from distinct primordial density and velocity power spectra as predicted by early-universe physics. We first demonstrate that such two-component simulations require an unconventional match between force and mass resolution (i.e. force softening on at least the mean particle separation scale). Otherwise, the growth on any scale is not correctly recovered because of a spurious coupling between the two species at the smallest scales. With these simulations, we then demonstrate how the primordial differences in the clustering of baryons and DM are progressively diminished over time. In particular, we explicitly show how the BAO signature is damped in the spatial distribution of baryons and imprinted in that of DM. This is a rapid process, yet it is still not fully completed at low redshifts. On large scales, the overall shape of the correlation function of baryons and DM differs by 2% at z = 9 and by 0.2% at z = 0. The differences in the amplitude of the BAO peak are approximately a factor of 5 larger: 10% at z = 9 and 1% at z = 0. These discrepancies are, however, smaller than effects expected to be introduced by galaxy formation physics in both the shape of the power spectrum and in the BAO peak, and are thus unlikely to be detected given the precision of the next generation of galaxy surveys. Hence, our results validate the standard practice of modelling the observed galaxy distribution using predictions for the total mass clustering in the Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Noiseless Gravitational Lensing Simulations

The microphysical properties of the DM particle can, in principle, be constrained by the properties and abundance of substructures in DM halos, as measured through strong gravitational lensing. Unfortunately, there is a lack of accurate theoretical predictions for the lensing signal of substructures, mainly because of the discreteness noise inherent to N-body simulations. Here we present Recursive-TCM, a method that is able to provide lensing predictions with an arbitrarily low discreteness noise, without any free parameters or smoothing scale. This solution is based on a novel way of interpreting the results of N-body simulations, where particles simply trace the evolution and distortion of Lagrangian phase-space volume elements. We discuss the advantages of this method over the widely used cloud-in-cells and adaptive-kernel smoothing density estimators. Applying the new method to a cluster-sized DM halo simulated in warm and cold DM scenarios, we show how the expected differences in their substructure population translate into differences in the convergence and magnification maps. We anticipate that our method will provide the high-precision theoretical predictions required to interpret and fully exploit strong gravitational lensing observations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Statistical and systematic errors in redshift-space distortion measurements from large surveys

We investigate the impact of statistical and systematic errors on measurements of linear redshift-space distortions (RSD) in future cosmological surveys, analyzing large catalogues of dark-matter halos from the BASICC simulation. These allow us to estimate the dependence of errors on typical survey properties, as volume, galaxy density and mass (i.e. bias factor) of the adopted tracer. We find that measures of the specific growth rate β=f/b using the Hamilton/Kaiser harmonic expansion of the redshift-space correlation function ξ(r_p,π) on scales larger than 3/h Mpc are typically under-estimated by up to 10% for galaxy sized halos. This is significantly larger than the corresponding statistical errors, which amount to a few percent, indicating the importance of non-linear improvements to the Kaiser model to obtain accurate measurements of the growth rate. We compare the statistical errors to predictions obtained with the Fisher information matrix, based on the usual FKP prescription for the errors on the power spectrum. We show that this produces parameter errors fairly similar to the standard deviations from the halo catalogues, but only if applied to strictly linear scales in Fourier space (k<0.2 h/Mpc). Finally, we present an accurate scaling formula describing the relative error on β as a function of the survey parameters, which closely matches the simulation results in all explored regimes. This provides a handy and plausibly more realistic alternative to the Fisher matrix approach, to quickly and accurately predict RSD statistical errors expected from future surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Dynamical State and Mass-Concentration Relation of Galaxy Clusters

We use the Millennium Simulation series to study how the dynamical state of dark matter halos affects the relation between mass and concentration. We find that a large fraction of massive systems are identified when they are substantially out of equilibrium and in a particular phase of their dynamical evolution: the more massive the halo, the more likely it is found at a transient stage of high concentration. This state reflects the recent assembly of massive halos and corresponds to the first pericentric passage of recently-accreted material when, before virialization, the kinetic and potential energies reach maximum and minimum values, respectively. This result explains the puzzling upturn in the mass-concentration relation reported in recent work for massive halos; indeed, the upturn disappears when only dynamically-relaxed systems are considered in the analysis. Our results warn against applying simple equilibrium models to describe the structure of rare, massive galaxy clusters and urges caution when extrapolating scaling laws calibrated on lower-mass systems, where such deviations from equilibrium are less common. The evolving dynamical state of galaxy clusters ought to be carefully taken into account if cluster studies are to provide precise cosmological constraints.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Mass Profile and Accretion History of Cold Dark Matter Halos

We use the Millennium Simulation series to study the relation between the accretion history (MAH) and mass profile of cold dark matter halos. We find that the mean density within the scale radius, r_{-2} (where the halo density profile has isothermal slope), is directly proportional to the critical density of the Universe at the time when the main progenitor's virial mass equals the mass enclosed within r_{-2}. Scaled to these characteristic values of mass and density, the mean MAH, expressed in terms of the critical density of the Universe, M(ρ_{crit}(z)), resembles that of the enclosed density profile, M(<ρ>), at z=0. Both follow closely the NFW profile, suggesting that the similarity of halo mass profiles originates from the mass-independence of halo MAHs. Support for this interpretation is provided by outlier halos whose accretion histories deviate from the NFW shape; their mass profiles show correlated deviations from NFW and are better approximated by Einasto profiles. Fitting both M(<ρ>) and M(ρ_{crit}) with either NFW or Einasto profiles yield concentration and shape parameters that are correlated, confirming and extending earlier work linking the concentration of a halo with its accretion history. These correlations also confirm that halo structure is insensitive to initial conditions: only halos whose accretion histories differ greatly from the NFW shape show noticeable deviations from NFW in their mass profiles. As a result, the NFW profile provides acceptable fits to hot dark matter halos, which do not form hierarchically, and for fluctuation power spectra other than CDM. Our findings, however, predict a subtle but systematic dependence of mass profile shape on accretion history which, if confirmed, would provide strong support for the link between accretion history and halo structure we propose here.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

We use the Millennium Simulation series to investigate the mass and redshift dependence of the concentration of equilibrium cold dark matter (CDM) halos. We extend earlier work on the relation between halo mass profiles and assembly histories to show how the latter may be used to predict concentrations for halos of all masses and at any redshift. Our results clarify the link between concentration and the ``collapse redshift'' of a halo as well as why concentration depends on mass and redshift solely through the dimensionless ``peak height'' mass parameter, $ν(M,z)=δ_{\rm crit}(z)/σ(M,z)$. We combine these results with analytic mass accretion histories to extrapolate the $c(M,z)$ relations to mass regimes difficult to reach through direct simulation. Our model predicts that, at given $z$, $c(M)$ should deviate systematically from a simple power law at high masses, where concentrations approach a constant value, and at low masses, where concentrations are substantially lower than expected from extrapolating published empirical fits. This correction may reduce the expected self-annihilation boost factor from substructure by about one order of magnitude. The model also reproduces the $c(M,z)$ dependence on cosmological parameters reported in earlier work, and thus provides a simple and robust account of the relation between cosmology and the mass-concentration-redshift relation of CDM halos.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Warm DM halo mass function below the cut-off scale

Warm Dark Matter (WDM) cosmologies are a viable alternative to the Cold Dark Matter (CDM) scenario. Unfortunately, an accurate scrutiny of the WDM predictions with N-body simulations has proven difficult due to numerical artefacts. Here, we report on cosmological simulations that, for the first time, are devoid of those problems, and thus, are able to accurately resolve the WDM halo mass function well below the cut-off. We discover a complex picture, with perturbations at different evolutionary stages populating different ranges in the halo mass function. On the smallest mass scales we can resolve, identified objects are typically centres of filaments that are starting to collapse. On intermediate mass scales, objects typically correspond to fluctuations that have collapsed and are in the process of relaxation, whereas the high mass end is dominated by objects similar to haloes identified in CDM simulations. We then explicitly show how the formation of low-mass haloes is suppressed, which translates into a strong cut-off in the halo mass function. This disfavours some analytic formulations that predict a halo mass function that would extend well below the free streaming mass. We argue for a more detailed exploration of the formation of the smallest structures expected to form in a given cosmology, which, we foresee, will advance our overall understanding of structure formation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Cosmic Structure and Dynamics of the Local Universe

We present a cosmography analysis of the Local Universe based on the recently released Two-Micron All-Sky Redshift Survey (2MRS). Our method is based on a Bayesian Networks Machine Learning algorithm (the Kigen-code) which self-consistently samples the initial density fluctuations compatible with the observed galaxy distribution and a structure formation model given by second order Lagrangian perturbation theory (2LPT). From the initial conditions we obtain an ensemble of reconstructed density and peculiar velocity fields which characterize the local cosmic structure with high accuracy unveiling nonlinear structures like filaments and voids in detail. Coherent redshift space distortions are consistently corrected within 2LPT. From the ensemble of cross-correlations between the reconstructions and the galaxy field and the variance of the recovered density fields we find that our method is extremely accurate up to k ~ 1 h Mpc^-1 and still yields reliable results down to scales of about 3-4 h^-1 Mpc. The motion of the local group we obtain within ~ 80 h^-1 Mpc (v_LG=522+-86 km s^-1, l_LG=291^o +- 16^o, b_LG=34^o+-8^o) is in good agreement with measurements derived from the CMB and from direct observations of peculiar motions and is consistent with the predictions of LambdaCDM.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Cosmology with clustering anisotropies: disentangling dynamic and geometric distortions in galaxy redshift surveys

We investigate the impact of different observational effects affecting a precise and accurate measurement of the growth rate of fluctuations from the anisotropy of clustering in galaxy redshift surveys. We focus on redshift measurement errors, on the reconstruction of the underlying real-space clustering and on the apparent degeneracy existing with the geometrical distortions induced by the cosmology-dependent conversion of redshifts into distances. We use a suite of mock catalogues extracted from large N-body simulations, focusing on the analysis of intermediate, mildly non-linear scales and apply the standard linear dispersion model to fit the anisotropy of the observed correlation function. We verify that redshift errors up to ~0.2% have a negligible impact on the precision with which the specific growth rate beta can be measured. Larger redshift errors introduce a positive systematic error, which can be alleviated by adopting a Gaussian distribution function of pairwise velocities. This is, in any case, smaller than the systematic error of up to 10% due to the limitations of the linear dispersion model, which is studied in a separate paper. We then show that 50% of the statistical error budget on beta depends on the deprojection procedure through which the real-space correlation function is obtained. Finally, we demonstrate that the degeneracy with geometric distortions can in fact be circumvented. This is obtained through a modified version of the Alcock-Paczynski test in redshift-space, which successfully recovers the correct cosmology by searching for the solution that optimizes the description of dynamical redshift distortions. For a flat cosmology, we obtain largely independent, robust constraints on beta and OmegaM. In a volume of 2.4(Gpc/h)^3, the correct OmegaM is obtained with ~12% error and negligible bias, once the real-space correlation function is properly reconstructed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Estimating cosmic velocity fields from density fields and tidal tensors

In this work we investigate the nonlinear and nonlocal relation between cosmological density and peculiar velocity fields. Our goal is to provide an algorithm for the recon- struction of the nonlinear velocity field from the fully nonlinear density. We find that including the gravitational tidal field tensor using second order Lagrangian perturba- tion theory (2LPT) based upon an estimate of the linear component of the nonlinear density field significantly improves the estimate of the cosmic flow in comparison to linear theory not only in the low density, but also and more dramatically in the high density regions. In particular we test two estimates of the linear component: the log- normal model and the iterative Lagrangian linearisation. The present approach relies on a rigorous higher order Lagrangian perturbation theory analysis which incorpo- rates a nonlocal relation. It does not require additional fitting from simulations being in this sense parameter free, it is independent of statistical-geometrical optimisation and it is straightforward and efficient to compute. The method is demonstrated to yield an unbiased estimator of the velocity field on scales ~> 5 Mpc/h with closely Gaussian distributed errors. Moreover, the statistics of the divergence of the peculiar velocity field is extremely well recovered showing a good agreement with the true one from N-body simulations. The typical errors of about 10 km/s (1 sigma confidence intervals) are reduced by more than 80% with respect to linear theory in the scale range between 5 and 10 Mpc/h in high density regions (δ > 2). We also find that iterative Lagrangian linearisation is significantly superior in the low density regime with respect to the lognormal model.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Semi-analytic galaxy formation in early dark energy cosmologies

We study the impact of early dark energy (EDE) cosmologies on galaxy properties by coupling high-resolution numerical simulations with semi-analytic modeling (SAM) of galaxy formation and evolution. EDE models are characterized by a non-vanishing high-redshift contribution of dark energy, producing an earlier growth of structures and a modification of large-scale structure evolution. They can be viewed as typical representatives of non-standard dark energy models in which only the expansion history is modified, and hence the impact on galaxy formation is indirect. We show that in EDE cosmologies the predicted space density of galaxies is enhanced at all scales with respect to the standard LCDM scenario, and the corresponding cosmic star formation history and stellar mass density is increased at high-redshift. We compare these results with a set of theoretical predictions obtained with alternative SAMs applied to our reference LCDM simulation, yielding a rough measure of the systematic uncertainty of the models. We find that the modifications in galaxy properties induced by EDE cosmologies are of the same order of magnitude as intra-SAM variations for a standard LCDM realization (unless rather extreme EDE models are considered), suggesting that is difficult to use such predictions alone to disentangle between different cosmological scenarios. However, when independent information on the underlying properties of host dark matter haloes is included, the SAM predictions on galaxy bias may provide important clues on the expansion history and the equation-of-state evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The effects of alignment and ellipticity on the clustering of galaxies

We investigate the effects of halo ellipticity and alignment with larger-scale structure on the galaxy correlation function. We base our analysis on the galaxy formation models of Guo et al. (2011), run on the Millennium Simulations. We quantify the importance of these properties of the galaxy distribution by randomizing the angular positions of satellite galaxies within haloes, either coherently or individually, while keeping the distance to their respective central galaxies fixed. We find that the effect of disrupting the alignment with larger-scale structure is a ~2 per cent decrease in the galaxy correlation function around r=1.8 Mpc/h. Sphericalizing the ellipsoidal distributions of galaxies within haloes decreases the correlation function by up to 20 per cent for r<1 Mpc/h. Similar results apply to power spectra and redshift-space correlation functions. Models such as those based on the Halo Occupation Distribution, which adopt a spherically averaged profile for the galaxy distributions within haloes, will therefore significantly underestimate the clustering on sub-Mpc scales.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The effects of halo alignment and shape on the clustering of galaxies

We investigate the effects of halo shape and its alignment with larger scale structure on the galaxy correlation function. We base our analysis on the galaxy formation models of Guo et al., run on the Millennium Simulations. We quantify the importance of these effects by randomizing the angular positions of satellite galaxies within haloes, either coherently or individually, while keeping the distance to their respective central galaxies fixed. We find that the effect of disrupting the alignment with larger scale structure is a ~2 per cent decrease in the galaxy correlation function around r=1.8 Mpc/h. We find that sphericalizing the ellipsoidal distributions of galaxies within haloes decreases the correlation function by up to 20 per cent for r<1 Mpc/h and increases it slightly at somewhat larger radii. Similar results apply to power spectra and redshift-space correlation functions. Models based on the Halo Occupation Distribution, which place galaxies spherically within haloes according to a mean radial profile, will therefore significantly underestimate the clustering on sub-Mpc scales. In addition, we find that halo assembly bias, in particular the dependence of clustering on halo shape, propagates to the clustering of galaxies. We predict that this aspect of assembly bias should be observable through the use of extensive group catalogues.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

One simulation to fit them all - changing the background parameters of a cosmological N-body simulation

We demonstrate that the output of a cosmological N-body simulation can, to remarkable accuracy, be scaled to represent the growth of large-scale structure in a cosmology with parameters similar to but different from those originally assumed. Our algorithm involves three steps: a reassignment of length, mass and velocity units, a relabelling of the time axis, and a rescaling of the amplitudes of individual large-scale fluctuation modes. We test it using two matched pairs of simulations. Within each pair, one simulation assumes parameters consistent with analyses of the first-year WMAP data. The other has lower matter and baryon densities and a 15% lower fluctuation amplitude, consistent with analyses of the three-year WMAP data. The pairs differ by a factor of a thousand in mass resolution, enabling performance tests on both linear and nonlinear scales. Our scaling reproduces the mass power spectra of the target cosmology to better than 0.5% on large scales (k < 0.1 h/Mpc) both in real and in redshift space. In particular, the BAO features of the original cosmology are removed and are correctly replaced by those of the target cosmology. Errors are still below 3% for k < 1 h/Mpc. Power spectra of the dark halo distribution are even more precisely reproduced, with errors below 1% on all scales tested. A halo-by-halo comparison shows that centre-of-mass positions and velocities are reproduced to better than 90 kpc/h and 5%, respectively. Halo masses, concentrations and spins are also reproduced at about the 10% level, although with small biases. Halo assembly histories are accurately reproduced, leading to central galaxy magnitudes with errors of about 0.25 magnitudes and a bias of about 0.13 magnitudes for a representative semi-analytic model.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Simulations of Quintessential Cold Dark Matter: beyond the cosmological constant

We study the nonlinear growth of cosmic structure in different dark energy models, using large volume N-body simulations. We consider a range of quintessence models which feature both rapidly and slowly varying dark energy equations of state, and compare the growth of structure to that in a universe with a cosmological constant. The adoption of a quintessence model changes the expansion history of the universe, the form of the linear theory power spectrum and can alter key observables, such as the horizon scale and the distance to last scattering. We incorporate these effects into our simulations in stages to isolate the impact of each on the growth of structure. The difference in structure formation can be explained to first order by the difference in growth factor at a given epoch; this scaling also accounts for the nonlinear growth at the 15% level. We find that quintessence models that are different from $Λ$CDM both today and at high redshifts $(z \sim 1000)$ and which feature late $(z<2)$, rapid transitions in the equation of state, can have identical baryonic acoustic oscillation (BAO) peak positions to those in $Λ$CDM. We find that these models have higher abundances of dark matter haloes at $z>0$ compared to $Λ$CDM and so measurements of the mass function should allow us to distinguish these quintessence models from a cosmological constant. However, we find that a second class of quintessence models, whose equation of state makes an early $(z>2)$ rapid transition to $w=-1$, cannot be distinguished from $Λ$CDM using measurements of the mass function or the BAO, even if these models have non-negligible amounts of dark energy at early times.

0 citations0 reviewsNo code linkedOpen access