Estimation of sources and factors affecting indoor VOC levels using basic numerical methods

1.   Introduction and statement of the main results

In his pioneering paper ^[35], Serrin refined Alexandrov's method of moving planes ^[6] and employed it to prove that a function defined in a bounded domain having constant nonzero Laplacian and satisfying constant Dirichlet and Neumann boundary condition is necessarily rotationally symmetric. In fact, in ^[35] it was shown that similar results hold more generally for a vast class of elliptic PDEs. Serrin's result and the method of moving planes have given rise to a prolific field of research on overdetermined boundary value problems. Nowadays Serrin's technique has been employed to study analogous problems in a number of works, among which we mention ^[8,14,23,27]. An alternative proof of Serrin's result was given by Weinberger ^[38] using a completely different method, essentially based on a comparison argument with the model solutions. Weinberger's method and in particular its use of a P-function has also experienced quite some success ^{[16,18,19,22,34]}. For further insights on these methods and their applications, see also ^[28] and references therein.

The above mentioned papers all gave rise to interesting characterizations of solutions supported in a ball. In this work, we are instead interested in characterizations of solutions supported in an annulus. Let us explain our setting in more detail. Let $E_i\Subset E_o\subset \mathbb{R}^2$ be two nonempty simply connected bounded domains, whose boundaries $\partial E_o = \Gamma_o$ and $\partial E_i = \Gamma_i$ are simple closed smooth curves, and let $\Omega = E_o\setminus \overline{E_i}$. By construction, $\Omega$ is a bounded domain that is not simply connected, with boundary

We are interested in studying pairs $(\Omega, u)$, where $\Omega$ has the above structure and $u:\Omega\to \mathbb{R}$ is a solution to the following problem

where $a, b, \alpha, \beta$ are constants and $\nu$ is the unit-normal pointing outside $\Omega$. This overdetermined boundary value problem has been well studied, especially in the case where $a > b$ and $\alpha\geq 0$. Under this assumption, the method of moving planes shows its strength and allows to prove some powerful rigidity results. A particularly nice one is the following, due to Reichel.

Theorem 1.1. ([31,Theorem 2]). Let $(\Omega, u)$ be a solution to problem $(1.1)$ such that $a > b$ and $b < u < a$ in $\Omega$. Then $\Omega$ is an annulus and $u$ is rotationally symmetric with $\partial u/ \partial |x| < 0$.

It is worth mentioning that Reichel's result is slightly stronger, as it also works for more general PDEs and it allows $E_i$ (and thus $\Gamma_i$) to be disconnected. In ^[36], Sirakov generalized previous results ^[1,5,40] and proved that Reichel's thesis remains true if one allows different values of $a$ and $\alpha$ on different connected components of $\Gamma_i$. He was also able to replace the hypothesis $u < a$ with the weaker $\alpha\geq 0$. We also remark that similar symmetry results have been proven to hold for more general families of elliptic operators. In ^[36] it is shown that the same result works for quasi-linear regular strongly elliptic operator. Similar results also hold for quasi-linear possibly degenerate elliptic operators ^[5], fully nonlinear operators ^[32] and for the fractional Laplacian ^[37]. For further symmetry results on this and related problems, see also ^{[7,16,17,29,33]}.

On the other hand, it seems that less attention has been brought to the case $a < b$. One possible explaination is that the moving plane method works well when the model solution are monotonically decreasing along the radial coordinate, but seems to be harder to employ when the model solutions are monotonically increasing. In this work, our aim is to discuss this case and to prove the analogue of Theorem 1.1, namely:

Theorem 1.2. Let $(\Omega, u)$ be a solution to problem $(1.1)$ such that $a < b$ and $a < u < b$ in $\Omega$. Then $\Omega$ is an annulus and $u$ is rotationally symmetric with $\partial u/ \partial |x| > 0$.

The strategy of the proof is based on the method developed in ^[2]. There, problem (1.1) has been studied in the case $a = b = 0$, using a new comparison strategy. Of course, under that hypothesis, monotonicity is lost and there are some complications, especially concerning the behaviour of the solution near the set of the maxima of the function $u$. However, we will show that the core of the method employed in ^[2] is perfectly suited to study the case discussed in this paper. In fact, we can adapt most of the argument in ^[2] to prove Theorem 1.2.

