The development of numerical N -body simulations have allowed to study formation process and evolution of galaxies at different scales. This paper presents the fundamental concepts of N-body systems applied to the cosmological evolution of the $\Lambda$-Cold Dark Matter ($\Lambda$CDM) model. In order to perform structure formation in the Universe, we provide an introduction to the basic equations and their implementation on the GADGET-2 software. We also present a simple guide to modify this code. First, we briefly describe the dark matter in the Universe as well as the theoretical and experimental basis of the $\Lambda$CDM model. Then, we focus on the simulation codes and provide the equations that govern most of the N-body simulations to model the dark matter. We describe the Smoothed Particle Hydrodynamics method used for simulating the gas, star dynamics and structure formation in these simulations. Then, cautiously, we guide the reader to the installation of GADGET-2 on a Linux-based computer, as well as to carry out a couple of examples to operate the code. Finally, by using a computational cluster, we show several results of a large structure simulation, analyse the outputs to display the matter power spectrum, and compare the outcome with theoretical predictions.

We study the cosmological constant ($\Lambda$) in the standard $\Lambda$CDM model by introducing the \textit{graduated dark energy} (gDE) characterised by a minimal dynamical deviation from the null inertial mass density of the $\Lambda$ in the form $\rho_{\rm inert}\propto \rho^{\lambda}<0$ with $\lambda<1$ being a ratio of two odd integers, for which its energy density $\rho$ dynamically takes negative values in the finite past. For large negative values of $\lambda$, it creates a phenomenological model described by a smooth function that approximately describes the $\Lambda$ spontaneously switching sign in the late universe to become positive today. We confront the model with the latest combined observational data sets of PLK+BAO+SN+$H$. It is striking that the data predict bimodal posterior probability distributions for the parameters of the model along with large negative $\lambda$ values; the new maximum significantly excludes the $\Lambda$, and the old maximum contains the $\Lambda$. The improvement in the goodness of fit for the $\Lambda$ reaches highly significant levels, $\Delta\chi_{\rm min}^2=6.4$ for the new maxima, while it remains at insignificant levels, $\Delta\chi_{\rm min}^2\lesssim0.02$, for the old maxima. We show that, in contrast to the old maxima, which do not distinguish from the $\Lambda$, the new maxima agree with the model-independent $H_0$ measurements, high-precision Ly-$\alpha$ data, and model-independent $Omh^2$ diagnostic estimates. Our results provide strong hints of a spontaneous sign switch in the cosmological constant and lead us to conjecture that the universe has transitioned from AdS vacua to dS vacua, at a redshift $z\approx 2.32$ and triggered the late-time acceleration, and suggests looking for such mechanisms in string theory constructions.

We present an explicit detailed theoretical and observational investigation of an anisotropic massive Brans-Dicke (BD) gravity extension of the standard $\Lambda$CDM model, wherein the extension is characterized by two additional degrees of freedom; the BD parameter, $\omega$, and the present day density parameter corresponding to the shear scalar, $\Omega_{\sigma^2,0}$. The BD parameter, determining the deviation from general relativity (GR), by alone characterizes both the dynamics of the effective dark energy (DE) and the redshift dependence of the shear scalar. These two affect each other depending on $\omega$, namely, the shear scalar contributes to the dynamics of the effective DE, and its anisotropic stress --which does not exist in scalar field models of DE within GR-- controls the dynamics of the shear scalar deviating from the usual $\propto(1+z)^6$ form in GR. We mainly confine the current work to non-negative $\omega$ values as it is the right sign --theoretically and observationally-- for investigating the model as a correction to the $\Lambda$CDM. By considering the current cosmological observations, we find that $\omega\gtrsim 250$, $\Omega_{\sigma^2,0}\lesssim 10^{-23}$ and the contribution of the anisotropy of the effective DE to this value is insignificant. We conclude that the simplest anisotropic massive BD gravity extension of the standard $\Lambda$CDM model exhibits no significant deviations from it all the way to the Big Bang Nucleosynthesis. We also point out the interesting features of the model in the case of negative $\omega$ values; for instance, the constraints on $\Omega_{\sigma^2,0}$ could be relaxed considerably, the values of $\omega\sim-1$ (relevant to string theories) predict dramatically different dynamics for the expansion anisotropy.

