IPHT-T14/088 CRM-3338 RIKEN-MP-90

## Towards U( N M | ) knot invariant from ABJM theory

Bertrand Eynard a,b † and Taro Kimura a,c ‡

a Institut de Physique Th´orique, CEA Saclay, F-91191 Gif-sur-Yvette, France e b Centre de Recherches Math´matiques, Universit´ de Montr´al, Montr´al, QC, Canada, e e e e c Mathematical Physics Laboratory, RIKEN Nishina Center, Saitama 351-0198, Japan

## Abstract

We study U( N M | ) character expectation value with the supermatrix Chern-Simons theory, known as the ABJM matrix model, with emphasis on its connection to the knot invariant. This average just gives the half-BPS circular Wilson loop expectation value in ABJM theory, which shall correspond to the unknot invariant. We derive the determinantal formula, which gives U( N M | ) character expectation values in terms of U(1 1) averages for a particular type of character representations. | This means that the U(1 1) character expectation value is a building block for the U( | N M | ) averages, and also, by an appropriate limit, for the U( N ) invariants. In addition to the original model, we introduce another supermatrix model obtained through the symplectic transform, which is motivated by the torus knot Chern-Simons matrix model. We obtain the Rosso-Jonestype formula and the spectral curve for this case.

† E-mail address: bertrand.eynard@cea.fr

‡

Present address: Department of Physics, Keio University; E-mail address: taro.kimura@keio.jp Keywords: knot invariant, ABJM theory, Chern-Simons theory, Matrix model MSC: 81T60, 81T45, 57M27, 14H81

## Contents

| 1 Introduction   | 1 Introduction                                   | 1 Introduction                                        |   1 |
|------------------|--------------------------------------------------|-------------------------------------------------------|-----|
| 2                | Unknot matrix model                              | Unknot matrix model                                   |   4 |
|                  | 2.1                                              | U(1 | 1) theory . . . . . . . . . . . . . . . . . . . |   6 |
|                  | 2.2                                              | U( N | N ) theory . . . . . . . . . . . . . . . . . . |   7 |
| 3                | Torus knot matrix model                          | Torus knot matrix model                               |   9 |
| 4                | U( N | M ) theory                                | U( N | M ) theory                                     |  11 |
| 5                | Spectral curve                                   | Spectral curve                                        |  14 |
|                  | 5.1                                              | Symplectic transformation . . . . . . . . . . . .     |  15 |
|                  | 5.2                                              | Saddle point analysis . . . . . . . . . . . . . . .   |  17 |
|                  | 5.3                                              | Asymptotic expansion and topological recursion        |  23 |
| 6                | Comments on related topics                       | Comments on related topics                            |  23 |
|                  | 6.1                                              | Topological A-model . . . . . . . . . . . . . . .     |  23 |
|                  | 6.2                                              | Topological B-model . . . . . . . . . . . . . . .     |  25 |
|                  | 6.3                                              | Matrix models . . . . . . . . . . . . . . . . . .     |  26 |
| 7                | Discussion                                       | Discussion                                            |  28 |
| A                | Mirror of the torus knot ABJM partition function | Mirror of the torus knot ABJM partition function      |  31 |

## 1 Introduction

Since it was shown by Witten [1] that a knot invariant is realized by using the Wilson loop operators in Chern-Simons gauge theory, knot theory has been providing various kinds of interesting topics not only for mathematicians, but also for physicists. In particular the most important example of the knot invariant, which is called the Jones polynomial, is obtained from Chern-Simons theory with the Wilson loop in a fundamental representation

$$J ( K ; q ) = \left \langle W _ { \Box } ( K ; q ) \right \rangle / \left \langle W _ { \Box } ( \bigcirc ; q ) \right \rangle,$$

where the Wilson loop in a representation R is given by

$$W _ { R } ( K ; q ) = \text{Tr} _ { R } \, \text{P} \exp \left ( \oint _ { K } A \right ) \,,$$

and its expectation value is taken with respect to Chern-Simons theory with SU(2) gauge symmetry on a three-sphere S 3

$$S _ { \text{CS} } [ A ] = \frac { k } { 4 \pi } \int _ { S ^ { 3 } } \text{Tr} \left ( A \wedge d A + \frac { 2 } { 3 } A \wedge A \wedge A \right ) \,.$$

In this case the parameter q is associated with the level of Chern-Simons theory as q = exp(2 πi/ k ( +2)). This prescription to derive the knot invariant is quite general: when the fundamental representation R = glyph[square] is replaced by a generic representation, one obtains the colored Jones polynomial, and SU( N ) and SO( N / ) Sp( N ) generalizations provide the HOMFLY and Kauffman polynomials, respectively.

The knot polynomial is usually defined by the Skein relation with a proper normalization condition for the unknot invariant. Although it is in principle computable for any knots based on that definition, their expressions get much complicated as the number of crossings in knots increases, or the representation of the knot polynomial becomes highly involved. Thus it is not yet achieved to write down such a knot invariant in a closed form. On the other hand, for a particular class of knots, the unknot and also torus knots, there is a useful integral representation of the knot invariant [2, 3, 4, 5], which is applicable to a generic gauge group and representation. This is based on the matrix integral formula for the partition function of Chern-Simons theory, especially defined on a three-sphere S 3 , and then given as the expectation value of the character in the corresponding group.

In this paper we consider a supergroup character average with the supermatrix ChernSimons theory, known as the ABJM matrix model [6], towards a supersymmetric generalization of the knot invariant. This supermatrix model is derived from N = 6 superconformal Chern-Simons-matter theory with gauge group U( N ) k × U( N ) -k , which is so-called ABJM theory [7], by implementing the localization technique for the path integral. In this theory the circular Wilson loop, in particular the half-BPS operator, is described by a holonomy with a superconnection taking a value in U( N N | ), which is written as a supergroup character [8]. Thus it is expressed in terms of the supersymmetric Schur function [9, 10]. In this sense the average we compute here is expected to be an unknot invariant for U( N N | ) theory.

The knot invariant and the character average have an analogous structure to a matrix integral in the presence of an external field, which is dual to a correlation function of characteristic polynomials. Especially in the case of supermatrix models, the corresponding correlation function is that for the characteristic polynomial ratio. For this correlation function, there is an interesting determinantal formula which consists of a single pair correlation function as a kernel [11, 12, 13]. In this way we shall expect that a similar determinantal structure can be found in U( N N | ) theory, and the U(1 1) expectation value plays a role as | the kernel function there. We will show in this paper such a determinantal formula for the average with a particular type of the character representation.

According to the general scheme of the topological recursion [14], one can determine all order perturbation series for the correlation function from the spectral curve. In the case of the knot invariant, there are basically two kinds of perturbative expansions. The first is based on the spectral curve, which is obtained from the A-polynomial, and its expansion is from the large representation limit [15, 16]. Although this A-polynomial was originally introduced to the Jones polynomial, namely SU(2) Chern-Simons theory, it can be now extended to more

generic theories. See for example [17]. The other expansion comes from the spectral curve arising in the large rank limit of the knot invariant. This kind of spectral curves is quite analogous to that discussed in random matrix theory, because for some kinds of knots we have matrix integral-like expressions for the knot invariant, and in this case the matrix size N corresponds to the rank of Chern-Simons gauge group SU( N ). This large N limit plays an important role in topological string theory, because it describes the geometric transition of the corresponding Calabi-Yau threefold [18]. Such a duality is still available for the situation in the presence of a knot [19], which gives a brane on a proper Lagrangian submanifold of the Calabi-Yau threefold. In the sense of topological strings, the corresponding spectral curve provides the mirror Calabi-Yau threefold. In this paper we will discuss the large N spectral curve for the supermatrix Chern-Simons theories.

The supergroup character average discussed in this paper is based on ABJM theory. As pointed out in [20], it is perturbatively equivalent to Chern-Simons theory on the lens space L (2 , 1) = S / 3 Z 2 ∼ = RP 3 , which is dual to the topological string on the local P 1 × P 1 geometry, and the associated spectral curve becomes a genus-one curve. Therefore the unknot spectral curve is just given by the mirror curve of the local P 1 × P 1 . In the case of the torus knot, the spectral curve is obtained by applying the symplectic transformation of the unknot [5]. As well as the ordinary HOMFLY polynomial for SU( N ) Chern-Simons theory, we will introduce the ( P, Q )-deformed supermatrix model through the symplectic transform of the original matrix model, which is motivated by the torus knot Chern-Simons matrix model. Since the Adams operation works well even for the Schur function associated with the supergroup U( N N | ), we can derive the Rosso-Jones formula for the torus knot average. This means that the torus knot character average can be represented as a linear combination of the fractionally framed unknot averages. We will also derive the spectral curve from the saddle point equations of the ( P, Q )-deformed supermatrix model, and then obtain a consistent result with the symplectic transform of the unknot curve.

This paper is organized as follows. In Sec. 2 we study the supergroup character average with the ABJM matrix model, which is just given as the half-BPS circular Wilson loop operator in ABJM theory. We especially focus on the representation, corresponding to the partition such that the number of its diagonal components is given by N . We will show the determinantal formula factorizing the U( N N | ) character average into the U(1 1) expectation | values. In Sec. 3 we then consider the ( P, Q )-deformation of the supermatrix model, which is obtained by the SL(2 , Z ) transform of the original ABJM matrix model. We will show that the Rosso-Jones formula holds even for the U( N N | ) theoy, and the U(1 1) average plays a | role of a building block for the ( P, Q )-deformed U( N N | ) character expectation value with a particular kind of representations. In Sec. 4 we extend the argument to U( N M | ) theory, and we will see that the determinantal formula for the character average can be derived even in this case. We obtain the expression which interpolates U( N ) and U( N N | ) theories. In Sec. 5 we discuss the spectral curve for the ( P, Q )-deformed matrix model. We show that there

are two consistent ways of obtaining the spectral curve: One is the symplectic transform of the original matrix model and the other is the saddle point analysis of the ( P, Q )-deformed matrix model itself. In Sec. 6 we give some comments on the relations to topological string and random matrix theory. Sec. 7 is devoted to a summary and discussions.

## 2 Unknot matrix model

Let us start with the matrix model description of Chern-Simons theory at level k defined on a three-sphere S 3 . When we take the gauge group G = U( N ), the partition function is given as the matrix model-like integral

$$\mathcal { Z } _ { \text{CS} } ( S ^ { 3 } ; q ) = \frac { 1 } { N! } \int \prod _ { i = 1 } ^ { N } \frac { d x _ { i } } { 2 \pi } \, e ^ { - \frac { 1 } { 2 g _ { s } } x _ { i } ^ { 2 } } \, \prod _ { i < j } ^ { N } \left ( 2 \sinh \frac { x _ { i } - x _ { j } } { 2 } \right ) ^ { 2 } \,,$$

where the parameter q is related to the coupling constant as q = exp g s with

$$g _ { s } & = \frac { 2 \pi i } { k + N } \,.$$

We then consider the expectation value of the circular Wilson loop operator with representation R , which corresponds to the unknot invariant,

$$\left \langle W _ { R } ( \Box ) \right \rangle = \frac { 1 } { \mathcal { Z } _ { \text{CS} } } \frac { 1 } { N! } \int _ { i = 1 } ^ { N } \frac { d x _ { i } } { 2 \pi } \, e ^ { - \frac { 1 } { 2 g s } x _ { i } ^ { 2 } } \, \prod _ { i < j } ^ { N } \left ( 2 \sinh \frac { x _ { i } - x _ { j } } { 2 } \right ) ^ { 2 } \, \text{Tr} _ { R } \, U ( x ) \,. \quad \, ( 2. 3 )$$

The holonomy matrix U x ( ) is given by

$$U ( x ) = \left ( \begin{array} { c c c } e ^ { x _ { 1 } } & & \\ & \ddots & \\ & & e ^ { x _ { N } } \end{array} \right ) \,,$$

and Tr R U is indeed the character of G = U( N ) in the representation R . This character is written as a Schur polynomial with the corresponding partition λ to the representation R

$$\text{Tr} _ { R } U = s _ { \lambda } ( e ^ { x _ { 1 } }, \cdots, e ^ { x _ { N } } ) \,.$$

Although we can deal with only the unknot Wilson loop based on this matrix model, it is possible to obtain the torus knot invariants by applying a slightly different matrix model, as discussed in Sec. 3.

We now consider a supersymmetric extension of this Chern-Simons theory. In this paper we especially apply the supermatrix generalization of the Chern-Simons matrix model (2.1), namely the ABJM matrix model [8]. It was shown in [6] that ABJM theory [7], which is the three-dimensional superconformal Chern-Simons-matter theory with gauge group U( N ) k ×

U( N ) -k , can be similarly reduced to the matrix model-like integral

$$\mathcal { Z } _ { \text{ABJM} } ( S ^ { 3 } ; q ) \ = \ \frac { 1 } { N! ^ { 2 } } \int _ { i = 1 } ^ { N } \frac { d x _ { i } } { 2 \pi } \frac { d y _ { i } } { 2 \pi } \, e ^ { - \frac { 1 } { 2 g _ { s } } ( x _ { i } ^ { 2 } - y _ { i } ^ { 2 } ) } \, \prod _ { i, j } ^ { N } \left ( 2 \cosh \frac { x _ { i } - y _ { j } } { 2 } \right ) ^ { - 2 } \\ \times \prod _ { i < j } ^ { N } \left ( 2 \sinh \frac { x _ { i } - x _ { j } } { 2 } \right ) ^ { 2 } \left ( 2 \sinh \frac { y _ { i } - y _ { j } } { 2 } \right ) ^ { 2 } \,. \ ( 2. 6 )$$

In this case there is no level shift in the coupling constant

$$g _ { s } & = \frac { 2 \pi i } { k } \,.$$

We can insert the Wilson loop operator into this matrix model as well as Chern-Simons theory with the classical group. Although there are some possibilities for the operators in this case, the relevant choice to this study is the half-BPS circular Wilson loop operator, which is given by a character of the supergroup U( N N | ) [8, 20],

$$\left \langle W _ { R } ( \mathcal { O } ) \right \rangle \ = \ \frac { 1 } { Z _ { \text{ABJM} } } \frac { 1 } { N! ^ { 2 } } \int _ { i = 1 } ^ { N } \frac { d x _ { i } } { 2 \pi } \, \frac { d y _ { i } } { 2 \pi } \, e ^ { - \frac { 1 } { 2 g s } ( x _ { i } ^ { 2 } - y _ { i } ^ { 2 } ) } \, \prod _ { i, j } ^ { N } \left ( 2 \cosh \frac { x _ { i } - y _ { j } } { 2 } \right ) ^ { - 2 } \\ \times \prod _ { i < j } ^ { N } \left ( 2 \sinh \frac { x _ { i } - x _ { j } } { 2 } \right ) ^ { 2 } \left ( 2 \sinh \frac { y _ { i } - y _ { j } } { 2 } \right ) ^ { 2 } \, \text{Str} _ { R } \, U ( x ; y ) \,, \, ( 2. 8 )$$

where the matrix U x y ( ; ) is of the size 2 N × 2 N

$$U ( x ; y ) = \left ( \begin{array} { c c } U ( x ) & \\ & - U ( y ) \end{array} \right ) \,.$$

Let us call this expectation value the unknot Wilson loop average as an analogy with the knot invariant. The supergroup character for U( N N | ) is obtained by replacing the power sum polynomial Tr U n in the U( N + N ) character with the supertrace Str U n [9].

As well as the U( N ) representation theory, the supergroup character can be also expressed as the Schur polynomial, but with a prescribed symmetry [10]

$$\text{Str} _ { R } U ( x ; y ) = s _ { \lambda } ( e ^ { x } ; e ^ { y } ) \,.$$

For U( N N | ) theory, we have a useful determinantal formula in terms of the Frobenius coordinate of the partition λ = ( α , 1 · · · , α d λ ( ) | β , 1 · · · , β d λ ( ) ) with α i = λ i -i and β i = λ t i -i [21]

