A Radioactive Dice Decay Simulation is a classic educational model that demonstrates the probabilistic nature of nuclear half-life. Since the exact moment an unstable atomic nucleus will decay is entirely random, we can model a large group of atoms using a pile of six-sided dice.
In this simulation, each die represents an unstable nucleus, a specific roll result (like a 6) represents a decay event, and each round of rolling represents one unit of time. Once a die “decays,” it is permanently removed from the sample. Key Mathematical Concepts
Decay Constant (λ): The probability that an individual atom will decay in a given time interval. For a standard 6-sided die where rolling a 6 means decay,
Decay Formula: The expected number of remaining undecayed nuclei N(t) after t time steps (rolls) follows an exponential decay trend:
N(t)=N0(1−λ)tcap N open paren t close paren equals cap N sub 0 open paren 1 minus lambda close paren to the t-th power Where N₀ is the initial sample size. Theoretical Half-Life ( T1/2cap T sub 1 / 2 end-sub
): The time it takes for half of the sample to decay. It can be derived by setting , which simplifies to:
T1/2=ln(2)−ln(1−λ)=ln(2)−ln(⁄6)≈3.8 stepscap T sub 1 / 2 end-sub equals the fraction with numerator l n 2 and denominator negative l n open paren 1 minus lambda close paren end-fraction equals the fraction with numerator l n 2 and denominator negative l n open paren 5 / 6 close paren end-fraction is approximately equal to 3.8 steps Step-by-Step Python Implementation Guide
Instead of rolling hundreds of real dice by hand, you can build a programmatic simulation in Python. Below is the step-by-step breakdown to code it. 1. Import Libraries and Initialize Parameters
Start by importing Python’s built-in random module to handle the random number generation. Radioactive Dice Simulation Lab | PDF – Scribd
Leave a Reply