In this work, we obtain measurements of the Hubble constant in the context of modified gravity theories. We set up our theoretical framework by considering viable cosmological $f(R)$ and $f(T)$ models, and we analyzed them through the use of geometrical data sets obtained in a model-independent way, namely, gravitationally lensed quasars with measured time delays, standard clocks from cosmic chronometers, and standard candles from the Pantheon Supernovae Ia sample. We find $H_0=(72.4\pm 1.4)$ km s$^{-1}$ Mpc$^{-1}$ and $H_0=(71.5\pm 1.3)$ km s$^{-1}$ Mpc$^{-1}$ for the $f(R)$ and $f(T)$ models, respectively. Our results represent 1.9\% and 1.8\% measurements of the Hubble constant, which are fully consistent with the local estimate of $H_0$ by the Hubble Space Telescope. We do not find significant departures from general relativity, as our study shows that the characteristic parameters of the extensions of gravity beyond general relativity are compatible with the $\Lambda$CDM cosmology. Moreover, within the standard cosmological framework, our full joint analysis suggests that it is possible to measure the dark energy equation of state parameter at 1.2\% accuracy, although we find no statistical evidence for deviations from the cosmological constant case.

Cosmography becomes non-predictive when cosmic data span beyond the red shift limit $z\simeq1 $. This leads to a \emph{strong convergence issue} that jeopardizes its viability. In this work, we critically compare the two main solutions of the convergence problem, i.e. the $y$-parametrizations of the redshift and the alternatives to Taylor expansions based on Pad\'e series. In particular, among several possibilities, we consider two widely adopted parametrizations, namely $y_1=1-a$ and $y_2=\arctan(a^{-1}-1)$, being $a$ the scale factor of the Universe. We find that the $y_2$-parametrization performs relatively better than the $y_1$-parametrization over the whole redshift domain. Even though $y_2$ overcomes the issues of $y_1$, we get that the most viable approximations of the luminosity distance $d_L(z)$ are given in terms of Pad\'e approximations. In order to check this result by means of cosmic data, we analyze the Pad\'e approximations up to the fifth order, and compare these series with the corresponding $y$-variables of the same orders. We investigate two distinct domains involving Monte Carlo analysis on the Pantheon Superovae Ia data, $H(z)$ and shift parameter measurements. We conclude that the (2,1) Pad\'e approximation is statistically the optimal approach to explain low and high-redshift data, together with the fifth-order $y_2$-parametrization. At high redshifts, the (3,2) Pad\'e approximation cannot be fully excluded, while the (2,2) Pad\'e one is essentially ruled out.

The observation of GW170817 binary neutron star (BNS) merger event has imposed strong bounds on the speed of gravitational waves (GWs) locally, inferring that the speed of GWs propagation is equal to the speed of light. Current GW detectors in operation will not be able to observe BNS merger to long cosmological distance, where possible cosmological corrections on the cosmic expansion history are expected to play an important role, specially for investigating possible deviations from general relativity. Future GW detectors designer projects will be able to detect many coalescences of BNS at high $z$, such as the third generation of the ground GW detector called Einstein Telescope (ET) and the space-based detector deci-hertz interferometer gravitational wave observatory (DECIGO). In this paper, we relax the condition $c_T/c = 1$ to investigate modified GW propagation where the speed of GWs propagation is not necessarily equal to the speed of light. Also, we consider the possibility for the running of the Planck mass corrections on modified GW propagation. We parametrize both corrections in terms of an effective GW luminosity distance and we perform a forecast analysis using standard siren events from BNS mergers, within the sensitivity predicted for the ET and DECIGO. We find at high $z$ very strong forecast bounds on the running of the Planck mass, namely $\mathcal{O}(10^{-1})$ and $\mathcal{O}(10^{-2})$ from ET and DECIGO, respectively. Possible anomalies on GW propagation are bound to $|c_T/c - 1| \leq 10^{-2} \,\,\, (10^{-2})$ from ET (DECIGO), respectively. We finally discuss the consequences of our results on modified gravity phenomenology.

