Search

The population of satellite galaxies in a host galaxy is a combination of the cumulative accretion of subhaloes and their associated star formation efficiencies, therefore, the luminosity distribution of satellites provides valuable information of both dark matter properties and star formation physics. Recently, the luminosity function of satellites in nearby Milky Way-mass galaxies has been well measured to satellites as faint as Leo I with $M_{V} \sim -8$. In addition to the finding of the diversity in the satellite luminosity functions, it has been noticed that there is a big gap among the magnitude of satellites in some host galaxies, such as M101, where the gap is around 5 in magnitude, noticeably larger than the prediction from the halo abundance matching method. The reason of this gap is still unknown. In this paper, we use a semi-analytical model of galaxy formation, combined with high-resolution N-body simulation, to investigate the probability and origin of such big gap in M101-alike galaxies. We found that, although M101 analogues are very rare with probability of \sim 0.1%-0.2% in the local universe, their formation is a natural outcome of the CDM model. The gap in magnitude is mainly due to the mass of the accreted subhaloes, not from the stochastic star formation in them. We also found that the gap is correlated with the total satellite mass and host halo mass. By tracing the formation history of M101 type galaxies, we find that they likely formed after $z \sim 1$ due to the newly accreted bright satellites. The gap is not in a stable state, and it will disappear in ~7 Gyr due to mergers of bright satellites with the central galaxy., Comment: 8 pages, 11 figures, accepted by MNRAS

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., 22 pages, 18 figures, 2 tables

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 ($\lambda_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., Comment: Accepted for publication in MNRAS. 15 pages, 10 figures

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., Comment: Accepted by MNRAS

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 $, Comment: 15 pages, 13 figures, accepted to ApJ

Observational work have shown that the two brightest satellite galaxies of the Milky Way (MW), the Large Magellanic Cloud (LMC) and the Small Magellanic Cloud (SMC), are rare amongst MW analogues. It is then interesting to know whether the presence of massive satellite has any effect on the whole satellite population in MW analogues. In this article, we investigate this problem using a semi-analytical model combined with the Millennium-II Simulation. MW-analogous galaxies are defined to have similar stellar mass or dark matter halo mass to the MW. We find that, in the first case, the halo mass is larger and there are, on average, twice as many satellites in Milky Way analogs if there is a massive satellite galaxy in the system. This is mainly from the halo formation bias. The difference is smaller if MW analogues are selected using halo mass. We also find that the satellites distribution is slightly asymmetric, being more concentrated on the line connecting the central galaxy and the massive satellite and that, on average, LMC have brought in 14.7 satellite galaxies with $M_{r}, 9 pages, 10 figures, published in MNRAS

Both simulation and observational data have shown that the spin and shape of dark matter halos are correlated with their nearby large-scale environment. As structure formation on different scales is strongly coupled, it is trick to disentangle the formation of halo with the large-scale environment, making it difficult to infer which is the driving force for the correlation between halo spin/shape with the large-scale structure. In this paper, we use N-body simulation to produce twin Universes that share the same initial conditions on small scales but different on large scales. This is achieved by changing the random seeds for the phase of those k modes smaller than a given scale in the initial conditions. In this way, we are able to disentangle the formation of halo and large-scale structure, making it possible to investigate how halo spin and shape correspond to the change of environment on large scales. We identify matching halo pairs in the twin simulations as those sharing the maximum number of identical particles within each other. Using these matched halo pairs, we study the cross match of halo spin and their correlation with the large-scale structure. It is found that when the large-scale environment changes (eigenvector) between the twin simulations, the halo spin has to rotate accordingly, although not significantly, to maintain the universal correlation seen in each simulation. Our results suggest that the large-scale structure is the main factor to drive the correlation between halo properties and their environment., Comment: 12 pages, 8 figures, accepted for publication in ApJ

Most cosmological structures in the universe spin. Although structures in the universe form on a wide variety of scales from small dwarf galaxies to large super clusters, the generation of angular momentum across these scales is poorly understood. We have investigated the possibility that filaments of galaxies - cylindrical tendrils of matter hundreds of millions of light-years across, are themselves spinning. By stacking thousands of filaments together and examining the velocity of galaxies perpendicular to the filament's axis (via their red and blue shift), we have found that these objects too display motion consistent with rotation making them the largest objects known to have angular momentum. The strength of the rotation signal is directly dependent on the viewing angle and the dynamical state of the filament. Just as it is easiest to measure rotation in a spinning disk galaxy viewed edge on, so too is filament rotation clearly detected under similar geometric alignment. Furthermore, the mass of the haloes that sit at either end of the filaments also increases the spin speed. The more massive the haloes, the more rotation is detected. These results signify that angular momentum can be generated on unprecedented scales., 12 pages, 6 figures, to appear in the 14 June issue of Nature Astronomy, a version after one round of reviews. The final authenticated version is available online at https://www.nature.com/articles/s41550-021-01380-6

Satellites are not randomly distributed around their central galaxies but show polar and planar structures. In this paper, we investigate the axis-asymmetry or lopsidedness of satellite galaxy distributions around isolated galaxies in a hydrodynamic cosmological simulation. We find a statistically significant lopsided signal by studying the angular distribution of the satellite galaxies' projected positions around isolated central galaxies in a two-dimensional plane. The signal is dependent on galaxy mass, color and large-scale environment. Satellites that inhabit low-mass blue hosts, or located further from the hosts show the most lopsided signal. Galaxy systems with massive neighbors exhibit stronger lopsidedness. This satellite axis-asymmetry signal also decreases as the universe evolves. Our findings are in agreement with recent observational results, and they provide a useful perspective for studying galaxy evolution, especially on the satellite accretion, internal evolution and interaction with the cosmic large-scale structure., Comment: 13 pages, 8 figures, accepted by ApJ

The Hot Universe Baryon Surveyor (HUBS) mission is proposed to study "missing" baryons in the universe. Unlike dark matter, baryonic matter is made of elements in the periodic table, and can be directly observed through the electromagnetic signals that it produces. Stars contain only a tiny fraction of the baryonic matter known to be present in the universe. Additional baryons are found to be in diffuse (gaseous) form, in or between galaxies, but a significant fraction has not yet been seen. The latter (missing baryons) are thought to be hiding in low-density warm-hot ionized medium (WHIM), based on results from theoretical studies and recent observations, and be distributed in the vicinity of galaxies (i.e., circum-galactic medium) and between galaxies (i.e., intergalactic medium). Such gas would radiate mainly in the soft X-ray band and the emission would be very weak, due to its very low density. HUBS is optimized to detect the X-ray emission from the hot baryons in the circum-galactic medium, and thus fill a void in observational astronomy. The goal is not only to detect the missing baryons, but to characterize their physical and chemical properties, as well as to measure their spatial distribution. The results would establish the boundary conditions for understanding galaxy evolution. Though highly challenging, detecting missing baryons in the intergalactic medium could be attempted, perhaps in the outskirts of galaxy clusters, and could shed significant light on the large-scale structures of the universe. The current design of HUBS will be presented, along with the status of technology development., Comment: 12 pages, 7 figures, Space Telescopes and Instrumentation 2020: Ultraviolet to Gamma Ray, eds. J.A. den Herder, S. Nikzad, & K. Nakazawa, Proc. of SPIE Vol. 11444, 114442S

We use hydrodynamical cosmological simulations to show that it is possible to create, via tidal interactions, galaxies lacking dark matter in a dark matter dominated universe. We select dwarf galaxies from the NIHAO project, obtained in the standard Cold Dark Matter model and use them as initial conditions for simulations of satellite-central interactions. After just one pericentric passage on an orbit with a strong radial component, NIHAO dwarf galaxies can lose up to 80 per~cent of their dark matter content, but, most interestingly, their central ($\approx 8$~kpc) dark matter to stellar ratio changes from a value of ${\sim}25$, as expected from numerical simulations and abundance matching techniques, to roughly unity as reported for NGC1052-DF2 and NGC1054-DF4. The stellar velocity dispersion drops from ${\sim}30$ ${\rm km\,s^{-1}}$ before infall to values as low as $6\pm 2$~ ${\rm km\,s^{-1}}$. These, and the half light radius around 3 kpc, are in good agreement with observations from van Dokkum and collaborators. Our study shows that it is possible to create a galaxy "without" dark matter starting from typical dwarf galaxies formed in a dark matter dominated universe, provided they live in a dense environment., 8 pages, 11 figures. Two extra figures added, accepted for publication in MNRAS

We have used the NIHAO simulations to explore how supernovae (SNe) affect star formation in galaxies. We find that SN feedback operates on all scales from the interstellar medium (ISM) to several virial radii. SNe regulate star formation by preventing condensation of HI into H$_2$ and by moving cold neutral gas to the hot HII phase. The first effect explains why the cold neutral gas in dwarf galaxies forms stars inefficiently. The second maintains the hot ISM of massive galaxies (HII vents out at lower masses). At $v_{\rm vir}\simeq 67{\rm\,km\,s}^{-1}$, the outflow rate follows the relation: $\dot{M}_{\rm out}=23\,(v_{\rm vir}/67{\rm\,km\,s}^{-1})^{-4.6}\,{\rm SFR}$. $20\%$ to $70\%$ of the gas expelled from galaxies escapes from the halo (ejective feedback) but outflows are dominated by cold swept-up gas, most of which falls back onto the galaxy on a $\sim 1\,$Gyr timescale. This `fountain feedback' reduces the masses of galaxies by a factor of two to four, since gas spends half to three quarter of its time in the fountain. Less than $10\%$ of the ejected gas mixes with the hot circumgalactic medium and this gas is usually not reaccreted. On scales as large as $6r_{\rm vir}$, galactic winds divert the incoming gas from cosmic filaments and prevent if from accreting onto galaxies (pre-emptive feedback). This process is the main reason for the low baryon content of ultradwarves., Comment: Submitted for publication in MNRAS

International audience; It has been shown, both in simulations and observationally, that the tidal field of a large galaxy can torque its satellites such that the major axis of satellite galaxies points towards their hosts. This so-called ‘shape alignment’ has been observed in isolated Milky Way-like galaxies but not in ‘Local Group’-like pairs. In this study, we investigate the shape alignment of satellite galaxies in galaxy pairs similar to the Local Group identified in the Sloan Digital Sky Survey Data Release 13 (SDSS DR13). By stacking tens of thousands of satellite galaxies around primary galaxy pairs, we find two statistically strong alignment signals. (1) The major axes of satellite galaxies located in the (projected) area between two primaries (the facing region) tend to be perpendicular to the line connecting the satellite to its host (tangential alignment), while (2) the major axes of satellite galaxies located in regions away from the other host (the away region) tend to be aligned with the line connecting the satellite to its host (radial alignment). These alignments are confirmed at ∼5σ levels. The alignment signal increases with increasing primary brightness, decreasing pair separation, and decreasing satellite distance. The alignment signal is also found to be stronger in filamentary environments. These findings will shed light on understanding the mechanisms of how satellite galaxies are affected by the tidal field in galaxy pairs and will be useful for investigating galaxy intrinsic alignment in the analyses of weak gravitational lensing.

Both simulations and observations have confirmed that the spin of haloes/galaxies is correlated with the large scale structure (LSS) with a mass dependence such that the spin of low-mass haloes/galaxies tend to be parallel with the LSS, while that of massive haloes/galaxies tend to be perpendicular with the LSS. It is still unclear how this mass dependence is built up over time. We use N-body simulations to trace the evolution of the halo spin-LSS correlation and find that at early times the spin of all halo progenitors is parallel with the LSS. As time goes on, mass collapsing around massive halo is more isotropic, especially the recent mass accretion along the slowest collapsing direction is significant and it brings the halo spin to be perpendicular with the LSS. Adopting the $fractional$ $anisotropy$ (FA) parameter to describe the degree of anisotropy of the large-scale environment, we find that the spin-LSS correlation is a strong function of the environment such that a higher FA (more anisotropic environment) leads to an aligned signal, and a lower anisotropy leads to a misaligned signal. In general, our results show that the spin-LSS correlation is a combined consequence of mass flow and halo growth within the cosmic web. Our predicted environmental dependence between spin and large-scale structure can be further tested using galaxy surveys., Comment: 9 pages, 7 figures, 2 tables, Accepted for publication in MNRAS

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., 10 pages, 7 figures, accepted by ApJ

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., Comment: 5 papes, 3 figures, accepted for publication in MNRAS Letters

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., Published in New Astronomy, https://doi.org/10.1016/j.newast.2020.101405

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 ($\alpha$) is fairly constant ($\Delta\alpha \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)., Comment: 5 pages, 4 figures, accepted for publication in MNRAS Letters

We use a sample of 83 high-resolution cosmological zoom-in simulations and a semi-analytic model (SAM) to study the stochasticity of galaxy formation in haloes ranging from dwarf to Milky Way masses. Our simulated galaxies reproduce the observed inefficiency of galaxy formation as expressed through the stellar, gas and baryonic Tully-Fisher relations. For HI velocities in the range (70 less than or similar to V less than or similar to 220 km s(-1)), the scatter is just 0.08 to 0.14 dex, consistent with the observed intrinsic scatter at these scales. At low velocities (20 less than or similar to V less than or similar to 70 km s(-1)), the simulated scatter is 0.2-0.25 dex, which could be tested with future observations. The scatter in the stellar mass versus dark halo velocity relation is constant for 30 less than or similar to V less than or similar to 180 km s(-1), and smaller (similar or equal to 0.17 dex) when using the maximum circular velocity of the dark-matter-only simulation, V-max(DMO), compared to the virial velocity (V-200 or V-200(DMO)). The scatter in stellar mass is correlated with halo concentration, and is minimized when using a circular velocity at a fixed fraction of the virial radius similar or equal to 0.4R(200) or with V-alpha = V-200(DMO)(V-max(DMO) / V-200(DMO))(alpha) with alpha similar or equal to 0.7, consistent with constraints from halo clustering. Using the SAM we show the correlation between halo formation time and concentration is essential in order to reproduce this result. This uniformity in galaxy formation efficiency we see in our hydrodynamical simulations and an SAM proves the simplicity and self-regulating nature of galaxy formation in a Lambda cold dark matter (Lambda CDM) universe.

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., Comment: 17 pages, 12 figures, MNRAS in press

Following Cui et al. 2018 (hereafter Paper I) on the classification of large-scale environments (LSE) at z = 0, we push our analysis to higher redshifts and study the evolution of LSE and the baryon distributions in them. Our aim is to investigate how baryons affect the LSE as a function of redshift. In agreement with Paper I, the baryon models have negligible effect on the LSE over all investigated redshifts. We further validate the conclusion obtained in Paper I that the gas web is an unbiased tracer of total matter -- even better at high redshifts. By separating the gas mainly by temperature, we find that about 40 per cent of gas is in the so-called warm-hot intergalactic medium (WHIM). This fraction of gas mass in the WHIM decreases with redshift, especially from z = 1 (29 per cent) to z = 2.1 (10 per cent). By separating the whole WHIM gas mass into the four large-scale environments (i.e. voids, sheets, filaments, and knots), we find that about half of the WHIM gas is located in filaments. Although the total gas mass in WHIM decreases with redshift, the WHIM mass fractions in the different LSE seem unchanged.

We analyze 88 independent high-resolution cosmological zoom-in simulations of disk galaxies in the NIHAO simulations suite to explore the connection between the atomic gas fraction and angular momentum of baryons throughout cosmic time. The study is motivated by the analytic model of \citet{obreschkow16}, which predicts a relation between the atomic gas fraction $f_{\rm atm}$ and the global atomic stability parameter $q \equiv j\sigma / (GM)$, where $M$ and $j$ are the mass and specific angular momentum of the galaxy (stars+cold gas) and $\sigma$ is the velocity dispersion of the atomic gas. We show that the simulated galaxies follow this relation from their formation ($z\simeq4$) to present within $\sim 0.5$ dex. To explain this behavior, we explore the evolution of the local Toomre stability and find that $90\%$--$100\%$ of the atomic gas in all simulated galaxies is stable at any time. In other words, throughout the entire epoch of peak star formation until today, the timescale for accretion is longer than the timescale to reach equilibrium, thus resulting in a quasi-static equilibrium of atomic gas at any time. Hence, the evolution of $f_{\rm atm}$ depends on the complex hierarchical growth history primarily via the evolution of $q$. An exception are galaxies subject to strong environmental effects., Comment: 12 pages, 7 figures; accepted to ApJ

We develop a new method StarGO (Stars' Galactic Origin) to identify the galactic origins of halo stars using their kinematics. Our method is based on self-organizing map (SOM), which is one of the most popular unsupervised learning algorithms. StarGO combines SOM with a novel adaptive group identification algorithm with essentially no free parameters. In order to evaluate our model, we build a synthetic stellar halo from mergers of nine satellites in the Milky Way. We construct the mock catalogue by extracting a heliocentric volume of 10 kpc from our simulations and assigning expected observational uncertainties corresponding to bright stars from Gaia DR2 and LAMOST DR5. We compare the results from StarGO against that from a Friends-of-Friends (FoF) based method in the space of orbital energy and angular momentum. We show that StarGO is able to systematically identify more satellites and achieve higher number fraction of identified stars for most of the satellites within the extracted volume. When applied to data from Gaia DR2, StarGO will enable us to reveal the origins of the inner stellar halo in unprecedented detail., 11 pages, 7 figures, Accepted for publication in ApJ

The convergence peaks, constructed from galaxy shape measurement in weak lensing, is a powerful probe of cosmology as the peaks can be connected with the underlined dark matter haloes. However the capability of convergence peak statistic is affected by the noise in galaxy shape measurement, signal to noise ratio as well as the contribution from the projected mass distribution from the large-scale structures along the line of sight (LOS). In this paper we use the ray-tracing simulation on a curved sky to investigate the correspondence between the convergence peak and the dark matter haloes at the LOS. We find that, in case of no noise and for source galaxies at $z_{\rm s}=1$, more than $65\%$ peaks with $\text{SNR} \geq 3$ (signal to noise ratio) are related to more than one massive haloes with mass larger than $10^{13} {\rm M}_{\odot}$. Those massive haloes contribute $87.2\%$ to high peaks ($\text{SNR} \geq 5$) with the remaining contributions are from the large-scale structures. On the other hand, the peaks distribution is skewed by the noise in galaxy shape measurement, especially for lower SNR peaks. In the noisy field where the shape noise is modelled as a Gaussian distribution, about $60\%$ high peaks ($\text{SNR} \geq 5$) are true peaks and the fraction decreases to $20\%$ for lower peaks ($ 3 \leq \text{SNR} < 5$). Furthermore, we find that high peaks ($\text{SNR} \geq 5$) are dominated by very massive haloes larger than $10^{14} {\rm M}_{\odot}$., 13 pages, 11 figures, 4 tables, accepted for publication in MNRAS. Our mock galaxy catalog is available upon request by email to the author (chengliangwei@pmo.ac.cn)

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., 15 papes, 11 figures, accepted for publication in MNRAS

We use the NIHAO simulations to investigate the effects of baryonic physics on the time evolution of Dark Matter central density profiles. The sample is made of $\approx 70$ independent high resolution hydrodynamical simulations of galaxy formation and covers a wide mass range: 1e10< Mhalo, Comment: 11 pages, 13 figures. Corrected typo in table 2 (middle row) with respect to the version published in MNRAS

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., 14 pages, 8 figures, accepted for publication in MNRAS

Using a series of cosmological simulations that includes one dark-matter-only (DM-only) run, one gas cooling-star formation-supernovae feedback (CSF) run and one that additionally includes feedback from active galactic nuclei (AGNs), we classify the large-scale structures with both a velocity-shear-tensor code (Vweb) and a tidal-tensor code (Pweb). We find that the baryonic processes have almost no impact on large-scale structures -- at least not when classified using aforementioned techniques. More importantly, our results confirm that the gas component alone can be used to infer the filamentary structure of the Universe practically un-biased, which could be applied to cosmology constrains. In addition, the gas filaments are classified with its velocity (Vweb) and density (Pweb) fields, which can theoretically connect to the radio observations, such as HI surveys. This will help us to bias-freely link the radio observations with DM distributions at large scale., 13 pages, 8 figures, matched to the MNRAS accepted version (V2). Matched to the the version after the Erratum. No conclusion is changed (V3). Corrections for the references papers (V4)

The correlation between the spins of dark matter halos and the large-scale structure (LSS) has been studied in great detail over a large redshift range, while investigations of galaxies are still incomplete. Motivated by this point, we use the state-of-the-art hydrodynamic simulation, Illustris-1, to investigate mainly the spin--LSS correlation of galaxies at redshift of $z=0$. We mainly find that the spins of low-mass, blue, oblate galaxies are preferentially aligned with the slowest collapsing direction ($e_3$) of the large-scale tidal field, while massive, red, prolate galaxy spins tend to be perpendicular to $e_3$. The transition from a parallel to a perpendicular trend occurs at $\sim10^{9.4} M_{\odot}/h$ in the stellar mass, $\sim0.62$ in the g-r color, and $\sim0.4$ in triaxiality. The transition stellar mass decreases with increasing redshifts. The alignment was found to be primarily correlated with the galaxy stellar mass. Our results are consistent with previous studies both in N-body simulations and observations. Our study also fills the vacancy in the study of the galaxy spin--LSS correlation at $z=0$ using hydrodynamical simulations and also provides important insight to understand the formation and evolution of galaxy angular momentum., Comment: 9 pages, 6 figures, 1 table. Accepted for publication in ApJ, match the proof version

Both simulations and observations have found that the spin of halo/galaxy is correlated with the large scale environment, and particularly the spin of halo flips in filament. A consistent picture of halo spin evolution in different environments is still lacked. Using N-body simulation we find that halo spin with its environment evolves continuously from sheet to cluster, and the flip of halo spin happens both in filament and nodes. For the flip in filament can be explained by halo formation time and the migrating time when its environment changes from sheet to filament. For low-mass haloes, they form first in sheets and migrate into filaments later, so their mass and spin growth inside filament are lower, and the original spin is still parallel to filament. For massive haloes, they migrate into filaments first, and most of their mass and spin growth are obtained in filaments, so the resulted spin is perpendicular to filament. Our results well explain the overall evolution of cosmic web in the cold dark matter model and can be tested using high-redshift data. The scenario can also be tested against alternative models of dark matter, such as warm/hot dark matter, where the structure formation will proceed in a different way., 5 pages, 2 figures, Accepted by MNRAS Letters

We use ROSAT All Sky Survey broad-band X-ray images and the optical clusters identified from Sloan Digital Sky Survey Data Release 7 to estimate the X-ray luminosities around similar to 65 000 candidate clusters with masses greater than or similar to 10(13) h(-1) M-circle dot based on an optical to X-ray (OTX) code we develop. We obtain a catalogue with X-ray luminosity for each cluster. This catalogue contains 817 clusters (473 at redshift z 3 in X-ray detection. We find about 65 per cent of these X-ray clusters have their most massive member located near the X-ray flux peak; for the rest 35 per cent, the most massive galaxy is separated from the X-ray peak, with the separation following a distribution expected from a Navarro-Frenk-White 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 haloes with masses greater than or similar to 10(14) h(-1) M-circle dot contain a larger fraction of red satellite galaxies when they are brighter in X-ray. An opposite trend is found in central galaxies in relative low-mass haloes with masses less than or similar to 10(14) h(-1) M-circle dot where X-ray brighter clusters have smaller fraction of red central galaxies. Clusters with masses greater than or similar to 10(14) h(-1) M-circle dot that are strong X-ray emitters contain many more low-mass satellite galaxies than weak X-ray emitters. These results are also confirmed by checking X-ray clusters of similar X-ray luminosities but having different characteristic stellar masses. A cluster catalogue containing the optical properties of member galaxies and the X-ray luminosity is available at http://gax.shao.ac.cn/data/Group.html.

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., Comment: 5 pages, 3 figures, accepted for publication in MNRAS Letters

We examine the quenched fraction of central and satellite galaxies as a function of galaxy stellar mass, halo mass, and the matter density of their large scale environment. Matter densities are inferred from our ELUCID simulation, a constrained simulation of local Universe sampled by SDSS, while halo masses and central/satellite classification are taken from the galaxy group catalog of Yang et al. The quenched fraction for the total population increases systematically with the three quantities. We find that the `environmental quenching efficiency', which quantifies the quenched fraction as function of halo mass, is independent of stellar mass. And this independence is the origin of the stellar mass-independence of density-based quenching efficiency, found in previous studies. Considering centrals and satellites separately, we find that the two populations follow similar correlations of quenching efficiency with halo mass and stellar mass, suggesting that they have experienced similar quenching processes in their host halo. We demonstrate that satellite quenching alone cannot account for the environmental quenching efficiency of the total galaxy population and the difference between the two populations found previously mainly arises from the fact that centrals and satellites of the same stellar mass reside, on average, in halos of different mass. After removing these halo-mass and stellar-mass effects, there remains a weak, but significant, residual dependence on environmental density, which is eliminated when halo assembly bias is taken into account. Our results therefore indicate that halo mass is the prime environmental parameter that regulates the quenching of both centrals and satellites., Comment: 21 pages, 16 figures, submitted to ApJ

Semi-analytic models of galaxy formation are powerful tools to study the evolution of galaxy population in a cosmological context. However, most models over-predict the number of low-mass galaxies at high redshifts and the color of model galaxies are not right in the sense that low-mass satellite galaxies are too red and centrals are too blue. The recent version of the L-Galaxies model by Henriques et al.(H15) is a step forward to solve these problems by reproducing the evolution of stellar mass function and the overall fraction of red galaxies. In this paper we compare the two model predictions of L-Galaxies (the other is Guo et al. , G13) to the SDSS data in detail. We find that in the H15 model the red fraction of central galaxies now agrees with the data due to their implementation of strong AGN feedback, but the stellar mass of centrals in massive haloes is now slightly lower than the data. For satellite galaxies, the red fraction of low-mass galaxies ($\log M_{*}/M_{\odot} < 10$) also agrees with the data, but the color of massive satellites ($10 < \log M_{*}/M_{\odot} < 11$) is slightly bluer. The correct color of centrals and bluer color of massive satellites indicate that the quenching in massive satellites are not strong enough. We also find that there are too much red spirals and less bulge-dominated galaxies in both H15 and G13 models. Our results suggest that additional mechanisms, such as more minor merger or disk instability, are needed to slightly increase the stellar mass of central galaxy in massive galaxies, mainly in the bulge component, and the bulge dominated galaxy will be quenched or then be quenched by any other mechanisms., 12 pages, 6 figures, accepted for publication in RAA

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., 19 pages, 17 figures, analysis strongly improved, main conclusions unchanged, accepted for publication in MNRAS

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., 7 pages, 4 figures, minor changes to match the version accepted on MNRAS Letters

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})/(\Omega_b/\Omega_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., Comment: 20 pages, 15 figures, accepted to MNRAS

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, Comment: 5 pages, 3 figures, accepted for publication in MNRAS

Using N-body simulations and galaxy formation models, we study the galaxy stellar mass correlation and the two-point auto-correlation. The simulations are run with cosmological parameters from the WMAP first, third and seven year results, which mainly differ in the perturbation amplitude of \sigma_{8}. The stellar mass of galaxies are determined using either a semi-analytical galaxy formation model or a simple empirical abundance matching method. Compared to the SDSS DR7 data at z=0 and the DEEP2 results at z=1, we find that the predicted galaxy clusterings from the semi-analytical model are higher than the data at small scales, regardless of the adopted cosmology. Conversely, the abundance matching method predicts good agreement with the data at both z=0 and z=1 for high \sigma_8 cosmologies (WMAP1 & WMAP7), but the predictions from a low \sigma_8 cosmology (WMAP3) are significantly lower than the data at z=0. We find that the excess clustering at small-scales in the semi-analytical model mainly arises from satellites in massive haloes, indicating that either the star formation is too efficient in low-mass haloes or tidal stripping is too inefficient at high redshift. Our results show that galaxy clustering is strongly affected by the models for galaxy formation, thus can be used to constrain the baryonic physics. The weak dependence of galaxy clustering on cosmological parameters makes it difficult to constrain the WMAP1 and WMAP7 cosmologies.

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.

We propose a new strategy of finding strongly-lensed supernovae (SNe) by monitoring known galaxy-scale strong-lens systems. Strongly lensed SNe are potentially powerful tools for the study of cosmology, galaxy evolution, and stellar populations, but they are extremely rare. By targeting known strongly lensed starforming galaxies, our strategy significantly boosts the detection efficiency for lensed SNe compared to a blind search. As a reference sample, we compile the 128 galaxy-galaxy strong-lens systems from the Sloan Lens ACS Survey (SLACS), the SLACS for the Masses Survey, and the Baryon Oscillation Spectroscopic Survey Emission-Line Lens Survey. Within this sample, we estimate the rates of strongly-lensed Type Ia SN (SNIa) and core-collapse SN (CCSN) to be $1.23 \pm 0.12$ and $10.4 \pm 1.1$ events per year, respectively. The lensed SN images are expected to be widely separated with a median separation of 2 arcsec. Assuming a conservative fiducial lensing magnification factor of 5 for the most highly magnified SN image, we forecast that a monitoring program with a single-visit depth of 24.7 mag (5$\sigma$ point source, $r$ band) and a cadence of 5 days can detect 0.49 strongly-lensed SNIa event and 2.1 strongly-lensed CCSN events per year within this sample. Our proposed targeted-search strategy is particularly useful for prompt and efficient identifications and follow-up observations of strongly-lensed SN candidates. It also allows telescopes with small field of views and limited time to efficiently discover strongly-lensed SNe with a pencil-beam scanning strategy., Comment: 14 pages, 5 figures, ApJ in press

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).

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

