Central Potential Scattering and Phase Shifts: Difference between revisions

From PhyWiki
Jump to navigation Jump to search
No edit summary
Line 27: Line 27:
When <math>r\to\infty,\!</math> we use the asymptotic approximations of the spherical Bessel functions and the spherical Neumann functions, obtaining
When <math>r\to\infty,\!</math> we use the asymptotic approximations of the spherical Bessel functions and the spherical Neumann functions, obtaining


:<math>\frac{u(r)}{r}\approx A_l \frac{\sin(kr-l\frac{\pi}{2})}{kr} -B_l \frac{\cos(kr-l\frac{\pi}{2})}{kr}.</math>
<math>\frac{u(r)}{r}\approx A_l \frac{\sin(kr-l\frac{\pi}{2})}{kr} -B_l \frac{\cos(kr-l\frac{\pi}{2})}{kr}.</math>
   
   
Let us now define
Let us now define
Line 35: Line 35:
The angle <math>\delta_l\!</math> is known as the phase shift of the <math>\ l^{\text{th}}</math> wave and it is the phase shift induced by scattering from the potential in the radial part of the wave function. Note that, in the absence of a scattering potential, the boundary condition that the wave function must be finite at the origin causes <math>B_l</math> to vanish for all values of <math>l.</math>  Therefore, the magnitude of <math>B_l</math> compared to <math>A_l</math> is a meausre of the intensity of the scattering.  We may rewrite the above expression as
The angle <math>\delta_l\!</math> is known as the phase shift of the <math>\ l^{\text{th}}</math> wave and it is the phase shift induced by scattering from the potential in the radial part of the wave function. Note that, in the absence of a scattering potential, the boundary condition that the wave function must be finite at the origin causes <math>B_l</math> to vanish for all values of <math>l.</math>  Therefore, the magnitude of <math>B_l</math> compared to <math>A_l</math> is a meausre of the intensity of the scattering.  We may rewrite the above expression as


<math>\frac{u_l(r)}{r}=A_l\frac{\sin(kr-l\frac{\pi}{2} + \delta_l )}{kr} .</math>
<math>\frac{u_l(r)}{r}=A_l\frac{\sin(kr-l\frac{\pi}{2} + \delta_l )}{kr}.</math>


Physically, we expect <math>\delta_l < 0\!</math> for repulsive potentials and <math>\delta_l > 0\!</math> for attractive potentials.  Also, if <math>l/k \gg d,\!</math> then the classical impact parameter is much larger than the range of the potential and in this case we expect <math>\delta_l\!</math> to be small.
Physically, we expect <math>\delta_l < 0\!</math> for repulsive potentials and <math>\delta_l > 0\!</math> for attractive potentials.  Also, if <math>l/k \gg d,\!</math> then the classical impact parameter is much larger than the range of the potential and in this case we expect <math>\delta_l\!</math> to be small.
Line 158: Line 158:
In the region <math>r>a,\!</math> <math>u_l\!</math> is
In the region <math>r>a,\!</math> <math>u_l\!</math> is


:<math>u_l(r)= B j_l (kr)+C n_l(kr).\!</math>
<math>u_l(r)= B j_l (kr)+C n_l(kr).\!</math>
   
   
Here, <math>A,\!</math> <math>B,\!</math> and <math>C\!</math> are arbitrary constants and <math>k^2 =  \frac{2m}{\hbar^2}E.\!</math>
Here, <math>A,\!</math> <math>B,\!</math> and <math>C\!</math> are arbitrary constants and <math>k^2 =  \frac{2m}{\hbar^2}E.\!</math>

Revision as of 23:47, 2 September 2013

