Source author record

H. J. Mo

H. J. Mo 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

59works
3topics
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

59 published item(s)

preprint2022arXiv

ELUCID VII: Using Constrained Hydro Simulations to Explore the Gas Component of the Cosmic Web

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$α$ 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.

preprint2022arXiv

Groups and protocluster candidates in the CLAUDS and HSC-SSP joint deep surveys

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<3$, $3\leq z<4$ and $z \geq 4$, respectively. In general, these catalogs of galaxy groups and protocluster candidates will provide useful environmental information in probing galaxy evolution along the cosmic time.

preprint2022arXiv

Massive Star-Forming Galaxies Have Converted Most of Their Halo Gas into Stars

In the local Universe, the efficiency for converting baryonic gas into stars is very low. In dark matter halos where galaxies form and evolve, the average efficiency varies with galaxy stellar mass and has a maximum of about twenty percent for Milky-Way-like galaxies. The low efficiency at higher mass is believed to be produced by some quenching processes, such as the feedback from active galactic nuclei. We perform an analysis of weak lensing and satellite kinematics for SDSS central galaxies. Our results reveal that the efficiency is much higher, more than sixty percent, for a large population of massive star-forming galaxies around $10^{11}M_{\odot}$. This suggests that these galaxies acquired most of the gas in their halos and converted it into stars without being affected significantly by quenching processes. This population of galaxies is not reproduced in current galaxy formation models, indicating that our understanding of galaxy formation is incomplete. The implications of our results on circumgalactic media, star formation quenching and disc galaxy rotation curves are discussed. We also examine systematic uncertainties in halo-mass and stellar-mass measurements that might influence our results.

preprint2020arXiv

Detection of missing baryons in galaxy groups with kinetic Sunyaev-Zel'dovich effect

We present the detection of the kinetic Sunyaev-Zel'dovich effect (kSZE) signals from groups of galaxies as a function of halo mass down to $\log (M_{500}/{\rm M_\odot}) \sim 12.3$, using the {\it Planck} CMB maps and stacking about $40,000$ galaxy systems with known positions, halo masses, and peculiar velocities. The signals from groups of different mass are constrained simultaneously to take care of projection effects of nearby halos. The total kSZE flux within halos estimated implies that the gas fraction in halos is about the universal baryon fraction, even in low-mass halos, indicating that the `missing baryons' are found. Various tests performed show that our results are robust against systematic effects, such as contamination by infrared/radio sources and background variations, beam-size effects and contributions from halo exteriors. Combined with the thermal Sunyaev-Zel'dovich effect, our results indicate that the `missing baryons' associated with galaxy groups are contained in warm-hot media with temperatures between $10^5$ and $10^6\,{\rm K}$.

preprint2020arXiv

Relating the structure of dark matter halos to their assembly and environment

We use a large $N$-body simulation to study the relation of the structural properties of dark matter halos to their assembly history and environment. The complexity of individual halo assembly histories can be well described by a small number of principal components (PCs), which, compared to formation times, provide a more complete description of halo assembly histories and have a stronger correlation with halo structural properties. Using decision trees built with the random ensemble method, we find that about $60\%$, $10\%$, and $20\%$ of the variances in halo concentration, axis ratio, and spin, respectively, can be explained by combining four dominating predictors: the first PC of the assembly history, halo mass, and two environment parameters. Halo concentration is dominated by halo assembly. The local environment is found to be important for the axis ratio and spin but is degenerate with halo assembly. The small percentages of the variance in the axis ratio and spin that are explained by known assembly and environmental factors suggest that the variance is produced by many nuanced factors and should be modeled as such. The relations between halo intrinsic properties and environment are weak compared to their variances, with the anisotropy of the local tidal field having the strongest correlation with halo properties. Our method of dimension reduction and regression can help simplify the characterization of the halo population and clarify the degeneracy among halo properties.

preprint2016arXiv

Alignments of dark matter halos with large-scale tidal fields: mass and redshift dependence

Large scale tidal field estimated directly from the distribution of dark matter halos is used to investigate how halo shapes and spin vectors are aligned with the cosmic web. The major, intermediate and minor axes of halos are aligned with the corresponding tidal axes, and halo spin axes tend to be parallel with the intermediate axes and perpendicular to the major axes of tidal field. The strengths of these alignments generally increase with halo mass and redshift, but the dependencies are only through the peak height, ν. The scaling relations of the alignment strengths with the value of ν indicate that the alignment strengths remain roughly constant when the structures within which the halos reside are still in quasi-linear regime, but decreases as nonlinear evolution becomes more important. We also calculate the alignments in projection so that our results can be compared directly with observations. Finally, we investigate the alignments of tidal tensors on large scales, and use the results to understand alignments of halo pairs separated at various distances. Our results suggest coherent structure of the tidal field is the underlying reason for the alignments of halos and galaxies seen in numerical simulations and in observations.

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.

preprint2016arXiv

Galaxy groups in the 2MASS Redshift Survey

