A VARIATIONAL ANALYSIS OF EINSTEIN SCALAR FIELD LICHNEROWICZ EQUATIONS ON COMPACT RIEMANNIAN MANIFOLDS

Size: px
Start display at page:

Download "A VARIATIONAL ANALYSIS OF EINSTEIN SCALAR FIELD LICHNEROWICZ EQUATIONS ON COMPACT RIEMANNIAN MANIFOLDS"

Transcription

1 A VARIATIONAL ANALYSIS OF EINSTEIN SCALAR FIELD LICHNEROWICZ EQUATIONS ON COPACT RIEANNIAN ANIFOLDS EANUEL HEBEY, FRANK PACARD, AND DANIEL POLLACK Abstract. We establis new existence and non-existence results for positive solutions of te Einstein scalar field Licnerowicz equation on compact manifolds. Tis equation arises from te Hamiltonian constraint equation for te Einstein scalar field system in general relativity. Our analysis introduces variational tecniques, in te form of te mountain pass lemma, to te analysis of te Hamiltonian constraint equation, wic as been previously studied by oter metods. 1. Introduction One of te foundations in te matematical analysis of te Einstein field equations of general relativity is te rigorous formulation of te Caucy problem. Te basic local existence result of Coquet-Bruat [10], and te important extension of tis due to Coquet-Bruat and Geroc [5], allows one to approac te study of globally yperbolic spacetimes via te analysis of initial data sets. Te Gauss and Codazzi equations impose constraints on te coices of initial data in general relativity, and tese constraints are expressed by te Einstein constraint equations. Tis perspective, originally studied in te context of vacuum spacetimes, as also been successfully employed in te study of many non-vacuum models obtained by minimally coupling gravity to many of te classical matter and field sources, suc as electromagnetism via te axwell equations), Yang-ills fields, fluids, and oters [8, 11, 12]. One of te simplest non-vacuum systems is te Einstein-scalar field system wic arises in coupling gravity to a scalar field satisfying a linear or nonlinear wave equation wit respect to te Lorentz metric describing te gravitational field. Te Einstein-scalar field system, wen posed in tis generality, includes as special cases te massless or massive) Einstein-Klein-Gordon equations as well as te vacuum Einstein equations wit a positive or negative) cosmological constant. Einstein scalar field teories ave been te subject of interesting developments in recent years. Among tese are te recent attempts to use suc teories to explain te observed acceleration of te expansion of te universe [15, 16, 17, 18]. Using te conformal metod, Coquet-Bruat, Isenberg, and Pollack [6, 7] reformulated te constraint equations for te Einstien scalar field system as a determined system of nonlinear partial differential equations. Te equations are semi-decoupled in te constant mean curvature CC) setting. One of tese equations, te conformally formulated momentum constraint, is a linear elliptic equation and its solvability is easy to address. Te oter one, te conformally formulated Hamiltonian constraint, is a nonlinear elliptic equation te Einstein scalar field Licnerowicz equation) as Date: February 2,

2 2 EANUEL HEBEY, FRANK PACARD, AND DANIEL POLLACK in 1.1) below see [3] for a survey on te constraint equations, and in particular, te conformal metod). Tis nonlinear equation, wic contains bot a positive critical Sobolev nonlinearity and a negative power nonlinearity, turns out to be of great matematical interest. In tis paper we provide a variational analysis of tis equation under certain conditions on its coefficients. Te analysis of te Licnerowicz equations wic arise as te conformally formulated Hamiltonian constraint equations in bot vacuum and non-vacuum settings as, in te past, been conducted primarily by eiter te metod of sub and supersolutions i.e. a barrier metod) or by perturbation or fixed point metods. Tis approac as been sufficient to allow for a complete understanding of solvability in, for example, te case of constant mean curvature vacuum initial data on compact manifolds [11]. In [7] tis metod was applied to constant mean curvature initial data for te Einstein-scalar field system on compact manifolds. In a number of cases, te metod of sub and supersolutions was sown to be sufficient to completely analyze te solvability of te Einstein scalar field Licnerowicz equation. In oter cases, te limitations of tis metod were exposed and only partial results were obtained. We establis ere two general teorems concerning non-existence and existence respectively, of positive solutions to te Einstein scalar field Licnerowicz equation 1.1). Tese results are of interest due bot to teir application to questions of existence and non-existence of solutions of te Einstein scalar field constraint equations, as well as, more generally, te introduction of variational tecniques to te analysis of te constraint equations. We expect tat similar variational tecniques will be of use in resolving oter open questions concerning initial data for te Caucy problem in general relativity. In wat follows we let, g) be a smoot compact Riemannian manifold of dimension n 3. We let also H 1 ) be te Sobolev space of functions in L 2 ) wit one derivative in L 2 ). Te H 1 norm on H 1 ) is given by u H 1 = u 2 + u 2 ) dv g. Let 2 = 2n n 2, so tat 2 is te critical Sobolev exponent for te embedding of H 1 into Lebesgue s spaces. Let also, A, and B be smoot functions on. We consider te following Einstein scalar field Licnerowicz type equations g u + u = Bu A u 2 +1, 1.1) were g = div g is te Laplace-Beltrami operator, and u > 0. Unless oterwise stated, solutions are always required to be smoot and positive. Te relationsip between te coefficients in 1.1) and initial data for te Einstein scalar field system are as follows see [7] for more details). We first note tat te sign convention for te Laplace-Beltrami operator wic we use ere is te opposite of te one used in [7]. Te conformal initial data for te purely gravitational portion of te Einstein scalar field system consists of a background Riemannian metric g indicating a coice of conformal class for te pysical metric) togeter wit a symmetric 0, 2)-tensor σ wic is divergence-free and trace-free wit respect to g so tat σ is wat is commonly referred to as a transverse-traceless, or TT-tensor) and a scalar function τ representing te mean curvature of te Caucy surface in te spacetime development of te initial data set. Te initial data for te scalar field

3 EINSTEIN SCALAR FIELD LICHNEROWICZ EQUATIONS 3 consists of two functions, ψ and π on, representing respectively te initial value for te scalar field and its normalized time derivative. Wit respect to tis set of conformal initial data, te constraint equations for te Einstein scalar field system can be realized as a determined elliptic system wose unknowns consist of a positive scalar function φ and a vector field W on. As previously remarked, in te CC case wen τ is constant) tis system becomes semi-decoupled. Tis means tat te portion of it corresponding to te momentum constraint equation is a linear, elliptic, vector equation for W in wic te unknown φ does not appear. Tis equation as a unique solution wen, g) as no conformal Killing vector fields. Te solution, W, of tis conformally formulated momentum constraint equation ten appears in te one of te coefficients of te conformally formulated Hamiltonian constraint equation wic is wat we refer to as te Einstein scalar field Licnerowicz equation. A positive solution φ of te Einstein scalar field Licnerowicz equation is ten used wit te vector field W to transform te conformal initial data set g, σ, τ, ψ, π) into a pysical initial data set satisfying te Einstein scalar field constraint equations see [7]). In terms of te conformal initial data set and te vector field W satisfying te conformally formulated momentum constraint equation) te coefficients of te Einstein scalar field Licnerowicz equation 1.1) are = c n Rg) ψ 2 g ), A = cn σ + DW 2 g + π 2) and ) n 1 B = c n n τ 2 4V ψ) were c n = n 2 4n 1), Rg) is te scalar curvature, is te covariant derivative for g, V ) is te potential in te wave equation for te scalar field, and te operator D is te conformal Killing operator relative to g, defined by DW ) ab := a W b + b W a 2 n g ab m W m. Te kernel of D consists of te conformal Killing fields on, g). Note tat relative to te notation of [7], we ave = R g,ψ, B = B τ,ψ and A = A g,w,π. We assume in wat follows tat A 0 in. Tis assumption implies no pysical restrictions since we always ave tat A 0 in te original Einstein scalar field teory. One of te results of [7] is te definition of a conformal invariant, te Yamabe scalar field conformal invariant, wose sign can be used, troug a judicious coice of te background metric g, to control te sign of. We prove two type of results in tis paper. Te first one, in Section 2, establises a set of sufficient conditions to guarantee te nonexistence of positive solutions of 1.1). Te second one, in Section 3, is concerned wit te existence of positive solutions of 1.1). Our existence result corresponds to but generalizes) te case of initial data wit a positive Yamabe scalar field conformal invariant considered in [7]. ore specifically te results presented ere sould be contrasted wit te partial results indicated in te tird row of Table 2 of [7], and specifically wit Teorems 4 and 5 in of [7]. Te results presented ere apply, for example, wen considering initial data for te Einstein massive Klein Gordon system wit small relative to te mass), or zero, values of te mean curvature. Te basic variational metod employed ere is to use te mountain pass lemma [1] to solve a family of ε-approximated equations, and let ten ε 0 to obtain a solution of

