1. Introduction
A classical problem arising in geometric analysis consists in prescribing certain geometric quantities on Riemannian manifolds via a conformal change of metric. It dates back to [Reference Berger9, Reference Kazdan36], where the authors proposed the following problem: given a smooth function $K$ defined on compact surface $(M^2,\,g)$
, can $K$
be achieved as the Gaussian curvature of $M$
with respect to some conformal metric $\tilde g=e^ug$
?
Analytically, this reduces to solve the following equation in $u$:

where $\Delta _g$ is the Laplace–Beltrami operator and $k_g$
denotes the Gaussian curvature of $M$
relative to $g$
. Over the last few decades, this equation has received significant attention and it is not possible to give here a comprehensive list of references; a collection of results can be found in [Reference Aubin4].
If $M$ has a boundary, then boundary conditions are in order. Here we consider a non-linear boundary condition corresponding to conformally prescribing the boundary geodesic curvature $H$
, for some given function $H$
defined on $\partial M$
. In this case, we are led to the boundary value problem:

where $\nu$ is the exterior normal vector to $\partial M$
and $h_g$
its initial geodesic curvature. Equation (1.1) has been considered in particular situations; the case of constant $K$
and $H$
has been studied by Brendle in [Reference Brendle10] by means of a parabolic flow. In this situation, some classification results are available when $M$
is an annulus or the half-space, see [Reference Jiménez35, Reference Li and Zhu39, Reference Zhang47]. The case of non-constant curvatures was first addressed by Cherrier in [Reference Cherrier17], but his results are obstructed by the presence of Lagrange multipliers. Recently, the general case with $K<0$
in surfaces topologically different from the disk has been studied in [Reference López-Soriano, Malchiodi and Ruiz41], and a blow-up analysis has been performed.
Generally speaking, the case of a disk is challenging due to the non-compact nature of the conformal map group acting on it, similarly to the Nirenberg problem on ${\mathbb {S}}^2$. The problem becomes

Integrating (1.2) and applying the Gauss–Bonnet theorem, we obtain

from which we see that $K$ or $H$
needs to be somewhere positive. Some partial results are available for the case in which one of the curvatures is zero, see [Reference Chang and Liu13, Reference Chang and Yang14, Reference Da Lio, Martinazzi and Riviére23, Reference Guo and Liu31, Reference Li and Liu38, Reference Liu and Huang40]. However, up to our knowledge, there are few results available for the case of non-constant functions $K$
and $H$
. In [Reference Cruz-Blázquez and Ruiz20], the problem is posed in a new variational setting and existence of solutions in the form of global minimizers is obtained for non-negative and symmetric curvatures. Existence results for not necessarily symmetric, non-negative curvatures are found in [Reference Ruiz45] via a Leray–Schauder degree argument. In [Reference Battaglia, Medina and Pistoia6], the authors construct blowing-up solutions for (1.2) under certain non-degeneracy assumptions on $K$
and $H$
using a Ljapunov–Schmidt reduction.
Concerning the blow-up behaviour of sequences of solutions, a rather exhaustive study is given in [Reference Jevnikar, López-Soriano, Medina and Ruiz34] (see also [Reference Da Lio, Martinazzi and Riviére23, Reference Guo and Liu31]). In particular, it is shown that if $K<0$, then problem (1.2) only admits blow-ups at boundary points where the scaling-invariant function ${\mathfrak {D}}_2:{\mathbb {S}}^1\to \mathbb {R}$
defined as

is greater or equal than one (see [Reference Jevnikar, López-Soriano, Medina and Ruiz34, Theorem 1.1] and [Reference López-Soriano, Malchiodi and Ruiz41, Theorem 1.4]).
The natural analogue of this question in higher dimensions is the prescription of the scalar curvature of a manifold and the mean curvature of the boundary, and has received more attention.
More precisely, if $(M^n,\,g)$ is a Riemannian manifold of dimension $n\ge 3$
with boundary, and $K:\overline M\to \mathbb {R},\,H:\partial M\to \mathbb {R}$
are given smooth functions, it consists of finding positive solution for the boundary problem

