Testing and selecting cosmological models with ultra-compact radio quasars

In this paper, we place constraints on four alternative cosmological models under the assumption of the spatial flatness of the Universe: CPL, EDE, GCG and MPC. A new compilation of 120 compact radio quasars observed by very-long-baseline interferometry, which represents a type of new cosmological standard rulers, are used to test these cosmological models. Our results show that the fits on CPL obtained from the quasar sample are well consistent with those obtained from BAO. For other cosmological models considered, quasars provide constraints in agreement with those derived with other standard probes at $1\sigma$ confidence level. Moreover, the results obtained from other statistical methods including Figure of Merit, $Om(z)$ and statefinder diagnostics indicate that: (1) Radio quasar standard ruler could provide better statistical constraints than BAO for all cosmological models considered, which suggests its potential to act as a powerful complementary probe to BAO and galaxy clusters. (2) Turning to $Om(z)$ diagnostics, CPL, GCG and EDE models can not be distinguished from each other at the present epoch. (3) In the framework of statefinder diagnostics, MPC and EDE will deviate from $\rm{\Lambda}$CDM model in the near future, while GCG model cannot be distinguished from $\rm{\Lambda}$CDM model unless much higher precision observations are available.


Introduction
It is strongly indicated that the Universe has entered a stage of accelerated expansion, which was confirmed by a lot of recent observations including Supernova Ia (SN Ia) [1][2][3], baryon acoustic oscillation (BAO) [4], and precise measurements of the spectrum of cosmic microwave background (CMB) [5][6][7][8]. However, it remains a big puzzle in modern cosmology about the origin to the current cosmic acceleration, which gives birth to a large class of cosmological models mathematically explaining this phenomenon. In general, these cosmological scenarios are mainly split into two large categories, the first of which adheres to General Relativity and drives the current accelerated expansion through a dark energy component, while the second focus on gravitational modifications without the inclusion of exotic dark energy.
In the first scenario, the most simple candidate for dark energy is the cosmological constant Λ, in which the equation of state (EoS) of dark energy is equal to −1. This model, the so-called ΛCDM, provides excellent agreements with a wide range of astronomical data so far [9,10]. a e-mail: caoshuo@bnu.edu.cn However, it is confronted with some theoretical problems including the well-known fine tuning problem and coincidence problem [11]. Other models allowing any constant dark energy equation of state (quintessence [12,13], phantom models [14], etc.), as well as models in which the dark energy equation of state is allowed to evolve with time (the well-known CPL parametrization [15,16]) have been extensively studied with various astrophysical probes in the literature [17]. Meanwhile, considering the possible interaction between dark energy and dark matter, the socalled interacting dark energy model [18][19][20] could also contribute to the alleviation of the coincidence problem. Originating from different aspects of new physics, many other dynamical dark energy models such as the Chaplygin gas [21,22] and the holographic dark energy models [23][24][25] have been explored by cosmologists for a long time. In the second scenario, significant interest in modifications to general relativity has also gained a lot attention [26][27][28], with the aim of explaining the acceleration of the Universe without introducing dark energy.
In the face of so many competing cosmological models, many authors turned to various observational probes such as Supernovae Ia acting as standard candles (z ∼ 1.40), strong gravitational lensing systems [29,30], BAO, and CMB (z ∼ 1000) acting as standard rulers, with the aim to determine which one is most favored by the observational data [10,17]. However, in order to achieve this difficult goal of model filtration, it is still necessary to acquire high precision data set and develop new complementary techniques bridging the redshift gap of current data. In the past decades, different types of radio sources have been proposed as possible candidates for standard rulers in cosmological studies [31][32][33]. For instance, the size of the line emitting regions was employed as a standard ruler to study the local expansion history of the Universe, using which Ref. [34] actually derived the distance to Active Galactic Nuclei (AGN). Another new type of useful cosmological ruler is Super-Eddington accreting quasar, the properties of which have been extensively studied in Refs. [35,36]. Recent studies have used the nonlinear relation between the ultraviolet and X-ray luminosity of quasars to place constraints on cosmological parameters [37]. In this paper, we highlight the usefulness of ultracompact structure in radio quasars as a reliable cosmic standard probe to assess some popular cosmological models. In the more recent work [38], a sub-sample of 120 intermediate-luminosity quasars in the redshift range of 0.46 < z < 2.8 was extracted from 613 milliarcsecond ultra-compact radio sources observed by the very-longbaseline interferometry (VLBI) all-sky survey. A pioneer work using this data was also performed in Ref. [39] to probe a flat ΛCDM model and XCDM model, in which the linear size of this standard ruler was calibrated as l m = 11.03 pc through a cosmological-model-independent method. As complementary to other cosmological standard rods, such as BAO and galaxy clusters [40], quasars are promising objects for studying the expansion rate of the Universe at much higher redshift, thus have become an effective probe in cosmology and astrophysics [20,39]. As an extension of the previous work [39], the aim of this analysis is to test alternative cosmological models using the quasar sample and investigate its possibility to provide additional information of model discrimination compared with that provided by other standard ruler data (BAO and galaxy clusters). Two model diagnostics, the Om diagnostic [41] and the statefinder diagnostic [42], are also applied to our work. This paper is organized as follows: In Section 2, we briefly introduce the observational data sets used. The details of the cosmological models to be considered are presented in Section 3. In Section 4, we describe the methods used to obtain the constraints for each data set. In Section 5, we perform a Markov Chain Monte Carlo (MCMC) analysis, and furthermore apply model diagnostics in Section 6. The main results are summarized in Section 7.

