Decay Chain Parent Daughter Interactive Calculator

Parent Decay

A_P(t) = A_P,0 e^-λ_Pt

Daughter Activity (Bateman Equation)

A_D(t) = A_D,0 e^-λ_Dt + (λ_D/(λ_D - λ_P)) A_P,0 (e^-λ_Pt - e^-λ_Dt)

Activity Ratio at Equilibrium

A_D/A_P = λ_D/(λ_D - λ_P)

(Valid for transient and secular equilibrium)

Decay Constants

λ = ln(2)/t_1/2 = 0.693147/t_1/2

Effective Half-Life

t_1/2,eff = ln(2)/(λ_P + λ_D)

(For combined parent-daughter decay)

Scenario: Hospital Radiopharmacy Daily Operations

Elena, a nuclear medicine radiopharmacist at a 400-bed teaching hospital, manages daily ⁹⁹Mo/^99mTc generator elutions for cardiac stress tests, bone scans, and thyroid imaging procedures. Her generator was calibrated Monday morning at 6:00 AM with 50 GBq of ⁹⁹Mo activity. To optimize patient scheduling and minimize radioactive waste, Elena must calculate the optimal elution time to maximize ^99mTc yield. Using the decay chain calculator with t_1/2,P = 66 hours and t_1/2,D = 6.01 hours, she determines that maximum ^99mTc activity (approximately 44 GBq) occurs at 22.9 hours post-elution, corresponding to 5:00 AM Tuesday morning. This timing allows her to prepare radiopharmaceuticals for the day's first procedures starting at 7:00 AM, with each dose calibrated for administration 2-3 hours later when patients arrive. The transient equilibrium activity ratio of 1.10 confirms that daughter activity will briefly exceed parent activity, validating her elution schedule. This precise timing optimization reduces wasted radioisotope by 15% compared to arbitrary 24-hour elution cycles, saving the hospital approximately $12,000 annually in generator costs while ensuring adequate activity for all scheduled procedures.

Scenario: Environmental Radon Assessment in Home Inspection

Marcus, a certified home inspector specializing in radon testing, conducts a 5-day alpha track detector measurement in a basement of a 1970s split-level home scheduled for sale. His detector measures 7.8 pCi/L (288 Bq/m³) of radon-222, exceeding the EPA action level of 4 pCi/L. To advise the homebuyer on mitigation costs and health risks, Marcus needs to calculate the equilibrium factor between radon gas and its short-lived progeny (²¹⁸Po, ²¹⁴Pb, ²¹⁴Bi, ²¹⁴Po), which deliver the actual lung dose. Using the decay chain calculator with ²²⁶Ra parent half-life (1,600 years) and ²²²Rn daughter half-life (3.82 days), he confirms secular equilibrium exists, meaning radon concentration remains essentially constant over the measurement period. For the progeny, he calculates time to 99% equilibrium using the shortest half-life (²¹⁸Po at 3.1 minutes): t_eq = 6.64 × 3.1 min = 20.6 minutes. This confirms that radon progeny reach full equilibrium within 30 minutes of radon entry, validating the standard assumption that progeny contribute the full equilibrium dose. Marcus's report recommends a sub-slab depressurization system (estimated cost $1,200-1,800) and explains that without mitigation, the homeowner faces a lifetime lung cancer risk increase of approximately 2.7% based on EPA dose-response models, justifying the mitigation investment for the concerned buyers.

Scenario: Nuclear Reactor Post-Shutdown Dose Rate Prediction

Dr. Yuki Tanaka, a reactor physicist at a research reactor facility, prepares a detailed maintenance schedule requiring personnel entry into the reactor pool area 72 hours after shutdown from 4 MW steady-state operation. To ensure worker doses remain below 2 mSv per maintenance session, she must calculate the gamma dose rate from fission product decay chains, particularly ¹⁴⁰Ba/¹⁴⁰La (barium-lanthanum, with half-lives of 12.75 days and 1.678 days respectively). At shutdown, her ORIGEN-S burnup calculation predicts 3.2 × 10¹⁵ Bq of ¹⁴⁰Ba activity. Using the decay chain calculator in transient equilibrium mode with initial daughter activity of zero, Dr. Tanaka calculates that ¹⁴⁰La activity reaches 95% of its equilibrium value (activity ratio = 1.15 times parent activity) within 11.3 hours post-shutdown. At t = 72 hours, parent ¹⁴⁰Ba activity has decayed to 2.48 × 10¹⁵ Bq, while daughter ¹⁴⁰La has reached transient equilibrium at 2.85 × 10¹⁵ Bq. Combined with shielding calculations accounting for the 1.596 MeV gamma from ¹⁴⁰La, she determines the dose rate at the planned work location will be 1.7 mSv/hr, requiring work restrictions to 70-minute maximum shifts to remain within dose limits. This calculation leads her to recommend delaying maintenance until 120 hours post-shutdown, when dose rates drop to 0.85 mSv/hr, allowing standard 2-hour work periods and completing the maintenance campaign more efficiently with lower worker doses.

