The peculiar velocity field offers a unique way to probe dark matter density field on large scales at low redshifts. In this work, we have compiled a new sample of 465 peculiar velocities from low redshift $(z < 0.067)$ Type Ia supernovae. We compare the reconstructed velocity field derived from the 2M++ galaxy redshift compilation to the supernovae, the SFI++ and the 2MTF Tully-Fisher distance catalogues. We used a forward method to jointly infer the distances and the velocities of distance indicators by comparing the observations to the reconstruction. Comparison of the reconstructed peculiar velocity fields to observations allows us to infer the cosmological parameter combination $f\sigma_8$, and the bulk flow velocity arising from outside the survey volume. The residual bulk flow arising from outside the 2M++ volume is inferred to be $159^{+12}_{-11}$ km s$^{-1}$ in the direction $l=259^{\circ} \pm 4^{\circ}$ and $b= 5^{\circ}\pm 3^{\circ}$. We obtain $f\sigma_{8} = 0.401 \pm 0.017$, equivalent to $\sigma_8(\Omega_m/0.3)^{0.55}=0.775 \pm 0.033$, a $\lesssim 4\%\,$ statistical uncertainty on the value of $f\sigma_8$. Our inferred value is consistent with other low redshift results in the literature.

Most measurements of mass in Astronomy that use kinematics of stars or gas rely on assumptions of equilibrium that are often hard to verify. Instead, we develop a novel idea that uses the clustering in action space, as a probe of underlying gravitational potential: the correct potential should maximize small-scale clustering in the action space. We provide a first-principle derivation of likelihood using the two-point correlation function in action space, and test it against simulations of stellar streams. We then apply this method to the 2nd data release of Gaia, and use it to measure the radial force fraction $f_h$ and logarithmic slope $\alpha$ of dark matter halo profile. We investigate stars within 9-11 kpc and 11.5-15 kpc from Galactic centre, and find $(f_h,\alpha)= (0.391\pm 0.009, 1.835\pm 0.092) $ and $(0.351\pm 0.012,1.687\pm 0.079)$, respectively. We also confirm that the set of parameters that maximize the likelihood function do correspond to the most clustering in the action space. The best-fit circular velocity curve for Milky Way potential is consistent with past measurements (although it is $\sim$ 5-10\% lower than previous methods that use masers or globular clusters). Our work provides a clear demonstration of the full statistical power that lies in the full phase space information, relieving the need for {\it ad hoc} assumptions such as virial equilibrium, circular motion, or steam-finding algorithms.

An ambitious goal in cosmology is to forward-model the observed distribution of galaxies in the nearby Universe today from the initial conditions of large-scale structures. For practical reasons, the spatial resolution at which this can be done is necessarily limited. Consequently, one needs a mapping between the density of dark matter averaged over ~Mpc scales, and the distribution of dark matter halos (used as a proxy for galaxies) in the same region. Here we demonstrate a method for determining the halo mass distribution function by learning the tracer bias between density fields and halo catalogues using a neural bias model. The method is based on the Bayesian analysis of simple, physically motivated, neural network-like architectures, which we denote as neural physical engines, and neural density estimation. As a result, we are able to sample the initial phases of the dark matter density field whilst inferring the parameters describing the halo mass distribution function, providing a fully Bayesian interpretation of both the initial dark matter density distribution and the neural bias model. We successfully run an upgraded BORG inference using our new likelihood and neural bias model with halo catalogues derived from full N-body simulations. We notice orders of magnitude improvement in modelling compared to classical biasing techniques.

This paper presents the first derivation of the quadratic action for curvature perturbations, $\zeta$, within the framework of cuscuton gravity. We study the scalar cosmological perturbations sourced by a canonical single scalar field in the presence of cuscuton field. We identify $\zeta$ as comoving curvature with respect to the source field and we show that it retains its conservation characteristic on super horizon scales. The result provides an explicit proof that cuscuton modification of gravity around Friedmann-Lemaitre-Robertson-Walker (FLRW) metric is ghost free. We also investigate the potential development of other instabilities in cuscuton models. We find that in a large class of these models, there is no generic instability problem. However, depending on the details of slow-roll parameters, specific models may display gradient instabilities.