Orthogonal polynomials for a class of measures with discrete rotational symmetries in the complex plane

We obtain the strong asymptotics of polynomials $p_n(\lambda)$, $\lambda\in\mathbb{C}$, orthogonal with respect to measures in the complex plane of the form $$e^{-N(|\lambda|^{2s}-t\lambda^s-\overline{t\lambda}^s)}dA(\lambda),$$ where $s$ is a positive integer, $t$ is a complex parameter and $dA$ stands for the area measure in the plane. Such problem has its origin from normal matrix models. We study the asymptotic behaviour of $p_n(\lambda)$ in the limit $n,N\to\infty$ in such a way that $n/N\to T$ constant. Such asymptotic behaviour has two distinguished regimes according to the topology of the limiting support of the eigenvalue distribution of the normal matrix model. If $0<|t|^2<T/s$, the eigenvalue distribution support is a simply connected compact set of the complex plane, while for $|t|^2>T/s$ the eigenvalue distribution support consists of $s$ connected components. Correspondingly the support of the limiting zero distribution of the orthogonal polynomials consists of a closed contour contained in each connected component. Our asymptotic analysis is obtained by reducing the planar orthogonality conditions of the polynomials to an equivalent system of contour integral orthogonality conditions. The strong asymptotics for the orthogonal polynomials is obtained from the corresponding Riemann--Hilbert problem by the Deift--Zhou nonlinear steepest descent method.