We revisit the longstanding question of whether first brightest cluster galaxies are statistically drawn from the same distribution as other cluster galaxies or are "special", using the new non-parametric, empirically based model presented in Vale&Ostriker (2006) for associating galaxy luminosity with halo/subhalo masses. We introduce scatter in galaxy luminosity at fixed halo mass into this model, building a conditional luminosity function (CLF) by considering two possible models: a simple lognormal and a model based on the distribution of concentration in haloes of a given mass. We show that this model naturally allows an identification of halo/subhalo systems with groups and clusters of galaxies, giving rise to a clear central/satellite galaxy distinction. We then use these results to build up the dependence of brightest cluster galaxy (BCG) magnitudes on cluster luminosity, focusing on two statistical indicators, the dispersion in BCG magnitude and the magnitude difference between first and second brightest galaxies. We compare our results with two simple models for BCGs: a statistical hypothesis that the BCGs are drawn from a universal distribution, and a cannibalism scenario merging two galaxies from this distribution. The statistical model is known to fail from work as far back as Tremaine&Richstone (1977). We show that neither the statistical model nor the simplest possibility of cannibalism provide a good match for observations. Our CLF models both give similar results, in good agreement with observations. Specifically, we find <m_1> between -25 and -25.5 in the K-band, sigma(m_1)~0.25 and <Delta_12> between 0.6 and 0.8, for cluster luminosities in the range of 10^12 to 10^13 h^-2 L_sun.

We present a non-parametric, empirically based, model for associating galaxy luminosities with halo/subhalo masses, based on a self-consistent treatment of subhalo mass loss and the subhalo mass function. We find that, at high mass, the mass-luminosity relation is almost independent of the actual luminosity function considered, when luminosity is scaled by the characteristic luminosity L*. Additionally, the shape of the total halo luminosity depends on the slope of the subhalo mass function. For these high mass, cluster sized haloes, we find that total luminosity scales as L_tot ~ M^0.88, while the luminosity of the first brightest galaxy has a much weaker dependence on halo mass, L_1 ~ M^0.28, in good agreement with observations and previous results. At low mass, the resulting slope of the mass-luminosity relation depends strongly of the faint end slope of the luminosity function, and we obtain a steep relation, with approximately L ~ M^4.5 in the K-band. The average number of galaxies per halo/cluster is also in very good agreement with observations, scaling as M^0.9. In general, we obtain a good agreement with several independent sets of observational data. We find that, when comparing with observations and for a flat cosmology, the model tends to prefer lower values for Omega_m and sigma_8. Within the WMAP+SDSS concordance plane of Tegmark et al. (2004), we find best agreement around Omega_m=0.25 and sigma_8=0.8, also in very good agreement with the results of the CMB+2dF study of Sanchez et al. (2005). We also check on possible corrections for observed mass based on a comparison of the equivalent number of haloes/clusters. Additionally, we include further checks on the model results based on the mass to light ratio, the occupation number, the group luminosity function and the multiplicity function. (abridged)

In this paper we present a new, essentially empirical, model for the relation between the mass of a dark matter halo/subhalo and the luminosity of a galaxy hosted in it. To estimate this, we replace the assumption of linearity between light and mass fluctuations with the assumption of monotonicity between galaxy light and halo or subhalo mass. We are enabled to proceed with this less restrictive ansatz by the availability of new, very high resolution dark matter simulations and more detailed and comprehensive global galactic luminosity functions. We find that the relation between halo/subhalo mass and hosted galaxy luminosity is fairly well fit by a double power law, while that between halo mass and group luminosity is fairly well fit by a two branch function, with both branches double power laws. Both relations asymptote to $L\propto M^4$ at low $M$, while at high mass the former follows $L\propto M^{0.28}$ and the latter $L\propto M^{0.9}$. We also derive results for the occupation number, luminosity function of cluster galaxies, group luminosity function and multiplicity function. Then, using some additional assumptions, we further derive results for biasing between mass and light and mass and galaxy number, light distribution function and the void probability distribution. Our results for the most part seem to match well with observations and previous expectations. We feel this is a potentially powerful way of modelling this relation, since the main inputs are readily testable against dark matter simulation results and galaxy surveys. (abridged)

There are now evidences that the cosmological constant $\Lambda$ has a non-zero positive value. Alternative scenarios to a pure cosmological constant model are provided by quintessence, an effective negative pressure fluid permeating the universe. Recent results indicate that the energy density $\rho$ and the pressure $p$ of this fluid are constrained by $-\rho\leq p \buildrel<\over\sim-0.6 \rho$. Since $p=-\rho$ is equivalent to the pure cosmological constant model, it is appropriate to analyze this particular, but important, case further. We study the linear theory of perturbations in a Friedmann-Robertson-Walker universe with a cosmological constant. We obtain the equations for the evolution of the perturbations in the fully relativistic case, for which we analyze the single-fluid and the two-fluid cases. We obtain solutions to these equations in appropriate limits. We also study the Newtonian approximation. We find that for a positive cosmological constant universe (i) the perturbations will grow slower in the relativistic regime for a two-fluid composed of dark matter and radiation, and (ii) in the Newtonian regime, the perturbations stop growing.