Search

One of the main challenges in galaxy formation that has emerged recently is the early assembly of massive galaxies. The observed number density and the maximum stellar mass ($M_{\star}$) of massive galaxies in the early Universe appear to be higher than model predictions, which may pose a serious problem to the LCDM cosmology. A major limitation in many previous studies is the large uncertainty in estimating $M_{\star}$ due to the lack of constraints in the rest-frame near-infrared part of the spectral energy distribution, which is critical to determining $M_{\star}$ accurately. Here we use data from a large JWST/MIRI survey in the PRIMER program to carry out a systematic analysis of massive galaxies at $z \sim 3-8$, leveraging photometric constraints at rest-frame $\gtrsim 1 \mu$m. We find a significant reduction in the number and mass densities of massive galaxies at $z > 5$ compared to earlier results that did not use the MIRI photometry. Within the standard $\Lambda$CDM cosmology, our results require a moderate increase in the baryon-to-star conversion efficiency ($\epsilon$) towards higher redshifts and higher $M_{\star}$. For the most massive galaxies at $z\sim 8$, the required $\epsilon$ is $\sim 0.3$, in comparison to $\epsilon \sim 0.14$ for typical low-redshift galaxies. Our findings are consistent with models assuming suppressed stellar feedback due to the high gas density and the associated short free-fall time expected for massive halos at high redshift., Comment: 23 pages, 10 figures, submitted

Using data from ALFALFA, xGASS, HI-MaNGA and the Sloan Digital Sky Survey (SDSS), we identify a sample of 47 "red but HI-rich"(RR) galaxies with $NUV-r > 5$ and unusually high HI-to-stellar mass ratios. We compare the optical properties and local environments between the RR galaxies and a control sample of "red and HI-normal"(RN) galaxies that are matched in stellar mass and color. The two samples are similar in the optical properties typical of massive red (quenched) galaxies in the local Universe. The RR sample tends to be associated with slightly lower-density environments and has lower clustering amplitudes and smaller neighbor counts at scales from several kiloparsecs to a few Megaparsecs. The results are consistent with the RR galaxies preferentially being located at the center of low-mass halos, with a median halo mass $\sim 10^{12}h^{-1}M_{\odot}$ compared to $\sim 10^{12.5}h^{-1}M_{\odot}$ for the RN sample. This result is confirmed by the SDSS group catalog which reveals a central fraction of 89% for the RR sample, compared to $\sim 60\%$ for the RN sample. If assumed to follow the HI size-mass relation of normal galaxies, the RR galaxies have an average HI-to-optical radius ratio of $R_{HI}/R_{90}\sim 4$, four times the average ratio for the RN sample. We compare our RR sample with similar samples in previous studies, and quantify the population of RR galaxies using the SDSS complete sample. We conclude that the RR galaxies form a unique but rare population, accounting for only a small fraction of the massive quiescent galaxy population. We discuss the formation scenarios of the RR galaxies., Comment: 25 pages, 12 figures, accepted for publication in ApJ, proof completed. Correct some typos in the text. Correct a typo in figure 2

We propose an efficient and robust method to estimate the halo concentration based on the first moment of the density distribution, which is $R_1\equiv \int_0^{r_{\rm vir}}4\pi r^3\rho(r)dr/M_{\rm vir}/r_{\rm vir}$. We find that $R_1$ has a monotonic relation with the concentration parameter of the NFW profile, and that a cubic polynomial function can fit the relation with an error $\lesssim 3\%$. Tests on ideal NFW halos show that the conventional NFW profile fitting method and the $V_{\rm max}/V_{\rm vir}$ method produce biased halo concentration estimation by $\approx 10\%$ and $\approx 30\%$, respectively, for halos with 100 particles. In contrast, the systematic error for our $R_1$ method is smaller than $0.5\%$ even for halos containing only 100 particles. Convergence tests on realistic halos in $N$-body simulations show that the NFW profile fitting method underestimates the concentration parameter for halos with $\lesssim 300$ particles by $\gtrsim 20\%$, while the error for the $R_1$ method is $\lesssim 8\%$. We also show other applications of $R_1$, including estimating $V_{\rm max}$ and the Einasto concentration $c_{\rm e}\equiv r_{\rm vir}/r_{-2}$. The calculation of $R_1$ is efficient and robust, and we recommend including it as one of the halo properties in halo catalogs of cosmological simulations., Comment: 17 pages, 9 + 9 figures. Main figure: Fig. 3. MNRAS accepted

