This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

We study nonlinear structure formation in the Bound Dark Energy model (BDE), where dark energy (DE) corresponds to a light scalar meson particle $\phi$ dynamically formed at a condensation energy scale $\Lambda_c$. The evolution of this dark-energy meson is determined by the potential $V(\phi)=\Lambda_c^{4+2/3}\phi^{-2/3}$, with a distinguishing phenomenology from other quintessence scenarios. Particularly, the expansion rate of the universe is affected not only at late times, but also when the condensation of $\phi$ occurs, which in linear theory leads to an enhancement (with respect to standard $\Lambda$CDM) of matter perturbations on small scales. We study how much of this signature is still present at late times as well as the properties of dark matter halos in the nonlinear regime through N-body simulations. Our results show that nonlinear corrections wash out this feature from the matter power spectrum even before DE becomes dominant. There is, however, a small but clear suppression of the BDE spectrum of $2\%$ today on the largest scales due to the distinct late-time dynamics of DE. The differences on the clustering power between BDE and $\Lambda$CDM are reflected in the halo mass function, where small halos are more abundant in BDE as opposed to large heavy structures, whose formation is delayed because of the expansion history of the universe. This result is well captured by the semi-analytical Sheth-Tormen formula. However, despite these differences, the halo concentration parameter is essentially the same in both models, which suggest that clustering inside the halos decouple from the general expansion once the halos form.

We analyze the anisotropic clustering of the Sloan Digital Sky Survey-IV Extended Baryon Oscillation Spectroscopic Survey (eBOSS) Luminous Red Galaxy Data Release 14 (DR14) sample combined with Baryon Oscillation Spectroscopic Survey (BOSS) CMASS sample of galaxies in the redshift range 0.6$<z<$1.0, which consists of 80,118 galaxies from eBOSS and 46,439 galaxies from the BOSS-CMASS sample. The eBOSS-CMASS Luminous Red Galaxy sample has a sky coverage of 1,844 deg$^2$, with an effective volume of 0.9 Gpc$^3$. The analysis was made in configuration space using a Legendre multipole expansion. The Redshift Space Distortion signal is modeled as a combination of the Convolution Lagrangian Perturbation Model and the Gaussian Streaming Model. We constrain the logarithmic growth of structure times the amplitude of dark matter density fluctuations, $f (z_{\rm eff})\sigma_8(z_{\rm eff})=0.454 \pm0.139 $, and the Alcock-Paczynski dilation scales which constraints the angular diameter distance $D_A(z_{eff})=1466.5 \pm 136.6 (r_s/r_s^{\rm fid})$ and $H(z_{\rm eff})=105.8 \pm 16 (r_s^{\rm fid}/r_s) \mathrm{km\,s^{-1}\,Mpc^{-1}}$, where $r_s$ is the sound horizon at the end of the baryon drag epoch and $r_s^{\rm fid}$ is its value in the fiducial cosmology at an effective redshift $z_{\rm eff}=0.72$. These results are in full agreement with the current $\Lambda$-Cold Dark Matter ($\Lambda$-CDM) cosmological model inferred from Planck measurements. This study is the first eBOSS LRG full-shape analysis i.e. including Redshift-Space Distortions (RSD) simultaneously with the Alcock-Paczynski (AP) effect and the Baryon Acoustic Oscillation (BAO) scale.

Cosmological observations such as structure formation, CMB, and cosmic distance ladder set tight constraints to the amount and nature of dark matter (DM). In particular structure formation strongly constraints not only the amount of energy density, but also the time when DM became non-relativistic, anr. Standard cold and thermic warm DM particles have a smooth transition from being relativistic at high energies to a non-relativistic regime since the mass of these particles is constant while the velocity redshifts with the expansion of the universe. However, here we explore the possibility that the DM particle acquires a non-perturbative mass at a phase transition scale and a scale factor ac. This transition acquired a more fundamental meaning for the Bound Dark Matter (BDM) model because they describe a particle getting its mass through a non-perturbative process. These BDM particles may go from being relativistic to non-relativistic at ac, which implies an abrupt transition of the velocity of the particles vc at that time, affecting the equation of state and the cosmological evolution. Here we study the cosmological impact of the values of vc and the scale transition ac and they by reducing the free-streaming scale and therefore the small scale structure. Using CMB Plank, Supernovae SNIa, and BAO data, we constrain the valid region in the parameter space ac-vc putting upper bounds to ac but not restricting vc. For instance, the transition must be ac<$2.66\times10^{-6}$ for vc = 0. We also find that the free-streaming and the Jeans mass of the dark matter particle is highly influenced by the velocity vc, for example, a 3 keV WDM have a Jeans mass of Mfs=$1.97\times10^7$ Msun/h^3 but an equivalent BDM with the same anr but an abrupt transition in the velocity, vc~0 would have a Jeans mass of Mfs=$2.08x10^4$ Msun/h^3 which significantly changes the large scale structure formation.

We study the nonlinear regime of large scale structure formation considering a dynamical dark energy (DE) component determined by a Steep Equation of State parametrization (SEoS) $w(z)=w_0+w_i\frac{(z/z_T)^q}{1+(z/z_T)^q}$. In order to perform the model exploration at low computational cost, we modified the public code L-PICOLA. We incorporate the DE model by means of the first and second-order matter perturbations in the Lagrangian frame and the expansion parameter. We analyze deviations of SEoS models with respect to $\Lambda$CDM in the non-linear matter power spectrum ($P_k$), the halo mass function (HMF), and the two-point correlation function (2PCF). On quantifying the nature of steep (SEoS-I) and smooth transitions in DE field (CPL-lim), no signature of steep transition is observed, rather found the overall impact of DE behaviors in $P_k$ at level of $\sim 2-3\%$ and $\sim 3-4\%$ differences w.r.t $\Lambda$CDM at $z=0$ respectively. HMF shows the possibility to distinguish between the models at the high mass ends. The best-fitted model assuming only background and linear perturbations dubbed as SEoS-II largely deviates from $\Lambda$CDM and current observations on studying the nonlinear growth. This large deviation in SEoS-II also quantified the combined effect of the dynamical DE and the larger amount of matter contained, $\Omega_{m0}$ and $H_{0}$ accordingly. 2PCF results are relatively robust with $\sim 1-2 \%$ deviation for SEoS-I and CPL-lim and a significant deviation for SEoS-II throughout $r$ from $\Lambda$CDM. Finally, we conclude that the search for viable DE models (like the SEoS) must include non-linear growth constraints.

We present the full evolution of the velocity for a massive particle, along with the equation of state we are able to analytically compute the energy density and pressure evolution. Therefore, is easy to compute the perturbation equations for any massive decoupled particle, i.e. warm dark matter (WDM) or neutrinos, treated as fluid. Using this approach we analytical calculate the moment when the WDM stop being relativistic, a_nr, which is just 15% difference with respect to the exact Boltzmann solution. Using the fluid approximation the matter power spectrum is computed fast and with great accuracy, the cut-off in structure formation due to the free-streaming, lambda_fs, of the particle, characteristic for a WDM particle, is replicated in both matter power spectrum and halo mass function, in which we found up to 30% correction on the Jeans mass. We also show that the matter power can be computed using a fitting formula that involves only the cut-off scale, k_fs. This formulation can be integrated in comprehensive numerical modeling reasonable increasing the performance in the calculations.