The alignment between satellites and central galaxies has been studied in detail both in observational and theoretical works. The widely accepted fact is that the satellites preferentially reside along the major axis of their central galaxy. However, the origin and large-scale environment dependence of this alignment are still unknown. In an attempt to figure out those, we use data constructed from SDSS DR7 to investigate the large-scale environmental dependence of this alignment with emphasis on examining the alignments' dependence on the colour of the central galaxy. We find a very strong large-scale environmental dependence of the satellite-central alignment in groups with blue centrals. Satellites of blue centrals in knots are preferentially located perpendicular to the major axis of the centrals, and the alignment angle decreases with environment namely when going from knots to voids. The alignment angle strongly depend on the ${}^{0.1}(g-r)$ colour of centrals. We suggest that the satellite-central alignment is the result of a competition between satellite accretion within large scale-structure and galaxy evolution inside host haloes. For groups containing red central galaxies, the satellite-central alignment is mainly determined by the evolution effect, while for blue central dominated groups, the effect of large-scale structure plays a more important role, especially in knots. Our results provide an explanation for how the satellite-central alignment forms within different large-scale environments. The perpendicular case in groups and knots with blue centrals may also provide insight into understanding similar polar arrangements such the formation of the Milky Way and Centaurus A's satellite system., 10 pages, 5 figures and 2 tables, accepted for publication in ApJ