We propose a novel method to quantify the assembly histories of dark matter halos with the redshift evolution of the mass-weighted spatial variance of their progenitor halos, i.e. the protohalo size history. We find that the protohalo size history for each individual halo at z~0 can be described by a double power-law function. The amplitude of the fitting function strongly correlates to the central-to-total stellar mass ratios of descendant halos. The variation of the amplitude of the protohalo size history can induce a strong halo assembly bias effect for massive halos. This effect is detectable in observation using the central-to-total stellar mass ratio as a proxy of the protohalo size. The correlation to the descendant central-to-total stellar mass ratio and the halo assembly bias effect seen in the protohalo size are much stronger than that seen in the commonly adopted half-mass formation time derived from the mass accretion history. This indicates that the information loss caused by the compression of halo merger trees to mass accretion histories can be captured by the protohalo size history. Protohalo size thus provides a useful quantity to connect protoclusters across cosmic time and to link protoclusters with their descendant clusters in observations., Comment: 20 pages, 12 + 10 figures, MNRAS accepted!

A stellar system is dynamically hot when its kinetic energy is dominated by random motion represented by the velocity dispersion $\sigma_{\rm hot} (M_*)$. We use MaNGA data to obtain inner and outer dispersion of a galaxy, $\sigma_{\rm in}$ and $\sigma_{\rm out}$, to characterize its dynamical status and study its connection with star formation quenching and the growth of supermassive black hole (SMBH). We divide galaxies into fully quenched (FQGs), partially quenched (PQGs) and fully star-forming (FSGs) populations, and identify quenched central cores (QCCs) in PQGs. The galaxy distribution in $\sigma_{\rm in}/\sigma_{\rm hot}$-$\sigma_{\rm out}/\sigma_{\rm hot}$ diagram is L-shaped, consisting of a horizontal sequence ($\sigma_{\rm out}/\sigma_{\rm hot}\sim0$) and a vertical sequence ($\sigma_{\rm in}/\sigma_{\rm hot}\sim1$). FQGs and QCCs are located at the top of vertical sequence, $\sigma_{\rm out}/\sigma_{\rm hot}\sim1$, therefore they are dynamically hot over their entire bodies. PQGs reside along vertical sequence, so they have hot center but cold outskirt. FSGs are diverse and can be found in both sequences. Galaxy structural properties, star formation and AGN activities make a transition along horizontal sequence at $\log(\sigma_{\rm in}/\sigma_{\rm hot})\sim-0.3$, and along vertical sequence at $\log(\sigma_{\rm out}/\sigma_{\rm hot})\sim-0.3$. The fractions of optical AGNs and barred galaxies increase rapidly in the first transition and decline rapidly in the second; radio galaxies are located at the top of vertical sequence. Our results demonstrate that star formation quenching and SMBH growth are effective only in dynamically hot systems. A simple model along this line can reproduce the observed SMBH scaling relations. We discuss how secular processes and strong interactions can make a system dynamically hot, and lead to the SMBH growth and star formation quenching., Comment: 24 pages, 19 figures, 1 table, accepted by ApJ

The presence of massive galaxies at high $z$ as recently observed by JWST appears to contradict the current $\Lambda$CDM cosmology. Here we aim to alleviate this tension by incorporating uncertainties from three sources in counting galaxies: cosmic variance, error in stellar mass estimation, and backsplash enhancement. Each of these factors significantly increases the cumulative stellar mass density $\rho_*(>M_*)$ at the high-mass end, and their combined effect can boost the density by more than one order of magnitude. Assuming a star formation efficiency of $\epsilon_* \sim 0.5$, cosmic variance alone reduces the tension to a $2\sigma$ level, except for the most massive galaxy at $z=8$. Additionally, incorporating a 0.3 dex lognormal dispersion in the stellar mass estimation brings the observed $\rho_*(>M_*)$ at $z \sim 7 - 10$ within $2\sigma$. The tension is completely eliminated when we account for the gas stripped from backsplash halos. These results highlight the importance of fully modeling uncertainties when interpreting observational data of rare objects. We use the constrained simulation, ELUCID, to investigate the descendants of high-$z$ massive galaxies. Our findings reveal that a significant portion of these galaxies ultimately reside in massive halos at $z=0$ with $M_{\rm halo} > 10^{13} h^{-1}M_\odot$. Moreover, a large fraction of local central galaxies in $M_{\rm halo} \geqslant 10^{14.5} h^{-1}M_\odot$ halos are predicted to contain substantial amounts of ancient stars formed in massive galaxies at $z \sim 8$. This prediction can be tested by studying the structure and stellar population of central galaxies in present-day massive clusters., Comment: 19 pages, 15 figures. MNRAS accepted

