Number Density Profile Of Galaxies Research Articles

The FLAMINGO simulation view of cluster progenitors observed in the epoch of reionization with JWST

ABSTRACT Motivated by the recent JWST discovery of galaxy overdensities during the Epoch of Reionzation, we examine the physical properties of high-z protoclusters and their evolution using the Full-hydro Large-scale structure simulations with All-sky Mapping for the Interpretation of Next Generation Observations (FLAMINGO) simulation suite. We investigate the impact of the apertures used to define protoclusters, because the heterogeneous apertures used in the literature have limited our understanding of the population. Our results are insensitive to the uncertainties of the subgrid models at a given resolution, whereas further investigation into the dependence on numerical resolution is needed. When considering galaxies more massive than $M_\ast \, {\simeq }\, 10^8\, {\rm M_\odot }$, the FLAMINGO simulations predict a dominant contribution from progenitors similar to those of the Coma cluster to the cosmic star formation rate density during the reionization epoch. Our results indicate the onset of suppression of star formation in the protocluster environments as early as $z\, {\simeq }\, 5$. The galaxy number density profiles are similar to NFW (Navarro–Frenk–White profile) at $z\, {\lesssim }\, 1$ while showing a steeper slope at earlier times before the formation of the core. Different from most previous simulations, the predicted star formation history for individual protoclusters is in good agreement with observations. We demonstrate that, depending on the aperture, the integrated physical properties including the total (dark matter and baryonic) mass can be biased by a factor of 2 to 5 at $z\, {=}\, 5.5$–7, and by an order of magnitude at $z\, {\lesssim }\, 4$. This correction suffices to remove the ${\simeq }\, 3\, \sigma$ tensions with the number density of structures found in recent JWST observations.

The boundary of cosmic filaments

ABSTRACT For decades, the boundary of cosmic filaments has been a subject of debate. In this work, we determine the physically motivated radii of filaments by constructing stacked galaxy number density profiles around the filament spines. We find that the slope of the profile changes with distance to the filament spine, reaching its minimum at approximately 1 Mpc at $z=0$ in both state-of-the-art hydrodynamical simulations and observational data. This can be taken as the average value of the filament radius. Furthermore, we note that the average filament radius rapidly decreases from $z=4$ to 1, and then slightly increases. Moreover, we find that the radius of the filament depends on the length of the filament, the distance from the connected clusters, and the masses of the clusters. These results suggest a two-phase formation scenario of cosmic filaments. The filaments experienced rapid contraction before $z=1$, but their density distribution has remained roughly stable since then. The subsequent mass transport along the filaments to the connected clusters is likely to have contributed to the formation of the clusters themselves.

The Galaxy Number Density Profile of Halos

More precise measurements of galaxy clustering will be provided by the next generation of galaxy surveys, such as DESI, WALLABY, and the Square Kilometre Array. To utilize this information to improve our understanding of the Universe, we need to accurately model the distribution of galaxies in their host dark matter halos. In this work, we present a new galaxy number density profile of halos, which makes predictions for the positions of galaxies in the host halo, different to the widely adopted Navarro–Frenk–White (NFW) profile, since galaxies tend to be found more in the outskirts of halos (nearer the virial radius) than an NFW profile. The parameterized galaxy number density profile model of halos is fit and tested using the Dark Sage semi-analytic model of galaxy formation. We find that our galaxy number density profile model of halos can accurately reproduce the halo occupation distribution and galaxy two-point correlation function of the Dark Sage simulation. We also derive the analytic expressions for the circular velocity and gravitational potential energy for this profile model. We use the SDSS Data Release 10 galaxy group catalog to validate this galaxy number density profile model of halos. Compared to the NFW profile, we find that our model more accurately predicts the positions of galaxies in their host halo and the galaxy two-point correlation function.

The Uchuu-universe machine data set: galaxies in and around clusters

