Singular integral operators on closed curves explained

In mathematics, singular integral operators on closed curves arise in problems in analysis, in particular complex analysis and harmonic analysis. The two main singular integral operators, the Hilbert transform and the Cauchy transform, can be defined for any smooth Jordan curve in the complex plane and are related by a simple algebraic formula. In the special case of Fourier series for the unit circle, the operators become the classical Cauchy transform, the orthogonal projection onto Hardy space, and the Hilbert transform a real orthogonal linear complex structure. In general the Cauchy transform is a non-self-adjoint idempotent and the Hilbert transform a non-orthogonal complex structure. The range of the Cauchy transform is the Hardy space of the bounded region enclosed by the Jordan curve. The theory for the original curve can be deduced from that of the unit circle, where, because of rotational symmetry, both operators are classical singular integral operators of convolution type. The Hilbert transform satisfies the jump relations of Plemelj and Sokhotski, which express the original function as the difference between the boundary values of holomorphic functions on the region and its complement. Singular integral operators have been studied on various classes of functions, including Hölder spaces, L^p spaces and Sobolev spaces. In the case of L² spaces—the case treated in detail below—other operators associated with the closed curve, such as the Szegő projection onto Hardy space and the Neumann–Poincaré operator, can be expressed in terms of the Cauchy transform and its adjoint.

Operators on the unit circle

See main article: Singular integral operators of convolution type. If f is in L²(T), then it has a Fourier series expansion

\displaystyle{f(\theta) = \sum _{n\in} } a_n e^.}

Hardy space H²(T) consists of the functions for which the negative coefficients vanish, a_n = 0 for n < 0. These are precisely the square-integrable functions that arise as boundary values of holomorphic functions in the unit disk |z| < 1. Indeed, f is the boundary value of the function

\displaystyle{F(z) = \sum _{n\ge} a _{n} z ^{n,}}

in the sense that the functions

\displaystyle{f ^{i\theta} _{r(\theta) = F(re} )},

defined by the restriction of F to the concentric circles |z| = r, satisfy

\displaystyle{\|f _{r-f\|} _{2} → 0 } as \displaystyle{r → 1 } .

The orthogonal projection P of L²(T) onto H²(T) is called the Szegő projection. It is a bounded operator on L²(T) with operator norm 1.

By Cauchy's theorem

\displaystyle{F(z) = { 1 \over 2 \pi i} \int _{|\zeta|=1} {f(\zeta)\over \zeta - z} d\zeta = { 1 \over 2 \pi} \int ^{\pi} _{- \pi} {f(\theta) \over 1- e ^{-i\theta} z} d\theta.}

Thus

\displaystyle{F(re ^{i\varphi} ) = { 1 \over 2 \pi} \int ^{\pi} _{- \pi} {f(\varphi - \theta) \over 1- re ^{i\theta} } \, d\theta.}

When r equals 1, the integrand on the right hand side has a singularity at θ = 0. The truncated Hilbert transform is defined by

\displaystyle{H ^{\varepsilon} f(\varphi) = {i\over \pi} \int _{\varepsilon\le} {f(\varphi - \theta) \over 1- e ^{i\theta} } \, d\theta= \int_

\ge \delta

\, d\zeta,}

where δ = |1 – e^iε|. Since it is defined as convolution with a bounded function, it is a bounded operator on L²(T). Now

\displaystyle{H ^{\pi} _{\varepsilon} 2 \Re ( 1- e ^{i\theta} ) ^{-1} d\theta = {i\over\pi}\int ^{\pi} _{\varepsilon} 1 d\theta = i - {i \varepsilon\over \pi}.}

If f is a polynomial in z then

\displaystyle{H ^{\varepsilon} f(z) - {i( 1- \varepsilon)\over \pi} f(z) = { 1 \over \pi i} \int _{|\zeta} {f(\zeta) - f(z)\over \zeta - z} d\zeta.}

