Matrix resolvent and the discrete KdV hierarchy

TL;DR

This paper introduces a matrix-resolvent approach to the discrete KdV hierarchy, defining tau-functions and deriving explicit formulas, with applications to ribbon graph enumeration and cubic Hodge integrals.

Contribution

It develops a novel matrix-resolvent method for the discrete KdV hierarchy and connects tau-functions with combinatorial and geometric applications.

Findings

Explicit formulas for tau-function derivatives are derived.

Connections between tau-functions and ribbon graph enumeration are established.

Applications to cubic Hodge integrals are demonstrated.

Abstract

Based on the matrix-resolvent approach, for an arbitrary solution to the discrete KdV hierarchy, we define the tau-function of the solution, and compare it with another tau-function of the solution defined via reduction of the Toda lattice hierarchy. Explicit formulae for generating series of logarithmic derivatives of the tau-functions are then obtained, and applications to enumeration of ribbon graphs with even valencies and to the special cubic Hodge integrals are considered.

Equations

$$P(n):=\Lambda +w_n{\Lambda }^{-1},$$

$$A_{\ell}\;:=\;\bigl{(}P^{\ell+1}\bigr{)}_{+}\,,\quad\ell\geq 0\,.$$

$$\frac{\partial P}{\partial s_{j}}\;=\;\bigl{[}A_{2j-1}\,,\,P\bigr{]}\,,\qquad j\geq 1\,.$$

$$\frac{\partial w_n}{\partial s_1}=w_n(w_{n+1}-w_{n-1}),$$

$$L:=P^2={\Lambda }^2+w_{n+1}+w_n+w_n{w}_{n-1}{\Lambda }^{-2}.$$

$$\frac{\partial L}{\partial s_{j}}\;=\;\bigl{[}A_{2j-1}\,,\,L\bigr{]}\,,\qquad j\geq 1\,.$$

$$\frac{\partial (w_{n+1}+w_n)}{\partial s_1}=w_{n+2}{w}_{n+1}-w_n{w}_{n-1},$$

$$\frac{\partial (w_n{w}_{n-1})}{\partial s_1}=(w_{n+1}+w_n-w_{n-1}-w_{n-2})w_n{w}_{n-1}.$$

$$L{\psi }_n=\lambda {\psi }_n,\text{mbox}i.e.{\psi }_{n+2}+(w_{n+1}+w_n-\lambda ){\psi }_n+w_n{w}_{n-1}{\psi }_{n-2}=0,$$

$$\frac{\partial {\psi }_n}{\partial s_j}=A_{2j-1}\psi .$$

$$\tilde{L}={\Lambda }^2+a_1(n)\Lambda +a_2(n)+a_3(n){\Lambda }^{-1}+a_4(n){\Lambda }^{-2},$$

$$\mathcal{L}:=\left(\begin{array}{cc}{\Lambda }^2& 0\\ 0& {\Lambda }^2\end{array}\right)+U_n,U_n:=\left(\begin{array}{cc}{w}_{n+1}+w_n-\lambda & w_n{w}_{n-1}\\ -1& 0\end{array}\right).$$

$$\mathcal{L}\left(\begin{array}{c}{\psi }_n\\ {\psi }_{n-2}\end{array}\right)=0.$$

$$R_{n+2}{U}_n-U_n{R}_n=0.$$

$$\displaystyle R_{n}\;=\;\begin{pmatrix}1&0\\ 0&0\\ \end{pmatrix}+\mathcal{O}\bigl{(}\lambda^{-1}\bigr{)}\,,$$

$$\mathrm{tr}{R}_n=1,\det R_n=0.$$

$$R_n(\lambda )=\left(\begin{array}{cc}1+{\alpha }_n(\lambda )& {\beta }_n(\lambda )\\ {\gamma }_n(\lambda )& -{\alpha }_n(\lambda )\end{array}\right).$$

$${\beta }_n=-w_n{w}_{n-1}{\gamma }_{n+2}$$

$${\alpha }_{n+2}+{\alpha }_n+1=(\lambda -w_{n+1}-w_n){\gamma }_{n+2}$$

$$(\lambda -w_{n+1}-w_n)({\alpha }_n-{\alpha }_{n+2})=w_n{w}_{n-1}{\gamma }_n-w_{n+2}{w}_{n+1}{\gamma }_{n+4}$$

$${\alpha }_n+{\alpha }_n^2+{\beta }_n{\gamma }_n=0$$

$$\displaystyle\sum_{i,j\geq 1}\frac{\partial^{2}\log\tau_{n}^{\textsc{\tiny\rm dKdV}}({\bf s})}{\partial s_{i}\partial s_{j}}\lambda^{-i-1}\mu^{-j-1}\;=\;\frac{{\rm tr}\,\bigl{(}R_{n}(\lambda)R_{n}(\mu)\bigr{)}-1}{(\lambda-\mu)^{2}}\,,$$

$$\displaystyle\frac{1}{\lambda}+\sum_{i\geq 1}\frac{1}{\lambda^{i+1}}\frac{\partial}{\partial s_{i}}\log\frac{\tau^{\textsc{\tiny\rm dKdV}}_{n+2}}{\tau^{\textsc{\tiny\rm dKdV}}_{n}}\;=\;\bigl{[}R_{n+2}(\lambda)\bigr{]}_{21}\,,$$

$$\frac{{\tau }_{n+2}^{\text{textsc}\mathrm{dKdV}}{\tau }_{n-1}^{\text{textsc}\mathrm{dKdV}}}{{\tau }_{n+1}^{\text{textsc}\mathrm{dKdV}}{\tau }_n^{\text{textsc}\mathrm{dKdV}}}=w_n.$$

$$e^{\alpha n+{\beta }_0+{\sum }_{k\ge 1}{\beta }_k{s}_k},$$

$$\sum_{j_{1},\dots,j_{k}=1}^{\infty}\frac{1}{\lambda_{1}^{j_{1}+1}\cdots\lambda_{k}^{j_{k}+1}}\frac{\partial^{k}\log\tau_{n}^{\textsc{\tiny\rm dKdV}}({\bf s})}{\partial s_{j_{1}}\dots\partial s_{j_{k}}}\;=\;-\frac{1}{k}\sum_{\sigma\in S_{k}}\frac{{\rm tr}\,\bigl{(}R_{n}(\lambda_{\sigma_{1}})\cdots R_{n}(\lambda_{\sigma_{k}})\bigr{)}}{\prod_{i=1}^{k}(\lambda_{\sigma_{i}}-\lambda_{\sigma_{i}+1})}\,,$$

$$w_n(\mathbf{s})=\frac{{\tau }_{n+1}(\mathbf{s}){\tau }_{n-1}(\mathbf{s})}{{\tau }_n^2(\mathbf{s})}.$$

