We examine the counts-in-cells probability distribution functions that describe dark matter halos in the Dark Energy Universe Simulations (DEUS). We describe the measurements between redshifts $z=0$ to $z=4$ on both linear and non-linear scales. The best-fits of the gravitational quasi-equilibrium distribution (GQED), the negative binomial distribution (NBD), the Poisson-Lognormal distribution (PLN), and the Poisson-Lognormal distribution with a bias parameter (PLNB) are compared to simulations. The fits agree reasonably consistently over a range of redshifts and scales. To distinguish quintessence (RPCDM) and phantom ($w$CDM) dark energy from $\Lambda$ dark energy, we present a new method that compares the model parameters of the counts-in-cells probability distribution functions. We find that the mean and variance of the halo counts-in-cells on $2-25h^{-1}$Mpc scales between redshifts $0.65<z<4$ show significant percentage differences for different dark energy cosmologies. On $15-25h^{-1}$Mpc scales, the $g$ parameter in NBD, $\omega$ parameter in PLN, $b$ and $C_b$ parameters in PLNB show larger percentage differences for different dark energy cosmologies than on smaller scales. On $2-6h^{-1}$Mpc scales, kurtosis and the $b$ parameter in the GQED show larger percentage differences for different dark energy cosmologies than on larger scales. For cosmologies explored in the DEUS simulations, the percentage differences between these statistics for the RPCDM and $w$CDM dark energy cosmologies relative to $\Lambda$CDM generally increases with redshift from a few percent to significantly larger percentages at $z=4$. Applying our method to simulations and galaxy surveys can provide a useful way to distinguish among dark energy models and cosmologies in general.

We identify low redshift clusters and groups in the Sloan Digital Sky Survey (SDSS) and estimate their kinetic and correlation potential energies. We compare the distribution of these energies to the predictions by Yang and Saslaw (2012) and in the process estimate a measure of an average 3-dimensional velocity and spatial anisotropy of a sample of clusters. We find that the inferred velocity anisotropy is correlated with the inferred spatial anisotropy. We also find that the general shape of the energy distribution agrees with theory over a wide range of scales from small groups to superclusters once the uncertainties and fluctuations in the estimated energies are included.

We develop a general theory for estimating the probability that a galaxy cluster of a given shape exists. The theory is based on the observed result that the distribution of galaxies is very close to quasi-equilibrium, in both its linear and nonlinear regimes. This places constraints on the spatial configuration of a cluster of galaxies in quasi-equilibrium. In particular, we show that that a cluster of galaxies may be described as a collection of nearly virialized subclusters of approximately the same mass. Clusters that contain more than 10 subclusters are very likely to be completely virialized. Using our theory, we develop a method for comparing probabilities of different spatial configurations of subclusters. As an illustrative example, we show that a cluster of galaxies arranged in a line is more likely to occur than a cluster of galaxies arranged in a ring.

We determine the galaxy counts-in-cells distribution from the Sloan Digital Sky Survey (SDSS) for 3D spherical cells in redshift space as well as for 2D projected cells. We find that cosmic variance in the SDSS causes the counts-in-cells distributions in different quadrants to differ from each other by up to 20%. We also find that within this cosmic variance, the overall galaxy counts-in-cells distribution agrees with both the gravitational quasi-equilibrium distribution and the negative binomial distribution. We also find that brighter galaxies are more strongly clustered than if they were randomly selected from a larger complete sample that includes galaxies of all luminosities. The results suggest that bright galaxies could be in dark matter haloes separated by less than ~10 Mpc/h.

We examine the evolution of the spatial counts-in-cells distribution of galaxies and show that the form of the galaxy distribution function does not change significantly as galaxies merge and evolve. In particular, bound merging pairs follow a similar distribution to that of individual galaxies. From the adiabatic expansion of the universe we show how clustering, expansion and galaxy mergers affect the clustering parameter b. We also predict the evolution of b with respect to redshift.

This is the first exploration of the galaxy distribution function at redshifts greater than about 0.1. Redshifts are based on the North and South GOODS Catalogs. In each catalog we examine clustering in the two redshift bands 0.47 < z < 0.8 and 0.9 < z < 1.5. The mean redshifts of the samples in these bands are about 0.6 and 1.1. Our main result is that at these redshifts the galaxy spatial distribution function f_V(N) has the form predicted by gravitational quasi-equilibrium dynamics for cosmological many-body systems. This constrains related processes such as galaxy merging and the role of dark matter in the range of these redshifts.

The cosmological many-body problem is effectively an infinite system of gravitationally interacting masses in an expanding universe. Despite the interactions' long-range nature, an analytical theory of statistical mechanics describes the spatial and velocity distribution functions which arise in the quasi-equilibrium conditions that apply to many cosmologies. Consequences of this theory agree well with the observed distribution of galaxies. Further consequences such as thermodynamics provide insights into the physical properties of this system, including its robustness to mergers, and its transition from a grand canonical ensemble to a collection of microcanonical ensembles with negative specific heat.

We determine the spatial distribution function of galaxies from a wide range of samples in the 2MASS survey. The results agree very well with the form of the distribution predicted by the theory of cosmological gravitational many-body galaxy clustering. On large scales we find a value of the clustering parameter b = 0.867 +/- 0.026, in agreement with b = 0.83 +/- 0.05 found previously for the Pisces-Perseus supercluster. We measure b(theta) as a function of scale, since this is a powerful test of the applicability of computer simulations. The results suggest that when galaxies clustered they were usually surrounded by individual, rather than by communal haloes.

We examine the peculiar velocity distribution function of galaxies in cosmological many-body gravitational clustering. Our statistical mechanical approach derives a previous basic assumption and generalizes earlier results to galaxies with haloes. Comparison with the observed peculiar velocity distributions indicates that individual massive galaxies are usually surrounded by their own haloes, rather than being embedded in common haloes. We then derive the density of energy states, giving the probability that a randomly chosen configuration of N galaxies in space is bound and virialized. Gravitational clustering is very efficient. The results agree well with the observed probabilities for finding nearby groups containing N galaxies. A consequence is that our local relatively low mass group is quite typical, and the observed small departures from the local Hubble flow beyond our group are highly probable.

We give the first determination of the observed peculiar velocity distribution function for a representative sample of galaxies (within 50 Mpc (H=100) of the Local Group), which includes a wide range of clustering properties. We explore in detail the effects of uncertainties in sampling and in distance measures on the estimated distribution function. The observed distribution function is consistent with an earlier prediction of gravitational clustering, over the entire range of peculiar velocities, from field galaxies to rich clusters, on scales up to 50 Mpc (H=100). In the simplest consistent model, most of the inhomogeneous mass of the Universe is in galaxies or their halos. We estimate the ``Mach Number'' for the bulk flow within 50 Mpc (H=100) from us to be approximately 0.8, which includes the effect of high-dispersion galaxies in clusters.