The aim of this thesis is to question some of the basic assumptions that go into building the $\Lambda$CDM model of our universe. The assumptions we focus on are the initial conditions of the universe, the fundamental forces in the universe on large scales and the approximations made in analysing cosmological data. For each of the assumptions we outline the theoretical understanding behind them, the current methods used to study them and how they can be improved and finally we also perform numerical analysis to quantify the novel solutions/methods we propose to extend the previous assumptions.

In the standard big bang model the universe starts in a radiation dominated era, where the gravitational perturbations are described by second order differential equations, which will generally have two orthogonal set of solutions. One is the so called {\it growing(cosine)} mode and the other is the {\it decaying(sine)} mode, where the nomenclature is derived from their behaviour on super-horizon(sub-horizon) scales. The decaying mode is qualitatively different to the growing mode of adiabatic perturbations as it evolves with time on \emph{super-horizon} scales. The time dependence of this mode on super-horizon scales is analysed in both the synchronous gauge and the Newtonian gauge to understand the true gauge invariant behaviour of these modes. We then explore constraints on the amplitude of this mode on scales between $k \sim 10^{-5}$ Mpc$^{-1}$ and $k \sim 10^{-1}$ Mpc$^{-1}$ using the temperature and polarization anisotropies from the cosmic microwave background, by computing the Fisher information. Binning the primordial power non-parametrically into 100 bins, we find that the decaying modes are constrained at comparable variance as the growing modes on scales smaller than the horizon today using temperature anisotropies. Adding polrisation data makes the decaying mode more constrained. The decaying mode amplitude is thus constrained by $\sim 1/l$ of the growing mode. On super-horizon scales, the growing mode is poorly constrained, while the decaying mode cannot substantially exceed the scale-invariant amplitude. This interpretation differs substantially from the past literature, where the constraints were quoted in gauge-dependent variables, and resulted in illusionary tight super-horizon decaying mode constraints. The results presented here can generally be used to non-parametrically constrain any model of the early universe.

Many models of the early universe predict that there should be primordial tensor perturbations. These leave an imprint into the temperature and polarisation anisotropies of the cosmic microwave background (CMB). The differential equation describing the primordial tensor perturbations is a second order differential equation and thus has two solutions. Canonically, the decaying solution of this equation in radiation domination is dropped as it diverges at early times and on superhorizon scales while it is then suppressed at late times. Furthermore, if there is an inflationary phase prior to the radiation domination phase, the amplitude of the decaying mode will also be highly suppressed as it enters the radiation phase, thus its effect will be negligible. In this study we remain agnostic to the early universe models describing pre-radiation domination physics and allow this mode to be present and see what effect it has on the CMB anisotropies. We find that the decaying mode, if normalised at the same time on subhorizon scales as the growing mode leaves an imprint on the CMB anisotropies that is identical to the growing mode. Contrary to expectation, on large scales both modes are poorly constrained for a scale invariant spectrum, and the apparent divergence of the decaying mode does not lead to a divergent physical observable. Quantitatively, the decaying mode can be more constrained both from temperature and polarisation anisotropies. We use a model independent, non-parametric, approach to constrain both of these primordial tensor perturbations using the temperature and polarisation anisotropies. We find that both modes are best constrained at the reionisation and recombination bumps and crucially, at the reionisation bump the decaying mode can be distinguished from the growing mode.

Cross-correlating two different types of galaxy gives rise to parity breaking in the correlation function that derives from differences in the galaxies' properties and environments. This is typically associated with a difference in galaxy bias, describing the relation between galaxy number density and dark matter density, although observational effects such as magnification bias also play a role. In this paper we show that the presence of a screened fifth force adds additional degrees of freedom to the correlation function, describing the effective coupling of the force to the two galaxy populations. These are also properties of the galaxies' environments, but with different dependence in general to galaxy bias. We derive the parity-breaking correlation function analytically as a function of fifth-force strength and the two populations' fifth-force charges, and explore the result numerically using Hu-Sawicki $f(R)$ as a toy model of chameleon screening. We find that screening gives rise to an octopole, which, in the absence of magnification bias, is not present in any gravity theory without screening and thus is a qualitatively distinct signature. The modification to the dipole and octopole can be $\mathcal{O}(10\%)$ and $\mathcal{O}(100\%)$ respectively at redshift $z \gtrsim 0.5$ due to screening, but decreases towards lower redshift. The change in the background power spectrum in $f(R)$ theories induces a change in the dipole of roughly the same size, but dominant to the effect of screening at low $z$. While current data is insufficient to measure the parity-breaking dipole or octopole to the precision required to test these models, future surveys such as DESI, Euclid and SKA have the potential to probe screened fifth forces through the dipole.

Cosmological large-scale structure analyses based on two-point correlation functions often assume a Gaussian likelihood function with a fixed covariance matrix. We study the impact on cosmological parameter estimation of ignoring the parameter dependence of this covariance matrix, focusing on the particular case of joint weak-lensing and galaxy clustering analyses. Using a Fisher matrix formalism (calibrated against exact likelihood evaluation in particular simple cases), we quantify the effect of using a parameter dependent covariance matrix on both the bias and variance of the parameters. We confirm that the approximation of a parameter-independent covariance matrix is exceptionally good in all realistic scenarios. The information content in the covariance matrix (in comparison with the two point functions themselves) does not change with the fractional sky coverage. Therefore the increase in information due to the parameter dependent covariance matrix becomes negligible as the number of modes increases. Even for surveys covering less than $1\%$ of the sky, this effect only causes a bias of up to ${\cal O}(10\%)$ of the statistical uncertainties, with a misestimation of the parameter uncertainties at the same level or lower. The effect will only be smaller with future large-area surveys. Thus for most analyses the effect of a parameter-dependent covariance matrix can be ignored both in terms of the accuracy and precision of the recovered cosmological constraints.

When a supernova explodes, most of its energy is released in a shell of relativistic neutirnos which changes the surrounding geometry. We calculate the potentially observable responses to such a change in both pulsar scintillation and conventional interferometers. In both cases, the responses are permanent changes due to such a transient event. This is by-definition a memory effect. In addition to the transverse component in the usual gravitational memory (Christodolou effect) effect, it also has a longitudinal component. Furthermore it is different from the Christodolou effect as the transverse component of this memory effect also has a term that grows with time.

We show that when a supernova explodes, a nearby pulsar signal goes through a very specific change. The observed period first changes smoothly, then is followed by a sudden change in the time derivative. A stable millisecond pulsar can allow us to measure such an effect. This may improve our measurement of the total energy released in neutrinos and also the orientation of the supernova-pulsar system.