Spatial Distribution Of Galaxies Research Articles

The Faber–Jackson relation and Fundamental Plane from halo abundance matching

The Fundamental Plane (FP) describes the relation between the stellar mass, size, and velocity dispersion of elliptical galaxies; the Faber-Jackson relation (FJR) is its projection onto {mass, velocity} space. In this work we redeploy and expand the framework of Desmond & Wechsler (2015) to ask whether abundance matching-based LCDM models that have shown success in matching the spatial distribution of galaxies are also capable of explaining key properties of the FJR and FP, including their scatter. Within our framework, agreement with the normalisation of the FJR requires haloes to expand in response to disc formation. We find that the tilt of the FP may be explained by a combination of the observed non-homology in galaxy structure and the variation in mass-to-light ratio produced by abundance matching with a universal initial mass function (IMF), provided that the anisotropy of stellar motions is taken into account. However, the predicted scatter around the FP is considerably increased by situating galaxies in cosmologically-motivated haloes due to variations in halo properties at fixed stellar mass, and appears to exceed that of the data. This implies that additional correlations between galaxy and halo variables may be required to fully reconcile these models with elliptical galaxy scaling relations.

A unified model for the spatial and mass distribution of subhaloes

N-body simulations suggest that the substructures that survive inside dark matter haloes follow universal distributions in mass and radial number density. We demonstrate that a simple analytical model can explain these subhalo distributions as resulting from tidal stripping which increasingly reduces the mass of subhaloes with decreasing halo-centric distance. As a starting point, the spatial distribution of subhaloes of any given infall mass is shown to be largely indistinguishable from the overall mass distribution of the host halo. Using a physically motivated statistical description of the amount of mass stripped from individual subhaloes, the model fully describes the joint distribution of subhaloes in final mass, infall mass and radius. As a result, it can be used to predict several derived distributions involving combinations of these quantities including, but not limited to, the universal subhalo mass function, the subhalo spatial distribution, the gravitational lensing profile, the dark matter annihilation radiation profile and boost factor. This model clarifies a common confusion when comparing the spatial distributions of galaxies and subhaloes, the so called "anti-bias", as a simple selection effect. We provide a Python code SubGen for populating haloes with subhaloes at http://icc.dur.ac.uk/data/

Optical-near-IR analysis of globular clusters in the IKN dwarf spheroidal: a complex star formation history

Age, metallicity and spatial distribution of globular clusters (GCs) provide a powerful tool to reconstruct major star-formation episodes in galaxies. IKN is a faint dwarf spheroidal (dSph) in the M81 group of galaxies. It contains five old GCs, which makes it the galaxy with the highest known specific frequency (SN=126). We estimate the photometric age, metallicity and spatial distribution of the poorly studied IKN GCs. We search SDSS for GC candidates beyond the HST field of view, which covers half of IKN. To break the age-metallicity degeneracy in the V-I colour we use WHT/LIRIS Ks-band photometry and derive photometric ages and metallicities by comparison with SSP models in the V,I,Ks colour space. IKN GCs' VIKs colours are consistent with old ages ($\geq\!8$ Gyr) and a metallicity distribution with a higher mean than typical for such a dSph ([Fe/H$]\!\simeq\!-1.4_{-0.2}^{+0.6}$ dex). Their photometric masses range ($0.5 <{\cal M_{\rm GC}}<4\times10^5M_\odot$) implies a high mass ratio between GCs and field stars, of $10.6\%$. Mixture model analysis of the RGB field stars' metallicity suggests that 72\% of the stars may have formed together with the GCs. Using the most massive GC-SFR relation we calculate a SFR of $\sim\!10M_\odot/$yr during its formation epoch. We note that the more massive GCs are closer to the galaxy photometric centre. IKN GCs also appear spatially aligned along a line close to the IKN major-axis and nearly orthogonal to the plane of spatial distribution of galaxies in the M81 group. We identify one new IKN GC candidate based on colour and PSF analysis of the SDSS data. The evidence towards i) broad and high metallicity distribution of the field IKN RGB stars and its GCs, ii) high fraction and iii), spatial alignment of IKN GCs, supports a scenario for tidally triggered complex IKN's SFH in the context of interactions with galaxies in the M81 group.

