\section*{\refname
       \@mkboth{\MakeUppercase\refname}{\MakeUppercase\refname}}%
     \list{\@biblabel{\@arabic\c@enumiv}}%
          {\settowidth\labelwidth{\@biblabel{#1}}%
           \leftmargin\labelwidth
           \advance\leftmargin\labelsep
           \itemsep\bibspace
           \parsep\z@skip     %
           \@openbib@code
           \usecounter{enumiv}%
           \let\p@enumiv\@empty
           \renewcommand\theenumiv{\@arabic\c@enumiv}}%
     \sloppy\clubpenalty4000\widowpenalty4000%
     \sfcode`\.\@m}
    {\def\@noitemerr
      {\@latex@warning{Empty `thebibliography' environment}}%
     \endlist}

\makeatother

\makeatletter

\renewcommand{\rmdefault}{ptm}

\newtheorem{thm}{Theorem}[section]
\newtheorem{lem}[thm]{Lemma}
\newtheorem{prop}[thm]{Proposition}
\newtheorem{defn}[thm]{Definition}
\newtheorem{ex}[thm]{Example}
\newtheorem{cor}[thm]{Corollary}
\newtheorem{rem}[thm]{Remark}
\newtheorem{rems}[thm]{Remarks}
\newtheorem{question}[thm]{Question}

\def\Xint#1{\mathchoice
  {\XXint\displaystyle\textstyle{#1}}%
  {\XXint\textstyle\scriptstyle{#1}}%
  {\XXint\scriptstyle\scriptscriptstyle{#1}}%
  {\XXint\scriptscriptstyle\scriptscriptstyle{#1}}%
  \!\int}
\def\XXint#1#2#3{{\setbox0=\hbox{$#1{#2#3}{\int}$}
  \vcenter{\hbox{$#2#3$}}\kern-.5\wd0}}
\def\Xint={\Xint=}
\def\Xint-{\Xint-}

\newcommand{\alpha}                \newcommand{\lda}{\lambda}{\alpha}                \newcommand{\lda}{\lambda}
\newcommand{\Omega}                \newcommand{\pa}{\partial}{\Omega}                \newcommand{\pa}{\partial}
\newcommand{\varepsilon}           \newcommand{\ud}{\mathrm{d}}{\varepsilon}           \newcommand{\ud}{\mathrm{d}}
\newcommand{\begin{equation}}      \newcommand{\ee}{\end{equation}}{\begin{equation}}      \newcommand{\ee}{\end{equation}}
\newcommand{\w}{\omega}
\newcommand{\Lambda}              \newcommand{\B}{\mathcal{B}}{\Lambda}              \newcommand{\B}{\mathcal{B}}
\newcommand{\mathbb{R}}              \newcommand{\Sn}{\mathbb{S}^n}{\mathbb{R}}              \newcommand{\Sn}{\mathbb{S}^n}
\newcommand{\mathcal{D}^{\sigma,2}}  \newcommand{\M}{\mathscr{M}}{\mathcal{D}^{\sigma,2}}  \newcommand{\M}{\mathscr{M}}
\newcommand{\displaystyle\int}{\displaystyle\int}
\newcommand{\displaystyle\lim}{\displaystyle\lim}
\newcommand{\displaystyle\sup}{\displaystyle\sup}
\newcommand{\displaystyle\min}{\displaystyle\min}
\newcommand{\displaystyle\max}{\displaystyle\max}
\newcommand{\displaystyle\inf}{\displaystyle\inf}
\newcommand{\displaystyle\sum}{\displaystyle\sum}
\newcommand{\abs}[1]{\lvert#1\rvert}

\begin{document}

\title{\textbf{Bubbling and extinction for some fast diffusion equations in bounded domains}
\bigskip}

\author{\medskip  Tianling Jin\footnote{T. Jin is partially supported by Hong Kong RGC grants GRF 16302217 and GRF 16306320.}, \  \
Jingang Xiong\footnote{J. Xiong is partially supported by NSFC 11922104 and 11631002.}}

\date{August 4, 2020}

\maketitle

\begin{abstract}
We study a Sobolev critical fast diffusion equation in bounded domains with the Brezis-Nirenberg effect. We obtain extinction profiles of its positive solutions, and show that the convergence rates of the relative error in regular norms are at least polynomial. Exponential decay rates are proved for generic domains. Our proof makes use of its regularity estimates, a curvature type evolution equation, as well as blow up analysis. Results for Sobolev subcritical fast diffusion equations are also obtained.

\end{abstract}

\section{Introduction}

Let $\Omega}                \newcommand{\pa}{\partial$ be a bounded domain in $\mathbb{R}}              \newcommand{\Sn}{\mathbb{S}^n^n$, $n\ge 3$, with smooth boundary $\pa \Omega}                \newcommand{\pa}{\partial$. We consider the Cauchy-Dirichlet  problem for the  fast diffusion equation with the  Sobolev critical exponent
\begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:BN-1}
\begin{split}
\frac{\pa }{\pa t}u^{\frac{n+2}{n-2}} &= \Delta u+ b u\quad  \mbox{in }\Omega}                \newcommand{\pa}{\partial\times (0,\infty),\\
u&=0 \quad\quad\quad\quad \mbox{on }\pa\Omega}                \newcommand{\pa}{\partial \times (0,\infty),\\
u(\cdot,0)&=u_0\ge 0\quad\ \ \   \mbox{in }\Omega}                \newcommand{\pa}{\partial,
\end{split}
\ee
where $\Delta=\sum_{i=1}^n \frac{\pa^2 }{\pa x_i^2}$ is the Laplace operator, $u_0$ is not identically zero, and
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:b}
b\in[0,\lda_1)\mbox{ is a constant}
\ee
with $\lda_1$ being the first eigenvalue of $-\Delta$ in $\Omega$ with zero Dirichlet boundary condition. Hence, the operator $-\Delta -b$ is coercive on the Sobolev space $H_0^1(\Omega}                \newcommand{\pa}{\partial)$. The fast diffusion equations arise in the modelling of gas-kinetics, plasmas, thin liquid film dynamics driven by Van der Waals forces, and etc. If $b=0$, this Sobolev critical equation \eqref{eq:BN-1} can be viewed a unnormalized Yamabe flow with metrics degenerate on the boundary.

The theory of existence and uniqueness of solutions to \eqref{eq:BN-1} is well understood, see V\'azquez \cite{Vaz,Vaz1}. If $u_0\in L^q(\Omega}                \newcommand{\pa}{\partial)$ for some $q>\frac{2n}{n-2}$, then the solution will become instantaneously positive in $\Omega$ and globally bounded. Moreover, the solution  will vanish in a finite time $T^*>0$. If we assume that $u_0\in H^1_0(\Omega}                \newcommand{\pa}{\partial)\cap L^q(\Omega}                \newcommand{\pa}{\partial)$ for some $q>\frac{2n}{n-2}$, then it follows from the work of Chen-DiBenedetto \cite{CDi}, DiBenedetto-Kwong-Vespri  \cite{DKV} and Jin-Xiong \cite{JX19} that the solution is of $C^{3,2}_{x,t}(\overline\Omega}                \newcommand{\pa}{\partial\times(0,T^*))$. In particular, the solutions are classical. Therefore, when we investigate the asymptotic behavior of nonnegative  solutions  to \eqref{eq:BN-1} as $t$ approaching to the extinction time $T^*$, there is no loss of generality to  consider classical (up to the boundary) solutions  to \eqref{eq:BN-1}.

When $\frac{n+2}{n-2}$ is replaced by $p\in (1,\frac{n+2}{n-2})$ if $n\ge 3$, or $p\in(1,\infty)$ if $n=1,2$, which is a Sobolev subcritical exponent, the extinction behavior of solutions to the fast diffusion equation
\begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:subcritical}
\begin{split}
\frac{\pa }{\pa t}u^{p} &= \Delta u\quad  \mbox{in }\Omega}                \newcommand{\pa}{\partial\times (0,\infty),\\
u&=0 \quad\quad \mbox{on }\pa\Omega}                \newcommand{\pa}{\partial \times (0,\infty)
\end{split}
\ee
has been well-studied.  By the scaling
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:subchanging_variables}
v(x,t)=\Big (\frac{p}{(p-1)(T^*-\tau)}\Big)^{\frac{1}{p-1}} u(x,\tau), \quad t =\frac{p}{p-1}   \ln \Big(\frac{T^*}{T^*-\tau}\Big),
\ee
where $T^*$ is the extinction time, the equation \eqref{eq:subcritical} becomes
 \begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:subcritical2}
 \begin{split}
\frac{\pa }{\pa t}v^{p} &= \Delta v+v^{p} \quad \mbox{in }\Omega}                \newcommand{\pa}{\partial\times (0,\infty),\\
 v&=0 \quad \mbox{on }\pa\Omega}                \newcommand{\pa}{\partial \times (0,\infty).
\end{split}
\ee
Berryman-Holland \cite{BH} proved that the solution of \eqref{eq:subcritical2} converges to a stationary solution $v_\infty$ in $H^1_0(\Omega)$ along a sequence of times. Feireisl-Simondon \cite{FS} proved the full convergence in the $C^0(\overline\Omega)$ topology. Bonforte-Grillo-V\'azquez \cite{BGV}  proved that the relative error $v(\cdot,t)/v_\infty$ converges to $1$ in $L^\infty(\Omega)$.
Recently, Bonforte-Figalli \cite{BFig} proved the exponential convergence of the relative error for generic domains $\Omega$, which means that the domains $\Omega$ satisfy
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{condition2}
\begin{split}
&\mbox{For every nonnegative } H^1_0\ \mbox{ solution } v \mbox{ of }  -\Delta v- v^{p}=0 \mbox{ in } \Omega, \mbox{the linearized} \\[-1mm]
&\mbox{operator at }v,
\mbox{that is } L_v:=-\Delta - p v^{p-1}, \mbox{ has a trivial kernel in } H^1_0(\Omega).
\end{split}
\ee
The set of  smooth domains satisfying \eqref{condition2}  has generic properties, see Saut-Temam \cite{ST}.

The main advantage of the subcritical regime is the upper bound of solutions $u$ to \eqref{eq:subcritical} proved in DiBenedetto-Kwong-Vespri  \cite{DKV}
\begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:idealbound}
u(x,t) \le C d(x)(T^*-t)^{\frac{1}{p-1}} \quad \mbox{for }t<T^*,
\ee
where $d(x)=dist(x,\pa \Omega}                \newcommand{\pa}{\partial)$. The estimate \eqref{eq:idealbound} implies that the function $v$ defined by \eqref{eq:subchanging_variables}, which satisfies \eqref{eq:subcritical2},  is uniformly bounded as $t\to\infty$, and consequently, has uniform regularity estimates up to the boundary $\partial\Omega$ by the work of \cite{CDi,DKV,JX19}.

However, this uniform bound in general does not hold for \eqref{eq:subcritical2} if $p= \frac{n+2}{n-2}$. For instance, it is the case if $\Omega}                \newcommand{\pa}{\partial$ is star-shaped, since there is no stationary solution of \eqref{eq:subcritical2} due to the Pohozaev identity. In this paper, we will show that the uniform boundedness still holds for the equation \eqref{eq:BN-1} assuming $b>0$ and $n\ge 4$. The role of the positivity of $b$ when $n\ge 4$ was first discovered  in the seminal paper Brezis-Nirenberg \cite{BN}, and is similar to the role that the non-vanishing Weyl tensor and the positive mass theorem play in the resolution of the Yamabe problem on  compact manifolds by Aubin \cite{Aubin} and Schoen \cite{Schoen}.

Under the scaling
 \begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:changing_variables}
 v(x,t)=\left(\frac{n+2}{4(T^*-\tau)}\right)^{\frac{n-2}{4}} u(x,\tau ), \quad t=\frac{n+2}{4} \ln \left(\frac{T^*}{T^*-\tau}\right),
 \ee
the equation \eqref{eq:BN-1} becomes
 \begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:BN-3}
 \begin{split}
\frac{\pa }{\pa t}v^{\frac{n+2}{n-2}} &= \Delta v+ b v +v^{\frac{n+2}{n-2}} \quad \mbox{in }\Omega}                \newcommand{\pa}{\partial\times (0,\infty),\\
 v&=0 \quad \mbox{on }\pa\Omega}                \newcommand{\pa}{\partial \times (0,\infty).
\end{split}
\ee
We will show that every solution of \eqref{eq:BN-3}  converges to a stationary solution, that is a solution of
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:BNstationary}
\Delta v+ b v +v^{\frac{n+2}{n-2}}=0 \quad \mbox{in }\Omega}                \newcommand{\pa}{\partial,\quad v=0\quad \mbox{on }\partial\Omega}                \newcommand{\pa}{\partial,
\ee
with at least polynomial rates. Moreover, the convergence rate will be  exponential if the domain $\Omega$ satisfies the following condition:
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{condition}
\begin{split}
&\mbox{For every nonnegative } H^1_0\ \mbox{ solution } v \mbox{ of }  -\Delta v- b v -v^{\frac{n+2}{n-2}}=0 \mbox{ in } \Omega, \mbox{the linearized } \\[-1mm]
&\mbox{operator at }v, \mbox{that is } L_v:=-\Delta - b  -\frac{n+2}{n-2}v^{\frac{4}{n-2}}, \mbox{ has a trivial kernel in } H^1_0(\Omega).
\end{split}
\ee
The set of  smooth domains satisfying \eqref{condition} also has generic properties, see Saut-Temam \cite{ST}.

\begin{thm}\label{thm:critical}
Let $n\ge 4$, and $b>0$ satisfy \eqref{eq:b}. Let $u$ be a classical nonnegative solution of \eqref{eq:BN-1} with extinction time $T^*>0$. Let $v$ be defined by \eqref{eq:changing_variables}. Then there is a nonzero stationary solution $v_\infty$ of \eqref{eq:BN-3}, and two positive constants $\theta$ and $C$ such that
\[
\left\|\frac{v(\cdot,t)}{v_\infty}-1\right\|_{C^2(\overline\Omega}                \newcommand{\pa}{\partial)}\le Ct^{-\theta}\quad\mbox{for all }t\ge 1.
\]
If $\Omega$ satisfies \eqref{condition}, then there exist  two positive constants $\gamma$ and $C$ such that
\[
\left\|\frac{v(\cdot,t)}{v_\infty}-1\right\|_{C^2(\overline\Omega}                \newcommand{\pa}{\partial)}\le Ce^{-\gamma t}\quad\mbox{for all }t\ge 1.
\]

\end{thm}

When $n=3$, it was shown in Brezis-Nirenberg \cite{BN} that the situation for the stationary equation \eqref{eq:BNstationary} changes drastically from dimensions $n\ge 4$. The positivity of $b$ is not sufficient to give a minimal energy solution of \eqref{eq:BNstationary}. Druet \cite{Druet} showed that the necessary and sufficient condition is the positivity of the regular part of the Green's function of $-\Delta-b$ at a diagonal point. There should be similar changes for the parabolic equation \eqref{eq:BN-3} as well.

When $b=0$,  Sire-Wei-Zheng \cite{SWZ}  recently proved the existence of some initial data such that the solution of \eqref{eq:BN-3} blows up at finitely many points with an explicit blow up rate as  $t\to\infty$, using the gluing method for parabolic equations in the spirit of Cort\'azar-del Pino-Musso \cite{CdM} and D\'avila-del Pino-Wei \cite{DdW}.  This generalizes and provides rigorous proof of a result of Galaktionov-King \cite{GKing} for the radially symmetric case, where the solution blows up at one point. A class of type II ancient solutions to the Yamabe flow, which are rotationally
symmetric and converge to a tower of spheres as $t\to-\infty$, was constructed by  Daskalopoulos-del Pino-Sesum \cite{DdS}.  Bubble tower solutions for the energy critical heat equation were constructed in  del Pino-Monica-Wei \cite{dMW}.  It is conjectured in Sire-Wei-Zheng \cite{SWZ} that bubble tower solutions to \eqref{eq:BN-3} with $b=0$ also exist. Nevertheless, if it is the global case ($\Omega$ replaced by $\mathbb{R}}              \newcommand{\Sn}{\mathbb{S}^n^n$), then it has been proved by del Pino-S\'aez \cite{delSaez} that the solution of \eqref{eq:BN-3} for $b=0$ with fast decay initial data will converge to a nontrivial stationary solution, which is in fact a standard bubble.

To prove Theorem \ref{thm:critical},  we will adapt the blow up analysis of Struwe \cite{St}, Bahri-Coron \cite{Bahri-C}, Schwetlick-Struwe \cite{SS} and Brendle \cite{Br05}.  See also Chen-Xu \cite{CX} and Mayer \cite{Mayer} for similar analysis of scalar curvature flows. Here, we define a curvature type quantity $\mathcal{R}$, and derive its equation along the parabolic equation \eqref{eq:BN-3}. Due to the lack of information of  $\mathcal{R}$ on the boundary $\partial\Omega$, extra work is needed to obtain estimates for $\mathcal{R}$. Here the optimal boundary regularity proved in our previous paper \cite{JX19} is crucial.   Part  of the blow up  analysis in this paper remains valid  when $b=0$ or $n=3$.  The condition  $n\ge 4$ and $b>0$ is used in the final step to rule out bubbles.

Our proof of the polynomial decay rates in Theorem \ref{thm:critical} can be applied to prove the polynomial rate of the convergence of the relative error for the Sobolev subcritical fast diffusion equation \eqref{eq:subcritical2} in all smooth domains. We also provide an alternative proof the exponential convergence result of Bonforte-Figalli \cite{BFig} for $\Omega$ satisfying \eqref{condition2}.
\begin{thm}\label{thm:subcritical}
Suppose $p\in (1,\frac{n+2}{n-2})$ if $n\ge 3$, and $p\in(1,\infty)$ if $n=1,2$. Let $u$ be a classical nonnegative solution of \eqref{eq:subcritical} with extinction time $T^*>0$. Let $v$ be defined by \eqref{eq:subchanging_variables}. Then there is a stationary solution $v_\infty$ of \eqref{eq:subcritical2}, and two positive constants $\theta$ and $C$ such that
\[
\left\|\frac{v(\cdot,t)}{v_\infty}-1\right\|_{C^2(\overline\Omega}                \newcommand{\pa}{\partial)}\le Ct^{-\theta}\quad\mbox{for all }t\ge 1.
\]
If $\Omega$ satisfies \eqref{condition2}, then there exist  two positive constants $\gamma$ and $C$ such that
\[
\left\|\frac{v(\cdot,t)}{v_\infty}-1\right\|_{C^2(\overline\Omega}                \newcommand{\pa}{\partial)}\le Ce^{-\gamma t}\quad\mbox{for all }t\ge 1.
\]
\end{thm}
 We know from the work of Carlotto-Chodosh-Rubinstein \cite{CCR}  that  there exists a Yamabe flow on $\mathbb{S}^1(1/\sqrt{n-2})\times\mathbb{S}^{n-1}(1)$ such that it converges exactly at a polynomial rate. In a forthcoming paper \cite{JXslow}, we will show examples of polynomial rate convergence of the relative error for the solutions of the fast diffusion equations \eqref{eq:subcritical2} and \eqref{eq:BN-3}.

This paper is organized as follows. Sections \ref{sec:integralbound}--\ref{sec:convergence} deal with the critical equation \eqref{eq:BN-3}. We first obtain certain integral bounds for solutions of this critical equation in section \ref{sec:integralbound}. Sections \ref{sec:concentration} is for the possible concentration phenomenon for its solutions. In section \ref{sec:blowup}, we use blow up analysis to rule out such possible concentration phenomenon. Section \ref{sec:convergence} is devoted to the proof of the uniform boundedness and convergence results in Theorem \ref{thm:critical}. In section \ref{sec:subcritical}, we consider  the subcritical equation \eqref{eq:subcritical2} and prove Theorem \ref{thm:subcritical}.

\bigskip

\noindent \textbf{Acknowledgement:} Part of this work was completed while the second named author was visiting the Hong Kong University of Science and Technology and Rutgers University, to which he is grateful for providing  very stimulating research environments and supports. Both authors would like to thank Professors Ha\"im Brezis and YanYan Li for their interests and useful comments.

\section{Integral bounds}\label{sec:integralbound}

 For an open set $\Omega}                \newcommand{\pa}{\partial$, let $H_0^1(\Omega}                \newcommand{\pa}{\partial)$ be the closure of $C_c^\infty(\Omega}                \newcommand{\pa}{\partial)$ under the norm
\[
\|u\|_{H_0^1(\Omega}                \newcommand{\pa}{\partial)}:=\left(\int_\Omega}                \newcommand{\pa}{\partial |\nabla u|^2 \,\ud x\right)^{1/2}.
\]
For convenience, we define
\begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:norm}
\|u\|:=\left(\int_\Omega}                \newcommand{\pa}{\partial |\nabla u|^2-bu^2 \,\ud x\right)^{1/2}
\ee
and
\begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:inner-product}
\langle u,v\rangle= \int_{\Omega}                \newcommand{\pa}{\partial} (\nabla u\nabla v-buv)\,\ud x
\ee
be the associated inner product.
Since $b<\lda_1$, there exists a constant $K_b>0$ such that
\begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:sobolev-b}
\|u\|_{L^{\frac{2n}{n-2}}} \le K_b^{1/2} \|u\| \quad \mbox{for any }u\in H_0^1(\Omega}                \newcommand{\pa}{\partial).
\ee

Recall that from \cite{JX19}, we know that the solution $u(x,t)$ of \eqref{eq:BN-1} is smooth in $t\in (0,T^*)$ for every $x\in\overline\Omega$, and $ \pa_t^l u(\cdot,t) \in C^{\frac{3n-2}{n-2}}(\overline \Omega}                \newcommand{\pa}{\partial)$ for all $l\ge 0$ and all $t\in(0,T^*)$.

\begin{lem}\label{lem:basic} Let $u$ be a solution of \eqref{eq:BN-1}, and $T^*$ be the extinction time of $u$. Then for every $0<t<T^*$,
 \[
\frac{1}{C}(T^*-t)^{\frac{n}{2}} \le  \int_\Omega}                \newcommand{\pa}{\partial u(x,t)^{\frac{2n}{n-2}}\,\ud x \le (1-t/T^*)^{\frac{n}{2}} \int_\Omega}                \newcommand{\pa}{\partial u_0(x)^{\frac{2n}{n-2}}\,\ud x,
 \]
 where $C$ is a positive constant depending only on $n$ and $\Omega}                \newcommand{\pa}{\partial$.

 \end{lem}

 \begin{proof} If $b=0$, the lemma was proved by \cite{BH} (noting that by our regularity result in \cite{JX19}, the regularity assumptions in \cite{BH} are satisfied).   The same proof applies if $b\in[0,\lambda_1)$.
 \end{proof}

 Let $v$ be as in \eqref{eq:changing_variables}. By Lemma \ref{lem:basic}, we have
\begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:BN-5}
\frac{1}{C} \le \int_{\Omega}                \newcommand{\pa}{\partial} v(x,t)^{\frac{2n}{n-2}}\,\ud x\le C
\ee
for some positive constant $C$ depending only on $n,\Omega$ and $u_0$.  Define
\begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:F-funct}
F(v(t))= \int_{\Omega}                \newcommand{\pa}{\partial } \Big( |\nabla v(x,t)|^2 -b v(x,t)^2 -\frac{n-2}{n} v(x,t)^{\frac{2n}{n-2}}\Big)\,\ud x.
\ee
It follows that $F(v(t))$ is bounded from below. By the equation of $v$ and integrating by parts,
\begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:F-dcr}
\frac{\ud}{\ud t} F(v(t))= - 2\int_{\Omega}                \newcommand{\pa}{\partial}(\Delta v +bv+v^{\frac{n+2}{n-2}}) \pa_t v \,\ud x= -\frac{2(n+2)}{n-2} \int_{\Omega}                \newcommand{\pa}{\partial }v^{\frac{4}{n-2}} |\pa_t v|^2 \,\ud x \le 0.
\ee
Hence, $F(v(t))$ is non-increasing in $t$, and thus, together with \eqref{eq:BN-5}, we have $\|v(\cdot,t)\|_{H^1_0(\Omega)}$ is uniformly bounded. Moreover, there exists some constant $F_\infty$ such that
\begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:f_infty}
\lim_{t\to \infty} F(v(t)) = F_\infty.
\ee

