adding critical radius formula + its Python bindings