THE SPATIAL DISTRIBUTION OF SATELLITE GALAXIES WITHIN HALOS: MEASURING THE VERY SMALL SCALE ANGULAR CLUSTERING OF SDSS GALAXIES

We measure the angular clustering of galaxies from the Sloan Digital Sky Survey (SDSS) Data Release 7 in order to probe the spatial distribution of satellite galaxies within their dark matter halos. Specifically, we measure the angular correlation function on very small scales (7??320?) in a range of luminosity threshold samples (absolute r-band magnitudes from ?18 up to ?21) that are constructed from the subset of SDSS that has been spectroscopically observed more than once (the so-called plate overlap region). We choose to measure angular clustering in this reduced survey footprint in order to minimize the effects of fiber collision incompleteness, which are otherwise substantial on these small scales, and we discuss the possible impact that fiber collisions have on our measurements. We model our clustering measurements using a fully numerical halo model that populates dark matter halos in N-body simulations to create realistic mock galaxy catalogs. The model has free parameters that specify both the number and spatial distribution of galaxies within their host halos. We adopt a flexible density profile for the spatial distribution of satellite galaxies that is similar to the dark matter Navarro?Frenk?White (NFW) profile, except that the inner slope is allowed to vary. We find that the angular clustering of our most luminous samples ( < ?20 and ?21) suggests that luminous satellite galaxies have substantially steeper inner density profiles than NFW. Lower-luminosity samples are less constraining, however, and are consistent with satellite galaxies having shallow density profiles. Our results confirm the findings of Watson et al. while using different clustering measurements and modeling methodology.

Finding and characterising WHIM structures using the luminosity density method

AbstractWe have developed a new method to approach the missing baryons problem. We assume that the missing baryons reside in a form of Warm Hot Intergalactic Medium, i.e. the WHIM. Our method consists of (a) detecting the coherent large scale structure in the spatial distribution of galaxies that traces the Cosmic Web and that in hydrodynamical simulations is associated to the WHIM, (b) mapping its luminosity into a galaxy luminosity density field, (c) using numerical simulations to relate the luminosity density to the density of the WHIM, (d) applying this relation to real data to trace the WHIM using the observed galaxy luminosities in the Sloan Digital Sky Survey and 2dF redshift surveys. In our application we find evidence for the WHIM along the line of sight to the Sculptor Wall, at redshifts consistent with the recently reported X-ray absorption line detections. Our indirect WHIM detection technique complements the standard method based on the detection of characteristic X-ray absorption lines, showing that the galaxy luminosity density is a reliable signpost for the WHIM. For this reason, our method could be applied to current galaxy surveys to optimise the observational strategies for detecting and studying the WHIM and its properties. Our estimates of the WHIM hydrogen column density NH in Sculptor agree with those obtained via the X-ray analysis. Due to the additional NH estimate, our method has potential for improving the constrains of the physical parameters of the WHIM as derived with X-ray absorption, and thus for improving the understanding of the missing baryons problem.

Determination of the abundance of cosmic matter via the cell count moments of the galaxy distribution

