# Topological open/closed string dualities: matrix models and wave functions

Sujay K. Ashok$^a$ and Jan Troost$^b$

E-mail: sashok@imsc.res.in, jan.troost@ens.fr

**ABSTRACT:** We sharpen the duality between open and closed topological string partition functions for topological gravity coupled to matter. The closed string partition function is a generalized Kontsevich matrix model in the large dimension limit. We integrate out off-diagonal degrees of freedom associated to one source eigenvalue, and find an open/closed topological string partition function, thus proving open/closed duality. We match the resulting open partition function to the generating function of intersection numbers on moduli spaces of Riemann surfaces with boundaries and boundary insertions. Moreover, we connect our work to the literature on a wave function of the KP integrable hierarchy and clarify the role of the extended Virasoro generators that include all time variables as well as the coupling to the open string observable.

**KEYWORDS:** D-branes, Integrable Hierarchies, Matrix Models, 2D Gravity

ARXIV ePRINT: 1907.02410

<footnote>$^a$Institute of Mathematical Sciences, Homi Bhabha National Institute (HBNI), IV Cross Road, C.I.T. Campus, Taramani, Chennai, 600113 Tamil Nadu, India</footnote>

<footnote>$^b$Département de Physique, Ecole Normale Supérieure, 24 rue Lhomond, Paris, France</footnote>
# Contents

<table><tr><td>1 Introduction</td><td>1</td></tr><tr><td>2 The open/closed string duality with matter</td><td>3</td></tr><tr><td>2.1 A brief overview</td><td>3</td></tr><tr><td>2.2 The generalized Kontsevich model extended</td><td>3</td></tr><tr><td>2.3 The closed and open/closed string partition functions</td><td>6</td></tr><tr><td>2.3.1 The closed string partition function</td><td>6</td></tr><tr><td>2.3.2 The wave potential and the open string partition function</td><td>7</td></tr><tr><td>3 The open/closed string Virasoro constraints</td><td>10</td></tr><tr><td>3.1 The open/closed Virasoro algebra</td><td>10</td></tr><tr><td>3.2 From closed to open Virasoro</td><td>13</td></tr><tr><td>3.3 The integrable hierarchy and the relation to geometry</td><td>14</td></tr><tr><td>4 Conclusions</td><td>16</td></tr><tr><td>A Properties of the Fourier transform</td><td>17</td></tr><tr><td>B An equivalent Virasoro algebra</td><td>18</td></tr><tr><td>C Another pure duality</td><td>20</td></tr></table>

## 1 Introduction

Two-dimensional gravity is an important tool in the study of quantum gravity. It is a simpler analogue of four-dimensional gravity, and it arises in the diffeomorphism invariant theory living on a string theory world sheet. Two-dimensional quantum gravity has been solved to a large extent, in three ways: via double-scaled matrix models [1–3], using conformal field theory techniques [4–6] and in its topological formulation [7–9]. The three approaches were proven to be equivalent in many instances. See e.g. [10, 11] for reviews of the two-dimensional theory of quantum gravity.

While the topological formulation of two-dimensional gravity on Riemann surfaces without boundary was put on a firm footing a while back [7, 8], the theory on Riemann surfaces with boundary was only recently understood rigorously [12]. This gave rise to a flurry of mathematical activity which made earlier observations in the physics literature [13–17] precise [18–21]. See [22] for a partial review peppered with physical insight and [23–25] for more recent results.
Of particular interest to us here is the open/closed duality which was understood fairly well in the string theory literature [13–17], and which was rigorously derived, including important additional details, in the more recent mathematical physics literature [18–21] in the case of pure two-dimensional gravity. This duality relates an open/closed string partition function to a purely closed string partition function; to be precise, the addition of a D-brane in topological string theory is transmuted into a shift of the background, a renormalization of the partition function, an operator insertion and an integral transform. While the conceptual framework goes back to [15–17], precise formulas were given more recently in [19]. Only the latter allow to make direct contact with the rigorous analysis in algebraic geometry [12, 19–21].

The main mathematical tool in the derivation of the precise form of open/closed string duality [19] was the Kontsevich matrix model of pure two-dimensional gravity [9]. This matrix model depends only on the closed string sources, and can be viewed as a closed string matrix model in that sense. The closed versus open/closed duality in the case of pure gravity was derived in a particularly clear manner by integrating out off-diagonal degrees of freedom of a $(N+1) \times (N+1)$ matrix model to obtain a $N \times N$ matrix model depending on one extra eigenvalue which is the integral transform of an open string coupling (or D-brane modulus) [19]. In the large $N$ limit, this then gives rise to a duality between a closed string theory and an open/closed string theory. See also [18, 26] for related results.

In this paper, we find an alternative to the original mathematical derivation that allows for a generalization to the case of two-dimensional topological gravity coupled to matter. The open/closed duality in this case was understood to a degree in [27]. We provide the mathematically precise relation between the open/closed and closed string partition function using the technique described above of integrating out off-diagonal degrees of freedom at finite $N$, and then taking the infinite $N$ limit. This will lead to considerable additional insight. Indeed, we naturally find a determinant insertion in the Kontsevich matrix model. We moreover find that the closed string integrable hierarchy is intuitively extended to include the missing times. Importantly, by matching to the integrable system literature, we are able to obtain a precise algebra-geometric understanding of the resulting open/closed and extended partition functions. We also improve our understanding of the Virasoro algebras that govern the generating functions.

The paper is structured as follows. In section 2, we first present the result of partially integrating out off-diagonal matrix elements in the Kontsevich matrix model in an elementary fashion. We then interpret the resulting equations in terms of concepts in the integrable systems literature, as well as in string theory. We moreover build the bridge to the open/closed topological intersection numbers studied in algebraic geometry. In section 3 we provide a guide to the various Virasoro algebras that appear in the context of these integrable systems, and clarify how to interpret them. The knowledge gained is used to prove a crucial statement in section 2. In section 4, we conclude with a summary and some open problems. In appendix A we provide properties of a generalized Fourier transform, in appendix B the result of conjugating Virasoro generators to a more familiar form, and in appendix C, we review an alternative viewpoint on open/closed duality in the case of pure gravity.
## 2 The open/closed string duality with matter

In this section, we first derive a central technical result in an elementary manner. We start from the generalized Kontsevich matrix model [9, 28–30] with matrix integration variables of dimension $(N+1) \times (N+1)$ and integrate out the off-diagonal matrix elements compared to $N \times N$ and $1 \times 1$ blocks around the diagonal. As we will discuss, the resulting equation can be interpreted in terms of standard concepts in integrable systems, and has interesting conceptual consequences in topological gravity and low-dimensional string theory.

### 2.1 A brief overview

The generalized Kontsevich model is a matrix generalization of the higher Airy function [9]. When the matrix model integration variable tends to infinite size, it becomes a generating function for the p-spin intersection numbers on moduli spaces of Riemann surfaces [9, 28–33]. Equivalently, it is the partition function of topological gravity coupled to matter of type p, or a closed string theory of gravity plus matter in dimension smaller than one [34–36].

Our starting point is the study of the generalized Kontsevich model at finite $N$. In particular, we will perform the integration over $2N$ off-diagonal degrees of freedom in the $(N+1) \times (N+1)$ matrix integration variables, and will be left with a $N \times N$ matrix model and an effective action for one more diagonal degree of freedom. In the large $N$ limit, both the $(N+1) \times (N+1)$ and the $N \times N$ dimensional model will correspond to closed string matrix models. The extra eigenvalue will eventually be related to a D-brane or open string modulus. The equation that results from the integration will thus allow for an interpretation in terms of open/closed string duality. Equivalently, it allows for an interpretation of intersection numbers of Riemann surfaces with boundaries and bulk and boundary insertions, in terms of intersection numbers of Riemann surfaces with only bulk insertions.

Our approach has been inspired by a number of sources. Firstly, the idea of open/closed string duality can at least be traced back to the advent of D-branes in string theory [37]. The fact that it takes a particularly simple form in the case of pure topological gravity was understood in [16]. These ideas were extended to the case of topological gravity with matter in [27]. At a technical level, these references have a different approach. Secondly, the more rigorous reference [19] follows almost the same technical route we described above to render the physical intuition in [16] mathematically precise, but there are important differences that we will highlight.

### 2.2 The generalized Kontsevich model extended

A topological closed string partition function $\tau(\Lambda)$ corresponding to topological gravity in two dimensions coupled to topological matter is captured by the generalized Kontsevich model [9, 28–30]:

$$ \tau(\Lambda) = \frac{N(\Lambda)}{D(\Lambda)}, \qquad (2.1) $$
where the numerator $N(\Lambda)$ and the denominator $D(\Lambda)$ read

$$ N(\Lambda) = \int [dM]_N e^{-\alpha \frac{1}{p+1} \mathrm{Tr} [(M+\Lambda)^{p+1}]_{\ge 2}} \\ D(\Lambda) = \int [dM]_N e^{-\frac{\alpha}{2} \sum_{k=0}^{p-1} \mathrm{Tr} M \Lambda^k M \Lambda^{p-1-k}} \qquad (2.2) $$

The square brackets with lower index $\ge 2$ indicate that one should consider only terms that are of order two or higher in the $N \times N$ hermitian matrix integration variable $M$. The coupling constant $\alpha$ can be identified with one over the string coupling, $\alpha = -1/g_s$. The partition function of the closed string is a function of the matrix source $\Lambda$ which codes the values of all couplings to the matter primaries and gravitational descendants. When the couplings go to zero, the matrix source $\Lambda$ goes to infinity, and the $\tau$ function approaches one. The matter content is specified by the order $p+1$ of the potential term. We have that $p=2$ for pure gravity.<fnref n="1" />

