The Wide Field InfraRed Survey Telescope (WFIRST) was the highest ranked large space-based mission of the 2010 New Worlds, New Horizons decadal survey. It is now a NASA mission in formulation with a planned launch in the mid-2020s. A primary mission objective is to precisely constrain the nature of dark energy through multiple probes, including Type Ia supernovae (SNe Ia). Here, we present the first realistic simulations of the WFIRST SN survey based on current hardware specifications and using open-source tools. We simulate SN light curves and spectra as viewed by the WFIRST wide-field channel (WFC) imager and integral-field channel (IFC) spectrometer, respectively. We examine 11 survey strategies with different time allocations between the WFC and IFC, two of which are based upon the strategy described by the WFIRST Science Definition Team, which measures SN distances exclusively from IFC data. We propagate statistical and, crucially, systematic uncertainties to predict the Dark Energy Task Force figure of merit (FoM) for each strategy. Of the strategies investigated, we find the most successful to be WFC-focused. However, further work in constraining systematics is required to fully optimize the use of the IFC. Even without improvements to other cosmological probes, the WFIRST SN survey has the potential to increase the FoM by more than an order of magnitude from the current values. Although the survey strategies presented here have not been fully optimized, these initial investigations are an important step in the development of the final hardware design and implementation of the WFIRST mission.

The light curves of Type Ia supernovae are routinely used to constrain cosmology models. Driven by radioactive decay of 56Ni, the light curves steadily decline over time, but >150 days past explosion, the near-infrared portion is poorly characterized. We report a year-long plateau in the near-infrared light curve at 150-500 days, followed by a second decline phase accompanied by a possible appearance of [Fe I] emission lines. This near-infrared plateau contrasts sharply with Type IIP plateaus and requires a new physical mechanism. We suggest a such as masking of the "near-infrared catastrophe," a predicted yet unobserved sharp light-curve decline, by scattering of ultraviolet photons to longer wavelengths. The transition off the plateau could be due to a change in the dominant ionization state of the supernova ejecta. Our results shed new light on the complex radiative transfer processes that take place in Type Ia supernovae and enhance their use as "standard candles."

We describe the simulated data sample for the "Photometric LSST Astronomical Time Series Classification Challenge" (PLAsTiCC), a publicly available challenge to classify transient and variable events that will be observed by the Large Synoptic Survey Telescope (LSST), a new facility expected to start in the early 2020s. The challenge was hosted by Kaggle, ran from 2018 September 28 to 2018 December 17, and included 1,094 teams competing for prizes. Here we provide details of the 18 transient and variable source models, which were not revealed until after the challenge, and release the model libraries at https://doi.org/10.5281/zenodo.2612896. We describe the LSST Operations Simulator used to predict realistic observing conditions, and we describe the publicly available SNANA simulation code used to transform the models into observed fluxes and uncertainties in the LSST passbands (ugrizy). Although PLAsTiCC has finished, the publicly available models and simulation tools are being used within the astronomy community to further improve classification, and to study contamination in photometrically identified samples of type Ia supernova used to measure properties of dark energy. Our simulation framework will continue serving as a platform to improve the PLAsTiCC models, and to develop new models.

We present the analysis underpinning the measurement of cosmological parameters from 207 spectroscopically classified type Ia supernovae (SNe Ia) from the first three years of the Dark Energy Survey Supernova Program (DES-SN), spanning a redshift range of 0.017<$z$<0.849. We combine the DES-SN sample with an external sample of 122 low-redshift ($z$<0.1) SNe Ia, resulting in a "DES-SN3YR" sample of 329 SNe Ia. Our cosmological analyses are blinded: after combining our DES-SN3YR distances with constraints from the Cosmic Microwave Background (CMB; Planck Collaboration 2016), our uncertainties in the measurement of the dark energy equation-of-state parameter, $w$, are .042 (stat) and .059 (stat+syst) at 68% confidence. We provide a detailed systematic uncertainty budget, which has nearly equal contributions from photometric calibration, astrophysical bias corrections, and instrumental bias corrections. We also include several new sources of systematic uncertainty. While our sample is <1/3 the size of the Pantheon sample, our constraints on $w$ are only larger by 1.4$\times$, showing the impact of the DES SN Ia light curve quality. We find that the traditional stretch and color standardization parameters of the DES SNe Ia are in agreement with earlier SN Ia samples such as Pan-STARRS1 and the Supernova Legacy Survey. However, we find smaller intrinsic scatter about the Hubble diagram (0.077 mag). Interestingly, we find no evidence for a Hubble residual step ( 0.007 $\pm$ 0.018 mag) as a function of host galaxy mass for the DES subset, in 2.4$\sigma$ tension with previous measurements. We also present novel validation methods of our sample using simulated SNe Ia inserted in DECam images and using large catalog-level simulations to test for biases in our analysis pipelines.

