Radioactive Equilibrium

In physics of nuclear decays, a radioactive equilibrium exists when a radioactive nuclide is decaying at the same rate at which it is being produced. The disintegrating nucleus is usually referred to as the parent nucleus and the nucleus remaining after the event as the daughter nucleus. The daughter nucleus can either be stable or radioactive. If it is radioactive, then it decays into a daughter nucleus and so on. Thus, each radioactive parent nucleus can initiate a series of decays, with each decay-product having its own characteristic decay constant.

Concentration of daughter nuclei in the radioactive equilibrium depends primarily on proportions of half-lives (or decay constants) of parent and daughter nuclei. Since the production rate and decay rate are equal, the number of atoms present remains constant over time. In any case, a radioactive equilibrium is not established immediately, but it only takes place after a transition period. This period is of the order of few half-lifes of the longest-lived nucleus in the decay chain. In case of radioactive decay chains, a radioactive equilibrium may be established between each member of the decay chain.

radioactive equilibrium - equation

As was written, proportionality of half-lives is a key parameter, which determines type of radioactive equilibrium:

  • Radioactive equilibrium is not established when a half-life of the parent nucleus is shorter than a half-life of the daughter nucleus. In this case the production rate and decay rate of certain member of decay chain cannot be equal.
  • Secular radioactive equilibrium exists when the parent nucleus has an extremely long half-life. This type of equilibrium is particularly important in nature. Over the 4.5 billion years of the Earth’s history, especially uranium 238, uranium 235 and thorium 232 and members of their decay chains have reached radioactive equilibria between the parent nucleus and the various descendants.
  • Transient radioactive equilibrium exists when a half-life of the parent nucleus is longer than a half-life of the daughter nucleus. In this case, the parent nuclide and the daughter nuclide decay at essentially the same rate.

Secular Equilibrium

Secular radioactive equilibrium exists when the parent nucleus has an extremely long half-life. Secular equilibrium is typical for natural radioactive series, such as the thorium series or the uranium series. For the uranium series with uranium-238 (with a half-life of 4.47 billion years), where all of the elements in the chain are in secular equilibrium, each of the descendants has built up to an equilibrium amount and all decay at the rate set by the original parent. The only exception is the final stable element (lead-206) on the end of the chain. Its number of atoms is constantly increasing. In any case, a radioactive equilibrium is not established immediately, but it only takes place after a transition period. This period is of the order of few half-lifes of the longest-lived nucleus in the decay chain (234U for the uranium series; 231Pa for the actinium series). In case of radioactive decay chains, a radioactive equilibrium may be established between each member of the decay chain.

Secular Equilibrium

As can be seen,  the secular equilibrium is particularly important in nature. Over the 4.5 billion years of the Earth’s history, especially uranium 238, uranium 235 and thorium 232 and members of their decay chains have reached radioactive equilibria between the parent nucleus and the various descendants. The half-lives of all their descendants are all extremely variable, and it is difficult to represent a range of timescales going from individual seconds to millions of years. Since daughter radioisotopes have different half-lives then secular equilibrium is reached after some time. In the long decay chain for a naturally radioactive element, such as uranium-238, where all of the elements in the chain are in secular equilibrium, each of the descendants has built up to an equilibrium amount and all decay at the rate set by the original parent. If and when equilibrium is achieved, each successive daughter isotope is present in direct proportion to its half-life (or to its decay constant).

radioactive equilibrium - equation

Since its activity is inversely proportional to its half-life, each nuclide in the decay chain finally contributes as many individual transformations as the head of the chain. Radioactive cascade in the secular equilibrium significantly influences radioactivity (disintegrations per second) of natural samples and natural materials. All the descendants are present, at least transiently, in any natural sample, whether metal, compound, or mineral. For example, pure uranium-238 is weakly radioactive (proportional to its long half-life), but a uranium ore is about 13 times more radioactive (approximately 13 descendants in the decay chain) than the pure uranium-238 metal because of its daughter isotopes (e.g. radon, radium etc.) it contains. Not only are unstable radium isotopes significant radioactivity emitters, but as the next stage in the decay chain they also generate radon, a heavy, inert, naturally occurring radioactive gas. Moreover the decay heat of uranium and its decay products (e.g. radon, radium etc.) contributes to heating of Earth’s core. Together with thorium and potassium-40 in the Earth’s mantle is thought that these elements are the main source of heat that keeps the Earth’s core liquid.

