Source author record

Xi Kang

Xi Kang 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

44works
4topics
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

44 published item(s)

preprint2026arXiv

Non-Parametric Attenuation Curves in Local Star-Forming Galaxies: Geometry Effect, Dust Evolution, and ISS

We introduce a non-parametric approach, the Stellar Population Synthesis with Equivalent Widths (SEW) method, to reconstruct spectrally-resolved attenuation curves for 169,568 star-forming galaxies from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7). Composite attenuation curves, stacked by stellar mass and inclination, reveal systematic trends: a higher stellar mass correlates with steeper slopes (lower $R_V$), while edge-on galaxies exhibit flatter curves due to geometric saturation effects. This flattening occurs because, as optical depth increases along the line of sight, the observed light becomes increasingly dominated by emission from the outer, less obscured layers of the galaxy. Using a simplified radiative transfer treatment based on a uniform dust-star mixture, we find the inclination-dependent slope variations are consistent with geometric effects, whereas the mass-dependent slope steepening indicates evolution in intrinsic dust properties, suggesting feedback-driven grain fragmentation in massive galaxies. Additionally, intermediate-scale structures (ISSs) are tentatively identified in the attenuation curves at approximately 4870, 6370, and 7690 Å. These results illustrate how the interplay among dust-star geometry, grain size evolution, and the galactic environment shapes attenuation curves.

preprint2025arXiv

Introduction to the Chinese Space Station Survey Telescope (CSST)

The Chinese Space Station Survey Telescope (CSST) is an upcoming Stage-IV sky survey telescope, distinguished by its large field of view (FoV), high image quality, and multi-band observation capabilities. It can simultaneously conduct precise measurements of the Universe by performing multi-color photometric imaging and slitless spectroscopic surveys. The CSST is equipped with five scientific instruments, i.e. Multi-band Imaging and Slitless Spectroscopy Survey Camera (SC), Multi-Channel Imager (MCI), Integral Field Spectrograph (IFS), Cool Planet Imaging Coronagraph (CPI-C), and THz Spectrometer (TS). Using these instruments, CSST is expected to make significant contributions and discoveries across various astronomical fields, including cosmology, galaxies and active galactic nuclei (AGN), the Milky Way and nearby galaxies, stars, exoplanets, Solar System objects, astrometry, and transients and variable sources. This review aims to provide a comprehensive overview of the CSST instruments, observational capabilities, data products, and scientific potential.

preprint2022arXiv

A stochastic model to reproduce the star-formation history of individual galaxies in hydrodynamic simulations

The star formation history (SFH) of galaxies is critical for understanding galaxy evolution. Hydrodynamical simulations enable us to precisely reconstruct the SFH of galaxies and establish a link to the underlying physical processes. In this work, we present a model to describe individual galaxies' SFHs from three simulations: TheThreeHundred, Illustris-1 and TNG100-1. This model divides the galaxy SFH into two distinct components: the "main sequence" and the "variation". The "main sequence" part is generated by tracing the history of the $SFR-M_*$ main sequence of galaxies across time. The "variation" part consists of the scatter around the main sequence, which is reproduced by fractional Brownian motions. We find that: 1) The evolution of the main sequence varies between simulations; 2) fractional Brownian motions can reproduce many features of SFHs, however, discrepancies still exist; 3) The variations and mass-loss rate are crucial for reconstructing the SFHs of the simulations. This model provides a fair description of the SFHs in simulations. On the other hand, by correlating the fractional Brownian motion model to simulation data, we provide a 'standard' against which to compare simulations.

preprint2022arXiv

About One-point Statistics of the Ratio of Two Fourier-transformed Cosmic Fields and an Application

The Fourier transformation is an effective and efficient operation of Gaussianization at the one-point level. Using a set of N-body simulation data, we verified that the one-point distribution functions of the dark matter momentum divergence and density fields closely follow complex Gaussian distributions. The one-point distribution function of the quotient of two complex Gaussian variables is introduced and studied. Statistical theories are then applied to model one-point statistics about the growth of individual Fourier mode of the dark matter density field, which can be obtained by the ratio of two Fourier transformed cosmic fields. Our simulation results proved that the models based on the Gaussian approximation are impressively accurate, and our analysis revealed many interesting aspects about the growth of dark matter's density fluctuation in Fourier space.

preprint2022arXiv

Cold-mode and hot-mode accretion in galaxy formation: an entropy approach

We have analysed two cosmological zoom simulations with $M_{\rm vir}\sim 10^{12}{\rm\,M}_\odot$ from the NIHAO series, both with and without feedback. We show that an entropy criterion based on the equation of state of the intergalactic medium can successfully separate cold- and hot-mode accretion. The shock-heated gas has non-negligible turbulent support and cools inefficiently. In the simulations without feedback, only a small fraction ($\sim 20$ per cent) of the stellar mass comes from baryons that have been in the hot circumgalactic medium, although quantitative conclusions should be taken with caution due to our small-number statistics. With feedback, the fraction is larger because of the reaccretion of gas heated by supernovae, which has lower entropies and shorter cooling times than the gas heated by accretion shocks. We have compared the results of NIHAO to predictions of the GalICS 2.1 semianalytic model of galaxy formation. The shock-stability criterion implemented in GalICS 2.1 successfully reproduces the transition from cold- to hot-mode accretion.

preprint2022arXiv

NIHAO XXVII: Crossing the green valley

The transition of high-mass galaxies from being blue and star forming to being red and dead is a crucial step in galaxy evolution, yet not fully understood. In this work, we use the NIHAO suite of galaxy simulations to investigate the relation between the transition time through the green valley and other galaxy properties. The typical green valley crossing time of our galaxies is approximately 400 Myr, somewhat shorter than observational estimates. The crossing of the green valley is triggered by the onset of AGN feedback and the subsequent shut down of star formation. Interestingly the time spent in the green valley is not related to any other galaxy properties, such as stellar age or metallicity, or the time at which the star formation quenching takes place. The crossing time is set by two main contributions: the ageing of the current stellar population and the residual star formation in the green valley. These effects are of comparable magnitude, while major and minor mergers have a negligible contribution. Most interestingly, we find the time that a galaxy spends to travel through the green valley is twice the $e$-folding time of the star formation quenching. This result is stable against galaxy properties and the exact numerical implementation of AGN feedback in the simulation. Assuming a typical crossing time of about one Gyr inferred from observations, our results imply that any mechanism or process aiming to quench star formation, must do it on a typical timescale of 500 Myr.

