Dark matter comprised of axion-like particles (ALPs) generated by the realignment mechanism in the post-inflationary scenario leads to primordial isocurvature fluctuations. The power spectrum of these fluctuations is flat for small wave numbers, extending to scales accessible with cosmological surveys. We use the latest measurements of Cosmic Microwave Background (CMB) primary anisotropies together with CMB lensing, Baryonic Acoustic Oscillations (BAO) and Sunyaev Zel'dovich (SZ) cluster counts to measure the amplitude and tilt of the isocurvature component. We find preference for a white-noise isocurvature component in the CMB primary anisotropies; this conclusion is, however, weakened by current large-scale structure (LSS) data. Interpreting the result as a conservative upper limit on the isocurvature component, the combined bound on the ALP mass from all probes is $m_{a} \gtrsim 10^{-19}$ eV, with some dependence on how $m_{a}$ evolves with temperature. The expected sensitivity of cosmic shear and galaxy clustering from future LSS experiments and CMB lensing suggests improved bounds of $m_{a} \gtrsim 10^{-18}$-$10^{-13}$ eV, depending on scale cuts used to avoid non-linearities and the ALP mass-temperature dependence.

We assume that dark matter is comprised of axion-like particles (ALPs) generated by the realignment mechanism in the post-inflationary scenario. This leads to isocurvature fluctuations with an amplitude of order one for scales comparable to the horizon at the time when the ALP field starts oscillating. The power spectrum of these fluctuations is flat for small wave numbers, extending to scales relevant for cosmological observables. Denoting the relative isocurvature amplitude at $k_*$ = 0.05 Mpc$^{-1}$ by $f_{\rm iso}$, Planck observations of the cosmic microwave background (CMB) yield $f_{\rm iso}$ < 0.31 at the 2$\sigma$-level. This excludes the hypothesis of post-inflationary ALP dark matter with masses $m_{a}$ < 10$^{-20}-$10$^{-16}$ eV, where the range is due to details of the ALP mass-temperature dependence. Future CMB stage IV and 21-cm intensity mapping experiments may improve these limits by 1$-$2 orders of magnitude in $m_{a}$.

We study the possibility of using correlations between spatial modulations in the observed luminosity distribution of galaxies and the underlying density field as a cosmological probe. Considering redshift ranges, where magnification effects due to gravitational lensing may be neglected, we argue that the dipole part of such luminosity-density correlations traces the corresponding velocity-density signal which may thus be measured from a given galaxy redshift catalogue. Assuming an SDSS-like survey with mean density $\overline{n}$ = 0.01 ($h^{-1}$ Mpc)$^{-3}$ and effective volume $V_{\rm eff}$ = 0.2 ($h^{-1}$ Gpc)$^{3}$ at a fiducial redshift z = 0.1, we estimate that the velocity-density correlation function can be constrained with high signal-to-noise ratio $\gtrsim$ 10 on scales 10-100 Mpc. Similar conclusions apply to the monopole which is sensitive to the environmental dependence of galaxy luminosities and relevant to models of galaxy formation.

Adopting Schwarzschild's orbit-superposition technique, we construct a series of self-consistent galaxy models, embedded in the external field of galaxy clusters in the framework of Milgrom's MOdified Newtonian Dynamics. These models represent relatively massive ellipticals with a Hernquist radial profile at various distances from the cluster centre. Using $N$-body simulations, we perform a first analysis of these models and their evolution. We find that self-gravitating axisymmetric density models, even under a weak external field, lose their symmetry by instability and generally evolve to triaxial configurations. A kinematic analysis suggests that the instability originates from both box and non-classified orbits with low angular momentum. We also consider a self-consistent isolated system which is then placed in a strong external field and allowed to evolve freely. This model, just as the corresponding equilibrium model in the same external field, eventually settles to a triaxial equilibrium as well, but has a higher velocity radial anisotropy and is rounder. The presence of an external field in MOND universe generically predicts some lopsidedness of galaxy shapes.

