--- title: Nonlocal Double Phase Operator url: https://www.emergentmind.com/topics/nonlocal-double-phase-operator type: topic --- # Nonlocal Double Phase Operator The nonlocal double phase operator is a nonlinear integro-differential operator that combines two fractional interaction mechanisms with different growth exponents and, in general, different differentiability orders. In its prototype form, it is written as \[ \mathcal L u(x)=\operatorname{P.V.}\int_{\mathbb R^n}|u(x)-u(y)|^{p-2}(u(x)-u(y))K_{sp}(x,y)\,dy+\operatorname{P.V.}\int_{\mathbb R^n}a(x,y)|u(x)-u(y)|^{q-2}(u(x)-u(y))K_{tq}(x,y)\,dy, \] or, in the notation of the nonhomogeneous theory, \[ L_a u(x)=2\,\operatorname{P.V.}\int_{\mathbb R^N}|u(x)-u(y)|^{p-2}(u(x)-u(y))K_{s,p}(x,y)\,dy+2\,\operatorname{P.V.}\int_{\mathbb R^N}a(x,y)|u(x)-u(y)|^{q-2}(u(x)-u(y))K_{t,q}(x,y)\,dy. \] The first integral is the \(p\)-phase, the second is the \(q\)-phase weighted by a modulating coefficient \(a(x,y)\ge 0\). This structure is the nonlocal analogue of the classical local double phase operator, and it includes the fractional \(p\)-Laplacian as the special case \(a\equiv 0\) [2505.16461] [2106.04412]. ## 1. Structural form and phase interaction The defining feature of the operator is the coexistence of two distinct fractional elliptic phases. One phase has \(p\)-growth and fractional order \(s\); the other has \(q\)-growth and fractional order \(t\), and is activated or suppressed by the coefficient \(a(x,y)\). The terminology “double phase” refers precisely to this switching between a \(p\)-phase and a \(q\)-phase. In the nonlocal setting the switch occurs at the level of pairwise interactions \((x,y)\), rather than pointwise in \(x\) alone [1901.05864] [2011.11466]. The kernels are typically assumed symmetric and comparable to fractional kernels. Representative hypotheses are \[ \Lambda^{-1}|x-y|^{-n-sp}\le K_{sp}(x,y)\le \Lambda |x-y|^{-n-sp}, \qquad \Lambda^{-1}|x-y|^{-n-tq}\le K_{tq}(x,y)\le \Lambda |x-y|^{-n-tq}, \] or, in the notation of the nonhomogeneous paper, \[ \frac{1}{\Lambda_1|x-y|^{N+sp}}\le K_{s,p}(x,y)\le \frac{\Lambda_1}{|x-y|^{N+sp}}, \qquad \frac{1}{\Lambda_2|x-y|^{N+tq}}\le K_{t,q}(x,y)\le \frac{\Lambda_2}{|x-y|^{N+tq}}. \] Symmetry, continuity away from the diagonal, and, in several papers, translation invariance are also imposed [2505.16461] [2106.04412]. Two parameter regimes recur in the literature. One is the “lower-order perturbation” regime \(tq\le sp\), which appears in weak/viscosity equivalence, self-improvement, parabolic Hölder theory, and Harnack theory [2505.16461] [2106.04412] [2112.04287] [2509.07433]. The other is the harder regime \(s\le t\), where the second phase has the same or higher nonlocal order than the first; in that case Hölder continuity of weak solutions is proved under stronger compatibility conditions such as \[ tq\le sp+\alpha \] when \(a\) is Hölder continuous, and gradient Hölder regularity is obtained under \[ tq\le sp+\min\{1,q-p\} \] when \(a\) is Lipschitz [2108.09623] [2604.22206]. This split is central to the theory. It shows that “nonlocal double phase operator” is not a single fixed model, but a family of mixed-growth fractional operators whose analytical behavior depends sharply on the relation between the two phases. ## 2. Variational, weak, and tail-based formulations The operator is naturally associated with mixed fractional energies of \((p,q)\)-type. In the weak theory, one works in fractional Sobolev spaces together with tail spaces that encode the influence of values outside the localization region. A standard requirement is \[ u\in W^{s,p}(\mathbb R^N)\cap L_{s,p}^{p-1}(\mathbb R^N)\cap L_{a,t,q}^{q-1}(\mathbb R^N), \] or analogous variants such as \(L^{p-1}_{sp}(\mathbb R^n)\) and \(L^{q-1}_{a,tq}(\mathbb R^n)\) [2505.16461] [2106.04412]. For the operator \(L_a\), the weak formulation is expressed through \[ \begin{aligned} H_a(u,v) &=\iint_{\mathbb{R}^N\times\mathbb{R}^N} h_p(u(x)-u(y))(v(x)-v(y))K_{s,p}(x,y)\,dx\,dy \\ &\quad +\iint_{\mathbb{R}^N\times\mathbb{R}^N} a(x,y)\,h_q(u(x)-u(y))(v(x)-v(y))K_{t,q}(x,y)\,dx\,dy, \end{aligned} \] where \(h_\ell(t)=|t|^{\ell-2}t\). A weak supersolution satisfies \[ H_a(u,v)\ge \int_\Omega f(x,u(x),D_s^p u(x),D_{a,t}^q u(x))\,v(x)\,dx \] for every nonnegative \(v\in C_c^\infty(\Omega)\), with equality in the solution case [2505.16461]. The nonhomogeneous theory introduces the nonlocal energy densities \[ D_s^p u(x)=\int_{\mathbb R^N}\frac{|u(x)-u(y)|^p}{|x-y|^{N+sp}}\,dy, \qquad D_{a,t}^q u(x)=\int_{\mathbb R^N}a(x,y)\frac{|u(x)-u(y)|^q}{|x-y|^{N+tq}}\,dy, \] so that the right-hand side may depend on \(x\), \(u\), and both phase-dependent nonlocal energies [2505.16461]. A distinctive analytic feature is the tail. The theory uses quantities such as \(\operatorname{Tail}_{s,p}(v;x,r)\) and \(\operatorname{Tail}_{a,t,q}(v;x,r)\), or their \(sp\)- and \(tq\)-variants, to measure the contribution of \(\mathbb R^n\setminus B_r(x)\) to the operator. This tail dependence is not auxiliary: it enters Caccioppoli estimates, logarithmic lemmas, local boundedness, Hölder estimates, and Harnack inequalities throughout the subject [2505.16461] [2106.04412] [2509.07433]. The functional framework also admits rougher generalizations. One paper replaces the prototype kernels by measurable kernels \(K_{sp},K_{tq}\) satisfying two-sided ellipticity bounds and nonlinearities \(\phi_p,\phi_q\) satisfying \[ |\phi_r(z)|\le \Lambda^+|z|^{r-1},\qquad \phi_r(z)z\ge |z|^r,\qquad r\in\{p,q\}, \] and proves the same self-improving conclusions for this broader class [2011.11466]. ## 3. Viscosity solutions and equivalence with weak solutions The viscosity framework adapts the standard touching-test paradigm to a nonlocal mixed-growth operator. In the homogeneous setting, a viscosity solution is tested with a \(C^2\) function patched to the original solution outside a contact ball; in the nonhomogeneous setting the test inequality must also evaluate the phase-dependent quantities \(D_s^p\) and \(D_{a,t}^q\) on the test function [2106.04412] [2505.16461]. For the nonhomogeneous operator, a viscosity supersolution is lower semicontinuous, satisfies the appropriate tail integrability, and obeys \[ L_a\psi(x_0)\ge f\bigl(x_0,u(x_0),D_s^p\psi(x_0),D_{a,t}^q\psi(x_0)\bigr) \] whenever \(\psi\in C^2(B_r(x_0))\) touches \(u\) from below at \(x_0\), under one of the regularity regimes required in the paper: either \(p>\frac{2}{2-s}\), or \(\nabla\psi(x_0)\neq 0\), or \(p\le \frac{2}{2-s}\) with \(x_0\) an isolated critical point and \(\psi\in C_\beta^2\) for \(\beta>\frac{sp}{p-1}\) [2505.16461]. A foundational pointwise compatibility statement was proved for bounded viscosity solutions under the natural energy finiteness assumption \[ \int_{\mathbb R^n}\int_{\mathbb R^n}\left(\frac{|u(x)-u(x+y)|^p}{|y|^{n+sp}}+a(x,y)\frac{|u(x)-u(x+y)|^q}{|y|^{n+tq}}\right)\,dx\,dy<\infty. \] Under the structural assumptions of that paper, if \(\mathcal L u\le C\) in viscosity sense and a \(C^2\) function touches from above, then \(\mathcal L u(x_0)\) exists pointwise and satisfies the same inequality at the contact point [1901.05864]. The first systematic weak/viscosity bridge for the nonlocal double phase equation was established in the homogeneous setting. Bounded weak supersolutions were shown to be viscosity supersolutions, using local Hölder continuity, a comparison principle, continuity of the operator on smooth functions, and the fact that affine functions are annihilated under the symmetry and translation-invariance assumptions [2106.04412]. The nonhomogeneous theory then established a two-way equivalence. If \(u\) is a continuous weak supersolution, a suitable comparison principle holds, and \(f\) is continuous in \(x\) and \(t\) and Lipschitz in the nonlocal derivative variables, then \(u\) is a viscosity supersolution. Conversely, if \(u\) is a bounded viscosity supersolution and \(f\) is uniformly continuous in \(\Omega\times\mathbb R^3\), Lipschitz in the nonlocal derivative variables, nonincreasing in \(t\), and satisfies the stated growth condition, then \(u\) is a weak supersolution. The second implication uses infimal convolution \(u_\epsilon\), an approximating source \(f_\epsilon\), convergence of the nonlocal energy terms, and passage to the limit in the weak formulation [2505.16461]. A persistent misconception is that weak and viscosity notions are automatically interchangeable for this class of operators. The published theory shows instead that equivalence is a theorem with nontrivial hypotheses, not a formal identity. In particular, one paper explicitly noted that the reverse implication viscosity \(\Rightarrow\) weak remained open in its setting, and a later paper resolved it for a nonhomogeneous class under additional structural conditions on \(f\) [2106.04412] [2505.16461]. ## 4. Regularity theory The regularity theory begins with Hölder continuity. The first Hölder result for bounded viscosity solutions proved that if \[ \mathcal L u=f \quad \text{in } B_2, \qquad f\in L^\infty(B_2), \] then \(u\in C^{0,\gamma}(B_1)\) for some \(\gamma=\gamma(\text{data})\in(0,1)\), assuming only boundedness of the modulating coefficient \[ 0\le a(x,y)\le M \] and suitable exponent restrictions. The proof follows a nonlocal oscillation-decay scheme in the style of Krylov–Safonov rather than a De Giorgi–Nash–Moser energy argument [1901.05864]. The weak-solution theory in the \(tq\le sp\) regime was then developed by De Giorgi–Nash–Moser methods. Under the structural assumptions of the homogeneous equation, bounded weak solutions are locally Hölder continuous, and the oscillation satisfies \[ \operatorname{osc}_{B_\rho(x_0)}u \le C\left(\frac{\rho}{r}\right)^\alpha \|u\|_{L^\infty(\mathbb R^n)} \] for \(B_{2r}(x_0)\Subset\Omega\) and \(0<\rho\le r\) [2106.04412]. A different phenomenon appears in the higher-order regime \(s\le t\). For weak solutions in that setting, local boundedness holds under the growth restriction \[ p\le q\le \frac{np}{n-sp}\quad \text{when } sps_2\) [2112.04287]. The current highest regularity result concerns gradients. Interior \(C^{1,\alpha}\) regularity of viscosity solutions has been proved for the degenerate case, where \(a\) may vanish, assuming \(a\) is symmetric, translation invariant, bounded, and Lipschitz, and that \[ tq\le sp+\min\{1,q-p\}. \] That work explicitly states that it resolves the higher regularity issue raised by De Filippis–Palatucci [2604.22206]. A common oversimplification is to say that boundedness of \(a\) is either always sufficient or always insufficient for regularity. The published results show a more nuanced picture. Bounded measurable \(a\) suffices for Hölder continuity of bounded viscosity solutions in one regime [1901.05864], whereas Hölder continuity of \(a\) is required for Hölder regularity of weak solutions in the \(s\le t\) regime [2108.09623], and Lipschitz continuity of \(a\) is used for \(C^{1,\alpha}\) regularity [2604.22206]. This suggests that the threshold on the coefficient depends both on the phase interaction and on the level of regularity being sought. ## 5. Self-improvement, De Giorgi classes, and Harnack theory Beyond baseline Hölder continuity, the operator exhibits self-improving regularity. For bounded weak solutions under \[ p\le q,\qquad t\le s,\qquad \frac1{p'}\le \frac{tq}{sp}\le 1,\qquad sp0\). The proof uses a dual-pair reformulation \((U,\nu)\), a reverse Hölder inequality on diagonal sets, and a fractional Gehring lemma adapted from Kuusi–Mingione–Sire [2011.11466]. A related self-improving theory has been proved for nonlocal double phase equations with VMO coefficients and solution-dependent principal coefficient \(a(x,y,u(x),u(y))\). Under \(2\le p\le q\le ps/t\), weak solutions satisfy a local higher differentiability and integrability property, \[ u\in W^{s+\delta,\;p(1+\delta)}_{\mathrm{loc}}(\Omega), \] and under stronger assumptions they satisfy \[ u\in C^{0,\alpha}_{\mathrm{loc}}(\Omega)\qquad \text{for every }\alpha\in(0,\Theta), \] where \[ \Theta=\min\left\{\frac{ps-\frac N\gamma}{p-1},\ \frac{qt}{q-p},\ 1\right\}. \] The mechanism there is perturbative freezing of VMO coefficients together with Campanato iteration [2303.07749]. The nonhomogeneous theory also develops a double phase De Giorgi class \(\mathrm{DG}_a^\pm\) for the simpler equation \[ L_a u=f(x,u), \qquad |f(x,t)|\le c_1+c_2|t|^{l-1},\quad 1