We present the first cosmological parameter constraints using measurements of type Ia supernovae (SNe Ia) from the Dark Energy Survey Supernova Program (DES-SN). The analysis uses a subsample of 207 spectroscopically confirmed SNe Ia from the first three years of DES-SN, combined with a low-redshift sample of 122 SNe from the literature. Our "DES-SN3YR" result from these 329 SNe Ia is based on a series of companion analyses and improvements covering SN Ia discovery, spectroscopic selection, photometry, calibration, distance bias corrections, and evaluation of systematic uncertainties. For a flat LCDM model we find a matter density Omega_m = 0.331 +_ 0.038. For a flat wCDM model, and combining our SN Ia constraints with those from the cosmic microwave background (CMB), we find a dark energy equation of state w = -0.978 +_ 0.059, and Omega_m = 0.321 +_ 0.018. For a flat w0waCDM model, and combining probes from SN Ia, CMB and baryon acoustic oscillations, we find w0 = -0.885 +_ 0.114 and wa = -0.387 +_ 0.430. These results are in agreement with a cosmological constant and with previous constraints using SNe Ia (Pantheon, JLA).

We describe catalog-level simulations of Type Ia supernova (SN~Ia) light curves in the Dark Energy Survey Supernova Program (DES-SN), and in low-redshift samples from the Center for Astrophysics (CfA) and the Carnegie Supernova Project (CSP). These simulations are used to model biases from selection effects and light curve analysis, and to determine bias corrections for SN~Ia distance moduli that are used to measure cosmological parameters. To generate realistic light curves, the simulation uses a detailed SN~Ia model, incorporates information from observations (PSF, sky noise, zero point), and uses summary information (e.g., detection efficiency vs. signal to noise ratio) based on 10,000 fake SN light curves whose fluxes were overlaid on images and processed with our analysis pipelines. The quality of the simulation is illustrated by predicting distributions observed in the data. Averaging within redshift bins, we find distance modulus biases up to 0.05~mag over the redshift ranges of the low-z and DES-SN samples. For individual events, particularly those with extreme red or blue color, distance biases can reach 0.4~mag. Therefore, accurately determining bias corrections is critical for precision measurements of cosmological parameters. Files used to make these corrections are available at https://des.ncsa.illinois.edu/releases/sn.

The combination of multiple observational probes has long been advocated as a powerful technique to constrain cosmological parameters, in particular dark energy. The Dark Energy Survey has measured 207 spectroscopically--confirmed Type Ia supernova lightcurves; the baryon acoustic oscillation feature; weak gravitational lensing; and galaxy clustering. Here we present combined results from these probes, deriving constraints on the equation of state, $w$, of dark energy and its energy density in the Universe. Independently of other experiments, such as those that measure the cosmic microwave background, the probes from this single photometric survey rule out a Universe with no dark energy, finding $w=-0.80^{+0.09}_{-0.11}$. The geometry is shown to be consistent with a spatially flat Universe, and we obtain a constraint on the baryon density of $\Omega_b=0.069^{+0.009}_{-0.012}$ that is independent of early Universe measurements. These results demonstrate the potential power of large multi-probe photometric surveys and pave the way for order of magnitude advances in our constraints on properties of dark energy and cosmology over the next decade.

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.

