Home Download Help Forum Resources Extensions FAQ NetLogo Publications Contact Us Donate Models: Library Community Modeling Commons Beginners Interactive NetLogo Dictionary (BIND) NetLogo Dictionary User Manuals: Web Printable Chinese Czech Farsi / Persian Japanese Spanish

NetLogo Models Library: 
If you download the NetLogo application, this model is included. You can also Try running it in NetLogo Web 
This model simulates the spontaneous decay of a collection of radioactive nuclei. As they decay and become stable, the plot of the number that are still radioactive demonstrates the notion of "halflife".
At each time tick, each undecayed (light blue) nucleus has a certain chance of decaying. When a nucleus decays, it briefly flashes bright yellow (as if giving off radiation), then turns dark blue. Eventually, if you wait long enough, all of the nuclei will have decayed and the model will stop.
Set the initial number of nuclei (NUMBERNUCLEI slider) and the likelihood of decay during each time interval (DECAYCHANCE slider). Then push the SETUP button. Push the GO button to run the model.
The number of radioactive nuclei that remain is shown in the "Radioactive Nuclei" plot. Each time the number of nuclei is reduced by half, red and green lines appear on the plot to mark the place where each halfway mark was reached. The "Decay Rate" plot shows the number of decays that occur during each clock tick.
What is the shape of the decay curve (Radioactive Nuclei)? How is this affected by the initial conditions?
Why is the decay curve this shape? Is it the same as the decay curve shown in books?
The time between red lines is called the halflife. What is its physical meaning? Is it constant as the number of nuclei decreases? Is it affected by the initial number of nuclei or the decaychance? Do you think it's a useful way to characterize a radioactive material?
How long does it take for all the nuclei to decay?
Watch one nucleus carefully. When will it decay?
Radioactivity depends on the number of decays per unit time (Decay Rate), because each decay event gives off radiation. What happens to the radioactivity of this sample over time? What is the shape of the decay rate curve? How is it related to the shape of the decay curve?
Examine the standard equation for nuclear decay:
> N = No (e exp (T/tau))
No is the initial number of nuclei, N is the number at a later time T, and tau is the "mean lifetime". Compare its behavior to what you see in this model. The corresponding equation for radiation is:
> R = Ro (e exp (T/tau))
Why are these two equations so similar?
Try the extremes of the initial conditions: many or few nuclei, high or low decaychance. How does this affect the "jaggedness" of the decay rate plot?
What does the plot do when very few nuclei are left?
What instructions would you give each turtle be to make it behave like an unstable nucleus? Check the code in the Code tab and see if it's what you thought it should be.
Carbon14 dating involves comparing the ratio of a stable (C12) to an unstable (C14) nucleus. Explain that method in terms of this model.
It is often the case in nature that two nuclei with different decay rates are present in the same sample. Modify this model to have two or more types of nuclei. What happens to the radiation curve?
In this model, the nuclei don't affect their neighbors when they decay  they just disappear. Get them to emit particles that in turn cause reactions in other nuclei. See if you can model a chain reaction. (See the Reactor Top Down and Reactor XSection models.)
If you mention this model or the NetLogo software in a publication, we ask that you include the citations below.
For the model itself:
Please cite the NetLogo software as:
Copyright 1997 Uri Wilensky.
This work is licensed under the Creative Commons AttributionNonCommercialShareAlike 3.0 License. To view a copy of this license, visit https://creativecommons.org/licenses/byncsa/3.0/ or send a letter to Creative Commons, 559 Nathan Abbott Way, Stanford, California 94305, USA.
Commercial licenses are also available. To inquire about commercial licenses, please contact Uri Wilensky at uri@northwestern.edu.
This model was created as part of the project: CONNECTED MATHEMATICS: MAKING SENSE OF COMPLEX PHENOMENA THROUGH BUILDING OBJECTBASED PARALLEL MODELS (OBPML). The project gratefully acknowledges the support of the National Science Foundation (Applications of Advanced Technologies Program)  grant numbers RED #9552950 and REC #9632612.
This model was converted to NetLogo as part of the projects: PARTICIPATORY SIMULATIONS: NETWORKBASED DESIGN FOR SYSTEMS LEARNING IN CLASSROOMS and/or INTEGRATED SIMULATION AND MODELING ENVIRONMENT. The project gratefully acknowledges the support of the National Science Foundation (REPP & ROLE programs)  grant numbers REC #9814682 and REC0126227. Converted from StarLogoT to NetLogo, 2001.
(back to the NetLogo Models Library)