Showing 1–50 of 240 results for author: Mo, H

From Halos to Galaxies. VI. Improved halo mass estimation for SDSS groups and measurement of the halo mass function

Authors: Dingyi Zhao, Yingjie Peng, Yipeng Jing, Xiaohu Yang, Luis C. Ho, Alvio Renzini, Anna R. Gallazzi, Cheqiu Lyu, Roberto Maiolino, Jing Dou, Zeyu Gao, Qiusheng Gu, Filippo Mannucci, Houjun Mo, Bitao Wang, Enci Wang, Kai Wang, Yu-Chen Wang, Bingxiao Xu, Feng Yuan, Xingye Zhu

Abstract: In $Λ$CDM cosmology, galaxies form and evolve in their host dark matter (DM) halos. Halo mass is crucial for understanding the halo-galaxy connection. The abundance matching (AM) technique has been widely used to derive the halo masses of galaxy groups. However, quenching of the central galaxy can decouple the coevolution of its stellar mass and DM halo mass. Different halo assembly histories can… ▽ More In $Λ$CDM cosmology, galaxies form and evolve in their host dark matter (DM) halos. Halo mass is crucial for understanding the halo-galaxy connection. The abundance matching (AM) technique has been widely used to derive the halo masses of galaxy groups. However, quenching of the central galaxy can decouple the coevolution of its stellar mass and DM halo mass. Different halo assembly histories can also result in significantly different final stellar mass of the central galaxies. These processes can introduce substantial uncertainties in the halo masses derived from the AM method, particularly leading to a systematic bias between groups with star-forming centrals (blue groups) and passive centrals (red groups). To improve, we developed a new machine learning (ML) algorithm that accounts for these effects and is trained on simulations. Our results show that the ML method eliminates the systematic bias in the derived halo masses for blue and red groups and is, on average, $\sim1/3$ more accurate than the AM method. With careful calibration of observable quantities from simulations and observations from SDSS, we apply our ML model to the SDSS Yang et al. groups to derive their halo masses down to $10^{11.5}\mathrm{M_\odot}$ or even lower. The derived SDSS group halo mass function agrees well with the theoretical predictions, and the derived stellar-to-halo mass relations for both red and blue groups matches well with those obtained from direct weak lensing measurements. These new halo mass estimates enable more accurate investigation of the galaxy-halo connection and the role of the halos in galaxy evolution. △ Less

Submitted 22 August, 2024; originally announced August 2024.

Comments: 21 pages, 9 figures. Submitted to ApJ, comments welcome!

On the kinematic nature of apparent discs at high redshifts: Local counterparts are not dominated by ordered rotation but by tangentially anisotropic random motion

Authors: Bitao Wang, Yingjie Peng, Michele Cappellari, Hua Gao, Houjun Mo

Abstract: It is not straightforward to physically interpret the apparent morphology of galaxies. Recent observations by James Webb Space Telescope (JWST) revealed a dominant galaxy population at high redshifts ($z>2$) that were visually classified as discs for their flattened shapes and/or exponential light profiles. The extensively accepted interpretation is that they are dynamically cold discs supported b… ▽ More It is not straightforward to physically interpret the apparent morphology of galaxies. Recent observations by James Webb Space Telescope (JWST) revealed a dominant galaxy population at high redshifts ($z>2$) that were visually classified as discs for their flattened shapes and/or exponential light profiles. The extensively accepted interpretation is that they are dynamically cold discs supported by bulk rotation. However, it is long known that flattened shapes and exponential profiles are not exclusive for rotating disc structure. To break degeneracy and assess the rotational support of typical high-$z$ galaxies in the JWST samples, those with active star formation and stellar masses $\mathrm{lg}(\mathcal{M}_{\star}/\mathcal{M}_{\odot})\sim9$, we study the kinematics of their equal-mass counterparts at $z=0$. While these local star-forming low-mass galaxies are photometrically similar to real dynamically cold discs, they are not supported by ordered rotation but primarily by random motion, and their flattened shapes result largely from tangential orbital anisotropy. Given the empirical and theoretical evidence that young galaxies are dynamically hotter at higher redshifts, our results suggest that the high-$z$ JWST galaxies may not be cold discs but are dynamically warm/hot galaxies with flattened shapes driven by anisotropy. While both having low rotational support, local low-mass galaxies possess oblate shapes, contrasting the prolate shapes (i.e. cigar-like) of low-mass systems at high redshifts. Such shape transition (prolate$\Rightarrow$oblate) indicates an associated change in orbital anisotropy (radial$\Rightarrow$tangential), with roots likely in the assembly of their host dark matter halos. △ Less

Submitted 20 August, 2024; originally announced August 2024.

Comments: 11 pages, 4 figures. Submitted to ApJL and revised according to referee's report

From Halos to Galaxies. X: Decoding Galaxy SEDs with Physical Priors and Accurate Star Formation History Reconstruction

Authors: Zeyu Gao, Yingjie Peng, Kai Wang, Luis C. Ho, Alvio Renzini, Anna R. Gallazzi, Filippo Mannucci, Houjun Mo, Yipeng Jing, Xiaohu Yang, Enci Wang, Dingyi Zhao, Jing Dou, Qiusheng Gu, Cheqiu Lyu, Roberto Maiolino, Bitao Wang, Yu-Chen Wang, Bingxiao Xu, Feng Yuan, Xingye Zhu

Abstract: The spectral energy distribution (SED) of galaxies is essential for deriving fundamental properties like stellar mass and star formation history (SFH). However, conventional methods, including both parametric and non-parametric approaches, often fail to accurately recover the observed cosmic star formation rate (SFR) density due to oversimplified or unrealistic assumptions about SFH and their inab… ▽ More The spectral energy distribution (SED) of galaxies is essential for deriving fundamental properties like stellar mass and star formation history (SFH). However, conventional methods, including both parametric and non-parametric approaches, often fail to accurately recover the observed cosmic star formation rate (SFR) density due to oversimplified or unrealistic assumptions about SFH and their inability to account for the complex SFH variations across different galaxy populations. To address this issue, we introduce a novel approach that improves galaxy broad-band SED analysis by incorporating physical priors derived from hydrodynamical simulations. Tests using IllustrisTNG simulations demonstrate that our method can reliably determine galaxy physical properties from broad-band photometry, including stellar mass within 0.05 dex, current SFR within 0.3 dex, and fractional stellar formation time within 0.2 dex, with a negligible fraction of catastrophic failures. When applied to the SDSS main photometric galaxy sample with spectroscopic redshift, our estimates of stellar mass and SFR are consistent with the widely-used MPA-JHU and GSWLC catalogs. Notably, using the derived SFHs of individual SDSS galaxies, we estimate the cosmic SFR density and stellar mass density with remarkable consistency to direct observations up to $z \sim 6$. This marks the first time SFHs derived from SEDs can accurately match observations. Consequently, our method can reliably recover observed spectral indices such as $\rm D_{\rm n}(4000)$ and $\rm Hδ_{\rm A}$ by synthesizing the full spectra of galaxies using the estimated SFHs and metal enrichment histories, relying solely on broad-band photometry as input. Furthermore, this method is extremely computationally efficient compared to conventional approaches. △ Less

Submitted 14 August, 2024; originally announced August 2024.

Comments: 22 pages, 12 figures. Submitted to ApJ, comments welcome!

From Halos to Galaxies. IX. Estimate of Halo Assembly History for SDSS Galaxy Groups

Authors: Cheqiu Lyu, Yingjie Peng, Yipeng Jing, Xiaohu Yang, Luis C. Ho, Alvio Renzini, Dingyi Zhao, Filippo Mannucci, Houjun Mo, Kai Wang, Bitao Wang, Bingxiao Xu, Jing Dou, Anna R. Gallazzi, Qiusheng Gu, Roberto Maiolino, Enci Wang, Feng Yuan

Abstract: The properties of the galaxies are tightly connected to their host halo mass and halo assembly history. Accurate measurement of the halo assembly history in observation is challenging but crucial to the understanding of galaxy formation and evolution. The stellar-to-halo mass ratio ($M_*/M_{\mathrm{h}}$) for the centrals has often been used to indicate the halo assembly time $t_{\mathrm{h,50}}$ of… ▽ More The properties of the galaxies are tightly connected to their host halo mass and halo assembly history. Accurate measurement of the halo assembly history in observation is challenging but crucial to the understanding of galaxy formation and evolution. The stellar-to-halo mass ratio ($M_*/M_{\mathrm{h}}$) for the centrals has often been used to indicate the halo assembly time $t_{\mathrm{h,50}}$ of the group, where $t_{\mathrm{h,50}}$ is the lookback time at which a halo has assembled half of its present-day virial mass. Using mock data from the semi-analytic models, we find that $M_*/M_{\mathrm{h}}$ shows a significant scatter with $t_{\mathrm{h,50}}$, with a strong systematic difference between the group with a star-forming central (blue group) and passive central (red group). To improve the accuracy, we develop machine-learning models to estimate $t_{\mathrm{h,50}}$ for galaxy groups using only observable quantities in the mocks. Since star-formation quenching will decouple the co-growth of the dark matter and baryon, we train our models separately for blue and red groups. Our models have successfully recovered $t_{\mathrm{h,50}}$, within an accuracy of $\sim$ 1.09 Gyr. With careful calibrations of individual observable quantities in the mocks with SDSS observations, we apply the trained models to the SDSS Yang et al. groups and derive the $t_{\mathrm{h,50}}$ for each group for the first time. The derived SDSS $t_{\mathrm{h,50}}$ distributions are in good agreement with that in the mocks, in particular for blue groups. The derived halo assembly history, together with the halo mass, make an important step forward in studying the halo-galaxy connections in observation. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Comments: 18 pages, 7 figures. Accepted by ApJ

