Search

In their recent study, Labb\'e et al. used multi-band infrared images captured by the James Webb Space Telescope (JWST) to discover a population of red massive galaxies that formed approximately 600 million years after the Big Bang. The authors reported an extraordinarily large density of these galaxies, with stellar masses exceeding $10^{10}$ solar masses, which, if confirmed, challenges the standard cosmological model as suggested by recent studies. However, this conclusion is disputed. We contend that during the early epochs of the universe the stellar mass-to-light ratio could not have reached the values reported by Labb\'e et al. A model of galaxy formation based on standard cosmology provides support for this hypothesis, predicting the formation of massive galaxies with higher ultraviolet (UV) luminosity, which produce several hundred solar masses of stars per year and containing significant dust. These forecasts are consistent with the abundance of JWST/HST galaxies selected photometrically in the rest-frame UV wavelengths and with the properties of the recent spectroscopically-confirmed JWST/HST galaxies formed during that era. Discrepancies with Labb\'e et al. may arise from overestimation of the stellar masses, systematic uncertainties, absence of JWST/MIRI data, heavy dust extinction affecting UV luminosities, or misidentification of faint red AGN galaxies at closer redshifts. The current JWST/HST results, combined with a realistic galaxy formation model, provide strong confirmation of the standard cosmology., Comment: Submitted to Nature, matters arising

Measurements of the luminosity function of active galactic nuclei (AGN) at high redshift ($z\gtrsim 6$) are expected to suffer from field-to-field variance, including cosmic and Poisson variances. Future surveys, such as those from the Euclid telescope and James Webb Space Telescope (JWST), will also be affected by field variance. We use the Uchuu simulation, a state-of-the-art cosmological $N$-body simulation with 2.1 trillion particles in a volume of $25.7~\mathrm{Gpc}^3$, that has sufficient mass resolution to resolve dwarf-size systems, combined with a semi-analytic galaxy and AGN formation model, to generate the Uchuu-$\nu^2$GC catalog, publicly available, that allows us to investigate the field-to-field variance of the luminosity function of AGN. With this Uchuu-$\nu^2$GC model, we quantify the cosmic variance as a function of survey area, AGN luminosity, and redshift. In general, cosmic variance decreases with increasing survey area and decreasing redshift. We find that at $z\sim6-7$, the cosmic variance depends weakly on AGN luminosity, in particular for small survey areas (0.01 and 0.1 deg$^2$). This is because the typical mass of dark matter haloes in which AGN reside does not significantly depend on luminosity. Due to the rarity of AGN, Poisson variance dominates the total field-to-field variance, especially for bright AGN. We discuss uncertainties present in the estimation of the faint-end of the AGN luminosity function from recent observations, and extend this to make predictions for the expected number of AGN and their variance for upcoming observations with Euclid, JWST, and the Legacy Survey of Space and Time (LSST). In particular, we predict that the Euclid deep survey will find 120--240 (16--80) AGN -- depending on the model -- with rest-frame UV absolute magnitude brighter than $-20$ ($-20.5$) at $z=6.3$ ($z=7$) in the Euclid H-band deep survey (abridged)., Comment: 17 pages, 17 figures, submitted to MNRAS, corrected the following link, All Uchuu-$\nu^2$GC data used for this paper is publicly available in our Skies & Universes website http://skiesanduniverses.org/Simulations/Uchuu/

The cluster mass-richness relation (MRR) is an observationally efficient and potentially powerful cosmological tool for constraining the mean matter density of the universe and the amplitude of fluctuations using the cluster abundance technique. We derive the MRR relation using GalWCat19, a publicly available galaxy cluster catalog we created from the Sloan Digital Sky Survey-DR13 spectroscopic dataset. The MRR shows a tail at the low-richness end. Using the Illustris-TNG and mini-Uchuu cosmological numerical simulations, we demonstrate that this tail is caused by systematical uncertainties. We show that, by means of a judicious cut, identified by the use of the Hinge function, it is possible to determine a richness threshold above which the MRR is linear i.e., where cluster mass scales with richness as logM_200 = alpha + beta logN_200. We derive the MRR and show it is consistent with both sets of simulations with a slope of beta ~ 1. We use our MRR to estimate cluster masses from the GalWCat19 catalog which we then use to set constraints on omega_m and sigma_8. Utilizing the all-member MRR, we obtain constraints of omega_m = 0.31 (+0.04-0.03) and sigma_8 = 0.82 (+0.05-0.04), and utilizing the red-member MRR, we obtain omega_m = 0.31 (+0.04-0.03) and sigma_8 = 0.81 (+0.05-0.04). Our constraints on omega_m and sigma_8 are consistent and very competitive with the Planck 2018 results.

We derive cosmological constraints on the matter density, \om, and the amplitude of fluctuations, \sig, using $\mathtt{GalWCat19}$, a catalog of 1800 galaxy clusters we identified in the Sloan Digital Sky Survey-DR13 spectroscopic data set using our GalWeight technique to determine cluster membership \citep{Abdullah18,Abdullah19}. By analyzing a subsample of 756 clusters in a redshift range of $0.045\leq z \leq 0.125$ and virial masses of $M\geq 0.8\times10^{14}$ \hm ~with mean redshift of $z = 0.085$, we obtain \om ~$=0.310^{+0.023}_{-0.027} \pm 0.041$ (systematic) and \sig ~$=0.810^{+0.031}_{-0.036}\pm 0.035$ (systematic), with a cluster normalization relation of $\sigma_8= 0.43 \Omega_m^{-0.55}$. There are several unique aspects to our approach: we use the largest spectroscopic data set currently available, and we assign membership using the GalWeight technique which we have shown to be very effective at simultaneously maximizing the number of {\it{bona fide}} cluster members while minimizing the number of contaminating interlopers. Moreover, rather than employing scaling relations, we calculate cluster masses individually using the virial mass estimator. Since $\mathtt{GalWCat19}$ is a low-redshift cluster catalog we do not need to make any assumptions about evolution either in cosmological parameters or in the properties of the clusters themselves. Our constraints on \om ~and \sig ~are consistent and very competitive with those obtained from non-cluster abundance cosmological probes such as Cosmic Microwave Background (CMB), Baryonic Acoustic Oscillation (BAO), and supernovae (SNe). The joint analysis of our cluster data with Planck18+BAO+Pantheon gives \om ~$=0.315^{+0.013}_{-0.011}$ and \sig ~$=0.810^{+0.011}_{-0.010}$., Comment: 14 pages, 10 figures, accepted for publication in the ApJ