Quantum Mechanics A
SchrodEq.png
Schrödinger Equation
The most fundamental equation of quantum mechanics; given a Hamiltonian Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathcal{H}} , it describes how a state Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle |\Psi\rangle} evolves in time.
Basic Concepts and Theory of Motion
UV Catastrophe (Black-Body Radiation)
Photoelectric Effect
Stability of Matter
Double Slit Experiment
Stern-Gerlach Experiment
The Principle of Complementarity
The Correspondence Principle
The Philosophy of Quantum Theory
Brief Derivation of Schrödinger Equation
Relation Between the Wave Function and Probability Density
Stationary States
Heisenberg Uncertainty Principle
Some Consequences of the Uncertainty Principle
Linear Vector Spaces and Operators
Commutation Relations and Simultaneous Eigenvalues
The Schrödinger Equation in Dirac Notation
Transformations of Operators and Symmetry
Time Evolution of Expectation Values and Ehrenfest's Theorem
One-Dimensional Bound States
Oscillation Theorem
The Dirac Delta Function Potential
Scattering States, Transmission and Reflection
Motion in a Periodic Potential
Summary of One-Dimensional Systems
Harmonic Oscillator Spectrum and Eigenstates
Analytical Method for Solving the Simple Harmonic Oscillator
Coherent States
Charged Particles in an Electromagnetic Field
WKB Approximation
The Heisenberg Picture: Equations of Motion for Operators
The Interaction Picture
The Virial Theorem
Commutation Relations
Angular Momentum as a Generator of Rotations in 3D
Spherical Coordinates
Eigenvalue Quantization
Orbital Angular Momentum Eigenfunctions
General Formalism
Free Particle in Spherical Coordinates
Spherical Well
Isotropic Harmonic Oscillator
Hydrogen Atom
WKB in Spherical Coordinates
Feynman Path Integrals
The Free-Particle Propagator
Propagator for the Harmonic Oscillator
Differential Cross Section and the Green's Function Formulation of Scattering
Central Potential Scattering and Phase Shifts
Coulomb Potential Scattering

We will now discuss scattering from a central potential in a different way. Recall that the wave function for an incident and scattered wave for a central potential is given by

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \psi_k(\mathbf{r})=\psi_k^{(0)}(\mathbf{r})=\psi_k^{(0)}(\mathbf{r})+f_k(\theta)\frac{e^{ikr}}{r},}

where Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \psi_k^{(0)}(\mathbf{r})\!} is the incoming wave and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle f_k(\theta)\!} is the scattering amplitude.

To determine Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle f_k(\theta),\!} we start with the Schrödinger equation,

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \left( -\frac{\hbar^2 }{2m}\nabla^2+V(r) \right) \psi=\frac{\hbar^2 k^2 }{2m}\psi .}

As before, this equation may be reduced to an effective one-dimensional equation,

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \left( -\frac{\hbar^2 }{2m}\frac{d^2 }{dr^2 }+\frac{\hbar^2 l(l+1)}{2mr^2}+V(|\mathbf{r}|) \right) u_l(r) =\frac{\hbar^2 k^2 }{2m}u_l(r),}

with the full wave function given by Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \psi(\mathbf{r})=\frac{u_l(r)}{r}Y_l^m(\theta,\phi).}

For a potential with a finite range Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle d,\!} we know that, for Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle r \gg d,\!} the problem reduces to that of a free particle, and thus

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \left( -\frac{\hbar^2 }{2m}\frac{d^2 }{dr^2 }+\frac{\hbar^2 l(l+1)}{2mr^2} \right) u_l(r) =\frac{\hbar^2 k^2 }{2m}u_l(r).}

The solution of this equation is a linear combination of the spherical Bessel functions and the spherical Neumann functions,

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle u_l(r)=A_l j_l(kr) +B_l n_l(kr).\!}

When Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle r\to\infty,\!} we use the asymptotic approximations of the spherical Bessel functions and the spherical Neumann functions, obtaining

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{u(r)}{r}\approx A_l \frac{\sin(kr-l\frac{\pi}{2})}{kr} -B_l \frac{\cos(kr-l\frac{\pi}{2})}{kr}.}

Let us now define

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{B_l }{A_l }=-\tan\delta_l.}

