Two-Phase Elliptic Measure

• Two-phase elliptic measure is the study of harmonic and elliptic measures on complementary domains, emphasizing the Radon–Nikodym derivative and transmission conditions on the common boundary.
• It employs monotonicity formulas, such as Almgren’s frequency and the ACF functional, combined with blow-up analysis to establish free-boundary regularity under specific geometric and log density conditions.
• By extending to variable-coefficient operators and VMO regimes, the theory bridges harmonic analysis, geometric measure theory, and multi-operator free-boundary problems.

Two-phase elliptic measure concerns a pair of complementary domains, typically denoted $\Omega^+$ and $\Omega^-$, together with the corresponding interior and exterior harmonic or elliptic measures $\omega^+$ and $\omega^-$. The central datum is the Radon–Nikodym derivative

$$h(Q)=\frac{d\omega^-}{d\omega^+}(Q),\qquad Q\in \partial\Omega,$$

or, more generally, the mutual absolute continuity of $\omega^+$ and $\omega^-$. In the harmonic case this leads to a two-phase free-boundary problem with transmission condition on the common boundary $\Gamma=\partial\Omega$; in the variable-coefficient setting it leads to structural questions for the boundary via blow-ups, tangent measures, and geometric measure theory (Engelstein, 2014, Prats et al., 2019, Goering et al., 4 Sep 2025).

1. Core setting and principal objects

In the basic harmonic-measure formulation, $\Omega^+\subset \mathbb R^n$ is an open 2-sided non-tangentially accessible (NTA) domain, and

$$\Omega^-:=\operatorname{Int}(\mathbb R^n\setminus \overline{\Omega^+}),$$

with $\Omega^-$0 as their common boundary. One writes $\Omega^-$1 and $\Omega^-$2 for the harmonic measures on $\Omega^-$3 seen from $\Omega^-$4 and $\Omega^-$5, assumes mutual absolute continuity $\Omega^-$6, and defines the Poisson-kernel ratio $\Omega^-$7 (Engelstein, 2014). In the lecture-note formulation of Engelstein’s work, the Green’s functions $\Omega^-$8 and $\Omega^-$9 are taken with poles at $\omega^+$0, equivalently on unbounded NTA domains, and the boundary condition is expressed in transmission form as

$\omega^+$1

where $\omega^+$2 denotes the normal derivative pointing into $\omega^+$3 (Engelstein, 2014).

The geometric hypotheses are formulated in terms of NTA connectivity and flatness. For $\omega^+$4, the NTA condition is given by the corkscrew condition and the Harnack-chain condition in the sense of Jerison–Kenig; if $\omega^+$5 is also NTA, one has a two-sided NTA pair $\omega^+$6 (Prats et al., 2019). Reifenberg flatness is quantified by

$\omega^+$7

with $\omega^+$8, and vanishing Reifenberg flatness means

$\omega^+$9

(Prats et al., 2019).

In the elliptic setting, one replaces the Laplacian by divergence-form operators

$\omega^-$0

where $\omega^-$1 are uniformly elliptic, symmetric, matrix-valued coefficients, each 2-quasicontinuous, satisfying

$\omega^-$2

for almost every $\omega^-$3 and all $\omega^-$4 (Goering et al., 4 Sep 2025). For fixed poles $\omega^-$5, the corresponding elliptic measures $\omega^-$6 are supported on $\omega^-$7 and solve the Dirichlet problem via the Riesz representation (Goering et al., 4 Sep 2025).

2. Two-phase free-boundary formulation for harmonic measure

The sharp regularity theory in the harmonic case is organized around the assumption

$\omega^-$8

with $\omega^-$9 and $$h(Q)=\frac{d\omega^-}{d\omega^+}(Q),\qquad Q\in \partial\Omega,$$0 (Engelstein, 2014). The main theorem stated in Engelstein’s exposition is that if $$h(Q)=\frac{d\omega^-}{d\omega^+}(Q),\qquad Q\in \partial\Omega,$$1 is a 2-sided NTA domain, or a Lipschitz domain, and $$h(Q)=\frac{d\omega^-}{d\omega^+}(Q),\qquad Q\in \partial\Omega,$$2, then:

