Despite the great observational success of the standard cosmological model some discrepancies in the inferred parameter constraints have manifested among a number of cosmological data sets. These include a tension between the expansion rate of our Cosmos as inferred from the cosmic microwave background (CMB) and as found from local measurements, the preference for an enhanced amplitude of CMB lensing, a somewhat low quadrupole moment of the CMB fluctuations as well as a preference for a lower amplitude of matter fluctuations in large-scale structure surveys than inferred from the CMB. We analyse these observational tensions under the addition of spatial curvature and a free CMB background temperature that may deviate from its locally measured value. With inclusion of these parameters, we observe a trend in the parameter constraints from cosmic microwave background and baryon acoustic oscillation data towards an open and hotter universe with larger current expansion rate, standard CMB lensing amplitudes, lower amplitude of matter fluctuations, and lower CMB quadrupole moment, consistently reducing the individual tensions among the cosmological data sets. Combining this data, we find a preference for an open and hotter universe. Finally, we briefly discuss a local void as a possible source for a deviation of the locally measured CMB temperature from its background value and as mimic of negative spatial curvature for CMB photons, which when correctly implemented may further ease parameter tensions.

To effectively exploit large-scale structure surveys, we depend on accurate and reliable predictions of non-linear cosmological structure formation. Tools for efficient and comprehensive computational modelling are therefore essential to perform cosmological parameter inference analyses. We present the public software package ReACT, demonstrating its capability for the fast and accurate calculation of non-linear power spectra from non-standard physics. We showcase ReACT through a series of forecast constraints on the DGP and $f(R)$ gravity models, adopting LSST-like cosmic shear power spectra. Accurate non-linear modelling with ReACT more than doubles LSST's constraining power on the $f(R)$ parameter, in contrast to an analysis that is limited to the quasi-linear regime. We find that ReACT is sufficiently robust for the inference of consistent constraints on theories beyond $\Lambda$CDM for current and ongoing surveys. With further improvement, particularly in terms of the accuracy of the non-linear $\Lambda$CDM power spectrum, ReACT can, in principle, meet the accuracy requirements for future surveys such as Euclid and LSST.

Accurate modelling of non-linear scales in galaxy clustering will be crucial for data analysis of Stage IV galaxy surveys. A selection of competing non-linear models must be made based on validation studies. We provide a comprehensive set of forecasts of two different models for the halo redshift space power spectrum, namely the commonly applied TNS model and an effective field theory of large scale structure (EFTofLSS) inspired model. Using simulation data and a least-$\chi^2$ analysis, we determine ranges of validity for the models. We then conduct an exploratory Fisher analysis using the full anisotropic power spectrum to investigate parameter degeneracies. We proceed to perform an MCMC analysis utilising the monopole, quadrupole, and hexadecapole spectra, with a restricted range of scales for the latter in order to avoid biasing our growth rate, $f$, constraint. We find that the TNS model with a Lorentzian damping and standard Eulerian perturbative modelling outperforms other variants of the TNS model. Our MCMC analysis finds that the EFTofLSS-based model may provide tighter marginalised constraints on $f$ at $z=0.5$ and $z=1$ than the TNS model, despite having additional nuisance parameters. However this depends on the range of scales used as well as the fiducial values and priors on the EFT nuisance parameters. Finally, we extend previous work to provide a consistent comparison between the Fisher matrix and MCMC forecasts using the multipole expansion formalism, and find good agreement between them.

Future large-scale structure surveys will measure three-point statistics with high statistical significance. This will offer significant improvements on our understanding of gravity, provided we can model these statistics accurately. We assess the performance of several schemes for theoretical modelling of the matter bispectrum, including halo-model based approaches and fitting formulae. We compare the model predictions against N-body simulations, considering scales up to $k_{\rm max} = 4 h/{\rm Mpc}$, well into non-linear regime of structure formation. Focusing on the equilateral configuration, we conduct this analysis for three theories of gravity: general relativity, $f(R)$ gravity, and the DGP braneworld model. Additionally, we compute the lensing convergence bispectrum for these models. We find that all current modelling prescriptions in modified gravity, in particular for theories with scale-dependent linear growth, fail to attain the accuracy required by the precision of the Stage IV surveys such as \emph{Euclid}. Among these models, we find that a halo-model corrected fitting formula achieves the best overall performance.

