Source author record

Y. P. Jing

Y. P. Jing appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

43works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

43 published item(s)

preprint2022arXiv

Construct the emission line galaxy-host halo connection through auto and cross correlations

We investigate the [O\,II] emission line galaxy (ELG)-host halo connection via auto and cross correlations, and propose a concise and effective method to populate ELGs in dark matter halos without assuming a parameterized halo occupation distribution (HOD) model. Using the observational data from VIMOS Public Extragalactic Redshift Survey (VIPERS), we measure the auto and cross correlation functions between ELGs selected by [O\,II] luminosity and normal galaxies selected by stellar mass. Combining the stellar-halo mass relation (SHMR) derived for the normal galaxies and the fraction of ELGs observed in the normal galaxy population, we demonstrate that we can establish an accurate ELG-halo connection. With the ELG-halo connection, we can accurately reproduce the auto and cross correlation functions of ELGs and normal galaxies both in real-space and in redshift-space, once the satellite fraction is properly reduced. Our method provides a novel strategy to generate ELG mock catalogs for ongoing and upcoming galaxy redshift surveys. We also provide a simple description for the HOD of ELGs.

preprint2020arXiv

Constraining the Milky Way Mass Profile with Phase-Space Distribution of Satellite Galaxies

We estimate the Milky Way (MW) halo properties using satellite kinematic data including the latest measurements from Gaia DR2. With a simulation-based 6D phase-space distribution function (DF) of satellite kinematics, we can infer halo properties efficiently and without bias, and handle the selection function and measurement errors rigorously in the Bayesian framework. Applying our DF from the EAGLE simulation to 28 satellites, we obtain an MW halo mass of $M=1.23_{-0.18}^{+0.21}\times 10^{12} M_\odot$ and a concentration of $c=9.4_{ -2.1}^{ +2.8}$ with the prior based on the $M$-$c$ relation. The inferred mass profile is consistent with previous measurements but with better precision and reliability due to the improved methodology and data. Potential improvement is illustrated by combining satellite data and stellar rotation curves. Using our EAGLE DF and best-fit MW potential, we provide much more precise estimates of kinematics for those satellites with uncertain measurements. Compared to the EAGLE DF, which matches the observed satellite kinematics very well, the DF from the semi-analytical model based on the dark-matter-only simulation Millennium II (SAM-MII) over-represents satellites with small radii and velocities. We attribute this difference to less disruption of satellites with small pericenter distances in the SAM-MII simulation. By varying the disruption rate of such satellites in this simulation, we estimate a $\sim 5\%$ scatter in the inferred MW halo mass among hydrodynamics-based simulations.

preprint2020arXiv

Orbital distribution of infalling satellite halos across cosmic time

The initial orbits of infalling subhalos largely determine the subsequent evolution of the subhalos and satellite galaxies therein and shed light on the assembly of their hosts. Using a large set of cosmological simulations of various resolutions, we quantify the orbital distribution of subhalos at infall time and its mass and redshift dependence in a large dynamic range. We further provide a unified and accurate model validated across cosmic time, which can serve as the initial condition for semi-analytic models. We find that the infall velocity $v$ follows a nearly universal distribution peaked near the host virial velocity $V_{\mathrm{h}}$ for any subhalo mass or redshift, while the infall orbit is most radially biased when $v\sim V_{\mathrm{h}}$. Moreover, subhalos that have a higher host mass or a higher sub-to-host ratio tend to move along a more radial direction with a relatively smaller angular momentum than their low host mass or low sub-to-host ratio counterparts, though they share the same normalized orbital energy. These relations are nearly independent of the redshift when using the density peak height as the proxy for host halo mass. The above trends are consistent with the scenario where the dynamical environment is relatively colder for more massive structures because their own gravity is more likely to dominate the local potentials. Based on this understanding, the more massive or isolated halos are expected to have higher velocity anisotropy.

preprint2020arXiv

Universal conditional distribution function of [OII] luminosity of galaxies, and prediction for the [OII] luminosity function at redshift $z<3$

The star-forming emission line galaxies (ELGs) with strong [OII] doublet are one of the main spectroscopic targets for the ongoing and upcoming fourth generation galaxy redshift surveys. In this work, we measure the [OII] luminosity $L_{\mathrm{[OII]}}$ and the near-ultraviolet band absolute magnitude $M_{\mathrm{NUV}}$ for a large sample of galaxies in the redshift range of $0.6\leq z <1.45$ from the Public Data Release 2 (PDR-2) of the VIMOS Public Extragalactic Redshift Survey (VIPERS). We aim to construct the intrinsic relationship between the $L_{\mathrm{[OII]}}$ and $M_{\mathrm{NUV}}$ through Bayesian analysis. In particular, we develop two different methods to properly correct for the incompleteness effect and observational errors in the [OII] emission line measurement. Our results indicate that the conditional distribution of $L_{\mathrm{[OII]}}$ at a given $M_{\mathrm{NUV}}$ can be well described by a universal probability distribution function (PDF), which is independent of $M_{\mathrm{NUV}}$ or redshift. Convolving the $L_{\mathrm{[OII]}}$ conditional PDF with the NUV Luminosity function (LF) available in the literature, we make a prediction for [OII] LFs at $z<3$. The predicted [OII] LFs are in good agreement with the observational results from the literature. Finally, we utilize the predicted [OII] LFs to estimate the number counts of [OII] emitters for the Subaru Prime Focus Spectrograph (PFS) survey. This universal conditional PDF of $L_{\mathrm{[OII]}}$ provides a novel way to optimize the source targeting strategy for [OII] emitters in future galaxy redshift surveys, and to model [OII] emitters in theories of galaxy formation.

preprint2016arXiv

Clustering properties and halo masses for central galaxies in the local Universe

We investigate the clustering and dark matter halo mass for a sample of $\sim$16,000 central galaxies selected from the SDSS/DR7 group catalog. We select subsamples of central galaxies on three two-dimensional planes, each formed by stellar mass ($M_{\star}$) and one of the other properties including optical color (g-r) , surface stellar mass density ($μ_{\star}$) and central stellar velocity dispersion ($σ_{\star}$). For each subsample we measure both the projected cross-correlation function ($w_{p}(r_{p})$) relative to a reference galaxy sample, and an average mass of the host dark matter halos ($M_{halo}$) . For comparison we have also estimated the $w_{p}(r_{p})$ for the full galaxy population and the subset of satellite galaxies. We find that, for central galaxies, both $w_{p}(r_{p})$ and $M_{halo}$ show strongest dependence on $M_{\star}$, and there is no clear dependence on other properties when stellar mass is fixed. This result provides strong support to the previously-adopted assumption that, for central galaxies, stellar mass is the galaxy property that is best indicative of the host dark halo mass. The full galaxy population and the subset of satellites show similar clustering properties in all cases. However, they are similar to the centrals only at high masses ($M_{\star} \geq 10^{11} M_{\odot}$). At low-mass ($M_{\star} \leq 10^{11} M_{\odot}$), the results are different: the $w_{p}(r_{p})$ increases with both $σ$ and g-r when $M_{\star}$ is fixed, and depends very weakly on $M_{\star}$ when $σ_{\star}$ or g-r is fixed. At fixed $M_{star}$, the $μ_{\star}$ shows weak correlations with clustering amplitude (and halo mass for central galaxies). Our results suggest that it is necessary to consider central and satellite galaxies separately when studying the link between galaxies and dark matter halos. (Abridged)