$${\tau }_n(\mathbf{s})=e^{\alpha n+{\beta }_0+{\sum }_{k\ge 1}{\beta }_k{s}_k}{\tau }_n^{\text{textsc}\mathrm{dKdV}}(\mathbf{s}){\tau }_{n+1}^{\text{textsc}\mathrm{dKdV}}(\mathbf{s}).$$

$$\det\bigl{(}\mu_{i+j-2}\bigr{)}_{i,j=1}^{n}\;=\;\det\bigl{(}\mu_{2i+2j-2}\bigr{)}_{i,j=1}^{[n/2]}\det\bigl{(}\mu_{2i+2j-4}\bigr{)}_{i,j=1}^{[(n+1)/2]}\,.$$

$$\sum _{j_1,\dots ,j_k=1}^{\infty }\frac{\frac{{\partial }^k\log {\tau }_n(\mathbf{s})}{\partial s_{j_1}\dots \partial s_{j_k}}}{{\lambda }_1^{j_1+1}\cdots {\lambda }_k^{j_k+1}}=$$

$$\displaystyle\bigl{\langle}{\rm tr}\,M^{2j_{1}}\,\cdots\,{\rm tr}\,M^{2j_{k}}\bigr{\rangle}_{c}\;:=\;k!\,\sum_{0\leq g\leq\frac{|j|}{2}-\frac{k}{2}+\frac{1}{2}}\,n^{2-2g-k+|j|}\,a_{g}(2j_{1},\dots,2j_{k})\,,$$

