<figure><img src="image_1.png" /></figure>

HAL
open science

# Principal eigenvalues and eigenfunctions for fully nonlinear equations in punctured balls

Isabeau Birindelli, Françoise Demengel, Fabiana Leoni

► To cite this version:

Isabeau Birindelli, Françoise Demengel, Fabiana Leoni. Principal eigenvalues and eigenfunctions for fully nonlinear equations in punctured balls. Journal de Mathématiques Pures et Appliquées, 2024, 186, pp.74-102. 10.1016/j.matpur.2024.04.004. hal-04596489

HAL Id: hal-04596489

https://hal.science/hal-04596489

Submitted on 31 May 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

Public Domain
# Principal eigenvalues and eigenfunctions for fully nonlinear equations in punctured balls

Isabeau Birindelli, Françoise Demengel, Fabiana Leoni

► To cite this version:

Isabeau Birindelli, Françoise Demengel, Fabiana Leoni. Principal eigenvalues and eigenfunctions for fully nonlinear equations in punctured balls. Journal de Mathématiques Pures et Appliquées, 2024, 186, pp.74-102. 10.1016/j.matpur.2024.04.004. hal-04596489

HAL Id: hal-04596489

https://hal.science/hal-04596489v1

Submitted on 31 May 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

Public Domain
# Principal eigenvalues and eigenfunctions for fully nonlinear equations in punctured balls

Isabeau Birindelli¹

Sapienza Università di Roma, P.le Aldo Moro 2, 00185 Roma, Italy

Françoise Demengel

CY Paris University, 33, Boulevard du Port 95011 Cergy-Pontoise Cedex, France

Fabiana Leoni

Sapienza Università di Roma, P.le Aldo Moro 2, 00185 Roma, Italy

## Abstract

This paper is devoted to the proof of the existence of the principal eigenvalue and related eigenfunctions for fully nonlinear uniformly elliptic equations posed in a punctured ball, in presence of a singular potential. More precisely, we analyze existence, uniqueness and regularity of solutions $(\bar{\lambda}_{\gamma}, u_{\gamma})$ of the equation

$$F(D^2 u_{\gamma}) + \bar{\lambda}_{\gamma} \frac{u_{\gamma}}{r^{\gamma}} = 0 \text{ in } B(0, 1) \setminus \{0\}, \quad u_{\gamma} = 0 \text{ on } \partial B(0, 1)$$

where $u_{\gamma} > 0$ in $B(0, 1) \setminus \{0\}$ and $\gamma > 0$. We prove existence of radial solutions which are continuous on $\bar{B}(0, 1)$ in the case $\gamma < 2$, existence of unbounded solutions in the case $\gamma = 2$ and a non existence result for $\gamma > 2$. We also give, in the case of Pucci's operators, the explicit value of $\bar{\lambda}_2$, which generalizes the Hardy-Sobolev constant for the Laplacian.

## Résumé

Dans cet article nous nous intéressons à l'existence de valeur propre principale et de fonctions propres associées pour des opérateurs complètement non linéaires dans un domaine épointé, en présence d'un potentiel singulier.

Plus pécisément, nous analysons l'existence, l'unicité et la régularité de solutions de l'équation

$$F(D^2 u_{\gamma}) + \bar{\lambda}_{\gamma} \frac{u_{\gamma}}{r^{\gamma}} = 0 \text{ dans } B(0, 1) \setminus \{0\}, \quad u_{\gamma} = 0 \text{ sur } \partial B(0, 1)$$

où $u_{\gamma} > 0$ est définie sur $B(0, 1) \setminus \{0\}$ et $\gamma > 0$. On montre l'existence de solutions radiales qui sont continuées sur $\bar{B}(0, 1)$ dans le cas $\gamma < 2$, l'existence de solutions non bornées dans le cas $\gamma = 2$ et un résultat de non existence dans le cas $\gamma > 2$. Nous donnons aussi dans
Principal eigenvalues and eigenfunctions for fully nonlinear equations
in punctured balls

Isabeau Birindelli¹

Sapienza Università di Roma, P.le Aldo Moro 2, 00185 Roma, Italy

Françoise Demengel

CY Paris University, 33, Boulevard du Port 95011 Cergy-Pontoise Cedex, France

Fabiana Leoni

Sapienza Università di Roma, P.le Aldo Moro 2, 00185 Roma, Italy

Abstract

This paper is devoted to the proof of the existence of the principal eigenvalue and related eigenfunctions for fully nonlinear uniformly elliptic equations posed in a punctured ball, in presence of a singular potential. More precisely, we analyze existence, uniqueness and regularity of solutions $(\bar{\lambda}_{\gamma}, u_{\gamma})$ of the equation

$$F(D^2u_{\gamma}) + \bar{\lambda}_{\gamma} \frac{u_{\gamma}}{r^{\gamma}} = 0 \text{ in } B(0,1) \setminus \{0\}, u_{\gamma} = 0 \text{ on } \partial B(0,1)$$

where $u_{\gamma} > 0$ in $B(0,1) \setminus \{0\}$ and $\gamma > 0$. We prove existence of radial solutions which are continuous on $\bar{B}(0,1)$ in the case $\gamma < 2$, existence of unbounded solutions in the case $\gamma = 2$ and a non existence result for $\gamma > 2$. We also give, in the case of Pucci's operators, the explicit value of $\bar{\lambda}_2$, which generalizes the Hardy-Sobolev constant for the Laplacian.

Résumé

Dans cet article nous nous intéressons à l'existence de valeur propre principale et de fonctions propres associées pour des opérateurs complètement non linéaires dans un domaine épointé, en présence d'un potentiel singulier.

Plus pécisément, nous analysons l'existence, l'unicité et la régularité de solutions de l'équation

$$F(D^2u_{\gamma}) + \bar{\lambda}_{\gamma} \frac{u_{\gamma}}{r^{\gamma}} = 0 \text{ dans } B(0,1) \setminus \{0\}, u_{\gamma} = 0 \text{ sur } \partial B(0,1)$$

où $u_{\gamma} > 0$ est définie sur $B(0,1) \setminus \{0\}$ et $\gamma > 0$. On montre l'existence de solutions radiales qui sont continuées sur $\bar{B}(0,1)$ dans le cas $\gamma < 2$, l'existence de solutions non bornées dans le cas $\gamma = 2$ et un résultat de non existence dans le cas $\gamma > 2$. Nous donnons aussi dans
le cas des opérateurs de Pucci la valeur explicite de $\bar{\lambda}_2$, ce qui généralise la constante de Hardy-Sobolev dans le cas du Laplacien.

Keywords: Fully nonlinear elliptic equations, singular potential, principal eigenvalues, regularity of eigenfunctions

2020 MSC: 35J60, 35J75, 35P15, 35P30.

# 1. Introduction

In this paper we will study radial eigenvalues and related positive radial eigenfunctions for the Dirichlet problem associated with the equation

$$F(D^2 u) + \mu r^{-\gamma} u = 0 \text{ in } B(0, 1) \setminus \{0\}$$

when $\gamma > 0$ and $F$ is a second order fully nonlinear uniformly elliptic operator. By radial eigenvalue and radial eigenfunction we mean a real value $\lambda_\gamma$ and a radial nontrivial function $u_\gamma$ respectively, satisfying the equation

