NuclearProperties Forest NucPhys I

From New IAC Wiki
Revision as of 03:25, 8 April 2009 by Oborn (talk | contribs) (→‎Fermi Momentum)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

Nuclear Properties

Quantum Chromodynamics (QCD) is the fundamental quantum field theory within the standard model that is used to describe the Strong interaction of the fundamental particles known as quarks and gluons, the constituents of a nucleon, in terms of their color. While the Strong force acts directly on the elementary quark and gluon particles, a residual of the force is observed acting between nucleons within an atomic nucleus that is referred to as the nuclear force. The degrees of freedom needed to describe the elementary particles within an average Atomic nucleus with A=50 makes a solution difficult. As a result, a phenomenological approach is used (i.e.: models) to describe nuclear physics processes. Below are some of the static properties of a nucleus that a quantum field theory or model would need to predict.

  1. nuclear charge and radius
  2. mass and binding energy
  3. angular momentum and parity
  4. magnetic dipole and electric quadrupole moments
  5. excited energy levels

These properties are explored further in the sections below.

Nuclear Charge and Radius

Once you know the Isotope then you know how many protons are in the nucleus of interest and therefore the charge.

The interest however is in how that charge is distributed inside the nucleus.

The density of nucleons within the nucleus tends to be uniform over a short distance and then rapidly goes to zero. There are two quantities which are used to characterize nuclear size.

  1. mean radius: The radius of the nucleus in which its density is half of its central value
  2. skin thickness: the distance over which the density of the nucleus drops from its max to its min.

Because the nucleus is charged one could also quantify the nuclear shape in term of

  1. electromagnetic multipole moments: charge = 1st moment, mag dipole = 2nd moment, electric quadrupole= 3rd moment

[math]\alpha[/math] probe

One of the first ways to measure the nuclear radius was to look for the failure of Rutherford scattering. Since Rutherford scattering only considers the coulomb repulsion of two charged object, one could assert that the physics changes when measurements of the scattered particles fail to obey the Rutherford prediction.

An alpha particle scattered from a Gold nucleus deviates from Rutherford scattering at incident alpha particle energies of about 27 MeV. If you assume that this is the energy at which the incident alpha particle starts to probe the inside of the nucleus then you can estimate the size of a uniformly charged spherical nucleus with a hard boundary as:

[math]\frac{1}{2} mv^2 = \frac{1}{4 \pi \epsilon_0} \frac{Q_1 Q_2}{r}[/math]

[math]\Rightarrow r= \frac{1}{4 \pi \epsilon_0}\frac{Q_1 Q_2}{\frac{1}{2} mv^2} = 8.3 fm[/math]

The modern value for the gold nucleus is 7.2 fm.

RutherfordScatteringBreakdown.jpg

optical measure of radius

Quantum measure of radius:

One direct method used to determine the size of an object in to alter the size of a probe until interference patterns emerge (ie [math]\lambda \lt D[/math]).

CircularDiskDiffraction.jpg DiffractionOpaqueSphere.jpg

[math]\sin(\theta) = \frac{1.22 \lambda}{D}[/math]
Note
The pattern broadens as [math]\lambda \rightarrow D[/math].

By virtue particle wave duality, the elastic scattering of an electron from a nucleus can behave in a similar fashion to the scattering of light by an opaque target.


Figure 3.1 from Krane's book.


DiffractionElectronsByO16andC12.jpg

Although 3-D objects reveal diffraction patterns that are similar to those from a two dimensional disk, such comparisons are really only rough estimates. As we saw from the above, the scattering of photons from O-16 revealed a diffraction pattern. A similar diffraction pattern can be seen if you elastically scatter electrons from a Pb-208 nucleus.

Distribution of Nuclear Charge

For our current interest we would like a means to measure the distribution of charge in the nucleus of an atom. This can be accomplished by elastically scattering a probe off of the nucleus which is sensitive to charge. Elastically scattering an electron off of a nucleus is one such probe. A description of this can be found using Fermi's Golden Rule where the transition amplitude

[math]| M_{i,f}| ^2 \equiv \int \psi_f^{*}(\vec{r}) H_{int} \psi_i(\vec{r}) dr^3[/math]

is the main term in Fermi's Golden Rule which describes the interaction that takes place.

In our current example we have an incident electron plane wave in some initial state which gets elastically scattered by some charge distribution (a nucleus) to a final plane wave state by means of a coulomb interaction.

[math]\psi_i(\vec{r}) = e^{i \vec{k}_i \cdot \vec{r}}[/math]
[math]\psi_f^{*}(\vec{r}) = e^{-i \vec{k}_f \cdot \vec{r}}[/math]
[math]H_{int} = \frac{-1}{4 \pi \epsilon_0} \frac{Q}{|\vec{r}- \vec{R}|} = \frac{-ze^2}{4 \pi \epsilon_0} \int \int e^{i \vec{q} \cdot \vec{r}} \frac{\rho(\vec{R})}{|\vec{r} - \vec{R}|} dV dV^{\prime}[/math]

