Hover over the menu bar to pick a physics animation.

Crystal Target Scattering Angular Distribution

This animation is similar to the Rutherford scattering experiment. The Rutherford scattering experiment (RSE) is famous because it established that atoms of a solid have a very compact nucleus. In the case of the RSE, the solid was a very thin gold foil and the incident particles were doubly ionized helium atoms (nucleii) which are also called alpha (`alpha`) particles. The molecule weight of gold is 197 and that of helium is 4. Therefore, if almost all the mass of the gold atom is concentrated in its small nucleus, when the alpha particle hits the gold nucleus, only a small fraction of its energy will be transferred to the nucleus. If the alpha particle hits the gold nucleus dead center, then the alpha particle will bounce back 180 degrees from its original direction. It turned out that that is what happened. It led to Rutherford making the now famous comment:

"It was quite the most incredible event that has ever happened to me in my life. It was almost as incredible as if you fired a 15-inch shell at a piece of tissue paper and it came back and hit you."

This result would not have been the case if the nucleus were not extremely compact. It turns out that the nucleus mean radius is of the order `r_(n)=10^(-14) ` meters. Gold has a face centered cubic (FCC) crystal structure where the cube edge is `a_(gold)=4.07cdot10^(-10)` meters. So viewed from the side, as seen in Canvas 0, the spacing between each gold atom occupies an area of 1/2 the area of a square of edge length `a_(gold)` so each gold atom occupies an area `A_(gold)=(16.56/2)cdot10^(-20) meters^2` (see Spreadsheet 0) . This means that it will be extremely unlikely that a small alpha particle would have a hard collision with a helium nucleus. In fact, the probability, P, is about

`P=(pi r_n^2)/A_(gold)=3.5 cdot 10^(-9)`

and yet a significant fraction of the alpha particles were deflected at other than small angles. Part of the reason for that is that the alpha particle did not have to actually contact the nucleus of gold since both have fairly large positive charge. The charge of gold is `Z_(gold)=79` while that of the alpha particles is `Z_(alpha)=2` where the Z values represent the number of oppositely signed electron charges, `e`, where `e=1.6cdot10^(-19)` coulombs.

With that introduction, we are now ready to work out the details of this animation. The illustration in Canvas 0 shows only a single layer of gold nucleii. An actual gold foil will have at least `n_L=200` layers so the probability of a direct hit is increased by `n_L` or to `P=0.7 cdot 10^(-6)` if `n_L=200`.

Spreadsheet 0: Gold Foil Target Spreadsheet
Canvas 1: Diagram Showing Impact Parameter
Impact Radius 50
Sphere Radius 50
Canvas 0: Diagram of Gold Crystal Geometry
Number of Unit Cell Columns or Rows
Canvas 2: Particle Scattering Animation
Beam-Target Opacity: 0.5
Single Step
Canvas 3: Plots of Scatter Vs Angle .

Equations for the Animation and Plots

Below I show that, when the target particle is much more massive than the beam particle as is the case for alpha particles incident on gold nucleii, the change in angle of the scattered particle is just twice the incident angle. The incident angle is defined as the angle between the incident velocity and the vector between centers of the beam and target particles when they have just achieved contact. See a diagram of these angles in Canvas 4. Then, using equations similar to bouncing off of a flat wall, the final beam particle velocity vector `bb v_f` becomes:

`bb v_f=bb v_i-2 hat bb"u" ( bb v_i cdot hat bb"u") `

where `bb v_i` is the incident velocity and `hat bb"u"` is a unit vector pointing from the target particle to the beam particle. Note that this equation conserves energy

`(bb v_f cdot bb v_f) =(bb v_i cdot bb v_i)`
` -4 (hat bb"u" cdot bb v_i) ( bb v_i cdot hat bb"u")`
`+4 (bb v_i cdot hat bb"u" ) (bb v_i cdot hat bb"u") =(bb v_i cdot bb v_i)`

