The cosmological coincidences between the matter and radiation energy densities at recombination as well as between the densities of matter and the cosmological constant at present time are well known. We point out that moreover the third intersection between the energy densities of radiation and the cosmological constant coincides with the reionization epoch. To quantify the statistical relevance of this concurrence, we compute the Bayes factor between the concordance cosmology with free Thomson scattering optical depth and a model for which this parameter is inferred from imposing a match between the time of density equality and the epoch of reionization. This is to characterize the potential explanatory gain if one were to find a parameter-free physical connection. We find a very strong preference for such a concurrence on the Jeffreys scale from current cosmological observations. We furthermore discuss the effect of choice of priors, changes in reionization history, and free sum of neutrino masses. We also estimate the impact of adding intermediate polarization data from the Planck High Frequency Instrument and prospects for future 21 cm surveys. In the first case, preference for the correlation remains substantial, whereas future data may give results more decisive in pro or substantial in contra. Finally, we provide a discussion on different interpretations of these findings. In particular, we show how a connection between the star-formation history and the cosmological background dynamics can give rise to this concurrence.

The Pseudo-Nambu-Goldstone Boson (PNGB) potential, defined through the amplitude $M^4$ and width $f$ of its characteristic potential $V(\phi) = M^4[1 + \cos(\phi/f)]$, is one of the best-suited models for the study of thawing quintessence. We analyze its present observational constraints by direct numerical solution of the scalar field equation of motion. Observational bounds are obtained using Supernovae data, cosmic microwave background temperature, polarization and lensing data from {\it Planck}, direct Hubble constant constraints, and baryon acoustic oscillations data. We find the parameter ranges for which PNGB quintessence gives a viable theory for dark energy. This exact approach is contrasted with the use of an approximate equation-of-state parametrization for thawing theories. We also discuss other possible parameterization choices, as well as commenting on the accuracy of the constraints imposed by {\it Planck} alone. Overall our analysis highlights a significant prior dependence to the outcome coming from the choice of modelling methodology, which current data are not sufficient to override.

Most of the information on our cosmos stems from either late-time observations or the imprint of early-time inhomogeneities on the cosmic microwave background. We explore to what extent early modifications of gravity, which become significant after recombination but then decay towards the present, can be constrained by current cosmological observations. For the evolution of the gravitational modification, we adopt the decaying mode of a hybrid-metric Palatini $f(\mathcal{R})$ gravity model which is designed to reproduce the standard cosmological background expansion history and due to the decay of the modification is naturally compatible with Solar-System tests. We embed the model in the effective field theory description of Horndeski scalar-tensor gravity with an early-time decoupling of the gravitational modification. Since the quasistatic approximation for the perturbations in the model breaks down at high redshifts, where modifications remain relevant, we introduce a computationally efficient correction to describe the evolution of the scalar field fluctuation in this regime. We compare the decaying early-time modification against geometric probes and recent Planck measurements and find no evidence for such effects in the observations. Current data constrains the scalar field value at $|f_{\mathcal{R}}(z=z_{\rm on})| \lesssim 10^{-2}$ for modifications introduced at redshifts $z_{\rm on}\sim(500-1000)$ with present-day value $|f_{\mathcal{R}0}|\lesssim10^{-8}$. Finally, we comment on constraints that will be achievable with future 21~cm surveys and gravitational wave experiments.

In this work, we introduce two models of the hybrid metric-Palatini theory of gravitation. We explore their background evolution, showing explicitly that one recovers standard General Relativity with an effective Cosmological Constant at late times. This happens because the Palatini Ricci scalar evolves towards and asymptotically settles at the minimum of its effective potential during cosmological evolution. We then use a combination of cosmic microwave background, supernovae and baryonic accoustic oscillations background data to constrain the models' free parameters. For both models, we are able to constrain the maximum deviation from the gravitational constant $G$ one can have at early times to be around $1\%$.