Data
Three types of standard rulers currently available are used to place constraint on different cosmological models: the compact radio quasars (QSO) data from VLBI, the angular-diameter distance (D A ) measurements derived from galaxy clusters, and the baryonic acoustic oscillations (BAO) data from large-scale structure (LSS) observations. QSO. In our analysis, we use the angular size measurements of 120 radio quasars covering the redshift range of 0.46 < z < 2.76 [39]. The linear sizes of compact structure in intermediate-luminosity radio quasars (10 27 W/Hz < L < 10 28 W/Hz) display negligible dependent on luminosity and redshift [38]. We refer the reader to [38] for the detailed selection methodology to obtain the final sample of radio quasars explicitly presented in [39], which could serve as standard cosmological rods with intrinsic linear size calibrated to l m = 11.03 ± 0.25 pc.
Galaxy cluster. X-ray observation of intracluster medium and radio observations of Sunyaev-Zeldovich effect allow us to estimate the angular diameter distance of galaxy clusters. In this paper, we will use the D A measurements of 38 galaxy clusters in the redshift range of 0.16 < z < 0.89. The final statistical sample with all necessary information can be found in Refs. [40].
BAO. The third astrophysical probe applied to the joint cosmological analysis is BAO, which measures the angular-diameter distance through the clustering of galaxies perpendicular to the line of sight. The acoustic peak in the galaxy correlation function has been detected over a redshift range of 0.1 < z < 0.7 with large scale surveys. The determination of the BAO scale at lower redshift, z = 0.106, was made in the 6dFGS survey [43], while the other four measurements of acoustic scale at higher redshifts were respectively obtained by SDSS-DR7 [44], SDSS-DR9 [45], and the WiggleZ galaxy survey [46]. These data extensively used in nine-year WMAP analysis were summarized in Table 1 of Ref. [7].
We remark here that, in order to test the cosmological models at higher redshifts, it is very necessary to turn to distance indicators located in the deeper universe. Higherredshift radio-loud quasars are valuable additions to standard rulers used for cosmological tests, since the predic-tions of cosmological models can be radically different. Fig. 1 shows the redshift coverage of different standard ruler data. One can see that inclusion of quasars could result in a fair coverage of redshifts, which enables QSO an excellent complement to other observational probes at lower redshifts.