A two-phase model of galaxy formation: III. The formation of globular clusters

Authors: Yangyao Chen, Houjun Mo, Huiyuan Wang

Abstract: We develop a model of globular cluster (GC) formation within the cosmological hierarchy of structure formation. The model is rooted in the `two-phase' scenario of galaxy formation developed in Paper-I, where the fast accretion of dark matter halos at high redshift leads to the formation of self-gravitating, turbulent gas clouds that subsequently fragment into dynamically hot systems of dense sub-c… ▽ More We develop a model of globular cluster (GC) formation within the cosmological hierarchy of structure formation. The model is rooted in the `two-phase' scenario of galaxy formation developed in Paper-I, where the fast accretion of dark matter halos at high redshift leads to the formation of self-gravitating, turbulent gas clouds that subsequently fragment into dynamically hot systems of dense sub-clouds with masses $\sim 10^6$-$10^7 M_\odot$. Here we elaborate on the formation, evolution, and fate of these sub-clouds, and show that some of the sub-clouds can be compactified via two distinctive channels into a 'supernova-free' regime to form two distinct populations of GCs. The model is simple, characterized by a small number of free parameters underpinned by physical considerations, and can be efficiently implemented into cosmological N-body simulations to generate a coherent sample of halos, galaxies, and GCs. Our model can reproduce a range of observations on GCs, including the mass function, the size-mass relation, the frequency per unit host galaxy/halo mass, the bimodal metallicity distribution, and the spatial profile. Predictions for GCs are made for both the local Universe and for redshift up to $z \approx 10$, and can be tested by upcoming observations. △ Less

Submitted 28 May, 2024; originally announced May 2024.

Comments: 35 pages, 19 figures, 1 table. Submitted to MNRAS

Gas-rich Ultra-diffuse Galaxies Are Originated from High Specific Angular Momentum

Authors: Yu Rong, Huijie Hu, Min He, Wei Du, Qi Guo, Hui-Yuan Wang, Hong-Xin Zhang, Houjun Mo

Abstract: Ultra-diffuse galaxies, characterized by comparable effective radii to the Milky Way but possessing 100-1,000 times fewer stars, offer a unique opportunity to garner novel insights into the mechanisms governing galaxy formation. Nevertheless, the existing corpus of observational and simulation studies has not yet yielded a definitive constraint or comprehensive consensus on the formation mechanism… ▽ More Ultra-diffuse galaxies, characterized by comparable effective radii to the Milky Way but possessing 100-1,000 times fewer stars, offer a unique opportunity to garner novel insights into the mechanisms governing galaxy formation. Nevertheless, the existing corpus of observational and simulation studies has not yet yielded a definitive constraint or comprehensive consensus on the formation mechanisms underlying ultra-diffuse galaxies. In this study, we delve into the properties of ultra-diffuse galaxies enriched with neutral hydrogen using a semi-analytic method, with the explicit aim of constraining existing ultra-diffuse galaxy formation models. We find that the gas-rich ultra-diffuse galaxies are statistically not failed $L^{\star}$ galaxies nor dark matter deficient galaxies. In statistical terms, these ultra-diffuse galaxies exhibit comparable halo concentration, but higher baryonic mass fraction, as well as higher stellar and gas specific angular momentum, in comparison to typical dwarf galaxy counterparts. Our analysis unveils that higher gas specific angular momentum serves as the underlying factor elucidating the observed heightened baryonic mass fractions, diminished star formation efficiency, expanded stellar disk sizes, and reduced stellar densities in ultra-diffuse galaxies. Our findings make significant contributions to advancing our knowledge of ultra-diffuse galaxy formation and shed light on the intricate interplay between gas dynamics and the evolution of galaxies. △ Less

Submitted 31 March, 2024; originally announced April 2024.

Comments: comments welcome

Connection between galaxy morphology and dark-matter halo structure I: a running threshold for thin discs and size predictors from the dark sector

Authors: Jinning Liang, Fangzhou Jiang, Houjun Mo, Andrew Benson, Avishai Dekel, Noa Tavron, Philip F. Hopkins, Luis C. Ho

Abstract: We present a series of studies on the connection between galaxy morphology and the structure of host dark-matter (DM) haloes using cosmological simulations. In this work, we introduce a new kinematic decomposition scheme that features physical identification of morphological components, enabling robust separation of thin and thick discs; and measure a wide range of halo properties, including their… ▽ More We present a series of studies on the connection between galaxy morphology and the structure of host dark-matter (DM) haloes using cosmological simulations. In this work, we introduce a new kinematic decomposition scheme that features physical identification of morphological components, enabling robust separation of thin and thick discs; and measure a wide range of halo properties, including their locations in the cosmic web, internal structures, and assembly histories. Our analysis of the TNG50 simulation reveals that the orbital-circularity threshold for disc differentiation varies across galaxies, with systematic trends in mass and redshift, so the widely used decomposition method with constant circularity cuts is oversimplified and underestimates thin disc at JWST redshifts. The energy threshold between the stellar halo and the inner galaxy is also a function of mass and redshift, minimizing at the sub-Galactic halo mass, where the circularity threshold peaks. Revisiting the issue of galaxy size predictor, we show that disc sizes in TNG50 exhibit correlations with three structural parameters besides virial mass and redshift: 1) a positive correlation with halo spin $λ$ across redshifts -- stronger than previously reported for zoom-in simulations but still weaker than the simple $r_{1/2}/R_{\rm vir} \propto λ$ scaling; 2) an anti-correlation with DM concentration $c$ that is well described by $r_{1/2}/R_{\rm vir} \propto c^{-0.7}$ even when $c$ is measured in the DM only run; 3) more actively accreting haloes having slightly larger discs, as well as more significant stellar haloes and lower thin-to-thick ratio. Disc mass fraction is higher in rounder haloes and in cosmic knots and filaments, implying that disc development needs both stable halo conditions and continuous material supply. Our methodology is public and adaptable to other simulations. △ Less

Submitted 21 March, 2024; originally announced March 2024.

Comments: 20 pages, 17 figures, submitted to MNRAS

CSST large-scale structure analysis pipeline: I. constructing reference mock galaxy redshift surveys

Authors: Yizhou Gu, Xiaohu Yang, Jiaxin Han, Yirong Wang, Qingyang Li, Zhenlin Tan, Wenkang Jiang, Yaru Wang, Jiaqi Wang, Antonios Katsianis, Xiaoju Xu, Haojie Xu, Wensheng Hong, Houjun Mo, Run Wen, Xianzhong Zheng, Feng Shi, Pengjie Zhang, Zhongxu Zhai, Chengze Liu, Wenting Wang, Ying Zu, Hong Guo, Youcai Zhang, Yi Lu , et al. (7 additional authors not shown)

Abstract: In this paper, we set out to construct a set of reference mock galaxy redshift surveys (MGRSs) for the future Chinese Space-station Survey Telescope (CSST) observation, where subsequent survey selection effects can be added and evaluated. This set of MGRSs is generated using the dark matter subhalos extracted from a high-resolution Jiutian $N$-body simulation of the standard $Λ$CDM cosmogony with… ▽ More In this paper, we set out to construct a set of reference mock galaxy redshift surveys (MGRSs) for the future Chinese Space-station Survey Telescope (CSST) observation, where subsequent survey selection effects can be added and evaluated. This set of MGRSs is generated using the dark matter subhalos extracted from a high-resolution Jiutian $N$-body simulation of the standard $Λ$CDM cosmogony with $Ω_m=0.3111$, $Ω_Λ=0.6889$, and $σ_8=0.8102$. The simulation has a boxsize of $1~h^{-1} {\rm Gpc}$, and consists of $6144^3$ particles with mass resolution $3.723 \times 10^{8} h^{-1} M_\odot $. In order to take into account the effect of redshift evolution, we first use all 128 snapshots in the Jiutian simulation to generate a light-cone halo/subhalo catalog. Next, galaxy luminosities are assigned to the main and subhalo populations using the subhalo abundance matching (SHAM) method with the DESI $z$-band luminosity functions at different redshifts. Multi-band photometries, as well as images, are then assigned to each mock galaxy using a 3-dimensional parameter space nearest neighbor sampling of the DESI LS observational galaxies and groups. Finally, the CSST and DESI LS survey geometry and magnitude limit cuts are applied to generate the required MGRSs. As we have checked, this set of MGRSs can generally reproduce the observed galaxy luminosity/mass functions within 0.1 dex for galaxies with $L > 10^8 L_\odot$ (or $M_* > 10^{8.5} M_\odot$) and within 1-$σ$ level for galaxies with $L < 10^8L_\odot$ (or $M_* < 10^{8.5} M_\odot$). Together with the CSST slitless spectra and redshifts for our DESI LS seed galaxies that are under construction, we will set out to test various slitless observational selection effects in subsequent probes. △ Less

