On generalized Macdonald polynomials

Generalized Macdonald polynomials (GMP) are eigenfunctions of specificallydeformed Ruijsenaars Hamiltonians and are built as triangular polylinear combinations of Macdonald polynomials. They are orthogonal with respect to a modified scalar product, which could be constructed with the help of an increasingly important triangular perturbation theory, showing up in a variety of applications. A peculiar feature of GMP is that denominators in this expansion are fully factorized, which is a consequence of a hidden symmetry resulting from the special choice of the Hamiltonian deformation. We introduce also a simplified but deformed version of GMP, which we call generalized Schur functions. Our basic examples are bilinear in Macdonald polynomials.


Introduction
Macdonald polynomials are getting increasingly important for practical calculations in string theory. The reason is their role in representation theory of the Ding-Iohara-Miki symmetry, which underlines dynamics of background brane networks in 6d super-Yang-Mills models. At the same time, their properties remain under-investigated, especially for the purposes of physical applications. This paper is one in the series of recent attempts to cure this situation. This time we concentrate on the subject of generalized Macdonald polynomials (GMP).
Macdonald polynomials M R labelled by Young diagrams R are [1] graded symmetric polynomials of variables x i , i = 1, . . . , N (in this case we use the notation M R (x i )), or graded polynomials of time-variables p k := N i=1 x k i (in this case we use the notation M R {p k }). In this paper, we will be mostly interested in the second point of view. They can be defined in many different ways, in particular, within the context of quantum toroidal algebras [2,3] (see also [4,5] and references therein). There are two other natural definitions: using their triangle structure and orthogonality condition w.r.t. a scalar product, JHEP01(2020)110 or using a Hamiltonian structure behind them. We will briefly review in section 2 both of these possibilities in order to demonstrate that these definitions are really effective.
Unfortunately, despite their undisputable importance, these GMP can be effectively described only using the quantum toroidal algebra behind them [47,54,55]. That is, the GMP are defined as common eigenfunctions of a deformed Calogero-Ruijsenaars Hamiltonian (generalized cut-and-join operators of [56][57][58]) that is read off from the quantum toroidal algebra [47,54,55]. The goal of the present paper is to look for a definition independent on this algebraic approach. We discuss the two possibilities: how the deformed Hamiltonian can be independently obtained and how the GMP can be defined through the Gauss decomposition of an appropriate scalar product.

Triangular structure
Let us fix the scalar product on time-variables and continue it to any polynomials by linearity. From no on, we use the notation {x} := x − x −1 . Here the Young diagram ∆ = δ 1 ≥ δ 2 ≥ . . . ≥ δ l ∆ , and p ∆ = l ∆ i=1 p δ i . The

JHEP01(2020)110
combinatorial factor z ∆ is best defined in the dual parametrization of the Young diagram, ∆ = . . . , 2 m 2 , 1 m 1 , then z ∆ = k k m k · m k !. Then, the ordinary Macdonald polynomials can be defined as a lower triangular combination of the Schur polynomials S R {p} orthogonal w.r.t. this scalar product (2.1): 3) The Young diagrams are ordered lexicographically:

JHEP01(2020)110
Ruijsenaars Hamiltonians in these variables as difference operators acting on the functions of N variables x i aŝ The first Hamiltonian in terms of time-variables. The Ruijsenaars Hamiltonians can be rewritten in terms of time-variables. The results reads as follows (we slightly redefine the Ruijsenaars Hamiltonians here to make formulas simpler). Introduce an operator Then, the first Hamiltonian isĤ with the eigenvalues when acting on the eigenfunction M R λ (1) In fact, in order to construct the Macdonald polynomials, one does not need to know all the Hamiltonians, it is sufficient to use only the first one: its polynomial solutions are unambiguously determined. The reason is the integrable structures behind the system. Nevertheless, in order to have the complete picture, we construct all higher Hamiltonians now.

JHEP01(2020)110
The general Hamiltonian looks likê (2.14) Symbolically, it can be rewritten aŝ 15) This expression becomes non-symbolic in the case of q = t: in this case, the pair product in this expression is a result of normal ordering so that one can re-writê This is a reflection of the fact that any 0 dz/z ·V k (z) is a Hamiltonian. This is not surprising, since the eigenfunctions in this case do not depend on t: these are the Schur functions. Hence, t is a free parameter and such is t k , expansion in any of these parameters generate the corresponding Hamiltonians having the Schur functions as their eigenfunctions. Unfortunately, in the general case of different q and t, the factor coming from the normal ordering is different from that in (2.14)-(2.15). For instance, Eigenvalues. One would better choose another normalization of the Hamiltonians: The eigenvalues in symmetric representations are: Similarly, in representation [n 1 , n 2 ] the eigenvalues are where