preprint2022arXiv

NIHAO XXVIII: Collateral effects of AGN on dark matter concentration and stellar kinematics

Although active galactic nuclei (AGN) feedback is required in simulations of galaxies to regulate star formation, further downstream effects on the dark matter distribution of the halo and stellar kinematics of the central galaxy can be expected. We combine simulations of galaxies with and without AGN physics from the Numerical Investigation of a Hundred Astrophysical Objects (NIHAO) to investigate the effect of AGN on the dark matter profile and central stellar rotation of the host galaxies. Specifically, we study how the concentration-halo mass ($c-M$) relation and the stellar spin parameter ($λ_R$) are affected by AGN feedback. We find that AGN physics is crucial to reduce the central density of simulated massive ($\gtrsim 10^{12}$ M$_\odot$) galaxies and bring their concentration to agreement with results from the Spitzer Photometry & Accurate Rotation Curves (SPARC) sample. Similarly, AGN feedback has a key role in reproducing the dichotomy between slow and fast rotators as observed by the ATLAS$^{3\text{D}}$ survey. Without star formation suppression due to AGN feedback, the number of fast rotators strongly exceeds the observational constraints. Our study shows that there are several collateral effects that support the importance of AGN feedback in galaxy formation, and these effects can be used to constrain its implementation in numerical simulations.

preprint2022arXiv

Quenching of Massive Disk Galaxies in the IllustrisTNG Simulation

A rare population of massive disk galaxies have been found to invade the red sequence dominated by early-type galaxies. These red/quenched massive disk galaxies have recently gained great interest into their formation and origins. The usually proposed quenching mechanisms, such as bar quenching and environment quenching, seem not suitable for those bulge-less quenched disks in low-density environment. In this paper, we use the IllustrisTNG-300 simulation to investigate the formation of massive quenched central disk galaxies. It is found that these galaxies contain less gas and harbor giant supermassive black holes(SMBHs) (above $ 10^{8}M_{\odot}$) than their star forming counterparts. By tracing their formation history, we found that quenched disk galaxies formed early and preserved disk morphology for cosmological time scales. They have experienced less than one major merger on average and it is mainly mini-mergers (mass ratio $<$1/10) that contribute to the growth of their SMBHs. In the Illustris-TNG simulation the black hole feedback mode switches from thermal to kinetic feedback when the black hole mass is more massive than $\sim 10^{8}M_{\odot}$, which is more efficient to eject gas outside of the galaxy and to suppress further cooling of hot gaseous halo. We conclude that kinetic AGN feedback in massive red/quenched disk galaxy is the dominant quenching mechanism.

preprint2022arXiv

Spin conservation of cosmic filaments

Cosmic filaments are the largest collapsing structure in the Universe. Recently both observations and simulations inferred that cosmic filaments have coherent angular momenta (spins). Here we use filament finders to identify the filamentary structures in cosmological simulations and study their physical origins, which are well described by the primordial tidal torque of their Lagrangian counterpart regions -- protofilaments. This initial angular momenta statistically preserve their directions to low redshifts. We further show that a spin reconstruction method can predict the spins of filaments and potentially relate their spins to the initial conditions of the Universe. This correlation provides a new way of constraining and obtaining additional information of the initial perturbations of the Universe.

preprint2022arXiv

Starduster: A multi-wavelength SED model based on radiative transfer simulations and deep learning

We present Starduster, a supervised deep learning model that predicts the multi-wavelength SED from galaxy geometry parameters and star formation history by emulating dust radiative transfer simulations. The model is comprised of three specifically designed neural networks, which take into account the features of dust attenuation and emission. We utilise the Skirt radiative transfer simulation to produce data for the training data of neural networks. Each neural network can be trained using ~ 4000 - 5000 samples. Compared with the direct results of the Skirt simulation, our deep learning model produces ~ 0.005 mag and ~ 0.1 - 0.2 mag errors for dust attenuation and emission respectively. As an application, we fit our model to the observed SEDs of IC4225 and NGC5166. Our model can reproduce the observations, and provide reasonable measurements of the inclination angle and stellar mass. However, some predicted geometry parameters are different from an image-fitting study. Our analysis implies that including a constraint at (rest-frame) ~ 40 micron could alleviate the degeneracy in the parameter space for both IC4225 and NGC5166, leading to broadly consistent results with the image-fitting predictions. Our SED code is publicly available and can be applied to both SED-fitting and SED-modelling of galaxies from semi-analytic models.

preprint2020arXiv

A robust determination of halo environment in the cosmic field

A number of methods for studying the large-scale cosmic matter distribution exist in the literature. One particularly common method employed to define the cosmic web is to examine the density, velocity or potential field. Such methods are advantageous since a Hessian matrix can be constructed whose eigenvectors (and eigenvalues) indicate the principal directions (and strength) of local collapse or expansion. Technically this is achieved by diagonalizing the Hessian matrix using a fixed finite grid. The resultant large-scale structure quantification is thus inherently limited by the grid's finite resolution. Here, we overcome the obstacle of finite grid resolution by introducing a new method to determine halo environment using an adaptive interpolation which is more robust to resolution than the typical "Nearest Grid Point" (NGP) method. Essentially instead of computing and diagonalizing the Hessian matrix once for the entire grid, we suggest doing so once for each halo or galaxy in question. We examine how the eigenvalues and eigenvector direction's computed using our algorithm and the NGP method converge for different grid resolutions, finding that our new method is convergent faster. Namely changes of resolution have a much smaller effect than in the NGP method. We therefore suggest this method for future use by the community.

preprint2020arXiv

NIHAO-XXIII: Dark Matter density shaped by Black Hole feedback

