Persistent tension between low-redshift observations and the Cosmic Microwave Background radiation (CMB), in terms of two fundamental distance scales set by the sound horizon $r_d$ and the Hubble constant $H_0$, suggests new physics beyond the Standard Model or residual systematics. We examine recently updated distance calibrations from Cepheids, gravitational lensing time-delay observations, and the Tip of the Red Giant Branch. Calibrating the Baryon Acoustic Oscillations (BAO) and Type Ia supernovae with combinations of the distance indicators, we obtain a joint and self-consistent measurement of $H_0$ and $r_d$ at low redshift, independent of cosmological models and CMB inference. In an attempt to alleviate the tension between late-time and CMB-based measurements, we consider four extensions of the standard $\Lambda$CDM model. The sound horizon from our different measurements is $r_d=(137\pm3^{stat.}\pm2^{syst.})$~Mpc. Depending on the adopted distance indicators, the $combined$ tension in $H_0$ and $r_d$ ranges between 2.3 and 5.1 $\sigma$. We find that modifications of $\Lambda$CDM that change the physics after recombination fail to solve the problem, for the reason that they only resolve the tension in $H_0$, while the tension in $r_d$ remains unchanged. Pre-recombination extensions (with early dark energy or the effective number of neutrinos $\rm{N}_{\rm{eff}}=3.24 \pm 0.16$) are allowed by the data, unless the calibration from Cepheids is included. Results from time-delay lenses are consistent with those from distance-ladder calibrations and point to a discrepancy between absolute distance scales measured from the CMB (assuming the standard cosmological model) and late-time observations. New proposals to resolve this tension should be examined with respect to reconciling not only the Hubble constant but also the sound horizon derived from the CMB and other cosmological probes.

We present an algorithm for inferring the dynamical mass of galaxy clusters directly from their respective phase-space distributions, i.e. the observed line-of-sight velocities and projected distances of galaxies from the cluster centre. Our method employs neural flows, a deep neural network capable of learning arbitrary high-dimensional probability distributions, and inherently accounts, to an adequate extent, for the presence of interloper galaxies which are not bounded to a given cluster, the primary contaminant of dynamical mass measurements. We validate and showcase the performance of our neural flow approach to robustly infer the dynamical mass of clusters from a realistic mock cluster catalogue. A key aspect of our novel algorithm is that it yields the probability density function of the mass of a particular cluster, thereby providing a principled way of quantifying uncertainties, in contrast to conventional machine learning approaches. The neural network mass predictions, when applied to a contaminated catalogue with interlopers, have a logarithmic residual scatter which goes down to 0.043 dex for the most massive clusters. This is nearly an order of magnitude improvement over the classical cluster mass scaling relation with the velocity dispersion, and outperforms recently proposed machine learning approaches. We also apply our neural flow mass estimator to a compilation of galaxy observations of some well-studied clusters with robust dynamical mass estimates, further substantiating the efficacy of our algorithm. Such sophisticated approaches would undoubtedly be relevant for robust and efficient dynamical mass inference from upcoming surveys covering unprecedented volumes of the sky.

The matter sound horizon can be inferred from the cosmic microwave background within the Standard Model. Independent direct measurements of the sound horizon are then a probe of possible deviations from the Standard Model. We aim at measuring the sound horizon $r_s$ from low-redshift indicators, which are completely independent of CMB inference. We used the measured product $H(z)r_s$ from baryon acoustic oscillations (BAO) together with supernovae~\textsc{I}a to constrain $H(z)/H_{0}$ and time-delay lenses analysed by the H0LiCOW collaboration to anchor cosmological distances ($\propto H_{0}^{-1}$). {Additionally, we investigated the influence of adding a sample of quasars with higher redshift with standardisable UV-Xray luminosity distances. We adopted polynomial expansions in $H(z)$ or in comoving distances} so that our inference was completely independent of any cosmological model on which the expansion history might be based. Our measurements are independent of Cepheids and systematics from peculiar motions {to within percent-level accuracy.} The inferred sound horizon $r_s$ varies between $(133 \pm 8)$~Mpc and $(138 \pm 5)$~Mpc across different models. The discrepancy with CMB measurements is robust against model choice. Statistical uncertainties are comparable to systematics. The combination of time-delay lenses, supernovae, and BAO yields a distance ladder that is independent of cosmology (and of Cepheid calibration) and a measurement of $r_s $ that is independent of the CMB. These cosmographic measurements are then a competitive test of the Standard Model, regardless of the hypotheses on which the cosmology is based.

We present the discovery and optical follow-up of the faintest supernova-like transient known. The event (SN 2019gsc) was discovered in a star-forming host at 53\,Mpc by ATLAS. A detailed multi-colour light curve was gathered with Pan-STARRS1 and follow-up spectroscopy was obtained with the NOT and Gemini-North. The spectra near maximum light show narrow features at low velocities of 3000 to 4000 km s$^{-1}$, similar to the extremely low luminosity SNe 2010ae and 2008ha, and the light curve displays a similar fast decline (\dmr $0.91 \pm 0.10$ mag). SNe 2010ae and 2008ha have been classified as type Iax supernovae, and together the three either make up a distinct physical class of their own or are at the extreme low luminosity end of this diverse supernova population. The bolometric light curve is consistent with a low kinetic energy of explosion ($E_{\rm k} \sim 10^{49}$ erg s$^{-1}$), a modest ejected mass ($M_{\rm ej} \sim 0.2$ \msol) and radioactive powering by $^{56}$Ni ($M_{\rm Ni} \sim 2 \times 10^{-3}$ \msol). The spectra are quite well reproduced with radiative transfer models (TARDIS) and a composition dominated by carbon, oxygen, magnesium, silicon and sulphur. Remarkably, all three of these extreme Iax events are in similar low-metallicity star-forming environments. The combination of the observational constraints for all three may be best explained by deflagrations of near $M_{\rm Ch}$ hybrid carbon-oxygen-neon white dwarfs which have short evolutionary pathways to formation.

We study the dynamics of small inhomogeneities in an expanding universe collapsing to form bound structures using full solutions of the Einstein-Vlasov (N-body) equations. We compare these to standard Newtonian N-body solutions using quantities defined with respect to fiducial observers in order to bound relativistic effects. We focus on simplified initial conditions containing a limited range of length scales, but vary the inhomogeneities from small magnitude, where the Newtonian and general-relativistic calculations agree quite well, to large magnitude, where the background metric receives an order one correction. For large inhomogeneities, we find that the collapse of overdensities tends to happen faster in Newtonian calculations relative to fully general-relativistic ones. Even in this extreme regime, the differences in the spacetime evolution outside the regions of large gravitational potential and velocity are small. For standard cosmological values, we corroborate the robustness of Newtonian N-body simulations to model large scale perturbations and the related cosmic variance in the local expansion rate.

