Functional calculus and harmonic analysis in geometry

TL;DR

This survey highlights recent advances in solving geometric problems by integrating functional calculus and harmonic analysis, aiming to inspire further research in these interconnected mathematical fields.

Contribution

It provides an overview of how functional calculus and harmonic analysis are combined to address complex geometric problems, encouraging cross-disciplinary development.

Findings

Resolution of specific geometric problems using harmonic analysis techniques

Demonstration of the interplay between functional calculus and geometry

Potential for broader applications in geometric analysis

Abstract

In this short survey article, we showcase a number of non-trivial geometric problems that have recently been resolved by marrying methods from functional calculus and real-variable harmonic analysis. We give a brief description of these methods as well as their interplay. This is a succinct survey that hopes to inspire geometers and analysts alike to study these methods so that they can be further developed to be potentially applied to a broader range of questions.

Equations

$$\mathrm{dom}(\sqrt{-\mathrm{div}B\nabla })={\mathrm{H}}^1({\mathbb{R}}^n)=\left\{u\in {\mathrm{L}}^2:\nabla u\in {\mathrm{L}}^2\right\}$$

$$\mathrm{Re}{⟨B(x)u,u⟩}_{{\mathbb{R}}^n}\ge \kappa |u{|}^2$$

$${\Pi }_B=\left(\begin{array}{cc}0& -\mathrm{div}B\\ \nabla & 0\end{array}\right)$$

$$\zeta \mapsto \frac{\sqrt{{\zeta }^2}}{\zeta }\text{and}\zeta \mapsto \frac{\zeta }{\sqrt{{\zeta }^2}}$$

$$u(\theta )=\sum _{n=-\infty }^{\infty }{a}_n{\mathrm{e}}^{ın\theta },$$

$$({\Delta }_{{\mathrm{S}}^1}u)(\theta )=\sum _{n=-\infty }^{\infty }{n}^2{a}_n{\mathrm{e}}^{ın\theta }.$$

$$(f({\Delta }_{{\mathrm{S}}^1})u)(\theta ):=\sum _{n=-\infty }^{\infty }f(n^2)a_n{\mathrm{e}}^{ın\theta }.$$

$$fg\mapsto f({\Delta }_{{\mathrm{S}}^1})g({\Delta }_{{\mathrm{S}}^1}),(f\circ g)\mapsto f(g({\Delta }_{{\mathrm{S}}^1})),\text{and}1\mapsto \mathrm{I}.$$

$${\partial }_{t}u(t,\theta )+{\Delta }_{{\mathrm{S}}^1}u(t,\theta )=0,\underset{t\to 0}{\lim }u(t,\theta )=u_0(\theta ).$$

$$f(\mathrm{D})u:=\sum _{n}f({\lambda }_n)a_n{\psi }_n.$$

$$|\psi (x)|\le C\min \left\{|x{|}^{\alpha },|x{|}^{-\alpha }\right\}$$

$$\Vert {\chi }_{[1,2]}(t\mathrm{D})u{\Vert }^2$$

$$=\sum _n{a}_n^2\Vert {\chi }_{[1,2]}(t{\lambda }_n){\psi }_n{\Vert }^2=\sum _n{a}_n^2|{\chi }_{[1,2]}(t{\lambda }_n){|}^2\Vert {\psi }_n{\Vert }^2.$$

$${\int }_0^{\infty }\Vert {\chi }_{[1,2]}(t\mathrm{D}){\psi }_n{\Vert }^2\frac{dt}{t}=\Vert {\psi }_n{\Vert }^2{\int }_0^{\infty }|{\chi }_{[1,2]}(t{\lambda }_n){|}^2\frac{dt}{t}.$$

$${\int }_0^{\infty }|{\chi }_{[1,2]}(t{\lambda }_n){|}^2\frac{dt}{t}$$

$$={\int }_0^{\infty }|{\chi }_{[1,2](s)}{|}^2\frac{ds}{s}={\int }_1^2\frac{1}{s}ds=\log (2).$$

$${\int }_0^{\infty }\Vert {\chi }_{[1,2]}(t\mathrm{D})u{\Vert }^2\frac{dt}{t}=\sum _n{a}_n\log (2)\Vert u_n{\Vert }^2=\log (2)\Vert u{\Vert }^2.$$

$${\int }_0^{\infty }\Vert \psi (t\mathrm{D})u_n{\Vert }^2\frac{dt}{t}=\Vert u_n{\Vert }^2{\int }_0^{\infty }|\psi (t{\lambda }_n){|}^2\frac{dt}{t}\simeq \Vert u_n{\Vert }^2$$

$${\int }_0^{\infty }\Vert \psi (t\mathrm{D})u{\Vert }^2\frac{dt}{t}\simeq \Vert u{\Vert }^2.$$

$${\mathrm{S}}_{\alpha }^{\mathrm{o}}:=\left\{\zeta \in \mathbb{C}\setminus \left\{0\right\}:|\arg (\zeta )|<\alpha \text{or}|\arg (-\zeta )|<\alpha \right\},$$

$$|\psi (\zeta )|\le C\max \left\{|\zeta {|}^{\alpha },|\zeta {|}^{-\alpha }\right\}$$

$$\zeta \mapsto \frac{\zeta }{1+{\zeta }^2}$$

$$\psi (T)u:=\frac{1}{2\pi ı}{\oint }_{\gamma }\psi (\zeta )(\zeta -T{)}^{-1}ud\zeta ,$$

\gamma=\left\{te^{\imath\nu}:\ \infty>t>0\right\}+\left\{-te^{-\imath\nu}:\ 0<t<\infty\right\}\\ +\left\{te^{\imath\nu}:\ \infty>t>0\right\}+\left\{te^{-\imath\nu}:\ 0<t<\infty\right\}

$$\Vert \psi (T)u\Vert \le C_{\psi ,T}\Vert u\Vert .$$

$$\mathcal{H}=\ker (T)\oplus \overline{\mathrm{ran}(T)}.$$

$$\Vert \psi (T)u\Vert \le C\Vert \psi {\Vert }_{\infty }\Vert u\Vert$$

