The energy of the system is expanded as a sum of 2-body and 3-body contributions using cubic B-splines, which combine the beneficial properties of smoothness and differentiability with the advantage of a compact support. Hence, the number of basis functions that need to be evaluated in every energy computation step is strongly limited, as a maximum of four functions can be non-zero in every segment. The low number of basis functions directly relates to a high computational efficiency .