Here $k_g$ and $h_g$
denote the scalar and boundary mean curvatures of $M$
with respect to $g$
. If $u>0$
solves (1.4), then the metric $\tilde g=u^\frac 4{n-2}g$
satisfies $k_{\tilde g}=K$
and $h_{\tilde g}=H$
.
In the literature, we can find many partial results for this equation. The case of prescribing a scalar flat metric with constant boundary mean curvature is known as the Escobar problem, in strong analogy with the Yamabe problem. Its study was initiated by Escobar in [Reference Escobar25, Reference Escobar26, Reference Escobar and García28], with later contributions in [Reference Almaraz2, Reference Marques42–Reference Mayer and Ndiaye44]. Different settings with constant curvatures are considered in [Reference Brendle and Chen11, Reference Chen, Ruan and Sun16, Reference Escobar27, Reference Han and Li32, Reference Han and Li33]. Some results are available for the case of non-constant functions when one of them is equal to zero. Existence results for the scalar flat problem are given in [Reference Abdelhedi, Chtioui and Ahmedou1, Reference Chang, Xu and Yang12, Reference Djadli, Malchiodi and Ahmedou24, Reference Xu and Zhang46], while the works [Reference Ben Ayed, El Mehdi and Ahmedou7, Reference Ben Ayed, El Mehdi and Ahmedou8, Reference Li37] concern the case with minimal boundaries.
On the other hand, the problem with non-constant functions $K$ and $H$
has received comparatively little study. In this regard, we highlight [Reference Ambrosetti, Li and Malchiodi3], which contains perturbative results about nearly constant positive curvature functions on the unit ball of $\mathbb {R}^n$
.
The case of non-constant $K>0$ and $H$
of arbitrary sign was also considered in [Reference Djadli, Malchiodi and Ahmedou22] in the half sphere of $\mathbb {R}^3$
, and a blow-up analysis was carried out. As for negative $K$
, in [Reference Chen, Ho and Sun15] the authors study equation (1.4) with $K<0$
and $H<0$
by means of a geometric flow, in the spirit of [Reference Brendle10], but solutions are obtained up to Lagrange multipliers. Finally, in the recent work [Reference Cruz-Blázquez, Malchiodi and Ruiz19], the case with non-constant functions $K<0$
and $H$
of arbitrary sign is treated on manifolds of non-positive Yamabe invariant. Similarly to the two-dimensional case, it is shown that the nature of the problem changes greatly depending on whether the function ${\mathfrak {D}}_n:\partial M\to \mathbb {R}$
given by

is less than one over the entire boundary or not. When ${\mathfrak {D}}_n<1$ the energy functional becomes coercive and a global minimizer can be found. However, if ${\mathfrak {D}}_n\ge 1$
somewhere on $\partial M$
, a min–max argument and a careful blow-up analysis are needed to recover the existence of solutions, although only in dimension three.
In this paper, we will focus on the following perturbative version of (1.2) and (1.4):

and if $n\ge 3$

where ${\mathfrak {D}}_n>1$ is defined in (1.3) and (1.5), $K:\mathbb {B}^n\to \mathbb {R},\,H:\mathbb {S}^{n-1}\to \mathbb {R}$
are smooth with bounded derivatives and the parameter $\varepsilon \in \mathbb {R}$
small.
Our main result for problem ($P^2_\varepsilon$) reads as follows:
Theorem 1.1 Assume ${\mathfrak {D}}_2\ne \frac 2{\sqrt 3},$ let $\psi :{\mathbb {S}}^1\to \mathbb {R}$
be defined by

and $\Phi _1:{\mathbb {S}}^1\to \mathbb {R}$ be defined by

and $\Phi _m$ be defined as in definition 2.1. If one of the following holds true:
(1) For any global maximum $\xi$
of $\psi$
there exists $m=m(\xi )\ge 1$
such that $\Phi _j(\xi )=0>\Phi _m(\xi )$
for any $j< m$
;
(2) For any global minimum $\xi$
of $\psi$
there exists $m=m(\xi )\ge 1$
such that $\Phi _j(\xi )=0<\Phi _m(\xi )$
for any $j< m$
;
(3) For any critical point $\xi$
of $\psi$
there exists $m=m(\xi )\ge 1$
such that $\Phi _j(\xi )=0\ne \Phi _m(\xi )$
for any $j< m$
, $\psi$
is Morse and
\begin{align*} \sum_{\{\xi:\nabla\psi(\xi)=0,\,\Phi_m(\xi)<0\}}({-}1)^{\operatorname{ind}_\xi\nabla\psi}\ne1; \end{align*}
then, problem ($P^2_\varepsilon$) has a solution for $|\varepsilon |$
small enough.
Our main result for problem ($P^n_\varepsilon$) reads as follows:
Theorem 1.2 Let $\psi :{\mathbb {S}^{n-1}}\to \mathbb {R}$ be defined by

with $\mathtt a({\mathfrak {D}}_n),\,\mathtt b({\mathfrak {D}}_n)$ as in (2.3), $\Phi _1:{\mathbb {S}^{n-1}}\to \mathbb {R}$
be defined by