We measure the imprint of primordial baryon acoustic oscillations (BAO) in the correlation function of Ly$\alpha$ absorption in quasar spectra from the Baryon Oscillation Spectroscopic Survey (BOSS) and the extended BOSS (eBOSS) in Data Release 14 (DR14) of the Sloan Digital Sky Survey (SDSS)-IV. In addition to 179,965 spectra with absorption in the Lyman-$\alpha$ (Ly$\alpha$) region, we use, for the first time, Ly$\alpha$ absorption in the Lyman-$\beta$ region of 56,154 spectra. We measure the Hubble distance, $D_H$, and the comoving angular diameter distance, $D_M$, relative to the sound horizon at the drag epoch $r_d$ at an effective redshift $z=2.34$. Using a physical model of the correlation function outside the BAO peak, we find $D_H(2.34)/r_d=8.86\pm 0.29$ and $D_M(2.34)/r_d=37.41\pm 1.86$, within 1$\sigma$ from the flat-$\Lambda$CDM model consistent with CMB anisotropy measurements. With the addition of polynomial "broadband" terms, the results remain within one standard deviation of the CMB-inspired model. Combined with the quasar-Ly$\alpha$ cross-correlation measurement presented in a companion paper Blomqvist19, the BAO measurements at $z=2.35$ are within 1.7$\sigma$ of the predictions of this model.

We present the status of the Dark Energy Spectroscopic Instrument (DESI) and its plans and opportunities for the coming decade. DESI construction and its initial five years of operations are an approved experiment of the US Department of Energy and is summarized here as context for the Astro2020 panel. Beyond 2025, DESI will require new funding to continue operations. We expect that DESI will remain one of the world's best facilities for wide-field spectroscopy throughout the decade. More about the DESI instrument and survey can be found at https://www.desi.lbl.gov.

We study the cosmological properties of a dynamical of dark energy (DE) component determined by a Steep Equation of State (SEoS) $w(z)=w_0+w_i\frac{(z/z_T)^q}{1+(z/z_T)^q}$. The SEoS has a transition at $z_T$ between two pivotal values ($w_i, w_0$) which can be taken as an early time and present day values of $w$ and the steepness is given by $q$. We describe the impact of this dynamical DE at background and perturbative level. The steepness of the transition has a better cosmological fit than a conventional CPL model with $w=w_0+w_a(1-a)$. Furthermore, we analyze the impact of steepness of the transition in the growth of matter perturbations and structure formation. This is manifest in the linear matter power spectrum, $P(k)$, the logarithmic growth function, $f\sigma_8(z)$, and the differential mass function $dn/d\log M(z=0)$. The differences in these last three quantities is at a percent-level using the same cosmological baseline parameters in our SEoS and a $\Lambda CDM$ model. However, we find an increase in the power spectrum, producing a bump at $k\approx k_T$ with $k_T\equiv a_TH(a_T)$ the mode associated to the time of the steep transition ($a_T = 1/(1+z_T)$). Different dynamics of DE lead to a different amount of DM at present time which has an impact in Power Spectrum and accordingly in structure formation.

We present constraints on local primordial non-Gaussianity (PNG), parametrized through $f^{\rm loc}_{\rm NL}$, using the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey Data Release 14 quasar sample. We measure and analyze the anisotropic clustering of the quasars in Fourier space, testing for the scale-dependent bias introduced by primordial non-Gaussianity on large scales. We derive and employ a power spectrum estimator using optimal weights that account for the redshift evolution of the PNG signal. We find constraints of $-51<f^{\rm loc}_{\rm NL}<21$ at 95% confidence level. These are amont the tightest constraints from Large Scale Structure (LSS) data. Our redshift weighting improves the error bar by 15% in comparison to the unweighted case. If quasars have lower response to PNG, the constraint degrades to $-81<f^{\rm loc}_{\rm NL}<26$, with a 40% improvement over the standard approach. We forecast that the full eBOSS dataset could reach $\sigma_{f^{\rm loc}_{\rm NL}}\simeq 5\text{-}8$ using optimal methods and full range of scales.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

We present a complete analysis of the observational constraints and cosmological implications of our Bound Dark Energy (BDE) model aimed to explain the late-time cosmic acceleration of the universe. BDE is derived from particle physics and corresponds to the lightest meson field $\phi$ dynamically formed at low energies due to the strong gauge coupling constant. The evolution of the dark energy is determined by the scalar potential $V(\phi)=\Lambda_c^{4+2/3}\phi^{-2/3}$ arising from non-perturbative effects at a condensation scale $\Lambda_c$ and scale factor $a_c$, related each other by $a_c\Lambda_c/\mathrm{eV}=1.0934\times 10^{-4}$. We present the full background and perturbation evolution at a linear level. Using current observational data, we obtain the constraints $a_c=(2.48 \pm 0.02)\times10^{-6}$ and $\Lambda_c=(44.09 \pm 0.28) \textrm{ eV}$, which is in complete agreement with our theoretical prediction $\Lambda_c^{th}=34^{+16}_{-11}\textrm{ eV}$. The bounds on the equation of state today, the dark energy density and the expansion rate are $w_\mathrm{BDE 0}=-0.929\pm 0.007$, $\Omega_\mathrm{BDE0}=0.696\pm0.007$ and $H_0=67.82\pm 0.05$ km s$^{-1}$Mpc, respectively. Even though the constraints on the six Planck base parameters are consistent at the 1$\sigma$ level between BDE and the concordance $\Lambda$CDM model, BDE improves the likelihood ratio by 2.1 of the Baryon Acoustic Oscillations (BAO) measurements with respect to $\Lambda$CDM and has an equivalent fit for type Ia supernovae and the Cosmic Microwave Background data. We present the constraints on the different cosmological parameters, and particularly we show the tension between BDE and $\Lambda$CDM in the BAO distance ratio $r_\mathrm{BAO}$ vs $H_\mathrm{0}$ and the growth index $\gamma$ at different redshifts, as well as the dark matter density at present time $\Omega_ch^2$ vs $H_0$.

We develop a new method, which is based on the optimal redshift weighting scheme, to extract the maximal tomographic information of baryonic acoustic oscillations (BAO) and redshift space distortions (RSD) from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 14 quasar (DR14Q) survey. We validate our method using the EZ mocks, and apply our pipeline to the eBOSS DR14Q sample in the redshift range of $0.8<z<2.2$. We report a joint measurement of $f\sigma_8$ and two-dimensional BAO parameters $D_{\rm A}$ and $H$ at four effective redshifts of $z_{\rm eff}=0.98, 1.23, 1.52$ and $1.94$, and provide the full data covariance matrix. Using our measurement combined with BOSS DR12, MGS and 6dFGS BAO measurements, we find that the existence of dark energy is supported by observations at a $7.4\sigma$ significance level. Combining our measurement with BOSS DR12 and Planck observations, we constrain the gravitational growth index to be $\gamma=0.580\pm0.082$, which is fully consistent with the prediction of general relativity. This paper is part of a set that analyses the eBOSS DR14 quasar sample.

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.

We present a Dark Energy (DE) model with a sound derivation as a natural extension of the Standard Model of particle physics with no free parameters and an excellent fit with current cosmological data improving by 21% the $\Lambda$CDM fit of the Baryon Acoustic Oscillations (BAO) measurements, specially designed to determine the dynamics of DE. DE corresponds to the lightest bound state scalar particle $\phi$ with a potential $V=\Lambda_c^{4+2/3}\phi^{-2/3}$ dynamically formed at the condensation energy scale $\Lambda_c$ and scale factor $a_c$. The value of $\Lambda_c$, the exponent $n=2/3$, and the initial conditions of $\phi$ are all derived quantities. We obtain an exact constraint $a_c\Lambda_c/\textrm{eV}=1.0939\times 10^{-4}$ and a theoretical prediction $\Lambda_c=34 ^{+16}_{-11} \textrm{ eV}$, consistent with the best fit $\Lambda_c=44.08\pm 0.27 \textrm{ eV}$. We test our model constraint on $a_c\Lambda_c$ by allowing $a_c$ and $\Lambda_c$ to vary independently and remarkably our prediction has a relative difference of only 0.2% with the best fit value. Unlike a cosmological constant $\Lambda$, our DE model predicts the amount of DE and leaves detectable cosmological imprints at different times and scales at a background and perturbation level.