We derive cosmological constraints on the matter density, \om, and the amplitude of fluctuations, \sig, using $\mathtt{GalWCat19}$, a catalog of 1800 galaxy clusters we identified in the Sloan Digital Sky Survey-DR13 spectroscopic data set using our GalWeight technique to determine cluster membership \citep{Abdullah18,Abdullah19}. By analyzing a subsample of 756 clusters in a redshift range of $0.045\leq z \leq 0.125$ and virial masses of $M\geq 0.8\times10^{14}$ \hm ~with mean redshift of $z = 0.085$, we obtain \om ~$=0.310^{+0.023}_{-0.027} \pm 0.041$ (systematic) and \sig ~$=0.810^{+0.031}_{-0.036}\pm 0.035$ (systematic), with a cluster normalization relation of $��_8= 0.43 ��_m^{-0.55}$. There are several unique aspects to our approach: we use the largest spectroscopic data set currently available, and we assign membership using the GalWeight technique which we have shown to be very effective at simultaneously maximizing the number of {\it{bona fide}} cluster members while minimizing the number of contaminating interlopers. Moreover, rather than employing scaling relations, we calculate cluster masses individually using the virial mass estimator. Since $\mathtt{GalWCat19}$ is a low-redshift cluster catalog we do not need to make any assumptions about evolution either in cosmological parameters or in the properties of the clusters themselves. Our constraints on \om ~and \sig ~are consistent and very competitive with those obtained from non-cluster abundance cosmological probes such as Cosmic Microwave Background (CMB), Baryonic Acoustic Oscillation (BAO), and supernovae (SNe). The joint analysis of our cluster data with Planck18+BAO+Pantheon gives \om ~$=0.315^{+0.013}_{-0.011}$ and \sig ~$=0.810^{+0.011}_{-0.010}$., 14 pages, 10 figures, accepted for publication in the ApJ

We present the first galaxy-OVI absorption kinematic study for 20 absorption systems (EW>0.1~Å) associated with isolated galaxies (0.15$33^{\circ}$). OVI absorption along the galaxy major axis is not correlated with galaxy rotation kinematics, with only 1/10 systems that could be explained with rotation/accretion models. This is in contrast to co-rotation commonly observed for MgII absorption. OVI along the minor axis could be modeled by accelerating outflows but only for small opening angles, while the majority of the OVI is decelerating. Along both axes, stacked OVI profiles reside at the galaxy systemic velocity with the absorption kinematics spanning the entire dynamical range of their galaxies. The OVI found in AMR cosmological simulations exists within filaments and in halos of ~50 kpc surrounding galaxies. Simulations show that major axis OVI gas inflows along filaments and decelerates as it approaches the galaxy while increasing in its level of co-rotation. Minor axis outflows in the simulations are effective within 50-75 kpc beyond that they decelerate and fall back onto the galaxy. Although the simulations show clear OVI kinematic signatures they are not directly comparable to observations. When we compare kinematic signatures integrated through the entire simulated galaxy halo we find that these signatures are washed out due to full velocity distribution of OVI throughout the halo. We conclude that OVI alone does not serve as a useful kinematic indicator of gas accretion, outflows or star-formation and likely best probes the halo virial temperature., 24 pages, 21 figures, 4 tables. Accepted to ApJ on November 14, 2018

We analyze the spectra of 300,000 luminous red galaxies (LRGs) with stellar masses $M_* \gtrsim 10^{11} M_{\odot}$ from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). By studying their star-formation histories, we find two main evolutionary paths converging into the same quiescent galaxy population at $z\sim0.55$. Fast-growing LRGs assemble $80\%$ of their stellar mass very early on ($z\sim5$), whereas slow-growing LRGs reach the same evolutionary state at $z\sim1.5$. Further investigation reveals that their clustering properties on scales of $\sim$1-30 Mpc are, at a high level of significance, also different. Fast-growing LRGs are found to be more strongly clustered and reside in overall denser large-scale structure environments than slow-growing systems, for a given stellar-mass threshold. Our results imply a dependence of clustering on stellar-mass assembly history (naturally connected to the mass-formation history of the corresponding halos) for a homogeneous population of similar mass and color, which constitutes a strong observational evidence of galaxy assembly bias., Comment: 6 pages, 5 figures. Submitted to ApJ Letters

Numerical simulations play a key important role in modern cosmology. Examples are plenty including the cosmic web - large scale structure of the distribution of galaxies in space - which was first observed in N-body simulations and later discovered in observations. The cuspy dark matter halo profiles, the overabundance of satellites, the Too-Big-Too-Fail problem are other examples of theoretical predictions that have a dramatic impact on recent developments in cosmology and galaxy formation. Large observational surveys such as e.g. SDSS, Euclid, and LSST are intimately connected with extensive cosmological simulations that provide statistical errors and tests for systematics. Accurate predictions for baryonic acoustic oscillations and redshift space distortions from high-resolution and large-volume cosmological simulations are required for interpretation of these large-scale galaxy/qso surveys. However, most of the results from extensive computer simulations, that would be greatly beneficial if publicly available, are still in hands of few research groups. Even when the simulation data is available, sharing vast amounts of data can be overwhelming. We argue that there is an effective and simple path to expand the data access and dissemination of numerous results from different cosmological models. Here we demonstrate that public access can be effectively provided with relatively modest resources. Among different results, we release for the astronomical community terabytes of raw data of th popular Bolshoi and MultiDark simulations. We also provide numerous results that are focused on mimicking observational data and galaxy surveys for major projects. Skies and Universes is a community effort: data are produced and shared by many research groups. We offer to other cosmologists and astronomers to host their data products in the skiesanduniverses.org space., Comment: 4 pages

International audience; 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 (ELGs) to probe cosmological models by mapping the large-scale structure of the Universe in the redshift range 0.6 \textless z \textless 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 \textless z \textless 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 on to 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 similar to 0.8 live in haloes of (1 +/- 0.5) x 10(12) h(-1)M(circle dot) and 22.5 +/- 2.5 per cent 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.

