Constraint on the equation of state parameter ($\omega$) in non-minimally coupled $f(Q)$ gravity

We study observational constraints on the modified symmetric teleparallel gravity, the non-metricity $f(Q)$ gravity, which reproduces background expansion of the universe. For this purpose, we use Hubble measurements, Baryonic Acoustic Oscillations (BAO), 1048 Pantheon supernovae type Ia data sample which integrate SuperNova Legacy Survey (SNLS), Sloan Digital Sky Survey (SDSS), Hubble Space Telescope (HST) survey, Panoramic Survey Telescope and Rapid Response System (Pan-STARRS1). We confront our cosmological model against observational samples to set constraints on the parameters using Markov Chain Monte Carlo (MCMC) methods. We find the equation of state parameter $\omega=-0.853^{+0.015}_{-0.020}$ and $\omega= -0.796^{+0.049}_{-0.074}$ for Hubble and Pantheon samples, respectively. As a result, the $f(Q)$ model shows the quintessence behavior and deviates from $\Lambda$CDM.


I. INTRODUCTION
The current cosmic expansion of the universe motivates the scientific community to understand its fundamental properties. The present status of the universe can comprehend by testing the Λ cold dark matter (ΛCDM) model, cosmological models, and any deviation from it. Moreover, General Relativity (GR) is a wellestablished gravitational theory, but its extensions and modifications are seeking more interest due to their successful properties in describing the universe's accelerated expansion [1,2]. Hence, we aim to develop gravitational theories which have GR as a specific limit. But, in general, that formulation must include extra degree(s) of freedom to fulfill the above prerequisites.
This report will highlight the cosmological model based on the recently proposed extension of symmetric teleparallel gravity, so-called f (Q) gravity, where the non-metricity Q describes the gravitational interaction [3]. Investigations on f (Q) gravity have developed rapidly and lead to interesting applications [4][5][6][7][8][9][10][11][12][13][14][15][16][17][18][19]. If we look at the universe's expansion history, one can see that some cosmological parameters play an important role in designating the cosmological model's cosmic evolution. And, it is well known that the equation of state parameter (ω) predicts various fluid descriptions of the universe. Therefore, it is interesting to constraint this parameter using observational data. For this purpose, we use Hubble measurements, Baryonic Acoustic Oscillations (BAO), Pantheon supernovae type Ia integrates SuperNova Legacy Survey (SNLS), Sloan Digital Sky Survey (SDSS), Hubble Space Telescope (HST) survey, Panoramic Survey Telescope and Rapid Response System (Pan-STARRS1). The MCMC methods use to do the numerically analysis.
In this report, the ideas presented in the following sections. In section II, we discuss the basic setup for the f (Q) gravity. In section III, we discuss the cosmological application in FRW space-time. In section IV, we discuss the various type of observational datasets, constraint the parameters using the MCMC method, and deliberate our results. Finally, gathering all the information, we conclude in section V.