$$a_g(2j_1,\dots ,2j_k):=\sum _{\Gamma }\frac{1}{\#\mathrm{Sym}\Gamma }.$$

$$E_{k}(n;\lambda_{1},\dots,\lambda_{k})\;:=\;\sum_{j_{1},\dots,j_{k}=1}^{\infty}\frac{\bigl{\langle}{\rm tr}\,M^{2j_{1}}\cdots{\rm tr}\,M^{2j_{k}}\bigr{\rangle}_{c}}{\lambda_{1}^{j_{1}+1}\cdots\lambda_{k}^{j_{k}+1}}\,.$$

$$R_n(\lambda ):=\left(\begin{array}{cc}1& 0\\ 0& 0\end{array}\right)+\sum _{j=0}^{\infty }\frac{(2j-1)!!}{{\lambda }^{j+1}}\left(\begin{array}{cc}(2j+1)A_{n,j}-(n-1)B_{n,j}& (n-n^2)B_{n+2,j}\\ & \\ B_{n,j}& (n-1)B_{n,j}-(2j+1)A_{n,j}\end{array}\right)$$

$$A_{n,j}:=(n-1){}_2{F}_1(-j,2-n;2;2),$$

Full text

*†*Deceased on March 19, 2019.

Keywords. Volterra lattice, matrix resolvent, tau-structure, modified GUE correlator, Hodge integral.

Mathematics Subject Classification (2010). Primary 37K10; Secondary 60B20, 39B05, 14N35.

1. 1 Introduction
2. 2 Basic formulation
3. 3 Tau-structure for the discrete KdV hierarchy
4. 4 Proof of Theorem 1
5. 5 Proofs of Theorem 2 and Corollary 2
6. A On consequence of the Hodge–GUE correspondence

1. Introduction

The discrete KdV equation (aka the Volterra lattice equation) is an integrable Hamiltonian equation in (1+1) dimensions, i.e. one discrete space variable and one continuous time variable, which extends to a commuting system of Hamiltonian equations, called the discrete KdV integrable hierarchy. This integrable hierarchy has important applications in algebraic geometry and symplectic geometry (in particular in the theory of Riemann surfaces) (see e.g. [22]). Significance of the discrete KdV hierarchy was further pointed out by E. Witten [39] in the study of the GUE partition function with even couplings — the “matrix gravity”, and was recently addressed also in the study of the special cubic Hodge partition function [15, 19, 20] — the topological gravity in the sense of [15, 20]. The explicit relationship between the two gravities, called the Hodge–GUE correspondence, has been established in [15, 20]. In this paper, by using the matrix-resolvent (MR) approach recently introduced and developed in [1, 2, 3, 14, 18] we study the tau-structure for the discrete KdV hierarchy, and apply it to studying the above mentioned enumerative problems.

1.1. The discrete KdV hierarchy

Let $P(n)$ be the following difference operator

[TABLE]

where $\Lambda$ denotes the shift operator $\Lambda:f_{n}\mapsto f_{n+1}$. Introduce

[TABLE]

Here, for an operator $Q$ of the form $Q=\sum_{k\in\mathbb{Z}}Q_{k}\,\Lambda^{k}$, the positive part $Q_{+}:=\sum_{k\geq 0}Q_{k}\,\Lambda^{k}$. The discrete KdV hierarchy is defined as the following system of commuting flows:

[TABLE]

For example, the $s_{1}$-flow reads

[TABLE]

which is the discrete KdV equation. The commutativity implies that equations (3) for all $j\geq 1$ can be solved together, yielding solutions of the form $w_{n}=w_{n}({\bf s})$, ${\bf s}:=(s_{1},s_{2},s_{3},\dots)$.

Let us introduce

[TABLE]

Then A_{2j-1}=\bigl{(}P^{2j}\bigr{)}_{+}=\bigl{(}L^{j}\bigr{)}_{+} .

Lemma 1**.**

The discrete KdV hierarchy (3) can be equivalently written as

[TABLE]

The proof will be given in Section 2. For the particular case $j=1$, we have

[TABLE]

It can be shown that equations (7)–(8) are equivalent to equation (4); the details for this equivalence are in Section 2.2.

Observe that equations (6) are the compatibility conditions of the following scalar Lax pairs:

[TABLE]

We want to write the spectral problem (9) into a matrix form. The scalar Lax operator $L$, defined in (5), could be viewed as a reduction of

[TABLE]

which is the Lax operator of a bigraded Toda hierarchy. However, observe that $L$ contains $\Lambda^{{\rm even}}$ only (with ${\rm even}=-2,0,2$). So, instead of considering a $4\times 4$ matrix-valued Lax operator, a $2\times 2$ matrix-valued operator will be sufficient. Indeed, introduce

[TABLE]

Then the spectral problem (9) reads

[TABLE]

1.2. The MR approach to tau-functions

In this subsection, we apply the MR approach to study further some basics in the theory of the discrete KdV hierarchy (in particular about tau-function), and will arrive at a formula for computing logarithm of the tau-function. Denote by ${\mathbb{Z}}[{\bf w}]$ the ring of polynomials with integer coefficients in the variables ${\bf w}:=(w_{n+i})|_{i\in{\mathbb{Z}}}$.

Definition 1**.**

An element R_{n}\in{\rm Mat}\bigl{(}2,{\mathbb{Z}}[{\bf w}](\hskip-1.42271pt(\lambda^{-1})\hskip-1.42271pt)\bigr{)} is called a matrix resolvent (MR) of ${\mathcal{L}}$, if

[TABLE]

Definition 2**.**

The basic (matrix) resolvent $R_{n}$ is defined as the MR of ${\mathcal{L}}$ satisfying

[TABLE]

The basic resolvent $R_{n}$ exists and is unique. See in Section 3 for the proof. Write

[TABLE]

Then Definition 2 for $R_{n}(\lambda)$ is equivalent to the following set of equations

[TABLE]

together with equation (14). These equations give recursive relations and initial values for the coefficients of $\alpha_{n},\,\beta_{n},\,\gamma_{n}$ (see (63)–(65) below), which will be called the MR recursive relations.

Lemma 2**.**

For an arbitrary solution $w_{n}({\bf s})$ to the discrete KdV hierarchy, let $R_{n}(\lambda)$ denote the basic resolvent of $\mathcal{L}$ evaluated at $w_{n}=w_{n}({\bf s})$. There exists a function $\tau_{n}^{\textsc{\tiny\rm dKdV}}({\bf s})$ satisfying

[TABLE]

Moreover, the function $\tau_{n}^{\textsc{\tiny\rm dKdV}}({\bf s})$ is uniquely determined by $w_{n}({\bf s})$ up to a factor of the form

[TABLE]

where $\alpha,\beta_{0},\beta_{1},\beta_{2},\cdots$ are arbitrary constants that are independent of $n,{\bf s}$.

We call $\tau_{n}^{\textsc{\tiny\rm dKdV}}({\bf s})$ the tau-function of the solution $w_{n}=w_{n}({\bf s})$ to the discrete KdV hierarchy.

The matrix-resolvent method then allows to compute logarithmic derivatives of $\tau_{n}^{\textsc{\tiny\rm dKdV}}({\bf s})$, which is achieved via the following proposition.

Proposition 1**.**

For any $k\geq 3$, the generating series of the $k_{th}$-order logarithmic derivatives of $\tau_{n}^{\textsc{\tiny\rm dKdV}}({\bf s})$ has the following expression:

[TABLE]

where it is understood that $\sigma_{k+1}=\sigma_{1}$.

The proof of this proposition is in Section 3.5.

1.3. The factorization formula

In [14] we gave the definition of tau-function for the Toda lattice using the MR approach. Observe that the discrete KdV hierarchy (3) is a reduction of the Toda lattice hierarchy. Therefore, for the arbitrary solution $w_{n}({\bf s})$ to the discrete KdV, we can also associate another tau-function $\tau_{n}({\bf s})$ of the solution $w_{n}({\bf s})$ obtained via the reduction (see Section 4.2 for the precise definition). In particular, this tau-function satisfies that

[TABLE]

It turns out that the $\tau_{n}({\bf s})$ factorizes into a product of two as given by the following theorem.

Theorem 1**.**

There exist constants $\alpha,\beta_{0},\beta_{1},\beta_{2},\cdots$ such that

[TABLE]

The proof of this theorem is in Section 4.

Remark*.*

Identity (25) echoes an identity between Hankel determinants. Indeed, let $d\mu(\lambda)$ be a measure with even moments on $\mathbb{R}$. Denote $\mu_{j}=\int\lambda^{j}d\mu(\lambda)$, $j\geq 0$. ($\mu_{\rm odd}=0$.) We know that

[TABLE]

If we deform the measure $d\mu(\lambda)$ to be $d\mu(\lambda;{\bf t})=e^{-\sum_{j\geq 1}t_{j-1}\lambda^{j}}d\mu(\lambda)$, then the LHS$\cdot(2\pi)^{-n}$ becomes a Toda tau-function (cf. the formula (3.9) of [10] and the references therein; cf. also [10, 34, 14]; the $(2\pi)^{-n}$ is a normalization factor for convenience that does not affect the fact that the LHS is already a Toda tau-function). If all the even Toda times are zero, then the ${\bf t}$-deformed measure remains even and the factorization (26) holds identically in ${\bf t}=(0,s_{1},0,s_{2},\cdots)$. Moreover, note that the RHS of (26) with deformation consists of two determinants which can be identified with the Hankel determinants associated with certain ${\bf s}$-deformed measures on $\mathbb{R}_{+}$, where ${\bf s}=(s_{1},s_{2},\cdots)$. Then to see (25) from (26), at least for special cases, one needs to further show that each of the two determinants is a tau-function for the discrete KdV hierarchy. The more precise statement for a special case together with the detailed proofs can be seen from the recent arXiv preprint by Massimo Gisonni, Giulio Ruzza and Tamara Grava [26] regarding Laguerre Unitary Ensemble (LUE) with the consideration of the parameters $\alpha=-1/2$ and $\alpha=1/2$, respectively in the notations of [26] (cf. also [8, 9]).

The next corollary follows from Proposition 1 and Theorem 1.

Corollary 1**.**

Fix $k\geq 2$. Let $w_{n}=w_{n}({\bf s})$ be an arbitrary solution to the discrete KdV hierarchy, and $\tau_{n}$ the tau-function reduced from the Toda lattice hierarchy of $w_{n}({\bf s})$. The following formula holds true:

[TABLE]

In practice, the two tau-functions $\tau_{n}$ and $\tau_{n}^{\textsc{\tiny\rm dKdV}}$ of some solution for the discrete KdV hierarchy may both have geometric/enumerative meanings; this is the case for the Hodge–GUE (see below).

Remark*.*

As we shall see from Section 4.2 that the above mentioned reduction does not mean that $v_{n}^{\textsc{\tiny\rm Toda}}$, $w_{n}^{\textsc{\tiny\rm Toda}}$ (see Section 4.2) are independent of the even Toda times $t_{0},t_{2}$, $\cdots$. The reduction means the $v_{n}^{\textsc{\tiny\rm Toda}}(0,t_{1},0,t_{3},\cdots)\equiv 0$; but the usage of the MR of Toda in the way of [14] would compute also the correlators containing the correspondence to $t_{0},t_{2}$, $\cdots$. The introductions of the MR of the discrete KdV hierarchy and of the operator $1+\Lambda$ are essential that surprisingly solve the problem in a simple form.

1.4. Application

We will first apply Corollary 1 to some counting problem. Then by using the Hodge–GUE correspondence [20, 15] we compute some combinations of Hodge integrals.

I. Enumeration of ribbon graphs with even valencies. Enumeration of ribbon graphs is closely related to the random matrix theory [4, 27, 30, 11, 34, 7]: e.g. to the Gaussian Unitary Ensembles (GUE) correlators; the partition function with coupling constants in a random matrix theory is often a tau-function of some integrable system. Given $k\geq 1$ and $j_{1},\dots,j_{k}\geq 1$, denote

[TABLE]

Here, $|j|=j_{1}+\dots+j_{k}$, and $\sum_{\Gamma}$ denotes summation over connected ribbon graphs $\Gamma$ with labelled half edges and unlabelled vertices of genus $g$ with $k$ vertices of valencies $2j_{1},\dots,2j_{k}$, and $\#\,{\rm Sym}\,\Gamma$ is the order of the symmetry group of $\Gamma$ generated by permuting the vertices.111The number $a_{g}(2j_{1},\dots,2j_{k})$ has the alternative expression $a_{g}(2j_{1},\dots,2j_{k})=\sum_{G}\frac{\prod_{\ell=1}^{k}(2j_{\ell})}{\#\,{\rm Sym}\,G}$, where $\sum_{G}$ denotes summation over connected ribbon graphs $G$ with unlabelled half-edges and unlabelled vertices of genus $g$ with $k$ vertices of valencies $2j_{1},\dots,2j_{k}$. The notation \bigl{\langle}{\rm tr}\,M^{2j_{1}}\,\cdots\,{\rm tr}\,M^{2j_{k}}\bigr{\rangle}_{c} is borrowed from the literature of random matrices, where it is often called a connected Gaussian Unitary Ensemble (GUE) correlator. For every $k\geq 1$, denote

[TABLE]

Definition 3**.**

Define a $2\times 2$ matrix-valued series R_{n}(\lambda)\in{\rm Mat}\bigl{(}2,\mathbb{Z}[n][[\lambda^{-1}]]\bigr{)} by

[TABLE]

with

[TABLE]

Theorem 2**.**

The following formulae hold true:

[TABLE]

where $R_{n}(\lambda)$ is defined in Definition 3, and it is understood that $\sigma_{k+1}=\sigma_{1}$.

In the above formulae

[TABLE]

is the Gauss hypergeometric function. Recall that it truncates to a polynomial if $a$ or $b$ are non-positive integers. In particular,

[TABLE]

The proof of Theorem 2 is in Section 5.

II. Combinations of certain special cubic Hodge integrals. The particular solution to the discrete KdV hierarchy considered here will be actually the same as in I. Denote by $\overline{\mathcal{M}}_{g,k}$ the Deligne–Mumford moduli space of stable algebraic curves of genus $g$ with $k$ distinct marked points, by $\mathcal{L}_{i}$ the $i_{\rm th}$ tautological line bundle on $\overline{\mathcal{M}}_{g,k}$, and $\mathbb{E}_{g,k}$ the Hodge bundle. Denote

[TABLE]

The Hodge integrals are some rational numbers defined by

[TABLE]

These numbers are zero unless the degree-dimension matching is satisfied

[TABLE]

We are particularly interested in the following special cubic Hodge integrals:

[TABLE]

where $\Lambda_{g}(z):=\sum_{j=0}^{g}\lambda_{j}\,z^{j}$ denotes the Chern polynomial of the Hodge bundle $\mathbb{E}_{g,k}$. Significance of these Hodge integrals is manifested by the Gopakumar–Mariño–Vafa conjecture [25, 33] regarding the Chern–Simons/string duality; see e.g. [37] and the references therein.

Notations. $\mathbb{Y}$ denotes the set of partitions. For a partition $\lambda$, denote by $\ell(\lambda)$ the length of $\lambda$ and by $|\lambda|$ the weight of $\lambda$. Denote $m(\lambda):=\prod_{i=1}^{\infty}m_{i}(\lambda)$ with $m_{i}(\lambda)$ being the multiplicity of $i$ in $\lambda$.

Definition 4**.**

For given $g,k\geq 0$ and an arbitrary set of integers $i_{1},\dots,i_{k}\geq 0$, define

[TABLE]

where $|i|:=i_{1}+\dots+i_{k}$, $\tau_{I}:=\tau_{i_{1}}\cdots\tau_{i_{k}}$, and $\tau_{\lambda+1}:=\tau_{\lambda_{1}+1}\cdots\tau_{\lambda_{\ell(\lambda)}+1}$.

It should be noted that according to (37), “$\sum_{\lambda\in\mathbb{Y}}$” in (38) is a finite sum.

The following lemma will be proved in Section 5.2.

Lemma 3**.**

The number $H_{g,i_{1},\dots,i_{k}}$ vanishes unless $|i|\leq 3g-3+k$.

Corollary 2**.**

The numbers $H_{g,i_{1},\dots,i_{k}}$ satisfy

i) For $k=0$,