It has been pointed out that the Galactic satellites all have a common mass around 1e7 Msun within 300 pc (M03), while they span almost four order of magnitudes in luminosity (Mateo et al. 1993, Strigari et al. 2008). It is argued that this may reflect a specific scale for galaxy formation or a scale for dark matter clustering. Here we use numerical simulations coupled with a semi-analytic model for galaxy formation, to predict the central mass and luminosity of galactic satellites in the LCDM model. We show that this common mass scale can be explained within the Cold Dark Matter scenario when the physics of galaxy formation is taken into account. The narrow range of M03 comes from the narrow distribution of circular velocities at time of accretion (peaking around 20 km/s) for satellites able to form stars and the not tight correlation between halo concentration and circular velocity. The wide range of satellite luminosities is due to a combination of the mass at time of accretion and the broad distribution of accretion redshifts for a given mass. This causes the satellites baryonic content to be suppressed by photo-ionization to very different extents. Our results favor the argument that the common mass M03 reflects a specific scale (circular velocity ~ 20 km/s) for star formation.

Detailed studies of galaxy formation require clear definitions of the structural components of galaxies. Precisely defined components also enable better comparisons between observations and simulations. We use a subsample of eighteen cosmological zoom-in simulations from the NIHAO project to derive a robust method for defining stellar kinematic discs in galaxies. Our method uses Gaussian Mixture Models in a 3D space of dynamical variables. The NIHAO galaxies have the right stellar mass for their halo mass, and their angular momenta and S��rsic indices match observations. While the photometric disc-to-total ratios are close to 1 for all the simulated galaxies, the kinematic ratios are around ~0.5. Thus, exponential structure does not imply a cold kinematic disc. Above log(M*)~9.5, the decomposition leads to thin discs and spheroids that have clearly different properties, in terms of angular momentum, rotational support, ellipticity, [Fe/H] and [O/Fe]. At log(M*), Accepted for publication in MNRAS

