Generalized Tachyonic Teleparallel cosmology

TL;DR

This paper introduces a new dark energy model within teleparallel gravity, utilizing a generalized non-minimal coupling of a tachyonic scalar field, and analyzes its stability and cosmological implications through phase space analysis.

Contribution

It proposes a novel tachyonic dark energy model with generalized boundary coupling in teleparallel gravity, exploring its dynamical properties and critical points.

Findings

Model exhibits accelerated expansion and scaling solutions.

Phase space analysis reveals critical points related to different cosmological behaviors.

Conditions for the existence of these solutions are identified.

Abstract

In this paper we propose a new dark energy model in the teleparallel alternative of general relativity, by considering a generalized non--minimal coupling of a tachyonic scalar field with the teleparallel boundary term. Within the framework of teleparallel gravity, the boundary coupling term is associated with the divergence of the torsion vector. Considering the linear stability technique for various potentials and couplings, we have analyzed the dynamical properties of the present tachyonic dark energy model in the phase space, uncovering the corresponding essential dynamical features. Our study of the phase space structure revealed that for a specific class of potential energy, this model exhibits various critical points which are related to different cosmological behaviors, such as accelerated expansion and scaling solutions, determining the existence conditions and the corresponding physical features.

Figures

Tables

Table 1. Table 1: The critical points in the case of exponential coupling and inverse hyperbolic sine potential. The auxiliary variable Ξ Ξ \Xi used in the description of the eigenvalues for the various critical points is equal to: Ξ = ω 1 4 ω 2 4 ( − w m ) ( w m + 1 ) ( 17 ω 1 2 ω 2 2 w m 2 + 14 ω 1 2 ω 2 2 w m − 96 ( − w m ) + 3 / 2 48 ( − w m ) + 5 / 2 48 − w m + ω 1 2 ω 2 2 ) ω 1 3 ​ ω 2 3 ​ − w m ​ ( w m + 1 ) \Xi=\frac{\sqrt{\omega_{1}^{4}\omega_{2}^{4}\left(-w_{\rm m}\right)\left(w_{\rm m}+1\right)\left(17\omega_{1}^{2}\omega_{2}^{2}w_{\rm m}^{2}+14\omega_{1}^{2}\omega_{2}^{2}w_{\rm m}-96\left(-w_{\rm m}\right){}^{3/2}+48\left(-w_{\rm m}\right){}^{5/2}+48\sqrt{-w_{\rm m}}+\omega_{1}^{2}\omega_{2}^{2}\right)}}{\omega_{1}^{3}\omega_{2}^{3}\sqrt{-w_{\rm m}\left(w_{\rm m}+1\right)}}

Point	x	y	u	$\lambda$	${\Omega }_{\varphi }$	$w_{\mathrm{eff}}$	Eigenvalues
A	0	1	u	6u	1	-1	$0,-3\left(w_m+1\right),-\frac{3}{2}\pm \frac{\sqrt{3}\sqrt{\left(6u^2+1\right){\omega }_1\left({\omega }_1\left(18u^2-24\alpha u+3\right)-144u^2+4{\omega }_1^2{\omega }_2^2\right)}}{2\left(6u^2+1\right){\omega }_1}$
$B_{+}^{-}$	$\sqrt{w_m+1}$	$-\frac{\sqrt{3}\sqrt{w_m+1}}{{\omega }_1{\omega }_2}$	0	$-{\omega }_1{\omega }_2$	$\frac{3\left(w_m+1\right)}{{\omega }_1^2{\omega }_2^2\sqrt{-w_m}}$	$w_m$	$-\frac{3\alpha \left(w_m+1\right)}{{\omega }_1{\omega }_2},-\frac{6\left(w_m+1\right)}{{\omega }_1},\frac{3}{4}\left(w_m-1\pm \Xi \right)$
$B_{-}^{+}$	$-\sqrt{w_m+1}$	$\frac{\sqrt{3}\sqrt{w_m+1}}{{\omega }_1{\omega }_2}$	0	$-{\omega }_1{\omega }_2$	$\frac{3\left(w_m+1\right)}{{\omega }_1^2{\omega }_2^2\sqrt{-w_m}}$	$w_m$	$-\frac{3\alpha \left(w_m+1\right)}{{\omega }_1{\omega }_2},-\frac{6\left(w_m+1\right)}{{\omega }_1},\frac{3}{4}\left(w_m-1\pm \Xi \right)$
$B_{-}^{-}$	$-\sqrt{w_m+1}$	$-\frac{\sqrt{3}\sqrt{w_m+1}}{{\omega }_1{\omega }_2}$	0	${\omega }_1{\omega }_2$	$\frac{3\left(w_m+1\right)}{{\omega }_1^2{\omega }_2^2\sqrt{-w_m}}$	$w_m$	$\frac{3\alpha \left(w_m+1\right)}{{\omega }_1{\omega }_2},-\frac{6\left(w_m+1\right)}{{\omega }_1},\frac{3}{4}\left(w_m-1\pm \Xi \right)$
$B_{+}^{+}$	$\sqrt{w_m+1}$	$\frac{\sqrt{3}\sqrt{w_m+1}}{{\omega }_1{\omega }_2}$	0	${\omega }_1{\omega }_2$	$\frac{3\left(w_m+1\right)}{{\omega }_1^2{\omega }_2^2\sqrt{-w_m}}$	$w_m$	$\frac{3\alpha \left(w_m+1\right)}{{\omega }_1{\omega }_2},-\frac{6\left(w_m+1\right)}{{\omega }_1},\frac{3}{4}\left(w_m-1\pm \Xi \right)$

