The $A_{\infty}$ condition, $\varepsilon$ -approximators, and Varopoulos extensions in uniform domains

S. Bortz, B. Poggi, O. Tapiola, and X. Tolsa Simon Bortz, Department of Mathematics
University of Alabama
Tuscaloosa, AL, 35487, USA sbortz@ua.edu Bruno Poggi, Departament de Matemàtiques
Universitat Autònoma de Barcelona
Barcelona, Catalonia. bgpoggi.math@gmail.com Olli Tapiola, Departament de Matemàtiques
Universitat Autònoma de Barcelona
Barcelona, Catalonia. olli.m.tapiola@gmail.com Xavier Tolsa, ICREA
Barcelona
Departament de Matemàtiques
Universitat Autònoma de Barcelona and Centre de Recerca Matemàtica
Barcelona
Catalonia. xtolsa@mat.uab.cat

Abstract.

Suppose that $\Omega\subset\mathbb{R}^{n+1}$ , $n\geq 1$ , is a uniform domain with $n$ -Ahlfors regular boundary and $L$ is a (not necessarily symmetric) divergence form elliptic, real, bounded operator in $\Omega$ . We show that the corresponding elliptic measure $\omega_{L}$ is quantitatively absolutely continuous with respect to surface measure of $\partial\Omega$ in the sense that $\omega_{L}\in A_{\infty}(\sigma)$ if and only if any bounded solution $u$ to $Lu=0$ in $\Omega$ is $\varepsilon$ -approximable for any $\varepsilon\in(0,1)$ . By $\varepsilon$ -approximability of $u$ we mean that there exists a function $\Phi=\Phi^{\varepsilon}$ such that $\|u-\Phi\|_{L^{\infty}(\Omega)}\leq\varepsilon\|u\|_{L^{\infty}(\Omega)}$ and the measure $\widetilde{\mu}_{\Phi}$ with $d\widetilde{\mu}=|\nabla\Phi(Y)|\,dY$ is a Carleson measure with $L^{\infty}$ control over the Carleson norm.

As a consequence of this approximability result, we show that boundary $\operatorname{BMO}$ functions with compact support can have Varopoulos-type extensions even in some sets with unrectifiable boundaries, that is, smooth extensions that converge non-tangentially back to the original data and that satisfy $L^{1}$ -type Carleson measure estimates with $\operatorname{BMO}$ control over the Carleson norm. Our result complements the recent work of Hofmann and the third named author who showed the existence of these types of extensions in the presence of a quantitative rectifiability hypothesis.

S.B. was supported by the Simons foundation grant “Travel support for Mathematicians” (grant number 959861). B.P., O.T. and X.T. were supported by the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement 101018680). X.T. is also partially supported by MICINN (Spain) under the grant PID2020-114167GB-I00, the María de Maeztu Program for units of excellence (Spain) (CEX2020-001084-M), and 2021-SGR-00071 (Catalonia).

1. Introduction

Carleson measures (see Definition 2.7) and their connections to geometry, harmonic analysis and partial differential equations (PDE) have been studied actively over the previous decades; see [Mat21] for a recent survey related to many key developments and other relevant research. Carleson-type estimates are particularly powerful for measures $\mu_{F}$ such that $d\mu_{F}=|\nabla F(Y)|^{2}\operatorname{dist}(Y,\partial\Omega)\,dY$ and $F$ is a solution to an elliptic PDE in the set $\Omega$ : they can be used to, for example, characterize functions of bounded mean oscillation (BMO) [FS72], quantitative boundary geometry [HMM16, GMT18] and absolute continuity properties of harmonic measure [HL18]. As it is discussed in [Gar07, Chapters VI and VIII], Carleson-type estimates for measures $\widetilde{\mu}_{F}$ such that $d\widetilde{\mu}_{F}=|\nabla F(Y)|\,dY$ would be very powerful but they fail even for harmonic functions in the unit disk. To circumvent this problem, Varopoulos [Var78] introduced a way to approximate harmonic functions in $L^{\infty}$ sense by other functions satisfying these estimates. This $\varepsilon$ -approximability theory has been studied from many points of view in the past years [Gar07, Dah80, KKPT00, HKMP15, HMM16, HR18, GMT18, AGMT22, HT20, BT19, BH20, Gar22].

The initial motivation behind $\varepsilon$ -approximability theory was to prove an extension theorem for BMO functions inspired by Carleson’s Corona Theorem [Car62]. Varopoulos [Var77, Var78] showed that any compactly supported BMO function $f$ in ${\mathbb{R}}^{n}$ has a smooth extension $V_{f}$ to the upper half-space such that the extension converges non-tangentially (see Definition 2.8) back to $f$ and $|\nabla V_{f}(Y)|\,dY$ defines a Carleson measure. Inspired by the power of these estimates, Hofmann and the third named author [HT21] recently showed that uniform rectifiability (see Definition 2.6) of the boundary is enough to guarantee the existence of these Varopoulos-type extensions. Since Carleson-type estimates for harmonic functions can be used to characterize uniform rectifiability [HMM16, GMT18] or even stronger geometric properties [AHMMT20], it is natural to ask if the existence of Varopoulos-type extensions (which satisfy better Carleson-type estimates than harmonic functions) characterizes some quantitative geometric properties for the boundary.

In this paper, we show that the existence of Varopoulos-type extensions does not characterize uniform rectifiability but they can exist even in sets with unrectifiable boundaries. Our main result is the following $\varepsilon$ -approximability result which can be used to build Varopoulos-type extensions. Throughout, set $\sigma\coloneqq\mathcal{H}^{n}|_{\partial\Omega}$ .

Theorem 1.1.

Let $\Omega\subset\mathbb{R}^{n+1}$ , $n\geq 1$ , be a uniform domain (see Definition 2.3) with $n$ -Ahlfors regular boundary (see Definition 2.4). Let $L=-\operatorname{div}A\nabla$ be a real, not necessarily symmetric, bounded elliptic operator in $\Omega$ such that the corresponding elliptic measure $\omega_{L}$ satisfies $\omega_{L}\in A_{\infty}(\sigma)$ (see Definition 2.31), and let $\varepsilon\in(0,1)$ . Then any solution $u\in W^{1,2}_{\operatorname{loc}}(\Omega)\cap L^{\infty}(\Omega)$ to $Lu=0$ in $\Omega$ is $\varepsilon$ -approximable: there exists a constant $C_{\varepsilon}$ and a function $\Phi=\Phi^{\varepsilon}\in C^{\infty}(\Omega)$ such that

i)

$\|u-\Phi\|_{L^{\infty}(\Omega)}\leq\varepsilon\|u\|_{L^{\infty}(\Omega)}$ ,

ii)

$\Phi$ satisfies a quantitative $L^{1}$ -type Carleson measure estimate

\displaystyle\sup_{x\in\partial\Omega,r>0}\frac{1}{r^{n}}\iint_{B(x,r)\cap\Omega}|\nabla\Phi(Y)|\,dY\leq C_{\varepsilon}\|u\|_{L^{\infty}(\Omega)},

iii)

$|\nabla\Phi(X)|\leq\tfrac{C_{\varepsilon}\|u\|_{L^{\infty}(\Omega)}}{\delta(X)}$ for every $X\in\Omega$ ,
iv)

if $|X-Y|\ll\operatorname{dist}(X,\partial\Omega)$ , then $|\Phi(X)-\Phi(Y)|\leq\tfrac{C_{\varepsilon}\|u\|_{L^{\infty}(\Omega)}}{\delta(X)}|X-Y|$ ,
v)

there exists a function $\varphi\in L^{\infty}(\partial\Omega)$ such that

$\displaystyle\lim_{Y\to x,\,\operatorname{n.t.}}\Phi(Y)=\varphi(x)\text{ for }\sigma\text{-a.e. }x\in\partial\Omega.$

The notation $\lim_{Y\to x,\,\operatorname{n.t.}}$ means non-tangential convergence (see Definition 2.8) and $\delta(\cdot)\coloneqq\operatorname{dist}(\cdot,\partial\Omega)$ . Here, $C_{\varepsilon}$ depends on $\varepsilon$ , the structural constants related to $\Omega$ and $\partial\Omega$ , ellipticity and the $\omega_{L}\in A_{\infty}(\sigma)$ constants.

The proof of Theorem 1.1 borrows some ideas from the proof of [HMM16, Theorem 1.3] (which is an adaptation of the classical construction in [Gar07, Chapter VIII, Theorem 6.1]), but very quickly our argument must differ significantly. In [HMM16], the authors constructed the approximators for harmonic functions in the presence of a quantitative rectifiability hypothesis. This allowed them to construct approximating chord-arc domains which, in turn, allowed them to use an “ $N\lesssim S$ ” estimate for harmonic functions. This was the most delicate part of their argument and our main challenges are strongly related to overcoming the fact that we cannot use the same tools due to our geometry (our boundary may be purely unrectifiable) and our operator $L$ (no control on its structure or smoothness).

We also obtain the following converse to Theorem 1.1.

Theorem 1.2.

Let $\Omega\subset\mathbb{R}^{n+1}$ , $n\geq 1$ , be a uniform domain (see Definition 2.3) with $n$ -Ahlfors regular boundary (see Definition 2.4), and let $L=-\operatorname{div}A\nabla$ be a real, not necessarily symmetric, bounded elliptic operator in $\Omega$ . Suppose also that every solution $u\in W^{1,2}_{\operatorname{loc}}(\Omega)\cap L^{\infty}(\Omega)$ to $Lu=0$ is $\varepsilon$ -approximable for every $\varepsilon\in(0,1)$ in the sense of Theorem 1.1. Then $\omega_{L}\in A_{\infty}(\sigma)$ (see Definition 2.31).

In particular, by combining Theorem 1.1 and Theorem 1.2 with [CHMT20, Theorem 1.1], we get a new characterization of the $A_{\infty}$ property of elliptic measure on uniform domains.

Corollary 1.3.

(a)

$\omega_{L}\in A_{\infty}(\sigma)$ ,
(b)

every solution $u\in W^{1,2}_{\operatorname{loc}}(\Omega)\cap L^{\infty}(\Omega)$ to $Lu=0$ in $\Omega$ is $\varepsilon$ -approximable for any $\varepsilon\in(0,1)$ in the sense of Theorem 1.1,

(c)

every solution $u\in W^{1,2}_{\operatorname{loc}}(\Omega)\cap L^{\infty}(\Omega)$ to $Lu=0$ in $\Omega$ satisfies an $L^{2}$ -type Carleson measure estimate with $L^{\infty}$ control over the Carleson norm: there exists $C\geq 1$ such that

\displaystyle\sup_{\begin{subarray}{c}x\in\partial\Omega\\ r\in(0,\operatorname{diam}(\partial\Omega))\end{subarray}}\frac{1}{r^{n}}\iint_{B(x,r)\cap\Omega}|\nabla u(X)|^{2}\operatorname{dist}(X,\partial\Omega)\,dX\leq C\|u\|_{L^{\infty}{\Omega}}^{2}.

Only the implications “(a) $\implies$ (b)” and “(b) $\implies$ (a)” in Corollary 1.3 are new. The equivalence “(a) $\iff$ (c)” was already shown in [CHMT20] for $n\geq 2$ , while the case $n=1$ follows as a particular case of a more general result in [FP22], with a similar method of proof.

In the setting of a uniform domain with Ahlfors regular boundary, the conditions (a), (b) and (c) in Corollary 1.3 are known to be equivalent with uniform rectifiability of $\partial\Omega$ for the special case $L=-\Delta$ [HM14, HMU14, HMM16, GMT18], or for $L=-\operatorname{div}A\nabla$ where $A$ is a locally Lipschitz symmetric matrix such that $|\nabla A|$ satisfies an $L^{1}$ -type Carleson measure condition [CHMT20, HMT17, AGMT22]. Even for the aforementioned operators with nice structure, all the available proofs in the literature of the implications “(a) $\implies$ (b)” or “(c) $\implies$ (b)” in Corollary 1.3 rely on uniform rectifiability techniques in a decisive fashion, and do not extend to domains with rougher boundaries. The main novelty of this manuscript is that we succeed in proving the implication “(a) $\implies$ (b)” without appealing to uniform rectifiability theory. This allows us to establish the equivalence “(a) $\iff$ (b)” for arbitrary elliptic operators in settings that are beyond chord-arc domains (see Definition 2.5). For further characterizations of the $\omega_{L}\in A_{\infty}(\sigma)$ property for arbitrary real divergence form elliptic operators, see [CDMT22] and [MPT].

Let us see an immediate corollary of our new characterization of the $\omega_{L}\in A_{\infty}(\sigma)$ property. We say that $L$ is a Dahlberg–Kenig–Pipher operator if $A\in\operatorname{Lip}_{\operatorname{loc}}(\Omega)$ with $|\nabla A|\operatorname{dist}(\cdot,\partial\Omega)\in L^{\infty}(\Omega)$ , and the measure $\mu_{A}$ such that $d\mu_{A}=|\nabla A(X)|^{2}\operatorname{dist}(X,\partial\Omega)\,dX$ is a Carleson measure. For these operators, the conditions (a) and (c) in Corollary 1.3 are equivalent with uniform rectifiability of $\partial\Omega$ when $\Omega$ is a uniform domain with Ahlfors regular boundary [HMMTZ21]. Combining the main result of [HMMTZ21] with Corollary 1.3 gives us a new result for Dahlberg–Kenig–Pipher operators:

Corollary 1.4.

Let $\Omega\subset\mathbb{R}^{n+1}$ , $n\geq 2$ , be a uniform domain (see Definition 2.3) with $n$ -Ahlfors regular boundary (see Definition 2.4), and let $L=-\operatorname{div}A\nabla$ be a not necessarily symmetric Dahlberg–Kenig–Pipher operator in $\Omega$ . The following are equivalent:

(a)

$\partial\Omega$ is uniformly rectifiable (see Definition 2.6).
(b)

every solution $u\in W^{1,2}_{\operatorname{loc}}(\Omega)\cap L^{\infty}(\Omega)$ to $Lu=0$ in $\Omega$ is $\varepsilon$ -approximable for any $\varepsilon\in(0,1)$ in the sense of Theorem 1.1.

As a consequence of Theorem 1.1 and the techniques in [HT21], we get the following generalization of the Varopoulos extension theorem [Var77, Var78]:

Theorem 1.5.

Let $\Omega\subset\mathbb{R}^{n+1}$ , $n\geq 1$ , be a uniform domain (see Definition 2.3) with Ahflors regular boundary (see Definition 2.4). If there exists a divergence form elliptic operator $L=-\mathop{\operatorname{div}}\nolimits A\nabla$ such that the corresponding elliptic measure $\omega_{L}$ satisfies $\omega_{L}\in A_{\infty}(\sigma)$ (see Definition 2.31), then every $f\in\operatorname{BMO}_{\operatorname{c}}(\partial\Omega)$ has a Varopoulos extension in $\Omega$ . That is, there exists $F$ with the following properties:

(1)

$F\in C^{\infty}(\Omega)$ and $|\nabla F(X)|\lesssim\tfrac{\|f\|_{\operatorname{BMO}(\partial\Omega)}}{\delta(X)}$ , for all $X\in\Omega$ ,
(2)

$\lim_{Y\to x,\,\operatorname{n.t.}}F(Y)=f(x)$ for $\sigma$ -a.e. $x\in\partial\Omega$ , and

(3)

$|\nabla F(Y)|$ is the density of a Carleson measure in the sense that

\sup_{r>0,x\in\partial\Omega}\frac{1}{r^{n}}\iint_{B(x,r)\cap\Omega}|\nabla F(Y)|\,dY\leq C\|f\|_{\operatorname{BMO}(\partial\Omega)}.

The notation $\lim_{Y\to x,\,\operatorname{n.t.}}$ means non-tangential limit and $\delta(\cdot)\coloneqq\operatorname{dist}(\cdot,\partial\Omega)$ . Here, $C$ depends on the structural constants related to $\Omega$ and $\partial\Omega$ , and the $\omega_{L}\in A_{\infty}(\sigma)$ constants.

Theorem 1.5 is not (and is not meant to be) a generalization of the main result in [HT21] where an extension theorem of this type was proven in the presence of a quantitative rectifiability hypothesis for the boundary. The novelty of Theorem 1.5 is that its assumptions hold for some sets with very rough boundaries. In particular, recently David and Mayboroda [DM21] showed that the key hypothesis of Theorem 1.5 holds for the exterior of the $4$ -corner Cantor set:

Theorem 1.6 ([DM21, Section 4]).

Let $\Omega$ be the complement of the $4$ -corner Cantor set in $\mathbb{R}^{2}$ (see Section 7). There exists a divergence form elliptic operator $L=-\mathop{\operatorname{div}}\nolimits A\nabla$ in $\Omega$ such that $\omega_{L}\in A_{\infty}(\sigma)$ (see Definition 2.31). More precisely, one can take the matrix $A$ to be diagonal and equal to the identity outside a ball of radius $1$ concentric with the Cantor set and so that the $L$ -elliptic measure with pole at $\infty$ equals $\sigma/\sigma(\partial\Omega)$ .

The remarkable thing about Theorem 1.6 is that since the complement of the $4$ -corner Cantor set is an unbounded uniform domain with unrectifiable $1$ -Ahlfors regular boundary, we know that harmonic measure for this set cannot satisfy the $A_{\infty}(\sigma)$ condition (see, for example, [AHMMT20]). Thus, constructing operators like this is highly non-trivial. In Section 7, we take the David–Mayboroda example and use it to build an example of a similar operator in $\mathbb{R}^{3}$ (see Proposition 7.6).

By combining Theorems 1.5 and 1.6 and Proposition 7.6, we get the following:

Corollary 1.7.

There exist uniform domains $\Omega$ with unrectifiable Ahlfors regular boundaries $\partial\Omega$ in $\mathbb{R}^{2}$ and $\mathbb{R}^{3}$ such that every function $f\in\operatorname{BMO}(\partial\Omega)$ with compact support has a Varopoulos extension in $\Omega$ . In particular, the existence of Varopoulos extensions does not imply rectifiability for the boundary.

By Corollary 1.7, $L^{1}$ -type Carleson measure estimates are simultaneously too strong and too weak from the point of view of the David–Semmes theory: harmonic functions fail these estimates even in the unit disk but the existence of non-harmonic extensions that satisfy these estimates does not imply even qualitative rectifiability for the boundary.

Finally, we want to mention an upcoming paper that is closely related to our results. As we were finishing this work, we were informed¹¹1Personal communication. about an upcoming manuscript by M. Mourgoglou and T. Zacharopoulos, where the authors construct Varopoulos extensions in corkscrew domains with Ahlfors regular boundaries satisfying a mild quantitative connectivity hypothesis by different methods.

The paper is organized as follows. In Section 2, we discuss basic definitions and consider some key tools from dyadic analysis and elliptic PDE theory. In Section 3, we prove some important preliminary estimates for Theorem 1.1, and in Section 4, we prove Theorem 1.1. In Section 5, we prove Theorem 1.2 (and hence, Corollary 1.3). Finally, in Section 6 we sketch the proof of Theorem 1.5 and in Section 7 we construct a David–Mayboroda-type example in $\mathbb{R}^{3}$ (which completes the proof of Corollary 1.7).

2. Preliminaries

Throughout, we let $\Omega\subset\mathbb{R}^{n+1}$ be an open set with $n\geq 1$ . We say that $\Omega$ is a domain if it is also connected.

Usually, we use capital letters $X,Y,Z$ , and so on to denote points in $\Omega$ , and lowercase letters $x,y,z$ , and so on to denote points in $\partial\Omega$ . For $X\in\mathbb{R}^{n+1}$ and $r>0$ , we let $B(X,r)$ be the Euclidean open ball of radius $r$ centered at $X$ . The letters $c$ and $C$ and their obvious variations denote constants that depend only on dimension, $n$ -Ahlfors regularity constant (see Definition 2.4), corkscrew constant (see Definition 2.1), Harnack chain constants, ellipticity constants (see Section 2.5), and so on. We call these kinds of constants structural constants. We write $a\lesssim b$ if $a\leq Cb$ for a structural constant $C$ and $a\approx b$ if $C_{1}b\leq a\leq C_{2}b$ for structural constants $C_{1}$ and $C_{2}$ .

2.1. Uniform domains, chord-arc domains, Ahlfors regularity and uniform rectifiability

Definition 2.1 (Corkscrew condition).

We say that a domain $\Omega\subset\mathbb{R}^{n+1}$ satisfies the corkscrew condition if there exists a constant $\gamma>0$ such that for every $x\in\partial\Omega$ and $r\in(0,\operatorname{diam}(\Omega))$ there exists $Y_{x,r}$ such that

B(Y_{x,r},\gamma r)\subset B(x,r)\cap\Omega.

We call $Y_{x,r}$ a corkscrew point relative to $x$ at scale $r$ .

Definition 2.2 (Harnack chain Condition).

We say a domain $\Omega\subset\mathbb{R}^{n+1}$ satisfies the Harnack chain condition if there exists a uniform constant $C$ such that for every $\rho>0$ and $\Lambda\geq 0$ and $X,X^{\prime}\in\Omega$ with $\operatorname{dist}(X,\partial\Omega),\operatorname{dist}(X^{\prime},\partial\Omega)\geq\rho$ and $|X-X^{\prime}|\leq\Lambda\rho$ there exists a chain of open balls $B_{1},\dots,B_{J}$ with $J\leq N(\Lambda)$ with $X\in B_{1}$ , $X^{\prime}\in B_{J}$ , $B_{j}\cap B_{j+1}\neq\mbox{{\O}}$ and $C^{-1}\operatorname{diam}(B_{j})\leq\operatorname{dist}(B_{j},\partial\Omega)\leq C\operatorname{diam}(B_{j})$ .

Definition 2.3 (Uniform domain).

We say that a domain $\Omega\subset\mathbb{R}^{n+1}$ is uniform if it satisfies the corkscrew and Harnack chain conditions.

Definition 2.4 (Ahlfors regularity).

We say $\Sigma\subset\mathbb{R}^{n+1}$ is $n$ -Ahlfors regular (or simply Ahlfors regular) if there exists $C$ such that

C^{-1}r^{n}\leq\mathcal{H}^{n}(B(x,r)\cap\Sigma)\leq Cr^{n},\quad\text{for each }x\in r\in(0,\operatorname{diam}(\Sigma)).

Here and below $\mathcal{H}^{n}$ denotes the $n$ -dimensional Hausdorff measure.

Definition 2.5 (Chord-arc domain).

We say that a domain $\Omega\subset\mathbb{R}^{n+1}$ is a chord-arc domain if $\Omega$ satisfies the Harnack chain condition, both $\Omega$ and $\text{int}\,\Omega^{\operatorname{c}}$ satisfy the corkscrew condition, and the boundary $\partial\Omega$ is $n$ -Ahlfors regular.

Definition 2.6 (Uniform rectifiability).

Following [DS91], we say that an $n$ -Ahlfors regular set $E\subset\mathbb{R}^{n+1}$ is uniformly rectifiable if it contains “big pieces of Lipschitz images” of $\mathbb{R}^{n}$ : there exist constants $\theta,M>0$ such that for every $x\in E$ and $r\in(0,\operatorname{diam}(E))$ there is a Lipschitz mapping $\rho=\rho_{x,r}\colon\mathbb{R}^{n}\to\mathbb{R}^{n+1}$ , with Lipschitz norm no larger that $M$ , such that

\displaystyle\mathcal{H}^{n}\big{(}E\cap B(x,r)\cap\rho(\{y\in\mathbb{R}^{n}\colon|y|<r\})\big{)}\geq\theta r^{n}.

2.2. Carleson measures, non-tangential convergence, BMO and local BV

Given a domain $\Omega\subset\mathbb{R}^{n+1}$ , we set $\sigma\coloneqq\mathcal{H}^{n}|_{\partial\Omega}$ and $\delta(X)\coloneqq\operatorname{dist}(X,\partial\Omega)$ for $X\in\Omega$ .

Definition 2.7 (Carleson measures).

We say that a Borel measure $\mu$ in $\Omega$ is a Carleson measure (with respect to $\partial\Omega$ ) if we have

\displaystyle C_{\mu}\coloneqq\sup_{x\in\partial\Omega,r>0}\frac{\mu(B(x,r)\cap\Omega)}{r^{n}}<\infty.

We call $C_{\mu}$ the Carleson norm of $\mu$ .

Definition 2.8 (Cones and non-tangential convergence).

Suppose that $m>1$ . For every $x\in\partial\Omega$ , the cone of $m$ -aperture at $x$ is the set

(2.9)

\displaystyle\widetilde{\Gamma}(x)\coloneqq\widetilde{\Gamma}^{m}(x)\coloneqq\{Z\in\Omega\colon\operatorname{dist}(Z,x)<m\operatorname{dist}(Z,\partial\Omega)\}.

Let $G$ be a function defined in $\Omega$ and $g$ be a function defined on $\partial\Omega$ . We say that $G$ converges non-tangentially to $g$ at $x\in\partial\Omega$ if there exists $m>1$ such that we have $\lim_{k\to\infty}G(Y_{k})=g(x)$ for every sequence $(Y_{k})$ in $\widetilde{\Gamma}^{m}(x)$ such that $\lim_{k\to\infty}Y_{k}=x$ . We denote this by $\lim_{Y\to x,\,\operatorname{n.t.}}G(Y)=g(x)$ .

Definition 2.10 (Non-tangential maximal operator).

We denote the non-tangential maximal operator by $N_{*}$ , that is, for a function $u\in L^{\infty}(\Omega)$ , the function $N_{*}u\colon\partial\Omega\to\mathbb{R}$ is defined as

\displaystyle N_{*}u(x)=\sup_{X\in\widetilde{\Gamma}(x)}|u(X)|.

We call $N_{*}u$ the non-tangential maximal function of $u$ .

Remark 2.11.

The aperture constant $m$ in Definition 2.8 does not play a big role in this paper and therefore we do not analyze it in detail for our results; we simply have that the results hold for some uniform aperture constant. Naturally, the aperture constant affects the values of the non-tangential maximal function but since the $L^{p}$ norms of non-tangential maximal functions that are defined with different aperture constants are comparable (with the comparability constant depending on the aperture constants) (see [FS72, Lemma 1] and [HT20, Lemma 1.10]), this is not important for us. In most computations, it is convenient to use dyadic cones (see (2.17)) instead of cones of the previous type.

Definition 2.12 (BMO).

The space $\operatorname{BMO}(\partial\Omega)$ (bounded mean oscillation) consists of $f\in L^{1}_{\operatorname{loc}}(\partial\Omega)$ with

\displaystyle\|f\|_{\operatorname{BMO}}\coloneqq\sup_{\Delta}\fint_{\Delta}|f(y)-\langle f\rangle_{\Delta}|\,d\sigma(y)<\infty,

where the supremum is taken over all surface balls $\Delta=\Delta(x,r)\coloneqq B(x,r)\cap\partial\Omega$ . We denote $f\in\operatorname{BMO}_{\operatorname{c}}(\partial\Omega)$ if $f$ is a $\operatorname{BMO}$ function with compact support.

Definition 2.13 (Local BV).

We say that locally integrable function $f$ has locally bounded variation in $\Omega$ (denote $f\in\operatorname{BV}_{\operatorname{loc}}(\Omega)$ ) if for any open relatively compact set $\Omega^{\prime}\subset\Omega$ the total variation over $\Omega^{\prime}$ is finite:

\displaystyle\iint_{\Omega^{\prime}}|\nabla f(Y)|\,dY\coloneqq\sup_{\begin{subarray}{c}\overrightarrow{\Psi}\in C_{0}^{1}(\Omega^{\prime})\\ \|\overrightarrow{\Psi}\|_{L^{\infty}(\Omega^{\prime})\leq 1}\end{subarray}}\iint_{\Omega^{\prime}}f(Y)\,\text{div}\overrightarrow{\Psi}(Y)\,dY<\infty,

where $C_{0}^{1}(\Omega^{\prime})$ is the class of compactly supported continuously differentiable vector fields in $\Omega^{\prime}$ .

2.3. Dyadic cubes, Whitney regions and approximating domains

An $n$ -Ahlfors regular set $E\subset\mathbb{R}^{n+1}$ equipped with the Euclidean distance and surface measure can be viewed as a space of homogeneous type of Coifman and Weiss [CW71], with ambient dimension $n+1$ . All such sets can be decomposed dyadically in the following sense:

Lemma 2.14 ([Chr90, DS91, HK12]).

Assume that $E\subset\mathbb{R}^{n+1}$ is $n$ -Ahlfors regular. Then $E$ admits a dyadic decomposition in the sense that there exist constants $a_{1}\geq a_{0}>0$ such that for each $k\in\mathbb{Z}$ there exists a collection of Borel sets, $\mathbb{D}_{k}$ , which we will call (dyadic) cubes, such that

\mathbb{D}_{k}\coloneqq\{Q_{j}^{k}\subset E\colon j\in\mathfrak{I}_{k}\},

where $\mathfrak{I}_{k}$ denotes a countable index set depending on $k$ , satisfying

(i)

for each fixed $k\in\mathbb{Z}$ , the sets $Q_{j}^{k}$ are disjoint and $E=\cup_{j}Q_{j}^{k}\,\,$ ,
(ii)