4 4 EANUEL HEBEY, FRANK PACARD, AND DANIEL POLLACK 1.1). Finally, Section 4 contains a brief discussion of a class of sligtly more general equations wic arise wen considering te Einstein axwell scalar field teory. 2. Nonnexistence of smoot positive solutions Examples of nonexistence results involving pointwise conditions on, A, and B are easy to get. Let u be a smoot positive solution of 1.1), and x 0 be a point were u is minimum. Ten g ux 0 ) 0 and we get tat x 0 )ux 0 ) Bx 0 )ux 0 ) Ax 0 )ux 0 ) 2 1. Let us assume tat bot A and B are positive functions. We ave x 0 ) Bx 0 ) X + Ax 0 ) X 1 n, 2.1) were we ave set X = ux 0 ) 4 n 2. Studying te least value of te rigt and side of 2.1) considered as a function of X), we get tat 1.1) does not possess a smoot positive solution if n n + ) n ) > max n 1) n 1 A B n ) It also follows from 2.1) tat ) n 2 4n 1) A ux) ux 0 ) min + for all x. Te idea of getting suc a bound will be used again in Section 3 wen proving Teorem 3.1. We now obtain a nonexistence result involving te Lebesgue norm of te functions A, B and. Teorem 2.1. Let, g) be a smoot compact Riemannian manifold of dimension n 3. Let also, A, and B be smoot functions on wit A 0 in. If B > 0 in, and n n n 1) n 1 ) n+2 4n A n+2 3n 2 4n B 4n dvg > + ) n+2 4 B 2 n 4 dv g, 2.3) were + = max0, ), ten te Einstein scalar field Licnerowicz equation 1.1) does not possess any smoot positive solution. Proof. We assume first tat B > 0. Let u be a smoot positive solution of 1.1). Integrating 1.1) over we get tat Bu 2 1 Adv g dv g + u 2 +1 = udv g. 2.4) By Hölder s inequality, ) 4 ) n 2 udv g + ) n+2 4 B 2 n n+2 4 dv g B u 2 1 n+2 dv g. Again by using Hölder s inequality, ) 3n 2 ) n+2 A n+2 3n 2 4n B 4n dvg B u 2 1 4n 4n Adv g dv g u Collecting tese inequalities and using 2.4), we get ) 4n ) 4 X + A n+2 3n 2 n+2 4n B 4n dvg X 1 n + ) n+2 4 B 2 n n+2 4 dv g 2.5)

5 EINSTEIN SCALAR FIELD LICHNEROWICZ EQUATIONS 5 were we ave set X = B u 2 1 dv g ) 4 n+2 Te study of te minimal value of te function of X wic appears on te left and side of 2.5) implies tat n n ) 4n ) 4n n 1) n 1 A n+2 3n 2 n+2 4n B 4n dvg + ) n+2 4 B 2 n n+2 4 dv g Tis completes te proof of te teorem. any more restrictive nonexistence conditions can be obtained easily from 2.3). For example, replacing B by min B in te two integrals in 2.3), we get tat if n n ) n+2 4n n 1) n 1 A n+2 4n dvg > + ) n+2 4 dv g min B) n 1)n+2) 4n is fulfilled, ten 2.3) olds true and te Einstein scalar field Licnerowicz equation 1.1) does not possess any smoot positive solution. In te same spirit, note tat condition 2.2) is more restricitive tan 2.3) since, for any triple of functions satisfying 2.2) we ave raising tis to te power n+2 4n n n 3n 2 A B n+2 n 1) n 1 > + ) n B n2 n) n+2 and integrating te result over yields 2.3). In wat follows we let S = S, g), S > 0, be te Sobolev constant of, g) defined as te smallest S > 0 suc tat u 2 dv g S u 2 + u 2) ) 2 2 dv g 2.6) for all u H 1 ). Explicit upper bounds for S can be given in special geometries, like, see Ilias [13], wen te Ricci curvature of te manifold is positive. Concerning lower bounds, it is well-known tat S Kn 2, were K n is te sarp Sobolev constant in te n-dimensional Euclidean space for te Sobolev inequality u L 2 K n u L 2. By letting u = 1 in 2.6) we also get tat S V 2 /n g, were V g is te volume of wit respect to g. Using tis, we prove some nonexistence result for solutions wit bound an a priori bound on teir H 1 energy. Teorem 2.2. Let, g) be a smoot compact Riemannian manifold of dimension n 3. Let also, A, and B be smoot functions on wit A 0 in. If B is arbitrary, not necessarily positive, and A 1 2 dvg > SΛ max 2 B + max 1, max + ) SΛ 4 n 2 ) ) for some Λ > 0, were B = max0, B) and S is as in 2.6), ten te Einstein scalar field Licnerowicz equation 1.1) does not possess smoot positive solutions of energy u H 1 Λ. oreover, 2.7) is sarp in te sense tat te power p = 1 2 in te left and side of 2.7) cannot be improved, and tat te bound on te energy cannot be removed.

6 6 EANUEL HEBEY, FRANK PACARD, AND DANIEL POLLACK Proof. We prove ere tat 2.7) proibits te existence of positive solutions of 1.1). Te discussion on te sarpness of tis condition is postponed after te proof. Let u be a smoot positive solution of 1.1) suc tat u H 1 Λ, Λ > 0. Let C = max 1, max + ), were + = max0, ). Ten, u 2 + u 2) dv g C u 2 + u 2) dv g. 2.8) ultiplying 1.1) by u, and integrating over, we get by 2.8) tat Adv Bu 2 g dv g + C Λ ) u 2 By te Sobolev inequality 2.6) we can write tat Bu 2 dv g max B ) SΛ 2, 2.10) were B = max 0, B). Ten, by combining 2.9) 2.10) we get tat Adv g C Λ 2 + max u B ) SΛ ) 2 Now, Hölder inequality yields A 1 2 dvg ) 1 Adv ) g u 2 dv g. 2.12) u 2 By combining tis inequality wit 2.11), and by te Sobolev inequality 2.6), we get tat A 1 2 dvg SΛ max 2 B + C ) 1 2. SΛ 4 n 2 Tis proves te teorem. We now discuss te sarpness of 2.7) in Teorem 2.2. Te Yamabe equation on a Riemannian manifold, g) may be written as g u + n 2 Rg)u = u2 1, 2.13) 4n 1) were Rg) is te scalar curvature of g. A positive solution u > 0 of 2.13) corresponds to a conformally related metric g = u 2 2 g wit constant positive scalar curvature R g) = 4n 1) n 2. Now, any solution of 2.13) is a solution of 1.1) wen we let = n 2 4n 1) Rg), B = α, and A = 1 α)u22 for some α R. Tis provides a transformation rule for rewriting equations like 2.13) into equations like 1.1). On te unit spere S n, g), for wic Rg) = nn 1), we know see, for instance, Aubin [2]) tat tere exist families u ε ) ε of solutions of 2.13), ε > 0, suc tat u ε H 1 = Kn n + o1) for all ε > 0, and u ε L p + as ε 0 for all p > 2, were K n is te sarp Sobolev constant in te n-dimensional Euclidean space for te Sobolev inequality u L 2 K n u L 2. Letting α = 1 2, te above transformation rule 2.13) 1.1) provides a family of Einstein scalar field Licnerowicz type equations indexed by ε > 0, wit and B independent of ε, suc tat any equation in te family possesses a solution of energy less tan or equal to 2Kn n, and for wic Ap εdv g + as ε 0 for all p > 1 2. Tis proves tat te power p = 1 2 in te left and side of 2.7) cannot be improved. Tis example can be modified in