where

[math]dV^{\prime}=[/math] integral over the volume of the nucleus
[math]dV =[/math] integral over all space
[math]\rho(\vec{R})=[/math] charge density as a function of R from the center of the nucleus
[math]\vec{r} =[/math] distance of probe from the center of the nucleus
[math]\vec{q} = \vec{k}_i - \vec{k}_f =[/math] momentum transfered to the charge distribution.


Let

[math]\vec{s} = \vec{r} - \vec{R}[/math]
[math]\mu = \frac{\vec{q}\cdot \vec{s}}{|\vec{q}||\vec{s}|} = \cos (\theta)[/math]
[math]\Rightarrow dV = 2 \pi s^2 ds d\mu[/math]

Then

[math]H_{int} = \frac{-ze^2}{4 \pi \epsilon_0} \int \int e^{i \vec{q} \cdot (\vec{s}+ \vec{R})} \frac{\rho(\vec{R})}{|s|} 2 \pi s^2 ds d \mu dV^{\prime}[/math]
[math] = \frac{-ze^2}{4 \pi \epsilon_0} \int e^{i \vec{q} \cdot\vec{s}} 2 \pi s ds d \mu \int e^{i \vec{q} \cdot \vec{R}}\rho(\vec{R})dV^{\prime}[/math]
[math] I = \int e^{i qs \mu} 2 \pi s ds d \mu =[/math] Some integral
[math] F(q) \equiv \int e^{i \vec{q} \cdot \vec{R}}\rho(\vec{R})dV^{\prime}=[/math] Physics of interest = Fourier transform of [math]\rho[/math]


The inverse Fourier Transform would be

[math]\rho(\vec{R}) = \frac{1}{(2 \pi)^3} \int d^3q F(q) e^{-i \vec{q} \cdot \vec{r}}[/math]


To find the charge density you would measure the transition rate as a function of the momentum transfer. The charge density will then be the inverse Fourier transform of that data.

The data below are the result of such measurements.

RadialChargeDistribution C-12 Pb-208.jpg

F(q) general form
[math]F(q) = \int e^{i \vec{q} \cdot \vec{r}} \rho(r) dV[/math]
[math]= \int e^{i q r \cos(\theta)} \rho(r) r^2 \sin(\theta) d \theta d \phi dr[/math]
[math]= \frac{2 \pi}{q} \int \rho(r) r dr \int_{0}^{\pi} e^{i q r \cos(\theta)} q r sin(\theta) d \theta d \phi[/math]

Let

[math]u = qr \cos(\theta)[/math]
[math]du = -qr \sin(\theta) d \theta[/math]
[math]F(q) = \frac{2 \pi}{q} \int \rho(r) r dr \int_{qr}^{-qr} e^{iu} (-du)[/math]
[math]= \frac{2 \pi}{q} \int \rho(r) r dr \left . \frac{e^{iu}}{i} \right |_{-qr}^{qr}[/math]
[math]= \frac{2 \pi}{q} \int \rho(r) r dr \frac{e^{iu} - e^{-iqr}}{i} [/math]
[math]= \frac{4 \pi}{q} \int \rho(r) r dr \sin(qr)[/math]
[math]= \frac{4 \pi}{q} \int \sin(qr) \rho(r) r dr [/math]

Charge Radius

A measure of [math]F(q)[/math] can also be used to extract the charge radius

[math]F(q) = \frac{4 \pi}{q} \int \sin(qr) \rho(r) r dr [/math]
[math]\sin(\theta) = \theta - \frac{\theta^3}{3!} + ...[/math]


[math]F(q)= \frac{4 \pi}{q} \int \left ( qr - \frac{(qr)^3}{3!} \right ) \rho(r) r dr[/math]
[math]= \frac{4 \pi}{q} \left [ \int q \rho(r) r^2 dr - \int \frac{(qr)^3}{3!} \rho(r) r dr\right ][/math]
[math]\int_0^{\infty} q \rho(r) r^2 dr = \frac{1}{4 \pi}[/math] : Density is normalized
[math]\int_0^{\infty} \frac{(qr)^3}{3!} \rho(r) r dr = \frac{q^3}{6} \int_0^{\infty} r^2 \rho(r) r^2 dr = \frac{q^3}{6} \frac{\lt r^2\gt }{4 \pi}[/math]
[math]\Rightarrow F(q) = 1 - \frac{q^2}{6} \lt r^2\gt [/math]

