## Abstract

The Gibbs distribution for n particles of the Coulomb gas on the unit circle at inverse temperature β is given by $Enβ(f)=(1/Zn,β)∫···∫f(eiθ1,…,eiθn)|Δ(eiθ1,…,eiθn) |β(dθ1/2π)…dθn/2π$ for any symmetric function f, where Δ denotes the Vandermonde determinant and Zn the normalization constant. We will describe an ensemble of (sparse) random matrices whose eigenvalues follow this distribution. Our approach combines elements from the theory of orthogonal polynomials on the unit circle with ideas from a recent work of Dumitriu and Edelman. In particular, we resolve a question left open by them: finding a tridiagonal model for the Jacobi ensemble.