Adopting the Tolman-Oppenheimer-Volkoff formalism, we propose a new analytical solution for a static and spherically symmetric black hole where the cosmological constant is generalized to a non-isotropic Anton-Schmidt fluid acting as a unified dark energy - dark matter source. Our novel result can describe a black hole in energetic equilibrium with the surrounding universe. We thus investigate the interplay between the mass of the black hole and the parameters of the Anton-Schmidt equation of state, and physically interpret the solution as a family of space-time metrics that may describe both Schwarzschild-de Sitter black holes and naked singularities. The result opens a new window for further constraining the physical properties of the cosmic fluid using the Event Horizon Telescope results, which complement the analysis relying on classical cosmological observations.

Standard sirens are the gravitational wave (GW) analog of the astronomical standard candles, and can provide powerful information about the dynamics of the Universe. In this work, we simulate a catalog with 1000 standard siren events from binary neutron star mergers, within the sensitivity predicted for the third generation of the ground GW detector called Einstein telescope. After correctly modifying the propagation of GWs as input to generate the catalog, we apply our mock data set on scalar-tensor theories where the speed of GW propagation is equal to the speed of light. As a first application, we find new observational bounds on the running of the Planck mass, when considering appropriate values within the stability condition of the theory, and we discuss some consequences on the amplitude of the running of the Planck mass. In the second part, we combine our simulated standard sirens catalog with other geometric cosmological tests (Supernovae Ia and cosmic chronometers measurements) to constrain the Hu-Sawicki $f(R)$ gravity model. We thus find new and non-null deviations from the standard $\Lambda$CDM model, showing that in the future the $f(R)$ gravity can be tested up to 95\% confidence level. The results obtained here show that the statistical accuracy achievable by future ground based GW observations, mainly with the ET detector (and planed detectors with a similar sensitivity), can provide strong observational bounds on modified gravity theories.

We here propose a new model-independent technique to overcome the circularity problem affecting the use of Gamma-Ray Bursts (GRBs) as distance indicators through the use of $E_{\rm p}$--$E_{\rm iso}$ correlation. We calibrate the $E_{\rm p}$--$E_{\rm iso}$ correlation and find the GRB distance moduli that can be used to constrain dark energy models. We use observational Hubble data to approximate the cosmic evolution through B\'ezier parametric curve obtained through the linear combination of Bernstein basis polynomials. In so doing, we build up a new data set consisting of 193 GRB distance moduli. We combine this sample with the supernova JLA data set to test the standard $\Lambda$CDM model and its $w$CDM extension. We place observational constraints on the cosmological parameters through Markov Chain Monte Carlo numerical technique. Moreover, we compare the theoretical scenarios by performing the AIC and DIC statistics. For the $\Lambda$CDM model we find $\Omega_m=0.397^{+0.040}_{-0.039}$ at the $2\sigma$ level, while for the $w$CDM model we obtain $\Omega_m=0.34^{+0.13}_{-0.15}$ and $w=-0.86^{+0.36}_{-0.38}$ at the $2\sigma$ level. Our analysis suggests that $\Lambda$CDM model is statistically favoured over the $w$CDM scenario. No evidence for extension of the $\Lambda$CDM model is found.

We apply the holographic principle in the cosmological context through the nonadditive Tsallis entropy, used to describe the thermodynamic properties of nonstandard statistical systems such as the gravitational ones. Assuming the future event horizon as the infrared cutoff, we build a dark energy model free from cosmological inconsistencies, which includes standard thermodynamics and standard holographic dark energy as a limiting case. We thus describe the dynamics of Tsallis holographic dark energy in a flat FLRW background. Hence, we investigate cosmological perturbations in the linear regime on sub-horizon scales. We study the growth of matter fluctuations in the case of clustering dark matter and a homogeneous dark energy component. Furthermore, we employ the most recent late-time cosmic data to test the observational viability of our theoretical scenario. We thus obtain constraints on the free parameters of the model by means of Monte Carlo numerical method. We also used Bayesian selection criteria to estimate the statistical preference for Tsallis holographic dark energy compared to the concordance $\Lambda$CDM paradigm. Our results show deviations from standard holographic dark energy within the $2\sigma$ confidence level. Finally, the analysis of the dark energy equation of state indicates a quintessence-like behaviour with no evidence for phantom-divide crossing at the $1\sigma$ level.

