Search

We identify voids as maximal non-overlapping spheres within the haloes of the Uchuu simulation and three smaller halo simulation boxes with smaller volume and different $\sigma_{8}$ values, and galaxies with redshift in the range $0.02

We measure the two point correlation function (CF) of 1357 galaxy clusters with a mass of $\log_{10}{M_{200}}\geq 13.6$~\hm~and at a redshift of $z \leq 0.125$. This work differs from previous analyses in that it utilizes a spectroscopic cluster catalogue, $\mathtt{SDSS-GalWCat}$, to measure the CF and detect the baryon acoustic oscillation (BAO) signal. Unlike previous studies which use statistical techniques, we compute covariance errors directly by generating a set of 1086 galaxy cluster lightcones from the GLAM $N$-body simulation. Fitting the CF with a power-law model of the form $\xi(s) = (s/s_0)^{-\gamma}$, we determine the best-fit correlation length and power-law index at three mass thresholds. We find that the correlation length increases with increasing the mass threshold while the power-law index is almost constant. For $\log_{10}{M_{200}}\geq 13.6$~\hm, we find $s_0 = 14.54\pm0.87$~\h~and $\gamma=1.97\pm0.11$. We detect the BAO signal at $s = 100$~\h~with a significance of $1.60 \sigma$. Fitting the CF with a $\Lambda$CDM model, we find $D_\mathrm{V}(z = 0.089)\mathrm{r}^{fid}_d/\mathrm{r}_d = 267.62 \pm 26$ \h, consistent with Planck 2015 cosmology. We present a set of 108 high-fidelity simulated galaxy cluster lightcones from the high-resolution \U~N-body simulation, employed for methodological validation. We find $D_\mathrm{V}(z = 0.089)/r_d = 2.666 \pm 0.129$, indicating that our method does not introduce any bias in the parameter estimation for this small sample of galaxy clusters., Comment: 7 pages, 5 figures

This study investigates the clustering and bias of Luminous Red Galaxies (LRG) in the BOSS-LOWZ, -CMASS, -COMB, and eBOSS samples, using two types of simulated lightcones: (i) high-fidelity lightcones from Uchuu $N$-body simulation, employing SHAM technique to assign LRG to (sub)halos, and (ii) 16000 covariance lightcones from GLAM-Uchuu $N$-body simulations, including LRG using HOD data from Uchuu. Our results indicate that Uchuu and GLAM lightcones closely replicate BOSS/eBOSS data, reproducing correlation function and power spectrum across scales from redshifts $0.2$ to $1.0$, from $2$ to $150\,h^{-1}\mathrm{Mpc}$ in configuration space, from $0.005$ to $0.7\,h\mathrm{Mpc}^{-1}$ in Fourier space, and across different LRG stellar masses. Furthermore, comparing with existing MD-Patchy and EZmock BOSS/eBOSS lightcones based on approximate methods, our GLAM-Uchuu lightcones provide more precise clustering estimates. We identify significant deviations from observations within $20\,h^{-1}\mathrm{Mpc}$ scales in MD-Patchy and EZmock, with our covariance matrices indicating that these methods underestimate errors by between $10\%$ and $60\%$. Lastly, we explore the impact of cosmology on galaxy clustering. Our findings suggest that, given the current level of uncertainties in BOSS/eBOSS data, distinguishing models with and without massive neutrino effects on LSS is challenging. This paper highlights the Uchuu and GLAM-Uchuu simulations' robustness in verifying the accuracy of Planck cosmological parameters, providing a strong foundation for enhancing lightcone construction in future LSS surveys. We also demonstrate that generating thousands of galaxy lightcones is feasible using $N$-body simulations with adequate mass and force resolution.

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

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 present the public data release of the Uchuu-UM galaxy catalogues by applying the UniverseMachine algorithm to assign galaxies to the dark matter halos in the Uchuu $N$-body cosmological simulation. It includes a variety of baryonic properties for all galaxies down to $\sim 5\times10^8 M_{\odot}$ with halos in a mass range of $10^{10}

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$, 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. 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 also examine the effect of parameters related to galaxy formation physics on the field variance. 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)., Comment: 18 pages, 16 figures, accepted for publication in MNRAS, Figs. 3, 4, 11, 13, and 14 revised, Figs. 14 and 16 in the original manuscript moved to the Appendix, All Uchuu-$\nu^2$GC data used in this paper are publicly available on our Skies & Universes website https://www.skiesanduniverses.org/Simulations/Uchuu/GalaxyCatalogues/

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., Comment: 54 pages, 14 figures; this paper is the twin paper of arXiv:2110.00566 submitted to arXiv on the same day, and the two papers have identical texts in part

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., Comment: 48 pages, 12 figures. This paper is the twin paper of arxiv:2110.00328 submitted to arXiv on the same day, and the two papers have identical texts in part

