Physics:Radioactive decay

Short description: Emissions from unstable atomic nuclei

Alpha decay by a nucleus emits an alpha particle of two neutrons and two protons; also known as a helium nucleus.

Radioactive decay (also known as nuclear decay, radioactivity, radioactive disintegration, or nuclear disintegration) is the process by which an unstable atomic nucleus loses energy by radiation. A material containing unstable nuclei is considered radioactive. Three of the most common types of decay are alpha, beta, and gamma decay. The weak force is the mechanism that is responsible for beta decay, while the other two are governed by the electromagnetism and nuclear force.

Radioactive decay is a stochastic (i.e., random) process at the level of single atoms. According to quantum theory, it is impossible to predict when a particular atom will decay, regardless of how long the atom has existed. However, for a significant number of identical atoms, the overall decay rate can be expressed as a decay constant or as half-life. The half-lives of radioactive atoms have a huge range; from nearly instantaneous to far longer than the age of the universe.

The decaying nucleus is called the parent radionuclide (or parent radioisotope), and the process produces at least one daughter nuclide. Except for gamma decay or internal conversion from a nuclear excited state, the decay is a nuclear transmutation resulting in a daughter containing a different number of protons or neutrons (or both). When the number of protons changes, an atom of a different chemical element is created.

There are 28 naturally occurring chemical elements on Earth that are radioactive, consisting of 34 radionuclides (six elements have two different radionuclides) that date before the time of formation of the Solar System. These 34 are known as primordial nuclides. Well-known examples are uranium and thorium, but also included are naturally occurring long-lived radioisotopes, such as potassium-40.

History of discovery

Pierre and Marie Curie in their Paris laboratory, before 1907

Radioactivity was discovered in 1896 by scientists Henri Becquerel and Marie Skłodowska-Curie, while working with phosphorescent materials. These materials glow in the dark after exposure to light, and Becquerel suspected that the glow produced in cathode ray tubes by X-rays might be associated with phosphorescence. He wrapped a photographic plate in black paper and placed various phosphorescent salts on it. All results were negative until he used uranium salts. The uranium salts caused a blackening of the plate in spite of the plate being wrapped in black paper. These radiations were given the name "Becquerel Rays".

It soon became clear that the blackening of the plate had nothing to do with phosphorescence, as the blackening was also produced by non-phosphorescent salts of uranium and by metallic uranium. It became clear from these experiments that there was a form of invisible radiation that could pass through paper and was causing the plate to react as if exposed to light.

At first, it seemed as though the new radiation was similar to the then recently discovered X-rays. Further research by Becquerel, Ernest Rutherford, Paul Villard, Pierre Curie, Marie Curie, and others showed that this form of radioactivity was significantly more complicated. Rutherford was the first to realize that all such elements decay in accordance with the same mathematical exponential formula. Rutherford and his student Frederick Soddy were the first to realize that many decay processes resulted in the transmutation of one element to another. Subsequently, the radioactive displacement law of Fajans and Soddy was formulated to describe the products of alpha and beta decay.

The early researchers also discovered that many other chemical elements, besides uranium, have radioactive isotopes. A systematic search for the total radioactivity in uranium ores also guided Pierre and Marie Curie to isolate two new elements: polonium and radium. Except for the radioactivity of radium, the chemical similarity of radium to barium made these two elements difficult to distinguish.

Marie and Pierre Curie's study of radioactivity is an important factor in science and medicine. After their research on Becquerel's rays led them to the discovery of both radium and polonium, they coined the term "radioactivity" to define the emission of ionizing radiation by some heavy elements. (Later the term was generalized to all elements.) Their research on the penetrating rays in uranium and the discovery of radium launched an era of using radium for the treatment of cancer. Their exploration of radium could be seen as the first peaceful use of nuclear energy and the start of modern nuclear medicine.

Early health dangers

Taking an X-ray image with early Crookes tube apparatus in 1896. The Crookes tube is visible in the centre. The standing man is viewing his hand with a fluoroscope screen; this was a common way of setting up the tube. No precautions against radiation exposure are being taken; its hazards were not known at the time.

The dangers of ionizing radiation due to radioactivity and X-rays were not immediately recognized.

X-rays

The discovery of X‑rays by Wilhelm Röntgen in 1895 led to widespread experimentation by scientists, physicians, and inventors. Many people began recounting stories of burns, hair loss and worse in technical journals as early as 1896. In February of that year, Professor Daniel and Dr. Dudley of Vanderbilt University performed an experiment involving X-raying Dudley's head that resulted in his hair loss. A report by Dr. H.D. Hawks, of his suffering severe hand and chest burns in an X-ray demonstration, was the first of many other reports in Electrical Review.

Other experimenters, including Elihu Thomson and Nikola Tesla, also reported burns. Thomson deliberately exposed a finger to an X-ray tube over a period of time and suffered pain, swelling, and blistering. Other effects, including ultraviolet rays and ozone, were sometimes blamed for the damage, and many physicians still claimed that there were no effects from X-ray exposure at all.

Despite this, there were some early systematic hazard investigations, and as early as 1902 William Herbert Rollins wrote almost despairingly that his warnings about the dangers involved in the careless use of X-rays were not being heeded, either by industry or by his colleagues. By this time, Rollins had proved that X-rays could kill experimental animals, could cause a pregnant guinea pig to abort, and that they could kill a foetus. He also stressed that "animals vary in susceptibility to the external action of X-light" and warned that these differences be considered when patients were treated by means of X-rays.

Radioactive substances