$$s _ { \lambda } ( u ; v ) = \det _ { 1 \leq i, j \leq d ( \lambda ) } \left ( \sum _ { k, l = 1 } ^ { N } u _ { k } ^ { \alpha _ { i } } \left ( C ^ { - 1 } \right ) _ { k l } v _ { l } ^ { \beta _ { j } } \right ) \,,$$

where the matrix C -1 is the inverse of the Cauchy matrix

$$C = \left ( \frac { 1 } { u _ { k } + v _ { l } } \right ) _ { 1 \leq k, l \leq N }.$$

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0010v2-figure-2261cd56-4650-447e-85f8-419c6e9a53df.jpg)

Figure 1: The partition λ = (12 9 8 6 5 5 4 3 2 2), , , , , , , , , , which is also represented as λ = (11 7 5 2 0 9 8 5 3 1) in the Frobenius coordinate with , , , , | , , , , d λ ( ) = N (= 5). We obtain sub diagrams µ = (7 4 3 1 0) and , , , , ν t = (5 4 3 2 2) involved in this partition. , , , ,

We remark that the supersymmetric Schur polynomial is identically zero when d λ ( ) &gt; N , or equivalently λ N +1 &gt; N .

The formula (2.11) also implies that it can be written only in terms of the hook representations

$$s _ { \lambda } ( u ; v ) = \det _ { 1 \leq i, j \leq d ( \lambda ) } s _ { ( \alpha _ { i } | \beta _ { j } ) } ( u ; v ) \,.$$

This is just a supersymmetric version of the Giambelli formula. Actually this relation is useful to study the Wilson loop operators with various representations in ABJM theory [22].

Especially for the most generic situation with d λ ( ) = N , on which we focus in this paper, the supersymmetric Schur polynomial is decomposed into the ordinary ones [21]

$$s _ { \lambda } ( u ; v ) \ & = \ \det _ { 1 \leq i, j \leq N } u _ { i } ^ { \alpha _ { j } } \cdot \det _ { 1 \leq i, j \leq N } v _ { i } ^ { \beta _ { j } } \cdot \det _ { 1 \leq i, j \leq N } C _ { i j } ^ { - 1 } \\ & = \ s _ { \mu } ( u ) \ s _ { \nu } ( v ) \prod _ { i, j = 1 } ^ { N } ( u _ { i } + v _ { j } ) \,,$$

where the partitions µ and ν are given by

$$\mu _ { i } = \lambda _ { i } - N \,, \quad \nu _ { i } = \lambda _ { i } ^ { t } - N \,, \quad i = 1, \cdots, N \,,$$

or equivalently, µ t i = λ t i + N , ν t i = λ i + N , as shown in Fig. 1. The determinant of the Cauchy matrix is given by

$$\det C = \Delta ( u ) \Delta ( v ) \prod _ { i, j = 1 } ^ { N } ( u _ { i } + v _ { j } ) ^ { - 1 } \,,$$

with the Vandermonde determinant

$$\Delta ( u ) = \prod _ { i < j } ^ { N } ( u _ { i } - u _ { j } ) \,.$$

## 2.1 U(1 1) | theory

Let us consider the simplest example with the supergroup U(1 1), which plays a fundamental | role in this study. In this case only the hook representation is possible, which is written as

λ = ( α β | ) with the Frobenius coordinate. The corresponding Schur function reads

$$s _ { ( \alpha | \beta ) } ( e ^ { x } ; e ^ { y } ) = ( e ^ { x } + e ^ { y } ) \ e ^ { \alpha x + \beta y } \,.$$

The unknot Wilson loop expectation value is given by

$$\left \langle W _ { ( \alpha | \beta ) } ( \mathcal { O } ) \right \rangle = \frac { 1 } { \mathcal { Z } _ { \text{ABJM} } } \int \frac { d x } { 2 \pi } \frac { d y } { 2 \pi } \, e ^ { \frac { i k } { 3 \pi } ( x ^ { 2 } - y ^ { 2 } ) } \left ( 2 \cosh \frac { x - y } { 2 } \right ) ^ { - 2 } ( e ^ { x } + e ^ { y } ) \, e ^ { \alpha x + \beta y } \,. \quad ( 2. 1 9 )$$

We can compute this integral explicitly by applying the Fourier transform formula

$$\frac { 1 } { 2 \cosh w } = \int \frac { d z } { 2 \pi } \, \frac { e ^ { 2 i w z / \pi } } { \cosh z } \,.$$

Thus we have

$$\int \frac { d x } { 2 \pi } \frac { d y } { 2 \pi } \frac { d z } { 2 \pi } \, \frac { 1 } { \cosh z } \, e ^ { \frac { i k } { 4 \pi } ( x ^ { 2 } - y ^ { 2 } ) + ( \alpha + \frac { 1 } { 2 } ) x + ( \beta + \frac { 1 } { 2 } ) y + \frac { 1 } { \pi } ( x - y ) z } = \frac { 1 } { k } \, \frac { q ^ { \frac { 1 } { 2 } ( \alpha + \beta + 1 ) ( \alpha - \beta ) } } { q ^ { \frac { 1 } { 2 } ( \alpha + \beta + 1 ) } + q ^ { - \frac { 1 } { 2 } ( \alpha + \beta + 1 ) } } \,. \ ( 2. 2 1 )$$

Since the partition function becomes Z ABJM = (4 k ) -1 for N = 1, we obtain the expectation value as follows,

$$\ddot { \ w \sigma }, & \left \langle W _ { ( \alpha | \beta ) } ( \mathrm O ) \right \rangle = \frac { 4 \, q ^ { \frac { 1 } { 2 } ( \alpha + \beta + 1 ) ( \alpha - \beta ) } } { q ^ { \frac { 1 } { 2 } ( \alpha + \beta + 1 ) } + q ^ { - \frac { 1 } { 2 } ( \alpha + \beta + 1 ) } } \,. & & ( 2. 2 2 )$$

This depends only on the total and relative lengths of the partition, given by α + β +1 and α -β , respectively. We will see that the numerator can be seen as the framing factor in the following section.

## 2.2 U( N N | ) theory

We then consider the character expectation value for U( N N | ) theory, in particular with a representation with d λ ( ) = N , λ = ( α , 1 · · · , α N | β , 1 · · · , β N ). Let us first rewrite the partition function (2.6)

$$\mathcal { Z } _ { \text{ABJM} } ( S ^ { 3 } ; q ) = \frac { 1 } { N! ^ { 2 } } \int [ d x ] ^ { N } [ d y ] ^ { N } \, \det \left ( \frac { 1 } { 2 \cosh \frac { x _ { i } - y _ { j } } { 2 } } \right ) ^ { 2 } \,,$$

with shorthand notations

$$[ d x ] = \frac { d x } { 2 \pi } \, e ^ { - \frac { 1 } { g _ { s } } x ^ { 2 } } \,, \quad [ d y ] = \frac { d y } { 2 \pi } \, e ^ { \frac { 1 } { g _ { s } } y ^ { 2 } } \,,$$

where we have used the Cauchy formula

$$\det _ { 1 \leq i, j \leq N } \left ( \frac { 1 } { 2 \cosh \frac { x _ { i } - y _ { j } } { 2 } } \right ) = \prod _ { i, j = 1 } ^ { N } \left ( 2 \cosh \frac { x _ { i } - y _ { j } } { 2 } \right ) ^ { - 1 } \prod _ { i < j } ^ { N } \left ( 2 \sinh \frac { x _ { i } - x _ { j } } { 2 } \right ) \left ( 2 \sinh \frac { y _ { i } - y _ { j } } { 2 } \right ) \,.$$

In this case the Schur function has a simple expression as shown in (2.14). Therefore the unnormalized unknot expectation value is now given by

$$\left \langle W _ { R } ( \bigcirc ) \right \rangle = \int [ d x ] ^ { N } [ d y ] ^ { N } \det _ { N } \left ( \frac { 1 } { 2 \cosh \frac { x _ { i } - y _ { j } } { 2 } } \right ) \det _ { N } e ^ { x _ { i } \xi _ { j } } \det _ { N } e ^ { y _ { i } \eta _ { j } } \,,$$

where the parameters ξ i and η i , defined as ξ i = λ i -i +1 2 = / α i +1 2 and / η i = λ t i -i +1 2 = / β i +1 2, play a similar role to external fields in matrix models, as discussed in Sec. 6. / Since all the x i and y i are not distinguishable, this integral can be expressed as a size N determinant

$$\det _ { 1 \leq i, j \leq N } \left [ \int \frac { d x } { 2 \pi } \frac { d y } { 2 \pi } \, e ^ { \frac { i k } { 4 \pi } ( x ^ { 2 } - y ^ { 2 } ) + x \xi _ { i } + y m _ { j } } \left ( 2 \cosh \frac { x - y } { 2 } \right ) ^ { - 1 } \right ] \,.$$

At this moment the computation is almost reduced to that for U(1 1) theory. | Again using the formula (2.20), we obtain the determinantal formula for the character expectation value

$$\left \langle W _ { R } ( \bigcirc ) \right \rangle = k ^ { - N } \prod _ { i = 1 } ^ { N } q ^ { \frac { 1 } { 2 } ( \xi _ { i } ^ { 2 } - \eta _ { i } ^ { 2 } ) } \det _ { 1 \leq i, j \leq N } \left ( \frac { 1 } { q ^ { \frac { 1 } { 2 } ( \xi _ { i } + \eta _ { j } ) } + q ^ { - \frac { 1 } { 2 } ( \xi _ { i } + \eta _ { j } ) } } \right ) \,.$$

This shows that the U( N N | ) unknot character average is factorized into that for U(1 1) | theory (2.22), and thus the measure of this matrix integral is Giambelli compatible in the sense of [12].

To see that the U( N N | ) theory contains the U( N ) knot invariant, it is convenient to rewrite the expression (2.28) by applying the Cauchy formula,

$$\left \langle W _ { R } ( \bigcirc ) \right \rangle = k ^ { - N } q ^ { \frac { 1 } { 2 } ( C _ { 2 } ( \mu ) - C _ { 2 } ( \nu ) ) } \prod _ { i, j = 1 } ^ { N } \left ( q ^ { \frac { 1 } { 2 } ( \alpha _ { i } + \beta _ { j } + 1 ) } + q ^ { - \frac { 1 } { 2 } ( \alpha _ { i } + \beta _ { j } + 1 ) } \right ) ^ { - 1 } \\ \times \prod _ { i < j } ^ { N } \left ( q ^ { \frac { 1 } { 2 } ( \alpha _ { i } - \alpha _ { j } ) } - q ^ { - \frac { 1 } { 2 } ( \alpha _ { i } - \alpha _ { j } ) } \right ) \left ( q ^ { \frac { 1 } { 2 } ( \beta _ { i } - \beta _ { j } ) } - q ^ { - \frac { 1 } { 2 } ( \beta _ { i } - \beta _ { j } ) } \right ). \\. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad.$$

where the 2nd Casimir operator is defined by

$$C _ { 2 } ( \lambda ) = \sum _ { i = 1 } ^ { \infty } \left ( \left ( \lambda _ { i } - i + \frac { 1 } { 2 } \right ) ^ { 2 } - \left ( - i + \frac { 1 } { 2 } \right ) ^ { 2 } \right ) \,.$$

Thus in this case we see that the standard framing factor is given by q 1 2 ( C 2 ( µ ) -C 2 ( ν )) . Up to the normalization constants, the factors in the second line coincide with the Wilson loop expectation value for U( N ) theory, which is given by the quantum dimension of the representation R ,

$$\left \langle W _ { R } ( \bigcirc ) \right \rangle _ { \mathrm U ( N ) } \ & = \ \prod _ { i < j } ^ { N } \frac { [ \lambda _ { i } - \lambda _ { j } - i + j ] _ { q } } { [ - i + j ] _ { q } } \\ & \equiv \ \dim _ { q } R \,,$$

with

$$[ x ] _ { q } = q ^ { x / 2 } - q ^ { - x / 2 } \,.$$

The denominator in (2.31) corresponds to the partition function of U( N ) Chern-Simons theory