We present a systematic analysis of the reaction of dark matter distribution to galaxy formation across more than eight orders of magnitude in stellar mass. We extend the previous work presented in the NIHAO-IV paper (Tollet et al.) by adding 46 new high resolution simulations of massive galaxies performed with the inclusion of Black Hole feedback. We show that outflows generated by the AGN are able to partially counteract the dark matter contraction due to the large central stellar component in massive haloes. The net effect is to relax the central dark matter distribution that moves to a less cuspy density profiles at halo masses larger than $\approx 3 \times 10^{12}$ M$_{\odot}$. The scatter around the mean value of the density profile slope ($α$) is fairly constant ($Δα\approx 0.3$), with the exception of galaxies with halo masses around $10^{12}$ M$_{\odot}$, at the transition from stellar to AGN feedback dominated systems, where the scatter increases by almost a factor three. We provide useful fitting formulas for the slope of the dark matter density profiles at few percent of the virial radius for the whole stellar mass range: $10^5-10^{12}$ M$_{\odot}$ ($2 \times 10^9-5 \times 10^{13}$ M$_{\odot}$ in halo mass).

preprint2020arXiv

The alignment of satellite systems with cosmic filaments in the SDSS DR12

Galaxies, as well as their satellites, are known to form within the cosmic web: the large, multi-scale distribution of matter in the universe. It is known that the surrounding large scale structure (LSS) can impact and influence the formation of galaxies, e.g. the spin and shape of haloes or galaxies are correlated with the LSS and the correlation depends on halo mass or galaxy morphology. In this work, we use group and filament catalogues constructed from the SDSS DR12 to investigate the correlation between satellite systems and the large scale filaments they are located in. We find that the distribution of satellites is significantly correlated with filaments, namely the major axis of the satellite systems are preferentially aligned with the spine of the closest filament. Stronger alignment signals are found for the cases where the system away from the filament spine, while systems close to the filament spine show significantly weaker alignment. Our results suggest that satellites are accreted along filaments, which agrees with previous works. The case of which away from the filament spine may help us to understand how the filament forms as well as the peculiar satellite distribution in the Local Universe.

preprint2020arXiv

The formation of blue cluster in local Universe

It is well known from the Butcher-Oemler effect that galaxies in dense environment are mostly red with little star formation and the fraction of blue galaxies in galaxy groups/clusters also declines rapidly with redshifts. A recent work by Hashimoto et al. reported a local 'blue cluster' with high fraction of blue galaxies ($\sim 0.57$), higher than the model predictions. They ascribed this blue cluster to the feeding of gas along a filamentary structure around the cluster. In this work we use group catalog from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7) and the state-of-art of semi-analytic model (SAM) to investigate the formation of blue clusters in local Universe. In total, we find four blue clusters with halo mass $\sim 10^{14}M_{\odot}$ at $0.02 < z < 0.082$, while only the one found by Hashimoto et al. is in a filamentary structure. The SAM predicts that blue clusters have later formation time and most blue satellite galaxies are recently accreted. We conclude that the formation of blue clusters is mainly governed by newly accreted blue satellites, rather than the effect of large-scale environment.

preprint2020arXiv

What has quenched the massive spiral galaxies?

Quenched massive spiral galaxies have attracted great attention recently, as more data is available to constrain their environment and cold gas content. However, the quenching mechanism is still uncertain, as it depends on the mass range and baryon budget of the galaxy. In this letter, we report the identification of a rare population of very massive, quenched spiral galaxies with stellar mass $\gtrsim10^{11}{\rm~M_\odot}$ and halo mass $\gtrsim10^{13}{\rm~M_\odot}$ from the Sloan Digital Sky Survey at redshift $z\sim0.1$. Our CO observations using the IRAM-30m telescope show that these galaxies contain only a small amount of molecular gas. Similar galaxies are also seen in the state-of-the-art semi-analytical models and hydro-dynamical simulations. It is found from these theoretical models that these quenched spiral galaxies harbor massive black holes, suggesting that feedback from the central black holes has quenched these spiral galaxies. This quenching mechanism seems to challenge the popular scenario of the co-evolution between massive black holes and massive bulges.

preprint2019arXiv

Measuring the HI mass function below the detection threshold

We present a Bayesian Stacking technique to directly measure the HI mass function (HIMF) and its evolution with redshift using galaxies formally below the nominal detection threshold. We generate galaxy samples over several sky areas given an assumed HIMF described by a Schechter function and simulate the HI emission lines with different levels of background noise to test the technique. We use Multinest to constrain the parameters of the HIMF in a broad redshift bin, demonstrating that the HIMF can be accurately reconstructed, using the simulated spectral cube far below the HI mass limit determined by the $5σ$ flux-density limit, i.e. down to $M_{\rm HI} = 10^{7.5}$ M$_{\odot}$ over the redshift range $0 < z < 0.55$ for this particular simulation, with a noise level similar to that expected for the MIGHTEE survey. We also find that the constraints on the parameters of the Schechter function, $ϕ_{\star}$, $M_\star$ and $α$ can be reliably fit, becoming tighter as the background noise decreases as expected, although the constraints on the redshift evolution are not significantly affected. All the parameters become better constrained as the survey area increases. In summary, we provide an optimal method for estimating the HI mass at cosmological distances that allows us to constrain the HI mass function below the detection threshold in forthcoming HI surveys. This study is a first step towards the measurement of the HIMF at high ($z>0.1$) redshifts.

preprint2019arXiv

Warm dark matter model with a few keV mass is bad for the too-big-to-fail problem

Theoretical studying of the very inner structure of faint satellite galaxy requires very high-resolution hydro-dynamical simulations with realistic models for star formation, which are beginning to emerge recently. In this work we present an analytical description to model the inner kinematic of satellites in the Milky Way (MW). We use a Monte-Carlo method to produce merger trees for MW mass halo and analytical models to produce stellar mass in the satellite galaxies. We consider two important processes which can significantly modify the inner mass distribution in satellite galaxy. The first is baryonic feedback which can induce a flat inner profile depending on the star formation efficiency in the galaxy. The second is the tidal stripping to reduce and re-distribute the mass inside satellite. We apply this model to MW satellite galaxies in both CDM and thermal relic WDM models. It is found that tidal heating must be effective to produce a relatively flat distribution of the satellite circular velocities, to agree with the data. The constraint on WDM mass depends on the host halo mass. For a MW halo with dark matter mass lower than $2\times 10^{12}M_{\odot}$, a 2 keV WDM model can be safely excluded as the predicted satellite circular velocities are systematically lower than the data. For WDM with mass of 3.5 keV, it requires the MW halo mass to be larger than $1.5\times 10^{12}M_{\odot}$, otherwise the 3.5 Kev model can also be excluded. Our current model can not exclude the WDM model with mass larger than 10 Kev.

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

