# Hardy space

In complex analysis, the **Hardy spaces** (or **Hardy classes**) *H ^{p}* are certain spaces of holomorphic functions on the unit disk or upper half plane. They were introduced by Frigyes Riesz (Riesz 1923), who named them after G. H. Hardy, because of the paper (Hardy 1915). In real analysis

**Hardy spaces**are certain spaces of distributions on the real line, which are (in the sense of distributions) boundary values of the holomorphic functions of the complex Hardy spaces, and are related to the

*L*spaces of functional analysis. For 1 ≤

^{p}*p*≤ ∞ these real Hardy spaces

*H*are certain subsets of

^{p}*L*, while for

^{p}*p*< 1 the

*L*spaces have some undesirable properties, and the Hardy spaces are much better behaved.

^{p}There are also higher-dimensional generalizations, consisting of certain holomorphic functions on tube domains in the complex case, or certain spaces of distributions on **R**^{n} in the real case.

Hardy spaces have a number of applications in mathematical analysis itself, as well as in control theory (such as *H*^{∞} methods) and in scattering theory.

For spaces of holomorphic functions on the open unit disk, the Hardy space *H*^{2} consists of the functions *f* whose mean square value on the circle of radius *r* remains bounded as *r* → 1 from below.

More generally, the Hardy space *H ^{p}* for 0 <

*p*< ∞ is the class of holomorphic functions

*f*on the open unit disk satisfying

The space *H*^{∞} is defined as the vector space of bounded holomorphic functions on the disk, with the norm

For 0 < p ≤ q ≤ ∞, the class *H ^{q}* is a subset of

*H*, and the

^{p}*H*-norm is increasing with

^{p}*p*(it is a consequence of Hölder's inequality that the

*L*-norm is increasing for probability measures, i.e. measures with total mass 1).

^{p}The Hardy spaces defined in the preceding section can also be viewed as certain closed vector subspaces of the complex *L ^{p}* spaces on the unit circle. This connection is provided by the following theorem (Katznelson 1976, Thm 3.8): Given

*f*∈

*H*, with

^{p}*p*≥ 0,

^{[clarification needed]}the radial limit

where the *ĝ*(*n*) are the Fourier coefficients of a function *g* integrable on the unit circle,

The space *H ^{p}*(

**T**) is a closed subspace of

*L*(

^{p}**T**). Since

*L*(

^{p}**T**) is a Banach space (for 1 ≤

*p*≤ ∞), so is

*H*(

^{p}**T**).

In applications, those functions with vanishing negative Fourier coefficients are commonly interpreted as the causal solutions.^{[clarification needed]} Thus, the space *H*^{2} is seen to sit naturally inside *L*^{2} space, and is represented by infinite sequences indexed by **N**; whereas *L*^{2} consists of bi-infinite sequences indexed by **Z**.

When 1 ≤ *p* < ∞, the *real Hardy spaces* *H ^{p}* discussed further down

^{[clarification needed]}in this article are easy to describe in the present context. A real function

*f*on the unit circle belongs to the real Hardy space

*H*(

^{p}**T**) if it is the real part of a function in

*H*(

^{p}**T**), and a complex function

*f*belongs to the real Hardy space iff Re(

*f*) and Im(

*f*) belong to the space (see the section on real Hardy spaces below). Thus for 1 ≤

*p*< ∞, the real Hardy space contains the Hardy space, but is much bigger, since no relationship is imposed between the real and imaginary part of the function.

For 0 < *p* < 1, such tools as Fourier coefficients, Poisson integral, conjugate function, are no longer valid. For example, consider the function

exists for a.e. *θ* and is in *H ^{p}*(

**T**), but Re(

*f*) is 0 almost everywhere, so it is no longer possible to recover

*F*from Re(

*f*). As a consequence of this example, one sees that for 0 <

*p*< 1, one cannot characterize the real-

*H*(

^{p}**T**) (defined below) in the simple way given above,

^{[clarification needed]}but must use the actual definition using maximal functions, which is given further along somewhere below.

For the same function *F*, let *f _{r}*(e

^{iθ}) =

*F*(

*re*

^{iθ}). The limit when

*r*→ 1 of Re(

*f*),

_{r}*in the sense of*

*distributions*on the circle, is a non-zero multiple of the Dirac distribution at

*z*= 1. The Dirac distribution at a point of the unit circle belongs to real-