Radioactivity is characteristic of elements with large atomic numbers. Elements with at least one stable isotope are shown in light blue. Green shows elements of which the most stable isotope has a half-life measured in millions of years. Yellow and orange are progressively less stable, with half-lives in thousands or hundreds of years, down toward one day. Red and purple show highly and extremely radioactive elements where the most stable isotopes exhibit half-lives measured on the order of one day and much less.

However, the biological effects of radiation due to radioactive substances were less easy to gauge. This gave the opportunity for many physicians and corporations to market radioactive substances as patent medicines. Examples were radium enema treatments, and radium-containing waters to be drunk as tonics. Marie Curie protested against this sort of treatment, warning that "radium is dangerous in untrained hands". Curie later died from aplastic anaemia, likely caused by exposure to ionizing radiation. By the 1930s, after a number of cases of bone necrosis and death of radium treatment enthusiasts, radium-containing medicinal products had been largely removed from the market (radioactive quackery).

Radiation protection

Only a year after Röntgen's discovery of X rays, the American engineer Wolfram Fuchs (1896) gave what is probably the first protection advice, but it was not until 1925 that the first International Congress of Radiology (ICR) was held and considered establishing international protection standards. The effects of radiation on genes, including the effect of cancer risk, were recognized much later. In 1927, Hermann Joseph Muller published research showing genetic effects and, in 1946, was awarded the Nobel Prize in Physiology or Medicine for his findings.

The second ICR was held in Stockholm in 1928 and proposed the adoption of the röntgen unit, and the International X-ray and Radium Protection Committee (IXRPC) was formed. Rolf Sievert was named chairman, but a driving force was George Kaye of the British National Physical Laboratory. The committee met in 1931, 1934, and 1937.

After World War II, the increased range and quantity of radioactive substances being handled as a result of military and civil nuclear programs led to large groups of occupational workers and the public being potentially exposed to harmful levels of ionising radiation. This was considered at the first post-war ICR convened in London in 1950, when the present International Commission on Radiological Protection (ICRP) was born. Since then the ICRP has developed the present international system of radiation protection, covering all aspects of radiation hazards.

In 2020, Hauptmann and other 15 international researchers of eight nations, among which: Institutes of Biostatistics, Registry Research, Centers of Cancer Epidemiology, Radiation Epidemiology, and then also U.S. National Cancer Institute (NCI), International Agency for Research on Cancer (IARC) and Radiation Effects Research Foundation of Hiroshima studied definitively through meta-analysis the damage resulting from the "low doses" that have afflicted the populations of survivors of the atomic bombings of Hiroshima and Nagasaki and also in numerous accidents of nuclear plants that have occurred in the world. These scientists reported, in JNCI Monographs: Epidemiological Studies of Low Dose Ionizing Radiation and Cancer Risk, that the new epidemiological studies directly support excess cancer risks from low-dose ionizing radiation. In 2021, Italian researcher Venturi reported the first correlations between radio-caesium and pancreatic cancer with the role of caesium in biology and in pancreatitis and in diabetes of pancreatic origin.

Units

Graphic showing relationships between radioactivity and detected ionizing radiation

The International System of Units (SI) unit of radioactive activity is the becquerel (Bq), named in honor of the scientist Henri Becquerel. One Bq is defined as one transformation (or decay or disintegration) per second.

An older unit of radioactivity is the curie, Ci, which was originally defined as "the quantity or mass of radium emanation in equilibrium with one gram of radium (element)". Today, the curie is defined as 3.7×10¹⁰ disintegrations per second, so that 1 curie (Ci) = 3.7×10¹⁰ Bq. For radiological protection purposes, although the United States Nuclear Regulatory Commission permits the use of the unit curie alongside SI units, the European Union European units of measurement directives required that its use for "public health ... purposes" be phased out by 31 December 1985.

The effects of ionizing radiation are often measured in units of gray for mechanical or sievert for damage to tissue.

Types

Radioactive decay results in a reduction of summed rest mass, once the released energy (the disintegration energy) has escaped in some way. Although decay energy is sometimes defined as associated with the difference between the mass of the parent nuclide products and the mass of the decay products, this is true only of rest mass measurements, where some energy has been removed from the product system. This is true because the decay energy must always carry mass with it, wherever it appears (see mass in special relativity) according to the formula E = mc². The decay energy is initially released as the energy of emitted photons plus the kinetic energy of massive emitted particles (that is, particles that have rest mass). If these particles come to thermal equilibrium with their surroundings and photons are absorbed, then the decay energy is transformed to thermal energy, which retains its mass.

Decay energy, therefore, remains associated with a certain measure of the mass of the decay system, called invariant mass, which does not change during the decay, even though the energy of decay is distributed among decay particles. The energy of photons, the kinetic energy of emitted particles, and, later, the thermal energy of the surrounding matter, all contribute to the invariant mass of the system. Thus, while the sum of the rest masses of the particles is not conserved in radioactive decay, the system mass and system invariant mass (and also the system total energy) is conserved throughout any decay process. This is a restatement of the equivalent laws of conservation of energy and conservation of mass.

Alpha, beta and gamma decay

Alpha particles may be completely stopped by a sheet of paper, beta particles by aluminium shielding. Gamma rays can only be reduced by much more substantial mass, such as a very thick layer of lead.