Define
\begin{equation}\label{eq:Q}
\mathcal{R}= v^{-\frac{n+2}{n-2}}(-\Delta v-bv)
\end{equation}
and
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:definitionMq}
M_q(t)= \int_{\Omega}                \newcommand{\pa}{\partial} |\mathcal{R} -1|^q v^{\frac{2n}{n-2}} \,\ud x, \quad 	\ge 1.
\ee
In \cite{JX19}, we proved that $\mathcal{R}=1-\frac{n+2}{n-2} \frac{\pa_t v}{v}$ is $C^{2}$ up to the boundary $\partial\Omega$. However, all the estimates there for solutions of \eqref{eq:BN-3} are only locally uniform in $t\in (0,\infty)$. We shall prove some uniform estimates for all $t\in [1,\infty)$ and $M_q(t) \to \infty$ as $t\to \infty$.

To do this, we will first use Moser's iteration to obtain a uniform lower bound of $\mathcal{R}$ as an intermediate step. So we need the following evolution equation of $\mathcal{R}$ and integration by parts formula.
\begin{lem} \label{lem:flow-properties} Let $g=u^{\frac{4}{n-2}} g_{flat}$. Then
\begin{itemize}
\item[(i).]
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:volumn}
\pa_t v^{\frac{2n}{n-2}}= -\frac{2n}{n+2}(\mathcal{R}-1) v^{\frac{2n}{n-2}}.
\ee
\item[(ii).]
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:R-1}
\pa_t (\mathcal{R}-1)=\frac{n-2}{n+2} \Delta_g (\mathcal{R}-1) +\frac{4}{n+2} (\mathcal{R}-1)^2 + \frac{4}{n+2} (\mathcal{R}-1),
\ee
where $\Delta_g $ is the Laplace-Beltrami operator of $g$.
\item[(iii).] For any $f\in H^2(\Omega}                \newcommand{\pa}{\partial)$ and $h\in  H^1(\Omega}                \newcommand{\pa}{\partial)$,
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:integrationbyparts}
\int_{\Omega}                \newcommand{\pa}{\partial}h \Delta_g f\,\ud vol_g = -\int_{\Omega}                \newcommand{\pa}{\partial} \langle \nabla_g f, \nabla_g h\rangle_g \,\ud vol_g.
\ee
\end{itemize}

\end{lem}

\begin{proof} The equation \eqref{eq:volumn} follows immediately from \eqref{eq:BN-3} and \eqref{eq:Q}. We also have $\pa_t v= \frac{n-2}{n+2} v (1-\mathcal{R})$.

By the definition of $\mathcal{R}$, we have
\begin{align*}
\pa_t (\mathcal{R}-1)&= \frac{n+2}{n-2} v^{-\frac{2n}{n-2}} \pa_t v (\Delta+ b) v- v^{-\frac{n+2}{n-2}}  (\Delta +b) \pa_t v\\&
= v^{-\frac{n+2}{n-2}} (1-\mathcal{R}) (\Delta+ b) v  - \frac{n-2}{n+2} v^{-\frac{n+2}{n-2}}  (\Delta +b) (v(1-\mathcal{R}) )\\&
=(\mathcal{R}-1)\mathcal{R} - \frac{n-2}{n+2} v^{-\frac{n+2}{n-2}}  (\Delta +b) (v(1-\mathcal{R}) ).
\end{align*}
Let $L_g= \Delta _g -\frac{n-2}{4(n-1)}R_g$ be the conformal Laplacian of $g$, where $R_g$ is the the scalar curvature of $g$. By the conformal transformation law
\[
L_g (v^{-1}\varphi)=v^{-\frac{n+2}{n-2}}\Delta\varphi,\quad\forall \varphi\in C^2(\Omega),
\]
we have
\[
\frac{n-2}{4(n-1)}R_g=-L_g(1)=-v^{-\frac{n+2}{n-2}}\Delta v=\mathcal{R} +b v^{-\frac{4}{n-2}}
\]
and
\begin{align*}
v^{-\frac{n+2}{n-2}}  (\Delta +b) (v(1-\mathcal{R}) )&=  L_g (1-\mathcal{R}) + bv^{-\frac{4}{n-2}} (1-\mathcal{R})\\&
= \Delta _g (1-\mathcal{R}) -\frac{n-2}{4(n-1)} R_g (1-\mathcal{R}) + bv^{-\frac{4}{n-2}} (1-\mathcal{R}) \\&
=  \Delta _g (1-\mathcal{R}) -\mathcal{R} (1-\mathcal{R}).
\end{align*}
Then, \eqref{eq:R-1} follows.

Finally,
\begin{align*}
\int_{\Omega}                \newcommand{\pa}{\partial}h \Delta_g f\,\ud vol_g &= \int_\Omega h  v^{-\frac{2n}{n-2}} \pa_i (v^{\frac{2n}{n-2}}  v^{-\frac{4}{n-2}} \pa_i  f) v^{\frac{2n}{n-2}}\,\ud x\\&
=  \int_\Omega h \pa_i (v^2 \pa_if) \,\ud x= - \int_\Omega v^2 \pa_i f \pa_i h\,\ud x= -\int_\Omega \langle \nabla_g f, \nabla_g h\rangle_g \,\ud vol_g,
\end{align*}
where we used $v=0$ on $\partial\Omega$ in the third equality.
\end{proof}

We have the following Sobolev inequality regarding the metric $g=v^{\frac{4}{n-2}} g_{flat}$:
\begin{lem} \label{lem:sob-elliptic}
There holds
\[
\left(\int_{\Omega}                \newcommand{\pa}{\partial} |f|^{\frac{2n}{n-2}}\,\ud vol_g \right)^{\frac{n-2}{n}}\le K_b\int_{\Omega}                \newcommand{\pa}{\partial}( |\nabla_g f|_g^2+\mathcal{R}f^2) \,\ud vol_g
\]
for any $f\in H^1(\overline \Omega}                \newcommand{\pa}{\partial)$, where $K_b$ is the constant in \eqref{eq:sobolev-b}.
\end{lem}
\begin{proof}
Note that
\[
|\nabla (fv)|= v^2 |\nabla f|^2 +f^2|\nabla v|^2 +2 vf\nabla v \cdot \nabla f,
\]
\begin{align*}
\int_\Omega (f^2|\nabla v|^2 +2 v\nabla v f \nabla f)\,\ud x&=  \int_\Omega (f^2|\nabla v|^2 +v\nabla v  \nabla f^2 ) \,\ud x\\&=
-\int_\Omega vf^2\Delta v \,\ud x\\&
= \int_\Omega (\mathcal{R} f^2 v^{\frac{2n}{n-2}}+ b v^2 f^2) \,\ud x.
\end{align*}
Hence,
\[
\int_{\Omega}                \newcommand{\pa}{\partial}( |\nabla_g f|_g^2+\mathcal{R}f^2) \,\ud vol_g=\int_{\Omega}                \newcommand{\pa}{\partial}( v^2|\nabla f|+\mathcal{R}f^2v^{\frac{2n}{n-2}}) \,\ud x = \int_\Omega}                \newcommand{\pa}{\partial (|\nabla (fv)|- b(fv)^2 )\,\ud x.
\]
Therefore, the lemma follows from \eqref{eq:sobolev-b}.
\end{proof}

For any $t_0\ge 0$ and $T> 0$, let
\[
V^1(\Omega}                \newcommand{\pa}{\partial \times (t_0,t_0+T)) = C^0((t_0,t_0+T); L^2(\Omega}                \newcommand{\pa}{\partial)) \cap L^2 ((t_0,t_0+T); H^1(\Omega}                \newcommand{\pa}{\partial)),
\]
equipped with the norm
\[
\|f\|_{V^1(\Omega}                \newcommand{\pa}{\partial \times (t_0,t_0+T))}^2= \sup_{t_0<t<t_0+T} \int_{\Omega}                \newcommand{\pa}{\partial } f(x,t)^2 \,\ud vol_g +\int_{t_0}^{t_0+T} \int_{\Omega}                \newcommand{\pa}{\partial} ( |\nabla_g f|_g^2+\mathcal{R}f^2) \,\ud vol_g \ud t.
\]

We have the following parabolic version of Sobolev inequality.
\begin{lem}\label{lem:sob-parabolic} For any $f\in V^1 (\Omega}                \newcommand{\pa}{\partial \times (t_0,t_0+T))$, we have
\[
\left(\int_{t_0}^{t_0+T} \int_{\Omega}                \newcommand{\pa}{\partial} |f|^{\frac{2(n+2)}{n}}\,\ud vol_g \ud t\right)^{\frac{n}{n+2}}\le K_b^{\frac{n}{n+2}} \|f\|_{V^1(\Omega}                \newcommand{\pa}{\partial \times (t_0,t_0+T))}^2.
\]
\end{lem}

\begin{proof} By H\"older's inequality and Lemma \ref{lem:sob-elliptic}, we have
\begin{align*}
 \int_{\Omega}                \newcommand{\pa}{\partial} |f|^{\frac{2(n+2)}{n}}\,\ud vol_g & =  \int_{\Omega}                \newcommand{\pa}{\partial} |f|^2|f|^{\frac{4}{n}}\,\ud vol_g \\&
 \le \Big( \int_{\Omega}                \newcommand{\pa}{\partial} |f|^{\frac{2n}{n-2}}\,\ud vol_g\Big )^{\frac{n-2}{n}} \Big( \int_{\Omega}                \newcommand{\pa}{\partial} |f|^{2}\,\ud vol_g\Big)^{\frac{2}{n}}\\&
 \le K_b\int_{\Omega}                \newcommand{\pa}{\partial}( |\nabla_g f|_g^2+\mathcal{R}f^2) \,\ud vol_g \Big( \int_{\Omega}                \newcommand{\pa}{\partial} |f|^{2}\,\ud vol_g\Big)^{\frac{2}{n}}.
\end{align*}
Hence, by  Young's inequality
\begin{align*}
&\left(\int_{t_0}^{t_0+T} \int_{\Omega}                \newcommand{\pa}{\partial} |f|^{\frac{2(n+2)}{n}}\,\ud vol_g \ud t\right)^{\frac{n}{n+2}} \\&\le K_b^{\frac{n}{n+2}} \left(\int_{\Omega}                \newcommand{\pa}{\partial}( |\nabla_g f|_g^2+\mathcal{R}f^2) \,\ud vol_g \right)^{\frac{n}{n+2}} \left( \sup_{t_0<t<t_0+T} \int_{\Omega}                \newcommand{\pa}{\partial } f(x,t)^2 \,\ud vol_g \right)^{\frac{2}{n+2}}\\&
\le K_b^{\frac{n}{n+2}}  \|f\|_{V^1(\Omega}                \newcommand{\pa}{\partial \times (t_0,t_0+T))}^2.
\end{align*}
Therefore, we complete the proof.
\end{proof}

With the Sobolev inequality in Lemma \ref{lem:sob-parabolic}, we will apply Moser's iterations to the equation \eqref{eq:R-1} to obtain a uniform lower bound of $\mathcal{R}$.
\begin{lem}\label{lem:Rlowerbound} For $t\ge 1$, we have
\[
\mathcal{R}-1\ge -C,
\]
where $C$ is a constant depending only on $\Omega}                \newcommand{\pa}{\partial,n,b$ and $v_0$.

\end{lem}

\begin{proof}
 Let $T>2$, $\frac 12\le T_2<T_1\le1$, $\eta(t)$ be a smooth cut-off function so that $\eta(t)=0$ for all $t<T_2$, $0\le \eta(t)\le 1$ for $t\in [T_2,T_1]$, $\eta(t)=1$ for all $t>T_1$, and $|\eta'(t)|\le \frac{2}{T_1-T_2}$. Denote $\phi= (1-\mathcal{R})^+$.
By \eqref{eq:R-1}, we have
\[
\pa_t (1-\mathcal{R})= \frac{n-2}{n+2} \Delta_g (1-\mathcal{R})-\frac{4}{n+2} (1-\mathcal{R})^2 +\frac{4}{n+2} (1-\mathcal{R}).
\]
Let $k\ge \frac{n}{2}-1$ be a real number. Multiplying both sides of the inequality by $\eta^2 \phi^{1+k}$ and integrating by parts, we see that, for any $0< s<T$,
\begin{align*}
&\frac{1}{2+k}\int_0^s  \int_\Omega}                \newcommand{\pa}{\partial \eta^2 \pa_t \phi^{2+k}  \,\ud vol_g\ud t +  \frac{4(n-2)(k+1)}{(n+2)(k+2)^2} \int_0^s \int_\Omega}                \newcommand{\pa}{\partial \eta^2 |\nabla_g \phi^{\frac{k+2}{2}}|_g^2 \,\ud vol_g\ud t  \\&
\le -\frac{4}{n+2}\int_0^s \int_\Omega}                \newcommand{\pa}{\partial \phi^{3+k} \eta ^2 \,\ud vol_g\ud t  +\frac{4}{n+2} \int_0^s \int_\Omega}                \newcommand{\pa}{\partial \phi^{2+k} \eta ^2 \,\ud vol_g\ud t.
\end{align*}
Note that using \eqref{eq:volumn},  we have
\begin{align*}
&\frac{1}{2+k}\int_0^s \int_\Omega}                \newcommand{\pa}{\partial \eta^2 \pa_t \phi^{2+k} \,\ud vol_g\ud t  \\
&= \frac{1}{2+k} \int_\Omega}                \newcommand{\pa}{\partial \phi^{2+k} \eta^2 \,\ud vol_g \Big|_{t=s} - \frac{1}{2+k} \int_0^s \int_\Omega}                \newcommand{\pa}{\partial \phi^{2+k} \left(2\eta\pa_t \eta +\frac{2n}{n+2} (1-\mathcal{R})\eta^2\right)  \,\ud vol_g\ud t\\
&= \frac{1}{2+k} \int_\Omega}                \newcommand{\pa}{\partial \phi^{2+k} \eta^2 \,\ud vol_g \Big|_{t=s} - \frac{1}{2+k} \int_0^s \int_\Omega}                \newcommand{\pa}{\partial  \left(2\phi^{2+k}\eta\pa_t \eta +\frac{2n}{n+2} \phi^{3+k}\eta^2\right)  \,\ud vol_g\ud t.
\end{align*}
Since  $k\ge \frac{n}{2}-1$, $\frac{1}{2+k}\frac{2n}{n+2}< \frac{4}{n+2}$. Furthermore,
\[
\int_{\Omega}                \newcommand{\pa}{\partial} \mathcal{R} \phi^{2+k}\eta^2 \,\ud vol_g= - \int_{\Omega}                \newcommand{\pa}{\partial} (1-\mathcal{R}) \phi^{2+k}\eta^2\,\ud vol_g + \int_\Omega}                \newcommand{\pa}{\partial \phi^{2+k}\eta^2 \,\ud vol_g\le  \int_\Omega}                \newcommand{\pa}{\partial \phi^{2+k}\eta^2\,\ud vol_g.
\] It follows that
\[
\|\eta \phi^{\frac{2+k}{2}}\|_{V^1(\Omega}                \newcommand{\pa}{\partial\times (0,T))}^2 \le C(2+k)  \int_0^T \int_\Omega}                \newcommand{\pa}{\partial \phi^{2+k} (\eta ^2+|\pa_t \eta| \eta) \,\ud vol_g\ud t,
\]
where $C>0$ depends only on $n$. Making use of Lemma  \ref{lem:sob-parabolic}, we have for all $\gamma:=k+2\ge \frac{n+2}{2}$ that
\[
\left(\int_{T_1}^T \int_\Omega}                \newcommand{\pa}{\partial \phi^{\frac{\gamma(n+2)}{n}} \,\ud vol_g\ud t\right)^{\frac{n}{\gamma(n+2)}} \le \left(\frac{C\gamma}{T_1-T_2}\right)^{\frac{1}{\gamma}}  \left(\int_{T_2}^T \int_\Omega}                \newcommand{\pa}{\partial \phi^{\gamma} \,\ud vol_g\ud t\right)^{\frac{1}{\gamma}}  ,
\]
By the standard Moser's iteration argument,  we have
\[
\sup_{\Omega}                \newcommand{\pa}{\partial\times [1,T]} \phi \le C(n,K_b) \Big(\int_{1/2}^T\int_\Omega}                \newcommand{\pa}{\partial \phi^{\frac{n+2}{2}}\,\ud vol_g\ud t \Big)^{\frac{2}{n+2}},
\]
where $C(n,K_b)>0$ depending only on $n$ and $K_b$. Thus,
\begin{align*}
\sup_{\Omega}                \newcommand{\pa}{\partial\times [1,T]} \phi &\le C(n,K_b) \Big(\int_{1/2}^1\int_\Omega}                \newcommand{\pa}{\partial \phi^{\frac{n+2}{2}}\,\ud vol_g\ud t \Big)^{\frac{2}{n+2}}+C(n,K_b) \Big(\int_{1}^T\int_\Omega}                \newcommand{\pa}{\partial \phi^{\frac{n+2}{2}}\,\ud vol_g\ud t \Big)^{\frac{2}{n+2}}\\
&\le C(n,K_b) \|\mathcal{R}-1\|_{L^\infty (\Omega}                \newcommand{\pa}{\partial \times (1/2,1))}+C(n,K_b) \Big(\int_{1}^T\int_\Omega}                \newcommand{\pa}{\partial \phi^{\frac{n+2}{2}}\,\ud vol_g\ud t \Big)^{\frac{2}{n+2}},
\end{align*}
By Young's inequality, we have
\begin{align*}
\Big(\int_1^T\int_\Omega}                \newcommand{\pa}{\partial \phi^{\frac{n+2}{2}}\,\ud vol_g\ud t \Big)^{\frac{2}{n+2}} &\le (\sup_{\Omega}                \newcommand{\pa}{\partial\times [1,T]} \phi)^{\frac{n-2}{n+2}} \Big(\int_1^T\int_\Omega}                \newcommand{\pa}{\partial \phi^{2}\,\ud vol_g\ud t \Big)^{\frac{2}{n+2}}\\&
\le \varepsilon}           \newcommand{\ud}{\mathrm{d} \sup_{\Omega}                \newcommand{\pa}{\partial\times [1,T]} \phi +C(\varepsilon}           \newcommand{\ud}{\mathrm{d}) \Big(\int_1^T\int_\Omega}                \newcommand{\pa}{\partial \phi^{2}\,\ud vol_g\ud t \Big)^{\frac{1}{2}},
\end{align*}
for any small constant $\varepsilon}           \newcommand{\ud}{\mathrm{d}$. Therefore, by choosing a small $\varepsilon}           \newcommand{\ud}{\mathrm{d}$, we have
\begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:finite-T}
\sup_{\Omega}                \newcommand{\pa}{\partial\times [1,T]} \phi \le C(n,K_b)\left\{  \|\mathcal{R}-1\|_{L^\infty (\Omega}                \newcommand{\pa}{\partial \times (1/2,1))}+ \Big(\int_{1}^T M_2\,\ud t \Big)^{\frac{1}{2}} . \right\}
\ee
By \eqref{eq:F-dcr} and the definition of $\mathcal{R}$, we have
\[
\frac{\ud }{\ud t}F(v(t))= -\frac{2(n-2)}{n+2}M_2(t).
\]
It follows that
\begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:m-2}
\int_0^\infty M_2(t)\,\ud t \le \frac{n+2}{2(n-2)}(F(v(0))- F_\infty)<\infty.
\ee
Moreover, it was proved in \cite{JX19} that   $\|\mathcal{R}-1\|_{L^\infty (\Omega}                \newcommand{\pa}{\partial \times (1/2,1))} \le C$.  Sending $T\to \infty$ in \eqref{eq:finite-T},  we have
\[
\sup_{\Omega}                \newcommand{\pa}{\partial\times [1,\infty)} (1-\mathcal{R})^+ \le C.
\]
Therefore, we complete the proof.
\end{proof}

Using this uniform lower bound of $\mathcal{R}$, we can derive some useful differential inequalities for $M_q$ defined in \eqref{eq:definitionMq}.

For $q>1$,  using Lemma \ref{lem:flow-properties},
 we have
\begin{align*}
\frac{\ud M_q}{\ud t}&= \int_\Omega q |\mathcal{R}-1|^{q-2} (\mathcal{R}-1 )\frac{\pa }{\pa t} (\mathcal{R}-1) \,\ud vol_g -\int_\Omega (\mathcal{R}-1)^{q} \frac{\partial}{\partial t}v^{\frac{2n}{n-2}}\,\ud x   \\&
=q\frac{n-2}{n+2} \int_\Omega |\mathcal{R}-1|^{q-2} ( \mathcal{R}-1) \Delta_g ( \mathcal{R}-1)\,\ud vol_g \\
&\quad + \frac{4q}{n+2}\int_\Omega   |\mathcal{R}-1|^{q} \,\ud vol_g  +\frac{4}{n+2}(q-\frac{n}{2}) \int_\Omega   |\mathcal{R}-1|^{q} ( \mathcal{R}-1)\,\ud vol_g.
\end{align*}
Using Lemma \ref{lem:Rlowerbound},  we have for $t\ge 1$ that
\[
\left|\int_\Omega   |\mathcal{R}-1|^{q} ( \mathcal{R}-1)\,\ud vol_g-\int_\Omega   |\mathcal{R}-1|^{q+1}\,\ud vol_g\right|= 2\int_\Omega   |\mathcal{R}-1|^{q} ( \mathcal{R}-1)^-\,\ud vol_g\le CM_q.
\]
Using Lemma \ref{lem:flow-properties}, we have
\begin{align*}
\int_\Omega |\mathcal{R}-1|^{q-2} ( \mathcal{R}-1) \Delta_g ( \mathcal{R}-1)\,\ud vol_g= -\frac{4(q-1)}{q^2}\int_\Omega  |\nabla_g |\mathcal{R}-1|^{\frac{q}{2}}|_g^2\,\ud vol_g\le 0.
\end{align*}
Therefore, for $q\le \frac{n}{2}$ we have,
\begin{equation}\label{eq:diff-ineq-2}
\frac{\ud M_q}{\ud t} + \frac{4}{n+2}\left(\frac{n}{2}-q\right) M_{q+1}\le CM_q\quad\mbox{for }t\ge 1,
\end{equation}
where $C>0$ is a constant depending on $q$.