if $m\geq k$ then either $Q_{i}^{m}\subset Q_{j}^{k}$ or $Q_{i}^{m}\cap Q_{j}^{k}=\mbox{{\O}}$ ,
(iii)

for each $k\in\mathbb{Z}$ , $j\in\mathfrak{I}_{k}$ and $m<k$ , there is a unique $i\in\mathfrak{I}_{m}$ such that $Q_{j}^{k}\subset Q_{i}^{m}$ ,
(iv)

$\operatorname{diam}(Q_{j}^{k})\leq a_{1}2^{-k}$ ,
(v)

for each $Q_{j}^{k}$ , there exists a point $z_{j}^{k}\in Q_{j}^{k}$ such that $E\cap B(z^{k}_{j},a_{0}2^{-k})\subset Q_{j}^{k}\subset E\cap B(z^{k}_{j},a_{1}2^{-k})$ .

We denote by $\mathbb{D}=\mathbb{D}(E)$ the collection of all cubes $Q^{k}_{j}$ , that is,

\mathbb{D}\coloneqq\cup_{k}\mathbb{D}_{k}.

If $E$ is bounded, we ignore cubes where $2^{-k}\gtrsim\operatorname{diam}(\partial\Omega)$ (in particular, where $a_{0}2^{-k}\geq\operatorname{diam}(\partial\Omega)$ ). Given a cube $Q=Q_{j}^{k}\in\mathbb{D}$ , we define the side-length of $Q$ as $\ell(Q)\coloneqq 2^{-k}$ . By Ahlfors regularity and property (v) in Lemma 2.14, we know that $\ell(Q)\approx\operatorname{diam}(Q)$ and $\mathcal{H}^{n}(Q)\approx\ell(Q)^{n}$ . Given $Q\in\mathbb{D}$ and $m\in\mathbb{Z}$ , we set

\mathbb{D}_{Q}\coloneqq\left\{Q^{\prime}\in\mathbb{D}\colon Q^{\prime}\subseteq Q\right\},\qquad\mathbb{D}_{m,Q}\coloneqq\left\{Q^{\prime}\in\mathbb{D}_{Q}\colon\ell(Q^{\prime})=2^{-m}\ell(Q)\right\}.

We call the cubes in the collection $\mathbb{D}_{1,Q}$ the children of $Q$ . Notice that by Ahlfors regularity and property (v) in Lemma 2.14, each cube has a uniformly bounded number of children.

Given a cube $Q=Q_{j}^{k}\in\mathbb{D}$ , we call the point $z_{j}^{k}\in Q$ in property (v) in Lemma 2.14 the center of $Q$ , denote $x_{Q}\coloneqq z_{j}^{k}$ , and set

\Delta_{Q}\coloneqq E\cap B(z_{j}^{k},a_{0}\ell(Q)).

We use Lemma 2.14 to decompose $\partial\Omega$ , so that $\mathbb{D}(E)=\mathbb{D}(\partial\Omega)\eqqcolon\mathbb{D}$ . For each $Q\in\mathbb{D}$ , we let $X_{Q}$ be the corkscrew point relative to $x_{Q}$ at scale $10^{-5}a_{0}\ell(Q)$ . We have $B(X_{Q},\gamma 10^{-5}a_{0}\ell(Q))\subset B(x_{Q},10^{-5}a_{0}\ell(Q))\cap\Omega$ , where $\gamma$ is the corkscrew constant in Definition 2.1.

For many of our techniques, it is important that we show that some collections of dyadic cubes are quantitatively small in the following sense:

Definition 2.15 (Carleson packing condition).

Let $\mathbb{D}$ be a dyadic system on $\partial\Omega$ and let $\mathcal{A}\subset\mathbb{D}$ . We say that $\mathcal{A}$ satisfies a Carleson packing condition if there exists a constant $C\geq 1$ such that for any $Q_{0}\in\mathbb{D}$ we have

\displaystyle\sum_{Q\in\mathcal{A},Q\subset Q_{0}}\sigma(Q)\leq C\sigma(Q_{0}).

We denote the smallest such constant $C$ by $\mathcal{C}_{\mathcal{A}}$ .

Next, we use a standard decomposition of $\Omega$ into Whitney cubes (see e.g. [Ste70, Chapter VI]), and then associate a collection of such Whitney cubes to each boundary cube to construct suitable Whitney-type regions. These Whitney regions are modeled after regions of the type $Q\times(\ell(Q)/2,\ell(Q))$ (that is, the upper halves of Carleson boxes) in the simpler geometry of the upper half-space. For this, we recall the construction found in [HM14] noting that we make some changes to the notation therein (following the notation of more recent papers, e.g. [HMM16]). We let $\mathcal{W}=\{I\}_{I}$ denote a Whitney decomposition of $\Omega$ , with the properties that each $I$ is a closed $(n+1)$ -dimensional cube satisfying

\displaystyle 4\operatorname{diam}(I)\leq\operatorname{dist}(4I,\partial\Omega)\leq\operatorname{dist}(I,\partial\Omega)\leq 40\operatorname{diam}(I),

where $4I$ is the standard concentric Euclidean dilate of a cube; the interiors of the cubes $I$ are disjoint, and for all $I_{1},I_{2}\in\mathcal{W}$ with $I_{1}\cap I_{2}\neq\mbox{{\O}}$ we have

\displaystyle\frac{1}{4}\operatorname{diam}(I_{1})\leq\operatorname{diam}(I_{2})\leq 4\operatorname{diam}(I_{1}).

For $I\in\mathcal{W}$ we let $\ell(I)$ denote the side length of $I$ .

For each cube $Q\in\mathbb{D}$ and constant $K\geq K_{0}$ , with $K_{0}$ to be described momentarily, we associate an initial collection of Whitney cubes

W_{Q}(K)\coloneqq\{I\in\mathcal{W}\colon K^{-1}\ell(I)\leq\ell(Q)\leq K\ell(I),\operatorname{dist}(I,Q)\leq K\ell(Q)\}.

We choose $K_{0}$ depending on the constants in the corkscrew condition and the Ahlfors regularity condition, insisting on two conditions being met:

(1)

If $X\in\Omega$ with $\operatorname{dist}(X,\partial\Omega)\leq 10^{5}\operatorname{diam}(\partial\Omega)$ then $X\in I\in W_{Q}(K)$ for some $Q\in\mathbb{D}$ .
(2)

For any $Q\in\mathbb{D}$ , we have $B(X_{Q},\operatorname{dist}(X_{Q},\partial\Omega)/2)\subseteq\cup_{I\in W_{Q}(K)}I$ , and if $Q^{\prime}\in\mathbb{D}$ is another cube such that $Q^{\prime}\subset Q$ with $\ell(Q^{\prime})=\tfrac{1}{2}\ell(Q)$ , then we also have $B(X_{Q^{\prime}},\operatorname{dist}(X_{Q^{\prime}},\partial\Omega)/2)\subseteq\cup_{I\in W_{Q}(K)}I$ .

Of course, condition (1) above is automatically satisfied if $\operatorname{diam}(\partial\Omega)=\infty$ .

Following [HM14, Section 3], we augment the collection $\mathcal{W}_{Q}(K)$ as follows. For each $I\in W_{Q}(K)$ , we take a Harnack chain $H(I)$ from the center of $I$ to the corkscrew point $X_{Q}$ , and we let $W_{Q,I}(K)$ be be the collection of Whitney cubes in $\mathcal{W}$ that meet at least one ball in the chain $H(I)$ . We then set $\mathcal{W}^{*}_{Q}(K)\coloneqq\bigcup_{I\in W_{Q}(K)}W_{Q,I}(K)$ . Finally, for a small dimensional parameter $\tau$ (this is the parameter $\lambda$ in [HM14, Section 3]), we define the Whitney region relative to $Q$ as

\displaystyle U_{Q}=U_{Q}(K_{0})\coloneqq\bigcup_{I\in\mathcal{W}^{*}_{Q}(K_{0})}(1+\tau)I.

By construction, we know that if $X\in U_{Q}$ , then

(2.16)

K_{1}^{-1}\ell(Q)\leq\operatorname{dist}(X,\partial\Omega)\leq\operatorname{dist}(X,Q)\leq K_{1}\ell(Q),

where $K_{1}$ depends on $K_{0}$ , the dimension and the Harnack chain condition. For $\kappa\gg K_{0}$ to be chosen, we also define the following fattened versions of the Whitney regions:

U_{Q}^{*}=U_{Q}(\kappa)=\bigcup_{I\in\mathcal{W}^{*}_{Q}(\kappa)}(1+\tau)I,\qquad U_{Q}^{**}=\bigcup_{I\in\mathcal{W}^{*}_{Q}(\kappa)}(1+2\tau)I,

that is, $U_{Q}^{*}$ is constructed the same way as $U_{Q}$ but we replace the constant $K_{0}$ by $\kappa$ (similarly for $U_{Q}^{**}$ ). We describe the reasoning and choice of $\kappa$ in the next subsection. We note that for $\tau$ small enough the regions $U_{Q}$ , $U^{*}_{Q}$ and $U_{Q}^{**}$ have bounded overlaps, that is, for a collection of dyadic cubes $\mathscr{Q}$ and the $(n+1)$ -dimensional Lebesgue measure $|\cdot|$ we have $|\bigcup_{Q\in\mathscr{Q}}U_{Q}^{**}|\approx\sum_{Q\in\mathscr{Q}}|U_{Q}^{**}|$ .

Using the Whitney regions above, we can now define objects like sawtooth regions, Carleson boxes and dyadic cones. Let $Q_{0}\in\mathbb{D}$ be a fixed cube and $\mathcal{F}\subset\mathbb{D}_{Q_{0}}$ a collection of pairwise disjoint cubes. We set

\mathbb{D}_{\mathcal{F},Q_{0}}\coloneqq\mathbb{D}_{Q_{0}}\setminus\cup_{Q\in\mathcal{F}}\mathbb{D}_{Q}.

We then define the local sawtooth relative to $\mathcal{F}$ (and its fattened version) as

\Omega_{\mathcal{F},Q_{0}}\coloneqq\operatorname{int}\Big{(}\bigcup_{Q\in\mathbb{D}_{\mathcal{F},Q_{0}}}U_{Q}\Big{)},\qquad\Omega^{*}_{\mathcal{F},Q_{0}}\coloneqq\operatorname{int}\Big{(}\bigcup_{Q\in\mathbb{D}_{\mathcal{F},Q_{0}}}U^{*}_{Q}\Big{)}.

In the special case where $\mathcal{F}=\mbox{{\O}}$ , we write $T_{Q_{0}}=\Omega_{\mathcal{F},Q_{0}}$ and $T^{*}_{Q_{0}}=\Omega^{*}_{\mathcal{F},Q}$ , that is,

T_{Q_{0}}\coloneqq\operatorname{int}\Big{(}\bigcup_{Q\in\mathbb{D}_{Q_{0}}}U_{Q}\Big{)},\qquad T_{Q_{0}}^{*}\coloneqq\operatorname{int}\Big{(}\bigcup_{Q\in\mathbb{D}_{Q_{0}}}U^{*}_{Q}\Big{)},\qquad T_{Q_{0}}^{**}\coloneqq\operatorname{int}\Big{(}\bigcup_{Q\in\mathbb{D}_{Q_{0}}}U^{**}_{Q}\Big{)}

We call $T_{Q_{0}}$ the Carleson box relative to $Q_{0}$ and $T_{Q_{0}}^{*}$ and $T_{Q_{0}}^{**}$ its fattened versions. Given a cube $Q_{0}\in\mathbb{D}$ and a point $x\in Q_{0}$ , we also define the (truncated) dyadic cone at $x\in\partial\Omega$ , $\Gamma(x)$ , by setting

(2.17)

\displaystyle\Gamma(x)\coloneqq\Gamma_{Q_{0}}(x)\coloneqq\operatorname{int}\Big{(}\bigcup_{Q\in\mathbb{D}_{Q_{0}}\colon x\in Q}U_{Q}\Big{)}.

Notice that $\Gamma(x)=\Omega_{\mathcal{F}_{x},Q_{0}}$ , where $\mathcal{F}_{x}$ is the collection of maximal²²2Since the sizes of the cubes in $\{Q\in\mathbb{D}_{Q_{0}}\colon x\notin Q\}$ are bounded from above, we can always choose the maximal cubes. In the case of non-truncated dyadic cones (that is, cones of the form $\operatorname{int}\left(\bigcup_{Q\in\mathbb{D}\colon x\in Q}U_{Q}\right)$ ), we can choose the maximal cubes if the dyadic system $\mathbb{D}$ forms a tree-like structure. This structure can be achieved by choosing the center points of the dyadic cubes in a suitable way (see, for example, [HT14, Theorem 2.4]). (and hence, disjoint) cubes in the collection $\{Q\in\mathbb{D}_{Q_{0}}\colon x\notin Q\}$ . It is straightforward to verify that there exists uniform constants $m_{1},m_{2}>1$ such that $\Gamma(x)$ contains a truncated version of the cone $\widetilde{\Gamma}^{m_{1}}(x)$ and it is contained in a truncated version of the cone $\widetilde{\Gamma}^{m_{2}}(x)$ , where $\widetilde{\Gamma}^{m_{1}}(x)$ and $\widetilde{\Gamma}^{m_{2}}(x)$ are cones of the type (2.9). Thus, since the aperture of the cones is not important for us, we mostly use dyadic cones when studying non-tangential convergence.

By the following lemma, the Whitney regions, sawtooth regions, Carleson boxes, and truncated dyadic cones inherit many quantitative geometric properties from $\Omega$ :

Lemma 2.18 ([HM14, Lemma 3.61]).

Suppose that $\Omega\subset\mathbb{R}^{n+1}$ is a uniform domain with $n$ -Ahlfors regular boundary. Let $Q_{0}\in\mathbb{D}(\partial\Omega)$ be a cube and $\mathcal{F}\subset\mathbb{D}_{Q_{0}}$ be a collection of pairwise disjoint cubes. Then $\Omega_{\mathcal{F},Q}$ and $\Omega^{*}_{\mathcal{F},Q}$ are uniform domains with $n$ -Ahlfors regular boundary whose structural constants depend only on the dimension, the structural constants of $\Omega$ , and the constant $\kappa$ . In particular, the Whitney regions $U_{Q}$ , $U_{Q}^{*}$ and $U_{Q}^{**}$ , the Carleson boxes $T_{Q}$ , $T_{Q}^{*}$ and $T_{Q}^{**}$ and the truncated dyadic cones $\Gamma(x)$ are uniform domains with $n$ -Ahlfors regular boundaries, with structural constants depending only on the dimension and the structural constants of $\Omega$ and the constant $\kappa$ .

2.4. The choice of the parameter $\kappa$

In contrast to the setting of the upper half space, we do not define the sawtooths by removing Whitney regions. This is due to the overlaps of the regions $U_{Q}$ : we may encounter situations where for $Q_{0}\in\mathbb{D}(\partial\Omega)$ and a collection of pairwise disjoint cubes $\mathcal{F}\subset\mathbb{D}_{Q_{0}}$ there exists a cube $Q\in\mathbb{D}_{\mathcal{F},Q_{0}}$ such that $\overline{U_{Q}}$ does not contribute to the boundary of $\Omega_{\mathcal{F},Q}$ . That being said, if $\kappa$ is chosen large enough, then the fattened Whitney region $U_{Q}^{*}$ meets the boundary of the unfattened region $\Omega_{\mathcal{F},Q}$ on a portion roughly the measure of $Q$ . We will consider this in Section 4 where it will be convenient for us, but we prove the technical estimates that give this property below.

Let us fix a cube $Q\in\mathbb{D}$ . Recall that $X_{Q}$ is a corkscrew point relative to $x_{Q}\in Q$ at scale $r_{Q}\coloneqq 10^{-5}a_{0}\ell(Q)$ and $\Delta_{Q}=B(x_{Q},10^{5}r_{Q})\cap\partial\Omega\subset Q$ is the surface ball associated to $Q$ . We let $\hat{x}_{Q}\in\partial\Omega$ denote a touching point for $X_{Q}$ , that is, a point such that $|X_{Q}-\hat{x}_{Q}|=\delta(X_{Q})$ . By triangle inequality, $|x_{Q}-\hat{x}_{Q}|<2r_{Q}$ , and thus,

(2.19)

\hat{\Delta}_{Q}=B(\hat{x}_{Q},10^{3}r_{Q})\cap\partial\Omega\subset\Delta_{Q}\subseteq Q.

For every $\theta\in(0,1)$ , we let

\displaystyle P_{Q}(\theta)\coloneqq\hat{x}_{Q}+\theta(X_{Q}-\hat{x}_{Q})

be the “ $\theta$ -point” on the directed line segment from $\hat{x}_{Q}$ to $X_{Q}$ . Then, by definitions,

(2.20)

\gamma\theta r_{Q}\leq|P_{Q}(\theta)-\hat{x}_{Q}|=\operatorname{dist}(P_{Q}(\theta),\partial\Omega)\leq\theta r_{Q}.

Lemma 2.21.

There exists $\theta_{0}\in(0,1)$ depending on $K_{0}$ and the structural constants such that if for some $Q^{\prime}\in\mathbb{D}$ we have that

\displaystyle B\big{(}P_{Q}(\theta_{0}),\tfrac{\gamma\theta_{0}}{10}r_{Q}\big{)}\cap U_{Q^{\prime}}\neq\mbox{{\O}},

then $Q^{\prime}\subseteq Q$ and $\ell(Q^{\prime})<\ell(Q)$ .

Proof.

Suppose that $\theta_{0}\leq 1/4C(K_{1})^{2}$ for a large structural constant $C\geq 1$ and

(2.22)

\displaystyle X\in B\big{(}P_{Q}(\theta_{0}),\tfrac{\gamma\theta_{0}}{10}r_{Q}\big{)}\cap U_{Q^{\prime}}

for some $Q^{\prime}$ , where $K_{1}$ is the constant in (2.16). By (2.16), (2.22) and (2.20), it holds that

	$\displaystyle\ell(Q^{\prime})$	$\displaystyle\leq K_{1}\operatorname{dist}(X,\partial\Omega)$
(2.23)			$\displaystyle\leq K_{1}\bigg{(}\operatorname{dist}(P_{Q}(\theta_{0}),\partial\Omega)+\frac{\gamma\theta_{0}}{10}r_{Q}\bigg{)}\leq 2K_{1}\theta_{0}r_{Q}=2K_{1}\theta_{0}10^{-5}a_{0}\ell(Q).$

In particular, we have $\ell(Q^{\prime})<\ell(Q)$ .

To show that $Q^{\prime}\subset Q$ , we first notice that we have

\displaystyle|X-x_{Q^{\prime}}|\leq CK_{1}\ell(Q^{\prime})

for a structural constant $C\geq 1$ by (2.16) and the fact that $\operatorname{diam}(Q^{\prime})\approx\ell(Q^{\prime})$ . This and (2.23) then give us

\displaystyle|X-x_{Q^{\prime}}|\leq 2C(K_{1})^{2}\theta_{0}10^{-5}a_{0}\ell(Q).

Thus, by (2.20) and (2.22), it holds that

\displaystyle|X-\hat{x}_{Q}|\leq 2\theta_{0}r_{Q}=2\theta_{0}10^{-5}a_{0}\ell(Q).

Combining the previous two inequalities then gives us

(2.24)

\displaystyle|\hat{x}_{Q}-x_{Q^{\prime}}|\leq 4C(K_{1})^{2}\theta_{0}10^{-5}a_{0}\ell(Q)=4C(K_{1})^{2}\theta_{0}r_{Q}.

In particular,

\displaystyle|x_{Q^{\prime}}-x_{Q}|\leq|x_{Q^{\prime}}-\hat{x}_{Q}|+|\hat{x}_{Q}-x_{Q}|\leq 4C(K_{1})^{2}\theta_{0}r_{Q}+2r_{Q}<3r_{Q}<a_{0}\ell(Q),

by (2.24), the fact that $|\hat{x}_{Q}-x_{Q}|<2r_{Q}$ , and the choice $\theta_{0}\leq 1/4C(K_{1})^{2}$ . Thus, $x_{Q^{\prime}}\in\Delta_{Q}\subset Q$ . Since $Q^{\prime}\cap Q\neq\mbox{{\O}}$ and $\ell(Q^{\prime})<\ell(Q)$ , we know that $Q^{\prime}\subset Q$ , which is what we wanted. ∎

Let us then fix $\theta_{0}$ so that Lemma 2.21 holds. For $Q\in\mathbb{D}$ , we set

\displaystyle\Xi_{Q}\coloneqq\bigcup_{\theta\in[\theta_{0},1]}B\big{(}P_{Q}(\theta),\tfrac{\gamma\theta_{0}}{10}r_{Q}\big{)},

which is a cylinder-like object. We get the following straightforward lemma:

Lemma 2.25.

Let $Q\in\mathbb{D}$ be a fixed cube and let $\mathscr{Q}_{Q}$ be the collection of cubes that share the same dyadic parent as $Q$ , that is,

\displaystyle\mathscr{Q}_{Q}\coloneqq\{P\in\mathbb{D}\colon P,Q\subset Q_{0}\text{ for a cube }Q_{0}\in\mathbb{D}\text{ such that }\ell(P)=\ell(Q)=\tfrac{1}{2}\ell(Q_{0})\}.

Let $\kappa\gg\max\{K_{0},(\theta_{0})^{-1}\}$ and $X\in\Xi_{P}$ for some $P\in\mathscr{Q}_{Q}$ . Then $X\in U_{Q}^{*}$ .

Proof.

Let $\kappa\gg\max\{K_{0},(\theta_{0})^{-1}\}$ and $X\in\Xi_{P}$ for some $P\in\mathscr{Q}_{Q}$ . By the Whitney decomposition, there exists a Whitney cube $I\in\mathcal{W}$ such that $X\in I$ . By the definition of $U_{Q}^{*}$ , it is enough to show that $I\in\mathcal{W}_{Q}^{*}(\kappa)$ .

By (2.20) and the definitions, we first notice that

\displaystyle\ell(I)\approx\operatorname{dist}(X,\partial\Omega)\approx\gamma\theta_{0}r_{P}\approx\theta_{0}\ell(P)=\theta_{0}\ell(Q)

with uniformly bounded implicit constants. In particular, since $\theta_{0}\gg\tfrac{1}{\kappa}$ and $\kappa\gg 1$ , we get

\displaystyle\frac{1}{\kappa}\ell(I)\leq\ell(Q)\leq\kappa\ell(I).

On the other hand, since $P\in\mathscr{Q}_{Q}$ , we know that

\displaystyle\operatorname{dist}(I,Q)\lesssim\operatorname{dist}(I,P)+\ell(Q),

and by (2.20), (2.19) and the fact that $X\in\Xi_{P}$ , we know that

\displaystyle\operatorname{dist}(I,P)\leq\operatorname{dist}(X,P)\leq|X-\hat{x}_{P}|\leq\frac{\gamma\theta_{0}}{10}r_{P}+|X_{P}-\hat{x}_{P}|\leq 2r_{P}\leq 2\ell(P)=2\ell(Q).

In particular, since $\kappa\gg 1$ , we have

\displaystyle\operatorname{dist}(I,Q)\leq\kappa\ell(Q).

Thus, $I\in\mathcal{W}_{Q}^{*}(\kappa)$ , which proves the claim. ∎

Let us also record the following simple lemma for future use:

Lemma 2.26.

Let $Q\in\mathbb{D}$ and let $X_{Q}\in\Omega$ be a corkscrew point, $\hat{x}_{Q}\in\partial\Omega$ be a touching point and $r_{Q}=10^{-5}a_{0}\ell(Q)$ as above. If $Y\in B(P_{Q}(\theta),r_{Q})$ for some $\theta\in[0,1]$ and $\hat{y}\in\partial\Omega$ is a point such that $|Y-\hat{y}|=\operatorname{dist}(Y,\partial\Omega)$ , then $\hat{y}\in\hat{\Delta}_{Q}\subseteq\Delta_{Q}\subseteq Q$ .

Proof.

By definitions and using (2.20) several times, we get

|\hat{y}-\hat{x}_{Q}|\leq|\hat{y}-Y|+|Y-P_{Q}(\theta)|+|P_{Q}(\theta)-\hat{x}_{Q}|<\operatorname{dist}(Y,\partial\Omega)+2r_{Q}\\ \leq|Y-P_{Q}(\theta)|+\operatorname{dist}(P_{Q}(\theta),\partial\Omega)+2r_{Q}<r_{Q}+r_{Q}+2r_{Q}=4r_{Q},

and thus, $\hat{y}\in\hat{\Delta}_{Q}\subset\Delta_{Q}\subset Q$ by (2.19). ∎

2.5. Elliptic PDE estimates

Here we collect some of the standard estimates for divergence form elliptic operators with real coefficients that will be used throughout the paper. In this section, $\Omega$ always denotes a uniform domain in $\mathbb{R}^{n+1}$ , $n\geq 1$ , with $n$ -Ahlfors regular boundary. We recall that a divergence form elliptic operator is of the form

\displaystyle L(\cdot)\coloneqq-\mathop{\operatorname{div}}\nolimits(A\nabla\cdot),

viewed in the weak sense, where $A$ is a uniformly elliptic matrix, that is, $A=(a_{i,j})_{i,j=1}^{n+1}$ is an $(n+1)\times(n+1)$ matrix-valued function on $\mathbb{R}^{n+1}$ and there exists a constant $\Lambda$ , the ellipticity parameter, such that

\displaystyle\Lambda^{-1}|\xi|^{2}\leq A(X)\xi\cdot\xi,\qquad\text{and}\qquad\|a_{i,j}\|_{L^{\infty}(\mathbb{R}^{n+1})}\leq\Lambda,

for all $\xi,\zeta\in\mathbb{R}^{n+1}$ and almost every $X\in\Omega$ . We say that a constant depends on ellipticity if it depends on $\Lambda$ . Given an open set $\mathcal{O}\subset\mathbb{R}^{n+1}$ we say a function $u\in W^{1,2}_{\operatorname{loc}}(\mathcal{O})$ is a solution to $Lu=0$ in $\mathcal{O}$ if

\displaystyle\iint_{\mathcal{O}}A\nabla u\cdot\nabla\varphi\,dX=0,\quad\text{ for every }\varphi\in C_{\operatorname{c}}^{\infty}(\mathcal{O}).

The most fundamental estimate for solutions to divergence form elliptic equations is the following local energy inequality.

Lemma 2.27 (Caccioppoli Inequality).

Let $L=-\mathop{\operatorname{div}}\nolimits A\nabla$ be a divergence form elliptic operator and $u$ a solution to $Lu=0$ in an open set $\mathcal{O}$ . If $a>0$ and $B$ is a ball such that $(1+a)B\subset\mathcal{O}$ then

\displaystyle\iint_{B}|\nabla u|\,dX\lesssim r^{-2}\iint_{(1+a)B}u^{2}\,dX,

where the implicit constant depends only on $a$ , dimension and ellipticity.

Solutions to divergence form elliptic equations are locally Hölder continuous.

Lemma 2.28 (Hölder continuity of solutions, [DG57, Nas58]).

Let $L=-\mathop{\operatorname{div}}\nolimits A\nabla$ be a divergence form elliptic operator and $u$ a non-negative solution to $Lu=0$ in an open set $\mathcal{O}$ . Suppose that $B=B(X_{0},R)$ is a ball such that $\lambda B\coloneqq B(X_{0},2\lambda R)\subset\mathcal{O}$ for $\lambda>1$ . Then we have

