Host Halo Mass Research Articles

Constraining the Low-mass End of the Stellar-to-halo Mass Relation with Surveys of Satellite Galaxies

The abundance of satellite galaxies is set by the hierarchical assembly of their host halo. We leverage this to investigate the low-mass end (M H < 1011 M ⊙) of the stellar-to-halo mass relation (SHMR), which is key to constraining theories of galaxy formation and cosmology. We argue that recent analyses of satellite galaxies in the Local Group environment have not adequately modeled the dominant source of scatter in satellite stellar mass functions: the variance in accretion histories for a fixed host halo mass. We present a novel inference framework that not only properly accounts for this halo-to-halo variance but also naturally identifies the amount of host halo mass mixing, which is generally unknown. Specifically, we use the semianalytical SatGen model to construct mock satellite galaxy populations consistent with the third data release of the Satellites Around Galactic Analogs survey. We demonstrate that even under the most idealized circumstances, the halo-to-halo variance makes it virtually impossible to put any meaningful constraints on the scatter in the SHMR. Even a satellite galaxy survey made up 100 hosts can at best only place an upper limit of ∼0.5 dex on the scatter (at the 95% confidence level). This is because the large variance in halo assembly histories dominates over the scatter in the SHMR. This problem can be overcome by increasing the sample size of the survey by an order of magnitude (∼1000 host galaxies), something that should be fairly straightforward with forthcoming spectroscopic surveys.

The VST ATLAS Quasar Survey – III. Halo mass function via quasar clustering and quasar-CMB lensing cross-clustering

ABSTRACT We exploit the VST ATLAS quasar (QSO) catalogue to perform three measurements of the quasar halo mass profile. First, we make a new estimate of the angular autocorrelation function of ≈230 000 ATLAS quasars with $z_{\rm photo}\lesssim 2.5$ and $17 &amp;lt; g &amp;lt; 22$. By comparing with the $\Lambda$CDM mass clustering correlation function, we measure the quasar bias to be $b_{\rm Q}\approx 2.1$, implying a quasar halo mass of $M_{\rm halo} \approx 8.5\times 10^{11}\,h^{-1}\,{\rm M}_\odot$. Second, we cross-correlate these $z\approx 1.7$ ATLAS quasars with the Planck cosmic microwave background (CMB) lensing maps, detecting a somewhat stronger signal at $4\,{\rm arcmin} &amp;lt; \theta &amp;lt; 60\,{\rm arcmin}$ than previous authors. Scaling these authors’ model fit to our data, we estimate a quasar host halo mass of $M_{\rm halo}\approx 8.3\times 10^{11}\,h^{-1}\,{\rm M}_\odot$. Third, we fit halo occupation sistribution (HOD) model parameters to our quasar autocorrelation function and from the derived halo mass function, we estimate a quasar halo mass of $M_{\rm halo}\approx 2.5\times 10^{12}\,h^{-1}\,{\rm M}_\odot$. We then compare our HOD model prediction to our quasar-CMB lensing result, confirming their consistency. We find that most (≈2/3) QSOs have halo masses within a factor of ≈3 of this average mass. An analysis based on the probability of X-ray detections of AGN in galaxies and the galaxy stellar mass function gives a similarly small mass range. Finally, we compare the quasar halo mass and luminosity functions and suggest that gravitational growth may produce the constant space density with redshift seen in the quasar luminosity function.

The Mass and Redshift Dependence of Halo Star Clustering