Methods. We perform numerical simulations of the evolution of the cosmic web for the conventional LCDM model. The simulations cover a wide range of box sizes L = 256 - 4000 Mpc/h, mass and force resolutions and epochs from very early moments z = 30 to the present moment z = 0. We calculate density fields with various smoothing lengths to find the dependence of the density field on smoothing scale. We calculate PDF and its moments - variance, skewness and kurtosis. Results. We focus on the third (skewness S) and fourth (kurtosis K) moments of the distribution functions: their dependence on the smoothing scale, the amplitude of fluctuations and the redshift. During the evolution the reduced skewness $S_3= S/\sigma$ and reduced kurtosis $S_4=K/\sigma^2$ present a complex behaviour: at a fixed redshift curves of $S_3(\sigma)$ and $S_4(\sigma)$ steeply increase with $\sigma$ at $\sigma\le 1$ and then flatten out and become constant at $\sigma\ge2$. If we fix the smoothing scale $R_t$, then after reaching the maximum at $\sigma\approx 2$, the curves at large $\sigma$ start to gradually decline. We provide accurate fits for the evolution of $S_{3,4}(\sigma,z)$. Skewness and kurtosis approach at early epochs constant levels, depending on smoothing length: $S_3(\sigma) \approx 3$ and $S_4(\sigma) \approx 15$. Conclusions. Most of statistics of dark matter clustering (e.g., halo mass function or concentration-mass relation) are nearly universal: they mostly depend on the $\sigma$ with the relatively modest correction to explicit dependence on the redshift. We find just the opposite for skewness and kurtosis: the dependence of moments on evolutionary epoch $z$ and smoothing length $R_t$ is very different, together they determine the evolution of $S_{3,4}(\sigma)$ uniquely. The evolution of $S_3$ and $S_4$ cannot be described by current theoretical approximations., Comment: 17 pages, 14 figures, revised version accepted by Astronomy and Astrophysics

The eROSITA X-ray telescope on board the Spectrum-Roentgen-Gamma (SRG) mission will measure the position and properties of about 100,000 clusters of galaxies and 3 million active galactic nuclei over the full sky. To study the statistical properties of this ongoing survey, it is key to estimate the selection function accurately. We create a set of full sky light-cones using the MultiDark and UNIT dark matter only N-body simulations. We present a novel method to predict the X-ray emission of galaxy clusters. Given a set of dark matter halo properties (mass, redshift, ellipticity, offset parameter), we construct an X-ray emissivity profile and image for each halo in the light-cone. We follow the eROSITA scanning strategy to produce a list of X-ray photons on the full sky. We predict scaling relations for the model clusters, which are in good agreement with the literature. The predicted number density of clusters as a function of flux also agrees with previous measurements. Finally, we obtain a scatter of 0.21 (0.07, 0.25) for the X-ray luminosity -- mass (temperature -- mass, luminosity -- temperature) model scaling relations. We provide catalogues with the model photons emitted by clusters and active galactic nuclei. These catalogues will aid the eROSITA end to end simulation flow analysis and in particular the source detection process and cataloguing methods., Comment: 15 pages, 11 figures, accepted in the open journal of astrophysics

