As a physical and sufficient compression of the full CMB data, the CMB distance priors, or shift parameters, have been widely used and provide a convenient way to include CMB data when obtaining cosmological constraints. In this paper, we revisit this data vector and examine its stability under different cosmological models. We find that the CMB distance priors are an accurate substitute for the full CMB data when probing dark energy dynamics. This is true when the primordial power spectrum model is directly generalized from the power spectrum of the model used in the derivation of the distance priors from the CMB data. We discover a difference when a non-flat model with the untilted primordial inflation power spectrum is used to measure the distance priors. This power spectrum is a radical change from the more conventional tilted primordial power spectrum and violates fundamental assumptions for the reliability of the CMB shift parameters. We also investigate the performance of CMB distance priors when the sum of neutrino masses $\sum m_{\nu}$ and the effective number of relativistic species $N_{\text{eff}}$ are allowed to vary. Our findings are consistent with earlier results: the neutrino parameters can change the measurement of the sound horizon from CMB data, and thus the CMB distance priors. We find that when the neutrino model is allowed to vary, the cold dark matter density $\omega_{c}$ and $N_{\text{eff}}$ need to be included in the set of parameters that summarize CMB data, in order to reproduce the constraints from the full CMB data. We present an updated and expanded set of CMB distance priors which can reproduce constraints from the full CMB data within $1\sigma$, and are applicable to models with massive neutrinos, as well as non-standard cosmologies.

We perform a simulation with Galacticus, a semi-analytical galaxy formation model, to predict the number counts of H$\alpha$ and [OIII] emitting galaxies. With a state-of-the-art N-body simulation, UNIT, we first calibrate Galacticus with the current observation of H$\alpha$ luminosity function. The resulting model coupled with a dust attenuation model, can reproduce the current observations, including the H$\alpha$ luminosity function from HiZELS and number density from WISP. We extrapolate the model prediction to higher redshift and the result is found to be consistent with previous investigations. We then use the same galaxy formation model to predict the number counts for [OIII] emitting galaxies. The result provides further validation of our galaxy formation model and dust model. We present number counts of H$\alpha$ and [OIII] emission line galaxies for three different line flux limits: $5\times10^{-17}$erg/s/cm$^{2}$, $1\times10^{-16}$erg/s/cm$^{2}$ (6.5$\sigma$ nominal depth for WFIRST GRS), and $2\times10^{-16}$erg/s/cm$^{2}$ (3.5$\sigma$ depth of Euclid GRS). At redshift $2<z<3$, our model predicts that WFIRST can observe hundreds of [OIII] emission line galaxies per square degree with a line flux limit of $1\times10^{-16}$erg/s/cm$^{2}$. This will provide accurate measurement of large scale structure to probe dark energy over a huge cosmic volume to an unprecedented high redshift. Finally, we compare the flux ratio of H$\alpha$/[OIII] within the redshift range of $0<z<3$. Our results show the known trend of increasing H$\alpha$/[O III] flux ratio with H$\alpha$ flux at low redshift, which becomes a weaker trend at higher redshifts.

