Search

Using reconstructed initial conditions in the SDSS survey volume, we carry out constrained hydrodynamic simulations in three regions representing different types of the cosmic web: the Coma cluster of galaxies; the SDSS great wall; and a large low-density region at $z\sim 0.05$. These simulations, which include star formation and stellar feedback but no AGN formation and feedback, are used to investigate the properties and evolution of intergalactic and intra-cluster media. About half of the warm-hot intergalactic gas is associated with filaments in the local cosmic web. Gas in the outskirts of massive filaments and halos can be heated significantly by accretion shocks generated by mergers of filaments and halos, respectively, and there is a tight correlation between gas temperature and the strength of the local tidal field. The simulations also predict some discontinuities associated with shock fronts and contact edges, which can be tested using observations of the thermal SZ effect and X-rays. A large fraction of the sky is covered by Ly$\alpha$ and OVI absorption systems, and most of the OVI systems and low-column density HI systems are associated with filaments in the cosmic web. The constrained simulations, which follow the formation and heating history of the observed cosmic web, provide an important avenue to interpret observational data. With full information about the origin and location of the cosmic gas to be observed, such simulations can also be used to develop observational strategies., Comment: 21 pages, 15 figures, 1 table

Using the extended halo-based group finder developed by Yang et al. (2021), which is able to deal with galaxies via spectroscopic and photometric redshifts simultaneously, we construct galaxy group and candidate protocluster catalogs in a wide redshift range ($0 < z < 6$) from the joint CFHT Large Area $U$-band Deep Survey (CLAUDS) and Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) deep data set. Based on a selection of 5,607,052 galaxies with $i$-band magnitude $m_{i} < 26$ and a sky coverage of $34.41\ {\rm deg}^2$, we identify a total of 2,232,134 groups, within which 402,947 groups have at least three member galaxies. We have visually checked and discussed the general properties of those richest groups at redshift $z>2.0$. By checking the galaxy number distributions within a $5-7\ h^{-1}\mathrm{Mpc}$ projected separation and a redshift difference $Δz \le 0.1$ around those richest groups at redshift $z>2$, we identified a list of 761, 343 and 43 protocluster candidates in the redshift bins $2\leq z, 25 pages, 16 figures, 2 tables, accepted for publication in ApJ

Low-density points (LDPs), obtained by removing high-density regions of observed galaxies, can trace the large-scale structures (LSSs) of the universe. In particular, it offers an intriguing opportunity to detect weak gravitational lensing from low-density regions. In this work, we investigate the tomographic cross-correlation between Planck cosmic microwave background (CMB) lensing maps and LDP-traced LSSs, where LDPs are constructed from the DR8 data release of the DESI legacy imaging survey, with about 106–107 galaxies. We find that, due to the large sky coverage (20,000 deg2) and large redshift depth (z ≤ 1.2), a significant detection (10σ–30σ) of the CMB lensing–LDP cross-correlation in all six redshift bins can be achieved, with a total significance of ∼53σ over ℓ ≤ 1024. Moreover, the measurements are in good agreement with a theoretical template constructed from our numerical simulation in the WMAP 9 yr ΛCDM cosmology. A scaling factor for the lensing amplitude A lens is constrained to A lens = 1 ± 0.12 for z < 0.2, A lens = 1.07 ± 0.07 for 0.2 < z < 0.4, and A lens = 1.07 ± 0.05 for 0.4 < z < 0.6, with the r-band absolute magnitude cut of −21.5 for LDP selection. A variety of tests have been performed to check the detection reliability against variations in LDP samples and galaxy magnitude cuts, masks, CMB lensing maps, multipole ℓ cuts, sky regions, and photo-z bias. We also perform a cross-correlation measurement between CMB lensing and galaxy number density, which is consistent with the CMB lensing–LDP cross-correlation. This work therefore further convincingly demonstrates that LDP is a competitive tracer of LSS.

