The Core Cosmology Library (CCL) provides routines to compute basic cosmological observables to a high degree of accuracy, which have been verified with an extensive suite of validation tests. Predictions are provided for many cosmological quantities, including distances, angular power spectra, correlation functions, halo bias and the halo mass function through state-of-the-art modeling prescriptions available in the literature. Fiducial specifications for the expected galaxy distributions for the Large Synoptic Survey Telescope (LSST) are also included, together with the capability of computing redshift distributions for a user-defined photometric redshift model. A rigorous validation procedure, based on comparisons between CCL and independent software packages, allows us to establish a well-defined numerical accuracy for each predicted quantity. As a result, predictions for correlation functions of galaxy clustering, galaxy-galaxy lensing and cosmic shear are demonstrated to be within a fraction of the expected statistical uncertainty of the observables for the models and in the range of scales of interest to LSST. CCL is an open source software package written in C, with a python interface and publicly available at https://github.com/LSSTDESC/CCL.

In this paper we investigate a time-varying neutrino mass model, motivated by the mild tension between cosmic microwave background (CMB) measurements of the matter fluctuations and those obtained from low-redshift data. We modify the minimal case of the model proposed by Dvali and Funcke (2016) that predicts late neutrino mass generation in a post-recombination cosmic phase transition, by assuming that neutrino asymmetries allow for the presence of relic neutrinos in the late-time Universe. We show that, if the transition is supercooled, current cosmological data (including CMB temperature, polarization and lensing, baryon acoustic oscillations, and Type Ia supernovae) prefer the scale factor $a_s$ of the phase transition to be very large, peaking at $a_s\sim 1$, and therefore supporting a cosmological scenario in which neutrinos are almost massless until very recent times. We find that in this scenario the cosmological bound on the total sum of the neutrino masses today is significantly weakened compared to the standard case of constant-mass neutrinos, with $\sum m_\nu<4.8$~eV at 95\% confidence, and in agreement with the model predictions. The main reason for this weaker bound is a large correlation arising between the dark energy and neutrino components in the presence of false vacuum energy that converts into the non-zero neutrino masses after the transition. This result provides new targets for the coming KATRIN and PTOLEMY experiments. We also show that the time-varying neutrino mass model considered here does not provide a clear explanation to the existing cosmological $\Omega_m$-$\sigma_8$ discrepancies.

Future surveys will access large volumes of space and hence very long wavelength fluctuations of the matter density and gravitational field. It has been argued that the set of secondary effects that affect the galaxy distribution, relativistic in nature, will bring new, complementary cosmological constraints. We study this claim in detail by focusing on a subset of wide-area future surveys: Stage-4 cosmic microwave background experiments and photometric redshift surveys. In particular, we look at the magnification lensing contribution to galaxy clustering and general relativistic corrections to all observables. We quantify the amount of information encoded in these effects in terms of the tightening of the final cosmological constraints as well as the potential bias in inferred parameters associated with neglecting them. We do so for a wide range of cosmological parameters, covering neutrino masses, standard dark-energy parametrizations and scalar-tensor gravity theories. Our results show that, while the effect of lensing magnification to number counts does not contain a significant amount of information when galaxy clustering is combined with cosmic shear measurements, this contribution does play a significant role in biasing estimates on a host of parameter families if unaccounted for. Since the amplitude of the magnification term is controlled by the slope of the source number counts with apparent magnitude, $s(z)$, we also estimate the accuracy to which this quantity must be known to avoid systematic parameter biases, finding that future surveys will need to determine $s(z)$ to the $\sim$5-10\% level. On the contrary, large-scale general-relativistic corrections are irrelevant both in terms of information content and parameter bias for most cosmological parameters, but significant for the level of primordial non-Gaussianity.

We study the correlations between parameters characterizing neutrino physics and the evolution of dark energy. Using a fluid approach, we show that time-varying dark energy models exhibit degeneracies with the cosmic neutrino background over extended periods of the cosmic history, leading to a degraded estimation of the total mass and number of species of neutrinos. We investigate how to break degeneracies and combine multiple probes across cosmic time to anchor the behaviour of the two components. We use Planck CMB data and BAO measurements from the BOSS, SDSS and 6dF surveys to present current limits on the model parameters, and then forecast the future reach from the CMB Stage-4 and DESI experiments. We show that a multi-probe analysis of current data provides only marginal improvement on the determination of the individual parameters and no reduction of the correlations. Future observations will better distinguish the neutrino mass and preserve the current sensitivity to the number of species even in case of a time-varying dark energy component.