What is the difference between secular and transient equilibrium in decay chains? +

Secular equilibrium occurs when the parent half-life vastly exceeds the daughter half-life (typically by a factor of 10 or more), such that the parent decay rate remains essentially constant over several daughter half-lives. In secular equilibrium, the daughter activity approaches equality with parent activity after approximately 7 daughter half-lives, and both subsequently decay at the parent's rate. Examples include ²²⁶Ra/²²²Rn (1,600 years/3.82 days) and ²³⁸U/²³⁴Th (4.47 billion years/24.1 days). Transient equilibrium develops when the parent half-life exceeds the daughter half-life but by less than an order of magnitude. Here, daughter activity builds to a maximum exceeding parent activity by the ratio λ_D/(λ_D - λ_P), then declines at the parent's decay rate. The ⁹⁹Mo/^99mTc system (66 hours/6.01 hours) exemplifies transient equilibrium with a maximum activity ratio of 1.10. The mathematical distinction lies in whether λ_P can be neglected compared to λ_D in the Bateman equation denominator—when negligible, secular equilibrium approximations apply; when significant but smaller, transient equilibrium governs the behavior.

How do I handle decay chains where the parent and daughter have very similar half-lives? +

When parent and daughter half-lives are nearly identical (λ_P ≈ λ_D), the standard Bateman equation becomes numerically unstable due to near-zero division in the denominator (λ_D - λ_P) and near-equal exponential terms that subtract to produce large relative errors. For cases where |λ_P - λ_D| / λ_D falls below approximately 0.01 (1% difference), use the limiting form derived by L'Hôpital's rule: A_D(t) = A_D,0 e^-λt + A_P,0 λt e^-λt, where λ represents the common decay constant. This formulation contains a term linear in time, indicating that daughter activity initially increases linearly before exponentially decaying. This situation occasionally occurs in isomeric transition decay chains where nuclear excited states and ground states have similar lifetimes. Practically, if your calculator produces wildly fluctuating results or negative activities when parent and daughter half-lives differ by less than 1%, switch to this special-case equation or use higher-precision numerical methods (such as arbitrary-precision arithmetic) to minimize catastrophic cancellation errors in floating-point subtraction operations.

Why does my Mo-99/Tc-99m generator produce less activity after several elutions? +

Generator yield decreases over time due to parent (⁹⁹Mo) radioactive decay, which proceeds with a 66-hour half-life regardless of daughter elutions. Each elution removes daughter activity but does not affect the parent, so subsequent daughter regeneration depends on the remaining parent activity according to the Bateman equation. After one parent half-life (66 hours or approximately 2.75 days), parent activity drops to 50% of initial calibration, correspondingly reducing maximum daughter yield to 50%. After three parent half-lives (about 8 days), only 12.5% of original parent activity remains, producing insufficient ^99mTc for most clinical applications. Most generators are replaced after 7-10 days of use when parent decay reduces yields below practical thresholds. Additionally, repetitive elutions may leave trace residual daughter activity (incomplete elution), which slightly reduces net activity available in subsequent elutions—this effect typically accounts for 2-5% yield reduction and depends on elution technique and column chemistry efficiency. Breakthrough of parent ⁹⁹Mo into elution (alumina column deterioration) can also occur but is typically maintained below 0.1% by proper generator design and regular quality control testing.

How long must I wait for radon progeny to reach equilibrium when measuring indoor radon? +