$$\overset { \nu J } { \mathcal { Z } } _ { \text{CS} } ( S ^ { 3 } ; q ) = \frac { e ^ { \frac { \pi i } { 8 } N ( N - 1 ) } } { N ^ { \frac {$$

Let us now comment on the situation such that the representation of the character obeys d λ ( ) &lt; N . In this case we do not obtain a simple determinantal formula, since the Schur function with such a representation is not simply factorized any more. For example, the character expectation value for the hook representation, corresponding to the simplest situation d λ ( ) = 1, is given by

$$\left \langle W _ { ( \alpha | \beta ) } ( \mathcal { O } ) \right \rangle = \frac { 1 } { \mathcal { Z } _ { \text{ABJM} } } \frac { 1 } { N! ^ { 2 }$$

where we have

$$\tilde { C } = \left ( \frac { 1 } { 2 \cosh \frac { x _ { i } - y _ { j } } { 2 } } \right ) _ { 1 \leq i, j \leq N }. \quad \quad \quad$$

Although it is difficult to find an explicit formula for this integral, we expect that its asymptotic behavior is obtained from the determinantal formula (2.28) by taking the limit of α , β 1 1 →∞ . Nevertheless, let us comment that we can apply the Fermi gas method even to this case. Actually if we consider the grand canonical partition function, it turns out to be written as a Fredholm determinant, due to the Giambelli formula. See, for example, [22].

## 3 Torus knot matrix model

In addition to the unknot invariant, there is a similar integral formula for the torus knot Wilson loop based on Chern-Simons theory [2, 3, 4, 5],

$$\left \langle W _ { R } ( K _ { P, Q } ) \right \rangle = \frac { 1 } { \mathcal { Z } _ { \text{CS} } ^ { ( P, Q ) } } \frac { 1 } { N! } \int$$

where the coupling constant is now rescaled ˆ g s = PQg s , and the corresponding partition function is then given by

$$\mathcal { Z } _ { \text{CS} } ^ { ( P, Q ) } ( S ^ { 3 } ; q ) = \frac { 1 } { N! } \int \prod _ { i = 1 } ^ { N } \frac { d x _ { i }$$

This is obtained from the ordinary matrix model by applying the SL(2 , Z ) transformation [5]. Note that it is also seen as the biorthogonal generalization of the Chern-Simons matrix model [23].

For the torus knot invariant there is a useful formula, which is called the Rosso-Jones formula [24]

$$\left \langle W _ { R } ( K _ { P, Q } ) \right \rangle = \sum _ { V } c _ { R, Q } ^ { V } \left \langle W _ { V } ( K _ { 1, f } ) \right \rangle,$$

with f = P/Q . This means that the ( P, Q ) torus knot invariant can be expressed as a linear combination of the fractionally framed unknot invariant. This formula is easily derived from the integral formula (3.1) by using the Adams operation

$$s _ { \lambda } ( u ^ { Q } ) = \sum _ { \mu } c _ { \lambda, Q } ^ { \mu } \, s _ { \mu } ( u ) \,.$$

The coefficient c µ λ,Q can be determined by the Frobenius formula for the Schur and power sum polynomials,

$$s _ { \lambda } = \sum _ { \mu } \frac { 1 } { z _ { \mu } } \chi _ { \lambda } ( C _ { \mu } ) \, p _ { \mu } \,, \quad p _ { \mu } = \sum _ { \nu } \chi _ { \nu } ( C _ { \mu } ) \, s _ { \nu } \,,$$

where χ λ and C µ are the character and the conjugacy class for the symmetric group, and the coefficient z µ is given by z µ = ∏ j µ j ! j µ j . Thus we have

$$c _ { \lambda, Q } ^ { \mu } = \sum _ { \nu } \frac { 1 } { z _ { \nu } } \chi _ { \lambda } ( C _ { \nu } ) \chi _ { \mu } ( C _ { Q \nu } ) \,.$$

Since the Rosso-Jones formula (3.3) is obtained from the representation theoretical point of view, it is natural to think that there is a similar formula even for supergroup theories. Actually supersymmetric polynomials obey the same Frobenius formula (3.5) by definition [9]: Let U( x, y ) be a diagonal U( N + M ) matrix, and given an expansion of the U( N + M ) character in terms of the power sum polynomial Tr U x,y ( ) n as

$$\text{Tr} _ { R } U ( x, y ) = \sum _ { R ^ { \prime } } C _ { R } ^ { R ^ { \prime } } \prod _ { i = 1 } ^ { \ell ( R ^ { \prime } ) } \text{Tr} U ( x, y ) ^ { R ^ { \prime } _ { i } }$$

with glyph[lscript] ( R ) to be the number of non-zero entries in the Young diagram R , the supergroup character is just obtained by replacing Tr U x,y ( ) with Str U x y ( ; ),

$$\text{Str} _ { R } U ( x ; y ) = \sum _ { R ^ { \prime } } C _ { R } ^ { R ^ { \prime } } \prod _ { i = 1 } ^ { \ell ( R ^ { \prime } ) } \text{Str} U ( x ; y ) ^ { R ^ { \prime } _ { i } }$$

Therefore we obtain the Adams operation for U( N M | ) theory with the same coefficient (3.6),

$$s _ { \lambda } ( u ^ { Q } ; v ^ { Q } ) = \sum _ { \mu } c _ { \lambda, Q } ^ { \mu } \, s _ { \mu } ( u ; v ) \,.$$

Following the above discussions, we now introduce a supermatrix version of the torus knot matrix model (3.2)

$$\mathcal { Z } _ { \text{ABJM} } ^ { ( P, Q ) } = \frac { 1 } { N! ^ { 2 } } \int [ d x ] ^ { N } [ d y ] ^ { N } \, \det \left ( \frac { 1 } { 2 \cosh \frac { x _ { i } - y _ { j } } { 2 P } } \right ) \det \left ( \frac { 1 } { 2 \cosh \frac { x _ { i } - y _ { j } } { 2 Q } } \right ) \,,$$

with the rescaled coupling constant ˆ g s = PQg s . 1 Then we consider the character expectation value with respect to this partition function

$$\left \langle W _ { R } ( K _ { P, Q } ) \right \rangle = \frac { 1 } { Z _ { \text{ABJM} } ^ { ( P, Q ) } } \, \frac { 1 } { N! ^ { 2 } } \int [ d x ] ^ { N } [ d y ] ^ { N } \, \det \left ( \frac { 1 } { 2 \cosh \frac { x _ { i } - y _ { j } } { 2 P } } \right ) \det \left ( \frac { 1 } { 2 \cosh \frac { x _ { i } - y _ { j } } { 2 Q } } \right ) \, s _ { \lambda } ( e ^ { x } ; e ^ { y } ) \,.$$

Let us call this expectation value the torus knot character average. We can easily show that this ( P, Q )-deformed U( N N | ) character average also satisfies the Rosso-Jones formula

1 We can derive the so-called mirror description for this partition function, as well as the ordinary ABJM matrix model [25]. It depends on the parameters ( P, Q ) in a trivial way (A.3). See Appendix A for details.

(3.3) by applying the Adams operation for U( N N | ) theory (3.9). In this sense it is enough to compute the framed unknot average to obtain the torus knot average. In fact, when we start with a generic partition λ = ( α , 1 · · · , α N | β , 1 · · · , β N ) for a particular torus knot, representations appearing in the expansion (3.3) only provide partitions satisfying d λ ( ) = N . 2 Therefore we now focus on the torus knot and framed unknot with a representation with d λ ( ) = N . In this case we can show that the framed unknot average has the same expression as (2.29), up to the framing factor,

$$& \text{(2.29), up to the training factor,} \\ & \left \langle W _ { R } ( K _ { 1, f } ) \right \rangle = \frac { k ^ { - N } } { \mathcal { Z } ^ { ( 1, 1 ) } \text{$q^{frac{f}{2}(C_{2}(\mu)-C_{2}(\nu))$} } } \det _ { 1 \leq i, j \leq N } \left ( \frac { 1 } { q ^ { \frac { 1 } { 2 } ( \alpha _ { i } + \beta _ { j } + 1 ) } + q ^ { - \frac { 1 } { 2 } ( \alpha _ { i } + \beta _ { j } + 1 ) } } \right ) \\ & \quad = \frac { k ^ { - N } } { \mathcal { Z } ^ { ( 1, 1 ) } \text{$q^{frac{f}{2}(C_{2}(\mu)-C_{2}(\nu))$} } } \prod _ { i, j = 1 } ^ { N } \left ( q ^ { \frac { 1 } { 2 } ( \alpha _ { i } + \beta _ { j } + 1 ) } + q ^ { - \frac { 1 } { 2 } ( \alpha _ { i } + \beta _ { j } + 1 ) } \right ) ^ { - 1 } \\ & \quad \times \prod _ { i < j } ^ { N } \left ( q ^ { \frac { 1 } { 2 } ( \alpha _ { i } - \alpha _ { j } ) } - q ^ { - \frac { 1 } { 2 } ( \alpha _ { i } - \alpha _ { j } ) } \right ) \left ( q ^ { \frac { 1 } { 2 } ( \beta _ { i } - \beta _ { j } ) } - q ^ { - \frac { 1 } { 2 } ( \beta _ { i } - \beta _ { j } ) } \right ) \,. \quad ( 3. 1 2 ) \\ & \text{This means that the II/1/1)\,execution value nlavs a role of the building block for the terms}$$

This means that the U(1 1) expectation value plays a role of the building block for the torus | knot average at least with this kind of representations.

For U(1 1) theory we can compute the character expectation value (3.11) explicitly, as | well as the unknot average (2.22). In this case we have

$$\left \langle W _ { R } ( K _ { P, Q } ) \right \rangle = \frac { 4 \, q ^ { \frac { P Q } { 2 } } ( \alpha + \beta + 1 ) ( \alpha - \beta ) \left ( q ^ { \frac { P Q } { 2 } } ( \alpha + \beta + 1 ) + q ^ { - \frac { P Q } { 2 } } ( \alpha + \beta + 1 ) \right ) } { \left ( q ^ { \frac { P } { 2 } } ( \alpha + \beta + 1 ) + q ^ { - \frac { P } { 2 } } ( \alpha + \beta + 1 ) \right ) } \left ( q ^ { \frac { Q } { 2 } } ( \alpha + \beta + 1 ) + q ^ { - \frac { Q } { 2 } } ( \alpha + \beta + 1 ) \right ) } \,. \quad \text{ (3.13)}$$

Then we obtain 'U(1 1) knot invariant' for the torus knot from this expectation value by | implementing the normalization with the unknot contribution (2.22), as shown in (1.1). Removing the framing factor, we have

$$J _ { R } ( K _ { P, Q } ) = \frac { \left ( q ^ { \frac { 1 } { 2 } ( \alpha + \beta + 1 ) } + q ^ { - \frac { 1 } { 2 } ( \alpha + \beta + 1 ) } \right ) \left ( q ^ { \frac { P Q } { 2 } ( \alpha + \beta + 1 ) } + q ^ { - \frac { P Q } { 2 } ( \alpha + \beta + 1 ) } \right ) } { \left ( q ^ { \frac { P } { 2 } ( \alpha + \beta + 1 ) } + q ^ { - \frac { P } { 2 } ( \alpha + \beta + 1 ) } \right ) \left ( q ^ { \frac { Q } { 2 } ( \alpha + \beta + 1 ) } + q ^ { - \frac { Q } { 2 } ( \alpha + \beta + 1 ) } \right ) } \,.$$

This is a generic formula for the ( P, Q ) torus knot with the representation λ = ( α β | ). This expression is not a polynomial of q and q -1 in general, but it is manifestly invariant under the exchange of P ↔ Q , and also the inversion q ↔ q -1 .

## 4 U( N M | ) theory

The argument shown above can be straightforwardly extended to U( N M | ) theory. The U( N M | ) supermatrix Chern-Simons model is obtained from the Chern-Simons-matter the-

2 Although we do not have an explicit proof of this statement, we check it with a number of examples by numerical calculations.

Figure 2: The partition λ = (14 11 11 9 8 6 5 5 4 3 2 2) satisfying , , , , , , , , , , , λ N ≥ M with N = 7 and M = 5, which includes µ = (9 6 6 4 3 1 0) and , , , , , , ν t = (5 4 3 2 2). , , , ,

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0010v2-figure-4a01ffde-785f-4488-98ba-05a2e2316e01.jpg)

ory with gauge group U( N ) k × U( M ) -k , which is called ABJ theory [26],

$$\mathcal { Z } _ { \text{ABJ} } ( S ^ { 3 } ; q ) = \frac { 1 } { N! M! } \int _ { \real & 1 \, j = 1 } & [ d x ] ^ { N } [ d y ] ^ { M } \prod _ { i = 1 \, j = 1 } ^ { N } \int _ { \real & 1 \, j = 1 } ^ { M } \left ( 2 \cosh \frac { x _ { i } - y _ { j } } { 2 } \right ) ^ { - 2 } \\ & \times \prod _ { i < j } ^ { N } \left ( 2 \sinh \frac { x _ { i } - x _ { j } } { 2 } \right ) ^ { 2 } \prod _ { i < j } ^ { M } \left ( 2 \sinh \frac { y _ { i } - y _ { j } } { 2 } \right ) ^ { 2 }. \\ - \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot$$

The matrix measure in this model can be also expressed as a determinant using the generalized Cauchy determinant formula [27]

$$\Delta _ { N } ( u ) \Delta _ { M } ( v ) \prod _ { i = 1 } ^ { N } \prod _ { j = 1 } ^ { M } ( u _ { i } + v _ { j } ) ^ { - 1 } = \det \left ( \begin{array} { c } u _ { i } ^ { k - 1 } \\ ( u _ { i } + v _ { j } ) ^ { - 1 } \end{array} \right ) \,,$$

where we assume N ≥ M , and the indices run as i = 1 , · · · , N , j = 1 , · · · , M and k = 1 , · · · , N -M . Thus we have

$$& \prod _ { i = 1 } ^ { N } \prod _ { j = 1 } ^ { M } \left ( 2 \cosh \frac { x _ { i } - y _ { j } } { 2 } \right ) ^ { - 1 } \prod _ { i < j } ^ { N } \left ( 2 \sinh \frac { x _ { i } - x _ { j } } { 2 } \right ) \prod _ { i < j } ^ { M } \left ( 2 \sinh \frac { y _ { i } - y _ { j } } { 2 } \right ) \\ & = \prod _ { i = 1 } ^ { N } e ^ { \frac { - N + M + 1 } { 2 } x _ { i } } \prod _ { j = 1 } ^ { M } e ^ { \frac { N - M + 1 } { 2 } y _ { j } } \det \left ( \begin{array} { c } e ^ { x _ { i } ( k - 1 ) } \\ ( e ^ { x _ { i } } + e ^ { y _ { j } } ) ^ { - 1 } \end{array} \right ).$$

We then consider the expectation value of the Wilson loop operator with respect to this partition function. As in the case of U( N N | ) theory, if the partition, corresponding to the representation of the Wilson loop, satisfies λ N ≥ M as shown in Fig. 2, the Schur function is factorized into the ordinary ones [21]

$$s _ { \lambda } ( u ; v ) & = s _ { \mu } ( u ) \, s _ { \nu } ( v ) \prod _ { i = 1 } ^ { N } \prod _ { j = 1 } ^ { M } ( u _ { i } + v _ { j } ) \\ & = \prod _ { i = 1 } ^ { N } u _ { i } ^ { \lambda _ { i } + N - M - i } \prod _ { j = 1 } ^ { M } v _ { j } ^ { \lambda _ { j } ^ { t } + M - N - j } \det \left ( \begin{array} { c } u _ { i } ^ { k - 1 } \\ ( u _ { i } + v _ { j } ) ^ { - 1 } \end{array} \right ) \,,$$

where the partitions µ and ν are defined as µ i = λ i -M for i = 1 , · · · , N and ν j = λ t j -N for j = 1 , · · · , M , or equivalently µ t i = λ t i + M and ν t j = λ j + N . Therefore the unknot character expectation value is now given by

$$& \exp e c { \text{aauon your way is now given } } \upsilon y \\ & \left \langle W _ { R } ( \bigcirc ) \right \rangle = \frac { 1 } { Z _ { \text{ABJ} } } \frac { 1 } { N! M! } \int [ d x ] ^ { N } [ d y ] ^ { M } \det \left ( \begin{array} { c } e ^ { x _ { i } ( k - 1 ) } \\ ( e ^ { x _ { i } } + e ^ { y _ { j } } ) ^ { - 1 } \end{array} \right ) \det e ^ { x _ { i } ( \lambda _ { j } - j + 1 ) } \det e ^ { y _ { i } ( \lambda _ { j } ^ { t } - j + 1 ) } \\ & \quad = \frac { 1 } { Z _ { \text{ABJ} } } \det \left ( \begin{array} { c } \int \frac { d x } { 2 \pi } \, e ^ { - \frac { 1 } { 2 g _ { s } } x ^ { 2 } + x ( \xi _ { i } + k - \frac { 1 } { 2 } ) } \\ \int \frac { d x } { 2 \pi } \, d y \\ \int \frac { d x } { 2 \pi } \, e ^ { - \frac { 1 } { 2 g _ { s } } ( x ^ { 2 } - y ^ { 2 } ) + x \xi _ { i } + y _ { n j } } \left ( 2 \cosh \frac { x - y } { 2 } \right ) ^ { - 1 } \end{array} \right ), \quad ( 4. 5 ) \\ \quad \det \psi _ { \ } c. = \, \vee. \quad \psi _ { \ } t + 1 \, \frac { \epsilon _ { \ } w } { \epsilon _ { \ } w } \, \colon = 1 \, \det \pi _ { \ } w \, \colon = 1 \, \det \pi _ { \ } t \, \det \pi _ { \ } t \, \det \pi _ { \ } t \, \det \pi _ { \ } t \, \det \pi _ { \ } t \, \det \pi _ { \ } t \, \det \pi _ { \ } t \, \det \pi _ { \ } t \, \det \pi _ { \ } t \, \det \pi _ { \ } t \, \det \pi _ { \ } t \, \det \pi _ { \ } t \, \det \pi _ { \ } t \, \det \pi _ { \ } t \, \det \pi _ { \ } t \, \det \pi _ { \ } t \, \det \pi _ { \ } t \, \det \pi _ { \ } t \, \det \pi _ { \ } t \end{array}$$

with ξ i = λ i -i + 1 2 for i = 1 , · · · , N and η j = λ t j -j + 1 2 for j = 1 , · · · , M . This yields the determinantal formula for U( N M | ) theory

$$\left \langle W _ { R } ( \bigcirc ) \right \rangle = \frac { 1 } { Z _ { \text{ABJ} } } \frac { i \frac { N - M } { 2 } } { k \frac { N + M } { 2 } } \prod _ { i = 1 } ^ { N } q ^ { \frac { 1 } { 2 } ( \epsilon _ { i } ^ { 2 } + \xi _ { i } ) } \prod _ { j = 1 } ^ { M } q ^ { \frac { 1 } { 2 } ( - \eta _ { j } ^ { 2 } + \eta _ { j } ) } \prod _ { k = 1 } ^ { N - M } q ^ { \frac { 1 } { 2 } ( k - \frac { 1 } { 2 } ) ^ { 2 } } \det \left ( \begin{array} { c } q ^ { \xi _ { i } ( k - 1 ) } \\ ( q ^ { \xi _ { i } } + q ^ { \eta _ { j } } ) ^ { - 1 } \end{array} \right ).$$

This formula is also represented as follows,

$$& \text{This formula is also represented as follows,} \\ & \left \langle W _ { R } ( \mathcal { O } ) \right \rangle \\ & = \frac { 1 } { \mathcal { Z } _ { A B J } } \frac { \frac { N - M } { 2 } } { \frac { N + M } { 2 } q ^ { \frac { 1 } { 6 } } ( N - M ) ( N - M + \frac { 1 } { 2 } ) ( N - M - \frac { 1 } { 2 } ) } \prod _ { i = 1 } ^ { N } q ^ { \frac { 1 } { 6 } } ( \mathcal { E } _ { i } ^ { 2 } + ( N - M ) \xi _ { i } ) } \prod _ { j = 1 } ^ { M } q ^ { - \frac { 1 } { 6 } } ( \mathcal { E } _ { j } ^ { 2 } + ( N - M ) \eta _ { j } ) \\ & \times \prod _ { i < j } \left ( q ^ { \frac { 1 } { 6 } ( \xi _ { i } - \xi _ { j } ) } - q ^ { - \frac { 1 } { 6 } ( \xi _ { i } - \xi _ { j } ) } \right ) \prod _ { i < j } ^ { M } \left ( q ^ { \frac { 1 } { 6 } ( \eta _ { i } - \eta _ { j } ) } - q ^ { - \frac { 1 } { 6 } ( \eta _ { i } - \eta _ { j } ) } \right ) \prod _ { i = 1 } ^ { M } \prod _ { j = 1 } ^ { M } \left ( q ^ { \frac { 1 } { 6 } ( \xi _ { i } + \eta _ { j } ) } + q ^ { - \frac { 1 } { 6 } ( \xi _ { i } + \eta _ { j } ) } \right ) ^ { - 1 }. \\ & \text{It is easy to see that this expression is reduced to the U ( N | N ) average (2.29) and the U ( N )$$

It is easy to see that this expression is reduced to the U( N N | ) average (2.29) and the U( N ) invariant (2.31) by taking N = M and M = 0, respectively.

We can similarly introduce the U( N M | ) supermatrix model for torus knots

$$\text{we can similarly introduce use } \cup ( \text{v} | \text{$m$} ) \sup e m a l i x { \text{in} \colon } \text{for two two} } \\ \mathcal { Z } _ { \text{ABJ} } ^ { ( P, Q ) } ( S ^ { 3 } ; q ) & = \frac { 1 } { N! M! } \int [ d x ] ^ { N } [ d y ] ^ { M } \prod _ { i = 1 } ^ { N } \prod _ { j = 1 } ^ { M } \left ( 2 \cosh \frac { x _ { i } - y _ { j } } { 2 P } \right ) ^ { - 1 } \left ( 2 \cosh \frac { x _ { i } - y _ { j } } { 2 Q } \right ) ^ { - 1 } \\ & \quad \times \prod _ { i < j } ^ { N } \left ( 2 \sinh \frac { x _ { i } - x _ { j } } { 2 P } \right ) \left ( 2 \sinh \frac { x _ { i } - x _ { j } } { 2 Q } \right ) \prod _ { i < j } ^ { M } \left ( 2 \sinh \frac { y _ { i } - y _ { j } } { 2 P } \right ) \left ( 2 \sinh \frac { y _ { i } - y _ { j } } { 2 Q } \right ). \\ \text{$m$} \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \colon } \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \, \text{$m$} / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / n / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m /m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / s / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / M / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / k / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / l m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m / m.$$