We present an algorithm to extend subhalo merger trees in a low-resolution dark-matter-only simulation by conditionally matching them to those in a high-resolution simulation. The algorithm is general and can be applied to simulation data with different resolutions using different target variables. We instantiate the algorithm by a case in which trees from ELUCID, a constrained simulation of $(500h^{-1}{\rm Mpc})^3$ volume of the local universe, are extended by matching trees from TNGDark, a simulation with much higher resolution. Our tests show that the extended trees are statistically equivalent to the high-resolution trees in the joint distribution of subhalo quantities and in important summary statistics relevant to modeling galaxy formation and evolution in halos. The extended trees preserve certain information of individual systems in the target simulation, including properties of resolved satellite subhalos, and shapes and orientations of their host halos. With the extension, subhalo merger trees in a cosmological scale simulation are extrapolated to a mass resolution comparable to that in a higher-resolution simulation carried out in a smaller volume, which can be used as the input for (sub)halo-based models of galaxy formation. The source code of the algorithm, and halo merger trees extended to a mass resolution of $\sim 2 \times 10^8 h^{-1}M_\odot$ in the entire ELUCID simulation, are available., Comment: 24 pages, 14 figures, 3 tables. MNRAS accepted

Using the IllustrisTNG simulation, we study the interaction of large-scale shocks with the circumgalactic medium (CGM) and interstellar medium (ISM) of star-forming (SF) satellite galaxies in galaxy clusters. These shocks are usually produced by mergers and massive accretion. Our visual inspection shows that approximately half of SF satellites have encountered shocks in their host clusters at $z\leq0.11$. After a satellite crosses a shock front and enters the postshock region, the ram pressure on it is boosted significantly. Both the CGM and ISM can be severely impacted, either by striping or compression. The stripping of the ISM is particularly important for low-mass galaxies with $\log (M_{*}/M_{\odot})<10$ and can occur even in the outskirts of galaxy clusters. In comparison, satellites that do not interact with shocks lose their ISM only in the inner regions of clusters. About half of the ISM is stripped within about 0.6 Gyr after it crosses the shock front. Our results show that shock-induced stripping plays an important role in quenching satellite galaxies in clusters., Comment: 14 pages,9 figures. Accepted for publication in ApJ

We present a new HI mass estimator which relates the HI-to-stellar mass ratio to four galaxy properties: stellar surface mass density, color index $u-r$, stellar mass and concentration index, with the scatter of individual galaxies around the mean HI mass modeled with a Gaussian distribution. We calibrate the estimator using the xGASS sample, including both HI detection and non-detection, and constrain the model parameters through Bayesian inferences. Tests with mock catalogs demonstrate that our estimator provides unbiased HI masses for optical samples like the SDSS, thus suitable for statistical studies of HI gas contents in galaxies and dark matter halos. We apply our estimator to the SDSS spectroscopic sample to estimate the local HI mass function (HIMF), the conditional HI mass function (CHIMF) in galaxy groups and the HI-halo mass (HIHM) relation. Our HIMF agrees with the ALFALFA measurements at $M_{HI}\gtrsim 5\times 10^9M_{\odot}$, but with higher amplitude and a steeper slope at lower masses. We show that this discrepancy is caused primarily by the cosmic variance which is corrected for the SDSS sample but not for the ALFALFA. The CHIMFs for all halo masses can be described by a single Schechter function, and this is true for red, blue and satellite galaxies. For central galaxies the CHIMFs show a double-Gaussian profile, with the two components contributed by the red and blue galaxies, respectively. The total HI mass in a group increases monotonically with halo mass. The HI mass of central galaxies in galaxy groups increases rapidly with halo mass only at $M_h\lesssim10^{12}M_{\odot}$, while the mass dependence becomes much weaker at higher halo masses. The observed HI-halo mass relation is not reproduced by current hydrodynamic simulations and semi-analytic models of galaxy formation., Comment: 10 figures, 2 tables, published in ApJ. $\mathbf{Note}$ : The version published in ApJ has a typo. In the last paragraph of section 3.2, the maximum posterior value of c_a should be c_a = 0.10 \pm 0.08, not c_a = 0.16 \pm 0.10