ABSTRACT We present the public data release of the Uchuu-UM galaxy catalogues by applying the UniverseMachine algorithm to assign galaxies to the dark matter haloes in the Uchuu N-body cosmological simulation. It includes a variety of baryonic properties for all galaxies down to ∼5 × 108 M⊙ with haloes in a mass range of 1010 < Mhalo/M⊙ < 5 × 1015 up to redshift z = 10. Uchuu-UM includes more than 104 cluster-size haloes in a volume of 8(h−1Gpc)3, reproducing observed stellar mass functions across the redshift range of z = 0−7, galaxy quenched fractions, and clustering statistics at low redshifts. Compared to the previous largest UM catalogue, the Uchuu-UM catalogue includes significantly more massive galaxies hosted by large-mass dark matter haloes. Overall, the number density profile of galaxies in dark matter haloes follows the dark matter profile, with the profile becoming steeper around the splashback radius and flattening at larger radii. The number density profile of galaxies tends to be steeper for larger stellar masses and depends on the colour of galaxies, with red galaxies having steeper slopes at all radii than blue galaxies. The quenched fraction exhibits a strong dependence on the stellar mass and increases towards the inner regions of clusters. The publicly available Uchuu-UM galaxy catalogue presented here can serve to model ongoing and upcoming large galaxy surveys.

Probing Galaxy Evolution in Massive Clusters Using ACT and DES: Splashback as a Cosmic Clock

Abstract We measure the projected number density profiles of galaxies and the splashback feature in clusters selected by the Sunyaev–Zel’dovich effect from the Advanced Atacama Cosmology Telescope (AdvACT) survey using galaxies observed by the Dark Energy Survey (DES). The splashback radius is consistent with CDM-only simulations and is located at 2.4 − 0.4 + 0.3 Mpc h − 1 . We split the galaxies on color and find significant differences in their profile shapes. Red and green-valley galaxies show a splashback-like minimum in their slope profile consistent with theory, while the bluest galaxies show a weak feature at a smaller radius. We develop a mapping of galaxies to subhalos in simulations and assign colors based on infall time onto their hosts. We find that the shift in location of the steepest slope and different profile shapes can be mapped to the average time of infall of galaxies of different colors. The steepest slope traces a discontinuity in the phase space of dark matter halos. By relating spatial profiles to infall time, we can use splashback as a clock to understand galaxy quenching. We find that red galaxies have on average been in clusters over 3.2 Gyr, green galaxies about 2.2 Gyr, while blue galaxies have been accreted most recently and have not reached apocenter. Using the full radial profiles, we fit a simple quenching model and find that the onset of galaxy quenching occurs after a delay of about a gigayear and that galaxies quench rapidly thereafter with an exponential timescale of 0.6 Gyr.

The Santiago–Harvard–Edinburgh–Durham void comparison – I. SHEDding light on chameleon gravity tests

We present a systematic comparison of several existing and new void finding algorithms, focusing on their potential power to test a particular class of modified gravity models - chameleon $f(R)$ gravity. These models deviate from standard General Relativity (GR) more strongly in low-density regions and thus voids are a promising venue to test them. We use Halo Occupation Distribution (HOD) prescriptions to populate haloes with galaxies, and tune the HOD parameters such that the galaxy two-point correlation functions are the same in both f(R) and GR models. We identify both 3D voids as well as 2D underdensities in the plane-of-the-sky to find the same void abundance and void galaxy number density profiles across all models, which suggests that they do not contain much information beyond galaxy clustering. However, the underlying void dark matter density profiles are significantly different, with f(R) voids being more underdense than GR ones, which leads to f(R) voids having a larger tangential shear signal than their GR analogues. We investigate the potential of each void finder to test f(R) models with near-future lensing surveys such as EUCLID and LSST. The 2D voids have the largest power to probe f(R) gravity, with a LSST analysis of tunnel (which is a new type of 2D underdensity introduced here) lensing distinguishing at 80 and 11$\sigma$ (statistical error) f(R) models with $|f_{R0}|=10^{-5}$ and $10^{-6}$ from GR.

THE QUENCHING TIMESCALE AND QUENCHING RATE OF GALAXIES