The torus knot average is obtained from this matrix model by inserting the U( N M | ) character, and thus satisfies the Rosso-Jones formula (3.3) thanks to the Adams operation (3.9). We remark that in the limit M → 0 this matrix model reproduces the integral formula for the ordinary U( N ) torus knot invariant, which is the HOMFLY polynomial [2, 3, 4, 5]. This supports a topological nature of the matrix model presented here.

In the next section we will discuss the spectral curve of the matrix model (4.8) arising in the large N limit.

Figure 3: Two ways of obtaining the ( P, Q ) torus knot ABJM (supermatrix) model.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0010v2-figure-66e6ddba-bf61-4c08-b70d-27a898b6495d.jpg)

## 5 Spectral curve

We then provide the spectral curve for the ABJ(M) matrix model for ( P, Q ) torus knots, which has been introduced in Sec. 3 and Sec. 4. In order to obtain the spectral curve, we have to solve the saddle point equations arising in the large N limit of the corresponding matrix model. It is well known that the ABJM matrix model is perturbatively equivalent to the Chern-Simons theory on the lens space L (2 1) = , S / 3 Z 2 ∼ = RP 3 , which is regarded as the two-cut solution of the Chern-Simons matrix model [20]. As well as the ordinary ABJM matrix model, we can obtain the spectral curve for the torus knot ABJM model from the ( P, Q )-modified Chern-Simons theory on the lens space L (2 , 1)

$$\Lambda ^ { - 1 } \, \infty / \, \cfrac { \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \end{array}$$

This is interpreted as the two-cut matrix model of the original Chern-Simons torus knot model (3.2), whose gauge symmetry is broken as as U( N + N ) → U( N ) × U( N ), or in general, U( N + M ) → U( N ) × U( M ).

As shown in Fig. 3, we have two ways of obtaining the spectral curve for the torus knot ABJM matrix model. The first is to apply the SL(2 , Z ) transformation to the spectral curve for ABJM theory, which is given by that for the lens space Chern-Simons theory, through the analytic continuation [28, 29]. We note that this kind of symplectic transformation is

considered to discuss the knot invariant for the lens space [30, 31]. The second is directly solving the saddle point equation for the matrix model (5.1). In this case its saddle point analysis is similar to the ordinary one-cut solution, which is discussed in [5]. We will show that these two methods provide a consistent result.

## 5.1 Symplectic transformation

Let us start with the spectral curve for the Chern-Simons theory on the lens space S / 3 Z 2 ∼ = RP 3 , which is essentially equivalent to that for ABJM theory. The spectral curve C is defined as the zero locus of the two-parameter function

$$\mathcal { C } = \left \{ ( U, V ) \in \mathbb { C } ^ { * } \times \mathbb { C } ^ { * } \, \Big | \, H ( U, V ) = 0 \right \},$$

where the function H U,V ( ) for the lens space L (2 1) is now given by [29] ,

$$H ( U, V ) = c \left ( U + \frac { V ^ { 2 } } { U } \right ) - V ^ { 2 } + \zeta V - 1 \,.$$

The parameter ζ is to be determined, and c is related to the 't Hooft coupling constant as c = exp g s ( N + M / ) 2 for the two-cut model with U( N ) × U( M ) symmetry. 3

Although the algebraic curve characterized by the function (5.3) is naturally obtained from the saddle point analysis of the matrix model, we replace the variable ( U, V ) → ( U , c 2 -1 2 U 1 2 V ) to obtain a homogeneous expression in U and V ,

$$H ( U, V ) = c U ^ { 2 } + V ^ { 2 } - c ^ { - 1 } V ^ { 2 } U ^ { 2 } + c ^ { - \frac { 1 } { 2 } } \zeta V U - 1 = 0 \,.$$

This solves

$$V = \frac { 1 } { 1 - c ^ { - 1 } U ^ { 2 } } \left ( - \frac { \zeta } { 2 c ^ { \frac { 1 } { 2 } } } U \pm \sqrt { \frac { \zeta ^ { 2 } } { 4 c } U ^ { 2 } + ( 1 - c U ^ { 2 } ) ( 1 - c ^ { - 1 } U ^ { 2 } ) } \right ) \,,$$

which is a more convenient expression to perform the contour integral as discussed later. By taking the limit ζ → 0, this solution is just reduced to the one-cut solution, which reproduces the previous result [5], up to a proper replacement of the variable V → V 2 ,

$$V \ \stackrel { \zeta \to 0 } { \longrightarrow } \ \frac { 1 - c U ^ { 2 } } { 1 - c ^ { - 1 } U ^ { 2 } } \,.$$

3 The spectral curve for Chern-Simons theory on the lens space L r, ( 1) = S / 3 Z r (the r -cut Chern-Simons matrix model) is given by [29]

$$\Phi _ { \Phi _ { 1 } \Re \nu _ { 1 } } \\ H ( U, V ) = c \left ( U + \frac { V ^ { r } } { U } \right ) - p _ { r } ( V ) = 0 \,,$$

where p r ( V ) is a degree r polynomial such that the coefficients of V r and V 0 are given by one, p r ( V ) = V r + · · · +1. When the Chern-Simons gauge group is broken as U( N 1 + · · · + N r ) → U( N 1 ) ×··· × U( N r ), the parameter c corresponds to the total 't Hooft coupling c = exp g s ( N 1 + · · · + N r ) / 2. Since the polynomial p r ( V ) has r -1 parameters, the total number of the parameters becomes 1 + ( r -1) = r , which is consistent with that of the subgroups, U( N i ) with i = 1 , · · · , r .

If we write the discriminant of (5.6) as

$$\frac { \zeta ^ { 2 } } { 4 c } U ^ { 2 } + ( 1 - c U ^ { 2 } ) ( 1 - c ^ { - 1 } U ^ { 2 } ) = \left ( U ^ { 2 } - \alpha \right ) \left ( U ^ { 2 } - \frac { 1 } { \alpha } \right ),$$

the end point of the cut is determined by

$$\alpha + \frac { 1 } { \alpha } = c + \frac { 1 } { c } - \frac { \zeta ^ { 2 } } { 4 c } \,.$$

Thus the parameter ζ is seen as the blow-up parameter for the spectral curve, which is determined by requiring that the filling fractions of the cuts are:

$$\frac { 1 } { 2 \pi i } \oint _ { [ \alpha ^ { - 1 / 2 }, \alpha ^ { 1 / 2 } ] } \log V \, \frac { d U } { U } = g _ { s } N \,,$$

$$\frac { 1 } { 2 \pi i } \oint _ { [ - \alpha ^ { 1 / 2 }, - \alpha ^ { - 1 / 2 } ] } ^ { \cdot } \log V \, \frac { d U } { U } = g _ { s } M \,,$$

for U( N ) × U( M ) theory where the integration contour surrounds the corresponding segments counterclockwise.

We can obtain the spectral curve for the ( P, Q ) torus knot from (5.6) through the symplectic transformation, which is characterized by the SL(2 , Z ) matrix [5]

$$M _ { P, Q } = \left ( \begin{array} { c c } Q & P \\ \gamma & \delta \end{array} \right ) \,,$$

where these integer entries satisfy the condition

$$Q \delta - P \gamma = 1 \,.$$

We then apply the following choice of variables corresponding to the matrix (5.12),

$$X & = U ^ { Q } V ^ { P } \,,$$

$$Y = U ^ { \gamma } V ^ { \delta } \,.$$

In this case by substituting the original expression (5.6) and rescaling the variable U 2 → c P Q U 2 , we obtain the spectral curve for the ( P, Q ) torus knot Chern-Simons theory on the lens space

$$X = \frac { U ^ { Q } } { \left ( 1 - c ^ { \frac { P } { Q } - 1 } U ^ { 2 } \right ) ^ { P } } \left ( - \frac { \zeta } { 2 } c ^ { ( \frac { P } { Q } - 1 ) / 2 } U \pm \sqrt { \frac { \zeta ^ { 2 } } { 4 } c ^ { \frac { P } { Q } - 1 } U ^ { 2 } + \left ( 1 - c ^ { \frac { P } { Q } - 1 } U ^ { 2 } \right ) \left ( 1 - c ^ { \frac { P } { Q } + 1 } U ^ { 2 } \right ) } \right ) ^ { P } \,,$$

(5.16)

$$V = \frac { 1 } { 1 - c ^ { \frac { P } { Q } - 1 } U ^ { 2 } } \left ( - \frac { \zeta } { 2 } c ^ { ( \frac { P } { Q } - 1 ) / 2 } U \pm \sqrt { \frac { \zeta ^ { 2 } } { 4 } c ^ { \frac { P } { Q } - 1 } U ^ { 2 } + \left ( 1 - c ^ { \frac { P } { Q } - 1 } U ^ { 2 } \right ) \left ( 1 - c ^ { \frac { P } { Q } + 1 } U ^ { 2 } \right ) } \right ) \,.$$

$$( 5. 1 7 )$$

It is shown in [30, 31] that the new curve obtained through this kind of symplectic transformation gives the topological invariants for torus knots in the lens space S / 3 Z 2 , which is dual to the topological string on the local P 1 × P 1 geometry.

## 5.2 Saddle point analysis

We study the torus knot spectral curve from the large N limit of the matrix model (5.1), and then check its consistency with the result obtained through the symplectic transformation, (5.16) and (5.17). We now rewrite the matrix integral (5.1) with another set of variables, u i = e x / PQ i ( ) and v i = e y / PQ i ( ) ,

$$u _ { i } = e ^ { x _ { i } / ( r \varPsi ) } \text{ and } v _ { i } = e ^ { y _ { i } / ( r \varPsi ) }, \\ \mathcal { Z } _ { \text{CS} } ^ { ( P, Q ) } ( L ( 2, 1 ) ; q ) & = \frac { 1 } { N! ^ { 2 } } \int \frac { d ^ { N } u } { ( 2 \pi ) ^ { N } } \frac { d ^ { N } v } { ( 2 \pi ) ^ { N } } \exp \left [ \sum _ { i, j = 1 } ^ { N } \left ( \log ( u _ { i } ^ { P } + v _ { j } ^ { P } ) + \log ( u _ { i } ^ { Q } + v _ { j } ^ { Q } ) \right ) \\ & \quad + \sum _ { i < j } ^ { N } \left ( \log ( u _ { i } ^ { P } - u _ { j } ^ { P } ) + \log ( u _ { i } ^ { Q } - u _ { j } ^ { Q } ) + \log ( v _ { i } ^ { P } - v _ { j } ^ { P } ) + \log ( v _ { i } ^ { Q } - v _ { j } ^ { Q } ) \right ) \\ & \quad - \sum _ { i = 1 } ^ { N } \left ( \frac { P Q } { 2 g _ { s } } \left ( \log ^ { 2 } u _ { i } + \log ^ { 2 } v _ { i } \right ) + \left ( \frac { P + Q } { 2 } ( 2 N - 1 ) + 1 \right ) \left ( \log u _ { i } + \log v _ { i } \right ) \right ) \right ]. \\ \text{The matrix integral has a convex notential (see [32]\). and this implies that the integrand has}$$

The matrix integral has a convex potential (see [32]), and this implies that the integrand has a unique minimum.

In this case we have two saddle point equations with u i and v i variables, glyph[negationslash]

$$\sum _ { j ( \ne i ) } ^ { N } \left [ \frac { P u _ { i } ^ { P } } { u _ { i } ^ { P } - u _ { j } ^ { P } } + \frac { Q u _ { i } ^ { Q } } { u _ { i } ^ { Q } - u _ { j } ^ { Q } } \right ] + \sum _ { j = 1 } ^ { N } \left [ \frac { P u _ { i } ^ { P } } { u _ { i } ^ { P } + v _ { j } ^ { P } } + \frac { Q u _ { i } ^ { Q } } { u _ { i } ^ { Q } + v _ { j } ^ { Q } } \right ] = \frac { P Q } { g _ { s } } \log u _ { i } + \frac { P + Q } { 2 } ( 2 N - 1 ) + 1 \,, \\ \dots \, - \, - \, - \, - \, - \, - \, - \, - \, - \, - \, -$$

glyph[negationslash]

$$\sum _ { j ( \ne i ) } ^ { N } \left [ \frac { P v _ { i } ^ { P } } { v _ { i } ^ { P } - v _ { j } ^ { P } } + \frac { Q v _ { i } ^ { Q } } { v _ { i } ^ { Q } - v _ { j } ^ { Q } } \right ] + \sum _ { j = 1 } ^ { N } \left [ \frac { P v _ { i } ^ { P } } { v _ { i } ^ { P } + u _ { j } ^ { P } } + \frac { Q v _ { i } ^ { Q } } { v _ { i } ^ { Q } + u _ { j } ^ { Q } } \right ] = \frac { P Q } { g _ { s } } \log v _ { i } + \frac { P + Q } { 2 } ( 2 N - 1 ) + 1 \,.$$

Let us now rewrite these equations in terms of the resolvent [5]. Using the formula

$$\frac { P x ^ { P - 1 } } { x ^ { P } - y ^ { P } } = \sum _ { \substack { k = 0 \\ \sim \. } } ^ { P - 1 } \frac { 1 } { x - \omega ^ { - k Q } y } \,,$$

$$\frac { P x ^ { P - 1 } } { x ^ { P } + y ^ { P } } = \sum _ { k = 0 } ^ { p - 1 } \frac { 1 } { x - \omega ^ { - ( k + \frac { 1 } { 2 } ) Q } y } \,,$$

with the primitive PQ -th root of unity ω = exp2 πi/ PQ ( ), we have glyph[negationslash]

glyph[negationslash]

$$& \sum _ { j ( \neq i ) } ^ { N } \frac { P u _ { i } ^ { P } } { u _ { i } ^ { P } - u _ { j } ^ { P } } + \sum _ { j = 1 } ^ { N } \frac { P u _ { i } ^ { P } } { u _ { i } ^ { P } + v _ { j } ^ { P } } \\ & = \sum _ { j ( \neq i ) } ^ { N } \frac { u _ { i } } { u _ { i } - u _ { j } } + \sum _ { k = 1 } ^ { P - 1 } \left ( \frac { 1 } { g _ { s } } W _ { 0 } ^ { ( 1 ) } ( u _ { i } \omega ^ { k Q } ) - \frac { 1 } { 1 - \omega ^ { - k Q } } \right ) + \sum _ { k = 0 } ^ { P - 1 } \frac { 1 } { g _ { s } } W _ { 0 } ^ { ( 2 ) } ( u _ { i } \omega ^ { ( k + \frac { 1 } { 2 } ) Q ) } ) \,, \quad ( 5. 2 3 )$$

