Lecture 16

pp-periodic distributions

Let p>0p>0, a distribution TT is pp-periodic if τpT=T\tau_p T = T. Similar to the 1-periodic case, by choosing a partition of unity φn\varphi_n subordinate to the intervals [np/2,np/2][-np/2,np/2] we get

(FT,f)=(T,Ff)=(T,φ0nZFf(x+np)).(\mathcal FT,f) = (T,\mathcal Ff) = (T,\varphi_0 \sum_{n\in \mathbb{Z}}\mathcal Ff(x+np)).

This time nZFf(x+np)=(ωp(y),Ff(y+x))=(Fτxωp(y),f)=(e2πixy1pω1p(y),f(y))=1pnZf(np)e2πpinx\sum_{n\in \mathbb{Z}}\mathcal Ff(x+np) = (\omega_p(y),\mathcal Ff(y+x)) = (\mathcal F \tau_x\omega_p(y),f) =( e^{-2\pi i xy}\frac{1}{p}\omega_{\frac{1}{p}}(y),f(y)) = \frac{1}{p}\sum_{n\in \mathbb{Z}}f(\frac{n}{p})e^{-\frac{2\pi}{p}inx}. Of course, one can also derive this by applying the Poisson summation formula directly.

This implies FT=nZ(T(y),φ0e2πpiny)1pδnp\mathcal FT = \sum_{n\in \mathbb{Z}}(T(y),\varphi_0e^{-\frac{2\pi }{p}iny})\frac{1}{p}\delta_{\frac{n}{p}}. If we write T^(n)=1p(T(y),φ0e2πpiny)\hat{T}(n) = \frac{1}{p}(T(y),\varphi_0e^{-\frac{2\pi}{p}iny}), then by taking inverse transform (recall that F1δnp=e2πpinx\mathcal F^{-1}\delta_{\frac{n}{p}} = e^{\frac{2\pi}{p}inx}), we can write T=nZT^(n)e2πpinxT = \sum_{n\in \mathbb{Z}}\hat{T}(n)e^{\frac{2\pi }{p}inx}, this is the Fourier expansion of TT, which is parallel with the Fourier series theory.

Fundamental solution of D+sD+s

Fundamental solution of a linear differential operator

The reason why fundamental solution is useful is because the following lemma and the fact that δ is the unit of convolution.

Lemma

ddx(fT)=fT\frac{d}{dx}(f*T) = f*T'. Consequently, P(fT)=fPTP(f*T) = f*PT.

Proof. ((fT),g)=(fT,g)=(T,fg)=(T,(fg))=(T,fg)=(fT,g)((f*T)',g) = -(f*T,g') = -(T,f^-*g') = -(T,(f^-*g)')=(T',f^-*g) = (f*T',g). \square

It follows that if we can solve PE=δPE = \delta, then P(gE)=(gPE)=gδ=gP(g*E) = (g*PE) = g*\delta = g, in other words, gEg*E is solution of Pf=gPf = g in sense of distribution!

The case of D+sD+s

Let D:=iddxD:=-i\frac{d}{dx} be the Dirac operator and ss be a complex number.

In problem set 3 we gave solution of the equation (D+s)f=g(D+s)f = g by explicitedly calculating the inverse of D+sD+s using Fourier series and showed that the solution is given by f=gGf = g*G, where GG is given by Fourier series nZ1n+seinx\sum_{n\in \mathbb{Z}}\frac{1}{n+s}e^{inx}.

But this method may not work for general differential operators. In this note we’ll revisit the problem using common method of solving differential equations from distribution theory. From the derivative theorem we see that (D+s)G(D+s) G has Fourier coefficients 1, so as a 2π2\pi-periodic distribution, (D+s)G(D+s)G is equal to 2πω2π2\pi\omega_{2\pi}. From the above point of view, we can interprate GG as a fundamental solution of 2π2\pi-periodic functions since ω2π\omega_{2\pi} is the unit of convolution for 2π2\pi-periodic functions. (Recall that convolving with ω2π\omega_{2\pi} = periodization to period 2π2\pi, but if it’s already 2π2\pi-periodic then periodization does nothing). You may wonder why (D+s)=2πω2π(D+s)=2\pi\omega_{2\pi} instead of ω2π\omega_{2\pi}, this is because when we defined convolution at that time, we made a normalization by deviding 2π2\pi.