We present a measurement of the anisotropic and isotropic Baryon Acoustic Oscillations (BAO) from the extended Baryon Oscillation Spectroscopic Survey Data Release 14 quasar sample with optimal redshift weights. Applying the redshift weights improves the constraint on the BAO dilation parameter $\alpha(z_{\rm eff})$ by 17\%. We reconstruct the evolution history of the BAO distance indicators in the redshift range of $0.8<z<2.2$. This paper is part of a set that analyses the eBOSS DR14 quasar sample.

We present an anisotropic analysis of Baryon Acoustic Oscillation (BAO) signal from the SDSS-IV extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 14 (DR14) quasar sample. The sample consists of 147,000 quasars distributed over a redshift range of $0.8 < z < 2.2$. We apply the redshift weights technique to the clustering of quasars in this sample and achieve a 4.6 per cent measurement of the angular distance measurement $D_M$ at $z = 2.2$ and Hubble parameter $H$ at $z=0.8$. We parameterize the distance-redshift relation, relative to a fiducial model, as a quadratic expansion. The coefficients of this expansion are used to reconstruct the distance-redshift relation and obtain distance and Hubble parameter measurements at all redshifts within the redshift range of the sample. Reporting the result at two characteristic redshifts, we determine $D_M(z=1) = 3405\pm305 \ (r_{\rm d} / r_{\rm d, fid}) \ {\rm Mpc}$, $H(z=1) = 120.7\pm 7.3 \ (r_{\rm d,fid} / r_{\rm d}) \ {\rm km} \ {\rm s}^{-1}{\rm Mpc}^{-1}$ and $D_M(z=2) = 5325\pm249 \ (r_{\rm d} / r_{\rm d, fid}) \ {\rm Mpc}$, $H(z=2) = 189.9\pm 32.9 \ (r_{\rm d,fid} / r_{\rm d}) \ {\rm km} \ {\rm s}^{-1}{\rm Mpc}^{-1}$. These measurements are highly correlated. We assess the outlook of BAO analysis from the final quasar sample by testing the method on a set of mocks that mimic the noise level in the final sample. We demonstrate on these mocks that redshift weighting shrinks the measurement error by over 25 per cent on average. We conclude redshift weighting can bring us closer to the cosmological goal of the final quasar sample.

We present measurements of the Baryon Acoustic Oscillation (BAO) scale in redshift-space using the clustering of quasars. We consider a sample of 147,000 quasars from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) distributed over 2044 square degrees with redshifts $0.8 < z < 2.2$ and measure their spherically-averaged clustering in both configuration and Fourier space. Our observational dataset and the 1400 simulated realizations of the dataset allow us to detect a preference for BAO that is greater than 2.8$\sigma$. We determine the spherically averaged BAO distance to $z = 1.52$ to 3.8 per cent precision: $D_V(z=1.52)=3843\pm147 \left(r_{\rm d}/r_{\rm d, fid}\right)\ $Mpc. This is the first time the location of the BAO feature has been measured between redshifts 1 and 2. Our result is fully consistent with the prediction obtained by extrapolating the Planck flat $\Lambda$CDM best-fit cosmology. All of our results are consistent with basic large-scale structure (LSS) theory, confirming quasars to be a reliable tracer of LSS, and provide a starting point for numerous cosmological tests to be performed with eBOSS quasar samples. We combine our result with previous, independent, BAO distance measurements to construct an updated BAO distance-ladder. Using these BAO data alone and marginalizing over the length of the standard ruler, we find $\Omega_{\Lambda} > 0$ at 6.6$\sigma$ significance when testing a $\Lambda$CDM model with free curvature.

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.

We present a parametrization for the Dark Energy Equation of State "EoS" which has a rich structure, performing a transition at pivotal redshift $z_T$ between the present day value $w_0$ to an early time $w_i=w_a+w_0\equiv w(z\gg0)$ with a steepness given in terms of $q$ parameter. The proposed parametrization is $w=w_0+w_a(z/z_T)^q/(1+(z/z_T))^q$, with $w_0$, $w_i$, $q $ and $z_T$ constant parameters. It reduces to the widely used EoS $w=w_0+w_a(1-a)$ for $z_T=q=1$. This transition is motivated by scalar field dynamics such as for example quintessence models. We study if a late time transition is favored by BAO measurements combined with local determination of $H_0$ and information from the CMB. According to our results, an EoS with a present value of $w_0 = -0.92$ and a high redshift value $w_i =-0.99$, featuring a transition at $z_T = 0.28$ with an exponent $q = 9.97$ was favored by data coming from local dynamics of the Universe (BAO combined with $H_0$ determination). We find that a dynamical DE model allows to simultaneously fit $H_0$ from local determinations and Planck CMB measurements, alleviating the tension obtained in a $\Lambda$CDM model. Additionally to this analysis we solved numerically the evolution of matter over-densities in the presence of dark energy both at background level and when its perturbations were considered. We show that the presence of a steep transition in the DE EoS gets imprinted into the evolution of matter overdensities and that the addition of an effective sound speed term does not erase such feature.

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.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

We present a parametrization for the Dark Energy Equation of State "EoS" which has a rich structure. Our EoS has a transition at pivotal redshift $z_T$ between the present day value $w_0$ to an early time $w_i=w_a+w_0\equiv w(z>>0)$ and the steepness of this transition is given in terms of the $q$ parameter. The proposed parametrization is $w=w_0+w_a(z/z_T)^q/(1+(z/z_T))^q$, with $w_0$, $w_i$, $q $ and $z_T$ constant parameters. This transition is motivated by scalar field dynamics such as for example quintessence models. Our parametrization reduces to the widely used EoS $w=w_0+w_a(1-a)$ for $z_T=q=1$. We study if a late time transition is favored by BAO measurements and Planck priors. According to our results, an EoS with a present value of $w_0 = -0.91$ and a high redshifts value $w_i =-0.62$, featuring a transition at a redshift of $z_T = 1.16$ with an exponent $q = 9.95$ is a good fit to the observational data. We found good agreement between the model and the data reported by the different surveys. A "thawing" dynamics is preferred by the use of BAO data alone (including Lymman-$\alpha$ forest measurements) and a "freezing" evolution of the EoS is preferred when we include the priors from Planck. The constraints imposed by the available BAO measurements (\cite{Beutler:2011hx, Ross:2014qpa, Anderson:2013oza, Kazin:2014qga, Font-Ribera:2013wce, Delubac:2014aqe, Gong:2015tta}) and its physical behavior are discussed.

We present a new selection technique of producing spectroscopic target catalogues for massive spectroscopic surveys for cosmology. This work was conducted in the context of the extended Baryon Oscillation Spectroscopic Survey (eBOSS), which will use ~200 000 emission line galaxies (ELGs) at 0.6<zspec<1.0 to obtain a precise baryon acoustic oscillation measurement. Our proposed selection technique is based on optical and near-infrared broad-band filter photometry. We used a training sample to define a quantity, the Fisher discriminant (linear combination of colours), which correlates best with the desired properties of the target: redshift and [OII] flux. The proposed selections are simply done by applying a cut on magnitudes and this Fisher discriminant. We used public data and dedicated SDSS spectroscopy to quantify the redshift distribution and [OII] flux of our ELG target selections. We demonstrate that two of our selections fulfil the initial eBOSS/ELG redshift requirements: for a target density of 180 deg^2, ~70% of the selected objects have 0.6<zspec<1.0 and only ~1% of those galaxies in the range 0.6<zspec<1.0 are expected to have a catastrophic zspec estimate. Additionally, the stacked spectra and stacked deep images for those two selections show characteristic features of star-forming galaxies. The proposed approach using the Fisher discriminant could, however, be used to efficiently select other galaxy populations, based on multi-band photometry, providing that spectroscopic information is available. This technique could thus be useful for other future massive spectroscopic surveys such as PFS, DESI, and 4MOST.

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct novel cosmological observations using the BOSS spectrograph at Apache Point Observatory. Observations will be simultaneous with the Time Domain Spectroscopic Survey (TDSS) designed for variability studies and the Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for studies of X-ray sources. eBOSS will use four different tracers to measure the distance-redshift relation with baryon acoustic oscillations (BAO). Using more than 250,000 new, spectroscopically confirmed luminous red galaxies at a median redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000 spectroscopically-confirmed quasars will provide the first BAO distance measurements over the redshift range 0.9<z<2.2, with expected precision of 2.8% and 4.2% on $d_A(z)$ and H(z), respectively. Finally, with 60,000 new quasars and re-observation of 60,000 quasars known from BOSS, we will obtain new Lyman-alpha forest measurements at redshifts z>2.1; these new data will enhance the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS. Furthermore, eBOSS will provide improved tests of General Relativity on cosmological scales through redshift-space distortion measurements, improved tests for non-Gaussianity in the primordial density field, and new constraints on the summed mass of all neutrino species. Here, we provide an overview of the cosmological goals, spectroscopic target sample, demonstration of spectral quality from early data, and projected cosmological constraints from eBOSS.