NIHAO IX: the role of gas inflows and outflows in driving the contraction and expansion of cold dark matter haloes

We use ~100 cosmological galaxy formation zoom-in simulations using the smoothed particle hydrodynamics code {\sc gasoline} to study the effect of baryonic processes on the mass profiles of cold dark matter haloes. The haloes in our study range from dwarf (M_{200}~10^{10}Msun) to Milky Way (M_{200}~10^{12}Msun) masses. Our simulations exhibit a wide range of halo responses, primarily varying with mass, from expansion to contraction, with up to factor ~10 changes in the enclosed dark matter mass at one per cent of the virial radius. Confirming previous studies, the halo response is correlated with the integrated efficiency of star formation: e_SF=(M_{star}/M_{200})/(Ω_b/Ω_m). In addition we report a new correlation with the compactness of the stellar system: e_R=r_{1/2}/R_{200}. We provide an analytic formula depending on e_SF and e_R for the response of cold dark matter haloes to baryonic processes. An observationally testable prediction is that, at fixed mass, larger galaxies experience more halo expansion, while the smaller galaxies more halo contraction. This diversity of dark halo response is captured by a toy model consisting of cycles of adiabatic inflow (causing contraction) and impulsive gas outflow (causing expansion). For net outflow, or equal inflow and outflow fractions, f, the overall effect is expansion, with more expansion with larger f. For net inflow, contraction occurs for small f (large radii), while expansion occurs for large f (small radii), recovering the phenomenology seen in our simulations. These regularities in the galaxy formation process provide a step towards a fully predictive model for the structure of cold dark matter haloes.

preprint2016arXiv

NIHAO project II: Halo shape, phase-space density and velocity distribution of dark matter in galaxy formation simulations

We use the NIHAO (Numerical Investigation of Hundred Astrophysical Objects) cosmological simulations to study the effects of galaxy formation on key properties of dark matter (DM) haloes. NIHAO consists of $\simeq 90$ high-resolution SPH simulations that include (metal-line) cooling, star formation, and feedback from massive stars and SuperNovae, and cover a wide stellar and halo mass range: $10^6 < M_* / M_{\odot} < 10^{11}$ ( $10^{9.5} < M_{\rm halo} / M_{\odot} < 10^{12.5}$). When compared to DM-only simulations, the NIHAO haloes have similar shapes at the virial radius, R_{\rm vir}, but are substantially rounder inside $\simeq 0.1R_{\rm vir}$. In NIHAO simulations $c/a$ increases with halo mass and integrated star formation efficiency, reaching $\sim 0.8$ at the Milky Way mass (compared to 0.5 in DM-only), providing a plausible solution to the long-standing conflict between observations and DM-only simulations. The radial profile of the phase-space $Q$ parameter ($ρ/σ^3$) is best fit with a single power law in DM-only simulations, but shows a flattening within $\simeq 0.1R_{\rm vir}$ for NIHAO for total masses $M>10^{11} M_{\odot}$. Finally, the global velocity distribution of DM is similar in both DM-only and NIHAO simulations, but in the solar neighborhood, NIHAO galaxies deviate substantially from Maxwellian. The distribution is more symmetric, roughly Gaussian, with a peak that shifts to higher velocities for Milky Way mass haloes. We provide the distribution parameters which can be used for predictions for direct DM detection experiments. Our results underline the ability of the galaxy formation processes to modify the properties of dark matter haloes.

preprint2016arXiv

NIHAO X: Reconciling the local galaxy velocity function with Cold Dark Matter via mock HI observations

We used 87 high resolution hydrodynamical cosmological simulations from the NIHAO suite to investigate the relation between the maximum circular velocity (Vmax) of a dark matter halo in a collisionless simulation and the velocity width of the HI gas in the same halo in the hydrodynamical simulation. These two quantities are normally used to compare theoretical and observational velocity functions and have led to a possible discrepancy between observations and predictions based on the Cold Dark Matter (CDM) model. We show that below 100 km/s, there is clear bias between HI based velocities and Vmax, that leads to an underestimation of the actual circular velocity of the halo. When this bias is taken into account the CDM model has no trouble in reproducing the observed velocity function and no lack of low velocity galaxies is actually present. Our simulations also reproduce the linewidth - stellar mass (Tully-Fisher) relation and HI sizes, indicating that the HI gas in our simulations is as extended as observed. The physical reason for the lower than expected linewidths is that, in contrast to high mass galaxies, low mass galaxies no longer have extended thin HI rotating disks, as is commonly assumed.

preprint2016arXiv

Resolution-independent modeling of environmental effects in semi-analytic models of galaxy formation that include ram-pressure stripping of both hot and cold gas

The quenching of star formation in satellite galaxies is observed over a wide range of dark matter halo masses and galaxy environments. In the recent Guo et al (2011) and Fu et al (2013) semi-analytic + N-body models, the gaseous environment of the satellite galaxy is governed by the properties of the dark matter subhalo in which it resides. This quantity depends of the resolution of the N-body simulation, leading to a divergent fraction of quenched satellites in high- and low-resolution simulations. Here, we incorporate an analytic model to trace the subhaloes below the resolution limit. We demonstrate that we then obtain better converged results between the Millennium I and II simulations, especially for the satellites in the massive haloes ($\rm log M_{halo}=[14,15]$). We also include a new physical model for the ram-pressure stripping of cold gas in satellite galaxies. However, we find very clear discrepancies with observed trends in quenched satellite galaxy fractions as a function of stellar mass at fixed halo mass. At fixed halo mass, the quenched fraction of satellites does not depend on stellar mass in the models, but increases strongly with mass in the data. In addition to the over-prediction of low-mass passive satellites, the models also predict too few quenched central galaxies with low stellar masses, so the problems in reproducing quenched fractions are not purely of environmental origin. Further improvements to the treatment of the gas-physical processes regulating the star formation histories of galaxies are clearly necessary to resolve these problems.

preprint2016arXiv

The gap of stellar mass in galaxy groups: another perspective of the Too-big-To-Fail problem in the Milky Way