and $\Phi _m:{\mathbb {S}^{n-1}}\to \mathbb {R}$ be defined as in definition 2.1. If one of the following holds true:
(1) For any global maximum $\xi$
of $\psi$
there exists $m=m(\xi )\ge 1$
such that $\Phi _j(\xi )=0>\Phi _m(\xi )$
for any $j< m$
;
(2) For any global minimum $\xi$
of $\psi$
there exists $m=m(\xi )\ge 1$
such that $\Phi _j(\xi )=0<\Phi _m(\xi )$
for any $j< m$
;
(3) For any critical point $\xi$
of $\psi$
there exists $m=m(\xi )\ge 1$
such that $\Phi _j(\xi )=0\ne \Phi _m(\xi )$
for any $j< m$
, $\psi$
is Morse and
\begin{align*} \sum_{\{\xi:\nabla\psi(\xi)=0,\,\Phi_m(\xi)<0\}}({-}1)^{\operatorname{ind}_\xi\nabla\psi}\ne1; \end{align*}
then, problem ($P^2_\varepsilon$) has a solution for $|\varepsilon |$
small enough.
Problems ($P^2_\varepsilon$) and ($P^n_\varepsilon$
) share many similarities, not only for their geometric importance, but also from an analytic point of view.
In fact, they both have critical terms in the interior and in the boundary non-linearities: exponential non-linearities in ($P^2_\varepsilon$) are critical in view of the Moser–Trudinger inequalities, whereas in ($P^n_\varepsilon$
) we have the critical Sobolev exponent and the critical trace exponent

Moreover, since we are prescribing a negative curvature in the interior and a positive curvature in the boundary, the two non-linear terms have different signs and are therefore in competition.
Theorem 1.1 seems to be the first result of prescribing both nearly constant curvatures on a disk. Similar results were recently obtained in [Reference Battaglia, Cozzi, Fernández and Pistoia5] in the case of zero curvature in the interior and in [Reference Grossi and Prashanth30] for the sphere. Theorem 1.2 is the counterpart of the result obtained in [Reference Ambrosetti, Li and Malchiodi3], where the authors perturb the positive constant curvature on the unit ball of $\mathbb {R}^n$.
We also provide higher-order expansions of the reduced energy functional, which permits to consider also some cases of degenerate critical points. This is the case when the functionals $\Phi _m(\xi )$ play a role in theorems 1.1 and 1.2.
Such expansions require sharper estimates (see proposition 3.4 and appendix A) and both derivatives of $K,\,H$ and non-local terms appear. In particular, non-local terms are present only if the order of the expansion is high enough, depending on the dimension. At the first order, we only get the fractional Laplacian in the two-dimensional case, which is why $\Phi _1(\xi )$
is defined differently in theorems 1.1 and 1.2.
The definition of $\Phi _m$ for $m\ge 2$
is rather involved and it is therefore postponed to definition 2.1.
Finally, we point out that, in theorem 1.1, $\Phi _1$ can be seen as the normal derivative (up to a constant) of the functional $\psi$
, which can be naturally extended from the circle to the closed disk. More precisely, for $\xi \in \overline {\mathbb {B}^2}$
, we set

therefore, in view of the Dirichlet-to-Neumann characterization of the fractional Laplacian, we have $\Phi _1(\xi )=\partial _\nu \Psi (\xi )$ for any $\xi \in {\mathbb {S}}^1$
.
Quite interestingly, this fact has no higher-dimensional counterpart in theorem 1.2.
The assumption ${\mathfrak {D}}_2\not =\frac 2{\sqrt 3}$ (i.e. $\alpha _{\mathfrak {D}}\not =0$
in proposition 4.3) allows to apply the degree argument to the function which also depend on the extra parameter that only appears in the 2D case (see (2) of proposition 4.3). It would be interesting to understand whether this is a mere technical assumption or not and also whether it has some geometrical meaning.
The plan of the paper is as follows.
In § 2, we introduce some notation and preliminaries which we will use in the following; in § 3, we study the energy functional associated to the system and show some of its crucial properties; in § 4, we apply the Ljapunov–Schmidt finite dimensional reduction; in § 5, we study the existence of critical points to the reduced energy functional; finally, in the appendix, we prove some crucial asymptotic estimates.
2. Notation and preliminaries
We remind that ${\mathbb {B}^n}$ will denote the unit ball of $\mathbb {R}^n$
, for $n\ge 2$
. We consider the well-known inversion map $\mathscr {I}:\mathbb {R}^n_+\to {\mathbb {B}^n}$
defined by

Straightforward computations show that $\mathscr {I}\circ \mathscr {I}=\hbox {Id}$, therefore $\mathscr {I}^{-1}$
has the same expression. For more details about this map, see for instance [Reference Escobar and García28], Section 2.
We point out that, up to the sign of the last coordinate, $\mathscr {I}$ extends the stereographic projection from $\partial \mathbb {R}^n_+$
to ${\mathbb {S}^{n-1}}$
and, in dimension $2$
, it coincides with the Riemann map from the half-plane to the disk. In particular, $\mathscr {I}$
is a conformal map and satisfies

For convenience, we define $\rho :\mathbb {R}^n_+\to \mathbb {R}_+$ by

