Blog Archive

Saturday, September 18, 2021

09-18-2021-1229 - Bateman equation

In nuclear physics, the Bateman equation is a mathematical model describing abundances and activities in a decay chain as a function of time, based on the decay rates and initial abundances. The model was formulated by Ernest Rutherford in 1905[1] and the analytical solution was provided by Harry Bateman in 1910.[2]

If, at time t, there are  atoms of isotope  that decays into isotope  at the rate , the amounts of isotopes in the k-step decay chain evolves as:

(this can be adapted to handle decay branches). While this can be solved explicitly for i = 2, the formulas quickly become cumbersome for longer chains.[3] The Bateman equation is a classical master 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 the Laplace 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]

Quantity calculation with the Bateman-Function for 241Pu

While the Bateman formula can be implemented in a computer code, if  for some isotope pair, cancellation can lead to computational errors. Therefore, other methods such as numerical integration or the matrix exponential method are also in use.[5]

For example, for the simple case of a chain of three isotopes the corresponding Bateman equation reduces to



https://en.wikipedia.org/wiki/Bateman_equation

 



No comments:

Post a Comment