Use Fourier transform to solve D+sD+s

One of the applications of Fourier transform is that it can help us to find solution of partial differential equations in a direct way. As an example we consider the simple case D+sD+s. Later we’ll use the Fourier transform to solve the wave equation and heat equation.

TBA.

The Cauchy kernel: A complex analysis interlude

Let DD be an open set with C1\mathscr{C}^1 bundary on C\mathbb{C}. We adapt the identification C\mathbb{C} with R2\mathbb{R}^2 via zz with x+iyx+iy, where i=1i = \sqrt{-1}. Then for f(z)f(z) which is C1\mathscr{C}^1 on a neighbourhood of DD, we have the following Cauchy integral formula

The Cauchy integral formula

f(ζ)=12πi(Dzˉf(z)zζdzdzˉ+Df(z)zζdz)f(\zeta) = \frac{1}{2\pi i}(\int_{D}\frac{\partial_{\bar{z}}f(z)}{z-\zeta}dz\wedge d\bar{z}+\int_{\partial D}\frac{f(z)}{z-\zeta}dz).

In the real notation it is f(ζ)=12πiDf(z)zζdz1πDzˉf(z)zζdλ(z)f(\zeta) = \frac{1}{2\pi i }\int_{\partial D}\frac{f(z)}{z-\zeta}dz -\frac{1}{\pi}\int_{D}\frac{\partial_{\bar{z}}f(z)}{z-\zeta} d\lambda(z), where dλ(z)=dxdyd\lambda(z) = dxdy stands for the Lebesgue measure on C=R2\mathbb{C}=\mathbb{R}^2.

Proof. We give a proof of this using a distribution version of divergence theorem (or Green formula in 2 dimensional case), see Problem set 8 problem 3. The result is

We’ll use the complex version of this result. Choose a parametrization γ(t)=(x(t),y(t))\gamma(t) = (x(t),y(t)) of D\partial D, such that the tangent is given by (x(t),y(t))(x'(t),y'(t)) (or γ(t)=x(t)+iy(t)\gamma'(t) = x'(t)+iy'(t)) and invards normal is given by (y(t),x(t))(-y'(t),x'(t)) (or iγ(t)=i(x(t)+iy(t))=y(t)+ix(t)i\gamma'(t) = i(x'(t)+iy'(t)) = -y'(t)+ix'(t)). Then ds=γ(t)dtds =|\gamma'(t)|dt, n=1γ(t)(y(t),x(t))n = \frac{1}{|\gamma'(t)|}(-y'(t),x'(t)) so n1ds=y(t)dtn_1ds = -y'(t)dt and n2ds=x(t)dtn_2 ds = x'(t)dt. Then (x[χD],f)=Dfxdxdy=(n1ds,f)=f(x(t),y(t))y(t)dt=Dfdy(-\partial_x[\chi_D],f) =\int_{D}\frac{\partial f}{\partial x}dxdy = -(n_1ds,f) = \int f(x(t),y(t))y'(t)dt = \int_{\partial D}fdy. Similarly (y[χD],f)=Dfdx(-\partial_y[\chi_D],f) = -\int_{\partial D}fdx. Then we get

Df(z)dz=Dfdx+iDfdy=((yix)[χD],f)=2i(zˉ[χD],f)=2iDzˉfdxdy\int_{\partial D}f(z)dz = \int_{\partial D}fdx+i\int_{\partial D}fdy = ((\partial_y-i\partial_x)[\chi_D],f) = -2i(\partial_{\bar{z}}[\chi_D],f) = 2i\int_D \partial_{\bar{z}}f dxdy

or zˉ[χD]=12idzD\partial_{\bar z}[\chi_D] = \frac{1}{2i}dz|_{\partial D}. One can easily translate the proof using green formula or Stokes theorem. Let χD\chi_D be the characteristic function of DD, χBϵ\chi_{B_\epsilon} be the characteristic function of a small ball of radius ε centered at ζ.