We study a new model of Energy-Momentum Squared Gravity (EMSG), called Energy-Momentum Log Gravity (EMLG), constructed by the addition of the term $f(T_{\mu\nu}T^{\mu\nu})=\alpha \ln(\lambda\,T_{\mu\nu}T^{\mu\nu})$, envisaged as a correction, to the Einstein-Hilbert action with cosmological constant $\Lambda$. The choice of this modification is made as a specific way of including new terms in the right-hand side of the Einstein field equations, resulting in constant effective inertial mass density and, importantly, leading to an explicit exact solution of the matter energy density in terms of redshift. We look for viable cosmologies, in particular, an extension of the standard $\Lambda$CDM model. EMLG provides an effective dynamical dark energy passing below zero at large redshifts, accommodating a mechanism for screening $\Lambda$ in this region, in line with suggestions for alleviating some of the tensions that arise between observational data sets within the standard $\Lambda$CDM model. We present a detailed theoretical investigation of the model and then constrain the free parameter $\alpha'$, a normalisation of $\alpha$, using the latest observational data. The data does not rule out the $\Lambda$CDM limit of our model ($\alpha'= 0$), but prefers slightly negative values of the EMLG model parameter ($\alpha'= -0.032\pm 0.043$), which leads to the screening of $\Lambda$. We also discuss how EMLG relaxes the persistent tension that appears in the measurements of $H_0$ within the standard $\Lambda$CDM model.

The current concordance model of cosmology is dominated by two mysterious ingredients: dark matter and dark energy. In this paper, we explore the possibility that, in fact, there exist two dark-energy components: the cosmological constant $\Lambda$, with equation-of-state parameter $w_\Lambda=-1$, and a `missing matter' component $X$ with $w_X=-2/3$, which we introduce here to allow the evolution of the universal scale factor as a function of conformal time to exhibit a symmetry that relates the big bang to the future conformal singularity, such as in Penrose's conformal cyclic cosmology. Using recent cosmological observations, we constrain the present-day energy density of missing matter to be $\Omega_{X,0}=-0.034 \pm 0.075$. This is consistent with the standard $\Lambda$CDM model, but constraints on the energy densities of all the components are considerably broadened by the introduction of missing matter; significant relative probability exists even for $\Omega_{X,0} \sim 0.1$, and so the presence of a missing matter component cannot be ruled out. As a result, a Bayesian model selection analysis only slightly disfavours its introduction by 1.1 log-units of evidence. Foregoing our symmetry requirement on the conformal time evolution of the universe, we extend our analysis by allowing $w_X$ to be a free parameter. For this more generic `double dark energy' model, we find $w_X = -1.01 \pm 0.16$ and $\Omega_{X,0} = -0.10 \pm 0.56$, which is again consistent with the standard $\Lambda$CDM model, although once more the posterior distributions are sufficiently broad that the existence of a second dark-energy component cannot be ruled out. The model including the second dark energy component also has an equivalent Bayesian evidence to $\Lambda$CDM, within the estimation error, and is indistinguishable according to the Jeffreys guideline.

The dark energy component of the universe still remains as a mystery, however, several papers based on observational data have shown that its equation of state may have an oscillatory behaviour. In this paper, we provide a general description for the dark-energy equation-of-state $w(z)$ in the form of Fourier series. This description generalises some previous dynamical dark energy models and is in agreement with the $w(z)$ reconstructions. We make use of a modified version of a simple and fast Markov Chain Monte Carlo code to constraint the model parameters. For the analysis we use data from supernovae type-Ia , baryon acoustic oscillations, $H(z)$ measurements and cosmic microwave background. We provide a comparison of the proposed model with $\Lambda$CDM, $w$CDM and the standard Taylor approximation. The Fourier series expansion of $w(z)$ is preferred from $\Lambda$CDM at more than $3\sigma$ significance level based on the improvement in the fit alone. We use the Akaike criteria to perform the model comparison and found that, even though there are extra parameters, there is a slight preference of the Fourier series compared with the $\Lambda$CDM model. The preferred shape of $w(z)$ found here puts in jeopardy the single scalar field models, as they as they cannot reproduce the crossing the phantom divide line $w=-1$.

Bayesian statistics and MCMC (Markov Chain Monte Carlo) algorithms have found their place in the field of Cosmology. They have become important mathematical and numerical tools, especially in parameter estimation and model comparison. In this paper we review some of the fundamental concepts to understand Bayesian statistics, to then introduce the MCMC algorithms and samplers that allow us to perform the parameter inference procedure. We also provide a general description of the standard cosmological model, known as the $\Lambda$CDM model, along with several alternatives to it; and current datasets coming from astrophysical and cosmological observations. Finally, with the tools acquired we use a MCMC algorithm implemented in python -called SimpleMC- to test the cosmological models and find out the combination of parameters that best describes the universe.

The aim of the Air Microwave Yield (AMY) experiment is to investigate the Molecular Bremsstrahlung Radiation (MBR) emitted from an electron beam induced air-shower. The measurements have been performed with a 510 MeV electron beam at the Beam Test Facility (BTF) of Frascati INFN National Laboratories in a wide frequency range between 1 and 20 GHz. We present the experimental apparatus and the first results of the measurements. Contrary to what have been reported in a previous similar experiment~\cite{Gorham-SLAC}, we have found that the intensity of the emission is strongly influenced by the particular time structure of the accelerator beam. This makes very difficult the interpretation of the emission process and a realistic extrapolation of the emission yield to the plasma generated during the development of an atmospheric shower.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

A flat Friedman-Roberson-Walker universe dominated by a cosmological constant ($\Lambda$) and cold dark matter (CDM) has been the working model preferred by cosmologists since the discovery of cosmic acceleration. However, tensions of various degrees of significance are known to be present among existing datasets within the $\Lambda$CDM framework. In particular, the Lyman-$\alpha$ forest measurement of the Baryon Acoustic Oscillations (BAO) by the Baryon Oscillation Spectroscopic Survey (BOSS) prefers a smaller value of the matter density fraction $\Omega_{\rm M}$ compared to the value preferred by cosmic microwave background (CMB). Also, the recently measured value of the Hubble constant, $H_0=73.24\pm1.74 \ {\rm km}\ {\rm s}^{-1} \ {\rm Mpc}^{-1}$, is $3.4\sigma$ higher than $66.93\pm0.62 \ {\rm km}\ {\rm s}^{-1} \ {\rm Mpc}^{-1}$ inferred from the Planck CMB data. In this work, we investigate if these tensions can be interpreted as evidence for a non-constant dynamical dark energy (DE). Using the Kullback-Leibler (KL) divergence to quantify the tension between datasets, we find that the tensions are relieved by an evolving DE, with the dynamical DE model preferred at a $3.5\sigma$ significance level based on the improvement in the fit alone. While, at present, the Bayesian evidence for the dynamical DE is insufficient to favour it over $\Lambda$CDM, we show that, if the current best fit DE happened to be the true model, it would be decisively detected by the upcoming DESI survey.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

We use flux-transmission correlations in \Lya forests to measure the imprint of baryon acoustic oscillations (BAO). The study uses spectra of 157,783 quasars in the redshift range $2.1\le z \le 3.5$ from the Sloan Digital Sky Survey (SDSS) Data Release 12 (DR12). Besides the statistical improvements on our previous studies using SDSS DR9 and DR11, we have implemented numerous improvements in the analysis procedure, allowing us to construct a physical model of the correlation function and to investigate potential systematic errors in the determination of the BAO peak position. The Hubble distance, $\DHub=c/H(z)$, relative to the sound horizon is $\DHub(z=2.33)/r_d=9.07 \pm 0.31$. The best-determined combination of comoving angular-diameter distance, $\DM$, and the Hubble distance is found to be $\DHub^{0.7}\DM^{0.3}/r_d=13.94\pm0.35$. This value is $1.028\pm0.026$ times the prediction of the flat-\lcdm model consistent with the cosmic microwave background (CMB) anisotropy spectrum. The errors include marginalization over the effects of unidentified high-density absorption systems and fluctuations in ultraviolet ionizing radiation. Independently of the CMB measurements, the combination of our results and other BAO observations determine the open-\lcdm density parameters to be $\om=0.296 \pm 0.029$, $\ol=0.699 \pm 0.100$ and $\Omega_k = -0.002 \pm 0.119$.

We analyse the Baryon Acoustic Oscillation (BAO) signal of the final Baryon Oscillation Spectroscopic Survey (BOSS) data release (DR12). Our analysis is performed in Fourier-space, using the power spectrum monopole and quadrupole. The dataset includes $1\,198\,006$ galaxies over the redshift range $0.2 < z < 0.75$. We divide this dataset into three (overlapping) redshift bins with the effective redshifts $\zeff = 0.38$, $0.51$ and $0.61$. We demonstrate the reliability of our analysis pipeline using N-body simulations as well as $\sim 1000$ MultiDark-Patchy mock catalogues, which mimic the BOSS-DR12 target selection. We apply density field reconstruction to enhance the BAO signal-to-noise ratio. By including the power spectrum quadrupole we can separate the line-of-sight and angular modes, which allows us to constrain the angular diameter distance $D_A(z)$ and the Hubble parameter $H(z)$ separately. We obtain two independent $1.6\%$ and $1.5\%$ constraints on $D_A(z)$ and $2.9\%$ and $2.3\%$ constraints on $H(z)$ for the low ($\zeff=0.38$) and high ($\zeff=0.61$) redshift bin, respectively. We obtain two independent $1\%$ and $0.9\%$ constraints on the angular averaged distance $D_V(z)$, when ignoring the Alcock-Paczynski effect. The detection significance of the BAO signal is of the order of $8\sigma$ (post-reconstruction) for each of the three redshift bins. Our results are in good agreement with the Planck prediction within $\Lambda$CDM. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are combined with others in~\citet{Alam2016} to produce the final cosmological constraints from BOSS.

We present cosmological results from the final galaxy clustering data set of the Baryon Oscillation Spectroscopic Survey, part of the Sloan Digital Sky Survey III. Our combined galaxy sample comprises 1.2 million massive galaxies over an effective area of 9329 deg^2 and volume of 18.7 Gpc^3, divided into three partially overlapping redshift slices centred at effective redshifts 0.38, 0.51, and 0.61. We measure the angular diameter distance DM and Hubble parameter H from the baryon acoustic oscillation (BAO) method after applying reconstruction to reduce non-linear effects on the BAO feature. Using the anisotropic clustering of the pre-reconstruction density field, we measure the product DM*H from the Alcock-Paczynski (AP) effect and the growth of structure, quantified by f{\sigma}8(z), from redshift-space distortions (RSD). We combine measurements presented in seven companion papers into a set of consensus values and likelihoods, obtaining constraints that are tighter and more robust than those from any one method. Combined with Planck 2015 cosmic microwave background measurements, our distance scale measurements simultaneously imply curvature {\Omega}_K =0.0003+/-0.0026 and a dark energy equation of state parameter w = -1.01+/-0.06, in strong affirmation of the spatially flat cold dark matter model with a cosmological constant ({\Lambda}CDM). Our RSD measurements of f{\sigma}_8, at 6 per cent precision, are similarly consistent with this model. When combined with supernova Ia data, we find H0 = 67.3+/-1.0 km/s/Mpc even for our most general dark energy model, in tension with some direct measurements. Adding extra relativistic species as a degree of freedom loosens the constraint only slightly, to H0 = 67.8+/-1.2 km/s/Mpc. Assuming flat {\Lambda}CDM we find {\Omega}_m = 0.310+/-0.005 and H0 = 67.6+/-0.5 km/s/Mpc, and we find a 95% upper limit of 0.16 eV/c^2 on the neutrino mass sum.

We analyse the broad-range shape of the monopole and quadrupole correlation functions of the BOSS Data Release 12 (DR12) CMASS and LOWZ galaxy sample to obtain constraints on the Hubble expansion rate $H(z)$, the angular-diameter distance $D_A(z)$, the normalised growth rate $f(z)\sigma_8(z)$, and the physical matter density $\Omega_mh^2$. We adopt wide and flat priors on all model parameters in order to ensure the results are those of a `single-probe' galaxy clustering analysis. We also marginalise over three nuisance terms that account for potential observational systematics affecting the measured monopole. However, such Monte Carlo Markov Chain analysis is computationally expensive for advanced theoretical models, thus we develop a new methodology to speed up our analysis. We obtain $\{D_A(z)r_{s,fid}/r_s$Mpc, $H(z)r_s/r_{s,fid}$kms$^{-1}$Mpc$^{-1}$, $f(z)\sigma_8(z)$, $\Omega_m h^2\}$ = $\{956\pm28$ , $75.0\pm4.0$ , $0.397 \pm 0.073$, $0.143\pm0.017\}$ at $z=0.32$ and $\{1421\pm23$, $96.7\pm2.7$ , $0.497 \pm 0.058$, $0.137\pm0.015\}$ at $z=0.59$ where $r_s$ is the comoving sound horizon at the drag epoch and $r_{s,fid}=147.66$Mpc for the fiducial cosmology in this study. In addition, we divide the galaxy sample into four redshift bins to increase the sensitivity of redshift evolution. However, we do not find improvements in terms of constraining dark energy model parameters. Combining our measurements with Planck data, we obtain $\Omega_m=0.306\pm0.009$, $H_0=67.9\pm0.7$kms$^{-1}$Mpc$^{-1}$, and $\sigma_8=0.815\pm0.009$ assuming $\Lambda$CDM; $\Omega_k=0.000\pm0.003$ assuming oCDM; $w=-1.01\pm0.06$ assuming $w$CDM; and $w_0=-0.95\pm0.22$ and $w_a=-0.22\pm0.63$ assuming $w_0w_a$CDM. Our results show no tension with the flat $\Lambda$CDM cosmological paradigm. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS.

We develop a new methodology called double-probe analysis with the aim of minimizing informative priors in the estimation of cosmological parameters. We extract the dark-energy-model-independent cosmological constraints from the joint data sets of Baryon Oscillation Spectroscopic Survey (BOSS) galaxy sample and Planck cosmic microwave background (CMB) measurement. We measure the mean values and covariance matrix of $\{R$, $l_a$, $\Omega_b h^2$, $n_s$, $log(A_s)$, $\Omega_k$, $H(z)$, $D_A(z)$, $f(z)\sigma_8(z)\}$, which give an efficient summary of Planck data and 2-point statistics from BOSS galaxy sample, where $R=\sqrt{\Omega_m H_0^2}\,r(z_*)$, and $l_a=\pi r(z_*)/r_s(z_*)$, $z_*$ is the redshift at the last scattering surface, and $r(z_*)$ and $r_s(z_*)$ denote our comoving distance to $z_*$ and sound horizon at $z_*$ respectively. The advantage of this method is that we do not need to put informative priors on the cosmological parameters that galaxy clustering is not able to constrain well, i.e. $\Omega_b h^2$ and $n_s$. Using our double-probe results, we obtain $\Omega_m=0.304\pm0.009$, $H_0=68.2\pm0.7$, and $\sigma_8=0.806\pm0.014$ assuming $\Lambda$CDM; and $\Omega_k=0.002\pm0.003$ and $w=-1.00\pm0.07$ assuming o$w$CDM. The results show no tension with the flat $\Lambda$CDM cosmological paradigm. By comparing with the full-likelihood analyses with fixed dark energy models, we demonstrate that the double-probe method provides robust cosmological parameter constraints which can be conveniently used to study dark energy models. We extend our study to measure the sum of neutrino mass and obtain $\Sigma m_\nu<0.10/0.22$ (68\%/95\%) assuming $\Lambda$CDM and $\Sigma m_\nu<0.26/0.52$ (68\%/95\%) assuming $w$CDM. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS.

The cosmological information contained in anisotropic galaxy clustering measurements can often be compressed into a small number of parameters whose posterior distribution is well described by a Gaussian. We present a general methodology to combine these estimates into a single set of consensus constraints that encode the total information of the individual measurements, taking into account the full covariance between the different methods. We illustrate this technique by applying it to combine the results obtained from different clustering analyses, including measurements of the signature of baryon acoustic oscillations (BAO) and redshift-space distortions (RSD), based on a set of mock catalogues of the final SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). Our results show that the region of the parameter space allowed by the consensus constraints is smaller than that of the individual methods, highlighting the importance of performing multiple analyses on galaxy surveys even when the measurements are highly correlated. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The methodology presented here is used in Alam et al. (2016) to produce the final cosmological constraints from BOSS.

We construct two simple effective field theory versions of {\it Hybrid Natural Inflation (HNI)} that illustrate the range of its phenomenological implications. The resulting inflationary sector potential, $V=\Delta^4(1+a\cos(\phi/f))$, arises naturally, with the inflaton field a pseudo-Nambu-Goldstone boson. The end of inflation is triggered by a waterfall field and the conditions for this to happen are determined. Also of interest is the fact that the slow-roll parameter $\epsilon$ (and hence the tensor $r$) is a non-monotonic function of the field with a maximum where observables take universal values that determines the maximum possible tensor to scalar ratio $r$. In one of the models the inflationary scale can be as low as the electroweak scale. We explore in detail the associated HNI phenomenology, taking account of the constraints from Black Hole production, and perform a detailed fit to the Planck 2015 temperature and polarisation data.

In recent years, the Lyman-$\alpha$ absorption observed in the spectra of high-redshift quasars has been used as a tracer of large-scale structure by means of the three-dimensional Lyman-$\alpha$ forest auto-correlation function at redshift $z\simeq 2.3$, but the need to fit the quasar continuum in every absorption spectrum introduces a broadband distortion that is difficult to correct and causes a systematic error for measuring any broadband properties. We describe a $k$-space model for this broadband distortion based on a multiplicative correction to the power spectrum of the transmitted flux fraction that suppresses power on scales corresponding to the typical length of a Lyman-$\alpha$ forest spectrum. Implementing the distortion model in fits for the baryon acoustic oscillation (BAO) peak position in the Lyman-$\alpha$ forest auto-correlation, we find that the fitting method recovers the input values of the linear bias parameter $b_{F}$ and the redshift-space distortion parameter $\beta_{F}$ for mock data sets with a systematic error of less than 0.5\%. Applied to the auto-correlation measured for BOSS Data Release 11, our method improves on the previous treatment of broadband distortions in BAO fitting by providing a better fit to the data using fewer parameters and reducing the statistical errors on $\beta_{F}$ and the combination $b_{F}(1+\beta_{F})$ by more than a factor of seven. The measured values at redshift $z=2.3$ are $\beta_{F}=1.39^{+0.11\ +0.24\ +0.38}_{-0.10\ -0.19\ -0.28}$ and $b_{F}(1+\beta_{F})=-0.374^{+0.007\ +0.013\ +0.020}_{-0.007\ -0.014\ -0.022}$ (1$\sigma$, 2$\sigma$ and 3$\sigma$ statistical errors). Our fitting software and the input files needed to reproduce our main results are publicly available.

We derive constraints on cosmological parameters and tests of dark energy models from the combination of baryon acoustic oscillation (BAO) measurements with cosmic microwave background (CMB) and Type Ia supernova (SN) data. We take advantage of high-precision BAO measurements from galaxy clustering and the Ly-alpha forest (LyaF) in the BOSS survey of SDSS-III. BAO data alone yield a high confidence detection of dark energy, and in combination with the CMB angular acoustic scale they further imply a nearly flat universe. Combining BAO and SN data into an "inverse distance ladder" yields a 1.7% measurement of $H_0=67.3 \pm1.1$ km/s/Mpc. This measurement assumes standard pre-recombination physics but is insensitive to assumptions about dark energy or space curvature, so agreement with CMB-based estimates that assume a flat LCDM cosmology is an important corroboration of this minimal cosmological model. For open LCDM, our BAO+SN+CMB combination yields $\Omega_m=0.301 \pm 0.008$ and curvature $\Omega_k=-0.003 \pm 0.003$. When we allow more general forms of evolving dark energy, the BAO+SN+CMB parameter constraints remain consistent with flat LCDM. While the overall $\chi^2$ of model fits is satisfactory, the LyaF BAO measurements are in moderate (2-2.5 sigma) tension with model predictions. Models with early dark energy that tracks the dominant energy component at high redshifts remain consistent with our constraints. Expansion history alone yields an upper limit of 0.56 eV on the summed mass of neutrino species, improving to 0.26 eV if we include Planck CMB lensing. Standard dark energy models constrained by our data predict a level of matter clustering that is high compared to most, but not all, observational estimates. (Abridged)

We introduce a new parametrization for the dark energy, led by the same idea to the linear expansion of the equation of state in scale factor $a$ and in redshift $z$, which diverges neither in the past nor future and contains the same number of degrees of freedom with the former two. We present constraints of the cosmological parameters using the most updated baryon acoustic oscillation (BAO) measurements along with cosmic microwave background (CMB) data and a recent reanalysis of Type Ia supernova (SN) data. This new parametrization allowed us to carry out successive observational analyses by decreasing its degrees of freedom systematically until ending up with a dynamical dark energy model that has the same number of parameters with $\Lambda$CDM. We found that the dark energy source with a dynamical equation of state parameter equal $-2/3$ at the early universe and $-1$ today fits the data slightly better than $\Lambda$.

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

(Abridged) We detect the large-scale structure of Lya emission in the Universe at redshifts z=2-3.5 by measuring the cross-correlation of Lya surface brightness with quasars in SDSS/BOSS. We use a million spectra targeting Luminous Red Galaxies at z<0.8, after subtracting a best fit model galaxy spectrum from each one, as an estimate of the high-redshift Lya surface brightness. The quasar-Lya emission cross-correlation we detect has a shape consistent with a LambdaCDM model with Omega_M =0.30^+0.10-0.07. The predicted amplitude of this cross-correlation is proportional to the product of the mean Lya surface brightness, <mu_alpha>, the amplitude of mass fluctuations, and the quasar and Lya emission bias factors. Using known values, we infer <mu_alpha>(b_alpha/3) = (3.9 +/- 0.9) x 10^-21 erg/s cm^-2 A^-1 arcsec^-2, where b_alpha is the Lya emission bias factor. If the dominant sources of Lya emission are star forming galaxies, we infer rho_SFR = (0.28 +/- 0.07) (3/b_alpha) /yr/Mpc^3 at z=2-3.5. For b_alpha=3, this value is a factor of 21-35 above previous estimates from individually detected Lya emitters, although consistent with the total rho_SFR derived from dust-corrected, continuum UV surveys. 97% of the Lya emission in the Universe at these redshifts is therefore undetected in previous surveys of Lya emitters. Our measurement is much greater than seen from stacking analyses of faint halos surrounding previously detected Lya emitters, but we speculate that it arises from similar Lya halos surrounding all luminous star-forming galaxies. We also detect redshift space anisotropy of the quasar-Lya emission cross-correlation, finding evidence at the 3.0 sigma level that it is radially elongated, consistent with distortions caused by radiative-transfer effects (Zheng et al. (2011)). Our measurements represent the first application of the intensity mapping technique to optical observations.

We study the Hybrid Natural Inflation (HNI) model and some of its realisations in the light of recent CMB observations, mainly Planck temperature and WMAP-9 polarization, and compare with the recent release of BICEP2 dataset. The inflationary sector of HNI is essentially given by the potential $V(\phi) = V_0(1+a\cos (\frac{\phi}{f} ) )$, where $a$ is a positive constant smaller or equal to one and $f$ is the scale of (pseudo Nambu-Goldstone) symmetry breaking. We show that to describe the HNI model realisations we only need two observables; the spectral index $n_s$, the tensor-to-scalar ratio, and a free parameter in the amplitude of the cosine function $a$. We find that in order to make the HNI model compatible with the BICEP2 observations, we require a large positive running of the spectra. We find that this could over-produce primordial black holes in the most consistent case of the model. This situation could be aleviated if, as recently argued, the BICEP2 data do not correspond to primordial gravitational waves.

The Pierre Auger Observatory is the world's largest cosmic ray observatory. Our current exposure reaches nearly 40,000 km$^2$ str and provides us with an unprecedented quality data set. The performance and stability of the detectors and their enhancements are described. Data analyses have led to a number of major breakthroughs. Among these we discuss the energy spectrum and the searches for large-scale anisotropies. We present analyses of our X$_{max}$ data and show how it can be interpreted in terms of mass composition. We also describe some new analyses that extract mass sensitive parameters from the 100% duty cycle SD data. A coherent interpretation of all these recent results opens new directions. The consequences regarding the cosmic ray composition and the properties of UHECR sources are briefly discussed.

The AMY experiment aims to measure the microwave bremsstrahlung radiation (MBR) emitted by air-showers secondary electrons accelerating in collisions with neutral molecules of the atmosphere. The measurements are performed using a beam of 510 MeV electrons at the Beam Test Facility (BTF) of Frascati INFN National Laboratories. The goal of the AMY experiment is to measure in laboratory conditions the yield and the spectrum of the GHz emission in the frequency range between 1 and 20 GHz. The final purpose is to characterise the process to be used in a next generation detectors of ultra-high energy cosmic rays. A description of the experimental setup and the first results are presented.

Joint contributions of the Pierre Auger and Telescope Array Collaborations to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013: cross-calibration of the fluorescence telescopes, large scale anisotropies and mass composition.

Recent cosmological observations hint at a deviation from the simple power-law form of the primordial spectrum of curvature perturbations. In this paper we show that in the presence of a tensor component, a turn-over in the initial spectrum is preferred by current observations, and hence non-power-law models ought to be considered. For instance, for a power-law parameterisation with both a tensor component and running parameter, current data show a preference for a negative running at more than $2.5\sigma$ C.L. As a consequence of this deviation from a power-law, constraints on the tensor-to-scalar ratio $r$ are slightly broader. We also present constraints on the inflationary parameters for a model-independent reconstruction and the Lasenby & Doran (LD) model. In particular, the constraints on the tensor-to-scalar ratio from the LD model are: $r_{\rm LD}=0.11\pm{0.024}$. In addition to current data, we show expected constraints from Planck-like and CMB-Pol sensitivity experiments by using Markov-Chain-Monte-Carlo sampling chains. For all the models, we have included the Bayesian Evidence to perform a model selection analysis. The Bayes factor, using current observations, shows a strong preference for the LD model over the standard power-law parameterisation, and provides an insight into the accuracy of differentiating models through future surveys.

The fluorescence yield data used for shower reconstruction in the Auger, HiRes and TA experiments are different, not only in the overall absolute value but also in the wavelength spectrum and the various atmospheric dependencies. The effect on the energy reconstruction of using different fluorescence yield parameterizations is discussed. In addition, the impact of a change in the fluorescence spectrum depends on the optical efficiency of the telescopes. A simple analytical procedure allows us to evaluate the combined effect of fluorescence yield and optical efficiency showing a non-negligible deviation between the energy scales of TA and Auger. However no relevant effect is found in the comparison between HiRes and Auger. Finally we show that a similar procedure could also be applied with real data.

One of the main challenges of modern cosmology is to investigate the nature of dark energy in our Universe. The properties of such a component are normally summarised as a perfect fluid with a (potentially) time-dependent equation-of-state parameter $w(z)$. We investigate the evolution of this parameter with redshift by performing a Bayesian analysis of current cosmological observations. We model the temporal evolution as piecewise linear in redshift between `nodes', whose $w$-values and redshifts are allowed to vary. The optimal number of nodes is chosen by the Bayesian evidence. In this way, we can both determine the complexity supported by current data and locate any features present in $w(z)$. We compare this node-based reconstruction with some previously well-studied parameterisations: the Chevallier-Polarski-Linder (CPL), the Jassal-Bagla-Padmanabhan (JBP) and the Felice-Nesseris-Tsujikawa (FNT). By comparing the Bayesian evidence for all of these models we find an indication towards possible time-dependence in the dark energy equation-of-state. It is also worth noting that the CPL and JBP models are strongly disfavoured, whilst the FNT is just significantly disfavoured, when compared to a simple cosmological constant $w=-1$. We find that our node-based reconstruction model is slightly disfavoured with respect to the $\Lambda$CDM model.

The preferred shape for the primordial spectrum of curvature perturbations is determined by performing a Bayesian model selection analysis of cosmological observations. We first reconstruct the spectrum modelled as piecewise linear in \log k between nodes in k-space whose amplitudes and positions are allowed to vary. The number of nodes together with their positions are chosen by the Bayesian evidence, so that we can both determine the complexity supported by the data and locate any features present in the spectrum. In addition to the node-based reconstruction, we consider a set of parameterised models for the primordial spectrum: the standard power-law parameterisation, the spectrum produced from the Lasenby & Doran (LD) model and a simple variant parameterisation. By comparing the Bayesian evidence for different classes of spectra, we find the power-law parameterisation is significantly disfavoured by current cosmological observations, which show a preference for the LD model.

We constrain the shape of the primordial power spectrum using recent measurements of the cosmic microwave background (CMB) from the Wilkinson Microwave Anisotropy Probe (WMAP) 7-year data and other high-resolution CMB experiments. We also include observations of the matter power spectrum from the luminous red galaxy (LRG) subset DR7 of the Sloan Digital Sky Survey (SDSS). We consider two different models of the primordial power spectrum. The first is the standard nearly scale-invariant spectrum in the form of a generalised power-law parameterised in terms of the spectral amplitude $A_{\rm s}$, the spectral index $n_{\rm s}$ and (possibly) the running parameter $n_{\rm run}$. The second spectrum is derived from the Lasenby and Doran (LD) model. The LD model is based on the restriction of the total conformal time available in a closed Universe and the predicted primordial power spectrum depends upon just two parameters. An important feature of the LD spectrum is that it naturally incorporates an exponential fall-off on large scales, which might provide a possible explanation for the lower-than-expected power observed at low multipoles in the CMB. In addition to parameter estimation, we compare both models using Bayesian model selection. We find there is a significant preference for the LD model over a simple power-law spectrum for a CMB-only dataset, and over models with an equal number of parameters for all the datasets considered.

An accurate knowledge of the fluorescence yield and its dependence on atmospheric properties such as pressure, temperature or humidity is essential to obtain a reliable measurement of the primary energy of cosmic rays in experiments using the fluorescence technique. In this work, several sets of fluorescence yield data (i.e. absolute value and quenching parameters) are described and compared. A simple procedure to study the effect of the assumed fluorescence yield on the reconstructed shower parameters (energy and shower maximum depth) as a function of the primary features has been developed. As an application, the effect of water vapor and temperature dependence of the collisional cross section on the fluorescence yield and its impact on the reconstruction of primary energy and shower maximum depth has been studied.

The dynamics of a cosmological model fueled by scalar field dark matter with a cosh-like potential plus a cosmological constant is investigated in detail. It is revealed that the late-time attractor is always the de Sitter solution, and that, depending on the values of the free parameters, the oscillating solution of the scalar field -- modeling cold dark matter -- mediates between some early stage (say, the radiation-dominated solution) and the accelerating de Sitter attractor.