The secular equilibrium can occasionally be disrupted when one of the intermediary nuclei leaves the sample (e.g. radon released from soil) where its ancestors are confined. These local disruptions are important to consider in the use of dating techniques. Disruption of the secular equilibrium also influences radioactivity of fresh nuclear fuel. Most of PWRs use the uranium fuel, which is in the form of uranium dioxide. But this uranium dioxide must be chemically purified and most of its decay chain is not present in the fresh nuclear fuel.

Radioactive Equilibrium of Uranium-234

Isotope of uranium-234 is a member of the uranium series. This isotope has the half-life of only 2.46×105 years and therefore it do not belong to primordial nuclides (unlike 235U and 238U). On the other hand, this isotope is still present in the Earth’s crust, but this is due to the fact 234U is an indirect decay product of 238U. 238U decays via alpha decay into 234U. 234U decays via alpha decay into 230Th, except very small fraction (on the order of ppm) of nuclei which decays by spontaneous fission.

In a natural sample of uranium, these nuclei are present in the unalterable proportions of the radioactive equilibrium of the 238U filiation at a ratio of one atom of 234U for about 18 500 nuclei of 238U. As a result of this equilibrium these two isotopes (238U and 234U) contribute equally to the radioactivity of natural uranium.

Transient Radioactive Equilibrium

The transient radioactive equilibrium exists when a half-life of the parent nucleus is longer than a half-life of the daughter nucleus, but the concentration of parent nuclei significantly decreases in time. In this case, the parent nuclide and the daughter nuclide may decay at essentially the same rate, but both concentrations of nuclides decreases as the concentration of parent nuclei decreases. Contrary to secular equilibrium, the half-life of the daughter nuclei is not negligible compared to parent’s half-life.

Transient Equilibrium - Radioactive Equilibrium

An example of this type of compound decay process is a Technetium-99m generator producing technetium-99m for nuclear medicine diagnostic procedures from molybdenum-99. Technetium-99m’s short half-life of 6 hours makes storage impossible and would make transport very expensive. Therefere, for medical purposes, molybdenum-99 is used to produce technetium-99m. These two isotopes are in the transient equilibrium. The decay constant for molybdenum-99 is considerably smaller than the decay constant for technetium-99m. Although the decay constant for molybdenum-99 is smaller, the actual rate of decay is initially larger than that of molybdenum-99 because of the great difference in their initial concentrations. As the concentration of the daughter increases, the rate of decay of the daughter will approach and eventually match the decay rate of the parent. When this occurs, they are said to be in the transient equilibrium. In case of  Technetium-99m generator, transient equilibrium occurs after about four half-lives. Today, technetium-99m is the most utilized element in nuclear medicine and is employed in a wide variety of nuclear medicine imaging studies.

Also the transient equilibrium can occasionally be disrupted when one of the intermediary nuclei leaves the sample where its ancestors are confined.

Transient Radioactive Equilibrium with Source – Example

A special example of radioactive equilibrium are concentrations of iodine-135 and xenon-135 in a nuclear reactor, but in this case, the xenon burnup must be taken into account. Note that, in this special case, the half-life of the parent nucleus is shorter than the half-life of the daughter nucleus. The production and removal of xenon can be characterized by the following differential equations:

xenon equilibrium - iodine equilibriumWhen the rate of production of iodine equals the rate of removal of iodine, equilibrium exists. This equilibrium also known as “xenon 135 reservoir”, since all of this iodine must undergo a decay to xenon. In equilibrium, the iodine concentration remains constant and is designated NI(eq). The following equation for the iodine equilibrium concentration can be determined from the preceding equation by setting the dNI/dt =0. Since the equilibrium iodine concentration is proportional to the fission reaction rate, it is also proportional to reactor power level.

When the rate of production of xenon 135 equals the rate of removal, equilibrium exists also for xenon. The xenon concentration remains constant and is designated NXe(eq). The following equation (1) for the xenon equilibrium concentration can be determined from the preceding equation by setting the dNXe/dt =0. For xenon 135 to be in equilibrium, iodine 135 must also be in equilibrium. Substituting the expression for equilibrium iodine 135 concentration into the equation for equilibrium xenon (1) results in the following (2).

