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

Anatoly Klypin

Anatoly Klypin 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-phastro-ph.HEastro-ph.IMDatabaseshep-ph

Catalog footprint

What is connected

35works
7topics
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

35 published item(s)

preprint2022arXiv

Fast full N-body simulations of generic modified gravity: conformal coupling models

We present MG-GLAM, a code developed for the very fast production of full $N$-body cosmological simulations in modified gravity (MG) models. We describe the implementation, numerical tests and first results of a large suite of cosmological simulations for three classes of MG models with conformal coupling terms: the $f(R)$ gravity, symmetron and coupled quintessence models. Derived from the parallel particle-mesh code GLAM, MG-GLAM incorporates an efficient multigrid relaxation technique to solve the characteristic nonlinear partial differential equations of these models. For $f(R)$ gravity, we have included new variants to diversify the model behaviour, and we have tailored the relaxation algorithms to these to maintain high computational efficiency. In a companion paper, we describe versions of this code developed for derivative coupling MG models, including the Vainshtein- and K-mouflage-type models. MG-GLAM can model the prototypes for most MG models of interest, and is broad and versatile. The code is highly optimised, with a tremendous speedup of a factor of more than a hundred compared with earlier $N$-body codes, while still giving accurate predictions of the matter power spectrum and dark matter halo abundance. MG-GLAM is ideal for the generation of large numbers of MG simulations that can be used in the construction of mock galaxy catalogues and the production of accurate emulators for ongoing and future galaxy surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Fast full $N$-body simulations of generic modified gravity: derivative coupling models

We present MG-GLAM, a code developed for the very fast production of full $N$-body cosmological simulations in modified gravity (MG) models. We describe the implementation, numerical tests and first results of a large suite of cosmological simulations for two broad classes of MG models with derivative coupling terms -- the Vainshtein- and Kmouflage-type models -- which respectively features the Vainshtein and Kmouflage screening mechanism. Derived from the parallel particle-mesh code GLAM, MG-GLAM incorporates an efficient multigrid relaxation technique to solve the characteristic nonlinear partial differential equations of these models. For Kmouflage, we have proposed a new algorithm for the relaxation solver, and run the first simulations of the model to understand its cosmological behaviour. In a companion paper, we describe versions of this code developed for conformally-coupled MG models, including several variants of $f(R)$ gravity, the symmetron model and coupled quintessence. Altogether, MG-GLAM has so far implemented the prototypes for most MG models of interest, and is broad and versatile. The code is highly optimised, with a tremendous (over two orders of magnitude) speedup when comparing its running time with earlier $N$-body codes, while still giving accurate predictions of the matter power spectrum and dark matter halo abundance. MG-GLAM is ideal for the generation of large numbers of MG simulations that can be used in the construction of mock galaxy catalogues and accurate emulators for ongoing and future galaxy surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Building a digital twin of a luminous red galaxy spectroscopic survey: galaxy properties and clustering covariance

Luminous red galaxies (LRGs) are one of the key tracers of the large-scale structure of the Universe used by galaxy surveys. Hence, it is important to make accurate predictions for their properties and clustering, including the errors on these statistics. Here, we describe a novel technique which uses the semi-analytical model of galaxy formation {\sc Galform}, embedded in the high-resolution $N$-body Planck-Millennium simulation, to populate a thousand halo catalogues generated using the Parallel-PM $N$-body {\sc glam} code. Our hybrid scheme allows us to make clustering predictions on scales that cannot be modelled in the original $N$-body simulation. LRGs are selected in the redshift range $z=0.6-1$ from the {\sc Galform} output using similar colour-magnitude cuts in the $r$, $z$ and $W1$ bands to those that will be applied in the Dark Energy Spectroscopic Instrument (DESI) survey. We find that the LRG-halo connection is non-trivial, leading to the prediction of a non-standard halo occupation distribution; in particular, the occupation of central galaxies does not reach unity for the most massive haloes, and drops with increasing mass. The {\sc glam} catalogues reproduce the abundance and clustering of the LRGs predicted by {\sc Galform}, and show good agreement with recent measurements of the clustering of DESI-like LRGs using photometric redshifts. We use the \glam{} mocks to compute the covariance matrices for the two-point correlation function and power spectrum of the LRGs and their background dark matter density field, revealing important differences. We also make predictions for the linear-growth rate and the baryon acoustic oscillations distances at $z=0.6$, $0.74$ and $0.93$. All DESI-like LRG catalogues are made publicly available.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Completed SDSS-IV Extended Baryon Oscillation Spectroscopic Survey: GLAM-QPM mock galaxy catalogs for the Emission Line Galaxy Sample

We present 2000 mock galaxy catalogs for the analysis of baryon acoustic oscillations in the Emission Line Galaxy (ELG) sample of the Extended Baryon Oscillation Spectroscopic Survey Data Release 16 (eBOSS DR16). Each mock catalog has a number density of $6.7 \times 10^{-4} h^3 \rm Mpc^{-3}$, covering a redshift range from 0.6 to 1.1. The mocks are calibrated to small-scale eBOSS ELG clustering measurements at scales of around 10 $h^{-1}$Mpc. The mock catalogs are generated using a combination of GaLAxy Mocks (GLAM) simulations and the Quick Particle-Mesh (QPM) method. GLAM simulations are used to generate the density field, which is then assigned dark matter halos using the QPM method. Halos are populated with galaxies using a halo occupation distribution (HOD). The resulting mocks match the survey geometry and selection function of the data, and have slightly higher number density which allows room for systematic analysis. The large-scale clustering of mocks at the baryon acoustic oscillation (BAO) scale is consistent with data and we present the correlation matrix of the mocks.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

GalWeight Application: A publicly-available catalog of dynamical parameters of 1,800 galaxy clusters from SDSS-DR13, ($\mathtt{GalWCat19}$)