We use a large N-body simulation to study the characteristic scales in the density gradient profiles in and around halos with masses ranging from $10^{12}$ to $10^{15} h^{-1}{\rm M_\odot}$. We investigate the profiles separately along the major (T_1) and minor (T_3) axes of the local tidal tensor and how the characteristic scales depend on halo mass, formation time, and environment. We find two kinds of prominent characteristic features in the gradient profiles, a deep `valley' and a prominent `peak'. We use the Gaussian Process Regression to fit the gradient profiles and identify the local extrema to determine the scales associate with these features. Around the valley, we identify three types of distinct local minima, corresponding to caustics of particles orbiting around halos. The appearance and depth of the three caustics depend significantly on the direction defined by the local tidal field, formation time and environment of halos. The first caustic is located at a radius r>0.8R_{200}, corresponding to the splashback feature, and is dominated by particles at their first apocenter after infall. The second and third caustics, around 0.6R_{200} and 0.4R_{200} respectively, can be determined reliably only for old halos. The first caustic is always the most prominent feature along T_3, but may not be the case along T_1 or in azimuthally-averaged profiles, suggesting that caution must be taken when using averaged profiles to investigate the splashback radius. We find that the splashback feature is approximately isotropic when proper separations are made between the first and the other caustics. We also identify a peak feature located at $\sim$ 2.5R_{200} in the density gradient profile. This feature is the most prominent along T_1 and is produced by mass accumulations from the structure outside halos. We also discuss the origins of these features and their observational implications., Comment: 15 pages, 12 figures, accepted by A&A

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

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

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

We develop a model to establish the interconnection between galaxies and their dark matter halos. We use Principal Component Analysis (PCA) to reduce the dimensionality of both the mass assembly histories of halos/subhalos and the star formation histories of galaxies, and Gradient Boosted Decision Trees (GBDT) to transform halo/subhalo properties into galaxy properties. We use two sets of hydrodynamic simulations to motivate our model architecture and to train the transformation. We then apply the two sets of trained models to dark matter only (DMO) simulations to show that the transformation is reliable and statistically accurate. The model trained by a high-resolution hydrodynamic simulation, or by a set of such simulations implementing the same physics of galaxy formation, can thus be applied to large DMO simulations to make "mock" copies of the hydrodynamic simulation. The model is both flexible and interpretable, which paves the way for future applications in which we will constrain the model using observations at different redshifts simultaneously and explore how galaxies form and evolve in dark matter halos empirically., Comment: 21 pages, 16 figures, 2 tables, MNRAS accepted

We develop a method to identify proto-clusters based on dark matter halos represented by galaxy groups selected from surveys of galaxies at high redshift. We test the performance of this method on halos in N-body simulations, and find it can correctly identify more than 85% of the true proto-clusters with $\gtrsim$ 95% purity and with mass estimates typical within 0.25 dex from their true values. We show how the information provided by the proto-clusters can be used to link galaxies in present-day clusters of galaxies with their high redshift progenitors. Our tests show that the proto-clusters identified by our method can recover reliably the progenitor stellar mass distribution of galaxies, thereby providing an avenue to investigate the formation and evolution of present-day galaxy clusters and their member galaxies., Comment: 15 pages, 10 + 4 figures, MNRAS accepted, Comments are welcome

Using two sets of large $N$-body simulations, we study the origin of the correlations of halo assembly time ($z_{\rm f}$), concentration ($v_{\rm max}/v_{\rm 200}$) and spin ($\lambda$) with the large-scale evolved density field at given halo mass, i.e. the secondary bias. We find that the secondary bias is the secondary effect of the correlations of halo properties with the linear density estimated at the same comoving scale. Using the linear density on different scales, we find two types of correlations. The internal correlation, which reflects the correlation of halo properties with the mean linear over-density $\delta_{\rm L}$ within the halo Lagrangian radius $R_{\rm L}$, is positive for both $z_{\rm f}$ and $v_{\rm max}/v_{\rm 200}$, and negative for $\lambda$. The external correlation, which describes the correlation of halo properties with linear overdensity at $R>R_{\rm L}$ for given $\delta_{\rm L}$, shows trends opposite to the internal correlation. Both of the external and internal correlations depend only weakly on halo mass, indicating a similar origin for halos of different masses. Our findings offer a transparent perspective on the origin of the secondary bias. The secondary bias can be largely explained by the competition of the external and internal correlations together with the correlation of the linear density field on different scales. These two types of correlations combined can establish the complex halo-mass dependence of the secondary bias observed in the simulations., Comment: 14 pages,14 figures, accepted by A&A