The Milky Way presents the too-big-to-fail (TBTF) problem that there are two observed satellite galaxies with maximum circular velocity larger than 55km/s, and others have velocity less than 25km/s, but the cold dark matter model predicts there should be more than 10 subhaloes with velocity larger than 25km/s. Those massive subhaloes with $25km/s<V_{max}<55km/s$ should not have failed to form stars. The TBTF problem severely challenges the CDM model. Most efforts are seeking the effects of baryonic feedback, decreasing the mass of the Milky Way, changing the properties of dark matter, so as to assign the observed low-velocity satellites into the massive subhaloes found in simulations. However, the TBTF problem can be avoided if the MW have not accreted subhaloes with velocity between $25km/s<V_{max}<55km/s$ although the probability of such a gap is lower as $\sim 1\%$ and can not be tested against observations. In this work we study the gap in stellar mass of satellite galaxies using the SDSS group catalogue and a semi-analytical model. We find that there are 1-2\% of galaxy groups with a large gap in the stellar mass of their satellites. These 'big gap' groups have accreted less massive subhaloes in their formation history and naturally display a gap between their satellite galaxies. If extrapolating our results to the Milky Way is appropriate, we conclude that it is very likely that our Milky Way has not accreted enough massive subhaloes to host those low-velocity satellites, and the TBTF problem is naturally avoided.

preprint2015arXiv

A study on the multicolour evolution of Red Sequence galaxy populations: insights from hydrodynamical simulations and semi-analytical models

By means of our own cosmological-hydrodynamical simulation and semi-analytical model we studied galaxy population properties in clusters and groups, spanning over 10 different bands from UV to NIR, and their evolution since redshift z=2. We compare our results in terms of galaxy red/blue fractions and luminous-to-faint ratio (LFR) on the Red Sequence (RS) with recent observational data reaching beyond z=1.5. Different selection criteria were tested in order to retrieve galaxies belonging to the RS: either by their quiescence degree measured from their specific SFR ("Dead Sequence"), or by their position in a colour-colour plane which is also a function of sSFR. In both cases, the colour cut and the limiting magnitude threshold were let evolving with redshift, in order to follow the natural shift of the characteristic luminosity in the LF. We find that the Butcher-Oemler effect is wavelength-dependent, with the fraction of blue galaxies increasing steeper in optical colours than in NIR. Besides, only when applying a lower limit in terms of fixed absolute magnitude, a steep BO effect can be reproduced, while the blue fraction results less evolving when selecting samples by stellar mass or an evolving magnitude limit. We then find that also the RS-LFR behaviour, highly debated in the literature, is strongly dependent on the galaxy selection function: in particular its very mild evolution recovered when measured in terms of stellar mass, is in agreement with values reported for some of the highest redshift confirmed (proto)clusters. As to differences through environments, we find that normal groups and (to a lesser extent) cluster outskirts present the highest values of both star forming fraction and LFR at low z, while fossil groups and cluster cores the lowest: this separation among groups begins after z~0.5, while earlier all group star forming properties are undistinguishable.

preprint2015arXiv

NIHAO project I: Reproducing the inefficiency of galaxy formation across cosmic time with a large sample of cosmological hydrodynamical simulations

We introduce project NIHAO (Numerical Investigation of a Hundred Astrophysical Objects), a set of 100 cosmological zoom-in hydrodynamical simulations performed using the gasoline code, with an improved implementation of the SPH algorithm. The haloes in our study range from dwarf to Milky Way masses, and represent an unbiased sampling of merger histories, concentrations and spin parameters. The particle masses and force softenings are chosen to resolve the mass profile to below 1% of the virial radius at all masses, ensuring that galaxy half-light radii are well resolved. Using the same treatment of star formation and stellar feedback for every object, the simulated galaxies reproduce the observed inefficiency of galaxy formation across cosmic time as expressed through the stellar mass vs halo mass relation, and the star formation rate vs stellar mass relation. We thus conclude that stellar feedback is the chief piece of physics required to limit the efficiency of star formation in galaxies less massive than the Milky Way.

preprint2015arXiv

NIHAO V: Too big doesn't fail -- reconciling the conflict between LCDM predictions and the circular velocities of nearby field galaxies

We compare the half-light circular velocities, V_{1/2}, of dwarf galaxies in the Local Group to the predicted circular velocity curves of galaxies in the NIHAO suite of LCDM simulations. We use a subset of 34 simulations in which the central galaxy has a stellar luminosity in the range 0.5 x 10^5 < L_V < 2 x 10^8 L_{sun}. The NIHAO galaxy simulations reproduce the relation between stellar mass and halo mass from abundance matching, as well as the observed half-light size vs luminosity relation. The corresponding dissipationless simulations over-predict the V_{1/2}, recovering the problem known as too big to fail (TBTF). By contrast, the NIHAO simulations have expanded dark matter haloes, and provide an excellent match to the distribution of V_{1/2} for galaxies with L_V > 2 x 10^6 L_{sun}. For lower luminosities our simulations predict very little halo response, and tend to over predict the observed circular velocities. In the context of LCDM, this could signal the increased stochasticity of star formation in haloes below M_{halo} \sim 10^{10} M_{sun}, or the role of environmental effects. Thus, haloes that are "too big to fail", do not fail LCDM, but haloes that are "too small to pass" (the galaxy formation threshold) provide a future test of LCDM.

preprint2015arXiv

Semi-Analytic Model Predictions of Mass Segregation from Groups to Clusters

Taking advantage of a high-resolution simulation coupled with a state-of-art semi-analytic model of galaxy formation, we probe the mass segregation of galaxies in groups and clusters, focusing on which physical mechanisms are driving it. We find evidence of mass segregation in groups and clusters up to the virial radius, both looking at the galaxy stellar mass and subhalo mass. The physical mechanism responsible for that is consistent with dynamical friction, a drag-force that brings more massive galaxies faster towards the innermost regions of the halo. At odds with observational results, we do not find the inclusion of low-mass galaxies in the samples, down to stellar mass $M_* = 10^9 \, M_{\odot}$, to change the overall trend shown by intermediate and massive galaxies. Moreover, stellar stripping as well as the growth of galaxies after their accretion, do not contribute either in shaping mass segregation or mixing the radial mass distribution. Beyond the virial radius we find an "anti-mass segregation" in groups that progressively weakens in clusters. The continuous accretion of new objects and recent merger events play a different role depending on the halo mass onto which accreting material is falling.

