We investigate primordial magnetogenesis and the evolution of the electromagnetic field through a quantum bounce, in a model that starts in the far past from a contracting phase where only dust is present and the electromagnetic field is in the adiabatic quantum vacuum state. By including a coupling between curvature and electromagnetism of the form $RF_{\mu \nu}F^{\mu \nu}$, we find acceptable magnetic field seeds within the current observational constraints at 1 mega-parsec (Mpc), and that the magnetic power spectrum evolves as a power-law with spectral index $n_B=6$. It is also shown that the electromagnetic backreaction is not an issue in the model under scrutiny.

An old question surrounding bouncing models concerns their stability under vector perturbations. Considering perfect fluids or scalar fields, vector perturbations evolve kinematically as $a^{-2}$, where $a$ is the scale factor. Consequently, a definite answer concerning the bounce stability depends on an arbitrary constant, therefore, there is no definitive answer. In this paper, we consider a more general situation where the primeval material medium is a non-ideal fluid, and its shear viscosity is capable of producing torque oscillations, which can create and dynamically sustain vector perturbations along cosmic evolution. In this framework, one can set that vector perturbations have a quantum mechanical origin, coming from quantum vacuum fluctuations in the far past of the bouncing model, as it is done with scalar and tensor perturbations. Under this prescription, one can calculate their evolution during the whole history of the bouncing model, and precisely infer the conditions under which they remain linear before the expanding phase. It is shown that such linearity conditions impose constraints on the free parameters of bouncing models, which are mild, although not trivial, allowing a large class of possibilities. Such conditions impose that vector perturbations are also not observationally relevant in the expanding phase. The conclusion is that bouncing models are generally stable under vector perturbations. As they are also stable under scalar and tensor perturbations, we conclude that bouncing models are generally stable under perturbations originated from quantum vacuum perturbations in the far past of their contracting phase.

The large amount of cosmological data already available (and in the near future) makes necessary the development of efficient numerical codes. Many software products have been implemented to perform cosmological analyses considering one or few probes. The need of multi-task software is rapidly increasing, in order to combine numerous cosmological probes along with their specificity (e.g., astrophysical descriptions and systematic errors). In this work we mention some of these libraries, bringing out some challenges they will face in the few-percent error era (on the cosmological parameters). We review some concepts of the standard cosmological model, and examine some specific topics on their implementation, bringing, for example, the discussion on how some quantities are numerically defined in different codes. We also consider implementation differences between public codes, mentioning their advantages/disadvantages.

We investigate a set of cosmological models for which the primordial power spectrum has a large-scale power deficit. The standard power-law spectrum is subject to long-wavelength modifications described by some new parameters, resulting in corrections to the anisotropies in the cosmic microwave background. The new parameters are fitted to different data sets: temperature only, temperature and polarization, the low-redshift determination of $H_0$, and baryonic acoustic oscillations. We discuss the statistical significance of the modified spectra, from both frequentist and Bayesian perspectives. Our analysis suggests motivations for considering models that break scalar-tensor consistency, or models with negligible power in the far super-Hubble limit. We present what appears to be substantial evidence for a new scale around 350 Mpc above which the primordial (scalar) power spectrum is sharply reduced by about 20%.

The attempt to describe the recent accelerated expansion of the Universe includes different propositions for dark energy models and modified gravity theories. Establish their features in order to discriminate and even rule out part of these models using observational data is a fundamental issue of cosmology. In the present work we consider a class of extended theories of gravity (ETG) that are minimally coupled to the ordinary matter fields. In this context, and assuming a homogeneous and isotropic spacetime, we derive the energy conditions for this ETG class. We then put constraints on these conditions using a model-independent approach to reconstruct the deceleration function along with the Joint Light-curve Analysis (JLA) supernova sample, 11 baryon acoustic oscillation and 22 cosmic-chronometer measurements. We also consider an additional condition imposing the strong energy condition only on the ordinary matter. This is to guarantee the presence of only attractive matter in the energy-momentum tensor, at least in the redshift range of the observations, i.e., the recent accelerated expansion of the Universe is due solely to the modifications in the gravity theory. The main result of this work is a general reconstruction of the energy conditions valid for every considered ETG.