For $q\ge \frac{n}{2}$, we first obtain from Lemma \ref{lem:sob-elliptic} that
\begin{align*}
&\int_\Omega |\mathcal{R}-1|^{q-2} ( \mathcal{R}-1) \Delta_g ( \mathcal{R}-1)\,\ud vol_g\\
&= -\frac{4(q-1)}{q^2}\int_\Omega  |\nabla_g |\mathcal{R}-1|^{\frac{q}{2}}|_g^2\,\ud vol_g\\
&\le -\beta M_{\frac{qn}{n-2}}^{\frac{n-2}{n}}  +\frac{4(q-1)}{q^2} \int_\Omega \mathcal{R} |\mathcal{R}-1|^q\,\ud vol_g\\
&\le -\beta M_{\frac{qn}{n-2}}^{\frac{n-2}{n}}  +\frac{4(q-1)}{q^2} \int_\Omega (\mathcal{R}-1) |\mathcal{R}-1|^q\,\ud vol_g +\frac{4(q-1)}{q^2} M_q,
\end{align*}
where $\beta>0$ is a constant depending on $K_b$ and $q$. Then, we have
\[
\frac{\ud M_q}{\ud t} +  \beta M_{\frac{qn}{n-2}}^{\frac{n-2}{n}} \le  \frac{4}{n+2}\left(q-\frac{n}{2}+\frac{(n-2)(q-1)}{q}\right) M_{q+1}+CM_q.
\]
By the interpolation inequality and Young's inequality we have
\begin{align*}
M_{q+1} \le M_{\frac{qn}{n-2}}^{\frac{n-2}{2q}} M_{q}^{\frac{2(q+1)-n}{2q}} \le \varepsilon}           \newcommand{\ud}{\mathrm{d} M_{q(p+1)/2}^{\frac{n-2}{n}} +C(\varepsilon}           \newcommand{\ud}{\mathrm{d}) M_{q}^{\frac{2(q+1)-n}{2q-n}}.
\end{align*}
By choosing a small $\varepsilon}           \newcommand{\ud}{\mathrm{d}$, we obtain
\begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:diff-ineq-3}
\begin{split}
\frac{\ud }{\ud t} M_q(t)  +\beta M_{\frac{qn}{n-2}}(t)^{\frac{n-2}{n}} &\le C\Big(M_q(t)+  M_{q}(t)^{1+\frac{2}{2q-n}}\Big) \quad\mbox{for }t\ge 1
\end{split}
\ee
for $q>\frac{n}{2}$, where $\beta$ and $C$ are positive constants depending on $q$.

The differential inequalities \eqref{eq:diff-ineq-2} and \eqref{eq:diff-ineq-3} will be used to prove the decay of $M_q$.
\begin{prop}\label{prop:Mq} For every $1\le q<\infty$, we have
\[
\lim_{t\to \infty} M_q(t) =0.
\]
\end{prop}

\begin{proof}
By H\"older's inequality and \eqref{eq:BN-5}, we only need to consider $q\ge 2$.

Let us assume $n\ge 4$ first.

Case 1. $2\le q\le \frac{n}{2}$.

Since $M_2\in L^1(0,\infty)$, we can pick $t_j\to \infty$ such that $M_2(t_j) \to 0$ as $j\to \infty$.  By \eqref{eq:diff-ineq-2} we have
\[
\frac{\ud }{\ud t} M_2(t) \le CM_2(t).
\]
Integrating the above inequality we have
\[
M_2(t)\le M_2(t_j) +C \int_{t_j}^\infty M_2(s)\,\ud s \quad \mbox{for }t\ge t_j.
\]
Hence,  $\lim_{t\to \infty} M_2(t)=0$.  If $2<q\le \frac{n}{2}$, \eqref{eq:diff-ineq-2} we have
\[
\int_{1}^\infty M_3(t)\,\ud t \le C \Big(\int_1^\infty M_2(t)\,\ud t +M_2(1)\Big) <\infty.
\]
For any $2<q\le \min\{3,\frac{n}{2}\}$, we have $M_q(t)\le M_2(t)+M_3(t)$. Hence, $\int_{1}^\infty M_q(t)\,\ud t <\infty$. We can repeat the argument for $M_2$ to show that $M_q(t) \to 0$ as $t\to \infty$.  If $3<\frac{n}{2}$, we can show that $\int_1^\infty M_4<\infty$ and $M_{q}(t)\to 0$ as $t\to \infty$ for all $3< q\le  \min\{4,\frac{n}{2}\}$. Repeating this argument in finite times, we then have $M_q \in L^1(1,\infty)$ and $M_{q}(t)\to 0$ as $t\to \infty$ for all $2\le q\le \frac{n}{2}$.

Case 2. $q> \max\{2, \frac{n}{2}\}$.

By \eqref{eq:diff-ineq-2} with $q=n/2$, we have
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:q0}
\int_1^\infty M_{\frac{n^2}{2(n-2)}}(t)^{\frac{n-2}{n}}\,\ud t <\infty.
\ee
Using \eqref{eq:diff-ineq-3} to have
\[
\frac{\ud }{\ud t} M_{q}(t) \le CM_{q}(t)^{\frac{n-2}{n}}\Big(M_{q}(t)^{\frac{2}{n}}+  M_{q}(t)^{\frac{2}{n}+\frac{2}{2q-n}}\Big) \quad\mbox{for }t\ge 1
\]
Hence,
\[
H(M_q(t))\le H(M_{q}(T)) +C\int_T^\infty M_{q}^{\frac{n-2}{n}}\,\ud t \quad \mbox{for } 1\le T<t<\infty
\]
where $$H(\rho)=\int_0^\rho \frac{1}{s^{\frac{2}{n}}+s^{\frac{2}{n}+\frac{2}{2q-n}}}\,\ud s.$$
Let
\[
q_0=\frac{n^2}{2(n-2)},\quad q_k=\frac{n}{n-2}q_{k-1}, \quad k=1,2,\cdots.
\]
Note that $\frac{2}{n}+\frac{2}{2q_0-n}=1$ and $\frac{2}{n}+\frac{2}{2q_k-n}<1$ for all $k\ge 1$. Hence, starting with \eqref{eq:q0} that $\displaystyle\int_1^\infty M_{q_{0}}(t)^{\frac{n-2}{n}}\,\ud t <\infty$, using similar arguments to those in Case 1, we can recursively prove in the order of $k=0,1,2,\cdots$ that $M_{q_k}(t)\to 0$ as $t\to\infty$,  $\displaystyle\int_1^\infty \Big(M_{q_k}(t)+  M_{q_k}(t)^{1+\frac{2}{2q-n}}\Big)<\infty$, $\displaystyle\int_1^\infty M_{q_{k+1}}(t)^{\frac{n-2}{n}}\,\ud t <\infty$, and $M_{q_{k+1}}(t)\to 0$ as $t\to\infty$. Hence, $M_q(t)\to 0$ for any $q\ge \frac{n^2}{2(n-2)}$.

Finally, let us consider $n=3$. By \eqref{eq:diff-ineq-3}, we have
\[
\frac{\ud }{\ud t} M_2(t) \le CM_2(t)(1+M_2(t)^{2}).
\]
Using \eqref{eq:m-2}, we can pick $t_i\to \infty$ such that $M_2(t_i) \to 0$. Hence,
\[
\arctan M_2(t) \le \arctan M_2(t_i) +C \int_{t_i}^\infty M_2(t)\,\ud t.
\] It follows that $\lim_{t\to \infty}\arctan M_2(t) =0$ and thus $\lim_{t\to \infty} M_2(t)=0$. Hence, $\displaystyle\int_1^\infty M_{6}(t)^{\frac{1}{3}}\,\ud t <\infty$. Since $6>q_0$ when $n=3$, we can use the argument of those in Case 2 to show that $M_q(t)\to 0$ for all $q\ge 6$. By H\"older inequality, we conclude that $M_q(t)\to 0$ for all $q\ge 1$.
\end{proof}

\begin{cor}
We have
\[
\lim_{t\to\infty}\|\mathcal{R}-1\|_{L^\infty(\Omega)}=0.
\]
\end{cor}
\begin{proof}
We know from the proof of Proposition \ref{prop:Mq} that there exists a sufficiently large $q>1$ such that
\[
\int_1^\infty M_q(t)\,\ud t<\infty.
\]
This means that $\mathcal{R}-1$ has very high integrability in space-time. Then we can apply the Moser's iteration as in the proof of Lemma \ref{lem:Rlowerbound} to obtain
\[
\|\mathcal{R}-1\|_{L^\infty(\Omega\times(T,\infty))} \le C \left(\int_{T-1}^\infty M_q(t)\,\ud t\right)^{\frac 1q}
\]
for all large $T$. Hence, the corollary follows.
\end{proof}

\section{Concentration compactness}\label{sec:concentration}

The solution of \eqref{eq:BN-3} may blow up as $t\to \infty$ because of the critical exponent $\frac{n+2}{n-2}$. Nevertheless, we also know how the solutions may blow up.

\begin{prop}\label{prop:c-c} Let $v$ be a solution of \eqref{eq:BN-3}. For any $t_\nu \to \infty$, $\nu \to \infty$, $v_\nu=v(\cdot,t_\nu)$ is  a Palais-Smale sequence of the functional $F$ given by \eqref{eq:F-funct} in $H_0^1(\Omega}                \newcommand{\pa}{\partial)$.

\end{prop}

\begin{proof}  We have already proved that $v_\nu$ is bounded in $H^1_0(\Omega}                \newcommand{\pa}{\partial)$ and $F(v_\nu)\to F_\infty$ as $\nu \to \infty$. It remains to show the derivative of $F$ at $v_\nu$ tends to zero. Indeed, for any $\varphi\in H^1_0(\Omega}                \newcommand{\pa}{\partial)$, we have
{\allowdisplaybreaks
\begin{align*}
\langle \ud F(v_\nu ), \varphi \rangle&= 2\int_{\Omega}                \newcommand{\pa}{\partial} (-\Delta v_\nu-bv_\nu- v_\nu ^{\frac{n+2}{n-2}}) \varphi \,\ud x \\&
= 2\int_{\Omega}                \newcommand{\pa}{\partial} (\mathcal{R} -1 ) v_\nu ^{\frac{n+2}{n-2}} \varphi \,\ud x\\&
\le 2 \left(\int_{\Omega}                \newcommand{\pa}{\partial} |\mathcal{R} -1 |^{\frac{2n}{n+2}} v_\nu ^{\frac{2n}{n-2}} \,\ud x\right)^{\frac{n+2}{2n}} \left(\int_{\Omega}                \newcommand{\pa}{\partial} | \varphi|^{\frac{2n}{n-2}} \,\ud x \right)^{\frac{n-2}{2n}} \\&
\le C(n)M_{\frac{2n}{n+2}}(t_\nu)^{\frac{n+2}{2n} } \|\varphi \|_{H^1_0(\Omega}                \newcommand{\pa}{\partial)},
\end{align*}
}
where we used H\"older's inequality  and the Sobolev inequality. It follows from Proposition \ref{prop:Mq} that $\ud F(v_\nu)$ strongly converges to $0$ in $H^{-1}(\Omega}                \newcommand{\pa}{\partial)$.

Therefore, we complete the proof.
\end{proof}

The next proposition shows that the blow up points, if exist, will stay uniformly away from the boundary $\partial\Omega$.

\begin{prop} \label{prop:away} There exist two positive constants  $\delta_0$ and $C$, depending on $v(\cdot,1)$,   such that for all $x\in \Omega}                \newcommand{\pa}{\partial$ with $d(x):=\mbox{dist}(x,\partial\Omega)<\delta_0$ and $t\ge 1$,
\[
v(x,t)\le C d(x).
\]
\end{prop}

\begin{proof} We are going to use the moving plane method as Han \cite{H} did for the elliptic case.  By the Hopf Lemma, there exist $\rho_0>0$ and $\alpha}                \newcommand{\lda}{\lambda_0>0$ such that $v(z-\rho e, 1)$ is nondecreasing for $0<\rho<\rho_0$, where $z\in \pa \Omega}                \newcommand{\pa}{\partial$, $e\in \mathbb{R}}              \newcommand{\Sn}{\mathbb{S}^n^n $ with $|e|=1$, and $(e,\nu(z))\ge \alpha}                \newcommand{\lda}{\lambda_0$ with $\nu(z)$ the unit out normal to $\pa \Omega}                \newcommand{\pa}{\partial $ at $z$. If $\Omega}                \newcommand{\pa}{\partial$ is strictly convex,  using the moving plane method we can conclude that  $v(z-\rho e, t)$ is nondecreasing for $0<\rho<\rho_0$, and for all $t\ge 1$. Therefore, we can find $\gamma>0 $ and $\delta>0$ such that for any fixed $t\ge 1$, and any $x\in\Omega$ satisfying $0<d(x)<\delta$, there exists a measurable set $\Gamma_x$  with (i) $meas(\Gamma_x )\ge \gamma$, (ii) $\Gamma_x \subset \{z: d(z)\ge \delta/2 \}$, and (iii) $v(y,t)\ge v(x,t)$ for any $y\in \Gamma_x$. Actually, $\Gamma_x$ can be taken to a piece of cone with vertex at $x$. It follows that for any $x\in \{z:0<d(z)<\delta\}$, we have
\[
v(x,t) \le \frac{1}{meas(\Gamma_x )} \int_{\Gamma_x}   v(y, t)\,\ud y \le \frac{C}{\gamma},
\]
where we used \eqref{eq:BN-5} and H\"older's inequality. Namely, $v(x,t)\le C$ for $(x,t)\in \{z:0<d(z)<\delta\} \times [1,\infty)$.  By the proof of Theorem 4.1 in \cite{DKV}, we have $v(x,t)\le Cd(x)$ for $(x,t)\in \{z:0<d(z)<\delta\} \times [1,\infty)$.

For a general domain, one can first use a Kelvin transform near each boundary point, and then apply the moving plane method. Pick any point $P\in \pa \Omega}                \newcommand{\pa}{\partial$ for instance. Since we assume the boundary of the domain $\Omega}                \newcommand{\pa}{\partial$ is smooth, we may assume, without loss of generality, that the unit ball $B_1$ contacts $P$ from the exterior of $\Omega}                \newcommand{\pa}{\partial$. Let $w(x,t)$ be the Kelvin transform of $v$:
\[
w(x,t)= |x|^{2-n} v\left(\frac{x}{|x|^2}, t\right).
\]
Then
\[
\begin{cases}
\pa_t w^{\frac{n+2}{n-2}} =\Delta w+b|x|^{-4} w +w^{\frac{n+2}{n-2}} \quad \mbox{in }\Omega}                \newcommand{\pa}{\partial_P \times (0,\infty)\\
w=0 \quad \mbox{on }\pa \Omega}                \newcommand{\pa}{\partial_P \times (0,\infty),
\end{cases}
\]
where $\Omega}                \newcommand{\pa}{\partial_P$ is the image of $\Omega}                \newcommand{\pa}{\partial$ under the Kelvin transform. Since $b\ge 0$,  $b|x|^{-4}$ is nondecreasing along the $-P$ direction. Applying the moving plane method we have that  $w(\cdot,t)$ is nondecreasing along the $-P$ direction in a neighborhood (uniform in $t$) of $P$. Since the $L^{\frac{2n}{n-2}}$ norm is invariant under the Kelvin transform, using the above argument in the case of strictly convex domains, we conclude that $w(\cdot, t)$ is bounded in a neighborhood of $P$ independent of $t$ and so is $v(\cdot, t)$.  It follows that $v(x,t)\le Cd(x)$ for $(x,t)\in \{z:0<d(z)<\delta\} \times [1,\infty)$ for some $\delta>0$.

Therefore, we complete the proof.
\end{proof}

 For $a\in \mathbb{R}}              \newcommand{\Sn}{\mathbb{S}^n^n$ and $\lda\in (0,\infty)$, let
\begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:bubble-def}
\bar \xi_{a,\lda}(x)=c_0 \left(\frac{\lda}{1+\lda^2|x-a|^2}\right)^{\frac{n-2}{2}}
\ee
with $c_0=(n(n-2))^{\frac{n-2}{4}}$. Then we have
\[
-\Delta \bar \xi_{a,\lda} = \bar \xi_{a,\lda}^{\frac{n+2}{n-2}} \quad \mbox{in }\mathbb{R}}              \newcommand{\Sn}{\mathbb{S}^n^n
\]
and
\[
\int_{\mathbb{R}}              \newcommand{\Sn}{\mathbb{S}^n^n} \bar \xi_{a,\lda} ^{\frac{2n}{n-2}}= Y(\Sn)^{\frac{n}{2}},
\]
where $\Sn$ is the standard unit sphere in $\mathbb{R}}              \newcommand{\Sn}{\mathbb{S}^n^{n+1}$,
\[
Y(\Sn)= \frac{n(n-2)}{4} |\Sn|^{\frac{2}{n}}=\inf_{u\in H^{1}(\mathbb{S}^n )} \frac{\int_{\Sn } |\nabla u|^2 +\frac{n(n-2)}{4} u^2\ud vol_{g_{\mathbb{S}^n}}}{(\int_{\Sn} | u|^{\frac{2n}{n-2}} vol_{g_{\mathbb{S}^n}} )^{\frac{n-2}{n}}}
\]
and $|\Sn|$ is the area of $\Sn$. Define
 \begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:xi}
\xi_{a,\lda}(x)=  \bar \xi_{a,\lda}(x)-h_{a,\lda}(x),
\ee
where $\Delta h_{a,\lda}(x) = 0$ in $\Omega}                \newcommand{\pa}{\partial$ and $h_{a,\lda}= \bar \xi_{a,\lda} $ on $\pa \Omega}                \newcommand{\pa}{\partial$.  By the maximum principle, $\xi_{a,\lda}>0$ in $\Omega$ and $h_{a,\lda}>0$ in $\overline\Omega$.

\begin{prop}\label{prop:s-b-c} Let $v$ be a solution of \eqref{eq:BN-3}.  For any $t_\nu \to \infty$, $\nu \to \infty$,  after passing to a subsequence if necessary, $v_\nu$ weakly converges to $v_\infty$ in $H_0^1(\Omega}                \newcommand{\pa}{\partial)$ and  we can find a non-negative integer $m$ and a sequence of m-tuplets $(x^*_{k,\nu}, \lda_{k,\nu}^*)_{1\le k\le m}$, $(x^*_{k,\nu}, \lda_{k,\nu}^*) \in \Omega}                \newcommand{\pa}{\partial \times (0,\infty)$, with the following properities.
\begin{itemize}

 \item[-] The function $v_\infty\in H_0^1(\Omega}                \newcommand{\pa}{\partial)$ satisfies the equation $-\Delta v_\infty-bv_\infty= v_\infty^{\frac{n+2}{n-2}} $ in $\Omega}                \newcommand{\pa}{\partial$.

\item[-] There hold, for all $i\neq j$,
\[
\frac{ \lda_{i, \nu}^* }{\lda_{j, \nu}^* }+\frac{ \lda_{j, \nu}^* }{\lda_{i, \nu}^* }+\lda_{i, \nu}^* \lda_{j, \nu}^* |x^*_{i,\nu}-x^*_{j,\nu} |^2 \to \infty,
\]
and for all $k$
\[
\lda_{k, \nu}^* d(x^*_{k,\nu}) \to \infty
\]
as $\nu \to \infty$.

\item[-] We have
\[
\left\|v_\nu- v_\infty-\sum_{k=1}^m \xi_{x^*_{k,\nu}, \lda_{k,\nu}^*}\right\| \to 0
\]
as $\nu\to \infty$.

\item[-] We have
\[
F(v_\nu) = F(v_\infty) +\frac{2m}{n} Y(\mathbb{S}^n)^{n/2} +o(1),
\]
where $o(1)\to 0$ as $\nu\to \infty$.

\end{itemize}

\end{prop}

\begin{proof} In view of Propositions \ref{prop:c-c} and \ref{prop:away}, this proposition follows from Struwe \cite{St} and Bahri-Coron \cite{Bahri-C}.
\end{proof}
A similar result for the harmonic map heat flow was proved by  Qing-Tian \cite{QT}.
The correction term $h_{a,\lambda}$ in \eqref{eq:xi} is small and can be controlled.

\begin{lem} \label{lem:bubble-derivative}  Let $\xi_{a,\lda}$ and $h_{a,\lda}$ be defined as in \eqref{eq:xi}. Suppose $a\in \Omega}                \newcommand{\pa}{\partial$ with $d(a)>\delta>0$ and $\lda>1$. Then we have
\[
|h_{a,\lda}(x)|+|\nabla_{a} h_{a,\lda}(x)|+\lda |\nabla_{\lda} h_{a,\lda}(x)| \le C(n, \Omega}                \newcommand{\pa}{\partial,\delta) \lda^{-\frac{n-2}{2}}, \quad x\in \Omega}                \newcommand{\pa}{\partial,
\]
\[
\nabla_{a} \xi_{a,\lda} (x)=(n-2)\xi_{a,\lda} \frac{\lda^2(x-a)}{1+\lda^2|x-a|^2} +O(\lda^{-\frac{n-2}{2}}),
\]
and
\[
\nabla_{\lda} \xi_{a,\lda} (x)=\frac{(n-2)}{2\lda } \xi_{a,\lda} \frac{1-\lda^2|x-a|^2}{1+\lda^2 |x-a|^2} +O(\lda^{-\frac{n}{2}}),
\]
where $|O(\lda^{-\frac{n-2}{2}}) |\le C \lda^{-\frac{n-2}{2}}$ for some $C$ depending only on $n, \Omega}                \newcommand{\pa}{\partial$ and $\delta$.
\end{lem}

\begin{proof} It suffices to prove the estimate of $h_{a,\lda}$, which follows from Poisson formula for the Laplace equation.
\end{proof}

\section{Refined blow up analysis}\label{sec:blowup}

We continue from Proposition \ref{prop:s-b-c}. By the strong maximum principle, the nonnegative limit $v_\infty$ either is positive in $\Omega}                \newcommand{\pa}{\partial$ or identically equals to zero. We will adapt the refined blow up analysis in Brendle \cite{Br05} to the equation \eqref{eq:BN-3}.

\subsection{The case $v_\infty \equiv 0$}\label{sec: zero}

For every $\nu$, let $\mathcal{A}_\nu$ be the closed set of all $m$-tuplets $(x_k,\lda_k, \alpha}                \newcommand{\lda}{\lambda_k)_{1\le k\le m}$ satisfying $(x_k,\lda_k, \alpha}                \newcommand{\lda}{\lambda_k) \in \overline B_{\frac{1}{\lda_{k,\nu}^*}}(x_{k,\nu}^*) \times [\frac{\lda_{k,\nu}^*}{2}, \frac{3\lda_{k,\nu}^*}{2} ] \times [\frac{1}{2},\frac32]$.  Choose an $m$-tuplet $(x_{k,\nu},\lda_{k,\nu}, \alpha}                \newcommand{\lda}{\lambda_{k,\nu})_{1\le k\le m}\in \mathcal{A}_\nu$ such that
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:finite-minimum}
\left\|v_\nu- \sum_{k=0}^m \alpha}                \newcommand{\lda}{\lambda_{k,\nu} \xi_{x_{k,\nu}, \lda_{k,\nu}}\right\| =\inf_{(x_k,\lda_k, \alpha}                \newcommand{\lda}{\lambda_k)_{1\le k\le m} \in \mathcal{A}_{\nu}} \left\|v_\nu- \sum_{k=0}^m \alpha}                \newcommand{\lda}{\lambda_{k} \xi_{x_{k}, \lda_{k}}\right\|.
\ee
By Proposition \ref{prop:s-b-c}, Proposition \ref{prop:away} and the definition of $(x_{k,\nu},\lda_{k,\nu}, \alpha}                \newcommand{\lda}{\lambda_{k,\nu})_{1\le k\le m}$, we have, for all $i\neq j$,
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:bubbles-1}
\frac{ \lda_{i, \nu} }{\lda_{j, \nu}}+\frac{ \lda_{j, \nu} }{\lda_{i, \nu}}+\lda_{i, \nu} \lda_{j, \nu} |x_{i,\nu}-x_{j,\nu} |^2 \to \infty,
\ee
and for all $k$
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:bubbles-2}
\lda_{k, \nu} d(x_{k,\nu}) \to \infty
\ee
as $\nu \to \infty$.
In addition, $d(x_{k,\nu})>\delta/2$, and
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:bubbles-3}
\left\|v_\nu-\sum_{k=1}^m  \alpha}                \newcommand{\lda}{\lambda_k\xi_{x_{k,\nu}, \lda_{k,\nu}}\right\| \to 0
\ee
as $\nu\to \infty$.

By the triangle inequality,
\begin{align*}
&\left\|\sum_{k=1}^m  \alpha}                \newcommand{\lda}{\lambda_k\xi_{x_{k,\nu}, \lda_{k,\nu}}- \sum_{k=1}^m \xi_{x^*_{k,\nu}, \lda_{k,\nu}^*} \right\| \\& \le \left\|v_\nu-\sum_{k=1}^m  \alpha}                \newcommand{\lda}{\lambda_k\xi_{x_{k,\nu}, \lda_{k,\nu}} \right\|+\left\|v_\nu- \sum_{k=1}^m \xi_{x^*_{k,\nu}, \lda_{k,\nu}^*} \right\|=o(1).
\end{align*}
It follows that, for all $1\le k\le m$,
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:bubbles-4}
|x_{k,\nu}- x_{k,\nu}^{*}|=o(1) \frac{1}{\lda_{k,\nu}^*},\quad  \frac{\lda_{k,\nu}}{\lda_{k,\nu}^*}=1+o(1), \quad \alpha}                \newcommand{\lda}{\lambda_{k,\nu}=1+o(1).
\ee
In particular, $(x_{k,\nu},\lda_{k,\nu}, \alpha}                \newcommand{\lda}{\lambda_{k,\nu})_{1\le k\le m}$ is an interior point of $\mathcal{A}_{\nu}$.

In the sequel, we assume
\begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:order}
\lda_{1,\nu}\ge \lda_{2,\nu} \ge \dots \ge \lda_{m,\nu}.
\ee
 Let
\begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:decomposition-1}
U_\nu= \sum_{k=1}^m  \alpha}                \newcommand{\lda}{\lambda_{k,\nu}\xi_{x_{k,\nu}, \lda_{k,\nu}}, \quad w_\nu=v_{\nu}-U_{\nu}.
\ee

\begin{lem} \label{lem:proj-1} We have for $1\le k\le m$,
\begin{align*}
&\Big|\int_\Omega}                \newcommand{\pa}{\partial \xi_{x_{k,\nu}, \lda_{k,\nu}}^{\frac{n+2}{n-2}}  w_\nu \,\ud x\Big|+\Big|\int_\Omega}                \newcommand{\pa}{\partial \xi_{x_{k,\nu}, \lda_{k,\nu}}^{\frac{n+2}{n-2}}  \frac{1-\lda^2|x-x_{k,\nu}|^2}{1+\lda^2|x-x_{k,\nu}|^2} w_\nu \,\ud x\Big |\\[2mm] &\quad  +\Big|\int_\Omega}                \newcommand{\pa}{\partial \xi_{x_{k,\nu}, \lda_{k,\nu}}^{\frac{n+2}{n-2}}  \frac{\lda^2(x-x_{k,\nu})}{1+\lda^2|x-x_{k,\nu}|^2} w_\nu \,\ud x\Big | \le o(1) \Big(\int_{\Omega}                \newcommand{\pa}{\partial} |w_\nu|^{\frac{2n}{n-2}}\,\ud x\Big)^{\frac{n-2}{2n}}
\end{align*}
\end{lem}

\begin{proof} By  the finite dimensional variational problem  \eqref{eq:finite-minimum} and \eqref{eq:bubbles-4}, taking derivatives in $\mathcal{A}_\nu$, we have
\[
\int_{\Omega}                \newcommand{\pa}{\partial}\Big[\nabla (\nabla_{a,\lda}\xi_{x_{k,\nu},\lda_{k,\nu}}) \nabla w_\nu-b\nabla_{a,\lda}\xi_{x_{k,\nu},\lda_{k,\nu}} w_\nu\Big]\,\ud x=0
\]
and
\[
\int_{\Omega}                \newcommand{\pa}{\partial}\Big[\nabla \xi_{x_{k,\nu},\lda_{k,\nu}} \nabla w_\nu-b\xi_{x_{k,\nu},\lda_{k,\nu}} w_\nu\Big]\,\ud x=0,
\]
where $\nabla_{a,\lda}\xi_{x_{k,\nu},\lda_{k,\nu}}= \nabla_{a,\lda}\xi_{a,\lda}\Big|_{(a,\lda)=(x_{k,\nu},\lda_{k,\nu})}$.  Integrating by parts, using the equation of $\bar \xi_{a,\lda}$, H\"older's inequality and Lemma \ref{lem:bubble-derivative}, the lemma follows.
\end{proof}

Note that the bubbles are non-degenerate, since we have the following well known lemma (see (3.14) in Rey \cite{Rey}).

\begin{lem}\label{lem:coercive} Let $\bar \xi_{a,\lda}$ be defined in \eqref{eq:bubble-def}. Then there exists a constant $c_1>0$ such that
\[
(1-c_1)\int_{\mathbb{R}}              \newcommand{\Sn}{\mathbb{S}^n^n} |\nabla \varphi|^2 \ge \frac{n+2}{n-2}\int_{\mathbb{R}}              \newcommand{\Sn}{\mathbb{S}^n^n} \bar \xi_{0,1}^\frac{4}{n-2} \varphi^2
\]
for any $\varphi\in H^{1}_0(\mathbb{R}}              \newcommand{\Sn}{\mathbb{S}^n^n)$ satisfying
\[
\int_{\mathbb{R}}              \newcommand{\Sn}{\mathbb{S}^n^n}  \bar \xi_{0,1}^\frac{4}{n-2} (\nabla_{a,\lda} \bar \xi_{0,1}) \varphi\,\ud x=0.
\]
\end{lem}
We have the following non-degeneracy estimates of the second variation of $F$ for $w_\nu$.
\begin{lem} \label{lem:2nd-variational} For large $\nu$, we have
\[
\frac{n+2}{n-2} \int_\Omega}                \newcommand{\pa}{\partial \sum_{k=1}^m  \xi_{x_{k,\nu},\lda_{k,\nu}}  ^{\frac{4}{n-2}} w_\nu^2   \le (1-c) \int_{\Omega}                \newcommand{\pa}{\partial} (|\nabla w_\nu|^2-bw_\nu^2)\,\ud x,
\]
where $c>0$ is independent of $\nu$.

\end{lem}

\begin{proof} We assume $w_\nu$ is not zero, otherwise there is nothing to prove. Define $\tilde w_\nu= \frac{w_\nu}{\|w_\nu\|}$. Suppose the lemma is not true.  Then we can find a subsequence of $\{\tilde w_\nu\}$ (still denoted by $\{\tilde w_\nu\}$) satisfying
\begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:coer-1}
\lim_{\nu\to \infty}\frac{n+2}{n-2}  \int_\Omega}                \newcommand{\pa}{\partial   \sum_{k=1}^m \xi_{x_{k,\nu},\lda_{k,\nu}}  ^{\frac{4}{n-2}} \tilde w_\nu^2  \ge 1.
\ee
By \eqref{eq:sobolev-b},
\begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:coer-2}
\int_{\Omega}                \newcommand{\pa}{\partial} |\tilde w_\nu|^{\frac{2n}{n-2}} \le K_b^{\frac{n}{n-2}} \|\tilde w_\nu\|= K_b^{\frac{n}{n-2}}.
\ee
By \eqref{eq:bubbles-1} and \eqref{eq:order}, we can find $R_\nu\to \infty$, $R_\nu \lda_{j,\nu}^{-1} \to 0 $ for all $1\le j\le m$, and
\begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:coer-3}
\frac{\lda_{i,\nu}}{R_\nu} (\lda_{j,\nu}^{-1} +|x_{i,\nu}-x_{j,\nu}|) \to \infty
\ee
for all $i<j$. Set
\[
\Omega}                \newcommand{\pa}{\partial_{j,\nu}= B_{R_\nu \lda_{j,\nu}^{-1}}(x_{j,\nu})\setminus \bigcup_{i=1}^{j-1} B_{R_\nu \lda_{i,\nu}^{-1}}(x_{i,\nu}).
\]
By \eqref{eq:coer-1} and $\|\tilde w_\nu\|=1$, we can find $1\le j\le m$ such that
\[
\lim_{\nu\to \infty } \int_\Omega}                \newcommand{\pa}{\partial    \xi_{x_{j,\nu},\lda_{j,\nu}}  ^{\frac{4}{n-2}} \tilde w_\nu^2 >0
\]
and
\[
\lim_{\nu\to \infty  } \int_{\Omega}                \newcommand{\pa}{\partial_{j,\nu}} (|\nabla \tilde w_\nu|^2- b\tilde w_\nu^2 ) \le \lim_{\nu\to \infty  }\frac{n+2}{n-2}  \int_\Omega}                \newcommand{\pa}{\partial    \xi_{x_{j,\nu},\lda_{j,\nu}}  ^{\frac{4}{n-2}} \tilde w_\nu^2.
\]
Let $\hat w_\nu(x)= \lda_{j,\nu}^{-\frac{n-2}{2}} \tilde w_\nu(x_{j,\nu} +\lda_{j,\nu}^{-1} x)$. Under this scaling, by using \eqref{eq:coer-3}, we know that either $B_{R_\nu \lda_{j,\nu}^{-1}}(x_{j,\nu})$ will be disjoint with $\bigcup_{i=1}^{j-1} B_{R_\nu \lda_{i,\nu}^{-1}}(x_{i,\nu})$, or $B_{R_\nu \lda_{i,\nu}^{-1}}(x_{i,\nu})$ will shrink to a point for every $1\le i\le j-1$.  By passing to a weak limit in $H^1_{loc}(\mathbb{R}}              \newcommand{\Sn}{\mathbb{S}^n^n)$, and using the above two inequalities and Lemma \ref{lem:proj-1}, we then obtain a contradiction to Lemma \ref{lem:coercive}.

Therefore, Lemma \ref{lem:2nd-variational} is proved.
\end{proof}

\begin{cor}\label{cor:coercive}
For large $\nu$, we have
\[
\frac{n+2}{n-2} \int_\Omega}                \newcommand{\pa}{\partial U_\nu ^{\frac{4}{n-2}} w_\nu^2   \le (1-c) \int_{\Omega}                \newcommand{\pa}{\partial} (|\nabla w_\nu|^2-bw_\nu^2)\,\ud x,
\]
where $c>0$ is independent of $\nu$.
\end{cor}
\begin{proof}
It follows from Lemma \ref{lem:2nd-variational}, H\"older's inequality, the Sobolev inequality \eqref{eq:sobolev-b} and the fact that
\[
\int_\Omega \Big|U_\nu ^{\frac{4}{n-2}} - \sum_{k=1}^m  \xi_{x_{k,\nu},\lda_{k,\nu}}  ^{\frac{4}{n-2}}\Big|^{\frac{n}{2}}=o(1).
\]
\end{proof}

Now we can have an expansion of the Hamitonian $F$ defined in \eqref{eq:F-funct}.
\begin{prop} \label{prop:bubble-energy} When $n\ge 4$ and  $\nu$ is sufficiently large,  we have
\[
F(U_\nu) \le  \sum_{k=1}^m F( \xi_{x_{k,\nu},\lda_{k,\nu}}  ) +o(1) \sum_{k=1}^m  \int_{\Omega}                \newcommand{\pa}{\partial} \xi_{x_{k,\nu},\lda_{k,\nu}}^2   +C \sum_{k=1}^m \lda_{k,\nu}^{2-n}.
\]
\end{prop}

\begin{proof}  We shall need the following inequality
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:basic}
\left(\sum_{k=1}^m a_k \right)^{\frac{2n}{n-2}} \ge \sum_{k=1}^m a_k ^{\frac{2n}{n-2}}+\frac{2n}{n-2}\sum_{k< l} a_k^{\frac{n+2}{n-2}} a_l + c_{n,m}\sum_{k< l} (a_k \lor a_l) ^{\frac{4}{n-2}} (a_k\land a_l)^{2}
\ee for  any $a_1,\dots, a_m\ge 0$, where  $c_{n,m}>0$ is a constant, and $a_k \lor a_l=\max(a_k,a_l)$ and $a_k\land a_l=\min(a_k,a_l)$. This inequality can be proved using Lemma \ref{lem:basic-1} and induction.

