Meixner multiple orthogonal polynomials on interlacing lattices (Q6569646)

The authors consider the polynomials \(P_{n_1,n_2}\), which are orthogonal with respect to two discrete measures \N\[\N\mu_j(y)=\sum_{x\in\gamma_j+\mathbb Z_{+}}R(x)\delta(y-x), \qquad j=1,2, \N\]\Nwhere the weight function \(R\) is defined on lattices \(\gamma_j+\mathbb Z_{+}\) as a product of two classical Meixner measures, namely \N\[\NR(x):=R_{\gamma_1,\gamma_2}^{\alpha_1,\alpha_2,b}(x) =b^x \frac{\Gamma(x-\gamma_1+\alpha_1)\Gamma(x-\gamma_2+\alpha_2)}{\Gamma(x-\gamma_1+1)\Gamma(x-\gamma_2+1)}.\N\]\NHere \N\[\N0<b<1, \quad -1<\gamma_2-\gamma_1<1, \quad -\alpha_1<\gamma_2-\gamma_1<\alpha_2. \N\]\NThe main goal of this work is to prove that the system \(\{\mu_1, \mu_2\}\) of positive measures is perfect, i.e., the polynomial \(P_{n_1,n_2}\) is, up to multiplication by a constant, the unique polynomial of multiple orthogonality with respect to the pair of measures \({\mu_1, \mu_2}\) for the index \((n_1,n_2)\).