This site uses cookies. To see our cookie policy and how to control cookies click hereCookies

Probability and differential equations

Some atomic nuclei are energetically unstable and can undergo spontaneous radioactive decay to form a more stable nucleus. The process of radioactive decay is purely random and, in a given fixed time period, a nucleus has a fixed probability of decaying. The process is also memoryless. This means that it doesn't matter whether the nucleus has gone one time period without decay or one hundred time periods, the probability of decay in the next time period remains the same.A time period of particular use is the time for which the probability of decay is 0.5. Such a probability means that over the time period, half of the nuclei will decay and half will not. This time period is called the

The number of nuclei remaining afer \(N\) half-lives, \(A_N\) is given mathematically by $$ A_N = \frac{A_0}{2^N},$$ where \(A_0\) is the initial number. If you are only considering the nuclide that is decaying the above equation also applies to the mass. Care must be taken with mass however! Consider 100 g of polonium-210, \(^{210}Po\). After one half-life (138.38 days) there will be 50 g left. However, there will not be 50 g of its stable decay daughter lead-206, \(^{206}Pb\), as each nucleus of lead has approximately 2% less mass than each nucleus of polonium, leading to there being about 49 g of lead. The missing mass is largely accounted for by the alpha particles (a helium nucleus consisting of two protons and two neutrons) emitted during the decay process.

The interactive element below simulates the radioactive decay of red nuclei to yellow nuclei. Three different half-life values can be selected and line markers are applied at every half-life interval. Notice how in each interval, the number of red nuclei halves, following the equation given above.

The plot in the interactive element above is drawn using time intervals that are much smaller than the half-life so a smooth curve is drawn. This has an impact on the probability of decay value that should be used. To see why, consider the following equation for how the number of nuclei changes over a time period: $$A_{t} = A_0 - pA_0 = A_0 \left( 1-P \; \right),$$ where \(P \;\) is the probability of decay over the time period. If the time period is one half-life then, by definition, \(P = 0.5\). Now consider taking two time steps, each of half a half-life. A first guess is that as the time-step is half as long, the probability of decay should be half as large, so that after the first half time-step the number of nuclei is $$A_{\frac{t}{2}} = A_0 \left( 1-\frac{P}{2} \; \right),$$ and after the second half time-step the final value is $$A_{t} = A_{\frac{t}{2}} \left( 1-\frac{P}{2} \; \right) = A_0 \left( 1-\frac{P}{2} \; \right) ^ 2.$$ Using \(P = 0.5\) in the above will give the wrong answer - it will not yield a reduction by half over the half-life. To achieve this a value of \(P = 0.586\) (3 d.p.) is required. This exercise can be extended to ever smaller sub-divisions of the half-life. When this is done the value of \(P \;\) converges to the natural logarithm of 2 (about 0.693).