We demonstrate that accurate and precise information about the matter content of the universe can be retrieved via a simple cell count analysis of the 3D spatial distribution of galaxies. A new clustering statistic, the {\it galaxy clustering ratio} $\eta$, is the key in this process. This is defined as the ratio between one- and two-point second-order moments of the smoothed galaxy density distribution. The distinguishing feature of this statistic is its universality: on large cosmic scales both galaxies (in redshift space) and mass (in real space) display the same $\eta$ amplitude. This quantity, in addition, does not evolve as a function of redshift. As a consequence, the $\eta$ statistic provides insight into characteristic parameters of the real-space power spectrum of mass density fluctuations without the need to specify the galaxy biasing function, a model for galaxy redshift distortions nor the growing mode of density ripples. We demonstrate the method with the luminous red galaxy (LRG) sample extracted from the spectroscopic Sloan Digital Sky Survey (SDSS) data release 7 (DR7) catalogue. Taking weak (flat) priors of the curvature of the universe ($\Omega_k$) and of the constant value of the dark energy equation of state ($w$), and strong (gaussian) priors of the physical baryon density $\Omega_bh^2$, of the Hubble constant $H_0$ and of the spectral index of primordial density perturbations $n_s$, we estimate the abundance of matter with a relative error of 8% ($\Omega_m= 0.283\pm0.023$). We expect that this approach will be instrumental in searching for evidence of new physics beyond the standard model of cosmology and in planning future redshift surveys such as BigBOSS or EUCLID.

Anisotropy in the Hubble constant as modeled by density gradients

The all-sky maps of the observed variation of the Hubble constant can be reproduced from a theoretical point of view by introducing an intergalactic plasma with a variable number density of electrons. The observed averaged value and variance of the Hubble constant are reproduced by adopting a rim model, an auto-gravitating model, and a Voronoi diagrams model as the backbone for an auto-gravitating medium. We also analyze an astronomer's model based on the 3D spatial distribution of galaxies as given by the 2MASS Redshift Survey and an auto-gravitating Lane--Emden ($n=5$) profile of the electrons. The simulation which involves the Voronoi diagrams is done in a cubic box with sides of 100 Mpc. The simulation which involves the 2MASS covers the range of redshift smaller than 0.05.

Astronomical Signatures of Dark Matter

Several independent astronomical observations in different wavelength bands reveal the existence of much larger quantities of matter than what we would deduce from assuming a solar mass to light ratio. They are very high velocities of individual galaxies within clusters of galaxies, higher than expected rotation rates of stars in the outer regions of galaxies, 21 cm line studies indicative of increasing mass to light ratios with radius in the halos of spiral galaxies, hot gaseous X-ray emitting halos around many elliptical galaxies, and clusters of galaxies requiring a much larger component of unseen mass for the hot gas to be bound. The level of gravitational attraction needed for the spatial distribution of galaxies to evolve from the small perturbations implied by the very slightly anisotropic cosmic microwave background radiation to its current web-like configuration requires much more mass than is observed across the entire electromagnetic spectrum. Distorted shapes of galaxies and other features created by gravitational lensing in the images of many astronomical objects require an amount of dark matter consistent with other estimates. The unambiguous detection of dark matter and more recently evidence for dark energy has positioned astronomy at the frontier of fundamental physics as it was in the 17th century.

Groups in the Millennium Simulation and in SDSS DR7

The Millennium N-body simulation and the Sloan Digital Sky Survey seventh data release (SDSS DR7) galaxy and galaxy group catalogues are compared to study the properties of galaxy groups and the distribution of galaxies in groups. We construct mock galaxy group catalogues for a Millennium semi-analytical galaxy catalogue by using the same friends-of-friends method, which was used by Tago et al to analyse the SDSS data. We analyse in detail the group luminosities, group richnesses, virial radii, sizes of groups and their rms velocities for four volume-limited samples from observations and simulations. Our results show that the spatial densities of groups agree within one order of magnitude in all samples with a rather good agreement between the mock catalogues and observations. All group property distributions have similar shapes and amplitudes for richer groups. For galaxy pairs and small groups, the group properties for observations and simulations are clearly different. In addition, the spatial distribution of galaxies in small groups is different: at the outskirts of the groups the galaxy number distributions do not agree, although the agreement is relatively good in the inner regions. Differences in the distributions are mainly due to the observational limitations in the SDSS sample and to the problems in the semi-analytical methods that produce too compact and luminous groups.

