Math Theory of Partial Differential Equations Lecture 3-2: Spectral properties of the Laplacian. Bessel functions.

Similar documents
MATH 311 Topics in Applied Mathematics Lecture 25: Bessel functions (continued).

The heat and wave equations in 2D and 3D

Final Examination Linear Partial Differential Equations. Matthew J. Hancock. Feb. 3, 2006

Eigenvalue comparisons in graph theory

Last Update: April 7, 201 0

FFTs in Graphics and Vision. The Laplace Operator

SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS

REMARKS ON THE MEMBRANE AND BUCKLING EIGENVALUES FOR PLANAR DOMAINS. Leonid Friedlander

MATH FALL 2014 HOMEWORK 10 SOLUTIONS

In what follows, we examine the two-dimensional wave equation, since it leads to some interesting and quite visualizable solutions.

Asymptotic distribution of eigenvalues of Laplace operator

Math 4263 Homework Set 1

24 Solving planar heat and wave equations in polar coordinates

MATH 304 Linear Algebra Lecture 20: The Gram-Schmidt process (continued). Eigenvalues and eigenvectors.

CHAPTER VIII HILBERT SPACES

Math 300 Winter 2011 Advanced Boundary Value Problems I. Bessel s Equation and Bessel Functions

Vectors in Function Spaces

AMS 212A Applied Mathematical Methods I Appendices of Lecture 06 Copyright by Hongyun Wang, UCSC. ( ) cos2

Section Consider the wave equation: ρu tt = T 0 u xx + αu + βu t

Notes on Special Functions

# Points Score Total 100

Isoperimetric inequalities and variations on Schwarz s Lemma

MATH 220: INNER PRODUCT SPACES, SYMMETRIC OPERATORS, ORTHOGONALITY

Spectral Gap for Complete Graphs: Upper and Lower Estimates

Quantum Computing Lecture 2. Review of Linear Algebra

SOLUTIONS TO SELECTED PROBLEMS FROM ASSIGNMENTS 3, 4

The spectral zeta function

SOLUTIONS TO THE FINAL EXAM. December 14, 2010, 9:00am-12:00 (3 hours)


Recitation 1 (Sep. 15, 2017)

Sturm-Liouville operators have form (given p(x) > 0, q(x)) + q(x), (notation means Lf = (pf ) + qf ) dx

HW7, Math 322, Fall 2016

MATH 205 HOMEWORK #3 OFFICIAL SOLUTION. Problem 1: Find all eigenvalues and eigenvectors of the following linear transformations. (a) F = R, V = R 3,

Definition and basic properties of heat kernels I, An introduction

1 A complete Fourier series solution

Method of Green s Functions

Lecture IX. Definition 1 A non-singular Sturm 1 -Liouville 2 problem consists of a second order linear differential equation of the form.

Essential Spectra of complete manifolds

Ordinary Differential Equations II

Gradient, Divergence and Curl in Curvilinear Coordinates

Partial Differential Equations

On the spectrum of the Laplacian

MATH325 - QUANTUM MECHANICS - SOLUTION SHEET 11

Math 3363 Examination I Solutions

1 Distributions (due January 22, 2009)

Electromagnetism HW 1 math review

Green s Theorem. MATH 311, Calculus III. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Green s Theorem

Kevin James. MTHSC 206 Section 16.4 Green s Theorem

Solutions: Problem Set 3 Math 201B, Winter 2007

Eigenvalues and Eigenfunctions of the Laplacian

Representation theory and quantum mechanics tutorial Spin and the hydrogen atom

f(s) e -i n π s/l d s

Physics 6303 Lecture 15 October 10, Reminder of general solution in 3-dimensional cylindrical coordinates. sinh. sin

MATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators.

Composite Membranes. 8th August 2006

Preparation for the Final

ENGI 4430 Line Integrals; Green s Theorem Page 8.01

ADIABATIC PHASES IN QUANTUM MECHANICS

ENGI 4430 Line Integrals; Green s Theorem Page 8.01

MATH MEASURE THEORY AND FOURIER ANALYSIS. Contents

MATH 124B: INTRODUCTION TO PDES AND FOURIER SERIES

MATH 590: Meshfree Methods

PHYS 404 Lecture 1: Legendre Functions

21 Laplace s Equation and Harmonic Functions

Waves on 2 and 3 dimensional domains

Quantum Ergodicity and Benjamini-Schramm convergence of hyperbolic surfaces

Nodal lines of Laplace eigenfunctions

Separation of Variables in Linear PDE: One-Dimensional Problems

Linear Algebra and Dirac Notation, Pt. 1

Before you begin read these instructions carefully:

The Postulates. What is a postulate? Jerry Gilfoyle The Rules of the Quantum Game 1 / 21

Sturm-Liouville Theory

Problem set 3: Solutions Math 207B, Winter Suppose that u(x) is a non-zero solution of the eigenvalue problem. (u ) 2 dx, u 2 dx.

Functional Analysis Review

Math October 20, 2006 Exam 2: Solutions

Physics 505 Homework No. 1 Solutions S1-1

MATH 31CH SPRING 2017 MIDTERM 2 SOLUTIONS

