New observational constraints on $f(T)$ cosmology from radio quasars

Using a new recently compiled milliarcsecond compact radio data set of 120 intermediate-luminosity quasars in the redshift range $0.46<z<2.76$, whose statistical linear sizes show negligible dependence on redshifts and intrinsic luminosity and thus represent standard rulers in cosmology, we constrain three viable and most popular $f(T)$ gravity models, where $T$ is the torsion scalar in teleparallel gravity. Our analysis reveals that constraining power of the quasars data (N=120) is comparable to the Union2.1 SN Ia data (N=580) for all three $f(T)$ models. Together with other standard ruler probes such as Cosmic Microwave Background and Baryon Acoustic Oscillation distance measurements, the present value of the matter density parameter $\Omega_m$ obtained by quasars is much lager than that derived from other observations. For two of the models considered ($f_1$CDM and $f_2$CDM) a small but noticeable deviation from $\Lambda$CDM cosmology is present, while in the framework of $f_3$CDM the effective equation of state may cross the phantom divide line at lower redshifts. These results indicate that intermediate-luminosity quasars could provide an effective observational probe comparable to SN Ia at much higher redsifts, and $f(T)$ gravity is a reasonable candidate for the modified gravity theory.


Introduction
The current cosmic acceleration has been supported by many independent astrophysical observations, including type Ia supernovae (SN Ia) [1], large scale structure [2], cosmic microwave background (CMB) anisotropy [3], etc. A mysterious component with negative pressure, dubbed as dark energy, has been proposed to explain this phenomenon in the framework of Einstein's general relativity, which gave birth a corresponding author e-mail: caoshuo@bnu.edu.cn to a large number of dark energy models including the cosmological constant (Λ CDM), scalar field theory [4][5][6], and dynamical dark energy models [7][8][9][10][11]. The other direction one could follow in search for solution of the accelerating cosmic expansion enigma is to construct modified theories of gravity instead of invoking exotic dark energy. Large majority of works in this direction concentrated on the braneworld Dvali-Gabadadze-Porrati (DGP) model [12], f (R) gravity [13], and Gauss-Bonnet gravity [14].
Equally well, one can also modify the gravity according to the scenario described by the so-called f (T ) theory [15], which was proposedin the framework of the Teleparallel Equivalent of General Relativity (also known as Teleparallel Gravity). In this approach, the Levi-Civita connection used in Einstein's general relativity is replaced by the Weitzenböck connection with torsion, while the Lagrangian density of this theory is the torsion scalar T . Compared with the f (R) theory leading to the fourth order equations, the field equations of the f (T ) theory are in the form of second order differential equations, which provides an important advantage of this approach. In addition, if certain conditions are satisfied, the behavior of f (T ) cosmologies is similar to several popular dark energy models, such as quintessence [16], phantom [17], DGP model [18] and transient acceleration [19]. Due to the above mentioned property, f (T ) theory and its cosmological applications has gained a lot interest in the literature. The detailed introduction to the f (T ) theory could be found in [18,20].
In this paper, we focus on using the currently released quasar data [21] to provide the constraints on various f (T ) gravity models. Recently, the angular size of compact structure in radio quasars versus redshift data from the very-longbaseline interferometry (VLBI) observations have become an effective probe in cosmology. Reliable standard rulers and standard candles at cosmological scales are crucial for measuring cosmic distances at different redshifts. For in-stance, the type Ia supernovae are regarded as standard candles, while the BAO peak location is commonly recognized as a fixed comoving ruler. The increasing observational material concerning these two distance indicators has been widely used in various cosmological studies. In the past, there were controversial discussions about whether the compact radio sources could act as standard rulers [22][23][24][25][26]. The difficulty lies in the fact that the linear sizes l m of compact radio sources might not be constant, i.e., its value is dependent on both redshifts and some intrinsic properties of the source (luminosity, for example). Based on a 2.29 GHz VLBI all-sky survey of 613 milliarcsecond ultra-compact radio sources [27,28], Cao et al. [21] presented a method to divide the full sample into different sub-samples, according to their optical counterparts and luminosity (low-luminosity quasars, intermediate-luminosity quasars, and high-luminosity quasars). The final results indicated that intermediate-luminosity quasars show negligible dependence on both redshifts z and intrinsic luminosity L, which makes them a fixed comoving-length standard ruler. More recently, based on a cosmological-modelindependent method to calibrate the linear sizes l m of intermediateluminosity quasars, Cao et al. [29] investigated the cosmological application of this data set and obtained stringent constraints on both the matter density Ω m and the Hubble constant H 0 , which agree very with the recent Planck results. The advantage of this data set, compared with other standard rulers: BAO [30][31][32], clusters [33], strong lensing systems [34][35][36]), is that quasars are observed at much higher redshifts (z ∼ 3.0). Therefore, it may be rewarding to test the f (T ) theory with this newly revised quasar data. In this paper, we examine constraints on the viable f (T ) cosmological models imposed by the quasars. We compare them with analogous results obtained with the newly revised Union2.1 set -the largest published and spectroscopically confirmed SN Ia sample to date. We expect that different systematics and sensitivities of these two different probes (rulers vs. candles) can give complementary results on the f (T ) theory. This paper is organized as follows: In section 2 we briefly introduce the f (T ) gravity and its cosmological consequences. In section 3 we present the latest data sets for our analysis and perform a Markov chain Monte Carlo analysis using different data sets. Finally, we summarize the main conclusions in Section 4.