where W ( ) i 0 ( u ) is the leading contribution to the resolvents

$$W ^ { ( 1 ) } ( u ) = \left \langle g _ { s } \sum _ { i = 1 } ^ { N } \frac { u } { u - u _ { i } } \right \rangle = \sum _ { g = 0 } ^ { \infty } g _ { s } ^ { 2 g - 1 } W _ { g } ^ { ( 1 ) } ( u ) \,,$$

$$W ^ { ( 2 ) } ( u ) = \left \langle g _ { s } \sum _ { i = 1 } ^ { N } \frac { u } { u - v _ { i } } \right \rangle = \sum _ { g = 0 } ^ { \infty } g _ { s } ^ { 2 g - 1 } W _ { g } ^ { ( 2 ) } ( u ) \,.$$

These resolvents are analytic in the complex plane except for a finite set of cuts. We now assume that each resolvent has only one cut in the complex plane C ∗ . We can check later that this assumption is correct, using unicity of the extremum guaranteed by [32].

We introduce the 't Hooft coupling constants for this matrix model. If we consider a generic situation with the gauge group U( N 1 + N 2 ) → U( N 1 ) × U( N 2 ), we have two constants

$$t ^ { ( i ) } = g _ { s } N _ { i } \,, \quad i = 1, \, 2 \,.$$

The summation of them is denoted by t = t (1) + t (2) . ABJM theory corresponds to the situation such that N 2 is analytically continuated as N 2 →-N 2 , and then the gauge group ranks are chosen to be N 1 2 , → N . This means that the total 't Hooft coupling becomes zero t = 0. Then the boundary conditions for the resolvents are given by

$$W _ { 0 } ^ { ( i ) } ( u ) \ & \longrightarrow \ \begin{cases} 0 & ( u \to 0 ) \\ t ^ { ( i ) } & ( u \to \infty ) \end{cases},$$

and the saddle point equations in the 't Hooft limit with N →∞ and g s → 0 yield

$$P Q \log u + \frac { P + Q } { 2 } t = W _ { 0 } ^ { ( 1 ) } ( u + i 0 ) + W _ { 0 } ^ { ( 1 ) } ( u - i 0 ) + \sum _ { k = 1 } ^ { P - 1 } W _ { 0 } ^ { ( 1 ) } ( u \omega ^ { k Q } ) + \sum _ { k = 1 } ^ { Q - 1 } W _ { 0 } ^ { ( 1 ) } ( u \omega ^ { k P } ) \\ + \sum _ { k = 0 } ^ { P - 1 } W _ { 0 } ^ { ( 2 ) } ( u \omega ^ { ( k + \frac { 1 } { 2 } ) Q } ) + \sum _ { k = 0 } ^ { Q - 1 } W _ { 0 } ^ { ( 2 ) } ( u \omega ^ { ( k + \frac { 1 } { 2 } ) P } ) \,,$$

$$P Q \log u + \frac { P + Q } { 2 } t = W _ { 0 } ^ { ( 2 ) } ( u + i 0 ) + W _ { 0 } ^ { ( 2 ) } ( u - i 0 ) + \sum _ { k = 1 } ^ { P - 1 } W _ { 0 } ^ { ( 2 ) } ( u \omega ^ { k Q } ) + \sum _ { k = 1 } ^ { Q - 1 } W _ { 0 } ^ { ( 2 ) } ( u \omega ^ { k P } ) \\ + \sum _ { k = 0 } ^ { P - 1 } W _ { 0 } ^ { ( 1 ) } ( u \omega ^ { ( k + \frac { 1 } { 2 } ) Q } ) + \sum _ { k = 0 } ^ { Q - 1 } W _ { 0 } ^ { ( 1 ) } ( u \omega ^ { ( k + \frac { 1 } { 2 } ) P } ) \,. \quad \quad ( 5. 2 9 ) \\ \cdot \quad \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \end{cases}$$

In order to deal with these equations, it is convenient to introduce the exponentiated resolvents

$$y ^ { ( a ) } ( u ) = - u \exp \frac { P + Q } { P Q } \left ( \frac { t } { 2 } - W _ { 0 } ^ { ( 1 ) } ( u ) - W _ { 0 } ^ { ( 2 ) } ( u \, \omega ^ { \frac { 1 } { 2 } Q } ) \right ) \,,$$

$$y ^ { ( b ) } ( u ) = - u \exp \frac { P + Q } { P Q } \left ( \frac { t } { 2 } - W _ { 0 } ^ { ( 1 ) } ( u ) - W _ { 0 } ^ { ( 2 ) } ( u \, \omega ^ { \frac { 1 } { 2 } \, P } ) \right ) \,,$$

′

1

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0010v2-figure-def2ce00-7072-48a9-88a5-bb0ecd70a8f2.jpg)

| F 0           | F 1                 | F 2                 | F 3                           | F 4                           |
|---------------|---------------------|---------------------|-------------------------------|-------------------------------|
| F F ′ 3 F ′ 4 | F 3 F 4 F ′ 3 F ′ 4 | F 3 F 4 F 3 F ′ 4   | F 0 F 1 F 2 F ′ 0 F ′ 1 F ′ 2 | F 0 F 1 F 2 F ′ 0 F F ′ 2     |
| F ′ 0         | F ′ 1               | F ′ 2               | F ′ 3                         | F ′ 4                         |
| F ′ 4 F 3 F 4 | F ′ 4 F ′ 3 F 3 F 4 | F ′ 4 F ′ 3 F 3 F 4 | F ′ 2 F ′ 0 F ′ 1 F 0 F 1 F 2 | F ′ 2 F ′ 0 F ′ 1 F 0 F 1 F 2 |

Table 1: The cuts of the functions F k ( u ) and F Q l + ( u ) for ( P, Q ) = (2 , 3), where F k = F k ( u ), F ′ k = F k ( uω 1 2 ( P -Q ) ) and F ′ Q l + = F Q l + ( uω 1 2 ( Q P -) ). The solid and dotted lines correspond to the cuts from the first resolvent W (1) ( u ) and the second resolvent W (2) ( u ). For example, we can see F 0 = F 3 , F 4 under crossing the corresponding cut of W (1) ( u ), and F 0 = F ′ 3 , F ′ 4 through the cut from W (2) ( u ).

with the boundary behavior

$$y ^ { ( a, b ) } ( u ) \ \longrightarrow \ \begin{cases} - u \, e ^ { \frac { P + Q } { 2 P Q } t } & ( u \to 0 ) \\ - u \, e ^ { - \frac { P + Q } { 2 P Q } t } & ( u \to \infty ) \end{cases}.$$

The saddle point equations are now written as follows,

$$y ^ { ( a ) } ( u + i 0 ) y ^ { ( b ) } ( u - i 0 ) \prod _ { k = 1 } ^ { P - 1 } y ^ { ( a ) } ( u \, \omega ^ { k Q } ) \prod _ { k = 1 } ^ { Q - 1 } y ^ { ( b ) } ( u \, \omega ^ { k P } ) = 1 \,, \quad ( 5. 3 3 )$$

$$y ^ { ( a ) } ( ( u + i 0 ) \omega ^ { \frac { 1 } { 2 } Q } ) y ^ { ( b ) } ( ( u - i 0 ) \omega ^ { \frac { 1 } { 2 } P } ) \prod _ { k = 1 } ^ { P - 1 } y ^ { ( a ) } ( u \omega ^ { ( k + \frac { 1 } { 2 } ) Q } ) \prod _ { k = 1 } ^ { Q - 1 } y ^ { ( b ) } ( u \omega ^ { ( k + \frac { 1 } { 2 } ) P } ) = 1 \,. \quad ( 5. 3 4 )$$

Because it is converted to each other under the exchange of W (1) ( u ) and W (2) ( u ), these equations imply an equivalent condition for the resolvents, which is essentially the same as the one-cut Chern-Simons matrix model for the ( P, Q ) torus knot [5]. Thus we can apply a similar approach to solve these equations.

We consider the products of the resolvents

$$F _ { k } ( u ) = \prod _ { \substack { l = 0 \\ \wedge \,. } } ^ { P - 1 } y ^ { ( a ) } ( u \omega ^ { k P + l Q } ) \,, \quad 0 \leq k \leq Q - 1 \,,$$

$$F _ { Q + l } ( u ) = \prod _ { k = 0 } ^ { Q - 1 } \frac { 1 } { y ^ { ( b ) } ( u \omega ^ { k P + l Q } ) } \,, \quad 0 \leq l \leq P - 1 \,.$$

Since we have assumed that the original resolvents W ( ) i ( u ) have a single cut, these functions F k ( u ) and F Q l + ( u ) have 2 P and 2 Q cuts, obtained by rotating the original one with integer multiple angles of 2 π/ PQ ( ). The total number of the cuts is thus 2 PQ . Due to the saddle point equations they satisfy

$$F _ { k } ( u - i 0 ) \ & = \ \ F _ { Q + l } ( u + i 0 ) \quad \text{ for } \ W ^ { ( 1 ) } ( u ) \,, \\ F _ { k } ( u - i 0 ) \ & = \ \ F _ { Q + l } ( ( u + i 0 ) \omega ^ { \frac { 1 } { 2 } ( Q - P ) } ) \quad \text{for } \ W ^ { ( 2 ) } ( u ) \,, \\ F _ { k } ( ( u - i 0 ) \omega ^ { \frac { 1 } { 2 } ( P - Q ) } ) \ & = \ \ F _ { Q + l } ( ( u + i 0 ) \omega ^ { \frac { 1 } { 2 } ( Q - P ) } ) \quad \text{for } \ W ^ { ( 1 ) } ( u ) \,, \\ F _ { k } ( ( u - i 0 ) \omega ^ { \frac { 1 } { 2 } ( P - Q ) } ) \ & = \ \ F _ { Q + l } ( u + i 0 ) \quad \text{for } \ W ^ { ( 2 ) } ( u ) \,.$$

This means that F k ( u -i 0) = F Q l + ( u + 0) under crossing the cut from the first resolvent i W (1) ( u ), F k ( u -i 0) = F Q l + (( u + 0) i ω 1 2 ( Q P -) ) for the cut from the second W (2) ( u ), and so on. See Table 1 for the case with ( P, Q ) = (2 , 3).

Using these functions we define a function

$$\mathcal { S } ( u, f ) = \mathcal { S } _ { 1 } ( u, f ) \, \mathcal { S } _ { 2 } ( u, f ) \,,$$

where

$$\mathcal { S } _ { 1 } ( u, f ) = \prod _ { \underset { \sim } { \overset { Q - 1 } { \wedge } } } ^ { Q - 1 } \left ( f - F _ { k } ( u ) \right ) \prod _ { l = 0 } ^ { P - 1 } \left ( f - F _ { Q + l } ( u ) \right ),$$

$$\mathcal { S } _ { 2 } ( u, f ) = \prod _ { k = 0 } ^ { Q - 1 } \left ( f - F _ { k } ( u \omega ^ { \frac { 1 } { 2 } ( P - Q ) } ) \right ) \prod _ { l = 0 } ^ { P - 1 } \left ( f - F _ { Q + l } ( u \omega ^ { \frac { 1 } { 2 } ( Q - P ) } ) \right ) \,. \quad \quad ( 5. 4 0 )$$

This function has no cut in the complex plane

$$\mathcal { S } ( u - i 0, f ) = \mathcal { S } ( u + i 0, f ) \,,$$

and the only singularities at u = 0 or u = ∞ as poles. This implies that S u, f ( ) is an entire function of u in C ∗ with polynomial behavior at 0 and at ∞ , therefore it must be a Laurent polynomial of u . Moreover since this function satisfies S ( u, f ) = S ( uω,f ), it depends only on u PQ , and must be a Laurent polynomial of u PQ . We remark that S 1 ( u, f ) and S 2 ( u, f ) still have the cuts, and thus they are not analytic functions.

By definition, S ( u, f ) vanishes when f = F k ( u ), and thus the spectral curve is the algebraic equation:

$$\mathcal { S } ( u, f ) = 0 \,.$$

As shown in [5] we can determine coefficients of the polynomial S ( u, f ) by the asymptotic behavior of F k ( u ) and F Q l + ( u ),

$$F _ { k } ( u ) \, \longrightarrow \, \begin{cases} - \omega ^ { k P ^ { 2 } } \, e ^ { \frac { ( P + Q ) } { 2 Q } t } \, u ^ { P } & ( u \to 0 ) \\ - \omega ^ { k P ^ { 2 } } \, e ^ { - \frac { ( P + Q ) } { 2 Q } t } \, u ^ { P } & ( u \to \infty ) \end{cases},$$

$$F _ { Q + l } ( u ) \, \longrightarrow \, \begin{cases} - \omega ^ { - l Q ^ { 2 } } \, e ^ { - \frac { ( P + Q ) } { 2 P } t } \, u ^ { - Q } & ( u \to 0 ) \\ - \omega ^ { - l Q ^ { 2 } } \, e ^ { \frac { ( P + Q ) } { 2 P } t } \, u ^ { - Q } & ( u \to \infty ) \end{cases}.$$

$$\i e \\ \mathcal { S } _ { i } ( u, f ) = \sum _ { k = 0 } ^ { P + Q } ( - 1 ) ^ { k } \mathcal { S } _ { i, k } ( u ) \, f ^ { P + Q - k } \,, \ \ i = 1, \, 2 \,,$$

this behavior implies