preprint2016arXiv

ELUCID - Exploring the Local Universe with reConstructed Initial Density field III: Constrained Simulation in the SDSS Volume

A method we developed recently for the reconstruction of the initial density field in the nearby Universe is applied to the Sloan Digital Sky Survey Data Release 7. A high-resolution N-body constrained simulation (CS) of the reconstructed initial condition, with $3072^3$ particles evolved in a 500 Mpc/h box, is carried out and analyzed in terms of the statistical properties of the final density field and its relation with the distribution of SDSS galaxies. We find that the statistical properties of the cosmic web and the halo populations are accurately reproduced in the CS. The galaxy density field is strongly correlated with the CS density field, with a bias that depend on both galaxy luminosity and color. Our further investigations show that the CS provides robust quantities describing the environments within which the observed galaxies and galaxy systems reside. Cosmic variance is greatly reduced in the CS so that the statistical uncertainties can be controlled effectively even for samples of small volumes.

preprint2015arXiv

A Chemical Evolution Model for the Fornax Dwarf Spheroidal Galaxy

Fornax is the brightest Milky Way (MW) dwarf spheroidal galaxy and its star formation history (SFH) has been derived from observations. We estimate the time evolution of its gas mass and net inflow and outflow rates from the SFH using a simple star formation law that relates the star formation rate to the gas mass. We present a chemical evolution model on a 2D mass grid with supernovae (SNe) as sources of metal enrichment. We find that a key parameter controlling the enrichment is the mass M_x of the gas to mix with the ejecta from each SN. The choice of M_x depends on the evolution of SN remnants and on the global gas dynamics. It differs between the two types of SNe involved and between the periods before and after Fornax became an MW satellite at time t = t_sat . Our results indicate that due to the global gas outflow at t > t_sat , part of the ejecta from each SN may directly escape from Fornax. Sample results from our model are presented and compared with data.

preprint2015arXiv

An extended view of the Pisces Overdensity from the SCUSS survey

SCUSS is a u-band photometric survey covering about 4000 square degree of the South Galactic Cap, reaching depths of up to 23 mag. By extending around 1.5 mag deeper than SDSS single-epoch u data, SCUSS is able to probe much a larger volume of the outer halo, i.e. with SCUSS data blue horizontal branch (BHB) stars can trace the outer halo of the Milky Way as far as 100-150 kpc. Utilizing this advantage we combine SCUSS u band with SDSS DR9 gri photometric bands to identify BHB stars and explore halo substructures. We confirm the existence of the Pisces overdensity, which is a structure in the outer halo (at around 80 kpc) that was discovered using RR Lyrae stars. For the first time we are able to determine its spatial extent, finding that it appears to be part of a stream with a clear distance gradient. The stream, which is ~5 degrees wide and stretches along ~25 degrees, consists of 20-30 BHBs with a total significance of around 6sigma over the background. Assuming we have detected the entire stream and that the progenitor has fully disrupted, then the number of BHBs suggests the original system was similar to smaller classical or a larger ultra-faint dwarf galaxy. On the other hand, if the progenitor still exists, it can be hunted for by reconstructing its orbit from the distance gradient of the stream. This new picture of the Pisces overdensity sheds new light on the origin of this intriguing system.

preprint2015arXiv

Estimating the evolution of gas in the Fornax dwarf spheroidal galaxy from its star formation history: an illustrative example

We propose that detailed data on the star formation history of a dwarf spheroidal galaxy (dSph) may be used to estimate the evolution of the total mass M_g(t) for cold gas in its star-forming disk. Using Fornax as an illustrative example, we estimate its M_g(t) and the corresponding net gas flow rate Delta F(t) assuming a global star formation rate psi(t) = lambda_*(t)[M_g(t)/M_sun]^alpha consistent with observations of nearby galaxies. We identify the onset of the transition in Delta F(t) from a net inflow to a net outflow as the time t_sat at which the Fornax halo became a Milky Way satellite and estimate the evolution of its total mass M_h(t) at t < t_sat using the median halo growth history in the LambdaCDM cosmology and its present mass within the half-light radius derived from observations. We examine three different cases of alpha = 1, 1.5, and 2, and justify the corresponding lambda_*(t) by comparing the gas mass fraction f_g(t) = M_g(t)/M_h(t) at t < t_sat with results from simulations of gas accretion by halos in a reionized universe. We find that the Fornax halo grew to M_h(t_sat) ~ 2x10^9 M_sun at t_sat ~ 5 or 8 Gyr, in broad agreement with previous studies using data on its stellar kinematics and its orbital motion. We describe qualitatively the evolution of Fornax as a satellite and discuss potential extension of our approach to other dSphs.

preprint2015arXiv

Modelling The Redshift-Space Three-Point Correlation Function in SDSS-III

We present the measurements of the redshift-space three-point correlation function (3PCF) for z~0.5 luminous red galaxies of the CMASS sample in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Data Release 11. The 3PCF measurements are interpreted within the halo occupation distribution framework using high-resolution N-body simulations, and the model successfully reproduces the 3PCF on scales larger than 1Mpc/h. As with the case for the redshift-space two-point correlation functions, we find that the redshift-space 3PCF measurements also favour the inclusion of galaxy velocity bias in the model. In particular, the central galaxy in a halo is on average in motion with respect to the core of the halo. We discuss the potential of the small-scale 3PCF to tighten the constraints on the relation between galaxies and dark matter haloes and on the phase-space distribution of galaxies.

preprint2015arXiv

The stellar-to-halo mass relations of local galaxies segregates by color

We derive the stellar-to-halo mass relations, SHMR, of local blue and red central galaxies separately, as well as the fraction of halos hosting blue/red central galaxies. We find that: 1) the SHMR of central galaxies is segregated by color, with blue centrals having a SHMR above the one of red centrals; at logMh~12, the Ms/Mh ratio of the blue centrals is ~0.05, which is ~1.7 times larger than the value of red centrals. 2) The intrinsic scatters of the SHMRs of red and blue centrals are ~0.14 and ~0.11dex, respectively. The intrinsic scatter of the average SHMR of all central galaxies changes from ~0.20dex to ~0.14dex in the 11.3<logMh<15 range. 3) The fraction of halos hosting blue centrals at Mh=1E11Msun is 87%, but at 2x1E12Msun decays to ~20%, approaching to a few per cents at higher masses. The characteristic mass at which this fraction is the same for blue and red galaxies is Mh~7x1E11Msun. Our results suggest that the SHMR of central galaxies at large masses is shaped by halo mass quenching (likely through shock virial heating and AGN feedback), but group richness also plays an important role: central galaxies living in less dense environments quenched their star formation later or did not quench it yet. At low masses, processes that delay star formation without invoking too strong supernova-driven outflows could explain the high Ms/Mh ratios of blue centrals as compared to those of the scarce red centrals.