Take the derivative of F(q) with respect to Q then you get a slope of the function. The slope of this function near the origin(q=0) tells you the mean charge radius squared. Making these electron scattering measurements on several different nuclei has revealed a cubed root relationship between the radius and the atomic number A.

[math] R = R_o A^{\frac{1}{3}} \approx 1.23 A^{\frac{1}{3}} fm[/math]

NuclearChargeRadiusFromElectronScattering.jpg

K X-ray isotope shift

Another method employed to measure the charge radius of a nucleus involves a measurement of the K X-rays produced by 2 different isotopes of the same atom.

A "K X-ray" is the photon given off when and electron undergoes a transition from the 2P orbit to the 1S.

Let

[math]E_K(A) = E_{2P}(A) - E_{1S}(A)[/math] = K X-ray energy from isotope [math]A[/math] (ie: A=208, Z=82, Pb-208)
[math]E_K(A^{\prime}) = E_{2P}(A^{\prime}) - E_{1S}(A^{\prime})[/math] = K X-ray from isotope A^{\prime} (ie: [math]A^{\prime}[/math] = 207, Z=82 , Pb-207)

Assume

[math]E_{2P}(A) = E_{2P}(A^{\prime})[/math]

This assumes that the energy difference [math]= E_{1S}(A) - E_{1S}(A^{\prime})[/math] is a lot bigger than the difference for the 2P electrons. Since the energy difference we are investigating is caused by the 1S electron spending some time inside the nucleus sampling the nucleus' charge distribution, one can assert that the higher energy 2P electrons state spends substantially less time sampling the nuclear charge distribution.


[math]\Delta E \equiv E_K(A) -E_K(A^{\prime}) = \left [ E_{2P}(A) - E_{1S}(A) \right ] - \left [ E_{2P}(A^{\prime}) - E_{1S}(A^{\prime}) \right ][/math]
[math]= E_{1S}(A^{\prime}) - E_{1S}(A)[/math]


Let

[math]= E_{1S}(A^{\prime}) = E_{1s_p}(A^{\prime})- \Delta E_{1S}(A^{\prime})[/math]

where

[math]E_{1s_p} =[/math] energy level for an electron with a point nucleus

similarly

[math]= E_{1S}(A) = E_{1s_p}(A)- \Delta E_{1S}(A)[/math]


[math]\Rightarrow \Delta E = \Delta E_{1S}(A^{\prime}) -\Delta E_{1S}(A)[/math]
[math]\Delta E_{1S}(A) = [/math] change in the electron energy eigen value when the point like assumption for the nucleus is removed


Electron energy correction

Use the coulomb potential to calculate the change [math]\Delta E_{1S}(A)[/math] in the electron energy eigen value when the nucleus is given a finite instead of being point like.

An electron which has a finite probability of existing within a finite size nucleus, of radius [math]R[/math], will feel a different coulomb potential when it is inside than outside. This possibility of having an electron spend some time inside the nucleus effectively changes its energy eigenvalue by an amount [math]\Delta E[/math]. This change in energy is

[math]\Delta E = \lt V(r \le R)\gt - \lt V(r \gt R)\gt [/math]
[math]= \int_0^R \Psi^*_N \hat{V}(r \le R) \Psi_N dV - \int_0^R \Psi^*_N \hat{V}(r \gt R) \Psi_N dV[/math]

When the electron is outside the nucleus

[math]V(r \gt R) = - \frac{Ze^2}{4 \pi \epsilon_0 r}[/math]
[math]V(r\le R) = - \frac{Ze^2}{4 \pi \epsilon_0 } \left ( \frac{3}{2R} - \frac{r^2}{2R^3}\right )[/math]

substitution:

[math]\Delta E =- \frac{Ze^2}{4 \pi \epsilon_0 } \int_0^R \Psi^*_N \left ( \frac{3}{2R} - \frac{r^2}{2R^3} - \frac{1}{r}\right ) \Psi_N dV[/math]

To make the calculation easy lets assume that the electrons wave function for a finite size nucleus is the same as the wave function for the point like nucleus. The radial wave function for Hydrogen is

[math]\Rightarrow \Psi_N = 2 \left ( \frac{Z}{a_0} \right )^{3/2} e^{\frac{-Zr}{a_0}}[/math]
[math]a_0 \equiv[/math] Bohr radius

Because the operator has no angular dependence we only need to do the radial part of the integral. The angular part of the integral is normalized to unity. So