The standard model of cosmology, {\Lambda}CDM, is the simplest model that matches the current observations, but it relies on two hypothetical components, to wit, dark matter and dark energy. Future galaxy surveys and cosmic microwave background (CMB) experiments will independently shed light on these components, but a joint analysis that includes cross-correlations will be necessary to extract as much information as possible from the observations. In this paper, we carry out a multi-probe analysis based on pseudo-spectra and test it on publicly available data sets. We use CMB temperature anisotropies and CMB lensing observations from Planck as well as the spectroscopic galaxy and quasar samples of SDSS-III/BOSS, taking advantage of the large areas covered by these surveys. We build a likelihood to simultaneously analyse the auto and cross spectra of CMB lensing and tracer overdensity maps before running Monte-Carlo Markov Chains (MCMC) to assess the constraining power of the combined analysis. We then add the CMB temperature anisotropies likelihood and obtain constraints on cosmological parameters ($H_0$, $\omega_b$, $\omega_c$, ${\ln10^{10}A_s}$, $n_s$ and $z_{re}$) and galaxy biases. We demonstrate that the joint analysis can additionally constrain the total mass of neutrinos ${\Sigma m_{\nu}}$ as well as the dark energy equation of state $w$ at once (for a total of eight cosmological parameters), which is impossible with either of the data sets considered separately. Finally, we discuss limitations of the analysis related to, e.g., the theoretical precision of the models, particularly in the non-linear regime.

Among the various possibilities to probe the theory behind the recent accelerated expansion of the universe, the energy conditions (ECs) are of particular interest, since it is possible to confront and constrain the many models, including different theories of gravity, with observational data. In this context, we use the ECs to probe any alternative theory whose extra term acts as a cosmological constant. For this purpose, we apply a model-independent approach to reconstruct the recent expansion of the universe. Using Type Ia supernova, baryon acoustic oscillations and cosmic-chronometer data, we perform a Markov Chain Monte Carlo analysis to put constraints on the effective cosmological constant $\Omega^0_{\rm eff}$. By imposing that the cosmological constant is the only component that possibly violates the ECs, we derive lower and upper bounds for its value. For instance, we obtain that $0.59 < \Omega^0_{\rm eff} < 0.91$ and $0.40 < \Omega^0_{\rm eff} < 0.93$ within, respectively, $1\sigma$ and $3\sigma$ confidence levels. In addition, about 30\% of the posterior distribution is incompatible with a cosmological constant, showing that this method can potentially rule it out as a mechanism for the accelerated expansion. We also study the consequence of these constraints for two particular formulations of the bimetric massive gravity. Namely, we consider the Visser's theory and the Hassan and Roses's massive gravity by choosing a background metric such that both theories mimic General Relativity with a cosmological constant. Using the $\Omega^0_{\rm eff}$ observational bounds along with the upper bounds on the graviton mass we obtain constraints on the parameter spaces of both theories.

Approximate Bayesian Computation (ABC) enables parameter inference for complex physical systems in cases where the true likelihood function is unknown, unavailable, or computationally too expensive. It relies on the forward simulation of mock data and comparison between observed and synthetic catalogues. Here we present cosmoabc, a Python ABC sampler featuring a Population Monte Carlo (PMC) variation of the original ABC algorithm, which uses an adaptive importance sampling scheme. The code is very flexible and can be easily coupled to an external simulator, while allowing to incorporate arbitrary distance and prior functions. As an example of practical application, we coupled cosmoabc with the numcosmo library and demonstrate how it can be used to estimate posterior probability distributions over cosmological parameters based on measurements of galaxy clusters number counts without computing the likelihood function. cosmoabc is published under the GPLv3 license on PyPI and GitHub and documentation is available at http://goo.gl/SmB8EX