The angle Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \delta_l\!} is known as the phase shift of the Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \ l^{\text{th}}} wave and it is the phase shift induced by scattering from the potential in the radial part of the wave function. Note that, in the absence of a scattering potential, the boundary condition that the wave function must be finite at the origin causes Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle B_l} to vanish for all values of Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle l.} Therefore, the magnitude of Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle B_l} compared to Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle A_l} is a meausre of the intensity of the scattering. We may rewrite the above expression as

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{u_l(r)}{r}=A_l\frac{\sin(kr-l\frac{\pi}{2} + \delta_l )}{kr}.}

Physically, we expect Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \delta_l < 0\!} for repulsive potentials and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \delta_l > 0\!} for attractive potentials. Also, if Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle l/k \gg d,\!} then the classical impact parameter is much larger than the range of the potential and in this case we expect Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \delta_l\!} to be small.

Because the scattering amplitude has azimuthal symmetry (i.e., it is independent of Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \phi\!} ), we can write the full solution of the Schrödinger equation as a superposition of Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle m=0\!} spherical harmonics only:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \psi(\mathbf{r})=\sum_{l=0}^{\mathop{ \infty}}a'_l(k)P_l(\cos\theta) \frac{u_l(r)}{r}=\sum_{l=0}^{\mathop{ \infty}}a_l(k)P_l(\cos\theta) \frac{\sin(kr-l\frac{\pi}{2}+\delta_l )}{kr}.}

We now determine the coefficients Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle a_l(k)\!} by substituting in the wave function in terms of the scattering amplitude on the left-hand side:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle e^{ikr\cos\theta}+f_k(\theta)\frac{e^{ikr}}{r}=\sum_{l=0}^{\mathop{ \infty}}a_l(k)P_l(\cos\theta) \frac{\sin(kr-l\frac{\pi}{2}+\delta_l )}{kr}}

Here, we assume that the incident wave propagates along the Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle z\!} direction. This must hold for large Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle r.\!} We may show that

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle e^{ikr\cos\theta}=\sum_{l=0}^{\mathop{ \infty}}(2l+1)i^l P_l(\cos\theta) \frac{\sin(kr-l\frac{\pi}{2} ) }{kr},}

so that

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \sum_{l=0}^{\mathop{ \infty}}(2l+1)i^l P_l(\cos\theta) \frac{\sin(kr-l\frac{\pi}{2} ) }{kr}+f_k(\theta)\frac{e^{ikr}}{r}=\sum_{l=0}^{\mathop{ \infty}}a_l(k)P_l(\cos\theta) \frac{\sin(kr-l\frac{\pi}{2}+\delta_l )}{kr} .}

Using the fact that

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \sin (x)=\frac{e^{ix}-e^{-ix}}{2i},\!}

we may rewrite this as

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{align} \sum_{l=0}^{\mathop{ \infty}}(2l+1)i^l P_l(\cos\theta) \frac{1}{kr} \frac{1}{2i} \left( e^{i\left( kr-l\frac{\pi}{2} \right)} - e^{-i\left( kr-l\frac{\pi}{2} \right)} \right) +f_k(\theta)\frac{e^{ikr}}{r} \\ = \sum_{l=0}^{\mathop{ \infty}}a_l(k)P_l(\cos\theta) \frac{1}{kr} \frac{1}{2i} \left( e^{i\left(kr-l\frac{\pi}{2}+\delta_l \right)} - e^{-i\left(kr-l\frac{\pi}{2}+\delta_l \right)}\right). \end{align} }

By matching the coefficients of Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle e^{-ikr}\!} , we get

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle a_l(k)=(2l+1)i^le^{i\delta_l},}

and doing the same for Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle e^{ikr}\!} yields

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle f_k(\theta)=\frac{1}{k}\sum_{l=0}^{\infty}(2l+1)e^{i\delta_l }\sin\delta_lP_l(\cos\theta) .}