We adopt the two-point correlation function (2PCF) as a statistical tool to quantify the spatial clustering of halo stars for galaxy systems spanning a wide range in host halo virial mass ( 11.25<log10M200c/M⊙<15 ) and redshifts (0 < z < 1.5) from the IllustrisTNG simulations. Consistent with a previous study, we identify clear correlations between the strength of the 2PCF signals and galaxy formation redshifts, but over a much wider mass range. We find that such correlations are slightly stronger at higher redshifts and get weakened with the increase of host halo mass. We demonstrate that the spatial clustering of halo stars is affected by two factors: (1) the clustering gets gradually weakened as time passes (phase mixing); (2) newly accreted stars at more recent times would increase the clustering. For more massive galaxy systems, they assemble late, and the newly accreted stars would increase the clustering. The late assembly of massive systems may also help to explain the weaker correlations between the 2PCF signals and the galaxy formation redshifts in massive halos, as their 2PCFs are affected more by recently accreted stars, while formation redshift characterizes mass accretion on a much longer timescale. We find that the orbits of satellite galaxies in more massive halos maintain larger radial anisotropy, reflecting the more active accretion state of their hosts while also contributing to their stronger mass loss rates.

The SAGA Survey. IV. The Star Formation Properties of 101 Satellite Systems around Milky Way–mass Galaxies

We present the star-forming properties of 378 satellite galaxies around 101 Milky Way analogs in the Satellites Around Galactic Analogs (SAGA) Survey, focusing on the environmental processes that suppress or quench star formation. In the SAGA stellar mass range of 106−10 M ⊙, we present quenched fractions, star-forming rates, gas-phase metallicities, and gas content. The fraction of SAGA satellites that are quenched increases with decreasing stellar mass and shows significant system-to-system scatter. SAGA satellite quenched fractions are highest in the central 100 kpc of their hosts and decline out to the virial radius. Splitting by specific star formation rate (sSFR), the least star-forming satellite quartile follows the radial trend of the quenched population. The median sSFR of star-forming satellites increases with decreasing stellar mass and is roughly constant with projected radius. Star-forming SAGA satellites are consistent with the star formation rate–stellar mass relationship determined in the Local Volume, while the median gas-phase metallicity is higher and median H i gas mass is lower at all stellar masses. We investigate the dependence of the satellite quenched fraction on host properties. Quenched fractions are higher in systems with larger host halo mass, but this trend is only seen in the inner 100 kpc; we do not see significant trends with host color or star formation rate. Our results suggest that lower-mass satellites and satellites inside 100 kpc are more efficiently quenched in a Milky Way–like environment, with these processes acting sufficiently slowly to preserve a population of star-forming satellites at all stellar masses and projected radii.

Weak Gravitational Lensing around Low Surface Brightness Galaxies in the DES Year 3 Data

We present galaxy-galaxy lensing measurements using a sample of low surface brightness galaxies (LSBGs) drawn from the Dark Energy Survey Year 3 (Y3) data as lenses. LSBGs are diffuse galaxies with a surface brightness dimmer than the ambient night sky. These dark-matter-dominated objects are intriguing due to potentially unusual formation channels that lead to their diffuse stellar component. Given the faintness of LSBGs, using standard observational techniques to characterize their total masses proves challenging. Weak gravitational lensing, which is less sensitive to the stellar component of galaxies, could be a promising avenue to estimate the masses of LSBGs. Our LSBG sample consists of 23,790 galaxies separated into red and blue color types at g−i≥0.60 and g−i&lt;0.60, respectively. Combined with the DES Y3 shear catalog, we measure the tangential shear around these LSBGs and find signal-to-noise ratios of 6.67 for the red sample, 2.17 for the blue sample, and 5.30 for the full sample. We use the clustering redshifts method to obtain redshift distributions for the red and blue LSBG samples. Assuming all red LSBGs are satellites, we fit a simple model to the measurements and estimate the host halo mass of these LSBGs to be . We place a 95% upper bound on the subhalo mass at . By contrast, we assume the blue LSBGs are centrals, and place a 95% upper bound on the halo mass at log(Mhost/M⊙)&lt;11.84. We find that the stellar-to-halo mass ratio of the LSBG samples is consistent with that of the general galaxy population. This work illustrates the viability of using weak gravitational lensing to constrain the halo masses of LSBGs.

A unified model for the clustering of quasars and galaxies at z ≈ 6

