Combining our knowledge of halo structure and internal kinematics from cosmological dark matter simulations and the distribution of halo interlopers in projected phase space measured in cosmological galaxy simulations, we develop MAGGIE, a prior- and halo-based, probabilistic, abundance matching (AM) grouping algorithm for doubly complete subsamples (in distance and luminosity) of flux-limited samples. We test MAGGIE-L and MAGGIE-M (in which group masses are derived from AM applied to the group luminosities and stellar masses, respectively) on groups of at least three galaxies extracted from a mock Sloan Digital Sky Survey Legacy redshift survey, incorporating realistic observational errors on galaxy luminosities and stellar masses. In comparison with the optimal Friends-of-Friends (FoF) group finder, groups extracted with MAGGIE are much less likely to be secondary fragments of true groups; in primary fragments, its galaxy memberships (relative to the virial sphere of the real-space group) are much more complete and usually more reliable, and its masses are much less biased and usually with less scatter, as are its group luminosities and stellar masses (computed in MAGGIE using the membership probabilities as weights). FoF outperforms MAGGIE only for high-mass clusters: for the reliability of the galaxy population and the dispersion of its total mass. In comparison with our implementation of the Yang et al. group finder, MAGGIE reaches much higher completeness and slightly lower group fragmentation and dispersion on group total masses, luminosities and stellar masses, but slightly greater bias in the latter two and lower reliabilities. MAGGIE should therefore lead to sharper trends of environmental effects on galaxies and more accurate mass/orbit modelling.

The velocity distribution function (VDF) of dark matter (DM) halos in $\Lambda$CDM dissipationless cosmological simulations, which must be non-separable in its radial and tangential components, is still poorly known. We present the first single-parameter, non-separable, anisotropic model for the VDF in $\Lambda$CDM halos, built from an isotropic $q$-Gaussian (Tsallis) VDF of the isotropic set of dimensionless spherical velocity components (after subtraction of streaming motions), normalized by the respective velocity dispersions. We test our VDF on 90 cluster-mass halos of a dissipationless cosmological simulation. Beyond the virial radius, $r_{\rm vir}$, our model VDF adequately reproduces that measured in the simulated halos, but no $q$-Gaussian model can adequately represent the VDF within $r_{\rm vir}$, as the speed distribution function is then flatter-topped than any $q$-Gaussian can allow. Nevertheless, our VDF fits significantly better the simulations than the commonly used Maxwellian (Gaussian) distribution, at virtually all radii within $5\,r_{\rm vir}$. Within 0.4 (1) $r_{\rm vir}$, the non-Gaussianity index $q$ is (roughly) linearly related to the slope of the density profile and also to the velocity anisotropy profile. We provide a parametrization of the modulation of $q$ with radius for both the median fits and the fit of the stacked halo. At radii of a few percent of $r_{\rm vir}$, corresponding to the Solar position in the Milky Way, our best-fit VDF, although fitting better the simulations than the Gaussian one, overproduces significantly the fraction of high velocity objects, indicating that one should not blindly use these $q$-Gaussian fits to make predictions on the direct detection rate of DM particles.

We propose a fresh look at the Main Galaxy Sample of the Sloan Digital Sky Survey by packing the galaxies in stellar mass and redshift bins. We show how important it is to consider the emission-line equivalent widths, in addition to the commonly used emission-line ratios, to properly identify retired galaxies (i.e. galaxies that have stopped forming stars and are ionized by their old stellar populations) and not mistake them for galaxies with low-level nuclear activity. We find that the proportion of star-forming galaxies decreases with decreasing redshift in each mass bin, while that of retired galaxies increases. Galaxies with $M_\star > 10^{11.5} M_\odot$ have formed all their stars at redshift larger than 0.4. The population of AGN hosts is never dominant for galaxy masses larger than $10^{10} M_\odot$. We warn about the effects of stacking galaxy spectra to discuss galaxy properties. We estimate the lifetimes of active galactic nuclei (AGN) relying entirely on demographic arguments --- i.e. without any assumption on the AGN radiative properties. We find upper-limit lifetimes of about 1--5 Gyr for detectable AGN in galaxies with masses between $10^{10}$--$10^{12} M_\odot$. The lifetimes of the AGN-dominated phases are a few $10^8$ yr. Finally, we compare the star-formation histories of star-forming, AGN and retired galaxies as obtained by the spectral synthesis code STARLIGHT. Once the AGN is turned on it inhibits star formation for the next $\sim$ 0.1 Gyr in galaxies with masses around $10^{10} M_\odot$, $\sim$ 1 Gyr in galaxies with masses around $10^{11} M_\odot$.

We test the Friends-of-Friends (FoF) grouping algorithm, which depends on two linking lengths (LLs), plane-of-sky and line-of-sight (LOS), normalized to the mean nearest neighbor separation of field galaxies. We argue, on theoretical grounds, that LLs should be $b_\perp \simeq 0.11$, and $b_\parallel \simeq 1.3$ to recover 95% of all galaxies with projected radii within the virial radius $r_{200}$ and 95% of the galaxies along the LOS. We then predict that 80 to 90% of the galaxies in FoF groups should lie within their parent real-space groups (RSGs), defined within their virial spheres. We test the FoF extraction for 16x16 pairs of LLs, using subsamples of galaxies, doubly complete in distance and luminosity, of a flux-limited mock SDSS galaxy catalog. We find that massive RSGs are more prone to fragmentation, while the fragments typically have low estimated mass, with typically 30% of groups of low and intermediate estimated mass being fragments. Group merging rises drastically with estimated mass. For groups of 3 or more galaxies, galaxy completeness and reliability are both typically better than 80% (after discarding the fragments). Estimated masses of extracted groups are biased low, by up to a factor 4 at low richness, while the inefficiency of mass estimation improves from 0.85 dex to 0.2 dex when moving from low to high multiplicity groups. The optimal LLs depend on the scientific goal for the group catalog. We propose $b_\perp \simeq 0.07$, with $b_\parallel \simeq 1.1$ for studies of environmental effects, $b_\parallel \simeq 2.5$ for cosmographic studies and $b_\parallel \simeq 5$ for followups of individual groups.