A galaxy group catalog is constructed from the 2MASS Redshift Survey (2MRS) with the use of a halo-based group finder. The halo mass associated with a group is estimated using a `GAP' method based on the luminosity of the central galaxy and its gap with other member galaxies. Tests using mock samples shows that this method is reliable, particularly for poor systems containing only a few members. On average 80% of all the groups have completeness >0.8, and about 65% of the groups have zero contamination. Halo masses are estimated with a typical uncertainty $\sim 0.35\,{\rm dex}$. The application of the group finder to the 2MRS gives 29,904 groups from a total of 43,246 galaxies at $z \leq 0.08$, with 5,286 groups having two or more members. Some basic properties of this group catalog is presented, and comparisons are made with other groups catalogs in overlap regions. With a depth to $z\sim 0.08$ and uniformly covering about 91% of the whole sky, this group catalog provides a useful data base to study galaxies in the local cosmic web, and to reconstruct the mass distribution in the local Universe.

preprint2016arXiv

Mapping the real space distributions of galaxies in SDSS DR7: I. Two Point Correlation Functions

Using a method to correct redshift space distortion (RSD) for individual galaxies, we mapped the real space distributions of galaxies in the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7). We use an ensemble of mock catalogs to demonstrate the reliability of our method. Here as the first paper in a series, we mainly focus on the two point correlation function (2PCF) of galaxies. Overall the 2PCF measured in the reconstructed real space for galaxies brighter than $^{0.1}{\rm M}_r-5\log h=-19.0$ agrees with the direct measurement to an accuracy better than the measurement error due to cosmic variance, if the reconstruction uses the correct cosmology. Applying the method to the SDSS DR7, we construct a real space version of the main galaxy catalog, which contains 396,068 galaxies in the North Galactic Cap with redshifts in the range $0.01 \leq z \leq 0.12$. The Sloan Great Wall, the largest known structure in the nearby Universe, is not as dominant an over-dense structure as appears to be in redshift space. We measure the 2PCFs in reconstructed real space for galaxies of different luminosities and colors. All of them show clear deviations from single power-law forms, and reveal clear transitions from 1-halo to 2-halo terms. A comparison with the corresponding 2PCFs in redshift space nicely demonstrates how RSDs boost the clustering power on large scales (by about $40-50\%$ at scales $\sim 10 h^{-1}{\rm {Mpc}}$) and suppress it on small scales (by about $70-80\%$ at a scale of $0.3 h^{-1}{\rm {Mpc}}$).

preprint2015arXiv

Assessing Colour-dependent Occupation Statistics Inferred from Galaxy Group Catalogues

We investigate the ability of current implementations of galaxy group finders to recover colour-dependent halo occupation statistics. To test the fidelity of group catalogue inferred statistics, we run three different group finders used in the literature over a mock that includes galaxy colours in a realistic manner. Overall, the resulting mock group catalogues are remarkably similar, and most colour-dependent statistics are recovered with reasonable accuracy. However, it is also clear that certain systematic errors arise as a consequence of correlated errors in group membership determination, central/satellite designation, and halo mass assignment. We introduce a new statistic, the halo transition probability (HTP), which captures the combined impact of all these errors. As a rule of thumb, errors tend to equalize the properties of distinct galaxy populations (i.e. red vs. blue galaxies or centrals vs. satellites), and to result in inferred occupation statistics that are more accurate for red galaxies than for blue galaxies. A statistic that is particularly poorly recovered from the group catalogues is the red fraction of central galaxies as function of halo mass. Group finders do a good job in recovering galactic conformity, but also have a tendency to introduce weak conformity when none is present. We conclude that proper inference of colour-dependent statistics from group catalogues is best achieved using forward modelling (i.e., running group finders over mock data), or by implementing a correction scheme based on the HTP, as long as the latter is not too strongly model-dependent.

preprint2015arXiv

Classical bulges, supermassive blackholes and AGN feedback: Extension to low-mass galaxies

The empirical model of Lu et al. 2014a for the relation between star formation rate and halo mass growth is adopted to predict the classical bulge mass ($M_{\rm cb}$) - total stellar mass ($M_\star$) relation for central galaxies. The assumption that the supermassive black hole (SMBH) mass ($M_{\rm BH}$) is directly proportional to the classical bulge mass, with the proportionality given by that for massive galaxies, predicts a $M_{\rm BH}$ - $M_\star$ relation that matches well the observed relation for different types of galaxies. In particular, the model reproduces the strong transition at $M_\star=10^{10.5}$ - $10^{11}M_{\odot}$, below which $M_{\rm BH}$ drops rapidly with decreasing $M_\star$. Our model predicts a new sequence at $M_\star <10^{10.5}M_{\odot}$, where $M_{\rm BH} \propto M_\star$ but the amplitude is a factor of $\sim 50$ lower than the amplitude of the sequence at $M_\star>10^{11}M_{\odot}$. If all SMBH grow through similar quasar modes with a feedback efficiency of a few percent, the energy produced in low-mass galaxies at redshift $z\gtrsim 2$ can heat the circum-galactic medium up to a specific entropy level that is required to prevent excessive star formation in low-mass dark matter halos.

preprint2015arXiv

Galaxy Ecosystems: gas contents, inflows and outflows

We use a set of observational data for galaxy cold gas mass fraction and gas phase metallicity to constrain the content, inflow and outflow of gas in central galaxies hosted by halos with masses between $10^{11} M_{\odot}$ to $10^{12} M_{\odot}$. The gas contents in high redshift galaxies are obtained by combining the empirical star formation histories of Lu et al. (2014) and star formation models that relate star formation rate with the cold gas mass in galaxies. We find that the total baryon mass in low-mass galaxies is always much less than the universal baryon mass fraction since $z = 2$, regardless of star formation model adopted. The data for the evolution of the gas phase metallicity require net metal outflow at $z\lesssim 2$, and the metal loading factor is constrained to be about $0.01$, or about $60\%$ of the metal yield. Based on the assumption that galactic outflow is more enriched in metal than both the interstellar medium and the material ejected at earlier epochs, we are able to put stringent constraints on the upper limits for both the net accretion rate and the net mass outflow rate. The upper limits strongly suggest that the evolution of the gas phase metallicity and gas mass fraction for low-mass galaxies at $z < 2$ is not compatible with strong outflow. We speculate that the low star formation efficiency of low-mass galaxies is owing to some preventative processes that prevent gas from accreting into galaxies in the first place.

preprint2015arXiv

Gas and Metal Contents of Galaxies and Gaseous Halos: Preventive versus Ejective Feedback

Using a semi-analytical approach we investigate the characteristics of predictions for the masses and metallicities of the baryonic matter in and around galaxies made by three galaxy formation models. These models represent three different feedback scenarios: one model with purely ejective feedback, one model with ejective feedback with reincorporation of ejected gas, and one preventative model. We find that, when the model parameters are adjusted to predict the correct stellar masses for a range of halo masses between 10^{10} to 10^{12}Msun, these three scenarios have very different predictions for the masses and metallicities of the interstellar and circum-galactic media. Compared with current observational data, the model implementing preventative feedback has a large freedom to match a broad range of observational data, while the ejective models have difficulties to match a number of observational constraints simultaneously, independent of how the ejection and reincorporation are implemented. Our results suggest that the feedback process which regulates the amounts of stars and cold gas in low-mass galaxies is preventative in nature.

preprint2015arXiv

Star Formation and Quenching of Satellite Galaxies

We study the quenching of satellite galaxies by gradual depletion of gas due to star formation, by ram-pressure striping and by tidally triggered starburst. Using progenitors constrained by the empirical model of Lu et al., in which outflow loading factor is low, we do not find an over-quenching problem in satellites even if there is no further cold gas supply from the cooling of the halo gas after a galaxy is accreted by its host. Gradual depletion alone predicts a unimodal distribution in specific star formation, in contrast to the bimodal distribution observed, and under-predicts the quenched fraction in low mass halos. Ram-pressure stripping nicely reproduces the bimodal distribution but under-predicts the quenched fraction in low-mass halos. Starbursts in gas-rich satellites triggered by tidal interactions with central galaxies can nicely reproduce the quenched satellite population in low-mass halos, but become unimportant for low-mass satellites in massive halos. The combined processes, together with the constrained progenitors, can reproduce the observed star formation properties of satellites in halos of different masses.

preprint2015arXiv

Star Formation and Stellar Mass Assembly in Dark Matter Halos: From Giants to Dwarfs

The empirical model of Lu et al. 2014 is updated with recent data and used to study galaxy star formation and assembly histories. At $z > 2$, the predicted galaxy stellar mass functions are steep, and a significant amount of star formation is hosted by low-mass haloes that may be missed in current observations. Most of the stars in cluster centrals formed earlier than $z\approx 2$ but have been assembled much later. Milky Way mass galaxies have had on-going star formation without significant mergers since $z\approx 2$, and are thus free of significant (classic) bulges produced by major mergers. In massive clusters, stars bound in galaxies and scattered in the halo form a homogeneous population that is old and with solar metallicity. In contrast, in Milky Way mass systems the two components form two distinct populations, with halo stars being older and poorer in metals by a factor of $\approx 3$. Dwarf galaxies in haloes with $M_{\rm h} < 10^{11}h^{-1}M_{\odot}$ have experienced a star formation burst accompanied by major mergers at $z > 2$, followed by a nearly constant star formation rate after $z = 1$. The early burst leaves a significant old stellar population that is distributed in spheroids.

preprint2014arXiv

Bayesian inferences of galaxy formation from the K-band luminosity and HI mass functions of galaxies: constraining star formation and feedback

We infer mechanisms of galaxy formation for a broad family of semi-analytic models (SAMs) constrained by the K-band luminosity function and HI mass function of local galaxies using tools of Bayesian analysis. Even with a broad search in parameter space the whole model family fails to match to constraining data. In the best fitting models, the star formation and feedback parameters in low-mass haloes are tightly constrained by the two data sets, and the analysis reveals several generic failures of models that similarly apply to other existing SAMs. First, based on the assumption that baryon accretion follows the dark matter accretion, large mass-loading factors are required for haloes with circular velocities lower than 200 km/s, and most of the wind mass must be expelled from the haloes. Second, assuming that the feedback is powered by Type-II supernovae with a Chabrier IMF, the outflow requires more than 25% of the available SN kinetic energy. Finally, the posterior predictive distributions for the star formation history are dramatically inconsistent with observations for masses similar to or smaller than the Milky-Way mass. The inferences suggest that the current model family is still missing some key physical processes that regulate the gas accretion and star formation in galaxies with masses below that of the Milky Way.

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.

preprint2014arXiv

Formation of disk galaxies in preheated media: a preventative feedback model

We introduce a semi-analytic galaxy formation model implementing a self-consistent treatment for the hot halo gas configuration and the assembly of central disks. Using the model, we explore a preventative feedback model, in which the circum-halo medium is assumed to be preheated up to a certain entropy level by early starbursts or other processes, and compare it with an ejective feedback model, in which baryons are first accreted into dark matter halos and subsequently ejected out by feedback. The model demonstrates that when the medium is preheated to an entropy comparable to the halo virial entropy the baryon accretion can be largely reduced and delayed. In addition, the preheated medium can establish an extended low density gaseous halo when it accretes into the dark matter halos, and result in a specific angular momentum of the cooling gas large enough to form central disks as extended as those observed. Combined with simulated halo assembly histories, the preventative feedback model can reproduce remarkably well a number of observational scaling relations. These include the cold baryon (stellar plus cold gas) mass fraction-halo mass relations, star formation histories, disk size-stellar mass relation and its evolution, and the number density of low-mass galaxies as a function of redshift. In contrast, the conventional ejective feedback model fails to reproduce these observational trends. Using the model, we demonstrate that the properties of disk galaxies are closely tied to the thermal state of hot halo gas and even possibly the circum-halo medium, which suggests that observational data for the disk properties and circum-galactic hot/warm medium may jointly provide interesting constraints for galaxy formation models.

preprint2014arXiv

Spin alignments of spiral galaxies within the large-scale structure from SDSS DR7

Using a sample of spiral galaxies selected from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7) and Galaxy Zoo 2 (GZ2), we investigate the alignment of spin axes of spiral galaxies with their surrounding large scale structure, which is characterized by the large-scale tidal field reconstructed from the data using galaxy groups above a certain mass threshold. We find that the spin axes of only have weak tendency to be aligned with (or perpendicular to) the intermediate (or minor) axis of the local tidal tensor. The signal is the strongest in a \cluster environment where all the three eigenvalues of the local tidal tensor are positive. Compared to the alignments between halo spins and local tidal field obtained in N-body simulations, the above observational results are in best agreement with those for the spins of inner regions of halos, suggesting that the disk material traces the angular momentum of dark matter halos in the inner regions.

preprint2013arXiv

Alignments of galaxies within cosmic filaments from SDSS DR7

Using a sample of galaxy groups selected from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7), we examine the alignment between the orientation of galaxies and their surrounding large scale structure in the context of the cosmic web. The latter is quantified using the large-scale tidal field, reconstructed from the data using galaxy groups above a certain mass threshold. We find that the major axes of galaxies in filaments tend to be preferentially aligned with the directions of the filaments, while galaxies in sheets have their major axes preferentially aligned parallel to the plane of the sheets. The strength of this alignment signal is strongest for red, central galaxies, and in good agreement with that of dark matter halos in N-body simulations. This suggests that red, central galaxies are well aligned with their host halos, in quantitative agreement with previous studies based on the spatial distribution of satellite galaxies. There is a luminosity and mass dependence that brighter and more massive galaxies in filaments and sheets have stronger alignment signals. We also find that the orientation of galaxies is aligned with the eigenvector associated with the smallest eigenvalue of the tidal tensor. These observational results indicate that galaxy formation is affected by large-scale environments, and strongly suggests that galaxies are aligned with each other over scales comparable to those of sheets and filaments in the cosmic web.

preprint2013arXiv

Constraining the Star Formation Histories in Dark Matter Halos: I. Central Galaxies

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]

preprint2013arXiv

Measuring the X-ray luminosities of SDSS DR7 clusters from RASS

We use ROSAT All Sky Survey (RASS) broadband X-ray images and the optical clusters identified from SDSS DR7 to estimate the X-ray luminosities around $\sim 65,000$ candidate clusters with masses $\ga 10^{13}\msunh$ based on an Optical to X-ray (OTX) code we develop. We obtain a catalogue with X-ray luminosity for each cluster. This catalog contains 817 clusters (473 at redshift $z\le 0.12$) with $S/N> 3$ in X-ray detection. We find about $65\%$ of these X-ray clusters have their most massive member located near the X-ray flux peak; for the rest $35\%$, the most massive galaxy is separated from the X-ray peak, with the separation following a distribution expected from a NFW profile. We investigate a number of correlations between the optical and X-ray properties of these X-ray clusters, and find that: the cluster X-ray luminosity is correlated with the stellar mass (luminosity) of the clusters, as well as with the stellar mass (luminosity) of the central galaxy and the mass of the halo, but the scatter in these correlations is large. Comparing the properties of X-ray clusters of similar halo masses but having different X-ray luminosities, we find that massive halos with masses $\ga 10^{14}\msunh$ contain a larger fraction of red satellite galaxies when they are brighter in X-ray. ... A cluster catalog containing the optical properties of member galaxies and the X-ray luminosity is available at {\it http://gax.shao.ac.cn/data/Group.html}.

preprint2013arXiv

Reconstructing the Initial Density Field of the Local Universe: Method and Test with Mock Catalogs

Our research objective in this paper is to reconstruct an initial linear density field, which follows the multivariate Gaussian distribution with variances given by the linear power spectrum of the current CDM model and evolves through gravitational instability to the present-day density field in the local Universe. For this purpose, we develop a Hamiltonian Markov Chain Monte Carlo method to obtain the linear density field from a posterior probability function that consists of two components: a prior of a Gaussian density field with a given linear spectrum, and a likelihood term that is given by the current density field. The present-day density field can be reconstructed from galaxy groups using the method developed in Wang et al. (2009a). Using a realistic mock SDSS DR7, obtained by populating dark matter haloes in the Millennium simulation with galaxies, we show that our method can effectively and accurately recover both the amplitudes and phases of the initial, linear density field. To examine the accuracy of our method, we use $N$-body simulations to evolve these reconstructed initial conditions to the present day. The resimulated density field thus obtained accurately matches the original density field of the Millennium simulation in the density range 0.3 <= rho/rho_mean <= 20 without any significant bias. Especially, the Fourier phases of the resimulated density fields are tightly correlated with those of the original simulation down to a scale corresponding to a wavenumber of ~ 1 h/Mpc, much smaller than the translinear scale, which corresponds to a wavenumber of ~ 0.15 h\Mpc.

preprint2012arXiv

Bayesian inference of galaxy formation from the K-band luminosity function of galaxies: tensions between theory and observation

We conduct Bayesian model inferences from the observed K-band luminosity function of galaxies in the local Universe, using the semi-analytic model (SAM) of galaxy formation introduced in Lu et al (2011). The prior distributions for the 14 free parameters include a large range of possible models. We find that some of the free parameters, e.g. the characteristic scales for quenching star formation in both high-mass and low-mass halos, are already tightly constrained by the single data set. The posterior distribution includes the model parameters adopted in other SAMs. By marginalising over the posterior distribution, we make predictions that include the full inferential uncertainties for the colour-magnitude relation, the Tully-Fisher relation, the conditional stellar mass function of galaxies in halos of different masses, the HI mass function, the redshift evolution of the stellar mass function of galaxies, and the global star formation history. Using posterior predictive checking with the available observational results, we find that the model family (i) predicts a Tully-Fisher relation that is curved; (ii) significantly over predicts the satellite fraction; (iii) vastly over predicts the HI mass function; (iv) predicts high-z stellar mass functions that have too many low mass galaxies and too few high mass ones. and (v) predicts a redshift evolution of the stellar mass density and the star formation history that are in moderate disagreement. These results suggest that some important processes are still missing in the current model family and we discuss a number of possible solutions to solve the discrepancies, such as interactions between galaxies and dark matter halos, tidal stripping, the bimodal accretion of gas, preheating, and a redshift-dependent initial mass function.

preprint2012arXiv

Cross identification between X-ray and Optical Clusters of Galaxies in the SDSS DR7 Field

We use the ROSAT all sky survey X-ray cluster catalogs and the optical SDSS DR7 galaxy and group catalogs to cross-identify X-ray clusters with their optical counterparts, resulting in a sample of 201 X-ray clusters in the sky coverage of SDSS DR7. We investigate various correlations between the optical and X-ray properties of these X-ray clusters, and find that the following optical properties are correlated with the X-ray luminosity: the central galaxy luminosity, the central galaxy mass, the characteristic group luminosity ($\propto \Lx^{0.43}$), the group stellar mass ($\propto \Lx^{0.46}$), with typical 1-$σ$ scatter of $\sim 0.67$ in $\log \Lx$. Using the observed number distribution of X-ray clusters, we obtain an unbiased scaling relation between the X-ray luminosity, the central galaxy stellar mass and the characteristic satellite stellar mass as ${\log L_X} = -0.26 + 2.90 [\log (M_{\ast, c} + 0.26 M_{\rm sat}) -12.0]$ (and in terms of luminosities, as ${\log L_X} = -0.15 + 2.38 [\log (L_{c} + 0.72 L_{\rm sat}) -12.0]$). We find that the systematic difference between different halo mass estimations, e.g., using the ranking of characteristic group stellar mass or using the X-ray luminosity scaling relation can be used to constrain cosmology. Comparing the properties of groups of similar stellar mass (or optical luminosities) and redshift that are X-ray luminous or under-luminous, we find that X-ray luminous groups have more faint satellite galaxies and higher red fraction in their satellites. The cross-identified X-ray clusters together with their optical properties are provided in Appendix B.

preprint2012arXiv

Evolution of the Galaxy - Dark Matter Connection and the Assembly of Galaxies in Dark Matter Halos

We present a new model to describe the galaxy-dark matter connection across cosmic time, which unlike the popular subhalo abundance matching technique is self-consistent in that it takes account of the facts that (i) subhalos are accreted at different times, and (ii) the properties of satellite galaxies may evolve after accretion. Using observations of galaxy stellar mass functions out to $z \sim 4$, the conditional stellar mass function at $z\sim 0.1$ obtained from SDSS galaxy group catalogues, and the two-point correlation function (2PCF) of galaxies at $z \sim 0.1$ as function of stellar mass, we constrain the relation between galaxies and dark matter halos over the entire cosmic history from $z \sim 4$ to the present. This relation is then used to predict the median assembly histories of different stellar mass components within dark matter halos (central galaxies, satellite galaxies, and halo stars). We also make predictions for the 2PCFs of high-$z$ galaxies as function of stellar mass. Our main findings are the following: (i) Our model reasonably fits all data within the observational uncertainties, indicating that the $Λ$CDM concordance cosmology is consistent with a wide variety of data regarding the galaxy population across cosmic time. (ii) ... [abridged]

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.

preprint2011arXiv

A Bayesian approach to the semi-analytic model of galaxy formation: methodology

We believe that a wide range of physical processes conspire to shape the observed galaxy population but we remain unsure of their detailed interactions. The semi-analytic model (SAM) of galaxy formation uses multi-dimensional parameterisations of the physical processes of galaxy formation and provides a tool to constrain these underlying physical interactions. Because of the high dimensionality, the parametric problem of galaxy formation may be profitably tackled with a Bayesian-inference based approach, which allows one to constrain theory with data in a statistically rigorous way. In this paper we develop a SAM in the framework of Bayesian inference. We show that, with a parallel implementation of an advanced Markov-Chain Monte-Carlo algorithm, it is now possible to rigorously sample the posterior distribution of the high-dimensional parameter space of typical SAMs. As an example, we characterise galaxy formation in the current $Λ$CDM cosmology using the stellar mass function of galaxies as an observational constraint. We find that the posterior probability distribution is both topologically complex and degenerate in some important model parameters, suggesting that thorough explorations of the parameter space are needed to understand the models. We also demonstrate that because of the model degeneracy, adopting a narrow prior strongly restricts the model. Therefore, the inferences based on SAMs are conditional to the model adopted. Using synthetic data to mimic systematic errors in the stellar mass function, we demonstrate that an accurate observational error model is essential to meaningful inference.

preprint2011arXiv

An analytical model for the accretion of dark matter subhalos

An analytical model is developed for the mass function of cold dark matter subhalos at the time of accretion and for the distribution of their accretion times. Our model is based on the model of Zhao et al. (2009) for the median assembly histories of dark matter halos, combined with a simple log-normal distribution to describe the scatter in the main-branch mass at a given time for halos of the same final mass. Our model is simple, and can be used to predict the un-evolved subhalo mass function, the mass function of subhalos accreted at a given time, the accretion-time distribution of subhalos of a given initial mass, and the frequency of major mergers as a function of time. We test our model using high-resolution cosmological $N$-body simulations, and find that our model predictions match the simulation results remarkably well. Finally, we discuss the implications of our model for the evolution of subhalos in their hosts and for the construction of a self-consistent model to link galaxies and dark matter halos at different cosmic times.

preprint2011arXiv

On the algorithms of radiative cooling in semi-analytic models

We study the behaviour of multiple radiative cooling algorithms implemented in seven Semi-Analytic Models (SAMs) of galaxy formation, including a new model we propose in this paper. We use versions of the models without feedback and apply them to dark matter haloes growing in a cosmological context, which have final masses that range from 10^{11}Msun to 10^{14}Msun. First, using simplified smoothly-growing halo models, we demonstrate that the different algorithms predict cooling rates and final cold gas masses that differ by a factor of ~5 for massive haloes (>10^{12}Msun). The algorithms are in better agreement for less massive haloes because they cool efficiently and, therefore, their cooling rates are largely limited by the halo accretion rate. However, for less massive haloes, all the SAMs predict less cooling than corresponding 1D hydrodynamic models. Second, we study the gas accretion history of the central galaxies of dark matter haloes using merger trees. The inclusion of mergers alters the cooling history of haloes by locking up gas in galaxies within small haloes at early times. For realistic halo models, the dispersion in the cold gas mass predicted by the algorithms is 0.5 dex for high mass haloes and 0.1 dex for low mass haloes, while the dispersion in the accretion rate is about two times larger. Comparing to cosmological SPH simulations, we find that most SAMs systematically under-predict the gas accretion rates for low-mass haloes but over-predict the gas accretion rates for massive haloes. Although the models all include both "rapid" and "slow" mode accretion, the transition between the two accretion modes varies between models and also differs from the simulations. Finally, we construct a new model that explicitly incorporates cold halo gas to illustrate that such a class of models can better match the results from cosmological hydrodynamic simulations.

preprint2011arXiv

Probing Hot Gas in Galaxy Groups through the Sunyaev-Zeldovich Effect

We investigate the potential of exploiting the Sunyaev-Zeldovich effect (SZE) to study the properties of hot gas in galaxy groups. It is shown that, with upcoming SZE surveys, one can stack SZE maps around galaxy groups of similar halo masses selected from large galaxy redshift surveys to study the hot gas in halos represented by galaxy groups. We use various models for the hot halo gas to study how the expected SZE signals are affected by gas fraction, equation of state, halo concentration, and cosmology. Comparing the model predictions with the sensitivities expected from the SPT, ACT and Planck surveys shows that a SPT-like survey can provide stringent constraints on the hot gas properties for halos with masses M ~> 10^{13} h^{-1}Msun. We also explore the idea of using the cross correlation between hot gas and galaxies of different luminosity to probe the hot gas in dark matter halos without identifying galaxy groups to represent dark halos. Our results show that, with a galaxy survey as large as the Sloan Digital Sky Survey and with the help of the conditional luminosity function (CLF) model, one can obtain stringent constraints on the hot gas properties in halos with masses down to 10^{13} h^{-1}Msun. Thus, the upcoming SZE surveys should provide a very promising avenue to probe the hot gas in relatively low-mass halos where the majority of L*-galaxies reside.

preprint2011arXiv

Reconstructing the Cosmic Velocity and Tidal Fields with Galaxy Groups Selected from the Sloan Digital Sky Survey

[abridge]Cosmic velocity and tidal fields are important for the understanding of the cosmic web and the environments of galaxies, and can also be used to constrain cosmology. In this paper, we reconstruct these two fields in SDSS volume from dark matter halos represented by galaxy groups. Detailed mock catalogues are used to test the reliability of our method against uncertainties arising from redshift distortions, survey boundaries, and false identifications of groups by our group finder. We find that both the velocity and tidal fields, smoothed on a scale of ~2Mpc/h, can be reliably reconstructed in the inner region (~66%) of the survey volume. The reconstructed tidal field is used to split the cosmic web into clusters, filaments, sheets, and voids, depending on the sign of the eigenvalues of tidal tensor. The reconstructed velocity field nicely shows how the flows are diverging from the centers of voids, and converging onto clusters, while sheets and filaments have flows that are convergent along one and two directions, respectively. We use the reconstructed velocity field and the Zel'dovich approximation to predict the mass density field in the SDSS volume as function of redshift, and find that the mass distribution closely follows the galaxy distribution even on small scales. We find a large-scale bulk flow of about 117km/s in a very large volume, equivalent to a sphere with a radius of ~170Mpc/h, which seems to be produced by the massive structures associated with the SDSS Great Wall. Finally, we discuss potential applications of our reconstruction to study the environmental effects of galaxy formation, to generate initial conditions for simulations of the local Universe, and to constrain cosmological models. The velocity, tidal and density fields in the SDSS volume, specified on a Cartesian grid with a spatial resolution of ~700kpc/h, are available from the authors upon request.

preprint2010arXiv

Satellite Kinematics III: Halo Masses of Central Galaxies in SDSS

We use the kinematics of satellite galaxies that orbit around the central galaxy in a dark matter halo to infer the scaling relations between halo mass and central galaxy properties. Using galaxies from the Sloan Digital Sky Survey, we investigate the halo mass-luminosity relation (MLR) and the halo mass-stellar mass relation (MSR) of central galaxies. In particular, we focus on the dependence of these scaling relations on the colour of the central galaxy. We find that red central galaxies on average occupy more massive haloes than blue central galaxies of the same luminosity. However, at fixed stellar mass there is no appreciable difference in the average halo mass of red and blue centrals, especially for M* $\lsim$ 10^{10.5} h^{-2} Msun. This indicates that stellar mass is a better indicator of halo mass than luminosity. Nevertheless, we find that the scatter in halo masses at fixed stellar mass is non-negligible for both red and blue centrals. It increases as a function of stellar mass for red centrals but shows a fairly constant behaviour for blue centrals. We compare the scaling relations obtained in this paper with results from other independent studies of satellite kinematics, with results from a SDSS galaxy group catalog, from galaxy-galaxy weak lensing measurements, and from subhalo abundance matching studies. Overall, these different techniques yield MLRs and MSRs in fairly good agreement with each other (typically within a factor of two), indicating that we are converging on an accurate and reliable description of the galaxy-dark matter connection. We briefly discuss some of the remaining discrepancies among the various methods.

preprint2010arXiv

The Stellar Mass Components of Galaxies: Comparing Semi-Analytical Models with Observation

We compare the stellar masses of central and satellite galaxies predicted by three independent semianalytical models with observational results obtained from a large galaxy group catalogue constructed from the Sloan Digital Sky Survey. In particular, we compare the stellar mass functions of centrals and satellites, the relation between total stellar mass and halo mass, and the conditional stellar mass functions, which specify the average number of galaxies of stellar mass M_* that reside in a halo of mass M_h. The semi-analytical models only predict the correct stellar masses of central galaxies within a limited mass range and all models fail to reproduce the sharp decline of stellar mass with decreasing halo mass observed at the low mass end. In addition, all models over-predict the number of satellite galaxies by roughly a factor of two. The predicted stellar mass in satellite galaxies can be made to match the data by assuming that a significant fraction of satellite galaxies are tidally stripped and disrupted, giving rise to a population of intra-cluster stars in their host halos. However, the amount of intra-cluster stars thus predicted is too large compared to observation. This suggests that current galaxy formation models still have serious problems in modeling star formation in low-mass halos.

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

Galaxy Groups in the SDSS DR4: III. the luminosity and stellar mass functions

Using a large galaxy group catalogue constructed from the Sloan Digital Sky Survey Data Release 4 (SDSS DR4) with an adaptive halo-based group finder, we investigate the luminosity and stellar mass functions for different populations of galaxies (central versus satellite; red versus blue; and galaxies in groups of different masses) and for groups themselves. The conditional stellar mass function (CSMF), which describes the stellar distribution of galaxies in halos of a given mass for central and satellite galaxies can be well modeled with a log-normal distribution and a modified Schechter form, respectively. On average, there are about 3 times as many central galaxies as satellites. Among the satellite population, there are in general more red galaxies than blue ones. For the central population, the luminosity function is dominated by red galaxies at the massive end, and by blue galaxies at the low mass end. At the very low-mass end ($M_\ast \la 10^9 h^{-2}\Msun$), however, there is a marked increase in the number of red centrals. We speculate that these galaxies are located close to large halos so that their star formation is truncated by the large-scale environments. The stellar-mass function of galaxy groups is well described by a double power law, with a characteristic stellar mass at $\sim 4\times 10^{10}h^{-2}\Msun$. Finally, we use the observed stellar mass function of central galaxies to constrain the stellar mass - halo mass relation for low mass halos, and obtain $M_{\ast, c}\propto M_h^{4.9}$ for $M_h \ll 10^{11} \msunh$.

preprint2009arXiv

Stellar Ages and Metallicities of Central and Satellite Galaxies: Implications for Galaxy Formation and Evolution

Using a large SDSS galaxy group catalogue, we study how the stellar ages and metallicities of central and satellite galaxies depend on stellar mass and halo mass. We find that satellites are older and metal-richer than centrals of the same stellar mass. In addition, the slopes of the age-stellar mass and metallicity-stellar mass relations are found to become shallower in denser environments. This is due to the fact that the average age and metallicity of low mass satellite galaxies increase with the mass of the halo in which they reside. A comparison with the semi-analytical model of Wang et al. (2008) shows that it succesfully reproduces the fact that satellites are older than centrals of the same stellar mass and that the age difference increases with the halo mass of the satellite. This is a consequence of strangulation, which leaves the stellar populations of satellites to evolve passively, while the prolonged star formation activity of centrals keeps their average ages younger. The resulting age offset is larger in more massive environments because their satellites were accreted earlier. The model fails, however, in reproducing the halo mass dependence of the metallicities of low mass satellites, yields metallicity-stellar mass and age-stellar mass relations that are too shallow, and predicts that satellite galaxies have the same metallicities as centrals of the same stellar mass, in disagreement with the data. We argue that these discrepancies are likely to indicate the need to (i) modify the recipes of both supernova feedback and AGN feedback, (ii) use a more realistic description of strangulation, and (iii) include a proper treatment of the tidal stripping, heating and destruction of satellite galaxies. [Abridged]

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.

preprint2009arXiv

The nature of red dwarf galaxies

Using dark matter halos traced by galaxy groups selected from the Sloan Digital Sky Survey Data Release 4, we find that about 1/4 of the faint galaxies ($\rmag >-17.05$, hereafter dwarfs) that are the central galaxies in their own halo are not blue and star forming, as expected in standard models of galaxy formation, but are red. In contrast, this fraction is about 1/2 for dwarf satellite galaxies. Many red dwarf galaxies are physically associated with more massive halos. In total, about $\sim 45$% of red dwarf galaxies reside in massive halos as satellites, while another $\sim 25$% have a spatial distribution that is much more concentrated towards their nearest massive haloes than other dwarf galaxies. We use mock catalogs to show that the reddest population of non-satellite dwarf galaxies are distributed within about 3 times the virial radii of their nearest massive halos. We suggest that this population of dwarf galaxies are hosted by low-mass halos that have passed through their massive neighbors, and that the same environmental effects that cause satellite galaxies to become red are also responsible for the red colors of this population of galaxies. We do not find any significant radial dependence of the population of dwarf galaxies with the highest concentrations, suggesting that the mechanisms operating on these galaxies affect color more than structure. However, over 30% of dwarf galaxies are red and isolated and their origin remains unknown.

preprint2008arXiv

Feedback from galactic stellar bulges and hot gaseous haloes of galaxies

We demonstrate that the feedback from stellar bulges can play an essential role in shaping the halo gas of galaxies with substantial bulge components by conducting 1-D hydrodynamical simulations. The feedback model we consider consists of two distinct phases: 1) an early starburst during the bulge formation and 2) a subsequent long-lasting mass and energy injection from stellar winds of low-mass stars and Type Ia SNe. An energetic outward blastwave is initiated by the starburst and is maintained and enhanced by the long-lasting stellar feedback. For a MW-like galactic bulge, this blastwave sweeps up the halo gas in the proto-galaxy and heats up the surrounding medium to a scale much beyond the virial radius of the halo, thus the accretion of the halo hot gas can be completely stopped. In addition, the long-lasting feedback in the later phase powers a galactic bulge wind that is reverse-shocked at a large radius in the presence of surrounding intergalactic medium and hence maintains a hot gaseous halo. As the mass and energy injection decreases with time, the feedback evolves to a subsonic and quasi-stable outflow, which is enough to prevent halo gas from cooling. The two phases of the feedback thus re-enforce each-other's impact on the gas dynamics. The simulation results demonstrate that the stellar bulge feedback may provide a plausible solution to the long-standing problems in understanding the MW type galaxies, such as the "missing stellar feedback" problem and the "over-cooling" problem. The simulations also show that the properties of the hot gas in the subsonic outflow state depend sensitively on the environment and the formation history of the bulge. This dependence and variance may explain the large dispersion in the X-ray to B-band luminosity ratio of the low $L_X/L_B$ Es.

preprint2008arXiv

The subhalo - satellite connection and the fate of disrupted satellite galaxies

In the standard paradigm, satellite galaxies are believed to be associated with the population of dark matter subhalos. In this paper, we use the conditional stellar mass functions of {\it satellite galaxies} obtained from a large galaxy group catalogue together with models of the subhalo mass functions to explore the fraction and fate of stripped stars from satellites in galaxy groups and clusters of different masses. The majority of the stripped stars in massive halos are predicted to end up as intra-cluster stars, and the predicted amounts of the intra-cluster component as a function of the velocity dispersion of galaxy system match well the observational results obtained by Gonzalez et al. (2007). The fraction of the mass in the stripped stars to that remain bound in the central and satellite galaxies is the highest ($\sim 40%$ of the total stellar mass) in halos with masses $M_h\sim 10^{14}\msunh$. If all these stars end up in the intra-cluster component (Max), or maximum of them are accreted into the central galaxy (Min), then we can predict that a maximum $\sim 19%$ and a minimum $\sim 5%$ of the total stars in the whole universe are in terms of the diffused intra-cluster component. In the former case, in massive halos with $M_h \sim 10^{15} \msunh$, the stellar mass of the intra-cluster component is roughly 6 times as large as that of the central galaxy. This factor decreases to $\sim 2$, 1 and 0.1 in halos with $M_h \sim 10^{14}$, $10^{13}$, and $10^{12} \msunh$, respectively. The total amount of stars stripped from satellite galaxies is insufficient to build up the central galaxies in halos with masses $\la 10^{12.5}\msunh$, and so the quenching of star formation must occur in halos with higher masses. Abridged.

preprint2007arXiv

The cross-correlation between galaxies of different luminosities and Colors

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<M_r < -18.0),but for red galaxies the dependence is weaker and becomes insignificant for luminosities below L^*. To examine whether a significant stochastic/nonlinear component exists in the bias relation, we study the ratio R_ij= W_{ii}W_{jj}/W_{ij}^2, where W_{ij} is the projected correlation between subsample i and j. We find that the values of R_ij are all consistent with 1 for all-all, red-red and blue-blue samples, however significantly larger than 1 for red-blue samples. For faint red - faint blue samples the values of R_{ij} are as high as ~ 2 on small scales r_p < 1 \mpch and decrease with increasing r_p.

preprint2006arXiv

Inclination-Dependent Luminosity Function of Spiral Galaxies in the Sloan Digital Sky Survey: Implication for Dust Extinction

Using a samples of 61506 spiral galaxies selected from the SDSS DR2, we examine the luminosity function (LF) of spiral galaxies with different inclination angles. We find that the characteristic luminosity of the LF, $L^*$, decreases with increasing inclination, while the faint-end slope, $α$, depends only weakly on it. The inclination-dependence of the LF is consistent with that expected from a simple model where the optical depth is proportional to the cosine of the inclination angle, and we use a likelihood method to recover both the coefficient in front of the cosine, $γ$, and the LF for galaxies viewed face-on. The value of $γ$ is quite independent of galaxy luminosity in a given band, and the values of $γ$ obtained in this way for the 5 SDSS bands give an extinction curve which is a power law of wavelength ($τ\proptoλ^{-n}$), with a power index $n=0.96\pm0.04$. Using the dust extinction for galaxies obtained by Kauffmann et al. (2003), we derive an `extinction-corrected' luminosity function for spiral galaxies. Dust extinction makes $M^*$ dimmer by about 0.5 magnitudes in the $z$-band, and about 1.2 magnitudes in the $u$- band. Since our analysis is based on a sample where selection effects are well under control, the dimming of edge-on galaxies relative to face-on galaxies is best explained by assuming that galaxy disks are optically thick in dust absorptions.