7 EINSTEIN SCALAR FIELD LICHNEROWICZ EQUATIONS 7 different ways wit te constructions given in Brendle [4] and in Druet and Hebey [9]. We prove next tat te bound on te energy in Teorem 2.2 cannot be removed. By Druet and Hebey [9] we know tat on te unit spere in dimension n 6, or on any quotient, g) of te unit spere in dimension n 6, tere exist families ε ) ε of smoot functions, suc tat ε nn 2) 4 in C 1 ), and families u ε ) ε of smoot positive functions suc tat, for any ε > 0, u ε solves te Yamabe type equation g0 u ε + ε u ε = u 2 1 ε, 2.14) and suc tat u ε H 1 + as ε 0. Rewriting 2.14) wit te transformation rule 2.13) 1.1), we see tat te u ε s solve 1.1) wit = ε, B = α, and A = 1 α)u 22 ε for some α R. Letting α = 1 2, we get families of Einstein scalar field Licnerowicz type equations indexed by ε > 0 suc tat any equation in te family possesses a solution, B is independent of ε, te ε s converge in te C 1 -topology to a positive constant function, and A1/2 ε dv g + as ε 0. In particular, we cannot ope to get tat tere exists C = Cn,, B), depending on te manifold and continuously on and B in te C 0 -topology, like tis is te case for te constant in 2.7) wen Λ is fixed, suc tat if A1/2 dv g C, ten te Einstein scalar field Licnerowicz type equation 1.1) does not possess a smoot positive solution. Tis proves tat te bound on te energy in Teorem 2.2 cannot be removed. In te same circle of ideas, we mention tat if B > 0 in, ten we can give anoter form to 2.7) were te constant appears as CΛ 2. In order to get tis dependancy in Λ 2 we may proceed as in te proof of Teorem 2.2, but now getting bounds from te estimate 2.9). By 2.9), since we assumed tat B > 0 in, we can write tat u 2 dv g C Λ 2 min B and Adv g C Λ ) u 2 Ten, by 2.12) as in te proof of te second part of Teorem 2.2, we get from 2.15) tat 1.1) does not possess a smoot positive solution if A 1 2 dvg > max 1, max + ) Λ ) min B) 1 2 Condition 2.16) is complementary to te condition in Teorem 2.2. For large Λ s, 2.16) is better tan 2.7) since it involves te energy Λ 2 and not Λ 2n 1)/n 2). 3. Existence of a smoot positive solution In tis section we use te mountain pass lemma [1], to get existence results tat complement te nonexistence results presented in Teorem 2.2. ore precisely, we prove tat if Adv g is sufficiently small, and A > 0 in, ten 1.1) possesses a solution. Wen A 0, 1.1) is te prescribed scalar curvature equation and we know from Kazdan and Warner [14] tat tere are situations in wic te equation does not possess a solution.

8 8 EANUEL HEBEY, FRANK PACARD, AND DANIEL POLLACK In te sequel we assume tat te function is cosen so tat g + is coercive. Tis amounts to say tat tere exists a constant K = K, g, ) > 0, suc tat u 2 dv g K u 2 + u 2) dv g for all u H 1 ). It will be convenient to define u H 1 = u 2 + u 2) ) 1 2 dv g. 3.1) We also denote by S = S, g, ) > 0, te Sobolev constant defined to be te smallest constant S > 0 suc tat u 2 dv g S u 2 + u 2) ) 2 2 dv g 3.2) for all u H 1 ). Observe tat, if > 0 in, ten g + is coercive and conversely coercivity implies tat dv g > 0, and tus tat max > 0. Also observe tat if A, B 0, A + B > 0, and if 1.1) possesses a smoot positive solution, ten g + is coercive. Indeed, in tat case, tere exists a function u > 0 suc tat g u+u > 0 everywere in, and te existence of suc an u implies te coercivity of g +. Finally, as already mentioned, wen > 0 in, ten g + is coercive and we ave te bound ) S max 1, S. min were S = S, g) > 0 is te Sobolev constant defined in 2.6). We prove ere tat te following existence result olds true. Teorem 3.1. Let, g) be a smoot compact Riemannian manifold of dimension n 3. Let, A, and B be smoot functions on for wic g + is coercive, A > 0 in, and max B > 0. Tere exists a constant C = Cn), C > 0 depending only on n, suc tat if and ϕ 2 H 1 A ϕ 2 dv g C S max B ) n 1 3.3) Bϕ 2 dv g > 0 for some smoot positive function ϕ > 0 in, were H 1 is as in 3.1) and S is as in 3.2), ten te Einstein scalar field Licnerowicz equation 1.1) possesses a smoot positive solution. Proof. Preliminary computations We define I 1) : H 1 ) R by I 1) u) = 1 u 2 + u 2) dv g Bu + ) 2 dv g, 3.4) and if we fix ε > 0 we define I 2) ε : H 1 ) R by I ε 2) u) = 1 2 Adv g, 3.5) ε + u + ) 2 2 )

9 EINSTEIN SCALAR FIELD LICHNEROWICZ EQUATIONS 9 were 2 = 2 2. Obviously, for any u H 1 ) we can write Φ u H 1 ) I 1) u) Ψ u H 1 ) 3.6) if te functions Φ, Ψ : [0, + ) R are defined by and Φt) = 1 2 t2 max B 2 S t 2 Ψt) = 1 2 t2 + max B 2 S t 2 for t R, were S > 0 and H 1 are as in 3.1) and 3.2). Let t 0 > 0 be given by t 0 = 1 S max B ) n ) 3.8) 3.9) so tat Φ is increasing in [0, t 0 ], and decreasing in [t 0, + ). We define θ > 0 suc tat θ 2 1 = 2n 1) and t 1 = θ t 0 for t 0 as in 3.9). It is easy to ceck tat Ψt 1 ) θ Φt 0) 1 2 Φt 0), 3.10) were Φ and Ψ are as in 3.7) and 3.8). Finally, we define te functional I ε = I 1) + I 2) ε, 3.11) were I 1) and I ε 2) are as in 3.4) and 3.5). Let ϕ C ), ϕ > 0 in, be te function in te statement of te teorem. In particular Bϕ 2 dv g > ) and, witout loss of generality, we can assume tat ϕ H 1 = 1. Now, provided te constant C in 3.3) is cosen to be C = θ 2 2 2, 4 we find tat 3.3) precisely translates into 1 A 2 dv g 1 t 1 ϕ) 2 2 Φt 0) 3.13) and by 3.6), 3.10), and 3.13) we get tat Finally, 3.12) implies tat I ε t 1 ϕ) Φt 0 ) < I ε t 0 ϕ) 3.14) lim I εt ϕ) =. +

10 10 EANUEL HEBEY, FRANK PACARD, AND DANIEL POLLACK Hence we can coose t 2 > t 0 suc tat were I ε is te functional in 3.11). I ε t 2 ϕ) < 0, 3.15) Application of te ountain Pass Lemma By 3.14) and 3.15), we can apply te mountain pass lemma [1] to te functional I ε. Let c ε = inf γ Γ max u γ I εu), 3.16) were Γ stands for te set of continuous pats joining u 1 = t 1 ϕ to u 2 = t 2 ϕ. Observe tat c ε > Φt 0 ) and, taking te pat γt) = t ϕ, for t [t 1, t 2 ], we see tat c ε is bounded uniformly as ε tends to 0. We will keep in mind, for furter use tat for all ε small enoug, were c > 0 is independent of ε. Φt 0 ) < c ε c 3.17) By te mountain pass lemma we get tat tere exists a sequence u k ) k in H 1 ) suc tat I ε u k ) c and I εu k ) ) as k +. By 3.18), = u k ϕ)dv g + u k ϕdv g Bu + 1 k ϕdv )2 g Au + k ϕdv ) 3.19) g ε + u + k )2 ) + o ϕ 2 +1 H 1 for all ϕ H 1 ), were u k ϕ) stands for te pointwise scalar product of u k and ϕ wit respect to g, and 1 uk 2 + u 2 ) k dvg Bu + k dv )2 g ) Adv g 2 ε + u + k )2 ) = c ε + o1). 2 Combining 3.19) wit ϕ = u k, and 3.20), we get tat 1 Bu + k dv n )2 g + 1 Au + k )2 dv g 2 ε + u + k )2 ) Adv ) g 2 ε + u + k )2 ) = c ε + o u k H + o1), ) and it follows from 3.21) tat for k sufficiently large, 1 ) Bu + k dv n )2 g 2c ε + o u k H. 3.22) 1 By 3.20) and 3.22) we ten get tat for k sufficiently large uk 2 + u 2 k) dvg n 2 Bu + k dv n )2 g + 4 c ε ) 2n c ε + o u k H )