The Hubble constant $H_{0}$ and matter density $\Omega_{m}$ of the Universe are measured using the latest $\gamma$-ray attenuation results from Fermi-LAT and Cherenkov telescopes. This methodology is based upon the fact that the extragalactic background light supplies opacity for very high energy photons via photon-photon interaction. The amount of $\gamma$-ray attenuation along the line of sight depends on the expansion rate and matter content of the Universe. This novel strategy results in a value of $H_{0}=67.4_{-6.2}^{+6.0}$~km~s$^{-1}$~Mpc$^{-1}$ and $\Omega_{m}=0.14_{-0.07}^{+0.06}$. These estimates are independent and complementary to those based on the distance ladder, cosmic microwave background (CMB), clustering with weak lensing, and strong lensing data. We also produce a joint likelihood analysis of our results from $\gamma$ rays and these from more mature methodologies, excluding the CMB, yielding a combined value of $H_{0}=66.6\pm 1.6$~km~s$^{-1}$~Mpc$^{-1}$ and $\Omega_{m}=0.29\pm 0.02$.

We revisit the claimed tension, or lack thereof, of measured values of the Hubble-Lema\^{i}tre parameter $H_0$ from Cosmic Microwave Background (CMB) data and low-redshift indicators. Baryon Acoustic Oscillations (BAO) rely on the scale of the sound horizon at recombination $r_s$ to convert angular measurements into angular-diameter distances, so fixing $r_s$ from CMB measurements already constrains $H_0.$ If departures from concordance cosmology are to be constrained, truly independent measurements of $H_0$ are needed. We use the angular-diameter distances to three time-delay lenses from the H0LiCOW collaboration to calibrate the distance ladder, combine them with relative distances from Supernovae Ia and BAO, leaving $r_s$ completely free, and provide the inferred coefficients ($q_{0},j_{0},s_{0}$) in the polynomial expansion of H(z). We obtain $H_{0}r_{s}=(9895\pm161)$km/s and $H_0=(72\pm7)$km/s/Mpc. Combined with $H_0$ from H0LiCOW, then $r_s=(137\pm4.5)$Mpc is consistent with previous work and systematically lower than the CMB-inferred value. Our results are independent of the adopted cosmology, and removing Supernovae with z<0.1 has a negligible effect.

Strongly lensed supernovae can be detected as multiply imaged or highly magnified transients. In order to compare the performances of these two observational strategies, we calculate expected discovery rates as a function of survey depth in five grizy filters and for different classes of supernovae (Ia, IIP, IIL, Ibc and IIn). We find that detections via magnification is the only effective strategy for relatively shallow pre-LSST surveys. For survey depths about the LSST capacity, both strategies yield comparable numbers of lensed supernovae. Supernova samples from the two methods are to a large extent independent and combining them increases detection rates by about 50 per cent. While the number of lensed supernovae detectable via magnification saturates at the limiting magnitudes of LSST, detection rates of multiply imaged supernova still go up drastically at increasing survey depth. Comparing potential discovery spaces, we find that lensed supernovae found via image multiplicity exhibit longer time delays and larger image separations making them more suitable for cosmological constraints than their counterparts found via magnification. We find that the ZTF will find about 2 type Ia and 4 core-collapse lensed supernovae per year at a limiting magnitude of 20.6 in the r band. Applying a hybrid method which combines searching for highly magnified or multiply imaged transients, we find that LSST will detect 89 type Ia and 254 core-collapse lensed supernovae per year. In all cases, lensed core-collapsed supernovae will be dominated by type IIn supernovae contributing to 80 per cent of the total counts, although this prediction relies quite strongly on the adopted spectral templates for this class of supernovae. Revisiting the case of the lensed supernova iPTF16geu, we find that it is consistent within the 2\sigma contours of predicted redshifts and magnifications for the iPTF survey.

The primary difficulty in measuring dynamical masses of galaxy clusters from galaxy data lies in the separation between true cluster members from interloping galaxies along the line of sight. We study the impact of membership contamination and incompleteness on cluster mass estimates obtained with 25 commonly used techniques applied to nearly 1000 mock clusters. We show that all methods overestimate or underestimate cluster masses when applied to contaminated or incomplete galaxy samples respectively. This appears to be the main source of the intrinsic scatter in the mass scaling relation. Applying corrections based on a prior knowledge of contamination and incompleteness can reduce the scatter to the level of shot noise expected for poorly sampled clusters. We establish an empirical model quantifying the effect of imperfect membership on cluster mass estimation and discuss its universal and method-dependent features. We find that both imperfect membership and the response of the mass estimators depend on cluster mass, effectively causing a flattening of the estimated - true mass relation. Imperfect membership thus alters cluster counts determined from spectroscopic surveys, hence the cosmological parameters that depend on such counts.

In the standard approach to studying cosmological structure formation, the overall expansion of the Universe is assumed to be homogeneous, with the gravitational effect of inhomogeneities encoded entirely in a Newtonian potential. A topic of ongoing debate is to what degree this fully captures the dynamics dictated by general relativity, especially in the era of precision cosmology. To quantitatively assess this, we directly compare standard N-body Newtonian calculations to full numerical solutions of the Einstein equations, for cold matter with various magnitude initial inhomogeneities on scales comparable to the Hubble horizon. We analyze the differences in the evolution of density, luminosity distance, and other quantities defined with respect to fiducial observers. This is carried out by reconstructing the effective spacetime and matter fields dictated by the Newtonian quantities, and by taking care to distinguish effects of numerical resolution. We find that the fully general relativistic and Newtonian calculations show excellent agreement, even well into the nonlinear regime. They only notably differ in regions where the weak gravity assumption breaks down, which arise when considering extreme cases with perturbations exceeding standard values.

The historic detection of gravitational waves from a binary neutron star merger (GW170817) and its electromagnetic counterpart led to the first accurate (sub-arcsecond) localization of a gravitational-wave event. The transient was found to be $\sim$10" from the nucleus of the S0 galaxy NGC 4993. We report here the luminosity distance to this galaxy using two independent methods. (1) Based on our MUSE/VLT measurement of the heliocentric redshift ($z_{\rm helio}=0.009783\pm0.000023$) we infer the systemic recession velocity of the NGC 4993 group of galaxies in the cosmic microwave background (CMB) frame to be $v_{\rm CMB}=3231 \pm 53$ km s$^{-1}$. Using constrained cosmological simulations we estimate the line-of-sight peculiar velocity to be $v_{\rm pec}=307 \pm 230$ km s$^{-1}$, resulting in a cosmic velocity of $v_{\rm cosmic}=2924 \pm 236$ km s$^{-1}$ ($z_{\rm cosmic}=0.00980\pm 0.00079$) and a distance of $D_z=40.4\pm 3.4$ Mpc assuming a local Hubble constant of $H_0=73.24\pm 1.74$ km s$^{-1}$ Mpc$^{-1}$. (2) Using Hubble Space Telescope measurements of the effective radius (15.5" $\pm$ 1.5") and contained intensity and MUSE/VLT measurements of the velocity dispersion, we place NGC 4993 on the Fundamental Plane (FP) of E and S0 galaxies. Comparing to a frame of 10 clusters containing 226 galaxies, this yields a distance estimate of $D_{\rm FP}=44.0\pm 7.5$ Mpc. The combined redshift and FP distance is $D_{\rm NGC 4993}= 41.0\pm 3.1$ Mpc. This 'electromagnetic' distance estimate is consistent with the independent measurement of the distance to GW170817 as obtained from the gravitational-wave signal ($D_{\rm GW}= 43.8^{+2.9}_{-6.9}$ Mpc) and confirms that GW170817 occurred in NGC 4993.