Photometric Effects and Voronoi-Diagrams as a Mixed Model for the Spatial Distribution of Galaxies

We review the model of the Voronoi Diagrams which allows to reproduce the large-scale structures of our universe as given by the astronomical catalogs. The observed number of galaxies in a given solid angle with a chosen flux/magnitude versus the redshift presents a maximum that is a function of the flux/magnitude ; it can be explained by a detailed analysis of the standard luminosity function for galaxies as well by two new luminosity function for galaxies. The current status of the research on the statistics of the Voronoi Diagrams is reviewed.

The scale of cosmic isotropy

The most fundamental premise to the standard model of the universe states that the large-scale properties of the universe are the same in all directions and at all comoving positions. Demonstrating this hypothesis has proven to be a formidable challenge.The cross-over scale Riso above which the galaxy distribution becomes statistically isotropic is vaguely defined and poorly (if not at all) quantified. Here we report on a formalism that allows us to provide an unambiguous operational definition and an estimate of Riso. We apply the method to galaxies in the Sloan Digital Sky Survey (SDSS) Data Release 7, finding that Riso ∼ 150h−1Mpc. Besides providing a consistency test of the Copernican principle, this result is in agreement with predictions based on numerical simulations of the spatial distribution of galaxies in cold dark matter dominated cosmological models.

A model of the anisotropic correlation function ξ (rp, π) in redshift space including redshift errors

With the advent of very large volume, wide-angle photometric redshift surveys like e.g. Pan-STARRS, DES, or PAU, which aim at using the spatial distribution of galaxies as a means to constrain the equation of state parameter of dark energy, w_DE, it has become extremely important to understand the influence of redshift inaccuracies on the measurement. We have developed a new model for the anisotropic two point large-scale (r > 64 h^-1 Mpc) correlation function xi(rp,pi), in which nonlinear structure growth and nonlinear coherent infall velocities are taken into account, and photometric redshift errors can easily be incorporated. In order to test its validity and investigate the effects of photometric redshifts, we compare our model with the correlation function computed from a suite of 50 large-volume, moderate-resolution numerical N-body simulation boxes, where we can perform the analysis not only in real- and redshift space, but also simulate the influence of a gaussian redshift error distribution with an absolute rms of sigma_z= 0.015, 0.03, 0.06, and 0.12, respectively. We conclude that for the given volume (V_box =2.4 h^-3 Gpc^3) and number density (n ~ 1.25 10^-4) of objects the full shape of xi(rp,pi) is modeled accurately enough to use it to derive unbiased constraints on the equation of state parameter of dark energy w_DE and the linear bias b, even in the presence of redshift errors of the order of sigma_z = 0.06.

THE EXTREME SMALL SCALES: DO SATELLITE GALAXIES TRACE DARK MATTER?

We investigate the radial distribution of galaxies within their host dark matter halos by modeling their small-scale clustering, as measured in the Sloan Digital Sky Survey. Specifically, we model the Jiang et al. (2011) measurements of the galaxy two-point correlation function down to very small projected separations (10 < r < 400 kpc/h), in a wide range of luminosity threshold samples (absolute r-band magnitudes of -18 up to -23). We use a halo occupation distribution (HOD) framework with free parameters that specify both the number and spatial distribution of galaxies within their host dark matter halos. We assume that the first galaxy in each halo lives at the halo center and that additional satellite galaxies follow a radial density profile similar to the dark matter Navarro-Frenk-White (NFW) profile, except that the concentration and inner slope are allowed to vary. We find that in low luminosity samples, satellite galaxies have radial profiles that are consistent with NFW. M_r < -20 and brighter satellite galaxies have radial profiles with significantly steeper inner slopes than NFW (we find inner logarithmic slopes ranging from -1.6 to -2.1, as opposed to -1 for NFW). We define a useful metric of concentration, M_(1/10), which is the fraction of satellite galaxies (or mass) that are enclosed within one tenth of the virial radius of a halo. We find that M_(1/10) for low luminosity satellite galaxies agrees with NFW, whereas for luminous galaxies it is 2.5-4 times higher, demonstrating that these galaxies are substantially more centrally concentrated within their dark matter halos than the dark matter itself. Our results therefore suggest that the processes that govern the spatial distribution of galaxies, once they have merged into larger halos, must be luminosity dependent, such that luminous galaxies become poor tracers of the underlying dark matter.