preprint2004arXiv

Galaxy formation in pre-processed dark halos

Recent N-body simulations show that the formation of a present-day, galaxy sized dark matter halo in the CDM cosmogony in general consists of an early fast collapse phase, during which the potential associated with a halo is established, followed by a slow accretion phase, during which mass is added rather gently in the outer region. Here we outline a scenario in which the fast collapse phase is accompanied with rapid formation of cold clouds and with starbursts that can eject a large amount of gas from the halo center.The loss of orbital energy of the cold clouds to the dark matter and the ejection of gas from halo center by starburst can significantly reduce the halo concentration. The outflow from the starburst can also heat the gas in the protogalaxy region. Subsequent formation of galaxies in the slow accretion regime is therefore in halos that have been pre-processed by these processes and may have properties different from that given by N-body simulations. This scenario can help to solve several outstanding problems in the standard Lambda CDM model of galaxy formation without compromising its success in allowing structure formation at high redshift. The predicted rotation curves are quite flat,alleviating the discrepancy of the Tully-Fisher relation predicted in the Lambda CDM model with observations. The flattened galaxy halos allow accreted minihalos to survive in their central regions longer, which may be helpful in producing the flux anomalies observed in some gravitational lensing systems. The preheating by the early starbursts effectively reduces the amount of gas that can be accreted into galaxy halos, which may explain why the baryon fraction in a spiral galaxy is in general much lower than the universal baryon fraction, f_B~0.16, in the standard Lambda CDM model.