Constraints on gravity and cosmology will greatly benefit from performing joint clustering and weak lensing analyses on large-scale structure data sets. Utilising non-linear information coming from small physical scales can greatly enhance these constraints. At the heart of these analyses is the matter power spectrum. Here we employ a simple method, dubbed "Hybrid $P_\ell(k)$", based on the Gaussian streaming model to calculate the quasi non-linear redshift space matter power spectrum multipoles from a non-linear prescription for the matter power spectrum, which remains agnostic about gravity or dark energy. We test this approach against comoving Lagrangian acceleration simulation measurements performed in GR, DGP and $f(R)$ gravity and find that our method performs comparably or better to the dark matter TNS redshift space power spectrum model. When comparing the redshift space multipoles for halos, we find that a linear bias implementation is insufficient to be competitive with the TNS-based model used in the BOSS clustering analysis, which however has additional degrees of freedom. Our approach offers many avenues for improvement in accuracy as well as further unification under the halo model.

We assess the performance of the Taruya, Nishimichi and Saito (TNS) model for the halo redshift space power spectrum, focusing on utilising mildly non-linear scales to constrain the growth rate of structure f. Using simulations with volume and number density typical of forthcoming Stage IV galaxy surveys, we determine ranges of validity for the model at redshifts z = 0.5 and z = 1. We proceed to perform a Bayesian MCMC analysis utilising the monopole, quadrupole, and hexadecapole spectra, followed by an exploratory Fisher matrix analysis. As previously noted in other forecasts as well as in real data analyses, we find that including the hexadecapole can significantly improve the constraints. However, a restricted range of scales is required for the hexadecapole in order for the growth parameter estimation to remain unbiased, limiting the improvement. We consistently quantify these effects by employing the multipole expansion formalism in both our Fisher and MCMC forecasts.

Recent studies have demonstrated that {\em secondary} non-Gaussianity induced by gravity will be detected with a high signal-to-noise (S/N) by future and even by on-going weak lensing surveys. One way to characterise such non-Gaussianity is through the detection of a non-zero three-point correlation function of the lensing convergence field, or of its harmonic transform, the bispectrum. A recent study analysed the properties of the squeezed configuration of the bispectrum, when two wavenumbers are much larger than the third one. We extend this work by estimating the amplitude of the (reduced) bispectrum in four generic configurations, i.e., {\em squeezed, equilateral, isosceles} and {\em folded}, and for four different source redshifts $z_s=0.5,1.0,1.5,2.0$, by using an ensemble of all-sky high-resolution simulations. We compare these results against theoretical predictions. We find that, while the theoretical expectations based on widely used fitting functions can predict the general trends of the reduced bispectra, a more accurate theoretical modelling will be required to analyse the next generation of all-sky weak lensing surveys. The disagreement is particularly pronounced in the squeezed limit.

We provide in depth MCMC comparisons of two different models for the halo redshift space power spectrum, namely a variant of the commonly applied Taruya-Nishimichi-Saito (TNS) model and an effective field theory of large scale structure (EFTofLSS) inspired model. Using many simulation realisations and Stage IV survey-like specifications for the covariance matrix, we check each model's range of validity by testing for bias in the recovery of the fiducial growth rate of structure formation. The robustness of the determined range of validity is then tested by performing additional MCMC analyses using higher order multipoles, a larger survey volume and a more highly biased tracer catalogue. We find that under all tests, the TNS model's range of validity remains robust and is found to be much higher than previous estimates. The EFTofLSS model fails to capture the spectra for highly biased tracers as well as becoming biased at higher wavenumbers when considering a very large survey volume. Further, we find that the marginalised constraints on $f$ for all analyses are stronger when using the TNS model.