ABSTRACT The average star formation rate (SFR) in galaxies has been declining since the redshift of 2. A fraction of galaxies quench and become quiescent. We constrain two key properties of the quenching process: the quenching timescale and the quenching rate among galaxies. We achieve this by analyzing the galaxy number density profile in NUV−u color space and the distribution in NUV−u versus u − i color–color diagram with a simple toy-model framework. We focus on galaxies in three mass bins between 1010 and 1010.6 M ⊙. In the NUV−u versus u − i color–color diagram, the red u − i galaxies exhibit a different slope from the slope traced by the star-forming galaxies. This angled distribution and the number density profile of galaxies in NUV−u space strongly suggest that the decline of the SFR in galaxies has to accelerate before they turn quiescent. We model this color–color distribution with a two-phase exponential decline star formation history. The models with an e-folding time in the second phase (the quenching phase) of 0.5 Gyr best fit the data. We further use the NUV−u number density profile to constrain the quenching rate among star-forming galaxies as a function of mass. Adopting an e-folding time of 0.5 Gyr in the second phase (or the quenching phase), we found the quenching rate to be 19%/Gyr, 25%/Gyr and 33%/Gyr for the three mass bins. These are upper limits of the quenching rate as the transition zone could also be populated by rejuvenated red-sequence galaxies.

Observing dynamical friction in galaxy clusters

We present a novel method to detect the effects of dynamical friction in observed galaxy clusters. Following accretion into clusters, massive satellite galaxies will backsplash to systematically smaller radii than less massive satellites, an effect that may be detected by stacking the number density profiles of galaxies around clusters. We show that this effect may be understood using a simple toy model which reproduces the trends with halo properties observed in simulations. We search for this effect using SDSS redMaPPer clusters with richness 010 < λ < 2, and find that bright (Mi < −21.5) satellites have smaller splashback radii than fainter (0Mi > −2) satellites at 99% confidence.

Mass, velocity anisotropy, and pseudo phase-space density profiles of Abell 2142(Corrigendum)

The pseudo phase-space density profile of Abell 2142, defined with the galaxy number density profile ν instead of the cluster mass density profile ρ, was shown in Fig. 12 of Munari et al. (2014). This graph was erroneous (because we had incorrectly considered the projected number density profile instead of the 3D profile). The correct figure is shown below. The values in Table 4 are changed, and the correct ones are reported in the table below. While Fig. 12 of Munari et al. (2014) indicated that the PPSDs computed with the number density profile are significantly shallower than the theoretical relation of Dehnen & McLaughlin (2005), the corrected version of the figure shown here indicates that the PPSDs computed with the number density profile are now either consistent with the relation of Dehnen & McLaughlin (Q(r) for BLUE sample) or only slightly shallower, but not less consistent with that relation than found for the PPSDs computed with the mass density profile. Therefore, the statement in Munari et al. (2014) that the mass density profile represents the PPSD and β − γ relations better is no longer correct. Indeed, Figs. 11 and 13 of Munari et al. show that the β − γ relations obtained using the mass density or tracer number density are indistinguishable for the RED and ALL samples. And therefore, the entire discussion of the greater relevance of the mass density profile relative to the galaxy number density profile must be dismissed. Fig. 12. Same as Fig. 10 of Munari et al. (2014), but now using the radial profiles of galaxy number density instead of total mass density to estimate the PPSD.

Introducing the Illustris project: the evolution of galaxy populations across cosmic time

We present an overview of galaxy evolution across cosmic time in the Illustris Simulation. Illustris is an N-body/hydrodynamical simulation that evolves 2*1820^3 resolution elements in a (106.5Mpc)^3 box from cosmological initial conditions down to z=0 using the AREPO moving-mesh code. The simulation uses a state-of-the-art set of physical models for galaxy formation that was tuned to reproduce the z=0 stellar mass function and the history of the cosmic star-formation rate density. We find that Illustris successfully reproduces a plethora of observations of galaxy populations at various redshifts, for which no tuning was performed, and provide predictions for future observations. In particular, we discuss (a) the buildup of galactic mass, showing stellar mass functions and the relations between stellar mass and halo mass from z=7 to z=0, (b) galaxy number density profiles around massive central galaxies out to z=4, (c) the gas and total baryon content of both galaxies and their halos for different redshifts, and as a function of mass and radius, and (d) the evolution of galaxy specific star-formation rates up to z=8. In addition, we (i) present a qualitative analysis of galaxy morphologies from z=5 to z=0, for the stellar as well as the gaseous components, and their appearance in HST mock observations, (ii) follow galaxies selected at z=2 to their z=0 descendants, and quantify their growth and merger histories, and (iii) track massive z=0 galaxies to high redshift and study their joint evolution in star-formation activity and compactness. We conclude with a discussion of several disagreements with observations, and lay out possible directions for future research.

