Difference between revisions of "Talk:Neutron scattering lengths"

From GISAXS
Jump to: navigation, search
(Potential well)
(Blanked the page)
 
(5 intermediate revisions by the same user not shown)
Line 1: Line 1:
==Origin of the scattering lengths==
 
The following description is adapted from [http://www.ncnr.nist.gov/programs/sans/pdf/polymer_tut.pdf 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<sup>−27</sup> kg). A thermal neutron (~100°C) would have energy of:
 
:<math>
 
\mathrm{KE} = \frac{1}{2} m v^2 = \frac{3}{2} kT = 7 \times 10^{-21} \, \mathrm{J} = 48 \,\mathrm{meV}
 
</math>
 
The velocity of such neutrons is ~3000 m/s, and the momentum is <math>p=mv=5\times10^{-24} \, \mathrm{kgm/s}</math>. Finally, the deBroglie wavelength would be:
 
:<math>
 
\lambda = \frac{h}{p} = 1.3 \, \AA
 
</math>
 
A cold neutron (~18 K) would have energy of 4×10<sup>−22</sup> J = 2 meV, velocity of ~660 m/s, and wavelength of 6 Å.
 
===Potential well===
 
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>
 
\left[ - \frac{h^2}{8 \pi^2 m}\nabla^2  + V(r) \right] \psi(r) = E \psi(r)
 
</math>
 
  
Outside of the square-well (<math>|r|>R</math>), <math>V(r)=0</math>, and so the equation is solved as simply:
 
:<math>
 
\psi_{s,\mathrm{out}} = \frac{\sin(kr)}{kr} - b \frac{e^{ikr}}{r}
 
</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:
 
:<math>
 
\psi_{s,\mathrm{in}} = A \frac{\sin(qr)}{qr}
 
</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>:
 
:<math>\begin{alignat}{2}
 
\psi_{s,\mathrm{out}} (r=R) & = \psi_{s,\mathrm{in}} (r=R) \\
 
\frac{\mathrm{d}}{\mathrm{d}r} \psi_{s,\mathrm{out}} (r=R) & = \frac{\mathrm{d}}{\mathrm{d}r}  \psi_{s,\mathrm{in}} (r=R)
 
\end{alignat}
 
</math>
 
Note that <math>kR = \sqrt{2 m E_i} R 2\pi/h \ll 1</math>; because of the small neutron mass and energy (see above), as well as the small size of a nucleus (femtometers). Therefore:
 
:<math>
 
\begin{alignat}{2}
 
\psi_{s,\mathrm{out}} & = \frac{\sin(kR)}{kR} - b \frac{e^{ikR}}{r} \\
 
    & \approx 1 - b/r
 
\end{alignat}
 
</math>
 
and so:
 
:<math>
 
\begin{alignat}{2}
 
\psi_{s,\mathrm{out}} (r=R) & = \psi_{s,\mathrm{in}} (r=R) \\
 
1 - \frac{b}{R} & = A \frac{\sin(qR)}{qR} \\
 
R - b & = A \frac{\sin(qR)}{q}
 
\end{alignat}
 
</math>
 
And from equating the derivatives:
 
:<math>
 
\begin{alignat}{2}
 
1 & = A \cos(qR) \\
 
A & = \frac{1}{ \cos(qR)  }
 
\end{alignat}
 
</math>
 

Latest revision as of 11:41, 6 June 2014