Early researchers found that an electric or magnetic field could split radioactive emissions into three types of beams. The rays were given the names alpha, beta, and gamma, in increasing order of their ability to penetrate matter. Alpha decay is observed only in heavier elements of atomic number 52 (tellurium) and greater, with the exception of beryllium-8 (which decays to two alpha particles). The other two types of decay are observed in all the elements. Lead, atomic number 82, is the heaviest element to have any isotopes stable (to the limit of measurement) to radioactive decay. Radioactive decay is seen in all isotopes of all elements of atomic number 83 (bismuth) or greater. Bismuth-209, however, is only very slightly radioactive, with a half-life greater than the age of the universe; radioisotopes with extremely long half-lives are considered effectively stable for practical purposes.

Transition diagram for decay modes of a radionuclide, with neutron number N and atomic number Z (shown are α, β^±, p⁺, and n⁰ emissions, EC denotes electron capture).

In analysing the nature of the decay products, it was obvious from the direction of the electromagnetic forces applied to the radiations by external magnetic and electric fields that alpha particles carried a positive charge, beta particles carried a negative charge, and gamma rays were neutral. From the magnitude of deflection, it was clear that alpha particles were much more massive than beta particles. Passing alpha particles through a very thin glass window and trapping them in a discharge tube allowed researchers to study the emission spectrum of the captured particles, and ultimately proved that alpha particles are helium nuclei. Other experiments showed beta radiation, resulting from decay and cathode rays, were high-speed electrons. Likewise, gamma radiation and X-rays were found to be high-energy electromagnetic radiation.

The relationship between the types of decays also began to be examined: For example, gamma decay was almost always found to be associated with other types of decay, and occurred at about the same time, or afterwards. Gamma decay as a separate phenomenon, with its own half-life (now termed isomeric transition), was found in natural radioactivity to be a result of the gamma decay of excited metastable nuclear isomers, which were in turn created from other types of decay. Although alpha, beta, and gamma radiations were most commonly found, other types of emission were eventually discovered. Shortly after the discovery of the positron in cosmic ray products, it was realized that the same process that operates in classical beta decay can also produce positrons (positron emission), along with neutrinos (classical beta decay produces antineutrinos).

Electron capture

In electron capture, some proton-rich nuclides were found to capture their own atomic electrons instead of emitting positrons, and subsequently, these nuclides emit only a neutrino and a gamma ray from the excited nucleus (and often also Auger electrons and characteristic X-rays, as a result of the re-ordering of electrons to fill the place of the missing captured electron). These types of decay involve the nuclear capture of electrons or emission of electrons or positrons, and thus acts to move a nucleus toward the ratio of neutrons to protons that has the least energy for a given total number of nucleons. This consequently produces a more stable (lower energy) nucleus.

A hypothetical process of positron capture, analogous to electron capture, is theoretically possible in antimatter atoms, but has not been observed, as complex antimatter atoms beyond antihelium are not experimentally available. Such a decay would require antimatter atoms at least as complex as beryllium-7, which is the lightest known isotope of normal matter to undergo decay by electron capture.

Nucleon emission

Shortly after the discovery of the neutron in 1932, Enrico Fermi realized that certain rare beta-decay reactions immediately yield neutrons as an additional decay particle, so called beta-delayed neutron emission. Neutron emission usually happens from nuclei that are in an excited state, such as the excited ¹⁷O* produced from the beta decay of ¹⁷N. The neutron emission process itself is controlled by the nuclear force and therefore is extremely fast, sometimes referred to as "nearly instantaneous". Isolated proton emission was eventually observed in some elements. It was also found that some heavy elements may undergo spontaneous fission into products that vary in composition. In a phenomenon called cluster decay, specific combinations of neutrons and protons other than alpha particles (helium nuclei) were found to be spontaneously emitted from atoms.

More exotic types of decay

Other types of radioactive decay were found to emit previously seen particles but via different mechanisms. An example is internal conversion, which results in an initial electron emission, and then often further characteristic X-rays and Auger electrons emissions, although the internal conversion process involves neither beta nor gamma decay. A neutrino is not emitted, and none of the electron(s) and photon(s) emitted originate in the nucleus, even though the energy to emit all of them does originate there. Internal conversion decay, like isomeric transition gamma decay and neutron emission, involves the release of energy by an excited nuclide, without the transmutation of one element into another.

Rare events that involve a combination of two beta-decay-type events happening simultaneously are known (see below). Any decay process that does not violate the conservation of energy or momentum laws (and perhaps other particle conservation laws) is permitted to happen, although not all have been detected. An interesting example discussed in a final section, is bound state beta decay of rhenium-187. In this process, the beta electron-decay of the parent nuclide is not accompanied by beta electron emission, because the beta particle has been captured into the K-shell of the emitting atom. An antineutrino is emitted, as in all negative beta decays.

If energy circumstances are favorable, a given radionuclide may undergo many competing types of decay, with some atoms decaying by one route, and others decaying by another. An example is copper-64, which has 29 protons, and 35 neutrons, which decays with a half-life of 12.7004(13) hours. This isotope has one unpaired proton and one unpaired neutron, so either the proton or the neutron can decay to the other particle, which has opposite isospin. This particular nuclide (though not all nuclides in this situation) is more likely to decay through beta plus decay (61.52(26)%) than through electron capture (38.48(26)%). The excited energy states resulting from these decays which fail to end in a ground energy state, also produce later internal conversion and gamma decay in almost 0.5% of the time.

Decay chains and multiple modes

The decay chain of 237-neptunium

The daughter nuclide of a decay event may also be unstable (radioactive). In this case, it too will decay, producing radiation. The resulting second daughter nuclide may also be radioactive. This can lead to a sequence of several decay events called a decay chain (see this article for specific details of important natural decay chains). Eventually, a stable nuclide is produced. Any decay daughters that are the result of an alpha decay will also result in helium atoms being created.

