Decay Chain Calculator
Free Decay chain Calculator for nuclear physics. Enter variables to compute results with formulas and detailed steps.
Formula
N2(t) = N0 * (lambda1/(lambda2-lambda1)) * (exp(-lambda1*t) - exp(-lambda2*t))
Where N2 is the daughter population at time t, N0 is the initial parent atoms, lambda1 and lambda2 are the parent and daughter decay constants (ln2/half-life), and t is the elapsed time. This is the Bateman equation for a two-member chain.
Worked Examples
Example 1: Radium-226 to Radon-222 Chain
Problem: Starting with 1e20 atoms of Ra-226 (t1/2 = 1600 years), calculate the Rn-222 (t1/2 = 3.82 days) population after 30 days.
Solution: lambda1 = ln2 / 1600 yr = 4.332e-4 yr^-1\nlambda2 = ln2 / (3.82/365.25) yr = 66.25 yr^-1\nt = 30/365.25 = 0.0821 yr\nN_Ra = 1e20 * exp(-4.332e-4 * 0.0821) = 9.9996e19 (essentially unchanged)\nN_Rn = 1e20 * (4.332e-4 / (66.25 - 4.332e-4)) * (exp(-4.332e-4*0.0821) - exp(-66.25*0.0821))\n= 1e20 * 6.539e-6 * (0.9999 - 0.00432) = 6.51e14 atoms
Result: Ra-226: ~1.00e20 atoms | Rn-222: ~6.51e14 atoms | Secular equilibrium nearly reached
Example 2: Mo-99 to Tc-99m Medical Generator
Problem: A Mo-99 (t1/2 = 66 hours) generator starts with 1e18 atoms. Find Tc-99m (t1/2 = 6.01 hours) population after 24 hours.
Solution: lambda1 = ln2 / 66h = 0.01050 h^-1\nlambda2 = ln2 / 6.01h = 0.11534 h^-1\nt = 24 h\nN_Mo = 1e18 * exp(-0.01050 * 24) = 7.77e17\nN_Tc = 1e18 * (0.01050/(0.11534-0.01050)) * (exp(-0.01050*24) - exp(-0.11534*24))\n= 1e18 * 0.1001 * (0.7772 - 0.0633) = 7.15e16 atoms
Result: Mo-99: 7.77e17 atoms | Tc-99m: 7.15e16 atoms | Transient equilibrium
Frequently Asked Questions
What is the Bateman equation and how does it describe decay chain kinetics?
The Bateman equations are a set of coupled first-order differential equations that describe the number of atoms of each member in a radioactive decay chain as a function of time. For a simple two-member chain, the daughter population equals N0 times lambda1 divided by (lambda2 minus lambda1) times the difference of two exponentials. For longer chains, the general Bateman solution involves sums of exponential terms with coefficients determined by all the decay constants in the chain. These equations assume that each decay produces exactly one daughter atom and that branching ratios are unity. The Bateman equations are the mathematical foundation for Decay Chain Calculator and enable prediction of the complete time evolution of every isotope in the chain.
What is secular equilibrium in a decay chain and when does it occur?
Secular equilibrium occurs when the half-life of the parent isotope is much longer (at least 100 times) than the half-life of any daughter in the chain. Under these conditions, the parent activity remains essentially constant over many daughter half-lives, and after a transient buildup period, each daughter reaches a steady-state population where its activity equals the parent activity. In secular equilibrium, the number of daughter atoms equals N_parent times lambda_parent divided by lambda_daughter. A classic example is radium-226 (half-life 1,600 years) in secular equilibrium with radon-222 (half-life 3.82 days). This principle is widely used in radiometric dating and environmental radiation monitoring.
How do you calculate the time of maximum daughter activity in a decay chain?
The time at which the daughter population reaches its maximum value can be calculated analytically from the Bateman equations. Taking the derivative of the daughter population equation with respect to time and setting it equal to zero gives t_max equals ln(lambda2/lambda1) divided by (lambda2 minus lambda1). At this time, the daughter decay rate exactly equals the daughter production rate from parent decay. Before t_max, the daughter is building up faster than it decays; after t_max, it decays faster than it is produced. This maximum time depends only on the two decay constants, not on the initial parent population. Knowledge of t_max is important for medical isotope production and nuclear waste management planning.
How does branching decay affect the calculations in a decay chain?
Branching decay occurs when a radioactive isotope can decay through more than one mode, such as both alpha and beta decay, producing different daughter products. When branching occurs, the decay constant for each branch equals the total decay constant multiplied by the branching ratio (the fraction decaying through that branch). This means the total half-life remains unchanged, but the effective production rate of each daughter is reduced by the branching fraction. In complex decay chains like the uranium series, several isotopes exhibit branching, creating parallel paths that eventually converge. Decay Chain Calculator assumes a single decay path for simplicity, but the input decay constants can be adjusted to represent the effective constants for a specific branch.
What practical applications rely on decay chain calculations?
Decay chain calculations have numerous practical applications across science and industry. In nuclear medicine, production of technetium-99m from molybdenum-99 generators requires precise knowledge of the Mo-99 to Tc-99m decay chain dynamics to optimize elution schedules. In geological dating, the uranium-lead method exploits the complete U-238 decay chain to date rocks billions of years old. Nuclear power plants must track the buildup of fission product decay chains to predict decay heat after reactor shutdown. Environmental monitoring uses radon-222 measurements (a member of the U-238 chain) to assess radioactive contamination. Nuclear forensics analyzes isotope ratios along decay chains to determine the age and origin of nuclear materials.
How does the initial condition affect decay chain evolution?
The initial conditions profoundly affect how a decay chain evolves over time. If only the parent isotope is present initially (a freshly purified sample), all daughters must build up from zero, and the system takes several daughter half-lives to approach equilibrium. If the system starts in secular equilibrium (as in an undisturbed ore body), all activities are equal and remain so as long as the parent activity does not change significantly. If daughters are selectively removed (as in a molybdenum-technetium generator), the daughter repopulation follows the Bateman equations with modified initial conditions. Understanding these different starting scenarios is crucial for correctly interpreting measurements and planning experiments in nuclear science.