We study the galaxy clustering dependence on the $\left[\mathrm{O\scriptsize{II}}\right]$ emission line luminosity in the SDSS DR7 Main galaxy sample at mean redshift $z\sim0.1$. We select volume-limited samples of galaxies with different $\left[\mathrm{O\scriptsize{II}}\right]$ luminosity thresholds and measure their projected, monopole and quadrupole two-point correlation functions. We model these observations using the 1$h^{-1}\rm{Gpc}$ MultiDark Planck cosmological simulation and generate light-cones with the SUrvey GenerAtoR algorithm. To interpret our results, we adopt a modified (Sub)Halo Abundance Matching scheme, accounting for the stellar mass incompleteness of the emission line galaxies. The satellite fraction constitutes an extra parameter in this model and allows to optimize the clustering fit on both small and intermediate scales (i.e. $r_p\lesssim 30h^{-1}\rm{Mpc})$, with no need of any velocity bias correction. We find that, in the local Universe, the $\left[\mathrm{O\scriptsize{II}}\right]$ luminosity correlates with all the clustering statistics explored and with the galaxy bias. This latter quantity correlates more strongly with the SDSS $r$-band magnitude than $\left[\mathrm{O\scriptsize{II}}\right]$ luminosity. In conclusion, we propose a straightforward method to produce reliable clustering models, entirely built on the simulation products, which provides robust predictions of the typical ELG host halo masses and satellite fraction values. The SDSS galaxy data, MultiDark mock catalogues and clustering results are made publicly available., Comment: 9 pages, 7 figures, 1 table

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., 11 pages, 7 figures, 2 tables

This article has been accepted for publication in Monthly Notices of the Royal Astronomical Society © 2011 RAS © 2011 The Authors Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved., We make detailed theoretical predictions for the assembly properties of the Local Group (LG) in the standard Λ cold dark matter cosmological model. We use three cosmological N-body dark matter simulations from the Constrained Local Universe Simulations project, which are designed to reproduce the main dynamical features of the matter distribution down to the scale of a few Mpc around the LG. Additionally, we use the results of an unconstrained simulation with a 60 times larger volume to calibrate the influence of cosmic variance. We characterize the mass aggregation history (MAH) for each halo by three characteristic times: the formation, assembly and last major merger times. A major merger is defined by a minimal mass ratio of 10: 1. We find that the three LGs share a similar MAH with formation and last major merger epochs placed on average ≈10-12 Gyr ago. Between 12 and 17 per cent of the haloes in the mass range 5 × 10 11 < M h < 5 × 10 12h -1M ⊙ have a similar MAH. In a set of pairs of haloes within the same mass range, a fraction of 1-3 per cent share similar formation properties as both haloes in the simulated LG. An unsolved question posed by our results is the dynamical origin of the MAH of the LGs. The isolation criteria commonly used to define LG-like haloes in unconstrained simulations do not narrow down the halo population into a set with quiet MAHs, nor does a further constraint to reside in a low-density environment. The quiet MAH of the LGs provides a favourable environment for the formation of disc galaxies like the Milky Way and M31. The timing for the beginning of the last major merger in the Milky Way dark matter halo matches with the gas-rich merger origin for the thick component in the galactic disc. Our results support the view that the specific large- and mid-scale environments around the LG play a critical role in shaping its MAH and hence its baryonic structure at present, GY would like to thank the MICINN (Spain) for financial support under project numbers FPA 2009-08958, AYA 2009-13875-C03 and the SyeC Consolider project CSD 2007-0050

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., 28 pages, 9 figures, submitted to New Astronomy

The condensation of gas and stars in the inner regions of dark matter halos leads to a more concentrated dark matter distribution. While this effect is based on simple gravitational physics, the question of its validity in hierarchical galaxy formation has led to an active debate in the literature. We use a collection of several state-of-the-art cosmological hydrodynamic simulations to study the halo contraction effect in systems ranging from dwarf galaxies to clusters of galaxies, at high and low redshift. The simulations are run by different groups with different codes and include hierarchical merging, gas cooling, star formation, and stellar feedback. We show that in all our cases the inner dark matter density increases relative to the matching simulation without baryon dissipation, at least by a factor of several. The strength of the contraction effect varies from system to system and cannot be reduced to a simple prescription. We present a revised analytical model that describes the contracted mass profile to an rms accuracy of about 10%. The model can be used to effectively bracket the response of the dark matter halo to baryon dissipation. The halo contraction effect is real and must be included in modeling of the mass distribution of galaxies and galaxy clusters., 11 pages, submitted to ApJ

We study the concentration of dark matter halos and its evolution in N-body simulations of the standard LCDM cosmology. The results presented in this paper are based on 4 large N-body simulations with about 10 billion particles each: the Millennium-I and II, Bolshoi, and MultiDark simulations. The MultiDark (or BigBolshoi) simulation is introduced in this paper. This suite of simulations with high mass resolution over a large volume allows us to compute with unprecedented accuracy the concentration over a large range of scales (about six orders of magnitude in mass), which constitutes the state-of-the-art of our current knowledge on this basic property of dark matter halos in the LCDM cosmology. We find that there is consistency among the different simulation data sets. We confirm a novel feature for halo concentrations at high redshifts: a flattening and upturn with increasing mass. The concentration c(M,z) as a function of mass and the redshift and for different cosmological parameters shows a remarkably complex pattern. However, when expressed in terms of the linear rms fluctuation of the density field sigma(M,z), the halo concentration c(sigma) shows a nearly-universal simple U-shaped behaviour with a minimum at a well defined scale at sigma=0.71. Yet, some small dependences with redshift and cosmology still remain. At the high-mass end (sigma < 1) the median halo kinematic profiles show large signatures of infall and highly radial orbits. This c-sigma(M,z) relation can be accurately parametrized and provides an analytical model for the dependence of concentration on halo mass. When applied to galaxy clusters, our estimates of concentrations are substantially larger -- by a factor up to 1.5 -- than previous results from smaller simulations, and are in much better agreement with results of observations. (abridged), Submitted to MNRAS

