Hostname: page-component-586b7cd67f-tf8b9 Total loading time: 0 Render date: 2024-11-26T13:08:26.464Z Has data issue: false hasContentIssue false

Local regularity for nonlinear elliptic and parabolic equations with anisotropic weights

Published online by Cambridge University Press:  28 April 2023

Changxing Miao
Affiliation:
Beijing Computational Science Research Center, Beijing 100193, China Institute of Applied Physics and Computational Mathematics, P.O. Box 8009, Beijing 100088, China ([email protected])
Zhiwen Zhao
Affiliation:
Beijing Computational Science Research Center, Beijing 100193, China ([email protected])
Rights & Permissions [Opens in a new window]

Abstract

The main purpose of this paper is to capture the asymptotic behaviour for solutions to a class of nonlinear elliptic and parabolic equations with the anisotropic weights consisting of two power-type weights of different dimensions near the degenerate or singular point, especially covering the weighted p-Laplace equations and weighted fast diffusion equations. As a consequence, we also establish the local Hölder estimates for their solutions in the presence of single power-type weights.

Type
Research Article
Copyright
© The Author(s), 2023. Published by Cambridge University Press on Behalf of The Edinburgh Mathematical Society.

1. Introduction and main results

In this paper, we will use the De Giorgi truncation method [Reference De Giorgi16] to study the local behaviour of solutions for a class of nonlinear elliptic and parabolic equations with the weights comprising two power-type weights of different dimensions. For that purpose, we should first establish the corresponding anisotropic weighted Sobolev embedding theorems and Poincaré’s inequality, which are fundamental tools to investigate relevant Sobolev spaces and partial differential equations. The former has recently been established by Li and Yan [Reference Li and Yan30], whose results improve and extend the classical Caffarelli–Kohn–Nirenberg type inequalities in [Reference Caffarelli, Kohn and Nirenberg10]. With regard to the latter, we prepare to prove that this type of anisotropic weights belongs to the Muckenhoupt class A q under certain conditions, and then the anisotropic weighted Poincaré inequality is obtained by utilizing the theories of A q-weights, $1 \lt q \lt \infty$, see Section 2 below for the finer details. This is another major novelty of this paper besides the regularity results with anisotropic weights. For more relevant investigations on weighted Sobolev and Poincaré inequalities, see [Reference Badiale and Tarantello2, Reference Bahouri, Chemin and Gallagher3, Reference Cabré and Ros-Oton8, Reference Cabré, Ros-Oton and Serra9, Reference Chanillo and Wheeden11, Reference Chua and Wheeden14, Reference Lin31, Reference Nguyen and Squassina33, Reference Nguyen and Squassina34] and the references therein.

The anisotropy of the weights considered in this paper comes from two power-type weights of different dimensions. This complex weighted form will bring great difficulties of analyses, computations and discussions in the following proofs, especially the findings for regular indices which make this type of anisotropic weights become A q-weights. The mathematical formulations and main results for the considered nonlinear elliptic and parabolic problems with anisotropic weights are, respectively, presented as follows.

1.1. The nonlinear elliptic equations with anisotropic weights

Consider a bounded smooth domain $\Omega\subset\mathbb{R}^{n}$ with $0\in\Omega$ and $n\geq2$. With regard to the weighted elliptic equations, we mainly study the local regularity of solution to the following problem

(1.1)\begin{align} \begin{cases} \mathrm{div}(Aw|\nabla u|^{p-2}\nabla u)=0,& \mathrm{in}\ \Omega,\\ 0\leq u\leq\overline{M} \lt \infty,&\mathrm{in}\ \Omega, \end{cases} \end{align}

where $w=|x^{\prime}|^{\theta_{1}}|x|^{\theta_{2}}$, the values of θ 1 and θ 2 are assumed in the following theorems, $1 \lt p \lt n+\theta_{1}+\theta_{2}$, $\overline{M}$ is a given positive constant, $A(x)=(a_{ij}(x))_{n\times n}$ is symmetric and satisfies

(1.2)\begin{align} \lambda^{-1}|\xi|^{2}\leq\sum^{n}_{i,j=1}a_{ij}(x)\xi_{i}\xi_{j}\leq\lambda|\xi|^{2},\quad\lambda\geq1\;\; \text{for a.e. } x\in\Omega \text{ and all } \xi\in\mathbb{R}^{n}. \end{align}

Here and throughout this paper, we use superscript prime to denote $(n-1)$-dimensional variables and domains, such as xʹ and Bʹ. Moreover, in the following, we simplify the notations $B_{R}(0)$ and $B_{R}^{\prime}(0^{\prime})$ as B R and $B^{\prime}_{R}$, respectively, where R > 0. The prototype equation is the anisotropic weighted p-Laplacian, that is, the equation in the case when A = I in Equation (1.1). Remark that the origin can be called the degenerate or singular point of the weight. For example, if $\theta_{1} \gt 0$, $\theta_{2} \gt 0$, then the weight $w=|x^{\prime}|^{\theta_{1}}|x|^{\theta_{2}}\rightarrow0$ as $|x|\rightarrow0$, while for $\theta_{1} \lt 0$, $\theta_{2} \lt 0$, it blows up as $|x|$ tends to zero. For the former, the origin is called the degenerate point of the weight, while it is called the singular point for the latter.

For the weighted elliptic problem (1.1), Fabes et al. [Reference Fabes, Kenig and Serapioni22] established the local Hölder regularity of weak solutions under the case of $\theta_{1}=0,\,\theta_{2} \gt -n$ and p = 2. However, the value of Hölder index α obtained in [Reference Fabes, Kenig and Serapioni22] is not explicit. Recently, Dong et al. [Reference Dong, Li and Yang21] utilized spherical harmonic expansion to find the exact value of index α for the solution near the degenerate point of the weight. To be precise, for problem (1.1) with Ω replaced by B R, R > 0, let $n\geq2$, $\theta_{1}=0$, $\theta_{2}=p=2$ and $A=\kappa(x)I$, where κ satisfies that $\lambda^{-1}\leq\kappa\leq\lambda$ in B R and $\int_{\mathbb{S}^{n-1}}\kappa x_{i}=0$, $i=1,2,\ldots,n$. Based on these assumed conditions, they derived

\begin{align*} u(x)=u(0)+O(1)|x|^{\alpha},\quad\alpha=\frac{-n+\sqrt{n^{2}+4\tilde{\lambda}_{1}}}{2}\quad\mathrm{in}\;B_{R/2}, \end{align*}

where O(1) represents some quantity such that $|O(1)|\leq C=C(n,\lambda,\overline{M}),$ $\tilde{\lambda}_{1}\leq n-1$ is the first non-zero eigenvalue of the following eigenvalue problem:

\begin{align*} -\mathrm{div}_{\mathbb{S}^{n-1}}(\kappa(\xi)\nabla_{\mathbb{S}^{n-1}}u(\xi))=\tilde{\lambda} \kappa(\xi) u(\xi),\qquad\xi\in\mathbb{S}^{n-1}. \end{align*}

In particular, $\tilde{\lambda}_{1}=n-1$ if A = I. See Lemmas 2.2 and 5.1 in [Reference Dong, Li and Yang21] for more details. By finding the explicit exponent α, they succeeded in solving the optimal gradient blow-up rate for solution to the insulated conductivity problem in dimensions greater than two, which has been previously regarded as a challenging problem. By their investigations in [Reference Dong, Li and Yang20, Reference Dong, Li and Yang21], we realize that the Hölder regularity for solutions to the weighted elliptic problem (1.1) is in close touch with the insulated conductivity problem arising from composite materials. Then the study on the regularity for weighted elliptic problem (1.1) is a topic of theoretical interest and also of great relevance to applications for the insulated composites. It is worth emphasizing that when p > 2, the exact value of index α still remains open. In addition, with regard to the Hölder regularity for nonlinear degenerate elliptic equations without weights, we refer to [Reference Lindqvist32, Reference Serrin36, Reference Trudinger38] and the references therein.

Before stating the definition of weak solution to problem (1.1), we first introduce some notations. Throughout this paper, we will use $L^{p}(\Omega,w)$ and $W^{1,p}(\Omega,w)$ to represent weighted L p space and weighted Sobolev space with their norms, respectively, written as

\begin{align*} \begin{cases} \|u\|_{L^{p}(\Omega,w)}=\left(\int_{\Omega}|u|^{p}w\,{\rm d}x\right)^{\frac{1}{p}},\\ \|u\|_{W^{1,p}(\Omega,w)}=\left(\int_{\Omega}|u|^{p}w\,{\rm d}x\right)^{\frac{1}{p}}+\left(\int_{\Omega}|\nabla u|^{p}wdx\right)^{\frac{1}{p}}. \end{cases} \end{align*}

We say that $u\in W^{1,p}(\Omega,w)$ is a weak solution of problem (1.1) if

\begin{align*} \int_{\Omega}Aw|\nabla u|^{p-2}\nabla u\cdot\nabla\varphi\, {\rm d}x=0,\quad\forall\varphi\in W_{0}^{1,p}(\Omega,w). \end{align*}

For later use, we introduce the following indexing sets:

\begin{align*} \begin{cases} \mathcal{A}=\{(a,b): a \gt -(n-1),\,b\geq0\},\\ \mathcal{B}=\{(a,b):a \gt -(n-1),\,b \lt 0,\,a+b \gt -n\},\\ \mathcal{C}_{q}=\{(a,b):a \lt (n-1)(q-1),\,b\leq0\},\quad q \gt 1,\\ \mathcal{D}_{q}=\{(a,b):a \lt (n-1)(q-1),\,b \gt 0,\,a+b \lt n(q-1)\},\quad q \gt 1,\\ \mathcal{F}=\{(a,b):a+b \gt -(n-1)\}. \end{cases} \end{align*}

The local behaviour of solution to problem (1.1) near the degenerate or singular point of the anisotropic weight is captured as follows.

Theorem 1.1. For $n\geq2$, $(\theta_{1},\theta_{2})\in(\mathcal{A}\cup\mathcal{B})\cap(\mathcal{C}_{q}\cup\mathcal{D}_{q})\cap\mathcal{F}$, $1 \lt q \lt p \lt n+\theta_{1}+\theta_{2}$, let u be a weak solution of problem (1.1) with $\Omega=B_{1}$. Then there exists a constant $0 \lt \alpha \lt 1$ depending only on $n,p,q,\theta_{1},\theta_{2},\lambda,$ such that

(1.3)\begin{align} u(x)=u(0)+O(1)|x|^{\alpha}\quad\text{for all}\;x\in B_{1/2}, \end{align}

where O(1) denotes some quantity satisfying that $|O(1)|\leq C=C(n,p,q,\theta_{1},\theta_{2},\lambda,\overline{M}).$

Remark 1.2. If the considered domain B 1 is replaced with $B_{R_{0}}$ for any given $R_{0} \gt 0$ in Theorems 1.1 and 1.4, then by applying their proofs with minor modification, we obtain that Equations (1.3) and (1.4) also hold with $B_{1/2}$ replaced by $B_{R_{0}/2}$. In this case, the constant C will depend on R 0, but the index α depends not on it.

Remark 1.3. The result in Theorem 1.1 can be extended to general degenerate elliptic equations as follows:

\begin{align*} \begin{cases} \mathrm{div}(\mathcal{G}(x,\nabla u))=0,& \mathrm{in}\;B_{R_{0}},\\ 0\leq u\leq\overline{M} \lt \infty,&\mathrm{in}\;B_{R_{0}}, \end{cases} \end{align*}

where $R_{0} \gt 0$, $\mathcal{G}:B_{R_{0}}\times\mathbb{R}^{n}\rightarrow\mathbb{R}^{n}$ is a Carathéodory function such that for a.e. $x\in B_{R_{0}}$ and any $\xi\in\mathbb{R}^{n}$, there holds

\begin{align*} \lambda^{-1}w(x)|\xi|^{p}\leq\mathcal{G}(x,\xi)\cdot\xi,\quad|\mathcal{G}(x,\xi)|\leq\lambda w(x)|\xi|^{p-1},\quad w(x)=|x'|^{\theta_{1}}|x|^{\theta_{2}}. \end{align*}

Here $\lambda\geq1$, $(\theta_{1},\theta_{2})\in(\mathcal{A}\cup\mathcal{B})\cap(\mathcal{C}_{q}\cup\mathcal{D}_{q})\cap\mathcal{F}$, $1 \lt q \lt p \lt n+\theta_{1}+\theta_{2}$. In fact, it only needs to slightly modify the proof of Lemma 3.1 below for the purpose of achieving this generalization.

When $\theta_{1}=0$, the above weight becomes a single power-type weight. In this case, we establish the Hölder estimates as follows.

Theorem 1.4. For $n\geq2$, $\theta_{1}=0,\,\theta_{2} \gt -(n-1)$, $1 \lt p \lt n+\theta_{2}$, let u be a bounded weak solution of problem (1.1) with $\Omega=B_{1}$. Then there exist a small constant $0 \lt \alpha=\alpha(n,p,\theta_{2},\lambda) \lt 1$ and a large constant $0 \lt C=C(n,p,\theta_{2},\lambda,\overline{M})$ such that

(1.4)\begin{align} |u(x)-u(y)|\leq C|x-y|^{\alpha}\quad\text{for all}\;x,y\in B_{1/2}. \end{align}

Observe that when $\theta_{1}=0$, Equation (1.1) will become degenerate elliptic equation in any domain away from the origin, then we can directly establish its Hölder regularity in these regions by using the interior Hölder estimates for degenerate elliptic equation. This, in combination with Remark 1.2 and Theorem 1.4, gives the following corollary.

Corollary 1.5. For $n\geq2$, $\theta_{1}=0,\,\theta_{2} \gt -(n-1)$, $1 \lt p \lt n+\theta_{2}$, let u be a weak solution of problem (1.1). Then u is locally Hölder continuous in Ω, that is, for any compact subset $K\subset\Omega$, there exists two constants $0 \lt \alpha=\alpha(n,p,\theta_{2},\lambda) \lt 1$ and $C=C(\mathrm{dist}(K,\partial\Omega),\mathrm{dist}(0,\partial\Omega),n,p,\theta_{2},\lambda,\overline{M}) \gt 0$ such that Equation (1.4) holds with $B_{1/2}$ replaced by K.

1.2. The nonlinear parabolic equations with anisotropic weights

Let $0\in\Omega\subset\mathbb{R}^{n}$, $n\geq2$, be defined as above. The second problem of interest is concerned with studying the local regularity of solution to the weighted nonlinear parabolic equation as follows:

(1.5)\begin{align} \begin{cases} w_{1}\partial_{t}u^{p}-\mathrm{div}(Aw_{2}\nabla u)=0,& \mathrm{in}\;\Omega_{T},\\ 0 \lt \overline{m}\leq u\leq\overline{M} \lt \infty,&\mathrm{in}\;\Omega_{T}, \end{cases} \end{align}

where $\Omega_{T}=\Omega\times(-T,0]$, T > 0, $w_{1}=|x^{\prime}|^{\theta_{1}}|x|^{\theta_{2}}$, $w_{2}=|x^{\prime}|^{\theta_{3}}|x|^{\theta_{4}}$, $p\geq1$, the ranges of θ i, $i=1,2,3,4$ are prescribed in the following theorems, $\overline{m}$ and $\overline{M}$ are two given positive constants, the symmetric matrix $A=(a_{ij}(x))_{n\times n}$ satisfies the uniformly elliptic condition in Equation (1.2). When $\theta_{i}=0$, $i=1,2,3,4$, and A = I, Equation (1.5) becomes fast diffusion equation, whose relevant mathematical problem is modeled by

(1.6)\begin{align} \begin{cases} \partial_{t}u^{p}-\Delta u=0,&\mathrm{in}\;\Omega\times(0,\infty),\\ u=0,&\mathrm{on}\;\partial\Omega\times(0,\infty),\\ u(x,0)=u_{0}(x)\geq0.& \end{cases} \end{align}

In physics, Equation (1.6) can be used to describe fast diffusion phenomena occurring in gas kinetics, plasmas and thin liquid film dynamics. For more related applications and physical explanations, see [Reference Daskalopoulos and Kenig15, Reference Vázquez39] and the references therein.

For problem (1.6), it is well known that when $u_{0}(x)\not\equiv0$, there exists a finite extinction time $T^{\ast} \gt 0$ such that $u(\cdot,t)\equiv0$ in Ω if $t\in[T^{\ast},\infty)$ and $u(\cdot,t) \gt 0$ in Ω if $t\in(0,T^{\ast})$. This, together with the continuity of u (see Chen-DiBenedetto [Reference Chen and DiBenedetto12]), indicates that for any $U\subset\subset\Omega\times(0,T^{\ast})$, there exist two positive constants $\overline{m}$ and $\overline{M}$ such that $0 \lt \overline{m}\leq u\leq\overline{M} \lt \infty$ for $(x,t)\in U$. This fact motivates our investigation on the local regularity of weak solution for the corresponding weighted problem (1.5). In particular, it can be called the weighted fast diffusion equation when A = I in Equation (1.5). For the fast diffusion problem (1.6), the regularity of solution and its asymptotic behaviour near extinction time have been extensively studied, for example, see [Reference Bonforte and Vázquez7, Reference Chen and DiBenedetto12, Reference Daskalopoulos and Kenig15, Reference DiBenedetto, Gianazza and Vespri17Reference DiBenedetto, Kwong and Vespri19, Reference Jin and Xiong27Reference Kwong29, Reference Sacks35] for the regularity and [Reference Akagi1, Reference Berryman and Holland4Reference Bonforte, Grillo and Vázquez6, Reference Feireisl and Simondon23] for the extinction behaviour, respectively. In particular, Jin and Xiong recently established a priori Hölder estimates for the solution to a weighted nonlinear parabolic equation in Theorem 3.1 of [Reference Jin and Xiong28], which is critical to the establishment of optimal global regularity for fast diffusion equation with any $1 \lt p \lt \infty$. Their results especially answer the regularity problem proposed by Berryman and Holland [Reference Berryman and Holland4]. It is worth pointing out that the degeneracy of weight in [Reference Jin and Xiong28] is located at the boundary. By contrast, the degeneracy or singularity of the weights considered in this paper lies in the interior. This will lead to some distinct differences in terms of the establishments of Hölder estimates under these two cases. Moreover, since the weights considered in this paper take more sophisticated forms comprising two power-type weights of different dimensions, it greatly increases the difficulties of analyses and calculations. With regard to the regularity for weighted parabolic problem in the case when p = 1 in Equation (1.5), we refer to [Reference Chiarenza and Serapioni13, Reference Gutiérrez and Wheeden25, Reference Surnachev37] and the references therein.

The weighted L p space and weighted Sobolev spaces with respect to space variable have been defined above. Similarly, for a weight w, let $W^{1,p}(\Omega_{T},w)$ represent the corresponding weighted Sobolev spaces in (x, t) with its norm given by

\begin{align*} \|u\|_{W^{1,p}(\Omega_{T},w)}=\left(\int_{\Omega_{T}}|u|^{p}w\,{\rm d}x\,{\rm d}t\right)^{\frac{1}{p}}+\left(\int_{\Omega_{T}}\left(\left|\partial_{t}u\right|^{p}+\left|\nabla u\right|^{p}\right)w\,{\rm d}x\,{\rm d}t\right)^{\frac{1}{p}}. \end{align*}

We say that $u\in W^{1,2}(\Omega_{T},w_{2})$ is a weak solution of problem (1.5) if

\begin{align*} \int^{t_{2}}_{t_{1}}\int_{\Omega}\left(w_{1}\partial_{t}u^{p}\varphi+Aw_{2}\nabla u\nabla\varphi\right)\,{\rm d}x\,{\rm d}t=0 \end{align*}

for any $-T\leq t_{1} \lt t_{2}\leq0$ and $\varphi\in C^{1}(\Omega_{T})$, which vanishes on $\partial\Omega\times(-T,0)$.

Introduce the following index conditions:

  • (S1) let $n\geq4$ and $1+2/(n-1) \lt q \lt 2$, if $(\theta_{1},\theta_{2})\in(\mathcal{A}\cup\mathcal{B})\cap\mathcal{C}_{q}$;

  • (S2) let $n\geq3$ and $1+2/n \lt q \lt 2$, if $(\theta_{1},\theta_{2})\in\mathcal{A}\cap\mathcal{D}_{q}$.

For the local behaviour of solution to problem (1.5), we have

Theorem 1.6. Suppose that $p\geq1$, $(\theta_{1},\theta_{2})$ satisfies condition $\mathrm{(}\mathbf{S}1\mathrm{)}$ or $\mathrm{(}\mathbf{S}2\mathrm{)}$, $(\theta_{3},\theta_{4})\in\mathcal{A}\cup\mathcal{B},$ $\theta_{1}+\theta_{2}\geq\theta_{3}+\theta_{4}=2$, $\theta_{1}/\theta_{3}=\theta_{2}/\theta_{4}$, $\theta_{3},\theta_{4}\neq0$. Let u be a weak solution of problem (1.5) with $\Omega\times(-T,0]=B_{1}\times(-1,0]$. Then there exists a small constant $0 \lt \alpha=\alpha(n,p,q,\theta_{1},\theta_{2},\theta_{3},\lambda,\overline{m},\overline{M}) \lt 1$ such that for any $t_{0}\in (-1/4,0)$,

(1.7)\begin{align} u(x,t)=u(0,t_{0})+O(1)\big(|x|+\sqrt[\theta_{1}+\theta_{2}]{|t-t_{0}|}\big)^{\alpha},\quad \forall (x,t)\in B_{1/2}\times(-1/4,t_{0}], \end{align}

where O(1) satisfies that $|O(1)|\leq C=C(n,p,q,\theta_{1},\theta_{2},\theta_{3},\lambda,\overline{m},\overline{M}).$

Remark 1.7. We provide here explanations for the index conditions $\mathrm{(}\mathbf{S}1\mathrm{)}$ and $\mathrm{(}\mathbf{S}2\mathrm{)}$. Observe that if $(\theta_{1},\theta_{2})\in C_{q}$, $1 \lt q \lt 2$ and $\theta_{1}+\theta_{2}\geq2$, then we have $(n-1)(q-1) \gt 2$, which requires that $n\geq4$ and $q \gt 2/(n-1)+1$. Similarly, if $(\theta_{1},\theta_{2})\in D_{q}$, $1 \lt q \lt 2$ and $\theta_{1}+\theta_{2}\geq2$, it requires that $n\geq3$ and $q \gt 2/n+1$.

Remark 1.8. For any fixed $R_{0} \gt 0$, let $B_{R_{0}}\times(-R_{0}^{\theta_{1}+\theta_{2}},0]$ substitute for $B_{1}\times(-1,0]$ in Theorems 1.6 and 1.9. Then applying their proofs with a slight modification, we derive that Equations (1.7)–(1.8) still hold with $t_{0}\in(-1/4,0)$, $B_{1/2}\times(-1/4,t_{0}]$ and $B_{1/2}\times(-1/4,0)$ replaced by $t_{0}\in(-R^{\theta_{1}+\theta_{2}}_{0}/4,0)$, $B_{R_{0}/2}\times(-R^{\theta_{1}+\theta_{2}}_{0}/4,t_{0}]$ and $B_{R_{0}/2}\times(-R^{\theta_{2}}_{0}/4,0)$, respectively. A difference lies in that the constant C will depend on R 0, but not on α.

