On multivariate orthogonal polynomials and elementary symmetric functions

We study families of multivariate orthogonal polynomials with respect to the symmetric weight function in d variables B γ ( x ) = ∏ i = 1 d ω ( x i ) ∏ i < j | x i - x j | 2 γ + 1 , x ∈ ( a , b ) d , for γ > - 1 , where ω ( t ) is an univariate weight function in t ∈ ( a , b ) and x = ( x 1 , x 2 , … , x d ) with x i ∈ ( a , b ) . Applying the change of variables x i , i = 1 , 2 , … , d , into u r , r = 1 , 2 , … , d , where u r is the r -th elementary symmetric function, we obtain the domain region in terms of the discriminant of the polynomials having x i , i = 1 , 2 , … , d , as its zeros and in terms of the corresponding Sturm sequence. Choosing the univariate weight function as the Hermite, Laguerre, and Jacobi weight functions, we obtain the representation in terms of the variables u r for the partial differential operators such that the respective Hermite, Laguerre, and Jacobi generalized multivariate orthogonal polynomials are the eigenfunctions. Finally, we present explicitly the partial differential operators for Hermite, Laguerre, and Jacobi generalized polynomials, for d = 2 and d = 3 variables.