*H*(

^{p}**T**) for every

*p*< 1 (see below).

For 0 < *p* ≤ ∞, every non-zero function *f* in *H ^{p}* can be written as the product

*f*=

*Gh*where

*G*is an

*outer function*and

*h*is an

*inner function*, as defined below (Rudin 1987, Thm 17.17). This "Beurling factorization" allows the Hardy space to be completely characterized by the spaces of inner and outer functions.

^{[1]}

^{[2]}

One says that *G*(*z*)^{[clarification needed]} is an **outer (exterior) function** if it takes the form

One says that *h* is an **inner (interior) function** if and only if |*h*| ≤ 1 on the unit disk and the limit

exists for almost all θ and its modulus is equal to 1 a.e. In particular, *h* is in *H*^{∞}.^{[clarification needed]} The inner function can be further factored into a form involving a Blaschke product.

The function *f*, decomposed as *f* = *Gh*,^{[clarification needed]} is in *H ^{p}* if and only if φ belongs to

*L*(

^{p}**T**), where φ is the positive function in the representation of the outer function

*G*.

Let *G* be an outer function represented as above from a function φ on the circle. Replacing φ by φ^{α}, α > 0, a family (*G*_{α}) of outer functions is obtained, with the properties:

It follows that whenever 0 < *p*, *q*, *r* < ∞ and 1/*r* = 1/*p* + 1/*q*, every function *f* in *H ^{r}* can be expressed as the product of a function in

*H*and a function in

^{p}*H*. For example: every function in

^{q}*H*

^{1}is the product of two functions in

*H*

^{2}; every function in

*H*,

^{p}*p*< 1, can be expressed as product of several functions in some

*H*,

^{q}*q*> 1.

Real-variable techniques, mainly associated to the study of *real Hardy spaces* defined on **R**^{n} (see below), are also used in the simpler framework of the circle. It is a common practice to allow for complex functions (or distributions) in these "real" spaces. The definition that follows does not distinguish between real or complex case.

Let *P _{r}* denote the Poisson kernel on the unit circle

**T**. For a distribution

*f*on the unit circle, set

where the *star* indicates convolution between the distribution *f* and the function e^{iθ} → *P _{r}*(θ) on the circle. Namely, (

*f*∗

*P*)(e

_{r}^{iθ}) is the result of the action of

*f*on the

*C*

^{∞}-function defined on the unit circle by

For 0 < *p* < ∞, the *real Hardy space* *H ^{p}*(

**T**) consists of distributions

*f*such that

*M f*is in

*L*(

^{p}**T**).

The function *F* defined on the unit disk by *F*(*re*^{iθ}) = (*f* ∗ *P _{r}*)(e

^{iθ}) is harmonic, and

*M f*is the

*radial maximal function*of

*F*. When

*M f*belongs to

*L*(

^{p}**T**) and

*p*≥ 1, the distribution

*f*"

*is*" a function in

*L*(

^{p}**T**), namely the boundary value of

*F*. For

*p*≥ 1, the

*real Hardy space*

*H*(

^{p}**T**) is a subset of

*L*(

^{p}**T**).

To every real trigonometric polynomial *u* on the unit circle, one associates the real *conjugate polynomial* *v* such that *u* + i*v* extends to a holomorphic function in the unit disk,

This mapping *u* → *v* extends to a bounded linear operator *H* on *L ^{p}*(

**T**), when 1 <

*p*< ∞ (up to a scalar multiple, it is the Hilbert transform on the unit circle), and

*H*also maps

*L*

^{1}(

**T**) to weak-

*L*

^{1}(

**T**). When 1 ≤

*p*< ∞, the following are equivalent for a

*real valued*integrable function

*f*on the unit circle:

When 1 < *p* < ∞, *H(f)* belongs to *L ^{p}*(

**T**) when

*f*∈

*L*(

^{p}**T**), hence the real Hardy space

*H*(

^{p}**T**) coincides with

*L*(

^{p}**T**) in this case. For

*p*= 1, the real Hardy space

*H*

^{1}(

**T**) is a proper subspace of

*L*

^{1}(

**T**).

The case of *p* = ∞ was excluded from the definition of real Hardy spaces, because the maximal function *M f* of an *L*^{∞} function is always bounded, and because it is not desirable that real-*H*^{∞} be equal to *L*^{∞}. However, the two following properties are equivalent for a real valued function *f*

When 0 < *p* < 1, a function *F* in *H ^{p}* cannot be reconstructed from the real part of its boundary limit

*function*on the circle, because of the lack of convexity of

*L*in this case. Convexity fails but a kind of "

^{p}*complex convexity*" remains, namely the fact that

*z*→ |

*z*|

^{q}is subharmonic for every

*q*> 0. As a consequence, if

is in *H ^{p}*, it can be shown that

*c*= O(

_{n}*n*

^{1/p–1}). It follows that the Fourier series

converges in the sense of distributions to a distribution *f* on the unit circle, and *F*(*re*^{iθ}) =(*f* ∗ *P _{r}*)(θ). The function

*F*∈

*H*can be reconstructed from the real distribution Re(

^{p}*f*) on the circle, because the Taylor coefficients

*c*of

_{n}*F*can be computed from the Fourier coefficients of Re(

*f*).

Distributions on the circle are general enough for handling Hardy spaces when *p* < 1. Distributions that are not functions do occur^{[where?]}, as is seen with functions *F*(*z*) = (1−*z*)^{−N} (for |*z*| < 1), that belong to *H ^{p}* when 0 <

*N*

*p*< 1 (and

*N*an integer ≥ 1).

A real distribution on the circle belongs to real-*H ^{p}*(

**T**) iff it is the boundary value of the real part of some

*F*∈

*H*. A Dirac distribution δ

^{p}_{x}, at any point

*x*of the unit circle, belongs to real-

*H*(

^{p}**T**) for every

*p*< 1; derivatives δ′

_{x}belong when

*p*< 1/2, second derivatives δ′′

_{x}when

*p*< 1/3, and so on.

It is possible to define Hardy spaces on other domains than the disc, and in many applications Hardy spaces on a complex half-plane (usually the right half-plane or upper half-plane) are used.

The Hardy space *H ^{p}*(

**H**) on the upper half-plane

**H**is defined to be the space of holomorphic functions

*f*on

**H**with bounded (quasi-)norm, the norm being given by

The corresponding *H*^{∞}(**H**) is defined as functions of bounded norm, with the norm given by

Although the unit disk **D** and the upper half-plane **H** can be mapped to one another by means of Möbius transformations, they are not interchangeable^{[clarification needed]} as domains for Hardy spaces. Contributing to this difference is the fact that the unit circle has finite (one-dimensional) Lebesgue measure while the real line does not. However, for *H*^{2}, one has the following theorem: if *m* : **D** → **H** denotes the Möbius transformation

In analysis on the real vector space **R**^{n}, the Hardy space^{[clarification needed]} *H ^{p}* (for 0 <

*p*≤ ∞) consists of tempered distributions

^{[clarification needed]}

*f*such that for some Schwartz function Φ with ∫Φ = 1, the maximal function

is in *L ^{p}*(

**R**

^{n}),

^{[clarification needed]}where ∗ is convolution and Φ

_{t }(

*x*) =

*t*

^{ −n}Φ(

*x*/