JHEP01(2020)110
In the general case, with the symmetric variables and the time-variables Another construction of Hamiltonians. There is another construction of the Hamiltonians that can be conveniently extended to the GMP. We borrow the formulas of this paragraph from [63]. Let us define an operatorP m := z=0 dz z m+1V 1 (z) (2.25) Then, the new first Hamiltonian coincides with the old one and is equal just toĤ 1 ∼Ĥ 1 = P 0 , the second one is the commutatorĤ 2 ∼ [P −1 , P 1 ] and all remaining Hamiltonians are repeated commutators 1Ĥ where we introduced the notation g k := {q k }{t −k }{t k /q k } and a non-unit pre-factor is chosen to simplify further formulas. The eigenvalues of these Hamiltonians are Λ (k) It can be understood within the context of the quantum toroidal algebra, since according to [2,3], these Hamiltonians are spectral dual (see also [64][65][66][67][68][69]) to ψ + k generators of the algebra and, hence, are commuting.

JHEP01(2020)110
where The generating function of these eigenvalues can be written in the form The Hamiltonians (2.14) are related with these by the triangle transformation One can also make another triangle transformation in order to construct new Hamiltonians H k , their generating function being The eigenvalues of these new Hamiltonians are L (k) k . This means that they celebrate the property Note that one can choose various generating functions for the Hamiltonians (2.26), e.g. that of the exponential type This generating function does not includeĤ 1 ∼P 0 .

Macdonald polynomials: summary
Thus, we explained that the Macdonald polynomials can be unambiguously defined either by the triangle expansion (2.2) with the requirement of the orthogonality condition w.r.t. the scalar product (2.1), or by the requirement that they are eigenfunctions of the Hamiltonian (2.10), which comes from the Ruijsenaars-Schneider integrable system. However, triangularity is not yet apparent from the shape of the Hamiltonian and this relation requires better understanding, see [62] for a possible approach to this problem.

JHEP01(2020)110
Because of integrability, in this system there are also higher Hamiltonians. No canonical way is still available to construct them from the first principles. We mention just three possible approaches.
One of the possibilities is the Hamiltonians (2.14), which have their eigenvalues (2.22) expressed through linear combinations of the antisymmetric Schur (or, which is the same, of antisymmetric Macdonald polynomials), the latter being symmetric polynomials of the variables (2.23). Let us note that the shifted (or interpolating) antisymmetric Macdonald polynomials [70] are linear combinations very similar to (2.22): The difference is only in an additional product in the denominator of the summand. Another way of choosing the Hamiltonians is (2.26), they have the eigenvalues (2.27), which are just antisymmetric Schur or Macdonald polynomials (without the shift). These polynomials are functions of four sets of time-variables HereR denotes the representation of SL N conjugate to R. The third approach is mentioned in [62]: the first Hamiltonian converts Schur functions into bilinear sums restricted to single-hook Young diagrams X and Y of same size. In higher Hamiltonians restriction is weakened to k-hook diagrams. This approach has chances of revealing the origins of triangularity and could allow generalization from the Schur polynomials not only to the Macdonald ones, but also to the Kerov functions, but there is long way to work all this out.

JHEP01(2020)110 3 Introduction to GMP
Let us now discuss if one can define the GMP using one of the two possibilities that we discussed in the previous section for the ordinary Macdonald polynomials.