We investigate scalar particle creation in a set of bouncing models where the bounce occurs due to quantum cosmological effects described by the Wheeler-DeWitt equation. The scalar field can be either conformally or minimally coupled to gravity, and it can be massive or massless, without self interaction. The analysis is made for models containing a single radiation fluid, and for the more realistic case of models containing the usual observed radiation and dust fluids, which can fit most of the observed features of our Universe, including an almost scale invariant power spectrum of scalar cosmological perturbations. In the conformal coupling case, the particle production is negligible. In the minimal coupling case, for massive particles, the results point to the same physical conclusion within observational constraints: particle production is most important at the bounce energy scale, and it is not sensitive neither to its mass nor whether there is dust in the background model. The only caveat is the case where the particle mass is larger than the bounce energy scale. On the other hand, the energy density of produced massive particles depend on their masses and the energy scale of the bounce. For very large masses and deep bounces, this energy density may overcome that of the background. In the case of massless particles, the energy density of produced particles can become comparable to the background energy density only for bounces occurring at energy scales comparable to the Planck scale or above, which lies beyond the scope of this paper: we expect that the simple Wheeler-DeWitt approach we are using should be valid only at scales some few orders of magnitude below the Planck energy. Nevertheless, in the case in which dust is present, there is an infrared divergence, which becomes important only for scales much larger than today's Hubble radius.

The formalism to treat quantization and evolution of cosmological perturbations of multiple fluids is described. We first construct the Lagrangian for both the gravitational and matter parts, providing the necessary relevant variables and momenta leading to the quadratic Hamiltonian describing linear perturbations. The final Hamiltonian is obtained without assuming any equations of motions for the background variables. This general formalism is applied to the special case of two fluids, having in mind the usual radiation and matter mix which made most of our current Universe history. Quantization is achieved using an adiabatic expansion of the basis functions. This allows for an unambiguous definition of a vacuum state up to the given adiabatic order. Using this basis, we show that particle creation is well defined for a suitable choice of vacuum and canonical variables, so that the time evolution of the corresponding quantum fields is unitary. This provides constraints for setting initial conditions for an arbitrary number of fluids and background time evolution. We also show that the common choice of variables for quantization can lead to an ill-defined vacuum definition. Our formalism is not restricted to the case where the coupling between fields is small, but is only required to vary adiabatically with respect to the ultraviolet modes, thus paving the way to consistent descriptions of general models not restricted to single-field (or fluid).

Distance measurements are currently the most powerful tool to study the expansion history of the universe without specifying its matter content nor any theory of gravitation. Assuming only an isotropic, homogeneous and flat universe, in this work we introduce a model-independent method to reconstruct directly the deceleration function via a piecewise function. Including a penalty factor, we are able to vary continuously the complexity of the deceleration function from a linear case to an arbitrary $(n+1)$-knots spline interpolation. We carry out a Monte Carlo (MC) analysis to determine the best penalty factor, evaluating the bias-variance trade-off, given the uncertainties of the SDSS-II and SNLS supernova combined sample (JLA), compilations of baryon acoustic oscillation (BAO) and $H(z)$ data. The bias-variance analysis is done for three fiducial models with different features in the deceleration curve. For each fiducial model, we test different reconstructions using, in each case, more than $10^4$ catalogs in a total of about $5\times 10^5$. This investigation proved to be essential in determining the best reconstruction to study these data. We show that, evaluating a single fiducial model, the conclusions about the bias-variance ratio are misleading. We determine the reconstruction method in which the bias represents at most $10\%$ of the total uncertainty. Finally, we apply the Ensemble Sampler Markov Chain Monte Carlo (ESMCMC) method to explore the posterior of the deceleration function up to redshift $1.3$ (using only JLA) and $2.3$ (JLA+BAO+$H(z)$). We obtain that the standard cosmological model agrees within $3\sigma$ level with the reconstructed results in the whole studied redshift intervals. Since our method is calibrated to minimize the bias, the error bars of the reconstructed functions are a good approximation for the total uncertainty.

Quantization of arbitrary free scalar fields in spatially homogeneous and isotropic space-times is considered. The quantum representation allowing a unitary evolution for the fields is taken as a requirement for the theory. Studying the group of linear canonical transformations, we show the relations between unitary evolution and choice of canonical variables. From these relations we obtain the conditions on the Hamiltonian such that there are canonical variables for which the field has unitary evolution. We then compute the linear transformation leading to these variables, also proving that they are unique. We obtain these results by developing the asymptotic analysis of the fields using the action angle variables, which proves to be a generalization of the usual Wentzel-Kramers-Brillouin approximation. These tools allow us to re-frame the adiabatic vacuum condition in a extensible format by using the action angle variables to relate these vacuum choices to those where the particle number density does not depend on angle (fast) variables. Finally, we develop a larger set of canonical variables relating the adiabatic vacuum conditions with the smearing of the quantum fields. This set of canonical variables also connects the adiabatic vacuum conditions with the instantaneous Hamiltonian diagonalization vacuum choice.