\displaystyle|u(X)-u(Y)|\leq C\left(\frac{|X-Y|}{\lambda R}\right)^{\alpha}\left(\mathchoice{{\vbox{\hbox{$\textstyle\raisebox{-8.61108pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-7.09656pt}}{{\vbox{\hbox{$\scriptstyle\raisebox{-6.02777pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-6.09258pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-4.29826pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-3.54826pt}}\!\iint_{2\lambda B}|u|^{2}\,dY\right)^{1/2}\quad\text{ for all }X,Y\in B,

where $\alpha$ and $C$ depend only on dimension and ellipticity.

Another celebrated result is Moser’s Harnack inequality for non-negative solutions.

Lemma 2.29 (Harnack inequality [Mos61]).

Let $L=-\mathop{\operatorname{div}}\nolimits A\nabla$ be a divergence form elliptic operator and $u$ a non-negative solution to $Lu=0$ in an open set $\mathcal{O}$ . If $B$ is a ball such that $2B\subset\mathcal{O}$ then $\sup_{B}u\leq C\inf_{B}u$ , where $C$ depends only on dimension and ellipticity.

We now turn our attention to the elliptic measure, for which we borrow the setting of [AGMT22]. Consider the compactified space $\overline{\mathbb{R}}^{n+1}=\mathbb{R}^{n+1}\cup\{\infty\}$ ; following [HKM93, Section 9], we will understand all topological notions with respect to this space. Hence, for instance, if $\Omega$ is unbounded, then $\infty\in\partial\Omega$ , and the functions in the space $C(\partial\Omega)$ are assumed continuous and real-valued, so that all functions in $C(\partial\Omega)$ lie in $L^{\infty}(\partial\Omega)$ even if $\partial\Omega$ is unbounded.

Given a domain $\Omega$ and a divergence form elliptic operator $L$ , we let $\omega_{L,\Omega}^{X}$ denote the elliptic measure with pole at $X\in\Omega$ . That is, by the Riesz representation theorem, for every $X\in\Omega$ there exists a probability measure $\omega_{L,\Omega}^{X}$ such that if $f\in C_{\operatorname{c}}(\partial\Omega)$ , then the solution to $Lu=0$ in $\Omega$ with $u\in C(\overline{\Omega})$ and $u=f$ on $\partial\Omega$ , constructed via Perron’s method, satisfies

(2.30)

u(X)=u_{f}(X)=\int_{\partial\Omega}f(y)\,d\omega^{X}_{L,\Omega}(y).

When the context is clear, we simply denote $\omega^{X}\coloneqq\omega^{X}_{L,\Omega}$ and, with slight abuse of terminology, call the family of elliptic measures $\omega=\omega_{L}=\omega_{L,\Omega}\coloneqq\{\omega^{X}\}_{X}$ just the elliptic measure.

Our main results consider characterizations and implications given by quantitative absolute continuity of elliptic measure in the sense of Muckenhoupt’s $A_{\infty}$ condition [Muc72, CF74]:

Definition 2.31 ( $A_{\infty}$ for elliptic measure).

Let $L$ be a divergence form elliptic operator in $\Omega$ . We say that the elliptic measure $\omega=\omega_{L,\Omega}$ satisfies the $A_{\infty}$ condition with respect to surface measure (denote $\omega\in A_{\infty}(\sigma)$ ) if there exist constants $C\geq 1$ and $s>0$ such that if $B\coloneqq B(x,r)$ with $x\in\partial\Omega$ and $r\in(0,\operatorname{diam}(\partial\Omega)/4)$ and $A\subset\Delta\coloneqq B\cap\partial\Omega$ is a Borel set, then

\displaystyle\omega^{Y}(A)\leq C\left(\frac{\sigma(A)}{\sigma(\Delta)}\right)^{s}\omega^{Y}(\Delta),\qquad\text{for every }Y\in\Omega\setminus 4B.

We refer to $C$ and $s$ here together as the “ $\omega_{L}\in A_{\infty}(\sigma)$ constants”.

Next we discuss the Green’s function and its properties.

Definition 2.32 (Green’s function).

Let $L=-\operatorname{div}A\nabla$ be a not necessarily symmetric divergence form elliptic operator with bounded measurable coefficients. There exists a unique non-negative function $G=G_{L}=G_{L,\Omega}:\Omega\times\Omega\rightarrow\mathbb{R}$ , called Green’s function for $L$ , satisfying the following properties:

(i)

For each $X,Y\in\Omega$ ,

(2.33)

0\leq G(X,Y)\lesssim\left\{\begin{matrix}|X-Y|^{1-n},&\qquad n\geq 2,&\qquad X\neq Y,\\[2.84526pt] 1,&\qquad n=1,&\qquad|X-Y|\geq\delta(Y)/10,\\[2.84526pt] \log\big{(}\frac{\delta(Y)}{|X-Y|}\big{)},&\qquad n=1,&\qquad|X-Y|\leq\delta(Y)/2.\end{matrix}\right.

(ii)

For every $a\in(0,1)$ there exists $c_{a}$ such that

(2.34)

G(X,Y)\geq\left\{\begin{matrix}c_{a}|X-Y|^{1-n},&\qquad n\geq 2,&\qquad|X-Y|\leq a\delta(Y),\\[2.84526pt] c_{a}\log\big{(}\frac{\delta(Y)}{|X-Y|}\big{)},&\qquad n=1,&\qquad|X-Y|\leq a\delta(Y).\end{matrix}\right.

(iii)

For each $Y\in\Omega$ , $G(\cdot,Y)\in C(\overline{\Omega}\backslash\{Y\})\cap W^{1,2}_{\operatorname{loc}}(\Omega\backslash\{Y\})$ and $G(\cdot,Y)|_{\partial\Omega}\equiv 0$ .

(iv)

For each $X\in\Omega$ , the identity $LG(\cdot,X)=\delta_{X}$ holds in the distributional sense; that is,

\int_{\Omega}A(Y)\nabla_{Y}G(Y,X)\nabla\Phi(Y)\,dY=\Phi(X),\qquad\text{for any }\Phi\in C_{c}^{\infty}(\Omega).

(v)

For each $X,Y\in\Omega$ with $X\neq Y$ , if $L^{*}=-\operatorname{div}A^{T}\nabla$ , then

$G_{L}(X,Y)=G_{L^{*}}(Y,X).$

If $n\geq 2$ , then it has been known for a long time that a non-negative Green’s function exists for any domain, without any further regularity assumptions on the geometry [GW82, HK07]. If $n=1$ , the situation has been more challenging; for instance, key Sobolev embeddings, available when $n\geq 2$ , fail when $n=1$ , and the fundamental solution changes sign when $n=1$ [KN85]; nevertheless, the paper [DK09] shows the construction of a Green’s function for any domain with either finite volume or finite width, and also, for the domain above a Lipschitz graph, improving on the result of [DM95] (but non-negativity is not shown in these works). For our setting of uniform domains $\Omega\subset\mathbb{R}^{n+1}$ , $n\geq 1$ , with $n$ -Ahlfors regular boundary, the unified (for $n=1$ and $n\geq 2$ ) existence of the non-negative Green’s function for arbitrary divergence form elliptic operators $L$ of merely bounded measurable coefficients with the properties stated above follows from the much more general, recent construction in [DFM, Theorem 14.60 and Lemma 14.78].

The Green’s function and the elliptic measure are related through the following Riesz formula: For every $F\in C_{\operatorname{c}}^{\infty}(\mathbb{R}^{n+1})$ , one has that

(2.35)

\int_{\partial\Omega}F(y)\,d\omega^{X}(y)=F(X)-\iint_{\Omega}A^{T}\nabla_{Y}G_{L}(X,Y)\cdot\nabla_{Y}F(Y)\,dY.

We need several estimates from the literature for the elliptic measure and Green’s function in our proofs and we list these estimates below. Although these have appeared in several works in the literature [CFMS81, HKM93, AGMT22], we cite [DFM] for their unified consideration of the cases $n=1$ and $n\geq 2$ and arbitrary elliptic operators on uniform domains with Ahlfors regular boundary. The first estimate is a non-degeneracy estimate for elliptic measure.

Lemma 2.36 (Bourgain’s estimate, [DFM, Lemma 15.1]).

Let $x\in\partial\Omega$ and $r\in(0,\operatorname{diam}(\partial\Omega)]$ and let $Y_{x,r}$ be a corkscrew point relative to $x$ at scale $r$ . We have

\displaystyle\omega^{Y_{x,r}}(\Delta(x,r))\geq c,

where $c$ depends only on dimension, ellipticity, and the Ahlfors regularity constant. Here and below $\Delta(x,r)=B(x,r)\cap\partial\Omega$ .

The elliptic measure is locally doubling in the following sense.

Lemma 2.37 ((Local) doubling property, [DFM, Lemma 15.43]).

Let $x\in\partial\Omega$ and $r\in(0,\operatorname{diam}(\Omega)]$ . If $Y\in\Omega\setminus B(x,4r)$ then $\omega^{Y}(\Delta(x,2r))\leq C\omega^{Y}(\Delta(x,r))$ , where $C$ depends on dimension, ellipticity, Harnack chain, corkscrew and Ahlfors regularity constants.

The following estimate allows us to connect the Green’s function and elliptic measure in a quantitative way:

Lemma 2.38 (CFMS estimate, [DFM, Lemma 15.28]).

If $x\in\partial\Omega$ and $r\in(0,\operatorname{diam}(\partial\Omega)]$ then

\displaystyle\frac{G_{L}(X,Y_{x,r})}{r}\approx\frac{\omega^{X}_{L}(\Delta(x,r))}{r^{n}}

for every $X\in\Omega\setminus B(x,2r)$ , where the implicit constants depend on dimension, ellipticity, Harnack chain, corkscrew and Ahlfors regularity constants, and $Y_{x,r}$ is any corkscrew point relative to $x$ at scale $r$ .

Non-negative solutions $u$ to $Lu=0$ that vanish on an open subset of the boundary of a uniform domain must do so at the same rate:

Lemma 2.39 (Boundary Harnack Principle, [DFM, Theorem 15.64]).

Let $x\in\partial\Omega$ and $r\in(0,\operatorname{diam}(\partial\Omega)]$ , and let $u$ and $v$ be positive functions such that $Lu=Lv=0$ in $\Omega\cap B(x,4r)$ that vanish continuously on $\partial\Omega\cap B(x,4r)$ . Then

\displaystyle\frac{u(X)}{v(X)}\approx\frac{u(Y_{x,r})}{v(Y_{x,r})},\qquad\text{for all }X\in B(x,r)\cap\Omega,

where $Y_{x,r}$ is a corkscrew point relative to $x$ at scale $r$ . The implicit constants depend on dimension, ellipticity, Harnack chain, corkscrew and Ahlfors regularity constants.

We have the following standard consequence of the boundary Harnack Principle.

Lemma 2.40 (Change of Poles, [DFM, Lemma 15.61]).

Let $x\in\partial\Omega$ , $r\in(0,\operatorname{diam}(\partial\Omega))$ , $Y_{x,r}$ a corkscrew point relative to $x$ at scale $r$ , and $E\subset\Delta(x,r)$ a Borel set. Then

\omega^{Y_{x,r}}(E)\approx\frac{\omega^{X}(E)}{\omega^{X}(\Delta(x,r))},\qquad\text{for each }X\in\Omega\backslash B(x,2r),

where the implicit constants depend on dimension, ellipticity, and structural constants of $\Omega$ .

Finally, solutions $u$ to $Lu=0$ that vanish continuously on the boundary do so at a Hölder rate:

Lemma 2.41 ([DFM, Lemma 11.32]).

Let $x\in\partial\Omega$ and $r\in(0,\operatorname{diam}(\partial\Omega)]$ , and let $u$ be a solution to $Lu=0$ in $\Omega\cap B(x,4r)$ that vanishes continuously on $\partial\Omega\cap B(x,4r)$ . We have the bound

\displaystyle|u(X)|\leq C\left(\frac{\operatorname{dist}(X,\partial\Omega)}{r}\right)^{\alpha}\sup_{Y\in B(x,4r)}|u(Y)|

for every $X\in B(x,r)\cap\Omega$ , where $\alpha$ and $C$ depend on dimension, the Harnack chain, corkscrew and Ahlfors regularity constants.

3. Set-up for Theorem 1.1

In this section we provide some of the preliminary estimates and observations required to prove Theorem 1.1. We divide this section into two subsections. The first records Carleson measure estimates and non-tangential convergence in our setting. The second subsection contains a few lemmas, which are roughly adapted from ideas in [DJK84] and play a crucial role in our analysis. Throughout this section, we suppose that the assumptions of Theorem 1.1 hold; that is, $\Omega$ is a uniform domain with Ahlfors regular boundary and $L$ is a not necessarily symmetric divergence form elliptic operator such that $\omega_{L}\in A_{\infty}(\sigma)$ . For the sequel, recall that $\delta(\cdot)\coloneqq\operatorname{dist}(\cdot,\partial\Omega)$ .

3.1. CME and non-tangential convergence

One of the key tools in most of the constructions of $\varepsilon$ -approximators (see, for example, [HMM16]) is $L^{2}$ -type Carleson measure estimates (CME), that is, Carleson properties (see Definition 2.7) of measures $\mu_{u}$ such that $d\mu_{u}=|\nabla u|^{2}\delta(Y)\,dY$ for a solution $u$ to $Lu=0$ . Under the hypotheses of Theorem 1.1, we have the following “classical” Carleson measure estimate and $L^{p}$ -solvability of the Dirichlet problem for $L$ for some $p$ [AHMT19, CHMT20, FP22]:

Lemma 3.1 ([FP22, Corollary 1.32]).

Suppose that $\Omega$ is a uniform domain in $\mathbb{R}^{n+1}$ , $n\geq 1$ , with $n$ -Ahlfors regular boundary and $L$ is a divergence form elliptic operator such that $\omega_{L}\in A_{\infty}(\sigma)$ . There exists a constant $C\geq 1$ such that if $u\in W^{1,2}_{\operatorname{loc}}(\Omega)\cap L^{\infty}(\Omega)$ is a solution to $Lu=0$ then

\displaystyle\sup_{r>0}\sup_{x\in\partial\Omega}\frac{1}{r^{n}}\iint_{B(x,r)\cap\Omega}|\nabla u(Y)|^{2}\delta(Y)\,dY\lesssim C\|u\|_{L^{\infty}(\Omega)}^{2}.

The constant $C$ depends on the structural constants and the $\omega_{L}\in A_{\infty}(\sigma)$ constants.

Lemma 3.2 ([HL18, Theorem 1.3]).

Suppose that $\Omega\subset\mathbb{R}^{n+1}$ , $n\geq 1$ , is a uniform domain with Ahlfors regular boundary and $L$ is a divergence form elliptic operator in $\Omega$ such that $\omega_{L}\in A_{\infty}(\sigma)$ . There exist $1<p<\infty$ and $C\geq 1$ such that for every $f\in L^{p}(\partial\Omega)$ there exists a solution $u$ to the boundary value problem

\displaystyle\begin{cases}Lu=0&\text{ in }\Omega,\\ \|Nu\|_{L^{p}(\partial\Omega)}\lesssim\|f\|_{L^{p}(\partial\Omega)},\\ \lim\limits_{Y\to x,\,\operatorname{n.t.}}u(Y)=f(x),&\text{ for $\sigma$-a.e.}x\in\partial\Omega.\end{cases}

Moreover, the solution is of the form

\displaystyle u(X)=u_{f}(X)=\int_{\partial\Omega}f(y)\,d\omega^{X}(y).

Using the $L^{p}$ result of Lemma 3.2 gives us the following non-tangential convergence result for $L^{\infty}$ data:

Lemma 3.3.

\displaystyle u_{f}(X)=\int_{\partial\Omega}f(y)\,d\omega^{X}(y)

converges non-tangentially to $f$ ; that is,

(3.4)

\lim\limits_{Y\to x,\,\operatorname{n.t.}}u_{f}(Y)=f(x),\quad\text{ for $\sigma$-a.e. }x\in\partial\Omega.

Equivalently,

(3.5)

\lim\limits_{Y\to x,\,\operatorname{n.t.}}u_{f}(Y)=f(x),\quad\text{ for $\omega^{Y}$-a.e. }x\in\partial\Omega.

Proof.

Let $p$ be as in Lemma 3.2. Let us note that $\omega^{Y}$ and $\omega^{Y^{\prime}}$ are mutually absolutely continuous for all $Y,Y^{\prime}\in\Omega$ by the Harnack chain condition and the Harnack inequality (see Lemma 2.29). Moreover, $\omega^{Y}$ is mutually absolutely continuous with respect to $\sigma$ by the $A_{\infty}(\sigma)$ condition [CF74, Lemma 5] and hence, (3.4) and (3.5) are equivalent. In addition, we have $f\in L^{\infty}(\partial\Omega,d\sigma)$ if and only if $f\in L^{\infty}(\partial\Omega,d\omega^{Y})$ for all $Y\in\Omega$ , with $\|f\|_{L^{\infty}(\partial\Omega,d\sigma)}=\|f\|_{L^{\infty}(\partial\Omega,d\omega^{Y})}$ .

If $\operatorname{diam}(\partial\Omega)<\infty$ , then $f\in L^{p}(\partial\Omega)$ and the claim follows from Lemma 3.2. Thus, we may assume that $\operatorname{diam}(\partial\Omega)=\infty$ . We show that the claim holds in any ball $B$ centered at $\partial\Omega$ . Let us fix $x_{0}\in\partial\Omega$ and $r>0$ . We write $f=g+h$ for $g=f\mathbbm{1}_{\Delta(x_{0},100r)}$ , where $\Delta(x_{0},100r)\coloneqq B(x_{0},100r)\cap\partial\Omega$ . By linearity, we have that $u_{f}(X)=u_{g}(X)+u_{h}(X)$ . Since $f\in L^{\infty}(\partial\Omega)$ , we know that $g\in L^{p}(\partial\Omega)$ and thus, Lemma 3.2 gives us

\displaystyle\lim\limits_{Y\to x,\,\operatorname{n.t.}}u_{g}(Y)=g(x)=f(x),\quad\text{ for $\sigma$-a.e. }x\in\partial\Omega\cap B(x_{0},r).

Therefore it suffices to show

(3.6)

\lim\limits_{Y\to x,\,\operatorname{n.t.}}u_{h}(Y)=0,\quad\text{ for $\sigma$-a.e. }x\in\partial\Omega\cap B(x_{0},r).

We now have

(3.7)

|u_{h}(X)|\leq\|f\|_{L^{\infty}(\partial\Omega)}\omega^{X}(\partial\Omega\setminus\Delta(x_{0},100r))

for any $X\in\Omega$ . Let $\phi\in C_{\operatorname{c}}(\Delta(x_{0},100r))$ be a non-negative function such that $\phi(x)\leq 1$ for every $x\in\partial\Omega$ and $\phi\equiv 1$ on $\Delta(x_{0},50r)$ . Then, since $\omega^{X}$ is a probability measure, we get

(3.8)

\displaystyle\omega^{X}(\partial\Omega\setminus\Delta(x_{0},100r))\leq 1-v(X)\coloneqq 1-\int_{\partial\Omega}\phi(y)\,d\omega^{X}(y).

Since $\phi$ is a compactly supported continuous function on $\partial\Omega$ , we know that $v$ is a continuous bounded solution to $Lv=0$ in $\overline{\Omega}$ such that $v=\phi\equiv 1$ on $\Delta(x_{0},50r)$ . In particular, we have

(3.9)

\displaystyle\lim_{Y\to x,\,\operatorname{n.t.}}v(Y)=1

for every $x\in\Delta(x_{0},50r)$ . Thus, combining (3.7) and (3.8) gives us $|u_{h}(X)|\leq\|f\|_{L^{\infty}(\partial\Omega)}\left(1-v(X)\right)$ for every $X\in\Omega$ and (3.6) follows then from (3.9). This completes the proof. ∎

3.2. A few important lemmas

In this subsection, we prove some lemmas that will be important for the proof of Theorem 1.1. The first three lemmas were inspired by the ideas in [DJK84].

Lemma 3.10.

Let $\Omega$ be a domain in $\mathbb{R}^{n+1}$ , $n\geq 1$ , and let $L$ be a not necessarily symmetric divergence form elliptic operator. Suppose that $u\in W^{1,2}_{\operatorname{loc}}(\Omega)$ is a weak solution to $Lu=0$ in $\Omega$ , that $\Omega^{\prime}\subset\Omega$ is a Wiener-regular domain which is compactly contained in $\Omega$ , and fix $X_{*}\in\Omega^{\prime}$ . If $n=1$ , assume in addition that $\Omega^{\prime}$ is a uniform domain with $n$ -Ahlfors regular boundary. Then $|\nabla u|^{2}G_{L,\Omega^{\prime}}(X_{*},\cdot)\in L^{1}(\Omega^{\prime})$ , and

(3.11)

\iint_{\Omega^{\prime}}|\nabla u(Y)|^{2}G_{L,\Omega^{\prime}}(X_{*},Y)\,dY\approx\int_{\partial\Omega^{\prime}}(u(y)-u(X_{*}))^{2}\,d\omega_{L,\Omega^{\prime}}^{X_{*}}(y),

where the implicit constants depend only on ellipticity.

Proof.

Throughout this argument we fix $X_{*}\in\Omega^{\prime}$ , we let $r=\operatorname{dist}(X_{*},\partial\Omega^{\prime})/8$ , and we write $G(Y)\coloneqq G_{\Omega^{\prime}}(X_{*},Y)$ . First, we show the finiteness of the integral in the left-hand side of (3.11). If $n=1$ , then it is trivial by (2.33) and the Meyers reverse Hölder estimate for gradients of solutions [Mey63], so suppose that $n\geq 2$ . We write

\iint_{\Omega^{\prime}}|\nabla u|^{2}G\,dY=\iint_{\Omega^{\prime}\backslash B(X_{*},r)}|\nabla u|^{2}G\,dY\\ +\iint_{B(X_{*},r)}|\nabla u|^{2}G\,dY\eqqcolon T_{1}+T_{2}.

Since $G\in L^{\infty}(\Omega^{\prime}\backslash B_{*}(X_{*},r))$ and $\nabla u\in L^{2}(\Omega^{\prime})$ , it is clear that $T_{1}<\infty$ . As for $T_{2}$ , for each $k\in\mathbb{N}_{0}$ , let $A_{k}\coloneqq B(X_{*},2^{-k}r)\backslash B(X_{*},2^{-k-1}r)$ , and using (2.33), Lemma 2.27, and Lemma 2.28, we see that

(3.12)

T_{2}=\sum_{k=0}^{\infty}\iint_{A_{k}}|\nabla u|^{2}G\,dY\lesssim\sum_{k=0}^{\infty}(2^{-k}r)^{-n+1}\iint_{B(X_{*},2^{-k}r)}|\nabla(u(Y)-u(X_{*}))|^{2}\,dY\\ \lesssim\sum_{k=0}^{\infty}(2^{-k}r)^{-n-1}\iint_{B(X_{*},2^{-k+1}r)}|u(Y)-u(X_{*})|^{2}\,dY\lesssim\sum_{k=0}^{\infty}2^{-2\alpha k}\|u\|_{L^{\infty}(B(X_{*},4r))}<\infty.

We turn to the proof of (3.11). By the ellipticity of $A$ and the product rule, we see that

(3.13)

\iint_{\Omega^{\prime}}|\nabla u|^{2}G\,dY\approx\iint_{\Omega^{\prime}}A\nabla u\cdot(\nabla u)G\,dY\\ =\iint_{\Omega^{\prime}}A\nabla u\cdot\nabla(uG)\,dY-\iint_{\Omega^{\prime}}A\nabla u\cdot(\nabla G)u\,dY\eqqcolon T_{3}+T_{4},

provided that the last two integrals are finite, which we now show. First, we claim that $|u||\nabla u||\nabla G|\in L^{1}(\Omega^{\prime})$ . As before, it is enough to show that $|u||\nabla u||\nabla G|\in L^{1}(B(X_{*},r))$ . Let $A_{k}$ as above, and for each $k\in\mathbb{N}_{0}$ , let $\{B_{k}^{j}\}_{j=1}^{J_{k}}$ be a cover of $A_{k}$ by balls centered on $A_{k}$ of radius $2^{-k-4}r$ with uniformly bounded overlap. Then $J_{k}\lesssim 1$ , and we have that

\iint_{B(X_{*},r)}|u||\nabla u||\nabla G|\,dY\leq C\sum_{k=0}^{\infty}\max_{j}\iint_{B_{k}^{j}}|\nabla u||\nabla G|\,dY\leq\left\{\begin{matrix}C\sum_{k=0}^{\infty}2^{-k\alpha},&\quad n\geq 2,\\[2.84526pt] C\sum_{k=0}^{\infty}k2^{-k\alpha},&\quad n=1,\end{matrix}\right.

where the last estimate follows from using Lemma 2.27 for both $u$ and $G$ , then (2.33) and Lemma 2.28. This proves the claim. By the product rule and the triangle inequality, we have also shown that $|\nabla u||\nabla(uG)|\in L^{1}(\Omega^{\prime})$ . Finally, by boundedness of $A$ we conclude that $A\nabla u\nabla(uG)$ and $A\nabla u(\nabla G)u$ belong to $L^{1}(\Omega^{\prime})$ , as desired.

The next step is to show that $T_{3}=0$ . For each $M\in\mathbb{N}$ large enough, let $\psi_{M}\in C^{\infty}(\Omega)$ satisfy $\psi_{M}\geq 0$ , $\psi_{M}\equiv 1$ in $\Omega\backslash B(X_{*},\frac{2}{M})$ , $\psi_{M}\equiv 0$ in $B(X_{*},\frac{1}{M})$ , and $|\nabla\psi_{M}|\lesssim M$ . We claim that

(3.14)

\iint_{\Omega^{\prime}}A\nabla u\nabla(uG)\psi_{M}\,dY\rightarrow 0,\qquad\text{as }M\rightarrow\infty.

Fix $M\in\mathbb{N}$ large enough. By the product rule and the fact that $uG\psi_{M}\in W_{0}^{1,2}(\Omega^{\prime})$ , since $Lu=0$ in $\Omega^{\prime}$ , we have that

(3.15)

\Big{|}\iint_{\Omega^{\prime}}A\nabla u\nabla(uG)\psi_{M}\,dY\Big{|}=\Big{|}\iint_{\Omega^{\prime}}A\nabla u(\nabla\psi_{M})uG\,dY\Big{|}\\ \leq CM\iint_{B(X_{*},\frac{2}{M})\backslash B(X_{*},\frac{1}{M})}|\nabla u|G\,dY\leq\left\{\begin{matrix}CM^{-\alpha},&\quad n\geq 2,\\[2.84526pt] CM^{-\alpha}\log M,&\quad n=1,\end{matrix}\right.

where in the last estimate we once again used (2.33), Lemma 2.27, and Lemma 2.28. Thus we have shown (3.14). Since it is also true that $A\nabla u\nabla(uG)\psi_{M}\rightarrow A\nabla u\nabla(uG)$ pointwise a.e. in $\Omega^{\prime}$ and since we have already proved that $|\nabla u||\nabla(uG)|\in L^{1}(\Omega^{\prime})$ , then by the Lebesgue Dominated Convergence Theorem we conclude that $T_{3}=0$ .

We proceed with the proof of (3.11) as follows:

T_{4}=-\frac{1}{2}\iint_{\Omega^{\prime}}A^{T}\nabla G\nabla(u^{2})\,dY=\frac{1}{2}\Big{(}\int_{\partial\Omega^{\prime}}u(y)^{2}\,d\omega_{L,\Omega^{\prime}}^{X_{*}}(y)-u(X_{*})^{2}\Big{)}\\ =\frac{1}{2}\int_{\partial\Omega^{\prime}}(u(y)-u(X_{*}))^{2}\,d\omega_{L,\Omega^{\prime}}^{X_{*}}(y),

where we used the Riesz formula (2.35), and in the last identity we used that $\omega_{L,\Omega^{\prime}}^{X_{*}}$ is a probability measure, and that $Lu=0$ in $\Omega^{\prime}$ so that

\displaystyle\int_{\partial\Omega^{\prime}}u(y)\,d\omega_{L,\Omega^{\prime}}^{X_{*}}(y)=u(X_{*}),

and hence

\displaystyle\int_{\partial\Omega^{\prime}}u(X_{*})u(y)\,d\omega_{L,\Omega^{\prime}}^{X_{*}}(y)=u(X_{*})\int_{\partial\Omega^{\prime}}u(y)\,d\omega_{L,\Omega^{\prime}}^{X_{*}}(y)=u(X_{*})^{2}.

This finishes the proof. ∎

The following result is a direct consequence of the maximum principle and the DeGiorgi–Nash–Moser theory (see for instance the proof of [ADFJM19, Proposition 3.2]).

Lemma 3.16.

Let $\Omega_{1}$ and $\Omega_{2}$ be Wiener regular domains such that $\Omega_{1}\subseteq\Omega_{2}$ . If $G_{i}(X,Y)$ is the Green function for $\Omega_{i}$ , $i=1,2$ then

\displaystyle G_{1}(X,Y)\leq G_{2}(X,Y)\quad\text{ for every }(X,Y)\in\Omega_{1}\times\Omega_{1}\setminus\{X=Y\}.

We will also use a corona-type decomposition of the elliptic measure $\omega_{L}$ [HLMN17, GMT18, AGMT22]. To formulate this decomposition, we recall the coherency and semicoherency of subcollections of $\mathbb{D}$ :

Definition 3.17.

We say that a subcollection $\mathcal{S}\subset\mathbb{D}$ is coherent if the following three conditions hold.

(a)

There exists a maximal element $Q(\mathcal{S})\in\mathcal{S}$ such that $Q\subset Q(\mathcal{S})$ for every $Q\in\mathcal{S}$ .
(b)

If $Q\in\mathcal{S}$ and $P\in\mathbb{D}$ is a cube such that $Q\subset P\subset Q(\mathcal{S})$ , then also $P\in\mathcal{S}$ .
(c)

If $Q\in\mathcal{S}$ , then either all children of $Q$ belong to $\mathcal{S}$ or none of them do.

If $\mathcal{S}$ satisfies only conditions (a) and (b), then we say that $\mathcal{S}$ is semicoherent.

We are ready to present the corona decomposition that we shall use in the sequel.

Lemma 3.18 ([AGMT22, Proposition 3.1]).

Suppose that $\Omega$ is a uniform domain and $L=-\mathop{\operatorname{div}}\nolimits A\nabla$ is a divergence form elliptic operator such that $\omega_{L}\in A_{\infty}(\sigma)$ . Then there exist constants, $C,M>1$ and a decomposition of the dyadic system $\mathbb{D}=\mathbb{D}(\partial\Omega)$ , with the following properties.

(i)

The dyadic grid breaks into a disjoint decomposition $\mathbb{D}=\mathcal{G}\cup\mathcal{B}$ , the good cubes and bad cubes respectively.
(ii)

The family $\mathcal{G}$ has a disjoint decomposition $\mathcal{G}=\cup\mathcal{S}$ where each $\mathcal{S}$ is a coherent stopping time regime with maximal cube $Q(\mathcal{S})$ .

(iii)

The maximal cubes and bad cubes satisfy a Carleson packing condition:

\displaystyle\sum_{\begin{subarray}{c}Q\in\mathcal{B}\\ Q\subseteq R\end{subarray}}\sigma(Q)+\sum_{\mathcal{S}\colon Q(\mathcal{S})\subseteq R}\sigma(Q(\mathcal{S}))\leq C\sigma(R)\quad\text{ for every }R\in\mathbb{D}.

(iv)

On each stopping time $\mathcal{S}$ , the elliptic measure ‘acts like surface measure’ in the sense that if $Q\in\mathcal{S}$ , then

\displaystyle M^{-1}\frac{\sigma(Q)}{\sigma(Q(\mathcal{S}))}\leq\frac{\omega^{X_{Q(\mathcal{S})}}(Q)}{\omega^{X_{Q(\mathcal{S})}}(Q(\mathcal{S}))}\leq M\frac{\sigma(Q)}{\sigma(Q(\mathcal{S}))}.

Proof.

Since $\omega_{L}\in A_{\infty}(\sigma)$ , then by Lemma 3.1 we have that the hypothesis (b) of [AGMT22, Proposition 3.1] is satisfied. By [AGMT22, Proposition 3.1], this yields a decomposition $\mathbb{D}=\mathcal{G}\cup\mathcal{B}$ like the one described above (in fact we have that $\mathcal{B}=\varnothing$ ), but with the caveat that the stopping time regimes need only be semicoherent. Then, by a standard mechanism (see for instance [DS93, pp.56–57], and [CHM, Remark 2.13]), we can modify the stopping time regimes so that they are coherent, while the rest of the properties are still satisfied. ∎

It is straightforward to check that for each semicoherent stopping time regime $\mathcal{S}$ , there exists a collection of pairwise disjoint cubes $\mathcal{F}_{\mathcal{S}}$ such that $\mathcal{S}=\mathbb{D}_{\mathcal{F}_{\mathcal{S}},Q(\mathcal{S})}$ . By Lemma 2.38 and the Harnack inequality, we get the following:

Lemma 3.19.

Suppose that $\Omega$ is a uniform domain in $\mathbb{R}^{n+1}$ , $n\geq 1$ , with $n$ -Ahlfors regular boundary, that $L=-\mathop{\operatorname{div}}\nolimits A\nabla$ is a divergence form elliptic operator such that $\omega_{L}\in A_{\infty}(\sigma)$ , and that $\mathcal{S}$ is a semicoherent stopping time regime satisfying the property in Lemma 3.18 (iv). Then

\displaystyle G_{L}(X_{Q},Y)\sigma(Q)\approx\delta(Y)\qquad\text{ for every }Q\in\mathcal{S}\text{ and }Y\in\Omega^{*}_{\mathcal{F}_{\mathcal{S}},Q},

where the implicit constants depend only on dimension, ellipticity, $\kappa$ , Harnack chain, corkscrew, and Ahlfors regularity constants, as well as the constant $M$ in Lemma 3.18 (iv).

Proof.

Fix $Q\in\mathcal{S}$ and³³3For relevant notation in this proof, see Section 2.3. $Y\in\Omega_{\mathcal{F}_{\mathcal{S}},Q}^{*}$ . By definition, there exists $P\in\mathbb{D}_{\mathcal{F}_{S},Q}$ so that $Y\in U_{P}^{*}$ , and thus $P\in\mathcal{S}$ , $P\subseteq Q$ , and $\delta(Y)\approx\ell(P)$ . Let $B_{P}^{\prime}\coloneqq B(x_{P},\eta a_{0}\ell(P))$ with $\eta\in(0,1)$ small, and let $X_{P}^{\prime}$ be the corkscrew point for $B_{P}^{\prime}$ . Define $B_{Q}^{\prime}$ and $X_{Q}^{\prime}$ analogously. We may guarantee that $X_{Q}\in\Omega\backslash 2B_{P}^{\prime}$ if $\eta$ is chosen small enough depending only on the corkscrew constant. Then, since $\delta(Y)\approx\delta(X_{P}^{\prime})$ , and $|Y-X_{P}^{\prime}|\lesssim\delta(Y)$ , by the Harnack inequality (for $L$ and $L^{*}$ ) and Harnack chains, Lemma 2.38, the doubling property of elliptic measure, and Ahlfors regularity, we have that

\frac{G_{L}(X_{Q},Y)}{\delta(Y)}\approx\frac{G_{L}(X_{Q},X_{P}^{\prime})}{\ell(P)}\approx\frac{\omega^{X_{Q}}(2B_{P}^{\prime})}{\ell(P)^{n}}\approx\frac{\omega^{X_{Q}^{\prime}}(B_{P}^{\prime})}{\sigma(P)}.

Now, since $X_{Q(\mathcal{S})}\in\Omega\backslash 2B_{Q}^{\prime}$ and $P,Q\in\mathcal{S}$ , then by Lemma 2.40, the doubling property of elliptic measure, Harnack chains, Harnack inequality, Bourgain’s estimate, and the property (iv) in Lemma 3.18, it follows that

\frac{\omega^{X_{Q}^{\prime}}(B_{P})}{\sigma(P)}\approx\frac{\omega^{X_{Q(\mathcal{S})}}(B_{P})}{\sigma(P)}\frac{1}{\omega^{X_{Q(\mathcal{S})}}(B_{Q}^{\prime})}\approx\frac{\omega^{X_{Q(\mathcal{S})}}(P)}{\sigma(P)}\frac{1}{\omega^{X_{Q(\mathcal{S})}}(Q)}\\ \approx\frac{1}{\sigma(Q(\mathcal{S}))}\frac{\sigma(Q(\mathcal{S}))}{\sigma(Q)}=\frac{1}{\sigma(Q)},

which completes the proof. ∎

4. Existence of $\varepsilon$ -approximators: proof of Theorem 1.1

In this section we prove Theorem 1.1. The key step in the proof is the construction of $\operatorname{BV}_{\operatorname{loc}}$ approximators since the existence of smooth approximators with more delicate properties follows then by using almost black box regularization arguments. More precisely, the heart of Theorem 1.1 is the following result:

Theorem 4.1.

Let $\Omega\subset\mathbb{R}^{n+1}$ , $n\geq 1$ , be a uniform domain with Ahflors regular boundary. Let $L=-\mathop{\operatorname{div}}\nolimits A\nabla$ be a (not necessarily symmetric) divergence form elliptic operator satisfying that $\omega_{L}\in A_{\infty}(\sigma)$ . Then, for any $\varepsilon\in(0,1)$ there exists a constant $C_{\varepsilon}$ such that if $u\in W^{1,2}_{\operatorname{loc}}(\Omega)\cap L^{\infty}(\Omega)$ is a solution to $Lu=0$ in $\Omega$ , then there exists $\Phi=\Phi^{\varepsilon}\in\operatorname{BV}_{\operatorname{loc}}(\Omega)$ satisfying

\displaystyle\|u-\Phi\|_{L^{\infty}(\Omega)}\leq\varepsilon\|u\|_{L^{\infty}(\Omega)}\quad\text{ and }\quad\sup_{x\in\partial\Omega,r>0}\frac{1}{r^{n}}\iint_{B(x,r)\cap\Omega}|\nabla\Phi(Y)|\,dY\leq C_{\varepsilon}\|u\|_{L^{\infty}(\Omega)},

where $C_{\varepsilon}$ depends on $\varepsilon$ , structural constants, ellipticity and the $\omega_{L}\in A_{\infty}(\sigma)$ constants.

As mentioned in the introduction (see the paragraph after Theorem 1.1), the proof of Theorem 4.1 is based on the construction of $\varepsilon$ -approximators in [HMM16, Theorem 1.3], with important differences. We will point out when we have to diverge from the strategy in [HMM16].

4.1. Set-up for the proof of Theorem 4.1

We start by making some preliminary considerations and fixing some notation and terminology. Let $u$ be a bounded solution to $Lu=0$ , and let $\varepsilon\in(0,1)$ . Without loss of generality, we may assume that $\|u\|_{L^{\infty}(\Omega)}=1$ , since otherwise we have $u=0$ or we can replace $u$ with $u/\|u\|_{L^{\infty}}$ . We fix a stopping time regime $\mathcal{S}$ from Lemma 3.18, with maximal cube $Q(\mathcal{S})$ . Following [Gar07, HMM16], we label each cube $Q\in\mathbb{D}$ depending on how much the the function $u$ oscillates in the corresponding fattened Whitney region $U_{Q}^{*}$ . To be more precise, we say that

\displaystyle Q\in\mathbb{D}\text{ is {\bf red} if }\quad\sup_{X,Y\in U_{Q}^{*}}|u(X)-u(Y)|\geq\frac{\varepsilon}{1000},

and

\displaystyle Q\in\mathbb{D}\text{ is {\bf blue} if }\quad\sup_{X,Y\in U_{Q}^{*}}|u(X)-u(Y)|<\frac{\varepsilon}{1000}.

We note that these conditions differ from the, ones in [HMM16] in two ways: in [HMM16], the oscillation threshold was $\varepsilon/10$ and the level of oscillation was measured in the unfattened region $U_{Q}$ . With our conditions we have a little bit more control on the blue cubes, which will be useful for us later. In [HMM16], it was enough to use these two labels but for our analysis it is important to take into consideration the minimal cubes of $\mathcal{S}$ , which we label as follows:

\displaystyle Q\in\mathcal{S}\text{ is {\bf yellow} if }\quad Q\text{ has a child $Q^{\prime}$ such that $Q^{\prime}\not\in\mathcal{S}$}.

Recall that $Q^{\prime}$ is a child of $Q$ if $Q^{\prime}\subset Q$ and $\ell(Q^{\prime})=\ell(Q)/2$ . Notice that yellow cubes are not a separate collection from the red and blue cubes but each yellow cube is also red or blue. We denote the collections of red and yellow cubes by

\mathcal{R}\coloneqq\{Q\in\mathbb{D}\colon Q\text{ is red}\},\qquad\mathcal{Y}=\mathcal{Y}(\mathcal{S})\coloneqq\{Q\in\mathcal{S}\colon Q\text{ is yellow}\}.

The rough idea of the construction of the $\varepsilon$ -approximator $\Phi$ of $u$ is to first construct $\Phi$ inside a Carleson box $T_{Q_{0}}$ for a fixed cube $Q_{0}$ and then use a “local to global”-type argument to define a global approximator. Working in a Carleson box allows us to reduce many of the challenging estimates to working with just one stopping time regime $\mathcal{S}$ given by Lemma 3.18. We then set $\Phi\equiv u$ in the regions $U_{Q}$ such that $Q\in\mathcal{R}\cup\mathcal{Y}$ and break up the blue cubes into smaller stopping time regimes where $u$ does not vary by more than $\varepsilon/100$ . In these new regimes, we set $\Phi=u(X_{0})$ in the union of $U_{Q}$ , where $X_{0}$ is any point in the union. The $L^{\infty}$ approximation property follows then just from the way we defined $\Phi$ but verifying the $L^{1}$ -type Carleson measure estimate for $\Phi$ is more challenging. For this, one of the key steps is to show that the collections of red, yellow and maximal cubes from the new stopping time regimes satisfy Carleson packing conditions. For the first two collections, this follows in a straightforward way:

Lemma 4.2.

The collections $\mathcal{R}$ and $\bigcup_{\mathcal{S}}\mathcal{Y}(\mathcal{S})$ satisfy the following Carleson packing conditions: for any $P\in\mathbb{D}$ we have

\displaystyle\sum_{Q\in\mathcal{R},Q\subset P}\sigma(Q)\lesssim\frac{1}{\varepsilon^{2}}\sigma(P)

and

\displaystyle\sum_{\mathcal{S}}\sum_{Q\in\mathcal{Y}(\mathcal{S}),Q\subset P}\sigma(Q)\leq(C+1)\sigma(P).

The constant $C$ is the same constant as $C$ in Lemma 3.18.

Proof.

The packing condition for $\bigcup_{\mathcal{S}}\mathcal{Y}(\mathcal{S})$ follows from the definition and the facts that the stopping time regimes are coherent and their maximal cubes satisfy a Carleson packing condition. Indeed, each yellow cube $Q\in\mathcal{Y}(\mathcal{S})$ is contained in $Q(\mathcal{S})$ and by the coherency of $\mathcal{S}$ , no yellow cube can contain a smaller yellow cube $\widetilde{Q}\in\mathcal{Y}(\mathcal{S})$ . Thus, the cubes in $\mathcal{Y}(\mathcal{S})$ are disjoint and we get

\sum_{\mathcal{S}}\sum_{Q\in\mathcal{Y}(\mathcal{S}),Q\subset P}\sigma(Q)\leq\sum_{\begin{subarray}{c}\mathcal{S}:Q\in\mathcal{Y}(\mathcal{S}),Q\subset P,\\ Q(\mathcal{S})\subset P\end{subarray}}\sigma(Q)+\sum_{\begin{subarray}{c}\mathcal{S}:Q\in\mathcal{Y}(\mathcal{S}),Q\subset P,\\ P\in\mathcal{S}\end{subarray}}\sigma(Q)\\ \leq\sum_{\begin{subarray}{c}\mathcal{S}:Q(\mathcal{S})\subset P\end{subarray}}\sigma(Q(\mathcal{S}))+\sum_{\begin{subarray}{c}\mathcal{S}:P\in\mathcal{S}\end{subarray}}\sigma(P).

Since the collection $\{Q(\mathcal{S})\}_{\mathcal{S}}$ satisfies a Carleson packing condition by Lemma 3.18, we know that the first sum on the right-hand side is bounded by $C\sigma(P)$ . In addition, since the stopping time regimes ${\bf S}$ are disjoint, we have $P\in\mathcal{S}$ for at most one stopping time regime $\mathcal{S}$ . Thus, the second sum on the right-hand side is bounded by $\sigma(P)$ . The desired bound follows from combining these two estimates.

Let us then prove the Carleson packing condition for $\mathcal{R}$ . If $Q\in\mathcal{R}$ , then there exist points $X_{1},X_{2}\in U_{Q}^{*}$ such that $|u(X_{1})-u(X_{2})|>\varepsilon$ . Since $X_{1},X_{2}\in U_{Q}^{*}$ , there exist Whitney cubes $I_{1},I_{2}\subset U_{Q}^{*}$ such that $X_{1}\in(1+\tau)I_{1}$ , $X_{2}\in(1+\tau)I_{2}$ and $\ell(Q)\approx\ell(I_{1})\approx\ell(I_{2})\approx\delta(X_{1})\approx\delta(X_{2})$ . In particular, we have $\operatorname{dist}(X_{1},\partial U_{Q}^{**})\approx\operatorname{dist}(X_{2},\partial U_{Q}^{**})\approx\ell(Q)$ for the twice-fattened region $U_{Q}^{**}$ (see Section 2.3), where the implicit constants depend on the dilation parameter $\tau$ . Thus, since $U_{Q}^{**}$ satisfies the Harnack chain condition by Lemma 2.18, there exists a uniformly bounded number of balls $B_{1},B_{2},\ldots,B_{N}$ and points $Y_{1},Y_{2},\ldots,Y_{N-1}$ such that

$\bullet$

$X_{1}\in B_{1}$ , $X_{2}\in B_{N}$ and $Y_{i}\in B_{i}\cap B_{i+1}$ for every $i=1,2,\ldots,N-1$ ,
$\bullet$

for a constant $\lambda>1$ (depending on $\tau$ ), we have $2\lambda B_{i}\subset U_{Q}^{**}$ for every $i=1,2,\ldots,N$ , and
$\bullet$

$|2\lambda B_{i}|\approx|U_{Q}^{**}|\approx\ell(Q)^{n+1}$ for every $i=1,2,\ldots,N$ .

These properties combined with the triangle inequality, the local Hölder continuity (that is, Lemma 2.28) and the Poincaré inequality applied for solution $v$ , $v(X)\coloneqq u(X)-\fint_{U_{Q}^{**}}u(Z)\,dZ$ , then give us

\varepsilon\leq|v(X_{1})-v(X_{2})|\leq|v(X_{1})-v(Y_{1})|+\sum_{i=1}^{N-2}|v(Y_{i})-v(Y_{i+1})|+|v(Y_{N})-v(X_{2})|\\ \lesssim\mathchoice{{\vbox{\hbox{$\textstyle\raisebox{-8.61108pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-7.09656pt}}{{\vbox{\hbox{$\scriptstyle\raisebox{-6.02777pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-6.09258pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-4.29826pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-3.54826pt}}\!\iint_{2\lambda B_{1}}|v|\,dX+\sum_{i=1}^{N-2}\mathchoice{{\vbox{\hbox{$\textstyle\raisebox{-8.61108pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-7.09656pt}}{{\vbox{\hbox{$\scriptstyle\raisebox{-6.02777pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-6.09258pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-4.29826pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-3.54826pt}}\!\iint_{2\lambda B_{i+1}}|v|\,dX+\mathchoice{{\vbox{\hbox{$\textstyle\raisebox{-8.61108pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-7.09656pt}}{{\vbox{\hbox{$\scriptstyle\raisebox{-6.02777pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-6.09258pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-4.29826pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-3.54826pt}}\!\iint_{2\lambda B_{N}}|v|\,dX\lesssim\mathchoice{{\vbox{\hbox{$\textstyle\raisebox{-8.61108pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-7.09656pt}}{{\vbox{\hbox{$\scriptstyle\raisebox{-6.02777pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-6.09258pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-4.29826pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-3.54826pt}}\!\iint_{U_{Q}^{**}}|v|\,dX\lesssim\ell(Q)\Big{(}\mathchoice{{\vbox{\hbox{$\textstyle\raisebox{-8.61108pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-7.09656pt}}{{\vbox{\hbox{$\scriptstyle\raisebox{-6.02777pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-6.09258pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-4.29826pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-3.54826pt}}\!\iint_{U_{Q}^{**}}|\nabla u|^{2}\,dX\Big{)}^{1/2}.

Thus, since $\delta(X)\approx\ell(Q)$ for every $X\in U_{Q}^{**}$ , we have

\varepsilon^{2}\sigma(Q)\approx\ell(Q)^{n}\varepsilon^{2}\lesssim\ell(Q)^{n}\ell(Q)^{2}\mathchoice{{\vbox{\hbox{$\textstyle\raisebox{-8.61108pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-7.09656pt}}{{\vbox{\hbox{$\scriptstyle\raisebox{-6.02777pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-6.09258pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-4.29826pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-3.54826pt}}\!\iint_{U_{Q}^{**}}|\nabla u(X)|^{2}\,dX\\ \lesssim\iint_{U_{Q}^{**}}|\nabla u(X)|^{2}\ell(Q)\,dX\approx\iint_{U_{Q}^{**}}|\nabla u(X)|^{2}\delta(X)\,dX.

By construction, we know that the regions $U_{P}^{**}$ have bounded overlaps and any twice-fattened Carleson box $T_{P}^{**}$ satisfies $T_{P}^{**}\subset B(x_{P},R_{P})\cap\Omega$ the center $x_{P}$ of $P$ and for a radius $R_{P}\approx\ell(P)$ . Hence, for any $Q_{0}\in\mathbb{D}$ , these facts, the previous estimate and the Carleson measure estimate (3.1) give us

\sum_{\begin{subarray}{c}P\in\mathcal{R}\\ P\subset Q_{0}\end{subarray}}\sigma(P)\lesssim\varepsilon^{-2}\sum_{\begin{subarray}{c}P\in\mathcal{R}\\ P\subset Q_{0}\end{subarray}}\iint_{U_{P}^{**}}|\nabla u(X)|^{2}\delta(X)\,dX\\ \lesssim\varepsilon^{-2}\iint_{T^{**}_{Q_{0}}}|\nabla u(X)|^{2}\delta(X)\,dX\lesssim\varepsilon^{-2}\ell(Q_{0})^{n}\approx\varepsilon^{-2}\sigma(Q_{0}),

which proves the claim. ∎

4.2. A stopping time decomposition for the family of blue cubes

We now move to decomposing the collection of blue cubes into more manageable subcollections using a stopping time procedure. A similar idea is utilized in [HMM16, p. 2360] but, due to our geometry, we use different stopping time conditions and different analysis of the subcollections. Set $\mathcal{L}=\mathcal{L}(\mathcal{S})$ to be the collection of blue cubes in $\mathcal{S}$ . We first take the largest blue cube in $\mathcal{S}$ with respect to side length (if there is more than one such cube, we just pick one) and denote this cube by $Q({\bf S}_{1})$ . The cube $Q({\bf S}_{1})$ will be the maximal cube in our first refined stopping time regime ${\bf S}_{1}$ . We let $\mathcal{F}_{{\bf S}_{1}}$ be the collection of cubes $Q\in\mathcal{S}\cap\mathbb{D}_{Q({\bf S}_{1})}\setminus\{Q({\bf S}_{1})\}$ such that $Q$ is a maximal cube with respect to having one of the following three properties:

(1)

$Q$ or one of its siblings⁴⁴4Here $Q^{\prime}$ is a sibling of $Q$ if $Q,Q^{\prime}\in\mathbb{D}_{k}$ and $Q,Q^{\prime}\subset Q^{*}\in\mathbb{D}_{k-1}$ , that is, $Q$ and $Q^{\prime}$ have the same parent. For simplicity, we consider $Q$ to be its own sibling. is red.
(2)

$Q$ and all of its siblings are blue, but for some $Q^{\prime}$ that is either $Q$ or a sibling of $Q$ it holds that

$|u(X_{Q^{\prime}})-u(X_{Q({\bf S}_{1})})|>\varepsilon/100.$
(3)

$Q$ is yellow.

Recall that for every $P\in\mathbb{D}$ , the point $X_{P}$ is a corkscrew point relative to the center $x_{P}$ at scale $10^{-5}a_{0}\ell(P)$ , as defined in Section 2.3. We now set ${\bf S}_{1}\coloneqq\mathbb{D}_{\mathcal{F}_{{\bf S}_{1}},Q({\bf S}_{1})}$ . By construction, ${\bf S}_{1}$ is a coherent stopping time regime in the sense of Definition 3.17. Since the cubes in ${\bf S}_{1}$ are blue and none of them satisfy the stopping condition (2) above, we know that

(4.3)

|u(X)-u(X_{Q({\bf S}_{1})})|\leq\varepsilon/50\quad\text{ for every }X\in\Omega_{\mathcal{F}_{{\bf S}_{1}},Q({\bf S}_{1})}^{*}.

We now express $\mathcal{F}_{{\bf S}_{1}}$ as a union of three collections,

\displaystyle\mathcal{F}_{{\bf S}_{1}}=\mathcal{F}_{{\bf S}_{1}}^{\operatorname{R}}\cup\mathcal{F}_{{\bf S}_{1}}^{\operatorname{SB}}\cup\mathcal{F}_{{\bf S}_{1}}^{\operatorname{Y}},

where $\mathcal{F}_{{\bf S}_{1}}^{\operatorname{R}}$ contains the cubes for which (1) holds, $\mathcal{F}_{{\bf S}_{1}}^{\operatorname{SB}}$ contains the cubes for which (2) holds and $\mathcal{F}_{{\bf S}_{1}}^{\operatorname{Y}}$ contains the cubes for which (3) holds. The superscripts stand for “red”, “stopping blue” and “yellow”. We note that the collections $\mathcal{F}_{{\bf S}_{1}}^{\operatorname{R}}$ and $\mathcal{F}_{{\bf S}_{1}}^{\operatorname{SB}}$ are disjoint but $\mathcal{F}_{{\bf S}_{1}}^{\operatorname{Y}}$ may overlap with both of them.

We now continue this way: we let $Q({\bf S}_{2})$ be the largest blue cube in $\mathcal{S}\setminus{\bf S}_{1}$ with respect to side length, we extract the collection of maximal stopping cubes $\mathcal{F}_{Q({\bf S}_{2})}$ (with an updated stopping condition (2)), we define the coherent stopping regime ${\bf S}_{2}$ and the collections $\mathcal{F}_{{\bf S}_{2}}^{\operatorname{R}}$ , $\mathcal{F}_{{\bf S}_{2}}^{\operatorname{SB}}$ and $\mathcal{F}_{{\bf S}_{2}}^{\operatorname{Y}}$ , choose the largest blue cube $Q({\bf S}_{3})$ in $\mathcal{S}\setminus{\bf S}_{1}\cup{\bf S}_{2}$ , and so on. Since each ${\bf S}_{i}$ contains at least the cube $Q({\bf S}_{i})$ , we know that this procedure exhausts $\mathcal{L}$ and gives us a disjoint decomposition $\mathcal{L}=\cup_{j}{\bf S}_{j}$ where each ${\bf S}_{j}$ is a coherent stopping time regime. Just like (4.3), we have the oscillation estimate

(4.4)

|u(X)-u(X_{Q({\bf S}_{j})})|\leq\varepsilon/50\quad\text{ for every }X\in\Omega_{\mathcal{F}_{{\bf S}_{j}},Q({\bf S}_{j})}^{*}

for every $j$ . We also get the collections $\mathcal{F}_{{\bf S}_{j}}^{\operatorname{R}}$ , $\mathcal{F}_{{\bf S}_{j}}^{\operatorname{SB}}$ and $\mathcal{F}_{{\bf S}_{j}}^{\operatorname{Y}}$ for each $j$ .

Our next goal is to show that the maximal cubes $\{Q({\bf S}_{j})\}_{j}$ satisfy a Carleson packing condition. This goal is an analog of [HMM16, Lemma 5.16] but since the proof of this lemma is based on the use of an “ $N\lesssim S$ ” estimate in sawtooth regions (which is possible in the presence of uniform rectifiability of $\partial\Omega)$ , this is the part where we significantly depart from [HMM16]. Following an idea in [DS91], we let $\lambda\in(0,10^{-10})$ be a small parameter (to be chosen) and break the stopping times into four groups. We say that ${\bf S}_{j}$ is of

$\bullet$

Type 1 (T1) if $\sigma(Q({\bf S}_{j})\setminus\cup_{Q\in\mathcal{F}_{{\bf S}_{j}}}Q)\geq\lambda\sigma(Q({\bf S}_{j}))$ .
$\bullet$

Type 2 (T2) if $\sigma(\cup_{Q\in\mathcal{F}_{{\bf S}_{j}}^{\operatorname{R}}}Q)\geq\lambda\sigma(Q({\bf S}_{j}))$ .
$\bullet$

Type 3 (T3) if $\sigma(\cup_{Q\in\mathcal{F}_{{\bf S}_{j}}^{\operatorname{Y}}}Q)\geq\lambda\sigma(Q({\bf S}_{j}))$ .
$\bullet$

Type 4 (T4) if ${\bf S}_{j}$ is not type 1, 2 or 3.

The Carleson packing condition for the cubes $Q({\bf S}_{j})$ follows in a straightforward way when ${\bf S}_{j}$ is of Type 1, 2 or 3:

Lemma 4.5.

We have the following Carleson packing conditions for the maximal cubes of the subregimes ${\bf S}_{j}$ in the decomposition $\mathcal{L}(\mathcal{S})=\cup_{j}{\bf S}_{j}$ : for any $P\in\mathbb{D}$ , we have

	$\displaystyle\sum_{\begin{subarray}{c}j\colon{\bf S}_{j}\text{ is T1}\\ Q({\bf S}_{j})\subset P\end{subarray}}\sigma(Q({\bf S}_{j}))+\sum_{\begin{subarray}{c}j\colon{\bf S}_{j}\text{ is T3}\\ Q({\bf S}_{j})\subset P\end{subarray}}\sigma(Q({\bf S}_{j}))$	$\displaystyle\lesssim\frac{1}{\lambda}\sigma(P),\qquad\text{and,}$
	$\displaystyle\sum_{\begin{subarray}{c}j\colon{\bf S}_{j}\text{ is T2}\\ Q({\bf S}_{j})\subset P\end{subarray}}\sigma(Q({\bf S}_{j}))$	$\displaystyle\lesssim\frac{1}{\lambda\varepsilon^{2}}\sigma(P).$

Proof.

For the regimes of Type 1, we first notice that if $Q({\bf S}_{i})\subsetneq Q({\bf S}_{j})$ , then $Q({\bf S}_{i})\subset Q$ for some $Q\in\mathcal{F}_{{\bf S}_{j}}$ . In particular, the sets $Q({\bf S}_{j})\setminus\cup_{Q\in\mathcal{F}_{{\bf S}_{j}}}Q$ are pairwise disjoint. Thus, by the definition of Type 1, we get

\displaystyle\sum_{\begin{subarray}{c}j\colon{\bf S}_{j}\text{ is T1}\\ Q({\bf S}_{j})\subset P\end{subarray}}\sigma(Q({\bf S}_{j}))\leq\frac{1}{\lambda}\sum_{\begin{subarray}{c}j\colon{\bf S}_{j}\text{ is T1}\\ Q({\bf S}_{j})\subset P\end{subarray}}\sigma(Q({\bf S}_{j})\setminus\cup_{Q\in\mathcal{F}_{{\bf S}_{j}}}Q)\leq\frac{1}{\lambda}\sigma(P).

For Type 3 regimes, since the cubes $Q\in\mathcal{F}_{{\bf S}_{j}}^{\operatorname{Y}}\subset\mathcal{Y}(\mathcal{S})$ are yellow cubes in $\mathcal{S}$ , they are disjoint. Thus, the claim for the regimes of Type 3 follows immediately from definition.

For the regimes of Type 2, we recall that if $Q\in\mathcal{F}_{{\bf S}_{j}}^{\operatorname{R}}$ , then $Q$ is red or one of its siblings is red. In particular, each $Q\in\mathcal{F}_{{\bf S}_{j}}^{\operatorname{R}}$ has approximately the same measure as some red sibling $R_{Q}\subset Q({\bf S}_{j})$ of $Q$ . If there is more than one red sibling, we just choose one of them for each $Q\in\mathcal{F}_{{\bf S}_{j}}^{\operatorname{R}}$ . On the other hand, since each cube has only a uniformly bounded number of siblings, for each $Q\in\bigcup_{j}\mathcal{F}_{{\bf S}_{j}}^{\operatorname{R}}$ we can have $R_{Q}=R_{Q^{\prime}}$ only for a uniformly bounded number of cubes $Q^{\prime}$ . Thus, we get

\sum_{\begin{subarray}{c}j\colon{\bf S}_{j}\text{ is T2}\\ Q({\bf S}_{j})\subset P\end{subarray}}\sigma(Q({\bf S}_{j}))\leq\frac{1}{\lambda}\sum_{\begin{subarray}{c}j\colon{\bf S}_{j}\text{ is T2}\\ Q({\bf S}_{j})\subset P\end{subarray}}\sum_{Q\in\mathcal{F}_{{\bf S}_{j}}^{\operatorname{R}}}\sigma(Q)\approx\frac{1}{\lambda}\sum_{\begin{subarray}{c}j\colon{\bf S}_{j}\text{ is T2}\\ Q({\bf S}_{j})\subset P\end{subarray}}\sum_{Q\in\mathcal{F}_{{\bf S}_{j}}^{\operatorname{R}}}\sigma(R_{Q})\\ \lesssim\frac{1}{\lambda}\sum_{R\in\mathcal{R},R\subset P}\sigma(R)\lesssim\frac{1}{\lambda\varepsilon^{2}}\sigma(P),

where we used Lemma 4.2 in the final estimate. ∎

By Lemma 4.5, to show the Carleson packing condition for the collection $\{Q({\bf S}_{j})\}_{j}$ it remains only to consider the regimes ${\bf S}_{j}$ of Type 4.

Lemma 4.6.

There exists $\lambda_{0}>0$ depending only on structural constants, ellipticity, and the $\omega_{L,\Omega}\in A_{\infty}(\sigma)$ constants, such that for any $\lambda\in(0,\lambda_{0})$ , there exist constants $C_{1},C_{2}\geq 1$ depending only on structural constants, ellipticity, and the $\omega_{L,\Omega}\in A_{\infty}(\sigma)$ constants (and independent of $\varepsilon$ , $\lambda$ , and $\mathcal{S}$ ), so that

\displaystyle\sum_{\begin{subarray}{c}j\colon{\bf S}_{j}\text{ is T4}\\ Q({\bf S}_{j})\subset P\end{subarray}}\sigma(Q({\bf S}_{j}))

\displaystyle\leq\frac{C_{1}}{\varepsilon^{2}}\sigma(P).

for every $P\in\mathbb{D}$ . In particular, we have

\displaystyle\sum_{j\colon Q({\bf S}_{j})\subset P}\sigma(Q({\bf S}_{j}))

\displaystyle\leq\frac{C_{2}}{\lambda\varepsilon^{2}}\sigma(P).

Proving Lemma 4.6 is much more delicate than proving Lemma 4.5 and we do this in several steps. The key idea is to reduce the proof to proving estimates for which we can use lemmas from Section 3.

Let us fix a regime ${\bf S}_{j}$ that is of Type 4. Since ${\bf S}_{j}$ is not of Type 1, 2 or 3, we know that, roughly speaking, at the “bottom” of the sawtooth domain $\Omega_{\mathcal{F}_{{\bf S}_{j}},Q({\bf S}_{j})}$ there is a large region where $u$ has some (uniform) oscillation from the value $u(X_{Q({\bf S}_{j})})$ . Let us be more precise. Since ${\bf S}_{j}$ is not of Type 1, we know that $\sigma(Q({\bf S}_{j})\setminus\cup_{Q\in\mathcal{F}_{{\bf S}_{j}}}Q)<\lambda\sigma(Q({\bf S}_{j})$ , and since ${\bf S}_{j}$ is not of Type 2 or 3, we have

\displaystyle\sigma\big{(}(\cup_{Q\in\mathcal{F}_{{\bf S}_{j}}^{\operatorname{R}}}Q)\cup(\cup_{Q\in\mathcal{F}_{{\bf S}_{j}}^{\operatorname{Y}}}Q)\big{)}<2\lambda\sigma(Q({\bf S}_{j})).

Thus, since $\mathcal{F}_{{\bf S}_{j}}=\mathcal{F}_{{\bf S}_{j}}^{\operatorname{R}}\cup\mathcal{F}_{{\bf S}_{j}}^{\operatorname{SB}}\cup\mathcal{F}_{{\bf S}_{j}}^{\operatorname{Y}}$ , it holds that

\displaystyle\sigma\big{(}\cup_{Q\in\mathcal{F}_{{\bf S}_{j}}^{\operatorname{SB}}}Q\big{)}\geq(1-3\lambda)\sigma(Q({\bf S}_{j})).

Let $N=N(j)$ be so large that the subcollection

\displaystyle\mathcal{F}_{N,j}^{\operatorname{SB}}\coloneqq\{Q\in\mathcal{F}_{{\bf S}_{j}}^{\operatorname{SB}}\colon\ell(Q)>2^{-N}\ell(Q({\bf S}_{j}))\}

satisfies

(4.7)

\sigma\big{(}\cup_{Q\in\mathcal{F}_{N,j}^{\operatorname{SB}}}Q\big{)}\geq(1-4\lambda)\sigma(Q({\bf S}_{j})).

Our estimates will not depend on $N$ but we work with the subcollection $\mathcal{F}_{N,j}^{\operatorname{SB}}$ to avoid dealing with estimates on the boundary. We let $\mathcal{F}_{N,j}$ denote the collection of maximal cubes in the collection $\mathcal{F}_{{\bf S}_{j}}\cup\mathbb{D}_{N,Q({\bf S}_{j})}$ , where $\mathbb{D}_{N,Q}=\{Q^{\prime}\in\mathbb{D}_{Q}\colon\ell(Q^{\prime})=2^{-N}\ell(Q)\}$ as earlier. We note that $\Omega_{\mathcal{F}_{N,j},Q({\bf S}_{j})}\subseteq\Omega_{\mathcal{F}_{{\bf S}_{j}},Q({\bf S}_{j})}$ . Then

\displaystyle\mathcal{F}_{N,j}=\mathcal{F}_{N,j}^{\operatorname{SB}}\cup\mathcal{F}_{N,j}^{\operatorname{O}},\ \text{ where }\ \mathcal{F}_{N,j}^{\operatorname{O}}\coloneqq\mathcal{F}_{N,j}\setminus\mathcal{F}_{N,j}^{\operatorname{SB}},

where the superscript O in $\mathcal{F}_{N,j}^{\operatorname{O}}$ stands for “other” cubes. By (4.7) we have

(4.8)

\displaystyle\sigma\big{(}\cup_{Q\in\mathcal{F}_{N,j}^{\operatorname{O}}}Q\big{)}=\sum_{Q\in\mathcal{F}_{N,j}^{\operatorname{O}}}\sigma(Q)\leq(4\lambda)\sigma(Q({\bf S}_{j})).

With the notation above, we can formulate our key estimate for the proof of Lemma 4.6:

Lemma 4.9.

Suppose that ${\bf S}_{j}$ is a stopping time regime of Type 4. There exists $\lambda_{0}>0$ depending only on structural constants, ellipticity, and the $\omega_{L,\Omega}\in A_{\infty}(\sigma)$ constants, such that for any $\lambda\in(0,\lambda_{0})$ , there exists a constant $C_{3}\geq 1$ depending only on structural constants, ellipticity, and the $\omega_{L,\Omega}\in A_{\infty}(\sigma)$ constants (and independent of $\varepsilon$ , $\lambda$ , $N$ , $j$ , and $\mathcal{S}$ ), so that the following estimate holds:

(4.10)

\displaystyle\sigma(Q({\bf S}_{j}))\leq\frac{C_{3}}{\varepsilon^{2}}\iint_{\Omega_{\mathcal{F}_{N,j},Q({\bf S}_{j})}}|\nabla u(Y)|^{2}\delta(Y)\,dY.

Taking Lemma 4.9 for granted momentarily, we can prove Lemma 4.6 in a straightforward way:

Proof of Lemma 4.6.

Fix $P\in\mathbb{D}$ . Then we have

\sum_{\begin{subarray}{c}j\colon{\bf S}_{j}\text{ is T4}\\ Q({\bf S}_{j})\subset P\end{subarray}}\sigma(Q({\bf S}_{j}))\leq\frac{C_{3}}{\varepsilon^{2}}\sum_{\begin{subarray}{c}j\colon{\bf S}_{j}\text{ is T4}\\ Q({\bf S}_{j})\subset P\end{subarray}}\iint_{\Omega_{\mathcal{F}_{N,j},Q({\bf S}_{j})}}|\nabla u(Y)|^{2}\delta(Y)\,dY\\ \lesssim\frac{C_{3}}{\varepsilon^{2}}\iint_{T_{P}}|\nabla u(Y)|^{2}\delta(Y)\,dY\lesssim\frac{C_{3}}{\varepsilon^{2}}\sigma(P),

where we used Lemma 4.9, the fact that the bounded overlap of the regions $U_{Q}$ and the disjointness of the collections ${\bf S}_{j}$ imply that the regions $\Omega_{\mathcal{F}_{N(j),j},Q({\bf S}_{j})}$ have bounded overlaps, the fact that $T_{P}\subset B(x_{P},C\ell(P))$ and Lemma 3.1. The rest of the claim follows now from Lemma 4.5. ∎

4.3. Proof of Lemma 4.9: A high oscillation estimate

Let us then start processing the estimate (4.10). Let ${\bf S}_{j}$ be a fixed stopping time regime of Type 4. To relax the notation, we denote

	$\displaystyle{\bf S}_{}\coloneqq\mathbb{D}_{\mathcal{F}_{N,j},Q({\bf S}_{j})}\subseteq{\bf S}_{j},\qquad Q({\bf S}_{})\coloneqq Q({\bf S}_{j}),\qquad X_{*}\coloneqq X_{Q({\bf S}_{j})},$
	$\displaystyle\Omega_{*}\coloneqq\Omega_{\mathcal{F}_{N,j},Q({\bf S}_{j})},\qquad\mathcal{F}\coloneqq\mathcal{F}_{N,j},\qquad\mathcal{F}^{\operatorname{SB}}\coloneqq\mathcal{F}_{N,j}^{\operatorname{SB}},\quad\text{ and}\quad\mathcal{F}^{\operatorname{O}}\coloneqq\mathcal{F}_{N,j}^{\operatorname{O}}.$

Recall that $Q({\bf S}_{*})\in\mathcal{S}$ and that $\mathcal{S}$ is a coherent stopping time regime satisfying the property (iv) in Lemma 3.18. Thus, by using Lemma 3.19, Lemma 3.16 and Lemma 3.10 in this order, we get

(4.11)

\iint_{\Omega_{*}}|\nabla u(Y)|^{2}\delta(Y)\,dY\approx\sigma(Q({\bf S}_{*}))\iint_{\Omega_{*}}|\nabla u(Y)|^{2}G_{L,\Omega}(X_{*},Y)\,dY\\ \geq\sigma(Q({\bf S}_{*}))\iint_{\Omega_{*}}|\nabla u(Y)|^{2}G_{L,\Omega_{*}}(X_{*},Y)\,dY\approx\sigma(Q({\bf S}_{*}))\int_{\partial\Omega_{*}}(u(y)-u(X_{*}))^{2}\,d\omega_{*}(y),

where

$\omega_{*}$ is the elliptic measure for $L$ in $\Omega_{*}$ with pole at $X_{*}$ .

Thus, Lemma 4.9 follows immediately from the following estimate. Recall that $\lambda$ is the parameter we used when we defined the Types 1–4 for the stopping time regimes.

Lemma 4.12.

There exists $\lambda_{0}>0$ depending on structural constants, ellipticity, and the $\omega_{L,\Omega}\in A_{\infty}(\sigma)$ constants, such that for any $\lambda\in(0,\lambda_{0})$ , there exists a constant $c_{4}>0$ depending only on structural constants, ellipticity, and the $\omega_{L,\Omega}\in A_{\infty}(\sigma)$ constants (and independent of $\varepsilon$ , $\lambda$ , $N$ , $j$ , and $\mathcal{S}$ ), so that the following estimate holds:

\int_{\partial\Omega_{*}}(u(y)-u(X_{*}))^{2}\,d\omega_{*}(y)\geq c_{4}\varepsilon^{2}.

For the the proof of Lemma 4.12, we need some auxiliary constructions and estimates. Recall that $X_{*}=X_{Q({\bf S}_{*})}$ is a corkscrew point relative to $Q({\bf S}_{*})$ at scale $r_{*}\coloneqq 10^{-5}a_{0}\ell(Q({\bf S}_{*}))$ (see Section 2.4). Let $\hat{x}_{*}\in\partial\Omega$ be a touching point for $X_{*}$ on $\partial\Omega$ , that is, $|\hat{x}_{*}-X_{*}|=\operatorname{dist}(X_{*},\partial\Omega)$ . For $\xi\in[0,1]$ , consider the points $X(\xi)\coloneqq\hat{x}_{*}+\xi(X_{*}-\hat{x}_{*})$ which lie on the line segment from $\hat{x}_{*}$ to $X_{*}$ . Since we know that $\hat{x}_{*}\not\in\overline{\Omega}_{*}$ and $X_{*}\in\Omega_{*}$ , there exists $\xi_{0}\in(0,1)$ such that $X(\xi_{0})\in\partial\Omega_{*}$ . Now we set $X_{**}=X(\xi_{0})$ , and $\Delta_{*}\coloneqq B(X_{**},r_{*})\cap\partial\Omega_{*}$ . Since we are working with the truncated collection of stopping cubes $\mathcal{F}=\mathcal{F}_{N,j}$ , we know that $\partial\Omega\cap\partial\Omega_{*}=\mbox{{\O}}$ . In particular, $\Delta_{*}\subset\Omega$ . By Lemma 2.26, we know that

(4.13)

\displaystyle\text{if }y\in\Delta_{*},\ \ \text{ then }\hat{y}\in\Delta_{Q({\bf S}_{*})}\subseteq Q({\bf S}_{*}),

where $\hat{y}$ is a touching point for $y$ in $\Omega$ , that is, $|y-\hat{y}|=\operatorname{dist}(y,\partial\Omega)$ .

By Lemma 2.18, we know that $\Omega_{*}$ is also a uniform domain with Ahlfors regular boundary.⁵⁵5By Lemma 2.18, the structural constants of $\Omega_{*}=\Omega_{\mathcal{F}_{N,j},Q({\bf S}_{j})}$ do not depend on $j$ or the truncation parameter $N=N(j)$ . Thus, by Lemma 2.36, we have

(4.14)

\displaystyle\omega^{\widetilde{X}}_{L,\Omega_{*}}(\Delta_{*})\gtrsim 1,

where $\widetilde{X}$ is a corkscrew point relative to $X_{**}$ at scale $r_{*}$ in the domain $\Omega_{*}$ . Recall that by the construction in Section 2.3, we know that $B(X_{*},\delta(X_{*})/2)\subset\Omega_{*}$ and $\operatorname{diam}(\Omega_{*})\approx\ell(Q({\bf S}_{*}))$ , and we have $\delta(X_{*})\approx\ell(Q({\bf S}_{*}))\approx r_{*}\approx\delta(\widetilde{X})$ by the definition of $X_{*}$ and $\widetilde{X}$ . Thus, by the Harnack chain property of $\Omega_{*}$ , there exists a Harnack chain of uniformly bounded length from $\widetilde{X}$ to $X_{*}$ inside $\Omega_{*}$ . Thus, by (4.14), formula (2.30) and Lemma 2.29 (that is, Harnack inequality), there exists a constant $c_{*}>0$ that depends only on structural constants such that

(4.15)

\displaystyle\omega_{*}(\Delta_{*})=\omega^{X_{*}}_{L,\Omega_{*}}(\Delta_{*})>c_{*}.

Next, we will construct a cover of $\Delta_{*}$ that consists of dilated surface balls on $\partial\Omega_{*}$ associated to the cubes $Q\in\mathcal{F}$ . We will construct the cover in such a way that there is oscillation of $u$ on the balls associated to cubes in $\mathcal{F}^{\operatorname{SB}}$ and the balls associated to cubes in $\mathcal{F}^{\operatorname{O}}$ do not have much $\omega_{*}$ -mass, provided $\lambda$ is sufficiently small. Given $Q\in\mathbb{D}$ , denote by $\hat{x}_{Q}$ a touching point for the corkscrew point $X_{Q}$ . For $\theta\in[0,1]$ , recall that we denote $P_{Q}(\theta)=\hat{x}_{Q}+\theta(X_{Q}-\hat{x}_{Q})$ , and by Lemma 2.21 we showed that there exists $\theta_{0}\in(0,1)$ such that if for some $Q^{\prime}\in\mathbb{D}$ we have $B(P_{Q}(\theta_{0}),\tfrac{\gamma\theta_{0}}{10}r_{Q})\cap U_{Q^{\prime}}\neq\mbox{{\O}}$ , then $Q^{\prime}\subset Q$ and $\ell(Q^{\prime})<\ell(Q)$ , where $\gamma$ is the corkscrew constant in Definition 2.1.

Fix $Q\in\mathcal{F}$ . Then its parent $\widetilde{Q}$ satisfies $\widetilde{Q}\in\mathbb{D}_{\mathcal{F},Q({\bf S}_{*})}$ and hence, by the construction of the Whitney regions in Section 2.3, we have $P_{Q}(1)=X_{Q}\in U_{\widetilde{Q}}\subset\Omega_{*}$ . By Lemma 2.21, we also know that $P_{Q}(\theta_{0})\not\in\Omega_{*}$ . Indeed, otherwise $P_{Q}(\theta_{0})\in U_{Q^{\prime}}$ for a cube $Q^{\prime}\in\mathbb{D}_{\mathcal{F},Q({\bf S}_{*})}$ , but by Lemma 2.21 we would then have $Q^{\prime}\subset Q$ with $\ell(Q^{\prime})<\ell(Q)$ . This is impossible since $Q\in\mathcal{F}$ . Thus, there exists $\theta^{\prime}\in(\theta_{0},1)$ such that

\displaystyle X^{*}_{Q}\coloneqq P_{Q}(\theta^{\prime})\in\partial\Omega_{*}.

We set

\displaystyle\Delta^{*}_{Q}\coloneqq B(X^{*}_{Q},\tfrac{\lambda\theta_{0}}{2}r_{Q})\cap\partial\Omega_{*}\quad\text{ and }\quad M\Delta^{*}_{Q}\coloneqq B(X^{*}_{Q},M\tfrac{\lambda\theta_{0}}{2}r_{Q})\cap\partial\Omega_{*}

for a constant $M\geq 1$ to be chosen momentarily.

Let us describe some of the properties of $\Delta^{*}_{Q}$ . First, by Lemma 2.25 and definition of $\Delta_{Q}^{*}$ , it holds that

(4.16)

\Delta^{*}_{Q}\subseteq\Xi_{Q}\subseteq U_{Q^{\prime}}^{*}

whenever $Q^{\prime}$ is a sibling of $Q$ , where

\displaystyle\Xi_{Q^{\prime}}=\bigcup_{\theta\in[\theta_{0},1]}B\big{(}P_{Q^{\prime}}(\theta),\tfrac{\gamma\theta_{0}}{10}r_{Q^{\prime}}\big{)}.

Next, let us observe that if $Q\in\mathcal{F}^{\operatorname{SB}}$ , then there exists a sibling $Q^{\prime}$ of $Q$ such that $Q^{\prime}$ is blue and $|u(X_{Q^{\prime}})-u(X_{Q({\bf S}_{*})})|>\varepsilon/100$ . Thus, for every $Q\in\mathcal{F}^{\operatorname{SB}}$ , it holds that

(4.17)

|u(X)-u(X_{*})|\geq|u(X_{Q^{\prime}})-u(X_{Q({\bf S}_{*})})|-|u(X_{Q^{\prime}})-u(X)|\\ \geq\varepsilon/100-\varepsilon/1000\geq\varepsilon/200,

for every $X\in\Delta^{*}_{Q}$ since $Q^{\prime}$ is blue and $\Delta^{*}_{Q}\subset U_{Q^{\prime}}^{*}$ by (4.16).

Lemma 4.18.

There exists $M\geq 1$ , depending only on structural constants, such that

\displaystyle\Delta_{*}=B(X_{**},r_{*})\cap\partial\Omega_{*}\subset\bigcup_{Q\in\mathcal{F}}M\Delta^{*}_{Q}.

Proof.

Fix $y\in\Delta_{*}\subset\partial\Omega_{*}$ and a touching point $\hat{y}\in\partial\Omega$ for $y$ . Then by (4.13) we have $\hat{y}\in\Delta_{Q({\bf S}_{*})}$ . Since $X_{**}$ lies on the line segment from the corkscrew point $X_{*}$ relative to $x_{Q({\bf S}_{*})}$ at scale $r_{*}$ to its touching point $\hat{x}_{*}$ , we have $|\hat{x}_{*}-X_{**}|=\operatorname{dist}(X_{**},\partial\Omega)\leq\operatorname{dist}(X_{*},\partial\Omega)\leq r_{*}=10^{-5}a_{0}\ell(Q({\bf S}_{*}))$ . Thus, $y\in B(\hat{x}_{*},2r_{*})=B(\hat{x}_{Q({\bf S}_{*})},2(10)^{-5}a_{0}\ell(Q({\bf S}_{*})))$ . By (4.13), we have

(4.19)

|y-\hat{y}|=\operatorname{dist}(y,\partial\Omega)=\operatorname{dist}(y,Q({\bf S}))\leq\ell(Q({\bf S}_{*})),

and by the definition of $\mathcal{W}_{Q^{\prime}}(K_{0})$ (which we used in the construction of the Whitney regions), for any cube $Q^{\prime}$ it holds that

(4.20)

\displaystyle\text{if }\hat{y}\in Q^{\prime}\text{ and }C_{\tau}^{-1}K_{0}^{-1}|y-\hat{y}|\leq\ell(Q^{\prime})\leq C_{\tau}K_{0}|y-\hat{y}|,\quad\text{ then }y\in\operatorname{int}(U_{Q^{\prime}}),

where $\tau$ is the dilation parameter in the definition of $U_{Q}$ . Since ${\bf S}_{*}$ and $\mathcal{F}$ are the truncated collections, we have $Q({\bf S}_{*})=\cup_{Q\in\mathcal{F}}Q$ . In particular, by Lemma 2.26, we have $\hat{y}\in Q({\bf S}_{*})=\cup_{Q\in\mathcal{F}}Q$ . Let $Q_{\hat{y}}\in\mathcal{F}$ be the cube such that $\hat{y}\in Q_{\hat{y}}$ . We now have $\ell(Q_{\hat{y}})\geq C_{\tau}^{-1}K_{0}^{-1}|y-\hat{y}|$ for the same constant $C_{\tau}$ as in (4.20), since otherwise there exists a cube $Q^{\prime}$ such that $Q_{\hat{y}}\subset Q^{\prime}\subseteq Q({\bf S}_{*})$ with $C^{-1}_{\tau}K_{0}^{-1}|y-\hat{y}|\leq\ell(Q^{\prime})\leq C_{\tau}K_{0}|y-\hat{y}|$ . This is not possible since (4.20) and the fact that $Q_{\hat{y}}\in\mathcal{F}$ would then imply that $y\in\text{int}\,(U_{Q^{\prime}})\subset\Omega_{*}$ , but we know that $y\in\partial\Omega_{*}$ . Thus, it holds that

(4.21)

\displaystyle\operatorname{dist}(Q_{\hat{y}},y)=|\hat{y}-y|\leq C_{\tau}K_{0}\ell(Q_{\hat{y}}).

Recall that $x_{Q_{\hat{y}}}$ is the center of $Q_{\hat{y}}$ and $X_{Q_{\hat{y}}}^{*}$ lies on a line segment from a corkscrew point $X_{Q_{\hat{y}}}$ to its touching point $\hat{x}_{Q_{\hat{y}}}$ . By Lemma 2.26, we know that $\hat{x}_{Q_{\hat{y}}}\in Q_{\hat{y}}$ . Thus, (4.21), the definitions of the points and the fact that $\hat{y},\hat{x}_{Q_{\hat{y}}}\in Q_{\hat{y}}$ give us

\displaystyle|y-X^{*}_{Q_{\hat{y}}}|\leq|y-\hat{y}|+|\hat{y}-\hat{x}_{Q_{\hat{y}}}|+|\hat{x}_{Q_{\hat{y}}}-X^{*}_{Q_{\hat{y}}}|\lesssim K_{0}\ell(Q_{\hat{y}})+\operatorname{diam}(Q_{\hat{y}})+r_{Q_{\hat{y}}}\approx r_{Q_{\hat{y}}}.

In particular, there exists $M\geq 1$ such that

\displaystyle y\in M\Delta^{*}_{Q_{\hat{y}}}\subseteq\bigcup_{Q\in\mathcal{F}}M\Delta^{*}_{Q},

which is what we wanted. ∎

Let us then fix $M\geq 1$ as in Lemma 4.18. By (4.15), it holds that

(4.22)

\displaystyle\omega_{*}\Big{(}\bigcup_{Q\in\mathcal{F}}M\Delta^{*}_{Q}\Big{)}\geq c_{*}.

Our next goal is to analyze how much the cubes $Q\in\mathcal{F}^{\operatorname{O}}$ contribute to (4.22) and then limit this contribution by choosing $\lambda$ in a suitable way. For this, we prove the following bound:

Lemma 4.23.

For any $Q\in\mathcal{F}$ , we have

(4.24)

\displaystyle\omega_{*}(\Delta^{*}_{Q})\lesssim\frac{\sigma(Q)}{\sigma(Q({\bf S}_{*}))}

for an implicit constant depending only on structural constants, ellipticity, and the $\omega_{L,\Omega}\in A_{\infty}(\sigma)$ constants (and independent of $\lambda$ , $\varepsilon$ , $N$ , $j$ , and $\mathcal{S}$ ).

Proof.

Let $Q\in\mathcal{F}$ . Recall that $\Delta^{*}_{Q}=B(X^{*}_{Q},\tfrac{\lambda\theta_{0}}{2}r_{Q})\cap\partial\Omega_{*}$ is a surface ball on $\partial\Omega_{*}$ with $X^{*}_{Q}\in\partial\Omega_{*}$ . Since $\Omega_{*}$ is a uniform domain, it satisfies the corkscrew condition. Let $\widetilde{X}_{Q}$ be a corkscrew point in $\Omega_{*}$ relative to $X^{*}_{Q}$ at scale approximately $r_{Q}^{\prime}=\frac{\lambda\theta_{0}}{2}r_{Q}\approx r_{Q}$ . By perhaps insisting that $\theta_{0}$ is smaller we have that this corkscrew point $\widetilde{X}_{Q}$ is far from $X_{*}$ . Then by connecting $\widetilde{X}_{Q}$ to $X_{Q}$ with a Harnack chain⁶⁶6Here we use that $\operatorname{dist}(X^{*}_{Q},\partial\Omega)\approx\operatorname{dist}(X^{*}_{Q},Q)\approx\ell(Q)\approx r_{Q}$ ; with very crude bounds this can be seen from (4.16). Then we also use that $\Omega_{*}\subset\Omega$ , so that $\ell(Q)\approx r_{Q}^{\prime}\approx\operatorname{dist}(\widetilde{X}_{Q},\partial\Omega_{*})\leq\operatorname{dist}(\widetilde{X}_{Q},\partial\Omega)$ . (of uniformly bounded length) in $\Omega$ and using Lemma 2.38,

G_{L,\Omega}(X_{*},\widetilde{X}_{Q})\approx G_{L,\Omega}(X_{*},X_{Q})\lesssim\ell(Q)\frac{\omega^{X_{*}}_{L,\Omega}(Q)}{\ell(Q)^{n}}.

Now, by 3.16 we have that $G_{L,\Omega_{*}}(X_{*},\widetilde{X}_{Q})\leq G_{L,\Omega}(X_{*},\widetilde{X}_{Q})$ , and then by Lemma 2.38 in $\Omega_{*}$ we conclude that

\frac{\omega_{*}(\Delta^{*}_{Q})}{(r^{\prime}_{Q})^{n}}r^{\prime}_{Q}=\frac{\omega^{X_{*}}_{L,\Omega_{*}}(\Delta^{*}_{Q})}{(r^{\prime}_{Q})^{n}}r^{\prime}_{Q}\lesssim G_{L,\Omega_{*}}(X_{*},\widetilde{X}_{Q})\leq G_{L,\Omega}(X_{*},\widetilde{X}_{Q}).

Combining the two previously displayed inequalities and using that $\ell(Q)\approx r^{\prime}_{Q}$ we have $\omega_{*}(\Delta^{*}_{Q})\lesssim\omega^{X_{*}}_{L,\Omega}(Q)$ . By Lemma 2.40 and Lemma 3.18 applied twice for both $Q,Q({\bf S}_{*})\in\mathcal{S}$ , it holds that

\omega^{X_{*}}_{L,\Omega}(Q)=\omega^{X_{Q({\bf S}_{*})}}_{L,\Omega}(Q)\approx\frac{\omega_{L,\Omega}^{X_{Q(\mathcal{S})}}(Q)}{\omega_{L,\Omega}^{X_{Q(\mathcal{S})}}(Q({\bf S}_{*}))}\approx\frac{\sigma(Q)}{\sigma(Q(\mathcal{S}))}\frac{\sigma(Q(\mathcal{S}))}{\sigma(Q({\bf S}_{*}))}=\frac{\sigma(Q)}{\sigma(Q({\bf S}_{*}))},

which ends the proof of (4.24). ∎

Now we are ready to conclude the proof of Lemma 4.12. Using (4.24) and the doubling property of $\omega_{*}$ we find that

(4.25)

\sum_{Q\in\mathcal{F}^{\operatorname{O}}}\omega_{*}\left(M\Delta^{*}_{Q}\right)\leq C\sum_{Q\in\mathcal{F}^{\operatorname{O}}}\omega_{*}\left(\Delta^{*}_{Q}\right)\leq C\sum_{Q\in\mathcal{F}^{\operatorname{O}}}\frac{\sigma(Q)}{\sigma(Q({\bf S}_{*}))}\leq\hat{C}\lambda,

where we used (4.8) in the last inequality. Now we choose $\lambda>0$ so that $\hat{C}\lambda<c_{*}/2$ and use (4.25) and (4.22) to deduce

(4.26)

\omega_{*}\Big{(}\bigcup_{Q\in\mathcal{F}^{\operatorname{SB}}}M\Delta^{*}_{Q}\Big{)}\geq c_{*}/2,

where we used $\mathcal{F}=\mathcal{F}^{\operatorname{SB}}\cup\mathcal{F}^{O}$ . Now we use the $5R$ -covering lemma [Mat95] to produce a countable collection of disjoint surface balls $\{M\Delta^{*}_{k}\}\coloneqq\{M\Delta^{*}_{Q_{k}}\}$ where each $Q_{k}$ is in $\mathcal{F}^{\operatorname{SB}}$ and such that

\bigcup_{Q\in\mathcal{F}^{\operatorname{SB}}}M\Delta^{*}_{Q}\subset\cup_{k}5M\Delta^{*}_{k}.

Then using (4.26) and the doubling property of $\omega_{*}$ it holds

\omega_{*}(\cup_{k}\Delta^{*}_{k})=\sum_{k}\omega_{*}(\Delta^{*}_{k})\gtrsim\sum_{k}\omega_{*}(5M\Delta^{*}_{k})\gtrsim\omega_{*}\Big{(}\bigcup_{Q\in\mathcal{F}^{\operatorname{SB}}}M\Delta^{*}_{Q}\Big{)}\geq c_{*}/2,

where we used that $\Delta^{*}_{k}$ are disjoint. To summarize we have produced a sequence of surface balls $\Delta^{*}_{k}=\Delta^{*}_{Q_{k}}$ with $Q_{k}\in\mathcal{F}^{\operatorname{SB}}$ such that

(4.27)

\omega_{*}(\cup_{k}\Delta^{*}_{k})\geq c_{**},

where $c_{**}$ depends on dimension, ellipticity, the Ahlfors regularity constant for $\partial\Omega$ , the corkscrew and Harnack Chain constants for $\Omega$ , and the $\omega_{L}\in A_{\infty}(\sigma)$ constants. Thus, using (4.17) we have

\int_{\partial\Omega_{*}}(u(y)-u(X_{*}))^{2}\,d\omega_{*}(y)\geq\omega_{*}(\cup_{k}\Delta^{*}_{k})\inf\limits_{y\in\cup_{k}\Delta^{*}_{k}}(u(y)-u(X_{*}))^{2}\\ \geq c_{**}(\varepsilon/200)^{2}\eqqcolon c_{4}\varepsilon^{2},

which proves Lemma 4.12.

As we had reduced the proof of the packing of the Type 4 maximal cubes to Lemma 4.12, this completes the proof of Lemma 4.6.

4.4. Construction of $\varepsilon$ -approximators

With the help of the previous constructions and estimates, we can prove the existence of $\operatorname{BV}_{\operatorname{loc}}$ $\varepsilon$ -approximators in a similar way as in [HMM16]. For the convenience of the reader, we recall the key steps of the construction below. For some of the details, we follow the construction of $L^{p}$ -type approximators in [HT20] which are an adaptation of the arguments in [HMM16]. Recall that we denote the collection of blue cubes in the disjoint stopping time regimes $\mathcal{S}$ in Lemma 3.18 by $\mathcal{L}=\mathcal{L}(\mathcal{S})$ , and each of these collections has a decomposition $\mathcal{L}=\cup_{j}{\bf S}_{j}$ .

Proof of Theorem 4.1.

Let us fix a dyadic cube $Q_{0}\in\mathbb{D}(\partial\Omega)$ and construct an $\varepsilon$ -approximator $\Phi_{Q_{0}}=\Phi_{Q_{0}}^{\varepsilon}$ first in the Carleson box $T_{Q_{0}}$ . We start by dividing the Carleson box $T_{Q_{0}}$ into a few types of different regions where we define the approximator differently. Let us choose the largest good (in the sense of the corona decomposition from Lemma 3.18), blue subcube $Q_{1}\subset Q_{0}$ which may be $Q_{0}$ itself; if there are several such cubes with the largest side length, we choose just one of them. Since $Q_{1}$ is a good blue cube, there exists a stopping time regime $\mathcal{S}_{Q_{1}}$ in Lemma 3.18 and subregime ${\bf S}_{Q_{1}}\subset\mathcal{L}(\mathcal{S}_{Q_{1}})$ such that $Q_{1}$ is the maximal element of the regime ${\bf S}^{1}\coloneqq{\bf S}_{Q_{1}}\cap\mathbb{D}_{Q_{0}}$ . We then choose the largest good blue cube $Q_{2}$ from the collection $\mathbb{D}_{Q_{0}}\setminus{\bf S}^{1}$ . Similarly, $Q_{2}$ is the maximal element of the regime ${\bf S}^{2}\coloneqq{\bf S}_{Q_{2}}\cap\mathbb{D}_{Q_{0}}$ . We then choose the largest good blue cube $Q_{3}\in\mathbb{D}_{Q_{0}}\setminus({\bf S}^{1}\cup{\bf S}^{2})$ , and continue like this. This gives us a sequence of good blue cubes $Q_{1},Q_{2},\ldots$ such that $\ell(Q_{1})\geq\ell(Q_{2})\geq\ldots$ , each cube $Q_{i}$ is a maximal element of a regime ${\bf S}^{i}$ and the collection $\cup_{i}{\bf S}^{i}$ contains all the good blue cubes in $\mathbb{D}_{Q_{0}}$ . The cubes $Q_{i}$ are “mostly” of the form $Q({\bf S})$ as in the decomposition of the collections $\mathcal{L}(\mathcal{S})$ earlier in the sense that there exists a collection of pairwise disjoint cubes $\{P_{k}\}_{k}\subset\mathbb{D}_{Q_{0}}$ (that may be empty) such that every cube in the collection $\{Q_{i}\}_{i}\setminus\{P_{k}\}_{k}$ is of the form of $Q({\bf S})$ for some ${\bf S}$ . This is because the cube $Q_{0}$ is arbitrary and hence, it may be a bad cube or a red cube. For each $i$ , we define the “bottom” cubes of ${\bf S}^{i}$ in the obvious way: we set $\mathcal{F}_{{\bf S}^{i}}\coloneqq\mathcal{F}_{{\bf S}_{Q_{i}}}$ , that is, $\mathcal{F}_{{\bf S}^{i}}$ is the collection of the stopping cubes associated to the unique regime ${\bf S}_{Q_{i}}$ that contains $Q_{i}$ .

For each $i$ , we define the regions $A_{i}$ recursively the following way:

\displaystyle A_{1}\coloneqq\Omega_{\mathcal{F}_{{\bf S}^{1}},Q_{1}},\qquad A_{i}\coloneqq\Omega_{\mathcal{F}_{{\bf S}^{i},Q_{i}}}\setminus\bigcup_{k=1}^{i-1}A_{k}\quad\text{ for }i\geq 2.

By construction, the regions $A_{i}$ are pairwise disjoint. We also set $\Omega_{0}\coloneqq\bigcup_{i}A_{i}$ , and we define the function $\Phi_{0}$ on $\Omega_{0}$ as

\displaystyle\Phi_{0}\coloneqq\sum_{i}u(X_{Q({\bf S}_{Q_{i}})})\mathbbm{1}_{A_{i}},

where $X_{Q({\bf S}_{Q_{i}})}$ is the corkscrew point we used in the stopping conditions in the definition of ${\bf S}_{Q_{i}}$ . In particular, for any $X\in A_{i}$ we have $|u(X)-u(X_{Q({\bf S}_{Q_{i}})})|\leq\varepsilon/100<\varepsilon$ . Furthermore, by the disjointness of the regions $A_{i}$ , we have $\|u-\Phi_{0}\|_{L^{\infty}(\Omega_{0})}<\varepsilon$ .

Let us then consider the cubes in $\mathbb{D}_{Q_{0}}\setminus\cup_{i}{\bf S}^{i}$ . Let us fix some enumeration $\{R_{j}\}_{j}$ for the cubes $\mathbb{D}_{Q_{0}}\setminus\cup_{i}{\bf S}^{i}$ . The cubes $R_{j}$ are red cubes or bad blue cubes. For each $j$ , we define the regions $V_{j}$ recursively the following way:

\displaystyle V_{1}\coloneqq U_{R_{1}},\qquad V_{j}\coloneqq U_{R_{j}}\setminus\bigcup_{k=1}^{j-1}V_{k}\quad\text{ for }j\geq 2.

By construction, the regions $V_{j}$ are pairwise disjoint. We also set $\Omega_{1}\coloneqq\bigcup_{j}V_{j}$ , and we define the function $\Phi_{1}$ on $\Omega_{1}$ as

\displaystyle\Phi_{1}(X)\coloneqq\left\{\begin{array}[]{cl}u(X),&\text{ if }X\in V_{k}\text{ for a red cube }R_{k},\\ u(X_{k}),&\text{ if }X\in V_{k}\text{ for a blue cube }R_{k},\end{array}\right.

where $X_{k}$ is any fixed point on $U_{R_{k}}$ . By the definitions, we have $\|u-\Phi_{1}\|_{L^{\infty}(\Omega_{1})}<\varepsilon/1000<\varepsilon$ .

We define the $\varepsilon$ -approximator $\Phi_{Q_{0}}$ of $u$ in the Carleson box $T_{Q_{0}}$ as

\displaystyle\Phi_{Q_{0}}(X)\coloneqq\left\{\begin{array}[]{cl}\Phi_{0}(X),&\text{ if }X\in\Omega_{0},\\ \Phi_{1}(X),&\text{ if }X\in T_{Q_{0}}\setminus\Omega_{0}.\end{array}\right.

By the construction, we have $\|u-\Phi_{Q_{0}}\|_{L^{\infty}(T_{Q_{0}})}<\varepsilon$ . The $L^{1}$ -type Carleson measure estimate for $\Phi_{Q_{0}}$ in $T_{Q_{0}}$ can be proven as in [HMM16] with small but quite obvious changes. Using a covering argument, the claim can be reduced to proving the estimate on Carleson boxes $T_{Q^{\prime}}$ , and since $u\in L^{\infty}(\Omega)$ , the core challenge is to handle the jumps across the boundaries of the sets $A_{i}$ and $V_{j}$ that contribute to the total variation of $\Phi_{Q_{0}}$ inside $T_{Q_{0}}$ . Since the boundaries of the sawtooth regions, Whitney regions and Carleson boxes are Ahlfors regular by Lemma 2.18, the estimates reduce to using the Carleson packing conditions in Lemma 3.18, Lemma 4.2, Lemma 4.5 and Lemma 4.6. The Carleson norm of the measure $\mu_{\Phi_{Q_{0}}}$ such that $d\mu_{\Phi_{Q_{0}}}(Y)=|\nabla\Phi_{Q_{0}}(Y)|\,dY$ is given (up to a structural constant) by the sizes of the Carleson packing norms in these results. We omit the details.

Using these kinds of local approximators, we build the global approximator of $u$ . If $\operatorname{diam}(\Omega)<\infty$ , it is enough to build a local approximator for a Carleson box that covers the whole space $\Omega$ . Thus, we may assume that $\operatorname{diam}(\Omega)=\infty$ . Suppose first that $\operatorname{diam}(\partial\Omega)<\infty$ . Then there exists a dyadic cube $Q_{0}$ that covers the whole boundary $\partial\Omega$ . We build the local approximator $\Phi_{Q_{0}}$ on $T_{Q_{0}}$ , extend it to whole $\Omega$ by setting it to be $0$ outside $T_{Q_{0}}$ and define the global approximator as $\Phi=\mathbbm{1}_{T_{Q_{0}}}\Phi_{Q_{0}}+\mathbbm{1}_{\Omega\setminus T_{Q_{0}}}u$ . The $L^{1}$ -type Carleson measure estimate follows from the same arguments as with the local approximators.

Finally, suppose that $\operatorname{diam}(\partial\Omega)=\infty$ . Fix a sequence of dyadic cubes $P_{k}$ such that $P_{1}\subset P_{2}\subset\cdots$ , $\ell(P_{1})<\ell(P_{2})<\cdots$ and $\partial\Omega=\cup_{k}P_{k}$ . This type of sequence of cubes does not exist in every dyadic system, but we can always construct a system where it exists (see, for example, [HT14]). We build a local approximator $\Phi_{P_{k}}$ in $T_{P_{k}}$ for every $k$ and extend the approximators to whole $\Omega$ by setting each of them to be $0$ outside $T_{P_{k}}$ . We then define the global approximator as $\Phi=\mathbbm{1}_{T_{P_{1}}}\Phi_{P_{1}}+\sum_{k=2}^{\infty}\mathbbm{1}_{T_{P_{k}}\setminus T_{P_{k-1}}}\Phi_{P_{k}}$ . The $L^{1}$ -type Carleson measure estimate follows from the Carleson measure estimates of the local approximators and the fact that the collection $\{P_{k}\}_{k}$ satisfies a Carleson packing condition with a uniformly bounded Carleson packing norm depending only on structural constants. Again, we omit the details. ∎

To finish the proof of Theorem 1.1, we regularize the approximators in Theorem 4.1. This regularization makes the constant $C_{\varepsilon}$ significantly larger but since the size of this constant is not important for our results, we do not track its size.

Lemma 4.28.

Let $\varepsilon\in(0,1)$ . There exists a unifomly bounded constant $\widetilde{C}_{\varepsilon}\geq 1$ such that we can choose the $\varepsilon$ -approximator $\Phi=\Phi^{\varepsilon}$ for the solution $u\in W^{1,2}(\Omega)\cap L^{\infty}(\Omega)$ to $Lu=0$ in Theorem 4.1 so that

i)

$\|u-\Phi\|_{L^{\infty}(\Omega)}\leq 2\varepsilon\|u\|_{L^{\infty}(\Omega)}$ ,
ii)

$\sup_{x\in\partial\Omega,r>0}\frac{1}{r^{n}}\iint_{B(x,r)\cap\Omega}|\nabla\Phi(Y)|\,dY\leq\widetilde{C}_{\varepsilon}\|u\|_{L^{\infty}(\Omega)},$
iii)

$\Phi\in C^{\infty}(\Omega)$ ,
iv)

$|\nabla\Phi(Y)|\leq\tfrac{\widetilde{C}_{\varepsilon}}{\delta(Y)}$ for every $Y\in\Omega$ ,
v)

if $|X-Y|\ll\delta(X)$ , then $|\Phi(X)-\Phi(Y)|\leq\tfrac{\widetilde{C}_{\varepsilon}|X-Y|}{\delta(X)}$ ,
vi)

there exists a function $\varphi\in L^{\infty}(\partial\Omega)$ such that

$\displaystyle\lim_{Y\to x,\,\operatorname{n.t.}}\Phi(Y)=\varphi(x)\text{ for }\sigma\text{-a.e. }x\in\partial\Omega.$

The constant $\widetilde{C}_{\varepsilon}$ depends on $\varepsilon$ , the structural constants of $\Omega$ , the constant $C_{\varepsilon}$ in Theorem 4.1 and the Hölder continuity constants $C$ and $\alpha$ in Lemma 2.28.

Proof.

The proof uses tweaked mollifier techniques combined with a regularized distance function. The properties follow mostly from [HT21, Section 3] but for the convenience of the reader, we define the core objects and give some explicit details below.

Let $\beta$ be a regularized version of the distance function $\delta=\operatorname{dist}(\cdot,\partial\Omega)$ , that is, a smooth function in $\Omega$ such that $\beta\approx\delta$ (see [Ste70, Theorem 2, p. 171]). Let $\zeta\geq 0$ be a smooth non-negative function supported on $B(0,\tfrac{1}{m})$ for a suitable constant $m>0$ (depending on the implicit constants in $\delta\approx\beta$ ), satisfying $\zeta\leq 1$ and $\int\zeta=1$ . For a constant $\xi_{\varepsilon}>0$ to be chosen momentarily, we set

\displaystyle\Lambda_{\xi_{\varepsilon}}(X,Y)\coloneqq\zeta_{\xi_{\varepsilon}\beta(X)}(X-Y)=\frac{1}{(\xi_{\varepsilon}\beta(X))^{n+1}}\zeta\Big{(}\frac{X-Y}{\xi_{\varepsilon}\beta(X)}\Big{)}.

For a suitable choice of $m$ , we have $\text{supp}\,\Lambda(X,\cdot)\subset B(X,\xi_{\varepsilon}\delta(X)/2)$ . Given the non-smooth $\varepsilon$ -approximator $\Phi_{0}$ of the solution $u\in W^{1,2}_{\operatorname{loc}}(\Omega)\cap L^{\infty}(\Omega)$ to $Lu=0$ in Theorem 4.1, we set

\displaystyle\Phi(X)\coloneqq\iint\Lambda_{\varepsilon}(X,Y)\Phi_{0}(Y)\,dY.

The property iii) follows from a standard modification of the case $\Omega=\mathbb{R}^{n+1}$ (for example, see [EG15, Theorem 1, p. 123]) and properties ii), iv) and v) are formulated explicitly in [HT21, Section 3]. Property i) follows from the local Hölder continuity of $u$ (that is, Lemma 2.28) and the fact that $\Phi_{0}$ is an $\varepsilon$ -approximator of $u$ : for almost every $X\in\Omega$ , we get

	$\displaystyle\|\Phi(X)-u(X)\|$	$\displaystyle=\Big{\|}\iint\Lambda_{\varepsilon}(X,Y)(\Phi_{0}(Y)-u(X))\,dY\Big{\|}$
		$\displaystyle<\varepsilon\\|u\\|_{L^{\infty}(\Omega)}+\iint\Lambda_{\varepsilon}(X,Y)\left\|u(Y)-u(X)\right\|\,dY$
		$\displaystyle\leq\varepsilon\\|u\\|_{L^{\infty}(\Omega)}+\iint\Lambda_{\varepsilon}(X,Y)\,C\Big{(}\frac{\|X-Y\|}{\tfrac{1}{4}\delta(X)}\Big{)}^{\alpha}\Big{(}\fint_{B(X,\tfrac{1}{2}\delta(X))}\|u(Z)\|^{2}\,dZ\Big{)}^{1/2}\,dY$
		$\displaystyle\leq\varepsilon\\|u\\|_{L^{\infty}(\Omega)}+C\left(2\xi_{\varepsilon}\right)^{\alpha}\\|u\\|_{L^{\infty}(\Omega)}\iint\Lambda_{\varepsilon}(X,Y)\,dY\leq 2\varepsilon\\|u\\|_{L^{\infty}(\Omega)}$

as long as we choose $\xi_{\varepsilon}\leq\tfrac{1}{2}\big{(}\tfrac{\varepsilon}{C}\big{)}^{\frac{1}{\alpha}}$ , where $C$ and $\alpha$ are the Hölder continuity constants in Lemma 2.28.

Property vi) follows from the same argument that is used in the proof of [HT21, Lemma 4.14] after some small additional considerations. The proof of [HT21, Lemma 4.14] is based on showing that almost every cone on a codimension $1$ uniformly rectifiable set has locally exactly two components, these local components satisfy the Harnack chain condition and the Harnack chain condition combined with the $L^{1}$ -type Carleson measure estimate ii) implies the existence of the a.e. non-tangential trace $\varphi$ . We do not assume that $\partial\Omega$ is uniformly rectifiable but by the definition of dyadic cones (2.17) and Lemma 2.18 we know that any truncated cone on $\partial\Omega$ has exactly one component inside $\Omega$ and this component satisfies the Harnack chain condition. Thus, the argument in the proof of [HT21, Lemma 4.14] works also for us. In particular, we can choose $\Phi$ in such a way that all the properties i) – vi) hold. This completes the proofs of Lemma 4.28 and Theorem 1.1. ∎

5. Proof of Theorem 1.2

In this section, we prove Theorem 1.2, that is, we prove that $\varepsilon$ -approximability of solutions $u$ to $Lu=0$ implies that $\omega_{L}\in A_{\infty}(\sigma)$ . To be more precise, we prove the following seemingly stronger result:

Theorem 5.1.

Let $\Omega\subset\mathbb{R}^{n+1}$ , $n\geq 1$ , be a uniform domain with Ahflors regular boundary, and let $L$ be a divergence form elliptic operator $L=-\mathop{\operatorname{div}}\nolimits A\nabla$ in $\Omega$ . Suppose also that for every bounded Borel set $S\subset\partial\Omega$ the solution $u=u_{S}$ to $Lu=0$ such that $u(X)=\omega_{L}^{X}(S)$ is $\varepsilon$ -approximable for every $\varepsilon\in(0,1)$ in the sense of Theorem 1.1 with the $\varepsilon$ -approximability constants depending only on structural constants and $\varepsilon$ . Then $\omega_{L}\in A_{\infty}(\sigma)$ .

In particular, by Theorem 1.1 and Theorem 5.1, $\varepsilon$ -approximability of the subclass of solutions $u$ to $Lu=0$ in Theorem 5.1 is equivalent with $\varepsilon$ -approximability of all solutions $u$ to $Lu=0$ (and hence, it is equivalent with the other conditions in Corollary 1.3).

The proof of Theorem 5.1 is based on the proof of [CHMT20, Theorem 1.1], which itself is based on the techniques used in [KKPT00] and [KKPT16]. The key idea is the following. We fix a cube $Q_{0}$ and a Borel set $F\subset Q_{0}$ and we build a suitable solution $u=u_{F}$ associated to $F$ such that $u$ oscillates a significant amount. We then use this oscillation to control the $L^{1}$ norm of the gradient of an $\varepsilon$ -approximator of $u$ from below near $Q_{0}$ for a very small $\varepsilon$ . The $L^{1}$ -type Carleson measure estimate of the $\varepsilon$ -approximator then allows us to verify the $\omega_{L}\in A_{\infty}(\sigma)$ condition.

For the proof of Theorem 5.1, we need some definitions and notation from [CHMT20]. For clarity, we adopt most of the notation as it is from [CHMT20]. We define the following fattened version of the Whitney regions $U_{Q}$ and a “wider” version of the truncated dyadic cone:

\displaystyle U_{Q,\eta^{3}}\coloneqq\bigcup_{\begin{subarray}{c}Q^{\prime}\in\mathbb{D}_{Q}\\ \ell(Q^{\prime})>\eta^{3}\ell(Q)\end{subarray}}U_{Q^{\prime}},\qquad\text{ and }\qquad\Gamma_{Q_{0}}^{\eta}(x)\coloneqq\bigcup_{\begin{subarray}{c}Q\in\mathbb{D}_{Q_{0}}\\ Q\ni x\end{subarray}}U_{Q,\eta^{3}}.

The main difference between our approach and the proof of the implication “CME $\implies$ $A_{\infty}$ ” in [CHMT20] is the following lemma which is a modification of [CHMT20, Lemma 3.10]:

Lemma 5.2.

There exist $\varepsilon\in(0,1)$ , $0<\eta\ll 1$ , depending only on structural constants, and $\alpha_{0}\in(0,1)$ , $C_{\eta}\geq 1$ , both depending on structural constants and on $\eta$ , such that for each $Q_{0}\in\mathbb{D}$ , for every $\alpha\in(0,\alpha_{0})$ , and for every $F\subset Q_{0}$ satisfying $\omega_{L}^{X_{Q_{0}}}(F)\leq\alpha\omega_{L}^{X_{Q_{0}}}(Q_{0})$ , there exists a Borel set $S\subset Q_{0}$ such that if $u(X)=\omega_{L}^{X}(S)$ and $\Phi=\Phi^{\varepsilon}$ is an $\varepsilon$ -approximator of $u$ , then

\displaystyle\iint_{\Gamma_{Q_{0}}^{\eta}(y)}|\nabla\Phi(Y)|\delta(Y)^{-n}\,dY\geq C_{\eta}^{-1}\log(\alpha^{-1}),\qquad\text{for each }y\in F.

Before proving Lemma 5.2, let us see how it gives us Theorem 5.1.

Proof of Theorem 5.1.

Following [CHMT20]⁷⁷7Although the results of [CHMT20] are stated only for $n\geq 2$ , the proof of “(a) $\implies$ (b)” in [CHMT20, Theorem 1.1] works for a large part also for $n=1$ ; see [FP22]., we show that for each $\beta\in(0,1)$ , there exists $\alpha\in(0,1)$ such that for every $Q_{0}\in\mathbb{D}$ and every Borel set $F\subset Q_{0}$ , we have that

(5.3)

\frac{\omega_{L}^{X_{Q_{0}}}(F)}{\omega_{L}^{X_{Q_{0}}}(Q_{0})}\leq\alpha\quad\implies\quad\frac{\sigma(F)}{\sigma(Q_{0})}\leq\beta,

where the constants $\alpha$ and $\beta$ are independent of the choice of the dyadic system $\mathbb{D}$ . This is a dyadic version of the $A_{\infty}(\sigma)$ condition. Although it looks different than Definition 2.31, the conditions are equivalent since we consider doubling measures (see, for example, [GR85, Chapter IV, Theorem 2.11]) and the constants are independent of the system $\mathbb{D}$ (see [CHMT20, pp. 16] or use adjacent dyadic techniques [HK12, HT14]).

Fix $\beta\in(0,1)$ and $Q_{0}\in\mathbb{D}$ . Moreover, fix $\eta\in(0,1)$ small enough, and constants $\varepsilon$ , $\alpha_{0}$ , and $C_{\eta}$ as in Lemma 5.2. Let $F\subset Q_{0}$ satisfy $\omega_{L}^{X_{0}}(F)\leq\alpha\omega_{L}^{X_{0}}(Q_{0})$ with $\alpha\in(0,\alpha_{0})$ . We now use Lemma 5.2 to see that there exists $S\subset Q_{0}$ such that if $\Phi^{\varepsilon}$ is an $\varepsilon$ -approximator of $u(X)=\omega_{L}^{X}(S)$ , then

C_{\eta}^{-1}\log(\alpha^{-1})\sigma(F)\leq\int_{F}\iint_{\Gamma_{Q_{0}}^{\eta}(y)}|\nabla\Phi^{\varepsilon}(Y)|\delta(Y)^{-n}\,dY\,d\sigma(y)\leq C\eta^{-3n}\sigma(Q_{0}),

where the last estimate follows by Fubini’s theorem and property ii) of $\varepsilon$ -approximability in Theorem 1.1 (see [CHMT20, pp. 15–16] for more details). We may then choose $\alpha$ small enough depending on $\beta$ so that (5.3) holds. ∎

We turn to the proof of Lemma 5.2. For this, we need the following machinery:

Definition 5.4.

Let $Q_{0}\in\mathbb{D}$ be a dyadic cube, $\mu$ a regular Borel measure on $Q_{0}$ , $F\subset Q_{0}$ be a Borel set, $\epsilon_{0}>0$ , and $k\in\mathbb{N}$ be fixed. We say that a collection of nested Borel subsets $\{\mathcal{O}_{\ell}\}_{\ell=1}^{k}$ of $Q_{0}$ is a good $\epsilon_{0}$ -cover of $F$ of length $k$ for $\mu$ if

(a)

$F\subset\mathcal{O}_{k}\subset\mathcal{O}_{k-1}\subset\cdots\subset\mathcal{O}_{2}\subset\mathcal{O}_{1}\subset Q_{0}$ ,
(b)

$\mathcal{O}_{\ell}=\bigcup_{i}Q_{i}^{\ell}$ for disjoint subcubes $Q_{i}^{\ell}\in\mathbb{D}_{Q_{0}}$ ,
(c)

$\mu(\mathcal{O}_{\ell}\cap Q_{i}^{\ell-1})\leq\epsilon_{0}\mu(Q_{i}^{\ell-1})$ for every $i$ and every $2\leq\ell\leq k$ .

Lemma 5.5 ([CHMT20, Lemma 3.5]).

Fix $Q_{0}\in\mathbb{D}$ and suppose that $\epsilon_{0}\in(0,\tfrac{1}{e})$ . If $F\subset Q_{0}$ is a Borel set such that $\omega_{L}^{X_{Q_{0}}}(F)\leq\alpha\omega_{L}^{X_{Q_{0}}}(Q_{0})$ for $\alpha\in(0,\epsilon_{0}^{2}/C^{\prime})$ , then there exists a good $\epsilon_{0}$ -cover of $F$ of length $k\approx\tfrac{\log\alpha^{-1}}{\log\epsilon_{0}^{-1}}$ for $\omega_{L}^{X_{Q_{0}}}$ . Here $C^{\prime}$ depends only on the constant $C$ of Lemma 2.37.

Let $Q_{0}\in\mathbb{D}$ be a fixed cube and $F\subset Q_{0}$ be a fixed subset such that $\omega_{L}^{X_{0}}(F)\leq\alpha\omega_{L}^{X_{0}}(Q_{0})$ for $\alpha$ small enough. For each $Q\in\mathbb{D}$ , we let $\widetilde{Q}\in\mathbb{D}$ be the unique dyadic cube such that $x_{Q}\in\widetilde{Q}$ and $\ell(\widetilde{Q})=\eta\ell(Q)$ for $\eta>0$ a small enough parameter to be determined later. Fix $\epsilon_{0}>0$ small enough. Let $\{\mathcal{O}_{l}\}_{l=1}^{k}$ be a good $\epsilon_{0}$ -cover of $F$ given by Lemma 5.5. We set

\displaystyle\widetilde{\mathcal{O}}_{j}\coloneqq\bigcup_{i}\widetilde{Q}_{i}^{l}\quad\text{ and }\quad S\coloneqq\bigcup_{j=2}^{k}\widetilde{\mathcal{O}}_{j-1}\setminus\mathcal{O}_{j}

By the nestedness of the sets $\mathcal{O}_{j}$ , we have $\mathbbm{1}_{S}=\sum_{j=2}^{k}\mathbbm{1}_{\widetilde{\mathcal{O}}_{j-1}\setminus\mathcal{O}_{j}}$ . We define the nonnegative solution $u\coloneqq u_{F}\colon\Omega\to\mathbb{R}$ to $Lu=0$ as

\displaystyle u(X)=\int_{\partial\Omega}\mathbbm{1}_{S}(y)\,d\omega_{L}^{X}(y)=\omega_{L}^{X}(S)=\sum_{j=2}^{k}\omega_{L}^{X}(\widetilde{\mathcal{O}}_{j-1}\setminus\mathcal{O}_{j}).

We have the following lower oscillation bound:

Lemma 5.6 ([CHMT20, Lemma 3.24]).

There exists a structural constant $c_{0}>0$ such that if $\eta$ and $\epsilon_{0}=\epsilon_{0}(\eta,c_{0})$ are small enough, then for any $y\in F$ and any $1\leq l\leq k-1$ , there exist dyadic cubes $Q_{i}^{l}$ and $P_{i}^{l}$ such that

\displaystyle\big{|}u(\hat{X}_{\widetilde{Q}_{i}^{l}})-u(\hat{X}_{\widetilde{P}_{i}^{l}})\big{|}\geq c_{0},

where $Q_{i}^{l}$ is the unique cube (in $\mathcal{O}_{l}$ ) such that $y\in Q_{i}^{l}$ , $P_{i}^{l}\in\mathbb{D}_{Q_{i}^{l}}$ is the unique cube such that $y\in P_{i}^{l}$ and $\ell(P_{i}^{l})=\eta\ell(Q_{i}^{l})$ , $\widetilde{Q}_{i}^{l}$ and $\widetilde{P}_{i}^{l}$ are defined as we did after Lemma 5.5 and $\hat{X}_{\widetilde{Q}_{i}^{l}}$ and $\hat{X}_{\widetilde{P}_{i}^{l}}$ are corkscrew points relative to $x_{\widetilde{Q}_{i}^{l}}$ at scale $a_{1}\ell(\widetilde{Q}_{i}^{l})$ and relative to $x_{\widetilde{P}_{i}^{l}}$ at scale $a_{1}\ell(\widetilde{P}_{i}^{l})$ , respectively.

With the help of Lemma 5.6, we can now prove Lemma 5.2:

Proof of Lemma 5.2.

The proof is a straightforward modification of the proofs of [CHMT20, Lemma 3.10] and [HT21, Lemma 4.14]. Let us fix $y\in F$ and $l\in\{1,2,\ldots,k-1\}$ . We borrow the notation from Lemma 5.6, and set $\varepsilon=c_{0}/4$ . By adjusting the construction parameters for the Whitney regions in Section 2.3, we may assume that there exist Whitney cubes $I_{\widetilde{Q}_{i}^{l}}$ and $I_{\widetilde{P}_{i}^{l}}$ such that⁸⁸8In Section 2.3, we constructed the Whitney regions $U_{Q}$ in such a way that we know that a corkscrew point relative to $x_{Q}$ at scale $10^{-5}a_{0}\ell(Q)$ belongs to $U_{Q}$ . Corkscrew points of the type $\hat{X}_{Q}$ are also relative to $x_{Q}$ but they are at scale $a_{1}\ell(Q)$ .

(5.7)

\displaystyle\hat{X}_{\widetilde{Q}_{i}^{l}}\in I_{\widetilde{Q}_{i}^{l}}\subset(1+\tau)I_{\widetilde{Q}_{i}^{l}}\subset U_{\widetilde{Q}_{i}^{l}}\quad\text{ and }\quad\hat{X}_{\widetilde{P}_{i}^{l}}\in I_{\widetilde{P}_{i}^{l}}\subset(1+\tau)I_{\widetilde{P}_{i}^{l}}\subset U_{\widetilde{P}_{i}^{l}}.

Since $\ell(\widetilde{Q}_{i}^{l})=\eta\ell(Q_{i}^{l})$ and $\ell(\widetilde{P}_{i}^{l})=\eta\ell(P_{i}^{l})=\eta^{2}\ell(Q_{i}^{l})$ . In particular, since $\eta<1$ , we have $U_{\widetilde{Q}_{i}^{l}},U_{\widetilde{P}_{i}^{l}}\subset U_{Q_{i}^{l},\eta^{3}}$ .

Since $\Phi$ is a $\tfrac{c_{0}}{4}$ -approximator of $u$ , Lemma 5.6 gives us

c_{0}\leq|u(\hat{X}_{\widetilde{Q}_{i}^{l}})-u(\hat{X}_{\widetilde{P}_{i}^{l}})|\leq|u(\hat{X}_{\widetilde{Q}_{i}^{l}})-\Phi(\hat{X}_{\widetilde{Q}_{i}^{l}})|+|\Phi(\hat{X}_{\widetilde{Q}_{i}^{l}})-\Phi(\hat{X}_{\widetilde{P}_{i}^{l}})|+|\Phi(\hat{X}_{\widetilde{P}_{i}^{l}})-u(\hat{X}_{\widetilde{P}_{i}^{l}})|\\ \leq\frac{c_{0}}{2}+|\Phi(\hat{X}_{\widetilde{Q}_{i}^{l}})-\Phi(\hat{X}_{\widetilde{P}_{i}^{l}})|.

By (5.7), we know that the points $\hat{X}_{\widetilde{Q}_{i}^{l}}$ and $\hat{X}_{\widetilde{P}_{i}^{l}}$ are well inside $U_{Q_{i}^{l},\eta^{3}}$ in the sense that $\operatorname{dist}(\hat{X}_{\widetilde{P}_{i}^{l}},\partial U_{Q_{i}^{l},\eta^{3}})\approx\operatorname{dist}(\hat{X}_{\widetilde{Q}_{i}^{l}},\partial U_{Q_{i}^{l},\eta^{3}})\approx\ell(Q_{i}^{l})$ . Since $\operatorname{diam}(U_{Q_{i}^{l},\eta^{3}})\approx_{\eta}\ell(Q_{i}^{l})$ and $U_{Q_{i}^{l},\eta^{3}}$ satisfies the Harnack chain condition by Lemma 2.18, there exists a chain of $N=N(\eta)$ balls $B_{1},B_{2},\ldots,B_{N}$ such that

(i)

$\hat{X}_{\widetilde{Q}_{i}^{l}}\in B_{1}$ , $\hat{X}_{\widetilde{P}_{i}^{l}}\in B_{N}$ and $B_{j}\cap B_{j+1}\neq\mbox{{\O}}$ for every $j=1,2,\ldots,N-1$ ,
(ii)

the radii of the balls $B_{j}$ are comparable to $\ell(Q_{i}^{l})$ , depending on $\eta$ and $c_{0}$ ,
(iii)

$|\Phi(\hat{X}_{\widetilde{Q}_{i}^{l}})-\Phi(X)|<\tfrac{c_{0}}{8}$ for every $X\in B_{1}$ and $|\Phi(\hat{X}_{\widetilde{P}_{i}^{l}})-\Phi(Y)|<\tfrac{c_{0}}{8}$ for every $Y\in B_{N}$ ,
(iv)

for each $j=1,2,\ldots,N-1$ , there exists a cylinder $\mathfrak{C}_{j}$ connecting $B_{j}$ to $B_{j+1}$ such that $B_{j}\cup B_{j+1}\subset\mathfrak{C}_{j}\subset U_{Q_{i}^{l},\eta^{3}}$ , the cylinders have bounded overlaps and $|\mathfrak{C}_{j}|\approx_{\eta}\ell(Q_{i}^{l})^{n+1}$ ,

where the bound (iii) follows from properties of $\varepsilon$ -approximators (see property iii) in Theorem 1.1). Then

|\Phi(\hat{X}_{\widetilde{Q}_{i}^{l}})-\Phi(\hat{X}_{\widetilde{P}_{i}^{l}})|\\ \leq\mathchoice{{\vbox{\hbox{$\textstyle\raisebox{-8.61108pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-7.09656pt}}{{\vbox{\hbox{$\scriptstyle\raisebox{-6.02777pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-6.09258pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-4.29826pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-3.54826pt}}\!\iint_{B_{1}}|\Phi(\hat{X}_{\widetilde{Q}_{i}^{l}})-\Phi(X)|\,dX+\sum_{j=1}^{N-1}\Big{|}\mathchoice{{\vbox{\hbox{$\textstyle\raisebox{-8.61108pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-7.09656pt}}{{\vbox{\hbox{$\scriptstyle\raisebox{-6.02777pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-6.09258pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-4.29826pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-3.54826pt}}\!\iint_{B_{j}}\Phi(X)\,dX-\mathchoice{{\vbox{\hbox{$\textstyle\raisebox{-8.61108pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-7.09656pt}}{{\vbox{\hbox{$\scriptstyle\raisebox{-6.02777pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-6.09258pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-4.29826pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-3.54826pt}}\!\iint_{B_{j+1}}\Phi(X)\,dX\Big{|}+\mathchoice{{\vbox{\hbox{$\textstyle\raisebox{-8.61108pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-7.09656pt}}{{\vbox{\hbox{$\scriptstyle\raisebox{-6.02777pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-6.09258pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-4.29826pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-3.54826pt}}\!\iint_{B_{N}}|\Phi(X)-\Phi(\hat{X}_{\widetilde{P}_{i}^{l}})|\,dX\\ \leq\frac{c_{0}}{4}+\sum_{j=1}^{N-1}\Big{|}\mathchoice{{\vbox{\hbox{$\textstyle\raisebox{-8.61108pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-7.09656pt}}{{\vbox{\hbox{$\scriptstyle\raisebox{-6.02777pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-6.09258pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-4.29826pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-3.54826pt}}\!\iint_{B_{j}}\Phi(X)\,dX-\mathchoice{{\vbox{\hbox{$\textstyle\raisebox{-8.61108pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-7.09656pt}}{{\vbox{\hbox{$\scriptstyle\raisebox{-6.02777pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-6.09258pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-4.29826pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-3.54826pt}}\!\iint_{B_{j+1}}\Phi(X)\,dX\Big{|}

and furthermore, by Poincaré inequality and properties of the Whitney regions and cylinders $\mathfrak{C}_{j}$ , we have

\Big{|}\mathchoice{{\vbox{\hbox{$\textstyle\raisebox{-8.61108pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-7.09656pt}}{{\vbox{\hbox{$\scriptstyle\raisebox{-6.02777pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-6.09258pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-4.29826pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-3.54826pt}}\!\iint_{B_{j}}\Phi(X)\,dX-\mathchoice{{\vbox{\hbox{$\textstyle\raisebox{-8.61108pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-7.09656pt}}{{\vbox{\hbox{$\scriptstyle\raisebox{-6.02777pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-6.09258pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-4.29826pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-3.54826pt}}\!\iint_{B_{j+1}}\Phi(X)\,dX\Big{|}=\Big{|}\mathchoice{{\vbox{\hbox{$\textstyle\raisebox{-8.61108pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-7.09656pt}}{{\vbox{\hbox{$\scriptstyle\raisebox{-6.02777pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-6.09258pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-4.29826pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-3.54826pt}}\!\iint_{B_{j}}\Phi(X)\,dX-\mathchoice{{\vbox{\hbox{$\textstyle\raisebox{-8.61108pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-7.09656pt}}{{\vbox{\hbox{$\scriptstyle\raisebox{-6.02777pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-6.09258pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-4.29826pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-3.54826pt}}\!\iint_{\mathfrak{C}_{j}}\Phi(X)\,dX+\mathchoice{{\vbox{\hbox{$\textstyle\raisebox{-8.61108pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-7.09656pt}}{{\vbox{\hbox{$\scriptstyle\raisebox{-6.02777pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-6.09258pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-4.29826pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-3.54826pt}}\!\iint_{\mathfrak{C}_{j}}\Phi(X)\,dX-\mathchoice{{\vbox{\hbox{$\textstyle\raisebox{-8.61108pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-7.09656pt}}{{\vbox{\hbox{$\scriptstyle\raisebox{-6.02777pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-6.09258pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-4.29826pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-3.54826pt}}\!\iint_{B_{j+1}}\Phi(X)\,dX\Big{|}\\ \lesssim_{\eta}\mathchoice{{\vbox{\hbox{$\textstyle\raisebox{-8.61108pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-7.09656pt}}{{\vbox{\hbox{$\scriptstyle\raisebox{-6.02777pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-6.09258pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-4.29826pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-3.54826pt}}\!\iint_{\mathfrak{C}_{j}}\Big{|}\Phi(X)-\mathchoice{{\vbox{\hbox{$\textstyle\raisebox{-8.61108pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-7.09656pt}}{{\vbox{\hbox{$\scriptstyle\raisebox{-6.02777pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-6.09258pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-4.29826pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-3.54826pt}}\!\iint_{\mathfrak{C}_{j}}\Phi(Y)\,dY\Big{|}\,dX\lesssim_{\eta}\frac{\ell(Q_{i}^{l})}{|\mathfrak{C}_{j}|}\iint_{\mathfrak{C}_{j}}|\nabla\Phi(X)|\,dX\lesssim_{\eta}\iint_{\mathfrak{C}_{j}}|\nabla\Phi(X)|\delta(X)^{-n}\,dX

for every $j=1,2,\ldots,N-1$ . Combining the previous estimates and using the bounded overlaps of the cylinders then gives us

\displaystyle\frac{c_{0}}{4}\leq\sum_{j=1}^{N-1}\Big{|}\mathchoice{{\vbox{\hbox{$\textstyle\raisebox{-8.61108pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-7.09656pt}}{{\vbox{\hbox{$\scriptstyle\raisebox{-6.02777pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-6.09258pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-4.29826pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-3.54826pt}}\!\iint_{B_{j}}\Phi(X)\,dX-\mathchoice{{\vbox{\hbox{$\textstyle\raisebox{-8.61108pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-7.09656pt}}{{\vbox{\hbox{$\scriptstyle\raisebox{-6.02777pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-6.09258pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-4.29826pt}}{{\vbox{\hbox{$\scriptscriptstyle\raisebox{-4.30554pt}{\rotatebox[origin={c}]{20.0}{${-}\mkern-2.0mu{-}$}}$}}\kern-3.54826pt}}\!\iint_{B_{j+1}}\Phi(X)\,dX\Big{|}

\displaystyle\lesssim_{\eta}\iint_{U_{Q_{i}^{l},\eta^{3}}}|\nabla\Phi(X)|\delta(X)^{-n}\,dX.

Finally, we sum over $l$ and use Lemma 5.5, the bounded overlaps of the regions $U_{Q_{i}^{l},\eta^{3}}$ and the structure of $\Gamma_{Q_{i}^{l}}^{\eta}$ to get

\frac{c_{0}}{4}\frac{\log\alpha^{-1}}{\log\epsilon_{0}^{-1}}\approx\frac{c_{0}}{4}(k-1)\lesssim_{\eta}\sum_{l=1}^{k-1}\iint_{U_{Q_{i}^{l},\eta^{3}}}|\nabla\Phi(X)|\delta(X)^{-n}\,dX\lesssim_{\eta}\iint_{\Gamma_{Q_{i}^{l}}^{\eta}(y)}|\nabla\Phi(X)|\delta(X)^{-n}\,dX,

which proves the claim. ∎

6. Proof of Theorem 1.5

With the help of the tools from the previous sections, we can now prove Theorem 1.5. The proof is an adaptation of the corresponding proof from [HT21] (which itself uses some core ideas of Varopoulos [Var77, Var78] and Garnett [Gar07]).

Proof of Theorem 1.5.

Let $f\in\operatorname{BMO}_{\operatorname{c}}(\partial\Omega)$ and $\mathbb{D}$ be a dyadic system on $\partial\Omega$ . By Lemma 10.1 and Remark 10.3 in [HT21], we know that

(6.1)

\displaystyle f=f_{0}+g,

where $f_{0}\in L^{\infty}(\partial\Omega)$ with $\|f_{0}\|_{L^{\infty}(\partial\Omega)}\lesssim\|f\|_{\operatorname{BMO}}$ and $g=\sum_{j}\alpha_{j}\mathbbm{1}_{Q_{j}}$ , for a collection of dyadic cubes $\mathcal{A}\coloneqq\{Q_{j}\}_{j}\subset\mathbb{D}$ satisfying a Carleson packing condition with $\mathcal{C}_{\mathcal{A}}\lesssim 1$ and coefficients $\alpha_{j}$ satisfying $\sup_{j}|\alpha_{j}|\lesssim\|f\|_{\operatorname{BMO}(\partial\Omega)}$ . By [HT21, Proposition 1.3], there exists an extension $G$ of $g$ in $\Omega$ that satisfies the corresponding versions of the properties (1) – (3) in Theorem 1.5. Thus, it is enough to construct the extension for the function $f_{0}$ in the decomposition (6.1).

Construction of the extension for $f_{0}$ follows the proof of [HT21, Theorem 1.2]. We repeatedly use Lemma 3.3 to solve boundary value problems for $L$ with updated data and Theorem 1.1 to take smooth $\tfrac{1}{2}$ -approximators (that is, $\varepsilon$ -approximators for $\varepsilon=\tfrac{1}{2}$ ) for the corresponding solutions and to take their a.e. non-tangential boundary traces. The idea is that a suitable approximator is a Varopoulos-type extension plus an error term and we can keep on halving the size of the error term through iteration.

We start by taking the solution $u_{0}$ to the boundary value problem with data $f_{0}$ , satisfying $u_{0}(X)=\int_{\partial\Omega}f_{0}\,d\omega^{X}$ . We then take the $\tfrac{1}{2}$ -approximator $\Phi_{0}$ to $u_{0}$ . This approximator has an a.e. non-tangential boundary trace $\varphi_{0}$ . By Theorem 1.1, these functions satisfy $\|u_{0}-\Phi_{0}\|_{L^{\infty}(\Omega)}\leq\tfrac{1}{2}\|u_{0}\|_{L^{\infty}(\Omega)}\leq\tfrac{1}{2}\|f_{0}\|_{L^{\infty}(\partial\Omega)}$ , and

\displaystyle\sup_{x\in\partial\Omega,r>0}\frac{1}{r^{n}}\iint_{B(x,r)\cap\Omega}|\nabla\Phi_{0}(Y)|\,dY\leq\widetilde{C}\|u_{0}\|_{L^{\infty}(\Omega)}\leq\widetilde{C}\|f_{0}\|_{L^{\infty}(\partial\Omega)}.

We set $f_{1}\coloneqq f_{0}-\varphi_{0}$ which is the a.e. non-tangential boundary trace of $u_{0}-\Phi_{0}$ . We take the solution $u_{1}$ with the data $f_{1}$ , satisfying $u_{1}(X)=\int_{\partial\Omega}f_{1}\,d\omega^{X}$ . This solution satisfies

\displaystyle\|u_{1}\|_{L^{\infty}(\Omega)}\leq\|f_{1}\|_{L^{\infty}(\partial\Omega)}\leq\|u_{0}-\Phi_{0}\|_{L^{\infty}(\Omega)}\leq\tfrac{1}{2}\|u_{0}\|_{L^{\infty}(\Omega)}\leq\tfrac{1}{2}\|f_{0}\|_{L^{\infty}(\partial\Omega)},

and thus, we can take a $\tfrac{1}{2}$ -approximator $\Phi_{1}$ of $u_{1}$ , with boundary trace $\varphi_{1}$ . We get

\displaystyle\|u_{1}-\Phi_{1}\|_{L^{\infty}(\Omega)}\leq\tfrac{1}{2}\|u_{1}\|_{L^{\infty}(\Omega)}\leq\tfrac{1}{4}\|u_{0}\|_{L^{\infty}(\Omega)}

and

\displaystyle\sup_{x\in\partial\Omega,r>0}\frac{1}{r^{n}}\iint_{B(x,r)\cap\Omega}|\nabla\Phi_{1}(Y)|\,dY\leq\widetilde{C}\|u_{1}\|_{L^{\infty}(\Omega)}\leq\frac{\widetilde{C}}{2}\|f_{0}\|_{L^{\infty}(\partial\Omega)}.

We set $f_{2}\coloneqq f_{1}-\varphi_{1}=f_{0}-\varphi_{0}-\varphi_{1}$ , and continue in the previous way. This gives us a sequences of solutions $u_{k}$ and their $\tfrac{1}{2}$ -approximators $\Phi_{k}$ . The functions $u_{k}$ and $\Phi_{k}$ have a.e. non-tangential boundary traces $f_{k}$ and $\varphi_{k}$ , respectively, and we have

(i)

$f_{k+1}=f_{0}-\sum_{i=0}^{k}\varphi_{i}$ ,
(ii)

$\|f_{k+1}\|_{L^{\infty}(\partial\Omega)}\leq\|u_{k}-\Phi_{k}\|_{L^{\infty}(\Omega)}\leq 2^{-k-1}\|u_{0}\|_{L^{\infty}(\Omega)}\leq 2^{-k-1}\|f_{0}\|_{L^{\infty}(\partial\Omega)}$ ,
(iii)

$\|u_{k}\|_{L^{\infty}(\Omega)}\leq\|f_{k}\|_{L^{\infty}(\partial\Omega)}\leq 2^{-k}\|u_{0}\|_{L^{\infty}(\Omega)}$ ,
(iv)

$\sup_{x\in\partial\Omega,r>0}\frac{1}{r^{n}}\iint_{B(x,r)\cap\Omega}|\nabla\Phi_{k}(Y)|\,dY\leq\widetilde{C}\|u_{k}\|_{L^{\infty}(\Omega)}\leq\widetilde{C}2^{-k}\|u_{0}\|_{L^{\infty}(\Omega)}$ , and
(v)

$\|\Phi_{k}\|_{L^{\infty}(\Omega)}\lesssim 2^{-k}\|u_{0}\|_{L^{\infty}(\Omega)}$ .

Property (v) follows from the properties (ii) and (iii) combined with the fact that $\Phi_{k}$ is a $\tfrac{1}{2}$ -approximator of $u_{k}$ . By this property, we can define a uniformly convergent series

\displaystyle\Phi(X)\coloneqq\sum_{k=0}^{\infty}\Phi_{k}(X)

for $X\in\Omega$ . Let $F_{0}$ be a version of $\Phi$ that has been smoothened using similar convolution techniques as in the proof of Lemma 4.28 (see [HT21, Section 3] for details). By the arguments in the proof of [HT21, Theorem 1.1], $F_{0}$ is an extension of $f_{0}$ that satisfies the properties (1)–(3) in Theorem 1.5. Thus, by the decomposition (6.1), we may define the extensions $F$ in Theorem 1.5 by setting $F\coloneqq F_{0}+G$ . This completes the proof. ∎

7. An example in $\mathbb{R}^{3}$

In this section we construct a three-dimensional version of the example provided in Corollary 1.7. That is, we show there is a domain $\Omega$ in $\mathbb{R}^{3}$ whose boundary is not rectifiable and such that every function $f\in\operatorname{BMO}(\partial\Omega)$ has a Varopoulos extension in $\Omega$ .

To define $\Omega$ , denote by $E$ the $4$ -corner Cantor set in $\mathbb{R}^{2}$ . That is $E=\bigcap_{k=0}^{\infty}E_{k}$ , where $E_{k}$ equals the union of $4^{k}$ closed squares $Q_{i}^{k}$ of side length $4^{-k}$ located in the corners of the squares $Q_{j}^{k-1}$ of the previous generation (see [DM21, Section 3] for the precise definition). We assume that the center of $E_{0}=Q_{1}^{0}$ coincides with the origin in $\mathbb{R}^{2}$ . We consider the half-plane $\Pi=\{(x,y)\in\mathbb{R}^{2}:y>-2\}$ and we set $V=\Pi\setminus E$ . We also write $\mathbb{L}=\partial\Pi=\{(x,y)\in\mathbb{R}^{2}:y=-2\}$ . Then we define $\Omega=V\times\mathbb{R}$ . It is straightforward to check that both $V$ and $\Omega$ are uniform domains.

Notice that $\partial V=E\cup\mathbb{L}$ is $1$ -Ahlfors regular, while $\partial\Omega=(E\cup\mathbb{L})\times\mathbb{R}$ is $2$ -Ahlfors regular. In fact, the purpose of introducing the half plane $\Pi$ and the line $\mathbb{L}$ in this previous construction is to ensure the $2$ -Ahlfors regularity of $\partial\Omega$ . It is also clear that $E\times\mathbb{R}$ is purely $2$ -unrectifiable (since this set has no approximate tangent planes at any point).

Let $A$ be the $2\times 2$ matrix in the David–Mayboroda example in Theorem 1.6, let $L=-\operatorname{div}A\nabla$ be the associated elliptic operator in $\mathbb{R}^{2}$ , and let

\hat{A}=\left(\begin{array}[]{ll}A&0\\ 0&1\end{array}\right).

Set $\hat{L}=-\mathop{\operatorname{div}}\nolimits\hat{A}\,\nabla$ in $\mathbb{R}^{3}$ . Below we will show that $\omega_{\hat{L},\Omega}\in A_{\infty}(\sigma)$ , where $\sigma$ is the surface measure on $\partial\Omega$ . Consequently, by Theorem 1.5, every function $f\in\operatorname{BMO}(\partial\Omega)$ has a Varopoulos extension in $\Omega$ .

First, we prove the following.

Lemma 7.1.

We have that $\omega_{L,V}\in A_{\infty}(\mathcal{H}^{1}|_{\partial V})$ . Further, there is a constant $C>0$ such that for any surface ball $\Delta\subset\partial V$ and any corkscrew point $p\in V$ for $\Delta$ ,

(7.2)

\frac{d\omega_{L,V}^{p}}{d\mathcal{H}^{1}|_{\partial V}}(x)\leq\frac{C}{\mathcal{H}^{1}(\Delta)},\qquad\text{for each }x\in\Delta.

Proof.

It is clear that the estimate (7.2) implies the local $A_{\infty}$ condition of $\omega_{L}\equiv\omega_{L,V}$ . First we will show that, for any $p\in\Pi$ ,

(7.3)

\frac{d\omega_{L,\Pi}^{p}}{d\mathcal{H}^{1}|_{\mathbb{L}}}(x)\leq C\,\frac{d\omega_{-\Delta,\Pi}^{p}}{d\mathcal{H}^{1}|_{\mathbb{L}}}(x),

where $\omega_{L,\Pi}$ stands for the $L$ -elliptic measure for the domain $\Pi$ , and $\omega_{-\Delta,\Pi}$ is the harmonic measure (i.e. for the Laplacian) for the domain $\Pi$ . To this end, we consider the auxiliary domain $U=\Pi\setminus\overline{B_{1}}$ , where we denoted $B_{r}=B(0,r)$ . Observe that $A$ is the identity matrix on $U$ , and thus $\omega_{L,U}^{X}\equiv\omega_{-\Delta,U}^{X}$ for all $X\in U$ . Let $F\subset\mathbb{L}$ be an arbitrary closed subset, and let $q\in\partial B_{3/2}$ be such that

\omega_{L,\Pi}^{q}(F)=\max_{X\in\partial B_{3/2}}\omega_{L,\Pi}^{X}(F).

Since $\omega_{L,\Pi}^{x}(F)$ is a harmonic function of $x$ in $U$ , we have

(7.4)

\omega_{L,\Pi}^{q}(F)=\int_{\partial U}\omega^{X}_{L,\Pi}(F)\,d\omega^{q}_{-\Delta,U}(X)=\int_{\mathbb{L}}\omega_{L,\Pi}^{X}(F)\,d\omega^{q}_{-\Delta,U}(X)+\int_{\partial B_{1}}\omega_{L,\Pi}^{X}(F)\,d\omega^{q}_{-\Delta,U}(X)\\ \leq\omega_{-\Delta,U}^{q}(F)+\sup_{X\in\partial B_{1}}\omega_{L,\Pi}^{X}(F)\,\omega_{-\Delta,U}^{q}(\partial B_{1}).

By the maximum principle and the definition of $q$ , we have

\sup_{X\in\partial B_{1}}\omega_{L,\Pi}^{X}(F)\leq\sup_{X\in\partial B_{3/2}}\omega_{L,\Pi}^{X}(F)=\omega_{L,\Pi}^{q}(F).

Also, it is immediate that $c_{B}\coloneqq\omega_{-\Delta,U}^{q}(\partial B_{1})<1$ . Hence,

\omega_{L,\Pi}^{q}(F)\leq\omega_{-\Delta,U}^{q}(F)+c_{B}\,\omega_{L,\Pi}^{q}(F),

or equivalently, $\omega_{L,\Pi}^{q}(F)\leq(1-c_{B})^{-1}\omega_{-\Delta,U}^{q}(F)$ . By a Harnack chain argument we deduce that $\omega_{L,\Pi}^{p}(F)\lesssim\omega_{-\Delta,U}^{p}(F)=\omega_{L,U}^{p}(F)$ for all $p\in\partial B_{3/2}$ , and then by the maximum principle, it follows that the same estimate is valid for all $p\in\Pi\setminus\overline{B}_{3/2}$ . Using again the maximum principle, we get, for all $p\in\Pi\setminus\overline{B}_{3/2}$ ,

\omega_{L,\Pi}^{p}(F)\lesssim\omega_{-\Delta,U}^{p}(F)\leq\omega_{-\Delta,\Pi}^{p}(F).

Since $\omega_{L,\Pi}^{p}(F)$ and $\omega_{-\Delta,U}^{p}(F)$ are, respectively, elliptic and harmonic in $B_{3/2}$ , by a Harnack chain argument it follows that the estimate above also holds for all $p\in B_{3/2}$ , possibly with a different implicit constant. This is equivalent to (7.3).

To prove (7.2), let $B$ be a ball centered in $\partial V$ such that $\Delta=B\cap\partial\Omega$ and consider an arbitrary closed set $F\subset\Delta$ . Let $p\in B\cap V$ be a corkscrew point for $\Delta$ . By the maximum principle and by (7.3),

\omega_{L,V}^{p}(F\cap\mathbb{L})\leq\omega_{L,\Pi}^{p}(F\cap\mathbb{L})\lesssim\omega_{-\Delta,\Pi}^{p}(F\cap\mathbb{L})\lesssim\frac{\mathcal{H}^{1}(F\cap\mathbb{L})}{\mathcal{H}^{1}(\Delta)},

where in the last estimate we used that $\frac{d\omega_{-\Delta,\Pi}^{p}}{d\mathcal{H}^{1}|_{\mathbb{L}}}(x)\lesssim\sigma(\Delta)^{-1}$ for all $x\in\Delta$ . Indeed, if $F\cap\mathbb{L}=\varnothing$ , then there is nothing to show; if $B$ is centered on $\mathbb{L}$ , this follows from Ahlfors regularity of $\partial V$ and classical properties of the harmonic measure, and if $B$ is centered on $E$ and there exists $z\in F\cap\mathbb{L}$ , then if $B^{\prime}=B(z,\operatorname{rad}(B))$ and $p^{\prime}$ is a corkscrew point for $B^{\prime}$ in $V$ , then by Harnack chains we have that $\omega^{p}_{-\Delta,\Pi}(F\cap\mathbb{L})\approx\omega^{p^{\prime}}_{-\Delta,\Pi}(F\cap\mathbb{L})$ , and the claim follows as in the previous case.

To estimate $\omega_{L,V}^{p}(F\cap E)$ , we assume that $B\cap E\neq\varnothing$ and we distinguish two cases. Suppose first that $r(B)$ , the radius of $B$ , satisfies $r(B)\leq 1$ . Denote $X_{E}=(0,2)$ and let $B^{\prime}$ be a ball centered in $E$ containing $B$ with radius at most $2r(B)$ . By the maximum principle, the change of poles formula, and the fact that $\frac{d\omega_{L,E^{c}}^{X_{E}}}{d\mathcal{H}^{1}|_{E}}(x)\approx 1$ (by Theorem 1.6), we obtain

\omega_{L,V}^{p}(F\cap E)\leq\omega_{L,E^{c}}^{p}(F\cap E)\approx\frac{\omega_{L,E^{c}}^{X_{E}}(F\cap E)}{\omega_{L,E^{c}}^{X_{E}}(B^{\prime})}\approx\frac{\mathcal{H}^{1}(F\cap E)}{\mathcal{H}^{1}(B^{\prime}\cap E)}\approx\frac{\mathcal{H}^{1}(F\cap E)}{\mathcal{H}^{1}(\Delta)}.

In the case $r(B)>1$ , recall that $B_{1}=B(0,1)$ , and then using the change of poles formula and the maximum principle, we write

\omega_{L,V}^{p}(F\cap E)\approx\omega_{L,V}^{X_{E}}(F\cap E)\,\omega_{L,V}^{p}(B_{1})\leq\omega_{L,E^{c}}^{X_{E}}(F\cap E)\,\omega_{L,V}^{p}(B_{1})\approx\mathcal{H}^{1}(F\cap E)\,\omega_{L,V}^{p}(B_{1}).

To estimate $\omega_{L,V}^{p}(B_{1})$ , consider a ball $B_{1}^{\prime}$ of radius $1$ centered in $\mathbb{L}$ , disjoint from $E$ , and contained in $4B_{1}$ . We claim that

(7.5)

\omega_{L,V}^{p}(B_{1})\approx\omega_{L,V}^{p}(B_{1}^{\prime}).

Indeed, if $p\notin 8B_{1}$ , then (7.5) follows directly from Lemma 2.37. On the other hand, if $p\in 8B_{1}$ , then denote $p^{\prime}=(0,16)$ and note that $\delta(p)\gtrsim 1$ , $\delta(p^{\prime})\geq 8$ , and $|p-p^{\prime}|\leq 24$ . Then, by the Harnack inequality, Harnack chains, and Lemma 2.37, the estimate (7.5) follows. Next, by the maximum principle, (7.3), and Ahlfors regularity, we get

\omega_{L,V}^{p}(B_{1}^{\prime})\leq\omega_{L,\Pi}^{p}(B_{1}^{\prime})\lesssim\omega_{-\Delta,\Pi}^{p}(B_{1}^{\prime})\lesssim\frac{\mathcal{H}^{1}(B_{1}^{\prime}\cap\mathbb{L})}{\mathcal{H}^{1}(\Delta)}\approx\frac{1}{\mathcal{H}^{1}(\Delta)}.

Therefore, again we derive

\omega_{L,V}^{p}(F\cap E)\lesssim\frac{\mathcal{H}^{1}(F\cap E)}{\mathcal{H}^{1}(\Delta)}.

Altogether, we deduce that

\omega_{L,V}^{p}(F)=\omega_{L,V}^{p}(F\cap\mathbb{L})+\omega_{L,V}^{p}(F\cap E)\lesssim\frac{\mathcal{H}^{1}(F)}{\mathcal{H}^{1}(\Delta)}

for any closed set $F\subset\Delta$ , which is equivalent to the statement in the lemma, by the Lebesgue-Radon-Nykodim Theorem. ∎

Now we are ready to prove the $A_{\infty}$ property of the elliptic measure $\omega_{\hat{L}}$ for the three-dimensional domain $\Omega$ :

Proposition 7.6.

The elliptic measure $\omega_{\hat{L}}$ for the domain $\Omega\subset\mathbb{R}^{3}$ defined above satisfies the $A_{\infty}$ condition with respect to $\mathcal{H}^{2}|_{\partial\Omega}$ . Further, there is a constant $C>0$ such that for any surface ball $\Delta\subset\partial\Omega$ and any corkscrew point $p\in\Omega$ for $\Delta$ ,

(7.7)

\frac{d\omega_{\hat{L}}^{p}}{d\mathcal{H}^{2}|_{\partial\Omega}}(x)\leq\frac{C}{\mathcal{H}^{2}(\Delta)},\qquad\text{for each }x\in\Delta.

Proof.

Let $B$ be a ball with radius $r(B)$ centered in $\partial\Omega$ such that $\Delta=B\cap\partial\Omega$ . Clearly, to prove the $A_{\infty}$ condition for $\omega_{\hat{L}}$ , it suffices to show (7.7). In turn, since $\hat{L}$ is symmetric, by a direct application of Lemma 2.38 and the Lebesgue-Radon-Nykodim Theorem, it is enough to prove that

(7.8)

G_{\hat{L}}(X,p)\lesssim\frac{\operatorname{dist}(X,\partial\Omega)}{r(B)^{2}}\quad\mbox{ for all $X\in B\cap\Omega\setminus B(p,\frac{1}{2}\operatorname{dist}(p,\partial\Omega))$,}

where $G_{\hat{L}}$ is the $\hat{L}$ -Green function for $\Omega$ .

Denote by $P$ the orthogonal projection of $\mathbb{R}^{3}$ onto $\mathbb{R}^{2}\equiv\mathbb{R}^{2}\times\{0\}$ . Let $p_{0}=P(p)$ and consider the function $u:\overline{\Omega}\setminus P^{-1}(\{p_{0}\})\to\mathbb{R}$ defined by

u(X)=G_{L}(P(X),P(p))=G_{L}(P(X),p_{0}),

where $G_{L}$ is the Green’s function for $L$ in the domain $V=P(\Omega)$ . It is immediate to check that $u$ is $\hat{L}$ -elliptic in $\Omega\setminus P^{-1}(\{p_{0}\})$ , and clearly it can be extended continuously by zero to the whole $\partial\Omega$ . Thus, by the boundary Harnack principle, choosing $p^{\prime}\in\partial B(p,\frac{1}{2}\operatorname{dist}(p,\partial\Omega))$ such that $P(p^{\prime})\in\partial B(p_{0},\frac{1}{2}\operatorname{dist}(p,\partial\Omega))\cap\mathbb{R}^{2}$ , we have

\frac{G_{\hat{L}}(X,p)}{G_{\hat{L}}(p^{\prime},p)}\approx\frac{u(X)}{u(p^{\prime})}=\frac{G_{L}(P(X),p_{0})}{G_{L}(P(p^{\prime}),p_{0})}\quad\mbox{ for all $X\in B\cap\Omega\setminus B(p,\frac{1}{2}\operatorname{dist}(p,\partial\Omega))$.}

Thus, for such points $X$ and by (2.33) and (2.34) applied both to $G_{\hat{L}}$ and $G_{L}$ ,

\frac{G_{\hat{L}}(X,p)}{|p^{\prime}-p|^{-1}}\approx\frac{G_{L}(P(X),p_{0})}{1}.

Thus, by Lemma 2.38 applied to $G_{L,V}$ , Harnack chains, and the Harnack inequality, we see that

(7.9)

G_{\hat{L}}(X,p)\approx\frac{G_{L,V}(P(X),p_{0})}{r(B)}\approx\frac{\omega_{L,V}^{p_{0}}(\Delta_{V,X})}{r(B)},

where $\Delta_{V,X}=B(P(X),2\operatorname{dist}(P(X),\partial V))$ . From (7.2) we infer that

(7.10)

\omega_{L,V}^{p_{0}}(\Delta_{V,X})\lesssim\frac{\mathcal{H}^{1}(\Delta_{V,X})}{\mathcal{H}^{1}(P(B)\cap\partial V)}\approx\frac{\operatorname{dist}(X,\partial\Omega)}{r(B)}.

From (7.9) and (7.10), we deduce (7.8), which concludes the proof. ∎

References

[AHMT19] M. Akman, S. Hofmann, J. M. Martell, and T. Toro. Perturbation of elliptic operators in 1-sided NTA domains satisfying the capacity density condition, 2019. Preprint, arxiv.1901.08261.
[ADFJM19] D. N. Arnold, G. David, M. Filoche, D. Jerison, and S. Mayboroda. Localization of eigenfunctions via an effective potential. Comm. Partial Differential Equations, 44(11):1186–1216, 2019.
[AGMT22] J. Azzam, J. Garnett, M. Mourgoglou, and X. Tolsa. Uniform rectifiability, elliptic measure, square functions, and $\varepsilon$ -approximability via an acf monotonicity formula. Int. Math. Res. Not. IMRN, 06 2022. rnab095.
[AHMMT20] J. Azzam, S. Hofmann, J. M. Martell, M. Mourgoglou, and X. Tolsa. Harmonic measure and quantitative connectivity: geometric characterization of the $L^{p}$ -solvability of the Dirichlet problem. Invent. Math., 222(3):881–993, 2020.
[BH20] S. Bortz and S. Hofmann. Quantitative Fatou theorems and uniform rectifiability. Potential Anal., 53(1):329–355, 2020.
[BT19] S. Bortz and O. Tapiola. $\varepsilon$ -approximability of harmonic functions in $L^{p}$ implies uniform rectifiability. Proc. Amer. Math. Soc., 147(5):2107–2121, 2019.
[CFMS81] L. Caffarelli, E. Fabes, S. Mortola, and S. Salsa. Boundary behavior of nonnegative solutions of elliptic operators in divergence form. Indiana Univ. Math. J., 30(4):621–640, 1981.
[CDMT22] M. Cao, O. Domínguez, J. M. Martell, and P. Tradacete. On the $A_{\infty}$ condition for elliptic operators in 1-sided nontangentially accessible domains satisfying the capacity density condition. Forum Math. Sigma, 10:Paper No. e59, 57, 2022.
[CHM] M. Cao, P. Hidalgo-Palencia, and J. Martell. Carleson measure estimates, corona decompositions, and perturbation of elliptic operators without connectivity. Preprint. arXiv:2202.06363.
[Car62] L. Carleson. Interpolations by bounded analytic functions and the corona problem. Ann. of Math. (2), 76:547–559, 1962.
[CHMT20] J. Cavero, S. Hofmann, J. M. Martell, and T. Toro. Perturbations of elliptic operators in 1-sided chord-arc domains. Part II: non-symmetric operators and Carleson measure estimates. Trans. Amer. Math. Soc., 373(11):7901–7935, 2020.
[Chr90] M. Christ. A $T(b)$ theorem with remarks on analytic capacity and the Cauchy integral. Colloq. Math., 60/61(2):601–628, 1990.
[CF74] R. R. Coifman and C. Fefferman. Weighted norm inequalities for maximal functions and singular integrals. Studia Math., 51:241–250, 1974.
[CW71] R. R. Coifman and G. Weiss. Analyse harmonique non-commutative sur certains espaces homogènes. Lecture Notes in Mathematics, Vol. 242. Springer-Verlag, Berlin-New York, 1971. Étude de certaines intégrales singulières.
[Dah80] B. E. J. Dahlberg. Approximation of harmonic functions. Ann. Inst. Fourier (Grenoble), 30(2):vi, 97–107, 1980.
[DJK84] B. E. J. Dahlberg, D. S. Jerison, and C. E. Kenig. Area integral estimates for elliptic differential operators with nonsmooth coefficients. Ark. Mat., 22(1):97–108, 1984.
[DS91] G. David and S. Semmes. Singular integrals and rectifiable sets in ${\bf R}^{n}$ : Beyond Lipschitz graphs. Astérisque, 193:152, 1991.
[DFM] G. David, J. Feneuil, and S. Mayboroda. Elliptic theory in domains with boundaries of mixed dimension. Preprint. May 2020. arXiv:2003.09037.
[DM21] G. David and S. Mayboroda. Good elliptic operators on Cantor sets. Adv. Math., 383:Paper No. 107687, 21, 2021.
[DS93] G. David and S. Semmes. Analysis of and on uniformly rectifiable sets, volume 38 of Mathematical Surveys and Monographs. American Mathematical Society, Providence, RI, 1993.
[DG57] E. De Giorgi. Sulla differenziabilità e l’analiticità delle estremali degli integrali multipli regolari. Mem. Accad. Sci. Torino. Cl. Sci. Fis. Mat. Nat. (3), 3:25–43, 1957.
[DM95] G. Dolzmann and S. Müller. Estimates for Green’s matrices of elliptic systems by $L^{p}$ theory. Manuscripta Math., 88(2):261–273, 1995.
[DK09] H. Dong and S. Kim. Green’s matrices of second order elliptic systems with measurable coefficients in two dimensional domains. Trans. Amer. Math. Soc., 361(6):3303–3323, 2009.
[EG15] L. C. Evans and R. F. Gariepy. Measure theory and fine properties of functions. Textbooks in Mathematics. CRC Press, Boca Raton, FL, revised edition, 2015.
[FS72] C. Fefferman and E. M. Stein. $H^{p}$ spaces of several variables. Acta Math., 129(3-4):137–193, 1972.
[FP22] J. Feneuil and B. Poggi. Generalized Carleson perturbations of elliptic operators and applications. Trans. Amer. Math. Soc., 375(11):7553–7599, 2022.
[GR85] J. García-Cuerva and J. L. Rubio de Francia. Weighted norm inequalities and related topics, volume 116 of North-Holland Mathematics Studies. North-Holland Publishing Co., Amsterdam, 1985. Notas de Matemática [Mathematical Notes], 104.
[GMT18] J. Garnett, M. Mourgoglou, and X. Tolsa. Uniform rectifiability from Carleson measure estimates and $\varepsilon$ -approximability of bounded harmonic functions. Duke Math. J., 167(8):1473–1524, 2018.
[Gar07] J. B. Garnett. Bounded analytic functions, volume 236 of Graduate Texts in Mathematics. Springer, New York, first edition, 2007.
[Gar22] J. B. Garnett. Carleson measure estimates and $\varepsilon$ -approximation for bounded harmonic functions, without Ahlfors regularity assumptions. Rev. Mat. Iberoam., 38(1):323–351, 2022.
[GW82] M. Grüter and K.-O. Widman. The Green function for uniformly elliptic equations. Manuscripta Math., 37(3):303–342, 1982.
[HKM93] J. Heinonen, T. Kilpeläinen, and O. Martio. Nonlinear potential theory of degenerate elliptic equations. Oxford Mathematical Monographs. The Clarendon Press, Oxford University Press, New York, 1993. Oxford Science Publications.
[HKMP15] S. Hofmann, C. Kenig, S. Mayboroda, and J. Pipher. Square function/non-tangential maximal function estimates and the Dirichlet problem for non-symmetric elliptic operators. J. Amer. Math. Soc., 28(2):483–529, 2015.
[HK07] S. Hofmann and S. Kim. The Green function estimates for strongly elliptic systems of second order. Manuscripta Math., 124(2):139–172, 2007.
[HL18] S. Hofmann and P. Le. BMO solvability and absolute continuity of harmonic measure. J. Geom. Anal., 28(4):3278–3299, 2018.
[HLMN17] S. Hofmann, P. Le, J. M. Martell, and K. Nyström. The weak- $A_{\infty}$ property of harmonic and $p$ -harmonic measures implies uniform rectifiability. Anal. PDE, 10(3):513–558, 2017.
[HMMTZ21] S. Hofmann, J. M. Martell, , S. Mayboroda, T. Toro, and Z. Zhao. Uniform rectifiability and elliptic operators satisfying a carleson measure condition. Geom. Funct. Anal., 31:325–401, 2021.
[HM14] S. Hofmann and J. M. Martell. Uniform rectifiability and harmonic measure I: Uniform rectifiability implies Poisson kernels in $L^{p}$ . Ann. Sci. Éc. Norm. Supér. (4), 47(3):577–654, 2014.
[HMM16] S. Hofmann, J. M. Martell, and S. Mayboroda. Uniform rectifiability, Carleson measure estimates, and approximation of harmonic functions. Duke Math. J., 165(12):2331–2389, 2016.
[HMT17] S. Hofmann, J. M. Martell, and T. Toro. $A_{\infty}$ implies NTA for a class of variable coefficient elliptic operators. J. Differ. Equations, 263(10):6147–6188, 2017.
[HMU14] S. Hofmann, J. M. Martell, and I. Uriarte-Tuero. Uniform rectifiability and harmonic measure, II: Poisson kernels in $L^{p}$ imply uniform rectifiability. Duke Math. J., 163(8):1601–1654, 2014.
[HT20] S. Hofmann and O. Tapiola. Uniform rectifability and $\varepsilon$ -approximability of harmonic functions in $L^{p}$ . Ann. Inst. Fourier (Grenoble), 70(4):1595–1638, 2020.
[HT21] S. Hofmann and O. Tapiola. Uniform rectifiability implies Varopoulos extensions. Adv. Math., 390:Paper No. 107961, 53, 2021.
[HK12] T. Hytönen and A. Kairema. Systems of dyadic cubes in a doubling metric space. Colloq. Math., 126(1):1–33, 2012.
[HR18] T. Hytönen and A. Rosén. Bounded variation approximation of $L_{p}$ dyadic martingales and solutions to elliptic equations. J. Eur. Math. Soc. (JEMS), 20(8):1819–1850, 2018.
[HT14] T. Hytönen and O. Tapiola. Almost Lipschitz-continuous wavelets in metric spaces via a new randomization of dyadic cubes. J. Approx. Theory, 185:12–30, 2014.
[KKPT16] C. Kenig, B. Kirchheim, J. Pipher, and T. Toro. Square functions and the $A_{\infty}$ property of elliptic measures. J. Geom. Anal., 26(3):2383–2410, 2016.
[KKPT00] C. Kenig, H. Koch, J. Pipher, and T. Toro. A new approach to absolute continuity of elliptic measure, with applications to non-symmetric equations. Adv. Math., 153(2):231–298, 2000.
[KN85] C. E. Kenig and W.-M. Ni. On the elliptic equation $Lu-k+K\,{\rm exp}[2u]=0$ . Ann. Scuola Norm. Sup. Pisa Cl. Sci. (4), 12(2):191–224, 1985.
[Mat95] P. Mattila. Geometry of sets and measures in Euclidean spaces, volume 44 of Cambridge Studies in Advanced Mathematics. Cambridge University Press, Cambridge, 1995. Fractals and rectifiability.
[Mat21] P. Mattila. Rectifiability; a survey. preprint, arXiv:2112.00540, 2021.
[Mey63] N. G. Meyers. An $L^{p}$ -estimate for the gradient of solutions of second order elliptic divergence equations. Ann. Scuola Norm. Sup. Pisa Cl. Sci. (3), 17:189–206, 1963.
[Mos61] J. Moser. On Harnack’s theorem for elliptic differential equations. Comm. Pure Appl. Math., 14:577–591, 1961.
[MPT] M. Mourgoglou, B. Poggi, and X. Tolsa. Solvability of the Poisson-Dirichlet problem with interior data in ${L}^{p^{\prime}}$ -carleson spaces and its applications to the ${L}^{p}$ -regularity problem. Preprint. January 2023 (v2). July 2022 (v1). arXiv: 2207.10554.
[Muc72] B. Muckenhoupt. Weighted norm inequalities for the Hardy maximal function. Trans. Amer. Math. Soc., 165:207–226, 1972.
[Nas58] J. Nash. Continuity of solutions of parabolic and elliptic equations. Amer. J. Math., 80:931–954, 1958.
[Ste70] E. M. Stein. Singular integrals and differentiability properties of functions. Princeton Mathematical Series, No. 30. Princeton University Press, Princeton, N.J., 1970.
[Var77] N. T. Varopoulos. BMO functions and the $\overline{\partial}$ -equation. Pacific J. Math., 71(1):221–273, 1977.
[Var78] N. T. Varopoulos. A remark on functions of bounded mean oscillation and bounded harmonic functions. Addendum to: “BMO functions and the $\overline{\partial}$ -equation” (Pacific J. Math. 71 (1977), no. 1, 221–273). Pacific J. Math., 74(1):257–259, 1978.

	$\displaystyle\|\Phi(X)-u(X)\|$	$\displaystyle=\Big{\|}\iint\Lambda_{\varepsilon}(X,Y)(\Phi_{0}(Y)-u(X))\,dY\Big{\|}$
		$\displaystyle<\varepsilon\\|u\\|_{L^{\infty}(\Omega)}+\iint\Lambda_{\varepsilon}(X,Y)\left\|u(Y)-u(X)\right\|\,dY$
		$\displaystyle\leq\varepsilon\\|u\\|_{L^{\infty}(\Omega)}+\iint\Lambda_{\varepsilon}(X,Y)\,C\Big{(}\frac{\|X-Y\|}{\tfrac{1}{4}\delta(X)}\Big{)}^{\alpha}\Big{(}\fint_{B(X,\tfrac{1}{2}\delta(X))}\|u(Z)\|^{2}\,dZ\Big{)}^{1/2}\,dY$
		$\displaystyle\leq\varepsilon\\|u\\|_{L^{\infty}(\Omega)}+C\left(2\xi_{\varepsilon}\right)^{\alpha}\\|u\\|_{L^{\infty}(\Omega)}\iint\Lambda_{\varepsilon}(X,Y)\,dY\leq 2\varepsilon\\|u\\|_{L^{\infty}(\Omega)}$

The A∞A_{\infty} condition, ε\varepsilon-approximators, and Varopoulos extensions in uniform domains

Abstract.

1. Introduction

Theorem 1.1.

Theorem 1.2.

Corollary 1.3.

Corollary 1.4.

Theorem 1.5.

Theorem 1.6 ([DM21, Section 4]).

Corollary 1.7.

2. Preliminaries

2.1. Uniform domains, chord-arc domains, Ahlfors regularity and uniform rectifiability

Definition 2.1 (Corkscrew condition).

Definition 2.2 (Harnack chain Condition).

Definition 2.3 (Uniform domain).

Definition 2.4 (Ahlfors regularity).

Definition 2.5 (Chord-arc domain).

Definition 2.6 (Uniform rectifiability).

2.2. Carleson measures, non-tangential convergence, BMO and local BV

Definition 2.7 (Carleson measures).

Definition 2.8 (Cones and non-tangential convergence).

Definition 2.10 (Non-tangential maximal operator).

Remark 2.11.

Definition 2.12 (BMO).

Definition 2.13 (Local BV).

2.3. Dyadic cubes, Whitney regions and approximating domains

Lemma 2.14 ([Chr90, DS91, HK12]).

Definition 2.15 (Carleson packing condition).

Lemma 2.18 ([HM14, Lemma 3.61]).

2.4. The choice of the parameter κ\kappa

Lemma 2.21.

Proof.

Lemma 2.25.

Proof.

Lemma 2.26.

Proof.

2.5. Elliptic PDE estimates

Lemma 2.27 (Caccioppoli Inequality).

Lemma 2.28 (Hölder continuity of solutions, [DG57, Nas58]).

Lemma 2.29 (Harnack inequality [Mos61]).

Definition 2.31 (A∞A_{\infty} for elliptic measure).

Definition 2.32 (Green’s function).

Lemma 2.36 (Bourgain’s estimate, [DFM, Lemma 15.1]).

Lemma 2.37 ((Local) doubling property, [DFM, Lemma 15.43]).

Lemma 2.38 (CFMS estimate, [DFM, Lemma 15.28]).

Lemma 2.39 (Boundary Harnack Principle, [DFM, Theorem 15.64]).

Lemma 2.40 (Change of Poles, [DFM, Lemma 15.61]).

Lemma 2.41 ([DFM, Lemma 11.32]).

3. Set-up for Theorem 1.1

3.1. CME and non-tangential convergence

Lemma 3.1 ([FP22, Corollary 1.32]).

Lemma 3.2 ([HL18, Theorem 1.3]).

Lemma 3.3.

Proof.

3.2. A few important lemmas

Lemma 3.10.

Proof.

Lemma 3.16.

Definition 3.17.

Lemma 3.18 ([AGMT22, Proposition 3.1]).

Proof.

Lemma 3.19.

Proof.

4. Existence of ε\varepsilon-approximators: proof of Theorem 1.1

Theorem 4.1.

4.1. Set-up for the proof of Theorem 4.1

Lemma 4.2.

Proof.

4.2. A stopping time decomposition for the family of blue cubes

Lemma 4.5.

Proof.

Lemma 4.6.

Lemma 4.9.

Proof of Lemma 4.6.

4.3. Proof of Lemma 4.9: A high oscillation estimate

Lemma 4.12.

Lemma 4.18.

Proof.

Lemma 4.23.

Proof.

The $A_{\infty}$ condition, $\varepsilon$ -approximators, and Varopoulos extensions in uniform domains

2.4. The choice of the parameter $\kappa$

Definition 2.31 ( $A_{\infty}$ for elliptic measure).

4. Existence of $\varepsilon$ -approximators: proof of Theorem 1.1

4.4. Construction of $\varepsilon$ -approximators

7. An example in $\mathbb{R}^{3}$