$$\text{ins oceanive} \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{  } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{\ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text
 } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text { } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text` \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{} \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text
}$$

$$\begin{cases} \mathcal { S } _ { i, k } ( u ) = \mathcal { O } ( u ^ { k P } ) \,, & 1 \leq k \leq Q - 1 \,, \\ \mathcal { S } _ { 1, k } ( u ) = ( - 1 ) ^ { P + Q + P Q } e ^ { - \frac { P + Q } { 2 } t } u ^ { P Q } \left ( 1 + \mathcal { O } ( 1 / u ) \right ) \,, & k = Q \,, \\ \mathcal { S } _ { 2, k } ( u ) = ( - 1 ) ^ { P Q } e ^ { - \frac { P + Q } { 2 } t } u ^ { P Q } \left ( 1 + \mathcal { O } ( 1 / u ) \right ) \,, & k = Q \,, \\ \mathcal { S } _ { i, k } ( u ) = \mathcal { O } ( u ^ { P Q } u ^ { - ( k - Q ) Q } ) \,, & Q + 1 \leq k \leq P + Q - 1 \,, \end{cases}$$

(5.47)

Thus we obtain that S ( u, f ) is a polynomial of f and of u PQ

$$\mathcal { S } ( u, f ) = f ^ { 2 ( P + Q ) } + 1 + ( - 1 ) ^ { P + Q } e ^ { - ( P + Q ) t } \left ( u ^ { - 2 P Q } f ^ { 2 Q } + u ^ { 2 P Q } f ^ { 2 P } \right ) + \cdots.$$

where · · · means terms within the Newton polygon, and which are not determined by asymptotic behaviors. The spectral curve for the two-cut matrix model is now given by the following polynomial relation,

$$\mathcal { S } ( u, f ) = 0 \,.$$

It is shown in [32, 33] that if the potential has convexity property, then the solution of the saddle point analysis is unique. This is the case here. We are looking for a two-cut solution, i.e. a genus-one spectral curve. A generic S ( u, f ) with the given Newton's polygon (5.48), would have genus the number of interior points of the Newton's polygon, i.e. 4( P + Q ) -3. However we are here looking for a genus-one curve, which implies that all but one coefficient

If we write

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0010v2-figure-d8e22b7b-78e2-4406-a7f6-5eb7f4b1e24b.jpg)

X

2

1

Figure 4: Newton's polygons for the spectral curve. The polygon in the left panel corresponds to the curve (5.53). There are cetain lattice points inside the polygon, which are fixed by the condition that the curve should have genus-one, and by the filling fraction condition. The right panel shows its singular limit, where the two resolvents W (1) ( u ) and W (2) ( u ) coincide with each other. In this case there is no lattice point inside the polygon.

of S ( u, f ) must be fixed by vanishing of some discriminants. The last coefficient is fixed by the filling fraction condition for U( N M | ), corresponding to (5.10) and (5.11):

$$\frac { 1 } { 2 \pi i } \oint _ { \mathcal { A } } W ( u ) \, d u = g _ { s } N \,, \quad \frac { 1 } { 2 \pi i } \oint _ { \mathcal { B } } W ( u ) \, d u = g _ { s } M \,,$$

where W u ( ) is the sum of the resolvents, and A and B are contours around the first cut and the second cut. One could determine S ( u, f ) by solving all vanishing discriminant equations, so as to impose that the curve have genus-one. However there is a short-cut: one can exhibit an S ( u, f ) with the correct Newton's polygon, the correct filling fraction condition and the correct asymptotic behaviors at 0 and ∞ and which is guaranteed to be genus-one. By unicity according to [32, 33] it must be the correct spectral curve.

The S ( u, f ) which has genus-one and the correct asymptotic behavior is simply the symplectic transform of the unknot of section 5.1, which we write parametrically as [5]

$$f = \frac { V } { U } \,, \quad u ^ { - P Q } = X = U ^ { Q } V ^ { P } \,,$$

we obtain another expression of the spectral curve (5.49) in terms of U and V ,

$$\mathcal { S } ( U, V ) = U ^ { 2 ( P + Q ) } + V ^ { 2 ( P + Q ) } + ( - 1 ) ^ { P + Q } e ^ { - ( P + Q ) t } ( U V ) ^ { 2 ( P + Q ) } + ( - 1 ) ^ { P + Q } e ^ { - ( P + Q ) t } + \cdots = 0 \,,$$

which is equivalent to

$$\mathcal { S } ( X, V ) = V ^ { 2 ( P + Q ) } + X ^ { 2 } + ( - 1 ) ^ { P Q } e ^ { - ( P + Q ) t } X ^ { 2 } V ^ { 2 Q } + ( - 1 ) ^ { P Q } e ^ { - ( P + Q ) t } V ^ { 2 P } + \cdots = 0 \,. \, \text{ (5.53)}$$

P

+

Q

Q

P

V

2

We depict the Newton's polygon for this spectral curve in Fig. 4.

Let us then comment on the singular limit of the spectral curve (5.53), which is realized when the two resolvents W (1) ( u ) and W (2) ( u ) coincide with each other. In this case the analytic function obeys S ( u, f ) = S ( uω ,f 1 2 ). This implies it depends only on u 2 PQ . Therefore the spectral curve (5.53) has to be written only in terms of X 2 and V 2 . The Newton's polygon corresponding to this situation is shown in the right panel of Fig. 4. Since there is no lattice point inside the polygon, the spectral curve has no free parameter. This means that it is just reduced to the genus-zero curve, which corresponds to the limit ζ → 0 of the curve (5.3).

## 5.3 Asymptotic expansion and topological recursion

The matrix integral (3.10) is of the type discussed in [32], therefore it guarantees that it has an asymptotic expansion of the type

$$\log \mathcal { Z } _ { \text{ABJM} } ^ { ( P, Q ) } = \sum _ { g = 0 } ^ { \infty } g _ { s } ^ { 2 g - 2 } F _ { g }$$

which obeys the topological recursion of [14], i.e. F g is the g th symplectic invariant of the spectral curve as defined by the topological recursion [14]. Similarly, all expectation values have a g s -expansion, whose coefficients are given by the topological recursion.

The character expectation values 〈 Tr R U 〉 can be decomposed on the basis of power sums. We have:

$$\left \langle \prod _ { i = 1 } ^ { n } \text{Tr} \, U ^ { p _ { i } } \right \rangle = \text{Res} _ { x _ { 1 } \rightarrow 0 } \dots \text{Res} _ { x _ { n } \rightarrow 0 } W _ { n } ( x _ { 1 }, \dots, x _ { n } ) \, x _ { 1 } ^ { p _ { 1 } } \dots x _ { n } ^ { p _ { n } } \prod _ { i = 1 } ^ { n } \frac { d x _ { i } } { x _ { i } } \, \quad \ \ ( 5. 5 5 )$$

and W n has a topological expansion:

$$W _ { n } ( x _ { 1 }, \dots, x _ { n } ) = \sum _ { g = 0 } ^ { \infty } g _ { s } ^ { 2 g - 2 + n } W _ { g, n } ( x _ { 1 }, \dots, x _ { n } )$$

and W g,n is computed by the topological recursion.

## 6 Comments on related topics

## 6.1 Topological A-model

Let us comment on the realization of the knot invariant in topological string theory, especially in the topological A-model. It was shown in [19] that the knot invariant for K is obtained in the topological string by adding a brane with a proper Lagrangian submanifold of the Calabi-Yau threefold L K . This insertion of a brane gives rise to the expectation value of the

characteristic polynomial with respect to Chern-Simons theory

$$\mathcal { Z } _ { \text{top} } ( K ; x ) & = \left \langle \det \left ( 1 \otimes 1 - U \otimes e ^ { - x } \right ) \right \rangle _ { \text{CS} } \\ & = \sum _ { n = 0 } ^ { \infty } \left \langle \text{Tr} _ { R _ { n } } U \right \rangle _ { \text{CS} } e ^ { - n x } \,,$$

where the matrix U is the holonomy along the knot K U , = Pexp (∮ K A ) , and R n is the totally symmetric representation with n boxes. This means that this topological string partition function is the discrete Fourier (Laplace) transform of the HOMFLY polynomial, since the expectation value of the holonomy 〈 Tr R n U 〉 just gives the knot invariant. If we consider a multi-point correlator of the characteristic polynomials, the knot invariant with more generic representations is obtained. Note that this knot invariant is analogous to the matrix integral with the external source, as discussed in Sec. 6.3. This kind of relation between the characteristic polynomial and the external source is naturally interpreted from the viewpoint of the topological expansion of spectral curves.

In our case it is natural to consider a supermatrix version of the partition function (6.1), corresponding to ABJM theory

$$\mathcal { Z } _ { \text{top} } ( K ; x, y ) = \left \langle \text{Sdet} \left ( 1 \otimes 1 - U \otimes \left ( \begin{array} { c c } e ^ { - x } \\ & e ^ { - y } \end{array} \right ) \right ) \right \rangle _ { \text{ABJM} } \,.$$

Actually one can obtain this partition function by applying both of bosonic and fermionic modes, describing strings stretching between the three-sphere S 3 and the Lagrangian L K

$$\mathcal { Z } _ { \text{top} } ( K ; x, y ) = \left \langle \frac { \det \left ( 1 \otimes 1 - U ^ { \prime } \otimes e ^ { - x } \right ) } { \det \left ( 1 \otimes 1 - U ^ { \prime \prime } \otimes e ^ { - y } \right ) } \right \rangle$$

where U = diag( U , U ′ ′′ ). Whether it comes in denominator (bosonic) or numerator (fermionic) depends on the relative orientation of S 3 and L K [34]. See also the argument on the characteristic polynomial ratio presented in Sec. 6.3.

We can similarly expand (6.2) with the corresponding expectation value 〈 Str R n U 〉 ABJM . This partition function is also expanded with the string coupling constant in the sense of the WKB expansion [19]. Using the identity

$$\text{Sdet} \left ( 1 \otimes 1 - U \otimes \left ( \begin{array} { c } e ^ { - x } \\ & e ^ { - y } \end{array} \right ) \right ) = \exp \left [ \text{Str log} \left ( 1 \otimes 1 - U \otimes \left ( \begin{array} { c } e ^ { - x } \\ & e ^ { - y } \end{array} \right ) \right ) \right ] \,,$$

we have

$$\mathcal { Z } _ { \text{top} } ( K ; x, y ) \, \sim \, \exp \left ( \frac { 1 } { g _ { s } } \int _ { y } ^ { x } p ( x ) d x \right ) \,,$$

where the integrand is given by

$$p ( x ) = \lim _ { g _ { s } \to 0 } \sum _ { n = 0 } ^ { \infty } g _ { s } \left \langle \text{Str} \, U ^ { n } \right \rangle _ { \text{ABJM} } e ^ { - n x } \,.$$

In this way we can show that its leading contribution is given by the disc amplitude as well as the ordinary knot invariant. In this case, since the ABJM matrix model is obtained from Chern-Simons theory on the lens space L (2 1) through the analytic continuation, the , mirror curve, on which the one-form is defined, is replaced with that for the local P 1 × P 1 geometry. Actually the Wilson loop expectation value is evaluated based on this spectral curve [20, 35]. Furthermore, in (6.5), both of the initial and end points of the integral have physical meanings, as positions of brane and anti-brane. Thus the partition function describes pair creation of branes in the topological string. If we take the limit y → ∞ , it goes back to the usual one, including either of bosonic or fermionic modes.

## 6.2 Topological B-model

In the B-model description of the topological strings, the n -point function defined in (5.56) plays an important role. As well as the one-point function W 1 = W which is given as the resolvent, we also have a series expansion of the multi-point function with respect to the coupling constant. If once a spectral curve is obtained, one can determine higher order terms by using the topological recursion [14]. In this study the spectral curve is given by the genus-one mirror curve for the local P 1 × P 1 geometry and its symplectic transform. Indeed this expansion has a natural interpretation in terms of the B-model topological strings. The multi-point correlation function corresponds to multiple insertion of the Wilson loop operators into the Chern-Simons matrix model. Thus a set of variables in the multi-point function (5.56) provides the boundary condition for topological strings. This implies that it computes the open string sector of the B-model, and we can obtain the corresponding open Gromov-Witten invariants based on the mirror symmetry in a perturbative way [36, 37, 38].

As in the case of the topological A-model, the knot invariant can be investigated also in the B-model through the mirror symmetry. In addition to the unknot invariants [39, 40], the torus knot invariant is formulated by using the symplectic transformation acting on the Bmodel open string moduli [5]. 4 Moreover this prescription for giving the torus knot invariant is now applied to the knot invariant not only in the three-sphere S 3 , but also in the lens space L (2 , 1) [30, 31]. Since the ABJM matrix model is perturbatively equivalent to Chern-Simons theory on the lens space L (2 1), the perturbative analysis of the knot invariant for , L (2 1), , which is based on the topological recursion, is apparently relevant to our supermatrix model for the torus knot introduced in Sec. 3. Although we can obtain a systematic expansion of the correlation functions, we have to also take into account non-perturbative contributions, which play an important role in knot theory [15, 16] and ABJM theory [42]. The study of such a non-perturbative effect on the torus knots is an interesting and important issue to be clarified in the future.

4 Let us note that another kind of approach to the B-model description, which is in principle applicable to any knots, is discussed based on the A-polynomial [41].

## 6.3 Matrix models

We comment on the idea underlying this work, which comes from random matrices. In random matrix theory, the wave functions are expectation values of characteristic polynomials

$$\psi ( x ) = \left \langle \det ( x - M ) \right \rangle,$$

where the expectation value is taken with respect to a certain measure, which is specified later. The Hamiltonian associated with this wave function is given by the non-commutative Riemann surface, which is obtained through quantization of the spectral curve. In fact, a more important object is the kernel

$$\mathcal { K } ( x ; y ) = \frac { 1 } { y - x } \left \langle \frac { \det ( x - M ) } { \det ( y - M ) } \right \rangle \,.$$

Out of that kernel one can reconstruct every other observables. For instance, if the matrix size is given by N , the wave function is obtained by sending y →∞ ,

$$\psi ( x ) = \lim _ { y \to \infty } y ^ { N + 1 } \, \mathcal { K } ( x ; y ) \,.$$

Every other correlations of characteristic polynomials is obtained by the Fay identity [11] (See also [12, 13])

$$\left \langle \prod _ { i = 1 } ^ { k } \frac { \det ( x _ { i } - M ) } { \det ( y _ { i } - M ) } \right \rangle = \frac { \prod _ { i, j } ( y _ { i } - x _ { j } ) } { \prod _ { i < j } ( x _ { i } - x _ { j } ) \, ( y _ { i } - y _ { j } ) } \, \det _ { 1 \leq i, j \leq k } \mathcal { K } ( x _ { i } ; y _ { j } ) \,,$$

which is also seen as Pl¨cker or Hirota equation. u The factor in front of the determinant of the kernel can be written as the Cauchy determinant (2.16).

Wave functions and kernels can be seen themselves as partition functions. Indeed let

$$\mathcal { Z } = \int d \mu ( M )$$

be a matrix integral with some measure dµ M ( ) depending on a coupling constant g s . In many cases, there is a small g s expansion of the type

$$\log \mathcal { Z } = \sum _ { g = 0 } ^ { \infty } g _ { s } ^ { 2 g - 2 } \, F _ { g } ( \mathcal { C } )$$

where C is the spectral curve associated to the measure dµ , i.e. the g s → 0 limit of the eigenvalue distribution.

An expectation value of characteristic polynomials can be written

$$\mathcal { K } ( x _ { 1 }, \dots, x _ { k } ; y _ { 1 }, \dots, y _ { k } ) = \frac { 1 } { \mathcal { Z } } \int d \mu ( M ) \prod _ { i = 1 } ^ { k } \frac { \det ( x _ { i } - M ) } { \det ( y _ { i } - M ) } \,,$$

i.e. it is the partition function with a new measure

$$d \mu _ { x ; y } ( M ) = \exp \sum _ { i = 1 } ^ { k } \left ( \, \text{Tr} \log ( x _ { i } - M ) - \text{Tr} \log ( y _ { i } - M ) \right ) d \mu ( M ) \,. \quad \quad \, ( 6. 1 4 )$$

There is also a g s expansion with a new spectral curve

$$\log \mathcal { K } ( x _ { 1 }, \dots, x _ { k } ; y _ { 1 }, \dots, y _ { k } ) = \sum _ { g = 0 } ^ { \infty } g _ { s } ^ { 2 g - 2 } \left ( F _ { g } ( \mathcal { C } _ { x ; y } ) - F _ { g } ( \mathcal { C } ) \right ).$$

The way to find the new spectral curve is as follows. The spectral curve C is a Riemann surface equipped with two analytic functions u and v as C = { ( u, v ) ∈ C × C | H u,v ( ) = 0 } . The new spectral curve C x y ; is a Riemann surface such that vdu has additional simple poles at the x i 's (residue +1) and at the y i 's (residue -1).

Instead of characteristic polynomials one may also be interested in external fields interactions

$$\langle \mathrm e ^ { \text{Tr} \, M A } \rangle = \frac { 1 } { \mathcal { Z } } \int d \tilde { \mu } ( M ) \, \mathrm e ^ { \text{Tr} \, M A } \,.$$

Again, this has (under some assumptions) a g s expansion with a new spectral curve C A

$$\log \left \langle e ^ { \text{Tr} \, M A } \right \rangle = \sum _ { g = 0 } ^ { \infty } g _ { s } ^ { 2 g - 2 } \, \left ( F _ { g } ( \mathcal { C } _ { A } ) - F _ { g } ( \mathcal { C } ) \right ) \,.$$

The way to find the new spectral curve is as follows. The spectral curve C is a Riemann surface equipped with two analytic functions u and v . The new spectral curve C A is a Riemann surface such that udv has additional simple poles at the a i 's (the eigenvalues of A ) with a residue equal to the multiplicity of a i .

We see that under the exchange u ↔ v and by identifying the a i 's with the x i 's (multiplicity α i = 1) and the y i 's (multiplicity α i = -1), this would be the same spectral curve and thus

$$\left \langle e ^ { \text{Tr} \, M A } \right \rangle _ { \tilde { \mu } } \in \left \langle \prod _ { i = 1 } ^ { k } \det ( x _ { i } - M ) ^ { \alpha _ { i } } \right \rangle _ { \mu }$$

where the measures µ and ˜ are obtained by exchanging the role of µ u and v in the spectral curves, and A is a matrix with eigenvalues x i with multiplicity α i [14, 43, 44, 45, 46]. This raises the question of how can a multiplicity be negative? This is why supermatrix models are needed. Negative multiplicities correspond to the fermionic side of the supermatrix [47]. This duality between external field and expectation values of characteristic polynomials appeared in [48].

As shown in (2.26), the half-BPS Wilson loop expectation value in ABJM theory has a quite similar expression to that of matrix integral with the external source. The Gaussian matrix model with the external source is given as follows,

$$\left \langle e ^ { \text{Tr} \, M A } \right \rangle & = \frac { 1 } { \mathcal { Z } } \int d M \, e ^ { - \frac { 1 } { 2 g _ { s } } \, \text{Tr} \, M ^ { 2 } + \text{Tr} \, M A } \\ & = \frac { 1 } { \mathcal { Z } } \frac { 1 } { \Delta ( a ) } \int \prod _ { i = 1 } ^ { N } \frac { d x _ { i } } { 2 \pi } \, e ^ { - \frac { 1 } { 2 g _ { s } } x _ { i } ^ { 2 } } \det _ { 1 \leq i, j \leq N } e ^ { x _ { i } a _ { j } } \prod _ { i < j } ^ { N } ( x _ { i } - x _ { j } ) \,,$$

where the U( N ) part is integrated out using the Harish-Chandra-Itzykson-Zuber formula [49, 50]. Writing the U( N ) character in terms of the Schur function

$$s _ { \lambda } ( e ^ { x } ) = \frac { 1 } { \Delta ( e ^ { x } ) } \det _ { 1 \leq i, j \leq N } e ^ { x _ { i } ( \lambda _ { j } + N - j ) } \,,$$

we obtain the Wilson loop expectation value with respect to the U( N ) Chern-Simons theory

$$\left \langle W _ { R } ( \mathcal { O } ) \right \rangle _ { \mathrm U ( N ) } = \frac { 1 } { \mathcal { Z } _ { \mathrm C S } } \int \prod _ { i = 1 } ^ { N } \frac { d x _ { i } } { 2 \pi } \, e ^ { - \frac { 1 } { 2 g _ { s } } x _ { i } ^ { 2 } } \det _ { 1 \leq i, j \leq N } e ^ { x _ { i } \tilde { \xi } _ { j } } \prod _ { i < j } ^ { N } \left ( 2 \sinh \frac { x _ { i } - x _ { j } } { 2 } \right ) \,,$$

with ˜ ξ i = λ i -i + N +1 2 . In this sense insertion of the Wilson loop operator corresponds to the external fields for the matrix integral.

This kind of interpretation can be possible even for the supergroup character average (2.26), at least as long as the partition λ , corresponding to the representation R , satisfies d λ ( ) = N . 5 All this means that characters of U( N N | ) are dual to expectation values of characteristic polynomials of another matrix model with some measure µ ,

$$\left \langle W _ { R } ( K ) \right \rangle _ { U ( N | N ) } = \left \langle \prod _ { i = 1 } ^ { N } \frac { \det ( \alpha _ { i } - M ) } { \det ( \beta _ { i } - M ) } \right \rangle _ { \mu } \,,$$

and where α , β i i are the Frobenius coordinates for the representation R . Then the Fay identity says that

$$\left \langle W _ { R } ( K ) \right \rangle _ { U ( N | N ) } = \det _ { 1 \leq i, j \leq N } \left \langle W _ { ( \alpha _ { i } | \beta _ { j } ) } ( K ) \right \rangle _ { U ( 1 | 1 ) }.$$

This is what we have checked in this article, especially as the determinantal formula for the unknot average shown in (2.28). This means that those supergroup character expectation values are Giambelli compatible [12].

## 7 Discussion

In this paper we have considered the supergroup character expectation value based on the ABJM matrix model as the supermatrix Chern-Simons theory, with emphasis on its connection to the knot invariant. We have explicitly computed the U(1 1) expectation value for |

5 When the representation does not satisfy this condition, there is no simple analogy between the matrix integral with the external fields and the Wilson loop average (2.34), because such an external field has to consist of N + N parameters for U( N N | ) theory. On the other hand, this analogy holds for arbitrary representations in the ordinary U( N ) Chern-Simons theory as (6.21). It is because even if the number of non-zero elements in the partition is less than N at the first place, it can be made N by the constant shift, since the average (2.31) is invariant under the shift of the partition

$$( \lambda _ { 1 }, \lambda _ { 2 }, \cdots, \lambda _ { N } ) \ \longrightarrow \ ( \lambda _ { 1 } + c, \lambda _ { 2 } + c, \cdots, \lambda _ { N } + c ) \,.$$

It is obvious that the U( N N | ) invariant (2.26) is not invariant under such a constant shift of the partition, since the corresponding external field is characterized by the Frobenius coordinates of the partition.

the unknot and torus knot matrix models. We have obtained the determinantal formula, where the U(1 1) average plays a role of a building block for U( | N N | ) theory, and shown the Rosso-Jones-type formula for the supergroup character average. We have also discussed the U( N M | ) theory, and found its determinantal formula which interpolates U( N N | ) and U( N ) theories. We have derived the spectral curve for the torus knot as the symplectic transform of the unknot curve, and by analyzing the saddle point equations for the torus knot matrix model itself. We have shown these two methods to obtain the spectral curve are consistent with each other. We have then commented on how to realize the knot invariant in topological string theory, and the underlying idea coming from random matrix theory.

Let us comment on some open issues to be investigated in the future. The most attractive one is to check whether the supergroup character average can be a knot invariant or not. We have certain evidence on this point. First of all, the U( N N | ) character expectation value contains U( N ) part as shown for the unknot Wilson loop (2.29). Furthermore, in the limit M → 0, the U( N M | ) matrix model (4.8) is explicitly reduced to the U( N ) torus knot matrix model, which is proven to be a topological invariant of the torus knot. This implies that, at least for the torus knot, the U( N M | ) average would play a role of the knot invariant which contains much information than the ordinary U( N ) invariant.

Secondly the ABJM matrix model is obtained from the lens space Chern-Simons theory through the analytic continuation. Since the knot invariant in the lens space is given as the character average with the corresponding matrix model, we can expect the character expectation value with the supermatrix model is also a knot invariant, at least up to some non-perturbative contributions. Moreover, the construction of the knot invariant shown in Sec. 6.1 seems natural from the viewpoint of topological string, and possibly applied to an arbitrary knot. These supporting facts suggest that we can obtain a knot invariant from ABJM theory. For this purpose, for example, it is interesting to see whether the torus knot supermatrix model can be directly obtained from ABJM theory. In the case of classical group theory, it is shown by [51] that the torus knot matrix model is obtained from N = 2 Chern-Simons theory on the ellipsoid-type squashed three-sphere S 3 b . This suggests that the supergroup torus knot character average is given by the half-BPS Wilson loop for the N = 6 ABJM theory on S 3 b . However the integral formula (3.10) is not obtained naively using the localization method for S 3 b , because the matter contribution cannot be written as a simple cosh function for such a case [52]. In addition, the mirror description of the torus knot model (3.10), presented in (A.3), is different from that on S 3 b . Although it is not yet obvious that such a way to obtain the torus knot invariant using the ellipsoid S 3 b is applicable to the present case, this discrepancy would reflect that ABJM theory itself is not topological, and we may need to consider, for example, topological twist to obtain the knot invariant for the Seifert manifold [2, 3, 4, 5]. We will return to this issue in the future.

If we can have the supergroup knot invariant, it is also interesting to provide another definition of the U( N N | ) knot invariant, for example, based on the Skein relation, which

enable us to compute the invariant for any knots in principle. It is naively expected that such a relation can be related to the supergroup WZW model. However, it is known that the U(1 1) WZW model describes the Skein relation for the Alexander-Conway polynomial [53, | 54, 55], while the torus knot average obtained in this paper (3.14) is not consistent with that. This reflects the fact that ABJM theory is not just Chern-Simons theory with the supergroup. We have to explore other possibility of the two-dimensional conformal field theory model, which shall live on the boundary of ABJM theory, if it exists. Furthermore, as the recent progress in knot theory, it is definitely interesting to study the volume conjecture and its generalization [56, 57, 58], the AJ conjecture [59], and also the knot homology [ ? , 60, 61, 62, 63] corresponding to the supergroup knot invariant. In order to discuss the volume conjecture, we have to deal with the hyperbolic knot and the volume of its complement in S 3 . Therefore our construction, which is so far available only for the unknot and torus knots, is not yet enough to study this conjecture. Also from this point of view, the definition of the knot invariant based on the Skein relation is highly desirable, as commented above. The AJ conjecture claims that the knot invariant satisfies some integrable equations, which are obtained through quantization of the A-polynomial. In other words, this implies that the knot invariant is associated with the corresponding τ -function. As pointed out in Sec. 6, the knot invariant is closely related to the matrix integral with the external fields, which can be seen as a certain τ -function. For example, the determinantal formula for the character expectation value given in this paper can be one of the supporting results for such a suggestive relation. The determinantal formula also gives an insight into the knot homology. Using the Jacobi identity for determinants, one can obtain some relations between the knot invariants for different rank groups. Such a relation could provide a natural differential on the knot homology.

## Acknowledgements

We would like to thank G. Borot and M. Mari˜o for fruitful discussions. We are also graten ful to A. Brini for carefully reading the manuscript and giving useful comments. BE thanks Centre de Recherches Math´ ematiques de Montr´al, the FQRNT grant from the Qu´bec gove e ernment, Piotr Suglyph[suppress]kowski l and the ERC starting grant Fields-Knots. The work of TK is supported in part by Keio Gijuku Academic Development Funds, MEXT-Supported Program for the Strategic Research Foundation at Private Universities 'Topological Science' (No. S1511006), and JSPS Grant-in-Aid for Scientific Research on Innovative Areas 'Topological Materials Science' (No. JP15H05855).

## A Mirror of the torus knot ABJM partition function

In this appendix we derive the mirror description of the torus knot ABJM partition function [25]. We first expand the determinants in (3.10) as summation over permutations

$$\mathcal { Z } _ { \text{ABJM} } ^ { ( P, Q ) } = \sum _ { \sigma, \sigma ^ { \prime } \in \mathfrak { S } _ { N } } ( - 1 ) ^ { \sigma + \sigma ^ { \prime } } \frac { 1 } { N! ^ { 2 } } \int _ { [ d x ] ^ { N } [ d y ] ^ { N } } \prod _ { i = 1 } ^ { N } \left ( 2 \cosh \frac { x _ { i } - y _ { \sigma ( i ) } } { 2 P } \, 2 \cosh \frac { x _ { i } - y _ { \sigma ^ { \prime } ( i ) } } { 2 Q } \right ) ^ { - 1 }.$$

Applying the formula (2.20), we have

$$\text{pying use formula } ( z _ { \sigma } \varnothing ), \, \text{we move} \\ \sum _ { \sigma, \sigma ^ { \prime } \in \mathfrak { S } _ { N } } ( - 1 ) ^ { \sigma + \sigma ^ { \prime } } \frac { 1 } { N! ^ { 2 } } \int _ { \frac { d ^ { N } x } { ( 2 \pi ) ^ { N } } } \frac { d ^ { N } y } { ( 2 \pi ) ^ { N } } \frac { d ^ { N } z } { ( 2 \pi ) ^ { N } } \, \frac { d ^ { N } w } { ( 2 \pi ) ^ { N } } \, \prod _ { i = 1 } ^ { N } \left ( \cosh z _ { i } \, \cosh w _ { i } \right ) ^ { - 1 } \\ \times \exp \left [ \frac { i k } { 4 P Q \pi } \sum _ { i = 1 } ^ { N } ( x _ { i } ^ { 2 } - y _ { i } ^ { 2 } ) + \frac { i } { \pi } \sum _ { i = 1 } ^ { N } \left ( x _ { i } \left ( \frac { z _ { i } } { P } + \frac { w _ { i } } { Q } \right ) - y _ { i } \left ( \frac { z _ { \sigma } - 1 ( i ) } { P } + \frac { w _ { \sigma ^ { \prime } - 1 ( i ) } } { Q } \right ) \right ) \right ] \\ = \sum _ { \sigma, \sigma ^ { \prime } \in \mathfrak { S } _ { N } } ( - 1 ) ^ { \sigma + \sigma ^ { \prime } } \frac { ( P Q ) ^ { N } } { N! ^ { 2 } \, k ^ { N } } \int _ { \frac { d ^ { N } z } { ( 2 \pi ) ^ { N } } } \frac { d ^ { N } w } { ( 2 \pi ) ^ { N } } \prod _ { i = 1 } ^ { N } \frac { \exp \left [ - \frac { 2 i } { k \pi } \left ( z _ { i } w _ { i } - z _ { \sigma } - 1 ( i ) w _ { \sigma ^ { \prime } - 1 ( i ) } \right ) \right ] } { \cosh z _ { i } \, \cosh w _ { i } } \,. \ \ ( A. 2 )$$

At this moment it is obvious that the partition function depends only on the composition of permutations σ · σ ′-1 . Thus, by fixing either of them σ ′ as the trivial permutation, we obtain

$$\text{ obtain} \\ \mathcal { Z } _ { \text{ABJM} } ^ { ( P Q ) } & = \sum _ { \sigma \in \mathfrak { S } _ { N } } ( - 1 ) ^ { N } \frac { ( P Q ) ^ { N } } { N! k ^ { N } } \int _ { \frac { d ^ { N } z } { ( 2 \pi ) ^ { N } } } \frac { d ^ { N } w } { ( 2 \pi ) ^ { N } } \prod _ { i = 1 } ^ { N } \frac { \exp \left [ \frac { - \frac { 2 i } { k \pi } } { \left ( z _ { i } - z _ { \sigma ( i ) } \right ) w _ { i } } \right ] } { \cosh z _ { i } \cosh w _ { i } } \\ & = \sum _ { \sigma \in \mathfrak { S } _ { N } } ( - 1 ) ^ { N } \frac { ( P Q ) ^ { N } } { N! k ^ { N } } \int _ { \frac { d ^ { N } z } { ( 2 \pi ) ^ { N } } } \prod _ { i = 1 } ^ { N } \left ( \cosh z _ { i } \cdot 2 \cosh \frac { z _ { i } - z _ { \sigma ( i ) } } { k } \right ) ^ { - 1 } \\ & = \frac { ( P Q ) ^ { N } } { N! ( 2 k ) ^ { N } } \int _ { \frac { d ^ { N } z } { ( 2 \pi ) ^ { N } } } \prod _ { i < j } ^ { N } \left ( \tanh \frac { z _ { i } - z _ { j } } { 2 k } \right ) ^ { 2 } \prod _ { i = 1 } ^ { N } \left ( 2 \cosh \frac { z _ { i } } { 2 } \right ) ^ { - 1 } \\ & = ( P Q ) ^ { N } \mathcal { Z } _ { \text{ABJM} } ^ { ( 1, 1 ) } \cdot \\ \text{This is the mirror expression of the partition function (3.10).  Especially for $k=1$, the mirror}$$

This is the mirror expression of the partition function (3.10). Especially for k = 1, the mirror theory turns out to be N = 4 SYM theory with a single fundamental and a single adjoint hypermultiplet. The dependence on the parameters ( P, Q ) becomes trivial in this mirror representation.

## References

- [1] E. Witten, 'Quantum field theory and the Jones polynomial,' Commun. Math. Phys. 121 (1989) 351-399.
- [2] R. Lawrence and L. Rozansky, 'Witten-Reshetikhin-Turaev Invariants of Seifert Manifolds,' Commun. Math. Phys. 205 (1999) 287-314.

| [3]   | C. Beasley, 'Localization for Wilson Loops in Chern-Simons Theory,' Adv. Theor. Math. Phys. 17 (2013) 1-240, arXiv:0911.2687 [hep-th] .                                                                   |
|-------|-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [4]   | J. K¨ll´n, a e 'Cohomological localization of Chern-Simons theory,' JHEP 1108 (2011) 008, arXiv:1104.5353 [hep-th] .                                                                                      |
| [5]   | A. Brini, B. Eynard, and M. Mari˜o, n 'Torus knots and mirror symmetry,' Ann. Henri Poincar´ e 13 (2012) 1873-1910, arXiv:1105.2012 [hep-th] .                                                            |
| [6]   | A. Kapustin, B. Willett, and I. Yaakov, 'Exact Results for Wilson Loops in Superconformal Chern-Simons Theories with Matter,' JHEP 1003 (2010) 089, arXiv:0909.4559 [hep-th] .                            |
| [7]   | O. Aharony, O. Bergman, D. L. Jafferis, and J. Maldacena, ' N = 6 superconformal Chern-Simons-matter theories, M2-branes and their gravity duals,' JHEP 0810 (2008) 091, arXiv:0806.1218 [hep-th] .       |
| [8]   | N. Drukker and D. Trancanelli, 'A Supermatrix model for N = 6 super Chern-Simons-matter theory,' JHEP 1002 (2010) 058, arXiv:0912.3006 [hep-th] .                                                         |
| [9]   | I. Bars, 'Supergroups and Their Representations,' Lectures Appl. Math. 21 (1983) 17.                                                                                                                      |
| [10]  | A. Berele and A. Regev, 'Hook Young-diagrams with applications to combinatorics and to representations of Lie superalgebras,' Adv. Math. 64 (1987) 118-175.                                               |
| [11]  | J. D. Fay, Theta Functions on Riemann Surfaces , vol. 352 of Lecture Notes in Mathematics . Springer Berlin Heidelberg, 1973.                                                                             |
| [12]  | A. Borodin, G. Olshanski, and E. Strahov, 'Giambelli compatible point processes,' Adv. Appl. Math. 37 (2006) 209-248, arXiv:math-ph/0505021 .                                                             |
| [13]  | M. Berg`re e and B. Eynard, 'Determinantal formulae and loop equations,' arXiv:0901.3273 [math-ph] .                                                                                                      |
| [14]  | B. Eynard and N. Orantin, 'Invariants of algebraic curves and topological expansion,' Commun. Num. Theor. Phys. 1 (2007) 347-452, arXiv:math-ph/0702045 [math-ph]                                         |
| [15]  | R. Dijkgraaf, H. Fuji, and M. Manabe, 'The Volume Conjecture, Perturbative Knot Invariants, and Recursion Relations for Topological Strings,' Nucl. Phys. B849 (2011) 166-211, arXiv:1010.4542 [hep-th] . |
| [16]  | G. Borot and B. Eynard, 'All-order asymptotics of hyperbolic knot invariants from non-perturbative topological recursion of A-polynomials,' Quantum Top. 6 (2015) 39-138, arXiv:1205.2261 [math-ph] .     |

