The growth rate of large-scale structure provides a powerful consistency test of the standard cosmological model and a probe of possible deviations from general relativity. We use a Fisher analysis to forecast constraints on the growth rate from a combination of next-generation spectroscopic surveys. In the overlap survey volumes, we use a multi-tracer analysis to significantly reduce the effect of cosmic variance. The non-overlap individual survey volumes are included in the Fisher analysis in order to utilise the entire volume. We use the observed angular power spectrum, which naturally includes all wide-angle and lensing effects and circumvents the need for an Alcock-Paczynski correction. Cross correlations between redshift bins are included by using a novel technique to avoid computation of the sub-dominant contributions. Marginalising over the standard cosmological parameters, as well as the clustering bias in each redshift bin, we find that the precision on $\gamma$ improves on the best single-tracer precision by up to $\sim$50\%.

A new generation of cosmological experiments will spectroscopically detect the H$\alpha$ line from emission-line galaxies (ELGs) at optical/near-infrared frequencies. Other emission lines will also be present, which may come from the same H$\alpha$ sample or constitute a new galaxy sample altogether. Our goal is to assess the value, for cosmological investigation, of galaxies at $z\gtrsim2$ present in H$\alpha$ galaxy surveys and identifiable by the highly redshifted ultra-violet and optical lines - namely the OII line and the OIII doublet in combination with the H$\beta$ line. We use state-of-the-art models of luminosity functions of astrophysical spectral lines to estimate the volumetric number density of OIII+H$\beta$ and OII ELGs. We focus on a wavelength range which will be covered by planned cosmological surveys. We study the constraining power of these high-redshift galaxy samples on cosmological parameters such as the amplitude of baryon acoustic oscillations, $H(z)$, $D_A(z)$, $f\sigma_8(z)$, and $b\sigma_8(z)$ for different survey designs. We present a strong science case for extracting the OIII+H$\beta$ sample, which we consider as an independent probe of the Universe in the redshift range 2 to 3. Moreover, we show that the OII sample can be used to measure the baryon acoustic oscillations and the growth of structures above $z=3$; albeit it may be shot-noise dominated, it will nonetheless provide valuable tomographic information. Summarising, we discuss the scientific potential of a sample of galaxies which, so far, has been mainly considered as a contaminant in H$\alpha$ galaxy surveys. Our findings indicate that planed H$\alpha$ surveys should include the extraction of these oxygen-line samples in their pipeline, to enhance their scientific impact on cosmology.

Tighter constraints on measurements of primordial non-Gaussianity will allow the differentiation of inflationary scenarios. The cosmic microwave background bispectrum-the standard method of measuring the local non-Gaussianity-is limited by cosmic variance. Therefore, it is sensible to investigate measurements of non-Gaussianity using the large-scale structure. This can be done by investigating the effects of non-Gaussianity on the power spectrum on large scales. In this study we forecast the constraints on the local primordial non-Gaussianity parameter $f_{\rm NL}$ that can be obtained with future radio surveys. We utilize the multi-tracer method which reduces the effect of cosmic variance and takes advantage of the multiple radio galaxy populations which are differently biased tracers of the same underlying dark matter distribution. Improvements on previous work include the use of observational bias and halo mass estimates, updated simulations and realistic photometric redshift expectations, thus producing more realistic forecasts. Combinations of SKA simulations and radio observations were used as well as different redshift ranges and redshift bin sizes. It was found that in the most realistic case the 1 - $\sigma$ error on $f_{\rm NL}$ falls within the range 4.07 and 6.58, rivalling the tightest constraints currently available.