` Also this equation is correct even when the beam particle makes multiple scatterings from the target particles even when the initial velocity is not along the x axis.
This equation gives the final velocity for each time a beam particle contacts a target particle. Then the beam particle position is iterated using this new velocity.

`x_p=x_p+v_xdt`
`y_p=y_p+v_ydt`

where `dt` is the iteration time increment. When the particle reaches the circular chamber wall, it is stopped. Then particle numbers that have reached the wall are placed in bins Vs their angle on the wall. A plot of bin occupation Vs angle is shown in Canvas 3. The same bin accumulation data is shown in Canvas 2 in a circular plot where the highest bin count is shown at the center of the scatter chamber. The learner has the option of adjusting the target opacity. This depends on both the target particle radius, `r_t`, and particle spacing `s_t` and the beam particle radius, `r_b` via the expression:

`s_t=2.5/(opacity)*(r_t+r_b)`;

where the 2.5 factor is just a parameter that seems to make the ratio of forward scatter to backscatter agree with the numerical results.

Discussion of the Animation Results

First I should say that this is a 2D animation while the Rutherford experiment was a 3D experiment. The 2D here means that the target nucleii are really cylinders along the z axis and perpenducilar to the x and y axes. Also the "gold" nucleii that we use here are much larger compared to their spacing than the nucleii in a real gold foil. This is necessary in order to show at least some backscattering in a finite time. We expect, with target opacity of 50%, to have about half of the beam particles go right through the target without being deflected. As you can see by looking at Canvas 3 after the Run has is complete, that is fairly reasonably valid.
Canvas 4: Diagram Showing Incident and Scattering Angles
Impact Radius 50

Deriving the Equations for Final Velocities After Collision

In the following equations, please be aware that bold letters denote vectors which have either (x,y) or (x,y,z) components depending on whether we are considering 2D or 3D calculations. For hard sphere collisions one must consider both conservation of momentum and energy. If the collision was mediated by a potential between the particles, one would only have to consider conservation of kinetic and potential energy. Conservation of momentum requires that the final momentum be the same as the initial momentum. For a collision between particles 1 and 2:

`m_1 bb"v"_{1f}+m_2 bb"v"_{2f}=m_1 bb"v"_{1i}+m_2 bb"v"_{2i} (1)`

where the subscript f stands for final value and the subscript i stands for initial value and the arrowed bold letter ` bb"v"` stands for the velocity vector. The conservation of kinetic energy equation is as follows:

`m_1 bb"v"_{1f} cdot bb"v"_{1f} +m_2 bb"v"_{2f} cdot bb"v"_{2f}`

`=m_1 bb"v"_{1i} cdot bb"v"_{1i} +m_2 bb"v"_{2i} cdot bb"v"_{2i} (2)`

where the `cdot` between the velocity vectors stands for the inner or dot product of the two vectors. When we combine equations 1 and 2 we get equations that relate the components of the velocity vectors.

The most important concept in these calculations is that of the unit vector ` hat bb"u"` between the centers of the two particles that are colliding

` hat bb"u"=[(x_2-x_1)hat bb"x"+(y_2-y_1)hat bb"y"+(z_2-z_1)hat bb"z"]/r_{12}`

where the x, y, and z that have hats are unit vectors along their respective axes and

`r_{12}=sqrt[(x_2-x_1)^2+(y_2-y_1)^2+(z_2-z_1)^2]`

Note that the unit vector points from sphere 1 to sphere 2 so, if `x_2` is greater than `x_1`, the unit vector's x component is positive.
Applying equation 1, we know that the changes of momentum of the spheres are equal and opposite and are along the unit vector `hat bb"u"` so we can write:

`m_1 delta bb"v"_1=M delta v hat bb"u"=-m_2 delta bb"v"_2 (3)`

where `M` has units of mass and `M` and the scalar speed increment `delta v` are still to be determined. Now we may re-write equation 2 using equation 3 and obtain:

`((m_1 bb"v"_1+M delta v hat bb"u") cdot (m_1 bb"v"_1+M delta v hat bb"u"))/(m_1)+ ((m_2 bb"v"_1-M delta vhat bb"u") cdot (m_2 bb"v"_2-M delta vhat bb"u"))/(m_2) =m_1 bb"v"_{1} cdot bb"v"_{1} +m_2 bb"v"_{2} cdot bb"v"_{2}(3)`

After cancellations of the terms on its right side, equation 4 simplifies to:

`M delta v=(2m_1m_2)/(m_1+m_2)( bb"v"_1- bb"v"_2) cdot hat bb"u"`

We can now make the identifications:

`M=(2m_1m_2)/(m_1+m_2)`

and

`delta v=( bb"v"_1- bb"v"_2) cdot hat bb"u"`

where M is known as the "reduced mass" Finally using equation 3 again we can write the expressions for the final velocity vectors:

` bb"v"_{1f}= bb"v"_{1i}+M/m_1 [( bb"v"_1- bb"v"_2) cdot hat bb"u"] hat bb"u" (5)`
` bb"v"_{2f}= bb"v"_{2i}-M/m_2 [( bb"v"_1- bb"v"_2) cdot hat bb"u"] hat bb"u"`

When `m_2/m_1 gtgt 1` then `M/m_1` becomes 2. We use equation 5 with `M/m_1=2` to compute the scatter direction of the beam particles.
As always, the learner is encouraged to press F12 (Windows) and select "Sources" to examine the HTML and javascript code that was used to create this page.