Some radionuclides may have several different paths of decay. For example, 35.94(6)% of bismuth-212 decays, through alpha-emission, to thallium-208 while 64.06(6)% of bismuth-212 decays, through beta-emission, to polonium-212. Both thallium-208 and polonium-212 are radioactive daughter products of bismuth-212, and both decay directly to stable lead-208.

Occurrence and applications

According to the Big Bang theory, stable isotopes of the lightest three elements (H, He, and traces of Li) were produced very shortly after the emergence of the universe, in a process called Big Bang nucleosynthesis. These lightest stable nuclides (including deuterium) survive to today, but any radioactive isotopes of the light elements produced in the Big Bang (such as tritium) have long since decayed. Isotopes of elements heavier than boron were not produced at all in the Big Bang, and these first five elements do not have any long-lived radioisotopes. Thus, all radioactive nuclei are, therefore, relatively young with respect to the birth of the universe, having formed later in various other types of nucleosynthesis in stars (in particular, supernovae), and also during ongoing interactions between stable isotopes and energetic particles. For example, carbon-14, a radioactive nuclide with a half-life of only 5700(30) years, is constantly produced in Earth's upper atmosphere due to interactions between cosmic rays and nitrogen.

Nuclides that are produced by radioactive decay are called radiogenic nuclides, whether they themselves are stable or not. There exist stable radiogenic nuclides that were formed from short-lived extinct radionuclides in the early Solar System. The extra presence of these stable radiogenic nuclides (such as xenon-129 from extinct iodine-129) against the background of primordial stable nuclides can be inferred by various means.

Radioactive decay has been put to use in the technique of radioisotopic labeling, which is used to track the passage of a chemical substance through a complex system (such as a living organism). A sample of the substance is synthesized with a high concentration of unstable atoms. The presence of the substance in one or another part of the system is determined by detecting the locations of decay events.

On the premise that radioactive decay is truly random (rather than merely chaotic), it has been used in hardware random-number generators. Because the process is not thought to vary significantly in mechanism over time, it is also a valuable tool in estimating the absolute ages of certain materials. For geological materials, the radioisotopes and some of their decay products become trapped when a rock solidifies, and can then later be used (subject to many well-known qualifications) to estimate the date of the solidification. These include checking the results of several simultaneous processes and their products against each other, within the same sample. In a similar fashion, and also subject to qualification, the rate of formation of carbon-14 in various eras, the date of formation of organic matter within a certain period related to the isotope's half-life may be estimated, because the carbon-14 becomes trapped when the organic matter grows and incorporates the new carbon-14 from the air. Thereafter, the amount of carbon-14 in organic matter decreases according to decay processes that may also be independently cross-checked by other means (such as checking the carbon-14 in individual tree rings, for example).

Szilard–Chalmers effect

The Szilard–Chalmers effect is the breaking of a chemical bond as a result of a kinetic energy imparted from radioactive decay. It operates by the absorption of neutrons by an atom and subsequent emission of gamma rays, often with significant amounts of kinetic energy. This kinetic energy, by Newton's third law, pushes back on the decaying atom, which causes it to move with enough speed to break a chemical bond. This effect can be used to separate isotopes by chemical means.

The Szilard–Chalmers effect was discovered in 1934 by Leó Szilárd and Thomas A. Chalmers. They observed that after bombardment by neutrons, the breaking of a bond in liquid ethyl iodide allowed radioactive iodine to be removed.

Origins of radioactive nuclides

Radioactive primordial nuclides found in the Earth are residues from ancient supernova explosions that occurred before the formation of the Solar System. They are the fraction of radionuclides that survived from that time, through the formation of the primordial solar nebula, through planet accretion, and up to the present time. The naturally occurring short-lived radiogenic radionuclides found in today's rocks, are the daughters of those radioactive primordial nuclides. Another minor source of naturally occurring radioactive nuclides are cosmogenic nuclides, that are formed by cosmic ray bombardment of material in the Earth's atmosphere or crust. The decay of the radionuclides in rocks of the Earth's mantle and crust contribute significantly to Earth's internal heat budget.

Rates

The decay rate, or activity, of a radioactive substance is characterized by the following time-independent parameters:

• The half-life, t_1/2, is the time taken for the activity of a given amount of a radioactive substance to decay to half of its initial value.
• The decay constant, λ "lambda", the reciprocal of the mean lifetime (in s⁻¹), sometimes referred to as simply decay rate.
• The mean lifetime, τ "tau", the average lifetime (1/e life) of a radioactive particle before decay.

Although these are constants, they are associated with the statistical behavior of populations of atoms. In consequence, predictions using these constants are less accurate for minuscule samples of atoms.

In principle a half-life, a third-life, or even a (1/√2)-life, can be used in exactly the same way as half-life; but the mean life and half-life t_1/2 have been adopted as standard times associated with exponential decay.

Those parameters can be related to the following time-dependent parameters:

• Total activity (or just activity), A, is the number of decays per unit time of a radioactive sample.
• Number of particles, N, in the sample.
• Specific activity, a, is the number of decays per unit time per amount of substance of the sample at time set to zero (t = 0). "Amount of substance" can be the mass, volume or moles of the initial sample.

These are related as follows:

[math]\displaystyle{ \begin{align} t_{1/2} &= \frac{\ln(2)}{\lambda} = \tau \ln(2) \\[2pt] A &= - \frac{\mathrm{d}N}{\mathrm{d}t} = \lambda N = \frac{\ln(2)}{t_{1/2}}N \\[2pt] S_A a_0 &= - \frac{\mathrm{d}N}{\mathrm{d}t}\bigg|_{t=0} = \lambda N_0 \end{align} }[/math]

where N₀ is the initial amount of active substance — substance that has the same percentage of unstable particles as when the substance was formed.

Mathematics

For the mathematical details of exponential decay in general context, see exponential decay.

For related derivations with some further details, see half-life.

For the analogous mathematics in 1st order chemical reactions, see consecutive reactions.

Universal law

The mathematics of radioactive decay depend on a key assumption that a nucleus of a radionuclide has no "memory" or way of translating its history into its present behavior. A nucleus does not "age" with the passage of time. Thus, the probability of its breaking down does not increase with time but stays constant, no matter how long the nucleus has existed. This constant probability may differ greatly between one type of nucleus and another, leading to the many different observed decay rates. However, whatever the probability is, it does not change over time. This is in marked contrast to complex objects that do show aging, such as automobiles and humans. These aging systems do have a chance of breakdown per unit of time that increases from the moment they begin their existence.

Aggregate processes, like the radioactive decay of a lump of atoms, for which the single-event probability of realization is very small but in which the number of time-slices is so large that there is nevertheless a reasonable rate of events, are modelled by the Poisson distribution, which is discrete. Radioactive decay and nuclear particle reactions are two examples of such aggregate processes. The mathematics of Poisson processes reduce to the law of exponential decay, which describes the statistical behaviour of a large number of nuclei, rather than one individual nucleus. In the following formalism, the number of nuclei or the nuclei population N, is of course a discrete variable (a natural number)—but for any physical sample N is so large that it can be treated as a continuous variable. Differential calculus is used to model the behaviour of nuclear decay.

One-decay process

Consider the case of a nuclide A that decays into another B by some process A → B (emission of other particles, like electron neutrinos ν_e and electrons e⁻ as in beta decay, are irrelevant in what follows). The decay of an unstable nucleus is entirely random in time so it is impossible to predict when a particular atom will decay. However, it is equally likely to decay at any instant in time. Therefore, given a sample of a particular radioisotope, the number of decay events −dN expected to occur in a small interval of time dt is proportional to the number of atoms present N, that is

[math]\displaystyle{ - \frac{\mathrm{d}N}{\mathrm{d}t} \propto N }[/math]

Particular radionuclides decay at different rates, so each has its own decay constant λ. The expected decay −dN/N is proportional to an increment of time, dt:

[math]\displaystyle{ -\frac{\mathrm{d}N}{N} = \lambda \mathrm{d}t }[/math]

The negative sign indicates that N decreases as time increases, as the decay events follow one after another. The solution to this first-order differential equation is the function:

[math]\displaystyle{ N(t) = N_0\,e^{-{\lambda}t} }[/math]

where N₀ is the value of N at time t = 0, with the decay constant expressed as λ

We have for all time t:

[math]\displaystyle{ N_A + N_B = N_\text{total} = N_{A0}, }[/math]

where N_total is the constant number of particles throughout the decay process, which is equal to the initial number of A nuclides since this is the initial substance.

If the number of non-decayed A nuclei is:

[math]\displaystyle{ N_A = N_{A0}e^{-\lambda t} }[/math]

then the number of nuclei of B (i.e. the number of decayed A nuclei) is

[math]\displaystyle{ N_B = N_{A0} - N_A = N_{A0} - N_{A0}e^{-\lambda t} = N_{A0} \left( 1 - e^{-\lambda t} \right). }[/math]

The number of decays observed over a given interval obeys Poisson statistics. If the average number of decays is ⟨N⟩, the probability of a given number of decays N is

[math]\displaystyle{ P(N) = \frac{\langle N \rangle^N \exp(-\langle N\rangle)}{N!} . }[/math]

Chain-decay processes

Chain of two decays

Now consider the case of a chain of two decays: one nuclide A decaying into another B by one process, then B decaying into another C by a second process, i.e. A → B → C. The previous equation cannot be applied to the decay chain, but can be generalized as follows. Since A decays into B, then B decays into C, the activity of A adds to the total number of B nuclides in the present sample, before those B nuclides decay and reduce the number of nuclides leading to the later sample. In other words, the number of second generation nuclei B increases as a result of the first generation nuclei decay of A, and decreases as a result of its own decay into the third generation nuclei C. The sum of these two terms gives the law for a decay chain for two nuclides:

[math]\displaystyle{ \frac{\mathrm{d}N_B}{\mathrm{d}t} = -\lambda_B N_B + \lambda_A N_A. }[/math]

The rate of change of N_B, that is dN_B/dt, is related to the changes in the amounts of A and B, N_B can increase as B is produced from A and decrease as B produces C.

Re-writing using the previous results:

[math]\displaystyle{ \frac{\mathrm{d}N_B}{\mathrm{d}t} = - \lambda_B N_B + \lambda_A N_{A0} e^{-\lambda_A t} }[/math]

The subscripts simply refer to the respective nuclides, i.e. N_A is the number of nuclides of type A; N_A0 is the initial number of nuclides of type A; λ_A is the decay constant for A – and similarly for nuclide B. Solving this equation for N_B gives:

[math]\displaystyle{ N_B = \frac{N_{A0}\lambda_A}{\lambda_B - \lambda_A} \left( e^{-\lambda_A t} - e^{-\lambda_B t}\right) . }[/math]