preprint2003arXiv

An analytic model for the galactic winds and mass outflows

Galactic winds and mass outflows are observed both in nearby starburst galaxies and in high-redshift star-forming galaxies. In this paper we develop a simple analytic model to understand the observed superwind phenomenon. Our model is built upon the model of McKee & Ostriker (1977) for the interstellar medium. It allows one to predict how properties of a superwind, such as wind velocity and mass outflow rate, are related to properties of its star-forming host galaxy, such as size, gas density and star formation rate. The model predicts a threshold of star formation rate density for the generation of observable galactic winds. Galaxies with more concentrated star formation produce superwinds with higher velocities. The predicted mass outflow rates are comparable to (or slightly larger than) the corresponding star formation rates. We apply our model to both local starburst galaxies and high-redshift Lyman break galaxies, and find its predictions to be in good agreement with current observations. Our model is simple, and so can be easily incorporated into numerical simulations and semi-analytical models of galaxy formation.

preprint2002arXiv

Galaxy Formation in Preheated Intergalactic Media

We outline a scenario of galaxy formation in which the gas in galaxy-forming regions was preheated to high entropy by vigorous energy feedback associated with the formation of stars in old ellipticals and bulges and with AGN activity. Such preheating likely occurred at redshifts z ~ 2-3, and can produce the entropy excess observed today in low-mass clusters of galaxies without destroying the bulk of the Lyman alpha forest. Subsequent galaxy formation is affected by the preheating, because the gas no longer follows the dark matter on galaxy scales. The hot gas around galaxy haloes has very shallow profiles and emits only weakly in the X-ray. Cooling in a preheated halo is not inside-out, because the cooling efficiency does not change significantly with radius. Only part of the gas in a protogalaxy region can cool and be accreted into the final galaxy halo. The accreted gas is likely in diffuse clouds and so does not lose angular momentum to the dark matter. Cluster ellipticals are produced by mergers of stellar systems formed prior to the preheating, while large galaxy disks form in low-density environments where gas accretion can continue to the present time.

