Modified cutoff behavior; decay goes with erfc(r/cutoff), and neighbors go up to cutoff * 2