Difference between revisions of "Talk:Neutron scattering lengths"

From GISAXS
Jump to: navigation, search
(Origin of the scattering lengths)
(Potential well)
Line 14: Line 14:
 
Consider a neutron of energy <math>E_i</math> interacting with a nucleus, which exhibits an attractive square well of depth <math>-V_0</math> and width <math>2R</math>; where <math>V_0 \gg E_i</math>. The [http://en.wikipedia.org/wiki/Schr%C3%B6dinger_equation Schrödinger equation] is:
 
Consider a neutron of energy <math>E_i</math> interacting with a nucleus, which exhibits an attractive square well of depth <math>-V_0</math> and width <math>2R</math>; where <math>V_0 \gg E_i</math>. The [http://en.wikipedia.org/wiki/Schr%C3%B6dinger_equation Schrödinger equation] is:
 
:<math>
 
:<math>
\left[ - \frac{h^2}{8 \pi^2 m}\nabla^2  + V(r) \right] \psi(r) = E \psi(r)
+
\left[ - \frac{\hbar^2}{2 m}\nabla^2  + V(r) \right] \psi(r) = E \psi(r)
 
</math>
 
</math>
  
Line 21: Line 21:
 
\psi_{s,\mathrm{out}} = \frac{\sin(kr)}{kr} - b \frac{e^{ikr}}{r}
 
\psi_{s,\mathrm{out}} = \frac{\sin(kr)}{kr} - b \frac{e^{ikr}}{r}
 
</math>
 
</math>
where <math>k=\sqrt{2mE_i} 2 \pi/h</math>. Inside the square-well (<math>|r|<R</math>), the potential is <math>V(r)=-V_0</math>, and the solution becomes:
+
where <math>k=\sqrt{2mE_i} /\hbar</math>. Inside the square-well (<math>|r|<R</math>), the potential is <math>V(r)=-V_0</math>, and the solution becomes:
 
:<math>
 
:<math>
 
\psi_{s,\mathrm{in}} = A \frac{\sin(qr)}{qr}
 
\psi_{s,\mathrm{in}} = A \frac{\sin(qr)}{qr}
 
</math>
 
</math>
where <math>q=\sqrt{2m(E_i+V_0)} 2 \pi/h</math>. The two solutions are subject to a continuity boundary condition at <math>|r|=R</math>:
+
where <math>q=\sqrt{2m(E_i+V_0)} /\hbar</math>. The two solutions are subject to a continuity boundary condition at <math>|r|=R</math>:
 
:<math>\begin{alignat}{2}
 
:<math>\begin{alignat}{2}
 
\psi_{s,\mathrm{out}} (r=R) & = \psi_{s,\mathrm{in}} (r=R) \\
 
\psi_{s,\mathrm{out}} (r=R) & = \psi_{s,\mathrm{in}} (r=R) \\

Revision as of 12:17, 6 June 2014

Origin of the scattering lengths

The following description is adapted from Boualem Hammouda's (NCNR) SANS tutorial.

Neutron energy

Consider first the energies of neutrons used in scattering experiments (recall the neutron mass is 1.67×10−27 kg). A thermal neutron (~100°C) would have energy of:

Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle \mathrm {KE} ={\frac {1}{2}}mv^{2}={\frac {3}{2}}kT=7\times 10^{-21}\,\mathrm {J} =48\,\mathrm {meV} }

The velocity of such neutrons is ~3000 m/s, and the momentum is Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle p=mv=5\times 10^{-24}\,\mathrm {kgm/s} } . Finally, the deBroglie wavelength would be:

A cold neutron (~18 K) would have energy of 4×10−22 J = 2 meV, velocity of ~660 m/s, and wavelength of 6 Å.

Potential well

Consider a neutron of energy interacting with a nucleus, which exhibits an attractive square well of depth Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle -V_{0}} and width ; where . The Schrödinger equation is:

Outside of the square-well (Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle |r|>R} ), Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle V(r)=0} , and so the equation is solved as simply:

where . Inside the square-well (), the potential is , and the solution becomes:

where . The two solutions are subject to a continuity boundary condition at Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle |r|=R} :

Note that ; because of the small neutron mass and energy (see above), as well as the small size of a nucleus (femtometers). Therefore:

Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle {\begin{alignedat}{2}\psi _{s,\mathrm {out} }&={\frac {\sin(kR)}{kR}}-b{\frac {e^{ikR}}{r}}\\&\approx 1-b/r\end{alignedat}}}

and so:

And from equating the derivatives:

Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle {\begin{alignedat}{2}1&=A\cos(qR)\\A&={\frac {1}{\cos(qR)}}\end{alignedat}}}

Combining the two results yields:

Failed to parse (Conversion error. Server ("https://wikimedia.org/api/rest_") reported: "Cannot get mml. Server problem."): {\displaystyle {\begin{alignedat}{2}R-b&=\left({\frac {1}{\cos(qR)}}\right){\frac {\sin(qR)}{q}}\\{\frac {R}{R}}-{\frac {b}{R}}&={\frac {\tan(qR)}{q}}{\frac {1}{R}}\\{\frac {b}{R}}&=1-{\frac {\tan(qR)}{qR}}\end{alignedat}}}

This final equation gives a first-order estimate for the scattering length, b, given the radius of the nucleus (R ~ 10−15 m) and the depth of the potential well (V0 ~ MeV).

TBD: Graph

The extreme variation of b/R as a function of qR means that with each nucleon added, the scattering length jumps to a very different value. This also demonstrates why the scattering length can be negative (indicative of negative phase shift). This model should only be taken qualitatively, of course (e.g. we have neglected absorption, as well as the detailed form of the nuclear potential).