Based on the spectroscopic and shear catalogs for SDSS galaxies in the local Universe, we compare optically-selected active galactic nuclei (AGNs) with control star-forming and quiescent galaxies on galactic, inter-halo and larger scales. We find that AGNs are preferentially found in two specific stages of galaxy evolution: star-burst and `green valley' phases, and that the stellar population of their host galaxies is quite independent of stellar mass, different from normal galaxies. Combining galaxy-galaxy lensing and galaxy clustering on large scales, we measure the mass of AGN host halos. The typical halo mass is about $10^{12}h^{-1}\rm M_{\odot}$, similar to the characteristic mass in the stellar mass-halo mass relation (SHMR). For given stellar mass, AGN host galaxies and star-forming galaxies share the same SHMR, while quiescent galaxies have more massive halos. Clustering analysis on halo scales reveals that AGNs are surrounded by a larger number of satellites (with stellar mass down to 1/1000 of the mass of the central galaxy) than star-forming galaxies, and that galaxies with larger stellar velocity dispersion have more satellites. The number of satellites also increase with halo mass, reaching unity around $10^{12}h^{-1}\rm M_{\odot}$. Our results suggest a scenario, in which the interaction of the central galaxy with the satellites triggers an early episode of star burst and AGN activities, followed by multiple AGN cycles driven by the non-axisymmetric structure produced by the interaction. The feedback from the starburst and AGN reduces the amount of cold gas for fueling the central black hole, producing a characteristic halo mass scale, $\sim 10^{12}h^{-1}\rm M_{\odot}$, where the AGN fraction peaks., Comment: accepted by A&A, 17 pages, 11 figures

We use TNG and EAGLE hydrodynamic simulations to investigate the central galaxy - dark matter halo relations that are needed for a halo-based empirical model of star formation in galaxies. Using a linear dimension reduction algorithm and a model ensemble method, we find that for both star-forming and quenched galaxies, the star formation history (SFH) is tightly related to the halo mass assembly history (MAH). The quenching of a low-mass galaxy is mainly due to the infall-ejection process related to a nearby massive halo, while the quenching of a high-mass galaxy is closely related to the formation of a massive progenitor in its host halo. The classification of star-forming and quenched populations based solely on halo properties contains contamination produced by sample imbalance and overlapping distributions of the two populations. Guided by the results from hydrodynamic simulations, we build an empirical model to predict the SFH of central galaxies based on the MAH of their host halos, and we model the star-forming and quenched populations separately. Our model is based on the idea of adopting star formation templates from hydrodynamic simulations to reduce model complexity. We use various tests to demonstrate that the model can recover star formation histories of individual galaxies, and can statistically reproduce the galaxy bimodal distribution, stellar mass - halo mass and star formation rate - halo mass relations from low to high redshift, and assembly bias. Our study provides a framework of using hydrodynamic simulations to discover, and to motivate the use of, key ingredients to model galaxy formation using halo properties., Comment: 20 pages, 17 figures, 2 tables. MNRAS accepted. Comments are welcomed

We present a comparison of the physical properties of the ionized gas in the circumgalactic (CGM) and intergalactic (IGM) media at $z\sim0$ between observations and four cosmological hydrodynamical simulations: Illustris, TNG300 of the IllustrisTNG project, EAGLE, and one of the Magneticum simulations. For the observational data, we use the gas properties that are inferred from cross-correlating the Sunyaev-Zel'dovich effect (SZE) from the {\it Planck} CMB maps with the haloes and the large-scale structure reconstructed from Sloan Digital Sky Survey data. Both the observational and simulation results indicate that the integrated gas pressure in haloes deviates from the self-similar case, showing that feedback impacts haloes with $M_{500}\sim 10^{12-13}\,{\rm M_\odot}$. The simulations predict that more than half the baryons are displaced from haloes, while the gas fraction inferred from our observational data roughly equals the cosmic baryon fraction throughout the $M_{500}\sim 10^{12-14.5}\,{\rm M_\odot}$ halo mass range. All simulations tested here predict that the mean gas temperature in haloes is about the virial temperature, while that inferred from the SZE is up to one order of magnitude lower than that from the simulations (and also from X-ray observations). While a remarkable agreement is found for the average properties of the IGM between the observation and some simulations, we show that their dependence on the large-scale tidal field can break the degeneracy between models that show similar predictions otherwise. Finally, we show that the gas pressure and the electron density profiles from simulations are not well described by a generalized NFW (GNFW) profile. Instead, we present a new model with a mass-dependent shape that fits the profiles accurately., Comment: 15 pages, 6 figures, accepted for publication in MNRAS

Identifying galaxy groups from redshift surveys of galaxies plays an important role in connecting galaxies with the underlying dark matter distribution. Current and future high-$z$ spectroscopic surveys, usually incomplete in redshift sampling, present both opportunities and challenges to identifying groups in the high-$z$ Universe. We develop a group finder that is based on incomplete redshift samples combined with photometric data, using a machine learning method to assign halo masses to identified groups. Test using realistic mock catalogs shows that $\gtrsim 90\%$ of true groups with halo masses $\rm M_h \gtrsim 10^{12} M_{\odot}/h$ are successfully identified, and that the fraction of contaminants is smaller than $10\%$. The standard deviation in the halo mass estimation is smaller than 0.25 dex at all masses. We apply our group finder to zCOSMOS-bright and describe basic properties of the group catalog obtained., Comment: 18 pages, 12 figures, 6 table. MNRAS accepted. Comments are welcomed

A clear transition feature of galaxy quenching is identified in the multi-parameter space of stellar mass ($M_*$), bulge to total mass ratio ($B/T_{\rm m}$), halo mass ($M_{\rm h}$) and halo-centric distance ($r/r_{180}$). For given halo mass, the characteristic stellar mass ($M_{*, \rm ch}$) for the transition is about one-fifth of that of the corresponding central galaxy, and almost independent of $B/T_{\rm m}$. Once $B/T_{\rm m}$ is fixed, the quenched fraction of galaxies with $M_{*} < M_{*, \rm ch}$ increases with $M_{\rm h}$, but decreases with $M_*$ in the inner part of halos ($r/r_{180} < 0.5$). In the outer part ($r/r_{180} > 0.5$), the trend with $M_{\rm h}$ remains but the correlation with $M_*$ is absent or becomes positive. For galaxies above $M_{\rm *, ch}$ and with $B/T_{\rm m}$ fixed, the quenched fraction increases with $M_{\rm *}$, but depends only weakly on $M_{\rm h}$ in both the inner and outer regions. At fixed $B/T_{\rm m}$ and $M_*$, the quenched fraction increases with decreasing $r/r_{180}$ for galaxies with $M_{*} < M_{*, \rm ch}$, and depends only weakly on $r/r_{180}$ for galaxies with $M_{*} > M_{*, \rm ch}$. Our finding provides a physically-motivated way to classify galaxies in halos into two classes based on their quenching properties: an `upper class' with $M_{*} > M_{\rm *,ch}$ and a `lower class' with $M_{*} < M_{\rm *,ch}$. Environmental quenching is important for `lower class' galaxies, while internal quenching plays the dominating role for the `upper class'., Comment: 11 pages, 7 figures, accepted by ApJ