In the case of $\theta_{1}=\theta_{3}=0$, w 1 and w 2 become single power-type weight. Then we have

Theorem 1.9. For $p\geq1$, $n\geq2$, $\theta_{1}=\theta_{3}=0$, $\theta_{2}\geq\theta_{4}=2$, let u be a weak solution of problem (1.5) with $\Omega\times(-T,0]=B_{1}\times(-1,0]$. Then there exist two constants $0 \lt \alpha \lt 1$ and C > 0, both depending only on $n,p,\theta_{2},\lambda,\overline{m},\overline{M},$ such that

(1.8)\begin{align} |u(x,t)-u(y,s)|\leq C\big(|x-y|+\sqrt[\theta_{2}]{|t-s|}\big)^{\alpha} \end{align}

for any $(x,t),(y,s)\in B_{1/2}\times(-1/4,0).$

When $\theta_{1}=\theta_{3}=0$ and $\theta_{2}\geq\theta_{4}=2$, Equation (1.5) will be uniformly parabolic in any domain away from the origin. Then we can directly use the interior Hölder estimates for uniformly parabolic equation to obtain its Hölder regularity in these regions. This, together with Remark 1.8 and Theorem 1.9, leads to the following corollary.

Corollary 1.10. For $p\geq1$, $n\geq2$, $\theta_{1}=\theta_{3}=0$, $\theta_{2}\geq\theta_{4}=2$, let u be a weak solution of problem (1.5). Then u is locally Hölder continuous in $\Omega\times(-T,0)$, that is, for any compact subset $K\subset\Omega\times(-T,0)$, there exist a small constant $0 \lt \alpha=\alpha(n,p,\theta_{2},\lambda,\overline{m},\overline{M}) \lt 1$ and a large constant $C=C(\mathrm{dist}(K,\partial(\Omega\times(-T,0))),\mathrm{dist}(0,\partial\Omega),n,p,\theta_{2},\lambda,\overline{m},\overline{M}) \gt 0$ such that Equation (1.8) holds with K substituting for $B_{1/2}\times(-1/4,0)$.

The rest of this paper is organized as follows. In $\S$ 2, we establish the anisotropic weighted Poincaré type inequality and its corresponding isoperimetric inequality. Then we make use of the De Giorgi truncation method [Reference De Giorgi16] to study the local regularity for solutions to the nonlinear elliptic and parabolic equations with anisotropic weights in $\S$ 3 and $\S$ 4, respectively.

2. Anisotropic weighted Poincaré inequality and its application to the isoperimetric inequality

As pointed out in the introduction, this section is devoted to establishing the anisotropic weighted Poincaré-type inequality. It will be achieved by proving that $w=|x^{\prime}|^{\theta_{1}}|x|^{\theta_{2}}$ is an A q-weight under the condition $(\theta_{1},\theta_{2})\in(\mathcal{A}\cup\mathcal{B})\cap(\mathcal{C}_{q}\cup\mathcal{D}_{q})$, see Theorem 2.6 and Corollary 2.8 below. As a consequence, we derive the isoperimetric inequality of De Giorgi type, which is critical to application for the De Giorgi truncation method in [Reference De Giorgi16].

Denote by ω n the volume of unit ball in $\mathbb{R}^{n}$. In this section, we employ the notation $a\sim b$ to represent that there exists a constant $C=C(n,\theta_{1},\theta_{2}) \gt 0$ such that $\frac{1}{C}b\leq a\leq Cb.$ To begin with, we have

Lemma 2.1. $d\mu:=wdx=|x^{\prime}|^{\theta_{1}}|x|^{\theta_{2}}\,{\rm d}x$ is a Radon measure if $(\theta_{1},\theta_{2})\in\mathcal{A}\cup\mathcal{B}$. Moreover, $\mu(B_{R})\sim R^{n+\theta_{1}+\theta_{2}}$ for any R > 0.

Proof. It suffices to verify that the weight w is a locally integrable function in $\mathbb{R}^{n}$ under these above cases. Specifically, it only needs to prove that for any r > 0, there holds $\mu(B_{r}) \lt \infty$. Recall the following elemental inequalities: for $a,b\geq0$,

(2.1)\begin{align} \begin{cases} \frac{a^{q}+b^{q}}{2}\leq(a+b)^{q}\leq a^{q}+b^{q},&0\leq q\leq1,\\ a^{q}+b^{q}\leq(a+b)^{q}\leq2^{q-1}(a^{q}+b^{q}),&q \gt 1. \end{cases} \end{align}

Step 1. Consider the case when $\theta_{2}\geq0$. Then we have from Equation (2.1) that

\begin{align*} \mu(B_{r})=&2\int_{B_{r}\cap\{x_{n} \gt 0\}}\left|x^{\prime}\right|^{\theta_{1}}\left|x\right|^{\theta_{2}}\,{\rm d}x\sim\int_{B_{r}^{\prime}}\left|x^{\prime}\right|^{\theta_{1}}\,{\rm d}x^{\prime}\int_{0}^{\sqrt{r^{2}-|x^{\prime}|^{2}}}\left(\left|x^{\prime}\right|^{\theta_{2}}+x_{n}^{\theta_{2}}\right)\,{\rm d}x_{n}\notag\\ \sim&\int^{r}_{0}\left(s^{n+\theta_{1}+\theta_{2}-2}\left(r^{2}-s^{2}\right)^{\frac{1}{2}}+s^{n+\theta_{1}-2}\left(r^{2}-s^{2}\right)^{\frac{\theta_{2}+1}{2}}\right)\,{\rm d}s\notag\\ \sim&\,r^{n+\theta_{1}+\theta_{2}}\int^{1}_{0}\left(s^{n+\theta_{1}+\theta_{2}-2}\left(1-s\right)^{\frac{1}{2}}+s^{n+\theta_{1}-2}\left(1-s\right)^{\frac{\theta_{2}+1}{2}}\right)\,{\rm d}s. \end{align*}

Observe that this type of integration is called Beta function. It makes sense if and only if $n+\theta_{1}+\theta_{2}-1 \gt 0$ and $n+\theta_{1}-1 \gt 0$. Then the conclusion is proved in the case when $(\theta_{1},\theta_{2})\in\mathcal{A}$.

Step 2. Consider $\theta_{2} \lt 0$. Then it follows from Equation (2.1) that

(2.2)\begin{align} \mu(B_{r})=&2\int_{B_{r}\cap\{x_{n} \gt 0\}}|x^{\prime}|^{\theta_{1}}|x|^{\theta_{2}}\,{\rm d}x\notag\\ \sim&\int_{B_{r}^{\prime}}|x^{\prime}|^{\theta_{1}}\,{\rm d}x^{\prime}\int^{\sqrt{r^{2}-|x^{\prime}|^{2}}}_{0}\frac{1}{|x^{\prime}|^{-\theta_{2}}+x_{n}^{-\theta_{2}}}\,{\rm d}x_{n}. \end{align}

For the last integration term in Equation (2.2), we further split it as follows:

\begin{align*} I_{1}=&\int_{B^{\prime}_{\frac{1}{\sqrt{2}}r}}|x^{\prime}|^{\theta_{1}}\,{\rm d}x^{\prime}\int^{|x^{\prime}|}_{0}\frac{1}{|x^{\prime}|^{-\theta_{2}}+x_{n}^{-\theta_{2}}}\,{\rm d}x_{n},\notag\\ I_{2}=&\int_{B^{\prime}_{\frac{1}{\sqrt{2}}r}}|x^{\prime}|^{\theta_{1}}\,{\rm d}x^{\prime}\int^{\sqrt{r^{2}-|x^{\prime}|^{2}}}_{|x^{\prime}|}\frac{1}{|x^{\prime}|^{-\theta_{2}}+x_{n}^{-\theta_{2}}}\,{\rm d}x_{n},\notag\\ I_{3}=&\int_{B_{r}^{\prime}\setminus B^{\prime}_{\frac{1}{\sqrt{2}}r}}|x^{\prime}|^{\theta_{1}}\,{\rm d}x^{\prime}\int^{\sqrt{r^{2}-|x^{\prime}|^{2}}}_{0}\frac{1}{|x^{\prime}|^{-\theta_{2}}+x_{n}^{-\theta_{2}}}\,{\rm d}x_{n}. \end{align*}

Observe that $|x^{\prime}|^{-\theta_{2}}\leq|x^{\prime}|^{-\theta_{2}}+x_{n}^{-\theta_{2}}\leq2|x^{\prime}|^{-\theta_{2}}$ if $0\leq x_{n}\leq |x^{\prime}|$. Then for the first term I 1, we have

\begin{align*} I_{1}\sim&\int_{B^{\prime}_{\frac{1}{\sqrt{2}}r}}|x^{\prime}|^{\theta_{1}}\,{\rm d}x^{\prime}\int^{|x^{\prime}|}_{0}|x^{\prime}|^{\theta_{2}}\,{\rm d}x_{n}\sim\int^{\frac{r}{\sqrt{2}}}_{0}s^{n+\theta_{1}+\theta_{2}-1}\,{\rm d}s\notag\\ \sim&r^{n+\theta_{1}+\theta_{2}}\int^{1}_{0}s^{n+\theta_{1}+\theta_{2}-1}\,{\rm d}s. \end{align*}

This integration makes sense iff $n+\theta_{1}+\theta_{2} \gt 0$.

With regard to the second term I 2, we divide it into two cases to discuss as follow.

Case 1. If $\theta_{2}=-1$, then

\begin{align*} I_{2}\sim&\int_{B^{\prime}_{\frac{1}{\sqrt{2}}r}}|x^{\prime}|^{\theta_{1}}\,{\rm d}x^{\prime}\int^{\sqrt{r^{2}-|x^{\prime}|^{2}}}_{|x^{\prime}|}x_{n}^{-1}\,{\rm d}x_{n}\sim\int^{\frac{r}{\sqrt{2}}}_{0}s^{n+\theta_{1}-2}\,\ln\frac{\sqrt{r^{2}-s^{2}}}{s}\,{\rm d}s\notag\\ \sim&\,r^{n+\theta_{1}-1}\int^{\frac{1}{\sqrt{2}}}_{0}s^{n+\theta_{1}-2}\,\ln\frac{\sqrt{1-s^{2}}}{s}\,{\rm d}s\sim-r^{n+\theta_{1}-1}\int^{\frac{1}{\sqrt{2}}}_{0}s^{n+\theta_{1}-2}\,\ln s\,{\rm d}s\notag\\ \sim&-r^{n+\theta_{1}+\theta_{2}}\bigg(s^{n+\theta_{1}+\theta_{2}}\,\ln s\Big|^{\frac{1}{\sqrt{2}}}_{0}-\int^{\frac{1}{\sqrt{2}}}_{0}s^{n+\theta_{1}+\theta_{2}-1}\,{\rm d}s\bigg). \end{align*}

It makes sense iff $n+\theta_{1}+\theta_{2} \gt 0$.

Case 2. If $\theta_{2}\neq-1$, then

\begin{align*} I_{2}\sim&\int_{B^{\prime}_{\frac{1}{\sqrt{2}}r}}|x^{\prime}|^{\theta_{1}}\,{\rm d}x^{\prime}\int^{\sqrt{r^{2}-|x^{\prime}|^{2}}}_{|x^{\prime}|}x_{n}^{\theta_{2}}\,{\rm d}x_{n}\notag\\ \sim&\int^{\frac{1}{\sqrt{2}}r}_{0}s^{n+\theta_{1}-2}\left(\left(r^{2}-s^{2}\right)^{\frac{\theta_{2}+1}{2}}-s^{\theta_{2}+1}\right)\,{\rm d}s\notag\\ \sim&\,r^{n+\theta_{1}+\theta_{2}}\int^{\frac{1}{\sqrt{2}}}_{0}s^{n+\theta_{1}-2}\left(\left(1-s^{2}\right)^{\frac{\theta_{2}+1}{2}}-s^{\theta_{2}+1}\right)\,{\rm d}s. \end{align*}

Note that $\min\{1,2^{-\frac{\theta_{2}+1}{2}}\}\leq|(1-s^{2})^{\frac{\theta_{2}+1}{2}}|\leq\max\{1,2^{-\frac{\theta_{2}+1}{2}}\}$ in $[0,\frac{1}{\sqrt{2}}].$ Then it makes sense iff $\theta_{1} \gt -(n-1)$ and $n+\theta_{1}+\theta_{2} \gt 0$.

The last term I 3 remains to be analyzed. Note that $|x^{\prime}|\geq\sqrt{r^{2}-|x^{\prime}|^{2}}\geq x_{n}$ if $\frac{1}{\sqrt{2}}r\leq|x^{\prime}|\leq r$ and $0\leq x_{n}\leq\sqrt{r^{2}-|x^{\prime}|^{2}}$. Then we deduce

\begin{align*} I_{3}\sim&\int_{B_{r}^{\prime}\setminus B^{\prime}_{\frac{1}{\sqrt{2}}r}}|x^{\prime}|^{\theta_{1}}\,{\rm d}x^{\prime}\int^{\sqrt{r^{2}-|x^{\prime}|^{2}}}_{0}|x^{\prime}|^{\theta_{2}}\,{\rm d}x_{n}\notag\\ \sim&\int_{\frac{1}{\sqrt{2}}r}^{r}s^{n+\theta_{1}+\theta_{2}-2}\sqrt{r^{2}-s^{2}}\,{\rm d}s\sim r^{n+\theta_{1}+\theta_{2}}\int_{\frac{1}{\sqrt{2}}}^{1}s^{n+\theta_{1}+\theta_{2}-2}\sqrt{1-s^{2}}\,{\rm d}s\notag\\ \sim&r^{n+\theta_{1}+\theta_{2}}, \end{align*}

where we used the fact that the integrand $s^{n+\theta_{1}+\theta_{2}}\sqrt{1-s^{2}}$ has no singular point in $[\frac{1}{\sqrt{2}},1].$ Consequently, combining these above facts, we obtain that if $(\theta_{1},\theta_{2})\in\mathcal{B}$, then ${\rm d}\mu$ is a Radon measure. The proof is complete.

Definition 2.2. A Radon measure ${\rm d}\mu$ is called doubling if there exists some constant $0 \lt C \lt \infty$ such that $\mu(B_{2R}(\bar{x}))\leq C\mu(B_{R}(\bar{x}))$ for any $\bar{x}\in\mathbb{R}^{n}$ and R > 0.

Theorem 2.3. The Radon measure ${\rm d}\mu=wdx$ is doubling if $(\theta_{1},\theta_{2})\in\mathcal{A}\cup\mathcal{B}$.

Remark 2.4. $w=|x^{\prime}|^{\theta_{1}}|x|^{\theta_{2}}$ degenerates to be an isotropic weight if $\theta_{1}=0$. In this case, it is doubling if $\theta_{2} \gt -n$. Its proof is simple and direct, see pages 505–506 in [Reference Grafakos24] for more details. By contrast, it will involve complex analyses, computations and discussions if $\theta_{1}\neq0$.

Proof. For any $\bar{x}\in\mathbb{R}^{n}$ and R > 0, we divide all balls $B_{R}(\bar{x})$ into two types as follows: the first type satisfies $|\bar{x}|\geq3R$ and the second type satisfies $|\bar{x}| \lt 3R$.

Step 1. Consider the case when $|\bar{x}|\geq 3R$. Observe that

(2.3)\begin{align} \int_{B_{R}(\bar{x})}|x^{\prime}|^{\theta_{1}}|x|^{\theta_{2}}\,{\rm d}x\geq \begin{cases} (|\bar{x}|-R)^{\theta_{2}}\int_{B_{R}(\bar{x})}|x^{\prime}|^{\theta_{1}}\,{\rm d}x,&\text{if }\theta_{2}\geq0,\\ (|\bar{x}|+R)^{\theta_{2}}\int_{B_{R}(\bar{x})}|x^{\prime}|^{\theta_{1}}\,{\rm d}x,&\text{if }\theta_{2} \lt 0, \end{cases} \end{align}

and

(2.4)\begin{align} \int_{B_{2R}(\bar{x})}|x^{\prime}|^{\theta_{1}}|x|^{\theta_{2}}\,{\rm d}x\leq \begin{cases} (|\bar{x}|+2R)^{\theta_{2}}\int_{B_{2R}(\bar{x})}|x^{\prime}|^{\theta_{1}}\,{\rm d}x,&\text{if }\theta_{2}\geq0,\\ (|\bar{x}|-2R)^{\theta_{2}}\int_{B_{2R}(\bar{x})}|x^{\prime}|^{\theta_{1}}\,{\rm d}x,&\text{if }\theta_{2} \lt 0. \end{cases} \end{align}

On one hand, we have

\begin{align*} \int_{B_{R}(\bar{x})}|x^{\prime}|^{\theta_{1}}\,{\rm d}x=&2\int_{B_{R}^{\prime}(\bar{x}^{\prime})}|x^{\prime}|^{\theta_{1}}\sqrt{R^{2}-|x^{\prime}-\bar{x}^{\prime}|^{2}}\,{\rm d}x^{\prime}\geq\sqrt{3}R\int_{B^{\prime}_{R/2}(\bar{x}^{\prime})}|x^{\prime}|^{\theta_{1}}\,{\rm d}x^{\prime}. \end{align*}

Observe that $|x^{\prime}|^{\theta_{1}}$ increases radially if $\theta_{1}\geq0$, while it decreases radially for $\theta_{1} \lt 0$. Then we obtain that

  1. (i) for $|\bar{x}^{\prime}|\geq\frac{3}{2}R$, then

    \begin{align*} \int_{B^{\prime}_{R/2}(\bar{x}^{\prime})}|x^{\prime}|^{\theta_{1}}\,{\rm d}x^{\prime}\geq&\omega_{n-1}\left(\frac{R}{2}\right)^{n-1} \begin{cases} \left(|\bar{x}^{\prime}|-R/2\right)^{\theta_{1}},&\text{if }\theta_{1}\geq0,\\ \left(|\bar{x}^{\prime}|+R/2\right)^{\theta_{1}},&\text{if }\theta_{1} \lt 0; \end{cases} \end{align*}
  2. (ii) for $|\bar{x}^{\prime}| \lt \frac{3}{2}R,$ then

    \begin{align*} \int_{B^{\prime}_{R/2}(\bar{x}^{\prime})}|x^{\prime}|^{\theta_{1}}\,{\rm d}x^{\prime}\geq& \begin{cases} \int_{B^{\prime}_{R/2}(0^{\prime})}|x^{\prime}|^{\theta_{1}}\,{\rm d}x^{\prime},&\text{if }\theta_{1}\geq0,\\ \int_{B^{\prime}_{R/2}\left(\frac{3}{2}R\frac{\bar{x}^{\prime}}{|\bar{x}^{\prime}|}\right)}|x^{\prime}|^{\theta_{1}}\,{\rm d}x^{\prime},&\text{if }\theta_{1} \lt 0 \end{cases}\notag\\ \geq&\omega_{n-1}R^{n-1+\theta_{1}} \begin{cases} \frac{n-1}{2^{n-1+\theta_{1}}(n-1+\theta_{1})},&\text{if }\theta_{1}\geq0,\\ \frac{1}{2^{n-1}},&\text{if }\theta_{1} \lt 0. \end{cases} \end{align*}

On the other hand, we have

(2.5)\begin{align} \int_{B_{2R}(\bar{x})}|x^{\prime}|^{\theta_{1}}\,{\rm d}x=&2\int_{B_{2R}^{\prime}(\bar{x}^{\prime})}|x^{\prime}|^{\theta_{1}}\sqrt{4R^{2}-|x^{\prime}-\bar{x}^{\prime}|^{2}}\,{\rm d}x^{\prime}\notag\\ \leq&4R\int_{B^{\prime}_{2R}(\bar{x}^{\prime})}|x^{\prime}|^{\theta_{1}}\,{\rm d}x^{\prime}. \end{align}

By the same argument as above, we deduce that

  1. (1) for $|\bar{x}^{\prime}|\geq3R$, then

    (2.6)\begin{align} \int_{B^{\prime}_{2R}(\bar{x}^{\prime})}|x^{\prime}|^{\theta_{1}}\,{\rm d}x^{\prime}\leq&\omega_{n-1}(2R)^{n-1} \begin{cases} \left(|\bar{x}^{\prime}|+2R\right)^{\theta_{1}},&\text{if }\theta_{1}\geq0,\\ \left(|\bar{x}^{\prime}|-2R\right)^{\theta_{1}},&\text{if }\theta_{1} \lt 0; \end{cases} \end{align}
  2. (2) for $|\bar{x}^{\prime}| \lt 3R$, then

    (2.7)\begin{align} \int_{B^{\prime}_{2R}(\bar{x}^{\prime})}|x^{\prime}|^{\theta_{1}}\,{\rm d}x^{\prime}\leq&\int_{B^{\prime}_{5R}(0^{\prime})}|x^{\prime}|^{\theta_{1}}\,{\rm d}x^{\prime}=\frac{(n-1)\omega_{n-1}}{n-1+\theta_{1}}(5R)^{n-1+\theta_{1}}. \end{align}

    Note that if $\theta_{1}\geq0$, then

    \begin{align*} \begin{cases} |\bar{x}^{\prime}|+2R\leq2(|\bar{x}^{\prime}|-R/2),&\text{for }|\bar{x}^{\prime}|\geq3R,\\ R \lt |\bar{x}^{\prime}|-R/2 \lt 5R/2,&\text{for }3R/2 \lt |\bar{x}^{\prime}| \lt 3R, \end{cases} \end{align*}
    while, if $\theta_{1} \lt 0$,
    \begin{align*} \begin{cases} |\bar{x}^{\prime}|-2R\geq\frac{2}{7}(|\bar{x}^{\prime}|+R/2),&\text{for }|\bar{x}^{\prime}|\geq3R,\\ 2R \lt |\bar{x}^{\prime}|+R/2 \lt 7R/2,&\text{for }3R/2 \lt |\bar{x}^{\prime}| \lt 3R. \end{cases} \end{align*}

    Then combining these above facts, we obtain

    (2.8)\begin{align} \int_{B_{2R}(\bar{x})}|x^{\prime}|^{\theta_{1}}\,{\rm d}x\leq C(n,\theta_{1},\theta_{2})\int_{B_{R}(\bar{x})}|x^{\prime}|^{\theta_{1}}\,{\rm d}x. \end{align}

    Since $|\bar{x}|\geq3R$, then $|\bar{x}|+2R\leq4(|\bar{x}-R|)$ and $|\bar{x}|+R\leq4(|\bar{x}|-2R)$. This, in combination with Equations (2.3)–(2.8), reads that for $|\bar{x}|\geq3R,$

    (2.9)\begin{align} \int_{B_{2R}(\bar{x})}|x^{\prime}|^{\theta_{1}}|x|^{\theta_{2}}\,{\rm d}x\leq C(n,\theta_{1},\theta_{2})\int_{B_{R}(\bar{x})}|x^{\prime}|^{\theta_{1}}|x|^{\theta_{2}}\,{\rm d}x. \end{align}

Step 2. Let $|\bar{x}| \lt 3R$. Then we have