Over the past century, major advances in astronomy and astrophysics have been largely driven by improvements in instrumentation and data collection. With the amassing of high quality data from new telescopes, and especially with the advent of deep and large astronomical surveys, it is becoming clear that future advances will also rely heavily on how those data are analyzed and interpreted. New methodologies derived from advances in statistics, computer science, and machine learning are beginning to be employed in sophisticated investigations that are not only bringing forth new discoveries, but are placing them on a solid footing. Progress in wide-field sky surveys, interferometric imaging, precision cosmology, exoplanet detection and characterization, and many subfields of stellar, Galactic and extragalactic astronomy, has resulted in complex data analysis challenges that must be solved to perform scientific inference. Research in astrostatistics and astroinformatics will be necessary to develop the state-of-the-art methodology needed in astronomy. Overcoming these challenges requires dedicated, interdisciplinary research. We recommend: (1) increasing funding for interdisciplinary projects in astrostatistics and astroinformatics; (2) dedicating space and time at conferences for interdisciplinary research and promotion; (3) developing sustainable funding for long-term astrostatisics appointments; and (4) funding infrastructure development for data archives and archive support, state-of-the-art algorithms, and efficient computing.

While Type Ia Supernovae (SNe Ia) are one of the most mature cosmological probes, the next era promises to be extremely exciting in the number of different ways SNe Ia are used to measure various cosmological parameters. Here we review the experiments in the 2020s that will yield orders of magnitudes more SNe Ia, and the new understandings and capabilities to constrain systematic uncertainties at a level to match these statistics. We then discuss five different cosmological probes with SNe Ia: the conventional Hubble diagram for measuring dark energy properties, the distance ladder for measuring the Hubble constant, peculiar velocities and weak lensing for measuring sigma8 and strong-lens measurements of H0 and other cosmological parameters. For each of these probes, we discuss the experiments that will provide the best measurements and also the SN Ia-related systematics that affect each one.

We analyze a set of 89 Type Ia supernovae (SN Ia) that have both optical and near-infrared (NIR) photometry to derive distances and construct low redshift ($z < 0.04$) Hubble diagrams. We construct mean light curve (LC) templates using a hierarchical Bayesian model. We explore both Gaussian process (GP) and template methods for fitting the LCs and estimating distances, while including peculiar velocity and photometric uncertainties. For the 56 SN Ia with both optical and NIR observations near maximum light, the GP method yields a NIR-only Hubble-diagram with a RMS of $0.117 \pm 0.014$ mag when referenced to the NIR maxima. For each NIR band, a comparable GP method RMS is obtained when referencing to NIR-max or B-max. Using NIR LC templates referenced to B-max yields a larger RMS value of $0.138 \pm 0.014$ mag. Fitting the corresponding optical data using standard LC fitters that use LC shape and color corrections yields larger RMS values of $0.179 \pm 0.018$ mag with SALT2 and $0.174 \pm 0.021$ mag with SNooPy. Applying our GP method to subsets of SN Ia NIR LCs at NIR maximum light, even without corrections for LC shape, color, or host-galaxy dust reddening, provides smaller RMS in the inferred distances, at the $\sim 2.3 - 4.1\sigma$ level, than standard optical methods that do correct for those effects. Our ongoing RAISIN program on the Hubble Space Telescope will exploit this promising infrared approach to limit systematic errors when measuring the expansion history of the universe to constrain dark energy.

We present spectroscopy from the first three seasons of the Dark Energy Survey Supernova Program (DES-SN). We describe the supernova spectroscopic program in full: strategy, observations, data reduction, and classification. We have spectroscopically confirmed 307 supernovae, including 251 type Ia supernovae (SNe Ia) over a redshift range of $0.017 < z < 0.85$. We determine the effective spectroscopic selection function for our sample, and use it to investigate the redshift-dependent bias on the distance moduli of SNe Ia we have classified. We also provide a full overview of the strategy, observations, and data products of DES-SN, which has discovered 12,015 likely supernovae during these first three seasons. The data presented here are used for the first cosmology analysis by DES-SN ('DES-SN3YR'), the results of which are given in DES Collaboration (2018a).

Empirical models of supernova (SN) spectral energy distributions (SEDs) are widely used for SN survey simulations and photometric classifications. The existing library of SED models has excellent optical templates but limited, poorly constrained coverage of ultraviolet (UV) and infrared (IR) wavelengths. However, both regimes are critical for the design and operation of future SN surveys, particularly at IR wavelengths that will be accessible with the James Webb Space Telescope (JWST) and the Wide-Field Infrared Survey Telescope (WFIRST). We create a public repository of improved empirical SED templates using a sampling of Type Ia and core-collapse (CC) photometric light curves to extend the Type Ia parameterized SALT2 model and a set of SN Ib, SN Ic, and SN II SED templates into the UV and near-IR. We apply this new repository of extrapolated SN SED models to examine how future surveys can discriminate between CC and Type Ia SNe at UV and IR wavelengths, and present an open-source software package written in Python, SNSEDextend, that enables a user to generate their own extrapolated SEDs.