Equations

$$R=-T+B,$$

$$S=\int \left[\frac{T}{2{\kappa }^2}+\frac{1}{2}f(\varphi )T+\frac{1}{2}g(\varphi )B-V(\varphi )\sqrt{1-\frac{2X}{V(\varphi )}}+L_{\mathrm{m}}\right]e{d}^{4}x,$$

$$X=\frac{1}{2}({\partial }_{\mu }\varphi )({\partial }^{\mu }\varphi )$$

$$T_{\mu \nu }^a={\partial }_{\mu }{e}_{\nu }^a-{\partial }_{\nu }{e}_{\mu }^a+{\omega }_{b\mu }^a{e}_{\nu }^b-{\omega }_{b\nu }^a{e}_{\mu }^b,$$

$$T=T_{\mu \nu }^a{S}_a^{\mu \nu },$$

$$S_a^{\mu \nu }=\frac{1}{2}\left(T_a^{\mu \nu }+T_a^{\nu \mu }-T_a^{\mu \nu }\right)-E_{\nu }^a{T}_{\alpha }^{\alpha \mu }+E_{\mu }^a{T}_{\alpha }^{\alpha \nu }.$$

$$\displaystyle 2\Big{(}\frac{1}{\kappa^{2}}+f(\phi)\Big{)}\left[e^{-1}\partial_{\mu}(eS_{a}{}^{\mu\nu})-E_{a}^{\lambda}T^{\rho}{}_{\mu\lambda}S_{\rho}{}^{\nu\mu}-\frac{1}{4}E^{\nu}_{a}T\right]-E_{a}^{\nu}V(\phi)\sqrt{1-\frac{2X}{V(\phi)}}$$

$$-\frac{1}{\sqrt{1-\frac{2X}{V(\varphi )}}}{E}_a^{\nu }{\partial }_{\mu }\varphi {\partial }^{\mu }\varphi +2({\partial }_{\mu }f(\varphi )+{\partial }_{\mu }g(\varphi ))E_a^{\rho }{S}_{\rho }{}^{\mu \nu }+E_a^{\nu }\square g(\varphi )-E_a^{\mu }{\nabla }^{\nu }{\nabla }_{\mu }g(\varphi )={\mathcal{T}}_a^{\nu }.$$

$$\displaystyle\frac{V^{\prime}(\phi)}{2V(\phi)\sqrt{1-\frac{2X}{V(\phi)}}}(\partial_{\mu}\phi)(\partial^{\mu}\phi)+\partial_{\mu}\Big{(}\frac{\partial^{\mu}\phi}{\sqrt{1-\frac{2X}{V(\phi)}}}\Big{)}+V^{\prime}(\phi)\sqrt{1-\frac{2X}{V(\phi)}}$$

$$d{s}^2=d{t}^2-a(t{)}^2(d{x}^2+d{y}^2+d{z}^2),$$

$$e_{\mu }^a=\text{diag}(1,a(t),a(t),a(t)),$$

$$3H^2$$

$$3H^2+2\dot{H}$$

$${\dot{\rho }}_{\mathrm{m}}+3H({\rho }_{\mathrm{m}}+p_{\mathrm{m}})=0,$$

$$\ddot{\varphi }+3H\dot{\varphi }\left(1-\frac{{\dot{\varphi }}^2}{V(\varphi )}\right)+\left(1-\frac{3{\dot{\varphi }}^2}{2V(\varphi )}\right)V^{\prime }(\varphi )$$

$$3H^2={\rho }_{\mathrm{eff}},3H^2+2\dot{H}=-p_{\mathrm{eff}},$$

$${\rho }_{\varphi }$$

$$p_{\varphi }$$