We measure the star formation histories (SFH) of a sample of low-mass galaxies with $M_\ast<10^9M_\odot$ from the SDSS-IV MaNGA survey. The large number of IFU spectra for each galaxy are either combined to reach a high signal to noise ratio or used to investigate spatial variations. We use Bayesian inferences based on full spectrum fitting. Our analysis based on Bayesian evidence ratio indicates a strong preferencefor a model that allows the presence of an old stellar population, and that an improper model for the SFH can significantly underestimate the old population in these galaxies. The addition of NIR photometry to the constraining data can further distinguish between different SFH model families and significantly tighten the constraints on the mass fraction in the old population. On average more than half of the stellar mass in present-day low-mass galaxies formed at least 8 Gyrs ago, while about 30\% within the past 4 Gyrs. Satellite galaxies on average have formed their stellar mass earlier than central galaxies. The radial dependence of the SFH is quite weak. Our results suggest that most of the low-mass galaxies have an early episode of active star formation that produces a large fraction of their present stellar mass., Comment: 21 pages, 22 figures, MNRAS accepted version

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

We present the detection of the kinetic Sunyaev-Zel'dovich effect (kSZE) signals from groups of galaxies as a function of halo mass down to $\log (M_{500}/{\rm M_\odot}) \sim 12.3$, using the {\it Planck} CMB maps and stacking about $40,000$ galaxy systems with known positions, halo masses, and peculiar velocities. The signals from groups of different mass are constrained simultaneously to take care of projection effects of nearby halos. The total kSZE flux within halos estimated implies that the gas fraction in halos is about the universal baryon fraction, even in low-mass halos, indicating that the `missing baryons' are found. Various tests performed show that our results are robust against systematic effects, such as contamination by infrared/radio sources and background variations, beam-size effects and contributions from halo exteriors. Combined with the thermal Sunyaev-Zel'dovich effect, our results indicate that the `missing baryons' associated with galaxy groups are contained in warm-hot media with temperatures between $10^5$ and $10^6\,{\rm K}$., Comment: We replace the submission of the old version of the manuscript (arXiv:1712.08619; same paper, basically) with the manuscript of this new version, and remove this new submission to avoid confusion, which we newly submitted only to make it noticed that there was a very significant revision to the old version although they are still not a separate paper, but same paper