Radon-222 progeny (the decay sequence ²¹⁸Po → ²¹⁴Pb → ²¹⁴Bi → ²¹⁴Po → ²¹⁰Pb) reach radioactive equilibrium with radon gas very rapidly, governed by the shortest half-life in the chain. ²¹⁸Po (polonium-218), the immediate daughter of ²²²Rn, has a half-life of 3.1 minutes. To reach 99% equilibrium requires approximately 6.64 half-lives, or about 20.6 minutes for ²¹⁸Po. The subsequent daughters ²¹⁴Pb (26.8 minutes) and ²¹⁴Bi (19.9 minutes) reach their respective equilibria within 3 hours. However, physical factors significantly affect apparent equilibrium in real environments: radon progeny are charged ions that plate out onto surfaces, attach to airborne particles, or are removed by ventilation at rates competing with radioactive production. The equilibrium factor F (ratio of actual progeny concentration to theoretical equilibrium concentration) typically ranges from 0.3 to 0.7 in homes, averaging around 0.4. For radon measurements using alpha track detectors (measuring radon gas directly) or continuous radon monitors, equilibration time refers to detector response time (typically hours to days) rather than radioactive equilibrium. For progeny-measuring devices like working level meters, allow 3-4 hours in closed-room conditions for physical-radioactive equilibrium establishment.

Can I use decay chain calculations for three-generation decay series? +

Yes, the Bateman equation extends to arbitrary-length decay chains through recursive application, though analytical solutions become algebraically complex beyond two generations. For a three-member chain (parent → daughter → granddaughter), the granddaughter activity equation contains three exponential terms with coefficients involving products of decay constants and differences. The general form is A₃(t) = C₁e^-λ₁t + C₂e^-λ₂t + C₃e^-λ₃t, where coefficients C_i depend on all three decay constants and initial activities. For practical calculations involving long chains (such as natural uranium and thorium decay series with 14+ members), numerical methods using matrix exponential techniques or stepwise integration provide more reliable results than extended analytical formulas. Modern reactor physics codes like ORIGEN employ these matrix methods to track hundreds of interconnected decay chains simultaneously. For manual calculations, two-generation (parent-daughter) equations often suffice because equilibrium establishment in early chain members makes subsequent generations proportional to their immediate parents, allowing segmented analysis. If you need three-generation calculations, specialized nuclear engineering software or iterative application of two-generation equations (treating each daughter as a parent for the next generation) provides accurate results with manageable computational complexity.

What units should I use for half-lives when parent and daughter have vastly different timescales? +

Always convert all time quantities (half-lives, elapsed times) to a common unit before performing decay chain calculations to avoid unit mismatch errors that produce incorrect results by orders of magnitude. For parent-daughter systems with disparate half-lives, choose a time unit that keeps both decay constants within a reasonable numerical range (neither extremely large nor extremely small) to minimize floating-point precision loss. For example, in ²³⁸U/²³⁴Th calculations (4.47 × 10⁹ years versus 24.1 days), expressing both in years yields λ_U = 1.55 × 10^-10 yr^-1 and λ_Th = 10.5 yr^-1, spanning 11 orders of magnitude and risking computational issues. Instead, convert to days: λ_U = 4.25 × 10^-13 day^-1 and λ_Th = 0.0288 day^-1, still spanning wide range but more computationally stable. For ⁹⁹Mo/^99mTc (66 hours/6.01 hours), hours work well for both. When in doubt, use seconds as the fundamental SI unit—all radioactive decay calculations inherently involve exponential functions that are unitless, so proper unit conversion in decay constant calculation automatically ensures correct results regardless of input time unit choice. Modern calculators handle unit conversion internally, but manual calculations require explicit attention to consistent time units throughout all equation terms.

Free Engineering Calculators

Explore our complete library of free engineering and physics calculators.

🔗 Explore More Free Engineering Calculators

• Criticality K Effective Interactive Calculator
• Fuel Burnup Mwd Mtu Interactive Calculator
• Steam Generator Sizing Nuclear Interactive Calculator

Browse All Engineering Calculators →

About the Author

Robbie Dickson — Chief Engineer & Founder, FIRGELLI Automations

Robbie Dickson brings over two decades of engineering expertise to FIRGELLI Automations. With a distinguished career at Rolls-Royce, BMW, and Ford, he has deep expertise in mechanical systems, actuator technology, and precision engineering.

Wikipedia · Full Bio