$$w_{\mathrm{eff}}=\frac{p_{\mathrm{eff}}}{{\rho }_{\mathrm{eff}}}=\frac{p_{\mathrm{m}}+p_{\varphi }}{{\rho }_{\mathrm{m}}+{\rho }_{\varphi }},$$

$${\Omega }_{\mathrm{m}}=\frac{{\rho }_{\mathrm{m}}}{3H^2},{\Omega }_{\varphi }=\frac{{\rho }_{\varphi }}{3H^2},$$

$$x=\frac{\dot{\varphi }}{\sqrt{V(\varphi )}},y=\frac{\sqrt{V(\varphi )}}{\sqrt{3}H},u=\frac{1}{2}{g}^{\prime }(\varphi ),\lambda =-\frac{V^{\prime }(\varphi )}{V(\varphi )},$$

$$0\le {\Omega }_{\mathrm{m}}=1-2\sqrt{3}uxy-\frac{y^2}{\sqrt{1-x^2}}\le 1,$$

$$g(\varphi )=\chi {\varphi }^p,$$

$$\Gamma =\frac{V(\varphi )V^{\prime \prime }(\varphi )}{V^{\prime }(\varphi )V^{\prime }(\varphi )}.$$

$$\frac{dx}{dN}=\frac{\ddot{\varphi }}{\sqrt{3}{H}^{2}y}+\frac{1}{2}\sqrt{3}\lambda x^{2}y,$$

$$\frac{dy}{dN}=-\frac{y\dot{H}}{H^2}-\frac{1}{2}\sqrt{3}\lambda x{y}^2,$$

$$\frac{du}{dN}=\sqrt{3}\cdot 2^{\frac{1}{1-p}}(p-1)p\chi xy{\left(\frac{u}{p\chi }\right)}^{\frac{p-2}{p-1}},$$

$$\frac{d\lambda }{dN}=-\sqrt{3}(\Gamma -1){\lambda }^{2}xy,$$

$$\ddot{\varphi }=3H^2\lambda \left(1-\frac{3x^2}{2}\right)y^2-3\sqrt{3}{H}^{2}x\left(1-x^2\right)y-6u{\left(1-x^2\right)}^{3/2}\left(3H^2+\dot{H}\right),$$

$$\dot{H}=\frac{3H^{2}}{12u^{2}\left(1-x^{2}\right)^{3/2}+2}\cdot\Bigg{[}w_{m}\left(2\sqrt{3}uxy+\frac{y^{2}}{\sqrt{1-x^{2}}}-1\right)+2^{\frac{p-2}{p-1}}(p-1)p\chi x^{2}y^{2}\left(\frac{u}{p\chi}\right)^{\frac{p-2}{p-1}}\\ -12u^{2}\left(1-x^{2}\right)^{3/2}-\lambda u\left(3x^{2}-2\right)y^{2}+2\sqrt{3}ux\left(x^{2}-1\right)y+\sqrt{1-x^{2}}y^{2}-1\Bigg{]}.$$

$$\frac{dx}{dN}=\frac{1}{y\left(6u^{2}\left(x^{2}-1\right)^{2}+\sqrt{1-x^{2}}\right)}\Big{(}-18u^{2}x^{5}yw_{\rm m}+36u^{2}x^{3}yw_{\rm m}-18u^{2}xyw_{\rm m}+3\sqrt{3}ux^{4}w_{\rm m}+3u\sqrt{3-3x^{2}}x^{2}y^{2}w_{\rm m}\\ -3u\sqrt{3-3x^{2}}y^{2}w_{\rm m}-6\sqrt{3}ux^{2}w_{\rm m}+3\sqrt{3}uw_{\rm m}+3\sqrt{3}\lambda u^{2}x^{6}y^{2}-6\sqrt{3}\lambda u^{2}x^{4}y^{2}+3\sqrt{3}\lambda u^{2}x^{2}y^{2}\\ -6\sqrt{3}u\chi x^{6}y^{2}+12\sqrt{3}u\chi x^{4}y^{2}-3\sqrt{3}ux^{4}-6\sqrt{3}u\chi x^{2}y^{2}+6\sqrt{3}u\sqrt{1-x^{2}}x^{2}y^{2}-3\sqrt{3}u\sqrt{1-x^{2}}y^{2}\\ +6\sqrt{3}ux^{2}-3\sqrt{3}u\sqrt{1-x^{2}}x^{4}y^{2}-3\sqrt{3}u-\sqrt{3}\lambda\sqrt{1-x^{2}}x^{2}y^{2}+\sqrt{3}\lambda\sqrt{1-x^{2}}y^{2}-3\sqrt{1-x^{2}}xy+3\sqrt{1-x^{2}}x^{3}y\Big{)}\,,$$