Cosmological models
Four cosmological models are considered with the data sets described above: Chevallier-Ploarski-Linder parametrization (CPL), Entropy Dark energy model (EDE), Generalized Chaplygin Gas Model (GCG) and modified polytropic Cardassian model (MPC). A flatness of Friedmann Robertson Walker (FRW) metric is assumed in our analysis, which is strongly supported by the recent observations [7,47]. Under this assumption, the angular diameter distance can be expressed as where E(z) = H(z)/H 0 and H 0 is the Hubble constant, which is fixed at 67.8 ± 0.9 km s −1 Mpc −1 based on recent Planck results [8]. Moreover, in order to obtain stringent constraints on other cosmological parameters, for the CPL, EDE, and MPC models, we use the prior on the matter density parameter Ω m from Planck [8]. It should be noted that, although the above priors always influence the cosmological analysis, conclusions should not be significantly affected concerning the core of this paper, i.e., the comparison of the confidence regions derived from different standard ruler data.

CPL model
A simple extension of the ΛCDM model is the XCDM model with constant equation of state. However, it would be natural to consider the equation of state varying with redshifts, i.e., it could be an arbitrary function of the redshift, w = w(z). One of the most popular functions is the CPL parametrization [15,16] where w 0 and w a are the two parameters to be fitted by the observational data. Note that the ΛCDM model can be always recovered by taking w 0 = −1 and w a = 0. Then the Hubble function can be expressed as where Ω DE = 1 − Ω m , and the cosmological parameters in this cosmological model are p = (w 0 , w a ).

EDE model
Recently, the Entropy Dark Energy model was proposed on the base of the theory of entropic gravity [48,49]. In the framework of entropic gravity theory, the gravity force can be explained as a kind of entropic force related to the change of entropy, while the field of equation of gravity is obtained with the second law of thermodynamics. The EDE in the entropic gravity model arises from the surface term in the Einstein-Hilbert's action. In the previous work of Refs. [48,49], a positive term C H H 2 + CḢḢ (the overdot means a derivative with time) was added to the surface part in the action, where C H and CḢ are the model parameters respectively falling into the range of 3/2π ≤ C H ≤ 1 and 0 ≤ CḢ ≤ 3/2π. In our work, no interaction between the DE and other cosmic components (especially matter) is assumed and one can derive the evolution of Hubble parameter as It is straightforward to show that η = Ω m when C H = 3CḢ/2, and the parameters to be considered in this model are p = (C H , CḢ ).

GCG model
The so-called General Chaplygin Gas model (GCG) has been widely studied to explain the accelerating universe [21,22,50]. In the GCG model, the dark sectors in the Universe, i.e., dark energy and dark matter, can be unified through an exotic equation of state. More specifically, the GCG background fluid with its energy density ρ GCG and pressure p GCG can be related with the equation of state [21] is the unified energy density of dark energy and dark matter. The Universe is filled with two components, the GCG component and the baryonic matter component, i.e., ρ = ρ GCG + ρ b . Under the assumption of flat FRW metric, the Hubble parameter of this model can be expressed as where Ω GCG = 1 − Ω b and A s = A ρ 1+α . In our analysis, the baryonic density parameter is fixed at Ω b = 0.0484 based on the recent Planck results [8], and the two cosmological parameters in this model are p = (A s , α).

MPC model
In order to explain the cosmic acceleration from a different perspective, Freese & Lewis (2002) put forward a Cardassian model without the introduction of dark energy [26], for which the Friedmann equation is modified as where ρ m is the total matter density. We emphasize here that, in order to lead to the cosmic acceleration in this parameterization, the value of the parameter n should always be n < 2/3. Then a simple generalized case of the Cardassian model was proposed in Ref. [51], in which an additional exponent q was introduced. The Hubble parameter with this generalization can be written as where the parameters to be constrained are p = (q, n). This MPC model will reduce to ΛCDM model with q = 1 and n = 0.