The evolution of K* and the halo occupation distribution since z= 1.5: observations versus simulations

We study the evolution of the K-band luminosity function (LF) and the Halo Occupation Distribution (HOD) using Subaru observations of 15 X-ray clusters at z=0.8-1.5 and compare the results with mock clusters (0<z<1.3) extracted from the Millennium Simulation and populated with galaxies using the semi-analytic model (SAM) of Bower et al., matched in mass to our observed sample. We find that the characteristic luminosity K* defined by a Shechter LF is consistent with SAM predictions, which mimic well the evolution of K* in z>1 rich clusters. However, we cannot distinguish between this model and a simple stellar population synthesis model invoking passive evolution with a formation redshift z~5 - consistent with the presence of an old red galaxy population ubiquitous in rich clusters at z=1.5. We also see a small difference (\Delta K*~0.5) between our clusters and studies of the field population at similar redshifts, suggesting only a weak dependence of the luminous (L>L*) part of the LF on environment. Turning to our HOD study, we find that within R_{500}, high-z clusters tend to host smaller numbers of galaxies to a magnitude K*+2 compared to their low-z counterparts. This behavior is also seen in the mock samples and is relatively insensitive to the average mass of the cluster haloes. In particular, we find significant correlations of the observed number of member cluster galaxies (N) with both z and cluster mass: $N(M,z)=(53\pm1)(1+z)^{-0.61^{+0.18}_{-0.20}}(M/10^{14.3})^{0.86\pm0.05}$. Finally, we examine the spatial distribution of galaxies and provide a new estimate of the concentration parameter for clusters at high z ($c_{g}=2.8^{+1.0}_{-0.8}$). Our result is consistent with predictions from both our SAM mock clusters and literature's predictions for dark matter haloes. The mock sample predictions rise slowly with decreasing redshift reaching $c_{g}=6.3^{+0.39}_{-0.36}$ at z=0.

Spatial galaxy distribution function for a three-component system

We derive the distribution function and the allied thermodynamic quantities for a system of galaxies with three mass species. A new clustering parameter b3 that inherently takes into account the masses and the number of galaxies of each kind, emerges directly from the calculations. Our general conclusion is that the inclusion of the third component does not significantly effect the overall features of the distribution function.

The non-uniform distribution of galaxies from data of the SDSS DR7 survey

We have analyzed the spatial distribution of galaxies from the release of the Sloan Digital Sky Survey of galactic redshifts (SDSS DR7), applying the complete correlation function (conditional density), two-point conditional density (cylinder), and radial density methods. Our analysis demonstrates that the conditional density has a power-law form for scales lengths 0.5-30 Mpc/h, with the power-law corresponding to the fractal dimension D = 2.2+-0.2; for scale lengths in excess of 30 Mpc/h, it enters an essentially flat regime, as is expected for a uniform distribution of galaxies. However, in the analysis applying the cylinder method, the power-law character with D = 2.0+-0.3 persists to scale lengths of 70 Mpc/h. The radial density method reveals inhomogeneities in the spatial distribution of galaxies on scales of 200 Mpc/h with a density contrast of two, confirming that translation invariance is violated in the distribution of galaxies to 300 Mpc/h, with the sampling depth of the SDSS galaxies being 600 Mpc/h.

From Dwarf Spheroidals to cDs: Simulating the Full Galaxy Population in a LCDM Cosmology