Apply the identity (2) to DBϵD\setminus B_{\epsilon} and acts to the function f(z)zζ\frac{f(z)}{z-\zeta} we get Df(z)zζdzBϵf(z)zζdz=2i(DzˉfzζdxdyBϵzˉfzζdxdy)\int_{\partial D}\frac{f(z)}{z-\zeta}dz - \int_{\partial B_\epsilon}\frac{f(z)}{z-\zeta}dz = 2i(\int_D \frac{\partial_{\bar{z}}f}{z-\zeta}dxdy - \int_{B_\epsilon}\frac{\partial_{\bar{z}}f}{z-\zeta}dxdy). Note that by parametrizing Bϵ\partial B_{\epsilon} with γ(t)=ϵeit+ζ\gamma(t) = \epsilon e^{it}+\zeta, we get Bϵf(z)zζdz=i02πf(ζ+ϵeiθ)dθ\int_{\partial B_\epsilon}\frac{f(z)}{z-\zeta}dz = i\int_{0}^{2\pi}f(\zeta+\epsilon e^{i\theta})d\theta which goes to 2πf(ζ)2\pi f(\zeta) when ϵ0\epsilon \to 0, by dominated convergence theorem.

Since Bϵzˉfzζdxdy\int_{B_\epsilon}|\frac{\partial_{\bar z}f}{z-\zeta}|dxdy supBϵzfBϵ1rrdrdθ=2πϵsupBϵzˉf0 \leq \sup_{B_{\epsilon}}|\partial_z f|\cdot \int_{B_{\epsilon}} \frac{1}{r}rdrd\theta = 2\pi\epsilon \sup_{B_\epsilon}|\partial_{\bar z}f|\to 0 when ϵ0\epsilon \to 0, we get 2πif(ζ)=Df(z)zζdz2iDzˉfzζdxdy2\pi i f(\zeta) = \int_{\partial D}\frac{f(z)}{z-\zeta}dz -2i\int_{D}\frac{\partial_{\bar{z}}f}{z-\zeta}dxdy. From the notation dzdzˉ=(dx+idy)(dxidy)=2idxdydz\wedge d\bar{z} = (dx+idy)\wedge (dx-idy)=-2idx\wedge dy, we get f(ζ)=12πi(Df(z)zζdz+Dzˉf(z)zζdzdzˉ)f(\zeta) = \frac{1}{2\pi i }(\int_{\partial D}\frac{f(z)}{z-\zeta}dz + \int_D \frac{\partial_{\bar{z}}f(z)}{z-\zeta}dz\wedge d\bar{z}).

Consequences of the Cauchy integral formula

Holomorphic case

If zˉf=0\partial_{\bar{z}}f = 0 on DD, then f(ζ)=12πiDf(z)zζdzf(\zeta) = \frac{1}{2\pi i }\int_{\partial D}\frac{f(z)}{z-\zeta}dz.

This is the usual Cauchy integral formula. Note that this holds for general simple closed curves than C1\mathscr{C}^1 case, see Ahlfors’ book.

Cauchy kernel

1π1zζ\frac{1}{\pi}\frac{1}{z-\zeta} is the fundamental solution of zˉ\partial_{\bar{z}}.

Proof. Let φ be a smooth function of compact support on C\mathbb{C}. Let DD be a domain with C1\mathscr{C}^1 boundary which contains the support of φ. Then by Cauchy integral formula φ(ζ)=12πiDφ(z)zζdz1πDzˉφ(z)zζdλ(z)\varphi(\zeta) = \frac{1}{2\pi i }\int_{\partial D}\frac{\varphi(z)}{z-\zeta}dz - \frac{1}{\pi}\int_{D}\frac{\partial_{\bar z}\varphi(z)}{z-\zeta}d\lambda(z) =01πDzˉφzζ  dλ(z)= 0 -\frac{1}{\pi}\int_D \frac{\partial_{\bar z}\varphi}{z-\zeta}\;d\lambda(z) because the support of φ is inside DD and hence φ0\varphi\equiv 0 on D\partial D.

It follows that (zˉ[1π1zζ],φ)=1πDzˉφzζ  dλ(z)=φ(ζ)=(δζ,φ)(\partial_{\bar z} [\frac{1}{\pi}\frac{1}{z-\zeta}],\varphi) = -\frac{1}{\pi}\int_D \frac{\partial_{\bar z}\varphi}{z-\zeta}\;d\lambda(z) = \varphi(\zeta) = (\delta_\zeta,\varphi).