GALAXY ENVIRONMENTS OVER COSMIC TIME: THE NON-EVOLVING RADIAL GALAXY DISTRIBUTIONS AROUND MASSIVE GALAXIES SINCEz= 1.6

We present a statistical study of the environments of massive galaxies in four redshift bins between z=0.04 and z=1.6, using data from the Sloan Digital Sky Survey (SDSS) and the NEWFIRM Medium Band Survey (NMBS). We measure the projected radial distribution of galaxies in cylinders around a constant number density selected sample of massive galaxies and utilize a statistical subtraction of contaminating sources. Our analysis shows that massive primary galaxies typically live in group halos and are surrounded by 2 to 3 satellites with masses more than one-tenth of the primary galaxy mass. The cumulative stellar mass in these satellites roughly equals the mass of the primary galaxy itself. We further find that the radial number density profile of galaxies around massive primaries has not evolved significantly in either slope or overall normalization in the past 9.5 Gyr. A simplistic interpretation of this result can be taken as evidence for a lack of mergers in the studied groups and as support for a static evolution model of halos containing massive primaries. Alternatively, there exists a tight balance between mergers and accretion of new satellites such that the overall distribution of galaxies in and around the halo is preserved. The latter interpretation is supported by a comparison to a semi-analytic model, which shows a similar constant average satellite distribution over the same redshift range.

Effects of superstructure environment on galaxy groups

We analyse properties of galaxy groups and their dependence on the large scale environment as defined by superstructures. We find that group galaxy cross correlations depend only on group properties regardless the groups reside in superstructures. This indicates that the total galaxy density profile around groups is independent of the global environment. At a given global luminosity, a proxy to group total mass, groups have a larger stellar mass content by a factor 1.3, a relative excess independent of the group luminosity. Groups in superstructures have 40 per cent higher velocity dispersions and systematically larger minimal enclosing radii. We also find that the stellar population of galaxies in groups in superstructures is systematically older as infered from the galaxy spectra Dn4000 parameter. Although the galaxy number density profile of groups is independent of environment, the star formation rate and stellar mass profile of the groups residing in superstructures differs from groups elsewhere. For groups residing in superstructures, the combination of a larger stellar mass content and star formation rate produces a larger time scale for star formation regardless the distance to the group center. Our results provide evidence that groups in superstructures formed earlier than elsewhere, as expected in the assembly bias scenario.

The radial distribution of galaxies in groups and clusters

We present a new catalogue of 55,121 groups and clusters centred on Luminous Red Galaxies from SDSS DR7 in the redshift range 0.15<z<0.4. We provide halo mass estimates for each of these groups derived from a calibration between the optical richness of bright galaxies (M_r<-20.5) within 1 Mpc, and X-ray-derived mass for a small subset of 129 groups and clusters with X-ray measurements. We derive the mean (stacked) surface number density profiles of galaxies as a function of total halo mass in different mass bins. We find that derived profiles can be well-described by a projected NFW profile with a concentration parameter (<c>~2.6) that is approximately a factor of two lower than that of the dark matter (as predicted by N-body cosmological simulations) and nearly independent of halo mass. Interestingly, in spite of the difference in shape between the galaxy and dark matter radial distributions, both exhibit a high degree of self-similarity. A self-consistent comparison to several recent semi-analytic models of galaxy formation indicates that: (1) beyond ~0.3 r_500 current models are able to reproduce both the shape and normalisation of the satellite profiles; and (2) within ~0.3 r_500 the predicted profiles are sensitive to the details of the satellite-BCG merger timescale calculation. The former is a direct result of the models being tuned to match the global galaxy luminosity function combined with the assumption that the satellite galaxies do not suffer significant tidal stripping, even though their surrounding dark matter haloes can be removed through this process. Combining our results with measurements of the intracluster light should provide a way to inform theoretical models on the efficacy of the tidal stripping and merging processes.

THE STRUCTURE OF 2MASS GALAXY CLUSTERS