The large-scale structure growth index $\gamma$ provides a consistency test of the standard cosmology and is a potential indicator of modified gravity. We investigate the constraints on $\gamma$ from next-generation spectroscopic surveys, using the power spectrum that is observed in redshift space, i.e., the angular power spectrum. The angular power spectrum avoids the need for an Alcock-Packzynski correction. It also naturally incorporates cosmic evolution and wide-angle effects, without any approximation. We include the cross-correlations between redshift bins, using a hybrid approximation when the total number of bins is computationally unfeasible. We show that the signal-to-noise on $\gamma$ increases as the redshift bin-width is decreased. Noise per bin also increases -- but this is compensated by the increased number of auto- and cross-spectra. In our forecasts, we marginalise over the amplitude of primordial fluctuations and other standard cosmological parameters, including the dark energy equation of state parameter, as well as the clustering bias. Neglecting cross-bin correlations increases the errors by $\sim40 - 150\%$. Using only linear scales, we find that a DESI-like BGS survey and an HI intensity mapping survey with the SKA1 precursor MeerKAT deliver similar errors of $\sim4-6\%$, while a Euclid-like survey and an SKA1 intensity mapping survey give $\sim3\%$ errors.

Intensity mapping has been attracting increasing interest as a way to study galaxy evolution and the large scale structure of the Universe. Instead of detecting individual galaxies, we measure the integrated emission from a volume of the universe. Contrary to galaxy surveys, it includes light from the faintest objects in the cosmic web. We will start this overview by introducing line intensity mapping and review its current observational status. We will enumerate the most prominent emission lines for intensity mapping and discuss how to model their observed signal. The prospects of using this technique in the near future are very good with a wide range of experiments in different parts of the electromagnetic spectrum. As for any other observational tool, intensity mapping is not free of systematic uncertainties, some of which are intrinsic to it. We will briefly discuss them as well as the means to deal with foreground and background contamination. Here, we will focus on the science cases for intensity mapping in the post-reionization universe. These range from the history of star formation to measuring the acoustic scale and the fundamental physics.

The Square Kilometre Array (SKA) is a planned large radio interferometer designed to operate over a wide range of frequencies, and with an order of magnitude greater sensitivity and survey speed than any current radio telescope. The SKA will address many important topics in astronomy, ranging from planet formation to distant galaxies. However, in this work, we consider the perspective of the SKA as a facility for studying physics. We review four areas in which the SKA is expected to make major contributions to our understanding of fundamental physics: cosmic dawn and reionisation; gravity and gravitational radiation; cosmology and dark energy; and dark matter and astroparticle physics. These discussions demonstrate that the SKA will be a spectacular physics machine, which will provide many new breakthroughs and novel insights on matter, energy and spacetime.

The angular power spectrum is a gauge-independent observable that is in principle the natural tool for analysing galaxy number counts. In practice, the problem is that the computational requirements for next-generation spectroscopic surveys such as Euclid and the Square Kilometre Array are currently unfeasible. We propose a new method to save computational time for spectroscopic angular power spectra. This hybrid method is modelled on the Fourier power spectrum approach of treating relatively thick redshift bins (redshift width ~0.1) as separate surveys. In the hybrid method, each thick bin is further subdivided into thin bins (redshift width ~0.01); all the correlations within each thick bin are computed, while cross-bin correlations beyond the thick bins are neglected. Constraints on cosmological parameters from the hybrid method are comparable to those from the standard galaxy power spectrum analysis - but they have the advantage that cosmic evolution, wide-angle and lensing effects are naturally included, while no Alcock-Paczynski correction is needed. The hybrid method delivers much tighter constraints than a 2D tomographic approach that is typical for photometric surveys, which considers only thick bins and the correlations between them. Furthermore, for standard cosmological parameters our method is not biased by neglecting the effects of lensing on number counts, while the tomographic method is strongly biased.

We present a detailed overview of the cosmological surveys that will be carried out with Phase 1 of the Square Kilometre Array (SKA1), and the science that they will enable. We highlight three main surveys: a medium-deep continuum weak lensing and low-redshift spectroscopic HI galaxy survey over 5,000 sqdeg; a wide and deep continuum galaxy and HI intensity mapping survey over 20,000 sqdeg from z = 0.35 - 3; and a deep, high-redshift HI intensity mapping survey over 100 sqdeg from z = 3 - 6. Taken together, these surveys will achieve an array of important scientific goals: measuring the equation of state of dark energy out to z ~ 3 with percent-level precision measurements of the cosmic expansion rate; constraining possible deviations from General Relativity on cosmological scales by measuring the growth rate of structure through multiple independent methods; mapping the structure of the Universe on the largest accessible scales, thus constraining fundamental properties such as isotropy, homogeneity, and non-Gaussianity; and measuring the HI density and bias out to z = 6. These surveys will also provide highly complementary clustering and weak lensing measurements that have independent systematic uncertainties to those of optical surveys like LSST and Euclid, leading to a multitude of synergies that can improve constraints significantly beyond what optical or radio surveys can achieve on their own. This document, the 2018 Red Book, provides reference technical specifications, cosmological parameter forecasts, and an overview of relevant systematic effects for the three key surveys, and will be regularly updated by the Cosmology Science Working Group in the run up to start of operations and the Key Science Programme of SKA1.