Method
In the following, we consider the observational constraints on the cosmological models from observational data.

QSO
If taking milliarcsecond structure in radio quasars as individual standard rulers, the angular sizes at redshift z can be written as where D A (z) is the corresponding angular diameter distance mentioned above and l m is the intrinsic length of milliarcsecond structure in radio quasars. In this work, we take the typical value of l m = 11.03 ± 0.25 pc calibrated with cosmic chronometers, which was obtained through a new cosmology-independent calibration technique [39]. Following the general classification of Active Galactic Nuclei (AGN), it is powered by the accretion of mass onto black holes in the center of galaxies and will produce jets of relativistic plasma in the central regions. There are two main physical meaning related to the linear size of this standard ruler: on the one hand, there is almost no stellar contribution when the distance from the AGN center approaches 10 pc, which is also the position at which AGN jets are typically generated [52]; On the other hand, according to the recent analysis of the correlation between the black hole's mass accretion and the star-formation rate, 10 pc represents the typical radius within which the two rates are almost equal to each other, a conclusion supported by the findings from both recent observations and simulations of AGN [53,54]. More importantly, the value of l m estimated from single-frequency VLBI measurements agrees very well with that obtained from multifrequency VLBI imaging observations [55]. Such consistency could also be seen from the comparison between the cosmological fits derived from two types of VLBI observations at different observing frequencies [56]. The data points of the 120 QSOs are given in terms of the angular sizes, θ obs i . One can then constrain cosmological parameters by minimizing the χ 2 function given by where θ(z i ; p) is the theoretical value of the angular size at redshift z (which is defined in Eq. (9)) and p represents the cosmological parameters of interest (which are specifically introduced in Section 3). θ obs i is the observed counterpart of the angular size for the ith quasar. Following the work of Ref. [39], in our analysis the total uncertainty expresses as σ 2 i = σ 2 sta,i + σ 2 sys,i , which includes the statistical error of observations in θ obs i and an additional 10% systematical uncertainty accounting for the intrinsic spread in the linear size. See Table 1 of Ref. [39] for details of the quasar data and reference to the source papers.

Galaxy clusters
We can obtain the angular diameter distances by using the Sunyaev-Zeldovich effect together with x-ray emission of galaxy clusters, which can be directly used to estimate the cosmological parameters by minimizing the corresponding χ 2 as Here D th A (z i ; p) is the theoretical angular diameter distance at redshift z i , which is defined in Eq. (1). D obs A (z i ) is the observed angular diameter distance of the ith galaxy cluster with total uncertainty defined as σ 2 DA,i = σ 2 mod + σ 2 stat + σ 2 sys , where the modeling error (σ mod ), statistical error (σ stat ) and systematical error (σ sys ) are explicitly shown in Table 2-3 in Ref. [40].

BAO
Compared with the previous works involving BAO as a standard ruler [9,57,58], we use the measurement of distance ratio r s (z d )/D V (z) or D V (z)/r s from the BAO peak in the distribution of SDSS luminous red galaxies, which contains the main information of the observations of LSS. Here r s (z d ) is the comoving sound horizon at the drag epoch, where the redshift z d at the baryonic drag epoch is fitted with the formula proposed in Ref. [59]. D V (z) is the effective distance given by where D A is the angular diameter distance and H(z) is Hubble parameter. The χ 2 function for BAO is defined as is the corresponding inverse covariance matrix [7].

Joint analysis
We will present a combined analysis of the above three tests to fit theoretical models to observational data. Meanwhile, in order to make a good comparison with the quasar sample, a joint analysis with galaxy clusters and BAO data sets is also performed in this analysis. The χ 2 function of the above two combined analysis are respectively expressed as and

