The statistical power of Lyman-${\alpha}$ forest Baryon Acoustic Oscillation (BAO) measurements is set to increase significantly in the coming years as new instruments such as the Dark Energy Spectroscopic Instrument deliver progressively more constraining data. Generating mock datasets for such measurements will be important for validating analysis pipelines and evaluating the effects of systematics. With such studies in mind, we present LyaCoLoRe: a package for producing synthetic Lyman-${\alpha}$ forest survey datasets for BAO analyses. LyaCoLoRe transforms initial Gaussian random field skewers into skewers of transmitted flux fraction via a number of fast approximations. In this work we explain the methods of producing mock datasets used in LyaCoLoRe, and then measure correlation functions on a suite of realisations of such data. We demonstrate that we are able to recover the correct BAO signal, as well as large-scale bias parameters similar to literature values. Finally, we briefly describe methods to add further astrophysical effects to our skewers - high column density systems and metal absorbers - which act as potential complications for BAO analyses.

We utilise mock catalogues from high-accuracy cosmological $N$-body simulations to quantify shifts in the recovery of the acoustic scale that could potentially result from galaxy clustering bias. The relationship between galaxies and dark matter halos presents a complicated source of systematic errors in modern redshift surveys, particularly when aiming to make cosmological measurements to sub-percent precision. Apart from a scalar, linear bias parameter accounting for the density contrast ratio between matter tracers and the true matter distribution, other types of galaxy bias, such as assembly and velocity biases, may also significantly alter clustering signals from small to large scales. We create mocks based on generalised halo occupation populations of 36 periodic boxes from the \abacuscosmos release with. In a total volume of $48 \, h^{-3} \mathrm{Gpc}^3$, we test various biased models along with an unbiased base case. Two reconstruction methods are applied to galaxy samples and the apparent acoustic scale is derived by fitting the two-point correlation function multipoles. With respect to the baseline, we find a $0.3\%$ shift in the line-of-sight acoustic scale for one variation in the satellite galaxy population, and we find an $0.7\%$ shift for an extreme level of velocity bias of the central galaxies. All other bias models are consistent with zero shift at the $0.2\%$ level after reconstruction. We note that the bias models explored are relatively large variations, producing sizeable and likely distinguishable changes in small-scale clustering, the modelling of which would further calibrate the BAO standard ruler.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later.