Triangular structure of GMP
As reviewed in detail in recent [73], the basic property of all symmetric polynomials naturally emerging in the course of study of non-perturbative quantum field theory, is that they are triangular combinations of Schur functions. This property is a bonus from existence of the natural lexicographical ordering for Young diagrams, and actually it breaks down when one goes up from ordinary to plane partitions [74][75][76], but it is still true for finite sets of Young diagrams, and thus for the Generalized Macdonald and Kerov functions.
. . , A n } depends on an ordered sequence of n Young diagrams, on n infinite sets of time variables and on n "deformation" parameters A i , sometimes additionally constrained by the condition n i=1 A i = 1: since nothing especially interesting depends on n, we mostly consider the case of n = 2 in order to simplify formulas, with a single parameter A (usually substituted by k =p k , i.e. our typical notation will be M Similarly to the ordinary Macdonald polynomial case, we present the GMP as a triangular linear combination of the products of the corresponding Schur functions: Let us discuss the ordering in this sum. At each level L = |R 1 | + |R 2 |, the pairs of Young diagrams are ordered in such a way that R 2 is ordered first. The ordering starts from diagrams of smaller sizes and is lexicographical for the same size diagrams. Then, one similarly orders R 1 within diagrams with coinciding R 2 . For instance, at level 2: One could fix the ordering for the same size diagrams to be the transposed lexicographical one (reftranspolex), which would give essentially different results in Kerov generalizations [72] (see [73] for details), however, in the Macdonald case, the difference between these orderings is inessential as we explained in the previous section, and it remains inessential for the GMP. Moreover, it is sufficient to consider the partial ordering and, in all cases when the diagrams can not be ordered, the corresponding Kostka-Mac-

JHEP01(2020)110
Note that, from the triangular structure of the usual Macdonald polynomials, it follows that there is also an expansion (3.4) For the purposes of the present paper the most convenient is this expansion in the Macdonald polynomials, not in the Schur ones.
The GMP are defined so that at the beginning of the lexicographic sequence, when R 2 = ∅, they do not depend on Q = A 2 and coincide with the ordinary Macdonald polynomials, ∅,R {p,p} is a full-fledged function of both time variables and A, and has not much to do with M R {p}). Another reduction appears in the undeformed case, for A = ∞: The scalar product. Now one has to specify the transform, i.e. to explain how the coefficients C are defined. In [73] for the ordinary Macdonald and Kerov functions, we imposed an orthogonality condition, but for the GMP this is not immediate: one should know what the relevant scalar product is. Our goal is to construct such a scalar product, and thus to put the theory of GMP into the general context of [73], but this requires some work and insights. Before doing this, let us note that there is a standard scalar product for the GMP, however, the system of GMP's is not orthogonal, but bi-orthogonal w.r.t. it: one has to introduce a dual set of the GMP The standard product is an independent product of those for p k andp k , (2.1): The ordinary Macdonald functions are orthogonal in this scalar product, (3.4).) The problem is that the freedom is not fixed by such orthogonality requirement, and there is no guiding principle to select the needed triangular transformation.
Let us note that the transformation from the Schur to Macdonald polynomials is lower triangular both in the normal and dual sectors, while the further transformation to the GMP continues to be lower triangular in normal sector, but becomes upper triangular in the dual one. Thus the net transform from the Schur polynomials to the dual GMP's is not really triangular.

JHEP01(2020)110
One can see in this example that the GMP themselves can be obtained by a triangular transformation from the Schur bilinears, but it is not the case for the dual GMP. This is obvious already for the very first one: in variance with M[1, 1], ∅{p,p} = S [1,1] {p}, the dual M ∅, [2] does not coincide with the corresponding Schur M [2] {p}. As already mentioned, the reason for this difference is that the GMP's are obtained by composing two upper triangular transformation, but the dual GMP's are compositions of lower and upper triangular transformations (see comment after formula (3.10)), which are no longer triangular themselves.

GMP as eigenfunctions of deformed Hamiltonians
Since the GMP are not orthogonal in the standard metric, the Hamiltonians that have the GMP's their eigenfunctions are not Hermitian in the standard metric, though the deformation is triangular. The typical form is p ∂ ∂p , so that the ordinary Macdonald polynomials M R,∅ {p,p} = M R {p} are not affected.
In fact, it is important to add some additional Hermitian pieces, and they bring into the game the deformation parameter(s). For example, has eigenfunctions p andp − η u−ū · p with eigenvalues u andū. It becomes Hermitian if the scalar product is deformed to In general, u is a coefficient in front of np n ∂ ∂pn (or higher Macdonald Hamiltonians), and this term just shifts the eigenvalue of M R {p}, while the mixing term ∼ ǫ n ·n 2 p n ∂ ∂pn annihilates it.
In fact, it is sufficient to consider only one Hamiltonian, it fully determines the deformation. All the higher Hamiltonians are deformed in a consistent way.
Perturbation theory. The main point is that perturbation theory is greatly simplified for triangular perturbations. If the deformed Hamiltonian isĤ +V andĤψ i = λ i ψ i , then The underlined sums are triangular, if such is the perturbation matrix: V k i = 0 only for k < i. In this case, the eigenvalues do not change, and the eigenfunctions are finite sums with V entering the expression for the i-th eigenfunctions at most in the i-th power.
In application to the GMP, the role of ψ i is played by M RMQ , and the operatorV written symbolically asV = k V k p k ∂ ∂p k decreases the size of Q and enlarges that of R: which brings (R, Q) down in the ordering. The new eigenfunctions, i.e. exactly the GMP, are given by The eigenvalues actually have the form Perturbed eigenfunctions are no longer orthogonal in the original scalar product, where ψ i ψ j = δ i,j . However, one can introduce a new one, where they are, and triangularity JHEP01(2020)110 allows one to reformulate it as a simple recursion relation: In particular, for the lowest eigenfunction ψ 0 , which remains un-deformed, (3.24) This expansion is very similar to above one for the eigenfunction itself, the only differences are alternated signs and the "common eigenvalue" in denominators: it is λ i for the eigenfunctions and λ 0 for the deformed scalar product. In the case of GMP, the role of ψ 0 is played by any of the ordinaryp-independent Macdonald functions M R {p}.