By Cauchy's theorem the right hand side tends to 0 uniformly as ε, and hence δ, tends to 0. So

\displaystyle{H ^{\varepsilon} f → if}

uniformly for polynomials. On the other hand, if u(z) = z it is immediate that

\displaystyle{\overline{H ^{\varepsilon} f} = - u ^{-1} H ^{\varepsilon(} u \overline{f}).}

Thus if f is a polynomial in z⁻¹ without constant term

\displaystyle{H ^{\varepsilon} f → - i f} uniformly.

Define the Hilbert transform on the circle by

\displaystyle{H = i( 2 P - I).}

Thus if f is a trigonometric polynomial

\displaystyle{H ^{\varepsilon} f → Hf} uniformly.

It follows that if f is any L² function

\displaystyle{H ^{\varepsilon} f → Hf} in the L² norm.

This is a consequence of the result for trigonometric polynomials since the H_ε are uniformly bounded in operator norm: indeed their Fourier coefficients are uniformly bounded.

It also follows that, for a continuous function f on the circle, H_εf converges uniformly to Hf, so in particular pointwise. The pointwise limit is a Cauchy principal value, written

\displaystyle{Hf = P.V. { 1 \over \pi} \int {f(\zeta)\over \zeta - e ^{i\varphi} }\, d\zeta.}

The Hilbert transform has a natural compatibility with orientation-preserving diffeomorphisms of the circle.^[1] Thus if H is a diffeomorphism of the circle with

\displaystyle{H(e ^{i\theta} ) = e ^{ih(\theta)} , h(\theta +2 \pi) = h(\theta) +2 \pi,}

then the operators

\displaystyle{H ^{\varepsilon} _{h} f(e ^{i\varphi} ) = { 1 \over \pi}\int _{|e ^{ih(\theta)} - e ^{ih(\varphi)} |\ge \varepsilon} {f(e ^{i\theta} )\over e ^{i\theta} - e ^{i\varphi} } \, e^\, d\theta,}

are uniformly bounded and tend in the strong operator topology to H. Moreover, if Vf(z) = f(H(z)), then VHV⁻¹ – H is an operator with smooth kernel, so a Hilbert–Schmidt operator.

Hardy spaces

See also: Hardy space. The Hardy space on the unit circle can be generalized to any multiply connected bounded domain Ω with smooth boundary ∂Ω. The Hardy space H²(∂Ω) can be defined in a number of equivalent ways. The simplest way to define it is as the closure in L²(∂Ω) of the space of holomorphic functions on Ω which extend continuously to smooth functions on the closure of Ω. As Walsh proved, in a result that was a precursor of Mergelyan's theorem, any holomorphic function on Ω that extends continuously to the closure can be approximated in the uniform norm by a rational function with poles in the complementary region Ω^c. If Ω is simply connected, then the rational function can be taken to be a polynomial. There is a counterpart of this theorem on the boundary, the Hartogs–Rosenthal theorem, which states that any continuous function ∂Ω can be approximated in the uniform norm by rational functions with poles in the complement of ∂Ω. It follows that for a simply connected domain when ∂Ω is a simple closed curve, H²(∂Ω) is just the closure of the polynomials; in general it is the closure of the space of rational functions with poles lying off ∂Ω.^[2]

On the unit circle an L² function f with Fourier series expansion

\displaystyle{f(e ^{i\theta} ) = \sum a _{n} e ^{in\theta} }

has a unique extension to a harmonic function in the unit disk given by the Poisson integral

\displaystyle{f(re ^{i\theta} ) = P _{r} f(e ^{i\theta} ) = \sum a _{n} r ^{|n|} e ^{in\theta} .}

In particular

\displaystyle{\|P ^{2= \sum} _{rf\|} |a ^{2} _{n|} r ^{2|n|} ,}

so that the norms increase to the value at r = 1, the norm of f. A similar in the complement of the unit disk where the harmonic extension is given by