Models of the spatial distribution of dark matter halos must achieve new levels of precision and accuracy in order to satisfy the requirements of upcoming experiments. In this work, we present a halo bias emulator for modeling the clustering of halos on large scales. It incorporates the cosmological dependence of the bias beyond the mapping of halo mass to peak height. The emulator makes substantial improvements in accuracy compared to the widely used Tinker et al. (2010) model. Halos in this work are defined using an overdensity criteria of 200 relative to the mean background density. Halo catalogs are produced for 40 N-body simulations as part of the Aemulus project at snapshots from z=3 to z=0. The emulator is trained over the mass range $6\times10^{12}-7\times10^{15}\ h^{-1}M_{\odot}$. Using an additional suite of 35 simulations, we determine that the precision of the emulator is redshift dependent, achieving sub-percent levels for a majority of the redshift range. Two additional simulation suites are used to test the ability of the emulator to extrapolate to higher and lower masses. Our high-resolution simulation suite is used to develop an extrapolation scheme in which the emulator asymptotes to the Tinker et al. (2010) model at low mass, achieving ~3% accuracy down to $10^{11}\ h^{-1}M_{\odot}$. Finally, we present a method to propagate emulator modeling uncertainty into an error budget. Our emulator is made publicly available at \url{https://github.com/AemulusProject/bias_emulator}.

Two decades after its discovery, cosmic acceleration remains the most profound mystery in cosmology and arguably in all of physics. Either the Universe is dominated by a form of dark energy with exotic physical properties not predicted by standard model physics, or General Relativity is not an adequate description of gravity over cosmic distances. WFIRST emerged as a top priority of Astro2010 in part because of its ability to address the mystery of cosmic acceleration through both high precision measurements of the cosmic expansion history and the growth of cosmic structures with multiple and redundant probes. We illustrate in this white paper how mission design changes since Astro2010 have made WFIRST an even more powerful dark energy facility and have improved the ability of WFIRST to respond to changes in the experimental landscape. WFIRST is the space-based probe of DE the community needs in the mid-2020s.

Weak lensing of Type Ia supernovae (SNe Ia) is a systematic uncertainty in the use of SNe Ia as standard candles, as well as an independent cosmological probe, if the corresponding magnification distribution can be extracted from data. We study the peak brightness distribution of SNe Ia in the Pantheon sample, and find that the high $z$ sub-sample shows distinct weak lensing signatures compared to the low $z$ subsample: a long tail at the bright end due to high magnifications and a shift of the peak brightness toward the faint end, consistent with findings from earlier work. We have developed a technique to reconstruct the weak lensing magnification distribution of SNe Ia, $p(\mu)$, from the measured SN Ia flux distribution, and applied it to the Pantheon sample. We find that $p(\mu)$ can be reconstructed at a significance better than 2$\sigma$ for the subsample of SNe Ia at $z>0.7$ (124 SNe Ia), and at a lower significance for the SNe Ia at $z>0.9$ (49 SNe Ia), due to the small number of SNe Ia at high redshifts. The large number of $z>1$ SNe Ia from future surveys will enable the use of $p(\mu)$ reconstructed from SNe Ia as an independent cosmological probe.

We present a systematic analysis of the cosmological constraints from the "Pantheon Sample" of 1048 Type Ia Supernovae (SNe Ia) in the redshift range $0.01<z<2.3$ compiled by Scolnic et al. (2018). Applying the flux-averaging method for detecting unknown systematic effects, we find that the "Pantheon" sample has been well calibrated and the bias caused by unknown systematic errors has been minimized. We present the estimate of distances measured from SNe Ia and reconstruct the expansion history of the Universe. The results are in agreement with a simple cosmological constant model and reveals the possible improvements that future SN Ia observations from WFIRST and LSST can target. We have derived distance priors using the CMB data from the Planck 2018 final data release, and combine them with SNe Ia and BAO data, to explore the impact from the systematic errors of SNe Ia on the combined cosmological parameter constraints. Using the combined data set of SNe Ia, BAO, and CMB distance priors, we measure the dark energy density function $X(z)=\rho_X(z)/\rho_X(0)$ as a free function (defined as a cubic spline of its values at $z=0.33, 0.67, 1.0$), along with the cosmological parameters ($\Omega_k$, $\Omega_m$, $\Omega_b$, $H_0$). We find no deviation from a flat Universe dominated by a cosmological constant ($X(z)=1$), and $H_0=68.4\pm 0.9\text{km s}^{-1}\text{Mpc}^{-1}$, straddling the Planck team's measurement of $H_0=67.4\pm 0.5\text{km s}^{-1}\text{Mpc}^{-1}$, and Riess et al. (2018) measurement of $H_0=73.52\pm 1.62\text{km s}^{-1}\text{Mpc}^{-1}$. Adding $H_0=73.52\pm 1.62\text{km s}^{-1}\text{Mpc}^{-1}$ as a prior to the combined data set leads to the time dependence of the dark energy density at $z\sim 0.33$ at 68\% confidence level. Not including the systematic errors on SNe Ia has a similar but larger effect on the dark energy density measurement.

Using the N-body simulations of the AEMULUS Project, we construct an emulator for the non-linear clustering of galaxies in real and redshift space. We construct our model of galaxy bias using the halo occupation framework, accounting for possible velocity bias. The model includes 15 parameters, including both cosmological and galaxy bias parameters. We demonstrate that our emulator achieves ~ 1% precision at the scales of interest, 0.1<r<10 h^{-1} Mpc, and recovers the true cosmology when tested against independent simulations. Our primary parameters of interest are related to the growth rate of structure, f, and its degenerate combination fsigma_8. Using this emulator, we show that the constraining power on these parameters monotonically increases as smaller scales are included in the analysis, all the way down to 0.1 h^{-1} Mpc. For a BOSS-like survey, the constraints on fsigma_8 from r<30 h^{-1} Mpc scales alone are more than a factor of two tighter than those from the fiducial BOSS analysis of redshift-space clustering using perturbation theory at larger scales. The combination of real- and redshift-space clustering allows us to break the degeneracy between f and sigma_8, yielding a 9% constraint on f alone for a BOSS-like analysis. The current AEMULUS simulations limit this model to surveys of massive galaxies. Future simulations will allow this framework to be extended to all galaxy target types, including emission-line galaxies.

The rapidly growing statistical precision of galaxy surveys has lead to a need for ever-more precise predictions of the observables used to constrain cosmological and galaxy formation models. The primary avenue through which such predictions will be obtained is suites of numerical simulations. These simulations must span the relevant model parameter spaces, be large enough to obtain the precision demanded by upcoming data, and be thoroughly validated in order to ensure accuracy. In this paper we present one such suite of simulations, forming the basis for the AEMULUS Project, a collaboration devoted to precision emulation of galaxy survey observables. We have run a set of 75 (1.05 h^-1 Gpc)^3 simulations with mass resolution and force softening of 3.51\times 10^10 (Omega_m / 0.3) ~ h^-1 M_sun and 20 ~ h^-1 kpc respectively in 47 different wCDM cosmologies spanning the range of parameter space allowed by the combination of recent Cosmic Microwave Background, Baryon Acoustic Oscillation and Type Ia Supernovae results. We present convergence tests of several observables including spherical overdensity halo mass functions, galaxy projected correlation functions, galaxy clustering in redshift space, and matter and halo correlation functions and power spectra. We show that these statistics are converged to 1% (2%) for halos with more than 500 (200) particles respectively and scales of r>200 ~ h^-1 kpc in real space or k ~ 3 h Mpc^-1 in harmonic space for z\le 1. We find that the dominant source of uncertainty comes from varying the particle loading of the simulations. This leads to large systematic errors for statistics using halos with fewer than 200 particles and scales smaller than k ~ 4 h^-1 Mpc. We provide the halo catalogs and snapshots detailed in this work to the community at https://AemulusProject.github.io.

Existing models for the dependence of the halo mass function on cosmological parameters will become a limiting source of systematic uncertainty for cluster cosmology in the near future. We present a halo mass function emulator and demonstrate improved accuracy relative to state-of-the-art analytic models. In this work, mass is defined using an overdensity criteria of 200 relative to the mean background density. Our emulator is constructed from the AEMULUS simulations, a suite of 40 N-body simulations with snapshots from z=3 to z=0. These simulations cover the flat wCDM parameter space allowed by recent Cosmic Microwave Background, Baryon Acoustic Oscillation and Type Ia Supernovae results, varying the parameters w, Omega_m, Omega_b, sigma_8, N_{eff}, n_s, and H_0. We validate our emulator using five realizations of seven different cosmologies, for a total of 35 test simulations. These test simulations were not used in constructing the emulator, and were run with fully independent initial conditions. We use our test simulations to characterize the modeling uncertainty of the emulator, and introduce a novel way of marginalizing over the associated systematic uncertainty. We confirm non-universality in our halo mass function emulator as a function of both cosmological parameters and redshift. Our emulator achieves better than 1% precision over much of the relevant parameter space, and we demonstrate that the systematic uncertainty in our emulator will remain a negligible source of error for cluster abundance studies through at least the LSST Year 1 data set.

The extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 14 sample includes 80,118 Luminous Red Galaxies. By combining these galaxies with the high-redshift tail of the BOSS galaxy sample, we form a sample of LRGs at an effective redshift $z=0.72$, covering an effective volume of 0.9~Gpc$^3$. We introduce new techniques to account for spurious fluctuations caused by targeting and by redshift failures which were validated on a set of mock catalogs. This analysis is sufficient to provide a $2.6$\% measurement of spherically averaged BAO, $D_V(z=0.72) = 2353^{+63}_{-61} (r_d/r_{d,\rm{fid}}) h^{-1}$Mpc, at 2.8$\sigma$ of significance. Together with the recent quasar-based BAO measurement at $z=1.5$, and forthcoming Emission Line Galaxy-based measurements, this measurement demonstrates that eBOSS is fulfilling its remit of extending the range of redshifts covered by such measurements, laying the ground work for forthcoming surveys such as the Dark Energy Spectroscopic Survey and Euclid.

We compare a large suite of theoretical cosmological models to observational data from the cosmic microwave background, baryon acoustic oscillation measurements of expansion, Type Ia SNe measurements of expansion, redshift space distortion measurements of the growth of structure, and the local Hubble constant. Our theoretical models include parametrizations of dark energy as well as physical models of dark energy and modified gravity. We determine the constraints on the model parameters, incorporating the redshift space distortion data directly in the analysis. To determine whether models can be ruled out, we evaluate the $p$ value (the probability under the model of obtaining data as bad or worse than the observed data). In our comparison, we find the well known tension of H$_0$ with the other data; no model resolves this tension successfully. Among the models we consider, the large scale growth of structure data does not affect the modified gravity models as a category particularly differently than dark energy models; it matters for some modified gravity models but not others, and the same is true for dark energy models. We compute predicted observables for each model under current observational constraints, and identify models for which future observational constraints will be particularly informative.

We investigate the possibility of detecting in redshift surveys a hemispherical power asymmetry similar to that first reported in CMB observations. We assume the hemispherical asymmetry arises from a linear gradient in comoving coordinates in the perturbation amplitude. We predict the resulting clustering of galaxy or galaxy cluster tracers using an excursion set approach; doing so accounts for the variation of both the underlying clustering and the tracer bias. Based on the predicted variation of the clustering of tracers, we perform a Fisher matrix forecast of the galaxy clustering amplitude and calculate the statistical significance for ideal surveys and planned surveys. The results indicate that the DESI galaxy survey would be able to detect this signal with higher than $3\sigma$ significance if the asymmetry does exist. We also investigate the amplitude and scale dependence of the above result. The DESI galaxy survey can probe the dipole amplitude higher than 0.04, which correspond to a $\pm4\%$ difference of the temperature fluctuation along and opposite the dipole direction, at least at the $2\sigma$ level. Additionally, we investigate a modulation of the power spectrum that exhibits asymmetry only for large scales. This modulation is potentially detectable. For Milky Way galaxy mass tracers, the scale-dependent modulation yields a larger change in the large scale power spectrum than does a scale-independent modulation, because the former does not alter the bias.

We present the first scientific results from the luminous red galaxy sample (LRG) of the extended Baryon Oscillation Spectroscopic Survey (eBOSS). We measure the small and intermediate scale clustering from a sample of more than 61,000 galaxies in the redshift range $0.6 < z < 0.9$. We interpret these measurements in the framework of the Halo Occupation Distribution. The bias of eBOSS LRGs is $2.30 \pm 0.03$, with a satellite fraction of $13\pm3$\% and a mean halo mass of $2.5\times10^{13}h^{-1}M_{\odot}$. These results are consistent with expectations, demonstrating that eBOSS galaxies will be reliable tracers of large scale structure at $z\sim 0.7$. The eBOSS galaxy bias implies a scatter of luminosity at fixed halo mass, $\sigma_{\log L}$, of 0.19 dex. Using the clustering of massive galaxies from BOSS-CMASS, BOSS-LOWZ, and SDSS, we find that $\sigma_{\log L}=0.19$ is consistent with observations over the full redshift range that these samples cover. The addition of eBOSS to previous surveys allows investigation of the evolution of massive galaxies over the past $\sim 7$ Gyr.

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

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.

Compared with the plentiful researches of the Hubble parameter and deceleration factor, the third time derivative of the scale factor $a(t)$ in the FRW cosmology, namely, the jerk parameter $j$ still lacks attention. In order to study the properties of $j$, we propose several kinds of parameterizations of $j(z)$ as a function of the redshift $z$. By setting the standard $\Lambda$CDM model as the fiducial model, we constrain the jerk models with the observational Hubble parameter data (OHD) and Type Ia Supernovae (SNe) observations. We find that the perturbation of $j(z)$ favors a value of nearly zero and the $\Lambda$CDM is well accommodated by the jerk reconstruction. We also compare the powers of OHD and SNe in constraining the jerk models in detail, and find that the newly released OHD measurement at $z=2.3$ can improve the constraint significantly, even tighter than the SNe one. Furthermore, we analyze the jerk models by calculating the Hubble parameter, equation of state, the deceleration factor and $Om(z)$ diagnostic. Our results show that the universe is indeed undergoing an accelerated expansion phase following the matter-dominated one, which is consistent with the standard model by observations.

We study the holographic and agegraphic dark energy models without interaction using the latest observational Hubble parameter data (OHD), the Union2.1 compilation of type Ia supernovae (SNIa), and the energy conditions. Scenarios of dark energy are distinguished by the cut-off of cosmic age, conformal time, and event horizon. The best-fit value of matter density for the three scenarios almost steadily located at $\Omega_{m0}=0.26$ by the joint constraint. For the agegraphic models, they can be recovered to the standard cosmological model when the constant $c$ which presents the fraction of dark energy approaches to infinity. Absence of upper limit of $c$ by the joint constraint demonstrates the recovery possibility. Using the fitted result, we also reconstruct the current equation of state of dark energy at different scenarios, respectively. Employing the model criteria $\chi^2_{\textrm{min}}/dof$, we find that conformal time model is the worst, but they can not be distinguished clearly. Comparing with the observational constraints, we find that SEC is fulfilled at redshift $0.2 \lesssim z \lesssim 0.3$ with $1\sigma$ confidence level. We also find that NEC gives a meaningful constraint for the event horizon cut-off model, especially compared with OHD only. We note that the energy condition maybe could play an important role in the interacting models because of different degeneracy between $\Omega_m$ and constant $c$.

We present the constraints on the Quintessence scalar field model from the observational data of the variation of the fine structure constant obtained from Keck and VLT telescopes. Within the theoretical frame proposed by Bekenstein, the constraints on the parameters of the Quintessence scalar field model are obtained. By the consideration of the prior of $\Omega_{m0}$ as WMAP 7 suggests, we obtain various results of the different samples. Based on these results, we also calculate the probability density function of the coupling constant $\zeta$. The best-fit values show a consistent relationship between $\zeta$ and the different experimental results. In our work, we test two different potential models, namely, the inverse power law potential and the exponential potential. The results show that both the large value of the parameters in the potential and the strong coupling can cause the variation of fine structure constant.

We investigate the interacting agegraphic dark energy in Brans-Dicke theory and introduce a new series general forms of dark sector coupling. As examples, we select three cases involving a linear interaction form (Model I) and two nonlinear interaction form (Model II and Model III). Our conclusions show that the accelerated scaling attractor solutions do exist in these models. We also find that these interacting agegraphic dark energy modes are consistent with the observational data. The difference in these models is that nonlinear interaction forms give more approached evolution to the standard $\Lambda$CDM model than the linear one. Our work implies that the nonlinear interaction forms should be payed more attention.

The newly released observational $H(z)$ data (OHD) is used to constrain $\Lambda(t)$CDM models as holographic and agegraphic dark energy. By the use of the length scale and time scale as the IR cut-off including Hubble horizon (HH), future event horizon (FEH), age of the universe (AU), and conformal time (CT), we achieve four different $\Lambda(t)$CDM models which can describe the present cosmological acceleration respectively. In order to get a comparison between such $\Lambda(t)$CDM models and standard $\Lambda$CDM model, we use the information criteria (IC), $Om(z)$ diagnostic, and statefinder diagnostic to measure the deviations. Furthermore, by simulating a larger Hubble parameter data sample in the redshift range of $0.1<z<2.0$, we get the improved constraints and more sufficient comparison. We show that OHD is not only able to play almost the same role in constraining cosmological parameters as SNe Ia does but also provides the effective measurement of the deviation of the DE models from standard $\Lambda$CDM model. In the holographic and agegraphic scenarios, the results indicate that the FEH is more preferable than HH scenario. However, both two time scenarios show better approximations to $\Lambda$CDM model than the length scenarios.

We use the Radial Baryon Acoustic Oscillation (RBAO) measurements, distant type Ia supernovae (SNe Ia), the observational $H(z)$ data (OHD) and the Cosmic Microwave Background (CMB) shift parameter data to constrain cosmological parameters of $\Lambda$CDM and XCDM cosmologies and further examine the role of OHD and SNe Ia data in cosmological constraints. We marginalize the likelihood function over $h$ by integrating the probability density $P\propto e^{-\chi^{2}/2}$ to obtain the best fitting results and the confidence regions in the $\Omega_{m}-\Omega_{\Lambda}$ plane.With the combination analysis for both of the {\rm $\Lambda$}CDM and XCDM models, we find that the confidence regions of 68.3%, 95.4% and 99.7% levels using OHD+RBAO+CMB data are in good agreement with that of SNe Ia+RBAO+CMB data which is consistent with the result of Lin et al's work. With more data of OHD, we can probably constrain the cosmological parameters using OHD data instead of SNe Ia data in the future.