\begin{align*} \int_{B_{2R}(\bar{x})}|x^{\prime}|^{\theta_{1}}|x|^{\theta_{2}}\,{\rm d}x\leq \int_{B_{5R}(0)}|x^{\prime}|^{\theta_{1}}|x|^{\theta_{2}}\,{\rm d}x\leq C(n,\theta_{1},\theta_{2})R^{n+\theta_{1}+\theta_{2}}. \end{align*}

First, if $\theta_{1} \lt 0$, then

\begin{align*} \int_{B_{R}(\bar{x})}|x^{\prime}|^{\theta_{1}}|x|^{\theta_{2}}\geq&\int_{B_{R}(\bar{x})}|x|^{\theta_{1}+\theta_{2}}\,{\rm d}x\notag\\ \geq& \begin{cases} \int_{B_{R}(0)}|x|^{\theta_{1}+\theta_{2}}\,{\rm d}x,&\text{if }\theta_{2}\geq-\theta_{1},\\ \int_{B_{R}\left(3R\frac{\bar{x}}{|\bar{x}|}\right)}|x|^{\theta_{1}+\theta_{2}}\,{\rm d}x,&\text{if } \theta_{2} \lt -\theta_{1} \end{cases}\notag\\ \geq&R^{n+\theta_{1}+\theta_{2}} \begin{cases} \frac{n\omega_{n}}{n+\theta_{1}+\theta_{2}},&\text{if }\theta_{2}\geq-\theta_{1},\\ 2^{\theta_{1}+\theta_{2}}\omega_{n},&\text{if } \theta_{2} \lt -\theta_{1}, \end{cases} \end{align*}

where we utilized the fact that $|x|^{\theta_{1}+\theta_{2}}$ is radially increasing if $\theta_{1}+\theta_{2}\geq0$ and radially decreasing if $\theta_{1}+\theta_{2} \lt 0$.

Second, if $\theta_{1}\geq0$, we discuss as follows:

  1. (i) for $\theta_{2}\geq0$, similarly as before, we have

    \begin{align*} \int_{B_{R}(\bar{x})}|x^{\prime}|^{\theta_{1}}|x|^{\theta_{2}}\,{\rm d}x\geq&\int_{B_{R}(\bar{x})}|x^{\prime}|^{\theta_{1}+\theta_{2}}\,{\rm d}x\geq\sqrt{3}R\int_{B^{\prime}_{R}(\bar{x}^{\prime})}|x^{\prime}|^{\theta_{1}+\theta_{2}}\,{\rm d}x^{\prime}\notag\\ \geq&\int_{B_{R}^{\prime}(0^{\prime})}|x^{\prime}|^{\theta_{1}+\theta_{2}}\,{\rm d}x^{\prime}=\frac{(n-1)\omega_{n-1}}{n+\theta_{1}+\theta_{2}-1}R^{n+\theta_{1}+\theta_{2}}; \end{align*}
  2. (ii) for $\theta_{2} \lt 0$, then

    \begin{align*} \int_{B_{R}(\bar{x})}|x^{\prime}|^{\theta_{1}}|x|^{\theta_{2}}\,{\rm d}x=&\int_{B_{R}(\bar{x})}|x^{\prime}|^{\theta_{1}+\theta_{2}}\left(\frac{|x^{\prime}|}{|x|}\right)^{-\theta_{2}}\,{\rm d}x\notag\\ \geq&\frac{8^{\theta_{2}}}{2^{\theta_{1}+\theta_{2}}}R^{\theta_{1}+\theta_{2}}\int_{B_{R}(\bar{x})\cap\{|x^{\prime}| \gt R/2\}}\,{\rm d}x\notag\\ \geq&\frac{8^{\theta_{2}}}{2^{\theta_{1}+\theta_{2}}}R^{\theta_{1}+\theta_{2}}\int_{B_{R}(0)\cap\{|x^{\prime}| \gt R/2\}}\,{\rm d}x\notag\\ \geq&\frac{8^{\theta_{2}}(\omega_{n}-2^{2-n}\omega_{n-1})}{2^{\theta_{1}+\theta_{2}}}R^{n+\theta_{1}+\theta_{2}}. \end{align*}

    Then combining these aforementioned facts, we obtain that Equation (2.9) also holds if $|\bar{x}| \lt 3R.$ The proof is complete.

Definition 2.5. Let $1 \lt q \lt \infty$. We say that w is an A q-weight, if there is a positive constant $C=C(n,q,w)$ such that

\begin{align*} \rlap{\mbox{--}}\kern-3.5pt\int_{B}w\,{\rm d}x\left(\rlap{\mbox{--}}\kern-3.5pt\int_{B}w^{-\frac{1}{q-1}}\,{\rm d}x\right)^{q-1}\leq C(n,q,w),\quad\mathrm{with}\;\rlap{\mbox{--}}\kern-3.5pt\int_{B}=\frac{1}{|B|}\int_{B} \end{align*}

for any ball B in $\mathbb{R}^{n}$.

Theorem 2.6. Let $1 \lt q \lt \infty$. If $(\theta_{1},\theta_{2})\in(\mathcal{A}\cup\mathcal{B})\cap(\mathcal{C}_{q}\cup\mathcal{D}_{q})$, then $w=|x^{\prime}|^{\theta_{1}}|x|^{\theta_{2}}$ is an A q-weight.

Remark 2.7. From Theorems 2.3 and 2.6, we see that the Radon measure ${\rm d}\mu=w\,{\rm d}x=|x^{\prime}|^{\theta_{1}}|x|^{\theta_{2}}\,{\rm d}x$ is doubling on a larger range $(\theta_{1},\theta_{2})\in\mathcal{A}\cup\mathcal{B}$. This implies that when $(\theta_{1},\theta_{2})\in(\mathcal{A}\cup\mathcal{B})\setminus(\mathcal{C}_{q}\cup\mathcal{D}_{q})$, the weight $w=|x^{\prime}|^{\theta_{1}}|x|^{\theta_{2}}$ provides an example of a doubling measure but is not in A q.

Proof. For $1 \lt q \lt \infty,$ according to the definition of A q-weight, it needs to verify the following inequality:

(2.10)\begin{align} \rlap{\mbox{--}}\kern-3.5pt\int_{B}w\,{\rm d}x\left(\rlap{\mbox{--}}\kern-3.5pt\int_{B}w^{-\frac{1}{q-1}}\,{\rm d}x\right)^{q-1}\leq C(n,q,\theta_{1},\theta_{2}) \end{align}

for any ball $B\subset\mathbb{R}^{n}$. For any R > 0 and $\bar{x}\in\mathbb{R}^{n}$, the ball $B_{R}(\bar{x})$ must belong to one of the following two types: $|\bar{x}|\geq3R$ and $|\bar{x}| \lt 3R$. On the one hand, if $|\bar{x}|\geq3R$, then we have

(2.11)\begin{align} \frac{2}{3}|\bar{x}|\leq|\bar{x}|-R\leq|x|\leq|\bar{x}|+R\leq\frac{4}{3}|\bar{x}|, \quad\text{for }x\in B_{R}(\bar{x}). \end{align}

Applying Equations (2.4)–(2.7) with $B_{2R}(\bar{x})$ and $B^{\prime}_{2R}(\bar{x}^{\prime})$ replaced by $B_{R}(\bar{x})$ and $B^{\prime}_{R}(\bar{x}^{\prime})$, it follows from Equation (2.11) that

\begin{align*} \int_{B_{R}(\bar{x})}|x^{\prime}|^{\theta_{1}}|x|^{\theta_{2}}\,{\rm d}x\leq&C(\theta_{2})|\bar{x}|^{\theta_{2}}\int_{B_{R}(\bar{x})}|x^{\prime}|^{\theta_{1}}\,{\rm d}x\leq C(\theta_{2})|\bar{x}|^{\theta_{2}+1}\int_{B_{R}^{\prime}(\bar{x}^{\prime})}|x^{\prime}|^{\theta_{1}}\,{\rm d}x^{\prime}\notag\\ \leq&C(n,\theta_{1},\theta_{2})|\bar{x}|^{n+\theta_{1}+\theta_{2}} \end{align*}

and

\begin{align*} \int_{B_{R}(\bar{x})}|x^{\prime}|^{-\frac{\theta_{1}}{q-1}}|x|^{-\frac{\theta_{2}}{q-1}}\,{\rm d}x\leq&C(q,\theta_{2})|\bar{x}|^{-\frac{\theta_{2}}{q-1}}\int_{B_{R}(\bar{x})}|x^{\prime}|^{-\frac{\theta_{1}}{q-1}}\,{\rm d}x\notag\\ \leq& C(q,\theta_{2})|\bar{x}|^{-\frac{\theta_{2}}{q-1}+1}\int_{B_{R}^{\prime}(\bar{x}^{\prime})}|x^{\prime}|^{-\frac{\theta_{1}}{q-1}}\,{\rm d}x^{\prime}\notag\\ \leq& C(n,q,\theta_{1},\theta_{2})|\bar{x}|^{n-\frac{\theta_{1}+\theta_{2}}{p-1}}, \end{align*}

where we require that $-(n-1) \lt \theta_{1} \lt (n-1)(q-1)$ and $\theta_{2}\in\mathbb{R}$. Combining these two relations, we obtain that Equation (2.10) holds in the case of $|\bar{x}|\geq3R$.

On the other hand, if $|\bar{x}| \lt 3R$, we have $|x|\leq4R$ for $x\in B_{R}(\bar{x})$. Therefore, it follows from Lemma 2.1 that

\begin{align*} \rlap{\mbox{--}}\kern-3.5pt\int_{B_{R}(\bar{x})}|x^{\prime}|^{\theta_{1}}|x|^{\theta_{2}}\,{\rm d}x\leq 4^{n}\rlap{\mbox{--}}\kern-3.5pt\int_{B_{4R}(0)}|x^{\prime}|^{\theta_{1}}|x|^{\theta_{2}}\,{\rm d}x\leq C(n,\theta_{1},\theta_{2})R^{\theta_{1}+\theta_{2}} \end{align*}

and

\begin{align*} \left(\rlap{\mbox{--}}\kern-3.5pt\int_{B_{R}(\bar{x})}|x^{\prime}|^{-\frac{\theta_{1}}{q-1}}|x|^{-\frac{\theta_{2}}{q-1}}\,{\rm d}x\right)^{q-1}\leq& 4^{n(q-1)}\left(\rlap{\mbox{--}}\kern-3.5pt\int_{B_{4R}(0)}|x^{\prime}|^{-\frac{\theta_{1}}{q-1}}|x|^{-\frac{\theta_{2}}{q-1}}\,{\rm d}x\right)^{q-1}\notag\\ \leq& C(n,q,\theta_{1},\theta_{2})R^{-\theta_{1}-\theta_{2}}, \end{align*}

where these two relations hold if $(\theta_{1},\theta_{2}),(-\frac{\theta_{1}}{q-1},-\frac{\theta_{2}}{q-1})\in\mathcal{A}\cup\mathcal{B}$, that is, $(\theta_{1},\theta_{2})\in(\mathcal{A}\cup\mathcal{B})\cap(\mathcal{C}_{q}\cup\mathcal{D}_{q})$. Therefore, Equation (2.11) holds for any $B\subset\mathbb{R}^{n}$. The proof is complete.

Denote ${\rm d}\mu:=w\,{\rm d}x=|x^{\prime}|^{\theta_{1}}|x|^{\theta_{2}}\,{\rm d}x.$ Combining Theorem 15.21 and Corollary 15.35 in [Reference Heinonen, Kilpeläinen and Martio26] and Theorem 2.6 above, we obtain the following anisotropic weighted Poincaré inequality.

Corollary 2.8. For $n\geq2$ and $1 \lt q \lt \infty$, let $(\theta_{1},\theta_{2})\in[(\mathcal{A}\cup\mathcal{B})\cap(C_{q}\cup\mathcal{D}_{q})]\cup\{\theta_{1}=0,\,\theta_{2} \gt -n\}$. Then for any $B:=B_{R}(\bar{x})\subset\mathbb{R}^{n}$, R > 0 and $\varphi\in W^{1,q}(B,w)$,

(2.12)\begin{align} \int_{B}|\varphi-\varphi_{B}|^{q}\,{\rm d}\mu\leq C(n,q,\theta_{1},\theta_{2})R^{q}\int_{B}|\nabla\varphi|^{q}\,{\rm d}\mu, \end{align}

where $\varphi_{B}=\frac{1}{\mu(B)}\int_{B}\varphi \,{\rm d}\mu.$

Remark 2.9. It is worth emphasizing that according to Corollary 15.35 in [Reference Heinonen, Kilpeläinen and Martio26], Equation (2.12) holds for any $(\theta_{1},\theta_{2})\in\{\theta_{1}=0,\,\theta_{2} \gt -n\}$ and $1 \lt q \lt \infty$. This conclusion is very strong, which is achieved by combining the theories of A q-weights and quasiconformal mappings; see Chapter 15 of [Reference Heinonen, Kilpeläinen and Martio26] for further details.

Making use of the anisotropic weighted Poincaré inequality in Corollary 2.8, we can establish the corresponding weighted isoperimetric inequality of De Giorgi type as follows.

Proposition 2.10. For $n\geq2$ and $1 \lt q \lt \infty$, let $(\theta_{1},\theta_{2})\in[(\mathcal{A}\cup\mathcal{B})\cap(C_{q}\cup\mathcal{D}_{q})]\cup\{\theta_{1}=0,\,\theta_{2} \gt -n\}$. Then for any R > 0, l > k and $u\in W^{1,q}(B_{R},w)$,

(2.13)\begin{align} &(l-k)^{q}\left(\int_{\{u\geq l\}\cap B_{R}}{\rm d}\mu\right)^{q}\int_{\{u\leq k\}\cap B_{R}}{\rm d}\mu\notag\\ &\leq C(n,q,\theta_{1},\theta_{2})R^{q(n+\theta_{1}+\theta_{2}+1)}\int_{\{k \lt u \lt l\}\cap B_{R}}|\nabla u|^{q}\,{\rm d}\mu \end{align}

and

(2.14)\begin{align} &(l-k)^{q}\left(\int_{\{u\leq k\}\cap B_{R}}{\rm d}\mu\right)^{q}\int_{\{u\geq l\}\cap B_{R}}{\rm d}\mu\notag\\ &\leq C(n,q,\theta_{1},\theta_{2})R^{q(n+\theta_{1}+\theta_{2}+1)}\int_{\{k \lt u \lt l\}\cap B_{R}}|\nabla u|^{q}\,{\rm d}\mu, \end{align}

where ${\rm d}\mu=w\,{\rm d}x=|x^{\prime}|^{\theta_{1}}|x|^{\theta_{2}}\,{\rm d}x.$

Remark 2.11. Since the index q > 1, we have to establish two isoperimetric inequalities in Proposition 2.10, which are used to capture the decaying rates of the distribution function in Lemmas 3.4 and 4.5 below. Meanwhile, it also causes more complex calculations in the proofs of Lemmas 3.4 and 4.5.

Proof. Step 1. Set

\begin{align*} u_{1}=\inf\{u,l\}-\inf\{u,k\},\quad\bar{u}_{1}=\frac{1}{\mu(B_{R})}\int_{B_{R}}u_{1}\,{\rm d}\mu. \end{align*}

First, we have

\begin{align*} \int_{\{u_{1}=0\}\cap B_{R}}\bar{u}_{1}^{q}\,{\rm d}\mu=&\frac{1}{(\mu(B_{R}))^{q}}\left(\int_{B_{R}}u_{1}\,{\rm d}\mu\right)^{q}\int_{\{u\leq k\}\cap B_{R}}{\rm d}\mu\notag\\ \geq&\frac{C(l-k)^{q}}{R^{q(n+\theta_{1}+\theta_{2})}}\left(\int_{\{u\geq l\}\cap B_{R}}{\rm d}\mu\right)^{q}\int_{\{u\leq k\}\cap B_{R}}{\rm d}\mu. \end{align*}

Second, it follows from Corollary 2.8 that

\begin{align*} \int_{\{u_{1}=0\}\cap B_{R}}|\bar{u}_{1}|^{q}\,{\rm d}\mu\leq&\int_{B_{R}}|u_{1}-\bar{u}_{1}|^{q}\,{\rm d}\mu\leq CR^{q}\int_{B_{R}}|\nabla u_{1}|^{q}\,{\rm d}\mu\notag\\ =&CR^{q}\int_{\{k \lt u \lt l\}\cap B_{R}}|\nabla u|^{q}\,{\rm d}\mu. \end{align*}

The proof of Equation (2.13) is finished.

Step 2. Denote

\begin{align*} u_{2}=\sup\{u,l\}-\sup\{u,k\},\quad\bar{u}_{2}=\frac{\int_{B_{R}}u_{2}\,{\rm d}\mu}{|B_{R}|_{\mu}}. \end{align*}

By the same argument as before, we have

\begin{align*} \int_{\{u_{2}=0\}\cap B_{R}}\bar{u}_{2}^{q}\,{\rm d}\mu=&\frac{1}{(\mu(B_{R}))^{q}}\left(\int_{B_{R}}u_{2}\,{\rm d}\mu\right)^{q}\int_{\{u\geq l\}\cap B_{R}}{\rm d}\mu\notag\\ \geq&\frac{C(l-k)^{q}}{R^{q(n+\theta_{1}+\theta_{2})}}\left(\int_{\{u\leq k\}\cap B_{R}}{\rm d}\mu\right)^{q}\int_{\{u\geq l\}\cap B_{R}}{\rm d}\mu \end{align*}

and

\begin{align*} \int_{\{u_{2}=0\}\cap B_{R}}|\bar{u}_{2}|^{q}\,{\rm d}\mu\leq&\int_{B_{R}}|u_{2}-\bar{u}_{2}|^{q}\,{\rm d}\mu\leq CR^{q}\int_{B_{R}}|\nabla u_{2}|^{q}\,{\rm d}\mu\notag\\ =&CR^{q}\int_{\{k \lt u \lt l\}\cap B_{R}}|\nabla u|^{q}\,{\rm d}\mu. \end{align*}

The proof is complete.

3. Regularity for solutions to degenerate elliptic equations with anisotropic weights

Throughout this section, denote ${\rm d}\mu:=w\,{\rm d}x=|x^{\prime}|^{\theta_{1}}|x|^{\theta_{2}}\,{\rm d}x.$ The first step is to establish a Caccioppoli inequality for the truncated solution.

Lemma 3.1. Let u be the solution of problem (1.1). Then for any non-negative $\eta\in C^{\infty}_{0}(B_{R}(x_{0}))$ with any $B_{R}(x_{0})\subset B_{1}$,

\begin{align*} \int_{B_{R}(x_{0})}|\nabla(v\eta)|^{p}w\,{\rm d}x\leq C(n,p,\lambda)\int_{B_{R}(x_{0})}|\nabla\eta|^{p}|v|^{p}w\,{\rm d}x, \end{align*}

where $v=(u-k)^{+}$ or $(u-k)^{-}$ with $k\geq0$.

Proof. First, pick test function $\varphi=v\eta^{p}$ if $v=(u-k)^{+}$. Since

\begin{align*} 0=\int_{B_{R}(x_{0})}Aw|\nabla u|^{p-2}\nabla u\cdot\nabla\varphi=\int_{B_{R}(x_{0})}Aw|\nabla v|^{p-2}\nabla v\cdot\nabla\varphi, \end{align*}

then it follows from Young’s inequality that

(3.1)\begin{align} &\lambda\int_{B_{R}(x_{0})}|\nabla v|^{p}\eta^{p}w\,{\rm d}x\notag\\ &\leq\int_{B_{R}(x_{0})}Aw|\nabla v|^{p}\eta^{p}=-p\int_{B_{R}(x_{0})}Aw|\nabla v|^{p-2}\nabla v\cdot\nabla\eta v\eta^{p-1}\notag\\ &\leq p\lambda\int_{B_{R}(x_{0})}|\eta\nabla v|^{p-1}|v\nabla\eta|w\,{\rm d}x\notag\\ &\leq\frac{\lambda}{2}\int_{B_{R}(x_{0})}|\nabla v|^{p}\eta^{p}w\,{\rm d}x+C\int_{B_{R}(x_{0})}|v|^{p}|\nabla\eta|^{p}w\,{\rm d}x, \end{align}

which yields that

\begin{align*} \int_{B_{R}(x_{0})}|\nabla(v\eta)|^{p}w\,{\rm d}x\leq&2^{p-1}\int_{B_{R}(x_{0})}\big(|\eta\nabla v|^{p}+|v\nabla\eta|^{p}\big)w\,{\rm d}x\notag\\ \leq& C\int_{B_{R}(x_{0})}|\nabla\eta|^{p}|v|^{p}w\,{\rm d}x. \end{align*}

Second, choose test function $\varphi=-v\eta^{2}$ if $v=(u-k)^{-}$. Then we have

\begin{align*} 0=\int_{B_{R}(x_{0})}Aw|\nabla u|^{p-2}\nabla u\cdot\nabla\varphi=-\int_{B_{R}(x_{0})}Aw|\nabla v|^{p-2}\nabla v\cdot\nabla\varphi. \end{align*}

Therefore, in exactly the same way to Equation (3.1), we obtain that Lemma 3.1 holds.

We now improve the oscillation of the solution u in a small domain provided that u is small on a large portion of a larger domain.

Lemma 3.2. Assume that $n\geq2$, $(\theta_{1},\theta_{2})\in(\mathcal{A}\cup\mathcal{B})\cap\mathcal{F}$, $1 \lt p \lt n+\theta_{1}+\theta_{2}$. For $R\in(0,1)$, let $0\leq m\leq\inf\limits_{B_{R}}u\leq\sup\limits_{B_{R}}u\leq M\leq\overline{M}$. Then there exists a small constant $0 \lt \tau_{0}=\tau_{0}(n,p,\theta_{1},\theta_{2},\lambda) \lt 1$ such that for any ɛ > 0 and $0 \lt \tau \lt \tau_{0}$,

  1. (a) if

    (3.2)\begin{align} |\{x\in B_{R}: u \gt M-\varepsilon\}|_{\mu}\leq\tau|B_{R}|_{\mu}, \end{align}
    then
    (3.3)\begin{align} u\leq M-\frac{\varepsilon}{2},\quad\text{for}\quad x\in B_{R/2}; \end{align}
  2. (b) if

    \begin{align*} |\{x\in B_{R}: u \lt m+\varepsilon\}|_{\mu}\leq\tau|B_{R}|_{\mu}, \end{align*}
    then
    (3.4)\begin{align} u\geq m-\frac{\varepsilon}{2},\quad\text{for}\quad x\in B_{R/2}. \end{align}

Remark 3.3. The assumed condition in Equation (3.2) is natural, since

\begin{align*} |\{x\in B_{R}:\,u \gt M-\varepsilon\}|_{\mu}\rightarrow0,\quad\mathrm{as}\; \varepsilon\rightarrow0. \end{align*}

This fact also implies that the value of τ can be chosen to satisfy that $\tau\rightarrow0$, as $\varepsilon\rightarrow0$. Then the key to applying Lemma 3.2 lies in making clear the dependency between τ and ɛ in condition (3.2). The purpose will be achieved by establishing the explicit decaying estimates in terms of the distribution function of u in Lemma 3.4 below.

Proof. Step 1. For ɛ > 0 and $i=0,1,2,\ldots,$ let