1. if $$h(Q)=\frac{d\omega^-}{d\omega^+}(Q),\qquad Q\in \partial\Omega,$$3, $$h(Q)=\frac{d\omega^-}{d\omega^+}(Q),\qquad Q\in \partial\Omega,$$4 is locally given by the graph of a $$h(Q)=\frac{d\omega^-}{d\omega^+}(Q),\qquad Q\in \partial\Omega,$$5 function;
2. if $$h(Q)=\frac{d\omega^-}{d\omega^+}(Q),\qquad Q\in \partial\Omega,$$6, then under an a priori small Reifenberg-flatness or Lipschitz assumption, the same conclusion holds.

Moreover, if $$h(Q)=\frac{d\omega^-}{d\omega^+}(Q),\qquad Q\in \partial\Omega,$$7 then $$h(Q)=\frac{d\omega^-}{d\omega^+}(Q),\qquad Q\in \partial\Omega,$$8 is $$h(Q)=\frac{d\omega^-}{d\omega^+}(Q),\qquad Q\in \partial\Omega,$$9, and if $\omega^+$0 is real-analytic then $\omega^+$1 is analytic. In the basic case $\omega^+$2 one obtains $\omega^+$3 (Engelstein, 2014).

A central feature of this theory is that there is no a priori non-degeneracy in the free boundary condition. The regularity argument therefore does not begin from a standard quantitative separation of phases; instead, non-degeneracy must be established from monotonicity formulae (Engelstein, 2014). This distinguishes the two-phase harmonic-measure problem from one-phase free-boundary results of Alt–Caffarelli, Jerison, and De Silva–Ferrari–Salsa, which are described in the source exposition as requiring an a priori non-degeneracy in the free-boundary condition (Engelstein, 2014).

This formulation places the Poisson-kernel ratio at the center of the analysis. The boundary regularity of $\omega^+$4 is not prescribed directly; rather, it is recovered from the regularity of $\omega^+$5 and the structure of the harmonic measures themselves. A plausible implication is that two-phase elliptic measure serves simultaneously as boundary data and as a geometric probe of the interface.

3. Monotonicity, blow-ups, and quantitative regularity

The basic harmonic-measure argument uses three nearly-monotone functionals for

$\omega^+$6

The first is Almgren’s frequency,

$\omega^+$7

when $\omega^+$8. For a truly harmonic $\omega^+$9, one has $\omega^-$0; in the two-phase setting $\omega^-$1 is “almost harmonic,” and the source states that

$\omega^-$2

which is integrable at $\omega^-$3, so that $\omega^-$4 as $\omega^-$5 (Engelstein, 2014).

The second is the Alt–Caffarelli–Friedman functional

$\omega^-$6

In the truly harmonic two-phase setting $\omega^-$7 is monotone; in the present setting it is “almost” monotone, and this yields the pointwise lower bound

$\omega^-$8

(Engelstein, 2014). The third is Monneau’s functional

$\omega^-$9

with $\Gamma=\partial\Omega$0 a 1-homogeneous polynomial. By relating $\Gamma=\partial\Omega$1 to $\Gamma=\partial\Omega$2, one deduces that $\Gamma=\partial\Omega$3 exists; together with $\Gamma=\partial\Omega$4, this forces uniqueness of the tangent half-plane $\Gamma=\partial\Omega$5 and rules out degenerate free boundary behaviour (Engelstein, 2014).

These monotonicity arguments yield the radial density

$\Gamma=\partial\Omega$6

which exists and satisfies $\Gamma=\partial\Omega$7 for each $\Gamma=\partial\Omega$8. Consequently, $\Gamma=\partial\Omega$9 can be decomposed into a regular set where the tangent is a plane and a negligible remainder (Engelstein, 2014).

The proof strategy then proceeds through blow-ups and compactness. After obtaining an initial Lipschitz bound for $\Omega^+\subset \mathbb R^n$0 from the ACF functional, one fixes $\Omega^+\subset \mathbb R^n$1 and rescales

$\Omega^+\subset \mathbb R^n$2