Submitted 15 March, 2024; originally announced March 2024.

Comments: 13 pages, 9 figures, accepted for publication in MNRAS

Properties of a Fading AGN from SDSS-IV MaNGA

Authors: Hao Mo, Yan-Mei Chen, Zhi-Yun Zhang, Alexei Moiseev, Dmitry Bizyaev, Yong Shi, Qiu-Sheng Gu, Min Bao, Xiao Cao, Song-Lin Li

Abstract: We identify a fading AGN SDSS J220141.64+115124.3 from the internal Product Launch-11 (MPL-11) in Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. The central region with a projected radius of $\sim$2.4 kpc is characterized as LINER-like line ratios while the outskirts extended to $\sim$15 kpc show Seyfert-like line ratios. The [OIII]$λ$5007 luminosity of the Seyfert regions is… ▽ More We identify a fading AGN SDSS J220141.64+115124.3 from the internal Product Launch-11 (MPL-11) in Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. The central region with a projected radius of $\sim$2.4 kpc is characterized as LINER-like line ratios while the outskirts extended to $\sim$15 kpc show Seyfert-like line ratios. The [OIII]$λ$5007 luminosity of the Seyfert regions is a factor of 37 (2) higher than the LINER regions without (with) dust attenuation correction, suggesting that the AGN activity decreases at least $\sim$8 $\times$ 10$^3$ yrs ($\sim$2.4 kpc/light-speed) ago. We model the emission line spectra in the central region with double Gaussian components (a narrow core and a broad wing) and analyze the properties of each component. The narrow core component mostly co-rotates with the stellar disc, whereas the broad wing component with a median of the velocity dispersion $\sim$300 km s$^{-1}$ is related to a wind outflow. The kinematic position angle (PA) of the ionized gas shows a $\sim$20° twist from the galaxy center to 1.5 effective radius. The median of the PA difference between the gas and stellar components is as large as $\sim$50° within 0.4 effective radius. The tidal feature in DESI image and star-gas misalignment suggest this galaxy is a merger remnant. Combining all these observational results as well as public available X-ray and MIR luminosities, we confirm this is a fading AGN, the merger process kick-started the central engine to quasar phase which ionized gas composed of tidal debris, and now the activity of the central black hole decreases. The discontinuity in [OIII]$λ$5007 flux and EQW maps is due to multiple AGN outbursts triggered by merger remnant gas inflows. △ Less

Submitted 14 March, 2024; originally announced March 2024.

Comments: Accepted for publication in MNRAS. 12 pages, 10 figures, 1 table

The true number density of massive galaxies in the early Universe revealed by JWST/MIRI

Authors: Tao Wang, Hanwen Sun, Luwenjia Zhou, Ke Xu, Cheng Cheng, Zhaozhou Li, Yangyao Chen, H. J. Mo, Avishai Dekel, Xianzhong Zheng, Zheng Cai, Tiacheng Yang, Y. -S. Dai, David Elbaz, J. -S. Huang

Abstract: 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… ▽ More 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 μ$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 $Λ$CDM cosmology, our results require a moderate increase in the baryon-to-star conversion efficiency ($ε$) towards higher redshifts and higher $M_{\star}$. For the most massive galaxies at $z\sim 8$, the required $ε$ is $\sim 0.3$, in comparison to $ε\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. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: 23 pages, 10 figures, submitted

ELUCID VIII: Simulating the Coma Galaxy Cluster to Calibrate Model and Understand Feedback

Authors: Xiong Luo, Huiyuan Wang, Weiguang Cui, Houjun Mo, RenJie Li, Yipeng Jing, Neal Katz, Romeel Davé, Xiaohu Yang, Yangyao Chen, Hao Li, Shuiyao Huang

Abstract: We conducted an investigation of the Coma cluster of galaxies by running a series of constrained hydrodynamic simulations with GIZMO-SIMBA and GADGET-3, based on initial conditions reconstructed from the SDSS survey volume in the ELUCID project. We compared simulation predictions and observations for galaxies, ICM and IGM in and around the Coma cluster to constrain galaxy formation physics. Our re… ▽ More We conducted an investigation of the Coma cluster of galaxies by running a series of constrained hydrodynamic simulations with GIZMO-SIMBA and GADGET-3, based on initial conditions reconstructed from the SDSS survey volume in the ELUCID project. We compared simulation predictions and observations for galaxies, ICM and IGM in and around the Coma cluster to constrain galaxy formation physics. Our results demonstrate that this type of constrained investigation allows us to probe in more detail the implemented physical processes, because the comparison between simulations and observations is free of cosmic variance and hence can be conducted in a ''one-to-one'' manner. We found that an increase in the earlier star formation rate and the supernova feedback of the original GIZMO-SIMBA model is needed to match observational data on stellar, ISM and ICM metallicity. The simulations without AGN feedback can well reproduce the observational ICM electron density, temperature, and entropy profiles, ICM substructures, and the IGM temperature-density relation, while the ones with AGN feedback usually fail. However, one requires something like AGN feedback to reproduce a sufficiently large population of quiescent galaxies, particularly in low-density regions. The constrained simulations of the Coma cluster thus provide a test bed to understand processes that drive galaxy formation and evolution. △ Less

Submitted 26 March, 2024; v1 submitted 26 January, 2024; originally announced January 2024.

Comments: 21 pages, 9 figures

Post-starburst galaxies in SDSS-IV MaNGA: Two broad categories of evolutionary pathways

Authors: Zhuo Cheng, Cheng Li, Niu Li, Renbin Yan, Houjun Mo

Abstract: We study the size-mass relation (SMR) and recent star formation history (SFH) of post-starburst (PSB) galaxies in the local Universe, using spatially resolved spectroscopy from the final data release of MaNGA. Our sample includes 489 PSB galaxies: 94 cPSB galaxies with central PSB regions, 85 rPSB galaxies with ring-like PSB regions and 310 iPSB galaxies with irregular PSB regions. When compared t… ▽ More We study the size-mass relation (SMR) and recent star formation history (SFH) of post-starburst (PSB) galaxies in the local Universe, using spatially resolved spectroscopy from the final data release of MaNGA. Our sample includes 489 PSB galaxies: 94 cPSB galaxies with central PSB regions, 85 rPSB galaxies with ring-like PSB regions and 310 iPSB galaxies with irregular PSB regions. When compared to control galaxies of similar SFR, redshift and mass, a similar SMR is found for all types of PSB samples except the cPSB galaxies which have smaller sizes at intermediate masses ($9.5\lesssim \log_{10}(\rm M_\ast/M_\odot)\lesssim 10.5$). The iPSB galaxies in the star-forming sequence (iPSB-SF) show no/weak gradients in $\textrm{D}_{n}(4000)$, $\textrm{EW}(\textrm{H}δ_{A})$ and $\textrm{EW}(\textrm{H}α)$, consistent with the global star-forming status of this type of galaxies, while the quiescent iPSB (iPSB-Q) sample shows negative gradients in $\textrm{D}_{n}(4000)$ and positive gradients in $\textrm{EW}(\textrm{H}δ_{A})$, indicating older stellar populations in the inner regions. Both cPSB and rPSB samples show positive gradients in $\textrm{D}_{n}(4000)$ and negative gradients in $\textrm{EW}(\textrm{H}δ_{A})$, indicating younger stellar populations in the inner regions. These results imply that the four types of PSB galaxies can be broadly divided into two distinct categories in terms of evolutionary pathway: (1) iPSB-SF and iPSB-Q which have SMRs and SFHs similar to control galaxies, preferring an inside-out quenching process, (2) rPSB and cPSB which appear to be different stages of the same event, likely to follow the outside-in quenching process driven by disruption events such as mergers that result in a more compact structure as quenching proceeds. △ Less

Submitted 28 January, 2024; v1 submitted 6 December, 2023; originally announced December 2023.

Comments: 20 pages, 10 figures, accepted for publication in ApJ

On the existence, rareness and uniqueness of quenched HI-rich galaxies in the local Universe

Authors: Xiao Li, Cheng Li, H. J. Mo, Jianhong Hu, Jing Wang, Ting Xiao

Abstract: 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… ▽ More 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. △ Less

Submitted 18 February, 2024; v1 submitted 6 December, 2023; originally announced December 2023.

Comments: 25 pages, 12 figures, accepted for publication in ApJ, proof completed. Correct some typos in the text. Correct a typo in figure 2

A two-phase model of galaxy formation: II. The size-mass relation of dynamically hot galaxies

Authors: Yangyao Chen, Houjun Mo, Huiyuan Wang