[math]dV \rightarrow r^2 dr[/math]
[math]\Delta E =- \frac{Ze^2}{4 \pi \epsilon_0 } \int_0^R 2 \left ( \frac{Z}{a_0} \right )^{3/2} e^{\frac{-Zr}{a_0}}\left ( \frac{3}{2R} - \frac{r^2}{2R^3} - \frac{1}{r}\right ) 2 \left ( \frac{Z}{a_0} \right )^{3/2} e^{\frac{-Zr}{a_0}} r^2dr [/math]
[math]= - \frac{Z^4e^2}{\pi \epsilon_0 a_0^3 } \int_0^R e^{\frac{-2Zr}{a_0}} \left ( \frac{3}{2R} - \frac{r^2}{2R^3} - \frac{1}{r}\right ) r^2 dr [/math]

Since

[math]\frac{R}{a_0} \sim 10^{-15} \Rightarrow e^{\frac{-2Zr}{a_0}} \sim e^{0} = 1[/math]


[math]\Delta E = - \frac{Z^4e^2}{\pi \epsilon_0 a_0^3 } \int_0^R \left ( \frac{3}{2R} - \frac{r^2}{2R^3} - \frac{1}{r}\right ) r^2dr[/math]
[math]= -\frac{Z^4e^2}{\pi \epsilon_0 a_0^3 } \left . \left ( \frac{3}{2R} \frac{r^3}{3} - \frac{1}{2R^3} \frac{r^5}{5} - \frac{r^2}{2}\right ) \right |_0^R[/math]
[math]= -\frac{Z^4e^2}{\pi \epsilon_0 a_0^3 } \left ( \frac{R^2}{2} - \frac{R^2}{10} - \frac{R^2}{2}\right )[/math]
[math]= -\frac{Z^4e^2}{\pi \epsilon_0 a_0^3 } \left ( - \frac{R^2}{10} \right )[/math]
[math]\Delta E = \frac{Z^4e^2}{\pi \epsilon_0 a_0^3 } \frac{R^2}{10} = \frac{Z^4e^2}{4\pi \epsilon_0 a_0^3 } \frac{2R^2}{5} [/math]
[math]\Rightarrow \Delta E_{1S}(A) = \frac{Z^4e^2}{4\pi \epsilon_0 } \frac{2R^2}{5a_0^3}[/math]


The K X-ray isotope shift is:

[math] \Delta E = \Delta E_{1S}(A^{\prime}) -\Delta E_{1S}(A)[/math]
[math]= \frac{Z^4e^2}{4\pi \epsilon_0 } \frac{2}{5a_0^3} \left [ \left ( R^{\prime}\right )^2 - R^2 \right ][/math]

From electron scattering experiments we found that

[math] R = R_o A^{\frac{1}{3}} \approx 1.23 A^{\frac{1}{3}} fm[/math]
[math]\Rightarrow \Delta E \equiv E_K(A) -E_K(A^{\prime}) = \frac{Z^4e^2}{4\pi \epsilon_0 } \frac{2R_o^2}{5a_0^3} \left [ \left ( A^{\prime}\right)^{2/3} - A^{2/3} \right ][/math]

If you measure the isotope shift you can infer [math]R_o[/math]

K-X-ray E shift Hg.jpg

Notice
The slope of the fit shown (solid line for even nuclei) is about 10% higher than what is predicted

Using wave functions which account for relativistic effect moves the prediction into agreement with experiment.

The experiment has been improved by using muonic atoms. A muonic atom has a muon in place of the electron. Because of its higher mass the muon has a tighter "orbit" , the muon isn't really orbiting the nucleus, and has a higher probability of being inside the nucleus thereby being more sensitive to the charge distribution. The X-rays emmited by muons cascading down the energy levels are in the MeV range while electrons are in the keV range.

Muon experiments have measured [math]R_0 = 1.25[/math].

NuclearRadiusMuonKxrayEnergyShiftMeasurement.jpg

Coulomb Energy

Yet another method to determine the radius of a nucleus considers the impact of the coulomb force on the binding energy of Mirror nuclei. Mirror nuclei are two nuclei with the same number of nucleons but different number of protons and neutrons such that the number of protons (Z) in one nucleus is equal to the number of neutrons(N) in the other nucleus. For example He-3 (Z=2,N=1) and H-3 (Z=1, N=2). The difference in the coulomb binding energy of such mirror nuclei can be used to determine the radius of the nuclei.

If we assume the nucleus is a uniformly charge non-conducting sphere, then the "self energy", or the energy needed to assemble the charge distribution can be found using the work energy theorem for conservative forces.


[math]W = -\Delta U =[/math] Work -Energy Theorem :work needed to assemble a charge distribution which corresponds to a potential energy U = qV

The potential energy for the assembled charge distribution would be