Starting from a very accurate model for density-in-cells statistics of dark matter based on large deviation theory, a bias model for the tracer density in spheres is formulated. It adopts a mean bias relation based on a quadratic bias model to relate the log-densities of dark matter to those of mass-weighted dark haloes in real and redshift space. The validity of the parametrised bias model is established using a parametrisation-independent extraction of the bias function. This average bias model is then combined with the dark matter PDF, neglecting any scatter around it: it nevertheless yields an excellent model for densities-in-cells statistics of mass tracers that is parametrised in terms of the underlying dark matter variance and three bias parameters. The procedure is validated on measurements of both the one and two point statistics of subhalo densities in the state-of-the-art Horizon Run 4 simulation showing excellent agreement for measured dark matter variance and bias parameters. Finally, it is demonstrated that this formalism allows for a joint estimation of the nonlinear dark matter variance and the bias parameters using solely the statistics of subhaloes. Having verified that galaxy counts in hydrodynamical simulations sampled on a scale of 10 Mpc/h closely resemble those of subhaloes, this work provides important steps towards making theoretical predictions for density-in-cells statistics applicable to upcoming galaxy surveys like Euclid or WFIRST.

Observed galaxy luminosities (derived from redshifts) hold information on the large-scale peculiar velocity field in the form of spatially correlated scatter, which allows for bounds on bulk flows and the growth rate of matter density perturbations using large galaxy redshift surveys. We apply this luminosity approach to galaxies from the recent SDSS Data Release 13. Our goal is twofold. First, we take advantage of the recalibrated photometry to identify possible systematic errors relevant to our previous analysis of earlier data. Second, we seek improved constraints on the bulk flow and the normalized growth rate f$\sigma_{8}$ at z ~ 0.1. Our results confirm the robustness of our method. Bulk flow amplitudes, estimated in two redshift bins with 0.02 < z$_{1}$ < 0.07 < z$_{2}$ < 0.22, are generally smaller than in previous measurements, consistent with both the updated photometry and expectations for the $\Lambda$CDM model. The obtained growth rate, f$\sigma_{8}$ = 0.48 +/- 0.16, is larger than, but still compatible with, its previous estimate, and closer to the reference value of Planck. Rather than precision, the importance of these results is due to the fact that they follow from an independent method that relies on accurate photometry, which is a top requirement for next-generation photometric catalogs.

We assess the effect of the local large scale structure on the estimation of two-point statistics of the observed radial peculiar velocities of galaxies. A large N-body simulation is used to examine these statistics from the perspective of random observers as well as "Local Group (LG)-like" observers conditioned to reside in an environment resembling the observed universe within 20 Mpc. The local environment systematically distorts the shape and amplitude of velocity statistics with respect to ensemble-averaged measurements made by a Copernican (random) observer. The Virgo cluster has the most significant impact, introducing large systematic deviations in all the statistics. For a simple "top-hat" selection function, an idealized survey extending to $\sim 160h^{-1}\,{\rm Mpc}$ or deeper is needed to completely mitigate the effects of the local environment. Using shallower catalogues leads to systematic deviations of the order of $50$ to $200\%$ depending on the scale considered. For a flat redshift distribution similar to the one of the CosmicFlows-3 survey, the deviations are even more prominent in both the shape and amplitude at all separations considered $({\stackrel{<}{{}_\sim}} 100h^{-1}\,{\rm Mpc})$. Conclusions based on statistics calculated without taking into account the impact of the local environment should be revisited.

I revisit cosmological perturbations in Bekenstein's tensor-vector-scalar theory (TeVeS). Considering only scalar modes in the conformal Newtonian gauge, the extra degrees of freedom are expressed in a way suitable for studying modifications at the level of the metric potentials. Assuming a universe in the matter-dominated phase, I discuss the mechanism responsible for boosting structure growth, and confirm the vector field as its key ingredient. Using a semi-analytic approach, I further characterize the evolution of density perturbations and the potentials on sub- and superhorizon scales.