The f (T ) theory
In this section we brief review the f (T ) gravity in the framework of cosmology, and then present three specific f (T ) models to be analyzed in this work.

The f (T ) cosmology
We use the vierbein fields e i (x µ ) (i = 0, 1, 2, 3), which is an orthonormal basis for the tangent space at each point x µ of the manifold e i · e j = η i j , and whose components are e µ i (µ = 0, 1, 2, 3) (here Latin indices stand for the tangent space and Greek indices refer to the manifold). Its dual vierbein gives the metric tensor g µν (x) = η i j e i µ (x) e j ν (x). In f (T ) theory, instead of the torsionless Levi-Civita connection in Einstein's General Relativity, the curvatureless Weitzenböck connection is considered, and hence the torsion tensor describing the gravitational field is The Lagrangian of teleparallel gravity is constructed by the torsion scalar as [15] T where and the contorsion tensor K µν ρ is given by In the f (T ) theory, the Lagrangian density is a function of T [15], and the action reads where e = det(e i µ ) = √ −g. The corresponding field equation is where µν ρ , and T µν is the matter energy-momentum tensor. Considering a flat homogeneous and isotropic FRW universe, we have where a(t) is the cosmological scale factor. By substituting Eqs. (7), (1), (3) and (4) into Eq. (2), one could obtain the torsion scalar as [15] T ≡ S ρ µν T ρ µν = −6H 2 , where H is the Hubble parameter H =ȧ/a. The dot represents the first derivative with respect to the cosmic time. Substituting Eq. (7) into (6), one can obtain the corresponding Friedmann equations where ρ and p are the total energy density and pressure, respectively. By defining the effective energy density density ρ eff , pressure p eff and effective equation of state (EoS) parameter w eff as The Friedmann equations could be rewrite as 3 1 Therefore, the cosmic acceleration could be driven by the torsion instead of dark energy. In this cosmological framework, the corresponding normalized Hubble parameter is where T 0 = −6H 2 0 (the subscript "0" denotes the current value). Here we consider the matter and radiation in the Universe -the components whose energy density evolves with redshift z as ρ m = ρ m0 (1 + z) 3 , ρ r = ρ r0 (1 + z) 4 , respectively. And then, Eq. (16) could be expressed as [37,38] where Ω F = 1 − Ω m − Ω r, and Ω i = k 2 ρ i0 3H 2 0 . In this way, a specific form of f (T ) is embodied in the function y(z, p), whose expression is where p stands for the parameters in different forms of f (T ) theory.

Specific f (T ) models
In this subsection we briefly review three specific f (T ) models, which have passed basic observational tests [37] and will be further investigated in this paper.
(1) The power-law model [15] (hereafter f 1 CDM) assumes that the Lagrangian density f (T ) of the theory is the following: where α and b are two model parameters. The distortion parameter b quantifies deviation from the Λ CDM model, whereas the parameter α can be expressed through the Hubble constant and density parameter Ω F0 by inserting Eq. (19) into Eq. (17) with the boundary condition E(z = 0) = 1 : Now Eq. (18) may be rewritten as Depending on the choice of parameter b, this f (T ) model can be connected with some popular dark energy models. For b = 0, it reduces to the Λ CDM, while it can mimic the Dvali-Gabadadze-Porrati (DGP) model when b = 1/2.
(2) The exponential model [39] (hereafter f 2 CDM) is characterized by where α and p are two dimensionless parameters. Similarly the expressions for α and y(z, p) can also be obtained as This model reduces to the Λ CDM in the limit p → +∞. By setting b = 1/p, Eq. (24) is rewritten as and Λ CDM is recovered when b → 0 + .
(3) Motivated by the exponential f (R) gravity, the hyperbolictangent model [17] (hereafter f 3 CDM) arises from the ansatz where α and n are the two model parameters. We obtain the expressions for α and y(z, p) as respectively. Compared with two previous f (T ) theories, this f (T ) model cannot be reduced to the Λ CDM for any value of its parameters. In addition, in order to have a positive value for ρ e f f , the parameter n must be greater than 3/2 [17].