\begin{align*} r_{i}=\frac{R}{2}+\frac{R}{2^{i+1}},\quad k_{i}=M-\varepsilon+\frac{\varepsilon}{2}(1-2^{-i}). \end{align*}

Take a cutoff function $\eta_{i}\in C_{0}^{\infty}(B_{r_{i}})$, satisfying that $\eta_{i}=1$ in $B_{r_{i+1}}$, $0\leq\eta_{i}\leq1$, $|\nabla\eta_{i}|\leq C(r_{i}-r_{i+1})^{-1}$ in $B_{r_{i}}$. For $k\in[m,M]$ and $\rho\in(0,R]$, write $v_{i}=(u-k_{i})^{+}$ and $A(k,\rho)=\{x\in B_{\rho}: u \gt k\}$. By Theorem 1.1 in [Reference Li and Yan30], we have the following anisotropic Caffarelli–Kohn–Nirenberg type inequality:

\begin{align*} \|u\|_{L^{\frac{(n+\theta_{1}+\theta_{2})p}{n+\theta_{1}+\theta_{2}-p}}(B_{R},w)}\leq C(n,p,\theta_{1},\theta_{2})\|\nabla u\|_{L^{p}(B_{R},w)},\quad \forall u\in W_{0}^{1,p}(B_{R},w), \end{align*}

which, together with Lemma 3.1, reads that

\begin{align*} \left(\int_{B_{R}}|\eta_{i}v_{i}|^{p\chi}w\,{\rm d}x\right)^{\frac{1}{\chi}}\leq C\int_{B_{R}}|v_{i}|^{p}|\nabla\eta_{i}|^{p}w\,{\rm d}x,\quad\chi=\frac{n+\theta_{1}+\theta_{2}}{n+\theta_{1}+\theta_{2}-p}. \end{align*}

Since

\begin{align*} \int_{B_{R}}|v_{i}|^{p}|\nabla\eta_{i}|^{p}w\,{\rm d}x\leq\frac{C(M-k_{i})^{p}}{(r_{i}-r_{i+1})^{p}}|A(k_{i},r_{i})|_{\mu} \end{align*}

and

\begin{align*} \left(\int_{B_{R}}|\eta_{i}v_{i}|^{p\chi}w\,{\rm d}x\right)^{\frac{1}{\chi}}\geq(k_{i+1}-k_{i})^{p}|A(k_{i+1},r_{i+1})|_{\mu}^{\frac{1}{\chi}}, \end{align*}

it then follows that there exists a positive constant $i_{0}=i_{0}(n,p,\theta_{1},\theta_{2}) \gt 0$ such that for $i\geq i_{0}$,

\begin{align*} \frac{|A(k_{i+1},r_{i+1})|_{\mu}}{|B_{R}|_{\mu}}\leq&(C2^{2p(i+2)})^{\chi}R^{(\theta_{1}+\theta_{2}+n)(\chi-1)-p\chi}\left(\frac{|A(k_{i},r_{i})|_{\mu}}{|B_{R}|_{\mu}}\right)^{\chi}\notag\\ =&(C2^{2p(i+2)})^{\chi}\left(\frac{|A(k_{i},r_{i})|_{\mu}}{|B_{R}|_{\mu}}\right)^{\chi}\notag\\ \leq&\prod\limits^{i}_{s=0}(C2^{2p(i+2-s)})^{\chi^{s+1}}\left(\frac{|A(k_{0},r_{0})|_{\mu}}{|B_{R}|_{\mu}}\right)^{\chi^{i+1}}\notag\\ \leq&(C^{\ast})^{i+1}\left(\frac{|A(k_{0},r_{0})|_{\mu}}{|Q_{R}|_{\mu}}\right)^{\chi(i+1)}. \end{align*}

Fix $\tau_{0}=(C^{\ast})^{-\chi}$. Then we deduce that for any ɛ > 0 and $0 \lt \tau \lt \tau_{0}$, if Equation (3.2) holds, then

\begin{align*} \frac{|A(k_{i+1},r_{i+1})|_{\mu}}{|B_{R}|_{\mu}}\leq (C^{\ast}\tau^{\chi})^{i+1}=\left(\frac{\tau}{\tau_{0}}\right)^{\chi(i+1)}\rightarrow0,\quad\mathrm{as}\;i\rightarrow\infty. \end{align*}

Hence, Equation (3.3) is proved.

Step 2. Similarly as above, set

\begin{align*} r_{i}=\frac{R}{2}+\frac{R}{2^{i+1}},\quad \tilde{k}_{i}=m+\varepsilon-\frac{\varepsilon}{2}(1-2^{-i}),\quad i\geq0. \end{align*}

For $k\in[m,M]$ and $\rho\in(0,R]$, let $\tilde{v}_{i}=(u-\tilde{k}_{i})^{-}$ and $\tilde{A}(k,\rho)=\{x\in B_{\rho}: u \lt k\}$. Then applying the proof of Equation (3.3) with minor modification, we obtain that Equation (3.4) also holds. The proof is finished.

The decaying estimates for the distribution function of the solution u are established as follows.

Lemma 3.4. Suppose that $n\geq2$, $(\theta_{1},\theta_{2})\in[(\mathcal{A}\cup\mathcal{B})\cap(C_{q}\cup\mathcal{D}_{q})\cap\mathcal{F}]\cup\{\theta_{1}=0,\theta_{2} \gt -(n-1)\}$, $1 \lt q \lt p \lt n+\theta_{1}+\theta_{2}$, $0 \lt \gamma \lt 1$, $0 \lt R \lt \frac{1}{2}$ and $0\leq m\leq\inf\limits_{B_{2R}}u\leq\sup\limits_{B_{2R}}u\leq M$. Then for any ɛ > 0,

  1. (a) if

    (3.5)\begin{align} \frac{|\{x\in B_{R}:u \gt M-\varepsilon\}|_{\mu}}{|B_{R}|_{\mu}}\leq1-\gamma, \end{align}
    then for any $j\geq1$,
    (3.6)\begin{align} \frac{|\{x\in B_{R}:u \gt M-\frac{\varepsilon}{2^{j}}\}|_{\mu}}{|B_{R}|_{\mu}}\leq\frac{C}{\sqrt[q]{\gamma}j^{\frac{p-q}{pq}}}; \end{align}
  2. (b) if

    (3.7)\begin{align} \frac{|\{x\in B_{R}:u \lt m+\varepsilon\}|_{\mu}}{|B_{R}|_{\mu}}\leq1-\gamma, \end{align}
    then for any $j\geq1$,
    \begin{align*} \frac{|\{x\in B_{R}:u \lt m+\frac{\varepsilon}{2^{j}}\}|_{\mu}}{|B_{R}|_{\mu}}\leq\frac{C}{\sqrt[q]{\gamma}j^{\frac{p-q}{pq}}}, \end{align*}
    where $C=C(n,p,q,\theta_{1},\theta_{2},\lambda)$.

Proof. Step 1. For $i\geq0$, let $k_{i}=M-\frac{\varepsilon}{2^{i}}$ and $A(k_{i},R)=B_{R}\cap\{u \gt k_{i}\}.$ From Equation (2.13), we know that for $q \gt 1,$

(3.8)\begin{align} &(k_{i+1}-k_{i})^{q}|A(k_{i+1},R)|_{\mu}^{q}|B_{R}\setminus A(k_{i},R)|_{\mu}\notag\\ &\leq CR^{q(n+\theta_{1}+\theta_{2}+1)}\int_{A(k_{i},R)\setminus A(k_{i+1},R)}|\nabla u|^{q}w\,{\rm d}x. \end{align}

Using Equation (3.5), we have

\begin{align*} |B_{R}\setminus A(k_{i},R)|_{\mu}\geq\gamma|B_{R}|_{\mu}=C(n,\theta_{1},\theta_{2})\gamma R^{n+\theta_{1}+\theta_{2}}. \end{align*}

This, together with Equation (3.8), shows that

\begin{align*} &|A(k_{i+1},R)|_{\mu}\leq\frac{C2^{i+1}}{\varepsilon\sqrt[q]{\gamma}}R^{\frac{(n+\theta_{1}+\theta_{2})(q-1)}{q}+1}\left(\int_{A(k_{i},R)\setminus A(k_{i+1},R)}|\nabla u|^{q}w\,{\rm d}x\right)^{\frac{1}{q}}. \end{align*}

Since $1 \lt q \lt p \lt n+\theta_{1}+\theta_{2}$, we then have from Hölder’s inequality that

\begin{align*} &\int_{A(k_{i},R)\setminus A(k_{i+1},R)}|\nabla u|^{q}w\,{\rm d}x\notag\\ &\leq\left(\int_{A(k_{i},R)\setminus A(k_{i+1},R)}|\nabla u|^{p}w\,{\rm d}x\right)^{\frac{q}{p}}\left(\int_{A(k_{i},R)\setminus A(k_{i+1},R)}w\,{\rm d}x\right)^{\frac{p-q}{p}}\notag\\ &\leq\left(\int_{B_{R}}|\nabla(u-k_{i})^{+}|^{p}w\,{\rm d}x\right)^{\frac{q}{p}}\left(\int_{A(k_{i},R)\setminus A(k_{i+1},R)}w\,{\rm d}x\right)^{\frac{p-q}{p}}. \end{align*}

Choose a cutoff function $\eta\in C_{0}^{\infty}(B_{2R})$ satisfying that

(3.9)\begin{align} \eta=1\;\mathrm{in}\;B_{R},\quad0\leq\eta\leq1,\;|\nabla\eta|\leq\frac{C(n)}{R}\;\,\mathrm{in}\;B_{2R}. \end{align}

It then follows from Lemma 3.1 that

\begin{align*} &\left(\int_{B_{R}}|\nabla (u-k_{i})^{+}|^{p}w\,{\rm d}x\right)^{\frac{1}{p}}\leq C\left(\int_{B_{2R}}|(u-k_{i})^{+}|^{p}|\nabla\eta|^{p}w\,{\rm d}x\right)^{\frac{1}{p}}\leq\frac{C\varepsilon}{2^{i}}R^{\frac{n+\theta_{1}+\theta_{2}-p}{p}}. \end{align*}

A combination of these above facts shows that

\begin{align*} |A(k_{i+1},R)|_{\mu}\leq&\frac{C}{\sqrt[q]{\gamma}}R^{(n+\theta_{1}+\theta_{2})(1-\frac{p-q}{pq})}|A(k_{i},R)\setminus A(k_{i+1},R)|_{\mu}^{\frac{p-q}{pq}}. \end{align*}

This leads to that for $j\geq1$,

\begin{align*} j|A(k_{j},R)|_{\nu}^{\frac{pq}{p-q}}\leq&\sum^{j-1}_{i=0}|A(k_{i+1},R)|_{\mu}^{\frac{pq}{p-q}}\leq\frac{C}{\gamma^{\frac{p}{p-q}}}R^{(n+\theta_{1}+\theta_{2})\left(\frac{pq}{p-q}-1\right)}|B_{R}|_{\mu}\notag\\ \leq&\frac{C}{\gamma^{\frac{p}{p-q}}}|B_{R}|_{\mu}^{\frac{pq}{p-q}}. \end{align*}

Then Equation (3.6) is proved.

Step 2. For $i\geq0$, denote $\tilde{k}_{i}=m+\frac{\varepsilon}{2^{i}}$ and $\tilde{A}(k_{i},R)=B_{R}\cap\{u \lt k_{i}\}.$ In light of Equation (2.14), we see that for $q \gt 1,$

\begin{align*} &(\tilde{k}_{i}-\tilde{k}_{i+1})^{q}|\tilde{A}(\tilde{k}_{i+1},R)|_{\mu}^{q}|B_{R}\setminus \tilde{A}(\tilde{k}_{i},R)|_{\mu}\notag\\ &\leq CR^{q(n+\theta_{1}+\theta_{2}+1)}\int_{\tilde{A}(\tilde{k}_{i},R)\setminus \tilde{A}(\tilde{k}_{i+1},R)}|\nabla u|^{q}w\,{\rm d}x. \end{align*}

From Equation (3.7), we have

\begin{align*} |B_{R}\setminus \tilde{A}(\tilde{k}_{i},R)|_{\mu}\geq\gamma|B_{R}|_{\mu}=C(n,\theta_{1},\theta_{2})\gamma R^{n+\theta_{1}+\theta_{2}}. \end{align*}

Hence, we obtain

\begin{align*} &|\tilde{A}(\tilde{k}_{i+1},R)|_{\mu}\leq\frac{C2^{i+1}}{\varepsilon\sqrt[q]{\gamma}}R^{\frac{(n+\theta_{1}+\theta_{2})(q-1)}{q}+1}\left(\int_{\tilde{A}(\tilde{k}_{i},R)\setminus \tilde{A}(\tilde{k}_{i+1},R)}|\nabla u|^{q}w\,{\rm d}x\right)^{\frac{1}{q}}. \end{align*}

Analogously as before, it follows from Hölder’s inequality and Lemma 3.1 that for $1 \lt q \lt p \lt n+\theta_{1}+\theta_{2}$

\begin{align*} &\left(\int_{\tilde{A}(\tilde{k}_{i},R)\setminus \tilde{A}(\tilde{k}_{i+1},R)}|\nabla u|^{q}w\,{\rm d}x\right)^{\frac{1}{q}}\notag\\ &\leq\left(\int_{\tilde{A}(\tilde{k}_{i},R)\setminus \tilde{A}(\tilde{k}_{i+1},R)}|\nabla u|^{p}w\,{\rm d}x\right)^{\frac{1}{p}}\left(\int_{\tilde{A}(\tilde{k}_{i},R)\setminus \tilde{A}(\tilde{k}_{i+1},R)}w\,{\rm d}x\right)^{\frac{p-q}{pq}}\notag\\ &\leq\left(\int_{B_{R}}|\nabla(u-\tilde{k}_{i})^{-}|^{p}w\,{\rm d}x\right)^{\frac{1}{p}}\left(\int_{\tilde{A}(\tilde{k}_{i},R)\setminus \tilde{A}(\tilde{k}_{i+1},R)}w\,{\rm d}x\right)^{\frac{p-q}{pq}}\notag\\ &\leq C\left(\int_{B_{2R}}|(u-\tilde{k}_{i})^{-}|^{p}|\nabla\eta|^{p}w\,{\rm d}x\right)^{\frac{1}{p}}|\tilde{A}(\tilde{k}_{i},R)\setminus \tilde{A}(\tilde{k}_{i+1},R)|_{\mu}^{\frac{p-q}{pq}}\notag\\ &\leq\frac{C\varepsilon}{2^{i}}R^{\frac{n+\theta_{1}+\theta_{2}-p}{p}}|\tilde{A}(\tilde{k}_{i},R)\setminus \tilde{A}(\tilde{k}_{i+1},R)|_{\mu}^{\frac{p-q}{pq}}, \end{align*}

where η is given by Equation (3.9). Then we obtain

\begin{align*} |\tilde{A}(\tilde{k}_{i+1},R)|_{\mu}\leq&\frac{C}{\sqrt[q]{\gamma}}R^{(n+\theta_{1}+\theta_{2})(1-\frac{p-q}{pq})}|\tilde{A}(\tilde{k}_{i},R)\setminus \tilde{A}(\tilde{k}_{i+1},R)|_{\mu}^{\frac{p-q}{pq}}, \end{align*}

and thus,

\begin{align*} j|\tilde{A}(\tilde{k}_{j},R)|_{\nu}^{\frac{pq}{p-q}}\leq&\sum^{j-1}_{i=0}|\tilde{A}(\tilde{k}_{i+1},R)|_{\mu}^{\frac{pq}{p-q}}\leq\frac{C}{\gamma^{\frac{p}{p-q}}}R^{(n+\theta_{1}+\theta_{2})\left(\frac{pq}{p-q}-1\right)}|B_{R}|_{\mu}\notag\\ \leq&\frac{C}{\gamma^{\frac{p}{p-q}}}|B_{R}|_{\mu}^{\frac{pq}{p-q}}\quad \text{for }j\geq1. \end{align*}

The proof is complete.

A combination of Lemmas 3.2 and 3.4 yields the following improvement on oscillation of u in a small domain.

Corollary 3.5. Assume that $n\geq2$, $(\theta_{1},\theta_{2})\in[(\mathcal{A}\cup\mathcal{B})\cap(C_{q}\cup\mathcal{D}_{q})\cap\mathcal{F}]\cup\{\theta_{1}=0,\,\theta_{2} \gt -(n-1)\}$, $1 \lt q \lt p \lt n+\theta_{1}+\theta_{2}$, $0 \lt \gamma \lt 1$, $0 \lt R \lt \frac{1}{2}$ and $0\leq m\leq\inf\limits_{B_{2R}}u\leq\sup\limits_{B_{2R}}u\leq M$. Then there exists a large constant $k_{0} \gt 1$ depending only on $n,p,q,\theta_{1},\theta_{2},\lambda,\gamma$ such that for any ɛ > 0,

  1. (i) if

    \begin{align*} \frac{|\{x\in B_{R}: u \gt M-\varepsilon\}|_{\mu}}{|B_{R}|_{\mu}}\leq1-\gamma, \end{align*}
    then
    \begin{align*} \sup\limits_{B_{R/2}}u\leq M-\frac{\varepsilon}{2^{k_{0}}}; \end{align*}
  2. (ii) if

    \begin{align*} \frac{|\{x\in B_{R}: u \lt m+\varepsilon\}|_{\mu}}{|B_{R}|_{\mu}}\leq1-\gamma, \end{align*}
    then
    \begin{align*} \inf\limits_{B_{R/2}}u\geq m+\frac{\varepsilon}{2^{k_{0}}}. \end{align*}

Proof. Applying Lemmas 3.2 and 3.4, we obtain that Corollary 3.5 holds. In particular, in the case of $\theta_{1}=0$, $\theta_{2} \gt -(n-1)$, $1 \lt p \lt n+\theta_{2}$, we fix $q=\frac{p+1}{2}$ in Lemma 3.4.

We are now ready to prove Theorems 1.1 and 1.4, respectively.

Proof of Theorem 1.1

For $0 \lt R\leq\frac{1}{2}$, denote

\begin{align*} \overline{\mu}(R)=\sup_{B_{R}}u,\quad\underline{\mu}(R)=\inf_{B_{R}}u,\quad \omega(R)=\overline{\mu}(R)-\underline{\mu}(R). \end{align*}

Note that one of the following two statements must hold: either

(3.10)\begin{align} |\{x\in B_{R}:u \gt \overline{\mu}(R)-2^{-1}\omega(R)\}|_{\mu}\leq\frac{1}{2}|B_{R}|_{\mu} \end{align}

or

(3.11)\begin{align} |\{x\in B_{R}:u \lt \underline{\mu}(R)+2^{-1}\omega(R)\}|_{\mu}\leq\frac{1}{2}|B_{R}|_{\mu}. \end{align}

Using Corollary 3.5 with $\gamma=\frac{1}{2}$, we derive that there is a large constant $k_{0} \gt 1$ such that

\begin{align*} \overline{\mu}(R/2)\leq\overline{\mu}(R)-\frac{\omega(R)}{2^{k_{0}+1}},\quad\text{when Equation~(3.10) holds} \end{align*}

and

\begin{align*} \underline{\mu}(R/2)\geq\underline{\mu}(R)+\frac{\omega(R)}{2^{k_{0}+1}},\quad\text{when Equation~(3.11) holds.} \end{align*}

In either case, we have

\begin{align*} \omega(R/2)\leq\left(1-\frac{1}{2^{k_{0}+1}}\right)\omega(R)=\frac{1}{2^{\alpha}}\omega(R),\quad\mathrm{with}\;\alpha=-\frac{\ln\big(1-\frac{1}{2^{k_{0}+1}}\big)}{\ln2}. \end{align*}

Observe that for each $0 \lt R\leq\frac{1}{2}$, there exists an integer $l\geq1$ such that $2^{-(l+1)} \lt R\leq2^{-l}$. Since $\omega(R)$ is non-decreasing with respect to R, we then have

\begin{align*} \omega(R)\leq\omega(2^{-l})\leq 2^{-(l-1)\alpha}\omega(2^{-1})=4^{\alpha}2^{-(l+1)\alpha}\omega(2^{-1})\leq CR^{\alpha}, \end{align*}

where $C=C(n,p,q,\theta_{1},\theta_{2},\lambda,\overline{M}).$ The proof is complete.

Proof of Theorem 1.4

First, by applying the proof of Theorem 1.1 with a slight modification, we also obtain that there exist a small constant $0 \lt \alpha=\alpha(n,p,\theta_{2},\lambda) \lt 1$ and a large constant $C=C(n,p,\theta_{2},\lambda,\overline{M}) \gt 0$ such that

(3.12)\begin{align} |u(x)-u(0)|\leq C|x|^{\alpha}\quad\text{for all}\;x\in B_{1/2}. \end{align}

For $R\in(0,1/2)$, $y\in Q_{1/R}$, denote

\begin{align*} u_{R}(y)=u(Ry),\qquad A_{R}(y)=A(Ry). \end{align*}

Hence, u R is the solution of

\begin{align*} \mathrm{div}(A_{R}|y|^{\theta_{2}}\nabla u_{R})=0\quad\mathrm{in}\;Q_{1/R}. \end{align*}

After the change of variables, we see that this equation becomes degenerate elliptic equation in $B_{1/2}(\bar{y})$ for any $\bar{y}\in\partial B_{1}$. For any two given points $x,\tilde{x}\in B_{1/2},$ let $|\tilde{x}|\leq|x|$ without loss of generality. Denote $R=|x|$. By the interior Hölder estimate for degenerate elliptic equation, we derive that there exist two constants $0 \lt \beta=\beta(n,p,\theta_{2},\lambda) \lt 1$ and $0 \lt C=C(n,p,\theta_{2},\lambda,\overline{M})$ such that for any $\bar{y}\in\partial B_{1}$,

(3.13)\begin{align} |u_{R}(y)-u_{R}(\bar{y})|\leq C|y-\bar{y}|^{\beta},\quad\forall y\in B_{1/2}(\bar{y}). \end{align}

Consequently, for $|x-\tilde{x}|\leq R^{2}$, we have from Equation (3.13) that

\begin{align*} |u(x)-u(\tilde{x})|=&\left|u_{R}(x/R)-u_{R}(\tilde{x}/R)\right|\leq C|(x-\tilde{x})/R|^{\beta}\leq C|x-\tilde{x}|^{\beta/2}, \end{align*}

while, for $|x-\tilde{x}| \gt R^{2}$, we deduce from Equation (3.12) that

\begin{align*} |u(x)-u(\tilde{x})|\leq&|u(x)-u(0)|+|u(0)-u(\tilde{x})|\leq C\big(R^{\alpha}+|\tilde{x}|^{\alpha}\big)\leq CR^{\alpha}\leq C|x-\tilde{x}|^{\frac{\alpha}{2}}. \end{align*}

Therefore, the proof of Theorem 1.4 is complete.

4. Regularity for solutions to nonlinear parabolic equations with anisotropic weights

Let $n\geq2$, R > 0 and $-T\leq t_{1} \lt t_{2}\leq0$. For $u\in C((t_{1},t_{2});L^{2}(B_{R},w_{1}))\cap L^{2}((t_{1},t_{2});W_{0}^{1,2}(B_{R},w_{2}))$, denote