Spatial variations in the distribution of galaxy luminosities, estimated from redshifts as distance proxies, are correlated with the peculiar velocity field. Comparing these variations with the peculiar velocities inferred from galaxy redshift surveys is a powerful test of gravity and dark energy theories on cosmological scales. Using ~ 2 $\times$ 10$^{5}$ galaxies from the SDSS Data Release 7, we perform this test in the framework of gravitational instability to estimate the normalized growth rate of density perturbations f$\sigma_{8}$ = 0.37 +/- 0.13 at z ~ 0.1, which is in agreement with the $\Lambda$CDM scenario. This unique measurement is complementary to those obtained with more traditional methods, including clustering analysis. The estimated accuracy at z ~ 0.1 is competitive with other methods when applied to similar datasets.

Large-scale peculiar motion modulates the observed luminosity distribution of galaxies. Using about half a million SDSS galaxies, this can be harnessed to obtain bounds on peculiar velocity moments, the amplitude of the linear matter power spectrum, $\sigma_{8}$, and the growth rate of density perturbations at z ~ 0.1. Results obtained from this approach agree well with the predictions of the $\Lambda$CDM model and are consistent with the reported ~ 1% zero-point tilt in the SDSS photometry.

Peculiar motion introduces systematic variations in the observed luminosity distribution of galaxies. This allows one to constrain the cosmic peculiar velocity field from large galaxy redshift surveys. Using around half a million galaxies from the SDSS Data Release 7 at z ~ 0.1, we demonstrate the applicability of this approach to large datasets and obtain bounds on peculiar velocity moments and $\sigma_{8}$, the amplitude of the linear matter power spectrum. Our results are in good agreement with the $\Lambda$CDM model and consistent with the previously reported ~ 1% zero-point tilt in the SDSS photometry. Finally, we discuss the prospects of constraining the growth rate of density perturbations by reconstructing the full linear velocity field from the observed galaxy clustering in redshift space.

Spatial modulations in the distribution of observed luminosities (computed using redshifts) of ~ 5 $\times$ 10$^5$ galaxies from the SDSS Data Release 7, probe the cosmic peculiar velocity field out to z ~ 0.1. Allowing for luminosity evolution, the r-band luminosity function, determined via a spline-based estimator, is well represented by a Schechter form with M$^{\star}$(z) - 5log$_{10}$h = -20.52 - 1.6(z - 0.1) $\pm$ 0.05 and $\alpha^{\star}$ = -1.1 $\pm$ 0.03. Bulk flows and higher velocity moments in two redshift bins, 0.02 < z < 0.07 and 0.07 < z < 0.22, agree with the predictions of the $\Lambda$CDM model, as obtained from mock galaxy catalogs designed to match the observations. Assuming a $\Lambda$CDM model, we estimate $\sigma_{8}$ $\approx$ 1.1 $\pm$ 0.4 for the amplitude of the linear matter power spectrum, where the low accuracy is due to the limited number of galaxies. While the low-z bin is robust against coherent photometric uncertainties, the bias of results from the second bin is consistent with the ~ 1% magnitude tilt reported by the SDSS collaboration. The systematics are expected to have a significantly lower impact in future datasets with larger sky coverage and better photometric calibration.

The observed density field in redshift space is directly affected by the radial motions generated from mass fluctuations outside the volume occupied by a given galaxy redshift survey. These motions introduce redshift space anisotropies which are more pronounced at larger distances from the survey's center, thus offering clues to the nature of mass fluctuations on super-survey scales. Furthermore, we note that all estimates of the growth factor derived from redshift space distortions are based on relations which explicitly assume that the velocity field is generated by mass fluctuations inside the survey volume. This may cause uncertainties in these estimates which are on the order of a few percent.

