Decay Chain Calculator
Free Decay chain Calculator for nuclear physics. Enter variables to compute results with formulas and detailed steps.
Calculator
Adjust values & calculateFormula
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.
Last reviewed: December 2025
Worked Examples
Example 1: Radium-226 to Radon-222 Chain
Example 2: Mo-99 to Tc-99m Medical Generator
Background & Theory
The Decay Chain Calculator applies the following established principles and formulas. Physics is the fundamental natural science concerned with matter, energy, and the interactions between them. Classical mechanics, founded on Newton's three laws of motion, provides the framework for analyzing the motion of objects. The first law states that an object remains at rest or in uniform motion unless acted upon by a net external force. The second law quantifies this relationship: F = ma, where force equals mass times acceleration in SI units of newtons (N = kgยทm/sยฒ). The third law establishes that every action produces an equal and opposite reaction. Kinematics describes motion without reference to its causes. The four fundamental equations relate displacement s, initial velocity u, final velocity v, acceleration a, and time t: v = u + at, s = ut + ยฝatยฒ, vยฒ = uยฒ + 2as, and s = ยฝ(u + v)t. These assume constant acceleration and are foundational for solving projectile motion, free fall, and linear dynamics problems. Energy conservation underpins much of physics. Kinetic energy is KE = ยฝmvยฒ, where m is mass in kilograms and v is speed in meters per second. Gravitational potential energy is PE = mgh, where g โ 9.81 m/sยฒ near Earth's surface and h is height in meters. The work-energy theorem states that the net work done on an object equals its change in kinetic energy: W = ฮKE. Electricity and circuits rely on Ohm's law: V = IR, where voltage V is in volts, current I in amperes, and resistance R in ohms. Electrical power is P = IV = IยฒR = Vยฒ/R, measured in watts. Wave mechanics connects frequency f, wave speed v, and wavelength ฮป through f = v/ฮป, with frequency in hertz (Hz). Pressure is defined as force per unit area, P = F/A, in pascals (Pa = N/mยฒ). The ideal gas law PV = nRT links pressure, volume, moles n, the gas constant R = 8.314 J/(molยทK), and absolute temperature in kelvin. Gravitational force between two masses follows Newton's law of universal gravitation: F = Gmโmโ/rยฒ, where G = 6.674ร10โปยนยน Nยทmยฒ/kgยฒ is the gravitational constant.
History
The history behind the Decay Chain Calculator traces back through the following developments. The history of physics spans over two millennia, beginning with the natural philosophy of ancient Greece. Aristotle (384โ322 BCE) proposed that all matter consisted of four elements and that objects moved toward their natural place, with heavier objects falling faster than lighter ones. While largely incorrect, his systematic approach to explaining nature dominated Western thought for nearly 2,000 years. The Scientific Revolution overturned Aristotelian physics. Galileo Galilei (1564โ1642) performed groundbreaking experiments on inclined planes and falling bodies, demonstrating that all objects fall with the same acceleration regardless of mass, and established the principle of inertia. His use of mathematics to describe motion was revolutionary. Isaac Newton synthesized these developments in his landmark Principia Mathematica (1687), laying out the three laws of motion and the law of universal gravitation. Newton's framework unified terrestrial and celestial mechanics, explaining planetary orbits with the same equations governing a falling apple. His calculus provided the mathematical language for expressing rates of change. The 19th century brought two major theoretical achievements. James Clerk Maxwell formulated his equations of electromagnetism between 1861 and 1862, unifying electricity, magnetism, and optics, and predicting the existence of electromagnetic waves traveling at the speed of light. Thermodynamics was developed by Carnot, Clausius, and Kelvin, establishing the laws governing heat, work, and entropy. The 20th century produced two revolutions that fundamentally altered the classical picture. Albert Einstein published the special theory of relativity in 1905, showing that space and time are not absolute but relative to the observer, and that mass and energy are equivalent via E = mcยฒ. His general theory of relativity in 1915 reinterpreted gravity as the curvature of spacetime. Simultaneously, quantum mechanics emerged from the work of Planck, Bohr, Heisenberg, and Schrรถdinger, revealing that at atomic scales energy is quantized and particles exhibit wave-particle duality. These developments culminated in the Standard Model of particle physics, which describes all known fundamental particles and three of the four fundamental forces.
Frequently Asked Questions
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.
References
Reviewed by Manoj Kumar, Mathematics Educator ยท Editorial policy