Note that Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \delta_l\!} is a function of Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle k\!} and therefore a function of the incident energy. If Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \delta_l(k)\!} is known, then we can reconstruct the entire scattering amplitude and consequently the differential cross section. The phase shifts themselves must be determined by solving the Schrödinger equation.

The differential scattering cross section is

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{d\sigma}{d\Omega}=|f_k(\theta) |^2=\frac{1}{k^2 }\left|\sum_{l=0}^{\infty}(2l+1)e^{i\delta_l }\sin\delta_lP_l(\cos\theta) \right|^2.}

By integrating Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{d\sigma}{d\Omega}\!} over the solid angle Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Omega, \!} we obtain the total scattering cross section:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{align} \sigma_{tot} &= \int \frac{d\sigma}{d\Omega}\,d\Omega \\ &= \sum_{l=0}^{\infty}\sum_{l'=0}^{\infty}(2l+1)(2l'+1)e^{i\delta_l }e^{-i\delta_{l'}}\sin\delta_l\sin\delta_{l'} \int_{0}^{2\pi} d\phi\,\int_{0}^{\pi}d\theta\,\sin\theta P_l(\cos\theta)P_{l'}(\cos\theta) \\ &= \frac{4\pi}{k^2 }\sum_{l=0}^{\mathop{ \infty}}(2l+1)\sin^2\delta_l \end{align} }

The last equality follows from the orthogonality of the Legendre polynomials,

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \int_{-1}^{1}dxP_l(x) P_{l'}(x)=\frac{2}{(2l+1)}\delta_{ll'}.}

Finally, note that since for all we obtain

If we take the imaginary part of the scattering amplitude,then

Therefore,

This relationship is known as the optical theorem. The optical theorem is a general law of wave scattering theory that relates the forward scattering amplitude to the total cross section of the scattering. It was originally discovered independently by Sellmeier and Lord Rayleigh in 1871.

Referring back to the formula for the scattering amplitude, one more important quantity can be discussed:

This quantity, for now referred to as the partial scattering for angular momentum is the ratio of the coefficients of the outgoing and incoming waves for a wave scattered on a potential of finite range

These ratios can simplify the problem of evaluating the continuity of the waveform at the boundary In general, if the interior wave function is known to be smoothly continuous across the boundary at then the phase shifts can be expressed in terms of the logarithmic derivatives evaluated at the boundary

Using the above equations for the form of beyond the region of scattering, the following relation is found:

Thus, after some algebra,

Note that, if which corresponds to then only the first portion of this expression survives. This is a special quantity corresponding to hard sphere scattering; we may define the phase angles known as the hard sphere phase shifts, as

Note that these phase shifts are present for any potential, not just that of a hard sphere.

Scattering by Square Well potential

Consider a beam of point particles of mass scattering from a finite spherical attractive well of depth and radius

The effective Schrödinger equation for is

Its solution is

where

In the region is

Here, and are arbitrary constants and

For large

where

We now apply the boundary conditions at which are continuity of and of its logarithmic derivative. We obtain

Let us now consider two limiting cases:

(a) and In this case, we find that, with some simplification, This behavior is a result of the centrifugal barrier that keeps waves of energy far below the barrier from feeling the effect of the potential.

(b) When the phase shift is the partial wave cross section is maximized. This is known as a resonant scattering.

From (a), we see that the phase shift is small for small. However, when changes and passes through the resonance condition, the phase shift rises rapidly and has a sharp peak at resonant energy . This can be represented as

so that the partial wave cross section is

This is the Breit-Wigner formula for a resonant cross section.

Problems

(1) Consider the scattering of a particle from a real spherically symmetric potential. If is the differential cross section and is the total cross section, then show that

for a general central potential using the partial-wave expansion of the scattering amplitude and the cross section.

Solution

(2) Consider an attractive delta shell potential () of the form,

(a) Derive the equation for the phase shift caused by this potential for arbitrary angular momentum.

(b) Obtain the expression for the wave phase shift.

(c) Obtain the wave scattering amplitude.

Solution