We use full sky simulations, including the effects of foreground contamination and removal, to explore multi-tracer synergies between a SKA-like 21cm intensity mapping survey and a LSST-like photometric galaxy redshift survey. In particular we study ratios of auto and cross-correlations between the two tracers as estimators of the ratio of their biases, a quantity that should benefit considerably from the cosmic variance cancellation of the multi-tracer approach. We show how well we should be able to measure the bias ratio on very large scales (down to $\ell \sim 3$), which is crucial to measure primordial non-Gaussianity and general relativistic effects on large scale structure. We find that, in the absence of foregrounds but with realistic noise levels of such surveys, the multi-tracer estimators are able to improve on the sensitivity of a cosmic-variance contaminated measurement by a factor of $2-4$. When foregrounds are included, estimators using the 21cm auto-correlation become biased. However, we show that cross-correlation estimators are immune to this and do not incur in any significant penalty in terms of sensitivity from discarding the auto-correlation data. However, the loss of long-wavelength radial modes caused by foreground removal in combination with the low redshift resolution of photometric surveys, reduces the sensitivity of the multi-tracer estimator, albeit still better than the cosmic variance contaminated scenario even in the noise free case. Finally we explore different alternative avenues to avoid this problem.

We study synergies between HI 21cm and H$\alpha$ intensity map observations, focusing on SKA1-like and SPHEREx-like surveys. We forecast how well such a combination can measure features in the angular power spectrum on the largest scales, that arise from primordial non-Gaussianity and from general relativistic effects. For the first time we consider Doppler, Sachs-Wolfe and integrated SW effects separately. We confirm that the single-tracer surveys on their own cannot detect general relativistic effects and can constrain the non-Gaussianity parameter $f_{\rm NL}$ only slightly better than Planck. Using the multi-tracer technique, constraints on $f_{\rm NL}$ can be pushed down to $\sim1$. Amongst the general relativistic effects, the Doppler term is detectable with the multi-tracer. The Sachs-Wolfe terms and the integrated SW effect are still not detectable.

Supernova cosmology without spectra will be an important component of future surveys such as LSST. This lack of supernova spectra results in uncertainty in the redshifts which, if ignored, leads to significantly biased estimates of cosmological parameters. Here we present a hierarchical Bayesian formalism -- zBEAMS -- that addresses this problem by marginalising over the unknown or uncertain supernova redshifts to produce unbiased cosmological estimates that are competitive with supernova data with spectroscopically confirmed redshifts. zBEAMS provides a unified treatment of both photometric redshifts and host galaxy misidentification (occurring due to chance galaxy alignments or faint hosts), effectively correcting the inevitable contamination in the Hubble diagram. Like its predecessor BEAMS, our formalism also takes care of non-Ia supernova contamination by marginalising over the unknown supernova type. We illustrate this technique with simulations of supernovae with photometric redshifts and host galaxy misidentification. A novel feature of the photometric redshift case is the important role played by the redshift distribution of the supernovae.