preprint2014arXiv

ELUCID - Exploring the Local Universe with reConstructed Initial Density field I: Hamiltonian Markov Chain Monte Carlo Method with Particle Mesh Dynamics

Simulating the evolution of the local universe is important for studying galaxies and the intergalactic medium in a way free of cosmic variance. Here we present a method to reconstruct the initial linear density field from an input non-linear density field, employing the Hamiltonian Markov Chain Monte Carlo (HMC) algorithm combined with Particle Mesh (PM) dynamics. The HMC+PM method is applied to cosmological simulations, and the reconstructed linear density fields are then evolved to the present day with N-body simulations. The constrained simulations so obtained accurately reproduce both the amplitudes and phases of the input simulations at various $z$. Using a PM model with a grid cell size of 0.75 Mpc/h and 40 time-steps in the HMC can recover more than half of the phase information down to a scale k~0.85 h/Mpc at high z and to k~3.4 h/Mpc at z=0, which represents a significant improvement over similar reconstruction models in the literature, and indicates that our model can reconstruct the formation histories of cosmic structures over a large dynamical range. Adopting PM models with higher spatial and temporal resolutions yields even better reconstructions, suggesting that our method is limited more by the availability of computer resource than by principle. Dynamic models of structure evolution adopted in many earlier investigations can induce non-Gaussianity in the reconstructed linear density field, which in turn can cause large systematic deviations in the predicted halo mass function. Such deviations are greatly reduced or absent in our reconstruction.

preprint2013arXiv

Detection of the large scale alignment of massive galaxies at z~0.6

We report on the detection of the alignment between galaxies and large-scale structure at z~0.6 based on the CMASS galaxy sample from the Baryon Oscillation Spectroscopy Survey data release 9. We use two statistics to quantify the alignment signal: 1) the alignment two-point correlation function which probes the dependence of galaxy clustering at a given separation in redshift space on the projected angle (theta_p) between the orientation of galaxies and the line connecting to other galaxies, and 2) the cos(2theta)-statistic which estimates the average of cos(2theta_p) for all correlated pairs at given separation. We find significant alignment signal out to about 70 Mpc/h in both statistics. Applications of the same statistics to dark matter halos of mass above 10^12 M_sun/h in a large cosmological simulation show similar scale-dependent alignment signals to the observation, but with higher amplitudes at all scales probed. We show that this discrepancy may be partially explained by a misalignment angle between central galaxies and their host halos, though detailed modeling is needed in order to better understand the link between the orientations of galaxies and host halos. In addition, we find systematic trends of the alignment statistics with the stellar mass of the CMASS galaxies, in the sense that more massive galaxies are more strongly aligned with the large-scale structure.

preprint2013arXiv

Photometric properties and luminosity function of nearby massive early-type galaxies

We perform photometric analyses for a bright early-type galaxy (ETG) sample with 2949 galaxies ($M_{\rm r}<-22.5$ mag) in the redshift range of 0.05 to 0.15, drawn from the SDSS DR7 with morphological classification from Galaxy Zoo 1. We measure the Petrosian and isophotal magnitudes, as well as the corresponding half-light radius for each galaxy. We find that for brightest galaxies ($M_{\rm r}<-23$ mag), our Petrosian magnitudes, and isophotal magnitudes to 25 ${\rm mag/arcsec^2}$ and 1\% of the sky brightness are on average 0.16 mag, 0.20 mag, and 0.26 mag brighter than the SDSS Petrosian values, respectively. In the first case the underestimations are caused by overestimations in the sky background by the SDSS PHOTO algorithm, while the latter two are also due to deeper photometry. Similarly, the typical half-light radii ($r_{50}$) measured by the SDSS algorithm are smaller than our measurements. As a result, the bright-end of the $r$-band luminosity function is found to decline more slowly than previous works. Our measured luminosity densities at the bright end are more than one order of magnitude higher than those of Blanton et al. (2003), and the stellar mass densities at $M_{\ast}\sim 5\times10^{11} M_{\odot}$ and $M_{\ast}\sim 10^{12} M_{\odot}$ are a few tenths and a factor of few higher than those of Bernardi et al. (2010). These results may significantly alleviate the tension in the assembly of massive galaxies between observations and predictions of the hierarchical structure formation model.

preprint2013arXiv

Stellar Mass and color dependence of the three-point correlation function of galaxies in the local universe

The three-point correlation function (3PCF) for galaxies provides an opportunity to measure the non-Gaussianity generated from nonlinear structure formation and also probes information about galaxy formation and evolution that is generally not available from the two-point correlation function (2PCF). We measure the 3PCF of the Sloan Digital Sky Survey DR7 main sample galaxies in both redshift and projected spaces on scales up to 40Mpc/h. We explore the dependence of the 3PCF on galaxy stellar mass and color in order to constrain the formation and evolution for galaxies of different properties. The study of the dependence on these properties also helps better constrain the relation between galaxy stellar mass and color and the properties of their hosting dark-matter halos. We focus on the study of the reduced 3PCF, Q, defined as the ratio between the 3PCF and the sum of the products of the 2PCFs. We find a very weak stellar mass dependence of Q in both redshift and projected spaces. On small scales, more massive galaxies tend to have slightly higher amplitudes of Q. The shape dependence of Q is also weak on these small scales, regardless of stellar mass and color. The reduced 3PCF has a strong color dependence for the low-mass galaxies, while no significant dependence on color is found for the high-mass galaxies. Low-mass red galaxies have higher amplitudes and stronger shape dependence of the reduced 3PCF than the blue galaxies, implying that these low-mass red galaxies tend to populate filamentary structures. The linear galaxy bias model fails to interpret the color dependence of Q, emphasizing the importance of a nonvanishing quadratic bias parameter in the correct modeling of the galaxy color distribution.

preprint2012arXiv

Internal kinematics of groups of galaxies in the Sloan Digital Sky Survey data release 7

We present measurements of the velocity dispersion profile (VDP) for galaxy groups in the final data release of the Sloan Digital Sky Survey (SDSS). For groups of given mass we estimate the redshift-space cross-correlation function (CCF) with respect to a reference galaxy sample, xi(r_p, pi), the projected CCF, w_p(r_p), and the real-space CCF, xi(r). The VDP is then extracted from the redshift distortion in xi(r_p, pi), by comparing xi(r_p, pi) with xi(r). We find that the velocity dispersion (VD) within virial radius (R_200) shows a roughly flat profile, with a slight increase at radii below ~0.3 R_200 for high mass systems. The average VD within the virial radius, sigma_v, is a strongly increasing function of central galaxy mass. We apply the same methodology to N-body simulations with the concordance Lambda cold dark matter cosmology but different values of the density fluctuation parameter sigma_8, and we compare the results to the SDSS results. We show that the sigma_v-M_* relation from the data provides stringent constraints on both sigma_8 and sigma_ms, the dispersion in log M_* of central galaxies at fixed halo mass. Our best-fitting model suggests sigma_8 = 0.86 +/- 0.03 and sigma_ms = 0.16 +/- 0.03. The slightly higher value of sigma_8 compared to the WMAP7 result might be due to a smaller matter density parameter assumed in our simulations. Our VD measurements also provide a direct measure of the dark matter halo mass for central galaxies of different luminosities and masses, in good agreement with the results obtained by Mandelbaum et al. (2006) from stacking the gravitational lensing signals of the SDSS galaxies.

