If you're defining the potential, then presumably you know for Be
and the cutoff you are using, what is the max value of rho that
an atom can experience from its near neighbors. Its a physical
value based on how much "rho" each atom defines in
its near field. Since you define the rho() function, you also
define this max value, don't you?