Using the  inequality \eqref{eq:basic},  we have
\begin{align}
&\int_\Omega (| \nabla U_\nu|^2 -bU_\nu^2) -\frac{n-2}{n} \int_\Omega U_\nu^{\frac{2n}{n-2}} \nonumber\\
&\le \sum_{k} \alpha}                \newcommand{\lda}{\lambda_{k,\nu}^2 \int_\Omega ( |\nabla \xi_{x_{k,\nu},\lda_{k,\nu}} |^2-b\xi_{x_{k,\nu},\lda_{k,\nu}} ^2)- \sum_{k} \alpha}                \newcommand{\lda}{\lambda_{k,\nu}^{\frac{2n}{n-2}} \frac{n-2}{n}\int_\Omega\xi_{x_{k,\nu},\lda_{k,\nu}}^{\frac{2n}{n-2}} \nonumber\\
&\quad +2 \sum_{i<j}  \alpha}                \newcommand{\lda}{\lambda_{j,\nu} \left[\alpha}                \newcommand{\lda}{\lambda_{i,\nu} \int_\Omega (\nabla \xi_{x_{i,\nu},\lda_{i,\nu}} \nabla \xi_{x_{j,\nu},\lda_{j,\nu}} -b\xi_{x_{i,\nu},\lda_{i,\nu}} \xi_{x_{j,\nu},\lda_{j,\nu}} )-  \alpha}                \newcommand{\lda}{\lambda_{i,\nu}^{\frac{n+2}{n-2}}  \int_\Omega   \xi_{x_{i,\nu},\lda_{i,\nu}} ^{\frac{n+2}{n-2}}\xi_{x_{j,\nu},\lda_{j,\nu}}\right] \nonumber\\
&\quad-c_{n,m} \sum_{i<j}  \int_\Omega   (\xi_{x_{i,\nu},\lda_{i,\nu}} \lor \xi_{x_{j,\nu},\lda_{j,\nu}})^{\frac{4}{n-2}} (\xi_{x_{i,\nu},\lda_{i,\nu}} \land \xi_{x_{j,\nu},\lda_{j,\nu}})^{2}.\label{eq:computeFU1}
\end{align}
By the equation of $\bar \xi_{x_{k,\nu},\lda_{k,\nu}}$ and the definition of $\xi_{x_{k,\nu},\lda_{k,\nu}}$, we have
\begin{align}
&\alpha}                \newcommand{\lda}{\lambda_{k,\nu}^2 \int_\Omega ( |\nabla \xi_{x_{k,\nu},\lda_{k,\nu}} |^2-b\xi_{x_{k,\nu},\lda_{k,\nu}} ^2)- \alpha}                \newcommand{\lda}{\lambda_{k,\nu}^{\frac{2n}{n-2}} \frac{n-2}{n}\int_\Omega\xi_{x_{k,\nu},\lda_{k,\nu}}^{\frac{2n}{n-2}}  \nonumber\\
&\le \alpha}                \newcommand{\lda}{\lambda_{k,\nu}^2 \int_{\mathbb{R}}              \newcommand{\Sn}{\mathbb{S}^n^n}  |\nabla \bar \xi_{x_{k,\nu},\lda_{k,\nu}} |^2- \alpha}                \newcommand{\lda}{\lambda_{k,\nu}^{\frac{2n}{n-2}} \frac{n-2}{n}\int_{\mathbb{R}}              \newcommand{\Sn}{\mathbb{S}^n^n} \bar \xi_{x_{k,\nu},\lda_{k,\nu}}^{\frac{2n}{n-2}}  - b \alpha}                \newcommand{\lda}{\lambda_{k,\nu}^2 \int_\Omega  \xi_{x_{k,\nu},\lda_{k,\nu}} ^2  +C \lda_{k,\nu}^{2-n}\nonumber\\
&\le \int_{\mathbb{R}}              \newcommand{\Sn}{\mathbb{S}^n^n}  |\nabla \bar \xi_{x_{k,\nu},\lda_{k,\nu}} |^2 - \frac{n-2}{n}\int_{\mathbb{R}}              \newcommand{\Sn}{\mathbb{S}^n^n}  \bar \xi_{x_{k,\nu},\lda_{k,\nu}}^{\frac{2n}{n-2}}  -\frac{4}{n-2}(\alpha}                \newcommand{\lda}{\lambda_{k,\nu}-1)^2 Y(\Sn)^{\frac{n}{2}} \nonumber\\
&\quad -b \alpha}                \newcommand{\lda}{\lambda_{k,\nu}^2 \int_\Omega  \xi_{x_{k,\nu},\lda_{k,\nu}} ^2 +C \lda_{k,\nu}^{2-n}\nonumber\\
&\le  F( \xi_{x_{k,\nu} ,\lda_{k,\nu}}) -\frac{4}{n-2}(\alpha}                \newcommand{\lda}{\lambda_{k,\nu}-1)^2Y(\Sn)^{\frac{n}{2}}+o(1) \int_\Omega  \xi_{x_{k,\nu},\lda_{k,\nu}} ^2 +C\lda_{k,\nu}^{2-n},\label{eq:computeFU2}
\end{align}
where we used $\alpha}                \newcommand{\lda}{\lambda_{k,\nu}=1+o(1)$, $\alpha}                \newcommand{\lda}{\lambda_{k,\nu}^2-\frac{n-2}{n} \alpha}                \newcommand{\lda}{\lambda_{k,\nu}^{\frac{2n}{n-2}}\le  \frac{2}{n}-\frac{4}{n-2}(\alpha}                \newcommand{\lda}{\lambda_{k,\nu}-1)^2 $ and
\[
\int_{\mathbb{R}}              \newcommand{\Sn}{\mathbb{S}^n^n} |\nabla \bar \xi_{x_{k,\nu},\lda_{k,\nu}} |^2= \int_{\mathbb{R}}              \newcommand{\Sn}{\mathbb{S}^n^n} \bar  \xi_{x_{k,\nu},\lda_{k,\nu}} ^{\frac{2n}{n-2}}=Y(\Sn)^{\frac{n}{2}}.
\]
In addition,
\begin{align}
& \alpha}                \newcommand{\lda}{\lambda_{i,\nu}  \int_\Omega (\nabla \xi_{x_{i,\nu},\lda_{i,\nu}} \nabla \xi_{x_{j,\nu},\lda_{j,\nu}} -b\xi_{x_{i,\nu},\lda_{i,\nu}} \xi_{x_{j,\nu},\lda_{j,\nu}} ) -\alpha}                \newcommand{\lda}{\lambda_{i,\nu}^{\frac{n+2}{n-2}}  \int_\Omega   \xi_{x_{i,\nu},\lda_{i,\nu}} ^{\frac{n+2}{n-2}}\xi_{x_{j,\nu},\lda_{j ,\nu}} \nonumber\\
& = \alpha}                \newcommand{\lda}{\lambda_{i,\nu} \int_\Omega (-\Delta \xi_{x_{i,\nu},\lda_{i,\nu}}  -b\xi_{x_{i,\nu},\lda_{i,\nu}} -\alpha}                \newcommand{\lda}{\lambda_{i,\nu}^{\frac{4}{n-2}}  \xi_{x_{i,\nu},\lda_{i,\nu}} ^{\frac{n+2}{n-2}} )\xi_{x_{j,\nu},\lda_{j,\nu}} \nonumber\\&
\le C|  \alpha}                \newcommand{\lda}{\lambda_{i,\nu} -1| \int_\Omega  \xi_{x_{i,\nu},\lda_{i,\nu}} ^{\frac{n+2}{n-2}} \xi_{x_{j,\nu},\lda_{j,\nu}}  + \int_\Omega  \Big(\bar\xi_{x_{i,\nu},\lda_{i,\nu}} ^{\frac{n+2}{n-2}}- \xi_{x_{i,\nu},\lda_{i,\nu}} ^{\frac{n+2}{n-2}}\Big) \xi_{x_{j,\nu},\lda_{j,\nu}} \nonumber\\&
\le C|  \alpha}                \newcommand{\lda}{\lambda_{i,\nu} -1| \int_\Omega  \xi_{x_{i,\nu},\lda_{i,\nu}} ^{\frac{n+2}{n-2}} \xi_{x_{j,\nu},\lda_{j,\nu}}  +O(\lda_{i,\nu} ^{2-n}+\lda_{j,\nu}^{2-n}) \nonumber\\
& \le \frac{2}{n-2}(\alpha}                \newcommand{\lda}{\lambda_{k,\nu}-1)^2  Y(\Sn)^{\frac{n}{2}}+C\Big (\int_\Omega  \xi_{x_{i,\nu},\lda_{i,\nu}} ^{\frac{n+2}{n-2}} \xi_{x_{j,\nu},\lda_{j,\nu}}\Big)^2 +O(\lda_{i,\nu} ^{2-n}+\lda_{j,\nu}^{2-n}),\label{eq:computeFU3}
\end{align}
where $C>0$ is independent of $\nu$, and in the second inequality we used
\begin{align*}
\int_\Omega  \Big(\bar\xi_{x_{i,\nu},\lda_{i,\nu}} ^{\frac{n+2}{n-2}}- \xi_{x_{i,\nu},\lda_{i,\nu}} ^{\frac{n+2}{n-2}}\Big) \xi_{x_{j,\nu},\lda_{j,\nu}}
&\le C \int_\Omega  \bar\xi_{x_{i,\nu},\lda_{i,\nu}} ^{\frac{4}{n-2}}|h_{x_{i,\nu},\lda_{i,\nu}} |  \bar\xi_{x_{j,\nu},\lda_{j,\nu}} \\
&\le C\lda_{i,\nu}^{\frac{2-n}{2}} \left(\int_\Omega  \bar\xi_{x_{i,\nu},\lda_{i,\nu}} ^{\frac{n+2}{n-2}}\right)^{\frac{4}{n+2}} \left(\int_\Omega\bar\xi_{x_{j,\nu},\lda_{j,\nu}} ^{\frac{n+2}{n-2}}\right)^{\frac{n-2}{n+2}}\\
&\le C\lda_{i,\nu}^{\frac{2-n}{2}}\lda_{i,\nu}^{\frac{2-n}{2}\frac{4}{n+2}}\lda_{j,\nu}^{\frac{2-n}{2}\frac{n-2}{n+2}}\\
&\le C\lda_{i,\nu}^{\frac{2-n}{2}}(\lda_{i,\nu}^{\frac{2-n}{2}}+\lda_{j,\nu}^{\frac{2-n}{2}})\\
&\le C(\lda_{i,\nu}^{2-n}+\lda_{j,\nu}^{2-n}).
\end{align*}
Combining \eqref{eq:computeFU1}, \eqref{eq:computeFU2} and  \eqref{eq:computeFU3}, we have
\begin{align*}
&F(U_\nu) \\
&\le  \sum_{k=1}^m F( \xi_{x_{k,\nu},\lda_{k,\nu}}  ) +o(1) \sum_{k=1}^m  \int_{\Omega}                \newcommand{\pa}{\partial} \xi_{x_{k,\nu},\lda_{k,\nu}}^2   +C \sum_{k=1}^m \lda_{k,\nu}^{2-n}\\
& \ \ + \sum_{i<j} \left[C\Big (\int_\Omega  \xi_{x_{i,\nu},\lda_{i,\nu}} ^{\frac{n+2}{n-2}} \xi_{x_{j,\nu},\lda_{j,\nu}}\Big)^2 - c_{n,m}  \int_\Omega   (\xi_{x_{i,\nu},\lda_{i,\nu}} \lor \xi_{x_{j,\nu},\lda_{j,\nu}})^{\frac{4}{n-2}} (\xi_{x_{i,\nu},\lda_{i,\nu}} \land \xi_{x_{j,\nu},\lda_{j,\nu}})^{2}\right].
\end{align*}
Meanwhile, we have
\[
\begin{split}
&C\Big (\int_\Omega  \xi_{x_{i,\nu},\lda_{i,\nu}} ^{\frac{n+2}{n-2}} \xi_{x_{j,\nu},\lda_{j,\nu}}\Big)^2 - c_{n,m}  \int_\Omega   (\xi_{x_{i,\nu},\lda_{i,\nu}} \lor \xi_{x_{j,\nu},\lda_{j,\nu}})^{\frac{4}{n-2}} (\xi_{x_{i,\nu},\lda_{i,\nu}} \land \xi_{x_{j,\nu},\lda_{j,\nu}})^{2} \\
& \le C\Big (\int_{\mathbb{R}}              \newcommand{\Sn}{\mathbb{S}^n^n} \bar \xi_{x_{i,\nu},\lda_{i,\nu}} ^{\frac{n+2}{n-2}} \bar\xi_{x_{j,\nu},\lda_{j,\nu}}\Big)^2 - c_{n,m}  \int_{\mathbb{R}}              \newcommand{\Sn}{\mathbb{S}^n^n}   (\bar\xi_{x_{i,\nu},\lda_{i,\nu}} \lor \bar\xi_{x_{j,\nu},\lda_{j,\nu}})^{\frac{4}{n-2}} (\bar\xi_{x_{i,\nu},\lda_{i,\nu}} \land \bar\xi_{x_{j,\nu},\lda_{j,\nu}})^{2} \\
&\quad+C(\lda_{i,\nu}^{2-n}+\lda_{j,\nu}^{2-n})\\
&\le C(\lda_{i,\nu}^{2-n}+\lda_{j,\nu}^{2-n})\quad\mbox{for all large }\nu,
\end{split}
\]
where we used \eqref{eq:I1I2} in the last inequality.

Therefore, we complete the proof.
\end{proof}

\begin{cor}\label{cor:pure-singular}
If $n\ge 4$ and $b>0$ satisfying \eqref{eq:b}, we have, for large $\nu$,
\[
F(U_\nu) \le \frac{2m}{n} Y(\Sn)^{\frac{n}{2}}.
\]

\end{cor}

\begin{proof} Note that
\[
\int_{\Omega}                \newcommand{\pa}{\partial}\xi_{x_{j,\nu},\lda_{j,\nu}}^{2}\ge\begin{cases}  \frac{1}{C} \lda_{j,\nu}^{-\frac{4}{n-2}},& \quad \mbox{if }n\neq 4, \\
 \frac{1}{C}  \lda_{j,\nu}^{-2} \ln \lda_{j,\nu},& \quad \mbox{if }n=4.
\end{cases}
\]
Hence, if $n\ge 4$, for any large constant $N$ we can find $j_N>0$ such that for all $j\ge j_N$ there holds   $\int_{\Omega}                \newcommand{\pa}{\partial}\xi_{x_{j,\nu},\lda_{j,\nu}}^{2} \ge N \lda_{j,\nu}^{2-n}$. The corollary follows immediately from Proposition \ref{prop:bubble-energy}.
\end{proof}

\subsection{The case $v_\infty >0$}

Let $\mathcal{L}^2(\Omega):=\big\{f: \int_{\Omega}f^2v_\infty^{\frac{4}{n-2}}<\infty\big\}$ with the inner product $\langle f,g\rangle=\int_{\Omega}fgv_\infty^{\frac{4}{n-2}}$. Then the operator
\[
f \longmapsto \left[v_\infty^{-\frac{4}{n-2}}(-\Delta-b)\right]^{-1}f
\]
is a bounded linear compact symmetric operator mapping $\mathcal{L}^2(\Omega)$ into itself. Using the spectral theorem, there exists a sequence of $H^1_0(\Omega)$ functions $\{\phi_l: l\in\mathbb{N}\}$ and a sequence of positive real numbers $\{\mu_l:l\in\mathbb{N}\}$ such that $0<\mu_1<\mu_2\le \mu_3 \le \cdots \to\infty$,
\[
-\Delta \phi_l -b\phi_l =\mu_l v_\infty^{\frac{4}{n-2}} \phi_l \quad \mbox{in }\Omega}                \newcommand{\pa}{\partial, \quad \phi_l=0 \quad \mbox{on }\pa \Omega}                \newcommand{\pa}{\partial,
\]
 and $\{\phi_l: l\in\mathbb{N}\}$ forms an orthonormal  basis of $\mathcal{L}^2(\Omega)$. In particular,
\[
\int_{\Omega}                \newcommand{\pa}{\partial} v_\infty^{\frac{4}{n-2}} \phi_i \phi_j=\begin{cases}
&1 \quad \mbox{for }i=j, \\
&0 \quad \mbox{for }i\neq j.
\end{cases}
\] By the regularity theory of linear elliptic equations, $\phi_l\in C^{2+\frac{4}{n-2}}(\bar \Omega}                \newcommand{\pa}{\partial) \cap C^\infty(\Omega}                \newcommand{\pa}{\partial)$ for every $l$.  By the equation of $v_\infty$ and the positivity of  $v_\infty$, we know that $\mu_1=1$ and $\phi_1= v_\infty(\int_{\Omega}                \newcommand{\pa}{\partial} v_\infty^{\frac{2n}{n-2}})^{-1/2}$. It is easy to check that $\{\frac{1}{\sqrt{\mu_l}}\phi_l\}$ is also an orthonormal basis of $H_0^1(\Omega}                \newcommand{\pa}{\partial)$ with respect to the inner product \eqref{eq:inner-product}.

Let $L$ be the largest number such that
\[
\mu_l\le \frac{n+2}{n-2} \quad \mbox{for all } l\le L.
\]
For $f\in L^p(\Omega)$, $p\ge 1$, we denote by $\Pi$ the projection operator
\[
\Pi
= f-\sum_{i=1}^L \left(\int_{\Omega}                \newcommand{\pa}{\partial }f \phi_i\,\ud x\right) v_\infty ^{\frac{4}{n-2}} \phi_i.
\]
It is clear that $\Pi( L^p(\Omega}                \newcommand{\pa}{\partial))=\{f\in L^p(\mathbb{R}}              \newcommand{\Sn}{\mathbb{S}^n^n): \int_\Omega f\phi_i=0, \ i=1,2,\cdots,L\}.$ Hence, $\Pi( L^p(\Omega}                \newcommand{\pa}{\partial))$ is a closed subspace of $L^p(\Omega}                \newcommand{\pa}{\partial)$, and thus, is a Banach space with the inherited $L^p$ norm.

We have several estimates regarding this projection.

\begin{lem}\label{lem:Lpprojection} For every $1\le p<\infty$, we can find a constant $C$ such that
\[
\|f\|_{L^p(\Omega}                \newcommand{\pa}{\partial)} \le C \Big \|\Delta f+bf+\frac{n+2}{n-2} v_\infty^{\frac{4}{n-2}} f \Big\|_{L^p(\Omega}                \newcommand{\pa}{\partial)}+C \sup_{1\le l\le L} \Big|\int_\Omega}                \newcommand{\pa}{\partial v_\infty^{\frac{4}{n-2}} \phi_l f \Big|
\]
for all $f\in W^{2,p}(\Omega)\cap W^{1,p}_0(\Omega)$.
\end{lem}

\begin{proof}
Suppose that this is not true. Then there exists a sequence of functions $f_k\in W^{2,p}(\Omega)\cap W^{1,p}_0(\Omega)$ such that $\|f_k\|_{L^p(\Omega}                \newcommand{\pa}{\partial)}=1$ for all $k$, and
\[
\lim_{k\to\infty}\Big \|\Delta f_k+bf_k+\frac{n+2}{n-2} v_\infty^{\frac{4}{n-2}} f_k \Big\|_{L^p(\Omega}                \newcommand{\pa}{\partial)}+ \lim_{k\to\infty}\sup_{1\le l\le L} \Big|\int_\Omega}                \newcommand{\pa}{\partial v_\infty^{\frac{4}{n-2}} \phi_l f_k \Big|=0.
\]
If $p>1$, then by the $W^{2,p}$ estimates, we have $\|f_k\|_{W^{2,p}(\Omega}                \newcommand{\pa}{\partial)}\le C$. If $p=1$, by the estimates of Brezis-Strauss \cite{BrezisS}, $\|f_k\|_{W^{1,q}(\Omega}                \newcommand{\pa}{\partial)}\le C$ for some $q>1$. Therefore, by the compactness, we obtain an $f$ such that $\|f\|_{L^p(\Omega}                \newcommand{\pa}{\partial)}=1$,
$\int_\Omega}                \newcommand{\pa}{\partial v_\infty^{\frac{4}{n-2}} \phi_l f_k =0$ for all $1\le l\le L$, and
\[
\Delta f+bf+\frac{n+2}{n-2} v_\infty^{\frac{4}{n-2}} f=0
\]
in the distribution sense. Multiplying $\phi_l$ and integrating by parts, we have
\[
\left(\mu_l- \frac{n+2}{n-2} \right)\int_\Omega}                \newcommand{\pa}{\partial v_\infty^{\frac{4}{n-2}} \phi_l f_k =0.
\]
Hence, $\int_\Omega}                \newcommand{\pa}{\partial v_\infty^{\frac{4}{n-2}} \phi_l f_k =0$ for all $l> L$.
Meanwhile, from the elliptic regularity, we know that $f\in L^\infty(\Omega)$. Hence, $f\in\mathcal{L}^2(\Omega)$, and thus, $f\equiv 0$, which is a contradiction.
\end{proof}

\begin{lem} \label{lem:projectionest1}  There exists a constant $C$ such that
\begin{itemize}
\item[(i)]
\[
\|f\|_{L^{\frac{n+2}{n-2}}(\Omega}                \newcommand{\pa}{\partial)}\le C \Big \| \Pi(\Delta f+bf+\frac{n+2}{n-2} v_\infty^{\frac{4}{n-2}} f )\Big\|_{L^{\frac{n(n+2)}{n^2+4}}(\Omega}                \newcommand{\pa}{\partial)}+C \sup_{1\le l\le L} \Big|\int_\Omega}                \newcommand{\pa}{\partial v_\infty^{\frac{4}{n-2}} \phi_l f \Big|
\]
for all $f\in W^{2,\frac{n(n+2)}{n^2+4}}(\Omega)\cap W^{1,\frac{n(n+2)}{n^2+4}}_0(\Omega)$.
\item[(ii)]
\[
\|f\|_{L^{1}(\Omega}                \newcommand{\pa}{\partial)}\le C \Big \| \Pi(\Delta f+bf+\frac{n+2}{n-2} v_\infty^{\frac{4}{n-2}} f )\Big\|_{L^{1}(\Omega}                \newcommand{\pa}{\partial)}+C \sup_{1\le l\le L} \Big|\int_\Omega}                \newcommand{\pa}{\partial v_\infty^{\frac{4}{n-2}} \phi_l f \Big|
\]
for all $f\in W^{2,1}(\Omega)\cap W^{1,1}_0(\Omega)$.
\end{itemize}

\end{lem}

\begin{proof}
Given Lemma \ref{lem:Lpprojection}, the proof is the same as that of Lemma 6.3 in \cite{Br05}.
\end{proof}

\begin{lem} \label{lem:ift} There exists $\delta_1>0$ such that for every $z=(z_1,\dots,z_L)\in \mathbb{R}}              \newcommand{\Sn}{\mathbb{S}^n^L$ with $|z|\le \delta_1$, there exists $\xi_z \in C_0^{\frac{3n-2}{n-2}}(\overline \Omega}                \newcommand{\pa}{\partial)$ satisfying $1/2\le \xi_z/v_\infty\le 2$ in $\Omega}                \newcommand{\pa}{\partial$,
\[
\int_{\Omega}                \newcommand{\pa}{\partial} v_\infty^{\frac{4}{n-2}} (\xi_z- v_\infty) \phi_l\,\ud x= z_l, \quad l=1,\dots, L,
\]
and
\begin{equation}\label{eq:projection0}
\Pi(\Delta \xi_z+b\xi_z+ \xi_z ^{\frac{n+2}{n-2}})=0.
\end{equation}
Furthermore, the map $z\mapsto \xi_z$ is real analytic and $\frac{\pa }{\pa z_1}\xi_z(0)=v_\infty$, $\frac{\pa }{\pa z_l}\xi_z(0)=\phi_l$ for $2\le l\le L$.
\end{lem}

\begin{proof} Let $\xi_{z}= (1+z_1) v_\infty +\sum_{l=2}^L z_l \phi_l+h$, where $$h\in \mathcal{H}:=span\{\phi_1,\dots, \phi_L\}^{\perp}$$ and the ``$\perp$" is with respect to the inner product \eqref{eq:inner-product}. By a direct computation,
\begin{align*}
&\Pi(\Delta \xi_z+b\xi_z+ \xi_z ^{\frac{n+2}{n-2}}) \\&= (\Delta +b)\xi_z - \sum_{l=1}^L\Big( \int_\Omega (\Delta +b)\xi_z \phi_l\Big) v_\infty^{\frac{4}{n-2}} \phi_l
+ \xi_z^{\frac{n+2}{n-2}}- \sum_{l=1}^L\Big( \int_\Omega \xi_z^{\frac{n+2}{n-2}} \phi_l\Big) v_\infty^{\frac{4}{n-2}} \phi_l\\&
= (\Delta+b)h +\xi_z^{\frac{n+2}{n-2}}- \sum_{l=1}^L\Big( \int_\Omega \xi_z^{\frac{n+2}{n-2}} \phi_l\Big) v_\infty^{\frac{4}{n-2}} \phi_l=:G(z,h).
\end{align*}
For any $p>n$, we claim that there exists a small constant $\delta>0$ such that
\[
G: \{|z|<\delta \}\times \{\|h\|_{W^{2,p}(\Omega}                \newcommand{\pa}{\partial) \cap W^{1,p}_0 (\Omega}                \newcommand{\pa}{\partial)}<\delta\} \to \Pi( L^p(\Omega}                \newcommand{\pa}{\partial))
\]
is analytic. Indeed, let $\Phi(z,h)=\xi_{z}$, $\mathcal{L}u= \Delta u +bu+u^{\frac{n+2}{n-2}}$. Then we have $G=\Pi \circ \mathcal{L} \circ \Phi $. Obviously, the linear maps $\Phi$ and $\Pi$ are analytic. By Lemma 5.3 of Feireisl-Simondon \cite{FS}, $\mathcal{L}$ is also analytic in some small neighborhood of $v_\infty$ in $W^{2,p}(\Omega}                \newcommand{\pa}{\partial) \cap W^{1,p}_0 (\Omega}                \newcommand{\pa}{\partial)$.

Note that $G(0,0)=0$ and
\[
G_{h}(0,0)\varphi= (\Delta+b) \varphi +\frac{n+2}{n-2}\left( v_\infty^{\frac{4}{n-2}} \varphi -\sum_{l=1}^L \Big(\int_\Omega   v_\infty^{\frac{4}{n-2}} \varphi \phi_l\Big)  v_\infty^{\frac{4}{n-2}} \phi_l\right).
\]
Since $G_h(0,0)$ is coercive on $H^1_0(\Omega}                \newcommand{\pa}{\partial)\cap\mathcal{H}$, then $G_h(0,0): W^{2,p}(\Omega}                \newcommand{\pa}{\partial) \cap W^{1,p}_0 (\Omega}                \newcommand{\pa}{\partial)\cap\mathcal{H} \to  \Pi( L^p(\Omega}                \newcommand{\pa}{\partial))$ is invertible, and both $G_h(0,0)$ and $(G_h(0,0))^{-1}$ are continuous. By the Implicit Function Theorem we can find $h(z)\in W^{2,p}(\Omega}                \newcommand{\pa}{\partial) \cap W^{1,p}_0 (\Omega}                \newcommand{\pa}{\partial)\cap\mathcal{H}$ such that $G(z,h(z))=0$ and $h$ is analytic in $z$.  The regularity of $(h(z))(\cdot)$ follows from linear elliptic equations theory.  Since $h(0)=0$, and $0=G_z(0,0)+G_h(0,0)\pa_z h(0)= G_h(0,0)\pa_z h(0)$, we have $\pa_z h(0)=0$. It follows that $\frac{\pa }{\pa z_1}\xi_z(0)=v_\infty$ and $\frac{\pa }{\pa z_l}\xi_z(0)=\phi_l$ for $2\le l\le L$.  Therefore, we complete the proof.
\end{proof}

The difference of the energy at $\xi_z$ and $v_\infty$ can be controlled as follows.
\begin{lem} \label{lem:lojasiewiczgradient} There exists a real number $0<\gamma<1$ such that
\[
F(\xi_z)-F(v_\infty) \le 2 \sup_{1\le l\le L} \Big|  \int_{\Omega}                \newcommand{\pa}{\partial} (\Delta \xi_z+b\xi_z +\xi_z^{\frac{n+2}{n-2}})  \phi_l\,\ud x \Big|^{1+\gamma}
\]
if $z$ is sufficiently small.
\end{lem}

\begin{proof} By Lemma \ref{lem:ift} above and Lemma 5.3 of \cite{FS}, the function $z\mapsto F(\xi_z)$ is real analytic. Using the Lojasiewicz inequality, we have
\[
|F(\xi_z)-F(v_\infty)| \le \sup_{l} \Big|\frac{\pa }{\pa z_l} F(\xi_z)\Big|^{1+\gamma}
\]
if $z$ is sufficiently small. By a direct computation, we have
\begin{align*}
\frac{\pa }{\pa z_l} F(\xi_z) =-2 \int_{\Omega}                \newcommand{\pa}{\partial} (\Delta \xi_z+b\xi_z +\xi_z^{\frac{n+2}{n-2}})  \frac{\pa \xi_z }{\pa z_l}\,\ud x=-2  \int_{\Omega}                \newcommand{\pa}{\partial} (\Delta \xi_z+b\xi_z +\xi_z^{\frac{n+2}{n-2}})  \phi_l\,\ud x.
\end{align*}
Therefore, we complete the proof.
\end{proof}

For every $\nu$, as in the beginning of Section \ref{sec: zero}, let $\mathcal{A}_\nu$ be the closed set of all $m$-tuplets $(x_k,\lda_k, \alpha}                \newcommand{\lda}{\lambda_k)_{1\le k\le m}$ satisfying $(x_k,\lda_k, \alpha}                \newcommand{\lda}{\lambda_k) \in \overline B_{\frac{1}{\lda_{k,\nu}^*}}(x_{k,\nu}^*) \times [\frac{\lda_{k,\nu}^*}{2}, \frac{3\lda_{k,\nu}^*}{2} ] \times [\frac{1}{2},\frac32]$.  Let  $\delta_1>0$ be the constant in Lemma \ref{lem:ift} and $\overline B_{\delta_1}^L$ is the open ball in $\mathbb{R}}              \newcommand{\Sn}{\mathbb{S}^n^L$ centered at origin with radius $\delta_1$.
Choose an element  $(z_\nu, (x_{k,\nu},\lda_{k,\nu}, \alpha}                \newcommand{\lda}{\lambda_{k,\nu})_{1\le k\le m}) \in  \overline B_{\delta_1}^L\times \mathcal{A}_\nu$ such that
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:finite-minimum-2}
\left\|v_\nu- \xi_{z_\nu}-\sum_{k=0}^m \alpha}                \newcommand{\lda}{\lambda_{k,\nu} \xi_{x_{k,\nu}, \lda_{k,\nu}}\right\| =\inf_{(z, (x_k,\lda_k, \alpha}                \newcommand{\lda}{\lambda_k)_{1\le k\le m}) \in \mathcal{\bar A}_{\nu}} \left\|v_\nu- \xi_z-\sum_{k=0}^m \alpha}                \newcommand{\lda}{\lambda_{k} \xi_{x_{k}, \lda_{k}}\right\|.
\ee
Similar to \eqref{eq:bubbles-1} - \eqref{eq:bubbles-3}, we have
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:bubbles-1'}
\frac{ \lda_{i, \nu} }{\lda_{j, \nu}}+\frac{ \lda_{j, \nu} }{\lda_{i, \nu}}+\lda_{i, \nu} \lda_{j, \nu} |x_{i,\nu}-x_{j,\nu} |^2 \to \infty,
\ee
and for all $k$
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:bubbles-2'}
\lda_{k, \nu} d(x_{k,\nu}) \to \infty
\ee
as $\nu \to \infty$.
In addition, $d(x_{k,\nu})>\delta/2$, and
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:bubbles-3'}
\left\|v_\nu-\xi_{z_\nu}-\sum_{k=1}^m  \alpha}                \newcommand{\lda}{\lambda_k\xi_{x_{k,\nu}, \lda_{k,\nu}}\right\| \to 0
\ee
as $\nu\to \infty$.

By the triangle inequality,
\begin{align*}
&\left\| \xi_{z_\nu} -v_\infty+\sum_{k=1}^m  \alpha}                \newcommand{\lda}{\lambda_k\xi_{x_{k,\nu}, \lda_{k,\nu}}- \sum_{k=1}^m \xi_{x^*_{k,\nu}, \lda_{k,\nu}^*} \right\| \\& \le \left\|v_\nu- \xi_{z_\nu}-\sum_{k=1}^m  \alpha}                \newcommand{\lda}{\lambda_k\xi_{x_{k,\nu}, \lda_{k,\nu}} \right\|+\left\|v_\nu-v_\infty- \sum_{k=1}^m \xi_{x^*_{k,\nu}, \lda_{k,\nu}^*} \right\|=o(1).
\end{align*}
It follows that, for all $1\le k\le m$,
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:bubbles-4'}
|z_{\nu}|=o(1), \quad |x_{k,\nu}- x_{k,\nu}^{*}|=o(1) \frac{1}{\lda_{k,\nu}^*},\quad  \frac{\lda_{k,\nu}}{\lda_{k,\nu}^*}=1+o(1), \quad \alpha}                \newcommand{\lda}{\lambda_{k,\nu}=1+o(1).
\ee
In particular, $(z_\nu, (x_{k,\nu},\lda_{k,\nu}, \alpha}                \newcommand{\lda}{\lambda_{k,\nu})_{1\le k\le m}) \in  \overline B_{\delta_1}^L\times \mathcal{A}_\nu$ is an interior point.

In the sequel, we assume
\begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:order'}
\lda_{1,\nu}\ge \lda_{2,\nu} \ge \dots \ge \lda_{m,\nu}.
\ee
 Let
\begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:decomposition-1'}
U_\nu=\xi_{z_\nu}+ \sum_{k=1}^m  \alpha}                \newcommand{\lda}{\lambda_{k,\nu}\xi_{x_{k,\nu}, \lda_{k,\nu}}, \quad w_\nu=v_{\nu}-U_{\nu}.
\ee

\begin{lem} \label{lem:proj-2}  We have for $1\le l \le L$,
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:proj-2}
\Big|\int_\Omega}                \newcommand{\pa}{\partial v_\infty^{\frac{4}{n-2}} \phi_l w_\nu\,\ud x \Big|\le o(1) \int_{\Omega}                \newcommand{\pa}{\partial} |w_\nu|\,\ud x,
\ee
and for $1\le k\le m$,
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:proj-3}
\begin{split}
&\Big|\int_\Omega}                \newcommand{\pa}{\partial \xi_{(x_{k,\nu}, \lda_{k,\nu})}^{\frac{n+2}{n-2}}  w_\nu \,\ud x\Big|+\Big|\int_\Omega}                \newcommand{\pa}{\partial \xi_{(x_{k,\nu}, \lda_{k,\nu})}^{\frac{n+2}{n-2}}  \frac{1-\lda^2|x-x_{k,\nu}|^2}{1+\lda^2|x-x_{k,\nu}|^2} w_\nu \,\ud x\Big |\\[2mm] &\quad  +\Big|\int_\Omega}                \newcommand{\pa}{\partial \xi_{(x_{k,\nu}, \lda_{k,\nu})}^{\frac{n+2}{n-2}}  \frac{\lda^2(x-x_{k,\nu})}{1+\lda^2|x-x_{k,\nu}|^2} w_\nu \,\ud x\Big | \le o(1) \Big(\int_{\Omega}                \newcommand{\pa}{\partial} |w_\nu|^{\frac{2n}{n-2}}\,\ud x\Big)^{\frac{n-2}{2n}}.
\end{split}
\ee
\end{lem}
\begin{proof} Let $\tilde \phi_l=\frac{\pa }{\pa z_l} \xi_z$. By \eqref{eq:bubbles-4'}, we have $\|\tilde \phi_1 -v_\infty\|_{C^2(\Omega}                \newcommand{\pa}{\partial)}=o(1)$ and $\|\tilde \phi_l -\phi_l\|_{C^2(\Omega}                \newcommand{\pa}{\partial)}=o(1)$ for $l=2,\dots, L$.  By the definition of $(z_\nu, (x_{k,\nu},\lda_{k,\nu}, \alpha}                \newcommand{\lda}{\lambda_{k,\nu})_{1\le k\le m})$, we have
\[
\int \nabla \tilde \phi_l \nabla w_\nu- b\tilde \phi_l w_\nu=0.
\]
Hence,
\begin{align*}
\mu_l\int_\Omega v_\infty^\frac{4}{n-2} \phi_lw_\nu\,\ud x&=\int_\Omega \Big(-\Delta\phi_l-b\phi_l \Big)w_\nu\,\ud x\\
&=\int_\Omega \Big(\Delta(\tilde\phi_l- \phi_l)+b(\tilde\phi_l- \phi_l) \Big)w_\nu\,\ud x.
\end{align*}
Since $\mu_l>0$, then we can conclude  \eqref{eq:proj-2}. The proof of \eqref{eq:proj-3} is the same as that of Lemma \ref{lem:proj-1}.
\end{proof}

Now we can show the  non-degeneracy estimates of the second variation of $F$ for $w_\nu$.
\begin{lem} \label{lem:2nd-variational-2} For large $\nu$, we have
\[
\frac{n+2}{n-2} \int_\Omega}                \newcommand{\pa}{\partial \big(v_\infty^{\frac{4}{n-2}}+ \sum_{k=1}^m  \xi_{x_{k,\nu},\lda_{k,\nu}}  ^{\frac{4}{n-2}}\big) w_\nu^2   \le (1-c) \int_{\Omega}                \newcommand{\pa}{\partial} (|\nabla w_\nu|^2-bw_\nu^2)\,\ud x,
\]
where $c>0$ is independent of $\nu$.
\end{lem}

\begin{proof} We assume $w_\nu$ is not zero, otherwise there is nothing to prove. Define $\tilde w_\nu= \frac{w_\nu}{\|w_\nu\|}$. Suppose the lemma is not true.  Then we can find a subsequence of $\{\tilde w_\nu\}$ (still denoted by $\{\tilde w_\nu\}$) satisfying
\begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:coer-1b}
\lim_{\nu\to \infty}\frac{n+2}{n-2}  \int_\Omega}                \newcommand{\pa}{\partial  \big(v_\infty^{\frac{4}{n-2}}+ \sum_{k=1}^m  \xi_{x_{k,\nu},\lda_{k,\nu}}  ^{\frac{4}{n-2}}\big)  \tilde w_\nu^2  \ge 1.
\ee
By \eqref{eq:sobolev-b},
\begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:coer-2b}
\int_{\Omega}                \newcommand{\pa}{\partial} |\tilde w_\nu|^{\frac{2n}{n-2}} \le K_b^{\frac{n}{n-2}} \|\tilde w_\nu\|= K_b^{\frac{n}{n-2}}.
\ee
By \eqref{eq:bubbles-1'} and \eqref{eq:order'}, we can find $R_\nu\to \infty$, $R_\nu \lda_{j,\nu}^{-1} \to 0 $ for all $1\le j\le m$, and
\begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:coer-3b}
\frac{\lda_{i,\nu}}{R_\nu} (\lda_{j,\nu}^{-1} +|x_{i,\nu}-x_{j,\nu}|) \to \infty
\ee
for all $i<j$. Set
\[
\Omega}                \newcommand{\pa}{\partial_{j,\nu}= B_{R_\nu \lda_{j,\nu}^{-1}}(x_{j,\nu})\setminus \bigcup_{i=1}^{j-1} B_{R_\nu \lda_{i,\nu}^{-1}}(x_{i,\nu}).
\]
By \eqref{eq:coer-1b} and $\|\tilde w_\nu\|=1$, there are two cases:
\begin{itemize}
\item[(i).] We can find $1\le j\le m$ such that
\[
\lim_{\nu\to \infty } \int_\Omega}                \newcommand{\pa}{\partial    \xi_{x_{j,\nu},\lda_{j,\nu}}  ^{\frac{4}{n-2}} \tilde w_\nu^2 >0
\]
and
\[
\lim_{\nu\to \infty  } \int_{\Omega}                \newcommand{\pa}{\partial_{j,\nu}} (|\nabla \tilde w_\nu|^2- b\tilde w_\nu^2 ) \le \frac{n+2}{n-2}  \int_\Omega}                \newcommand{\pa}{\partial    \xi_{x_{j,\nu},\lda_{j,\nu}}  ^{\frac{4}{n-2}} \tilde w_\nu^2.
\]
\item[(ii).]
\[
\lim_{\nu\to \infty } \int_\Omega}                \newcommand{\pa}{\partial    v_\infty^{\frac{4}{n-2}} \tilde w_\nu^2 >0
\]
and
\[
\lim_{\nu\to \infty  } \int_{\Omega}                \newcommand{\pa}{\partial\setminus \cup_{j}\Omega}                \newcommand{\pa}{\partial_{j,\nu}} (|\nabla \tilde w_\nu|^2- b\tilde w_\nu^2 ) \le \frac{n+2}{n-2}  \int_\Omega}                \newcommand{\pa}{\partial    v_\infty  ^{\frac{4}{n-2}} \tilde w_\nu^2.
\]
\end{itemize}