extracting a subsequence with $\Omega^+\subset \mathbb R^n$3 in $\Omega^+\subset \mathbb R^n$4, where $\Omega^+\subset \mathbb R^n$5 is a 1-homogeneous harmonic polynomial with $\Omega^+\subset \mathbb R^n$6 a hyperplane. The almost-monotonicity of $\Omega^+\subset \mathbb R^n$7 and Monneau’s functional shows that $\Omega^+\subset \mathbb R^n$8 is unique, so $\Omega^+\subset \mathbb R^n$9 has a unique tangent plane at every $$\Omega^-:=\operatorname{Int}(\mathbb R^n\setminus \overline{\Omega^+}),$$0 (Engelstein, 2014).

The passage from $$\Omega^-:=\operatorname{Int}(\mathbb R^n\setminus \overline{\Omega^+}),$$1 to $$\Omega^-:=\operatorname{Int}(\mathbb R^n\setminus \overline{\Omega^+}),$$2 uses an improvement-of-flatness scheme. Once $$\Omega^-:=\operatorname{Int}(\mathbb R^n\setminus \overline{\Omega^+}),$$3 and the unique tangent plane $$\Omega^-:=\operatorname{Int}(\mathbb R^n\setminus \overline{\Omega^+}),$$4 are known, the boundary lies within $$\Omega^-:=\operatorname{Int}(\mathbb R^n\setminus \overline{\Omega^+}),$$5 of $$\Omega^-:=\operatorname{Int}(\mathbb R^n\setminus \overline{\Omega^+}),$$6 in every sufficiently small ball $$\Omega^-:=\operatorname{Int}(\mathbb R^n\setminus \overline{\Omega^+}),$$7, with $$\Omega^-:=\operatorname{Int}(\mathbb R^n\setminus \overline{\Omega^+}),$$8 uniformly in $$\Omega^-:=\operatorname{Int}(\mathbb R^n\setminus \overline{\Omega^+}),$$9; this implies that $\Omega^-$00 is locally a $\Omega^-$01 graph (Engelstein, 2014). The one-sided and two-sided Harnack lemmata then allow one to push a gap from the boundary of a ball into its center when the solution is trapped between rotated two-plane solutions. In the key iterative lemma, if in $\Omega^-$02 one has

$\Omega^-$03

and $\Omega^-$04 varies by at most $\Omega^-$05, then for any $\Omega^-$06 one finds a new direction $\Omega^-$07, a new height $\Omega^-$08 with $\Omega^-$09, and a smaller error $\Omega^-$10 so that

$\Omega^-$11

Iterating yields $\Omega^-$12 regularity for some $\Omega^-$13 depending only on dimension and $\Omega^-$14, and a second iteration with the sharp Hölder control on $\Omega^-$15 yields exactly $\Omega^-$16 (Engelstein, 2014).

Higher regularity is obtained by a partial hodograph transform in each phase $\Omega^-$17, which converts the free-boundary PDE into a fully elliptic system for two height-functions $\Omega^-$18 over the common boundary. The system is checked to be elliptic with coercive boundary conditions, and a nonlinear Schauder theory of Agmon–Douglis–Nirenberg / Kinderlehrer–Stampacchia promotes $\Omega^-$19 from $\Omega^-$20 to $\Omega^-$21, yielding $\Omega^-$22 regularity of $\Omega^-$23 by iteration (Engelstein, 2014).

4. VMO regularity and geometric equivalences

The VMO regime replaces pointwise Hölder control of $\Omega^-$24 by vanishing mean oscillation relative to harmonic measure. For a doubling measure $\Omega^-$25, the source defines $\Omega^-$26 by

$\Omega^-$27

and $\Omega^-$28 precisely when

$\Omega^-$29

(Prats et al., 2019). In the two-phase setting one asks that

$\Omega^-$30

Prats and Tolsa prove a two-phase analogue of the Kenig–Toro one-phase characterization. If $\Omega^-$31 is a bounded NTA domain, $\Omega^-$32 is also NTA, and $\Omega^-$33 is $\Omega^-$34-Reifenberg flat for $\Omega^-$35 sufficiently small, then the following are equivalent:

• $\Omega^-$36 and $\Omega^-$37 are mutually absolutely continuous and