Following the first two annual intensity mapping workshops at Stanford in March 2016 and Johns Hopkins in June 2017, we report on the recent advances in theory, instrumentation and observation that were presented in these meetings and some of the opportunities and challenges that were identified looking forward. With preliminary detections of CO, [CII], Lya and low-redshift 21cm, and a host of experiments set to go online in the next few years, the field is rapidly progressing on all fronts, with great anticipation for a flood of new exciting results. This current snapshot provides an efficient reference for experts in related fields and a useful resource for nonspecialists. We begin by introducing the concept of line-intensity mapping and then discuss the broad array of science goals that will be enabled, ranging from the history of star formation, reionization and galaxy evolution to measuring baryon acoustic oscillations at high redshift and constraining theories of dark matter, modified gravity and dark energy. After reviewing the first detections reported to date, we survey the experimental landscape, presenting the parameters and capabilities of relevant instruments such as COMAP, mmIMe, AIM-CO, CCAT-p, TIME, CONCERTO, CHIME, HIRAX, HERA, STARFIRE, MeerKAT/SKA and SPHEREx. Finally, we describe recent theoretical advances: different approaches to modeling line luminosity functions, several techniques to separate the desired signal from foregrounds, statistical methods to analyze the data, and frameworks to generate realistic intensity map simulations.

We discuss the ground-breaking science that will be possible with a wide area survey, using the MeerKAT telescope, known as MeerKLASS (MeerKAT Large Area Synoptic Survey). The current specifications of MeerKAT make it a great fit for science applications that require large survey speeds but not necessarily high angular resolutions. In particular, for cosmology, a large survey over $\sim 4,000 \, {\rm deg}^2$ for $\sim 4,000$ hours will potentially provide the first ever measurements of the baryon acoustic oscillations using the 21cm intensity mapping technique, with enough accuracy to impose constraints on the nature of dark energy. The combination with multi-wavelength data will give unique additional information, such as exquisite constraints on primordial non-Gaussianity using the multi-tracer technique, as well as a better handle on foregrounds and systematics. Such a wide survey with MeerKAT is also a great match for HI galaxy studies, providing unrivalled statistics in the pre-SKA era for galaxies resolved in the HI emission line beyond local structures at z > 0.01. It will also produce a large continuum galaxy sample down to a depth of about 5\,$\mu$Jy in L-band, which is quite unique over such large areas and will allow studies of the large-scale structure of the Universe out to high redshifts, complementing the galaxy HI survey to form a transformational multi-wavelength approach to study galaxy dynamics and evolution. Finally, the same survey will supply unique information for a range of other science applications, including a large statistical investigation of galaxy clusters as well as produce a rotation measure map across a huge swathe of the sky. The MeerKLASS survey will be a crucial step on the road to using SKA1-MID for cosmological applications and other commensal surveys, as described in the top priority SKA key science projects (abridged).

It is usually assumed that we will need to wait until next-generation surveys like Euclid, LSST and SKA, in order to improve on the current best constraints on primordial non-Gaussianity from the Planck experiment. We show that two contemporary surveys, with the SKA precursor MeerKAT and the Dark Energy Survey (DES), can be combined using the multi-tracer technique to deliver an accuracy on measurement of $f_{\rm NL}$ that is up to three times better than Planck.

We present a systematic study of the intensity mapping technique using updated models for the different emission lines from galaxies and identify which ones are more promising for cosmological studies of the post reionization epoch. We consider the emission of ${\rm Ly\alpha}$, ${\rm H\alpha}$, H$\beta$, optical and infrared oxygen lines, nitrogen lines, CII and the CO rotational lines. We then identify that ${\rm Ly\alpha}$, ${\rm H\alpha}$, OII, CII and the lowest rotational CO lines are the best candidates to be used as IM probes. These lines form a complementary set of probes of the galaxies emission spectra. We then use reasonable experimental setups from current, planned or proposed experiments to access the detectability of the power spectrum of each emission line. Intensity mapping of ${\rm Ly\alpha}$ emission from $z=2$ to 3 will be possible in the near future with HETDEX, while far-infrared lines require new dedicated experiments. We also show that the proposed SPHEREx satellite can use OII and ${\rm H\alpha}$ IM to study the large-scale distribution of matter in intermediate redshifts of 1 to 4. We found that submilimeter experiments with bolometers can have similar performances at intermediate redshifts using CII and CO(3-2).