preprint2012arXiv

Resolving Subhaloes' Lives with the Hierarchical Bound-Tracing Algorithm

We develop a new code, the Hierarchical Bound-Tracing (HBT for short) code, to find and trace dark matter subhaloes in simulations based on the merger hierarchy of dark matter haloes. Application of this code to a recent benchmark test of finding subhaloes demonstrates that HBT stands as one of the best codes to trace the evolutionary history of subhaloes. The success of the code lies in its careful treatment of the complex physical processes associated with the evolution of subhaloes and in its robust unbinding algorithm with an adaptive source subhalo management. We keep a full record of the merger hierarchy of haloes and subhaloes, and allow growth of satellite subhaloes through accretion from its "satellite-of-satellites", hence allowing mergers among satellites. Local accretion of background mass is omitted, while rebinding of stripped mass is allowed. The justification of these treatments is provided by case studies of the lives of individual subhaloes and by the success in finding the complete subhalo catalogue. We compare our result to other popular subhalo finders and show that HBT is able to well resolve subhaloes in high density environment and keep strict physical track of subhaloes' merger history. This code is fully parallelized and freely available upon request to the authors.

preprint2012arXiv

Stellar mass versus stellar velocity dispersion: which is better for linking galaxies to their dark matter halos?

It was recently suggested that, compared to its stellar mass (M*), the central stellar velocity dispersion (sigma*) of a galaxy might be a better indicator for its host dark matter halo mass. Here we test this hypothesis by estimating the dark matter halo mass for central alaxies in groups as function of M* and sigma*. For this we have estimated the redshift-space cross-correlation function (CCF) between the central galaxies at given M* and sigma* and a reference galaxy sample, from which we determine both the projected CCF, w_p(r_p), and the velocity dispersion profile (VDP) of satellites around the centrals. A halo mass is then obtained from the average velocity dispersion within the virial radius. At fixed M*, we find very weak or no correlation between halo mass and sigma*. In contrast, strong mass dependence is clearly seen even when sigma* is limited to a narrow range. Our results thus firmly demonstrate that the stellar mass of central galaxies is still a good (if not the best) indicator for dark matter halo mass, better than the stellar velocity dispersion. The dependence of galaxy clustering on sigma* fixed M*, as recently discovered by Wake et al. (2012), may be attributed to satellite galaxies, for which the tidal stripping occurring within halos has stronger effect on stellar mass than on central stellar velocity dispersion.

preprint2011arXiv

Autocorrelations of stellar light and mass at z~0 and ~1: From SDSS to DEEP2

We present measurements of projected autocorrelation functions w_p(r_p) for the stellar mass of galaxies and for their light in the U, B and V bands, using data from the third data release of the DEEP2 Galaxy Redshift Survey and the final data release of the Sloan Digital Sky Survey (SDSS). We investigate the clustering bias of stellar mass and light by comparing these to projected autocorrelations of dark matter estimated from the Millennium Simulations (MS) at z=1 and 0.07, the median redshifts of our galaxy samples. All of the autocorrelation and bias functions show systematic trends with spatial scale and waveband which are impressively similar at the two redshifts. This shows that the well-established environmental dependence of stellar populations in the local Universe is already in place at z=1. The recent MS-based galaxy formation simulation of Guo et al. (2011) reproduces the scale-dependent clustering of luminosity to an accuracy better than 30% in all bands and at both redshifts, but substantially overpredicts mass autocorrelations at separations below about 2 Mpc. Further comparison of the shapes of our stellar mass bias functions with those predicted by the model suggests that both the SDSS and DEEP2 data prefer a fluctuation amplitude of sigma_8 0.8 rather than the sigma_8=0.9 assumed by the MS.

preprint2011arXiv

Galaxy clustering and projected density profiles as traced by satellites in photometric surveys: Methodology and luminosity dependence

We develop a new method which measures the projected density distribution w_p(r_p)n of photometric galaxies surrounding a set of spectroscopically-identified galaxies, and simultaneously the projected correlation function w_p(r_p) between the two populations. In this method we are able to divide the photometric galaxies into subsamples in luminosity intervals when redshift information is unavailable, enabling us to measure w_p(r_p)n and w_p(r_p) as a function of not only the luminosity of the spectroscopic galaxy, but also that of the photometric galaxy. Extensive tests show that our method can measure w_p(r_p) in a statistically unbiased way. The accuracy of the measurement depends on the validity of the assumption in the method that the foreground/background galaxies are randomly distributed and thus uncorrelated with those galaxies of interest. Therefore, our method can be applied to the cases where foreground/background galaxies are distributed in large volumes, which is usually valid in real observations. We applied our method to data from SDSS including a sample of 10^5 LRGs at z~0.4 and a sample of about half a million galaxies at z~0.1, both of which are cross-correlated with a deep photometric sample drawn from the SDSS. On large scales, the relative bias factor of galaxies measured from w_p(r_p) at z~0.4 depends on luminosity in a manner similar to what is found at z~0.1, which are usually probed by autocorrelations of spectroscopic samples. On scales smaller than a few Mpc and at both z~0.4 and z~0.1, the photometric galaxies of different luminosities exhibit similar density profiles around spectroscopic galaxies at fixed luminosity and redshift. This provides clear support for the assumption commonly-adopted in HOD models that satellite galaxies of different luminosities are distributed in a similar way, following the dark matter distribution within their host halos.

preprint2010arXiv

Systematic Effects on Determination of the Growth Factor from Redshift-space Distortions

The linear growth factor of density perturbations is believed to be a powerful observable of future redshift surveys to probe physical properties of dark energy and to distinguish among gravity theories. We investigate systematic effects on determination of the growth factor f from a measurement of redshift-space distortions. Using N-body simulations we identify dark matter halos over a broad mass range. We compute the power spectra and correlation functions for the halos and then examine how well the redshift distortion parameter beta=f/b can be reconstructed as a function of halo mass. We find that beta measured for a fixed halo mass is generally a function of scale even on large scales, in contrast with the common expectation that beta approaches a constant described by Kaiser's formula on such scales. The scale dependence depends on the halo mass, being stronger for smaller halos. It also cannot be easily explained with the well-known distribution function of the halo peculiar velocities. We demonstrate that the biasing for smaller halos has larger nonlinearity and stochasticity, thus the linear bias assumption becomes worse for smaller halos. Only for massive halos with b>1.5, beta approaches the linear theory prediction on scales of r or pi/k>30Mpc/h. Luminous red galaxies (LRG), targeted by the SDSS-III's BOSS survey, tend to reside in very massive halos. Our results indicate that if the LRG is used for the measurement of redshift distortions, f can be measured unbiasedly. On the other hand, if one considers to use emission line galaxies, which are targeted by the BigBOSS survey and inhabited in halos of a broad mass range, the scale dependence of beta must be taken into account carefully; otherwise one might give incorrect constraints on dark energy or modified gravity theories. We also find that beta reconstructed in Fourier space behaves better than that in configuration space.