$$\frac{dy}{dN}=\frac{1}{2\left(6u^{2}\left(x^{2}-1\right)^{2}+\sqrt{1-x^{2}}\right)}\Big{(}-6u\sqrt{3-3x^{2}}xy^{2}w_{\rm m}+3\sqrt{1-x^{2}}yw_{\rm m}-3y^{3}w_{\rm m}-6\sqrt{3}\lambda u^{2}x^{5}y^{2}+36u^{2}x^{4}y\\ +12\sqrt{3}\lambda u^{2}x^{3}y^{2}-72u^{2}x^{2}y-6\sqrt{3}\lambda u^{2}xy^{2}+36u^{2}y+9\lambda u\sqrt{1-x^{2}}x^{2}y^{3}-6\lambda u\sqrt{1-x^{2}}y^{3}\\ +6\sqrt{3}u\sqrt{1-x^{2}}xy^{2}-6\sqrt{3}u\sqrt{1-x^{2}}x^{3}y^{2}-6\chi\sqrt{1-x^{2}}x^{2}y^{3}+3x^{2}y^{3}-\lambda\sqrt{3-3x^{2}}xy^{2}+3\sqrt{1-x^{2}}y-3y^{3}\Big{)}\,,$$

Full text

I Introduction

Einstein’s theory of General Relativity (GR) has experienced unprecedented success in its power to explain astrophysical phenomena ranging from Solar System tests to strong field gravitational wave physics Misner:1974qy . However, this theory of gravity has required important modifications due to observational realities which have arisen over the past few decades. In terms of the energy budget of the Universe, the first modification comes from observations of galaxies and their dynamical structure, which is only possible with the addition of approximately purely gravitational interacting particles, namely dark matter, which may potentially be beyond the standard model of particle physics Rubin:1970zza ; Navarro:1995iw . The second and larger contribution to the modification of GR comes from the relatively recent observation of the accelerating expansion of the Universe Riess:1998cb ; Perlmutter:1998np which is an observational fact, called dark energy. This can be accounted for in GR through the introduction of the cosmological constant however this poses its own problems Weinberg:1988cp . The $\Lambda$CDM is the most successful model which can explain the current accelerated expansion and the evolution of the observable Universe at the level of background dynamics, involving the superposition between the dark matter fluid and the cosmological constant. On the other hand, the early period of the Universe also features several facets that need remedy. Most prominently, for $\Lambda$CDM to correctly produce our current picture of the Universe, a period of cosmological inflation must of taken place Guth:1980zm ; Linde:1981mu which would allow for a natural solution to the horizon problem. However, this may also necessitate further particles beyond the standard model Lyth:1998xn . Cosmologically, the time that should be best described by the $\Lambda$CDM model is the present or late-time period of the Universe. However, recent releases by the Planck collaboration have revealed a growing tension in the local and global measurements of $H_{0}$ and $f\sigma_{8}$ Aghanim:2018eyx .

Now, it may be the case that the fundamental and observational problems surrounding $\Lambda$CDM may be resolved in the coming years, or it may be the case that $\Lambda$CDM needs to be changed in some way. Over the previous decades there has been concerted efforts in extending GR to account for certain elements of these problems Capozziello:2011et . However, it may also be the case that a new paradigm is needed to confront the growing requirements of constructing a viable theory of gravity. One such treatment is the teleparallel gravity approach where the Levi-Civita connection is replaced with the Weitzenböck connection AP . The connection plays a crucial role in gravitational physics in that the expression of curvature, torsion, or nonmetricity is not a property of the manifold itself but of the connection which relates the elemental tangent spaces of the manifold Heisenberg:2018vsk . In this way, one can choose to consider gravitation in terms of the curvatureless Weitzenböck connection which also observes the metricity property.

In teleparallel gravity, the gravitation is characterized by the torsion tensor, $T^{\rho}_{\phantom{\rho}\mu\nu}$, instead of the Riemann tensor in GR and its extensions. As in the GR framework, a Lagrangian can be constructed to represent the gravitational field. Of particular interest, in teleparallel gravity, is that a Lagrangian can be constructed such that it is equivalent to the Einstein-Hilbert Lagrangian up to a total divergence or boundary term, $B$, that is

[TABLE]