preprint2002arXiv

Modelling the IGM and the Lyalpha forest at high redshift from the dark matter distribution

A variety of approximate schemes for modelling the low-density Intergalactic Medium (IGM) in the high-redshift Universe is compared to the results of a large high-resolution hydro-dynamical simulation. These schemes use either an analytical description of the dark matter distribution and the IGM or numerical simulations of the DM distributions combined with different approximate relations between dark matter field and the gas distribution. Schemes based on a filtering of the dark matter distribution with a global Jeans scale result in a rather poor description of the gas distribution. An adaptive filtering which takes into account the density/temperature dependence of the Jeans scale is required. A reasonable description of the gas distribution can be achieved using a fit of the mean relation between the dark matter and gas densities in the hydro-dynamical simulation to relate dark matter and gas distribution. In the hydro-dynamical simulations deviations from this mean relation are correlated with gradients in the dark matter peculiar velocity field indicative of shocks in the gas component. A scheme which takes into account this correlation results in a further improved gas distribution. Such adaptive filtering schemes applied to dark matter simulations will be very well suited for studies of statistical properties of the Lyalpha forest which investigate the IGM and the underlying dark matter distribution and require a large dynamic range and/or an extensive parameter study.

preprint2001arXiv

Observational signatures of feedback in QSO absorption spectra