We use the combination of the very large Millennium and Millennium-II simulations to follow galaxy formation over a stellar mass range of more than five orders of magnitude throughout representative cosmological volumes. These are the first simulations to address the abundance and spatial distribution of galaxies ranging from cD’s down to dSph’s. We demonstrate that stronger feedback is needed than in previous models, epically at low mass, if the SDSS stellar mass function is to be matched in a LCDM universe. Models which satisfy this constraint also reproduce the observed abundance of faint dwarfs around the Milky Way and M 31. Only the very faintest objects are significantly affected by the treatment of reionisation. Our model can reproduce most galaxy properties in the local Universe. Problems remain with respect to the stellar age distributions and kinematics of dwarfs.

COUNTS-IN-CYLINDERS IN THE SLOAN DIGITAL SKY SURVEY WITH COMPARISONS TON-BODY SIMULATIONS

Environmental statistics provide a necessary means of comparing the properties of galaxies in different environments and a vital test of models of galaxy formation within the prevailing, hierarchical cosmological model. We explore counts-in-cylinders, a common statistic defined as the number of companions of a particular galaxy found within a given projected radius and redshift interval. Galaxy distributions with the same two-point correlation functions do not necessarily have the same companion count distributions. We use this statistic to examine the environments of galaxies in the Sloan Digital Sky Survey, Data Release 4. We also make preliminary comparisons to four models for the spatial distributions of galaxies, based on N-body simulations, and data from SDSS DR4 to study the utility of the counts-in-cylinders statistic. There is a very large scatter between the number of companions a galaxy has and the mass of its parent dark matter halo and the halo occupation, limiting the utility of this statistic for certain kinds of environmental studies. We also show that prevalent, empirical models of galaxy clustering that match observed two- and three-point clustering statistics well fail to reproduce some aspects of the observed distribution of counts-in-cylinders on 1, 3 and 6-Mpc/h scales. All models that we explore underpredict the fraction of galaxies with few or no companions in 3 and 6-Mpc/h cylinders. Roughly 7% of galaxies in the real universe are significantly more isolated within a 6 Mpc/h cylinder than the galaxies in any of the models we use. Simple, phenomenological models that map galaxies to dark matter halos fail to reproduce high-order clustering statistics in low-density environments.

The ergodicity bias in the observed galaxy distribution

The spatial distribution of galaxies we observed is subject to the given condition thatwe, human beings are sitting right in a galaxy — the Milky Way. Thus theergodicity assumption is questionable in interpretation of the observedgalaxy distribution. The resultant difference between observed statistics (volume average) and the true cosmic value (ensemble average) istermed as the ergodicity bias. We perform explicit numerical investigation ofthe effect for a set of galaxy survey depths and near-end distance cuts. It is found that the ergodicity bias in observed two- and three-point correlationfunctions in most cases is insignificant for modern analysis of samples fromgalaxy surveys and thus close a loophole in precision cosmology. However, itmay become non-negligible in certain circumstances, such as thoseapplications involving three-point correlation function at large scales oflocal galaxy samples. Thus one is reminded to take extra care in galaxysample construction and interpretation of the statistics of the sample, especially when the characteristic redshift is low.

Method for analyzing the spatial distribution of galaxies on gigaparsec scales. II. Application to a grid of the HUDF-FDF-COSMOS-HDF surveys

Using the deep fields of COSMOS, FDF, HUDF, and HDF-N as an example, we discuss the prospects for and limitations on the method for searching for superlarge structures in the spatial distribution of galaxies proposed in the preceding article of this series. An analysis of the distribution N(z) of photometric redshifts in a grid of the deep fields of HUDF-FDF-COSMOS-HDFN reveals the possible existence of superlarge structures with a contrast dN/N∼50% and tangential and radial dimensions of about 1000 Mpc. The reality of the detected candidate superlarge structures in the universe can be verified by further observations with a finer grid of deep fields. The influence of systematic errors can be reduced by observing the same deep fields with several 3-10 meter telescopes and utilizing different methods for determining the photometric redshifts.