\begin{align*} \|u\|_{V^{1}_{0}(B_{R}\times(t_{1},t_{2}))}=\sqrt{\sup\limits_{t\in(t_{1},t_{2})}\int_{B_{R}}|u|^{2}w_{1}\,{\rm d}x+\int^{t_{2}}_{t_{1}}\int_{B_{R}}|\nabla u|^{2}w_{2}\,{\rm d}x\,{\rm d}t}, \end{align*}

where the anisotropic weights w 1 and w 2 are defined in Equation (1.5). The parabolic Sobolev inequality with anisotropic weights is now given as follows.

Proposition 4.1. For $n\geq2$, R > 0, $\theta_{1}+\theta_{2} \gt -(n-2)$ and $-T\leq t_{1} \lt t_{2}\leq0$, let $u\in C((t_{1},t_{2});L^{2}(B_{R},w_{1}))\cap L^{2}((t_{1},t_{2});W_{0}^{1,2}(B_{R},w_{2}))$. Then

\begin{align*} \|u\|_{L^{2\chi}(B_{R}\times(t_{1},t_{2}),w_{2})}\leq C\|u\|_{V^{1}_{0}(B_{R}\times(t_{1},t_{2}))},\quad\mathrm{with}\;\chi=\frac{n+\theta_{1}+\theta_{2}+2}{n+\theta_{1}+\theta_{2}}, \end{align*}

where $C=C(n,\theta_{1},\theta_{2}).$

Proof. Applying the anisotropic version of the Caffarelli–Kohn–Nirenberg inequality in [Reference Li and Yan30], we obtain that for any $u\in W_{0}^{1,2}(B_{R},w_{2})$,

\begin{align*} &\left(\int_{B_{R}}|u|^{\frac{2(n+\theta_{1}+\theta_{2})}{n+\theta_{1}+\theta_{2}-2}}|x^{\prime}|^{\frac{\theta_{3}(n+\theta_{1}+\theta_{2})-2\theta_{1}}{n+\theta_{1}+\theta_{2}-2}}|x|^{\frac{\theta_{4}(n+\theta_{1}+\theta_{2})-2\theta_{2}}{n+\theta_{1}+\theta_{2}-2}}\,{\rm d}x\right)^{\frac{n+\theta_{1}+\theta_{2}-2}{n+\theta_{1}+\theta_{2}}}\notag\\ &\leq C\int_{B_{R}}|\nabla u|^{2}|x^{\prime}|^{\theta_{3}}|x|^{\theta_{4}}\,{\rm d}x. \end{align*}

This, in combination with the Hölder’s inequality, leads to

(4.1)\begin{align} &\int_{B_{R}}|u|^{2\chi}|x^{\prime}|^{\theta_{3}}|x|^{\theta_{4}}\,{\rm d}x\notag\\ &=\int_{B_{R}}|u|^{2}|x^{\prime}|^{\theta_{3}-\theta_{1}(\chi-1)}|x^{\prime}|^{\theta_{4}-\theta_{2}(\chi-1)}|u|^{2(\chi-1)}|x^{\prime}|^{\theta_{1}(\chi-1)}|x|^{\theta_{2}(\chi-1)}\,{\rm d}x\notag\\ &\leq\left(\int_{B_{R}}|u|^{\frac{2}{2-\chi}}|x^{\prime}|^{\frac{\theta_{3}-\theta_{1}(\chi-1)}{2-\chi}}|x|^{\frac{\theta_{4}-\theta_{2}(\chi-1)}{2-\chi}}\,{\rm d}x\right)^{2-\chi}\left(\int_{B_{R}}|u|^{2}|x^{\prime}|^{\theta_{1}}|x|^{\theta_{2}}\,{\rm d}x\right)^{\chi-1}\notag\\ &\leq C\int_{B_{R}}|\nabla u|^{2}|x^{\prime}|^{\theta_{3}}|x|^{\theta_{4}}\,{\rm d}x\left(\int_{B_{R}}|u|^{2}|x^{\prime}|^{\theta_{1}}|x|^{\theta_{2}}\,{\rm d}x\right)^{\chi-1}. \end{align}

Then integrating Equation (4.1) from t 1 to t 2, it follows from Young’s inequality that

\begin{align*} \left(\int^{t_{2}}_{t_{1}}\int_{B_{R}}|u|^{2\chi}w_{2}\right)^{\frac{1}{\chi}}\leq& C\left(\sup\limits_{t\in(t_{1},t_{2})}\int_{B_{R}}|u|^{2}w_{1}\,{\rm d}x\right)^{\frac{\chi-1}{\chi}}\left(\int_{B_{R}}|\nabla u|^{2}w_{2}\,{\rm d}x\right)^{\frac{1}{\chi}}\notag\\ \leq&C\left(\sup\limits_{t\in(t_{1},t_{2})}\int_{B_{R}}|u|^{2}w_{1}\,{\rm d}x+\int^{t_{2}}_{t_{1}}\int_{B_{R}}|\nabla u|^{2}w_{2}\,{\rm d}x\,{\rm d}t\right). \end{align*}

The proof is complete.

The Caccioppoli inequality for the truncated solution is given as follows.

Lemma 4.2. Set $\overline{m}\leq k\leq\overline{M}$. Then for any $B_{R}(x_{0})\subset B_{1}$ and non-negative $\eta\in C^{\infty}(B_{R}(x_{0})\times(-1,0))$, which vanishes on $\partial B_{R}(x_{0})\times(-1,0)$, we obtain that for $-1 \lt t_{1} \lt t_{2} \lt 0$,

\begin{align*} &\sup\limits_{t\in(t_{1},t_{2})}\int_{B_{R}(x_{0})}(v\eta)^{2}w_{1}\,{\rm d}x+\int_{B_{R}(x_{0})\times(t_{1},t_{2})}|\nabla(v\eta)|^{2}w_{2}\,{\rm d}x\,{\rm d}t\notag\\ &\leq\int_{B_{R}(x_{0})}(v^{2}+C_{0}v^{3})\eta^{2}w_{1}\,{\rm d}x\Big|_{t_{1}}+C_{0}\int_{B_{R}(x_{0})\times(t_{1},t_{2})}(\eta|\partial_{t}\eta|w_{1}+|\nabla\eta|^{2}w_{2})v^{2}\,{\rm d}x\,{\rm d}t \end{align*}

and

\begin{align*} &\sup\limits_{t\in(t_{1},t_{2})}\int_{B_{R}(x_{0})}(\tilde{v}^{2}-C_{0}\tilde{v}^{3})\eta^{2}w_{1}\,{\rm d}x+\int_{B_{R}(x_{0})}|\nabla(\tilde{v}\eta)|^{2}w_{2}\,{\rm d}x\,{\rm d}t\notag\\ &\leq\int_{B_{R}(x_{0})}(\tilde{v}\eta)^{2}w_{1}\,{\rm d}x\Big|_{t_{1}}+C_{0}\int_{B_{R}(x_{0})\times(t_{1},t_{2})}\left(\eta|\partial_{t}\eta|w_{1}+|\nabla\eta|^{2}w_{2}\right)\tilde{v}^{2}\,{\rm d}x\,{\rm d}t, \end{align*}

where $C_{0}=C_{0}(n,p,\lambda,\overline{m},\overline{M})$, $v=(u-k)^{+}$, $\tilde{v}=(u-k)^{-}$, u is the solution of problem (1.5).

Proof. Choose test function $\varphi=v\eta^{2}$. By denseness, we obtain that for $t_{1}\leq s\leq t_{2}$,

\begin{align*} \int^{s}_{t_{1}}\int_{B_{R}(x_{0})}w_{1}\partial_{t}u^{p}v\eta^{2}\,{\rm d}x\,{\rm d}t+\int^{s}_{t_{1}}\int_{B_{R}(x_{0})}Aw_{2}\nabla u\nabla(v\eta^{2})\,{\rm d}x\,{\rm d}t=0. \end{align*}

On the one hand,

\begin{align*} \int^{s}_{t_{1}}\int_{B_{R}(x_{0})}w_{1}\partial_{t}u^{p}v\eta^{2}dxdt=&p \int^{s}_{t_{1}}\int_{B_{R}(x_{0})}w_{1}v\eta^{2}u^{p-1}\partial_{t}v\,dxdt\notag\\ =&p \int^{s}_{t_{1}}\int_{B_{R}(x_{0})}w_{1}\eta^{2}[(v+k)^{p}-k(v+k)^{p-1}]\partial_{t}v\,dxdt\notag\\ =&p \int^{s}_{t_{1}}\int_{B_{R}(x_{0})}w_{1}\eta^{2}\partial_{t}\mathcal{H}, \end{align*}

where

\begin{align*} \mathcal{H}:=\frac{(v+k)^{p+1}}{p+1}-\frac{k(v+k)^{p}}{p}+\frac{k^{p+1}}{p(p+1)}. \end{align*}

Remark that the last term $\frac{k^{p+1}}{p(p+1)}$ in $\mathcal{H}$ is added to keep it non-negative. In fact, by Taylor expansion, we obtain

\begin{align*} \frac{(v+k)^{p+1}}{p+1}=&\frac{k^{p+1}}{p+1}\left(1+\frac{v}{k}\right)^{p+1}\notag\\ =&\frac{k^{p+1}}{p+1}\left(1+(p+1)\frac{v}{k}+\frac{p(p+1)}{2}\frac{v^{2}}{k^{2}}+\frac{(p-1)p(p+1)}{6}\frac{v^{3}}{k^{3}}+O\left(\frac{v^{4}}{k^{4}}\right)\right), \end{align*}

and

\begin{align*} \frac{k(v+k)^{p}}{p}=&\frac{k^{p+1}}{p}\left(1+\frac{v}{k}\right)^{p}\notag\\ =&\frac{k^{p+1}}{p}\left(1+p\frac{v}{k}+\frac{p(p-1)}{2}\frac{v^{2}}{k^{2}}+\frac{p(p-1)(p-2)}{6}\frac{v^{3}}{k^{3}}+O\left(\frac{v^{4}}{k^{4}}\right)\right). \end{align*}

A consequence of these two relations shows that

\begin{align*} \frac{(v+k)^{p+1}}{p+1}-\frac{k(v+k)^{p}}{p}=-\frac{k^{p+1}}{p(p+1)}+\frac{1}{2}k^{p-1}v^{2}+\frac{p-1}{3}k^{p-2}v^{3}+O\left(\frac{v^{4}}{k^{4}}\right), \end{align*}

which yields that

(4.2)\begin{align} 0\leq\mathcal{H}-\frac{1}{2}k^{p-1}v^{2}\leq C(p,\overline{m},\overline{M})v^{3}. \end{align}

In light of Equation (4.2), it follows from integration by parts that

\begin{align*} &\int^{s}_{t_{1}}\int_{B_{R}(x_{0})}w_{1}\partial_{t}u^{p}v\eta^{2}\,{\rm d}x\,{\rm d}t\notag\\ &\geq\frac{p}{2}k^{p-1}\int_{B_{R}(x_{0})}w_{1}\eta^{2}(x,s)v^{2}(x,s)\,{\rm d}x-\frac{p}{2}k^{p-1}\int_{B_{R}(x_{0})}w_{1}\eta^{2}(x,t_{1})v^{2}(x,t_{1})\,{\rm d}x\notag\\ &\quad-C\int_{B_{R}(x_{0})}w_{1}\eta^{2}(x,t_{1})v^{3}(x,t_{1})\,{\rm d}x-C\int^{s}_{t_{1}}\int_{B_{R}(x_{0})}w_{1}\eta|\partial_{t}\eta|v^{2}\,{\rm d}x\,{\rm d}t. \end{align*}

On the other hand, utilizing Young’s inequality, we have

\begin{align*} &\int_{B_{R}(x_{0})\times(t_{1},s)}Aw_{2}\nabla u\nabla(v\eta^{2})\,{\rm d}x\,{\rm d}t\notag\\ &=\int_{B_{R}(x_{0})\times(t_{1},s)}Aw_{2}\nabla v(\eta^{2}\nabla v+2v\eta\nabla\eta)\,{\rm d}x\,{\rm d}t\notag\\ &\geq\frac{1}{\lambda}\int_{B_{R}(x_{0})\times(t_{1},s)}|\eta\nabla v|^{2}w_{2}\,{\rm d}x\,{\rm d}t-C\int_{B_{R}(x_{0})\times(t_{1},s)}v\eta|\nabla\eta||\nabla v|w_{2}\,{\rm d}x\,{\rm d}t\notag\\ &\geq\frac{1}{2\lambda}\int_{B_{R}(x_{0})\times(t_{1},s)}|\eta\nabla v|^{2}w_{2}\,{\rm d}x\,{\rm d}t-C\int_{B_{R}(x_{0})\times(t_{1},s)}v^{2}|\nabla\eta|^{2}w_{2}\,{\rm d}x\,{\rm d}t\notag\\ &\geq\frac{1}{4\lambda}\int_{B_{R}(x_{0})\times(t_{1},s)}|\nabla(\eta v)|^{2}w_{2}\,{\rm d}x\,{\rm d}t-C\int_{B_{R}(x_{0})\times(t_{1},s)}v^{2}|\nabla\eta|^{2}w_{2}\,{\rm d}x\,{\rm d}t, \end{align*}

where in the last inequality, we used the following elementary inequality:

\begin{align*} |a-b|^{2}\geq\frac{1}{2}|a|^{2}-|b|^{2}\quad\text{for any } a,b\in\mathbb{R}^{n}. \end{align*}

Therefore, the first inequality in Lemma 4.2 holds.

The proof of the second inequality in Lemma 4.2 is analogous by picking test function $\varphi=-\tilde{v}\eta^{2}$. Then we obtain

\begin{align*} -\int^{s}_{t_{1}}\int_{B_{R}(x_{0})}w_{1}\partial_{t}u^{p}\tilde{v}\eta^{2}=&p \int^{s}_{t_{1}}\int_{B_{R}(x_{0})}w_{1}\tilde{v}\eta^{2}u^{p-1}\,\partial_{t}\tilde{v}\notag\\ =&p \int^{s}_{t_{1}}\int_{B_{R}(x_{0})}w_{1}\eta^{2}[-(k-\tilde{v})^{p}+k(k-\tilde{v})^{p-1}]\,\partial_{t}\tilde{v}\notag\\ =&p \int^{s}_{t_{1}}\int_{B_{R}(x_{0})}w_{1}\eta^{2}\,\partial_{t}\widetilde{\mathcal{H}}, \end{align*}

where

\begin{align*} \widetilde{\mathcal{H}}=\frac{(k-\tilde{v})^{p+1}}{p+1}-\frac{k(k-\tilde{v})^{p}}{p}+\frac{k^{p+1}}{p(p+1)}. \end{align*}

Similarly as before, it follows from Taylor expansion that

\begin{align*} -C(p,\overline{m},\overline{M})\tilde{v}^{3}\leq\widetilde{\mathcal{H}}-\frac{1}{2}k^{p-1}\tilde{v}^{2}\leq0, \end{align*}

which reads that

\begin{align*} &\int^{s}_{t_{1}}\int_{B_{R}(x_{0})}w_{1}\partial_{t}u^{p}v\eta^{2}\,{\rm d}x\,{\rm d}t\notag\\ &\geq\frac{p}{2}k^{p-1}\int_{B_{R}(x_{0})}w_{1}\eta^{2}(x,s)\tilde{v}^{2}(x,s)\,{\rm d}x-C\int_{B_{R}(x_{0})}w_{1}\eta^{2}(x,s)\tilde{v}^{3}(x,s)\,{\rm d}x\notag\\ &\quad-\frac{p}{2}k^{q-1}\int_{B_{R}(x_{0})}w_{1}\eta^{2}(x,t_{1})\tilde{v}^{2}(x,t_{1})\,{\rm d}x-C\int^{s}_{t_{1}}\int_{B_{R}(x_{0})}w_{1}\eta|\partial_{t}\eta|\tilde{v}^{2}\,{\rm d}x\,{\rm d}t. \end{align*}

By the same argument as before, we have

\begin{align*} &-\int_{B_{R}(x_{0})\times(t_{1},s)}Aw_{2}\nabla u\nabla(\tilde{v}\eta^{2})\,{\rm d}x\,{\rm d}t=\int_{B_{R}(x_{0})\times(t_{1},s)}Aw_{2}\nabla\tilde{v}\nabla(\tilde{v}\eta^{2})\,{\rm d}x\,{\rm d}t\notag\\ &\geq\frac{1}{4\lambda}\int_{B_{R}(x_{0})\times(t_{1},s)}|\nabla(\eta\tilde{v})|^{2}w_{2}\,{\rm d}x\,{\rm d}t-C\int_{B_{R}(x_{0})\times(t_{1},s)}\tilde{v}^{2}|\nabla\eta|^{2}w_{2}\,{\rm d}x\,{\rm d}t. \end{align*}

The proof is complete.

For R > 0 and $(x_{0},t_{0})\in B_{1}\times[-1+R^{\theta_{1}+\theta_{2}},0]$, denote

\begin{align*} Q_{R}(x_{0},t_{0}):=B_{R}(x_{0})\times(t_{0}-R^{\theta_{1}+\theta_{2}},t_{0}]. \end{align*}

For brevity, we use Q R to represent $Q_{R}(0,0)$ in the following. Introduce two Radon measures associated with the weights w 1 and w 2 as follows:

\begin{align*} {\rm d}\mu_{w_{i}}=w_{i}\,{\rm d}x,\qquad {\rm d}\nu_{w_{i}}=w_{i}\,{\rm d}x\,{\rm d}t,\quad i=1,2, \end{align*}

satisfying that for $E\subset B_{1}$ and $\widetilde{E}\subset Q_{1}$,

\begin{align*} |E|_{\mu_{w_{i}}}=\int_{E}w_{i}\,{\rm d}x,\qquad |\widetilde{E}|_{\nu_{w_{i}}}=\int_{\widetilde{E}}w_{i}\,{\rm d}x\,{\rm d}t. \end{align*}

Observe that by Hölder’s inequality, we know that for $\widetilde{E}\subset Q_{R}$,

(4.3)\begin{align} \frac{|\widetilde{E}|_{\nu_{w_{2}}}}{|Q_{R}|_{\nu_{w_{2}}}}\leq&\leq\frac{C|\widetilde{E}|_{\nu_{w_{1}}}^{\frac{\theta_{3}}{\theta_{1}}}R^{\frac{(n+\theta_{1}+\theta_{2})(\theta_{1}-\theta_{3})}{\theta_{1}}}}{R^{n+\theta_{1}+\theta_{2}+\theta_{3}+\theta_{4}}}\leq C\left(\frac{|\widetilde{E}|_{\nu_{w_{1}}}}{|Q_{R}|_{\nu_{w_{1}}}}\right)^{\frac{\theta_{3}}{\theta_{1}}}, \end{align}

where $C=C(n,\theta_{1},\theta_{2},\theta_{3},\theta_{4}).$ Here we used the assumed condition that $\theta_{1}/\theta_{3}=\theta_{2}/\theta_{4}$, $\theta_{3},\theta_{4}\neq0.$

Similar to Lemma 3.2, we improve the oscillation of the solution u in a small region as follows.

Lemma 4.3. Assume as in Theorem 1.6 or Theorem 1.9. For $R\in(0,1)$ and $t_{0}\in[-1+R^{\theta_{1}+\theta_{2}},0]$, let $0 \lt \overline{m}\leq m\leq\inf\limits_{Q_{R}(0,t_{0})}u\leq\sup\limits_{Q_{R}(0,t_{0})}u\leq M\leq\overline{M}$. Then

  1. (a) there exists a small constant $0 \lt \tau_{0}=\tau_{0}(n,p,\theta_{1},\theta_{2},\theta_{3},\lambda,\overline{m},\overline{M}) \lt 1$ such that for any ɛ > 0 and $0 \lt \tau \lt \tau_{0}$, if

    (4.4)\begin{align} |\{(x,t)\in Q_{R}(0,t_{0}): u(x,t) \gt M-\varepsilon\}|_{\nu_{w_{1}}}\leq\tau|Q_{R}(0,t_{0})|_{\nu_{w_{1}}}, \end{align}
    then we have
    (4.5)\begin{align} u(x,t)\leq M-\frac{\varepsilon}{2}\quad\text{for}\quad(x,t)\in Q_{R/2}(0,t_{0}); \end{align}
  2. (b) there exist two small constant $0 \lt \varepsilon_{0}=\varepsilon_{0}(n,p,\lambda,\overline{m},\overline{M}) \lt 1$ and $0 \lt \tau_{0}=\tau_{0}(n,p,\theta_{1},\theta_{2},\theta_{3},\lambda,\overline{m},\overline{M}) \lt 1$ such that for any $0 \lt \varepsilon\leq\varepsilon_{0}$ and $0 \lt \tau \lt \tau_{0}$, if

    \begin{align*} |\{(x,t)\in Q_{R}(0,t_{0}): u(x,t) \lt m+\varepsilon\}|_{\nu_{w_{1}}}\leq\tau|Q_{R}(0,t_{0})|_{\nu_{w_{1}}}, \end{align*}
    then we have
    (4.6)\begin{align} u(x,t)\geq m-\frac{\varepsilon}{2}\quad\text{for}\quad(x,t)\in Q_{R/2}(0,t_{0}). \end{align}

Remark 4.4. From the proof of Lemma 4.3 below, we see that the value of $\theta_{3}+\theta_{4}$ has to be restricted to 2 and thus affect our final regular index in Theorems 1.6 and 1.9.

Proof. Without loss of generality, let $t_{0}=0$.

Step 1. For ɛ > 0 and $i=0,1,2,\ldots,$ set

\begin{align*} r_{i}=\frac{R}{2}+\frac{R}{2^{i+1}},\qquad k_{i}=M-\varepsilon+\frac{\varepsilon}{2}(1-2^{-i}). \end{align*}

Choose a cutoff function $\eta_{i}\in C_{0}^{\infty}(Q_{r_{i}})$ such that

\begin{align*} \eta_{i}=1\;\mathrm{in}\;Q_{r_{i+1}},\quad 0\leq\eta_{i}\leq1,\;|\nabla\eta_{i}|\leq\frac{C}{r_{i}-r_{i+1}},\;|\partial_{t}\eta_{i}|\leq\frac{C}{r_{i}^{\theta_{1}+\theta_{2}}-r_{i+1}^{\theta_{1}+\theta_{2}}}\;\,\mathrm{in}\;Q_{r_{i}}. \end{align*}

Denote $v_{i}=(u-k_{i})^{+}$ and $A(k,\rho)=\{(x,t)\in Q_{\rho}: u \gt k\}$ for $k\in[m,M]$ and $\rho\in(0,R]$. Then combining Proposition 4.1 and Lemma 4.2, we deduce

(4.7)\begin{align} \left(\int_{Q_{R}}|\eta_{i}v_{i}|^{2\chi}w_{2}\right)^{\frac{1}{\chi}}\leq C\int_{Q_{R}}\left(|\nabla\eta_{i}|^{2}+|\partial_{t}\eta_{i}||x^{\prime}|^{\theta_{1}-\theta_{3}}|x|^{\theta_{2}-\theta_{4}}\right)v_{i}^{2}w_{2}, \end{align}