Models for the formation of galaxies and clusters of galaxies require strong feedback in order to explain the observed properties of these systems. We investigate whether such feedback has observational consequences for the intergalactic medium, as probed in absorption towards background quasars. A typical quasar sight-line intersects one proto-cluster per unit redshift, and significant feedback from forming galaxies or AGN, heating the proto-cluster gas, will result in a large clearing of reduced absorption in the Lyman alpha-forest. Such a gap could be detected at redshift > 3 when the mean opacity is high. Feedback from Lyman-break galaxies in proto-clusters can be probed by the absorption lines produced in their winds. Strong feedback from galaxies has a major impact on the number and properties of absorption lines with neutral hydrogen column densities of 10^16/ cm^2. This feedback can be probed with multiple sight-lines and by studying the unsaturated higher-order lines of the Lyman series. Galactic winds from dwarf galaxies should break-up into clouds, in order not to over produce the number of absorption lines. These clouds can then coast to large distances.

preprint2000arXiv

Tully-Fisher Relation and its Implications for Halo Density Profile and Self-interacting Dark Matter

We show that the Tully-Fisher relation observed for spiral galaxies can be explained in the current scenario of galaxy formation without invoking subtle assumptions, provided that galactic-sized dark haloes have shallow, core-like central profiles, with a core radius proportional to halo circular velocity. In such a system, both disk and halo contribute significantly to the maximum rotation of the disk, and the interaction between the disk and halo components acts to reduce the scatter in the Tully-Fisher relation. With model parameters chosen in plausible ranges, the model can well accommodate the zero-point, slope, scatter of the observed Tully-Fisher relation, as well as the large range of disk sizes. The model predicts that LSB disks obey a Tully-Fisher relation similar to that of normal disks, if disk mass-to-light ratio is properly taken into account. The halo profile required by the Tully-Fisher relation is as shallow as that required by the rotation curves of faint disks, but much shallower than that predicted by conventional CDM models. Our results cannot be explained by some of the recent proposals for resolving the conflict between conventional CDM models and the rotation-curve shapes of faint galaxies. If dark matter self-interaction (either scattering or annihilation) is responsible for the shallow profile, the observed Tully-Fisher relation requires the interaction cross section σ_X to satisfy <σ_{X}|v|>/m_{X}~10^{-16} cm^3/s/GeV.

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

