Singularity in Rutherford cross section

which means at [tex]\vartheta[/tex]=0 the differential cross section is infinite, which is ok.
My question is, given that the differential cross section is proportional to the probability per unit solid angle [tex]\frac{dP}{d\Omega}[/tex], which is proportional to the expected rate of scattered particles, why/how does the expected rate go to infinity at [tex]\vartheta[/tex]=0?

You are correct. After integrating over solid angle for a deflection greater than θ, the scattering cross section is

σ(>θ) = (πb2/4) cot2(θ/2)

where b= "collision diameter" = zZe2/(½mv2)

For small scattered angles, the impact parameter x** (distance of closest approach) is very large. But the Coulomb field does not extend to infinity, due to screening of the nuclear Coulomb field by atomic electrons. Thus at small angles, the cross section cuts off, and the total cross section is finite.

** x = (b/2) cot(θ/2)

I hope this helps.

Bob S

[added] Very roughly, the maximum impact parameter x should be set to about the "radius" of the 1s atomic electrons, just as electron screening begins to shield the nucleus.

When you do a Rutherford scattering experiment using an alpha source and a gold foil, you will see the counting rate drop by several orders of magnitude as θ increases, following the Rutherford scattering formula. Then at about θ = ~90 degrees the differential cross section starts deviating remarkably from the Rutherford cross section (Coulomb scattering). You can use this to establish the impact parameter of the gold nucleus.

The problem with Bob's answer is that you still have a singularity in the scattering with a bare nucleon. This singularity is real.

What you have to remember is that the Rutherford cross-section is a theoretical calculation. Like all calculations, how well it matches the real world depends on how well the assumptions match the real world.

The assumption that matters is that the interaction between the scattered particle and the scattering particle is instantaneous and depends pretty much only on the closest point of approach of the particles. This is a good approximation for wide angle scatters: a free particle comes in close, has a sudden interaction, and then a free particles moves out.

For interactions where the two particles never get close, this is a poor approximation - and that's exactly where the Rutherford cross-section becomes singular. So the lesson you should take from the singularity is not that you can ignore it because it's difficult to prepare a system that's not under the influence of nearby electrons: it's that you should recognize that this is telling you that (at least) one of the assumptions in the derivation is no longer valid.

Coulomb scattering (the Rutherford scattering model) is for a charged particle scattering in a 1/r2 central force (i.e., the Coulomb field). The calculated trajectory is hyperbolic, and is exact for a simple scattering of a charged particle in a central Coulomb field. As you may realize, this is exactly the same as the scattering of an extra-solar object in the gravitational field of the Sun. The main differences are that for Coulomb scattering of a nucleus, 1) the force may either attractive (electron, negative muon) or repulsive (alpha, proton, positive muon), 2) the force at large distances is shielded (screened) by atomic electrons, and 3) the force is non-Coulomb at small distances (due to strong interactions, radiative interactions, vacuum polarization, etc.).

The infinite cross section for very small-angle scattering is correct for Coulomb scattering off a bare nucleus, because the Coulomb field extends to infinite distances. But for neutral atoms, the Coulomb field is cut off due to screening by atomic electrons. A good electron-screening model is the Thomas-Fermi model, which can be approximated by a cutoff maximum impact parameter of x = 0.529·Z-1/3 x 10-8 cm. Moliere scattering theory also addresses the small scattering angle, large impact parameter regime.

It ca't be singularity simply from mathematical point of view. Your Coulomb wave functions are finite for any angle. But you have singular cross section. What is it? Then try to calculate total cross section or transport, which is extremely important for calculating kinetic characteristics in plasma or semiconductors!

When you have a semiconductor, lets say extrinsic. There you have screening radius Rs. or in angles 1/(kRs), k is a wave vector of your particle(p=hk). then when kRs^2/r >1 then you don't have screening ( r=n^{1/3}, n-concentration)? What are you doing in this case?

You mathematically calculate your cross section incorrectly.

I know that cross section is finite(Not actually the differential cross section itself, but sigma(\te)*sin(\te) remains finite for any angle!), I calculated it.

It ca't be singularity simply from mathematical point of view. Your Coulomb wave functions are finite for any angle. But you have singular cross section. What is it? Then try to calculate total cross section or transport, which is extremely important for calculating kinetic characteristics in plasma or semiconductors!

When you have a semiconductor, lets say extrinsic. There you have screening radius Rs. or in angles 1/(kRs), k is a wave vector of your particle(p=hk). then when kRs^2/r >1 then you don't have screening ( r=n^{1/3}, n-concentration)? What are you doing in this case?

You mathematically calculate your cross section incorrectly.

I know that cross section is finite(Not actually the differential cross section itself, but sigma(\te)*sin(\te) remains finite for any angle!), I calculated it.
.

The infinite cross section for very small-angle scattering is correct for Coulomb scattering off a bare nucleus, because the Coulomb field extends to infinite distances. But for neutral atoms, the Coulomb field is cut off due to screening by atomic electrons. A good electron-screening model is the Thomas-Fermi model, which can be approximated by a cutoff maximum impact parameter of x = 0.529·Z-1/3 x 10-8 cm. Moliere scattering theory also addresses the small scattering angle, large impact parameter regime.

tupos, i think you've completely understood what it is that was asked.

I was curious to why the Rutherford cross section tends to infinity. I think we're all aware that the its not the complete picture, but i wanted to know why it broke down so catastrophically in the small angle region.

The problem with vanadium calculation is that he use asymptotic calculation of differential cross section. If you open any book on quantum mechanics, you'll find out the integral equation for scattering theory. it is a sum of plain wave and integral of \int{g(r-r')\psi(r')V(r')dr'}, you know the green function exp[ik|r-r'|]/|r-r'|, then you expand your green function and obtain asymptotic representation exp[i\vec k \vec r] + f(\te) exp{ikr}/r, from here you make a conclusion that your d\sigma is proportional to |f(\te)|^2!

But if you carefully think for the Coulomb case you can't expand the green function, because it is mathematically incorrect. And your asymptotic representation can't be used. That is why your Rutherford cross section calculated wrong.

Except of this you should use exact solution of Schroedinger equation in the Coulomb field

exp{i \vec k \vec r}F(i\xi,1,ikr(1-cos(\te))), which is a sum of two functions U_1 and U_2,

U_1 has spherical wave asymptotic and U_2 has plain wave.(See Mors and Feshbah for more details). After this you should calculate flux in a usual way

\psi^*\nabla \psi - \psi\nabla\psi* using U_1 function as a \psi in this formula.

And only after this you should calculate the cross section using standard definition

\si = j_{sc} / j_0 * r^2 d\Omega.

If you do all this you find out the your cross section is proportional to function |G_1(\te)|^2*sin(\te). G_1 behaves on small angles like ln(x), (Mors,Feshbah). From here you can see that differential cross section is finite because |ln(x)|^2 * x = 0 when x->0.

From here you can calculate total cross section which is of course not infinite and has the following form

\sigma_total = 2\pi r /k * I(\xi),

\xi here and upper in F is the same and characterized your incident beam, \xi = \alpha /(hv), where v is a velocity v = hk/m, alpha = Ze^2. And I is some function which depends only from \xi.

Maybe now you understand it better. It is a problem of every book of quantum mechanics. People simply calculated this incorrectly for the whole angle region.

Regards.

P.S. asymptotic calculation are correct in angle region \te>\te_0

where \te_0 describe the "near" zone, a zone where you can't distinguish plain and spherical wave from each other.