\displaystyle{F ^{i\theta} _{R(e} ) = F(Re ^{i\theta} ) = \sum a _{n} R ^{-|n|} a ^{in\theta} _{ne} .}

In this case the norms increase from the value at R = ∞ to the norm of f, the value at R = 1.

A similar result holds for a harmonic function f on a simply connected region with smooth boundary provided the L² norms are taken over the level curves in a tubular neighbourhood of the boundary. Using vector notation v(t) = (x(t), y(t)) to parametrize the boundary curve by arc length, the following classical formulas hold:

• \displaystyle{ v ⋅ • v =1 , \ddot{ v }\cdot \dot = 0.}

Thus the unit tangent vector t(t) at t and oriented normal vector n(t) are given by

\displaystyle{ t = • v , n = ( - • y , • x ).}

The constant relating the acceleration vector to the normal vector is the curvature of the curve:

\displaystyle{\ddot{ v }=\kappa(t)\,\mathbf(t),\,\,\,\,\, \kappa(t) =\ddot\cdot \mathbf=\ddot\dot - \ddot\dot.}

There are two further formulas of Frenet:

• \displaystyle{ n =- \kappa t , \ddot{ n }= \dot\mathbf -\kappa^2 \mathbf.}

A tubular neighbourhood of the boundary is given by

\displaystyle{ v _{s(t)=v(t)+sn(t),}}

so that the level curves ∂Ω_s with s constant bound domains Ω_s. Moreover

• \displaystyle{ v _{s(t)} = ( 1- s\kappa) t ,\partial _{s} | • v _{s|=-\kappa.}}

Hence differentiating the integral means with respect to s, the derivative in the direction of the inward pointing normal, gives

\displaystyle{\partial _{s} \int _{\partial\Omega _{s}} |f| ^{2} =- \int _{\partial\Omega _{s}} (\partial _{n} f \overline{f} + f \overline{\partial _{n} f}) - \int _{\partial\Omega _{s}} \kappa( 1- \kappa s) ^{-1} |f| ^{2} = -2 \iint _{\Omega _{s}} |\nabla f| ^{2} - \int _{\partial\Omega _{s}} \kappa( 1- \kappa s) ^{-1} |f| ^{2,}}

using Green's theorem. Thus for s small

\displaystyle{\partial _{s} \|f| _{\partial \Omega _{s}} \| ^{2} \le M \|f| _{\partial \Omega _{s}} \| ^{2,}}

for some constant M independent of f. This implies that

\displaystyle{\partial _{s} e ^{-Ms} \|f| _{\partial \Omega _{s}} \| ^{2} \le 0 ,}

so that, on integrating this inequality, the norms are bounded near the boundary:

\displaystyle{\|f| _{\partial \Omega _{s}} \| \le e ^{Ms/2} \|f| _{\partial} \|.}

This inequality shows that a function in the L² Hardy space H²(Ω) leads, via the Cauchy integral operator C, to a holomorphic function on Ω satisfying the classical condition that the integral means

\displaystyle{ \int _{\partial\Omega _{s}} |f| ^{2}}

are bounded. Furthermore, the restrictions f_s of f to ∂Ω_s, which can be naturally identified with ∂Ω, tend in L² to the original function in Hardy space. In fact H²(Ω) has been defined as the closure in L²(Ω) of rational functions (which can be taken to be polynomials if Ω is simply connected). Any rational function with poles only in Ω^c can be recovered inside Ω from its boundary value g by Cauchy's integral formula

\displaystyle{Cg(a) = { 1 \over 2 \pi i}\int _{\partial\Omega} {g(z)\over z - a} dz.}

The estimates above show that the functions Cg|_∂Ωs depend continuously on Cg|_∂Ω. Moreover, in this case the functions tend uniformly to the boundary value and hence also in L², using the natural identification of the spaces L²(∂Ω_s) with L²(∂Ω). Since Ch can be defined for any L² function as a holomorphic function on Ω since h is integrable on ∂Ω. Since h is a limit in L² of rational functions g, the same results hold for h and Ch, with the same inequalities for the integral means. Equally well h is the limit in L²(∂Ω) of the functions Ch|_∂Ωs.