With the advent of wide-field cosmological surveys, we are approaching samples of hundreds of thousands of galaxy clusters. While such large numbers will help reduce statistical uncertainties, the control of systematics in cluster masses becomes ever more crucial. Here we examine the effects of an important source of systematic uncertainty in galaxy-based cluster mass estimation techniques: the presence of significant dynamical substructure. Dynamical substructure manifests as dynamically distinct subgroups in phase-space, indicating an 'unrelaxed' state. This issue affects around a quarter of clusters in a generally selected sample. We employ a set of mock clusters whose masses have been measured homogeneously with commonly-used galaxy-based mass estimation techniques (kinematic, richness, caustic, radial methods). We use these to study how the relation between observationally estimated and true cluster mass depends on the presence of substructure, as identified by various popular diagnostics. We find that the scatter for an ensemble of clusters does not increase dramatically for clusters with dynamical substructure. However, we find a systematic bias for all methods, such that clusters with significant substructure have higher measured masses than their relaxed counterparts. This bias depends on cluster mass: the most massive clusters are largely unaffected by the presence of significant substructure, but masses are significantly overestimated for lower mass clusters, by $\sim10\%$ at $10^{14}$ and $\geq20\%$ for $\leq10^{13.5}$. The use of cluster samples with different levels of substructure can, therefore, bias certain cosmological parameters up to a level comparable to the typical uncertainties in current cosmological studies.

The standard relation between the cosmological redshift and cosmic scale factor underlies cosmological inference from virtually all kinds of cosmological observations, leading to the emergence of the LambdaCDM cosmological model. This relation is not a fundamental theory and thus observational determination of this function (redshift remapping) should be regarded as an insightful alternative to holding its standard form in analyses of cosmological data. Here we present non-parametric reconstructions of redshift remapping in dark-matter-dominated models and constraints on cosmological parameters from a joint analysis of all primary cosmological probes including the local measurement of the Hubble constant, Type Ia supernovae, baryonic acoustic oscillations (BAO), Planck observations of the cosmic microwave background (CMB) radiation (temperature power spectrum) and cosmic chronometers. The reconstructed redshift remapping points to an additional boost of redshift operating in late epoch of cosmic evolution, but affecting both low-redshift observations and the CMB. The model predicts a significant difference between the actual Hubble constant, h=0.48+/-0.02, and its local determination, h_obs=0.73+/-0.02. The ratio of these two values coincides closely with the maximum expansion rate inside voids formed in the corresponding open cosmological model with Omega_m=0.87+/-0.03, whereas the actual value of the Hubble constant implies the age of the Universe that is compatible with the Planck LambdaCDM cosmology. The new dark-matter-dominated model with redshift remapping provides excellent fits to all data and eliminates recently reported tensions between the Planck LambdaCDM cosmology, the local determination of the Hubble constant and the BAO measurements from the Ly-alpha forest of high-redshift quasars.

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 study evolution of voids in cosmological simulations using a new method for tracing voids over cosmic time. The method is based on tracking watershed basins (contiguous regions around density minima) of well developed voids at low redshift, on a regular grid of density field. It enables us to construct a robust and continuous mapping between voids at different redshifts, from initial conditions to the present time. We discuss how the new approach eliminates strong spurious effects of numerical origin when voids evolution is traced by matching voids between successive snapshots (by analogy to halo merger trees). We apply the new method to a cosmological simulation of a standard LambdaCDM cosmological model and study evolution of basic properties of typical voids (with effective radii between 6Mpc/h and 20Mpc/h at redshift z=0) such as volumes, shapes, matter density distributions and relative alignments. The final voids at low redshifts appear to retain a significant part of the configuration acquired in initial conditions. Shapes of voids evolve in a collective way which barely modifies the overall distribution of the axial ratios. The evolution appears to have a weak impact on mutual alignments of voids implying that the present state is in large part set up by the primordial density field. We present evolution of dark matter density profiles computed on iso-density surfaces which comply with the actual shapes of voids. Unlike spherical density profiles, this approach enables us to demonstrate development of theoretically predicted bucket-like shape of the final density profiles indicating a wide flat core and a sharp transition to high-density void walls.

The canonical redshift-scale factor relation, 1/a=1+z, is a key element in the standard LambdaCDM model of the big bang cosmology. Despite its fundamental role, this relation has not yet undergone any observational tests since Lemaitre and Hubble established the expansion of the Universe. It is strictly based on the assumption of the Friedmann-Lemaitre-Robertson-Walker metric describing a locally homogeneous and isotropic universe and that photons move on null geodesics of the metric. Thus any violation of this assumption, within general relativity or modified gravity, can yield a different mapping between the model redshift z=1/a-1 and the actually observed redshift z_obs, i.e. z_obs neq z. Here we perform a simple test of consistency for the standard redshift-scale factor relation by determining simultaneous observational constraints on the concordance LambdaCDM cosmological parameters and a generalized redshift mapping z=f(z_obs). Using current baryon acoustic oscillations (BAO) and Type Ia supernova (SN) data we demonstrate that the generalized redshift mapping is strongly degenerated with dark energy. Marginalization over a class of monotonic functions f(z_obs) changes substantially degeneracy between matter and dark energy density: the density parameters become anti correlated with nearly vertical axis of degeneracy. Furthermore, we show that current SN and BAO data, analysed in a framework with the generalized redshift mapping, do not constrain dark energy unless the BAO data include the measurements from the Ly-alpha forest of high-redshift quasars.

We study the internal kinematics of galaxy clusters in the region beyond the sphere of virialization. Galaxies around a virialized cluster are infalling towards the cluster centre with a non-zero mean radial velocity. We develop a new formalism for describing the dynamical state of clusters, by generalizing the standard Jeans formalism with the inclusion of the peculiar infall motions of galaxies and the Hubble expansion as well as the contributions due to background cosmology. Using empirical fits to the radial profiles of density, mean radial velocity and velocity anisotropy of both a stacked cluster-mass halo and two isolated halos of a cosmological dark matter only simulation, we verify that our generalized Jeans equation correctly predicts the radial velocity dispersion out to 4 virial radii. We find that the radial velocity dispersion inferred from the standard Jeans equation is accurate up to 2 virial radii, but overestimated by \approx 20% for the stacked halo and by \approx 40% for the isolated halos, in the range \approx 2-3 virial radii. Our model depends on the logarithmic growth rate of the virial radius (function of halo mass or concentration), which we estimate in 7 different ways, and on the departure from self-similarity of the evolution of the peculiar velocity profile in virial units.