where $\chi=(n+\theta_{1}+\theta_{2}+2)(n+\theta_{1}+\theta_{2})^{-1}$. Note that

(4.8)\begin{align} |\partial_{t}\eta_{i}|R^{\theta_{1}+\theta_{2}-\theta_{3}-\theta_{4}}\leq\frac{CR^{2-\theta_{3}-\theta_{4}}}{(r_{i}-r_{i+1})^{2}}=\frac{C}{(r_{i}-r_{i+1})^{2}}. \end{align}

Therefore, we have

\begin{align*} \int_{Q_{R}}\left(|\nabla\eta_{i}|^{2}+|\partial_{t}\eta_{i}||x^{\prime}|^{\theta_{1}-\theta_{3}}|x|^{\theta_{2}-\theta_{4}}\right)v_{i}^{2}w_{2}\leq\frac{C(M-k_{i})^{2}}{(r_{i}-r_{i+1})^{2}}|A(k_{i},r_{i})|_{\nu_{w_{2}}} \end{align*}

and

\begin{align*} \left(\int_{Q_{R}}|\eta_{i}v_{i}|^{2\chi}w_{2}\,{\rm d}x\,{\rm d}t\right)^{\frac{1}{\chi}}\geq(k_{i+1}-k_{i})^{2}|A(k_{i+1},r_{i+1})|_{\nu_{w_{2}}}^{\frac{1}{\chi}}. \end{align*}

Define

\begin{align*} F_{i}:=\frac{|A(k_{i},r_{i})|_{\nu_{w_{2}}}}{|Q_{R}|_{\nu_{w_{2}}}}. \end{align*}

Then we have

(4.9)\begin{align} F_{i+1}\leq&(C2^{4(i+2)})^{\chi}R^{(n+\theta_{1}+\theta_{2}+\theta_{3}+\theta_{4})(\chi-1)-2\chi}F_{i}^{\chi}\notag\\ =&(C2^{4(i+2)})^{\chi}R^{\frac{2(\theta_{3}+\theta_{4}-2)}{n+\theta_{1}+\theta_{2}}}F_{i}^{\chi}=(C2^{4(i+2)})^{\chi}F_{i}^{\chi}\notag\\ \leq&\prod\limits^{i}_{s=0}(C2^{4(i+2-s)})^{\chi^{s+1}}F_{0}^{\chi^{i+1}}. \end{align}

Observe from Equations (4.8)–(4.9) that the value of $\theta_{3}+\theta_{4}$ must be chosen to be 2. A consequence of Equations (4.3) and (4.9) shows that there exists a constant $i_{0}=i_{0}(n,\theta_{1},\theta_{2},\theta_{3}) \gt 0$ such that if $i\geq i_{0}$,

\begin{align*} F_{i+1}\leq (C^{\ast})^{i+1}\left(\frac{|A(k_{0},r_{0})|_{\nu_{w_{1}}}}{|Q_{R}|_{\nu_{w_{1}}}}\right)^{\frac{\theta_{3}}{\theta_{1}}\chi^{i+1}}\leq (C^{\ast})^{i+1}\left(\frac{|A(k_{0},r_{0})|_{\nu_{w_{1}}}}{|Q_{R}|_{\nu_{w_{1}}}}\right)^{\frac{\theta_{3}}{\theta_{1}}\chi(i+1)}. \end{align*}

By taking $\tau_{0}=(C^{\ast})^{-\frac{\theta_{1}}{\theta_{3}}\chi}$, we obtain that for any ɛ > 0 and $0 \lt \tau \lt \tau_{0}$, if Equation (4.4) holds, then

\begin{align*} F_{i+1}\leq \left(C^{\ast}\tau^{\frac{\theta_{3}}{\theta_{1}}\chi}\right)^{i+1}=\left(\frac{\tau}{\tau_{0}}\right)^{\frac{\theta_{3}}{\theta_{1}}\chi(i+1)}\rightarrow0\quad\mathrm{as}\;i\rightarrow\infty. \end{align*}

That is, Equation (4.5) holds.

Step 2. Analogously as before, pick

\begin{align*} r_{i}=\frac{R}{2}+\frac{R}{2^{i+1}},\quad \tilde{k}_{i}=m+\varepsilon-\frac{\varepsilon}{2}(1-2^{-i}),\quad i\geq0. \end{align*}

Let $\varepsilon_{0}=\frac{1}{C_{0}}$, where C 0 is given in Lemma 4.2. Denote $\tilde{v}_{i}=(u-\tilde{k}_{i})^{-}$. Then we obtain that for any $0 \lt \varepsilon\leq\varepsilon_{0}$,

\begin{align*} \tilde{v}_{i}^{2}-C_{0}\tilde{v}_{i}^{3}\geq(1-C_{0}\varepsilon)\tilde{v}_{i}^{2}\geq0, \end{align*}

which implies that Equation (4.7) holds with v i replaced by $\tilde{v}_{i}$. Then following the left proof of Equation (4.5) above, we deduce that Equation (4.6) holds. The proof is complete.

The decaying estimates for the distribution function of u are stated as follows.

Lemma 4.5. Let the values of $n,p,q,\theta_{i},i=1,2,3,4$ be assumed in Theorem 1.6 or Theorem 1.9 with $\theta_{3}+\theta_{4}=2$ replaced by $0\leq\theta_{3}+\theta_{4}\leq2$. Suppose that $0 \lt \gamma \lt 1$, $0 \lt R \lt \frac{1}{2}$, $0 \lt a\leq1$, $-\frac{1}{2} \lt t_{0}\leq-aR^{\theta_{1}+\theta_{2}}$ and $\overline{m}\leq m_{a}\leq\inf\limits_{B_{2R}\times[t_{0},t_{0}+a R^{\theta_{1}+\theta_{2}}]}u\leq\sup\limits_{B_{2R}\times[t_{0},t_{0}+a R^{\theta_{1}+\theta_{2}}]}u\leq M_{a}\leq\overline{M}$. Then

  1. (a) for any ɛ > 0, if

    \begin{align*} \frac{|\{x\in B_{R}:u(x,t) \gt M_{a}-\varepsilon\}|_{\mu_{w_{1}}}}{|B_{R}|_{\mu_{w_{1}}}}\leq1-\gamma,\quad\forall t\in[t_{0},t_{0}+aR^{\theta_{1}+\theta_{2}}], \end{align*}
    then for any $j\geq1$,
    (4.10)\begin{align} \frac{|\{(x,t)\in B_{R}\times[t_{0},t_{0}+aR^{\theta_{1}+\theta_{2}}]:u(x,t) \gt M_{a}-\frac{\varepsilon}{2^{j}}\}|_{\nu_{w_{1}}}}{|B_{R}\times[t_{0},t_{0}+aR^{\theta_{1}+\theta_{2}}]|_{\nu_{w_{1}}}}\leq\frac{C}{\sqrt[q]{\gamma}\sqrt{a}j^{\frac{2-q}{2q}}}; \end{align}
  2. (b) for any $0 \lt \varepsilon\leq\varepsilon_{0}=C_{0}^{-1}$ with $C_{0}=C_{0}(n,p,\lambda,\overline{m},\overline{M})$ given by Lemma 4.2, if

    \begin{align*} \frac{|\{x\in B_{R}:u(x,t) \lt m_{a}+\varepsilon\}|_{\mu_{w_{1}}}}{|B_{R}|_{\mu_{w_{1}}}}\leq1-\gamma,\quad\forall t\in[t_{0},t_{0}+aR^{\theta_{1}+\theta_{2}}], \end{align*}
    then for any $j\geq1$,
    \begin{align*} \frac{|\{(x,t)\in B_{R}\times[t_{0},t_{0}+aR^{\theta_{1}+\theta_{2}}]:u(x,t) \lt m_{a}+\frac{\varepsilon}{2^{j}}\}|_{\nu_{w_{1}}}}{|B_{R}\times[t_{0},t_{0}+aR^{\theta_{1}+\theta_{2}}]|_{\nu_{w_{1}}}}\leq\frac{C}{\sqrt[q]{\gamma}\sqrt{a}j^{\frac{2-q}{2q}}}, \end{align*}
    where $C=C(n,p,q,\theta_{1},\theta_{2},\theta_{3},\lambda,\overline{m},\overline{M})$.

Remark 4.6. Since the proof of Lemma 4.5 only uses the aforementioned Proposition 2.10 and Lemma 4.2 instead of Lemma 4.3, we can obtain a larger range for the value of $\theta_{3}+\theta_{4}$ than that in Lemma 4.3.

Proof. Step 1. For $i\geq0$, denote $k_{i}=M_{a}-\frac{\varepsilon}{2^{i}}$ and

\begin{align*} A(k_{i},R;t)=B_{R}\cap\{u(\cdot,t) \gt k_{i}\},\qquad A(k_{i},R)=(B_{R}\times[t_{0},t_{0}+aR^{\theta_{1}+\theta_{2}}])\cap\{u \gt k_{i}\}. \end{align*}

It then follows from Equation (2.13) that for $1 \lt q \lt 2,$

(4.11)\begin{align} &(k_{i+1}-k_{i})^{q}|A(k_{i+1},R;t)|_{\mu_{w_{1}}}^{q}|B_{R}\setminus A(k_{i},R;t)|_{\mu_{w_{1}}}\notag\\ &\leq CR^{q(n+\theta_{1}+\theta_{2}+1)}\int_{A(k_{i},R;t)\setminus A(k_{i+1},R;t)}|\nabla u|^{q}w_{1}\,{\rm d}x. \end{align}

From the assumed condition, we have

\begin{align*} |B_{R}\setminus A(k_{i},R;t)|_{\mu_{w_{1}}}\geq\gamma|B_{R}|_{\mu_{w_{1}}}=C(n,\theta_{1})\gamma R^{n+\theta_{1}+\theta_{2}}. \end{align*}

Substituting this into Equation (4.11) and integrating from t 0 to $t_{0}+a R^{\theta_{1}+\theta_{2}}$, we deduce from Hölder’s inequality that

\begin{align*} &\int^{t_{0}+aR^{\theta_{1}+\theta_{2}}}_{t_{0}}|A(k_{i+1},R;t)|_{\mu_{w_{1}}}\,{\rm d}t\notag\\ &\leq\frac{C2^{i+1}}{\varepsilon\sqrt[q]{\gamma}}R^{\frac{(n+\theta_{1}+\theta_{2})(q-1)}{q}+1}\int_{t_{0}}^{t_{0}+a R^{\theta_{1}+\theta_{2}}}\left(\int_{A(k_{i},R;t)\setminus A(k_{i+1},R;t)}|\nabla u|^{q}w_{1}\,{\rm d}x\right)^{\frac{1}{q}}\,{\rm d}t\notag\\ &\leq\frac{C2^{i+1}a^{\frac{q-1}{q}}}{\varepsilon\sqrt[q]{\gamma}}R^{\frac{(n+2\theta_{1}+2\theta_{2})(q-1)}{q}+1}\left(\int_{A(k_{i},R)\setminus A(k_{i+1},R)}|\nabla u|^{q}w_{1}\,{\rm d}x\,{\rm d}t\right)^{\frac{1}{q}}. \end{align*}

In light of $1 \lt q \lt 2$, it follows from Hölder’s inequality again that

