We have described a new, open-source implementation of the straightforward MBE(n) method without any kind of electrostatic embedding. The asceticism is by design, as we aim to circumvent complexities arising from self-consistent charge embedding including overly complicated analytic gradients, 66 instead using high-order nbody expansions and low-level supersystem corrections to capture cooperative polarization. Numerical problems arising from unfavorable combinatorics for n > 3, which we have previously documented in large-scale MBE(n) calculations, are avoided here by means of energy screening. Stability in the presence of diffuse basis func-tions is achieved by iterating the n-body self-consistent field calculations to convergence.