Innuclear physics, theBateman equation is amathematical model describing abundances and activities in adecay chain as a function of time, based on thedecay rates and initial abundances. The model was formulated byErnest Rutherford in 1905[1] and the analytical solution was provided byHarry Bateman in 1910.[2]
If, at timet, there are atoms of isotope that decays into isotope at the rate, the amounts of isotopes in thek-step decay chain evolves as:
(this can be adapted to handle decay branches). While this can be solved explicitly fori = 2, the formulas quickly become cumbersome for longer chains.[3] The Bateman equation is a classicalmaster equation where the transition rates are only allowed from one species (i) to the next (i+1) but never in the reverse sense (i+1 to i is forbidden).
Bateman found a general explicit formula for the amounts by taking theLaplace transform of the variables.
(it can also be expanded with source terms, if more atoms of isotope i are provided externally at a constant rate).[4]

While the Bateman formula can be implemented in a computer code, if for some isotope pair,catastrophic cancellation can lead to computational errors. Therefore, other methods such asnumerical integration or thematrix exponential method are also in use.[5][6]
For example, for the simple case of a chain of three isotopes the corresponding Bateman equation reduces to
Which gives the following formula for activity of isotope (by substituting)
{{cite web}}: CS1 maint: archived copy as title (link)