preprint2015arXiv

The accretion of dark matter subhaloes within the cosmic web: primordial anisotropic distribution and its universality

The distribution of galaxies displays anisotropy on different scales and it is often referred as galaxy alignment. To understand the origin of galaxy alignments on small scales, one must investigate how galaxies were accreted in the early universe and quantify their primordial anisotropic at the time of accretion. In this paper we use N-body simulations to investigate the accretion of dark matter subhaloes, focusing on their alignment with the host halo shape and the orientation of mass distribution on large scale, defined using the hessian matrix of the density field. The large/small (e1/e3) eigenvalues of the hessian matrix define the fast/slow collapse direction of dark matter on large scale. We find that: 1) the halo major axis is well aligned with the e3 (slow collapse) direction, and it is stronger for massive haloes; 2) subhaloes are predominately accreted along the major axis of the host halo, and the alignment increases with the host halo mass. Most importantly, this alignment is universal; 3) accretion of subhaloes with respect to the e3 direction is not universal. In massive haloes, subhaloes are accreted along the e3 (even stronger than the alignment with the halo major axis), but in low-mass haloes subhaloes are accreted perpendicular to the e3. The transit mass is lower at high redshift. The last result well explains the puzzled correlation (both in recent observations and simulations) that massive galaxies/haloes have their spin perpendicular to the filament, and the spin of low-mass galaxies/haloes is slightly aligned with the filament, under the assumption that the orbital angular momentum of subhaloes is converted to halo spin.

preprint2014arXiv

Cold stream stability during minor mergers

We use high-resolution Eulerian simulations to study the stability of cold gas flows in a galaxy size dark matter halo (10^12 Msun) at redshift z=2. Our simulations show that a cold stream penetrating a hot gaseous halo is stable against thermal convection and Kelvin-Helmholtz instability. We then investigate the effect of a satellite orbiting the main halo in the plane of the stream. The satellite is able to perturb the stream and to inhibit cold gas accretion towards the center of the halo for 0.5 Gyr. However, if the supply of cold gas at large distances is kept constant, the cold stream is able to re-establish itself after 0.3 Gyr. We conclude that cold streams are very stable against a large variety of internal and external perturbations.

preprint2014arXiv

Satellite Alignment: I. Distribution of Substructures and Their Dependence On Assembly History From N-Body Simulations

Observations have shown that the spatial distribution of satellite galaxies is not random, but aligned with the major axes of central galaxies. This alignment is dependent on galaxy properties, such that red satellites are more strongly aligned than blue satellites. Theoretical work done to interpret this phenomena has found that it is due to the non-spherical nature of dark matter halos. However, most studies over-predict the alignment signal under the assumption that the central galaxy shape follows the shape of the host halo. It is also not clear whether the color dependence of alignment is due to an assembly bias or an evolution effect. In this paper we study these problems using a cosmological N-body simulation. Subhalos are used to trace the positions of satellite galaxies. It is found that the shape of dark matter halos are mis-aligned at different radii. If the central galaxy shares the same shape as the inner host halo, then the alignment effect is weaker and agrees with observational data. However, it predicts almost no dependence of alignment on the color of satellite galaxies, though the late accreted subhalos show stronger alignment with the outer layer of the host halo than their early accreted counterparts. We find that this is due to the limitation of pure N-body simulations that satellites galaxies without associated subhalos ('orphan galaxies') are not resolved. These orphan (mostly red) satellites often reside in the inner region of host halos and should follow the shape of the host halo in the inner region.

preprint2014arXiv

The Distribution of Satellites Around Central Galaxies in a Cosmological Hydrodynamical Simulation

Observations have shown that the spatial distribution of satellite galaxies is not random, but rather is aligned with the major axes of central galaxies (CGs). The strength of the alignment is dependent on the properties of both the satellites and centrals. Theoretical studies using dissipationless N-body simulations are limited by their inability to directly predict the shape of CGs. Using hydrodynamical simulations including gas cooling, star formation, and feedback, we carry out a study of galaxy alignment and its dependence on the galaxy properties predicted directly from the simulations.We found that the observed alignment signal is well produced, as is the color dependence: red satellites and red centrals both show stronger alignments than their blue counterparts. The reason for the stronger alignment of red satellites is that most of them stay in the inner region of the dark matter halo where the shape of the CG better traces the dark matter distribution. The dependence of alignment on the color of CGs arises from the halo mass dependence, since the alignment between the shape of the central stellar component and the inner halo increases with halo mass. We also find that the alignment of satellites is most strongly dependent on their metallicity, suggesting that the metallicity of satellites, rather than color, is a better tracer of galaxy alignment on small scales. This could be tested in future observational studies.

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

Resolving the problem of galaxy clustering on small scales: any new physics needed?

Galaxy clustering sets strong constraints on the physics governing galaxy formation and evolution. However, most current models fail to reproduce the clustering of low-mass galaxies on small scales ($r<1Mpc/h$). In this paper we study the galaxy clusterings predicted from a few semi-analytical models. We firstly compare two Munich versions, Guo et al. (2011, Guo11) and De Lucia \& Blazoit (2007, DLB07). The Guo11 model well reproduces the galaxy stellar mass function, but over-predicts the clustering of low-mass galaxies on small scales. The DLB07 model provides a better fit to the clustering on small scales, but over-predicts the stellar mass function. These seem to be puzzling. We find that there is slightly more fraction of satellite galaxies residing in massive haloes in the Guo11 model, which is the dominant contribution to the clustering discrepancy between the two models. However, both models still over-predict the clustering at $0.1Mpc/h<r<10Mpc/h$ for low mass galaxies. This is because both models over-predict the number of satellites by $30\%$ in massive halos than the data. Actually, the better agreement of DLB07 model with the data on small scales comes as a coincidence as it predicts too many low-mass central galaxies which are less clustered and thus bring down the total clustering. Finally, we show the predictions from the semi-analytical of Kang et al. (2012). We find that this model can simultaneously fit the stellar mass function and galaxy clustering if the supernova feedback in satellite galaxies is stronger. We conclude that semi-analytical models are now able to solve the small-scales clustering problem, without invoking of any other new physics or changing the dark matter properties, such as the recent favored warm dark matter.

