CONICAL SQUARE FUNCTION ESTIMATES AND FUNCTIONAL CALCULI FOR PERTURBED HODGE-DIRAC OPERATORS IN Lp DOROTHEE FREY, ALAN MCINTOSH, AND PIERRE PORTAL Abstract. Perturbed Hodge-Dirac operators and their holomorphic functional calculi, as investigated in the papers by Axelsson, Keith and the second author, provided insight into the solution of the Kato square-root problem for elliptic operators in L2 spaces, and allowed for an extension of these estimates to other systems with applications to nonsmooth boundary value problems. In this paper, we determine conditions under which such operators satisfy conical square function estimates in a range of Lp spaces, thus allowing us to apply the theory of Hardy spaces associated with an operator, to prove that they have a bounded holomorphic functional calculus in those Lp spaces. We also obtain functional calculi results for restrictions to certain subspaces, for a larger range of p. This provides a framework for obtaining Lp results on perturbed Hodge Laplacians, generalising known Riesz transform bounds for an elliptic operator L with bounded measurable coefficients, one Sobolev exponent below the Hodge exponent, and Lp bounds on the square-root of L by the gradient, two Sobolev exponents below the Hodge exponent. Our proof shows that the heart of the harmonic analysis in L2 extends to Lp for all p ∈ (1, ∞), while the restrictions in p come from the operator-theoretic part of the L2 proof. In the course of our work, we obtain some results of independent interest about singular integral operators on tent spaces, and about the relationship between conical and vertical square functions. Mathematics Subject Classification (2010): 47A60, 47F05, 42B30, 42B37

Contents 1. Introduction 1.1. Acknowledgments 2. Preliminaries 2.1. Notation 2.2. Holomorphic functional calculus 2.3. Off-diagonal bounds 2.4. Tent spaces 2.5. Hardy spaces associated with bisectorial operators 2.6. Hodge-Dirac operators 3. Main results 4. Consequences 4.1. Differential forms. 4.2. Second order elliptic operators 4.3. First order systems of the form DA Date: August 13, 2014. 1

2 4 4 4 5 6 7 8 9 14 19 19 20 22

2

DOROTHEE FREY, ALAN MCINTOSH, AND PIERRE PORTAL

5. 6. 7. 8. 9. 10.

Low frequency estimates: The Carleson measure argument High frequency estimates for p ∈ (2∗ , 2] Lp -L2 off-diagonal bounds High frequency estimates for p ∈ (max{1, (pH )∗ }, 2] Technical tools 1: Schur estimates in tent spaces Technical tools 2: Estimating conical square functions by vertical square functions References

23 26 27 30 33 36 41

1. Introduction In [18], Axelsson, Keith, and the second author introduced a general framework to study various harmonic analytic problems, such as boundedness of Riesz transforms or the construction of solutions to boundary value problems, through the holomorphic functional calculus of certain first order differential operators that generalise the Hodge-Dirac operator d + d∗ (where d is the exterior derivative) of Riemannian geometry. By proving that such Hodge-Dirac operators have a bounded holomorphic functional calculus in L2 , they recover, in particular, the solution of Kato’s square root problem obtained by Auscher, Hofmann, Lacey, McIntosh and Tchamitchian in [9]. Their results also provide the harmonic analytic foundation to new approaches to problems in PDE (see e.g. [5–7]) and geometry (see e.g. [15]). The main result in [18] is of a perturbative nature. Informally speaking, it states that the functional calculus of the standard Hodge-Dirac operator in L2 is stable under perturbation by rough coefficients. It is natural, and important in applications, to know whether or not such a result also holds in Lp for p ∈ (1, ∞). There are two main approaches to this question. The first one uses the extrapolation method pioneered by Blunck and Kunstmann in [20], and developed by Auscher in [3] to show that the relevant L2 bounds remain valid in certain intervals (p− , p+ ) about 2 which depend on the operator involved. This approach has been mostly developed to study second order differential operators, but has also been adapted to first order operators by Ajiev [1] and by Auscher and Stahlhut in [16,17]. The other approach to Lp estimates for the holomorphic functional calculus of Hodge-Dirac operators consists in adapting the entire machinery of [18] to Lp . This was done in the series of papers [31–33] by the second and third authors, together with Hytönen, using ideas from (UMD) Banach space valued harmonic analysis. At the technical level, all these results are fundamentally perturbation results for square function estimates. In L2 , the heart of [18] is an estimate of the form ˆ∞ ˆ ˆ∞ ˆ dxdt 1 2 −1 2 2 dxdt 12 ( |tΠB (I + t ΠB ) u(x)| ) .( |tΠ(I + t2 Π2 )−1 u(x)|2 )2 t t 0

Rn

0

∀u ∈ R(Γ),

Rn

where Π = Γ + Γ∗ is a first order differential (Hodge-Dirac) operator with constant coefficients, and ΠB = Γ + B1 Γ∗ B2 is a perturbation by L∞ coefficients B1 , B2 . See Section 2

CONICAL SQUARE FUNCTION ESTIMATES FOR DIRAC OPERATORS

3

for precise definitions. In Lp , the papers [31–33] establish analogues of the form ˆ∞ ˆ∞ dt 1 2 −1 2 2 dt 12 k( |tΠB (I + t ΠB ) u(.)| ) kp . k( |tΠ(I + t2 Π2 )−1 u(.)|2 ) 2 kp t t 0

∀u ∈ R(Γ).

0 p

While these (vertical) L square function estimates are traditionally used to establish the boundedness of the holomorphic functional calculus (see e.g. [24]), the same result could also be obtained using the conical Lp square function estimates: k(t, x) 7→ (tΠB (I+t2 ΠB 2 )−1 )M u(x)kT p,2 . k(t, x) 7→ (tΠ(I+t2 Π2 )−1 )M u(x)kT p,2

∀u ∈ R(Γ),

where M is a suitably large integer and T p,2 is one of Coifman-Meyer-Stein’s tent spaces (see [22] and Section 2 for precise definitions). This fact has been noticed in the development of a Hardy space theory associated with bisectorial operators (starting with [15,25,29], see also [34, Theorem 7.10]). In this paper, we prove such conical Lp square function estimates for the Hodge-Dirac operators introduced in [18]. This allows us to strengthen the results from [31–33] (in the scalar-valued setting) by eliminating the R-boundedness assumptions. Instead of relying on probabilistic/dyadic methods, we use the more flexible theory of Hardy spaces associated with operators, and recent results about integral operators on tent spaces. Our proof then exhibits an interesting phenomenon. As in [18] and other papers on functional calculus of Hodge-Dirac operators or Kato square root estimates, we consider separately the “high frequency" part of the estimate (involving k(t, x) 7→ (tΠB (I + t2 ΠB 2 )−1 )M (I + t2 Π2 )−M )u(x)kT p,2 ), and the “low frequency" part (involving k(t, x) 7→ (tΠB (I + t2 ΠB 2 )−1 )M (I − (I + t2 Π2 )−M )u(x)kT p,2 ). In L2 , the proof of the high frequency estimate is purely operator theoretic, while the low frequency requires the techniques from real analysis used in the solution of the Kato square root problem. In the approach to the Lp case given in [31–33], the same is true, but both the high and the low frequency estimate use an extra assumption: the R-bisectoriality of ΠB in Lp . With the approach through conical square function given here, we obtain the low frequency estimate for all p ∈ (1, ∞) without any assumption on the Lp behaviour of the operator ΠB . Restrictions in p, and appropriate assumptions (which are necessary, as can be seen in [3]), are needed for the high frequency part. We believe that this will be helpful in future projects, as the theory moves away from the Euclidean setting (see e.g. the work of Morris [40], Bandara and the second author [19]). Dealing with a specific Hodge-Dirac operator in a geometric context, one can hope to prove sharp high frequency estimates using methods specific to the context at hand, and combine them with the harmonic analytic machinery developed here to get the full square function estimates, and hence the functional calculus result. Another feature of the approach given here is that we obtain, from Lp assumptions, not just functional calculus results in Lp , but also functional calculus results on some subspaces of Lq for certain q < p. In particular, we obtain Riesz transform estimates for q ∈ (p∗ , 2], and reverse Riesz transform estimates for q ∈ (p∗∗ , 2]. Here p∗ and p∗∗ denote the first and second Sobolev exponents below p. This can also be relevant in geometric settings, where one expects the results to depend not only on the geometry, but on the

4

DOROTHEE FREY, ALAN MCINTOSH, AND PIERRE PORTAL

different levels of forms. The paper is organised as follows. In Section 2, we give the relevant definitions and recall the main results from the theories that this paper builds upon. In Section 3, we state our main results - relevant high and low frequency square function estimates - and establish their functional calculus consequences as corollaries in Section 4. In Section 5, we prove low frequency estimates by developing Lp conical square function versions of the tools used in [18]. In Section 6, we prove high frequency estimates for p ∈ (max(1, 2∗ ), 2]. In this range, the proof is straightforward, and does not require any Lp assumption. In dimensions 1 and 2 this already gives the result for all p ∈ (1, ∞). In Section 7 we establish the relevant Lp -L2 off-diagonal bounds for the resolvents of our Hodge-Dirac operator. In Section 8, we use these off-diagonal bounds to prove the high frequency estimates. This uses singular integral operator theory on tent spaces, and, in particular, results from the final two sections. The latter are devoted to proving some technical estimates required earlier in the paper. We think that the results proven there, including Schur-type extrapolation results for integral operator on tent spaces (Section 9), and a comparison of conical square functions by vertical square functions for bisectorial operators with appropriate off-diagonal bounds (Section 10), are of independent interest.

1.1. Acknowledgments. All three authors gratefully acknowledge support from the Australian Research Council through the Discovery Project DP120103692. This work is a key outcome of DP120103692. Frey and McIntosh also acknowledge support from ARC DP110102488. Portal is further supported by the ARC through the Future Fellowship FT130100607. The authors thank Pascal Auscher for stimulating discussions, and for keeping us aware of the progress of his student Sebastian Stahlhut on related questions. There is a connection between the results in [17] by Auscher and Stahlhut and our results, though the approaches are rather different because Auscher and Stahlhult rely on the results from [31–33] through [16], while one of our aims is to give an alternative approach to [31–33]. We remark that Auscher and Stahlhult apply their results to develop an extensive theory of a priori estimates for related non-smooth boundary value problems. 2. Preliminaries 2.1. Notation. Throughout the paper n and N denote two fixed positive natural numbers. We express inequalities “up to a constant" between two positive quantities a, b with the notation a . b. By this we mean that there exists a constant C > 0, independent of all relevant quantities in the statement, such that a ≤ Cb. If a . b and b . a, we write a ≈ b. We denote R∗ = R \ {0}. For a Banach space X, we write L(X) for the set of all bounded linear operators on X. For p ∈ (1, ∞) and an unbounded linear operator A on Lp (Rn ; CN ), we denote by Dp (A), Rp (A), Np (A) its domain, range and null space, respectively. We use upper and lower stars to denote Sobolev exponents: For p ∈ [1, ∞), we denote np np p∗ = n+p and p∗ = n−p , with the convention p∗ = ∞ for p ≥ n.

CONICAL SQUARE FUNCTION ESTIMATES FOR DIRAC OPERATORS

5

For a ball (resp. cube) B ⊆ Rn with radius (resp. side length) r > 0 and given α > 0, we write αB for the ball (resp. the cube) with the same centre and radius (resp. side length) αr. We define dyadic shells by S1 (B) := 4B and Sj (B) := 2j+1 B \ 2j B for j ≥ 2. 2.2. Holomorphic functional calculus. This paper deals with the holomorphic functional calculus of certain bisectorial first order differential operators. The fundamental results concerning this calculus have been developed in [14,24,36,38]. References for this theory include the lecture notes [2] and [37], and the book [27]. Definition 2.1. Let 0 ≤ ω < µ < π2 . Define closed and open sectors and double sectors in the complex plane by Sω+ := {z ∈ C : | arg z| ≤ ω} ∪ {0},

Sω− := −Sω+ ,

o Sµ+

o o := {z ∈ C : z 6= 0, | arg z| < µ}, Sµ− := −Sµ+ , o o Sω := Sω+ ∪ Sω− , Sµo := Sµ+ ∪ Sµ− .

Denote by H(Sµo ) the space of all holomorphic functions on Sµo . Let further H ∞ (Sµo ) := {ψ ∈ H(Sµo ) : kψkL∞ (Sµo ) < ∞}, Ψβα (Sµo ) := {ψ ∈ H(Sµo ) : ∃C > 0 : |ψ(z)| ≤ C|z|α (1 + |z|α+β )−1 ∀z ∈ Sµo } S for every α, β > 0, and set Ψ(Sµo ) := α,β>0 Ψβα (Sµo ). We say that ψ ∈ Ψ(Sµo ) is nono degenerate if neither of the restrictions ψ|Sµ± vanishes identically. Definition 2.2. Let 0 ≤ ω < π2 . A closed operator D acting on a Banach space X is called ω-bisectorial if σ(D) ⊂ Sω , and for all θ ∈ (ω, π2 ) there exists Cθ > 0 such that kλ(λI − D)−1 kL(X) ≤ Cθ

∀λ ∈ C \ Sθ .

We say that D is bisectorial if it is ω-bisectorial for some ω ∈ [0, π2 ). For D bisectorial with angle ω ∈ [0, π2 ) and ψ ∈ Ψ(Sµo ) for µ ∈ (ω, π2 ), we define ψ(D) through the Cauchy integral ˆ 1 ψ(D) = ψ(z)(zI − D)−1 dz, 2πi γ

where γ denotes the boundary of Sθ for some θ ∈ (ω, µ), oriented counter-clockwise. Definition 2.3. Let 0 ≤ ω < π2 and µ ∈ (ω, π2 ). An ω-bisectorial operator D, acting on a Banach space X, is said to have a bounded H ∞ functional calculus with angle µ if there exists C > 0 such that for all ψ ∈ Ψ(Sµ0 ) kψ(D)kL(X) ≤ Ckψk∞ . For such an operator, the functional calculus extends to a bounded algebra homomorphism from H ∞ (Sµo ) to L(X). More precisely, for all bounded functions f : Sµo ∪ {0} → C which are holomorphic on Sµo , one can define a bounded operator f (D) by f (D)u = f (0)PN (D) u + lim ψn (D)u, n→∞

u ∈ X,

6

DOROTHEE FREY, ALAN MCINTOSH, AND PIERRE PORTAL

where PN (D) denotes the bounded projection onto N (D) with null space R(D), and the functions ψn ∈ Ψ(Sµo ) are uniformly bounded and tend locally uniformly to f on Sµo ; see [2,24]. The definition is independent of the choice of the approximating sequence (ψn )n∈N . 2.3. Off-diagonal bounds. The operator theoretic property that captures the relevant aspect of the differential nature of our operators, is the following notion of off-diagonal bounds. This notion plays a central role in many current developments of singular integral operator theory. We refer to [3] for more information and references. Definition 2.4. Let p ∈ [1, 2]. A family of operators {Ut ; t > 0} ⊂ L2 (Rn ; CN ) is said to have Lp -L2 off-diagonal bounds of order M > 0 if there exists CM > 0 such that for all t > 0, all Borel sets E, F ⊆ Rn and all u ∈ Lp (Rn ; CN ) with supp u ⊆ F , we have  −M dist(E, F ) −n( p1 − 12 ) (2.1) kUt ukL2 (E) ≤ CM t 1+ kukLp (F ) , t where dist(E, F ) = inf{|x − y|; x ∈ E, y ∈ F }. The following properties of off-diagonal bounds with respect to composition and interpolation are essentially known (see [12]). We nonetheless include some proofs. Lemma 2.5. Let p ∈ (1, 2]. Let {Tt ; t > 0} ⊂ L(L2 (Rn ; CN )) have L2 -L2 off-diagonal bounds of every order, {Vt ; t > 0} ⊂ L(L2 (Rn ; CN )) have L2 -L2 off-diagonal bounds of order M > 0, and {St ; t > 0} ⊂ L(L2 (Rn ; CN )) have Lp -L2 off-diagonal bounds of order M , Then 1

1

(1) If sup ktn( p − 2 ) Tt kL(Lp ,L2 ) < ∞, then for all q ∈ (p, 2], {Tt ; t > 0} has Lq -L2 t>0

off-diagonal bounds of every order. (2) If {Tt ; t > 0} has Lp -Lq off-diagonal bounds of every order for some q ∈ [p, 2], then sup kTt kL(Lp ) < ∞. t>0

(3) {Vt St ; t > 0} has Lp -L2 off-diagonal bounds of order M . (4) For all t > 0, Tt extends to an operator Tt : L∞ (Rn ; CN ) → L2loc (Rn ; CN ) with n

kTt ukL2 (B(x0 ,t)) . t 2 kuk∞

∀u ∈ L∞ (Rn ; CN ), x0 ∈ Rn .