Dark-matter halos grown in cosmological simulations appear to have central NFW-like density cusps with mean values of $d\log\rho/d\log r \approx -1$, and some dispersion, which is generally parametrized by the varying index $\alpha$ in the Einasto density profile fitting function. Non-universality in profile shapes is also seen in observed galaxy clusters and possibly dwarf galaxies. Here we show that non-universality, at any given mass scale, is an intrinsic property of DARKexp, a theoretically derived model for collisionless self-gravitating systems. We demonstrate that DARKexp - which has only one shape parameter, $\phi_0$ - fits the dispersion in profile shapes of massive simulated halos as well as observed clusters very well. DARKexp also allows for cored dark-matter profiles, such as those found for dwarf spheroidal galaxies. We provide approximate analytical relations between DARKexp $\phi_0$, Einasto $\alpha$, or the central logarithmic slope in the Dehnen-Tremaine analytical $\gamma$-models. The range in halo parameters reflects a substantial variation in the binding energies per unit mass of dark-matter halos.

Measurements of cosmological parameters via the distance-redshift relation usually rely on models that assume a homogenous universe. It is commonly presumed that the large-scale structure evident in our Universe has a negligible impact on the measurement if distances probed in observations are sufficiently large (compared to the scale of inhomogeneities) and are averaged over different directions on the sky. This presumption does not hold when considering the effect of the gravitational redshift caused by our local gravitational potential, which alters light coming from all distances and directions in the same way. Despite its small magnitude, this local gravitational redshift gives rise to noticeable effects in cosmological inference using SN Ia data. Assuming conservative prior knowledge of the local potential given by sampling a range of gravitational potentials at locations of Milky-Way-like galaxies identified in cosmological simulations, we show that ignoring the gravitational redshift effect in a standard data analysis leads to an additional systematic error of ~1 per cent in the determination of density parameters and the dark energy equation of state. We conclude that our local gravitational field affects our cosmological inference at a level that is important in future observations aiming to achieve percent-level accuracy.

This article is the second in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our aim is to quantify the scatter, systematic bias and completeness of cluster masses derived from a diverse set of 25 galaxy-based methods using two contrasting mock galaxy catalogues based on a sophisticated halo occupation model and a semi-analytic model. Analysing 968 clusters, we find a wide range in the RMS errors in log M200c delivered by the different methods (0.18 to 1.08 dex, i.e., a factor of ~1.5 to 12), with abundance matching and richness methods providing the best results, irrespective of the input model assumptions. In addition, certain methods produce a significant number of catastrophic cases where the mass is under- or over-estimated by a factor greater than 10. Given the steeply falling high-mass end of the cluster mass function, we recommend that richness or abundance matching-based methods are used in conjunction with these methods as a sanity check for studies selecting high mass clusters. We see a stronger correlation of the recovered to input number of galaxies for both catalogues in comparison with the group/cluster mass, however, this does not guarantee that the correct member galaxies are being selected. We do not observe significantly higher scatter for either mock galaxy catalogues. Our results have implications for cosmological analyses that utilise the masses, richnesses, or abundances of clusters, which have different uncertainties when different methods are used.

The caustic technique for measuring mass profiles of galaxy clusters relies on the assumption of spherical symmetry. When applied to aspherical galaxy clusters, the method yields mass estimates affected by the cluster orientation. Here we employ mock redshift catalogues generated from cosmological simulations to study the effect of clusters intrinsic shape and surrounding filamentary structures on the caustic mass estimates. To this end, we develop a new method for removing perturbations from large-scale structures, modelled as the two-halo term, in a caustic analysis of stacked cluster data. We find that the cluster masses inferred from kinematical data of ~10^14 Msun clusters observed along the major axis are larger than masses from those observed along the minor axis by a factor of 1.7 within the virial radius, increasing to 1.8 within three virial radii. This discrepancy increases by 20% for the most massive clusters. In addition a smaller but still significant mass discrepancy arises when filamentary structures are present near a galaxy cluster. We find that the mean cluster mass from random sightlines is unbiased at all radii and their scatter ranges from 0.14 to 0.17 within one and three virial radii, with a 40% increase for the most massive clusters. We provide tables which estimate the caustic mass bias given observational constraints on the cluster orientation.

Using N-body simulations we study the phenomenon of radial orbit instability occurring in dark matter haloes of the size of a dwarf galaxy. We carried out simulations of seven spherical models, with the same standard NFW density profile but different anisotropy profiles of particle orbits. Four of them underwent instability: two with a constant positive anisotropy, one with an anisotropic core and an isotropic envelope and one with a very small isotropic core and an anisotropic envelope. Haloes affected by the instability become approximately axisymmetric and prolate, with the profile of the shortest-to-longest axis ratio increasing with radius. The lower limit for the central value of this axis ratio is 0.3 for an NFW halo. The density profiles of the haloes did not change significantly, whereas the velocity distributions became axisymmetric. The total angular momentum rose due to large-amplitude oscillations of its components perpendicular to the symmetry axis of the halo. We also studied orbits of individual particles assigning them to classical orbit families in triaxial potentials. We find that the membership of a given particle in a family depends on its initial total angular momentum and its components along the principal axes of matter distribution.

We apply two new state-of-the-art methods that model the distribution of observed tracers in projected phase space to lift the mass / velocity anisotropy (VA) degeneracy and deduce constraints on the mass profiles of galaxies, as well as their VA. We first show how a distribution function based method applied to the satellite kinematics of otherwise isolated SDSS galaxies shows convincing observational evidence of age matching: red galaxies have more concentrated dark matter (DM) halos than blue galaxies of the same stellar or halo mass. Then, applying the MAMPOSSt technique to M87 (traced by its red and blue globular clusters) we find that very cuspy DM is favored, unless we release priors on DM concentration or stellar mass (leading to unconstrained slope). For the Fornax dwarf spheroidal (traced by its metal-rich and metal-poor stars), the inner DM slope is unconstrained, with weak evidence for a core if the stellar mass is fixed. This highlights how priors are crucial for DM modeling. Finally, we find that blue GCs around M87 and metal-rich stars in Fornax have tangential outer VA.

We report estimates of the cosmological parameters $\Omega_m$ and $\Omega_{\Lambda}$ obtained using supernovae (SNe) associated with gamma-ray bursts (GRBs) at redshifts up to 0.606. Eight high-fidelity GRB-SNe with well-sampled light curves across the peak are used. We correct their peak magnitudes for a luminosity-decline rate relation to turn them into accurate standard candles with dispersion $\sigma = 0.18$ mag. We also estimate the peculiar velocity of the low-redshift host galaxy of SN 1998bw, using constrained cosmological simulations. In a flat universe, the resulting Hubble diagram leads to best-fit cosmological parameters of $(\Omega_m, \Omega_{\Lambda}) = (0.58^{+0.22}_{-0.25},0.42 ^{+0.25}_{-0.22})$. This exploratory study suggests that GRB-SNe can potentially be used as standardizable candles to high redshifts to measure distances in the universe and constrain cosmological parameters.

This paper is the first in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our primary aim is to test the performance of these cluster mass estimation techniques on a diverse set of models that will increase in complexity. We begin by providing participating methods with data from a simple model that delivers idealised clusters, enabling us to quantify the underlying scatter intrinsic to these mass estimation techniques. The mock catalogue is based on a Halo Occupation Distribution (HOD) model that assumes spherical Navarro, Frenk and White (NFW) haloes truncated at R_200, with no substructure nor colour segregation, and with isotropic, isothermal Maxwellian velocities. We find that, above 10^14 M_solar, recovered cluster masses are correlated with the true underlying cluster mass with an intrinsic scatter of typically a factor of two. Below 10^14 M_solar, the scatter rises as the number of member galaxies drops and rapidly approaches an order of magnitude. We find that richness-based methods deliver the lowest scatter, but it is not clear whether such accuracy may simply be the result of using an over-simplistic model to populate the galaxies in their haloes. Even when given the true cluster membership, large scatter is observed for the majority non-richness-based approaches, suggesting that mass reconstruction with a low number of dynamical tracers is inherently problematic.