Observational data and fitting method
In order to measure the angular diameter distance, we always turn to objects of known comoving size acting as "standard rulers". In this paper, we will consider a combination of three types of standard rulers using the most recent and significantly improved observations, i.e., the compact radio quasars data from VLBI, baryonic acoustic oscillations (BAO) from the large-scale structure, and the cosmic microwave background (CMB) measurements.

Quasars data
It is well known that the baryon acoustic oscillations (BAO) peak location is commonly recognized as a fixed comoving ruler of about 100 Mpc. Therefore it has already been used in cosmological studies [30][31][32]. In the similar spirit, as extensively discussed in the literature, compact radio sources (quasars, in particular) constitute another possible class of standard rulers of about 10 pc comoving length. Following the analysis of Gurvits [28], luminosity and redshift dependence of the linear sizes of quasars can be parametrized as where β and n are two parameters quantifying the "angular size -redshift" and "angular size -luminosity" relations, respectively. The parameter l is the linear size scaling factor representing the apparent distribution of radio brightness within the core. The data used in this paper were derived from an old 2.29 GHz VLBI survey undertaken by Preston et al. (1985), which contains 613 milliarcsecond ultracompact radio sources covering the redshift range 0.0035 < z < 3.787. More recently, Cao et al. [29] presented a method to identify a sub-sample which can serve as a certain class of individual standard rulers in the Universe. According to the optical counterparts and luminosities, the full sample could be divided into three sub-samples: low-luminosity quasars (L < 10 27 W/Hz), intermediate-luminosity quasars (10 27 W/Hz < L < 10 28 W/Hz) and high-luminosity quasars (L > 10 28 W/Hz).
The final results showed that only intermediate-luminosity quasars show negligible dependence (|n| ≃ 10 −3 , β ≃ 10 −4 ), and thus they could be a population of rulers once the characteristic length l is fixed. In our analysis, we will use the observations of 120 intermediate-luminosity quasars covering the redshift range 0.46 < z < 2.80, while the linear size of this standard ruler is calibrated to l = 11.03 pc through a new cosmology-independent technique [29]. The observable quantity in this data-set is the angular size of the compact structure in intermediate-luminosity radio quasars, whose theoretical (i.e. determined by the cosmological model) counterpart is where D A is the angular diameter distance at redshift z and the f (T ) model parameters p directly enter the angular diameter distance through where E(z ′ ; p) is the dimensionless Hubble parameter and h is the dimensionless Hubble constant. We estimate the f (T ) parameters by minimizing the corresponding χ 2 defined as where θ obs (z i ) is the observed value of the angular size and σ θ (z i ) is the corresponding uncertainty for the ith data point in the sample. In order to properly account for the intrinsic spread in linear sizes and systematics we have added in quadrature 10% uncertainties to the σ θ (z i ).

CMB and BAO data
In order to diminish the degeneracy between f (T ) model parameters we also used the accurate measurements of BAO and CMB.
The CMB experiments measure the temperature and polarization anisotropy of the cosmic radiation in the early epoch. In general, they are a very important tool for the inference of cosmological model parameters. In particular, the shift parameter R defined as: where z * = 1090.43 denotes the decoupling redshift, is a convenient quantity for a quick fitting of cosmological model parameters. The first-year data release of Planck reported its value of R = 1.7499 ± 0.0088 [41]. We estimate the model parameters by minimizing the corresponding χ 2 The measurements of Baryon acoustic oscillation (BAO) in the large-scale structure power spectrum and CMB angular power spectrum have also been widely used for cosmological applications. In this work we consider the measure- is the co-moving angular-diameter distance, and the dilation scale is given by The BAO data are shown in Table 1. Similarly, the corresponding χ 2 for the BAO probes is defined as i j is the inverse covariance matrix given by Ref. [40].  Table 1 Ratios of distances and the so called dilation scale D V (z BAO ) at different redshifts z BAO taken after [30][31][32] and [40].