preprint2010arXiv

The Dependence of the Occupation of Galaxies on the Halo Formation Time

We study the dependence of the galaxy contents within halos on the halo formation time using two galaxy formation models, one being a semianalytic model utilizing the halo assembly history from a high resolution N-body simulation and the other being a smoothed particle hydrodynamics simulation including radiative cooling, star formation, and energy feedback from galactic winds. We confirm the finding by Gao et al. that at fixed mass, the clustering of halos depends on the halo formation time, especially for low-mass halos. This age dependence of halo clustering makes it desirable to study the correlation between the occupation of galaxies within halos and the halo age. We find that, in halos of fixed mass, the number of satellite galaxies has a strong dependence on halo age, with fewer satellites in older halos. The youngest one-third of the halos can have an order of magnitude more satellites than the oldest one-third. For central galaxies, in halos that form earlier, they tend to have more stars and thus appear to be more luminous, and the dependence of their luminosity on halo age is not as strong as that of stellar mass. The results can be understood through the star formation history in halos and the merging of satellites onto central galaxies. The age dependence of the galaxy contents within halos would constitute an important ingredient in a more accurate halo-based model of galaxy clustering.

preprint2009arXiv

A determination of dark matter bispectrum with a large set of N-body simulations

We use a set of numerical N-body simulations to study the large-scale behavior of the reduced bispectrum of dark matter and compare the results with the second-order perturbation theory and the halo models for different halo mass functions. We find that the second-order perturbation theory (PT2) agrees with the simulations fairly well on large scales of k<0.05 h/Mpc, but it shows a signature of deviation as the scale goes down. Even on the largest scale where the bispectrum can be measured reasonably well in our simulations, the inconsistency between PT2 and the simulations appears for the colinear triangle shapes. For the halo model, we find that it can only serve as a qualitative method to help study the behavior of Q on large scales and also on relatively small scales. The failure of second-order perturbation theory will also affect the precise determination of the halo models, since they are connected through the 3-halo term in the halo model. The 2-halo term has too much contribution on the large scales, which is the main reason for the halo model to overpredict the bispectrum on the large scales. Since neither of the models can provide a satisfying description for the bispectrum on scales of about 0.1h/Mpc for the requirement of precision cosmology, we release the reduced bispectrum of dark matter on a large range of scales for future analytical modeling of the bispectrum.

preprint2009arXiv

Accurate universal models for the mass accretion histories and concentrations of dark matter halos

A large amount of observations have constrained cosmological parameters and the initial density fluctuation spectrum to a very high accuracy. However, cosmological parameters change with time and the power index of the power spectrum varies with mass scale dramatically in the so-called concordance Lambda CDM cosmology. Thus, any successful model for its structural evolution should work well simultaneously for various cosmological models and different power spectra. We use a large set of high-resolution N-body simulations of a variety of structure formation models (scale-free, standard CDM, open CDM, and Lambda CDM) to study the mass accretion histories (MAHs), the mass and redshift dependence of concentrations and the concentration evolution histories of dark matter halos. We find that there is significant disagreement between the much-used empirical models in the literature and our simulations. According to two simple but tight correlations we find from the simulation results, we develop new empirical models for both the MAHs and the concentration evolution histories of dark matter halos, and the latter can also be used to predict the mass and redshift dependence of halo concentrations. These models are accurate and universal: the same set of model parameters works well for different cosmological models and for halos of different masses at different redshifts and the model predictions are highly accurate even when the histories are traced to very high redshift. These models are also simple and easy to implement. A web calculator and a user-friendly code to make the relevant calculations are available from http://www.shao.ac.cn/dhzhao/mandc.html . We explain why Lambda CDM halos on nearly all mass scales show two distinct phases in their evolution histories.

preprint2009arXiv

Halo Occupation Distribution Modeling of Clustering of Luminous Red Galaxies

We perform Halo Occupation Distribution (HOD) modeling to interpret small-scale and intermediate-scale clustering of 35,000 luminous early-type galaxies and their cross-correlation with a reference imaging sample of normal L* galaxies in the Sloan Digital Sky Survey. The modeling results show that most of these luminous red galaxies (LRGs) are central galaxies residing in massive halos of typical mass M ~ a few times 10^13 to 10^14 Msun/h, while a few percent of them have to be satellites within halos in order to produce the strong auto-correlations exhibited on smaller scales. The mean luminosity Lc of central LRGs increases with the host halo mass, with a rough scaling relation of Lc \propto M^0.5. The halo mass required to host on average one satellite LRG above a luminosity threshold is found to be about 10 times higher than that required to host a central LRG above the same threshold. We find that in massive halos the distribution of L* galaxies roughly follows that of the dark matter and their mean occupation number scales with halo mass as M^1.5. The HOD modeling results also allows for an intuitive understanding of the scale-dependent luminosity dependence of the cross-correlation between LRGs and L_* galaxies. Constraints on the LRG HOD provide tests to models of formation and evolution of massive galaxies, and they are also useful for cosmological parameter investigations. In one of the appendices, we provide LRG HOD parameters with dependence on cosmology inferred from modeling the two-point auto-correlation functions of LRGs.

preprint2009arXiv

Impacts of Source Properties on Strong Lensing by Rich Galaxy Clusters

We use a high-resolution $N$-body simulation to investigate the influence of background galaxy properties, including redshift, size, shape and clustering, on the efficiency of forming giant arcs by gravitational lensing of rich galaxy clusters. Two large sets of ray-tracing simulations are carried out for 10 massive clusters at two redshifts, i.e. $z_{\rm l} \sim 0.2$ and 0.3. The virial mass ($M_{\rm vir}$) of the simulated lens clusters at $z\sim0.2$ ranges from $6.8\times10^{14} h^{-1} {M_{\odot}}$ to $1.1\times 10^{15} h^{-1} M_{\odot}$. The information of background galaxies brighter than 25 magnitude in the $I$-band is taken from Cosmological Evolution Survey (COSMOS) imaging data. Around $1.7\times 10^5$ strong lensing realizations with these images as background galaxies have been performed for each set. We find that the efficiency for forming giant arcs for $z_{\rm l}=0.2$ clusters is broadly consistent with observations. The efficiency of producing giant arcs by rich clusters is weakly dependent on the source size and clustering. Our principal finding is that a small proportion ($\sim 1/3$) of galaxies with elongated shapes (e.g. ellipticity $ε=1-b/a>0.5$) can boost the number of giant arcs substantially. Compared with recent studies where a uniform ellipticity distribution from 0 to 0.5 is used for the sources, the adoption of directly observed shape distribution increases the number of giant arcs by a factor of $\sim2$. Our results indicate that it is necessary to account for source information and survey parameters (such as point-spread-function, seeing) to make correct predictions of giant arcs and further to constrain the cosmological parameters.(abridged)