Future galaxy clustering surveys will probe small scales where non-linearities become important. Since the number of modes accessible on intermediate to small scales is very high, having a precise model at these scales is important especially in the context of discriminating alternative cosmological models from the standard one. In the mildly non-linear regime, such models typically differ from each other, and galaxy clustering data will become very precise on these scales in the near future. As the observable quantity is the angular power spectrum in redshift space, it is important to study the effects of non-linear density and redshift space distortion (RSD) in the angular power spectrum. We compute non-linear contributions to the angular power spectrum using a flat-sky approximation that we introduce in this work, and compare the results of different perturbative approaches with $N$-body simulations. We find that the TNS perturbative approach is significantly closer to the $N$-body result than Eulerian or Lagrangian 1-loop approximations, effective field theory of large scale structure or a halofit-inspired model. However, none of these prescriptions is accurate enough to model the angular power spectrum well into the non-linear regime. In addition, for narrow redshift bins, $\Delta z \lesssim 0.01$, the angular power spectrum acquires non-linear contributions on all scales, right down to $\ell=2$, and is hence not a reliable tool at this time. To overcome this problem, we need to model non-linear RSD terms, for example as TNS does, but for a matter power spectrum that remains reasonably accurate well into the deeply non-linear regime, such as halofit.

Cosmic microwave background (CMB) lensing is an integrated effect whose kernel is greater than half the peak value in the range $1<z<5$. Measuring this effect offers a powerful tool to probe the large-scale structure of the Universe at high redshifts. With the increasing precision of ongoing CMB surveys, other statistics than the lensing power spectrum, in particular the lensing bi-spectrum, will be measured at high statistical significance. This will provide ways to improve the constraints on cosmological models and lift degeneracies. Following on an earlier paper, we test analytical predictions of the CMB lensing bi-spectrum against full-sky lensing simulations, and discuss their validity and limitation in detail. The tree-level prediction of perturbation theory agrees with the simulation only up to $\ell\sim 200$, but the one-loop order allows capturing the simulation results up to $\ell\sim 600$. We also show that analytical predictions based on fitting formulas for the matter bi-spectrum agree reasonably well with simulation results, although the precision of the agreement depends on the configurations and scales considered. For instance, the agreement is at the $10\%$-level for the equilateral configuration at multipoles up to $\ell\sim2000$, but the difference in the squeezed limit raises to more than a factor of two at $\ell\sim2000$. This discrepancy appears to come from limitations in the fitting formula of the matter bi-spectrum. We also find that the analytical prediction for the post-Born correction to the bi-spectrum is in good agreement with the simulation. We conclude by discussing the bi-spectrum prediction in some theories of modified gravity.

Gravity-induced non-Gaussianity in the large-scale structure of the Universe, characterised by higher-order statistics such as the bispectrum (three-point cumulant), is expected to contain rich cosmological information. A measurement of the bispectrum will not only improve the cosmological constraints, but also give us the possibility to probe gravity on cosmological scales. In this paper, we present a framework to numerically calculate the one-loop matter bispectrum based on standard perturbation theory (SPT). This approach allows general modifications to the standard $\Lambda$CDM model to be easily implemented. We demonstrate the performance of the bispectrum calculation in three representative cases, namely the Vainshtein-screened Dvali-Gabadadze-Porrati (DGP) model, the chameleon-screened Hu-Sawicki $f(R)$ model and the phenomenological dark scattering (DS) momentum-exchange model. The predicted bispectra are then compared with measured results from a set of cosmological $N$-body simulations, and the impact of possible systematics arising from simplified or approximate treatments in the perturbative calculation is studied in detail. We find that the one-loop bispectrum calculation offers significantly more information on general screening and momentum exchange effects than the leading-order bispectrum calculation. Further, the accuracy of the one-loop prediction is shown to be comparable to non-linear fitting formulas over a wide range of wavenumbers ($k\lesssim0.3\,h$$\mbox{Mpc}^{-1}$) even at lower redshifts, $z\lesssim 1$.

We compare analytical computations with numerical simulations for dark-matter clustering, in general relativity and in the normal branch of DGP gravity (nDGP). Our analytical frameword is the Effective Field Theory of Large-Scale Structure (EFTofLSS), which we use to compute the one-loop dark-matter power spectrum, including the resummation of infrared bulk displacement effects. We compare this to a set of 20 COLA simulations at redshifts $z = 0$, $z=0.5$, and $z =1$, and fit the free parameter of the EFTofLSS, called the speed of sound, in both $\Lambda$CDM and nDGP at each redshift. At one-loop at $z = 0$, the reach of the EFTofLSS is $k_{\rm reach}\approx 0.14 \, h { \rm Mpc^{-1}}$ for both $\Lambda$CDM and nDGP. Along the way, we compare two different infrared resummation schemes and two different treatments of the time dependence of the perturbative expansion, concluding that they agree to approximately $1\%$ over the scales of interest. Finally, we use the ratio of the COLA power spectra to make a precision measurement of the difference between the speeds of sound in $\Lambda$CDM and nDGP, and verify that this is proportional to the modification of the linear coupling constant of the Poisson equation.