In the case where B is a stable nuclide (λ_B = 0), this equation reduces to the previous solution:

[math]\displaystyle{ \lim_{\lambda_B\rightarrow 0} \left[ \frac{N_{A0}\lambda_A}{\lambda_B - \lambda_A} \left( e^{-\lambda_A t} - e^{-\lambda_B t} \right) \right] = \frac{N_{A0}\lambda_A}{0 - \lambda_A} \left( e^{-\lambda_A t} - 1 \right) = N_{A0} \left( 1- e^{-\lambda_A t} \right), }[/math]

as shown above for one decay. The solution can be found by the integration factor method, where the integrating factor is e^λ_Bt. This case is perhaps the most useful since it can derive both the one-decay equation (above) and the equation for multi-decay chains (below) more directly.

Chain of any number of decays

For the general case of any number of consecutive decays in a decay chain, i.e. A₁ → A₂ ··· → A_i ··· → A_D, where D is the number of decays and i is a dummy index (i = _{1, 2, 3, ...D}), each nuclide population can be found in terms of the previous population. In this case N₂ = 0, N₃ = 0, ..., N_D = 0. Using the above result in a recursive form:

[math]\displaystyle{ \frac{\mathrm{d}N_j}{\mathrm{d}t} = - \lambda_j N_j + \lambda_{j-1} N_{(j-1)0} e^{-\lambda_{j-1} t}. }[/math]

The general solution to the recursive problem is given by Bateman's equations:

'Bateman's equations'

[math]\displaystyle{ \begin{align} N_D &= \frac{N_1(0)}{\lambda_D} \sum_{i=1}^D \lambda_i c_i e^{-\lambda_i t} \\[3pt] c_i &= \prod_{j=1, i\neq j}^D \frac{\lambda_j}{\lambda_j - \lambda_i} \end{align} }[/math]

Alternative modes

In all of the above examples, the initial nuclide decays into just one product. Consider the case of one initial nuclide that can decay into either of two products, that is A → B and A → C in parallel. For example, in a sample of potassium-40, 89.3% of the nuclei decay to calcium-40 and 10.7% to argon-40. We have for all time t:

[math]\displaystyle{ N = N_A + N_B + N_C }[/math]

which is constant, since the total number of nuclides remains constant. Differentiating with respect to time:

[math]\displaystyle{ \begin{align} \frac{\mathrm{d}N_A}{\mathrm{d}t} & = - \left(\frac{\mathrm{d}N_B}{\mathrm{d}t} + \frac{\mathrm{d}N_C}{\mathrm{d}t} \right) \\ - \lambda N_A & = - N_A \left ( \lambda_B + \lambda_C \right ) \\ \end{align} }[/math]

defining the total decay constant λ in terms of the sum of partial decay constants λ_B and λ_C:

[math]\displaystyle{ \lambda = \lambda_B + \lambda_C . }[/math]

Solving this equation for N_A:

[math]\displaystyle{ N_A = N_{A0} e^{-\lambda t} . }[/math]

where N_A0 is the initial number of nuclide A. When measuring the production of one nuclide, one can only observe the total decay constant λ. The decay constants λ_B and λ_C determine the probability for the decay to result in products B or C as follows:

[math]\displaystyle{ N_B = \frac{\lambda_B}{\lambda} N_{A0} \left ( 1 - e^{-\lambda t} \right ), }[/math]

[math]\displaystyle{ N_C = \frac{\lambda_C}{\lambda} N_{A0} \left ( 1 - e^{-\lambda t} \right ). }[/math]

because the fraction λ_B/λ of nuclei decay into B while the fraction λ_C/λ of nuclei decay into C.

Corollaries of laws

The above equations can also be written using quantities related to the number of nuclide particles N in a sample;

• The activity: A = λN.
• The amount of substance: n = N/N_A.
• The mass: m = Mn = MN/N_A.

where N_A = 6.02214076×10²³ mol⁻¹ is the Avogadro constant, M is the molar mass of the substance in kg/mol, and the amount of the substance n is in moles.

Decay timing: definitions and relations

Time constant and mean-life

For the one-decay solution A → B:

[math]\displaystyle{ N = N_0\,e^{-{\lambda}t} = N_0\,e^{-t/ \tau}, \,\! }[/math]

the equation indicates that the decay constant λ has units of t⁻¹, and can thus also be represented as 1/τ, where τ is a characteristic time of the process called the time constant.

In a radioactive decay process, this time constant is also the mean lifetime for decaying atoms. Each atom "lives" for a finite amount of time before it decays, and it may be shown that this mean lifetime is the arithmetic mean of all the atoms' lifetimes, and that it is τ, which again is related to the decay constant as follows:

[math]\displaystyle{ \tau = \frac{1}{\lambda}. }[/math]

This form is also true for two-decay processes simultaneously A → B + C, inserting the equivalent values of decay constants (as given above)

[math]\displaystyle{ \lambda = \lambda_B + \lambda_C \, }[/math]

into the decay solution leads to:

[math]\displaystyle{ \frac{1}{\tau} = \lambda = \lambda_B + \lambda_C = \frac{1}{\tau_B} + \frac{1}{\tau_C}\, }[/math]

Simulation of many identical atoms undergoing radioactive decay, starting with either 4 atoms (left) or 400 (right). The number at the top indicates how many half-lives have elapsed.

Half-life

A more commonly used parameter is the half-life T_1/2. Given a sample of a particular radionuclide, the half-life is the time taken for half the radionuclide's atoms to decay. For the case of one-decay nuclear reactions:

[math]\displaystyle{ N = N_0\,e^{-{\lambda}t} = N_0\,e^{-t/\tau}, \,\! }[/math]

the half-life is related to the decay constant as follows: set N = N₀/2 and t = T_1/2 to obtain

[math]\displaystyle{ t_{1/2} = \frac{\ln 2}{\lambda} = \tau \ln 2. }[/math]

This relationship between the half-life and the decay constant shows that highly radioactive substances are quickly spent, while those that radiate weakly endure longer. Half-lives of known radionuclides vary by almost 54 orders of magnitude, from more than 2.25(9)×10²⁴ years (6.9×10³¹ sec) for the very nearly stable nuclide ¹²⁸Te, to 8.6(6)×10⁻²³ seconds for the highly unstable nuclide ⁵H.

The factor of ln(2) in the above relations results from the fact that the concept of "half-life" is merely a way of selecting a different base other than the natural base e for the lifetime expression. The time constant τ is the e −1 -life, the time until only 1/e remains, about 36.8%, rather than the 50% in the half-life of a radionuclide. Thus, τ is longer than t_1/2. The following equation can be shown to be valid:

[math]\displaystyle{ N(t) = N_0\,e^{-t/\tau} = N_0\,2^{-t/t_{1/2}}. \,\! }[/math]

Since radioactive decay is exponential with a constant probability, each process could as easily be described with a different constant time period that (for example) gave its "(1/3)-life" (how long until only 1/3 is left) or "(1/10)-life" (a time period until only 10% is left), and so on. Thus, the choice of τ and t_1/2 for marker-times, are only for convenience, and from convention. They reflect a fundamental principle only in so much as they show that the same proportion of a given radioactive substance will decay, during any time-period that one chooses.

Mathematically, the n^th life for the above situation would be found in the same way as above—by setting N = N₀/n, t = T_1/n and substituting into the decay solution to obtain

[math]\displaystyle{ t_{1/n} = \frac{\ln n}{\lambda} = \tau \ln n. }[/math]

Example for carbon-14

Carbon-14 has a half-life of 5700(30) years and a decay rate of 14 disintegrations per minute (dpm) per gram of natural carbon.

If an artifact is found to have radioactivity of 4 dpm per gram of its present C, we can find the approximate age of the object using the above equation:

[math]\displaystyle{ N = N_0\,e^{-t/ \tau}, }[/math]

where:

[math]\displaystyle{ \begin{align} \frac{N}{N_0} &= 4/14 \approx 0.286, \\ \tau &= \frac{T_{1/2}}{\ln 2} \approx 8267\text{ years}, \\ t &= -\tau\,\ln\frac{N}{ N_0} \approx 10356\text{ years}. \end{align} }[/math]

Changing rates

The radioactive decay modes of electron capture and internal conversion are known to be slightly sensitive to chemical and environmental effects that change the electronic structure of the atom, which in turn affects the presence of 1s and 2s electrons that participate in the decay process. A small number of nuclides are affected. For example, chemical bonds can affect the rate of electron capture to a small degree (in general, less than 1%) depending on the proximity of electrons to the nucleus. In ⁷Be, a difference of 0.9% has been observed between half-lives in metallic and insulating environments. This relatively large effect is because beryllium is a small atom whose valence electrons are in 2s atomic orbitals, which are subject to electron capture in ⁷Be because (like all s atomic orbitals in all atoms) they naturally penetrate into the nucleus.

In 1992, Jung et al. of the Darmstadt Heavy-Ion Research group observed an accelerated β⁻ decay of ¹⁶³Dy⁶⁶⁺. Although neutral ¹⁶³Dy is a stable isotope, the fully ionized ¹⁶³Dy⁶⁶⁺ undergoes β⁻ decay into the K and L shells to ¹⁶³Ho⁶⁶⁺ with a half-life of 47 days.

Rhenium-187 is another spectacular example. ¹⁸⁷Re normally undergoes beta decay to ¹⁸⁷Os with a half-life of 41.6 × 10⁹ years, but studies using fully ionised ¹⁸⁷Re atoms (bare nuclei) have found that this can decrease to only 32.9 years. This is attributed to "bound-state β⁻ decay" of the fully ionised atom – the electron is emitted into the "K-shell" (1s atomic orbital), which cannot occur for neutral atoms in which all low-lying bound states are occupied.

Example of diurnal and seasonal variations in gamma ray detector response.

A number of experiments have found that decay rates of other modes of artificial and naturally occurring radioisotopes are, to a high degree of precision, unaffected by external conditions such as temperature, pressure, the chemical environment, and electric, magnetic, or gravitational fields. Comparison of laboratory experiments over the last century, studies of the Oklo natural nuclear reactor (which exemplified the effects of thermal neutrons on nuclear decay), and astrophysical observations of the luminosity decays of distant supernovae (which occurred far away so the light has taken a great deal of time to reach us), for example, strongly indicate that unperturbed decay rates have been constant (at least to within the limitations of small experimental errors) as a function of time as well.

Recent results suggest the possibility that decay rates might have a weak dependence on environmental factors. It has been suggested that measurements of decay rates of silicon-32, manganese-54, and radium-226 exhibit small seasonal variations (of the order of 0.1%). However, such measurements are highly susceptible to systematic errors, and a subsequent paper has found no evidence for such correlations in seven other isotopes (²²Na, ⁴⁴Ti, ¹⁰⁸Ag, ¹²¹Sn, ¹³³Ba, ²⁴¹Am, ²³⁸Pu), and sets upper limits on the size of any such effects. The decay of radon-222 was once reported to exhibit large 4% peak-to-peak seasonal variations (see plot), which were proposed to be related to either solar flare activity or the distance from the Sun, but detailed analysis of the experiment's design flaws, along with comparisons to other, much more stringent and systematically controlled, experiments refute this claim.

