Integral Operators

Note that any operator K on a function space U may be expressed as an integral,

K=IKI=abdxabdxw(x)w(x)|xx|K|xx|=abdxdxw(x)w(x)x|K|x|xx|=abdxdxw(x)w(x)K(x,x)|xx|

so that Kxf(x)x|K|f=x|dxdxw(x)w(x)K(x,x)|xx|f=dxdxw(x)w(x)K(x,x)δ(xx)w(x)f(x)=dxw(x)K(x,x)f(x)

If K(x,x) is just a function (not a distribution or generalized function) we K an integral operator.

We can think of K(x,x) as a matrix element wich x the row and x the column and f(x) the column vector. Then (as homework) K(x,x)=K(x,x), just like a matrix.

K(x,x) is called the integral kernel of K.

(Example the momentum operator as seen Monday)

Inverses of Differential Operators

We may expect it to involve integration (to undo the differentiation). L1 an integral operator. We get constants of integration, which can be fixed with boundary conditions. We may have to eliminate the nullspace, Differential opeators have non-empty nullspaces but sometimes the homogeneous solutions cannot satisfy the boundary conditions which would eliminate the nullspace from the domain of the operator.

So let’s procede, assume L=L is self adjoint.

Solve L|u=|s on N(L) i.e. the subspace spanned by eigenvectors with nonzero eigenvalues. Let’s call L1G. I.e. L1L=LL1=LG=I. When G exists, it’s called a Green’s operator. In general, there may be more than one Green’s operator depending on the boundary conditions.

Properties of Green’s Operator

x|x=δ(xx)w(x)=x|LG|x=x|L|Gx=Lxx|G|x=LxG(x,x)

So, LxG(x,x)=δ(xx)w(x) is a differential equation for the Green’s operator. I.e. the Green’s operator satisfies the differential equation. It is the same differential equation we had before L|u=|s but with a delta function source (e.g. point charge). The specific solution depends on the bounary conditions. Lets say we found G(x,x), then |u=G|s. Or, x|u=x|G|s=u(x)=abdxw(x)G(x,x)s(x). This gives the particular solution.

Finding Green Functions

Solve the differential Equation

LxG(x,x)=δ(xx)w(x) that satisfies the problems’ given boundary conditions. So, solves the homogeneous equation for all x not x, at x equal to x it is discontinuous in the last order derivative in L of the Green’s function. This is called a junction condition at x=x.

D&K does this for 2nd order linear operators and does example for a string.

For Self Adjoint Operators, Expand the Green Function in the Eigenfunctions of L

Fourier Transform

  • Example: Simple Harmonic Oscillator

    mx¨+mω2x=F(t) [d2dt2+ω2]x(t)=Tx(t)=1mF(t)=s(t). So, LtG(t,t)=δ(tt).

    The solution is x(t)=dt(G(t,t)s(t))+x0(t) and fixed with boundary conditions.

    Note coefficients in Lt=d2dt2=ω2 are independent of time. Thus, we can transform tt+α to no effect. Hence, G(t,t)=G(tt).

    Define ttt. So, LtG(t)=δ(t). Then, Lt12πdEexp(iEt)G~(E)=δ(t). So, Ltexp(iEt)=(E2+ω2)exp(iEt). Then, LtG(t)=12πdE(ω2E2)exp(iEt)G~(t)=δ(t).

    Then, dtexp(iEt)12πdE(ω2E2)exp(iEt)G~(t)=dtexp(iEt)δ(t).

    Then, 12πdE(ω2E2)dtexp(i(EE)t)G~(t)=exp(iEt).

    From a definition of the Dirac Delta Function, δ(t)=12πexp±iωt.

    Then, 2π(ω2E2)G~(E)=1.

    Hence, G~(E)=12π1ω2E2.

    G(t)=12πRdEexp(iEt)G~(E)=12πexp(iEt)E2ω2. We have an integral we will need to solve with Complex analysis.

    RdEexp(iEt)(Eω)(E+ω)=RdEg(E), simple poles at ±ω.

    This is a singular integral, so it is ambiguous. Thus, we must choose the contour.

    • Principal Value (Non-singular bit, limit as you go straight through the singularities) [This would be a solution without the delta function]
    • Small semicircles around the poles (4 possibilities) [This gives the solutions with the delta function]

    5 choices! Our boundary conditions are needed to distinguish which solution to use.

    Res g(ω)=limzωexp(izt)2z=exp(iωt)2ω, hence for our integral Res(ω)=exp(iωt)4πω

    Res g(ω)=limzωexp(izt)2z=exp(iωt)2ω, hence for our integral Res(ω)=exp(iωt)4πω

    Since exp(iEt)=exp(itRcosθ)exp(tRsinθ) we must choose the upper contour for positive t and the lower contour for negative t.

    Choose the contour that goes into the negative imaginary axis. If we use Γ+ then the contour encloses the poles. If we use Γ then the contour does not enclose the poles.

    t>0: G(t)=2πiResenclosed=2πi1ωsin(ωt) t<0: G(t)=0 Hence, G(t)=sinωtωΘ(t).

    This choice of contour gives rise to the retarded Green function. The opposite choice gives rise to the advanced Green function. The cross is the Feynman Green function.

Author: Christian Cunningham

Created: 2024-05-30 Thu 21:15

Validate