Cosmography can be considered as a sort of a model-independent approach to tackle the dark energy/modified gravity problem. In this review, the success and the shortcomings of the $\Lambda$CDM model, based on General Relativity and standard model of particles, are discussed in view of the most recent observational constraints. The motivations for considering extensions and modifications of General Relativity are taken into account, with particular attention to $f(R)$ and $f(T)$ theories of gravity where dynamics is represented by curvature or torsion field respectively. The features of $f(R)$ models are explored in metric and Palatini formalisms. We discuss the connection between $f(R)$ gravity and scalar-tensor theories highlighting the role of conformal transformations in the Einstein and Jordan frames. Cosmological dynamics of $f(R)$ models is investigated through the corresponding viability criteria. Afterwards, the equivalent formulation of General Relativity (Teleparallel Equivalent General Relativity) in terms of torsion and its extension to $f(T)$ gravity is considered. Finally, the cosmographic method is adopted to break the degeneracy among dark energy models. A novel approach, built upon rational Pad\'e and Chebyshev polynomials, is proposed to overcome limits of standard cosmography based on Taylor expansion. The approach provides accurate model-independent approximations of the Hubble flow. Numerical analyses, based on Monte Carlo Markov Chain integration of cosmic data, are presented to bound coefficients of the cosmographic series. These techniques are thus applied to reconstruct $f(R)$ and $f(T)$ functions and to frame the late-time expansion history of the universe with no \emph{a priori} assumptions on its equation of state. A comparison between the $\Lambda$CDM cosmological model with $f(R)$ and $f(T)$ models is reported.

In this paper we describe QUBIC, an experiment that will observe the polarized microwave sky with a novel approach, which combines the sensitivity of state-of-the art bolometric detectors with the systematic effects control typical of interferometers. QUBIC unique features are the so-called "self-calibration", a technique that allows us to clean the measured data from instrumental effects, and its spectral imaging power, i.e. the ability to separate the signal in various sub-bands within each frequency band. QUBIC will observe the sky in two main frequency bands: 150 GHz and 220 GHz. A technological demonstrator is currently under testing and will be deployed in Argentina during 2019, while the final instrument is expected to be installed during 2020.

The effective theory of the Anton-Schmidt cosmic fluid within the Debye approximation is investigated. In this picture, the universe is modeled out by means of a medium without cosmological constant. In particular, the Anton-Schmidt representation of matter describes the pressure of crystalline solids under deformations imposed by isotropic stresses. The approach scheme is related to the fact that the universe deforms under the action of the cosmic expansion itself. Thus, we frame the dark energy term as a function of scalar fields and obtain the corresponding dark energy potential $V(\varphi)$. Different epochs of the universe evolution are investigated in terms of the evolution of $\varphi$. We show how the Anton-Schmidt equation of state is capable of describing both late and early epochs of cosmic evolution. Finally, numerical bounds on the Anton-Schmidt model with $n=-1$ are derived through a Markov Chain Monte Carlo analysis on the combination of data coming from type Ia Supernovae, observations of Hubble parameter and baryon acoustic oscillations. Statistical comparison with the $\Lambda$CDM model is performed by the AIC and BIC selection criteria. Results are in excellent agreement with the low-redshift data. A further generalization of the model is presented to satisfy the theoretical predictions at early-stage cosmology.

We here study extended classes of logotropic fluids as \textit{unified dark energy models}. Under the hypothesis of the Anton-Schmidt scenario, we consider the universe obeying a single fluid whose pressure evolves through a logarithmic equation of state. This result is in analogy with crystals under isotropic stresses. Thus, we investigate thermodynamic and dynamical consequences by integrating the speed of sound to obtain the pressure in terms of the density, leading to an extended version of the Anton-Schmidt cosmic fluids. Within this picture, we get significant outcomes expanding the Anton-Schmidt pressure in the infrared regime. The low-energy case becomes relevant for the universe to accelerate without any cosmological constant. We therefore derive the effective representation of our fluid in terms of a Lagrangian $\mathcal{L}=\mathcal{L}(X)$, depending on the kinetic term $X$ only. We analyze both the relativistic and non-relativistic limits. In the non-relativistic limit we construct both the Hamiltonian and Lagrangian in terms of density $\rho$ and scalar field $\vartheta$, whereas in the relativistic case no analytical expression for the Lagrangian can be found. Thus, we obtain the potential as a function of $\rho$, under the hypothesis of irrotational perfect fluid. We demonstrate that the model represents a natural generalization of \emph{logotropic dark energy models}. Finally, we analyze an extended class of generalized Chaplygin gas models with one extra parameter $\beta$. Interestingly, we find that the Lagrangians of this scenario and the pure logotropic one coincide in the non-relativistic regime.