We use a sample of galaxies from the Two Micron All Sky Survey (2MASS) Extended Source Catalog to refine a matched filter method of finding galaxy clusters that takes into account each galaxy's position, magnitude, and redshift if available. The matched filter postulates a radial density profile, luminosity function, and line-of-sight velocity distribution for cluster galaxies. We use this method to search for clusters in the galaxy catalog, which is complete to an extinction-corrected K-band magnitude of 13.25 and has spectroscopic redshifts for roughly 40% of the galaxies, including nearly all brighter than K = 11.25. We then use a stacking analysis to determine the average luminosity function, radial distribution, and velocity distribution of cluster galaxies in several richness classes, and use the results to update the parameters of the matched filter before repeating the cluster search. We also investigate the correlations between a cluster's richness and its velocity dispersion and core radius, using these relations to refine priors that are applied during the cluster search process. After the second cluster search iteration, we repeat the stacking analysis. We find a cluster galaxy luminosity function that fits a Schechter form, with parameters M_K* - 5 log h = -23.64\pm0.04 and \alpha = -1.07\pm0.03. We can achieve a slightly better fit to our luminosity function by adding a Gaussian component on the bright end to represent the brightest cluster galaxy (BCG) population. The radial number density profile of galaxies closely matches a projected Navarro-Frenk-White (NFW) profile at intermediate radii, with deviations at small radii due to well-known cluster centering issues and outside the virial radius due to correlated structure. The velocity distributions are Gaussian in shape, with velocity dispersions that correlate strongly with richness.

QUANTIFYING THE COLLISIONLESS NATURE OF DARK MATTER AND GALAXIES IN A1689

