221B Lecture Notes Notes on Spherical Bessel Functions

Similar documents
221B Lecture Notes Notes on Spherical Bessel Functions

Legendre Polynomials - Lecture 8

David Eigen. MA112 Final Paper. May 10, 2002

Lecture 8 February 18, 2010

Physics 116C Helmholtz s and Laplace s Equations in Spherical Polar Coordinates: Spherical Harmonics and Spherical Bessel Functions

Week 6 Lectures, Math 6451, Tanveer

Math 124B January 17, 2012

Math 124B January 31, 2012

An approximate method for solving the inverse scattering problem with fixed-energy data

Radiation Fields. Lecture 12

Lecture Notes for Math 251: ODE and PDE. Lecture 32: 10.2 Fourier Series

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES

Physics 505 Fall 2007 Homework Assignment #5 Solutions. Textbook problems: Ch. 3: 3.13, 3.17, 3.26, 3.27

FFTs in Graphics and Vision. Spherical Convolution and Axial Symmetry Detection

SEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l

4 Separation of Variables

Separation of Variables and a Spherical Shell with Surface Charge

Physics 506 Winter 2006 Homework Assignment #6 Solutions

LECTURE NOTES 8 THE TRACELESS SYMMETRIC TENSOR EXPANSION AND STANDARD SPHERICAL HARMONICS

More Scattering: the Partial Wave Expansion

Physics 235 Chapter 8. Chapter 8 Central-Force Motion

$, (2.1) n="# #. (2.2)

Lecture Notes 4: Fourier Series and PDE s

2M2. Fourier Series Prof Bill Lionheart

Quantum Mechanical Models of Vibration and Rotation of Molecules Chapter 18

11.1 One-dimensional Helmholtz Equation

Higher dimensional PDEs and multidimensional eigenvalue problems

Fourier Series. 10 (D3.9) Find the Cesàro sum of the series. 11 (D3.9) Let a and b be real numbers. Under what conditions does a series of the form

6 Wave Equation on an Interval: Separation of Variables

APPENDIX B. Some special functions in low-frequency seismology. 2l +1 x j l (B.2) j l = j l 1 l +1 x j l (B.3) j l = l x j l j l+1

CS229 Lecture notes. Andrew Ng

Assignment 7 Due Tuessday, March 29, 2016

4 1-D Boundary Value Problems Heat Equation

Strauss PDEs 2e: Section Exercise 2 Page 1 of 12. For problem (1), complete the calculation of the series in case j(t) = 0 and h(t) = e t.

MAT 167: Advanced Linear Algebra

MA 201: Partial Differential Equations Lecture - 10

Volume 13, MAIN ARTICLES

Physics 505 Fall Homework Assignment #4 Solutions

Jackson 4.10 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell

Formulas for Angular-Momentum Barrier Factors Version II

Applied Nuclear Physics (Fall 2006) Lecture 7 (10/2/06) Overview of Cross Section Calculation

Notes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.

Data Search Algorithms based on Quantum Walk

PHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased

TRANSFORMATION OF REAL SPHERICAL HARMONICS UNDER ROTATIONS

Gaussian Curvature in a p-orbital, Hydrogen-like Atoms

14 Separation of Variables Method

Vector Spherical Harmonics

Lecture 9. Stability of Elastic Structures. Lecture 10. Advanced Topic in Column Buckling

LECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL HARMONICS

1D Heat Propagation Problems

Math 220B - Summer 2003 Homework 1 Solutions

17 Lecture 17: Recombination and Dark Matter Production

1.2 Partial Wave Analysis

MAS 315 Waves 1 of 8 Answers to Examples Sheet 1. To solve the three problems, we use the methods of 1.3 (with necessary changes in notation).

This shoud be true for any, so the equations are equivaent to the Schrodingerequation (H E) : () Another important property of the (2) functiona is, t

Lecture VIII : The pseudopotential

Section 6: Magnetostatics

Using the Green s Function to find the Solution to the Wave. Equation:

Introduction to LMTO method

arxiv: v1 [math.ca] 6 Mar 2017

ON THE REPRESENTATION OF OPERATORS IN BASES OF COMPACTLY SUPPORTED WAVELETS