Xenon Worth - different power levelsFrom this equation it can be seen that the equilibrium value for xenon 135 increases as power increases, because the numerator is proportional to the fission reaction rate. But the thermal flux is also in the denominator. Therefore, as the thermal flux exceeds some value, the xenon burnup begins to dominate, and at approximately 1015 neutrons.cm-2.s-1, the xenon-135 concentration approaches a limiting value. The equilibrium iodine 135 and xenon 135 concentrations as a function of neutron flux are illustrated in the following figure.

Bateman Equations

In physics, the Bateman equations are a set of first-order differential equations, which describe the time evolution of nuclide concentrations undergoing serial or linear decay chain. The model was formulated by Ernest Rutherford in 1905 and the analytical solution for the case of radioactive decay in a linear chain was provided by Harry Bateman in 1910. This model can be also used in nuclear depletion codes to solve nuclear transmutation and decay problems.

For example, ORIGEN is a computer code system for calculating the buildup, decay, and processing of radioactive materials. ORIGEN uses a matrix exponential method to solve a large system of coupled, linear, first-order ordinary differential equations (similar to the Bateman equations) with constant coefficients.

The Bateman equations for radioactive decay case of n – nuclide series in linear chain describing nuclide concentrations are as follows:

Bateman Equations

Bateman Equations for Nuclear Transmutation

Fuel Depletion - Isotopic Changes
Isotopic changes of 4% uranium-235 fuel as a function of fuel burnup.

As was written, this model can be also used in nuclear depletion codes to solve nuclear transmutation and decay problems as well. In case of transmutation the decay constants that govern Bateman equations for a decay case are substituted by transmutation constants. By the transmutation constant λi,j we understand probability of the ith nuclide production per time unit from the jth nuclide destruction, as a result of nuclear interaction with the whole spectrum of interacting particles or due to the natural nuclear decay.

These equations are usually used for the exact evolution of isotopic changes in the nuclear fuel during fuel depletion. Fuel depletion is usually modelled mathematically as a set of differential equations known as evolution equations.

Special Reference: Jerzy Cetnar, General solution of Bateman equations for nuclear transmutations. Annals of Nuclear Energy 33 (2006). January 2006.

References:

Radiation Protection:

  1. Knoll, Glenn F., Radiation Detection and Measurement 4th Edition, Wiley, 8/2010. ISBN-13: 978-0470131480.
  2. Stabin, Michael G., Radiation Protection and Dosimetry: An Introduction to Health Physics, Springer, 10/2010. ISBN-13: 978-1441923912.
  3. Martin, James E., Physics for Radiation Protection 3rd Edition, Wiley-VCH, 4/2013. ISBN-13: 978-3527411764.
  4. U.S.NRC, NUCLEAR REACTOR CONCEPTS
  5. U.S. Department of Energy, Nuclear Physics and Reactor Theory. DOE Fundamentals Handbook, Volume 1 and 2. January 1993.

Nuclear and Reactor Physics:

  1. J. R. Lamarsh, Introduction to Nuclear Reactor Theory, 2nd ed., Addison-Wesley, Reading, MA (1983).
  2. J. R. Lamarsh, A. J. Baratta, Introduction to Nuclear Engineering, 3d ed., Prentice-Hall, 2001, ISBN: 0-201-82498-1.
  3. W. M. Stacey, Nuclear Reactor Physics, John Wiley & Sons, 2001, ISBN: 0- 471-39127-1.
  4. Glasstone, Sesonske. Nuclear Reactor Engineering: Reactor Systems Engineering, Springer; 4th edition, 1994, ISBN: 978-0412985317
  5. W.S.C. Williams. Nuclear and Particle Physics. Clarendon Press; 1 edition, 1991, ISBN: 978-0198520467
  6. G.R.Keepin. Physics of Nuclear Kinetics. Addison-Wesley Pub. Co; 1st edition, 1965
  7. Robert Reed Burn, Introduction to Nuclear Reactor Operation, 1988.
  8. U.S. Department of Energy, Nuclear Physics and Reactor Theory. DOE Fundamentals Handbook, Volume 1 and 2. January 1993.
  9. Paul Reuss, Neutron Physics. EDP Sciences, 2008. ISBN: 978-2759800414.

See above:

Radioactive Decay