[TABLE]

ii) For $k=1$, $\forall\,j\geq 1$,

[TABLE]

where $A_{n,j}$ and $B_{n,j}$ are defined in (32)–(33).

iii) For $k\geq 2$,

[TABLE]

where $R_{n}(\lambda)$ is defined as in (31).

The proof, using the Hodge–GUE correspondence and Theorem 2, will be given in Section 5.2. We note that the sum $\sum_{i_{1},\dots,i_{k}\geq 0\atop|i|\leq 3g-3+k}$ appearing in the LHS of (42), (43) has the following alternative expression, which can be deduced from Appendix A:

[TABLE]

Organization of the paper. In Section 2, we derive several useful formulae. In Section 3, we study MR, and use it to describe the discrete KdV flows and the tau-structure. Section 4 is devoted to the proof of Theorem 1. Proofs of Theorem 2 and Corollary 2 are in Section 5.

Acknowledgements. We would like to thank the anonymous referee for valuable suggestions and constructive comments that improve a lot the presentation of the paper. One of the authors D.Y. is grateful to Youjin Zhang and Don Zagier for their advising, and to Giulio Ruzza for helpful discussions. Part of the work of D.Y. was done when he was a post-doc at MPIM, Bonn; he thanks MPIM for excellent working conditions and financial supports.

2. Basic formulation

In this section we will do some preparations for the later sections by reviewing the basics of the theory of the discrete KdV hierarchy.

2.1. Some useful identities

Recall that $P(n)\;:=\;\Lambda+w_{n}\,\Lambda^{-1}$, $L=P^{2}$. Denote

[TABLE]

where the coefficients $A_{\ell,k}(n)$ and $m_{j,k}(n)$, $k\in{\mathbb{Z}}$ belong to ${\mathbb{Z}}[{\bf w}]$. It is easy to see that if $k$ is odd, or if $|k|>2j$, then $m_{j,k}\equiv 0$. It is also easy to see that

[TABLE]

Lemma 4**.**

The following identities hold true

[TABLE]

**Proof  **Comparing the constant terms of the identity

[TABLE]

we obtain that

[TABLE]

This proves (47)–(48). Similarly, comparing the coefficients of $\Lambda^{-2}$ of (50) we obtain

[TABLE]

which implies identity (49). The lemma is proved.

Lemma 5**.**

The following identities hold true

[TABLE]

**Proof  **Identities (51)–(53) are contained in the Lemma 2.2.1 of [14] (see the proof therein). Identity (54) follows from comparing coefficients of $\Lambda$ on the both sides of the following identity:

[TABLE]

The lemma is proved.

Taking $\ell=2j-1$ in identity (54) and using (46) we obtain

[TABLE]