ABSTRACT Recent observations from the EIGER JWST program have measured for the first time the quasar–galaxy cross-correlation function at $z\approx 6$. The autocorrelation function of faint $z\approx 6$ quasars was also recently estimated. These measurements provide key insights into the properties of quasars and galaxies at high redshift and their relation with the host dark matter haloes. In this work, we interpret these data building upon an empirical quasar population model that has been applied successfully to quasar clustering and demographic measurements at $z\approx 2\!-\!4$. We use a new, large-volume N-body simulation with more than a trillion particles, FLAMINGO-10k, to model quasars and galaxies simultaneously. We successfully reproduce observations of $z\approx 6$ quasars and galaxies (i.e. their clustering properties and luminosity functions), and infer key quantities such as their luminosity–halo mass relation, the mass function of their host haloes, and their duty cycle/occupation fraction. Our key findings are (i) quasars reside on average in $\approx 10^{12.5}\, {\rm M}_{\odot }$ haloes (corresponding to $\approx 5\sigma$ fluctuations in the initial conditions of the linear density field), but the distribution of host halo masses is quite broad; (ii) the duty cycle of (UV-bright) quasar activity is relatively low ($\approx 1~{{\ \rm per\ cent}}$); (iii) galaxies (that are bright in [O iii]) live in much smaller haloes ($\approx 10^{10.9}\, {\rm M}_{\odot }$) and have a larger duty cycle (occupation fraction) of $\approx 13~{{\ \rm per\ cent}}$. Finally, we focus on the inferred properties of quasars and present a homogeneous analysis of their evolution with redshift. The picture that emerges reveals a strong evolution of the host halo mass and duty cycle of quasars at $z\approx 2\!-\!6$, and calls for new investigations of the role of quasar activity across cosmic time.

EIGER. VI. The Correlation Function, Host Halo Mass, and Duty Cycle of Luminous Quasars at z ≳ 6

We expect luminous (M 1450 ≲ −26.5) high-redshift quasars to trace the highest-density peaks in the early Universe. Here, we present observations of four z ≳ 6 quasar fields using JWST/NIRCam in the imaging and wide-field slitless spectroscopy mode and report a wide range in the number of detected [O iii]-emitting galaxies in the quasars’ environments, ranging between a density enhancement of δ ≈ 65 within a 2 cMpc radius—one of the largest protoclusters during the Epoch of Reionization discovered to date—to a density contrast consistent with zero, indicating the presence of a UV-luminous quasar in a region comparable to the average density of the Universe. By measuring the two-point cross-correlation function of quasars and their surrounding galaxies, as well as the galaxy autocorrelation function, we infer a correlation length of quasars at 〈z〉 = 6.25 of r0QQ=22.0−2.9+3.0cMpch−1 , while we obtain a correlation length of the [O iii]-emitting galaxies of r0GG=4.1±0.3cMpch−1 . By comparing the correlation functions to dark-matter-only simulations we estimate the minimum mass of the quasars’ host dark matter halos to be log10(Mhalo,min/M⊙)=12.43−0.15+0.13 (and log10(Mhalo,min[OIII]/M⊙)=10.56−0.03+0.05 for the [O iii] emitters), indicating that (a) luminous quasars do not necessarily reside within the most overdense regions in the early Universe, and that (b) the UV-luminous duty cycle of quasar activity at these redshifts is f duty ≪ 1. Such short quasar activity timescales challenge our understanding of early supermassive black hole growth and provide evidence for highly dust-obscured growth phases or episodic, radiatively inefficient accretion rates.

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

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 h) for the centrals has often been used to indicate the halo assembly time t h,50 of the group, where t 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 h shows a significant scatter with t 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 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 h,50, within an accuracy of ∼1.09 Gyr. With careful calibrations of individual observable quantities in the mocks with Sloan Digital Sky Survey (SDSS) observations, we apply the trained models to the SDSS Yang et al. groups and derive the t h,50 for each group for the first time. The derived SDSS t 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.