The cosmological peculiar velocity field (deviations from the pure Hubble flow) of matter carries significant information on dark energy, dark matter and the underlying theory of gravity on large scales. Peculiar motions of galaxies introduce systematic deviations between the observed galaxy redshifts z and the corresponding cosmological redshifts z_cos. A novel method for estimating the angular power spectrum of the peculiar velocity field based on observations of galaxy redshifts and apparent magnitudes m (or equivalently fluxes) is presented. This method exploits the fact that a mean relation between z_cos and m of galaxies can be derived from all galaxies in a redshift-magnitude survey. Given a galaxy magnitude, it is shown that the z_cos(m) relation yields its cosmological redshift with a 1-sigma error of sigma_z~0.3 for a survey like Euclid (~10^9 galaxies at z<~2), and can be used to constrain the angular power spectrum of z-z_cos(m) with a high signal-to-noise ratio. At large angular separations corresponding to l<~15, we obtain significant constraints on the power spectrum of the peculiar velocity field. At 15<~l<~60, magnitude shifts in the z_cos(m) relation caused by gravitational lensing magnification dominate, allowing us to probe the line-of-sight integral of the gravitational potential. Effects related to the environmental dependence in the luminosity function can easily be computed and their contamination removed from the estimated power spectra. The amplitude of the combined velocity and lensing power spectra at z~1 can be measured with <~5% accuracy.