We point out that $\rho$ satisfies (1.1) or (1.4) for some particular choices of the curvatures. More precisely, in dimension $n=2$

while in dimensions $n\ge 3$

For the reader's convenience, we collect here all the constants that appear in our computations. In the following, we agree that ${\mathfrak {D}}={\mathfrak {D}}_n$.
Definition 2.1 Let $n\ge 2$ and ${\mathfrak {D}}>1$
. We define

In particular, set

We define the functionals $\Phi _j(\xi )$ as follows.
For $j\le n-2$ we set:

For $j=n-1,\,n$ we set:

For $m\ge n$ we set:

The symbol $a\lesssim b$ will be used to mean $a\le cb$
with $c$
independent on the quantities.
We denote as $\partial _\nu$ the (outer) normal derivative of a function at a point on ${\mathbb {S}^{n-1}}$
and as $\Delta _\tau$
the tangential Laplacian.
For a multi-index $\alpha =(\alpha _1,\,\dots,\,\alpha _n)\in \mathbb {N}^n$ we denote:

2.1. Conformal metrics
Throughout this article, we will use the existing conforming equivalence between $\mathbb {R}^n_+$ and ${\mathbb {B}^n}$
via the inversion map (2.1), often without explicitly specifying it. Therefore, it is important to remember the conformal properties of the conformal Laplacian and conformal boundary operator.
If $n\ge 3$ and $\tilde g=\rho ^\frac 4{n-2}g$
is a conformal metric, then the conformal Laplacian and conformal boundary operators, defined by

are conformally invariant in the following sense:

If $n=2$, then the Laplace–Beltrami operator and the normal derivative satisfy the following conformal property: if $\tilde g=e^\rho g$
is a conformal metric, then

The following result establishes the conformal invariance of a certain geometric quantity that will be very much related to our energy functionals.
Lemma 2.2 Let $(M^n,\,g)$ be a compact Riemannian manifold of dimension $n\ge 3$
and $\tilde g=\varphi ^{\frac {4}{n-2}}g$
a conformal metric with $\varphi$
smooth and positive. If we set $\hat f=f\varphi ^{-1},$
then

Proof. We will use the following basic identities:

where ${2^*},\,{2^\sharp }$ are as in (1.6) and the relation between $k_{\tilde g}$
, $k_g$
, $h_{\tilde g}$
and $h_g$
given by (1.4). The first term in the left-hand side of (2.5) can be decomposed using the previous identities:

On the other hand, integrating by parts on $M$ and using (1.4):

2.2. Solutions of the unperturbed problems
By means of the inversion map and the classification results available for $\mathbb {R}^n_+$, we can give an $n$
-dimensional family of solutions for the problems ($P^2_\varepsilon$
) and ($P^n_\varepsilon$
) with $\varepsilon =0$
.
First, we consider the problem in ${\mathbb {B}^2}$:

By [Reference Zhang47], a family of solutions of the problem in the half space

is given by

for $\lambda >0$ and $x_0\in \mathbb {R}$
. Other classification results for solutions to (2.8) are given in [Reference Gálvez and Mira29, Reference Li and Zhu39].
Let us call $\hat U_{x_0,\lambda }=\left (U_{x_0,\lambda }-\rho \right )\circ \mathscr {I}^{-1}$. Taking into account equation (2.2) and the conformal properties of the Laplacian and normal derivative in $\mathbb {R}^2$
, it is clear that

Therefore, a family of solutions for ($P^2_0$) is given by

with

$\lambda >0$ and $x_0\in \mathbb {R}.$
Now, we address the unperturbed problem in ${\mathbb {B}^n}$ for $n\ge 3$
:

Consider $\mathbb {R}^n_+$ with its usual metric, and the problem

The results in [Reference Chipot, Fila and Shafrir18] imply that all solutions of (2.10) have the form

for any $x_0\in \partial \mathbb {R}^n_+$ and $\lambda >0$
.
Then, by (2.4), we can write (2.10) as:

If we call $\hat u=\left (\frac u\rho \right )\circ \mathscr {I}^{-1}$, it is clear that

which is exactly ($P^{n\,'}_0$). Hence, a family of solutions of ($P^{n\,'}_0$
) is given by

with

In view of formulae (2.9) and (2.11) we set:

with $\bar z,\,z_n$ as before, $x_0\in \mathbb {R}^{n-1}$
, $\lambda >0$
, and define

3. Properties of the energy functionals
We define the functionals $J^n_\varepsilon :H^1\left ({\mathbb {B}^n}\right )\to \mathbb {R}$ by


Observe that we can write

with

with $\alpha _n,\,\beta _n$ as in definition 2.1.
Let $V_{x_0,\lambda }$ be given by (2.12). We set