We introduce the Uchuu suite of large high-resolution cosmological $N$-body simulations. The largest simulation, named Uchuu, consists of 2.1 trillion ($12800^3$) dark matter particles in a box of side-length 2.0 Gpc/h, with particle mass $3.27 \times 10^{8}$ Msun/h. The highest resolution simulation, Shin-Uchuu, consists of 262 billion ($6400^3$) particles in a box of side-length 140 Mpc/h, with particle mass $8.97 \times 10^{5}$ Msun/h. Combining these simulations we can follow the evolution of dark matter halos and subhalos spanning those hosting dwarf galaxies to massive galaxy clusters across an unprecedented volume. In this first paper, we present basic statistics, dark matter power spectra, and the halo and subhalo mass functions, which demonstrate the wide dynamic range and superb statistics of the Uchuu suite. From an analysis of the evolution of the power spectra we conclude that our simulations remain accurate from the Baryon Acoustic Oscillation scale down to the very small. We also provide parameters of a mass-concentration model, which describes the evolution of halo concentration and reproduces our simulation data to within 5 per cent for halos with masses spanning nearly eight orders of magnitude at redshift 00.5, whereas no upturn is detected at z<0.5. We make publicly available various $N$-body products as part of Uchuu Data Release 1 on the Skies & Universes site. Future releases will include gravitational lensing maps and mock galaxy, X-ray cluster, and active galactic nuclei catalogues., Comment: 22 pages, 14 figures, accepted by MNRAS. We release various $N$-body products as data release 1 on http://skiesanduniverses.org/ such as subsets of simulation particles, matter power spectra, halo/subhalo catalogues, and their merger trees

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., Comment: 12 pages, 11 figures, submitted to MNRAS