where $T$ is called the torsion scalar and contains only second order terms, while the boundary term, $B$, encapsulates the higher order contributions to the Ricci scalar, $R$. This is the so-called teleparallel equivalent of general relativity (TEGR) which is equivalent to GR at the level of the field equations Cai:2015emx ; Krssak:2018ywd . The natural consequence of this realization is that every test of GR also becomes a test of TEGR with the difference that in TEGR gravity is again observed to act as a (Lorentz) force, and that the barrier with the quantum regime seems to have less tension AP . Moreover, due to the second order nature of the torsion tensor, even in the case where the theory is extended to an $f(T)$ Lagrangian, the resulting field equations remain second order which has important consequences for the gravitational wave polarization modes of the theory. In fact, $f(T)$ gravity continues to exhibit the equivalent polarization modes as in the GR and TEGR settings Farrugia:2018gyz .

Given the decomposition of the Ricci scalar into the second and fourth order terms expressed in Eq. (1), we consider the analysis of a tachyonic dark energy model nonminimally coupled to the aforementioned separate contributions through different functionals. In particular, we choose to study a tachyonic scalar field which has been shown to produce an inflationary epoch and late-time accelerating solutions that do not violate the strong energy condition Quiros:2009mz ; Aguirregabiria:2004xd . These models are partially inspired by string theory Bagla:2002yn ; Fang:2010zze ; Ortin:2015hya and $k$-essence theory Almeida:2016ixq . In the teleparallel setting, scalar fields have been investigated to a moderate degree, with various extensions having been investigated in the cosmological context Geng:2011aj ; Xu:2012jf ; Hohmann:2018rwf ; PhysRevD.97.084008 ; Otalora:2013tba ; Otalora:2014aoa . Furthermore, the effects of the boundary couplings in scalar tensor theories have been discussed in various papers Hohmann:2018ijr ; Bahamonde:2018miw ; Abedi:2018lkr ; Hohmann:2018rwf ; Bahamonde:2017wwk ; Marciu:2017sji ; Gecim:2017hmn ; Bahamonde:2017ifa ; Abedi:2017jqx ; Fazlpour:2018jzz ; Bahamonde:2016grb ; Bahamonde:2016cul ; MohseniSadjadi:2016ukp ; Bahamonde:2015zma . A recent review on various studies related to dynamical analysis in different cosmological constructions can be found in BAHAMONDE20181 . In Ref. Otalora:2013dsa , a tachyonic field is investigated for the modified teleparallel setting, late-time accelerating attractor solutions are found with a field equation of state that realistically tends to the current dark energy value. This approach was extended to more general models in Ref. Fazlpour:2014qaa where the attractor solution context is further clarified.

The work is divided as follows: in Sec. II the tachyonic approach to extended teleparallel theories of gravity is introduced with a focus on the cosmological consequences of the treatment. In Sec. III, the dynamical analysis of the system is undertaken for specific choices of tachyonic field. Finally, in Sec. IV the conclusions are summarized and discussed. Unless stated otherwise geometric units are used throughout the paper. In addition, $e^{a}_{\mu}$ and $E_{a}^{\mu}$ represent the tetrads and the inverse of the tetrads respectively and the $(+---)$ metric signature is used.

II Generalized Tachyonic Teleparallel theories of gravity

In this paper, we present a new teleparallel tachyonic model based on the following action

[TABLE]

where $\kappa^{2}=8\pi G$, $L_{\rm m}$ is a matter Lagrangian, $T$ is the scalar torsion, $B=(2/e)\partial_{\mu}(eT^{\mu})$ is the boundary term, $f(\phi)$ and $g(\phi)$ are scalar field dependent coupling functions, $V(\phi)$ is the potential and

[TABLE]

is the kinetic term. The torsion tensor is identified as

[TABLE]

where $e^{a}_{\nu}$ form a tetrad field of the gravitational system and represent coordinate transformations between the general manifold and the tangent space at any point, while $\omega^{a}_{\phantom{a}b\mu}$ form the spin connection components which are purely inertial and sustain the local Lorentz invariance of the theory Krssak:2015oua ; Cai:2015emx . The torsion scalar is then defined through the contraction

[TABLE]

where the superpotential is defined as

[TABLE]

This work considers an analogous generalization of other Tachyonic models studied in the literature Otalora:2013dsa . By taking, $f(\phi)=-g(\phi)$, one recovers a Tachyonic theory with a non-minimally coupling between the scalar field and the Ricci scalar $R$ due to the relation in Eq. (1). By taking $g(\phi)=0$, one recovers a teleparallel Tachyonic theory where the torsion scalar $T$ is non-minimally coupled with the scalar field. In Banijamali:2012nx , the authors found that this coupling allows the crossing of the phantom divide line. The new coupling between the scalar field and the boundary term is motivated from the scalar tensor theory studied in Bahamonde:2015hza ; Zubair:2016uhx , where the authors found that, without fine-tunning, the system evolves to a late-time accelerating attractor solution.

