×
近期发现有不法分子冒充我刊与作者联系,借此进行欺诈等不法行为,请广大作者加以鉴别,如遇诈骗行为,请第一时间与我刊编辑部联系确认(《中国物理C》(英文)编辑部电话:010-88235947,010-88236950),并作报警处理。
本刊再次郑重声明:
(1)本刊官方网址为cpc.ihep.ac.cn和https://iopscience.iop.org/journal/1674-1137
(2)本刊采编系统作者中心是投稿的唯一路径,该系统为ScholarOne远程稿件采编系统,仅在本刊投稿网网址(https://mc03.manuscriptcentral.com/cpc)设有登录入口。本刊不接受其他方式的投稿,如打印稿投稿、E-mail信箱投稿等,若以此种方式接收投稿均为假冒。
(3)所有投稿均需经过严格的同行评议、编辑加工后方可发表,本刊不存在所谓的“编辑部内部征稿”。如果有人以“编辑部内部人员”名义帮助作者发稿,并收取发表费用,均为假冒。
                  
《中国物理C》(英文)编辑部
2024年10月30日

Observational Hubble parameter data constraints on the interactive model of ${\boldsymbol f(\boldsymbol T) }$ gravity with particle creation

  • In this study, we consider an open system from the thermodynamic perspective for an adiabatic FRW universe model in which particle creation occurs within the system. In this case, the modified continuity equation is obtained, and then, we make it correspond to the continuity equation of $ f(T) $ gravity. Therefore, we take $ f(T) $ gravity with a viscous fluid in the flat-FRW metric, where T is the torsion scalar. We assume the contents of the universe to be dark matter and dark energy and consider an interaction term between them. An interesting point of this study is that we make the modified continuity equation resulting from particle creation equivalent to the matter continuity equation resulting from $ f(T) $ gravity. The result of this evaluation establishes a relationship between the number of particles and scale factor. In what follows, we write the corresponding cosmological parameters in terms of the number of particles and also reconstruct the number of particles in terms of the redshift parameter. We then parameterize the Hubble parameter derived from power-law cosmology with 51 data points from the Hubble observational parameter data. Next, we plot the corresponding cosmological parameters for dark energy in terms of the redshift to investigate the accelerated expansion of the universe. In addition, by using the sound speed parameter, we discuss the stability and instability analyses of the present model in different eras of the universe. Finally, we plot the density parameter values for dark energy and dark matter in terms of the redshift parameter.
      PCAS:
  • 加载中
  • [1] A. G. Riess, A. V. Filippenko, P. Challis et al., AJ 116(3), 1009 (1998) doi: 10.1086/300499
    [2] S. Perlmutter, G. Aldering, G. Goldhabe et al., AJ 517(2), 565 (1999) doi: 10.1086/307221
    [3] C. L. Bennett, M. Halpern, G. Hinshaw et al., ApJS 148(1), 1 (2003) doi: 10.1086/377253
    [4] M. Tegmark, M. A. Strauss, M. R. Blanton et al., Phys. Rev. D 69(10), 103501 (2004) doi: 10.1103/PhysRevD.69.103501
    [5] S. Weinberg, Rev. Mod. Phys. 61(1), 1 (1989) doi: 10.1103/RevModPhys.61.1
    [6] R. R. Caldwell, Phys. Lett. B 545(1), 23 (2002)
    [7] A. R. Amani, Int. J. Theor. Phys. 50(10), 3078 (2011) doi: 10.1007/s10773-011-0808-9
    [8] J. Sadeghi and A. R. Amani, Int. J. Theor. Phys. 48(1), 14 (2009) doi: 10.1007/s10773-008-9776-0
    [9] M. R. Setare, J. Sadeghi, and A. R. Amani, Phys. Lett. B 673(4), 241 (2009)
    [10] M. R. Setare, J. Sadeghi, and A. R. Amani, International Journal of Modern Physics D 18(08), 1291-1301 (2009) doi: 10.1142/S0218271809015096
    [11] R. A. Battye and F. Pace, Phys. Rev. D 94(6), 063513 (2016) doi: 10.1103/PhysRevD.94.063513
    [12] M. Li, T. Qiu, Y. Cai et al., Journal of Cosmology and Astroparticle Physics 2012(04), 003 (2012) doi: 10.1088/1475-7516/2012/04/003
    [13] B. Pourhassan and J. Naji, International Journal of Modern Physics D, 1750012 (2014)
    [14] A. R. Amani and S. L. Dehneshin, Canadian Journal of Physics 93, 1453-1459 (2015) doi: 10.1139/cjp-2015-0024
    [15] V. Faraoni, Physics of the Dark Universe 11, 11-15 (2016) doi: 10.1016/j.dark.2015.11.001
    [16] J. Sadeghi, B. Pourhassan, A. S. Kubeka et al., International Journal of Modern Physics D 25(07), 1650077 (2016) doi: 10.1142/S0218271816500772
    [17] H. Wei, Communications in Theoretical Physics 52(4), 743 (2009) doi: 10.1088/0253-6102/52/4/35
    [18] A. R. Amani and A. Samiee-Nouri, Communications in Theoretical Physics 64(4), 485 (2015) doi: 10.1088/0253-6102/64/4/485
    [19] S. Nojiri and S. D. Odintsov, Phys. Lett. B 657(4), 238-245 (2007)
    [20] M. Li, Phys. Lett. B 603(1), 1-5 (2004)
    [21] S. Del Campo, J. C. Fabris, R. Herrera et al., Phys. Rev. D 83(12), 123006 (2011) doi: 10.1103/PhysRevD.83.123006
    [22] Y. Hu, M. Li, N. Li et al., Journal of Cosmology and Astroparticle Physics 2015(08), 012 (2015) doi: 10.1088/1475-7516/2015/08/012
    [23] A. R. Amani, C.Escamilla-Rivera, and H. R. Faghani, Phys. Rev. D 88, 124008 (2013) doi: 10.1103/PhysRevD.88.124008
    [24] A. R. Amani and B. Pourhassan, International Journal of Geometric Methods in Modern Physics 11(08), 1450065 (2014) doi: 10.1142/S0219887814500650
    [25] J. Naji, B. Pourhassan, and A. R. Amani, International Journal of Modern Physics D 23(02), 1450020 (2014) doi: 10.1142/S0218271814500205
    [26] J. Morais, M. Bouhmadi-Lopez, K. Sravan Kumar et al., Physics of the Dark Universe 15, 7-30 (2017) doi: 10.1016/j.dark.2016.11.002
    [27] Y. Zhang, Physics of the Dark Universe 15, 82 (2017) doi: 10.1016/j.dark.2016.12.003
    [28] J. Sadeghi, B. Pourhassan, and Z. Abbaspour Moghaddam, Int. J. Theor. Phys. 53(1), 125-135 (2014) doi: 10.1007/s10773-013-1790-1
    [29] J. Sadeghi, F. Milani, and A. R. Amani, Mod. Phys. Lett. A 24(29), 2363-2376 (2009) doi: 10.1142/S0217732309030783
    [30] J. Sadeghi, M. R. Setare, A. R. Amani et al., Phys. Lett. B 685(4), 229-234 (2010)
    [31] A. R. Amani, International Journal of Modern Physics D 25(06), 1650071 (2016) doi: 10.1142/S0218271816500711
    [32] T. Singh, R. Chaubey, and A. Singh, Canadian Journal of Physics 94(7), 623-627 (2016) doi: 10.1139/cjp-2016-0081
    [33] V. Sahni, and Y. Shtanov, Journal of Cosmology and Astroparticle Physics 2003(11), 014 (2003) doi: 10.1088/1475-7516/2003/11/014
    [34] M. R. Setare, J. Sadeghi, and A. R. Amani, Phys. Lett. B 660(4), 299 (2008) doi: 10.1016/j.physletb.2008.01.024
    [35] G. P. de Brito, J. M. Hoff da Silva, P. Michel LT da Silva et al., International Journal of Modern Physics D 24(11), 1550089 (2015) doi: 10.1142/S0218271815500893
    [36] M. R. Setare, J. Sadeghi, and A. R. Amani., Phys. Lett. B 666(4), 288-298 (2008) doi: 10.1016/j.physletb.2008.07.087
    [37] A. R. Amani and H. Farahani, Int. J. Theor. Phys. 51(5), 1498-1502 (2012) doi: 10.1007/s10773-011-1025-2
    [38] A. R. Amani and H. Farahani, Int. J. Theor. Phys. 51(9), 2943-2949 (2012) doi: 10.1007/s10773-012-1174-y
    [39] A. R. Amani and B. Pourhassan, Int. J. Theor. Phys. 51(1), 49-54 (2012) doi: 10.1007/s10773-011-0876-x
    [40] S. R. Bhoyar, V. R. Chirde, and S. H. Shekh, Astrophys. 60(2), 259-272 (2017) doi: 10.1007/s10511-017-9480-y
    [41] V. R. Chirde and S. H. Shekh, J. Astrophys. Astron. 39(5), 56 (2018) doi: 10.1007/s12036-018-9555-0
    [42] J. K. Singh, K. Bamba, R. Nagpal et al., Phys. Rev. D 97(12), 123536 (2018) doi: 10.1103/PhysRevD.97.123536
    [43] R. Nagpal, S. K. J. Pacif, J. K. Singh et al., Eur. Phys. J. C 78(11), 946 (2018) doi: 10.1140/epjc/s10052-018-6403-y
    [44] A. R. Amani, J. Sadeghi, H. Farajollahi et al., Canadian Journal of Physics 90(1), 61 (2011)
    [45] H. Wei and R. G. Cai, Phys. Lett. B 663, 1 (2008) doi: 10.1016/j.physletb.2008.03.048
    [46] H. Wei and R.G. Cai, Phys. Lett. B 665, 1 (2007)
    [47] A. Einstein, Sitzungsberichte der Preussischen Akademie der Wissenschaften, Physikalisch-mathematische Klasse 17, 224227 (1928)
    [48] R. Weitzenbock, Noordhoff, Groningen (1923)
    [49] E. V. Linder, Phys. Rev. D 81(12), 127301 (2010) doi: 10.1103/PhysRevD.81.127301
    [50] R. Myrzakulov, Eur. Phys. J. C 71(9), 1-8 (2011)
    [51] R. Myrzakulov, Eur. Phys. J. C 72(11), 2203 (2012) doi: 10.1140/epjc/s10052-012-2203-y
    [52] B. Li, T. P. Sotiriou, and J. D. Barrow, Phys. Rev. D 83(6), 064035 (2011) doi: 10.1103/PhysRevD.83.064035
    [53] T. Harko, F. S. N. Lobo, S. Nojiri et al., Phys. Rev. D 84(2), 024020 (2011) doi: 10.1103/PhysRevD.84.024020
    [54] S. Nojiri and S. D. Odintsov, Phys. Lett. B 631, 1-6 (2005) doi: 10.1016/j.physletb.2005.10.010
    [55] Y. B. ZelDovich and A. A. Starobinskii, Soviet Journal of Experimental and Theoretical Physics 34, 1159 (1972)
    [56] I. Prigogine, Int. J. Theor. Phys. 28(9), 927-933 (1989) doi: 10.1007/BF00670337
    [57] L. Parker, Journal of Physics A: Mathematical and Theoretical 45(37), 374023 (2012) doi: 10.1088/1751-8113/45/37/374023
    [58] V. H. Cardenas, Eur. Phys. J. C 72(9), 1-6 (2012)
    [59] J. C. Fabris, J. A. de Freitas Pacheco, and O. F. Piattella, Journal of Cosmology and Astroparticle Physics 2014(06), 038 (2014) doi: 10.1088/1475-7516/2014/06/038
    [60] S. K. Biswas, W. Khyllep, J. Dutta et al., Phys. Rev. D 95(10), 103009 (2017) doi: 10.1103/PhysRevD.95.103009
    [61] S. Saha, and A. Mondal, Eur. Phys. J. C 77(3), 1-6 (2017)
    [62] F. Q. Tu, Y. X. Chen, and Q. H. Huang, Chin. Phys. C 44(6), 065103 (2020) doi: 10.1088/1674-1137/44/6/065103
    [63] S. Mandal and P. K. Sahoo, Mod. Phys. Lett. A 35(40), 2050328 (2020) doi: 10.1142/S0217732320503289
    [64] L. Parker, Phys. Rev. Lett. 21(8), 562 (1968) doi: 10.1103/PhysRevLett.21.562
    [65] S. W. Hawking, Commun. Math. Phys. 43, 199 (1975) doi: 10.1007/BF02345020
    [66] M. Sharif and S. Rani, Mod. Phys. Lett. A 28(27), 1350118 (2013) doi: 10.1142/S0217732313501186
    [67] D. Pavón, and B. Wang, Gen. Relativ. Gravit. 41(1), 1-5 (2009) doi: 10.1007/s10714-008-0656-y
    [68] W. Zimdahl, D. Pavón, and L. P. Chimento, Phys. Lett. B 521(3-4), 133-138 (2001) doi: 10.1016/S0370-2693(01)01174-1
    [69] J. Dutta, W. Khyllep, and N. Tamanini, Journal of Cosmology and Astroparticle Physics 2018(01), 038 (2018) doi: 10.1088/1475-7516/2018/01/038
    [70] N. Aghanim et al., Astronomy and Astrophysics 641, A6 (2020) doi: 10.1051/0004-6361/201833910
    [71] K. Bamba, C. Q. Geng, C. C. Lee et al., Journal of Cosmology and Astroparticle Physics 2011(01), 021 (2011) doi: 10.1088/1475-7516/2011/01/021
    [72] Y. F. Cai, S. Capozziello, M. De Laurentis et al., Reports on Progress in Physics 79(10), 106901 (2016) doi: 10.1088/0034-4885/79/10/106901
    [73] D. Shafer, Phys. Rev. D 91(10), 103516 (2015) doi: 10.1103/PhysRevD.91.103516
    [74] I. Tutusaus, B. Lamine, A. Blanchard et al., Phys. Rev. D 94(10), 103511 (2016) doi: 10.1103/PhysRevD.94.103511
    [75] C. Blake, S. Brough, M. Colless et al., Monthly Notices of the Royal Astronomical Society 425(1), 405-414 (2012) doi: 10.1111/j.1365-2966.2012.21473.x
    [76] A. Font-Ribera, D. Kirkby et al., 2014, JCAP, 1405, 027. Journal of Cosmology and Astroparticle Physics 2014, no. 05 (2014): 027
    [77] T. Delubac, J. E. Bautista, N. G. Busca et al., Astronomy and Astrophysics 574, A59 (2015) doi: 10.1051/0004-6361/201423969
    [78] U. Alam, S. Bag, and V. Sahni, Phys. Rev. D 95(2), 023524 (2017) doi: 10.1103/PhysRevD.95.023524
    [79] M. Moresco, L. Pozzetti, A. Cimatti et al., Journal of Cosmology and Astroparticle Physics 2016(05), 014 (2016) doi: 10.1088/1475-7516/2016/05/014
    [80] O. Farooq, F. R. Madiyar, S. Crandall et al., The Astrophysical Journal 835(1), 26 (2017) doi: 10.3847/1538-4357/835/1/26
    [81] S. K. J. Pacif, R. Myrzakulov, and S. Myrzakul, International Journal of Geometric Methods in Modern Physics 14(07), 1750111 (2017) doi: 10.1142/S0219887817501110
    [82] J. Magana, M. H. Amante, M. A. Garcia-Aspeitia et al., Monthly Notices of the Royal Astronomical Society 476(1), 1036-1049 (2018) doi: 10.1093/mnras/sty260
    [83] E. Mahichi, A. Amani, and M. A. Ramzanpour, arXiv: 2107.11827 (2021)
    [84] D. M. Scolnic, D. O. Jones, A. Rest et al., AJ 859(2), 101 (2018) doi: 10.3847/1538-4357/aab9bb
  • 加载中

Figures(5)

Get Citation
S. Ganjizadeh, Alireza Amani and M. A. Ramzanpour. Observational Hubble parameter data constraints on the interactive model of ${\boldsymbol f(\boldsymbol T) }$ gravity with particle creation[J]. Chinese Physics C. doi: 10.1088/1674-1137/ac8c22
S. Ganjizadeh, Alireza Amani and M. A. Ramzanpour. Observational Hubble parameter data constraints on the interactive model of ${\boldsymbol f(\boldsymbol T) }$ gravity with particle creation[J]. Chinese Physics C.  doi: 10.1088/1674-1137/ac8c22 shu
Milestone
Received: 2022-06-22
Article Metric

Article Views(1371)
PDF Downloads(39)
Cited by(0)
Policy on re-use
To reuse of subscription content published by CPC, the users need to request permission from CPC, unless the content was published under an Open Access license which automatically permits that type of reuse.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

Observational Hubble parameter data constraints on the interactive model of ${\boldsymbol f(\boldsymbol T) }$ gravity with particle creation

  • Department of Physics, Ayatollah Amoli Branch, Islamic Azad University, Amol, Iran

Abstract: In this study, we consider an open system from the thermodynamic perspective for an adiabatic FRW universe model in which particle creation occurs within the system. In this case, the modified continuity equation is obtained, and then, we make it correspond to the continuity equation of $ f(T) $ gravity. Therefore, we take $ f(T) $ gravity with a viscous fluid in the flat-FRW metric, where T is the torsion scalar. We assume the contents of the universe to be dark matter and dark energy and consider an interaction term between them. An interesting point of this study is that we make the modified continuity equation resulting from particle creation equivalent to the matter continuity equation resulting from $ f(T) $ gravity. The result of this evaluation establishes a relationship between the number of particles and scale factor. In what follows, we write the corresponding cosmological parameters in terms of the number of particles and also reconstruct the number of particles in terms of the redshift parameter. We then parameterize the Hubble parameter derived from power-law cosmology with 51 data points from the Hubble observational parameter data. Next, we plot the corresponding cosmological parameters for dark energy in terms of the redshift to investigate the accelerated expansion of the universe. In addition, by using the sound speed parameter, we discuss the stability and instability analyses of the present model in different eras of the universe. Finally, we plot the density parameter values for dark energy and dark matter in terms of the redshift parameter.

    HTML

    I.   INTRODUCTION
    • In the late twentieth century, by measuring cosmic distances and observing cosmic objects beyond the Milky Way, we realized that the universe is rapidly expanding, and this event occurs when there is a force greater than gravity but in the opposite direction in the universe. However, the origin of this mysterious force is introduced by an unknown energy called dark energy, which is supposed to surround the entire space and increase the expansion rate of the universe. In fact, this energy was first discovered in type-Ia supernovae and then reconfirmed by two approaches to the cosmic microwave background radiation and the large-scale structure of the universe [14]. The problem of dark energy originates from a fundamental theoretical framework such as string theory and quantum gravity. However, despite the existence of evidence for the fact that the universe is accelerating, these observations do not fit the standard model of cosmology within the framework of general relativity. To resolve this contradiction, we introduce a factor called negative pressure for the dark energy source, which can describe the observed accelerated expansion of the universe. Extensive studies of dark energy have been explored with models such as the cosmological constant, scalar fields, modified gravity, agegraphic, holography, and teleparallel gravity [546].

      In addition, another scenario called teleparallel gravity, first proposed by Einstein, was introduced for the unification of electromagnetism and gravity. In this theory, space-time is characterized by a linear connection without curvature by a metric tensor field called the dynamical tetrad field. Therefore, the geometry of the teleparallel theory uses the linear Weitzenböck connection to represent a torsion connection without curvature. However, the geometry of general relativity is based on the Levi-Civita connection to represent a curvature connection without torsion [4754].

      As mentioned above, current observations with different scenarios suggest that the universe is undergoing accelerated expansion, which leads to the crossing of the phantom divided-line (equation of state (EoS) is less than –1). Furthermore, the particle creation mechanism is an appropriate alternative to describe the current acceleration of the universe, so that it can successfully mimic standard cosmology. In that case, from the point of view of the first law of thermodynamics, the creation of a dynamical number of particles, N, leads to a modified continuity equation for a system called the adiabatic universe [5563]. Note that the corresponding continuity equation is also obtained via Einstein's equation in standard cosmology as $ \dot \rho + 3 H (\rho + p) = 0 $, in which ρ, p, and H are the energy density, pressure, and Hubble parameter, respectively. An important point that should be noted in the mechanism of particle creation is that its continuity equation, in addition to standard terms, also has the term of the rate of the number of particles as $({\dot N}/{N}) (\rho+p) $. A number of papers attribute the origin of particle number rate to the bulk viscosity. In this event, an open thermodynamic system in which the number of fluid particles is not conserved leads to a modified particle conservation equation, known as the equilibrium equation for particle flow. Hence, we consider that the corresponding particle flow is related to the energy flow between the dark components of the universe. Then, the particle flux resulting from particle creation depends on the energy transfer between dark energy and dark matter. In that case, we consider this issue as energy transfer, the so-called interaction between the dark components of the universe.

      It should be noted that the particle creation in a black hole and the expanding universe was also studied [64, 65]. For this purpose, in this work, we intend to implement the mechanism of the particle creation to describe the accelerated expansion of the universe. This motivated us to make a correspondence between the particle creation mechanism and the $ f(T) $ gravity model. Therefore, using this correspondence and a power-law model for the scale factor, we expect to examine the accelerated expansion of the universe with observational data.

      This paper is organized as follows:

      In Sec. II, we implement the fundamental formation of viscous $ f(T) $ gravity in the flat-FRW metric. In Sec. III, we consider that the universe is dominated by a bulk viscosity fluid, and also obtain the Friedmann equations and EoS in an interacting model. In Sec. IV, we first explain the mechanism of particle creation, and then using the interaction term, we correspond the particle creation with the viscous $ f(T) $ gravity model. In what follows, in Sec. V, we investigate the present model by observational constraints with Hubble data, and then the cosmological parameters are plotted in terms of the redshift parameter. In Sec. VI, we explore the stability conditions for the corresponding model. Finally, in Sec. VII, we present a conclusion for the present model.

    II.   FOUNDATION OF f(T) GRAVITY
    • In this section, we review teleparallel gravity. For this purpose, we can write the corresponding action for $ f(T) $ gravity in the following form:

      $ \begin{equation} S = \int e \left(\frac{f(T)}{2 \kappa^2} + \mathcal{L}_m \right) {\rm d}^4x, \end{equation} $

      (1)

      where $ \kappa^2=8 \pi G $, and $ e^i_{\,\,\mu} $ is a vierbein field from which the determinant of the vierbein field is obtained as $ e={\rm det}\left(e^i_{\,\,\mu}\right) $, and $ \mathcal{L}_m $ is the matter Lagrangian density. For a more complete explanation, we write the basis of the vierbein field $ e_i(x^\mu) $ as the tangent vectors at point $ x^\mu $ of a manifold as $ e_i. e_j = \eta_{ij} $, in which we consider the corresponding signature of Minkowski space-time as $ \eta_{ij}={\rm diag}(+1,-1,-1,-1) $. In this case, by using the relationship $ e_i = e_i\,^\mu \partial_\mu $, we can obtain metric tensor $ g_{\mu \nu} $ by the tetrad or vierbein field $ e^i_{\,\,\mu} $ as $ g_{\mu \nu}=\eta_{\mu \nu} e^i_{\,\,\mu} e^i_{\,\,\nu} $, which in turn leads to the relationships $ e_i\,^\mu e^i\,_\nu=\delta_\nu^\mu $ and $ e_i\,^\mu e^j\,_\mu=\delta_i^j $. It should be noted that the Greek and Latin letters denote components of space-time coordinates and tangent space-time coordinates, respectively.

      As we know, we use the Weitzenböck connection in teleparallel gravity instead of the Levi-Civita connection in general relativity as

      $ \begin{equation} \Gamma^\lambda_{\; \mu\nu} = e_i^{\; \lambda} \partial_\mu e^i_{\; \nu} = -e^i_{\; \mu} \partial_\nu e_i^{\; \lambda}. \end{equation} $

      (2)

      By using this connection, we can introduce the torsion tensor and the torsion scalar in the following form:

      $ \begin{eqnarray} &{T^\lambda }_{\mu \nu } = \Gamma^\lambda_{\; \mu\nu} - \Gamma^\lambda_{\; \nu \mu} = {e_a}^\lambda \,\left( {{\partial _\mu }{e^a}_\nu - {\partial _\nu }{e^a}_\mu } \right), \end{eqnarray} $

      (3)

      $ \begin{eqnarray} &T = {S_\lambda }^{\mu \nu }\,{T^\lambda }_{\mu \nu }, \end{eqnarray} $

      (4)

      where

      $ \begin{eqnarray} &S_\lambda ^{\mu \nu } = \frac{1}{2}\,\left( {K_\lambda ^{\mu \nu } + \delta _\lambda ^\mu \,T_\alpha ^{\alpha \nu } - \delta _\lambda ^\nu \,T_\alpha ^{\alpha \mu }} \right), \end{eqnarray} $

      (5)

      $ \begin{eqnarray} &K_\lambda ^{\mu \nu } = - \frac{1}{2}\,\left( {{T^{\mu \nu }}_\lambda - T_\lambda ^{\nu \mu } - T_\lambda ^{\mu \nu }} \right), \end{eqnarray} $

      (6)

      to be known as the antisymmetric tensor and the contortion tensor, respectively. We can also obtain the Ricci tensor and the Ricci scalar by the aforesaid relations in the following form:

      $ \begin{eqnarray} & R_{\mu \nu} = -\nabla^\rho S_{\nu \rho \mu} - g_{\mu \nu} \nabla^\rho T^\sigma_{\; \rho \sigma} - S^{\rho \sigma}_{\; \; \mu} K_{\sigma \rho \nu}, \end{eqnarray} $

      (7)

      $ \begin{eqnarray} &R = -T - 2 \nabla^\mu \left(T^\nu_{\; \mu \nu}\right). \end{eqnarray}$

      (8)

      By taking variation of action (1) with respect to the tetrad field, the field equations are obtained as

      $\begin{aligned}[b] & e^{-1} \partial_{\mu}\left(e\, e^\lambda_a\, S_{\lambda}^{\; \mu\nu}\right) f_{T} - e_{a}^{\rho} \, T^{\lambda}_ {\; \mu\rho} S_{\lambda}{}^{\nu\mu} f_{T} \\&\quad+ e^\lambda_a S_{\lambda}{}^{\mu\nu}\, \partial_{\mu}(T) f_{TT} +\frac{1}{4}e_{a}^{\nu} f = \frac{1}{2} \kappa^2 e_{a}^{\lambda} \mathcal{T}_{\lambda}{}^{\nu}, \end{aligned} $

      (9)

      where index T represents the derivative with respect to torsion scalar, and $ \mathcal{T}_{\lambda}{}^{\nu} $ is the energy-momentum tensor. Here, we consider the background manifold as a flat Friedmann-Robertson-Walker (FRW) metric as

      $ \begin{equation} {\rm d}s^{2} = {\rm d}t^{2} - a^{2}(t) \left({\rm d}x^2 +{\rm d}y^2 + {\rm d}z^2 \right), \end{equation} $

      (10)

      where $ a(t) $ is the scale factor. Then, one can immediately obtain the tetrad field as $ e_{\mu}^a=\mathrm{diag}(1,a,a,a) $ or $ e^{\mu}_a= \mathrm{diag}(1,a^{-1},a^{-1},a^{-1}) $, and the torsion scalar is obtained in the following form:

      $ \begin{eqnarray} T = -6 H^2, \end{eqnarray} $

      (11)

      where $ H = \dfrac{\dot{a}}{a} $ is introduced as the Hubble parameter.

    III.   VISCOUS f(T) GRAVITY
    • In this section, we intend to consider that the current universe is dominated by an anisotropic fluid with the view that the universe is in a more realistic form. This means that we would like to describe the corresponding universe by $ f(T) $ gravity in the presence of a bulk viscous fluid. Considering that the fluid viscosity indicates its resistance to flow as a dissipative phenomenon, the pressure within the universe becomes negative. Then, the presence of bulk viscosity causes the accelerated expansion of the universe [66]. In this case, we can write the energy-momentum tensor in the presence of bulk viscosity as follows:

      $ \begin{equation} \mathcal{T}_i^j=(\rho_{\rm tot} + p_{\rm tot} + p_{\rm b}) u_i u^j - \left(p_{\rm tot} + p_{\rm b}\right)\, \delta_i^j, \end{equation} $

      (12)

      where $ \rho_{\rm tot} $, $ p_{\rm tot} $, and $ p_{\rm b} $ are introduced as the total energy density, total pressure, and pressure of bulk viscosity of fluid within the universe, respectively. We use the $ 4 $-velocity $ u_\mu $ notation as $ u^i $ = (+1,0,0,0), which gives rise to $ u_i u^j $ = 1. Therefore, we can acquire the non-zero energy-momentum tensor elements as

      $ \begin{equation} \mathcal{T}_i^j ={\rm diag} \left( \rho_{\rm tot}, - p_{\rm tot} - p_{\rm b}, - p_{\rm tot} - p_{\rm b}, - p_{\rm tot} - p_{\rm b}\right), \end{equation} $

      (13)

      where by substituting the aforesaid energy-momentum tensor elements and FRW metric into Eq. (6), the Friedmann equations are obtained as

      $ \begin{eqnarray} & \kappa^2 \rho_{\rm tot} = 6 H^2 \partial_T f + \frac{1}{2} f, \end{eqnarray} $

      (14)

      $ \begin{eqnarray} & \kappa^2 \left({p}_{\rm tot} +p_{\rm b}\right) = -2 \left(\dot{H} + 3 H^2 \right) \partial_T f + 24 H^2 \dot{H} \partial_{T T} f - \frac{1}{2} f, \end{eqnarray} $

      (15)

      where the dot and index T indicate the derivative with respect to time evolution and the torsion scalar, respectively. Note that the Friedmann Eqs. (14) and (15) are written in the standard form in general relativity by applying condition $ f(T) = T $ in a perfect fluid ($ p_{\rm b} = 0 $). Therefore, we consider that the components of the universe include matter and dark energy in an anisotropic fluid ($ p_{\rm b} \neq 0 $). In this case, we write the total energy density and total pressure in terms of the components of the universe in the following form:

      $ \begin{eqnarray} & \rho_{\rm tot} = \rho + \rho_{\rm de}, \end{eqnarray} $

      (16)

      $ \begin{eqnarray} & {p}_{\rm tot} = {p} + p_{\rm de}, \end{eqnarray} $

      (17)

      where ρ and p are the energy density and pressure of matter, and $ \rho_{\rm de} $ and $ p_{\rm de} $ are the energy density and pressure of dark energy. However, from Eqs. (14), (15) and (16), (17), we obtain relations $ \rho_{\rm de} $ and $ p_{\rm de} $ in the following form:

      $ \begin{eqnarray} & \kappa^2 \rho_{\rm de} = 6 H^2 \partial_T f + \frac{1}{2} f - \kappa^2 \rho, \end{eqnarray} $

      (18)

      $ \begin{aligned}[b] \kappa^2 {p}_{\rm de} =& -2 \left(\dot{H} + 3 H^2 \right) \partial_T f + 24 H^2 \dot{H} \partial_{T T} f\\& - \frac{1}{2} f - \kappa^2 \rho \, \omega_m - 3 \kappa^2 \xi H, \end{aligned} $

      (19)

      where $ \omega_m = p/\rho $ is the matter equation of state (EoS), which is a constant, and the term $ 3 \xi H $ is substituted with the pressure of bulk viscosity, $ p_{\rm b} $, in which the coefficient ξ is a positive constant. Afterward, we can introduce the EoS of dark energy as follows:

      $ \begin{aligned}[b] \omega_{\rm de} =& \frac{p_{\rm de}}{\rho_{\rm de}} \\=& -1 - \frac{2 \dot{H} \partial_T f - 24 H^2 \dot{H} \partial_{T T} f + \kappa^2 \rho \, (1+ \omega_m) + 3 \kappa^2 \xi H}{6 H^2 \partial_T f + \dfrac{1}{2} f - \kappa^2 \rho}, \end{aligned} $

      (20)

      where the corresponding EoS of dark energy describes the dynamic of the universe in late time and depends on the modified gravity, matter, and viscous fluid.

      In what follows, we can obtain the total continuity equation for the total fluid inside the universe as

      $ \begin{equation} \dot{\rho}_{\rm tot}+3 H \left(\rho_{\rm tot} + {p}_{\rm tot} + p_{\rm b}\right)=0, \end{equation} $

      (21)

      which can be written separately in terms of the universe components in the presence of an interacting term in the form

      $ \begin{eqnarray} &\dot{\rho}+3 H \left(\rho + p \right)= \mathcal{Q}, \end{eqnarray}$

      (22)

      $ \begin{eqnarray} &\dot{\rho}_{\rm de}+3 H \left(\rho_{\rm de} + {p}_{\rm de} + p_{\rm b}\right)= -\mathcal{Q}, \end{eqnarray} $

      (23)

      where $ \mathcal{Q} $ is introduced as an interaction term. Note that the interaction term appears when the energy flow is transferred between the dark components of the universe. It is evident that the quantity $ \mathcal{Q} $ should be positive, that is, energy transfer from dark energy to dark matter occurs. Therefore, the positivity of the interaction term guarantees that the second law of thermodynamics is realized [67]. In that case, since the unit $ \mathcal{Q} $ is the inverse of the time evolution, it is natural to choose this value as the product of the Hubble parameter and energy density. Herein, we consider $ \mathcal{Q} = 3 b^2 H \rho $, in which $ b^2 $ is the coupling parameter or transmission intensity.

      In the next section, we discuss the thermodynamic study of the particle creation in an open system. Then, the corresponding open system is considered an adiabatic universe, so we see how the particle creation is related to the interaction term between the dark parts of the universe.

    IV.   CORRESPONDENCE BETWEEN MATTER CREATION AND VISCOUS f(T) GRAVITY
    • As we know, the first law of thermodynamics, or in other words, the Gibbs equation, with the dynamic number of particles N and enthalpy density $ h = \rho + p $ in an open system is written as follows:

      $ \begin{equation} {\rm d}E = {\rm d}Q - p {\rm d}V + \frac {h}{n} {\rm d}(nV), \end{equation} $

      (24)

      where $ \rho = \dfrac{E}{V} $ and $ n = \dfrac{N}{V} $ are the energy density and density of particles, respectively, where E and V are the internal energy and volume of the system, respectively.

      Now, we consider the open system as an adiabatic FRW universe model as $ {\rm d}Q = 0 $ with the creation or elimination of the number of particles in the corresponding system. In order to realize the creation of the particles, we consider the volume of the universe as a sphere with a radius of the scale factor $ a(t) $ in the form of the relationship $ V = \dfrac{4 \pi}{3}a^3 $, so that Eq. (24) yields

      $ \begin{equation} \dot \rho - (\rho + p) \frac{\dot{n}}{n} = 0. \end{equation} $

      (25)

      By taking the derivative from the relationship $ N = n V $ with respect to time evolution, we can clearly obtain the following relation:

      $ \begin{equation} \frac {\dot n}{n} = \frac {\dot N}{N} -3 H, \end{equation} $

      (26)

      where we will immediately have

      $ \begin{equation} \dot \rho + 3 H (\rho + p) - \frac {\dot N}{N} (\rho+p) = 0, \end{equation} $

      (27)

      and this equation is introduced as the modified continuity equation by particle creation. In this case, the third sentence provides negative pressure for matter inside the universe, which causes the present accelerated expansion. Then, we clearly obtain the solution of Eq. (27) as below:

      $ \begin{equation} \rho = C \left(\frac{N}{a^3}\right)^{1 + \omega_m} = C \, n^{1 + \omega_m}. \end{equation} $

      (28)

      where C is an integral constant.

      Now, we can consider that the matter component within the universe is dominated by the dynamical number of particles N; therefore, the matter creation is related to theories of gravity and cosmology. It is interesting to note that researchers consider the negative pressure for Eq. (27) by a viscous fluid, but in this study, we tackle this issue with the interacting model. The existence of interaction between dark energy and dark matter affects the formation of the structure and evolution of the universe. Therefore, the coupling between dark energy and dark matter leads to the evolution of a uniform distribution of baryons in the formation of the structure of the universe and instability in the dark parts from the early era to the dominance of dark energy in the late era [68, 69]. This means that particle creation is considered as an interaction term between the universe components, that is, particle creation is expressed by the conditions $ \mathcal{Q} > 0 $ and $ \mathcal{Q} < 0 $, for the source and sink of energy flow, respectively. In this case, by applying $ \mathcal{Q} > 0 $ and comparing Eqs. (22) and (27), we find the interaction term $ \mathcal{Q} $ with particle creation in the following form

      $ \begin{equation} \mathcal{Q} = \frac {\dot N}{N} \rho (1 + \omega_m), \end{equation} $

      (29)

      where the rate of matter production is immediately obtained as

      $ \begin{equation} \frac {\dot N}{N} = \frac{3 b^2}{1 + \omega_m} H, \end{equation} $

      (30)

      where N is given as follows:

      $ \begin{equation} N = c \, a^{{3 b^2}/({1 + \omega_m})}\; \; \; {\rm or}\; \; \; a = N^{({1 + \omega_m})/{3 b^2}}, \end{equation} $

      (31)

      where c is an integral constant, assumed to be equal to 1. By substituting Eq. (31) into Eq. (28), we obtain the energy density in terms of the scale factor as follows:

      $ \begin{equation} \rho = C \, a^{3 \left(b^2 - 1 - \omega_m \right)}, \end{equation} $

      (32)

      which shows that, by applying $ b = 0 $ and $ \omega_m = 0 $, we reach the dominated universe with dust-like matter. Note that Eq. (32) satisfies Eq. (22). However, Eq. (32) is directly derived from (22), but with the particle creation approach, we can obtain the number of particles, N, as an explicit function in the form of Eq. (31). In contrast, we write the density parameter, $ \Omega_m $, as follows:

      $ \begin{equation} \Omega_m = \frac{\rho}{\rho_{c}}, \end{equation} $

      (33)

      where $ \rho_c = 3 H^2 / \kappa^2 $ is the critical density. We describe the present density parameter as

      $ \begin{equation} \Omega_{m_0} = \frac{\kappa^2 \rho_0}{3 H_0^2}, \end{equation} $

      (34)

      where $ H_0 = 67.4 \pm 0.5 \; {\rm km\,s^{-1}\,Mpc^{-1}} $ is the current Hubble parameter, and here, we take $ \Omega_{m_0} = 0.315 \pm 0.007 $ [70]. The matter energy density, ρ, is given by Eqs. (32) and (34) as

      $ \begin{equation} \rho = \frac{3 H_0^2 \Omega_{m_0}}{\kappa^2} \left(\frac{a}{a_0}\right)^{3 \left(b^2 - 1 - \omega_m \right)}, \end{equation} $

      (35)

      where $ a_0 = 1 $ is the current scale factor.

    V.   OBSERVATIONAL CONSTRAINTS WITH HUBBLE PARAMETER DATA
    • As mentioned in the previous section, we considered an open system from the thermodynamic point of view for an adiabatic FRW universe model, where there is no heat exchange with the outside universe but the number of particles within the system is created. Therefore, the presence of particles within the universe leads to a modified continuity equation that is related to the interaction term compared to the continuity equation of $ f(T) $ gravity. With this perspective, we treat the number of particles N as a function of the scale factor in the form of Eq. (31). Therefore, in this work, the perspective of the particle creation from the point of view of $ f(T) $ gravity is considered.

      In what follows, we insert Eq. (32) into Eqs. (18), (19) and then we have

      $ \begin{eqnarray} & \kappa^2 \rho_{\rm de} = 6 H^2 \partial_T f + \frac{1}{2} f - 3 H_0^2\, \Omega_{m_0} \, a^{3\left(b^2 - 1 - \omega_m \right)}, \end{eqnarray} $

      (36)

      $ \begin{aligned}[b] \kappa^2 {p}_{\rm de} =& -2 \left(\dot{H} + 3 H^2 \right) \partial_T f + 24 H^2 \dot{H} \partial_{T T} f \\&- \frac{1}{2} f - 3 \omega_m\, H_0^2\, \Omega_{m_0} \, a^{3 \left(b^2 - 1 - \omega_m \right)} - 3 \kappa^2 \xi H. \end{aligned} $

      (37)

      In order to solve the aforesaid Friedmann equations, we need to find the form of the function $ f(T) $ in terms of the torsion scalar. For this purpose, if we compare the Friedmann Eq. (14) with the standard form of the Friedmann equation, we find function $ f(T) $ as $ f(T) = T - \alpha \sqrt{T} $. Since we are dealing with the modified Friedmann equation (14), the form of function $ f(T) $ is expected to be different from the form of function $ f(T) = T - \alpha \sqrt{T} $. For this purpose, we take the form of function $ f(T) $ as follows:

      $ \begin{equation} f(T) = T + \alpha T \left(1 - {\rm e}^{{\beta T_0}/{T}} \right) - \alpha T_0 \sqrt{\frac{T}{\beta T_0}} \, \ln\left(\frac{\beta T_0}{T}\right), \end{equation} $

      (38)

      where α and β are dimensionless parameters, and $ T_0 = -6 H_0^2 $ is the current amount of the torsion scalar which arises from Eq. (11). The above particular model is expected to demonstrate the crossing of the phantom divided-line, which can solve dark energy as a serious challenge in cosmology. It should be noted that this choice is mentioned in Refs. [71, 72].

      In order to solve the current system, we substitute Eqs. (11), (31), and (38) into Eqs. (36), (37) and then we have

      $ \begin{aligned}[b] \rho_{\rm de} =& \frac{3 \alpha}{\kappa^2}\left(2 \beta H_0^2 - \delta^2 \frac{\dot{N}^2}{N^2}\right) {\rm e}^{\sigma ({N^2}/{\dot{N}^2})} + \frac{3 \delta^2}{\kappa^2}\frac{\dot{N}^2}{N^2} \\&\times\left(\alpha + 1 + \frac{2 \alpha \beta \sqrt{\beta} H_0}{\delta}\frac{N}{\dot{N}}\right) - \frac{3 H_0^2 \Omega_{m_0}}{\kappa^2} N^{3 \delta \eta}, \end{aligned} $

      (39)

      $ \begin{aligned}[b]{p}_{\rm de} =& \frac{\alpha}{\kappa^2}\Bigg[\delta (3 \delta - 2) \frac{\dot{N}^2}{N^2} - 2 \beta H_0^2 \left(3 - \frac{1}{\delta}\right) \\&+ 2 \delta \frac{\ddot{N}}{N} - \frac{4 \beta^2 H_0^4}{\delta^3} \frac{N^2}{\dot{N}^2} \left(1 - \frac{N \ddot{N}}{\dot{N}^2}\right) \\&- \frac{2 \beta H_0^2}{\delta} \frac{N \ddot{N}}{\dot{N}^2}\Bigg] {\rm e}^{\sigma ({N^2}/{\dot{N}^2})} \\& - \frac{\delta (\alpha + 1)(3 \delta - 2)}{\kappa^2}\frac{\dot{N}^2}{N^2} - 3 \xi \delta \frac{\dot{N}}{N} - \frac{2 \delta (\alpha + 1)}{\kappa^2}\frac{\ddot{N}}{N} \\&- \frac{3 H_0^2 \omega_m \Omega_{m_0}}{\kappa^2} N^{3 \delta \eta} - \frac{2 \alpha H_0}{\kappa^2 \sqrt{\beta}} \left((3 \delta - 1) \frac{\dot{N}}{N} + \frac{\ddot{N}}{\dot{N}}\right), \end{aligned} $

      (40)

      where the dot symbol demonstrates the derivative with respect to time evolution, and we have $ \eta = b^{2} - 1 - \omega_{m} $, $ \delta = (1 + \omega_m)/3 b^2 $, and $ \sigma = \beta H_0^2/\delta^2 $. Here we can see that the Friedmann equations depend on the number of particles, N, and its corresponding derivatives with respect to time evolution.

      Now we intend to reconstruct the aforesaid cosmological quantities in terms of the redshift parameter, z, which is written in terms of the scale factor by the relation $ a(t) = 1/(1+z) $. We can write down the transformation between the derivative with respect to z and the derivative with respect to time evolution as $ {\rm d}/{\rm d}t = -H (1+z) {\rm d}/{\rm d}z $. Therefore, Eq. (31) is written in terms of the redshift parameter as

      $\begin{equation} N = \left(1+z\right)^{-\frac{1}{\delta}}, \end{equation} $

      (41)

      and we write $ \dot{N} $ and $ \ddot{N} $ in terms of the redshift parameter in the following form

      $ \begin{eqnarray} & \dot{N}(t) = - H (1+z) N', \end{eqnarray} $

      (42)

      $ \begin{eqnarray} & \ddot{N}(t) = H H' (1+z)^2 N' + H^2 (1+z) N' + H^2 (1+z)^2 N''(z), \end{eqnarray} $

      (43)

      where the prime symbol indicates the derivative with respect to the redshift parameter. We can see the variation in the number of particles in terms of the redshift parameter in Fig. 1. Figure 1 shows that the value of the number of particles increases from the early universe to the late universe.

      Figure 1.  Number of particles in terms of redshift parameter for $ \alpha = 0.001 $, $ \beta = 0.5 $, $ b = 0.5 $, and $ \omega_m = 0.01 $.

      In what follows, in order to solve the above scenario, we take a particular model called the power-law cosmology model as below:

      $ \begin{equation} a(t) = \left(\frac{t}{t_0}\right)^m, \end{equation} $

      (44)

      where $ t_0 $ is the present age of the universe, and m is a dimensionless positive coefficient [73, 74]. It should be noted that if $ t =t_0 $, we will have $ a(t_0) = a_0 =1 $, which has already been mentioned. The Hubble parameter yields

      $ \begin{equation} H = \frac{m}{t}, \end{equation} $

      (45)

      where the present age of the universe in terms of the present Hubble parameter, $ H_0 $, is written as

      $ \begin{equation} t_0 = \frac{m}{H_0}, \end{equation} $

      (46)

      where the coefficient m is introduced as a correction factor. In this case, the Hubble parameter is obtained in terms of the redshift parameter in the following form:

      $ \begin{equation} H = H_0 (1+z)^{{1}/{m}}, \end{equation} $

      (47)

      where only one free parameter is included, which in turn is less complicated and can be easily calculated with observational constraints.

      Since Eq. (47) is acquired by the power-law model in terms of the redshift parameter, we have to deliberate the correctness of using it in our model. For this purpose, we fit the power-law model with 51 supernova data points from observational Hubble parameter data [7582], and then, by using the likelihood analysis (the chi-square value, $ \chi_{\rm min}^2 $), we obtain the value of $ m = 0.95 $ for these data. In this case, by using Eq. (46), the age of the universe is obtained as $ t_0 = 13.78 \;{\rm Gyr} $ (note that the results presented in this paragraph are the same as those in Ref. [83], and we will not go into detail in writing).

      Now, we can obtain $ \dot{N} $ and $ \ddot{N} $ in terms of the pure redshift parameter in the form

      $ \begin{eqnarray} & \dot{N}(t) = \frac{H_0}{\delta} (1+z)^{({1}/{m})-({1}/{\delta})}, \end{eqnarray} $

      (48)

      $ \begin{eqnarray} & \ddot{N}(t) = \frac{H_0^2 (m - \delta)}{m \,\delta^2}(1+z)^{({2}/{m})-({1}/{\delta})}. \end{eqnarray}$

      (49)

      By inserting Eqs. (48), (49) into Eqs. (39), (40) we have

      $ \begin{aligned}[b] \rho_{\rm de} =& \frac{3 \alpha H_0^2}{\kappa^2}\left(2 \beta - (1+z)^{{2}/{m}}\right) {\rm e}^{\beta (1+z)^{-{2}/{m}}} \\&+ \frac{3 (\alpha+1) H_0^2}{\kappa^2} (1+z)^{{2}/{m}} + \frac{6 \alpha \beta \sqrt{\beta} H_0^2}{\kappa^2} (1+z)^{{1}/{m}}\\& - \frac{3 H_0^2 \Omega_{m_0}}{\kappa^2} (1+z)^{-3 \eta}, \end{aligned} $

      (50)

      $ \begin{aligned}[b] {p}_{\rm de} =& \frac{\alpha H_0^2}{\kappa^2 m} \Big(-4 \beta^2 \left(1+z \right)^{-{2}/{m}} + (3 m - 2) \left(1+z \right)^{{2}/{m}} \\&- 2 \beta (3 m - 1)\Big) {\rm e}^{\beta \left(1+z \right)^{-{2}/{m}}}\\& - \frac{3 H_0^2 \omega_m \Omega_{m_0}}{\kappa^2} \left(1+z \right)^{-3 \eta} - \frac{H_0}{\kappa^2 m \sqrt{\beta}} \Big(\kappa^2 m \sqrt{\beta}\, \xi \\& + 2 \alpha H_0 \left(3m - 1\right)\Big) \left(1+z \right)^{{1}/{m}}\\& - \frac{\left(3 m - 2\right) \left(\alpha +1\right) H_0^2}{\kappa^{2}m } \left(1+z \right)^{{2}/{m}}, \end{aligned} $

      (51)

      which shows that $ \rho_{\rm de} $ and $ p_{\rm de} $ are explicit functions of the redshift parameter. We note that these obtained relationships, i.e., those of the energy density and pressure of dark energy, depend on a few free parameters that arise from scenarios of the $ f(T) $gravity, matter, particle number, and viscous fluid.

      The variations of relationships $ \rho_{\rm de} $ and $ p_{\rm de} $ are plotted in terms of the redshift parameter as shown in Fig. 2. It should be noted that free parameters play an important role in drawing the corresponding graphs. For this purpose, we try to choose them by applying conditions such as $ \rho_{\rm de} > 0 $ and $ p_{\rm de} < 0 $, especially in the late time ($ z = 0 $). However, we consider the values of the free parameters for the corresponding graphs as $ \alpha = 0.001 $, $ \beta = 0.5 $, $ b = 0.5 $, $ \omega_m = 0.01 $, and $ \xi = 26, 30, 34 $. Fig. 2 shows that the energy density decreases from a very high positive value in the early time to a much lower positive value in the late time ($ z = 0 $), and the pressure decreases from a very large negative value in the early time to a much lower negative value in the late time ($ z = 0 $).

      Figure 2.  Energy density and pressure of dark energy in terms of redshift parameter for $ \alpha = 0.001 $, $ \beta = 0.5 $, $ b = 0.5 $, $ \omega_m = 0.01 $, $ \xi = 26 $ (line), $ \xi = 30 $ (dot), and $ \xi = 34 $ (dash).

      Since the dark energy EoS is one of the most practical cosmic parameters, it can describe the nature of dark energy and dark matter for different eras of the universe from the Big Bang to the late time. Now, using $ \rho_{\rm de} $ and $ p_{\rm de} $, we clearly find the dark energy EoS as $ \omega_{\rm de} = p_{\rm de}/\rho_{\rm de} $ by Eqs. (50), (51) (we avoid writing the corresponding EoS because it is long). In that case, we can plot the variation in the dark energy EoS, $ \omega_{\rm de} $, in terms of the redshift parameter as shown in Fig. 3. Figure 3 shows that the amount of EoS decreases from early time to late time, i.e., the evolution of the universe changes from the matter era ($ \omega_{\rm de} > -1 $) to the phantom era or dark energy era ($ \omega_{\rm de} < -1 $). Interestingly, the amount of EoS depends on the effect of the viscous fluid; hence, from Fig. 3, we can see that, when the amount of the viscosity coefficient increases, the dark energy period occurs with more growth. Therefore, this result indicates that the universe is undergoing an accelerated expansion phase, which is compatible with the obtained results in Ref. [84] for EoS.

      Figure 3.  EoS of dark energy in terms of redshift parameter for $ \alpha = 0.001 $, $ \beta = 0.5 $, $ b = 0.5 $, $ \omega_m = 0.01 $, $ \xi = 26 $ (line), $ \xi = 30 $ (dot), and $ \xi = 34 $ (dash).

    VI.   STABILITY ANALYSIS
    • In this section, we study the stability and instability of the present model from the thermodynamic perspective. In this case, we employ a useful function called the sound speed parameter, $ c_s^2 $. Since the universe is an adiabatic thermodynamic system, heat energy or mass is not transferred from the inside to outside of the universe. This means that the entropy perturbation is zero, and the pressure changes in terms of the energy density. Therefore, we introduce the pressure-to-density ratio with the below sound speed parameter:

      $ \begin{equation} c_s^2=\frac{\partial p_{\rm de}}{\partial \rho_{\rm de}} = \frac{\partial_z p_{\rm de}}{\partial_z \rho_{\rm de}}, \end{equation} $

      (52)

      where symbol $ \partial_z $ represents the derivative with respect to the redshift parameter. Now, the parameter $ c_s^2 $ helps us find the range of the redshift parameter wherein the universe is in a state of stability or instability. It should be noted that the stability and instability analyses are expressed with conditions $ c_s^2 > 0 $ and $ c_s^2 < 0 $, respectively. The variation in the sound speed parameter is plotted in terms of the redshift parameter as shown in Fig. 4. Figure 4 shows that the universe is in a phase of instability in all eras, even with the existence of viscous fluid. This instability indicates that the energy density of dark energy does not exhibit a controllable growth.

      Figure 4.  Sound speed parameter in terms of redshift parameter for $ \alpha = 0.001 $, $ \beta = 0.5 $, $ b = 0.5 $, $ \omega_m = 0.01 $, $ \xi = 26 $ (line), $ \xi = 30 $ (dot), and $ \xi = 34 $ (dash).

      Now, for a more complete description of our model, we attempt to calculate the present values of the density parameters for dark energy and dark matter. According to the density parameter for dark matter (33), the density parameter for dark energy is introduced in the form

      $ \begin{equation} \Omega_{\rm de} = \frac{\rho_{\rm de}}{\rho_{\rm c}}, \end{equation} $

      (53)

      where $ \rho_{\rm c} $ is the critical density. By inserting Eq. (35) into Eq. (33), the variation in the density parameter for dark matter is plotted in terms of the redshift parameter as shown in Fig. 5(right side). Similarly, we plot the variation in the density parameter for dark energy in terms of the redshift parameter by substituting Eq. (51) into Eq. (53) as shown in Fig. 5(left side). Figure 5 shows that the value of $ \Omega_m $ decreases from the early universe to the late universe, eventually reaching the current value of $ 0.315 $. In contrast, the value of $ \Omega_{\rm de} $ increases from the early time to the late time, reaching the present value of $ 0.689 $. Therefore, these calculations show that the proportion of dark energy is increasing, which means that the universe is under a phase of accelerated expansion.

      Figure 5.  Density parameters of dark energy and dark matter in terms of redshift parameter for $ \alpha = 0.001 $, $ \beta = 0.5 $, $ b = 0.5 $, and $ \omega_m = 0.01 $.

    VII.   CONCLUSION
    • In this study, we explored particle creation from the perspective of $ f(T) $ gravity using the flat-FRW metric in the universe, which contains a viscous fluid. First, to obtain the Einstein equation in $ f(T) $ gravity, we used the Weitzenböck connection in teleparallel gravity instead of the Levi-Civita connection in general relativity. The energy-momentum tensor was considered in terms of the total energy density, total pressure, and pressure of bulk viscosity of fluid within the universe. Next, we acquired the Friedmann equations and EoS of dark energy in terms of a function of torsion scalar T and its first and second derivatives. Then, the continuity equations for the universe components were obtained by considering the interaction model. The interaction term $ \mathcal{Q} $ was considered as $ \mathcal{Q} = 3 b^2 H \rho $, in which $ b^2 $ is the transfer strength from dark matter to dark energy.

      Another significant aspect of this paper is that the universe should be considered as an adiabatic system from the perspective of the first law of thermodynamics. For this purpose, we considered the creation of the number of particles in this universe, and thus, we obtained the modified continuity equation for particle creation. An interesting point is that the obtained continuity equation for matter creation has a negative pressure, which causes the present accelerated expansion. Next, we acquired the energy density of matter in terms of the number of particles, N, as shown in Eq. (28). Another interesting point is that Eq. (27) is equivalent to Eq. (22). This means that we consider particle creation to be equivalent to the exchange of energy between the components of the universe (same as the interaction term $ \mathcal{Q}) $. The result of this equivalency leads to obtaining the number of particles and energy density of matter in terms of the scale factor.

      Subsequently, we related the obtained energy density of matter to Friedmann equations of $ f(T) $ gravity with a specified function $ f(T) $ of Eq. (38), which is exponential and logarithmic. In this case, we obtained the energy density (39) and pressure (40) of dark energy in terms of the number of particles N. Then, we reconstructed the number of particles in terms of the redshift parameter; the benefit of this reconstruction is that the cosmological parameters can be evaluated for late time ($ z = 0 $). For evaluations in this study, we took the power-law for the scale factor, and thus, we immediately acquired the Hubble parameter (47) in terms of the redshift parameter. After that, we used the result of Ref. [83] to fit the Hubble parameter by the minimum chi-square value with $ 51 $ observational Hubble parameter data points. Then, we plotted the graphs of the cosmological parameters such as the energy density, pressure, and EoS of dark energy in terms of the redshift parameter. These graphs indicated that the universe is in a phase of accelerated expansion. Further, stability analysis of the present model was investigated using sound speed, and the corresponding graph showed that the universe is in an instability phase during late time, i.e., the growth of the energy density of dark energy is not controlled in the present universe. Finally, in order for the results to be consistent with the observational data, we calculated the density parameter values for dark energy and dark matter, which were $ 0.689 $ and $ 0.315 $, respectively, and also plotted them against the redshift parameter.

Reference (84)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return