We study the growth rate of matter perturbations in the context of teleparallel dark energy in a flat universe. We investigate the dynamics of different theoretical scenarios based on specific forms of the scalar field potential. Allowing for non-minimal coupling between torsion scalar and scalar field, we perform a phase-space analysis of the autonomous systems of equations through the study of critical points. We thus analyze the stability of the critical points, and discuss the cosmological implications searching for possible attractor solutions at late times. Furthermore, combing the growth rate data and the Hubble rate measurements, we place observational constraints on the cosmological parameters of the models through Monte Carlo numerical method. We find that the scenario with a non-minimal coupling is favoured with respect to the standard quintessence case. Adopting the best-fit results, we show that the dark energy equation of state parameter can cross the phantom divide. Finally, we compare our results with the predictions of the concordance $\Lambda$CDM paradigm by performing Bayesian model selection.

A kinematic treatment to trace out the form of $f(R)$ cosmology, within the Palatini formalism, is discussed by only postulating the universe homogeneity and isotropy. To figure this out we build model-independent approximations of the luminosity distance through rational expansions. These approximants extend the Taylor convergence radii computed for usual cosmographic series. We thus consider both Pad\'e and the rational Chebyshev polynomials. They can be used to accurately describe the universe late-time expansion history, providing further information on the thermal properties of all effective cosmic fluids entering the energy momentum tensor of Palatini's gravity. To perform our numerical analysis, we relate the Palatini's Ricci scalar with the Hubble parameter $H$ and thus we write down a single differential equation in terms of the redshift $z$. Therefore, to bound $f(R)$, we make use of the most recent outcomes over the cosmographic parameters obtained from combined data surveys. In particular our clue is to select two scenarios, i.e. $(2,2)$ Pad\'e and $(2,1)$ Chebyshev approximations, since they well approximate the luminosity distance at the lowest possible order. We find that best analytical matches to the numerical solutions lead to $f(R)=a+bR^n$ with free parameters given by the set $(a, b, n)=(-1.627, 0.866, 1.074)$ for $(2,2)$ Pad\'e approximation, whereas $f(R)=\alpha+\beta R^m$ with $(\alpha, \beta, m)=(-1.332, 0.749, 1.124)$ for $(2,1)$ rational Chebyshev approximation. Finally, our results are compared with the $\Lambda$CDM predictions and with previous studies in the literature. Slight departures from General Relativity are also discussed.

QUBIC, the QU Bolometric Interferometer for Cosmology, is a novel forthcoming instrument to measure the B-mode polarization anisotropy of the Cosmic Microwave Background. The detection of the B-mode signal will be extremely challenging; QUBIC has been designed to address this with a novel approach, namely bolometric interferometry. The receiver cryostat is exceptionally large and cools complex optical and detector stages to 40 K, 4 K, 1 K and 350 mK using two pulse tube coolers, a novel 4He sorption cooler and a double-stage 3He/4He sorption cooler. We discuss the thermal and mechanical design of the cryostat, modelling and thermal analysis, and laboratory cryogenic testing.

In the present study, we consider an extended form of teleparallel Lagrangian $f(T,\phi,X)$, as function of a scalar field $\phi$, its kinetic term $X$ and the torsion scalar $T$. We use linear perturbations to obtain the equation of matter density perturbations on sub-Hubble scales. The gravitational coupling is modified in scalar modes with respect to the one of General Relativity, albeit vector modes decay and do not show any significant effects. We thus extend these results by involving multiple scalar field models. Further, we study conformal transformations in teleparallel gravity and we obtain the coupling as the scalar field is non-minimally coupled to both torsion and boundary terms. Finally, we propose the specific model $f(T,\phi,X)=T + \partial_\mu \phi\ \partial^\mu \phi +\xi T \phi^2$. To check its goodness, we employ the observational Hubble data, constraining the coupling constant, $\xi$, through a Monte Carlo technique based on the Metropolis-Hastings algorithm. Hence, fixing $\xi$ to its best-fit value got from our numerical analysis, we calculate the growth rate of matter perturbations and we compare our outcomes with the latest measurements and the predictions of the $\Lambda$CDM model.