Our set-up is elementary. We take the source matrix as well as the integration variable to be an $(N+1) \times (N+1)$ matrix and wish to interpret the last diagonal entry of the source matrix as (the dual of) an open string (or boundary insertion) modulus. The matrix integral only depends on the eigenvalues of the source matrix $\Lambda_z$, which we take to be diagonal:

$$ \Lambda_z = \mathrm{diag}(\lambda_1, \dots, \lambda_N, z). \qquad (2.3) $$

We used the special notation $z$ for the last diagonal entry, which we plan to single out. We have the corresponding closed string partition function:

$$ \tau(\Lambda_z) = \frac{N(\Lambda_z)}{D(\Lambda_z)}. \qquad (2.4) $$

To perform the integration over the $2N$ off-diagonal degrees of freedom, we integrate over all $N(N+1)$ off-diagonal degrees of freedom, and then reinstate those that we wish to keep. To perform the integration, our main tool is the Harish-Chandra-Itzykson-Zuber integration formula [38, 39]. We integrate in two steps. Firstly, we concentrate on the numerator, and then we perform the Gaussian integration in the denominator.

The exponent in the numerator (2.2) of equation (2.4) can be simplified by shifting away the source $\Lambda_z$, to keep only a term which has power $p+1$, and a linear and constant term in the integration variable $M$. Indeed, we have

$$ [(M + \Lambda_z)^{p+1}]_{\ge 2} = (M + \Lambda_z)^{p+1} - \Lambda_z^{p+1} - (p+1)\Lambda_z^p M, \qquad (2.5) $$

and therefore upon shifting $M \to M - \Lambda_z$, we find that the numerator $N(\Lambda_z)$ takes the form:

$$ N(\Lambda_z) = \int [dM]_{N+1} e^{-\alpha \mathrm{Tr} \left( \frac{M^{p+1}}{p+1} - M \Lambda_z^p + \frac{p}{p+1} \Lambda_z^{p+1} \right)}. \qquad (2.6) $$

<footnote n="1">¹We remark that these matrix integrals should be thought off as consisting of a Gaussian integration plus an exponential of higher powers that are to be expanded as formal power series in the integration variable. For a careful treatment of these and other aspects see e.g. [9, 19].</footnote>
We can then integrate over the unitary, angular variables that serve to diagonalise the matrix M. To that end we parameterize the $(N+1) \times (N+1)$ hermitian matrix M in terms of a unitary matrix U and diagonal matrix $M_d$ as $M = UM_dU^{-1}$. The matrix integral factorizes:

$$ N(\Lambda_z) = \frac{\pi^{\frac{N(N+1)}{2}}}{\prod_{i=1}^{N+1} i!} \int \prod_{i=1}^{N+1} dm_i \Delta_{N+1}(m) e^{-\frac{\alpha}{p+1} \left( \sum_{i=1}^{N+1} m_i^{p+1} + p \operatorname{Tr} \Lambda_z^{p+1} \right)} \\ \times \int [dU]_{N+1} e^{\alpha \operatorname{Tr} \Lambda_z^p U M_d U^{-1}}, \qquad (2.7) $$

where the integration variables $m_i$ are the eigenvalues of the diagonal matrix $M_d$ and we defined the Vandermonde determinant measure factor

$$ \Delta_N(m) = \prod_{1 \le i < j \le N} (m_j - m_i)^2. \qquad (2.8) $$

We use the Harish-Chandra-Itzykson-Zuber formula [38, 39]

$$ \int [dU]_{N+1} e^{\alpha \operatorname{Tr} AUBU^{-1}} = \frac{\prod_{i=1}^{N} i!}{\alpha^{\frac{N(N+1)}{2}}} \frac{\det(e^{\alpha a_i b_j})}{\prod_{1 \le i < j \le N+1} (a_j - a_i)(b_j - b_i)}, \qquad (2.9) $$

to perform the integral over the unitary matrix U:

$$ N(\Lambda_z) = \frac{1}{(N+1)!} \left(\frac{\pi}{\alpha}\right)^{\frac{N(N+1)}{2}} \frac{1}{\prod_{1 \le i < j \le N+1} (\lambda_{z,j}^p - \lambda_{z,i}^p)} e^{-\frac{\alpha p}{p+1} \operatorname{Tr} \Lambda_z^{p+1}} \\ \times \int \prod_{i=1}^{N+1} dm_i e^{-\frac{\alpha}{p+1} \sum_{i=1}^{N+1} m_i^{p+1}} \prod_{1 \le i < j \le N+1} (m_j - m_i) \det \left( e^{\alpha \lambda_{z,i}^p m_j} \right). \qquad (2.10) $$

In the next step, we isolate the integral over the last diagonal variable $m_{N+1}$ and write the remaining factors as an $N \times N$ matrix integral that we wish to identify as a generalized Kontsevich model with an extra insertion. The first step towards this goal is to expand the determinant along the $(N+1)$st row of the matrix:

$$ \det_{N+1} (e^{\alpha \lambda_{z,j}^p m_i}) \prod_{i<j}^{N+1} (m_j - m_i) = \sum_{l=1}^{N+1} e^{\alpha z^l m_l} \det_N (e^{\alpha \lambda_j^p m_i})_{i \neq l} \prod_{i<j \neq l}^{N+1} (m_j - m_i) \prod_{i \neq l}^{N+1} (m_l - m_i). \qquad (2.11) $$

By suitably changing variables in each term of the sum and using the permutation invariance of the measure, one can show that each resulting integral contributes equally.<fnref n="2" /> Therefore, this sum may be written as $(N+1)$ times the contribution from the $l = N+1$ term. For easier writing, we rename the diagonal integration variable $m_{N+1} = s$ and find the numerator:

$$ N(\Lambda_z) = \frac{1}{N!} \left(\frac{\pi}{\alpha}\right)^{\frac{N(N+1)}{2}} \frac{1}{\prod_{i=1}^{N} (z^p - \lambda_i^p)} e^{-\alpha \frac{p}{p+1} (z^{p+1} + \operatorname{Tr} \Lambda^{p+1})} \\ \times \int ds \, e^{-\alpha \frac{s^{p+1}}{p+1} + \alpha z^p s} \int \prod_{i=1}^{N} dm_i \Delta_N(m_i) \frac{\det(e^{\alpha \lambda_i^p m_j}) \prod_{i=1}^{N} (s - m_i)}{\prod_{1 \le i < j \le N} (m_j - m_i)(\lambda_j^p - \lambda_i^p)}. \qquad (2.12) $$

<footnote n="2">²We thank Alexandr Buryak for clarifying analogous steps that arise in [19].</footnote>
At this point, we trace our steps backwards and use the Harish-Chandra-Itzykson-Zuber formula in reverse. Thus, we write the numerator in terms of the integration over a single variable s and a $N \times N$ hermitian matrix M:

$$ N(\Lambda_z) = \left(\frac{\pi}{\alpha}\right)^N \frac{1}{\det(z^p - \Lambda^p)} \int ds \, e^{-\alpha \left(\frac{s^{p+1}}{p+1} - z^p s + \frac{p}{p+1} z^{p+1}\right)} \\ \times \int [dM]_N e^{-\alpha \mathrm{Tr}\left(\frac{M^{p+1}}{p+1} - \Lambda^p M + \frac{p}{p+1} \Lambda^{p+1}\right)} \det(s - M). \quad (2.13) $$

This finishes the first and harder part of the calculation. Secondly, we must keep track of the denominator $D(\Lambda_z)$ that serves to anchor the partition function $\tau$ at 1 for large source matrix $\Lambda$. The integration in the denominator is Gaussian, and the $N \times N$ determinant that will serve to normalize the final $N \times N$ matrix integration can easily be factored out:

$$ \begin{aligned} D(\Lambda_z) &= \int [dM]_{N+1} e^{-\frac{\alpha}{2} \sum_{k=0}^{p-1} \mathrm{Tr} M \Lambda_z^k M \Lambda_z^{p-1-k}} \\ &= \prod_{i=1}^{N+1} \sqrt{\frac{2\pi}{\alpha}} \frac{1}{\sqrt{p\lambda_{z,i}^{p-1}}} \prod_{i \neq j}^{N+1} \sqrt{\frac{\pi}{\alpha}} \frac{1}{\sqrt{\lambda_{z,j}^p - \lambda_{z,i}^p}} \\ &= \sqrt{\frac{2\pi}{\alpha}} (pz^{p-1})^{-\frac{1}{2}} \prod_{i=1}^{N} \frac{\pi}{\alpha} \frac{\lambda_i - z}{\lambda_i^p - z^p} D(\Lambda) \\ &= \sqrt{\frac{2\pi}{\alpha}} \left(\frac{\pi}{\alpha}\right)^N \frac{\det(z - \Lambda)}{\det(z^p - \Lambda^p)} (pz^{p-1})^{-\frac{1}{2}} D(\Lambda). \end{aligned} \quad (2.14) $$

The integrating-out in both numerator and denominator provides us with a final formula for the tau-function at finite N:

$$ \tau(\Lambda_z) = \sqrt{\frac{\alpha p}{2\pi}} \frac{z^{\frac{p-1}{2}}}{\det(z-\Lambda)} e^{-\alpha \frac{p}{p+1} z^{p+1}} \int ds \, e^{-\frac{\alpha}{p+1} s^{p+1} + \alpha z^p s} \\ \times \frac{1}{D(\Lambda)} \int [dM]_N e^{-\alpha \frac{1}{(p+1)} \mathrm{Tr}[(M+\Lambda)^{p+1}]}_{\ge 2} \det(s-M-\Lambda). \quad (2.15) $$

