Interpolation Macdonald polynomials and Cauchy-type identities

Let Sym denote the algebra of symmetric functions and $P_\mu(\,\cdot\,;q,t)$ and $Q_\mu(\,\cdot\,;q,t)$ be the Macdonald symmetric functions (recall that they differ by scalar factors only). The $(q,t)$-Cauchy identity $$\sum_\mu P_\mu(x_1,x_2,\dots;q,t)Q_\mu(y_1,y_2,\dots;q,t)=\prod_{i,j=1}^\infty\frac{(x_iy_jt;q)_\infty}{(x_iy_j;q)_\infty}$$ expresses the fact that the $P_\mu(\,\cdot\,;q,t)$'s form an orthogonal basis in Sym with respect to a special scalar product $\langle\,\cdot\,,\,\cdot\,\rangle_{q,t}$. The present paper deals with the inhomogeneous \emph{interpolation} Macdonald symmetric functions $$I_\mu(x_1,x_2,\dots;q,t)=P_\mu(x_1,x_2,\dots;q,t)+\text{lower degree terms}.$$ These functions come from the $N$-variate interpolation Macdonald polynomials, extensively studied in the 90's by Knop, Okounkov, and Sahi. The goal of the paper is to construct symmetric functions $J_\mu(\,\cdot\,;q,t)$ with the biorthogonality property $$\langle I_\mu(\,\cdot\,;q,t), J_\nu(\,\cdot\,;q,t)\rangle_{q,t}=\delta_{\mu\nu}.$$ These new functions live in a natural completion of the algebra Sym. As a corollary one obtains a new Cauchy-type identity in which the interpolation Macdonald polynomials are paired with certain multivariate rational symmetric functions. The degeneration of this identity in the Jack limit is also described.