Clusters and groups of galaxies are highly aspherical, with shapes approximated by nearly prolate ellipsoids of revolution. An equally fundamental property is the shape of these objects in velocity space which is the anisotropy of the global velocity dispersion tensor. Here we make use of kinematical data comprising around 600 nearby clusters and rich groups of galaxies from the SDSS to place constraints on the phase-space shapes of these objects, i.e. their shapes in both position and velocity space. We show that the line of sight velocity dispersion normalised by a mass dependent velocity scale correlates with the apparent elongation, with circular (elongated) clusters exhibiting an excessive (decremental) normalised velocity dispersion. This correlation holds for dynamically young or old clusters and, therefore, it originates from projecting their intrinsic phase-space shapes rather than from dynamical evolution. It signifies that clusters are preferentially prolate not only in position space, but also in velocity space. The distribution of the axial ratios in position space is found to be well approximated by a Gaussian with a mean 0.66+/-0.01 and a dispersion 0.07+/-0.008. The velocity ellipsoids representing the shapes in velocity space are more spherical, with a mean axial ratio of 0.78+/-0.03. This finding has important implications for mass measurements based on the line of sight velocity dispersion profiles in individual clusters. For typical axial ratios of the velocity ellipsoids in the analysed cluster sample, systematic errors on the mass estimates inferred from the line of sight velocity dispersions become comparable to statistical uncertainties for galaxy clusters with as few as 40 spectroscopic redshifts.

It has been shown in previous work that DARKexp, which is a theoretically derived, maximum entropy, one shape parameter model for isotropic collisionless systems, provides very good fits to simulated and observed dark-matter halos. Specifically, it fits the energy distribution, N(E), and the density profiles, including the central cusp. Here, we extend DARKexp N(E) to include the distribution in angular momentum, L^2, for spherically symmetric systems. First, we argue, based on theoretical, semi-analytical, and simulation results, that while dark-matter halos are relaxed in energy, they are not nearly as relaxed in angular momentum, which precludes using maximum entropy to uniquely derive N(E,L^2). Instead, we require that when integrating N(E,L^2) over squared angular momenta one retrieves the DARKexp N(E). Starting with a general expression for N(E,L^2) we show how the distribution of particles in L^2 is related to the shape of the velocity distribution function, VDF, and velocity anisotropy profile, \beta(r). We then demonstrate that astrophysically realistic halos, as judged by the VDF shape and \beta(r), must have linear or convex distributions in L^2, for each separate energy bin. The distribution in energy of the most bound particles must be nearly flat, and become more tilted in favor of radial orbits for less bound particles. These results are consistent with numerical simulations and represent an important step towards deriving the full distribution function for spherically symmetric dark-matter halos.

The increasing precision in the determination of the Hubble parameter has reached a per cent level at which large-scale cosmic flows induced by inhomogeneities of the matter distribution become non-negligible. Here we use large-scale cosmological N-body simulations to study statistical properties of the local Hubble parameter as measured by local observers. We show that the distribution of the local Hubble parameter depends not only on the scale of inhomogeneities, but also on how one defines the positions of observers in the cosmic web and what reference frame is used. Observers located in random dark matter haloes measure on average lower expansion rates than those at random positions in space or in the centres of cosmic voids, and this effect is stronger from the halo rest frames compared to the CMB rest frame. We compare the predictions for the local Hubble parameter with observational constraints based on type Ia supernovae (SNIa) and CMB observations. Due to cosmic variance, for observers located in random haloes we show that the Hubble constant determined from nearby SNIa may differ from that measured from the CMB by 0.8 per cent at 1sigma statistical significance. This scatter is too small to significantly alleviate a recently claimed discrepancy between current measurements assuming a flat LCDM model. However, for observers located in the centres of the largest voids permitted by the standard LCDM model, we find that Hubble constant measurements from SNIa would be biased high by 5 per cent, rendering this tension inexistent in this extreme case.

The mass measurement of galaxy clusters is an important tool for the determination of cosmological parameters describing the matter and energy content of the Universe. However, the standard methods rely on various assumptions about the shape or the level of equilibrium of the cluster. We present a novel method of measuring cluster masses. It is complementary to most of the other methods, since it only uses kinematical information from outside the virialized cluster. Our method identifies objects, as galaxy sheets or filaments, in the cluster outer region, and infers the cluster mass by modeling how the massive cluster perturbs the motion of the structures from the Hubble flow. At the same time, this technique allows to constrain the three-dimensional orientation of the detected structures with a good accuracy. We use a cosmological numerical simulation to test the method. We then apply the method to the Coma cluster, where we find two galaxy sheets, and measure the mass of Coma to be Mvir=(9.2\pm2.4)10^{14} Msol, in good agreement with previous measurements obtained with the standard methods.

The velocity anisotropy of particles inside dark matter (DM) haloes is an important physical quantity, which is required for the accurate modelling of mass profiles of galaxies and clusters of galaxies. It is typically measured using the ratio of the radial-to-tangential velocity dispersions at a given distance from the halo centre. However, this measure is insufficient to describe the dynamics of realistic haloes, which are typically quite elongated. Studying the velocity distribution in massive DM haloes in cosmological simulations, we find that in the inner parts of the haloes the local velocity ellipsoids are strongly aligned with the major axis of the halo, the alignment being stronger for more relaxed haloes. In the outer regions of the haloes, the alignment becomes gradually weaker and the orientation is more random. These two distinct regions of different degree of the alignment coincide with two characteristic regimes of the DM density profile: shallower and steeper than \rho r^{-2}. This alignment of the local velocity ellipsoids requires reinterpretation of features found in measurements based on the spherically averaged ratio of the radial-to-tangential velocity dispersions. In particular, we show that the velocity distribution in the central halo regions is highly anisotropic. For cluster-size haloes with mass 10^{14}-10^{15} h^-1 Msun, the velocity anisotropy along the major axis is nearly independent of radius and is equal to \beta=1-\sigma^2_{perp}/\sigma^2_{radial}=0.4, which is significantly larger than the previously estimated spherically averaged velocity anisotropy. The alignment of density and velocity anisotropies, and the radial trends may also have some implications for the mass modelling based on kinematical data of such objects as galaxy clusters or dwarf spheroidals, where the orbital anisotropy is a key element in an unbiased mass inference.