This article has been accepted for publication in Monthly Notices of the Royal Astronomical Society © 2011 RAS © 2011 The authors Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved, 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, spherical-overdensity and phase-space-based algorithms. We further introduce a robust (and publicly available) suite of test scenarios that allow 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. Through 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). We further suggest to utilize the peak of the rotation curve, vmax, as a proxy for mass, given the arbitrariness in defining a proper halo edge, We are greatly indebted to the ASTROSIM network of the European Science Foundation (Science Meeting 2910) for financially supporting the workshop ‘Haloes going MAD’ held in Miraflores de la Sierra near Madrid in 2010 May where all of this work was initiated.AK is supported by the Spanish Ministerio de Ciencia e Innovación (MICINN) in Spain through the Ramón y Cajal programme as well as the grants AYA 2009-13875-C03-02, AYA2009- 12792-C03-03, CSD2009-00064 and CAM S2009/ESP-1496. SRK acknowledges the support by the MICINN under the Consolider-Ingenio, SyeC project CSD-2007-00050. SP and VQ have also been supported by the MICINN (grants AYA2010-21322-C03-01 and CONSOLIDER2007-00050) and the Generalitat Valenciana (grant PROMETEO-2009-103). MZ is supported by NSF grant AST-0708087. GY acknowledges financial support from MICINN (Spain) under project AYA 2009-13875-C03-02 and the ASTROMADRID project S2009/ESP-1496 financed by Comunidad de Madrid

We describe the basic ideas of MPI parallelization of the N-body Adaptive Refinement Tree (ART) code. The code uses self-adaptive domain decomposition where boundaries of the domains (parallelepipeds) constantly move -- with many degrees of freedom -- in the search of the minimum of CPU time. The actual CPU time spent by each MPI task on previous time-step is used to adjust boundaries for the next time-step. For a typical decomposition of 5^3 domains, the number of possible changes in boundaries is 3^{84}. We describe two algorithms of finding minimum of CPU time for configurations with a large number of domains. Each MPI task in our code solves the N-body problem where the large-scale distribution of matter outside of the boundaries of a domain is represented by relatively few temporary large particles created by other domains. At the beginning of a zero-level time-step, domains create and exchange large particles. Then each domain advances all its particles for many small time-steps. At the end of the large step, the domains decide where to place new boundaries and re-distribute particles. The scheme requires little communications between processors and is very efficient for large cosmological simulations., 15 pages, 5 figures, conference "High Performance Computing in Science and Engineering Garching/Munich 2007", Eds. S. Wagner et al, Springer, Berlin 2008

We present an approximation for the average density profile of dark matter haloes in the LambdaCDM cosmological model, which is accurate to within 10--15% even for large radial distances from 0.05R_{vir} up to 10R_{vir} for halo masses ranging from 10^{11.5} to 10^{15.0} h^{-1}M_{Sun}. We propose a modified form of the Navarro, Frenk & White (NFW) approximation: \rho(r)=\rho_{NFW}(r)+ A (r/R_{vir})^{-1}+B (1+r/R_{vir})^{-1}. This generalized expression, which is applicable to the external regions of dark matter haloes, only very slightly affects the density in the inner regions of haloes. The strong correlation among the different parameters in the model allows us to describe the profile in terms of just one parameter: the virial mass. We integrate our density profile to derive the enclosed mass in a sphere of a given radius and compare it with the NFW results. We find that the NFW underestimates the enclosed mass by more than 50% at 10R_{vir}, whereas our model reproduces the results from numerical simulations to within 2% accuracy even at this distance. We also use this new approximation to study the weak gravitational lensing and to obtain an analytic expression for the tangential shear. This allows us to quantify the contribution to the shear from the outer regions of the density profile. For the first time we calculate the difference between the tangential shear calculated via the NFW profile and the corresponding result when the external regions of haloes in cosmological simulations are taken into account. We find a 4% difference for all the mass ranges under study., Comment: 12 pages, 10 figures, submitted for publication in MNRAS

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., 5 pages 1 figure. To appear in proceedings of the conference "Galaxies in the Local Volume", Sydney, 8 to 13 July 2007

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 then 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. There are smaller halos in the voids, but they should not produce any luminous matter. Small voids look quite similar to their giant cousins: the density has a minimum at the center of a void and it increases as we get closer to the border. Small nonluminous halos inside the void form a web of tiny filaments. 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., Comment: 5 pages, 1 figure

We present a statistical analysis of voids in the 2dF galaxy redshift survey (2dFGRS). In order to detect the voids, we have developed two robust algorithms. We define voids as non-overlapping maximal spheres empty of halos or galaxies with mass or luminosity above a given one. We search for voids in cosmological $N$-Body simulations to test the performance of our void finders. We obtain and analyze the void statistics for several volume-limited samples for the North Galactic Strip (NGP) and the South Galactic Strip (SGP) constructed from the 2dFGRS full data release. We find that the results obtained from the NGP and the SGP are statistically compatible. From the results of several statistical tests we conclude that voids are essentially uncorrelated, with at most a mild anticorrelation and that there is a dependence of the void number density on redshift at least at the 99.5% confidence level. We develop a technique to correct the distortion caused by the fact that we use the redshift as the radial coordinate. We calibrate this technique with mock catalogues and find that the correction might be of some relevance to carry out accurate inferences from void statistics. We study the statistics of the galaxies inside nine nearby voids. We find that galaxies in voids are not randomly distributed: they form structures like filaments. We also obtain the galaxy number density profile in voids. This profile follow a similar but steeper trend to that follow by halos in voids., 16 pages, 10 figures, replaced to match the published version in MNRAS