The Massive and Distant Clusters of WISE Survey 2: A Stacking Analysis Investigating the Evolution of Star Formation Rates and Stellar Masses in Groups and Clusters

The evolution of galaxies depends on their masses and local environments; understanding when and how environmental quenching starts to operate remains a challenge. Furthermore, studies of the high-redshift regime have been limited to massive cluster members, owing to sensitivity limits or small fields of view when the sensitivity is sufficient, intrinsically biasing the picture of cluster evolution. In this work, we use stacking to investigate the average star formation history of more than 10,000 groups and clusters drawn from the Massive and Distant Clusters of WISE Survey 2. Our analysis covers near-ultraviolet to far-infrared wavelengths, for galaxy overdensities at 0.5 ≲ z ≲ 2.54. We employ spectral energy distribution fitting to measure the specific star formation rates (sSFRs) in four annular apertures with radii between 0 and 1000 kpc. At z ≳ 1.6, the average sSFR evolves similarly to the field in both the core and the cluster outskirts. Between z¯=1.60 and z¯=1.35 , the sSFR in the core drops sharply, and it continues to fall relative to the field sSFR at lower redshifts. We interpret this change as evidence that the impact of environmental quenching dramatically increases at z ∼ 1.5, with the short time span of the transition suggesting that the environmental quenching mechanism dominant at this redshift operates on a rapid timescale. We find indications that the sSFR may decrease with increasing host halo mass, but lower-scatter mass tracers than the signal-to-noise ratio are needed to confirm this relationship.

Impact of astrophysical scatter on the epoch of reionization [H i]21 bispectrum

It is believed that the first star-forming galaxies are the main drivers of cosmic reionization. It is usually assumed that there is a one-to-one relationship between the star formation rate (SFR) inside a galaxy and the host halo mass in semi-analytical/numerical modeling of large-scale ionization maps during the epoch of reionization. However, more accurate simulations and observations suggest that the SFR and ionizing luminosity in galaxies may vary considerably even if the host halo mass is the same. This astrophysical scatter can introduce an additional non-Gaussianity in the [H i]21cm signal, which might not be captured adequately in the power spectrum. In this work, we have studied the impact of the scatter on the [H i]21cm bispectrum using semi-numerical simulations. We find that the scatter primarily affects small ionized regions, whereas the large ionized bubbles remain largely unaffected. Although, the fractional change in the [H i]21cm bispectra due to the scatter is found to be more than a factor of 10 at large scales (k 1 ≲ 1 Mpc-1) for z=7.4, it is found to be statistically insignificant. However, at small scales (k 1 ~ 2.55 Mpc-1), we have found the impact due to the scatter to be high in magnitude (|〈Δ B 〉/B no-scatter| ~ 1) and statistically significant (|〈Δ B〉/σ ΔB| ≳ 5) at neutral fraction, x̅HI ~ 0.8 for z=7.4. The impact due to scatter is found to be even more prominent (|〈Δ B 〉/B no-scatter| ≳ 10) at small scales for z=10 and x̅HI ~ 0.8, but with reduced statistical significance to some extent (|〈Δ B〉/σ ΔB| ~ 3), compared to z=7.4 at the same neutral fraction. We have also found that in the most optimistic scenario, SKA1-Low might be able to detect these signatures of astrophysical scatter, at ~ 3σ and ~ 5σ detection significance for x̅HI ~ 0.8 and 0.9 respectively, for the equilateral [H i]21cm bispectrum at z=7.4.

The beyond-halo mass effects of the cosmic web environment on galaxies