11 EINSTEIN SCALAR FIELD LICHNEROWICZ EQUATIONS 11 In particular, by 3.22) and 3.23), uk 2 + u 2 ) k dvg 2n c ε + 1, and 4n n 2 c ε Bu + k )2 dv g 3nc ε 3.24) for k sufficiently large, were c ε is as in 3.16). By 3.24), te sequence u k ) k is bounded in H 1 ). Up to passing to a subsequence we may ten assume tat tere exists u ε H 1 ) suc tat u k u ε weakly in H 1 ), u k u ε strongly in L p ) for some p > 2, and u k u ε almost everywere in as k +. As a consequence, u + k )2 1 u + ε ) 2 1 weakly in L 2 /2 1) ), and u + k ε + u + k )2 ) q u + ε ε + u + ε ) 2 ) q strongly in L2 ) 3.25) for all q > 0, as k +. Indeed, by 3.24), te u + k 1 s are bounded in )2 L 2 /2 1) ). Since tey converge almost everywere to u + ε ) 2 1, te first equation in 3.25) follows from standard integration teory. By te Lebesgue s dominated convergence teorem we also ave tat ε + u + k )2 ) q ε + u + ε ) 2 ) q strongly in L p ) for all p 1 and all q > 0, and since u k u ε in L p ) for some p > 2, we easily get tat te second equation in 3.25) olds true. By 3.25), letting k + in 3.19), it follows tat u ε satisfies g u ε + u ε = Bu + ε ) Au + ε ε + u + ε ) 2 ) ) in te weak sense. Te weak maximum principle and 3.26) imply tat u ε 0. As a consequence, g u ε + u ε = Bu 2 1 Au ε ε ) ε + u 2 ε) 2 +1 in te weak sense. Regularity and positivity of te solution We may rewrite 3.27) as ) A g u ε + u ε + u 2 ε) 2 +1 ε = Bu 2 1 ε, and since A ε + u 2 ε) 2 +1 L ), te regularity arguments developed in Trudinger [19] apply to 3.27). It follows tat u ε L s ) for some s > 2. Since we ave tat Aε + u 2 ε) 2 +1 u ε L p ) if u ε L p ), and u ε L s ) for some s > 2, te standard bootstrap procedure, togeter wit regularity teory, gives tat u ε H 2,p ) for all p 1, were H 2,p is te Sobolev space of functions in L p wit two derivatives in L p. By te Sobolev embedding teorem we ten get tat te rigt and side in 3.27) is in C 0,α ) for α 0, 1), and by regularity teory it follows tat u ε C 2,α ) for α 0, 1). In particular, te strong maximum principle can be applied and we get tat eiter

12 12 EANUEL HEBEY, FRANK PACARD, AND DANIEL POLLACK u ε 0, or u ε > 0 in. Ten we easily get tat u ε C ) is smoot. By 3.24) and 3.25), letting k + in 3.21), we get tat 1 Adv g 2 ε + u 2 ε) 2 1)c. 3.28) 2 were c is te upper bound for c ε. If, for a sequence of ε j tending to 0, u εj were to be equal to 0, we would conclude tat 1 Adv g c 2 2 1)ε 2 j 3.29) wic is clearly impossible since we ave assumed tat A > 0. Terefore, for ε small enoug u ε 0. Ten, according to te above discussion, u ε is a smoot positive solution of 3.27). By 3.24), and standard properties of te weak limit, we also get tat uε 2 + u 2 ) ε dvg 2nc ε ) for all ε > 0 small enoug. Passing to te limit as ε tends to 0 In wat follows we let ε k ) k be a sequence of positive real numbers suc tat ε k 0 as k + and 3.29) olds true wit ε = ε k for all k, and let u k = u εk. Ten u k is a smoot positive function in suc tat g u k + u k = Bu 2 1 k + Au k ε k + u 2 k ) ) in wile, by 3.17) and 3.30), te sequence u k ) k is bounded in H 1 ). Let x k be a point were u k is minimum. Ten g u k x k ) 0 and we get wit 3.31) tat Let δ 0 > 0 be suc tat x k ) + B x k )u k x k ) 2 2 δ 22 +1) 0 max Ax k ) ε k + u k x k ) 2 ) ) ) + max 2 B )δ2 0 = min A 2 By 3.32) we obtain tat u k x k ) δ 0, and tus tat min u k δ ) wen k is sufficiently large. Since u k ) k is bounded in H 1 ) we may assume tat tere exists u H 1 ) suc tat, up to passing to a subsequence, u k u weakly in H 1 ), u k u strongly in L p ) for some p > 2, and u k u almost everywere in as k +. By 3.33), u δ 0 almost everywere in. Still by 3.33), we get wit similar arguments to tose used to prove 3.25) tat u 2 1 k u 2 1 weakly in L 2 /2 1) ), and u k ε k + u 1 2 k )2 +1 u 2 +1 strongly in L2 ). 3.34) as k +. By 3.31) and 3.34), letting k + in 3.31), we get tat u is a weak solution of te Einstein scalar field Licnerowicz equation 1.1). Rewriting 1.1) as g u + A ) u = Bu 2 1, u 2 +2

13 EINSTEIN SCALAR FIELD LICHNEROWICZ EQUATIONS 13 and since Au 2 2 L ), te regularity arguments developed in Trudinger [19] apply to 1.1). It follows tat u L s ) for some s > 2. Since u δ 0 almost everywere, and δ 0 > 0, te standard bootstrap procedure, togeter wit regularity teory, gives tat u is a smoot positive solution of 1.1). Tis ends te proof of te teorem. As a remark, te above proof provides an explicit expression for te dimensional constant C in 3.3). As anoter remark, it can be noted tat wen Bdv g > 0, ten we can take ϕ to be constant in 3.12). In particular, our existence result as te following Corollary. Corollary 3.1. Let, g) be a smoot compact Riemannian manifold of dimension n 3 and a smoot functions on for wic g + is coercive. Tere exists a constant C = Cn, ), C > 0, suc tat if A and B are smoot functions on, wit A > 0 in, max B > 0, and Bdv g > 0, and if we furter assume tat max B )n 1 A dv g Cn, ), 3.35) ten te Einstein scalar field Licnerowicz equation 1.1) possesses a smoot positive solution. Wen A > 0 and B > 0, we can also take ϕ = A n 2 4n result as te following Corollary. in 3.12), and our existence Corollary 3.2. Let, g) be a smoot compact Riemannian manifold of dimension n 3 and a smoot functions on for wic g + is coercive. Tere exists a constant C = Cn, ), C > 0, suc tat if A and B are smoot functions on, wit A > 0 and B > 0 in and if we furter assume tat max B )n 1 A n 2 4n 2 H 1 A 1 2 dvg Cn, ), 3.36) ten te Einstein scalar field Licnerowicz equation 1.1) possesses a smoot positive solution. Interestingly, Sobolev embedding implies tat A 1 2 dvg S A n 2 4n 2 H 1 and so, if A and B satisfy 3.36), ten max B )n 1 A 1 2 dvg ) 2 Cn, ) S wic is reminiscent of condition tat ensured te non existence of a solution wic was obtained in Teorem Einstein-axwell-scalar field teory Te metods employed in sections 2 and 3 are strong enoug to deal wit additional nonlinear negative power terms in te equation of te form Cu p for C 0 and p > 1. Suc terms arise, for example, in te Einstein-axwell-scalar field teory. Given, g) compact of dimension n 3, we let, A, B, and C be smoot functions in, and we briefly discuss in tis section equations of te form g u + u = Bu A u C u p, 4.1),

