Detect correctly radial bonds changes during integration step

[#55563212]
Perform chemical reactions calculations every 50 timesteps

[#56525850]