\begin{align*} &\left(\int_{A(k_{i},R)\setminus A(k_{i+1},R)}|\nabla u|^{q}w_{1}\right)^{\frac{1}{q}}\notag\\ &\leq\left(\int_{A(k_{i},R)\setminus A(k_{i+1},R)}|\nabla u|^{2}w_{2}\right)^{\frac{1}{2}}\left(\int_{A(k_{i},R)\setminus A(k_{i+1},R)}|x'|^{\frac{2\theta_{1}-q\theta_{3}}{2-q}}|x|^{\frac{2\theta_{2}-q\theta_{4}}{2-q}}\right)^{\frac{2-q}{2q}}\notag\\ &\leq R^{\frac{\theta_{1}+\theta_{2}-\theta_{3}-\theta_{4}}{2}}|A(k_{i},R)\setminus A(k_{i+1},R)|_{\nu_{w_{1}}}^{\frac{2-q}{2q}}\left(\int_{B_{R}\times[t_{0},t_{0}+aR^{\theta_{1}+\theta_{2}}]}|\nabla (u-k_{i})^{+}|^{2}w_{2}\right)^{\frac{1}{2}}. \end{align*}

Pick a cutoff function $\eta\in C_{0}^{\infty}(B_{2R})$ such that

(4.12)\begin{align} \eta=1\;\mathrm{in}\;B_{R},\quad0\leq\eta\leq1,\;|\nabla\eta|\leq\frac{C(n)}{R}\;\,\mathrm{in}\;B_{2R}. \end{align}

Then from Lemma 4.2, we deduce

\begin{align*} &\left(\int_{B_{R}\times[t_{0},t_{0}+aR^{\theta_{1}+\theta_{2}}]}|\nabla (u-k_{i})^{+}|^{2}w_{2}\right)^{\frac{1}{2}}\notag\\ &\leq C\left(\int_{B_{2R}}|(u-k_{i})^{+}(x,t_{0})|^{2}w_{1}+\frac{1}{R^{2}}\int_{B_{2R}\times[t_{0},t_{0}+aR^{\theta_{1}+\theta_{2}}]}|(u-k_{i})^{+}|^{2}w_{2}\right)^{\frac{1}{2}}\notag\\ &\leq\frac{C\varepsilon}{2^{i}}R^{\frac{n+\theta_{1}+\theta_{2}+\theta_{3}+\theta_{4}-2}{2}}, \end{align*}

where in the last inequality, we used the assumed condition that $0\leq\theta_{3}+\theta_{4}\leq2$. Therefore, combining these above facts, we obtain

\begin{align*} |A(k_{i+1},R)|_{\nu_{w_{1}}}\leq&\frac{Ca^{\frac{q-1}{q}}}{\sqrt[q]{\gamma}}R^{\frac{(n+2\theta_{1}+2\theta_{2})(3q-2)}{2q}}\left|A(k_{i},R)\setminus A(k_{i+1},R)\right|_{\nu_{w_{1}}}^{\frac{2-q}{2q}}, \end{align*}

which yields that for $j\geq1$,

\begin{align*} j|A(k_{j},R)|_{\nu_{w_{1}}}^{\frac{2q}{2-q}}\leq&\sum^{j-1}_{i=0}|A(k_{i+1},R)|_{\nu_{w_{1}}}^{\frac{2q}{2-q}}\notag\\ \leq&\frac{Ca^{\frac{2(q-1)}{2-q}}}{\gamma^{\frac{2}{2-q}}}R^{\frac{(n+2\theta_{1}+2\theta_{2})(3q-2)}{2-q}}|B_{R}\times[t_{0},t_{0}+aR^{\theta_{1}+\theta_{2}}]|_{\nu_{w_{1}}}\notag\\ \leq&\frac{C}{\gamma^{\frac{2}{2-q}}a^{\frac{q}{2-q}}}|B_{R}\times[t_{0},t_{0}+aR^{\theta_{1}+\theta_{2}}]|_{\nu_{w_{1}}}^{\frac{2q}{2-q}}. \end{align*}

Then Equation (4.10) holds.

Step 2. For $i\geq0$, write

\begin{align*} \tilde{k}_{i}=m_{a}+\frac{\varepsilon}{2^{i}} \end{align*}

and

\begin{align*} \tilde{A}(k_{i},R;t)=B_{R}\cap\{u(\cdot,t) \lt k_{i}\},\qquad \tilde{A}(k_{i},R)=(B_{R}\times[t_{0},t_{0}+aR^{\theta_{1}+\theta_{2}}])\cap\{u \lt k_{i}\}. \end{align*}

Using Equation (2.14), we obtain that for $1 \lt q \lt 2,$

\begin{align*} &(\tilde{k}_{i}-\tilde{k}_{i+1})^{q}|\tilde{A}(\tilde{k}_{i+1},R;t)|_{\mu_{w_{1}}}^{q}|B_{R}\setminus \tilde{A}(\tilde{k}_{i},R;t)|_{\mu_{w_{1}}}\notag\\ &\leq CR^{q(n+\theta_{1}+\theta_{2}+1)}\int_{\tilde{A}(\tilde{k}_{i},R;t)\setminus \tilde{A}(\tilde{k}_{i+1},R;t)}|\nabla u|^{q}w_{1}\,{\rm d}x. \end{align*}

Observe by the assumed condition that

\begin{align*} |B_{R}\setminus \tilde{A}(\tilde{k}_{i},R;t)|_{\mu_{w_{1}}}\geq\gamma|B_{R}|_{\mu_{w_{1}}}=C(n,\theta_{1},\theta_{2})\gamma R^{n+\theta_{1}+\theta_{2}}. \end{align*}

Analogously as above, integrating from t 0 to $t_{0}+a R^{\theta_{1}+\theta_{2}}$ and using Hölder’s inequality, we have

\begin{align*} &\int^{t_{0}+aR^{\theta_{1}+\theta_{2}}}_{t_{0}}|\tilde{A}(\tilde{k}_{i+1},R;t)|_{\mu_{w_{1}}}\,{\rm d}t\notag\\ &\leq\frac{C2^{i+1}a^{\frac{q-1}{q}}}{\varepsilon\sqrt[q]{\gamma}}R^{\frac{(n+2\theta_{1}+2\theta_{2})(q-1)}{q}+1}\left(\int_{\tilde{A}(\tilde{k}_{i},R)\setminus \tilde{A}(\tilde{k}_{i+1},R)}|\nabla u|^{q}w_{1}\,{\rm d}x\,{\rm d}t\right)^{\frac{1}{q}} \end{align*}

and

\begin{align*} &\int_{\tilde{A}(\tilde{k}_{i},R)\setminus A(\tilde{k}_{i+1},R)}|\nabla u|^{q}w_{1}\notag\\ &\leq\left(\int_{\tilde{A}(\tilde{k}_{i},R)\setminus \tilde{A}(\tilde{k}_{i+1},R)}|\nabla u|^{2}w_{2}\right)^{\frac{q}{2}}\left(\int_{\tilde{A}(\tilde{k}_{i},R)\setminus \tilde{A}(\tilde{k}_{i+1},R)}|x'|^{\frac{2\theta_{1}-q\theta_{3}}{2-q}}|x|^{\frac{2\theta_{2}-q\theta_{4}}{2-q}}\right)^{\frac{2-q}{2}}\notag\\ &\leq R^{\frac{q(\theta_{1}+\theta_{2}-\theta_{3}-\theta_{4})}{2}}|\tilde{A}(\tilde{k}_{i},R)\setminus \tilde{A}(\tilde{k}_{i+1},R)|_{\nu_{w_{1}}}^{\frac{2-q}{2q}}\left(\int_{B_{R}\times[t_{0},t_{0}+aR^{\theta_{1}+\theta_{2}}]}|\nabla (u-\tilde{k}_{i})^{-}|^{2}w_{2}\right)^{\frac{q}{2}}. \end{align*}

For any $0 \lt \varepsilon\leq\varepsilon_{0}=C_{0}^{-1}$ with C 0 given in Lemma 4.2, we know that

\begin{align*} [(u-\tilde{k}_{i})^{-}]^{2}-C_{0}[(u-\tilde{k}_{i})^{-}]^{3}\geq(1-C_{0}\varepsilon)[(u-\tilde{k}_{i})^{-}]^{2}\geq0. \end{align*}

Therefore, in view of $0\leq\theta_{3}+\theta_{4}\leq2$ and applying Lemma 4.2 with η defined by Equation (4.12), we derive

\begin{align*} &\int_{B_{R}\times[t_{0},t_{0}+aR^{\theta_{1}+\theta_{2}}]}|\nabla (u-\tilde{k}_{i})^{-}|^{2}w_{2}\notag\\ &\leq C\left(\int_{B_{2R}}|(u-\tilde{k}_{i})^{-}(x,t_{0})|^{2}w_{1}+\frac{1}{R^{2}}\int_{B_{R}\times[t_{0},t_{0}+aR^{\theta_{1}+\theta_{2}}]}|(u-\tilde{k}_{i})^{-}|^{2}w_{2}\right)\notag\\ &\leq\frac{C\varepsilon^{2}}{4^{i}}R^{n+\theta_{1}+\theta_{2}+\theta_{3}+\theta_{4}-2}. \end{align*}

Then we deduce

\begin{align*} |\tilde{A}(\tilde{k}_{i+1},R)|_{\nu_{w_{1}}}^{\frac{2q}{2-q}}\leq&\frac{Ca^{\frac{2(q-1)}{2-q}}}{\gamma^{\frac{2}{2-q}}}R^{\frac{(n+2\theta_{1}+2\theta_{2})(3q-2)}{2-q}}|\tilde{A}(\tilde{k}_{i},R)\setminus \tilde{A}(\tilde{k}_{i+1},R)|_{\nu_{w_{1}}}. \end{align*}

This leads to that for $j\geq1$,

\begin{align*} j|\tilde{A}(\tilde{k}_{j},R)|_{\nu_{w_{1}}}^{\frac{2q}{2-q}}\leq&\sum^{j-1}_{i=0}|\tilde{A}(\tilde{k}_{i+1},R)|_{\nu_{w_{1}}}^{\frac{2q}{2-q}}\notag\\ \leq&\frac{Ca^{\frac{2(q-1)}{2-q}}}{\gamma^{\frac{2}{2-q}}}R^{\frac{(n+2\theta_{1}+2\theta_{2})(3q-2)}{2-q}}|B_{R}\times[t_{0},t_{0}+aR^{\theta_{1}+\theta_{2}}]|_{\nu_{w_{1}}}\notag\\ \leq&\frac{C}{\gamma^{\frac{2}{2-q}}a^{\frac{q}{2-q}}}|B_{R}\times[t_{0},t_{0}+aR^{\theta_{1}+\theta_{2}}]|_{\nu_{w_{1}}}^{\frac{2q}{2-q}}. \end{align*}

The proof is complete.

We now give explicit estimates for the distribution function of u at each time slice from the starting time.

Lemma 4.7. Assume as in Theorem 1.6 or Theorem 1.9. Let $0 \lt \gamma \lt 1$, $0 \lt R \lt \frac{1}{2}$, $-\frac{1}{2} \lt t_{0}\leq-R^{\theta_{1}+\theta_{2}}$ and $\overline{m}\leq m_{1}\leq\inf\limits_{B_{2R}\times[t_{0},t_{0}+R^{\theta_{1}+\theta_{2}}]}u\leq\sup\limits_{B_{2R}\times[t_{0},t_{0}+R^{\theta_{1}+\theta_{2}}]}u\leq M_{1}\leq\overline{M}$. Then there exist a small constant $\bar{\varepsilon}_{0}=\bar{\varepsilon}_{0}(n,p,\theta_{1},\theta_{2},\lambda,\gamma,\overline{m},\overline{M}) \gt 0$ and a large constant $\bar{l}_{0}=\bar{l}_{0}(n,p,q,\theta_{1},\theta_{2},\theta_{3},\lambda,\gamma,\overline{m},\overline{M}) \gt 1$ such that

  1. (i) for every $0 \lt \varepsilon\leq\bar{\varepsilon}_{0}$, if

    (4.13)\begin{align} \frac{|\{x\in B_{R}: u(x,t_{0}) \gt M_{1}-\varepsilon\}|_{\mu_{w_{1}}}}{|B_{R}|_{\mu_{w_{1}}}}\leq1-\gamma, \end{align}
    then for any $t_{0}\leq t\leq t_{0}+R^{\theta_{1}+\theta_{2}},$
    (4.14)\begin{align} \frac{|\{x\in B_{R}: u(x,t) \gt M_{1}-2^{-l_{0}}\varepsilon\}|_{\mu_{w_{1}}}}{|B_{R}|_{\mu_{w_{1}}}}\leq1-\frac{\gamma}{2}; \end{align}
  2. (ii) for every $0 \lt \varepsilon\leq\bar{\varepsilon}_{0}$, if

    (4.15)\begin{align} \frac{|\{x\in B_{R}: u(x,t_{0}) \lt m_{1}+\varepsilon\}|_{\mu_{w_{1}}}}{|B_{R}|_{\mu_{w_{1}}}}\leq1-\gamma, \end{align}
    then for any $t_{0}\leq t\leq t_{0}+R^{\theta_{1}+\theta_{2}},$
    (4.16)\begin{align} \frac{|\{x\in B_{R}: u(x,t) \lt m_{1}+2^{-l_{0}}\varepsilon\}|_{\mu_{w_{1}}}}{|B_{R}|_{\mu_{w_{1}}}}\leq1-\frac{\gamma}{2}. \end{align}

Remark 4.8. It is worth emphasizing that the explicit values of $\bar{\varepsilon}_{0}$ and $\bar{l}_{0}$ are given by Equations (4.20) and (4.22) below.

Proof. Step 1. For $a\in(0,1]$ and $k\in[\overline{m},\overline{M}]$, define

\begin{align*} A^{a}(k,R)=(B_{R}\times[t_{0},t_{0}+a R^{\theta_{1}+\theta_{2}}])\cap\{u \gt k\}. \end{align*}

Take a smooth cutoff function $\eta\in C^{\infty}_{0}(B_{R})$ satisfying that η = 1 in $B_{\sigma R}$, where $\sigma\in(0,1)$ to be determined later. Set $k_{1} \gt 1$. Denote $v=(u-(M_{1}-\varepsilon))^{+}$. From Lemma 4.2, we obtain

(4.17)\begin{align} &\sup\limits_{t\in(t_{0},t_{0}+aR^{\theta_{1}+\theta_{2}})}\int_{B_{R}}v^{2}\eta^{2}w_{1}\,{\rm d}x\notag\\ &\leq\int_{B_{R}}(v^{2}+Cv^{3})\eta^{2}w_{1}\,{\rm d}x\big|_{t_{0}}+C\int_{B_{R}\times[t_{0},t_{0}+aR^{\theta_{1}+\theta_{2}}]}v^{2}|\nabla\eta|^{2}w_{2}\,{\rm d}x\,{\rm d}t. \end{align}

Observe that for $t\in[t_{0},t_{0}+aR^{\theta_{1}+\theta_{2}}]$,

\begin{align*} \int_{B_{R}}v^{2}\eta^{2}w_{1}\,{\rm d}x\big|_{t}\geq\varepsilon^{2}(1-2^{-k_{1}})^{2}|B_{\sigma R}\cap\{u(x,t) \gt M_{1}-2^{-k_{1}}\varepsilon\}|_{\mu_{w_{1}}}, \end{align*}

and by Equation (4.13),

\begin{align*} \int_{B_{R}}(v^{2}+Cv^{3})\eta^{2}w_{1}\,{\rm d}x\big|_{t_{0}}\leq&\varepsilon^{2}(1+C\varepsilon)|\{x\in B_{R}: u(x,t_{0}) \gt M_{1}-\varepsilon\}|_{\mu_{w_{1}}}\notag\\ \leq&\varepsilon^{2}(1+C\varepsilon)(1-\gamma)|B_{R}|_{\mu_{w_{1}}}, \end{align*}

and

\begin{align*} \int^{t_{0}+aR^{\theta_{1}+\theta_{2}}}_{t_{0}}\int_{B_{R}}v^{2}|\nabla\eta|^{2}w_{2}\,{\rm d}x\,{\rm d}t\leq&\frac{C\varepsilon^{2}}{(1-\sigma)^{2}R^{2}}|A^{a}(M_{1}-\varepsilon,R)|_{\nu_{w_{2}}}\notag\\ \leq&\frac{C\varepsilon^{2}R^{\theta_{3}+\theta_{4}-2}}{(1-\sigma)^{2}}|B_{R}|_{\mu_{w_{1}}}\frac{|A^{a}(M_{1}-\varepsilon,R)|_{\nu_{w_{2}}}}{|Q_{R}|_{\nu_{w_{2}}}}\notag\\ =&\frac{C\varepsilon^{2}}{(1-\sigma)^{2}}|B_{R}|_{\mu_{w_{1}}}\frac{|A^{a}(M_{1}-\varepsilon,R)|_{\nu_{w_{2}}}}{|Q_{R}|_{\nu_{w_{2}}}}, \end{align*}

where we utilized the assumed condition of $\theta_{3}+\theta_{4}=2$. A consequence of these facts gives that for $t\in[t_{0},t_{0}+aR^{\theta_{1}+\theta_{2}}]$,

\begin{align*} &|B_{\sigma R}\cap\{u(x,t) \gt M_{1}-2^{-k_{1}}\varepsilon\}|_{\mu_{w_{1}}}\notag\\ &\leq|B_{R}|_{\mu_{w_{1}}}\left(\frac{(1+C\varepsilon)(1-\gamma)}{(1-2^{-k_{1}})^{2}}+\frac{C}{(1-\sigma)^{2}}\frac{|A^{a}(M_{1}-\varepsilon,R)|_{\nu_{w_{2}}}}{|Q_{R}|_{\nu_{w_{2}}}}\right), \end{align*}

which, together with the fact that $|B_{R}\setminus B_{\sigma R}|_{\mu_{w_{1}}}\leq C(1-\sigma)|B_{R}|_{\mu_{w_{1}}}$, reads that

\begin{align*} &\frac{|B_{R}\cap\{u(x,t) \gt M_{1}-2^{-k_{1}}\varepsilon\}|_{\mu_{w_{1}}}}{|B_{R}|_{\mu_{w_{1}}}}\notag\\ &\leq \frac{(1+C\varepsilon)(1-\gamma)}{(1-2^{-k_{1}})^{2}}+\frac{C}{(1-\sigma)^{2}}\bigg((1-\sigma)^{3}+\frac{|A^{a}(M_{1}-\varepsilon,R)|_{\nu_{w_{2}}}}{|Q_{R}|_{\nu_{w_{2}}}}\bigg). \end{align*}

Pick σ such that

\begin{align*} (1-\sigma)^{3}=\frac{|A^{a}(M_{1}-\varepsilon,R)|_{\nu_{w_{2}}}}{|Q_{R}|_{\nu_{w_{2}}}}, \end{align*}

which yields that for $t\in[t_{0},t_{0}+aR^{\theta_{1}+\theta_{2}}],$

(4.18)\begin{align} &\frac{|B_{R}\cap\{u(x,t) \gt M_{1}-2^{-k_{1}}\varepsilon\}|_{\mu_{w_{1}}}}{|B_{R}|_{\mu_{w_{1}}}}\notag\\ &\leq \frac{(1+\overline{C}\varepsilon)(1-\gamma)}{(1-2^{-k_{1}})^{2}}+\overline{C}\bigg(\frac{|A^{a}(M_{1}-\varepsilon,R)|_{\nu_{w_{2}}}}{|Q_{R}|_{\nu_{w_{2}}}}\bigg)^{\frac{1}{3}}, \end{align}

where $\overline{C}=\overline{C}(n,p,\theta_{1},\theta_{2},\lambda,\overline{m},\overline{M}).$ Note that

\begin{align*} \frac{|A^{a}(M_{1}-\varepsilon,R)|_{\nu_{w_{2}}}}{|Q_{R}|_{\nu_{w_{2}}}}\leq a. \end{align*}

Take a small positive constant a such that a −1 is an integer and

\begin{align*} \overline{C}a^{\frac{1}{3}}\leq\frac{\gamma}{8}. \end{align*}

By fixing the value of a, we now divide the time interval $[t_{0},t_{0}+R^{\theta_{1}+\theta_{2}}]$ into finite small intervals. Denote $N=a^{-1}$ and $t_{i}=t_{0}+iaR^{\theta_{1}+\theta_{2}}$, $i=1,2,\ldots,N.$

Claim that there exist a small positive constant $\bar{\varepsilon}_{0}$ and a large positive constant $k_{0} \gt 1$ depending only on $n,p,\theta_{1},\theta_{2},\lambda,\gamma,\overline{m},\overline{M}$ such that for any $0 \lt \varepsilon\leq\bar{\varepsilon}_{0}$ and $k_{1}\geq k_{0}$,

(4.19)\begin{align} \frac{(1+\overline{C}\varepsilon)(1-\gamma)}{(1-2^{-k_{1}})^{2}}\leq1-\gamma+\frac{\gamma}{8N}. \end{align}

In fact, since $(1-t)^{-2}\leq(1+6t)$ for $t\in(0,\frac{1}{2})$, then

\begin{align*} \frac{1+\overline{C}\varepsilon}{(1-2^{-k_{1}})^{2}}\leq(1+\overline{C}\varepsilon)(1+6\cdot2^{-k_{1}}). \end{align*}

Let $\overline{C}\varepsilon=6\cdot2^{-k_{1}}$. Then we have

\begin{align*} \frac{1+\overline{C}\varepsilon}{(1-2^{-k_{1}})^{2}}\leq1+2\overline{C}\varepsilon+\overline{C}^{2}\varepsilon^{2}. \end{align*}

Pick

(4.20)\begin{align} \bar{\varepsilon}_{0}=\frac{-2\overline{C}+\sqrt{4\overline{C}^{2}+\frac{\gamma}{2N(1-\gamma)}}}{2\overline{C}^{2}},\qquad k_{0}=-\frac{\ln2}{\ln(\overline{C}\bar{\varepsilon}_{0})-\ln6}. \end{align}

Then we obtain that for any $0 \lt \varepsilon\leq\bar{\varepsilon}_{0}$ and $k_{1}\geq k_{0}$,

\begin{align*} \frac{1+\overline{C}\varepsilon}{(1-2^{-k_{1}})^{2}}\leq&1+2\overline{C}\bar{\varepsilon}_{0}+\overline{C}^{2}\bar{\varepsilon}^{2}_{0}=1+\frac{\gamma}{8N(1-\gamma)}. \end{align*}

That is, Equation (4.19) holds.

Consequently, it follows from Equation (4.18) that for $0 \lt \varepsilon\leq\bar{\varepsilon}_{0}$, $k_{1}\geq k_{0}$ and $t\in[t_{0},t_{1}]$,

\begin{align*} &\frac{|B_{R}\cap\{u(x,t) \gt M_{1}-2^{-k_{1}}\varepsilon\}|_{\mu_{w_{1}}}}{|B_{R}|_{\mu_{w_{1}}}}\leq1-\left(\frac{7}{8}-\frac{1}{8N}\right)\gamma. \end{align*}

Then applying Lemma 4.5, we deduce from Equation (4.3) that for any $k_{2} \gt k_{1}\geq k_{0}$,

\begin{align*} \frac{|A^{a}(M_{1}-2^{-k_{2}}\varepsilon,R)|_{\nu_{w_{2}}}}{|Q_{R}|_{\nu_{w_{2}}}}\leq& C\left(\frac{|A^{a}(M_{1}-2^{-k_{2}}\varepsilon,R)|_{\nu_{w_{1}}}}{|Q_{R}|_{\nu_{w_{1}}}}\right)^{\frac{\theta_{3}}{\theta_{1}}}\notag\\ \leq& \widehat{C}\left(\frac{\sqrt{a}}{\sqrt[q]{\gamma}(k_{2}-k_{1})^{\frac{2-q}{2q}}}\right)^{\frac{\theta_{3}}{\theta_{1}}}, \end{align*}

where $\widehat{C}=\widehat{C}(n,p,q,\theta_{1},\theta_{2},\theta_{3},\lambda,\overline{m},\overline{M}).$ Pick

\begin{align*} k_{2}=k_{1}+a^{\frac{q}{2-q}}\gamma^{-\frac{2}{2-q}}\left(\frac{\gamma}{8N\overline{C}\sqrt[3]{\widehat{C}}}\right)^{-\frac{6\theta_{1}q}{\theta_{3}(2-q)}}. \end{align*}

Then we have

(4.21)\begin{align} \overline{C}\left(\frac{|A^{a}(M_{1}-2^{-k_{2}}\varepsilon,R)|_{\nu_{w_{2}}}}{|Q_{R}|_{\nu_{w_{2}}}}\right)^{\frac{1}{3}}\leq\frac{\gamma}{8N}. \end{align}

Choose $k_{1}=k_{0}$ and $l_{1}=k_{1}+k_{2}$. By letting $2^{-k_{2}}\varepsilon$ substitute for ɛ in Equation (4.18), we have

\begin{align*} \sup\limits_{t\in[t_{0},t_{1}]}|B_{R}\cap\{u(x,t) \gt M_{1}-2^{-l_{1}}\varepsilon\}|_{\mu_{w_{1}}}\leq\Big(1-\gamma+\frac{\gamma}{4N}\Big)|B_{R}|_{\mu_{w_{1}}}. \end{align*}

Then it can be inductively proved that there exist a strictly increasing integer set $\{l_{i}\}_{i=1}^{N}$ such that for $i=1,2,\ldots,N$,

\begin{align*} \sup\limits_{t\in[t_{i-1},t_{i}]}|B_{R}\cap\{u(x,t) \gt M_{1}-2^{-l_{i}}\varepsilon\}|_{\mu_{w_{1}}}\leq\Big(1-\gamma+\frac{i\gamma}{4N}\Big)|B_{R}|_{\mu_{w_{1}}}. \end{align*}

In fact, let the above relation hold in interval $[t_{i-1},t_{i}]$ and then prove that it also holds in the next interval $[t_{i},t_{i+1}]$. For simplicity, denote $\varepsilon_{i}=2^{-l_{i}}\varepsilon$ and $\gamma_{i}=\gamma(1-\frac{i}{4N}).$ Then the assumption implies that

\begin{align*} |B_{R}\cap\{u(x,t_{i}) \gt M_{1}-\varepsilon_{i}\}|_{\mu_{w_{1}}}\leq(1-\gamma_{i})|B_{R}|_{\mu_{w_{1}}}. \end{align*}

By the same argument as in Equation (4.18), it follows from Equations (4.19)–(4.21) that for $\bar{k}_{1}\geq k_{0}$ and $t\in[t_{i},t_{i+1}],$

\begin{align*} &\frac{|B_{R}\cap\{u(x,t) \gt M_{1}-2^{-\bar{k}_{1}}\varepsilon_{i}\}|_{\mu_{w_{1}}}}{|B_{R}|_{\mu_{w_{1}}}}\notag\\ &\leq \frac{(1+\overline{C}\varepsilon_{i})(1-\gamma_{i})}{(1-2^{-\bar{k}_{1}})^{2}}+\overline{C}\bigg(\frac{|A^{a}(M_{1}-\varepsilon_{i},R)|_{\nu_{w_{2}}}}{|Q_{R}|_{\nu_{w_{2}}}}\bigg)^{\frac{1}{3}}\notag\\ &\leq1-\gamma_{i}+\frac{\gamma_{i}}{8N}+\overline{C}\bigg(\frac{|A^{a}(M_{1}-2^{-l_{i}}\varepsilon,R)|_{\nu_{w_{2}}}}{|Q_{R}|_{\nu_{w_{2}}}}\bigg)^{\frac{1}{3}}\notag\\ &\leq1-\gamma_{i}+\frac{\gamma_{i}}{8N}+\frac{\gamma}{8N}\notag\\ & \lt 1-\gamma+\frac{i+1}{4N}\gamma, \end{align*}

where $\overline{C}=\overline{C}(n,p,\theta_{1},\theta_{2},\lambda,\overline{m},\overline{M})$ is defined above and in the third inequality, we used the fact that $l_{i}\geq l_{1} \gt k_{2}$. By taking $\bar{k}_{1}=k_{0}$ and $l_{i+1}=l_{i}+\bar{k}_{1}$, we obtain

\begin{align*} \sup\limits_{t\in[t_{i},t_{i+1}]}|B_{R}\cap\{u(x,t) \gt M_{1}-2^{-l_{i+1}}\varepsilon\}|_{\mu_{w_{1}}}\leq\Big(1-\gamma+\frac{(i+1)\gamma}{4N}\Big)|B_{R}|_{\mu_{w_{1}}}. \end{align*}

Then picking

(4.22)\begin{align} \bar{l}_{0}:=&l_{N}=l_{1}+(N-1)k_{0}\notag\\ =&(N+1)k_{0}++a^{\frac{q}{2-q}}\gamma^{-\frac{2}{2-q}}\left(\frac{\gamma}{8N\overline{C}\sqrt[3]{\widehat{C}}}\right)^{-\frac{6\theta_{1}q}{\theta_{3}(2-q)}}, \end{align}

we obtain that Equation (4.14) holds.

Step 2. For $0 \lt a\leq1$ and $\overline{m}\leq k\leq\overline{M}$, let

\begin{align*} \tilde{A}^{a}(k,R)=(B_{R}\times[t_{0},t_{0}+a R^{\theta_{1}+\theta_{2}}])\cap\{u \lt k\}. \end{align*}

Define $\tilde{v}=(u-(m_{1}+\varepsilon))^{-}$. A direct application of Lemma 4.2 gives that

\begin{align*} &\sup\limits_{t\in(t_{0},t_{0}+aR^{\theta_{1}+\theta_{2}})}\int_{B_{R}}(\tilde{v}^{2}-C_{0}\tilde{v}^{3})\eta^{2}w_{1}\,{\rm d}x\notag\\ &\leq\int_{B_{R}}\tilde{v}^{2}\eta^{2}w_{1}\,{\rm d}x\big|_{t_{0}}+C_{0}\int_{B_{R}\times[t_{0},t_{0}+aR^{\theta_{1}+\theta_{2}}]}\tilde{v}^{2}|\nabla\eta|^{2}w_{2}\,{\rm d}x\,{\rm d}t, \end{align*}

where $C_{0}=C_{0}(n,p,\lambda,\overline{m},\overline{M})$ and η is defined in Equation (4.17). Pick a small constant $0 \lt \bar{\varepsilon}_{1}\leq(2C_{0})^{-1}$, which implies that $1-C_{0}\bar{\varepsilon}_{1}\geq\frac{1}{2}$. Then we obtain that for $t_{0} \lt t \lt t_{0}+aR^{\theta_{1}+\theta_{2}}$, $0 \lt \varepsilon\leq\bar{\varepsilon}_{0}$ and $k_{1} \gt 1$,

\begin{align*} \int_{B_{R}}(\tilde{v}^{2}-C_{0}\tilde{v}^{3})\eta^{2}w_{1}\,{\rm d}x\big|_{t}\geq(1-C_{0}\varepsilon)\varepsilon^{2}(1-2^{-k_{1}})^{2}|B_{\sigma R}\cap\{u(x,t) \lt m_{1}+2^{-k_{1}}\varepsilon\}|_{\mu_{w_{1}}}, \end{align*}

and in view of Equation (4.15),

\begin{align*} \int_{B_{R}}\tilde{v}^{2}\eta^{2}w_{1}\,{\rm d}x\big|_{t_{0}}\leq&\varepsilon^{2}|\{x\in B_{R}: u(x,t_{0}) \lt m_{1}+\varepsilon\}|_{\mu_{w_{1}}}\leq\varepsilon^{2}(1-\gamma)|B_{R}|_{\mu_{w_{1}}}, \end{align*}

and

\begin{align*} \int_{B_{R}\times[t_{0},t_{0}+aR^{\theta_{1}+\theta_{2}}]}\tilde{v}^{2}|\nabla\eta|^{2}w_{2}\,{\rm d}x\,{\rm d}t\leq&\frac{C\varepsilon^{2}}{(1-\sigma)^{2}R^{2}}|\tilde{A}^{a}(m_{1}+\varepsilon,R)|_{\nu_{w_{2}}}\notag\\ \leq&\frac{C\varepsilon^{2}R^{\theta_{3}+\theta_{4}-2}}{(1-\sigma)^{2}}|B_{R}|_{\mu_{w_{1}}}\frac{|\tilde{A}^{a}(m_{1}+\varepsilon,R)|_{\nu_{w_{2}}}}{|Q_{R}|_{\nu_{w_{2}}}}\notag\\ =&\frac{C\varepsilon^{2}}{(1-\sigma)^{2}}|B_{R}|_{\mu_{w_{1}}}\frac{|\tilde{A}^{a}(m_{1}+\varepsilon,R)|_{\nu_{w_{2}}}}{|Q_{R}|_{\nu_{w_{2}}}}. \end{align*}

Therefore, we deduce that for $t_{0}\leq t\leq t_{0}+aR^{\theta_{1}+\theta_{2}}$,

\begin{align*} &|B_{\sigma R}\cap\{u(x,t) \lt m_{1}+2^{-k_{1}}\varepsilon\}|_{\mu_{w_{1}}}\notag\\ &\leq|B_{R}|_{\mu_{w_{1}}}\left(\frac{1-\gamma}{(1-C\varepsilon)(1-2^{-k_{1}})^{2}}+\frac{C}{(1-\sigma)^{2}}\frac{|\tilde{A}^{a}(m_{1}+\varepsilon,R)|_{\nu_{w_{2}}}}{|Q_{R}|_{\nu_{w_{2}}}}\right)\notag\\ &\leq|B_{R}|_{\mu_{w_{1}}}\left(\frac{(1+C\varepsilon)(1-\gamma)}{(1-2^{-k_{1}})^{2}}+\frac{C}{(1-\sigma)^{2}}\frac{|\tilde{A}^{a}(m_{1}+\varepsilon,R)|_{\nu_{w_{2}}}}{|Q_{R}|_{\nu_{w_{2}}}}\right), \end{align*}

and thus,

\begin{align*} &\frac{|B_{R}\cap\{u(x,t) \lt m_{1}+2^{-k_{1}}\varepsilon\}|_{\mu_{w_{1}}}}{|B_{R}|_{\mu_{w_{1}}}}\notag\\ &\leq \frac{(1+C\varepsilon)(1-\gamma)}{(1-2^{-k_{1}})^{2}}+\frac{C}{(1-\sigma)^{2}}\bigg((1-\sigma)^{3}+\frac{|\tilde{A}^{a}(m_{1}+\varepsilon,R)|_{\nu_{w_{2}}}}{|Q_{R}|_{\nu_{w_{2}}}}\bigg). \end{align*}

Take σ such that

\begin{align*} (1-\sigma)^{3}=\frac{|\tilde{A}^{a}(m_{1}+\varepsilon,R)|_{\nu_{w_{2}}}}{|Q_{R}|_{\nu_{w_{2}}}}. \end{align*}

Then we obtain that for $t_{0}\leq t\leq t_{0}+aR^{\theta_{1}+\theta_{2}},$

\begin{align*} &\frac{|B_{R}\cap\{u(x,t) \lt m_{1}+2^{-k_{1}}\varepsilon\}|_{\mu_{w_{1}}}}{|B_{R}|_{\mu_{w_{1}}}}\notag\\ &\leq \frac{(1+\overline{C}\varepsilon)(1-\gamma)}{(1-2^{-k_{1}})^{2}}+\overline{C}\bigg(\frac{|\tilde{A}^{a}(m_{1}+\varepsilon,R)|_{\nu_{w_{2}}}}{|Q_{R}|_{\nu_{w_{2}}}}\bigg)^{\frac{1}{3}}, \end{align*}

where $\overline{C}=\overline{C}(n,p,\theta_{1},\theta_{2},\lambda,\overline{m},\overline{M}).$ Consequently, by the same argument as in the left proof of Equation (4.14) above, we deduce that Equation (4.16) holds. The proof is complete.

A consequence of Lemmas 4.3, 4.5 and 4.7 gives the improvement on oscillation of u in a small region.

Corollary 4.9. Assume as in Theorem 1.6 or Theorem 1.9. Let $0 \lt \gamma \lt 1$, $0 \lt R \lt \frac{1}{2}$, $-\frac{1}{4} \lt t_{0}\leq0$ and $\overline{m}\leq m\leq\inf\limits_{B_{2R}\times[t_{0}-R^{\theta_{1}+\theta_{2}},t_{0}]}u\leq\sup\limits_{B_{2R}\times[t_{0}-R^{\theta_{1}+\theta_{2}},t_{0}]}u\leq M\leq\overline{M}$. Then there exist a small constant $\tilde{\varepsilon}_{0}=\tilde{\varepsilon}_{0}(n,p,\theta_{1},\theta_{2},\lambda,\gamma,\overline{m},\overline{M}) \gt 0$ and a large constant $l_{0}=l_{0}(n,p,q,\theta_{1},\theta_{2},\theta_{3},\lambda,\gamma,\overline{m},\overline{M}) \gt 1$ such that for any $0 \lt \varepsilon\leq\tilde{\varepsilon}_{0}$,

  1. (i) if

    \begin{align*} \frac{|\{x\in B_{R}: u(x,t_{0}-R^{\theta_{1}+\theta_{2}}) \gt M-\varepsilon\}|_{\mu_{w_{1}}}}{|B_{R}|_{\mu_{w_{1}}}}\leq1-\gamma, \end{align*}
    then
    \begin{align*} \sup\limits_{Q_{R/2}(0,t_{0})}u\leq M-\frac{\varepsilon}{2^{l_{0}}}; \end{align*}
  2. (ii) if

    \begin{align*} \frac{|\{x\in B_{R}: u(x,t_{0}-R^{\theta_{1}+\theta_{2}}) \lt m+\varepsilon\}|_{\mu_{w_{1}}}}{|B_{R}|_{\mu_{w_{1}}}}\leq1-\gamma, \end{align*}
    then
    \begin{align*} \inf\limits_{Q_{R/2}(0,t_{0})}u\geq m+\frac{\varepsilon}{2^{l_{0}}}. \end{align*}

Proof. Applying Lemma 4.3, Lemma 4.5 with a = 1 and Lemma 4.7, we obtain that Corollary 4.9 holds. In particular, we fix $q=\frac{3}{2}$ under the assumed conditions in Theorem 1.9.

Based on these above facts, we now give the proofs of Theorems 1.6 and 1.9, respectively.

Proof of Theorem 1.6

Pick a sufficiently large constant $\kappa_{0}\geq2$ such that

\begin{align*} \frac{\overline{M}-\overline{m}}{\kappa_{0}} \lt \tilde{\varepsilon}_{0}, \end{align*}

where $\tilde{\varepsilon}_{0}$ is given by Corollary 4.9 with $\gamma=\frac{1}{2}$. For $0 \lt R\leq\frac{1}{2}$ and $-\frac{1}{4} \lt t_{0} \lt 0$, define

\begin{align*} \overline{\mu}(R)=\sup\limits_{(x,t)\in Q_{R}(0,t_{0})}u(x,t),\quad\underline{\mu}(R)=\inf\limits_{(x,t)\in Q_{R}(0,t_{0})}u(x,t),\quad \omega(R)=\overline{\mu}(R)-\underline{\mu}(R). \end{align*}

Observe that there is at least one inequality holding in terms of the following two inequalities:

(4.23)\begin{align} |\{x\in B_{R/2}:u(x,t_{0}-(R/2)^{\theta_{1}+\theta_{2}}) \gt \overline{\mu}(R)-\kappa_{0}^{-1}\omega(R)\}|_{\mu_{w_{1}}}\leq\frac{1}{2}|B_{R/2}|_{\mu_{w_{1}}}, \end{align}

and

(4.24)\begin{align} |\{x\in B_{R/2}:u(x,t_{0}-(R/2)^{\theta_{1}+\theta_{2}}) \lt \underline{\mu}(R)+\kappa_{0}^{-1}\omega(R)\}|_{\mu_{w_{1}}}\leq\frac{1}{2}|B_{R/2}|_{\mu_{w_{1}}}. \end{align}

From Corollary 4.9, it follows that there exists a large constant $l_{0} \gt 1$ such that

\begin{align*} \overline{\mu}(R/4)\leq\overline{\mu}(R)-\frac{\omega(R)}{\kappa_{0}2^{l_{0}}}\quad\text{if Equation~(4.23) holds,} \end{align*}

and

\begin{align*} \underline{\mu}(R/4)\geq\underline{\mu}(R)+\frac{\omega(R)}{\kappa_{0}2^{l_{0}}}\quad\text{if Equation~(4.24) holds.} \end{align*}

In both cases, we have

\begin{align*} \omega(R/4)\leq\left(1-\frac{1}{\kappa_{0}2^{l_{0}}}\right)\omega(R)=\frac{1}{4^{\alpha}}\omega(R),\quad\mathrm{with}\;\alpha=-\frac{\ln\big(1-\frac{1}{\kappa_{0}2^{l_{0}}}\big)}{\ln4}. \end{align*}

Note that for any $0 \lt R\leq\frac{1}{2}$, there is an integer k such that $4^{-(k+1)}\cdot2^{-1} \lt R\leq4^{-k}\cdot2^{-1}$. In light of the fact that $\omega(R)$ is non-decreasing in R, it follows that

\begin{align*} \omega(R)\leq\omega(4^{-k}\cdot2^{-1})\leq 4^{-k\alpha}\omega(2^{-1})=8^{\alpha}(4^{-(k+1)}\cdot2^{-1})^{\alpha}\omega(2^{-1})\leq CR^{\alpha}, \end{align*}

where $C=C(n,p,q,\theta_{1},\theta_{2},\theta_{3},\lambda,\overline{m},\overline{M})$. Therefore, for any $(x,t)\in B_{1/2}\times(-1/4,t_{0})$, we obtain that

  1. (i) if $|t-t_{0}|\leq2^{-(\theta_{1}+\theta_{2})}$, then

    \begin{align*} |u(x,t)-u(0,t_{0})|\leq&|u(x,t)-u(x,t_{0})|+|u(x,t_{0})-u(0,t_{0})|\\ &\leq C\left(|t-t_{0}|^{\frac{\alpha}{\theta_{1}+\theta_{2}}}+|x|^{\alpha}\right)\notag\\ \leq& C\left(|x|+|t-t_{0}|^{\frac{1}{\theta_{1}+\theta_{2}}}\right)^{\alpha}; \end{align*}
  2. (ii) if $|t-t_{0}| \gt 2^{-(\theta_{1}+\theta_{2})}$, there exists a set $\{t_{i}\}_{i=1}^{N}$ such that $t \lt t_{1}\leq\cdots\leq t_{N} \lt t_{0}$,

    \begin{align*} |u(x,t)-u(0,t_{0})|\leq&|u(x,t)-u(x,t_{1})|+|u(x,t_{1})-u(x,t_{0})|+|u(x,t_{0})-u(0,t_{0})|\notag\\ \leq&C\left(|t-t_{1}|^{\frac{\alpha}{\theta_{1}+\theta_{2}}}+|t_{1}-t_{0}|^{\frac{\alpha}{\theta_{1}+\theta_{2}}}+|x|^{\frac{\alpha}{\theta_{1}+\theta_{2}}}\right)\notag\\ \leq&C\left(|x|+|t-t_{0}|^{\frac{1}{\theta_{1}+\theta_{2}}}\right)^{\alpha}\quad\text{if}\;N=1, \end{align*}
    and
    \begin{align*} &|u(x,t)-u(0,t_{0})|\notag\\ &\leq|u(x,t)-u(x,t_{1})|+\sum^{N-1}_{i=1}|u(x,t_{i})-u(x,t_{i+1})|\notag\\ &\quad+|u(x,t_{N})-u(x,t_{0})|+|u(x,t_{0})-u(0,t_{0})|\notag\\ &\leq C\left(|t-t_{1}|^{\frac{\alpha}{\theta_{1}+\theta_{2}}}+\sum^{N-1}_{i=1}|t_{i}-t_{i+1}|^{\frac{\alpha}{\theta_{1}+\theta_{2}}}+|t_{N}-t_{0}|^{\frac{\alpha}{\theta_{1}+\theta_{2}}}+|x|^{\frac{\alpha}{\theta_{1}+\theta_{2}}}\right)\notag\\ &\leq C\left(|x|+|t-t_{0}|^{\frac{1}{\theta_{1}+\theta_{2}}}\right)^{\alpha}\quad \text{if}\;N\geq2. \end{align*}

    The proof is complete.

Proof of Theorem 1.9

To begin with, applying the aforementioned proof of Theorem 1.6 with minor modification, we also obtain that there exists a small constant $0 \lt \alpha \lt 1$ and a large constant C > 0, both depending only on $n,p,\theta_{2},\lambda,\overline{m},\overline{M},$ such that for any $t_{0}\in (-1/4,0)$,

(4.25)\begin{align} |u(x,t)-u(0,t_{0})|\leq C\left(|x|+\sqrt[\theta_{2}]{|t-t_{0}|}\right)^{\alpha},\quad\forall (x,t)\in B_{1/2}\times(-1/4,t_{0}]. \end{align}

For $R\in(0,1/2)$, $(y,s)\in Q_{1/R}$, define

\begin{align*} u_{R}(y,s)=u(Ry,R^{\theta_{2}}s),\qquad A_{R}(y)=A(Ry). \end{align*}

Therefore, u R verifies

\begin{align*} |y|^{\theta_{2}}\partial_{s}u^{q}_{R}-\mathrm{div}\left(A_{R}|y|^{2}\nabla u_{R}\right)=0\quad\mathrm{in}\;Q_{1/R}. \end{align*}

By the change of variables, we obtain that this equation keeps uniformly parabolic in $B_{1/2}(\bar{y})\times(-R^{-\theta_{2}},0)$ for any $\bar{y}\in\partial B_{1}$.

For any $(x,t),(\tilde{x},\tilde{t})\in B_{1/2}\times(-1/4,0),$ let $|\tilde{x}|\leq|x|$ without loss of generality. Write $R=|x|$. It then follows from the interior Hölder estimates for uniformly parabolic equations that there exist two constants $0 \lt \beta=\beta(n,p,\theta_{2},\lambda,\overline{m},\overline{M}) \lt 1$ and $0 \lt C=C(n,p,\theta_{2},\lambda,\overline{m},\overline{M})$ such that for any $\bar{y}\in\partial B_{1}$ and $\bar{s}\in(-4^{-1} R^{-\theta_{2}},0)$,

(4.26)\begin{align} |u_{R}(y,s)-u_{R}(\bar{y},\bar{s})|\leq C(|y-\bar{y}|+\sqrt{|s-\bar{s}|})^{\beta}, \end{align}

for any (y, s) satisfying that $|y-\bar{y}|+\sqrt{|s-\bar{s}|} \lt 1/2$.

Observe that for any $(x,t),(\tilde{x},\tilde{t})\in B_{1/2}\times(-1/4,0),$

\begin{align*} |u(x,t)-u(\tilde{x},\tilde{t})|\leq|u(x,t)-u(x,\tilde{t})|+|u(x,\tilde{t})-u(\tilde{x},\tilde{t})|. \end{align*}

On the one hand, if $|t-\tilde{t}|\leq R^{2\theta_{2}}$, then we deduce from Equation (4.26) that

\begin{align*} |u(x,t)-u(x,\tilde{t})|\leq&\left|u_{R}(x/R,t/R^{\theta_{2}})-u_{R}(x/R,\tilde{t}/R^{\theta_{2}})\right|\notag\\ \leq&C|(t-\tilde{t})/R^{\theta_{2}}|^{\beta/2}\leq C|t-\tilde{t}|^{\beta/4}, \end{align*}

while, if $|t-\tilde{t}| \gt R^{2\theta_{2}}$, then we have from Equation (4.25) that

\begin{align*} &|u(x,t)-u(x,\tilde{t})|\notag\\ &\leq|u(x,t)-u(0,t)|+|u(0,t)-u(0,\tilde{t})|+|u(0,\tilde{t})-u(x,\tilde{t})|\notag\\ &\leq C\big(R^{\alpha}+|t-\tilde{t}|^{\frac{\alpha}{\theta_{2}}}\big)\leq C|t-\tilde{t}|^{\frac{\alpha}{2\theta_{2}}}. \end{align*}

On the other hand, if $|x-\tilde{x}|\leq R^{2}$, then it follows from Equation (4.26) that

\begin{align*} |u(x,\tilde{t})-u(\tilde{x},\tilde{t})|=&\left|u_{R}(x/R,\tilde{t}/R^{\theta_{2}})-u_{R}(\tilde{x}/R,\tilde{t}/R^{\theta_{2}}\big)\right|\notag\\ \leq&C|(x-\tilde{x})/R|^{\beta}\leq C|x-\tilde{x}|^{\beta/2}, \end{align*}

while, if $|x-\tilde{x}| \gt R^{2}$, then we see from Equation (4.25) that

\begin{align*} |u(x,\tilde{t})-u(\tilde{x},\tilde{t})|\leq&|u(x,\tilde{t})-u(0,\tilde{t})|+|u(0,\tilde{t})-u(\tilde{x},\tilde{t})|\notag\\ \leq& C\big(R^{\alpha}+|\tilde{x}|^{\alpha}\big)\leq CR^{\alpha}\leq C|x-\tilde{x}|^{\frac{\alpha}{2}}. \end{align*}

Consequently, we complete the proof of Theorem 1.9.

Funding Statement

This work was supported in part by the National Key Research and Development Program of China (No. 2022YFA1005700 and 2020YFA0712903). C. Miao was partially supported by the National Natural Science Foundation of China (No. 12026407 and 12071043). Z. Zhao was partially supported by China Postdoctoral Science Foundation (No. 2021M700358).

References

Akagi, G., Rates of convergence to non-degenerate asymptotic profiles for fast diffusion via energy methods, Arch. Ration. Mech. Anal. 247(no. 2) (2023), .10.1007/s00205-023-01843-2CrossRefGoogle Scholar
Badiale, M. and Tarantello, G., A Sobolev-Hardy inequality with applications to a nonlinear elliptic equation arising in astrophysics, Arch. Ration. Mech. Anal. 163(4) (2002), 259293.10.1007/s002050200201CrossRefGoogle Scholar
Bahouri, H., Chemin, J. -Y. and Gallagher, I., Refined Hardy inequalities, Ann. Sc. Norm. Super. Pisa Cl. Sci. (5) 5(3) (2006), 375391.Google Scholar
Berryman, J. B. and Holland, C. J., Stability of the separable solution for fast diffusion, Arch. Ration. Mech. Anal. 74 (1980), 379388.10.1007/BF00249681CrossRefGoogle Scholar
Bonforte, M. and Figalli, A., Sharp extinction rates for fast diffusion equations on generic bounded domains, Comm. Pure Appl. Math. 74(4) (2021), 744789.10.1002/cpa.21887CrossRefGoogle Scholar
Bonforte, M., Grillo, G. and Vázquez, J. L., Behaviour near extinction for the fast diffusion equation on bounded domains, J. Math. Pures Appl. (9) 97(1) (2012), 138.10.1016/j.matpur.2011.03.002CrossRefGoogle Scholar
Bonforte, M. and Vázquez, J. L., Positivity, local smoothing, and Harnack inequalities for very fast diffusion equations, Adv. Math. 223(2) (2010), 529578.10.1016/j.aim.2009.08.021CrossRefGoogle Scholar
Cabré, X. and Ros-Oton, X., Sobolev and isoperimetric inequalities with monomial weights, J. Differential Equations 255(11) (2013), 43124336.10.1016/j.jde.2013.08.010CrossRefGoogle Scholar
Cabré, X., Ros-Oton, X. and Serra, J., Sharp isoperimetric inequalities via the ABP method, J. Eur. Math. Soc. (JEMS) 18(12) (2016), 29712998.10.4171/JEMS/659CrossRefGoogle Scholar
Caffarelli, L., Kohn, R. and Nirenberg, L., First order interpolation inequalities with weights, Compos. Math. 53 (1984), 259275.Google Scholar
Chanillo, S. and Wheeden, R. L., Weighted Poincaré and Sobolev inequalities and estimates for weighted Peano maximal functions, Amer. J. Math. 107(5) (1985), 11911226.10.2307/2374351CrossRefGoogle Scholar
Chen, Y. Z. and DiBenedetto, E., On the local behavior of solutions of singular parabolic equations, Arch. Ration. Mech. Anal. 103(4) (1988), 319345.Google Scholar
Chiarenza, F. and Serapioni, R. P., A Harnack inequality for degenerate parabolic equations, Comm. Partial Differential Equations 9(8) (1984), 719749.10.1080/03605308408820346CrossRefGoogle Scholar
Chua, S. -K. and Wheeden, R. L., Estimates of best constants for weighted Poincaré inequalities on convex domains, Proc. Lond. Math. Soc. 93(1) (2006), 197226.10.1017/S0024611506015826CrossRefGoogle Scholar
Daskalopoulos, P. and Kenig, C., Degenerate diffusions. Initial value problems and local regularity theory, EMS Tracts in Mathematics, Volume 1 (European Mathematical Society (EMS), Zürich, 2007).10.4171/033CrossRefGoogle Scholar
De Giorgi, E., Sulla differenziabilità e l’analiticità delle estremali degli integrali multipli regolari, Mem. Accad. Sci. Torino. Cl. Sci. Fis. Mat. Nat. (3) 3 (1957), 2543.Google Scholar
DiBenedetto, E., Gianazza, U. and Vespri, V., Harnack’s inequality for degenerate and singular parabolic equations, Springer Monographs in Mathematics (Springer, New York, 2012).10.1007/978-1-4614-1584-8CrossRefGoogle Scholar
DiBenedetto, E. and Kwong, Y. C., Harnack estimates and extinction profile for weak solution of certain singular parabolic equations, Trans. Amer. Math. Soc. 330(2) (1992), 783811.10.1090/S0002-9947-1992-1076615-7CrossRefGoogle Scholar
DiBenedetto, E., Kwong, Y. C. and Vespri, V., Local space-analyticity of solutions of certain singular parabolic equations, Indiana Univ. Math. J. 40(2) (1991), 741765.10.1512/iumj.1991.40.40033CrossRefGoogle Scholar
Dong, H. J., Li, Y. Y. and Yang, Z. L., Optimal gradient estimates of solutions to the insulated conductivity problem in dimension greater than two, arXiv:2110.11313.Google Scholar
Dong, H. J., Li, Y. Y. and Yang, Z. L., Gradient estimates for the insulated conductivity problem: the non-umbilical case, arXiv:2203.10081.Google Scholar
Fabes, E. B., Kenig, C. E. and Serapioni, R. P., The local regularity of solutions of degenerate elliptic equations, Comm. Partial Differential Equations 7(1) (1982), 77116.10.1080/03605308208820218CrossRefGoogle Scholar
Feireisl, E. and Simondon, F., Convergence for semilinear degenerate parabolic equations in several space dimension, J. Dynam. Differential Equations 12 (2000), 647673.10.1023/A:1026467729263CrossRefGoogle Scholar
Grafakos, L., Classical Fourier analysis, Graduate Texts in Mathematics, 3rd edn, Volume 249 (Springer, New York, 2014).10.1007/978-1-4939-1194-3CrossRefGoogle Scholar
Gutiérrez, C. E. and Wheeden, R. L., Harnack’s inequality for degenerate parabolic equations, Comm. Partial Differential Equations 16(4–5) (1991), 745770.10.1080/03605309108820776CrossRefGoogle Scholar
Heinonen, J., Kilpeläinen, T. and Martio, O., Nonlinear potential theory of degenerate elliptic equations, Unabridged Republication of the 1993 Original (Dover Publications, Inc., Mineola, NY, 2006).Google Scholar
Jin, T. L., and Xiong, J. G., Optimal boundary regularity for fast diffusion equations in bounded domains, Amer. J. Math. 145(no. 1) (2023), 151219.10.1353/ajm.2023.0003CrossRefGoogle Scholar
Jin, T. L. and Xiong, J. G., Regularity of solutions to the Dirichlet problem for fast diffusion equations, arXiv:2201.10091.Google Scholar
Kwong, Y. C., Interior and boundary regularity of solutions to a plasma type equation, Proc. Amer. Math. Soc. 104(2) (1988), 472478.10.1090/S0002-9939-1988-0962815-5CrossRefGoogle Scholar
Li, Y. Y., and Yan, X. K., Anisotropic Caffarelli–Kohn–Nirenberg type inequalities, Adv. Math. 419 (2023), .10.1016/j.aim.2023.108958CrossRefGoogle Scholar
Lin, C. S., Interpolation inequalities with weights, Comm. Partial Differential Equations 11(14) (1986), 15151538.Google Scholar
Lindqvist, P., Notes on the stationary p-Laplace equation, SpringerBriefs in Mathematics, (Springer, Cham, 2019).10.1007/978-3-030-14501-9CrossRefGoogle Scholar
Nguyen, H. and Squassina, M., Fractional Caffarelli–Kohn–Nirenberg inequalities, J. Funct. Anal. 274(9) (2018), 26612672.10.1016/j.jfa.2017.07.007CrossRefGoogle Scholar
Nguyen, H. and Squassina, M., On Hardy and Caffarelli–Kohn–Nirenberg inequalities, J. Anal. Math. 139(2) (2019), 773797.10.1007/s11854-025-0077-9CrossRefGoogle Scholar
Sacks, P. E., Continuity of solutions of a singular parabolic equation, Nonlinear Anal. 7(4) (1983), 387409.10.1016/0362-546X(83)90092-5CrossRefGoogle Scholar
Serrin, J., Local behavior of solutions of quasi-linear equations, Acta Math. 111 (1964), 247302.10.1007/BF02391014CrossRefGoogle Scholar
Surnachev, M., A Harnack inequality for weighted degenerate parabolic equations, J. Differential Equations 248(8) (2010), 20922129.10.1016/j.jde.2009.08.021CrossRefGoogle Scholar
Trudinger, N. S., On Harnack type inequalities and their application to quasilinear elliptic equations, Comm. Pure Appl. Math. 20 (1967), 721747.10.1002/cpa.3160200406CrossRefGoogle Scholar
Vázquez, J. L., The porous medium equation: mathematical theory, Oxford Mathematical Monographs (The Clarendon Press, Oxford University Press, Oxford, 2007).Google Scholar