preprint2013arXiv

The dependence of tidal stripping efficiency on the satellite and host galaxy morphology

In this paper we study the tidal stripping process for satellite galaxies orbiting around a massive host galaxy, and focus on its dependence on the morphology of both satellite and host galaxy. For this purpose, we use three different morphologies for the satellites: pure disc, pure bulge and a mixture bulge+disc. Two morphologies are used for the host galaxies: bulge+disc and pure bulge. We find that while the spheroidal stellar component experiences a constant power-law like mass removal, the disc is exposed to an exponential mass loss when the tidal radius of the satellite is of the same order of the disc scale length. This dramatic mass loss is able to completely remove the stellar component on time scale of 100 Myears. As a consequence two satellites with the same stellar and dark matter masses, on the same orbit could either retain considerable fraction of their stellar mass after 10 Gyrs or being completely destroyed, depending on their initial stellar morphology. We find that there are two characteristic time scales describing the beginning and the end of the disc removal, whose values are related to the size of the disc. This result can be easily incorporated in semi-analytical models. We also find that the host morphology and the orbital parameters also have an effect on the determining the mass removal, but they are of secondary importance with respect to satellite morphology. We conclude that satellite morphology has a very strong effect on the efficiency of stellar stripping and should be taken into account in modeling galaxy formation and evolution.

preprint2013arXiv

The effect of Warm Dark Matter on galaxy properties: constraints from the stellar mass function and the Tully-Fisher relation

In this paper we combine high resolution N-body simulations with a semi analytical model of galaxy formation to study the effects of a possible Warm Dark Matter (WDM) component on the observable properties of galaxies. We compare three WDM models with a dark matter mass of 0.5, 0.75 and 2.0 keV, with the standard Cold Dark Matter case. For a fixed set of parameters describing the baryonic physics the WDM models predict less galaxies at low (stellar) masses, as expected due to the suppression of power on small scales, while no substantial difference is found at the high mass end. However these differences in the stellar mass function, vanish when different set of parameters are used to describe the (largely unknown) galaxy formation processes. We show that is possible to break this degeneracy between DM properties and the parameterization of baryonic physics by combining observations on the stellar mass function with the Tully-Fisher relation (the relation between stellar mass and the rotation velocity at large galactic radii as probed by resolved HI rotation curves). WDM models with a too warm candidate (m<0.75 keV) cannot simultaneously reproduce the stellar mass function and the Tully-Fisher relation. We conclude that accurate measurements of the galaxy stellar mass function and the link between galaxies and dark matter haloes down to the very low-mass end can give very tight constraints on the nature of DM candidates.

preprint2012arXiv

Bulk flow of halos in ΛCDM simulation

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 ΛCDM model we used; some results deviate from expectation at a level of ~ 3σ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.

preprint2011arXiv

Quantifying Kinematic Substructure in the Milky Way's Stellar Halo

We present and analyze the positions, distances, and radial velocities for over 4000 blue horizontal-branch (BHB) stars in the Milky Way's halo, drawn from SDSS DR8. We search for position-velocity substructure in these data, a signature of the hierarchical assembly of the stellar halo. Using a cumulative "close pair distribution" (CPD) as a statistic in the 4-dimensional space of sky position, distance, and velocity, we quantify the presence of position-velocity substructure at high statistical significance among the BHB stars: pairs of BHB stars that are close in position on the sky tend to have more similar distances and radial velocities compared to a random sampling of these overall distributions. We make analogous mock-observations of 11 numerical halo formation simulations, in which the stellar halo is entirely composed of disrupted satellite debris, and find a level of substructure comparable to that seen in the actually observed BHB star sample. This result quantitatively confirms the hierarchical build-up of the stellar halo through a signature in phase (position-velocity) space. In detail, the structure present in the BHB stars is somewhat less prominent than that seen in most simulated halos, quite possibly because BHB stars represent an older sub-population. BHB stars located beyond 20 kpc from the Galactic center exhibit stronger substructure than at $\rm r_{gc} < 20$ kpc.

preprint2011arXiv

What Sets the Sizes of the Faintest Galaxies?

We provide a comprehensive description and offer an explanation for the sizes of the faintest known galaxies in the universe, the dwarf spheroidal (dSph) satellites of the Milky Way and Andromeda. After compiling a consistent data set of half-light radii (r_{1/2}) and luminosities, we describe the size-luminosity relation of dSphs by a log-normal distribution in r_{1/2} with a mean size that varies as a function of luminosity. Accounting for modest number statistics, measurement uncertainties and surface brightness limitations, we find that the size-luminosity relations of the Milky Way and Andromeda dSph populations are statistically indistinguishable, and also very similar: their mean sizes at a given stellar luminosity differ by no more than 30%. In addition, we find that the mean size, slope and scatter of this log-normal size description of Local Group dSphs matches onto the relation of more massive low-concentration galaxies. This suggests that the stellar sizes of dSphs are ultimately related to their overall initial baryonic angular momentum. To test this hypothesis we perform a series of high resolution N-body simulations that we couple with a semi-analytic model of galaxy formation. These predict the same mean size and slope as observed in dSph satellites. At the same time, these models predict that the size-luminosity distributions for satellite galaxies around similar host-halos must be similar providing a natural explanation as to why the size distributions of Milky Way and Andromeda satellites are similar. Although strong rotation is currently not observed in dSphs, this may well be consistent with our angular-momentum-based explanation for their sizes if the disks of these galaxies have become sufficiently stirred through tidal interaction.

preprint2010arXiv

An Improved Model for the Dynamical Evolution of Dark Matter Subhaloes