We present the interesting coincidence of cosmology and astrophysics that points toward a dimensionless age of the universe H_0*t_0 that is close to one. Despite cosmic deceleration for 9 Gyr and acceleration since then, we find H_0t*_0 = 0.96 +/- 0.01 for the LCDM model that fits SN Ia data from Pan-STARRS, CMB power spectra, and baryon acoustic oscillations. Similarly, astrophysical measures of stellar ages and the Hubble constant derived from redshifts and distances point to H_0*t ~ 1.0 +/- 0.1$. The wide range of possible values for H_0*t_0 realized during comic evolution means that we live at what appears to be a special time. This "synchronicity problem" is not precisely the same as the usual Coincidence problem because there are combinations of Omega_Matter and Omega_Lambda for which the usual coincidence problem holds but for which H_0*t_0 is not close to 1.

The braneworld model proposed by Dvali, Gabadadze and Porrati (DGP) leads to an accelerated universe without cosmological constant or other form of dark energy for the positive branch $(\epsilon =+1)$. For the negative branch $(\epsilon =-1)$ we have investigated the behavior of a model with an holographic Ricci-like dark energy and dark matter, where the IR cutoff takes the form $\alpha H^2 + \beta \dot{H}$, being $H$ the Hubble parameter and $\alpha$, $\beta$ positive constants of the model. We perform an analytical study of the model in the late-time dark energy dominated epoch, where we obtain a solution for $r_cH(z)$, where $r_c$ is the leakage scale of gravity into the bulk, and conditions for the negative branch on the holographic parameters $\alpha$ and $\beta$, in order to hold the conditions of weak energy and accelerated universe. On the other hand, we compare the model versus the late-time cosmological data using the latest type Ia supernova sample of the Joint Light-curve Analysis (JLA), in order to constraint the holographic parameters in the negative branch, as well as $r_cH_0$ in the positive branch, where $H_0$ is the Hubble constant. We find that the model has a good fit to the data and that the most likely values for $(r_cH_0, \alpha, \beta)$ lie in the permitted region found from an analytical solution in a dark energy dominated universe. We give a justification to use holographic cut-off in 4D for the dark energy in the 5 dimensional DGP model. Finally, using the Bayesian Information Criterion we find that this model it is disfavored compared with the flat $\Lambda$CDM model.

By means of a combined use of the type Ia supernovae and H(z) data tests, together with the study of the asymptotic properties in the equivalent phase space - through the use of the dynamical systems tools - we demonstrate that the bulk viscous matter-dominated scenario is not a good model to explain the accepted cosmological paradigm, at least, under the parametrization of bulk viscosity considered in this paper. The main objection against such scenarios is the absence of conventional radiation and matter-dominated critical points in the phase space of the model. This entails that radiation and matter dominance are not generic solutions of the cosmological equations, so that these stages can be implemented only by means of unique and very specific initial conditions, i. e., of very unstable particular solutions. Such a behavior is in marked contradiction with the accepted cosmological paradigm which requires of an earlier stage dominated by relativistic species, followed by a period of conventional non-relativistic matter domination, during which the cosmic structure we see was formed. Also, we found that the bulk viscosity is positive just until very late times in the cosmic evolution, around z < 1. For earlier epochs it is negative, been in tension with the local second law of thermodynamics.