The dynamics of the inner disk of many galaxies is dominated by bars. As a result, resonances between the bar and the disk become an important factor. In order to detect resonances, we measure angular and radial frequencies of individual orbits in N-body models of barred galaxies. In our simulations with several million particles, stellar disks are immersed in live halos of dark matter. A bar forms spontaneously and interacts with the disk through resonances. We detect 9 of them. The corotation (CR) and the inner Lindblad resonance (ILR) are the most significant. The spatial distributions of particles trapped at these resonances show that ILR is not localized at a given radius. Particles trapped at ILR lie along the bar spanning a wide range of radii. On the other hand, particles at the corotation resonance form a ring around the corotation radius. We study how orbits evolve near corotation and how they are trapped. When an orbit is getting trapped, it evolves strongly, but after that there is little evolution. Those trapped orbits circulate along the corotation ring or librate around a stable lagragian point. The same technique is applied to study resonances in the halo. The bar interacts strongly with the halo particles that stay close to the disk. We find the corotation resonance and the inner Lindblad resonance for halo particles., Comment: This paper contains the figures presented in the conference: "Dynamics of Galaxies: Baryons and Dark Matter'' (Las Vegas, NV, March 2005). The full talk is available in http://astronomy.nmsu.edu/dept/html/directory.grads.danielcv.shtml

If dark matter (DM), which is considered to constitute most of the mass of galaxies, is made of supersymmetric (SUSY) particles, the centers of galaxies should emit gamma-rays produced by their self-annihilation. We present accurate estimates of continuum gamma-ray fluxes due to neutralino annihilation in the central regions of the Milky Way. We use detailed models of our Galaxy, which satisfy available observational data, and include some important physical processes, which were previously neglected. Our models predict that spatially extended annihilation signal should be detected at high confidence levels by incoming experiments assuming that neutralinos make up most of the DM in the Universe and that they annihilate according to current SUSY models., 4 pages, submitted to Physical Review Letters

The theory of structure formation predicts that galaxies form within extended massive dark matter halos built from smaller pieces that collided and merged, resulting in the hierarchy of galaxies, groups, and clusters observed today. Here we present constrained numerical simulations designed to match the observed local universe as well as possible and to study the formation, evolution and present day properties of such dark matter halos in different environments. Simulations have been done with different versions of the Adaptive Refinement Tree (ART) N-body code, some of them including eulerian gasdynamics., Comment: 10 pages, 6 figures, to appear in the proceedings of the "First Joint HLRB and KONWIHR Result and Reviewing Workshop" (Munich, October 2002). Version with high resolution figures available at http://www.aip.de/People/SGottloeber/papers/h009za.ps.gz (gzipped 5.3 Mb)

[Abridged] The interaction rates of dark-matter halos and subhalos, are computed using high-resolution cosmological N-body simulations of the Lambda-CDM model. Although the number fraction of subhalos of mass $>2x10^{11}\hsolmass$ is only ~10%, we find that the interaction rate of such subhalos is relatively high because they reside in high density environments. At low redshift, the subhalo collisions dominate the total collision rate, and even at z=3 they are involved in more than 30% of all collisions. About 40% of the "major" mergers are between subhalos. Therefore subhalo interactions must be incorporated in models of structure formation. We find that the collision rate between halos in physical density units, is $\propto (1+z)^��$, with $��= 3-4$, in agreement with earlier simulations and most observational data. We test previous analytic estimates of the interaction rates of subhalos based on statistical models, which could be very inaccurate because of the small number of subhalos and the variation of conditions within small host halos. We find that, while such statistical estimates may severely overestimate the rate within hosts $< 10^{13} \hsolmass$, typical of high redshifts, they are valid for larger hosts regardless of the number of subhalos in them. We find the Makino & Hut (1997) estimate of the subhalo merger rate to be valid for hosts $ > 6x10^{11}\hsolmass$ at all redshifts. The collision rate between subhalos and the central object of their host halo is approximated relatively well using the timescale for dynamical friction in circular orbits. This approximation fails in ~40% of the cases, partly because of deviations from the assumption of circular orbits and partly because of the invalidity of the assumption that the subhalo mass is negligible., 18 pages, submitted to MNRAS

We give a summary of recent results on spatial and velocity biases in cosmological models. Progress in numerical techniques made it possible to simulate halos in large volumes with a such accuracy that halos survive in dense environments of groups and clusters of galaxies. Halos in simulations look like real galaxies, and, thus, can be used to study the biases -- differences between galaxies and the dark matter. The biases depend on scale, redshift, and circular velocities of selected halos. Two processes seem to define the evolution of the spatial bias: (1) statistical bias and (2) merger bias (merging of galaxies, which happens preferentially in groups, reduces the number of galaxies, but does not affect the clustering of the dark matter). There are two kinds of velocity bias. The pair-wise velocity bias is b_12=0.6-0.8 at r< 5Mpc/h, z=0. This bias mostly reflects the spatial bias and provides almost no information on the relative velocities of the galaxies and the dark matter. One-point velocity bias is a better measure of the velocities. Inside clusters the galaxies should move slightly faster (b_v =1.1-1.3) than the dark matter. Qualitatively this result can be understood using the Jeans equations of the stellar dynamics. For the standard LCDM model we find that the correlation function and the power spectrum of galaxy-size halos at z=0 are antibiased on scales r, Lecture at the Summer School "Relativistic Cosmology: Theory and Observations", Italy, Como, May 2000. AASTeX and natbib

This lecture gives a short introduction to different methods used in cosmology. The focus is on major features of N-body simulations: equations, main numerical techniques, effects of resolution, and methods of halo identification., Lecture at the Summer School "Relativistic Cosmology: Theory and Observations", Italy, Como, May 2000. AASTeX and natbib

We have simulated the evolution of a large sample of X-ray clusters in a Lambda-CDM universe at high spatial resolution using adaptive mesh refinement and placed the results in an online archive for public access. The Simulated Cluster Archive website http://sca.ncsa.uiuc.edu provides tools for interactive 2D and 3D analysis of gas and dark matter fields, X-ray and SZ imaging, and data export. We encourage community use and solicit their feedback., 2 pages, 0 figures. Poster paper presented at "Large-Scale Structure in the X-Ray Universe", Santorini, Greece 20-22 September 1999. Eds. M. Plionis & I. Georgantopoulos, (Atlantisciences, Athens). Visit the Simulated Cluster Archive at http://sca.ncsa.uiuc.edu