Perhaps surprisingly, our method proves to be less effective in the more studied case $b < a$. We will try to explain along the work why this is the case. Nevertheless, we will be able to employ our method to prove a slightly weaker version of Theorem 1.1 (namely, we will need to assume the additional condition $2a+\alpha^2\leq 2b+\beta^2$, see Theorem 4.1). As we will see, the proof is also less elementary than the one for Theorem 1.2. It does seem then that the moving plane method is still the one best suited to characterize rotationally symmetric solutions that are radially decreasing. Nevertheless, our alternative proof may be of interest: it is novel and it expands on our method in ways that may find further applications in the future.

Let us conclude with some final comments on the above results and on future applications. First of all, one natural question would be whether similar characterizations can be achieved for model solutions that are not monotonic along the radial coordinate. In this case, the situation is less clear as there are negative results: in ^[25] it is proven that there exist solutions to (1.1) with $\alpha = \beta$ that are not rotationally symmetric. In the case $a = b = 0$, there are both positive ([2,Theorem B and Theorem C]) and negative results ([2,Theorem A]). Another natural question concerns whether our method can be applied to more general elliptic equations (in the same way as the Weinberger's comparison method ^[38] was then developed in ^[18,19,22]) or to other problems in Riemannian geometry. In fact, the most part of our arguments (with the notable exception of the Pohozaev identity) do not rely on the structure of the Euclidean space at all. As a couple of examples of further applications, this method has proven to be quite successful to characterize static spacetimes in General Relativity ^{[10,11,12,13]} and somewhat similar techniques have been employed for other problems in General Relativity ^[4,21] and for $p$-harmonic functions in manifolds with nonnegative Ricci curvature in ^[3,20].

The paper is structured as follows. In Section 2 we show that for any acceptable choice of $a, b, \alpha, \beta$ there is a rotationally symmetric solution solving problem (1.1). This is necessary in order to start our method, as we need to select a model solution to compare with. Section 3 is the heart of the paper: a crucial gradient estimate (Theorem 3.1) is introduced and exploited to obtain some important length bounds for the boundary components (Proposition 3.2). In Subsection 3.3 these two results are then used, in combination with a divergence theorem argument, to prove Theorem 1.2. Unfortunately, the method described in Subsection 3.3 does not work to prove Theorem 1.1. In Section 4 we then develop a new alternative argument, based on a combination of the Pohozaev identity and the isoperimetric inequality, leading to the proof of Theorem 4.1, which is a weaker version of Theorem 1.1. Finally, in Section 5 we further comment on the complications that one has to deal with when studying the case $a > b$, trying to analyze where they come from and how one may try to overcome them. Ultimately, the purpose of this final section is that of understanding what are the limits of our method and how far they can be pushed.

2.   Setting up the comparison argument

Our first aim is that of understanding whether there are some relations that are always in place between the constants $a, \alpha, b, \beta$. The main tool that will help us in this regard is the well known Pohozaev identity ^[30]. In our framework, this formula has the following form

where $X = (x_1, x_2)$ is the position vector and $\nu$ is the exterior normal to $\partial\Omega$. We use $d\ell$, $dA$ to denote the line element and area element respectively, whereas $\langle\, \cdot\, |\, \cdot\, \rangle$ is the standard scalar product between vectors of $\mathbb{R}^2$. Recalling that $\operatorname{div}X = 2$ and using the divergence theorem, we deduce

From this formula we deduce some relations between $a, b, \alpha, \beta$, summarized in the following result.

Proposition 2.1. Let $(\Omega, u)$ be a solution to problem $(1.1)$.

● If $a > b$ and $b < u < a$ in $\Omega$, then $4a+\alpha^2 > 4b+\beta^2$, $\alpha\geq 0$, $\beta < 0$,

● If $a < b$ and $a < u < b$ in $\Omega$, then $4a+\alpha^2 > 4b+\beta^2$, $\alpha < 0$, $\beta\geq 0$.

Proof. We start by recalling that $\Omega = E_o\setminus\overline{E_i}$, so that in particular $|E_o| = |\Omega|+|E_i|$. As a consequence, we can rewrite (2.1) in the following two forms

If $b < u < a$ in $\Omega$, then the left hand side of the first equation is negative, whereas if $a < u < b$ in $\Omega$, then the left hand side of the second equation is negative. In both cases, we deduce $4a+\alpha^2 > 4b+\beta^2$.

The signs for $\alpha$ and $\beta$ follow immediately from the hypothesis that $a$ is the maximum (resp. minimum) of $u$ and $b$ is the minimum (resp. maximum) of $u$. The fact that $\beta\neq 0$ when $b < u < a$ and that $\alpha\neq 0$ when $a < u < b$ is a consequence of the Hopf Lemma (recall $\Delta u = -2 < 0$).

In order to start our comparison argument, we first need to take a close look at the model solutions. In particular, we first have to make sure that we always have a model solution to compare with, for any acceptable value of $a, b, \alpha, \beta$.

It can be easily checked that the rotationally symmetric solutions to $\Delta u = -2$ have the form

with $L, M\in \mathbb{R}$. When $M\leq 0$, these solutions are monotonically decreasing with respect to $|x|$. When $M > 0$, these solutions are monotonically increasing for $|x| < \sqrt{ M}$ and they are monotonically decreasing for $|x| > \sqrt{ M}$.

Remark 1. As in ^[2], one can fix the value of $L$ by means of a rescaling of the function and of the domain, ending up with the one-parameter family

However, for our purposes in this paper this is not necessary, as it is actually easier to work directly with $(2.2)$.

We expect that for any acceptable choice (according to Proposition 2.1) of $a, b, \alpha, \beta$, there should be a solution of the form (2.2) solving our problem (1.1). This is true in the case $a < b$, as guaranteed by the following lemma:

Lemma 2.2. If $a < b$, $\alpha < 0$, $\beta\geq 0$, $4a+\alpha^2 > 4b+\beta^2$, then there exist constants $L$ and $M > 0$ and two radii $0 < r_i < r_o\leq\sqrt{ M}$ so that the function $u\, = \, L-|x|^2/2+ M\log |x|$ satisfies problem $(1.1)$ in the annulus $\Omega = \{r_i < |x| < r_o\}$.

Proof. The proof of the lemma is more of a long exercise. We start by imposing the conditions we want on $L, M, r_o, r_i$, namely

and we employ the last two equations to obtain

Since we are assuming $\alpha < 0$ and $\beta\geq 0$, in order for $r_o, r_i$ to be positive, we need to choose the plus sign in (2.4) and we need $M > 0$.

Substituting in the first two equations, we then deduce that we have a solution if and only if we can find a zero for the function

Taking the limit at zero and infinity, we compute

It follows by continuity of $F$ that there is a value $0 < M < \infty$ such that $F(M) = 0$, as wished.

The case $a > b$ appears to be more complicated. In the next lemma we will show that we can still find a model to compare with as long as we assume the additional condition $2a+\alpha^2\leq 2b+\beta^2$. This is the first instance in which we notice complications in the case $a > b$. More specifically, it seems that negative values of $M$ are harder to deal with. In this respect, the hypothesis $2a+\alpha^2\leq 2b+\beta^2$ is helpful as it forces $M\geq 0$.

Lemma 2.3. Let $a > b$, $\alpha\geq 0$, $\beta < 0$, $4a+\alpha^2 > 4b+\beta^2$. If we further have $2a+\alpha^2\leq 2b+\beta^2$, then there exist constants $L$ and $M\geq 0$ and two radii $\sqrt{M}\leq r_i < r_o$ so that the function $u\, = \, L-|x|^2/2+ M\log |x|$ satisfies problem $(1.1)$ in the annulus $\Omega = \{r_i < |x| < r_o\}$.

Proof. Proceeding as in the proof of Lemma 2.2, we find again the expressions (2.4) for $r_o$ and $r_i$. Notice that, in contrast with the previous case, both signs for $\pm$ in (2.4) are possible. However, since we are ultimately interested in having $M\geq 0$, this again essentially forces us to choose the plus sign in the formulæ (2.4) for $r_i$ and $r_o$. Therefore, as in the previous lemma, we reduce ourselves to finding a zero for the function (2.5). We compute

Since we are assuming $4a+2\alpha^2-4b-2\beta^2\leq 0$, we conclude again that there is a value $0\leq M < \infty$ such that $F(M) = 0$.

3.   Exploiting the comparison argument

In the previous section we have developed a way to compare a general solution of (1.1) with a model solution (2.2). Building on that, in this section we develop our comparison technique. In Subsection 3.1 we prove a crucial gradient estimate, that is then exploited in Subsection 3.2 to prove length bounds for the boundary components $\Gamma_i$ and $\Gamma_o$. The results in both these subsection work in both the case $a < b$ and $b < a$. In Subsection 3.3 we will then specialize to the case $a < b$ and we will prove Theorem 1.2.

3.1.   Gradient estimate

Let $u$ be a solution to problem (1.1) for some fixed values of $a, \alpha, b, \beta$. In the case $a > b$, assume further that $2a+\alpha^2\leq 2b+\beta^2$. Then, Lemmata 2.2 and 2.3 tell us that there exist constants $L\in \mathbb{R}$, $M\geq 0$, $r_o > 0$, $r_i > 0$ such that the model solution

in the annulus $\{r_i < |x| < r_o\}$, solves problem (1.1) for the same values of $a, \alpha, b, \beta$.

From now on, the constants $L\in \mathbb{R}$, $M\geq 0$, $r_o > 0$, $r_i > 0$ will always be the ones prescribed by Lemmata 2.2 and 2.3. We are now ready to set up our comparison argument, in the spirit of ^[2]. We start by defining the pseudo-radial function $\Psi:\Omega\to \mathbb{R}$ implicitly via the following identity

It is easily seen that the function $G(\Psi) = L-\Psi^2/2+ M\log\Psi$ is monotone in $(0, \sqrt{ M}]$ and $[\sqrt{ M}, +\infty)$, hence the pseudo-radial function $\Psi = G^{-1}(u)$ is well defined as long as $\Psi^2\neq M$ inside $\Omega$. From now on we will always assume that $a < u < b$ or $b < u < a$, and from this it easily follows that the function $\Psi$ actually takes values between $r_i$ (attained at $\Gamma_i$) and $r_o$ (attained at $\Gamma_o$). On the other hand, Lemmata 2.2 and 2.3 tell us that either $r_i < r_o\leq \sqrt{ M}$ or $\sqrt{ M}\leq r_i < r_o$. Therefore, for all our purposes, the function $\Psi:\Omega\to[r_i, r_o]$ will always be well defined.

We also introduce the two functions

The function $W_0$ corresponds to the value that $| \nabla u|^2$ would have on the model solution. In fact, it is easily checked that, for the model solutions (2.2), it holds $W_0\equiv W$ pointwise. Our strategy is then that of comparing $W$ and $W_0$ in $\Omega$, ultimately trying to prove that they have to coincide.

First of all, notice that, by construction, $W$ and $W_0$ coincide on the boundary of $\Omega$, namely when $\Psi = r_o$ and when $\Psi = r_i$. We will now show that, under our hypotheses, $W$ is actually controlled by $W_0$ on the whole $\Omega$.

Theorem 3.1 (Gradient Estimates). Let $(\Omega, u)$ be a solution to problem $(1.1)$ with respect to some choice of $a, \alpha, b, \beta$. Suppose that $a < u < b$ or $b < u < a$. In the case $b < u < a$, assume further that $2a+\alpha^2\leq 2b+\beta^2$. Let $W$, $W_0$ be defined as above. Then it holds

Moreover, if $W = W_0$ at some point in the interior of $\Omega$, then $(\Omega, u)$ is rotationally symmetric.

Remark 2. This theorem is clearly inspired by [2,Theorem 3.1], however we mention that similar gradient estimates have already found applications in a number of problems. One of the most notable ones is the paper ^[9] and the subsequent developments in ^[15,26], where a gradient comparison argument very much resembling Theorem 3.1 has been obtained and exploited for static spacetimes in General Relativity. The introduction of the pseudo-radial function is instead more recent (it was exploited in the series of papers ^[11,12,13] for static spacetimes). This function is really helpful as it allows to have an explicit formula for $W_0$. This is crucial in the computations that follows: we will see that $\Psi$ will appear in all the important computations in this section.

Proof. The proof is essentially the same as the one in ^[2], but let us give some comments. Let $L\in \mathbb{R}$, $M\geq 0$, $r_o > 0$, $r_i > 0$ be the constants prescribed by Lemmata 2.2 and 2.3.

The crucial estimate is deduced starting from the following quantity

where $du\otimes du$ is the $2\times 2$ matrix with entries $\frac{ \partial u}{ \partial x^i}\frac{ \partial u}{ \partial x^j}$ and $g_{ \mathbb{R}^2}$ is the $2\times 2$ identity matrix. The heuristic behind the choice of the quantity above is essentially the fact that it can be computed to be zero on the model solution. Computing explicitly its square norm and using the fact that it is nonnegative, we obtain the following estimate for the hessian of $u$:

Starting now from the equality $\Delta| \nabla u|^2 = 2| \nabla^2 u|^2$ and plugging in this estimate, after some computations we obtain an elliptic inequality for the quantity $W-W_0$

Considering the function $F_\gamma = \gamma\, (W-W_0)$, where $\gamma = \gamma(\Psi) > 0$, one gets

where we have used $'$ to denote the differentiation with respect to $\Psi$. We now need to find a function $\gamma$ such that the coefficients of the zero order terms have the right sign. A good choice is to set

which corresponds to choosing

We know that $\Psi^2\neq M$ inside $\Omega$, hence $\gamma$ is well defined in the interior of $\Omega$. However, $\Psi^2$ may attain the value $M$ on the boundary $\partial\Omega$: this happens when $\alpha = 0$ (which implies from (2.4) that $r_i = \sqrt{ M}$) or when $\beta = 0$ (which implies from (2.4) that $r_o = \sqrt{ M}$). Let us for the moment assume that this does not happen, so that $F_\gamma$ is well defined on the whole $\Omega$. With this choice of $\gamma$, it is easily seen that $F_\gamma$ satisfies

Since $M\geq 0$, it follows that $F_\gamma$ satisfies the Maximum Principle in $\Omega$. Furthermore, $W = W_0$ on $\partial\Omega$ by construction, hence the Maximum Principle implies that $F_\gamma\leq 0$ (equivalently, $W\leq W_0$) on the whole $\Omega$. Furthermore, if the equality $W = W_0$ holds at one point $p$ in the interior of $N$, then, applying the Strong Maximum Principle in a neighborhood of $p$, we deduce that $W = W_0$ on the whole $\Omega$. It is then easy to conclude the desired rigidity statement with standard arguments, see for instance [10,Theorem 4.2].

It remains to discuss the case where $\Psi^2 = M$ on one of the boundary components, say $\Gamma_i$ (for $\Gamma_o$, the same argument applies). In this case, we consider $\varepsilon > 0$ and a small neighborhood $U_ \varepsilon = \{x\in\Omega\, :\, d(x, \Gamma_i) < \varepsilon\}$ of $\Gamma_i$, and we work on the domain $\Omega_ \varepsilon = \Omega-U_ \varepsilon$. Clearly $F_\gamma$ is well defined and satisfies the Maximum Principle in $\Omega_ \varepsilon$. We will then obtain the desired result by taking $\varepsilon\to 0$, provided we can prove that $F_\gamma\to 0$ as we approach $\Gamma_i$. It is easily seen that $F_\gamma$ goes to zero if and only if $W/\sqrt{W_0}$ goes to zero, which in turn is equivalent to $| \nabla u|^2/\sqrt{a-u}\to 0$, see [2,Lemma A.1]. The latter is granted by the Reverse Łojasiewicz Inequality [11,Theorem 2.2] (to be more precise, [11,Theorem 2.2] cannot be applied directly, as it is written for interior points only; however the proof extends without modifications to the case at hand).

3.2.   Length bounds

Following again ^[2], we can now exploit the gradient estimate proven above to obtain length bounds for both $\Gamma_i$ and $\Gamma_o$.

Proposition 3.2. Let $(\Omega, u)$ be a solution to problem $(1.1)$ with respect to some choice of $a, \alpha, b, \beta$. Suppose that $a < u < b$ or $b < u < a$. In the case $b < u < a$, assume further that $2a+\alpha^2\leq 2b+\beta^2$. Then

Proof. We will prove the length bound for $\Gamma_i$ only, as the one for $\Gamma_o$ is proven in the same way (up to some sign differences coming from the fact that $\nabla\Psi$ is pointing outside $\Omega$ on $\Gamma_o$). Let us also suppose for the moment that $\alpha\neq 0$, so that the gradient of $u$ does not vanish on $\Gamma_i$. Recalling that $\Psi = G^{-1}(u)$, where $G(x) = L-x^2/2+ M\log x$, we easily compute

In particular the norm $| \nabla\Psi| = \sqrt{W/W_0}$ is equal to $1$ on $\Gamma_i$ by construction and is less than $1$ inside $\Omega$ thanks to (3.3). In particular, the inequality $\langle \nabla| \nabla\Psi|^2\, |\, \nabla\Psi\rangle\leq 0$ holds on $\Gamma_i$ and the curvature $\kappa$ of $\Gamma_i$ at any point of $\Gamma_i$ satisfies

Integrating $\kappa\geq r_i$ on $\Gamma_i$, recalling that the total curvature of a simple curve is equal to $2\pi$, we get

This concludes the proof of the length bound for $\Gamma_i$ when $\alpha\neq 0$.

If instead $\alpha = 0$, we have $W = W_0 = 0$ on $\Gamma_i$, so $| \nabla\Psi| = \sqrt{W/W_0}$ is not well defined and we need to find a workaround. Let $p\in\Gamma_i$ and let us extend $u$ to a smooth function in a neighborhood $U$ of $p$. We can do it by means of Whitney Extension Theorem ^[39], that clearly applies to our case since $u$ is smooth up to the boundary (see for instance [24,Theorem 6.19]). With a slight abuse of notation, let us still call $u$ the smooth function defined on $U$ that coincides with the original $u$ in $\overline{\Omega}\cap U$. In the spirit of [11,Theorem 3.1], we consider the signed distance function $r$ to $\Gamma_i$. Namely, $|r(x)| = {{\rm{dist}}}(x, \Gamma_i)$ and $r > 0$ inside $\Omega$, $r < 0$ outside $\Omega$. If $U$ is taken to be small enough, the function $r$ is known to be smooth with $| \nabla r|\equiv 1$ in $U$. We can then proceed exactly as in [11,Theorem 3.1] to show that $u = a+(\Delta u/2) r^2+r^3 \omega$ in $U$, where $\omega$ is a smooth function. In particular, in $\overline{\Omega}\cap U$ we have

Recalling $\Delta u = -2$, from this expansion we easily compute that $\Delta r = 3\omega$ on $\Gamma_i$. Furthermore, recalling that $u = G(\Psi) = L-\Psi^2/2+r_i^2\log \Psi$ (notice that, if $\alpha = 0$, from (2.3) we have $M = r_i^2$), starting from the Taylor series for $G$ at $r_i$, with some computations one also gets an expansion for $\Psi$ in $\overline{\Omega}\cap U$:

Computing $| \nabla\Psi|$ from this expansion using $| \nabla r| = 1$ in $U$ and recalling that $| \nabla\Psi|\leq 1$ thanks to the gradient estimates, we deduce $3\omega\leq 1/r_i$. We now use the fact that $\Gamma_i = \{r = 0\}$ is a level set of $r$ to compute its curvature as follows

Integrating on $\Gamma_i$ we then obtain the desired length bound.

3.3.   Proof of Theorem 1.2

In this subsection, we focus on the case $a < b$ and we prove Theorem 1.2. For simplicity, let us start by discussing the case where $\nabla u$ does not vanish on the boundary of $\Omega$, namely for the moment let us assume $\beta\neq 0$ (notice that we know that $\alpha$ is strictly positive in the case $a < b$, so $\alpha$ cannot vanish). Using the divergence theorem and recalling that $\Delta u = -2$, we compute

where $\nu$ is the outward unit normal to $\partial\Omega = \Gamma_i\sqcup\Gamma_o$. Since we are in the case $a < b$, then $\nu = - \nabla u/| \nabla u|$ on $\Gamma_i$ and $\nu = \nabla u/| \nabla u|$ on $\Gamma_o$. Setting $W = | \nabla u|^2$ as usual, and recalling the expression (3.2) of $W_0$ in terms of $\Psi$, the identity above can be written as

Notice that the left hand side is nonnegative because of (3.3). Furthermore, since we are focusing on the case $a < b$, we have $\Psi < r_o < \sqrt{ M}$ hence

Since by construction $W = W_0$ on $\partial\Omega$, from (3.10) we immediately deduce

On the other hand, the length bounds in Proposition 3.2 give us the opposite inequality. Therefore the equality must hold in (3.11), and from (3.10) in particular we obtain that $W\equiv W_0$ in $\Omega$. The result now follows from the rigidity statement in Theorem 3.1.

It remains to discuss the case where $\beta$ is equal to zero. In this case, the boundary term in (3.9) is ill-defined, as one has $\Psi^2 = r_o^2 = M$ on $\Gamma_o$. To avoid this problem, we just need to apply the divergence theorem in $\Omega\setminus\{u > b- \varepsilon\}$, and then take the limit as $\varepsilon\to 0$. From (3.8) it follows easily that the quantity $\langle \nabla u\, |\, \nu\rangle/(M-\Psi^2)$ goes to $1/\Psi = 1/r_o$ as we approach $\Gamma_o$. The proof then proceeds exactly as in the case $\beta\neq 0$.

4.   Exploiting the Pohozaev identity

In this section we focus on the case $a > b$ and we prove Theorem 1.1 under the additional hypothesis that $2a+\alpha^2\leq 2b+\beta^2$. Let us give the precise statement for the convenience of the reader.

Theorem 4.1. Let $(\Omega, u)$ be a solution to problem $(1.1)$ such that $a > b$, $b < u < a$ in $\Omega$ and $2a+\alpha^2\leq 2b+\beta^2$. Then $\Omega$ is an annulus and $u$ is rotationally symmetric with $\partial u/ \partial |x| < 0$.

As discussed, all the arguments in Subsections 3.1 and 3.2 still work in this case. Unfortunately however, the strategy employed in Subsection 3.3 to complete the proof for $a < b$, does not work when $a > b$. In fact, one can check that in this case the inequality that one gets at the end has the opposite sign of (3.11), hence it does not combine with the length bounds from Proposition 3.2 anymore. We will then need to exploit a more delicate argument based on sharp estimates coming from the Pohozaev identity and the isoperimetric inequality.

Let us start by writing

where

for $k\in \mathbb{R}$ constant to be chosen later. The function $\varphi$ is constant on the level sets of $u$, hence we can consider the function $\dot \varphi$ representing the derivative of $\varphi$ with respect to $u$. More precisely, $\dot \varphi:\Omega\to \mathbb{R}$ is the function satisfying $\nabla \varphi = \dot \varphi \nabla u$. We then use the divergence theorem to write the following sequence of identities

where we have denoted by $\varphi_i$ and $\varphi_o$ the (constant) value of $\varphi$ on $\Gamma_i$, $\Gamma_o$, respectively.

The choice of $\varphi$ in (4.1) has been done so that it holds $(2 \varphi-4u)/\dot \varphi = W_0$. Combining (4.2) with the Pohozaev identity (2.1), we obtain

Since $|E_o| = |E_i|+|\Omega|$, we can write the above equation in the following way:

Recall from Proposition 2.1 that $4a+\alpha^2 > 4b+\beta^2$. Furthermore, combining the isoperimetric inequality with the length bound in Proposition 3.2 we get

Finally, we have the identity

Combining these insights together, we obtain the following inequality:

Recall that the constants appearing in the above formula are related as follows:

Using these identities and formula (4.1) for $\varphi$, with some computations we can rewrite the above inequality as

Let us now analyze $\dot \varphi$. Differentiating (4.1), we get

One can check that the quantity in square brakets is monotonically decreasing in $\Psi$. Since we know from Lemma 2.3 that $\Psi^2\geq r_i^2\geq M$, in order for $\dot \varphi$ to be positive in $\Omega$ it is then sufficient to choose

Choosing $k$ exactly equal to the above value, recalling from (3.3) that $W\leq W_0$, we then get from (4.3):

Since $r_i^2\geq M > 0$ and $4b+\beta^2-4a-\alpha^2 < 0$, this implies

But we have the opposite inequality from Proposition 3.2, therefore everything must be an equality. As a consequence, $W = W_0$ on the whole $\Omega$ and we conclude using the rigidity statement in Theorem 3.1.

5.   Further comments

In the previous sections we have shown that we are able to deal with the case $a > b$ and prove Theorem 1.1 if we assume the additional hypothesis $2a+\alpha^2\leq 2b+\beta^2$. While this hypothesis is somewhat restrictive, it still allows to get some nontrivial applications. As a particularly relevant example, we now argue that our weaker version of Theorem 1.1, namely Theorem 4.1, is still strong enough to deal with the cases of interest in ^[2]. There, Theorem 1.1 was invoked in the proof of [2,Theorem B] on a domain $\Omega_o$ (the outer domain, in the terminology of ^[2]) for a function $u$ satisfying problem (1.1) with $\alpha = 0$, $b = 0$ and $\beta^2/2a\geq 1$ (the latter inequality followed from [2,Theorem 2.1]). Under those hypotheses, the inequality $2a+\alpha^2\leq 2b+\beta^2$ is trivially satisfied. In other words, Theorem 4.1 is enough for the intended applications in ^[2].

The rest of this section is devoted to further comment on the hypothesis $2a+\alpha^2\leq 2b+\beta^2$. Specifically, we discuss what can go wrong when this hypothesis is not in place and what are the complications that one would need to overcome in order to prove Theorem 1.1 in full generality using our approach.

The first basic step that one needs in order to start our comparison argument is of course to have a model to compare with. In other words, given a solution $(\Omega, u)$ of problem (1.1) for some $a, b, \alpha, \beta$ with $a > b$, we need to have a model solution (2.2) solving the same problem. This is granted by Lemma 2.3 when $2a+\alpha^2\leq 2b+\beta^2$. Unfortunately, an analogous result seems harder to prove when $2a+\alpha^2 > 2b+\beta^2$.

In case we are able to show that there is a model solution to compare with, there is still another crucial complication, namely the proof of Theorem 3.1. In fact, in order to prove the gradient estimate $W\leq W_0$, we relied on the Maximum Principle applied to the elliptic inequality (3.5). Unfortunately, when $M < 0$ (this can happen when $2a+\alpha^2 > 2b+\beta^2$), the zeroth order term of (3.5) has the wrong sign, hence we cannot apply the Maximum Principle. It is possible with some work to find a workaround, at least under the additional hypothesis that $r_i^2\geq- M$.

Theorem 5.1 (Gradient Estimates, $M < 0$). Let $(\Omega, u)$ be a solution to problem $(1.1)$ and suppose that $b < u < a$. Suppose that there exist $L, M, r_o, r_i$ such that the corresponding model solution $(2.2)$ solves the problem for the same $a, b, \alpha, \beta$. Let $W$, $W_0$ be defined as in $(3.2)$. If $r^2_i\geq - M > 0$, then it holds

Moreover, if $W = W_0$ at some point in the interior of $\Omega$, then $(\Omega, u)$ is rotationally symmetric.

Proof. We start by observing that we can write (3.4) as

If we ask the quantity in the last square bracket to be equal to zero, we obtain the following formula for $\gamma$

where $k\in \mathbb{R}$ is a constant. In particular we get

Since $\Psi^2\geq r_i^2\geq - M > 0$, choosing $k$ to be big enough we have that the right hand side of (5.3) is negative at any point in $\Omega$. As a consequence, for such $k$, from (5.2) we get

We can then apply the Maximum Principle to find out that $F_\gamma$ is nonpositive on the whole $\Omega$. This concludes the proof.

The question still remains open of whether it is possible to remove the hypothesis $r_i^2\geq - M$ in the above theorem. Looking at the proof, it is clear that the freedom in the choice of $k$ should allow to weaken that hypothesis. However, it does seem hard to remove it entirely.

The rest of our arguments does not seem to depend heavily on the hypothesis $2a+\alpha^2\leq 2b+\beta^2$. In other words, even in the case $2a+\alpha^2 > 2b+\beta^2$, our method should allow to conclude the desired rigidity, provided we have a model solution to compare with and the gradient estimate $W\leq W_0$ is in force. In order to prove Theorem 1.1 in full generality using our approach, these seem to be the two issues that are left to address.

Acknowledgments

The author is a member of the Gruppo Nazionale per l'Analisi Matematica, la Probabilità e le loro Applicazioni (GNAMPA) of the Istituto Nazionale di Alta Matematica (INdAM).

Conflict of interest

The author declares no conflict of interest.