*t*). The

*H*-quasinorm ||

^{p}*f*||

_{Hp}of a distribution

*f*of

*H*is defined to be the

^{p}*L*norm of

^{p}*M*

_{Φ}

*f*(this depends on the choice of Φ, but different choices of Schwartz functions Φ give equivalent norms). The

*H*-quasinorm is a norm when

^{p}*p*≥ 1, but not when

*p*< 1.

If 1 < *p* < ∞, the Hardy space *H ^{p}* is the same vector space as

*L*, with equivalent norm. When

^{p}*p*= 1, the Hardy space

*H*

^{1}is a proper subspace of

*L*

^{1}. One can find sequences in

*H*

^{1}that are bounded in

*L*

^{1}but unbounded in

*H*

^{1}, for example on the line

The *L*^{1} and *H*^{1} norms are not equivalent on *H*^{1}, and *H*^{1} is not closed in *L*^{1}. The dual of *H*^{1} is the space *BMO* of functions of bounded mean oscillation. The space *BMO* contains unbounded functions (proving again that *H*^{1} is not closed in *L*^{1}).

If *p* < 1 then the Hardy space *H ^{p}* has elements that are not functions, and its dual

^{[clarification needed]}is the homogeneous Lipschitz space of order

*n*(1/

*p*− 1). When

*p*< 1, the