Results and discussion
For all cosmological models described in Section 3, we estimate the constraint ability of different angular diameter distance data, QSO, BAO, galaxy clusters, BAO+Cluster and BAO+Cluster+QSO, by minimizing the χ 2 function given in Sect. 4. Furthermore, the Figure of Merit (FoM) [60] is also applied to quantify the constraining power of each data (especially the quasar sample) on cosmological model parameters.

CPL
The best fits for CPL parameters w 0 , w a and the estimated χ 2 from different data sets are shown in Table  1. The 1σ, 2σ contours of the model parameters are presented in Fig.2. As can be seen from Table 1 and Fig.2, the fitting results from QSO are in good agreement with those from BAO, whereas in tension with the results from galaxy cluster data. Notice that concordance ΛCDM cosmology (w 0 = −1, w a = 0) is consistent with the quasar and BAO standard ruler data at less than 1σ level. More importantly, compared with the previous literature using other independent and precise experiments [30,58,61], the currently compiled quasar data may improve the constraints on model parameters significantly, in the framework of CPL parametrization. When adding QSO data set to the joint data set of BAO and galaxy cluster, we will get more precise assessment of w 0 , w a , which is consistent with that obtained from the recent Planck CMB data [8] as well as the combination of the CMB measurements from Atacama Cosmology Telescope (ACT) and the South Pole Telescope, BAO and H 0 measurements [7]. Table 1 shows the best fits of EDE parameters (CḢ , C H ) and the χ 2 derived from different observational data sets. Although the exact value of CḢ are not independently obtained with QSO or BAO, appreciable consistency between the same type of probes (standard rulers) is indeed revealed in our analysis. On the other hand, it is clear that the quasar data set could provide constraints on the other parameter C H comparable to the other two types of standard probes. These implications can be more clearly seen from the corresponding contours for EDE model, which are explicitly presented in Fig. 3. Fitting results from the joint angular diameter distance data of QSO+BAO+Cluster give the best-fit parameters CḢ > 0.469 and C H = 0.895 +0.040 −0.027 , which agree very well with the results yielded from the luminosity distance data in-

GCG
Working on the GCG model, we obtain the fitting results from different combinations of observational data, which are displayed in Table. 1 and Fig. 4. The excellent consistency between the three types of standard rulers can be clearly seen through the comparison of these plots. The joint analysis with standard rulers provides the best-fit parameters and the marginalized 1σ constraints as A s = 0.708 ± 0.040 and α = −0.05 +0.10 −0.16 . For comparison, it is necessary to refer to the previous results obtained with other astrophysical measurements. The results obtained with the combination analysis of the X-ray mass fractions of galaxy clusters, the dimensionless coordinate distance to SN I and FRIIb radio galaxies gave A s = 0.70 +0. 16 −0.17 and α = −0.09 +0. 54 −0.33 [63]. The former work done in Ref. [64] with the joint data of 9 Hubble parameters data points, 115 SN Ia and BAO peak at z = 0.35 showed that 0.67 ≤ A s ≤ 0.83 and −0.21 ≤ α ≤ 0.42. Recent work done in [50] indicated A s = 0.73 ± 0.06 and α = −0.09 +0. 15 −0.15 , which strengthens the indication that joint analysis of cosmic standard rulers (QSO+BAO+Cluster) could provide consistent but more stringent fitting results compared with these previous results.

