Abstract
we present new strategy to perform vibrational mean field and vibrational configuration interaction coalculations for large molecular systems. Thanks to an efficient pruning scheme based on the many body expansion of the potential energy surface, we show that variational bound state calculation is possible form large molecular systems