14 14 EANUEL HEBEY, FRANK PACARD, AND DANIEL POLLACK were A, C 0 and p > 1. In te case of te Einstein-axwell-scalar field teory in spatial) dimension n = 3 we ave p = 3 and C 0 represents te sum of te squares of te norms of te electric and magnetic fields on. Te approac we used to prove Teorem 2.2 deals wit inequalities resulting from te signs of te coefficients and te powers of te unknown function u and tus applies to 4.1). Let ˆp = 2 +p Ten, if we concentrate on getting nonexistence results of smoot positive solutions wit no a priori bound on te energy, te approac we used to prove Teorem 2.2 gives in particular tat 4.1) does not possess a smoot positive solution if B > 0 in, A, C 0 in, and eiter 2.3) olds true, or α + 1) α+1 α α ) 1 ˆp C 1ˆp B ˆp 1 ˆp dv g > + ) n+2 4 B 2 n 4 dv g, 4.2) were α = n 2)p + 1)/4. We also do get similar conditions to 4.2) for te nonexistence of solutions of 4.1) of energy bounded by Λ. Te metod we used to prove Teorem 3.1 applies to 4.1) as well. Assume g + is coercive, A, C 0 in, A + C > 0 in, and max B > 0. Following te proof of Teorem 3.1 we get tat tat tere exists Λ = Λn, p), Λ > 0 depending only on n and p, suc tat if and A ϕ 2 dv g Λ S max B ) n 1, Bϕ 2 dv g > 0 C ϕ p 1 dv Λ g S max B ) α 4.3) for some smoot positive function ϕ > 0 in suc tat ϕ H 1 = 1, were H 1 is as in 3.1), S is as in 3.2), and α is as in 4.2), ten 4.1) possesses a smoot positive solution. As for 3.3), te constant Λ in 4.3) can be made explicit. References [1] A. Ambrosetti and P. Rabinowitz, Dual variational metods in critical point teory and applications, J. Functional Analysis ), [2] T. Aubin, Nonlinear analysis on manifolds. onge-ampre equations. Grund. der at. Wissenscaften, 252 Springer-Verlag, New York, [3] R. Bartnik and J. Isenberg, Te Constraint Equations, in Te Einstein equations and te large scale beavior of gravitational fields P.T. Cruściel and H. Friedric, eds.), Birkäuser, Basel, 1 39, [4] S. Brendle, Blow-up penomena for te Yamabe PDE in ig dimensions, Preprint 2006). [5] Y. Coquet-Bruat and R. Geroc, Global aspects of te Caucy problem in general relativity, Comm. at. Pys ), [6] Y. Coquet-Bruat, J. Isenberg and D. Pollack Te Einstein scalar field constraints on asymptotically Euclidean manifolds, Cin. Ann. at., Vol 27, ser. B, no ), 31 52; ArXiv: gr-qc/ [7] Coquet-Bruat, Y., Isenberg, J., and Pollack, D., Te constraint equations for te Einstein scalar field system on compact manifolds, Class. Quantum Grav ) : ArXiv: gr-qc/ [8] Y. Coquet-Bruat and J. York, Te Caucy Problem, in General Relativity and Gravitation - Te Einstein Centenary, A. Held ed.), Plenum, , [9] O. Druet and E. Hebey, Blow-up examples for second order elliptic PDEs of critical Sobolev growt, Trans. Amer. at. Soc ) [10] Y. Foures-Bruat, Téorème d existence pour certains systèmes d équations aux dérivées partialles non linéaires, Acta. at ) [11] J. Isenberg, Constant mean curvature solutions of te Einstein constraint equations on closed manifolds, Class. Quantum Grav )

15 EINSTEIN SCALAR FIELD LICHNEROWICZ EQUATIONS 15 [12] J. Isenberg, D. axwell and D. Pollack, A gluing constructions for non-vacuum solutions of te Einstein constraint equations, Adv. Teor. at. Pys ), no. 1, ; ArXiv: gr-qc/ [13] S. Ilias, Constantes explicites pour les inégalités de Sobolev sur les variétés riemanniennes compactes, Ann. Inst. Fourier ) [14] J.L. Kazdan and F.W. Warner, Scalar curvature and conformal deformation of Riemannian structure, J. Differential Geometry ) [15] A. Rendall, Accelerated cosmological expansion due to a scalar field wose potential as a p ositive lower bound, Class. Quantum Grav. 21, ); ArXive: gr-qc/ [16] A. Rendall, atematical properties of cosmological models wit accelerated expansion., Analytical and numerical approaces to matematical relativity, , Lecture Notes in Pys., 692, Springer, Berlin, 2006; ArXive: gr-qc/ [17] A. Rendall, Intermediate inflation and te slow-roll approximation, Class. Quantum Grav. 22, ); ArXive: gr-qc/ [18] V. Sani, Dark matter and dark energy, In E. Papantonopoulos ed.) Pysics of te early universe. Springer, Berlin. 2005): ArXive: astro-p/ [19] N.S. Trudinger, Remarks concerning te conformal deformation of Riemannian structures on compact manifolds, Ann. Scuola Norm. Sup. Pisa ) Emmanuel Hebey, Université de Cergy-Pontoise, Département de atématiques, Site de Saint-artin, 2 avenue Adolpe Cauvin, Cergy-Pontoise cedex, France address: Emmanuel.Hebey@mat.u-cergy.fr Frank Pacard, Université Paris XII, Département de atématiques, 61 avenue du Général de Gaulle, Créteil cedex, France address: pacard@univ-paris12.fr Daniel Pollack, University of Wasington, Department of atematics, Box , Seattle, WA , USA address: pollack@mat.wasington.edu

Existence, stability and instability for Einstein-scalar field Lichnerowicz equations by Emmanuel Hebey

Existence, stability and instability for Einstein-scalar field Lichnerowicz equations by Emmanuel Hebey Existence, stability and instability for Einstein-scalar field Lichnerowicz equations by Emmanuel Hebey Joint works with Olivier Druet and with Frank Pacard and Dan Pollack Two hours lectures IAS, October

More information

Poisson Equation in Sobolev Spaces

Poisson Equation in Sobolev Spaces Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on

More information

Differentiation in higher dimensions

Differentiation in higher dimensions Capter 2 Differentiation in iger dimensions 2.1 Te Total Derivative Recall tat if f : R R is a 1-variable function, and a R, we say tat f is differentiable at x = a if and only if te ratio f(a+) f(a) tends

More information

Stationary Kirchhoff equations with powers by Emmanuel Hebey (Université de Cergy-Pontoise)

Stationary Kirchhoff equations with powers by Emmanuel Hebey (Université de Cergy-Pontoise) Stationary Kirchhoff equations with powers by Emmanuel Hebey (Université de Cergy-Pontoise) Lectures at the Riemann center at Varese, at the SNS Pise, at Paris 13 and at the university of Nice. June 2017

More information

RIGIDITY OF TIME-FLAT SURFACES IN THE MINKOWSKI SPACETIME

RIGIDITY OF TIME-FLAT SURFACES IN THE MINKOWSKI SPACETIME RIGIDITY OF TIME-FLAT SURFACES IN THE MINKOWSKI SPACETIME PO-NING CHEN, MU-TAO WANG, AND YE-KAI WANG Abstract. A time-flat condition on spacelike 2-surfaces in spacetime is considered ere. Tis condition

More information

4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.

4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these. Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra

More information

232 Calculus and Structures

232 Calculus and Structures 3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE

More information

Symmetry Labeling of Molecular Energies

Symmetry Labeling of Molecular Energies Capter 7. Symmetry Labeling of Molecular Energies Notes: Most of te material presented in tis capter is taken from Bunker and Jensen 1998, Cap. 6, and Bunker and Jensen 2005, Cap. 7. 7.1 Hamiltonian Symmetry

More information

Lecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator.

Lecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator. Lecture XVII Abstract We introduce te concept of directional derivative of a scalar function and discuss its relation wit te gradient operator. Directional derivative and gradient Te directional derivative

More information

How to Find the Derivative of a Function: Calculus 1

How to Find the Derivative of a Function: Calculus 1 Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te

More information

Mass Lumping for Constant Density Acoustics

Mass Lumping for Constant Density Acoustics Lumping 1 Mass Lumping for Constant Density Acoustics William W. Symes ABSTRACT Mass lumping provides an avenue for efficient time-stepping of time-dependent problems wit conforming finite element spatial

More information

Reflection Symmetries of q-bernoulli Polynomials

Reflection Symmetries of q-bernoulli Polynomials Journal of Nonlinear Matematical Pysics Volume 1, Supplement 1 005, 41 4 Birtday Issue Reflection Symmetries of q-bernoulli Polynomials Boris A KUPERSHMIDT Te University of Tennessee Space Institute Tullaoma,

More information

MA455 Manifolds Solutions 1 May 2008

MA455 Manifolds Solutions 1 May 2008 MA455 Manifolds Solutions 1 May 2008 1. (i) Given real numbers a < b, find a diffeomorpism (a, b) R. Solution: For example first map (a, b) to (0, π/2) and ten map (0, π/2) diffeomorpically to R using

More information

MIXED DISCONTINUOUS GALERKIN APPROXIMATION OF THE MAXWELL OPERATOR. SIAM J. Numer. Anal., Vol. 42 (2004), pp

MIXED DISCONTINUOUS GALERKIN APPROXIMATION OF THE MAXWELL OPERATOR. SIAM J. Numer. Anal., Vol. 42 (2004), pp MIXED DISCONTINUOUS GALERIN APPROXIMATION OF THE MAXWELL OPERATOR PAUL HOUSTON, ILARIA PERUGIA, AND DOMINI SCHÖTZAU SIAM J. Numer. Anal., Vol. 4 (004), pp. 434 459 Abstract. We introduce and analyze a

More information

A SHORT INTRODUCTION TO BANACH LATTICES AND