The Evolution of Correlation Functions and Power Spectra in Gravitational Clustering

Hamilton et al. (1991) proposed a simple formula relating the nonlinear autocorrelation function of the mass distribution to the primordial spectrum of density fluctuations for gravitational clustering in an $Ω=1$ universe. High resolution N-body simulations show this formula to work well for scale-free spectra $P(k) \propto k^n$ when the spectral index $n\sim 0$, but not when $n \lsim -1$. We show that a modified version of the formula can work well provided its form depends on $n$. This dependence can be derived from a simple physical model for collapse from Gaussian initial conditions. Our modified formula is easy to apply and is an excellent fit to N-body simulations with $0 \leq n \leq -2$. It can also be applied to non-power law initial spectra such as that of the standard Cold Dark Matter model by using the local spectral index at the current nonlinear scale as the effective value of $n$ at any given redshift. We give analytic expressions both for the nonlinear correlation function and for the nonlinear power spectrum.

preprint1995arXiv

An analytic model for the spatial clustering of dark matter haloes

We develop a simple analytic model for the gravitational clustering of dark matter haloes to understand how their spatial distribution is biased relative to that of the mass. The statistical distribution of dark haloes within the initial density field (assumed Gaussian) is determined by an extension of the Press-Schechter formalism. Modifications of this distribution caused by gravitationally induced motions are treated using a spherical collapse approximation. We test this model against results from a variety of N-body simulations, and find that it gives an accurate description of a bias function. This bias function is sufficient to calculate the cross-correlation between dark haloes and mass, and again we find excellent agreement between simulation results and analytic predictions. Because haloes are spatially exclusive, the variance in the count of objects within spheres of fixed radius and overdensity is significantly smaller than the Poisson value. This seriously complicates any analytic calculation of the autocorrelation function of dark halos. Our simulation results show that this autocorrelation function is proportional to that of the mass over a wide range in $R$, even including scales where both functions are significantly greater than unity. The constant of proportionality is very close to that predicted on large scales by the analytic model. This result permits an entirely analytic estimate of the autocorrelation function of dark haloes. We use our model to study how the distribution of galaxies may be biased with respect to that of the mass. In conjunction with other data these techniques should make it possible to measure the amplitude of cosmic mass fluctuations and the density of the Universe.