Course 2BA1, Section 11: Periodic Functions and Fourier Series

VTU-NPTEL-NMEICT Project

Midterm 2 Review. Drew Rollins

Lecture 11. Fourier transform

Massive complex scalar field in the Kerr Sen geometry: Exact solution of wave equation and Hawking radiation

Chapter 5. Wave equation. 5.1 Physical derivation

Motion in Spherically Symmetric Potentials

HYDROGEN ATOM SELECTION RULES TRANSITION RATES

Vibrations of Structures

C. Fourier Sine Series Overview

Wave Equation Dirichlet Boundary Conditions

Lecture Note 3: Stationary Iterative Methods

Physics 127c: Statistical Mechanics. Fermi Liquid Theory: Collective Modes. Boltzmann Equation. The quasiparticle energy including interactions

Term Test AER301F. Dynamics. 5 November The value of each question is indicated in the table opposite.

THINKING IN PYRAMIDS

Scattering of Particles by Potentials

arxiv: v1 [hep-lat] 23 Nov 2017

ODE Homework 2. Since M y N x, the equation is not exact. 2. Determine whether the following equation is exact. If it is exact, M y N x 1 x.

Lecture 17 - The Secrets we have Swept Under the Rug

Homework 5 Solutions

On the Number of Limit Cycles for Discontinuous Generalized Liénard Polynomial Differential Systems

Module 22: Simple Harmonic Oscillation and Torque

Electromagnetism in the Spherical-Wave Basis:

SVM: Terminology 1(6) SVM: Terminology 2(6)

Strauss PDEs 2e: Section Exercise 1 Page 1 of 7

BALANCING REGULAR MATRIX PENCILS

On Integrals Involving Universal Associated Legendre Polynomials and Powers of the Factor (1 x 2 ) and Their Byproducts

Identites and properties for associated Legendre functions

Introduction to Simulation - Lecture 13. Convergence of Multistep Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy

1 Heat Equation Dirichlet Boundary Conditions

arxiv: v1 [math.ca] 31 Aug 2017

Atom-Atom Interactions in Ultracold Quantum Gases

Algorithms to solve massively under-defined systems of multivariate quadratic equations

6.434J/16.391J Statistics for Engineers and Scientists May 4 MIT, Spring 2006 Handout #17. Solution 7

EXISTENCE OF SOLUTIONS FOR ONE-DIMENSIONAL WAVE EQUATIONS WITH NONLOCAL CONDITIONS

THE THREE POINT STEINER PROBLEM ON THE FLAT TORUS: THE MINIMAL LUNE CASE

Jost Function for Singular Potentials

hole h vs. e configurations: l l for N > 2 l + 1 J = H as example of localization, delocalization, tunneling ikx k

Transcription:

Definitions B Lecture Notes Notes on Spherica Besse Functions We woud ike to sove the free Schrödinger equation [ h d r R(r) = h k R(r). () m r dr r m R(r) is the radia wave function ψ( x) = R(r)Y m (θ, φ). By factoring out h /m and defining = kr, we find the equation [ d + R() = 0. () d The soutions to this equation are spherica Besse functions. Due to some reason, I don t see the integra representations I use beow in books on mathemtica formuae, but I beieve they are right. The behavior at the origin can be studied by power expansion. Assuming R n, and coecting terms of the owest power in, we get There are two soutions, n(n + ) = 0. (3) n = or. (4) The first soution gives a positive power, and hence a reguar soution at the origin, whie the second a negative power, and hence a singuar soution at the origin. It is easy to check that the foowing integra representations sove the above equation Eq. (): and () = (/) i e it ( t ) dt, (5)! + h () () = (/) i e it ( t ) dt. (6)!

By acting the derivatives in Eq. (), one finds [ d + d () = (/)! = (/)! i ± i i [ ( t ) + ± d dt ( + )it t + dt [ e it ( t ) + dt. (7) Therefore ony boundary vaues contribute, which vanish both at t = and t = i for = kr > 0. The same hods for h () (). One can aso easiy see that () = h () ( ) by taking the compex conjugate of the expression Eq. (5) and changing the variabe from t to t. The integra representation Eq. (5) can be expanded in powers of /. For instance, for, we change the variabe from t to x by t = + ix, and find () = (/) ( e i(+ix) x ( i) x idx! 0 i) Simiary, we find = i (/) e i ( i)! = i ei h () () = i e i ( i) k ( + k)! k k!( k)! C k 0 e x ( x i) k x dx. (8) k i k ( + k)! k k!( k)!. (9) k Therefore both h (,) are singuar at = 0 with power. The combination j () = ( + h () )/ is reguar at = 0. This can be seen easiy as foows. Because h () is an integra from t = to i, whie from t = + to i, the differencd between the two corresponds to an integra from t = to t = i and coming back to t = +. Because the integrand does not have a poe, this contour can be deformed to a straight integra from t = to +. Therefore, j () = (/) e it ( t ) dt. (0)!

In this expression, 0 can be taken without any probems in the integra and hence j, i.e., reguar. The other inear combination n = ( )/i is of course singuar at = 0. Note that h () is anaogous to It is usefu to see some exampes for ow. () = j () + i n () () e i = cos + i sin. () j 0 = sin, j = sin cos, j = 3 n 0 = cos, n = cos sin 0 = i ei, h() = i ( h () 0 = i e i, h() = i ( + i i sin 3 cos, 3, n = 3 cos 3 sin, ) 3 e i = i ( ) 3 3i 3 e i. ) e i h () = i ( ) 3 + 3i 3 e i. (3) Asymptotic Behavior Eqs. (8,9) give the asymptotic behaviors of for : By taking inear combinations, we aso find 3 Pane Wave Expansion i ei ( i) = i ei( π/). (4) j sin( π/), (5) n cos( π/). (6) The non-trivia ooking formua we used in the cass e ikz = ( + )i j (kr)p (cos θ) (7) =0 can be obtained quite easiy from the integra representation Eq. (0). The point is that one can keep integrating it in parts. By integrating e it factor 3

and differentiating ( t ) factor, the boundary terms at t = ± aways vanish up to -th time because of the ( t ) factor. Therefore, j = (/) (! (i) eit dt) d ( t ) dt. (8) Note that the definition of the Legendre poynomias is P (t) = d! dt (t ). (9) Using this definition, the spherica Besse function can be written as j = i e it P (t)dt. (0) Then we use the fact that the Legendre poynomias form a compete set of orthogona poynomias in the interva t [,. Noting the normaization P n (t)p m (t)dt = n + δ n,m, () the orthonorma basis is P n (t) (n + )/, and hence n=0 n + P n (t)p n (t ) = δ(t t ). () By mutipying Eq. (0) by P (t )( + )/ and summing over n, n= + P (t )j n () = i n e it P (t )P (t)dt = n=0 i n eit. (3) By setting = kr and t = cos θ, we prove Eq. (7). If the wave vector is pointing at other directions than the positive z- axis, the formua Eq. (7) needs to be generaized. Noting Y 0 (θ, φ) = ( + )/4π P (cos θ), we find e i k x = 4π i j (kr) =0 m= Y m (θ k, φ k )Y m (θ x, φ x ) (4) 4

4 Deta-Function Normaization An important consequence of the identity Eq. (4) is the innerproduct of two spherica Besse functions. We start with d xe i k x e i k x = (π) 3 δ( k k ). (5) Using Eq. (4) in the.h.s of this equation, we find d xe i k x e i k x = (4π),m,m = (4π),m dω x drr Y m (Ω k )Y m (Ω x )Y m (Ω x)y m (Ω k )j (kr)j (k r) drr j (kr)j (k r)y m (Ω k )Y m (Ω k ). (6) On the other hand, the r.h.s. of Eq. (5) is (π) 3 δ( k k ) = (π) 3 k δ(k k )δ(ω k Ω k ) Comparing Eq. (6) and (7) and noting,m = (π) 3 k sin θ δ(k k )δ(θ θ )δ(φ φ ). (7) Y m (Ω k )Y m (Ω k ) = δ(ω k Ω k ), (8) we find 0 drr j (kr)j (k r) = π k δ(k k ). (9) 5