$$\left\{ \begin{array}{ll} F(D^2 u_\gamma) + \lambda_\gamma r^{-\gamma} u_\gamma = 0 & \text{in } B(0, 1) \setminus \{0\} \\ u_\gamma = 0 & \text{on } \partial B(0, 1). \end{array} \right. \qquad (1.1)$$

In principle, eigenfunctions are required to satisfy the above eigenvalue problem in the viscosity sense, but, due to the radial symmetry, this is equivalent to consider classical solutions. We will focus on constant sign eigenfunctions, in particular positive eigenfunctions, thus referring to the so called principal eigenvalues. If necessary, in order to emphasize the dependence of the eigenvalue on the operator $F$, the potential $f(r)$ appearing in the zero order term and the domain $\Omega$ in which the equation is considered, we will use the notation $\lambda = \lambda(F, f(r), \Omega)$. Interestingly, we will see that for problem (1.1), as in the case when $F$ is the Laplace operator, $\gamma = 2$ is a critical value, in the sense that for $\gamma < 2$ there exists smooth eigenfunctions, for $\gamma > 2$ there are no eigenfunctions and, for $\gamma = 2$, the eigenfunctions are unbounded.

Let us recall some known results when $F$ is the Laplacian. In the case $\gamma = 2$, the equation is naturally linked to Hardy's inequality. Indeed, if $N > 2$ and $u \in H_0^1(B(0, 1))$ (respectively, $u \in H^1(\mathbb{R}^N)$), then $\frac{u(x)}{|x|}$ belongs to $L^2(B(0, 1))$ (respectively $\frac{u(x)}{|x|} \in L^2(\mathbb{R}^N)$), and there exists a positive constant $c$ such that

$$\int \left( \frac{|u(x)|}{|x|} \right)^2 \le c \int |\nabla u|^2.$$

<footnote>Email addresses: isabeau@mat.uniroma1.it (Isabeau Birindelli), francoise.demengel@cyu.fr (Françoise Demengel), leoni@mat.uniroma1.it (Fabiana Leoni)</footnote>

<footnote n="1">¹Corresponding author</footnote>
le cas des opérateurs de Pucci la valeur explicite de $\bar{\lambda}_2$, ce qui généralise la constante de Hardy-Sobolev dans le cas du Laplacien.

*Keywords:* Fully nonlinear elliptic equations, singular potential, principal eigenvalues, regularity of eigenfunctions

2020 MSC: 35J60, 35J75, 35P15, 35P30.

# 1. Introduction

In this paper we will study radial eigenvalues and related positive radial eigenfunctions for the Dirichlet problem associated with the equation

$$F(D^2u) + \mu r^{-\gamma}u = 0 \quad \text{in } B(0, 1) \setminus \{0\}$$

when $\gamma > 0$ and $F$ is a second order fully nonlinear uniformly elliptic operator. By radial eigenvalue and radial eigenfunction we mean a real value $\lambda_\gamma$ and a radial nontrivial function $u_\gamma$ respectively, satisfying the equation

$$\left\{ \begin{array}{ll} F(D^2u_\gamma) + \lambda_\gamma r^{-\gamma}u_\gamma = 0 & \text{in } B(0, 1) \setminus \{0\} \\ u_\gamma = 0 & \text{on } \partial B(0, 1). \end{array} \right. \qquad (1.1)$$

In principle, eigenfunctions are required to satisfy the above eigenvalue problem in the viscosity sense, but, due to the radial symmetry, this is equivalent to consider classical solutions. We will focus on constant sign eigenfunctions, in particular positive eigenfunctions, thus referring to the so called principal eigenvalues. If necessary, in order to emphasize the dependence of the eigenvalue on the operator $F$, the potential $f(r)$ appearing in the zero order term and the domain $\Omega$ in which the equation is considered, we will use the notation $\lambda = \lambda(F, f(r), \Omega)$. Interestingly, we will see that for problem (1.1), as in the case when $F$ is the Laplace operator, $\gamma = 2$ is a critical value, in the sense that for $\gamma < 2$ there exists smooth eigenfunctions, for $\gamma > 2$ there are no eigenfunctions and, for $\gamma = 2$, the eigenfunctions are unbounded.

Let us recall some known results when $F$ is the Laplacian. In the case $\gamma = 2$, the equation is naturally linked to Hardy's inequality. Indeed, if $N > 2$ and $u \in H_0^1(B(0, 1))$ (respectively, $u \in H^1(\mathbb{R}^N)$), then $\frac{u(x)}{|x|}$ belongs to $L^2(B(0, 1))$ (respectively $\frac{u(x)}{|x|} \in L^2(\mathbb{R}^N)$), and there exists a positive constant $c$ such that

$$\int \left( \frac{|u(x)|}{|x|} \right)^2 \le c \int |\nabla u|^2.$$

<footnote>Email addresses: isabeau@mat.uniroma1.it (Isabeau Birindelli), francoise.demengel@cyu.fr (Françoise Demengel), leoni@mat.uniroma1.it (Fabiana Leoni)</footnote>

<footnote n="1">¹Corresponding author</footnote>
Furthermore, the best constant $c = \frac{4}{(N-2)^2}$ is not achieved, in the sense that

$$ \inf_{u \in H_0^1(B(0,1)), \int_{B(0,1)} \left( \frac{|u(x)|}{|x|} \right)^2 dxdy = 1} \int_{B(0,1)} |\nabla u|^2 = \frac{(N-2)^2}{4} \quad (1.2) $$

but there is no $u \in H_0^1$ which realizes the infimum. By obvious arguments, $B(0,1)$ can be replaced by any bounded regular open set of $\mathbb{R}^N$ containing 0, and the optimal constant does not depend on the size of $\Omega$. Note that the right hand side of (1.2) coincides with the variational characterization of the first (or principal) eigenvalue for the equation

$$ -\Delta u = \lambda \frac{u}{|x|^2}. $$

For further knowledge on the Hardy–Sobolev inequality and for the case of the p-Laplacian, we refer to [1, 2, 3].

On the other hand, if the exponent $\gamma$ of the potential is strictly less than 2, since $H_0^1(B(0,1))$ is compactly embedded into the weighted space $L^2(B(0,1), \frac{1}{r^\gamma})$, then existence of minima in $H_0^1(B(0,1))$ can be obtained by standard arguments of the direct method in calculus of variations. In that case, denoting

$$ \bar{\lambda}_{\gamma} = \inf_{u \in H_0^1(B(0,1)), \int_{B(0,1)} \frac{|u(x)|^2}{|x|^{\gamma}} dx = 1} \int_{B(0,1)} |\nabla u|^2 \quad (1.3) $$

one sees that $\bar{\lambda}_{\gamma}$ is also the first eigenvalue for the equation

$$ \Delta u + \bar{\lambda}_{\gamma} \frac{u}{r^{\gamma}} = 0, $$

meaning that $\bar{\lambda}_{\gamma}$ is such that there exists $u > 0$ in $H_0^1(B(0,1))$ satisfying the equation. Note that, by its definition, $\bar{\lambda}_{\gamma}$ depends on the domain, since

$$ \bar{\lambda}_{\gamma}(B(0, t)) = \frac{1}{t^{2-\gamma}} \bar{\lambda}_{\gamma}(B(0, 1)). $$

If $\gamma > 2$ there is no embedding from $H_0^1(B(0,1))$ into $L^2(B(0,1), \frac{1}{r^\gamma})$. Indeed, as an example, the function

$$ u(r) = r^{-\frac{N-2}{2}} + \epsilon(-\log r) $$

with $0 < \epsilon < \frac{\gamma-2}{2}$, belongs to $H_0^1(B(0,1))$ and satisfies

$$ \int_{B(0,1)} \frac{u(|x|)^2}{|x|^{\gamma}} dx = +\infty $$
Furthermore, the best constant $c = \frac{4}{(N-2)^2}$ is not achieved, in the sense that

$$ \inf_{u \in H_0^1(B(0,1)), \int_{B(0,1)} \left( \frac{|u(x)|}{|x|} \right)^2 dxdy = 1} \int_{B(0,1)} |\nabla u|^2 = \frac{(N-2)^2}{4} \quad (1.2) $$

but there is no $u \in H_0^1$ which realizes the infimum. By obvious arguments, $B(0,1)$ can be replaced by any bounded regular open set of $\mathbb{R}^N$ containing 0, and the optimal constant does not depend on the size of $\Omega$. Note that the right hand side of (1.2) coincides with the variational characterization of the first (or principal) eigenvalue for the equation

$$ -\Delta u = \lambda \frac{u}{|x|^2}. $$

For further knowledge on the Hardy–Sobolev inequality and for the case of the p-Laplacian, we refer to [1, 2, 3].

On the other hand, if the exponent $\gamma$ of the potential is strictly less than 2, since $H_0^1(B(0,1))$ is compactly embedded into the weighted space $L^2(B(0,1), \frac{1}{r^\gamma})$, then existence of minima in $H_0^1(B(0,1))$ can be obtained by standard arguments of the direct method in calculus of variations. In that case, denoting

$$ \bar{\lambda}_{\gamma} = \inf_{u \in H_0^1(B(0,1)), \int_{B(0,1)} \frac{|u(x)|^2}{|x|^{\gamma}} dx = 1} \int_{B(0,1)} |\nabla u|^2 \quad (1.3) $$

one sees that $\bar{\lambda}_{\gamma}$ is also the first eigenvalue for the equation

$$ \Delta u + \bar{\lambda}_{\gamma} \frac{u}{r^{\gamma}} = 0, $$

meaning that $\bar{\lambda}_{\gamma}$ is such that there exists $u > 0$ in $H_0^1(B(0,1))$ satisfying the equation. Note that, by its definition, $\bar{\lambda}_{\gamma}$ depends on the domain, since

$$ \bar{\lambda}_{\gamma}(B(0, t)) = \frac{1}{t^{2-\gamma}} \bar{\lambda}_{\gamma}(B(0, 1)). $$

If $\gamma > 2$ there is no embedding from $H_0^1(B(0,1))$ into $L^2(B(0,1), \frac{1}{r^\gamma})$. Indeed, as an example, the function

$$ u(r) = r^{-\frac{N-2}{2}} + \epsilon(-\log r) $$

with $0 < \epsilon < \frac{\gamma-2}{2}$, belongs to $H_0^1(B(0,1))$ and satisfies

$$ \int_{B(0,1)} \frac{u(|x|)^2}{|x|^{\gamma}} dx = +\infty $$
For results in the variational linear case we refer to the works of many authors, but in particular we wish to mention the works of Cirstea and collaborators [4, 5, 6, 7].

Let us now focus on the case of concern of this paper i.e. when $F$ is a fully nonlinear uniformly elliptic operator, that is $F$ is a continuous function defined on the set $S_N$ of symmetric $N \times N$ matrices, and it satisfies, for positive constants $\Lambda \ge \lambda > 0$,

$$ \lambda \operatorname{tr}(M') \le F(M + M') - F(M) \le \Lambda \operatorname{tr}(M'), \quad (1.4) $$

for all $M, M' \in S_N$, with $M'$ positive semidefinite.

We suppose also that $F$ is rotationally invariant, that is

$$ F(O^t MO) = F(M) \quad (1.5) $$

for every orthogonal matrix $O$ and for all $M \in S_N$, and that $F$ is positively homogeneous of degree 1, i.e.

$$ F(tM) = tF(M) \quad (1.6) $$

for any $M \in S_N$ and for all $t > 0$. In this case we will see that, as in the regular case i.e. $\gamma = 0$, the first eigenvalue for problem (1.1) can be defined on the model of [8], i.e. by the optimization formula

$$ \begin{aligned} \bar{\lambda}_{\gamma} &= \bar{\lambda}_{\gamma}(F, r^{-\gamma}, B(0,1) \setminus \{0\}) \\ &= \sup\{\mu : \exists u \in C(B(0,1) \setminus \{0\}), u > 0 \text{ in } B(0,1) \setminus \{0\}, F(D^2u) + \mu \frac{u}{r^{\gamma}} \le 0\}, \end{aligned} \quad (1.7) $$

where the differential inequality is understood in the viscosity sense.

A first easy observation is that, by considering constant super-solutions, one always has $\bar{\lambda}_{\gamma} \ge 0$. One of the goals of the present paper is to show, in particular, that $\bar{\lambda}_{\gamma} > 0$ for $\gamma \le 2$. In case of smooth coefficients and regular domains, the principal eigenvalues and related eigenfunctions for fully nonlinear operators $F$ have been largely investigated. We refer to e.g. [9, 10, 11, 12, 13, 14].

In the class of fully nonlinear operators satisfying (1.4), (1.5) and (1.6), the Pucci's extremal operators play a crucial role and we will treat them in depth. Let us recall that, by decomposing each matrix $M \in S_N$ as $M = M^+ - M^-$, where $M^+$ and $M^-$ are positive semidefinite matrices satisfying $M^+M^- = O$, then Pucci's sup operator can be defined as

$$ \mathcal{M}_{\lambda, \Lambda}^{+}(M) = \Lambda \operatorname{tr}(M^{+}) - \lambda \operatorname{tr}(M^{-}), $$

as well as Pucci's inf operator is given by

$$ \mathcal{M}_{\lambda, \Lambda}^{-}(M) = \lambda \operatorname{tr}(M^{+}) - \Lambda \operatorname{tr}(M^{-}) = -\mathcal{M}_{\lambda, \Lambda}^{+}(-M). $$

As it is well known, see [15], under assumption (1.4) each operator $F$ satisfies

$$ \mathcal{M}_{\lambda, \Lambda}^{-}(M_1 - M_2) \le F(M_1) - F(M_2) \le \mathcal{M}_{\lambda, \Lambda}^{+}(M_1 - M_2) \quad (1.8) $$
For results in the variational linear case we refer to the works of many authors, but in particular we wish to mention the works of Cirstea and collaborators [4, 5, 6, 7].

Let us now focus on the case of concern of this paper i.e. when $F$ is a fully nonlinear uniformly elliptic operator, that is $F$ is a continuous function defined on the set $S_N$ of symmetric $N \times N$ matrices, and it satisfies, for positive constants $\Lambda \ge \lambda > 0$,

$$ \lambda \operatorname{tr}(M') \le F(M + M') - F(M) \le \Lambda \operatorname{tr}(M'), \quad (1.4) $$

for all $M, M' \in S_N$, with $M'$ positive semidefinite.

We suppose also that $F$ is rotationally invariant, that is

$$ F(O^t MO) = F(M) \quad (1.5) $$

for every orthogonal matrix $O$ and for all $M \in S_N$, and that $F$ is positively homogeneous of degree 1, i.e.

$$ F(tM) = tF(M) \quad (1.6) $$

for any $M \in S_N$ and for all $t > 0$. In this case we will see that, as in the regular case i.e. $\gamma = 0$, the first eigenvalue for problem (1.1) can be defined on the model of [8], i.e. by the optimization formula

$$ \begin{aligned} \bar{\lambda}_{\gamma} &= \bar{\lambda}_{\gamma}(F, r^{-\gamma}, B(0,1) \setminus \{0\}) \\ &= \sup\{\mu : \exists u \in C(B(0,1) \setminus \{0\}), u > 0 \text{ in } B(0,1) \setminus \{0\}, F(D^2u) + \mu \frac{u}{r^{\gamma}} \le 0\}, \end{aligned} \quad (1.7) $$

where the differential inequality is understood in the viscosity sense.

A first easy observation is that, by considering constant super-solutions, one always has $\bar{\lambda}_{\gamma} \ge 0$. One of the goals of the present paper is to show, in particular, that $\bar{\lambda}_{\gamma} > 0$ for $\gamma \le 2$. In case of smooth coefficients and regular domains, the principal eigenvalues and related eigenfunctions for fully nonlinear operators $F$ have been largely investigated. We refer to e.g. [9, 10, 11, 12, 13, 14].

In the class of fully nonlinear operators satisfying (1.4), (1.5) and (1.6), the Pucci's extremal operators play a crucial role and we will treat them in depth. Let us recall that, by decomposing each matrix $M \in S_N$ as $M = M^+ - M^-$, where $M^+$ and $M^-$ are positive semidefinite matrices satisfying $M^+M^- = O$, then Pucci's sup operator can be defined as

$$ \mathcal{M}_{\lambda, \Lambda}^{+}(M) = \Lambda \operatorname{tr}(M^{+}) - \lambda \operatorname{tr}(M^{-}), $$

as well as Pucci's inf operator is given by

$$ \mathcal{M}_{\lambda, \Lambda}^{-}(M) = \lambda \operatorname{tr}(M^{+}) - \Lambda \operatorname{tr}(M^{-}) = -\mathcal{M}_{\lambda, \Lambda}^{+}(-M). $$

As it is well known, see [15], under assumption (1.4) each operator $F$ satisfies

$$ \mathcal{M}_{\lambda, \Lambda}^{-}(M_1 - M_2) \le F(M_1) - F(M_2) \le \mathcal{M}_{\lambda, \Lambda}^{+}(M_1 - M_2) \quad (1.8) $$
for all $M_1, M_2 \in \mathcal{S}_N$, and, if (1.6) is also satisfied, then

$$ \mathcal{M}_{\lambda, \Lambda}^{-}(M) \leq F(M) \leq \mathcal{M}_{\lambda, \Lambda}^{+}(M), \quad \forall M \in \mathcal{S}_{N}. \qquad (1.9) $$

This shows in particular that Pucci's operators act as explicit extremal operators in the whole class of uniformly elliptic operators having the same ellipticity constants. In the sequel, we will omit in the notation the dependence on the ellipticity constants, which are fixed once for all.

We further recall that for a $C^2$ radial function $u(x) = u(|x|)$, one has

$$ D^2 u(x) = u''(r) \frac{x \otimes x}{r^2} + \frac{u'(r)}{r} \left( I - \frac{x \otimes x}{r^2} \right), $$

and, as a consequence,

$$ \mathcal{M}^{+}(D^{2}u) = \Lambda(N-1) \left(\frac{u'(r)}{r}\right)^{+} - \lambda(N-1) \left(\frac{u'(r)}{r}\right)^{-} + \Lambda(u''(r))^{+} - \lambda(u''(r))^{-}, $$

$$ \mathcal{M}^{-}(D^{2}u) = \lambda(N-1) \left(\frac{u'(r)}{r}\right)^{+} - \Lambda(N-1) \left(\frac{u'(r)}{r}\right)^{-} + \lambda(u''(r))^{+} - \Lambda(u''(r))^{-}. $$

Thus, the ODEs satisfied by radial solutions of Pucci's extremal equations have coefficients depending on the dimension like parameters, associated with $\mathcal{M}^+$ and $\mathcal{M}^-$ respectively, defined as

$$ \tilde{N}_{+} = \frac{\lambda}{\Lambda}(N-1) + 1, \quad \tilde{N}_{-} = \frac{\Lambda}{\lambda}(N-1) + 1. $$

Note that one has always

$$ \tilde{N}_{-} \geq N \geq \tilde{N}_{+}, $$

with equalities holding true if and only if $\Lambda = \lambda$. We will assume always that $\tilde{N}_{+} > 2$.

We can now state the main results of the paper. We will always assume that $F$ satisfies assumptions (1.4), (1.5) and (1.6). Let us start with the case $\gamma < 2$.

**Theorem 1.1.** Suppose that $\gamma < 2$. Then:

(i) $\bar{\lambda}_{\gamma}$ defined in (1.7) is positive and there exists a function $u$, continuous in $\overline{B(0, 1)}$, radial, strictly positive in $B(0, 1)$, such that

$$ \left\{ \begin{array}{ll} F(D^2u) + \bar{\lambda}_{\gamma} \frac{u}{r^{\gamma}} = 0 & \text{in } B(0,1) \setminus \{0\} \\ u = 0 & \text{on } \partial B(0,1) \end{array} \right. . $$

Furthermore $u$ is $C^2(B(0,1) \setminus \{0\})$ and it can be extended on $B(0,1)$ as a Lipschitz continuous function if $\gamma \le 1$, as a function of class $C^1(B(0,1))$ when $\gamma < 1$, and as an Hölder continuous function with exponent $2 - \gamma$ if $\gamma > 1$.
for all $M_1, M_2 \in \mathcal{S}_N$, and, if (1.6) is also satisfied, then

$$ \mathcal{M}_{\lambda, \Lambda}^{-}(M) \leq F(M) \leq \mathcal{M}_{\lambda, \Lambda}^{+}(M), \quad \forall M \in \mathcal{S}_{N}. \qquad (1.9) $$

This shows in particular that Pucci's operators act as explicit extremal operators in the whole class of uniformly elliptic operators having the same ellipticity constants. In the sequel, we will omit in the notation the dependence on the ellipticity constants, which are fixed once for all.

We further recall that for a $C^2$ radial function $u(x) = u(|x|)$, one has

$$ D^2 u(x) = u''(r) \frac{x \otimes x}{r^2} + \frac{u'(r)}{r} \left( I - \frac{x \otimes x}{r^2} \right), $$

and, as a consequence,

$$ \mathcal{M}^{+}(D^{2}u) = \Lambda(N-1) \left(\frac{u'(r)}{r}\right)^{+} - \lambda(N-1) \left(\frac{u'(r)}{r}\right)^{-} + \Lambda(u''(r))^{+} - \lambda(u''(r))^{-}, $$

$$ \mathcal{M}^{-}(D^{2}u) = \lambda(N-1) \left(\frac{u'(r)}{r}\right)^{+} - \Lambda(N-1) \left(\frac{u'(r)}{r}\right)^{-} + \lambda(u''(r))^{+} - \Lambda(u''(r))^{-}. $$

Thus, the ODEs satisfied by radial solutions of Pucci's extremal equations have coefficients depending on the dimension like parameters, associated with $\mathcal{M}^+$ and $\mathcal{M}^-$ respectively, defined as

$$ \tilde{N}_{+} = \frac{\lambda}{\Lambda}(N-1) + 1, \quad \tilde{N}_{-} = \frac{\Lambda}{\lambda}(N-1) + 1. $$

Note that one has always

$$ \tilde{N}_{-} \geq N \geq \tilde{N}_{+}, $$

with equalities holding true if and only if $\Lambda = \lambda$. We will assume always that $\tilde{N}_{+} > 2$.

We can now state the main results of the paper. We will always assume that $F$ satisfies assumptions (1.4), (1.5) and (1.6). Let us start with the case $\gamma < 2$.

**Theorem 1.1.** Suppose that $\gamma < 2$. Then:

(i) $\bar{\lambda}_{\gamma}$ defined in (1.7) is positive and there exists a function $u$, continuous in $\overline{B(0, 1)}$, radial, strictly positive in $B(0, 1)$, such that

$$ \left\{ \begin{array}{ll} F(D^2u) + \bar{\lambda}_{\gamma} \frac{u}{r^{\gamma}} = 0 & \text{in } B(0,1) \setminus \{0\} \\ u = 0 & \text{on } \partial B(0,1) \end{array} \right. . $$

Furthermore $u$ is $C^2(B(0,1) \setminus \{0\})$ and it can be extended on $B(0,1)$ as a Lipschitz continuous function if $\gamma \le 1$, as a function of class $C^1(B(0,1))$ when $\gamma < 1$, and as an Hölder continuous function with exponent $2 - \gamma$ if $\gamma > 1$.
(ii) $\bar{\lambda}_{\gamma}$ is stable under various regular approximations :

$$ \bar{\lambda}_{\gamma} = \lim_{\epsilon \to 0} \bar{\lambda}(F, \frac{1}{(r^2 + \epsilon^2)^{\frac{\gamma}{2}}}, B(0, 1)), $$

$$ \bar{\lambda}_{\gamma} = \lim_{\delta \to 0} \bar{\lambda}(F, \frac{1}{r^{\gamma}}, B(0, 1) \setminus \overline{B(0, \delta)}). $$

Statement (i) of the above theorem shows in particular that $\bar{\lambda}_{\gamma}$ is actually achieved on smooth radial eigenfunctions. Thus, if we define

$$ \bar{\lambda}'_{\gamma} := \sup\{\mu : \exists u \in C^2(B(0,1) \setminus \{0\}), u > 0 \text{ in } B(0,1) \setminus \{0\}, u \text{ radial}, \\ F(D^2u) + \mu \frac{u}{r^{\gamma}} \le 0\}, \quad (1.10) $$

it then follows that

$$ \bar{\lambda}_{\gamma} = \bar{\lambda}'_{\gamma}. $$

Actually, we will work initially with the smooth eigenvalue $\bar{\lambda}'_{\gamma}$, and we will finally prove that it coincides with $\bar{\lambda}_{\gamma}$. Note that, due to the lack of regularity of the coefficient function $\frac{1}{r^{\gamma}}$ we cannot employ directly the results of [16] which ensure that solutions of

$$ F(D^2u) + f(r)u = 0 $$

in a radial domain are radial when $f$ is non increasing. Nonetheless, we will prove that the two eigenvalues coincide for any $\gamma \le 2$. However, due to the singularity at zero, we cannot prove that any eigenfunction is radial. For symmetry results for solutions of equations having isolated singularities, we refer to [17], where the case of $C^1$ operators is considered.

Theorem 1.1 will be proved after several steps and intermediate results. In particular, we will prove a comparison theorem for smooth, bounded, radial sub- and super-solutions in the punctured ball, without assuming any order condition at the origin. Furthermore, we will show that for $\mu < \bar{\lambda}_{\gamma}$ the problem

$$ \left\{ \begin{array}{ll} F(D^2u) + \mu u r^{-\gamma} = f(r)r^{-\gamma} & \text{in } B(0,1) \setminus \{0\} \\ u = 0 & \text{on } \partial B(0,1) \end{array} \right. $$

admits a unique radial solution $u \in C^2(B(0,1)\setminus\{0\}) \cap C(\overline{B(0,1)})$ for any radial and continuous datum $f \in \overline{B(0,1)}$ satisfying $f \le 0$.

Next, for the case $\gamma = 2$ we have the following result, which gives explicit expressions for the eigenvalues and the eigenfunctions in case of Pucci's operators.

**Theorem 1.2.** Assume that $\gamma = 2$. Then:
(ii) $\bar{\lambda}_{\gamma}$ is stable under various regular approximations :

$$ \bar{\lambda}_{\gamma} = \lim_{\epsilon \to 0} \bar{\lambda}(F, \frac{1}{(r^2 + \epsilon^2)^{\frac{\gamma}{2}}}, B(0, 1)), $$

$$ \bar{\lambda}_{\gamma} = \lim_{\delta \to 0} \bar{\lambda}(F, \frac{1}{r^{\gamma}}, B(0, 1) \setminus \overline{B(0, \delta)}). $$

Statement (i) of the above theorem shows in particular that $\bar{\lambda}_{\gamma}$ is actually achieved on smooth radial eigenfunctions. Thus, if we define

$$ \bar{\lambda}'_{\gamma} := \sup\{\mu : \exists u \in C^2(B(0,1) \setminus \{0\}), u > 0 \text{ in } B(0,1) \setminus \{0\}, u \text{ radial}, \\ F(D^2u) + \mu \frac{u}{r^{\gamma}} \le 0\}, \quad (1.10) $$

it then follows that

$$ \bar{\lambda}_{\gamma} = \bar{\lambda}'_{\gamma}. $$

Actually, we will work initially with the smooth eigenvalue $\bar{\lambda}'_{\gamma}$, and we will finally prove that it coincides with $\bar{\lambda}_{\gamma}$. Note that, due to the lack of regularity of the coefficient function $\frac{1}{r^{\gamma}}$ we cannot employ directly the results of [16] which ensure that solutions of

$$ F(D^2u) + f(r)u = 0 $$

in a radial domain are radial when $f$ is non increasing. Nonetheless, we will prove that the two eigenvalues coincide for any $\gamma \le 2$. However, due to the singularity at zero, we cannot prove that any eigenfunction is radial. For symmetry results for solutions of equations having isolated singularities, we refer to [17], where the case of $C^1$ operators is considered.

Theorem 1.1 will be proved after several steps and intermediate results. In particular, we will prove a comparison theorem for smooth, bounded, radial sub- and super-solutions in the punctured ball, without assuming any order condition at the origin. Furthermore, we will show that for $\mu < \bar{\lambda}_{\gamma}$ the problem

$$ \left\{ \begin{array}{ll} F(D^2u) + \mu u r^{-\gamma} = f(r)r^{-\gamma} & \text{in } B(0,1) \setminus \{0\} \\ u = 0 & \text{on } \partial B(0,1) \end{array} \right. $$

admits a unique radial solution $u \in C^2(B(0,1)\setminus\{0\}) \cap C(\overline{B(0,1)})$ for any radial and continuous datum $f \in \overline{B(0,1)}$ satisfying $f \le 0$.

Next, for the case $\gamma = 2$ we have the following result, which gives explicit expressions for the eigenvalues and the eigenfunctions in case of Pucci's operators.

**Theorem 1.2.** Assume that $\gamma = 2$. Then:
(i) For the operator $M^+$ one has

$$ \bar{\lambda}_2(M^+) = \Lambda \frac{(\tilde{N}_+ - 2)^2}{4} $$

and the function $u(x) = r^{-\frac{\tilde{N}_+-2}{2}}(-\ln r)$ is an explicit solution of

$$ \begin{cases} M^+(D^2u) + \bar{\lambda}_2 \frac{u}{r^2} = 0 & \text{in } B(0,1) \setminus \{0\} \\ u = 0 & \text{on } \partial B(0,1) \end{cases} $$

Analogously, for the operator $M^-$ one has

$$ \bar{\lambda}_2(M^-) = \lambda \frac{(\tilde{N}_- - 2)^2}{4} $$

and the function $u(x) = r^{-\frac{\tilde{N}_--2}{2}}(-\ln r)$ is an explicit solution of

$$ \begin{cases} M^-(D^2u) + \bar{\lambda}_2 \frac{u}{r^2} = 0 & \text{in } B(0,1) \setminus \{0\} \\ u = 0 & \text{on } \partial B(0,1) \end{cases} $$

(ii) The eigenvalues $\bar{\lambda}_2(M^\pm)$ are stable under various regularization

$$ \bar{\lambda}_2(M^\pm) = \lim_{\gamma \to 2} \bar{\lambda}_\gamma(M^\pm) $$

$$ \bar{\lambda}_2(M^\pm) = \lim_{\delta \to 0} \bar{\lambda}(M^\pm, B(0,1) \setminus \overline{B(0, \delta)}) $$

$$ \bar{\lambda}_2(M^\pm) = \lim_{\epsilon \to 0} \bar{\lambda}(M^\pm, \frac{1}{r^2 + \epsilon^2}) $$

(iii) For any operator $F$ satisfying (1.4) and (1.6) one has

$$ \Lambda \frac{(\tilde{N}_+ - 2)^2}{4} \le \bar{\lambda}_2(F) \le \Lambda \frac{(\tilde{N}_- - 2)^2}{4} $$

We observe that we cannot prove the existence of eigenfunctions for a general operator $F$, but we can merely provide the estimate on the eigenvalue given by statement (iii) above. Theorem 1.2 will be obtained by using a variational approach adapted to the fully nonlinear radial framework. Indeed, we will define variational eigenvalues associated with the operators $M^\pm$ in an analogous way as in (1.3), taking advantage of the radial symmetry of solutions. Then, the full statements of Theorem 1.2 will follow as consequences of the properties established for $\bar{\lambda}_\gamma$ in the case $\gamma < 2$ and the stability of the variational formulation as $\gamma \to 2$.

Finally, for the case $\gamma > 2$, the singularity of the coefficient is too strong and it prevents the existence of positive smooth super-solutions, as stated by the following non existence result.
(i) For the operator $M^+$ one has

$$ \bar{\lambda}_2(M^+) = \Lambda \frac{(\tilde{N}_+ - 2)^2}{4} $$

and the function $u(x) = r^{-\frac{\tilde{N}_+ - 2}{2}}(-\ln r)$ is an explicit solution of

$$ \begin{cases} M^+(D^2u) + \bar{\lambda}_2 \frac{u}{r^2} = 0 & \text{in } B(0,1) \setminus \{0\} \\ u = 0 & \text{on } \partial B(0,1) \end{cases} $$

Analogously, for the operator $M^-$ one has

$$ \bar{\lambda}_2(M^-) = \lambda \frac{(\tilde{N}_- - 2)^2}{4} $$

and the function $u(x) = r^{-\frac{\tilde{N}_- - 2}{2}}(-\ln r)$ is an explicit solution of

$$ \begin{cases} M^-(D^2u) + \bar{\lambda}_2 \frac{u}{r^2} = 0 & \text{in } B(0,1) \setminus \{0\} \\ u = 0 & \text{on } \partial B(0,1) \end{cases} $$

(ii) The eigenvalues $\bar{\lambda}_2(M^\pm)$ are stable under various regularization

$$ \bar{\lambda}_2(M^\pm) = \lim_{\gamma \to 2} \bar{\lambda}_\gamma(M^\pm) $$

$$ \bar{\lambda}_2(M^\pm) = \lim_{\delta \to 0} \bar{\lambda}(M^\pm, B(0,1) \setminus \overline{B(0,\delta)}) $$

$$ \bar{\lambda}_2(M^\pm) = \lim_{\epsilon \to 0} \bar{\lambda}(M^\pm, \frac{1}{r^2 + \epsilon^2}) $$

(iii) For any operator $F$ satisfying (1.4) and (1.6) one has

$$ \Lambda \frac{(\tilde{N}_+ - 2)^2}{4} \le \bar{\lambda}_2(F) \le \Lambda \frac{(\tilde{N}_- - 2)^2}{4} $$

We observe that we cannot prove the existence of eigenfunctions for a general operator $F$, but we can merely provide the estimate on the eigenvalue given by statement (iii) above. Theorem 1.2 will be obtained by using a variational approach adapted to the fully nonlinear radial framework. Indeed, we will define variational eigenvalues associated with the operators $M^\pm$ in an analogous way as in (1.3), taking advantage of the radial symmetry of solutions. Then, the full statements of Theorem 1.2 will follow as consequences of the properties established for $\bar{\lambda}_\gamma$ in the case $\gamma < 2$ and the stability of the variational formulation as $\gamma \to 2$. Finally, for the case $\gamma > 2$, the singularity of the coefficient is too strong and it prevents the existence of positive smooth super-solutions, as stated by the following non existence result.
**Theorem 1.3.** If $\gamma > 2$, then the eigenvalue $\bar{\lambda}'_\gamma$ defined by (1.10) satisfies $\bar{\lambda}'_\gamma = 0$.

Let us observe that, symmetrically, one could define the eigenvalue associated with negative eigenfunctions, by setting

$$ \bar{\lambda}_{\gamma}^{-} = \sup\{\mu : \exists u \in C(B(0,1) \setminus \{0\}), u < 0 \text{ in } B(0,1), F(D^2u) + \mu \frac{u}{r^{\gamma}} \ge 0 \} = 0. $$

In this case, the results above can be extended to $\bar{\lambda}_{\gamma}^{-}$ with obvious modifications.

Let us conclude this introduction by observing that, in case of semilinear or quasilinear equations, many existence, non existence and classification results have been obtained in presence of zero order terms having Hardy's potential perturbed with additional sub- or superlinear terms. In particular, we refer to [18, 19, 4? ] for results related to Laplace operator, and to [1] for the p-Laplace operator.

The case where, in all directions above, $\Delta$ or $\Delta_p$ is replaced by a non variational fully nonlinear operator will be the object of future works.

## 2. The case $\gamma < 2$: proof of Theorem 1.1

Theorem 1.1 will be proved as a consequence of several classical steps: a comparison principle, existence and regularity results and a maximum principle "below" the first eigenvalue.

### 2.1. Maximum principles, existence and regularity results

The first result of the present section is a crucial technical lemma.

**Lemma 2.1.** Let $f \in C(B(0,1) \setminus \{0\})$ be a radial, bounded and positive function and assume that $u \in C^2(B(0,1) \setminus \{0\})$ is a radial, bounded function satisfying

$$ \mathcal{M}^+(D^2u) \geq f r^{-\gamma} \quad \text{in } B(0,1) \setminus \{0\}. \qquad (2.1) $$

Then

(i) $u' \ge 0$ in a right neighborhood of 0;

(ii) $\lim_{r \to 0} u'(r)r^{\tilde{N}-1} = 0$ and in a right neighborhood of 0 one has

$$ u'(r) \geq \frac{\inf f}{\Lambda(\tilde{N} - \gamma)} r^{1-\gamma}; $$
**Theorem 1.3.** If $\gamma > 2$, then the eigenvalue $\bar{\lambda}'_\gamma$ defined by (1.10) satisfies $\bar{\lambda}'_\gamma = 0$.

Let us observe that, symmetrically, one could define the eigenvalue associated with negative eigenfunctions, by setting

$$ \bar{\lambda}_{\gamma}^{-} = \sup\{\mu : \exists u \in C(B(0,1) \setminus \{0\}), u < 0 \text{ in } B(0,1), F(D^2u) + \mu \frac{u}{r^{\gamma}} \ge 0\} = 0. $$

In this case, the results above can be extended to $\bar{\lambda}_\gamma^-$ with obvious modifications.

Let us conclude this introduction by observing that, in case of semilinear or quasilinear equations, many existence, non existence and classification results have been obtained in presence of zero order terms having Hardy's potential perturbed with additional sub- or superlinear terms. In particular, we refer to [18, 19, 4? ] for results related to Laplace operator, and to [1] for the p-Laplace operator.

The case where, in all directions above, $\Delta$ or $\Delta_p$ is replaced by a non variational fully nonlinear operator will be the object of future works.

## 2. The case $\gamma < 2$: proof of Theorem 1.1

Theorem 1.1 will be proved as a consequence of several classical steps: a comparison principle, existence and regularity results and a maximum principle "below" the first eigenvalue.

### 2.1. Maximum principles, existence and regularity results

The first result of the present section is a crucial technical lemma.

**Lemma 2.1.** Let $f \in C(B(0,1) \setminus \{0\})$ be a radial, bounded and positive function and assume that $u \in C^2(B(0,1) \setminus \{0\})$ is a radial, bounded function satisfying

$$ \mathcal{M}^+(D^2u) \geq f r^{-\gamma} \quad \text{in } B(0,1) \setminus \{0\}. \qquad (2.1) $$

Then

(i) $u' \ge 0$ in a right neighborhood of 0;

(ii) $\lim_{r \to 0} u'(r)r^{\tilde{N}-1} = 0$ and in a right neighborhood of 0 one has

$$ u'(r) \geq \frac{\inf f}{\Lambda(\tilde{N} - \gamma)} r^{1-\gamma}; $$
(iii) if, furthermore, $M^+(D^2u) = fr^{-\gamma}$, then, in a right neighborhood of zero, one has also

$$u'(r) \leq \frac{\sup f}{\lambda(N - \gamma)} r^{1-\gamma}.$$

In particular, there exists a constant $c > 0$ such that, for $r$ sufficiently small,

$$|u'(r)| \leq cr^{1-\gamma}$$

and then $u$ is locally Lipschitz continuous in $B(0,1)$ if $\gamma \leq 1$, it belongs to $C^1(B(0,1))$ if $\gamma < 1$, and it is locally Hölder continuous in $B(0,1)$ with exponent $2 - \gamma$ if $\gamma > 1$.

*Proof.* Let us prove, by contradiction, that $u'$ does not change sign in a right neighborhood of 0. If not, there exists a decreasing sequence $\{r_n\}$ converging to 0, such that $u'(r_n) = 0$ for all $n$, and $u' \leq 0$ in $]r_{2n+1}, r_{2n}[$, $u' \geq 0$ in $]r_{2n+2}, r_{2n+1}[$. Since $u'(r_{2n}) = u'(r_{2n+1}) = 0$, there exists some $s_{2n} \in ]r_{2n+1}, r_{2n}]$ such that $u''(s_{2n}) = 0$. This yields the contradiction

$$0 \geq \lambda(N-1) \frac{u'(s_{2n})}{s_{2n}} = M^+(D^2u(s_{2n})) \geq f(s_{2n}) s_{2n}^{-\gamma} > 0.$$

Next, arguing again by contradiction, if $u'(r) \leq 0$ for $r$ sufficiently small, then, by (2.1), one has $u''(r) > 0$ and, in a right neighborhood of 0,

$$M^+(D^2u) = \Lambda u''(r) + \lambda(N-1) \frac{u'(r)}{r} \geq f(r)r^{-\gamma}.$$

Hence,

$$(u'(r)r^{\tilde{N}_+-1})' \geq \frac{fr^{\tilde{N}_+-1-\gamma}}{\Lambda} > 0$$

and, in particular, $u'r^{\tilde{N}_+-1}$ is increasing in a right neighborhood of 0. Thus, the limit $\lim_{r \to 0} u'(r)r^{\tilde{N}_+-1}$ exists and it is lesser than or equal to 0. If it was lesser than zero, then we would have, for some constant $l > 0$,

$$u'(r) \leq -lr^{1-\tilde{N}_+}$$

in a right neighborhood of 0, yielding a contradiction to the boundedness of $u$. This shows that $\lim_{r \to 0} u'(r)r^{\tilde{N}_+-1} = 0$ and, by monotonicity, $u'(r) > 0$ in a right neighborhood of 0. The reached contradiction proves statement (i).

In order to prove (ii), let us observe that, for $r$ sufficiently small, by (2.1) we have either

$$u''(r) + (N-1) \frac{u'(r)}{r} \geq \frac{f(r)r^{-\gamma}}{\Lambda}$$
(iii) if, furthermore, $M^+(D^2u) = fr^{-\gamma}$, then, in a right neighborhood of zero, one has also

$$u'(r) \leq \frac{\sup f}{\lambda(N - \gamma)} r^{1-\gamma}.$$

In particular, there exists a constant $c > 0$ such that, for $r$ sufficiently small,

$$|u'(r)| \leq cr^{1-\gamma}$$

and then $u$ is locally Lipschitz continuous in $B(0,1)$ if $\gamma \leq 1$, it belongs to $C^1(B(0,1))$ if $\gamma < 1$, and it is locally Hölder continuous in $B(0,1)$ with exponent $2 - \gamma$ if $\gamma > 1$.

*Proof.* Let us prove, by contradiction, that $u'$ does not change sign in a right neighborhood of 0. If not, there exists a decreasing sequence $\{r_n\}$ converging to 0, such that $u'(r_n) = 0$ for all $n$, and $u' \leq 0$ in $]r_{2n+1}, r_{2n}[$, $u' \geq 0$ in $]r_{2n+2}, r_{2n+1}[$. Since $u'(r_{2n}) = u'(r_{2n+1}) = 0$, there exists some $s_{2n} \in ]r_{2n+1}, r_{2n}]$ such that $u''(s_{2n}) = 0$. This yields the contradiction

$$0 \geq \lambda(N - 1) \frac{u'(s_{2n})}{s_{2n}} = M^+(D^2u(s_{2n})) \geq f(s_{2n}) s_{2n}^{-\gamma} > 0.$$

Next, arguing again by contradiction, if $u'(r) \leq 0$ for $r$ sufficiently small, then, by (2.1), one has $u''(r) > 0$ and, in a right neighborhood of 0,

$$M^+(D^2u) = \Lambda u''(r) + \lambda(N - 1) \frac{u'(r)}{r} \geq f(r)r^{-\gamma}.$$

Hence,

$$(u'(r)r^{\tilde{N}_+-1})' \geq \frac{fr^{\tilde{N}_+-1-\gamma}}{\Lambda} > 0$$

and, in particular, $u'r^{\tilde{N}_+-1}$ is increasing in a right neighborhood of 0. Thus, the limit $\lim_{r \to 0} u'(r)r^{\tilde{N}_+-1}$ exists and it is lesser than or equal to 0. If it was lesser than zero, then we would have, for some constant $l > 0$,

$$u'(r) \leq -lr^{1-\tilde{N}_+}$$

in a right neighborhood of 0, yielding a contradiction to the boundedness of $u$. This shows that $\lim_{r \to 0} u'(r)r^{\tilde{N}_+-1} = 0$ and, by monotonicity, $u'(r) > 0$ in a right neighborhood of 0. The reached contradiction proves statement (i).

In order to prove (ii), let us observe that, for $r$ sufficiently small, by (2.1) we have either

$$u''(r) + (N - 1) \frac{u'(r)}{r} \geq \frac{f(r)r^{-\gamma}}{\Lambda}$$
if $u''(r) \ge 0$, or

$$u''(r) + (\tilde{N}_- - 1) \frac{u'(r)}{r} \ge \frac{f(r)r^{-\gamma}}{\lambda}$$

if $u''(r) \le 0$. Since $u'(r) \ge 0$ and $\tilde{N}_- \ge N$, in both cases one has

$$u''(r) + (\tilde{N}_- - 1) \frac{u'(r)}{r} \ge \frac{f(r)r^{-\gamma}}{\Lambda},$$

that is

$$(u'r^{\tilde{N}_- - 1})' \ge \frac{f(r)r^{\tilde{N}_- - 1 - \gamma}}{\Lambda} \ge \frac{\inf f}{\Lambda} r^{\tilde{N}_- - 1 - \gamma}.$$

Arguing as above, we deduce that $u'r^{\tilde{N}_- - 1}$ is increasing in a right neighborhood of 0, hence it has a nonnegative limit as $r \to 0$, and such a limit must be 0, since $u$ is bounded. Moreover, by integrating the above inequality, we obtain

$$u'(r) \ge \frac{\inf f}{\Lambda(\tilde{N}_- - \gamma)} r^{1-\gamma}.$$

Let us finally prove (iii). Assuming that $M^+(D^2u) = f(r)r^{-\gamma}$ and using statement (i), it follows that, for every $r > 0$ sufficiently small, one has either

$$u''(r) + (N - 1) \frac{u'(r)}{r} = \frac{f(r)r^{-\gamma}}{\Lambda}$$

or

$$u''(r) + (\tilde{N}_- - 1) \frac{u'(r)}{r} = \frac{f(r)r^{-\gamma}}{\lambda}.$$

In both cases, we deduce

$$u''(r) + (N - 1) \frac{u'(r)}{r} \le \frac{f(r)r^{-\gamma}}{\lambda},$$

which yields

$$(u'(r)r^{N-1})' \le \frac{f(r)}{\lambda} r^{N-1-\gamma} \le \frac{\sup f}{\lambda} r^{N-1-\gamma}.$$

Hence, $u'(r)r^{N-1} - \frac{\sup f}{\lambda(N-\gamma)}r^{N-\gamma}$ is non increasing in a right neighborhood of 0 and it has a limit as $r \to 0$. This implies that $u'(r)r^{N-1}$ has a limit as $r \to 0$ as well, and such a limit is zero by the boundedness of $u$. By integrating the last inequality, we finally deduce

$$u'(r) \le \frac{\sup f}{\lambda(N-\gamma)} r^{1-\gamma}.$$

The regularity of $u$ at zero is then a consequence of the estimate $|u'(r)| \le cr^{1-\gamma}$. Elsewhere, it follows from the assumption $u \in C^2(B(0,1) \setminus \{0\})$. $\square$

<figure><img src="image_2.png" /></figure>
if $u''(r) \ge 0$, or

$$u''(r) + (\tilde{N}_- - 1) \frac{u'(r)}{r} \ge \frac{f(r)r^{-\gamma}}{\lambda}$$

if $u''(r) \le 0$. Since $u'(r) \ge 0$ and $\tilde{N}_- \ge N$, in both cases one has

$$u''(r) + (\tilde{N}_- - 1) \frac{u'(r)}{r} \ge \frac{f(r)r^{-\gamma}}{\Lambda},$$

that is

$$(u'r^{\tilde{N}_- - 1})' \ge \frac{f(r)r^{\tilde{N}_- - 1 - \gamma}}{\Lambda} \ge \frac{\inf f}{\Lambda} r^{\tilde{N}_- - 1 - \gamma}.$$

Arguing as above, we deduce that $u'r^{\tilde{N}_- - 1}$ is increasing in a right neighborhood of 0, hence it has a nonnegative limit as $r \to 0$, and such a limit must be 0, since $u$ is bounded. Moreover, by integrating the above inequality, we obtain

$$u'(r) \ge \frac{\inf f}{\Lambda(\tilde{N}_- - \gamma)} r^{1-\gamma}.$$

Let us finally prove (iii). Assuming that $M^+(D^2u) = f(r)r^{-\gamma}$ and using statement (i), it follows that, for every $r > 0$ sufficiently small, one has either

$$u''(r) + (N - 1) \frac{u'(r)}{r} = \frac{f(r)r^{-\gamma}}{\Lambda}$$

or

$$u''(r) + (\tilde{N}_- - 1) \frac{u'(r)}{r} = \frac{f(r)r^{-\gamma}}{\lambda}.$$

In both cases, we deduce

$$u''(r) + (N - 1) \frac{u'(r)}{r} \le \frac{f(r)r^{-\gamma}}{\lambda},$$

which yields

$$(u'(r)r^{N-1})' \le \frac{f(r)}{\lambda}r^{N-1-\gamma} \le \frac{\sup f}{\lambda}r^{N-1-\gamma}.$$

Hence, $u'(r)r^{N-1} - \frac{\sup f}{\lambda(N-\gamma)}r^{N-\gamma}$ is non increasing in a right neighborhood of 0 and it has a limit as $r \to 0$. This implies that $u'(r)r^{N-1}$ has a limit as $r \to 0$ as well, and such a limit is zero by the boundedness of $u$. By integrating the last inequality, we finally deduce

$$u'(r) \le \frac{\sup f}{\lambda(N-\gamma)} r^{1-\gamma}.$$

The regularity of $u$ at zero is then a consequence of the estimate $|u'(r)| \le cr^{1-\gamma}$. Elsewhere, it follows from the assumption $u \in C^2(B(0,1) \setminus \{0\})$. $\square$

<figure><img src="image_1.png" /></figure>
**Remark 2.2.** By using the change of variable $v = -u$, one gets that if $f$ belongs to $C(B(0,1) \setminus \{0\})$, it is radial, bounded and positive and $u \in C^2(B(0,1) \setminus \{0\})$ is a bounded radial function satisfying

$$ \mathcal{M}^{-}(D^{2}u) \leq -f r^{-\gamma} \quad \text{in } B(0,1) \setminus \{0\}, $$

then, for $r$ sufficiently small, $u'(r) \le 0$, $\lim_{r \to 0} u'(r) r^{\tilde{N}-1} = 0$ and

$$ u'(r) \leq - \frac{\inf f}{\Lambda(\tilde{N}_- - \gamma)} r^{1-\gamma}. $$

Moreover, if $\mathcal{M}^{-}(D^{2}u) = -f r^{-\gamma}$ in $B(0,1) \setminus \{0\}$, then $|u'(r)| \le cr^{1-\gamma}$ for a positive constant $c$. Hence $u$ is locally Lipschitz continuous in $B(0,1)$ for $\gamma \le 1$, it belongs to $C^{1}(B(0,1))$ if $\gamma < 1$, and it is locally Hölder continuous in $B(0,1)$ with exponent $2 - \gamma$ for $\gamma > 1$. Obviously, since $\mathcal{M}^{-} \le \mathcal{M}^{+}$, one gets an analogous conclusion when

$$ \mathcal{M}^{+}(D^{2}u) \leq -f r^{-\gamma}. $$

We can now prove a comparison principle for general radial fully nonlinear singular equations of the form

$$ F(D^2u) - \beta u r^{-\gamma} = f(r)r^{-\gamma} \quad \text{in } B(0,1) \setminus \{0\} $$

when no boundary condition at the origin is assumed.

**Theorem 2.3.** Let $f, g \in C(B(0,1))$ be radial functions and assume that $u, v \in C(\overline{B(0,1)}) \cap C^2(B(0,1) \setminus \{0\})$ are radial functions satisfying in $B(0,1) \setminus \{0\}$

$$ F(D^2u) - \beta u(r)r^{-\gamma} \geq f(r)r^{-\gamma} $$

$$ F(D^2v) - \beta v(r)r^{-\gamma} \leq g(r)r^{-\gamma} $$

with $\beta \ge 0$ and $f \ge g$ in $B(0,1)$. Then, $u \le v$ on $\partial B(0,1)$ implies $u \le v$ in $\overline{B(0,1)}$.

*Proof.* Let us first consider the case in which either $\beta > 0$ or $f > g$ in $B(0,1)$. We suppose by contradiction that

$$ \max_{B(0,1)}(u-v) > 0. $$

If the maximum is achieved at 0, then $(u-v)(0) > 0$ and, by the assumptions on $f, g$ and $\beta$, there exist $\delta > 0$ and a neighborhood on the right of 0 on which

$$ (\beta(u(r) - v(r)) + f(r) - g(r))r^{-\gamma} \geq \delta r^{-\gamma} > 0 $$

By (1.8), we then obtain, for $r$ sufficiently small,

$$ \mathcal{M}^{+}(D^{2}(u-v)) \geq \delta r^{-\gamma} $$
**Remark 2.2.** By using the change of variable $v = -u$, one gets that if $f$ belongs to $C(B(0,1) \setminus \{0\})$, it is radial, bounded and positive and $u \in C^2(B(0,1) \setminus \{0\})$ is a bounded radial function satisfying

$$ \mathcal{M}^{-}(D^{2}u) \leq -f r^{-\gamma} \quad \text{in } B(0,1) \setminus \{0\}, $$

then, for $r$ sufficiently small, $u'(r) \le 0$, $\lim_{r \to 0} u'(r) r^{\tilde{N}-1} = 0$ and

$$ u'(r) \leq - \frac{\inf f}{\Lambda(\tilde{N}_- - \gamma)} r^{1-\gamma}. $$

Moreover, if $\mathcal{M}^{-}(D^{2}u) = -f r^{-\gamma}$ in $B(0,1) \setminus \{0\}$, then $|u'(r)| \le cr^{1-\gamma}$ for a positive constant $c$. Hence $u$ is locally Lipschitz continuous in $B(0,1)$ for $\gamma \le 1$, it belongs to $C^{1}(B(0,1))$ if $\gamma < 1$, and it is locally Hölder continuous in $B(0,1)$ with exponent $2 - \gamma$ for $\gamma > 1$. Obviously, since $\mathcal{M}^{-} \le \mathcal{M}^{+}$, one gets an analogous conclusion when

$$ \mathcal{M}^{+}(D^{2}u) \leq -f r^{-\gamma}. $$

We can now prove a comparison principle for general radial fully nonlinear singular equations of the form

$$ F(D^2u) - \beta u r^{-\gamma} = f(r)r^{-\gamma} \quad \text{in } B(0,1) \setminus \{0\} $$

when no boundary condition at the origin is assumed.

**Theorem 2.3.** Let $f, g \in C(B(0,1))$ be radial functions and assume that $u, v \in C(\overline{B(0,1)}) \cap C^2(B(0,1) \setminus \{0\})$ are radial functions satisfying in $B(0,1) \setminus \{0\}$

$$ F(D^2u) - \beta u(r)r^{-\gamma} \geq f(r)r^{-\gamma} $$

$$ F(D^2v) - \beta v(r) r^{-\gamma} \leq g(r) r^{-\gamma} $$

with $\beta \ge 0$ and $f \ge g$ in $B(0,1)$. Then, $u \le v$ on $\partial B(0,1)$ implies $u \le v$ in $\overline{B(0,1)}$.

*Proof.* Let us first consider the case in which either $\beta > 0$ or $f > g$ in $B(0,1)$. We suppose by contradiction that

$$ \max_{B(0,1)}(u-v) > 0. $$

If the maximum is achieved at 0, then $(u-v)(0) > 0$ and, by the assumptions on $f, g$ and $\beta$, there exist $\delta > 0$ and a neighborhood on the right of 0 on which

$$ (\beta(u(r) - v(r)) + f(r) - g(r))r^{-\gamma} \geq \delta r^{-\gamma} > 0 $$

By (1.8), we then obtain, for $r$ sufficiently small,

$$ \mathcal{M}^{+}(D^{2}(u-v)) \geq \delta r^{-\gamma} $$
Using Lemma 2.1 for $u - v$, one gets that for some positive constant $c$, $(u - v)' \ge c\delta r^{1-\gamma}$, which contradicts the fact that $u - v$ attains its maximum at 0. Hence, there exists $0 < \bar{r} < 1$ such that $u(\bar{r}) - v(\bar{r}) = \max(u - v)$. Then, $(D^2u - D^2v)(\bar{r}) \le 0$ and, by ellipticity, we get the contradiction

$$ (f(\bar{r}) + \beta u(\bar{r}))r^{-\gamma} \le F(D^2u(\bar{r})) \le F(D^2v(\bar{r})) \le (g(\bar{r}) + \beta v(\bar{r}))r^{-\gamma}. $$

For the case $\beta = 0$ and $f \ge g$, let us introduce the radial function

$$ w(r) = 1 - r^{\tau} $$

with $0 < \tau \le 2 - \gamma$. A direct computation shows that

$$ \mathcal{M}^{+}(D^{2}w) \le \tau\Lambda(|\tau - 1| - (\tilde{N}_{+} - 1))r^{\tau-2}. $$

We observe that $|\tau - 1| < 1 < \tilde{N}_{+} - 1$, so that

$$ \mathcal{M}^{+}(D^{2}w) \le -Cr^{-\gamma} \quad \text{in } B(0,1) \setminus \{0\} $$

with $C = \tau\Lambda(\tilde{N}_{+} - 1 - |\tau - 1|) > 0$. Thus, for any $\epsilon > 0$, we have

$$ F(D^2(u - \epsilon w)) \ge F(D^2u) - \epsilon \mathcal{M}^+(D^2w) \ge (f + \epsilon C)r^{-\gamma}. $$

Since $f + \epsilon C > g$ and $u - \epsilon w = u \le v$ on $\partial B(0, 1)$, the previous argument proves that

$$ u - \epsilon w \le v \quad \text{in } \overline{B(0,1)} $$

and the conclusion follows by letting $\epsilon \to 0$. $\square$

**Remark 2.4.** The auxiliary function introduced in the proof of Theorem 2.3 shows that there exist a radial function $w \in C^2(B(0, 1) \setminus \{0\})$, strictly positive in $B(0, 1) \setminus \{0\}$, such that

$$ F(D^2w) \le -cwr^{-\gamma} \quad \text{in } B(0,1) \setminus \{0\} $$

for a constant $c > 0$. This proves that $\bar{\lambda}_{\gamma}(F) \ge \bar{\lambda}'_{\gamma}(F) \ge c > 0$.

Next, we have the following existence, uniqueness and regularity result.

**Theorem 2.5.** Let $f \in C(B(0, 1))$ be a radial, bounded function. For $\beta \ge 0$ and $b \in \mathbb{R}$ there exists a unique bounded radial function $u \in C(\overline{B}(0, 1) \setminus \{0\}) \cap C^2(B(0, 1) \setminus \{0\})$ satisfying

$$ \left\{ \begin{array}{ll} F(D^2u) - \beta u r^{-\gamma} = r^{-\gamma} f(r) & \text{in } B(0,1) \setminus \{0\} \\ u = b & \text{on } \partial B(0,1) \end{array} \right. \qquad (2.2) $$

Moreover, $u$ can be extended up to $\overline{B}(0, 1)$, and one has: $u \in C^1(\overline{B}(0, 1))$ if $\gamma < 1$, $u$ is Lipschitz continuous in $\overline{B}(0, 1)$ if $\gamma \le 1$, $u$ is Hölder continuous in $\overline{B}(0, 1)$ with exponent $2 - \gamma$ if $\gamma > 1$.
Using Lemma 2.1 for $u-v$, one gets that for some positive constant $c$, $(u-v)' \ge c\delta r^{1-\gamma}$, which contradicts the fact that $u-v$ attains its maximum at 0. Hence, there exists $0 < \bar{r} < 1$ such that $u(\bar{r}) - v(\bar{r}) = \max(u-v)$. Then, $(D^2u - D^2v)(\bar{r}) \le 0$ and, by ellipticity, we get the contradiction

$$ (f(\bar{r}) + \beta u(\bar{r}))r^{-\gamma} \le F(D^2u(\bar{r})) \le F(D^2v(\bar{r})) \le (g(\bar{r}) + \beta v(\bar{r}))r^{-\gamma}. $$

For the case $\beta = 0$ and $f \ge g$, let us introduce the radial function

$$ w(r) = 1 - r^{\tau} $$

with $0 < \tau \le 2 - \gamma$. A direct computation shows that

$$ \mathcal{M}^{+}(D^{2}w) \le \tau\Lambda(|\tau - 1| - (\tilde{N}_{+} - 1))r^{\tau-2}. $$

We observe that $|\tau - 1| < 1 < \tilde{N}_{+} - 1$, so that

$$ \mathcal{M}^{+}(D^{2}w) \le -Cr^{-\gamma} \quad \text{in } B(0,1) \setminus \{0\} $$

with $C = \tau\Lambda(\tilde{N}_{+} - 1 - |\tau - 1|) > 0$. Thus, for any $\epsilon > 0$, we have

$$ F(D^2(u - \epsilon w)) \ge F(D^2u) - \epsilon \mathcal{M}^+(D^2w) \ge (f + \epsilon C)r^{-\gamma}. $$

Since $f + \epsilon C > g$ and $u - \epsilon w = u \le v$ on $\partial B(0, 1)$, the previous argument proves that

$$ u - \epsilon w \le v \quad \text{in } \overline{B(0,1)} $$

and the conclusion follows by letting $\epsilon \to 0$. $\square$

**Remark 2.4.** The auxiliary function introduced in the proof of Theorem 2.3 shows that there exist a radial function $w \in C^2(B(0,1) \setminus \{0\})$, strictly positive in $B(0,1) \setminus \{0\}$, such that

$$ F(D^2w) \le -cwr^{-\gamma} \quad \text{in } B(0,1) \setminus \{0\} $$

for a constant $c > 0$. This proves that $\bar{\lambda}_{\gamma}(F) \ge \bar{\lambda}'_{\gamma}(F) \ge c > 0$.

Next, we have the following existence, uniqueness and regularity result.

**Theorem 2.5.** Let $f \in C(B(0,1))$ be a radial, bounded function. For $\beta \ge 0$ and $b \in \mathbb{R}$ there exists a unique bounded radial function $u \in C(\overline{B}(0,1) \setminus \{0\}) \cap C^2(B(0,1) \setminus \{0\})$ satisfying

$$ \left\{ \begin{array}{ll} F(D^2u) - \beta u r^{-\gamma} = r^{-\gamma} f(r) & \text{in } B(0,1) \setminus \{0\} \\ u = b & \text{on } \partial B(0,1) \end{array} \right. \qquad (2.2) $$

Moreover, $u$ can be extended up to $\overline{B}(0,1)$, and one has: $u \in C^1(\overline{B}(0,1))$ if $\gamma < 1$, $u$ is Lipschitz continuous in $\overline{B}(0,1)$ if $\gamma \le 1$, $u$ is Hölder continuous in $\overline{B}(0,1)$ with exponent $2 - \gamma$ if $\gamma > 1$.
Proof. For every $n \in \mathbb{N}$ let us introduce the regularized Dirichlet boundary value problem

$$ \begin{cases} F(D^2 u_n) - \beta u_n (r^2 + 1/n)^{-\gamma/2} = (r^2 + 1/n)^{-\gamma/2} f(r) & \text{in } B(0,1) \\ u_n = b & \text{on } \partial B(0,1) \end{cases} \quad (2.3) $$

which, by standard viscosity solutions theory, see [20], has a unique solution $u_n \in C(\overline{B(0,1)})$. By the symmetry results of [16], it follows that $u_n$ is radial, hence, as a solution of the associated ODE, $u_n$ belongs to $C^2(B(0,1))$.

For $0 < \tau \le 2 - \gamma$, let us consider the radial function

$$ w(r) = L(1 - r^{\tau}) + b, $$

where $L > 0$ is a constant to be suitably chosen. The same computation used in the proof of Theorem 2.3 yields

$$ \mathcal{M}^{+}(D^2 w) \le -L C r^{-\gamma} \le -L C (r^2 + 1/n)^{-\gamma/2}, $$

and then, by uniform ellipticity, it follows that

$$ F(D^2 w) - \beta w r^{-\gamma} \le (-L C - \beta b)(r^2 + 1/n)^{-\gamma/2} \le f(r)(r^2 + 1/n)^{-\gamma/2} $$

as soon as $C$ is chosen large enough. Analogously, for some convenient positive constant $L$, the function $-L(1 - r^{\tau}) + b$ is a radial sub-solution of problem (2.3).

The standard comparison principle then implies that the sequence $\{u_n\}$ is uniformly bounded in $C(\overline{B(0,1)})$. Hence, it is locally uniformly bounded in $C^2(B(0,1) \setminus \{0\})$ and, up to a subsequence, it is converging locally uniformly in $\overline{B(0,1)} \setminus \{0\}$ to a radial solution $u$ of problem (2.2), which is a globally bounded function belonging to $C^2(B(0,1) \setminus \{0\})$.

Let us now show that the constructed bounded radial solution $u$ is actually continuous in the whole ball $\overline{B(0,1)}$. Indeed, the same approximation argument used in order to prove the existence of $u$ can be applied, in particular, in order to show the existence of a radial bounded solution $\bar{w}$ of the Dirichlet problem

$$ \begin{cases} \mathcal{M}^{+}(D^2 \bar{w}) = -(\|f\|_{\infty} + B + 1)r^{-\gamma} & \text{in } B(0,1) \setminus \{0\} \\ \bar{w} = u & \text{on } \partial B(0,1) \end{cases} $$

where $B > 0$ is a constant such that $\beta|u| \le B$ in $B(0,1)$. Lemma 2.1 and Remark 2.2 applied to $\bar{w}$ yield that

$$ |\bar{w}'| \le cr^{1-\gamma} $$

for some $c > 0$. As a consequence, we have

$$ \mathcal{M}^{+}(D^2(u - \bar{w})) \ge F(D^2(u - \bar{w})) \ge F(D^2 u) - \mathcal{M}^{+}(D^2 \bar{w}) \ge r^{-\gamma} \quad \text{in } B(0,1) \setminus \{0\} $$
Proof. For every $n \in \mathbb{N}$ let us introduce the regularized Dirichlet boundary value problem

$$ \begin{cases} F(D^2 u_n) - \beta u_n (r^2 + 1/n)^{-\gamma/2} = (r^2 + 1/n)^{-\gamma/2} f(r) & \text{in } B(0,1) \\ u_n = b & \text{on } \partial B(0,1) \end{cases} \quad (2.3) $$

which, by standard viscosity solutions theory, see [20], has a unique solution $u_n \in C(\overline{B(0,1)})$. By the symmetry results of [16], it follows that $u_n$ is radial, hence, as a solution of the associated ODE, $u_n$ belongs to $C^2(B(0,1))$.

For $0 < \tau \le 2 - \gamma$, let us consider the radial function

$$ w(r) = L(1 - r^{\tau}) + b, $$

where $L > 0$ is a constant to be suitably chosen. The same computation used in the proof of Theorem 2.3 yields

$$ \mathcal{M}^{+}(D^{2}w) \leq -L C r^{-\gamma} \leq -L C (r^{2} + 1/n)^{-\gamma/2}, $$

and then, by uniform ellipticity, it follows that

$$ F(D^2 w) - \beta w r^{-\gamma} \leq (-L C - \beta b)(r^2 + 1/n)^{-\gamma/2} \leq f(r)(r^2 + 1/n)^{-\gamma/2} $$

as soon as $C$ is chosen large enough. Analogously, for some convenient positive constant $L$, the function $-L(1 - r^{\tau}) + b$ is a radial sub-solution of problem (2.3).

The standard comparison principle then implies that the sequence $\{u_n\}$ is uniformly bounded in $C(\overline{B}(0,1))$. Hence, it is locally uniformly bounded in $C^2(B(0,1) \setminus \{0\})$ and, up to a subsequence, it is converging locally uniformly in $\overline{B}(0,1) \setminus \{0\}$ to a radial solution $u$ of problem (2.2), which is a globally bounded function belonging to $C^2(B(0,1) \setminus \{0\})$.

Let us now show that the constructed bounded radial solution $u$ is actually continuous in the whole ball $\overline{B}(0,1)$. Indeed, the same approximation argument used in order to prove the existence of $u$ can be applied, in particular, in order to show the existence of a radial bounded solution $\bar{w}$ of the Dirichlet problem

$$ \begin{cases} \mathcal{M}^{+}(D^{2}\bar{w}) = -(\|f\|_{\infty} + B + 1)r^{-\gamma} & \text{in } B(0,1) \setminus \{0\} \\ \bar{w} = u & \text{on } \partial B(0,1) \end{cases} $$

where $B > 0$ is a constant such that $\beta|u| \le B$ in $B(0,1)$. Lemma 2.1 and Remark 2.2 applied to $\bar{w}$ yield that

$$ |\bar{w}'| \le cr^{1-\gamma} $$

for some $c > 0$. As a consequence, we have

$$ \mathcal{M}^{+}(D^{2}(u - \bar{w})) \geq F(D^{2}(u - \bar{w})) \geq F(D^{2}u) - \mathcal{M}^{+}(D^{2}\bar{w}) \geq r^{-\gamma} \quad \text{in } B(0,1) \setminus \{0\} $$
and then, by Lemma 2.1 (i), in a right neighborhood of zero one has $(u - \bar{w})'(r) \ge 0$. Hence,

$$u'(r) \ge -cr^{1-\gamma}$$

for $r$ small enough. Analogously, we have

$$\mathcal{M}^{-}(D^{2}(u + \bar{w})) \le F(D^{2}(u + \bar{w})) \le F(D^{2}u) + \mathcal{M}^{+}(D^{2}\bar{w}) \le -r^{-\gamma} \quad \text{in } B(0, 1) \setminus \{0\}$$

which implies, by Remark 2.2,

$$u'(r) \le -\bar{w}' \le cr^{1-\gamma}$$

for $r$ small enough. Arguing as in the proof of Lemma 2.1, from the estimate $|u'(r)| \le cr^{1-\gamma}$ for $r$ sufficiently small, we deduce that $u$ is Lipschitz continuous in $\overline{B}(0, 1)$ if $\gamma \le 1$, it belongs to $C^1(\overline{B}(0, 1))$ if $\gamma < 1$, and it is Hölder continuous in $\overline{B}(0, 1)$ with exponent $2 - \gamma$ if $\gamma > 1$. Let us observe that the argument above shows that any bounded radial solution of problem (2.2) is continuous in $\overline{B}(0, 1)$. This, jointly with Theorem 2.3, implies that problem (2.2) has a unique radial bounded solution. $\square$

The argument used in the above proof yields also the following compactness result, which we state separately for the sake of clarity.

**Theorem 2.6.** Let $\{u_n\}_n$ be a uniformly bounded sequence of radial functions belonging to $C^2(B(0, 1) \setminus \{0\})$ and satisfying

$$F(D^2 u_n) = f_n r^{-\gamma} \quad \text{in } B(0, 1) \setminus \{0\},$$

where $\{f_n\}_n$ are radial, bounded and continuous on $B(0, 1) \setminus \{0\}$. If $\{f_n\}$ is uniformly bounded, then $\{u_n\}_n$ are equicontinuous, thus uniformly converging in $\overline{B}(0, 1)$ up to a subsequence. If $\{f_n\}$ is uniformly converging to $f \in C(\overline{B}(0, 1))$, then, up to a subsequence, $\{u_n\}_n$ is uniformly converging to a radial solution $u \in C^2(B(0, 1) \setminus \{0\}) \cap C(\overline{B}(0, 1))$ of

$$F(D^2 u) = r^{-\gamma} f \quad \text{in } B(0, 1) \setminus \{0\}.$$

In the next results, we prove several properties of the “smooth” eigenvalue

$$\bar{\lambda}'_\gamma := \sup\{\mu : \exists u \in C^2(B(0, 1)\setminus\{0\}), u > 0 \text{ in } B(0, 1)\setminus\{0\}, u \text{ radial}, F(D^2u) + \mu \frac{u}{r^\gamma} \le 0\}.$$

In Subsection 2.3 we will prove in fact that $\bar{\lambda}'_\gamma = \bar{\lambda}_\gamma$.

Let us start by proving the validity of the maximum principle below the value $\bar{\lambda}'_\gamma$.

**Theorem 2.7.** Let $\mu < \bar{\lambda}'_\gamma$ and suppose that $u \in C(\overline{B(0, 1)}) \cap C^2(B(0, 1) \setminus \{0\})$ is a radial function satisfying

$$F(D^2 u) + \mu u r^{-\gamma} \ge 0 \quad \text{in } B(0, 1) \setminus \{0\}.$$

If $u(1) \le 0$, then $u \le 0$ in $\overline{B}(0, 1)$.
and then, by Lemma 2.1 (i), in a right neighborhood of zero one has $(u - \bar{w})'(r) \ge 0$. Hence,

$$u'(r) \ge -cr^{1-\gamma}$$

for $r$ small enough. Analogously, we have

$$\mathcal{M}^{-}(D^{2}(u + \bar{w})) \le F(D^{2}(u + \bar{w})) \le F(D^{2}u) + \mathcal{M}^{+}(D^{2}\bar{w}) \le -r^{-\gamma} \quad \text{in } B(0, 1) \setminus \{0\}$$

which implies, by Remark 2.2,

$$u'(r) \le -\bar{w}' \le cr^{1-\gamma}$$

for $r$ small enough. Arguing as in the proof of Lemma 2.1, from the estimate $|u'(r)| \le cr^{1-\gamma}$ for $r$ sufficiently small, we deduce that $u$ is Lipschitz continuous in $\overline{B}(0, 1)$ if $\gamma \le 1$, it belongs to $C^1(\overline{B}(0, 1))$ if $\gamma < 1$, and it is Hölder continuous in $\overline{B}(0, 1)$ with exponent $2 - \gamma$ if $\gamma > 1$. Let us observe that the argument above shows that any bounded radial solution of problem (2.2) is continuous in $\overline{B}(0, 1)$. This, jointly with Theorem 2.3, implies that problem (2.2) has a unique radial bounded solution. $\square$

The argument used in the above proof yields also the following compactness result, which we state separately for the sake of clarity.

**Theorem 2.6.** Let $\{u_n\}_n$ be a uniformly bounded sequence of radial functions belonging to $C^2(B(0, 1) \setminus \{0\})$ and satisfying

$$F(D^2 u_n) = f_n r^{-\gamma} \quad \text{in } B(0, 1) \setminus \{0\},$$

where $\{f_n\}_n$ are radial, bounded and continuous on $B(0, 1) \setminus \{0\}$. If $\{f_n\}$ is uniformly bounded, then $\{u_n\}_n$ are equicontinuous, thus uniformly converging in $\overline{B}(0, 1)$ up to a subsequence. If $\{f_n\}$ is uniformly converging to $f \in C(\overline{B}(0, 1))$, then, up to a subsequence, $\{u_n\}_n$ is uniformly converging to a radial solution $u \in C^2(B(0, 1) \setminus \{0\}) \cap C(\overline{B}(0, 1))$ of

$$F(D^2 u) = r^{-\gamma} f \quad \text{in } B(0, 1) \setminus \{0\}.$$

In the next results, we prove several properties of the “smooth” eigenvalue

$$\bar{\lambda}'_\gamma := \sup\{\mu : \exists u \in C^2(B(0, 1)\setminus\{0\}), u > 0 \text{ in } B(0, 1)\setminus\{0\}, u \text{ radial, } F(D^2u) + \mu \frac{u}{r^\gamma} \le 0\}.$$

In Subsection 2.3 we will prove in fact that $\bar{\lambda}'_\gamma = \bar{\lambda}_\gamma$.

Let us start by proving the validity of the maximum principle below the value $\bar{\lambda}'_\gamma$.

**Theorem 2.7.** Let $\mu < \bar{\lambda}'_\gamma$ and suppose that $u \in C(\overline{B(0, 1)}) \cap C^2(B(0, 1) \setminus \{0\})$ is a radial function satisfying

$$F(D^2 u) + \mu u r^{-\gamma} \ge 0 \quad \text{in } B(0, 1) \setminus \{0\}.$$

If $u(1) \le 0$, then $u \le 0$ in $\overline{B}(0, 1)$.
Proof. If $\mu < 0$, we just apply Theorem 2.3 with $v \equiv g \equiv f \equiv 0$. So, we can assume without loss of generality that $\mu \ge 0$.

For $\mu' \in ]\mu, \bar{\lambda}'_{\gamma}[$, let $v \in C^2(B(0, 1) \setminus \{0\})$ be a radial function satisfying

$$F(D^2v) + \mu'vr^{-\gamma} \le 0, \quad v > 0 \quad \text{in } B(0, 1) \setminus \{0\}.$$

We can assume without loss of generality that $v > 0$ on $\partial B(0, 1)$, e.g. by performing a dilation in $r$ (this may change a little $\mu'$ but we can still suppose by continuity that $\mu' \in ]\mu, \bar{\lambda}'_{\gamma}[$). Arguing as in the proof of Lemma 2.1, since $M^-(D^2v) + \mu'vr^{-\gamma} \le 0$, we easily obtain that $v'(r)$ has constant sign near zero. Assuming by contradiction that $v'(r) \ge 0$ in a neighborhood of zero, then $v$ is bounded in $B(0, 1) \setminus \{0\}$. Hence, Remark 2.2 applies and yields $v'(r) \le 0$ for $r$ small enough: a contradiction. This shows that $v'(r) \le 0$ in a neighborhood of 0.

Then, there are two possible cases: either $\lim_{r \to 0} v(r) = +\infty$ or $v$ can be extended as a continuous function on $\bar{B}(0, 1)$. In the first case, by applying the standard comparison principle, it is easy to prove that for all $\epsilon > 0$ one has $u \le \epsilon v$ in $\bar{B}(0, 1) \setminus \{0\}$. In this case, letting $\epsilon \to 0$, we get the conclusion. On the other hand, if $v$ is bounded and continuous on $\bar{B}(0, 1)$, we can argue by contradiction. Let us assume that $u$ is positive somewhere in $B(0, 1)$, so that $\frac{u}{v}$ has a positive maximum on $\bar{B}(0, 1)$, achieved at some point inside $B(0, 1)$. Up to a multiplicative constant for $v$, we can suppose that

$$\max_{B(0,1)} \frac{u}{v} = 1,$$

so that $u(r) \le v(r)$. If the maximum is achieved at 0, then one has $u(0) = v(0) > 0$. By continuity, for $r$ small enough one has

$$M^+(D^2(u-v)) \ge F(D^2u) - F(D^2v) \ge \frac{1}{2}(\mu'v(0) - \mu u(0))r^{-\gamma}.$$

Since $\mu'v(0) - \mu u(0) = (\mu' - \mu)v(0) > 0$, we can use Lemma 2.1 (ii) and we get that $(u-v)'(r) > 0$ in a right neighborhood of 0. This is a contradiction to $u - v$ has a maximum point at zero. Hence, we have that $1 = \max \frac{u}{v} > \frac{u(0)}{v(0)}$. Let us select $\eta < 1$ such that $\eta > \max\{\frac{\mu}{\mu'}, \frac{u(0)}{v(0)}\}$. Then, the function $u - \eta v$ has a positive maximum achieved at some point $0 < \bar{r} < 1$.

Since $D^2u(\bar{r}) \le \eta D^2v(\bar{r})$, by ellipticity and using (1.6), we get

$$-\mu v(\bar{r}) \bar{r}^{-\gamma} \le -\mu u(\bar{r}) \bar{r}^{-\gamma} \le F(D^2u(\bar{r})) \le F(\eta D^2v(\bar{r})) \le -\mu' \eta v(\bar{r}) \bar{r}^{-\gamma},$$

which gives the contradiction $\mu \ge \mu'\eta$.

□

The next result provides the existence, uniqueness and regularity of solutions below the eigenvalue $\bar{\lambda}'_{\gamma}$.
Proof. If $\mu < 0$, we just apply Theorem 2.3 with $v \equiv g \equiv f \equiv 0$. So, we can assume without loss of generality that $\mu \ge 0$.

For $\mu' \in ]\mu, \bar{\lambda}'_{\gamma}[$, let $v \in C^2(B(0, 1) \setminus \{0\})$ be a radial function satisfying

$$F(D^2v) + \mu'vr^{-\gamma} \le 0, \quad v > 0 \quad \text{in } B(0, 1) \setminus \{0\}.$$

We can assume without loss of generality that $v > 0$ on $\partial B(0, 1)$, e.g. by performing a dilation in $r$ (this may change a little $\mu'$ but we can still suppose by continuity that $\mu' \in ]\mu, \bar{\lambda}'_{\gamma}[$). Arguing as in the proof of Lemma 2.1, since $M^-(D^2v) + \mu'vr^{-\gamma} \le 0$, we easily obtain that $v'(r)$ has constant sign near zero. Assuming by contradiction that $v'(r) \ge 0$ in a neighborhood of zero, then $v$ is bounded in $B(0, 1) \setminus \{0\}$. Hence, Remark 2.2 applies and yields $v'(r) \le 0$ for $r$ small enough: a contradiction. This shows that $v'(r) \le 0$ in a neighborhood of 0.

Then, there are two possible cases: either $\lim_{r \to 0} v(r) = +\infty$ or $v$ can be extended as a continuous function on $\bar{B}(0, 1)$. In the first case, by applying the standard comparison principle, it is easy to prove that for all $\epsilon > 0$ one has $u \le \epsilon v$ in $\bar{B}(0, 1) \setminus \{0\}$. In this case, letting $\epsilon \to 0$, we get the conclusion. On the other hand, if $v$ is bounded and continuous on $\bar{B}(0, 1)$, we can argue by contradiction. Let us assume that $u$ is positive somewhere in $B(0, 1)$, so that $\frac{u}{v}$ has a positive maximum on $\bar{B}(0, 1)$, achieved at some point inside $B(0, 1)$. Up to a multiplicative constant for $v$, we can suppose that

$$\max_{B(0,1)} \frac{u}{v} = 1,$$

so that $u(r) \le v(r)$. If the maximum is achieved at 0, then one has $u(0) = v(0) > 0$. By continuity, for $r$ small enough one has

$$M^+(D^2(u-v)) \ge F(D^2u) - F(D^2v) \ge \frac{1}{2}(\mu'v(0) - \mu u(0))r^{-\gamma}.$$

Since $\mu'v(0) - \mu u(0) = (\mu' - \mu)v(0) > 0$, we can use Lemma 2.1 (ii) and we get that $(u-v)'(r) > 0$ in a right neighborhood of 0. This is a contradiction to $u - v$ has a maximum point at zero. Hence, we have that $1 = \max \frac{u}{v} > \frac{u(0)}{v(0)}$. Let us select $\eta < 1$ such that $\eta > \max\{\frac{\mu}{\mu'}, \frac{u(0)}{v(0)}\}$. Then, the function $u - \eta v$ has a positive maximum achieved at some point $0 < \bar{r} < 1$.

Since $D^2u(\bar{r}) \le \eta D^2v(\bar{r})$, by ellipticity and using (1.6), we get

$$-\mu v(\bar{r}) \bar{r}^{-\gamma} \le -\mu u(\bar{r}) \bar{r}^{-\gamma} \le F(D^2u(\bar{r})) \le F(\eta D^2v(\bar{r})) \le -\mu' \eta v(\bar{r}) \bar{r}^{-\gamma},$$

which gives the contradiction $\mu \ge \mu'\eta$.

□

The next result provides the existence, uniqueness and regularity of solutions below the eigenvalue $\bar{\lambda}'_{\gamma}$.
**Theorem 2.8.** Let $f$ be a radial and continuous function in $\overline{B}(0,1)$ satisfying $f \le 0$, with $f$ not identically zero. Then, for every $\mu < \bar{\lambda}'_\gamma$ there exists a unique, bounded, radial function $u \in C^2(B(0,1) \setminus \{0\})$ satisfying

$$ \begin{cases} F(D^2u) + \mu u r^{-\gamma} = f(r)r^{-\gamma} & \text{in } B(0,1) \setminus \{0\} \\ u = 0 & \text{on } \partial B(0,1) \end{cases} $$

Moreover, $u$ can be extended as a strictly positive continuous function in $B(0,1)$, Lipschitz continuous in $\overline{B}(0,1)$ if $\gamma \le 1$, $(2-\gamma)-\text{Hölder continuous in } \overline{B}(0,1)$ if $\gamma > 1$.

*Proof.* As in the proof of Theorem 2.7, we can assume without loss of generality that $\mu > 0$, otherwise the conclusion just follows from Theorem 2.5.

Let us first show the existence of a bounded, radial, smooth solution $u$. We recursively define a sequence $\{u_n\}_{n \ge 0}$ as follows: we set

$$ u_0 = 0, $$

and then, by using Theorem 2.5, we define $u_{n+1} \in C^2(B(0,1) \setminus \{0\}) \cap C(\overline{B(0,1)})$ as the unique bounded radial solution of

$$ \begin{cases} F(D^2u_{n+1}) = (f - \mu u_n)r^{-\gamma} & \text{in } B(0,1) \setminus \{0\} \\ u_{n+1} = 0 & \text{on } \partial B(0,1) \end{cases} $$

By Theorem 2.3, we have that $u_{n+1} \ge 0$, hence it is strictly positive in $B(0,1) \setminus \{0\}$ by the standard strong maximum principle, since $f$ is not identically zero. In particular, $u_n$ is not identically zero for all $n \ge 1$. By applying the comparison principle in Theorem 2.3 again, we deduce also that $u_{n+1} \ge u_n$. Let us prove that $\{u_n\}_n$ is uniformly bounded. If not, by setting $v_n = \|u_n\|_{\infty}^{-1}u_n$ and $k_n = \|u_{n+1}\|_{\infty}^{-1}\|u_n\|_{\infty} \le 1$, one gets that $v_{n+1}$ satisfies

$$ \begin{cases} F(D^2v_{n+1}) = \left( \frac{f(r)}{\|u_{n+1}\|_{\infty}} - \mu k_n v_n(r) \right) r^{-\gamma} & \text{in } B(0,1) \setminus \{0\} \\ v_{n+1} = 0 & \text{on } \partial B(0,1) \end{cases} $$

Since $\{v_n\}_n$ is uniformly bounded, by applying Theorem 2.6, we can extract a subsequence still denoted by $\{v_n\}_n$ uniformly converging to a function $v \ge 0$ satisfying

$$ \begin{cases} F(D^2v) + \mu k v r^{-\gamma} = 0 & \text{in } B(0,1) \setminus \{0\} \\ v = 0 & \text{on } \partial B(0,1) \end{cases} $$

where $k \le 1$ is the limit of some converging subsequence of $\{k_n\}_n$. Since $v$ is a radial solution, one has that $v \in C^2(B(0,1) \setminus \{0\})$ and, since $\mu k \le \mu < \bar{\lambda}'_\gamma$, Theorem 2.7 yields $v \le 0$. Hence, we get $v \equiv 0$, a contradiction with $\|v\|_{\infty} = 1$.

We have obtained that $\{u_n\}_n$ is bounded, and using once more Theorem 2.6, we deduce that $\{u_n\}$ uniformly converges to some $u$, which satisfies the desired equation. By the strong
**Theorem 2.8.** Let $f$ be a radial and continuous function in $\overline{B}(0,1)$ satisfying $f \le 0$, with $f$ not identically zero. Then, for every $\mu < \bar{\lambda}'_\gamma$ there exists a unique, bounded, radial function $u \in C^2(B(0,1) \setminus \{0\})$ satisfying

$$ \begin{cases} F(D^2u) + \mu u r^{-\gamma} = f(r)r^{-\gamma} & \text{in } B(0,1) \setminus \{0\} \\ u = 0 & \text{on } \partial B(0,1) \end{cases} $$

Moreover, $u$ can be extended as a strictly positive continuous function in $B(0,1)$, Lipschitz continuous in $\overline{B}(0,1)$ if $\gamma \le 1$, $(2-\gamma)-\text{Hölder continuous in } \overline{B}(0,1)$ if $\gamma > 1$.

*Proof.* As in the proof of Theorem 2.7, we can assume without loss of generality that $\mu > 0$, otherwise the conclusion just follows from Theorem 2.5.

Let us first show the existence of a bounded, radial, smooth solution $u$. We recursively define a sequence $\{u_n\}_{n \ge 0}$ as follows: we set

$$ u_0 \equiv 0, $$

and then, by using Theorem 2.5, we define $u_{n+1} \in C^2(B(0,1) \setminus \{0\}) \cap C(\overline{B(0,1)})$ as the unique bounded radial solution of

$$ \begin{cases} F(D^2u_{n+1}) = (f - \mu u_n)r^{-\gamma} & \text{in } B(0,1) \setminus \{0\} \\ u_{n+1} = 0 & \text{on } \partial B(0,1) \end{cases} $$

By Theorem 2.3, we have that $u_{n+1} \ge 0$, hence it is strictly positive in $B(0,1) \setminus \{0\}$ by the standard strong maximum principle, since $f$ is not identically zero. In particular, $u_n$ is not identically zero for all $n \ge 1$. By applying the comparison principle in Theorem 2.3 again, we deduce also that $u_{n+1} \ge u_n$. Let us prove that $\{u_n\}_n$ is uniformly bounded. If not, by setting $v_n = \|u_n\|_{\infty}^{-1}u_n$ and $k_n = \|u_{n+1}\|_{\infty}^{-1}\|u_n\|_{\infty} \le 1$, one gets that $v_{n+1}$ satisfies

$$ \begin{cases} F(D^2v_{n+1}) = \left( \frac{f(r)}{\|u_{n+1}\|_{\infty}} - \mu k_n v_n(r) \right) r^{-\gamma} & \text{in } B(0,1) \setminus \{0\} \\ v_{n+1} = 0 & \text{on } \partial B(0,1) \end{cases} $$

Since $\{v_n\}_n$ is uniformly bounded, by applying Theorem 2.6, we can extract a subsequence still denoted by $\{v_n\}_n$ uniformly converging to a function $v \ge 0$ satisfying

$$ \begin{cases} F(D^2v) + \mu k v r^{-\gamma} = 0 & \text{in } B(0,1) \setminus \{0\} \\ v = 0 & \text{on } \partial B(0,1) \end{cases} $$

where $k \le 1$ is the limit of some converging subsequence of $\{k_n\}_n$. Since $v$ is a radial solution, one has that $v \in C^2(B(0,1) \setminus \{0\})$ and, since $\mu k \le \mu < \bar{\lambda}'_\gamma$, Theorem 2.7 yields $v \le 0$. Hence, we get $v \equiv 0$, a contradiction with $\|v\|_{\infty} = 1$.

We have obtained that $\{u_n\}_n$ is bounded, and using once more Theorem 2.6, we deduce that $\{u_n\}$ uniformly converges to some $u$, which satisfies the desired equation. By the strong
maximum principle, we get that $u > 0$ in $B(0, 1) \setminus \{0\}$. Moreover, by Remark 2.2, we have that $u'(r) \le 0$ for $r > 0$ small enough, which implies $u(0) > 0$. Finally, the global regularity of $u$ follows from Theorem 2.5.

As far as uniqueness is concerned, let $u$ and $v$ be two positive, bounded, radial solutions of

$$F(D^2u) + \mu ur^{-\gamma} = fr^{-\gamma},$$

and let us define

$$\eta = \sup \frac{u}{v}.$$

We assume, by contradiction, that $\eta > 1$. Then $\eta f \le f$, $\eta f \ne f$ and $\eta v - u \ge 0$, $\inf(\eta v - u) = 0$. If $\eta v - u$ achieves its minimum at $\bar{r} \ne 0$, then, by the equation, one has

$$\mathcal{M}^{-}(D^{2}(v\eta - u)) + \mu(\eta v - u)r^{-\gamma} = (\eta - 1)fr^{-\gamma} \le 0.$$

By the strong maximum principle, one gets that $v\eta - u = 0$ for every $r \ne 0$, which yields a contradiction with the equation since $\eta f \ne f$. Hence, we deduce that $\eta v - u > 0$ for any $r \ne 0$, and $\eta v(0) - u(0) = 0$. We then have

$$\mathcal{M}^{-}(D^{2}(\eta v - u)) < gr^{-\gamma},$$

with $g := u - \eta v + f < 0$ for any $r \ne 0$. By Lemma 2.1, we obtain $(\eta v - u)' \le 0$, a contradiction with the fact that $(\eta v - u)$ achieves its minimum at zero. Thus, we have obtained that $u \le v$ and, by exchanging the roles of $u$ and $v$, we conclude that $u \equiv v$. $\square$

We can now prove that the smooth eigenvalue $\bar{\lambda}'_\gamma$ is actually achieved on smooth eigenfunctions.

**Theorem 2.9.** There exists $u \in C(\overline{B(0, 1)}) \cap C^2(B(0, 1) \setminus \{0\})$, radial, strictly positive in $B(0, 1)$ and satisfying

$$\begin{cases} F(D^2u) + \bar{\lambda}'_\gamma ur^{-\gamma} = 0 & \text{in } B(0, 1) \setminus \{0\} \\ u = 0 & \text{on } \partial B(0, 1) \end{cases}$$

Furthermore, in $\overline{B(0, 1)}$, $u$ is Lipschitz continuous when $\gamma \le 1$ and Hölder continuous with exponent $2 - \gamma$ if $\gamma > 1$.

*Proof.* We consider a sequence $\{\lambda_n\}$, with $\lambda_n \to \bar{\lambda}'_\gamma$ and $\lambda_n < \bar{\lambda}'_\gamma$ and, for all $n$, the solution $u_n \in C(\overline{B(0, 1)}) \cap C^2(B(0, 1) \setminus \{0\})$ provided by Theorem 2.8 of

$$F(D^2u_n) + \lambda_n u_n r^{-\gamma} = -r^{-\gamma}, \quad u_n(1) = 0.$$
maximum principle, we get that $u > 0$ in $B(0, 1) \setminus \{0\}$. Moreover, by Remark 2.2, we have that $u'(r) \le 0$ for $r > 0$ small enough, which implies $u(0) > 0$. Finally, the global regularity of $u$ follows from Theorem 2.5.

As far as uniqueness is concerned, let $u$ and $v$ be two positive, bounded, radial solutions of

$$F(D^2u) + \mu ur^{-\gamma} = fr^{-\gamma},$$

and let us define

$$\eta = \sup \frac{u}{v}.$$

We assume, by contradiction, that $\eta > 1$. Then $\eta f \le f$, $\eta f \ne f$ and $\eta v - u \ge 0$, $\inf(\eta v - u) = 0$. If $\eta v - u$ achieves its minimum at $\bar{r} \ne 0$, then, by the equation, one has

$$\mathcal{M}^{-}(D^{2}(v\eta - u)) + \mu(\eta v - u)r^{-\gamma} = (\eta - 1)fr^{-\gamma} \le 0.$$

By the strong maximum principle, one gets that $v\eta - u = 0$ for every $r \ne 0$, which yields a contradiction with the equation since $\eta f \ne f$. Hence, we deduce that $\eta v - u > 0$ for any $r \ne 0$, and $\eta v(0) - u(0) = 0$. We then have

$$\mathcal{M}^{-}(D^{2}(\eta v - u)) < gr^{-\gamma},$$

with $g := u - \eta v + f < 0$ for any $r \ne 0$. By Lemma 2.1, we obtain $(\eta v - u)' \le 0$, a contradiction with the fact that $(\eta v - u)$ achieves its minimum at zero. Thus, we have obtained that $u \le v$ and, by exchanging the roles of $u$ and $v$, we conclude that $u \equiv v$. $\square$

We can now prove that the smooth eigenvalue $\bar{\lambda}'_{\gamma}$ is actually achieved on smooth eigenfunctions.

**Theorem 2.9.** There exists $u \in C(\overline{B(0, 1)}) \cap C^2(B(0, 1) \setminus \{0\})$, radial, strictly positive in $B(0, 1)$ and satisfying

$$\begin{cases} F(D^2u) + \bar{\lambda}'_{\gamma}ur^{-\gamma} = 0 & \text{in } B(0, 1) \setminus \{0\} \\ u = 0 & \text{on } \partial B(0, 1) \end{cases}$$

Furthermore, in $\overline{B(0, 1)}$, $u$ is Lipschitz continuous when $\gamma \le 1$ and Hölder continuous with exponent $2 - \gamma$ if $\gamma > 1$.

*Proof.* We consider a sequence $\{\lambda_n\}$, with $\lambda_n \to \bar{\lambda}'_{\gamma}$ and $\lambda_n < \bar{\lambda}'_{\gamma}$ and, for all $n$, the solution $u_n \in C(\overline{B(0, 1)}) \cap C^2(B(0, 1) \setminus \{0\})$ provided by Theorem 2.8 of

$$F(D^2u_n) + \lambda_n u_n r^{-\gamma} = -r^{-\gamma}, \quad u_n(1) = 0.$$
We claim that the positive sequence ${\|u_n\|_{\infty}}_n$ is unbounded. Indeed, arguing by contradiction, if ${u_n}_n$ is uniformly bounded, then, by using Theorem 2.6 and considering a subsequence if necessary, we obtain that there exists a solution $u \in C(\overline{B(0,1)}) \cap C^2(B(0,1) \setminus \{0\})$, $u \ge 0$, of

$$F(D^2u) + \bar{\lambda}'_\gamma ur^{-\gamma} = -r^{-\gamma}, \quad u(1) = 0.$$

Then, arguing as in the proof of Theorem 2.8, we deduce that $u$ is strictly positive in $B(0,1)$, and by taking $0 < \epsilon < \frac{1}{\|u\|_{\infty}}$, we see that $u$ satisfies

$$F(D^2u) + (\bar{\lambda}'_\gamma + \epsilon)ur^{-\gamma} \le 0,$$

a contradiction to the definition of $\bar{\lambda}'_\gamma$. It then follows that the sequence ${u_n}_n$ is not uniformly bounded. Normalizing and considering a subsequence, letting $n \to \infty$ yields the existence of $u \in C(\overline{B(0,1)}) \cap C^2(B(0,1) \setminus \{0\})$ satisfying $\|u\|_{\infty} = 1$ and

$$F(D^2u) + \bar{\lambda}'_\gamma ur^{-\gamma} = 0 \text{ in } B(0,1) \setminus \{0\}, \quad u(1) = 0.$$

Finally, the strict positivity of $u$ in $B(0,1)$ and its global regularity in $\overline{B(0,1)}$ follow by arguing as in the proof of Theorem 2.8. $\square$

The uniqueness, up to positive multiplicative constants, of the solution given by Theorem 2.9 is provided by the last result of the present section.

**Proposition 2.10.** The eigenvalue $\bar{\lambda}'_\gamma$ is simple.

*Proof.* From Theorem 2.9, there exists a bounded eigenfunction $v$. Let $u$ be another eigenfunction.

We define $\eta := \sup\{t > 0, tv < u\}$, or, equivalently, $\frac{1}{\eta} = \sup_{B(0,1)\setminus\{0\}} \frac{v}{u}$, which is well defined by Hopf Lemma applied to $u$ and $v$ on $\partial B(0,1)$. Then, the function $u - \eta v$ is nonnegative in $\overline{B(0,1)} \setminus \{0\}$ and it satisfies

$$M^-(D^2(u - \eta v)) \le 0 \quad \text{in } B(0,1) \setminus \{0\}.$$

If, by contradiction, $u - \eta v > 0$ at some point in $B(0,1) \setminus \{0\}$, then, by the strong maximum principle, $u - \eta v > 0$ in the whole of $B(0,1) \setminus \{0\}$.

We now distinguish the cases $u(0)$ finite or infinite.

- In the case $u(0) = +\infty$, one necessarily has

$$\frac{1}{\eta} = \lim_{r \to 1} \frac{v}{u} = \frac{v'(1)}{u'(1)}$$

but this contradicts Hopf's Lemma.
We claim that the positive sequence ${\|u_n\|_\infty}_n$ is unbounded. Indeed, arguing by contradiction, if ${u_n}_n$ is uniformly bounded, then, by using Theorem 2.6 and considering a subsequence if necessary, we obtain that there exists a solution $u \in C(\overline{B(0,1)}) \cap C^2(B(0,1) \setminus \{0\})$, $u \ge 0$, of

$$F(D^2u) + \bar{\lambda}'_\gamma ur^{-\gamma} = -r^{-\gamma}, \quad u(1) = 0.$$

Then, arguing as in the proof of Theorem 2.8, we deduce that $u$ is strictly positive in $B(0,1)$, and by taking $0 < \epsilon < \frac{1}{\|u\|_\infty}$, we see that $u$ satisfies

$$F(D^2u) + (\bar{\lambda}'_\gamma + \epsilon)ur^{-\gamma} \le 0,$$

a contradiction to the definition of $\bar{\lambda}'_\gamma$. It then follows that the sequence ${u_n}_n$ is not uniformly bounded. Normalizing and considering a subsequence, letting $n \to \infty$ yields the existence of $u \in C(\overline{B(0,1)}) \cap C^2(B(0,1) \setminus \{0\})$ satisfying $\|u\|_\infty = 1$ and

$$F(D^2u) + \bar{\lambda}'_\gamma ur^{-\gamma} = 0 \text{ in } B(0,1) \setminus \{0\}, \quad u(1) = 0.$$

Finally, the strict positivity of $u$ in $B(0,1)$ and its global regularity in $\overline{B(0,1)}$ follow by arguing as in the proof of Theorem 2.8. $\square$

The uniqueness, up to positive multiplicative constants, of the solution given by Theorem 2.9 is provided by the last result of the present section.

**Proposition 2.10.** The eigenvalue $\bar{\lambda}'_\gamma$ is simple.

*Proof.* From Theorem 2.9, there exists a bounded eigenfunction $v$. Let $u$ be another eigenfunction.

We define $\eta := \sup\{t > 0, tv < u\}$, or, equivalently, $\frac{1}{\eta} = \sup_{B(0,1)\setminus\{0\}} \frac{v}{u}$, which is well defined by Hopf Lemma applied to $u$ and $v$ on $\partial B(0,1)$. Then, the function $u - \eta v$ is nonnegative in $\overline{B(0,1)} \setminus \{0\}$ and it satisfies

$$M^{-}(D^{2}(u - \eta v)) \le 0 \quad \text{in } B(0,1) \setminus \{0\}.$$

If, by contradiction, $u - \eta v > 0$ at some point in $B(0,1) \setminus \{0\}$, then, by the strong maximum principle, $u - \eta v > 0$ in the whole of $B(0,1) \setminus \{0\}$.

We now distinguish the cases $u(0)$ finite or infinite.

- In the case $u(0) = +\infty$, one necessarily has

$$\frac{1}{\eta} = \lim_{r \to 1} \frac{v}{u} = \frac{v'(1)}{u'(1)}$$

but this contradicts Hopf's Lemma.
- If $u(0) < \infty$, we have either $\frac{1}{\eta} = \lim_{r \to 1} \frac{v}{u} = \frac{v'(1)}{u'(1)}$ or $\frac{1}{\eta} = \frac{v(0)}{u(0)}$. Here, the contradiction follows either from Hopf's Lemma or from Remark 2.2, which gives $(u - \eta v)'(r) \le 0$ for $r > 0$ sufficiently small, so that $u - \eta v$ cannot have a strict minimum point at zero.
Thus, we have obtained that all the eigenfunctions are bounded and multiple of each others. $\square$

## 2.2. The eigenvalue inherited from some equation on $\mathbb{R}^+$

We suppose in this section that $F$ is one of Pucci's operators. We consider only the case $F = \mathcal{M}^+$, the changes to bring for $F = \mathcal{M}^-$ being obvious.

We present below an alternative proof of the existence of radial eigenfunctions related to the eigenvalue $\bar{\lambda}_{\gamma}$. Here, the idea is to show the existence of global solutions defined in $(0, +\infty)$ of the ODE associated with radial solutions of the equation

$$ \mathcal{M}^{+}(D^{2}u) = -\frac{u}{r^{\gamma}} \quad \text{in } \mathbb{R}^{N} \setminus \{0\}. $$

We recall that $u \in C^2(\mathbb{R}^N \setminus \{0\})$ is a radial solution of the above equation if and only if $u = u(r)$ is a $C^2((0, +\infty))$ solution of the second order ODE

$$ u'' = M_{+} \left( -\frac{(N-1)}{r} K_{+}(u') - \frac{u}{r^{\gamma}} \right) \quad \text{in } (0, +\infty), \qquad (2.4) $$

where

$$ K_{+}(s) = \begin{cases} \Lambda s & \text{if } s \ge 0 \\ \lambda s & \text{if } s < 0 \end{cases}, \quad M_{+}(s) = \begin{cases} \frac{1}{\Lambda} s & \text{if } s \ge 0 \\ \frac{1}{\lambda} s & \text{if } s < 0 \end{cases}. $$

**Theorem 2.11.** There exists a global solution $u \in C^2((0, +\infty))$ of equation (2.4), which extends as a continuous function on $[0, +\infty)$ satisfying $u(0) = 1$. Moreover, there exists $\bar{r} > 0$ such that $u(\bar{r}) = 0$ and $u(r) > 0$ for $0 \le r < \bar{r}$.

*Proof.* We begin by proving the local existence of $u$ near zero. We distinguish the cases $\gamma < 1$ and $\gamma \ge 1$.

### 1st case : $\gamma < 1$

For fixed $r_0 > 0$ to be conveniently chosen, let us define the function set

$$ V_{r_0} = \{u \in C([0, r_0]) : |u(r) - 1| \le \frac{1}{2}, u(0) = 1\}. $$

For $u \in V_{r_0}$, let us define

$$ T(u)(r) := 1 - \int_{0}^{r} \frac{1}{\lambda s^{N-1}} \int_{0}^{s} u(t) t^{N-1-\gamma} dt ds. $$
- If $u(0) < \infty$, we have either $\frac{1}{\eta} = \lim_{r \to 1} \frac{v}{u} = \frac{v'(1)}{u'(1)}$ or $\frac{1}{\eta} = \frac{v(0)}{u(0)}$. Here, the contradiction follows either from Hopf's Lemma or from Remark 2.2, which gives $(u - \eta v)'(r) \le 0$ for $r > 0$ sufficiently small, so that $u - \eta v$ cannot have a strict minimum point at zero.
Thus, we have obtained that all the eigenfunctions are bounded and multiple of each others. $\square$

## 2.2. The eigenvalue inherited from some equation on $\mathbb{R}^+$

We suppose in this section that $F$ is one of Pucci's operators. We consider only the case $F = \mathcal{M}^+$, the changes to bring for $F = \mathcal{M}^-$ being obvious.

We present below an alternative proof of the existence of radial eigenfunctions related to the eigenvalue $\bar{\lambda}_{\gamma}$. Here, the idea is to show the existence of global solutions defined in $(0, +\infty)$ of the ODE associated with radial solutions of the equation

$$ \mathcal{M}^{+}(D^{2}u) = -\frac{u}{r^{\gamma}} \quad \text{in } \mathbb{R}^{N} \setminus \{0\}. $$

We recall that $u \in C^2(\mathbb{R}^N \setminus \{0\})$ is a radial solution of the above equation if and only if $u = u(r)$ is a $C^2((0, +\infty))$ solution of the second order ODE

$$ u'' = M_{+}\left(-\frac{(N-1)}{r}K_{+}(u') - \frac{u}{r^{\gamma}}\right) \quad \text{in } (0, +\infty), \qquad (2.4) $$

where

$$ K_{+}(s) = \begin{cases} \Lambda s & \text{if } s \ge 0 \\ \lambda s & \text{if } s < 0 \end{cases}, \quad M_{+}(s) = \begin{cases} \frac{1}{\Lambda} s & \text{if } s \ge 0 \\ \frac{1}{\lambda} s & \text{if } s < 0 \end{cases}. $$

**Theorem 2.11.** There exists a global solution $u \in C^2((0, +\infty))$ of equation (2.4), which extends as a continuous function on $[0, +\infty)$ satisfying $u(0) = 1$. Moreover, there exists $\bar{r} > 0$ such that $u(\bar{r}) = 0$ and $u(r) > 0$ for $0 \le r < \bar{r}$.

*Proof.* We begin by proving the local existence of $u$ near zero. We distinguish the cases $\gamma < 1$ and $\gamma \ge 1$.

### 1st case : $\gamma < 1$

For fixed $r_0 > 0$ to be conveniently chosen, let us define the function set

$$ V_{r_0} = \{u \in C([0, r_0]) : |u(r) - 1| \le \frac{1}{2}, u(0) = 1\}. $$

For $u \in V_{r_0}$, let us define

$$ T(u)(r) := 1 - \int_{0}^{r} \frac{1}{\lambda s^{N-1}} \int_{0}^{s} u(t) t^{N-1-\gamma} dt ds. $$
We fix $r_0$ such that $r_0^{2-\gamma} \le \frac{\lambda(N-\gamma)(2-\gamma)}{3}$. Then, it is easy to verify that $T$ maps $V_{r_0}$ into itself and it is a contraction mapping on it. Let us denote by $u \in V_{r_0}$ its fixed point. It then follows that $u$ satisfies

$$u'(r) = -\frac{1}{\lambda r^{N-1}} \int_{0}^{r} u(t) t^{N-1-\gamma} dt$$

as well as

$$u'' = -\frac{u}{\lambda r^{\gamma}} - (N-1) \frac{u'}{r}.$$

Thus, we clearly have $u' \le 0$ and, if we prove that $u'' \le 0$ as well, then $u$ is a solution of (2.4) in $(0, r_0)$. We observe that $u' \le 0$ implies $u \le 1$ and, consequently,

$$u'(r) \ge -\frac{1}{\lambda r^{N-1}} \int_{0}^{r} t^{N-1-\gamma} dt = -\frac{r^{1-\gamma}}{\lambda(N-\gamma)}.$$

Hence

$$u(r) \ge 1 - \frac{r^{2-\gamma}}{\lambda(2-\gamma)(N-\gamma)}, \qquad (2.5)$$

which in turn implies

$$\begin{aligned} u'(r) &\le -\frac{1}{\lambda r^{N-1}} \int_{0}^{r} \left(1 - \frac{t^{2-\gamma}}{\lambda(2-\gamma)(N-\gamma)}\right) t^{N-1-\gamma} dt \\ &= -\frac{r^{1-\gamma}}{\lambda(N-\gamma)} + \frac{r^{3-2\gamma}}{\lambda^2(2-\gamma)(N-\gamma)(N+2-2\gamma)} \\ &\le -\frac{r^{1-\gamma}}{2\lambda(N-\gamma)} \end{aligned}$$

by the choice of $r_0$. Thus, we have proved that

$$-\frac{r^{1-\gamma}}{\lambda(N-\gamma)} \le u'(r) \le -\frac{r^{1-\gamma}}{2\lambda(N-\gamma)} \qquad (2.6)$$

From estimates (2.5) and (2.6), we further deduce

$$\begin{aligned} u'' &= -\frac{u}{\lambda r^{\gamma}} - (N-1) \frac{u'}{r} \\ &\le -\left(1 - \frac{r^{2-\gamma}}{\lambda(2-\gamma)(N-\gamma)}\right) \frac{r^{-\gamma}}{\lambda} + (N-1) \frac{r^{-\gamma}}{\lambda(N-\gamma)} \\ &= -\frac{(1-\gamma)r^{-\gamma}}{\lambda(N-\gamma)} + \frac{r^{2-2\gamma}}{\lambda^2(2-\gamma)(N-\gamma)} \\ &\le -\frac{(1-\gamma)r^{-\gamma}}{2\lambda(N-\gamma)} \end{aligned}$$
We fix $r_0$ such that $r_0^{2-\gamma} \le \frac{\lambda(N-\gamma)(2-\gamma)}{3}$. Then, it is easy to verify that $T$ maps $V_{r_0}$ into itself and it is a contraction mapping on it. Let us denote by $u \in V_{r_0}$ its fixed point. It then follows that $u$ satisfies

$$u'(r) = -\frac{1}{\lambda r^{N-1}} \int_{0}^{r} u(t) t^{N-1-\gamma} dt$$

as well as

$$u'' = -\frac{u}{\lambda r^{\gamma}} - (N-1) \frac{u'}{r}.$$

Thus, we clearly have $u' \le 0$ and, if we prove that $u'' \le 0$ as well, then $u$ is a solution of (2.4) in $(0, r_0)$. We observe that $u' \le 0$ implies $u \le 1$ and, consequently,

$$u'(r) \geq -\frac{1}{\lambda r^{N-1}} \int_{0}^{r} t^{N-1-\gamma} dt = -\frac{r^{1-\gamma}}{\lambda(N-\gamma)}.$$

Hence

$$u(r) \geq 1 - \frac{r^{2-\gamma}}{\lambda(2-\gamma)(N-\gamma)}, \qquad (2.5)$$

which in turn implies

$$\begin{align*}
u'(r) &\le -\frac{1}{\lambda r^{N-1}} \int_0^r \left(1 - \frac{t^{2-\gamma}}{\lambda(2-\gamma)(N-\gamma)}\right) t^{N-1-\gamma} dt \\
&= -\frac{r^{1-\gamma}}{\lambda(N-\gamma)} + \frac{r^{3-2\gamma}}{\lambda^2(2-\gamma)(N-\gamma)(N+2-2\gamma)} \\
&\le -\frac{r^{1-\gamma}}{2\lambda(N-\gamma)}
\end{align*}$$

by the choice of $r_0$. Thus, we have proved that

$$-\frac{r^{1-\gamma}}{\lambda(N-\gamma)} \le u'(r) \le -\frac{r^{1-\gamma}}{2\lambda(N-\gamma)} \qquad (2.6)$$

From estimates (2.5) and (2.6), we further deduce

$$\begin{align*}
u'' &= -\frac{u}{\lambda r^{\gamma}} - (N-1) \frac{u'}{r} \\
&\le -\left(1 - \frac{r^{2-\gamma}}{\lambda(2-\gamma)(N-\gamma)}\right) \frac{r^{-\gamma}}{\lambda} + (N-1) \frac{r^{-\gamma}}{\lambda(N-\gamma)} \\
&= -\frac{(1-\gamma)r^{-\gamma}}{\lambda(N-\gamma)} + \frac{r^{2-2\gamma}}{\lambda^2(2-\gamma)(N-\gamma)} \\
&\le -\frac{(1-\gamma)r^{-\gamma}}{2\lambda(N-\gamma)}
\end{align*}$$
if we fix $r_0$ by setting $r_0^{2-\gamma} = \frac{\lambda(1-\gamma)(2-\gamma)}{3}$. By this choice for $r_0$, we obtain that $u$ satisfies equation (2.4) in $(0, r_0]$ and, moreover, that $u \in C^1([0, r_0])$.

**2nd case :** $\gamma \ge 1$

In this case, the solution $u$ is expected to be locally convex near zero. Thus, for $v \in V_{r_0}$, we define the map

$$T(v)(r) = 1 - \int_{0}^{r} \frac{1}{\Lambda s^{\tilde{N}_{+}-1}} \int_{0}^{s} u(t) t^{\tilde{N}_{+}-1-\gamma} dt \, ds$$

As in the previous case, it is easy to check that if $r_0^{2-\gamma} \le \frac{\Lambda(\tilde{N}_+-\gamma)(2-\gamma)}{3}$, then $T$ is a contraction mapping on $V_{r_0}$. Let $u \in V_{r_0}$ be its fixed point. We then have

$$u'(r) = -\frac{1}{\Lambda r^{\tilde{N}_{+}-1}} \int_{0}^{r} u(t) t^{\tilde{N}_{+}-1-\gamma} dt$$

and

$$u'' = -\frac{(\tilde{N}_{+}-1)}{r} u' - \frac{u}{\Lambda r^{\gamma}}.$$

Thus, in order to prove that $u$ satisfies (2.4) in $(0, r_0]$, it is enough to show that $u'' \ge 0$. By using a bootstrap argument analogous to the one used in the first case, we deduce

$$-\frac{r^{1-\gamma}}{\Lambda(\tilde{N}_{+}-\gamma)} \le u' \le -\frac{r^{1-\gamma}}{\Lambda(\tilde{N}_{+}-\gamma)} + \frac{r^{3-2\gamma}}{\Lambda^2(\tilde{N}_{+}-\gamma)(2-\gamma)(\tilde{N}_{+}+2-2\gamma)}$$

as well as

$$1 - \frac{r^{2-\gamma}}{\Lambda(\tilde{N}_{+}-\gamma)(2-\gamma)} \le u \le 1 - \frac{r^{2-\gamma}}{\Lambda(\tilde{N}_{+}-\gamma)(2-\gamma)} + \frac{r^{4-2\gamma}}{2\Lambda^2(\tilde{N}_{+}-\gamma)(2-\gamma)^2(\tilde{N}_{+}+2-2\gamma)}.$$

We then have

$$\begin{aligned} u''(r) &= -\frac{u}{\Lambda r^{\gamma}} - \frac{(\tilde{N}_{+}-1)}{r} u' \\ &\ge \frac{\gamma-1}{\Lambda(\tilde{N}_{+}-\gamma)} r^{-\gamma} + \frac{3-2\gamma}{\Lambda^2(\tilde{N}_{+}-\gamma)(2-\gamma)(\tilde{N}_{+}+2-2\gamma)} r^{2-2\gamma} \\ &\quad - \frac{r^{4-3\gamma}}{2\Lambda^3(\tilde{N}_{+}-\gamma)(2-\gamma)^2(\tilde{N}_{+}+2-2\gamma)} \end{aligned}$$

By choosing $r_0$ such that

$$r_0^{2-\gamma} \le \Lambda(3-2\gamma)(2-\gamma)$$

we obtain in any case that $u'' \ge 0$ in $(0, r_0]$.

Thus, there exists a local solution $u$ of equation (2.4) in $(0, r_0]$, which is positive and satisfies the initial condition $u(0) = 1$. By observing the Lipschitz continuity of the functions $K_+$ and
if we fix $r_0$ by setting $r_0^{2-\gamma} = \frac{\lambda(1-\gamma)(2-\gamma)}{3}$. By this choice for $r_0$, we obtain that $u$ satisfies equation (2.4) in $(0, r_0]$ and, moreover, that $u \in C^1([0, r_0])$.

**2nd case :** $\gamma \ge 1$

In this case, the solution $u$ is expected to be locally convex near zero. Thus, for $v \in V_{r_0}$, we define the map

$$T(v)(r) = 1 - \int_{0}^{r} \frac{1}{\Lambda s^{\tilde{N}_{+}-1}} \int_{0}^{s} u(t) t^{\tilde{N}_{+}-1-\gamma} dt \, ds$$

As in the previous case, it is easy to check that if $r_0^{2-\gamma} \le \frac{\Lambda(\tilde{N}_+-\gamma)(2-\gamma)}{3}$, then $T$ is a contraction mapping on $V_{r_0}$. Let $u \in V_{r_0}$ be its fixed point. We then have

$$u'(r) = -\frac{1}{\Lambda r^{\tilde{N}_{+}-1}} \int_{0}^{r} u(t) t^{\tilde{N}_{+}-1-\gamma} dt$$

and

$$u'' = -\frac{(\tilde{N}_{+}-1)}{r}u' - \frac{u}{\Lambda r^{\gamma}}.$$

Thus, in order to prove that $u$ satisfies (2.4) in $(0, r_0]$, it is enough to show that $u'' \ge 0$. By using a bootstrap argument analogous to the one used in the first case, we deduce

$$-\frac{r^{1-\gamma}}{\Lambda(\tilde{N}_{+}-\gamma)} \le u' \le -\frac{r^{1-\gamma}}{\Lambda(\tilde{N}_{+}-\gamma)} + \frac{r^{3-2\gamma}}{\Lambda^2(\tilde{N}_{+}-\gamma)(2-\gamma)(\tilde{N}_{+}+2-2\gamma)}$$

as well as

$$1 - \frac{r^{2-\gamma}}{\Lambda(\tilde{N}_{+}-\gamma)(2-\gamma)} \le u \le 1 - \frac{r^{2-\gamma}}{\Lambda(\tilde{N}_{+}-\gamma)(2-\gamma)} + \frac{r^{4-2\gamma}}{2\Lambda^2(\tilde{N}_{+}-\gamma)(2-\gamma)^2(\tilde{N}_{+}+2-2\gamma)}.$$

We then have

$$\begin{aligned} u''(r) &= -\frac{u}{\Lambda r^{\gamma}} - \frac{(\tilde{N}_{+}-1)}{r}u' \\ &\ge \frac{\gamma-1}{\Lambda(\tilde{N}_{+}-\gamma)}r^{-\gamma} + \frac{3-2\gamma}{\Lambda^2(\tilde{N}_{+}-\gamma)(2-\gamma)(\tilde{N}_{+}+2-2\gamma)}r^{2-2\gamma} \\ &\quad - \frac{r^{4-3\gamma}}{2\Lambda^3(\tilde{N}_{+}-\gamma)(2-\gamma)^2(\tilde{N}_{+}+2-2\gamma)} \end{aligned}$$

By choosing $r_0$ such that

$$r_0^{2-\gamma} \le \Lambda(3-2\gamma)(2-\gamma)$$

we obtain in any case that $u'' \ge 0$ in $(0, r_0]$.

Thus, there exists a local solution $u$ of equation (2.4) in $(0, r_0]$, which is positive and satisfies the initial condition $u(0) = 1$. By observing the Lipschitz continuity of the functions $K_+$ and
$M_+$ and using Cauchy-Lipschitz Theorem, we can extend the solution $u$ as a global solution in $(0, +\infty)$.

It remains to prove that there exists a first point $\bar{r}$ such that $u(\bar{r}) = 0$ and $u(r) > 0$ for $0 \le r < \bar{r}$.

We argue by contradiction, and we assume that $u(r) > 0$ for all $r \ge 0$. Since $u'$ is initially negative, if there exists a first point $r_1 > 0$ such that $u'(r_1) = 0$, then $u''(r_1) \ge 0$. On the other hand, from equation (2.4), we deduce $u''(r_1) < 0$ since $u(r_1) > 0$. Hence, one has $u'(r) < 0$ for all $r > 0$. From (2.4), it then follows that, independently of the sign of $u''$, one has

$$u'' + \frac{N-1}{r}u' \le -\frac{ur^{-\gamma}}{\Lambda}.$$

Inspired by [21] and [22], let us introduce the function

$$y(r) = \frac{u'(r)}{u(r)}r^{N-1},$$

which is, then, negative and it satisfies

$$y' \le -\frac{r^{N-1-\gamma}}{\Lambda} - \frac{y^2}{r^{N-1}}.$$

By integrating between some $r_1 > 0$ and $r$, we obtain

$$y(r) + k(r) \le -c_1 r^{N-\gamma},$$

for some $c_1 > 0$ and $k(r) = \int_{r_1}^{r} \frac{y^2(t)}{t^{N-1}} dt$. This yields in particular $y(r) \le -c_1 r^{N-\gamma}$ and therefore, for $r$ sufficiently large,

$$k(r) \ge c_2 r^{N+2(1-\gamma)}. \qquad (2.7)$$

On the other hand, we also have

$$k(r) \le -y(r),$$

that is

$$k(r) \le \sqrt{k'(r)r^{N-1}}$$

which yields, after integration on $(r, +\infty)$,

$$k(r) \le (N-2)r^{N-2}. \qquad (2.8)$$

Being $N+2(1-\gamma) > N-2$, estimates (2.7) and (2.8) give a contradiction, showing that the constructed solution $u$ cannot be globally positive in $[0, +\infty)$. $\square$

<figure><img src="image_3.png" /></figure>
$M_+$ and using Cauchy-Lipschitz Theorem, we can extend the solution $u$ as a global solution in $(0, +\infty)$.

It remains to prove that there exists a first point $\bar{r}$ such that $u(\bar{r}) = 0$ and $u(r) > 0$ for $0 \le r < \bar{r}$.

We argue by contradiction, and we assume that $u(r) > 0$ for all $r \ge 0$. Since $u'$ is initially negative, if there exists a first point $r_1 > 0$ such that $u'(r_1) = 0$, then $u''(r_1) \ge 0$. On the other hand, from equation (2.4), we deduce $u''(r_1) < 0$ since $u(r_1) > 0$. Hence, one has $u'(r) < 0$ for all $r > 0$. From (2.4), it then follows that, independently of the sign of $u''$, one has

$$u'' + \frac{N-1}{r}u' \le -\frac{ur^{-\gamma}}{\Lambda}.$$

Inspired by [21] and [22], let us introduce the function

$$y(r) = \frac{u'(r)}{u(r)}r^{N-1},$$

which is, then, negative and it satisfies

$$y' \le -\frac{r^{N-1-\gamma}}{\Lambda} - \frac{y^2}{r^{N-1}}.$$

By integrating between some $r_1 > 0$ and $r$, we obtain

$$y(r) + k(r) \le -c_1 r^{N-\gamma},$$

for some $c_1 > 0$ and $k(r) = \int_{r_1}^{r} \frac{y^2(t)}{t^{N-1}} dt$. This yields in particular $y(r) \le -c_1 r^{N-\gamma}$ and therefore, for $r$ sufficiently large,

$$k(r) \ge c_2 r^{N+2(1-\gamma)}. \qquad (2.7)$$

On the other hand, we also have

$$k(r) \le -y(r),$$

that is

$$k(r) \le \sqrt{k'(r)r^{N-1}}$$

which yields, after integration on $(r, +\infty)$,

$$k(r) \le (N-2)r^{N-2}. \qquad (2.8)$$

Being $N+2(1-\gamma) > N-2$, estimates (2.7) and (2.8) give a contradiction, showing that the constructed solution $u$ cannot be globally positive in $[0, +\infty)$. $\square$

<figure><img src="image_2.png" /></figure>
As an immediate consequence of the above result and Theorem 2.7, we deduce the following

**Corollary 2.12.** Let $\bar{r}$ be defined as in Theorem 2.11. Then $\bar{\lambda}'_\gamma(B(0,1) \setminus \{0\}) = \bar{r}^{\gamma-2}$.

**Remark 2.13.** Let us observe that, as in the case of equations with continuous coefficients, one could prove the existence of a numerable set of radial eigenvalues, by proving the oscillatory behavior of the solution $u$ constructed in Theorem 2.11, see [21] and [22].

### 2.3. The stability of the principal eigenvalue and related eigenfunctions

The results of the present section give, as a corollary, the proof of Theorem 1.1.

Let us start by proving the stability with respect to the $\epsilon$-regularization of the singular potential. We recall that $r_\epsilon = (r^2 + \epsilon^2)^{\frac{1}{2}}$ and $\bar{\lambda}_\gamma^\epsilon = \bar{\lambda}(F, \frac{1}{r_\epsilon^\gamma}, B(0,1))$.

**Theorem 2.14.** One has

$$ \bar{\lambda}'_\gamma = \lim_{\epsilon \to 0} \bar{\lambda}_\gamma^\epsilon. $$

Furthermore, if $\{u_\epsilon\}$ is the sequence of the eigenfunctions associated with the eigenvalue $\bar{\lambda}_\gamma^\epsilon$ and satisfying $u_\epsilon(0) = 1$, then, one can extract from $\{u_\epsilon\}$ a subsequence uniformly converging on $\overline{B(0,1)}$ to the eigenfunction associated with $\bar{\lambda}'_\gamma$ which takes the value 1 at zero.

*Proof.* Let $\{u_\epsilon\}$ be the sequence as in the statement. Then, each $u_\epsilon$ is a smooth positive function in $B(0,1)$, radially symmetric by [16], satisfying in particular

$$ F(D^2 u_\epsilon) + \bar{\lambda}_\gamma^\epsilon \frac{u_\epsilon}{r^\gamma} \geq 0 \quad \text{in } B(0,1) \setminus \{0\}, $$

so that, by Theorem 2.7, one has $\bar{\lambda}_\gamma^\epsilon \geq \bar{\lambda}'_\gamma$. Moreover, the sequence $\{\bar{\lambda}_\gamma^\epsilon\}$ is monotone increasing with respect to $\epsilon$. Thus, we deduce

$$ \mu := \lim_{\epsilon \to 0} \bar{\lambda}_\gamma^\epsilon \geq \bar{\lambda}'_\gamma. $$

On the other hand, by the monotonicity properties of radially symmetric solutions of elliptic equations, we know that $u'_\epsilon(r) \leq 0$ for $r \in [0, 1]$. Since

$$ \mathcal{M}^+(D^2 u_\epsilon) \geq F(D^2 u_\epsilon) $$

we deduce that, independently of the sign of $u''_\epsilon(r)$, one has

$$ u''_\epsilon + (\tilde{N}_+ - 1) \frac{u'_\epsilon}{r} \geq - \frac{\bar{\lambda}_\gamma^\epsilon}{\lambda} u_\epsilon r_\epsilon^{-\gamma}. $$

This implies

$$ (u'_\epsilon r^{\tilde{N}_+ - 1})' \geq - \frac{\bar{\lambda}_\gamma^\epsilon}{\lambda} u_\epsilon \frac{r^{\tilde{N}_+ - 1}}{r_\epsilon^\gamma} \geq - \frac{\bar{\lambda}_\gamma^\epsilon}{\lambda} r^{\tilde{N}_+ - 1 - \gamma}, $$
As an immediate consequence of the above result and Theorem 2.7, we deduce the following

**Corollary 2.12.** Let $\bar{r}$ be defined as in Theorem 2.11. Then $\bar{\lambda}'_\gamma(B(0,1) \setminus \{0\}) = \bar{r}^{\gamma-2}$.

**Remark 2.13.** Let us observe that, as in the case of equations with continuous coefficients, one could prove the existence of a numerable set of radial eigenvalues, by proving the oscillatory behavior of the solution $u$ constructed in Theorem 2.11, see [21] and [22].

### 2.3. The stability of the principal eigenvalue and related eigenfunctions

The results of the present section give, as a corollary, the proof of Theorem 1.1.

Let us start by proving the stability with respect to the $\epsilon$-regularization of the singular potential. We recall that $r_\epsilon = (r^2 + \epsilon^2)^{\frac{1}{2}}$ and $\bar{\lambda}_\gamma^\epsilon = \bar{\lambda}(F, \frac{1}{r_\epsilon^\gamma}, B(0,1))$.

**Theorem 2.14.** One has

$$ \bar{\lambda}'_\gamma = \lim_{\epsilon \to 0} \bar{\lambda}_\gamma^\epsilon. $$

Furthermore, if $\{u_\epsilon\}$ is the sequence of the eigenfunctions associated with the eigenvalue $\bar{\lambda}_\gamma^\epsilon$ and satisfying $u_\epsilon(0) = 1$, then, one can extract from $\{u_\epsilon\}$ a subsequence uniformly converging on $\overline{B(0,1)}$ to the eigenfunction associated with $\bar{\lambda}'_\gamma$ which takes the value 1 at zero.

*Proof.* Let $\{u_\epsilon\}$ be the sequence as in the statement. Then, each $u_\epsilon$ is a smooth positive function in $B(0,1)$, radially symmetric by [16], satisfying in particular

$$ F(D^2u_\epsilon) + \bar{\lambda}_\gamma^\epsilon \frac{u_\epsilon}{r^\gamma} \geq 0 \quad \text{in } B(0,1) \setminus \{0\}, $$

so that, by Theorem 2.7, one has $\bar{\lambda}_\gamma^\epsilon \geq \bar{\lambda}'_\gamma$. Moreover, the sequence $\{\bar{\lambda}_\gamma^\epsilon\}$ is monotone increasing with respect to $\epsilon$. Thus, we deduce

$$ \mu := \lim_{\epsilon \to 0} \bar{\lambda}_\gamma^\epsilon \geq \bar{\lambda}'_\gamma. $$

On the other hand, by the monotonicity properties of radially symmetric solutions of elliptic equations, we know that $u'_\epsilon(r) \leq 0$ for $r \in [0,1]$. Since

$$ \mathcal{M}^+(D^2u_\epsilon) \geq F(D^2u_\epsilon) $$

we deduce that, independently of the sign of $u''_\epsilon(r)$, one has

$$ u''_\epsilon + (\tilde{N}_+ - 1) \frac{u'_\epsilon}{r} \geq - \frac{\bar{\lambda}_\gamma^\epsilon}{\lambda} u_\epsilon r_\epsilon^{-\gamma}. $$

This implies

$$ (u'_\epsilon r^{\tilde{N}_+ - 1})' \geq - \frac{\bar{\lambda}_\gamma^\epsilon}{\lambda} u_\epsilon \frac{r^{\tilde{N}_+ - 1}}{r_\epsilon^\gamma} \geq - \frac{\bar{\lambda}_\gamma^\epsilon}{\lambda} r^{\tilde{N}_+ - 1 - \gamma}, $$
and therefore, by integrating,

$$0 \ge u'_{\epsilon}(r) \ge - \frac{\bar{\lambda}_{\gamma}^{\epsilon}}{\lambda(\tilde{N}_{+} - \gamma)} r^{1-\gamma}.$$

Hence, on $\overline{B(0, 1)}$, the functions $u_{\epsilon}$ are uniformly Lipschitz continuous if $\gamma \le 1$, and uniformly $(2 - \gamma)$-Hölder continuous if $\gamma > 1$. In both cases, up to a subsequence, $\{u_{\epsilon}\}$ is uniformly converging to a continuous radial function $u \in C(\overline{B(0, 1)})$ which satisfies $u(0) = 1$ and

$$F(D^2 u) + \mu \frac{u}{r^{\gamma}} = 0.$$

Hence, $u$ is $C^2(B(0, 1) \setminus \{0\})$ and, by the standard strong maximum principle, $u$ is strictly positive in $B(0, 1)$. This yields, by definition, $\mu \le \bar{\lambda}'_{\gamma}$. Hence, $\mu = \bar{\lambda}'_{\gamma}$ and the conclusion follows from Proposition 2.10. $\square$

As a consequence of the previous theorem, we finally obtain the following

**Corollary 2.15.** One has

$$\bar{\lambda}_{\gamma} = \lim_{\delta \to 0} \bar{\lambda}_{\gamma} \left( B(0, 1) \setminus \overline{B(0, \delta)} \right) = \bar{\lambda}'_{\gamma}.$$

*Proof.* We observe that the function $\delta \mapsto \bar{\lambda}_{\gamma}(B(0, 1) \setminus B(0, \delta))$ is monotone increasing. Moreover, by their own definition, we have that

$$\bar{\lambda}'_{\gamma} \le \bar{\lambda}_{\gamma} \le \bar{\lambda}_{\gamma} (B(0, 1) \setminus B(0, \delta)) \quad \text{for all } \delta \ge 0.$$

On the other hand, by Theorem 2.14, for any $\eta > 0$ there exists $\epsilon_0 > 0$ such that

$$\bar{\lambda}_{\gamma}^{\epsilon_0} \le \bar{\lambda}'_{\gamma} + \frac{\eta}{2}.$$

Furthermore, by using the continuity of the principal eigenvalue with respect to the domain for equations with regular coefficients, see e.g. [14], there exists $\delta_0 > 0$ such that

$$\bar{\lambda}_{\gamma}^{\epsilon_0}(B(0, 1)) \setminus B(0, \delta_0) \le \bar{\lambda}_{\gamma}^{\epsilon_0} + \frac{\eta}{2} \le \bar{\lambda}'_{\gamma} + \eta.$$

Now, since $\epsilon \mapsto \bar{\lambda}_{\gamma}^{\epsilon}(B(0, 1)) \setminus B(0, \delta_0)$ decreases when $\epsilon$ decreases to zero, one gets

$$\bar{\lambda}_{\gamma}(B(0, 1) \setminus B(0, \delta_0)) \le \bar{\lambda}_{\gamma}^{\epsilon_0}(B(0, 1)) \setminus B(0, \delta_0) \le \bar{\lambda}'_{\gamma} + \eta,$$

which gives the conclusion. $\square$
and therefore, by integrating,

$$0 \geq u'_{\epsilon}(r) \geq - \frac{\bar{\lambda}_{\gamma}^{\epsilon}}{\lambda(\tilde{N}_{+} - \gamma)} r^{1-\gamma}.$$

Hence, on $\overline{B(0, 1)}$, the functions $u_{\epsilon}$ are uniformly Lipschitz continuous if $\gamma \leq 1$, and uniformly $(2 - \gamma)$-Hölder continuous if $\gamma > 1$. In both cases, up to a subsequence, $\{u_{\epsilon}\}$ is uniformly converging to a continuous radial function $u \in C(\overline{B(0, 1)})$ which satisfies $u(0) = 1$ and

$$F(D^2 u) + \mu \frac{u}{r^{\gamma}} = 0.$$

Hence, $u$ is $C^2(B(0, 1) \setminus \{0\})$ and, by the standard strong maximum principle, $u$ is strictly positive in $B(0, 1)$. This yields, by definition, $\mu \leq \bar{\lambda}'_{\gamma}$. Hence, $\mu = \bar{\lambda}'_{\gamma}$ and the conclusion follows from Proposition 2.10. $\square$

As a consequence of the previous theorem, we finally obtain the following

**Corollary 2.15.** One has

$$\bar{\lambda}_{\gamma} = \lim_{\delta \to 0} \bar{\lambda}_{\gamma} \left( B(0, 1) \setminus \overline{B(0, \delta)} \right) = \bar{\lambda}'_{\gamma}.$$

*Proof.* We observe that the function $\delta \mapsto \bar{\lambda}_{\gamma}(B(0, 1) \setminus B(0, \delta))$ is monotone increasing. Moreover, by their own definition, we have that

$$\bar{\lambda}'_{\gamma} \leq \bar{\lambda}_{\gamma} \leq \bar{\lambda}_{\gamma} (B(0, 1) \setminus B(0, \delta)) \quad \text{for all } \delta \geq 0.$$

On the other hand, by Theorem 2.14, for any $\eta > 0$ there exists $\epsilon_0 > 0$ such that

$$\bar{\lambda}_{\gamma}^{\epsilon_0} \leq \bar{\lambda}'_{\gamma} + \frac{\eta}{2}.$$

Furthermore, by using the continuity of the principal eigenvalue with respect to the domain for equations with regular coefficients, see e.g. [14], there exists $\delta_0 > 0$ such that

$$\bar{\lambda}_{\gamma}^{\epsilon_0}(B(0, 1)) \setminus B(0, \delta_0)) \leq \bar{\lambda}_{\gamma}^{\epsilon_0} + \frac{\eta}{2} \leq \bar{\lambda}'_{\gamma} + \eta.$$

Now, since $\epsilon \mapsto \bar{\lambda}_{\gamma}^{\epsilon}(B(0, 1)) \setminus B(0, \delta_0))$ decreases when $\epsilon$ decreases to zero, one gets

$$\bar{\lambda}_{\gamma}(B(0, 1) \setminus B(0, \delta_0)) \leq \bar{\lambda}_{\gamma}^{\epsilon_0}(B(0, 1)) \setminus B(0, \delta_0)) \leq \bar{\lambda}'_{\gamma} + \eta,$$

which gives the conclusion. $\square$
### 3. The case $\gamma = 2$: Proof of Theorem 1.2.

The aim of this section is to give the proof of Theorem 1.2.

Let us start by recalling that, in the semilinear case, the eigenvalue related to Laplace operator with an inverse quadratic potential can be defined by a variational approach, i.e. by considering the minimum problem

$$ \bar{\lambda}_2(\Delta) := \inf_{\substack{u \in H_0^1(B(0,1)) \\ \int_{B(0,1)} \frac{u^2}{|x|^2} dx = 1}} \int_{B(0,1)} |\nabla u|^2 dx. $$

In this case, one has

$$ \bar{\lambda}_2(\Delta) = \left( \frac{N-2}{2} \right)^2. $$

Indeed, on the one hand, by Hardy inequality, every function $u \in H_0^1(B(0,1))$ satisfies

$$ \int_{B(0,1)} \frac{u^2}{|x|^2} dx \le \left( \frac{2}{N-2} \right)^2 \int_{B(0,1)} |\nabla u|^2 dx $$

and then

$$ \bar{\lambda}_2(\Delta) \ge \left( \frac{N-2}{2} \right)^2. $$

On the other hand, for every $\epsilon > 0$, the function $u_\epsilon(r) = r^{-\frac{N-2}{2}+\epsilon}(-\log r)$ belongs to $H_0^1(B(0,1))$ and satisfies

$$ \int_{B(0,1)} |\nabla u_\epsilon|^2 = \left[ \left( \frac{N-2}{2} \right)^2 + \epsilon^2 \right] \int_{B(0,1)} \frac{|u_\epsilon|^2}{|x|^2} dx, $$

so that

$$ \bar{\lambda}_2(\Delta) \le \left( \frac{N-2}{2} \right)^2 + \epsilon^2, \quad \forall \epsilon > 0. $$

As it is well known, the infimum defining $\bar{\lambda}_2(\Delta)$ is not achieved, that is the Dirichlet problem

$$ \begin{cases} -\Delta u = \bar{\lambda}_2(\Delta) \frac{u}{r^2} & \text{in } B(0,1) \\ u = 0 & \text{on } \partial B(0,1) \end{cases} $$

has not finite energy solutions $u \in H_0^1(B(0,1))$.

A kind of variational approach is possible also in the fully nonlinear framework for the case of Pucci's operators. From now on, we consider the operator $M^+$, being obvious the changes to be made for the operator $M^-$.
### 3. The case $\gamma = 2$: Proof of Theorem 1.2.

The aim of this section is to give the proof of Theorem 1.2.

Let us start by recalling that, in the semilinear case, the eigenvalue related to Laplace operator with an inverse quadratic potential can be defined by a variational approach, i.e. by considering the minimum problem

$$ \bar{\lambda}_2(\Delta) := \inf_{\substack{u \in H_0^1(B(0,1)) \\ \int_{B(0,1)} \frac{u^2}{|x|^2} dx = 1}} \int_{B(0,1)} |\nabla u|^2 dx. $$

In this case, one has

$$ \bar{\lambda}_2(\Delta) = \left( \frac{N-2}{2} \right)^2. $$

Indeed, on the one hand, by Hardy inequality, every function $u \in H_0^1(B(0,1))$ satisfies

$$ \int_{B(0,1)} \frac{u^2}{|x|^2} dx \le \left( \frac{2}{N-2} \right)^2 \int_{B(0,1)} |\nabla u|^2 dx $$

and then

$$ \bar{\lambda}_2(\Delta) \ge \left( \frac{N-2}{2} \right)^2. $$

On the other hand, for every $\epsilon > 0$, the function $u_\epsilon(r) = r^{-\frac{N-2}{2}+\epsilon}(-\log r)$ belongs to $H_0^1(B(0,1))$ and satisfies

$$ \int_{B(0,1)} |\nabla u_\epsilon|^2 = \left[ \left( \frac{N-2}{2} \right)^2 + \epsilon^2 \right] \int_{B(0,1)} \frac{|u_\epsilon|^2}{|x|^2} dx, $$

so that

$$ \bar{\lambda}_2(\Delta) \le \left( \frac{N-2}{2} \right)^2 + \epsilon^2, \quad \forall \epsilon > 0. $$

As it is well known, the infimum defining $\bar{\lambda}_2(\Delta)$ is not achieved, that is the Dirichlet problem

$$ \begin{cases} -\Delta u = \bar{\lambda}_2(\Delta) \frac{u}{r^2} & \text{in } B(0,1) \\ u = 0 & \text{on } \partial B(0,1) \end{cases} $$

has not finite energy solutions $u \in H_0^1(B(0,1))$.

A kind of variational approach is possible also in the fully nonlinear framework for the case of Pucci's operators. From now on, we consider the operator $M^+$, being obvious the changes to be made for the operator $M^-$.
Let us introduce the space of functions

$$ \mathcal{V} = \{ u \in C^2([0, 1]) : u'(0) = 0, \text{supp}(u) \text{ compact in } [0, 1) \}, $$

endowed with the norm

$$ \|u\| = \left( \int_{0}^{1} |u'|^2 r^{\tilde{N}_{+}-1} dr \right)^{1/2}, $$

and let us denote by $\mathcal{H}_0^1$ the closure of $\mathcal{V}$. Then, for all $\gamma \le 2$, we can consider the minimum problem

$$ \bar{\lambda}_{\gamma, var} := \inf_{\substack{u \in \mathcal{H}_0^1 \\ \int_0^1 u^2 r^{\tilde{N}_+-1-\gamma} dr = 1}} \int_0^1 |u'|^2 r^{\tilde{N}_+-1} dr. $$

In the next results, we will relate the two values $\bar{\lambda}_{\gamma}(\mathcal{M}^+)$ and $\bar{\lambda}_{\gamma, var}$, and we will study their asymptotic behavior as $\gamma \to 2$.

**Theorem 3.1.** One has

$$ \bar{\lambda}_{2, var} = \left( \frac{\tilde{N}_+ - 2}{2} \right)^2. $$

*Proof.* For any $\epsilon > 0$, let $u_{\epsilon}(r) = r^{-\frac{\tilde{N}_+-2}{2}} + \epsilon(-\log r)$. Then, it is easy to check that $u_{\epsilon} \in \mathcal{H}_0^1$ and a direct computation shows that

$$ \int_0^1 |u_{\epsilon}'|^2 r^{\tilde{N}_+-1} dr = \left[ \left( \frac{\tilde{N}_+-2}{2} \right)^2 + \epsilon^2 \right] \int_0^1 u_{\epsilon}^2 r^{\tilde{N}_+-3} dr, $$

hence

$$ \bar{\lambda}_{2, var} \le \left( \frac{\tilde{N}_+ - 2}{2} \right)^2 + \epsilon^2, \quad \forall \epsilon > 0. $$

On the other hand, we observe that the function $u = r^{-\frac{\tilde{N}_+-2}{2}}(-\log r)$ satisfies, for $r > 0$,

$$ u'' + (\tilde{N}_+ - 1) \frac{u'}{r} = - \left( \frac{\tilde{N}_+ - 2}{2} \right)^2 u r^{-2}. $$

Let us multiply the above equation by $\frac{v^2}{u}r^{\tilde{N}_+-1}$, where $v \in \mathcal{V}$ is arbitrarily fixed. Since $\tilde{N}_+ > 2$, we have that $\frac{u'r^{\tilde{N}_+-1}}{u}$ tends to zero as $r \to 0$. As a consequence, integrating by parts, we get

$$ - \left( \frac{\tilde{N}_+ - 2}{2} \right)^2 \int_0^1 v^2 r^{\tilde{N}_+-3} dr = \int_0^1 \left( \frac{u'}{u} v - v' \right)^2 r^{\tilde{N}_+-1} dr - \int_0^1 (v')^2 r^{\tilde{N}_+-1} dr, $$
Let us introduce the space of functions

$$ \mathcal{V} = \{ u \in C^2([0, 1]) : u'(0) = 0, \text{supp}(u) \text{ compact in } [0, 1) \}, $$

endowed with the norm

$$ \|u\| = \left( \int_{0}^{1} |u'|^2 r^{\tilde{N}_{+}-1} dr \right)^{1/2}, $$

and let us denote by $\mathcal{H}_0^1$ the closure of $\mathcal{V}$. Then, for all $\gamma \le 2$, we can consider the minimum problem

$$ \bar{\lambda}_{\gamma, var} := \inf_{\substack{u \in \mathcal{H}_0^1 \\ \int_0^1 u^2 r^{\tilde{N}_+ - 1 - \gamma} dr = 1}} \int_0^1 |u'|^2 r^{\tilde{N}_+ - 1} dr. $$

In the next results, we will relate the two values $\bar{\lambda}_{\gamma}(\mathcal{M}^+)$ and $\bar{\lambda}_{\gamma, var}$, and we will study their asymptotic behavior as $\gamma \to 2$.

**Theorem 3.1.** One has

$$ \bar{\lambda}_{2, var} = \left( \frac{\tilde{N}_+ - 2}{2} \right)^2. $$

*Proof.* For any $\epsilon > 0$, let $u_{\epsilon}(r) = r^{-\frac{\tilde{N}_+ - 2}{2} + \epsilon}(-\log r)$. Then, it is easy to check that $u_{\epsilon} \in \mathcal{H}_0^1$ and a direct computation shows that

$$ \int_0^1 |u_{\epsilon}'|^2 r^{\tilde{N}_+ - 1} dr = \left[ \left( \frac{\tilde{N}_+ - 2}{2} \right)^2 + \epsilon^2 \right] \int_0^1 u_{\epsilon}^2 r^{\tilde{N}_+ - 3} dr, $$

hence

$$ \bar{\lambda}_{2, var} \le \left( \frac{\tilde{N}_+ - 2}{2} \right)^2 + \epsilon^2, \quad \forall \epsilon > 0. $$

On the other hand, we observe that the function $u = r^{-\frac{\tilde{N}_+ - 2}{2}}(-\log r)$ satisfies, for $r > 0$,

$$ u'' + (\tilde{N}_+ - 1) \frac{u'}{r} = - \left( \frac{\tilde{N}_+ - 2}{2} \right)^2 u r^{-2}. $$

Let us multiply the above equation by $\frac{v^2}{u} r^{\tilde{N}_+ - 1}$, where $v \in \mathcal{V}$ is arbitrarily fixed. Since $\tilde{N}_+ > 2$, we have that $\frac{u' r^{\tilde{N}_+ - 1}}{u}$ tends to zero as $r \to 0$. As a consequence, integrating by parts, we get

$$ - \left( \frac{\tilde{N}_+ - 2}{2} \right)^2 \int_0^1 v^2 r^{\tilde{N}_+ - 3} dr = \int_0^1 \left( \frac{u'}{u} v - v' \right)^2 r^{\tilde{N}_+ - 1} dr - \int_0^1 (v')^2 r^{\tilde{N}_+ - 1} dr, $$
which yields, by the arbitrariness of $v \in \mathcal{V}$,

$$ \bar{\lambda}_{2, var} \geq \left( \frac{\tilde{N}_{+} - 2}{2} \right)^{2}. $$

□

In order to establish the relationship between $\bar{\lambda}_{\gamma, var}$ and $\bar{\lambda}_{\gamma}$ we need to investigate on the monotonicity and convexity properties of the functions $u$ realizing the infimum in the definition of $\bar{\lambda}_{\gamma, var}$.

**Proposition 3.2.** Let $1 < \gamma < 2$ and assume that $u_{\gamma} \in \mathcal{H}_{0}^{1}$, with $u_{\gamma} \geq 0$, realizes the infimum defining $\bar{\lambda}_{\gamma, var}$. Then, $u_{\gamma} \in C^{2}((0, 1])$ is bounded, $u'_{\gamma} \leq 0$ and $u''_{\gamma} \geq 0$ in $(0, 1)$.

*Proof.* Since $u_{\gamma}$ is a minimum, for any $v \in \mathcal{H}_{0}^{1}$ one has

$$ \int_{0}^{1} u'_{\gamma} v' r^{\tilde{N}_{+}-1} = \bar{\lambda}_{\gamma, var} \int_{0}^{1} u_{\gamma} v r^{\tilde{N}_{+}-1-\gamma}. $$

In particular, $u_{\gamma}$ satisfies in the distributional sense

$$ -(u'_{\gamma} r^{\tilde{N}_{+}-1})' = \bar{\lambda}_{\gamma, var} u_{\gamma} r^{\tilde{N}_{+}-1-\gamma} \quad (3.1) $$

By regularity theory, this implies that $u_{\gamma}$ belongs to $C^{2}((0, 1])$, it is strictly positive in $(0, 1)$ and it satisfies $u_{\gamma}(1) = 0$. Let us prove that $u_{\gamma}$ is bounded and that it can be extended as a continuous function on $[0, 1]$. Indeed, by multiplying equation (3.1) by a smooth function $v \in C^{2}([0, 1])$, having compact support in $[0, 1)$ and satisfying $v(0) \neq 0$, and integrating on $[\epsilon, 1]$ for $\epsilon > 0$, one has

$$ \bar{\lambda}_{\gamma, var} \int_{\epsilon}^{1} u_{\gamma}(r) r^{\tilde{N}_{+}-1-\gamma} v(r) dr = \int_{\epsilon}^{1} u'_{\gamma}(r) r^{\tilde{N}_{+}-1} v'(r) dr + u'_{\gamma}(\epsilon) \epsilon^{\tilde{N}_{+}-1} v(\epsilon). $$

Letting $\epsilon$ go to zero, we deduce $\lim_{\epsilon \to 0} u'_{\gamma}(\epsilon) \epsilon^{\tilde{N}_{+}-1} = 0$. It then follows, again from (3.1), that $u'_{\gamma}(r) \leq 0$ and that there exists some positive constant $c_{0}$ such that

$$ u'_{\gamma}(r) \geq -c_{0} r^{1-\tilde{N}_{+}} \quad \text{for } r \in (0, 1]. $$

This implies that

$$ u_{\gamma}(r) = - \int_{r}^{1} u'_{\gamma}(s) ds \leq c_{0} \int_{r}^{1} s^{1-\tilde{N}_{+}} ds \leq d_{0} r^{2-\tilde{N}_{+}} $$
which yields, by the arbitrariness of $v \in \mathcal{V}$,

$$ \bar{\lambda}_{2, var} \geq \left( \frac{\tilde{N}_{+} - 2}{2} \right)^{2}. $$

□

In order to establish the relationship between $\bar{\lambda}_{\gamma, var}$ and $\bar{\lambda}_{\gamma}$ we need to investigate on the monotonicity and convexity properties of the functions $u$ realizing the infimum in the definition of $\bar{\lambda}_{\gamma, var}$.

**Proposition 3.2.** Let $1 < \gamma < 2$ and assume that $u_{\gamma} \in \mathcal{H}_{0}^{1}$, with $u_{\gamma} \geq 0$, realizes the infimum defining $\bar{\lambda}_{\gamma, var}$. Then, $u_{\gamma} \in C^{2}((0, 1])$ is bounded, $u'_{\gamma} \leq 0$ and $u''_{\gamma} \geq 0$ in $(0, 1)$.

*Proof.* Since $u_{\gamma}$ is a minimum, for any $v \in \mathcal{H}_{0}^{1}$ one has

$$ \int_{0}^{1} u'_{\gamma} v' r^{\tilde{N}_{+}-1} = \bar{\lambda}_{\gamma, var} \int_{0}^{1} u_{\gamma} v r^{\tilde{N}_{+}-1-\gamma}. $$

In particular, $u_{\gamma}$ satisfies in the distributional sense

$$ -(u'_{\gamma} r^{\tilde{N}_{+}-1})' = \bar{\lambda}_{\gamma, var} u_{\gamma} r^{\tilde{N}_{+}-1-\gamma} \quad (3.1) $$

By regularity theory, this implies that $u_{\gamma}$ belongs to $C^{2}((0, 1])$, it is strictly positive in $(0, 1)$ and it satisfies $u_{\gamma}(1) = 0$. Let us prove that $u_{\gamma}$ is bounded and that it can be extended as a continuous function on $[0, 1]$. Indeed, by multiplying equation (3.1) by a smooth function $v \in C^{2}([0, 1])$, having compact support in $[0, 1)$ and satisfying $v(0) \neq 0$, and integrating on $[\epsilon, 1]$ for $\epsilon > 0$, one has

$$ \bar{\lambda}_{\gamma, var} \int_{\epsilon}^{1} u_{\gamma}(r) r^{\tilde{N}_{+}-1-\gamma} v(r) dr = \int_{\epsilon}^{1} u'_{\gamma}(r) r^{\tilde{N}_{+}-1} v'(r) dr + u'_{\gamma}(\epsilon) \epsilon^{\tilde{N}_{+}-1} v(\epsilon). $$

Letting $\epsilon$ go to zero, we deduce $\lim_{\epsilon \to 0} u'_{\gamma}(\epsilon) \epsilon^{\tilde{N}_{+}-1} = 0$. It then follows, again from (3.1), that $u'_{\gamma}(r) \leq 0$ and that there exists some positive constant $c_{0}$ such that

$$ u'_{\gamma}(r) \geq -c_{0} r^{1-\tilde{N}_{+}} \quad \text{for } r \in (0, 1]. $$

This implies that

$$ u_{\gamma}(r) = - \int_{r}^{1} u'_{\gamma}(s) ds \leq c_{0} \int_{r}^{1} s^{1-\tilde{N}_{+}} ds \leq d_{0} r^{2-\tilde{N}_{+}} $$
with $d_0 = \frac{c_0}{\tilde{N}_+-2} > 0$, which, in turn, yields

$$u'_{\gamma}(r)r^{\tilde{N}_+-1} = -\bar{\lambda}_{\gamma,\text{var}} \int_0^r u_{\gamma}(s)s^{\tilde{N}_+-1-\gamma}ds \geq -\bar{\lambda}_{\gamma,\text{var}} d_0 \int_0^r s^{1-\gamma}ds = -c_1 r^{2-\gamma}.$$

Thus, we have

$$u'_{\gamma}(r) \geq -c_1 r^{1-\tilde{N}_++2-\gamma}$$

and, then,

$$u_{\gamma}(r) \leq d_1 r^{2-\tilde{N}_++2-\gamma}.$$

Iterating the above inequalities, we obtain that for all integers $j \geq 0$ such that $2 - \tilde{N}_+ + j(2 - \gamma) < 0$, there exist positive constants $c_j$ and $d_j$ satisfying

$$u'_{\gamma}(r) \geq -c_j r^{1-\tilde{N}_++j(2-\gamma)}, \quad u_{\gamma}(r) \leq d_j r^{2-\tilde{N}_++j(2-\gamma)}. \qquad (3.2)$$

Now, if there exists $j \in \mathbb{N}$ such that $2 - \tilde{N}_+ + j(2 - \gamma) = 0$, i.e. if $\frac{\tilde{N}_+-2}{2-\gamma} \in \mathbb{N}$, then, by integrating the estimates obtained at the $(j-1)$-th step, we obtain

$$u'_{\gamma}(r) \geq -c_j r^{-1}, \quad u_{\gamma}(r) \leq d_j (-\ln r).$$

Integrating once more, we finally deduce

$$u'_{\gamma}(r) \geq -c_{j+1}(-\ln r)r^{1-\gamma} \implies u_{\gamma}(r) \leq d_{j+1} = \frac{c_{j+1}}{(2-\gamma)^2}.$$

On the other hand, if $\frac{\tilde{N}_+-2}{2-\gamma}$ is not integer, by integrating estimates (3.2) for $j = \lfloor \frac{\tilde{N}_+-2}{2-\gamma} \rfloor$, we obtain

$$u'_{\gamma}(r) \geq -c_{j+1}r^{1-\tilde{N}_++(j+1)(2-\gamma)} \implies u_{\gamma}(r) \leq d_{j+1} = \frac{c_{j+1}}{2 - \tilde{N}_+ + (j+1)(2-\gamma)}.$$

This shows that, in any case, $u_{\gamma}$ is bounded.

Let us finally prove that $u''_{\gamma} \geq 0$. We introduce the function

$$y_{\gamma}(r) := (\tilde{N}_{+} - 1)u'_{\gamma}(r) + \bar{\lambda}_{\gamma,\text{var}}r^{1-\gamma}u_{\gamma}(r),$$

which verifies $y_{\gamma} = -ru''_{\gamma}$. Hence, we need to prove that $y_{\gamma}(r) \leq 0$ for $r \in (0,1]$. An easy computation shows that

$$y'_{\gamma}(r) + (\tilde{N}_{+} - 1) \frac{y_{\gamma}(r)}{r} = \bar{\lambda}_{\gamma,\text{var}} (u'_{\gamma}(r)r^{1-\gamma} + (1-\gamma)u_{\gamma}(r)r^{-\gamma}) \leq 0,$$

so that

$$(y_{\gamma}(r)r^{\tilde{N}_+-1})' \leq 0.$$

Since $u_{\gamma}$ is bounded and $-c_0 r^{1-\gamma} \leq u'_{\gamma}(r) \leq 0$, we deduce that $y_{\gamma}(r)r^{\tilde{N}_+-1} \to 0$ as $r \to 0$. Hence, $y_{\gamma}(r) \leq 0$ for $r > 0$.

□
with $d_0 = \frac{c_0}{\tilde{N}_+-2} > 0$, which, in turn, yields

$$u'_{\gamma}(r)r^{\tilde{N}_+-1} = -\bar{\lambda}_{\gamma,\text{var}} \int_0^r u_{\gamma}(s)s^{\tilde{N}_+-1-\gamma}ds \geq -\bar{\lambda}_{\gamma,\text{var}} d_0 \int_0^r s^{1-\gamma}ds = -c_1 r^{2-\gamma}.$$

Thus, we have

$$u'_{\gamma}(r) \geq -c_1 r^{1-\tilde{N}_++2-\gamma}$$

and, then,

$$u_{\gamma}(r) \leq d_1 r^{2-\tilde{N}_++2-\gamma}.$$

Iterating the above inequalities, we obtain that for all integers $j \geq 0$ such that $2 - \tilde{N}_+ + j(2 - \gamma) < 0$, there exist positive constants $c_j$ and $d_j$ satisfying

$$u'_{\gamma}(r) \geq -c_j r^{1-\tilde{N}_++j(2-\gamma)}, \quad u_{\gamma}(r) \leq d_j r^{2-\tilde{N}_++j(2-\gamma)}. \qquad (3.2)$$

Now, if there exists $j \in \mathbb{N}$ such that $2 - \tilde{N}_+ + j(2 - \gamma) = 0$, i.e. if $\frac{\tilde{N}_+-2}{2-\gamma} \in \mathbb{N}$, then, by integrating the estimates obtained at the $(j-1)$-th step, we obtain

$$u'_{\gamma}(r) \geq -c_j r^{-1}, \quad u_{\gamma}(r) \leq d_j (-\ln r).$$

Integrating once more, we finally deduce

$$u'_{\gamma}(r) \geq -c_{j+1}(-\ln r)r^{1-\gamma} \implies u_{\gamma}(r) \leq d_{j+1} = \frac{c_{j+1}}{(2-\gamma)^2}.$$

On the other hand, if $\frac{\tilde{N}_+-2}{2-\gamma}$ is not integer, by integrating estimates (3.2) for $j = \lfloor \frac{\tilde{N}_+-2}{2-\gamma} \rfloor$, we obtain

$$u'_{\gamma}(r) \geq -c_{j+1}r^{1-\tilde{N}_++(j+1)(2-\gamma)} \implies u_{\gamma}(r) \leq d_{j+1} = \frac{c_{j+1}}{2 - \tilde{N}_+ + (j+1)(2-\gamma)}.$$

This shows that, in any case, $u_{\gamma}$ is bounded.

Let us finally prove that $u''_{\gamma} \geq 0$. We introduce the function

$$y_{\gamma}(r) := (\tilde{N}_{+} - 1)u'_{\gamma}(r) + \bar{\lambda}_{\gamma,\text{var}}r^{1-\gamma}u_{\gamma}(r),$$

which verifies $y_{\gamma} = -ru''_{\gamma}$. Hence, we need to prove that $y_{\gamma}(r) \leq 0$ for $r \in (0,1]$. An easy computation shows that

$$y'_{\gamma}(r) + (\tilde{N}_{+} - 1) \frac{y_{\gamma}(r)}{r} = \bar{\lambda}_{\gamma,\text{var}} (u'_{\gamma}(r)r^{1-\gamma} + (1-\gamma)u_{\gamma}(r)r^{-\gamma}) \leq 0,$$

so that

$$(y_{\gamma}(r)r^{\tilde{N}_+-1})' \leq 0.$$

Since $u_{\gamma}$ is bounded and $-c_0 r^{1-\gamma} \leq u'_{\gamma}(r) \leq 0$, we deduce that $y_{\gamma}(r)r^{\tilde{N}_+-1} \to 0$ as $r \to 0$. Hence, $y_{\gamma}(r) \leq 0$ for $r > 0$.

□
**Corollary 3.3.** Let $\gamma \in ]1, 2[$. Then

$$ \bar{\lambda}_{\gamma}(\mathcal{M}^{+}) = \Lambda \bar{\lambda}_{\gamma, var}. $$

*Proof.* It is not difficult to prove that the infimum defining $\bar{\lambda}_{\gamma, var}$ is achieved for $1 < \gamma < 2$. Thus, there exists $v_{\gamma} \in \mathcal{H}_{0}^{1}$, which can be assumed to be positive in $[0, 1)$, and which satisfies in $(0, 1)$

$$ \Lambda v_{\gamma}'' + \lambda \frac{N-1}{r} v_{\gamma}' = -\Lambda \bar{\lambda}_{\gamma, var} v_{\gamma} r^{-\gamma}. $$

By Proposition 3.2, we have that $v_{\gamma}$ is bounded, $v_{\gamma}' \le 0$ and $v_{\gamma}'' \ge 0$, so that $v_{\gamma}$ satisfies

$$ \mathcal{M}^{+}(D^{2}v_{\gamma}) + \Lambda \bar{\lambda}_{\gamma, var}v_{\gamma}r^{-\gamma} = 0 \quad \text{in } B(0, 1) \setminus \{0\}. $$

By definition, it then follows that $\bar{\lambda}_{\gamma}(\mathcal{M}^{+}) = \bar{\lambda}_{\gamma}'(\mathcal{M}^{+}) \ge \Lambda \bar{\lambda}_{\gamma, var}$. Furthermore, analyzing the boundary condition, we get, by regularity, that $v_{\gamma}(1) = 0$ in the classical sense. If, by contradiction, $\Lambda \bar{\lambda}_{\gamma, var} < \bar{\lambda}_{\gamma}'(\mathcal{M}^{+})$, then Theorem 2.7 would give $v_{\gamma} \le 0$ in $B(0, 1)$, a contradiction. $\square$

**Corollary 3.4.** One has

$$ \lim_{\gamma \to 2} \bar{\lambda}_{\gamma}(\mathcal{M}^{+}) = \Lambda \left( \frac{\tilde{N}_{+} - 2}{2} \right)^{2}. $$

*Proof.* By Theorem 3.1 and Corollary 3.3, it is sufficient to prove that $\bar{\lambda}_{\gamma, var} \to \bar{\lambda}_{2, var}$ as $\gamma \to 2$.

We first observe that, by their own definition, $\bar{\lambda}_{2, var} \le \bar{\lambda}_{\gamma, var}$.

On the other hand, for any $\epsilon > 0$ there exists $v \in \mathcal{H}_{0}^{1}$ such that

$$ \int_{0}^{1} |v'|^{2} r^{\tilde{N}_{+}-1} dr \le (\bar{\lambda}_{2, var} + \epsilon) \int_{0}^{1} |v|^{2} r^{\tilde{N}_{+}-3} dr. $$

Moreover, there exists $\gamma_{0}$ sufficiently close to 2 in order that, for $\gamma \ge \gamma_{0}$,

$$ \int_{0}^{1} |v|^{2} r^{\tilde{N}_{+}-1-\gamma} dr \ge (1-\epsilon) \int_{0}^{1} |v|^{2} r^{\tilde{N}_{+}-3} dr. $$

Thus, one has

$$ \int_{0}^{1} |v'|^{2} r^{\tilde{N}_{+}-1} dr \le (\bar{\lambda}_{2, var} + \epsilon)(1-\epsilon)^{-1} \int_{0}^{1} |v|^{2} r^{\tilde{N}_{+}-1-\gamma} dr $$

which yields

$$ \bar{\lambda}_{\gamma, var} \le (\bar{\lambda}_{2, var} + \epsilon)(1 - \epsilon)^{-1}. $$

$\square$
**Corollary 3.3.** Let $\gamma \in ]1, 2[$. Then

$$ \bar{\lambda}_{\gamma}(\mathcal{M}^{+}) = \Lambda \bar{\lambda}_{\gamma, var}. $$

*Proof.* It is not difficult to prove that the infimum defining $\bar{\lambda}_{\gamma, var}$ is achieved for $1 < \gamma < 2$. Thus, there exists $v_{\gamma} \in \mathcal{H}_{0}^{1}$, which can be assumed to be positive in $[0, 1)$, and which satisfies in $(0, 1)$

$$ \Lambda v_{\gamma}'' + \lambda \frac{N-1}{r} v_{\gamma}' = -\Lambda \bar{\lambda}_{\gamma, var} v_{\gamma} r^{-\gamma}. $$

By Proposition 3.2, we have that $v_{\gamma}$ is bounded, $v_{\gamma}' \le 0$ and $v_{\gamma}'' \ge 0$, so that $v_{\gamma}$ satisfies

$$ \mathcal{M}^{+}(D^{2}v_{\gamma}) + \Lambda \bar{\lambda}_{\gamma, var}v_{\gamma}r^{-\gamma} = 0 \quad \text{in } B(0,1) \setminus \{0\}. $$

By definition, it then follows that $\bar{\lambda}_{\gamma}(\mathcal{M}^{+}) = \bar{\lambda}_{\gamma}'(\mathcal{M}^{+}) \ge \Lambda \bar{\lambda}_{\gamma, var}$. Furthermore, analyzing the boundary condition, we get, by regularity, that $v_{\gamma}(1) = 0$ in the classical sense. If, by contradiction, $\Lambda \bar{\lambda}_{\gamma, var} < \bar{\lambda}_{\gamma}'(\mathcal{M}^{+})$, then Theorem 2.7 would give $v_{\gamma} \le 0$ in $B(0,1)$, a contradiction. $\square$

**Corollary 3.4.** One has

$$ \lim_{\gamma \to 2} \bar{\lambda}_{\gamma}(\mathcal{M}^{+}) = \Lambda \left( \frac{\tilde{N}_{+} - 2}{2} \right)^{2}. $$

*Proof.* By Theorem 3.1 and Corollary 3.3, it is sufficient to prove that $\bar{\lambda}_{\gamma, var} \to \bar{\lambda}_{2, var}$ as $\gamma \to 2$.

We first observe that, by their own definition, $\bar{\lambda}_{2, var} \le \bar{\lambda}_{\gamma, var}$.

On the other hand, for any $\epsilon > 0$ there exists $v \in \mathcal{H}_{0}^{1}$ such that

$$ \int_{0}^{1} |v'|^{2} r^{\tilde{N}_{+}-1} dr \le (\bar{\lambda}_{2, var} + \epsilon) \int_{0}^{1} |v|^{2} r^{\tilde{N}_{+}-3} dr. $$

Moreover, there exists $\gamma_{0}$ sufficiently close to 2 in order that, for $\gamma \ge \gamma_{0}$,

$$ \int_{0}^{1} |v|^{2} r^{\tilde{N}_{+}-1-\gamma} dr \ge (1-\epsilon) \int_{0}^{1} |v|^{2} r^{\tilde{N}_{+}-3} dr. $$

Thus, one has

$$ \int_{0}^{1} |v'|^{2} r^{\tilde{N}_{+}-1} dr \le (\bar{\lambda}_{2, var} + \epsilon)(1-\epsilon)^{-1} \int_{0}^{1} |v|^{2} r^{\tilde{N}_{+}-1-\gamma} dr $$

which yields

$$ \bar{\lambda}_{\gamma, var} \le (\bar{\lambda}_{2, var} + \epsilon)(1 - \epsilon)^{-1}. $$

$\square$
We are now ready to prove statement (i) of Theorem 1.2.

**Theorem 3.5.** *One has*

$$ \bar{\lambda}_2(\mathcal{M}^+) = \Lambda \left( \frac{\tilde{N}_+ - 2}{2} \right)^2 $$

and the function $u(r) = r^{-\frac{\tilde{N}_+-2}{2}}(-\ln r)$ is an explicit solution of

$$ \begin{cases} \mathcal{M}^+(D^2u) + \bar{\lambda}_\gamma \frac{u}{r^\gamma} = 0 & \text{in } B(0,1) \setminus \{0\} \\ u = 0 & \text{on } \partial B(0,1) \end{cases} $$

*Proof.* For any positive constants $c_1$ and $c_2$, let us consider the function

$$ u(r) = r^{-\frac{\tilde{N}_+-2}{2}} (c_1(-\ln r) + c_2). \qquad (3.3) $$

An easy computation, analogous to the one made in the proof of Theorem 3.1, leads to

$$ \mathcal{M}^+(D^2u) + \Lambda \left( \frac{\tilde{N}_+ - 2}{2} \right)^2 \frac{u}{r^2} = 0 \quad \text{in } B(0,1) \setminus \{0\}. $$

This gives, by definition, that

$$ \bar{\lambda}_2(\mathcal{M}^+) \geq \Lambda \left( \frac{\tilde{N}_+ - 2}{2} \right)^2. $$

On the other hand, by observing that $\bar{\lambda}_\gamma \geq \bar{\lambda}_2$ for all $\gamma \leq 2$ and by using Corollary 3.4, we also have

$$ \bar{\lambda}_2(\mathcal{M}^+) \leq \lim_{\gamma \to 2} \bar{\lambda}_\gamma(\mathcal{M}^+) = \Lambda \left( \frac{\tilde{N}_+ - 2}{2} \right)^2. $$

□

As a consequence of Corollary 3.4 and Theorem 3.5, we immediately deduce the first stability property of $\bar{\lambda}_2(\mathcal{M}^+)$ stated in Theorem 1.2-(ii). The other ones are given by the following result.

**Corollary 3.6.** *For the operators $F = \mathcal{M}^{\pm}$, one has*

$$ \lim_{\delta \to 0} \bar{\lambda}_2(B(0,1) \setminus \overline{B(0,\delta)}) = \bar{\lambda}_2 = \lim_{\epsilon \to 0} \bar{\lambda}_2^\epsilon. $$
We are now ready to prove statement (i) of Theorem 1.2.

**Theorem 3.5.** *One has*

$$ \bar{\lambda}_2(\mathcal{M}^+) = \Lambda \left( \frac{\tilde{N}_+ - 2}{2} \right)^2 $$

and the function $u(r) = r^{-\frac{\tilde{N}_+-2}{2}}(-\ln r)$ is an explicit solution of

$$ \begin{cases} \mathcal{M}^+(D^2u) + \bar{\lambda}_{\gamma} \frac{u}{r^{\gamma}} = 0 & \text{in } B(0,1) \setminus \{0\} \\ u = 0 & \text{on } \partial B(0,1) \end{cases} $$

*Proof.* For any positive constants $c_1$ and $c_2$, let us consider the function

$$ u(r) = r^{-\frac{\tilde{N}_+-2}{2}} (c_1(-\ln r) + c_2). \qquad (3.3) $$

An easy computation, analogous to the one made in the proof of Theorem 3.1, leads to

$$ \mathcal{M}^+(D^2u) + \Lambda \left( \frac{\tilde{N}_+ - 2}{2} \right)^2 \frac{u}{r^2} = 0 \quad \text{in } B(0,1) \setminus \{0\}. $$

This gives, by definition, that

$$ \bar{\lambda}_2(\mathcal{M}^+) \geq \Lambda \left( \frac{\tilde{N}_+ - 2}{2} \right)^2 . $$

On the other hand, by observing that $\bar{\lambda}_{\gamma} \geq \bar{\lambda}_2$ for all $\gamma \leq 2$ and by using Corollary 3.4, we also have

$$ \bar{\lambda}_2(\mathcal{M}^+) \leq \lim_{\gamma \to 2} \bar{\lambda}_{\gamma}(\mathcal{M}^+) = \Lambda \left( \frac{\tilde{N}_+ - 2}{2} \right)^2 . $$

□

As a consequence of Corollary 3.4 and Theorem 3.5, we immediately deduce the first stability property of $\bar{\lambda}_2(\mathcal{M}^+)$ stated in Theorem 1.2-(ii). The other ones are given by the following result.

**Corollary 3.6.** *For the operators $F = \mathcal{M}^{\pm}$, one has*

$$ \lim_{\delta \to 0} \bar{\lambda}_2(B(0,1) \setminus \overline{B(0,\delta)}) = \bar{\lambda}_2 = \lim_{\epsilon \to 0} \bar{\lambda}_2^{\epsilon}. $$
Proof. We observe that

$$ \bar{\lambda}_2(B(0,1) \setminus \overline{B(0,\delta)}) \geq \bar{\lambda}_2 $$

and that $\bar{\lambda}_2(B(0,1) \setminus \overline{B(0,\delta)})$ is decreasing with respect to $\delta$. Hence,

$$ \lim_{\delta \to 0} \bar{\lambda}_2(B(0,1) \setminus \overline{B(0,\delta)}) \geq \bar{\lambda}_2. $$

In order to prove the reverse inequality, we can use Theorem 3.5 jointly with Corollary 3.4, as well as Corollary 2.15. Indeed, for any $\eta > 0$ let $\gamma_0 < 2$ such that, for all $\gamma_0 \leq \gamma < 2$, one has

$$ \bar{\lambda}_2 \geq \bar{\lambda}_{\gamma} - \eta. $$

Moreover, let $\delta_0 = \delta(\gamma, \eta)$ be such that, for any $\delta < \delta_0$,

$$ \bar{\lambda}_{\gamma}(B(0,1) \setminus \overline{B(0,\delta)}) \leq \bar{\lambda}_{\gamma} + \eta. $$

It then follows

$$ \bar{\lambda}_2 \geq \bar{\lambda}_{\gamma}(B(0,1) \setminus \overline{B(0,\delta)}) - 2\eta \geq \bar{\lambda}_2(B(0,1) \setminus \overline{B(0,\delta)}) - 2\eta. $$

The assertion concerning $\lim \bar{\lambda}_2^\epsilon$ can be proved in the same way by using Theorem 2.14. $\square$

In order to complete the proof of Theorem 1.2, it is enough to observe that statement (iii) immediately follows from statement (i), the definition of $\bar{\lambda}_2(F)$ and the ellipticity inequalities (1.9).

### 4. The case $\gamma > 2$: Proof of Theorem 1.3

This section is completely devoted to the proof of Theorem 1.3. Let us assume, by contradiction, that for $\gamma > 2$ there exists $u \in C^2(B(0,1) \setminus \{0\})$ positive and radial, satisfying

$$ M^{-}(D^{2}u) \leq -\mu u r^{-\gamma} \quad \text{in } B(0,1) \setminus \{0\}, $$

for some $\mu > 0$.

Then, arguing as in the proof of Lemma 2.1, it follows that $u'(r)$ has constant sign in a right neighborhood of zero. If $u'(r) \geq 0$ for $r$ small, then, by the equation, $u''(r) \leq 0$ and then we would have

$$ (u'r^{\bar{N}_+-1})' \leq -\frac{\mu u r^{\bar{N}_+-1-\gamma}}{\Lambda}. $$

This implies that $u'r^{\bar{N}_+-1}$ has a nonnegative limit for $r \to 0$, and if this limit is strictly positive, we get that $u$ becomes large negative as $r \to 0$, a contradiction. Then the limit is zero, and then from the inequality above we get $u'(r)r^{\bar{N}_+-1} < 0$, a contradiction again.
Proof. We observe that

$$ \bar{\lambda}_2(B(0,1) \setminus \overline{B(0,\delta)}) \geq \bar{\lambda}_2 $$

and that $\bar{\lambda}_2(B(0,1) \setminus \overline{B(0,\delta)})$ is decreasing with respect to $\delta$. Hence,

$$ \lim_{\delta \to 0} \bar{\lambda}_2(B(0,1) \setminus \overline{B(0,\delta)}) \geq \bar{\lambda}_2. $$

In order to prove the reverse inequality, we can use Theorem 3.5 jointly with Corollary 3.4, as well as Corollary 2.15. Indeed, for any $\eta > 0$ let $\gamma_0 < 2$ such that, for all $\gamma_0 \le \gamma < 2$, one has

$$ \bar{\lambda}_2 \ge \bar{\lambda}_\gamma - \eta. $$

Moreover, let $\delta_0 = \delta(\gamma, \eta)$ be such that, for any $\delta < \delta_0$,

$$ \bar{\lambda}_\gamma(B(0,1) \setminus \overline{B(0,\delta)}) \le \bar{\lambda}_\gamma + \eta. $$

It then follows

$$ \bar{\lambda}_2 \ge \bar{\lambda}_\gamma(B(0,1) \setminus \overline{B(0,\delta)}) - 2\eta \ge \bar{\lambda}_2(B(0,1) \setminus \overline{B(0,\delta)}) - 2\eta. $$

The assertion concerning $\lim \bar{\lambda}_2^\epsilon$ can be proved in the same way by using Theorem 2.14. $\square$

In order to complete the proof of Theorem 1.2, it is enough to observe that statement (iii) immediately follows from statement (i), the definition of $\bar{\lambda}_2(F)$ and the ellipticity inequalities (1.9).

### 4. The case $\gamma > 2$: Proof of Theorem 1.3

This section is completely devoted to the proof of Theorem 1.3. Let us assume, by contradiction, that for $\gamma > 2$ there exists $u \in C^2(B(0,1) \setminus \{0\})$ positive and radial, satisfying

$$ M^-(D^2u) \le -\mu u r^{-\gamma} \quad \text{in } B(0,1) \setminus \{0\}, $$

for some $\mu > 0$.

Then, arguing as in the proof of Lemma 2.1, it follows that $u'(r)$ has constant sign in a right neighborhood of zero. If $u'(r) \ge 0$ for $r$ small, then, by the equation, $u''(r) \le 0$ and then we would have

$$ (u'r^{\bar{N}_+-1})' \le -\frac{\mu u r^{\bar{N}_+-1-\gamma}}{\Lambda}. $$

This implies that $u'r^{\bar{N}_+-1}$ has a nonnegative limit for $r \to 0$, and if this limit is strictly positive, we get that $u$ becomes large negative as $r \to 0$, a contradiction. Then the limit is zero, and then from the inequality above we get $u'(r)r^{\bar{N}_+-1} < 0$, a contradiction again.
Therefore, we have $u'(r) \le 0$ for $r > 0$ sufficiently small. Then, we observe that, whatever is the sign of $u''$, one has

$$u'' + (\tilde{N}_{-} - 1) \frac{u'}{r} \le -\frac{\mu}{\Lambda} u r^{-\gamma}.$$

Thus,

$$(u'r^{\tilde{N}_{-}-1})' \le -\frac{\mu}{\Lambda} u r^{\tilde{N}_{-}-1-\gamma} < 0, \qquad (4.1)$$

and then $u'(r)r^{\tilde{N}_{-}-1}$ has a non positive limit as $r \to 0$. If the limit is strictly negative, then $u'(r)r^{\tilde{N}_{-}-1} \le -c < 0$ in a neighborhood of zero. Then, we get

$$u(r) \ge c_1 r^{2-\tilde{N}_{-}}$$

for $r$ small and a positive $c_1$. Hence, by integrating (4.1) between $r$ and $s > r$ sufficiently small, we deduce

$$-u'(s)s^{\tilde{N}_{-}-1} + u'(r)r^{\tilde{N}_{-}-1} \ge c_2 \int_r^s t^{1-\gamma} dt = c_2 \frac{s^{2-\gamma} - r^{2-\gamma}}{2-\gamma}$$

and, since $\gamma > 2$, this yields $u'(r) > 0$ for $r$ small enough: a contradiction.

Thus, one has $\lim_{r \to 0} u'r^{\tilde{N}_{-}-1} = 0$ and, by (4.1), $u'(r) < 0$ for $r > 0$.

Next, by an inductive argument analogous to the one used in the proof of Proposition 3.2, we prove that for all integer $j \ge 0$ such that $\tilde{N}_{-} - 2 + j(2 - \gamma) > 0$ and for $r$ sufficiently small, one has, for some $c_j > 0$,

$$u(r) \ge c_j r^{j(2-\gamma)}. \qquad (4.2)$$

Indeed, (4.2) holds true for $j = 0$, since $u' < 0$ and $u$ is positive. Let us suppose that (4.2) is true for $j$ and that $\tilde{N}_{-} - 2 + (j + 1)(2 - \gamma) > 0$. Then, by (4.1),

$$-u'(r)r^{\tilde{N}_{-}-1} \ge \frac{\mu}{\Lambda} c_j \int_0^r s^{\tilde{N}_{-}-1+j(2-\gamma)-\gamma} ds = \frac{\mu}{\Lambda} \frac{c_j}{\tilde{N}_{-}-2+(j+1)(2-\gamma)} r^{\tilde{N}_{-}-2+(j+1)(2-\gamma)},$$

which, by integration, yields (4.2) for $j + 1$.

Now, let us assume that $\frac{\tilde{N}_{-}-2}{\gamma-2}$ is not integer. Then, using estimate (4.2) with $j = \lfloor \frac{\tilde{N}_{-}-2}{\gamma-2} \rfloor$ jointly with (4.1), we deduce for $r_0 > r > 0$

$$-u'(r_0)r_0^{\tilde{N}_{-}-1} + u'(r)r^{\tilde{N}_{-}-1} \ge \frac{\mu c_j}{\Lambda(\tilde{N}_{-}-2+(j+1)(2-\gamma))} s^{\tilde{N}_{-}-2+(j+1)(2-\gamma)} \Big|_{r}^{r_0}.$$

Since $\tilde{N}_{-} - 2 + (j + 1)(2 - \gamma) < 0$, this yields the contradiction

$$\lim_{r \to 0} u'(r) r^{\tilde{N}_{-}-1} = +\infty.$$
Therefore, we have $u'(r) \le 0$ for $r > 0$ sufficiently small. Then, we observe that, whatever is the sign of $u''$, one has

$$u'' + (\tilde{N}_- - 1) \frac{u'}{r} \le -\frac{\mu}{\Lambda} u r^{-\gamma}.$$

Thus,

$$(u'r^{\tilde{N}_--1})' \le -\frac{\mu}{\Lambda} u r^{\tilde{N}_--1-\gamma} < 0, \quad (4.1)$$

and then $u'(r)r^{\tilde{N}_--1}$ has a non positive limit as $r \to 0$. If the limit is strictly negative, then $u'(r)r^{\tilde{N}_--1} \le -c < 0$ in a neighborhood of zero. Then, we get

$$u(r) \ge c_1 r^{2-\tilde{N}_-}$$

for $r$ small and a positive $c_1$. Hence, by integrating (4.1) between $r$ and $s > r$ sufficiently small, we deduce

$$-u'(s)s^{\tilde{N}_--1} + u'(r)r^{\tilde{N}_--1} \ge c_2 \int_r^s t^{1-\gamma} dt = c_2 \frac{s^{2-\gamma} - r^{2-\gamma}}{2-\gamma}$$

and, since $\gamma > 2$, this yields $u'(r) > 0$ for $r$ small enough: a contradiction.

Thus, one has $\lim_{r \to 0} u'r^{\tilde{N}_--1} = 0$ and, by (4.1), $u'(r) < 0$ for $r > 0$.

Next, by an inductive argument analogous to the one used in the proof of Proposition 3.2, we prove that for all integer $j \ge 0$ such that $\tilde{N}_- - 2 + j(2 - \gamma) > 0$ and for $r$ sufficiently small, one has, for some $c_j > 0$,

$$u(r) \ge c_j r^{j(2-\gamma)}. \quad (4.2)$$

Indeed, (4.2) holds true for $j = 0$, since $u' < 0$ and $u$ is positive. Let us suppose that (4.2) is true for $j$ and that $\tilde{N}_- - 2 + (j+1)(2-\gamma) > 0$. Then, by (4.1),

$$-u'(r)r^{\tilde{N}_--1} \ge \frac{\mu}{\Lambda} c_j \int_0^r s^{\tilde{N}_--1+j(2-\gamma)-\gamma} ds = \frac{\mu c_j}{\Lambda \tilde{N}_- - 2 + (j+1)(2-\gamma)} r^{\tilde{N}_--2+(j+1)(2-\gamma)},$$

which, by integration, yields (4.2) for $j + 1$.

Now, let us assume that $\frac{\tilde{N}_--2}{\gamma-2}$ is not integer. Then, using estimate (4.2) with $j = \lfloor \frac{\tilde{N}_--2}{\gamma-2} \rfloor$ jointly with (4.1), we deduce for $r_0 > r > 0$

$$-u'(r_0)r_0^{\tilde{N}_--1} + u'(r)r^{\tilde{N}_--1} \ge \frac{\mu c_j}{\Lambda (\tilde{N}_- - 2 + (j+1)(2-\gamma))} s^{\tilde{N}_--2+(j+1)(2-\gamma)} \Big|_{r}^{r_0}.$$

Since $\tilde{N}_- - 2 + (j+1)(2-\gamma) < 0$, this yields the contradiction

$$\lim_{r \to 0} u'(r) r^{\tilde{N}_--1} = +\infty.$$
On the other hand, if $\frac{\tilde{N}_- - 2}{\gamma - 2} = j + 1$ is integer, then $\tilde{N}_- - 2 + j(2 - \gamma) = \gamma - 2 > 0$, and from (4.2) it follows that

$$u(r) \geq c_j r^{\gamma - \tilde{N}_-},$$

hence

$$u(r) r^{\tilde{N}_- - 1 - \gamma} \geq c_j r^{-1}.$$

From (4.1) we then deduce, for $0 < r < r_0$,

$$-u'(r_0) r_0^{\tilde{N}_- - 1} + u'(r) r^{\tilde{N}_- - 1} \geq \frac{\mu c_j}{\Lambda} (\ln r_0 - \ln r)$$

and we reach, also in this case, the contradiction

$$\lim_{r \to 0} u'(r) r^{\tilde{N}_- - 1} = +\infty.$$

## Acknowledgement

This work has been partially supported by INdAM-GNAMPA Research Project CUP_E53C22001930001.

## References

[1] J. P. García Azorero, I. Peral, Hardy inequalities and some critical elliptic and parabolic problems, J. Differential Equations 144 (2) (1998) 441–476. doi:10.1006/jdeq.1997.3375. URL https://doi.org/10.1006/jdeq.1997.3375

[2] I. Peral, F. Soria, Elliptic and Parabolic Equations Involving the Hardy-Leray Potential, Series Nonlinear Analysis and Applications, De Gruyter, 2021.

[3] G. Stampacchia, Le problème de Dirichlet pour les équations elliptiques du second ordre à coefficients discontinus, Ann. Inst. Fourier (Grenoble) 15 (1965) 189–258. URL http://www.numdam.org/item?id=AIF19651511890

[4] F. C. Cirstea, A complete classification of the isolated singularities for nonlinear elliptic equations with inverse square potentials, Mem. Amer. Math. Soc. 227 (1068) (2014) vi+85.

[5] F. C. Cirstea, Y. Du, Isolated singularities for weighted quasilinear elliptic equations, J. Funct. Anal. 259 (1) (2010) 174–202. doi:10.1016/j.jfa.2010.03.015. URL https://doi.org/10.1016/j.jfa.2010.03.015
On the other hand, if $\frac{\tilde{N}_- - 2}{\gamma - 2} = j + 1$ is integer, then $\tilde{N}_- - 2 + j(2 - \gamma) = \gamma - 2 > 0$, and from (4.2) it follows that

$$u(r) \geq c_j r^{\gamma - \tilde{N}_-},$$

hence

$$u(r) r^{\tilde{N}_- - 1 - \gamma} \geq c_j r^{-1}.$$

From (4.1) we then deduce, for $0 < r < r_0$,

$$-u'(r_0) r_0^{\tilde{N}_- - 1} + u'(r) r^{\tilde{N}_- - 1} \geq \frac{\mu c_j}{\Lambda} (\ln r_0 - \ln r)$$

and we reach, also in this case, the contradiction

$$\lim_{r \to 0} u'(r) r^{\tilde{N}_- - 1} = +\infty.$$

## Acknowledgement

This work has been partially supported by INdAM-GNAMPA Research Project CUP_E53C22001930001.

## References

[1] J. P. García Azorero, I. Peral, Hardy inequalities and some critical elliptic and parabolic problems, J. Differential Equations 144 (2) (1998) 441–476. doi:10.1006/jdeq.1997.3375. URL https://doi.org/10.1006/jdeq.1997.3375

[2] I. Peral, F. Soria, Elliptic and Parabolic Equations Involving the Hardy-Leray Potential, Series Nonlinear Analysis and Applications, De Gruyter, 2021.

[3] G. Stampacchia, Le problème de Dirichlet pour les équations elliptiques du second ordre à coefficients discontinus, Ann. Inst. Fourier (Grenoble) 15 (1965) 189–258. URL http://www.numdam.org/item?id=AIF19651511890

[4] F. C. Cirstea, A complete classification of the isolated singularities for nonlinear elliptic equations with inverse square potentials, Mem. Amer. Math. Soc. 227 (1068) (2014) vi+85.

[5] F. C. Cirstea, Y. Du, Isolated singularities for weighted quasilinear elliptic equations, J. Funct. Anal. 259 (1) (2010) 174–202. doi:10.1016/j.jfa.2010.03.015. URL https://doi.org/10.1016/j.jfa.2010.03.015
[6] F. C. Cirstea, Y. Du, Asymptotic behavior of solutions of semilinear elliptic equations near an isolated singularity, J. Funct. Anal. 250 (2) (2007) 317-346.
doi:10.1016/j.jfa.2007.05.005.
URL https://doi.org/10.1016/j.jfa.2007.05.005

[7] N. Chaudhuri, F. C. Cirstea, On trichotomy of positive singular solutions associated with the Hardy-Sobolev operator, C. R. Math. Acad. Sci. Paris 347 (3-4) (2009) 153-158.
doi:10.1016/j.crma.2008.12.018.
URL https://doi.org/10.1016/j.crma.2008.12.018

[8] H. Berestycki, L. Nirenberg, S. R. S. Varadhan, The principal eigenvalue and maximum principle for second-order elliptic operators in general domains, Comm. Pure Appl. Math. 47 (1) (1994) 47-92.
doi:10.1002/cpa.3160470105.
URL https://doi.org/10.1002/cpa.3160470105

[9] P.-L. Lions, Bifurcation and optimal stochastic control, Nonlinear Anal. 7 (2) (1983) 177-207.
doi:10.1016/0362-546X(83)90081-0.
URL https://doi.org/10.1016/0362-546X(83)90081-0

[10] H. Ishii, Y. Yoshimura, Demi-eigenvalues for uniformly elliptic isaacs operators, unpublished Manuscript (2006).

[11] I. Birindelli, F. Demengel, Eigenvalue, maximum principle and regularity for fully non linear homogeneous operators, Commun. Pure Appl. Anal. 6 (2) (2007) 335-366.
doi:10.3934/cpaa.2007.6.335.
URL https://doi.org/10.3934/cpaa.2007.6.335

[12] I. Birindelli, F. Demengel, Regularity and uniqueness of the first eigenfunction for singular fully nonlinear operators, J. Differential Equations 249 (5) (2010) 1089-1110.
doi:10.1016/j.jde.2010.03.015.
URL https://doi.org/10.1016/j.jde.2010.03.015

[13] I. Birindelli, F. Leoni, F. Pacella, Symmetry and spectral properties for viscosity solutions of fully nonlinear equations, J. Math. Pures Appl. (9) 107 (4) (2017) 409-428.
doi:10.1016/j.matpur.2016.07.004.
URL https://doi.org/10.1016/j.matpur.2016.07.004

[14] A. Quaas, B. Sirakov, Principal eigenvalues and the Dirichlet problem for fully nonlinear elliptic operators, Adv. Math. 218 (1) (2008) 105-135.
doi:10.1016/j.aim.2007.12.002.
URL https://doi.org/10.1016/j.aim.2007.12.002

[15] L. A. Caffarelli, X. Cabré, Fully nonlinear elliptic equations, Vol. 43 of American Mathematical Society Colloquium Publications, American Mathematical Society, Providence,
[6] F. C. Cirstea, Y. Du, Asymptotic behavior of solutions of semilinear elliptic equations near an isolated singularity, J. Funct. Anal. 250 (2) (2007) 317-346.
doi:10.1016/j.jfa.2007.05.005.
URL https://doi.org/10.1016/j.jfa.2007.05.005

[7] N. Chaudhuri, F. C. Cirstea, On trichotomy of positive singular solutions associated with the Hardy-Sobolev operator, C. R. Math. Acad. Sci. Paris 347 (3-4) (2009) 153-158.
doi:10.1016/j.crma.2008.12.018.
URL https://doi.org/10.1016/j.crma.2008.12.018

[8] H. Berestycki, L. Nirenberg, S. R. S. Varadhan, The principal eigenvalue and maximum principle for second-order elliptic operators in general domains, Comm. Pure Appl. Math. 47 (1) (1994) 47-92.
doi:10.1002/cpa.3160470105.
URL https://doi.org/10.1002/cpa.3160470105

[9] P.-L. Lions, Bifurcation and optimal stochastic control, Nonlinear Anal. 7 (2) (1983) 177-207.
doi:10.1016/0362-546X(83)90081-0.
URL https://doi.org/10.1016/0362-546X(83)90081-0

[10] H. Ishii, Y. Yoshimura, Demi-eigenvalues for uniformly elliptic Isaacs operators, unpublished Manuscript (2006).

[11] I. Birindelli, F. Demengel, Eigenvalue, maximum principle and regularity for fully non linear homogeneous operators, Commun. Pure Appl. Anal. 6 (2) (2007) 335-366.
doi:10.3934/cpaa.2007.6.335.
URL https://doi.org/10.3934/cpaa.2007.6.335

[12] I. Birindelli, F. Demengel, Regularity and uniqueness of the first eigenfunction for singular fully nonlinear operators, J. Differential Equations 249 (5) (2010) 1089-1110.
doi:10.1016/j.jde.2010.03.015.
URL https://doi.org/10.1016/j.jde.2010.03.015

[13] I. Birindelli, F. Leoni, F. Pacella, Symmetry and spectral properties for viscosity solutions of fully nonlinear equations, J. Math. Pures Appl. (9) 107 (4) (2017) 409-428.
doi:10.1016/j.matpur.2016.07.004.
URL https://doi.org/10.1016/j.matpur.2016.07.004

[14] A. Quaas, B. Sirakov, Principal eigenvalues and the Dirichlet problem for fully nonlinear elliptic operators, Adv. Math. 218 (1) (2008) 105-135.
doi:10.1016/j.aim.2007.12.002.
URL https://doi.org/10.1016/j.aim.2007.12.002

[15] L. A. Caffarelli, X. Cabré, Fully nonlinear elliptic equations, Vol. 43 of American Mathematical Society Colloquium Publications, American Mathematical Society, Providence,
RI, 1995. doi:10.1090/coll/043.

URL https://doi.org/10.1090/coll/043

[16] F. Da Lio, B. Sirakov, Symmetry results for viscosity solutions of fully nonlinear uniformly elliptic equations, J. Eur. Math. Soc. (JEMS) 9 (2) (2007) 317–330.
doi:10.4171/JEMS/81.

URL https://doi.org/10.4171/JEMS/81

[17] L. Caffarelli, Y. Y. Li, L. Nirenberg, Some remarks on singular solutions of nonlinear elliptic equations. I, J. Fixed Point Theory Appl. 5 (2) (2009) 353–395.
doi:10.1007/s11784-009-0107-8.

URL https://doi.org/10.1007/s11784-009-0107-8

[18] L. Boccardo, L. Orsina, Sublinear elliptic equations with singular potentials, Adv. Nonlinear Stud. 12 (2) (2012) 187–198.
doi:10.1515/ans-2012-0201.

URL https://doi.org/10.1515/ans-2012-0201

[19] L. Boccardo, L. Orsina, I. Peral, A remark on existence and optimal summability of solutions of elliptic problems involving Hardy potential, Discrete Contin. Dyn. Syst. 16 (3) (2006) 513–523.
doi:10.3934/dcds.2006.16.513.

URL https://doi.org/10.3934/dcds.2006.16.513

[20] M. G. Crandall, H. Ishii, P.-L. Lions, User's guide to viscosity solutions of second order partial differential equations, Bull. Amer. Math. Soc. (N.S.) 27 (1) (1992) 1–67.
doi:10.1090/S0273-0979-1992-00266-5.

URL https://doi.org/10.1090/S0273-0979-1992-00266-5

[21] J. Busca, M. J. Esteban, A. Quaas, Nonlinear eigenvalues and bifurcation problems for Pucci's operators, Ann. Inst. H. Poincaré C Anal. Non Linéaire 22 (2) (2005) 187–206.
doi:10.1016/j.anihpc.2004.05.004.

URL https://doi.org/10.1016/j.anihpc.2004.05.004

[22] F. Demengel, Generalized eigenvalues for fully nonlinear singular or degenerate operators in the radial case, Adv. Differential Equations 14 (11-12) (2009) 1127–1154.
RI, 1995. doi:10.1090/coll/043.

URL https://doi.org/10.1090/coll/043

[16] F. Da Lio, B. Sirakov, Symmetry results for viscosity solutions of fully nonlinear uniformly elliptic equations, J. Eur. Math. Soc. (JEMS) 9 (2) (2007) 317–330.
doi:10.4171/JEMS/81.

URL https://doi.org/10.4171/JEMS/81

[17] L. Caffarelli, Y. Y. Li, L. Nirenberg, Some remarks on singular solutions of nonlinear elliptic equations. I, J. Fixed Point Theory Appl. 5 (2) (2009) 353–395.
doi:10.1007/s11784-009-0107-8.

URL https://doi.org/10.1007/s11784-009-0107-8

[18] L. Boccardo, L. Orsina, Sublinear elliptic equations with singular potentials, Adv. Nonlinear Stud. 12 (2) (2012) 187–198.
doi:10.1515/ans-2012-0201.

URL https://doi.org/10.1515/ans-2012-0201

[19] L. Boccardo, L. Orsina, I. Peral, A remark on existence and optimal summability of solutions of elliptic problems involving Hardy potential, Discrete Contin. Dyn. Syst. 16 (3) (2006) 513–523.
doi:10.3934/dcds.2006.16.513.

URL https://doi.org/10.3934/dcds.2006.16.513

[20] M. G. Crandall, H. Ishii, P.-L. Lions, User's guide to viscosity solutions of second order partial differential equations, Bull. Amer. Math. Soc. (N.S.) 27 (1) (1992) 1–67.
doi:10.1090/S0273-0979-1992-00266-5.

URL https://doi.org/10.1090/S0273-0979-1992-00266-5

[21] J. Busca, M. J. Esteban, A. Quaas, Nonlinear eigenvalues and bifurcation problems for Pucci's operators, Ann. Inst. H. Poincaré C Anal. Non Linéaire 22 (2) (2005) 187–206.
doi:10.1016/j.anihpc.2004.05.004.

URL https://doi.org/10.1016/j.anihpc.2004.05.004

[22] F. Demengel, Generalized eigenvalues for fully nonlinear singular or degenerate operators in the radial case, Adv. Differential Equations 14 (11-12) (2009) 1127–1154.