*H*-quasinorm is not a norm, as it is not subadditive. The

^{p}*p*th power ||

*f*||

_{Hp}

^{p}is subadditive for

*p*< 1 and so defines a metric on the Hardy space

*H*, which defines the topology and makes

^{p}*H*into a complete metric space.

^{p}When 0 < *p* ≤ 1, a bounded measurable function *f* of compact support is in the Hardy space *H ^{p}* if and only if all its moments

whose order *i*_{1}+ ... +*i _{n}* is at most

*n*(1/

*p*− 1), vanish. For example, the integral of

*f*must vanish in order that

*f*∈

*H*, 0 <

^{p}*p*≤ 1, and as long as

*p*>

*n*/ (

*n*+1) this is also sufficient.

If in addition *f* has support in some ball *B* and is bounded by |*B*|^{−1/p} then *f* is called an ** H^{p}-atom** (here |

*B*| denotes the Euclidean volume of

*B*in

**R**

^{n}). The

*H*-quasinorm of an arbitrary

^{p}*H*-atom is bounded by a constant depending only on

^{p}*p*and on the Schwartz function Φ.

When 0 < *p* ≤ 1, any element *f* of *H ^{p}* has an

**atomic decomposition**as a convergent infinite combination of

*H*-atoms,

^{p}On the line for example, the difference of Dirac distributions *f* = δ_{1}−δ_{0} can be represented as a series of Haar functions, convergent in *H ^{p}*-quasinorm when 1/2 <

*p*< 1 (on the circle, the corresponding representation is valid for 0 <

*p*< 1, but on the line, Haar functions do not belong to

*H*when

^{p}*p*≤ 1/2 because their maximal function is equivalent at infinity to

*a*

*x*

^{−2}for some

*a*≠ 0).

Let (*M _{n}*)

_{n≥0}be a martingale on some probability space (Ω, Σ,

*P*), with respect to an increasing sequence of σ-fields (Σ

_{n})

_{n≥0}. Assume for simplicity that Σ is equal to the σ-field generated by the sequence (Σ

_{n})

_{n≥0}. The

*maximal function*of the martingale is defined by

Let 1 ≤ *p* < ∞. The martingale (*M _{n}*)

_{n≥0}belongs to

*martingale*-

*H*when

^{p}*M**∈

*L*.

^{p}If *M** ∈ *L ^{p}*, the martingale (

*M*)

_{n}_{n≥0}is bounded in

*L*; hence it converges almost surely to some function

^{p}*f*by the martingale convergence theorem. Moreover,

*M*converges to

_{n}*f*in

*L*-norm by the dominated convergence theorem; hence

^{p}*M*can be expressed as conditional expectation of

_{n}*f*on Σ

_{n}. It is thus possible to identify martingale-

*H*with the subspace of

^{p}*L*(Ω, Σ,

^{p}*P*) consisting of those

*f*such that the martingale

Doob's maximal inequality implies that martingale-*H ^{p}* coincides with

*L*(Ω, Σ,

^{p}*P*) when 1 <

*p*< ∞. The interesting space is martingale-

*H*

^{1}, whose dual is martingale-BMO (Garsia 1973).

The Burkholder–Gundy inequalities (when *p* > 1) and the Burgess Davis inequality (when *p* = 1) relate the *L ^{p}*-norm of the maximal function to that of the

*square function*of the martingale

Martingales with continuous time parameter can also be considered. A direct link with the classical theory is obtained via the complex Brownian motion (*B _{t}*) in the complex plane, starting from the point

*z*= 0 at time

*t*= 0. Let τ denote the hitting time of the unit circle. For every holomorphic function

*F*in the unit disk,

is a martingale, that belongs to martingale-*H ^{p}* iff

*F*∈

*H*(Burkholder, Gundy & Silverstein 1971).

^{p}In this example, Ω = [0, 1] and Σ_{n} is the finite field generated by the dyadic partition of [0, 1] into 2^{n} intervals of length 2^{−n}, for every *n* ≥ 0. If a function *f* on [0, 1] is represented by its expansion on the Haar system (*h _{k}*)

then the martingale-*H*^{1} norm of *f* can be defined by the *L*^{1} norm of the square function

This space, sometimes denoted by *H*^{1}(δ), is isomorphic to the classical real *H*^{1} space on the circle (Müller 2005). The Haar system is an unconditional basis for *H*^{1}(δ).