In this paper we describe QUBIC, an experiment that takes up the challenge posed by the detection of primordial gravitational waves with a novel approach, that combines the sensitivity of state-of-the art bolometric detectors with the systematic effects control typical of interferometers. The so-called "self-calibration" is a technique deeply rooted in the interferometric nature of the instrument and allows us to clean the measured data from instrumental effects. The first module of QUBIC is a dual band instrument (150 GHz and 220 GHz) that will be deployed in Argentina during the Fall 2018.

We here propose a new class of barotropic factor for matter, motivated by properties of isotropic deformations of crystalline solids. Our approach is dubbed Anton-Schmidt's equation of state and provides a non-vanishing, albeit small, pressure term for matter. The corresponding pressure is thus proportional to the logarithm of universe's volume, i.e. to the density itself since $V\propto \rho^{-1}$. In the context of solid state physics, we demonstrate that by only invoking standard matter with such a property, we are able to frame the universe speed up in a suitable way, without invoking a dark energy term by hand. Our model extends a recent class of dark energy paradigms named \emph{logotropic} dark fluids and depends upon two free parameters, namely $n$ and $B$. Within the Debye approximation, we find that $n$ and $B$ are related to the Gr\"uneisen parameter and the bulk modulus of crystals. We thus show the main differences between our model and the logotropic scenario, and we highlight the most relevant properties of our new equation of state on the background cosmology. Discussions on both kinematics and dynamics of our new model have been presented. We demonstrate that the $\Lambda$CDM model is inside our approach, as limiting case. Comparisons with CPL parametrization have been also reported in the text. Finally, a Monte Carlo analysis on the most recent low-redshift cosmological data allowed us to place constraints on $n$ and $B$. In particular, we found $n=-0.147^{+0.113}_{-0.107}$ and $B=3.54 \times 10^{-3}$.

The limits of standard cosmography are here revised addressing the problem of error propagation during statistical analyses. To do so, we propose the use of Chebyshev polynomials to parameterize cosmic distances. In particular, we demonstrate that building up rational Chebyshev polynomials significantly reduces error propagations with respect to standard Taylor series. This technique provides unbiased estimations of the cosmographic parameters and performs significatively better than previous numerical approximations. To figure this out, we compare rational Chebyshev polynomials with Pad\'e series. In addition, we theoretically evaluate the convergence radius of (1,1) Chebyshev rational polynomial and we compare it with the convergence radii of Taylor and Pad\'e approximations. We thus focus on regions in which convergence of Chebyshev rational functions is better than standard approaches. With this recipe, as high-redshift data are employed, rational Chebyshev polynomials remain highly stable and enable one to derive highly accurate analytical approximations of Hubble's rate in terms of the cosmographic series. Finally, we check our theoretical predictions by setting bounds on cosmographic parameters through Monte Carlo integration techniques, based on the Metropolis-Hastings algorithm. We apply our technique to high-redshift cosmic data, using the JLA supernovae sample and the most recent versions of Hubble parameter and baryon acoustic oscillation measurements. We find that cosmography with Taylor series fails to be predictive with the aforementioned data sets, while turns out to be much more stable using the Chebyshev approach.

We propose a model-independent formalism to numerically solve the modified Friedmann equations in the framework of $f(T)$ teleparallel cosmology. Our strategy is to expand the Hubble parameter around the redshift $z=0$ up to a given order and to adopt cosmographic bounds as initial settings to determine the corresponding $f(z)\equiv f(T(H(z)))$ function. In this perspective, we distinguish two cases: the first expansion is up to the jerk parameter, the second expansion is up to the snap parameter. We show that inside the observed redshift domain $z\leq1$, only the net strength of $f(z)$ is modified passing from jerk to snap, whereas its functional behavior and shape turn out to be identical. As first step, we set the cosmographic parameters by means of the most recent observations. Afterwards, we calibrate our numerical solutions with the concordance $\Lambda$CDM model. In both cases, there is a good agreement with the cosmological standard model around $z\leq 1$, with severe discrepancies outer of this limit. We demonstrate that the effective dark energy term evolves following the test-function: $f(z)=\mathcal A+\mathcal Bz^2e^{\mathcal Cz}$. Bounds over the set $\left\{\mathcal A, \mathcal B, \mathcal C\right\}$ are also fixed by statistical considerations, comparing discrepancies between $f(z)$ with data. The approach opens the possibility to get a wide class of test-functions able to frame the dynamics of $f(T)$ without postulating any model \emph{a priori}. We thus re-obtain the $f(T)$ function through a back-scattering procedure once $f(z)$ is known. We figure out the properties of our $f(T)$ function at the level of background cosmology, to check the goodness of our numerical results. Finally, a comparison with previous cosmographic approaches is carried out giving results compatible with theoretical expectations.