preprint2009arXiv

Influence of baryonic physics on the merger time-scale of galaxies in N-body/hydrodynamical simulations

Following our previous work(Jiang et al.(2008)), in which we studied the merger time-scale of galaxies in a high-resolution cosmological hydro/N-body simulation, we investigate the potential influence of uncertainties in the numerical implementation of baryonic physics on the merger time-scale. The simulation used in the previous work suffers from the overcooling problem which causes the central galaxies of large halos too massive. This may result in a shortened merger time-scale compared to that in the real universe. We run a similar simulation, but the stellar mass is significantly reduced to model another extreme case of low stellar mass. Our result shows that the merger time-scale is little affected by the star formation recipes, except for the satellites in nearly radial orbits which show a 22 percent higher time-scale in the lower stellar mass case. Since the radial orbits only account for a small part of the satellites' orbits, the fitting formula in Jiang et al.(2008) is still applicable to a reasonable accuracy, nearly immune to the uncertainty in the baryonic physics.

preprint2009arXiv

Intrinsic Ellipticity Correlation of SDSS Luminous Red Galaxies and Misalignment with their Host Dark Matter Halos

We investigate the orientation correlation of giant elliptical galaxies by measuring the intrinsic ellipticity correlation function of 83,773 luminous red galaxies (LRGs) at redshifts 0.16 -- 0.47 from the Sloan Digital Sky Survey. We have accurately determined the correlation up to 30 $h^{-1}$Mpc. Luminosity dependence of the ellipticity correlation is also detected although the error bars are large, while no evidence is found for its redshift evolution between $z=0.2$ and $z=0.4$. Then we use a cosmological $N$-body simulation to examine misalignment between the central LRGs and their parent dark matter halos. Central and satellite galaxies are assigned to simulated halos by employing a halo occupation distribution model for the LRGs. The ellipticity correlation is predicted to have the same shape as but an amplitude about 4 times higher than our observation if the central LRGs are perfectly aligned with their host halos. This indicates that the central LRG galaxies are preferentially but not perfectly aligned with their host halos. With the assumption that there is a misalignment angle between a central LRG and its host halo which follows a Gaussian distribution with a zero mean and a width $σ_θ$, we obtain a tight constraint on the misalignment parameter, $σ_θ={35.4}^{+4.0}_{-3.3}$ deg. This type of intrinsic ellipticity correlation, if not corrected, can lead to contamination at 5% level to the shear power spectrum in weak lensing surveys of limiting magnitude $R_{AB}=24.5$ if the source central galaxies follow the same misalignment distribution as the LRGs.

preprint2009arXiv

Modelling galaxy stellar mass evolution from z~0.8 to today

We apply the empirical method built for z=0 in the previous work of Wang et al. to a higher redshift, to link galaxy stellar mass directly with its hosting dark matter halo mass at z~0.8. The relation of the galaxy stellar mass and the host halo mass M_infall is constrained by fitting both the stellar mass function and the correlation functions at different stellar mass intervals of the VVDS observation, where M_infall is the mass of the hosting halo at the time when the galaxy was last the central galaxy. We find that for low mass haloes, their residing central galaxies are less massive at high redshift than those at low redshift. For high mass haloes, central galaxies in these haloes at high redshift are a bit more massive than the galaxies at low redshift. Satellite galaxies are less massive at earlier times, for any given mass of hosting haloes. Fitting both the SDSS and VVDS observations simultaneously, we also propose a unified model of the M_stars-M_infall relation, which describes the evolution of central galaxy mass as a function of time. The stellar mass of a satellite galaxy is determined by the same M_stars-M_infall relation of central galaxies at the time when the galaxy is accreted. With these models, we study the amount of galaxy stellar mass increased from z~0.8 to the present day through galaxy mergers and star formation. Low mass galaxies gain their stellar masses from z~0.8 to z=0 mainly through star formation. For galaxies of higher mass, the increase of stellar mass solely through mergers from z=0.8 can make the massive galaxies a factor ~2 larger than observed at z=0. We can also predict stellar mass functions of redshifts up to z~3, and the results are consistent with the latest observations.

preprint2009arXiv

The distribution of ejected subhalos and its implication for halo assembly bias