By taking variation with respect to the tetrad, we find the corresponding gravitational field equations given by

[TABLE]

Here, $\nabla_{\mu}$ is the covariant derivative with respect to the Levi-Civita connection and $\Box=\nabla^{\mu}\nabla_{\mu}$. By taking variations with respect to the scalar field, one finds

[TABLE]

For the flat FLRW metric in Cartesian coordinates, the metric is given by

[TABLE]

where $a(t)$ is the cosmological scale factor. This can equivalently be described by the tetrad field

[TABLE]

which naturally induce vanishing spin connection components Krssak:2015oua ; Cai:2015emx . The ensuing Friedmann equations then turn out to be represented by

[TABLE]

where $\kappa=1$ was assumed, dots represent differentiation with respect to cosmic time and primes with respect to the scalar field. It should be noted that a standard perfect fluid with energy density $\rho_{\rm m}$ and pressure $p_{\rm m}$ is being assumed. It is straightforward to show that the standard conservation of the energy-momentum tensor gives

[TABLE]

which is the standard conservation equation for matter.

The scalar field relation in Eq. (8) takes the following form

[TABLE]

It can be shown that this equation can be also found directly from the modified FLRW equations (11)–(12), therefore, it is not an extra constraint equation.

The Friedmann equations can also equivalently be rewritten as their GR analogue with an additional effective fluid component so that

[TABLE]

where $\rho_{\rm eff}=\rho_{\rm m}+\rho_{\phi}$ and $p_{\rm eff}=p_{\rm m}+p_{\phi}$ with the new fluid quantities defined as follows

[TABLE]

It is then convenient to introduce the effective equation of state parameter

[TABLE]

which means that the density parameters of the contributing components take the form

[TABLE]

in such a way that $\Omega_{\rm m}+\Omega_{\phi}=1$ holds.

III Dynamical system of the model

In this section, we will study the dynamical system of our model which only has non-minimally couplings between the scalar field and the boundary terms, therefore we will assume that $f(\phi)=0$, and $g(\phi)\neq 0$. The dynamical system for the case $f(\phi)\neq 0$ and $g(\phi)=0$ was studied previously in Otalora:2013dsa ; Fazlpour:2014qaa . Let us first introduce the following dimensionless variables

[TABLE]

so that $y>0$. By using these variables in the first FLRW equation in Eq. (11) and the definition of the energy density parameter given by Eq. (19), one gets the constraint

[TABLE]

which gives the phase space of the dynamical system. One needs to choose $g(\phi)$ to write down the full dynamical system of the model. Therefore, in the next sections, we will study two different kind of couplings: power-law and exponential types. We will further assume standard a barotropic fluid given by $p_{\rm m}=w_{\rm m}\rho_{\rm m}$.

III.1 Power-law coupling and exponential potential

Assuming a power-law coupling between the boundary term and the scalar field means setting

[TABLE]

where $\chi$ and $p$ are both constants. It is possible to write down the dynamical equations of the system as a 4-dimensional one with a generic potential. This can be done if one introduces the variable

[TABLE]

By introducing the variable $N=\log(a)$, the dynamical system for this kind of couplings can be written as

[TABLE]

[TABLE]

[TABLE]

[TABLE]

where

[TABLE]

and

[TABLE]

This dynamical system assumes that $p\neq 1$. The case $p=1$ is a very special one which gives a 3-dimensional dynamical system of equations, but generically, does not give any interesting cosmological behaviour. One more realistic model is the one where $p=2$ which gives a coupling like $(1/2)\chi\phi^{2}B$ which was studied in Bahamonde:2015hza . In that case however, non-tachyonic scalar fields were considered. In what follows, we shall consider the case where we have a specific coupling, $g(\phi)=\chi\phi^{p}$, with $p=2$ and an exponential potential $V(\phi)=V_{0}e^{-\lambda\phi}$, where $\lambda$ is a positive constant. In this case, the dynamical evolution of the model can be described by a 3D autonomous system of differential equations

[TABLE]

[TABLE]

[TABLE]

In this scenario, the effective equation of state for the dark energy field can be written as:

[TABLE]

Next, the critical points are determined by considering that the RHS of Eqs. (25)–(27) are equal to zero, taking into account also the physical viability which requires: $0\leq\Omega_{m}=1-\Omega_{\phi}\leq 1$, $y\geq 0$, $1-x^{2}>0$ and the location of the corresponding critical points is the the real space. In the case of this specific power law coupling and for an exponential potential, the phase space structure has a reduced complexity, having only one critical point located at $O(x,y,u)=\left(0,1,\frac{\lambda}{6}\right)$, with the following eigenvalues

[TABLE]