We investigate the capabilities of perturbation theory in capturing non-linear effects of dark energy. We test constant and evolving $\omega$ models, as well as models involving momentum exchange between dark energy and dark matter. Specifically, we compare perturbative predictions at 1-loop level against N-body results for four non-standard equations of state as well as varying degrees of momentum exchange between dark energy and dark matter. The interaction is modelled phenomenologically using a time dependent drag term in the Euler equation. We make comparisons at the level of the matter power spectrum and the redshift space monopole and quadrupole. The multipoles are modelled using the Taruya, Nishimichi and Saito (TNS) redshift space spectrum. We find perturbation theory does very well in capturing non-linear effects coming from dark sector interaction. We isolate and quantify the 1-loop contribution coming from the interaction and from the non-standard equation of state. We find the interaction parameter $\xi$ amplifies scale dependent signatures in the range of scales considered. Non-standard equations of state also give scale dependent signatures within this same regime. In redshift space the match with N-body is improved at smaller scales by the addition of the TNS free parameter $\sigma_v$. To quantify the importance of modelling the interaction, we create mock data sets for varying values of $\xi$ using perturbation theory. This data is given errors typical of Stage IV surveys. We then perform a likelihood analysis using the first two multipoles on these sets and a $\xi=0$ modelling, ignoring the interaction. We find the fiducial growth parameter $f$ is generally recovered even for very large values of $\xi$ both at $z=0.5$ and $z=1$. The $\xi=0$ modelling is most biased in its estimation of $f$ for the phantom $\omega=-1.1$ case.

Modifications to gravity can provide attractive alternatives to the dark components of the standard model of cosmology. These modifications to general relativity (GR) must be hidden at small scales where theory is well tested, and so one naturally looks to the large scales in order to detect any deviations from GR. One particularly promising avenue in testing gravity at cosmological scales is within the anisotropy of galaxy clustering in redshift space. This thesis presents a framework for consistently constructing large scale structure observables in redshift space for gravitational theories that include an additional scalar degree of freedom, specifically, the Horndeski class of theories with a generalized potential term. The relevance of such a framework in the context of next generation spectroscopic surveys is then investigated using N-body simulations. The thesis concludes with ongoing and recently completed extensions to this framework, including interacting dark energy models and the effective field theory of large scale structure.

We extend our previous redshift space power spectrum code to the redshift space correlation function. Here we focus on the Gaussian Streaming Model (GSM). Again, the code accommodates a wide range of modified gravity and dark energy models. For the non-linear real space correlation function used in the GSM, we use the Fourier transform of the RegPT 1-loop matter power spectrum. We compare predictions of the GSM for a Vainshtein screened and Chameleon screened model as well as GR. These predictions are compared to the Fourier transform of the Taruya, Nishimichi and Saito (TNS) redshift space power spectrum model which is fit to N-Body data. We find very good agreement between the Fourier transform of the TNS model and the GSM predictions, with $\leq 6\%$ deviations in the first two correlation function multipoles for all models for separations in $50$Mpc$/h \leq s \leq 180$Mpc/$h$. Excellent agreement is found in the differences between the modified gravity and GR multipole predictions for both approaches to the redshift space correlation function, highlighting their matched ability in picking up deviations from GR. We elucidate the timeliness of such non-standard templates at the dawn of stage-IV surveys and discuss necessary preparations and extensions needed for upcoming high quality data.