The 21 cm intensity mapping experiments promise to obtain the large-scale distribution of HI gas at the post-reionization epoch. In order to reveal the underlying matter density fluctuations from the HI mapping, it is important to understand how HI gas traces the matter density distribution. Both nonlinear halo clustering and nonlinear effects modulating HI gas in halos may determine the scale below which the HI bias deviates from linearity. We employ three approaches to generate the mock HI density from a large-scale N-body simulation at low redshifts, and demonstrate that the assumption of HI linearity is valid at the scale corresponding to the first peak of baryon acoustic oscillations, but breaks down at $k \gtrsim 0.1\,h\, {\rm Mpc}^{-1}$. The nonlinear effects of halo clustering and HI content modulation counteract each other at small scales, and their competition results in a model-dependent "sweet-spot" redshift near $z$=1 where the HI bias is scale-independent down to small scales. We also find that the linear HI bias scales approximately linearly with redshift for $z\le 3$., Comment: 9 pages, 6 figures. Match the accepted version for publication at the ApJ

In the CDM paradigm, the halo mass function is a sensitive probe of the cosmic structure. In observations, halo mass is typically estimated from its relation with other observables. The resulting halo mass function is subject to systematic bias, such as the Eddington bias, due to the scatter or uncertainty in the observable - mass relation. Exact correction for the bias is not easy, as predictions for the observables are typically model-dependent in simulations. In this paper, we point out an interesting feature in the halo mass function of the concordence $\Lambda$CDM model: the total halo mass within each evenly-spaced logarithmic mass bin is approximately the same over a large mass range. We show that this property allows us to construct an almost bias-free halo mass function using only an observable (as a halo mass estimator) and stacked weak lensing measurements as long as the scatter between the true halo mass and the observable-inferred mass has a stable form in logarithmic units. The method is not sensitive to the form of the mass-observable relation. We test the idea using cosmological simulations, and show that the method performs very well for realistic observables., Comment: 10 pages,16 figures

We study the topology of cosmic large-scale structure through the genus statistics, using galaxy catalogues generated from the Millennium Simulation and observational data from the latest Sloan Digital Sky Survey Data Release (SDSS DR7). We introduce a new method for constructing galaxy density fields and for measuring the genus statistics of its isodensity surfaces. It is based on a Delaunay tessellation field estimation (DTFE) technique that allows the definition of a piece-wise continuous density field and the exact computation of the topology of its polygonal isodensity contours, without introducing any free numerical parameter. Besides this new approach, we also employ the traditional approaches of smoothing the galaxy distribution with a Gaussian of fixed width, or by adaptively smoothing with a kernel that encloses a constant number of neighboring galaxies. Our results show that the Delaunay-based method extracts the largest amount of topological information. Unlike the traditional approach for genus statistics, it is able to discriminate between the different theoretical galaxy catalogues analyzed here, both in real space and in redshift space, even though they are based on the same underlying simulation model. In particular, the DTFE approach detects with high confidence a discrepancy of one of the semi-analytic models studied here compared with the SDSS data, while the other models are found to be consistent., 14 pages, 9 figures, accepted by ApJ

As we demonstrated in Paper I, the quenched fractions of central and satellite galaxies as function of halo mass are extremely similar, as long as one controls for stellar mass. The same holds for the quenched fractions as a function of central velocity dispersion, which is tightly correlated with black hole mass, as long as one controls for both stellar and halo mass. Here we use mock galaxy catalogs constructed from the latest semi-analytic model, L-GALAXIES, and the state-of-the-art hydrodynamical simulation, EAGLE, to investigate whether these models can reproduce the trends seen in the data. We also check how the group finder used to identify centrals and satellites impacts our results. We find that L-GALAXIES fails to reproduce the trends. The predicted quenched fraction of central galaxies increases sharply with halo mass around $10^{12.5}h^{-1}M_{\odot}$ and with black hole mass around $\sim10^{6.5}M_{\odot}$, while the predicted quenched fraction of satellites increases with both halo and black hole masses gradually. In contrast, centrals and satellites in EAGLE follow almost the same trend as seen in the data. We discuss the implications of our results for how feedback processes regulate galaxy quenching., Comment: 19 pages, 12 figures; accepted for publication in ApJ

