The MBIS implementation is great! Would it be possible to add the ability to partition the charge into core and valence components like they do in the original paper? This is essential to use the "valence width" parameter that is currently calculated for any electrostatic potential or electrostatic interaction energy calculation that includes a charge penetration correction!