This is our first technical result. To clarify its significance, we repackage it in various ways, and then provide an interpretation.

## 2.3 The closed and open/closed string partition functions

To make contact with both the string theory and integrability literature, we slightly reshuf-fle the result (2.15):

$$ \det\left(1 - \frac{z}{\Lambda}\right) \tau(\Lambda_z) = \sqrt{\frac{\alpha p}{2\pi}} z^{\frac{p-1}{2}} e^{-\alpha \frac{p}{p+1} z^{p+1}} \int ds \, e^{\alpha z^p s} \quad (2.16) \\ \times \frac{1}{D(\Lambda)} \int [dM]_N e^{-\alpha \frac{1}{(p+1)} \mathrm{Tr}[(M+\Lambda)^{p+1}]}_{\ge 2} e^{-\alpha \frac{1}{p+1} s^{p+1}} \frac{\det(\Lambda + M - s)}{\det(\Lambda)}. $$

### 2.3.1 The closed string partition function

We wish to be more specific about the interpretation and meaning of both the left and the right hand side of the equality (2.16). Firstly, let us remind the reader that the closed
string partition function $\tau(\Lambda)$ is the $\tau$-function of a reduced KP integrable hierarchy. The times $t_n$ of the integrable hierarchy are defined in terms of the source matrix $\Lambda$ by

$$t_n = -\frac{1}{n} \mathrm{Tr} \Lambda^{-n}. \qquad (2.17)$$

The times play the role of closed string couplings in topological string theory. A well known and important result is the independence of the closed string partition function on the times $t_{np}$, where $n \in \mathbb{N}$ [9, 30]. We will make use of this repeatedly in what follows.

Secondly, we consider the closed string partition function $\tau(\Lambda_z)$ with the extra source eigenvalue $z$ (as in equation (2.3)). We can think of the addition of the extra source eigenvalue $z$ as a redefinition of the closed string couplings $t_n$. Indeed, we have:

$$\tilde{t}_n = -\frac{1}{n} \mathrm{Tr} \Lambda_z^{-n} = -\frac{1}{n} \mathrm{Tr} \Lambda^{-n} - \frac{1}{n z^n} = t_n - \frac{1}{n z^n}. \qquad (2.18)$$

Thus, if we view the $\tau$-function as a function of the time variables $t_n$, which are independent in the large $N$ limit, then we can define a shifted closed string partition function by the formula:

$$\tau(\Lambda_z) = \tau \left( t_n - \frac{1}{n z^n} \right). \qquad (2.19)$$

The extra source term redefines the closed string background. This elementary aspect of the duality formula (2.16) is well-understood in the literature.

Let us introduce an alternative symbol for the determinant on the left hand side of equation (2.16):

$$e^{\xi} = \det \left( 1 - \frac{z}{\Lambda} \right), \qquad (2.20)$$

which implies that

$$\xi = \sum_{n=1}^{\infty} z^n t_n. \qquad (2.21)$$

We then recognize on the left hand side of equation (2.16) a quantity from integrable systems [42], namely the extended partition function or wave potential $\tau_{\mathrm{ext}}(\Lambda, z)$ equal to:

$$\tau_{\mathrm{ext}}(\Lambda, z) = e^{\xi} \tau(\Lambda_z). \qquad (2.22)$$

It is important that the extended tau-function $\tau_{\mathrm{ext}}(\Lambda_z)$ (in contrast to the tau-function $\tau(\Lambda)$) depends on all times of the KP integrable hierarchy, and in particular on the times $t_{np}$. The factor $e^{\xi}$ introduces an exponential $t_{np}$ dependence. Thus the left-hand side $\tau_{\mathrm{ext}}(\Lambda_z)$ of the equality (2.16) is an extended and shifted closed string partition function, equal to the wave potential of the KP integrable hierarchy.

### 2.3.2 The wave potential and the open string partition function

Moreover, we also wish to compactly code and interpret the right hand side of equation (2.16). To that end we define an open/closed partition function:

$$\tau^{\mathrm{op+cl}}(\Lambda, s) = \frac{1}{D(\Lambda)} \int [dM] e^{-\alpha \frac{1}{(p+1)} \mathrm{Tr} [(M+\Lambda)^{p+1}]}_{\ge 2} e^{-\alpha \frac{1}{p+1} s^{p+1}} \frac{\det(\Lambda + M - s)}{\det(\Lambda)}, \qquad (2.23)$$
which is equal to the generalized Kontsevich matrix integral with a normalized determinant insertion labelled by a coupling constant $s$. Our result (2.16) then reads more compactly:

$$ \tau_{\text{ext}}(\Lambda, z) = \sqrt{\frac{\alpha p}{2\pi}} z^{\frac{p-1}{2}} e^{-\alpha \frac{p}{p+1} z^{p+1}} \int ds \, e^{\alpha z^p s} \tau^{\text{op+cl}}(\Lambda, s). \qquad (2.24) $$

If we moreover define a formal Fourier transform of a function (or rather a formal power series) $f(s)$ by:

$$ \Phi[f(s)](z) := \sqrt{\frac{\alpha p}{2\pi}} z^{\frac{p-1}{2}} \int ds \, e^{\alpha \left( \frac{z^{p+1}}{p+1} + z^p s \right)} f(s+z), \qquad (2.25) $$

we have that our duality succinctly reads:

$$ \tau_{\text{ext}}(\Lambda, z) = \Phi \left[ \tau^{\text{op+cl}}(\Lambda, s) \right] (z). \qquad (2.26) $$

We summarize that the (finite $N$) generalization of the) extended closed string partition function equals the integral transform of the open/closed string partition function. Importantly, we still need to argue that the correct interpretation of the right hand side of the result (2.26) is indeed as an open/closed string partition function. In order to do so, we need considerably more background.

Relatively recently, the following wave function was introduced in the integrable system literature [41, 42]. Given a $p$-reduced integrable KP hierarchy and the associated Lax operator $L$ — see [40] for a pedagogical introduction to the subject — satisfying the evolution equations:

$$ \frac{\partial L}{\partial t_m} = [[L^{\frac{m}{p}}]_{\ge 0}, L], \qquad (2.27) $$

one defines a wave function $\Psi(t_k)$ that satisfies the differential equations:

$$ \frac{\partial \Psi}{\partial t_m} = [L^{\frac{m}{p}}]_{\ge 0} \Psi, \qquad (2.28) $$

as well as the initial condition [41, 42]

$$ \Psi_{t_n \ge 0} = 1. \qquad (2.29) $$

The interest in the wave function $\Psi(t_k)$ originates in the advance that has been made in computing intersection numbers on moduli spaces of Riemann surfaces with boundaries [21, 41]. Indeed, in the presence of boundaries, the reduction of the integrable hierarchy no longer takes place, and the generating function of intersection numbers does depend on the times $t_{np}$. Moreover, there is an extra dependence on a coupling constant $s$ that counts open string or boundary insertions. The partition function of closed/open intersection numbers can be expressed in terms of the wave function $\Psi$ [20, 21]. Let us describe the final result of the algebro-geometric calculations. We introduce free energies that generate intersection numbers on Riemann surfaces with and without boundary:

$$ \tau = e^{F^{\text{cl}}} \\ \tau_{\text{geom}}^{\text{op+cl}} = e^{F^{\text{cl}}} + F^{\text{op}} = \tau e^{F^{\text{op}}}. \qquad (2.30) $$
Crucially, it is shown in [21] that the generator $F^{\text{open}}$ of geometric open/closed intersection numbers is given by the logarithm of the wave function $\Psi(t_k)$ after a substitution of variables:

$$F^{\text{op}} = \log \Psi(t_i, t_p \rightarrow t_p - \alpha s). \qquad (2.31)$$

It is important to note here that in the generating function, we have allowed for extended closed string amplitudes, namely, amplitudes where we have added boundaries but no explicit open string insertions. In this regard, see the useful intuitive remarks in [20].

Finally, we are ready to state the relation our matrix model has to the algebro-geometric open/closed intersection numbers. We propose that our matrix integral $\tau^{\text{op+cl}}$ is equal to its geometric counterpart:

$$\tau^{\text{op+cl}} = \tau_{\text{geom}}^{\text{op+cl}} = \tau(t) \Psi(t, t_p - \alpha s). \qquad (2.32)$$

The left hand side is our open/closed string partition function $\tau^{\text{op+cl}}$ defined in terms of the generalized Kontsevich integral with normalized determinant insertion (2.23), while the right hand side is the generator of geometric open/closed string intersection numbers, according to [21].

We will prove the identification (2.32) in the following. As a warm-up, let us first show that our $\tau^{\text{op+cl}}$ indeed only depends on the combination $t_p - \alpha s$ as expected from equation (2.32). We write the open/closed partition function $\tau^{\text{op+cl}}(\Lambda, s)$ as the inverse integral transform of the wave potential $\tau_{\text{ext}}(\Lambda, z):^3$

$$\tau^{\text{op+cl}}(\Lambda, s) = \sqrt{\frac{\alpha p}{2\pi}} \int dz \, z^{\frac{p-1}{2}} e^{\alpha \frac{p}{p+1} z^{p+1} - \alpha z^p s} \tau_{\text{ext}}(\Lambda, z). \qquad (2.34)$$