Using an analytical model, we study the evolution of subhalo, including its mass, angular momentum and merging time-scale. This model considers the dominant processes governing subhalo evolution, such as dynamical friction, tidal stripping and tidal heating. We find that in order to best match the evolution of angular momentum measured from N-body simulation, mass stripping by tidal force should become inefficient after subhalo has experienced a few passages of pericenter. It is also found that the often used Coulomb logarithm $\ln M/m$ has to be revised to best fit the merging time-scales from simulation. Combining the analytical model with the Extended Press-Schechter (EPS) based merger trees, we study the subhalo mass function, and their spatial distribution in a Milky-Way (MW) type halo. By tuning the tidal stripping efficiency, we can gain a better match to the subhalo mass function from simulation. The predicted distribution of subhaloes is found to agree with the distribution of MW satellites, but is more concentrated than the simulation results. The radial distribution of subhaloes depends weakly on subhaloes mass at both present day and the time of accretion, but strongly on the accretion time. Using the improved model, we measure the second moment of the subhalo occupation distribution, and it agrees well with the results of Kravtsov et al. (2004a) and Zheng et al. (2005).

preprint2010arXiv

Galaxy formation with cold gas accretion and evolving stellar initial mass function

The evolution of the galaxy stellar mass function is especially useful to test the current model of galaxy formation. Observational data have revealed a few inconsistencies with predictions from the $Λ{\rm CDM}$ model. For example, most massive galaxies have already been observed at very high redshifts, and they have experienced only mild evolution since then. In conflict with this, semi-analytical models of galaxy formation predict an insufficient number of massive galaxies at high redshift and a rapid evolution between redshift 1 and 0 . In addition, there is a strong correlation between star formation rate and stellar mass for star-forming galaxies, which can be roughly reproduced with the model, but with a normalization that is too low at high redshift. Furthermore, the stellar mass density obtained from the integral of the cosmic star formation history is higher than the measured one by a factor of 2. In this paper, we study these issues using a semi-analytical model that includes: 1) cold gas accretion in massive halos at high redshift; 2) tidal stripping of stellar mass from satellite galaxies; and 3) an evolving stellar initial mass function (bottom-light) with a higher gas recycle fraction. Our results show that the combined effects from 1) and 2) can predict sufficiently massive galaxies at high redshifts and reproduce their mild evolution at low redshift, While the combined effects of 1) and 3) can reproduce the correlation between star formation rate and stellar mass for star-forming galaxies across wide range of redshifts. A bottom-light/top-heavy stellar IMF could partly resolve the conflict between the stellar mass density and cosmic star formation history.

preprint2010arXiv

The dependence of AGN activity on stellar and halo mass in Semi-Analytic Models

AGN feedback is believed to play an important role in shaping a variety of observed galaxy properties, as well as the evolution of their stellar masses and star formation rates. In particular, in the current theoretical paradigm of galaxy formation, AGN feedback is believed to play a crucial role in regulating the levels of activity in galaxies, in relatively massive halos at low redshift. Only in recent years, however, has detailed statistical information on the dependence of galaxy activity on stellar mass, parent halo mass and hierarchy has become available. In this paper, we compare the fractions of galaxies belonging to different activity classes (star-forming, AGN and radio active) with predictions from four different and independently developed semi-analytical models. We adopt empirical relations to convert physical properties into observables (H_alpha emission lines, OIII line strength and radio power). We demonstrate that all models used in this study reproduce the overall distributions of galaxies belonging to different activity classes as a function of stellar mass and halo mass: star forming galaxies and the strongest radio sources are preferentially associated with low-mass and high-mass galaxies/halos respectively. However, model predictions differ from observational measurements in a number of ways. All models used in our study predict that almost every >1.e12 Msun dark matter halo and/or >1.e11 Msun galaxy should host a bright radio source, while only a small fraction of galaxies belong to this class in the data. In addition, radio brightness is expected to depend strongly on the mass of the parent halo mass in the models, while strong and weak radio galaxies are found in similar environments in data. Our results highlight that the distribution of AGN as a function of stellar mass provides one of the most promising discriminants between different gas accretion schemes.

preprint2010arXiv

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

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

preprint2009arXiv

Luminosity function and radial distribution of Milky Way Satellites in a LCDM Universe

We study the luminosity function and the radial distribution of satellite galaxies within Milky Way sized haloes as predicted in Cold Dark Matter based models of galaxy formation, making use of numerical N-body techniques as well as three different semi-analytic model (SAMs) galaxy formation codes. We extract merger trees from very high-resolution dissipationless simulations of four Galaxy-sized DM haloes, and use these as common input for the semi-analytic models. We present a detailed comparison of our predictions with the observational data recently obtained on the Milky Way satellite luminosity function (LF). We find that semi-analytic models with rather standard astrophysical ingredients are able to reproduce the observed luminosity function over six orders of magnitude in luminosity, down to magnitudes as faint as M_V=-2. We also perform a comparison with the actual observed number of satellites as a function of luminosity, by applying the selection criteria of the SDSS survey to our simulations instead of correcting the observations for incompleteness. Using this approach we again find good agreement for both the luminosity and radial distributions of MW satellites. We investigate which physical processes in our models are responsible for shaping the predicted satellite LF, and find that tidal destruction, suppression of gas infall by a photo-ionizing background, and supernova feedback all make important contributions. We conclude that the number and luminosity of Milky Way satellites can be naturally accounted for within the (Lambda) Cold Dark Matter paradigm, and this should no longer be considered a problem.

preprint2005arXiv

Reconstructing the Triaxial Shapes of Dark Matter Halos from the Anisotropic Spatial Distributions of their Substructures in the Concordance Cosmology

We develop an algorithm to reconstruct the triaxial shapes of dark matter halos from the anisotropic spatial distributions of their substructures for the concordance background cosmology. First, we construct an analytic model for the anisotropic spatial distribution of dark halo substructures under the assumption that the tidal field with non-zero trace in the triaxial mass distribution of the host halos generates the substructure bulk motions toward the major principal axes of the hosts. Our analytic model implies that the degree of anisotropy depends sensitively on the triaxiality of the host halos as well as the correlation between the substructure locations and the tidal shear field. Second, we set the axis ratios of the triaxial host halos as free parameters in the analytic model, fit the model to the numerical results from high-resolution N-body simulation of the concordance cosmology, and reconstruct the two axis-ratios of the host halos from the best-fit values of the free parameters. The comparison of the reconstructed axis ratios with the numerical results reveals a good agreement. Finally, we conclude that our analytic model may provide a physical understanding of the anisotropic spatial distribution of dark halo substructures and a new way to reconstruct in principle the triaxial shapes of dark matter halos from the observables.