GMP Hamiltonian. The actual Hamiltonian in the case of GMP is [50-55]
and perturbation is generated by the ǫ n -term, ǫ n := 1 − (t/q) 2n . Thus it looks more likê Explicit shifts of the arguments can be substituted by action of the operators like The eigenvalue of M R,S {p,p} remains un-deformed, i.e. is the same as it was for M R {p} · M S {p} with the un-deformed Hamiltonian (cf. with (2.11)): Higher Hamiltonians. Similarly to the case of ordinary Macdonald polynomials, the GMP's are unambiguously defined as graded polynomial eigenfunctions of the first Hamiltonian (3.25). This clearly implies an underlying integrable structure. It simultaneously implies that there exist higher Hamiltonians, and we again borrow formulas from [63] (see also [80, appendix B]

JHEP01(2020)110
in order to directly extended them to the GMP case: one just needs to substitute in all formulas P m witĥ with ǫ k := 1 − t 2k q 2k . Then, the Hamiltonians for the GMP are given by the same formulas (2.26), and the eigenvalues are given by formulas similar to (2.27), but depending on two Young diagrams:Ĥ (3.33) is expressed through the eigenvalues (2.27).
The generating function for these eigenvalues is given by An extension from the bilinear to multilinear GMP is evident. One can again introduce the HamiltoniansĤ One more notable thing is emergence of the triple of parameters (q 1 , q 2 , q 3 ) = (q, t −1 , tq −1 ), which is not seen at the level of the ordinary Macdonald functions. Thus the GMP are better suited for studying the apparent triality of the DIM algebra [77][78][79].

GMP: summary
Thus, we studied the Macdonald polynomial features that could be used for their definition and that can be extended to the GMP case. In particular, there is still a triangular structure of expansion, (3.1) in the Schur polynomials. Moreover, in this case, there is also a triangular expansion in Macdonald polynomials (3.4). However, there is no known scalar product so far that could fix the triangular expansion coefficients from the orthogonality condition.
There is a Hamiltonian (3.25) whose eigenfunctions are the GMP's, and it unambiguously defines them. However, this Hamiltonian at the moment comes only from an additional algebraic set-up of quantum toroidal algebras [47,54,55]. Moreover, in a complete analogy with the standard Macdonald case, one can construct higher Hamiltonians (3.31) with eigenvalues given by the same antisymmetric Schur or Macdonald polynomials. However, they are now polynomials of a sum of two sets of time-variables (3.33), each of them being associated its own Young diagram in accordance with the same Macdonald case formula (2.28).
Note that such a sum of times corresponding to two different Young diagrams emerges within the same context of refined topological vertex and Hopf hyperpolynomial as we discussed in section 2.3, however, in the case of composite representations (see [71,81], in particular, formula (32) in [71]).
In the next sections, we are going to modify the Macdonald scalar product so that the Q-deformed Hamiltonians become Hermitian. With this scalar product, the generalized Macdonald functions are themselves orthogonal, there is no need for the dual set of functions. If one could define (select) this deformed product by some appealing and simply formulated ansatz, this would fix the triangular transform and thus define the generalized Macdonald polynomials without any direct reference to the deformed Hamiltonians. We will look for this scalar product basing on the known set of GMP's, and then reverse the logic: we postulate it as a new definition of the GMP.
Our other goal in the forthcoming sections is to clarify reasons why the GMP Hamiltonian (3.25) is distinguished without references to the algebraic set-up. On this way, we will construct another set of symmetric functions of two sets of variables, which we call generalized Schur functions.