We elucidate the importance of the consistent treatment of gravity-model specific non-linearities when estimating the growth of cosmological structures from redshift space distortions (RSD). Within the context of standard perturbation theory (SPT), we compare the predictions of two theoretical templates with redshift space data from COLA (COmoving Lagrangian Acceleration) simulations in the normal branch of DGP gravity (nDGP) and General Relativity (GR). Using COLA for these comparisons is validated using a suite of full N-body simulations for the same theories. The two theoretical templates correspond to the standard general relativistic perturbation equations and those same equations modelled within nDGP. Gravitational clustering non-linear effects are accounted for by modelling the power spectrum up to one loop order and redshift space clustering anisotropy is modelled using the Taruya, Nishimichi and Saito (TNS) RSD model. Using this approach, we attempt to recover the simulation's fiducial logarithmic growth parameter $f$. By assigning the simulation data with errors representing an idealised survey with a volume of $10\mbox{Gpc}^3/h^3$, we find the GR template is unable to recover fiducial $f$ to within 1$\sigma$ at $z=1$ when we match the data up to $k_{\rm max}=0.195h$/Mpc. On the other hand, the DGP template recovers the fiducial value within $1\sigma$. Further, we conduct the same analysis for sets of mock data generated for generalised models of modified gravity using SPT, where again we analyse the GR template's ability to recover the fiducial value. We find that for models with enhanced gravitational non-linearity, the theoretical bias of the GR template becomes significant for stage IV surveys. Thus, we show that for the future large data volume galaxy surveys, the self-consistent modelling of non-GR gravity scenarios will be crucial in constraining theory parameters.

We use a set of N-body simulations employing a modified gravity (MG) model with Vainshtein screening to study matter and halo hierarchical clustering. As test-case scenarios we consider two normal branch Dvali-Gabadadze-Porrati (nDGP) gravity models with mild and strong growth rate enhancement. We study higher-order correlation functions $\xi_n(R)$ up to $n=9$ and associated hierarchical amplitudes $S_n(R)\equiv\xi_n(R)/\sigma(R)^{2n-2}$. We find that the matter PDFs are strongly affected by the fifth-force on scales up to $50h^{-1}$Mpc, and the deviations from GR are maximised at $z=0$. For reduced cumulants $S_n$, we find that at small scales $R\leq10h^{-1}$Mpc the MG is characterised by lower values, with the deviation growing from $7\%$ in the reduced skewness up to even $40\%$ in $S_5$. To study the halo clustering we use a simple abundance matching and divide haloes into thee fixed number density samples. The halo two-point functions are weakly affected, with a relative boost of the order of a few percent appearing only at the smallest pair separations ($r\leq 5h^{-1}$Mpc). In contrast, we find a strong MG signal in $S_n(R)$'s, which are enhanced compared to GR. The strong model exhibits a $>3\sigma$ level signal at various scales for all halo samples and in all cumulants. In this context, we find that the reduced kurtosis to be an especially promising cosmological probe of MG. Even the mild nDGP model leaves a $3\sigma$ imprint at small scales $R\leq3h^{-1}$Mpc, while the stronger model deviates from a GR-signature at nearly all scales with a significance of $>5\sigma$. Since the signal is persistent in all halo samples and over a range of scales, we advocate that the reduced kurtosis estimated from galaxy catalogues can potentially constitute a strong MG-model discriminatory as well as GR self-consistency test.

We develop a code to produce the power spectrum in redshift space based on standard perturbation theory (SPT) at 1-loop order. The code can be applied to a wide range of modified gravity and dark energy models using a recently proposed numerical method by A.Taruya. This includes Horndeski's theory with a general potential, which accommodates both chameleon and Vainshtein screening mechanisms and provides a non-linear extension of the effective theory of dark energy up to third order. Focus is on a recent non-linear model of the redshift space power spectrum which has been shown to model the anisotropy very well at relevant scales for the SPT framework, as well as capturing relevant non-linear effects typical of modified gravity theories. We provide consistency checks of the code against established results and elucidate it's application within the light of upcoming high precision RSD data.

Using a Gyrator-II Loop antenna tuned at 19.0Khz, we monitored the meteor shower during 17-24th November, 2002. We observe the primary peak at 3h58m (UT) on the 19th of November, 2002. We distinctly observed several `beadlike' and `exponentially dropping' signals. The `beadlike' signals were more in abundance on the 18th of November, 2002, one day prior to the actual encounter.

We report our first results of monitoring sudden ionospheric disturbances (SID). We present data taken continuously for two weeks during 20th Sept. 2002 and 4th Oct. 2002. We compare the effects of solar flares of the VLF signal with those obtained by GOES Satellite of NASA monitoring in X-rays and found excellent agreement.