Recent deep photometry of the dwarf spheroidal Ursa Major II's morphology, and spectroscopy of individual stars, have provided a number of new constraints on its properties. With a velocity dispersion $\sim$6 km s$^{-1}$, and under the assumption that the galaxy is virialised, the mass-to-light ratio is found to be approaching $\sim$2000 - apparently heavily dark matter dominated. Using N-Body simulations, we demonstrate that the observed luminosity, ellipticity, irregular morphology, velocity gradient, and the velocity dispersion can be well reproduced through processes associated with tidal mass loss, and in the absence of dark matter. These results highlight the considerable uncertainty that exists in measurements of the dark matter content of Ursa Major II. The dynamics of the inner tidal tails, and tidal stream, causes the observed velocity dispersion of stars to be boosted to values of $>$5 km s$^{-1}$ ($>$20 km s$^{-1}$ at times). This effect is responsible for raising the velocity dispersion of our model to the observed values in UMaII. We test an iterative rejection technique for removing unbound stars from samples of UMaII stars whose positions on the sky, and line-of-sight velocities, are provided. We find this technique is very effective at providing an accurate bound mass from this information, and only fails when the galaxy has a bound mass less than 10$%$ of its initial mass. However when $<2%$ mass remains bound, mass overestimation by $>$3 orders of magnitude are seen. Additionally we find that mass measurements are sensitive to measurement uncertainty in line-of-sight velocities. Measurement uncertainties of 1-4 km s$^{-1}$ result in mass overestimates by a factor of $\sim$1.3-5.7.

When measuring the mass profile of any given cosmological structure through internal kinematics, the distant background density is always ignored. This trick is often refereed to as the "Jeans Swindle". Without this trick a divergent term from the background density renders the mass profile undefined, however, this trick has no formal justification. We show that when one includes the expansion of the Universe in the Jeans equation, a term appears which exactly cancels the divergent term from the background. We thereby establish a formal justification for using the Jeans Swindle.

We study the kinematics of satellites around isolated galaxies selected from the Sloan Digital Sky Survey (SDSS) spectroscopic catalog. Using a model of the phase-space density previously measured for the halos of LCDM dark matter cosmological simulations, we determine the properties of the halo mass distribution and the orbital anisotropy of the satellites as a function of the colour-based morphological type and the stellar mass of the central host galaxy. We place constraints on the halo mass and the concentration parameter of dark matter and the satellite number density profiles. We obtain a concentration-mass relation for galactic dark matter haloes that is consistent with predictions of a standard LCDM cosmological model. At given halo or stellar mass, red galaxies have more concentrated halos than their blue counterparts. The fraction of dark matter within a few effective radii is minimal for 11.25<log M_star<11.5. The number density profile of the satellites appears to be shallower than of dark matter, with the scale radius typically 60 per cent larger than of dark matter. The orbital anisotropy around red hosts exhibits a mild excess of radial motions, in agreement with the typical anisotropy profiles found in cosmological simulations, whereas blue galaxies are found to be consistent with an isotropic velocity distribution. Our new constraints on the halo masses of galaxies are used to provide analytic approximations of the halo-to-stellar mass relation for red and blue galaxies.

Measurements of galaxy cluster kinematics are important in understanding the dynamical state and evolution of clusters of galaxies, as well as constraining cosmological models. While it is well established that clusters exhibit non-spherical geometries, evident in the distribution of galaxies on the sky, azimuthal variations of galaxy kinematics within clusters have yet to be observed. Here we measure the azimuthal dependence of the line-of-sight velocity dispersion profile in a stacked sample of 1743 galaxy clusters from the Sloan Digital Sky Survey (SDSS). The clusters are drawn from the SDSS DR8 redMaPPer catalog. We find that the line-of-sight velocity dispersion of galaxies lying along the major axis of the central galaxy is larger than those that lie along the minor axis. This is the first observational detection of anisotropic kinematics of galaxies in clusters. We show that the result is consistent with predictions from numerical simulations. Furthermore we find that the degree of projected anisotropy is strongly dependent on the line-of-sight orientation of the galaxy cluster, opening new possibilities for assessing systematics in optical cluster finding.

The theoretical framework of cosmology is mainly defined by gravity, of which general relativity is the current model. Recent tests of general relativity within the \Lambda Cold Dark Matter (CDM) model have found a concordance between predictions and the observations of the growth rate and clustering of the cosmic web. General relativity has not hitherto been tested on cosmological scales independent of the assumptions of the \Lambda CDM model. Here we report observation of the gravitational redshift of light coming from galaxies in clusters at the 99 per cent confidence level, based upon archival data. The measurement agrees with the predictions of general relativity and its modification created to explain cosmic acceleration without the need for dark energy (f(R) theory), but is inconsistent with alternative models designed to avoid the presence of dark matter.

We compare the DARKexp differential energy distribution, N(E) \propto \exp(\phi_0-E)-1, obtained from statistical mechanical considerations, to the results of N-body simulations of dark matter halos. We first demonstrate that if DARKexp halos had anisotropic velocity distributions similar to those of N-body simulated halos, their density and energy distributions could not be distinguished from those of isotropic DARKexp halos. We next carry out the comparison in two ways, using (1) the actual energy distribution extracted from simulations, and (2) N-body fitting formula for the density distribution as well as N(E) computed from the density using the isotropic Eddington formula. Both the methods independently agree that DARKexp N(E) with \phi_0\approx 4-5 is an excellent match to N-body N(E). Our results suggest (but do not prove) that statistical mechanical principles of maximum entropy can be used to explain the equilibrated final product of N-body simulations.

We analyze kinematic data of 41 nearby (z<0.1) relaxed galaxy clusters in terms of the projected phase-space density using a phenomenological, fully anisotropic model of the distribution function. We apply the Markov Chain Monte Carlo approach to place constraints on total mass distribution approximated by the universal NFW profile and the profile of the anisotropy of galaxy orbits. We find the normalization of the mean mass-concentration relation is c=6.9_{-0.7}^{+0.6} at the virial mass M_v=5x10^{14}M_sun. Assuming a one-to-one correspondence between sigma_8 and the normalization of the mass-concentration relation in the framework of the concordance model we estimate the normalization of the linear power spectrum to be sigma_8=0.91_{-0.08}^{+0.07}. Our constraints on the parameters of the mass profile are compared with estimates from other methods. We show that galaxy orbits are isotropic at the cluster centres (with the mean ratio of the radial-to-tangential velocity dispersions sigma_r/sigma_theta=0.97+/-0.04) and radially anisotropic at the virial sphere (with the mean ratio sigma_r/sigma_theta=1.75^{+0.23}_{-0.19}). Although the value of the central anisotropy appears to be universal, the anisotropy at the virial radius differs between clusters within the range 1<(sigma_r/sigma_theta)<2. Utilizing the Bautz-Morgan morphological classification and information on the prominence of a cool core we select two subsamples of galaxy clusters corresponding to less and more advanced evolutionary states. It is demonstrated that less evolved clusters have shallower mass profiles and their galaxy orbits are more radially biased at the virial sphere. This property is consistent with the expected evolution of the mass profiles as well as with the observed orbital segregation of late and early type galaxies.