[math]U = qV = \int V dq = \int V \rho 4 \pi r^2 dr = \int \frac{1}{4 \pi \epsilon_0 } \frac{Q}{r} \rho 4 \pi r^2 dr[/math]
[math]= \frac{1}{ \epsilon_0 } \int \rho \frac{4 \pi r^3}{3} \rho r dr = \frac{4 \pi\rho^2}{ 3\epsilon_0 } \int_0^R r^5 dr[/math]
[math]\Rightarrow U = \frac{4 \pi\rho^2}{ 3\epsilon_0 } \frac{R^5}{5} = \frac{4 \pi}{ 3\epsilon_0 } \left ( \frac{Q}{4 \pi R^3/3}\right )^2 \frac{R^5}{5} =\frac{1}{ 4 \pi \epsilon_0 } \frac{3}{5} \frac{Q^2}{R}[/math]

If you compare the coulomb energy difference between two mirror nuclei with one having Z protons and one have (Z-1) protons then the coulomb energy difference would be

[math]\Delta E_C = \frac{3}{5} \frac{e^2}{4 \pi \epsilon_0 R} \left [ Z^2 - \left ( Z-1 \right )^2 \right ] =\frac{3}{5} \frac{e^2}{4 \pi \epsilon_0 } \left [ \frac{2Z-1}{R}\right ] [/math]

For mirror nuclei : Z = N-1

[math]2Z - 1 = Z + (Z-1) = Z+N = A[/math] of the mirror nucleus

and

[math]R = R_o A^{1/3}[/math]
[math]\Delta E_C =\frac{3}{5} \frac{e^2}{4 \pi \epsilon_0 } \frac{A}{R_0 A^{1/3}} = \frac{3}{5} \frac{e^2}{4 \pi \epsilon_0 } \frac{A^{2/3}}{R_0 }= mA^{2/3}[/math]

One way to measure [math]\Delta E_C[/math] is to detect the [math]\beta[/math] decay of a mirror nucleus in which the proton changes into a neutron and emits a positron. The Max positron energy observed is [math]\Delta E_C[/math]. If you plot [math]\Delta E_C[/math] as a function of [math]A^{2/3}[/math] then the slope is proportional to [math]R_0[/math].

The graph below is from Krane's book

R 0 FromMirrorNucleiCoulombEnergyDiff.jpg

Matter radius

The above electron probe is a fine way to measure charge distribution in a nucleus. A nucleus however does contain uncharged nucleons known as neutrons. Neutrons do have a charge distribution, more positive core to a more negative surface charge. In order to measure the distribution of nuclear matter you will need to use a probe which depends more on the strong force and less on the electromagnetic.


Mesons (pions) are used in a manner similar to the use of muons by looking at the X-rays emmited as they cascade down the energy levels. A meson interacts with the nucleus through both the Strong and E&M forces.

Summary

There are 3 ways to characterize the nuclear shape

  1. mean radius:the density of nucleons within the nucleus drops to half its mean value
  2. skin thickness: The distance over which the density drops from a max to a min
  3. electromagnetic dipole moment: 1st moment = charge, 2nd moment = mag dipole, 3rd moment = electric quadrupole

Nuclear Binding Energy and Mass

Binding energy

The binding energy of a nucleus is defined as the mass difference between the constituents of a nucleus and the nucleus.

[math]B(Z,A) = \left [ Z m(H-1) + N M_n - m(X-A)\right ]c^2[/math]

If you make a plot of B/A -vs- A for ground state nuclei you would see something similar to the curve below.

BindingEnergyPerNucleon-vs-A.jpgRadialChargeDistribution C-12 Pb-208.jpg


The above plots may be described by some closed functional form which

Volume Term

As seen in the above graph, the ratio of B/A is almost a constant 8 MeV until you start to get to low A (A<9). As we saw previously, the density of nucleon inside the nucleus appears to be constant until you get to the edge.

[math]\frac{ \mbox{Num. of Nucleons}}{\mbox{Volume}} = \frac{A}{4 \pi r^3/3} = \mbox{Constant} \Rightarrow R =R_0 A^{1/3}[/math]

or in other words the atomic number [math]A[/math] is proportional to the volume of the nucleus. At first glance you can imaging the the nuclear force acting on a group of A nucleons would be proportional to [math]A^2[/math] thinking any given nucleon is experiencing a force from all the other nucleons. Based on the above graph you expect the binding energy to be proportional the the Atomic number, suggesting that the strong force is so short range that nucleons really only see their nearest neighbor.

So our first term to fit the above curve has the form

[math]B(Z,A) = \alpha_V A[/math]

Surface Term

The next parameter for the fit function is called the surface term [math](\alpha_S)[/math]. Experiments have shown that as you get near the edge of a nucleus the density of nucleons changes. The nucleons on the surface of the nucleus have fewer nearest neighbors and they are farther apart (Less bound) from each other. As a result the surface term should reduce the overestimate from the volume term. Furthermore, if the volume term is proportional to A then the radius is proportional to [math]A^{1/3}[/math] so the surface (Area) term should be proportional to [math]\left( A^{1/3} \right)^2[/math]