[Abridged] We analyze the exchange of dark matter between halos, subhalos, and their environments in a high-resolution cosmological N-body simulation of a Lambda CDM cosmology. At each analyzed redshift z we divide the dark matter particles into 4 components: (i) isolated galactic halos, (ii) subhalos, (iii) the diffuse medium of group and cluster halos, and (iv) the background outside of virialized halos. We follow the time evolution of the mass distribution and flows between these components and provide fitting functions for the exchange rates. We use our derived exchange rates to gauge the importance of metal redistribution in the universe due solely to gravity-induced interactions. The diffuse metallicity in clusters is predicted to be ~40% that in isolated galaxies (~55% of groups) at z=0, and should be lower only slightly by z=1, consistent with observations. The metallicity of the diffuse media in poor groups is expected to be lower by a factor of 5 by z~2, in agreement with the observed metallicity of damped Ly$\alpha$ systems. The metallicity of the background IGM is predicted to be (1-3)x10^{-4} that of z=0 clusters, also consistent with observations. The agreement of predicted and observed trends indicates that gravitational interaction alone may play an important role in metal enrichment of the intra-cluster and intergalactic media., Comment: 12 pages, submitted to MNRAS

Properties of dark matter halos are reviewed. Taken from different publications, we present results on (1) the mass and velocity functions, (2) density and velocity profiles, and (3) concentration of halos. In the range of radii r=(0.005-1)rvir the density profile for a quiet isolated halo is very accurately approximated by a fit suggested by Moore etal (1997): rho=1/x^1.5(1+x^1.5), where x=r/rs and rs is a characteristic radius. The fit suggested by Navarro et al (1995) rho= 1/x(1+x)^2, also gives a very satisfactory approximation with relative errors of about 10% for radii not smaller than 1% of the virial radius. The mass function of z=0 halos with mass below 10^{13}Msun/h is approximated by a power law with slope alpha =-1.85. The slope increases with the redshift. The velocity function of halos with Vmax< 500km/s is also a power law with the slope beta= -3.8-4. The power-law extends to halos at least down to 10km/s. It is also valid for halos inside larger virialized halos. The concentration of halos depends on mass (more massive halos are less concentrated) and environment, with isolated halos being less concentrated than halos of the same mass inside clusters. Halos have intrinsic scatter of concentration: at 1sigma level halos with the same mass have Delta(log{c_vir})=0.18 or, equivalently, DeltaVmax/Vmax =0.12. Velocity anisotropy for both subhalos and the dark matter is approximated by beta(r) =0.15+2x/[x^2+4], where x is radius in units of the virial radius., Lecture at the Summer School "Relativistic Cosmology: Theory and Observations", Italy, Como, May 2000. AASTeX and natbib

We study the distribution functions of mass and circular velocity for dark matter halos in N-body simulations of the $\Lambda$CDM cosmology, addressing redshift and environmental dependence. The dynamical range enables us to resolve subhalos and distinguish them from "distinct" halos. The mass function is compared to analytic models, and is used to derive the more observationally relevant circular velocity function. The distribution functions in the velocity range 100--500 km/s are well fit by a power-law with two parameters, slope and amplitude. We present the parameter dependence on redshift and provide useful fitting formulae. The amplitudes of the mass functions decrease with z, but, contrary to naive expectation, the comoving density of halos of a fixed velocity ~200 km/s actually increases out to z=5. This is because high-z halos are denser, so a fixed velocity corresponds to a smaller mass. The slope of the velocity function at z=0 is as steep as ~ -4, and the mass and velocity functions of distinct halos steepen with increasing z, while the functions of subhalos do not steepen with $z$, and become even flatter at z>2. A simple observable prediction is that the slope of the velocity function of isolated galaxies is steeper than that of galaxies in groups by as much as unity, reflecting the density biasing of high-velocity halos. We confirm that the Press-Schechter approximation typically overestimates the halo mass function by a factor of ~2, while modified approximations provide improved predictions., Comment: 26 pages, 9 figures. submitted to MNRAS

We give a summary of our recent studies of spatial and velocity biases of galaxy-size halos in cosmological models. Recent progress in numerical techniques made it possible to simulate halos in large volumes with a such accuracy that halos survive in dense environments of groups and clusters of galaxies. Dark matter halos in simulations look like real galaxies, and, thus, can be used to study the biases - differences between galaxies and the dark matter. For the standard LCDM model we find that the correlation function and the power spectrum of galaxy-size halos at z=0 are antibiased on scales r, 8 pages, 6 figures, Proceedings of "Cosmic Flows", Victoria, UVic, July 1999, paspconf and epsf

We discuss the impact of the cosmological environment on the evolution of dark matter halos using a high-resolution simulation within a spatially flat LCDM cosmology., 4 pages, 1 figure, contribution to the Morphology Conference in Johannesburg, September 1999

We use N-body simulations to study the velocity bias of dark matter halos, the difference in the velocity fields of dark matter and halos, in a flat low- density LCDM model. The high force, 2kpc/h, and mass, 10^9Msun/h, resolution allows dark matter halos to survive in very dense environments of groups and clusters making it possible to use halos as galaxy tracers. We find that the velocity bias pvb measured as a ratio of pairwise velocities of the halos to that of the dark matter evolves with time and depends on scale. At high redshifts (z ~5) halos move generally faster than the dark matter almost on all scales: pvb(r)~1.2, r>0.5Mpc/h. At later moments the bias decreases and gets below unity on scales less than r=5Mpc/h: pvb(r)~(0.6-0.8) at z=0. We find that the evolution of the pairwise velocity bias follows and probably is defined by the spatial antibias of the dark matter halos at small scales. One-point velocity bias b_v, defined as the ratio of the rms velocities of halos and dark matter, provides a more direct measure of the difference in velocities because it is less sensitive to the spatial bias. We analyze b_v in clusters of galaxies and find that halos are ``hotter'' than the dark matter: b_v=(1.2-1.3) for r=(0.2-0.8)r_vir, where r_vir is the virial radius. At larger radii, b_v decreases and approaches unity at r=(1-2)r_vir. We argue that dynamical friction may be responsible for this small positive velocity bias b_v>1 found in the central parts of clusters. We do not find significant difference in the velocity anisotropy of halos and the dark matter. The dark matter the velocity anisotropy can be approximated as beta(x)=0.15 +2x/(x^2+4), where x is measured in units of the virial radius., 13 pages, Latex, AASTeXv5 and natbib