We call this identity the key identity. It should be noted that the above identities (46)–(49), (51)–(54) hold in ${\mathbb{Z}}[{\bf w}]$ absolutely (namely, the validity does not require that $w_{n}$ is a solution of the discrete KdV hierarchy), because they are nothing but properties of the operators $P$ and $L$.

2.2. Proof of Lemma 1

Note that this lemma means the following: if $w_{n}=w_{n}({\bf s})$ satisfies (3), then it satisfies (6); vice versa. Firstly, let $w_{n}=w_{n}({\bf s})$ be an arbitrary solution to (3), i.e.,

[TABLE]

for all $j\geq 1$. Since $L=P^{2}$ we have

[TABLE]

Secondly, let $w_{n}=w_{n}({\bf s})$ be an arbitrary solution to (6), namely, it satisfies that

[TABLE]

Identity (56) implies that

[TABLE]

where we have used identity (55). Identity (57) implies that

[TABLE]

Combining the above two identities and assuming that $w_{n}\not\equiv w_{n+1}$ yields

[TABLE]

(One can see from (56) that solutions satisfying $w_{n}\equiv w_{n+1}$ are independent of ${\bf s}$. Therefore these trivial solutions also satisfy (3).) The proposition is proved.

2.3. Lax pairs in matrix form

In this subsection we write the scalar Lax pairs (9)–(10) into matrix form. The following lemma plays an important role.

Lemma 6**.**

The wave function $\psi_{n}$ satisfies that

[TABLE]

**Proof  **We have for any $j\geq 1$

[TABLE]

In the above derivations it is understood that $L=L(n)$ and $m_{j,k}=m_{j,k}(n)$. Therefore,

[TABLE]

The lemma is proved.

Lemma 7**.**

The vector-valued wave function \Psi_{n}=\bigl{(}\psi_{n},\psi_{n-2}\bigr{)}^{T} satisfies that

[TABLE]

where $V_{j}(n)$ are the following $2\times 2$ matrices

[TABLE]

**Proof  **Equation (60) follows straightforwardly from (59) and (9).

We therefore arrive at

Proposition 2**.**

The discrete KdV hierarchy are the compatibility conditions of (12) and (60):

[TABLE]

3. Tau-structure for the discrete KdV hierarchy

In this section, we use the MR method to study the tau-structure of the discrete KdV hierarchy; in particular, we will prove Proposition 1. The notations about the matrix-resolvents are the same as in the Introduction.

3.1. The MR recursive relations

Write

[TABLE]

Then we find that $a_{n,j},\,c_{n,j}$ satisfy

[TABLE]

as well as

[TABLE]

Lemma 8**.**

The basic resolvent of ${\mathcal{L}}$ exists and is unique.

**Proof  **Observe that multiplying (18) and (19) gives (20). This proves existence of $R_{n}$. Uniqueness follows directly from the MR recursive relations (63)–(65), as we can solve $a_{n,j},\,c_{n,j}$ uniquely in an algebraic way for all $j\geq 1$. The lemma is proved.

For the reader’s convenience we give in below the first few terms of the basic resolvent of $\mathcal{L}$:

[TABLE]

3.2. MR and the discrete KdV flows

In this subsection we use the basic MR to express the discrete KdV flows. (We would like to mention that the materials that we give in this subsection are rather standard.) Let $R_{n}$ be the basic matrix resolvent of $\mathcal{L}$.

Lemma 9**.**

The following formulae hold true:

[TABLE]

**Proof  **By identifying their recursive relations as well as the initial values of the recursions.

It follows from the above Lemma 9 that the matrices $V_{j}(n)$ defined in (61) have the following expressions:

[TABLE]

where $``+"$ means taking the polynomial part in $\lambda$ (including the constant term).

3.3. Loop operator

Introduce a linear operator $\nabla(\lambda)$ by

[TABLE]

It readily follows from equation (69) that

[TABLE]

where

[TABLE]

Lemma 10**.**

The following formula holds true:

[TABLE]

3.4. From MR to tau-function

The MR allows us to define tau-function of an arbitrary solution of the discrete KdV hierarchy. Recall that a family of elements $\Omega_{p;q}(n)\in{\mathbb{Z}}[{\bf w}]$, $p,q\geq 1$ are called a tau-structure of the discrete KdV hierarchy if

[TABLE]

and for an arbitrary solution $w_{n}=w_{n}({\bf s})$ of the discrete KdV hierarchy

[TABLE]

Definition 5**.**

Define $\Omega_{i;j}(n)$, $i,j\geq 1$ via the generating series

[TABLE]

Lemma 11**.**

The $\Omega_{i;j}(n)$, $i,j\geq 1$ (74) are well-defined, and live in ${\mathbb{Z}}[{\bf w}]$. Moreover, they form a tau-structure of the discrete KdV hierarchy.

**Proof  **The proof is almost identical with the one for the Toda lattice hierarchy [14] (or the one for the Drinfeld–Sokolov hierarchies [3]); details are omitted here.

**Proof  **of Lemma 2.  By Lemma 11, it suffices to prove the compatibility between (21)–(23).

Firstly, on one hand,

[TABLE]

On the other hand,

[TABLE]

Hence by using (18) we find that

[TABLE]

This proves the compatibility between (21) and (22).

Secondly, on one hand,

[TABLE]

On the other hand,

[TABLE]

Using (75) we find

[TABLE]

This proves compatibility between (21) and (23). Thirdly, the following identity

[TABLE]

shows the compatibility between (22) and (23). The proposition is proved.

3.5. Generating series of multi-point correlations functions

For an arbitrary solution $w_{n}=w_{n}({\bf s})$ to the discrete KdV hierarchy, let $\tau_{n}^{\textsc{\tiny\rm dKdV}}=\tau_{n}^{\textsc{\tiny\rm dKdV}}({\bf s})$ denote the tau-function of this solution. The logarithmic derivatives

[TABLE]

can be called the $k$-point correlation functions222We can say in a more accurate sense that the logarithmic derivatives are identified with the correlation functions, where the latter are defined as abstract differential polynomials; see for example [18] for the details. of the solution $w_{n}=w_{n}({\bf s})$.

**Proof  **of Proposition 1.  The proof can be achieved by the mathematical induction, as in [1]; we hence omit the details.

We see from Proposition 1 that the logarithmic derivatives $\frac{\partial^{k}\log\tau_{n}^{\textsc{\tiny\rm dKdV}}({\bf s})}{\partial s_{j_{1}}\dots\partial s_{j_{k}}}$ with $k\geq 2$ all live in ${\mathbb{Z}}[{\bf w}]$, as their generating series are expressed by MR via algebraic manipulations; this simple fact agrees with footnote 2 (and can be of course deduced from other techniques).

4. Proof of Theorem 1

The goal of this section is to prove Theorem 1.

4.1. Review of the MR approach to the Toda lattice hierarchy

Denote

[TABLE]

The Toda lattice hierarchy is defined as the following system of commuting flows

[TABLE]

