We propose a novel framework for approximations to intractable probabilistic models. The method is based on a free energy formulation of inference and allows for a simultaneous computation of marginal