We explore a cosmological model composed by a dark matter fluid interacting with a dark energy fluid. The interaction term has the non-linear lambda*(rho_m)^alpha * (rho_e)^beta form, where rho_m and rho_e are the energy densities of the dark matter and dark energy, respectively. The parameters alpha and beta are in principle not constrained to take any particular values, and were estimated from observations. We perform an analytical study of the evolution equations, finding the fixed points and their stability properties in order to characterize suitable physical regions in the phase space of the dark matter and dark energy densities. The constants (lambda, alpha, beta) as well as w_m and w_e of the EoS of dark matter and dark energy respectively, were estimated using the cosmological observations of the type Ia supernovae and the Hubble expansion rate H(z) data sets. We find that the best estimated values for the free parameters of the model correspond to a warm dark matter interacting with a phantom dark energy component, with a well goodness-of-fit to data. However, using the Bayesian Information Criterion (BIC) we find that this model is overcame by a warm dark matter -- phantom dark energy model without interaction, as well as by the LCDM model. We find also a large dispersion on the best estimated values of the (lambda, alpha, beta) parameters, so even if we are not able to set strong constraints on their values, given the goodness-of-fit to data of the model, we find that a large variety of theirs values are well compatible with the observational data used.

We revise the conditions for the physical viability of a cosmological model in which dark matter has bulk viscosity and also interacts with dark energy. We have also included radiation and baryonic matter components; all matter components are represented by perfect fluids, except the dark matter, that is treated as an imperfect fluid. We impose upon the model the condition of a complete cosmological dynamics that results in an either null or negative bulk viscosity, but the latter also disagrees with the Local Second Law of Thermodynamics. The model is also compared with cosmological observations at different redshifts: type Ia supernova, the shift parameter of CMB, the acoustic peak of BAO, and the Hubble parameter H(z). In general, observations consistently point out to a negative value of the bulk viscous coefficient, and in overall the fitting procedure shows no preference for the model over the standard LCDM model.

We explore the cosmological constraints on the parameter w_dm of the dark matter barotropic equation of state (EoS) to investigate the "warmness" of the dark matter fluid. The model is composed by the dark matter and dark energy fluids in addition to the radiation and baryon components. We constrain the values of w_dm using the latest cosmological observations that measure the expansion history of the Universe. When w_dm is estimated together with the parameter w_de of the barotropic EoS of dark energy we found that the cosmological data favor a value of w_dm = 0.006 +- 0.001, suggesting a -warm- dark matter, and w_de= -1.11 +- 0.03$ that corresponds to a phantom dark energy, instead of favoring a cold dark matter and a cosmological constant (w_dm = 0, w_de = -1). When w_dm is estimated alone but assuming w_de = -1, -1.1, -0.9, we found w_dm = 0.009 +- 0.002, 0.006 +- 0.002, 0.012 +- 0.002 respectively, where the errors are at 3 sigma (99.73%), i.e., w_dm > 0 with at least 99.73% of confidence level. When (w_dm, \Omega_dm0) are constrained together, the best fit to data corresponds to (w_dm=0.005 +- 0.001, \Omega_dm0 = 0.223 +- 0.008) and with the assumption of w_de = -1.1 instead of a cosmological constant (i.e., w_de = -1). With these results we found evidence of w_dm > 0 suggesting a -warm- dark matter, independent of the assumed value for w_{\rm de}, but where values w_de < -1 are preferred by the observations instead of the cosmological constant. These constraints on w_dm are consistent with perturbative analyses done in previous works.

We study a cosmological model composed of a dark energy fluid interacting with a viscous matter fluid in a spatially flat Universe. The matter component represents the baryon and dark matter and it is taken into account, through a bulk viscosity, the irreversible process that the matter fluid undergoes because of the accelerated expansion of the universe. The bulk viscous coefficient is assumed to be proportional to the Hubble parameter. The radiation component is also taken into account in the model. The model is constrained using the type Ia supernova observations, the shift parameter of the CMB, the acoustic peak of the BAO and the Hubble expansion rate, to constrain the values of the barotropic index of dark energy and the bulk viscous coefficient. It is found that the bulk viscosity is constrained to be negligible (around zero) from the observations and that the barotropic index for the dark energy to be negative and close to zero too, indicating a phantom energy.