| [17]   | S. Gukov and I. Saberi, 'Lectures on Knot Homology and Quantum Curves,' arXiv:1211.6075 [hep-th] .                                                        |
|--------|-----------------------------------------------------------------------------------------------------------------------------------------------------------|
| [18]   | R. Gopakumar and C. Vafa, 'On the gauge theory/geometry correspondence,' Adv. Theor. Math. Phys. 3 (1999) 1415-1443, arXiv:hep-th/9811131 [hep-th] .      |
| [19]   | H. Ooguri and C. Vafa, 'Knot invariants and topological strings,' Nucl. Phys. B577 (2000) 419-438, arXiv:hep-th/9912123 [hep-th] .                        |
| [20]   | M. Mari˜o n and P. Putrov, 'Exact Results in ABJM Theory from Topological Strings,' JHEP 1006 (2010) 011, arXiv:0912.3074 [hep-th] .                      |
| [21]   | E. Moens and J. Van der Jeugt, 'A Determinantal Formula for Supersymmetric Schur Polynomials,' J. Alg. Comb. 17 (2003) 283-307.                           |
| [22]   | Y. Hatsuda, M. Honda, S. Moriyama, and K. Okuyama, 'ABJM Wilson Loops in Arbitrary Representations,' JHEP 1310 (2013) 168, arXiv:1306.4297 [hep-th] .     |
| [23]   | Y. Dolivet and M. Tierz, 'Chern-Simons matrix models and Stieltjes-Wigert polynomials,' J. Math. Phys. 48 (2007) 023507, arXiv:hep-th/0609167 [hep-th] .  |
| [24]   | M. Rosso and V. Jones, 'On the invariants of torus knots derived from quantum groups,' J. Knot Theory Ramifications 2 (1993) 97-112.                      |
| [25]   | A. Kapustin, B. Willett, and I. Yaakov, 'Nonperturbative Tests of Three-Dimensional Dualities,' JHEP 1010 (2010) 013, arXiv:1003.5694 [hep-th] .          |
| [26]   | O. Aharony, O. Bergman, and D. L. Jafferis, 'Fractional M2-branes,' JHEP 0811 (2008) 043, arXiv:0807.4924 [hep-th] .                                      |
| [27]   | E. L. Basor and P. J. Forrester, 'Formulas for the Evaluation of Toeplitz Determinants with Rational Generating Functions,' Math. Nachr. 170 (1994) 5-18. |
| [28]   | M. Aganagic, A. Klemm, M. Mari˜o, n and C. Vafa, 'Matrix model as a mirror of Chern-Simons theory,' JHEP 0402 (2004) 010, hep-th/0211098 .                |
| [29]   | N. Halmagyi and V. Yasnov, 'The spectral curve of the lens space matrix model,' JHEP 0911 (2009) 104, arXiv:hep-th/0311117 .                              |
| [30]   | H. Jockers, A. Klemm, and M. Soroush, 'Torus Knots and the Topological Vertex,' Lett. Math. Phys. 104 (2014) 953-989, arXiv:1212.0321 [hep-th] .          |
| [31]   | S. Stevan, 'Torus Knots in Lens Spaces and Topological Strings,' Ann. Henri Poincar´ e 16 (2015) 1937-1967, arXiv:1308.5509 [hep-th] .                    |