Two special equations: Bessel s and Legendre s equations. p Fourier-Bessel and Fourier-Legendre series. p

Linear Operators, Eigenvalues, and Green s Operator

Math 265H: Calculus III Practice Midterm II: Fall 2014

The purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and engineering.

ON APPROXIMATION OF LAPLACIAN EIGENPROBLEM OVER A REGULAR HEXAGON WITH ZERO BOUNDARY CONDITIONS 1) 1. Introduction

( ) = 9φ 1, ( ) = 4φ 2.

Page 404. Lecture 22: Simple Harmonic Oscillator: Energy Basis Date Given: 2008/11/19 Date Revised: 2008/11/19

Analysis IV : Assignment 3 Solutions John Toth, Winter ,...). In particular for every fixed m N the sequence (u (n)

Math 46, Applied Math (Spring 2009): Final

TOEPLITZ OPERATORS. Toeplitz studied infinite matrices with NW-SE diagonals constant. f e C :

Notes on Integrable Functions and the Riesz Representation Theorem Math 8445, Winter 06, Professor J. Segert. f(x) = f + (x) + f (x).

QUASI-SEPARATION OF THE BIHARMONIC PARTIAL DIFFERENTIAL EQUATION

Differential Geometry and Lie Groups with Applications to Medical Imaging, Computer Vision and Geometric Modeling CIS610, Spring 2008

AN INTRODUCTION TO CURVILINEAR ORTHOGONAL COORDINATES

GENERALIZED STABILITY OF THE TWO-LAYER MODEL

Lecture 12: Detailed balance and Eigenfunction methods

4. Solvability of elliptic PDEs

Outline 1. Real and complex p orbitals (and for any l > 0 orbital) 2. Dirac Notation :Symbolic vs shorthand Hilbert Space Vectors,

3.3 Unsteady State Heat Conduction

ON THE EXISTENCE OF TRANSMISSION EIGENVALUES. Andreas Kirsch1

Infinite Series. 1 Introduction. 2 General discussion on convergence

ACM/CMS 107 Linear Analysis & Applications Fall 2017 Assignment 2: PDEs and Finite Element Methods Due: 7th November 2017

Line integrals of 1-forms on the Sierpinski gasket

Transcription:

Math 412-501 Theory of Partial ifferential Equations Lecture 3-2: Spectral properties of the Laplacian. Bessel functions.

Eigenvalue problem: 2 φ + λφ = 0 in, ( αφ + β φ ) n = 0, where α, β are piecewise continuous real functions on such that α + β 0 everywhere on. We assume that the boundary is piecewise smooth.

6 spectral properties of the Laplacian Property 1. All eigenvalues are real. Property 2. All eigenvalues can be arranged in the ascending order λ 1 < λ 2 <... < λ n < λ n+1 <... so that λ n as n. This means that: there are infinitely many eigenvalues; there is a smallest eigenvalue; on any finite interval, there are only finitely many eigenvalues.

Property 3. An eigenvalue λ n may be multiple but its multiplicity is finite. Moreover, the smallest eigenvalue λ 1 is simple, and the corresponding eigenfunction φ 1 has no zeros inside the domain. Property 4. Eigenfunctions corresponding to different eigenvalues are orthogonal relative to the inner product f, g = f (x, y)g(x, y) dx dy.

Property 5. Any eigenfunction φ can be related to its eigenvalue λ through the Rayleigh quotient: φ φ n ds + φ 2 dx dy λ =. φ 2 dx dy

Property 6. There exists a sequence φ 1, φ 2,... of pairwise orthogonal eigenfunctions that is complete in the Hilbert space L 2 (). Any square-integrable function f L 2 () is expanded into a series f (x, y) = n=1 c nφ n (x, y), that converges in the mean. The series is unique: c n = f, φ n φ n, φ n. If f is piecewise smooth then the series converges pointwise to f at points of continuity.

Rayleigh quotient Suppose that 2 φ = λφ in the domain. Multiply both sides by φ and integrate over : φ 2 φ dx dy = λ φ 2 dx dy. Green s formula: ψ 2 φ da = ψ φ n ds ψ φ da This is an analog of integration by parts. Now φ φ n ds φ 2 dx dy = λ φ 2 dx dy.

It follows that λ = φ φ n ds + φ 2 dx dy. φ 2 dx dy If φ satisfies the boundary condition φ = 0 or φ n = 0 (or mixed), then the one-dimensional integral vanishes. In particular, λ 0. If φ n + αφ = 0 on, then φ φ n ds = α φ 2 ds. In particular, if α 0 everywhere on, then λ 0.

Self-adjointness ψ 2 φ dx dy = ψ φ n ds ψ φ dx dy (Green s first identity) ( (φ 2 ψ ψ 2 φ) dx dy = φ ψ n ψ φ ) ds n (Green s second identity) If φ and ψ satisfy the same boundary condition ( αφ + β φ ) ( n = αψ + β ψ ) n = 0 then φ ψ n ψ φ n = 0 everywhere on.

If φ and ψ satisfy the same boundary condition then (φ 2 ψ ψ 2 φ) dx dy = 0. If φ and ψ are complex-valued functions then also (φ 2 ψ ψ 2 φ) dx dy = 0 (because 2 ψ = 2 ψ and ψ satisfies the same boundary condition as ψ). Thus 2 φ, ψ = φ, 2 ψ, where f, g = f (x, y)g(x, y) dx dy.

Eigenvalue problem: 2 φ + λφ = 0 in, ( αφ + β φ ) n = 0. The Laplacian 2 is self-adjoint in the subspace of functions satisfying the boundary condition. Suppose φ is an eigenfunction belonging to an eigenvalue λ. Let us show that λ R. Since 2 φ = λφ, we have that 2 φ, φ = λφ, φ = λ φ, φ, φ, 2 φ = φ, λφ = λ φ, φ. Now 2 φ, φ = φ, 2 φ and φ, φ > 0 imply λ R.

Suppose φ 1 and φ 2 are eigenfunctions belonging to different eigenvalues λ 1 and λ 2. Let us show that φ 1, φ 2 = 0. Since 2 φ 1 = λ 1 φ 1, 2 φ 2 = λ 2 φ 2, we have that 2 φ 1, φ 2 = λ 1 φ 1, φ 2 = λ 1 φ 1, φ 2, φ 1, 2 φ 2 = φ 1, λ 2 φ 2 = λ 2 φ 1, φ 2. But 2 φ 1, φ 2 = φ 1, 2 φ 2, hence λ 1 φ 1, φ 2 = λ 2 φ 1, φ 2. We already know that λ 2 = λ 2. Also, λ 1 λ 2. It follows that φ 1, φ 2 = 0.

The main purpose of the Rayleigh quotient Consider a functional (function on functions) φ φ n ds + φ 2 dx dy RQ[φ] =. φ 2 dx dy If φ is an eigenfunction of 2 in the domain with some boundary condition, then RQ[φ] is the corresponding eigenvalue. What if φ is not?

Let λ 1 < λ 2 λ 3... λ n λ n+1... be eigenvalues of a particular eigenvalue problem counted with multiplicities. That is, a simple eigenvalue appears once in this sequence, an eigenvalue of multiplicity two appears twice, and so on. There is a complete orthogonal system φ 1, φ 2,... in the Hilbert space L 2 () such that φ n is an eigenfunction belonging to λ n.

Theorem (i) λ 1 = min RQ[φ], where the minimum is taken over all nonzero functions φ which are differentiable in and satisfy the boundary condition. Moreover, if RQ[φ] = λ 1 then φ is an eigenfunction. (ii) λ n = min RQ[φ], where the minimum is taken over all nonzero functions φ which are differentiable in, satisfy the boundary condition, and such that φ, φ k = 0 for 1 k < n. Moreover, the minimum is attained only on eigenfunctions. Main idea of the proof: RQ[φ] = 2 φ, φ. φ, φ (see Haberman 5.6)

Spectral properties of the Laplacian in a circle Eigenvalue problem: 2 φ + λφ = 0 in = {(x, y) : x 2 + y 2 R 2 }, u = 0. In polar coordinates (r, θ): 2 φ r 2 + 1 r φ r + 1 2 φ r 2 θ + λφ = 0 2 (0 < r < R, π < θ < π), φ(r, θ) = 0 ( π < θ < π).

Additional boundary conditions: φ(r, π) = φ(r, π), φ(0, θ) < φ θ ( π < θ < π), φ (r, π) = θ (r, π) Separation of variables: φ(r, θ) = f (r)h(θ). Substitute this into the equation: (0 < r < R). f (r)h(θ) + r 1 f (r)h(θ) + r 2 f (r)h (θ) + λf (r)h(θ) = 0. ivide by f (r)h(θ) and multiply by r 2 : r 2 f (r) + r f (r) + λr 2 f (r) f (r) + h (θ) h(θ) = 0.

It follows that r 2 f (r) + r f (r) + λr 2 f (r) f (r) The variables have been separated: = h (θ) h(θ) = µ = const. r 2 f + rf + (λr 2 µ)f = 0, h = µh. Boundary conditions φ(r, θ) = 0 and φ(0, θ) < hold if f (R) = 0 and f (0) <. Boundary conditions φ(r, π) = φ(r, π) and φ φ θ (r, π) = θ (r, π) hold if h( π) = h(π) and h ( π) = h (π).

Eigenvalue problem: h = µh, h( π) = h(π), h ( π) = h (π). Eigenvalues: µ m = m 2, m = 0, 1, 2,.... µ 0 = 0 is simple, the others are of multiplicity 2. Eigenfunctions: h 0 = 1, h m (θ) = cos mθ and h m (θ) = sin mθ for m 1.

ependence on r: r 2 f + rf + (λr 2 µ)f = 0, f (R) = 0, f (0) <. We may assume that µ = m 2, m = 0, 1, 2,.... Also, we know that λ > 0 (Rayleigh quotient!). New variable z = λ r removes dependence on λ: z 2 d2 f dz 2 + z df dz + (z2 m 2 )f = 0. This is Bessel s differential equation of order m. Solutions are called Bessel functions of order m.