A SHORT INTRODUCTION TO BANACH LATTICES AND CHAPTER A SHORT INTRODUCTION TO BANACH LATTICES AND POSITIVE OPERATORS In tis capter we give a brief introduction to Banac lattices and positive operators. Most results of tis capter can be found, e.g.,

More information

Consider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.

Consider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx. Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions

More information

3 Parabolic Differential Equations

3 Parabolic Differential Equations 3 Parabolic Differential Equations 3.1 Classical solutions We consider existence and uniqueness results for initial-boundary value problems for te linear eat equation in Q := Ω (, T ), were Ω is a bounded

More information

CHAPTER 4. Elliptic PDEs

CHAPTER 4. Elliptic PDEs CHAPTER 4 Elliptic PDEs One of te main advantages of extending te class of solutions of a PDE from classical solutions wit continuous derivatives to weak solutions wit weak derivatives is tat it is easier

More information

INTRODUCTION AND MATHEMATICAL CONCEPTS

INTRODUCTION AND MATHEMATICAL CONCEPTS Capter 1 INTRODUCTION ND MTHEMTICL CONCEPTS PREVIEW Tis capter introduces you to te basic matematical tools for doing pysics. You will study units and converting between units, te trigonometric relationsips

More information

Volume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households

Volume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households Volume 29, Issue 3 Existence of competitive equilibrium in economies wit multi-member ouseolds Noriisa Sato Graduate Scool of Economics, Waseda University Abstract Tis paper focuses on te existence of

More information

Approximation of the Viability Kernel

Approximation of the Viability Kernel Approximation of te Viability Kernel Patrick Saint-Pierre CEREMADE, Université Paris-Daupine Place du Marécal de Lattre de Tassigny 75775 Paris cedex 16 26 october 1990 Abstract We study recursive inclusions

More information

3.1 Extreme Values of a Function

3.1 Extreme Values of a Function .1 Etreme Values of a Function Section.1 Notes Page 1 One application of te derivative is finding minimum and maimum values off a grap. In precalculus we were only able to do tis wit quadratics by find

More information

1 Calculus. 1.1 Gradients and the Derivative. Q f(x+h) f(x)

1 Calculus. 1.1 Gradients and the Derivative. Q f(x+h) f(x) Calculus. Gradients and te Derivative Q f(x+) δy P T δx R f(x) 0 x x+ Let P (x, f(x)) and Q(x+, f(x+)) denote two points on te curve of te function y = f(x) and let R denote te point of intersection of

More information

Order of Accuracy. ũ h u Ch p, (1)

Order of Accuracy. ũ h u Ch p, (1) Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical

More information

HOMEWORK HELP 2 FOR MATH 151

HOMEWORK HELP 2 FOR MATH 151 HOMEWORK HELP 2 FOR MATH 151 Here we go; te second round of omework elp. If tere are oters you would like to see, let me know! 2.4, 43 and 44 At wat points are te functions f(x) and g(x) = xf(x)continuous,

More information

COMPOSITION IN FRACTIONAL SOBOLEV SPACES

COMPOSITION IN FRACTIONAL SOBOLEV SPACES COMPOSITION IN FRACTIONAL SOBOLEV SPACES HAIM BREZIS (1)() AND PETRU MIRONESCU (3) 1. Introduction A classical result about composition in Sobolev spaces asserts tat if u W k,p (Ω) L (Ω) and Φ C k (R),

More information

COMPOSITION IN FRACTIONAL SOBOLEV SPACES. 1. Introduction. A classical result about composition in Sobolev spaces asserts

COMPOSITION IN FRACTIONAL SOBOLEV SPACES. 1. Introduction. A classical result about composition in Sobolev spaces asserts DISCRETE AND CONTINUOUS Website: ttp://mat.smsu.edu/journal DYNAMICAL SYSTEMS Volume 7, Number, April 001 pp. 41 46 COMPOSITION IN FRACTIONAL SOBOLEV SPACES HAIM BREZIS (1)() AND PETRU MIRONESCU (3) 1.

More information

Derivation Of The Schwarzschild Radius Without General Relativity

Derivation Of The Schwarzschild Radius Without General Relativity Derivation Of Te Scwarzscild Radius Witout General Relativity In tis paper I present an alternative metod of deriving te Scwarzscild radius of a black ole. Te metod uses tree of te Planck units formulas:

More information

Influence of the Stepsize on Hyers Ulam Stability of First-Order Homogeneous Linear Difference Equations

Influence of the Stepsize on Hyers Ulam Stability of First-Order Homogeneous Linear Difference Equations International Journal of Difference Equations ISSN 0973-6069, Volume 12, Number 2, pp. 281 302 (2017) ttp://campus.mst.edu/ijde Influence of te Stepsize on Hyers Ulam Stability of First-Order Homogeneous

More information

REVIEW LAB ANSWER KEY

REVIEW LAB ANSWER KEY REVIEW LAB ANSWER KEY. Witout using SN, find te derivative of eac of te following (you do not need to simplify your answers): a. f x 3x 3 5x x 6 f x 3 3x 5 x 0 b. g x 4 x x x notice te trick ere! x x g

More information

Parameter Fitted Scheme for Singularly Perturbed Delay Differential Equations

Parameter Fitted Scheme for Singularly Perturbed Delay Differential Equations International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department

More information

Chapter 5 FINITE DIFFERENCE METHOD (FDM)

Chapter 5 FINITE DIFFERENCE METHOD (FDM) MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential

More information

Gradient Descent etc.

Gradient Descent etc. 1 Gradient Descent etc EE 13: Networked estimation and control Prof Kan) I DERIVATIVE Consider f : R R x fx) Te derivative is defined as d fx) = lim dx fx + ) fx) Te cain rule states tat if d d f gx) )

More information

Weierstraß-Institut. im Forschungsverbund Berlin e.v. Preprint ISSN

Weierstraß-Institut. im Forschungsverbund Berlin e.v. Preprint ISSN Weierstraß-Institut für Angewandte Analysis und Stocastik im Forscungsverbund Berlin e.v. Preprint ISSN 0946 8633 Stability of infinite dimensional control problems wit pointwise state constraints Micael

More information

A Hybrid Mixed Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems

A Hybrid Mixed Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems A Hybrid Mixed Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems

More information

ERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*

ERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER* EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T

More information

INTRODUCTION AND MATHEMATICAL CONCEPTS

INTRODUCTION AND MATHEMATICAL CONCEPTS INTODUCTION ND MTHEMTICL CONCEPTS PEVIEW Tis capter introduces you to te basic matematical tools for doing pysics. You will study units and converting between units, te trigonometric relationsips of sine,

More information

FINITE ELEMENT EXTERIOR CALCULUS FOR PARABOLIC EVOLUTION PROBLEMS ON RIEMANNIAN HYPERSURFACES

FINITE ELEMENT EXTERIOR CALCULUS FOR PARABOLIC EVOLUTION PROBLEMS ON RIEMANNIAN HYPERSURFACES FINITE ELEMENT EXTERIOR CALCULUS FOR PARABOLIC EVOLUTION PROBLEMS ON RIEMANNIAN HYPERSURFACES MICHAEL HOLST AND CHRIS TIEE ABSTRACT. Over te last ten years, te Finite Element Exterior Calculus (FEEC) as

More information

Dynamics and Relativity

Dynamics and Relativity Dynamics and Relativity Stepen Siklos Lent term 2011 Hand-outs and examples seets, wic I will give out in lectures, are available from my web site www.damtp.cam.ac.uk/user/stcs/dynamics.tml Lecture notes,

More information

Dedicated to the 70th birthday of Professor Lin Qun

Dedicated to the 70th birthday of Professor Lin Qun Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,

More information

The derivative function

The derivative function Roberto s Notes on Differential Calculus Capter : Definition of derivative Section Te derivative function Wat you need to know already: f is at a point on its grap and ow to compute it. Wat te derivative

More information

Diffraction. S.M.Lea. Fall 1998

Diffraction. S.M.Lea. Fall 1998 Diffraction.M.Lea Fall 1998 Diffraction occurs wen EM waves approac an aperture (or an obstacle) wit dimension d > λ. We sall refer to te region containing te source of te waves as region I and te region

More information

5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems

5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems 5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we

More information

Numerical Differentiation

Numerical Differentiation Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function

More information

Decay of solutions of wave equations with memory

Decay of solutions of wave equations with memory Proceedings of te 14t International Conference on Computational and Matematical Metods in Science and Engineering, CMMSE 14 3 7July, 14. Decay of solutions of wave equations wit memory J. A. Ferreira 1,

More information

Discontinuous Galerkin Methods for Relativistic Vlasov-Maxwell System