MPC
All of the fitting results obtained with QSO, BAO, Cluster and the joint data are presented in Table 1. The 1σ, 2σ contours of the MPC model parameters β and n are also illustrated in Fig. 5. The results from QSO, BAO and galaxy cluster are consistent with each other within 1σ confidence level. Several authors have tested the MPC model using various data sets. For instance, Cosmic All-Sky Survey (CLASS) lensing sample [65] gave β = 0.05, n = −2.32, which is in tension with our results from QSO+BAO+Cluster. However, It is important to note that the shape of the constraint contours derived in our analysis are very similar to those shown in Refs. [65][66][67][68][69]. Moreover, the results obtained with standard rulers turned out to correspond well with previous works. Our results are similar to the results obtained with SN Ia+BAO+WMAP5 (β = 0.48 +2.020 −0.080 , n = −0.600 +0.980 −0.450 ) and SN Ia+BAO+CMB data sets (β = 1.098 +1.015 −0.465 , n = −0.041 +0.364 −0.964 ) at 1σ confidence level [66]. More recent works [69] have suggested the cutoffs of 0.45 < β < 1.05 and −0.8 < n < 0.05, which achieved a similar precision to our work.

Figure of Merit
As discussed above, the method of using the angular size measurements of compact structure in radio quasars distance could provide a complementary and effective probe in cosmological applications. However, in order to quantify the constraining power of the quasar sample, we introduce the Figure of Merit (FoM) [60,70], a useful statistical tool originally defined by the Dark Energy Task Force as the inverse of the area enclosed by the 95% confidence level contour of CPL parameters, (w 0 , w a ). Later there was a more general definition of FoM [71] where Cov(f 1 , f 2 , f 3 , ...) is the covariance matrix of the cosmological parameters f i . Note that a larger FoM corresponds to a smaller error ellipse, which therefore denotes a tighter constraints on the cosmological parameters.
We have calculated the FoM of the cosmological models for each data set analyzed, which is explicitly summarized in Table 1. A graphical representation of the FoM results is also provided in Fig. 6, which directly shows the results in the FoM test for each cosmological model. Out of all the candidate models considered, it is obvious that the QSO data could provide better statistically constraints on cosmological parameters than BAO. This could attribute to the large sample size and the higher redshift range covered by QSOs in comparison to other cosmological probes. When taking galaxy cluster observations into consideration, quasars perform better than galaxy clusters in the framework of two cosmological models, GCG and EDE models. On the other hand, when comparing the FoM of BAO+Cluster and that of all observations, we find the inclusion of the QSO sample will generate more stringent cosmological constraints.

Model diagnostics
In order to discriminate the four cosmological scenarios considered in this paper, it is important to find sensitive and robust diagnostics to illustrate the dynamic behavior of different cosmologies. As is well known, the expansion rate of the Universe can be expressed by the Hubble parameter H =ȧ/a, where a is the scale factor, while the rate of cosmic acceleration ia always quantified by the deceleration parameter However, it is very difficult for the Hubble parameter H and the deceleration parameter q to accurately distinguish cosmological models cause all the models considered will give similar results, e.g.,ä > 0 and H > 0 or q < 0, which encourages us to invoke some newer and more effective quantities to substitute the two original parameters. In this work, we will take the Om diagnostic and the statefinder diagnostic into consideration. It is well known that Om(z) is a combination of the Hubble parameter and the redshift, which provides a null test of dark energy being a cosmological constant at different stages for the ΛCDM model [41]. Therefore, this diagnostic, which has been extensively used to discriminate different cosmological models as well as ΛCDM model [72,73], can be defined as where E(z) = H(z)/H 0 . In the basic ΛCDM model neglecting the radiation at low redshifts, one can easily get The combination of Eq. (19) and Eq. (20) will lead to It is obvious that Om(z) should be constant and exactly equal to the present mass density parameter Ω m if the ΛCDM model is the true one, while for other cosmological models, the Om(z) diagnostic are expected to give different values. Applying the Om(z) diagnostic to the cosmological models considered in our work, we can get the relation between the redshift and Om(z) for different cosmological models, which is specifically presented in Fig. 7. The Om(z) for the CPL and GCG models cannot be distinguished from each other as well as from ΛCDM model. Moreover, the Om(z) of CPL, GCG, EDE models cannot be distinguished at present time unless high-precision observations are obtained and applied. Another impressive feature of Fig.7 is that the Om(z) for the MPC model, a cosmological candidate proposed without introducing dark energy in the Universe, absolutely deviates from the ΛCDM model and other cosmological models.
Apart from the Om(z), the statefinder diagnostic, which has been extensively applied to discriminate different cosmological models, involves the third derivative of the scale factor a as [42] and one can plot the corresponding trajectories in the r−s plane. For a certain cosmological model, the statefinder can be easily derived as and where E(z) = H(z)/H 0 and the deceleration parameter q(z) can be expressed as Applying the best fits from the joint analysis to each cosmological model, we obtain the evolution of the statefinder (r, s) and the deceleration parameter q. The evolution of the statefinder pair (r, s) for different cosmological models is shown in Fig. 8. The red point at (r, s) = (0, 1) indicates the statefinder of ΛCDM model and the diamond on each curve shows the present value of the statefinder (r, s) for each cosmological model. It is apparent that the CPL model can be distinguished from other cosmological models at present time, however, it will approach ΛCDM in the near future. Meanwhile, the MPC and EDE models, which are not distinguishable from each other by the statefinder, is deviating from ΛCDM at the present epoch. More importantly, the GCG model exhibits very similar evolution tendency to the concordance cosmological constant.
The evolution trajectories in the r−q plane are plotted in Fig. 9. For the four cosmological models considered in this analysis, we observe the signature flip from positive to negative in the value of q, which successfully explains the recent phase transition of these models. The diamond points on different curves in Fig. 9 denote the value of q and r at present time for different cosmological models. One can see the value of deceleration parameter q is very close to each other at present time, which is quite different from the behavior of r. As for the evolution of cosmological models, at the present epoch, the GCG model and ΛCDM model are not distinguishable and the MPC model can not be distinguished from EDE model. However, in the near future they will evolve diversely, which is in well consistent with the results obtained from the r − s plot.