$${\chi }^{\pm }(\zeta ):=\{\begin{array}{ll}1& \pm \mathrm{Re}\zeta >0,\\ 0& \text{otherwise},\end{array}$$

$$f(T)u:=f(0)P_{\ker (T),\overline{\mathrm{ran}(T)}}u+\underset{n\to \infty }{\lim }{f}_n(T)u.$$

$$\Vert f(T)\Vert \le C\Vert f{\Vert }_{\infty }.$$

$${\mathrm{L}}_{A,a}u:=-a\mathrm{div}A\nabla ,$$

Full text

1. 1 Introduction
2. 2 Functional Calculus
3. 3 Geometry
4. 4 Harmonic Analysis

1. Introduction

The purpose of this short survey article is to highlight some geometric questions that have recently been addressed via the use of functional calculus and harmonic analytic methods. In the literature, these phrases are umbrella terms that carry different meanings in different schools. Therefore, before we embark on considering the details, let us first note what we mean by these terms.

Functional calculus is the ability to take functions of operators and manipulate them as if they were functions. This is an extraordinarily powerful way to think and one motivation is that it provides a conceptual method to solve partial differential equations, particularly evolution equations. Harmonic analysis, on the other hand, is the art in which a mathematical object can be seen as a “signal”, and whose goal is to decompose this signal, typically in some scale invariant way, to simpler parts which are mathematically tractable. Perhaps the simplest of this kind is the Fourier series, where every square integrable function on the circle can be written as an infinite sum of trigonometric functions, where the properties of these latter functions are well known and understood. Geometry, to us, almost always means the shape of space. In practice, this is quantified in a multitude of ways, but here, we simply contrast it from the topology or structure of the space, which is seldom of concern to us.

At a first glance, perhaps it might seem peculiar to talk about a combination of functional calculus and harmonic analysis. However, this is an extremely fruitful marriage, and the late Alan McIntosh was one of the pioneers of this school of thought. His idea was to demonstrate the ability to construct a powerful and desirable functional calculus, namely the $\mathrm{H}^{\infty}$ functional calculus, via quadratic estimates. These estimates can be interpreted as the ability to reconstruct a signal in norm via operator adapted band-pass filters, where these filters only rely on the spectral properties of the operator. Taking a perspective of Fourier theory as a functional calculus for the Laplacian, the $\mathrm{H}^{\infty}$ functional calculus can be seen as a Fourier theory adapted to an operator. The quadratic estimates point of view of this functional calculus is extremely useful as it provides a bridge to methods emerging from real-variable harmonic analysis in order to compute these estimates.

An important motivating factor for the development of these methods was the treatment of the Kato square root problem, which remained open for almost half a century. In [27], Kato formulated the most abstract version of the problem. He first demonstrated that if $A$ is a maximally-accretive operator, meaning that the spectrum and numerical range sit in the complex plane with positive real part, then for any $\alpha<\frac{1}{2}$, the domains ${\mathrm{dom}}(A^{\alpha})={\mathrm{dom}}({{A}^{\ast}}^{\alpha})$. The equality was known to be invalid in general when $\alpha>\frac{1}{2}$ and Lions in [29] showed that that this equality even fails in general for maximally-accretive operators for the critical case $\alpha=\frac{1}{2}$. However, the question remained open as to whether the critical case was true for regularly-accretive operators, which are operators whose numerical range and spectrum are contained in a proper sector of the complex plane with positive real part. In 1972, McIntosh authored [31] which dealt a fatal blow to the abstract question by constructing a regularly-accretive operator $A$ for which ${\mathrm{dom}}(\sqrt{A})\neq{\mathrm{dom}}(\sqrt{{A}^{\ast}})$. Despite this fact, given that the main interest of Kato’s original question was in its applications to partial differential equations consisting of operators in divergence form, the Kato square root problem as we know it today became rephrased as to determine whether

[TABLE]

for coefficients $B\in{\rm L}^{\infty}(\mathbb{R}^{n})$ with a constant $\kappa>0$ such that

[TABLE]

almost-everywhere in $\mathbb{R}^{n}$. If the coefficients $B(x)$ are real and symmetric, then the desired equality follows simply from operator theoretic considerations alone. It is the complex, non-symmetric coefficient case that is of significance.

The resolution of this question for $n=1$ was the famous paper [21] due to Coifman, McIntosh and Meyer which they authored in 1981. It is a paper that has spawned a number of new directions in mathematics including wavelet theory. Moreover, this paper made significant contributions to the theory of singular integrals as well as the development of the so-called $T(b)$ theorems. Beyond their general usefulness, this was of paramount importance to the resolution of the Kato square root problem in arbitrary dimensions.

In one dimension, the Kato square root problem is actually equivalent to the boundedness of the Cauchy integral operator on a Lipschitz curve. However, it is the description of this problem from a differential operator point of view, where the functional calculus perspective takes centre stage, that was a key observation to its resolution. The general case was finally settled in 2002 by Auscher, Hofmann, Lacey, McIntosh and Tchamitchian in [8]. There are numerous surveys of the Kato square root problem and its resolution. Namely, the lecture series [7] by Auscher, as well as the survey article [26] by Hofmann and McIntosh give detailed accounts. See also references therein.

The general Kato square root problem was initially proved using second-order methods in [8]. In 2005, in [10], Axelsson (Rosén), Keith and McIntosh produced a first-order perspective of this problem. There, they constructed a larger first-order Dirac-type operator which encoded the problem. The idea was to rewrite the divergence form equation into the system

[TABLE]

and by showing that this operator has an $\mathrm{H}^{\infty}$ functional calculus, the desired conclusion was obtained. The required domain equality follows from showing that the functions

[TABLE]

applied to $\Pi_{B}$ define a bounded operator. In [10], the authors also demonstrate how to prove a similar problem for perturbations of the Hodge-Dirac operator of the form $D_{B}={\rm d}+B^{-1}{{\rm d}}^{\ast}B$. They further obtain corresponding results on compact manifolds. In fact, it this paper which provided the key conceptual point of view to make similar methods work in more geometric settings. It is important to also acknowledge Morris, where in his thesis [33], he developed many necessary technical results on metric spaces which were later adapted to attack geometric problems on vector bundles over manifolds.

Despite the fact that these methods were developed primarily for non-smooth Euclidean problems, their success in analysing differential operators on vector bundles has been for reasons similar in spirit to their success in the non-smooth setting. Pseudo-differential or Fourier analytic methods often rely on localisation and smoothness. The functional calculus that we will talk about here is sufficiently implicit and general to be able to treat non-smooth coefficients, which requires a global point of view. In a similar vein, bundles typically carry non-trivial geometry, and even if the coefficients of the problem are smooth, the tools need to be sufficiently abstract in order to capture the global nature without resorting to local descriptions. This is really the key indication that the methods borne out of the resolution of the Kato square root problem, and in particular its first-order reincarnation, might have useful applications to geometry, in particular to analysis on bundles.

Lastly, we remark that this survey is far from being comprehensive. The results we discuss here are almost exclusively contributions of the author and co-authors. It represents only a small subset of the successes of coupling functional calculus and harmonic analysis to be an effective force in solving problems in geometry. The author does not have sufficient expertise to comment on the general picture, and so the limited focus here is for the sake of conciseness and accuracy. Moreover, the author apologises in advance for the many omitted names and references. It is the author’s hope that geometers and analysts become better aware of these methods and embark on further developing these tools. In particular, there is a real need to adapt the methods arising from real-variable harmonic analysis, which are typically Euclidean in nature, to better account for the underlying geometry.

Acknowledgements

The author was supported by SPP2026 from the German Research Foundation (DFG). David Rule deserves mention as this paper was borne out of a Colloquium talk given at Linköping University upon his invitation. Andreas Rosén needs to acknowledged for useful feedback on an earlier version of this article. The anonymous referee also deserves a mention for their helpful suggestions. Furthermore, the author acknowledges the gracious hospitality of his auntie Harshya Perera in Sri Lanka who hosted him during the time in which this paper was written.

2. Functional Calculus

2.1. Fourier series

Let us begin our exposition by stating a beautiful and well known fact seen from a seldom mentioned perspective.

Recall that for a function $u\in{\rm L}^{2}(\mathrm{S}^{1})$, we can write its Fourier series as

[TABLE]

where the coefficients $a_{n}\in\mathbb{C}$ are uniquely determined by $u$. Now, consider the Laplacian $\Delta_{\mathrm{S}^{1}}$ on ${\rm L}^{2}(\mathrm{S}^{1})$, or equivalently, $-\frac{d^{2}}{d\theta^{2}}$ on $[0,1]$ with periodic boundary conditions on ${\rm L}^{2}([0,1])$. It is well known that the spectrum of this operator consists only of eigenvalues and it is precisely $\mathrm{spec}(\Delta_{\mathrm{S}^{1}})=\left\{\lambda_{n}=n^{2}\right\}_{n=0}^{\infty}$. The corresponding eigenfunctions are $\left\{\mathrm{e}^{\imath n\theta},\mathrm{e}^{-\imath n\theta}\right\}_{n=0}^{\infty}$.

When $u\in{\mathrm{dom}}(\Delta_{\mathrm{S}^{1}})$, the operator $\Delta_{\mathrm{S}^{1}}$ acts on $u$ as

[TABLE]

Therefore, we see the decomposition of $u$ as a Fourier series is, in fact, the decomposition of $f$ with respect to the frequencies of the Laplacian on $\mathrm{S}^{1}$.

Via (2.1), we can construct a rudimentary functional calculus for $\Delta_{\mathrm{S}^{1}}$. We call a function $f:\mathbb{Z}\to\mathbb{C}$ a symbol and define a function of $\Delta_{\mathrm{S}^{1}}$ by $f$ via:

[TABLE]

Here, the significance of spectral theory in the construction of this functional calculus is explicit. Indeed, this will be true for functional calculi we consider for more general operators. For an arbitrary function $f$, we would need to compute the domain of this operator, which is precisely the $u\in{\rm L}^{2}(\mathrm{S}^{1})$ for which the defining series converges. Moreover, it is easy to see that if $f$ is bounded, then $f(\Delta_{\mathrm{S}^{1}})$ is a bounded operator on ${\rm L}^{2}(\mathrm{S}^{1})$.

Note that this functional calculus, at least for bounded functions, enjoys a few desirable properties. Of significance is that it is an algebra homomorphism $f\mapsto f(\Delta_{\mathrm{S}^{1}})$, which satisfies

[TABLE]

These are essential features of the functional calculi we consider in this survey.

An important motivation for functional calculi comes from evolution equations, where solutions can be simply generated and represented through the functional calculus. For example, let us consider the following heat equation:

[TABLE]

The unique solution to this is given by $u(t,\theta)=(\mathrm{e}^{-t\Delta_{\mathrm{S}^{1}}}u_{0})(\theta)=\sum_{n=-\infty}^{\infty}e^{-tn^{2}}a_{n}e^{\imath n\theta}$.

2.2. Non-negative self-adjoint operators with discrete spectrum

Taking inspiration from this very concrete and classical example, let us instead consider the following abstract setup. Let $\mathscr{H}$ be a Hilbert space, ${\rm D}$ a non-negative self-adjoint operator on $\mathscr{H}$ with discrete spectrum $\mathrm{spec}({\rm D})=\left\{\lambda_{i}\geq 0\right\}$. We denote its eigenfunctions by $\left\{\psi_{i}\right\}$, which are orthogonal due to self-adjointness and as before, whenever $u\in\mathscr{H}$, we obtain coefficients $a_{n}\in\mathbb{C}$ such that $u=\sum_{n}a_{n}\psi_{n}$. In this context, a functional calculus for ${\rm D}$ can be generalised mutatis mutandis from the setting of $\mathrm{S}^{1}$ in the following way:

[TABLE]

As before, the $f$ needs to be chosen from an appropriate class in order to obtain a well-defined expression and to understand the mapping properties of $f({\rm D})$. For instance, it is easy to see that if $f:\mathbb{R}\to\mathbb{C}$ is bounded and continuous, the operator $f({\rm D})$ is well-defined and bounded on $\mathscr{H}$.

There are good reasons for considering this more abstract setup as the following example illustrates. Let $(\mathcal{M},\mathrm{g})$ be a smooth compact Riemannian manifold (without boundary), $\mathscr{H}={\rm L}^{2}(\mathcal{M})={\rm L}^{2}(\mathcal{M},\mathrm{g})$, and ${\rm D}=\Delta_{\mathcal{M},\mathrm{g}}=\nabla^{\ast,g}\nabla$, the Laplacian on $(\mathcal{M},\mathrm{g})$. Then, ${\rm id}:{\rm H}^{\rm 1}(\mathcal{M})\hookrightarrow{\rm L}^{2}(\mathcal{M})$ is a compact map and ${\mathrm{dom}}(\Delta_{\mathrm{g}})={\rm H}^{\rm 2}(\mathcal{M})\subset{\rm H}^{\rm 1}(\mathcal{M})$. The combination of these two facts yield that $\mathrm{spec}(\Delta_{\mathcal{M}})$ is discrete with no finite accumulation points. This exactly fits the abstract picture we have just painted. Of particular significance is the fact that we have not alluded to the exact value of the eigenvalues and their corresponding eigenfunctions, since unlike $\Delta_{\mathrm{S}^{1}}$ on $\mathrm{S}^{1}$, in general, they cannot be computed.

It is a maxim that the better you know the spectrum, the larger the class of functions that you can consider in order to construct a functional calculus. Since ${\rm D}$ is a self-adjoint operator, it further admits a functional calculus for Borel functions $f:\mathbb{R}\to\mathbb{C}$, which is obtained by extracting a spectral measure. We shall not indulge in the details of its construction but it can be found in [28] by Kato.

We now touch upon an important property of these functional calculi - the reconstruction of a signal in norm. Let us consider $u\in\mathscr{H}$ as a “signal”. We want to use the functional calculus of ${\rm D}$ to reconstruct $\|u\|$ but up to a constant. To keep our exposition simple, let us assume ${\mathrm{ker}}({\rm D})=0$. This is not a severe restriction: if ${\mathrm{ker}}({\rm D})\neq 0$, then we can split $\mathscr{H}={\mathrm{ker}}({\rm D})\oplus^{\perp}\overline{{\mathrm{ran}}({\rm D})}$ and instead work on $\overline{{\mathrm{ran}}({\rm D})}$ as the new space. Let $\varphi:\mathbb{R}\to\mathbb{R}$ be a piecewise smooth function such that $\psi\neq 0$ and for which there exists $\alpha>0$ and $C>0$ satisfying:

[TABLE]

for almost-every $x\in\mathbb{R}$. We want to justify the idea that $\psi(t{\rm D})$ is a band-pass filter for frequencies localised around $\frac{1}{t}$. Let us first compute with $\psi(x)=\chi_{[1,2]}(x)$ in order to make the exposition clear.

On writing $u=\sum_{n}a_{n}\psi_{n}$, note that

[TABLE]

Therefore,

[TABLE]

On substituting $s=t\lambda_{n}$ so that $dt=\frac{ds}{\lambda_{n}}$ and $\frac{1}{t}=\frac{\lambda_{n}}{s}$, we obtain

[TABLE]

Combining (2.3) with (2.4), we obtain

[TABLE]

If we were to repeat this exercise with a general $\psi$ satisfying (2.2), then by a similar calculation, we would obtain that

[TABLE]

so that

[TABLE]

We emphasise here that in (2.4), the cancellation of the spectral points $\lambda_{n}$ under a change of variables occurs precisely due to the fact that we compute with respect to the Haar measure $\frac{dt}{t}$ on $\mathbb{R}$. It is instructive to repeat this simple calculation with respect to another measure to better understand the significance of $\frac{dt}{t}$.

2.3. Bisectorial operators and the McIntosh Theorem

The key idea in the calculation in the previous section is that $\psi(t{\rm D})$ filters out all frequencies but those concentrated near $\frac{1}{t}$, and at least for the class of invertible self-adjoint operators with discrete spectrum that we considered, and despite the lack of ability to explicitly compute the spectrum, we were able to reconstruct the signal in norm up to a constant. We will see here that, in fact, this is a characterising condition for the existence of a more significant and useful functional calculus that can be applied to a broader class of operators.

Before we introduce this class of operators, for $\alpha<\pi/2$, let us define

[TABLE]

the open bisector of angle $\alpha$ in the complex plane, and $\mathrm{S}_{\alpha}:=\overline{\mathrm{S}^{\mathrm{o}}_{\alpha}}$, the closed bisector of angle $\alpha$. For $\omega<\pi/2$, an operator $T:{\mathrm{dom}}(T)\subset\mathscr{H}\to\mathscr{H}$ is called $\omega$-bisectorial if:

1. (i)

$T$ is closed, 2. (ii)

$\mathrm{spec}(T)\subset\mathrm{S}_{\omega}$, and 3. (iii)

for all $\mu\in(\omega,\pi/2)$, there exists $C_{\mu}>0$ such that the resolvent estimate $|\zeta|\|(\zeta-T)^{-1}\|\leq C_{\mu}$ holds for all $\zeta\in\mathbb{C}\setminus\mathrm{S}^{\mathrm{o}}_{\mu}$.

Operators that are $\omega$-bisectorial for which the half plane with non-positive real part $\left\{\zeta\in\mathbb{C}:\mathrm{Re}\ \zeta<0\right\}$ is free of spectrum are called $\omega$-sectorial. Every self-adjoint operator is [math]-bisectorial and non-negative self-adjoint operators are [math]-sectorial.

We remark that sectorial operators arise naturally in the study of holomorphic semigroup theory. In fact, in [28], there is a notion of sectoriality which further includes the condition on the numerical range, namely that $\mathrm{Re}\ \left\langle Tu,u\right\rangle\geq 0$. In modern language, these are called Kato sectorial operators.

Bisectorial operators automatically admit a rudimentary functional calculus. Namely, fixing $\mu\in(\omega,\pi/2)$, by $\Psi(\mathrm{S}^{\mathrm{o}}_{\mu})$, let us denote the set of holomorphic $\psi:\mathrm{S}^{\mathrm{o}}_{\mu}\to\mathbb{C}$ satisfying the estimate: there exists $\alpha>0$ and $C>0$ such that

[TABLE]

for all $\zeta\in\mathrm{S}^{\mathrm{o}}_{\mu}$. Note that this is reminiscent of the condition (2.2), and in light of this, functions in $\Psi(\mathrm{S}^{\mathrm{o}}_{\mu})$ can be thought of as complex valued functions that can be used to define operator adapted band-pass filters. For instance,

[TABLE]

is an example of a $\Psi(\mathrm{S}^{\mathrm{o}}_{\mu})$ function.

Now, a rudimentary functional calculus for $T$ can be defined for each $\psi\in\Psi(\mathrm{S}^{\mathrm{o}}_{\mu})$ via the expression

[TABLE]

where $\gamma$ is the unbounded contour which runs through zero and infinity as

[TABLE]

with $\nu\in(\omega,\mu)$. The reason we say it is rudimentary is because via this definition, it is not possible to construct $1(T)$. However, for $\psi\in\Psi(\mathrm{S}^{\mathrm{o}}_{\mu})$, the integral in (2.5) converges absolutely in the strong operator topology and defines a bounded operator satisfying: there exists a constant $C_{\psi,T}$ (possibly dependent on $\psi$ and $T$) such that

[TABLE]

Moreover, the bisectoriality of $T$ means that the kernel of $T$ is is accessible from the imaginary axis which, excluding the origin, is contained in the resolvent set with resolvent estimates iii. This induces a decomposition of $\mathscr{H}$ as:

[TABLE]

Note that in general, this is a topological sum and not an orthogonal one. As aforementioned, self-adjoint $T$ are [math]-bisectorial and in this case, this sum is indeed orthogonal.

Despite the ability to create such a functional calculus with many salient features, in applications, it is often necessary to control the norm $C_{\psi,T}$. If there exists a universal constant $C>0$ such that $C_{\psi,T}\leq C\|\psi\|_{\infty}$, i.e.,

[TABLE]

for all $\psi\in\Psi(\mathrm{S}^{\mathrm{o}}_{\mu})$, then we say that $T$ has an $\mathrm{H}^{\infty}$ functional calculus. This nomenclature will soon be justified. Let us define a larger class of functions $\mathrm{H}^{\infty}(\mathrm{S}^{\mathrm{o}}_{\mu})$, where we say that $f\in\mathrm{H}^{\infty}(\mathrm{S}^{\mathrm{o}}_{\mu})$ if:

1. (i)

$f:\mathrm{S}^{\mathrm{o}}_{\mu}\cup\left\{0\right\}\to\mathbb{C}$ is bounded, 2. (ii)

$f{{\lvert}}_{\mathrm{S}^{\mathrm{o}}_{\mu}}:\mathrm{S}^{\mathrm{o}}_{\mu}\to\mathbb{C}$ is holomorphic, and 3. (iii)

$\|f\|_{\infty}<\infty$.

Note that $f$ is allowed to jump across zero. In particular, $\operatorname{sgn}(\zeta):=\chi^{+}(\zeta)-\chi^{-}(\zeta)$, where

[TABLE]

is a $\mathrm{H}^{\infty}(\mathrm{S}^{\mathrm{o}}_{\mu})$ class function.

Before we extend the rudimentary functional calculus to $\mathrm{H}^{\infty}(\mathrm{S}^{\mathrm{o}}_{\mu})$ functions, let us note the following lemma due to McIntosh.

Lemma 2.1** (McIntosh convergence lemma).**

If $T$ has an $\mathrm{H}^{\infty}$ functional calculus, i.e. satisfies (2.7), then for each $f\in\mathrm{H}^{\infty}(\mathrm{S}^{\mathrm{o}}_{\mu})$, there exists $f_{n}\in\Psi(\mathrm{S}^{\mathrm{o}}_{\mu})$ such that $f_{n}\to f$ uniformly on compact subsets of $\mathrm{S}^{\mathrm{o}}_{\mu}$ and for each $u\in\mathscr{H}$, the sequence $\left\{f_{n}(T)u\right\}$ is Cauchy.

This lemma first appeared in [32] (page 218) by McIntosh. Armed with this, and writing $\operatorname{\mathbf{P}}_{{\mathrm{ker}}(T),\overline{{\mathrm{ran}}(T)}}:\mathscr{H}\to{\mathrm{ker}}(T)$ for the projector with kernel $\overline{{\mathrm{ran}}(T)}$ via the splitting (2.6), we define

[TABLE]

This defines a bounded operator $f(T)$ with norm

[TABLE]

It is readily checked that this limit is independent of the particular sequence $\left\{f_{n}(T)u\right\}$ as guaranteed by Lemma 2.1. Moreover, this justifies the language “$\mathrm{H}^{\infty}$ functional calculus”. In this setting, by definition, we are using the space $\Psi(\mathrm{S}^{\mathrm{o}}_{\mu})$ as a dense subset of $\mathrm{H}^{\infty}(\mathrm{S}^{\mathrm{o}}_{\mu})$. In a private communication, Alan McIntosh alluded to the construction of the $\mathrm{H}^{\infty}$ functional calculus as reminiscent to the construction of the Fourier transform in ${\rm L}^{2}(\mathbb{R}^{n})$, where the $\Psi(\mathrm{S}^{\mathrm{o}}_{\mu})$ can be thought of as “complex valued Schwartz class functions”.

Bisectorial operators and their accompanying functional calculi are very naturally motivated. A quintessential example is the first-order factorisation of the Kato square root problem. Let $x\mapsto A(x)\in{\rm L}^{\infty}(\mathrm{SymMat}(\mathbb{R}^{n}))$, real self-adjoint at almost-every $x$ and $a:\mathbb{R}\to[0,\infty]$ a measurable function.

We assume that both $A$ and $a$ are uniformly bounded above and below. Define the operator:

[TABLE]

which is seen to be densely-defined and closed via divergence-form methods. Of interest is the following nonlinear perturbation estimate:

[TABLE]

Note that here, only the ${\rm L}^{\infty}$ difference of the coefficients appear in the bound - there are no first-order differences.

In order to analyse this operator, the authors of [10] consider the operator

[TABLE]

where the matrix-valued map $x\mapsto B(x)$ and the function $x\mapsto b(x)$ are both complex valued and bounded. Although a moot point, we remark that $B(x)$ need not be symmetric. We also assume that there exists $\kappa>0$ such that $B$ and $b$ satisfy the following accretivity condition:

[TABLE]

The operator $\Pi_{B,b}$ is an $\omega$-bisectorial operator on $\mathscr{H}={\rm L}^{2}(\mathbb{R}^{n})\oplus^{\perp}{\rm L}^{2}(\mathbb{R}^{2n})$ (since $\nabla:{\rm H}^{\rm 1}(\mathbb{R}^{n})\subset{\rm L}^{2}(\mathbb{R}^{n})\to{\rm L}^{2}(\mathbb{R}^{2n})$) with $\omega<\pi/2$ dependent only on $\kappa$ and the ${\rm L}^{\infty}$ bound on $B$.

Then, note that its square is precisely

[TABLE]

and if $b=a$ and $B=A$, the operator (2.9) is accessed by projecting to the first entry.

Now, let us assume that $\Pi_{B,b}$ has an $\mathrm{H}^{\infty}$ functional calculus. Then, as aforementioned, $f(\Pi_{B,b})$ defined for $f\in\mathrm{H}^{\infty}(\mathrm{S}^{\mathrm{o}}_{\mu})$. Moreover, we have that $(B^{\prime},b^{\prime})\mapsto f(\Pi_{B^{\prime},b^{\prime}})$ is holomorphic in a small ball around $(B,b)$ in the ${\rm L}^{\infty}$ topology whose radius is determined by the $\kappa$ and the ${\rm L}^{\infty}$ bounds on $B$ and $b$. Consequently, we obtain the Lipschitz estimate

[TABLE]

for $(B^{\prime},b^{\prime})$ in this ball. In particular, on setting $f_{t}(\zeta)=\mathrm{e}^{-t\zeta^{2}}$ for $t>0$, we obtain

[TABLE]

Then, for $u\in{\rm L}^{2}(\mathbb{R}^{n})$,

[TABLE]

This is precisely the desired estimate in (2.10).

Although it is possible to obtain such deep consequences in a conceptual manner via the $\mathrm{H}^{\infty}$ functional calculus, proving that an operator actually possesses this calculus is not often an easy task. The criterion (2.7) is hardly quantitative. For that, let us return back to the general situation of the $\omega$-bisectorial operator $T$ on a Hilbert space $\mathscr{H}$. At this level of generality, a remarkable theorem due to McIntosh comes to our aid. It harnesses the perspective that we took in §2.2 of the reconstruction of the signal in norm as a criterion for the $\mathrm{H}^{\infty}$ functional calculus. The theorem we allude to is the following, although we note that in its original form, it was proved for $\omega$-sectorial operators.

Theorem 2.2** (McIntosh’s Theorem).**

Let $T$ be an $\omega$-bisectorial operator, $\omega\in[0,\pi/2)$, on a Hilbert space $\mathscr{H}$. Then, $\mathscr{H}={\mathrm{ker}}(T)\oplus\overline{{\mathrm{ran}}(T)}$ and for $\mu\in(\omega,\pi/2)$, the following are equivalent:

1. (i)

there exists $C>0$ such that for all $\psi\in\Psi(\mathrm{S}^{\mathrm{o}}_{\mu})$,

[TABLE] 2. (ii)

there exists $\psi\in\Psi(\mathrm{S}^{\mathrm{o}}_{\mu})$, not identically zero on either sector, such that

[TABLE]

for all $u\in\overline{{\mathrm{ran}}(T)}$, 3. (iii)

for all $\psi\in\Psi(\mathrm{S}^{\mathrm{o}}_{\mu})$ not identically zero on either sector, (2.11) holds for all $u\in\overline{{\mathrm{ran}}(T)}$.

See [32] (page 221) for a proof of this theorem. It is also a central theme of the book [25] by Haase.

Of significance are the estimates (2.11), which are aptly called quadratic estimates in the literature. This is an incredibly useful quantitative criterion that provides a method with which to assert that $T$ has an $\mathrm{H}^{\infty}$ functional calculus. In the same vein as before, $\psi(tT)$ can be through of as a band-pass filter, localised about frequencies with real part $\frac{1}{t}$. In fact, the terminology “band-pass filter” was initially introduced to the author by Andreas Rosén in a conversation regarding the conceptual foundations of quadratic estimates for bisectorial operators. Our earlier exposition for operators with discrete spectrum in §2.2 was borne out of the desire to verify this perspective in a more concrete setting. Moreover, note that at this level of generality, there is no restriction on type of spectrum that the operator may admit. In particular, it may have both residue and continuous spectrum in addition to point spectrum.

In the classical setting, via the Fourier transform, the celebrated Plancherel’s theorem asserts that the square norm integral of a function is equal to the square norm integral of its Fourier transform. Via the band-pass filter perspective, it can be seen that (2.11) is also of this form. Similarly, in the classical setting, the Calderón reproducing formulas allow the reconstruction of a function through the use of two radial convolutions. Such reproducing formulas also exist at this level of generality through the $\mathrm{H}^{\infty}$ functional calculus using specific pairs of $\Psi(\mathrm{S}^{\mathrm{o}}_{\mu})$ functions. This is more or less immediate from the quadratic estimate perspective of the $\mathrm{H}^{\infty}$ functional calculus.

It is also worth pointing out a Hilbert space geometric interpretation of quadratic estimates. The estimates (2.11) actually define a new, but comparable Hilbert space norm on $\overline{{\mathrm{ran}}(T)}$. Therefore, the ability to assert an $\mathrm{H}^{\infty}$ functional calculus is equivalent to saying that there are many comparable norms on $\mathscr{H}$ which are adapted to the operator. Given recent developments in synthetic notions of curvature in metric spaces, this interpretation might lend a way of accessing the $\mathrm{H}^{\infty}$ functional calculus via metric space geometric methods in future. For a more detailed description of recent developments in metric geometry, see the survey article [30] by Lott and references therein.

For a more detailed treatment of functional calculi for bisectorial and sectorial operators, we cite the survey paper [1] by Albrecht, Duong, and McIntosh, as well as [25]. There are also Banach space versions of the $\mathrm{H}^{\infty}$ functional calculus, initially developed by Cowling, Doust, McIntosh and Yagi in [22]. As this is a grossly inexhaustive list, we urge the reader to consult the references therein.

3. Geometry

3.1. Kato square root problem in the smooth setting

We begin our presentation of a selection of geometric problems resolved through the use of functional calculus by introducing the Kato square root problem in the smooth manifold setting. To that end, let $(\mathcal{M},\mathrm{g})$ be a smooth, complete Riemannian manifold. We emphasise that it need not be compact. The Laplacian with respect to $\mathrm{g}$ is then given by

[TABLE]

where $\operatorname{div}_{\mathrm{g}}:=-\nabla^{\ast,\mathrm{g}}$ and $\nabla:{\rm H}^{\rm 1}(\mathcal{M})\subset{\rm L}^{2}(\mathcal{M})\to{\rm L}^{2}({\rm T}^{\ast}\mathcal{M})$. Let $\mathrm{h}$ be another metric with $C\geq 1$ satisfying

[TABLE]

for all $u\in{\rm T}_{x}\mathcal{M}$.

Let us define an extended distance metric on the set of all such Riemannian metrics via

[TABLE]

It is easily verified that there exists $A\in{\rm C}^{\infty}\cap{\rm L}^{\infty}(\mathrm{Sym}({\rm T}^{\ast}\mathcal{M}\otimes{\rm T}\mathcal{M}))$ such that for all $u,v\in{\rm T}_{x}\mathcal{M}$,

[TABLE]

where $\theta(x):=\sqrt{\det A(x)}$. Then,

[TABLE]

See [11] where these claims are demonstrated in detail.

As before, let us write

[TABLE]

for accretive $B\in{\rm L}^{\infty}({\rm T}^{\ast}\mathcal{M}\otimes{\rm T}\mathcal{M})$ and accretive $b\in{\rm L}^{\infty}(\mathcal{M})$. Then, as before, from operator theoretic considerations alone, it is possible to prove that $\Pi_{B,b,\mathrm{g}}$ is an $\omega$-bisectorial operator. If $\Pi_{B,b,\mathrm{g}}$ has an $\mathrm{H}^{\infty}$ functional calculus, via a similar calculation as before, whenever $\mathrm{h}$ is another metric with with $\uprho_{\mathcal{M}}(g,h)\leq C$, where $C$ is fixed and depends on $(\mathcal{M},\mathrm{g})$, we obtain that

[TABLE]

In fact, more generally, for any $f\in\mathrm{H}^{\infty}(S_{\mu}^{o})$ with $\mu\in(\omega,\pi/2)$, we obtain that

[TABLE]

The significance here is that in (3.3), it is effectively the comparison of two heat kernels with respect to two different metrics, and this estimate says that they are continuous in an ${\rm L}^{\infty}$ sense under the change of metric. A direct, operator theoretic approach would yield that it is close in the ${\rm L}^{\infty}$ topology in first-derivatives of the change of metric.

The estimate (3.4) is obtained from the generalisation of the Kato square root problem to a geometric setting. They were initially established in the compact setting in [10] in 2002. More precisely, the authors of [10] prove the following.

Theorem 3.1** (Theorem 7.1 in [10]).**

Let $(\mathcal{M},\mathrm{g})$ be a smooth compact Riemannian manifold and suppose that $B\in{\rm L}^{\infty}({\rm T}\mathcal{M}\otimes{\rm T}^{\ast}\mathcal{M})$ and $b\in{\rm L}^{\infty}(\mathcal{M})$ such that there exists a $\kappa>0$ satisfying

[TABLE]

for almost-every $x\in\mathcal{M}$ and for all $v\in{\rm T}^{\ast}_{x}\mathcal{M}$. Then, the operator $\Pi_{B,b,\mathrm{g}}$ has an $\mathrm{H}^{\infty}$ functional calculus with the bound in the estimate depending only on the ${\rm L}^{\infty}$ norms of $B$, $b$ as well as $(\mathcal{M},\mathrm{g})$ and $\kappa$.

A similar problem, but with $\nabla$ replaced by

[TABLE]

with domain ${\mathrm{dom}}(S)={\rm H}^{\rm 1}(\mathcal{M})$, was considered by Morris in 2012 in [34]. For him, $\mathcal{M}\subset\mathbb{R}^{n+k}$ is an embedded submanifold. With this, he considers the inhomogeneous operator

[TABLE]

where here $\tilde{B}\in{\rm L}^{\infty}(\mathcal{M})\oplus{\rm L}^{\infty}({\rm T}\mathcal{M}\otimes{\rm T}^{\ast}\mathcal{M})$ and $\mathrm{g}$ is the induced metric on $\mathcal{M}$ from $\mathbb{R}^{n+k}$. Note that the nomenclature “inhomogeneous” arises from the fact that the first entry in the square of this matrix realises an operator

[TABLE]

from which its easy to see that this is a perturbation of the Laplacian through the introduction of lower order terms. He proves the following.

Theorem 3.2** (Theorem 1.1 in [34]).**

Let $\mathcal{M}\subset\mathbb{R}^{n+k}$ be a smooth, complete Riemannian submanifold and suppose that there exists $C<\infty$ such that the second fundamental form satisfies the bound $|\mathrm{II}(\mathcal{M})(x)|<C$. Moreover, suppose there exists $\kappa_{1},\kappa_{2}>0$ such that

[TABLE]

Then, $\tilde{\Pi}_{\tilde{B},\tilde{b},\mathrm{g}}$ has an $\mathrm{H}^{\infty}$ functional calculus with the implicit constant dependent on the ${\rm L}^{\infty}$ norms of $\tilde{B}$ and $\tilde{b}$ as well as $C$, $\kappa_{1}$, $\kappa_{2}$ and $\mathcal{M}$.

This theorem was subsequently improved to the intrinsic setting in [13] by McIntosh and the author to obtain the following.

Theorem 3.3** (Theorem 1.1 in [13]).**

Let $(\mathcal{M},\mathrm{g})$ be a smooth, complete Riemannian manifold and suppose there exists $C<\infty$ and $\kappa>0$ with

[TABLE]

Suppose that the following ellipticity condition holds: there exist $\kappa_{1},\kappa_{2}>0$ such that

[TABLE]

for all $u\in{\rm L}^{2}(\mathcal{M})$ and $v\in{\rm H}^{\rm 1}(\mathcal{M})$. Then $\tilde{\Pi}_{\tilde{B},\tilde{b},\mathrm{g}}$ has an $\mathrm{H}^{\infty}$ functional calculus with the implicit constant dependent on the ${\rm L}^{\infty}$ norms of $\tilde{B}$ and $\tilde{b}$ as well as $C$, $\kappa$, $\kappa_{1}$, $\kappa_{2}$ and $(\mathcal{M},\mathrm{g})$.

A remarkable geometric consequence of these theorems is that, given any $\tilde{\mathrm{g}}$ with $\uprho_{\mathcal{M}}(\mathrm{g},\tilde{\mathrm{g}})<\infty$, we are able to obtain the statement corresponding to (3.4) for $(\mathcal{M},\tilde{\mathrm{g}})$ in a sufficiently small ${\rm L}^{\infty}$ neighbourhood of $\tilde{\mathrm{g}}$. The size of this neighbourhood, in addition to the geometric assumptions on $(\mathcal{M},\mathrm{g})$, now also depends on $\uprho_{\mathcal{M}}(\mathrm{g},\tilde{\mathrm{g}})$. This observation was exploited in [11] to study the necessity of geometric assumptions needed to prove these theorems. In particular, it was shown that the curvature and injectivity radius bound assumptions in Theorem 3.3 can be dropped. This indicates that there is ample room for improving the methods used in the proofs.

As a concluding remark, we note that the last theorem, Theorem 3.3, is a special case of a more general Kato square root theorem on vector bundles in [13]. We have abstained from presenting the most general result, which requires some technicalities to facilitate its description, in order to make this exposition more accessible.

3.2. Rough metrics and a geometric flow tangential to the Ricci flow

In [11], a notion of a rough metric was defined. We first note that in order to facilitate the definition of this object, it is important to note that the notion of “measurable” and “null measure” on a manifold, even when it is noncompact, is independent of a particular choice of Riemannian metric. That is to say, we can say that a set is measurable if the intersection of the set with a coordinate chart is measurable against the pullback measure in that chart. It is easily verified that for any smooth or even continuous metric tensor $\mathrm{g}$, the notion of $\mu_{\mathrm{g}}$-measurable agrees with this notion of measurable. Similarly, a null measure set is where the intersected set is null via the pullback measure in each chart. This allows us to talk about measurable sections of tensorfields, and we define a rough metric to be a symmetric, measurable $(2,0)$-tensorfield for which at each $x\in\mathcal{M}$, there exists a coordinate chart $(\psi_{x},U_{x})$ with a constant $C_{x}=C(U_{x})\geq 1$ satisfying

[TABLE]

for almost-every $y\in U_{x}$.

It is unclear how to obtain a distance for this metric tensor since the length function is not well-defined. However, it induces a Borel-regular measure $\mu_{\mathrm{g}}$, defined via the usual expression $d\mu_{\mathrm{g}}(x)=\sqrt{\det\mathrm{g}(x)}\ d\mathscr{L}(x)$. As in the smooth setting, this is readily checked to be coordinate independent. As before, the notions of measurable that we remarked were independent of a particular Riemannian structure, despite the lack of regularity of $\mathrm{g}$, agrees with $\mu_{\mathrm{g}}$. Moreover, a set is of null measure in our sense if and only if it is null with respect to the $\mu_{\mathrm{g}}$ measure.

In this context, it is again possible to ask the question (3.3), and more generally (3.4), for $(\mathcal{M},\mathrm{g})$. As before, this follows from establishing an $\mathrm{H}^{\infty}$ functional calculus for $\Pi_{B,b,\mathrm{g}}$. From the functional calculus perspective, it can also be shown that if $\Pi_{B,b,\mathrm{g}}$ has an $\mathrm{H}^{\infty}$ functional calculus, then so does $\Pi_{B,b,\mathrm{h}}$ for any rough metric $\mathrm{h}$ that satisfies (3.1). This shows that such metric tensors are geometric invariances of the Kato square root problem, which further leads to questions about the existence of geometries which are counterexamples to this problem, and consequently, $\mathrm{H}^{\infty}$ functional calculus.

There are many natural geometries that are rough metrics. Take, for instance, a smooth metric $\mathrm{h}$ and a lipeomorphism (locally bi-Lipschitz map) $\varphi:\mathcal{M}\to\mathcal{M}$, and consider the metric $\mathrm{g}=\varphi^{\ast}\mathrm{h}$. The metric $\mathrm{g}$ is a rough metric, and in general it will only have measurable coefficients since the definition of this metric involves derivatives of the Lipschitz map $\varphi$. Note that even in $\mathbb{R}^{2}$, given a null measure set, there exists a Lipschitz function for which the non-differentiable points are contained in this null set. This demonstrates the existence of crude and non-trivial rough metrics even in the limited context of lipeomorphic pullbacks alone. In particular, the set of singularities can be a dense subset. From a geometric point of view, rough metrics are useful when a manifold can be given a smooth differentiable structure so that singularities can be represented purely in terms of a lack of regularity of the metric tensor. Such is the case, for instance, for a conical singularity with spherical cross sections.

Beyond being an object of interest for the reasons we have just outlined, rough metrics have become useful in the study of a geometric flow with non-smooth initial data. In [23], Gigli and Mantegazza defined a geometric flow on a smooth compact manifold $(\mathcal{M},\mathrm{g})$, which they ultimately generalise to measure-metric spaces satisfying a certain synthetic Ricci curvature bound. The flow in the smooth setting is obtained by solving for $\varphi_{t,x,v}$ in the continuity equation

[TABLE]

for each fixed $x\in\mathcal{M}$, $t>0$. Here, $\uprho^{\mathrm{g}}_{t}$ is the heat kernel of $\Delta_{\mathrm{g}}$, the operator $\operatorname{div}_{\mathrm{g},y}$ denotes the divergence operator acting on the variable $y$, the vector $v\in{\rm T}_{x}\mathcal{M}$, and ${\rm d}_{x}(\uprho^{\mathrm{g}}_{t}(x,y))(v)$ is the directional derivative of $\uprho^{\mathrm{g}}_{t}(x,y)$ in the variable $x$ in the direction $v$. They define a new family of metrics evolving in time by the expression

[TABLE]

Moreover, they demonstrate that this flow satisfies

[TABLE]

for almost-every $s$ along $\mathrm{g}$-geodesics $\gamma$. That is, this flow $t\mapsto\mathrm{g}_{t}$ is tangential (in this weak sense) to the Ricci flow. While the authors of [23] do not touch upon questions of regularity in their paper, they show that $x\mapsto\mathrm{g}_{t}(x)$ remains smooth for all $t>0$. It is of interest is to understand smoothing or non-smoothing properties of this flow in measure-metric settings. However, such a task is difficult in the general measure-metric space world given that we lack even a basic yardstick to make sense of reasonable notions of regularity. However, between the smooth Riemannian manifold setting and the metric space world, there is an entire jungle of non-smooth metric tensors against smooth differentiable structures. This is a fertile playground in which we can attempt to address regularity questions.

Given that (3.5) is in divergence form, it can certainly be understood in a weak sense. Starting with an initial rough metric, the best we can expect for the regularity of the heat kernel is ${\rm C}^{\alpha}$-regularity via the methods of Nash-Moser-de Giorgi. This is precisely what is done in [12]. However, to understand the continuity of $x\mapsto\mathrm{g}_{t}(x)$ on an open region where the heat kernel improves to ${\rm C}^{1}$, we are forced with a non-smooth perturbation problem with measurable coefficients. Note from (3.2), by choosing any smooth metric, equation (3.6) is equivalent to solving a divergence form equation with bounded measurable coefficients coming from the regularity of the initial metric $\mathrm{g}$, regardless of whether $\uprho^{\mathrm{g}}_{t}$ improves in regularity in an open region. In fact, in [12], via this method, the regularity question is reduced to showing the Kato square root problem for rough metrics. There, it is shown that the Kato square root problem can be solved for any rough metric on a compact manifold. More precisely, the following is proved in [12].

Theorem 3.4** (Theorem 3.4 and Theorem 4.3 in [12]).**

On a compact manifold $\mathcal{M}$ with a rough metric $\mathrm{g}$, the operator $\Pi_{B,b,\mathrm{g}}$ admits an $\mathrm{H}^{\infty}$ functional calculus. Consequently, if $\varnothing\neq\mathcal{N}\subset\mathcal{M}$ is an open set where the initial heat kernel satisfies $\uprho^{\mathrm{g}}_{t}\in{\rm C}^{1}(\mathcal{N}^{2})$, $\mathrm{g}_{t}$ as defined by (3.6) exists on $\mathcal{N}$ and it is continuous.

As far as the author is aware, this is the first and only known instance where the Kato square root problem has been applied to obtain regularity properties of a geometric flow.

3.3. Elliptic boundary value problems

Let $(\mathcal{M},\mathrm{g})$ be a smooth, compact Riemannian manifold with smooth boundary $\Sigma:=\partial\mathcal{M}$. Moreover, let $(\mathcal{E},\mathrm{h}^{\mathcal{E}}),(\mathcal{F},\mathrm{h}^{\mathcal{F}})\to\mathcal{M}$ be Hermitian vector bundles and ${\rm D}:{\rm C}^{\infty}(\mathcal{M};\mathcal{E})\to{\rm C}^{\infty}(\mathcal{M};\mathcal{F})$ a first-order elliptic differential operator.

Here, elliptic means that the principal symbol $\upsigma_{{\rm D}}(x,\xi):\mathcal{E}_{x}\to\mathcal{F}_{x}$ is invertible.

By ${\rm C}^{k}_{\rm c}(\mathcal{M};\mathcal{E})$, let us denote the $u\in{\rm C}^{k}(\mathcal{M};\mathcal{E})$ such that ${\rm spt}{\text{ }}u$ is compact.

In particular, we allow for the possibility that ${\rm spt}{\text{ }}u\cap\Sigma\neq\varnothing$. On the other hand, by ${\rm C}^{k}_{\rm cc}(\mathcal{M};\mathcal{E})$, we mean $u\in{\rm C}^{k}_{\rm c}(\mathcal{M};\mathcal{E})$ such that ${\rm spt}{\text{ }}u\subset\mathring{\mathcal{M}}$, i.e. ${\rm spt}{\text{ }}u\cap\Sigma=\varnothing$.

The restrictions $\mathcal{E}{{\lvert}}_{\Sigma}$ and $\mathcal{F}{{\lvert}}_{\Sigma}$ define smooth bundles on $\Sigma$. We remark that our discussion in this section also carries over to the noncompact setting (but still with compact boundary), but there, further assumptions need to be placed on the operator ${\rm D}$. We avoid conducting a conversation at this level of generality as it adds an unnecessary layer of complexity that detracts from the key features of this example.

By local considerations alone, there exists a formal adjoint ${\rm D}^{\dagger}:{\rm C}^{\infty}(\mathcal{M};\mathcal{F})\to{\rm C}^{\infty}(\mathcal{M};\mathcal{E})$. On writing ${\rm D}_{\mathrm{cc}}^{\dagger}={\rm D}^{\dagger}$ with ${\mathrm{dom}}({\rm D}_{\mathrm{cc}}^{\dagger})={\rm C}^{\infty}_{\rm cc}(\mathcal{M};\mathcal{F})$, we obtain a maximal operator ${\rm D}_{\max}$ and a minimal operator ${\rm D}_{\min}$ by

[TABLE]

It is easy to see that both are closed and densely-defined operators.

In order study boundary value problems for first-order differential operators, it is important to be able to define the boundary restriction map $u\mapsto u{{\lvert}}_{\Sigma}$, initially defined on ${\rm C}^{\infty}_{\rm c}(\mathcal{M};\mathcal{E})$, on the space ${\mathrm{dom}}({\rm D}_{\max})$. Moreover, it is essential that we understand the topology of the range of this map in order to obtain a bounded surjection. By abstract considerations, noting ${\mathrm{dom}}({\rm D}_{\min})$ is the kernel of this extended boundary restriction map, it is possible to show that this topology is unique. This allows for a complete description of all closed extensions of ${\mathrm{dom}}(D_{\min})$ as closed subspaces of the range of the boundary restriction map. An important motivation comes from index theory. As Atiyah-Patodi-Singer demonstrated in [2, 4, 3, 5], it is essential to be able to deal with non-local boundary conditions to phrase index theorems. It is also desirable to be able to understand how the index behaves under perturbations of boundary conditions. The panoramic perspective afforded through understanding the boundary restriction map on ${\mathrm{dom}}(D_{\max})$ is necessary to consider such questions.

In [16], Bär and Ballmann demonstrate how this can by done by computing the topology of the range of the boundary trace map on ${\mathrm{dom}}(D_{\max})$ via an associated elliptic first-order differential operator ${\rm A}$ on ${\rm L}^{2}(\Sigma;\mathcal{E})$. This adapted operator on the boundary has principal symbol

[TABLE]

where $\xi\in{\rm T}^{\ast}_{x}\Sigma$ and $x\mapsto\tau(x)$ is an interior pointing co-vectorfield on $\Sigma$. There are many operators that satisfy this condition, but in their treatise, Bär and Ballmann assume that the symbol $\upsigma_{{\rm A}}(x,\xi)$ is skew-adjoint. In this case, there are self-adjoint operators, unique up to the addition of a zeroth order operator, that satisfy (3.7). For such an operator, the range

[TABLE]

and $u\mapsto u{{\lvert}}_{\Sigma}$ is a bounded surjection in this topology. Boundary conditions are precisely closed subspaces $B\subset\check{\mathrm{H}}(A)$ and the associated closed operator ${\rm D}_{B,\max}$ to $B$ is ${\rm D}$ with domain

[TABLE]

Since the operator ${\rm A}$ is self-adjoint, elliptic, first-order, and the boundary $\Sigma$ is compact, it only admits discrete spectrum with orthogonal eigenspaces. Therefore, the analysis in [16] is carried out in the spirit of the Fourier series, heavily exploiting the orthogonality of eigenspaces. Although this condition on ${\rm A}$ seems restrictive, their results are very general and can be applied liberally to a wide class of operators, including all Dirac type operators.

A quintessential example is the operator $\not{{\rm D}}$, the Atiyah-Singer Dirac operator. This induces an operator ${\rm A}$ such that $\upsigma_{{\rm A}}(x,\tau(x))$ anti-commutes with ${\rm A}$. Then, the boundary condition

[TABLE]

is precisely the famed Atiyah-Patodi-Singer boundary condition used in the proof of their index theorem. In this situation, the operator $\not{{\rm D}}_{B_{\text{APS}}}$ is self-adjoint. Given $r\in\mathbb{R}$, we can also consider generalised APS conditions:

[TABLE]

When $r\neq 0$, the operator $\not{{\rm D}}_{B_{\text{APS}}(r)}$ is no longer self-adjoint.

A more exotic object of interest is the Rarita-Schwinger operator. Again, the manifold $\mathcal{M}$ is assumed to be Spin and this operator is defined as follows. Let $\mathcal{W}:={\rm T}^{\ast}\mathcal{M}\otimes\operatorname{\not{\Delta}}\mathcal{M}$ where $\operatorname{\not{\Delta}}\mathcal{M}$ is the Spin bundle and define the map $\iota:\operatorname{\not{\Delta}}\mathcal{M}\to\mathcal{W}$ via

[TABLE]

Then, let

[TABLE]

so that $\mathcal{W}=\iota(\operatorname{\not{\Delta}}\mathcal{M})\stackrel{{\scriptstyle\perp}}{{\oplus}}\operatorname{\not{\Delta}}_{\frac{3}{2}}\mathcal{M}.$ There, we obtain an induced twisted Dirac operator $\not{{\rm D}}_{\mathcal{W}}$ which can be written as the operator matrix

[TABLE]

The operator ${\rm D}_{\mathrm{RS}}$ is the Rarita-Schwinger operator. The principal symbol of any adapted boundary operator ${\rm A}_{\text{RS}}$ to ${\rm D}_{\mathrm{RS}}$, given by $\upsigma_{{\rm A}_{\text{RS}}}(x,\xi)=\upsigma_{{\rm D}_{\mathrm{RS}}}(x,\tau(x))^{-1}\,\circ\,\upsigma_{{\rm D}_{\mathrm{RS}}}(x,\xi)$, fails to be skew-symmetric in general dimensions. This is one operator that falls outside of the scope of the framework in [16].

In [17], Bär and the author demonstrate that any differential operator ${\rm A}$ obtained from a first-order elliptic differential operator ${\rm D}$ satisfying (3.7) is, up to the addition of a real constant, an invertible $\omega$-bisectorial operator. In this case, via the work of Grubb in [24], they obtain the spectral projectors $\chi^{\pm}({\rm A})$ to the left and right of the complex plane as classical pseudo-differential operators of order zero. With the aid of this, the operator

[TABLE]

is defined, where $\operatorname{sgn}(\zeta)=\chi^{+}(\zeta)-\chi^{-}(\zeta)$ and where $\chi^{\pm}$ were defined in (2.8). By the boundedness of $\chi^{\pm}({\rm A})$ and commutativity with ${\rm A}$, it is clear that ${\mathrm{dom}}(|{\rm A}|)={\mathrm{dom}}({\rm A})$. Moreover, the operator $|{\rm A}|$ is $\omega$-sectorial and invertible. As before, define

[TABLE]

Unfortunately, at this level of generality, the analysis is far more subtle than [16]. For instance, even the rudimentary task of demonstrating the boundedness of $u\mapsto u{{\lvert}}_{\Sigma}:{\mathrm{dom}}({\rm D}_{\max})\to\check{\mathrm{H}}({\rm A})$ can no longer be conducted as in [16] since the generalised eigenspaces associated to the spectrum of ${\rm A}$ may no longer be orthogonal. It turns out, however, that the analysis precisely reduces to establishing that $|{\rm A}|$ has an $\mathrm{H}^{\infty}$ functional calculus.

The criterion that McIntosh’s Theorem 2.2 provides is useful as it allows one to move between the $\mathrm{H}^{\infty}$ functional calculus and quadratic estimates. In the analysis in [17], it is really the quadratic estimates that are needed. These estimates are accessed via establishing the $\mathrm{H}^{\infty}$ functional calculus by other means. Namely, by a more recent work due to Auscher, Nahmod and McIntosh in [9], it is possible to assert the $\mathrm{H}^{\infty}$ functional calculus for an $\omega$-sectorial operator via interpolation methods. Via their Corollary 5.5 in [9], it suffices to show that there are two numbers $s,t>0$ such that ${\mathrm{dom}}(|{\rm A}|^{s})\subset{\mathrm{dom}}(|{\rm A}^{\ast}|^{s})$ and ${\mathrm{dom}}(|{\rm A}^{\ast}|^{t})\subset{\mathrm{dom}}(|{\rm A}|^{t})$. In our case, we can choose $s=t=1$ from the fact that ${\rm A}$ is a first-order elliptic differential operator on a smooth bundle $\mathcal{E}{{\lvert}}_{\Sigma}$ over a compact manifold $\Sigma$. More precisely, by pseudo-differential elliptic regularity theory, we obtain

[TABLE]

In fact, Corollary 5.5 in [9] can be thought of as a conduit between pseudo-differential methods and the $\mathrm{H}^{\infty}$ functional calculus. The following is an important result obtained in [17].

Theorem 3.5**.**

Let $(\mathcal{M},\mathrm{g})$ be a compact smooth manifold with compact boundary, and let ${\rm D}:{\rm C}^{\infty}(\mathcal{M};\mathcal{E})\to{\rm C}^{\infty}(\mathcal{M};\mathcal{F})$ be a first-order elliptic differential operator between Hermitian bundles $(\mathcal{E},\mathrm{h}^{\mathcal{E}})\to\mathcal{M}$ and $(\mathcal{F},\mathrm{h}^{\mathcal{F}})\to\mathcal{M}$. Then, the following hold:

1. (i)

${\rm C}^{\infty}_{\rm c}(\mathcal{M};\mathcal{E})$* is dense in ${\mathrm{dom}}({\rm D}_{\max})$ with respect to the corresponding graph norm.* 2. (ii)

The trace map ${\rm C}^{\infty}_{\rm c}(\mathcal{M};\mathcal{E})\to{\rm C}^{\infty}(\Sigma;\mathcal{E})$ given by $u\mapsto u{{\lvert}}_{\Sigma}$ extends uniquely to a surjective bounded linear map ${\mathrm{dom}}({\rm D}_{\max})\to\check{\mathrm{H}}({\rm A})$. 3. (iii)

The spaces

[TABLE] 4. (iv)

The ${\rm L}^{2}$ inner product $\left\langle\cdot,\cdot\right\rangle$ extends to a perfect paring $\left\langle\cdot,\cdot\right\rangle_{\check{\mathrm{H}}(A)\times\check{\mathrm{H}}(-A^{\ast})}$ between $\check{\mathrm{H}}(A)$ and $\check{\mathrm{H}}(-A^{\ast})$. For all $u\in{\mathrm{dom}}({\rm D}_{\max})$ and $v\in{\mathrm{dom}}(({\rm D}^{\dagger})_{\max})$,

[TABLE]

The corresponding statements hold for ${\rm D}^{\dagger}$ with the roles $\mathcal{E}$ and $\mathcal{F}$ interchanged and on replacing $A$ by $\tilde{A}$, a boundary adapted operator to ${\rm D}^{\dagger}$.

To illustrate through an example how the $\mathrm{H}^{\infty}$ functional calculus appears in these proofs, let us consider the following simple calculation. As we have aforementioned, these methods can be made to work for noncompact manifolds but with compact boundary. For the purpose of the calculation which we are about to perform, it is beneficial for us to consider $\mathcal{M}=\mathbb{R}_{+}\times\Sigma$ for some closed manifold $\Sigma$.

Despite that it is noncompact, it affords us with a greater simplicity than choosing a compact example. Moreover, further assume that ${\rm D}=\partial_{t}+{\rm A}$ over some bundle $\mathcal{E}=\mathcal{F}$ with ${\rm A}$ invertible $\omega$-bisectorial first-order differential operator. In this setting, the operator $|{\rm A}|$ has an $\mathrm{H}^{\infty}$ functional calculus.

Let $\eta\in{\rm C}^{\infty}_{\rm c}(\mathbb{R})$ be a cutoff function such that $\eta=1$ on $[0,1]$ and [math] on $[2,\infty)$. Then, for $v\in{\rm C}^{\infty}_{\rm c}(\Sigma)$, define

[TABLE]

By the properties of the semigroup, we have that $(\mathscr{E}v){{\lvert}}_{\Sigma}(x)=(\mathscr{E}v)(x,0)=v(x)$. Our goal is to show that

[TABLE]

which is a crucial step in proving ii of Theorem 3.5. So, first note that

[TABLE]

Moreover, write $v_{\pm}=\chi^{\pm}({\rm A})v$ so that $v=v_{-}+v_{+}$. This enables us to consider the cases $v_{\pm}$ separately from linearity of (3.9) and the boundedness of the projectors $\chi^{\pm}({\rm A})$.

Since $|{\rm A}|u={\rm A}u$, from (3.11), it follows that

[TABLE]

Thus, it follows easily from the fact that $\|\eta^{\prime}\|_{{\rm L}^{\infty}}<\infty$ that

[TABLE]

On setting $\psi(\zeta)=\zeta^{\frac{1}{2}}\mathrm{e}^{-\zeta}$, we have that $\psi\in\Psi(S_{\mu}^{o})$ for $\mu\in(0,\pi)$ and

[TABLE]

Therefore, since we assert that $|{\rm A}|$ has an $\mathrm{H}^{\infty}$ functional calculus, via McIntosh’s Theorem 2.2, we obtain that

[TABLE]

Note that the ultimate equivalence in this calculation follows from elliptic regularity theory.

For the remaining case, we have that $|{\rm A}|v_{-}=-{\rm A}v_{-}$ and from (3.11),

[TABLE]

Therefore,

[TABLE]

The first term is bounded similar to (3.12), but for the second term,

[TABLE]

The penultimate and crucial inequality again relies upon the $\mathrm{H}^{\infty}$ functional calculus of $|{\rm A}|$. On combining (3.12) and (3.13) yields (3.10).

As a concluding remark, let us emphasise that despite the discreteness of the spectrum of ${\rm A}$, in this more general setting, this does not allow us to reduce the problem to considering individual generalised eigenspaces of ${\rm A}$. One issue is their potential non-orthogonality. This apparent inconvenience is, in fact, a boon, as it has inspired the development of methods that are robust enough to be applied to even more general settings. For instance, the functional calculus perspective developed in understanding the general elliptic operator situation for compact boundary can, in fact, be applied to adapted operators that are symmetric on the boundary when the boundary is non-compact. Under a set of mild hypotheses, this yields a self-adjoint operator. The spectrum may now contain residue and continuous blobs in addition to points, but the $\mathrm{H}^{\infty}$ functional calculus methods are completely insensitive and are able to accommodate for this situation. This calculus exists for these operators since every self-adjoint operator enjoys an $\mathrm{H}^{\infty}$ functional calculus. Through careful analysis, it is possible to prove a theorem, similar in spirit to Theorem 3.5, in the non-compact boundary case.

3.4. Spectral flows and Riesz continuity for the Atiyah-Singer Dirac operator

In §3.1, we demonstrated that an $\mathrm{H}^{\infty}$ functional calculus for the operator $\Pi_{B,b,\mathrm{g}}$ yields Lipschitz estimates (3.4) for small perturbations of the coefficients $(B,b)$ in the ${\rm L}^{\infty}$ topology. This is obtained via establishing a much stronger result, namely that $(A,a)\mapsto f(\Pi_{A,a,\mathrm{g}})$ is holomorphic in an ${\rm L}^{\infty}$ neighbourhood of $(B,b)$. The proof of this exploits the fact that $\Pi_{B,b,\mathrm{g}}$ admits the structure

[TABLE]

where $\Gamma$ is a nilpotent, closed, densely-defined operator, and the $\tilde{B}_{1}$ and $\tilde{B_{2}}$ appropriately encode the coefficients $B$ and $b$. In particular, this affords us with a Hodge-decomposition in ${\rm L}^{2}$. Many desirable consequences follow from (3.14) and it is unclear how these could be obtained without such structure.

A geometric operator that is also of physical significance is the Atiyah-Singer Dirac operator $\not{{\rm D}}_{\mathrm{g}}$ over the Spin bundle $\operatorname{\not{\Delta}}\mathcal{M}$ when the manifold $\mathcal{M}$ is Spin with a Riemannian structure given by $\mathrm{g}$. There are some natural questions that one can ask. The first is in the case that $\mathcal{M}$ has compact boundary $\Sigma=\partial\mathcal{M}$, and where $\mathcal{B}$ is a self-adjoint boundary condition for $\not{{\rm D}}_{\mathrm{g}}$, yielding an operator $\not{{\rm D}}_{\mathcal{B},\mathrm{g}}$. It is interesting to know what happens to the operator when we perturb the boundary condition $\mathcal{B}$ in an appropriate sense. A primary complication is that the domain ${\mathrm{dom}}(\not{{\rm D}}_{\mathcal{B},\mathrm{g}})$ moves when we perturb $\mathcal{B}$. In fact, in a private communication with Alan Carey, he revealed that the late Krzysztof Wojciechowski was present at the Australian National University in 2004, and he was excited by the first-order factorisation of the Kato square root problem as well as the corresponding problem for the Hodge Dirac operator in [10] which naturally accounts for moving domains. He wondered whether these methods maybe applicable to study perturbations of boundary conditions for Dirac operators. His ambition was to understand the perturbation of such operators in the Riesz topology, which has connections to index theory and spectral flows.

A second related question, perhaps simpler in nature, is on a manifold without boundary. Here, it is desirable to to understand what happens to $\not{{\rm D}}_{\mathrm{g}}$ under perturbation of the metric $\mathrm{g}$. An added complication in this scenario that the Spin bundle $\operatorname{\not{\Delta}}\mathcal{M}$, outside of a set of exceptional perturbations, itself might change when the metric moves.

Wojciechowski’s excitement is indeed justified. It is possible to consider these problems via methods motivated by the resolution of the Kato square root problem. However, unlike this problem which can be seen via (3.14), it is unclear the kind of structural features that $\not{{\rm D}}$ might enjoy. Consequently, we are forced to consider the Lipschitz estimate directly, and for this, it is actually necessary to go beyond the regime of a single bisectorial operator and instead consider the functional calculus of two self-adjoint operators. We will provide a brief exposé of the details later. First, let us present some results from [14] by Rosén, McIntosh and the author where the metric perturbation question was addressed, and [15] by Rosén and the author where the boundary condition perturbation problem was studied.

Let us consider the question of metric perturbation. As we have aforementioned, the bundle itself can change here. However, if we are interested in understanding spectral properties, we are able to pull back through a similarity transform using a unitary operator. This allows us to compare the spectra of two operators despite the fact that they may act on two different Spin bundles. More precisely, if we have another metric $\mathrm{h}$ on $\mathcal{M}$, we can construct a unitary map $U:({\rm T}\mathcal{M},\mathrm{g})\to({\rm T}\mathcal{M},\mathrm{h})$, which naturally extends to a map between Spin bundles. In what is to follow, let

[TABLE]

The metric perturbation theorem is then the following.

Theorem 3.6** (Theorem 3.1 in [14]).**

Let $\mathcal{M}$ be a smooth Spin manifold (without boundary) with a smooth, complete metric $\mathrm{g}$ and Levi-Civita connection ${\nabla}^{\mathrm{g}}$. Fixing constant $C_{\mathrm{h}}>0$, let $\mathrm{h}$ be a ${\rm C}^{0,1}$ metric with $\uprho_{\mathcal{M}}(\mathrm{g},\mathrm{h})\leq 1$ satisfying:

1. (i)

there exists $\kappa>0$ such that $\operatorname{inj}(\mathcal{M},\mathrm{g})\geq\kappa$, 2. (ii)

there exists $C_{R}>0$ such that $|{\rm Ric}_{\mathrm{g}}|\leq C_{R}$ and $|{\nabla}^{\mathrm{g}}{\rm Ric}_{\mathrm{g}}|\leq C_{R}$, 3. (iii)

$|{\nabla}^{\mathrm{g}}\mathrm{h}|\leq C_{\mathrm{h}}$* almost-everywhere.*

Then, for $\omega\in(0,\pi/2)$, $\sigma>0$, whenever $f\in\mathrm{H}^{\infty}(\mathrm{S}^{\mathrm{o}}_{\omega,\sigma})$, we have the perturbation estimate

[TABLE]

The implicit constant depends on $\dim\mathcal{M}$ and the constants appearing in i-iii.

Returning to the original question posed by Wojciechowski, when the manifold $\mathcal{M}$ has boundary $\Sigma$, we can consider perturbations of local boundary conditions for the Atiyah-Singer Dirac operator using the framework in [16]. A local boundary condition is a space

[TABLE]

where $\mathcal{E}$ is a smooth subbundle. The operator $\not{{\rm D}}$ with boundary condition $\mathcal{B}$, denoted $\not{{\rm D}}_{\mathcal{B}}$, is the operator with domain

[TABLE]

To describe our perturbation result, we require two additional conditions on the local boundary condition $\mathcal{B}$:

1. (i)

Self-adjointness, which by §3.5 in [16] occurs if and only if $\upsigma_{\not{{\rm D}}}(\vec{\mathrm{n}}^{\flat})$ maps the ${\rm L}^{2}$ closure of $\mathcal{B}$ onto its orthogonal complement. Here $\vec{\mathrm{n}}^{\flat}$ the unit normal interior covectorfield. 2. (ii)

$\not{{\rm D}}$-ellipticity, which is defined in terms of a self-adjoint boundary operator $\not{{\rm A}}$ adapted to $\not{{\rm D}}$ with principal symbol $\upsigma_{\not{{\rm A}}}(\xi)=\upsigma_{\not{{\rm D}}}(\vec{\mathrm{n}}^{\flat})^{-1}\,\circ\,\upsigma_{\not{{\rm D}}}(\xi)$, and for which the operator

[TABLE]

is a Fredholm operator. Here, $\pi_{\mathcal{B}}:{\rm L}^{2}(\operatorname{\not{\Delta}}\Sigma)\to\mathcal{B}$ is projection induced from the fibrewise orthogonal projection $\pi_{\mathcal{E}}:\operatorname{\not{\Delta}}\Sigma\to\mathcal{E}$, and \raisebox{0.0pt}{\chi}_{[0,\infty)}(\not{{\rm A}}) is the projection onto the positive spectrum of the operator $\not{{\rm A}}$ (see Theorem 3.15 in [16]).

For two local boundary conditions $\mathcal{B}$ and $\tilde{\mathcal{B}}$, following §2 in Chapter IV in [28], we define the ${\rm L}^{\infty}$-gap between the subspaces $\mathcal{B}$ and $\tilde{\mathcal{B}}$ as

[TABLE]

where $\pi_{\mathcal{E}}$ and $\pi_{\tilde{\mathcal{E}}}$ are the orthogonal projections from $\operatorname{\not{\Delta}}\Sigma$ to $\mathcal{E}$ and $\tilde{\mathcal{E}}$ respectively. We let

[TABLE]

and similarly for $\tilde{\mathcal{B}}$. For a set $Z\subset\mathcal{M}$ which is a neighbourhood of $\Sigma$, given $r>0$, we write $Z_{r}=\left\{x\in\mathcal{M}:\uprho_{\mathrm{g}}(x,Z)<r\right\}$. By $Z_{r}\sqcup Z_{r}$, we denote the double of $Z_{r}$ by pasting along $\Sigma$.

Theorem 3.7** (Theorem 3.1 in [15]).**

Let $(\mathcal{M},\mathrm{g})$ be a smooth, Spin manifold with smooth, compact boundary $\Sigma=\partial\mathcal{M}$ that is complete as a metric space and suppose:

1. (i)

there exists a precompact open neighbourhood $Z$ of $\Sigma$ and $\kappa>0$ such that $\operatorname{inj}(\mathcal{M}\setminus Z,\mathrm{g})>\kappa$, 2. (ii)

there exists $C_{R}<\infty$ such that $|{\rm Ric}_{\mathrm{g}}|\leq C_{R}$ and $|{\nabla}{\rm Ric}_{\mathrm{g}}|\leq C_{R}$ on $\mathcal{M}\setminus Z$, and 3. (iii)

fix any smooth metric $\mathrm{g}_{Z}$ on the double $Z_{4}\sqcup Z_{4}$ obtained by pasting along $\Sigma$ and $C_{Z}<\infty$ and let $\kappa_{Z}>0$ with $|{\rm Ric}_{\mathrm{g}_{Z}}|\leq C_{Z}$ and $\operatorname{inj}(Z_{2}\sqcup Z_{2},\mathrm{g}_{Z})\geq\kappa_{Z}$.

Fixing $C_{B}<\infty$, let $\mathcal{B}$ and $\tilde{\mathcal{B}}$ be two local self-adjoint $\not{{\rm D}}$-elliptic boundary conditions which satisfy:

1. (iv)

$\|\mathcal{B}\|_{{\rm Lip}}+\|\tilde{\mathcal{B}}\|_{{\rm Lip}}\leq C_{B}$, and 2. (v)

$\not{{\rm D}}$-ellipticity constants of orders $1$ and $2$ for $\mathcal{B}$ in a given compact neighbourhood $K$ of the boundary.

Then, for $\omega\in(0,\pi/2)$ and $\sigma>0$, whenever we have $f\in\mathrm{H}^{\infty}(\mathrm{S}^{\mathrm{o}}_{\omega,\sigma})$, we have the perturbation estimate

[TABLE]

The implicit constant depends on $\dim\mathcal{M}$ and the constants appearing in (i)-(v).

The sole motivation for accounting for the function class $\mathrm{H}^{\infty}(\mathrm{S}^{\mathrm{o}}_{\omega,\sigma})$ is that it includes

[TABLE]

This allows us to measure the distance between the two operators in the Riesz topology. This is a desirable metric to work with in attempting to understand the spectral flow for these operators as it connects better with topological properties of these operators, as seen through K-theory. This was observed in [6] by Atiyah and Singer for bounded self-adjoint Fredholm operators. Moreover, note that for $\mathcal{M}$ compact, the geometric hypotheses in both Theorem 3.6 and Theorem 3.7 are automatically satisfied. A virtue of these results, even when applied to the compact setting, is that they quantify how these quantities enter the perturbation estimate.

In both these results, it is important to note that the continuity is in an ${\rm L}^{\infty}$ sense. That is, like for the perturbation estimate in the Kato square root problem, the continuity depends only in an ${\rm L}^{\infty}$ sense of the perturbation, although here, we are also required to perturb in a way in which we are uniformly ${\rm L}^{\infty}$ bounded in the gradient of the metric or the projector defining the local boundary condition. This is contrast to earlier results where an additive term of $\|\nabla^{\mathrm{g}}\mathrm{h}\|_{{\rm L}^{\infty}}$ or $\|\mathcal{B}\|_{{\rm Lip}}+\|\tilde{\mathcal{B}}\|$ would appear in the right side of the estimate.

A concrete example of a family of metrics $\mathrm{g}_{\varepsilon}$ on $\mathbb{R}^{2}$ that fits the hypothesis of Theorem 3.6 are:

[TABLE]

For $\varepsilon>0$, note that these metrics are not conformal to the standard Euclidean metric, but they are ${\rm C}^{0,1}$ with the property that $\|\nabla^{\mathbb{R}^{2}}\mathrm{g}_{\varepsilon}\|_{{\rm L}^{\infty}}\leq 1.$

As we have remarked earlier, in contrast to the resolution of the Kato square root problem, the analysis here is not to establish Lipschitz estimates through holomorphic dependency results for a single operator that encodes the perturbation. Instead, we are forced to consider the difference of functional calculus of two different operators. This is a complicated matter despite the fact that the operators are self-adjoint. Denoting the perturbation of either question by operators ${\rm D}_{1}$ and ${\rm D}_{2}$, and on drawing inspiration from methods arising from the resolution of the Kato square root problem, we are able to reduce the estimate of $\|f({\rm D}_{1})-f({\rm D}_{2})\|_{{\rm L}^{2}\to{\rm L}^{2}}$ to a local quadratic estimate of the form

[TABLE]

The coefficients $X$ either encode the difference between the metrics or the gap metric for the case of boundary condition perturbations. It is a bounded operator, a multiplication operator, or a special kind of singular integral operator. This reduction is obtained on drawing inspiration from McIntosh’s Theorem 2.2, where quadratic estimates are shown to imply (in fact, are equivalent to) the $\mathrm{H}^{\infty}$ functional calculus.

Note here that the upper limit of the integral is $1$ and not $\infty$ as in the case of the quadratic estimates (2.11) in Theorem 2.2. This is effectively because functions in $\mathrm{H}^{\infty}(\mathrm{S}^{\mathrm{o}}_{\omega,\sigma})$ are holomorphic on a neighbourhood of [math], and on recalling that $\psi(t{\rm D})$ can be thought of as a localisation of a signal to the spectrum of ${\rm D}$ with real part $\frac{1}{t}$, we can expect that there are only contributions from high frequencies.

4. Harmonic Analysis

In our exposition so far, we have already used the phrase “band-pass filter” in §2. This is suggestive, at least superficially, of potential links between functional calculus and harmonic analysis. Indeed, a particular significance of the quadratic estimates perspective of the $\mathrm{H}^{\infty}$ functional calculus is that it is a bridge to the realm of real-variable harmonic analysis from which these estimates can be calculated. In the two examples described in §3.1 and §3.2, we have seen that the the goal has been to obtain an $\mathrm{H}^{\infty}$ functional calculus via quadratic estimates. In §3.4, despite the fact that this problem is outside of the scope of the $\mathrm{H}^{\infty}$ functional calculus, there is a reduction to a “local” quadratic estimate. The aim of this section is to illustrate, using the philosophy which emerges from the resolution of the Kato square root problem, of how such an estimate might be computed.

Lets fix a scenario before we flesh out the intended outline. Fix $(\mathcal{M},\mathrm{g})$ to be a complete Riemannian manifold and suppose that ${\rm D}_{1}$ and ${\rm D}_{2}$ are first-order differential operators on a Hermitian bundle $(\mathcal{E},\mathrm{h}^{E})\to\mathcal{M}$. For simplicity, assume and that they are self-adjoint in ${\rm L}^{2}(\mathcal{M};\mathcal{E})$. We remark that, since we allow for vector bundles, the first-order nature of these operators are not a limitation. This can be seen in the setting of §3.1, where the problem is actually second-order for functions. Here, we have factorised this problem into a first-order system involving functions and vectorfields. Even in the situation of $(\mathcal{M},\mathrm{g})=(\mathrm{S}^{2},\mathrm{g}_{\text{round}})$, the two sphere with the round metric, the tangent and the cotangent bundles are not globally trivialisable. This is seen easily from the hairy ball theorem which states that any continuous vectorfield on $\mathrm{S}^{2}$ must vanish at some point. This simple example illustrates the reason which the vector bundle setting is unavoidable.

At a first glance, we have seemingly complicated the situation by allowing for vector bundles. Worse still, a factorisation of an elliptic higher-order problem might result in a non-elliptic first-order system. However, this tall cost is easily justified. A feature of a first-order operator ${\rm D}$ is that the principal symbol, or equivalently the commutator of ${\rm D}$ with functions, is a multiplication operator. Often, such operators are geometric, by which we mean that it arises from the geometry in some natural way. Quantitatively, this means that there is some $C>0$ such that

[TABLE]

This, along with an assumption $f{\mathrm{dom}}({\rm D})\subset{\mathrm{dom}}({\rm D})$ for $f\in{\rm C}^{\infty}_{\rm c}(\mathcal{M})$ and the self-adjointness of ${\rm D}$, is enough to show that this operator admits exponential off-diagonal estimates. That is, on fixing a $\psi\in\Psi(\mathrm{S}^{\mathrm{o}}_{\mu})$, there exists $C_{\psi}>0$ and for each $M>0$, there exists a constant $C_{\Delta,\psi,M}>0$ so that

[TABLE]

for every Borel set $E,\ F\subset\mathcal{M}$ and $u\in{\rm L}^{2}(\mathcal{M};\mathcal{E})$. These estimates can be found in Lemma 5.3 in [18] by Carbonaro, McIntosh and Morris for a specific choice of operator. A closer inspection of their proof reveals that the only necessary ingredients are the global boundedness of the commutator as in (4.1) and the fact that the domain of the operator is preserved under multiplication by compactly supported smooth functions.

The significance of such an estimate is that it allows us to localise computations. However, the functional calculus is built out of holomorphic functions and they admit tails when localised. In order to ensure convergence over a sum of norms of localisations, an estimate of the form (4.2) is required. More precisely, given a disjoint cover $\left\{B_{j}\right\}$, it is desirable to localise by decomposing as

[TABLE]

where $S$ is typically a multiplication operator. It is clear from the term on the right where off-diagonal estimates assist in the computation. These estimates can be thought of as a replacement for Gaussian estimates, the bread and butter of the scalar valued second-order world. However, as we have seen, for a natural class of first-order operators, exponential off-diagonal decay is obtained gratis. Moreover, this is also suggestive that our analysis might be carried out in spaces with exponential growth such as negatively curved manifolds. We shall touch upon this shortly.

Without loss of generality, let us assume that we are attempting to prove quadratic estimates of the form

[TABLE]

where either $\tau=1$ or $\tau=\infty$. Here, the operator

[TABLE]

whereas the operator ${\rm\bf Q}_{t}=\psi(t{\rm D}_{1})$ for some $\psi\in\Psi(\mathrm{S}^{\mathrm{o}}_{\mu})$. This satisfies

[TABLE]

as well as the exponential off-diagonal estimates (4.2). The quintessential example of such an operator is

[TABLE]

In order to proceed, we need to also understand how the underlying measure-metric geometry of the manifold $(\mathcal{M},\mathrm{g})$ affects the ability to perform harmonic analysis. This means we need to exploit some notion of scale invariance coupled to the geometry in a useful manner. In Euclidean space, this is captured through dyadic cubes. This is just a decomposition of $\mathbb{R}^{n}$ into a grid, where at each scale $j\in\mathbb{Z}$, the cubes of the grid are of length $2^{j}$ and are a nested subdivision in a natural way of the cubes of scale $j+1$ with length $2^{j+1}$.

Such structures exist in more general settings, and for us, the most significant result is due to Christ in [20]. There, he constructs a dyadic structure that naturally takes into account the measure-metric structure of $(\mathcal{M},\mathrm{g})$. As observed in [33], his result is easily generalised to spaces of exponential volume growth. More precisely, this means that there exists $c_{E}\geq 1$, $\kappa,c>0$ such that

[TABLE]

for every $t\geq 1,\ r>0$ and $x\in\mathcal{M}$. For instance, if there exists $\eta\in\mathbb{R}$ such that the Ricci curvature satisfies ${\rm Ric}_{\mathrm{g}}\geq\eta\mathrm{g}$, then $(\mathcal{M},\mathrm{g})$ satisfies (4.4). Under this assumption, the theorem is as follows.

Theorem 4.1** (Existence of a truncated dyadic structure).**

Suppose that $(\mathcal{M},\mathrm{g})$ satisfies (4.4). Then, there exist countably many index sets $I_{k}$, a countable collection of open subsets $\left\{Q_{\alpha}^{k}\subset\mathcal{M}:\alpha\in I_{k},\ k\in\mathbb{N}\right\}$, points $z_{\alpha}^{k}\in Q_{\alpha}^{k}$ (called the centre of $Q_{\alpha}^{k}$), and constants $\delta\in(0,1)$, $a_{0}>0$, $\eta>0$ and $C_{1},C_{2}<\infty$ satisfying:

1. (i)

for all $k\in\mathbb{N}$, $\mu(\mathcal{M}\setminus\cup_{\alpha}Q_{\alpha}^{k})=0$, 2. (ii)

if $l\geq k$, then either $Q_{\beta}^{l}\subset Q_{\alpha}^{k}$ or $Q_{\beta}^{l}\cap Q_{\alpha}^{k}=\varnothing$, 3. (iii)

for each $(k,\alpha)$ and each $l<k$ there exists a unique $\beta$ such that $Q_{\alpha}^{k}\subset Q_{\beta}^{l}$, 4. (iv)

$\operatorname{diam}Q_{\alpha}^{k}<C_{1}\delta^{k}$, 5. (v)

$\mathrm{B}(z_{\alpha}^{k},a_{0}\delta^{k})\subset Q_{\alpha}^{k}$, 6. (vi)

for all $k,\alpha$ and for all $t>0$, $\mu\left\{x\in Q_{\alpha}^{k}:d(x,\mathcal{M}\setminus Q_{\alpha}^{k})\leq t\delta^{k}\right\}\leq C_{2}t^{\eta}\mu(Q_{\alpha}^{k}).$

As in the Euclidean case, this theorem allows us to obtain a dyadic grid, suitability adjusted to the distance and volume structure of $(\mathcal{M},\mathrm{g})$, beneath some scale $\mathrm{t_{S}}>0$. Note this is in contrast to the Euclidean setting, where the dyadic grid is arbitrarily large. This is possible more generally in non-compact settings if the measure-metric structure of $(\mathcal{M},\mathrm{g})$ satisfies the stronger condition of doubling, rather than the weaker notion of exponential volume growth that we assume.

This dyadic decomposition is essential to the construction of a dyadic averaging operator,

[TABLE]

uniformly bounded for $t\leq\mathrm{t_{S}}$. Moreover, an essential feature of this operator is that, in an appropriate sense, it is “constant” for $x\in Q$ for each dyadic cube $Q$. This notion requires some mild bounded geometry assumptions on the bundle and the manifold, and we shall remark on this towards the end. Via this notion of constancy, along with the off-diagonal decay coupled to the dyadic decomposition, we can construct a certain principal part operator

[TABLE]

This is a kind of local smoothing operator and we will soon see that this is a key operator in the analysis. Using these operators, the underlying philosophy at the heart of the resolution of the Kato square root problem is to break up the integrand in (4.3) in the following way:

[TABLE]

The first term,

[TABLE]

aptly called the principal part term, is estimated by resorting to a Poincaré inequality, appropriately generalised to the bundle setting to match with the notion of constancy that we have aforementioned in passing. The second term,

[TABLE]

reduces to exploiting certain cancellation properties of the differential operator, and again, this is typically a mild term to estimate.

The key purpose of this decomposition is really to reduce the whole problem, through these two mild estimates, to a local Carleson measure estimate. While the first two estimates can be treated by “soft” methods, that is, more or less exploiting operator theoretic properties of the problem, it is in the local Carleson measure estimate where recent developments in modern real-variable harmonic analysis become of great importance.

To describe this in further detail, let us provide the definition of a local Carleson measure. A measure $\nu$ is called a local Carleson measure on $\mathcal{M}\times(0,t^{\prime}]$ (for some fixed $t^{\prime}\in(0,\mathrm{t_{S}}]$) if

[TABLE]

Here, $\mathrm{R}(Q)=Q\times(0,\operatorname{\ell}(Q))$ is the Carleson box over $Q$, $\operatorname{\ell}(Q)$ is the length of the cube $Q$, and $\mathscr{Q}_{t}$ are the dyadic cubes whose length $\operatorname{\ell}(Q)\sim t$. The norm $\|\nu\|_{\mathcal{C}}$ is the local Carleson norm of $\nu$. For a local Carleson measure $\nu$, Carleson’s inequality yields

[TABLE]

for all $u\in{\rm L}^{2}(\mathcal{M};\mathcal{E})$. A proof of this, in the case of functions, can be found as Theorem 4.2 in [34], but this proof is equally valid in our vector bundle setting.

Returning back to estimating the remaining term

[TABLE]

note that by what we have just mentioned, it suffices to prove that

[TABLE]

is a local Carleson measure. This typically involves heavy duty harmonic analysis machinery including non-tangential maximal functions, local $T(b)$ theorems, and other ideas arising from real-variable harmonic analysis. It is unfortunately beyond the scope of this article to give a description of these methods, but the books [36, 35] by Stein and [19] by Christ give detailed accounts of these ideas. The survey article [26] touches upon more recent developments of these ideas, particularly in connection with the Kato square root problem.

As we have aforementioned, this procedure for obtaining the quadratic estimates require mild bounded geometry assumptions on both the manifold and the bundle. For instance, in Theorem 3.3, Theorem 3.6, and Theorem 3.7, this is the reason that Ricci curvature and injectivity radius bounds appear as an assumption. We have already remarked that a bounded geometry assumption is required for the notion of constancy to which we alluded earlier. We also require this for a generalised local Poincaré inequality which is bootstrapped from functions. These are required to estimate (4.6) as well (4.7).

It is unclear how necessary these assumptions are to carry out quadratic estimates. These assumptions exist, in part, due to the fact that these methods are a generalisation of real-variable harmonic analysis machinery that has been developed in the Euclidean setting, where geometrically speaking, the bundles are flat and trivial. However, it is reasonable to expect, given that we now have an awareness of the geometric issues that arise in the analysis, that these methods may be better adapted to account for the geometry of the problem. Investigations along these lines may allow for the bounded geometry notion to be further weakened, or perhaps dispensed with entirely. While this is expected to be a gross undertaking that is technical in nature, this would be a worthwhile task as it would likely reveal deeper links between geometry and harmonic analysis through the medium of functional calculus.