We study the connection between the observed star formation rate-stellar mass (SFR-$M_*$) relation and the evolution of the stellar mass function (SMF) by means of a Subhalo Abundance Matching technique coupled to merger trees extracted from a N-body simulation. Our approach, which considers both galaxy mergers and stellar stripping, is to force the model to match the observed SMF at redshift $z>2$, and let it evolve down to the present time according to the observed (SFR-$M_*$) relation. In this study, we use two different sets of SMFs and two SFR-$M_*$ relations: a simple power law and a relation with a mass-dependent slope. Our analysis shows that the evolution of the SMF is more consistent with a SFR-$M_*$ relation with a mass-dependent slope, in agreement with predictions from other models of galaxy evolution and recent observations. In order to fully and realistically describe the evolution of the SMF, both mergers and stellar stripping must be considered, and we find that both have almost equal effects on the evolution of SMF at the massive end. Taking into account the systematic uncertainties in the observed data, the high-mass end of the SMF obtained by considering stellar stripping results in good agreement with recent observational data from the Sloan Digital Sky Survey (SDSS). At $\log M_* < 11.2$, our prediction at z=0.1 is close to \citet{li-white09} data, but the high-mass end ($\log M_* > 11.2$) is in better agreement with \citet{dsouza15} data which account for more massive galaxies., 16 pages, 9 figures. Accepted by ApJ