A significant fraction of mass in the universe is believed to be in the form of dark matter (DM). Due to gravitational instability, the DM collapses hierarchically into DM halos. In this contribution we present a study of the formation and evolution of such DM halos in a COBE-normalized spatially flat LCDM model (Omega_0=1-Omega_L=0.3; h=0.7) using high-resolution N-body simulations. The novelty of this study is use of the newly developed halo-finding algorithms to study the evolution of both isolated and satellite (located inside virial radii of larger group- and cluster-size systems) halos. The force and mass resolution required for a simulated halo to survive in the high-density environments typical of groups and clusters is high: ~1-3 kpc and \~10^9 M_sun, respectively. We use the high-resolution Adaptive Refinement Tree (ART) N-body code to follow the evolution of 256^3 dark matter particles with dynamic range in spatial resolution of 32,000 in a box of 60/h Mpc. We show that the correlation function of these halos is anti-biased with respect to the dark matter correlation function and is high and steeper than the correlation function of the isolated virialized objects. The correlation function evolves only mildly between z=3 and z=1. The mass evolution of isolated virialized objects determined from the simulation is in good agreement with prediction of semi-analytical models. The differences exist, however, if we include satellite halos in the halo catalogs., Contributions to the Proceedings of the MPA/ESO Conference "Evolution of LSS: from Recombination to Garching", Garching, August 2-7, 1998. LaTeX (uses tp.sty, psfig.sty), 5 pages, 5 figs

We study the evolution of the halo-halo correlation function and bias in a LCDM model using very high-resolution N-body simulations. The high force and mass resolution allows dark matter (DM) halos to survive in the tidal fields in high-density regions and thus prevents the ambiguities related with the ``overmerging problem.'' Numbers of galaxy-size halos in cluster-like objects in our simulation are similar to the numbers of galaxies observed in real clusters. This allows us to estimate for the first time the evolution of the correlation function and bias at small (down to ~100/h kpc) scales. We compare particle distribution, dark matter correlation function, density profiles, and halo mass function produced with our N-body code and corresponding results of the AP3M simulations. We find that at all epochs the 2-point correlation function of galaxy-size halos xihh is well approximated by a power-law with slope ~1.6-1.8. The difference between the shape of xihh and the shape of the correlation function of matter results in the scale-dependent bias at scales, Contributions to the Proceedings of the MPA/ESO Conference "Evolution of LSS: from Recombination to Garching", Garching, August 2-7, 1998. LaTeX (uses tp.sty, psfig.sty), 5 pages, 4 figs. Full-resolution version of Fig.2 is available at ftp://charon.nmsu.edu/pub/kravtsov/PAPERS/Garching98/

Formation of luminous matter in the Universe is a complicated process, which includes many processes and components. It is the vastly different scales involved in the process (from star formation on few parsec scales to galaxy clusters and superclusters on megaparsecs scales) and numerous ill-understood processes, which make the whole field a maze of unsolved, but exciting problems. We present new approximations for numerical treatment of multiphase ISM forming stars. The approximations were tested and calibrated using N-body+fluid numerical simulations. We specifically target issues related with effects of unresolved lumpinesses of the gas. We show that the degree of freedom is much smaller than naively expected because of self-regulating nature of the process of global star formation. One of the problems of numerical simulations is related with the uncertainties of approximation of the supernovae (SN) feedback. It is often assumed that the feedback is mainly due to momentum transferred by SN in to the ISM. We argue that this may not be true. We present a realistic example of gas actively forming stars with short cooling time for which the SF feedback is important, but the kinetic energy of the gas motion due to SN is negligibly small as compared with the thermal energy of gas., Comment: 5 pages, 2 PS figures, uses epsf.sty. Talk presented at the 12-th Postdam Cosmology Workshop, ``Large Scale Structure: Tracks and Traces'', World Scientific 1998. Available at http://astro.nmsu.edu/~aklypin/index.html

We present two matched sets of five simulations each, covering five presently favored simple modifications to the standard cold dark matter (CDM) scenario. One simulation suite, with a linear box size of 75 Mpc/h, is designed for high resolution and good statistics on the group/poor cluster scale, and the other, with a box size of 300 Mpc/h, is designed for good rich cluster statistics. All runs had 57 million cold particles, and models with massive neutrinos had an additional 113 million hot particles. We consider separately models with massive neutrinos, tilt, curvature, and a nonzero cosmological constant in addition to the standard CDM model. We find that our tilted Omega+Omega_Lambda=1 (TLCDM) model produces too much small-scale power by a factor of ~3, and our open Lambda=0 (OCDM) model also exceeds observed small-scale power by a factor of 2. In addition, we take advantage of the large dynamic range in detectable halo masses our simulations allow to check the shape of the Press-Schechter approximation. We find good fits at cluster masses for delta_c=1.27--1.35 for a Gaussian filter and delta_c=1.57--1.73 for a tophat filter. However, Press-Schechter overpredicts the number density of halos compared to the simulations in the high resolution suite by a weakly cosmology-dependent factor of 1.5--2 at galaxy and group masses, which cannot be fixed by adjusting delta_c within reasonable bounds. An appendix generalizes the spherical collapse model to any isotropic cosmology., Comment: 18 pages Latex using Monthly Notices style, with 13 inlined EPS figures. This version matches the one accepted by MNRAS. The appendix has been removed and may now be found instead at http://fozzie.gsfc.nasa.gov/thesis/appendixC.ps.gz

Particle-Mesh (PM) codes are still very useful tools for testing predictions of cosmological models in cases when extra high resolution is not very important. We release for public use a cosmological PM N-body code. We provide a complete package of routines needed to set initial conditions, to run the code, and to analyze the results. The package allows you to simulate models with numerous combinations of parameters: open/flat/closed background, with or without the cosmological constant, different values of the Hubble constant, with or without hot neutrinos, tilted or non-tilted initial spectra, different amount of baryons. Routines are included to measure the power spectrum and the density distribution function in your simulations, and a bound-density-maxima code for halo finding. We also provide results of test runs. A simulation with 256^3 mesh and 128^3 particles can be done in a couple of days on a typical workstation (70Mb of RAM are needed). To run simulations with 800^3 mesh and 256^3 particles one needs a computer with 1Gb memory and 1Gb disk space. The code has been successfully tested on an HP workstation and on a Sun workstation running Solaris. Most of the files (not tests) can be obtained from ftp://astro.nmsu.edu/pub/aklypin/PMCODE The package can be downloaded from http://astro.nmsu.edu/~aklypin/PM/pmcode/index.html We provide this tool as a service to the astronomical community, but we cannot guarantee results or publications., Comment: 18 pages, 3 Postscript figures, uses epsf.sty

