Holonomic systems of Gegenbauer type polynomials of matrix arguments related with Siegel modular forms

Differential operators on Siegel modular forms which behave well under the restriction of the domain are essentially intertwining operators of the tensor product of holomorphic discrete series to its irreducible components. These are characterized by polynomials in the tensor of pluriharmonic polynomials with some invariance properties. We give a concrete study of such polynomials in the case of the restriction from Siegel upper half space of degree 2n to the product of degree n. These generalize the Gegenbauer polynomials which appear for n = 1. We also describe their radial parts parametrization and differential equations which they satisfy, and show that these differential equations give holonomic systems of rank 2 ⁿ.