Looking for a scalar product
Let us discuss what could be the scalar product that gives rise to the correct GMP from the orthogonality relations. Let us denote

JHEP01(2020)110
Now we require the orthogonality of the GMP, In particular, 3) for all R and S at all levels. However, generally the proportionality coefficient in (4.2) is not known, though is unambiguously fixed by the triangular expansion (3.1) and by (4.2) at (R 1 , S 1 ) = (R 2 , S 2 ). Now let us proceed at the first two levels. At the first level, there is only one nondiagonal item of the symmetric matrix G R 1 ,R 2 |S 1 ,S 2 : This defines the scalar product completely. At level 2, there are more relations, but many more non-diagonal matrix elements of G R 1 ,R 2 |S 1 ,S 2 : and [1] Q−q 2 + +κ·G [2],∅ − κ·{t} q 2 {q}{q 2 }{qt}(Q−q 2 ) + t{q 2 /t} q 2 {q} 2 (Q−1) − q 8 t 2 +q 6 t 2 −q 4 t 4 −q 6 −q 4 t 2 +t 2 q 4 t{q}{q 2 }{qt}(Qq 2 −1) JHEP01(2020)110 The ugly-looking underlined terms can be eliminated or improved by choosing appropriate deformations of the norms: Then Now there are two natural versions of the other elements of the scalar product matrix Version A. Assume that the norms of un-deformed polynomials M R {p} also remain un-deformed: Then (4.9) However, the other components get very complicated. There is another choice, which provides simply looking expressions for all matrix elements. Now we describe it.
Version B. In fact, requiring the orthogonality of the generalized Macdonald polynomials, one obtains that the scalar product is defined up to arbitrary 3 components, say, , [1] [1], [1] (note that the components G ∅, [2] ∅, [2] , G ∅,[1,1] ∅, [1,1] always remain unfixed and are not taken into account). The choice with the simplest structure of zeroes and poles in Q looks as follows  [1,1] = −{q 2 /t}{q/t} 2 t 3 Y q 6 (Q−1) (4.12) Here we leave Y arbitrary, since the scalar product can be arbitrarily rescaled; one may require that G [1], [1] [1], [1] = {q} 2 {t} 2 as Q → ∞, which unambiguously fixes Y . Note that an additional pole at Q = 1 in the denominator emerge iff the second and the third representations coincide. This pole probably comes from the level 1 scalar product. The structure of poles is quite clear from these expressions. The only concern is the element G ∅, [2] ∅, [1,1] .
This search can be continued to higher levels, however, the calculations becomes much more involved, and the final answer is not available at the moment.

Understanding GMP Hamiltonians
We discuss now a specific property of the Hamiltonian (3.25) that makes its eigenfunctions much simpler, in particular, it leads to a cancellation of many denominators that are generally expected to emerge. One can say that the Hamiltonian is actually distinguished by this property.
Thus, we define non-trivial generalized Schur polynomials (GSP), which are a simpler version of GMP, but not their particular case unless one takes the limit (3.6). Looking at GSP can help to reveal and describe the essential properties of generic generalized polynomials, and this can be the first step to untie them from peculiarities of the Macdonald polynomials.
The Hamiltonian (6.3) at t = q unambiguously defines the GSP's, however, higher Hamiltonians are no longer given by the construction of section 3.2.

Cauchy formula vs. triangular expansion
The Cauchy formula. The Cauchy formula is always correct and is a direct corollary of the orthogonality relation. In this appendix, we demonstrate how it works for the GMP. We start, however, with the usual Macdonald polynomials. Then, the Cauchy formula has the form In order to prove this formula, one has to act with the operator on the both sides of this formula and then put all p k = 0. Then, the r.h.s. is trivially equal to M Q {p ′ }, while, at the l.h.s., one has to use the orthogonality condition (2.3) and realize the scalar product (2.1) by the differential operation in order to obtain The superscript ∨ here denotes the transposition of the Young diagram. This consideration demonstrates that, in order to obtain the Cauchy formula for the GMP's one has to use both M R 1 ,R 2 and M R 1 ,R 2 , due to the orthogonality relation (3.9). Then, the Cauchy formula acquires the form The transposition formula now looks like (7.6) and the Cauchy formula is The triangular structure. Let us now note that the Cauchy formula for the Schur polynomials because of the identity