Discontinuous Galerkin Methods for Relativistic Vlasov-Maxwell System Discontinuous Galerkin Metods for Relativistic Vlasov-Maxwell System He Yang and Fengyan Li December 1, 16 Abstract e relativistic Vlasov-Maxwell (RVM) system is a kinetic model tat describes te dynamics

More information

Superconvergence of energy-conserving discontinuous Galerkin methods for. linear hyperbolic equations. Abstract

Superconvergence of energy-conserving discontinuous Galerkin methods for. linear hyperbolic equations. Abstract Superconvergence of energy-conserving discontinuous Galerkin metods for linear yperbolic equations Yong Liu, Ci-Wang Su and Mengping Zang 3 Abstract In tis paper, we study superconvergence properties of

More information

3.4 Worksheet: Proof of the Chain Rule NAME

3.4 Worksheet: Proof of the Chain Rule NAME Mat 1170 3.4 Workseet: Proof of te Cain Rule NAME Te Cain Rule So far we are able to differentiate all types of functions. For example: polynomials, rational, root, and trigonometric functions. We are

More information

An approximation method using approximate approximations

An approximation method using approximate approximations Applicable Analysis: An International Journal Vol. 00, No. 00, September 2005, 1 13 An approximation metod using approximate approximations FRANK MÜLLER and WERNER VARNHORN, University of Kassel, Germany,

More information

Math 161 (33) - Final exam

Math 161 (33) - Final exam Name: Id #: Mat 161 (33) - Final exam Fall Quarter 2015 Wednesday December 9, 2015-10:30am to 12:30am Instructions: Prob. Points Score possible 1 25 2 25 3 25 4 25 TOTAL 75 (BEST 3) Read eac problem carefully.

More information

POLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY

POLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY APPLICATIONES MATHEMATICAE 36, (29), pp. 2 Zbigniew Ciesielski (Sopot) Ryszard Zieliński (Warszawa) POLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY Abstract. Dvoretzky

More information

THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Math 225

THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Math 225 THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Mat 225 As we ave seen, te definition of derivative for a Mat 111 function g : R R and for acurveγ : R E n are te same, except for interpretation:

More information

Global Existence of Classical Solutions for a Class Nonlinear Parabolic Equations

Global Existence of Classical Solutions for a Class Nonlinear Parabolic Equations Global Journal of Science Frontier Researc Matematics and Decision Sciences Volume 12 Issue 8 Version 1.0 Type : Double Blind Peer Reviewed International Researc Journal Publiser: Global Journals Inc.

More information

SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY

SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY (Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative

More information

Continuity and Differentiability of the Trigonometric Functions

Continuity and Differentiability of the Trigonometric Functions [Te basis for te following work will be te definition of te trigonometric functions as ratios of te sides of a triangle inscribed in a circle; in particular, te sine of an angle will be defined to be te

More information

NUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,

NUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example, NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing

More information

ERROR ESTIMATES FOR A FULLY DISCRETIZED SCHEME TO A CAHN-HILLIARD PHASE-FIELD MODEL FOR TWO-PHASE INCOMPRESSIBLE FLOWS

ERROR ESTIMATES FOR A FULLY DISCRETIZED SCHEME TO A CAHN-HILLIARD PHASE-FIELD MODEL FOR TWO-PHASE INCOMPRESSIBLE FLOWS ERROR ESTIMATES FOR A FULLY DISCRETIZED SCHEME TO A CAHN-HILLIARD PHASE-FIELD MODEL FOR TWO-PHASE INCOMPRESSIBLE FLOWS YONGYONG CAI, AND JIE SHEN Abstract. We carry out in tis paper a rigorous error analysis

More information

Math 31A Discussion Notes Week 4 October 20 and October 22, 2015

Math 31A Discussion Notes Week 4 October 20 and October 22, 2015 Mat 3A Discussion Notes Week 4 October 20 and October 22, 205 To prepare for te first midterm, we ll spend tis week working eamples resembling te various problems you ve seen so far tis term. In tese notes

More information

Mathematics 5 Worksheet 11 Geometry, Tangency, and the Derivative

Mathematics 5 Worksheet 11 Geometry, Tangency, and the Derivative Matematics 5 Workseet 11 Geometry, Tangency, and te Derivative Problem 1. Find te equation of a line wit slope m tat intersects te point (3, 9). Solution. Te equation for a line passing troug a point (x

More information

A Reconsideration of Matter Waves

A Reconsideration of Matter Waves A Reconsideration of Matter Waves by Roger Ellman Abstract Matter waves were discovered in te early 20t century from teir wavelengt, predicted by DeBroglie, Planck's constant divided by te particle's momentum,

More information

BOUNDARY REGULARITY FOR SOLUTIONS TO THE LINEARIZED MONGE-AMPÈRE EQUATIONS

BOUNDARY REGULARITY FOR SOLUTIONS TO THE LINEARIZED MONGE-AMPÈRE EQUATIONS BOUNDARY REGULARITY FOR SOLUTIONS TO THE LINEARIZED MONGE-AMPÈRE EQUATIONS N. Q. LE AND O. SAVIN Abstract. We obtain boundary Hölder gradient estimates and regularity for solutions to te linearized Monge-Ampère

More information

ANALYTIC SMOOTHING EFFECT FOR THE CUBIC HYPERBOLIC SCHRÖDINGER EQUATION IN TWO SPACE DIMENSIONS

ANALYTIC SMOOTHING EFFECT FOR THE CUBIC HYPERBOLIC SCHRÖDINGER EQUATION IN TWO SPACE DIMENSIONS Electronic Journal of Differential Equations, Vol. 2016 2016), No. 34, pp. 1 8. ISSN: 1072-6691. URL: ttp://ejde.mat.txstate.edu or ttp://ejde.mat.unt.edu ftp ejde.mat.txstate.edu ANALYTIC SMOOTHING EFFECT

More information

Continuity and Differentiability Worksheet

Continuity and Differentiability Worksheet Continuity and Differentiability Workseet (Be sure tat you can also do te grapical eercises from te tet- Tese were not included below! Typical problems are like problems -3, p. 6; -3, p. 7; 33-34, p. 7;

More information

Finite Difference Method

Finite Difference Method Capter 8 Finite Difference Metod 81 2nd order linear pde in two variables General 2nd order linear pde in two variables is given in te following form: L[u] = Au xx +2Bu xy +Cu yy +Du x +Eu y +Fu = G According

More information

Polynomial Interpolation

Polynomial Interpolation Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc

More information

A Finite Element Primer

A Finite Element Primer A Finite Element Primer David J. Silvester Scool of Matematics, University of Mancester d.silvester@mancester.ac.uk. Version.3 updated 4 October Contents A Model Diffusion Problem.................... x.

More information

Analysis of A Continuous Finite Element Method for H(curl, div)-elliptic Interface Problem

Analysis of A Continuous Finite Element Method for H(curl, div)-elliptic Interface Problem Analysis of A Continuous inite Element Metod for Hcurl, div)-elliptic Interface Problem Huoyuan Duan, Ping Lin, and Roger C. E. Tan Abstract In tis paper, we develop a continuous finite element metod for

More information

Finite Difference Methods Assignments

Finite Difference Methods Assignments Finite Difference Metods Assignments Anders Söberg and Aay Saxena, Micael Tuné, and Maria Westermarck Revised: Jarmo Rantakokko June 6, 1999 Teknisk databeandling Assignment 1: A one-dimensional eat equation

More information

A = h w (1) Error Analysis Physics 141

A = h w (1) Error Analysis Physics 141 Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.

More information

Lecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines

Lecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to

More information

New Streamfunction Approach for Magnetohydrodynamics

New Streamfunction Approach for Magnetohydrodynamics New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite

More information

MANY scientific and engineering problems can be

MANY scientific and engineering problems can be A Domain Decomposition Metod using Elliptical Arc Artificial Boundary for Exterior Problems Yajun Cen, and Qikui Du Abstract In tis paper, a Diriclet-Neumann alternating metod using elliptical arc artificial

More information

1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point

1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note

More information

Lyapunov characterization of input-to-state stability for semilinear control systems over Banach spaces

Lyapunov characterization of input-to-state stability for semilinear control systems over Banach spaces Lyapunov caracterization of input-to-state stability for semilinear control systems over Banac spaces Andrii Mironcenko a, Fabian Wirt a a Faculty of Computer Science and Matematics, University of Passau,

More information

A Mixed-Hybrid-Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems

A Mixed-Hybrid-Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems A Mixed-Hybrid-Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems

More information

Andrea Braides, Anneliese Defranceschi and Enrico Vitali. Introduction

Andrea Braides, Anneliese Defranceschi and Enrico Vitali. Introduction HOMOGENIZATION OF FREE DISCONTINUITY PROBLEMS Andrea Braides, Anneliese Defrancesci and Enrico Vitali Introduction Following Griffit s teory, yperelastic brittle media subject to fracture can be modeled

More information

MVT and Rolle s Theorem

MVT and Rolle s Theorem AP Calculus CHAPTER 4 WORKSHEET APPLICATIONS OF DIFFERENTIATION MVT and Rolle s Teorem Name Seat # Date UNLESS INDICATED, DO NOT USE YOUR CALCULATOR FOR ANY OF THESE QUESTIONS In problems 1 and, state

More information

Exam 1 Review Solutions

Exam 1 Review Solutions Exam Review Solutions Please also review te old quizzes, and be sure tat you understand te omework problems. General notes: () Always give an algebraic reason for your answer (graps are not sufficient),

More information

Strongly continuous semigroups

Strongly continuous semigroups Capter 2 Strongly continuous semigroups Te main application of te teory developed in tis capter is related to PDE systems. Tese systems can provide te strong continuity properties only. 2.1 Closed operators

More information

Solution. Solution. f (x) = (cos x)2 cos(2x) 2 sin(2x) 2 cos x ( sin x) (cos x) 4. f (π/4) = ( 2/2) ( 2/2) ( 2/2) ( 2/2) 4.

Solution. Solution. f (x) = (cos x)2 cos(2x) 2 sin(2x) 2 cos x ( sin x) (cos x) 4. f (π/4) = ( 2/2) ( 2/2) ( 2/2) ( 2/2) 4. December 09, 20 Calculus PracticeTest s Name: (4 points) Find te absolute extrema of f(x) = x 3 0 on te interval [0, 4] Te derivative of f(x) is f (x) = 3x 2, wic is zero only at x = 0 Tus we only need

More information

Material for Difference Quotient

Material for Difference Quotient Material for Difference Quotient Prepared by Stepanie Quintal, graduate student and Marvin Stick, professor Dept. of Matematical Sciences, UMass Lowell Summer 05 Preface Te following difference quotient

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL IFFERENTIATION FIRST ERIVATIVES Te simplest difference formulas are based on using a straigt line to interpolate te given data; tey use two data pints to estimate te derivative. We assume tat

More information

Stability properties of a family of chock capturing methods for hyperbolic conservation laws

Stability properties of a family of chock capturing methods for hyperbolic conservation laws Proceedings of te 3rd IASME/WSEAS Int. Conf. on FLUID DYNAMICS & AERODYNAMICS, Corfu, Greece, August 0-, 005 (pp48-5) Stability properties of a family of cock capturing metods for yperbolic conservation

More information

Click here to see an animation of the derivative

Click here to see an animation of the derivative Differentiation Massoud Malek Derivative Te concept of derivative is at te core of Calculus; It is a very powerful tool for understanding te beavior of matematical functions. It allows us to optimize functions,

More information

OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS. Sandra Pinelas and Julio G. Dix

OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS. Sandra Pinelas and Julio G. Dix Opuscula Mat. 37, no. 6 (2017), 887 898 ttp://dx.doi.org/10.7494/opmat.2017.37.6.887 Opuscula Matematica OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS Sandra

More information

Generic maximum nullity of a graph

Generic maximum nullity of a graph Generic maximum nullity of a grap Leslie Hogben Bryan Sader Marc 5, 2008 Abstract For a grap G of order n, te maximum nullity of G is defined to be te largest possible nullity over all real symmetric n

More information

Notes on wavefunctions II: momentum wavefunctions

Notes on wavefunctions II: momentum wavefunctions Notes on wavefunctions II: momentum wavefunctions and uncertainty Te state of a particle at any time is described by a wavefunction ψ(x). Tese wavefunction must cange wit time, since we know tat particles

More information

1 Introduction. FILOMAT (Niš) 16 (2002), GEODESIC MAPPINGS BETWEEN KÄHLERIAN SPACES. Josef Mikeš, Olga Pokorná 1 and Galina Starko

1 Introduction. FILOMAT (Niš) 16 (2002), GEODESIC MAPPINGS BETWEEN KÄHLERIAN SPACES. Josef Mikeš, Olga Pokorná 1 and Galina Starko FILOMAT (Niš) 16 (2002), 43 50 GEODESIC MAPPINGS BETWEEN KÄHLERIAN SPACES Josef Mikeš, Olga Pokorná 1 and Galina Starko Abstract Geodesic mappings from a Kälerian space K n onto a Kälerian space K n will

More information

A compactness theorem for Yamabe metrics

A compactness theorem for Yamabe metrics A compactness theorem for Yamabe metrics Heather acbeth November 6, 2012 A well-known corollary of Aubin s work on the Yamabe problem [Aub76a] is the fact that, in a conformal class other than the conformal

More information

University Mathematics 2

University Mathematics 2 University Matematics 2 1 Differentiability In tis section, we discuss te differentiability of functions. Definition 1.1 Differentiable function). Let f) be a function. We say tat f is differentiable at