ABSTRACT Galaxy properties primarily depend on their host halo mass. Halo mass, in turn, depends on the cosmic web environment. We explore if the effect of the cosmic web on galaxy properties is entirely transitive via host halo mass, or if the cosmic web has an effect independent of mass. The secondary galaxy bias, sometimes referred to as ‘galaxy assembly bias’, is the beyond-mass component of the galaxy–halo connection. We investigate the link between the cosmic web environment and the secondary galaxy bias in simulations. We measure the secondary galaxy bias through the following summary statistics: projected two-point correlation function, $w_{\mathrm{p}}(r_{\mathrm{p}})$, and counts-in-cylinders statistics, $P(N_{\mathrm{CIC}})$. First, we examine the extent to which the secondary galaxy bias can be accounted for with a measure of the environment as a secondary halo property. We find that the total secondary galaxy bias preferentially places galaxies in more strongly clustered haloes. In particular, haloes at fixed mass tend to host more galaxies when they are more strongly associated with nodes or filaments. This tendency accounts for a significant portion, but not the entirety, of the total secondary galaxy bias effect. Secondly, we quantify how the secondary galaxy bias behaves differently depending on the host halo proximity to nodes and filaments. We find that the total secondary galaxy bias is relatively stronger in haloes more associated with nodes or filaments. We emphasize the importance of removing halo mass effects when considering the cosmic web environment as a factor in the galaxy–halo connection.

LIGHTS. Survey Overview and a Search for Low Surface Brightness Satellite Galaxies

We present an overview of the LBT Imaging of Galactic Halos and Tidal Structures survey, which currently includes 25 nearby galaxies that are on average ∼1 mag fainter than the Milky Way, and a catalog of 54 low central surface brightness (24 < μ 0,g /mag arcsec−2 < 28) satellite galaxy candidates, most of which were previously uncatalogued. The depth of the imaging exceeds the full 10 yr depth of the Rubin Observatory’s Legacy Survey of Space and Time. We find, after applying completeness corrections, rising numbers of candidate satellites as we approach the limiting luminosity (M r ∼ −8 mag) and central surface brightness (μ 0,g ∼ 28 mag arcsec−2). Over the parameter range we explore, each host galaxy (excluding those that are in overdense regions, apparently groups) has nearly four such candidate satellites to a projected radius of ∼100 kpc. These objects are mostly just at or beyond the reach of spectroscopy unless they are H i rich or have ongoing star formation. We identify three, possibly four, ultra-diffuse satellite galaxies (effective radius >1.5 kpc). This incidence rate falls within expectations of the extrapolation of the published relationship between the number of ultra-diffuse satellite galaxies and host halo mass. Last, we visually identify 12 candidate satellites that host a nuclear star cluster (NSC). The NSC occupation fraction for the sample (12/54) matches that published for satellites of early-type galaxies, suggesting that the parent’s morphological type plays at most a limited role in determining the NSC occupation fraction.

Devouring the Centaurus A Satellites: Modeling Dwarf Galaxies with Galacticus

For the first time, systematic studies of dwarf galaxies are being conducted throughout the Local Volume, including the dwarf satellites of the nearby giant elliptical galaxy Centaurus A (NGC 5128). Given Centaurus A's mass (roughly 10 times larger than that of the Milky Way), AGN activity, and recent major mergers, investigating the dwarf galaxies of Centaurus A and their star formation physics is imperative. However, simulating the faintest dwarfs around a galaxy of Centaurus A's mass with sufficient resolution in a hydrodynamic simulation is computationally expensive and currently infeasible. In this study, we seek to reproduce the properties of Centaurus A dwarfs using the semianalytic model Galacticus to model dwarfs within a 700 kpc region around Centaurus A, corresponding approximately to its splashback radius. We investigate the effects of host halo mass and environment and predict observable properties of Centaurus A dwarfs using astrophysical prescriptions and parameters previously tuned to match properties of the Milky Way’s satellite galaxies. This approach allows us to approximately replicate cumulative luminosity functions, and luminosity–metallicity and luminosity–half-light-radii relations observed in the Centaurus A satellites. We provide predictions for the velocity dispersions, and star formation histories of Centaurus A dwarfs. The agreement between our predicted star formation histories for Centaurus A dwarfs and those of the Milky Way dwarfs implies the presence of universal processes governing star formation in dwarf galaxies. Overall, our findings shed light on the star formation physics of dwarf galaxies in the Centaurus A system, revealing insights into their properties and dependence on the host environment.