This critical point corresponds to a de Sitter evolution, acting as a cosmological constant $w_{\rm eff}=-1$, implying the full domination of the tachyonic dark energy field over the matter component, $\Omega_{m}=1-\Omega_{\phi}=0$. For this critical point, we show in Fig. 1, various regions for the model’s parameters which are connected to different cosmological scenarios, corresponding to a stable and stable spiral evolution, respectively.

III.2 Exponential coupling and inverse hyperbolic sine potential

Here, we investigate the case where the coupling functional is represented by an exponential, $g(\phi)=g_{0}e^{\alpha\phi}$, with $g_{0}$ and $\alpha$ constants. The potential energy associated with the present tachyonic dark energy model is considered to be an inverse hyperbolic sine, $V(\phi)=V_{0}\sinh^{-\omega_{1}}(\omega_{2}\phi)$, with $\omega_{1,2}$ constants. This type of potential is beyond the usual exponential type considered in many dynamical analysis and is motivated by the recent work in Ref. Roy:2017uvr , where the structure of the phase space for non–canonical fields with the potential energy beyond exponential type was investigated. The potential energy considered in this section, the inverse hyperbolic sine UrenaLopez:2000aj ; Sahni:1999gb represents a possible parameterization for the dark energy field which is associated to a second order polynomial in the dynamical equation for the dimensionless variable $\lambda$. As discussed in Ref. Roy:2017uvr , if we denote $f(\lambda)=\lambda^{2}(\Gamma-1)$ with $\Gamma$ defined in the relation in Eq. (23), then for this specific potential, the function $f(\lambda)$ obeys a second order polynomial parameterization $f(\lambda)=a\lambda^{2}+b\lambda+c$, with $a,b,c$ constant parameters. For the specific potential considered here, the inverse hyperbolic sine, we find that $f(\lambda)=\lambda^{2}/\omega_{1}-\omega_{1}\omega_{2}^{2}$.

In this case, the evolution of the dynamical system can be written as the following system of differential equations

[TABLE]

[TABLE]

[TABLE]

and finally the last equation for the dynamical system given by

[TABLE]

The total effective equation of state is:

[TABLE]

Next, as in the previous case, the critical points are obtained by considering that the RHS of the evolution equations in Eqs. (30)–(33) are equal to zero. For the potential energy of the inverse hyperbolic sine type, the critical points and the main physical properties are expressed in Table 1. It should be noted that for dust matter $w_{\rm m}=0$, only the critical point $A$ exists in the phase space. As can be noted from this table, one can observe two main classes of critical points. The first class, denoted as $A$ represents a critical line which corresponds to a cosmological constant behavior, having an interrelation between the auxiliary variable $u$ associated to the coupling function $g(\phi)$ and the strength of the potential energy, embedded into the non–constant $\lambda$ variable. This epoch corresponds to a de Sitter universe, a critical line where the dark energy field dominates in terms of density parameters. Analyzing the corresponding eigenvalues, due to the presence of a zero eigenvalue, the linear stability method fails to provide a viable theoretical framework for determining the stability properties. Hence, for this critical line we can only argue on the specific cases where we have a saddle cosmological behavior, due to the presence of eigenvalues with both positive and negative real parts. Concerning this critical line, we display in Fig. (2) various regions for the model’s parameters which correspond to a saddle cosmological behavior, due to the existence of at least one eigenvalue with negative real part, and at least an eigenvalue with positive real part.

