![Arxiv:2011.12917V3 [Hep-Th] 30 Jan 2021 Introduction 1 Mti Oeso H Eodkn) Ntebt Ae,Oecnpro Can One Cases, Both the in Kind first Kind)](http://data.docslib.org/img/8d93312fa0f516f0cdd88741589a0b60-1.webp)
arXiv:2011.12917v3 [hep-th] 30 Jan 2021 Introduction 1 mti oeso h eodkn) ntebt ae,oecnpro can one cases, both the In kind first kind). the second τ of the models all of (matrix of models action functions (matrix the in Generating potentials eigenvalues. arbitrary with or matrices over integrals nepii obntra oe u 4 ] hspoet smc iia t similar much is property This term. the 5]. [4, hence [6], sum systems power combinatorial explicit an hsppr u oli oetn omlto otescn id So kind. second the to formulation extend to is goal our paper, this ucin fK/oaitgal hois e 13 o eal,we details, for [1–3] see theories, integrable KP/Toda of functions ∗ † nteoiia ento,mti oesaedfie saeae ove averages as defined are models matrix definition, original the In oee,bsdsbigjs a just being besides However, h ueitgaiiyhsbe xlctyfruae ofrfrtem • the for far so formulated explicitly been has superintegrability The [email protected] [email protected] [email protected]; cu ucin pern ntedsrpino pnHurwit spin of description the in appearing functions Schur ale eut n[,1],wihaecaatr flna groups: linear of characters are which [4,5, 10]), expressed [9, are in model results complex earlier rectangular the in Correlators model complex Rectangular eeaiaint h ope esrmdl ntelte ca latter the In mo model. Kontsevich tensor cubic complex the the and to model generalization complex of version special udai n olm oetasaedsigihdaogo featur among distinguished stronger grability are a potentials Coulomb is and calculability quadratic explicit This averages. h matrix The n hi eeaiain,adsprnerblt ok l looks superintegrability and generalizations, their and ymti function symmetric a eew d w oeeape fpicpliprac,where importance, principal Hermit of of examples cases more important two most add the we in Here happen to known already r( tr vrYugdarm ttelevel the at diagrams polynomials. Young Schur over the about special is what explains formula the and spae yteegnauso h matrix the of eigenvalues the by played is n()adblw l h nerto esrsaenraie ounit to normalized are measures integration the all below, and (1) In mtnraiainitgasi eoiaost ipiyformulas. simplify to denominators in integrals normalization omit X ayegnau arxmdl oss euirbsso ob of basis peculiar a possess models matrix eigenvalue Many ueitgaiiyo oteihmti model matrix Kontsevich of Superintegrability X ¯ ) sapclaiyo arxmdl,terlvn ai sfor is basis relevant the models, matrix of peculiarity a As . k
e se per χ X R { c r( Tr eeis here nttt o nomto rnmsinPolm,Mso 1 Moscow Problems, Transmission Information for Institute r oeivle hntoeo hrcesadatal oti an contain actually and characters of those than involved more are X nriMironov Andrei X ¯ N χ ) k R 1 } ( × a
z i N := eee hsc nttt,Mso 191 Russia 119991, Moscow Institute, Physics Lebedev ftevariables the of ) τ 2 fnto,tepriinfnto fmti oe a eotnpre often be can model matrix of function partition the -function, Z etnua matrix, rectangular d χ IT ogpun,110,Russia 141701, Dolgoprudny, MIPT, | R b R { TP ocw171,Russia 117218, Moscow ITEP, | r( Tr o xmlso h orsodn ae-airfrua,se[ see formulas, Harer-Zagier corresponding the of examples for , X a,b,c X ¯ X ) Abstract z and , ∗ k i X · } ¯ rapolynomial a or , hl tter.h.s. the at while , 1 exp χ k property a ike R hnodnr nerblt,js ietecssof cases the like just integrability, ordinary than e { lxiMorozov Alexei − p ubr n te eae contexts. related other and numbers z e h eeatcaatr r h celebrated the are characters relevant the se, k hrcnrlptnil,adw alit call we and potentials, central ther Tr } e.I h omrcs,srihfradi a is straightforward case, former the In del. steShrplnma wihi endt be to defined is (which polynomial Schur the is X e yteShrplnmas(characters) polynomials Schur the by med a,uiay n ope arxmodels. matrix complex and unitary, ian, X h oe eed on depends model the ¯ eta h atto ucin r actually are functions partition the that ve evbe hthv xlctycalculable explicitly have that servables orltr r hndsrbdb integrals by described then are correlators ontne hmi h rsn text. present the in them need not do a h nw xmlsinclude: examples known the far r lentvl,wt akrudfields background with alternatively, or, ) d χ ,]truhteShrplnmas(see polynomials Schur the through 7,8] hrce > character < 2 R arxesmls fe ecie by described often ensembles, matrix r X ti oeso h rtkn,ad in and, kind, first the of models atrix { p k p = y, ueitgaiiyo mechanical of superintegrability o } k 79,Russia 27994, ftepwrsums power the of = < χ ttelhs 1 h oeof role the (1) l.h.s. the At R d,b,c 1 N { 1 χ N > , R † 1 ∼ =1 nohrwrs we words, other in 1; := N } { δ character 2 χ k, diinlsummation additional xenlfields external R or 1 } { N δ MIPT/TH-16/20 FIAN/TD-08/20 k, ITEP/TH-27/20 2 IITP/TH-18/20 } 1 superinte- vrgsof Averages . hsis This . p k := etdas sented a : Q P i z 11]. (1) i k z ), i For the quantity χR{δk,1}, there is a hook formula [13]
lR 1 i [[n]]s,a = n if n = a mod(s)and1otherwise (4) which will be also useful beyond Gaussian models. Also simple in terms of the box coordinates is the expression on the topological locus: χR{N} = χR{δk,1} (N + α − β) (5) (α,βY)∈R Using eq.(1) and the Cauchy formula pk k k exp Tr (XX¯) = χR{pk}· χR{Tr (XX¯) } (6) k ! Xk XR one can immediately obtain a combinatorial expression for the generating function of all correlators: p χ {p } χ {N } χ {N } exp −Tr XX¯ + k Tr (XX¯)k d2X = R k R 1 R 2 (7) k ! χR{δk,1} Z Xk XR • Gaussian Hermitian model Correlators in the Gaussian Hermitian model are [4,5,8] (cf. also with a Fourier expansion of [14, Eq.(2.18)] in terms of characters) 1 χ {N} χ {δ } χ {Tr Xk} := χ {Tr Xk}· exp − Tr X2 dX = R R k,2 (8) R R 2 χ {δ } Z R k,1 where X is an N × N matrix, and the generating function of all correlators is 1 p χ{p }χ {N} χ {δ } exp − Tr X2 + k Tr Xk dX = k R R k,2 (9) 2 k ! χR{δk,1} Z Xk XR • Trigonometric (unitary type) matrix model This model is given by the eigenvalue integral, and the correlators are [15] N N 2 m − m m κ hχ (emi )i := χ (emi ) sinh2 i j exp − i dm = A|R| · q2 R · χ {p∗} (10) R R 2 2g2 i R i 2 • Knot matrix model Further deformation of the measure to the knot eigenvalue model gives basically the same expression for the correlators [16]: N 2 m m m b/a i i mi − mj mi − mj − i |R| 2κ ∗ χ (e a ) := χ (e a ) sinh sinh e 2g2 dm = A · q R · χ {p } (11) R R 2a 2b i R i g2 with q = e 2ab . mi The correlators of hχR(e )i are more involved this time, but are related to the HOMFLY-PT polynomials Torusa,b of the torus knots HR (A, q) [16]: mi Torusa,b Torusa,b hχR(e )i = HR (A, q)= χR{N}· HR (A, q) (12) where κ m b|R| 2b Q Torusa,b i a a a ∗ HR (A, q)= cR,QχQ(e ) = A cR,Q · q · χQ{p } (13) * + QX⊢a|R| QX⊢a|R| where cR,Q are peculiar Adams coefficients [17]. Generalization to non-torus knots, for example, like that in [18], remains to be analyzed. • q,t deformed Gaussian Hermitian model The eigenvalue q,t-deformed Gaussian Hermitian model is associated with a deformation of Schur polyno- mials: all of them in the formula (8) are replaced by the corresponding Macdonald polynomials MR [19]: z Nk 1+(−1)k 1 i 1−t z ; q MR k MR k j ∞ (N−1) logq t 2 2 2 1−t 1−t hMR(zi)i := MR(zi) zi (q zi ; q )∞dqzi = (14) zi n o 1 n o −1 t ; q M k Z j6=i zj i R 1−t Y ∞ Y n o where ∞ k (z; q)∞ := (1 − q z) (15) kY=0 and the integral is defined to be the Jackson integral, 1 ∞ k k k dqzf(z) := (1 − q) q f(ξq )+ f(−ξq ) (16) −1 Z kX=0 2 2 2 1 2 The measure dµq,t(z) = (q z ; q )∞dqz gives rise to the Gaussian measure dµG(z) = exp − 2 z dz in the limit of q → 1. One can expect also a further extension of similar character identities to the elliptic q,t model [20]. • Monomial non-Gaussian models Superintegrability is in no way restricted to the Gaussian integration measures or their deformations, though beyond them the knowledge is still restricted. The most important example is provided by mono- mial non-Gaussian actions: k k −Tr Xs/s χR{Tr X } := ⊗N χR{Tr X }e dX = a Cs,a R = χR{δk,s}· (α,β)∈R[[N + α − β]]s,0 · [[N + α − β]]s,a (17) Q N =0 or a mod s 3 where we used the notation from (4). The new trick here is the use of a special star-like (closed) integration contour Cs,a: s ∞ s s F (x) e−x /s dx := e−2πi(a−1)b/s · F e2πib/sx e−x /s dx (18) ZCs,a b=1 Z0 X which picks up only the powers of x, which are equal to a − 1 mod s, in particular, s k +1 xke−x /sdx = δ(s) · Γ (19) k+1−a s ZCs,a (s) δk is defined to be 1 if k = 0 mod s and to vanish otherwise. This makes the answer depending on an additional parameter a = 0,...s − 1. The r.h.s. of (17) contains some factors N + j from χR{N}, (5), that is, those with N + i = 0,a mod(s), and, hence, its vanishing depends also on the value of N. Note that if the condition N =0 or a mod s is not satisfied, one can not define the correlator by the condition < 1 >= 1 because of zeroes in the denominator. Note also that the Vandermonde determinant in the integrand is independent of s and a. For more details, see the original paper [21]. In this paper, we extend these results to the simplest background field models: the matrix integral over N1 × N2 complex matrices X with the measure −tr AXBX¯ 2 µC (X) ∼ e d X (20) where N1 × N1 matrix A and N2 × N2 matrix B are fixed external matrices, and the Hermitian generalized Kontsevich model [22] in background matrix field Λ, i.e. the matrix integral over N × N Hermitian matrix X with the measure tr(W (X)+ΛX) µK (X) ∼ e dX (21) and Λ is an external N × N matrix. In this paper, we consider these two examples, the potential in the second X3 one being W (X)= − 3 , which is just the original Kontsevich model [23]. Remarkably, the difference between these two models is severely increased in this case: in particular, in the Hermitian case, only one background field can be easily handled. The question we address is what are the functions of X-variables that have simple and explicitly calcu- lable averages. As we demonstrate in this paper, in the complex case, these are still the Schur polynomials k χR{tr (XX¯) }, however, the r.h.s. of (1) now contain traces of A and B instead of N1 and N2: −k −k k χR{tr A }χR{tr B } χR{Tr (XX¯) } = (22) χR{δk,1} see sec.2 below for details, and for generalization to the complex tensor model. In the Hermitian (Kontsevich) case, the story makes a new twist: the relevant functions are more restricted Q Schur polynomials, see sec.3.3 for the correlators in this case, −k k QR/2{Tr Λ }QR/2{δk,1} QR{Tr X } = δR,even · (23) QR{δk,1} and sec.4.2 for the combinatorial expression for the Kontsevich model −k 1 QR{Tr Λ }QR{δk,1}Q2R{δk,3} ZK = |R| (24) 4 Q2R{δk,1} RX∈SP Contributing here are only diagrams R with all lines of even length. Eq.(24) can be considered as the main new result of this paper. Important feature of formulas (22)-(24) is that they depend on the matrix size only implicitly, through the traces of powers of background fields. This is the crucial feature, which allows one to forget about the matrix-integral origin/realization of these models, in particular to treat them in terms of τ-functions of integrable hierarchies [22]. 4 2 Rectangular complex model 2.1 Correlators in terms of permutations We start with calculating the correlators in the complex matrix model with the external matrix (see [24] for a square matrix example), and use the symmetry group technique worked out in [5]. In this section, we closely follow that paper. We consider the rectangular complex model 2 ZC := d X exp −Tr AXBX¯ (25) Z with correlators defined as hO(X)i := O(X)exp −Tr AXBX¯ d2X (26) Z where X is N1 ×N2 rectangular matrix, A and B are square matrices of sizes N1 ×N1 and N2 ×N2 respectively, and (notice that, with this definition, the matrix multiplication is defined by convolution of the first indices with each other, and of the second indices with each other) ¯ i j ¯ kl Tr AXBX := AkXij Bl X (27) Xijk In order to calculate the pair correlator is equal to kl −1 −1 Xij X¯ = A B (28) ik jl The 2m-point correlator can be labelled by a permutation of indices σ that belongs to the symmetric group Sm: l m ¯ mp ¯ aiασ(i) hOσi = Tr (XX) = Xaiαi X (29) *p=1 + *i=1 + Y Y In order to calculate this correlator, we apply the Wick theorem and use formula (28): m m bγ(i) βγ(i) ¯ biβi −1 −1 Xaiαi X = A B (30) ai αi *i=1 + i=1 Y γX∈Sm Y and obtain A B hOσi = pγ pγ◦σ (31) γX∈Sm where lγ A −k B −k A,B A,B pk := Tr A , pk := Tr B , pγ := pi(γ) (32) i=1 Y and lγ is the number of cycles in the permutation γ, with i(γ) denoting the i-th cycle in the permutation γ. Now we use the standard identity [13] pγ = ψR(γ)χR{p} (33) RX⊢|γ| Here the sum goes over all Young diagrams R with |γ| boxes, χR{p} is the Schur polynomial (the character of lγ the linear group GL(N)), pγ = i=1 pγi and ψR(γ) is the character of the symmetric group S|R|. Now using the orthogonality relation Q ψ (σ) ψ (γ)ψ (γ ◦ σ)= ψ (γ−1)ψ (γ ◦ σ)= R δ (34) R Q R Q d QR γ γ R X X where dR := χR{δk,1}, we finally come to A A χR{p }χR{p } hOσi = ψR(σ) (35) χR{δk,1} RX⊢m where the sizes of R and σ coincide. 5 2.2 Complex model in background fields Let us now use the Frobenius formula [13] 1 χ {p } = ψ (γ)p (36) R k |R|! R γ γ X where |γ| is the size of the Young diagram R, and the orthogonality relation (34) with the particular value |R| ψR(id)= ψR([1 ]) = |R|!dR. Then, one immediately obtains from (35) −k −k k χR{tr A }χR{tr B } χR{Tr (XX¯) } = (37) χR{δk,1} This generalizes the answer (1) to the case of A 6= I, B 6= I. Similarly to (7), we also can obtain the combinatorial expression for the generating function p χ {p }χ {tr A−k}χ {tr B−k} exp −Tr AXBX¯ + k Tr (XX¯)k d2X = R k R R (38) k ! χR{δk,1} Z Xk XR 2.3 Tensor model in background fields The results of this section for the complex matrix model can be straightforwardly extended to the Gaussian tensor model in background fields. To this end, we follow [25] and again apply the technique of [5]. The Gaussian tensor model is a model of complex r-tensors Xa1,...ar where each runs trough its interval ai =1,...,Ni, in the background of r external square matrices A(i) of the sizes Ni with the Gaussian action N1 Nr r ai ¯ b1,...br S := ... Xa1,...ar X A(i) (39) bi 1 1 r r i=1 a X,b =1 a X,b =1 Y Gauge invariant operators at the level m in any (not obligatory Gaussian) tensor model are the tensorial counterparts of “multi-trace” operators N1 Nr m a1 ,...,ar 1 r ¯ σ1(p) σr(p) Oσ1,...,σr = ... Xap,...ap X (40) 1 r p=1 ! ~aX=1 ~aX=1 Y labeled by the set of m elements σi of the permutation group Sm. In fact, the labeling is reduced to a double ⊗r coset Sm\Sm /Sm [26, 27], but we do not need these details here. There is also a distinguished set of operators called generalized characters in [25] that are defined as 1 χ (X, X¯)= ψ (σ ) ...ψ (σ ) · O (41) R1,...,Rr n! R1 1 Rr r σ1,...,σr σ1X,...,σr ∈Sn These generalized characters do not form a full basis in the space of all gauge invariant operators, instead they form an over-complete basis in the space of all gauge invariant operators with non-vanishing Gaussian averages [25]. b br b1 ,...br −1 1 −1 The pair correlator is now X 1 r X¯ p p = A ··· A and, in complete analogy with what ap,...ap (1) (r) a1 ar we were doing in the complex matrixD model case,E we obtain the averag es of these generalized characters CR1,...,Rr −k −k χR1,...,Rr = · χR1 Tr A(1) · ... · χRr Tr A(r) (42) dR1 · ... · dRr D E n o n o where r i=1 ψRi (∆) CR1,...,Rr := (43) z∆ ∆X⊢n Q mk and z∆ := k k mk! is the standard symmetric factor of the Young diagram (order of the automorphism). In the case of r = 3, CR ,R ,R are the Clebsch-Gordan coefficients of the three irreducible representations R1, R2, Q 1 2 3 R3 of the symmetric group. 6 3 Correlators in Hermitian model 3.1 Examples of correlators Let us calculate the correlators in the Gaussian Hermitian matrix model with external fields Since the matrices are square, we return to the standard convention of convolution of indices in the matrix multiplication. One could expect that, similarly to coming from the rectangular complex model to the Gaussian Hermitian one, the results will not be too much different from those in sec.2. It turns out, however, not to be the case, and both the calculations in the Hermitian case with the external field are much more tedious even in the case of one external matrix, and the results remarkably differ from those of sec.2. We consider here only the case of one external matrix, i.e. the correlators of the form hF (X)i := F (X)exp −Tr X2Λ dX (44) Z Since, the partition function is an invariant function of the external matrix Λ, i.e. depends only on its traces, −k one suffices to choose Λ diagonal. Similarly to (32), we define pk := Tr Λ . In this model, the pair correlator in terms of eigenvalues λi of the (diagonal) matrix Λ is equal to 2δilδjk hXij Xkli = (45) λk + λl Here is the difference with the case of two external matrices: in that case one can not work in terms of eigenvalues, and the pair correlator is much more involved. Now we consider examples of simple correlators in order to demonstrate the structure of answers. 3 3 First consider < Tr X ·Tr X >= P 2 < Tr X >=2ξ1 < Tr X · Tr X >= p1 4 1 < Tr X >=8 + p2 (λi + λj )(λi + λk) i,j,kX 3 2 < Tr X · Tr X >=3p1 2 2 2 < Tr X · Tr X >=4ξ1 +8ξ2 2 < Tr X · Tr X · Tr X >=2p1ξ1 +2p2 2 < Tr X · Tr X · Tr X · Tr X >=3p1 7 6 1 1 < Tr X >= 16 + 24 + p3 +6p1p2 (λi + λj )(λi + λk)(λi + λl) (λi + λj )(λi + λk)(λk + λl) i,j,k,lX i,j,k,lX 5 3 < Tr X · Tr X >=5p3 + 10p1 4 2 1 1 < Tr X · Tr X >=8p1 + 64 2 + p1p2 + p3 (λi + λj )(λi + λk) (λi + λj ) (λi + λk) i,j,kX i,j,kX 1 1 < Tr X4 · Tr X · Tr X >=8p +8 +9p p 1 (λ + λ )(λ + λ ) λ2(λ + λ ) 1 2 i j i k i,j i i j i,j,kX X 3 3 3 < Tr X · Tr X >=3p3 + 12p1 3 2 2 < Tr X · Tr X · Tr X >=6p1ξ1 + 12p1p2 3 3 < Tr X · Tr X · Tr X · Tr X >= 3(2p3 +3p1) 2 2 2 3 < Tr X · Tr X · Tr X >=8ξ1 + 48ξ1ξ2 + 64ξ3 2 2 2 < Tr X · Tr X · Tr X · Tr X >=4p1ξ1 +8p1ξ2 +8p2ξ1 +8p3 2 2 < Tr X · Tr X · Tr X · Tr X · Tr X >=6p1ξ1 + 12p1p2 3 < Tr X · Tr X · Tr X · Tr X · Tr X · Tr X >= 15p1 (49) and all odd level correlators evidently vanish. These examples demonstrate that the correlators that involve traces of only odd degrees of the matrix X (they are boxed in (49)) are expressed only through the times (32), moreover, through the odd times. We will return to other correlators elsewhere, and here we consider only this type of correlators. We will need a set of polynomials of odd times that gives a complete basis and forms a closed ring. 3.2 Q Schur polynomials Emergency of only the odd times immediately gives one a hint to use the Q Schur polynomials instead of the standard Schur polynomials, as we did in the previous sections. These polynomials1 depend only on odd time- variables p2k+1 and only on strict Young diagrams R = {r1 > r2 >...>rlR > 0} ∈ SP. They were introduced by I. Schur [29] in the study of projective representations of symmetric groups, and were later identified by I. Macdonald [30] with the Hall-Littlewood polynomials HLR{p} := MacR{p} at t = −1: q=0 l /2 2 R · HLR(t = −1) for R ∈ SP Q := (50) R 0 otherwise Macdonald’s observations were that HL R{p} for R ∈ SP depend only on odd time-variables p2k+1, and that R HLR{p} for R ∈ SP form a sub-ring, i.e. the Littlewood-Richardson coefficients NR1,R2 in the ring R HLR1 {p}· HLR2 {p} = NR1,R2 HLR{p} (51) R∈R1⊗R2 RX∈SP vanish for R/∈ SP, provided t = −1 and R1, R2 ∈ SP. Note that HLR{p} do not vanish for R/∈ SP, and then they can also depend on even p2k, thus the set of QR{p} is not the same as the set of HLR, it is a sub-set, and a sub-ring. There is also a manifest way to construct the Q Schur polynomial as a Pfaffian (see, e.g., [28, Eq.(74)]). The Q Schur polynomials form a system, which has very close properties to the standard Schur functions, only they form a basis in a subspace of time-variables. In particular, there is a counterpart of the Frobenius formula (36), for the Q Schur polynomials it states ΨR(∆) QR{pk} = p∆ (52) z∆ ∆X∈OP 1More detailed review of these polynomials can be found, e.g., in [28]. 8 where OP (odd partitions) is a set of Young diagrams with all lengths of lines odd, and ΨR(∆) are the characters of the Sergeev group [31,32]. As any characters, they satisfy the orthogonality conditions: ′ ΨR(∆)ΨR′ (∆) ΨR(∆)ΨR(∆ ) ′ ′ l = δRR , l = δ∆∆ (53) 2 ∆ z∆ 2 ∆ z∆ ∆X∈OP RX∈SP Actually relevant for the Q Schur polynomials is the restriction to odd times, i.e. the Young diagram ∆ in (52), l∆ which defines the monomial p∆ = i p∆i should have all the lines of odd length: ∆ ∈ OP. Therefore of crucial importance is the celebrated one-to-one correspondence between the sets of SP and OP. At last, the Q Schur polynomialsQ satisfy the Cauchy formula, which we write in a form similar to (6), because this is how it will be used in sec.4.2 below: 2k+1 k p2k+1Tr X QR{p}QR{Tr X } = exp (54) k +1/2 ! RX∈SP Xk 3.3 Combinatorial expression for correlators We are now ready to formulate a nice result for the correlators (44) in terms of the Q Schur polynomials. Let us denote through 2R the Young diagram produced from the diagram R by doubling all line lengths, and denote through R|2 the Young diagram with all even line lengths (in this case, R/2 denotes the Young diagram with all line lengths being half of those of R)2. Then (compare with [33]), −k Q {Tr Λ }Q {δk,1} R/2 R/2 if R|2 k QR{δk,1} QR{Tr X } = (55) 0 otherwise For the quantity QR{δk,1}, there is a counterpart of the standard hook formula (2): l R 1 Rk − Rm |R|− 2 QR{δk,1} =2 l (56) R R + R Rj !! k m j k p2k+1 p2l+1 = (k +1/2) · δk,l (58) D E the Q-polynomials are orthonormal: QR QR′ = δR,R′ (59) D E and the Cauchy formula (54) acquires an especially simple form with unit coefficients in the sum. 4 Kontsevich model Now we are ready to consider the Kontsevich model given by the integral [23] Tr X3 Z = exp − + TrΛ2X dX (60) 3 Z and it can be further generalized to non-cubic potentials [22], including the quadratic one, when the model is equivalent to the Hermitian matrix model [2]. In this paper, we concentrate on the original cubic case, and we denote the background field by Λ2 to simplify the formulas in sec.4.2. 2For instance, an example of R|2 is R = [6, 4, 2]. In this case, R/2= [3, 2, 1] and 2R = [12, 8, 4]. 9 4.1 Kontsevich model in the character phase The normalization of measure of the Kontsevich integral depends on the phase of the model [34]. In the character phase, the integral is understood as a formal series in positive powers of Tr Λk, and the generating function of the correlators is defined to be Tr X3 2 dX exp − 3 + TrΛ X Z = (61) ch X3 R dX exp −Tr 3 Since the cubic part does not depend on the angularR part of the matrix X, one can first perform the angular integration using the expansion [35] of the Itzykson-Zuber integral [36] 2 † χR{δk,1} 2k k exp Tr Λ UXU [dU]= · χR{Tr Λ }χR{Tr X } (62) χR{N} Z R X so that one remains with χR{δk,1} 2k k Zch = · χR{Tr Λ } χR{Tr X } (63) χR{N} R X k and then calculate the correlators χR{Tr X } for the star-like integration contours directly applying formulas like (17). Thus, in this phase, the expression in terms of characters is straightforward and nearly obvious. 4.2 Kontsevich model in the Kontsevich phase Much more interesting is the case of Kontsevich phase, where the character calculus is much less trivial. In this phase, the model is given by the integral Tr X3 2 2 dX exp − 3 + TrΛ X Z = exp − Tr Λ3 (64) K 3 dXexp(−Tr X2Λ) R and is understood as a formal series in TrΛ−k. TheR simplest way to deal with this integral is to shift the integration variable to the saddle point X → X + Λ so that the integral takes the form X3 2 dX exp −Tr 3 − Tr X Λ ZK = (65) R dXexp(−Tr X2Λ) and one calculates this integral perturbativelyR using the Cauchy formula (54), Tr X3 1 exp − = Q {−Tr Xk}· Q δ (66) 3 R R 2 k,3 RX∈SP Since, for even sizes of Young diagrams |R|, QR{−pk} = QR{pk} and averages are non-vanishing only for even sizes, one finally obtains the Kontsevich partition function: −k 1 k 1 QR{Tr Λ }QR{δk,1}Q2R{δk,3} ZK = |R| QR{Tr X } · QR {δk,3} = |R| (67) 2 4 Q2R{δk,1} RX∈SP RX∈SP where we used (55) for the average. Note that, because of the factor Q2R{δk,3}, this sum effectively runs over R of sizes divisible by 3. As expected in Kontsevich phase, this partition function depends on the matrix size N only through the variables Tr Λ−k. This is exactly the same phenomenon which we observed in (22). Derivation of the formula (67) could be one of the goals in [33], but at that time the Q Schur polynomials were even less known then now, and an explicit answer like this was unavailable. 10 4.3 Integrability and character expansion Note that expansion of the partition function in characters typically allows one to check if it is a τ-function of the KP hierarchy immediately by checking that the expansion coefficients satisfy the Pl¨ucker relations. For instance, this is the case for the partition functions (7) and (9). Indeed, any partition function of the form Zw = SR{p}SR{p¯}wR (68) XR with the function wR being the product wR = f(i − j) (69) i,jY∈R andp ¯k just arbitrary parameters, is a KP τ-function [7,37,38], since arbitrary Schur polynomial SR{p¯} satisfies the Pl¨ucker relations, and multiplying any solution to the Pl¨ucker relations by wR of the form (69) preserves the solution. Now using (5) and fixingp ¯k = N2, one obtains that (7) is of the form (68). Similarly, fixing pk = δk,1, one obtains that (9) is of the form (68). As we explained in [28, sec.8] (see also [39]), similarly looking at expansion in the Q Schur polynomials, one can expect that the partition function is a τ-function of the BKP hierarchy. Moreover, we proposed an example of expansion literally repeating formula (7) but with the Schur polynomials substituted by the Q Schur polynomials, though have not managed to construct a matrix model that possesses such an expansion. Hence, it is surprising to realize that the Kontsevich model, which is a τ-function of the KP hierarchy [22], not BKP have a similar character expansion in the Q Schur polynomials. This realization in terms of matrix model is especially important as the Q Schur polynomials are related to the spin Hurwitz numbers and to the corresponding cut-and-join operators [28]. Let us emphasize that, as soon as the Kontsevich partition function is a KdV τ-function [22], and the KdV τ-function depends on odd times only, this is quite natural to expand the Kontsevich partition function in the Q Schur polynomials, which, in variance with the ordinary Schur polynomials, depends only on odd times. One may ask why people usually consider expansions of KdV τ-functions into the ordinary Schur polynomials and not into the Q Schur polynomials (see, however, [40]), which could look more natural. The reason is that the standard technique of the KdV hierarchy as an embedding into the KP hierarchy uses the formalism of free fermions [41], hence the ordinary Schur polynomials, and it is only using the neutral fermions that gives rise to the Q Schur polynomials, but this type of fermions is associated with the BKP hierarchy [39, 42]. We will consider this issue in detail elsewhere. 5 Conclusion In this paper, we extended the superintegrability relation < character >∼ character to the second kind of matrix models, depending on background fields. This provides a very nice generalization for the complex matrix model, where dimensions of matrices are lifted to more general traces of background matrices. The best known example of the background-field model is, however, different: this is the Kontsevich model, which is rather analogous to the Hermitian, not complex model. In the Hermitian case, the simple formulas are known to get a little more involved: they include characters at a peculiar locus δk,s, where s is the power of the potential. Thus, it was an intriguing question what happens in background fields. Our result in this paper is for the ordinary Kontsevich model with cubic potential, and it is that the very set of characters is changed from the Schur polynomials to the very interesting set of Q Schur polynomials, and then their averages are given by a rather transparent formula (55). The formula is, however, somewhat non-trivial, because it requires a correlation between normalization of different Q Schur polynomials, in most other examples this does not matter because just a single character appears in the answer which is made from its values at three different loci. In the Kontsevich case, this is different, and this adds new colors to the yet-no-so-well-known subject of Q Schur polynomials. It is important that the superintegrability-related expansion (67) of the Kontsevich partition function in Q Schur polynomials is much simpler than the expansion through the ordinary Schur polynomials, which could seem natural from the point of view of KP integrability [3, 43]. Indeed, the coefficients of the latter expansion were found in [44, 45] and turned out to be quite complicated and do not expose any clear structure. Even the character nature of the coefficients is not evident in this KP induced expansion. Eq.(67) is obviously free from all such drawbacks, in full accordance with the expectations. It, however, remains a problem to lift (67) to generalized Kontsevich model [22], while the ordinary Schur expansion of them should be a direct lifting of [44, 45]. 11 Acknowledgements This work was supported by the Russian Science Foundation (Grant No.20-12-00195). References [1] A. Gerasimov, A. Marshakov, A. Mironov, A. Morozov, A. Orlov, Nucl.Phys. B357 (1991) 565-618 S. Kharchev, A. Marshakov, A. Mironov, A. Orlov, A. Zabrodin, Nucl.Phys. B366 (1991) 569-601 [2] S. Kharchev, A. Marshakov, A. Mironov, A. Morozov, Nucl.Phys. B397 (1993) 339-378, hep-th/9203043 [3] A. Morozov, Phys.Usp.(UFN) 37 (1994) 1; hep-th/9502091; hep-th/0502010 A. Mironov, Int.J.Mod.Phys. A9 (1994) 4355; Phys.Part.Nucl. 33 (2002) 537; hep-th/9409190 [4] A. Mironov, A. Morozov, Phys. Lett. B771 (2017) 503, arXiv:1705.00976 [5] A. Mironov, A. Morozov, Phys. Lett. B774 (2017) 210, arXiv:1706.03667 [6] A. Mishchenko, A. Fomenko, Funct. Anal. Appl. 12 (1978) 113 W. Miller Jr., S. Post, P. Winternitz, J. Phys. A46 (2013) 423001, arXiv:1309.2694 [7] A. Alexandrov, A. Mironov, A. Morozov, S. Natanzon, JHEP 11 (2014) 080, arXiv:1405.1395 [8] H. Itoyama, A. Mironov, A. Morozov, JHEP 1706 (2017) 115, arXiv:1704.08648 [9] C. Kristjansen, J. Plefka, G. W. Semenoff, M. Staudacher, Nucl.Phys. B643 (2002) 3-30, hep-th/0205033 [10] S. Corley, A. Jevicki, S. Ramgoolam, Adv.Theor.Math.Phys. 5 (2002) 809-839, hep-th/0111222 [11] J. Harer, D. Zagier, Invent.Math. 85 (1986) 457-485 C. Itzykson, J.-B. Zuber, Comm.Math.Phys. 134 (1990) 197-208 S.K. Lando, A.K. Zvonkin, Embedded graphs, Max-Plank-Institut f¨ur Mathematik, Preprint 2001 (63) [12] A. Morozov and S. Shakirov, JHEP 0904 (2009) 064, arXiv:0902.2627 [13] W. Fulton, Young tableaux: with applications to representation theory and geometry, LMS, 1997 [14] R. de Mello Koch, S. Ramgoolam, arXiv:1002.1634 [15] A. Mironov, A. Morozov, JHEP 1808 (2018) 163, arXiv:1807.02409 [16] M. Tierz, Mod. Phys. Lett. A19 (2004) 1365-1378, hep-th/0212128 A. Brini, B. Eynard, M. Mari˜no, Annales Henri Poincar´e. Vol. 13. No. 8. SP Birkh¨auser Verlag Basel, 2012, arXiv:1105.2012 [17] M. Rosso, V. F. R. Jones, J. Knot Theory Ramifications, 2 (1993) 97-112 X.-S. Lin, H. Zheng, Trans.Amer.Math.Soc. 362 (2010) 1-18 math/0601267 [18] A. Alexandrov, A. Mironov, A. Morozov, An. Morozov, JETP Letters 100 (2014) 271-278 (Pis’ma v ZhETF 100 (2014) 297-304), arXiv:1407.3754 A. Alexandrov, D. Melnikov, JETP Letters 101 (2015) 51–56, arXiv:1411.5698 [19] A. Morozov, A. Popolitov and S. Shakirov, Phys. Lett. B784 (2018) 342, arXiv:1803.11401 [20] A. Mironov, A. Morozov, arXiv:2011.02855 A. Mironov, A. Morozov, arXiv:2011.01762 [21] C. Cordova, B. Heidenreich, A. Popolitov, S. Shakirov, Commun. Math. Phys. 361 (2018) 1235, arXiv:1611.03142 [22] S. Kharchev, A. Marshakov, A. Mironov, A. Morozov, A. Zabrodin, Phys. Lett. B275 (1992) 311, hep-th/9111037 S. Kharchev, A. Marshakov, A. Mironov, A. Morozov, A. Zabrodin, Nucl. Phys. B380 (1992) 181, hep-th/9201013 12 [23] M. Kontsevich, Commun. Math. Phys. 147 (1992) 1 [24] E. Langmann, R. J. Szabo and K. Zarembo, JHEP 01 (2004) 017, arXiv:hep-th/0308043 [25] H. Itoyama, A. Mironov and A. Morozov, JHEP 1912 (2019) 127, arXiv:1909.06921 H. Itoyama, A. Mironov and A. Morozov, Phys. Lett. B802 (2020) 135237, arXiv:1910.03261 [26] J. Ben Geloun, S. Ramgoolam, arXiv:1307.6490 R. de Mello Koch, D. Gossman, L. Tribelhorn, JHEP, 2017 (2017) 011, arXiv:1707.01455 P. Diaz, S.J. Rey, arXiv:1706.02667 J. Ben Geloun, S. Ramgoolam, arXiv:1708.03524 [27] H. Itoyama, A. Mironov and A. Morozov, Nucl. Phys. B932 (2018) 52, arXiv:1710.10027 [28] A. Mironov, A. Morozov, S. Natanzon, Eur. Phys. J. C80 (2020) 97, arXiv:1904.11458 [29] I. Schur, J. Reine Angew. Math. 139 (1911) 155-250 [30] I.G. Macdonald, Symmetric functions and Hall polynomials, Second Edition, Oxford University Press, 1995 [31] A. Sergeev, Math. Sb. USSR, 51 (1985) 419–427 [32] M. Yamaguchi, J. Algebra 222 (1999) 301–327, math/9811090 A. Kleshchev, Linear and projective representations of symmetric groups, Cambridge Tracts in Mathematics 163, Cambridge Univ. Press (2005) [33] P. Di Francesco, C. Itzykson, J. B. Zuber, Commun. Math. Phys. 151 (1993) 193, hep-th/9206090 [34] A. Mironov, A. Morozov, G. W. Semenoff, Int. J. Mod. Phys. A11 (1996) 5031, hep-th/9404005 [35] A. Y. Morozov, Theor. Math. Phys. 162 (2010) 1 [Teor. Mat. Fiz. 161 (2010) 3], arXiv:0906.3518 [36] Harish-Chandra, Am. J. Math. 79 (1957) 87 C. Itzykson and J. B. Zuber, J. Math. Phys. 21 (1980) 411 [37] S. Kharchev, A. Marshakov, A. Mironov, A. Morozov, Int.J.Mod.Phys. A10 (1995) 2015, hep-th/9312210 [38] A. Orlov, D.M. Shcherbin, Theor.Math.Phys. 128 (2001) 906-926 A. Orlov, Theor.Math.Phys. 146 (2006) 183–206 [39] A. Orlov, Theor. Math. Phys. 137 (2003) 1574-1589, math-ph/0302011 J. Harnad, J.W. van de Leur, A.Yu. Orlov, Theor. Math. Phys. 168 (2011) 951-962, arXiv:1101.4216 J.W. van de Leur, A.Yu. Orlov, arXiv:1404.6076 ; arXiv:1611.04577 A.Yu. Orlov, T. Shiota, K. Takasaki, arXiv:1201.4518 ; arXiv:1611.02244 [40] J.J.C. Nimmo, A. Orlov, Glasgow Mathematical Journal 47 (A) (2005) 149-168, nlin/0405009 [41] E. Date, M. Jimbo, M. Kashiwara, T. Miwa, RIMS Symp. ”Non-linear integrable systems - classical theory and quantum theory” (World Scientific, Singapore, 1983) [42] E. Date, M. Jimbo, M. Kashiwara, T. Miwa, Physica D4 (1982) 343-365 M. Jimbo, T. Miwa, Publ. RIMS Kyoto Univ. 19 (1983) 943-1001 [43] Y.Ohta, J.Satsuma, D.Takahashi, T.Tokihiro, Prog.Theor.Phys.Suppl. 94 (1988) 210 [44] J. Zhou, arXiv:1306.5429 [45] F. Balogh and D. Yang, Lett. Math. Phys. 107 (2017) 1837, arXiv:1412.4419 13