Axion-like particles (ALPs) decaying into photons are known to affect a wide range of astrophysical and cosmological observables. In this study we focus on ALPs with masses in the keV–MeV range and lifetimes between 10⁴ and 10¹³ seconds, corresponding to decays between the end of Big Bang Nucleosynthesis and the formation of the Cosmic Microwave Background (CMB). Using the CosmoBit module of the global fitting framework GAMBIT, we combine state-of-the-art calculations of the irreducible ALP freeze-in abundance, primordial element abundances (including photodisintegration through ALP decays), CMB spectral distortions and anisotropies, and constraints from supernovae and stellar cooling. This approach makes it possible for the first time to perform a global analysis of the ALP parameter space while varying the parameters of ΛCDM as well as several nuisance parameters. We find a lower bound on the ALP mass of around m_a > 300 keV, which can only be evaded if ALPs are stable on cosmological timescales. Future observations of CMB spectral distortions with a PIXIE-like mission are expected to improve this bound by two orders of magnitude.

We introduce CosmoBit, a module within the open-source GAMBIT software framework for exploring connections between cosmology and particle physics with joint global fits. CosmoBit provides a flexible framework for studying various scenarios beyond ACDM, such as models of inflation, modifications of the effective number of relativistic degrees of freedom, exotic energy injection from annihilating or decaying dark matter, and variations of the properties of elementary particles such as neutrino masses and the lifetime of the neutron. Many observables and likelihoods in CosmoBit are computed via interfaces to AlterBBN, CLASS, DarkAges, MontePython, MultiModeCode, and plc. This makes it possible to apply a wide range of constraints from large-scale structure, Type Ia supernovae, Big Bang Nucleosynthesis and the cosmic microwave background. Parameter scans can be performed using the many different statistical sampling algorithms available within the GAMBIT framework, and results can be combined with calculations from other GAMBIT modules focused on particle physics and dark matter. We include extensive validation plots and a first application to scenarios with non-standard relativistic degrees of freedom and neutrino temperature, showing that the corresponding constraint on the sum of neutrino masses is much weaker than in the standard scenario.

The dark sectors of our Universe, dark matter and dark energy, together constitute about 96 % of the total energy content of the Universe. To date, we only have observational evidence for their existence. What is still lacking is a complete theoretical framework consistent with all observational data to embed a dark matter particle or component into the standard models of particle physics and cosmology, as well as an explanation for the nature or origin of dark energy.

Since the discovery of these dark components decades ago, a variety of different theories have been proposed to overcome the shortcomings of our current standard models. To assess the viability of these non-standard theories, they ideally should be tested against all relevant available datasets. In this thesis, I show two examples of how cosmological and astrophysical observables are used to constrain or even rule out non-standard cosmological models. Further, I present the first software tool that provides a general framework to test non-standard physics with global fits to data from particle physics and cosmology simultaneously.

The first example is minimally coupled covariant Galileons, a modification of General Relativity to explain dark energy without the need for a fine-tuned cosmological constant. I demonstrate how the combination of constraints arising from the integrated Sachs-Wolf effect and the propagation speed of gravitational waves can rule out all three branches of the theory.

The second example shows how the existence and parameter space of cosmic superstrings can be constrained. These are the hypothesised fundamental building blocks of Type IIb Superstring theory, stretched out to cosmological scales during the phase of inflation. The theory can be tested through the unique microlensing signature of cosmic superstrings when crossing the line of sight of an observer monitoring a point-like source. I show how, based on simulations, we can estimate the expected detection rates from observations of distant Type Ia Supernovae and stars in Andromeda; from these estimates I assess the implications for the theory.

Finally, I present CosmoBit, a new module for the Global and Modular Beyond-Standard Model Inference Tool (GAMBIT). \gambit allows the user to test a variety of extensions to the Standard Model of particle physics against data from, e.g. collider searches, dark matter direct and indirect detection experiments, as well as laboratory measurements of neutrino properties. CosmoBit augments this with the inclusion of cosmological likelihoods. This addition opens up the possibility to test a given model against data from, e.g. the Big Bang Nucleosynthesis proceeding minutes after the Big Bang, probes of the Cosmic Microwave Background ~ 380,000 years later, and (laboratory) measurements from the present day, 13.8 billion years after the Big Bang. Including measurements that span several different epochs and orders of magnitude in energy, the combination of CosmoBit with other GAMBIT modules provides a promising tool for shedding light on the dark sectors of the Universe.