The estimates above for the integral means near the boundary show that Cf lies in L²(Ω) and that its L² norm can be bounded in terms of that of f. Since Cf is also holomorphic, it lies in the Bergman space A²(Ω) of Ω. Thus the Cauchy integral operator C defines a natural mapping from the Hardy space of the boundary into the Bergman space of the interior.

The Hardy space H²(Ω) has a natural partner, namely the closure in L²(∂Ω) of boundary values of rational functions vanishing at ∞ with poles only in Ω. Denoting this subspace by H²₊(∂Ω) to distinguish it from the original Hardy space, which will also denoted by H²₋(∂Ω), the same reasoning as above can be applied. When applied to a function h in H²₊(∂Ω), the Cauchy integral operator defines a holomorphic function F in Ω^c vanishing at ∞ such that near the boundary the restriction of F to the level curves, each identified with the boundary, tend in L² to h. Unlike the case of the circle, H²₋(∂Ω) and H²₊(∂Ω) are not orthogonal spaces. By the Hartogs−Rosenthal theorem, their sum is dense in L²(∂Ω). As shown below, these are the ±i eigenspaces of the Hilbert transform on ∂Ω, so their sum is in fact direct and the whole of L²(∂Ω).

Hilbert transform on a closed curve

For a bounded simply connected domain Ω in the complex plane with smooth boundary ∂Ω, the theory of the Hilbert transform can be deduced by direct comparison with the Hilbert transform for the unit circle.^[3]

To define the Hilbert transform H_∂Ω on L²(∂Ω), take ∂Ω to be parametrized by arclength and thus a function z(t). The Hilbert transform is defined to be the limit in the strong operator topology of the truncated operators H_∂Ω^ε defined by

\displaystyle{H ^{\varepsilon} _{\partial\Omega} f(s) = { 1 \over \pi i} \int _{|s-t|\ge} {f(t)\over z(t) - z(s)} • z (t) dt.}

To make the comparison it will be convenient to apply a scaling transformation in C so that the length of ∂Ωis 2π. (This only changes the operators above by a fixed positive factor.) There is then a canonical unitary isomorphism of L²(∂Ω) onto L²(T), so the two spaces can be identified. The truncated operators H_∂Ω^ε can be compared directly with the truncated Hilbert transformH^ε:

\displaystyle{H ^{\varepsilon} _{\partial\Omega} g(s) - H ^{\varepsilon} g(s) = { 1 \over \pi i} \int K(s,t) ⋅ g(t) dt,}

where

• \displaystyle{K(u,v) = { z (t)\over z(t) - z(s)} - {ie ^{it} \over e ^{it} - e ^{is} }= \partial_t \log\left(\right).}

The kernel K is thus smooth on T × T, so the difference above tends in the strong topology to the Hilbert–Schmidt operator defined by the kernel. It follows that the truncated operators H_∂Ω^ε are uniformly bounded in norm and have a limit in the strong operator topology denoted H_∂Ω and called the Hilbert transform on ∂Ω.

\displaystyle{H _{\partial\Omega} g = \lim _{\varepsilon →} H ^{\varepsilon} _{\partial\Omega} g.}

Letting ε tend to 0 above yields

\displaystyle{H _{\partial\Omega} g(s) - H g(s) = { 1 \over \pi i} \int ^{2 \pi} _{0} K(s,t) ⋅ g(t) dt.}

Since H is skew-adjoint and H_∂Ω differs from H by a Hilbert–Schmidt operator with smooth kernel, it follows that H_∂Ω + H_∂Ω* is a Hilbert-Schmidt operator with smooth kernel. The kernel can also be computed explicitly using the truncated Hilbert transforms for ∂Ω:

\displaystyle{C(s,t) = { 1 \over \pi i}\left({ • z (t)\over z(t) - z(s)} - {\overline{ • z (s)}\over \overline{z(t)} - \overline{z(s)}} \right),}

and it can be verified directly that this is a smooth function on T × T.

Plemelj–Sokhotski relation

See also: Sokhotski–Plemelj theorem. Let C₋ and C₊ be the Cauchy integral operators for Ω and Ω^c. Then

\displaystyle{C _{\pm} \circ H = \pm i C _{\pm.}}

Since the operators C₋, C₊ and H are bounded, it suffices to check this on rational functions F with poles off ∂Ω and vanishing at ∞ by the Hartogs–Rosenthal theorem. The rational function can be written as a sum of functions F = F₋ + F₊ where F₋ has poles only in Ω^c and F₊ has poles only in Let f, f_± be the restrictions of f, f_± to ∂Ω. By Cauchy's integral formula

\displaystyle{C _{\pm} f _{\pm} = F _{\pm,} C _{\pm} f _{\mp} = 0 .}

On the other hand, it is straightforward to check that^[4]

\displaystyle{Hf _{\pm} = \pm i f _{\pm.}}

Indeed, by Cauchy's theorem, since F₋ is holomorphic in Ω,

\displaystyle{\int _{\partial\Omega,|z-w|\ge} {f _{-(z)\over} z - w} dz =- \int _{|z-w|=\varepsilon,z\in} {F _{-(z)\over} z - w} dz.}

As ε tends to 0, the latter integral tends to πi f₋(w) by the residue calculus. A similar argument applies to f₊, taking the circular contour on the right inside Ω^c.

By continuity it follows that H acts as multiplication by i on H²₋ and as multiplication by −i on H²₊. Since these spaces are closed and their sum dense, it follows that

\displaystyle{H ^{2=-I.}}

Moreover, H²₋ and H²₊ must be the ±i eigenspaces of H, so their sum is the whole of L²(∂Ω). The Plemelj–Sokhotski relation for f in L²(∂Ω) is the relation

\displaystyle{ - iHf = C _{-f|} _{\partial\Omega} - C _{+f|} _{\partial\Omega} .}

It has been verified for f in the Hardy spaces H²_±(∂Ω), so is true also for their sum. The Cauchy idempotent E is defined by

\displaystyle{H = i( 2 E - I).}

The range of E is thus H²₋(∂Ω) and that of I − E is H²₊(∂Ω). From the above

\displaystyle{Ef = C _{-f|} _{\partial\Omega} , (I - E)f = C _{+f|} _{\partial\Omega} .}

Operators on a closed curve

See also: Szegő kernel and Neumann–Poincaré operator. Two other operators defined on a closed curve ∂Ω can be expressed in terms of the Hilbert and Cauchy transforms H and E.^[5]

The Szegő projection P is defined to be the orthogonal projection onto Hardy space H²(∂Ω). Since E is an idempotent with range H²(∂Ω), P is given by the Kerzman–Stein formula:

\displaystyle{P = E(I + E - E ^{*)} ^{-1} .}

Indeed, since E − E* is skew-adjoint its spectrum is purely imaginary, so the operator I + E − E* is invertible. It is immediate that

\displaystyle{EP = P,PE = E.}

Hence PE* = P. So

\displaystyle{P(I + E - E ^{*)=} P + E - P = E.}

Since the operator H + H* is a Hilbert–Schmidt operator with smooth kernel, the same is true for E − E*.

Moreover, if J is the conjugate-linear operator of complex conjugation and U the operator of multiplication by the unit tangent vector:

\displaystyle{Jf(t) = \overline{f(t)}, Uf(t) = • z (t) ⋅ f(t).}

then the formula for the truncated Hilbert transform on ∂Ω immediately yields the following identity for adjoints