We study the cosmological implications of the Nambu-Jona-Lasinio (NJL model) when the coupling constant is field dependent. The NJL model has a four-fermion interaction describing two different phases due to quantum interaction effects and determined by the strength of the coupling constant g. It describes massless fermions for weak coupling and a massive fermions and strong coupling, where a fermion condensate is formed. In the original NJL model the coupling constant g is indeed constant, and in this work we consider a modified version of the NJL model by introducing a dynamical field dependent coupling motivated by string theory. The effective potential as a function of the varying coupling (aimed to implement a natural phase transition) is seen to develop a negative divergence, i.e. becomes a "bottomless well" in certain limit region. Although we explain how an lower unbounded potential is not necessarily unacceptable in a cosmological context, the divergence can be removed if we consider a mass term for the coupling-like field. We found that for a proper set of parameters, the total potential obtained has two minima, one located at the origin (the trivial solution, in which the fluid associated with the fields behave like matter); and the other related to the non-trivial solution. This last solution has three possibilities: 1) if the minimum is positive V_{min}>0, the system behave as a cosmological constant, thus leading eventually to an accelerated universe; 2) if the minimized potential vanishes V_{min}=0, then we have matter with no acceleration ; 3) finally a negative minimum V_{min}<0 leads an eventually collapsing universe, even though we have a flat geometry.Therefore, a possible interpretation as Dark Matter or Dark Energy is allowed among the behaviors implicated in the model.

We review the Nambu-Jona-Lasinio model (NJL), proposed long time ago, as a four-fermion interaction theory with chiral symmetry. The theory is not renormalizable and presents a symmetry breaking due to quantum corrections which depends on the strength of the coupling constant. We may associate a phase transition with this symmetry breaking, leading from a fermion massless states to fermion condensates. This condensates can be described effectively by a scalar field. We are interested in this paper in the cosmological dynamics of the NJL model, and in the possibility that it can be related to dark energy and/or dark matter, which form up to 95% of the energy content of the universe at present time. We consider exclusively gravitational interaction between the NJL and the SM particles.