As the second paper of a series on studying galaxy-galaxy lensing signals using the Sloan Digital Sky Survey Data Release 7 (SDSS DR7), we present our measurement and modelling of the lensing signals around groups of galaxies. We divide the groups into four halo mass bins, and measure the signals around four different halo-center tracers: brightest central galaxy (BCG), luminosity-weighted center, number-weighted center and X-ray peak position. For X-ray and SDSS DR7 cross identified groups, we further split the groups into low and high X-ray emission subsamples, both of which are assigned with two halo-center tracers, BCGs and X-ray peak positions. The galaxy-galaxy lensing signals show that BCGs, among the four candidates, are the best halo-center tracers. We model the lensing signals using a combination of four contributions: off-centered NFW host halo profile, sub-halo contribution, stellar contribution, and projected 2-halo term. We sample the posterior of 5 parameters i.e., halo mass, concentration, off-centering distance, sub halo mass, and fraction of subhalos via a MCMC package using the galaxy-galaxy lensing signals. After taking into account the sampling effects (e.g. Eddington bias), we found the best fit halo masses obtained from lensing signals are quite consistent with those obtained in the group catalog based on an abundance matching method, except in the lowest mass bin. Subject headings: (cosmology:) gravitational lensing, galaxies: clusters: general, 12 pages, 7 figures, submitted to ApJ

In a recent study, using the distribution of galaxies in the north galactic pole of SDSS DR7 region enclosed in a 500$\mpch$ box, we carried out our ELUCID simulation (Wang et al. 2016, ELUCID III). Here we {\it light} the dark matter halos and subhalos in the reconstructed region in the simulation with galaxies in the SDSS observations using a novel {\it neighborhood} abundance matching method. Before we make use of thus established galaxy-subhalo connections in the ELUCID simulation to evaluate galaxy formation models, we set out to explore the reliability of such a link. For this purpose, we focus on the following a few aspects of galaxies: (1) the central-subhalo luminosity and mass relations; (2) the satellite fraction of galaxies; (3) the conditional luminosity function (CLF) and conditional stellar mass function (CSMF) of galaxies; and (4) the cross correlation functions between galaxies and the dark matter particles, most of which are measured separately for all, red and blue galaxy populations. We find that our neighborhood abundance matching method accurately reproduces the central-subhalo relations, satellite fraction, the CLFs and CSMFs and the biases of galaxies. These features ensure that thus established galaxy-subhalo connections will be very useful in constraining galaxy formation processes. And we provide some suggestions on the three levels of using the galaxy-subhalo pairs for galaxy formation constraints. The galaxy-subhalo links and the subhalo merger trees in the SDSS DR7 region extracted from our ELUCID simulation are available upon request., 18 pages, 13 figures, ApJ accepted

The intrinsic alignment of galaxies is an important systematic effect in weak-lensing surveys, which can affect the derived cosmological parameters. One direct way to distinguish different alignment models and quantify their effects on the measurement is to produce mocked weak-lensing surveys. In this work, we use full-sky ray-tracing technique to produce mock images of galaxies from the ELUCID $N$-body simulation run with the WMAP9 cosmology. In our model we assume that the shape of central elliptical galaxy follows that of the dark matter halo, and spiral galaxy follows the halo spin. Using the mocked galaxy images, a combination of galaxy intrinsic shape and the gravitational shear, we compare the predicted tomographic shear correlations to the results of KiDS and DLS. It is found that our predictions stay between the KiDS and DLS results. We rule out a model in which the satellite galaxies are radially aligned with the center galaxy, otherwise the shear-correlations on small scales are too high. Most important, we find that although the intrinsic alignment of spiral galaxies is very weak, they induce a positive correlation between the gravitational shear signal and the intrinsic galaxy orientation (GI). This is because the spiral galaxy is tangentially aligned with the nearby large-scale overdensity, contrary to the radial alignment of elliptical galaxy. Our results explain the origin of detected positive GI term from the weak-lensing surveys. We conclude that in future analysis, the GI model must include the dependence on galaxy types in more detail., 23 pages, 13 figures, published in ApJ. Our mock galaxy catalog is available upon request by email to the author (Guoliang@pmo.ac.cn, Kangxi@pmo.ac.cn)