The first term of the energy is constant along our family of solutions:
Proposition 3.1 There exist constants $\mathtt E_{n,{\mathfrak {D}}}$, independent on $\lambda$
and $x_0$
, such that

Proof. Let us study the cases $n=2$ and $n\ge 3$
separately.
When $n=2$, integrating by parts and using ($P^2_0$
) and (2.2), we can see that:

Now,

Finally,

As for the case $n\ge 3$, from lemma 2.2 it follows

Moreover, by a direct change of variables, we obtain

Therefore,

By a change of variables and using the relations in § 2, we can move to $\mathbb {R}^n_+$ and write our function $\Gamma$
in a more suitable way.
Proposition 3.2 It holds

where $\tilde K=K\circ \mathscr {I},\,\tilde H=H\circ \mathscr {I}$.
We are interested in the behaviour of $\Gamma$ at infinity and when $\lambda \to 0$
.
Proposition 3.3 $\lim _{\left \vert x_0\right \vert +\lambda \to +\infty }\Gamma (x_0,\,\lambda )=\psi ((0,\,-1))$
Proof. First, notice that

With that in mind, we fix $\varepsilon >0$ small enough and write

Then, taking limits when $\lambda +\left \vert x_0\right \vert \to +\infty$,

The claim follows from taking limits when $\epsilon \to 0$.
The following result describes the behaviour of $\Gamma$ around $\lambda =0$
. Its proof will be postponed to appendix A.
Proposition 3.4 Define $\psi :{\mathbb {S}^{n-1}}\to \mathbb {R}$ by $\psi (\xi ):=\mathtt a({\mathfrak {D}}_n)K(\xi )-\mathtt b({\mathfrak {D}}_n)H(\xi )$
, and let us write $\xi =\mathscr {I}(x_0)\in {\mathbb {S}^{n-1}}$
. The following expansions hold, for any $m\in \mathbb {N}$
, when $\lambda \ll 1$
:
If $n=2$,

If $n\ge 3$,

Here $\mathtt a({\mathfrak {D}}_n),\,\mathtt b({\mathfrak {D}}_n),\,a_{n,0,1},\,\Phi _j(\xi )$ are given in definition 2.1.
4. The linear theory
In this section, we develop the technicalities of the Ljapunov–Schmidt finite dimensional reduction. Most of the results hereby presented are well-known in the literature of this argument, therefore details of the proofs will be skipped.
4.1. The $2$
-dimensional case
It is known (see [Reference Jevnikar, López-Soriano, Medina and Ruiz34]) that the solutions of the linear problem

are a linear combination of

given $\kappa >0$, set

Arguing as in theorem 3.3 of [Reference Battaglia, Medina and Pistoia6] we can prove that
Proposition 4.1 Fix $p>1$ and $\kappa >0$
. For any $(x_0,\,\lambda )\in \mathtt C_\kappa$
(see (4.1)) and $\mathfrak f\in L^p\left ({\mathbb {B}^2}\right )$
and $\mathfrak g\in L^p\left ({\mathbb {S}}^1\right )$
such that

there exists a unique $\phi \in H^1\left ({\mathbb {B}^2}\right )$ such that

which solves the problem

Furthermore

4.1.1. Rewriting the problem
We look for a solution of ($P^2_\varepsilon$) in the form

where $\phi$ satisfies the orthogonality condition (4.2). We shall rewrite problem ($P^2_\varepsilon$
) as a system

where $c_i$'s are real numbers.
The error that we are paying by using this approximating solution equals to

and the non-linear part is

Here we set

We have the following result
Proposition 4.2 Fix $\kappa >0$. There exists $\varepsilon _\kappa >0$
such that or any $(x_0,\,\lambda )\in \mathtt C_\kappa$
(see (4.1)) there exists a unique $\phi =\phi (\varepsilon,\,x_0,\,\lambda )\in H^1\left ({\mathbb {B}^2}\right )$
and $c_i\in \mathbb {R}$
which solve (4.3). Moreover, $(x_0,\,\lambda )\to \phi (\varepsilon,\,x_0,\,\lambda )$
is a $C^1-$
function and $\|\phi \|\lesssim \varepsilon.$
Proof. The proof is standard and relies on a contraction mapping argument combined with the linear theory developed in proposition 4.1 and the estimates for $p>1$

4.1.2. The reduced energy
Let us consider the energy functional $J_\varepsilon ^2$ defined in (3.1), whose critical points produce solutions of ($P^2_\varepsilon$
). We define the reduced energy

