Fractional radioactive decay law and Bateman equations

The aim of this work is to develop the fractional Bateman equations, which can model memory effects in successive isotopes transformations. Such memory effects have been previously reported in the alpha decay, which exhibits a non-Markovian behavior. Since there are radioactive decay series with consecutive alpha decays, it is convenient to include the mentioned memory effects, developing the fractional Bateman Equations, which can reproduce the standard ones when the fractional order is equal to one. The proposed fractional model preserves the mathematical shape and the symmetry of the standard equations, being the only difference the presence of the Mittag-Leffler function, instead of the exponential one. This last is a very important result, because allows the implementation of the proposed fractional model in burnup and activation codes in a straightforward way. Numerical experiments show that the proposed equations predict high decay rates for small time values, in comparison with the standard equations, which have high decay rates for large times. This work represents a novelty approach to the theory of successive transformations, and opens the possibility to study properties of the Bateman equation from a fractional approach.