Mayer cluster expansion is an important tool in statistical physics to evaluate grand canonical partition functions. It has recently been applied to the Nekrasov instanton partition function of N = 2 4d gauge theories. The associated canonical model involves coupled integrations that take the form of a generalized matrix model. It can be studied with the standard techniques of matrix models, in particular collective field theory and loop equations. In the first part of these notes, we explain how the results of collective field theory can be derived from the cluster expansion. The equalities between free energies at first orders is explained by the discrete Laplace transform relating canonical and grand canonical models. In a second part, we study the canonical loop equations and associate them to similar relations on the grand canonical side. It leads to relate the multi-point densities, fundamental objects of the matrix model, to the generating functions of multi-rooted clusters. Finally, a method is proposed to derive loop equations directly on the grand canonical model. It is crucial for our considerations that f is independent of . In this way, we exclude a class of models more relevant to the study of Nekrasov partition functions. For instance, setting = γ 2 , one recovers the model proposed by J. Hoppe in [27]. This model is a one-parameter version of the Nekrasov partition function that depends on two Ω-background equivariant deformation parameters 1 and 2 [28,29]. As → 0, it exhibits a phenomenon referred 1 Except for the first (planar) equation, which is a Ricatti equation, therefore non-linear. It is equivalent to a Schrödinger equation, i.e. a linear differential equation of second order.2 Such a structure should be related to the invariance of Nekrasov partition functions under transformations representing the SHc algebra uncovered in [23] (see also [24,25]). 3 In the case of real singularities, a prescription should be given to move away the poles from the contour by a small imaginary shift.