Coulomb Potential Scattering: Difference between revisions
Line 43: | Line 43: | ||
==Exact Coulomb Scattering Cross Section== | ==Exact Coulomb Scattering Cross Section== | ||
[[Image:ParabolicCoordinates.png|thumb|400px|left|'''Parabolic coordinates'''<br/><math>\phi\!</math> represents rotation about the z-axis, <math>\xi\!</math> represents the parabolas with their | [[Image:ParabolicCoordinates.png|thumb|400px|left|'''Parabolic coordinates'''<br/><math>\phi\!</math> represents rotation about the z-axis, <math>\xi\!</math> represents the parabolas with their vertices at a minimum, and <math>\eta\!</math> represents parabolas with their vertices at a maximum.]] | ||
When we are considering scattering due to the Coulomb potential, we can not neglect the effect of this potential at large distances because it is only a <math>\frac{1}{r}</math> potential. | When we are considering scattering due to the Coulomb potential, we can not neglect the effect of this potential at large distances because it is only a <math>\frac{1}{r}</math> potential. | ||
We will work in parabolic coordinates, which are related to Cartesian coordinates by | |||
<math>\xi=\sqrt{x^2+y^2+z^2}-z</math> | |||
<math>\eta=\sqrt{x^2+y^2+z^2}+z</math> | |||
<math>\phi=\tan^{-1}\left (\frac{y}{x}\right )</math> | |||
The [[Schrödinger Equation|Schrödinger equation]] in parabolic coordinates is | |||
<math> \left[ - \frac{\hbar^2}{2\mu}\frac{4}{\xi+\eta} \left( \frac{\partial}{\partial \xi}\xi\frac{\partial}{\partial \xi}+\frac{\partial}{\partial \eta}\eta\frac{\partial}{\partial \eta}+\frac{\xi+\eta}{4\xi\eta}\frac{\partial^2}{\partial \phi^2} \right) -\frac{2Ze^2}{\xi+\eta} \right] \psi=\frac{\hbar^2k^2}{2\mu}\psi.</math> | |||
Recall that, for a spherically symmetric potential, the scattering amplitude is a function of <math>\theta\!</math> only; therefore, we will seek solutions that are independent of <math>\phi.\!</math> | |||
We | We will look for solution of the form, | ||
<math>\psi(\xi,\eta,\phi)=e^{\frac{i}{2}k(\eta-\xi)}\Phi(\xi)=e^{ikz}\Phi(r-z).</math> | |||
The [[Schrödinger Equation|Schrödinger equation]] then becomes | |||
<math> \left[ \xi\frac{\partial^2}{\partial \xi^2}+(1-ik\xi)\frac{\partial}{\partial \xi}+\frac{Ze^2\mu}{k\hbar^2}k \right]\Phi(\xi)=0.</math> | |||
We now assume a series solution, | |||
<math>\Phi(\xi)=\sum_{n=0}^\infty a_n \xi^n.</math> | |||
The recursion relation for the coefficients is then | |||
<math>\frac{a_{n+1}}{a_{n}}=\frac{in-\lambda}{(n+1)^2}k=\frac{n+i\lambda}{(n+1)^2}ik,</math> | |||
where | |||
<math>\lambda=\frac{Ze^2\mu}{k\hbar^2}.</math> | |||
Recall that the confluent hypergeometric function <math>{}_1F_1</math> is given by | |||
<math>_1F_1(a,c,z)=1+\frac{a}{c}z+\frac{a(a+1)}{c(c+1)}\frac{z^2}{2!}+\dots+\frac{a(a+1)\dots(a+n-1)}{c(c+1)\dots(c+n-1)}\frac{z^n}{n!}+\dots</math> | |||
The recursion formula for its coefficients is | |||
<math>\dfrac{\alpha_{n+1}}{\alpha_n}= \left( \dfrac{a+n}{c+n} \right) \dfrac{1}{n+1}.</math> | |||
Comparing this to what we obtained earlier, we find that <math>\Phi(\xi)</math> is | |||
<math>\Phi(\xi)=A_1F_1(i\lambda,1,ik\xi),\!</math> | |||
where <math>A\!</math> is a c-number. | |||
The full wave function is thus | |||
<math>\psi(r,z)=A_1F_1(i\lambda,1,ik(r-z))e^{ikz}.\!</math> | |||
Now we should look at the limit where z is taken to go to infinity, and our confluent hypergeometric function is rewritten as: | Now we should look at the limit where z is taken to go to infinity, and our confluent hypergeometric function is rewritten as: |
Revision as of 13:11, 12 January 2014
We now consider the scattering of an electron from the Coulomb potential. This problem is important because it is relevant to the famous scattering experiment by Rutherford that showed that the atomic nucleus only makes up a very small fraction of the total size of an atom. We will first use the Born approximation to find the cross section for a Yukawa potential,
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 V(r)=V_0\frac{e^{-\alpha r}}{r},}
which reduces to the Coulomb potential 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 \alpha\to 0.} We will then find the exact cross section for the Coulomb potential, starting with the Schrödinger equation.
Scattering From a Yukawa Potential in the Born Approximation