Next-generation cosmological surveys will probe ever larger volumes of the Universe, including the largest scales, near and beyond the horizon. On these scales, the galaxy power spectrum carries signatures of local primordial non-Gaussianity (PNG) and horizon-scale general relativistic (GR) effects. However, cosmic variance limits the detection of horizon-scale effects. Combining different surveys via the multi-tracer method allows us to reduce the effect down cosmic variance. This method benefits from large bias differences between two tracers of the underlying dark matter distribution, which suggests a multi-wavelength combination of large volume surveys that are planned on a similar timescale. We show that the combination of two contemporaneous surveys, a large neutral hydrogen intensity mapping survey in SKA Phase1 and a Euclid-like photometric survey, will provide unprecedented constraints on PNG as well as detection of the GR effects. We forecast that the error on local PNG will break through the cosmic variance limit on cosmic microwave background surveys and achieve $\sigma(f_{NL})\simeq1.4-0.5$, depending on assumed priors, bias, and sky coverage. GR effects are more robust to changes in the assumed fiducial model, and we forecast that they can be detected with a signal-to-noise of about $14$.

A new dynamic control algorithm in order to direct the trajectory of a glider to a pre-assigned target point is proposed. The algorithms runs iteratively and the approach to the target point is self-correcting. The algorithm is applicable to any non-powered lift-enabled vehicle (glider) travelling in planetary atmospheres. As a proof of concept, we have applied the new algorithm to the command and control of the trajectory of the Space Shuttle during the Terminal Area Energy Management (TAEM) phase.

We characterise the primordial perturbations produced due to both inflaton and curvaton fluctuations in models where the curvaton has a quadratic, cosine or hyperbolic potential, and the inflaton potential is characterised by the usual slow-roll parameters. Isocurvature curvaton field perturbations can produce significant non-Gaussianity in the primordial density field, in contrast with adiabatic inflaton field perturbations which produce negligible non-Gaussianity for canonical scalar fields. A non-self-interacting curvaton with quadratic potential produces a local-type non-Gaussianity that is well described by the non-linearity parameter fNL, which may be scale-dependent when the inflaton perturbations dominate the power spectrum. We show how observational bounds on non-linearity parameters and the tensor-scalar ratio can be used to constrain curvaton and inflaton parameters. We find a consistency relation between the bispectrum and trispectrum parameters in a mixed inflaton-curvaton model for a quadratic curvaton potential. Self-interaction terms in the curvaton potential can lead to both a large trispectrum parameter, gNL, and scale-dependence of the non-linearity parameters.

We consider a simple model of cosmological collapse driven by canonical fields with exponential potentials. We generalise the two-field ekpyrotic collapse to consider non-orthogonal or tilted potentials and give the general condition for isocurvature field fluctuations to have a scale-invariant spectrum in this model. In particular we show that tilted potentials allow for a slightly red spectrum of perturbations as required by current observations. However a red spectrum of fluctuations implies that the two-field ekpyrotic phase must have a finite duration and requires a preceding phase which sets the initial conditions for what otherwise appears to be a fine-tuned trajectory in the phase space.

In this paper we consider how non-Gaussianity of the primordial density perturbation and the amplitude of gravitational waves from inflation can be used to determine parameters of the curvaton scenario for the origin of structure. We show that in the simplest quadratic model, where the curvaton evolves as a free scalar field, measurement of the bispectrum relative to the power spectrum, fNL, and the tensor-to-scalar ratio can determine both the expectation value of the curvaton field during inflation and its dimensionless decay rate relative to the curvaton mass. We show how these predictions are altered by the introduction of self-interactions, in models where higher-order corrections are determined by a characteristic mass scale and discuss how additional information about primordial non-Gaussianity and scale dependence may constrain curvaton interactions.

We estimate large-scale curvature perturbations from isocurvature fluctuations in the waterfall field during hybrid inflation, in addition to the usual inflaton field perturbations. The tachyonic instability at the end of inflation leads to an explosive growth of super-Hubble scale perturbations, but they retain the steep blue spectrum characteristic of vacuum fluctuations in a massive field during inflation. The power spectrum thus peaks around the Hubble-horizon scale at the end of inflation. We extend the usual delta-N formalism to include the essential role of these small fluctuations when estimating the large-scale curvature perturbation. The resulting curvature perturbation due to fluctuations in the waterfall field is second-order and the spectrum is expected to be of order 10^{-54} on cosmological scales.