Observational constraints
In this section, we determine the model parameters of three f (T ) cosmologies through the maximum likelihood method based on χ 2 introduced in previous section using the Markov Chain Monte Carlo (MCMC) method. Our code is based on CosmoMC [42] and we generated eight chains after setting R − 1 = 0.001 to guarantee the accuracy of this work. 4 In the case of the f (T ) theory based on f (T ) = α(−T ) b , different data sets and their combinations led to the marginalized 2D confidence contours presented in Fig. 1-2. The corresponding marginal 1σ error bars can also be seen in Table  2.
Left panel of the Fig. 1 shows the contours obtained from the quasars only and in combination with CMB and BAO. We remark that the quasar data only can not tightly constrain the model parameters. In order to clearly illustrate the constraint comparison between different data sets, a prior b > −1 is applied to the likelihood contours obtained from the quasar data. Quantitatively, the value of the distortion parameter b, which quantifies the deviation from the Λ CDM model varies over the interval [-3, 0.56] within 1σ confidence level. As it is well known, the main evidence for cosmic acceleration came from the other type of distance indicators in cosmology, those probing the luminosity distance, D L by observing the flux of type Ia supernovae (SN Ia). In order to compare our fits with the results obtained using SN Ia, likelihood contours obtained with the latest Union2.1 compilation [43] consisting of 580 SN Ia data points are also plotted in the right panel of the Fig. 1. It is clear that the quasar data could give more stringent constraints than SN Ia, and its constraining power becomes obvious when the large size difference between the samples is taken into consideration. This may happen due to the wider redshift range of the quasars data (0.46 < z < 2.8) compared with SN Ia (0.015 ≤ z ≤ 1.41). Moreover, one can clearly see from Fig. 1 that principal axes of confidence regions obtained with SN and quasars are inclined at higher angles, which sustains the hope that careful choice of the quasar sample would eventually provide a complementary probe breaking the degeneracy in the f (T ) model parameters. Fi-nally, our method based on the observations of intermediateluminosity quasars may also contribute to testing the consistency between luminosity and angular diameter distances [44][45][46].
With the combined standard ruler data sets of quasars, BAO and CMB, the best-fit value for the parameters are Ω m = 0.321 ± 0.012 and b = 0.080 ± 0.077 within 68.3% confidence level. For comparison, fitting results from SN+BAO+CMB are also given in Fig. 1. The best-fit value is Ω m = 0.297 +0.015 −0.017 and b = −0.12 +0. 17 −0.13 , which is in good agreement with that of the Quasar+BAO+CMB data. It is obvious that the quasar data, when combined to CMB and BAO observations, can give more stringent constraints on this f (T ) cosmology, which demonstrates the strong constraining power of BAO and CMB on the cosmological parameters. This situation has also been extensively discussed in the previous works investigating dark energy scenarios with other astrophysical observations [46][47][48][49][50]. Again, the constraining power of 120 quasar data is comparable to that of 580 SN Ia. On the one hand, the present value of the matter density parameter Ω m given by quasars is much lager than that derived from other observations. This has been noted by our previous analysis Cao et al. [29] and the first-year Planck results, in the framework of Λ CDM cosmology. Such a result indicates that quasars data at high redshifts may provide us a different understanding of the parameters describing the components of the Universe. On the other hand, the parameter b, which captures the deviation of f (T ) cosmology from the Λ CDM scenario, seems to be vanishing or slightly larger than 0 with the combined Quasar+BAO+CMB data. It is interesting to note that Λ CDM is not included at 1σ confidence level (b = 0.08 ± 0.077), this slight deviation from Λ CDM is also consistent with a similar conclusion obtained in Ref. [37] for this f 1 CDM model. This tendency can be more clearly seen from Fig. 3, which illustrates the comparison between the effective equation of state for f (T ) and the EoS for Λ CDM model at z ∼ 4, with the best-fitted value as well as the 1σ and 2σ uncertainties derived from the joint data of Quasars, BAO and CMB.
The contours constrained with the total combination of Quasars+SN Ia+BAO+CMB are presented in Fig. 2, and the best-fit value is Ω m = 0.317 ± 0.010 and b = 0.057 +0.091 −0.065 . The combined data give no stronger constraint, which indicates the constraint ability of quasars data is already very strong, while SN Ia do not play a leading role in the joint   Table 2 we also list alternative constraints obtained by the others using different probes.
Performing a similar analysis as before, this time with the other f (T ) model in which Λ CDM is also nested, namely, , we made the same comparison as f 1 CDM discussed above, i.e. Quasars vs. SN Ia and Quasars+BAO+CMB vs. SN Ia+BAO+CMB. The results are presented in Fig. 4 and the estimated cosmic parameters are briefly summarized in Table 3. It is apparent that the quasars data exhibit similar constraining power as in the case of f 1 CDM model, which implies that the constraint ability of 120 quasar data can be comparable to that of 580 SN Ia. By fitting the f 2 CDM model to Quasars+BAO+CMB, we obtain Ω m = 0.319 ± 0.011 and b < 0.268 (let us recall that here we introduced b = 1/p).
With the combined data set of Quasars+SN Ia+BAO+CMB, we also get the marginalized 1σ constraints of the parameters as Ω m = 0.319 +0.010 −0.011 and b < 0.224. The marginalized 1σ and 2σ contours of each parameter are presented in Fig. 6. In Table 3, the best-fit parameters and their 1σ uncertainties for three data sets are displayed. As previously the results from the others using different probes are shown for comparison. Obviously, the present matter density parameter Ω m fitted by quasars is lager than given by other ob-    servations. The parameter b quantifying the deviation from the Λ CDM scenario, tends to be zero for all of observations listed in Table 3, which results in that the exponential gravity is practically undistinguishable from Λ CDM. As can be seen from the results presented in Fig. 5, even at 2σ confidence level, the effective EoS of f 2 CDM model from the joint analysis of Quasars, BAO and CMB agrees very well with that of Λ CDM at z ∼ 4, which strongly indicates the consistency between the two types of cosmological models at much higher redshifts. Data Ref.   Now we will discuss the third f (T ) cosmology which is truly an alternative to the Λ CDM since the concordance cosmological model cannot be recovered as a limiting case of f 3 CDM model. Consequently, the parameter n does not characterize the deviation from Λ CDM. In Fig. 7 we presented contour plots of f 3 CDM model parameters fitted to four different probes, namely Quasars, SN Ia, Quasars+BAO+CMB, and SN Ia+BAO+CMB. As we can see, the quasar data provide more stringent constraints than SN Ia, which indicates that the constraining ability of quasar data can be comparable to or better than that of SN Ia at least in this particular model. In Fig. 8 we show the contour plots for the combination of all data sets Quasars+SN Ia+BAO+CMB. Additionally, in Table 4 we summarize the best-fit values for the three combined data sets respectively. The table 4 also includes the best-fit values and their 68% confidence levels for the previous results from the literature. Similar to the cases of f 1 CDM model and f 2 CDM model, the present matter density parameter Ω m implied by quasars is lager than that given by other observations. Concerning the value of the parameter n, its the value constrained by all of the current observations satisfies the condition n > 3/2, which is necessary to achieve the cosmic acceleration in the framework of f 3 CDM.
In Fig. 9 we show the evolution of the effective equation of state for f 3 CDM model as a function of redshift, concerning the best-fitted value with the 1σ and 2σ uncertainties from the joint data of Quasars, BAO and CMB. In particular, we find that the value of n obtained with quasars suggests that the effective equation of state crosses the phantom divide line at lower redshifts [17].