We propose a new Dark Energy parametrization based on the dynamics of a scalar field. We use an equation of state w=(x-1)/(x+1), with x=E_k/V, the ratio of kinetic energy E_k=\dotphi^2/2 and potential V. The equation of motion gives x=(L/6)(V/3H^2) and has a solution x=([(1+y)^2+2 L/3]^{1/2}-(1+y))/2 where y\equiv \rmm/V and L= (V'/V)^2 (1+q)^2, q=\ddotphi/V'. The resulting EoS is w=[6+ L- 6 \sqrt((1+y)^2+2L/3)]/(L+6y). Since the universe is accelerating at present time we use the slow roll approximation in which case we have |q|<< 1 and L\simeq (V'/V)^2. However, the derivation of w is exact and has no approximation. By choosing an appropriate ansatz for L we obtain a wide class of behavior for the evolution of Dark Energy without the need to specify the potential V. The EoS w can either grow and later decrease, or other way around, as a function of redshift and it is constraint between -1\leq w\leq 1 as for any canonical scalar field with only gravitational interaction. To determine the dynamics of Dark Energy we calculate the background evolution and its perturbations, since they are important to discriminate between different DE models. Our parametrization follows closely the dynamics of a scalar field scalar fields and the function L allow us to connect it with the potential V(phi) of the scalar field phi.

It has been argued that the small perturbations to the homogeneous and isotropic configurations of a canonical scalar field in an expanding universe do not grow. We show that this is not true in general, and clarify the root of the misunderstanding. We revisit a simple model in which the zero-mode of a free scalar field oscillates with high frequency around the minimum of the potential. Under this assumption the linear perturbations grow like those in the standard cold dark matter scenario, but with a Jeans length at the scale of the Compton wavelength of the scalar particle. Contrary to previous analyses in the literature our results do not rely on time-averages and/or fluid identifications, and instead we solve both analytically (in terms of a well-defined series expansion) and numerically the linearized Einstein-Klein-Gordon system. Also, we use gauge-invariant fields, which makes the physical analysis more transparent and simplifies the comparison with previous works carried out in different gauges. As a byproduct of this study we identify a time-dependent modulation of the different physical quantities associated to the background as well as the perturbations with potential observational consequences in dark matter models.

Large sets of objects with spectroscopic redshift measurements will be needed for imaging dark energy experiments to achieve their full potential, serving two goals:_training_, i.e., the use of objects with known redshift to develop and optimize photometric redshift algorithms; and_calibration_, i.e., the characterization of moments of redshift (or photo-z error) distributions. Better training makes cosmological constraints from a given experiment stronger, while highly-accurate calibration is needed for photo-z systematics not to dominate errors. In this white paper, we investigate the required scope of spectroscopic datasets which can serve both these purposes for ongoing and next-generation dark energy experiments, as well as the time required to obtain such data with instruments available in the next decade. Large time allocations on kilo-object spectrographs will be necessary, ideally augmented by infrared spectroscopy from space. Alternatively, precision calibrations could be obtained by measuring cross-correlation statistics using samples of bright objects from a large baryon acoustic oscillation experiment such as DESI. We also summarize the additional work on photometric redshift methods needed to prepare for ongoing and future dark energy experiments.

In order to unravel the nature of the dark matter (DM) we have proposed a particle-physics motivated model called Bound Dark Matter (BDM) that consist in DM massless particles above a threshold energy Ec that acquire mass below it due to nonperturbative methods. Therefore, the BDM model describes DM particles which are relativistic, hot dark matter (HDM) in the denser (inner) regions of galaxies and describes nonrelativistic, cold dark matter (CDM) where halo density is below rho_c = Ec^4. We test this model by fitting rotation curves from Low Surface Brightness (LSB) galaxies from The HI Nearby Galaxy Survey (THINGS). We use a particular DM cored profile that contains three parameters: a typical scale length (rs) and density (rho_0) of the halo, and a core radius (rc) stemming from the relativistic nature of the BDM model. Since the energy Ec parameterizes the phase transition due to the underlying particle physics model, it is independent on the details of galaxy and/or structure formation and therefore the DM profile parameters rs, rc, Ec are constrained, leaving only two free parameters. Through the results we agree with previous ones implying that cored profiles are preferred over the N-body motivated cuspy profiles. We also compute 2D likelihoods of the BDM parameters rc and Ec for the different galaxies and matter contents, and find an average galaxy core radius rc = 1.48kpc and a transition energy between hot and cold dark matter at Ec = 0.06 eV. The phase transition scale Ec is a new fundamental scale for our DM model well motivated theoretical origin from gauge group dynamics.

Recently a new -quantum motivated- theory of gravity has been proposed that modifies the standard Newtonian potential at large distances when spherical symmetry is considered. Accordingly, Newtonian gravity is altered by adding an extra Rindler acceleration term that has to be phenomenologically determined. Here we consider a standard and a power-law generalization of the Rindler modified Newtonian potential. The new terms in the gravitational potential are hypothesized to play the role of dark matter in galaxies. Our galactic model includes the mass of the integrated gas, and stars for which we consider three stellar mass functions (Kroupa, diet-Salpeter, and free mass model). We test this idea by fitting rotation curves of seventeen low surface brightness galaxies from The HI Nearby Galaxy Survey (THINGS). We find that the Rindler parameters do not perform a suitable fit to the rotation curves in comparison to standard dark matter profiles (Navarro-Frenk-White and Burkert) and, in addition, the computed parameters of the Rindler gravity show a high spread, posing the model as a nonacceptable alternative to dark matter.

This document summarizes the results of a community-based discussion of the potential science impact of the Mayall+BigBOSS highly multiplexed multi-object spectroscopic capability. The KPNO Mayall 4m telescope equipped with the DOE- and internationally-funded BigBOSS spectrograph offers one of the most cost-efficient ways of accomplishing many of the pressing scientific goals identified for this decade by the "New Worlds, New Horizons" report. The BigBOSS Key Project will place unprecedented constraints on cosmological parameters related to the expansion history of the universe. With the addition of an open (publicly funded) community access component, the scientific impact of BigBOSS can be extended to many important astrophysical questions related to the origin and evolution of galaxies, stars, and the IGM. Massive spectroscopy is the critical missing ingredient in numerous ongoing and planned ground- and space-based surveys, and BigBOSS is unique in its ability to provide this to the US community. BigBOSS data from community-led projects will play a vital role in the education and training of students and in maintaining US leadership in these fields of astrophysics. We urge the NSF-AST division to support community science with the BigBOSS multi-object spectrograph through the period of the BigBOSS survey in order to ensure public access to the extraordinary spectroscopic capability.

Recently a new particle physics model called Bound Dark Matter (BDM) has been proposed in which dark matter (DM) particles are massless above a threshold energy (Ec) and acquire mass below it due to nonperturbative methods. Therefore, the BDM model describes DM particles which are relativistic, hot dark matter, in the inner regions of galaxies and describes nonrelativistic, cold dark matter, where halo density is below rho_c = Ec^4. To realize this idea in galaxies we use a particular DM cored profile that contains three parameters: a scale length (rs) and density (rho_0) of the halo, and a core radius (rc) stemming from the relativistic nature of the BDM model. We test this model by fitting rotation curves of seventeen Low Surface Brightness galaxies from The HI Nearby Galaxy Survey (THINGS). Since the energy Ec parameterizes the phase transition due to the underlying particle physics model, it is independent on the details of galaxy or structure formation and therefore the DM profile parameters rs, rc, Ec are constrained, leaving only two free parameters. The high spatial and velocity resolution of this sample allows to derive the model parameters through the numerical implementation of the chi^2-goodness-of-fit test to the mass models. We compare the fittings with those of Navarro-Frenk-White (NFW), Burkert, and Pseudo-Isothermal (ISO) profiles. Through the results we conclude that the BDM profile fits better, or equally well, than NFW, Burkert, and ISO profiles and agree with previous results implying that cored profiles are preferred over the N-body motivated cuspy profile. We also compute 2D likelihoods of the BDM parameters rc and Ec for the different galaxies and matter contents, and find an average galaxy core radius rc=300 pc and a transition energy Ec = 0.11 eV when the DM halo is the only component. In Kroupa mass model, we obtain a core rc=1.48 kpc, and energy Ec=0.06 eV.

We present a model where inflation and Dark Matter takes place via a single scalar field phi. Without introducing any new parameters we are able unify inflation and Dark Matter using a scalar field phi that accounts for inflation at an early epoch while it gives a Dark Matter WIMP particle at low energies. After inflation our universe must be reheated and we must have a long period of radiation dominated before the epoch of Dark Matter. Typically the inflaton decays while it oscillates around the minimum of its potential. If the inflaton decay is not complete or sufficient then the remaining energy density of the inflaton after reheating must be fine tuned to give the correct amount of Dark Matter. An essential feature here, is that Dark Matter-Inflaton particle is produced at low energies without fine tuning or new parameters. This process uses the same coupling g as for the inflaton decay. Once the field phi becomes non-relativistic it will decouple as any WIMP particle, since n_phi is exponentially suppressed. The correct amount of Dark Matter determines the cross section and we have a constraint between the coupling $g$ and the mass $m_o$ of phi. The unification scheme we present here has four free parameters, two for the scalar potential V(phi) given by the inflation parameter lambda of the quartic term and the mass m_o. The other two parameters are the coupling $g$ between the inflaton phi and a scalar filed varphi and the coupling h between varphi with standard model particles psi or chi. These four parameters are already present in models of inflation and reheating process, without considering Dark Matter. Therefore, our unification scheme does not increase the number of parameters and it accomplishes the desired unification between the inflaton and Dark Matter for free.

We analyze the a set of seventeen rotation curves of Low Surface Brightness (LSB) galaxies from the The HI Nearby Galaxy Survey (THINGS) with different mass models to study the core structure and to determine a phase transition energy scale (E_c) between hot and cold dark matter, due to nonperturbative effects in the Bound Dark Matter (BDM) model. Our results agree with previous ones implying the cored profiles are preferred over the N-body motivated cuspy NFW profile. We find an average galactic core radius of r_c = 260 \times 10^{+/- 1.3} pc and a phase transition energy E_c = 0.11\times 10^{+/- 0.46} eV, that is of the same order of magnitude as the sum of the neutrino masses.

We propose a new Dark Energy parametrization based on the dynamics of a scalar field. We use an equation of state $w=(x-1)/(x+1)$, with $x=E_k/V$, the ratio of kinetic energy $E_k=\dot\phi^2/2$ and potential $V$. The eq. of motion gives $x=(L/6)(V/3H^2)$ and with a solution $x=([1+2 L/3(1+y)]^{1/2}-1)(1+y)/2$ where $y\equiv \rm/V$ and $L\equiv (V'/V)^2 (1+q)^2,\, q\equiv\ddot\p/V'$. Since the universe is accelerating at present time we use the slow roll approximation in which case we have $|q|\ll 1$ and $L\simeq (V'/V)^2$. However, the derivation of $L$ is exact and has no approximation. By choosing an appropriate ansatz for $L$ we obtain a wide class of behavior for the evolution of Dark Energy without the need to specify the potential $V$. In fact $w$ can either grow and later decrease, or other way around, as a function of redshift and it is constraint between $-1\leq w\leq 1$ as for any canonical scalar field with only gravitational interaction. Furthermore, we also calculate the perturbations of DE and since the evolution of DE is motivated by the dynamics of a scalar field the homogenous and its perturbations can be used to determine the form of the potential and the nature of Dark Energy. Since our parametrization is on $L$ we can easily connect it with the scalar potential $V(\phi)$.

BigBOSS is a Stage IV ground-based dark energy experiment to study baryon acoustic oscillations (BAO) and the growth of structure with a wide-area galaxy and quasar redshift survey over 14,000 square degrees. It has been conditionally accepted by NOAO in response to a call for major new instrumentation and a high-impact science program for the 4-m Mayall telescope at Kitt Peak. The BigBOSS instrument is a robotically-actuated, fiber-fed spectrograph capable of taking 5000 simultaneous spectra over a wavelength range from 340 nm to 1060 nm, with a resolution R = 3000-4800. Using data from imaging surveys that are already underway, spectroscopic targets are selected that trace the underlying dark matter distribution. In particular, targets include luminous red galaxies (LRGs) up to z = 1.0, extending the BOSS LRG survey in both redshift and survey area. To probe the universe out to even higher redshift, BigBOSS will target bright [OII] emission line galaxies (ELGs) up to z = 1.7. In total, 20 million galaxy redshifts are obtained to measure the BAO feature, trace the matter power spectrum at smaller scales, and detect redshift space distortions. BigBOSS will provide additional constraints on early dark energy and on the curvature of the universe by measuring the Ly-alpha forest in the spectra of over 600,000 2.2 < z < 3.5 quasars. BigBOSS galaxy BAO measurements combined with an analysis of the broadband power, including the Ly-alpha forest in BigBOSS quasar spectra, achieves a FOM of 395 with Planck plus Stage III priors. This FOM is based on conservative assumptions for the analysis of broad band power (kmax = 0.15), and could grow to over 600 if current work allows us to push the analysis to higher wave numbers (kmax = 0.3). BigBOSS will also place constraints on theories of modified gravity and inflation, and will measure the sum of neutrino masses to 0.024 eV accuracy.

We unify inflation and dark matter via a single scalar field phi. One of the main difficulties for this unification is that between inflation and dark matter one needs a successful reheating process and a long lasting period of radiation. Therefore the amount of energy density in the inflaton-dark matter field phi must be fine tune after reheating to give dark matter. Here we show an alternative scheme in which the inflaton decays completely, disappearing entirely from the spectrum. However, at low energies, before matter-radiation equality, the same interaction term that leads to the inflaton decay, regenerates phi. An essential feature is that the transition between the intermediate radiation dominated to the dark matter phase is related to a quantum generation of the scalar field \phi instead to purely classical dynamics. Thanks to this quantum transition the inflation-dark matter unification can take place naturally without fine tuning. The unification scheme presented here has three parameters, the mass of the dark matter particle m_o, the inflation parameter lambda and the coupling g for the inflaton interaction. Phenomenology fixes the value for lambda and gives a constraint between g and m_o, leaving only the mass of the dark matter particle m_o as a free parameter. These same three parameters lambda, m_o, g are present in models with inflation and a dark matter wimp particle but without unification. Therefore our unification scheme does not increase the number of parameters and it accomplishes the desired unification between inflaton and dark matter for free.

We present a type of dark energy models where the particles of dark energy phi are dynamically produced via a quantum transition at very low energies. The scale where the transition takes places depends on the strength g of the interaction between phi and a relativistic field varphi. We show that a g \simeq 10^{-12} gives a generation scale E_gen simeq eV with a cross section sigma simeq 1 pb close to the WIMPs cross section sigma_w simeq pb at decoupling. The number density n_phi of the \phi particles are a source term in the eq. of motion of phi and it generates the scalar potential v(phi) responsible for the late time acceleration of our universe. Since the appearance of phi may be at very low scales, close to present time, the cosmological coincidence problem can be explained simply due to the size of the coupling constant. We unify dark energy with inflation in terms of a single scalar field phi, and we use the same potential v for inflation and dark energy. However, after inflation phi decays completely and reheats the universe at a scale E_RH propto h^2 m_Pl, where h is the coupling between the SM particles and varphi. The field phi disappears from the spectrum during, from reheating until its re-generation, and therefore it does not interfere with the standard decelerating radiation/matter cosmological model allowing for a successful unification scheme. The same interaction term that gives rise to the inflaton decay accounts for the late time re-generation of phi. We present a simple model where the strength of the g and h couplings are set by the inflation scale E_I with g=h^2 propto E_I/m_Pl giving a reheating scale E_RH propto E_I=100TeV and phi-generation scale E_gen propto E_I^2/m_pl=eV << E_RH.

We present a new dark matter model BDM which is an hybrid between hot dark matter HDM and cold dark matter CDM, in which the BDM particles behave as HDM above the energy scale E_c and as CDM below this scale. Evolution of structure formation is similar to that of CDM model but BDM predicts a nonvanishing free streaming l_fs scale and a inner galaxy core radius r_core, both quantities determined in terms of a single parameter E_c, which corresponds to the phase transition energy scale of the subjacent elementary particle model. For energies above E_c or for a scale factor a smaller then a_c, with a<a_c<a_{eq}, the particles are massless and rho redshifts as radiation. However, once the energy becomes E\leq E_c or a>a_c then the BDM particles acquire a large mass through a non perturbative mechanism, as baryons do, and rho redshifts as matter with the particles having a vanishing velocity. Typical energies are E_c=O(10-100) eV giving a l_fs \propto E_c^{-4/3}\lesssim Mpc and m_fs\propto E_c^{-4}\lesssim 10^9 M\odot. A l_fs\neq 0, r_core\neq 0 help to resolve some of the shortcomings of CDM such as overabundance substructure in CDM halos and numerical fit to rotation curves in dwarf spheroidal and LSB galaxies. Finally, our BDM model and the phase transition scale E_c can be derived from particle physics.

We study a scalar phi field that unifies inflation and dark energy with a long period of a hot decelerating universe in between these two stages of inflation. A key feature is that the transition between the intermediate decelerated phase to the dark energy phase is related to a quantum regeneration of the scalar field phi instead to purely classical dynamics. The interaction V_int(phi,varphi) between phi and a second scalar field varphi allows not only for phi to decay at a high energy E_I, even though phi does not oscillate around the minimum of its potential, but it also regenerates phi at an energy E_{BD} \simeq E_I^2 close to present time. It predicts extra relativistic energy and an interacting dark energy which may account for a w<-1. In the model presented here we have E_I \simeq 100 TeV and E_{BD} \simeq 1 eV. The low value of the back decay energy E_{BD} gives an explanation of the coincidence problem.

We study cosmological solutions for the very early universe beginning at the Planck scale for a universe containing radiation, curvature and, as a simplification of a possible scalar field potential, a cosmological constant term. The solutions are the natural counterpart of the well known results for a post-inflationary universe of non-relativistic matter, curvature and a cosmological constant. Contrary to the common belief that inflation arises independently of the initial curvature we show that in the positive curvature case the universe collapses again into a Big Crunch without allowing the cosmological term to dominate and to produce inflation. There is a critical value for the cosmological constant which divide the regions where inflation is allowed from those where inflation cannot occur. One can also have loitering solutions where the scale factor remains almost constant growing to produce inflation (or decreasing to a Big Crunch) after a time which depends on the amount of energy above (or below) the critical energy. At the critical energy the solution approaches asymptotically a particular value for the scale factor (Einstein's static pre-inflationary universe). The cases where the cosmological term vanishes or becomes negative are also studied providing a complete discussion of Friedmann models.

We consider the spectrum of mesons for the gauge theory dual to a supergravity configuration of intersecting D3/D7 branes, and use the expression for the Lagrangian of the scalar mesons to compute explicitly the Lagrangian for the lightest states in the infrared limit. Assuming that the matter content of this gauge theory is part of a hidden sector, which interacts with the standard model only via gravity, we explore the cosmological consequences of these lightest scalar mesons for a FRW universe. We show that phantom fields may appear naturally in this kind of scenarios.

We study the cosmological evolution of two coupled scalar fields with an arbitrary interaction term V_T(\phi,\vp) in the presence of a barotropic fluid, which can be matter or radiation. The force between the barotropic fluid and the scalar fields is only gravitational. We show that the dynamics is completely determine by only three parameters L_i i=1,2,3. We determine all critical points and study their stability. We find six different attractor solutions depending on the values of L_i and we calculate the relevant cosmological parameters. We discuss the possibility of having one of the scalar fields as of dark energy while the other could be a scalar field redshifting as matter.

A dark energy component is responsible for the present stage of acceleration of our universe. If no fine tuning is assumed on the dark energy potential then it will end up dominating the universe at late times and the universe will not stop this stage of acceleration. On the other hand, the equation of state of dark energy seems to be smaller than -1 as suggested by the cosmological data. We take this as an indication that dark energy does indeed interact with another fluid (we consider fermion fields) and we determine the interaction through the cosmological data and extrapolate it into the future. We study the conditions under which a dark energy can dilute faster or decay into the fermion fields. We show that it is possible to live now in an accelerating epoch dominated by the dark energy and without introducing any fine tuning parameters the dark energy can either dilute faster or decaying into fermions in the future. The acceleration of the universe will then cease.

We consider the impact of the Heidelberg-Moscow claim for a detection of neutrino mass on the determination of the dark energy equation of state. By combining the Heidelberg-Moscow result with the WMAP 3-years data and other cosmological datasets we constrain the equation of state to -1.67< w <-1.05 at 95% c.l., While future data are certainly needed for a confirmation of the controversial Heildelberg-Moscow claim, our result shows that future laboratory searches for neutrino masses may play a crucial role in the determination of the dark energy properties.

We study the possibility that dark energy decays in the future and the universe stops accelerating. The fact thatthe cosmological observations prefer an equation of state of dark energy smaller than -1 can be a signal that dark energy will decay in the future. This conclusion is based in interpreting a w<-1 as a signal of dark energy interaction with another fluid. We determine the interaction through the cosmological data and extrapolate it into the future. The resulting energy density for dark energy becomes rho=a^{-3(1+w)}e^{-\beta(a-1)}, i.e. it has an exponential suppression for a >> a_o=1. In this scenario the universe ends up dominated by this other fluid, which could be matter, and the universe stops accelerating at some time in the near future.

Recently, the Heidelberg-Moscow double beta decay experiment has claimed a detection for a neutrino mass with high significance. Here we consider the impact of this measurement on the determination of the dark energy equation of state. By combining the Heidelberg-Moscow result with the WMAP 3-years data and other cosmological datasets we constrain the equation of state to -1.67< w <-1.05 at 95% c.l., ruling out a cosmological constant at more than 95% c.l.. Interestingly enough, coupled neutrino-dark energy models may be consistent with such equation of state. While future data are certainly needed for a confirmation of the controversial Heildelberg-Moscow claim, our result shows that future laboratory searches for neutrino masses may play a crucial role in the determination of the dark energy properties.

We construct a two-stage inflationary model which can accommodate early inflation at a scale $\Lambda_1$ as well as a second stage of inflation at $\Lambda_2$ with a single scalar field $\phi$. We use a symmetric potential, valid in a frictionless world, in which the two inflationary periods have exactly the same scale, i.e. $\Lambda_1=\Lambda_2$. However, we see today $\Lambda_1 \gg \Lambda_2$ due to the friction terms (expansion of the universe and interaction with matter). These type of models can be motivated from supergravity. Inflation occurs close to the maxima of the potential. As a consequence both inflations are necessarily finite. This opens the interesting possibility that the second inflation has already or is about to end. A first inflation is produced when fluctuations displace the inflaton field from its higher maximum rolling down the potential as in new inflation. Instead of rolling towards a global minimum the inflaton approaches a lower maximum where a second inflation takes place.

We study a supergravity model of inflation essentially depending on one parameter which can be identified with the slope of the potential at the origin. In this type of models the inflaton rolls at high energy from negative values and a positive curvature potential. At some point defined by the equation $\eta=1$ inflation starts. The potential curvature eventually changes to negative values and inflation ends when $\eta=-1$. No spontaneous symmetry breaking type mechanism for inflation of the new type to occur is here required. The model naturally gives a bounded $\it{total}$ number of e-folds which is typically close to the required number for observable inflation and it is independent of the initial conditions for the inflaton. The energy scale introduced is fixed by the amplitude of the anisotropies and is of the order of the supersymmetry breaking scale. The model can also accommodate an spectral index bigger or smaller than one without extreme fine tuning. We show that it is possible to obtain reasonable numbers for cosmological parameters and, as an example, we reproduce values obtained recently by Tegmark $\it{et. al.}$, from WMAP and SDSS data alone.

We present a realistic dark energy model derived from particle physics. Our model has essentially no free parameters and has an equivalent fit to the observational data (CMB, SN1a and LSS) as LCDM and a better fit than the best effective $w(z)$ model. With the lack of a clear determination of the cosmological parameters theoretical considerations should be taken seriously to distinguish between dark energy models.

The result from the SN1a projects suggest that the dark energy can be represented by a fluid with w<-1. However, it is commonly argued that a fluid with |w|>1 contradicts causality. Here, we will show that a fluid with |w|>1 does not contradict causality if $w$ is not constant. Scalar field are the most promising candidates for describing the dark energy and they do not have a constant equation of state parameter w. Scalar potentials may lead to regions where w is larger or smaller than one and even regions where the group velocity $d p/d\rho$ diverges. We study the evolution of scalar field perturbations and we show that the "sound speed" is always smaller than the speed of light independently of the value of w=p/\rho or d p/d\rho. In general, it is neither the phase velocity nor the group velocity that gives the "sound speed". In the analysis we include the special cases of $|w|>1$, tachyonic, phantom, k-essence scalar fields and curvature corrections. Our results show that scalar fields do not contradict causality as long as there is dispersion.

We study the cosmology of canonically normalized scalar fields that lead to an equation of state parameter of w_\phi=p_\phi/\rho_\phi<-1 without violating the weak energy condition: rho=\Sigma_i\rho_i \geq 0 and \rho_i+p_i\geq 0. This kind of behavior requires a negative scalar potential V, widely predicted in particle physics. We show that the energy density \rho_\phi=E_k+V takes negative values with an equation of state with w_\phi < -1. However, the net effect of the $\phi$ field on the scale factor is to decelerate it giving a total equation of state parameter w=p/\rho> w_b=p_b/\rho_b, where \rho_b stands for any kind of energy density with -1\leq w_b \leq 1. The fact that \rho_\phi<0 allows, at least in principle, to have a small cosmological constant or quintessence today as the cancellation of high energy scales such as the electroweak or susy breaking scale. While V is negative |\rho_\phi| is smaller than the sum of all other energy densities regardless of the functional form of the potential V. We show that the existence of a negative potential leads, inevitable, to a collapsing universe, i.e. to a would be "big crunch". In this picture we would still be living in the expanding universe.

We study the possibility that a dark group, a gauge group with particles interacting with the standard model particles only via gravity, is responsible for containing the dark energy and dark matter required by present day observations. We show that it is indeed possible and we determine the constrains for the dark group. The non-perturbative effects generated by a strong gauge coupling constant can de determined and a inverse power law scalar potential IPL for the dark meson fields is generated parameterizing the dark energy. On the other hand it is the massive particles, e.g. dark baryons, of the dark gauge group that give the corresponding dark matter. The mass of the dark particles is of the order of the condensation scale $\Lambda_c$ and the temperature is smaller then the photon's temperature. The dark matter is of the warm matter type. The only parameters of the model are the number of particles of the dark group. The allowed values of the different parameters are severely restricted. The dark group energy density at $\Lambda_c$ must be $ \Omdgc \leq 0.17$ and the evolution and acceptable values of dark matter and dark energy leads to a constrain of $\Lmc$ and the IPL parameter $n$ giving $\Lambda_c=O(1-10^3) eV$ and $0.28 \leq n \leq 1.04$.

We present model independent determination of the CMB from any kind of fluid that has an equation of state taking four different values. The first region has $w=1/3$, the second $w=1$, the third $w=-1$ while the last one has $-1<w=cte<-2/3$. This kind of dynamical $w$ contains as a limit the cosmological constant and tracker models. We derive the model independent evolution of $\wp$, for scalar fields, and we see that it remains most of the time in either of its three extremal values given by $\wp=1,-1,w_{tr}$. This "varying" $w$ is the generic behavior of scalar fields, quintessence, and we determine the size of the different regions by solving the dynamical equations in a model independent way. The dynamical $w$ models have a better fit to CMB data then the cosmological constant and the tracker models. We determine the effect of having the first two regions $w=1/3,1$ and depending on the size of these periods they can be observed in the CMB. In general, the CMB spectrum sets a lower limit to $\D N_T$ and to the phase transition scale $\Lm_c$. For smaller $\D N_T$ the CMB peaks are moved to the right of the spectrum and the hight increases considerably. For $\Ompi=0.1$ the CMB sets a lower limit to the phase transition scale $\Lm_c\geq 0.2 eV$.

We present a particle physics quintessence model which agrees well with existing cosmological data, including the position of the acoustic peaks. This model has an inverse power law potential (IPL) with $n=18/7\sim 2.57$ and it gives $\weff=-0.75$, an acoustic scale $l_A=307$ and a density contrast $\s8=0.95$. Models with $n>1$ have been said to be disfavored by the analysis of the acoustic peaks. However, the results are not correct. The main reason is that the tracker approximation has been used in deriving the IPL constrains and for $n<5$ the scalar field has not reached its tracker value by present day. The model can be derived from particle physics, using Affleck-Dine-Seiberg "ADS" superpotential, for a non-abelian gauge group with $N_c=8, N_f=1$. The advantage of having $N_f=1$ is that there is only one degree of freedom below the condensation scale given by the condensate (quintessence) $\phi^2=<Q\tilde Q>$ field. The condensation scale is at $1 GeV$ a very interesting scale since it connects the quintessence "Q" with the standard model "SM" scale. The similarity in energy scales between $Q$ and $SM$ scale gives an "explanation" to the coincidence problem. The fact that only recently the universe is accelerating is a natural consequence of the Q scale and the evolution of $\phi$.

We show that the condensates of a non-abelian gauge group, unified with the standard model gauge groups, can parameterize the present day cosmological constant and play the role of quintessence. The models agree with SN1a and recent CMB analysis. These models have {\it no free} parameters. Even the initial energy density at the unification scale and at the condensation scale are fixed by the number of degrees of freedom of the gauge group (i.e. by $N_c, N_f$). The values of $N_c, N_f$ are determined by imposing gauge coupling unification and the number of models is quite limited. Using Affleck-Dine-Seiberg superpotential one obtains a scalar potential $V=\Lm_c^{4+n}\phi^{-n}$. Models with $2<n<4.27$ or equivalently $2\times 10^{-2} GeV < \Lm_c < 6 \times 10^3 GeV$ do not satisfy the unification constrain. In fact, there are only three models and they have an inverse power potential with $6/11 \leq n \leq 2/3$. Imposing primordial nucleosynthesis bounds the preferred model has $N_c=3, N_f=6$, with $n=2/3$, a condensation scale $\Lm_c=4.2\times 10^{-8} GeV$ and $\wpo=-0.90$ with an average value $\weff=-0.93$. Notice that the tracker solution is not a good approximation since it has $w_{tr}=-\fr{2}{n+2}=-0.75$ for $n=2/3$. We study the evolution of all fields from the unification scale and we calculate the relevant cosmological quantities. We also discuss the supersymmetry breaking mechanism which is relevant for these models.

We study the cosmological evolution of scalar fields that arise from a phase transition at some energy scale $\Lm_c$. We focus on negative power potentials given by $V=c\Lm_c^{4+n}\phi^{-n}$ and restrict the cosmological viable values of $\Lm_c$ and $n$. We make a complete analysis of $V$ and impose $SN1a$ conditions on the different cosmological parameters. The cosmological observations ruled out models where the scalar field has reached its attractor solution. For models where this is not the case, the analytic approximated solutions are not good enough to determine whether a specific model is phenomenologically viable or not and the full differential equations must be numerically solved. The results are not fine tuned since a change of 45% on the initial conditions does not spoil the final results. We also determine the values of $N_c, N_f$ that give a condensation scale $\Lm_c$ consistent with gauge coupling unification, leaving only four models that satisfy unification and SN1a constraints.

We show that a positive accelerating universe can be obtained simply by the dynamics of a non-abelian gauge group. It is the condensates of the chiral fields that obtain a negative power potential, below the condensation scale, and allow for a quintessence interpretation of these fields. The only free parameters in this model are $N_c$ and $N_f$ and the number of dynamically gauge singlet bilinear fields $\phi$ generated below the condensation scale. We show that it is possible to have unification of all coupling constants, including the standard and non standard model couplings, while having an acceptable phenomenology of $\phi$ as the cosmological constant. This is done without any fine tuning of the initial conditions. The problem of coincidence (why the universe has only recently started an accelerating period) is not solved but it is put at the same level as what the particle content of the standard model is.

We show that the evolution of the quintessence energy density $\Omega_Q$ is model independent in an accelerating universe. The accelerating behaviour has lasted at most 0.5 e-folds of expansion assuming a present day value of $\Omega_Q=2/3$. The generic evolution differs from the exact solution by a percentage given by $r=2\gm_Q$, with $\gm_Q=(\rho_Q+p_Q)/\rho_Q$. For a small $\gm_Q$ the evolution remains a good approximation for a long period. Nucleosynthesis bounds on $\Omega_Q$ suggest that the model independent solution is valid for at least 12 e-folds of expansion which includes then the scale of radiation and matter equality. We can, therefore, establish model independent conditions on the cosmological parameters at some factor scale $a_i \ll a_o$. Finally, we discuss the relevance of this result to the cosmological coincidence problem.

We study the cosmological evolution of string/M moduli fields T. We use T-duality to fix the potential and show that the superpotential W is a function of the duality invariant function j(T) only. If W is given as a finite polynomial of j then the moduli fields {\it do not} give an accelerating universe, i.e. they {\it cannot} be used as quintessence. Furthermore, at T >>1 the potential is given by a double exponential potential V \simeq e^{-a e^{\sqrt{2} T}} leading to a fast decaying behaviour at large times. For moduli potentials with a finite v.e.v. of T the energy density redshift is model dependent but if T has a finite mass, m < \infty, then the moduli energy density redshifts faster or equal to matter. Only if the moduli mass is infinite can the moduli energy density dominate the universe independently of the initial conditions.

We study the cosmological evolution of scalar fields with arbitrary potentials in the presence of a barotropic fluid (matter or radiation) without making any assumption on which term dominates. We determine what kind of potentials V(phi) permits a quintessence interpretation of the scalar field phi and to obtain interesting cosmological results. We show that all model dependence is given in terms of lambda= - V'/V only and we study all possible asymptotic limits: lambda approaching zero, a finite constant or infinity. We determine the equation of state dynamically for each case. For the first class of potentials, the scalar field quickly dominates the universe behaviour, with an inflationary equation of state allowing for a quintessence interpretation. The second case gives the extensively studied exponential potential. While in the last case, when lambda approaches infinity, if it does not oscillate then the energy density redshifts faster than the barotropic fluid but if lambda oscillates then the energy density redshift depends on the specific potential.

We investigate the possibility of obtaining inflationary solutions of the slow roll type from a low energy Lagrangian coming from superstrings. The advantage of such an approach is that in these theories the scalar potential has only one free parameter (the Planck scale) and therefore no unnatural fine tuning may be accommodated. We find that in any viable scheme the dilaton and the moduli fields have to be stabilized and that before this happens, no other field may be used as the inflaton. Then inflation may occur due to chiral matter fields. Demanding that the potential terms associated with the chiral fields do not spoil the dilaton and moduli minimization leads to severe constraints on the magnitude of the density fluctuations.

We study the possibility of canceling the cosmological constant in supergravity string models. We show that with a suitable choice of superpotential the vacuum energy may vanish with the dilaton field at its minimum and supersymmetry broken with a large hierarchy. We derive the condition for which the introduction of a chiral potential, e.g. the inflaton potential, does not destabilize the dilaton field even in the region where the scalar potential takes positive values. This allows for an inflationary potential with the dilaton frozen at its minimum.