LCDM cosmological models with Early Dark Energy (EDE) have been proposed to resolve tensions between the Hubble constant H0 = 100h km/s/Mpc measured locally, giving h ~ 0.73, and H0 deduced from Planck cosmic microwave background (CMB) and other early universe measurements plus LCDM, giving h ~ 0.67. EDE models do this by adding a scalar field that temporarily adds dark energy equal to about 10% of the cosmological energy density at the end of the radiation-dominated era at redshift z ~ 3500. Here we compare linear and nonlinear predictions of a Planck-normalized LCDM model including EDE giving h = 0.728 with those of standard Planck-normalized LCDM with h = 0.678. We find that nonlinear evolution reduces the differences between power spectra of fluctuations at low redshifts. As a result, at z = 0 the halo mass functions on galactic scales are nearly the same, with differences only 1-2%. However, the differences dramatically increase at high redshifts. The EDE model predicts 50% more massive clusters at z = 1 and twice more galaxy-mass halos at z = 4. Even greater increases in abundances of galaxy-mass halos at higher redshifts may make it easier to reionize the universe with EDE. Predicted galaxy abundances and clustering will soon be tested by JWST observations. Positions of baryonic acoustic oscillations (BAOs) and correlation functions differ by about 2% between the models -- an effect that is not washed out by nonlinearities. Both standard LCDM and the EDE model studied here agree well with presently available acoustic-scale observations, but DESI and Euclid measurements will provide stringent new tests., Comment: 13 pages, 12 figures, submitted to MNRAS

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., Comment: 22 pages, 18 figure, submitted to MNRAS, the galaxy catalogues are available at http://www.skiesanduniverses.org/

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

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 $\Delta = 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(\sigma_{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 $\delta_{\log\sigma} = 0.06$. The catalogs are publicly available at the following website\footnote{\url{https://mohamed-elhashash-94.webself.net/galwcat/}}., Comment: 15 pages, 6 figures

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., Comment: Submitted to MNRAS

We present the first galaxy-OVI absorption kinematic study for 20 absorption systems (EW>0.1~{\AA}) 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., Comment: 24 pages, 21 figures, 4 tables. Accepted to ApJ on November 14, 2018

Halo bias is the main link between the matter distribution and dark matter halos. In its simplest form, halo bias is determined by halo mass, but there are known additional dependencies on other halo properties which are of consequence for accurate modeling of galaxy clustering. Here we present the most precise measurement of these secondary-bias dependencies on halo age, concentration, and spin, for a wide range of halo masses spanning from 10$^{10.7}$ to 10$^{14.7}$ $h^{-1}$ M$_{\odot}$. At the high-mass end, we find no strong evidence of assembly bias for masses above M$_{vir}$ $\sim10^{14}$ $h^{-1}$ M$_{\odot}$. Secondary bias exists, however, for halo concentration and spin, up to cluster-size halos, in agreement with previous findings. For halo spin, we report, for the first time, two different regimes: above M$_{vir}\sim$10$^{11.5}$ $h^{-1}$ M$_{\odot}$, halos with larger values of spin have larger bias, at fixed mass, with the effect reaching almost a factor 2. This trend reverses below this characteristic mass. In addition to these results, we test, for the first time, the performance of a multi-tracer method for the determination of the relative bias between different subsets of halos. We show that this method increases significantly the signal-to-noise of the secondary-bias measurement as compared to a traditional approach. This analysis serves as the basis for follow-up applications of our multi-tracer method to real data., Comment: 11 pages, 6 figures, submitted to MNRAS

In order to capture as much information as possible large galaxy surveys have been increasing their volume and redshift depth. To face this challenge theory has responded by making cosmological simulations of huge computational volumes with equally increasing the number of dark matter particles and supercomputing resources. Thus, it is taken for granted that the ideal situation is when a single computational box encompasses the whole effective volume of the observational survey, e.g., ~50 Gpch^3 for the DESI and Euclid surveys. Here we study the effects of missing long-waves in a finite volume using several relevant statistics: the abundance of dark matter halos, the PDF, the correlation function and power spectrum, and covariance matrices. Finite volume effects can substantially modify the results if the computational volumes are less than ~(500Mpch)^3. However, the effects become extremely small and practically can be ignored when the box-size exceeds ~1Gpch^3. We find that the average power spectra of dark matter fluctuations show remarkable lack of dependence on the computational box-size with less than 0.1% differences between 1Gpch and 4Gpch boxes. No measurable differences are expected for the halo mass functions for these volumes. The covariance matrices are scaled trivially with volume, and small corrections due to super-sample modes can be added. We conclude that there is no need to make those extremely large simulations when a box-size of 1-1.5Gpch is sufficient to fulfil most of the survey science requirements., Comment: 15 pages, 14 figures, accepted to MNRAS

We introduce GalWeight, a new technique for assigning galaxy cluster membership. This technique is specifically designed to simultaneously maximize the number of bona fide cluster members while minimizing the number of contaminating interlopers. The GalWeight technique can be applied to both massive galaxy clusters and poor galaxy groups. Moreover, it is effective in identifying members in both the virial and infall regions with high efficiency. We apply the GalWeight technique to MDPL2 \& Bolshoi N-body simulations, and find that it is $> 98\%$ accurate in correctly assigning cluster membership. We show that GalWeight compares very favorably against four well-known existing cluster membership techniques (shifting gapper, den Hartog, caustic, SIM). We also apply the GalWeight technique to a sample of twelve Abell clusters (including the Coma cluster) using observations from the Sloan Digital Sky Survey. We end by discussing GalWeight's potential for other astrophysical applications., Comment: Accepted for publication in APJ

Cosmological simulations are fundamental tools to study structure formation and the astrophysics of evolving structures, in particular clusters of galaxies. While hydrodynamical simulations cannot sample efficiently large volumes and explore different cosmologies at the same time, N-body simulations lack the baryonic physics that is crucial to determine the observed properties of clusters. One solution is to use (semi-)analytical models to implement the needed baryonic physics. In this way, we can generate the many mock universes that will be required to fully exploit future large sky surveys, such as that from the upcoming eROSITA X-ray telescope. We developed a phenomenological model based on observations of clusters to implement gas density and temperature information on the dark-matter-only halos of the MultiDark simulations. We generate several full-sky mock light-cones of clusters for the WMAP and Planck cosmologies, adopting different parameters in our phenomenological model of the intra-cluster medium. For one of these simulations and models, we also generate 100 light-cones corresponding to 100 random observers and explore the variance among them in several quantities. In this first paper on MultiDark mock galaxy cluster light-cones, we focus on presenting our methodology and discuss predictions for eROSITA, in particular, exploring the potential of angular power spectrum analyses of its detected (and undetected) cluster population to study X-ray scaling relations, the intra-cluster medium, and the composition of the cosmic X-ray background. We make publicly available on-line more than 400 GB of light-cones, which include the expected eROSITA count rate, on Skies & Universes (http://www.skiesanduniverses.org)., Comment: 21 pages, 13 figures, 2 tables. Accepted for publication in MNRAS

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

We present the public release of the MultiDark-Galaxies: three distinct galaxy catalogues derived from one of the Planck cosmology MultiDark simulations (i.e. MDPL2, with a volume of (1 Gpc/$h$)$^{3}$ and mass resolution of $1.5 \times 10^{9} M_{\odot}/h$) by applying the semi-analytic models GALACTICUS, SAG, and SAGE to it. We compare the three models and their conformity with observational data for a selection of fundamental properties of galaxies like stellar mass function, star formation rate, cold gas fractions, and metallicities - noting that they sometimes perform differently reflecting model designs and calibrations. We have further selected galaxy subsamples of the catalogues by number densities in stellar mass, cold gas mass, and star formation rate in order to study the clustering statistics of galaxies. We show that despite different treatment of orphan galaxies, i.e. galaxies that lost their dark-matter host halo due to the finite mass resolution of the N-body simulation or tidal stripping, the clustering signal is comparable, and reproduces the observations in all three models - in particular when selecting samples based upon stellar mass. Our catalogues provide a powerful tool to study galaxy formation within a volume comparable to those probed by on-going and future photometric and redshift surveys. All model data consisting of a range of galaxy properties - including broad-band SDSS magnitudes - are publicly available., Comment: 29 pages, 16 figures, 8 tables, accepted for publication in MNRAS. All data incl. the complete galaxy catalogues for all models are publicly available from the CosmoSim database (http://www.cosmosim.org); a selected set of galaxy properties is available via the Skies & Universes website (http://www.skiesanduniverses.org)

The one-point probability distribution function (PDF) of the matter density field in the universe is a fundamental property that plays an essential role in cosmology for estimates such as gravitational weak lensing, non-linear clustering, massive production of mock galaxy catalogs, and testing predictions of cosmological models. Here we make a comprehensive analysis of the dark matter PDF using a suite of 7000 N-body simulations that covers a wide range of numerical and cosmological parameters. We find that the PDF has a simple shape: it declines with density as a power-law P~rho**(-2), which is exponentially suppressed on both small and large densities. The proposed double-exponential approximation provides an accurate fit to all our N-body results for small filtering scales R< 5Mpc/h with rms density fluctuations sigma>1. In combination with the spherical infall model that works well for small fluctuations sigma<1, the PDF is now approximated with just few percent errors over the range of twelve orders of magnitude -- a remarkable example of precision cosmology. We find that at 5-10% level the PDF explicitly depends on redshift (at fixed sigma) and on cosmological density parameter Omega_m. We test different existing analytical approximations and find that the often used log-normal approximation is always 3-5 times less accurate than either the double-exponential approximation or the spherical infall model.

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

$N$-body cosmological simulations are an essential tool to understand the observed distribution of galaxies. We use the MultiDark simulation suite, run with the Planck cosmological parameters, to revisit the mass and velocity functions. At redshift $z=0$, the simulations cover four orders of magnitude in halo mass from $\sim10^{11}M_\odot$ with 8,783,874 distinct halos and 532,533 subhalos. The total volume used is $\sim$515 Gpc$^3$, more than 8 times larger than in previous studies. We measure and model the halo mass function, its covariance matrix w.r.t halo mass and the large scale halo bias. With the formalism of the excursion-set mass function, we explicit the tight interconnection between the covariance matrix, bias and halo mass function. We obtain a very accurate ($<2\%$ level) model of the distinct halo mass function. We also model the subhalo mass function and its relation to the distinct halo mass function. The set of models obtained provides a complete and precise framework for the description of halos in the concordance Planck cosmology. Finally, we provide precise analytical fits of the $V_{max}$ maximum velocity function up to redshift $z<2.3$ to push for the development of halo occupation distribution using $V_{max}$. The data and the analysis code are made publicly available in the \textit{Skies and Universes} database., Comment: Corresponding data is available at the Skies and Universes data base: http://projects.ift.uam-csic.es/skies-universes

Upcoming and existing large-scale surveys of galaxies require accurate theoretical predictions of the dark matter clustering statistics for thousands of mock galaxy catalogs. We demonstrate that this goal can be achieve with our new Parallel Particle-Mesh (PM) Nbody code (PPM-GLAM) at a very low computational cost. We run about 15,000 simulations with ~2 billion particles that provide ~1% accuracy of the dark matter power spectra P(k) for wave-numbers up to k~ 1h/Mpc. Using this large data-set we study the power spectrum covariance matrix, the stepping stone for producing mock catalogs. In contrast to many previous analytical and numerical results, we find that the covariance matrix normalised to the power spectrum C(k,k')/P(k)P(k') has a complex structure of non-diagonal components. It has an upturn at small k, followed by a minimum at k=0.1-0.2h/Mpc. It also has a maximum at k=0.5-0.6h/Mpc. The normalised covariance matrix strongly evolves with redshift: C(k,k')~delta(t)^alpha P(k)P(k'), where delta is the linear growth factor and alpha ~ 1-1.25, which indicates that the covariance matrix depends on cosmological parameters. We also show that waves longer than 1Gpc have very little impact on the power spectrum and covariance matrix. This significantly reduces the computational costs and complexity of theoretical predictions: relatively small volume ~ (1Gpc)^3 simulations capture the necessary properties of dark matter clustering statistics. All the power spectra obtained from many thousands of our simulations are publicly available., Comment: 20 pages, 13 figures

In current and future surveys, quasars play a key role. The new data will extend our knowledge of the Universe as it will be used to better constrain the cosmological model at redshift $z>1$ via baryon acoustic oscillation and redshift space distortion measurements. Here, we present the first clustering study of quasars observed by the extended Baryon Oscillation Spectroscopic Survey. We measure the clustering of $\sim 70,000$ quasars located in the redshift range $0.9

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

Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics 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 ~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., Comment: 28 pages, 35 figures, Accepted for publication in the Astrophysical Journal, Image resolution greatly reduced, High-resolution version of this article is available at http://www.jihoonkim.org/agora/AGORA_Paper4_draft.pdf, The first paper of the AGORA Initiative is at http://adsabs.harvard.edu/abs/2014ApJS..210...14K, More information on AGORA is at http://www.AGORAsimulations.org/

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., Comment: 8 pages, 5 figures. Accepted for publication in ApJ

We report and provide fitting functions for the abundance of dark matter halos and subhalos as a function of mass, circular velocity, and redshift from the new Bolshoi-Planck and MultiDark-Planck $\Lambda$CDM cosmological simulations, based on the Planck cosmological parameters. We also report the halo mass accretion rates, which may be connected with galaxy star formation rates. We show that the higher cosmological matter density of the Planck parameters compared with the WMAP parameters leads to higher abundance of massive halos at high redshifts. We find that the median halo spin parameter $\lambda_{\rm B} = J(2M_{\rm vir}R_{\rm vir}V_{\rm vir})^{-1}$ is nearly independent of redshift, leading to predicted evolution of galaxy sizes that is consistent with observations, while the significant decrease with redshift in median $\lambda_{\rm P} = J|E|^{-1/2}G^{-1}M^{-5/2}$ predicts more decrease in galaxy sizes than is observed. Using the Tully-Fisher and Faber-Jackson relations between galaxy velocity and mass, we show that a simple model of how galaxy velocity is related to halo maximum circular velocity leads to increasing overprediction of cosmic stellar mass density as redshift increases beyond redshifts $z\sim1$, implying that such velocity-mass relations must change at redshifts $z>1$. By making a realistic model of how observed galaxy velocities are related to halo circular velocity, we show that recent optical and radio observations of the abundance of galaxies are in good agreement with our $\Lambda$CDM simulations. Our halo demographics are based on updated versions of the \rockstar\ and \ctrees\ codes, and this paper includes appendices explaining all of their outputs. This paper is an introduction to a series of related papers presenting other analyses of the Bolshoi-Planck and MultiDark-Planck simulations., Comment: 27 pages, 26 figures. To match the published version at MNRAS. Rockstar and Consistent Trees outputs are available at http://hipacc.ucsc.edu/Bolshoi/MergerTrees.html

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 $\Lambda$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-$\sigma$ 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-$\sigma$. The BigMD-BOSS light-cone including properties of BOSS galaxies and halo properties is made publicly available., Comment: 16 pages, 16 figures. See companion papers that share the "The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey:" title

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 $\sigma$, 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., Comment: 20 pages, 15 figures; accepted for publication at MNRAS; adjusted to match the published version

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 $\chi^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., Comment: Accepted for publication in MNRAS

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

We present the measurements and modelling of the small-to-intermediate scale (0.1--25 Mpc/h) projected and three-dimensional (3D) redshift-space two-point correlation functions (2PCFs) of local galaxies in the Sloan Digital Sky Survey (SDSS) Data Release 7. We find a clear dependence of galaxy clustering on luminosity in both projected and redshift spaces, generally being stronger for more luminous samples. The measurements are successfully interpreted within the halo occupation distribution (HOD) framework with central and satellite velocity bias parameters to describe galaxy kinematics inside haloes and to model redshift-space distortion (RSD) effects. In agreement with previous studies, we find that more luminous galaxies reside in more massive haloes. Including the redshift-space 2PCFs helps tighten the HOD constraints. Moreover, we find that luminous central galaxies are not at rest at the halo centres, with the velocity dispersion about 30% that of the dark matter. Such a relative motion may reflect the consequence of galaxy and halo mergers, and we find that central galaxies in lower mass haloes tend to be more relaxed with respect to their host haloes. The motion of satellite galaxies in luminous samples is consistent with their following that of the dark matter. For faint samples, satellites tends to have slower motion, with velocity dispersion inside haloes about 85% that of the dark matter. We discuss possible applications of the velocity bias constraints on studying galaxy evolution and cosmology. In the appendix, we characterize the distribution of galaxy redshift measurement errors, which is well described by a Gaussian-convolved double exponential distribution., Comment: Fixed the Arxiv pdf generator bug for Figure 4. Figure 4 highlights the agreement between our models and the data

We present a comparison of major methodologies of fast generating mock halo or galaxy catalogues. The comparison is done for two-point and the three-point clustering statistics. The reference catalogues are drawn from the BigMultiDark N-body simulation. Both friend-of-friends (including distinct halos only) and spherical overdensity (including distinct halos and subhalos) catalogs have been used with the typical number density of a large-volume galaxy surveys. We demonstrate that a proper biasing model is essential for reproducing the power spectrum at quasilinear and even smaller scales. With respect to various clustering statistics a methodology based on perturbation theory and a realistic biasing model leads to very good agreement with N-body simulations. However, for the quadrupole of the correlation function or the power spectrum, only the method based on semi-N-body simulation could reach high accuracy (1% level) at small scales, i.e., r<25 Mpc/h or k>0.15 h/Mpc. Full N-body solutions will remain indispensable to produce reference catalogues. Nevertheless, we have demonstrated that the far more efficient approximate solvers can reach a few percent accuracy in terms of clustering statistics at the scales interesting for the large-scale structure analysis after calibration with a few reference N-body calculations. This makes them useful for massive production aimed at covariance studies, to scan large parameter spaces, and to estimate uncertainties in data analysis techniques, such as baryon acoustic oscillation reconstruction, redshift distortion measurements, etc., Comment: 16 pages, 15 figures; matches the version accepted by MNRAS; a bug in PINOCCHIO code has been fixed; no major modification from previous version

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., Comment: 22 pages, 23 figures, accepted by MNRAS