GSI anomaly

An unexpected series of experimental results for the rate of decay of heavy highly charged radioactive ions circulating in a storage ring has provoked theoretical activity in an effort to find a convincing explanation. The rates of weak decay of two radioactive species with half lives of about 40 s and 200 s are found to have a significant oscillatory modulation, with a period of about 7 s. The observed phenomenon is known as the GSI anomaly, as the storage ring is a facility at the GSI Helmholtz Centre for Heavy Ion Research in Darmstadt, Germany . As the decay process produces an electron neutrino, some of the proposed explanations for the observed rate oscillation invoke neutrino properties. Initial ideas related to flavour oscillation met with skepticism. A more recent proposal involves mass differences between neutrino mass eigenstates.

Theoretical basis

The neutrons and protons that constitute nuclei, as well as other particles that approach close enough to them, are governed by several interactions. The nuclear force (also known as residual strong force), not observed at the familiar macroscopic scale, is the most powerful force over subatomic distances. The electrostatic force is almost always significant, and, in the case of beta decay, the weak nuclear force is also involved.

The combined effects of these forces produces a number of different phenomena in which energy may be released by rearrangement of particles in the nucleus, or else the change of one type of particle into others. These rearrangements and transformations may be hindered energetically so that they do not occur immediately. In certain cases, random quantum vacuum fluctuations are theorized to promote relaxation to a lower energy state (the "decay") in a phenomenon known as quantum tunneling. Radioactive decay half-life of nuclides has been measured over timescales of 54 orders of magnitude, from 8.6(6)×10⁻²³ seconds (for hydrogen-5) to 7.10(28)×10³¹ seconds (for tellurium-128). The limits of these timescales are set by the sensitivity of instrumentation only, and there are no known natural limits to how brief or long a decay half-life for radioactive decay of a radionuclide may be.

The decay process, like all hindered energy transformations, may be analogized by a snowfield on a mountain. While friction between the ice crystals may be supporting the snow's weight, the system is inherently unstable with regard to a state of lower potential energy. A disturbance would thus facilitate the path to a state of greater entropy; the system will move towards the ground state, producing heat, and the total energy will be distributable over a larger number of quantum states thus resulting in an avalanche. The total energy does not change in this process, but, because of the second law of thermodynamics, avalanches have only been observed in one direction and that is toward the "ground state" — the state with the largest number of ways in which the available energy could be distributed.

Such a collapse (a gamma-ray decay event) requires a specific activation energy. For a snow avalanche, this energy comes as a disturbance from outside the system, although such disturbances can be arbitrarily small. In the case of an excited atomic nucleus decaying by gamma radiation in a spontaneous emission of electromagnetic radiation, the arbitrarily small disturbance is theorized to come from quantum vacuum fluctuations.

A radioactive nucleus (or any excited system in quantum mechanics) is unstable, and can, thus, spontaneously stabilize to a less-excited system. The resulting transformation alters the structure of the nucleus and results in the emission of either a photon or a high-velocity particle that has mass (such as an electron, alpha particle, or other type).

Hazard warning signs

• 

  The trefoil symbol used to warn of presence of radioactive material or ionising radiation

• 

  2007 ISO radioactivity hazard symbol intended for IAEA Category 1, 2 and 3 sources defined as dangerous sources capable of death or serious injury

• 

  One of several dangerous goods transport classification signs for radioactive materials

See also

• Actinides in the environment
• Background radiation
• Chernobyl disaster
• Crimes involving radioactive substances
• Decay chain
• Decay correction
• Fallout shelter
• Geiger counter
• Induced radioactivity
• Lists of nuclear disasters and radioactive incidents
• National Council on Radiation Protection and Measurements
• Nuclear engineering
• Nuclear pharmacy
• Nuclear physics
• Nuclear power
• Nuclear chain reaction
• Particle decay
• Poisson process
• Radiation therapy
• Radioactive contamination
• Radioactivity in biology
• Radiometric dating
• Transient equilibrium

References

General

• "Radioactivity", Encyclopædia Britannica. 2006. Encyclopædia Britannica Online. December 18, 2006
• Radio-activity by Ernest Rutherford Phd, Encyclopædia Britannica Eleventh Edition

External links

• The Lund/LBNL Nuclear Data Search – Contains tabulated information on radioactive decay types and energies.
• Nomenclature of nuclear chemistry
• Specific activity and related topics.
• The Live Chart of Nuclides – IAEA
• Interactive Chart of Nuclides
• Health Physics Society Public Education Website
•  Beach, Chandler B., ed (1914). "Becquerel Rays". The New Student's Reference Work. Chicago: F. E. Compton and Co. 
• Annotated bibliography for radioactivity from the Alsos Digital Library for Nuclear Issues
• Stochastic Java applet on the decay of radioactive atoms by Wolfgang Bauer
• Stochastic Flash simulation on the decay of radioactive atoms by David M. Harrison
• "Henri Becquerel: The Discovery of Radioactivity", Becquerel's 1896 articles online and analyzed on BibNum [click 'à télécharger' for English version].
• "Radioactive change", Rutherford & Soddy article (1903), online and analyzed on Bibnum [click 'à télécharger' for English version].