Let us briefly review part of the results of [14]. Introduce \mathcal{U}_{n}\;=\;\biggl{(}\begin{array}[]{cc}v_{n}^{\textsc{\tiny\rm Toda}}-\lambda&w_{n}^{\textsc{\tiny\rm Toda}}\\ -1&0\\ \end{array}\biggr{)}. The basic resolvent ${\mathcal{R}}_{n}$ associated to $\mathcal{P}^{\rm M}:=\Lambda+\mathcal{U}_{n}$ is defined as the unique solution in {\rm Mat}\bigl{(}2,{\mathbb{Z}}\bigl{[}{\bf v}^{\textsc{\tiny\rm Toda}},{\bf w}^{\textsc{\tiny\rm Toda}}\bigr{]}[[\lambda^{-1}]]\bigr{)} to the problem:

[TABLE]

Write

[TABLE]

Then $\mathcal{A}_{n},\mathcal{B}_{n},\mathcal{G}_{n}$ satisfy that

[TABLE]

The following lemma was proven in [14].

Lemma 12** ([14]).**

For an arbitrary solution $v_{n}^{\textsc{\tiny\rm Toda}}=v_{n}^{\textsc{\tiny\rm Toda}}({\bf t})$, $w_{n}^{\textsc{\tiny\rm Toda}}=w_{n}^{\textsc{\tiny\rm Toda}}({\bf t})$ to the Toda lattice hierarchy there exists a function $\tau_{n}^{\textsc{\tiny\rm Toda}}({\bf t})$ such that

[TABLE]

The function $\tau_{n}^{\textsc{\tiny\rm Toda}}({\bf t})$ is uniquely determined by the solution $v_{n}^{\textsc{\tiny\rm Toda}}({\bf t})$, $w_{n}^{\textsc{\tiny\rm Toda}}({\bf t})$ up to

[TABLE]

for some constants $a_{0}$, $a_{1}$, $b_{0}$, $b_{1}$, $b_{2}$, $\dots$..

In [14] the $\tau_{n}^{\textsc{\tiny\rm Toda}}({\bf t})$ is called the tau-function of the solution $v_{n}^{\textsc{\tiny\rm Toda}}({\bf t})$, $w_{n}^{\textsc{\tiny\rm Toda}}({\bf t})$ to the Toda lattice hierarchy. The logarithmic derivatives of $\tau_{n}^{\textsc{\tiny\rm Toda}}({\bf t})$

[TABLE]

can be called $k$-point correlations functions (cf. footnote 2) of the Toda lattice hierarchy. Define

[TABLE]

4.2. Reduction to the discrete KdV hierarchy

Now consider solutions to the Toda lattice hierarchy in the ring $\mathbb{C}[[t_{0},t_{1},\dots]]\otimes{\mathcal{V}}$, where ${\mathcal{V}}$ is any ring of functions of $n$, closed under $\Lambda$ and $\Lambda^{-1}$. These solutions can be specified by (i.e. are in 1-1 correspondence to) the initial value:

[TABLE]

Let us explain how a subset of solutions to the Toda lattice hierarchy be reduced to solutions of the discrete KdV hierarchy. On one hand, let $v_{n}^{\textsc{\tiny\rm Toda}}=v_{n}^{\textsc{\tiny\rm Toda}}({\bf t})$, $w_{n}^{\textsc{\tiny\rm Toda}}=w_{n}^{\textsc{\tiny\rm Toda}}({\bf t})$ be an arbitrary solution in $\mathbb{C}[[t_{0},t_{1},\dots]]\otimes{\mathcal{V}}$ of the Toda lattice hierarchy satisfying the following type of initial conditions

[TABLE]

It follows that

[TABLE]

This further implies that the commuting flows \frac{\partial w_{n}^{\textsc{\tiny\rm Toda}}({\bf t})}{\partial t_{2j-1}}\big{|}_{t_{0}=t_{2}=t_{4}=\dots=0}\,(j\geq 1) are decoupled, namely, there are no $v_{n}^{\textsc{\tiny\rm Toda}}$-dependence in these flows (of course when restricting to $t_{0}=t_{2}=t_{4}=\dots=0$). Moreover, these flows coincide with the discrete KdV hierarchy (3). Therefore if we define

[TABLE]

then $w_{n}=w_{n}({\bf s})$ is a solution to the discrete KdV hierarchy. On the other hand, let $w_{n}=w_{n}({\bf s})$ be an arbitrary solution to the discrete KdV hierarchy in the ring $\mathbb{C}[s_{1},s_{2},\dots]\otimes{\mathcal{V}}$. Let $g(n)$ denote its initial value, i.e. $g(n):=w_{n}({\bf s}={\bf 0})$. Define $v_{n}^{\textsc{\tiny\rm Toda}}({\bf t})\,,w_{n}^{\textsc{\tiny\rm Toda}}({\bf t})$ as the unique solution in $\mathbb{C}[[t_{0},t_{1},\dots]]\otimes{\mathcal{V}}$ to the Toda lattice hierarchy with $(f(n)\equiv 0,g(n))$ as the initial value. Then $w_{n}^{\textsc{\tiny\rm Toda}}({\bf t})|_{t_{2i-1}=s_{i},\,t_{2i-2}=0,\,i\geq 1}=w_{n}({\bf s})$.

Hence the correspondence between solutions of the discrete KdV hierarchy and a suitable subset of solutions of the Toda lattice hierarchy has been established.

For a solution \bigl{(}v_{n}^{\textsc{\tiny\rm Toda}}({\bf t}),w_{n}^{\textsc{\tiny\rm Toda}}({\bf t})\bigr{)} to the Toda lattice hierarchy satisfying $v_{n}^{\textsc{\tiny\rm Toda}}({\bf 0})\equiv 0~{}(\forall\,n)$, let $\tau_{n}^{\textsc{\tiny\rm Toda}}({\bf t})$ denote the tau-function of this solution. Define $w_{n}({\bf s})$ as in (91), and

[TABLE]

Then we know that the function $w_{n}=w_{n}({\bf s})$ satisfies the discrete KdV hierarchy (3), and that

[TABLE]

As indicated above, all solutions of the discrete KdV hierarchy can be obtained from this way.

Definition 6**.**

We call $\tau_{n}({\bf s})$ the tau-function reduced from the Toda lattice hierarchy of the solution $w_{n}=w_{n}({\bf s})$ to the discrete KdV hierarchy.

Introduce the notations:

[TABLE]

Clearly, $A_{n},~{}B_{n},~{}G_{n}$ belong to ${\mathbb{Z}}[{\bf w}][[\lambda^{-1}]]$. Note that definitions of $A_{n}(\lambda),~{}B_{n}(\lambda),~{}G_{n}(\lambda)$ are in the absolute sense, namely, they do not depend on whether $w_{n}$ is a solution or not.

Lemma 13**.**

The $A_{n}(\lambda)$ satisfies

[TABLE]

**Proof  **Following from (84) and (85) with $v_{n}^{\textsc{\tiny\rm Toda}}\equiv 0$.