In the former case, we can obtain a contradiction similar to that in the proof of Lemma \ref{lem:2nd-variational}.

In the later case, after passing to subsequence we suppose $ \tilde w_\nu \rightharpoonup \tilde w$ in $H_0^1$ as $\nu\to \infty$. It follows that
\begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:non-trivial}
\int_\Omega}                \newcommand{\pa}{\partial    v_\infty^{\frac{4}{n-2}} \tilde w^2 >0
\ee
and
\begin{equation}}      \newcommand{\ee}{\end{equation} \label{eq:lower-eigen-space}
 \int_{\Omega}                \newcommand{\pa}{\partial} (|\nabla \tilde w|^2- b\tilde w^2 ) \le \frac{n+2}{n-2}  \int_\Omega}                \newcommand{\pa}{\partial    v_\infty  ^{\frac{4}{n-2}} \tilde w^2.
\ee
By \eqref{eq:proj-2}, we further have
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:o-1}
\int_\Omega}                \newcommand{\pa}{\partial v_\infty^{\frac{4}{n-2}} \tilde w \phi_l=0 \quad \mbox{for }l=1,\dots, L.
\ee
Combing \eqref{eq:lower-eigen-space} and \eqref{eq:o-1}, $\tilde w$ has to be identically zero, which contradicts to \eqref{eq:non-trivial}.

Therefore, Lemma \ref{lem:2nd-variational-2} is proved.
\end{proof}

\begin{cor}\label{cor:coercive2}
For large $\nu$, we have
\[
\frac{n+2}{n-2} \int_\Omega}                \newcommand{\pa}{\partial U_\nu ^{\frac{4}{n-2}} w_\nu^2   \le (1-c) \int_{\Omega}                \newcommand{\pa}{\partial} (|\nabla w_\nu|^2-bw_\nu^2)\,\ud x,
\]
where $c>0$ is independent of $\nu$.
\end{cor}
\begin{proof}
It follows from Lemma \ref{lem:2nd-variational-2}, H\"older's inequality, the Sobolev inequality \eqref{eq:sobolev-b} and the fact that
\[
\int_\Omega \Big|U_\nu ^{\frac{4}{n-2}} -v_\infty^{\frac{4}{n-2}}- \sum_{k=1}^m  \xi_{x_{k,\nu},\lda_{k,\nu}}  ^{\frac{4}{n-2}}\Big|^{\frac{n}{2}}=o(1).
\]
\end{proof}

The following two lemmas are estimates of $v_\nu-\xi_{z_\nu}$ in $L^{\frac{n+2}{n-2}}(\Omega}                \newcommand{\pa}{\partial)$ and $L^{1}(\Omega}                \newcommand{\pa}{\partial)$, respectively.
\begin{lem}\label{lem:n+2norm}
For large $\nu$, we have
\[
\|v_\nu-\xi_{z_\nu}\|_{L^{\frac{n+2}{n-2}}(\Omega}                \newcommand{\pa}{\partial)}^{\frac{n+2}{n-2}} \le C \|v_\nu^{\frac{n+2}{n-2}}(\mathcal{R}(t_\nu)-1)\|_{L^{\frac{2n}{n+2}}(\Omega}                \newcommand{\pa}{\partial)}^{\frac{n+2}{n-2}}+C  \sum_{k=1}^m \lda_{k,\nu}^{\frac{2-n}{2}},
\]
where $C>0$ is independent of $\nu$.
\end{lem}
\begin{proof}
From \eqref{eq:Q}, we have
\[
\Delta v_\nu+bv_\nu+v_\nu^{\frac{n+2}{n-2}}=(1-\mathcal{R}(t_\nu)) v_\nu^{\frac{n+2}{n-2}}.
\]
Combining with \eqref{eq:projection0}, we obtain
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:projcorrection}
\begin{split}
&\Pi\Big(\Delta (v_\nu-\xi_{z_\nu})+b(v_\nu-\xi_{z_\nu})+ \frac{n+2}{n-2} v_\infty^\frac{4}{n-2}(v_\nu-\xi_{z_\nu})\Big)\\
&=\Pi\Big((1-\mathcal{R}(t_\nu))v_\nu^{\frac{n+2}{n-2}} -\frac{n+2}{n-2} (\xi_{z_\nu}^{\frac{4}{n-2}}-v_\infty^\frac{4}{n-2})(v_\nu-\xi_{z_\nu})\\
&\quad\quad\quad+ \xi_{z_\nu}^{\frac{n+2}{n-2}}+ \frac{n+2}{n-2}\xi_{z_\nu}^{\frac{4}{n-2}}(v_\nu-\xi_{z_\nu})-v_\nu^{\frac{n+2}{n-2}}\Big).
\end{split}
\ee
Apply (i) of Lemma \ref{lem:projectionest1}  to $v_\nu-\xi_{z_\nu}$, we obtain
\begin{align*}
&\left\|v_\nu-\xi_{z_\nu}\right\|_{L^\frac{n+2}{n-2}(\Omega)}\\
&\le C\left\|(1-\mathcal{R}(t_\nu))v_\nu^{\frac{n+2}{n-2}}\right\|_{L^\frac{n(n+2)}{n^2+4}(\Omega)} +C\left\|(\xi_{z_\nu}^{\frac{4}{n-2}}-v_\infty^\frac{4}{n-2})(v_\nu-\xi_{z_\nu})\right\|_{L^\frac{n(n+2)}{n^2+4}(\Omega)}\\
&\quad+ C\left\|\xi_{z_\nu}^{\frac{n+2}{n-2}}+ \frac{n+2}{n-2}\xi_{z_\nu}^{\frac{4}{n-2}}(v_\nu-\xi_{z_\nu})-v_\nu^{\frac{n+2}{n-2}}\right\|_{L^\frac{n(n+2)}{n^2+4}(\Omega)}+C\sup_{1\le l\le L}\left|\int_\Omega}                \newcommand{\pa}{\partial v_\infty^{\frac{4}{n-2}} \phi_l (v_\nu-\xi_{z_\nu})\right|.
\end{align*}
Using the estimates for all $a,b\ge 0$ that
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:pointwise3}
\left|a^{\frac{n+2}{n-2}} +\frac{n+2}{n-2} a^\frac{4}{n-2}(b-a)-b^\frac{n+2}{n-2}\right|\le C a^{\max(0,\frac{4}{n-2}-1)}|b-a|^{\min(\frac{n+2}{n-2},2)}+ C |b-a|^\frac{n+2}{n-2},
\ee
we obtain
\begin{align*}
&\left\|\xi_{z_\nu}^{\frac{n+2}{n-2}}+ \frac{n+2}{n-2}\xi_{z_\nu}^{\frac{4}{n-2}}(v_\nu-\xi_{z_\nu})-v_\nu^{\frac{n+2}{n-2}}\right\|_{L^\frac{n(n+2)}{n^2+4}(\Omega)}\\
&\le C \left\| |v_\nu-\xi_{z_\nu}|^{\min(\frac{n+2}{n-2},2)}+  |v_\nu-\xi_{z_\nu}|^\frac{n+2}{n-2}\right\|_{L^\frac{n(n+2)}{n^2+4}(\Omega)}\\
&\le C \left\| |v_\nu-\xi_{z_\nu}|^{\min(\frac{n+2}{n-2},2)}+  |v_\nu-\xi_{z_\nu}|^\frac{n+2}{n-2}\right\|_{L^\frac{n(n+2)}{n^2+4}(\cup_{k=1}^m B_{N/\lambda_{k,\nu}}(x_{k,\nu}))}\\
&\quad +C \left\| |v_\nu-\xi_{z_\nu}|^{\min(\frac{n+2}{n-2},2)}+  |v_\nu-\xi_{z_\nu}|^\frac{n+2}{n-2}\right\|_{L^\frac{n(n+2)}{n^2+4}(\Omega\setminus \cup_{k=1}^m B_{N/\lambda_{k,\nu}}(x_{k,\nu}))},
\end{align*}
where $N$ is a large real number to be chosen later. Using H\"older's inequality, we have
\begin{align*}
&\left\| |v_\nu-\xi_{z_\nu}|^{\min(\frac{n+2}{n-2},2)}+  |v_\nu-\xi_{z_\nu}|^\frac{n+2}{n-2}\right\|_{L^\frac{n(n+2)}{n^2+4}(\cup_{k=1}^m B_{N/\lambda_{k,\nu}}(x_{k,\nu}))}\\
&\le C \sum_{k=1}^m (N/\lambda_{k,\nu})^\frac{(n-2)^2}{2(n+2)} \left\| |v_\nu-\xi_{z_\nu}|^{\min(\frac{n+2}{n-2},2)}+  |v_\nu-\xi_{z_\nu}|^\frac{n+2}{n-2}\right\|_{L^\frac{2n}{n+2}(\Omega)}\\
& \le C \sum_{k=1}^m (N/\lambda_{k,\nu})^\frac{(n-2)^2}{2(n+2)}
\end{align*}
and
\begin{align*}
&\left\| |v_\nu-\xi_{z_\nu}|^{\min(\frac{n+2}{n-2},2)}+  |v_\nu-\xi_{z_\nu}|^\frac{n+2}{n-2}\right\|_{L^\frac{n(n+2)}{n^2+4}(\Omega\setminus \cup_{k=1}^m B_{N/\lambda_{k,\nu}}(x_{k,\nu}))}\\
&\le \left\| |v_\nu-\xi_{z_\nu}|^{\min(\frac{4}{n-2},1)}+  |v_\nu-\xi_{z_\nu}|^\frac{4}{n-2}\right\|_{L^\frac{n}{2}(\Omega\setminus \cup_{k=1}^m B_{N/\lambda_{k,\nu}}(x_{k,\nu}))}\cdot  \left\| v_\nu-\xi_{z_\nu}\right\|_{L^\frac{n+2}{n-2}(\Omega)}.
\end{align*}
Since
\begin{align*}
&\left\| v_\nu-\xi_{z_\nu}\right\|_{L^\frac{2n}{n-2}(\Omega\setminus \cup_{k=1}^m B_{N/\lambda_{k,\nu}}(x_{k,\nu}))}\\
&=  \left\|  \sum_{k=1}^m  \alpha}                \newcommand{\lda}{\lambda_{k,\nu}\xi_{x_{k,\nu}, \lda_{k,\nu}}+ w_\nu   \right\|_{L^\frac{2n}{n-2}(\Omega\setminus \cup_{k=1}^m B_{N/\lambda_{k,\nu}}(x_{k,\nu}))}\\
&\le \sum_{k=1}^m  \alpha}                \newcommand{\lda}{\lambda_{k,\nu}\left\|  \xi_{x_{k,\nu}, \lda_{k,\nu}} \right\|_{L^\frac{2n}{n-2}(\Omega\setminus B_{N/\lambda_{k,\nu}}(x_{k,\nu}))} + \left\|  w_\nu  \right\|_{L^\frac{2n}{n-2}(\Omega)} \\
&\le C N^{-\frac{n-2}{2}}+o(1),
\end{align*}
we have
\begin{align*}
&\left\|\xi_{z_\nu}^{\frac{n+2}{n-2}}+ \frac{n+2}{n-2}\xi_{z_\nu}^{\frac{4}{n-2}}(v_\nu-\xi_{z_\nu})-v_\nu^{\frac{n+2}{n-2}}\right\|_{L^\frac{n(n+2)}{n^2+4}(\Omega)}\\
&\le C \sum_{k=1}^m (N/\lambda_{k,\nu})^\frac{(n-2)^2}{2(n+2)} +C (N^{-\frac{n-2}{2}}+N^{-2}+o(1))\left\| v_\nu-\xi_{z_\nu}\right\|_{L^\frac{n+2}{n-2}(\Omega)}.
\end{align*}
Also,
\begin{align}
\sup_{1\le l\le L}\left|\int_\Omega}                \newcommand{\pa}{\partial v_\infty^{\frac{4}{n-2}} \phi_l (v_\nu-\xi_{z_\nu})\right|
&=\sup_{1\le l\le L}\left|\int_\Omega}                \newcommand{\pa}{\partial v_\infty^{\frac{4}{n-2}} \phi_l \Big( \sum_{k=1}^m  \alpha}                \newcommand{\lda}{\lambda_{k,\nu}\xi_{x_{k,\nu}, \lda_{k,\nu}}+ w_\nu \Big)\right|\nonumber\\
& \le C \sum_{k=1}^m  \lda_{k,\nu}^{\frac{2-n}{2}}+ o(1)\|w_\nu\|_{L^1(\Omega)}\nonumber\\
&\le C \sum_{k=1}^m  \lda_{k,\nu}^{\frac{2-n}{2}}+o(1)\left\|v_\nu -\xi_{z_\nu}-  \sum_{k=1}^m  \alpha}                \newcommand{\lda}{\lambda_{k,\nu}\xi_{x_{k,\nu}, \lda_{k,\nu}}\right\|_{L^1(\Omega)}\nonumber\\
&\le C \sum_{k=1}^m  \lda_{k,\nu}^{\frac{2-n}{2}}+o(1)\left\|v_\nu -\xi_{z_\nu}\right\|_{L^1(\Omega)}\label{eq:projectiondifference}.
\end{align}
Putting these facts together, we have
\begin{align*}
&\left\|v_\nu-\xi_{z_\nu}\right\|_{L^\frac{n+2}{n-2}(\Omega)}\\
&\le C\left\|(1-\mathcal{R}(t_\nu))v_\nu^{\frac{n+2}{n-2}}\right\|_{L^\frac{2n}{n+2}(\Omega)}+C \sum_{k=1}^m (N/\lambda_{k,\nu})^\frac{(n-2)^2}{2(n+2)}  \\
&+C (N^{-\frac{n-2}{2}}+N^{-2}+o(1))\left\| v_\nu-\xi_{z_\nu}\right\|_{L^\frac{n+2}{n-2}(\Omega)}+C \sum_{k=1}^m  \lda_{k,\nu}^{\frac{2-n}{2}}.
\end{align*}
By choosing $N$ sufficiently large, we obtain
\begin{align*}
\left\|v_\nu-\xi_{z_\nu}\right\|_{L^\frac{n+2}{n-2}(\Omega)}\le C\left\|(1-\mathcal{R}(t_\nu))v_\nu^{\frac{n+2}{n-2}}\right\|_{L^\frac{2n}{n+2}(\Omega)}+C \sum_{k=1}^m \lambda_{k,\nu}^{-\frac{(n-2)^2}{2(n+2)}},
\end{align*}
from which the conclusion follows.
\end{proof}

\begin{lem}\label{lem:1norm}
For large $\nu$, we have
\[
\|v_\nu-\xi_{z_\nu}\|_{L^{1}(\Omega}                \newcommand{\pa}{\partial)}\le C \|v_\nu^{\frac{n+2}{n-2}}(\mathcal{R}(t_\nu)-1)\|_{L^{\frac{2n}{n+2}}(\Omega}                \newcommand{\pa}{\partial)}^{\frac{n+2}{n-2}}+C  \sum_{k=1}^m \lda_{k,\nu}^{\frac{2-n}{2}},
\]
where $C>0$ is independent of $\nu$.
\end{lem}
\begin{proof}
Using \eqref{eq:projcorrection}, and applying (ii) of Lemma \ref{lem:projectionest1}  to $v_\nu-\xi_{z_\nu}$, we obtain
\begin{align*}
&\left\|v_\nu-\xi_{z_\nu}\right\|_{L^1(\Omega)}\\
&\le C\left\|(1-\mathcal{R}(t_\nu))v_\nu^{\frac{n+2}{n-2}}\right\|_{L^1(\Omega)} +C\left\|(\xi_{z_\nu}^{\frac{4}{n-2}}-v_\infty^\frac{4}{n-2})(v_\nu-\xi_{z_\nu})\right\|_{L^1(\Omega)}\\
&\quad+ C\left\|\xi_{z_\nu}^{\frac{n+2}{n-2}}+ \frac{n+2}{n-2}\xi_{z_\nu}^{\frac{4}{n-2}}(v_\nu-\xi_{z_\nu})-v_\nu^{\frac{n+2}{n-2}}\right\|_{L^1(\Omega)}+C\sup_{1\le l\le L}\left|\int_\Omega}                \newcommand{\pa}{\partial v_\infty^{\frac{4}{n-2}} \phi_l (v_\nu-\xi_{z_\nu})\right|.
\end{align*}
It follows from \eqref{eq:pointwise3} that
\begin{align*}
&\left\|\xi_{z_\nu}^{\frac{n+2}{n-2}}+ \frac{n+2}{n-2}\xi_{z_\nu}^{\frac{4}{n-2}}(v_\nu-\xi_{z_\nu})-v_\nu^{\frac{n+2}{n-2}}\right\|_{L^1(\Omega)}\\
&\le C \left\| |v_\nu-\xi_{z_\nu}|^{\min(\frac{n+2}{n-2},2)}+  |v_\nu-\xi_{z_\nu}|^\frac{n+2}{n-2}\right\|_{L^1(\Omega)}\\
&\le C \left\| v_\nu-\xi_{z_\nu}\right\|_{L^1(\Omega)}^{\max(0,1-\frac{n-2}{4})} \left\||v_\nu-\xi_{z_\nu}|^\frac{n+2}{n-2}\right\|_{L^1(\Omega)}^{\min(1,\frac{n-2}{4})} + C \left\| v_\nu-\xi_{z_\nu}\right\|_{L^{\frac{n+2}{n-2}}(\Omega)}^{\frac{n+2}{n-2}}\\
& \le  C \left\| v_\nu-\xi_{z_\nu}\right\|_{L^1(\Omega)}^{\max(0,1-\frac{n-2}{4})} \left\|v_\nu-\xi_{z_\nu}\right\|_{L^\frac{n+2}{n-2}(\Omega)}^{\frac{n+2}{n-2}\min(1,\frac{n-2}{4})} + C \left\| v_\nu-\xi_{z_\nu}\right\|_{L^{\frac{n+2}{n-2}}(\Omega)}^{\frac{n+2}{n-2}}\\
&\le  \frac{1}{2C} \left\| v_\nu-\xi_{z_\nu}\right\|_{L^1(\Omega)}+ C \left\| v_\nu-\xi_{z_\nu}\right\|_{L^{\frac{n+2}{n-2}}(\Omega)}^{\frac{n+2}{n-2}},
\end{align*}
where we used H\"older's inequality in the second inequality and the Young inequality in the last inequality. Combining \eqref{eq:projectiondifference}, we have
\begin{align*}
&\left\|v_\nu-\xi_{z_\nu}\right\|_{L^1(\Omega)}\\
&\le C\left\|(1-\mathcal{R}(t_\nu))v_\nu^{\frac{n+2}{n-2}}\right\|_{L^1(\Omega)}+o(1)\left\| v_\nu-\xi_{z_\nu}\right\|_{L^1(\Omega)}  \\
&+ \frac{1}{2} \left\| v_\nu-\xi_{z_\nu}\right\|_{L^1(\Omega)}+ C \left\| v_\nu-\xi_{z_\nu}\right\|_{L^{\frac{n+2}{n-2}}(\Omega)}^{\frac{n+2}{n-2}} +C \sum_{k=1}^m  \lda_{k,\nu}^{\frac{2-n}{2}}.
\end{align*}
Then the conclusion follows from Lemma \ref{lem:n+2norm}.
\end{proof}

Using the above two lemmas, we can continue to estimate $F(\xi_{z_\nu})- F(v_\infty)$ from  Lemma \ref{lem:lojasiewiczgradient}.

\begin{prop} \label{prop:regular-part}
For all large $\nu$, we have
\[
F(\xi_{z_\nu})- F(v_\infty) \le  C\left (\int_\Omega}                \newcommand{\pa}{\partial |\mathcal{R}(t_\nu)-1|^{\frac{2n}{n+2}} v_\nu^{\frac{2n}{n-2}}\right)^{\frac{n+2}{2n} (1+\gamma)} +C \sum_{k=1}^m \lda_{k,\nu}^{\frac{2-n}{2} (1+\gamma)},
\]
where $\gamma\in (0,1)$ is the one in Lemma \ref{lem:lojasiewiczgradient}.
\end{prop}
\begin{proof}
It follows from integration by parts that
\begin{align*}
& \int_{\Omega}                \newcommand{\pa}{\partial} (\Delta \xi_{z_\nu}+b\xi_{z_\nu} +\xi_{z_\nu}^{\frac{n+2}{n-2}})  \phi_l\,\ud x \\
 & =  \int_{\Omega}                \newcommand{\pa}{\partial} (\Delta v_\nu+bv_\nu +v_\nu^{\frac{n+2}{n-2}})  \phi_l\,\ud x + \mu_l \int_\Omega}                \newcommand{\pa}{\partial v_\infty^\frac{4}{n-2}\phi_l (v_\nu- \xi_{z_\nu})\,\ud x - \int_\Omega}                \newcommand{\pa}{\partial \phi_l (v_\nu^{\frac{n+2}{n-2}}- \xi_{z_\nu}^{\frac{n+2}{n-2}})\,\ud x\\
 & =\int_{\Omega}                \newcommand{\pa}{\partial} (1-\mathcal{R}(t_\nu))v_\nu^{\frac{n+2}{n-2}} \phi_l\,\ud x + \mu_l \int_\Omega}                \newcommand{\pa}{\partial v_\infty^\frac{4}{n-2}\phi_l (v_\nu- \xi_{z_\nu})\,\ud x - \int_\Omega}                \newcommand{\pa}{\partial \phi_l (v_\nu^{\frac{n+2}{n-2}}- \xi_{z_\nu}^{\frac{n+2}{n-2}})\,\ud x.
\end{align*}
Using the pointwise estimate
\[
|v_\nu^{\frac{n+2}{n-2}}- \xi_{z_\nu}^{\frac{n+2}{n-2}}| \le C \xi_{z_\nu}^{\frac{4}{n-2}} |v_\nu- \xi_{z_\nu}|+ C |v_\nu- \xi_{z_\nu}|^{\frac{n+2}{n-2}},
\]
we have
\begin{align*}
&\sup_{1\le l\le L} \Big|  \int_{\Omega}                \newcommand{\pa}{\partial} (\Delta \xi_z+b\xi_z +\xi_z^{\frac{n+2}{n-2}})  \phi_l\,\ud x \Big|\\
&\le C\left\|(1-\mathcal{R}(t_\nu))v_\nu^{\frac{n+2}{n-2}}\right\|_{L^{\frac{2n}{n+2}}(\Omega)} + C\left\| v_\nu-\xi_{z_\nu}\right\|_{L^1(\Omega)}+ C\left\| v_\nu-\xi_{z_\nu}\right\|_{L^\frac{n+2}{n-2}(\Omega)}^\frac{n+2}{n-2}\\
&\le C\left\|(1-\mathcal{R}(t_\nu))v_\nu^{\frac{n+2}{n-2}}\right\|_{L^{\frac{2n}{n+2}}(\Omega)} +C \left\|v_\nu^{\frac{n+2}{n-2}}(\mathcal{R}(t_\nu)-1)\right\|_{L^{\frac{2n}{n+2}}(\Omega}                \newcommand{\pa}{\partial)}^{\frac{n+2}{n-2}}+C  \sum_{k=1}^m \lda_{k,\nu}^{\frac{2-n}{2}}\\
&\le C\left\|(1-\mathcal{R}(t_\nu))v_\nu^{\frac{n+2}{n-2}}\right\|_{L^{\frac{2n}{n+2}}(\Omega)} +C  \sum_{k=1}^m \lda_{k,\nu}^{\frac{2-n}{2}},
\end{align*}
where we used Lemma \ref{lem:n+2norm} and Lemma \ref{lem:1norm} in the second inequality, and Proposition \ref{prop:Mq} in the last inequality.

Then the conclusion follows from Lemma \ref{lem:lojasiewiczgradient}.
\end{proof}

\begin{cor} \label{cor:regular+singular} If $n\ge 4$ and $b>0$ satisfying \eqref{eq:b}, we have
\[
F(U_\nu) \le F(v_\infty) +  \frac{2m}{n} Y(\mathbb{S}^n)^{n/2} + C \left(\int_\Omega}                \newcommand{\pa}{\partial |\mathcal{R}(t_\nu)-1|^{\frac{2n}{n+2}} v_\nu^{\frac{2n}{n-2}}\right)^{\frac{n+2}{2n} (1+\gamma)}.
\]
\end{cor}

\begin{proof} Let $\widetilde U_{\nu}= \sum_{k=1}^m \alpha}                \newcommand{\lda}{\lambda_k \xi_{x_{k,\nu}, \lda_{k,\nu}} $.
\begin{align*}
F(U_\nu)&= \int_\Omega}                \newcommand{\pa}{\partial |\nabla( \xi_{z_\nu}+\widetilde U_{\nu})|^2- b(\xi_{z_\nu}+\widetilde U_{\nu})^2 -\frac{n-2}{n} \int_\Omega}                \newcommand{\pa}{\partial (\xi_{z_\nu}+\widetilde U_{\nu})^{\frac{2n}{n-2}} \\
&= F(\xi_{z_\nu})+ F(\widetilde U_{\nu}) +2  \int_\Omega}                \newcommand{\pa}{\partial (\nabla  \xi_{z_\nu} \nabla \widetilde U_{\nu} -b  \xi_{z_\nu}  \widetilde U_{\nu} - \xi_{z_\nu} ^{\frac{n+2}{n-2}}  \widetilde U_{\nu})\\
&\quad- \frac{n-2}{n}\int_\Omega}                \newcommand{\pa}{\partial\left((\xi_{z_\nu}+ \widetilde U_{\nu})^{\frac{2n}{n-2}}-\frac{2n}{n-2}\xi_{z_\nu}^{\frac{n+2}{n-2}}\widetilde U_{\nu}- \xi_{z_\nu}^{\frac{2n}{n-2}}-\widetilde U^{\frac{2n}{n-2}} \right)
\end{align*}
We have
\begin{align*}
& \left|\int_\Omega}                \newcommand{\pa}{\partial (\nabla  \xi_{z_\nu} \nabla \widetilde U_{\nu} -b  \xi_{z_\nu}  \widetilde U_{\nu} - \xi_{z_\nu} ^{\frac{n+2}{n-2}}  \widetilde U_{\nu})\right|\\
 &= \left|\int_\Omega}                \newcommand{\pa}{\partial \Big(\Delta(\xi_{z_\nu}-v_\infty)  +b  (\xi_{z_\nu}-v_\infty)  + \xi_{z_\nu} ^{\frac{n+2}{n-2}} - v_\infty^{\frac{n+2}{n-2}} \Big) \widetilde U_{\nu})\right|\\
 & \le o(1)  \sum_{k=1}^m \lda_k^{\frac{2-n}{2}}.
 \end{align*}
