Nuclear Physics Tutorial 5 - Exponential Law of Decay

 

It is important to understand that radioactive decay is entirely a random and unpredictable process.  If we look at any one nucleus, it might decay in ten seconds or ten million years.  There is no way of telling when the decay will happen, and there is certainly no way of speeding up the process.  Remember that chemical reactions involve the outer shell electrons; radioactive decay involves the nucleus.

 

However when we have many millions of millions of nuclei, we can apply statistical models of probability.

 

The rate of disintegration of any nuclide at any time is directly proportional to the number of atoms left at the time:

This is calculus notation for rate of change in number.

 

The minus sign tells us that N decreases as time increases.

 

There is a constant of proportionality called the radioactive decay constant, which is given the physics code l (greek letter ‘l’, not to be confused with wavelength), and has the units s-1.  So we can write:

The radioactive decay constant is defined as the fraction of the total number of nuclei present that decays per unit time, provided that the time interval is small.  For nuclides with relatively long half-lives, the decay constant will have units “per second”.  However some very unstable nuclides decay in microseconds, so the decay constant would have to be “per nanosecond”.

 

Example

0.25 kg radon-226 emits alpha particles at a measured rate of 9 × 1012 s-1.  What is the decay constant of radium?  (No of atoms in a mole = 6 × 1023)

Work out the number of particles:

 

  1 mole of radon atoms has a mass of 0.226 kg.

 

0.250 × 6 × 1023 = 6.64 × 1023 atoms

0.226

We know that the rate of decay is 9 × 1012 s-1.  So we use DN/Dt = -lN

                        - 9 × 1012 s-1 = -l × 6.64 × 1023

l = 1.36 × 10-11 s-1

(The minus sign indicates a decay)

 

The unit Becquerel (Bq) is often used.  1 Bq = 1 count per second.

 

Question 1

A sample of living material contains carbon 14 with an activity of 260 Bq kg-1.  What is the decay constant? (The fraction that is made of carbon-14 is 1.4 × 10-12)

Answer

 

Over longer periods of time, the relationship above does not hold.  It can be shown by calculus methods (which you are not expected to know), that the decay follows the relationship:

[N – no of nuclei; N0 – original number of nuclei; e – exponential number, 2.718…; l - decay constant (s-1); t – time (s)]

 

This relationship is described as an exponential decay and the graph looks like this:

 

 

We should note the following:

  • The rate of decay represents the number of atoms remaining.  So we can use this graph for representation of the count rate, or ionisation current.

  • The rate of decay is the number of disintegrations per second, or the activity.  It is measured in Becquerels  (Bq).

  • 1 Bq = 1 disintegration per second.

 

Notice that the horizontal axis is calibrated in periods of time in which the decay goes from an initial value to half that value.  This period is called the half-life.  The formal definition of half-life is the time taken for the activity of a sample to decrease to half some initial value.  After 1 half life the activity is 50 %, between 1 and 2 half lives it falls from 50 % to 25 % and so on.  If we have a whole number of half-lives, we can do an iterative calculation of the fraction remaining, i.e. ½, ¼, 1/8, etc.

 

We can relate the half-life to the decay constant:

 

By definition

 

The formula is N = N0e-lt.  Rearranging gives us

 

T1/2 (pronounced "T-half") is the physics code for half-life.

 

 

This will work for any units of time, although we need to be consistent.

 

To work with these relationships, you must be familiar with the concept of logarithms (a way of expressing a number as a power of 10 or e).  You will see logarithms expressed as ‘lg’ or log10, which means log to the base 10, or ‘ln’ or loge, log to the base e.  The number e is the exponential number, 2.718... The latter are known as natural logarithms, which are at the heart of exponential functions.

 

Example

A radiographer has calculated that a patient is to be injected with 1 ´ 1018 atoms of iodine 131 to monitor thyroid activity.  The half-life is 8 days.  Calculate:

(a)    the radioactive decay constant

(b)   the initial activity

(c)    the number of undecayed atoms of iodine 131 after 24 days.

(d)   The total activity after 3 days.

(a) We need to use T1/2 = 0.693

                                               l

                  we need to convert the 8 days into seconds.

 

                        Ž l =    0.693   _ = 1.00 ´ 10-6 s-1

                                    8 ´ 86400

(b) Use DN = - lN = 1.00 ´ 10-6 s-1 ´ 1 ´ 1018 = 1 ´ 1012 Bq

                    Dt

         (c) 24 days is 3 half-lives.  Therefore the number atoms remaining undecayed is 1/8 of the original.

N  =  1 ´ 1018 ø 8 = 1.25 ´ 1017

(d) 3 is not so easy.  We use A = A0e-lt

Ž A = 1 ´ 1012 Bq ´ e-(1.00 ´ 10-6 s-1 ´ 3 ´ 86400s)

 

Ž A = 1 ´ 1012 Bq ´ e-(0.2592)

        =  1 ´ 1012 Bq ´ 0.772 = 7.72 ´ 1011 Bq.

[Alternative method] If you are not so confident in the use of the natural logarithm, you can work out the number of half-lives a particular time gives.  In this case 3 days = 3/8 of a half life.

 

Activity = 1 ´ 1012 ´ 1/23/8 = 1 ´ 1012 ´ 1/1.30 = 7.7 ´ 1011 Bq

 

This is a perfectly valid alternative approach.

 

Question 2

What is the half life of radon-226?  l = 1.36 × 10-11 s-1

Answer

 