We analyze the stellar initial mass functions (IMF) of a large sample of early type galaxies (ETGs) provided by MaNGA. The large number of IFU spectra of individual galaxies provide high signal-to-noise composite spectra that are essential for constraining IMF and to investigate possible radial gradients of the IMF within individual galaxies. The large sample of ETGs also make it possible to study how the IMF shape depends on various properties of galaxies. We adopt a novel approach to study IMF variations in ETGs, use Bayesian inferences based on full spectrum fitting. The Bayesian method provides a statistically rigorous way to explore potential degeneracy in spectrum fitting, and to distinguish different IMF models with Bayesian evidence. We find that the IMF slope depends systematically on galaxy velocity dispersion, in that galaxies of higher velocity dispersion prefer a more bottom-heavy IMF, but the dependence is almost entirely due to the change of metallicity, $Z$, with velocity dispersion. The IMF shape also depends on stellar age, $A$, but the dependence is completely degenerate with that on metallicity through a combination $AZ^{-1.42}$. Using independent age and metallicity estimates we find that the IMF variation is produced by metallicity instead of age. The IMF near the centers of massive ETGs appears more bottom-heavy than that in the outer parts, while a weak opposite trend is seen for low-mass ETGs. Uncertainties produced by star formation history, dust extinction, $\alpha$-element abundance enhancement and noise in the spectra are tested., Comment: 21 pages,20 figures, accepted for publication in MNRAS