We use extensive measurements of the cluster A1689 to assess the expected similarity in the dynamics of galaxies and dark matter (DM) in their motion as collisionless `particles' in the cluster gravitational potential. To do so we derive the radial profile of the specific kinetic energy of the cluster galaxies from the Jeans equation and observational data. Assuming that the specific kinetic energies of galaxies and DM are roughly equal, we obtain the mean value of the DM velocity anisotropy parameter, and the DM density profile. Since this deduced profile has a scale radius that is higher than inferred from lensing observations, we tested the validity of the assumption by repeating the analysis using results of simulations for the profile of the DM velocity anisotropy. Results of both analyses indicate a significant difference between the kinematics of galaxies and DM within $r \lesssim 0.3r_{\rm vir}$. This finding is reflected also in the shape of the galaxy number density profile, which flattens markedly with respect to the steadily rising DM profile at small radii. Thus, $r \sim 0.3r_{\rm vir}$ seems to be a transition region interior to which collisional effects significantly modify the dynamical properties of the galaxy population with respect to those of DM in A1689

DYNAMICAL STUDY OF A1689 FROM WIDE-FIELD VLT/VIMOS SPECTROSCOPY: MASS PROFILE, CONCENTRATION PARAMETER, AND VELOCITY ANISOTROPY

We examine the dynamics structure of the rich cluster A1689, combining VLT/VIMOS spectroscopy with Subaru/Suprime-Cam imaging. The radial velocity distribution of $\sim 500$ cluster members is bounded by a pair of clearly defined velocity caustics, with a maximum amplitude of $\sim|4000|$ km/s at $\simeq$ 300 h$^{-1}$ kpc, beyond which the amplitude steadily declines, approaching zero velocity at a limiting radius of $\sim$ 2 h$^{-1}$ Mpc. We derive the 3D velocity anisotropy and galaxy number density profiles using a model-independent method to solve the Jeans equation, simultaneously incorporating the observed velocity dispersion profile, the galaxy counts from deep Subaru imaging, and our previously derived cluster mass profile from a joint lensing and X-ray analysis. The velocity anisotropy is found to be predominantly radial at large radius, becoming increasingly tangential towards the center, in accord with expectations. We also analyze the galaxy data independently of our previous analysis using two different methods: The first is based on a solution of the Jeans equation assuming an NFW form for the mass distribution, whereas in the second method the caustic amplitude is used to determine the escape velocity. The cluster virial mass derived by both of these dynamical methods is in good agreement with results from our earlier lensing and X-ray analysis. We also confirm the high NFW concentration parameter, with results from both methods combined to yield $c_{\rm vir}>13$ (1$\sigma$). The inferred virial radius is consistent with the limiting radius where the caustics approach zero velocity and where the counts of cluster members drop off, suggesting that infall onto A1689 is currently not significant.

NGC 1600: Cluster or Field Elliptical?

A study of the galaxy distribution in the field of the elliptical galaxy NGC 1600 has been undertaken. Although this galaxy is often classified as a member of a loose group, all the neighboring galaxies are much fainter and could be taken as satellites of NGC 1600. The number density profile of galaxies in the field of this galaxy shows a decline with radius, with evidence of a background at approximately 1.3 Mpc. The density and number density profile are consistent with that found for other isolated early-type galaxies. NGC 1600 appears as an extended source in X-rays, and the center of the X-ray emission seems not to coincide with the center of the galaxy. The velocity distribution of neighboring galaxies has been measured from optical spectroscopic observations and shows that the mean radial velocity is approximately 85 km s?1 less than that of NGC 1600, indicating that the center of mass could lie outside the galaxy. The velocity dispersion of the group is estimated at -->429 ? 57 km s?1. The inferred mass of the system is therefore of the order of -->1014 -->M?, a value that corresponds to a large group. NGC 1600 therefore shares some similarities, but is not identical to, the fossil clusters detected in X-ray surveys. Implications of this result for studies of isolated early-type galaxies are briefly discussed.

Statistics of voids in the two-degree Field Galaxy Redshift Survey

We present a statistical analysis of voids in the two-degree Field Galaxy Redshift Survey (2dFGRS). In order to detect the voids, we have developed two robust algorithms. We define voids as non-overlapping maximal spheres empty of haloes or galaxies with mass or luminosity above a given value. We search for voids in cosmological N-Body simulations to test the performance of our void finders. We obtain and analyse the void statistics for several volume-limited samples for the North Galactic Pole (NGP) and the South Galactic Pole (SGP) constructed from the 2dFGRS full data release. We find that the results obtained from the NGP and the SGP are statistically compatible. From the results of several statistical tests we conclude that voids are essentially uncorrelated, with at most a mild anticorrelation and that at the 99.5 per cent confidence level there is a dependence of the void number density on redshift. We develop a technique to correct the distortion caused by the fact that we use the redshift as the radial coordinate. We calibrate this technique with mock catalogues and find that the correction might be of some relevance to carry out accurate inferences from void statistics. We study the statistics of the galaxies inside nine nearby voids. We find that galaxies in voids are not randomly distributed: they form structures like filaments. We also obtain the galaxy number density profile in voids. This profile follow a similar but steeper trend to that followed by haloes in voids.

Evolution of CDM halos and the Milky Way satellites

A cold dark matter halo big enough to host the Milky Way contains hundreds of subhalos massive enough to host dwarf galaxies. The difference to the much smaller observed number of satellite galaxies seems to be a problem for CDM. The galaxy number density profile and disk like configuration are also different form the total subhalo populations in CDM simulations. A number of different models of dwarf galaxy formation which are able to reproduce the right number of luminous subhalos have been proposed. Some of them also give the right radial distributions and make disk like configurations more probable. Additional information about the typical formation times and sites of dwarf galaxies can be found in the stellar halo of the Milky Way, i.e. from the stellar debris of tidally disrupted dwarfs. A stellar halo with a realistic concentration is obtained when most dwarfs form early (before redshift 10) in small halos (virial mass above which is needed for efficient atomic cooling.

Central matter distributions in rich clusters of galaxies from ${z\sim{0}}$ to ${z\sim {0.5}}$

We have analyzed the galaxy number density and luminosity density profiles of rich clusters of galaxies from redshifts to . We show that the luminosity profile computed with bright galaxies is significantly cusped in the center of the clusters, whatever the redshift. This is in agreement with the dark matter profiles predicted by numerical simulations. The galaxy number density profile for the bright galaxies is fitted equally well with a core model or a cusped model. In contrast, the luminosity and the galaxy number density profiles of the fainter galaxies are significantly better fitted by a core model. We did not detect any statistically significant different fits when applied to data in the range from to . The difference in profile between faint and bright galaxies may be due to the rapid (relative to the age of the universe at versus ) destruction of the faint galaxies by tidal forces and merging events in the denser central regions of the clusters. This process could erase the cusp by turning faint galaxies into diffuse light. In this case, the galaxies (with a cusp visible in the bright galaxy number density and mainly in luminosity profiles) would trace the total mass distribution.