![Arxiv:1110.1736V3 [Gr-Qc] 9 Aug 2012](http://data.docslib.org/img/94f76551f10ca168d69849108baa3cfa-1.webp)
arXiv:1110.1736v3 [gr-qc] 9 Aug 2012 ysvrlatos n[,9 ti hw htte(in- the that shown is conditions it wide allow 9] model [8, potential In power-law studied verse) been authors. have Universe several stage the dominant by of dynamics quadratic The the energy. during Fried- modifica- high the gravitational very at in for tions density responsible energy is total equation a mann the of 4- of appearance standard the term fact the quadratic In from [5–7]. different cosmology on essentially dimensional equations are field cosmological brane the the that well-known constant. is cosmological a It the negative in with embedded with bulk dimen- tension start four 5-dimensional positive a model with in hypersurface the confined sional model of standard inflationary setup the of the The particles because in [2–4]. reasons, impact other scenario cosmological among appreciable years, last its the in intensively studied been the have to [1], mechanism compactifications alternative Kaluza-Klein an as originally introduced § ‡ † ∗ Rand a in trapped fields quintessence of analysis space Phase lcrncades yoelsy.leyva@fisica.ugto.mx address: Electronic [email protected] address: Electronic [email protected] address: Electronic lcrncades [email protected] address: Electronic h adl-udu rn ftp Imdl(RS2), model II type of brane Randall-Sundrum The iiio eCeca neirı el nvria eGua de Universidad la Ingenier´ıa de e Ciencias Divisi´on de eatmnod ae´tc,UiesddCnrld a Vi Las de Central Matem´atica, Universidad de Departamento nttt eFıia otfii nvria eCat´olica d de Universidad F´ısica, Pontificia de Instituto ASnmes 42.q 42.v 42.b 45.d 11. 04.50.Kd, 04.20.Jb, 98.80.Jk 04.20.Cv, 98.80.Cq, 04.20.-q, numbers: PACS R gravity, modified system, dynamical Cosmology, KeyWords: type). saddle (of unstable is solution Sitter rnin rmriliflto.I h atclrcs fas a of sad case are particular domin (they the scalar-field In 5D-modifications the the inflation. with on primordial of sufficient attractors transient conditions stability time obtain obtain the late to We for not employed well solution. as is Sitter solutions Theory de Manifold of stability Center Randall-Sundrum Friedmann-Robertson-Walkerfluid. a isotropic in and trapped homogeneous an potential arbitrary with field epoeta for that prove we nti ae eivsiae rmtednmclssesper systems dynamical the from investigate, we paper this In eatmnod ´sc,UiesddCnrld a Villas Las de Central F´ısica, Universidad de Departamento .INTRODUCTION I. < χ eatmnod ´sc nvria eCmguy Cuba. Camag¨uey, de F´ısica Universidad de Departamento h eSte ouini smttclysal.However, stable. asymptotically is solution Sitter de the 0 Dtd coe 5 2018) 25, October (Dated: alsR Fadragas R. Carlos aoet Escobar Dagoberto endstudy refined a olyLeyva Yoelsy el Leon Genly xaso steMdfidCalgnGs[0.Frthis For [20]. Gas accelerated Chaplygin the our Modified explaining of the for feature is these approach expansion have with One deal framework, to RS2 universe. order on in based proposed been models, expansion Several accelerated an phase. experiences currently our that universe Cosmic confirm and strongly [16] [17–19] Background Structure Microwave Scale Large [13–15], pernovae nrydniyo oemte opnn rwa the as the grow if component modified [12] matter be proceeds some can that expansion of expansion is density cosmic scenario the energy this of of fate feature the recently and interesting rsneo efc udwt rirr aorpcin- barotropic arbitrary with fluid the perfect dex in a model scaling potential of power-law the presence of of case the stability for The solutions with not scenarios. contrast do favorable which in models present k-essence scenario and quintessence potential exponential successful a for y[1 oagnrlzdbackground generalized a to [11] by trary 1 aut,AP 5,310 eo,Gaaut,M´exico. Le´on, Guanajuato, 37150, 150, A.P. najuato, socr ihtecmo atrsucs unesnescal quintessence sources: etc. matter dust, common radiation, o expan density the field, cosmic energy with the the with occurs dilutes if as brane appear the in not content do matter the modification of kind This eea srpyia bevtossc sTp aSu- Ia Type as such observations astrophysical Several apr´s,Csla45,Valpara´ıso, Chile. 4950, Valpara´ıso, Casilla e ‡ 5-,1.5W,9.6+,9.0-,98.80.Bp, 98.80.-k, 95.36.+x, 11.25.Wx, 25.-w, § γ aa edwt potential with field calar n a eeoe y[0.Ti td a extended was study This [10]. by developed was tdslto.W rv h hr are there the prove We solution. ated ∗ † ls at lr P580 ua and Cuba. 54830, CP Clara Santa llas, teR2cs srcvr when recovery is case RS2 (the FW rn le lowt perfect a with also filled brane (FRW) sBaeol ftp I econsider We II. type of Braneworld ’s l-ie.Ti atcreae iha with correlates fact This dle-like). oeta o h tblt fscaling of stability the for potential andall-Sundrum at lr P580 Cuba. 54830, CP Clara Santa , pcie h vlto fa scalar an of evolution the spective, odtosfrteasymptotic the for conditions l-udu Braneworld: all-Sundrum 1 for V = > χ V 0 e H h de the 0 − χφ 2 ∝ n Λ + ρ ) Another 2). = T n o narbi- an for sion ar f 2 models was showed that the Universe follows a power law- cal discussion of the above results, while the conclusions expansion around the critical points. Another important are given in Section IV. approach consist in adding a self-interacting scalar field to the matter inventory in the brane [8, 21–24]. Scalar fields arise in a natural way in particle physics and they can act as a candidate for dark energy playing the roles II. DYNAMICAL SYSTEMS ANALYSIS OF THE such that quintessence, phantom, quintom, tachyons field FRW BRANE etc. [25] The dynamical behavior of scalar field coupled with a In this section we will focus our attention in a brane- barotropic fluid in a spatially flat Friedmann-Robertson world model where an scalar field, with arbitrary self- Walker universe has been studied by many authors, see interaction potential, is trapped on a RS2 brane. In the for instance the references [26–32]. A natural general- flat FRW metric, the field equations read [36–39]: ization of [26], is to include higher-dimensional behav- ior (RS2 scenario). This program was carried out in [33] where was found the corresponding scalar field potentials 1 ρT 2 H2 = ρ 1+ + U (1) which lead to attractor scaling solution for several energy 3 T 2λ λ density modifications to the Friedmann equation; for the ρT 2 4 2 2H˙ = (1 + )(φ˙ + γρm) U (2) RS2 framework the potential V cosech (Aφ) was found − λ − λ 2 ∝ . The dynamics of a scalar field with constant and expo- ρ˙m = 3γHρm (3) nential potentials was investigated in [34]. These results − φ¨ + ∂φV = 3H φ˙ (4) were extended to a wider class of self-interaction poten- − tial in [35] using a method proposed by [29] supporting the idea that this scenario modifies gravity only at very where we have used the Randall-Sundrum fine tunning high energy/short scales (UV modifications only) having condition i.e., we neglect the cosmological constant term an appreciable impact on primordial inflation but does (Λ4 = 0). ρT = ρφ + ρm, λ is the brane tension, γ not affecting the late-time dynamics of the Universe 3. is the barotropic index of the background fluid, V is C In this paper we make a step forward with respect to the the scalar field self-interaction potential. (t) = 4 U a(t) previous studies by exploring more deeply the dynam- is the dark radiation term which arises from a non- ics in the phase space associated to this scenario around vanishing bulk Weyl tensor being C a constant parameter both hyperbolic an non-hyperbolic critical points. The related with black hole mass in the bulk: if the bulk is last subject cannot be consistently studied with the help AdS Schwarzschild C = 0 [37]. When the black hole of linear analysis, but using the Center Manifold Theory. mass− vanishes, the bulk6 geometry reduces to AdS, and Our claim is that the more interesting solution are the C = 0 [7, 39]. In the following we will study the latter non-hyperbolic critical points, in particular the de Sitter case i. e., we will not consider here the dark radiation critical points. term . Here and throughout, we use ∂φ to denote deriva- In this paper we employ the Center Manifold Theory tive with respect to φ. to obtain sufficient conditions for the asymptotic stabil- From the Friedmann equation (1) it is deduced how ity of de Sitter solution and for proving that here are the brane effects modify the early time dynamics: at high not late time attractors with 5D-modifications. This fact energy (ρT >> λ) this equation reduces to H ρT . A correlates with a transient primordial inflation. Also we late times, due to the expansion rate, the energy∝ density provide conditions on the potential for the stability of of the matter trapped in the brane dilutes (ρT << λ), scaling solutions as well for the stability of the scalar- and the standard 4D TGR behavior is recovered, leading field dominated solution. to H √ρT . ∝ The paper is organized as follows. In Section II we give Having presented the cosmological equations, our pur- the essential details of the Randall-Sundrum Model and pose now is to define a dynamical system from (1)-(4) in deals with the dynamical system analysis, these include order to examine all possible cosmological behaviors. As the center manifold study. We explore the dynamics of we know dynamical systems techniques provides one of an scalar field with exponential potential plus a Cosmo- the better tools for obtaining useful information about logical Constant trapped on the brane using the previous the evolution of a wide class of cosmological models 4. In results in Section III. Section IV is devoted to the physi- order to take advantage from these tools, we introduce the Hubble-normalized variables
˙ 2 φ V ρT 2 x = y = Ωλ = , (5) At early times, where the quadratic energy term dominates, this √6H 3H2 6λH2 potential behave as a (inverse) power-law potential being consis- tent with previous analysis [8, 9]. 3 A word of caution: this claim is not in general true, specially if the energy density of the matter trapped in the brane increase at late times [12] 4 See, for instance, the seminal work [26], and [40, 41]. 3 the new temporal variable τ = Hdt, and the additional 5. Due to its classic nature, our model is not appropri- dynamical (non-compact) variable, s, given by ate to describing the dynamics near the initial singular- R ity, where quantum effects appear. However, from the s = ∂φ ln V (φ). (6) mathematical viewpoint, this region (Ωλ = 1) is reached − asymptotically. In fact, as some numerical integrations which is a function of the scalar field. corroborate, there exists an open set of orbits in the phase For the scalar potential treatment, we proceed follow- interior that tends to the boundary Ωλ =1 as τ . ing the reference [29]. Let be defined the scalar function Therefore, for mathematical motivations it is common→ −∞ to attach the boundary Ω = 1 to the phase space. On ′′ λ V V the other hand the points with (Ωλ = 0) are associated f =Γ 1, Γ= ′ . (7) − V 2 to the standard 4D behavior (ρT λ or λ ) and corresponds to the low energy regime.≪ → ∞ Since Γ is a function of the scalar field Γ(φ) (see definition From definition (5) and from the restriction (13), and (7)), also is the variable s = S(φ). Assuming that the taking into account the previous statements, it is enough inverse of S exists, we have φ = S−1(s). Thus, one can to investigate to the flow of (8)-(11) defined in the phase obtain de relation Γ = Γ(S−1(s)) and finally the scalar space field potential can be parameterized by a function f(s). 2 Thus, in general, it is possible the treatment of general Ψ= (x, y, Ωλ):0 x + y +Ωλ 1, 1 x 1, f { ≤ ≤ − ≤ ≤ classes of potentials using an “ -deviser”. In the table I 0 y 1, 0 Ωλ 1 s R . (15) are shown the functions f(s) for some usual quintessence ≤ ≤ ≤ ≤ }×{ ∈ } potentials. The cases (a)-(c) have been studied in RS2 Some cosmological parameters like the equation of pφ branes in [35]. state parameter of the scalar matter ωφ = , the decel- ρφ Bearing this in mind, and using the variables (5)-(6) H˙ eration parameter q = 1+ 2 can be re-expressed we deduce the following autonomous system of ordinary − H differential equations (ODE) as functions of the new variables as follows 2 x y 2 ′ 3 3(Ωλ + 1)(γ 2) 3 ω = − , Ω = x + y (16) x = sy 3x + − x + φ x2 + y φ 2 − 2(Ωλ 1) r − 3γ(Ωλ + 1)(y +Ωλ 1) + − x (8) 2(Ωλ 1) 1+Ωλ 2 3γ 2 − q = 3x + 1 x y Ωλ 1 (17) ′ 3y(Ωλ + 1) 2 1 Ωλ 2 − − − − y = (γ 2) x + γ (Ωλ + y 1) √6xys − (Ωλ 1) − − − − (9) A. Critical points ′ 2 Ω = 3Ωλ (γ 2) x + γ (Ωλ + y 1) (10) λ − − ′ 2 s = √6xs f(s), (11) The system (8)-(11), admits the curves of critical − ± points P1, P2, P3; the critical points P4 and P5; and the ± where the comma denotes derivative with respect τ. classes of critical points P6 , P7 and P8 parameterized From the Friedmann equation (1) follows the relation by s∗ satisfying f(s∗)=0. In Table II are displayed the location, existence conditions and some basic observables 2 6 ± Ωm =1 x y Ωλ (12) of these critical points . The critical points of P1 to P6 − − − ∗2 always exist; the point P7 exists for s 3γ, whereas, ∗2 ∗ ≥ Using (12), the energy condition 0 Ωm 1 can be P exists for s 6 with f(s )=0. ≤ ≤ 8 written as Now let us comment≤ on the stability of the first or- der perturbations of (8)-(11) near the critical points dis- 2 0 x + y +Ωλ 1. (13) played in table II. Let us comment briefly in their physical ≤ ≤ interpretation. From the definition of Ωλ and the Friedmann equation The line of critical points y =1 Ωλ called P represent − 1 we obtain the useful relation solutions with 5D-corrections, since, in general, Ωλ = 6 ρ 2Ω T = λ (14) λ 1 Ωλ − 5 See the references [30, 50] for a classical treatment of cosmological From (16) follows that the invariant set Ωλ = 1 cor- solutions near the initial singularity. responds to cosmological solutions where ρT λ (cor- 6 Strictly speaking the system admits one more curve of critical 2 ≫ x (γ 2) responding to the formal limit λ 0). Therefore, they point with coordinates x ∈ [−1, 1] ,y = − − , Ω = 1,s = 0, → γ λ are associate to high energy regions, i.e., to cosmologi- but since the energy condition (13) is not satisfied, we omit it cal solutions in a neighborhood of the initial singularity from the analysis. 4
TABLE I: Explicit forms of f(s) for some self-interaction potentials. To homogenize the notations we use units in which κ2 8πG = 1. ≡
Label Potential f(s) Reference
2 α 1 αχ (a) V = V0 sinh− χφ 2 [42, 43] α − s 2 α 1 αχ (b) V = V0[cosh(χφ) 1] + 2 [44] − − 2α 2s V0 1 α (c) V = β + [45] (η+e−αφ) β s χφ χ (d) V = V0e− + Λ 1 [46] − − s 2 2 2 eχφ s +8mχ+s√s +8mχ (e) V = V0 φm 2ms2 [47, 48] αφ βφ (s+α)(s+β) (f) V = V0 e + e 2 [49] − s
TABLE II: Location, existence conditions and some basic observables for the critical points of the system of equations (8)-(11).
Pi x y Ωλ s Existence ωφ Ωφ q
P1 0 1 Ωλ Ωλ [ 0, 1 [ 0 Always 1 1 Ωλ 1 − ∈ − − − 3γ P2 0 0 0 s R ” undefined 0 1 ∈ 2 −
P3 0 0 1 s R ” undefined 0 undefined ∈ P ± 10 00” 112 4 ±
P5 01 0 0” 1 1 1 − − P ± 1 0 0 s∗ ” 1 1 2 6 ± 3 √ 2 γ 3(γ 2)γ 2 3γ 3γ P7 ∗ −∗ 2 0 s∗ s∗ 3γ γ 1 ∗2 1 s − 2(s ) ≥ − s 2 − 2 ∗ ∗ ∗2 s (s ) 2 1 2 s P8 1 0 s∗ s∗ 6 s∗ 3 1 1 √6 − 6 ≤ 3 − 2 −
0. From the relationship between y and Ωλ follows that proving that this critical point is locally asymptotically this solution is dominated by the potential energy of the unstable. ρ V φ scalar field T = ( ); that is, it is de Sitter-like solution The point P2 represents a matter-dominated solution (ωφ = 1). In this case the Friedmann equation can be − (Ωm = 1). Although it is non-hyperbolic, it behaves like expressed as a saddle point in the space of phase of the RS model, since they have a nonempty stable and unstable manifolds (see V 3H2 = V 1+ (18) the table III) 7. 2λ The critical point P3 is located at the boundary Ωλ =1 In the early universe, where λ V, the expansion rate of the phase space region (15). From the physical view- of the universe for the RS model≪ differs from the general point, this solution represents the Big Bang singularity relativity predictions (ρT ). The eigenvalues for P3 are displayed in tables III. They→ ∞ were calculated for orbits contained completely H V in the invariant set x = y = 0 and by taking the limit RS = (19) HGR r2λ
P1 admits a 2D stable manifold, M2. Due the importance of de Sitter solutions in the cosmological context, in sec- 7 Strictly speaking, the concept of saddle point is not applicable tion II A 2 we calculate explicitly their center manifold to nonhyperbolic critical points. 5 as Ωλ 1. For orbits outside the above invariant set we The sufficient conditions for the asymptotic stability of → cannot make the above limit process since the system is the scalar-field-dominated solution (P8) are either not of class C1 at Ω = 1. However, several numerical λ iii) 0 γ < 2, √3γ
TABLE III: Eigenvalues for the critical points of the equations system (8)-(11). We use the notation β = ± 3 24γ2 γ γ ∗ 2 γ . 4 2 (2 ) (s ) 9 + 2 − ± r − −
Pi λ1 λ2 λ3 λ4
P1 0 0 3 3γ − − 3 P2 0 3γ 3γ (γ 2) − 2 −
P3 0 3(γ 1) 3γ 6γ − P ± 0 6 6 6 3γ 4 − −
P5 0 0 3 3γ − − 2 P ± 6 6 3γ 6 √6s∗ √6(s∗) f ′(s∗) 6 − − ∓ ∓
P7 3γ β β+ 3γs∗f ′(s∗) − − − 1 2 2 2 3 P8 s∗ 6 s∗ 3γ s∗ s∗ f ′(s∗) 2 − − − −
∂G ∂G y′ x′ x′ . Neglecting the error terms, and introducing the co- 2 1 2 =0 (29) 2 − ∂x1 − ∂x2 ordinate transformation u1 = x1, the system (34)-(35) reduces to the simpler form Since we have used Taylor expansions up to third order for obtaining the system (23)-(26) we must seek a solu- tion for (28)-(29) in the following form (see [32, 40, 41]): u′ = 2u2f(0) (36) 1 − 1 x′ = u x , (37) 3 2 2 2 2 − 1 2 F (x1, x2)= a1x1 + a2x1 + a3x1 + a4x1x2 + a5x1x2+ 3 2 + a6x2 + a7x2 + (4) (30) where the region of physical interest is u 0, x 0. O 1 ≥ 2 ≥ Observe that the dynamics on the center manifold, 3 2 2 2 governed by (36)-(37), depends on the value f(0). If ei- G(x1, x2)= b1x1 + b2x1 + b3x1 + b4x1x2 + b5x1x2+ ther f(0) = 0 or f is singular at the origin, the system + b x3 + b x2 + (4), (31) 6 2 7 2 O (36)-(37) does not represents correctly the dynamics of the center manifold. I such a case, we must incorporate as xi 0 where (4) is an error term of fourth order → O higher order terms in the scheme, increasing the problem in the vector norm. Substituting expressions (30) and complexity. Thus, we assume that f(0) is a real number (31) in the equations (28)-(29), and comparing terms of such that f(0) =0. equal powers, we obtain that the non-null coefficients in 6 the above expressions (30) and (31) are According to the Center Manifold Theorem, the sta- bility analysis of P5 is reduced to the analysis of the f(0) 1 1 1 stability of the origin of the system (36)-(37). For this a1 = , a5 = , b2 = , b3 = , (32) analysis we resort to numerical investigation. In figure 3√6 −√6 −6 3 1 are displayed several orbits contained the physical re- i.e., gion u1 0, x2 0. Observe that the axis are invariant sets. For≥f(0) >≥0 (see the panel (a) in figure 1), there 3 x1x2 x1f(0) is an open sets of orbits that converge to the origin as y1 = F (x1, x2)= + + (4), − √6 3√6 O time goes forward; thus, the origin is asymptotic stable for initial conditions in a vicinity of the origin whenever x2 x2x y = G(x , x )= 1 + 1 2 + (4) (33) f(0) > 0. From the asymptotic stability of the origin of 2 1 2 − 6 3 O (34)-(35) follows that, for f(0), the center manifold of P5 Thus, the dynamics on the center manifold, is given by is locally asymptotic stable, and hence, the solution P5 of the system(8)-(11) also is. Therefore, P5 with f(0) > 0 x′ = x3f(0) + (4) (34) corresponds to a late time de Sitter attractor. This result 1 − 1 O for RS2 brane cosmology is in a perfect agreement with x′ = x2x + (4). (35) 2 − 1 2 O the standard 4-dimensional TGR framework. 7
tion is of saddle type, we can correlate this behavior with a transient inflationary stage for the universe. In order to verify our claim, we appeal to the Center Manifold Theory. Let us consider an arbitrary critical point with coor- dinates (x = 0,y = 1 uc, Ωλ = uc,s = 0) located at − P1. In order to prepare the system (8)-(11) for the appli- cation of the Center Manifold Theorem we introduce the coordinate change
s(uc 1) u1 = − , u2 = Ωλ uc(y +Ωλ 2), − √6 − − − s(uc 1) v1 = (uc + 1)(y +Ωλ 1), v2 = − + x. (38) − √6 ′ ′ ′ ′ Then, we Taylor expand the system u1,u2, v1, v2 in a neighborhood of the origin with error of order (4). Accordingly to the Center Manifold theorem,O the local center manifold of the origin for the resulting vector field (a) is given by the graph: W c (0)= (u ,u , v , v ): v = F (u ,u ), loc { 1 2 1 2 1 1 1 2 2 2 v2 = G1(u1,u2),u1 + u2 <δ (39) for δ > 0 a small enough real value. The functions F1 and G1 in definition (39) are a so- lution of a system of quasi-linear differential equations analogous to (28)-(29). This system should be solved with an error of order (4), obtaining the functional de- pendence: O 2 2u1u2(uc + 1) 2 v1 = u1(uc + 1), uc 1 − 3 − 3 2 u1uc u1f(0) u1u2 v2 = − . (40) uc 1 − uc 1 − − Then, the dynamics on the center manifold is given by 3 ′ 6u1f(0) u1 = + (4) (41) uc 1 O − 2 ′ 2 6u2(1 3uc)u1 u2 =6ucu1 + − + (4). (42) uc 1 O − In the same way as for P , the dynamics of the sys- (b) 5 tem (41)-(42) depends on the values of f(0). We assume that f(0) R 0 . Otherwise it is required to include FIG. 1: Phase space of the system (36)(35) for: (a) f(0) = 1 higher order∈ terms\{ } in the Taylor expansion, increasing and (b) f(0) = 0.1. the numerical complexity. − In figure (2) are displayed some orbits in the phase space of the system (41)-(42) for the choices: (a) f(0) = 2 and uc = 0.5 and (b) f(0) = 2 and uc = 0.5. The 2. Dynamics of the center manifold of P1 origin of coordinates is locally− asymptotically unstable (of saddle type) irrespective the sign of f(0). Henceforth, In this section we investigate the stability of the curve the center manifold of P1 is locally asymptotic unstable of critical points P1 for 0 < Ωλ < 1. This solutions cor- (saddle type) for f(0) = 0; also is P1. respond to a de Sitter expansion with 5D-corrections. The physical interpretation6 of this result is that there According to the RS2 model this solution cannot behave are not late time attractors with 5D-modifications. This like a late time attractor since 5D-corrections are typical type of corrections are characteristic of the early universe. of the high energy (early universe) and not for low energy In this sense the cosmological solution associated to the (universe late) regimes. If we can prove that this solu- critical P1 correlates with the primordial inflation. 8
This potential have been widely investigated in the lit- erature. It was studied for quintessence models in [46] were it is considered a negative cosmological constant Λ. In our case of interest we assume Λ 0, to avoid dealing with negative values of y. However≥ we can ap- ply our procedure by permitting negative values for y for the case Λ < 0. The dark energy models with ex- ponential potential and negative cosmological constant were baptized as Quinstant Cosmologies. They were in- vestigated in [51] using an alternative compactification scheme. The asymptotic properties of a cosmological model with a scalar field with exponential potential have been investigated in the context of the General Relativ- ity by the authors of [26, 29], and in the context of the RS braneworlds by [34, 52]. In both cases it was studied the pure exponential potential (Λ = 0). Potentials of ex- ponential orders at infinity were studied in the context of Scalar-tensor theories and conformal F (R) theories by the authors of [30, 32]. We comment that the procedure introduced in previ- (a) ous sections is fairly general and can be applied to others potentials as those showed in table I. We remain in the exponential potential for simplicity. Also, we consider a pressureless (dust) background, i.e., γ =1. The function f(s) corresponding to the potential (43) is given by χ f(s)= 1 . (44) − − s The zero of this function is χ 1 s∗ = χ f ′(s∗)= = . (45) − s∗2 χ Observe that for the potential (43) s∗f ′(s∗) < 0. Thus, the only relevant late-time attractor should the de Sitter solution. In fact, the critical points P7 of the table II are reduced to the single point
3 1 3 P , , , χ , 7 = 2 0 (46) −r2 χ −2χ !
(b) which represents a saddle point in the phase space. The critical points P8 are reduced to the single one
FIG. 2: Phase space of the system (41)-(42) for the choices χ χ2 (a) f(0) = 2 and uc = 0.5 (b) f(0) = 2 and uc = 0.5. Ob- P8 = , 1 , 0, χ . (47) − −√6 − 6 − serve that the axis u1 is a line of points of critical asymptotic unstable in both cases for initial conditions in a vicinity of This point represents a scalar-field dominated solution the origin. The origin behaves as a saddle point. (Ωφ = 1). It is a saddle point in the phase space. Observe that all the trajectories in the phase space always emerge from the point (x, y, Ωλ)=(0, 0, 1). In the Fig 3, we present, for different choices of the III. EXPONENTIAL POTENTIAL free parameters, two numerical integrations which sug- gest that P5 that is a de Sitter late-time attractor with The objective of this section is to illustrate our analyt- an standard 4D behavior (Ωλ = 0). However, in order ical results for the exponential potential, the prove this claim we need to use the Center Manifold Theory. Although for the potential (43) the result in the V (φ)= V e−χφ +Λ. (43) appendix does not apply since f(0) = sgn(χ) is not a 0 − ∞ 9
y 1.0 defined in the phase plane x R, x [0, 1]. 0.5 1 ∈ 2 ∈ 0.0 In general, the system (8)-(11) is not invariant under 0.0 the change (s, x) ( s, x) unless f(s) is an even func- tion, i.e., f( s) =→f(−s). −However, in the particular case of the potential− (43) we observe from (44) that the sys- tem (8)-(11) is invariant under the discrete symmetry W0.5 (s,x,χ) ( s, x, χ). It is easy to show that for the exponential→ − case− the− symmetry (s,x,χ) ( s, x, χ) is induced by the discrete symmetry (φ,→ χ) − ( −φ, −χ). 1.0 → − − 0.4 On the other hand, for the potential (43) the function 0.2 0.0 s(φ) is given by x -0.2 -0.4 V χ s(φ)= 0 . V + eφχΛ (a) 0
Assuming V0 > 0 and Λ > 0 we have sgn(s) = sgn(χ). s φ y 1.0 As we see, the function ( ) is not defined for all the real 0.5 values, that is, the region of physical interest is restricted 0.0 to one semiplane s ≶ 0 (or in the degenerate case χ =0 0.0 to the hyperplane s = 0) depending of the sign of χ. That is, s is either positive, negative or zero, for χ either positive, negative or zero, provided V0 > 0 and Λ > 0. This fact limits the application of the above discrete W0.5 symmetry. In the figures (4) (a) and (4) (b) are displayed some orbits of the flow of (49)-(50) for the choice χ = +0.5, and χ = 0.5, respectively. In both cases, the x - − 2 1.0 axis is invariant so the orbits cannot cross throught it, 0.4 s 0.2 i.e., the orbits with (0) ≶ 0 at the intial time, remains in the respective semiplane all the time. Also, in both x 0.0 figures the flow is the same, thus, the discrete symmetry (a) (s,x,χ) ( s, x, χ) is preserved. However, from the physical→ view− point,− − for the choice χ = +0.5 the physical region is the right semiplane (x1 s > 0). Thus, the physical orbits depart form the origin≡ as the time goes
FIG. 3: Some orbits in the projection (x,y, Ωλ) of the phase forward. This means that the origin is unstable. For the χφ space for (8)-(11)) and potential V (φ) = V0e− + Λ for the choice χ = 0.5, the physical region is the left semiplane − choice (a) (γ,χ)=(1, 0.5), and (b) (γ,χ)=(1, 100). x1 s< 0, thus, the physical orbits approach the origin as the≡ time goes forward. Summarizing, for V0 > 0 and Λ > 0 we have sgn(s)= sgn(χ). Thus, for the potentials (43) with χ < 0 the de real number, we can use the same procedure as in section Sitter solution (P5) is asymptotically stable. This fact is III to obtain the center manifold of P5 by setting from self-consistent with the condition f(0) = sgn(χ) 0; the beginning the functional form of f(s) given by (44). the sufficient condition for the stability− of the de∞ Sitter ≫ For P5, the graph of the center manifold is given, up solution derived in section II A 1. However, for the po- to an error term (4), by O tentials (43) with χ> 0 the de Sitter solution is unstable (of saddle type) to perturbations in the s-direction. In 2 3 χ 1 x1 1 2 χ 2 x2x1 this case the late-time solution is given by the asymptotic y1 = − + x2χ x1 , 3√6 3 3 − 3√6 − √6 configuration x =0,y =1, Ωλ =0, φ . r ! → −∞ χx3 x 1 y = 1 + 2 x2 (48) 2 9 3 − 6 1 IV. RESULTS AND DISCUSSION where we have introduced the variables (22). The dynamics on the center manifold is governed by The main results of this paper can be summarized as the equations follows. The critical point P3 = (0, 0, 1) represents a Big Bang singularity. According to our numerical integra- χ2 tions in Figs. 3, we observe that all the trajectories in x′ = 1 x3 + χ(1 x )x2 + (4), (49) 1 − 3 1 − 2 1 O the phase space, but a measure zero set, emerge from the vicinity of this point. This result agrees with the previous x′ = x2x + (4), (50) 2 − 1 2 O results obtained in [35]. In [35] the authors use another 10
we have the following results. By an explicit computation of the center manifold of P1 and of P5 we prove that
P1 is locally asymptotic unstable (of saddle type) • irrespectively the sign of f(0) R 0 . This fea- ture is corroborated in the Figures∈ \{ 2. }
P5 is locally asymptotically stable for f(0) > 0 and • unstable (of saddle type) for f(0) < 0. This result is illustrated in the Figures 1 and 4.
The solutions dominated by the kinetic energy of the ± ± scalar field P4 and P6 behave like saddle-type solutions. This is a main difference with respect the standard 4D theory where this type of solutions are always past at- tractors. In this general case, the possible late-time attractors (a) are:
the standard 4D de Sitter solution P5 (ωφ = 1) • whenever f(0) > 0; −
the matter-scalar-field scaling solution P (Ωφ • 7 ∼ Ωm). The sufficient conditions for its asymptotic stability are s∗ < √3γ,f ′ (s∗) < 0 or s∗ > √3γ,f ′ (s∗) > 0; and−
the scalar-field-dominated solution P8 (Ωφ = 1). • The sufficient conditions for its asymptotic stabil- ity are √3γ
V. CONCLUSIONS
(b) In the present paper we have investigated the phase space of the Randall-Sundrum braneworlds models with a self-interacting scalar field trapped in the brane with FIG. 4: Some orbits of the flow of (49)-(50) for the choices: arbitrary potential. (a) χ = 0.5, and (b) χ = 0.5. In both cases the x2-axis is − From our numerical experiments we claim that P3 is invariant. In (a) the physical region is the semiplane x1 > 0; associated with the Big Bang singularity type. The nu- s > for orbits with (0) 0 the origin is unstable to perturbations merical investigation suggest that it is always the past in the s-direction. In (b) the physical region is the semiplane attractor in the phase space of the Randall-Sundrum cos- x1 < 0. Thus for orbits with s(0) < 0 the origin (hence P5) is asymptotically stable. mological models. Using the center manifold theory we have obtained suf- ficient conditions for the asymptotic stability of de Sitter solution. coordinate system, that is equivalent, except diffeomor- We have obtained conditions on the potential for the phisms, to the system used in this paper. stability of the scaling solutions as well for the stability In the particular case of a scalar field with poten- of the scalar-field dominated solution. −χφ tial V = V0e + Λ trapped in the brane, we have We have proved, using the center manifold theory proved that for χ < 0 the de Sitter solution (P5) is and numerical investigation, that there are not late time asymptotically stable. However, for χ > 0 the ori- attractors with 5D-modifications since they are always gin, i.e., P5 is unstable (of saddle type) and the late- saddle-like. This fact correlates with a transient primor- time solution is given by the asymptotic configuration dial inflation. x = 0,y = 1, Ωλ = 0, φ . This class of potentials In the particular case of a scalar field with potential → −∞ −χφ contains the previously studied potentials in [34] with V = V0e + Λ we have proved that for χ < 0 the Λ = 0. Thus our present results generalize those in [34]. de Sitter solution is asymptotically stable. However, for In the general case, for potentials satisfying f(0) R, χ> 0 the de Sitter solution is unstable (of saddle type). ∈ 11
Acknowledgments for partial financial support of this investigation. YL is grateful to the Departamento de F´ısica and the CA de This work was partially supported by PROMEP, Gravitaci´on y F´ısica Matem´atica for their kind hospital- DAIP, and by CONACyT, M´exico, under grant 167335; ity and their joint support for a postdoctoral fellowship. by MECESUP FSM0806, from Ministerio de Educaci´on The authors wish to thank two anonymous referees for de Chile; and by the National Basic Science Program their useful comments and for bringing our attention to (PNCB) and Territorial CITMA Project (No. 1115), several references. Cuba. DE, CRF and GL wish to thank the MES of Cuba
[1] Lisa Randall and Raman Sundrum. An Alternative to [17] N. Jarosik, C.L. Bennett, J. Dunkley, B. Gold, compactification. Phys.Rev.Lett., 83:4690–4693, 1999. M.R. Greason, et al. Seven-Year Wilkinson Mi- [2] Rachael M. Hawkins and James E. Lidsey. Inflation on crowave Anisotropy Probe (WMAP) Observations: Sky a single brane: Exact solutions. Phys.Rev., D63:041301, Maps, Systematic Errors, and Basic Results. Astro- 2001. phys.J.Suppl., 192:14, 2011. [3] Greg Huey and James E. Lidsey. Inflation, brane worlds [18] D. Larson, J. Dunkley, G. Hinshaw, E. Komatsu, and quintessence. Phys.Lett., B514:217–225, 2001. M.R. Nolta, et al. Seven-Year Wilkinson Microwave [4] Greg Huey and James E. Lidsey. Inflation and brane Anisotropy Probe (WMAP) Observations: Power Spec- worlds: Degeneracies and consistencies. Phys.Rev., tra and WMAP-Derived Parameters. Astrophys.J.Suppl., D66:043514, 2002. 192:16, 2011. [5] Pierre Binetruy, Cedric Deffayet, and David Lan- [19] E. Komatsu et al. Seven-Year Wilkinson Microwave glois. Nonconventional cosmology from a brane universe. Anisotropy Probe (WMAP) Observations: Cosmological Nucl.Phys., B565:269–287, 2000. Interpretation. Astrophys.J.Suppl., 192:18, 2011. [6] Pierre Binetruy, Cedric Deffayet, Ulrich Ellwanger, and [20] Prabir Rudra, Ritabrata Biswas, and Ujjal Debnath. David Langlois. Brane cosmological evolution in a bulk Dynamics of Modified Chaplygin Gas in Brane World with cosmological constant. Phys.Lett., B477:285–291, Scenario: Phase Plane Analysis. Astrophys.Space Sci., 2000. 339:53–64, 2012. [7] Peter Bowcock, Christos Charmousis, and Ruth Gregory. [21] Pedro F. Gonzalez-Diaz. Quintessence in brane cosmol- General brane cosmologies and their global space-time ogy. Phys.Lett., B481:353–359, 2000. structure. Class.Quant.Grav., 17:4745–4764, 2000. [22] A.S. Majumdar. From brane assisted inflation to [8] Kei-ichi Maeda. Brane quintessence. Phys.Rev., quintessence through a single scalar field. Phys.Rev., D64:123525, 2001. D64:083503, 2001. [9] Shuntaro Mizuno, Kei-ichi Maeda, and Kohta Ya- [23] N.J. Nunes and Edmund J. Copeland. Tracking mamoto. Dynamics of scalar field in a brane world. quintessential inflation from brane worlds. Phys.Rev., Phys.Rev., D67:023516, 2003. D66:043524, 2002. [10] N. Yu. Savchenko and A.V. Toporensky. Scaling solutions [24] M. Sami and N. Dadhich. Unifying brane world inflation on a brane. Class.Quant.Grav., 20:2553–2562, 2003. with quintessence. TSPU Vestnik, 44N7:25–36, 2004. [11] Shinji Tsujikawa and M. Sami. A Unified approach to [25] Edmund J. Copeland, M. Sami, and Shinji Tsujikawa. scaling solutions in a general cosmological background. Dynamics of dark energy. Int.J.Mod.Phys., D15:1753– Phys.Lett., B603:113–123, 2004. 1936, 2006. [12] Ricardo Garcia-Salcedo, Tame Gonzalez, Claudia [26] Edmund J. Copeland, Andrew R Liddle, and David Moreno, and Israel Quiros. Randall-Sundrum brane cos- Wands. Exponential potentials and cosmological scaling mology: modification of late-time cosmic dynamics by solutions. Phys.Rev., D57:4686–4690, 1998. exotic matter. Class.Quant.Grav., 28:105017, 2011. [27] J.M. Aguirregabiria and Ruth Lazkoz. Tracking solutions [13] Adam G. Riess, Louis-Gregory Strolger, Stefano Caser- in tachyon cosmology. Phys.Rev., D69:123502, 2004. tano, Henry C. Ferguson, Bahram Mobasher, et al. New [28] Ruth Lazkoz, Genly Leon, and Israel Quiros. Quin- Hubble Space Telescope Discoveries of Type Ia Super- tom cosmologies with arbitrary potentials. Phys.Lett., novae at z >= 1: Narrowing Constraints on the Early B649:103–110, 2007. Behavior of Dark Energy. Astrophys.J., 659:98–121, 2007. [29] Wei Fang, Ying Li, Kai Zhang, and Hui-Qing Lu. Exact [14] Tamara M. Davis, E. Mortsell, J. Sollerman, A.C. Becker, Analysis of Scaling and Dominant Attractors Beyond the S. Blondin, et al. Scrutinizing Exotic Cosmological Mod- Exponential Potential. Class.Quant.Grav., 26:155005, els Using ESSENCE Supernova Data Combined with 2009. Other Cosmological Probes. Astrophys.J., 666:716–725, [30] Genly Leon. On the Past Asymptotic Dynamics of Non- 2007. minimally Coupled Dark Energy. Class.Quant.Grav., [15] W. Michael Wood-Vasey et al. Observational Constraints 26:035008, 2009. on the Nature of the Dark Energy: First Cosmological [31] Genly Leon, Pavel Silveira, and Carlos R. Fadragas. Results from the ESSENCE Supernova Survey. Astro- Classical and Quantum Gravity: Theory, Analysis and phys.J., 666:694–715, 2007. Applications, chapter Phase-space of flat Friedmann- [16] Max Tegmark et al. Cosmological parameters from SDSS Robertson-Walker models with both a scalar field cou- and WMAP. Phys.Rev., D69:103501, 2004. pled to matter and radiation. Nova Science Publishers, 12
New York, USA, 2012. arXiv:1009.0689 [gr-qc]. D9:373–444, 2000. [32] Genly Leon and Carlos R. Fadragas. Cosmological Dy- [43] L. Arturo Urena-Lopez and Tonatiuh Matos. A New namical Systems. LAP Lambert Academic Publishing, cosmological tracker solution for quintessence. Phys.Rev., Saarbr¨ucken, Germany, 2011. D62:081302, 2000. [33] Edmund J. Copeland, Seung-Joo Lee, James E. Lidsey, [44] Varun Sahni and Li-Min Wang. A New cosmologi- and Shuntaro Mizuno. Generalised cosmological scaling cal model of quintessence and dark matter. Phys.Rev., solutions. Phys.Rev., D71:023526, 2005. D62:103517, 2000. [34] Tame Gonzalez, Tonatiuh Matos, Israel Quiros, and Al- [45] Shuang-Yong Zhou. A New Approach to Quintessence berto Vazquez-Gonzalez. Self-interacting Scalar Field and Solution of Multiple Attractors. Phys.Lett., B660:7– Trapped in a Randall-Sundrum Braneworld: The Dy- 12, 2008. namical Systems Perspective. Phys.Lett., B676:161–167, [46] Rolando Cardenas, Tame Gonzalez, Yoelsy Leiva, Osmel 2009. Martin, and Israel Quiros. A model of the universe in- [35] Yoelsy Leyva, Dania Gonzalez, Tame Gonzalez, Tonatiuh cluding dark energy accounted for by both a quintessence Matos, and Israel Quiros. Dynamics of a self-interacting field and a (negative) cosmological constant. Phys.Rev., scalar field trapped in the braneworld for a wide variety of D67:083501, 2003. self-interaction potentials. Phys.Rev., D80:044026, 2009. [47] Philippe Brax and Jerome Martin. Quintessence and su- [36] David Langlois. Brane cosmology: An Introduction. pergravity. Phys.Lett., B468:40–45, 1999. Prog.Theor.Phys.Suppl., 148:181–212, 2003. [48] Philippe Brax and Jerome Martin. The Robustness of [37] Philippe Brax and Carsten van de Bruck. Cosmology and quintessence. Phys.Rev., D61:103502, 2000. brane worlds: A Review. Class.Quant.Grav., 20:R201– [49] T. Barreiro, Edmund J. Copeland, and N.J. Nunes. R232, 2003. Quintessence arising from exponential potentials. [38] David Langlois. Cosmology of brane - worlds. 2004. Phys.Rev., D61:127301, 2000. arXiv: astro-ph/0403579 [50] Scott Foster. Scalar field cosmologies and the initial [39] Roy Maartens. Brane world gravity. Living Rev.Rel., 7:7, space-time singularity. Class.Quant.Grav., 15:3485–3504, 2004. 1998. [40] Alan Coley. Dynamical Systems and Cosmology, vol- [51] Genly Leon, Yoelsy Leyva, Emmanuel N. Saridakis, Os- ume 291 of Astrophysics and Space Science Library. mel Martin, and Rolando Cardenas. Dark Energy: Theo- Dordrecht-Kluwer Publishers, Dordrecht, The Nether- ries, Developments, and Implications, chapter Falsifying lands, 2003. Field-based Dark Energy Models. Nova Science Publish- [41] R Tavakol. Dynamical Systems in Cosmology, chapter ers, New York, USA, 2010. arXiv:0912.0542 [gr-qc]. Introduction to dynamical systems, page 360. Cambridge [52] Naureen Goheer and Peter K.S. Dunsby. Exponential University Press, Cambridge, UK, 2005. potentials on the brane. Phys.Rev., D67:103513, 2003. [42] Varun Sahni and Alexei A. Starobinsky. The Case for a positive cosmological Lambda term. Int.J.Mod.Phys.,