| [32]   | G. Borot, B. Eynard, and N. Orantin, 'Abstract loop equations, topological recursion and new applications,' Commun. Num. Theor. Phys. 9 (2015) 51-187, arXiv:1303.5808 [math-ph] .                                                  |
|--------|-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [33]   | G. Borot, A. Guionnet, and K. K. Kozlowski, 'Large- N asymptotic expansion for mean field models with Coulomb gas interaction,' Int. Math. Res. Not. 2015 (2015) 10451-10524, arXiv:1312.6664 [math-ph] .                           |
| [34]   | M. Aganagic, A. Klemm, M. Marino, and C. Vafa, 'The Topological vertex,' Commun. Math. Phys. 254 (2005) 425-478, arXiv:hep-th/0305132 [hep-th] .                                                                                    |
| [35]   | A. Klemm, M. Mari˜o, n M. Schiereck, and M. Soroush, 'ABJM Wilson loops in the Fermi gas approach,' Z. Naturforsch. A68 (2013) 178-209, arXiv:1207.0611 [hep-th] .                                                                  |
| [36]   | M. Mari˜o, n 'Open string amplitudes and large order behavior in topological string theory,' JHEP 0803 (2008) 060, arXiv:hep-th/0612127 [hep-th] .                                                                                  |
| [37]   | V. Bouchard, A. Klemm, M. Mari˜o, n and S. Pasquetti, 'Remodeling the B-model,' Commun. Math. Phys. 287 (2009) 117-178, arXiv:0709.1453 [hep-th] .                                                                                  |
| [38]   | B. Eynard and N. Orantin, 'Computation of open Gromov-Witten invariants for toric Calabi-Yau 3-folds by topological recursion, a proof of the BKMP conjecture,' Commun. Math. Phys. 337 (2015) 483-567, arXiv:1205.1103 [math-ph] . |
| [39]   | M. Aganagic and C. Vafa, 'Mirror symmetry, D-branes and counting holomorphic discs,' arXiv:hep-th/0012041 [hep-th] .                                                                                                                |
| [40]   | M. Aganagic, A. Klemm, and C. Vafa, 'Disk instantons, mirror symmetry and the duality web,' Z. Naturforsch. A57 (2002) 1-28, arXiv:hep-th/0105045 [hep-th] .                                                                        |
| [41]   | M. Aganagic and C. Vafa, 'Large N Duality, Mirror Symmetry, and a Q-deformed A-polynomial for Knots,' arXiv:1204.4709 [hep-th] .                                                                                                    |
| [42]   | Y. Hatsuda, M. Mari˜o, n S. Moriyama, and K. Okuyama, 'Non-perturbative effects and the refined topological string,' JHEP 1409 (2014) 168, arXiv:1306.1734 [hep-th] .                                                               |
| [43]   | B. Eynard and N. Orantin, 'Topological expansion of mixed correlations in the Hermitian 2-matrix model and x - y symmetry of the F g invariants,' J. Phys. A41 (2008) 015203, arXiv:0705.0958 [math-ph] .                           |
| [44]   | B. Eynard and N. Orantin, 'About the x - y symmetry of the F g algebraic invariants,' arXiv:1311.4993 [math-ph] .                                                                                                                   |
| [45]   | T. Kimura, 'Note on a duality of topological branes,' PTEP 2014 (2014) 103B04, arXiv:1401.0956 [hep-th] .                                                                                                                           |

| [46]   | T. Kimura, 'Duality and integrability of a supermatrix model with an external source,' PTEP 2014 (2014) 123A01, arXiv:1410.0680 [math-ph] .                                                                                                     |
|--------|-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [47]   | P. Desrosiers and B. Eynard, 'Supermatrix models, loop equations, and duality,' J. Math. Phys. 51 (2010) 123304, arXiv:0911.1762 [math-ph] .                                                                                                    |
| [48]   | E. Br´zin e and S. Hikami, 'Characteristic Polynomials of Random Matrices,' Commun. Math. Phys. 214 (2000) 111-135, arXiv:math-ph/9910005 .                                                                                                     |
| [49]   | Harish-Chandra, 'Differential Operators on a Semisimple Lie Algebra,' Amer. J. Math. 79 (1957) 87-120.                                                                                                                                          |
| [50]   | C. Itzykson and J.-B. Zuber, 'The Planar Approximation. II,' J. Math. Phys. 21 (1980) 411-421.                                                                                                                                                  |
| [51]   | A. Tanaka, 'Comments on knotted 1/2 BPS Wilson loops,' JHEP 1207 (2012) 097, arXiv:1204.5975 [hep-th] .                                                                                                                                         |
| [52]   | N. Hama, K. Hosomichi, and S. Lee, 'SUSY Gauge Theories on Squashed Three-Spheres,' JHEP 1105 (2011) 014, arXiv:1102.4716 [hep-th] .                                                                                                            |
| [53]   | T. Deguchi and Y. Akutsu, 'Graded Solutions of the Yang-Baxter Relation and Link Polynomials,' J. Phys. A23 (1990) 1861-1876.                                                                                                                   |
| [54]   | L. H. Kauffman and H. Saleur, 'Free fermions and the Alexander-Conway polynomial,' Commun. Math. Phys. 141 (1991) 293-327.                                                                                                                      |
| [55]   | L. Rozansky and H. Saleur, ' S and T matrices for the super U(1,1) WZW model: Application to surgery and three manifolds invariants based on the Alexander-Conway polynomial,' Nucl. Phys. B389 (1993) 365-423, arXiv:hep-th/9203069 [hep-th] . |
| [56]   | R. M. Kashaev, 'The Hyperbolic volume of knots from quantum dilogarithm,' Lett. Math. Phys. 39 (1997) 269-275, arXiv:q-alg/9601025 [math.QA] .                                                                                                  |
| [57]   | H. Murakami and J. Murakami, 'The colored Jones polynomials and the simplicial volume of a knot,' Acta Math. 186 (2001) 85-104, arXiv:math/9905075 [math.GT]                                                                                    |
| [58]   | S. Gukov, 'Three-dimensional quantum gravity, Chern-Simons theory, and the A-polynomial,' Commun. Math. Phys. 255 (2005) 577-627, arXiv:hep-th/0306165 [hep-th] .                                                                               |
| [59]   | S. Garoufalidis, 'On the characteristic and deformation varieties of a knot,' Geom. Topol. Monogr. 7 (2004) 291-309, arXiv:math/0306230 [math.GT] .                                                                                             |
| [60]   | M. Khovanov and L. Rozansky, 'Matrix factorizations and link homology,' Fund. Math. 199 (2008) 1-91, arXiv:math/0401268 [math.QA] .                                                                                                             |

- [61] M. Khovanov and L. Rozansky, 'Matrix factorizations and link homology II,' Geom. Topol. 12 (2008) 1387-1425, arXiv:math/0505056 [math.QA] .
- [62] S. Gukov, A. S. Schwarz, and C. Vafa, 'Khovanov-Rozansky homology and topological strings,' Lett. Math. Phys. 74 (2005) 53-74, arXiv:hep-th/0412243 [hep-th] .
- [63] N. M. Dunfield, S. Gukov, and J. Rasmussen, 'The Superpolynomial for Knot Homologies,' Experimental Math. 15 (2006) 129-159, arXiv:math/0505662 [math.GT] .
Bertrand Eynard Taro Kimura
Towards U(N|M) knot invariant from ABJM theory
We study U(N|M) character expectation value with the supermatrix Chern-Simons
theory, known as the ABJM matrix model, with emphasis on its connection to the
knot invariant. This average just gives the half BPS circular Wilson loop
expectation value in ABJM theory, which shall correspond to the unknot
invariant. We derive the determinantal formula, which gives U(N|M) character
expectation values in terms of U(1|1) averages for a particular type of
character representations. This means that the U(1|1) character expectation
value is a building block for all the U(N|M) averages, and in particular, by an
appropriate limit, for the U(N) invariants. In addition to the original model,
we introduce another supermatrix model obtained through the symplectic
transform, which is motivated by the torus knot Chern-Simons matrix model. We
obtain the Rosso-Jones-type formula and the spectral curve for this case.