Using a high-resolution cosmological $N$-body simulation, we identify the ejected population of subhalos, which are halos at redshift $z=0$ but were once contained in more massive `host' halos at high redshifts. The fraction of the ejected subhalos in the total halo population of the same mass ranges from 9% to 4% for halo masses from $\sim 10^{11}$ to $\sim 10^{12}\msun$. Most of the ejected subhalos are distributed within 4 times the virial radius of their hosts. These ejected subhalos have distinct velocity distribution around their hosts in comparison to normal halos. The number of subhalos ejected from a host of given mass increases with the assembly redshift of the host. Ejected subhalos in general reside in high-density regions, and have a much higher bias parameter than normal halos of the same mass. They also have earlier assembly times, so that they contribute to the assembly bias of dark matter halos seen in cosmological simulations. However, the assembly bias is {\it not} dominated by the ejected population, indicating that large-scale environmental effects on normal halos are the main source for the assembly bias.

preprint2006arXiv

Bispectrum and Nonlinear Biasing of Galaxies: Perturbation Analysis, Numerical Simulation and SDSS Galaxy Clustering

We consider nonlinear biasing models of galaxies with particular attention to a correlation between linear and quadratic biasing coefficients, b_1 and b_2. We first derive perturbative expressions for b_1 and b_2 in halo and peak biasing models. Then we compute power spectra and bispectra of dark matter particles and halos using N-body simulation data and of volume-limited subsamples of Sloan Digital Sky Survey (SDSS) galaxies, and determine their b_1 and b_2. We find that the values of those coefficients at linear regimes (k<0.2h/Mpc) are fairly insensitive to the redshift-space distortion and the survey volume shape. The resulting normalized amplitudes of bispectra, Q, for equilateral triangles, are insensitive to the values of b_1 implying that b_2 indeed correlates with b_1. The present results explain the previous finding of Kayo et al. (2004) for the hierarchical relation of three-point correlation functions of SDSS galaxies. While the relations between b_1 and b_2 are quantitatively different for specific biasing models, their approximately similar correlations indicate a fairly generic outcome of the biasing due to the gravity in primordial Gaussian density fields.

preprint2006arXiv

Smoothing Algorithms and High-order Singularities in Gravitational Lensing

We propose a new smoothing method for obtaining surface densities from discrete particle positions from numerical simulations. This is an essential step for many applications in gravitational lensing. This method is based on the ``scatter'' interpretation of the discrete density field in the Smoothed Particle Hydrodynamics. We use Monte Carlo simulations of uniform density fields and one isothermal ellipsoid to empirically derive the noise properties, and best smoothing parameters (such as the number of nearest neighbors used). A cluster from high-resolution simulations is then used to assess the reality of high-order singularities such as swallowtails and butterflies in caustics, which are important for the interpretation of substructures in gravitational lenses. We also compare our method with the Delaunay tesselation field estimator using the galaxy studied by Bradac et al. (2004), and find good agreements. We show that higher order singularities are not only connected with bound subhaloes but also with the satellite streams. However, the presence of high-order singularities are sensitive to not only the fluctuation amplitude of the surface density, but also the detailed form of the underlying smooth lensing potential (such as ellipticity and external shear).

preprint2002arXiv

Warm Dark Matter Model of Galaxy Formation

Cold Dark Matter (CDM) models of galaxy formation had been remarkably successful to explain a number of observations in the past decade. However, with both the theoretical modeling and the observations being improved, CDM models have been very recently shown to have excessive clustering on the sub-galactic scale. Here I discuss a solution, based on our high-resolution numerical simulations, to this outstanding problem by considering Warm Dark Matter (WDM). Our results show that the over-clustering problem on sub-galactic scales can be overcome by WDM models, and all the advantages of CDM models are preserved by WDM models. Therefore, the WDM model will become an interesting alternative to the well-studied CDM models.

preprint2001arXiv

The Fundamental Plane of Spiral Galaxies: Search from Observational Data

The fundamental plane of spiral galaxies was searched from observational data, which can be approximately represented in observational parameters by $ L \propto V^2 R $, where $L$, $R$, and $V$ are the luminosity, the linear size of galactic disk, and the rotation velocity. This plane exists at all optical bands in our samples of more than 500 spiral galaxies in total. It is more fundamental than the relationship between any two of the parameters, and can reduce the residual of the Tully--Fisher $L$--$V$ relations by about 50\%. Noticed that the power index of $V$ is doubled from that of $R$, which implies that the total mass and the mass distribution of galaxies plays an important role in forming the fundamental plane. Involving the third parameter of galactic size has a strong physical implication on galaxy formation and the dark--mass distribution.

preprint2000arXiv

Nonlinear Stochastic Biasing of Peaks and Halos: Scale-Dependence, Time-Evolution, and Redshift-Space Distortion from Cosmological N-body Simulations

We quantify the degree of nonlinearity and stochasticity of the clustering of biased objects, using cosmological N-body simulations. Adopting the peaks and the halos as representative biasing models, we focus on the two-point correlation of the biased objects, dark matter and their cross-correlation. Especially, we take account of the effect of redshift-space distortion and attempt to clarify the scale-dependence and the time-dependence by analyzing the biasing factor and the cross-correlation factor. On small scales, stochasticity and nonlinearity become appreciable and strongly object-dependent, especially in redshift space due to the pair-wise velocity dispersion of the biased objects. Nevertheless, an approximation of deterministic linear biasing $δ_{\rm obj}\simeq b_{\rm obj} δ_{\rm mass}$ works reasonably well even in the quasi-linear regime $r > 10 h^{-1}$ Mpc, and linear redshift-space distortion explains the clustering amplitudes in redshift space in this regime.

preprint1996arXiv

Analytical Approximations to the Low-Order Statistics of Dark Matter Distribution

We show that in a hierarchical clustering model the low-order statistics of the density and the peculiar velocity fields can all be modelled semianalytically for a given cosmology and an initial density perturbation power spectrum $P(k)$. We present such models for the two-point correlation function $ξ(r)$, the amplitude $Q$ of the three-point correlation function, the mean pairwise peculiar velocity $< v_{12}(r)> $, the pairwise peculiar velocity dispersion $< v_{12}^2(r)>$, and the one-point peculiar velocity dispersion $< v_1^2 >$. We test our models against results derived from N-body simulations. These models allow us to understand in detail how these statistics depend on $P(k)$ and cosmological parameters. They can also help to interpret, and maybe correct for, sampling effects when these statistics are estimated from observations. The dependence of the small-scale pairwise peculiar velocity dispersion on rich clusters in the sample, for instance, can be studied quantitatively. There are also significant implications for the reconstruction of the cosmic density field from measurements in redshift space.

preprint1996arXiv

High-Order Correlations of Peaks and Halos: a Step toward Understanding Galaxy Biasing

We develop an analytic model for the hierarchical correlation amplitudes S_j(R)= \bxi_j(R)/\bxi_2^{j-1}(R) of density peaks and dark matter halos in the quasi-linear regime. The statistical distribution of density peaks and dark matter halos within the initial density field are determined by the peak formalism and by an extension of the Press-Schechter formalism, respectively. Modifications of these distributions caused by gravitationally induced motions are treated using a spherical collapse model. We test our model against results from a variety of N-body simulations. The model works well for peaks and for halos that are identified earlier than the time when the moments are calculated. Because halos are spatially exclusive at the time of their identification, our model is only qualitatively correct for halos identified at the same time as the moments are calculated. The S_j depend only weakly on the bias parameter b for large b but increase rapidly with decreasing b at b\sim 1. Thus if galaxies are associated with peaks in the initial density field, or with dark halos formed at high redshifts, a measurement of S_j in the quasilinear regime should determine whether galaxies are significantly biased relative to the mass. We use our model to interpret the observed high order correlation functions of galaxies and clusters. We find that if the values of S_j for galaxies are as high as those given by the APM survey, then APM galaxies should not be significantly biased.

preprint1996arXiv

The Correlation Function of Clusters of Galaxies and the Amplitude of Mass Fluctuations in the Universe

We show that if a sample of galaxy clusters is complete above some mass threshold, then hierarchical clustering theories for structure formation predict its autocorrelation function to be determined purely by the cluster abundance and by the spectrum of linear density fluctuations. Thus if the shape of the initial fluctuation spectrum is known, its amplitude $σ_8$ can be estimated directly from the correlation length of a cluster sample in a way which is independent of the value of $Ω_0$. If the cluster mass corresponding to the sample threshold is also known, it provides an independent estimate of the quantity $σ_8Ω_0^{0.6}$. Thus cluster data should allow both $σ_8$ and $Ω_0$ to be determined observationally. We explore these questions using N-body simulations together with a simple but accurate analytical model based on extensions of Press-Schechter theory. Applying our results to currently available data we find that if the linear fluctuation spectrum has a shape similar to that suggested by the APM galaxy survey, then a correlation length $r_0$ in excess of $20\mpch$ for Abell clusters would require $σ_8>1$, while $r_0<15\mpch$ would require $σ_8<0.5$. With conventional estimates of the relevant mass threshold these imply $Ω_0\la 0.3$ and $Ω_0\ga 1$ respectively.

preprint1996arXiv

Three-point correlation function in the quasilinear regime

Using the second-order Eulerian perturbation theory (SEPT), we study the three-point correlation function $ζ$ in the quasilinear regime for the SCDM, LCDM and MDM models, with the interesting result that these three models have distinctive three-point correlation functions. We test this SEPT prediction using a large set of high-resolution N-body simulations. The N-body results show that the SEPT prediction for $ζ$ is not accurate even in the quasilinear regime ($r\la 10 \mpc$), in contrast to previous N-body tests on the skewness. However, similar to the perturbation theory, our N-body results still predict a strong dependence of the three-point correlation on the triangle shape which is observable in the distribution of galaxies if the galaxies trace the distribution of the underlying mass

preprint1995arXiv

The velocity dispersion profiles of clusters of galaxies: a cosmological test and the sampling effect

In this paper we investigate the velocity dispersion profiles of clusters of galaxies for seven cosmological models. One model is the SCDM model, and the others are six low-density models with the density parameter $Ω=0.1$, 0.2 or 0.3 and with or without a cosmological constant $Λ=1-Ω$. We find that the velocity dispersion profiles depend both on $Ω$ and on $Λ$. For $Λ=0$, the profiles are steeper in a lower-$Ω$ model than in a higher-$Ω$ one. The cosmological constant significantly weakens the dependence on $Ω$: the difference in the profile distributions between two flat models is much smaller than that between the two corresponding open models with the same $Ω$. These results in principle can be used to constrain the cosmological parameters when a large sample of the velocity dispersion profiles is available. Motivated by the practical situation that a sample of $\sim 100$ clusters with $\sim 100$ measured redshifts per cluster is still the best sample available in the foreseen future, we examine carefully to what degree the cosmological parameters can be constrained with the velocity dispersion profiles of such a sample of clusters. The limited sampling around clusters and the limited number of clusters seriously degrade the discriminative power of the velocity dispersion profiles among cosmological models. We find that the five models of $Ω\ge 0.2$ cannot be distinguished by this type of observation. Due to the limited sampling, one should be very cautious in extracting information about the density profile and/or the dynamics around single clusters from the diluted velocity dispersion profiles

preprint1995arXiv

Three-point correlation function of galaxy clusters in cosmological models: a strong dependence on triangle shape

In this paper, we use large $\pppm$ N-body simulations to study the three-point correlation function $ζ$ of clusters in two theoretical models. The first model (LCDM) is a low-density flat model of $Ω_0=0.3$, $Λ_0=0.7$ and $h=0.75$, and the second model (PIM) is an Einstein-de-Sitter model with its linear power spectrum obtained from observations. We found that the scaled function $Q(r,u,v)$, which is defined as the ratio of $ζ(r, ru, ru+rv)$ to the hierarchical sum $ξ(r)ξ(ru)+ ξ(ru) ξ(ru+rv) +ξ(ru+rv)ξ(r)$ (where $ξ$ is the two-point correlation function of clusters), depends weakly on $r$ and $u$, but very strongly on $v$. $Q(r,u,v)$ is about 0.2 at $v=0.1$ and 1.8 at $v=0.9$. A model of $Q(r,u,v)=Θ10^{1.3v^2}$ can fit the data of $ζ$ very nicely with $Θ\approx 0.14$. This model is found to be universal for the LCDM clusters and for the PIM clusters. Furthermore, $Q(r,u,v)$ is found to be insensitive to the cluster richness. We have compared our N-body results with simple analytical theories of cluster formation, like the peak theories or the local maxima theories. We found that these theories do not provide an adequate description for the three-point function of clusters. We have also examined the observational data of $ζ$ presently available, and have not found any contradiction between the observations and our model predictions. The $v$-dependence of $q$ in a projected catalogue of clusters is shown to be much weaker than the $v$-dependence of $Q$ found in the three-dimensional case. It would be important to search for the $v$-dependence of $Q$ in redshift samples of rich clusters.

preprint1994arXiv

Substructures and Density Profiles of Clusters in Models of Galaxy Formation

In this paper we investigate, using high resolution N-body simulations, the density profiles and the morphologies of galaxy clusters in seven models of structure formation. We show that these properties of clusters are closely related to the occurrence of a significant merging event in the recent past. The seven models are: (1) the standard CDM model (SCDM) with $Ω_0 = 1$, $Λ_0=0$ and $h=0.5$; (2) a low-density flat model (FL03) with $Ω_0=0.3$, $Λ_0=0.7$ and $h=0.75$; (3) an open model (OP03) with $Ω_0=0.3$, $Λ_0=0$ and $h=0.75$; (4) a low-density flat model (FL02) with $Ω_0=0.2$, $Λ_0=0.8$ and $h=1$; (5) an open model (OP02) with $Ω_0=0.2$, $Λ_0=0$ and $h=1$; (6) a low-density flat model (FL01) with $Ω_0=0.1$ and $Λ_0=0.9$; (7) an open model (OP01) with $Ω_0=0.1$ and $Λ_0=0$. We find that the density profiles and morphologies of clusters depend both on $Ω_0$ and on $Λ_0$. For $Λ_0=0$, these properties are a monotonic function of $Ω_0$. Clusters in OP01 have the steepest density profiles, their density contours are the roundest and show the smallest center shifts. The other extreme case is SCDM, where clusters show the least steep density profiles and the most elongated contours. For a given $Ω_0$ ($<1$), clusters in the flat model (i.e. with $Λ_0=1-Ω_0$) have flatter density profiles and less substructures than in the corresponding open model. In particular, our results show that low-density flat models with $Ω_0\sim 0.3$, which are currently considered as a successful alternative to SCDM, can produce a substantial fraction of clusters with substructures. This is in contrast to the conception that this kind of models may have serious problem in this aspect.

preprint1993arXiv

The large-scale structure in a universe dominated by cold plus hot dark matter

Using numerical simulations, we investigate the large-scale gravitational clustering in a flat universe dominated by cold plus hot dark matter (i.e., $Ω_0=\ocdm+\ohdm+\obaryon=1$). Primordial density fluctuation spectrum is taken to have the Zel'dovich-Harrison form. Three models are studied, with Model I having $\ocdm=0.69$, $\obaryon=0.01$, and $\ohdm=0.30$ in one flavor of neutrinos; Model II having $\ocdm=0.60$, $\obaryon=0.10$, and $\ohdm=0.30$ in one flavor of neutrinos; Model III having $\ocdm=0.69$, $\obaryon=0.01$, and $\ohdm=0.30$ in three flavors of neutrinos. The initial density spectra are normalized by the COBE quadrupole measurement, and galaxies are identified from the peaks of initial density fields above a certain threshold chosen, to match the observed two-point correlation on scales $\ls 10\mpc$. Thus the clustering properties of both the mass and the galaxies are completely specified. The biasing parameter (for the `galaxies') determined in this way is $b_g\approx 1.2$ for Model I, 1.5 for Model II and 1.6 for Model III. The clustering and motions of the simulated `galaxies' are compared with recent observations. The spatial distributions of galaxies in the hybrid models are very frothy; filaments, sheets, voids etc. of sizes 10 -- 50$\mpc$ are frequently seen in the simulations.