Model selection
In order to to make a good comparison between different models or decide which model is preferred by the observational data, we will use two standard information criteria, namely the Akaike Information Criterion (AIC) [53] and the Bayesian Information Criterion (BIC) [54] to study competing models. The above two information criteria are respectively defined as and where L = exp(−χ 2 min /2), k represents the number of free parameters in the model and N is the sample size used in the statistical analysis. In addition, we introduce the ratio of χ 2 min to the degrees of freedom (d.o.f), χ 2 min /d.o. f , to judge the quality of observational data set.
In Table 5, we list the values of AIC, BIC and χ 2 min /d.o. f for different models from the joint analysis Quasar+BAO+CMB and SN Ia+BAO+CMB. It is obvious that both of AIC and BIC criteria support Λ CDM to be the best cosmological model consistent with the available observations, since the IC value it yields is the smallest. Concerning the ranking of the three f (T ) models, AIC and BIC criteria tend to provide the same conclusions as follows. The f 2 CDM model performs the best in explaining the current data, which can be clearly seen from the similarity between f 2 CDM and Λ CDM shown in Fig. 5. Then next after f 2 CDM is the f 1 CDM model, which can also reduce to the Λ CDM model and its best-fit parameters indeed do so. The worst model according to the AIC and BIC criteria is f 3 CDM, which is unable to provide a good fit to the data and can not nest Λ CDM.

