Hölder stability for Serrin's overdetermined problem
In a bounded domain $\Omega$, we consider a positive solution of the problem $\Delta u+f(u)=0$ in $\Omega$, $u=0$ on $\partial\Omega$, where $f:\mathbb{R}\to\mathbb{R}$ is a locally Lipschitz continuous function. Under sufficient conditions on $\Omega$ (for instance, if $\Omega$ is convex), we show that $\partial\Omega$ is contained in a spherical annulus of radii $r_i<r_e$, where $r_e-r_i\leq C\,[u_\nu]_{\partial\Omega}^\alpha$ for some constants $C>0$ and $\alpha\in (0,1]$. Here, $[u_\nu]_{\partial\Omega}$ is the Lipschitz seminorm on $\partial\Omega$ of the normal derivative of $u$. This result improves to Hölder stability the logarithmic estimate obtained in [1] for Serrin's overdetermined problem. It also extends to a large class of semilinear equations the Hölder estimate obtained in [6] for the case of torsional rigidity ($f\equiv 1$) by means of integral identities. The proof hinges on ideas contained in [1] and uses Carleson-type estimates and improved Harnack inequalities in cones.
\section[Introduction]{Introduction}
Serrin's overdetermined problem has been the object of many investigations.
In its classical form, it involves a sufficiently smooth bounded domain $\Omega$ in ${\mathbb R^N}$ and a
classical solution of the set of equations:
\begin{eqnarray}
&\Delta u + f(u) = 0 \ \mbox{ and } \
u\ge 0 \ \mbox{ in } \ \Omega, \ u=0 \ \mbox{ on } \ \partial \Omega,
\label{serrin1}
\\
&u_\nu=c \ \mbox{ on } \ \partial \Omega.
\label{serrin2}
\end{eqnarray}
Here, $f:[0,+\infty)\to\mathbb R$ is a locally Lipschitz continuous function, $u_\nu$ denotes the {\it inward} normal derivative of $u$ and $c$ is a positive constant. Under these assumptions, Serrin \cite{Se} proved that $\Omega$ must be a ball and $u$ be radially symmetric. For his proof, he adapted and improved a method created by Aleksandrov to prove his {\it Soap Bubble Theorem} (see \cite{Al}). 
In fact, the radial symmetry of $\Omega$ and $u$ is obtained by the so-called {\it method of moving planes}
based upon the observation that the euclidean ball is the only bounded domain that is symmetric with respect to any hyperplane passing through its center of mass. There are other interesting proofs of this symmetry result, based on integral identities, that generally need severe restrictions on $f$ (see for instance \cite{We}, \cite{PS}, \cite{BNST}).
\par
Overdetermined problems like \eqref{serrin1}-\eqref{serrin2} arise in many physical and geometric situations; they can be seen as a prototype of inverse problems and often emerge in free boundary and shape optimization problems (see \cite{He}). 
\par
Despite the intense research that has been devoted to them for almost five decades, there are still many open problems. An important one  --- the focus of this paper --- concerns the study of the stability of the radial configuration.  
\par
The first contribution
in this direction is \cite{ABR}, where
it is proved that if one assumes that $u_\nu$ is \emph{almost} constant on $\partial \Omega$, then there exist two concentric balls $B_{r_i}$ and $B_{r_e}$, with 
\begin{equation}
\label{concentric balls}
B_{r_i} \subset \Omega \subset B_{r_e}
\end{equation} 
and such that 
$r_e-r_i$ can be bounded in terms of some measure of the deviation of $u_\nu$ from being a constant. More precisely, in \cite{ABR} it is proved the estimate
\begin{equation} 
\label{stability ABR}
r_e-r_i\le C\,|\log  \left\|u_\nu-d\|_{C^1(\partial \Omega)} \right|^{-1/N},
\end{equation}
where $d$ is some given constant, provided $\| u_\nu-d\|_{C^1(\partial \Omega)}$ is sufficiently small; here, $C$ is a constant depending on $N$, the regularity of $\partial\Omega$, the diameter of $\Omega$ and the Lipschtz constant of $f$.
The proof is based on a quantitative study of the method of moving planes and works for a general locally Lipschitz non-linearity $f$.
\par
In the case of the {\it torsional rigidity} problem, that is when $f\equiv 1$, \eqref{stability ABR} was improved in \cite{BNST2} (see also \cite{BNST3} for Monge-Amp\`ere equations). Indeed, the authors replace the logarithmic dependence at the right-hand side of \eqref{stability ABR} 
by a power law of H\"older type. Furthermore, they also give a stability estimate in terms of the $L^1$-norm of the deviation instead of its $C^1$-norm.
\par
In our main result, Theorem \ref{stability serrin}, we will show that the logarithmic estimate
\eqref{stability ABR} can be improved to obtain a stability of H\"older type. In order to avoid unnecessary technicalities, in this section we present our result in a relevant particular case. In what follows, we denote by $d_\Omega$ the {\it diameter} of $\Omega$, $r_\Omega$ the radius of the optimal interior touching sphere to $\partial\Omega$ (see Section \ref{sec:MovPlanes} for its precise definition), and use the following notation:
\begin{equation}
\label{seminorm normal derivative}
\big[u_\nu \big]_{\pa \Om} = \sup_{\substack{x,y \in \partial \Omega\\ \ x \neq y}} \frac{|u_\nu(x) - u_\nu(y)|}{|x-y|}.
\end{equation}



\begin{theorem}
\label{th:stability-convex}
Let $\Omega \subset \mathbb R^N$ be a convex domain with boundary of class $C^{2,\alpha}$
and let $f$ be a locally Lipschitz continuous function such that $f(0)\ge 0$.  
Let $u\in C^{2,\alpha}(\overline{\Omega})$ be a solution of \eqref{serrin1}.
\par
There exist two positive numbers $\varepsilon$, $C$, that depend on 
\begin{eqnarray*}
\displaystyle N, f, d_\Omega, r_\Omega, \mbox{the $C^{2,\alpha}$-regularity of } \Omega, \max_{\overline{\Omega}}u, \min_{\partial\Omega} u_\nu, 
\end{eqnarray*}
and a number $\tau\in (0,1)$ such that, if $[u_\nu]_{\partial \Omega} \leq \varepsilon$,
then 
\begin{equation*}
B_{r_i} \subset \Omega \subset B_{r_e},
\end{equation*}
where $B_{r_i}$ and $B_{r_e}$ are two concentric balls and their radii satisfy 
\begin{equation} 
\label{stability-improved}
r_e-r_i \leq C\,[u_\nu]_{\partial \Omega}^{\tau}.
\end{equation}
\end{theorem}
 The number $\tau$ can be explicitly determined (see Theorem \ref{stability serrin}, for details). 
Notice that the semi-norm \eqref{seminorm normal derivative} can be bounded in terms of the deviation used in \eqref{stability ABR}, if we choose $d$ as the minimum of $u_\nu$ on $\partial\Omega$.
\par
Therefore, \eqref{stability-improved} significantly improves the estimates of \cite{ABR} and \cite{BNST2}, since it enhances the stability from logarithmic, as in \eqref{stability ABR}, to that of H\"older type proved in \cite{BNST2}, but for {\it any locally Lipschitz non-linearity $f$} (that is, not only for the case $f\equiv 1$).
\par
The assumption on the convexity of $\Omega$ can be slightly relaxed (see Theorem \ref{stability serrin}), by requiring that $\Omega$ be what we call a $(\mathcal{C},\theta)${\it -domain}. Roughly speaking, we require that every maximal cap that comes about in employing the method of moving planes has a boundary with a Lipschitz constant bounded uniformly with respect to the direction of mirror reflection chosen (see Section \ref{sec:MovPlanes} for details). 
\par
As already mentioned, our results and the ones in \cite{BNST2} are obtained by using very different techniques. The starting point of \cite{BNST2} is the proof of symmetry given by the same authors in \cite{BNST} which makes use of information, such as Newton inequalities and Pohozaev identity, which is in a sense more global and avoids an extensive use of maximum principles needed to employ the method of moving planes. On the other hand, the method in \cite{BNST2} is more restrictive, since it seems to be suitable only for $f$ constant. 
\par
Instead, based on Serrin's original proof,  our approach is more flexible and allows to treat a {\it general 
Lipschitz continuous non-linearity $f$}; in fact, our study relies on quantitative versions of the maximum principle, such as (global) Harnack-type inequalities. Thus, the quality of the relevant stability estimate is affected by that of the Harnack inequality we employ. This is the reason why, at this stage, even if we can consider a general non-linearity, we need to put a restriction on the type of domain under study. In particular, for convex and $(\mathcal{C},\theta)${-domains}, we can use improved Harnack-type inequalities in every cap which is generated by the method of moving planes.
\par
Besides by those in \cite{ABR}, the techniques used in this paper are inspired by those employed in \cite{CMS2} (see also \cite{CMS}), where a quantitative study of the radially symmetric configuration was carried out for a related problem --- the {\it parallel surface problem} --- motivated by the remark, made in \cite{MSaihp} (see also \cite{MS1}, \cite{MSm2as}), that {\it time-invariant level surfaces} of solutions of certain nonlinear non-degenerate {\it fast diffusion} equations are parallel surfaces.
\par
As in \cite{ABR} and \cite{CMS2}, our approach consists in fixing a direction, defining an approximate set $X(\delta)$ --- built upon the so-called maximal cap (see Section \ref{sec:MovPlanes}) and its mirror-symmetric image in that direction --- which fits $\Omega$ well, as the parameter $\delta$ tends to $0$. This approximation process is controlled in terms of $\big[u_\nu \big]_{\pa \Om}$ and does not depend on the  particular direction chosen. 
\par
The application of Harnack's inequality and Carleson estimates in the maximal cap plays a crucial role in obtaining the new stability estimates. Since we are assuming that every maximal cap has Lipschitz regularity, the improvement in \eqref{stability-improved} is obtained by a refinement of Harnack's inequality in suitable cones.
\par
The paper is organized as follows. In Section \ref{sec:MovPlanes}, we recall the notations
and preliminaries necessary to use the method of moving planes. There, we also discuss on  
the definition of $(\mathcal{C},\theta)$-domains and present our main result, Theorem \ref{stability serrin}, 
of which Theorem \ref{th:stability-convex} is a straightforward corollary.
In Section \ref{sec:enhanced harnack} --- the core of the paper --- we begin the proof of Theorem \ref{stability serrin}, by producing the necessary enhanced Harnack estimates. Finally, in Section \ref{sec:stability}, we complete the proof of \eqref{stability-improved}.

\setcounter{equation}{0}

\section{Remarks on the method of moving planes and\\ statement of the main result} 
\label{sec:MovPlanes}
We consider a bounded domain $\Omega$ of class $C^{2,\alpha}$, $0<\alpha \leq 1$. 
In what follows we will often use the notation $C(N, \Omega, d_\Omega, r_\Omega, f, \dots)$ to denote a constant that
depends on the relevant parameters; in particular, the dependence
on $\Omega$ is meant to be only on the $C^{2,\alpha}$ regularity of $\partial\Omega$, as explained in \cite[Remark 1]{ABR}.
The assumed regularity of $\Omega$ implies that $\Omega$ satisfies a {\it uniform interior sphere condition}; for $x \in \partial \Omega$, we let $r(x)$ be the radius of the maximal ball $B \subset \Omega$ with $x \in \partial B$ and set
$$
r_\Omega = \min_{x \in \partial \Omega} r(x).
$$ 
\par
In the sequel, it will be useful to consider the {\it parallel set} of $\Omega$:
\begin{equation} \label{Om de}
\Omega(\delta) = \{x \in \Omega:\ \mbox{\rm dist}(x,\partial \Omega)>\delta\}.
\end{equation}
If $0<\delta<r_\Omega$, $\partial\Omega(\delta)$ is $C^{2,\alpha}$-smooth; also, $\Omega(\delta)$ satisfies a uniform interior sphere condition with optimal radius $r_\Omega-\delta$.
 \par
For a unit vector $\omega\in{\mathbb R^N}$ and a parameter $\mu\in\mathbb R$, we define the following objects:
\begin{equation}
\label{definitions}
\begin{array}{lll}
&\pi_{\mu}=\{ x\in{\mathbb R^N}: x\cdot\omega=\mu\}\ &\mbox{a hyperplane orthogonal to $\omega,$}\\
& \mathcal{H}_\mu=\{ x\in{\mathbb R^N}: x\cdot\omega>\mu\}\ &\mbox{the half-space \emph{on the right} of $\pi_\mu$,}\\
&\Omega_{\mu}=\{x\in A: x\cdot\omega>\mu\}\ &\mbox{the right-hand cap of $\Omega$},\\
&x^{\mu}=x-2(x\cdot\omega-\mu)\,\omega\ &\mbox{the reflected image of $x$ in $\pi_{\mu},$}\\
&\Omega^{\mu}=\{x\in{\mathbb R^N}:x^{\mu}\in \Omega_{\mu}\}\ &\mbox{the reflected cap in $\pi_{\mu}$.}
\end{array}
\end{equation}
Set $\Lambda=\sup\{x\cdot\omega: x\in \Omega\}$, the {\it extent} of $\Omega$ in direction $\omega$; if $\mu<\Lambda$ is close to $\Lambda$, the reflected cap $\Omega^\mu$ is contained in $\Omega$ (see \cite{Fr}). Set
\begin{equation}
\label{lambda def}
\lambda=\inf\{\mu: \Omega^{\mu'}\subset \Omega \mbox{ for all } \mu'\in(\mu,\Lambda)\}.
\end{equation}
Then at least one of the following cases occurs (see \cite{Se}, \cite{Fr}):
 \begin{enumerate}
\item[(S1)]
$\Omega^{\lambda}$ becomes tangent to $\partial \Omega$ at some
point $P^\lambda\in\partial \Omega^\lambda \setminus\pi_{\lambda}$, that is the reflected image of a point $P \in \partial \Omega_\lambda \setminus\pi_{\lambda}$;
\item[(S2)]
$\pi_{\lambda}$ is orthogonal to $\partial \Omega$ at some point $Q\in\partial \Omega\cap\pi_{\lambda}.$
\end{enumerate}
The cap $\Omega_\lambda$ will be called the {\it maximal cap}.
\par
As customary in the method of moving planes, we define the function
\begin{equation}\label{w def}
w(x)= u(x^\lambda)-u(x),\quad x\in \Omega_\lambda;
\end{equation}
$w$ satisfies the equation
\begin{equation} 
\label{eq w semilinear}
\Delta w+c(x)\,w=0 \  \mbox{ in } \Omega_\lambda,
\end{equation}
where for $x\in\Omega_\lambda$
\begin{equation*}
\label{defc}
c(x)=\left\{
\begin{array}{lll}
\displaystyle\frac{f(u(x^\lambda))-f(u(x))}{u(x^\lambda)-u(x)} &\mbox{ if } u(x^\lambda)\not= u(x),\\
\displaystyle 0 &\mbox{ if } u(x^\lambda) = u(x).
\end{array}
\right.
\end{equation*}
Notice that $c(x)$ is bounded by the {\it Lipschitz constant} $\frak{L}_f$ of $f$ in the interval
$[0,\max\limits_{\overline{\Omega}}u].$
\par
All the improved estimates in Section \ref{sec:enhanced harnack} concern $w$. As proved in \cite{Se} and refined in \cite{BNV} (see also \cite{Fr}), since $w\ge 0$ on $\partial\Omega_\lambda$, we can assume that $w\ge 0$ in $\Omega_\lambda$. Hence, a standard application of the strong maximum principle to the inequality $\Delta w-c^-(x)\, w\le 0$ with $c^-(x)=\max[-c(x),0]$ shows that either $w=0$ in $\Omega_\lambda$ (and $\Omega$ and $u$ are symmetric about $\pi_\lambda$) or 
$$
w>0 \ \mbox{ in } \ \Omega_\lambda.
$$

The following lemma ensures that the maximal cap always contains a half ball tangent to $\partial\Omega$
at either point $P$ or $Q$.



\begin{lm} \label{lemma ball ABR}
Let $P$ and $Q$ be as in case (S1) and (S2), respectively. Let $B_\rho(p)\subset\Omega$ be a ball with
$0<\rho\le r_\Omega$ and such that $P\in\partial B_\rho(p)$ or $Q\in\partial B_\rho(p)$.
\par
Then, $p\in \overline{\Omega}_\lambda$ and $B_\rho(p) \cap \mathcal{H}_\lambda \subset \Omega_\lambda$.
\end{lm}
\begin{proof}
The assertion is trivial for case (S2). 
\par
If case (S1) occurs, without loss of generality, we can assume that $\omega=e_1=(1, 0,\dots, 0)$ and $\lambda=0$. Since (S1) holds, the point $P^\lambda$ lies on $\partial\Omega$ and cannot fall inside $B_\rho(p)$, since $P\in \partial B_\rho (p)$ and $B_\rho (p)\subset\Omega$. Thus,  $|p-P^\lambda|\geq \rho=|p-P|$ and hence 
$|p_1+P_1| \geq |p_1-P_1|$, which implies that $p_1 \geq 0$, being $P_1>0$.
\end{proof}

As mentioned in the Introduction, the convexity assumption of Theorem \ref{th:stability-convex} can be relaxed; with this purpose, we introduce the class of $(\mathcal{C}, \theta)$-domains that, roughly speaking, have the property that every maximal cap $\Omega_\lambda$ has a boundary with a Lipschitz constant bounded by a number which does not
depend on the direction of reflection chosen.
To this aim, for a fixed $0<t<1/2$, by formula \eqref{Om de}, we define the set
\begin{equation}\label{G def}
G= \Omega(t\,r_\Omega);
\end{equation}
we know that $G$ is connected (see \cite[pp. 923--924]{ABR}).
Thus, we say that $\Omega$ is a {\it $(\mathcal{C}, \theta)$-domain} if 
there exists $\theta > 0$ such that for any direction $\omega$ and for any $x\in \Omega_\lambda\setminus\overline{G}_\lambda$ 
there exists $\zeta \in \partial G_\lambda \setminus \pi_\lambda$ such that $x$ and $\zeta$ belong to the axis of a (finite) right spherical cone $\mathcal{C} \subset \Omega_\lambda$ with vertex at $x$ and aperture $2\theta$.
\par
This property is not easy to check, since it relies on the knowledge of the position of the critical hyperplane (see Fig.\ref{fig1}). However, the class of $(\mathcal{C}, \theta)$-domains is not empty, as shown by the following proposition.

\begin{figure}[h]
\centering
 \subfigure[A $(\mathcal{C}, \theta)$-domain.]{\includegraphics[width=0.5\textwidth]{ok_Cth.png}}
 \hspace{2em}
 \subfigure[Not a $(\mathcal{C}, \theta)$-domain.]{\includegraphics[width=0.4\textwidth]{no_Cth.png}}
  \caption{Every maximal cap of a $(\mathcal{C}, \theta)$-domain has Lipschitz boundary.} \label{fig1}
\end{figure}


\begin{prop} 
\label{proposition cone}
Any convex domain $\Omega$ of class $C^1$ and satisfying the uniform interior sphere condition with
minimal radius $r_\Omega$ is a $(\mathcal{C}, \theta)$-domain with
\begin{equation}
\label{theta}
\theta\ge \arctan \frac{(1-t)\,r_\Omega}{2\,d_\Omega},
\end{equation}
where $t$ is the parameter appearing in \eqref{G def}.
\end{prop}


\begin{proof}
Since $\Omega$ is of class $C^1$, then the method of moving planes can be applied (see \cite{Fr}). As already observed, $G$ satisfies a uniform interior sphere condition with optimal radius $(1-t)\,r_\Omega$; also, 
\begin{equation*}
\Omega=\{ x+y: x\in G,\, |y|<t\, r_\Omega\},
\end{equation*}
and $G$, $G_\lambda$ and $\Omega_\lambda$ are all convex, since $\Omega$ is convex.

From Lemma \ref{lemma ball ABR} we have that $\Omega_\lambda$ contains a half-ball of radius $r_\Omega$ with center $p\in \overline{\Omega}_\lambda$, and $G_\lambda$ contains the half-ball of radius $(1-t)\,r_\Omega$ centered at $p$. The maximal ball $B$ contained in this half ball has radius $\overline{r}=(1-t)\,r_\Omega /2$ and is contained in $G_\lambda$; we denote its center by $y$.
\par
Let $x$ be any point in $\Omega_\lambda \setminus \overline{G}_\lambda$ and let $\mathcal{C}$ be the (finite) right circular cone with vertex at $x$ and based on the $(N-1)$-dimensional ball of radius $\overline{r}$ obtained by intersecting $B$ with the hyperplane perpendicular to the vector $x-y$ and passing through $y$. Since $\Omega_\lambda$ is convex, then $\mathcal{C} \subset \Omega_\lambda$. Moreover, being $B\subset G_\lambda$, the axis
of $\mathcal{C}$ intersects $\partial G_\lambda$ at a point $\zeta \notin\pi_\lambda$.
\par
It is clear that the aperture $2\theta$ of $\mathcal{C}$ is such that
$
\theta \geq  \arctan(\overline{r}/d_\Omega),
$
and hence \eqref{theta} holds.
\end{proof}

\begin{corollary}
Let $\Omega$ be a convex domain and $\partial \Omega$ be of class $C^1$. 
\par
Then any maximal cap $\Omega_\lambda$  has a Lipschitz continuous boundary  $\partial\Omega_\lambda$ with Lipschitz constant bounded by $2d_\Omega / r_\Omega$.
\end{corollary}
\par
We now state our main result.
\begin{theorem} 
\label{stability serrin}
Let $\Omega \subset \mathbb R^N$ be a $(\mathcal{C}, \theta)$-domain with boundary of class $C^{2,\alpha}$
and let $f$ be a locally Lipschitz continuous function with $f(0)\ge 0$.  
\par
There exist two positive constants $\varepsilon$ and $C$, that depend on 
$$
N, \, f, \, d_\Omega, \, r_\Omega, \, \mbox{the $C^{2,\alpha}$-regularity of } \Omega, \, \max_{\overline{\Omega}} u, \, \min_{\partial\Omega} u_\nu,
$$
and a number $\tau\in (0,1)$  such that, if $[u_\nu]_{\partial \Omega} \leq \varepsilon$,
then
\begin{equation*}
B_{r_i} \subset \Omega \subset B_{r_e},
\end{equation*}
where $B_{r_i}$ and $B_{r_e}$ are two concentric balls and their radii satisfy \eqref{stability-improved}:
\begin{equation*} 
r_e-r_i \leq C\,[u_\nu]_{\partial \Omega}^\tau.
\end{equation*}
\end{theorem}
\begin{remark}
(a) As it will be clear from the proof, it holds that $\tau=1/(1+\gamma)$ where $\gamma$ is given by \eqref{gamma}; $\gamma$ depends on the half-aperture $\theta$ of the cone $\mathcal{C}$ in the definition 
of $(\mathcal{C},\theta)$-domain and on the Harnack's constant for equation \eqref{eq w semilinear}.
\par
(b) It is clear that Theorem \ref{th:stability-convex} is an easy corollary of this theorem, by Proposition \ref{proposition cone}.
\par
(c) As explained in \cite{ABR}, in Theorem \ref{stability serrin} some precautions are in order. 
The dependence of $\varepsilon$ and $C$ on a lower bound for $u_\nu$ is needed in case $f(0)=0$.
In fact, if $f(0)>0$ a comparison argument in an interior ball touching $\partial\Omega$ shows that the minimum of $u_\nu$ on $\partial \Omega$ can be bounded from below by $f(0)$ times a constant that only depends on $N, \Omega, \| u\|_\infty$, and $f$. 
When $f(0)=0$, instead, the first inequality in \eqref{estimate-distance} below does not hold and the constant $\varepsilon$ (and hence $C$) must depend on a lower bound for $u_\nu$. This fact can be seen by considering any (positive) multiple of the first Dirichlet eigenfunction $\phi_1$ for $-\Delta$: in fact, for any $n\in\mathbb N$ the function $\phi_1/n$ satisfies \eqref{serrin1} with $f(u)=\lambda_1\,u$,
being $\lambda_1$ the first Dirichlet eigenvalue;
although $(\phi_1/n)_\nu\to 0$ on $\partial\Omega$ as $n\to\infty$, one cannot expect to derive any
information on the shape of $\Omega$.
\par
The question whether  Theorem \ref{stability serrin} or \cite[Theorem 1]{ABR} still hold if $f(0)<0$
remains open (see \cite[Section 5.1]{ABR}).
\end{remark}
\par
The proof of Theorem \ref{stability serrin} will be carried out in Sections \ref{sec:enhanced harnack} and 
\ref{sec:stability}.  Since it is quite long and complex, for the reader's convenience, we give an outline which consists of three main steps.
\begin{itemize}
\item[(I)]
We first improve \cite[Proposition 1]{ABR}. We fix a direction $\omega$ and apply the procedure of moving planes. For a sufficiently small $\delta>0$, we consider the set $\Omega(\delta)$ in \eqref{Om de}
and show that there exists a connected component $\Sigma_\delta$ of $\Omega(\delta) \cap \mathcal{H}_\lambda$ such that
$$
\|w\|_{L^\infty(\Sigma_\delta)} \leq \frac{C}{\delta^\gamma} \big[u_\nu \big]_{\pa \Om}\, ,
$$ 
where $C$ and $\gamma$ do not depend on $\delta$ and $\omega$ (see Proposition \ref{prop:new-prop-1} below).

\item[(II)]
Let $X(\delta)$ be the union of $\Sigma_\delta$ and its reflected image in the critical hyperplane $\pi_\lambda$. Since $u$ grows linearly near $\partial \Omega$, the smallness of $w$ in $\Sigma_\delta$ implies that $X(\delta)$ fits $\Omega$ well (see Lemma \ref{lemma 34 ABR} and Proposition \ref{prop 5 ABR} below). This step gives the approximate symmetry of $\Omega$ in the direction $\omega$.

\item[(III)]
Since the arguments in (I) and (II) do not depend on the chosen direction, we apply them for $N$ (mutually) orthogonal directions and obtain a point $\mathcal{O}$ as the intersection  of the corresponding $N$ critical hyperplanes. The point $\mathcal{O}$ can be chosen as an approximate center of symmetry. In fact, we are able to define two balls centered at $\mathcal{O}$ such that Theorem \ref{stability serrin} holds (see the proof of  Theorem \ref{stability serrin} in Section \ref{sec:stability} below). 
\end{itemize}

We notice that once one proves the first step, (I), the remaining ones, (II) and (III), follow from a well-established argument as in \cite{ABR} and \cite{CMS2}. 

\setcounter{equation}{0}

\section{Harnack's inequality in a cone}
\label{sec:enhanced harnack}

In this section, we prove some results based on Harnack's type inequality in a cone, which will be used in Section \ref{sec:stability} to prove step (I) (in particular Proposition \ref{prop:new-prop-1}). 

Let $0<a<1$ be fixed. It is well-known (see \cite[Theorem 8.20]{GT}) that a solution $w$ of \eqref{eq w semilinear} satisfies the following {\it Harnack's inequality}
\begin{equation} \label{harnack}
\sup_{B_{ar}} w \leq \frak{H}_a \inf_{B_{ar}} w ,
\end{equation}
for any ball $B_r \subset \Omega_\lambda$; the Harnack constant $\frak{H}_a$ can be bounded by the power 
$\sqrt{N} + \sqrt{r\,\| c\|_\infty}$ of a constant
only depending on $N$ and $a$ (see \cite{GT}). For instance, if $c(x) \equiv 0$, by the explicit Poisson's representation formula for harmonic functions, we have that
\begin{equation}\label{Harnack Harmonic}
\sup_{B_{ar}} w \leq \left( \frac{1+a}{1-a} \right)^N \inf_{B_{ar}} w ,
\end{equation}
for any $B_r \subset \Omega_\lambda$  (see \cite{GT}).



The following Lemma consists of an application of Harnack's inequality to a Harnack's chain of  balls 
contained in a cone. The result is well-known 
 (see \cite{Mo}, \cite{JK}, \cite{Ke}); however, since we are interested in a quantitative version of it, we provide our proof.


\begin{lm} 
\label{lemma harnack cone}
Fix a number $a\in (0,1)$. Pick $z \in \partial \Omega \cap \overline{\mathcal{H}}_\lambda$ and let $\mathcal{C}$ be any right spherical cone contained in $\Omega_\lambda$  and with vertex at $z$; let $2\theta$ be the aperture of $\mathcal{C}$.
\par
Let $w$ be given by \eqref{w def} and pick any two points $x$ and $\zeta$ on the axis of $\mathcal{C}$ such that
$$
|x-z|<|\zeta-z|.
$$
\par
Then we have that
\begin{equation}\label{Harnack x e y}
\frac{|x-z|^\gamma}{K}\, w (x) \leq w(\zeta)\, \leq \frac{K}{|x-z|^\gamma}\,w (x),
\end{equation}
where
\begin{eqnarray}
&& K=\frak{H}_a \,\left[\frac{|\zeta-z|(1-a\sin\theta)}{1-a}\right]^\gamma,\label{K} \\
&& \gamma = \log_\beta \frak{H}_a,  \ \ \beta= \frac{1+a\sin\theta}{1-a\sin\theta},\label{gamma}
\end{eqnarray}
and $\frak{H}_a$ is given by \eqref{harnack}.
\end{lm}



\begin{proof}
Here we prove the second inequality in \eqref{Harnack x e y}; the first one can be proved similarly.
\par
Let $\ell$ be the unit vector defining the axis of $\mathcal{C}$ through $z$, that is for instance
\begin{equation*}
    \ell=\frac{x-z}{|x-z|}.
\end{equation*}
We now construct a chain of balls $B_{r_i}(p_i)$, $i=0, 1, \dots, n$, joining $x$ to $y$ with the following
specifications:
\begin{itemize}
\item[(i)] the centers $p_0, p_1,\dots, p_n$ belong to the axis of $\mathcal{C}$;
\item[(ii)] $p_0=x$, $r_0=|x-z|$ and $\zeta\in B_{r_n}(p_n)$;
\item[(iii)] the balls $B_{r_1}(p_1), \dots, B_{r_n}(p_n)$ are all contained in $\mathcal{C}$ and tangent
to the lateral surface of $\mathcal{C}$;
\item[(iv)] the radii $r_1,\dots, r_n$ are chosen such that the balls $B_{ar_0}(p_0), \dots, B_{ar_n}(p_n)$ are pairwise disjoint and
\begin{equation*}
\overline{B_{ar_i}(p_i)} \cap \overline{B_{ar_{i+1}}(p_{i+1})} = \{p_i + a r_i \ell \},\ i=0,\ldots,n-1.
\end{equation*}
\end{itemize}
\par
A calculation shows that (i)-(iv) determine the $r_i$'s and $p_i$'s as
\begin{equation}
\label{ri&pi}
r_i=r_0\,\frac{(1-a) \sin\theta}{1-a\sin\theta}\,\beta^i, \quad p_i=z+r_0\,\frac{1-a}{1-a\sin\theta}\,\beta^i\,\ell , \ \ i=1, \dots, n.
\end{equation}
\par
Since $|p_{n-1}-z|\le |\zeta-z|$,  from the second formula in \eqref{ri&pi} we obtain a bound for $n$:
\begin{equation}
\label{bound-n}
n\le 1+\frac{\log_{\frak{H}_a}\left(\frac{|\zeta-z|}{r_0}\frac{1-a\sin\theta}{1-a}\right)}{\log_{\frak{H}_a}\beta}
\end{equation}
\par
As usual, the application of Harnack's inequality \eqref{harnack} to each ball $B_{r_i}(p_i)$ of the chain gives that
$
w(\zeta)\le\frak{H}_a^n\,w(x);
$
\eqref{Harnack x e y} then follows from \eqref{bound-n} and some algebraic manipulations.
\end{proof}


\begin{remark}
We notice that when $c(x) \equiv 0$ (and hence $\Delta w=0$) we have that
\begin{equation}
\label{gamma-torsion}
\gamma= N \frac{\log \frac{1+a}{1-a}}{\log\frac{1+a \sin\theta}{1-a \sin\theta}}.
\end{equation}
We observe that $\gamma \geq N$ and equality holds only if $\theta=\pi/2$.
\end{remark}

\setcounter{equation}{0}


\section{Stability for Serrin's problem: the proof}
\label{sec:stability}

In this section, thanks to the preparatory lemmas of Section  \ref{sec:enhanced harnack} and by following the outline announced in Section \ref{sec:MovPlanes}, we shall bring the proof of Theorem \ref{stability serrin} to an end.

\par
As already mentioned, the crucial step in the proof is item (I) in the outline; it is the analog of \cite[Proposition 1]{ABR}. We modify and improve the procedure used in \cite{ABR} at two salient points: (i) the definition of the approximating symmetric set $X(\delta)$; (ii)
the bound on the smallness of the function $w$ in \eqref{w def} in terms of the parameter $\delta$ and the semi-norm $\big[u_\nu \big]_{\pa \Om}$. 
We draw the reader's attention on the fact that, while (i) is due to a refinement,
based on Carleson-type estimates, of the technique used in \cite{ABR}, for (ii) the assumption that $\Omega$ is $(\mathcal{C}, \theta)$-domain is necessary --- and, so far, we are not able to remove it --- to treat the case of a general bounded domain with $C^{2,\alpha}$-smooth boundary.
More precisely, to every maximal cap that comes about in the method of moving planes, we want to apply Lemma \ref{lemma harnack cone} that, roughly speaking, requires that such maximal cap is Lipschitz continuous, which is essentially our definition of $(\mathcal{C}, \theta)$-domain, that is certainly fulfilled if $\Omega$ is convex, as shown in Proposition \ref{proposition cone}.  In other words, we are so far unable to exclude (even generically) that the situation depicted in Figure 1 (b) may occur.
\par
We now start the procedure.
For a fixed direction $\omega$, we let $\lambda$ be the number defined in \eqref{lambda def}; we then consider the connected component $\Sigma$ of $\Omega_\lambda$ which intersects the interior touching ball at the point $P$, if case (S1) occurs, or at the point $Q$, if (S2) occurs.
\par
For $\delta >0$, we consider the set $\Omega(\delta)$ defined in \eqref{Om de} and
define $\Sigma_\delta$ as the connected component of
$
\Omega(\delta) \cap \mathcal{H_\lambda}
$
contained in $\Sigma$ --- notice that our definition of $\Sigma_\delta$ differs from that in \cite{ABR}.
In addition, we fix the domain $G$ in \eqref{G def} by choosing $t=1/32$ and we set $a=1/2$ in Lemma \ref{lemma harnack cone}. 



\begin{prop}
\label{prop:new-prop-1}
Let $\Omega$ be a $(\mathcal{C}, \theta)$-domain. 
\par
Then, there is a constant $C=C(N, \Omega, d_\Omega, r_\Omega, f, \max_{\overline{\Omega}} u)$ such that
\begin{equation}
\label{w Linfty Sigma}
\|w\|_{L^\infty(\Sigma_\delta)} \leq C\, \delta^{-\gamma} [u_\nu]_{\partial \Omega}
\ \mbox{ for } \ 0<\delta \le r_\Omega/32, 
\end{equation}
where
$\gamma$ and $\beta$ are the numbers defined in \eqref{gamma}.
\end{prop}

\begin{proof}
Let $\widetilde{\Sigma}=\{x\in \Sigma:\ \mbox{\rm dist}(x,\partial \Sigma)> r_\Omega/64\}$. We apply  \cite[Proposition 1]{ABR} with $\delta=r_\Omega/64$ and obtain that 
\begin{equation*}
\|w\|_{L^\infty(\widetilde{\Sigma}_\delta)} \leq \widetilde{C}\,
[u_\nu]_{\partial \Omega},
\end{equation*}
for some constant $\widetilde{C}=\widetilde{C}(N, \Omega, d_\Omega, r_\Omega, f, \| u\|_\infty)$ . 
By using an argument analogous to the ones in the proofs of \cite[Lemmas 3.2 and 4.2]{CMS2} -- which are based on boundary Harnack's inequality -- we can extend the previous estimate to $G_\lambda$:
\begin{equation} \label{w Gm serrin}
\|w\|_{L^\infty(G_\lambda)} \leq M C\, [u_\nu]_{\partial \Omega},
\end{equation}
where $M$ is the constant appearing in the boundary Harnack inequality (see Theorem 1.3 in \cite{BCN}). It is important to notice that the bound in \eqref{w Gm serrin} does not depend on $\delta$.

Let $x$ be any point in $\Sigma_\delta \setminus \overline{G}_\lambda$, with $\delta < r_\Omega/32$. Since $\Omega$ is a $(\mathcal{C}, \theta)$-domain, from Lemma \ref{lemma harnack cone}, by choosing $\zeta$ equal to the point $\zeta\in \partial G_\lambda\setminus\pi_\lambda$ (in the definition of $(\mathcal{C}, \theta)$-domain)  
in Lemma \eqref{lemma harnack cone}, we obtain that
\begin{equation*}
    w(x) \leq K\, \delta^{-\gamma} w(\zeta);
\end{equation*}
\eqref{w Gm serrin} then yields \eqref{w Linfty Sigma} with $C=\widetilde{C} KM$.
\end{proof}

\begin{remark}
The statement of Proposition \ref{prop:new-prop-1} differs from that of \cite[Proposition 1]{ABR}
in three aspects: (i) the set $\Sigma_\delta$ we consider extends up to the hyperplane $\mathcal{H}_\lambda$; (ii) the dependence on $\delta$ in \eqref{w Linfty Sigma} greatly improves that of \cite[Eq. (8)]{ABR}, that blows up exponentially as $\delta\to 0$; (iii) the seminorm $[u_\nu]_{\partial\Omega}$ replaces the norm
$\| u_\nu-d\|_{C^1(\partial\Omega)}$. 
\end{remark}

Now, we define a symmetric open set as
\begin{equation}
\label{X def}
X(\delta)=\mbox{the interior of} \ \,\Sigma_\delta \cup \Sigma_\delta^{\lambda} \cup (\partial \Sigma_\delta \cap \pi_\lambda);
\end{equation}
we want to show that $X(\delta)$ \emph{fits $\Omega$ well}; how well is the main point of this paper.
\par
The main idea is to combine the estimate \eqref{w Linfty Sigma} on the smallness of $w$, together with the fact the $u$ grows linearly near $\partial \Omega$; as shown in  \cite[Proposition 4]{ABR} indeed, we know that
\begin{equation}
\label{estimate-distance}
\underline{K}\,\mbox{\rm dist}(x,\partial\Omega)\le u(x)\le \overline{K}\,\mbox{\rm dist}(x,\partial\Omega) \ \mbox{ for all} \ x\in\Omega,
\end{equation}
where 
$$
1/\underline{K}, \overline{K}\le C=C(\Omega,d_\Omega, \max_{\overline{\Omega}}u,f,\min_{\partial\Omega} u_\nu).
$$ 
In the following lemma we give our version  \cite[Eq.(34)]{ABR}.



\begin{lm} \label{lemma 34 ABR}
For $0<\sigma, \delta \leq r_\Omega/16$ let $\Omega(\sigma)$ and $X(\delta)$ be the sets defined in \eqref{Om de}
and \eqref{X def}, respectively.  Let $\gamma$ and $C$ be given by \eqref{gamma} and \eqref{w Linfty Sigma}, respectively. 
\par
If 
\begin{equation} \label{K si geq}
\underline{K}\,\sigma > C \delta^{-\gamma} [u_\nu]_{\partial \Omega} + \overline{K}\,\delta,
\end{equation}
then we have that
\begin{equation}
\label{Om sig subset X de}
\Omega(\sigma) \subset X(\delta) \subset \Omega.
\end{equation}
\end{lm}

\begin{proof}
We have that $X(\delta) \subset \Omega$ by construction.
\par
To show the first inclusion in \eqref{Om sig subset X de} we proceed by contradiction. Since the maximal cap $\Omega_\lambda$ contains a ball of radius $r_\Omega/4$, then $X(\delta)$ intersects $\Omega(\sigma)$. Assume that there exists a point $y \in \Omega(\sigma) \setminus X(\delta)$ and let $x$ be any point in $X(\delta) \cap \Omega(\sigma)$. Since $\Omega(\sigma)$ is connected, $x$ is joined to $y$ by a path contained in $\Omega(\sigma)$. Let $z$ be the first point on this path which falls outside $X(\delta)$. It is clear that $z \in \partial X(\delta) \cap \Omega(\sigma)$. We now consider two cases.
\par
If $z\cdot \omega < \lambda$ then the reflection $z^\lambda$ of $z$ in $\pi_\lambda$ is such that $z^\lambda \in \partial \Sigma_\delta$ and $\mbox{\rm dist}(z^\lambda,\partial \Omega) =\delta$. Since $u(z)=w(z^\lambda)+u(z^\lambda)$, from \eqref{w Linfty Sigma} and  \eqref{estimate-distance} we have that
\begin{equation} \label{34 by contrad}
\underline{K}\,\sigma\le u(z) \leq C \delta^{-\gamma} [u_\nu]_{\partial \Omega} + \overline{K}\,\delta,
\end{equation}
a contradiction.
\par
If $z\cdot \omega \geq \lambda$, then $z \in \partial \Sigma_\delta$ and $\mbox{\rm dist}(z,\partial \Omega) =\delta$. Hence, from \eqref{estimate-distance} we obtain that $u(z) \leq \overline{K}\,\delta$ and \eqref{34 by contrad} holds as well.
Since $z \in \Omega(\sigma)$ then from \eqref{estimate-distance} we have that $u(z) \geq \underline{K}\, \sigma$ which contradicts \eqref{K si geq} on account of \eqref{34 by contrad}.
\end{proof}

We draw the reader's attention on the differences in \eqref{K si geq} compared to \cite[Eq. (34)]{ABR}:
(i) thanks to Lemma \ref{lemma harnack cone}, the term $\delta^{-\gamma}$ replaces one in \cite[Eq. (34)]{ABR} that blows up exponentially; (ii) due to the different definition of the symmetrized set $X(\delta)$ (denoted by $X_\delta$ in \cite{ABR}), we
simplify the last summand in \eqref{K si geq}.

\begin{prop} 
\label{prop 5 ABR}
Let $\gamma$ be given by \eqref{gamma}. 
There exist positive numbers $C$ and $\varepsilon$, depending on $\Omega,d_\Omega, r_\Omega, \max_{\overline{\Omega}}u,f,\min_{\partial\Omega} u_\nu$, and $\sigma,\, \delta>0$ such that  \eqref{Om sig subset X de} holds with
\begin{equation} 
\label{de leq si}
\delta < \sigma < C\, [u_\nu]_{\partial \Omega}^\frac{1}{\gamma+1},
\end{equation}
provided that $[u_\nu]_{\partial \Omega}\le\varepsilon$.
\end{prop}

\begin{proof}
We must choose $\delta, \sigma\le r_\Omega/16$ that satisfy \eqref{K si geq}. 
If we let
\begin{equation*}
\delta = \left(\frac{C}{\overline{K}}\right)^{\frac{1}{\gamma+1}} [u_\nu]_{\partial \Omega}^{\frac{1}{\gamma+1}}
\ \mbox{ and } \ 
\sigma = \frac{4\overline{K}}{\underline{K}} \delta,
\end{equation*}
since $\overline{K}\,\delta=C\,\delta^{-\gamma}\,[u_n]_{\partial\Omega}$, then \eqref{K si geq} holds provided that
$\sigma\le r_\Omega/16$, that is for
\begin{equation*}
\varepsilon \leq \frac{\overline{K}}{C} \left( \frac{r_\Omega \underline{K}}{64\,\overline{K}} \right)^{\gamma +1},
\end{equation*}
The conclusion then follows from Lemma \ref{lemma 34 ABR}.
\end{proof}
\par
We are now ready to prove the main result of this paper.

\begin{proof}[Proof of Theorem \ref{stability serrin}]
The proof is analogous to that of \cite[Theorem 1]{ABR}; here, we use
\eqref{de leq si} in place of \cite[Eq. (33)]{ABR}. 
\par
Indeed, for any fixed direction $\omega\in{\mathbb S}^{N-1}$, 
Proposition \ref{prop 5 ABR} implies that, if $\sigma$ satisfies \eqref{de leq si},  then for any $x \in \partial \Omega$ there
exists $y_\omega\in \partial \Omega$ such that
\begin{equation} 
\label{x lam y}
|x^\lambda - y_\omega | \leq 2 \sigma,
\end{equation}
where $x^\lambda$ is the reflection of $x$ in the critical hyperplane $\pi_\lambda$ with $\lambda$ defined by 
\eqref{lambda def} (see also \cite[Corollary 2]{ABR}). 
\par
We now choose $N$
orthogonal directions, say $e_1,\ldots,e_N$, denote by
$\pi_1,\ldots,\pi_N$ the corresponding critical hyperplanes and we place the origin $0$ of ${\mathbb R^N}$
at the (unique) point in common to all the $\pi_j$'s. If we denote by
$R_j(x)$ the reflection of a point $x$ in $\pi_j$, we have that 
$$
-x= (\mathcal{R}_{N} \circ \cdots \circ \mathcal{R}_{1})(x);
$$
also, if we set $y_0=x$, $y_1=y_{e_1}$  and, for $j=2,\dots, N$,  $y_j$ as the point in $\partial\Omega$ such
that $|\mathcal{R}_j(y_{j-1})-y_j|\le 2 \sigma$ determined by \eqref{x lam y}, we obtain that
\begin{multline*}
|-x-y_N|\le \\
|(\mathcal{R}_{N} \circ \cdots \circ \mathcal{R}_{1})(x)-\mathcal{R}_{N}(y_{N-1})|+
|\mathcal{R}_{N}(y_{N-1})-y_N|\le \\
|(\mathcal{R}_{N-1} \circ \cdots \circ \mathcal{R}_{1})(x)-y_{N-1}|+2 \sigma\le \\
|(\mathcal{R}_{N-2} \circ \cdots \circ \mathcal{R}_{1})(x)-y_{N-2}|+4 \sigma\le \cdots \le 2N\sigma.
\end{multline*}
\par
Thus, we showed that there exists $y=y_N\in\partial\Omega$ such that
\begin{equation*} \label{Rx lam y}
|x +y | \leq 2 N \sigma.
\end{equation*}
This fact and \eqref{x lam y} imply that $0$ is an approximate
center of symmetry for $\Omega$, in the sense of \cite[Proposition 6]{ABR}, that is, for any direction $\omega$,  for the critical hyperplane $\pi$ in the direction $\omega$ we have that
\begin{equation} \label{bound prop 6 ABR}
\mbox{\rm dist}(0,\pi) \leq 4N\,[1+d_\Omega]\,\sigma.
\end{equation}
\par
By letting
\begin{equation*}
r_i=\min_{x \in \partial \Omega} |x|, \quad \textmd{and }\ \ r_e=\max_{x \in \partial
\Omega} |x|,
\end{equation*}
we obtain \eqref{stability-improved} thanks to \eqref{bound prop 6 ABR}, \cite[Proposition 7]{ABR} and \eqref{de leq si}.
\end{proof}
\par
Notice that Theorem \ref{stability serrin} was proved by fixing $a=1/2$ in Lemma \ref{lemma harnack cone}. However, an analog of Theorem \ref{stability serrin} for any fixed $a \in (0,1)$ can be proved by the same arguments and the number $\gamma$ in \eqref{gamma}  and the exponent in \eqref{stability-improved} may be optimized in terms of $a$. It is clear that the constants $\varepsilon$ and $C$ will change accordingly  (the smaller is $\gamma$, the larger is $C$ and smaller is $\varepsilon$). 
By using this plan, in case $\Omega$ is convex and $u$ is a solution of the torsional rigidity problem  
--- i.e. when $f(u)=1$ in \eqref{serrin1} --- we are able to give more explicit formulas
and partially compare our results to those in \cite{BNST2}. 




\begin{corollary}
Let $\Omega$ be a bounded convex domain with boundary of class $C^{2,\alpha}$. Let $u$ be a solution of 
$$
\Delta u = -1 \ \mbox{ in } \ \Omega, \quad u=0 \ \mbox{ on } \ \partial\Omega.
$$
\par
Then, for any $\eta>0$ there exist constants $C$ and $\varepsilon$, that
depend on $N, \Omega,d_\Omega, r_\Omega$ and $\eta$,
such that
\begin{equation} \label{stability convex torsion}
r_e-r_i \leq C\,[u_\nu]_{\partial \Omega}^{\frac{1}{\tau+\eta}},
\end{equation}
provided that $[u_\nu]_{\partial \Omega} \leq \varepsilon$, where
\begin{equation} \label{tau}
\tau = 1+N\,\sqrt{1+\left[\frac{2d_\Omega}{r_\Omega}\right]^2}.
\end{equation}
\end{corollary}

\begin{proof}
We recall that, in the case in hand, $\gamma$ is given by \eqref{gamma-torsion} and is increasing as $a$ 
grows; thus, the optimal exponent should be looked for when $a\to 0^+$.  Therefore, for
a $(\mathcal{C},\theta)$-domain, the exponent in \eqref{stability-improved} behaves as 
$$
\frac1{1+\gamma}= \frac{1}{1+N/\sin\theta+ o(1)} \ \mbox{ as } \ a\to 0^+.
$$ 
The conclusion then follows from Proposition \ref{proposition cone}.
\par
We finally notice that, in this case, the dependence of the constants $C$ and $\varepsilon$ appearing
in Theorem \ref{stability serrin} on the quantities $\|u\|_\infty$ and $\min_{\partial\Omega} u_\nu)$ can be removed, since these quantities can be bounded in terms of the $C^2$ regularity of $\Omega$, by standard barrier techniques.
\end{proof}
\begin{remark}
In \cite[Theorem 2]{BNST2} the authors prove an estimate similar to \eqref{stability convex torsion}. When $\Omega$ is convex, our estimate improves that in \cite[Theorem 2]{BNST2} if 
$$
\frac{d_\Omega}{r_\Omega} \leq \sqrt{\frac{2N^2+N-5/2}{N}}.
$$
\end{remark}

\section*{Acknowledgements} 
The authors wish to thank Paolo Salani for the useful discussions we had together.

The authors have been supported by the Gruppo Nazionale per l'Analisi Matematica, la Probabilit\`a e le loro
Applicazioni (GNAMPA) of the Istituto Nazionale di Alta Matematica (INdAM). 

The paper was completed while the author was visiting \lq\lq The Institute for Computational Engineering and Sciences\rq\rq (ICES) of The University of Texas at Austin, and he wishes to thank the Institute for hospitality and support. The author has been also supported by the NSF-DMS Grant 1361122 and the Firb project 2013 ``Geometrical and Qualitative aspects of PDE''.



2015-06-22T02:06:55.000
1410
1410.7791
en
https://arxiv.org/abs/1410.7791
In a bounded domain $\Omega$, we consider a positive solution of the problem $\Delta u+f(u)=0$ in $\Omega$, $u=0$ on $\partial\Omega$, where $f:\mathbb{R}\to\mathbb{R}$ is a locally Lipschitz continuous function. Under sufficient conditions on $\Omega$ (for instance, if $\Omega$ is convex), we show that $\partial\Omega$ is contained in a spherical annulus of radii $r_i<r_e$, where $r_e-r_i\leq C\,[u_\nu]_{\partial\Omega}^\alpha$ for some constants $C>0$ and $\alpha\in (0,1]$. Here, $[u_\nu]_{\partial\Omega}$ is the Lipschitz seminorm on $\partial\Omega$ of the normal derivative of $u$. This result improves to Hölder stability the logarithmic estimate obtained in [1] for Serrin's overdetermined problem. It also extends to a large class of semilinear equations the Hölder estimate obtained in [6] for the case of torsional rigidity ($f\equiv 1$) by means of integral identities. The proof hinges on ideas contained in [1] and uses Carleson-type estimates and improved Harnack inequalities in cones.
Analysis of PDEs (math.AP)
Hölder stability for Serrin's overdetermined problem
Qwen/Qwen-72B
1. YES
2. YES