$\Omega^-$38

• $\Omega^-$39 is vanishing Reifenberg flat, the inner unit normal $\Omega^-$40 to $\Omega^-$41 exists $\Omega^-$42-a.e. and belongs to $\Omega^-$43, and both Radon–Nikodym derivatives satisfy a reverse Hölder inequality of exponent $\Omega^-$44;
• $\Omega^-$45 is vanishing Reifenberg flat, $\Omega^-$46, and either $\Omega^-$47 or $\Omega^-$48;
• $\Omega^-$49 is vanishing Reifenberg flat, $\Omega^-$50 and $\Omega^-$51 have joint big pieces of chord-arc subdomains, and

$\Omega^-$52

where $\Omega^-$53 is the unit normal to the best-approximating plane over $\Omega^-$54 (Prats et al., 2019).

The proof uses several distinct ingredients. In the implication $\Omega^-$55, the mutual absolute continuity of $\Omega^-$56 and the capacity-density condition yield rectifiability and tangent $\Omega^-$57-planes $\Omega^-$58-a.e.; jump formulas for the Riesz transform relate boundary integrals of the normal $\Omega^-$59 to non-tangential limits of singular integrals against $\Omega^-$60; a Christ–David stopping-time cube decomposition isolates good regions where the blow-up of $\Omega^-$61 is small; and Carleson-measure or square-function estimates, combined with the John–Nirenberg/VMO condition, produce the vanishing oscillation estimate for $\Omega^-$62 (Prats et al., 2019).

The implication $\Omega^-$63 proceeds through approximating chord-arc subdomains obtained by bumping $\Omega^-$64 up and down on a Whitney-type covering of the bad set; on each approximator, the one-phase Kenig–Toro theorem gives that the logarithm of its Poisson kernel belongs to VMO of surface measure, and the estimate is transferred back by the maximum principle for harmonic measures (Prats et al., 2019). This equivalence makes explicit that VMO regularity in the two-phase logarithm encodes both vanishing flatness and quantitative geometric approximation by chord-arc pieces.

5. Multi-operator two-phase elliptic measure

The multi-operator theory replaces a single operator by a pair of divergence-form operators with distinct coefficients on the two sides of the boundary. Let

$\Omega^-$65

so that on $\Omega^-$66 the two measures are mutually absolutely continuous (Goering et al., 4 Sep 2025). The main structural theorem states that if $\Omega^-$67 are complementary NTA domains and $\Omega^-$68 are 2-quasicontinuous uniformly elliptic symmetric matrices, then there exists $\Omega^-$69 of full $\Omega^-$70-measure and a decomposition

$\Omega^-$71

such that:

1. for $\Omega^-$72,

$\Omega^-$73

one has $\Omega^-$74 on $\Omega^-$75, and $\Omega^-$76;

1. on $\Omega^-$77, one has $\Omega^-$78;
2. $\Omega^-$79 (Goering et al., 4 Sep 2025).

Here $\Omega^-$80 consists of the points where $\Omega^-$81, and $\Omega^-$82 those where the Radon–Nikodym derivative vanishes (Goering et al., 4 Sep 2025). The theorem therefore separates the boundary into a full-measure set with flat tangent measures, a singular part where the two elliptic measures are mutually singular, and a null remainder.

The blow-up construction is explicit. For $\Omega^-$83 and radii $\Omega^-$84, one rescales

$\Omega^-$85

where $\Omega^-$86 are the Green’s functions with fixed poles $\Omega^-$87 (Goering et al., 4 Sep 2025). After passing to a subsequence, one obtains locally Hausdorff convergence of the rescaled domains, local uniform convergence of $\Omega^-$88, and weak-* convergence of $\Omega^-$89 to a common limit measure $\Omega^-$90 (Goering et al., 4 Sep 2025).

In the limit, $\Omega^-$91 are unbounded complementary NTA domains with common boundary $\Omega^-$92, $\Omega^-$93 is their harmonic measure with pole at $\Omega^-$94, and

$\Omega^-$95

is a global $\Omega^-$96 solution of

$\Omega^-$97

with

$\Omega^-$98