where $\phi$ is given in proposition 4.2.
Proposition 4.3 The following are true:
(1) If $(x_0,\,\lambda )$
is a critical point of $\widetilde J_\varepsilon$
, then $U_{x_0,\lambda }+\phi$
is a solution to ($P^2_\varepsilon$
).
(2) The following expansion holds
\begin{align*} \widetilde J_\varepsilon(t,x_0,\lambda)=\mathtt E_{2,{\mathfrak{D}}}-\varepsilon\left(\alpha_{\mathfrak{D}} t^2+\Gamma(x_0,\lambda)\right)+o(\varepsilon) \end{align*}$C^1$-uniformly in compact sets of $\mathbb {R}\times (0,\,+\infty )\times \mathbb {R}$
.
Here $\mathtt E_{2,{\mathfrak {D}}}$
is a constant independent on $x_0$
, $t$
and $\lambda$
whose expression is given by (3.3), $\Gamma$
is defined in (3.5) and
\begin{align*} \alpha_{\mathfrak{D}}=\pi\left(\frac{\mathfrak{D}}{\sqrt{{\mathfrak{D}}^2-1}}-2\right). \end{align*}
Proof. We use the choice $\tau =t\sqrt \varepsilon$ and the fact that

4.2. The $n$
-dimensional case
Recently, in [Reference Cruz-Blázquez, Pistoia and Vaira21], it has been proved that all the solutions to the linearized problem

are a linear combination of the $n$ functions

Given $\kappa >0$ set

Arguing as in [Reference Cruz-Blázquez, Pistoia and Vaira21] we can prove that
Proposition 4.4 Fix $\kappa >0$. For any $(x_0,\,\lambda )\in \mathtt C_\kappa$
(see (4.5)) and $\mathfrak f\in L^\frac {2n}{n+2}\left (\mathbb {B}^n\right )$
and $\mathfrak g\in L^\frac {2(n-1)}n\left ({\mathbb {S}^{n-1}}\right )$
such that

there exists a unique $\phi \in H^1\left (\mathbb {B}^n\right )$ such that

which solves the problem

Furthermore

4.2.1. Rewriting the problem
We look for a positive solution of ($P^n_\varepsilon$) as

where $\phi$ satisfies (4.6). We rewrite problem ($P^n_\varepsilon$
) as a system

where the $c_i$ are real numbers. Moreover, the error is given by

and the non-linear part is

Here we set

We have the following result:
Proposition 4.5 Fix $\kappa >0$. There exists $\varepsilon _\kappa >0$
such that or any $(x_0,\,\lambda )\in \mathtt C_\kappa$
(see (4.5)) there exists a unique $\phi =\phi (\varepsilon,\,x_0,\,\lambda )\in H^1\left ({\mathbb {B}^2}\right )$
and $c_i\in \mathbb {R}$
which solve (4.7). Moreover, $(x_0,\,\lambda )\to \phi (\varepsilon,\,x_0,\,\lambda )$
is a $C^1-$
function and $\|\phi \|\lesssim \varepsilon.$
Proof. The proof is standard and relies on a contraction mapping argument combined with the linear theory developed in proposition 4.4 and the estimates

4.2.2. The reduced energy
We consider the functional $J^n_\varepsilon$ defined on (3.2). It is easy to see that its critical points are positive solutions to equation ($P^n_\varepsilon$
). Now, we introduce the reduced energy

where $\phi$ is given in proposition 4.5. It is quite standard to prove the following result
Proposition 4.6 The following assertions hold true
(1) If $(x_0,\,\lambda )$
is a critical point of $\widetilde J_\varepsilon$
, then $U_{x_0,\lambda }+\phi$
is a solution to ($P^n_\varepsilon$
).
(2) Moreover, we have the following expansion
\begin{align*} \widetilde J^n_\varepsilon(x_0,\lambda)=\mathtt E_{n,{\mathfrak{D}}}-\varepsilon\Gamma(x_0,\lambda)+o(\varepsilon) \end{align*}$C^1$-uniformly with respect to $(x_0,\,\lambda )$
in compact sets of $(0,\,+\infty )\times \mathbb {R}^{n-1}.$
Here $\mathtt E_{n,{\mathfrak {D}}}$
is a constant independent on $x_0$
and $\lambda$
, given by (3.4), and $\Gamma$
is the function defined on (3.5)
5. Existence of critical points of $\Gamma$
In this section, we are finally able to get critical points of the map $(x_0,\,\lambda )\mapsto \Gamma (x_0,\,\lambda )$, hence solutions to problems ($P^2_\varepsilon$
), ($P^n_\varepsilon$
).
We start with the following abstract result about critical points of maps defined on balls in dependence of the boundary behaviour.
Proposition 5.1 Let $f:{\mathbb {B}^n}\to \mathbb {R}$ be a $C^1$
map satisfying, as $\xi$
goes to ${\mathbb {S}^{n-1}}$
,