We perform full spectrum fitting stellar population analysis and Jeans Anisotropic modelling (JAM) of the stellar kinematics for about 2000 early-type galaxies (ETGs) and spiral galaxies from the MaNGA DR14 sample. Galaxies with different morphologies are found to be located on a remarkably tight mass plane which is close to the prediction of the virial theorem, extending previous results for ETGs. By examining an inclined projection (`the mass-size' plane), we find that spiral and early-type galaxies occupy different regions on the plane, and their stellar population properties (i.e. age, metallicity and stellar mass-to-light ratio) vary systematically along roughly the direction of velocity dispersion, which is a proxy for the bulge fraction. Galaxies with higher velocity dispersions have typically older ages, larger stellar mass-to-light ratios and are more metal rich, which indicates that galaxies increase their bulge fractions as their stellar populations age and become enriched chemically. The age and stellar mass-to-light ratio gradients for low-mass galaxies in our sample tend to be positive ($\rm centre

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

We extend the real-space mapping method developed in Shi et at. (2016) so that it can be applied to flux-limited galaxy samples. We use an ensemble of mock catalogs to demonstrate the reliability of this extension, showing that it allows for an accurate recovery of the real-space correlation functions and galaxy biases. We also demonstrate that, using an iterative method applied to intermediate-scale clustering data, we can obtain an unbiased estimate of the growth rate of structure $f\sigma_8$, which is related to the clustering amplitude of matter, to an accuracy of $\sim 10\%$. Applying this method to the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7), we construct a real-space galaxy catalog spanning the redshift range $0.01 \leq z \leq 0.2$, which contains 584,473 galaxies in the north Galactic cap (NGC). Using this data, we infer $\fss$ at a median redshift $z=0.1$, which is consistent with the WMAP9 cosmology at the $1\sigma$ level. By combining this measurement with the real-space clustering of galaxies and with galaxy-galaxy weak lensing measurements for the same sets of galaxies, we are able to break the degeneracy between $f$, $\sigma_8$, and $b$. From the SDSS DR7 data alone, we obtain the following cosmological constraints at redshift $z=0.1$: $f=$$0.464^{+0.040}_{-0.040}$, $\sigma_8=0.769^{+0.121}_{-0.089}$, and $b=1.910^{+0.234}_{-0.268}$, $1.449^{+0.194}_{-0.196}$, $1.301^{+0.170}_{-0.177}$, and $1.196^{+0.159}_{-0.161}~$ for galaxies within different absolute magnitude bins $^{0.1}{\rm M}_r-5\log h=[-23,0, -22.0], [-22,0, -21.0], [-21.0, -20.0]$ and $[-20.0, -19.0]$, respectively.

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., Comment: 18 pages, 13 figures, ApJ accepted

We use a 200 $h^{-1}Mpc$ a side N-body simulation to study the mass accretion history (MAH) of dark matter halos to be accreted by larger halos, which we call infall halos. We define a quantity $a_{\rm nf}\equiv (1+z_{\rm f})/(1+z_{\rm peak})$ to characterize the MAH of infall halos, where $z_{\rm peak}$ and $z_{\rm f}$ are the accretion and formation redshifts, respectively. We find that, at given $z_{\rm peak}$, their MAH is bimodal. Infall halos are dominated by a young population at high redshift and by an old population at low redshift. For the young population, the $a_{\rm nf}$ distribution is narrow and peaks at about $1.2$, independent of $z_{\rm peak}$, while for the old population, the peak position and width of the $a_{\rm nf}$ distribution both increases with decreasing $z_{\rm peak}$ and are both larger than those of the young population. This bimodal distribution is found to be closely connected to the two phases in the MAHs of halos. While members of the young population are still in the fast accretion phase at $z_{\rm peak}$, those of the old population have already entered the slow accretion phase at $z_{\rm peak}$. This bimodal distribution is not found for the whole halo population, nor is it seen in halo merger trees generated with the extended Press-Schechter formalism. The infall halo population at $z_{\rm peak}$ are, on average, younger than the whole halo population of similar masses identified at the same redshift. We discuss the implications of our findings in connection to the bimodal color distribution of observed galaxies and to the link between central and satellite galaxies., Comment: 21 pages, 7 figures, accepted for publication by ApJ

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

We study the dependence of galaxy clustering on atomic gas mass using a sample of $\sim$16,000 galaxies with redshift in the range of $0.002510^8M_{\odot}$, drawn from the 70% complete sample of the Arecibo Legacy Fast ALFA survey. We construct subsamples of galaxies with $M_{\rm HI}$ above different thresholds, and make volume-limited clustering measurements in terms of three statistics: the projected two-point correlation function, the projected cross-correlation function with respect to a reference sample selected from the Sloan Digital Sky Survey, and the redshift-space monopole moment. In contrast to previous studies, which found no/weak HI-mass dependence, we find both the clustering amplitude on scales above a few Mpc and the bias factors to increase significantly with increasing HI mass for subsamples with HI mass thresholds above $10^9M_{\odot}$. For HI mass thresholds below $10^9M_{\odot}$, while the measurements have large uncertainties caused by the limited survey volume and sample size, the inferred galaxy bias factors are systematically lower than the minimum halo bias factor from mass-selected halo samples. The simple halo model, in which galaxy content is only determined by halo mass, has difficulties in interpreting the clustering measurements of the HI-selected samples. We extend the simple model by including the halo formation time as an additional parameter. A model that puts HI-rich galaxies into halos that formed late can reproduce the clustering measurements reasonably well. We present the implications of our best-fitting model on the correlation of HI mass with halo mass and formation time, as well as the halo occupation distributions and HI mass functions for central and satellite galaxies. These results are compared with the predictions from semi-analytic galaxy formation models and hydrodynamic galaxy formation simulations., Comment: Accepted for publication in ApJ. The 2PCF measurements are available at http://sdss4.shao.ac.cn/guoh/

The ELUCID project aims to build a series of realistic cosmological simulations that reproduce the spatial and mass distribution of the galaxies as observed in the Sloan Digital Sky Survey (SDSS). This requires powerful reconstruction techniques to create constrained initial conditions. We test the reconstruction method by applying it to several $N$-body simulations. We use 2 medium resolution simulations from each of which three additional constrained $N$-body simulations were produced. We compare the resulting friend of friend catalogs by using the particle indexes as tracers, and quantify the quality of the reconstruction by varying the main smoothing parameter. The cross identification method we use proves to be efficient, and the results suggest that the most massive reconstructed halos are effectively traced from the same Lagrangian regions in the initial conditions. Preliminary time dependence analysis indicates that high mass end halos converge only at a redshift close to the reconstruction redshift. This suggests that, for earlier snapshots, only collections of progenitors may be effectively cross-identified.

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

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