[math]B(Z,A) = \alpha_V A - \alpha_S A^{2/3}[/math]

Coulomb Term

While the nuclear force is trying to bind the nucleons, the coulomb force is trying to push protons apart thereby making them less bound. As we saw previously, the self energy of a uniformly charged sphere suggests that

[math] U = \frac{1}{ 4 \pi \epsilon_0 } \frac{3}{5} \frac{Q^2}{R} = \frac{1}{ 4 \pi \epsilon_0 } \frac{3}{5} \frac{-Z(Z-1)}{R_0 A^{1/3}}[/math]

which leads a reduction in the binding energy

[math]B(Z,A) = \alpha_V A - \alpha_S A^{2/3} - \alpha_C \frac{Z(Z-1)}{A^{1/3}}[/math]

Asymmetry Term

Stable isotopes tend to have Z [math]\approx \frac{A}{2}[/math],

If Z is not half of A then you the nucleus is less bound (more unstable). So your expect a functional dependence like

[math]A-2Z[/math]

to appear such that

[math]B(Z,A) = \alpha_V A - \alpha_S A^{2/3} - \alpha_C \frac{Z(Z-1)}{A^{1/3}} - \alpha_{sym} \frac{(A-2Z)^2}{A} [/math]
Note
You may see the asymmetry term sometimes written as
[math]\frac{(N-Z)^2}{A} : A=N+Z[/math]

Pairing Term

it is also observed that nucleons like to pair up ( spin couple) inside the nucleus. If you had an even number of nucleons you would be able to spin couple such that you can occupy a lower energy state. Whereas if you had an odd number of nucleons, you would have some less bound un-paired nucleons within the nucleus.

[math]B(Z,A) = \alpha_V A - \alpha_S A^{2/3} - \alpha_C \frac{Z(Z-1)}{A^{1/3}} - \alpha_{sym} \frac{(A-2Z)^2}{A} - \delta[/math]

one functional form used for [math]\delta[/math] is