for some $f_i:{\mathbb {S}^{n-1}}\to \mathbb {R}$ with $f_0$
of class $C^1$
and some increasing $g_{\frac \xi {|\xi |}}:(0,\,1)\to (0,\,+\infty )$
such that $g_{\frac \xi {|\xi |}}(t)\underset {t\to 0}\to 0$
.
If one of the following holds true:
(1) $f_1(\xi )>0$
at any global maximum $\xi$
of $f_0$
;
(2) $f_1(\xi )<0$
at any global minimum $\xi$
of $f_0$
;
(3) $f_1(\xi )\ne 0$
at any critical point $\xi$
of $f_0$
, $f_0$
is Morse and
\begin{align*} \sum_{\{\xi:\nabla f_0(\xi)=0,\,f_1(\xi)>0\}}({-}1)^{\operatorname{ind}_\xi\nabla f_0}\ne1; \end{align*}
then, $f$ has at least a stable critical point.
Theorems 1.1 and 1.2 will follow without much difficulty from this proposition and proposition 3.4.
Proof of theorems 1.1 and 1.2 We only consider the case of theorem 1.1, since the same arguments also work for theorem 1.2.
Thanks to proposition 4.3, we get a solutions to the problem ($P^2_\varepsilon$) whenever $\frac {\mathfrak {D}}{\sqrt {{\mathfrak {D}}^2-1}}-2\ne 0$
, that is ${\mathfrak {D}}\ne \frac 2{\sqrt 3}$
, and $(x_0,\,\lambda )$
is a stable critical point of $\Gamma$
. After composing with $\mathscr {I}$
, this is equivalent to getting a critical point of the map $f(\xi )=\Gamma \left (\mathscr {I}^{-1}(\xi )\right )$
, which is well-defined and smooth in the whole $\overline {\mathbb {B}^n}$
thanks to proposition 3.3.
In view of proposition 3.4, $f$ satisfies the assumptions of proposition 5.1 with

with $m=m(\xi )$ as in theorem 1.1 (if $m(\xi )$
is not well-defined, as for minima of $\psi$
in case (1) or maxima of $\psi$
in case (2), one can just set $g_\xi (t)=t,\,f_1=-2\pi \Phi _1$
).
Here, we used that $\lambda =\frac {1-|\xi |}{1+\xi _n}+o(1-|\xi |)$ and that

hence the two definitions of $\psi$ given in theorem 1.1 and proposition 3.4 actually coincide.
Since $-2\pi <0$, then the assumptions on $K,\,H$
in theorem 4.3 are equivalent to the ones in proposition 5.1, hence they ensure existence of solutions.
To prove proposition 5.1, we will compute the Leray–Schauder degree of the map $f$.
Proof of proposition 5.1 First of all, $f$ can be extended up to ${\mathbb {S}^{n-1}}$
as $f_0$
. Since $g$
vanishes at $0$
, this extension is continuous.
Assume $(1)$ holds and take an absolute maximum point $\xi _0$
for $f$
on $\overline {\mathbb {B}^n}$
. To get a critical point for $f$
on ${\mathbb {B}^n}$
we suffice to show that $\xi _0\not \in {\mathbb {S}^{n-1}}$
.
If $\xi _0\in {\mathbb {S}^{n-1}}$, we would have $f_1(\xi _0)>0$
, therefore, for $0< t\ll 1$
we would have

contradicting the fact that $\xi _0$ is a maximum point.
If $(2)$ holds, then the same argument shows that the minimum of $f$
on $\overline {\mathbb {B}^n}$
lies in the interior of ${\mathbb {B}^n}$
, therefore it is a critical point of $f$
.
Assume now that $(3)$ holds. We consider the double of $\overline {\mathbb {B}^n}$
, namely the manifold obtained by gluing two copies of ${\mathbb {B}^n}$
along the boundary: $\frac {\overline {\mathbb {B}^n}\times \{0,\,1\}}\sim$
, where $(\xi,\,0)\sim (\xi,\,1)$
for $\xi \in {\mathbb {S}^{n-1}}$
. This manifold is clearly diffeomorphic to $\mathbb {S}^n$
, hence we will identify it as $\mathbb {S}^n$
.
$f$ can be naturally extended to $\tilde f:\mathbb {S}^n\to \mathbb {R}$
as $\tilde f(\xi,\,i)=f(\xi )$
for $i=0,\,1$
. The extension is continuous and, after a suitable rescalement of $g$
close to $0$
, of class $C^1$
(with vanishing normal derivative on the equator). Such a rescalement does not affect the presence of critical points to $\tilde f$
, $f$
and $f|_{{\mathbb {S}^{n-1}}}=f_0$
, which we will now investigate.
We use the Euler–Poincaré formula to compute the Leray–Schauder degree of $\tilde f$, which is a Morse function by assumption:

To deal with the critical points on ${\mathbb {S}^{n-1}}$, we notice that they are exactly the same critical points of $f_0$
, but their index may change, since each can be either a minimum or a maximum in the orthogonal direction; precisely:

and so