Clusters, filaments, sheets and voids are the building blocks of the cosmic web. In this study, we present and compare two distinct algorithms for finding cosmic filaments and sheets, a task which is far less well established than the identification of dark matter halos or voids. One method is based on the smoothed dark matter density field, the other uses the halo distributions directly. We apply both techniques to one high resolution N-body simulation and reconstruct the filamentary/sheet like network of the dark matter density field. We focus on investigating the properties of the dark matter halos inside these structures, in particular on the directions of their spins and the orientation of their shapes with respect to the directions of the filaments and sheets. We find that both the spin and the major axes of filament-halos with masses, 13 pages, 16 figures, accepted by ApJ

In measuring the power spectrum of the distribution of large numbers of dark matter particles in simulations, or galaxies in observations, one has to use Fast Fourier Transforms (FFT) for calculational efficiency. However, because of the required mass assignment onto grid points in this method, the measured power spectrum $\la |\delta^f(k)|^2\ra$ obtained with an FFT is not the true power spectrum $P(k)$ but instead one that is convolved with a window function $|W(\vec k)|^2$ in Fourier space. In a recent paper, Jing (2005) proposed an elegant algorithm to deconvolve the sampling effects of the window function and to extract the true power spectrum, and tests using N-body simulations show that this algorithm works very well for the three most commonly used mass assignment functions, i.e., the Nearest Grid Point (NGP), the Cloud In Cell (CIC) and the Triangular Shaped Cloud (TSC) methods. In this paper, rather than trying to deconvolve the sampling effects of the window function, we propose to select a particular function in performing the mass assignment that can minimize these effects. An ideal window function should fulfill the following criteria: (i) compact top-hat like support in Fourier space to minimize the sampling effects; (ii) compact support in real space to allow a fast and computationally feasible mass assignment onto grids. We find that the scale functions of Daubechies wavelet transformations are good candidates for such a purpose. Our tests using data from the Millennium Simulation show that the true power spectrum of dark matter can be accurately measured at a level better than 2% up to $k=0.7k_N$, without applying any deconvolution processes. The new scheme is especially valuable for measurements of higher order statistics, e.g. the bi-spectrum,........

We use a modified version of the halo-based group finder developed by Yang et al. to select galaxy groups from the Sloan Digital Sky Survey (SDSS DR4). In the first step, a combination of two methods is used to identify the centers of potential groups and to estimate their characteristic luminosity. Using an iterative approach, the adaptive group finder then uses the average mass-to-light ratios of groups, obtained from the previous iteration, to assign a tentative mass to each group. This mass is then used to estimate the size and velocity dispersion of the underlying halo that hosts the group, which in turn is used to determine group membership in redshift space. Finally, each individual group is assigned two different halo masses: one based on its characteristic luminosity, and the other based on its characteristic stellar mass. Applying the group finder to the SDSS DR4, we obtain 301237 groups in a broad dynamic range, including systems of isolated galaxies. We use detailed mock galaxy catalogues constructed for the SDSS DR4 to test the performance of our group finder in terms of completeness of true members, contamination by interlopers, and accuracy of the assigned masses. This paper is the first in a series and focuses on the selection procedure, tests of the reliability of the group finder, and the basic properties of the group catalogue (e.g. the mass-to-light ratios, the halo mass to stellar mass ratios, etc.). The group catalogues including the membership of the groups are available at http://gax.shao.ac.cn/data/Group.html and http://www.astro.umass.edu/~xhyang/Group.html, 19 pages, 12 figures, Accepted for publication in ApJ. Group catalogues are available at http://gax.shao.ac.cn/data/Group.html and http://www.astro.umass.edu/~xhyang/Group.html

