A multi-orthogonal polynomials' approach to bulk queueing theory

We consider a stationary Markov process that models certain queues with a bulk service of a fixed number $m$ of admitted customers. We find an integral expression of its transition probability function in terms of certain multi-orthogonal polynomials with respect to a system of distributions that contain measures supported on starlike subsets of the complex plane. We also study the corresponding steady state.