Utilizing the SDSS-DR13 spectroscopic dataset, we create a new publicly-available catalog of 1,800 galaxy clusters (GalWeight cluster catalog, $\mathtt{GalWCat19}$) and a corresponding catalog of 34,471 identified member galaxies. The clusters are identified from overdensities in redshift-phase space. The GalWeight technique introduced in Abdullah, Wilson and Klypin (AWK18) is then applied to identify cluster members. The completeness of the cluster catalog ($\mathtt{GalWCat19}$) and the procedure followed to determine cluster mass are tested on the Bolshoi N-body simulations. The 1,800 $\mathtt{GalWCat19}$ clusters range in redshift between $0.01 - 0.2$ and in mass between $(0.4 - 14) \times 10^{14}h^{-1}M_{\odot}$. The cluster catalog provides a large number of cluster parameters including sky position, redshift, membership, velocity dispersion, and mass at overdensities $Δ= 500, 200, 100, 5.5$. The 34,471 member galaxies are identified within the radius at which the density is 200 times the critical density of the Universe. The galaxy catalog provides the coordinates of each galaxy and the ID of the cluster that the galaxy belongs to. The cluster velocity dispersion scales with mass as $\log(σ_{200})=\log(946\pm52~ \mbox{km} ~ \mbox{s}^{-1}) +(0.349\pm0.142)\log\left[h(z) ~ M_{200}/10^{15}M_\odot\right]$ with scatter of $δ_{\logσ} = 0.06$. The catalogs are publicly available at the following website\footnote{\url{https://mohamed-elhashash-94.webself.net/galwcat/}}.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Suppressing cosmic variance with paired-and-fixed cosmological simulations: average properties and covariances of dark matter clustering statistics

Making cosmological inferences from the observed galaxy clustering requires accurate predictions for the mean clustering statistics and their covariances. Those are affected by cosmic variance -- the statistical noise due to the finite number of harmonics. The cosmic variance can be suppressed by fixing the amplitudes of the harmonics instead of drawing them from a Gaussian distribution predicted by the inflation models. Initial realizations also can be generated in pairs with 180 degrees flipped phases to further reduce the variance. Here, we compare the consequences of using paired-and-fixed vs Gaussian initial conditions on the average dark matter clustering and covariance matrices predicted from N-body simulations. As in previous studies, we find no measurable differences between paired-and-fixed and Gaussian simulations for the average density distribution function, power spectrum and bispectrum. Yet, the covariances from paired-and-fixed simulations are suppressed in a complicated scale- and redshift-dependent way. The situation is particularly problematic on the scales of Baryon Acoustic Oscillations where the covariance matrix of the power spectrum is lower by only 20% compared to the Gaussian realizations, implying that there is not much of a reduction of the cosmic variance. The non-trivial suppression, combined with the fact that paired-and-fixed covariances are noisier than from Gaussian simulations, suggests that there is no path towards obtaining accurate covariance matrices from paired-and-fixed simulations. Because the covariances are crucial for the observational estimates of galaxy clustering statistics and cosmological parameters, paired-and-fixed simulations, though useful for some applications, cannot be used for the production of mock galaxy catalogs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Abundance of Field Galaxies

We present new measurements of the abundance of galaxies with a given circular velocity in the Local Volume: a region centered on the Milky Way Galaxy and extending to distance 10Mpc. The sample of 750 mostly dwarf galaxies provides a unique opportunity to study the abundance and properties of galaxies down to absolute magnitudes MB= -10, and virial masses Mvir= 1e9Msun. We find that the standard LCDM model gives remarkably accurate estimates for the velocity function of galaxies with circular velocities V>60kms and corresponding virial masses Mvir> 3e10Msun, but it badly fails by over-predicting 5 times the abundance of large dwarfs with velocities V= 30-50kms. The Warm Dark Matter models cannot explain the data either, regardless of mass of the WDM particle. Just as in previous observational studies, we find a shallow asymptotic slope dN/dlog V = V**alpha, alpha =-1 of the velocity function, which is inconsistent with the standard LCDM model that predicts the slope alpha =-3. Though reminiscent to the known overabundance of satellites problem, the overabundance of field galaxies is a much more difficult problem. For the LCDM model to survive, in the 10Mpc radius of the Milky Way there should be 1000 dark galaxies with virial mass Mvir= 1e10Msun, extremely low surface brightness and no detectable HI gas. So far none of this type of galaxies have been discovered.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Clustering properties of $g$-selected galaxies at $z\sim0.8$

Current and future large redshift surveys, as the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey (SDSS-IV/eBOSS) or the Dark Energy Spectroscopic Instrument (DESI), will use emission-line galaxies (ELG) to probe cosmological models by mapping the large-scale structure of the Universe in the redshift range $0.6 < z < 1.7$. With current data, we explore the halo-galaxy connection by measuring three clustering properties of $g$-selected ELGs as matter tracers in the redshift range $0.6 < z < 1$: (i) the redshift-space two-point correlation function using spectroscopic redshifts from the BOSS ELG sample and VIPERS; (ii) the angular two-point correlation function on the footprint of the CFHT-LS; (iii) the galaxy-galaxy lensing signal around the ELGs using the CFHTLenS. We interpret these observations by mapping them onto the latest high-resolution MultiDark Planck N-body simulation, using a novel (Sub)Halo-Abundance Matching technique that accounts for the ELG incompleteness. ELGs at $z\sim0.8$ live in halos of $(1\pm 0.5)\times10^{12}\,h^{-1}$M$_{\odot}$ and 22.5$\pm2.5$% of them are satellites belonging to a larger halo. The halo occupation distribution of ELGs indicates that we are sampling the galaxies in which stars form in the most efficient way, according to their stellar-to-halo mass ratio.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Galaxy Three-Point Correlation Functions and Halo/Subhalo Models

We present the measurements of the luminosity-dependent redshift-space three-point correlation functions (3PCFs) for the Sloan Digital Sky Survey (SDSS) DR7 Main galaxy sample. We compare the 3PCF measurements to the predictions from three different halo and subhalo models. One is the halo occupation distribution (HOD) model and the other two are extensions of the subhalo abundance matching (SHAM) model by allowing the central and satellite galaxies to have different occupation distributions in the host halos and subhalos. Parameters in all the models are chosen to best describe the projected and redshift-space two-point correlation functions (2PCFs) of the same set of galaxies. All three model predictions agree well with the 3PCF measurements for the most luminous galaxy sample, while the HOD model better performs in matching the 3PCFs of fainter samples (with luminosity threshold below $L^*$), which is similar in trend to the case of fitting the 2PCFs. The decomposition of the model 3PCFs into contributions from different types of galaxy triplets shows that on small scales the dependence of the 3PCFs on triangle shape is driven by nonlinear redshift-space distortion (and not by the intrinsic halo shape) while on large scales it reflects the filamentary structure. The decomposition also reveals more detailed differences in the three models, which are related to the radial distribution, the mean occupation function, and the velocity distribution of satellite galaxies inside halos. The results suggest that galaxy 3PCFs can further help constrain the above galaxy-halo relation and test theoretical models.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Modelling galaxy clustering: halo occupation distribution versus subhalo matching

We model the luminosity-dependent projected and redshift-space two-point correlation functions (2PCFs) of the Sloan Digital Sky Survey (SDSS) DR7 Main galaxy sample, using the halo occupation distribution (HOD) model and the subhalo abundance matching (SHAM) model and its extension. All the models are built on the same high-resolution $N$-body simulations. We find that the HOD model generally provides the best performance in reproducing the clustering measurements in both projected and redshift spaces. The SHAM model with the same halo-galaxy relation for central and satellite galaxies (or distinct haloes and subhaloes), when including scatters, has a best-fitting $χ^2/\rm{dof}$ around $2$--$3$. We therefore extend the SHAM model to the subhalo clustering and abundance matching (SCAM) by allowing the central and satellite galaxies to have different galaxy--halo relations. We infer the corresponding halo/subhalo parameters by jointly fitting the galaxy 2PCFs and abundances and consider subhaloes selected based on three properties, the mass $M_{\rm acc}$ at the time of accretion, the maximum circular velocity $V_{\rm acc}$ at the time of accretion, and the peak maximum circular velocity $V_{\rm peak}$ over the history of the subhaloes. The three subhalo models work well for luminous galaxy samples (with luminosity above $L_*$). For low-luminosity samples, the $V_{\rm acc}$ model stands out in reproducing the data, with the $V_{\rm peak}$ model slightly worse, while the $M_{\rm acc}$ model fails to fit the data. We discuss the implications of the modeling results.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

MultiDark simulations: the story of dark matter halo concentrations and density profiles

Accurately predicting structural properties of dark matter halos is one of the fundamental goals of modern cosmology. We use the new suite of MultiDark cosmological simulations to study the evolution of dark matter halo density profiles, concentrations, and velocity anisotropies. The MultiDark simulations cover a large range of masses 1e10-1e15Msun and volumes upto 50Gpc**3. The total number of dark matter halos in all the simulations exceeds 60 billion. We find that in order to understand the structure of dark matter halos and to make ~1% accurate predictions for density profiles, one needs to realize that halo concentration is more complex than the traditional ratio of the virial radius to the core radius in the NFW profile. For massive halos the averge density profile is far from the NFW shape and the concentration is defined by both the core radius and the shape parameter alpha in the Einasto approximation. Combining results from different redshifts, masses and cosmologies, we show that halos progress through three stages of evolution. (1) They start as rare density peaks that experience very fast and nearly radial infall. This radial infall brings mass closer to the center producing a high concentrated halo. Here, the halo concentration increases with the increasing halo mass and the concentration is defined by the alpha parameter with nearly constant core radius. Later halos slide into (2) the plateau regime where the accretion becomes less radial, but frequent mergers still affect even the central region. Now the concentration does not depend on halo mass. (3) Once the rate of accretion slows down, halos move into the domain of declining concentration-mass relation because new accretion piles up mass close to the virial radius while the core radius is staying constant. We provide accurate analytical fits to the numerical results for halo density profiles and concentrations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: mock galaxy catalogues for the BOSS Final Data Release

We reproduce the galaxy clustering catalogue from the SDSS-III Baryon Oscillation Spectroscopic Survey Final Data Release (BOSS DR11&DR12) with high fidelity on all relevant scales in order to allow a robust analysis of baryon acoustic oscillations and redshift space distortions. We have generated (6,000) 12,288 MultiDark PATCHY BOSS (DR11) DR12 light-cones corresponding to an effective volume of $\sim192,000\,[h^{-1}\,{\rm Gpc}]^3$ (the largest ever simulated volume), including cosmic evolution in the redshift range from 0.15 to 0.75. The mocks have been calibrated using a reference galaxy catalogue based on the halo abundance matching modelling of the BOSS DR11&DR12 galaxy clustering data and on the data themselves. The production follows three steps. First, we apply the PATCHY code to generate a dark matter field and an object distribution including nonlinear stochastic galaxy bias. Secondly, we run the halo/stellar distribution reconstruction HADRON code to assign masses to the various objects. This step uses the mass distribution as a function of local density and non-local indicators (i.e., tidal field tensor eigenvalues and relative halo exclusion separation for massive objects) from the reference simulation applied to the corresponding PATCHY dark matter and galaxy distribution. Finally, we apply the SUGAR code to build the light cones. The resulting MultiDark PATCHY mock light cones reproduce the number density, selection function, survey geometry, and in general within 1 $σ$, for arbitrary stellar mass bins, the power spectrum up to $k=0.3\,h\,{\rm Mpc}^{-1}$, the two-point correlation functions down to a few Mpc scales, and the three-point statistics of the BOSS DR11&DR12 galaxy samples.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Modeling the clustering and halo occupation distribution of BOSS-CMASS galaxies in the Final Data Release

We present a study of the clustering and halo occupation distribution of BOSS CMASS galaxies in the redshift range 0.43 < z < 0.7 drawn from the Final SDSS-III Data Release. We compare the BOSS results with the predictions of a Halo Abundance Matching (HAM) clustering model that assigns galaxies to dark matter halos selected from the large BigMultiDark $N$-body simulation of a flat $Λ$CDM Planck cosmology. We compare the observational data with the simulated ones on a light-cone constructed from 20 subsequent outputs of the simulation. Observational effects such as incompleteness, geometry, veto masks and fiber collisions are included in the model, which reproduces within 1-$σ$ errors the observed monopole of the 2-point correlation function at all relevant scales: from the smallest scales, 0.5 $h^{-1}$ Mpc, up to scales beyond the Baryonic Acoustic Oscillation feature. This model also agrees remarkably well with the BOSS galaxy power spectrum (up to $k\sim1$ $h$ Mpc$^{-1}$), and the Three-point correlation function. The quadrupole of the correlation function presents some tensions with observations. We discuss possible causes that can explain this disagreement, including target selection effects. Overall, the standard HAM model describes remarkably well the clustering statistics of the CMASS sample. We compare the stellar to halo mass relation for the CMASS sample measured using weak lensing in the CFHT Stripe 82 Survey with the prediction of our clustering model, and find a good agreement within 1-$σ$. The BigMD-BOSS light-cone including properties of BOSS galaxies and halo properties is made publicly available.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Direct Insights into Observational Absorption Line Analysis Methods of the Circumgalactic Medium Using Cosmological Simulations

We study the circumgalactic medium (CGM) of a z=0.54 simulated dwarf galaxy using hydroART simulations. We present our analysis methods, which emulate observations, including objective absorption line detection, apparent optical depth (AOD) measurements, Voigt profile (VP) decomposition, and ionization modeling. By comparing the inferred CGM gas properties from the absorption lines directly to the gas selected by low ionization HI and MgII, and by higher ionization CIV and OVI absorption, we examine how well observational analysis methods recover the "true" properties of CGM gas. In this dwarf galaxy, low ionization gas arises in sub-kiloparsec "cloud" structures, but high ionization gas arises in multiple extended structures spread over 100 kpc; due to complex velocity fields, highly separated structures give rise to absorption at similar velocities. We show that AOD and VP analysis fails to accurately characterize the spatial, kinematic, and thermal conditions of high ionization gas. We find that HI absorption selected gas and OVI absorption gas arise in totally distinct physical gas structures, calling into question current observational techniques employed to infer metallicities and the total mass of "warm-hot" CGM gas. We present a method to determine whether CIV and OVI absorbing gas is photo or collisionally ionized and whether the assumption of ionization equilibrium is sound. As we discuss, these and additional findings have strong implications for how accurately currently employed observational absorption line methods recover the true gas properties, and ultimately, our ability to understand the CGM and its role in galaxy evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Low-mass galaxy assembly in simulations: regulation of early star formation by radiation from massive stars

Despite recent success in forming realistic present-day galaxies, simulations still form the bulk of their stars earlier than observations indicate. We investigate the process of stellar mass assembly in low-mass field galaxies, a dwarf and a typical spiral, focusing on the effects of radiation from young stellar clusters on the star formation (SF) histories. We implement a novel model of SF with a deterministic low efficiency per free-fall time, as observed in molecular clouds. Stellar feedback is based on observations of star-forming regions, and includes radiation pressure from massive stars, photoheating in H II regions, supernovae and stellar winds. We find that stellar radiation has a strong effect on the formation of low-mass galaxies, especially at z > 1, where it efficiently suppresses SF by dispersing cold and dense gas, preventing runaway growth of the stellar component. This behaviour is evident in a variety of observations but had so far eluded analytical and numerical models without radiation feedback. Compared to supernovae alone, radiation feedback reduces the SF rate by a factor of ~100 at z < 2, yielding rising SF histories which reproduce recent observations of Local Group dwarfs. Stellar radiation also produces bulgeless spiral galaxies and may be responsible for excess thickening of the stellar disc. The galaxies also feature rotation curves and baryon fractions in excellent agreement with current data. Lastly, the dwarf galaxy shows a very slow reduction of the central dark matter density caused by radiation feedback over the last ~7 Gyr of cosmic evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

What controls the ionized gas turbulent motions in dwarf galaxies?

Using 3D spectroscopy with a scanning Fabry-Perot interferometer, we study the ionized gas kinematics in 59 nearby dwarf galaxies. Combining our results with data from literature, we provide a global relation between the gas velocity dispersion (sigma) and the star formation rate (SFR) and Hαluminosity for galaxies in a very broad range of star formation rates SFR=0.001-300 Msun/yr. We find that the SFR-sigma relation for the combined sample of dwarf galaxies, star forming, local luminous, and ultra-luminous infrared galaxies can be fitted as sigma~ SFR^(5.3+-0.2). This implies that the slope of the L-sigma relation inferred from the sample of rotation supported disc galaxies (including mergers) is similar to the L-sigma relation of individual giant HII regions. We present arguments that the velocity dispersion of the ionized gas does not reflect the virial motions in the gravitational potential of dwarf galaxies, and instead is mainly determined by the energy injected into the interstellar medium by the ongoing star formation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Radiative feedback and the low efficiency of galaxy formation in low-mass haloes at high redshift

Any successful model of galaxy formation needs to explain the low rate of star formation in the small progenitors of today's galaxies. This inefficiency is necessary for reproducing the low stellar-to-virial mass fractions, suggested by current abundance matching models. A possible driver of this low efficiency is the radiation pressure exerted by ionizing photons from massive stars. The effect of radiation pressure in cosmological, zoom-in galaxy formation simulations is modeled as a non-thermal pressure that acts only in dense and optically thick star-forming regions. We also include photoionization and photoheating by massive stars. The full photoionization of hydrogen reduces the radiative cooling in the $10^{4-4.5}$ K regime. The main effect of radiation pressure is to regulate and limit the high values of gas density and the amount of gas available for star formation. This maintains a low star formation rate of $\sim 1 \ {\rm M_\odot} \ {\rm yr}^{-1}$ in halos with masses about $10^{11} \ {M_\odot}$ at $z\simeq3$. Infrared trapping and photoionization/photoheating processes are secondary effects in this mass range. The galaxies residing in these low-mass halos contain only $\sim0.6\%$ of the total virial mass in stars, roughly consistent with abundance matching. Radiative feedback maintains an extended galaxy with a rising circular velocity profile.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Constraints on WIMP Annihilation for Contracted Dark Matter in the Inner Galaxy with the Fermi-LAT

We derive constraints on parameters of generic dark matter candidates by comparing theoretical predictions with the gamma-ray emission observed by the Fermi-LAT from the region around the Galactic Center. Our analysis is conservative since it simply requires that the expected dark matter signal does not exceed the observed emission. The constraints obtained in the likely case that the collapse of baryons to the Galactic Center is accompanied by the contraction of the dark matter are strong. In particular, we find that for bb and τ+τ- or W+W- dark matter annihilation channels, the upper limits on the annihilation cross section imply that the thermal cross section is excluded for a Weakly Interacting Massive Particle (WIMP) mass smaller than about 700 and 500 GeV, respectively. For the μ+ μ- channel, where the effect of the inverse Compton scattering is important, depending on models of the Galactic magnetic field the exclusion of the thermal cross-section is for a WIMP mass smaller than about 150 to 400 GeV. The upper limits on the annihilation cross section of dark matter particles obtained are two orders of magnitude stronger than without contraction. In the latter case our results are compatible with the upper limits from the Galactic halo analysis reported by the Fermi-LAT collaboration for the case in which the same conservative approach without modeling of the astrophysical background is employed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Effects of baryon removal on the structure of dwarf spheroidal galaxies

Dwarf spheroidal galaxies (dSphs) are extremely gas-poor, dark matter-dominated galaxies, which make them ideal to test the predictions of the cold dark matter (CDM) model. We argue that the removal of the baryonic component from gas-rich dwarf irregular galaxies, the progenitors of dSphs, can substantially reduce their central density. Thus, it may play an important role in alleviating one of the problems of the CDM model related to the structure of relatively massive satellite galaxies of the Milky Way (MW). Traditionally, collisionless cosmological N-body simulations are used when confronting theoretical predictions with observations. However, these simulations assume that the baryon fraction everywhere in the Universe is equal to the cosmic mean, which can be incorrect for dSphs. We find that the combination of (i) the lower baryon fraction in dSphs compared to the cosmic mean and (ii) the concentration of baryons in the inner part of the MW halo can go a long way towards explaining the observed circular velocity profiles of dSphs. We find that the blowing away of baryons by ram pressure, when the dwarfs fall into larger galaxies, decreases the circular velocity profile of the satellite and reduces the density in the central ~200-500 pc by a factor of (1 - f_b)^4 ~ 0.5, where f_b is the cosmological fraction of baryons. Additionally, the enhanced baryonic mass in the central regions of the parent galaxy generates tidal forces, which are larger than those experienced by subhaloes in traditional N-body simulations. Increased tidal forces substantially alter circular velocity profiles for satellites with pericentres less than 50 kpc. We show that these two effects are strong enough to bring the predictions of subhaloes from CDM simulations into agreement with the observed structure of MW dSphs, regardless of the details of the baryonic processes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Halo Abundance Matching: accuracy and conditions for numerical convergence

Accurate predictions of the abundance and clustering of dark matter haloes play a key role in testing the standard cosmological model. Here, we investigate the accuracy of one of the leading methods of connecting the simulated dark matter haloes with observed galaxies -- the Halo Abundance Matching (HAM) technique. We show how to choose the optimal values of the mass and force resolution in large-volume N-body simulations so that they provide accurate estimates for correlation functions and circular velocities for haloes and their subhaloes -- crucial ingredients of the HAM method. At the 10% accuracy, results converge for 50 particles for haloes and 150 particles for progenitors subhaloes. In order to achieve this level of accuracy a number of conditions should be satisfied. The force resolution for the smallest resolved (sub)haloes should be in the range (0.1-0.3)rs, where rs is the scale radius of (sub)haloes. The number of particles for progenitors of subhaloes should be 150. We also demonstrate that the two-body scattering plays a minor role for the accuracy of N-body simulations thanks to the relatively small number of crossing-times of dark matter in haloes, and the limited force resolution of cosmological simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Mergers and Mass Accretion for Infalling Halos Both End Well Outside Cluster Virial Radii

We find that infalling dark matter halos (i.e., the progenitors of satellite halos) begin losing mass well outside the virial radius of their eventual host halos. The peak mass occurs at a range of clustercentric distances, with median and 68th percentile range of 1.8 +2.3/-1.0 R_(vir,host) for progenitors of z=0 satellites. The peak circular velocity for infalling halos occurs at significantly larger distances (3.7 +3.3/-2.2 R_(vir,host) at z=0). This difference arises because different physical processes set peak circular velocity (typically, ~1:5 and larger mergers which cause transient circular velocity spikes) and peak mass (typically, smooth accretion) for infalling halos. We find that infalling halos also stop having significant mergers well before they enter the virial radius of their eventual hosts. Mergers larger than a 1:40 ratio in halo mass end for infalling halos at similar clustercentric distances (~ 1.9 R_(vir,host)) as the end of overall mass accretion. However, mergers larger than 1:3 typically end for infalling halos at more than 4 virial radial away from their eventual hosts. This limits the ability of mergers to affect quenching and morphology changes in clusters. We also note that the transient spikes which set peak circular velocity may lead to issues with abundance matching on that parameter, including unphysical galaxy stellar mass growth profiles near clusters; we propose a simple observational test to check if a better halo proxy for galaxy stellar mass exists.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Orbital anisotropy in cosmological haloes revisited

The velocity anisotropy of particles inside dark matter (DM) haloes is an important physical quantity, which is required for the accurate modelling of mass profiles of galaxies and clusters of galaxies. It is typically measured using the ratio of the radial-to-tangential velocity dispersions at a given distance from the halo centre. However, this measure is insufficient to describe the dynamics of realistic haloes, which are typically quite elongated. Studying the velocity distribution in massive DM haloes in cosmological simulations, we find that in the inner parts of the haloes the local velocity ellipsoids are strongly aligned with the major axis of the halo, the alignment being stronger for more relaxed haloes. In the outer regions of the haloes, the alignment becomes gradually weaker and the orientation is more random. These two distinct regions of different degree of the alignment coincide with two characteristic regimes of the DM density profile: shallower and steeper than ρr^{-2}. This alignment of the local velocity ellipsoids requires reinterpretation of features found in measurements based on the spherically averaged ratio of the radial-to-tangential velocity dispersions. In particular, we show that the velocity distribution in the central halo regions is highly anisotropic. For cluster-size haloes with mass 10^{14}-10^{15} h^-1 Msun, the velocity anisotropy along the major axis is nearly independent of radius and is equal to β=1-σ^2_{perp}/σ^2_{radial}=0.4, which is significantly larger than the previously estimated spherically averaged velocity anisotropy. The alignment of density and velocity anisotropies, and the radial trends may also have some implications for the mass modelling based on kinematical data of such objects as galaxy clusters or dwarf spheroidals, where the orbital anisotropy is a key element in an unbiased mass inference.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Structure Finding in Cosmological Simulations: The State of Affairs

The ever increasing size and complexity of data coming from simulations of cosmic structure formation demands equally sophisticated tools for their analysis. During the past decade, the art of object finding in these simulations has hence developed into an important discipline itself. A multitude of codes based upon a huge variety of methods and techniques have been spawned yet the question remained as to whether or not they will provide the same (physical) information about the structures of interest. Here we summarize and extent previous work of the "halo finder comparison project": we investigate in detail the (possible) origin of any deviations across finders. To this extent we decipher and discuss differences in halo finding methods, clearly separating them from the disparity in definitions of halo properties. We observe that different codes not only find different numbers of objects leading to a scatter of up to 20 per cent in the halo mass and Vmax function, but also that the particulars of those objects that are identified by all finders differ. The strength of the variation, however, depends on the property studied, e.g. the scatter in position, bulk velocity, mass, and the peak value of the rotation curve is practically below a few per cent, whereas derived quantities such as spin and shape show larger deviations. Our study indicates that the prime contribution to differences in halo properties across codes stems from the distinct particle collection methods and -- to a minor extent -- the particular aspects of how the procedure for removing unbound particles is implemented. We close with a discussion of the relevance and implications of the scatter across different codes for other fields such as semi-analytical galaxy formation models, gravitational lensing, and observables in general.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The clustering of galaxies at z~0.5 in the SDSS-III Data Release 9 BOSS-CMASS sample: a test for the LCDM cosmology

We present results on the clustering of 282,068 galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS) sample of massive galaxies with redshifts 0.4<z<0.7 which is part of the Sloan Digital Sky Survey III project. Our results cover a large range of scales from ~0.5 to ~90 Mpc/h. We compare these estimates with the expectations of the flat LCDM cosmological model with parameters compatible with WMAP7 data. We use the MultiDark cosmological simulation together with a simple halo abundance matching technique, to estimate galaxy correlation functions, power spectra, abundance of subhaloes and galaxy biases. We find that the LCDM model gives a reasonable description to the observed correlation functions at z~0.5, which is a remarkably good agreement considering that the model, once matched to the observed abundance of BOSS galaxies, does not have any free parameters. However, we find a deviation (>~10%) in the correlation functions for scales less than ~1 Mpc/h and ~10-40 Mpc/h. A more realistic abundance matching model and better statistics from upcoming observations are needed to clarify the situation. We also estimate that about 12% of the "galaxies" in the abundance-matched sample are satellites inhabiting central haloes with mass M>~1e14 M_sun/h. Using the MultiDark simulation we also study the real space halo bias b(r) of the matched catalogue finding that b=2.00+/-0.07 at large scales, consistent with the one obtained using the measured BOSS projected correlation function. Furthermore, the linear large-scale bias depends on the number density n of the abundance-matched sample as b=-0.048-(0.594+/-0.02)*log(n/(h/Mpc)^3). Extrapolating these results to BAO scales we measure a scale-dependent damping of the acoustic signal produced by non-linear evolution that leads to ~2-4% dips at ~3 sigma level for wavenumbers k>~0.1 h/Mpc in the linear large-scale bias.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The velocity shear tensor: tracer of halo alignment

The alignment of DM halos and the surrounding large scale structure (LSS) is examined in the context of the cosmic web. Halo spin, shape and the orbital angular momentum of subhaloes is investigated relative to the LSS using the eigenvectors of the velocity shear tensor evaluated on a grid with a scale of 1 Mpc/h, deep within the non-linear regime. Knots, filaments, sheets and voids are associated with regions that are collapsing along 3, 2, 1 or 0 principal directions simultaneously. Each halo is tagged with a web classification (i.e. knot halo, filament halo, etc) according to the nature of the collapse at the halo's position. The full distribution of shear eigenvalues is found to be substantially different from that tagged to haloes, indicating that the observed velocity shear is significantly biased. We find that larger mass haloes live in regions where the shear is more isotropic, namely the expansion or collapse is more spherical. A correlation is found between the halo's shape and the eigenvectors of the shear tensor, with the longest (shortest) axis of the halo's shape being aligned with the slowest (fastest) collapsing eigenvector. This correlation is web independent, suggesting that the velocity shear is a fundamental tracer of the halo alignment. A similar result is found for the alignment of halo spin with the cosmic web. It has been shown that high mass haloes exhibit a spin flip with respect to the LSS: we find the mass at which this spin flip occurs is web dependent and not universal as suggested previously. Although weaker than haloes, subhalo orbits too exhibit an alignment with the LSS, providing a possible insight into the highly correlated co-rotation of the Milky Way's satellite system. The present study suggests that the velocity shear tensor constitutes the natural framework for studying the directional properties of the non-linear LSS and of halos and galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Galaxies in LCDM with Halo Abundance Matching: luminosity-velocity relation, baryonic mass-velocity relation, velocity function and clustering

It has long been regarded as difficult for a cosmological model to account simultaneously for the galaxy luminosity, mass, and velocity distributions. We revisit this issue using a modern compilation of observational data along with the best available large-scale cosmological simulation of dark matter. We find that the standard cosmological model, used in conjunction with halo abundance matching (HAM) and simple dynamical corrections, fits all basic statistics of galaxies with circular velocities Vcirc > 80 km/s. Our observational constraint is the luminosity-velocity relation which allows all types of galaxies to be included. We have compiled data for a variety of galaxies ranging from dwarf irregulars to giant ellipticals. The data present a clear monotonic luminosity-velocity relation from 50 km/s to 500 km/s, with a bend below 80 km/s and a systematic offset between late- and early-type galaxies. For comparison to theory, we employ our LCDM "Bolshoi" simulation of dark matter, which has unprecedented mass and force resolution. We use halo abundance matching to assign rank-ordered galaxy luminosities to the dark matter halos. The resulting predictions for the luminosity-velocity relation are in excellent agreement with the available data on both early-type and late-type galaxies for the luminosity range from Mr = -14-22. We also compare our predictions for the "cold" baryon mass (i.e., stars and cold gas) of galaxies as a function of circular velocity with the available observations, again finding a very good agreement. The predicted circular velocity function is in agreement with the galaxy velocity function for 80-400 km/s. However, we find that the dark matter halos with Vcirc < 80 km/s are much more abundant than observed galaxies with the same Vcirc . We find that the two-point correlation function of galaxies in our model matches very well the results from the SDSS.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Haloes gone MAD: The Halo-Finder Comparison Project

[abridged] We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including friends-of-friends (FOF), spherical-overdensity (SO) and phase-space based algorithms. We further introduce a robust (and publicly available) suite of test scenarios that allows halo finder developers to compare the performance of their codes against those presented here. This set includes mock haloes containing various levels and distributions of substructure at a range of resolutions as well as a cosmological simulation of the large-scale structure of the universe. All the halo finding codes tested could successfully recover the spatial location of our mock haloes. They further returned lists of particles (potentially) belonging to the object that led to coinciding values for the maximum of the circular velocity profile and the radius where it is reached. All the finders based in configuration space struggled to recover substructure that was located close to the centre of the host halo and the radial dependence of the mass recovered varies from finder to finder. Those finders based in phase space could resolve central substructure although they found difficulties in accurately recovering its properties. Via a resolution study we found that most of the finders could not reliably recover substructure containing fewer than 30-40 particles. However, also here the phase space finders excelled by resolving substructure down to 10-20 particles. By comparing the halo finders using a high resolution cosmological volume we found that they agree remarkably well on fundamental properties of astrophysical significance (e.g. mass, position, velocity, and peak of the rotation curve).

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Mass Distribution and Assembly of the Milky Way from the Properties of the Magellanic Clouds

We present a new measurement of the mass of the Milky Way (MW) based on observed properties of its largest satellite galaxies, the Magellanic Clouds (MCs), and an assumed prior of a ΛCDM universe. The large, high-resolution Bolshoi cosmological simulation of this universe provides a means to statistically sample the dynamical properties of bright satellite galaxies in a large population of dark matter halos. The observed properties of the MCs, including their circular velocity, distance from the center of the MW, and velocity within the MW halo, are used to evaluate the likelihood that a given halo would have each or all of these properties; the posterior PDF for any property of the MW system can thus be constructed. This method provides a constraint on the MW virial mass, 1.2 +0.7 -0.3(stat.) +0.3 -0.4 (sys.) x 10^12 M\odot (68% confidence), which is consistent with recent determinations that involve very different assumptions. In addition, we calculate the posterior PDF for the density profile of the MW and its satellite accretion history. Although typical satellites of 10^12 M\odot halos are accreted over a wide range of epochs over the last 10 Gyr, we find a \sim72% probability that the Magellanic Clouds were accreted within the last Gyr, and a 50% probability that they were accreted together.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The MultiDark Database: Release of the Bolshoi and MultiDark Cosmological Simulations

We present the online MultiDark Database -- a Virtual Observatory-oriented, relational database for hosting various cosmological simulations. The data is accessible via an SQL (Structured Query Language) query interface, which also allows users to directly pose scientific questions, as shown in a number of examples in this paper. Further examples for the usage of the database are given in its extensive online documentation (www.multidark.org). The database is based on the same technology as the Millennium Database, a fact that will greatly facilitate the usage of both suites of cosmological simulations. The first release of the MultiDark Database hosts two 8.6 billion particle cosmological N-body simulations: the Bolshoi (250/h Mpc simulation box, 1/h kpc resolution) and MultiDark Run1 simulation (MDR1, or BigBolshoi, 1000/h Mpc simulation box, 7/h kpc resolution). The extraction methods for halos/subhalos from the raw simulation data, and how this data is structured in the database are explained in this paper. With the first data release, users get full access to halo/subhalo catalogs, various profiles of the halos at redshifts z=0-15, and raw dark matter data for one time-step of the Bolshoi and four time-steps of the MultiDark simulation. Later releases will also include galaxy mock catalogs and additional merging trees for both simulations as well as new large volume simulations with high resolution. This project is further proof of the viability to store and present complex data using relational database technology. We encourage other simulators to publish their results in a similar manner.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Large Scale Bias of Dark Matter Halos: Numerical Calibration and Model Tests

We measure the clustering of dark matter halos in a large set of collisionless cosmological simulations of the flat LCDM cosmology. Halos are identified using the spherical overdensity algorithm, which finds the mass around isolated peaks in the density field such that the mean density is Delta times the background. We calibrate fitting functions for the large scale bias that are adaptable to any value of Delta we examine. We find a ~6% scatter about our best fit bias relation. Our fitting functions couple to the halo mass functions of Tinker et. al. (2008) such that bias of all dark matter is normalized to unity. We demonstrate that the bias of massive, rare halos is higher than that predicted in the modified ellipsoidal collapse model of Sheth, Mo, & Tormen (2001), and approaches the predictions of the spherical collapse model for the rarest halos. Halo bias results based on friends-of-friends halos identified with linking length 0.2 are systematically lower than for halos with the canonical Delta=200 overdensity by ~10%. In contrast to our previous results on the mass function, we find that the universal bias function evolves very weakly with redshift, if at all. We use our numerical results, both for the mass function and the bias relation, to test the peak-background split model for halo bias. We find that the peak-background split achieves a reasonable agreement with the numerical results, but ~20% residuals remain, both at high and low masses.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Toward a halo mass function for precision cosmology: the limits of universality

We measure the mass function of dark matter halos in a large set of collisionless cosmological simulations of flat LCDM cosmology and investigate its evolution at z<~2. Halos are identified as isolated density peaks, and their masses are measured within a series of radii enclosing specific overdensities. We argue that these spherical overdensity masses are more directly linked to cluster observables than masses measured using the friends-of-friends algorithm (FOF), and are therefore preferable for accurate forecasts of halo abundances. Our simulation set allows us to calibrate the mass function at z=0 for virial masses in the range 10^{11} Msol/h < M < 10^{15} Msol/h, to <~ 5%. We derive fitting functions for the halo mass function in this mass range for a wide range of overdensities, both at z=0 and earlier epochs. In addition to these formulae, which improve on previous approximations by 10-20%, our main finding is that the mass function cannot be represented by a universal fitting function at this level of accuracy. The amplitude of the "universal" function decreases monotonically by ~20-50%, depending on the mass definition, from z=0 to 2.5. We also find evidence for redshift evolution in the overall shape of the mass function.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Properties of voids in the Local Volume

Current explanation of the overabundance of dark matter subhalos in the Local Group (LG) indicates that there maybe a limit on mass of a halo, which can host a galaxy. This idea can be tested using voids in the distribution of galaxies: at some level small voids should not contain any (even dwarf) galaxies. We use observational samples complete to M_B=-12 with distances less than 8 Mpc to construct the void function (VF): the distribution of sizes of voids empty of any galaxies. There are ~ 30 voids with sizes ranging from 1 to 5 Mpc. We also study the distribution of dark matter halos in very high resolution simulations of the LCDM model. The theoretical VF matches the observations remarkably well only if we use halos with circular velocities larger than 45 +/- 10 km/s. This agrees with the Local Group predictions. Small voids look quite similar to heir giant cousins: the density has a minimum at the center of a void and it increases as we get closer to the border. Thus, both the Local Group data and the nearby voids indicate that isolated halos below 45 +/- 10 km/s must not host galaxies and that small (few Mpc) voids are truly dark.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Observing the dark matter density profile of isolated galaxies

Using the Sloan Digital Sky Survey (SDSS), we probe the halo mass distribution by studying the velocities of satellites orbiting isolated galaxies. In a subsample that covers 2500 sq. degrees on the sky, we detect about 3000 satellites with absolute blue magnitudes going down to M_B = -14; most of the satellites have M_B=-16 to -18, comparable to the magnitudes of M32 and the Magellanic Clouds. After a careful, model-independent removal of interlopers, we find that the line-of-sight velocity dispersion of satellites declines with distance to the primary. For an L* galaxy the r.m.s. line-of-sight velocity changes from ~120 km/s at 20 kpc to ~60 km/s at 350 kpc. This decline agrees remarkably well with theoretical expectations, as all modern cosmological models predict that the density of dark matter in the peripheral parts of galaxies declines as rho_DM propto r^{-3}. Thus, for the first time we find direct observational evidence of the density decline predicted by cosmological models; we also note that this result contradicts alternative theories of gravity such as MOND. We also find that the velocity dispersion of satellites within 100 kpc scales with the absolute magnitude of the central galaxy as sigma propto L^{0.3}; this is very close to the Tully-Fisher relation for normal spiral galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Filament and Shape Statistics: A Quantitative Comparison of Cold + Hot and Cold Dark Matter Cosmologies vs. CfA1 Data

A new class of geometric statistics for analyzing galaxy catalogs is presented. Filament statistics quantify filamentarity and planarity in large scale structure in a manner consistent with catalog visualizations. These statistics are based on sequences of spatial links which follow local high-density structures. From these link sequences we compute the discrete curvature, planarity, and torsion. Filament statistics are applied to CDM and CHDM ($Ω_ν= 0.3$) simulations of Klypin \etal (1996), the CfA1-like mock redshift catalogs of Nolthenius, Klypin and Primack (1994, 1996), and the CfA1 catalog. We also apply the moment-based shape statistics developed by Babul \& Starkman (1992), Luo \& Vishniac (1995), and Robinson \& Albrecht (1996) to these same catalogs, and compare their robustness and discriminatory power versus filament statistics. For 100 Mpc periodic simulation boxes ($H_0 = 50$ km s$^{-1}$ Mpc$^{-1}$), we find discrimination of $\sim 4σ$ (where $σ$ represents resampling errors) between CHDM and CDM for selected filament statistics and shape statistics, including variations in the galaxy identification scheme. Comparing the CfA1 data versus the models does not yield a conclusively favored model; no model is excluded at more than a $\sim 2σ$ level for any statistic, not including cosmic variance which could further degrade the discriminatory power. We find that CfA1 discriminates between models poorly mainly due to its sparseness and small number of galaxies, not due to redshift distortion, magnitude limiting, or geometrical effects. We anticipate that the proliferation of large redshift surveys and simulations will enable the statistics presented here to provide robust discrimination between large-scale structure in various cosmological models.

0 citations0 reviewsNo code linkedOpen access
preprint1993arXiv

Structure Formation with Cold + Hot Dark Matter

We report results from high-resolution particle-mesh (PM) N-body simulations of structure formation in an $Ω=1$ cosmological model with a mixture of Cold plus Hot Dark Matter (C+HDM) having $Ω_{\rm cold}=0.6$, $Ω_ν=0.3$, and $Ω_{\rm baryon}=0.1$. We present analytic fits to the C+HDM power spectra for both cold and hot ($ν$) components, which provide initial conditions for our nonlinear simulations. In order to sample the neutrino velocities adequately, these simulations included six times as many neutrino particles as cold particles. Our simulation boxes were 14, 50, and 200~Mpc cubes (with $H_0=50$ km s$^{-1}$ Mpc$^{-1}$); we also did comparison simulations for Cold Dark Matter (CDM) in a 50~Mpc box. C+HDM with linear bias factor $b=1.5$ is consistent both with the COBE data and with the galaxy correlations we calculate. We find the number of halos as a function of mass and redshift in our simulations; our results for both CDM and C+HDM are well fit by a Press-Schechter model. The number density of galaxy-mass halos is smaller than for CDM, especially at redshift $z>2$, but the numbers of cluster-mass halos are comparable. We also find that on galaxy scales the neutrino velocities and flatter power spectrum in C+HDM result in galaxy pairwise velocities that are in good agreement with the data, and about 30\% smaller than in CDM with the same biasing factor. On scales of several tens of Mpc, the C+HDM streaming velocities are considerably larger than CDM. Thus C+HDM looks promising as a model of structure formation.

0 citations0 reviewsNo code linkedOpen access