We study the cross-correlation between galaxies of different luminosities and colors, using a sample selected from the SDSS Dr 4. Galaxies are divided into 6 samples according to luminosity, and each of these samples is divided into red and blue subsamples. Projected auto-correlation and cross-correlation is estimated for these subsample. At projected separations r_p > 1\mpch, all correlation functions are roughly parallel, although the correlation amplitude depends systematically on luminosity and color. On r_p < 1\mpch, the auto- and cross-correlation functions of red galaxies are significantly enhanced relative to the corresponding power laws obtained on larger scales. Such enhancement is absent for blue galaxies and in the cross-correlation between red and blue galaxies. We esimate the relative bias factor on scales r > 1\mpch for each subsample using its auto-correlation function and cross-correlation functions. The relative bias factors obtained from different methods are similar. For blue galaxies the luminosity-dependence of the relative bias is strong over the luminosity range probed (-23.0, Comment: 25 pages, 18 figures, Accepted for publication in ApJ

Halo bias is the one of the key ingredients of the halo models. It was shown at a given redshift to be only dependent, to the first order, on the halo mass. In this study, four types of cosmic web environments: clusters, filaments, sheets and voids are defined within a state of the art high resolution $N$-body simulation. Within those environments, we use both halo-dark matter cross-correlation and halo-halo auto correlation functions to probe the clustering properties of halos. The nature of the halo bias differs strongly among the four different cosmic web environments we describe. With respect to the overall population, halos in clusters have significantly lower biases in the {$10^{11.0}\sim 10^{13.5}h^{-1}\rm M_\odot$} mass range. In other environments however, halos show extremely enhanced biases up to a factor 10 in voids for halos of mass {$\sim 10^{12.0}h^{-1}\rm M_\odot$}. Such a strong cosmic web environment dependence in the halo bias may play an important role in future cosmological and galaxy formation studies. Within this cosmic web framework, the age dependency of halo bias is found to be only significant in clusters and filaments for relatively small halos $\la 10^{12.5}\msunh$., Comment: 14 pages, 14 figures, ApJ accepted

In this paper we use high-resolution cosmological simulations to study halo intrinsic alignment and its dependence on mass, formation time and large-scale environment. In agreement with previous studies using N-body simulations, it is found that massive halos have stronger alignment. For given mass, older halos have stronger alignment than younger ones. By identifying the cosmic environment of halo using Hessian matrix, we find that for given mass, halos in cluster regions also have stronger alignment than those in filament. The existing theory has not addressed these dependencies explicitly. In this work we extend the linear alignment model with inclusion of halo bias and find that the halo alignment with its mass and formation time dependence can be explained by halo bias. However, the model can not account for the environment dependence, as it is found that halo bias is lower in cluster and higher in filament. Our results suggest that halo bias and environment are independent factors in determining halo alignment. We also study the halo alignment correlation function and find that halos are strongly clustered along their major axes and less clustered along the minor axes. The correlated halo alignment can extend to scale as large as $100h^{-1}$Mpc where its feature is mainly driven by the baryon acoustic oscillation effect., Comment: 11 pages, 10 figures, ApJ Accepted

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., submitted to ApJ, 19 pages, 22 figures. Please download the high-resolution version at http://staff.ustc.edu.cn/~whywang/paper/

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, Comment: Accepted for publication in ApJ (20/Nov/2014). 23 pages, 16 figures. This version has substantial changes after Referee's report. Parts of Section 5 from our previous version will be publish in a follow-up paper

Major mergers and interactions between gas-rich galaxies with comparable masses are thought to be the main triggers of starburst. In this work, we study, for a large stellar mass range, the interaction rate of the starburst galaxies in the local universe. We focus independently on central and satellite star forming galaxies extracted from the Sloan Digital Sky Survey. Here the starburst galaxies are selected in the star formation rate (SFR) stellar mass plane with SFR five times larger than the median value found for "star forming" galaxies of the same stellar mass. Through visual inspection of their images together with close companions determined using spectroscopic redshifts, we find that ~50% of the "starburst" populations show evident merger features, i.e., tidal tails, bridges between galaxies, double cores and close companions. In contrast, in the control sample we selected from the normal star forming galaxies, only ~19% of galaxies are associated with evident mergers. The interaction rates may increase by ~5% for the starburst sample and 2% for the control sample if close companions determined using photometric redshifts are considered. The contrast of the merger rate between the two samples strengthens the hypothesis that mergers and interactions are indeed the main causes of starburst., 7 pages, 4 figures, ApJL accepted

Using the self-consistent modeling of the conditional stellar mass functions across cosmic time by Yang et al. (2012), we make model predictions for the star formation histories (SFHs) of {\it central} galaxies in halos of different masses. The model requires the following two key ingredients: (i) mass assembly histories of central and satellite galaxies, and (ii) local observational constraints of the star formation rates of central galaxies as function of halo mass. We obtain a universal fitting formula that describes the (median) SFH of central galaxies as function of halo mass, galaxy stellar mass and redshift. We use this model to make predictions for various aspects of the star formation rates of central galaxies across cosmic time. Our main findings are the following. (1) The specific star formation rate (SSFR) at high $z$ increases rapidly with increasing redshift [$\propto (1+z)^{2.5}$] for halos of a given mass and only slowly with halo mass ($\propto M_h^{0.12}$) at a given $z$, in almost perfect agreement with the specific mass accretion rate of dark matter halos. (2) The ratio between the star formation rate (SFR) in the main-branch progenitor and the final stellar mass of a galaxy peaks roughly at a constant value, $\sim 10^{-9.3} h^2 {\rm yr}^{-1}$, independent of halo mass or the final stellar mass of the galaxy. However, the redshift at which the SFR peaks increases rapidly with halo mass. (3) More than half of the stars in the present-day Universe were formed in halos with $10^{11.1}\msunh < M_h < 10^{12.3}\msunh$ in the redshift range $0.4 < z < 1.9$. (4) ... [abridged], 15 figures, 22 pages, Accepted for publication in ApJ

We study the alignment signal between the distribution of the brightest satellite galaxies (BSGs) and the major axes of their host groups using the Sloan Digital Sky Survey group catalog constructed by Yang et al. After correcting for the effect of group ellipticity, a statistically significant (~5σ) major-axis alignment is detected and the alignment angle is found to be 43.°0 ± 0.°4. More massive and richer groups show a stronger BSG alignment. The BSG alignment around blue brightest central galaxies (BCGs) is slightly stronger than that around red BCGs. Red BSGs have a much stronger major-axis alignment than blue BSGs. Unlike BSGs, other satellites do not show very significant alignment with their group's major axis. We further explore BSG alignment using the semi-analytic model (SAM) constructed by Guo et al. In general, we found good agreement of the model with observations: BSGs in the SAM show a strong major-axis alignment that depends on group mass and richness in the same way as in observations and none of the other satellites exhibit prominent alignment. However, a discrepancy also exists in that the SAM shows a BSG color dependence opposite of that in observations, which is most probably induced by a missing large-scale environment ingredient in the SAM. The combination of two popular scenarios can explain the BSG alignment we detected. First, satellites merged into the group along the surrounding filaments, which are strongly aligned with the major axis of the group. Second, BSGs entered their host group more recently than other satellites, so they have preserved more information about their assembling history and major-axis alignment. In the SAM, we found positive evidence for the second scenario in the fact that BSGs merged into groups statistically more recently than other satellites. We also found that most of the BSGs (80%) were BCGs before they merged into groups and earlier merging BSGs tend to be closer to their BCGs than other BSGs. On the other hand, although the BSG color dependence in observations is the opposite of that in the SAM, it might also indicate that the first scenario is a good explanation.

Analysis of the Pangu N-body simulation validates that the bulk flow of halos follows a Maxwellian distribution which variance is consistent with the prediction of the linear theory of structure formation. We propose that the consistency between the observed bulk velocity and theories should be examined at the effective scale of the radius of a spherical top-hat window function yielding the same smoothed velocity variance in linear theory as the sample window function does. We compared some recently estimated bulk flows from observational samples with the prediction of the \Lambda CDM model we used; some results deviate from expectation at a level of ~ 3\sigma but the discrepancy is not as severe as previously claimed. We show that bulk flow is only weakly correlated with the dipole of the internal mass distribution, the alignment angle between the mass dipole and the bulk flow has a broad distribution peaked at ~ 30-50 deg., and also that the bulk flow shows little dependence on the mass of the halos used in the estimation. In a simulation of box size 1Gpc/h, for a cell of radius 100 Mpc/h the maximal bulk velocity is >500 km/s, dipoles of the environmental mass outside the cell are not tightly aligned with the bulk flow, but are rather located randomly around it with separation angles ~ 20-40 deg. In the fastest cell there is a slightly smaller number of low-mass halos; however halos inside are clustered more strongly at scales > ~ 20 Mpc/h, which might be a significant feature since the correlation between bulk flow and halo clustering actually increases in significance beyond such scales.