1db8542 354ac36 eb6554c c011cca
================================================================================
SHARP INTEGRAL FILLINGS IN CAT(0) SPACES
the isoperimetric inequality in spaces of nonpositive curvature,
at the best possible price
a complete account, readable by a high school student
================================================================================

This file is plain text on purpose. There is no LaTeX, no browser and no
lecturer needed. Every symbol is defined the first time it appears, and
section 10 is a glossary in case you meet a symbol again later and forget
it. Section 11 lists every result of this file in order. Section 12 is one
page that summarises the entire proof.

Results are numbered in the order they appear. [L3] is the third result and
it is a lemma; [P7] is the seventh result and it is a proposition; [T1] is a
theorem. The three main results of the file are [T1] (the sharp filling
inequality), [T2] (its geometric consequence, the Cartan-Hadamard domain
inequality) and [T3] (the induction statement that drives the proof).
Results taken from other people's papers are marked [lit] and are stated
exactly, with their hypotheses, in section 4.

A warning about size. The statement of the main theorem fits in five lines.
Its proof is long because it must work in a space whose only information is
a distance, and where surfaces have integer weights and may be wildly
singular. Nothing in the proof is optional: each section removes one
obstacle. Section 0.4 lists the obstacles so that you always know which one
is being attacked.


================================================================================
0.  THE SHORT VERSION
================================================================================

0.1  The question.

A soap film wants to be as small as possible. If you dip a wire loop into
soap water, the film that forms is (to very good approximation) the surface
of least area among all surfaces that have the wire as boundary. The oldest
question in the subject is: how small can that surface be?

Start in the flat plane. A closed curve of length L bounds a disk of area
L^2/(4 pi), and no surface with that boundary has less area, because a disk
has the least area among all sets with given perimeter (perimeter is the
two dimensional word for length of the boundary; area is the volume of the
region). For a curve that is not a circle the disk of the same length is
strictly bigger, so L^2/(4 pi) is an upper bound for the least filling area
of any closed curve of length L, and it is attained only by the round one.

The same question in every dimension. A closed n-dimensional surface of
area A bounds an (n+1)-dimensional region of volume V. Euler's constant
pi is the first of a family: write

    omega_j  = volume of the unit ball of R^j
    s_n      = area of the unit sphere S^n inside R^(n+1)
             = (n+1) * omega_(n+1)

so omega_1 = 2, omega_2 = pi, omega_3 = 4 pi / 3, and s_0 = 2, s_1 = 2 pi,
s_2 = 4 pi, s_3 = 2 pi^2. The Euclidean isoperimetric inequality says

    V  <=  A^((n+1)/n) / ( (n+1)^((n+1)/n) * omega_(n+1)^(1/n) ) .       (0.1)

The ball is the equality case: for the sphere of radius R we have
A = s_n R^n and V = omega_(n+1) R^(n+1), and (0.1) holds with equality. So
(0.1) cannot be improved, and the "sharp constant"

    c_n  =  1 / ( (n+1) * s_n^(1/n) )                                    (0.2)

is the best possible coefficient in the inequality V <= c_n A^((n+1)/n).

In numbers: c_2 = 1/(6 sqrt(pi)) = 0.09403..., and
c_3 = (128 pi^2)^(-1/3) = 0.09248...

Now let the room be curved. Our question becomes: if the ambient space is
not flat R^(n+1) but some curved space X, and a closed n-dimensional surface
inside X has area A, how much volume must a filling of it have? Equivalently
and more usefully: if the filling has n+1 dimensional volume V, how much
area must its boundary have?

0.2  The answer.

The main theorem of this file. The precise statement is [T1] in section 3;
here it is.

    Let X be a proper CAT(0) space and let n >= 2. Let T be an "integral
    n-cycle" of compact support in X, with mass A. Then T bounds an
    integral current S of compact support with

        boundary of S = T,        mass of S <= c_n A^((n+1)/n),          (0.3)

    where c_n is the Euclidean constant (0.2). There is no restriction on
    the dimension of X, on the topology or singularities of T, or on the
    integer multiplicities that T carries.

Two words need a promise behind them.

"CAT(0) space" means a metric space in which distance is measured along
geodesics and triangles are no fatter than their Euclidean comparison
triangles. The plane, all hyperbolic spaces, all trees, Hilbert space and
every convex subset of any of these are CAT(0). Nonpositive curvature in
the classical Riemannian sense implies CAT(0). Section 1.2 gives the exact
definition and section 1.3 gives the examples.

"Integral current" means an oriented n-dimensional surface with an integer
weight at each point, allowing pieces that are only Lipschitz (so, corners
and kinks are allowed) and allowing weights bigger than one (so a surface
can be covered twice). Its mass is the weighted area. Its boundary is the
oriented (n-1)-dimensional surface that is "the edge", and a cycle is a
current whose boundary is zero. Section 2 builds this up from nothing.

So the theorem says: in any CAT(0) space, however singular, and for cycles
of any dimension n >= 2, the Euclidean price for filling is still the true
price, even though there is no linear structure, no differentiable charts,
no normal bundle, and no smoothness anywhere.

0.3  Why the answer is not obvious.

Three separate reasons, and each one is a section of the proof.

(a) Curvature can make filling worse, and there is a cheap example showing
    the exponent can change completely. Take X to be an infinite tree with
    one branch point where three unit segments meet (a tripod). Take n = 1:
    a 1-cycle is a formal difference of two points, and its only fillings
    are the paths joining them, so filling two points at distance 2 costs 2.
    But the Euclidean bound for n = 1 would cost 2^2/(4 pi) = 0.32. False by
    a factor of six. So the theorem is genuinely about n >= 2; for n = 1 the
    truth is a linear law, not a (n+1)/n law. This is the honest content of
    the hypothesis n >= 2, and no proof can drop it.

    Positive curvature also breaks the inequality: on a round sphere of
    radius R, the equator has length 2 pi R and the cheapest thing it bounds
    is a hemisphere of area 2 pi R^2, while the Euclidean price would be
    (2 pi R)^2/(4 pi) = pi R^2, three times smaller. So nonpositive curvature
    is not decoration either.

(b) There are no charts. In a CAT(0) space we cannot integrate, cannot take
    derivatives, and cannot even write down a surface in coordinates. The
    theory of integral currents in metric spaces (Ambrosio-Kirchheim) gives
    the right definitions, but then all the classical machinery - the area
    formula, the coarea formula, the Euler-Lagrange equation - has to be
    re-derived from the distance alone. Sections 2 and 7 do that.

(c) The extremal configuration is never smooth. The sharp constant is
    attained by a round ball, so the natural proof has to look at what
    happens near a round ball. But the objects in the proof are extremizers
    of a variational problem in a singular space; they can have fractional
    exponents and can vanish on open sets. Sections 7.7 to 7.12 are the
    careful handling of exactly those two degeneracies.

0.4  The shape of the proof.

Here is the whole strategy in seven paragraphs. Every word here is made
precise later; read it once now and once at the end.

Step 1 (reduce to a compact room). A compactly supported cycle lives inside
some closed ball, and closed balls in a CAT(0) space are convex and, if the
space is proper, compact. So it is enough to prove the theorem inside a
compact CAT(0) space Y. [L4].

Step 2 (assume the theorem fails, get a champion). Inside Y, the least
filling mass of a cycle B is written V(B). Suppose the sharp bound fails for
some cycle. Then for a constant c slightly bigger than c_n, the functional

    B  |-->  V(B) - c * (mass of B)^((n+1)/n)                             (0.4)

is positive somewhere. A compactness theorem for currents gives a cycle A
that maximizes (0.4). Then A has a first variation property [L9]: for every
cycle B,

    V(A - B)  >=  c * ( m^((n+1)/n) - (mass of B)^((n+1)/n) ),            (0.5)

where m = mass of A. Step 3 (a rescaling of all distances) makes

    m < s_n ,        c * ((n+1)/n) * m^(1/n) = 1/n .                      (0.6)

So the mass of the champion is slightly less than the area of the unit
sphere in its own dimension. That numerical coincidence is the engine: the
proof will show that any cycle satisfying (0.5) and (0.6) must have
m >= s_n, and contradiction.

Step 3 (radial variation). Pull every point of A towards a fixed centre y at
a speed that depends on the point. Each point x is moved to the point at
fraction lambda(x) on the segment [y,x], where lambda = 1 - t g(x) and g is
a nonnegative Lipschitz function. Sliding A for time h sweeps out a
filling D_h of (A - B_h), where B_h is the moved copy of A. Comparing the
mass lost by A with the mass swept, and using the champion inequality (0.5),
gives the radial inequality [P1]:

    integral ( n g + r * (grad r . grad g) ) dmu
        <=  n * integral ( r g sqrt(1 - |grad r|^2) ) dmu                 (0.7)

for every centre y and every nonnegative Lipschitz g, where r(x) = d(y,x),
mu is the mass measure of A, and grad, dots and |.| are the chart
derivatives of section 2.5.

Step 4 (base case, n = 2). Feed (0.7) the cutoff g = 1/max(r, eps)^2 and
let eps go to 0. A Euclidean lower density bound at almost every point of
A turns that into m >= 4 pi = s_2, contradicting (0.6). So the theorem
holds for 2-cycles. [P3].

Step 5 (induction, the small set estimate). Now n > 2 and we assume the
theorem in dimension n-1. Cut the champion A by a level set of a Lipschitz
function. Each small piece has an (n-1)-dimensional boundary, and the
induction hypothesis fills that boundary cheaply. Feeding the result back
into the champion inequality (0.5) gives a perimeter bound: every small
piece of A has boundary mass at least (1 - eta) n omega_n^(1/n) * (its
mass)^((n-1)/n), the Euclidean sphere bound. [P4].

Step 6 (from sets to functions to functions). Two classical moves in metric
form. First, coarea plus rearrangement [P5, L12, L13] convert the small set
perimeter bound into an "almost Euclidean" Sobolev inequality [P6]:

    (S_* - eps) (integral |u|^p dmu)^(2/p)
       <=  4 beta integral |grad u|^2 dmu + C_eps integral u^2 dmu ,       (0.8)

where S_* = n s_n^(2/n), beta = 1/(n-2) and p = 2n/(n-2). Second, that
inequality makes a certain variational problem on A have a minimizer v with
[P8]

    W = integral v^p dmu  in  (0, s_n),      4 beta E(v) + n |v|_2^2 = n W.

The strict inequality W < s_n is inherited from (0.6) and is the whole
point: v is a slightly subcritical Sobolev bubble living on the current A.

Step 7 (chords force contradiction). Normalize dnu = v^p dmu / W, a
probability measure on Y, and consider for a centre y the squared weighted
chord

    s_y(x) = v(y)^beta d(y,x) v(x)^beta ,     L_y = integral s_y^2 dnu .

Testing the minimizer equation and the radial inequality (0.7) with the
functions built from the Euclidean bubble profile (1 + a s^2)^(-(n-2)/2)
gives the transform bound [P9]

    integral (1 + a s_y^2)^(-n) dnu  <=  (1 + 2 a L_y)^(-n/2) ,           (0.9)

which forces L_y <= 2 (W/s_n)^(2/n) < 2 for almost every centre y [L15].
On the other hand, test the minimizer equation once more with v^(1+2 beta)
and use Cauchy-Schwarz + the CAT(0) variance inequality to get the opposite
averaged bound [L16]:

    integral L_y dnu  >=  2 .                                            (0.10)

Since (W/s_n)^(2/n) < 1, (0.9) says the average is strictly less than 2 and
(0.10) says it is at least 2. Contradiction. So no champion exists, the
sharp inequality holds in dimension n, and induction gives the theorem in
every dimension. [T1].

Step 8 (the classical payoff). In a Cartan-Hadamard manifold there are no
compactly supported top-dimensional cycles other than zero [L17], so the
filling of the boundary of a domain must be the domain itself. Applying
[T1] to that filling gives the classical Euclidean domain inequality for
every bounded smooth domain in every dimension d >= 3. [T2].

0.5  Where this sits in the literature.

The Euclidean case is Almgren (1986): the sharp inequality (0.1) for
integral currents in R^m, all dimensions and codimensions. Doing the same in
a metric space needs currents in metric spaces, which is Ambrosio and
Kirchheim (2000). Wenger (2005) proved an isoperimetric inequality with the
Euclidean exponent in spaces with convex geodesic bicombings, with a
constant depending on the space and on the dimension, not sharp.

For the applications, the classical statement is: in a Cartan-Hadamard
manifold (complete, simply connected, sectional curvature <= 0), domains of
dimension d satisfy the Euclidean volume-perimeter inequality, namely

    area of the boundary  >=  d * omega_d^(1/d) * volume^((d-1)/d) .      (0.11)

That is "the Cartan-Hadamard conjecture". Dimension 2 is classical (Weil
1926; Beckenbach-Rado 1933). Dimension 3 is Kleiner (1992), dimension 4 is
Croke (1984); dimensions 3 to 9 were done by Chen-Ghomi-Wang (2026), and
various conditional and small-volume versions exist. In this file's
theorem, all dimensions simultaneously are consequences of one sharp
current inequality. The main result here is the filling inequality [T1]
itself: sharp coefficient, all cycle dimensions n >= 2, all proper CAT(0)
spaces, arbitrary multiplicities, no dimension bound on the ambient space.

Related but different: Lang-Stadler-Urech (2026) and Peteranderl (2026)
prove that in CAT(0) spaces of finite asymptotic Nagata dimension, fillings
can be much cheaper than Euclidean at very large mass (linear, or almost
linear). Those theorems are about the growth exponent at large mass and
their constants depend on the space. [T1] is about the exact coefficient at
every mass. There is no conflict: different regimes, different claims.

0.6  How to read this file.

Three reading routes.

  (i) For the statement and the idea: sections 0, 1, 3, then 12. About
      twenty minutes.

  (ii) For the actual proof: add section 2 (what a current is and how mass
      is computed), then 5, 6, 7 in order. Section 4 is a list of borrowed
      results; read it as a shopping list and go back to it when section 7
      says "by [lit 3]".

  (iii) For a full understanding including the technical layers: read
      everything, and then section 9, which explains honestly which steps are
      deep, which are fiddly, and why each fiddly step is unavoidable.

Exercises, when they appear, are short and are marked "Exercise:". They are
not part of the proof. If an exercise is hard, skip it; nothing later needs
it.

================================================================================
1.  THE SETTING
================================================================================

1.1  Distance, geodesics, triangles.

A metric space is a set X with a function d(x,y) >= 0 that is zero exactly
on the diagonal, is symmetric, and satisfies the triangle inequality
d(x,z) <= d(x,y) + d(y,z). That is all the structure a CAT(0) space uses.

A geodesic segment from x to y is a map gamma: [0,1] -> X with
gamma(0) = x, gamma(1) = y, and d(gamma(s), gamma(t)) = |s-t| d(x,y) for
all s,t in [0,1]. In words: it is a path that covers the distance exactly,
at constant speed. The point at fraction a is written [x,y]_a, so
[x,y]_0 = x and [x,y]_1 = y; this notation is used constantly below. X is a
geodesic space if every pair of points is joined by a geodesic segment.
Note that a geodesic segment need not be unique, but in a CAT(0) space it
always is (see 1.2); we then speak of "the" segment.

A geodesic triangle with vertices x_1,x_2,x_3 is a choice of three segments
joining them pairwise. Its Euclidean comparison triangle is a triangle in
the plane with the same three side lengths; it exists because three
positive numbers satisfying the three triangle inequalities are the side
lengths of a Euclidean triangle. Points on the comparison triangle
correspond to points at the same fractions along the sides.

1.2  CAT(0).

Definition. X is CAT(0) if X is a geodesic space and, for every geodesic
triangle with vertices x_1,x_2,x_3 and every pair of points p,q on two sides
of it, the corresponding points p',q' of the Euclidean comparison triangle
satisfy

    d(p,q)  <=  |p' - q'| .                                              (1.1)

Inequality (1.1) says: triangles in X are no fatter than Euclidean
triangles. It is the exact metric version of "sectional curvature <= 0", and
on a Riemannian manifold it is equivalent to it (Cartan-Alexandrov-Toponogov).
The "(0)" in CAT(0) is the curvature bound.