The host dark matter haloes of the first quasars

ABSTRACT If z &amp;gt; 6 quasars reside in rare, massive haloes, ΛCDM cosmology predicts they should be surrounded by an anomalously high number of bright companion galaxies. Here, I show that these companion galaxies should also move unusually fast. Using a new suite of cosmological, ‘zoom-in’ hydrodynamic simulations, I present predictions for the velocity distribution of quasar companion galaxies and its variation with quasar host halo mass at $z \, = \, 6$. Satellites accelerate as they approach the quasar host galaxy, producing a line-of-sight velocity profile that broadens with decreasing distance to the quasar host galaxy. This increase in velocity dispersion is particularly pronounced if the host halo mass is $\gtrsim 5 \times 10^{12} \, \rm M_\odot$. Typical line-of-sight speeds rise to $\approx 500 \, \rm km \, s^{-1}$ at projected radii $\sim 10 \, \rm kpc$. For about 10 per cent of satellites, they should exceed $800 \, \rm km \, s^{-1}$, with $\approx 5~{{\ \rm per\ cent}}$ of companions reaching line-of-sight speeds $\sim 1000 \, \rm km \, s^{-1}$. For lower host halo masses $\approx 5 \times 10^{11}{-}10^{12} \, \rm M_\odot$, the velocity profile of companion galaxies is significantly flatter. In this case, typical line-of-sight velocities are $\approx 250 \, \rm km \, s^{-1}$ and do not exceed $\approx 500 \, \rm km \, s^{-1}$. A comparison with existing Atacama Large Millimeter/submillimeter Array (ALMA), JWST and Multi-Unit Spectroscopic Explorer (MUSE) line-of-sight velocity measurements reveals that observed z &amp;gt; 6 quasar companions closely follow the velocity distribution expected for a host halo with mass $\gtrsim 5 \times 10^{12} \, \rm M_\odot$, ruling out a light host halo. Finally, through an estimate of ultraviolet and [$\rm O \, {\small III}$] luminosity functions, I show that the velocity distribution more reliably discriminates between halo mass than companion number counts, which are strongly affected by cosmic variance.

Zero Metallicity with Zero CPU Hours: Masses of the First Stars on the Laptop

We develop an analytic model for the mass of the first stars forming in the centers of primordial gas clouds as a function of host halo mass, redshift, and degree of rotation. The model is based on the estimation of key timescales determining the following three processes: the collapse of the gas cloud, the accretion onto the protostellar core, and the radiative feedback of the protostellar core. The final stellar mass is determined by the total mass accreted until the radiative feedback halts the accretion. The analytic estimation, motivated by the result of the full numerical simulations, leads to algebraic expressions allowing an extremely fast execution. Despite its simplicity, the model reproduces the stellar mass scale and its parameter dependencies observed in state-of-the-art cosmological zoom-in simulations. This work clarifies the basic physical principles undergirding such numerical treatments and provides a path to efficiently calibrating numerical predictions against eventual observations of the first stars.

Revisiting the extreme clustering of z ≈ 4 quasars with large volume cosmological simulations