We present a new method of constraining the mass and velocity anisotropy profiles of galaxy clusters from kinematic data. The method is based on a model of the phase space density which allows the anisotropy to vary with radius between two asymptotic values. The characteristic scale of transition between these asymptotes is fixed and tuned to a typical anisotropy profile resulting from cosmological simulations. The model is parametrized by two values of anisotropy, at the centre of the cluster and at infinity, and two parameters of the NFW density profile, the scale radius and the scale mass. In order to test the performance of the method in reconstructing the true cluster parameters we analyze mock kinematic data for 20 relaxed galaxy clusters generated from a cosmological simulation of the standard LCDM model. We use Bayesian methods of inference and the analysis is carried out following the Markov Chain Monte Carlo approach. The parameters of the mass profile are reproduced quite well, but we note that the mass is typically underestimated by 15 percent, probably due to the presence of small velocity substructures. The constraints on the anisotropy profile for a single cluster are in general barely conclusive. Although the central asymptotic value is determined accurately, the outer one is subject to significant systematic errors caused by substructures at large clustercentric distance. The anisotropy profile is much better constrained if one performs joint analysis of at least a few clusters. In this case it is possible to reproduce the radial variation of the anisotropy over two decades in radius inside the virial sphere.

We study the distribution function (DF) of dark matter particles in haloes of mass range 10^{14}--10^{15}\Msun. In the numerical part of this work we measure the DF for a sample of relaxed haloes formed in the simulation of a standard \LambdaCDM model. The DF is expressed as a function of energy E and the absolute value of the angular momentum L, a form suitable for comparison with theoretical models. By proper scaling we obtain the results that do not depend on the virial mass of the haloes. We demonstrate that the DF can be separated into energy and angular momentum components and propose a phenomenological model of the DF in the form f_{E}(E)[1+L^{2}/(2L_{0}^{2})]^{-\beta_{\infty}+\beta_{0}}L^{-2\beta_{0}}. This formulation involves three parameters describing the anisotropy profile in terms of its asymptotic values (\beta_{0} and \beta_{\infty}) and the scale of transition between them (L_{0}). The energy part f_{E}(E) is obtained via inversion of the integral for spatial density. We provide a straightforward numerical scheme for this procedure as well as a simple analytical approximation for a typical halo formed in the simulation. The DF model is extensively compared with the simulations: using the model parameters obtained from fitting the anisotropy profile, we recover the DF from the simulation as well as the profiles of the dispersion and kurtosis of radial and tangential velocities. Finally, we show that our DF model reproduces the power-law behaviour of phase space density Q=\rho(r)/\sigma^{3}(r).

Using an example of the Abell 2634 galaxy cluster we discuss the effect of contamination of kinematic data by interlopers and its impact on mass modelling. The cluster data reveal rich substructure along the line of sight. We demonstrate that it is necessary to apply a few independent methods of interloper removal in order to obtain a reliable sample of cluster members. We present results of three such methods which are commonly used in the literature and have been recently extensively tested on simulated data. Only two of them lead to consistent and reliable samples of cluster galaxies. For both of them we provide parameters of the best-fitting NFW density profile by fitting an isotropic solution of the Jeans equation to the velocity dispersion profiles.

We discuss how different approaches to selecting member stars in kinematic samples of dwarf spheroidal galaxies affect the estimates of their mass and anisotropy of stellar orbits. We demonstrate that the selection of members is an additional source of error compared to the usual uncertainties due to the sampling of velocity moments. As an example we use the kinematic data set for 202 stars in the Fornax dwarf galaxy for which we model the velocity dispersion profile and estimate the mass-to-light ratio and anisotropy assuming that mass follows light. We also show that stronger constraints on these parameters can be obtained if kurtosis of the velocity distribution is included in the analysis. Using the Besancon model of the Milky Way we demonstrate that the majority of contamination in Fornax probably comes from the Milky Way stars.

We summarize the method of mass modelling of galaxy clusters based on reproducing the dispersion and kurtosis of the projected velocity distribution of galaxies. The models are parametrized within the framework of the NFW density profile, characterized by the virial mass and concentration, together with the constant anisotropy of galaxy orbits. The use of velocity dispersion alone does not allow to constrain all the three parameters from kinematic data due to the mass-anisotropy degeneracy. The degeneracy is broken by introducing the fourth velocity moment, the kurtosis. We tested the method based on fitting both moments on mock data sets drawn from simulated dark matter haloes and showed it to reproduce reliably the properties of the haloes. The method has been applied to estimate the mass, concentration and anisotropy of more than 20 clusters which allowed us to confirm, for the first time using kinematic data, the mass-concentration relation found in N-body simulations.

A new revision of the gamma flux that we expect to detect in Imaging Atmospheric Cherenkov Telescopes (IACTs) from neutralino annihilation in the Draco dSph is presented in the context of the minimal supersymmetric standard models (MSSM) compatible with the present phenomenological and cosmological constraints, and using the dark matter (DM) density profiles compatible with the latest observations. This revision takes also into account the important effect of the Point Spread Function (PSF) of the telescope, and is valid not only for Draco but also for any other DM target. We show that this effect is crucial in the way we will observe and interpret a possible signal detection. Finally, we discuss the prospects to detect a possible gamma signal from Draco for current or planned gamma-ray experiments, i.e. MAGIC, GLAST and GAW. Even with the large astrophysical and particle physics uncertainties we find that the chances to detect a neutralino annihilation signal in Draco seem to be very scarce for current experiments. However, the prospects for future IACTs with upgraded performances (especially lower threshold energies and higher sensitivities) such as those offered by the CTA project, might be substantially better.

We study the effect of contamination by interlopers in kinematic samples of galaxy clusters. We demonstrate that without the proper removal of interlopers the inferred parameters of the mass distribution in the cluster are strongly biased towards higher mass and lower concentration. The interlopers are removed using two procedures previously shown to work most efficiently on simulated data. One is based on using the virial mass estimator and calculating the maximum velocity available to cluster members and the other relies on the ratio of the virial and projected mass estimators. We illustrate the performance of the methods in detail using the example of A576, a cluster with a strong uniform background contamination, and compare the case of A576 to 15 other clusters with different degree of contamination. We model the velocity dispersion and kurtosis profiles obtained for the cleaned data samples of these clusters solving the Jeans equations to estimate the mass, concentration and anisotropy parameter. We present the mass-concentration relation for the total sample of 22 clusters.

The aim of this paper is to study the efficiency of different approaches to interloper treatment in dynamical modelling of galaxy clusters. Using cosmological N-body simulation of standard LCDM model, we select 10 massive dark matter haloes and use their particles to emulate mock kinematic data in terms of projected galaxy positions and velocities as they would be measured by a distant observer. Taking advantage of the full 3D information available from the simulation, we select samples of interlopers defined with different criteria. The interlopers thus selected provide means to assess the efficiency of different interloper removal schemes found in the literature. We study direct methods of interloper removal based on dynamical or statistical restrictions imposed on ranges of positions and velocities available to cluster members. In determining these ranges, we use either the velocity dispersion criterion or a maximum velocity profile. We also generalize the common approaches taking into account both the position and velocity information. Another criterion is based on the dependence of the commonly used virial mass and projected mass estimators on the presence of interlopers. We find that the direct methods exclude on average 60-70 percent of unbound particles producing a sample with contamination as low as 2-4 percent. Next, we consider indirect methods of interloper treatment which are applied to the data stacked from many objects. In these approaches, interlopers are treated in a statistical way as a uniform background which modifies the distribution of cluster members. Using a Bayesian approach, we reproduce the properties of composite clusters and estimate the probability of finding an interloper as a function of distance from the object centre.