The half-life has important implications for the storage of radioactive waste.  Radioactive waste is some of the nastiest muck known to man, so it has to be stored carefully.  Isotopes with a short half-life have intense activities, whereas those with a long half-life have lower activities, but it takes much longer for the activity to decay to a reasonable level.  Either way, it’s not very nice.

 

Even when there have been several half-lives, there remains quite a considerable activity.  If in our example above the activity had decayed to 1/1000th of its original value, we would still have 109 disintegrations per second, which is quite a lot.

 

Another use of the decay equations is in radioactive dating.  Isotopes used are carbon-14, rubidium-87, and hydrogen 3.

 

Question 3

Strontium-90 is a beta emitter.  It is one of the radio-nuclides found in the fall out from an atomic bomb explosion.  It can be absorbed into the bone.  It emits beta particles and has a half life of 28 years.  What is the time needed for the activity to fall to 5 % of the original?

Answer
Question 4

A GM tube placed close to a radium source gives an initial average corrected count rate of 334 s-1

(a)     The GM tube detects 10 % of the radiation.  What is the initial activity?

(b)    Initially there were 1.5 ´ 109 nuclei in the sample.  What is the decay constant?

(c)     What is the half life of the radium in days?

Answer

 

Radioactive Dating

Carbon-14 has a half-life of 5700 years.  It decays by beta minus decay by this equation:

 

Carbon-14 is made by the collision of an energetic neutron with a nitrogen nucleus, which results in the ejection of a proton.  The energetic neutron is formed by cosmic ray activity.  Therefore carbon-14 is being made all the time. The rate of decay is balanced by the rate of formation, so the amount of carbon-14 remains almost constant.

 

The proportion of carbon-14 in living material is very low.  For each gram of carbon, the count-rate is about 140 min-1, 2.33 Bq.  The count rate is therefore 2300 Bq kg-1.  This is the reference point we use for 100 % activity.

 

If we know the activity of the sample, A, and the reference activity, A0, we can use the equation:

 

We can take natural logs to give us:

 

And this rearranges to:

-lt = ln A - ln A0

 

Since ln A0 is bigger than ln A, the right hand side of the equation is negative.  The left hand side is negative, so the two negative signs cancel.  This ensures that the time is positive.  We can only go forward in time.

 

To work out the decay constant, we need to use the equation:

 

In this case the half life is 5700 years, which we have to convert to seconds.  You do know how to do this, don't you?

 

Worked example

In an archaeological dig, remains are found of a warrior who was buried using a wooden boat for a coffin.  A 0.50 gram sample of the boat is found to have an activity of 3800 counts per hour.  What is the age of the boat?  The activity of carbon-14 in living material is 2300 Bq kg-1 and the half-life of carbon-14 is 5700 years.

Answer

Find out the activity of the sample in Becquerel.

A = 3800 ÷ 3600 = 1.056 Bq

 

Convert this to activity per unit mass.

 

A = 1.056 ÷ 0.5 × 10-3 = 2111 Bq kg-1

 

Now we need to work out the decay constant.

l = 0.693 ÷ (5700 × 365 × 24 × 3600) = 3.86 × 10-12 s-1

 

Now use:

-lt = ln A - ln A0

 

-3.86 × 10-12  × t = ln (2111) - ln (2300)

 

-3.86 × 10-12  × t =7.655 - 7.741 = -0.0857

Notice that the minus signs cancel:

t = 0.0857 ÷ 3.86 × 10-12 = 2.2 × 1010 s (700 years)

 

Normally you would leave your answer in seconds.  The answer is given to 2 significant figures as the "activity in living material" data item is to 2 significant figures.

 

The calculation above assumes that the level of carbon-14 is the same in the year 1300 as it is now.  If the cosmic ray activity were less then, the amount of carbon 14 would be reduced.  The warrior's remains would be more recent.

 

 

Rock Dating

Carbon-14 dating is only useful for objects under about 50000 years old.  It's no good for dating rocks.  The isotope potassium-40 is useful.  There are two modes of decay for potassium-40:

  1. Normal beta minus decay;

  2. Electron Capture.

For beta minus, the decay is:

The half life for this is about 1250 million years.

 

For electron capture, the electron is captured from the inner shell:

 

About 1 in 8 decays is of the latter type.  So out of 9 potassium-40 atoms, 8 will decay to calcium, while 1 decays to argon.  The problem with the decay to calcium is that calcium is as common as muck, so we wouldn't be able to tell which calcium atom came from the potassium decay, and which came from other sources.

 

However, Argon is an inert gas, and does not react with any crystals.  Argon atoms remain trapped in the rock, so if we knew how many argon atoms there were, we could work out how many potassium atoms had decayed.  So if we counted 12 × 106 argon atoms, we could say that these resulted from 108 × 106 potassium 40 atoms.  So if we have 500 × 106 potassium 40 atoms left, we can say that we would have started off with 108 × 106 + 500 × 106 = 608 × 106 atoms.

 

Question 5

If potassium-40 has a half-life of 1250 million years, what is the age of the sample of rock discussed in the paragraph above?

Answer

 

Summary

 

Radioactive decays are random.

 

Rate of decay depends on the number of atoms left

 

The probability of any one nucleus decaying in any one second is the decay constant

 

Decay constant is given the code l

 

DN = -lN

Dt

 

Over a longer period of time, decay is exponential.

 

N = No e-lt

 

Half life is the time taken for ½ the remaining atoms to decay

 

T1/2 = 0.693

              l