• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    New exact solutions to the Eddingtoninspired-Born-Infeld equations for a relativistic perfect fluid in a Bianchi type I spacetime

    2022-08-02 02:48:10CalvinTadmonandGuichardDjiodjoSeugmo
    Communications in Theoretical Physics 2022年6期

    Calvin Tadmonand Guichard Djiodjo-Seugmo

    1 Department of Mathematics and Computer Science,Faculty of Science,University of Dschang,PO Box 67,Dschang,Cameroon

    2 The Abdus Salam International Centre for Theoretical Physics,Strada Costiera 11,I-34151 Trieste,Italy

    Abstract We consider a Bianchi type I physical metric g,an auxiliary metric q and a density matter ρ in Eddington-inspired-Born-Infeld theory.We first derive a system of second order nonlinear ordinary differential equations.Then,by a suitable change of variables,we arrive at a system of first order nonlinear ordinary differential equations.Using both the solution-tube concept for the first order nonlinear ordinary differential equations and the nonlinear analysis tools such as the Arzelá–Ascoli theorem,we prove an existence result for the nonlinear system obtained.The resolution of this last system allows us to obtain new exact solutions for the model considered.Finally,by studying the asymptotic behaviour of the exact solutions obtained,we conclude that this solution is the counterpart of the Friedman–Lema?tre–Robertson–Walker spacetime in Eddington-inspired-Born-Infeld theory.

    Keywords:Eddington-inspired-Born-Infeld theory,Bianchi type I spacetime,solution-tube,Friedmann–Lema?tre–Robertson–Walker spacetime,density matter,ordinary differential equations

    1.Introduction

    The very successful theory of General Relativity (GR) [1,2]can be validly considered as the foundation of many other theories of gravity.One more proof to confirm once again this gigantic theory is the recent appearance of gravitational waves[3].But perfection is not of this world.Some limitations of GR can be mentioned like the inescapable existence of singularities and the lack of effective explanations of phenomena such as dark matter and dark energy problems.In order to provide answers to some of these problems,several researchers have built alternative theories which deviate from the GR in the distributions of matter.Therefore,apart from quantum gravity,we have several modified gravities.As modified gravity,we have the one inspired by Born-Infeld(BI)electrodynamics[4],where the infiniteness of the electric field at the location of a point charge is regularized.Then,Deser and Gibbons[5]used this BI structure to set up a new theory of gravity in the metric formulation.It should not be forgotten that the pioneer Eddington had proposed a determining form of gravity in this reformulation of GR for a de Sitter spacetime [6],but the coupling of matter remained a problem in his approach.A few years later,Vollick [7]set up the Palatini formulation of BI gravity and worked on several axes,namely the non-trivial couplings of matter[8,9].Subsequently,Banados and Ferreira[10]proposed to the scientific world a new formulation,called Eddington-inspired-Born-Infeld (EiBI) theory of gravity,whose material coupling was different and simple than that proposed by Vollick in[7].The EiBI theory is a theory of bigravity[11–14],which reduces to GR theory in a vacuum.The main characteristic of the EiBI bi-gravity theory is the existence of a physical metric g which is coupled to matter and of an auxiliary metric q which is not coupled to matter.Throughout the work,the Universe described by the bi-metric(g,q) will be referred to as a bi-universe.The action of the EiBI theory proposed by [10]is given by

    where k is the Eddington parameter linked to the cosmological constant Λ by the dimensionless constant λ=1+kΛ;∣g∣denotes the absolute value of the determinant of the matrix(gαβ);SM[g,ψM]is the action of field.The equations of motion are obtained by varying the action(1)with respect to the fields gαβand,respectively,where

    This yields the following equations:

    with qβδ,αstanding for.Here Rαβis the symmetric part of the Ricci tensor built from the connectionΓαγβ,given by

    (qαβ)is the inverse of(qαβ),and the energy-momentum tensor is given by

    In GR,the energy-momentum tensor Tαβfulfils the following conservation equations,

    where ?is the covariant derivative.

    Many research activities have been carried out on EiBI theory in recent years.Astrophysical scenarios have been discussed in [15–20,46,22–28],spherically symmetric solutions have been obtained in [14,28–33],and some interesting cosmological solutions have been derived in[9,13,14,46,29,31,33–35].In the present paper,we assume the following equation of state for the fluid,p=wρ,where p,w and ρ are the pressure of the fluid,a constant which defines the fluid and the energy density of the fluid,respectively [14,36].We study the evolution of an uncharged perfect fluid of pure radiation type,i.e.whenw=,the geometry is a Bianchi type I spacetime,which is in fact a generalization of Friedmann–Lema?tre–Robertson–Walker (FLRW) spacetime [14,31,33].The motivation for choosing a radiation fluid is that it naturally shows an upper bound on its energy density,which further induces an upper bound for scalar curvature equivalent of the inverse Eddington parameter>0[14].Indeed,our motivation for considering the Bianchi type I geometry is that,in Physics,the Universe as we observe it today is very well described by the cosmological and homogeneous models of FLRW.We assume that the Universe is initially anisotropic and becomes asymptotically isotropic.Here we consider six functions a,b,c,A,B and C,of the time variable t.These functions are determined by the equations(3)–(4)of the EiBI theory,with the material content of spacetime represented by a material stress tensor defined by the density of material ρ.Equations of the EiBI theory,coupled with the laws of conservation,turn out to be a differential system allowing the determination of the functions a,b,c,A,B,C and ρ.By setting suitable initial data,we arrive at the solution of the model considered.We study the asymptotic behaviours of the volume of the scalar factorV,the Hubble parameter H,the Anisotropic parameterA and the deceleration parameter Q.Indeed,the study of these different parameters makes it possible to foresee the cosmological applications of the EiBI bi-gravity theory for homogeneous and isotropic backgrounds.In [14],the authors obtain and discuss the stabilities of different interesting cosmological solutions for early universe,including loitering,quasi de Sitter and bouncing solutions.Therefore,the anomalies observed in the cosmic microwave background(CMB)might be due to the presence of small anisotropies.Moreover,by studying asymptotic behaviours,we show that these solutions tend toward the vacuum at the late time under appropriate conditions.

    The present work is organized as follows.In section 2 we formulate the model by using the Bianchi type I geometry and consider a metric tensor g,an auxiliary tensor q together with a relativistic perfect fluid of density ρ in the Eddingtoninspired-Born-Infeld theory.In section 3,we obtain explicit expressions of functions a,b,c,A,B and C,and study the asymptotic behaviour of the bi-universe described by the physical metric g and the auxiliary metric q.The conclusion is elaborated in section 4.

    2.The EiBI equations

    We consider a Bianchi type I bi-universe(g,q,R4)and denote by xα=(x0,xi) the usual coordinates inR4;g and q stand for physical metric and auxiliary metric of signature(-,+,+,+)which can be written respectively as[14,31]:

    where a>0,b>0,c>0,A>0,B>0,and C>0 are unknown functions of the single time variable t=x0.

    The general expression of the stress matter tensor of a relativistic perfect fluid is given by

    where ρ>0 and p>0 are unknown functions of t representing the matter density and the pressure,respectively.We consider a perfect fluid of pure radiation type,which means thatp=.The matter tensor then can be written as:

    whereu=(uα)is a unit vector tangent to the geodesics of g.

    In order to simplify,we consider a frame in which the fluid is spatially at rest.This implies ui=ui=0.Recall that indices are raised and lowered using the usual rules:Vα=gαβVβ;Vα=gαβVβ,(gαβ)standing for the inverse matrix of(gαβ).

    Solving EiBI equations (3)–(4) consists of determining,on one hand,the gravitational field represented by the physical metricg=(gαβ)and auxiliary metricq=(qαβ),through the six unknown functions a>0,b>0,c>0,A>0,B>0,and C>0,and,on the order hand,its source represented by the matter density ρ.

    From the expressions of g and q given respectively by equations (8) and (9) we deduce that

    Using equations(2)and(13)we obtain by a direct calculation the non-vanishing Christofell coefficients,which are defined byΓii0andΓii0,i=1,2,3.So,we have

    Using the equation (5),we obtain the non-vanishing components of the Ricci tensor Rαβ,given by

    Using equations (11) and (12) we obtain the following nonvanishing contravariant components of the energy-momentum tensor:

    Note that we will have to consider both the evolution and the constraints systems.Therefore,we have the following two Propositions.

    Proposition 1.Letk≠ 0.The Eddington-inspired-Born-Infeld evolution system in(g,q,ρ)can be written as the following nonlinear second order system

    where a,b and c are given by (8),A,B and C are given by (9).

    Proof.From equations (4),(12) (13) and (15),a direct calculation gives equations (17)–(20).We now prove equation (21).It holds that

    Using the conservation law (7),we obtain the following equation

    Differentiatinguβ uβ=-1 yieldsuβ?α uβ=0.The contracted multiplication of equation (22) byβugives

    This is equivalent to

    Using equation (14),we get the following equation

    Sinceuα uα=-1 andui=ui=0,we have(u0)2=(u0)2=1.So u0is different from 0.From equation (23) we readily obtain equation (21).This ends the proof of proposition 1.□

    Proposition 2.Assuming thatλ+kρ>0 andλ->0.The Eddington-inspired-Born-Infeld constraints system can be written as

    where a,b and c are given by (8),A,B and C are given by (9).

    Proof.From equations (3),(12) (13) and (16),we straightforwarly get the result.□

    Lemma1.Letk≠0,λ+kρ>0 andλ->0.The Eddington-inspired-Born-Infeld evolution system(17)-(20) can be written as the following nonlinear second order system

    where a,b and c are given by (8),A,B and C are given by (9).

    Proof.Using equations (24)–(27) we have

    whereλ+kρ>0 andλ->0.Moreover,we have

    Therefore,by dividing equation (18) by A,equation (19) by B,equation (20) by C and using equation (32),we obtain equations (29),(30) and (31) respectively.□

    Lemma 2.We assume thatλ+kρ>0 andλ->0.Then the Eddington-inspired-Born-Infeld constraints system can be written as an equation of degree 3 for ρ,as follows:

    Proof.From equations (32) and (33) we readily obtain equation (34) by simple computation.□

    We are looking for solutions a,b,c,A,B,C,ρ,of the EiBI equations satisfying the initial conditions

    where a0>0,b0>0,c0>0 A0>0,B0>0,C0>0 ρ0>0 and Ci,i=1,2,3,4,5,6,are given real constants.Our aim is to prove the existence of solutions on [0,∞) of the above Cauchy problem.Equation(21)shows that the matter density ρ will be determined by

    which shows that ρ is known once A,B and C are known.It is well known that the EiBI constraint equation(34),is satisfied over the entire domain of the solutions a,b,c,A,B,C,ρ if and only if(34)is satisfied for t=t0,i.e.the initial data a0,b0,c0,A0,B0,B0,ρ0in (35) satisfy the initial constraint [37]

    By setting the following change of variables

    we get

    Hence,the evolution equations (28)–(31) can be written in terms of U and V.

    Lemma 3.Letk≠ 0.The Eddington-inspired-Born-Infeld evolution system (28)-(31) can be written as the following nonlinear first order system

    where U,V and W are given by (38),U˙,V˙andW˙are given by (39).

    Proof.By substituting equations(38)and(39)in equations(28),(29) (30) and (31),we obtain system (40)–(43).□

    We look for solutions(U,V,W) to system (40)–(43)satisfying at t=t0,the following initial conditions U(t0)=U0,V(t0)=V0,W(t0)=W0with

    Notice that system (28)–(31) imply (40)–(43),and,conversely,the knowledge ofgives A,B and C by direct integration.

    Theorem 1.Assume thatk> 0 and λ satisfy the inequalityand λ satisfy the inequalityλ.Then the initial constraint equation (37) admits a unique positive real solution,given by

    Proof.By settingequation (37) becomes

    According to the Cardan method,the discriminant of this polynomial of degree 3 is given byTherefore,we get the value ofγ0,given by

    Then,we obtain the value ofρ0>0for values of k and λ fixed under the prescribed assumptions.□

    According to equation (37),we haveλ->0and λ+kρ0>0,this makes it possible to determine the domains of the Eddington parameter k and the dimensionless constant λ,given by the following equation

    Lemma 4.Letk>0 andork<0 andλ∈ (-kρ0,∞).The Eddington-inspired-Born-Infeld evolution system (40)-(43) can be written as the following nonlinear first order equation

    whereχ=U+V+Wandρ0is given by (44).

    Proof.By summing up equations (41)–(43) and using (40),we get

    Moreover (40) gives us

    By settingχ=U+V+Wand using (36) and (49),equation (50) gives (48).□

    To advance with our analysis,we will first establish an existence theorem of the solution χ for the nonlinear ODE (48).

    There is a variety of methods that can be used to solve equation (48).We can cite:(1) the solution-tube method,(2)the structure-preserving method,(3) the Lie group method,(4) the symmetry method,(5) the multiplication method,(6)the generalized multi-symplectic method,to name but the few.The solution-tube method has been used in [38]for solving systems of ordinary differential and time-scale equations.In [39],Hu and his collaborators used the preserving structure method to design a strategy for removing space debris by a tether tow.The Lie group method has been developed in [40]to efficiently solve a nonlinear modified Gardner type partial differential equation (PDE) and its timefractional form.In [41],the authors used the symmetry method to solve generalized KdV-Burgers-Kuramoto equations and its fractional version.The multiplication method has been developed in[42]to find analytical solutions of a (3+1)-dimensional sine-Gordon and a sinh-Gordon equations.In[43,44],the authors combined the Runge-Kutta numerical method to a generalized multi-symplectic method for solving some PDEs arising from physical phenomena such as wave propagation and signal processes.

    In this work,we will mainly use the solution-tube method.

    For this purpose,we need the following notations.

    (N-1)C1(R,R),the set of differentiable functions onR such thatχ˙ is continuous onR.This space is endowed with the normwhereis the space of measurable functions χ such that |χ| is integrable.

    We now solve the following problem

    where T>t0and F is the function defined by

    We will find the solution to problem (51) inC1(R,R).We now give the necessary mathematical tools,for the resolution of problem (51).

    Let’s introduce a notion of solution-tube for this problem.

    Definition 1.Let (v,M) ∈C1(R,R)×C1(R,[0,∞)).(v,M)is called a solution-tube of (51) if the following conditions are satisfied:

    (iii)∣v(T)-v(t0)∣≤M(T)-M(t0).

    The results below will help to establish the existence theorem for problem (51).

    Theorem 2.[38]Let W be an open set ofR and t be a real number.Ifg:R—R is differentiable in t and iff:W—R is differentiable ing(t) ∈W,thenf°gis differentiable in t and

    Theorem 2 gives the following example.

    Example 1.Let a functionχ:R—R differentiable int∈R.We know that the function∣·∣:R-{0} — (0,∞)is differentiable.If∣χ(t)∣>0,then by continuity,there existsδ>0 such that∣χ(s)∣>0 fors∈ (t-δ,t+δ).Moreover,by virtue of theorem 2 we have

    The remark below allows us to define explicitly the exponential function which is of capital utility for the proof of the existence theorem of the Problem (51).

    Remark 1.For∈>0,the exponential function e∈(·,t1):R —R can be defined as being the unique solution to the problem with the initial value

    More explicitly,the exponential function e∈(·,t1) can be given by the formula

    where forh≥ 0,we defineζ∈(h) by

    This remark also allows us to state the following Theorem.

    Theorem 3.[38]Ifg∈L1(R,R),the functionχ:R—R defined by

    is a solution to the problem

    Lemma 5.Let a functionL∈C1(R,R)such thatL˙ (t)>0 for allt∈ {s∈ R:L(s)>0}.IfL(t0)≥L(T),thenL(t) ≤0,for allt∈R.

    Proof.Suppose that there existst∈Rsuch thatL(t)>0,then there existst1∈R such thatL(t1)=maxt∈RL(t)>0.

    – Ift1<T,then there exists an interval[t1,t2]such thatL(t)>0 for allt∈[t1,t2].So,L(t2)-L(t1),which contradicts the maximality ofL(t1).

    – Ift1=T,then by hypothesis of lemma 5,it would be necessary thatL(t0)=L(T).

    – By takingt1=t0,by what precedes,we would find thatL(t0) ≤0.□

    Let us impose the following hypothesis:

    (H) There exists (v,M) a solution-tube of (51).

    In order to prove the existence Theorem,we will have recourse to the following modified problem

    where

    We define the operatorTP*:C(R,R) —C(R,R)by

    where the function e1(·,t0) is defined by (54).

    We can easily see that for k ≠0,λ->0and λ+kρ>0,the functionF:R×R—R defined in (52) is continuous.Therefore we state the following Proposition.

    Proposition 3.LetF:R×R—R be a continuous function.If the hypothesis(H) is satisfied,then the operatoris compact.

    Proof.This result is established in two main steps.

    Step 1.Let us first show the continuity of the operatorTP*.Let {χn}n∈?be a sequence inC(R,R)converging towards an elementχ∈C(R,R).

    Using equation (60) and after some calculations,we get the following inequalities

    Step 2.Let us now show that the setTP*(C(R,R))is relatively compact.Consider a sequence {yn}n∈?inTP*(C(R,R)).For alln∈?,there exitsχn∈C(R,R)such thatyn=TP*(χn).Using equation (60),we get the following inequality

    By definition,there existsR>0 such that∣(s)∣≤Rfor alls∈R and alln∈?.The function F being uniformly continuous overR×BR(0),we can deduce the existence of a constantA>0 such that∣F(s,(s))∣≤Afor alls∈Randn∈?.Thus,the sequence {yn}n∈?is uniformly bounded.Note also from the above fact that fort1,t2∈R,we have

    where B is a constant which can be chosen such that it is greater than

    Thus the sequence {yn}n∈?is also equicontinuous and by virtue of the Arzelà–Ascoli Theorem,{yn}n∈?has a convergent subsequence.

    Steps 1 and 2 show thatTP*is compact.□

    We can now state and prove the existence Theorem.

    Theorem 4.If the hypothesis(H) is satisfied,then problem(51) has a solutionχ∈C1(R,R) ∩T(v,M),

    where

    Proof.Step 1.We need to show that a.e.on the setA={t∈ R:∣χ(t)-v(t)∣>M(t)},we have

    This is done by using Example 1.

    Step 2.We will now show that almost everywhere on A,we have

    ?IfM(t)>0,then by the hypothesis of solution-tube,and by (98),almost everywhere,we have

    ?IfM(t)=0,by the hypothesis of the solution-tube and by (98),a.e.we have that

    Hence,almost everywhere on A,we have

    Step 3.By settingL(t)=∣χ(t)-v(t)∣-M(t),we have for allt∈ {t∈ R:L(t)>0},(t)>0.Moreover,by the hypothesis of the solution-tube,notice that

    Thus,the assumptions of Lemma 5 are satisfied.This proves Theorem 4.□

    As F is continuous,we also have the following result.

    Theorem 5.[38]LetF:R×R—R be a continuous function.If there exists non-negative constants α and K such that

    for all (t,p) ∈ R×R whereh:R—R is defined bywith

    Then the problem (51) has a solution.

    In particularF:{0,1,...,N,N+1}×R —R (N being a positive integer) is continuous,theorem 5 yields another important result.

    Theorem 6.[38]LetF:{0,1,...,N,N+1}×R —R be a continuous function.If there exists non-negative constants α and K such that

    for all (t,p) ∈{0,1,...,N,N+1}×R.Then the following system of finite difference equations

    has a solution.

    Since we have just shown that χ exists,then equations (41)–(43) of lemma 3 admit solutions given by the following Lemma.

    Lemma 6.Letk>0 andλ∈ (,∞)ork<0 andλ∈ (-kρ0,∞).The Eddington-inspired-Born-Infeld evolution system (41)-(43) admits a solution(U,V,W)given by

    whereC1=A0U0,C2=B0V0andC3=C0W0are given real constants,i s given by (44) andA0>0,B0>0andC0>0.

    Proof.Since χ exists,equations(41)–(43)can be rewritten as follows

    which is a linear system of ODEs with non-constant coefficients.Consequently,by using the constant variation method one arrives at the solution given by (64)–(66).□

    We have just proved that the auxiliary metric q exists.Indeed,according to the above discussion,we obtain the component functions A,B and C which are given respectively by

    where U,V and W are defined respectively by (64),(65)and (66).

    Now,we are going to prove the existence of the physical metric g.Using equation (32),the expressions of a,b and c are given by

    The following theorem gives a general summary of the existence of the physical metric g,the auxiliary metric q and the matter density ρ.

    Theorem 7.Letk>0 andλ∈ork<0 andλ∈ (-kρ0,∞).Assume thata0>0,b0>0andc0>0.Then,there exists a solution(g,q,ρ)given by

    where

    and U,V and W are defined respectively by (64),(65)and (66).

    We now study the asymptotic behaviour of the solutions obtained in Theorem 7.

    3.Asymptotic behaviour

    From values of the functions a,b,c,A,B,C and ρ,we now define important cosmological indicators in the directions of the physical metric g and the auxiliary metric q.

    The volumes of scalar factorsVqandVgare defined respectively by

    The Hubble’s functions Hqand Hgare defined respectively by

    where χ0=U0+V0+W0=χ(t0).

    It is of great cosmological importance to evaluate the anisotropy parameters Agrelated to the physical metric g and Aqrelated to the auxiliary metric q.Indeed,these different parameters are useful indicators to study the behaviour of anisotropic cosmological models.In GR,these parameters are finite for singular states [31].

    We have

    Similarly

    We presently define the decelerations parameters Qgand Qqin the directions of g and q,respectively.The role of these parameters in Cosmology is to know when the bi-universe described by metrics g and q inflate or not.We have

    Remark 2.

    (i) Assuming thatU(t)U0+V(t)V0+W(t)W0<we can easily show thatχ˙ (t)<0.ThenQq(t)≥0andQg(t)≥0.So,in this case,there is standard deceleration of the bi-universe in the direction of the auxiliary metric q and in the direction of the physical metric g.

    (ii) Assuming thatU(t)U0+V(t)V0+W(t)W0≥we can easily show thatχ˙ (t)≥0.ThenQq(t) <0andQg(t) <0.So the bi-universe inflate in the direction of the auxiliary metric q and in the direction of the physical metric g.

    (iii) By choosingχ(t)=withχ0=0,we readily getQq(t)=Qg(t)=-1,for allt∈R.Hence,one obtains de Sitter phase in both directions at early and at late time bi-universe.

    (iv) Fort=t0,notice that when(U0+V0+W0)2,we haveV0W0+U0V0+U0W0=0 and Aq(t)=2.Therefore the auxiliary metric q describes the Kasner geometry (see [31]).

    (v) Assuming thatA0=B0=C0,we haveA(t)=B(t)=Then the physical metric g and the auxiliary metric q are given by

    So in both directions of the bi-universe,the physical metric g and auxiliary metric q describe the FLRW metrics with cosmological expansion factorsandrespectively.

    (vi) Assuming that att=t0,we haveA0=B0=C0=1.Then the physical metric g and the auxiliary metric q are given by

    which describe Minkowski spacetimes.

    (vii) Fort=t0,the volumes of scalar factors satisfy the initial conditionsVq(t0)=A0B0C0≥8πkρ0andTherefore,in EiBI bi-gravity the stiff causal bi-universe starts its evolution from a non-singular state (see [31]).

    (viii) The asymptotic behaviours of Hg2,Hq2and the time derivativeof the Hubble rate in the EiBI regime are

    Similarly,we have

    So,the solutions obtained are unstable when cosmic time t tends to t0(see figure 37 in [14]).

    (ix) It is worth noting that whenλ=1,both physical and auxiliary metrics g and q coincide,and the bi-universe described by them is asymptotically flat [14].

    In theorem 8 below,we summarize the asymptotic behaviour of the physical metric g and the auxiliary metric q.

    Theorem 8.Letor(k<0 andλ∈ (-kρ0,∞)),A0>0,B0>0andC0>0.Assume thatρ(t)0.Then the bi-universe(R4,gαβ,qαβ,Tαβ)tends towards the vacuum at late time.

    Proof.We just have to show that

    According to the hypothesis of theorem 8 and using equation (74),we getA(t)∞,B(t)∞andC(t)∞.By direct calculation,we have

    We notice that(73)holds forα=β.But forα≠β,we haveT αβ(t)=0.So the result follows.□

    4.Conclusion

    In this work,we were in the quest of finding new exact solutions of Bianchi type I,comprising a relativistic perfect fluid,in Eddington-inspired-Born-Infeld theory of gravity.For this purpose,we first established a system of second order ordinary differential equations (Proposition 1),supplemented by a suitable constraints system(Proposition 2).Then,after an appropriate change of variables,combined with a thorough analysis of the systems obtained,we arrived at a nonlinear first order ordinary differential equation (see equation (48) of Lemma 4).To establish a solution to equation (48),we consider problem (51).By using the mathematical tools of nonlinear analysis and the notion of solution-tube for nonlinear first order differential equations,we succeed in establishing an existence result (Theorem 4) for problem (51).From this solution,we retrieve the physical metric g,the auxiliary metric q and the matter density ρ (Theorem 7).Finally,by dint of many meticulous calculations,we investigated the asymptotic behaviours of the physical metric g and the auxiliary metric q obtained.As a matter of fact,we provided the expressions of the volumes of scalar factors,Hubble’s functions,anisotropy parameters and deceleration parameters related to both metrics.These aforementioned functions are very important in Cosmology.We saw from Remark 4 that the auxiliary metric q is equivalent to the Kasner metric when the anisotropy parameter Aq(t)equals 2 for t=t0so thatU02+V02+W02=(U0+V0+W0)2.We also found from equations (89) and(90) that in both directions the solutions obtained are unstable when the cosmic time t tends to t0.We also noticed that when λ=1 both physical metric g and auxiliary metric q coincide.From Theorem 8,we concluded that,under appropriate conditionsor(k<0 andλ∈ (-kρ0,∞)),A0>0,B0>0 and C0>0,when t tends towards a late time the solution obtained in the EiBI theory tends towards a vacuum.

    It is worth stressing the fact that,in the quest for searching exact solutions to the model studied in this work,we have set the lapse function to 1 in(8)and(9).However,it would be interesting to consider the cases where one of the lapse function is 1 and the other one is unknown.

    Firstly,let’s consider the physical metric g given by equation(8)and the auxiliary metric q given by the following equation

    where N>0 represents the lapse function.In this case the non-vanishing components of the Ricci tensor are now given by

    After lengthy and tedious calculations,we arrive at the following equations which are counterparts of equations (17)–(20)

    It should be noted that equation (94) is the counterpart of equation (17) in this case.So there is no way to avoid this.Also,equations (94)–(97) are quite involved and not easy to handle with the approach used for solving equations (17)–(20) in the present work.The nonlinear second order ODEs system (94)–(97) is difficult to manipulate for the search of exact solutions.

    Secondly,we now consider the auxiliary metric defined by (9) and the physical metric defined as follows:

    where M>0 represents the lapse function.In this case,direct calculations show that only equation (17) of Proposition 1 is changed to its following analogue:

    where M verifies the following equation

    So,in all cases,equation (17) and its analogues cannot be avoided.We postpone the complete study of these two cases to a forthcoming paper.Another perspective pertaining to the present work is to take into consideration the spatial dependence of the physical and the auxiliary metrics.This would yield a more intricate system of partial differential equations to analyze mathematically and interpret physically.One could begin with the spherical symmetry assumption on g and q.

    Acknowledgments

    Calvin Tadmon thanks Professor Burin Gumjudpai for a discussion on EiBI gravity.

    ORCID iDs

    99九九线精品视频在线观看视频| 人妻少妇偷人精品九色| 国产黄a三级三级三级人| 看非洲黑人一级黄片| 国产精品福利在线免费观看| 国内精品美女久久久久久| 亚洲经典国产精华液单| 久久精品国产亚洲av涩爱| 久久精品国产亚洲av天美| 国产美女午夜福利| 国产免费一区二区三区四区乱码| 国产高清国产精品国产三级 | 三级国产精品欧美在线观看| 国产欧美亚洲国产| 搞女人的毛片| 精品一区二区三卡| 亚洲图色成人| 色视频www国产| 好男人视频免费观看在线| 听说在线观看完整版免费高清| 国产精品人妻久久久影院| 七月丁香在线播放| 日韩欧美精品v在线| 在线天堂最新版资源| 亚洲性久久影院| a级一级毛片免费在线观看| 在线观看一区二区三区激情| 亚洲不卡免费看| 精品人妻熟女av久视频| 熟女人妻精品中文字幕| 国产淫片久久久久久久久| 亚洲va在线va天堂va国产| 国模一区二区三区四区视频| 国产成人a∨麻豆精品| 一个人观看的视频www高清免费观看| 久久97久久精品| 69人妻影院| 国产av国产精品国产| 欧美日韩视频高清一区二区三区二| 日韩强制内射视频| 久久久精品94久久精品| 欧美另类一区| 五月玫瑰六月丁香| 国产一区亚洲一区在线观看| 七月丁香在线播放| 亚洲综合色惰| 久久久久久伊人网av| 亚洲欧美日韩无卡精品| 国产免费一级a男人的天堂| 大话2 男鬼变身卡| 寂寞人妻少妇视频99o| 日本-黄色视频高清免费观看| 国产免费福利视频在线观看| 成人毛片60女人毛片免费| 国产女主播在线喷水免费视频网站| 亚洲av欧美aⅴ国产| 国产乱人视频| 日日啪夜夜撸| 免费高清在线观看视频在线观看| 精品久久久久久久久亚洲| 草草在线视频免费看| 男女下面进入的视频免费午夜| www.色视频.com| 亚洲伊人久久精品综合| 日韩三级伦理在线观看| 777米奇影视久久| 国产黄a三级三级三级人| 在线观看国产h片| 欧美日韩视频高清一区二区三区二| 69av精品久久久久久| 国产免费一区二区三区四区乱码| 高清欧美精品videossex| 丝袜脚勾引网站| 天堂俺去俺来也www色官网| 国产免费一区二区三区四区乱码| 男插女下体视频免费在线播放| 男男h啪啪无遮挡| 能在线免费看毛片的网站| 亚洲国产最新在线播放| 六月丁香七月| 一级毛片我不卡| 老司机影院成人| 狂野欧美激情性bbbbbb| 国产中年淑女户外野战色| 99热全是精品| 亚洲精品456在线播放app| 亚洲美女搞黄在线观看| 免费黄频网站在线观看国产| 国产淫语在线视频| 18+在线观看网站| 亚洲,欧美,日韩| 在线精品无人区一区二区三 | 黄色欧美视频在线观看| 免费人成在线观看视频色| 尾随美女入室| 天天一区二区日本电影三级| 亚洲经典国产精华液单| 日本黄大片高清| freevideosex欧美| 国产午夜精品久久久久久一区二区三区| 国产精品偷伦视频观看了| 亚洲成人中文字幕在线播放| 国产精品秋霞免费鲁丝片| 18禁裸乳无遮挡动漫免费视频 | 国产高清三级在线| a级一级毛片免费在线观看| 麻豆成人av视频| 白带黄色成豆腐渣| 肉色欧美久久久久久久蜜桃 | 毛片女人毛片| 国产精品久久久久久精品电影小说 | 黄色一级大片看看| 美女脱内裤让男人舔精品视频| 国内精品美女久久久久久| 自拍偷自拍亚洲精品老妇| 国内揄拍国产精品人妻在线| 亚州av有码| 九色成人免费人妻av| 69人妻影院| 久久久久精品性色| 少妇熟女欧美另类| 51国产日韩欧美| 欧美 日韩 精品 国产| 丝瓜视频免费看黄片| 国产精品麻豆人妻色哟哟久久| 黄色欧美视频在线观看| 尾随美女入室| 欧美区成人在线视频| 国产黄色免费在线视频| 3wmmmm亚洲av在线观看| 七月丁香在线播放| av一本久久久久| 免费观看性生交大片5| 高清欧美精品videossex| 少妇 在线观看| 狠狠精品人妻久久久久久综合| 老师上课跳d突然被开到最大视频| 成人一区二区视频在线观看| 欧美人与善性xxx| 国产精品一区二区三区四区免费观看| 亚洲色图av天堂| 欧美精品一区二区大全| 久久国产乱子免费精品| 精品一区二区免费观看| 国产成人精品婷婷| 青春草亚洲视频在线观看| 久久精品久久精品一区二区三区| 成人无遮挡网站| 国产精品嫩草影院av在线观看| 三级国产精品片| 又爽又黄a免费视频| 久久这里有精品视频免费| 国产精品不卡视频一区二区| 人体艺术视频欧美日本| 亚洲一区二区三区欧美精品 | 韩国av在线不卡| 伊人久久国产一区二区| 欧美日本视频| 国产精品一区二区在线观看99| 亚洲欧美一区二区三区黑人 | 国产黄频视频在线观看| 最近最新中文字幕免费大全7| 国国产精品蜜臀av免费| 亚洲高清免费不卡视频| 国产 精品1| 亚洲国产精品国产精品| 午夜免费男女啪啪视频观看| 成人亚洲精品一区在线观看 | 国产精品国产av在线观看| 国产白丝娇喘喷水9色精品| 91久久精品电影网| 欧美日韩精品成人综合77777| 久久久国产一区二区| 国产精品av视频在线免费观看| 爱豆传媒免费全集在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 丝袜美腿在线中文| 国产 精品1| 国产精品麻豆人妻色哟哟久久| 亚洲性久久影院| 高清午夜精品一区二区三区| 2022亚洲国产成人精品| 草草在线视频免费看| 欧美bdsm另类| a级毛片免费高清观看在线播放| 校园人妻丝袜中文字幕| 99久久精品国产国产毛片| 少妇熟女欧美另类| 18禁动态无遮挡网站| 久久久久久久亚洲中文字幕| 色5月婷婷丁香| 肉色欧美久久久久久久蜜桃 | 老司机影院成人| 婷婷色av中文字幕| 真实男女啪啪啪动态图| 青春草视频在线免费观看| 搞女人的毛片| 国产女主播在线喷水免费视频网站| 国产黄频视频在线观看| 我的女老师完整版在线观看| av一本久久久久| 大片免费播放器 马上看| 亚洲欧美一区二区三区国产| 国产日韩欧美亚洲二区| 国产成人一区二区在线| 国产熟女欧美一区二区| 深夜a级毛片| 午夜激情久久久久久久| 啦啦啦中文免费视频观看日本| 国产精品一区www在线观看| 91精品国产九色| 99久久精品一区二区三区| 九九在线视频观看精品| 亚洲熟女精品中文字幕| 国产成人福利小说| 一二三四中文在线观看免费高清| 国产精品国产三级国产专区5o| 最近最新中文字幕大全电影3| 少妇高潮的动态图| 亚洲av中文av极速乱| 精品久久国产蜜桃| 久久久久久久午夜电影| 免费看a级黄色片| 午夜激情久久久久久久| 在线免费观看不下载黄p国产| 欧美变态另类bdsm刘玥| 亚洲激情五月婷婷啪啪| av在线天堂中文字幕| 亚洲av不卡在线观看| 一个人看视频在线观看www免费| 亚洲欧美精品专区久久| 国产成年人精品一区二区| 精品久久久精品久久久| 午夜精品一区二区三区免费看| 亚洲激情五月婷婷啪啪| 美女xxoo啪啪120秒动态图| 麻豆乱淫一区二区| 伊人久久精品亚洲午夜| 免费播放大片免费观看视频在线观看| 亚洲精品国产色婷婷电影| 国产黄色视频一区二区在线观看| 色5月婷婷丁香| 亚洲欧美日韩无卡精品| 亚洲精品影视一区二区三区av| 97在线人人人人妻| 纵有疾风起免费观看全集完整版| 狂野欧美激情性xxxx在线观看| 欧美xxxx性猛交bbbb| 日韩在线高清观看一区二区三区| 精品人妻视频免费看| 亚洲欧美日韩卡通动漫| 久久久久久久国产电影| 看十八女毛片水多多多| av播播在线观看一区| 久久99热这里只有精品18| 亚洲人成网站在线观看播放| 国产精品久久久久久久电影| 久久久久精品久久久久真实原创| 午夜亚洲福利在线播放| 久久久久国产精品人妻一区二区| 大又大粗又爽又黄少妇毛片口| 久久久久久久久久久免费av| 日本午夜av视频| 一级毛片电影观看| 最后的刺客免费高清国语| 久久久久九九精品影院| 内地一区二区视频在线| 91狼人影院| 2021天堂中文幕一二区在线观| 大片电影免费在线观看免费| 丰满少妇做爰视频| 蜜桃亚洲精品一区二区三区| 久久久欧美国产精品| 日韩av不卡免费在线播放| 高清av免费在线| 简卡轻食公司| 国产成人a区在线观看| 一个人看的www免费观看视频| 国产一区二区在线观看日韩| 日本爱情动作片www.在线观看| 国产精品一二三区在线看| 熟女av电影| 午夜福利在线在线| 99九九线精品视频在线观看视频| 少妇熟女欧美另类| 看黄色毛片网站| 纵有疾风起免费观看全集完整版| 高清av免费在线| 狂野欧美激情性bbbbbb| 少妇的逼水好多| 亚洲精品,欧美精品| 美女xxoo啪啪120秒动态图| 另类亚洲欧美激情| 日韩一区二区视频免费看| 麻豆久久精品国产亚洲av| 一本色道久久久久久精品综合| videossex国产| 国产在线一区二区三区精| 九九在线视频观看精品| 国产欧美另类精品又又久久亚洲欧美| 成人无遮挡网站| 国产大屁股一区二区在线视频| 免费av不卡在线播放| 又大又黄又爽视频免费| 18禁在线无遮挡免费观看视频| 亚洲美女搞黄在线观看| 国产av码专区亚洲av| 高清日韩中文字幕在线| 欧美一区二区亚洲| 男人狂女人下面高潮的视频| 久久久久久久久久久免费av| 久久久成人免费电影| 97热精品久久久久久| 日日啪夜夜爽| 少妇高潮的动态图| 免费不卡的大黄色大毛片视频在线观看| 高清欧美精品videossex| 精品视频人人做人人爽| 22中文网久久字幕| 亚洲av电影在线观看一区二区三区 | 一级毛片电影观看| 精品人妻熟女av久视频| 波多野结衣巨乳人妻| 2022亚洲国产成人精品| 久久人人爽av亚洲精品天堂 | 国产精品国产三级国产专区5o| 国产精品不卡视频一区二区| 久久久午夜欧美精品| 观看美女的网站| 中文字幕亚洲精品专区| 国产精品.久久久| 国产精品av视频在线免费观看| 最近手机中文字幕大全| 国产亚洲av片在线观看秒播厂| 三级男女做爰猛烈吃奶摸视频| 中文精品一卡2卡3卡4更新| 午夜爱爱视频在线播放| 成年女人在线观看亚洲视频 | 久久国内精品自在自线图片| 久久久a久久爽久久v久久| 高清午夜精品一区二区三区| 黑人高潮一二区| 最近的中文字幕免费完整| 免费av观看视频| 少妇被粗大猛烈的视频| 精品久久久久久久末码| 老女人水多毛片| 国产乱人视频| 国产黄片视频在线免费观看| 伊人久久国产一区二区| 交换朋友夫妻互换小说| 久久久久久国产a免费观看| 国产乱来视频区| 亚洲欧美日韩东京热| 亚洲欧美成人综合另类久久久| 毛片女人毛片| 99久久精品一区二区三区| 亚洲成人精品中文字幕电影| 国精品久久久久久国模美| 久久精品国产a三级三级三级| 欧美精品一区二区大全| 国产综合精华液| av黄色大香蕉| 男的添女的下面高潮视频| 在线观看美女被高潮喷水网站| 熟妇人妻不卡中文字幕| 国产午夜福利久久久久久| 欧美少妇被猛烈插入视频| 日本黄色片子视频| h日本视频在线播放| 爱豆传媒免费全集在线观看| 少妇 在线观看| 91精品一卡2卡3卡4卡| av免费观看日本| 禁无遮挡网站| 小蜜桃在线观看免费完整版高清| 五月伊人婷婷丁香| 国产亚洲av嫩草精品影院| 久久久欧美国产精品| 亚洲欧洲国产日韩| 最后的刺客免费高清国语| 小蜜桃在线观看免费完整版高清| 久久国产乱子免费精品| 最近手机中文字幕大全| 波多野结衣巨乳人妻| 欧美日韩在线观看h| 丰满少妇做爰视频| 一级毛片aaaaaa免费看小| 深夜a级毛片| 2021少妇久久久久久久久久久| 久久久精品免费免费高清| 日产精品乱码卡一卡2卡三| 好男人视频免费观看在线| 欧美变态另类bdsm刘玥| 美女主播在线视频| 能在线免费看毛片的网站| 成年女人在线观看亚洲视频 | 久久久久精品性色| 色5月婷婷丁香| 亚洲一级一片aⅴ在线观看| 国产亚洲一区二区精品| 黄色日韩在线| 亚洲精品国产成人久久av| 免费黄网站久久成人精品| 国产高潮美女av| 亚洲,一卡二卡三卡| 成人欧美大片| 在线观看一区二区三区| 日本欧美国产在线视频| 熟女av电影| 日韩三级伦理在线观看| 热re99久久精品国产66热6| av.在线天堂| 午夜免费鲁丝| av线在线观看网站| 久久热精品热| 成人无遮挡网站| 中文资源天堂在线| 精品一区二区免费观看| 久久久久性生活片| 99视频精品全部免费 在线| 精品久久国产蜜桃| 18禁动态无遮挡网站| 22中文网久久字幕| 97热精品久久久久久| 中国美白少妇内射xxxbb| 99久久中文字幕三级久久日本| 免费在线观看成人毛片| 一级毛片黄色毛片免费观看视频| 最近最新中文字幕免费大全7| 青青草视频在线视频观看| 欧美zozozo另类| 男女那种视频在线观看| av播播在线观看一区| 好男人在线观看高清免费视频| 少妇高潮的动态图| 噜噜噜噜噜久久久久久91| 真实男女啪啪啪动态图| 看非洲黑人一级黄片| 99热全是精品| 亚洲av日韩在线播放| 美女内射精品一级片tv| 欧美激情国产日韩精品一区| 亚洲欧洲日产国产| 国产日韩欧美亚洲二区| 日本wwww免费看| 国产成人一区二区在线| 97超视频在线观看视频| 亚洲国产成人一精品久久久| 边亲边吃奶的免费视频| 亚洲精品国产av成人精品| 黄色一级大片看看| 综合色丁香网| 久久久精品免费免费高清| 国产精品三级大全| 久久这里有精品视频免费| 欧美另类一区| 人人妻人人看人人澡| 成人特级av手机在线观看| 国产日韩欧美在线精品| 丰满人妻一区二区三区视频av| 亚洲精品一区蜜桃| 亚洲av成人精品一二三区| 欧美高清性xxxxhd video| 91久久精品电影网| 好男人视频免费观看在线| 不卡视频在线观看欧美| 男女国产视频网站| 国产成人精品一,二区| 亚洲国产av新网站| 老司机影院成人| 成人欧美大片| 国产午夜精品久久久久久一区二区三区| 一二三四中文在线观看免费高清| 精品99又大又爽又粗少妇毛片| 日本免费在线观看一区| 国产亚洲最大av| 国产精品伦人一区二区| 在线a可以看的网站| 国产男女内射视频| 亚洲精品中文字幕在线视频 | 日韩电影二区| 亚洲精品成人久久久久久| 日本欧美国产在线视频| av在线亚洲专区| 九九久久精品国产亚洲av麻豆| 久久久久网色| 高清视频免费观看一区二区| 一本色道久久久久久精品综合| 久久精品国产a三级三级三级| 国产真实伦视频高清在线观看| 一级黄片播放器| 日韩av不卡免费在线播放| 免费观看的影片在线观看| 久久女婷五月综合色啪小说 | 欧美成人午夜免费资源| 国产一区二区三区av在线| 一级黄片播放器| 久久久午夜欧美精品| 少妇熟女欧美另类| 只有这里有精品99| 最近的中文字幕免费完整| 欧美精品人与动牲交sv欧美| 天天躁日日操中文字幕| 五月开心婷婷网| 亚洲aⅴ乱码一区二区在线播放| 久久久久国产网址| 亚洲欧美日韩卡通动漫| 欧美激情国产日韩精品一区| 精品少妇久久久久久888优播| 麻豆成人av视频| 校园人妻丝袜中文字幕| 赤兔流量卡办理| 精品久久国产蜜桃| 久久久精品免费免费高清| 国产国拍精品亚洲av在线观看| av在线老鸭窝| 久久久成人免费电影| 精品亚洲乱码少妇综合久久| 草草在线视频免费看| 亚洲精品一区蜜桃| 高清日韩中文字幕在线| 久久久a久久爽久久v久久| av又黄又爽大尺度在线免费看| 国产 一区 欧美 日韩| 少妇人妻久久综合中文| 97在线人人人人妻| 国内揄拍国产精品人妻在线| 亚洲久久久久久中文字幕| 男女下面进入的视频免费午夜| 国产久久久一区二区三区| 成人国产麻豆网| 高清毛片免费看| 一区二区av电影网| 伊人久久国产一区二区| 美女国产视频在线观看| 涩涩av久久男人的天堂| 成人国产av品久久久| 亚洲精品成人久久久久久| 久久久久久久大尺度免费视频| 少妇人妻一区二区三区视频| 欧美成人精品欧美一级黄| 欧美日韩精品成人综合77777| 成年av动漫网址| 日本wwww免费看| 国产亚洲av片在线观看秒播厂| 在线观看三级黄色| 久久这里有精品视频免费| 激情 狠狠 欧美| 高清视频免费观看一区二区| 91在线精品国自产拍蜜月| 成人亚洲精品一区在线观看 | 国产黄片视频在线免费观看| 内射极品少妇av片p| 免费观看无遮挡的男女| 三级男女做爰猛烈吃奶摸视频| 午夜老司机福利剧场| 在线天堂最新版资源| 亚洲精品日本国产第一区| 天天躁夜夜躁狠狠久久av| 人妻系列 视频| 国产一区二区三区av在线| 亚洲精品国产av蜜桃| 日韩中字成人| 黄色视频在线播放观看不卡| 夫妻性生交免费视频一级片| 五月天丁香电影| 国产黄色免费在线视频| 美女高潮的动态| 日韩伦理黄色片| 九九爱精品视频在线观看| 国产极品天堂在线| 好男人在线观看高清免费视频| 国产精品精品国产色婷婷| kizo精华| 久久久久久久精品精品| 精品酒店卫生间| 一区二区三区乱码不卡18| 夜夜看夜夜爽夜夜摸| 在线免费十八禁| 国产真实伦视频高清在线观看| 涩涩av久久男人的天堂| 综合色av麻豆| 成年版毛片免费区| 国产老妇伦熟女老妇高清| 国产 一区精品| 国产日韩欧美亚洲二区| 亚洲精品国产成人久久av| 久久久久久久久久久免费av| 日本一二三区视频观看| 国产免费福利视频在线观看| 亚洲av免费高清在线观看| 欧美xxⅹ黑人| 狂野欧美白嫩少妇大欣赏| 欧美日韩视频精品一区| 晚上一个人看的免费电影| 久久久久久伊人网av| 一级毛片aaaaaa免费看小| 深爱激情五月婷婷| 成人免费观看视频高清| 毛片女人毛片| 亚洲成人中文字幕在线播放| 成人午夜精彩视频在线观看| 插阴视频在线观看视频| 日本色播在线视频| 午夜视频国产福利| 丰满少妇做爰视频| 97人妻精品一区二区三区麻豆| 成人午夜精彩视频在线观看| 日韩免费高清中文字幕av| 日韩欧美精品v在线| 精品一区二区三卡| 久久午夜福利片| 久久国产乱子免费精品| 国产探花在线观看一区二区|