Conclusions
In this paper, we place constraints on four alternative cosmological models under the assumption of the spatial flatness of the Universe: Chevallier-Ploarski-Linder A new compilation of 120 angular-size/redshift data compact radio quasars observed by very-long-baseline interferometry (VLBI), whose statistical linear sizes show negligible dependence on redshifts and intrinsic luminosity and thus represents standard rulers in cosmology, are used to test these cosmological models. Compared with BAO and galaxy clusters acting as cosmological standard rulers, higher-redshift radioloud quasars are valuable additions to standard rulers used for cosmological tests and the inclusion of quasars could result in a fair coverage of redshifts, which enable QSO to be an excellent complement to other observational probes at lower redshifts.
Our results show that the constraints on CPL obtained from the quasar sample are well consistent with that obtained from BAO but in tension with that from galaxy clusters. Note that the concordance ΛCDM cosmology (w 0 = −1, w a = 0) is consistent with the quasar and BAO standard ruler data at less than 1σ level. For other cosmological models considered, quasars provide fits in agreement with those obtained with other probes at 1σ level. Meanwhile, we have calculated the Figure of Merit for each cosmological model, which is explicitly summarized in Table 1. Out of all the candidate models considered, it is obvious that the QSO data could provide better statistically constraints on cosmological parameters than BAO. When taking the observations of galaxy cluster into consideration, quasars perform better than galaxy clusters in the framework of two cosmological models, GCG and EDE models. Based on the best-fits obtained with QSO+BAO+Cluster, we apply two model diagnostics, Om(z) and statefinder to differentiate the dynamical behavior of the four cosmological models. On the one hand, the results from the Om(z) diagnostic show that the CPL, GCG, EDE models cannot be distinguished at the present epoch. However, the MPC model, a cosmological candidate proposed without introducing dark energy in the Universe, absolutely deviates from the ΛCDM model and other cosmological models. On the other hand, in the framework of statefinder diagnostics, MPC and EDE are will deviate from ΛCDM model in the near future, while GCG model cannot be distinguished from ΛCDM model unless much higher precision observations are available.