4.3. Proof of Theorem 1

Firstly, on one hand, it follows from the Lemma 1.2.3 of [14] that

[TABLE]

On the other hand, from (22) and (67) we find

[TABLE]

Comparing the above two expressions we find

[TABLE]

where $S(n)$ is some function depending only on $n$. Equation (96) implies that

[TABLE]

where $\widetilde{S}(n)$ is some function depending only on $n$, and $f({\bf s})$ is some function depending only on ${\bf s}$.

Secondly, it follows from (23) and (92) that

[TABLE]

Substituting (97) in (98) we find that $\widetilde{S}(n)$ can only be an affine function of $n$, namely,

[TABLE]

where $\alpha$, $\alpha^{\prime}$ are some constants independent of $n$, ${\bf s}$.

Thirdly, on one hand, using (21) we find

[TABLE]

Therefore,

[TABLE]

So

[TABLE]

On the other hand, for $P(n)=\Lambda+w_{n}\Lambda^{-1}$, recall the notation

[TABLE]

Using Lemma 12 we have

[TABLE]

where

[TABLE]

Taking

[TABLE]

we have

[TABLE]

It follows from (84), (100), and respectively (18), that

[TABLE]

Lemma 14**.**

The following identities hold true:

[TABLE]

**Proof  **Identities (103), (104) are easy consequences of (101), (100).

Note that identity (55) implies that

[TABLE]

Applying Lemma 13 in this identity yields

[TABLE]

The lemma is proved.

Observe that $C_{2}(\lambda,\mu;n,{\bf s})$ satisfies the parity symmetries

[TABLE]

So

[TABLE]

Lemma 15**.**

The following identity hold true:

[TABLE]

**Proof  **Applying (103)–(105) and the parity symmetry

[TABLE]

we find that it suffices to prove the following equality

[TABLE]

Noting that

[TABLE]

and that

[TABLE]

we find

[TABLE]

where Lemma 13 is used. The lemma is proved.

End of proof of Theorem 1. It follows from Lemma 15 that

[TABLE]

Combining with (99) we find that

[TABLE]

where $\beta_{0},\beta_{1},\beta_{2},\cdots$ are constants (independent of $n$). The theorem is proved.

5. Proofs of Theorem 2 and Corollary 2

In this section, using Proposition 1, Corollary 1 and Theorem 1, we are going to prove Theorem 2 and Corollary 2.

5.1. Ribbon graphs with even valencies

In this subsection, we first prove Theorem 2, then we give a further study to the modified GUE partition function with even couplings.

**Proof  **of Theorem 2.  Define ${\mathcal{F}}_{n}({\bf s})$ and $Z_{n}({\bf s})$ by

[TABLE]

Here $B_{m}$ denotes the $m^{\rm th}$ Bernoulli number. Then $Z_{n}({\bf s})$ is a particular tau-function (of the discrete KdV hierarchy) reduced from the Toda lattice hierarchy [14]. The initial value of $w_{n}({\bf s}):=\frac{Z_{n+1}({\bf s})\,Z_{n-1}({\bf s})}{Z_{n}({\bf s})^{2}}$ is given by $w_{n}({\bf s}={\bf 0})=n$. The theorem then follows from Lemma 14, Corollary 1, as well as the Theorem 1.1.1 of [14].

Define a formal series $Z(x,{\bf s};\epsilon)$ by

[TABLE]

Here, $x$ is the t’Hooft coupling constant [28, 29]. Recall that we could view $Z(x,{\bf s};\epsilon)$ as a tau-function reduced from the Toda lattice of the discrete KdV hierarchy under the identification $n=x/\epsilon$ as well as the flow rescalings $\partial_{s_{j}}\mapsto\epsilon\,\partial_{s_{j}}$. More precisely, define

[TABLE]

then $w(x,{\bf s};\epsilon)$ is a particular solution to the discrete KdV hierarchy:

[TABLE]

with $L:=\Lambda^{2}\,+\,w(x+\epsilon)+w(x)\,+\,w(x)\,w(x-\epsilon)\,\Lambda^{-2}$, $A_{2j-1}:=L^{j}$, $\Lambda:=e^{\epsilon\partial_{x}}$. Validity of these statements can be found in the Appendix of [14]. The initial data of this solution is given by

[TABLE]

Let $Z^{\textsc{\tiny\rm dKdV}}(x,{\bf s};\epsilon)$ be the tau-function of the solution $w(x,{\bf s};\epsilon)$. The following corollary follows from Theorem 1.

Corollary 3**.**

There exist constants $\alpha,\beta_{0},\beta_{1},\beta_{2},\cdots$ such that

[TABLE]

Note that the constants $\alpha,\beta_{0},\beta_{1},\beta_{2},\cdots$ right above now can depend on $\epsilon$. In what follows, we fix the ambiguities simply by requiring $Z^{\textsc{\tiny\rm dKdV}}(x,{\bf s};\epsilon)$ to be the unique function satisfying

[TABLE]

Remark*.*

The following formal series of ${\bf s}$

[TABLE]

was introduced in [20] by Si-Qi Liu, Youjin Zhang and the authors of the present paper, called the modified GUE partition function with even couplings, which plays an important role in a proof of the Hodge–GUE correspondence [20]. Moreover, Liu, Zhang and the authors derived the Dubrovin–Zhang loop equation for $\log\widetilde{Z}$ from the corresponding Virasoro constraints, which also provides an algorithm for computing the modified GUE correlators of an arbitrary genus [20]. Very recently, Jian Zhou [42] derived the topological recursion of Chekhov–Eynard–Orantin type for the modified GUE correlators from the Virasoro constraints constructed in [20]; moreover, as a consequence of the topological recursion, an interesting formula between intersection numbers and $k$-point functions of modified GUE correlators was obtained by Zhou [42] (see the Theorem 3 in [42] for the details); it remains an open question to match the formula of Zhou with another interesting formula obtained by Gaëtan Borot and Elba Garcia-Failde [5] (see the Corollary 12.3 of [5]) as a consequence of the Hodge–GUE correspondence (or with a slightly different but equivalent consequence like (124) in below), which may lead to a new proof of the Hodge–GUE correspondence. Last but not least, as a corollary of Theorem 2, let us give a third algorithm of computing the modified GUE correlators based on the following full genera formulae:

[TABLE]

where $\langle\phi_{j_{1}}\cdots\phi_{j_{k}}\rangle(x;\epsilon)$ denote the modified GUE correlators with even couplings, defined by

[TABLE]

and $R_{n}(\lambda)$ is defined in Definition 3. We notice that the reason that one can talk about “genus” for $\log\widetilde{Z}$ is because $\log\widetilde{Z}$ is even in $\epsilon$ and so are $\langle\phi_{j_{1}}\cdots\phi_{j_{k}}\rangle(x;\epsilon)$. A concrete algorithm using the formulae of the form (114)–(115) for calculating the corresponding correlators including certain large genus asymptotics is given in [16].