A unique signature of the modified Newtonian dynamics (MOND) paradigm is its peculiar behavior in the vicinity of the points where the total Newtonian acceleration exactly cancels. In the Solar System, these are the saddle points of the gravitational potential near the planets. Typically, such points are embedded into low-acceleration bubbles where modified gravity theories a` la MOND predict significant deviations from Newton's laws. As has been pointed out recently, the Earth-Sun bubble may be visited by the LISA Pathfinder spacecraft in the near future, providing a unique occasion to put these theories to a direct test. In this work, we present a high-precision model of the Solar System's gravitational potential to determine accurate positions and motions of these saddle points and study the predicted dynamical anomalies within the framework of quasi-linear MOND. Considering the expected sensitivity of the LISA Pathfinder probe, we argue that interpolation functions which exhibit a "faster" transition between the two dynamical regimes have a good chance of surviving a null result. An example of such a function is the QMOND analog of the so-called simple interpolating function which agrees well with much of the extragalactic phenomenology. We have also discovered that several of Saturn's outermost satellites periodically intersect the Saturn-Sun bubble, providing the first example of Solar System objects that regularly undergo the MOND regime.

Certain covariant theories of the modified Newtonian dynamics paradigm seem to require an additional hot dark matter (HDM) component - in the form of either heavy ordinary neutrinos or more recently light sterile neutrinos (SNs) with a mass around 11eV - to be relieved of problems ranging from cosmological scales down to intermediate ones relevant for galaxy clusters. Here we suggest using gravitational lensing by galaxy clusters to test such a marriage of neutrino HDM and modified gravity, adopting the framework of tensor-vector-scalar theory (TeVeS). Unlike conventional cold dark matter (CDM), such HDM is subject to strong phase-space constraints, which allows one to check cluster lens models inferred within the modified framework for consistency. Since the considered HDM particles cannot collapse into arbitrarily dense clumps and only form structures well above the galactic scale, systems which indicate the need for dark substructure are of particular interest. As a first example, we study the cluster lens Abell 2390 and its impressive straight arc with the help of numerical simulations. Based on our results, we outline a general and systematic approach to model cluster lenses in TeVeS which significantly reduces the calculation complexity. We further consider a simple bimodal lens configuration, capable of producing the straight arc, to demonstrate our approach. We find that such a model is marginally consistent with the hypothesis of 11eV SNs. Future work including more detailed and realistic lens models may further constrain the necessary SN distribution and help to conclusively assess this point. Cluster lenses could therefore provide an interesting discriminator between CDM and such modified gravity scenarios supplemented by SNs or other choices of HDM.

We lay out the framework to numerically study nonlinear structure formation in the context of scalar-field-coupled cold dark matter models (phiCDM models) where the scalar field phi serves as dynamical dark energy. Adopting parameters for the scalar field which yield a realistic CMB spectrum, we generate the initial conditions for our Nbody simulations, which follow the spatial distributions of dark matter and the scalar field, by solving their equations of motion using the multilevel adaptive grid technique. We show that the spatial configuration of the scalar field tracks both the voids and clusters of dark matter. The phiCDM model differs from standard LCDM at small scales with observable modifications of, e.g., the mass function of halos as well as the matter power spectrum. Nevertheless, the predictions of both models for the Hubble expansion and the CMB spectrum are virtually indistinguishable. Hence, galaxy cluster counts and weak lensing observations, which probe structure formation at small scales, are needed to falsify this class of models.

Strong gravitational lensing by galaxies in MOdified Newtonian Dynamics (MOND) has until now been restricted to spherically symmetric models. These models were able to account for the size of the Einstein ring of observed lenses, but were unable to account for double-imaged systems with collinear images, as well as four-image lenses. Non-spherical models are generally cumbersome to compute numerically in MOND, but we present here a class of analytic non-spherical models that can be applied to fit double-imaged and quadruple-imaged systems. We use them to obtain a reasonable MOND fit to ten double-imaged systems, as well as to the quadruple-imaged system Q2237+030 which is an isolated bulge-disc lens producing an Einstein cross. However, we also find five double-imaged systems and three quadruple-imaged systems for which no reasonable MOND fit can be obtained with our models. We argue that this is mostly due to the intrinsic limitation of the analytic models, even though the presence of small amounts of additional dark mass on galaxy scales in MOND is also plausible.

Intercluster filaments negligibly contribute to the weak lensing signal in general relativity (GR), $\gamma_{N}\sim 10^{-4}-10^{-3}$. In the context of relativistic modified Newtonian dynamics (MOND) introduced by Bekenstein, however, a single filament inclined by $\approx 45^\circ$ from the line of sight can cause substantial distortion of background sources pointing towards the filament's axis ($\kappa=\gamma=(1-A^{-1})/2\sim 0.01$); this is rigorous for infinitely long uniform filaments, but also qualitatively true for short filaments ($\sim 30$Mpc), and even in regions where the projected matter density of the filament is equal to zero. Since galaxies and galaxy clusters are generally embedded in filaments or are projected on such structures, this contribution complicates the interpretation of the weak lensing shear map in the context of MOND. While our analysis is of mainly theoretical interest providing order-of-magnitude estimates only, it seems safe to conclude that when modeling systems with anomalous weak lensing signals, e.g. the "bullet cluster" of Clowe et al., the "cosmic train wreck" of Abell 520 from Mahdavi et al., and the "dark clusters" of Erben et al., filamentary structures might contribute in a significant and likely complex fashion. On the other hand, our predictions of a (conceptual) difference in the weak lensing signal could, in principle, be used to falsify MOND/TeVeS and its variations.

Aims: We explore the lensing properties of asymmetric matter density distributions in Bekenstein's Tensor-Vector-Scalar theory (TeVeS). Methods: Using an iterative Fourier-based solver for the resulting non-linear scalar field equation, we numerically calculate the total gravitational potential and derive the corresponding TeVeS lensing maps. Results: Considering variations on rather small scales, we show that the lensing properties significantly depend on the lens's extent along the line of sight. Furthermore, all simulated TeVeS convergence maps strongly track the dominant baryonic components, non-linear effects, being capable of counteracting this trend, turn out to be very small. Setting up a toy model for the cluster merger 1E0657-558, we infer that TeVeS cannot explain observations without assuming an additional dark mass component in both cluster centers, which is in accordance with previous work.