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 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 a new parameterization for the equation of state (EoS) $\omega_X=P_X/\rho_X$, which can reproduce a $f(R)$-like evolution with a precision between $[0.5\%-0.8\%]$ over the numerical solutions. Also, our proposal can render a variety of popular $f(R)$ models that are considered as viable candidates for the cosmic late time acceleration. By using observational data from baryonic acoustic oscillations, supernovae and cosmic chronometers we investigate the constraints on the new EoS parameters. This proposal set a EoS formulation which can be used in an efficient way and makes a good candidate to be implemented in a variety of surveys in order to test the $f(R)$ generic behaviour.

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 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.