arxiv: v2 [math.dg] 21 Nov 2013

Transcription

1 arxiv:11.579v [math.dg] 1 Nov 13 Fractional conformal Laplacians and fractional Yamabe problems María del Mar González Univ. Politècnica de Catalunya Abstract Jie Qing Univ. of California, Santa Cruz Based on the relations between scattering operators of asymptotically hyperbolic metrics and Dirichlet-to-Neumann operators of uniformly degenerate elliptic boundary value problems observed by Chang and González, we formulate fractional Yamabe problems that include the boundary Yamabe problem studied by Escobar. We observe an interesting Hopf type maximum principle together with interplays between analysis of weighted trace Sobolev inequalities and conformal structure of the underlying manifolds, which extend the phenomena displayed in the classic Yamabe problem and boundary Yamabe problem. 1 Introduction In this paper, based on the relations between scattering operators of asymptotically hyperbolic metrics and Dirichlet-to-Neumann operators of uniformly degenerate elliptic boundary value problems observed in [11], we formulated and solved fractional order Yamabe problems that include the boundary Yamabe problem studied by Escobar in [16]. Suppose that X n1 is a smooth manifold with smooth boundary M n for n 3. A function ρ is a defining function of the boundary M n in X n1 if ρ > in X n1, ρ = on M n, dρ on M n. We say that g is conformally compact if, for some defining function ρ, the metric ḡ = ρ g extends to X n1 so that ( X n1,ḡ) is a compact iemannian manifold. This induces a conformal class of metrics ĥ = ḡ TM n on Mn when defining functions vary. The conformal manifold (M n,[ĥ]) is called the conformal infinity of (Xn1,g ). A metric g is said to be asymptotically hyperbolic if it is conformally compact and the sectional curvature approaches 1 at infinity. In the recent work[7], Graham and Zworski introduced the meromorphic family of scattering operators S(s), which is a family of pseudo-differential operators, for a given asymptotically hyperbolic manifold (X n1, g ) and a choice of the representative ĥ of the conformal infinity (M n, [ĥ]). Often one instead considers the normalized scattering operators Γ(γ) ( n ) P γ [g,ĥ] = γ Γ( γ) S γ. The normalized scattering operators P γ [g,ĥ] are conformally covariant, P γ [g,w 4 n γ ĥ]φ = w nγ n γ Pγ [g,ĥ](wφ), González is supported by Spain Government project MTM C3-1 and GenCat 9SG345. esearch of Qing is partially supported by NSF 7535 and CNSF

