Radiochemistry inside a homogeneous volume element of a liquid phase irradiated with a steady electron beam is simulated under the assumption of an isotropic distribution of all contributing species. In this voxel, the chemical interplay is described by a set of coupled ordinary differential equations (ODEs) , describing the evolution of concentration c over time t of a reactant i depending on the concentrations of different reactants l and n (see the supporting information S1.1 for details):