In this paper we proposed to use the group of analysis of symmetries of the dynamical system to describe the evolution of the Universe. This methods is used in searching for the unknown equation of state. It is shown that group of symmetries enforce the form of the equation of state for noninteracting scaling multifluids. We showed that symmetries give rise the equation of state in the form $p=-\Lambda+w_{1}\rho(a)+w_{2}a^{\beta}+0$ and energy density $\rho=\Lambda+\rho_{01}a^{-3(1+w)}+\rho_{02}a^{\beta}+\rho_{03}a^{-3}$, which is commonly used in cosmology. The FRW model filled with scaling fluid (called homological) is confronted with the observations of distant type Ia supernovae. We found the class of model parameters admissible by the statistical analysis of SNIa data. We showed that the model with scaling fluid fits well to supernovae data. We found that $\Omega_{\text{m},0} \simeq 0.4$ and $n \simeq -1$ ($\beta = -3n$), which can correspond to (hyper) phantom fluid, and to a high density universe. However if we assume prior that $\Omega_{\text{m},0}=0.3$ then the favoured model is close to concordance $\Lambda$CDM model. Our results predict that in the considered model with scaling fluids distant type Ia supernovae should be brighter than in $\Lambda$CDM model, while intermediate distant SNIa should be fainter than in $\Lambda$CDM model. We also investigate whether the model with scaling fluid is actually preferred by data over $\Lambda$CDM model. As a result we find from the Akaike model selection criterion prefers the model with noninteracting scaling fluid.

We study the mass distribution in six nearby (z<0.06) relaxed Abell clusters of galaxies A0262, A0496, A1060, A2199, A3158 and A3558. Given the dominance of dark matter in galaxy clusters we approximate their total density distribution by the NFW formula characterized by virial mass and concentration. We also assume that the anisotropy of galactic orbits is reasonably well described by a constant and that galaxy distribution traces that of the total density. Using the velocity and position data for 120-420 galaxies per cluster we calculate, after removal of interlopers, the profiles of the lowest-order even velocity moments, dispersion and kurtosis. We then reproduce the velocity moments by jointly fitting the moments to the solutions of the Jeans equations. Including the kurtosis in the analysis allows us to break the degeneracy between the mass distribution and anisotropy and constrain the anisotropy as well as the virial mass and concentration. The method is tested in detail on mock data extracted from N-body simulations of dark matter haloes. We find that the best-fitting galactic orbits are remarkably close to isotropic in most clusters. Using the fitted pairs of mass and concentration parameters for the six clusters we conclude that the trend of decreasing concentration for higher masses found in cosmological N-body simulations is consistent with the data. By scaling the individual cluster data by mass we combine them to create a composite cluster with 1465 galaxies and perform a similar analysis on such sample. The estimated concentration parameter then lies in the range 1.5 < c < 14 and the anisotropy parameter in the range -1.1 < \beta < 0.5 at the 95 percent confidence level.

We investigate the possibility of testing cosmological models with interaction between matter and energy sector. We assume the standard FRW model while the so called energy conservation condition is interpreted locally in terms of energy transfer. We analyze two forms of dark energy sectors: the cosmological constant and phantom field. We find a simple exact solution of the models in which energy transfer is described by a Cardassian like term in the relation of $H^{2}(z)$, where $H$ is Hubble's function and $z$ is redshift. The considered models have two additional parameters $(\Omega_{\text{int}},n)$ (apart the parameters of the $\Lambda$CDM model) which can be tested using SNIa data. In the estimation of the model parameters Riess et al.'s sample is used. We also confront the quality of statistical fits for both the $\Lambda$CDM model and the interacting models with the help of the Akaike and Bayesian informative criteria. Our conclusion from standard best fit method is that the interacting models explains the acceleration of the Universe better but they give rise to a universe with high matter density. However, using the tools of information criteria we find that the two new parameters play an insufficient role in improving the fit to SNIa data and the standard $\Lambda$CDM model is still preferred. We conclude that high precision detection of high redshift supernovae could supply data capable of justifying adoption of new parameters.

The Abell 1689 galaxy cluster has recently become a subject of intensive study. Thanks to its intermediate redshift (z=0.183) its mass distribution can be reconstructed using numerous methods including gravitational lensing, galaxy kinematics and X-ray imaging. The methods used to yield conflicting mass estimates in the past and recently the cluster mass distribution has been claimed to be in conflict with standard CDM scenarios due to rather large concentration and steep mass profile obtained from detailed studies of Broadhurst et al. using lensing. By studying in detail the kinematics of about 200 galaxies with measured redshifts in the vicinity of the cluster we show that the cluster is probably surrounded by a few structures, quite distant from each other, but aligned along the line of sight. We support our arguments by referring to cosmological N-body simulations and showing explicitly that distant, non-interacting haloes can produce entangled multi-peak line-of-sight velocity distributions similar to that in A1689. We conclude that it is difficult to estimate the cluster mass reliably from galaxy kinematics, but the value we obtain after applying a simple cut-off in velocity agrees roughly with the mass estimated from lensing. The complicated mass distribution around the cluster may however increase the uncertainty in the determination of the density profile shape obtained with weak lensing.

Using cosmological N-body simulations we study the line-of-sight velocity distribution of dark matter haloes focusing on the lowest-order even moments, dispersion and kurtosis, and their application to estimate the mass profiles of cosmological structures. For each of the ten massive haloes selected from the simulation box we determine the virial mass, concentration and the anisotropy parameter. In order to emulate observations from each halo we choose randomly 300 particles and project their velocities and positions along the line of sight and on the surface of the sky, respectively. After removing interlopers we calculate the profiles of the line-of-sight velocity moments and fit them with the solutions of the Jeans equations. The estimates of virial mass, concentration parameter and velocity anisotropy obtained in this way are in good agreement with the values found from the full 3D analysis.

Using cosmological N-body simulations we study the radial velocity distribution in dark matter haloes focusing on the lowest-order even moments, dispersion and kurtosis. We determine the properties of ten massive haloes in the simulation box approximating their density distribution by the NFW formula characterized by the virial mass and concentration. We also calculate the velocity anisotropy parameter of the haloes and find it mildly radial and increasing with distance from the halo centre. The radial velocity dispersion of the haloes shows a characteristic profile with a maximum, while the radial kurtosis profile decreases with distance starting from a value close to Gaussian near the centre. We therefore confirm that dark matter haloes possess intrinsically non-Gaussian, flat-topped velocity distributions. We find that the radial velocity moments of the simulated haloes are quite well reproduced by the solutions of the Jeans equations obtained for the halo parameters with the anisotropy measured in the simulations. We also study the radial velocity moments for a composite cluster made of ten haloes out to ten virial radii. In this region the velocity dispersion decreases systematically to reach the value of the background, while kurtosis increases from below to above the Gaussian value of 3 signifying a transition from a flat-topped to a strongly peaked velocity distribution with respect to the Gaussian, which can be interpreted as the dominance of ordered flow with a small dispersion. We illustrate the transition by showing explicitly the velocity distribution of the composite cluster in a few radial bins.