More information

Nonlinear elliptic-parabolic problems

Nonlinear elliptic-parabolic problems Nonlinear elliptic-parabolic problems Inwon C. Kim and Norbert Požár Abstract We introduce a notion of viscosity solutions for a general class of elliptic-parabolic pase transition problems. Tese include

More information

Function Composition and Chain Rules

Function Composition and Chain Rules Function Composition and s James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 8, 2017 Outline 1 Function Composition and Continuity 2 Function

More information

Practice Problem Solutions: Exam 1

Practice Problem Solutions: Exam 1 Practice Problem Solutions: Exam 1 1. (a) Algebraic Solution: Te largest term in te numerator is 3x 2, wile te largest term in te denominator is 5x 2 3x 2 + 5. Tus lim x 5x 2 2x 3x 2 x 5x 2 = 3 5 Numerical

More information

arxiv: v1 [math.dg] 4 Feb 2015

arxiv: v1 [math.dg] 4 Feb 2015 CENTROID OF TRIANGLES ASSOCIATED WITH A CURVE arxiv:1502.01205v1 [mat.dg] 4 Feb 2015 Dong-Soo Kim and Dong Seo Kim Abstract. Arcimedes sowed tat te area between a parabola and any cord AB on te parabola

More information

Work and Energy. Introduction. Work. PHY energy - J. Hedberg

Work and Energy. Introduction. Work. PHY energy - J. Hedberg Work and Energy PHY 207 - energy - J. Hedberg - 2017 1. Introduction 2. Work 3. Kinetic Energy 4. Potential Energy 5. Conservation of Mecanical Energy 6. Ex: Te Loop te Loop 7. Conservative and Non-conservative

More information

Linearized Primal-Dual Methods for Linear Inverse Problems with Total Variation Regularization and Finite Element Discretization

Linearized Primal-Dual Methods for Linear Inverse Problems with Total Variation Regularization and Finite Element Discretization Linearized Primal-Dual Metods for Linear Inverse Problems wit Total Variation Regularization and Finite Element Discretization WENYI TIAN XIAOMING YUAN September 2, 26 Abstract. Linear inverse problems

More information

Applications of the van Trees inequality to non-parametric estimation.

Applications of the van Trees inequality to non-parametric estimation. Brno-06, Lecture 2, 16.05.06 D/Stat/Brno-06/2.tex www.mast.queensu.ca/ blevit/ Applications of te van Trees inequality to non-parametric estimation. Regular non-parametric problems. As an example of suc

More information

NONNEGATIVITY OF EXACT AND NUMERICAL SOLUTIONS OF SOME CHEMOTACTIC MODELS

NONNEGATIVITY OF EXACT AND NUMERICAL SOLUTIONS OF SOME CHEMOTACTIC MODELS NONNEGATIVITY OF EXACT AND NUMERICAL SOLUTIONS OF SOME CHEMOTACTIC MODELS PATRICK DE LEENHEER, JAY GOPALAKRISHNAN, AND ERICA ZUHR Abstract. We investigate nonnegativity of exact and numerical solutions

More information

Quasiperiodic phenomena in the Van der Pol - Mathieu equation

Quasiperiodic phenomena in the Van der Pol - Mathieu equation Quasiperiodic penomena in te Van der Pol - Matieu equation F. Veerman and F. Verulst Department of Matematics, Utrect University P.O. Box 80.010, 3508 TA Utrect Te Neterlands April 8, 009 Abstract Te Van

More information