In a previous work we obtained a set of necessary conditions for the linear approximation in cosmology. Here we discuss the relations of this approach with the so called covariant perturbations. It is often argued in the literature that one of the main advantages of the covariant approach to describe cosmological perturbations is that the Bardeen formalism is coordinate dependent. In this paper we will reformulate the Bardeen approach in a completely covariant manner. For that, we introduce the notion of pure and mixed tensors, which yields an adequate language to treat both perturbative approaches in a common framework. We then stress that in the referred covariant approach one necessarily introduces an additional hyper-surface choice to the problem. Using our mixed and pure tensors approach, we were able to construct a one-to-one map relating the usual gauge dependence of the Bardeen formalism with the hyper-surface dependence inherent to the covariant approach. Finally, through the use of this map, we define full non-linear tensors that at first order correspond to the three known gauge invariant variables $\Phi$, $\Psi$ and $\Xi$, which are simultaneously foliation and gauge invariant. We then stress that the use of the proposed mixed tensors allows one to construct simultaneously gauge and hyper-surface invariant variables at any order.

A recent paper by Kumar (2012) (hereafter K12) claimed that in a contracting model, described by perturbations around a collapsing Friedmann model containing dust or radiation, the perturbations can grow in such a way that the linearity conditions would become invalid. This conclusion is not correct due to the following facts: first, it is claimed that the linearity conditions are not satisfied, but nowhere in K12 the amplitudes of the perturbations were in fact estimated. Therefore, without such estimates, the only possible conclusion from this work is the well known fact that the perturbations indeed grow during contraction, which, per se, does not imply that the linearity conditions become invalid. Second, some evaluations of the linearity conditions are incorrect because third other terms, instead of the appropriate second order ones, are mistakenly compared with first order terms, yielding artificially fast growing conditions. Finally, it is claimed that the results of K12 are in sharp contrast with the results of the paper by Vitenti and Pinto-Neto (2012) (hereafter VPN12), because the former was obtained in a gauge invariant way. However, the author of K12 did not realized that the evolution of the perturbations were also calculated in a gauge invariant way in VPN12, but some of the linearity conditions which are necessary to be checked cannot be expressed in terms of gauge invariant quantities. In the present work, the incorrect or incomplete statements of K12 are clarified and completed, and it is shown that all other correct results of K12 were already present in VPN12, whose conclusions remain untouched, namely, that cosmological perturbations of quantum mechanical origin in a bouncing model can remain in the linear regime all along the contracting phase and at the bounce itself for a wide interval of energy scales of the bounce. (Abstract abridged)

In previous works, it was shown that the Lagrangians and Hamiltonians of cosmological linear scalar, vector and tensor perturbations of homogeneous and isotropic space-times with flat spatial sections containing a perfect fluid can be put in a simple form through the implementation of canonical transformations and redefinitions of the lapse function, without ever using the background classical equations of motion. In this paper, we generalize this result to general fluids, which includes entropy perturbations, and to arbitrary spacelike hyper-surfaces through a new method together with the Faddeev-Jackiw procedure for the constraint reduction. A simple second order Hamiltonian involving the Mukhanov-Sasaki variable is obtained, again without ever using the background equations of motion.

It has been shown that a contracting universe with a dust-like ($w \approx 0$) fluid may provide an almost scale invariant spectrum for the gravitational scalar perturbations. As the universe contracts, the amplitude of such perturbations are amplified. The gauge invariant variable $\Phi$ develops a growing mode which becomes much larger than the constant one around the bounce phase. The constant mode has its amplitude fixed by Cosmic Background Explorer (COBE) normalization, thus the amplitude of the growing mode can become much larger than 1. In this paper, we first show that this is a general feature of bouncing models, since we expect that general relativity should be valid in all scales away from the bounce. However, in the Newtonian gauge, the variable $\Phi$ gives the value of the metric perturbation $\phi$, raising doubts on the validity of the linear perturbative regime at the bounce. In order to address this issue, we obtain a set of necessary conditions for the perturbative series to be valid along the whole history of the model, and we show that there is a gauge in which all these conditions are satisfied, for a set of models, if the constant mode is fixed by COBE normalization. As a by-product of this analysis, we point out that there are sets of solutions for the perturbation variables where some gauge-fixing conditions are not well defined, turning these gauges prohibited for those solutions.

