This python driver intends to implement many body expansion techniques for calculating energies and properties of weakly solvated clusters using psi4 (and Gaussian) with proper accounting of basis set superposition errors using well known correction schemes lke mbcp, vmfc, etc.