Global eigenvalue distribution regime of random matrices with an anharmonic potential and an external source

We consider ensembles of random Hermitian matrices with a distribution measure determined by a polynomial potential perturbed by an external source. We find the genus-zero algebraic function describing the limit mean density of eigenvalues in the case of an anharmonic potential and a diagonal external source with two symmetric eigenvalues. We discuss critical regimes where the density support changes the connectivity or increases the genus of the algebraic function and consequently obtain local universal asymptotic representations for the density at interior and boundary points of its support (in the generic cases). The investigation technique is based on an analysis of the asymptotic properties of multiple orthogonal polynomials, equilibrium problems for vector potentials with interaction matrices and external fields, and the matrix Riemann-Hilbert boundary value problem.