Matrix models for $\beta$-ensembles from Nekrasov partition functions

We relate Nekrasov partition functions, with arbitrary values of $\epsilon_1,\epsilon_2$ parameters, to matrix models for $\beta$-ensembles. We find matrix models encoding the instanton part of Nekrasov partition functions, whose measure, to the leading order in $\epsilon_2$ expansion, is given by the Vandermonde determinant to the power $\beta=-\epsilon_1/\epsilon_2$. An additional, trigonometric deformation of the measure arises in five-dimensional theories. Matrix model potentials, to the leading order in $\epsilon_2$ expansion, are the same as in the $\beta=1$ case considered in 0810.4944 [hep-th]. We point out that potentials for massive hypermultiplets include multi-log, Penner-like terms. Inclusion of Chern-Simons terms in five-dimensional theories leads to multi-matrix models. The role of these matrix models in the context of the AGT conjecture is discussed.