Equivalently, after the linear change of variable $\Omega^-$99, one reduces to

$\omega^+$00

where

$\omega^+$01

with $\omega^+$02 on $\omega^+$03 and a matching normal-flux condition that becomes a two-plane-solution jump inequality in the viscosity formulation (Goering et al., 4 Sep 2025).

This reduction shows that multi-operator two-phase elliptic measure is controlled by a corresponding multi-operator free-boundary problem. The geometric content of the structural theorem is therefore mediated by the rigidity of the free-boundary blow-ups.

6. Tangent measures, rigidity, and unresolved questions

The tangent-measure framework follows Preiss. A nonzero Radon measure $\omega^+$04 is a tangent measure $\omega^+$05 if for some $\omega^+$06 and $\omega^+$07,

$\omega^+$08

weak-*, where $\omega^+$09 (Goering et al., 4 Sep 2025). A $\omega^+$10-cone is a family of measures closed under multiplication by positive constants and dilations, and Preiss’s connectedness lemma gives a dichotomy once the ambient cone has the requisite property $\omega^+$11 (Goering et al., 4 Sep 2025).

In the multi-operator application, the relevant cone is the cone of blow-up elliptic measures

$\omega^+$12

and the flat cone $\omega^+$13 consists of $\omega^+$14-plane measures (Goering et al., 4 Sep 2025). The source verifies property $\omega^+$15 by a Liouville-type rigidity statement: if a blow-up measure is sufficiently close to flat at all large scales, then a global Liouville theorem for the two-phase free boundary, obtained through Caffarelli’s “flat$\omega^+$16Lipschitz$\omega^+$17two-plane” program, forces the measure itself to be flat (Goering et al., 4 Sep 2025). Since touching corkscrew balls in $\omega^+$18 produce a two-plane tangent measure, the alternative $\omega^+$19 is excluded, and one concludes that all tangent measures are flat for $\omega^+$20-almost every $\omega^+$21 (Goering et al., 4 Sep 2025).

Several limitations are explicit. In the single-operator harmonic case, mutual absolute continuity $\omega^+$22 forces $\omega^+$23 and $\omega^+$24-rectifiability. In the multi-operator setting, one still has

$\omega^+$25

but the reverse inequality and full rectifiability remain open because an Alt–Caffarelli–Friedman monotonicity formula “in the case of coefficient jumps fails to characterize flatness” (Goering et al., 4 Sep 2025). Accordingly, the paper gives only a partial answer to Bishop’s question and leaves the multi-operator Oksendal conjecture open (Goering et al., 4 Sep 2025).

Open directions also appear in the harmonic and VMO theories. The lecture-note exposition of Engelstein’s work lists possible extensions to second-order elliptic operators with variable coefficients and corresponding elliptic measure, to two-phase problems for $\omega^+$26-harmonic measure or nonlinear PDEs, and to finer geometric analysis of the singular set where the tangent polynomial has degree at least $\omega^+$27 (Engelstein, 2014). Prats and Tolsa state that the same equivalences should hold for unbounded NTA domains with poles at $\omega^+$28, and they identify singular sets and stratification, as well as higher-order analogues such as biharmonic or polyharmonic measures and fully nonlinear elliptic operators, as active directions (Prats et al., 2019).

A recurrent misconception is that mutual absolute continuity alone automatically yields the full smooth or rectifiable picture. The cited results are more specific. In the Hölder regime, one obtains sharp $\omega^+$29 and then $\omega^+$30 regularity under the stated assumptions on $\omega^+$31 and the ambient geometry (Engelstein, 2014). In the VMO regime, one obtains an equivalence with vanishing Reifenberg flatness, VMO normal, and chord-arc approximation under the stated small-flatness hypotheses (Prats et al., 2019). In the multi-operator setting, mutual absolute continuity yields flat tangent measures almost everywhere and the upper bound $\omega^+$32, but not the missing lower bound on dimension (Goering et al., 4 Sep 2025). Together these results describe two-phase elliptic measure as a boundary invariant whose analytic regularity, free-boundary structure, and geometric consequences depend sensitively on the operator, the oscillation class of the logarithmic density, and the quantitative flatness available a priori.