\displaystyle{(H ^{\varepsilon)} ^{*=JUH} ^{\varepsilon} U ^{*J.}}

Letting ε tend to 0, it follows that

\displaystyle{H ^{*=JUH} U ^{*J}}

and hence

\displaystyle{I - E ^{*=JUEU} ^{*J.}}

The comparison with the Hilbert transform for the circle shows that commutators of H and E with diffeomorphisms of the circle are Hilbert–Schmidt operators. Similar their commutators with the multiplication operator corresponding to a smooth function f on the circle is also Hilbert–Schmidt operators. Up to a constant the kernel of the commutator with H is given by the smooth function

\displaystyle{A(s,t) = {f(s) - f(t)\over z(s) - z(t)}.}

The Neumann–Poincaré operator T is defined on real functions f as

\displaystyle{Tf(w) = { 1 \over 2 \pi}\int _{\partial\Omega} \partial _{n} ( log |z - w|) f(z) = { 1 \over 2 }\Re (Hf)(w).}

Writing h = f + ig,

\displaystyle{ 2 Th = \Re(Hf) + i\Re(Hg) = { 1 \over 2 } (Hf + JHf + iHg + iJHg) = { 1 \over 2 }(H + JHJ)h}

so that

\displaystyle{T = { 1 \over 4 } (H + JHJ) = { 1 \over 4 }(H + UH ^{*U} ^{*),}}

a Hilbert–Schmidt operator.

Classical definition of Hardy space

The classical definition of Hardy space is as the space of holomorphic functions F on Ω for which the functions F_s = F|_∂Ωs have bounded norm in L²(∂Ω). An argument based on the Carathéodory kernel theorem shows that this condition is satisfied whenever there is a family of Jordan curves in Ω, eventually containing any compact subset in their interior, on which the integral means of F are bounded.

To prove that the classical definition of Hardy space gives the space H²(∂Ω), take F as above. Some subsequence h_n = F_sn converges weakly in L²(∂Ω) to h say. It follows that Ch = F in Ω. In fact, if C_n is the Cauchy integral operator corresponding to Ω_sn, then

\displaystyle{Ch(a) - F(a) = Ch(a) - C _{nh} _{n(a)=C(h-h} _{n)(a)} + [(C - C _{n)h} _{n](a).}}

Since the first term on the right hand side is defined by pairing h − h_n with a fixed L² function, it tends to zero. If z_n(t) is the complex number corresponding to v_sn, then

\displaystyle{[(C - C _{n)h} _{n](a)={1\over} 2 \pi i} \int h _{n(t)} \left({ • z (t) \over z(t) - a} - { • z _{n(t)} \over z _{n(t)-a}\right)} dt.}

This integral tends to zero because the L² norms of h_n are uniformly bounded while the bracketed expression in the integrand tends to 0 uniformly and hence in L².

Thus F = Ch. On the other hand, if E is the Cauchy idempotent with range H²(∂Ω), then C ∘ E = C. Hence F =Ch = C (Eh). As already shown F_s tends to Ch in L²(∂Ω). But a subsequence tends weakly to h. Hence Ch = h and therefore the two definitions are equivalent.

Generalizations

The theory for multiply connected bounded domains with smooth boundary follows easily from the simply connected case.^[6] There are analogues of the operators H, E and P. On a given component of the boundary, the singular contributions to H and E come from the singular integral on that boundary component, so the technical parts of the theory are direct consequences of the simply connected case.

Singular integral operators on spaces of Hölder continuous functions are discussed in . Their action on L^p and Sobolev spaces is discussed in .

Notes and References

1. See:
2. See:
3. See:
4. See:
5. See:
6. See:

This article is licensed under the GNU Free Documentation License. It uses material from the Wikipedia article "Singular integral operators on closed curves".

Except where otherwise indicated, Everything.Explained.Today is © Copyright 2009-2024, A B Cryer, All Rights Reserved. Cookie policy.