Conclusions and discussions
As an interesting approach to modify gravity, f (T ) theory based on the concept of teleparallel gravity, was proposed to explain the accelerated expansion of the Universe without the need of dark energy. In this paper, we have used the recently-released sample of VLBI observations of the compact structure in 120 intermediate-luminosity quasars (0.46 < z < 2.80) to get the constraints on the viable and most popular f (T ) gravity models. The statistical linear sizes of these quasars observed at 2.29 GHz show negligible dependence on redshifts and intrinsic luminosity, and thus represent a fixed comoving-length of the standard ruler. Therefore, the other motivation of this work was to investigate the constraining ability of quasar data in the context of f (T ) models. In particular, we have considered three f (T ) models with two parameters, out of which two could nest the concordance Λ CDM model and we quantifed their deviation from Λ CDM cosmology through a single parameter b.
For the third f (T ) cosmology which can not be directly reduced to Λ CDM, we discussed the possibility for the effective equation of state to cross the phantom divide line.
In our investigation we have used (i) the very recently released "angular size -redshift" data sets of 120 intermediateluminosity quasars in the redshift range 0.46 < z < 2.76, (ii) the cosmic microwave background and baryon acoustic oscillation data points. Meanwhile, in order to compare our fits obtained with 120 quasars (standard rulers), to the similar constraints obtained with the Union 2.1 compilation consisting of 580 SN Ia data points (standard candles) we also carried out respective analysis based on SNIa data. Here we summarize our main conclusions in more detail: -For all of the three the f (T ) models, all of the fitting results show that the quasar data (N=120) could provide more stringent constraints than the Union2.    range covered by the quasar data (0.46 < z < 2.8) compared with SN Ia (0.015 ≤ z ≤ 1.41). The constraining power of the former becomes obvious when the large size difference between the samples is taken into consideration. Moreover, one can clearly see that principal axes of confidence regions obtained with SN and quasars are inclined at higher angles, which sustains the hope that careful choice of the quasar sample would eventually provide a complementary probe breaking the degeneracy in the f (T ) model parameters. Our method based on the observations of intermediate-luminosity quasars may also contribute to testing the consistency between luminosity and angular diameter distances. -The present value of the matter density parameter Ω m implied by quasars is much lager than that derived from other observations, which has been noted by our previous analysis and the first-year Planck results, in the framework of Λ CDM cosmology. Such result indicates that quasar data at high redshifts may provide us a different understanding of the components in the Universe. -For f 1 CDM and f 2 CDM models, deviation from Λ CDM cosmology is also allowed in the obtained confidence level, although the best-fit value is very close to its Λ CDM one. It is interesting in the present work to note that Λ CDM is not included at 1σ confidence level for the power-law model f 1 CDM model, this slight deviation from Λ CDM is also consistent with a similar conclusion obtained in the previous observational studies on f (T ) gravity. In the framework of f 3 CDM, the value of n constrained by all of the current observations satisfies the limit of n > 3/2, which is necessary to achieve the cosmic acceleration. Moreover, we find that the value of n obtained with quasars suggests that the effective equa-tion of state can cross the phantom divide line at lower redshifts . -The information criteria (AIC and BIC) demonstrate that, compared with other three f (T ) scenarios considered in this paper, the cosmological constant model is still the best cosmological model consistent with the available observations. Concerning the ranking of the f (T ) cosmologies, the f 2 CDM model performs the best in explaining the current data, while the f 3 CDM model gets the smallest support and can not nest the concordance Λ CDM model. -In summary, using for the recently released quasar data acting as a new source of standard rulers, we were able to set more stringent limits on the viable and most used f (T ) gravity models. Our results highlight the importance of quasar measurements to provide additional information of various candidates for modified gravity, especially the possible deviation from Λ CDM cosmology. More importantly, given the usefulness of this angular size data in pinning down parameter values, we also anticipate that near-future quasar observations will provide significantly more restrictive constraints on other torsional modified gravity theories [11,55,56].