Remark*.*

Very recently it was shown in [26] that $Z^{\textsc{\tiny\rm dKdV}}(x,{\bf s};\epsilon)$ and $Z^{\textsc{\tiny\rm dKdV}}(x+\epsilon,{\bf s};\epsilon)$ are identified with the LUE partition functions with $\alpha=-1/2$ and $\alpha=1/2$, respectively. One can obtain their $k$-point series by putting $x\rightarrow x\mp\frac{\epsilon}{2}$ in (114)–(115). An interesting genus expansion for $Z^{\textsc{\tiny\rm dKdV}}(x,{\bf s};\epsilon)$ was discovered in [9]. The interplay between $Z^{\textsc{\tiny\rm dKdV}}$ and $\widetilde{Z}$ suggests a Hurwitz/Hodge correspondence that deserves a further study.

Using the definitions of $Z(x,{\bf s};\epsilon)$ and $\widetilde{Z}(x,{\bf s};\epsilon)$ and using the expansion

[TABLE]

with $E_{k}$, $k\geq 0$ being the Euler numbers, we have the following formula:

[TABLE]

In other words, the modified GUE correlators with even couplings (116) have the expressions:

[TABLE]

where $k\geq 1$ and $j_{1},\dots,j_{k}\geq 1$. It should be noted that the $\langle\phi_{j_{1}}\cdots\phi_{j_{k}}\rangle(x;\epsilon)$ with $k\geq 1$, $j_{1},\dots,j_{k}\geq 1$ is a polynomial of $x$.

5.2. Combinations of certain special cubic Hodge integrals

Based on the Hodge–GUE correspondence and using Theorem 2, we compute in this subsection combinations of certain special cubic Hodge integrals. More precisely, we will prove Corollary 2.

The cubic Hodge free energy associated with \Lambda_{g}(-1)\,\Lambda_{g}(-1)\,\Lambda_{g}\bigl{(}\tfrac{1}{2}\bigr{)} is defined by

[TABLE]

Here, ${\bf t}=(t_{0},t_{1},\dots)$. (Warning: Avoid from confusing with the variables $t_{\ell}$, $\ell\geq 0$ of the Toda lattice hierarchy used in Section 4.) The Hodge–GUE correspondence connects $\mathcal{H}({\bf t};\epsilon)$ with the GUE partition function with even couplings, which is given by the following theorem.

Theorem A. ([20, 15]) The following identity holds true:

[TABLE]

where $\bar{s}_{j}:=\binom{2j}{j}\,s_{j}$ and

[TABLE]

Recall from (113) that the modified GUE partition function with even couplings $\widetilde{Z}$ is defined as the unique series of $x-1$ and ${\bf s}$ satisfying

[TABLE]

Combining (120) with (118) we obtain the following corollary.

Corollary 4**.**

The following formula holds true:

[TABLE]

Denote \Omega_{g}:=\Lambda_{g}(-1)\,\Lambda_{g}(-1)\,\Lambda_{g}\bigl{(}\frac{1}{2}\bigr{)} as in the introduction, and write

[TABLE]

It might be helpful to notice that for $g=1$, $\deg\,\Omega_{1}\leq 1$; for $g\geq 2$, $\deg\,\Omega_{g}\leq 3g-3$. Motivated by Theorem A, let us consider the following combination of Hodge integrals. For any given $k\geq 0$, $i_{1},\dots,i_{k}\geq 0$, define a formal series $H_{i_{1},\dots,i_{k}}(x;\epsilon)\in\epsilon^{-2}{\mathbb{Q}}[[x-1,\epsilon^{2}]]$ by

[TABLE]

Note that in the notation $\langle\,\dots\,\rangle_{g}$ , we omit the index $m$ from $\langle\,\dots\,\rangle_{g,m}$. For such an abbreviation, $m$ should be recovered from counting the number of $\tau$’s in $``\dots"$. Therefore, for each fixed $g,d$ and for each monomial in the Taylor expansion $e^{(x-1)\tau_{0}}=\sum_{r=0}^{\infty}\frac{1}{r!}(x-1)^{r}$, the above summation over partitions $\sum_{\lambda\in\mathbb{Y}}$ is a finite sum, i.e., the degree-dimension matching $|\lambda|=3g-3+k+r-d-|i|$ has to be hold. Lemma 3 also easily follows from this constrain with $r=0$ taken. The numbers $H_{g,i_{1},\dots,i_{k}}$ defined by (38) and the formal series $H_{i_{1},\dots,i_{k}}(x;\epsilon)$ are clearly related by

[TABLE]

Proposition 3**.**

For any $k\geq 0$ and $j_{1},\dots,j_{k}\geq 1$, the following formula holds true:

[TABLE]

**Proof  **Note that the {\mathcal{H}}\bigl{(}{\bf t};\sqrt{2}\epsilon\bigr{)} has the expression

[TABLE]

Formula (124) is then proved by substituting (119) and by using Corollary 4.

**Proof  **of Corollary 2. Note that the $k=0$ case is already given in [20]. By taking $x=1$ in (124) and using Theorem 2 we find (43). Formula (42) is then implied in a standard way by using the following linear equation (proven in [20])

[TABLE]

and the fact that

[TABLE]

Here $\widetilde{Z}=\widetilde{Z}(x,{\bf s};\epsilon)$ denotes the modified GUE partition function with even couplings. Note that the fact (126) can be obtained by taking the coefficients of $\lambda^{-1}$ on the both sides of (21). The corollary is proved.

Appendix A On consequence of the Hodge–GUE correspondence

In this appendix, we derive a consequence of the Hodge–GUE correspondence that has a similar flavour to formula (124). Note that

[TABLE]

Then by comparing the coefficients of $s_{p_{1}}\dots s_{p_{l}}$ of the both sides of (121) we get

[TABLE]

Here $\langle\sigma_{p_{1}}\dots\sigma_{p_{l}}\rangle(x;\epsilon)=:\sum_{g\geq 0}\epsilon^{2g-2}\langle\sigma_{p_{1}}\dots\sigma_{p_{l}}\rangle_{g}(x)$, and $\langle\sigma_{p_{1}}\dots\sigma_{p_{l}}\rangle(x;\epsilon)$ are the modified GUE correlators with even couplings defined in (116). Taking $x=1$ we find

[TABLE]

A further consideration to (128) was given in [5].

Combining (127) with (117) we find for any fixed $l\geq 1$, $p_{1},\dots,p_{l}\geq 1$ the following identities:

[TABLE]

Note that for any $g\geq 0$, the RHS is a priori a power series of $x-1$, but the LHS shows that it is actually a monomial of $x$ and so is also a polynomial of $x-1$. This subset of the identities deserve a further investigation. Moreover, the LHS vanishes when $g$ is sufficiently large, an so is the RHS; this provides another subset of the identities for the cubic Hodge integrals.