Abstract: In Paper-I we developed a two-phase model to connect dynamically hot galaxies (such as ellipticals and bulges) with the formation of self-gravitating gas clouds (SGCs) associated with the fast assembly of dark matter halos. Here we explore the implications of the model for the size-stellar mass relation of dynamically hot galaxies. Star-forming sub-clouds resulting from the fragmentation of the tu… ▽ More In Paper-I we developed a two-phase model to connect dynamically hot galaxies (such as ellipticals and bulges) with the formation of self-gravitating gas clouds (SGCs) associated with the fast assembly of dark matter halos. Here we explore the implications of the model for the size-stellar mass relation of dynamically hot galaxies. Star-forming sub-clouds resulting from the fragmentation of the turbulent SGC inherit its spatial structure and dynamical hotness, producing a `homologous' relation, $r_{\rm f}\approx\, 100 r_{\rm bulge}$, between the size of a dynamically hot galaxy ($r_{\rm bulge}$) and that of its host halo assembled in the fast regime ($r_{\rm f}$), independent of redshift and halo mass. This relation is preserved by the `dry' expansion driven by dynamical heating when a galaxy becomes gas-poor due to inefficient cooling, and is frozen due to the stop of bulge growth during the slow assembly regime of the halo. The size-stellar mass relation is thus a simple combination of the galaxy-halo homology and the non-linear stellar mass-halo mass relation. Using a set of halo assembly histories we reproduce all properties in the observed size-mass relation of dynamically hot galaxies, including the flattening in the low-mass end and the upturn in the massive end. The prediction matches observational data currently available to $z \approx 4$, and can be tested in the future at higher $z$. Our results indicate that the sizes of dynamically hot galaxies are produced by the dissipation and collapse of gas in halos to establish SGCs in which stars form. △ Less

Submitted 3 July, 2024; v1 submitted 20 November, 2023; originally announced November 2023.

Comments: 11 pages, 5 figures, 1 table; submitted to MNRAS

A two-phase model of galaxy formation: I. The growth of galaxies and supermassive black holes

Authors: Houjun Mo, Yangyao Chen, Huiyuan Wang

Abstract: We develop a model for galaxy formation and the growth of supermassive black holes (SMBHs), based on the fact that cold dark matter (CDM) halos form their gravitational potential wells through a fast phase with rapid change in the potential, and that the high universal baryon fraction makes cooled gas in halos self-gravitating and turbulent before it can form rotation-supported disks. Gas fragment… ▽ More We develop a model for galaxy formation and the growth of supermassive black holes (SMBHs), based on the fact that cold dark matter (CDM) halos form their gravitational potential wells through a fast phase with rapid change in the potential, and that the high universal baryon fraction makes cooled gas in halos self-gravitating and turbulent before it can form rotation-supported disks. Gas fragmentation produces sub-clouds so dense that cloud-cloud collision and drag on clouds are not significant, producing a dynamically hot system of sub-clouds that form stars and move ballistically to feed the central SMBH. Active galactic nucleus (AGN) and supernova (SN) feedback is effective only in the fast phase, and the cumulative effects are to regulate star formation and SMBH growth, as well as to reduce the amount of cold gas in halos to allow the formation of globally stable disks. Using a set of halo assembly histories, we demonstrate that the model can reproduce a number of observations, including correlations among SMBH mass, stellar mass of galaxies and halo mass, the number densities of galaxies and SMBH, as well as their evolution over the cosmic time. △ Less

Submitted 18 July, 2024; v1 submitted 8 November, 2023; originally announced November 2023.

Comments: 32 pages, 22 figures, 1 table. MNRAS accepted. Comments are welcomed

An efficient and robust method to estimate halo concentration

Authors: Kai Wang, H. J. Mo, Yangyao Chen, Joop Schaye