A second class of critical points displayed in Table 1 is represented by the $B_{i}^{j}$, $i,j=\left\{+,-\right\}$ critical points which have a scaling cosmological behavior, an epoch in which the dark energy field acts as a matter component and mimics a matter dominated epoch. This type of solutions can in principle solve the cosmic coincidence problem. For these solutions, we have shown in Table 1 the locations in the phase space and the corresponding physical properties associated. From a physical point of view, the validity of these solutions implies that various existence conditions are satisfied. These corresponds to the requirement that the critical points belong to the real space, and the corresponding density parameters are physically viable, $\Omega_{\phi}=1-\Omega_{m}\in\left[0,1\right]$. Moreover, due to the definition of the dimensionless variables within this section, one should add the requirement that the $y$ variable is real and positive and $x\in\left(-1,1\right)$ due to the form of the matter (dark energy) density parameter. From the expression of the $x$ variable presented in the table, it can be seen that the existence conditions imply that for the present tachyonic dark energy model $w_{\rm m}\in\left[-1,0\right)$. Assuming that the matter component is embedded into the dark matter fluid, this implies that the pressure associated to the dark matter fluid is negative, an exotic situation which is not excluded by different cosmological observations Thomas:2016iav ; Yang:2016dhx ; PhysRevD.88.127301 ; PhysRevD.71.047302 . Moreover, for the critical points $B_{\left[+,-\right]}^{-}$, the existence conditions imply that the product $\omega_{1}\omega_{2}$ is negative, while for the solutions $B_{\left[+,-\right]}^{+}$ we have an inverted situation, $\omega_{1}\omega_{2}>0$. Hence, the current tachyonic dark energy model might contribute to a solution of the cosmic coincidence problem due to the scaling solutions since it can recover matter and de Sitter cosmological epochs. The stability of the points $B_{\pm}^{\pm}$ depend on the sign of the eigenvalues. In general, the stability conditions for each point are very cumbersome since they depend on the parameters $\alpha,w_{\rm m}$ and $\omega_{1,2}$. Concerning the stability properties, we show in Fig. 3 various cases which corresponds to the stability associated to the $B_{-}^{-}$ and $B_{+}^{+}$ critical points, determining possible values of the $\omega_{1,2}$ parameters which result in stable scaling solutions. Point $A$ is a non-hyperbolic point and standard linear stability theory fails on describing any stability property of it. One can use other dynamical system techniques as centre manifold theory to study its stability (see BAHAMONDE20181 ; Leon:2015via ). However, due to limited physical effects associated to the $A$ critical point, we relied our analysis only on linear stability methods, exploring the specific conditions where the stability corresponds to a saddle dynamical behavior. The present discussion can be adapted also for various potential types beyond exponential which have been studied in Ref. Roy:2017uvr with compatible results.

IV Conclusions

In this paper, we have considered a new dark energy model in the teleparallel equivalent of general relativity, based on modifications due to tachyonic fields which are non–minimally coupled with the torsion scalar and its boundary term. In this approach, the boundary term is related to the divergence of the torsion vector and the torsion scalar reproduces the same theory as GR at the level of the field equations. After finding the corresponding field equations for this tachyonic dark energy model, we analyzed the effects of the non–minimal coupling by employing the linear stability theory. In the first scenario, we considered the case where the coupling function is represented by a power law dependence, and the potential energy term corresponds to an exponential. In this case, we observed that the structure of the phase space has a reduced complexity. The critical points which are present corresponds to a cosmological constant behavior. Hence, in this case one notices that the evolution of the dynamical system can explain the current accelerated expansion of the present Universe due to the cosmological constant behavior of the system at the critical points. However, due to the reduced complexity of the phase space and without the presence of the scaling solutions, the cosmic coincidence problem cannot be alleviated.

A second cosmological scenario was also considered by taking into account that the non–minimal coupling functional has an exponential type parameterization. Concerning the potential energy term in the action, the study considered that the potential is beyond the usual exponential type found in many dynamical constructions in scalar tensor theories. As can be noted from the previous section, the dynamical equation associated to the potential term involves a second order polynomial parameterization. In this specific case in the analysis, the potential corresponds to an inverse hyperbolic sine. For the second cosmological scenario the phase space have a 4-dimensional structure and a richer complexity. As can be noted from the previous section, we have shown that the second cosmological scenario analyzed can reproduce the known evolution of the Universe and solve the cosmic coincidence problem due to the existence of scaling solutions. In these critical points, the dark energy field mimics a matter era due to the specific form the corresponding effective equation of state. However, the existence conditions associated to the scaling solutions imply the existence of an exotic warm dark matter fluid having a limited negative pressure, a comportment not ruled out by present astrophysical observations Thomas:2016iav ; Aghanim:2018eyx . Furthermore, the rest of the critical points in the phase space corresponds to a cosmological constant behavior and can explain the current accelerated expansion of the Universe with a constant equation of state.

One can then say that the current dark energy model constructed in the teleparallel equivalent of general relativity modified by a tachyonic field non–minimally coupled with a boundary term represents a potentially realistic model in scalar tensor theory which might solve the cosmic coincidence problem and the nature of the dark energy phenomenon, a feasible tachyonic prototype. This offers one possible alternative avenue to tachyonic fields to compliment the curvature-based work in the literature Quiros:2009mz ; Aguirregabiria:2004xd ; Bagla:2002yn ; Fang:2010zze ; Otalora:2013tba ; Otalora:2014aoa .

Acknowledgements.

The JLS and SB would like to acknowledge networking support by the COST Action GWverse CA16104. This article is based upon work from CANTATA COST (European Cooperation in Science and Technology) action CA15117, EU Framework Programme Horizon 2020. SB is supported by Mobilitas Pluss N∘ MOBJD423 by the Estonian government. M. Marciu acknowledge partial support by the project 29/2016 ELI-RO from the Institute of Atomic Physics, Bucharest–Magurele.