From this expression, one verifies that the open/closed partition function is annihilated by the following operator:

$$\left(\frac{1}{\alpha} \frac{\partial}{\partial s} + \frac{\partial}{\partial t_p}\right) \tau^{\text{op+cl}}(\Lambda, s) = \sqrt{\frac{\alpha p}{2\pi}} \int dz \, z^{\frac{p-1}{2}} e^{\alpha \frac{p}{p+1} z^{p+1} - \alpha z^p s} \left(-z^p + \frac{\partial}{\partial t_p}\right) \tau_{\text{ext}}(\Lambda, z). \qquad (2.35)$$

In the second equality, we have used that the dependence of the wave potential on the time $t_p$ comes purely from the $e^\xi$ factor as the purely closed string partition function $\tau(\Lambda_z)$ is independent of the times $t_{np}$. This finishes our proof of the elementary property.

More importantly, we will prove that the wave function $\Psi(t_k)$ which is equal to the ratio of our matrix model integral $\tau^{\text{op+cl}}$ and the original closed string tau-function $\tau(\Lambda)$ does indeed satisfy the differential equations (2.28) and the initial conditions (2.29) in section 3.3. Since the solution to these equations is unique, that will prove our claim (2.32).

$$\frac{\alpha p}{2\pi} \int dz \, z^{p-1} e^{\alpha z^p (s-s')} = \frac{1}{2\pi} \int d(\alpha z^p) e^{\alpha z^p (s-s')} = \delta(s-s'). \qquad (2.33)$$

<footnote n="3">³To obtain this expression, the contour in the z-plane is engineered to give the equality</footnote>
In closing we remark that the matrix model perspective on both the algebraic and stringy duality naturally gives rise to the extended closed string partition function encountered in [20]. This is directly related to the wave function that is canonical from the perspective of the integrable system [41, 42], i.e. it is the wave function that governs the algebraic intersection numbers [21], including their extra time dependencies.

## 3 The open/closed string Virasoro constraints

The KdV integrable hierarchy combined with one Virasoro constraint — the string equation — is sufficient to determine all correlators of topological gravity coupled to matter. Equivalently, the W-algebra constraints on topological gravity correlators are sufficient to determine them all. A subset of the W-algebra constraints are the Virasoro equations, which are strong constraints on the partition function and which allow to make clear contact with the free fermion formulation of the integrable hierarchy, for instance. They will also allow us to make further contact with the geometric framework [21]. Crucially, they are instrumental in proving the identification between our matrix integral and the generator of geometric invariants. The identification of these Virasoro algebras goes back to the study of open/closed string matrix models [13, 14] and has also been touched upon in [16, 22], amongst many other places. We hope our treatment clarifies the various guises of the Virasoro algebra in the literature.

### 3.1 The open/closed Virasoro algebra

The closed string partition function is annihilated by half of a Virasoro algebra [29, 43–46]:

$$L_n \tau(\Lambda_z) = 0 \quad \text{for} \quad n \ge -1. \qquad (3.1)$$

The generators of the Virasoro algebra are given traditionally in terms of the times $t_k$ of the integrable system:

$$ \begin{aligned} L_{-1} &= \alpha \frac{\partial}{\partial t_1} + \sum_{k=p+1}^{\infty} \frac{k}{p} t_k \frac{\partial}{\partial t_{k-p}} + \frac{1}{2p} \sum_{k=1}^{p-1} k(p-k) t_k t_{p-k} \\ L_0 &= \alpha \frac{\partial}{\partial t_{p+1}} + \sum_{k=1}^{\infty} \frac{k}{p} t_k \frac{\partial}{\partial t_k} + \frac{p^2-1}{24p} \\ L_n &= \alpha \frac{\partial}{\partial t_{(n+1)p+1}} + \sum_{k=1}^{\infty} \frac{k}{p} t_k \frac{\partial}{\partial t_{np+k}} + \frac{1}{2p} \sum_{k=1}^{pn-1} \frac{\partial^2}{\partial t_k \partial t_{np-k}}. \end{aligned} \qquad (3.2) $$

The primed sums have indices in the reduced set of times $t_{k \notin pN}$. Here we study a set of extended Virasoro generators $L_n^{\text{ext}}$ that act on the wave potential $\tau_{\text{ext}}(\Lambda, z)$ in such a manner that they become differential operators in the spectral parameter $z$ only. Moreover, after Fourier transform, they become differential operators in the variable $s$. Thus, we will find the following equalities:

$$L_n^{\text{ext}}(t_k) \tau_{\text{ext}} = -L_n^{\text{ext}}(z) \Phi[\tau^{\text{op+cl}}(s)](z) = -\Phi[L_n^{\text{ext}}(s) \tau^{\text{op+cl}}]. \qquad (3.3)$$
Therefore, the wave potential $\tau_{\text{ext}}$ will be annihilated by the generators

$$L_n^{\text{cl}} = L_n^{\text{ext}}(t_k) + L_n^{\text{ext}}(z), \qquad (3.4)$$

while the open/closed partition function $\tau^{\text{op+cl}}$ will be killed by the Virasoro generators:

$$L_n^{\text{op}} = L_n^{\text{ext}}(t_k) + L_n^{\text{ext}}(z). \qquad (3.5)$$

To prove the equalities (3.3), we firstly propose that the extended Virasoro generators are given by the following expressions in terms of the times $t_k$:

$$ \begin{aligned} L_{-1}^{\text{ext}} &= \alpha \frac{\partial}{\partial t_1} + \sum_{k=p+1}^{\infty} \frac{k}{p} t_k \frac{\partial}{\partial t_{k-p}} + \frac{1}{2p} \sum_{k=1}^{p-1} k(p-k) t_k t_{p-k} + t_p \\ L_0^{\text{ext}} &= \alpha \frac{\partial}{\partial t_{p+1}} + \sum_{k=1}^{\infty} \frac{k}{p} t_k \frac{\partial}{\partial t_k} + \left( \frac{p^2-1}{24p} + \frac{1}{2p} \right) \\ L_n^{\text{ext}} &= \alpha \frac{\partial}{\partial t_{(n+1)p+1}} + \sum_{k=1}^{\infty} \frac{k}{p} t_k \frac{\partial}{\partial t_{np+k}} + \frac{1}{2p} \sum_{k=1}^{pn-1} \frac{\partial^2}{\partial t_k \partial t_{np-k}} + \frac{1}{p} \frac{\partial}{\partial t_{np}}. \end{aligned} \qquad (3.6) $$

Note that all times $t_{np}$ appear in these extended Virasoro generators. This is a natural change, since the wave potential depends on all times. Secondly, we observe further additional terms that can formally be argued as follows. If we were to introduce a time $t_0$ such that $t_0 = -\frac{1}{0}\text{Tr}\Lambda^0$ — compare to equation (2.17) —, then $0t_0$ would measure minus the dimension of the matrix $\Lambda$. Since the latter changes by one in our process of integrating out off-diagonal degrees of freedom, the combination $0t_0$ can be argued to change by one. This formally gives rise to the last term in the first and third lines of the proposal (3.6). The additional contribution to the constant term in $L_0$ can then be obtained by computing the $[L_1^{\text{ext}}, L_{-1}^{\text{ext}}]$ commutator. A more solid argument for the proposal is the long calculation that follows.

By explicitly acting with the extended Virasoro generators $L_n^{\text{ext}}(t_k)$ on the wave potential $\tau_{\text{ext}}$, we will manage to represent the extended Virasoro operators as differential operators in $z$ exclusively. Firstly, we calculate the action of the generators $L_{n>1}^{\text{ext}}$ on the wave potential:

$$ \begin{aligned} L_n^{\text{ext}} \tau_{\text{ext}} &= L_n^{\text{ext}} (e^{\xi} \tau(\Lambda_z)) \\ &= \tau(\Lambda_z) L_n^{\text{ext}} e^{\xi} + e^{\xi} L_n^{\text{ext}} \tau(\Lambda_z) + \frac{1}{p} \sum_{k=1}^{pn-1} \frac{\partial e^{\xi}}{\partial t_{np-k}} \frac{\partial \tau(\Lambda_z)}{\partial t_k}. \end{aligned} \qquad (3.7) $$

We have taken into account the cross term that arises due to the double-derivative term in the $L_{n>0}^{\text{ext}}$ generators. Let us consider each of these terms in turn. The first term is proportional to:

$$ L_n^{\text{ext}} e^{\xi} = \left( \alpha z^{np+p+1} + \frac{1}{p} \sum_{k=1}^{\infty} k t_k z^{np+k} + \frac{1}{2p} \sum_{k=1}^{np-1} z^{np} + \frac{1}{p} z^{np} \right) e^{\xi}. \qquad (3.8) $$
Observe that the second term in equation (3.8) can be rewritten as a z-derivative:

$$L_n^{\text{ext}} e^{\xi} = \left( \alpha z^{np+p+1} + \frac{z^{np+1}}{p} \frac{\partial}{\partial z} + \frac{np+1}{2p} z^{np} \right) e^{\xi}. \qquad (3.9)$$

Let us now consider the second term in equation (3.7). The important input here is that the closed string partition function $\tau(\Lambda_z)$ is annihilated by the action of the shifted Virasoro algebra where the shifted times $\tilde{t}_k$ are given by

$$\tilde{t}_k = t_k - \frac{1}{k z^k}. \qquad (3.10)$$

The shifted Virasoro generators take the simple form:

$$L_n^{\text{ext}}(t_k) = L_n^{\text{ext}}(\tilde{t}_k) + \frac{1}{p} \sum_{k=1}^{\infty} \frac{1}{z^k} \frac{\partial}{\partial t_{np+k}}. \qquad (3.11)$$

Acting on $\tau(\Lambda_z)$ the first term gives a zero contribution while the second term can be rewritten as:

$$L_n^{\text{ext}} \tau(\Lambda_z) = \frac{1}{p} z^{np+1} \frac{\partial \tau^{\text{cl}}(\Lambda_z)}{\partial z} - \frac{1}{p} z^{np+1} \sum_{k=1}^{np-1} z^{-k-1} \frac{\partial \tau(\Lambda_z)}{\partial t_k}. \qquad (3.12)$$

Finally, the last term in equation (3.7) can be evaluated to be:

$$\frac{e^{\xi}}{p} z^{np+1} \sum_{k=1}^{np-1} z^{-k-1} \frac{\partial \tau(\Lambda_z)}{\partial t_k}. \qquad (3.13)$$

We see that this contribution cancels the second term in equation (3.12) after multiplication by $e^{\xi}$. Furthermore, the coefficients of the $\partial/\partial z$ terms in formulas (3.9) and (3.12) are the same. Therefore, summing over all terms, we find the following differential operator as the action of the extended Virasoro generators on the wave potential:

$$L_n^{\text{ext}}(t_k) \tau_{\text{ext}} = \left( \alpha z^{np+p+1} + \frac{1}{p} z^{np+1} \frac{\partial}{\partial z} + \frac{np+1}{2p} z^{np} \right) \tau_{\text{ext}}. \qquad (3.14)$$

The differential operator on the right hand side is the negative of what we defined to be $L_n^{\text{ext}}(z)$ in equation (3.4) for the $n > 0$ Virasoro generators. We turn to the remaining two cases. The calculations for $L_0$ are similar to the ones we have already done and we find

$$L_0^{\text{ext}} \tau_{\text{ext}} = \left( \alpha z^{p+1} + \frac{z}{p} \frac{\partial}{\partial z} + \frac{1}{2p} \right) \tau_{\text{ext}}. \qquad (3.15)$$

Note that the constant $1/(2p)$ term survives in the first equality because the shifted Virasoro generator $L_0(\tilde{t}_k)$ that annihilates the closed string partition function $\tau(\Lambda_z)$ does not include this constant term.

The analysis for the operator $L_{-1}^{\text{ext}}$ is more involved due to the terms quadratic in the times. We use again the fact that the operator $L_{-1}^{\text{ext}}(\tilde{t}_k)$, with shifted times, annihilates the
closed string partition function and find the following relation:

$$L_{-1}^{\text{ext}}\tau_{\text{ext}} = \left( \alpha z + \sum_{k=p+1}^{\infty} \frac{k}{p} t_k z^{k-p} \right) \tau_{\text{ext}} + e^{\xi} \sum_{k=p+1}^{\infty} \frac{1}{p} \frac{1}{z^k} \frac{\partial \tau}{\partial t_{k-p}} \\ + \left( t_p + \frac{1}{2p} \sum_{k=1}^{p-1} k(p-k) \left( \frac{t_{p-k}}{kz^k} + \frac{t_k}{(p-k)z^{p-k}} - \frac{1}{k(p-k)z^p} \right) \right) \tau_{\text{ext}}. \quad (3.16)$$

We combine terms linear in the times $t_k$ and use the equations

$$\sum_{k=1}^{\infty} \frac{k}{p} t_k z^{k-1} e^{\xi} = \frac{1}{p} \frac{\partial e^{\xi}}{\partial z} \quad \text{and} \quad \sum_{k=1}^{\infty} z^{-k-1} \frac{\partial \tau(\tilde{t}_k)}{\partial t_k} = \frac{\partial \tau(\tilde{t}_k)}{\partial z}, \quad (3.17)$$

to finally obtain

$$L_{-1}^{\text{ext}} \tau_{\text{ext}} = \alpha \left( z + \frac{1}{\alpha p z^{p-1}} \frac{\partial}{\partial z} - \frac{p-1}{2\alpha p z^p} \right) \tau_{\text{ext}}. \quad (3.18)$$

We note that the $t_p$ term in the extended Virasoro generator (3.6) is crucial in order to obtain the sum over all times that leads to the $\partial_z$-derivative acting on the $e^\xi$ factor. In short, the form of the Virasoro generators (3.14) is valid for all $n$. This concludes our analysis of the extended closed Virasoro generators and how they are represented as differential operators in $z$ on the wave potential. We surmise that indeed the closed string Virasoro generator $L_n^{\text{cl}}$ in equation (3.4) annihilates the wave potential $\tau_{\text{ext}}$.

## 3.2 From closed to open Virasoro

Our next step in proving the equalities (3.3) is to Fourier transform the $z$-differential operator using the open/closed duality equation (2.26) to find the open string realization of the Virasoro algebra. To make things more transparent, we first rewrite the extended Virasoro algebra in (3.14) in the following manner:

$$L_n^{\text{ext}} \tau_{\text{ext}} = \alpha z^{p(n+1)} \left( z + \frac{1}{\alpha p z^{p-1}} \frac{\partial}{\partial z} - \frac{p-1}{2\alpha p z^p} \right) + \frac{n+1}{2} z^{pn}. \quad (3.19)$$

We now make use of the properties of the Fourier transform that are proven in appendix A (see equation (A.5)) that effectively show the equivalence between differential operators of $z$ acting on the closed string side and differential operators of $s$ acting on the open string side. We find the Virasoro generators:

$$\begin{aligned} L_n^{\text{ext}}(s) &= (-\alpha)^{-n} \left( \frac{\partial^{n+1}}{\partial s^{n+1}} \cdot s - \frac{1}{2}(n+1) \frac{\partial^n}{\partial s^n} \right) \\ &= (-\alpha)^{-n} \left( s \frac{\partial^{n+1}}{\partial s^{n+1}} + \frac{1}{2}(n+1) \frac{\partial^n}{\partial s^n} \right). \end{aligned} \quad (3.20)$$

Thus, we have shown that the operators $L_n^{\text{op}}$ defined in equation (3.5) as the sum of the extended Virasoro generators and the $s$ differential operators annihilate the open/closed partition function $\tau^{\text{op+cl}}$.
Let us briefly remark on the string coupling dependence of the Virasoro generators. The terms in equation (3.6) become more transparent when we rescale the times by a factor of $\alpha = -1/g_s$. The resulting string coupling dependence is then $g_s^{-2}$ for the term quadratic in times in $L_{-1}^{\text{ext}}$, arising from the sphere, and there is an extra term arising from the disk, proportional to $g_s^{-1}$. The term with the two derivatives in $L_n^{\text{ext}}$ is proportional to $g_s^2$, and the last, single derivative term is linear in the string coupling $g_s$. The terms above, in equation (3.20) are proportional to $g_s^n$. Each open string insertion comes with an extra factor of the string coupling $g_s$.

**An alternative point of view.** Let us add a remark on how the open/closed Virasoro generators are related to those found in the literature for $p = 2$. See e.g. [19, 22]. We show in appendix B that those Virasoro generators (which we also compute in the appendix for general $p$) are related to the ones we have by conjugation. Indeed, we can conjugate the Virasoro generators to eliminate the $t_{np}$ dependence of the wave potential $\tau_{\text{ext}}$ on which they act. The open/closed Virasoro generators then match the generators of references [19] and [22] that study the case $p = 2$. It should be remarked however (see appendix B) that this point of view introduces either an intricate operator insertion $\det(\partial_s/\alpha + \Lambda^p)^{-1/p}$ inside the matrix expression for the open string partition function, or the ratio of determinants insertion discussed in appendix C for $p = 2$ [19]. Since our matrix model perspective in the bulk of the paper naturally matches both the algebro-geometric considerations as well as the integrability literature, we have stuck to this point of view in section 2 as well as in the calculation of the open/closed Virasoro generators presented in subsections 3.1 and 3.2. Finally, we remind the reader that these Virasoro algebras have an interpretation in terms of a primary operator insertion that dates back to [13, 14].

### 3.3 The integrable hierarchy and the relation to geometry

In this subsection, we tie up several loose ends. We first make good use of the operators we found in our analysis to prove the statement made in subsection 2.3.2 that our open/closed partition function matches the algebro-geometric generating function. Secondly, we note that a subset of our Virasoro constraints correspond to geometric constraints derived in [21], and that our Virasoro algebra extends those constraints to an infinite family.

**The connection to algebraic geometry.** The first part of the proof of the fact that our matrix partition function matches the geometric one is based on a classic result in integrable systems, which says that the Baker-Akhiezer wave function $\psi(t_k, z)$ is given by:

$$ \psi(t_k, z) = \frac{\tau(\tilde{t}_k)}{\tau(t_k)} e^{\xi(t_k, z)} = \frac{\tau_{\text{ext}}}{\tau}, \qquad (3.21) $$

where $\tilde{t}_k$ are the shifted times (3.10). The Baker-Akhiezer wave function $\psi(t_k, z)$ satisfies the differential equations (2.28) as well as the eigenvalue equation

$$ L \psi(t_k, z) = z^p \psi(t_k, z), \qquad (3.22) $$

where $L$ is the Lax operator of the integrable hierarchy. Moreover, from the linear relation between the open/closed partition function $\tau^{\text{op+cl}}$ and the wave potential $\tau_{\text{ext}}$ in (2.34),
it is clear that the ratio $\tau^{\text{op+cl}}/\tau$ is related to the Baker-Akhiezer function $\psi(t_k, z)$ via the Fourier transform. Therefore the ratio $\tau^{\text{op+cl}}/\tau$ also satisfies the differential equations (2.28).

It remains to prove that the wave function $\Psi(t_k) = \tau^{\text{op+cl}}/\tau$ also satisfies the initial conditions (2.29). Our proof closely follows the proof in [42]. Previously, we proved that the wave potential satisfies the equation (see (3.18)):

$$L_{-1}^{\text{ext}} \cdot \tau_{\text{ext}} = \alpha S_z \cdot \tau_{\text{ext}}, \qquad (3.23)$$

where we define $S_z$ to be the differential operator:

$$S_z = \left( z + \frac{1}{\alpha p z^{p-1}} \frac{\partial}{\partial z} - \frac{p-1}{2\alpha p z^p} \right). \qquad (3.24)$$

We now set all times to zero except the time $t_1$ and study the reduced Baker-Akhiezer function $\psi(t_k, z)|_{t_{\ge 2}=0}$. When all times except time $t_1$ are zero, the closed string partition function $\tau(t)$ equals one. Also, recall that the operator $L_{-1}^{\text{ext}}$ coincides with the operator $\alpha \partial_{t_1}$ when all times but the first are zero.<fnref n="4" /> Combining this fact with equation (3.23) and the identification (3.21), we conclude that [42]

$$\left(\frac{\partial}{\partial t_1}\right)^k \psi(t_k, z)|_{t_{\ge 2}=0} = S_z^k \psi(t_k, z)|_{t_{\ge 2}=0} \qquad (3.25)$$

If we therefore define the function $A(z)$ to be the value that the wave function takes at $t_1 = 0$,

$$\psi(t_k, z)|_{t_{\ge 1}=0} = A(z), \qquad (3.26)$$

then we see that by Taylor expansion, we have

$$\psi(t_k, z)|_{t_{\ge 2}=0} = \sum_{n=0}^{\infty} \frac{1}{n!} S_z^n A(z) t_1^n. \qquad (3.27)$$

Moreover, one can check that the function $A(z)$ satisfies the differential constraint [42]:

$$S_z^p \cdot A(z) = z^p A(z). \qquad (3.28)$$

This follows from firstly, the property (3.22) of the Baker-Akhiezer function that it is an eigenvector of the Lax operator, secondly, the initial conditions that require the Lax operator $L$ to be $(\partial_{t_1})^p$ when all times are zero and thirdly, the $L_{-1}^{\text{ext}}$ Virasoro identity (3.23). We now claim that the function $A(z)$ that satisfies the differential constraint (3.28) is given by the following contour integral:

$$A(z) = \sqrt{\frac{\alpha p}{2\pi}} z^{\frac{p-1}{2}} e^{-\alpha \frac{p}{p+1} z^{p+1}} \int ds e^{-\alpha \frac{s^{p+1}}{p+1} + \alpha s z^p}. \qquad (3.29)$$

<footnote n="4">⁴These two statements have to be modified for $p=2$. For this value of $p$, the partition function at zero higher times equals $\tau(t) = \exp \frac{t_1^3}{6}$, and the operator $L_{-1}^{\text{ext}}$ has an extra term $t_1^2/2$. These two modifications cancel each other in the reasoning.</footnote>
We check by explicitly acting with the operator $S_z$:

$$S_z^p \cdot A(z) = \sqrt{\frac{\alpha p}{2\pi}} z^{\frac{p-1}{2}} e^{-\alpha \frac{p}{p+1} z^{p+1}} \int ds s^p e^{-\alpha \frac{s^{p+1}}{p+1} + \alpha s z^p}. \quad (3.30)$$

Within the $s$-integral, one adds and subtracts $z^p$. Then, the combination $(s^p - z^p)$ can be written as a total derivative, which integrates to zero. We have thus proven (3.28). We have chosen the normalization of the function (3.29) such that it is the inverse Fourier transform of the constant function equal to 1.

Recall that our wave function $\Psi(t_k) = \tau^{\text{op+cl}}/\tau$ is the formal Fourier transform of the Baker-Akhiezer function (3.27):

$$\Psi(t_k)|_{t_{\ge 2}=0} = \sqrt{\frac{\alpha p}{2\pi}} \int dz \, z^{\frac{p-1}{2}} e^{\alpha p \frac{z^{p+1}}{p+1}} - \alpha z^p t_p \psi(t_k, z)|_{t_{\ge 2}=0} \quad (3.31)$$

The key point is now that every term in the sum (3.27) for $n > 0$ vanishes inside the integral (3.31) [42]. This can be checked by integration by parts. Therefore the only term that contributes to $\Psi(t_k)|_{t_{\ge 2}=0}$ is the $n = 0$ term. By explicitly substituting the integral expression for the power series $A(z)$, we obtain

$$\Psi(t_k)|_{t_{\ge 2}=0} = \int ds \, e^{-\alpha \frac{s^{p+1}}{p+1}} \frac{\alpha p}{2\pi} \int dz \, z^{p-1} e^{\alpha z^p (s-t_p)} = \int ds \, e^{-\alpha \frac{s^{p+1}}{p+1}} \delta(s-t_p) = 1 \quad (3.32)$$

where the last equality follows from the fact that all higher times, amongst which is $t_p$, are set to zero. Thus, we have proven the initial condition. Therefore we have proven the identification (2.32) between the matrix and the geometric open/closed generating functions.

**Further constraints.** Finally, let us remark that our Virasoro constraints $L_{-1}^{\text{op}}$ and $L_0^{\text{op}}$ correspond to the string equation and dilaton equation of [21] (after using charge conservation, or the dimension constraint). Our Virasoro generators $L_n^{\text{op}}$ extend these geometrically proven constraints to an infinite set.

## 4 Conclusions

Topological gravity coupled to topological matter is well-studied [8]. Introducing D-branes and applying concepts from holography has allowed for a better understanding of these simple string theories, and has provided hands-on illustrations of profound concepts [15, 16]. Certainly, it has been very useful to underpin these achievements with the rigorous mathematical definition of topological gravity theories on Riemann surfaces with boundaries [12]. This has stimulated progress in the identification of open/closed string correlators on Riemann surfaces with boundaries, and their relation to integrable systems [21]. This has in turn allowed for a more rigorous understanding [20] of the notion that boundaries can be replaced by closed string insertions [37].

Amongst the many insights that string theory has provided into topological gravity is the fact that the integration variables in the Kontsevich matrix model correspond to open strings. Indeed, in [17], it was argued that the degrees of freedom in the Kontsevich matrix integral are mesons made of open strings stretching between extended and localized branes.
In string theory, it is often useful to split a given set of D-branes into a heavy stack, and a single “probe brane”. In this work we adopted such a probe brane analysis in treating the Kontsevich matrix model. We integrated out the open strings stretching between a large N stack, and one extra D-brane. In doing so we have naturally produced a matrix model realization of open/closed duality. The process of integrating out generates a determinant operator insertion corresponding to the addition of a D-brane. From the integrable systems perspective this was argued in the case of pure gravity in [18]. It is important to mention that the insertion is distinct from the more familiar determinant operator in the Kazakov matrix model before taking the double scaling limit [17, 27].

The hands-on treatment of the matrix model formulae allowed for a very precise treatment of the duality, the constraints the generating functions satisfy, and their relation to integrable systems. It is this relation to integrable systems that allowed us to precisely match our treatment to the rigorous algebro-geometric treatment of the p-spin intersection numbers on moduli spaces of Riemann surfaces with boundaries. We note that both the extension of the closed string partition function that depends on all times of the KP integrable hierarchy [20] as well as the wave function that appears in the integrable system [42] pop out of our matrix model analysis spontaneously.

One direction for future research amongst many is to carefully match the geometric analysis with a first principle string field theory or conformal field theory derivation of the amplitudes on Riemann surfaces with boundaries.

## Acknowledgments

We thank our colleagues for creating a stimulating research environment, and Alexandr Buryak for patient explanations of the results in [19]. SA would like to thank the École Normale Supérieure, Paris and the Università di Torino, Italy for their hospitality during the completion of this work.

## A Properties of the Fourier transform

We exhibit some basic properties of the generalized Fourier transform that prove useful in the analysis of the Virasoro generators. We recall the definition of the Fourier transform:

$$ \Phi[g(s)](z) := \sqrt{\frac{\alpha p}{2\pi}} z^{\frac{p-1}{2}} \int ds e^{\alpha \left( \frac{z^{p+1}}{p+1} + z^p s \right)} g(s+z). \quad (A.1) $$

We apply the transform to a function g that is an s-derivative, $g(s) = \partial_s f$:

$$ \begin{aligned} \Phi \left[ \frac{\partial f}{\partial s} \right] &= \sqrt{\frac{\alpha p}{2\pi}} z^{\frac{p-1}{2}} \int ds e^{\alpha \left( \frac{z^{p+1}}{p+1} + z^p s \right)} \frac{\partial f(s+z)}{\partial s} \\ &= -\alpha z^p \Phi[f(s)](z), \end{aligned} \quad (A.2) $$

where we have integrated by parts and used convergence properties of the functions under consideration.
Secondly, we begin by calculating the z-derivative of the Fourier transform:

$$ \frac{\partial}{\partial z} \Phi[f(s)](z) = \frac{p-1}{2z} \Phi[f(s)](z) + \sqrt{\frac{\alpha p}{2\pi}} z^{\frac{p-1}{2}} \int ds \left( p \alpha z^{p-1} s + \alpha z^p + \frac{\partial}{\partial z} \right) f(s+z). \quad (A.3) $$

The last two terms cancel against each other by first converting the z-derivative into an s-derivative and then using the relation derived in equation (A.2). In the remaining term, we add and subtract $p \alpha z^p$ within the integrand, and obtain the inverse Fourier transform of the function $s f(s)$ on the right hand side. Rearranging the terms suitably, we finally obtain

$$ \Phi[s f(s)](z) = \left( \frac{1}{\alpha p z^{p-1}} \frac{\partial}{\partial z} - \frac{(p-1)}{2 \alpha p z^p} + z \right) \Phi[f(s)](z). \quad (A.4) $$

We summarize these results schematically:

$$ \begin{gathered} -\alpha z^p \longleftrightarrow \frac{\partial}{\partial s} \\ \left( \frac{1}{\alpha p z^{p-1}} \frac{\partial}{\partial z} - \frac{(p-1)}{2 \alpha p z^p} + z \right) \longleftrightarrow s. \end{gathered} \quad (A.5) $$

## B An equivalent Virasoro algebra

In this appendix, we provide the formulas that accompany the point of view that the (non-extended) closed string partition function should be strictly dependent on the couplings $t_{k \notin pZ}$ only. In particular, this perspective gives rise to another form of the Virasoro algebra acting on the open/closed partition function dual to the pure closed partition function.

Recall that the extended Virasoro generators were obtained by conjugating the usual Virasoro generators by $e^\xi$ given in (2.21). Since this function involved all the times including the $t_{np}$, the extended Virasoro generators necessarily involved the times $t_{np}$ and associated derivatives. This is in contrast to the closed string partition function; defined via the generalized Kontsevich matrix model the partition function is necessarily independent of these times [9, 30]. It is therefore also instructive to define a new wave potential $\tilde{\tau}$ that is obtained by conjugation using an exponent $\tilde{\xi}$ that does not involve these times $t_{np}$, and investigate how the closed-open Virasoro operators are related in this picture. This will make contact with the existing literature [19, 22, 46].

We begin by defining

$$ \tilde{\tau}(z, \Lambda) = e^{\tilde{\xi}} \tau(\Lambda_z), \quad (B.1) $$

where

$$ \begin{aligned} \tilde{\xi} &= \log \det \left(1 - \frac{z}{\Lambda}\right) - \frac{1}{p} \log \det \left(1 - \frac{z^p}{\Lambda^p}\right) \\ &= - \sum_{\substack{k=1 \\ k \neq 0 \bmod p}}^{\infty} \frac{z^k}{k} \mathrm{Tr} \Lambda^{-k} = \sum_{\substack{k=1 \\ k \neq 0 \bmod p}}^{\infty} t_k z^k. \end{aligned} \quad (B.2) $$
The modified statement of open-closed duality is the following compact equation:

$$ \tilde{\tau}(z, \Lambda) = \Phi[\tilde{\tau}^{\text{op+cl}}(s)](z), \qquad (B.3) $$

where $\Phi$ is the same generalized Fourier transform defined previously and the open/closed partition functions are different from the ones obtained earlier due to the extra factor that eliminated the $t_{np}$ times. In particular, we have:

$$ \tilde{\tau}^{\text{op+cl}}(\Lambda, s) = \frac{1}{D(\Lambda)} \int [dM] e^{-\alpha \frac{1}{(p+1)} \text{Tr}[(M+\Lambda)^{p+1}]}_{\ge 2} \mathcal{O} \cdot e^{-\alpha \frac{1}{p+1} s^{p+1}} \det(M+\Lambda-s). \quad (B.4) $$

The operator $\mathcal{O}$ is given by the differential operator $\det(\Lambda^p + \frac{1}{\alpha} \frac{\partial}{\partial s})^{-\frac{1}{p}}$ obtained by translating the $z^p$ operator into an operator acting on the open string variables using the results in appendix A.

The open Virasoro generators one eventually obtains agree with those in the literature as we now show in brief. We begin with the closed string Virasoro generators in equation (3.2). The action of these generators on $\tilde{\xi}$ gives us:

$$ \begin{aligned} L_n e^{\tilde{\xi}} &= \left( \alpha z^{np+p+1} + \frac{1}{p} \sum_{\substack{k=1 \\ k \neq 0 \text{ mod } p}}^{\infty} k t_k z^{np+k} + \frac{1}{2p} \sum_{\substack{k=1 \\ k \neq 0 \text{ mod } p}}^{np-1} z^{np} \right) e^{\tilde{\xi}} \\ &= \left( \alpha z^{np+p+1} + \frac{z^{np+1}}{p} \frac{\partial}{\partial z} + \frac{n(p-1)}{2p} z^{np} \right) e^{\tilde{\xi}} \end{aligned} \qquad (B.5) $$

The main difference compared to the calculation in the main text is that the finite sum in the last term on the right hand side is over all integers from 1 to $np-1$, but not including the integers that are 0 mod $p$. This leads to a different coefficient for the $z^{np}$ term. The rest of the calculation proceeds as before and we obtain

$$ \begin{aligned} L_n \tilde{\tau}(z, \Lambda) &= \left( \alpha z^{np+p+1} + \frac{1}{p} z^{np+1} \frac{\partial}{\partial z} + \frac{n(p-1)}{2p} z^{np} \right) \tilde{\tau}(z, \Lambda) \qquad (B.6) \\ &= \alpha z^{(n+1)p} \left( z + \frac{1}{\alpha p z^{p-1}} \frac{\partial}{\partial z} - \frac{(p-1)}{2\alpha p z^p} \right) \tilde{\tau}(z, \Lambda) + \frac{(n+1)(p-1)}{2p} z^{np} \tilde{\tau}(z, \Lambda). \end{aligned} $$

For the generator $L_0$, we similarly find

$$ \begin{aligned} L_0 \tilde{\tau}(z, \Lambda) &= \left( \alpha z^{p+1} + \frac{1}{p} z \frac{\partial}{\partial z} \right) \tilde{\tau}(z, \Lambda) \\ &= \alpha z^p \left( z + \frac{1}{\alpha p z^{p-1}} \frac{\partial}{\partial z} - \frac{(p-1)}{2\alpha p z^p} \right) \tilde{\tau}(z, \Lambda) + \frac{p-1}{2p} \tilde{\tau}(z, \Lambda). \end{aligned} \qquad (B.7) $$

Lastly, for the $L_{-1}$ operator, we find the same result as before:

$$ L_{-1} \tilde{\tau}(z, \Lambda) = \alpha \left( z + \frac{1}{\alpha p z^{p-1}} \frac{\partial}{\partial z} - \frac{p-1}{2\alpha p z^p} \right) \tilde{\tau}(z, \Lambda). \qquad (B.8) $$
By utilizing the properties of the generalized Fourier transform (A.5), we finally obtain the expression for the open string Virasoro generators acting on the modified open/closed partition function:

$$L_n^{\text{op}}(s) = (-\alpha)^{-n} \left( \frac{\partial^{n+1}}{\partial s^{n+1}} \cdot s - \frac{(p-1)}{2p} (n+1) \frac{\partial^n}{\partial s^n} \right) \\ = (-\alpha)^{-n} \left( s \frac{\partial^{n+1}}{\partial s^{n+1}} + \frac{(p+1)}{2p} (n+1) \frac{\partial^n}{\partial s^n} \right) \qquad (B.9)$$

For $p = 2$, this agrees with the Virasoro generators obtained in [19, 22].

## C Another pure duality

In this appendix, we review an alternative expression of open/closed duality derived in the mathematical literature for the $p = 2$ case. The derivation has features in common with our derivation in section 2, but it also differs in important aspects. The rigorous mathematical treatment of this version of the duality is recent [19–21]. The string theoretic conceptual background is older [15–17]. We refer to the literature for much of the background, and merely draw attention to the alternative derivation, on the one hand, and on the other hand, the consequent differences in the final expression of the duality.

The starting point in [19] is also the Kontsevich matrix model. However the Harish-Chandra-Itzykson-Zuber formula is applied to the matrix integral with an exponent quadratic in the matrix variable of integration.<fnref n="5" /> This differs from the road we pursued in the bulk where we performed the Harish-Chandra-Itzykson-Zuber integration for an exponent linear in the matrix variable of integration for all $p$. This results in a surprisingly different formula for the open/closed string partition function. Firstly, the closed string partition function naturally becomes independent of the even times, after combining with the determinants, and the combination equals a formal Fourier transform of the open/closed partition function — we refer to [19] for the detailed formulas sketched in this appendix:

$$\sqrt{\frac{\det(1 - \frac{z}{\Lambda})}{\det(1 + \frac{z}{\Lambda})}} \tau(\Lambda_z) = \Phi_s^{\text{form}}[\tilde{\tau}^{\text{op+cl}}](\Lambda, z). \qquad (C.1)$$

The formal Fourier transform $\Phi_s^{\text{form}}$ is a close analogue of our transform $\Phi$ [19]. Note that the left hand side matches the function $\tilde{\tau}$ we defined in equation (B.1) when $p = 2$. The modified open string partition function, though, takes the form [19]:

$$\tilde{\tau}^{\text{op+cl}}(\Lambda, z) = \int d\mu(s) \int d\mu(M) \det \left( \frac{\Lambda + \sqrt{\Lambda^2 - z^2} - M + s - z}{\Lambda + \sqrt{\Lambda^2 - z^2} - M - s + z} \right). \qquad (C.2)$$

We observe an insertion of a ratio of determinants, corresponding to both a brane and an anti-brane insertion. It would be interesting to see whether the formula (C.2) allows for a generalization to larger $p$. Moreover, given the equality of the closed string sides in (B.3)

<footnote n="5">⁵This approach is only directly applicable to the case of pure gravity, namely $p = 2$.</footnote>
and (C.1), it would be important to understand the relation between the brane/anti-brane operator in equation (C.2) and the insertion we obtained in expression (B.4). While in the case of pure gravity, the relation can be established by retracing all steps in both lines of reasoning, a more direct physical understanding would be welcome.

**Open Access.** This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.

## References

[1] E. Brézin and V.A. Kazakov, *Exactly Solvable Field Theories of Closed Strings*, Phys. Lett. B 236 (1990) 144 [INSPIRE].

[2] M.R. Douglas and S.H. Shenker, *Strings in Less Than One-Dimension*, Nucl. Phys. B 335 (1990) 635 [INSPIRE].

[3] D.J. Gross and A.A. Migdal, *Nonperturbative Two-Dimensional Quantum Gravity*, Phys. Rev. Lett. 64 (1990) 127 [INSPIRE].

[4] V.G. Knizhnik, A.M. Polyakov and A.B. Zamolodchikov, *Fractal Structure of 2D Quantum Gravity*, Mod. Phys. Lett. A 3 (1988) 819 [INSPIRE].

[5] F. David, *Conformal Field Theories Coupled to 2D Gravity in the Conformal Gauge*, Mod. Phys. Lett. A 3 (1988) 1651 [INSPIRE].

[6] J. Distler and H. Kawai, *Conformal Field Theory and 2D Quantum Gravity*, Nucl. Phys. B 321 (1989) 509 [INSPIRE].

[7] E. Witten, *Topological Quantum Field Theory*, Commun. Math. Phys. 117 (1988) 353 [INSPIRE].

[8] E. Witten, *On the Structure of the Topological Phase of Two-dimensional Gravity*, Nucl. Phys. B 340 (1990) 281 [INSPIRE].

[9] M. Kontsevich, *Intersection theory on the moduli space of curves and the matrix Airy function*, Commun. Math. Phys. 147 (1992) 1 [INSPIRE].

[10] P.H. Ginsparg and G.W. Moore, *Lectures on 2-D gravity and 2-D string theory*, in *Proceedings, Theoretical Advanced Study Institute (TASI 92): From Black Holes and Strings to Particles: Boulder, U.S.A.*, June 1–26, 1992, pp. 277–469, 1993, hep-th/9304011 [INSPIRE].

[11] P. Di Francesco, P.H. Ginsparg and J. Zinn-Justin, *2-D Gravity and random matrices*, Phys. Rept. 254 (1995) 1 [hep-th/9306153] [INSPIRE].

[12] R. Pandharipande, J.P. Solomon and R.J. Tessler, *Intersection theory on moduli of disks*, open KdV and Virasoro, arXiv:1409.2191 [INSPIRE].

[13] S. Dalley, C.V. Johnson, T.R. Morris and A. Watterstam, *Unitary matrix models and 2-D quantum gravity*, Mod. Phys. Lett. A 7 (1992) 2753 [hep-th/9206060] [INSPIRE].

[14] C.V. Johnson, *On integrable $c < 1$ open string theory*, Nucl. Phys. B 414 (1994) 239 [hep-th/9301112] [INSPIRE].

[15] J. McGreevy and H.L. Verlinde, *Strings from tachyons: The $c = 1$ matrix reloaded*, JHEP 12 (2003) 054 [hep-th/0304224] [INSPIRE].
[16] D. Gaiotto and L. Rastelli, *A paradigm of open/closed duality: Liouville D-branes and the Kontsevich model*, JHEP 07 (2005) 053 [hep-th/0312196] [INSPIRE].

[17] J.M. Maldacena, G.W. Moore, N. Seiberg and D. Shih, *Exact vs. semiclassical target space of the minimal string*, JHEP 10 (2004) 020 [hep-th/0408039] [INSPIRE].

[18] A. Alexandrov, *Open intersection numbers, Kontsevich-Penner model and cut-and-join operators*, JHEP 08 (2015) 028 [arXiv:1412.3772] [INSPIRE].

[19] A. Buryak and R.J. Tessler, *Matrix Models and A Proof of the Open Analog of Witten's Conjecture*, Commun. Math. Phys. 353 (2017) 1299 [arXiv:1501.07888] [INSPIRE].

[20] A. Buryak, E. Clader and R.J. Tessler, *Closed extended r-spin theory and the Gelfand-Dickey wave function*, J. Geom. Phys. 137 (2019) 132 arXiv:1710.04829.

[21] A. Buryak, E. Clader and R.J. Tessler, *Open r-spin theory and the Gelfand-Dickey wave function*, J. Geom. Phys. 137 (2019) 132 [arXiv:1809.02536] [INSPIRE].

[22] R. Dijkgraaf and E. Witten, *Developments in Topological Gravity*, Int. J. Mod. Phys. A 33 (2018) 1830029 [arXiv:1804.03275] [INSPIRE].

[23] H. Muraki and C. Rim, *Open KdV hierarchy of 2d minimal gravity of Lee-Yang series*, arXiv:1808.07304 [INSPIRE].

[24] K. Aleshkin and V. Belavin, *Open minimal strings and open Gelfand-Dickey hierarchies*, JHEP 02 (2019) 043 [arXiv:1811.04066] [INSPIRE].

[25] A. Alexandrov, H. Muraki and C. Rim, *From minimal gravity to open intersection theory*, arXiv:1904.06885 [INSPIRE].

[26] E. Brézin and S. Hikami, *On an Airy matrix model with a logarithmic potential*, J. Phys. A 45 (2012) 045203 [arXiv:1108.1958] [INSPIRE].

[27] A. Hashimoto, M.-x. Huang, A. Klemm and D. Shih, *Open/closed string duality for topological gravity with matter*, JHEP 05 (2005) 007 [hep-th/0501141] [INSPIRE].

[28] E. Witten, *Algebraic geometry associated with matrix models of two-dimensional gravity*, proceedings of *Topological methods in modern mathematics*, Stony Brook, NY, U.S.A., 1991, (1993), p. 235.

[29] M. Adler and P. van Moerbeke, *A matrix integral solution to two-dimensional $W_p$-gravity*, Commun. Math. Phys. 147 (1992) 25 [INSPIRE].

[30] C. Itzykson and J.B. Zuber, *Combinatorics of the modular group. 2. The Kontsevich integrals*, Int. J. Mod. Phys. A 7 (1992) 5661 [hep-th/9201001] [INSPIRE].

[31] A. Chiodo, *Stable twisted curves and their r-spin structures*, Annales Inst. Fourier 58 (2008) 1635.

[32] C. Faber, S. Shadrin and D. Zvonkine, *Tautological relations and the r-spin Witten conjecture*, math/0612510.

[33] E. Brézin and S. Hikami, *The intersection numbers of the p-spin curves from random matrix theory*, JHEP 02 (2013) 035 [arXiv:1212.6096] [INSPIRE].

[34] R. Dijkgraaf and E. Witten, *Mean Field Theory, Topological Field Theory and Multimatrix Models*, Nucl. Phys. B 342 (1990) 486 [INSPIRE].

[35] K. Li, *Topological gravity with minimal matter*, Nucl. Phys. B 354 (1991) 711 [INSPIRE].
[36] R. Dijkgraaf, H.L. Verlinde and E.P. Verlinde, *Topological strings in $d < 1$, Nucl. Phys. B 352 (1991) 59* [INSPIRE].

[37] J. Polchinski, *Dirichlet Branes and Ramond-Ramond charges*, Phys. Rev. Lett. 75 (1995) 4724 [hep-th/9510017] [INSPIRE].

[38] Harish-Chandra, *Differential Operators on a Semisimple Lie Algebra*, Am. J. Math. 79 (1957) 87 [INSPIRE].

[39] C. Itzykson and J.B. Zuber, *The Planar Approximation. 2.*, J. Math. Phys. 21 (1980) 411 [INSPIRE].

[40] L.A. Dickey, *Soliton equations and Hamiltonian systems*, Adv. Ser. Math. Phys. 12 (1991) 1.

[41] A. Buryak, *Open intersection numbers and the wave function of the KdV hierarchy*, Moscow Math. J. 16 (2016) 27 [arXiv:1409.7957] [INSPIRE].

[42] M. Bertola and D. Yang, *The partition function of the extended $r$-reduced Kadomtsev-Petviashvili hierarchy*, J. Phys. A 48 (2015) 195205 [arXiv:1411.5717] [INSPIRE].

[43] M. Fukuma, H. Kawai and R. Nakayama, *Continuum Schwinger-dyson Equations and Universal Structures in Two-dimensional Quantum Gravity*, Int. J. Mod. Phys. A 6 (1991) 1385 [INSPIRE].

[44] R. Dijkgraaf, H.L. Verlinde and E.P. Verlinde, *Loop equations and Virasoro constraints in nonperturbative 2-D quantum gravity*, Nucl. Phys. B 348 (1991) 435 [INSPIRE].

[45] S. Kharchev, A. Marshakov, A. Mironov, A. Morozov and A. Zabrodin, *Towards unified theory of 2-D gravity*, Nucl. Phys. B 380 (1992) 181 [hep-th/9201013] [INSPIRE].

[46] R. Dijkgraaf, *Intersection theory, integrable hierarchies and topological field theory*, hep-th/9201003 [INSPIRE].