By Lemma \ref{lem:basic-1}, there exists $c>0$, depending only on $n$ such that
\begin{align*}
(\xi_{z_\nu}+ \widetilde U_{\nu})^{\frac{2n}{n-2}}-\frac{2n}{n-2}\xi_{z_\nu}^{\frac{n+2}{n-2}}\widetilde U_{\nu}- \xi_{z_\nu}^{\frac{2n}{n-2}}-\widetilde U^{\frac{2n}{n-2}}\ge
 \begin{cases}  c\xi_{z_\nu}^{\frac{4}{n-2}}(\widetilde U_{\nu})^2 ,& \quad \mbox{if } \xi_{z_\nu}
 \ge U_{\nu}', \\
 c\xi_{z_\nu}(\widetilde U_{\nu})^{\frac{n+2}{n-2}} , &\quad \mbox{if } \xi_{z_\nu} < U_{\nu}'.
 \end{cases}
\end{align*}
Since $v_\infty/2\le \xi_{z_\nu}\le 2 v_\infty$, and
\[
\int_{|x|<\sqrt{\lda^{-1}}}  \left(\frac{\lda}{1+\lda^2 |x|^2}\right)^{\frac{n+2}{2}} \ge \lda^{-\frac{n-2}{2}} \int_{|y|<1} \left(\frac{1}{1+ |y|^2}\right)^{\frac{n+2}{2}} \ge c \lda^{-\frac{n-2}{2}},
\]
we have
\[
\int_\Omega}                \newcommand{\pa}{\partial\left((\xi_{z_\nu}+ \widetilde U_{\nu})^{\frac{2n}{n-2}}-\frac{2n}{n-2}\xi_{z_\nu}^{\frac{n+2}{n-2}}\widetilde U_{\nu}- \xi_{z_\nu}^{\frac{2n}{n-2}}-\widetilde U^{\frac{2n}{n-2}} \right)\ge c \sum_{k=1}^m \lda_k^{\frac{2-n}{2}}.
\]
Then, the conclusion follows from Proposition \ref{prop:regular-part} and Corollary \ref{cor:pure-singular}.
\end{proof}

\section{Convergence}\label{sec:convergence}

Using the  estimates in Corollaries \ref{cor:coercive} and \ref{cor:coercive2}, we have the following estimate of $F(v_{\nu}) -F_\infty $ for any sequence of times $\{t_\nu: \nu\in\mathbb{N}\}$.
\begin{prop}\label{prob:subsequenceestimate}
Let $n\ge 4$, and $b>0$ satisfy \eqref{eq:b}. Let $\{t_\nu: \nu\in\mathbb{N}\}$ be a sequence of times such that $t_\nu\to \infty$ as $\nu\to\infty$. Then, we can find a real number $\gamma\in (0,1)$ and a
constant $C>0$ such that, after passing to a subsequence, we have
\[
F(v_{\nu}) -F_\infty \le  C \left(\int_\Omega}                \newcommand{\pa}{\partial |\mathcal{R}(t_\nu)-1|^{\frac{2n}{n+2}} v_\nu^{\frac{2n}{n-2}}\right)^{\frac{n+2}{2n} (1+\gamma)}
\]
for all integers $\nu$ in that subsequence, where $F_\infty$ is the one defined in \eqref{eq:f_infty}. Note that $\gamma$ and $C$ may depend
on the sequence $\{t_\nu: \nu\in\mathbb{N}\}$.
\end{prop}

\begin{proof}
It follows from \eqref{eq:f_infty} that $F(v_\infty)=F_\infty$. Recall that $U_\nu=v_\nu-w_\nu$.  We have
\begin{align*}
F(v_{\nu}) -F(U_\nu)&= \frac{n-2}{n}\int_\Omega}                \newcommand{\pa}{\partial (U_\nu^{\frac{2n}{n-2}}-v_\nu^{\frac{2n}{n-2}}) +2 \int_\Omega}                \newcommand{\pa}{\partial (\nabla v_\nu \nabla w_\nu -b v_\nu w_\nu)    -\int_\Omega}                \newcommand{\pa}{\partial (| \nabla w_\nu|^2 -b  w_\nu^2)  \\&
= \frac{n-2}{n}\int_\Omega}                \newcommand{\pa}{\partial (U_\nu^{\frac{2n}{n-2}}-v_\nu^{\frac{2n}{n-2}}) +2 \int_\Omega}                \newcommand{\pa}{\partial \mathcal{R} v_\nu^{\frac{n+2}{n-2}} w_\nu    -\int_\Omega}                \newcommand{\pa}{\partial (| \nabla w_\nu|^2 -b  w_\nu^2)\\&=
\frac{n-2}{n}\int_\Omega}                \newcommand{\pa}{\partial \left(U_\nu^{\frac{2n}{n-2}}-v_\nu^{\frac{2n}{n-2}} +\frac{2n}{n-2} v_\nu^{\frac{n+2}{n-2}} w_\nu -\frac{n(n+2)}{(n-2)^2} U_\nu^{\frac{4}{n-2}} w_\nu^2\right) \\& \quad  +2 \int_\Omega}                \newcommand{\pa}{\partial (\mathcal{R}-1)v_\nu^{\frac{n+2}{n-2}} w_\nu -\int_\Omega}                \newcommand{\pa}{\partial \left(| \nabla w_\nu|^2 -b  w_\nu^2 -\frac{n+2}{n-2} U_\nu^{\frac{4}{n-2}} w_\nu^2\right)
\end{align*}
Using the pointwise estimate
\begin{align*}
 &\left|U_\nu^{\frac{2n}{n-2}}-v_\nu^{\frac{2n}{n-2}} +\frac{2n}{n-2} v_\nu^{\frac{n+2}{n-2}} w_\nu -\frac{n(n+2)}{(n-2)^2} U_\nu^{\frac{4}{n-2}} w_\nu^2\right|\\
&=\left|U_\nu^{\frac{2n}{n-2}}- (w_\nu+U_\nu)^{\frac{2n}{n-2}}+\frac{2n}{n-2}(w_\nu+U_\nu)^{\frac{n+2}{n-2}} w_\nu -\frac{n(n+2)}{(n-2)^2} U_\nu^{\frac{4}{n-2}} w_\nu^2\right|
 \\& \le C U_\nu^{\max\{0,\frac{4}{N-2}-1\}} |w_\nu|^{\min\{\frac{2N}{N-2},3\}} +C |w_\nu |^{\frac{2N}{N-2}},
 \end{align*}
 it follows that
 \begin{align*}
 &\int_\Omega}                \newcommand{\pa}{\partial \left|U_\nu^{\frac{2n}{n-2}}-v_\nu^{\frac{2n}{n-2}} +\frac{2n}{n-2} v_\nu^{\frac{n+2}{n-2}} w_\nu -\frac{n(n+2)}{(n-2)^2} U_\nu^{\frac{4}{n-2}} w_\nu^2\right|  \\&
 \le C \int_\Omega}                \newcommand{\pa}{\partial U_\nu^{\max\{0,\frac{4}{n-2}-1\}} |w_\nu|^{\min\{\frac{2n}{n-2},3\}} +C \int_\Omega}                \newcommand{\pa}{\partial |w_\nu |^{\frac{2n}{n-2}} \\&
 \le C \left(\int_\Omega}                \newcommand{\pa}{\partial |w_\nu |^{\frac{2n}{n-2}} \right)^{\frac{n-2}{n}\min\{\frac{n}{n-2}, \frac32\}}.
 \end{align*}
 By H\"older inequality and Cauchy inequality, we have
 \begin{align*}
\left| \int_\Omega}                \newcommand{\pa}{\partial (\mathcal{R}-1)v_\nu^{\frac{n+2}{n-2}} w_\nu \right| & \le C\left(\int_\Omega}                \newcommand{\pa}{\partial |\mathcal{R}-1|^{\frac{2n}{n+2}} v_\nu^{\frac{2n}{n-2}}\right)^{\frac{n+2}{2n}} \left(\int_\Omega}                \newcommand{\pa}{\partial |w_\nu|^{\frac{2n}{n-2}} \right)^{\frac{n-2}{2n}} \\&
\le \varepsilon}           \newcommand{\ud}{\mathrm{d} \left(\int_\Omega}                \newcommand{\pa}{\partial |w_\nu|^{\frac{2n}{n-2}} \right)^{\frac{n-2}{n}} + C(\varepsilon}           \newcommand{\ud}{\mathrm{d})\left(\int_\Omega}                \newcommand{\pa}{\partial |\mathcal{R}-1|^{\frac{2n}{n+2}} v_\nu^{\frac{2n}{n-2}}\right)^{\frac{n+2}{n}} .
 \end{align*}
 Finally, by Corollaries \ref{cor:coercive} and \ref{cor:coercive2}, we have
 \[
 \int_\Omega}                \newcommand{\pa}{\partial |\nabla w_\nu|^2 -b w_\nu^2- \frac{n+2}{n-2} U_\nu^{\frac{4}{n-2}} w_\nu^2\ge
 \frac{1}{C} \left(\int_\Omega}                \newcommand{\pa}{\partial |w_\nu|^{\frac{2n}{n-2}} \right)^{\frac{n-2}{n}}.
 \]
 Since $\displaystyle\int_\Omega}                \newcommand{\pa}{\partial |w_\nu|^{\frac{2n}{n-2}}  \to 0$, we have, by choosing $\varepsilon}           \newcommand{\ud}{\mathrm{d} $ being small, that
 \begin{align*}
 F(v_{\nu}) -F(U_\nu) &\le C \left(\int |w_\nu |^{\frac{2n}{n-2}} \right)^{\frac{n-2}{n}\min\{\frac{n}{n-2}, \frac32\}} +C\left(\int |\mathcal{R}-1|^{\frac{2n}{n+2}} v_\nu^{\frac{2n}{n-2}}\right)^{\frac{n+2}{n}}  \\& \quad - \frac{1}{2C} \left(\int |w_\nu|^{\frac{2n}{n-2}} \right)^{\frac{n-2}{n}}\\&
 \le C\left(\int |\mathcal{R}-1|^{\frac{2n}{n+2}} v_\nu^{\frac{2n}{n-2}}\right)^{\frac{n+2}{n}}.
 \end{align*}
 By Corollary \ref{cor:regular+singular}, we complete the proof.
 \end{proof}

Then we can show the estimate for all large time.
\begin{cor}\label{cor:fullestimate}
There exist real numbers $\gamma\in (0,1)$ and $t_0>0$ such that
\[
F(v(t)) -F_\infty \le  \left(\int_\Omega}                \newcommand{\pa}{\partial |\mathcal{R}-1|^{\frac{2n}{n+2}} v(x,t)^{\frac{2n}{n-2}}\,\ud x\right)^{\frac{n+2}{2n} (1+\gamma)}
\]
for all $t\ge t_0$.
\end{cor}
\begin{proof}
Suppose this is not true. Then, there exists a sequence of
times $\{t_\nu: \nu\in\mathbb{N}\}$ such that $t_\nu>\nu$ and
\[
F(v(t_\nu)) -F_\infty \ge  \left(\int_\Omega}                \newcommand{\pa}{\partial |\mathcal{R}(t_\nu)-1|^{\frac{2n}{n+2}} v(x,t_\nu)^{\frac{2n}{n-2}}\,\ud x\right)^{\frac{n+2}{2n} (1+\frac{1}{\nu})}
\]
for all $\nu\in\mathbb{N}$. By applying Proposition \ref{prob:subsequenceestimate} to this sequence $\{t_\nu: \nu\in\mathbb{N}\}$, there exists an infinite subset $I\subset\mathbb{N}$, a real number
$\gamma\in (0,1)$ and $C>0$ such that
\[
F(v(t_\nu)) -F_\infty \le  C \left(\int_\Omega}                \newcommand{\pa}{\partial |\mathcal{R}(t_\nu)-1|^{\frac{2n}{n+2}} v(x,t_\nu)^{\frac{2n}{n-2}}\,\ud x\right)^{\frac{n+2}{2n} (1+\gamma)}
\]
for all $\nu\in I$. Thus, we have
\[
1\le C \left(\int_\Omega}                \newcommand{\pa}{\partial |\mathcal{R}(t_\nu)-1|^{\frac{2n}{n+2}} v(x,t_\nu)^{\frac{2n}{n-2}}\,\ud x\right)^{\frac{n+2}{2n} (\gamma-\frac{1}{\nu})}
\]
for all $\nu\in I$. However, from Proposition \ref{prop:Mq}, we have
\[
\lim_{\nu\to\infty}\left(\int_\Omega}                \newcommand{\pa}{\partial |\mathcal{R}(t_\nu)-1|^{\frac{2n}{n+2}} v(x,t_\nu)^{\frac{2n}{n-2}}\,\ud x\right)^{\frac{n+2}{2n} (\gamma-\frac{1}{\nu})}=0.
\]
We have reached a contradiction.
\end{proof}

Now we can use a differential inequality of $F$ to obtain a decay estimate.
\begin{prop} \label{prop:converge-1}
There exist $\theta\in (0,1)$ and $C>0$ such that for all $T>1$, there holds
\[
\int_{T}^\infty M_2(t)^{1/2} \,\ud t \le C T^{-\theta},
\]
where $M_2$ is defined in \eqref{eq:definitionMq},
\end{prop}

\begin{proof}
It follows from  Corollary \ref{cor:fullestimate}, H\"older's inequality, and \eqref{eq:BN-5} that
\[
0\le F(v(t))-F_\infty \le \left(\int_\Omega}                \newcommand{\pa}{\partial |\mathcal{R}-1|^{\frac{2n}{n+2}} v(x,t)^{\frac{2n}{n-2}}\,\ud x\right)^{\frac{n+2}{2n} (1+\gamma)}  \le C M_2(t)^{\frac{1+\gamma}{2}}.
\]
It follows that
\begin{align*}
\frac{\ud }{\ud t}(F(v(t))-F_\infty )= -\frac{2(n-2)}{n+2} M_2(t) \le -C (F(v(t))-F_\infty)^{\frac{2}{1+\gamma}}.
\end{align*}
Hence,
\[
\frac{\ud }{\ud t}(F(v(t))-F_\infty  )^{\frac{\gamma-1}{\gamma+1}} \ge C \frac{1-\gamma}{1+\gamma}>0.
\]
It follows that
\[
F(v(t))-F_\infty  \le Ct^{-\frac{1+\gamma}{1-\gamma}}
\]
for sufficiently large $t$.
Then we have
\begin{align*}
\left(\int_{T}^{2T} M_2(s)^{1/2}\,\ud s\right)^2&  \le T \int_{T }^{2T } M_2(s) \,\ud s
\\& \le \frac{n+2}{n-2}  T  (F(v(T))- F(v(2T))) \\&
\le  \frac{n+2}{n-2}  T  (F(v(T))- F_\infty)  \\&
 \le C T^{-\frac{2\gamma}{1-\gamma}},
\end{align*}
where we used the monotonicity of $F$. It follows that
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:M2rate}
\int_{T}^\infty M_2(t)^{1/2} \,\ud t  = \sum_{k=0}^{\infty}  \int_{2^kT}^{2^{k+1}T} M_2(t)^{1/2} \,\ud t \le C T^{-\frac{\gamma}{1-\gamma}}\sum_{k=1}^{\infty} 2^{-\frac{\gamma}{1-\gamma} k}\le CT^{-\frac{\gamma}{1-\gamma}}.
\ee
This finishes the proof.
\end{proof}