JHEP01(2020)110
(which is straightforwardly lifted to Kerov-Kostka matrix with arbitrary Kerov couplings g k ). In order to prove eq. (7.10), note that the scalar product (2.1) of the Schur functions is equal to Thus, the inverse Macdonald-Kostka matrix diagonalizes the matrix µ R,R ′ with the diagonal elements being ||M R || 2 . We write it in the matrix form where M is the diagonal matrix with elements ||M R (q, t)|| 2 and the tilde denotes the matrix transposition. At the same time, the inverse matrix is diagonalized by the Macdonald-Kostka matrix itself: Now it remains to notice that ||M R (q, t)|| 2 = ||M R ∨ (t, q)|| −2 (7.16) Let us denote the interchanging q and t in a matrix with a bar and introduce the matrix Λ := δ R,R ∨ so that this latter identity reads as Then one obtains from (7.14)-(7.17) which is exactly (7.10). This consideration directly extended to the case of GMP and of more general Kerov functions. To better understand the reason for additional matrices Λ emerging in (7.18) and somewhat strangely looking formula (7.10), it is instructive to look at a straightforward generalization of the Macdonald polynomials, Kerov functions [72,73].

JHEP01(2020)110
Kerov functions. One can consider an extension of the scalar product (2.1) to (7.19) which induces the Kerov functions, 2 which are defined to be triangular transforms from the Schur functions orthogonal with the scalar product (7.19). The Young diagrams are again ordered here lexicographically. However, in this case, there are two systems of symmetric functions, since the Kostka-Kerov matrices K (g) R,R ′ between diagrams placed differently in two different orderings are non-zero: Since the orderings begin to differ at level 6, the two system of the Kerov functions become different at level 6.
The Cauchy formula in this case looks like (7.21) since the rule that relates the Kerov functions of transposed Young diagrams (a counterpart of (7.3)) is The identity (7.10) becomes, in the Kerov case, R,R ′′ (t, q) = δ R ′ R ′′ (7.23) Notice that the definition of K in (7.20) differs from K. In the Macdonald case, K is related with K by conjugation with the matrix Λ, which explains emergence of this latter in (7.18). Formula (7.23) is proved in the same way as in s. 7 with taking into account that ||Ker (g) R ∨ || −2 (7.24)

Conclusion
The goal of this paper was to understand the notion of generalized symmetric polynomials.
In the literature they are actually known in just one particular example of the generalized JHEP01(2020)110 Macdonald polynomials, which appear in the theory of AGT relations, where q, t-deformed Selberg-Kadell integrals of GMP presumably reproduce Nekrasov functions. The question, however, is to find more algebraic and, most important, generalizable definitions. Usually, one uses two approaches, which we described and discussed in the text. Both are based on the triangularity property, which is now understood to be the basic one for the actual theory of symmetric functions and their application to integrability in general and supersymmetric gauge theories in particular. All the relevant special functions in these fields are triangular combinations of Schur functions: triangular with respect to one or another kind of lexicographical ordering of Young diagrams. However, the triangularity per se is not enough. The two approaches consist of imposing orthogonality conditions either directly or through a requirement of forming the eigenfunctions of Hermitian Hamiltonians.
We discussed advantages and drawbacks of these two approaches in application to GMP.
The problem with straightforward orthogonality is the complexity (and ambiguity) of the relevant scalar product, which we did not even manage to construct in the full form. For naive scalar products, the GMP's are not orthogonal: one needs to introduce a complementary set of "dual" GMP, and then the question is to impose relations between the dual and original GMP.
The Hamiltonian approach in the case of GMP is more successful, and we advanced it further to explicit description of higher Hamiltonians. However, even the triangularity property is somewhat obscure in this case, even for the ordinary Macdonald functions, and this makes the construction of Hamiltonians and the very fact of their existence a kind of mysterious art, preventing making integrability from becoming a clear deductive piece of science. This is a serious drawback, if one looks for generalizations, which are desperately needed, because most partition functions in quantum field theory are related to still unknown integrable systems. Moreover, even in pure mathematics, Hamiltonians are not yet known for an immediate generalization from the Macdonald polynomials to the Kerov functions.
We hope that our detailed presentation of these issues, of many problems which remain unsolved, and of many "miracles" which can provide the clues for their resolution, will attract an attention and help to advance the subject in the near future.