Here we discuss what are perhaps the two most popular variants of CDM that might agree with the data: \lcdm\ and CHDM. While the predictions of COBE-normalized \lcdm\ and CHDM both agree well with the available data on scales of $\sim 10$ to $100 \hMpc$, each has potential virtues and defects. \lcdm\ with $��_0 \sim 0.3$ has the possible virtue of allowing a higher expansion rate $H_0$ for a given cosmic age $t_0$, but the defect of predicting too much fluctuation power on small scales. CHDM has less power on small scales, so its predictions appear to be in good agreement with data on the galaxy distribution, but it remains to be seen whether it predicts early enough galaxy formation to be compatible with the latest high-redshift data. Also, two very recent observational results favor high cosmic density, and thus favor $��=1$ models such as CHDM over \lcdm\ --- (1) the positive deceleration parameter $q_0>0$ measured using high-redshift Type Ia supernovae, and (2) the low primordial deuterium/hydrogen ratio measured in two different quasar absorption spectra. We try to identify ``best'' variants of both \lcdm\ and CHDM, and discuss critical observational tests for both models., 9 pages, postscript only, includes 3 figures, to be published in Proceedings of the International Conference on Sources and Detection of Dark Matter in the Universe, UCLA, February 1996, D. Cline and D. Sanders, editors, Nucl. Phys. B Proc. Suppl., in press

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 ($\Omega_\nu = 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\sigma$ (where $\sigma$ 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\sigma$ 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., Comment: 17 pages, 12 figures, LaTex (uses mn.sty). Accepted by MNRAS

A new class of geometric statistics for analyzing galaxy catalogs is presented. {\it 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 ($\Omega_\nu = 0.3$) simulations of Klypin \etal (1995), the CfA1-like sky catalogs of Nolthenius, Klypin and Primack (1994, 1995), and the CfA1 catalog. For 100 Mpc periodic simulation boxes ($H_0 = 50$ km s$^{-1}$ Mpc$^{-1}$), we find robust discrimination of over 4$\sigma$ (where $\sigma$ represents resampling errors) between CHDM and CDM. The {\it reduced filament statistics} show that CfA1 data is intermediate between CHDM and CDM, but more consistent with the CHDM models. Filament statistics provide robust and discriminatory shape statistics with which to test cosmological simulations of various models against present and future redshift survey data., Comment: 29 pages total LaTex format. 6 compressed, tarred, uuencoded postscript figures. Or, PS figs available from Romeel Dav\'e at rad@lick.ucolick.org. Submitted to MNRAS

Cosmological models with a positive cosmological constant and $\Omega_0, Comment: 12 pages + 2 tables + 4 figures are available from http://charon.nmsu.edu/ftp/aklypin/LCDM

Cold + Hot Dark Matter (CHDM) $\Omega=1$ cosmological models require a total neutrino mass $\sim 5$ eV. Because recent data support the $\nu_\mu \to \nu_\tau$ oscillation explanation of the cosmic ray $\nu_\mu$ deficit, which requires that $m(\nu_\mu) \approx m(\nu_\tau)$, this suggests that $m(\nu_\mu) \approx m(\nu_\tau) \approx 2.4$ eV. The linear calculations and N-body simulation reported here indicate that an $\Omega=1$ CHDM model with two 2.4 eV neutrinos (designated C$\nu^2$DM) agrees remarkably well with all available observations, but only if the Hubble parameter $h \approx 0.5$. We also show that even one 2.4 eV neutrino raises serious difficulties for low-$\Omega$ flat CDM models.

Cold + Hot Dark Matter (CHDM) $��=1$ cosmological models require a total neutrino mass $\sim 5$ eV. Because recent data support the $��_��\to ��_��$ oscillation explanation of the cosmic ray $��_��$ deficit, which requires that $m(��_��) \approx m(��_��)$, this suggests that $m(��_��) \approx m(��_��) \approx 2.4$ eV. The linear calculations and N-body simulation reported here indicate that an $��=1$ CHDM model with two 2.4 eV neutrinos (designated C$��^2$DM) agrees remarkably well with all available observations, but only if the Hubble parameter $h \approx 0.5$. We also show that even one 2.4 eV neutrino raises serious difficulties for low-$��$ flat CDM models.

This letter presents results of new high resolution $\Omega=1$ Cold + Hot Dark Matter (CHDM) and Cold Dark Matter (CDM) simulations. Properties of groups in these simulations reflect the lower small-scale velocities and the greater tendency to form distinct filaments on both small and large scales in CHDM as compared to CDM. The fraction of galaxies in groups and the median group rms velocity are found to be powerful discriminators between models. We combine these two features into a very robust statistic, median group rms velocity $v_{\rm gr}(f_{\rm gr})$ as a function of the fraction $f_{\rm gr}$ of galaxies in groups. Using this statistic, we compare ``observed'' simulations to CfA data in redshift space in a careful and consistent way. We find that CHDM remains a promising model, with for example $v_{\rmgr}(0.45) \approx 125 \pm 25 \kms$ in agreement with the CfA data, while CDM with bias b=1.0 (COBE-compatible) or b=1.5, both giving $v_{\rm gr}(0.45) \approx 400 \pm 25 \kms$, can be virtually ruled out. Using median $M/L$, the observed value of $\Omega$ is $0.10$ (CHDM) to $0.38$ (CDM)., Comment: to appear in Ap J Letters, 12 pages including 3 figures, uuencoded compressed postscript, preprint SCIPP 93/46

Using an isolated Milky Way-mass galaxy simulation, we compare results from nine state-of-the-art gravitohydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e. g., radiative cooling and UV background by the standardized package GRACKLE) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly formed stellar clump mass functions show more significant variation (difference by up to a factor of similar to 3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low-density region, and between more diffusive and less diffusive schemes in the high-density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.