Abstract: 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πr^3ρ(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… ▽ More 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πr^3ρ(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. △ Less

Submitted 19 December, 2023; v1 submitted 29 September, 2023; originally announced October 2023.

Comments: 17 pages, 9 + 9 figures. Main figure: Fig. 3. MNRAS accepted

Characterizing the assembly of dark matter halos with protohalo size histories: I. Redshift evolution, relation to descendant halos, and halo assembly bias

Authors: Kai Wang, H. J. Mo, Yangyao Chen, Huiyuan Wang, Xiaohu Yang, Jiaqi Wang, Yingjie Peng, Zheng Cai

Abstract: 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 centra… ▽ More 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. △ Less

Submitted 12 January, 2024; v1 submitted 2 September, 2023; originally announced September 2023.

Comments: 20 pages, 12 + 10 figures, MNRAS accepted!

Halo mass-observable proxy scaling relations and their dependencies on galaxy and group properties

Authors: Ziwen Zhang, Huiyuan Wang, Wentao Luo, Houjun Mo, Jun Zhang, Xiaohu Yang, Hao Li, Qinxun Li

Abstract: Based on the DECaLS shear catalog, we study the scaling relations between halo mass($M_{\rm h}$) and various proxies for SDSS central galaxies, including stellar mass($M_*$), stellar velocity dispersion($σ_*$), abundance matching halo mass($M_{\rm AM}$) and satellite velocity dispersion($σ_{\rm s}$), and their dependencies on galaxy and group properties. In general, these proxies all have strong p… ▽ More Based on the DECaLS shear catalog, we study the scaling relations between halo mass($M_{\rm h}$) and various proxies for SDSS central galaxies, including stellar mass($M_*$), stellar velocity dispersion($σ_*$), abundance matching halo mass($M_{\rm AM}$) and satellite velocity dispersion($σ_{\rm s}$), and their dependencies on galaxy and group properties. In general, these proxies all have strong positive correlations with $M_{\rm h}$, consistent with previous studies. We find that the $M_{\rm h}$-$M_*$ and $M_{\rm h}$-$σ_*$ relations depend strongly on group richness($N_{\rm sat}$), while the $M_{\rm h}$-$M_{\rm AM}$ and $M_{\rm h}$-$σ_{\rm s}$ relations are independent of it. Moreover, the dependence on star formation rate(SFR) is rather weak in the $M_{\rm h}$-$σ_*$ and $M_{\rm h}$-$σ_{\rm s}$ relations, but very prominent in the other two. $σ_{\rm s}$ is thus the best proxy among them, and its scaling relation is in good agreement with hydro-dynamical simulations. However, estimating $σ_{\rm s}$ accurately for individual groups/clusters is challenging because of interlopers and the requirement for sufficient satellites. We construct new proxies by combining $M_*$, $σ_*$, and $M_{\rm AM}$, and find the proxy with 30\% contribution from $M_{\rm AM}$ and 70\% from $σ_*$ can minimize the dependence on $N_{\rm sat}$ and SFR. We obtain the $M_{\rm h}$-supermassive black hole(SMBH) mass relation via the SMBH scaling relation and find indications for rapid and linear growth phases for SMBH. We also find that correlations among $M_{\rm h}$, $M_*$ and $σ_*$ change with $M_*$, indicating that different processes drive the growth of galaxies and SMBH at different stages. △ Less

Submitted 3 November, 2023; v1 submitted 11 May, 2023; originally announced May 2023.

Comments: 22 pages, 12 figures and 2 tables, accepted for publication in ApJ

Dynamical hotness, star formation quenching and growth of supermassive black holes

Authors: Hui Hong, Huiyuan Wang, H. J. Mo, Ziwen Zhang, Guangwen Chen, Wentao Luo, Tinggui Wang, Pengfei Li, Renjie Li, Yao yao, Aoxiang Jiang

Abstract: A stellar system is dynamically hot when its kinetic energy is dominated by random motion represented by the velocity dispersion $σ_{\rm hot} (M_*)$. We use MaNGA data to obtain inner and outer dispersion of a galaxy, $σ_{\rm in}$ and $σ_{\rm out}$, to characterize its dynamical status and study its connection with star formation quenching and the growth of supermassive black hole (SMBH). We divid… ▽ More A stellar system is dynamically hot when its kinetic energy is dominated by random motion represented by the velocity dispersion $σ_{\rm hot} (M_*)$. We use MaNGA data to obtain inner and outer dispersion of a galaxy, $σ_{\rm in}$ and $σ_{\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 $σ_{\rm in}/σ_{\rm hot}$-$σ_{\rm out}/σ_{\rm hot}$ diagram is L-shaped, consisting of a horizontal sequence ($σ_{\rm out}/σ_{\rm hot}\sim0$) and a vertical sequence ($σ_{\rm in}/σ_{\rm hot}\sim1$). FQGs and QCCs are located at the top of vertical sequence, $σ_{\rm out}/σ_{\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(σ_{\rm in}/σ_{\rm hot})\sim-0.3$, and along vertical sequence at $\log(σ_{\rm out}/σ_{\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. △ Less

Submitted 19 July, 2023; v1 submitted 4 May, 2023; originally announced May 2023.

Comments: 24 pages, 19 figures, 1 table, accepted by ApJ

Massive Dark Matter Halos at High Redshift: Implications for Observations in the JWST Era

Authors: Yangyao Chen, H. J. Mo, Kai Wang

Abstract: The presence of massive galaxies at high $z$ as recently observed by JWST appears to contradict the current $Λ$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… ▽ More The presence of massive galaxies at high $z$ as recently observed by JWST appears to contradict the current $Λ$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 $ρ_*(>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 $ε_* \sim 0.5$, cosmic variance alone reduces the tension to a $2σ$ 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 $ρ_*(>M_*)$ at $z \sim 7 - 10$ within $2σ$. 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. △ Less

Submitted 18 October, 2023; v1 submitted 26 April, 2023; originally announced April 2023.

Comments: 19 pages, 15 figures. MNRAS accepted

A Conditional Abundance Matching Method of Extending Simulated Halo Merger Trees to Resolve Low-Mass Progenitors and Sub-halos

Authors: Yangyao Chen, H. J. Mo, Cheng Li, Kai Wang, Huiyuan Wang, Xiaohu Yang

Abstract: 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… ▽ More 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. △ Less

Submitted 18 October, 2023; v1 submitted 21 January, 2023; originally announced January 2023.

Comments: 24 pages, 14 figures, 3 tables. MNRAS accepted

Mapping dust attenuation and the 2175 Å bump at kpc scales in nearby galaxies

Authors: Shuang Zhou, Cheng Li, Niu Li, Houjun Mo, Renbin Yan, Michael Eracleous, Mallory Molina, Caryl Gronwall, Nikhil Ajgaonkar, Zhuo Cheng, Ruonan Guo

Abstract: We develop a novel approach to measure dust attenuation properties of galaxies, including the dust opacity, shape of the attenuation curve and the strength of the 2175Å absorption feature. From an observed spectrum, the method uses a model-independent approach to derive a relative attenuation curve, with absolute amplitude calibrated using NIR photometry. The dust-corrected spectrum is fitted with… ▽ More We develop a novel approach to measure dust attenuation properties of galaxies, including the dust opacity, shape of the attenuation curve and the strength of the 2175Å absorption feature. From an observed spectrum, the method uses a model-independent approach to derive a relative attenuation curve, with absolute amplitude calibrated using NIR photometry. The dust-corrected spectrum is fitted with stellar population models to derive the dust-free model spectrum, which is compared with the observed SED/spectrum from NUV to NIR to determine dust attenuation properties. We apply this method to investigate dust attenuation on kpc scales, using a sample of 134 galaxies with integral field spectroscopy from MaNGA, NIR imaging from 2MASS, and NUV imaging from Swift/UVOT. We find the attenuation curve slope and the 2175Å bump in both optical and NUV span a wide range at kpc scales. The slope is shallower at higher optical opacity, regardless of the specific star formation rate (sSFR), minor-to-major axis ratio (b/a) of galaxies and the location of spaxels within individual galaxies. The 2175Å bump presents a strong negative correlation with the sSFR, while the correlations with the optical opacity, b/a and the location within individual galaxies are all weak. All these trends appear to be independent of the stellar mass of galaxies. Our results support the scenario that the variation of the 2175Å bump is driven predominantly by processes related to star formation, such as the destruction of small dust grains by UV radiation in star-forming regions. △ Less

Submitted 19 September, 2023; v1 submitted 4 December, 2022; originally announced December 2022.

Comments: 21 pages,9 figures, APJ accepted version

Relating galaxies across different redshift to study galaxy evolution

Authors: Kai Wang, Houjun Mo, Cheng Li, Yangyao Chen

Abstract: We propose a general framework leveraging the halo-galaxy connection to link galaxies observed at different redshift in a statistical way, and use the link to infer the redshift evolution of the galaxy population. Our tests based on hydrodynamic simulations show that our method can accurately recover the stellar mass assembly histories up to $z\sim 3$ for present star-forming and quiescent galaxie… ▽ More We propose a general framework leveraging the halo-galaxy connection to link galaxies observed at different redshift in a statistical way, and use the link to infer the redshift evolution of the galaxy population. Our tests based on hydrodynamic simulations show that our method can accurately recover the stellar mass assembly histories up to $z\sim 3$ for present star-forming and quiescent galaxies down to $10^{10}h^{-1}M_{\odot}$. Applying the method to observational data shows that the stellar mass evolution of the main progenitors of galaxies depends strongly on the properties of descendants, such as stellar mass, halo mass, and star formation states. Galaxies hosted by low-mass groups/halos at the present time have since $z\sim 1.8$ grown their stellar mass $\sim 2.5$ times as fast as those hosted by massive clusters. This dependence on host halo mass becomes much weaker for descendant galaxies with similar star formation states. Star-forming galaxies grow about 2-4 times faster than their quiescent counterparts since $z\sim 1.8$. Both TNG and EAGLE simulations over-predict the progenitor stellar mass at $z>1$, particularly for low-mass descendants. △ Less

Submitted 8 February, 2023; v1 submitted 1 November, 2022; originally announced November 2022.

Comments: 16 pages, 10+4 figures, published on MNRAS

Galaxy populations in groups and clusters: evidence for a characteristic stellar mass scale at $M_\ast\sim 10^{9.5}M_\odot$

Authors: Jiacheng Meng, Cheng Li, Houjun Mo, Yangyao Chen, Zhen Jiang, Lizhi Xie

Abstract: We use the most recent data release (DR9) of the DESI legacy imaging survey and SDSS galaxy groups to measure the conditional luminosity function (CLF) for groups with halo mass $M_{\rm h}\ge 10^{12}M_{\odot}$ and redshift $0.01\le z\le 0.08$, down to a limiting $r$-band magnitude of $M_{\rm r}=-10\sim-12$. For a given halo mass we measure the CLF for the total satellite population, as well as sep… ▽ More We use the most recent data release (DR9) of the DESI legacy imaging survey and SDSS galaxy groups to measure the conditional luminosity function (CLF) for groups with halo mass $M_{\rm h}\ge 10^{12}M_{\odot}$ and redshift $0.01\le z\le 0.08$, down to a limiting $r$-band magnitude of $M_{\rm r}=-10\sim-12$. For a given halo mass we measure the CLF for the total satellite population, as well as separately for the red and blue populations classified using the $(g-z)$ color. We find a clear faint-end upturn in the CLF of red satellites, with a slope $α\approx-1.8$ which is almost independent of halo mass. This faint-end upturn is not seen for blue satellites and for the total population. Our stellar population synthesis modeling shows that the $(g-z)$ color provides a clean red/blue division, and that group galaxies in the red population defined by $(g-z)$ are all dominated by old stellar populations. The fraction of old galaxies as a function of galaxy luminosity shows a minimum at a luminosity $M_{\rm r}\sim-18$, corresponding to a stellar mass $M_\ast\sim10^{9.5}M_\odot$. This mass scale is independent of halo mass and is comparable to the characteristic luminosity at which galaxies show a dichotomy in surface brightness and size, suggesting that the dichotomy in the old fraction and in galaxy structure may have a common origin. The rising of the old fraction at the faint end for Milky Way (MW)-sized halos found here is in good agreement with the quenched fraction measured both for the MW/M31 system and from the ELVES survey. We discuss the implications of our results for the formation and evolution of low-mass galaxies, and for the stellar mass functions of low-mass galaxies to be observed at high redshift. △ Less

Submitted 21 December, 2022; v1 submitted 31 October, 2022; originally announced October 2022.

Comments: 26 pages, 13 figures, accepted by ApJ

Shock-induced stripping of satellite ISM/CGM in IllustrisTNG clusters at $z\sim0$

Authors: Hao Li, Huiyuan Wang, H. J. Mo, Yuan Wang, Xiong Luo, Renjie Li

Abstract: 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… ▽ More 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. △ Less

Submitted 29 November, 2022; v1 submitted 16 September, 2022; originally announced September 2022.

Comments: 14 pages,9 figures. Accepted for publication in ApJ

Conditional HI mass functions and the HI-to-halo mass relation in the local Universe

Authors: Xiao Li, Cheng Li, H. J. Mo, Ting Xiao, Jing Wang

Abstract: 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 cons… ▽ More 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. △ Less

Submitted 13 December, 2022; v1 submitted 15 September, 2022; originally announced September 2022.

Comments: 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

Anisotropy and characteristic scales in halo density gradient profiles

Authors: X. Wang, H. Wang, H. J. Mo

Abstract: 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 ki… ▽ More 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. △ Less

Submitted 8 September, 2022; v1 submitted 24 June, 2022; originally announced June 2022.

Comments: 15 pages, 12 figures, accepted by A&A

Journal ref: A&A 667, A99 (2022)

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

Authors: Renjie Li, Huiyuan Wang, H. J. Mo, Shuiyao Huang, Neal Katz, Xiong Luo, Weiguang Cui, Hao Li, Xiaohu Yang, Ning Jiang, Yuning Zhang

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

Submitted 7 September, 2022; v1 submitted 16 June, 2022; originally announced June 2022.

Comments: 21 pages, 15 figures, 1 table

Journal ref: APJ, 2022

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

Authors: Qingyang Li, Xiaohu Yang, Chengze Liu, Yipeng Jing, Min He, Jiasheng Huang, Y. Sophia Dai, Marcin Sawicki, Stephane Arnouts, Stephen Gwyn, Thibaud Moutard, H. J. Mo, Kai Wang, Antonios Katsianis, Weiguang Cui, Jiaxin Han, I-Non Chiu, Yizhou Gu, Haojie Xu

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

Submitted 11 May, 2022; originally announced May 2022.

Comments: 25 pages, 16 figures, 2 tables, accepted for publication in ApJ

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

Authors: Ziwen Zhang, Huiyuan Wang, Wentao Luo, Jun Zhang, H. J. Mo, YiPeng Jing, Xiaohu Yang, Hao Li

Abstract: 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… ▽ More 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. △ Less

Submitted 19 May, 2022; v1 submitted 9 December, 2021; originally announced December 2021.

Comments: Accepted for publication in A&A; 15 pages, 13 figures

Journal ref: A&A 663, A85 (2022)

MAHGIC: A Model Adapter for the Halo-Galaxy Inter-Connection

Authors: Yangyao Chen, H. J. Mo, Cheng Li, Kai Wang, Huiyuan Wang, Xiaohu Yang, Youcai Zhang, Neal Katz

Abstract: 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 hydrodyn… ▽ More 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. △ Less

Submitted 13 August, 2021; v1 submitted 7 June, 2021; originally announced June 2021.

Comments: 21 pages, 16 figures, 2 tables, MNRAS accepted

Finding proto-clusters to trace galaxy evolution: I. The finder and its performance

Authors: Kai Wang, H. J. Mo, Cheng Li, Yangyao Chen

Abstract: 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.… ▽ More 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. △ Less

Submitted 3 June, 2021; v1 submitted 25 April, 2021; originally announced April 2021.

Comments: 15 pages, 10 + 4 figures, MNRAS accepted, Comments are welcome

Evaluating the origins of the secondary bias based on the correlation of halo properties with the linear density field

Authors: Xiaoyu Wang, Huiyuan Wang, H. J. Mo, JingJing Shi, Yipeng Jing

Abstract: 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 ($λ$) 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… ▽ More 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 ($λ$) 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 $δ_{\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 $λ$. The external correlation, which describes the correlation of halo properties with linear overdensity at $R>R_{\rm L}$ for given $δ_{\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. △ Less

Submitted 27 September, 2021; v1 submitted 20 April, 2021; originally announced April 2021.

Comments: 14 pages,14 figures, accepted by A&A

Journal ref: A&A 654, A67 (2021)

Estimating dust attenuation from galactic spectra. II. Stellar and gas attenuation in star-forming and diffuse ionized gas regions in MaNGA

Authors: Niu Li, Cheng Li, Houjun Mo, Shuang Zhou, Fu-heng Liang, Médéric Boquien, Niv Drory, José G. Fernández-Trincado, Michael Greener, Rogério Riffel

Abstract: We investigate the dust attenuation in both stellar populations and ionized gas in kpc-scale regions in nearby galaxies, using integral field spectroscopy data from MaNGA MPL-9. We identify star-forming (HII) and diffuse ionized gas (DIG) regions from MaNGA datacubes. From the stacked spectrum of each region, we measure the stellar attenuation, $E(B-V)_{\rm star}$, using the technique developed by… ▽ More We investigate the dust attenuation in both stellar populations and ionized gas in kpc-scale regions in nearby galaxies, using integral field spectroscopy data from MaNGA MPL-9. We identify star-forming (HII) and diffuse ionized gas (DIG) regions from MaNGA datacubes. From the stacked spectrum of each region, we measure the stellar attenuation, $E(B-V)_{\rm star}$, using the technique developed by Li et al.(2020), as well as the gas attenuation, $E(B-V)_{\rm gas}$, from the Balmer decrement. We then examine the correlation of $E(B-V)_{\rm star}$, $E(B-V)_{\rm gas}$, $E(B-V)_{\rm gas}-E(B-V)_{\rm star}$ and $E(B-V)_{\rm star}/E(B-V)_{\rm gas}$ with 16 regional/global properties, and for regions with different $\rm Hα$ surface brightnesses ($Σ_{\rm Hα}$). We find a stronger correlation between $E(B-V)_{\rm star}$ and $E(B-V)_{\rm gas}$ in regions of higher $Σ_{\rm Hα}$. Luminosity-weighted age ($t_L$) is found to be the property that is the most strongly correlated with $E(B-V)_{\rm star}$, and consequently with $E(B-V)_{\rm gas}-E(B-V)_{\rm star}$ and $E(B-V)_{\rm star}/E(B-V)_{\rm gas}$. At fixed $Σ_{\rm Hα}$, $\log_{10}t_L$ is linearly and negatively correlated with $E(B-V)_{\rm star}/E(B-V)_{\rm gas}$ at all ages. Gas-phase metallicity and ionization level are important for the attenuation in the gas. Our results indicate that the ionizing source for DIG regions is likely distributed in the outer-skirt of galaxies, while for HII regions our results can be well explained by the two-component dust model of Charlot & Fall (2000). △ Less

Submitted 7 June, 2021; v1 submitted 28 February, 2021; originally announced March 2021.

Comments: accepted by ApJ

SDSS-IV MaNGA: the physical origin of off-galaxy H$α$ blobs in the local Universe

Authors: Xihan Ji, Cheng Li, Renbin Yan, Houjun Mo, Lihwai Lin, Hu Zou, Jianhui Lian, David V. Stark, Rogemar A. Riffel, Hsi-An Pan, Dmitry Bizyaev, Kevin Bundy

Abstract: H$α$ blobs are off-galaxy emission-line regions with weak or no optical counterparts. They are mostly visible in H$α$ line, appearing as concentrated blobs. Such unusual objects have been rarely observed and studied, and their physical origin is still unclear. We have identified 13 H$α$ blobs in the public data of MaNGA survey, by visually inspecting both the optical images and the spatially resol… ▽ More H$α$ blobs are off-galaxy emission-line regions with weak or no optical counterparts. They are mostly visible in H$α$ line, appearing as concentrated blobs. Such unusual objects have been rarely observed and studied, and their physical origin is still unclear. We have identified 13 H$α$ blobs in the public data of MaNGA survey, by visually inspecting both the optical images and the spatially resolved maps of H$α$ line for $\sim 4600$ galaxy systems. Among the 13 H$α$ blobs, 2 were reported in previously MaNGA-based studies and 11 are newly discovered. This sample, though still small in size, is by far the largest sample with both deep imaging and integral field spectroscopy. Therefore, for the first time we are able to perform statistical studies to investigate the physical origin of H$α$ blobs. We examine the physical properties of these H$α$ blobs and their associated galaxies, including their morphology, environments, gas-phase metallicity, kinematics of ionized gas, and ionizing sources. We find that the H$α$ blobs in our sample can be broadly divided into two groups. One is associated with interacting/merging galaxy systems, of which the ionization is dominated by shocks or diffuse ionized gas. It is likely that these H$α$ blobs used to be part of their nearby galaxies, but were stripped away at some point due to tidal interactions. The other group is found in gas-rich systems, appearing as low-metallicity star-forming regions that are visually detached from the main galaxy. These H$α$ blobs could be associated with faint disks, spiral arms, or dwarf galaxies. △ Less

Submitted 24 September, 2021; v1 submitted 18 January, 2021; originally announced January 2021.

Comments: 26 pages, 18 figures, accepted for publication in MNRAS

Wolf-Rayet galaxies in SDSS-IV MaNGA. II. Metallicity dependence of the high-mass slope of the stellar initial mass function

Authors: Fu-Heng Liang, Cheng Li, Niu Li, Shuang Zhou, Renbin Yan, Houjun Mo, Wei Zhang

Abstract: As hosts of living high-mass stars, Wolf-Rayet (WR) regions or WR galaxies are ideal objects for constraining the high-mass end of the stellar initial mass function (IMF). We construct a large sample of 910 WR galaxies/regions that cover a wide range of stellar metallicity (from Z~0.001 up to Z~0.03), by combining three catalogs of WR galaxies/regions previously selected from the SDSS and SDSS-IV/… ▽ More As hosts of living high-mass stars, Wolf-Rayet (WR) regions or WR galaxies are ideal objects for constraining the high-mass end of the stellar initial mass function (IMF). We construct a large sample of 910 WR galaxies/regions that cover a wide range of stellar metallicity (from Z~0.001 up to Z~0.03), by combining three catalogs of WR galaxies/regions previously selected from the SDSS and SDSS-IV/MaNGA surveys. We measure the equivalent widths of the WR blue bump at ~4650 A for each spectrum. They are compared with predictions from stellar evolutionary models Starburst99 and BPASS, with different IMF assumptions (high-mass slope α of the IMF ranging from 1.0 up to 3.3). Both singular evolution and binary evolution are considered. We also use a Bayesian inference code to perform full spectral fitting to WR spectra with stellar population spectra from BPASS as fitting templates. We then make model selection among different α assumptions based on Bayesian evidence. These analyses have consistently led to a positive correlation of IMF high-mass slope α with stellar metallicity Z, i.e. with steeper IMF (more bottom-heavy) at higher metallicities. Specifically, an IMF with α=1.00 is preferred at the lowest metallicity (Z~0.001), and a Salpeter or even steeper IMF is preferred at the highest metallicity (Z~0.03). These conclusions hold even when binary population models are adopted. △ Less

Submitted 2 October, 2021; v1 submitted 8 January, 2021; originally announced January 2021.

Comments: 22 pages, 11 figures, accepted by ApJ

An Extended Halo-based Group/Cluster finder: application to the DESI legacy imaging surveys DR8

Authors: Xiaohu Yang, Haojie Xu, Min He, Yizhou Gu, Antonios Katsianis, Jiacheng Meng, Feng Shi, Hu Zou, Youcai Zhang, Chengze Liu, Zhaoyu Wang, Fuyu Dong, Yi Lu, Qingyang Li, Yangyao Chen, Huiyuan Wang, Houjun Mo, Jian Fu, Hong Guo, Alexie Leauthaud, Yu Luo, Jun Zhang, Ying Zu

Abstract: We extend the halo-based group finder developed by \citet[][]{Yang2005a} to use data {\it simultaneously} with either photometric or spectroscopic redshifts. A mock galaxy redshift survey constructed from a high-resolution N-body simulation is used to evaluate the performance of this extended group finder. For galaxies with magnitude ${\rm z\le 21}$ and redshift $0<z\le 1.0$ in the DESI legacy ima… ▽ More We extend the halo-based group finder developed by \citet[][]{Yang2005a} to use data {\it simultaneously} with either photometric or spectroscopic redshifts. A mock galaxy redshift survey constructed from a high-resolution N-body simulation is used to evaluate the performance of this extended group finder. For galaxies with magnitude ${\rm z\le 21}$ and redshift $0<z\le 1.0$ in the DESI legacy imaging surveys (the Legacy Surveys), our group finder successfully identifies more than 60\% of the members in about $90\%$ of halos with mass $\ga 10^{12.5}\msunh$. Detected groups with mass $\ga 10^{12.0}\msunh$ have a purity (the fraction of true groups) greater than 90\%. The halo mass assigned to each group has an uncertainty of about 0.2 dex at the high mass end $\ga 10^{13.5}\msunh$ and 0.40 dex at the low mass end. Groups with more than 10 members have a redshift accuracy of $\sim 0.008$. We apply this group finder to the Legacy Surveys DR8 and find 5.2 Million groups with at least 3 members. About 387,000 of these groups have at least 10 members. The resulting catalog containing 3D coordinates, richness, halo masses, and total group luminosities, is made publicly available. △ Less

Submitted 7 April, 2022; v1 submitted 29 December, 2020; originally announced December 2020.

Comments: Due to a bug in our code, our calculation of the completeness and interloper fraction of the group members was not correct, which resulted in our choice of a low background B=2.5. In this new version of the paper, we updated results with the theoretical fiducial background B=10. Most of the changes in the figures are virtually too small to notice. This version supersedes the one published in ApJ

Hosts and triggers of AGNs in the Local Universe

Authors: Ziwen Zhang, Huiyuan Wang, Wentao Luo, H. J. Mo, Zhixiong Liang, Ran Li, Xiaohu Yang, Tinggui Wang, Hongxin Zhang, Hui Hong, Xiaoyu Wang, Enci Wang, Pengfei Li, JingJing Shi

Abstract: 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 t… ▽ More 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. △ Less

Submitted 30 April, 2021; v1 submitted 19 December, 2020; originally announced December 2020.

Comments: accepted by A&A, 17 pages, 11 figures

Journal ref: A&A 650, A155 (2021)

Star formation histories of massive red spiral galaxies in the local universe

Authors: Shuang Zhou, Cheng Li, Cai-Na Hao, Rui Guo, Houjun Mo, Xiaoyang Xia

Abstract: We investigate the star formation histories (SFHs) of massive red spiral galaxies with stellar mass $M_\ast>10^{10.5}M_\odot$, and make comparisons with blue spirals and red ellipticals of similar masses. We make use of the integral field spectroscopy from the SDSS-IV/DR15 MaNGA sample, and estimate spatially resolved SFHs and stellar population properties of each galaxy by applying a Bayesian spe… ▽ More We investigate the star formation histories (SFHs) of massive red spiral galaxies with stellar mass $M_\ast>10^{10.5}M_\odot$, and make comparisons with blue spirals and red ellipticals of similar masses. We make use of the integral field spectroscopy from the SDSS-IV/DR15 MaNGA sample, and estimate spatially resolved SFHs and stellar population properties of each galaxy by applying a Bayesian spectral fitting code to the MaNGA spectra. We find that both red spirals and red ellipticals have experienced only one major star formation episode at early times, and the result is independent of the adopted SFH model. On average, more than half of their stellar masses were formed $>$10 Gyrs ago, and more than 90\% were formed $>6$ Gyrs ago. The two types of galaxies show similarly flat profiles in a variety of stellar population parameters: old stellar ages indicated by $D4000$ (the spectral break at around 4000Å), high stellar metallicities, large Mgb/Fe ratios indicating fast formation, and little stellar dust attenuation. In contrast, although blue spirals also formed their central regions $>$10 Gyrs ago, both their central regions and outer disks continuously form stars over a long timescale. Our results imply that, massive red spirals are likely to share some common processes of formation (and possibly quenching) with massive red ellipticals in the sense that both types were formed at $z > 2$ through a fast formation process.Possible mechanisms for the formation and quenching of massive red spirals are discussed. △ Less

Submitted 30 May, 2021; v1 submitted 27 November, 2020; originally announced November 2020.

Comments: 17 pages, 6 figures, APJ accepted

Elucidating Galaxy Assembly Bias in SDSS

Authors: Andrés N. Salcedo, Ying Zu, Youcai Zhang, Huiyuan Wang, Xiaohu Yang, Yiheng Wu, Yipeng Jing, Houjun Mo, David H. Weinberg

Abstract: We investigate the level of galaxy assembly bias in the Sloan Digital Sky Survey (SDSS) main galaxy sample using ELUCID, a state-of-the-art constrained simulation that accurately reconstructed the initial density perturbations within the SDSS volume. On top of the ELUCID haloes, we develop an extended HOD model that includes the assembly bias of central and satellite galaxies, parameterized as… ▽ More We investigate the level of galaxy assembly bias in the Sloan Digital Sky Survey (SDSS) main galaxy sample using ELUCID, a state-of-the-art constrained simulation that accurately reconstructed the initial density perturbations within the SDSS volume. On top of the ELUCID haloes, we develop an extended HOD model that includes the assembly bias of central and satellite galaxies, parameterized as $\mathcal{Q}_\mathrm{cen}$ and $\mathcal{Q}_\mathrm{sat}$, respectively, to predict a suite of one- and two-point observables. In particular, our fiducial constraint employs the probability distribution of the galaxy number counts measured on $8\,\mathrm{Mpc}\,h^{-1}$ scales $N_8^g$ and the projected cross-correlation functions of quintiles of galaxies selected by $N_8^g$ with our entire galaxy sample. We perform extensive tests of the efficacy of our method by fitting the same observables to mock data using both constrained and non-constrained simulations. We discover that in many cases the level of cosmic variance between the two simulations can produce biased constraints that lead to an erroneous detection of galaxy assembly bias if the non-constrained simulation is used. When applying our method to the SDSS data, the ELUCID reconstruction effectively removes an otherwise strong degeneracy between cosmic variance and galaxy assembly bias in SDSS, enabling us to derive an accurate and stringent constraint on the latter. Our fiducial ELUCID constraint, for galaxies above a stellar mass threshold $M_*{=}10^{10.2}\,h^{-2}\,M_\odot$, is $\mathcal{Q}_\mathrm{cen}{=}{-}0.09\pm{0.05}$ and $\mathcal{Q}_\mathrm{sat}{=}0.09\pm{0.10}$, indicating no evidence for a significant~($>2σ$) galaxy assembly bias in the local Universe probed by SDSS. Finally, our method provides a promising path to the robust modelling of the galaxy-halo connection within future surveys like DESI and PFS. △ Less

Submitted 7 September, 2022; v1 submitted 8 October, 2020; originally announced October 2020.

Comments: 25 pages, 7 main figures, and 2 appendix figures; accepted by SCPMA (https://www.sciengine.com/SCPMA/doi/10.1007/s11433-022-1955-7); highlight review (https://www.sciengine.com/SCPMA/doi/10.1007/s11433-022-1967-6)

Journal ref: Sci. China Phys. Mech. Astron. 65, 109811 (2022)

How to empirically model star formation in dark matter halos: I. Inferences about central galaxies from numerical simulations

Authors: Yangyao Chen, H. J. Mo, Cheng Li, Kai Wang

Abstract: 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 h… ▽ More 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. △ Less

Submitted 3 March, 2021; v1 submitted 25 September, 2020; originally announced September 2020.

Comments: 20 pages, 17 figures, 2 tables. MNRAS accepted. Comments are welcomed

Measuring galaxy abundance and clustering at high redshift from incomplete spectroscopic data: Tests on mock catalogs

Authors: Jiacheng Meng, Cheng Li, Houjun Mo, Yangyao Chen, Kai Wang

Abstract: The number density and correlation function of galaxies are two key quantities to characterize the distribution of the observed galaxy population. High-$z$ spectroscopic surveys, which usually involve complex target selection and are incomplete in redshift sampling, present both opportunities and challenges to measure these quantities reliably in the high-$z$ Universe. Using realistic mock catalog… ▽ More The number density and correlation function of galaxies are two key quantities to characterize the distribution of the observed galaxy population. High-$z$ spectroscopic surveys, which usually involve complex target selection and are incomplete in redshift sampling, present both opportunities and challenges to measure these quantities reliably in the high-$z$ Universe. Using realistic mock catalogs we show that target selection and redshift incompleteness can lead to significantly biased results, especially due to the flux limit selection criteria. We develop a new method to correct the flux limit effect, using information provided by the parent photometric data from which the spectroscopic sample is constructed. Our tests using realistic mock samples show that the method is able to reproduce the true stellar mass function and correlation function reliably. Mock catalogs are constructed for the existing zCOSMOS and VIPERS surveys, as well as for the forthcoming PFS galaxy evolution survey. The same set of mock samples are used to quantify the total variance expected for different sample sizes. We find that the total variance decreases very slowly when the survey area reaches about 4 deg$^2$ for the abundance and about 8 deg$^2$ for the clustering, indicating that the cosmic variance is no longer the dominant source of error for PFS-like surveys. We also quantify improvements expected in the PFS-like galaxy survey relative to zCOSMOS and VIPERS surveys. △ Less

Submitted 2 February, 2024; v1 submitted 31 August, 2020; originally announced August 2020.

Comments: The paper is significantly revised and the flux limit effect is emphasized. The method to correct the flux limit effect on clustering measurement is also improved. 24 pages, 15 figures, and this version is accepted by ApJ

Properties of the ionized CGM and IGM: a comparison of galaxy formation models via the Sunyaev-Zel'dovich effect

Authors: S. H. Lim, David Barnes, Mark Vogelsberger, H. J. Mo, Dylan Nelson, Annalisa Pillepich, Klaus Dolag, Federico Marinacci

Abstract: 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 t… ▽ More 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. △ Less

Submitted 16 April, 2021; v1 submitted 22 July, 2020; originally announced July 2020.

Comments: 15 pages, 6 figures, accepted for publication in MNRAS

Constraints on circum-galactic media from Sunyaev-Zel'dovich effects and X-ray data

Authors: Xuanyi Wu, Houjun Mo, Cheng Li, Seunghwan Lim

Abstract: We use observational measurements of thermal and kinetic Sunyaev-Zel'dovich effects, as well as soft X-ray emission associated with galaxy groups to constrain the gas density and temperature in the circumgalactic medium (CGM) for dark matter halos with masses above $10^{12.5}M_{\odot}$. A number of generic models are used together with a Bayesian scheme to make model inferences. We find that gas w… ▽ More We use observational measurements of thermal and kinetic Sunyaev-Zel'dovich effects, as well as soft X-ray emission associated with galaxy groups to constrain the gas density and temperature in the circumgalactic medium (CGM) for dark matter halos with masses above $10^{12.5}M_{\odot}$. A number of generic models are used together with a Bayesian scheme to make model inferences. We find that gas with a single temperature component cannot provide a consistent model to match the observational data. A simple two-phase model assuming a hot component and an ionized warm component can accommodate all the three observations. The total amount of the gas in individual halos is inferred to be comparable to the universal baryon fraction corresponding to the halo mass. The inferred temperature of the hot component is comparable to the halo virial temperature. The fraction of the hot component increases from $(15-40)\%$ for $10^{12.5}{M}_\odot$ halos to $(40-60)\%$ for $10^{14.5}{M}_\odot$ halos, where the ranges reflect uncertainties in the assumed gas density profile. Our results suggest that a significant fraction of the halo gas is in a non-thermalized component with temperature much lower than the virial temperature. △ Less

Submitted 13 September, 2020; v1 submitted 10 July, 2020; originally announced July 2020.

Comments: 10pages, 3 figures

Identifying galaxy groups at high redshift from incomplete spectroscopic data: I. The group finder and application to zCOSMOS

Authors: Kai Wang, H. J. Mo, Cheng Li, Jiacheng Meng, Yangyao Chen

Abstract: 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 sa… ▽ More 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. △ Less

Submitted 12 September, 2020; v1 submitted 9 June, 2020; originally announced June 2020.

Comments: 18 pages, 12 figures, 6 table. MNRAS accepted. Comments are welcomed

Characteristic mass in galaxy quenching: environmental versus internal effects

Authors: Pengfei Li, Huiyuan Wang, H. J. Mo, Enci Wang, Hui Hong

Abstract: 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… ▽ More 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'. △ Less

Submitted 12 September, 2020; v1 submitted 21 March, 2020; originally announced March 2020.

Comments: 11 pages, 7 figures, accepted by ApJ

SDSS-IV MaNGA: Bayesian analysis of the star formation history of low-mass galaxies in the local Universe

Authors: Shuang Zhou, H. J. Mo, Cheng Li, Médéric Boquien, Graziano Rossi

Abstract: 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… ▽ More 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. △ Less

Submitted 30 September, 2020; v1 submitted 19 March, 2020; originally announced March 2020.

Comments: 21 pages, 22 figures, MNRAS accepted version

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

Authors: Yangyao Chen, H. J. Mo, Cheng Li, Huiyuan Wang, Xiaohu Yang, Youcai Zhang, Kai Wang

Abstract: 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 correla… ▽ More 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. △ Less

Submitted 28 July, 2020; v1 submitted 11 March, 2020; originally announced March 2020.

Comments: 22 pages, 11 figures, 2 tables. ApJ accepted. Comments are welcomed

The Formation History of Subhalos and the Evolution of Satellite Galaxies

Authors: Jingjing Shi, Huiyuan Wang, Houjun Mo, Mark Vogelsberger, Luis C. Ho, Min Du, Dylan Nelson, Annalisa Pillepich, Lars Hernquist

Abstract: Satellites constitute an important fraction of the overall galaxy population and are believed to form in dark matter subhalos. Here we use the cosmological hydrodynamic simulation TNG100 to investigate how the formation histories of subhalos affect the properties and evolution of their host galaxies. We use a scaled formation time ($a_{\rm nf}$) to characterize the mass assembly histories of the s… ▽ More Satellites constitute an important fraction of the overall galaxy population and are believed to form in dark matter subhalos. Here we use the cosmological hydrodynamic simulation TNG100 to investigate how the formation histories of subhalos affect the properties and evolution of their host galaxies. We use a scaled formation time ($a_{\rm nf}$) to characterize the mass assembly histories of the subhalos before they are accreted by massive host halos. We find that satellite galaxies in young subhalos (low $a_{\rm nf}$) are less massive and more gas rich, and have stronger star formation and a higher fraction of ex situ stellar mass than satellites in old subhalos (high $a_{\rm nf}$). Furthermore, these low $a_{\rm nf}$ satellites require longer timescales to be quenched as a population than the high $a_{\rm nf}$ counterparts. We find very different merger histories between satellites in fast accretion (FA, $a_{\rm nf}<1.3$) and slow accretion (SA, $a_{\rm nf}>1.3$) subhalos. For FA satellites, the galaxy merger frequency dramatically increases just after accretion, which enhances the star formation at accretion. While, for SA satellites, the mergers occur smoothly and continuously across the accretion time. Moreover, mergers with FA satellites happen mainly after accretion, while a contrary trend is found for SA satellites. Our results provide insight into the evolution and star formation quenching of the satellite population. △ Less

Submitted 28 March, 2020; v1 submitted 13 January, 2020; originally announced January 2020.

Comments: 21 pages, 10 figures, accepted for publication in ApJ

Estimating dust attenuation from galactic spectra. I. methodology and tests

Authors: Niu Li, Cheng Li, Houjun Mo, Jian Hu, Shuang Zhou, Cheng Du

Abstract: We develop a method to estimate the dust attenuation curve of galaxies from full spectral fitting of their optical spectra. Motivated from previous studies, we separate the small-scale features from the large-scale spectral shape, by performing a moving average method to both the observed spectrum and the simple stellar population model spectra. The intrinsic dust-free model spectrum is then deriv… ▽ More We develop a method to estimate the dust attenuation curve of galaxies from full spectral fitting of their optical spectra. Motivated from previous studies, we separate the small-scale features from the large-scale spectral shape, by performing a moving average method to both the observed spectrum and the simple stellar population model spectra. The intrinsic dust-free model spectrum is then derived by fitting the observed ratio of the small-scale to large-scale (S/L) components with the S/L ratios of the SSP models. The selective dust attenuation curve is then determined by comparing the observed spectrum with the dust-free model spectrum. One important advantage of this method is that the estimated dust attenuation curve is independent of the shape of theoretical dust attenuation curves. We have done a series of tests on a set of mock spectra covering wide ranges of stellar age and metallicity. We show that our method is able to recover the input dust attenuation curve accurately, although the accuracy depends slightly on signal-to-noise ratio of the spectra. We have applied our method to a number of edge-on galaxies with obvious dust lanes from the ongoing MaNGA survey, deriving their dust attenuation curves and $E(B-V)$ maps, as well as dust-free images in $g$, $r$, and $i$ bands. These galaxies show obvious dust lane features in their original images, which largely disappear after we have corrected the effect of dust attenuation. The vertical brightness profiles of these galaxies become axis-symmetric and can well be fitted by a simple model proposed for the disk vertical structure. Comparing the estimated dust attenuation curve with the three commonly-adopted model curves, we find that the Calzetti curve provides the best description of the estimated curves for the inner region of galaxies, while the Milky Way and SMC curves work better for the outer region. △ Less

Submitted 11 May, 2020; v1 submitted 8 January, 2020; originally announced January 2020.

Comments: 27 pages, 16 figures, Accepted for publication in ApJ