2 with principal symbol σ(p γ [g,ĥ]) = σ(( ĥ )γ ). Hence they may be considered to be conformal fractional Laplacians for γ (,1) for a given asymptotically hyperbolic metric g. As proven in [7], [19], when g is Poincaré-Einstein, P 1 is the conformal Laplacian, P is the Paneitz operator, and in general P k for k N are the conformal powers of the Laplacian discovered in [6]. When g is a fixed asymptotically hyperbolic metric we may simply denote Pĥγ := P γ [g,ĥ]. We will consider the associated fractional order curvature and the normalized total curvature Qĥγ = Pĥγ (1), I γ [ĥ] = M n Qĥγdvĥ ( M n dvĥ When a background metric ĥ is fixed, we may write I γ [w,ĥ] = I γ[w 4 n γ ĥ] = )n γ n. M n wpĥγ wdvĥ ( M n w n n γ dvĥ ) n γ n This functional I γ [ĥ] is clearly an analogue to the Yamabe functional. Hence one may ask if there is a metric which is the minimizer of I γ among metrics in the class [ĥ] and whose curvature Q γ is a constant. We will refer to that problem as a fractional Yamabe problem when γ (,1). For the original Yamabe problem readers are refereed to [3], [4]. A similar question was studied in [39] for γ > 1 and g being a Poincaré-Einstein metric. Because of the lack of a maximum principle these generalized Yamabe problems in general are difficult to solve. Yet this new window to the analytic aspects of conformal geometry remains fascinating. For example, it was proven in [8] that the location of the first scattering pole is dictated by the sign of the Yamabe constant and the Green s function of Pĥγ is positive for γ (,1) when the Yamabe constant is positive, at least in the case where g is conformally compact Einstein. It turns out that one may use the relations of scattering operators and the Dirichletto-Neumann operators to reformulate the above fractional Yamabe problems as degenerate elliptic boundary value problems. The correspondence between pseudo-differential equations and degenerate elliptic boundary value problems is inspired by the works in [1]. Interestingly, the corresponding degenerate elliptic boundary value problem is a natural extension of the boundary Yamabe problem raised and studied in [16]. ecall from [11] that, given an asymptotically hyperbolic manifold (X n1, g ) and a representative ĥ of the conformal infinity (Mn, [ĥ]), one can find a geodesic defining function ρ such that the compactified metric can be written as ḡ := ρ g = dρ h ρ = dρ ĥh(1) ρh () ρ o(ρ ) near infinity. One may consider the degenerate elliptic boundary value problem of ḡ as follows: { div(ρ a U)E(ρ)U = in (X n1, ḡ), U ρ= = f on M n (1.1), where s = n γ, and a = 1 γ. E(ρ) = ρ 1 s( g s(n s) ) ρ n s,.

3 Lemma 1.1 (Chang and Gonzalez [11]). Let (X n1, g ) be an asymptotically hyperbolic manifold. Suppose that U is the solution to the boundary value problem (1.1). Then 1. For γ (, 1 ) and n 4 γ not an L -eigenvalue for the Laplacian of g, P γ [g,ĥ]f = d γ lim ρ ρ a ρ U, (1.) where d γ = γ 1 Γ(γ) γγ( γ). (1.3). For γ = 1, P1[g,ĥ]f = lim ρu n 1 ρ Hf, where H := 1 n Tr ĥ (h(1) ) is the mean curvature of M. 3. For γ ( 1,1), (1.) still holds if H =. In light of Lemma 1.1, consider, for γ (,1), Iγ [U,ḡ] = d γ (ρ a U E(ρ)U X n1 )dvḡ. U M n n n γ dvĥ It is then a very natural variational problem for I γ. For instance, right away one sees that a minimizer of I γ is automatically nonnegative, which was a huge issue for the functional I γ. OnekeyingredientinourworkhereisthefollowingHopftypemaximumprinciple. Wedrew inspiration from some version of Hopf s lemma for the Euclidean half space case (Proposition 4.11 in [9]). Proposition 1.. Let γ (,1). Suppose that U is a nonnegative solution to (1.1) in X n1. Let p M n = X n1 and B r be a geodesic ball of radius r centered at p in M n. Then, for sufficiently small r, if U(q ) = for q B r \B1 and U > on r B1 r, then y a y U q >. (1.4) It seems weaker than the original one, but it suffices for our purposes. A nice and immediate consequence of the above maximum principle is that the first eigenfunction of the fractional conformal Laplacian Pĥγ is always positive, which has been a rather challenging question in general for the pseudo-differential operators Pĥγ (cf. [8]). Hence one can produce a metric in the class [ĥ] that has positive, negative, or zero Q γ curvature when the first eigenvalue is positive, negative, or zero respectively. Our approach to solve the γ-yamabe problem is very similar to the one taken in [16], where one of the crucial steps is the understanding of a trace inequality. In our case, the relevant sharp weighted trace Sobolev inequality appeared in the works [31], [13], [37]: Proposition 1.3. Let γ (,1) and a = 1 γ. Suppose that U W 1, ( n1,ya ) with trace TU = w. Then, for some constant S(n,γ), w L ( n ) S(n,γ) n1 where = n n γ. Moreover the equality holds if and only if w(x) = c ( µ x x µ y a U dxdy, (1.5) )n γ, x n, for c, µ > and x n fixed, and U is its Poisson extension of w as given in (.13).

4 As in the case of original Yamabe problem, one can define the γ-yamabe constant n Λ γ (M,[ĥ]) = inf I γ [h]. h [ĥ] It is then easily seen that Λ γ (S n,[g c ]) = d γ S(n, γ) where [g c ] is the canonical conformal class of metrics on the sphere S n. Analogous to the cases of the original Yamabe problem we obtain Theorem 1.4. Suppose that (X n1,g ) is an asymptotically hyperbolic manifold. Suppose, in addition, that H = when γ ( 1,1). Then, if then the γ-yamabe problem is solvable for γ (,1). < Λ γ (M,[ĥ]) < Λ γ(s n,[g c ]), (1.6) emark. It is easily seen that Λ γ (M,[ĥ]) > in the light of (1.4) in Theorem 1.1 and Theorem 1. in [9] when γ (, 1 ] or if some additional assumptions in Theorem 1. in [9] hold. Based on computations similar to ones in [16], we have Theorem 1.5. Suppose that (X n1, g ) is an asymptotically hyperbolic manifold and that ρ ( [g ] ic[g ](ρ ρ )n ) as ρ. (1.7) If X n1 has a non-umbilic point on X n1 and then na 3 Γ(γ) γ1 1 a Γ( γ) n 1a <, (1.8) a1 Λ γ (M,[ĥ]) < Λ γ(s n,[g c ]) and hence the γ-yamabe problem is solvable for γ (,1). We remark now that the 1 -Yamabe problem introduced in here reduces back to the boundary Yamabe problem consider in [16] in this way. Notice that, in this case, we have I 1[U, φn 1ḡ] 4 = I 1 [Uφ, ḡ] (1.9) for any positive function φ on X n1 and therefore (1.7) is no longer needed. Also notice that the condition (1.8) becomes n > 5 when γ = 1, which agrees with the conclusion in [16]. Suppose we start with a compact iemannian manifold (X n1, ḡ) and its boundary (M n, ĥ). Then one can construct an asymptotically hyperbolic manifold (Xn1, g ) which is conformal to (X n1, ḡ). For example, as observed in [11], one may require according to the works in [33], [] that [g ] = n(n1). (1.1) Then the induced degenerate equation becomes div(ρ a U) n 1a 4n [ḡ]ρa U = in (X n1,ḡ) (1.11) whose associated variational functional becomes F[U] = ρ a U ḡ dv ḡ n 1a 4n X X [ḡ]ρ a U dvḡ. (1.1)

5 In section we recall the work from [11] to make possible the passage from pseudodifferential equations to second order elliptic boundary value problems as in [1]. In Section 3 we study regularity (L and Schauder estimates) for degenerate elliptic boundary value problems. And more importantly we establish the Hopf type maximum principle. In Section 4 we formulate the fractional Yamabe problem and obtain some properties for the fractional case that are analogous to the original Yamabe problem with the help of the Hopf-type maximum principle. In Section 5 we analyze sharp weighted Sobolev trace inequalities. We define, on any conformal manifold, the fractional Yamabe constant associated with an asymptotically hyperbolic metric and show that the one of the standard round spheres associated to the standard hyperbolic metric is the largest. In Section 6 we take a subcritical approximation and prove our Theorem 1.4. In the last section we adopt the calculation from [16] and prove our Theorem 1.5 by choosing a suitable test function. We finally mention the two related works [4, 41] on nonlinearities with critical exponents for the fractional Laplacian. Conformal fractional Laplacians In this section we introduce the recent works in [11] to relate two equivalent definitions of conformal fractional Laplacians. Conformal fractional Laplacians are defined via scattering theory on asymptotically hyperbolic manifolds in [7], [19]. We also have seen fractional Laplacians defined as Dirichlet-to-Neumann operators for degenerate equations on compact manifolds with boundary in [1]. It turns out in some way these two fractional Laplacians are the same. Let X n1 be a smooth manifold of dimension n1 with compact boundary X = M n. A function ρ is a defining function of X in X if ρ > in X, ρ = on X, dρ on X. We say that g is conformally compact if the metric ḡ = ρ g extends to X n1 for a defining function ρ so that ( X n1,ḡ) is a compact iemannian manifold. This induces a conformal class of metrics ĥ = ḡ TM n on Mn when the defining function varies, which is called the conformal infinity of (X n1,g ). A metric g is said to be asymptotically hyperbolic if it is conformally compact and the sectional curvature approaches to 1 at infinity. Given an asymptotically hyperbolic manifold (X n1,g ) and a representative ĥ of the n conformal infinity (M,[ĥ]), there is a uniquely geodesic defining function ρ such that, on a neighborhood M (,δ) in X, g has the normal form where h ρ is a one parameter family of metrics on M such that g = ρ (dρ h ρ ) (.1) h ρ = ĥh(1) ρo(ρ ). (.) From [34], [7] it follows that, given f C (M), e(s) > n and s(n s) is not a L - eigenvalue for g, the generalized eigenvalue problem has a solution of the form g u s(n s)u =, in X (.3) u = Fρ n s Gρ s, F,G C ( X), F ρ= = f. (.4) The scattering operator on M is then defined as S(s)f = G M.

6 It is shown in [7] that, by a meromorphic continuation, S(s) is a meromorphic family of pseudo-differential operators in the whole complex plane. Instead, it is often useful to consider the normalized scattering operators P γ [g,ĥ] defined as: ( n ) P γ [g,ĥ] := d γs γ γ Γ(γ), d γ = Γ( γ). (.5) Note that s = n γ. With this regularization the principal symbol of P γ[g,ĥ] is exactly the principal symbol of the fractional Laplacian ( ĥ) γ. Hence we will call (assuming implicitly the dependence on the extension metric g ) Pĥγ := P γ[g,ĥ] a conformal fractional Laplacian for each γ (, 1) which is not a pole of the scattering n operator, i.e. 4 γ is not a L -eigenvalue for g. It is a conformally covariant operator, in the sense that it behaves like for a conformal change of metric ĥw = w 4 n γ ĥ. We will call Pĥw γ ϕ = w nγ n γ P ĥ γ (wϕ) (.6) Qĥγ = Pĥγ (1) the fractional scalar curvature associated to the conformal fractional Laplacian Pĥγ. From the above (.6) we have Pĥγ (w) = Qĥw γ w nγ n γ. (.7) The familiar case is γ = 1, where Pĥ1 = ĥ n 4(n 1) [ĥ] becomes the conformal Laplacian and the associated curvature is the scalar curvature Qĥ1 = [ĥ] of the metric ĥ which undergoes the change n 4(n 1) Pĥ1 w = n n [ĥw]w n 4(n 1) when taking conformal change of metrics, provided that (X n1, g ) is a Poincaré-Einstein as established in [7], [19]. The conformal fractional Laplacians and fractional scalar curvatures should also be compared to the higher order generalization of the conformal Laplacian and scalar curvature: the Paneitz operator Pĥ and its associated Q-curvature (see [38], [6], [39]). It was observed by Chang and González in[11] that the generalized eigenvalue problem(.3) on a non-compact manifold (X n1,g ) is equivalent to a linear degenerate elliptic problem on the compact manifold ( X n1,ḡ), for ḡ = ρ g. Hence Chang and González reconciled the definition of the fractional Laplacians given in the above as normalized scattering operators and the one given in the spirit of the Dirichlet-to-Neumann operators by Caffarelli and Silvestre in [1]. This observation in [11] plays a fundamental role in this paper and provides an alternative way to study the fractional partial differential equation (.7). First, we know by the conformal covariance that P g 1 u = ρn3 P ḡ 1 (ρ n 1 u). Let a = 1 γ ( 1,1), s = n γ, and U = ρs n u. Then we may write the equation (.3) as div(ρ a ḡu)e(ρ)u =, in (X n1, ḡ),

7 where or writing everything back in the metric g, E(ρ) := ρ a P ḡ 1 ρa ( s(n s) n 1 4n [g ] ) ρ a, (.8) E(ρ) = ρ 1 s( g s(n s) ) ρ n s. (.9) Notice that, in a neighborhood M (,δ) where the metric g is in the normal form, [ E(ρ) = n 1a 4n [ḡ] (n(n1)[g ])ρ ] ρ a in M (,δ). (.1) Proposition.1 (Chang and González [11]). Let (X n1, g ) be an asymptotically hyperbolic manifold. Then, given f C (M), the generalized eigenvalue problem (.3)-(.4) is equivalent to { div(ρ a U)E(ρ)U = in (X,ḡ), (.11) U ρ= = f on M, where U = ρ n s u and U is the unique minimizer of the energy F[V] = ρ a V ḡ dv ḡ E(ρ) V dvḡ X among all the functions V W 1, (X,ρ a ) with fixed trace V ρ= = f. Moreover, X 1. For γ (, 1 ), Pĥ γf = d γ lim ρ a ρ U, (.1) ρ where the constant d γ is given in (1.3).. For γ = 1, we have an extra term Pĥ1 f = lim ρ U n 1 ρ Hf, where H := 1 n Tr ĥ (h(1) ) is the mean curvature of M. 3. For γ ( 1,1), (.1) still holds if and only if H =. emark. It should be noted here that there are many asymptotically hyperbolic manifolds (X n1, g ) whose conformal infinity is prescribed as (M n, [ĥ]). If one insists (Xn1, g ) to be Poincaré-Einstein, then the normalized scattering operators Pĥγ are a bit more intrinsic, at least at positive integers as observed in [7], [19]. It should also be noted that one can simply start with a compact iemannian manifold ( X n1, ḡ) with boundary (M n, ĥ) and easily build an asymptotically hyperbolic manifold whose conformal infinity is given by (M n, [ĥ]). Please see the details of this observation in [11]. The simplest example of a conformally compact Einstein manifold is the hyperbolic space (H n1,g H ). It can be characterizedasthe upper half-space (with coordinatesx n, y ), endowed with the metric: g = dy dx y. Then (.11) with Dirichlet condition w reduces to { div(y a U) = in n1, U y= = w on n, and the fractional Laplacian at the boundary n is just P dx γ w = ( dx ) γ w = d γ lim y (y a y U).

8 This is precisely the Caffarelli-Silvestre extension [1]. Note that this extension U can be written in terms of the Poisson kernel K γ as follows: U(x,y) = K γ x w = C n,γ n y 1 a ( x ξ y ) n1 a w(ξ) dξ, (.13) for some constant C n,γ. Moreover, given w H γ ( n ), U is the minimizer of the functional: F[V] = y a V dxdy n1 among all the possible extensions in the set { } V : n1 : y a V dxdy <, V(,) = w. n1 Based on (.9) it is observed in [11] that one may use as a defining function, where v solves ρ = v 1 n s g v s(n s)v = and ρ s n v = 1 on M, to eliminate E(ρ ) from equation (.11). It suffices to show that v is strictly positive in the interior. But this is true because, away from the boundary, it is the solution of an uniformly elliptic equation in divergence form, thus it cannot have a non-positive minimum. Hence we arrive at an improvement of Proposition.1 as follows: Proposition.. The function ρ is a defining function of M in X such that E(ρ ). Hence U = (ρ ) s n u solves { div((ρ ) a U) = in (X,ḡ ), (.14) U = w on M, with respect to the metric ḡ = (ρ ) g and U is the unique minimizer of the energy F[V] = (ρ ) a V ḡ dv ḡ (.15) among all the extensions V W 1, (X,(ρ ) a ) satisfying V M = w. Moreover, [ ] ρ Qĥγ (ρ) = ρ 1 ρ γ O(ρ ) (n s)d γ near the infinity and provided that H = when γ ( 1,1). X Pĥγ w = d γ lim ρ (ρ ) a ρ U wqĥγ, (.16) We will sometimes use the defining function ρ, denoted by y unless explicitly stated otherwise, because it allowsus to workwith a pure divergenceequation with no lowerorderterms. We end this section by discussing the assumption that H = for an asymptotically hyperbolic metric g. It turns out that this indeed is an intrinsic condition.

9 Lemma.3. Suppose that (X n1, g ) is an asymptotically hyperbolic manifold and that ρ and ρ are the geodesic defining functions of M in X associated with representatives ĥ and h of the conformal infinity (M n, [ĥ]) respectively. Hence g = ρ (dρ h ρ ) = ρ (d ρ h ρ ) where and near the infinity. Then In particular h ρ = ĥρh(1) O(ρ ) h ρ = h ρ h (1) O( ρ ) h (1) = h (1) on M. H = ρ ρ H on M. ρ= Proof. This simply follows from the equations that define the geodesic defining functions. Let ρ = e w ρ near the infinity. Then 1 = d(e w ρ) e w ρ g = dρ ρ g ρ dw,dρ ρ g ρ dw ρ g, which implies w ρ ρ [ ( w ρ ) w h ρ ] =. Hence it is rather obvious that w ρ = at ρ =. Therefore the proof is complete in the light of the fact that g = ρ g = e w ρ g = e w ḡ. 3 Uniformly degenerate elliptic equations Considering the fractional powers of the Laplacian as Dirichlet-to-Neumann operators in Proposition. allows to relate the properties of non-local operators to those of uniformly degenerate elliptic equations in one more dimension. The same strategy has been used, for instance, in the recent work of Cabré-Sire [9]. Fix γ (,1). Let y = ρ be the special defining function given in Proposition. and set ḡ = y g. We are concerned with the uniformly degenerate elliptic equation { div(y a U) = in (X,ḡ ), (3.1) U = w on M. For our purpose we will concentrate on the local behaviors of the solutions to (3.1) near the boundary. First, wewrite our equationin localcoordinatesnearafixed boundarypoint (p,). More precisely, for some >, we set B = {(x,y) n1 : y >, (x,y) < }, Γ = {(x,) n1 : x < }, Γ = {(x,y) n1 : y, (x,y) = }.

10 In local coordinates on Γ the metric ĥ is of the form dx (1 O( x )), where x(p ) =. Consider the matrix A(x,y) = detḡ y a (ḡ ) 1. Then the equation (3.1) is equivalent to n1 i,j=1 i (A ij j U) =. (3.) Moreover we know that 1 c ya I A cy a I. (3.3) This shows that (3.) is a uniformly degenerate elliptic equation. For instance, the weight ψ(y) = y a is an A weight in the sense of [36]. Equation (3.) has been well understood in a series of papers by Fabes, Jerison, Kenig, Serapioni ([18], [17]). Let us state a regularity result that is relevant to us. We will concentrate on problems of the form { Div(A(DU)) = in B, y a y U = F, on Γ, (3.4) where, for the rest of the section, A satisfies the ellipticity condition (3.3) for a ( 1,1), the derivatives are Euclidean, that is, D := ( x1,..., xn,y), and Div(A(DU)) := n1 i,j=1 i (A ij j U). Definition 3.1. Given > and a function F L 1 (Γ ), we call U a weak solution of (3.4) if U satisfies (DU) t A(DU) L 1 (B ) and B (Dφ) t A(DU)dxdy Fφdx = Γ for all φ C 1 (B ) such that φ on Γ and (Dξ)t A(Dφ) L 1 (B ). Hölder regularity for weak solutions was shown in [18], Lemma.3.1, for any A satisfying (3.3). Using this main result, regularity of weak solutions up to the boundary was carefully shown in [9], Lemma 4.3, at least when A = y a I. However, their proof only depends on the divergence structure of the equation and the behavior of the weight. Hence we have Proposition 3.. Let γ (,1), γ = 1 a and β (,min{1,1 a}). Let > and U L (B ) W 1, (B,ya ) be a weak solution of { Div(A(DU)) = in B, y a y U = F(U) on Γ, (3.5) for A satisfying (3.3). If F C 1,β, then U C, β(b ) and x i U C, β(b ), i = 1,...,n, for some β (,1). Particularly, when F(x,t) = α(x)tβ(x)t nγ n γ, to get smoothness it is necessary to know the local boundedness of weak solutions U on B. To get this local boundedness for weak solutions we employ the usual Moser s iteration scheme adapted to boundary valued problems (see Theorem 3.4 below). However, a new idea is required: we will perform two coupled iterations, one in the interior and one at the boundary, that need to be handled simultaneously. Note that in the linear case when F, local boundedness was shown in [18, Corollary.3.4], using the weighted Sobolev embeddings in the interior described in Proposition 3.3. However,

11 when a non-linearity F(U) is present at the boundary term, instead we need to use weighted trace Sobolev embeddings. First, we recall a weighted Sobolev embedding theorem in the interior (c.f. [18, Theorem 1.3], see also [1]): Proposition 3.3. Let Ω be an open bounded set in n1. Take 1 < p <. There exist positive constants C Ω and δ such that for all u C (Ω) and all k satisfying 1 k n1 n δ, u L kp (Ω,y a ) C Ω u L p (Ω,y a ). C Ω maybe taken to depend only on n, p, a and the diameter of Ω. Now we can state the theorem. Note that we actually prove it in the flat case but it is straightforward to generalize it to the manifold setting: Theorem 3.4. Let U be a weak solution of the problem { div(y a U) = in B, where F(z) satisfies F(z) = O y a y U = F(U) on Γ, ( z β 1), when z, for some < β <. (3.6) Assume, in addition, that Γ r U dx =: V <. Then for each p > 1, there exists a constant C p = C( p,v) > such that sup B [ ( 1 U sup U C p Γ n1a ) 1/ p ( ) ] 1/ p 1 U L p(b,ya) n U L p(γ ). Proof. Let p X. Note that we can work with normal coordinates x 1,...,x n n, y > near p. Without loss of generality, assume that = 1. Then the general case is obtained by rescaling. Let η = η(r), r = ( x y ) 1/, be a smooth cutoff function such that η = 1 if r < 1, η = if r, η 1 if r (1,). Next, by working with U := max{u,}, U := max{ U,} separately, we can assume that U is positive. A good reference for Moser iteration arguments in divergence structure equations is [, chapter 8]. We generalize this method, considering a double iteration: one at the boundary, using Sobolev trace inequalities to handle the non-linear term F(U), the other in the interior domain. The first step is to use that U is a weak solution of (3.6) by finding a good test function. Formally we can write the following: multiply equation (3.6) by η U α and integrate by parts: = y a ηu α η U dxdy α y a η U α 1 U dxdy η U α F(U)dx. (3.7) B B This implies, using Hölder estimates to handle the crossed term, y a η U α 1 U dxdy η U α F(U)dx 4 α α y a η U α1 dxdy. (3.8) B On the other hand, again using Hölder inequality, we have y a (ηu δ ) dxdy δ y a η U (δ 1) U dxdy B B Γ B Γ B y a U δ η dxdy. If we insert formula (3.8) into the inequality above, for the choice α = δ 1, we obtain J := y a (ηu δ ) dxdy B ( 1 ( α1 α =: I 1 I. ) ) B y a η U δ dxdy (α1) α Γ η U α F(U)dx (3.9)

12 For the left hand side above, recall the trace Sobolev embedding (Corollary 5.3): J = B ( y a (ηu δ ) dxdy Γ ) (ηu δ ) dx, (3.1) and the standard weighted Sobolev embedding from Proposition 3.3. J = B y a (ηu δ ) dxdy ( B y a (ηu δ ) k ) k (3.11) for some 1 < k < n1 n. Next, we estimate from above the terms I 1,I in (3.9). I 1 can be easily handled since η C: I 1 = y a η U δ dxdy y a U δ dxdy. (3.1) B Now we consider the second term. To estimate I, if we write U δ β = U β U δ, then using Hölder inequality with p = β, 1 p 1 q = 1, we obtain Γ η U δ 1 F(U)dx [ Γ U dx ] 1 p [ Γ B η q U δq dx ] 1 q V 1 p [ Γ η q U δq dx ] 1 q. (3.13) Thislastintegralcanbehandledasfollows. Callχ =, forsimplicity. Becauseourhypothesis on β, we know that q (1,χ). Then, there exists λ (,1) such that q = λ(1 λ)χ, and an interpolation inequality gives: [ f q ]1 q [ f ]λ [ ]1 λ [ q f χ q = ] 1 ( f χ χ [ f][ ] 1)λ q f χ χ. (3.14) Since λ q < 1, Young s inequality reads z λ q Cǫ z ǫ, for ǫ small. If we substitute z = [ f ][ f χ] 1 χ above, together with (3.14), we arrive at [ f q ] 1 q ǫ [ f χ ] 1 χ Cǫ Then from (3.13) it follows that ( ) I V 1 p ǫ (ηu δ ) dx C ǫ η U δ dx Γ Γ, (3.15) where ǫ will be chosen later and will depend on the value of α,δ. We go back now to the main iteration formula (3.9). It is clear from (3.1), that the first integral of the right hand side of the formula for I (3.15) can be absorbed into the left hand side of (3.9), and using (3.11) and (3.1) we get that ( Γ 1 U δ dx ) ( B 1 U kδ dxdy ) 1 k C(δ) for some suitable choice of ǫ. Or switching notation from δ to δ, ( Γ 1 U δχ dx ) 1 χ ( B 1 U kδ dxdy ) 1 k C(δ) [ [ Γ Γ f. U δ dx B U δ dxdy ] U δ dx U δ dxdy. (3.16) B ],

13 Next, because we will always have δ > 1, we can use that so from (3.16) we get that For simplicity, we set and C 1 (a 1 δ b 1 δ ) (ab) 1 δ C (a 1 δ b 1 δ ), U L χδ (Γ 1 ) U L kδ (B 1,ya ) U L δ (Γ ) U L δ (B,ya ). Φ(δ,) := θ := min{χ,k} > 1, ( ) 1 ( 1 δ U L n δ (Γ 1 ) 1 n1a )1 δ U L δ (B 1,ya ). Then, after explicitly writing all the constants involved, formula (3.16) simply reduces to Φ(θδ,1) [C(1δ) σ ] δ Φ(δ,), for some positive number σ. It is clear that the same proof works if we replace B 1,B by B 1, B. The only difference is in (3.1), where we need to estimate η C( 1 ) 1. Thus we would obtain [ ] C(1δ) σ δ Φ(θδ, 1 ) Φ(δ, ). (3.17) 1 Now we iterate equation (3.17): set m = 1 1 m and θ m = θ m p. Then m 1 Φ(θ m,1) Φ(θ m, m ) (c 1 θ) c i= for some constant C because the series Finally, note that so that (3.18) is telling us that i= i θ i is convergent. i θ i Φ( p,) CΦ( p,), (3.18) supu = lim U L Γ δ δ (Γ 1 ), supu = lim U 1 B δ L δ (B 1,ya ), 1 sup B 1 U supu C Γ 1 [ U L p(b,y a) U L p(γ ) ]. escaling to a ball of radius concludes the proof of the theorem. The next main ingredient is the proof of the positivity of a solution to (3.5). We observed that a Hopf lemma, some version of which was known for the Euclidean half space case (Proposition 4.1 in [9]), can be obtained for the uniformly degenerate elliptic equation (3.1). This nice Hopf s lemma turns out to be one of the keys for us in this paper. It is interesting to observe a different behavior between the cases γ (,1/) and γ [1/,1) in our proof - this dichotomy does not seem to appear in the flat case in [9]. We continue to use the setting as in Proposition.. Let p X and (x,y) be the local coordinate at p for X with x(p ) =, where x is the normal coordinate at p with respect to the metric ĥ on the boundary Mn. Theorem 3.5. Suppose that U is a nonnegative solution to (3.1) in X n1. Then, for sufficiently small r, if U(q ) = for q Γ r \Γ 1 and U > on Γ on the boundary r 1 Mn, r then y a y U q >. (3.19)

14 Proof. First we assume that γ [1/,1), i.e., a ( 1,]. We consider a positive function W = y a (y Ay )(e B x e Br ). (3.) To calculate div(y a W) in the metric ḡ we first calculate from Proposition. that ḡ = (1α 1 y)dy (1α y)ĥo(y) for some constants α 1,α and detḡ = detĥ(1α 3y)o(y), for some constant α 3. Then div(y a W) = I 1 I I 3 I 4, where I 1 = 1 ) detḡ y( detḡ (ḡ ) yy ((1 a)( a)ya)(e B x e Br ) = (α 4 ( a)ao(1))(e B x e Br ), 1 ( ) I = detḡ x k detḡ (ḡ ) ky ((1 a)( a)ya)(e B x e Br ) for some constant α 4, I 3 = = o(1)(e B x e Br )o(y)be Br, 1 ) detḡ y( detḡ (ḡ ) yk (y y A) x k(e B x e Br ) = o(y)be Br, and I 4 = y y A ( ) detḡ x k detḡ (ḡ ) kj x j(e B x e Br ) Thus = y y A detḡ x k ( detḡ (ḡ ) kj ( x j r Be Br)) = yb e Br o(y)b e Br yb o(r )e Br o(y)be Br. div(y a W) = (α 4 ( a)ao(1))(e B x e Br ) (B o(1)b)ye Br. We remark here that all constants α s can be explicit, but it would not be any more use. Take r sufficiently small and A and B sufficiently large so that provided that a. Now we know div(y a W) div(y a (U ǫw)) ( ) in Γ r \Γ 1 (,r ) for all ǫ >, and moreover r U ǫw {( ) } on Γ r \Γ 1 (,r ), provided we choose ǫ appropriately small. Therefore, due to r the maximum principle we know that U ǫw >.

15 ( ) in Γ r \Γ 1 (,r ). Thus, when U(x(q ),) =, we have r which implies y a y (U ǫw) (x(q),), y a y U (x(q),) ǫy a y W (x(q),) = ǫ(1 a)(e B x(q) e Br ) >, as desired. When a (,1), or equivalently, γ (, 1 ), we instead use the function W = y a (y Ay a )(e B x e Br ). Then a similar calculation will prove that the conclusion still holds. Positivity of solutions for (3.1) is now clear: Corollary 3.6. Suppose that U C (X) C( X) is a nonnegative solution to the equation { div(y a U) = in (X, ḡ ), y a y U = F(U) on M, where F() =. Then U > on X unless U. Proof. First, U > in X, and U is not identically zero on the boundary if it is not identically zero on X. Then, on the boundary, the set where U is positive is nonempty and open. Hence, if the set where U vanishes is not empty, then, for any small number r, there always exist points p and q as given in the assumptions of Theorem 3.5. Thus we would arrive at the contradiction from Theorem The γ-yamabe problem Now we are ready to set up the fractional Yamabe problem for γ (,1). On the conformal infinity (M n, [ĥ]) of an asymptotically hyperbolic manifold (Xn1, g ), we consider a scalefree functional on metrics in the class [ĥ] given by I γ [ĥ] = M Qĥγ dv ĥ ( M dv ĥ )n γ n Or, if we set a base metric ĥ and write a conformal metric then I γ [w,ĥ] = ĥ w = w 4 n γ ĥ, M wpĥγ (w)dv ĥ ) ( M w dvĥ. (4.1) (4.) where = n n γ. We will call I γ the γ-yamabe functional. The γ-yamabe problem is to find a metric in the conformal class [ĥ] that minimizes the γ-yamabe functional I γ. It is clear that a metric ĥw, where w is a minimizer of I γ [w,ĥ], has a constant fractional scalar curvature Qĥw γ, that is, for some constant c on M. nγ Pĥγ (w) = cw n γ, w >, (4.3)

16 This suggests that we define the γ-yamabe constant Λ γ (M,[ĥ]) = inf { I γ [h] : h [ĥ] }. (4.4) It is then apparent that Λ γ (M,[ĥ]) is an invariant on the conformal class [ĥ] when g is fixed. In the mean time, based on Proposition.1, we set Iγ [U,ḡ] = d γ X ρa U ḡ dv ḡ X E(ρ) U dvḡ ( ), (4.5) M U dvĥ or similarly, using Proposition., we may set Iγ[U,ḡ ] = d γ X ya U ḡ dv ḡ M Qĥγ U dvĥ ( ). (4.6) M U dvĥ It is obvious that it is equivalent to solve the minimizing problems for I γ and Iγ. But a very pleasant surprising is that this immediately tells us that Λ γ (X, [ĥ]) = inf{ I γ [U,ḡ] : U W1, (X,y a ) } (4.7) (please see the definitions and discussions of the weighted Sobolev spaces in Section 5). Note that one has that I γ[ U ] I γ[u], to handle positivity issues. Therefore we have Lemma 4.1. Suppose that U is a minimizer of the functional Iγ [,ḡ] in the weighted Sobolev space W 1, (X,y a ) with M TU dvĥ = 1. Then its trace w = TU H γ (M) solves the equation Pĥγ (w) = Λ nγ γ(x, [ĥ])w n γ. To resolve the γ-yamabe problem is to verify I γ has a minimizer w, which is positive and smooth. But before launching our resolution to the γ-yamabe problem we are first due to discuss the sign of the γ-yamabe constant. These statements are familiar and easy ones for the Yamabe problem but not so easy at all for the γ-yamabe problem, where the conformal fractional Laplacians are just pseudo-differential operators. One knows that eigenvalues and eigenfunctions of the conformal fractional Laplacians are even more difficult to study than the differential operators. There are some affirmative results analogous to the conformal Laplacian proven in [8] when the Yamabe constant of the conformal infinity is assumed to be positive. Here we will take the advantage of our Hopf s Lemma and the interpretation of the conformal fractional Laplacians through extensions provided in Proposition.. For each γ (,1) we know that each conformal fractional Laplacian is self-adjoint (cf. [7], []). Hence we may look for the first eigenvalue λ 1 by minimizing the quotient M wpĥγ wdvĥ. (4.8) M w dvĥ Moreover, again in the light of Proposition., it is equivalent to minimizing d γ X ya U ḡ dv ḡ M Qĥγ U dvĥ. (4.9) M U dvĥ We arrive at the eigenvalue equation: Pĥγ w = λ 1w, on M. Or, equivalently, div(y a U) = in (X, ḡ ), d γ lim y y a y U QĥγU = λ 1 U on M, As a consequence of Proposition. and Theorem 3.5 we have: (4.1)

17 Theorem 4.. Suppose that (X n1, g ) is an asymptotically hyperbolic manifold. For each γ (,1) there is a smooth, positive first eigenfunction for Pĥγ and the first eigenspace is of dimension one, provided H = when γ ( 1,1). Proof. We use the variational characterization (4.9) of the first eigenvalue. We first observe that one may always assume there is a nonnegative minimizer for (4.9). Then regularity and the maximum principle in Section 3 insure that such a first eigenfunction is smooth and positive. To show that the first eigenspace is of dimension 1, we suppose that φ and ψ are positive first eigenfunctions for Pĥγ. Then ψ Pĥφ γ φ = φ nγ n γ P ĥ γ ψ = λ 1 φ nγ = (φ nγ n γ P ĥ γ φ) ψ φ ψ = Qĥφ γ φ, where ĥφ = φ 4 n γ ĥ. That is, there is a function U satisfying lim n γ ψ div(y a φ U) = in (X, ḡ φ), y φ ya φ U U = on M, y φ and U = ψ φ on M, where y φ and ḡ φ are associated with ĥφ as y and ḡ are associated with ĥ in Proposition. respectively. eplace U by U U m for U m = min X U and apply Theorem 3.5 and Corollary 3.6 to conclude that U has to be a constant. Consequently, we get the following. Corollary 4.3. Suppose that (X n1, g ) is an asymptotically hyperbolic manifold. Assume that γ (,1) and that H = when γ ( 1,1). Then there are three mutually exclusive possibilities for the conformal infinity (M n, [ĥ]): 1. The first eigenvalue of Pĥγ is positive, the γ-yamabe constant is positive, and M admits a metric in [ĥ] that has pointwise positive fractional scalar curvature.. The first eigenvalue of Pĥγ is negative, the γ-yamabe constant is negative, and M admits a metric in [ĥ] that has pointwise negative fractional scalar curvature. 3. The first eigenvalue of Pĥγ is zero, the γ-yamabe constant is zero, and M admits a metric in [ĥ] that has vanishing fractional scalar curvature. Proof. First of all it is obvious that the sign of the first eigenvalue of the conformal fractional Laplacian Pĥγ does not change within the conformal class due to the conformal covariance property of the conformal fractional Laplacian. The three possibilities are distinguished by the sign of the first eigenvalue λ 1 of the conformal fractional Laplacian Pĥγ. Because, if φ is the positive first eigenfunction of Pĥγ, then where ĥφ = φ 4 n γ ĥ. Qĥφ γ = λĥ1φ 4γ n γ

18 5 Weighted Sobolev trace inequalities Let us continue in the setting provided by Proposition.. On the compact manifold M n, for γ (,1), we recall the fractional order Sobolev space H γ (M), with its usual norm w H γ (M) := w L (M) w( ĥ) γ wdvĥ. An equivalent norm on this space is M w H γ (M) := A w L (M) M wpĥγ wdv ĥ, for some appropriately large number A, since Pĥγ is an elliptic pseudo-differential operator of order γ with its principal symbol being the same as that of ( ĥ) γ. Note that in n, this Sobolev norm can be easily written in terms of the Fourier transform as w H γ ( n ) = n (1 ξ ) γ ŵ (ξ)dξ. (5.1) We would also like to recall the definition of the weighted Sobolev spaces. For γ (,1) and a = 1 γ, consider the norm U W 1, (X,y a ) = y a U ḡ dv ḡ y a U dvḡ. The following is then known. Lemma 5.1. There exists a unique linear bounded operator X T : W 1, (X,y a ) H γ (M) such that TU = U M for all U C ( X), which is called the trace operator. Lemma 5.1 was explored by Nekvinda [37] in the case when X is a subset of n1 and M n a piece of its boundary; see also [3]. It then takes some standard argument to derive the Lemma 5.1 from, for instance, [37]. The classical Sobolev trace inequality on Euclidean space is well known (see, for instance, Escobar [15]), and reads: ( n Tu n n 1 dx ) n 1 n C(n) ( n1 X u dxdy )1 (5.) where the constant C(n) is sharp and the equality case is completely characterized. This corresponds to a = for our cases. The same result is true for any other real a ( 1,1). Indeed there are general Weighted Sobolev trace inequalities. Let us first recall the well known fractional Sobolev inequalities. They were considered first in the remarkable paper by Lieb [31] (see also the more recent [1], [13], or the survey [14]): Lemma 5.. Let < γ < n/, = n n γ. Then, for all w Hγ ( n ) we have w L ( n ) S(n,γ) ( )γ w Hγ ( n ) = S(n,γ) n w( ) γ wdx, (5.3) where ) [ S(n,γ) = γ π γγ( n γ Γ ( Γ(n) ) nγ Γ ( ) n We have equality in (5.3) if and only if ( µ w(x) = c x x µ for c, µ > and x n fixed. ]γ n )n γ = Γ( ) n γ Γ ( nγ) vol(s n ) γ n., x n,

19 Note that we may interpret the above inequality as a calculation of the best γ-yamabe constant on the standard sphere as the conformal infinity of the Hyperbolic space. Namely, if g c is the standard round metric on the unit sphere, w L (S n ) S(n,γ) S n wp gc γ wdv g c. (5.4) Such an inequality for the sphere case was also considered independently by Beckner [5], Branson [6], and Morpurgo [35], in the setting of interwining operators. Indeed, we have the following explicit expression for P Sn P Sn It is clear from (5.4) that γ : γ = Γ( ) B γ 1 Γ ( ), B γ 1 where B := S n ( ) n 1. Λ γ (S n,[g c ]) = 1 S(n,γ). (5.5) Sobolev trace inequalities can be obtained by the composition of the trace theorem and the Sobolev embedding theorem above. There have been some related works that deal with these types of energy inequalities, for instance, Nekvinda [37], González [3], and Cabré-Cinti [7]. In particular, in the light of the work of Caffarelli and Silvestre [1] and Lemma 5., we easily see the more general form of (5.) as follows: Corollary 5.3. Let w H γ ( n ), γ (,1), a = 1 γ, and U W 1, ( n1,ya ) with trace TU = w. Then w L ( n ) S(n,γ) y a U dxdy, (5.6) n1 where Equality holds if and only if w(x) = c ( S(n,γ) := d γs(n,γ). (5.7) µ x x µ )n γ, x n, for c, µ > and x n fixed, and U is its Poisson extension of w as given in (.13). In the following lines we take a closer look at the extremal functions that attain the best constant in the inequality above. On n we fix w µ (x) := ( µ x µ )n γ these correspond to the conformal diffeomorphisms of the sphere. We set as given in (.13). Then we have the equality w µ L ( n ) = S(n,γ) It is clear that w µ (x) = 1 ( ) x w µ n γ 1 µ, (5.8) U µ = K γ x w µ (5.9) n1 y a U µ dxdy., and U µ (x,y) = 1 µ n γ U 1 ( x µ, y µ ). (5.1)

20 Moreover, U µ is the (unique) solution of the problem div(y a U µ ) = in n1, lim y a y U µ = c n,γ (w µ ) nγ n γ on n, y (5.11) On the other hand, if we multiply equation (5.11) by U µ and integrate by parts, y a U µ dxdy = c n,γ (w µ ) dx. (5.1) n n1 Now we compare (5.1) with (5.6). Using (5.5) we arrive at Λ(S n,[g c ]) = c n,γ d γ [ n (w µ ) ] γ n dx. (5.13) Before the end of this section we calculate the general upper bound of the γ-yamabe constants. Indeed there is a complete analogue to the case of the usual Yamabe problem (cf. [3], [3]). Namely, the following. Proposition 5.4. Let γ (,1). Then Λ γ (M,[ĥ]) Λ γ(s n,[g c ]). Proof. First of all we will instead use the functional (4.6) to estimate the γ-yamabe constant for a good reason. The approach is rather the standard method of gluing a bubble (5.8) to the manifold M (see, for instance, [3], Lemma 3.4). For any fixed ǫ >, let B ǫ be the ball of radius ǫ centered at the origin in n1 and B ǫ be the half ball of radius ǫ in n1. Choose a smooth radial cutoff function η, η 1 supported on B ǫ, and satisfying η 1 on B ǫ. Then, consider the function V = ηu µ with its trace v = ηw µ on n. We have that n1 y a V dxdy (1ǫ) n1 y a U µ dxdy C(ǫ) B ǫ \B ǫ U µdxdy. (5.14) Note that w µ = O(µ n γ x γ n ) in the annulus ǫ x ǫ and U µ is O(µ n γ ) in the annulus B ǫ \B ǫ. This allows to estimate the second term in right hand side of (5.14) by O ( µ n γ) as µ, for ǫ fixed. For the first term in the right hand side of (5.14) we first use the fact that w µ attains the best constant in the Sobolev inequality, so ( S(n, γ) y a U µ dxdy = n1 w n µ dx ) ( v n ) dx O(µ n ). (5.15) Now we need to transplant the function V to the manifold ( X,ḡ ). Fix a point on the boundary M and use normal coordinates {x 1,...,x n,y} around it, in a half ball B ǫ where V is supported. Two things must be modified: when ǫ, V ḡ = V (1O(ǫ)), and so that dvḡ = (1O(ǫ))dxdy, I ǫ,µ := d γ y a V ḡ dv ḡ B ǫ Qĥγv dvĥ x ǫ ( (1O(ǫ)) y a V dxdy C B ǫ x <ǫ ) v dx.

21 It is easily seen that x <ǫ w µdx = o(1). This is a small computation that can be found in Lemma 3.5 of [3]. Then, from (5.15), fixing ǫ small and then µ small, we can get that ( ) 1 I ǫ,µ (1Cǫ) S(n,γ) v L (M) Cµ which implies Λ γ (M,[ĥ]) 1 S(n,γ) = Λ γ(s n,[g c ]). We end this section by remarking that, although most of the results mentioned here were already known in different contexts, it is certainly very interesting to put all the analysis and geometry together in the context of conformal fractional Laplacians and the associated γ-yamabe problems in a way that is analogous to what has been done on the subject of the Yamabe problem, which becomes fundamental to the development of geometric analysis. 6 Subcritical approximations In this section we take a well known subcritical approximation method to solve the γ-yamabe problem and prove Theorem 1.4. There does not seem to be any more difficulty than usual after our discussions in previous sections. But, for the convenience of the readers, we present a brief sketch of the proof. Similar to the case of the usual Yamabe problem we consider the following subcritical approximations to the functionals I γ and I γ respectively. Set and I β [w] = ( M wpĥγ wdvĥ ) β M wβ dvĥ Iβ[U] = d γ X ya U ḡ dvḡ M Qĥγ U dvĥ (. M Uβ dvĥ for β [, ), where = n n γ and γ (,1). These are subcritical problems and can be solved through standard variational methods. For clarity we state the following: Proposition 6.1. For each β <, there exists a smooth positive minimizer U β for Iβ [U] in W 1, (X,y a ), which satisfies the equations div(y a U β ) = in (X,ḡ ), ) β d γ lim y y a U β QĥγU β = c β U β 1 β on M, where the derivatives are taken with respect to the metric ḡ in X and c β = I β [U β] = mini β. And the boundary value w β of U β, which is a positive smooth minimizer for I β [w] in H γ (M), satisfies Pĥγw β = c β w β 1 β. Using a similar argument as in the proof of Lemma 4.3 in [3] (see also [3]) we have the following. Lemma 6.. If vol(m,ĥ) = 1, then c β is non-increasing as a function of β [, ]; and if Λ γ (M,[ĥ]), then c β is continuous from the left at β =.

22 We now start the proof of Theorem 1.4. eaders are referred to [16], [3], [4] for more details. Instead of applying the standard Sobolev embedding in the Yamabe problem we apply the weighted trace ones discussed in the previous section. To ensure that U β as β produces a minimizer for the γ-yamabe problem, we want to establish the a priori estimates for U β. In the light of the discussions in Section 3, we only need to have a uniform L bound for w β. We will establish the L bound for w β by the so-called blow-up method. Otherwise, assume there exist sequences β k, w k := w βk and U k := U βk, x k M such that w k (x k ) = max M {w k } = m k and x k x M as k. Take a normal coordinate system centered at x, and rescale with the boundary value V k (x,y) = m 1 k U k(δ k xx k,δ k y), v k (x) = m 1 k w k(δ k xx k ), 1 β k γ where δ k = mk. Then V k is defined in a half ball of radius k = 1 x k δ k and is a solution of div(ρ a V k ) = in B k, d γ lim y ya y V k (Qĥγ ) kv k = c k v β 1 k on B k, with respect to the metric ḡ (δ k xx k,δ k y), where (Qĥγ) k = δ 1 a k Qĥγ(δ k xx k ). Due to, for example, C,α a priori estimates for the rescaled solutions V k, to extract a subsequence if necessary, we have V k V in C loc,α. Moreoverthe metrics ḡ (δ k xx k,δ k y) converge to the Euclidean metric. Hence V is a non-trivial, non-negative solution of Let v = TV. It is easily seen that div(y a V ) = in n1, (6.1) d γ lim y a nγ n γ y V = c V y on n,. (6.) v n (x)dx 1. (6.3) Theorem 3.5 and Corollary 3.6 then assure that V > on n1. Therefore we can obtain y a V dxdy = c d γ (x)dx. (6.4) n1 n v It is then obvious that c >, that is, c = Λ γ (M,[ĥ]) in the light of Lemma 6.. Moreover, by the trace inequalities from Lemma 5.3, we have ( S(n,γ) y a V dxdy. (6.5) n v (x)dx ) Then (6.3), (6.4) and (6.5), together with the definition of Λ γ (S n,[g c ]) in (5.5) contradict the initial hypothesis (1.6). Once we have a uniform L estimate, by the regularity theorems in Section 3 we may extract a subsequence if necessary and pass to a limit U, whose boundary value w satisfies n1 Pĥγw = Λw 1, I γ [w ] = Λ, Λ = limc β. (6.6) Theorem 3.5 and Corollary 3.6 also ensure that w > on M. It remains to check that Λ = Λ γ (M,[ĥ]). However, this is a direct consequence of Lemma 6. when Λγ(M,[ĥ]). Meanwhile it is easily seen that by the definition of the γ-yamabe constants and (6.6) that Λ can not be less than Λ γ (M,[ĥ]). Hence it is also implied that Λ = Λγ(M,[ĥ]) by Lemma 6. when Λ γ (M,[ĥ]) <. Thus, in any case, w is a minimizer of I γ, as desired.

23 7 A sufficient condition In this section we give the proof of Theorem 1.5, which provides a sufficient condition for the resolution of the γ-yamabe problem. Here the precise structure of the metric will play a crucial role since a careful computation of the asymptotics is required, following the calculation in [16]. The section is divided into two parts: the first contains the necessary estimates on the Euclidean case, while in the second we go back to the geometry setting and finish the proof of the theorem. 7.1 Some preliminary results on n1 Here we consider the divergence equation (.11) on n1, as understood in [1], [3]. The main point is that by using the Fourier transform, a solution to this problem can be written in terms on its trace value on n and the well known Bessel functions. Indeed, let U be a solution of { div(y a U) = in n1, U(x,) = w on n (7.1) {}, or equivalently, U = K γ x w, where K γ is the Poisson kernel as given in (.13). The main idea is to reduce (7.1) to an ODE by taking Fourier transform in x. We obtain ξ û(ξ,y) a yûy(ξ,y)û yy (ξ,y) =, Û(ξ,) = ŵ(ξ), that is an ODE for each fixed value of ξ. On the other hand, consider the solution ϕ : [, ) of the problem ϕ(y) a y ϕ y(y)ϕ yy (y) =, (7.) subject to the conditions ϕ() = 1 and lim ϕ(t) =. This is a Bessel function and its t properties are summarized in Lemma 7.1. Then we have that Û(ξ,y) = ŵ(ξ)ϕ( ξ y). (7.3) For a review of Bessel functions (see, for instance, Lemma 5.1 in [3], or section in [1]): Lemma 7.1. Consider the following ODE in the variable y > : ϕ(y) a y ϕ y(y)ϕ yy (y) =, with boundary conditions ϕ() = 1, ϕ( ) =. Its solution can be written in terms of Bessel functions: ϕ(y) = c 1 y γ K γ (y), where K γ is the modified Bessel function of the second kind that has asymptotic behavior K γ (y) Γ(γ) K γ (y) ( ) γ, when y, y π y e y, when y, for a constant c 1 = 1 γ Γ(γ).

24 Now we are ready to prove the main technical lemmas in the proof of Theorem 1.5. More precisely, we will explicitly compute several energy terms through Fourier transforms, thanks to expression(7.3). Such precise computation is needed in order to obtain the exact value of the constant(1.8). Forthe restofthe section, wedenote U = ( x1 U)...( xn U) ( y U), and x U = ( x1 U)...( xn U). Lemma 7.. Given w H γ ( n ), let U = K γ w defined on n1. Then A 1 (w) := y a U dxdy = d 1 ŵ(ξ) ξ (γ 1) dξ, (7.4) n where A (w) := A 3 (w) := n1 n1 n1 y a x U dxdy = d n ŵ(ξ) ξ (γ 1) dξ, (7.5) y a U dxdy = d 3 n ŵ(ξ) ξ (γ 1) dξ, (7.6) d = a3 d 1, d 3 = 1 6 a1 d 1. Proof. We write A i := A i (w), i = 1,,3, for simplicity. Note that the integrals in the right hand side of (7.4), (7.5), (7.6) are finite because w H γ ( n ) H γ 1 ( n ), and because of the definition of the Sobolev norm (5.1). Thanks to (7.3) we can easily compute, using the properties of the Fourier transform, A 1 : = y a U dxdy = y a( x U y U ) dxdy n n = y a( ξ Û y Û ) dydξ n = y a ŵ(ξ) ξ ( ϕ( ξ y) ϕ ( ξ y) ) dydξ n = ŵ(ξ) ξ 1 a t a( ϕ(t) ϕ (t) ) dtdξ n = d 1 n ŵ(ξ) ξ 1 a dξ (7.7) for a constant Similarly, for d 1 := t a( ϕ(t) ϕ (t) ) dt. (7.8) A : = y a x U dxdy = y a ξ Û dydξ n n = y a ŵ(ξ) ξ ϕ( ξ y) dydξ n = ŵ(ξ) ξ 1 a t a ϕ(t) dtdξ n = d n ŵ(ξ) ξ 1 a dξ d := t a ϕ(t) dt. (7.9)

25 And finally, A 3 : = for n1 y a U dxdy = = ŵ(ξ) ξ 1 a n n y a Û dydξ = n y a ŵ(ξ) ϕ( ξ y) dydξ t a ϕ(t) dtdξ = d 3 n ŵ(ξ) ξ 1 a dξ, d 3 = t a ϕ(t) dt. (7.1) In the next step, we find the relation between the constants d 1,d,d 3. All the integrals will be evaluated between zero and infinity in the following. Multiply (7.) by ϕ t t a3 and integrate by parts: ϕϕ t t a3 a ϕ tt a ϕ tt ϕ t t a3 =. (7.11) In the above formula, we estimate the first term by t a3 ϕϕ t = 1 t a3 ( t ϕ ) = a3 t a ϕ, and the last one by t a3 ϕ tt ϕ t = 1 t a3 t ( ϕ t ) = a3 t a ϕ t, so from (7.11) we obtain (a3) t a ϕ = ( a3) t a ϕ t. Together with (7.8) and (7.9) this gives d 1 = 6 a3 d, as desired. Now, multiply equation (7.) by ϕt a and integrate: t a ϕϕ t a t a1 ϕ t t a ϕ tt ϕ =. (7.1) The third term above is computed as t a ϕ tt ϕ = t a ϕ t (a) t a1 ϕ t ϕ, so (7.1) becomes d 1 = t a1 ϕ t ϕ = (a1) t a ϕ = (a1)d 3. (7.13) This completes the proof of the lemma. In the following, we continue the estimates of the different error terms, although now we only need the asymptotic behavior and not the precise constant. Lemma 7.3. Let w be defined on n and U = K γ x w. Then 1. For each k N, if w H γ k/ ( n ), E k := n1 y ak U dxdy <. (7.14)