Therefore, applying again the Euler–Poincaré formula, this time to $f_0$ on ${\mathbb {S}^{n-1}}$
, we get:

By summing the previous equalities we get:

The latter quantity is non-zero by assumptions, therefore the set of critical points of $f$ on ${\mathbb {B}^n}$
, on which we are taking the first sum, cannot be empty.
Acknowledgements
S. C. acknowledges financial support from the Spanish Ministry of Universities and Next Generation EU funds, through a Margarita Salas grant from the University of Granada, by the FEDER-MINECO Grant PID2021-122122NB-I00 and by J. Andalucia (FQM-116). This work was carried out during his long visit to the University ‘Sapienza Universitá di Roma’, to which he is grateful. The three authors are partially supported by the group GNAMPA of the Istituto Nazionale di Alta Matematica (INdAM). In particular, the first and second authors are funded by the project ‘Fenomeni di blow-up per equazioni nonlineari’, project code CUP_E55F22000270001.
Appendix A. Proof of proposition 3.4
By introducing a rotation in ${\mathbb {B}^n}$ and moving to $\mathbb {R}^n_+$
via $\mathscr {I}$
we can give an expression for $\Gamma$
which is more convenient for our computation.
Lemma A.1 Let $A:{\mathbb {B}^n}\to {\mathbb {B}^n}$ be the rotation corresponding, via the $\mathscr {I}$
, to the translation of $T:x\to x+x_0$
on the half-place, namely $A=\mathscr {I}\circ T\circ \mathscr {I}^{-1}$
. There holds:

where $\tilde K_A=K\circ A\circ \mathscr {I},\,\tilde H_A=H\circ A\circ \mathscr {I}$.
Proof. By doing a change of variables, we observe that

Here we are using that $A$ is a rotation and its very definition, and we set $K_A\!=\!K\circ A, H_A=H\circ A$
. Finally, changing variables twice and using the definitions in § 2:

Proof of proposition 3.4 We start by estimating the boundary term, where some cancellations occur due to symmetry. We expand $\tilde H_A(\lambda y)$ in $\lambda$
up to order $n-2$
:

where we used the formula

and the vanishing, due to symmetry, of integrals of homogeneous polynomials of odd degree or of degree $2j$ which are $j$
-harmonic.
Moreover, in view of the conformal properties of the Laplacian, one has

hence the $j^\mathrm {th}$ term in the expansion equals

In the $j^\mathrm {th}$ order expansion, the remainder is actually $o\left (\lambda ^j\right )$
because we get

In order to deal with higher order terms, we need another argument, since this would get non-converging integrals.
We split the cases $n$ even and $n$
odd.
If $n$ is even, the main order term in the denominator of $I$
is of odd order, hence its integral vanishes. Therefore,

where we used again (A.1); one easily verifies that, due to the behaviours at $0$ at infinity, all the integrals are converging, hence everything is well defined.
After changing variables, the main terms are now

where we used the fact that $\frac 1{|\bar x|^2}=\frac {|z+\xi |^2}{|z-\xi |^2}$ and again (A.2). The small $o$
term contains some new quantities arising when the terms of order $\lambda ^{n-1}$
are transformed into each other.
The smallness of the remainder can be shown similarly as before, here and in the following.
Due to the asymptotic behaviour of both factors, $I'$ can be dealt with similarly as $I$
and one can iterate the argument. In particular, using the series expansion

we get, for any even $m>n$, the following $m^\mathrm {th}$
order term:

In particular, we point out that if $n=2$ this is the main order term in the boundary estimates, and it equals

Let us now consider the case $n$ odd. Here, the first term does not vanish and it gives rise to a logarithmic term. In fact,

iterating, for any odd $m>n$ we get

The argument to estimate the interior terms is similar. We expand $\tilde K(\lambda y)$ up to order $n-1$
, which is the highest power that can be integrated against $P_{x_0,\lambda }^\frac n2$
.

where we wrote the derivation in $\bar x$ and $x_n$
as

and used again cancellation by symmetry and (A.1). In the last step, we used (A.2) (in $\bar x$) and that

In the case $n=2$, since $\int _{\mathbb {R}^2_+}\frac {y_2}{\left (\bar y^2+(y_n+{\mathfrak {D}})^2-1\right )^n}{\rm d}\bar y{\rm d}y_2=\frac \pi 2\left (\sqrt {{\mathfrak {D}}^2-1}-{\mathfrak {D}}\right )$
, putting together with (A.3) we get the first-order expansion

whereas when $n\ge 3$ the first-order expansion contains only the interior term:

As for $I''$, we get local terms involving derivatives of $\tilde K_A$
and non-local terms similarly as before:

In order to iterate and find the next order terms, we again need a series expansion: we get

which in turn comes from

Therefore, for $m>n$ we get:

The proof is now complete, since all the quantities are the same as in definition 2.1.