We explore the viability of a bulk viscous matter-dominated Universe to explain the present accelerated expansion of the Universe. The model is composed by a pressureless fluid with bulk viscosity of the form \zeta = \zeta_0 + \zeta_1 * H where \zeta_0 and \zeta_1 are constants and H is the Hubble parameter. The pressureless fluid characterizes both the baryon and dark matter components. We study the behavior of the Universe according to this model analyzing the scale factor as well as some curvature scalars and the matter density. On the other hand, we compute the best estimated values of \zeta_0 and \zeta_1 using the type Ia Supernovae (SNe Ia) probe. We find that from all the possible scenarios for the Universe, the preferred one by the best estimated values of (\zeta_0, \zeta_1) is that of an expanding Universe beginning with a Big- Bang, followed by a decelerated expansion at early times, and with a smooth transition in recent times to an accelerated expansion epoch that is going to continue forever. The predicted age of the Universe is a little smaller than the mean value of the observational constraint coming from the oldest globular clusters but it is still inside of the confidence interval of this constraint. A drawback of the model is the violation of the local second law of thermodynamics in redshifts z >= 1. However, when we assume \zeta_1 = 0, the simple model \zeta = \zeta_0 evaluated at the best estimated value for \zeta_0 satisfies the local second law of thermodynamics, the age of the Universe is in perfect agreement with the constraint of globular clusters, and it also has a Big-Bang, followed by a decelerated expansion with the smooth transition to an accelerated expansion epoch in late times, that is going to continue forever.

We test a cosmological model which the only component is a pressureless fluid with a constant bulk viscosity as an explanation for the present accelerated expansion of the universe. We classify all the possible scenarios for the universe predicted by the model according to their past, present and future evolution and we test its viability performing a Bayesian statistical analysis using the SCP ``Union'' data set (307 SNe Ia), imposing the second law of thermodynamics on the dimensionless constant bulk viscous coefficient \zeta and comparing the predicted age of the universe by the model with the constraints coming from the oldest globular clusters. The best estimated values found for \zeta and the Hubble constant Ho are: \zeta=1.922 \pm 0.089 and Ho=69.62 \pm 0.59 km/s/Mpc with a \chi^2=314. The age of the universe is found to be 14.95 \pm 0.42 Gyr. We see that the estimated value of Ho as well as of \chi^2 are very similar to those obtained from LCDM model using the same SNe Ia data set. The estimated age of the universe is in agreement with the constraints coming from the oldest globular clusters. Moreover, the estimated value of \zeta is positive in agreement with the second law of thermodynamics (SLT). On the other hand, we perform different forms of marginalization over the parameter Ho in order to study the sensibility of the results to the way how Ho is marginalized. We found that it is almost negligible the dependence between the best estimated values of the free parameters of this model and the way how Ho is marginalized in the present work. Therefore, this simple model might be a viable candidate to explain the present acceleration in the expansion of the universe.

We present and constrain a cosmological model where the only component is a pressureless fluid with bulk viscosity as an explanation for the present accelerated expansion of the universe. We study the particular model of a bulk viscosity coefficient proportional to the Hubble parameter. The model is constrained using the SNe Ia Gold 2006 sample, the Cosmic Microwave Background (CMB) shift parameter R, the Baryon Acoustic Oscillation (BAO) peak A and the Second Law of Thermodynamics (SLT). It was found that this model is in agreement with the SLT using only the SNe Ia test. However when the model is constrained using the three cosmological tests together (SNe+CMB+BAO) we found: 1.- The model violates the SLT, 2.- It predicts a value of H_0 \approx 53 km sec^{-1} Mpc^{-1} for the Hubble constant, and 3.- We obtain a bad fit to data with a \chi^2_{min} \approx 532. These results indicate that this model is viable just if the bulk viscosity is triggered in recent times.

We present and constrain a cosmological model which component is a pressureless fluid with bulk viscosity as an explanation for the present accelerated expansion of the universe. We study the particular model of a constant bulk viscosity coefficient \zeta_m. The possible values of \zeta_m are constrained using the cosmological tests of SNe Ia Gold 2006 sample, the CMB shift parameter R from the three-year WMAP observations, the Baryon Acoustic Oscillation (BAO) peak A from the Sloan Digital Sky Survey (SDSS) and the Second Law of Thermodynamics (SLT). It was found that this model is in agreement with the SLT using only the SNe Ia test. However when the model is submitted to the three cosmological tests together (SNe+CMB+BAO) the results are: 1.- the model violates the SLT, 2.- predicts a value of H_0 \approx 53 km sec^{-1} Mpc^{-1} for the Hubble constant, and 3.- we obtain a bad fit to data with a \chi^2_{min} \approx 400 (\chi^2_{d.o.f.} \approx 2.2). These results indicate that this model is ruled out by the observations.