Proof. For (1), use Stein’s interpolation [42, Theorem 1] for the analytic family of opera1 1 ) M 0 (1−z) tors {tn( p − 2 )z (1+ dist(E,F ) 1E Tt 1F ; z ∈ S}, where S = {z ∈ C ; Re(z) ∈ [0, 1]} and t M 0 ∈ N. For q ∈ [p, 2], this gives Lq -L2 off-diagonal bounds of order max(0, M 0 (1 − which implies the result by choosing M 0 large enough.

1 − 12 q 1 − 12 p

)),

We refer to [3, Lemma 3.3] for a proof of (2). We now turn to (3). Let E, F ⊂ Rn be two Borel sets, and t > 0. Set δ = dist(E, F ) and G = {x ∈ Rn ; dist(x, F ) < 2δ }. Then dist(E, G) ≥ 2δ and dist(Rn \ G, F ) ≥ 2δ . Observe that the assumptions on Vt and St in particular imply that supt>0 kVt kL(L2 ) < ∞ 1 1 and supt>0 ktn( p − 2 ) St kL(Lp ,L2 ) < ∞ (taking E = F = Rn in the definition of off-diagonal

CONICAL SQUARE FUNCTION ESTIMATES FOR DIRAC OPERATORS

7

bounds). We have the following for all u ∈ Lp : k1E Vt St 1F uk2 ≤ k1E Vt 1G St 1F uk2 + k1E Vt 1Rn \G St 1F uk2 dist(E, G) −M ) k1G St 1F uk2 + k1Rn \G St 1F uk2 t 1 1 dist(E, G) −M dist(Rn \ G, F ) −M ) + (1 + ) )k1F ukp . t−n( p − 2 ) ((1 + t t 1 1 dist(E, F ) −M . t−n( p − 2 ) (1 + ) k1F ukp . t . (1 +

This proves (3). (4) The extension Tt : L∞ (Rn ; CN ) → L2loc (Rn ; CN ) can be defined as X 1Q (Tt u) = lim 1Q (Tt (1R u)), ρ→∞

R∈∆t dist(Q,R)<ρ

where u ∈ L∞ (Rn ; CN ), and Q ∈ ∆t a dyadic cube in Rn (see the beginning of Section 5 for a definition of ∆t ). It is shown in [18, Corollary 5.3] that the limit exists and the extension is well-defined.  The next lemma was shown in [34, Lemma 7.3] (as in [15, Lemma 3.6]). Lemma 2.6. Suppose M > 0, 0 ≤ ω < θ < µ < π2 . Let D be an ω-bisectorial operator in L2 (Rn ; CN ) such that {z(zI − D)−1 ; z ∈ C \ Sθ } has L2 -L2 off diagonal bounds of order M in the sense that kz(zI − D)−1 ukL2 (E) ≤ CM (1 + |z| dist(E, F ))−M kuk2 for all z ∈ C \ Sθ , all Borel subsets E, F ⊆ Rn , and all u ∈ L2 (Rn ; CN ) with supp u ⊆ F . If ψ ∈ Ψαβ (Sµo ) for α > 0, β > M , then {ψ(tD) ; t > 0} has L2 -L2 off-diagonal bounds of order M . 2.4. Tent spaces. Recall that the tent space T p,2 (Rn+1 + ), first introduced by Coifman, n+1 ∞ Meyer, and Stein in [22], is the completion of Cc (R+ ) with respect to the norm kF kT p,2

ˆ ˆ∞ ˆ =( ( Rn

|F (t, y)|2

1 dydt p 2 dx) p ) tn+1

0 B(x,t)

for p ∈ [1, ∞), and with respect to the norm ˆ rˆ −n sup (r kF kT ∞,2 = (r,x)∈R+ ×Rn

0

B(x,r)

|F (t, y)|2

dydt 1 )2 t

for p = ∞. The tent spaces interpolate by the complex method, in the sense that [T p0 ,2 , T p1 ,2 ]θ = T pθ ,2 + pθ1 . We recall a basic result about tent spaces, and another for θ ∈ [0, 1] and p1θ = 1−θ p0 about operators acting on them.

8

DOROTHEE FREY, ALAN MCINTOSH, AND PIERRE PORTAL

n+1 ∞ Lemma 2.7. [4] Let p ∈ [1, ∞), α ≥ 1 and Tαp,2 (Rn+1 + ) denote the completion of Cc (R+ ) with respect to the norm ˆ ˆ∞ ˆ 1 dydt p kF kTαp,2 = ( ( |F (t, y)|2 n+1 ) 2 dx) p . t Rn

0 B(x,αt)

p,2 Then Tαp,2 (Rn+1 (Rn+1 + ) = T + ) with the equivalence of norms n

kF kT p,2 ≤ kF kTαp,2 . α min{p,2} kF kT p,2

∀F ∈ T p,2 (Rn+1 + ).

Lemma 2.8. [34, Theorem 5.2] Let p ∈ (1, ∞). Let {Tt }t>0 be a family of operators n acting on L2 (Rn ) with L2 -L2 off-diagonal bounds of order M > min{p,2} . Then there exists n+1 C > 0 such that for all F ∈ T p,2 (R+ ) k(t, x) 7→ Tt F (t, . )(x)kT p,2 ≤ CkF kT p,2 . 2.5. Hardy spaces associated with bisectorial operators. We consider Hardy spaces associated with bisectorial operators. We refer to [15,25,28–30,34] and the references therein for more details about such spaces, and just recall here the main definition and result. Let 0 ≤ ω < µ < π2 , and D be an ω-bisectorial operator in L2 (Rn ; CN ) such that {(I + itD)−1 ; t ∈ R \ {0}} has L2 -L2 off-diagonal bounds of order M > n2 . Assume further that D has a bounded H ∞ functional calculus with angle θ ∈ (ω, µ). Given u ∈ L2 (Rn ; CN ) and ψ ∈ Ψ(Sµo ), write Qψ u(x, t) := ψ(tD)u(x),

x ∈ Rn , t > 0.

Definition 2.9. Let p ∈ [1, ∞), let ψ ∈ Ψ(Sµo ) be non-degenerate. The Hardy space p HD,ψ (Rn ; CN ) associated with D and ψ is the completion of the space N {u ∈ R2 (D) : Qψ u ∈ T p,2 (Rn+1 + ; C )}

with respect to the norm p kukHD,ψ := kQψ ukT p,2 .

Let us also recall [34, Theorem 7.10]: n

+ε Theorem 2.10. Let ε > 0. Let p ∈ (1, 2] and ψ, ψ˜ ∈ Ψε2 (Sµo ), or p ∈ [2, ∞) and ψ, ψ˜ ∈ Ψεn +ε (Sµo ), where µ > ω and both ψ and ψ˜ are non-degenerate. Then 2

p p p (1) HD,ψ (Rn ; CN ) = HD, (Rn ; CN ) =: HD (Rn ; CN ); ψ˜ p (2) For all u ∈ HD (Rn ; CN ), and all f ∈ Ψ(Sµo ), we have

k(t, x) 7→ ψ(tD)f (D)u(x)kT p,2 . kf k∞ kukHDp . p In particular, D has a bounded H ∞ functional calculus on HD (Rn ; CN ).

CONICAL SQUARE FUNCTION ESTIMATES FOR DIRAC OPERATORS

9

2.6. Hodge-Dirac operators. Throughout the paper, we work with the following class of Hodge-Dirac operators. It is a slight modification of the classes considered in [18] and [33]. Definition 2.11. A Hodge-DiracPoperator with constant coefficients is an operator of the ˆ j ∂j is a Fourier multiplier with symbol defined by form Π = Γ + Γ∗ , where Γ = −i nj=1 Γ ˆ = Γ(ξ) ˆ Γ =

n X

ˆ j ξj Γ

∀ξ ∈ Rn ,

j=1

ˆ j ∈ L (CN ), the operator Γ is nilpotent, i.e. Γ(ξ) ˆ 2 = 0 for all ξ ∈ Rn , and there with Γ exists κ > 0 such that (Π1)

ˆ κ|ξ||w| ≤ |Π(ξ)w|

ˆ ∀w ∈ R(Π(ξ)), ∀ξ ∈ Rn .

We list some results about these operators. Proposition 2.12. Suppose p ∈ (1, ∞). (1) The operator identity Π = Γ + Γ∗ holds in Lp (Rn ; CN ), in the sense that Dp (Π) = Dp (Γ) ∩ Dp (Γ∗ ) and Πu = Γu + Γ∗ u for all u ∈ Dp (Π). (2) There holds Np (Π) = Np (Γ) ∩ Np (Γ∗ ). (3) Π Hodge decomposes Lp (Rn ; CN ) in the sense that Lp (Rn ; CN ) = Np (Π) ⊕ Rp (Γ) ⊕ Rp (Γ∗ ) , or equivalently, Lp = Np (Γ) ⊕ Rp (Γ∗ ) and Lp = Np (Γ∗ ) ⊕ Rp (Γ). (4) Np (Γ), Np (Γ∗ ), Rp (Γ) and Rp (Γ∗ ) each form complex interpolation scales, p ∈ (1, ∞). (5) Hodge-Dirac operators with constant coefficients have a bounded H ∞ functional calculus in Lp (Rn ; CN ). (6) There holds k∇ ⊗ uk2 . kΠuk2 for all u ∈ D2 (Π) ∩ R2 (Π). Proof. See [33], Lemma 5.3, Proposition 5.4. For (4), see [31]. Part (5) is proven in [33, Theorem 3.6]. Part (6) is a consequence of (Π1), as shown in [33, Proposition 5.2].  We now consider perturbed Hodge-Dirac operators. Definition 2.13. A perturbed Hodge-Dirac operator is an operator of the form ΠB := Γ + Γ∗B := Γ + B1 Γ∗ B2 , where Π = Γ + Γ∗ is a Hodge-Dirac operator with constant coefficients, and B1 , B2 are multiplication operators by L∞ (Rn ; L (CN )) functions which satisfy Γ∗ B2 B1 Γ∗ = 0 in the sense that R2 (B2 B1 Γ∗ ) ⊂ N2 (Γ∗ ); ΓB1 B2 Γ = 0 in the sense that R2 (B1 B2 Γ) ⊂ N2 (Γ); Re(B1 Γ∗ u, Γ∗ u) ≥ κ1 kΓ∗ uk22 , Re(B2 Γu, Γu) ≥ κ2 kΓuk22 ,

∀u ∈ D2 (Γ∗ ) ∀u ∈ D2 (Γ)

and

10

DOROTHEE FREY, ALAN MCINTOSH, AND PIERRE PORTAL

for some κ1 , κ2 > 0. Let the angles of accretivity be π , 2 u∈R(Γ∗ )\{0} π ω2 := sup | arg(B2 u, u)| < , 2 u∈R(Γ)\{0}

ω1 :=

sup

| arg(B1 u, u)| <

and set ω := 12 (ω1 + ω2 ). Such operators satisfy the invertibility properties (denoting (ΠB (p))

kukp ≤ Cp kB1 ukp

1 p0

= 1 − p1 )

∀u ∈ Rp (Γ∗ ) and kvkp0 ≤ Cp0 kB2 ∗ vkp0

∀v ∈ Rp0 (Γ)

when p = 2. In many cases they satisfy (ΠB (p)) for all p ∈ (1, ∞), for example if B1 and B2 are invertible in L∞ , though in general all we can say is that the set of p for which (ΠB (p)) holds is open in (1, ∞). This follows on applying the extrapolation result of Kalton and Mitrea ([35], Theorem 2.5) to the interpolation families B1 : Rp (Γ∗ ) → Lp (Rn ) and 0 B2 ∗ : Rp0 (Γ) → Lp (Rn ). As noted in [33], it is a consequence of (ΠB (p)) that Γ∗B is a closed operator in Lp with ad0 joint (Γ∗B )∗ = B2 ∗ ΓB1 ∗ acting in Lp , that Rp (Γ∗B ) = B1 Rp (Γ∗ ), and that Rp0 (B2 ∗ ΓB1 ∗ ) = B2 ∗ Rp0 (Γ). Moreover, if (ΠB (p)) holds for all p in a subinterval of (1, ∞), then the spaces Rp (Γ∗B ) interpolate for those p also. Definition 2.14. A perturbed Hodge-Dirac operator ΠB Hodge decomposes Lp (Rn ; CN ) for some p ∈ (1, ∞), if (ΠB (p)) holds and there is a splitting into complemented subspaces Lp (Rn ; CN ) = Np (ΠB ) ⊕ Rp (ΠB ) = Np (ΠB ) ⊕ Rp (Γ) ⊕ Rp (Γ∗B ). It is proved in [18, Proposition 2.2] that ΠB Hodge decomposes L2 (Rn ; CN ). In investigating the property of Hodge Decomposition, let Pq denote the bounded projection of Lq (Rn ; CN ) onto Rq (Γ∗ ) with nullspace Nq (Γ), and let Qq denote the bounded projection of Lq (Rn ; CN ) onto Rq (Γ) with nullspace Nq (Γ∗ ) (1 < q < ∞). When ΠB Hodge decomposes Lp (Rn ; CN ), we denote by PRp (ΠB ) the projection of Lp (Rn ; CN ) onto Rp (ΠB ) with nullspace Np (ΠB ). Proposition 2.15. Let ΠB be a perturbed Hodge-Dirac operator, and let p ∈ (1, ∞). Then (i) ΠB Hodge decomposes Lp (Rn ; CN ) if and only if both (A) and (B) hold, where (A)

kukp . kB1 ukp

(B)

kvkp0 . kB2 ∗ vkp0

∀u ∈ Rp (Γ∗ ) and Lp (Rn ; CN ) = Np (Γ) ⊕ B1 Rp (Γ∗ ) ; 0

∀v ∈ Rp0 (Γ) and Lp (Rn ; CN ) = Np0 (Γ∗ ) ⊕ B2 ∗ Rp0 (Γ) .

(ii) Moreover (A) is equivalent to (A’), and (B) is equivalent to (B’) where (A’)

Pp B1 : Rp (Γ∗ ) → Rp (Γ∗ ) is an isomophism ;

CONICAL SQUARE FUNCTION ESTIMATES FOR DIRAC OPERATORS

(B’)

11

Qp0 B2 ∗ : Rp0 (Γ) → Rp0 (Γ) is an isomophism .

Proof. (i) Under the invertibility assumption (ΠB (p)), ΠB Hodge decomposes Lp (Rn ; CN ) if and only if both Lp (Rn ; CN ) = Np (Γ)⊕Rp (Γ∗B ) and Lp (Rn ; CN ) = Np (Γ∗B )⊕Rp (Γ) hold, 0 i.e. if and only if Lp (Rn ; CN ) = Np (Γ)⊕B1 Rp (Γ∗ ) and Lp (Rn ; CN ) = Np0 (Γ∗ )⊕B2 ∗ Rp0 (Γ) [33, Lemmas 6.1, 6.2]. This gives the proof of (i). (ii) (A) implies (A’): Let u ∈ Rp (Γ∗ ). Then Pp B1 u = −(I − Pp )B1 u + B1 u with, by (A), kPp B1 ukp ≈ k(I − Pp )B1 ukp + kB1 ukp , so that kukp . kB1 ukp . kPp B1 ukp . It remains to prove surjectivity. Let v ∈ Rp (Γ∗ ). By (A), there exist w ∈ Np (Γ) and u ∈ Rp (Γ∗ ) such that v = w + B1 u, and hence v = Pp v = Pp B1 u as claimed. (A’) implies (A): First we have that if u ∈ Rp (Γ∗ ), then kukp . kPp B1 ukp . kB1 ukp . Next we show that Np (Γ) ∩ B1 Rp (Γ∗ ) = {0}. Indeed if u ∈ Np (Γ), and u = B1 v with v ∈ Rp (Γ∗ ), then Pp B1 v = Pp u = 0, so by (A’), v = 0 and thus u = 0. Now we show that every element u ∈ Lp (Rn ; CN ) can be decomposed as stated. Let u ∈ Lp (Rn ; CN ). Then u = (I − Pp )u + Pp u for some v ∈ Rp (Γ∗ ) (by (A’))

= (I − Pp )u + Pp B1 v

= (I − Pp )(u − B1 v) + B1 v ∈ Np (Γ)

+

B1 Rp (Γ∗ )

with kB1 vkp . kvkp . kPp ukp . kukp . This gives the claimed direct sum decomposition. The proof that (B) is equivalent to (B’) follows the same lines, with p, Γ, B1 replaced by p0 , Γ∗ , B2 ∗ .  Proposition 2.16. The set of all p for which ΠB Hodge decomposes Lp (Rn ; CN ), is an open interval (pH , pH ), where 1 ≤ pH < 2 < pH ≤ ∞. Proof. By the interpolation properties of Rp (Γ∗ ), the set of p for which (A’) holds, is an open interval which contains 2, and the same can be said about (B’). So the set of all p for which ΠB Hodge decomposes Lp (Rn ; CN ) is the intersection of these two intervals, and thus is itself an open interval which we denote by (pH , pH ), with 1 ≤ pH < 2 < pH ≤ ∞.  An investigation of ΠB involves the related operator ΠB = Γ∗ + B2 ΓB1 , which is also a perturbed Hodge-Dirac operator with (Γ, Γ∗ , B1 , B2 ) replaced by (Γ∗ , Γ, B2 , B1 ), and for it we need the invertibility properties (ΠB (p))

kukp ≤ Cp kB2 ukp

∀u ∈ Rp (Γ) and kvkp0 ≤ Cp0 kB1 ∗ vkp0

∀v ∈ Rp0 (Γ∗ ) .

The formulae connecting ΠB and ΠB are, for θ ∈ (ω, π2 ), f ∈ H ∞ (Sθo ) and u ∈ D2 (Γ∗ ), f (ΠB )(Γ∗ u) = B2 f (ΠB )(B1 Γ∗ u), (2.2)





B1 g(ΠB )(Γ u) = g(ΠB )(B1 Γ u),

when f is odd, when g is even.

Proposition 2.17. Suppose ΠB is a perturbed Hodge-Dirac operator which Hodge decomposes Lp (Rn ; CN ) for all p ∈ (pH , pH ). Then:

12

DOROTHEE FREY, ALAN MCINTOSH, AND PIERRE PORTAL

(1) ΠB ∗ = Γ∗ + B2 ∗ ΓB1 ∗ is a perturbed Hodge-Dirac operator which Hodge decomposes Lq (Rn ; CN ) for all q ∈ ((pH )0 , (pH )0 ), i.e. (ΠB (q)) holds and Lq (Rn ; CN ) = Nq (ΠB ∗ ) ⊕ Rq (Γ∗ ) ⊕ Rq (B2 ∗ ΓB1∗ ) . (2) The perturbed Hodge-Dirac operator ΠB Hodge decomposes Lp (Rn ; CN ) for all p ∈ (pH , pH ), i.e. (ΠB (p)) holds and Lp (Rn ; CN ) = Np (ΠB ) ⊕ Rp (Γ∗ ) ⊕ Rp (B2 ΓB1 ) . (3) If, for some p ∈ (pH , pH ), ΠB is ω-bisectorial in Lp (Rn ; CN ), then ΠB is also ω-bisectorial in Lp (Rn ; CN ). Proof. (1) First note that the invertibility condition (ΠB (p)) for ΠB is the same as the invertibility condition (ΠB ∗ (p0 )) for ΠB ∗ = Γ∗ + B2 ∗ ΓB1 ∗ . Using this, it is proved in [33], 0 Lemma 6.3 that the Hodge decomposition for ΠB ∗ in Lp (Rn ; CN ) is equivalent to the Hodge decomposition for ΠB in Lp (Rn ; CN ). (2) On applying Proposition 2.15, ΠB Hodge decomposes Lp (Rn ; CN ) if and only if (A”)

Qp B2 : Rp (Γ) → Rp (Γ) is an isomophism and

(B”)

Pp0 B1 ∗ : Rp0 (Γ∗ ) → Rp0 (Γ∗ ) is an isomophism .

Using the Hodge decompositions for the unperturbed operators to identify the dual of Rp (Γ∗ ) with Rp0 (Γ∗ ), we find by duality that (A’) is equivalent to (B”) and (B’) is equivalent to (A”). This proves (2). (3) This is essentially proved in [33], Lemma 6.4.  Remark 2.18. We are not saying that (ΠB (p)) is equivalent to (ΠB (p)) for general p. We now define the operators RtB := (I + itΠB )−1 ,

t ∈ R, 1 B B ) = RtB R−t , t > 0, PtB := (I + t2 ΠB 2 )−1 = (RtB + R−t 2 1 2 −1 2 B QB = (−RtB + R−t ), t > 0. t := tΠB (I + t ΠB ) 2i In the unperturbed case B1 = B2 = I, we write Rt , Pt and Qt for RtB , PtB and QB t , B B B B B respectively. If we replace ΠB by ΠB , we replace Rt , Pt and Qt by Rt , Pt and QB t , respectively. We state some basic results for the unperturbed operator Π, noting that when we apply [33], we do not make use of the probabilistic/dyadic methods developed there. Proposition 2.19. Let M > n2 . (1) For all p ∈ (1, 2], {Rs M PRp (Π) ; s ∈ R∗ } has Lp -L2 off-diagonal bounds of every order. (2) For all p ∈ (1, ∞), k(s, x) 7→ Qs M u(x)kT p,2 ≈ kukp

∀u ∈ Rp (Π).

CONICAL SQUARE FUNCTION ESTIMATES FOR DIRAC OPERATORS

13

(3) For all p ∈ (max{2∗ , 1}, 2], k(s, x) 7→ Qs u(x)kT p,2 ≈ kukp

∀u ∈ Rp (Π).

Proof. Let s > 0. By [33, Proposition 4.8], Rs M PRp (Π) is a Fourier multiplier with bounded symbol ξ 7→ m(sξ). We also have that ΠM −1 R1 M is a Fourier multiplier with bounded b M −1 m(ξ). Since |ξ|M −1 |m(ξ)w| . |m(ξ)w| symbol m ˜ : ξ 7→ Π(ξ) ˜ for every ξ ∈ Rn and b every w ∈ R(Π(ξ)) (by (Π1), we have that kmk ˜ ∞ sup |m(ξ)| . M −1 . |ξ| ξ∈Rn \{0} For u ∈ L2 ∩ L1 , this implies n

n

bk2 . kmk2 kb uk∞ . kuk1 . ks 2 Rs M PRp (Π) uk2 . ks 2 m(s .) u Using Lemma 2.5 to interpolate this uniform bound with the L2 -L2 off-diagonal bounds gives (1). Using Corollary 10.2, Proposition 2.12, and Theorem 2.10, we then get k(s, x) 7→ Qs M u(x)kT p,2 . kukp

∀u ∈ Rp (Π),

for all p ∈ (1, ∞). The equivalence (2) then follows by duality. To prove (3), we remark that, for p ∈ (max(1, 2∗ ), ∞) and u ∈ Lp ∩ L2 , ˆ∞ ds k(t, x) 7→ Qt u(x)kT p,2 . k(t, x) 7→ Qt Qs Qs M u(x) kT p,2 . s 0 ( s (I − Pt )Ps if 0 < s ≤ t, Noting that Qt Qs = tt we consider the integral operator de(I − Ps )Pt if 0 < t ≤ s, s fined by ˆ∞ t s ds TK F (t, x) = min( , )K(t, s)F (s, x) ∀t > 0 ∀x ∈ Rn , s t s 0 ( (I − Pt )Ps if 0 < s ≤ t, for F ∈ T 2,2 and K(t, s) = Since, for every ε > 0, the integral (I − Ps )Pt if 0 < t ≤ s. ´∞ operator defined by T˜K F (t, x) = min( st , st )ε K(t, s)F (s, x) ds , for F ∈ T 2,2 and all t > 0, s 0

x ∈ Rn , is bounded on T 2,2 by Schur’s lemma, the result follows by Corollary 9.2 and (2).  We conclude the section by recalling the main result of Axelsson, Keith and the second author in [18]. Note that perturbed Hodge-Dirac operators satisfy the assumptions of [18] and [33]. In particular, k∇ ⊗ uk2 . kΠuk2 for all u ∈ D2 (Π) ∩ R2 (Π) as stated in Proposition 2.12 (6). Theorem 2.20. Suppose ΠB is a perturbed Hodge-Dirac operator with angles of accretivity as specified in Definition 2.13. Then: (1) ΠB is an ω-bisectorial operator in L2 (Rn ; CN ).

14

DOROTHEE FREY, ALAN MCINTOSH, AND PIERRE PORTAL

(2) The family {RtB ; t ∈ R} has L2 -L2 off-diagonal bounds of every order. (3) ΠB satisfies the quadratic estimate k(t, x) 7→ QB t u(x)kT 2,2 ≈ kukL2 for all u ∈ R2 (ΠB ) ⊆ L2 (Rn ; CN ). (4) For all µ > ω, ΠB has a bounded H ∞ functional calculus with angle µ in L2 (Rn ; CN ). 3. Main results Our main results are conical square function estimates on the range of Γ. Combining these estimates, and using the structure of Hodge-Dirac operators, we obtain functional calculus results as corollaries. Theorem 3.1. Suppose ΠB is a perturbed Hodge-Dirac operator. Given p ∈ (max{1, (pH )∗ }, ∞), we have M k(t, x) 7→ (QB t ) u(x)kT p,2 ≤ Cp kukp

∀u ∈ R2 (Γ) ∩ Lp (Rn ; CN ) and

M k(t, x) 7→ (QB t ) u(x)kT p,2 ≤ Cp kukp

∀u ∈ R2 (Γ∗ ) ∩ Lp (Rn ; CN ) ,

where M ∈ N if p ≥ 2, and M ∈ N with M >

n 2

if p < 2.

The proof of Theorem 3.1 consists of two parts: a low-frequency estimate and a highfrequency estimate, stated below in Theorem 3.8 and Theorem 3.9, respectively, and proven in the subsequent sections. We show that Theorem 3.1 is a consequence of Theorems 3.8 and 3.9 after their statements. M In the above theorem, when we consider a function of the form (t, x) 7→ (QB t ) u(x) in a p,2 B 2 tent space T , we are considering Qt as a bounded operator on L for each t. When p p ∈ (pH , pH ), QB t does in fact extend to a bounded operator on L .

As a first consequence, we obtain equivalence of the Hardy space HΠp B (Rn ; CN ) with the Lp closure of Rp (ΠB ) whenever p ∈ (pH , pH ), and corresponding results restricted to the ranges of Γ and Γ∗B for p below pH . We recall that (ΠB (p)) always holds for p ∈ (pH , pH ). Corollary 3.2. Suppose ΠB is a perturbed Hodge-Dirac operator. Suppose that µ ∈ (ω, π2 ), and that ψ ∈ Ψβα (Sµo ) is non-degenerate with either p ∈ (1, 2], and α > 0, β > n2 ;

or p ∈ [2, ∞), and α > n2 , β > 0.

(1) Let p ∈ (max{1, (pH )∗ }, pH ). Then k(t, x) 7→ ψ(tΠB )u(x)kT p,2 ≈ kukp

∀u ∈ Rp (Γ) and

k(t, x) 7→ ψ(tΠB )u(x)kT p,2 ≈ kukp

∀u ∈ Rp (Γ∗ ).

In particular, Rp (Γ) ⊆ HΠp B (Rn ; CN ) and Rp (Γ∗ ) ⊆ HΠp B (Rn ; CN ). (2) Let p ∈ (max{1, (pH )∗ }, pH ) and suppose that (ΠB (p)) holds. Then k(t, x) 7→ ψ(tΠB )u(x)kT p,2 ≈ kukp In particular, Rp (Γ∗B ) ⊆ HΠp B (Rn ; CN ).

∀u ∈ Rp (Γ∗B ).

CONICAL SQUARE FUNCTION ESTIMATES FOR DIRAC OPERATORS

15

(3) Let p ∈ (pH , pH ). Then k(t, x) 7→ ψ(tΠB )u(x)kT p,2 ≈ kukp

∀u ∈ Rp (ΠB ).

In particular, HΠp B (Rn ; CN ) = Rp (ΠB ). Remark 3.3. In Corollary 3.2 (2), one has in fact HΠp B (Rn ; CN ) = Rp (Γ) ⊕ Rp (Γ∗B ), for p ∈ (max{1, (pH )∗ }, pH ). This follows from the fact that the Hodge projections preserve Hardy spaces, as can be seen by considering their actions on HΠ1 B molecules (as defined in [13]). Remark 3.4. An inspection of our proof shows that we are actually proving that kukHΠp ≈ kukHΠp B

∀u ∈ R2 (Γ) ∩ HΠp .

When p > 1, we then use that kukHΠp ≈ kukLp . The proof still works if (pH )∗ < 1 and p = 1. In this case we get that kukHΠ1 ≈ kukHΠ1 B

∀u ∈ R2 (Γ) ∩ HΠ1 .

As Π is a Fourier multiplier one can then relate the HΠ1 norm to the classical H 1 norm: ∀u ∈ HΠ1 .

kukHΠ1 ≈ kukH 1 (Rn ,CN )

This can be done, for instance, by using the molecular theory presented in [13]. As a second consequence, we obtain functional calculus results for ΠB . Corollary 3.5. Suppose ΠB is a perturbed Hodge-Dirac operator. Suppose µ ∈ (ω, π2 ). (1) Let p ∈ (max{1, (pH )∗ }, pH ). Then for all f ∈ Ψ(Sµo ), kf (ΠB )ukp ≤ Cp kf k∞ kukp

∀u ∈ Rp (Γ).

(2) Let p ∈ (max{1, (pH )∗ }, pH ) and suppose (ΠB (p)) holds. Then for all f ∈ Ψ(Sµo ), kf (ΠB )ukp ≤ Cp kf k∞ kukp

∀u ∈ Rp (Γ∗B ).

(3) Let p ∈ (pH , pH ). Then ΠB is ω-bisectorial, and has a bounded H ∞ functional calculus with angle µ in Lp (Rn ; CN ). For the proofs, we use the following result, that establishes the reverse square function estimates when p < 2. Proposition 3.6. Suppose ΠB is a perturbed Hodge-Dirac operator. For all p ∈ [2, ∞] and all M ∈ N, we have (3.1)

M k(t, x) 7→ (QB t ) u(x)kT p,2 ≤ Cp kukp

∀u ∈ L2 (Rn ; CN ) ∩ Lp (Rn ; CN ).

Consequently, for all p ∈ (1, 2] and all M ∈ N, we have M kukp ≤ Cp k(t, x) 7→ (QB t ) u(x)kT p,2

∀u ∈ R2 (ΠB ) ∩ Lp (Rn ; CN ).

16

DOROTHEE FREY, ALAN MCINTOSH, AND PIERRE PORTAL

M Proof. The result for p = 2 holds by Theorem 2.20. We show that u 7→ (QB t ) u maps L∞ to T ∞,2 . The claim for p ∈ (2, ∞) then follows by interpolation. The argument goes back to Fefferman and Stein [26], and was used in a similar context in e.g. [10, Section 3.2]. Fix a cube Q in Rn and split u ∈ L∞ (Rn ; CN ) into u = u14Q + u1(4Q)c . Recall Sj (Q) = 2j+1 Q\2j Q for all j ≥ 2. Theorem 2.20 gives ˆ l(Q) ˆ 1 1 M 2 dxdt 21 ( |(QB ) . |Q|− 2 k14Q uk2 . kuk∞ . t ) 14Q u(x)| |Q| 0 t Q M of order N 0 > n2 yield On the other hand, L2 -L2 off-diagonal bounds for (QB t ) ˆ l(Q) ˆ 1 M 2 dxdt 12 ( |(QB ) t ) 1(4Q)c u(x)| |Q| 0 t Q ˆ l(Q) ˆ ∞ X 1 M 2 dxdt 21 |(QB ) . ( t ) 1Sj (Q) u(x)| |Q| 0 t Q j=2 ˆ l(Q) ∞ X t 2N 0 dt 1 −jN 0 1 . 2 ( ) k12j+1 Q uk22 ) 2 . kuk∞ . ( |Q| 0 l(Q) t j=2 0

Consider now p ∈ (1, 2). Let u ∈ R2 (ΠB ) ∩ Lp (Rn ; CN ) and v ∈ Lp (Rn ; CN ) ∩ L2 (Rn ; CN ). 0 We apply the above result to Π∗B in Lp (Rn ; CN ), noting that Π∗B Hodge decomposes Lq for all q ∈ ((pH )0 , (pH )0 ) by Proposition 2.17. By Calderón reproducing formula, tent space duality and the argument above, we have that ˆ ˆ ˆ∞ dt M B ∗ M |(QB | u(x).v(x)dx| . dx t ) u(x)||((Qt ) ) v(x)| t Rn

. k(t, x) 7→ . k(t, x) 7→

Rn 0 M (QB t ) u(x)kT p,2 k(t, x) M (QB t ) u(x)kT p,2 kvkp0 .

∗ M 7→ ((QB t ) ) v(x)kT p0 ,2



This gives the assertion.

M 2 Remark 3.7. Note that for the proof of (3.1), we only use that ((QB t ) )t>0 satisfies L n 2 0 2 2,2 L off-diagonal bounds of order N > 2 , and defines a bounded mapping from L to T . In particular, we do not use any assumptions on ΠB in Lp for p 6= 2. The proof gives a way to define a bounded extension from Lp to T p,2 of this mapping. In the case p = ∞, M 2 dxdt the above result shows that for every u ∈ L∞ (Rn ; CN ), |(QB is a Carleson t ) u(x)| t measure. We will make use of this fact in Proposition 5.5 below.

We next show that Corollaries 3.2 and 3.5 follow from Theorem 3.1 and Proposition 3.6. Proof of Corollary 3.2. By Theorem 2.10, it suffices to show M k(t, x) 7→ (QB t ) u(x)kT p,2 ≈ kukp , M and the corresponding equivalence for (QB in case (1), for M = 10n and u as given t ) in (1), (2) or (3). First suppose p ∈ (max{1, (pH )∗ }, 2]. Combining Theorem 3.1 and M Proposition 3.6 gives the equivalence for (QB and all u ∈ R2 (Γ) ∩ Lp (Rn ; CN ). The t ) M same reasoning applies to (QB and u ∈ R2 (Γ∗ ) ∩ Lp (Rn ; CN ). Now note that Np (Γ∗ ) ∩ t )

CONICAL SQUARE FUNCTION ESTIMATES FOR DIRAC OPERATORS

17

L2 (Rn ; CN ) ⊆ N2 (Γ∗ ), and the same holds with p and 2 interchanged. Using the Hodge decomposition for the unperturbed operator Π, we therefore have Lp (Rn ; CN ) ∩ L2 (Rn ; CN ) = [Np (Γ∗ ) ∩ N2 (Γ∗ )] ⊕ [Rp (Γ) ∩ R2 (Γ)]. Since this space is dense in Lp (Rn ; CN ), the space R2 (Γ) ∩ Lp (Rn ; CN ) is dense in Rp (Γ). M on Rp (Γ∗ ). As This gives the above equivalence on Rp (Γ), and, similarly, for (QB t ) stated before Definition 2.14, we have Rp (Γ∗B ) = B1 Rp (Γ∗ ) under (ΠB (p)). Using that M is even, the identity (2.2), kB1 k∞ < ∞, and that (ΠB (p)) holds by assumption, we therefore deduce from the above that, for u = B1 Γ∗ B2 v ∈ Rp (Γ∗B ), M B M ∗ k(t, x) 7→ (QB t ) u(x)kT p,2 = k(t, x) 7→ (Qt ) B1 Γ B2 v(x)kT p,2 M ∗ B M ∗ = k(t, x) 7→ B1 (QB t ) Γ B2 v(x)kT p,2 . k(t, x) 7→ (Qt ) Γ B2 v(x)kT p,2

. kΓ∗ B2 vkp . kukp . This gives (2). In the case p ∈ (pH , 2], ΠB Hodge decomposes Lp . This yields the result  on Rp (ΠB ). The case p ∈ [2, pH ) follows by duality, cf. the proof of Corollary 3.5. Proof of Corollary 3.5. First suppose p ∈ (max{1, (pH )∗ }, 2]. Let M = 10n, and µ ∈ (ω, π2 ). Let f ∈ Ψ(Sµo ) and u ∈ Rp (Γ). Using Corollary 3.2 and Theorem 2.10, we have that M kf (ΠB )ukp . k(t, x) 7→ (QB t ) f (ΠB )u(x)kT p,2 M . kf k∞ k(t, x) 7→ (QB t ) u(x)kT p,2 . kf k∞ kukp .

The same reasoning applies to u ∈ Rp (Γ∗B ), assuming (ΠB (p)). Now let p ∈ (pH , 2]. Since ΠB Hodge decomposes Lp , we have that, for all f ∈ Ψ(Sµo ), kf (ΠB )ukp . kf k∞ kukp

∀u ∈ Lp (Rn ; CN ).

This implies that ΠB is ω-bisectorial and has a bounded H ∞ functional calculus in Lp . Finally, we consider the case p ∈ [2, pH ). We apply the above result to Π∗B , which Hodge decomposes Lq for all q ∈ ((pH )0 , (pH )0 ) by Proposition 2.17. Hence, Π∗B has a bounded H ∞ functional calculus in Lq for all q ∈ ((pH )0 , 2]. By duality, ΠB has a bounded H ∞ functional calculus in Lp for all p ∈ [2, pH ).  The first conical square function estimate is an Lp version of the low frequency estimate in the main result of [18], Theorem 2.7 (and hence captures the harmonic analytic part of the proof of the Kato square root problem). The separation into low and high frequency is ˜ ˜ ˜ is a large natural number. Throughout done via the operators Pt N and I − Pt N where N ˜ = 10n. the paper we fix N Theorem 3.8. Suppose ΠB is a perturbed Hodge-Dirac operator. Suppose M ∈ N and p ∈ (1, ∞). Then ˜

N M k(t, x) 7→ (QB t ) Pt u(x)kT p,2 ≤ Cp kukp

∀u ∈ R2 (Π) ∩ Lp (Rn ; CN ).

18

DOROTHEE FREY, ALAN MCINTOSH, AND PIERRE PORTAL

This result is proven in Section 5. The second conical square function estimate is an Lp version of the high frequency estimate [18, Proposition 4.8, part (i)]. Note that this operator theoretic part of the proof in the case p = 2, is the part that does not necessarily hold for all p ∈ (1, ∞). Theorem 3.9. Suppose ΠB is a perturbed Hodge-Dirac operator. Suppose M = 10n and p ∈ (max{1, (pH )∗ }, 2]. Then ˜

N M k(t, x) 7→ (QB t ) (I − Pt )u(x)kT p,2 ≤ Cp kukp

∀u ∈ R2 (Γ) ∩ Lp (Rn ; CN ).

This result is proven in Sections 6, 7 and 8. We now show how to prove Theorem 3.1 from Theorems 3.8 and 3.9. Notice that the large (and somewhat arbitrary) value of M appearing in Theorem 3.9 is appropriately reduced as part of this proof. Proof of Theorem 3.1. For p ∈ (2, ∞), the claim has been shown in Proposition 3.6. From now on, suppose p ∈ (max{1, (pH )∗ }, 2]. Without loss of generality, we can assume that M = 10n. Indeed, the result for M > n2 will then follow by Theorem 2.10. Combining Theorem 3.8 and Theorem 3.9, we have that M k(t, x) 7→ (QB t ) u(x)kT p,2 . kukp

∀u ∈ R2 (Γ) ∩ Lp (Rn ; CN ).

Applying the same results to ΠB gives M k(t, x) 7→ (QB t ) u(x)kT p,2 . kukp

∀u ∈ R2 (Γ∗ ) ∩ Lp (Rn ; CN ). 

We conclude this section by showing that in certain situations the results can be improved when restricted to subspaces of the form Lp (Rn ; W ), where W is a subspace of CN . The proof given depends on Corollaries 7.3 and 8.3, as well as the preceding material. Theorem 3.10. Suppose ΠB is a perturbed Hodge-Dirac operator. Let W be a subspace c∗ (ξ)Γ(ξ) b b Γ c∗ (ξ) for all ξ ∈ Rn . of CN that is stable under Γ and Γ(ξ) (1) Suppose further that p ∈ (max{1, r∗ }, 2] for some r ∈ (1, 2), M ∈ N with M ≥ M 10n, and such that ΠB 2 is sectorial on Lr (Rn ; W ) and {(PtB ) 2 −1 PRr (ΠB 2 ) ; t ∈ R∗ } has Lr (Rn ; W ) − L2 (Rn ; CN ) off-diagonal bounds of every order. Then, for µ ∈ (ω, π2 ), we have (i)

kf (ΠB )Γukp . kf k∞ kΓukp

∀f ∈ H ∞ (Sµo ) ∀u ∈ Dp (Γ) ∩ L2 (Rn ; W ).

Moreover, Rp (Γ|Lp (Rn ;W ) ) ⊆ HΠp B (Rn ; CN ) with (ii)

kvkHΠp ≈ kvkp B

∀v ∈ Rp (Γ|Lp (Rn ;W ) ).

(2) If L2 (Rn ; W ) ⊂ R2 (Γ∗B ), (pH )∗ > 1, and (ΠB (r)) holds for all r ∈ ((pH )∗ , 2], then the hypotheses and conclusions of (1) hold for all p ∈ (max{1, (pH )∗∗ }, 2]. In particular, (iii)

k(ΠB 2 )1/2 ukp . kΓukp

∀u ∈ Dp (Γ) ∩ L2 (Rn ; W ).

CONICAL SQUARE FUNCTION ESTIMATES FOR DIRAC OPERATORS

19

Proof. (1) (i) The hypotheses of Corollary 8.3 are satisfied by assumption. We therefore obtain, for all p ∈ (max{1, r∗ }, 2], ˜

N M k(t, x) 7→ (QB t ) (I − Pt )Γv(x)kT p,2 . kΓvkp

∀v ∈ Dp (Γ) ∩ Lp (Rn ; W ).

Combined with Theorem 3.8, this gives M k(t, x) 7→ (QB t ) Γv(x)kT p,2 . kΓvkp

(3.2)

∀v ∈ Dp (Γ) ∩ Lp (Rn ; W ).

Therefore we have, for all f ∈ Ψ(Sµo ), v ∈ Dp (Γ) ∩ Lp (Rn ; W ), that (i) holds: M p kf (ΠB )Γvkp . k(t, x) 7→ (QB t ) f (ΠB )Γv(x)kT p,2 ≈ kf (ΠB )ΓvkHΠ

B

M . kΓvkHΠp ≈ k(t, x) 7→ (QB t ) Γv(x)kT p,2 . kΓvkp , B

where we have used Proposition 3.6, Theorem 2.10, and (3.2). The estimate holds for all f ∈ H ∞ (Sµo ) on taking limits as usual. (ii) This follows from (3.2) and the reverse inequality shown in Proposition 3.6. B (2) Let q > pH with q∗ > 1 and r ∈ (q∗ , q]. By Corollary 7.3, the identity PtB = RtB R−t and the fact that L2 (Rn ; W ) ⊂ R2 (Γ∗B ) by assumption, we have that 1

1

M

ktn( q∗ − 2 ) (PtB ) 2 −1 uk2 . kukLq∗ (Rn ;W )

∀t > 0 ∀u ∈ L2 (Rn ; W ) ∩ Lq∗ (Rn ; W ). M

Interpolating with L2 -L2 off-diagonal bounds as in Lemma 2.5(1), we get that {(PtB ) 2 −1 ; t > 0} has Lr (Rn ; W )-L2 off-diagonal bounds and that {PtB ; t > 0} has Lr (Rn ; W )-Lr˜ offdiagonal bounds for some r˜ ∈ (r, 2]. From Lemma 2.5(2), we then get that ΠB 2 is sectorial on Lr (Rn ; W ). This yields the hypotheses and hence the conclusions of (1). To obtain (iii), apply (i) with f (z) = sgn(z): k(ΠB 2 )1/2 ukp = k sgn(ΠB )ΠB ukp = k sgn(ΠB )Γukp ≤ Cp kΓukp ∀u ∈ Dp (Γ) ∩ L2 (Rn ; W ), noting that u ∈ R2 (Γ∗B ) by assumption.



4. Consequences 4.1. Differential forms. The motivating example for our formalism is perturbed differential forms, where CN = Λ = ⊕nk=0 Λk = ∧C Rn , the complex exterior algebra over Rn , and Γ = d, the exterior derivative, acting in Lp (Rn ; Λ) = ⊕nk=0 Lp (Rn ; Λk ). If the multiplication operators B1 , B2 satisfy the conditions of Definition 2.13, then ΠB = d + B1 d∗ B2 is a perturbed Hodge-Dirac operator, and it is from here that it gets its name. The Lp results stated in Section 3 all apply to this operator. Typically, but not necessarily, the operators Bj , j = 1, 2 split as Bj = Bj0 ⊕· · ·⊕Bjn , where Bjk ∈ L∞ (Rn ; L(Λk )), in which case Corollary 3.5 has a converse in the following sense (cf. [15, Theorem 5.14] for an analogous result for Hodge-Dirac operators on Riemannian manifolds). Proposition 4.1. Suppose ΠB = d+B1 d∗ B2 is a perturbed Hodge-Dirac operator as above, with Bj , j = 1, 2 splitting as Bj = Bj0 ⊕· · ·⊕Bjn , where Bjk ∈ L∞ (Rn ; L(Λk )), k = 0, . . . , n. Suppose that for some p ∈ (1, ∞), (ΠB (p)) holds and ΠB is an ω-bisectorial operator in Lp (Rn ; Λ) with a bounded H ∞ functional calculus in Lp (Rn ; Λ). Then p ∈ (pH , pH ).

20

DOROTHEE FREY, ALAN MCINTOSH, AND PIERRE PORTAL

We do not know if this converse holds for all perturbed Hodge-Dirac operators. It does, however, hold for all examples given in this section. See also the discussion before Corollary 10.2. Proof. We need to show that ΠB Hodge decomposes Lp (Rn ; Λ), i.e. Lp (Rn ; Λ) = Np (ΠB )⊕ Rp (Γ) ⊕ Rp (Γ∗B ), where Γ = d and Γ∗B = B1 d∗ B2 . Since ΠB is bisectorial in Lp (Rn ; Λ), we know that Lp (Rn ; Λ) = Np (ΠB ) ⊕ Rp (ΠB ). Therefore, it suffices to show that kΓukp + kΓ∗B ukp ≈ kΠB ukp

∀u ∈ Dp (ΠB ) = Dp (Γ) ∩ Dp (Γ∗B ).

For k = 0, . . . , n and u ∈ Lp (Rn ; Λ), denote by u(k) ∈ Lp (Rn ; Λk ) the k-th component of u. Note that Γ : Lp (Rn ; Λk ) → Lp (Rn ; Λk+1 ), Γ∗B : Lp (Rn ; Λk+1 ) → Lp (Rn ; Λk ), k = 0, . . . , n−1, and ΠB 2 : Lp (Rn ; Λk ) → Lp (Rn ; Λk ), k = 0, . . . , n. Using that sgn(ΠB ), where ( 1, if Re z > 0, sgn(z) = ∀z ∈ Sµ \ {0} and sgn(0) = 0, −1, if Re z < 0, is bounded in Lp (Rn ; Λ) since ΠB has a bounded H ∞ calculus, we therefore get for u ∈ Dp (ΠB ): kΓukp + kΓ∗B ukp ≈ ≈ ≈

n X j=0 n X k=0 n X

k(Γu)(j) kp +

X

k(Γ∗B u)(j) kp

j (k)

(kΓu kp +

kΓ∗B u(k) kp )



n X

(k)

kΠB u kp ≈

k=0

n X

k(ΠB 2 )1/2 u(k) kp

k=0

k((ΠB 2 )1/2 u)(k) kp ≈ k(ΠB 2 )1/2 ukp ≈ kΠB ukp .

k=0

 4.2. Second order elliptic operators. Let L denote the uniformly elliptic second order operator defined by n X Lf = −a div A∇f = −a ∂j (Aj,k ∂k f ) j,k=1 ∞

where a ∈ L (R ) with Re(a(x)) ≥ κ1 > 0 a.e. and A ∈ L∞ (Rn ; L(Cn )) with Re(A(x)) ≥ κ2 I > 0 a.e. Associated with L is the Hodge-Dirac operator   L2 (Rn ) 0 −a div A ⊕ acting in L2 (Rn ; (C1+n )) = ΠB = Γ+Γ∗B = Γ+B1 Γ∗ B2 = ∇ 0 2 L (Rn ; Cn ) n

where  Γ=

0 0 ∇ 0







, Γ =

0 − div 0 0



 , B1 =

a 0 0 0



 , B2 =

so that 2

ΠB =



L 0 ˜ 0 L



˜ = −∇a div A) . (where L

0 0 0 A

 ,

CONICAL SQUARE FUNCTION ESTIMATES FOR DIRAC OPERATORS

21

As shown in [18] (and recalled in Theorem 2.20), ΠB is an ω-bisectorial operator with an H ∞ functional calculus in L2 , so that in particular sgn(ΠB ) is a bounded operator on L2 (Rn ; C1+n ). Using the expression   0 −L−1/2 a div A 2 −1/2 sgn(ΠB ) = (ΠB ) ΠB = , ∇L−1/2 0 on D(ΠB ), and the fact that (sgn(ΠB ))2 u = u for all u ∈ R2 (ΠB ) = L2 (Rn ) ⊕ R2 (∇), we find that k∇L−1/2 gk2 ≈ kgk2 for all g ∈ R(L1/2 ), i.e. k∇f k2 ≈ kL1/2 f k2 for all f ∈ D(L1/2 ) = W 1,2 (Rn ), this being the Kato conjecture, previously solved in [9] (when a = 1). Turning now to Lp , we see that by our hypotheses, (ΠB (p)) holds for all p ∈ (1, ∞), and that {0} {0} Lp (Rn ) ⊕ , Rp (Γ∗B ) = ⊕ ⊕ , Rp (Γ) = . Np (ΠB ) = Np (div A) {0} Rp (∇) So p ∈ (pH , pH ), i.e. the Hodge decomposition Lp (Rn ; (C1+n )) = Np (ΠB )⊕Rp (Γ)⊕Rp (Γ∗B ) holds, if and only if Lp (Rn ; Cn ) = Np (div A) ⊕ Rp (∇). Turning briefly to Hardy space theory, we have HΠ2 B

L2 (Rn ) HL2 ⊕ = R2 (ΠB ) = = ⊕ HL2˜ R2 (∇)

and

HΠp B

=

HΠp B 2

HLp = ⊕ HLp˜

for all p ∈ (1, ∞). We remark that L has a bounded H ∞ functional calculus in HLp , and that sgn(ΠB ) is an isomorphism interchanging HLp and HLp˜ . We now state how the results of this section apply to ΠB , and have as consequences for L and its Riesz transform, results which are known, at least when a = 1 (see [3] and [30, Section 5]). Corollary 4.2. Let L = −a div A∇ be a uniformly elliptic operator as above. Then the following hold: (1) If pH < p < pH , then ΠB is an ω-bisectorial operator in Lp (Rn ; C1+n ) with a bounded H ∞ functional calculus. (2) If max{1, (pH )∗ } < p < pH , then HΠp B = Rp (ΠB ) and ΠB is an ω-bisectorial operator in Rp (ΠB ) with a bounded H ∞ functional calculus, so that L has a bounded H ∞ functional calculus in Lp (Rn ), and Dp (L1/2 ) = W 1,p (Rn ) with kL1/2 f kp ≈ k∇f kp . (3) If r ∈ (1, 2], max{1, min{r∗ , (pH )∗∗ }} < p < pH , M ∈ N with M ≥ 10n, and M {(I + t2 L)−( 2 −1) ; t > 0} has Lr (Rn )-L2 (Rn ) off-diagonal bounds of every order, then HLp˜ = Rp (∇) and kL1/2 f kp . k∇f kp for all f ∈ W 1,p (Rn ). Also g ∈ HLp if and only if ∇L−1/2 g ∈ Lp (Rn ; Cn ), with kgkHLp ≈ k∇L−1/2 gkp .

22

DOROTHEE FREY, ALAN MCINTOSH, AND PIERRE PORTAL

We remark that the hypotheses of (3) can also be stated in terms of off-diagonal bounds for the semigroup (e−tL )t>0 . Proof. As described above, ΠB is a perturbed Hodge-Dirac operator. (1) follows from Corollary 3.5 (3). (2) follows from Corollary 3.2 (1) and (2), noting that in our situation, the decomposition Rp (ΠB ) = Rp (Γ) ⊕ Rp (Γ∗B ) holds for all p ∈ (1, ∞). (3) Set W = C. As stated before, L2 (Rn ; W ) ⊆ R2 (Γ∗B ) = L2 (Rn ) ⊕ {0}. For w ∈ W and ξ ∈ Cn , we n P c∗ (ξ)Γ(ξ)w b b Γ c∗ (ξ)w = 0, so W is stable under Γ c∗ (ξ)Γ(ξ) b have that Γ = ( |ξj |2 )w and Γ(ξ) j=1

b Γ c∗ (ξ). If (pH )∗ > 1, we can therefore apply Theorem 3.10, which gives (3). If and Γ(ξ) (pH )∗ ≤ 1, (3) follows from (2).  Remark 4.3. If A = I, one has (pH , pH ) = (1, ∞), so the estimates in Corollary 4.2 hold for all p ∈ (1, ∞), in agreement with the results of [39] concerning L = −a∆. 4.3. First order systems of the form DA. Results for operators of the form DA or AD, used in studying boundary value problems as in [7], can be obtained in a similar way to those in this paper, building on the L2 theory in [8]. However they can also be obtained as consequences of the results for ΠB , as was shown in Section 3 of [18] when p = 2. Let us briefly summarise this in the Lp case. Let D be a first order system which is self-adjoint in L2 (Rn ; CN ), and A ∈ L∞ (Rn ; L(CN )) with Re(ADu, DU ) ≥ κkDuk2 2 for all u ∈ D2 (D). Set   L2 (Rn ; CN ) 0 ADA ∗ ∗ 2 n 2N ⊕ ΠB = Γ + ΓB = Γ + B1 Γ B2 = acting in L (R ; C ) = D 0 2 L (Rn ; CN ) where  Γ=

0 0 D 0





, Γ =



0 D 0 0



 , B1 =

A 0 0 0



 , B2 =

0 0 0 A

 .

Then ΠB is a Hodge-Dirac operator, and so, by [18], has a bounded H ∞ functional calculus in L2 (Rn ; C2N ). Turning to p ∈ (1, ∞), we find that (ΠB (p)) holds if and only if kukp . kAukp kvkp0 . kA∗ vkp0

∀u ∈ Rp (D)

and

∀v ∈ Rp0 (D) ,

and that (ΠB (p)) is the same. Assuming this (in particular if A is invertible in L∞ ), we find that {0} Np (D) Rp (AD) ∗ ⊕ ⊕ Np (ΠB ) = , Rp (Γ) = , Rp (ΓB ) = . ⊕ Np (DA) Rp (D) {0} and hence that ΠB Hodge decomposes Lp (Rn ; C2N ), i.e. p ∈ (pH , pH ), if and only if (4.1)

Lp (Rn ; CN ) = Np (DA) ⊕ Rp (D).

CONICAL SQUARE FUNCTION ESTIMATES FOR DIRAC OPERATORS

23

This can be seen following the arguments in Proposition 2.15: Under (ΠB (p)), (4.1) 0 holds if and only if Lp (Rn ; CN ) = Np0 (D) ⊕ Rp0 (A∗ D), i.e. if and only if Lp (Rn ; CN ) = Np (D) ⊕ Rp (AD). As in [18, Proof of Theorem 3.1] and [33, Corollary 8.17], we compute that, when defined,     A f (DA)u = 0 I f (ΠB ) u, I so that results concerning DA having a bounded H ∞ functional calculus in Lp (Rn ; CN ) can be obtained from our results for ΠB in Lp (Rn ; C2N ). Moreover results concerning bounds on f (DA)u when u ∈ Rp (D) can be obtained from our results on f (ΠB )v when v ∈ Rp (Γ) and on f (ΠB )w when w ∈ Rp (Γ∗B ). We leave further details to the reader, as well as consideration of AD. 5. Low frequency estimates: The Carleson measure argument In this section, we prove the low frequency estimate, Theorem 3.8. It suffices to show the result for M = 1, as for arbitrary M ∈ N, Lemma 2.8 and Theorem 2.20 yield ˜

˜

N N M B k(t, x) 7→ (QB t ) Pt u(x)kT p,2 . k(t, x) 7→ Qt Pt u(x)kT p,2

∀u ∈ R2 (Π) ∩ Lp (Rn ; CN ).

According to Theorem 2.20 and Lemma 2.5, the operator QB t extends to an operator ∞ n N 2 n N QB : L (R ; C ) → L (R ; C ) with t loc (5.1)

n

∀u ∈ L∞ (Rn ; CN ), x0 ∈ Rn , t > 0.

2 kQB t ukL2 (B(x0 ,t)) . t kuk∞

We can therefore define γt (x)w := (QB t w)(x)

(5.2)

∀w ∈ CN , x ∈ Rn ,

where, on the right-hand side, w is considered as the constant function defined by w(x) = w for all x ∈ Rn . Note that the definition of γt is different from the one in [18, Definition 5.1]. We use the splitting ˜

˜

˜

˜

N N N N B QB t Pt u = [Qt Pt u − γt At Pt u] + γt At Pt u, ˜

˜

˜

N N and refer to γt At Pt N u as the principal part, and [QB t Pt u − γt At Pt u] as the principal part approximation. S We use the following dyadic decomposition of Rn . Let ∆ = ∞ j=−∞ ∆2j , where ∆2j := j n n {2 (k + (0, 1] ) : k ∈ Z }. For a dyadic cube Q ∈ ∆2j , denote by l(Q) = 2j its sidelength, by |Q| = 2jn its volume. We set ∆t = ∆2j , if 2j−1 < t ≤ 2j . The dyadic averaging operator At : L2 (Rn ; CN ) → L2 (Rn ; CN ) is defined by ˆ 1 At u(x) := u(y) dy =: huiQx,t ∀u ∈ L2 (Rn ; CN ), x ∈ Rn , t > 0, |Qx,t | Qx,t

where Qx,t is the unique dyadic cube in ∆t that contains x.

24

DOROTHEE FREY, ALAN MCINTOSH, AND PIERRE PORTAL

Let us make the following simple observation: for all ε > 0, there exists a constant C > 0 such that for all t > 0 −(n+ε) X dist(Q, R) sup 1+ ≤ C. t Q∈∆t R∈∆ t

We first consider the principal part approximation, similar to [18, Proposition 5.5]. Proposition 5.1. Suppose ΠB is a perturbed Hodge-Dirac operator. Suppose p ∈ (1, ∞). Then ˜

˜

N N k(t, x) 7→ QB t Pt u(x) − γt (x)At Pt u(x)kT p,2 ≤ Cp kukp

∀u ∈ R2 (Π) ∩ Lp (Rn ; CN ).

Proof. Fix x ∈ Rn . For t > 0, we cover the ball B(x, t) by a finite number of cubes 2 2 Q ∈ ∆t . According to Theorem 2.20, QB t has L -L off-diagonal bounds of every order N 0 > 0. This, together with the Cauchy-Schwarz inequality and the Poincaré inequality (see [18, Lemma 5.4]), yields the following for Q ∈ ∆t : ˆ ∞ˆ ˜ ˜ N N 2 dydt 12 |QB ( ) t Pt u(y) − γt (y)At Pt u(y)| tn+1 0 Q ˆ ∞ˆ ˜ ˜ N N 2 dydt 12 |QB =( ) t (Pt u − hPt uiQ )(y)| tn+1 0 Q ˆ ∞ X 1 ˜ ˜ N N 2 dt ( kQB ≤( )2 t 1R (Pt u − hPt uiQ )kL2 (Q) ) n+1 t 0 R∈∆t ˆ ∞ X dist(Q, R) −N 0 dt 1 ˜ ˜ (1 + .( ( ) kPt N u − hPt N uiQ kL2 (R) )2 n+1 ) 2 t t 0 R∈∆t ˆ ∞ X dist(Q, R) −N 0 dt 1 ˜ ˜ .( (1 + ) kPt N u − hPt N uiQ k2L2 (R) n+1 ) 2 t t 0 R∈∆t ˆ ∞ˆ dist(Q, y) −N 0 +2n dydt 1 ˜ (1 + .( ) |t∇Pt N u(y)|2 n+1 ) 2 t t 0 Rn ˆ ˆ ∞ X ∞ dydt 1 0 ˜ . ( 2−j(N −2n) |t∇Pt N u(y)|2 n+1 ) 2 . t 0 Sj (Q) j=0 By change of angle in tent spaces, see Lemma 2.7, we thus get ∞ X j 0 ˜ ˜ ˜ N N B k(t, x) 7→ Qt Pt u(x) − γt (x)At Pt u(x)kT p,2 . 2− 2 (N −2n) k(t, x) 7→ t∇Pt N u(x)kT p,2 2j

j=0

.

∞ X

j

0

n

˜

˜

2− 2 (N −2n) 2j min{p,2} k(t, x) 7→ t∇Pt N u(x)kT p,2 . k(t, x) 7→ t∇Pt N u(x)kT p,2 ,

j=0 ˜

2n choosing N 0 > 2n + min{p,2} . Since Pt N is a Fourier multiplier, we have that, for u = Πv with v ∈ D2 (Π), and all j = 1, ..., n: ˜ ˜ t (∂x v) t∂xj Pt N u = Q j

CONICAL SQUARE FUNCTION ESTIMATES FOR DIRAC OPERATORS

25

˜ t = tΠPt N˜ . Therefore, by Proposition 2.19 and Proposition 2.12 (6), we have that with Q ˜ ˜ t (∂x v)(x)kT p,2 . max k∂x vkp . kukp , k(t, x) 7→ t∇Pt N u(x)kT p,2 . max k(t, x) 7→ Q j j j=1,...,n

j=1,...,n



which concludes the proof.

We now show that {γt At }t>0 defines a bounded operator on T p,2 for all p ∈ (1, ∞). This is an analogue of [18, Proposition 5.7]. Lemma 5.2. Suppose ΠB is a perturbed Hodge-Dirac operator. Suppose p ∈ (1, ∞). Then k(t, x) 7→ γt (x)At F (t, . )(x)kT p,2 ≤ Cp kF kT p,2

N ∀F ∈ T p,2 (Rn+1 + ; C ).

Proof. First observe that, given x ∈ Rn and t > 0, kAt F (t, . )kL∞ (B(x,t)) = sup |At F (t, y)| = sup

sup

|Q| |

F (t, z) dz|

Q∈∆t B(x,t)∩Q6=∅

y∈B(x,t)



ˆ −1

Q

1

n

|Q|− 2 kF (t, . )kL2 (B(x,5t)) . t− 2 kF (t, . )kL2 (B(x,5t)) .

Q∈∆t B(x,t)∩Q6=∅ n

According to (5.1), we have on the other hand kγt kL2 (B(x,t)) . t 2 , and consequently ˆ ∞ˆ ˆ ∞ˆ dydt 1 2 dydt 12 |γt (y) · At F (t, y)| n+1 ) . ( ( |F (t, y)|2 n+1 ) 2 . t t 0 B(x,t) 0 B(x,5t) Taking the Lp norm with respect to x ∈ Rn then yields the assertion.



˜

The corresponding estimate for the principal part γt (x)At Pt N relies on the following factorisation result for tent spaces: N Theorem 5.3 ([21], Theorem 1.1). Let p, q ∈ (1, ∞). If F ∈ T p,∞ (Rn+1 + ; C ) and G ∈ N p,q N T ∞,q (Rn+1 (Rn+1 + ; C ), then F G ∈ T + ; C ) and

kF · GkT p,q ≤ CkF kT p,∞ kGkT ∞,q , with a constant C which is independent of F and G. This plays the role of the Lp vertical square function version of Carleson’s inequality proven in [32, Lemma 8.1]. Note that this conical version is substantially simpler than its vertical counterpart. We also use the following conical maximal function estimate for operators with Lq -L2 off-diagonal bounds. Lemma 5.4. Let q ∈ [1, 2] and p ∈ (1, ∞) with q < p. Let {Tt }t>0 be a family of operators acting on L2 (Rn ; CN ) with Lq -L2 off-diagonal bounds of order N 0 > nq . Then k(t, x) 7→ At Tt u(x)kT p,∞ ≤ Cp kukp

∀u ∈ L2 (Rn ; CN ) ∩ Lp (Rn ; CN ).

26

DOROTHEE FREY, ALAN MCINTOSH, AND PIERRE PORTAL

Proof. Let u ∈ L2 (Rn ; CN )∩Lp (Rn ; CN ). Using Hölder’s inequality and Lq -L2 off-diagonal bounds for Tt , we obtain, given x ∈ Rn , the pointwise estimate ˆ 1 −n |Tt u(z)|2 dz) 2 sup |At Tt u(y)| . sup (t (y,t)∈Γ(x)

(y,t)∈Γ(x)



. .

sup

Qy,t ∞ X

(y,t)∈Γ(x) j=0 ∞ X

sup

(y,t)∈Γ(x) j=0 ∞ X

sup

ˆ

1

−n

|Tt 1Sj (Qy,t ) u(z)|2 dz) 2

(t

Qy,t n

t− q (1 +

2j t −N 0 ) kukLq (2j Qy,t ) t

n

0

n

2−j(N − q ) (2j t)− q kukLq (2j Qy,t )

(y,t)∈Γ(x) j=0

1 . sup ( B ball |B|

ˆ

1

|u(z)|q dz) q =: Mq u(x). B p

Since q < p, the boundedness of the Hardy-Littlewood maximal operator in L q implies that the maximal operator Mq is bounded in Lp (Rn ; CN ). Thus, k(t, x) 7→ At Tt u(x)kT p,∞ . kMq ukp . kukp .  The estimate for the principal part is a direct consequence of the two results above, together with the Carleson measure estimate for |γt (x)|2 dtdx . t Proposition 5.5. Suppose ΠB is a perturbed Hodge-Dirac operator. Let (t, x) 7→ γt (x) be defined as in (5.2). Suppose p ∈ (1, ∞). Then ˜

k(t, x) 7→ γt (x)At Pt N u(x)kT p,2 ≤ Cp kukp

∀u ∈ R2 (Π) ∩ Lp (Rn ; CN ).

Proof. Since A2t = At , Theorem 5.3 yields ˜

˜

k(t, x) 7→ γt (x)At Pt N u(x)kT p,2 . k(t, x) 7→ At Pt N u(x)kT p,∞ · k(t, x) 7→ γt (x)kT ∞,2 . The boundedness of the last factor is shown in Proposition 3.6 and noted in Remark 3.7, as a consequence of the L2 theory for ΠB established in [18], cf. Theorem 2.20. The first factor is bounded by a constant times kukp as an application of Lemma 5.4: take ˜ ˜ Tt := Pt N and notice that Pt N satisfies Lq -L2 off-diagonal bounds of every order for every q ∈ (1, 2] by Proposition 2.19.  6. High frequency estimates for p ∈ (2∗ , 2] In this section, we give a proof of Theorem 3.9 for the case 2∗ < pH < 2. In particular, this gives a proof for n ∈ {1, 2}, a case we have to exclude in Section 7 below for technical reasons. The proof is similar to the corresponding proof in L2 in [18], and is less technically involved than the case pH ≤ 2∗ considered in the next sections.

CONICAL SQUARE FUNCTION ESTIMATES FOR DIRAC OPERATORS

27

Proposition 6.1. Suppose ΠB is a perturbed Hodge-Dirac operator. Suppose M ∈ N and p ∈ (2∗ , 2]. Then ˜

N M k(t, x) 7→ (QB t ) (I − Pt )u(x)kT p,2 ≤ Cp kukp

∀u ∈ R2 (Γ) ∩ Lp (Rn ; CN ).

Proof of Proposition 6.1. Let p ∈ (2∗ , 2], M ∈ N, and u ∈ R2 (Γ) ∩ Lp (Rn ; CN ). Lemma 2.8 and Lemma 6.2 below yield k(t, x) 7→

M (QB t ) (I

˜ N

− Pt )u(x)kT p,2 = k(t, x) 7→

˜ −1 N X B Qt tΓ( Pt k )Qt u(x)kT p,2 k=0

. k(t, x) 7→ Qt u(x)kT p,2 . The assertion then follows from Proposition 2.19.



We use the following lemma in the proof of Proposition 6.1 above. The result and its proof are a slight modification of [18, Proposition 5.2]. 2 2 B Lemma 6.2. The families {tΓ∗B QB t ; t ∈ R} and {tΓQt ; t ∈ R} have L -L off-diagonal bounds of every order. B Proof. We prove the result for {tΓ∗B QB t ; t ∈ R}. The result for {tΓQt ; t ∈ R} then follows, given that for all t ∈ R, B ∗ B tΓQB t = (I − Pt ) − tΓB Qt ,

and {PtB ; t ∈ R} has L2 -L2 off-diagonal bounds of every order by Theorem 2.20. By B Theorem 2.20, we also have that the family {tΓ∗B QB t ; t ∈ R} = {PR(Γ∗B ) (I − Pt ) ; t ∈ R} is uniformly bounded in L2 . Let E, F ⊂ Rn be two Borel sets, u ∈ L2 (Rn ; CN ), and t ∈ R. As in [18, Proposition 5.2], let η be a Lipschitz function supported in E˜ = {x ∈ Rn ; dist(x, E) < 12 dist(x, F )}, constantly equal to 1 on E, and such that k∇ηk∞ ≤ 4 . We have the following: dist(E,F ) ∗ B ∗ B ∗ B ktΓ∗B QB t ukL2 (E) ≤ kηtΓB Qt uk2 ≤ k[ηI, tΓB ]Qt uk2 + ktΓB ηQt uk2 .

To estimate the first term, we use that [ηI, tΓ∗B ] = tB1 [ηI, Γ∗ ]B2 is a multiplication operator with norm bounded by tk∇ηk∞ , together with the off-diagonal bounds for QB t . For 2 the second term, observe that, since ΠB Hodge decomposes L according to Proposition 2.16, we have that B B B ktΓ∗B ηQB t uk2 . ktΠB ηQt uk2 ≤ k[ηI, tΠB ]Qt uk2 + kηtΠB Qt uk2 .

Here, we use that the commutator in the first part of the sum is again a multiplication B 2 2 operator. For the second part, we use that tΠB QB t = I − Pt , which satisfies L -L off-diagonal bounds.  7. Lp -L2 off-diagonal bounds We assume n ≥ 3 throughout this section. In this section, we prove Theorem 3.9 in the case p > pH . In a first step, we use an induction argument to establish bisectoriality in Lp , as well as Lp -L2 off-diagonal bounds.

28

DOROTHEE FREY, ALAN MCINTOSH, AND PIERRE PORTAL

We then apply the main result of Section 10, which gives a comparison of conical and vertical square functions under Lp -L2 off-diagonal bounds, and reduces the high freqency estimate to a vertical square function estimate on the unperturbed operator. The assertion of Theorem 3.9 then follows from the boundedness of the H ∞ functional calculus of the unperturbed operator. The first lemma shows how to deduce Lp∗ -Lp bounds from Lp bisectoriality via a Sobolev inequality, and serves as an induction step in the following. Lemma 7.1. Suppose n ≥ 3. Suppose ΠB is a perturbed Hodge-Dirac operator. Suppose p ∈ (1, 2] with p∗ > 1, and assume that p > pH and that ΠB is bisectorial in Lp (Rn ; CN ). Then sup ktRtB ukp . kukp∗

(7.1)

∀u ∈ R2 (Γ) ∩ Lp∗ (Rn ; CN ).

t∈R

Moreover, if (ΠB (p∗ )) holds, then sup ktRtB ukp . kukp∗

(7.2)

t∈R

∀u ∈ R2 (Γ∗B ) ∩ Lp∗ (Rn ; CN ).

Proof. We first show (7.1) for u = Γv, and v ∈ S(Rn ; CN ). We use the assumptions and a Sobolev inequality to obtain ktRtB Γvkp = ktΓ∗B RtB vkp . kvkp . kvkW˙ 1,p∗ . Using [33, Proposition 5.2], which gives that kvkW˙ 1,p∗ . kΠvkp∗ for all v ∈ Dp∗ (Π) ∩ Rp∗ (Π), we then get ktRtB Γvkp . kΠvkp∗ = kΓvkp∗ . The estimate (7.1) thus follows by density of S(Rn ; CN ) in W 1,p∗ . In order to obtain the analogous estimate (7.2), we apply (7.1) to the related operator ΠB , which is bisectorial in Lp and Hodge decomposes Lp by Proposition 2.17. We first prove the result with RtB n N ∗ ∗ ∗ replaced by QB t . Let v ∈ S(R ; C ), and consider u = B1 Γ v. Observe that Γ B2 B1 Γ = 0 B ∗ B ∗ implies Pt B1 Γ v = B1 Pt Γ v. Using (ΠB (p∗ )), (ΠB (p)), and (7.1), we obtain ∗ kt2 ΓPtB B1 Γ∗ vkp = kt2 ΓB1 PtB Γ∗ vkp . kt2 B2 ΓB1 PtB Γ∗ vkp = ktQB t Γ vkp B . kt(RtB − R−t )Γ∗ vkp . kΓ∗ vkp∗ . kB1 Γ∗ vkp∗ .

We thus get ∗ 2 B ∗ ∗ ktQB t B1 Γ vkp = kt ΓPt B1 Γ vkp . kB1 Γ vkp∗ .

In the exact same way, we also get ktPtB B1 Γ∗ vkp . kB1 Γ∗ vkp∗ , and thus ktRtB B1 Γ∗ vkp . kB1 Γ∗ vkp∗ , since RtB = PtB − iQB t . The result follows by approximating functions of the form B2 w with w ∈ D2 (Γ∗B ) ∩ Lp∗ (Rn ; CN ) by functions in the Schwartz class.  ( p∗(k) = (p∗(k−1) )∗ ∀k ∈ N, We use the following induction argument in which and p∗(0) = p, Ms (p) is the smallest natural number such that p∗(Ms (p)) ≥ 2. A simple induction argument np for all M ∈ N, so that Ms (p) ≥ n( p1 − 21 ). gives p∗(M ) = n−pM

CONICAL SQUARE FUNCTION ESTIMATES FOR DIRAC OPERATORS

29

Proposition 7.2. Suppose n ≥ 3. Suppose ΠB is a perturbed Hodge-Dirac operator. Suppose p ∈ (pH , 2] with p∗ > 1. Assume that ΠB is bisectorial in Lp (Rn ; CN ). Assume further that for all M ∈ N such that M ≥ Ms (p) and all r ∈ (p, 2] (with r = 2 if p = 2), {(RtB )M PRr (ΠB ) ; t ∈ R∗ } has Lr -L2 off-diagonal bounds of every order. (1) Given q ∈ (p∗ , 2] and M ∈ N such that M ≥ Ms (p), we have 1

1

sup ktn( q − 2 ) (RtB )M uk2 . kukq

∀u ∈ R2 (Γ) ∩ Lq (Rn ; CN ).

t∈R

Moreover, assuming the (ΠB (p∗ )) holds if q < pH , we have that 1

1

sup ktn( q − 2 ) (RtB )M uk2 . kukq t∈R

∀u ∈ R2 (Γ∗B ) ∩ Lq (Rn ; CN ).

(2) Given q ∈ (max{pH , p∗ }, 2] and M ∈ N such that M ≥ Ms (p), we have 1

1

sup ktn( q − 2 ) (RtB )M PRq (ΠB ) uk2 . kukq

∀u ∈ Lq (Rn ; CN ),

t∈R

and {(RtB )M PRq (ΠB ) ; t ∈ R∗ } has Lq -L2 off-diagonal bounds of every order. (3) For all q ∈ (max{pH , p∗ }, 2], ΠB is bisectorial in Lq (Rn ; CN ). Proof. (1) We have for all r ∈ (p, 2] and all M ∈ N such that M ≥ Ms (p) 1

1

sup ktn( r − 2 ) (RtB )M PRr (ΠB ) uk2 . kukr

∀u ∈ Lr (Rn ; CN )

t∈R

as a consequence of the Lr -L2 off-diagonal bounds for (RtB )M (take E = F = Rn in (2.1)). Combining this with Lemma 7.1 gives the assertion for q = r∗ . (2) For q > pH , ΠB Hodge decomposes Lq (Rn ; CN ) by assumption. We therefore get the first estimate in (2) as a direct consequence of (1). By Lemma 2.5 and Theorem 2.20, this implies that, for all q˜ ∈ (q, 2], {(RtB )M PRq˜(ΠB ) ; t ∈ R∗ } has Lq˜-L2 off-diagonal bounds of every order. (3) Interpolating between the uniform bounds from Lemma 7.1 and the L2 −L2 off-diagonal bounds as in Lemma 2.5 (1), and using Lemma 2.5 (2), gives that {RtB PRq (ΠB ) ; t ∈ R∗ } is uniformly bounded in Lq (Rn CN ). The result follows, since we have that RtB = RtB PRq (ΠB ) + RtB PNq (ΠB ) = RtB PRq (ΠB ) + PNq (ΠB ) for all t ∈ R∗ and all q > pH .  Corollary 7.3. Suppose n ≥ 3. Suppose ΠB is a perturbed Hodge-Dirac operator. Then, given q ∈ (pH , 2], the assertions (2) and (3) of Proposition 7.2 are satisfied for every M ∈ N such that M ≥ Ms (q). Moreover, if (ΠB ((pH )∗ )) holds, then given q ∈ (max{1, (pH )∗ }, 2], assertion (1) is satisfied for every M ∈ N such that M ≥ Ms (q). Proof. We deduce the statement inductively from Proposition 7.2. For p = 2, the assumptions of Proposition 7.2 are satisfied according to Theorem 2.20 (note that 2∗ > 1 for n ≥ 3). This yields (2) and (3) for all q ∈ (max{pH , 2∗ }, 2], and (1) for all q ∈ (2∗ , 2]. We now apply Proposition 7.2 inductively with exponent p. In each step, if p > pH and p∗ > 1, we know from the previous step that the assumptions of the proposition are satisfied. We run the induction until either pH > p or p∗ ≤ 1. In both cases, the claim of the corollary then follows from the penultimate induction step. The process is finite.  We now prove Theorem 3.9 in the case p > pH , which we restate here.

30

DOROTHEE FREY, ALAN MCINTOSH, AND PIERRE PORTAL

Corollary 7.4. Suppose n ≥ 3. Suppose ΠB is a perturbed Hodge-Dirac operator. Suppose M ∈ N with M ≥ 10n and p ∈ (pH , 2]. Then ˜

N M k(t, x) 7→ (QB t ) (I − Pt )u(x)kT p,2 ≤ Cp kukp

∀u ∈ R2 (Γ) ∩ Lp (Rn ; CN ).

Proof. Let p ∈ (pH , 2]. For u ∈ R2 (Γ) ∩ Lp , write ˜

N M B 2 B M −2 ˜ (QB Qt u, t ) (I − Pt )u = tΓ(Qt ) (Qt ) ˜ −1 N 2 2 2 ˜ t = ( P Ptk )Qt . Theorem 2.20 and Lemma 2.5 yield that tΓ(QB with Q t ) satisfies L -L k=0

off-diagonal bounds of every order. According to Corollary 7.3, ΠB is bisectorial in Lp , and {(RtB )M −2 PRp (ΠB ) ; t ∈ R∗ } satisfies Lp -L2 off-diagonal bounds of every order. Thus, Lemma 2.8, Lemma 2.5 and Theorem 10.1 yield for all u ∈ R2 (Γ) ∩ Lp (Rn ; CN ) ˜ N M B M −2 ˜ k(t, x) 7→ (QB Qt u(x)kT p,2 t ) (I − Pt )u(x)kT p,2 . k(t, x) 7→ (Qt ) ˜ t u(x)kLp (L2 ) . . k(t, x) 7→ Q



Now Proposition 2.19 gives the result. 8. High frequency estimates for p ∈ (max{1, (pH )∗ }, 2]

In this section, we prove Theorem 3.9 in the case p ∈ (max{1, (pH )∗ }, pH ]. We remark that our proof also provides an alternative (though less direct) approach to Corollary 7.4 in the case when pH < p ≤ 2. The idea of the proof is to show, using Corollary 9.2, that the integral operator defined by ˆ ∞ ds TK (F )(t, . ) := K(t, s)F (s, . ) s 0 ˜

˜

N 2N M with K(t, s) := (QB and M sufficiently large, extends to a bounded t ) (I − Pt )Qs operator on tent spaces. The square function estimate of Theorem 3.9 is then reduced to the square function estimate for the unperturbed operator shown in Proposition 2.19.

Theorem 8.1. Suppose ΠB is a perturbed Hodge-Dirac operator. Let M be even and such that M ≥ 10n, and let p ∈ (1, 2). Moreover, assume that ΠB 2 is sectorial in Lp (Rn ; CN ) M and that {(PtB ) 2 −1 PRp (ΠB 2 ) ; t ∈ R∗ } has Lp -L2 off-diagonal bounds of every order. Then, for all q ∈ (max{1, p∗ }, 2], ˜

N M k(t, x) 7→ (QB t ) (I − Pt )u(x)kT q,2 ≤ Cq kukq

∀u ∈ R2 (Γ) ∩ Lq (Rn ; CN ).

Proof. Let p ∈ (1, ∞). Let u ∈ Lp . By Theorem 2.20 and Lemma 2.8 we have that ˆ∞ ˆ∞ ds ds ˜ ˜ ˜ N 2N M k(t, x) 7→ (QB u(x) kT p,2 . k(t, x) 7→ Qs 2N u(x) kT p,2 t ) (I − Pt )Qs s s t

t

= k(t, x) 7→ ψ(tΠ)u(x)kT p,2 ,

CONICAL SQUARE FUNCTION ESTIMATES FOR DIRAC OPERATORS

31

´∞ zs 2N˜ ds for ψ(z) = ( 1+(zs) . Therefore 2) s 1

ˆ∞ ds ˜ ˜ N 2N M k(t, x) 7→ (QB u(x) kT p,2 . kukp , t ) (I − Pt )Qs s t

by Theorem 2.10 and Proposition 2.19. Now let q ∈ (max{1, p∗ }, 2]. For u ∈ R2 (Γ) ∩ Lq (Rn ; CN ), we have that ˆt

˜

˜

N 2N M u(x) (QB t ) (I − Pt )Qs

k(t, x) 7→

ds kT q,2 s

0

ˆt = k(t, x) 7→

ds s ˜ ˜ −1 ˜ N N M Qs N u(x) kT q,2 ( )tΓ(QB t ) (I − Pt )Ps Qs t s

0

ˆt . k(t, x) 7→

s ds ˜ ˜ ˜ M −2 ( )(QB (I − Pt N )Ps Qs N −1 Qs N u(x) kT q,2 , t ) t s

0

where in the last step we have used Lemma 6.2 and Lemma 2.8. We now consider the integral operator TK with kernel ˜

˜

M −2 K(t, s) = 1(0,∞) (t − s)(QB (I − Pt N )Ps Qs N −1 . t )

Using the results of Section 9, we aim to show that TK1+ extends to a bounded operator on T q,2 . The result then follows from Proposition 2.19. M From our assumption, we have that {(PtB ) 2 −1 PRp (ΠB 2 ) ; t ∈ R∗ } has Lp -L2 off-diagonal M

M

M −2 = (I − PtB ) 2 −1 .(PtB ) 2 −1 and {PtB ; t ∈ R} has bounds of every order. Since (QB t ) 2 2 M −2 L -L off-diagonal bounds of every order, we have that {(QB ; t ∈ R∗ } has Lp -L2 t ) off-diagonal bounds of every order. Using the bisectoriality of the unperturbed operator in Lq (see Proposition 2.19) and Lemma 2.5, we thus get that K satisfies (9.1) with max{t, s} = t for all r ∈ (p, 2]. To conclude the proof using Corollary 9.2, we thus only have to show that + supkTKε+iγ kL(T r,2 ) < ∞ ∀ε > 0 ∀r ∈ (p, 2].

γ∈R

To do so we use Lemma 2.8, Lemma 2.6, and Theorem 10.1, and obtain the following, for ε > 0, r ∈ (p, 2], F ∈ T r,2 , and γ ∈ R (with implicit constants independent of F and γ): ˆt + kTKε+iγ F kT r,2 . k(t, x) 7→

s ds ˜ ˜ ( )ε+iγ (I − Pt N )Ps Qs N −1 F (s, x) kLr (Rn ;L2 ((0,∞), dt )) t t s

0

= kTKe + F kLr (Rn ;L2 ((0,∞), dt )) , ε+iγ

˜

t

˜

e s) = (I − Pt N )Ps Qs N −1 . Since the unperturbed operator Π has a bounded where K(t, ˜ s) ; t, s > 0} is R-bounded in Lr by H ∞ functional calculus in Lr , the family {K(t,

32

DOROTHEE FREY, ALAN MCINTOSH, AND PIERRE PORTAL

[36, Theorem 5.3]. Therefore, Lemma 9.3 gives + F kT r,2 . kF kLr (Rn ;L2 ((0,∞), dt )) . kTKε+iγ t

We conclude the proof using [10, Proposition 2.1] to get that kF kLr (Rn ;L2 ((0,∞), dt )) . t kF kT r,2 .  B Combining this theorem with Corollary 7.3, and noting that PtB = RtB R−t , we have the following:

Corollary 8.2. Suppose ΠB is a perturbed Hodge-Dirac operator. Suppose p ∈ (max{1, (pH )∗ }, 2] and M ∈ N with M ≥ 10n. Then ˜

N M k(t, x) 7→ (QB t ) (I − Pt )u(x)kT p,2 ≤ Cp kukp

∀u ∈ R2 (Γ) ∩ Lp (Rn ; CN ).

If we restrict the off-diagonal bound assumptions to certain subspaces, the following restricted version of the theorem remains valid. Corollary 8.3. Suppose ΠB is a perturbed Hodge-Dirac operator. Let W be a subspace c∗ (ξ)Γ(ξ) b b Γ c∗ (ξ) for all ξ ∈ Rn . Let p ∈ (1, 2), M ∈ N of CN that is stable under Γ and Γ(ξ) be even and such that M ≥ 10n, and assume that ΠB 2 is sectorial on Lp (Rn ; W ) and M {(PtB ) 2 −1 PRp (ΠB 2 ) ; t ∈ R∗ } has Lp (Rn ; W ) − L2 (Rn ; CN ) off-diagonal bounds of every order. Then, for all q ∈ (max{1, p∗ }, 2], ˜

N M k(t, x) 7→ (QB t ) (I − Pt )Γv(x)kT q,2 ≤ Cq kΓvkq

∀v ∈ Dq (Γ) ∩ Lq (Rn ; W ).

Proof. Let p ∈ (1, 2), v ∈ S(Rn ; W ), and s > 0. Notice that, for all ξ ∈ Rn , −1 c2 (ξ))−1 Γ(ξ)b b b v (ξ) = sΓ c∗ (ξ)Γ(ξ)(I b c∗ (ξ)Γ(ξ)) b sΠ(ξ)(I + s2 Π + s2 Γ vb(ξ) ∈ W,

b c∗ (ξ)Γ(ξ). b since Γ(ξ) is nilpotent, and W is stable under Γ We thus have that Qs Γv belongs ˜ ˜ p n to L (R ; W ). The same reasoning also gives that (I − Pt N )Ps Q2sN −1 Γv ∈ Lp (Rn ; W ) for all t, s > 0. Therefore, denoting by PW the projection from Lp (Rn ; CN ) onto Lp (Rn ; W ), we have that, for all t > 0 and x ∈ Rn , ˆt

s ds ˜ ˜ ˜ M −2 ( )(QB (I − Pt N )Ps Qs N −1 Qs N Γv(x) t ) t s

0

ˆt =

s ds ˜ ˜ ˜ M −2 ( )(QB PW (I − Pt N )Ps Qs N −1 Qs N Γv(x) . t ) t s

0

This allows us to use the proof of Theorem 8.1. The only change required is to replace ˜ ˜ M −2 e s) = 1(0,∞) (t − s)(QB the kernel K by K(t, PW (I − Pt N )Ps Qs N −1 , for all t, s > 0, t ) M and remark that {(PtB ) 2 −1 PRp (ΠB 2 ) PW ; t > 0} has Lp (Rn ; CN )-L2 (Rn ; CN ) off-diagonal bounds of every order. 

CONICAL SQUARE FUNCTION ESTIMATES FOR DIRAC OPERATORS

33

9. Technical tools 1: Schur estimates in tent spaces We need boundedness criteria for integral operators of the form ˆ∞ ds N TK F (t, x) = K(t, s)F (s, .)(x) ∀F ∈ T p,2 (Rn+1 + ; C ), s 0

where {K(t, s) ; t, s > 0} is a uniformly bounded family of bounded linear operators acting on L2 (Rn ; CN ). We are interested here in Schur type estimates, i.e. estimates for integral operators with kernels satisfying size conditions of the form kK(t, s)k . min( st , st )α for some α > 0. The proofs are similar to those developed in [11] to treat singular integral operators with kernels satisfying size conditions of the form kK(t, s)k . |t − s|−1 . The appropriate off-diagonal bound assumptions are as follows. Let p ∈ [1, 2]. Let {K(t, s), s, t > 0} be a uniformly bounded family of bounded linear operators acting on L2 (Rn ; CN ) that satisfies Lp -L2 off-diagonal bounds of the following form: there exists C > 0, N 0 > 0, such that for all Borel sets E, F ⊆ Rn and all s, t > 0 −N 0  dist(E, F ) −n( p1 − 12 ) (9.1) k1E K(t, s)1F kLp →L2 ≤ C max{t, s} . 1+ max{t, s} Given a kernel K, we also consider s Kz+ (t, s) = 1(0,∞) (t − s)( )z K(t, s), t t ∀t, s ∈ (0, ∞), ∀z ∈ C. Kz− (t, s) = 1(0,∞) (s − t)( )z K(t, s), s N We then obtain the following result on T 1,2 (Rn+1 + ; C ), which is a refined version of the arguments in [15, Theorem 4.9]. Proposition 9.1. Suppose K satisfies (9.1) for some N 0 > n2 and p ∈ [1, 2]. Then the following holds. N (1) Given α ∈ (0, ∞), we have TKα− ∈ L(T 1,2 (Rn+1 + ; C )). N (2) Given β ∈ ( pn0 , ∞) and γ ∈ R, we have TK + ∈ L(T 1,2 (Rn+1 + ; C )) with β+iγ

sup kTK + kL(T 1,2 ) < ∞. γ∈R

β+iγ

Corollary 9.2. Suppose K satisfies (9.1) for some N 0 > n2 and p ∈ [1, 2]. Suppose + kL(T p,2 ) < ∞, then T − q ∈ [1, p], α ∈ (0, ∞), and β ∈ (n( 1q − p1 ), ∞). If sup kTKiγ Kα +K + ∈ β

γ∈R

N L(T q,2 (Rn+1 + ; C )).

Proof of Corollary 9.2. This follows from Proposition 9.1 by applying Stein’s interpolation [42, Theorem 1] to the analytic family of operators {TK −n ; Re(z) ∈ [0, 1]}. We choose z p0

p,2

1,2

1 q

the spaces T and T as endpoints, and set = 1 0 1 θ = p ( q − p ) and thus the condition β > pn0 θ = n( 1q −

1−θ p 1 ). p

+ θ for θ ∈ [0, 1]. This gives 

We now turn to the proof of Proposition 9.1, which follows the one of [15, Theorem 4.9].

34

DOROTHEE FREY, ALAN MCINTOSH, AND PIERRE PORTAL

Proof of Proposition 9.1. Let α > 0, β >

n . p0

It suffices to show that

kTKα− +K + F kT 1,2 ≤ C β+iγ

uniformly for all atoms F in T 1,2 and all γ ∈ R. Let F be a T 1,2 atom associated with a ball B ⊆ Rn of radius r > 0. Then ¨ |F (s, x)|2 T (B)

dxds ≤ |B|−1 , s

˜ = Kα− + K + , F˜ := T ˜ (F ), and F˜1 := F˜ 1T (4B) , F˜k := where T (B) = (0, r) × B. Set K K β+iγ ˜ F 1T (2k+1 B)\T (2k B) , k ≥ 2. We show that there exists δ > 0, independent of γ, such that ¨ |F˜k (t, x)|2

dxdt . 2−kδ |2k+1 B|−1 . t

ε ´∞ Let k = 1. Observe that for every ε > 0, 0 min st , st ds ≤ C, uniformly in t > 0. Using s Minkowski’s inequality and the assumption on K, we obtain ¨

dxdt |F˜ (t, x)|2 ≤ t T (4B)

ˆ

l(4B)

ˆ

ˆ



( 0

ˆ

4B

l(4B)

ˆ (

0

s t ds dxdt min( , )min(α,β) |K(t, s)F (s, . )(x)| )2 t s s t



s t ds dt min( , )min(α,β) kK(t, s)F (s, . )kL2 (4B) )2 t s s t ˆ0 ∞ ˆ ∞0 s t ds dt . ( min( , )min(α,β) kF (s, . )kL2 (B) )2 t s s t ˆ0 ∞ 0 ds . kF (s, . )k2L2 (B) ≤ |B|−1 . s 0 ≤

For k ≥ 2, we cover T (2k+1 B) \ T (2k B) with the two parts (0, 2k−1 r) × 2k+1 B \ 2k−1 B and (2k−1 r, 2k+1 r) × 2k+1 B. Via Minkowski’s inequality, we have ¨ (

|F˜k (t, x)|2

T (2k+1 B)\T (2k B)

ˆ

r



ˆ (

0

ˆ

+ 0

2k−1 r

0

r

ˆ (

s t dt 1 ds min( , )min(α,β) kK(t, s)F (s, . )k2L2 (2k+1 B\2k−1 B) ) 2 t s t s

2k+1 r

2k−1 r

=: I1 + I2 .

dxdt 1 )2 t

s dt 1 ds ( )β kK(t, s)F (s, . )k2L2 (2k+1 B) ) 2 t t s

CONICAL SQUARE FUNCTION ESTIMATES FOR DIRAC OPERATORS

35

For I2 , the fact that s < t and the assumed Lp -L2 boundedness of K(t, s) yield ˆ

r

ˆ (

r

ˆ (

I2 ≤ 0

ˆ

2k+1 r

2k−1 r 2k+1 r

s dt 1 ds ( )β kK(t, s)F (s, . )k2L2 (Rn ) ) 2 t t s

1 1 s dt 1 ds (( )β t−n( p − 2 ) kF (s, . )kLp (B) )2 ) 2 t t s k−1 ˆ0 r 2 r 1 1 1 1 s ds ( k )β (2k r)−n( p − 2 ) rn( p − 2 ) kF (s, . )kL2 (B) . s 0 2 r ˆ r ˆ r 1 1 ds 1 s ds 1 kF (s, . )k2L2 (B) ) 2 . 2−k(β+n( p − 2 )) ( ( )2β ) 2 ( s s 0 r 0

.

1

1

n

1

. 2−k(β+n( p − 2 )− 2 ) |2k B|− 2 .

Since, by assumption, β > n(1 − p1 ), this yields the desired estimate for I2 . We split the term I1 into the two parts ˆ

ˆ s t dt 1 ds I1 ≤ ( ( )α kK(t, s)F (s, . )k2L2 (2k+1 B\2k−1 B) ) 2 t s 0 0 s ˆ r ˆ 2k−1 r s dt 1 ds + ( ( )β kK(t, s)F (s, . )k2L2 (2k+1 B\2k−1 B) ) 2 t t s 0 s =: I1,1 + I1,2 . r

For I1,1 , we have t < s. The assumed Lp -L2 off-diagonal bounds and the fact that N 0 > n( p1 − 21 ) yield ˆ

ˆ

r

s

1 1 t dist(B, 2k−1 B) −N 0 dt 1 ds (( )α s−n( p − 2 ) (1 + ) kF (s, . )kLp (B) )2 ) 2 s s t s ˆ0 r 0 ˆ s 1 1 s s t dt 1 ds 0 . kF (s, . )kL2 (B) ( )−n( p − 2 ) ( k )N ( ( )2α ) 2 r 2 r t s 0 0 s ˆ r ˆ r 1 1 ds 1 ds 1 s 0 0 . 2−kN ( kF (s, . )k2L2 (B) ) 2 ( ( )2N −2n( p − 2 ) ) 2 s s 0 0 r

I1,1 .

(

0

1

0

n

1

. 2−kN |B|− 2 . 2−k(N − 2 ) |2k B|− 2 .

Since N 0 > n2 , this yields the assertion for I1,1 . ˜ > 0 with For I1,2 , we have s < t. According to our assumptions, there exists N

n 2

˜<
36

DOROTHEE FREY, ALAN MCINTOSH, AND PIERRE PORTAL

max(N 0 , β + n( p1 − 12 )). Using the Lp -L2 off-diagonal bounds, we get ˆ r ˆ 2k−1 r 1 1 dt 1 ds s 2k r −N 0 I1,2 . ( (( )β t−n( p − 2 ) (1 + ) kF (s, . )kLp (B) )2 ) 2 t t t s 0 s ˆ r ˆ r 1 1 dt 1 ds 0 s t 0 kF (s, . )kL2 (B) ( ( )2β ( )2N −2n( p − 2 ) ) 2 . 2−kN t r t s s 0 ˆ 2k−1 r ˆ r 1 1 dt 1 ds s t ˜ ˜ ( )2β ( )2N −2n( p − 2 ) ) 2 kF (s, . )kL2 (B) ( + 2−kN t r t s r ˆ r0 s ds ˜ ( )β kF (s, . )kL2 (B) . 2−kN r s ˆ r ˆ0 r 1 ds 1 ds 1 s ˜ n ˜ kF (s, . )k2L2 (B) ) 2 . 2−k(N − 2 ) |2k B|− 2 , . 2−kN ( ( )2β ) 2 ( s s 0 0 r where again we use the assumptions N 0 >

n 2

and β > n(1 − p1 ).



We conclude this section by pointing out that such estimates are much simpler in the context of vertical, rather than conical, square functions. In particular we have the following lemma (see [41, Section 5] for the relevant information regarding R-boundedness). Lemma 9.3. Suppose p ∈ (1, ∞) and ε > 0. If {K(t, s) ; t, s > 0} is a R-bounded family of bounded operators on Lp (Rn ), then TKε+ +Kε− ∈ L(Lp (Rn ; L2 (R+ ; dtt ))). Proof. Let F ∈ Lp (Rn ; L2 (R+ ; dtt )). By Kalton-Weis’ γ-multiplier theorem (see [41, Theorem 5.2]), we have the following: ˆ1 sε k(x, t) 7→ K(t, ts)F (x, ts)kLp (Rn ;L2 (R+ ; dt ))

kTKε+ F kLp (Rn ;L2 (R+ ; dt )) ≤ t

t

ds s

0

ˆ1 sε k(x, t) 7→ F (x, ts)kLp (Rn ;L2 (R+ ; dt ))

.

t

ds s

0

ˆ1 sε k(x, t) 7→ F (x, t)kLp (Rn ;L2 (R+ ; dt ))

=

t

ds . kF kLp (Rn ;L2 (R+ ; dt )) . t s

0

The same reasoning applies to TKε− .



10. Technical tools 2: Estimating conical square functions by vertical square functions While vertical and conical square functions look similar, the conical square functions are applied quite differently here compared with the way the vertical square functions are used in [32]. Nevertheless, as is the case classically (see e.g. [43]), there are relationships between conical and vertical square functions, as Auscher, Hofmann, and Martell have already pointed out in [10]. Here we prove a new comparison theorem that exploits Lp -L2 off-diagonal bounds, and H ∞ functional calculus. The proof is based on some unpublished work of Auscher, Duong and the second author, where a similar result was

CONICAL SQUARE FUNCTION ESTIMATES FOR DIRAC OPERATORS

37

obtained for operators with pointwise Gaussian bounds. Recall that a sectorial operator D2 in Lp (Rn ; CN ) induces a decomposition Lp (Rn ; CN ) = Np (D2 ) ⊕ Rp (D2 ). We denote by PRp (D2 ) the corresponding projection onto Rp (D2 ). Theorem 10.1. Suppose p ∈ (1, 2] and M ∈ N even with M > n + 4. Suppose D is a bisectorial operator in L2 (Rn ; CN ) with a bounded H ∞ functional calculus in L2 (Rn ; CN ), and suppose that {(I + itD)−1 ; t ∈ R} has L2 -L2 off-diagonal bounds of every order. M Assume further that D2 is sectorial in Lp (Rn ; CN ), and that {(I +t2 D2 )−( 2 −1) PRp (D2 ) ; t ∈ R∗ } has Lp -L2 off-diagonal bounds of every order. Then k(t, x) 7→ (tD(I + t2 D2 )−1 )M F (t, .)(x)kT p,2

ˆ∞ dt 1 ≤ Cp k( |F (t, .)|2 ) 2 kp t 0

∀F ∈ Cc ((0, ∞) × Rn ; CN ). Proof. We use a variant of the Blunck-Kunstmann extrapolation method established by Auscher in [3, Theorem 1.1]. As pointed out in [3] the extrapolation theorem holds in Hilbert space valued Lp spaces. Let us therefore consider H1 = L2 (R+ , dtt ; CN ), H2 = N L2 (R+ × Rn , tdtdx n+1 ; C ), and an operator T defined by T (F ) : x 7→ [(t, y) 7→ 1B(y,t) (x)(tD(I + t2 D2 )−1 )M PRp (D2 ) F (t, .)(y)], for F ∈ L2 (Rn ; H1 ). Since D has a bounded H ∞ functional calculus in L2 , T is a bounded operator from L2 (Rn ; H1 ) to L2 (Rn ; H2 ). Define ψt (D) = (tD(I + t2 D2 )−1 )M ,

φt (D) = I − (I − (I + t2 D2 )−M )M

∀t > 0.

Notice that by assumption, both {ψt (D) ; t > 0} and {φt (D)PRp (D2 ) ; t ∈ R∗ } have Lp -L2 M

off-diagonal bounds of every order. This follows from the fact that {(I+t2 D2 )− 2 PRp (D2 ) ; t ∈ R∗ } has Lp -L2 off-diagonal bounds of every order, and {(I + t2 D2 )−1 ; t ∈ R∗ } has L2 -L2 off-diagonal bounds of every order. To check the two off-diagonal bound conditions of [3, Theorem 1.1], let B ⊂ Rn be a ball of radius r > 0, F ∈ L2 (Rn ; H1 ) be supported in B, and j ∈ N. The two conditions are, in our setting, ˆ ˆ 1 1 1 2 −1 2 (10.1) ( j+1 kF (., y)kpH1 dy) p , kφr (D)PRp (D2 ) F (., y)kH1 dy) . g(j)(|B| |2 B| B

Sj (B)

for j ≥ 1, and (10.2) (

1

ˆ kT (I −

|2j+1 B| Sj (B)

1 φr (D)PRp (D2 ) )F (., y)k2H2 dy) 2

ˆ

1

kF (., y)kpH1 dy) p

−1

. g(j)(|B|

B

38

DOROTHEE FREY, ALAN MCINTOSH, AND PIERRE PORTAL

P for j ≥ 2, with g(j) satisfying j∈N g(j)2nj < ∞. Using the Lp -L2 off-diagonal bounds of {φt (D)PRp (D2 ) ; t ∈ R∗ }, we have that (

1

ˆ

1 kφr (D)PRp (D2 ) F (., y)k2H1 dy) 2

|2j+1 B|

−j n −n 2 2

≈2

r

ˆ∞ ˆ (

|φr (D)PRp (D2 ) F (t, y)|2

dydt 1 )2 t

0 Sj (B)

Sj (B) −j(M + n ) −n p 2

.2

r

ˆ∞ (

dt 1 kF (t, .)k2p ) 2 t

−j(M + n ) 2

.2

− p1

|B|

0 −j(M + n ) 2

.2

−1

(|B|

ˆ

1

|F (., y)|p dyk p2

k

p

B

ˆ ˆ∞ ˆ 1 1 n 2 dt p2 −j(M + ) −1 2 (|B| ( |F (t, y)| ) dy) p = 2 kF (., y)kpH1 dy) p . t 0

B

B

This establishes (10.1). The proof will thus be complete once we have established (10.2). To do so, we first use the straightforward integration lemma [22, Lemma 1], and obtain that for j ≥ 2 ˆ 1 1 kT (I − φr (D)PRp (D2 ) )F (., y)k2H2 dy) 2 ( j+1 |2 B| Sj (B)

.(

ˆ∞

1

ˆ

|2j+1 B|

|ψt (D)(I − φr (D))F (t, y)|2

dydt 1 )2 t

|ψt (D)(I − φr (D))F (t, y)|2

dydt 1 )2 t

0 R(Sj (B))

≤(

ˆ∞

1

ˆ

|2j+1 B| 0 (2j−1 B)c

+

j−2 X

(

k=0

where R(Sj (B)) =

ˆ∞

1

ˆ

2 dydt

|ψt (D)(I − φr (D))F (t, y)|

|2j+1 B|

t

1 2

) =: I +

Ik ,

k=0

(2j−1 −2k )r Sk (B)

S

j−2 X

{(t, y) ∈ R+ × Rn ; |y − x| ≤ t}. Let k ∈ {0, ..., j − 2}, and let

x∈Sj (B)

us first estimate Ik . Notice that r ˜ ψt (D)(I − φr (D)) = ( )M (t2 D2 )M (I + t2 D2 )−M φ(rD) ∀t, r > 0, t   2 )M −1 M ˜ for φ(z) = (1+z , and that {(t2 D2 )M (I + t2 D2 )−M ; t > 0} has L2 -L2 offz(1+z 2 )M ˜ diagonal bounds. In order to show that {φ(rD) ; r > 0} has Lp -L2 off-diagonal bounds M M ˜ ˜ of order N 0 > n2 , write φ(z) = (1 + z 2 ) 2 −1 φ(z)(1 + z 2 )−( 2 −1) . By Lemma 2.6, {(1 + M ˜ r2 D2 ) 2 −1 φ(rD) ; r > 0} has L2 -L2 off-diagonal estimates of order N 0 > n . On the M

2

other hand, {(I + r2 D2 )−( 2 −1) PRp (D2 ) ; r ∈ R∗ } has Lp -L2 off-diagonal estimates of every order. Combining the two families of operators gives the statement (see Lemma 2.5). 1 1 ˜ Thus {( rt )n( p − 2 ) (t2 D2 )M (I + t2 D2 )−M φ(rD) ; t > 0} has Lp -L2 off-diagonal bounds of n 0 order N > 2 , again by Lemma 2.5, whenever 0 < r < t. In particular, kψt (D)(I −

CONICAL SQUARE FUNCTION ESTIMATES FOR DIRAC OPERATORS

39

1 1 ˜ ˜ := M − n( 1 − 1 ). This gives φr (D))kL(Lp ,L2 ) . ( rt )M t−n( p − 2 ) , for M p 2

Ik = (

ˆ∞

1

ˆ |ψt (D)(I − φr (D))F (t, y)|2

|2j+1 B|

dydt 1 )2 t

(2j−1 −2k )r Sk (B)

.(

ˆ∞

1

1 1 r ˜ dt 1 (( )M t−n( p − 2 ) kF (t, .)kp )2 ) 2 t t

|2j+1 B| (2j−1 −2k )r

ˆ∞

1 1 n r dt 1 ˜ ((2j r)− 2 ( )M t−n( p − 2 ) kF (t, .)kp )2 ) 2 . t t

.( (2j−1 −2k )r

Taking into account that 2j−1 − 2k ≥ 2j−2 , we can estimate the above by ∞ 1 ´ ´ 1 n ´ n 1 ˜ n ˜ n ˜ n 1 kF (., y)kpH1 dy) p . 2−j(M + p ) r− p ( kF (t, .)k2p dtt ) 2 . 2−j(M + p ) r− p k |F (., y)|p dyk p2 . 2−j(M + p ) ( |B| 2j r

p

B

Summing up over k then yields

j−2 P

˜ +n) −j(M p

Ik . j2

k=0

I ≤ J1 + J2 := (

ˆr

1

1 ( |B|

´

B

1 p

kF (., y)kpH1 dy) . We now estimate

B

ˆ |ψt (D)(I − φr (D))F (t, y)|2

|2j+1 B|

dydt 1 )2 t

0 (2j−1 B)c

+(

1

ˆ∞

ˆ |ψt (D)(I − φr (D))F (t, y)|2

|2j+1 B|

dydt 1 )2 . t

r (2j−1 B)c

For J1 , we use that for 0 < t < r, {ψt (D) ; t > 0} and {ψt (D)φr (D) ; r > 0} have Lp -L2 off-diagonal bounds of every order. Thus, (

ˆr

1

ˆ |ψt (D)(I − φr (D))F (t, y)|2

|2j+1 B|

dydt 1 )2 t

0 (2j−1 B)c

≤(

1

ˆr

ˆ

|ψt (D)F (t, y)|

|2j+1 B| 0 (2j−1 B)c

.(

1

2 dydt

ˆr

−n +n p 2

(t

|2j+1 B| 0

) +(

1

ˆr

ˆ |ψt (D)φr (D)F (t, y)|2

|2j+1 B|

B

dydt 1 )2 t

0 (2j−1 B)c

2j r −N 0 dt 1 1 (1 + ) kF (t, .)kp )2 ) 2 + ( j+1 t t |2 B|

ˆr

n

n

0

(r− p + 2 2−jN kF (t, .)kp )2 0

ˆr t N 0 − np + n2 dt 1 0) − n 2 dt 21 −j( n +N .2 r ( ( ) kF (t, .)kp ) ) + 2 2 r p ( kF (t, .)kp )2 ) 2 r t t 0 0 ˆ 1 n 1 0 . 2−j( 2 +N ) ( kF (., y)kpH1 dy) p . |B| −j( n +N 0 ) − n p 2

ˆr

t

1 2

dt 1 )2 t

40

DOROTHEE FREY, ALAN MCINTOSH, AND PIERRE PORTAL 1

1

˜ Turning to J2 , we now use that {( rt )n( p − 2 ) (t2 D2 )M (I + t2 D2 )−M φ(rD) ; t > 0} has Lp -L2 n 0 off-diagonal bounds of order N > 2 , which gives (

ˆ∞

1

ˆ |ψt (D)(I − φr (D))F (t, y)|2

|2j+1 B|

dydt 1 )2 t

r (2j−1 B)c

ˆ

2j r

.( r

1 1 n r dt 1 2j r −N 0 ˜ ((2j r)− 2 ( )M t−n( p − 2 ) ( ) kF (t, .)kp )2 ) 2 t t t

ˆ∞ 1 1 n r dt 1 ˜ + ( ((2j r)− 2 ( )M t−n( p − 2 ) kF (t, .)kp )2 ) 2 t t 2j r

ˆ ˆ2j r 1 r M˜ 2j r −(N 0 −n( p1 − 21 )) ˜ +n) 1 −j(M 2 dt 12 p ) kF (., y)kpH1 dy) p .2 r ( (( ) ( kF (t, .)kp ) ) + 2 ( t t t |B| r B ˆ 1 1 0 n ˜ n kF (., y)kpH1 dy) p . . (2−j(N + 2 ) + 2−j(M + p ) )( |B| −j n p

−n p

B

Combining all the estimates gives ˆ ˆ 1 1 1 2 −j min{M + n ,N 0 + n } −1 2 2 2 ( j+1 kT (I−φr (D))F (., y)kH2 dy) . j2 (|B| kF (., y)kpH1 dy) p , |2 B| B

Sj (B)

which shows (10.2), given that M > 1.1], and conclude the proof.

n 2

and N 0 > n2 . This allows us to apply [3, Theorem 

As a Corollary, we obtain a slight generalisation of a result of P. Auscher [3, Corollary 6.10] (see also M. Uhl’s thesis [44, Theorem 4.19]). This result is interesting in view of the following observation. Given p ∈ (1, ∞), a bisectorial operator D has a bounded H ∞ functional calculus in Lp (Rn ; CN ) if and only if k(t, x) 7→ tD(I + t2 D2 )−1 ukLp (Rn ;L2 (R+ , dt )) ≈ kukp t

∀u ∈ Rp (D).

On the other hand, for a bisectorial operator D with a bounded H ∞ functional calculus in L2 (Rn ; CN ), that is bisectorial in Lp , and is such that {(I + itD)−1 ; t ∈ R} has L2 -L2 off-diagonal bounds of order N 0 > n2 , the conical square function estimate k(t, x) 7→ tD(I + t2 D2 )−M ukT p,2 ≈ kukp

∀u ∈ Rp (D)

for M > n2 implies that D has a bounded H ∞ functional calculus in Lp (Rn ; CN ). Indeed, bisectoriality gives the decomposition Lp (Rn ; CN ) = Np (D) ⊕ Rp (D), and D has a p functional calculus on HD (Rn ; CN ) = Rp (D) by Theorem 2.10. The corollary below thus provides a partial converse to this result. We do not know, however, if the Lp -L2 offdiagonal bounds assumption can be weakened to L2 -L2 off-diagonal bounds, which would give the full converse. As shown in Proposition 4.1, the converse is known for specific

CONICAL SQUARE FUNCTION ESTIMATES FOR DIRAC OPERATORS

41

p operators (note that p ∈ (pH , pH ) is equivalent to HD (Rn ; CN ) = Rp (D), since D always p n N Hodge decomposes HD (R ; C )).

Corollary 10.2. Let D be a bisectorial operator with a bounded H ∞ functional calculus in L2 (Rn ; CN ), and suppose that {(I + itD)−1 ; t ∈ R} has L2 -L2 off-diagonal bounds of every order. Let p ∈ (1, 2] and M ∈ N even with M > n + 4. Assume further that D is a bisectorial operator with a bounded H ∞ functional calculus in Lp (Rn ; CN ), and that M {(I + t2 D2 )−( 2 −1) PRp (D2 ) ; t ∈ R∗ } has Lp -L2 off-diagonal bounds of every order. Then k(t, x) 7→ (tD(I + t2 D2 )−1 )M +1 ukT p,2 ≈ kukp

∀u ∈ Rp (D).

p In particular, HD (Rn ; CN ) = Rp (D).

Proof. Let u ∈ Lp (Rn ; CN ). Applying Theorem 10.1 with F (t, x) = tD(I + t2 D2 )−1 u(x) for all (t, x) ∈ (0, ∞) × Rn , we get that ˆ∞ dt 1 2 2 −1 M +1 k(t, x) 7→ (tD(I + t D ) ) ukT p,2 . k( |tD(I + t2 D2 )−1 u|2 ) 2 kp . t 0

Since D has bounded H



p

functional calculus in L (Rn ; CN ), this gives

k(t, x) 7→ (tD(I + t2 D2 )−1 )M +1 ukT p,2 . kukp . The reverse inequality for u ∈ Rp (D) is proven exactly as in Proposition 3.6, using the L2 -L2 off-diagonal bounds, and the fact that D has a bounded H ∞ functional calculus in L2 (Rn ; CN ).  References [1] S. Ajiev. Extrapolation of the functional calculus of generalized Dirac operators and related embedding and Littlewood-Paley-type theorems. J. Aust. Math. Soc. 83 (2007) 297–326. [2] D. Albrecht, X. Duong, A. McIntosh. Operator theory and harmonic analysis. In: Instructional Workshop on Analysis and Geometry, Part III (Canberra, 1995), Proc. Centre Math. Appl. Austral. Nat. Univ. 34 (1996) 77–136. [3] P. Auscher. On necessary and sufficient conditions for Lp -estimates of Riesz transforms associated to elliptic operators on Rn and related estimates. Mem. Amer. Math. Soc. 871 (2007). [4] P. Auscher. Change of angle in tent spaces. C. R. Math. Acad. Sci. Paris 349 (2011) 297–301. [5] P. Auscher, A. Axelsson. Weighted maximal regularity estimates and solvability of non-smooth elliptic systems I. Invent. Math. 184 (2011) 47–115. [6] P. Auscher, A. Axelsson, S. Hofmann. Functional calculus of Dirac operators and complex perturbations of Neumann and Dirichlet problems. J. Funct. Anal. 255 (2008) 374–448. [7] P. Auscher, A. Axelsson, A. McIntosh. Solvability of elliptic systems with square integrable boundary data. Ark. Mat. 48 (2010) 253–287. [8] P. Auscher, A. Axelsson, A. McIntosh. On a quadratic estimate related to the Kato conjecture and boundary value problems. Harmonic analysis and partial differential equations, 105-129, Contemp. Math., 505, Amer. Math. Soc., Providence, RI (2010). [9] P. Auscher, S. Hofmann, M. Lacey, A. McIntosh, P. Tchamitchian. The solution of the Kato square root problem for second order elliptic operators on Rn . Ann. of Math. 156 (2002) 633–654. [10] P. Auscher, S. Hofmann, J.M. Martell. Vertical versus conical square functions. Trans. Amer. Math. Soc. 364 (2012) 5469–5489. [11] P. Auscher, C. Kriegler, S. Monniaux, P. Portal. Singular integral operators on tent spaces. J. Evol. Equ. 12 (2012) 741–765.

42

DOROTHEE FREY, ALAN MCINTOSH, AND PIERRE PORTAL

[12] P. Auscher, J.M. Martell. Weighted norm inequalities, off-diagonal estimates and elliptic operators Part II: Off-diagonal estimates on spaces of homogeneous type. J. Evol. Equ. 7 (2007) 265–316. [13] P. Auscher, A. McIntosh, A. Morris. Calderón reproducing formulas and applications to Hardy spaces. Preprint arXiv:1304.0168. [14] P. Auscher, A. McIntosh, A. Nahmod. Holomorphic functional calculi of operators, quadratic estimates and interpolation. Indiana Univ. Math. J. 46 (1997) 375–403. [15] P. Auscher, A. McIntosh, E. Russ. Hardy Spaces of Differential Forms on Riemannian Manifolds. J. Geom. Anal. 18 (2008) 192–248. [16] P. Auscher, S. Stahlhut. Remarks on functional calculus for perturbed first order Dirac operator. To appear in Operator Theory in Harmonic and Non Commutative Analysis, Oper. Theory Adv. Appl. 240 Birkhäuser/Springer Basel AG. [17] P. Auscher, S. Stahlhut. A priori estimates for boundary value elliptic problems via first order systems. arXiv preprint arXiv:1403.5367. [18] A. Axelsson, S. Keith, A. McIntosh. Quadratic estimates and functional calculi of perturbed Dirac operators. Invent. math. 163 (2006) 455–497. [19] L. Bandara, A. McIntosh. The Kato square root problem on vector bundles with generalised bounded geometry. arXiv preprint arXiv:1203.0373. [20] S. Blunck, P. Kunstmann. Calderón-Zygmund theory for non-integral operators and the H ∞ functional calculus. Rev. Mat. Iberoamericana. 19 (2003) 919–942. [21] W. Cohn, I. Verbitsky. Factorization of tent spaces and Hankel operators. J. Funct. Anal. 175 (2000) 308–329. [22] R. Coifman, Y. Meyer, E. Stein. Some new function spaces and their applications to harmonic analysis. J. Funct. Anal. 62 (1985) 304–335. [23] R. Coifman, G. Weiss. Analyse harmonique sur certains espaces homogènes. Lecture Notes in Mathematics. 242. Springer-Verlag (1971). [24] M. Cowling, I. Doust, A. McIntosh, A. Yagi. Banach space operators with a bounded H ∞ functional calculus. J. Austral. Math. Soc. Ser. A 60 (1996) 51–89. [25] X.T. Duong, L. Yan. Duality of Hardy and BMO spaces associated with operators with heat kernel bounds. J. Amer. Math. Soc. 18 (2005) 943–973. [26] C. Fefferman, E. Stein. H p spaces of several variables. Acta Math. 129 (1972) 137–193. [27] M. Haase. The functional calculus for sectorial operators. Birkhäuser Verlag, Basel (2006). [28] S. Hofmann, G. Lu, D. Mitrea, M. Mitrea, L. Yan. Hardy spaces associated to non-negative selfadjoint operators satisfying Davies-Gaffney estimates. Mem. Amer. Math. Soc. 214 (2011). [29] S. Hofmann, S. Mayboroda. Hardy and BM O spaces associated to divergence form elliptic operators. Math. Ann. 344 (2009) 37–116. [30] S. Hofmann, S. Mayboroda, A. McIntosh. Second order elliptic operators with complex bounded measurable coefficients in Lp , Sobolev and Hardy spaces. Ann. Sci. Éc. Norm. Supér. (4) 44 (2011) 723–800. [31] T. Hytönen, A. McIntosh. Stability in p of the H ∞ -calculus of first-order systems in Lp . Proc. Centre Math. Appl. Austral. Nat. Univ. 44 (2010) 167–181. [32] T. Hytönen, A. McIntosh, P. Portal. Kato’s square root problem in Banach spaces. J. Funct. Anal. 254 (2008) 675–726. [33] T. Hytönen, A. McIntosh, P. Portal. Holomorphic functional calculus of Hodge-Dirac operators in Lp . J. Evol. Equ. 11 (2011) 71–105. [34] T. Hytönen, J. van Neerven, P. Portal. Conical square function estimates in U M D Banach spaces and applications to H ∞ -functional calculi. J. Anal. Math. 106 (2008) 317–351. [35] N. Kalton, M. Mitrea. Stability results on interpolation scales of quasi-Banach spaces and applications. Trans. Amer. Math. Soc. 350 (1998) 3903–3922. [36] N. J. Kalton, L. Weis. The H ∞ -calculus and sums of closed operators. Math. Ann. 321 (2001) 319– 345. [37] P. C. Kunstmann, L. Weis. Maximal Lp regularity for parabolic problems, Fourier multiplier theorems and H ∞ -functional calculus, Lect. Notes in Math. 1855. Springer-Verlag (2004).

CONICAL SQUARE FUNCTION ESTIMATES FOR DIRAC OPERATORS

43

[38] A. McIntosh. Operators which have an H ∞ functional calculus. Proc. Centre Math. Appl. Austral. Nat. Univ. 14 (1986) 210–231. [39] A. McIntosh, A. Nahmod. Heat kernel estimates and functional calculi of −b∆. Math. Scand. 87, no. 2 (2000) 287–319. [40] A.J. Morris. The Kato square root problem on submanifolds. J. Lond. Math. Soc. (2) 86 (2012) 879–910. [41] J. van Neerven. γ-radonifying operators–a survey. Proc. Centre Math. Appl. Austral. Nat. Univ. 44 (2010) 1–61. [42] E. Stein. Interpolation of linear operators. Trans. Amer. Math. Soc. 83 (1956) 482–492. [43] E. M. Stein. Singular Integrals and Differentiability Properties of Functions. Princeton University Press, Princeton, NJ, 1970. [44] M. Uhl. Spectral multiplier theorems of Hörmander type via generalized Gaussian estimates. PhDThesis, Karlsruher Institut für Technologie. http://digbib.ubka.uni-karlsruhe.de/volltexte/1000025107 (2011).

Dorothee Frey Australian National University, Mathematical Sciences Institute, John Dedman Building, Acton ACT 0200, Australia. [email protected] Alan McIntosh Australian National University, Mathematical Sciences Institute, John Dedman Building, Acton ACT 0200, Australia. [email protected] Pierre Portal Australian National University, Mathematical Sciences Institute, John Dedman Building, Acton ACT 0200, Australia. [email protected]

PerturbedDirac-final-arxiv-v2.pdf

is one of Coifman-Meyer-Stein's tent spaces. (see [22] and Section 2 for precise definitions). This fact has been noticed in the de- velopment of a Hardy space ...

641KB Sizes 1 Downloads 142 Views

Recommend Documents

No documents