Using a semi-analytical model for galaxy formation, combined with a large N-body simulation, we investigate the origin of the dichotomy among early-type galaxies. We find that boxy galaxies originate from mergers with a progenitor mass ratio $n < 2$ and with a combined cold gas mass fraction $F_{\rm cold} < 0.1$. Our model accurately reproduces the observed fraction of boxy systems as a function of luminosity and halo mass, for both central galaxies and satellites. After correcting for the stellar mass dependence, the properties of the last major merger of early-type galaxies are independent of their halo mass. This provides theoretical support for the conjecture of Pasquali et al (2007) that the stellar mass of an early-type galaxy is the main parameter that governs its isophotal shape. We argue that the observed dichotomy of early-type galaxies has a natural explanation within hierarchical structure formation, and does not require AGN feedback. Rather, we argue that it owes to the fact that more massive systems (i) have more massive progenitors, (ii) assemble later, and (iii) have a larger fraction of early-type progenitors. Each of these three trends causes the cold gas mass fraction of the progenitors of more massive early-types to be lower, so that their last major merger was dryer. Finally, our model predicts that (i) less than 10 percent of all early-type galaxies form in major mergers that involve two early-type progenitors, (ii) more than 95 percent of all boxy early-type galaxies with $M_* < 2 \times 10^{10} h^{-1} \Msun$ are satellite galaxies, and (iii) about 70 percent of all low mass early-types do not form a supermassive black hole binary at their last major merger. The latter may help to explain why low mass early-types have central cusps, while their massive counterparts have cores.