II. ACTION AND FIELD EQUATIONS
Here, we consider the action for matter coupling in f (Q) gravity is given by [3] where g is the determinant of metric, f 1 (Q) & f 2 (Q) are the arbitrary functions of the non-metricity Q, and L M is the Lagrangian for the matter fields. The nonmetricity tensor and its traces are such that Moreover, the superpotential as a function of nonmetricity tensor is given by where the trace of nonmetricity tensor [3] reads To simplify the formulation, let us introduce the following notations where primes (') represent the derivatives of functions The energy-momentum tensor for the fluid description of the spacetime can be written by its definition By varying action (1) with respect to metric tensor, we can write the gravitational field equation, which is given by Now, one can use (9) to explore the cosmological applications in f (Q) modified gravity.

III. THE f (Q) COSMOLOGY
To explore several cosmological applications, we presume the homogeneous, isotropic and spatially flat line element given by where N(t) is the lapse function and for the usual time reparametrization freedom, we can take N = 1 at any time. δ ij is the Kronecker delta and i, j run over spatial components. The expansion rate and dilation rate can be written as respectively. For this line element the non-metricity read as Q = 6(H/N) 2 . We shall work on the perfect fluid matter distribution, for which the energy momentum tensor (8) become diagonal. The gravitational equations (9) in this case generalized to two Friedman equations: (13) respectively. Here, ρ is the energy density and p is the pressure of the fluid content of the spacetime. It is easy to verify that, for f 1 = −Q and f 2 = 1 = −F, the above Friedman equations reduce to standard one [7]. From the above field equations, the continuity equation for matter field can be derived aṡ From (14), one can recover the standard continuity equation by imposing L M = −ρ aṡ This is compatible with the isotropic and homogeneous background of the universe (see more details about the continuity equation in f (Q) gravity [7]).

A. Cosmological Model
In this subsection, we shall proceed with N = 1. Since, we are working on the Friedman-Robertson-Walker (FRM) framework, the non-metricity Q and dilation rate T reduce to Moreover, the modified Friedman equations (12) and (13) can be rewritten as Now, we have two equations with five unknown such as H, ρ, p, f 1 , f 2 . To proceed further, we consider the relation p = ωρ and from equation (15), we find the following relation where ρ 0 is the proportionality constant and ω is the equation of state parameter. The relation between scale factor a(t) and redshift z read Without loss of generality, we adopt ρ 0 = a 0 = 1. In modified theories of gravity the cosmological scenarios can be discussed through the properties of cosmological models. For this purpose, we introduce two Lagrangian functions f 1 (Q) and f 2 (Q) as where α and n = 1 are the arbitrary constants. Using (19), (20), (21) in (17), we find the following expression for H(z) Now, our aim is to put constraint on the parameters α, n, ω using the astronomical observation data.

IV. DATA, METHODOLOGY AND RESULTS
This section deals with the various observational datasets to constraint the parameters α, n, ω. For this, we have adopted some statistical analysis to perform the numerical analysis. Specifically, we employ a Markov Chain Monte Carlo (MCMC) method to obtain the posterior distributions of the parameters, using the standard Bayesian technique. This stimulation is done by using the Hubble measurements ( i.e., Hubble data) and SNe Ia data. The best fits of the parameetrs are maximized by using the probability function where χ 2 is the pseudo chi-squared function [20]. The χ 2 functions for various dataset are discussed below.

A. Hubble Dataset
Recently, a list of 57 data points of Hubble parameter in the redshift range 0.07 ≤ z ≤ 2.41 were compiled by Sharov and Vasiliev [21]. This H(z) dataset was measured from the line-of-sight BAO data [22][23][24][25][26] and the differential ages ∆t of galaxies [27][28][29][30]. The complete list of datasets is presented in [21]. To estimate the model parameters, we use the Chi-square function which is given by where H obs (z i ) represents the observed Hubble parameter values, H th (p s , z i ) represents the Hubble parameter with the model parameters, and σ 2 H(z i ) is the standard deviation. In Fig. 2, the profile of our model against Hubble data shown. The marginalized constraining results are displayed in Fig. 4.

B. Pantheon Dataset
Supernovae type Ia is a powerful distance indicator to explore the background evolution of the universe. Therefore, to constraint the above parameters, we take the recent Pantheon supernovae type Ia sample, which is collectively 1048 SNe Ia data points from various SN Ia samples in the redshift-range z ∈ [0.01, 2.3] such as SDSS, SNLS, Pan-STARRS1, low-redshift survey, and HST surveys [31]. The χ 2 SN function is given to do the statistical analysis for this sample [31] by where p j represents the assumed model's free parameters and C SN is the covariance metric [31], and µ represents the distance moduli is given by;

C. Results
In Table I, we show constraints at 68% C.L. of the cosmological parameter ω and model parameters α, n for the cosmological model. Figures 1, 3 present the contour plots of the parameters at 68% and 95% C.L. for Hubble and Pantheon samples, respectively. For reference, we compare the f (Q) model with the ΛCDM model with constraints values of parameters in Figures  2 and 1. It is observed that the f (Q) model is perfectly fitting with the observational data and deviates slightly from the ΛCDM. Moreover, it is well-known that the present scenario of the universe, i.e., accelerated expan-sion, can be discussed with the presence of additional constant Λ in Einstein's field equations or by modifying the fundamental formulation of gravity for the evolution of the universe. Besides this, the equation of state parameter (ω) plays a vital role for a cosmological model to predict its' different phases of evolution. The present scenario of the universe can predict by either quintessence behavior of ω (i.e., −1 < ω < −1/3) or phantom behavior of ω (i.e., ω < −1). For our model, we found ω = −0.853 +0.015 −0.020 for Hubble dataset and ω = −0.796 +0.049 −0.074 for Pantheon dataset at 1σ confidence level. One can clearly see that the f (Q) model shows the quintessence behavior and it is near to ΛCDM.

V. CONCLUSION
The rising concern in the current scenarios of the universe motivates us to go beyond the standard formulation of gravity. In this context, we have worked on the modified f (Q) gravity to obtain observational constraints for the background candidates of the accelerated expansion of the universe. For this, we have used a wide variety of observational samples such as Hubble data and Pantheon data (which includes SDSS, SNLS, Pan-STARRS1, low-redshift survey, and HST surveys). The well-known parametrization technique is adopted to obtain the expression for H(z). The best fit ranges of the parameters are obtained by applying the Bayesian method in MCMC simulation. The constraint values of the equation of state parameter (ω) suggest that the f (Q) model shows quintessence behavior. In addition, we have depicted the profiles of Hubble parameter with the constraint values of parameters for both the dataset in Figures 2 and 4, which helps us to compare our model with the ΛCDM.
Moreover, one can compare our outcomes with the existing results to discuss the current scenario of the universe. Particularly, the cosmographic tool and observational constraints are widely used to discuss such scenarios in modified theories of gravity. For instance, one can see some of the interesting studies used these ideas, such as the cosmographic idea is used to test the observational viability of various modified theories to discuss the universe's evolution. For example, the deceleration parameter q is used to check the universe's status [32], and with the present values of cosmographic parameter are used checked the viability of class of f (T) gravity [34], class of f (Q) gravity [6]. Without a prior assumption on the equation of state parameter ω and with the cosmographic parameters the cosmological evolution studied through the profile of ω in f (R) and f (T) gravity [33]. Also, some of the observational studies in f (Q) gravity have been done in last few years [4,15,16,18]. In most of these studies, authors have presumed differ- ent values of the equation of state parameter ω such as ω = 1/3, ω = 0 to include the radiation and matter content in the fluid description of the spacetime, respectively. Whereas in our study, such assumptions on ω are avoided, and the value of ω constraint against the observational data. And, this is the advantage of our work over such types of investigations. The main goal to describe the accelerated expansion of the universe is successfully archived.
In conclusion, our findings could motivate further research into the f (Q) gravity because it is one of the alternatives to the coherence model that, aside from being preferred by the data. The significance of this model is that it does not face the cosmological constant problem because it does not comprise any additional constant inside the Lagrangian f (Q) form. In a further study, it would be interesting to explore these types of models using weak lensing data, full CMB and LSS spectra, and other datasets. Some of these tests will be addressed in the near future, and we hope to report on them.