The immediate observational consequence of a non-trivial spatial topology of the Universe is that an observer could potentially detect multiple images of radiating sources. In particular, a non-trivial topology will generate pairs of correlated circles of temperature fluctuations in the anisotropies maps of the cosmic microwave background (CMB), the so-called circles-in-the-sky. In this way, a detectable non-trivial spatial topology may be seen as an observable attribute, which can be probed through the circles-in-the-sky for all locally homogeneous and isotropic universes with no assumptions on the cosmological dark energy (DE) equation of state (EOS) parameters. We show that the knowledge of the spatial topology through the circles-in-the-sky offers an effective way of reducing the degeneracies in the DE EOS parameters. We concretely illustrate the topological role by assuming, as an exanple, a Poincar\'{e} dodecahedral space topology and reanalyzing the constraints on the parameters of a specific EOS which arise from the supernovae type Ia, baryon acoustic oscillations and the CMB plus the statistical topological contribution.

In the standard Friedmann-Lemaitre-Robertson-Walker (FLRW) approach to model the Universe the violation of the so-called energy conditions is related to some important properties of the Universe as, for example, the current and the inflationary accelerating expansion phases. The energy conditions are also necessary in the formulation and proofs of Hawking-Penrose singularity theorems. In two recent articles we have derived bounds from energy conditions and made confrontations of these bounds with supernovae data. Here, we extend these results in following way: first, by using our most recent statistical procedure for calculating new q(z) estimates from the \emph{gold} and \emph{combined} type Ia supernovae samples; second, we use these estimates to obtain a new picture of the energy conditions fulfillment and violation for the recent past ($z\leq 1 $) in the context of the standard cosmology.

The energy conditions play an important role in the description of some important properties of the Universe, including the current accelerating expansion phase and the possible recent phase of super-acceleration. In a recent work we have provided a detailed study of the energy conditions for the recent past by deriving bounds from energy conditions and by making the confrontation of the bounds with supernovae data. Here, we extend and update these results in two different ways. First, by carrying out a new statistical analysis for $q(z)$ estimates needed for the confrontation between the bounds and supernovae data. Second, by providing a new picture of the energy conditions fulfillment and violation in the light of the recently compiled Union set of 307 type Ia supernovae and by using two different statistical approaches.

The energy conditions play an important role in the understanding of several properties of the Universe, including the current accelerating expansion phase and the possible existence of the so-called phantom fields. We show that the integrated bounds provided by the energy conditions on cosmological observables such as the distance modulus $\mu(z)$ and the lookback time $t_L(z)$ are not sufficient (nor necessary) to ensure the local fulfillment of the energy conditions, making explicit the limitation of these bounds in the confrontation with observational data. We recast the energy conditions as bounds on the deceleration and normalized Hubble parameters, obtaining new bounds which are necessary and sufficient for the local fulfillment of the energy conditions. A statistical confrontation, with $1\sigma-3\sigma$ confidence levels, between our bounds and supernovae data from the gold and combined samples is made for the recent past. Our analyses indicate, with $3\sigma$ confidence levels, the fulfillment of both the weak energy condition (WEC) and dominant energy condition (DEC) for $z \leq 1$ and $z \lesssim 0.8$, respectively. In addition, they suggest a possible recent violation of the null energy condition (NEC) with $3\sigma$, i.e. a very recent phase of super-acceleration. Our analyses also show the possibility of violation of the strong energy condition (\textbf{SEC}) with $3\sigma$ in the recent past ($z \leq 1$), but interestingly the $q(z)$-best-fit curve crosses the SEC-fulfillment divider at $z \simeq 0.67$, which is a value very close to the beginning of the epoch of cosmic acceleration predicted by the standard concordance flat $\Lambda$CDM scenario.