We are ready to show the uniform boundedness, and uniform higher order estimates.

\begin{prop}\label{prop:boundedness}
For any $\varepsilon}           \newcommand{\ud}{\mathrm{d}>0$, there exists $T_0>0$ such that
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:closedness}
\|v(\cdot,t)-v(\cdot, T_0)\|_{L^{\frac{2n}{n-2}}(\Omega}                \newcommand{\pa}{\partial)} <\varepsilon}           \newcommand{\ud}{\mathrm{d}\quad\mbox{for all }t>T_0.
\ee
Consequently, there exists $C>0$ such that
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:closedness2}
v(x, t)\le C\quad\mbox{ in }\ \Omega \mbox{ for all }t>1.
\ee
\end{prop}

\begin{proof}
For $b>a>1$, using the pointwise estimate
\[
 |v(x,b)-v(x,a)|^{\frac{n}{n-2}}\le |v(x,b)^{\frac{n}{n-2}} -v(x,a) ^{\frac{n}{n-2}}|,
\]
we have
\begin{equation}  \label{eq:stable-5}
\begin{split}
\left(\int_{\Omega}                \newcommand{\pa}{\partial} |v(x,b)-v(x,a)|^{\frac{2n}{n-2}}\,\ud x\right)^{1/2} &\le\left( \int_{\Omega}                \newcommand{\pa}{\partial} |v(x,b)^{\frac{n}{n-2}} -v(x,a) ^{\frac{n}{n-2}} |^2\,\ud x\right)^{1/2} \\&
\le\left(\int_{\Omega}                \newcommand{\pa}{\partial} \Big(\int_{a}^{b} |\pa_t (v(x,t)^{\frac{n}{n-2}})| \,\ud t\Big) ^2 \,\ud x \right)^{1/2}\\&
\le \int_{a}^{b} \left(\int_{\Omega}                \newcommand{\pa}{\partial} | \pa_t (v(x,t)^{\frac{n}{n-2}})| ^2 \,\ud x\right)^{1/2}\,\ud t\\&
\le C \int_{a}^{\infty} M_2(t)^{1/2} \,\ud t\\
&\le C a^{-\theta},
\end{split}
\end{equation}
where we used Minkowski's integral inequality in the third inequality, and  Proposition \ref{prop:converge-1} in the last inequality. Hence, for any $\varepsilon}           \newcommand{\ud}{\mathrm{d}>0$, there exists $T_0>0$ such that \eqref{eq:closedness} holds.

Furthermore, it follows from \eqref{eq:Q} that
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:ellipticsplit}
-\Delta v-bv=\Big[(\mathcal{R}-1)v^{\frac{4}{n-2}}+v^{\frac{4}{n-2}}\Big]v.
\ee
Using \eqref{eq:closedness} and the fact that $M_{n/2}(t) \to 0$ as $t \to \infty$, we can apply the standard Brezis-Kato estimates \cite{BrezisKato} to the equation \eqref{eq:ellipticsplit} to obtain that there exists $p>\frac{2n}{n-2}$ and $C>0$ such that
\[
\|v(t)\|_{L^p(\Omega}                \newcommand{\pa}{\partial)} <C\quad\mbox{for all }t>T_0.
\]
Then we can further use Proposition \ref{prop:Mq} with large $q$ and apply the standard Moser iteration to the equation \eqref{eq:ellipticsplit} to obtain \eqref{eq:closedness2}.
\end{proof}

\begin{thm}\label{thm:smoothbound}
There exists  $C>0$ such that
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:smoothestimate}
\|v(\cdot,t)\|_{C^{2+\frac{n+2}{n-2}}(\overline\Omega)}\le C\quad\mbox{for all }t>1.
\ee
\end{thm}
\begin{proof}
By using \eqref{eq:closedness2}, it follows from Proposition 6.2 in \cite{DKV} (more precisely, its proof) that
\[
\frac{1}{C} d(x)\le v(\cdot,t)\le C d(x)\quad\mbox{for all }x\in\Omega,\ t>1,
\]
where $d(x):=\mbox{dist}(x,\partial\Omega)$. Then \eqref{eq:smoothestimate} follows from Theorem 5.1 in \cite{JX19}.
\end{proof}

Let us conclude this section with the proof of Theorem \ref{thm:critical}.

\begin{proof}[Proof of Theorem \ref{thm:critical}]
It follows from Proposition \ref{prop:boundedness},  Theorem \ref{thm:smoothbound} and \eqref{eq:BN-5}  that there exists a nonzero stationary solution $v_\infty$ of \eqref{eq:BN-3} such that
\[
\lim_{t\to\infty}\|v(\cdot,t)-v_\infty\|_{C^3(\overline\Omega)}=0.
\]
From \eqref{eq:stable-5}, we know that there exist $C>0$ and $\theta>0$ such that
\[
\|v(\cdot,t)-v_\infty\|_{L^\frac{2n}{n-2}(\Omega)}\le C t^{-\theta}\quad\mbox{for all }t>1.
\]
Using \eqref{eq:smoothestimate} and interpolation inequalities (Gagliardo--Nirenberg interpolation inequality and Morrey inequality), we have
\[
\|v(\cdot,t)-v_\infty\|_{C^1(\overline\Omega)}\le C t^{-\theta}\quad\mbox{for all }t>1.
\]
Since $v(\cdot,t)\equiv v_\infty\equiv 0$ on $\partial\Omega$, we have for all $x\in\Omega$ that
\[
\left|\frac{v(x,t)-v_\infty(x)}{v_\infty(x)}\right| \le \frac{\|v(\cdot,t)-v_\infty\|_{C^1(\overline\Omega)} \ d(x)}{d(x)/C} \le C \|v(\cdot,t)-v_\infty\|_{C^1(\overline\Omega)} \le C t^{-\theta}\quad\mbox{for all }t>1.
\]
That is,
\[
\left\|\frac{v(\cdot,t)}{v_\infty}-1\right\|_{L^\infty(\Omega}                \newcommand{\pa}{\partial)}\le Ct^{-\theta}\quad\mbox{for all }t\ge 1.
\]
Using \eqref{eq:smoothestimate} again, we have
\[
\left\|\frac{v(\cdot,t)}{v_\infty}\right\|_{C^{1+\frac{n+2}{n-2}}(\overline\Omega)}\le C\quad\mbox{for all }t>1.
\]
Then by interpolation inequalities, we have
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:relativealgebraic}
\left\|\frac{v(\cdot,t)}{v_\infty}-1\right\|_{C^2(\overline\Omega}                \newcommand{\pa}{\partial)}\le Ct^{-\theta}\quad\mbox{for all }t\ge 1.
\ee

Now let us assume that $\Omega$ satisfies the condition \eqref{condition}. Then there exists $C>0$ such that for all $\varphi\in H^1_0$ satisfying
\[
-\Delta \varphi- b \varphi -\frac{n+2}{n-2}v_\infty^{\frac{4}{n-2}}\varphi=f\quad\mbox{in }\Omega,
\]
there holds
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:invertibleestimate}
\|\varphi\|_{H^1_0(\Omega)}\le C \|f\|_{L^2(\Omega)}.
\ee

On one hand, using the equation of $v_\infty$, we have
\begin{align*}
&F(v(t))-F_\infty
\\&= \int_{\Omega}                \newcommand{\pa}{\partial } \Big( |\nabla v(\cdot,t)|^2 -b v(\cdot,t)^2 -\frac{n-2}{n} v(\cdot,t)^{\frac{2n}{n-2}}\Big) - \int_{\Omega}                \newcommand{\pa}{\partial } \Big( |\nabla v_\infty|^2 -b v_\infty^2 -\frac{n-2}{n} v_\infty^{\frac{2n}{n-2}}\Big)\\
&\quad-2\int_\Omega}                \newcommand{\pa}{\partial \Big(\nabla v_\infty\nabla (v(\cdot,t)-v_\infty) -b v_\infty (v(\cdot,t)-v_\infty)- v_\infty^{\frac{n+2}{n-2}} (v(\cdot,t)-v_\infty)\Big)\\
&= \int_{\Omega}                \newcommand{\pa}{\partial }  |\nabla (v(\cdot,t)- v_\infty)|^2 -b (v(\cdot,t)- v_\infty)^2 -\frac{n-2}{n} \Big(v(\cdot,t)^{\frac{2n}{n-2}}-v_\infty^{\frac{2n}{n-2}}-\frac{2n}{n-2} (v(\cdot,t)- v_\infty)\Big)\\
&\le C \|v(\cdot,t)-v_\infty\|_{H^1_0(\Omega)}^2.
\end{align*}
On the other hand,  using the equation of $v_\infty$, we have
\begin{align*}
&\|\Delta v(\cdot, t)+b v(\cdot,t)+v(\cdot,t)^{\frac{n+2}{n-2}} \|_{L^2(\Omega)}\\
&=\|\Delta (v(\cdot, t)-v_\infty)+b (v(\cdot,t)-v_\infty)+v(\cdot,t)^{\frac{n+2}{n-2}}-v_\infty^{\frac{n+2}{n-2}} \|_{L^2(\Omega)}\\
&\ge \left\|\Delta (v(\cdot, t)-v_\infty)+b (v(\cdot,t)-v_\infty)+\frac{n+2}{n-2} v_\infty^\frac{4}{n-2}(v(\cdot,t)-v_\infty)\right\|_{L^2(\Omega)}\\
&\quad- \left\|v(\cdot,t)^{\frac{n+2}{n-2}}-v_\infty^{\frac{n+2}{n-2}}-\frac{n+2}{n-2} v_\infty^\frac{4}{n-2}(v(\cdot,t)-v_\infty)\right\|_{L^2(\Omega)}\\
&\ge \frac 1C \|v(\cdot, t)-v_\infty\|_{H^1_0(\Omega)}- C \left\| v_\infty ^{\max(0, \frac{4}{n-2}-1)}|v(\cdot,t)-v_\infty|^{\min(2, \frac{n+2}{n-2})}\right\|_{L^2(\Omega)}\\
&\ge \frac 1C \|v(\cdot, t)-v_\infty\|_{H^1_0(\Omega)} -o(1) \|v(\cdot, t)-v_\infty\|_{H^1_0(\Omega)} \quad \mbox{for large }t\\
&\ge \frac 1C \|v(\cdot, t)-v_\infty\|_{H^1_0(\Omega)}\quad \mbox{for large }t,
\end{align*}
where we used \eqref{eq:invertibleestimate} in the second inequality. Combining these two inequalities together, we have
\[
F(v(t))-F_\infty\le \|\Delta v(\cdot, t)+b v(\cdot,t)+v(\cdot,t)^{\frac{n+2}{n-2}} \|_{L^2(\Omega)}^2=CM_2(t).
\]
It follows that
\begin{align*}
\frac{\ud }{\ud t}(F(v(t))-F_\infty )= -\frac{2(n-2)}{n+2} M_2(t) \le -C (F(v(t))-F_\infty),
\end{align*}
and thus,
\[
F(v(t))-F_\infty  \le Ce^{-\gamma t}
\]
for some $C>0,\gamma>0$.
Hence, the proof of \eqref{eq:M2rate} will give
\[
\int_{T}^\infty M_2(t)^{1/2} \,\ud t \le Ce^{-\gamma t}.
\]
From \eqref{eq:stable-5}, we obtain that
\[
\|v(\cdot,t)-v_\infty\|_{L^\frac{2n}{n-2}(\Omega)} \le Ce^{-\gamma t}.
\]
Then the proof of \eqref{eq:relativealgebraic} gives
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:relativeexp}
\left\|\frac{v(\cdot,t)}{v_\infty}-1\right\|_{C^2(\overline\Omega}                \newcommand{\pa}{\partial)}\le Ce^{-\gamma t}\quad\mbox{for all }t\ge 1.
\ee
This finishes the proof of Theorem \ref{thm:critical}.
\end{proof}

\section{Subcritical case}\label{sec:subcritical}
In this last section, we consider the Sobolev subcritical case \eqref{eq:subcritical2}, and prove Theorem \ref{thm:subcritical}. The proof is similar to that of Theorem \ref{thm:critical}.
\begin{proof}[Proof of Theorem \ref{thm:subcritical}]
First, we know from Proposition 6.2 in \cite{DKV} that
\[
\frac{1}{C} d(x)\le v(\cdot,t)\le C d(x)\quad\mbox{for all }x\in\Omega,\ t>1,
\]
Secondly, it follows from Theorem 1.1 in \cite{FS} and Theorem 5.1 in \cite{JX19}  that there exists a nonzero stationary solution $v_\infty$ of \eqref{eq:subcritical2} such that
\[
\lim_{t\to\infty}\|v(\cdot,t)-v_\infty\|_{C^3(\overline\Omega)}=0.
\]
Let
\begin{align*}
F(v(t))&= \int_{\Omega}                \newcommand{\pa}{\partial}\left( |\nabla v(x,t)|^2 -\frac{2}{p+1} v(x,t)^{p+1}\right)\,\ud x, \\
\mathcal{R}&= -v^{-p} \Delta v=1-\frac{p\partial_t v}{v},\\
M_q(t)&=\int_{\Omega}                \newcommand{\pa}{\partial} |\mathcal{R}-1|^q v^{p+1} \,\ud x,
\end{align*}
where $q\ge 1$. Note that
\[
\frac{\ud }{\ud t} F(v(t))= -2 \int_\Omega}                \newcommand{\pa}{\partial (\Delta v+v^p) \pa_t v= -2p\int_\Omega}                \newcommand{\pa}{\partial |\pa_t v|^2 v^{p-1}=-\frac{2}{p }M_2(t).
\]
Hence $F(v(\cdot,t))$ deceases to $F(v_\infty)$ as $t\to\infty$. Furthermore,
\[
\ud  F(v)=-2\Delta v-2v^p=2(\mathcal{R}-1)v^p.
\]
Hence,
\[
\|\ud  F(v(\cdot,t))\|_{L^2(\Omega)} \le C M_{2}(t)^{1/2}.
\]
Then it follows from Proposition 6.1 in \cite{FS} that there exist $C>0, T_0>0, \gamma>0$ such that for all $t>T_0$, we have
\[
F(v(t))-F(v_\infty)\le C\|\ud  F(v(\cdot,t))\|_{L^2(\Omega)}^{1+\gamma}\le C M_{2}(t)^{\frac{1+\gamma}{2}}.
\]
Therefore, similar to the proof of Proposition \ref{prop:converge-1}, there exist $\theta\in (0,1)$ and $C>0$ such that for all $T>1$, there holds
\[
\int_{T}^\infty M_2(t)^{1/2} \,\ud t \le C T^{-\theta},
\]
Then it follows from the proof of \eqref{eq:stable-5} that
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:subM2exp}
\|v(\cdot,t)-v_\infty\|_{L^{p+1}(\Omega}                \newcommand{\pa}{\partial)}\le  C \int_t^\infty M_2(s)^{1/2}\,\ud s\le C t^{-\theta},
\ee
Using Theorem 5.1 in \cite{JX19} (which is the regularity estimate) and interpolation inequalities, we have
\[
\|v(\cdot,t)-v_\infty\|_{C^1(\overline\Omega)}\le C t^{-\theta}\quad\mbox{for all }t>1.
\]
Since $v(\cdot,t)\equiv v_\infty\equiv 0$ on $\partial\Omega$, we have for all $x\in\Omega$ that
\[
\left|\frac{v(x,t)-v_\infty(x)}{v_\infty(x)}\right| \le C \|v(\cdot,t)-v_\infty\|_{C^1(\overline\Omega)} \le C t^{-\theta}\quad\mbox{for all }t>1.
\]
That is,
\[
\left\|\frac{v(\cdot,t)}{v_\infty}-1\right\|_{L^\infty(\Omega}                \newcommand{\pa}{\partial)}\le Ct^{-\theta}\quad\mbox{for all }t\ge 1.
\]
Using Theorem 5.1 in \cite{JX19} again, we have
\[
\left\|\frac{v(\cdot,t)}{v_\infty}\right\|_{C^{1+p}(\overline\Omega)}\le C\quad\mbox{for all }t>1.
\]
Then by interpolation inequalities, we have
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:relativealgebraicsub}
\left\|\frac{v(\cdot,t)}{v_\infty}-1\right\|_{C^2(\overline\Omega}                \newcommand{\pa}{\partial)}\le Ct^{-\theta}\quad\mbox{for all }t\ge 1.
\ee

Now let us assume that $\Omega$ satisfies the condition \eqref{condition2}. Then there exists $C>0$ such that for all $\varphi\in H^1_0$ satisfying
\[
-\Delta \varphi-pv_\infty^{p-1}\varphi=f\quad\mbox{in }\Omega,
\]
there holds
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:invertibleestimate2}
\|\varphi\|_{H^1_0(\Omega)}\le C \|f\|_{L^2(\Omega)}.
\ee

As before, on one hand, using the equation of $v_\infty$, we have
\begin{align*}
F(v(t))-F_\infty \le C \|v(\cdot,t)-v_\infty\|_{H^1_0(\Omega)}^2.
\end{align*}
On the other hand,  using the equation of $v_\infty$, we have
\begin{align*}
&\|\Delta v(\cdot, t)+v(\cdot,t)^{p} \|_{L^2(\Omega)}\\
&\ge \left\|\Delta (v(\cdot, t)-v_\infty)+p v_\infty^{p-1}(v(\cdot,t)-v_\infty)\right\|_{L^2(\Omega)}\\
&\quad- \left\|v(\cdot,t)^p-v_\infty^p-p v_\infty^{p-1}(v(\cdot,t)-v_\infty)\right\|_{L^2(\Omega)}\\
&\ge \frac 1C \|v(\cdot, t)-v_\infty\|_{H^1_0(\Omega)}- C \left\| v_\infty ^{\max(0, p-2)}|v(\cdot,t)-v_\infty|^{\min(2, p)}\right\|_{L^2(\Omega)}\\
&\ge \frac 1C \|v(\cdot, t)-v_\infty\|_{H^1_0(\Omega)} -o(1) \|v(\cdot, t)-v_\infty\|_{H^1_0(\Omega)} \quad \mbox{for large }t\\
&\ge \frac 1C \|v(\cdot, t)-v_\infty\|_{H^1_0(\Omega)}\quad \mbox{for large }t,
\end{align*}
where we used \eqref{eq:invertibleestimate2} in the second inequality. Combining these two inequalities together, we have
\[
F(v(t))-F_\infty\le C\|\Delta v(\cdot, t)+v(\cdot,t)^{p} \|_{L^2(\Omega)}^2\le CM_2(t).
\]
It follows that
\begin{align*}
\frac{\ud }{\ud t}(F(v(t))-F_\infty )= -\frac{1}{p} M_2(t) \le -C (F(v(t))-F_\infty),
\end{align*}
and thus,
\[
F(v(t))-F_\infty  \le Ce^{-\gamma t}
\]
for some $C>0,\gamma>0$.
Hence, the proof of \eqref{eq:M2rate} will give
\[
\int_{T}^\infty M_2(t)^{1/2} \,\ud t \le Ce^{-\gamma t}.
\]
From \eqref{eq:subM2exp}, we obtain that
\[
\|v(\cdot,t)-v_\infty\|_{L^{p+1}(\Omega)} \le Ce^{-\gamma t}.
\]
Then the proof of \eqref{eq:relativealgebraicsub} gives
\begin{equation}}      \newcommand{\ee}{\end{equation}\label{eq:relativeexpsub}
\left\|\frac{v(\cdot,t)}{v_\infty}-1\right\|_{C^2(\overline\Omega}                \newcommand{\pa}{\partial)}\le Ce^{-\gamma t}\quad\mbox{for all }t\ge 1.
\ee

This finishes the proof of Theorem \ref{thm:subcritical}.
\end{proof}