ABSTRACT Observations from wide-field quasar surveys indicate that the quasar autocorrelation length increases dramatically from z ≈ 2.5 to ≈ 4. This large clustering amplitude at z ≈ 4 has proven hard to interpret theoretically, as it implies that quasars are hosted by the most massive dark matter haloes residing in the most extreme environments at that redshift. In this work, we present a model that simultaneously reproduces both the observed quasar autocorrelation and quasar luminosity functions. The spatial distribution of haloes and their relative abundance are obtained via a novel method that computes the halo mass and halo cross-correlation functions by combining multiple large-volume dark-matter-only cosmological simulations with different box sizes and resolutions. Armed with these halo properties, our model exploits the conditional luminosity function framework to describe the stochastic relationship between quasar luminosity, L, and halo mass, M. Assuming a simple power-law relation L ∝ Mγ with lognormal scatter, σ, we are able to reproduce observations at z ∼ 4 and find that: (i) the quasar luminosity–halo mass relation is highly non-linear (γ ≳ 2), with very little scatter (σ ≲ 0.3 dex); (ii) luminous quasars ($\log _{10} L/{\rm erg}\, {\rm s}^{-1}\gtrsim 46.5-47$) are hosted by haloes with mass log10M/M⊙ ≳ 13–13.5; and (iii) the implied duty cycle for quasar activity approaches unity ($\varepsilon _{\rm DC}\approx 10\,\mathrm{ per}\,\mathrm{ cent}-60~{{\ \rm per\ cent}}$). We also consider observations at z ≈ 2.5 and find that the quasar luminosity–halo mass relation evolves significantly with cosmic time, implying a rapid change in quasar host halo masses and duty cycles, which in turn suggests concurrent evolution in black hole scaling relations and/or accretion efficiency.

Modelling the orbital histories of satellites of Milky Way-mass galaxies: testing static host potentials against cosmological simulations

ABSTRACT Understanding the evolution of satellite galaxies of the Milky Way (MW) and M31 requires modelling their orbital histories across cosmic time. Many works that model satellite orbits incorrectly assume or approximate that the host halo gravitational potential is fixed in time and is spherically symmetric or axisymmetric. We rigorously benchmark the accuracy of such models against the FIRE-2 cosmological baryonic simulations of MW/M31-mass haloes. When a typical surviving satellite fell in ($3.4\!-\!9.7\, \rm {Gyr}$ ago), the host halo mass and radius were typically 26–86 per cent of their values today, respectively. Most of this mass growth of the host occurred at small distances, $r\lesssim 50\, \rm {kpc}$, opposite to dark matter only simulations, which experience almost no growth at small radii. We fit a near-exact axisymmetric gravitational potential to each host at z = 0 and backward integrate the orbits of satellites in this static potential, comparing against the true orbit histories in the simulations. Orbital energy and angular momentum are not well conserved throughout an orbital history, varying by 25 per cent from their current values already $1.6\!-\!4.7\, \rm {Gyr}$ ago. Most orbital properties are minimally biased, ≲10 per cent, when averaged across the satellite population as a whole. However, for a single satellite, the uncertainties are large: recent orbital properties, like the most recent pericentre distance, typically are ≈20 per cent uncertain, while earlier events, like the minimum pericentre or the infall time, are ≈40–80 per cent uncertain. Furthermore, these biases and uncertainties are lower limits, given that we use near-exact host mass profiles at z = 0.

The DESI One-Percent Survey: exploring a generalized SHAM for multiple tracers with the UNIT simulation

ABSTRACT We perform SubHalo Abundance Matching (SHAM) studies on UNIT simulations with {σ, Vceil, vsmear}-SHAM and {σ, Vceil, fsat}-SHAM. They are designed to reproduce the clustering on 5–30 $\, {\, h^{-1}\, {\rm Mpc}}$ of luminous red galaxies (LRGs), emission-line galaxies (ELGs), and quasi-stellar objects (QSOs) at 0.4 &amp;lt; z &amp;lt; 3.5 from DESI (Dark Energy Spectroscopic Instrument) One Percent Survey. Vceil is the incompleteness of the massive host (sub)haloes and is the key to the generalized SHAM. vsmear models the clustering effect of redshift uncertainties, providing measurements consistent with those from repeat observations. A free satellite fraction fsat is necessary to reproduce the clustering of ELGs. We find ELGs present a more complex galaxy–halo mass relation than LRGs reflected in their weak constraints on σ. LRGs, QSOs, and ELGs show increasing Vceil values, corresponding to the massive galaxy incompleteness of LRGs, the quenched star formation of ELGs and the quenched black hole accretion of QSOs. For LRGs, a Gaussian vsmear presents a better profile for subsamples at redshift bins than a Lorentzian profile used for other tracers. The impact of the statistical redshift uncertainty on ELG clustering is negligible. The best-fitting satellite fraction for DESI ELGs is around 4 per cent, lower than previous estimations for ELGs. The mean halo mass log10(〈Mvir〉) in ${{\, h^{-1}\, \mbox{M}_\odot }}{}$ for LRGs, ELGs, and QSOs are 13.16 ± 0.01, 11.90 ± 0.06, and 12.66 ± 0.45, respectively. Our generalized SHAM algorithms facilitate the production of multitracer galaxy mocks for cosmological tests.