By comparing semi-analytic galaxy catalogues with data from the Sloan Digital Sky Survey (SDSS), we show that current galaxy formation models reproduce qualitatively the dependence of galaxy clustering and pairwise peculiar velocities on luminosity, but some subtle discrepancies with the data still remain. The comparisons are carried out by constructing a large set of mock galaxy redshift surveys that have the same selection function as the SDSS Data Release Four (DR4). The mock surveys are based on two sets of semi-analytic catalogues presented by Croton et al. and Kang et al. From the mock catalogues, we measure the redshift-space projected two-point correlation function, the power spectrum, and the pairwise velocity dispersion (PVD) in Fourier space and in configuration space, for galaxies in different luminosity intervals. We then compare these theoretical predictions with the measurements derived from the SDSS DR4. On large scales and for galaxies brighter than L*, both sets of mock catalogues agree well with the data. For fainter galaxies, however, both models predict stronger clustering and higher pairwise velocities than observed. We demonstrate that this problem can be resolved if the fraction of faint satellite galaxies in massive haloes is reduced by ~30% compared to the model predictions. A direct look into the model galaxy catalogues reveals that a signifcant fraction (15%) of faint galaxies ($-1810^{13}\msun$, and this population is predominantly red in colour. These faint red galaxies are responsible for the high PVD values of low-luminosity galaxies on small scales., Comment: 16 pages, 12 figures, accepted for publication in Monthly Notices, typos corrected

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., 6 pages, 3 figures, 2 tables. Accepted for publication in MNRAS letters