Three consequences are used all the time. All three follow from (1.1) in a
few lines.

  (a) Geodesics are unique, and the space is "Busemann convex": if x -> x_t
      and y -> y_t are geodesics with the same parameter, then
      t |--> d(x_t, y_t) is convex. In particular the function
      x |--> d(y,x) is convex along every geodesic.

  (b) Balls are convex: for any centre o and radius R, the closed ball is a
      geodesic space in its own right, with the restricted distance.

  (c) Shortening a segment by a fixed fraction is a contraction: for fixed
      t in [0,1], the map x |--> [o,x]_t is t-Lipschitz. In symbols,

          d([o,x]_t, [o,x']_t)  <=  t d(x,x') .                          (1.2)

      Both (a) and (c) are immediate from (1.1) applied to the triangle
      o,x,x' and its comparison triangle, where the corresponding map is
      the honest Euclidean contraction towards o. We shall also use the
      "equal fraction" comparison, valid for two different fractions
      lambda, lambda':

          d([o,x]_lambda, [o,x']_lambda')^2
             <= (lambda r - lambda' r')^2
                + lambda lambda' ( d(x,x')^2 - (r - r')^2 )              (1.3)

      where r = d(o,x) and r' = d(o,x'). This is exactly the squared
      distance of the corresponding two points in the Euclidean comparison
      triangle, so it is a consequence of (1.1). It is the key estimate of
      section 6, where lambda and lambda' differ because they depend on x
      and on x' through a function g.

Finally, the four point inequality. For any four points a,b,c,d of a CAT(0)
space,

    d(a,c)^2 + d(b,d)^2
       <= d(a,b)^2 + d(b,c)^2 + d(c,d)^2 + d(d,a)^2 .                    (1.4)

Proof (short, and worth reading once, because it is used to prove that
chart metrics are Euclidean in section 2.5). Let m be the midpoint of
[a,c]. Midpoint comparison in the triangles a,b,c and a,d,c says

    2 d(b,m)^2 + d(a,c)^2/2 <= d(a,b)^2 + d(b,c)^2 ,
    2 d(d,m)^2 + d(a,c)^2/2 <= d(a,d)^2 + d(c,d)^2 .

Add them and use d(b,d)^2 <= 2 d(b,m)^2 + 2 d(d,m)^2, which is the triangle
inequality squared. The two d(a,c)^2/2 terms cancel against the common
factor 2 d(...)^2 and (1.4) follows.

Inequality (1.4) is a parallelogram law with an inequality sign. In
Euclidean space it is an equality for the four points of a parallelogram.
It says that CAT(0) spaces are, in a weak integral sense, no more curved
than flat ones; the analytic payoff is that a "metric derivative" of a
Lipschitz map out of R^k, which is a priori only a seminorm, satisfies the
parallelogram law and hence comes from an honest inner product. That is the
content of [L5] in section 2.5, and it is the reason the whole proof can
compute determinants and Jacobians at all.

1.3  What CAT(0) admits, and what it forbids.

The examples. R^m with its Euclidean distance, all hyperbolic spaces H^m,
any Hilbert space, all trees with the path metric, all CAT(0) cubical
complexes, all Euclidean buildings of noncompact type, all Gromov-Hausdorff
limits of CAT(0) spaces, and every convex subset of a CAT(0) space. The
class is closed under products with the Euclidean factor and under passing
to convex closed subsets, which is what step 1 of the proof does.

What it forbids. The space must be geodesic and its geodesics must be
unique, so, for instance, R^2 with the l^1 distance (where the segment from
(0,0) to (1,1) is not unique) is not CAT(0), even though it is flat in a
weaker sense. CAT(0) also forbids positive curvature, so a round sphere is
not CAT(0); a ball of radius larger than a quarter of the circumference in a
sphere is not CAT(0) either.

Why nobody objects to the singular examples. Trees, buildings and limits are
exactly where the classical differential-geometric methods fail, and exactly
where this theorem still gives a sharp answer. That is the point of proving
it with currents rather than with the calculus of variations on manifolds.

1.4  Proper spaces and the compact reduction.

X is proper if every closed bounded subset of X is compact. R^m, trees with
locally finite branching, and complete Riemannian manifolds are proper;
infinite dimensional Hilbert space is not.

Properness is used for one thing only: to make closed balls compact so that
compactness theorems for currents apply. Since a compactly supported current
lives in a bounded set, and closed balls in a CAT(0) space are convex, we
can always move the whole problem into a compact CAT(0) space: [L4] in
section 5 states and proves that reduction. From section 5 onward the
ambient space is a compact CAT(0) space named Y, and the letter X is only
used again in the final assembly.

1.5  Balls, spheres, and the constants.

Constants. Let omega_j be the volume of the unit ball of R^j (so
omega_1 = 2, omega_2 = pi, omega_3 = 4 pi/3, omega_4 = pi^2/2), and let s_n
be the area of the unit sphere S^n in R^(n+1). The relation is
s_n = (n+1) omega_(n+1): differentiate the volume of the ball of radius R,
which is omega_(n+1) R^(n+1), with respect to R and set R = 1. Thus
s_1 = 2 pi, s_2 = 4 pi, s_3 = 2 pi^2.

Three constants appear in the sharp statements:

    s_n           area of the unit sphere in dimension n
    c_n = 1/((n+1) s_n^(1/n))                    the sharp filling constant
    S_* = n s_n^(2/n)                            the sharp Sobolev constant

and two more in the induction, beta = 1/(n-2) and p = 2n/(n-2), which are
the classical Sobolev exponents in dimension n (so p > 2 exactly when
n > 2, which is why n = 2 is done separately).

For orientation: c_2 = 1/(6 sqrt(pi)) = 0.09403..., c_3 = (128 pi^2)^(-1/3)
= 0.09248..., and c_n decreases slowly to zero, roughly like a constant
times n^(-3/2).

Exercise. Check that c_n = omega_(n+1)^(1/(n+1)) / (n+1) * something is NOT
the formula; the right formula is c_n = 1/((n+1) s_n^(1/n)). Then check with
n = 2 that 1/(3 * sqrt(4 pi)) = 1/(6 sqrt(pi)).

1.6  Worked check: the ball is optimal.

Take X = R^(n+1), let T be the boundary of the ball of radius R, oriented.
Then

    mass of T = s_n R^n =: A ,
    mass of the ball = omega_(n+1) R^(n+1) =: V ,

and by (0.2),

    c_n A^((n+1)/n) = s_n R^(n+1) / (n+1) = omega_(n+1) R^(n+1) = V .

So the bound in [T1] is attained for every R. By [L17] (top dimensional
constancy in a noncompact manifold, section 8) the ball is the ONLY
compactly supported filling of T in R^(n+1). Hence the constant c_n cannot
be improved even in flat space, which settles optimality; and it also shows
the theorem is exactly the statement "every CAT(0) space is at least as
expensive as flat space", which is what nonpositive curvature is supposed to
mean.


================================================================================
2.  CURRENTS: THE RIGHT NOTION OF SURFACE
================================================================================

2.1  Why ordinary surfaces are not enough.

A soap film is not a smooth surface. It meets itself along curves, three
sheets meet at 120 degrees along a curve, and a minimizing sequence of
smooth surfaces can collapse to something with kinks and crossings. So if we
define "surface" as "smooth embedded manifold", two things go wrong at once.

First, no compactness. The set of smooth surfaces with a given boundary and
bounded area is not closed: take a wavy sheet and let the waves get finer
and finer with the same amplitude. The areas stay bounded but the surfaces
have no limit in any reasonable sense. Since the whole proof needs to look
at a surface of least area, and least area is defined by an infimum over a
sequence, this is fatal.

Second, no cancellation. If you run a surface one way and then the same
surface the other way, the union is a legitimate "filling" of the boundary
zero, and it should count as nothing. Smooth embedded surfaces cannot
express that; oriented surfaces can, if you count multiplicity with signs.

The theory of currents fixes both problems at once. A k-current is an
object for which two things make sense: "integrate a bounded Lipschitz
weight against the Jacobian of k coordinate functions" and "take the
oriented boundary". It contains every oriented Lipschitz k-surface with
integer multiplicity, it is closed under weak limits, and the limit of a
minimizing sequence exists and is again a current. Everything hard in this
file is a consequence of working in that completed world.

2.2  The definition, and what each piece means.

Standing notation for this file: X is a complete metric space (later a
compact CAT(0) space).
  Lip(X)    = real valued functions on X with finite Lipschitz constant.
  Lip_b(X)  = the bounded ones.
For a Lipschitz function pi, its Lipschitz constant is written Lip(pi).

Definition (Ambrosio-Kirchheim). Let k >= 0 be an integer. A k-current on X
is a map

    C : Lip_b(X) x Lip(X) x ... x Lip(X)  ->  R     (k copies of Lip(X))

which is linear in each of the k+1 slots separately, antisymmetric under
permuting the last k slots, and continuous under the following weak
convergence: if b_j -> b pointwise with a common bound on |b_j| and on
Lip(b_j), and if pi^j_a -> pi_a pointwise with a common bound on Lip(pi^j_a),
then C(b_j, pi^j_1, ..., pi^j_k) -> C(b, pi_1, ..., pi_k).

Reading the definition. Take the model case: an oriented surface S given by
a bi-Lipschitz parametrization phi: E -> X, E a Borel subset of R^k, with an
integer weight theta. Then the corresponding current is

    C(b, pi_1, ..., pi_k) = int_E theta(z) b(phi(z))
                             det D( (pi_1,...,pi_k) o phi )(z) dz .       (2.1)

That is: integrate the weight b against the Jacobian determinant of the k
coordinate functions along the surface. If you choose pi_1,...,pi_k to be
honest coordinates and b = 1, this is exactly the k-dimensional volume of
the piece. The point of the definition is that (2.1) makes sense for every
tuple of Lipschitz tests, and knowing all those numbers determines C.

Mass. For a finite Borel measure eta on X, say that eta controls C if

    |C(b, pi_1, ..., pi_k)|  <=  int |b| deta    for all b, pi_a with
                                                       Lip(pi_a) <= 1.   (2.2)

The mass of C is the smallest such eta, written ||C||; it exists (the
supremum of the left sides over admissible tuples is additive over disjoint
Borel sets, so it defines a measure). The total mass is M(C) = ||C||(X), and
for the model current (2.1) built from a surface,

    ||C|| = (phi)_# ( |theta| * J dz ),   J = sqrt(det G) the area density,

so M(C) is exactly "area counted with the absolute value of the weight". The
formula (2.2) with the sup over all tuples is the honest definition of mass
even for currents that are not surfaces; it says that mass measures all the
ways the current can be tested.

Why the definition is written with k coordinate functions. The determinant
in (2.1) is the volume distortion in the k coordinate directions. In a
metric space we have no coordinates, but we do have real valued Lipschitz
functions, and k of them together measure a k-dimensional distortion. All
of the smooth theory is recovered by taking the pi_a to be coordinates.

Boundary. The boundary of a k-current C is the (k-1)-current

    (dC)(b, pi_1, ..., pi_{k-1})  =  C(1, b, pi_1, ..., pi_{k-1}) .      (2.3)

For a curve (k = 1) this reads (dC)(b) = C(1,b) = integral of b over the
boundary of the curve, i.e. b(end) - b(start). So (2.3) is the algebraic
version of "the edge of the surface". One checks d(dC) = 0 whenever the
objects involved have finite mass; that algebraic identity is why slicing,
cutting and gluing work in this theory without any smoothness.

Cycles, fillings, filling volume. C is a cycle if dC = 0 (zero current).
C is normal if M(C) + M(dC) < infinity, and integral if in addition it has
the representation

    C = sum_i (phi_i)_# ( theta_i [[E_i]] ),                            (2.4)

with phi_i bi-Lipschitz from a bounded Borel set E_i onto its image, the
images pairwise disjoint Borel sets, theta_i integer valued and integrable
against the area density, and the sum of masses finite. Write I_k(X) for
the integral k-currents. For B in I_k(X) with dB = 0, put

    V(B) = inf { M(S) : S in I_(k+1)(X), dS = B } .                      (2.5)

V(B) is the least filling mass of B, and its finiteness is a theorem, not a
definition (see [L5] for an explicit filling, and [lit6] for the Euclidean
exponent bound). [T1] is precisely the statement that
V(B) <= c_k M(B)^((k+1)/k)
for cycles in a proper CAT(0) space of any dimension, with k >= 2.

The whole point of integral currents, rather than arbitrary normal ones, is
the integrality of the weight: it will give a lower density bound of section 2.7
which is what forces the constant s_n in the base case, and which is the
only place the integers theta are used quantitatively. Running the proof
with real weights would give a strictly weaker inequality.

2.3  Charts.

The representation (2.4) is the metric replacement for "a surface is a union
of coordinate patches". Two facts make it usable.

  (i) By subdividing and using the disjointified images, we may always take
      the chart domains E_i bounded and the images Z_i = phi_i(E_i)
      pairwise disjoint. The mass is then a sum of absolutely continuous
      measures on the charts, and the charts are "almost flat at almost
      every point" (that is the next subsection).

  (ii) Each phi_i is bi-Lipschitz, so its inverse is defined on Z_i and has
      a Lipschitz constant; by McShane's extension theorem [lit 11] every
      Lipschitz function on a subset of a metric space extends to the whole
      space with the same constant. This is what lets us test a current on
      "coordinates of the chart" even though the chart is only a subset of
      X: we extend and test on the extension; the value depends only on the
      values on the chart, because the current is supported there.

Fact (i) is algebraic bookkeeping. Fact (ii) is used twice with real force:
once to compute mass exactly ([L3]) and once to close the gradient in
section 7.7. Both times, the point is the same: a chart which is bi-Lipschitz
can be "unrolled" by Lipschitz functions that the current can be evaluated
on, with no loss in the constants that matter.

2.4  The metric derivative, and why chart metrics are Euclidean.

Classical area formula: for an injective C^1 map phi: R^k -> R^m, the
Jacobian factor is sqrt(det (D phi)^T (D phi)). In a metric space there is no
derivative, but there is a replacement.

Definition (metric differential). Let phi: E -> X be Lipschitz, E a Borel
subset of R^k, and let z be a density point of E. A seminorm N on R^k is the
metric differential of phi at z if the two-point formula

    d( phi(z+w), phi(z+w') )  =  N(w - w') + o(|w| + |w'|)               (2.6)

holds as the two domain points tend to z inside E.

Kirchheim's theorem [lit 2] says that at almost every density point of E a
metric differential exists (it is a seminorm, not necessarily a norm, and
not necessarily unique, but almost every point has one and any two agree).
In a general metric space (2.6) is the best structure available; in
particular the "norm" can be a weird polyhedral thing. In CAT(0) space it
cannot: it must be an honest Euclidean norm. That is our first real result.

[L1] (Euclidean chart metrics). Let X be CAT(0), phi: E -> X bi-Lipschitz,
and let z be a density point of E where the metric differential N exists.
Then
   (a) N is a norm, that is N(w) > 0 for w != 0;
   (b) N comes from an inner product: there is a symmetric positive
       definite k x k matrix G with N(w)^2 = w^T G w for all w in R^k;
   (c) with this G, the centered form of (2.6) holds:

           d(phi(z+w), phi(z))  =  |w|_G + o(|w|) ,   |w|_G^2 = w^T G w;
   (d) the entries of G are measurable functions of z on E, and J = sqrt(det G)
       is a measurable area density with
           Lip(phi^(-1))^(-k) <= J(z) <= Lip(phi)^k                    (2.7)
       at almost every z in E.

Proof. (a) Nondegeneracy. phi is bi-Lipschitz: there is c > 0 with
d(phi(z),phi(z')) >= c |z - z'|. Fix a direction and a scale; a
density-one point has the property that any fixed finite configuration in
R^k is approximated, to relative error o(1), by the rescaled set (E - z)/t
for small t. (If a point of the configuration stayed a positive distance
away from (E-z)/t along a sequence of scales, a ball of radius comparable to
t, sitting at distance O(t) from z, would be missing from E; that contradicts
density one.) Applying this to two points z, z+w in E and using the
bi-Lipschitz lower bound gives N(w) >= c|w|. Together with the upper bound
N(w) <= Lip(phi)|w| (which is just (2.6) and the Lipschitz bound), this makes
N a norm.

(b) Parallelogram law. Take four domain points 0, u, u+v, v in R^k and
approximate them by rescaled points of E - z, as in (a). Apply the four point
inequality (1.4) of section 1.2 to the four images, divide by t^2, and use
(2.6). The result is

    N(u+v)^2 + N(u-v)^2  <=  2 N(u)^2 + 2 N(v)^2 .

The reverse inequality is the same inequality applied to the pair
(u+v)/2, (u-v)/2, together with homogeneity of N. So

    N(u+v)^2 + N(u-v)^2  =  2 N(u)^2 + 2 N(v)^2 ,

a parallelogram law. Polarize: define the bilinear form

    <u,v> = ( N(u+v)^2 - N(u-v)^2 ) / 4 .

The parallelogram law is exactly the statement that <.,.> is additive in each
slot, and homogeneity gives bilinearity. It is positive definite because N is
a norm. Hence N(w)^2 = w^T G w with G the matrix of <.,.>.

(c) Put w' = 0 in (2.6) and use that N(w) = |w|_G.

(d) Each entry of G is a finite combination of squared metric differentials
in the fixed coordinate directions; more precisely, G_ab = (1/2)( N(e_a+e_b)^2
- N(e_a)^2 - N(e_b)^2 ), so measurability follows from measurability of
z |--> N_z(w) for fixed rational w, which follows from (2.6) and continuity
of the distance. The bounds (2.7) are consequences of the bi-Lipschitz
constants: eigenvalues of G lie between Lip(phi^(-1))^(-2) and Lip(phi)^2.
End of proof of [L1].

Two remarks, both used later.

  (i) Equation (2.6) with N Euclidean means the chart looks like a small
      copy of Euclidean space, up to relative error o(1), for at almost every
      point. So all the classical computations in Euclidean space - the
      determinant bounds, the coarea formula, the divergence theorem applied
      to a homotopy - transfer to a CAT(0) chart with no dimensional loss.
      This is the entire reason the proof is possible.

  (ii) The matrix G depends on the chart and on the point, and is only
      measurable. There is no uniform bound on Lip(phi) across the (countable)
      chart family. That is why all the estimates below are pointwise in the
      chart and then summed/integrated, and never uniform.

2.5  Mass from a quadratic majorant: Hadamard's inequality does the work.

The next two results are the ones that let us compute mass inside a CAT(0)
space with no loss. The first is an upper bound for mass, given a quadratic
upper bound for the distance distortion. Intuitively: if a chart distorts
small displacements by no more than a quadratic form P, then the area
density of that chart is no more than sqrt(det P), because a linear map that
sends the unit ball of P into the unit ball has determinant at most one.

Setting for [L2]. Let C be a finite-mass d-current with a representation by
countably many determinant integrals:

    C(b, pi_1, ..., pi_d) = sum_i int_{F_i} vartheta_i(z) b(psi_i(z))
                       det D( (pi_1,...,pi_d) o psi_i )(z) dz ,            (2.8)

where psi_i: F_i -> Y are Lipschitz, F_i in R^d Borel, and vartheta_i are
real signed weights (not yet assumed integer). Assume that at almost every
point carrying these weights there is a measurable nonnegative quadratic
form P_i(z) on R^d with

    d(psi_i(z+w), psi_i(z))^2  <=  P_i(z)[w,w] + o(|w|^2) ,  z+w in F_i . (2.9)

[L2] (quadratic majorant for mass). In the setting above, put

    eta  =  sum_i (psi_i)_# ( |vartheta_i| sqrt(det P_i) dz ) .            (2.10)

If eta(Y) < infinity, then eta controls C in the sense of (2.2), and
therefore M(C) <= eta(Y). No normality of the individual chart currents, and
no assumption that the psi_i are injective, is required.

Proof. Fix d scalar 1-Lipschitz functions pi_1,...,pi_d on Y. Almost
everywhere in each chart, the scalar compositions pi_a o psi_i are
differentiable in the sense of scalar Lipschitz extensions (that is,
pi_a o psi_i(z+w) - pi_a o psi_i(z) = ell_a(w) + o(|w|) with ell_a a
bounded linear functional). Write ell_a for the derivative row of pi_a o
psi_i at z, and let E be the d x d matrix with rows ell_1,...,ell_d.

Step 1: each row is bounded by P_i. Take any direction w in R^d. Choosing
domain points z + t w_j with w_j -> w and using (2.9) together with the
1-Lipschitz property of pi_a gives

    |ell_a(w)|  <=  sqrt( P_i(z)[w,w] ) .                                 (2.11)

Indeed the left side is the derivative of pi_a o psi_i in the direction w,
and the increment of pi_a is bounded by the distance increment of psi_i,
which is bounded by sqrt(P_i[w,w]) t + o(t).

Step 2: Hadamard's inequality in the P-geometry. Suppose first that P_i is
positive definite and put Q = P_i(z)^{1/2}. By (2.11), the row ell_a is the
composition of Q with a row vector of Euclidean length at most 1. Hence
E = S Q with each row of S of length at most 1, and

    |det E|  =  |det S| sqrt(det P_i)  <=  sqrt(det P_i) ,                (2.12)

because a matrix whose rows all have length at most 1 is the matrix of a
linear map that sends the unit cube into the unit ball, so it cannot increase
volume. (That last sentence is Hadamard's inequality, and it is the whole
content of this step.) If P_i is singular, then by (2.11) every row vanishes
on the kernel of P_i, so E has rank less than d as soon as the kernel is
nonzero; then both det E and det P_i are zero and (2.12) still holds.

Step 3: conclude. In the representation (2.8), the determinant is det E, so
each chart integral is bounded in absolute value by the chart integral of
|vartheta_i| sqrt(det P_i) against the test function |b|. Summing over the
charts, the finiteness of eta makes everything summable (the integrals of the
absolute values converge to the integral of the sum by monotone convergence)
and gives (2.2) with the measure eta. Minimality of the mass measure
finishes the proof.

One point of bookkeeping, because it is the kind of thing that is easy to
get wrong. The set where a chart is not differentiable depends on the chosen
tuple (pi_1,...,pi_d), so the exceptional set moves with the tuple. That is
harmless: mass is defined as the smallest measure that controls EVERY tuple,
and the definition of eta does not refer to a tuple at all. So one fixed
measure eta controls all tuples separately, which is exactly what (2.2)
demands. End of proof of [L2].

2.6  Exact chart mass.

Now the reverse inequality, which says that the chart formula for mass is an
identity and not merely an upper bound. This is what lets us convert any
integral over the current into an ordinary integral over a chart domain, and
it is used in every section after this one. In particular it identifies the
null sets: a set is null for the mass measure exactly when its inverse image
in the charts has Lebesgue measure zero.

[L3] (exact chart mass). Let C in I_d(Y) be given by the pairwise disjoint
charts of (2.4), with integer weights theta_i, metric matrices G_i and area
densities J_i = sqrt(det G_i) from [L1]. Then

    ||C||  =  sum_i (phi_i)_# ( |theta_i| J_i dz ) ,                      (2.13)

as measures on Y. Consequently, for a Borel set B in Y,

    sum_i int_{phi_i^(-1)(B)} |theta_i| J_i dz  =  ||C||(B) ,             (2.14)

and the sum on the left converges absolutely.

Proof. Let eta denote the right side of (2.13), allowing eta(Y) = infinity
for the moment.

First, eta <= ||C||. Fix 0 < eps < 1. We claim that each chart domain, up to
a Lebesgue-null set, is a disjoint countable union of Borel pieces E on which
the chart is (1 +- eps)-almost isometric to a fixed linear model, meaning:
there is a positive definite matrix Q (from a fixed countable dense family)
with

    (1-eps)|z - z'|_Q  <=  d(phi_i(z), phi_i(z'))  <=  (1+eps)|z - z'|_Q
                                                     for all z,z' in E .  (2.15)

Construction of the pieces. At a density point z, the metric differential
G_i(z) is close to some Q of the dense family (relative error at most eps/4,
say). The centered estimate [L1](c) then provides a radius 1/l, depending on
z, on which (2.15) holds with that Q against every other point of the domain.
Group points by the matrix Q and by the integer l, then subdivide into pieces
of diameter < 1/l, and disjointify. Borel measurability: for a fixed Q and
radius, the inequalities in (2.15) against all points of the domain can be
checked against a countable dense subset of the domain, using continuity of
the distance in the second variable, so each piece is Borel.

Now fix one such piece E of one chart. At almost every point of E,
differentiating (2.15) in the domain direction gives J_i(z) <= (1+eps)^k
sqrt(det Q), and the inverse map Q^{1/2} phi_i^(-1) is (1-eps)^(-1)-Lipschitz
on phi_i(E) by the lower inequality in (2.15). Extend each of its k coordinate
functions to all of Y by McShane (same Lipschitz constant) and call the
resulting tuple pi.

Use these tests in (2.1) with the first slot b chosen to be a bounded Borel
function, supported in phi_i(D) for a Borel D in E, equal there to
sgn(theta_i) o phi_i^(-1). The determinant in the chart is exactly
det(Q^{1/2}) = sqrt(det Q) at almost every point of E, because the extended
coordinates agree with Q^{1/2} z there. All other charts give zero, because
their images are disjoint from phi_i(D) and b is supported in phi_i(D).
Hence

    int_D |theta_i| sqrt(det Q) dz
       = C(b, pi_1, ..., pi_k)
       <= (1-eps)^(-k) ||C||(phi_i(D)) .                                  (2.16)

The last inequality is the mass bound (2.2) applied to these k functions,
each of Lipschitz constant at most (1-eps)^(-1). Combining (2.16) with
J_i <= (1+eps)^k sqrt(det Q) gives

    int_D |theta_i| J_i dz  <=  ((1+eps)/(1-eps))^k ||C||(phi_i(D)) .     (2.17)

Now let B be Borel in Y, take D = E intersected with phi_i^(-1)(B), sum over
all pieces and all charts: the images are disjoint so the right side adds up
to a constant times ||C||(B), and Lebesgue-null sets contribute nothing to
eta. Therefore eta(B) <= ((1+eps)/(1-eps))^k ||C||(B), first with eta(Y) =
infinity allowed and then, taking B = Y, as a finite bound. Let eps go to
zero to get eta <= ||C||.

Second, ||C|| <= eta. This is now immediate from [L2]: the quadratic
majorant (2.9) holds with P_i = G_i by the centered form of [L1](c), the
measure (2.10) is exactly eta, and eta is finite by the first part. So
||C|| <= eta. End of proof of [L3].

Two consequences we will use without comment.

  (i) The area density J_i is positive almost everywhere (by the bi-Lipschitz
      bound (2.7)), so the inverse images of null sets are null and the
      charts really do carry the mass measure faithfully.

  (ii) Integrating a Lipschitz function u against ||C|| is the same as
      integrating u o phi_i times |theta_i| J_i over the chart domains. All
      the Sobolev theory of section 7 is built this way, with the "chart
      gradient" Du computed in the metric G_i of [L1].

2.7  Euclidean lower density.

Here the integers theta enter for the first and only time in a way that
affects a constant.

[L4] (Euclidean lower density). Let C in I_k(Y). Then

    liminf_{rho -> 0} rho^(-k) ||C||( B(y,rho) )  >=  omega_k            (2.18)

for ||C||-almost every y. (Balls are open; omega_k is the volume of the unit
ball of R^k.)

Proof. By (2.13) it suffices to consider points y = phi_i(z) where z is a
density-one point of E_i, the metric differential exists at z, z is a
Lebesgue point of the density rho_i = |theta_i| J_i (extended by zero off
E_i), and |theta_i(z)| >= 1. These conditions hold at almost every point
carrying chart mass, because the weight is an integer and nonzero there.

Fix eps > 0. For all sufficiently small rho, every domain point

    z + w   with   |w|_{G_i(z)} < rho/(1+eps)

lies in E_i and satisfies d(y, phi_i(z+w)) < rho, by the centered estimate
[L1](c). The set of such w is an ellipsoid of Lebesgue volume

    omega_k rho^k / ( (1+eps)^k J_i(z) ) .

Lebesgue differentiation at z on dilates of this fixed ellipsoid says that
the chart contribution to ||C||(B(y,rho)) is at least

    ( |theta_i(z)| J_i(z) + o(1) ) * omega_k rho^k / ( (1+eps)^k J_i(z) ) .

Dividing by rho^k, letting rho -> 0, then eps -> 0, and using
|theta_i(z)| >= 1, gives (2.18). Other charts contribute nonnegative mass,
so they can only help. End of proof of [L4].

Read (2.18) as: an integral current cannot be thinner than k-dimensional
Lebesgue measure with coefficient omega_k. That is a purely Euclidean
statement about the model space, but it applies at almost every point of a
current in a wildly curved space, because at that point the chart is
metrically Euclidean to first order by [L1].

2.8  Products, cones, and the two coarse bounds on V.

Two constructions are needed before the variational argument can start: an
explicit filling of any cycle (a cone), and the compactness/continuity
properties of the filling volume V. We record them together.

Products. Let C in I_k(Y) and h > 0. The whole-current product of the
interval [0,h] with C, written [0,h] x C, is the (k+1)-current on
[0,h] x Y with the product metric. It is integral [lit 3], and its boundary
is

    d([0,h] x C) = [h] x C - [0] x C - [0,h] x (dC) .                    (2.19)

The signs in (2.19) are the ones that make d(d(.)) = 0; the third term is the
only place the boundary of C appears, and for a cycle it vanishes. The action
of the product on test functions is a sum of k+1 terms, one for each
coordinate slot, each of them an integral over t in [0,h] of an evaluation of
C with the time derivative in the first slot [lit 3]. We use this formula
exactly twice, in sections 6 and 7.6; both times it is applied to a homotopy
F: [0,h] x Y -> Y, so what is really used is the composition rule for
pushforwards.

[L5] (cone filling, and the linear bound). Let Y be a compact CAT(0) space,
B in I_k(Y) a cycle, and D = diam(spt B) (zero if B = 0). Then V(B) is
finite, and

    V(B)  <=  (k+1) D M(B) .                                             (2.20)

If B = 0 then V(0) = 0. Also, for every cycle B in Y,

    V(B)  <=  K_k M(B)^((k+1)/k)   and   V(B) <= L_k M(B) ,              (2.21)

with constants K_k, L_k depending only on k and on the diameter of Y, and
every cycle has a mass minimizing integral filling.

Proof. The zero case is trivial. Otherwise choose o in spt B, and let
H: [0,1] x Y -> Y be H(t,x) = [o,x]_t, the point at fraction t on the
segment from o to x; it is well defined because geodesics in a CAT(0) space
are unique. Two estimates, both from section 1.2:

    d(H(t,x), H(t,x')) <= t d(x,x') ,        (contraction, (1.2))
    d(H(t,x), H(s,x))  = |t-s| d(o,x) ,      (constant speed)

and d(o,x) <= D on spt B. Hence H is jointly Lipschitz on [0,1] x Y, with
spatial constant 1 and time constant D.

Now let S = H_#([0,1] x B). By (2.19) its boundary is
(H_1)_# B - (H_0)_# B - H_#([0,1] x dB). The first map is the identity, the
second is the constant map to o (whose pushforward of a positive dimensional
current is zero), and dB = 0. So dS = B: the cone over B is a filling.

Mass of the cone. Write the action of S as a sum of k+1 terms, as in the
product formula, with test functions psi_a(t,x) = pi_a(H(t,x)). In each of
those k+1 terms, one factor is a time derivative of a 1-Lipschitz scalar
composed with H, hence bounded by D; the remaining k spatial derivatives are
at most 1-Lipschitz. Applying the mass bound (2.2) for B to each term and
adding the k+1 of them yields

    M(S)  <=  (k+1) D M(B) ,                                             (2.22)

which proves (2.20).

The second bound in (2.21) is (2.20) with L_k = (k+1) diam(Y). The first
bound in (2.21) is Wenger's isoperimetric theorem [lit 6]: a compact CAT(0)
space has a convex geodesic bicombing (its unique geodesics), and then
V(B) <= K_k M(B)^((k+1)/k) with K_k depending on k and the space. We never
use the value of K_k, only its existence.

Existence of minimizers. Take a minimizing sequence S_j for V(B): bounded
masses, common boundary B. The compactness theorem for integral currents in a
compact metric space [lit 4] gives a weakly convergent subsequence with
integral limit S, and dS = B because boundary commutes with weak convergence.
Mass is lower semicontinuous under weak convergence [lit 4], so M(S) <=
liminf M(S_j) = V(B), and S is a minimizer. End of proof of [L5].

Finally, the continuity of filling volume, which is what makes the maximizing
argument of section 5 work. Suppose B_j, B are integral k-cycles in Y with
sup_j M(B_j) < infinity and B_j -> B weakly. Then [lit 7]

    V(B_j - B)  ->  0   and   V(B_j) -> V(B) .                           (2.23)

The second limit follows from the first and subadditivity,
|V(B_j) - V(B)| <= V(B_j - B), which is just the observation that a filling
of the difference fills the difference. The first limit is Wenger's flat
convergence theorem [lit 7], and it needs the local cone inequalities: the
cone over a piece of a cycle has mass comparable to its mass, which is exactly
(2.20) applied inside small balls. That is why we bothered to prove (2.20)
by hand.


================================================================================
3.  THE THEOREM, THE COROLLARY, AND THE INDUCTION
================================================================================

3.1  The main theorem.

[T1] (sharp integral filling in CAT(0) spaces). Let X be a proper CAT(0)
space, let n >= 2, and let T in I_n(X) have compact support and dT = 0.
Then there exists S in I_(n+1)(X) with compact support such that

    dS  =  T ,                                                          (3.1)

    M(S)  <=  c_n M(T)^((n+1)/n)
          =  M(T)^((n+1)/n) / ( (n+1)^((n+1)/n) omega_(n+1)^(1/n) ) ,   (3.2)

where c_n = 1/((n+1) s_n^(1/n)) as in (0.2). There is no hypothesis on the
dimension of X, and no hypothesis on the multiplicity function of T: it may
take any integer values, with either sign. The coefficient c_n is optimal,
by the ball computation of section 1.6.

The rest of this file is the proof. The proof structure, stated as three
separate results so that nothing is hidden:

    [T1] is implied by [P2] (the case n = 2), by [T3] (the induction step)
    and by [L6] (the reduction from proper X to compact Y).

Concretely: [L6] moves a compactly supported cycle into a closed ball;
[P2] handles 2-cycles in a compact CAT(0) space; [T3] handles n-cycles
given that the theorem is known for (n-1)-cycles. Induction over n gives
[T1] in compact spaces; [L6] gives it back in proper spaces.

3.2  What the theorem says, in one sentence, and what it does not say.

One sentence: in a space of nonpositive curvature, filling a closed
n-surface of area A never costs more than the flat-space optimum
c_n A^((n+1)/n), and the optimum is achieved by round spheres in flat space.

It does not say that the minimizer is smooth, unique, or a piece of a
sphere. It does not classify the equality case in a curved space. In a
general CAT(0) space the least filling of a cycle need not have any
regularity beyond integrality, and in fact the round balls of flat space are
the only known equality cases in this generality. The statement is an
inequality for the infimum (2.5), which is attained by [L5].

It also does not give a construction. The proof is by contradiction, and
produces a filling only through the compactness theorem: the infimum is
attained, and the argument shows the infimum is at most c_n A^((n+1)/n).

A different question, which this file does not answer, is Gromov's: given a
map from a fixed polyhedral complex into X, how much extra volume does one
need to extend it to a map from the cone over the complex? That question
fixes the topology and the parametrization of the filling. [T1] allows the
filling to be any integral current with the right boundary, which is why the
two questions are different. See section 0.5 for the literature on that.

3.3  The Cartan-Hadamard consequence.

A Cartan-Hadamard manifold is a complete, simply connected Riemannian
manifold whose sectional curvature is at most zero. By the Cartan-Hadamard
theorem plus Hopf-Rinow such a manifold is a proper CAT(0) space and is
diffeomorphic to R^d. So [T1] applies to it, and the top-dimensional case of
[T1] becomes the classical geometric inequality.

[T2] (Cartan-Hadamard domain inequality). Let M be a Cartan-Hadamard
manifold of dimension d >= 3, and let Omega be a bounded domain with smooth
boundary. Then

    area(dOmega)  >=  d * omega_d^(1/d) * vol(Omega)^((d-1)/d) .        (3.3)

The same inequality holds for relatively compact sets of finite perimeter,
with the perimeter of the set in place of area(dOmega).

The proof is three lines once the lemma "there are no compactly supported
top-dimensional cycles in a connected noncompact manifold" is available: the
domain is itself a filling of its boundary, and in a noncompact connected
oriented manifold the only compactly supported d-cycle inside I_d is zero, so
the minimizer of [T1] must BE the domain. Hence mass of the minimizer =
vol(Omega) is at most the right side of (3.2) with n = d-1, and a little
algebra produces (3.3). That lemma, and the algebra, are section 8; the
formula check is: with n = d-1 we have s_n = d omega_d, so

    c_(d-1)^(-(d-1)/d) = ( d s_(d-1)^(1/(d-1)) )^((d-1)/d)
                       = d^((d-1)/d) (d omega_d)^(1/d)
                       = d omega_d^(1/d) ,

so vol <= c_(d-1) area^(d/(d-1)) is the same as (3.3).

For d = 2 (surfaces) the same inequality is classical, going back to Weil
and Beckenbach-Rado, and it is not a consequence of [T1] because [T1] starts
at n = 2, that is d = 3. Section 8.3 explains the two-dimensional argument:
cut the domain into disks, apply the classical disk inequality (which follows
from the Gauss-Bonnet formula plus the general isoperimetric inequality), and
sum.

Dimensions 3 and 4 of (3.3) were known (Kleiner 1992, Croke 1984) and
dimensions 3 to 9 were known (Chen-Ghomi-Wang 2026) before this theorem. What
is new here is not a single dimension: it is that one sharp current
inequality in a general CAT(0) space implies all dimensions at once, without
curvature hypotheses and without a smoothness or comparison argument.

3.4  The induction step, stated as a theorem.

[T3] (induction step). Let n >= 3, let Y be a compact CAT(0) space, and
assume:

    (H_(n-1))   for every compact CAT(0) space Z and every integral
                (n-1)-cycle B in Z, there is an integral n-current R in Z
                with dR = B and M(R) <= c_(n-1) M(B)^(n/(n-1)) .

Then for every integral n-cycle T in Y there is an integral (n+1)-current S
in Y with dS = T and M(S) <= c_n M(T)^((n+1)/n).

Three remarks about [T3].

  (a) It is an implication between statements in consecutive dimensions, and
      it never mentions X. Together with the base case [P2] (n = 2) and the
      reduction [L6], it proves [T1] by induction on n. The induction is
      over the dimension of the CYCLE, not over the dimension of the space;
      the space is arbitrary in every step.

  (b) The hypothesis (H_(n-1)) is used in exactly one place: to fill the
      boundaries of small pieces of the champion in the small-set
      inequality [P3], section 7.1. Everywhere else the proof of [T3] is
      self-contained.

  (c) The proof of [T3] is a contradiction argument, and it is long because
      it has to pass through an analytic inequality ([P4]) and a
      variational problem ([P7]). The pieces, in order, are:
        section 5   the champion A and its normalization (uses no induction)
        section 6   the radial inequality [P1] (uses no induction)
        section 7.1 the small-set inequality [P3] (uses H_(n-1))
        section 7.2 coarea and countable mass tests ([L10], [L11])
        section 7.3 rearrangement and the radial Sobolev inequality
                     ([L12], [L13])
        section 7.4 the almost sharp Sobolev inequality [P4]
        section 7.5 the closed gradient and compactness ([P5], [P6])
        section 7.6 the minimizer [P7] and its powers ([L14], [L15])
        section 7.7 the chord transform [P8] and its two consequences
        section 7.8 the contradiction, which closes the induction

3.5  Why the theorem starts at n = 2, and why n = 2 is separate.

The hypothesis n >= 2 in [T1] is not a convenience. Section 0.3 gave the
example: in a tree, a 1-cycle is a pair of points and the cheapest filling is
the segment between them, whose length is the distance, not the square of the
distance over 4 pi. So the (n+1)/n exponent is false for n = 1 in a CAT(0)
space, and no argument can rescue it. Fortunately, 2-cycles cannot see trees:
a CAT(0) space of dimension 1 carries no nonzero 2-current at all, so the
theorem is vacuous there and the first interesting dimension is 2.

Why dimension 2 is proved separately rather than by the general argument.
The general argument of section 7 uses the Sobolev exponents

    beta = 1/(n-2) ,       p = 2n/(n-2) ,      S_* = n s_n^(2/n) ,

and at n = 2 they blow up: 1/(n-2) is infinite, p is infinite. There is no
Sobolev inequality in the critical dimension n = 2; the correct tool there is
not analysis but a direct computation with a 1/r^2 cutoff. That computation
is short, exact, and is [P2] in section 6.6. It also has a useful side effect:
it shows exactly where the integer weights and the constant s_2 = 4 pi come
from, which makes the general argument easier to believe.


================================================================================
4.  THE TOOLBOX: WHAT IS BORROWED, STATED EXACTLY
================================================================================

Everything in this section is someone else's theorem. Each item is stated in
the form in which this file uses it, with the hypotheses it needs. Nothing
here is proved; if an item is used in a delicate spot, the spot is flagged.

[lit 1] Ambrosio-Kirchheim, currents in metric spaces (2000).
  Content used: (a) the definition of k-current, mass measure and boundary
  of section 2.2, with mass having the minimality property (2.2);
  (b) the rectifiable representation (2.4) of an integral current by
  countably many bi-Lipschitz charts with integer weights;
  (c) pushforward by a Lipschitz map f: X -> Z, defined by
  (f_#C)(b, pi_1,...,pi_k) = C(b o f, pi_1 o f,...,pi_k o f), which sends
  integral currents to integral currents, commutes with boundary
  (d(f_#C) = f_#(dC)), and satisfies M(f_#C) <= Lip(f)^k M(C);
  (d) restriction to a Borel set E, defined by inserting the coefficient
  1_E in the first slot, with the identity ||C|_E|| = ||C|||_E and the bound
  |C(b, pi)| <= int |b| d||C|| for every admissible tuple;
  (e) for a locally Lipschitz (or merely Lipschitz) u and an integral
  current C, the level sets: for almost every real t the restriction
  C|{u > t} is an integral current and its boundary is an integral
  (k-1)-cycle. Also: the differential of a Lipschitz scalar function
  vanishes almost everywhere on each of its level sets, in the chart sense
  of section 2.5.
  Needed in: every section from 2 on. In particular (d) is the inequality
  that makes [L2] and [L3] work, and (e) is what makes slicing legal in
  section 7.2.

[lit 2] Kirchheim, rectifiable metric spaces (1994), in the two point form
  of Duda (2007). Existence, at almost every density point of a Borel set
  E in R^k, of a seminorm N (the metric differential) with
      d(phi(z+w), phi(z+w')) = N(w - w') + o(|w| + |w'|)
  for a Lipschitz phi. In this file it is used only through [L1], which
  upgrades N to a Euclidean norm when the target is CAT(0).
  Needed in: section 2.4.

[lit 3] Wenger (2007), product currents. For C in I_k and h > 0 there is an
  integral current [0,h] x C on the product of the interval with the space,
  with
      d([0,h] x C) = [h] x C - [0] x C - [0,h] x (dC)
  and with the action formula (sum of k+1 terms, time derivative in the
  first slot). The product is a "whole-current" product in the sense that no
  boundary is given to the individual pieces; that is what makes it usable
  in a metric space.
  Needed in: section 2.8 (statement), section 6.2 and section 7.7 (used on
  the homotopy F(t,x) = [y,x]_(lambda(t,x))).

[lit 4] Ambrosio-Kirchheim compactness and closure. If C_j in I_k(Y) with
  sup_j ( M(C_j) + M(dC_j) ) < infinity and Y compact, then a subsequence
  converges weakly to an integral current C; dC_j converges weakly to dC;
  and M(C) <= liminf M(C_j). ("Weak" means: evaluation against every
  admissible tuple converges.)
  Needed in: section 2.8 (existence of minimizers), section 5 (existence of
  the champion), section 7.5 (compactness of the Sobolev embedding).

[lit 5] Ambrosio-Kirchheim slicing. For C in I_k and Lipschitz u, the slices
  <C, u, t> are integral (k-1)-currents for almost every t, the
  codimension-one coarea formula holds, and the fundamental theorem
      d(C|{u > t}) = (dC)|{u > t} + <C, u, t>
  holds for almost every t (with the orientation conventions of [lit 1]).
  Needed in: section 7.2, where it is used in integrated form: the integral
  over an interval of levels of the boundary evaluations of the slices
  equals an evaluation of the original current on a function of u.

[lit 6] Wenger (2005), isoperimetric inequality of Euclidean type. In a
  complete metric space which is geodesic and has a convex geodesic
  bicombing (in particular any complete CAT(0) space), every integral
  k-cycle B has, for every k >= 1,
      V(B) <= K_k M(B)^((k+1)/k)
  with K_k depending only on k and on the local geometry constants.
  Needed in: section 2.8, as the first bound of (2.21), which makes the
  champion's functional bounded above.

[lit 7] Wenger (2007), flat convergence and filling volume continuity. If
  B_j, B are integral k-cycles in a compact metric space with bounded masses
  and B_j -> B weakly, and the space has the local cone inequalities, then
  V(B_j - B) -> 0, and consequently |V(B_j) - V(B)| -> 0.
  Needed in: section 5.2, to show the maximizing functional is upper
  semicontinuous along the minimizing/compactness sequence.

[lit 8] Brezis-Lieb (1983). If u_j -> u almost everywhere with
  sup_j ||u_j||_p < infinity, then
      ||u_j||_p^p - ||u_j - u||_p^p - ||u||_p^p  ->  0 .
  The version used here is the elementary consequence, with a proof, in
  section 7.6: for every eps there is C_eps with
      | |a+b|^p - |a|^p |  <=  eps |a|^p + C_eps |b|^p ,
  which together with dominated convergence gives the same conclusion.
  Needed in: section 7.6 (no mass is lost in the limit).

[lit 9] Aubin (1976) and Talenti (1976). The sharp Sobolev inequality in
  R^n, n >= 3:
      S_* ||f||_p^2 <= 4 beta ||grad f||_2^2 ,  S_* = n s_n^(2/n),
      beta = 1/(n-2) ,  p = 2n/(n-2) ,
  for smooth compactly supported f, with equality for the bubbles
  (a^2 + |x|^2)^(-(n-2)/2). In this file we need only the RADIAL case, and
  we prove it, by mass transport, as [L13]. The literature reference is
  used for context and to reassure the reader that the constant is the
  classical one.
  Needed in: section 7.3 (as orientation; the actual theorem is [L13]).

[lit 10] Cordero-Erausquin, Nazaret, Villani (2004). The mass transport
  method for sharp Sobolev inequalities: monotone rearrangement transports
  the level sets of f to the level sets of a radial profile, and the
  determinant comparison plus an integration by parts gives the sharp
  constant. We run exactly this argument for radial functions, in full, in
  [L13].
  Needed in: section 7.3 (method; details local).

[lit 11] McShane extension, in the quantitative form of Rieffel (2006). A
  real valued L-Lipschitz function on a subset A of a metric space X extends
  to an L-Lipschitz function on X; for a tuple of scalar functions, extend
  them one at a time. (Rieffel's paper identifies the best constant
  available for such extensions, which is what the "no loss" claims in
  [L3] and section 7.5 refer to.)
  Needed in: section 2.3, section 2.6, section 7.5.

[lit 12] Sturm (2003), variance inequality in CAT(0) spaces. For a finite
  measure sigma on a compact CAT(0) space and
      I(y) = int d(y,x)^2 dsigma(x) ,
  the function I is continuous, has a minimizer o (a barycentre), and
      I(y) >= I(o) + M d(o,y)^2   for every y,  M = sigma(X).
  The proof is the segment computation reproduced in section 7.8.
  Needed in: section 7.8, to double the chord contribution.

[lit 13] Almgren (1986), optimal isoperimetric inequalities. The sharp
  Euclidean inequality (0.1) for integral currents in R^m, in all
  dimensions and codimensions, with the ball as the equality case. This file
  does not use it as a black box in the proof; it is the theorem that
  [T1] generalizes, and the reason the constant c_n is the right target.
  Needed in: section 0.5 (context), section 8.4 (sharpness).

[lit 14] Simon (2014), lecture notes on geometric measure theory. Recipe for
  compactness of a family of finite measures from a uniform bound on the
  distributional derivatives: translate, average against a mollifier, use
  Ascoli-Arzela on the smoothed family, then let the mollifier shrink. This
  is the classical BV compactness argument; we run it in section 7.5 for the
  projected measures of the weighted current, in the exact form needed.
  Needed in: section 7.5.

[lit 15] Brezis-Kato (1979) and Brezis (1986). Truncated-power testing.
  Multiplying a weak solution by a power of itself times a cutoff, dropping
  the nonnegative term, and absorbing the "critical" potential with a large
  level set, gives L^q membership for all q. The absorption device is
  reproduced in section 7.6 as [L14].
  Needed in: section 7.6.

Two warnings about how to use this list.

First, the items are not interchangeable with their smooth analogues. For
instance, [lit 5] says that the slices are integral for almost every level,
not for every level, and that "the" boundary identity holds only almost
everywhere in the level. Every place where this file passes from a current
to its slices, it either avoids discussing the exceptional levels or sums
over all levels. The reader should watch for that: it is the most common
place to make a silent error in this subject.

Second, the constants. K_k of [lit 6], Lip(phi_i) of the charts, and the
constants C_eps of section 7.4 are all allowed to depend on the space, the
chart or eps, and NONE of them is ever compared with the sharp constants
c_n, S_*. The whole point of the proof is that the sharp constants survive,
so any step that would trade a sharp constant for a loose one is a red flag.
There are exactly two such steps, both in the base case and both harmless:
the bound 2b - 2b^2 <= 1/2 in [P2], and the absorption constant in [L14].


================================================================================
5.  REDUCTION, AND THE CHAMPION OF A FAILED INEQUALITY
================================================================================

This section does two things. First it removes the word "proper" from the
problem, so that from now on the ambient space is a single compact CAT(0)
space Y. Second it turns a failure of the inequality into an object A with a
very rigid property, the "champion inequality" (5.4), and it normalizes A so
that its mass is slightly below s_n. That normalization is the only place
where the sharp constant c_n is used numerically.

5.1  Every cycle lives in a compact convex ball.

[L6] (compact reduction). Fix n >= 2. Suppose that for every compact CAT(0)
space Y and every integral n-cycle T in Y with dT = 0 there is an integral
(n+1)-current S in Y with dS = T and M(S) <= c_n M(T)^((n+1)/n). Then the
same statement holds in every proper CAT(0) space X, with the extra
information that S can be taken to have compact support.

Proof. Let X be proper CAT(0) and let T in I_n(X) be a cycle with compact
support. If T = 0 take S = 0. Otherwise pick a point o and a radius R so
large that spt T is contained in the open ball B(o,R). Let
Y = closure of B(o,R). Y is closed and bounded, hence compact (properness),
and it is convex (section 1.2(b)), so Y with the restricted distance is a
compact CAT(0) space.

Define the radial projection P: X -> Y by

    P(x)  =  [o,x]_lambda ,     lambda = min{ 1, R / d(o,x) } ,          (5.1)

that is, walk from o towards x and stop at distance R from o. On Y, P is the
identity. On all of X, P is 1-Lipschitz, and here is the proof, which is the
only place where the equal-fraction comparison (1.3) is used with different
factors at the same time.

Let x, x' in X, put r = d(o,x), r' = d(o,x'), and let lambda, lambda' be the
corresponding fractions from (5.1). By (1.3),

    d(P x, P x')^2
       <= ( lambda r - lambda' r' )^2
          + lambda lambda' ( d(x,x')^2 - (r - r')^2 ) .                  (5.2)

Now lambda r = min{r,R} and lambda' r' = min{r',R}, so
|lambda r - lambda' r'| <= |r - r'| because s |--> min{s,R} is 1-Lipschitz.
Also 0 <= lambda lambda' <= 1. Substitute into (5.2):

    d(P x, P x')^2
       <= (r - r')^2 + lambda lambda' ( d(x,x')^2 - (r - r')^2 )
       =  lambda lambda' d(x,x')^2 + (1 - lambda lambda')(r - r')^2
       <= lambda lambda' d(x,x')^2 + (1 - lambda lambda') d(x,x')^2
       =  d(x,x')^2 ,

where the last inequality is the triangle inequality |r - r'| <= d(x,x').
So P is 1-Lipschitz.

Consequences. P_#T is an integral n-cycle in Y with M(P_#T) <= M(T) (use
[lit 1](c) with Lip(P) <= 1). Apply the hypothesis in Y: there is an
integral (n+1)-current S in Y with dS = P_#T and
M(S) <= c_n M(P_#T)^((n+1)/n) <= c_n M(T)^((n+1)/n). Let iota: Y -> X be the
inclusion and put S' = iota_#S. Then M(S') = M(S) because iota is an
isometric embedding, dS' = iota_#(P_#T) = (iota o P)_# T, and iota o P is the
identity on the open set B(o,R) which contains spt T. By locality of
pushforward [lit 1](c), (iota o P)_# T = T. So dS' = T and spt S' is a
compact subset of Y. End of proof of [L6].

Remark. The lemma is stated as an implication between theorems in one fixed
dimension n. In the proof of [T1] it is applied at the very end, after the
compact case is settled for all n, so nothing circular happens.

5.2  A failed inequality produces a champion.

From now on: Y is a compact CAT(0) space, n >= 2, q_0 = (n+1)/n, V is the
filling volume (2.5) of section 2.2, and c > 0 is a constant. Define

    Phi(B)  =  V(B) - c M(B)^(q_0) ,                                   (5.3)

initially on integral n-cycles B in Y. The point of Phi is that it is
scale-invariant up to a power: if all distances are multiplied by a, both
terms get multiplied by a^(n+1) (see 5.3 below). So the champion of Phi is a
cycle whose shape, not whose size, is extremal.

[L7] (positive champion, and its first variation). Suppose Phi is positive
at some integral n-cycle. Then:
  (a) Phi attains a positive maximum on the set of integral n-cycles;
  (b) every maximizer A satisfies m := M(A) > 0 and the champion
      inequality

          V(A - B)  >=  c ( m^(q_0) - M(B)^(q_0) )
                                 for every integral n-cycle B .          (5.4)

Proof. (a) Boundedness above. By (2.21), V(B) <= L_n M(B) for a constant L_n
depending only on n and Y. With s = M(B),

    Phi(B) <= L_n s - c s^(q_0) =: psi(s) ,

and psi is bounded above on [0, infinity) because q_0 = (n+1)/n > 1. So the
supremum of Phi is finite; call it Sup > 0 by hypothesis.

Boundedness of the maximizing sequence. If Phi(B) > 0 then
c s^(q_0) < V(B) <= L_n s, hence s^(q_0 - 1) < L_n / c, hence

    M(B) < ( L_n / c )^n ,      since   q_0 - 1 = 1/n .                 (5.5)

So every cycle with Phi(B) > 0 has mass below a fixed number, and a
maximizing sequence may be chosen inside that set.

Compactness. Let B_j be a maximizing sequence: all masses below (L_n/c)^n,
all boundaries zero. By [lit 4] (compactness in a compact space) a
subsequence, again written B_j, converges weakly to an integral current A
with dA = 0. By [lit 7] the filling volume is continuous along this
convergence, so V(B_j) -> V(A); and by [lit 4] mass is lower
semicontinuous, so M(A) <= liminf M(B_j), hence

    M(A)^(q_0) <= liminf_j M(B_j)^(q_0) .

Combining, and using that s |--> -c s^(q_0) is monotone decreasing,

    Phi(A)  =  V(A) - c M(A)^(q_0)
            >=  limsup_j ( V(B_j) - c M(B_j)^(q_0) )
            =   Sup .

So A attains the supremum, and in particular M(A) > 0 because Sup > 0.

(b) First variation. Let B be any integral n-cycle. Maximality says
Phi(A) >= Phi(B), that is

    V(A) - V(B)  >=  c ( M(A)^(q_0) - M(B)^(q_0) ) .

Now A - B is a cycle, and we need subadditivity of the filling volume: if S
fills A - B and R fills B, then S + R is an integral current whose boundary
is (A - B) + B = A, so M(S) + M(R) >= V(A). Take the infimum over S and over
R separately:

    V(A)  <=  V(A - B) + V(B) .                                        (*)

Combine (*) with the maximality inequality above:

    c ( m^(q_0) - M(B)^(q_0) )  <=  V(A) - V(B)  <=  V(A - B) ,

which is (5.4). End of proof of [L7].

How (5.4) is used. Take B to be a small perturbation of A, for example a
copy of A pushed inwards along radial lines. Then the left side of (5.4) is
approximately c q_0 m^(1/n) times the amount of mass that A loses, and the
right side is the mass swept out by the homotopy. That is exactly the
comparison performed in section 6, and the normalization below makes the
coefficient c q_0 m^(1/n) equal to 1/n, which is the coefficient appearing
in the radial inequality [P1].

5.3  Rescaling to the subcritical mass.

Distances in a metric space can be multiplied by a positive number a; call
the resulting space aY, and the correspondingly scaled copy of a current
aB. Two facts, both immediate from the definitions:

    (i) the mass of a j-dimensional current is multiplied by a^j;
    (ii) the filling volume V is multiplied by a^(n+1).

Fact (i) deserves a line. Let phi: Y -> aY be the identity map between the
two copies of the same set; it is Lipschitz with constant a from Y to aY, and
with constant 1/a from aY to Y. By [lit 1](c), M(phi_# C) <= a^j M(C) and
M(C) = M(phi^(-1)_# phi_# C) <= a^(-j) M(phi_#C). Squeeze. Fact (ii) follows
because the correspondence S <--> aS is a mass-preserving bijection between
fillings of B and fillings of aB, with masses scaled by a^(n+1).

[L8] (normalization). Let A be a champion as in [L7], with c > c_n and
m = M(A) > 0. Then there is a factor a > 0 such that the rescaled champion
A' = aA satisfies

    m' := M(A')  =  ( (n+1) c )^(-n)  <  s_n                         (5.6)

and, automatically,

    c q_0 ( m' )^(1/n)  =  1/n .                                     (5.7)

Moreover, after the rescaling, A' is still a maximizer of the
corresponding functional in the rescaled space, so (5.4) holds for A' with
m' in place of m.

Proof. Choose a = ( (n+1) c m^(1/n) )^(-1). Then
m' = a^n m = m / ( (n+1) c m^(1/n) )^n = m / ( (n+1)^n c^n m ) =
( (n+1)c )^(-n), which is (5.6). The inequality m' < s_n is equivalent to
((n+1)c)^(-n) < ((n+1)c_n)^(-n), that is c > c_n, which is our hypothesis;
recall from section 1.5 that ((n+1)c_n)^(-n) = s_n exactly. For (5.7):

    c q_0 (m')^(1/n) = c ((n+1)/n) * ( (n+1)c )^(-1) = 1/n .

Finally, by (i) and (ii), Phi is multiplied by a^(n+1) under rescaling, so a
maximizer stays a maximizer. End of proof of [L8].

We now relabel: the ambient space is still written Y, the champion is
written A, and m = M(A) satisfies (5.6) and (5.7). The mass measure of A is

    mu  =  ||A|| ,      mu(Y) = m < s_n .                             (5.8)

This is the running notation for the rest of the proof. Note what (5.6) says
numerically: the champion is a cycle whose total weighted area is strictly
less than the area of the unit sphere in the same dimension. Everything that
follows is an attempt to prove m >= s_n, which will be a contradiction.

Exercise (why c > c_n is available). Suppose the sharp inequality (3.2)
fails for some cycle T, so that V(T) > c_n M(T)^((n+1)/n). Show that then
there is c with c_n < c < V(T) / M(T)^((n+1)/n), and that for this c the
functional (5.3) is positive at T. (This is the only way the proof of [T3]
starts.)


================================================================================
6.  RADIAL VARIATION: THE INEQUALITY [P1]
================================================================================

This section contains the one computation that uses the CAT(0) hypothesis
directly, through the comparison inequality (1.3). Everything else in the
proof is analysis (Sobolev inequalities, minimizers) or algebra
(determinants). So if you want to see why curvature <= 0 is the right
hypothesis, read 6.2.

Setting for the whole section: Y is a compact CAT(0) space, n >= 2, A is the
champion of [L7] normalized by [L8], so that

    mu = ||A|| ,      m = mu(Y) < s_n ,      c q_0 m^(1/n) = 1/n ,      (6.1)

with q_0 = (n+1)/n. Fix a centre y in Y and put

    r(x) = d(y,x) .                                                     (6.2)

The function r is 1-Lipschitz, so its chart gradient Dr exists mu-almost
everywhere in the sense of section 2.5 and satisfies |Dr| <= 1 there. (The
null set may depend on y; that is harmless and is mentioned again in 7.8.)

6.1  The deformation, and why it is Lipschitz.

Fix a nonnegative Lipschitz function g on Y, and choose h > 0 so small that

    h M_g <= 1/2 ,      M_g := sup g .                                  (6.3)

For 0 <= t <= h define

    lambda(t,x) = 1 - t g(x) ,      F(t,x) = [y,x]_{lambda(t,x)} ,      (6.4)

so that 1/2 <= lambda <= 1, F(0,x) = x, and each point x slides towards y as
t grows, by a fraction which depends on x. Write F_t(x) = F(t,x), and define
the moved cycle and the swept current

    B_h = (F_h)_# A ,       D_h = F_# ( [0,h] x A ) .                   (6.5)

Both are integral currents by [lit 1](c) and [lit 3]. Their boundaries are,
by (2.19) and dA = 0,

    dB_h = 0 ,       dD_h = B_h - A .                                   (6.6)

So D_h fills A - B_h, and this is the filling that will be compared with the
champion inequality (5.4).

F is jointly Lipschitz on [0,h] x Y. Here is the estimate, and it is needed
to justify everything after it. Let R = diam(Y), let t,t' in [0,h] and
x,x' in Y, and set lambda = 1 - t g(x), lambda' = 1 - t' g(x'). By the
two-point comparison (1.3),

    d(F(t,x), F(t',x'))^2
       <= ( lambda r - lambda' r' )^2
          + lambda lambda' ( d(x,x')^2 - (r - r')^2 ) ,                 (6.7)

where r = d(y,x), r' = d(y,x'). Now

    |lambda r - lambda' r'| <= |r - r'| + R |t g(x) - t' g(x')|
                            <= |r - r'| + R h Lip(g) d(x,x') + R M_g |t-t'|,

and 0 < lambda lambda' <= 1. Substituting in (6.7) and taking square roots,

    d(F(t,x), F(t',x'))
       <= 2 ( 1 + R h Lip(g) ) d(x,x') + 2 R M_g |t - t'| .             (6.8)

So F is Lipschitz with spatial constant 2(1 + R h Lip(g)) and time constant
2 R M_g. (The factor 2 is not important; what matters is that the constant
depends only on g, R and h, and not on the chart or on the point.)

6.2  The comparison forms.

We now look at the derivative of F. There are no charts in Y, so we work in
one chart at a time: let phi_i: E_i -> Z_i be the charts of (2.4), with
metric matrices G_i from [L1] and values r, g evaluated at phi_i(z), and let
alpha and gamma be the chart covectors

    alpha = D(r o phi_i)(z) ,      gamma = D(g o phi_i)(z) .            (6.9)

By (2.7) analogue for covectors, |alpha|_{G_i^(-1)} <= 1 and
|gamma|_{G_i^(-1)} <= Lip(g). Write lambda = 1 - t g, and set

    P_tot = lambda^2 ( G_i - alpha(x)alpha ) + u u^T ,
            u = lambda alpha - t r gamma - r g dt ,                       (6.10)
    P_sp  = lambda^2 ( G_i - alpha(x)alpha )
            + ( lambda alpha - t r gamma )( lambda alpha - t r gamma )^T . (6.11)

Here "alpha alpha" means the rank-one matrix alpha_i alpha_j, and in P_tot
the first summand acts only on the spatial directions w, while u has one
time component and one spatial component. Both forms are nonnegative:

    G_i - alpha alpha  >=  0   because G_i^{-1}[alpha,alpha] = |alpha|^2 <= 1,

and a sum of a nonnegative form and a square of a linear functional is
nonnegative.

The point of these forms is the following differential bound. For increments
(s,w) with (t+s, z+w) in [0,h] x E_i,

    d( F_i(t+s, z+w), F_i(t,z) )^2
        <= P_tot[(s,w),(s,w)] + o( |s|^2 + |w|^2 ) ,                     (6.12)

where F_i = F(t, phi_i(z)), and similarly with s = 0 for P_sp. Let us check
the derivation, because this is exactly where the CAT(0) hypothesis enters
and there is nothing else like it in the proof.

Let A_lam(t,z) = lambda(t,z) r(z), a scalar function on [0,h] x E_i, where
r(z) means r(phi_i(z)). Its differential is

    dA_lam(s,w) = ( lambda alpha - t r gamma ) . w - ( r g ) s + o(|s|+|w|) ,

because d(lambda) = -s g - t (gamma . w) and d r = alpha . w, and
lambda r = (1-tg) r. Now apply (6.7) with x = phi_i(z), x' = phi_i(z+w),
t replaced by t and t+s:

    d(F_i, F_i)^2 <= ( (1+o(1)) dA_lam(s,w) )^2
                     + lambda (lambda + O(s) ) ( d(phi(z+w),phi(z))^2
                                              - (r(z+w)-r(z))^2 ) .

Using the centered metric differential [L1](c), the bracket becomes
(G_i w . w) - (alpha . w)^2 + o(|w|^2). Also lambda(lambda+O(s)) = lambda^2 +
O(s). Since all terms are quadratic, the O(s) correction contributes only to
order o(|s|^2+|w|^2) times a bounded coefficient (the coefficients are
bounded because Lip(g), R, Lip(phi_i), |alpha| are). Thus

    d^2 <= ( dA_lam(s,w) )^2 + lambda^2 ( G_i - alpha alpha )[w,w]
            + o(|s|^2 + |w|^2) ,

which is (6.12). Setting s = 0 gives the same statement with
(lambda alpha - t r gamma) in place of u.

Determinant of the spatial form. Work in coordinates in which G_i(z) is the
identity, and remember that at the end we divide by J_i = sqrt(det G_i).
From (6.11),

    P_sp = lambda^2 I - lambda t r ( alpha gamma^T + gamma alpha^T )
             + t^2 r^2 gamma gamma^T
         = I - t ( 2 g I + r ( alpha gamma^T + gamma alpha^T ) )
             + O(t^2) ,

where we substituted lambda = 1 - tg and collected: lambda^2 I = I - 2tgI +
O(t^2), and the t^2 term of lambda^2 combined with the last term is O(t^2).
Now use that the derivative at the identity of H |--> sqrt(det H) is
H' |--> (1/2) trace(H'). The trace of (alpha gamma^T + gamma alpha^T) is
2 alpha . gamma = 2 Dr . Dg in the chart metric. Hence

    sqrt(det P_sp) / J_i
       = 1 - t ( n g + r Dr . Dg ) + O(t^2) .                           (6.13)

The remainder is bounded by C t^2 with C depending only on n, R, M_g and
Lip(g), uniformly in the chart and the point: after G-orthonormalisation
every entry of the matrices involved is controlled by those four numbers,
because |alpha| <= 1, |gamma| <= Lip(g), r <= R, g <= M_g. That uniformity
is used in 6.3 and is the reason the O(h^2) there is honest.

Determinant of the total form. In time-first coordinates, (6.10) reads

    P_tot = [[0, 0], [0, H]] + u u^T ,    H = lambda^2 ( I - alpha alpha^T ) ,
    u = ( -r g , lambda alpha - t r gamma ) .

Write q = lambda alpha - t r gamma, a = -r g, so that

    P_tot = [[ a^2 , a q^T ] , [ a q , H + q q^T ]] .

If a != 0, the Schur complement of the block a^2 gives

    det P_tot = a^2 det( H + q q^T - q q^T ) = a^2 det H .

If a = 0 both sides are zero, because the first row and column vanish and
the matrix has dimension n+1 >= 3. Finally
det H = lambda^(2n) det(I - alpha alpha^T) = lambda^(2n)(1 - |alpha|^2).
Taking square roots,

    sqrt(det P_tot) / J_i
       = r g lambda^n sqrt( 1 - |Dr|^2 ) .                              (6.14)

This is an exact identity, not an estimate: the comparison form has the
sphere-like factor sqrt(1 - |Dr|^2), which is the classical "transverse
component" of the radial direction flat-space computation. Note that (6.14)
is valid also when r g = 0 (then both sides are 0) and when |Dr| = 1.

Figure-free version, in words: the spatial form loses mass at rate
n g + r Dr.Dg because the map shrinks by the factor lambda; the total form
sweeps mass at rate r g sqrt(1-|Dr|^2) because a radial motion has volume
weight equal to the component of the velocity transverse to the current.
Both rates are exactly what one would compute in R^n, and that is the
content of (6.13) and (6.14).

6.3  From the forms to mass.

We now apply [L2] twice. First write the action formulas in a chart, which
is just the chart formula (2.1) composed with F. For bounded Lipschitz b and
Lipschitz pi_1,...,pi_n,

    B_h(b, pi) = sum_i int_{E_i} theta_i(z) b(F_i(h,z))
                 det D_z( (pi o F_i(h,.)) )(z) dz ,                     (6.15)

and for the homotopy,

    D_h(b, pi_1,...,pi_{n+1})
       = sum_i int_0^h int_{E_i} theta_i(z) b(F_i(t,z))
         det D_{(t,z)}( (pi_1 o F_i, ..., pi_{n+1} o F_i) )(t,z) dz dt . (6.16)

These two formulas are the chart form of [lit 3] applied to the whole
current A, with the pushforward rules of [lit 1]. Two comments, both from the
source of this argument and both worth knowing:

  (a) The interchanges of sum and integral in (6.16) are legitimate: for a
      fixed test tuple, every cofactor in the expansion of the determinant is
      bounded by a constant times J_i, by the Euclidean determinant bound; so
      the total absolute integral is at most C h m by (2.13), and Fubini
      applies.

  (b) In (6.16), the separate time and space derivatives agree almost
      everywhere with the joint derivatives of the Lipschitz compositions, so
      the formula is the honest determinant of the joint derivative.

Now (6.12) says precisely that the quadratic majorant hypothesis (2.9) of
[L2] holds, for B_h with P_i = P_sp(h,z), and for D_h with
P_i = P_tot(t,z). The two majorant measures (2.10) are finite:

    by (6.13),  sqrt(det P_sp) <= (1 + C h) J_i ,
    by (6.14),  sqrt(det P_tot) <= R M_g J_i ,

so their integrals are bounded by (1+Ch) m and h R M_g m respectively, using
(2.13). [L2] therefore gives

    M(B_h)  <=  m - h int_Y ( n g + r Dr . Dg ) dmu + C h^2 ,           (6.17)

    M(D_h)  <=  h int_Y r g sqrt(1 - |Dr|^2) dmu .                      (6.18)

Both inequalities use only (6.13), (6.14) and the exact chart mass (2.13).
Note what (6.17) says: the moved cycle lost mass at the rate
n g + r Dr.Dg. Note what (6.18) says: the swept current has mass at most the
integral of r g sqrt(1-|Dr|^2).

6.4  The radial inequality [P1].

Everything is now in place for the champion comparison.

[P1] (radial first variation). Let Y be a compact CAT(0) space, A the
champion of (6.1) with mass measure mu and mass m < s_n, y in Y, and
r(x) = d(y,x). Then for every nonnegative Lipschitz g on Y,

    int_Y ( n g + r Dr . Dg ) dmu
        <=  n int_Y r g sqrt(1 - |Dr|^2) dmu .                           (6.19)

The square root is defined mu-almost everywhere, since |Dr| <= 1.

Proof. Abbreviate

    I_g = int ( n g + r Dr.Dg ) dmu ,      H_g = int r g sqrt(1-|Dr|^2) dmu .

By (6.17) there is C = C(g,y) with M(B_h) <= U_h := m - h I_g + C h^2 for
all small h; since m > 0, also U_h > 0 for small h. The current -D_h fills
A - B_h (because dD_h = B_h - A). So the champion inequality (5.4) applies
with B = B_h and gives

    c ( m^(q_0) - M(B_h)^(q_0) )  <=  V(A - B_h)  <=  M(D_h)  <=  h H_g .

Since U_h >= M(B_h) and s |--> s^(q_0) is increasing,

    c ( m^(q_0) - U_h^(q_0) )  <=  c ( m^(q_0) - M(B_h)^(q_0) )  <=  h H_g .

Divide by h > 0 and let h -> 0. By (6.17), ( m^(q_0) - U_h^(q_0) ) / h ->
q_0 m^(q_0 - 1) I_g = q_0 m^(1/n) I_g (using q_0 - 1 = 1/n). Hence

    c q_0 m^(1/n) I_g  <=  H_g .

The normalization (6.1) says c q_0 m^(1/n) = 1/n, so I_g <= n H_g, which is
(6.19). End of proof of [P1].

Three comments.

  (a) [P1] is an inequality between two integrals over the current, with no
      reference to the space's dimension and no reference to charts. It is the
      complete output of the geometric part of the proof, together with the
      base case below.

  (b) The champion inequality enters only through the numeric identity
      c q_0 m^(1/n) = 1/n of (6.1); that is why the rescaling [L8] was worth
      doing. Without it the inequality would read
      c q_0 m^(1/n) I_g <= H_g with an unhelpful coefficient.

  (c) The function g is an arbitrary nonnegative Lipschitz function, and the
      centre y is arbitrary. In section 6.5 g is a cutoff 1/max(r,eps)^2; in
      section 7.7 it is a more elaborate function built from the minimizer.
      The freedom to choose g is what makes one inequality enough for two
      very different jobs.

Exercise. Put g = 1 identically and y any point. Show that (6.19) reads
n m + int r Dr.D(1) <= n int r sqrt(1-|Dr|^2) dmu, that is
m <= int r sqrt(1-|Dr|^2) dmu. Conclude m <= R m trivially, i.e. this
choice of g gives no information. (Moral: g must concentrate mass.)

6.5  The base case: 2-cycles, by an explicit cutoff.

This is the one place where the famous constant 4 pi = s_2 appears from a
line-by-line computation, and it is worth reading carefully because it shows
what the integers in "integral current" are for.

[P2] (sharp filling in cycle dimension two). Let Z be a compact CAT(0)
space and let T in I_2(Z) be a cycle. Then there exists S in I_3(Z) with

    dS = T ,        M(S) <= c_2 M(T)^(3/2) = M(T)^(3/2) / (6 sqrt(pi)) .  (6.20)

Proof. Suppose not. Choose c > c_2 so that Phi = V - c M^(3/2) is positive
somewhere, and let A be the champion of [L7], rescaled by [L8] so that
m = M(A) < s_2 = 4 pi and c (3/2) m^(1/2) = 1/2. Fix a Lebesgue point y of
the mass measure mu = ||A||, in the sense of [L4] with k = 2:

    liminf_{eps -> 0} eps^(-2) mu(B(y,eps))  >=  omega_2 = pi .         (6.21)

Such a y exists because [L4] holds mu-almost everywhere and m > 0. Put r(x) =
d(y,x) and, for 0 < eps < 1,

    g(x) = max{ r(x), eps }^(-2) .                                      (6.22)

This is Lipschitz and nonnegative, so it is admissible in [P1]. Two
computations.

Inside the ball, r <= eps, g = eps^(-2) is constant, so Dg = 0 and

    n g + r Dr.Dg = 2 eps^(-2) ,      r g sqrt(1 - |Dr|^2) = 0 ,        (6.23)

the second because |Dr| = 1 mu-almost everywhere in the ball (r is a distance
function, so it moves at unit speed in its own gradient direction; this is
[lit 1](e) applied to the level sets of r).

Outside the ball, r > eps, g = r^(-2) and the chain rule gives
Dg = -2 r^(-3) Dr. Then with b = r^(-1) sqrt(1 - |Dr|^2) >= 0,

    2 g + r Dr.Dg = 2 r^(-2) ( 1 - |Dr|^2 ) = 2 b^2 ,                   (6.24)

    r g sqrt(1-|Dr|^2) = r^(-1) sqrt(1-|Dr|^2) = b .                    (6.25)

Since the level set {r = eps} has mu-measure zero for almost every eps (there
are only countably many levels of positive measure under a finite measure),
we may assume it for our eps. Inserting (6.23)-(6.25) into the radial
inequality [P1] gives

    2 eps^(-2) mu(B(y,eps))  <=  2 int_{r>eps} ( b - b^2 ) dmu .         (6.26)

Now 2(b - b^2) = 1/2 - 2(b - 1/2)^2 <= 1/2 pointwise, so the right side of
(6.26) is at most m/2. Hence

    mu(B(y,eps))   <=  eps^2 m / 4 .                                    (6.27)

Combine (6.21) with (6.27): divide by eps^2 and let eps -> 0 to get
pi <= m/4, that is m >= 4 pi = s_2. That contradicts m < s_2. So no champion
exists, and the filling inequality holds in dimension two with the sharp
constant c_2. End of proof of [P2].

Two remarks. First, the factor 4 is forced by the pointwise inequality
b - b^2 <= 1/4, with equality exactly at b = 1/2; that is the value a flat
two-dimensional sheet takes at distance 2r from the centre. The base case
therefore detects the flat model and nothing else, which is why it delivers
the sharp constant rather than just some constant. Second, the algebra:
(6.21) says mu(B(y,eps)) >= (pi + o(1)) eps^2 and (6.27) says
mu(B(y,eps)) <= eps^2 m/4; divide by eps^2 and let eps -> 0.

================================================================================
7.  THE INDUCTION STEP: ANALYSIS ON THE CHAMPION
================================================================================

From now on n >= 3, the space Y is compact CAT(0), and A is the champion of
(6.1) with mass measure mu, mass m < s_n, and normalization
c q_0 m^(1/n) = 1/n. We assume (H_(n-1)), the sharp filling inequality in
cycle dimension n-1, and we aim for a contradiction.

The plan, in one sentence per step.

    7.1  Cut A into small pieces. Each piece has small mass, so by (H_(n-1))
         its (n-1)-boundary is filled cheaply; feeding that into the champion
         inequality gives a sharp perimeter bound for small pieces. [P3]
    7.2  Convert that bound into a Sobolev inequality: coarea turns perimeter
         into a derivative, rearrangement turns the current into a Euclidean
         radial function. The two bookkeeping lemmas are [L9] and [L10].
    7.3  On R^n the radial Sobolev inequality is sharp by mass transport
         [L12]; the rearrangement does not lose energy. [L11], [L12]
    7.4  Result: the almost sharp Sobolev inequality on Y. [P4]
    7.5  Analysis set-up: the closed gradient [P5] and its compactness [P6].
    7.6  Variational set-up: a minimizer v exists and is subcritical, and
         its powers are admissible tests. [P7], [L13], [L14]
    7.7  Testing with a Euclidean bubble profile gives the chord transform
         bound [P8]. [L15], [L16]
    7.8  Tangent density bounds the chord moments above [L17]; variance
         bounds their average below [L18]. Contradiction. [T3]

7.1  Small pieces of A have Euclidean perimeter.

The only place where (H_(n-1)) enters.

[P3] (uniform small-set isoperimetry). For every eta in (0,1) there is
delta in (0, m/2), depending only on eta, n, m and the space, such that the
following holds. Let u: Y -> R be Lipschitz and put

    a(t) = mu{ u > t } ,    U_t = A|{u > t} ,    P(t) = M( d U_t ) ,     (7.1)

where U_t and its boundary are integral for almost every t (by [lit 1](e)).
Then for almost every level t with a(t) <= delta,

    P(t)  >=  (1 - eta) n omega_n^(1/n) a(t)^((n-1)/n) .                (7.2)

Proof. Fix a level with U_t integral, write a = a(t), P = P(t). By
(H_(n-1)) applied to the integral (n-1)-cycle dU_t, there is an integral
n-current R in Y with

    dR = dU_t ,      d_R := M(R) <= ( P / (n omega_n^(1/n)) )^(n/(n-1)) (7.3)

This is exactly the sharp filling inequality in dimension n-1, solved for
P in terms of the filling mass. Note that d_R is a monotone function of P,
which is why (7.3) is the right way round.

Now B := A - U_t + R is a cycle, and so is A - B = U_t - R. Masses: by
restriction (2.14) the mass of A splits exactly, so

    M(B) <= m - a + d_R ,      M(A - B) <= a + d_R ,                    (7.4)

and no disjointness of R from U_t is needed for these bounds. The champion
inequality (5.4) applies to B, and the coarse linear bound (2.21) applies to
A - B:

    c q_0 (m/2)^(1/n) (a - d_R)  <=  c ( m^(q_0) - M(B)^(q_0) )
                                <=  V(A - B)  <=  L_n M(A-B)^(q_0)
                                <=  L_n (2a)^(q_0) .                    (7.5)

The first inequality used the mean value theorem for s^(q_0) on the interval
[m - a + d_R, m] and the bound s^(1/n) >= (m/2)^(1/n) there, which holds
because a + d_R <= m (we will have this for small a). Solving (7.5) for
d_R gives

    d_R  >=  a ( 1 - C a^(1/n) ) ,      C = n K_n' 2^(1+2/n) ,          (7.6)

with a constant C depending only on n and the space; if d_R >= a already,
(7.6) holds trivially. Now combine (7.6) with (7.3), solved for P, and choose
delta so small that C delta^(1/n) <= 1 - (1-eta)^(n/(n-1)). Then for
a <= delta,

    P  >=  n omega_n^(1/n) a^((n-1)/n) ( 1 - C a^(1/n) )^((n-1)/n)
       >=  (1 - eta) n omega_n^(1/n) a^((n-1)/n) ,

which is (7.2). The case a = 0 is immediate. End of proof of [P3].

Read (7.2) as: small pieces of the champion cannot be cheap. Their boundary
must have at least the Euclidean sphere area, up to the factor 1-eta. The
loss eta is the price of the fact that we only know the induction hypothesis
in dimension n-1 and not a sharp value for how R sits relative to U_t.

7.2  Two bookkeeping lemmas about slices.

Coarea has to be done with the chart gradient, not with the Lipschitz
constant, otherwise the Dirichlet energy is lost. To do that, one needs to
recover the total mass of a slice from countably many scalar evaluations, and
that is [L9].

[L9] (countable family of mass tests). For a compact metric space Z and
k >= 1 there is a countable family T_k(Z) of finite lists of pairs
((b_1,pi^1),...,(b_N,pi^N)), where the b_j are Lipschitz with sum |b_j| <= 1
and each pi^j is a k-tuple of 1-Lipschitz functions, such that every
finite-mass k-current C on Z satisfies

    M(C)  =  sup_{T in T_k(Z)} sum_{(b,pi) in T} C(b, pi) .             (7.7)

Sketch. Normalize the tuples by requiring pi(z_0) = 0 for a fixed base
point; the normalized 1-Lipschitz functions form a compact metric space in
the uniform topology (Arzela-Ascoli), so it has a countable dense subset.
Let C have mass measure lambda, and let rho_pi be the density of the signed
measure b |--> C(b, pi) with respect to lambda. These densities satisfy
|rho_pi| <= 1 and, because the dense family approximates every tuple,
lambda is controlled by the measure rho lambda with rho = sup_j |rho_{pi^j}|;
minimality of lambda forces rho = 1 almost everywhere. Selecting, at each
point, the index of the largest |rho_{pi^j}| and its sign gives coefficients
b_j with sum |b_j| <= 1 and sum_j C(b_j,pi^j) = int max_j |rho_{pi^j}| d
lambda. Letting the finite families grow gives (7.7). The coefficients may be
taken Lipschitz by approximating Borel functions in L^1(lambda) and
normalizing the tuple. End of sketch of [L9].

[L10] (intrinsic coarea). With the notation (7.1), set

    e(t) = int_{{u > t}} |Du|^2 dmu .

Then for almost every t,

    P(t)^2  <=  ( -a'(t) ) ( -e'(t) ) ,      P(t) <= Lip(u) ( -a'(t) ) ,  (7.8)

where the derivatives are the densities of the absolutely continuous parts of
the distribution functions a and e; neither a nor e is assumed to be
absolutely continuous.

Sketch. For a bounded interval I and H(s) = int_I 1{t < s} dt, the slice
identity [lit 5] integrated over I, together with the cycle identity dA = 0,
gives

    int_I dU_t(b, pi) dt = A( H(u), b, pi ) = -A( b, H(u), pi ) .        (7.9)

The chain rule in charts gives D(H(u)) = 1_{u in I} Du, and the differential
of a Lipschitz function vanishes on its own level sets, so no endpoint
corrections appear. For each admissible list from [L9], (7.9) and the
determinant bound give

    | int_I sum C(b_j,pi^j) dt |  <=  int_{{u in I}} |Du| dmu .

Taking the supremum over the countable family ([L9]) and differentiating
with respect to the interval gives P(t) <= (d/dt) gamma_ac(t), where
gamma = u_#(|Du| mu). Then Cauchy-Schwarz on each level piece:
gamma(I)^2 <= lambda(I) zeta(I) with lambda = u_#mu, zeta = u_#(|Du|^2 mu),
and dividing by the interval length squared and differentiating gives the
second inequality of (7.8). End of sketch of [L10].

7.3  Rearrangement, and the sharp radial Sobolev inequality.

The pieces of A are now known to have Euclidean perimeter. Two more transfers
are needed to turn that into a Sobolev inequality; both are classical, and
both have to be done carefully because the mass measure is not Lebesgue
measure and the level function a(t) is not continuous.

[L11] (rearranged energy). Fix eta, let delta come from [P3], and let u >= 0
be Lipschitz with mu{u > 0} <= delta. Then there is a nonnegative, radially
nonincreasing, compactly supported Lipschitz function u* on R^n such that

    int_{R^n} (u*)^q dx = int_Y u^q dmu  for all q > 0 ,   and            (7.10)

    int_{R^n} |grad u*|^2 dx  <=  (1 - eta)^(-2) int_Y |Du|^2 dmu .     (7.11)

Sketch. For each level t, define R(t) by omega_n R(t)^n = a(t): the radius of
the Euclidean ball of the same mass. [P3] and [L10] give

    -R'(t)  >=  (1 - eta) / Lip(u) ,

which also controls the singular part of R. Let h be the generalized inverse
of R and set u*(x) = h(|x|); then the superlevel sets of u* are balls of
volume a(t), so (7.10) is the layer cake formula. The energy computation is

    int |grad u*|^2 dx
       = int_0^U ( n omega_n R(t)^(n-1) )^2 / ( -a'(t) ) dt
       <= (1-eta)^(-2) int_0^U ( -e'(t) ) dt
       <= (1-eta)^(-2) int |Du|^2 dmu ,

where the middle inequality is (7.8). The substitution is valid even when a
has atoms or a singular part: an atom produces a central plateau of u*, and a
singular part contributes no energy, because the substitution uses the
absolutely continuous measure |h'| d rho. End of sketch of [L11].

[L12] (sharp radial Sobolev). For every nonnegative, radially nonincreasing,
compactly supported Lipschitz f on R^n with n >= 3,

    n s_n^(2/n) ( int f^p dx )^(2/p)  <=  ( 4/(n-2) ) int |grad f|^2 dx ,  (7.12)

with p = 2n/(n-2). Equality is approached by the bubbles
(1 + |x|^2)^(-(n-2)/2).

Sketch (this is the classical computation, and it is short). Normalize
int f^p = 1. Let g(x) = A (1+|x|^2)^(-(n-2)/2) be the bubble, normalized
the same way. Match the cumulative masses of f^p and g^p by a radial
transport T: [0,r_0) -> (0,infinity). Then

    dT := t'(r) ( t(r)/r )^(n-1)  =  f(r)^p / g(t(r))^p ,

and the arithmetic-geometric mean inequality on the n eigenvalues of DT,
t' and (n-1) copies of t/r, together with an integration by parts, gives

    n int g^ell  <=  ( int |grad f|^2 )^(1/2) ( int |x|^2 g^p )^(1/2) ,  (7.13)

where ell = p(1 - 1/n) = 2(n-1)/(n-2). For the bubble g the two sides of
(7.13) are both proportional to the same finite constant. Differentiating
-g'(r) = K r g(r)^(p/2) with K = (n-2) A^(-2/(n-2)) and integrating by parts
against the vector field x gives the two identities

    n int g^ell dx  =  ell K int |x|^2 g^p dx ,
    int |grad g|^2 dx  =  K^2 int |x|^2 g^p dx .

Dividing, the constant K and the moment cancel, and (7.13) becomes

    int |grad g|^2 dx  <=  int |grad f|^2 dx .

So the bubble has the least energy among all admissible f with the same L^p
mass. Its energy is n(n-2)/4 * s_n^(2/n), and int g^p = 1, so this is
exactly (7.12). The exact numbers come from the stereographic projection identity

    int_{R^n} (1 + |x|^2)^(-n) dx  =  2^(-n) s_n ,                       (7.14)

and from -Delta w = n(n-2)(1+|x|^2)^(-(n+2)/2) for
w = (1+|x|^2)^(-(n-2)/2). End of sketch of [L12].

7.4  Globalization: the almost sharp Sobolev inequality.

[P4] (almost sharp Sobolev). For every eps in (0, S_*) there is a finite
constant C_eps such that every real Lipschitz u on Y satisfies

    ( S_* - eps ) ( int |u|^p dmu )^(2/p)
       <=  4 beta int |Du|^2 dmu + C_eps int u^2 dmu .                  (7.15)

Here S_* = n s_n^(2/n), beta = 1/(n-2), p = 2n/(n-2).

Sketch. Combine [L11] and [L12]: for nonnegative z with mu{z > 0} <= delta,

    S_* (1 - eta)^2 ( int z^p dmu )^(2/p)  <=  4 beta int |Dz|^2 dmu .  (7.16)

For general Lipschitz u with ||u||_2 > 0 put tau = ||u||_2 / sqrt(delta) and
z = (|u| - tau)_+. By Chebyshev, mu{z > 0} <= delta; the scalars |u| and z
differ by at most tau, so || |u| - z ||_p <= m^(1/p) tau; and
E(z) <= E(u) because z is a scalar Lipschitz composition of |u| (the gradient
vanishes on the corner level set, so the usual chain rule holds almost
everywhere). Applying (7.15) to z and using the elementary inequality
(b+d)^2 <= (1+xi) b^2 + (1+1/xi) d^2 gives

    S_* (1-eta)^2 / (1+xi) ||u||_p^2
       <= 4 beta E(u) + ( S_* (1-eta)^2 m^(2/p) / (xi delta) ) ||u||_2^2 .

Choose eta, xi > 0 so small that S_* (1-eta)^2 / (1+xi) >= S_* - eps, and
set C_eps = S_* (1-eta)^2 m^(2/p) / (xi delta). That is (7.15). End of
sketch of [P4].

The constant C_eps is not small and is not meant to be: it can be absorbed
later because the quadratic term is subcritical. What matters is that the
coefficient of the critical term is S_* - eps, with S_* the sharp Euclidean
constant.

7.5  The closed gradient and its compactness.

The Sobolev inequality (7.15) is for Lipschitz functions. To make a
variational argument we need a Hilbert space and a minimizer, and that means
closing the gradient.

[P5] (Sobolev space and calculus). Identify real Lipschitz functions that
agree mu-almost everywhere, and complete them in

    ||u||_H^2  =  ||u||_2^2 + E(u) ,     E(u) = int |Du|^2 dmu .         (7.17)

Then:
  (a) the completion H is a Hilbert space, and its gradient D is a
      well-defined closed operator;
  (b) H embeds continuously in L^p(mu), and (7.15) holds on H;
  (c) H is stable under C^1 composition with bounded first derivatives;
  (d) Du = 0 mu-almost everywhere on each level set {u = t}; consequently
      the usual chain rule holds for piecewise C^1 scalars with finitely
      many corners;
  (e) u |--> u_+ is continuous from H to H, and every nonnegative element of
      H is approximated in H by nonnegative Lipschitz functions;
  (f) for each fixed y, the radial inequality [P1] holds on H, that is, for
      every nonnegative g in H.

The only point that is not standard is (a), and here is why it works and why
it is not obvious: two Lipschitz functions that agree almost everywhere could
in principle have different charts of different gradients, and nothing in the
definition forbids that.

Proof of (a) (graph closure). Suppose u_j are Lipschitz,
u_j -> 0 in L^2(mu), and Du_j -> L in L^2(mu) (the covariant Hilbert space
of section 2.5). For a Lipschitz b and a Lipshitz (n-1)-tuple pi, the cycle
identity dA = 0 and the product rule give

    A( b, u_j, pi )  =  -A( u_j, b, pi )  ->  0 ,                       (7.18)

because |A(u_j,b,pi)| <= const ||u_j||_1 -> 0. The left side converges to the
integral of b against the signed chart measure

    nu_pi(B) = sum_i int_{phi_i^(-1)(B)} theta_i
                  det( L_i, D(pi o phi_i) ) dz ,

which by the determinant bound is a finite signed measure with total
variation controlled by ||L||_1. So nu_pi = 0 for every admissible (b,pi).

Now fix a chart i, and let F^1,...,F^n be the components of phi_i^(-1),
extended to all of Y by McShane [lit 11]. Then at almost every point of the
chart, D(F^a o phi_i) = e_a^* (the standard unit covector in the domain). Use
as pi the n tuples obtained by deleting one component of (F^1,...,F^n); the
corresponding measures nu_pi vanish, so restricting them to Z_i gives
theta_i L_i = 0 almost everywhere on E_i. The multiplicity is nonzero, so
L_i = 0. There are countably many charts and finitely many tuples per chart,
so L = 0.

Thus the graph of D is closed, and its completion is a Hilbert space of
functions in L^2. The embedding (b) and the calculus (c)-(e) are standard
once the graph is closed; the level-set statement (d) is proved by
approximating the indicator of a level set by smooth cutoffs and using the
chain rule (i), so the corner cases follow by approximation. Finally (f) is
proved by approximating g in H by nonnegative Lipschitz functions (part e)
and passing to the limit in [P1]: every term of (6.19) is continuous in g for
the H-norm, because r and Dr are bounded. End of proof of [P5].

[P6] (compact embedding). The inclusion H -> L^2(mu) is compact.

Sketch. Let u_j be bounded in H, hence also in L^p by (b). Fix a chart i and
a Lipschitz cutoff chi on Y. Define the projected signed measure

    lambda_{j,chi}(B)  =  A( (1_B o F) chi u_j, F^1, ..., F^n ) ,
    B in R^n ,

where F is the chart's inverse coordinate map. These measures are supported
in the compact set F(Y), and their total variation is bounded uniformly in j.
The cycle identity gives a bound on the distributional derivatives of
lambda_{j,chi}: for smooth compactly supported h,

    | int partial_a h d lambda_{j,chi} |

is bounded by a constant times ||h||_infinity (the constant depends on chi
but not on j), because moving the factor h to the first slot and applying
[lit 1](d) gives a bound in terms of ||D(chi u_j)||_1 <= const.

Apply the BV compactness recipe of [lit 14]: translate, average against a
mollifier of scale eps, apply Ascoli-Arzela to the smoothed family (which has
common bounds on values and first derivatives), then let eps shrink. This
shows that {lambda_{j,chi} : j >= 1} is precompact in total variation. Now
approximate 1_{Z_i} by Lipschitz cutoffs chi_k (possible by regularity of mu),
and use the uniform bound

    || lambda_{j,chi_k} - lambda_{j, 1_{Z_i}} ||_TV  <=  C ||chi_k - 1_{Z_i}||_2

to conclude that {lambda_{j, 1_{Z_i}}} is also precompact in TV. On this
isolated chart the density is exactly 1_{E_i}(z) theta_i(z) u_j(phi_i(z)),
so

    ||u_j - u_k||_{L^1(Z_i,mu)}  <=  Lip(phi_i)^n
                                     || lambda_{j,1_{Z_i}}
                                         - lambda_{k,1_{Z_i}} ||_TV .

Take a diagonal subsequence that is Cauchy in L^1(Z_i, mu) for every i; the
tail sum is small because mu is finite. The subsequence is then Cauchy in
L^1(mu). It is also bounded in L^p (mu), and p > 2, so interpolation gives

    ||u_j - u_k||_2  <=  ||u_j - u_k||_1^alpha ||u_j - u_k||_p^{1-alpha} ,

with alpha = (p-2)/(2(p-1)) > 0, hence Cauchy in L^2. End of sketch of [P6].

7.6  The minimizer, and why its powers are tests.

The almost sharp inequality (7.15) has a strict gap: the quadratic term is
irrelevant at scale of the minimizer because the small mass m < s_n pushes
the critical coefficient below S_*. Concretely:

[P7] (subcritical minimizer). Define

    Q(u)  =  4 beta E(u) + n ||u||_2^2 ,     Lambda  =  inf Q(u) / ||u||_p^2 ,
    over nonzero u in H .                                                 (7.19)

Then 0 < Lambda < S_*. There is a nonzero nonnegative v in H minimizing
Q(u)/||u||_p^2, and

    0 < W := int v^p dmu  =  ( Lambda / n )^(n/2)  <  s_n ,             (7.20)

    Q(v)  =  n W ,                                                       (7.21)

and for every psi in H,

    4 beta int <Dv, Dpsi> dmu + n int v psi dmu  =  n int v^(p-1) psi dmu . (7.22)

Proof. Positivity of Lambda follows from (7.15) with eps = S_*/2:
(S_*/2) ||u||_p^2 <= 4 beta E(u) + C ||u||_2^2 <= C' Q(u). Upper bound: the
constant function u = m^(-1/p) has ||u||_p = 1 and Q(u) = n m^(1-2/p) =
n m^(2/n), which is < n s_n^(2/n) = S_* because m < s_n.

Let u_j be an L^p-normalized minimizing sequence. It is bounded in H (Q is a
positive definite quadratic form equivalent to ||.||_H^2). By Hilbert space
compactness and [P6], a subsequence converges weakly in H and strongly in
L^2, and almost everywhere, to an element u. Put t = int |u|^p.

If t < 1, consider z = u_j - u. Brézis-Lieb [lit 8] (or the elementary
inequality in section 4) gives int |z|^p -> 1 - t. Quadratic splitting gives
Q(u_j) = Q(u) + Q(z) + o(1), and by definition of Lambda,
Q(u) >= Lambda t^(2/p) and Q(z) >= (S_* - eps)(1-t)^(2/p) - C_eps||z||_2^2.
Since ||z||_2 -> 0,

    Lambda  >=  Lambda t^(2/p) + S_* (1 - t)^(2/p) ,

and the right side is strictly bigger than Lambda whenever t < 1 and
S_* > Lambda (use a^(2/p) >= a for a in [0,1]). So t = 1, u attains the
minimum, and ||u||_p = 1. Replacing by |u| preserves both ||.||_p and E(u) by
[P5](e), so we take u >= 0.

Differentiate Q(u+s psi)/||u+s psi||_p^2 at s = 0: by [P5], the gradient is
continuous in H, and the L^p-derivative is dominated by p(|u|+|psi|)^(p-1)
|psi|, which is integrable. Since ||u||_p = 1, this gives

    4 beta int <Du, Dpsi> dmu + n int u psi dmu
       = Lambda int u^(p-1) psi dmu  for all psi in H .

Set v = (Lambda/n)^(1/(p-2)) u. Then v also minimizes the quotient, and its
Euler equation becomes (7.22), because Lambda * (Lambda/n)^(2/(p-2)) =
(Lambda/n) * (Lambda/n)^(2/(p-2)) * n ... more simply: multiplying (7.22's
predecessor) by (Lambda/n)^((p-1)/(p-2)) gives (7.22). Its p-th moment is

    W  =  int v^p dmu  =  (Lambda/n)^(p/(p-2))  =  (Lambda/n)^(n/2) ,

where we used p/(p-2) = n/2. This is < s_n because Lambda < S_* =
n s_n^(2/n). Finally, testing (7.22) with psi = v gives
4 beta E(v) + n ||v||_2^2 = n W, which is (7.21). End of proof of [P7].

[L13] (powers of the minimizer). For every real gamma >= 1,
v^gamma in H and D(v^gamma) = gamma v^(gamma-1) Dv. Moreover,
int v^q dmu < infinity for every finite q > 0, and
int v^q |Dv|^2 dmu < infinity for every q >= 0.

Sketch. Assume temporarily v in L^(2 gamma). Put w_N = v min(v,N)^(gamma-1)
and psi_N = v min(v,N)^(2 gamma-2). Both are in H by [P5](c). For their
gradients,

    <Dv, Dpsi_N> >= d_gamma |Dw_N|^2 ,      d_gamma = (2 gamma - 1) / gamma^2 :

there is equality on {v < N}, and on {v > N} the ratio of coefficients is 1
and 0 < d_gamma <= 1; the gradient vanishes on the level {v = N} by [P5](d).
Test (7.22) with psi_N, drop the nonnegative term n int v psi_N, and get

    4 beta d_gamma E(w_N)  <=  n int v^(p-2) w_N^2 dmu .                 (7.23)

Now apply (7.15) with eps = S_*/2 to w_N and use (7.23):

    ( d_gamma S_* / 2 ) ||w_N||_p^2
       <=  n int v^(p-2) w_N^2 dmu + d_gamma C ||w_N||_2^2 .

Hölder on the "critical potential":

    int_{v > K} v^(p-2) w_N^2 dmu
       <=  (int_{v > K} v^p dmu)^((p-2)/p) ||w_N||_p^2 .

Since v^p is integrable, choose K so large that the coefficient n (int_{v >
K} v^p)^((p-2)/p) is at most d_gamma S_*/4. What is left of the potential,
on {v <= K}, is bounded by n K^(p-2) ||w_N||_2^2. Absorb and get

    ( d_gamma S_* / 4 ) ||w_N||_p^2  <=  C_K ||w_N||_2^2 .

With the assumed L^(2 gamma) moment the right side is bounded uniformly in N.
So ||w_N||_p is bounded, and (7.23) bounds E(w_N) too. On {v < N},
Dw_N = gamma v^(gamma-1) Dv; the functions w_N and Dw_N converge monotonically
and are dominated by v^gamma and gamma v^(gamma-1)|Dv|, so both converge in
L^2. Closedness of D gives v^gamma in H with the asserted gradient.

Start with v in L^p (from [P7]), set l_j = p (p/2)^j, and iterate the
conditional statement with gamma = l_j / 2, j = 0, 1, 2, ...: each step gains
v in L^(l_(j+1)), and l_j -> infinity since p/2 > 1. Finite measure then
gives every finite moment, and every q >= 0 has its weighted energy from the
gradient of v^(1 + q/2). End of sketch of [L13].

[L14] (compositions with polynomial growth). Let r in H with 0 <= r <= R <
infty almost everywhere and |Dr| in L^infty(mu). Let F be C^1 on a
neighbourhood of [0,R] x [0, infty) with

    |F(a,t)| + |partial_1 F(a,t)| + |partial_2 F(a,t)|  <=  C (1+t)^M

for some finite C, M. Then F(r,v) in H and

    D(F(r,v))  =  partial_1 F(r,v) Dr + partial_2 F(r,v) Dv .             (7.24)

Sketch. Let v_N = min(v,N); by [P5] it is in H with gradient
1_{v<N} Dv. Multiply F by a compactly supported cutoff equal to 1 on
[0,R] x [0,N]; the resulting bounded-slope function has, by [P5](c), the
expected chain rule for F(r,v_N). The growth hypothesis shows that
F(r,v_N), D(F(r,v_N)) are dominated by C' (1+v)^M (1+|Dv|), whose square is
integrable by [L13]. Dominated convergence and closedness give the result.
End of sketch of [L14].

7.7  Chord tests and the conformal transform.

The next several results use the following objects. Normalize

    nu  =  v^p dmu / W ,      a probability measure.                    (7.25)

For a centre y in Y with v(y) > 0, put

    kappa = v(y)^beta ,       r(x) = d(y,x) ,                           (7.26)

    s(x)  =  kappa * r(x) * v(x)^beta ,     L_y = int s^2 dnu ,         (7.27)

    J_y(a)  =  int (1 + a s^2)^(-n) dnu                      (a >= 0) . (7.28)

J_y is a transform of the distribution of squared chords; the problem is to
control it in terms of its second moment L_y.

[L15] (admissible chord tests). Fix a > 0, let w(t) = (1 + a t^2)^(-(n-2)/2),
and let s as in (7.27). Then:

  (a) v w(s), v w(s)^2 in H, with

          D(v w(s))     =  w(s) Dv + w'(s) xi ,
          D(v w(s)^2)   =  w(s)^2 Dv + 2 w(s) w'(s) xi ,                 (7.29)

      where xi = kappa ( v^(1+beta) Dr + beta r v^beta Dv ) is a
      square-integrable chart covector field;

  (b) for every smooth f: R -> [0, infty) with bounded value and bounded
      derivative, the functions

          g = kappa^2 v^(2 + 2 beta) f(s) ,
          psi = beta kappa^2 r^2 v^(1 + 2 beta) f(s)                     (7.30)

      are in H. In particular g is admissible in [P1] (for the current mu),
      and psi is admissible in (7.22).

Sketch. The only obstruction is that v^beta is not differentiable at v = 0
when beta < 1 (and beta = 1/(n-2) < 1 for n >= 4, so the obstruction is
real). The fix is to write the functions as F(r,v) and use [L14]. For
instance v w(s)^j = F_j(r,v) with t = kappa r v^beta and
F_j(r,z) = z w(t)^j. Computation:

    F_j(r,z) = z + O(z^(1+2 beta)) ,    partial_z F_j -> 1 ,
    partial_r F_j -> 0  as z -> 0 ,

so extension by F_j(r,z) = z on z <= 0 makes F_j jointly C^1. The functions
F_g and F_psi of (7.30) behave similarly, with their derivatives going to
zero at z = 0 since the smallest power of z in their derivative formulas is
z^(2 beta) = z^(2/(n-2)) > 0. So [L14] applies and gives the chain rules
(7.29). End of sketch of [L15].

[L16] (completed radial inequality). For every y in Y and every nonnegative
g in H,

    int ( n g |Dr|^2 + r <Dr, Dg> ) dmu  <=  (n/4) int r^2 g dmu .        (7.31)

Sketch. [P1] says int (n g + r <Dr, Dg>) dmu <= n int r g sqrt(1-|Dr|^2) dmu.
Consider the pointwise identity, with b = sqrt(1-|Dr|^2):

    1 - r b - |Dr|^2 + r^2/4  =  (b - r/2)^2  >= 0 .

Multiply by n g >= 0 and integrate:

    n int g - n int r g b  - n int g |Dr|^2  + (n/4) int r^2 g  >= 0 .

Substituting [P1] (n int g + int r <Dr, Dg>) <= n int r g b then gives
(7.31) at once. End of sketch of [L16].

[P8] (chord transform bound). For every centre y with v(y) > 0 and every
a >= 0,

    J_y(a)  <=  (1 + 2 a L_y)^(-n/2) .                                  (7.32)

Sketch. Fix y with v(y) > 0 and abbreviate kappa = v(y)^beta, r = d(y,.),
s = kappa r v^beta, xi = kappa (v^(1+beta) Dr + beta r v^beta Dv).

Step 1: a ground state identity. Test the Euler equation (7.22) with
psih = v w(s)^2, which is in H by [L15]. Using the product rules (7.29) and
noting that the first-order term in <D(v w(s)), D(v w(s))>,
2 w w' <Dv, xi>, equals the term in <Dv, D(v w(s)^2)>, one gets

    Q(v w(s))  =  n int v^p w(s)^2 dmu  +  4 beta int (w'(s))^2 |xi|^2 dmu . (7.33)

Step 2: quotient minimality. By the definition of v as a minimizer,

    Q(v w(s))  >=  (Q(v)/||v||_p^2) ||v w(s)||_p^2
              =  n W J(a)^((n-2)/n) ,                                (7.34)

since int v^p w(s)^p = int v^p (1 + a s^2)^(-n) dmu/W * W = W J(a), and
p (n-2)/2 = n.

Step 3: control the |xi|^2 term. Use [L16] plus the Euler equation divided by
4 beta with test psi from (7.30), the two adding to cancel the undifferentiated
terms and giving

    int ( n f(s) + s f'(s) ) |xi|^2 dmu
       <=  (n/4) int s^2 f(s) v^p dmu .                            (7.35)

Now choose f so that n f(t) + t f'(t) = (w'(t))^2, namely

    f(t)  =  int_0^1 tau^(n-1) (w'(tau t))^2 dtau ,

which is smooth, bounded, with bounded derivative, and satisfies the
identity by integration of d/d tau (tau^n (w'(tau t))^2). Combining (7.33),
(7.34) and (7.35):

    J(a)^((n-2)/n)
       <=  int ( w(s)^2 + beta s^2 f(s) ) dnu .                    (7.36)

Step 4: scalar reduction. With q = (n-2)/2 > 0 and b = a s^2,

    w(s)^2 + beta s^2 f(s)  =  q int_0^1 u^(q-1) / (1 + a s^2 u)^n du ,

by an exact derivative computation. In particular, the right side of (7.36)
is in (0,1]. So (7.36) becomes

    J(a)^((n-2)/n)  <=  G(a)  :=  q int_0^1 J(au) u^(q-1) du .          (7.37)

Step 5: solve the ODE. With h = G^(1/q) = G^(-2/(n-2)), (7.37) gives
ah' >= h - 1, hence the quotient (h(a) - 1)/a is nondecreasing. Its limit at
a -> 0+ is 2 L_y (differentiate under the integral using the finite second
moment), so h(a) >= 1 + 2 a L_y. Combining with (7.37):

    J(a)  <=  G(a)^(n/(n-2))  =  h(a)^(-n/2)  <=  (1 + 2 a L_y)^(-n/2) .

Which is (7.32). End of sketch of [P8].

7.8  The two opposite bounds on chord moments, and the contradiction.

[L17] (tangent density). For nu-almost every centre y,

    liminf_{a -> infty} a^(n/2) J_y(a)  >=  s_n / (2^n W) .             (7.38)

Consequently,

    L_y  <=  2 ( W / s_n )^(2/n)  <  2                          (7.39)

for nu-almost every y.

Sketch. Choose a chart and a point y = phi_i(z) where z is a density point of
E_i, the metric differential is the Euclidean norm |.|_{G_0}, and u_i = v o
phi_i, rho_i = |theta_i| sqrt(det G_i) are Lebesgue points with values
v_0 > 0, rho_0 > 0. Such points carry full nu-mass by [L2] and [L3].

Let a_j -> infinity realize the lower limit l of a^(n/2) J_y(a), and put
eps_j = a_j^(-1/2). After passing to a diagonal subsequence in R^n, the
rescaled functions u_i(z + eps_j h), rho_i(z + eps_j h), 1_{E_i}(z + eps_j h)
converge almost everywhere to v_0, rho_0, 1; and

    eps_j^(-1) d(y, phi_i(z + eps_j h))  ->  |h|_{G_0} .

Restricting the defining integral of J_y to the chart, changing variables,
and using Fatou (with a_i^(n/2) cancelling the Jacobian eps_j^(-n)):

    W l  >=  v_0^p rho_0 int_{R^n} (1 + v_0^(4 beta) |h|_{G_0}^2)^(-n) dh .

The linear substitution H = v_0^(2 beta) G_0^(1/2) h has Jacobian
v_0^(2n beta) sqrt(det G_0) = v_0^(p - 2) sqrt(det G_0). The integrand under
this substitution is exactly (1 + |H|^2)^(-n), and using the stereographic
integral (7.14):

    W l  >=  |theta_i(z)| s_n / 2^n  >=  s_n / 2^n ,

where the last step uses |theta_i(z)| >= 1 (integrality). This is (7.38). To
deduce (7.39), combine (7.32) with (7.38), solve for L_y:

    s_n / (2^n W)  <=  (2 L_y)^(-n/2)  =>  L_y  <=  2 (W/s_n)^(2/n) ,

and 2(W/s_n)^(2/n) < 2 because W < s_n. The case L_y = 0 is immediate.
End of sketch of [L17].

[L18] (opposite averaged bound). Define

    B_* = int v^(p + 2 beta) dmu ,     I(y) = int d(y,x)^2 v(x)^(p + 2 beta) dmu(x) .

Then 0 < B_* < infinity, I is finite and continuous, and

    W^2  <=  B_* I(y)  for every y in Y .                           (7.40)

Consequently,

    int L_y dnu  >=  2 .                                                 (7.41)

Sketch. (7.40). Fix y, put r(x) = d(y,x). By [L13], g = v^p is in H with
Dg = p v^(p-1) Dv, so [P1] applies to g. After dividing by n and using
p/n = 2 beta, [P1] reads

    W  <=  int r ( v^p sqrt(1-|Dr|^2) - 2 beta v^(p-1) <Dr, Dv> ) dmu .  (7.42)

Pointwise, in the direct sum of the chart cotangent space with R, the vector
(sqrt(1-|Dr|^2), Dr) has norm one. So the bracket is the inner product of
that unit vector with (v^(1+beta), -2 beta v^beta Dv), multiplied by
v^(p/2 + beta), which is at most v^(p/2+beta) ( v^(2+2 beta) + 4 beta^2
v^(2 beta) |Dv|^2 )^(1/2). Applying Cauchy-Schwarz to (7.42):

    W^2  <=  ( int d(y,x)^2 v^(p+2 beta) dmu(x) ) *
             ( int ( v^(2+2 beta) + 4 beta^2 v^(2 beta) |Dv|^2 ) dmu ) .

The first factor is I(y). The second is computed from (7.22) with
psi = v^(1+2 beta), whose gradient is (1+2 beta) v^(2 beta) Dv. Then

    4 beta (1+2 beta) int v^(2 beta)|Dv|^2 + n int v^(2+2 beta)
       =  n int v^(p+2 beta) ,

i.e. the second factor equals B_* by the identity 1+2 beta = n beta. This
gives (7.40).

(7.41). Use the variance inequality [lit 12]: the function I has a minimizer
o (a barycentre) and satisfies I(y) >= I(o) + B_* d(o,y)^2 for every y.
Now

    int L_y dnu  =  (1/W^2) int I(y) v^(p+2 beta) dmu(y)
                >=  (1/W^2) ( B_* I(o) + B_* int d(o,y)^2 v^(p+2beta) dmu )
                =  (1/W^2) ( B_* I(o) + B_*^2 ... cancel that ... )

Hmm, careful: the second term of the inequality is
B_* int d(o,y)^2 v^(p+2beta) dmu(y) = B_* I(o)? Not exactly; I(o) is the
integral of d(o,x)^2 v(x)^(p+2beta) dmu(x), which is the same as
int d(o,y)^2 v(y)^(p+2beta) dmu(y). Yes, both are I(o), by renaming. So

    int L_y dnu  >=  ( 1/W^2 ) ( B_* I(o) + B_* I(o) )  =  2 B_* I(o)/W^2 .

By (7.40) at y = o, W^2 <= B_* I(o), and the right side is >= 2. End of
sketch of [L18].

[T3] (induction step). Let n >= 3, let Y be a compact CAT(0) space, and
assume (H_(n-1)). Then every integral n-cycle T in Y has an integral filling
S with dS = T and M(S) <= c_n M(T)^((n+1)/n).

Proof. Suppose not. By [L7] and [L8] there is a champion A with
m = M(A) < s_n and normalization (6.1). Build v from [P7] and nu from
(7.25).

Bound above: by (7.32) and (7.39), J_y(a) <= (1+2aL_y)^{-n/2}, so

    int (1 + a s_y^2)^(-n) dnu  <=  (1 + 2 a L_y)^(-n/2)
                                 <=  (1 + 2 a * 2(W/s_n)^(2/n))^(-n/2)
                                 =  (1 + 4 a (W/s_n)^(2/n))^(-n/2) ,

which sends a -> infinity. That is not the useful direction; the useful one
is:

Bound on the integral itself. Integrate (7.39) over nu: by the definition of
L_y,

    int L_y dnu  =  int int (v(y)^beta d(y,x) v(x)^beta)^2 dnu(x) dnu(y)
                 <=  2 (W/s_n)^(2/n) .                              (7.43)

Bound on the integral from below. (7.41) says int L_y dnu >= 2. Combine with
(7.43) and use W < s_n:

    2  <=  int L_y dnu  <=  2 (W/s_n)^(2/n)  <  2 ,

a contradiction. So no champion exists, and the sharp filling inequality
holds in dimension n. End of proof of [T3].