Tracking the evolution of satellite galaxies: mass stripping and dark-matter deficient galaxies

ABSTRACT Satellite galaxies undergo a variety of physical processes when they are accreted by groups and clusters, often resulting in the loss of baryonic and dark matter (DM) mass. In this work, we evaluate the predictions from the IllustrisTNG hydrodynamical simulation regarding the evolution of the matter content of satellites, focusing on a population that are accreted at z &amp;gt; 1 and retain their identity as satellites down to z = 0. At fixed host halo mass, the amount of DM and stellar mass stripped depends mostly on the pericentric distance, dperi, here normalized by host halo virial radius. The closest encounters result in significant loss of DM, with subhaloes retaining between 20 and a few per cent of their z = 1 mass. At fixed dperi, DM mass stripping seems more severe in lower mass haloes. Conversely, the average satellite in higher mass haloes has its stellar mass growth halted earlier, having lost a higher fraction of stellar mass by z = 0. We also show that mass stripping has a strong impact on the quenched fractions. The IllustrisTNG boxes are qualitatively consistent in these predictions, with quantitative differences mostly originating from the distinct subhalo mass ranges covered by the boxes. Finally, we have identified DM-deficient systems in all TNG boxes. These objects are preferentially found in massive clusters (Mhost ≳ 1013 M⊙), had very close encounters with their central galaxies ($d_{\rm peri}\simeq 0.05\, R_{\rm vir}$), and were accreted at high redshift (zinfall ≳ 1.4), reinforcing the notion that tidal stripping is responsible for their remarkable lack of DM.

Resolving the physics of quasar Ly α nebulae (RePhyNe): I. Constraining quasar host halo masses through circumgalactic medium kinematics

ABSTRACT Ly α nebulae ubiquitously found around z &amp;gt; 2 quasars can supply unique constraints on the properties of the circumgalactic medium, such as its density distribution, provided the quasar halo mass is known. We present a new method to constrain quasar halo masses based on the line-of-sight velocity dispersion maps of Ly α nebulae. By using MUSE-like mock observations obtained from cosmological hydrodynamic simulations under the assumption of maximal quasar fluorescence, we show that the velocity dispersion radial profiles of Ly α emitting gas are strongly determined by gravity and that they are thus self-similar with respect to halo mass when rescaled by the virial radius. Through simple analytical arguments and by exploiting the kinematics of He ii1640 Å emission for a set of observed nebulae, we show that Ly α radiative transfer effects plausibly do not change the shape of the velocity dispersion profiles but only their normalization without breaking their self-similarity. Taking advantage of these results, we define the variable $\eta ^{140-200}_{40-100}$ as the ratio of the median velocity dispersion in two specifically selected annuli and derive an analytical relation between $\eta ^{140-200}_{40-100}$ and the halo mass which can be directly applied to observations. We apply our method to 37 observed quasar Ly α nebulae at 3 &amp;lt; z &amp;lt; 4.7 and find that their associated quasars are typically hosted by ∼1012.16 ± 0.14M⊙ haloes independent of redshift within the explored range. This measurement, which is completely independent of clustering methods, is consistent with the lowest mass estimates based on quasar autocorrelation clustering at z∼3 and with quasar-galaxies cross-correlation results.