Here we compare the differential cross sections for the Yukawa, Gaussian and exponential potential. We can see that for small 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 q\alpha \!} all differential cross sections are similar. The differential cross section of the exponential potential decreases faster than the others and that of the Yukawa potential decreases more slowly than the others.
As alluded to earlier, let us first consider the Yukawa potential,
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 V(r)=V_0\frac{e^{-\alpha r}}{r}.}
Within the Born approximation, the scattering amplitude 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 \begin{align} f(\theta) &= -\frac{m}{2\pi\hbar^2}2\pi\int_0^\infty dr'\,r'^2 \frac{e^{-i|\mathbf{k}'-\mathbf{k}|r'}-e^{i|\mathbf{k}-\mathbf{k}'|r'}}{-i|\mathbf{k}'-\mathbf{k}|r'}\frac{V_0e^{-\alpha r'}}{r'} \\ &= -\frac{2mV_0}{\hbar^2}\int_0^\infty dr' r'^2 \frac{\sin\left(|\mathbf{k}'-\mathbf{k}|r'\right)}{|\mathbf{k}'-\mathbf{k}|r'}\frac{e^{-\alpha r'}}{r'} \\ &= -\frac{2mV_0}{\hbar^2}\frac{1}{\left(\left|\mathbf{k}'-\mathbf{k}\right|\right)^2+\alpha^2} \end{align} }
For the elastic scattering 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| \mathbf{k} \right| = \left| \mathbf{k'} \right| = k \!} . Therefore,
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|\mathbf{k}'-\mathbf{k}\right|)^2 = 2k\sin\frac{\theta}{2} }
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 f = -\frac{2mV_0}{\hbar^2}\frac{1}{4k^2\sin^2\frac{\theta}{2}+\alpha^2}. }
The differential cross section is then
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}= \left( \frac{2mV_0}{\hbar^2} \right) ^2 \left( \frac{1}{4k^2\sin^2\frac{\theta}{2}+\alpha^2}\right)^2.}
Exact Coulomb 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 \phi\!} represents rotation about the z-axis, 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 \xi\!} represents the parabolas with their vertices at a minimum, 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 \eta\!} represents parabolas with their vertices at a maximum.
When we are considering scattering due to the Coulomb potential, we can not neglect the effect of this potential at large distances because it is only a 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{1}{r}} potential.
We will work in parabolic coordinates, which are related to Cartesian coordinates 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 \xi=\sqrt{x^2+y^2+z^2}-z}
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 \eta=\sqrt{x^2+y^2+z^2}+z}
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=\tan^{-1}\left (\frac{y}{x}\right )}
The Schrödinger equation in parabolic coordinates 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 \left[ - \frac{\hbar^2}{2\mu}\frac{4}{\xi+\eta} \left( \frac{\partial}{\partial \xi}\xi\frac{\partial}{\partial \xi}+\frac{\partial}{\partial \eta}\eta\frac{\partial}{\partial \eta}+\frac{\xi+\eta}{4\xi\eta}\frac{\partial^2}{\partial \phi^2} \right) -\frac{2Ze^2}{\xi+\eta} \right] \psi=\frac{\hbar^2k^2}{2\mu}\psi.}
Recall that, for a spherically symmetric potential, the scattering amplitude 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 \theta\!} only; therefore, we will seek solutions that are 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 will look for solution of the form,
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(\xi,\eta,\phi)=e^{\frac{i}{2}k(\eta-\xi)}\Phi(\xi)=e^{ikz}\Phi(r-z).}
The Schrödinger equation then becomes
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[ \xi\frac{\partial^2}{\partial \xi^2}+(1-ik\xi)\frac{\partial}{\partial \xi}+\frac{Ze^2\mu}{k\hbar^2}k \right]\Phi(\xi)=0.}
We now assume a series solution,
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(\xi)=\sum_{n=0}^\infty a_n \xi^n.}
The recursion relation for the coefficients is then
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{a_{n+1}}{a_{n}}=\frac{in-\lambda}{(n+1)^2}k=\frac{n+i\lambda}{(n+1)^2}ik,}
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 \lambda=\frac{Ze^2\mu}{k\hbar^2}.}
Recall that the confluent hypergeometric function 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 {}_1F_1} 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 _1F_1(a,c,z)=1+\frac{a}{c}z+\frac{a(a+1)}{c(c+1)}\frac{z^2}{2!}+\dots+\frac{a(a+1)\dots(a+n-1)}{c(c+1)\dots(c+n-1)}\frac{z^n}{n!}+\dots}
The recursion formula for its coefficients 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 \dfrac{\alpha_{n+1}}{\alpha_n}= \left( \dfrac{a+n}{c+n} \right) \dfrac{1}{n+1}.}
Comparing this to what we obtained earlier, we find 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 \Phi(\xi)} 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 \Phi(\xi)=A_1F_1(i\lambda,1,ik\xi),\!}
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 A\!} is a c-number.
The full wave function is 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 \psi(r,z)=A_1F_1(i\lambda,1,ik(r-z))e^{ikz}.\!}
Now we should look at the limit where z is taken to go to infinity, and our confluent hypergeometric function is rewritten 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 _1F_1(a,c,z)=\frac{\Gamma(c)(-z)^a}{\Gamma(c-a)} \left[ 1-\frac{a(a-c+1)}{z} \right] +\frac{\Gamma(c)}{\Gamma(a)}e^zz^{a-c}}
Now we can use this to rewrite our equation 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 \Phi\!} 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 \xi\!} :
- 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(\xi)=Ae^{\frac{-\pi}{2}\lambda} \left[ \frac{e^{i\lambda \ln(k\xi)}}{\Gamma(1-i\lambda)} \left( 1-\frac{\lambda^2}{ik\xi} \right) +\frac{i\lambda}{ik\xi}\frac{e^{ik\xi+i\lambda \ln(k\xi)}}{\Gamma(1+i\lambda)} \right] }
Rewriting our wavefunction 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\!} :
- 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(r,\theta)=\frac{Ce^{\frac{-\pi\lambda}{2}}}{\Gamma(1-i\lambda)} \left[ \left( 1-\frac{\lambda^2}{2ikr}\frac{1}{\sin^2\frac{\theta}{2}} \right) e^{ikz}e^{i\lambda \ln(k(r-z))}+\frac{f(\theta)}{r}e^{ikr+i\lambda \ln(2kr)} \right] }
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 f(\theta)\!} is the following:
- 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(\theta)=\frac{\lambda\Gamma(1-i\lambda)}{2k\Gamma(1+i\lambda)}\frac{1}{\sin^2\frac{\theta}{2}}e^{i\lambda \ln(\sin^2\frac{\theta}{2})}=\frac{\lambda\Gamma(1-i\lambda)}{2k\Gamma(1+i\lambda)} \left( \sin^2\frac{\theta}{2} \right) ^{i\lambda-1}}
We can then get our differential cross section from that by squaring it:
- 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(\theta)|^2=\frac{\lambda^2}{4k^2\sin^4(\frac{\theta}{2})}=\frac{(Ze^2)^2}{16E^2}\frac{1}{\sin^4(\frac{\theta}{2})}}
If we normalize the wavefunction to give unit flux at large distances, we must take the following for the constant 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 C\!} :
- 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 C=\sqrt{\frac{\mu}{\hbar k}}\Gamma(1-i\lambda)e^{\frac{\pi\lambda}{2}}}
So the wavefunction at large distances is given by the following:
- 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(0)|^2=|C|^2=\frac{\mu}{\hbar k}|\Gamma(1-i\lambda)|^2e^{\pi\lambda}}
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 |\Gamma(1-i\lambda)|^2=\frac{2\pi|\lambda|e^{-\pi\lambda}}{e^{2\pi\lambda}-1}}
Plugging this in for our wavefunction squared:
- 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(0)|^2=\frac{2\pi|\lambda|}{\frac{\hbar k}{\mu}|1-e^{2\pi\lambda}|}}
Now let's use the following quantity to represent the velocity:
- 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{\hbar k}{\mu}=v}
For small incident velocities, we can write:
- 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(0)|^2=\frac{2\pi Ze^2}{\hbar^2 v^2}}
- 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(0)|^2=\frac{2\pi Ze^2}{\hbar^2 v^2}e^{-\frac{2\pi Ze^2}{\hbar v}}}
where the first equation is for an attractive Coulomb potential, and the second equation is for a repulsive Coulomb potential.
The factor , the exponent of , is known as the Gamow factor. The Gamow Factor or Gamow-Sommerfeld Factor, named after its discoverer George Gamow, is a probability factor for two nuclear particles' chance of overcoming the Coulomb barrier in order to undergo nuclear reactions. For example, nuclear fusion. By classical physics, there is almost no possibility for protons to fuse by crossing each other's Coulomb barrier, but when George Gamow instead applied quantum mechanics to the problem, he found that there was a significant chance for the fusion due to tunneling.
The Gamow factor for the Coulomb potential can also be obtained via the WKB method. Recall that
- and
For a broad and high barrier, , and
- .
In this case
Using our previous potential and we get
So, as above, the transmission coefficient for small v is