preprint1994arXiv

Lyman alpha absorption at low redshifts and hot gas in galactic haloes

Motivated by recent observation of Lanzetta et al. that most luminous galaxies at low redshifts produce $\lya$ absorption at impact parameter $l\lsim 160 \kpch$, we propose that these absorbers are clouds confined by the pressure of ambient hot gas in galactic haloes. We determine the properties of this hot gas and of the absorption systems on the basis of observational and theoretical constraints. The absorbing clouds need to be replenished on about one orbital time ($\sim 10^9$ yrs) in the galactic halo. The pressure and temperature of the gas at radius $r\sim 100\kpc$ are $P=(10-100){\rm cm^{-3} K}$, $T=10^{(5.5-6.5)}{\rm K}$. The model requires that about 10 per cent of the gas in low-redshift galactic haloes is in the hot phase. Such gas in galactic haloes emits x-ray with bolometric luminosity of the order $10^{37-40}{\rm erg\,s^{-1}}$. The plausibility for such gas to exist in current models of galaxy formation is discussed.

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

Spatial Distribution of Low Surface Brightness Galaxies

Using redshift samples, we calculate the cross-correlation functions of LSB galaxies with normal galaxies in complete samples (i.e. CfA and IRAS), which enables us to compare directly the amplitudes and shapes of the correlation functions. For pair separations $r\gs 2\mpc$, we find that the shape is in agreement with that of the correlation functions for other galaxies. The amplitudes ($A$) of $\xiab (r)$ are lower than those of the autocorrelation functions for the CfA and IRAS samples, with $A_{\rm LSB-CfA}:A_{\rm CfA-CfA}\approx 0.4$ and $A_{\rm LSB-IRAS}:A_{\rm IRAS-IRAS}\approx 0.6$. These results suggest that LSB galaxies are imbedded in the same large scale structure as other galaxies, but are less strongly clustered. This offers the hope that LSB galaxies may be unbiased tracers of the mass density on large scales. For $r\ls 2\mpc$, the cross-correlation functions are significantly lower than that expected from the extrapolation of $\xiab$ on larger scales, showing that the formation and survival of LSB galaxies may be inhibited by interaction with neighboring galaxies. The model which matches these observations suggests that strong luminosity segregation in galaxy clustering is not a necessary consequence of biased galaxy formation, unless the effect of surface brightness is taken in to account. It is also implies a significant mass density in LSB galaxies.

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.