QUBIC is an instrument aiming at measuring the B mode polarisation anisotropies at medium scales angular scales (30-200 multipoles). The search for the primordial CMB B-mode polarization signal is challenging, because of many difficulties: smallness of the expected signal, instrumental systematics that could possibly induce polarization leakage from the large E signal into B, brighter than anticipated polarized foregrounds (dust) reducing to zero the initial hope of finding sky regions clean enough to have a direct primordial B-modes observation. The QUBIC instrument is designed to address all aspects of this challenge with a novel kind of instrument, a Bolometric Interferometer, combining the background-limited sensitivity of Transition-Edge-Sensors and the control of systematics allowed by the observation of interference fringe patterns, while operating at two frequencies to disentangle polarized foregrounds from primordial B mode polarization. Its characteristics are described in details in this Technological Design Report.

A recent analysis of the Supernova Ia data claims a 'marginal' ($\sim3\sigma$) evidence for a cosmic acceleration. This result has been complemented with a non-accelerating $R_{h}=ct$ cosmology, which was presented as a valid alternative to the $\Lambda$CDM model. In this paper, we use the same analysis to show that a non-marginal evidence for acceleration is actually found. We compare the standard Friedmann models to the $R_{h}=ct$ cosmology by complementing SN Ia data with the Baryon Acoustic Oscillations, Gamma Ray Bursts and Observational Hubble datasets. We also study the power-law model which is a functional generalisation of $R_{h}=ct$. We find that the evidence for late-time acceleration is beyond refutable at a 4.56$\sigma$ confidence level from SN Ia data alone, and at an even stronger confidence level ($5.38\sigma$) from our joint analysis. Also, the non-accelerating $R_{h}=ct$ model fails to statistically compare with the $\Lambda$CDM having a $\Delta(\text{AIC})\sim30$.

We test the theoretical predictions of several cosmological models against different observables to compare the indirect estimates of the current expansion rate of the Universe determined from model fitting with the direct measurements based on Cepheids data published recently. We perform a statistical analysis of type Ia supernova (SN Ia), Hubble parameter, and baryon acoustic oscillation data. A joint analysis of these datasets allows us to better constrain cosmological parameters, but also to break the degeneracy that appears in the distance modulus definition between $H_0$ and the absolute B-band magnitude of SN Ia, $M_0$. From the theoretical side, we considered spatially flat and curvature-free $\Lambda$CDM, $w$CDM, and inhomogeneous Lema\^{i}tre-Tolman-Bondi (LTB) models. To analyse SN Ia we took into account the distributions of SN Ia intrinsic parameters. For the $\Lambda$CDM model we find that $\Omega_m=0.35\pm0.02$, $H_0=(67.8\pm1.0)\,$km$\,$s$^{-1}/$Mpc, while the corrected SN absolute magnitude has a normal distribution ${\cal N}(19.13,0.11)$. The $w$CDM model provides the same value for $\Omega_m$, while $H_0=(66.5\pm1.8)\,$km$\,$s$^{-1}/$Mpc and $w=-0.93\pm0.07$. When an inhomogeneous LTB model is considered, the combined fit provides $H_0=(64.2\pm1.9)\,$km$\,$s$^{-1}/$Mpc. Both the Akaike information criterion and the Bayes factor analysis cannot clearly distinguish between $\Lambda$CDM and $w$CDM cosmologies, while they clearly disfavour the LTB model. For the $\Lambda$CDM, our joint analysis of the SN Ia, the Hubble parameter, and the baryon acoustic oscillation datasets provides $H_0$ values that are consistent with CMB-only Planck measurements, but they differ by $2.5\sigma$ from the value based on Cepheids data.