The existence of cosmic superstrings may be probed by astronomical time domain surveys. When crossing the line of sight to point-like sources, strings produce a distinctivemicrolensing signature. We consider two avenues to hunt for a relic population of superstring loops: frequent monitoring of (1) stars in Andromeda, lensed by loops in the haloes of the Milky-Way and Andromeda and (2) supernovae at cosmological distances, lensed by loops in the intergalactic medium. We assess the potential of such experiments to detect and/or constrain strings with a range of tensions, 10(-15) less than or similar to G mu/c(2) less than or similar to 10(-6). The practical sensitivity is tied to cadence of observations which we explore in detail. We forecast that high-cadence monitoring of similar to 10(5) stars on the far side of Andromeda over a year-long period will detect microlensing events if G mu/c(2) similar to 10(-13), while similar to 10(6) stars will detect events if 10(-13.5) < G mu/c(2) < 10(-11.5); the upper and lower bounds of the accessible tension range continue to expand as the number of stars rises. We also analyse the ability to reject models in the absence of fluctuations. While challenging, these studies are within reach of forthcoming time-domain surveys. Supernova observations can hypothetically constrain models with 10(-12) < G mu/c(2) < 10(-6) without any optimization of the survey cadence. However, the event rate forecast suggests it will be difficult to reject models of interest. As a demonstration, we use observations from the Pantheon Type Ia supernova cosmology data set to place modest constraints on the number density of cosmic superstrings in a poorly tested region of the parameter space.

Cosmological models with Galileon gravity are an alternative to the standard ACDM paradigm with testable predictions at the level of its self-accelerating solutions for the expansion history, as well as large-scale structure formation. Here, we place constraints on the full parameter space of these models using data from the cosmic microwave background (CMB) (including lensing), baryonic acoustic oscillations (BAO) and the Integrated Sachs Wolfe (ISW) effect. We pay special attention to the ISW effect for which we use the cross spectra, C-l(Tg), of CMB temperature maps and foreground galaxies from the WISE survey. The sign of C-l(Tg) is set by the time evolution of the lensing potential in the redshift range of the galaxy sample: it is positive if the potential decays (like in ACDM), negative if it deepens. We constrain three subsets of Galileon gravity separately known as the Cubic, Quartic and Quintic Galileons. The cubic Galileon model predicts a negative C-l(Tg) and exhibits a 7.8 sigma tension with the data, which effectively rules it out. For the quartic and quintic models the ISW data also rule out a significant portion of the parameter space but permit regions where the goodness-of-fit is comparable to ACDM. The data prefers a non zero sum of the neutrino masses (Sigma m(v) approximate to 0.5eV) with similar to 5 sigma significance in these models. The best-fitting models have values of Ho consistent with local determinations, thereby avoiding the tension that exists in ACDM. We also identify and discuss a similar to 2 sigma tension that Galileon gravity exhibits with recent BAO measurements. Our analysis shows overall that Galileon cosmologies cannot be ruled out by current data but future lensing, BAO and ISW data hold strong potential to do so.

We address the impact of consistent modifications of gravity on the largest observable scales, focusing on relativistic effects in galaxy number counts and the cross-correlation between the matter large scale structure (LSS) distribution and the cosmic microwave background (CMB). Our analysis applies to a very broad class of general scalar-tensor theories encoded in the Horndeski Lagrangian and is fully consistent on linear scales, retaining the full dynamics of the scalar field and not assuming quasi-static evolution. As particular examples we consider self-accelerating Covariant Galileons, Brans-Dicke theory and parameterizations based on the effective field theory of dark energy, using the hi_class code to address the impact of these models on relativistic corrections to LSS observables. We find that especially effects which involve integrals along the line of sight (lensing convergence, time delay and the integrated Sachs-Wolfe effect- ISW) can be considerably modified, and even lead to O(1000%) deviations from General Relativity in the case of the ISW effect for Galileon models, for which standard probes such as the growth function only vary by O(10%). These effects become dominant when correlating galaxy number counts at different redshifts and can lead to similar to 50% deviations in the total signal that might be observable by future LSS surveys. Because of their integrated nature, these deep-redshift cross-correlations are sensitive to modifications of gravity even when probing eras much before dark energy domination. We further isolate the ISW effect using the cross-correlation between LSS and CMB temperature anisotropies and use current data to further constrain Horndeski models. Forthcoming large-volume galaxy surveys using multiple-tracers will search for all these effects, opening a new window to probe gravity and cosmic acceleration at the largest scales available in our universe.