[math]\delta = \alpha_p \frac{1}{A^{1/2}}[/math]
[math]\alpha_p =\left \{ {0 \;\;\;\; \mbox {A odd} \atop \alpha_p \;\;\;\; \mbox {A even} } \right .[/math]

The sign of [math]\alpha_p[/math] changes depending on if [math]Z[/math] and [math]N[/math] are both even or both odd.

Semiemirical mass formula constants

Parameter Krane Wapstra Rohlf
[math]\alpha_V[/math] 15.5 14.1 15.75
[math]\alpha_S[/math] 16.8 13 17.8
[math]\alpha_C[/math] 0.72 0.595 0.711
[math]\alpha_{sym}[/math] 23 19 23.7
[math]\alpha_p[/math] [math]\pm[/math]34 [math]\pm[/math] 33.5 [math]\mp[/math] 11.18

Except for the Rohlf measurements the "+" sign is used when N & Z are even and the "-" is used when N & Z are both odd integers.

Ref:

Wapstra: Atomic Masses of Nuclides, A. H. Wapstra, Springer, 1958

Rohlf: Modern Physics from a to Z0, James William Rohlf, Wiley, 1994

Liquid Drop Model

Semi-empirical (Weizacker) mass formula

The semi-empirical mass formula takes the definition of binding energy B(Z,A), solves for the mass of the nucleus, and then inserts the binding energy fit equation.

take the binding energy equation

[math]B(Z,A) = \left [ Z m({1 \atop 1 }H) + N M_n - m({A \atop Z }X_{N})\right ]c^2[/math]


and solve for the nucleus mass

[math]m({A \atop Z }X_{N}) = \left [ Z m({1 \atop 1 }H) + N M_n - \frac{B(Z,A)}{c^2}\right ][/math]

Examples

BE(26,59)

[math]B({59 \atop 26} Fe_{33}) = \left [ 26(1.007825 u) + 33 (1.008665 u) - 58.93488 u\right ] \frac{931.5 MeV}{u} = 516.53 MeV[/math] [math]\frac{B({59 \atop 26} Fe_{33})}{59} = 8.755 MeV[/math]

Nucleon Separation Energies

Neutron Separation Energy=[math]S_n = B({A \atop Z} X_{N})-B({{A-1} \atop Z} X_{N-1})[/math]= Energy needed to liberate a neutron from the Nucleus Proton Separation Energy=[math]S_p = B({A \atop Z} X_{N})-B({{A-1} \atop {Z-1}} X_{N})[/math]= Energy needed to liberate a proton from the Nucleus

Mass decay chains

A=125 decay chain

BetaDecayChainA-125.png

Isotope B(MeV) N Z t_{1/2}
Praseodymium (Pr-59) 1020.67 even odd 3.3 sec (inverse Beta)
Cerium (Ce-58) 1030.1 odd even 9 sec
Lanthanum (La-57) 1038.21 even odd 76 sec (inverse Beta)
Barium (Ba-56) 1044.63 odd even 3.5 min
Cesium (Cs-55) 1049.98 even odd 45 min (inverse Beta)
Xenon (Xe-54) 1053.86 odd even 16.9 hours
Iodine (I-53) 1056.29 even odd 59.41 days inverse beta decay
Tellurium (Te-52) 1057.26 odd even Stable
Antimony (Sb-51) 1057.28 even odd 2.76 years
Tin (Sn-50) 1055.7 odd even 9.64 days
Indium (In-49) 1051.06 even odd 12.2 sec
Cadmium (Cd-48) 1044.72 odd even 0.65 sec Beta decay
A=128 decay chain
Isotope B(MeV) N Z t_{1/2}
Lanthanum (La-57) 1067.29 odd odd 5 min
Barium (Ba-56) 1074.73 even even 2.43 days
Cesium (Cs-55) 1076.03 odd odd 3.66 min electron capture
Xenon (Xe-54) 1080.74 even even Stable
Iodine (I-53) 1079.41 odd odd 25 min Beta decay
Tellurium (Te-52) 1081.44 even even 10^24 years Stable for our lifetime
Antimony (Sb-51) 1077.84 odd odd 10.4 min
Tin (Sn-50) 1077.35 even even 59.07 min
Indium (In-49) 1069.15 odd odd 0.72,0.84 sec
Cadmium (Cd-48) 1062.87 even even 0.34 sec

Angular Momentum and Parity

Angular Momentum (j)

Consider the case of a proton orbiting around a nuclear shell

[math]\ell[/math] = orbital angular momentum
[math]s[/math] = spin
[math]\vec{j} = \vec{\ell} + \vec{s}[/math] = total angular momentum
[math]\vec{I} = \sum_i \vec{j}_i[/math] = theoretical spin of the nucleus

Often though a single valence nucleon determines the total angular momentum of a nucleus. In some cases the total angular momentum is determine by two valence nucleons.

At other times the total angular momentum is given by the combination

[math]\vec{I} = \vec{j}_{valence} + \vec{j}_{core}[/math]

Magneton

Bohr Magnton
[math]\frac{e \hbar}{2 m_e} = 5.7884 \times 10^{-5} \frac{eV}{T} = \mu_B[/math]
Nuclear Magnton
[math]\frac{e \hbar}{2 m_p} = 3.1525 \times 10^{-8} \frac{eV}{T} = \mu_N[/math]

Nuclear Magnetic moment

Magnetic Dipole Moment

Classical magnetic dipole moment was found by looking at the torque on a current carrying wire of area A carrying a current i and immersed in an external magetic field

[math]\tau = iAB = \mu B[/math]

If you consider a particle of charge e moving in a circular orbit of radius r then

[math]i = \frac{e}{2 \pi r/v}[/math]


[math]|\mu| = \frac{e \hbar}{2 m_p} \ell = g_{\ell} \mu_N \ell[/math] : nuclear model
[math]g_{\ell} =\left \{ {1 \;\;\;\; proton \atop 0 \;\;\;\; neutron} \right .[/math]

Bulk effects involving magnetism are usually determined by electrons (atomic magnetism). Only in special cases (NMR/MRI) can you see nuclear magnetism.

Spin Magnetic moment
[math]|\mu| = g_{s} s \mu_N [/math]
[math]g_{s} =\left \{ {5.5856912 \pm 0.0000022 \;\;\;\; proton \atop -3.8260837 \pm 0.0000018 \;\;\;\; neutron} \right .[/math]

Dirac Equation [math]\Rightarrow[/math] [math]g_s = 2[/math] for point particles.

[math]g_s[/math] (electron) = 2.0023

[math]\Rightarrow[/math] Perhaps the proton and neutron are not point particles!

Fermi Momentum

Fermi Momentum/Energy represents the amount of energy a nucleon (proton or neutron) has when confined to a nucleus. One way to determine this energy is to approximate the nucleus as a 3-D infinite potential well in the form of a cube of side L.

From the 1-D infinite potential well problem you have the eigenfunctions (wave functions) and eigenvalues:

[math]\psi_n(x) = \sqrt{\frac{2}{L}} \sin\left(\frac{n \pi x}{L}\right) [/math]
[math]E_n = \frac{n^2\hbar^2 \pi ^2}{2mL^2} = \frac{n^2 h^2}{8mL^2} [/math]

with [math]n[/math] being a positive integer of quantization.


The three dimensions are independent and thus separable allowing the total wave function and energy eigenvalues to be written as:


[math]\psi_n(x) = A \sin\left(\frac{n_x \pi x}{L}\right)\sin\left(\frac{n_y \pi y}{L}\right)\sin\left(\frac{n_z \pi z}{L}\right) [/math]
[math]E_n = \frac{\hbar^2 \pi ^2}{2mL^2} |\vec{n}|^2[/math]

where

[math]|\vec{n}|^2 = n_x^2 + n_y^2 + n_z^2[/math]

Now add N non-interacting, spin 1/2 fermions into this box.

Note
An atom with atomic number [math]A[/math] identified by the Chemical symbol [math]X[/math] with [math]Z[/math] protons has [math]A-Z[/math] neutrons is denoted as [math]{A \atop Z} X[/math]
in this case N = A = number of nucleons (fermions) in the system (box).


The nucleons in the system occupy occupy the lowest energy state ([math]E_{f}[/math]) which forms a sphere in momentum space of radius [math]n_f[/math].

In this ground state the number of nucleons ([math]N[/math]) is equal to the number of states within the sphere of radius [math]n_f[/math] in momentum space. In momentum space, the quantum numbers (n_x, n_y, and n_z) are all positive integers. This means that only 1/8 of the sphere is occupied. Each state can hold 2 spin 1/2 nucleons and still obey the Pauli principle.


[math]\frac{1}{8}[/math] Volume of sphere = [math]\frac{1}{8} \frac{4}{3} \pi n_f^2 = \frac{N}{2}[/math]

or

[math]n_f = \left ( \frac{3 N}{\pi} \right)^{1/3}[/math]

so the Fermi energy is given by

[math]E_f =\frac{\hbar^2 \pi^2}{2m L^2} n_f^2 [/math]
[math] = \frac{\hbar^2 \pi^2}{2m L^2} \left( \frac{3 N}{\pi} \right)^{2/3}[/math]


Which results in a relationship between the fermi energy and the number of particles per volume (when we replace L2 with V2/3):

[math]E_f = \frac{\hbar^2}{2m} \left( \frac{3 \pi^2 N}{V} \right)^{2/3} \,[/math]

Fermi Momentum

The total energy of a Fermi sphere of [math]A[/math] fermions is given by

[math]E_t = {\int_0}^{A} E_f(N) dN = {\int_0}^{A} \frac{\hbar^2 \pi^2}{2m L^2} \left( \frac{3 N}{\pi} \right)^{\frac{2}{3}} dN[/math]
[math] = \frac{3^{\frac{2}{3}} \pi^{\frac{4}{3}} \hbar^2}{2m L^2} {\int_0}^{A} N^{\frac{2}{3}} dN = \frac{3^{\frac{2}{3}} \pi^{\frac{4}{3}} \hbar^2}{2m L^2} \left( \frac{3}{5} A^{\frac{5}{3}} \right) = \frac{3^{\frac{5}{3}} \pi^{\frac{4}{3}} \hbar^2}{10 m L^2} A^{\frac{5}{3}}[/math]
[math]= \frac{\hbar^2}{2m_p} \left [ \frac{3A}{5} \left( \frac{3 \pi^2 A}{V} \right )^{2/3} \right ][/math]
[math]=\frac{3 A E_f}{5}[/math]
[math]L^2 = V^{\frac{2}{3}}[/math]

The Nuclear size/radius of the nucleus is roughly:

[math]R = \left(1.25 \times 10^{-15} \right) \times A^{1/3}[/math] m
where A is the number of nucleons.

The number density of nucleons in a nucleus is therefore:

[math]V = \frac{4}{3} \pi R^3 \approx 8.2 \times 10^{-45} m^{3} A[/math]

So the fermi energy of a nucleus is about:

[math]m= 938.2 MeV/c^22[/math]
[math]\hbar = 6.582 \times 10^{-16} eV s[/math]


[math]E_t = = \frac{\hbar^2}{2m_p} \left [ \frac{3A}{5} \left( \frac{3 \pi^2 A}{V} \right )^{2/3} \right ][/math]
[math]= \frac{(6.6 \times 10^{-22} MeV s)^2}{2 \times 938 MeV/c^2} \frac{(3 \times 10^8 m/s)^2}{c^2}\left [ \frac{3A}{5} \left( \frac{3 \pi^2 A}{8.2 \times 10^{-45} m^{-3} A} \right )^{2/3} \right ][/math]
[math]= 30 MeV \times A [/math]

or

[math]E_f = \frac{5 E_{tot}}{3 A} = 50 MeV[/math]


Note: in reality the neutrons and protons are independent systems so I have over counted the total energy when integrating over A.

[math]E_f = 25 MeV[/math]


Forest_NucPhys_I