Complete factorization in minimal N=4 Chern-Simons-matter theory

We investigate an N=4 U(N)_k x U(N+M)_{-k} Chern-Simons theory coupled to one bifundamental hypermultiplet by employing its partition function, which is given by 2N+M dimensional integration via localization. Surprisingly, by performing the integration explicitly we find that the partition function completely factorizes into that of the pure Chern-Simons theory for two gauge groups and an analogous contribution for the bifundamental hypermultiplet. Using the factorized form of the partition function we argue the level/rank duality, which is also expected from the Hanany-Witten transition in the type IIB brane realization. We also present the all order 't Hooft expansion of the partition function and comment on the connection to the higher-spin theory.


Introduction
It has been more than a decade since a three dimensional maximally superconformal Chern-Simons-matter theory has been discovered. The construction of such a theory was highly important since they are expected to describe multiple M2-branes in M-theory [1,2,3]. The first example was obtained by generalizing the usual Lie algebra to the three bracket to describe "Chan-Paton" degrees of freedom associated with the M2-branes [4,5,6,7] (see also [8]), which turned out to describe at most two M2-branes in a countable case. Subsequently superconformal Chern-Simons theories with eight supercharges were constructed as linear quiver gauge theories [9] associated with type IIB brane configuration [10], without relying on the three bracket. Then it was insightfully pointed out that by changing the minimal linear quiver gauge theory to the circular one the system achieves supersymmetry enhancement and describes interaction of arbitrary number of membranes [11]. 1 This opened up the study of superconformal Chern-Simons-matter theories from AdS/CFT correspondence [17,18,19]. 2 One of the non-trivial checks for the AdS/CFT correspondence was given by the exact computation of the partition function with the help of the supersymmetric localization technique [24,25,26]. The localization technique reduces the path integral for the partition function to a finite dimensional matrix integration [27,28,29]. These matrix models can be evaluated by the saddle point approximation in the large N limit [30,31], and found to be consistent with the classical approximation in the gravity side. It also provides non-trivial checks for various field theory dualities as some simple integration identities [32,33,34,35]. Moreover, for some special class of N = 4 quiver superconformal Chern-Simons theories the matrix model can be regarded as the partition function of N free fermions in one dimensional space [36], which enables us to evaluate the all order 1/N corrections explicitly [37,38,39]. This lead us to the discovery of the coincidence between the large N instanton effects in the partition function and a non-perturbative completion of the topological string free energy [40], which was not clear a priori from neither the construction of the theory nor the gravity dual.
In this situation we revisit the superconformal Chern-Simons theories characterized by linear quiver diagrams, which are also called Gaiotto-Witten theories [9], from its exact partition function determined by the supersymmetric localization. Particularly in this paper we intensively study the partition function of the minimal N = 4 Chern-Simons theory with U(N) k × U(N + M) −k gauge group interacting with a bifundamental hypermultiplet. Interestingly, we find that the partition function can be decomposed into a product of the pure Chern-Simons partition functions for the two gauge nodes and an analogous factor for the bifundamental representation in U(N)×U(N +M). We shall refer to this decomposition as complete factorization.
Our result allows us to argue the level/rank duality of the U(N) k ×U(N +M) −k linear quiver theory, whose evidence was provided in [41] in the vector model limit. We can provide the dictionary of particular observables under the duality by turning on the FI mass deformations. Once we obtain the partition function in a closed form expression we can also evaluate it in various limits directly. In particular, written in the form of complete factorization, the all order 't Hooft expansion for the free energy can be obtained easily. Since the vector model limit corresponds to a boundary of the parameter space of the two 't Hooft couplings, we can discuss the relation of the U(N) k × U(N + M) −k Gaiotto-Witten theory to the dual higher spin theory, supergravity and their interpolation.
The rest of this paper is organized as follows. In the next section we review the basic properties of the U(N) k ×U(N +M) −k linear quiver theory. In §3, starting from the localization formula, we write down the partition function of this theory with mass deformation in a closed form expression for general k, N and M. We observe two novel structures of the partition function: the complete factorization (3.12) and the singularity in the massless limit (see §3.2). The complete factorization enables us various subsequent analysis. In §4 we argue the level/rank duality from the symmetry of the partition function, which, regarding the mass parameter as fugacity, also provides the dictionary of some observables under the duality. In §5 we evaluate the free energy in the 't Hooft expansion N, M, k → ∞ with λ 1 = N/k and λ 2 = (N + M)/k kept finite. We also comment on the vector model limit, which is realized as λ 1 = 0. §6 is devoted for discussion and comments on further directions. This paper also contains two appendices in order to help the derivation of the complete factorization. §A is a collection of determinant formulas used in §3. In §B we explain a part of the computation (3.9) in detail.

Minimal N = 4 Chern-Simons theory
In this paper we study the minimal N = 4 Chern-Simons theory, that is, U(N) k × U(N + M) −k Chern-Simons theory interacting with one bifundamental hypermultiplet. This theory was first constructed in [9] as the simplest example in a class of N = 4 linear quiver Chern-Simons matter theory. This theory has SO(4) R-symmetry and possesses the parity invariance accompanied with the exchange of the two gauge fields when the ranks of the gauge groups are the same. For general levels and ranks, the system has the generalized parity symmetry, under which the levels as well as the ranks are exchanged. This system is realized by the following type IIB brane configuration. (2.1) The D3-branes are extending in the x 6 -direction, ending on the two NS5-branes extending in the x 7 x 8 x 9 -directions. Between the NS5-branes the D3-branes are also intersecting with a (1, k) 5-brane which is extending in the x i x i+4 -planes (i = 3, 4, 5) with angle arctan k from each x i+4 axis. In the figure we have omitted x 0 x 1 x 2 -directions, where all the branes are extending as well.
It is known that this system admits the mass term keeping N = 4 supersymmetry as well as SO(4) R-symmetry [13]. In this class of N = 4 theories the mass term is equivalent to FI term. In what follows we study the partition function of this system including the N = 4 mass term induced by FI term. The generalization including the mass term is computationally not difficult but conceptually important. Indeed, in §4 we find a non-trivial coincidence of a particular pair of the partition functions by turning on the continuous deformation parameter as well. By differentiating the partition function by the deformation parameter in arbitrary times, we obtain the coincidence for many observables.

Exact partition function
In this section we study the S 3 partition function of the N = 4 U(N) k × U(N + M) −k Gaiotto-Witten theory with FI deformation. With the help of the supersymmetric localization we can reduce the partition function down to the following 2N + M dimensional integration [32,28,29] It is easy to see that this partition function is non-negative when M = 0, ζ 1 = −ζ 2 and satisfies the generalized parity, Below we first obtain the closed form expression for Z ζ 1 ,ζ 2 (k, N, N + M) for general k, N, and M by performing the integrations explicitly. Interestingly, the final expression (3.12) have completely the same structure as the one-loop determinant factors in the original integrand (3.1). That is, the partition function "factorizes" into those of the pure Chern-Simons gauge fields and the contributions from the bifundamental matter fields. We also argue the massless limit ζ 1 + ζ 2 → 0 of the partition function, where we encounter divergence if the ranks are greater than some critical value determined by k.
Note that in the matrix model (3.1) the one-loop determinant factor is not bounded in a generic direction at infinity for general k,N,M, which is in contrast to the case of ABJ(M) matrix model. Hence the convergence of (3.1) is not trivial. For the theories without Chern-Simons couplings, it was pointed out [32] that the convergence property of the partition function gives the same classification as that obtained by using monopole operators [9]. It would be interesting to examine the convergence of (3.1) rigorously, which will provide a generalization of such classification for the theories including Chern-Simons terms. In this section, however, we simply assume the convergence and perform the integrations formally to obtain the exact expression for the partition function Z ζ 1 ,ζ 2 (k, N, N +M). Nevertheless the convergence property may be partly reflected in the pole structure of Z ζ 1 ,ζ 2 (k, N, N + M).
Even if the partition function is divergent in some parameter regime, we may be able to extract a physical quantity with the help of an appropriate regularization, as is usual in a general quantum field theory. For example, we can follow the prescription given in [32]. That is, we can start with the matrix model Z ζ 1 ,ζ 2 (k 1 , k 2 , N, N + M) for the partition function of the U(N) k 1 × U(N + M) k 2 theory with two independent levels k 1 and k 2 , which is absolutely convergent when Im[k 1 ], Im[k 2 ] > 0. Then the desired partition function Z ζ 1 ,ζ 2 (k, N, N + M) can be defined as lim k 1 →k,k 2 →−k Z ζ 1 ,ζ 2 (k 1 , k 2 , N, N + M), the analytic continuation of Z ζ 1 ,ζ 2 (k 1 , k 2 , N, N + M) as a function of k 1 and k 2 . Though it is not trivial whether this scheme results in the same expression as (3.12), determining the explicit expression of Z ζ 1 ,ζ 2 (k 1 , k 2 , N, N + M) is beyond the scope of this paper hence we shall leave it as a future work.

Complete factorization
In order to perform the integration (3.1) first we rewrite the one-loop determinant by using the Cauchy-Vandermonde determinant formula (A.1) (with appropriate substitutions of x a = e µa , etc.) as follows where we have used different characters for each type of indices to encode their ranges: a, b, · · · = 1, 2, · · · , N; ℓ, m, · · · = 1, 2, · · · , M; i, j, · · · = 1, 2, · · · , N + M. By using these the partition function can be rewritten as Then we can combine the determinants by using the formula (A.2), which allows us to move the integration into each matrix element inside a whole determinant. We finally obtain Here the matrix elements of each block are given as which are directly computed as follows. The second block is easily computed by the gaussian integration as To compute the first block, we change the integration variables to µ ± = i(µ ± ν) (3.8) Thanks to the fact that the Chern-Simons levels add up to zero, the µ + integration just gives the delta function and we can easily perform the µ − integration as Substituting these two expressions to (3.5) we obtain Here we have factored out some phases in (3.9) and (3.7), which depend only on one of the row/column indices and hence can be interpreted as the multiplication of a diagonal matrix.
Notice that the determinant in the last expression (3.10) is of the form of the right-hand side of (3.3) with the following identification of the parameters 3 which allows us to decompose the determinant into the product of hyperbolic functions again.
Collecting the overall factors together, we finally find a surprisingly simple expression for the partition function Here P is a phase factor given by and Z (CS) ζ (k, N) is the partition function of pure Chern-Simons theory with FI deformation That is, the partition function factorizes into the contributions from the two vector multiplets and the one from the bifundamental hypermultiplet completely. This is non-trivial since in the original theory the bifundamental hypermultiplet couples to both vector multiplets nontrivially. Note that, although we have introduced two FI parameters ζ 1 , ζ 2 , the "hypermultiplet contribution" Z (mat) ζ (k, N, N + M) depends only on a single FI parameter ζ with ζ = ζ 1 + ζ 2 . Indeed this dependence corresponds to the mass of the hypermultiplet (ζ 1 + ζ 2 )/k, which also motivate us to call the ratio Z ζ 1 ,ζ 2 (k, N, N + M)/(P Z (CS) . It is not difficult to check that this expression of partition function is non-negative when M = 0, ζ 1 = −ζ 2 and satisfies the generalized parity (3.2).
For later convenience we present a different expression for the matter partition function: In the above computation, we have naively used the delta function (3.9) to perform the integrations in M 1,a,j (3.6). Since the kernels of the delta function are not located on the real axis, however, this would require a more careful argument on the deformation of the integration contours, which we explain in §B in detail.
Unfortunately the explanation in §B works only in the regime N + M/2 − 1 < |k|/2. 4 In the following sections, however, we shall tentatively assume that our result (3.15) is valid beyond the region N + M/2 − 1 < |k|/2 by the analytic continuation in k, and argue the level/rank duality relation in §4.1 from the partition function. Performing the integration beyond this region is left as a future work.

Poles in massless limit
Lastly, let us consider the massless limit ζ 1 + ζ 2 → 0 of the partition function. In this limit, the partition function diverges for sufficiently large values of ranks N, M due to the cosine factors in the denominator of Z Let us clarify the concrete upper bound for the ranks where the partition function diverges at the saturation of the bound. The divergence occurs when a pair of indices (a, j) satisfies the following condition for some m ∈ Z. From this expression we immediately find that there are no divergence for odd k + N 2 − N 1 . In the case of even k + N 2 − N 1 , to avoid the divergence the range of the left-hand side of (3.17) has to be shorter than |k|. As the left-hand side of (3.17) is maximized at (a, j) = (N 1 , 1) and minimized at (a, j) = (1, N 2 ), we can easily identify this condition as the following inequality Note that the upper bound N + M/2 − 1 = |k|/2 is below the zero of the pure Chern-Simons partition functions, where the brane configuration (2.1) violates the s-rule [10] and hence breaks supersymmetry. Therefore it seems that the partition function can actually be singular inside the regime where the matrix model (3.1) obtained by supersymmetric localization physically makes sense. 5 We hope to clarify the interpretation for this divergence in future work [47].

Level/rank duality
In this section we argue the level/rank duality of the N = 4 U(N) k × U(N + M) −k superconformal Chern-Simons theory. We provide a strong evidence for this from the partition function exactly computed as (3.12) by verifying its invariance under the level/rank duality transformation, and from the type IIB brane configuration by moving 5-branes taking into account the Hanany-Witten effect [10].

From partition function
Let us show the level/rank duality from the partition function. Since our partition function is completely factorized (3.12), it is sufficient to show the duality for each factor. We shall ignore the phase factor P independent of the FI parameters since it is irrelevant for the observables generated by ∂ m ζ 1 ∂ n ζ 2 log Z ζ 1 ,ζ 2 (k, N, N + M). It is famous that the pure Chern-Simons partition function enjoys the level/rank duality with a phase factor Q = e πi(k−N) 2 2 independent of ζ. We can show that the matter contribution Z (mat) ζ (k, N 1 , N 2 ) also satisfy a similar self-dual relation with a ζ-independent phase factor R = (−1) . This relation follows from the multiplication formula Indeed, by filling one product by the complement N j=1 = k j=1 / k j=N 2 +1 we can apply this formula to the product with full range k j=1 , which generates a factor with a trivial dependence on the other index a. Therefore, up to overall phases we can replace N 2 → k − N 2 by moving the product from the denominator to the numerator. Since we have two products, after the operation twice the product of cosine factors finally come back to the denominator together with the replacement where we have renamed the indices as i = k − j + 1 in the second line and b = k − a + 1 in the third line respectively.
Substituting these results into the factorized expression for the partition function (3.12), we obtain the following self-dual relation up to a ζ-independent phase. We shall call this the level/rank duality for the N = 4 U(N) k × U(N + M) −k theory. Now let us consider the massless limit and argue the relation among the observables, regarding both sides of (4.5) as the generating function. First of all we notice that the prefactor in the massless limit 2 sin π((k +M)/2) k−2N −M is zero or infinity for even k +M. This is indeed consistent with the observation in §3.2, as either Z ζ 1 ,ζ 2 (k, N, N + M) or Z ζ 2 ,ζ 1 (k, k − N − M, k − N) is singular. Here we shall focus on the case with odd k + M. Since sin π((k + M)/2 + i(ζ 1 + ζ 2 )) ∝ cosh π(ζ 1 + ζ 2 ) up to a phase independent of ζ 1 , ζ 2 in this case, we obtain (m + n ≥ 1) . (4.6)

From brane realization
The level/rank duality of this system can be understood from the type IIB brane realization as shown in other N = 3 or N = 4 superconformal Chern-Simons matter theories [48,33].
For simplicity, let us start with the massless case, whose brane configuration is given as (2.1). To consider the Hanany-Witten transition involving (1, k)5-brane, we shall follow [48]. 6 First we decompose the (1, k) 5-brane into the NS5-brane extending in the x 7 x 8 x 9 -directions and k D5-branes extending in the x 3 x 4 x 5 -directions, where the D5-branes are slightly on the left of the NS5-brane as and move the rightmost NS5-brane to the leftmost. When the 5-branes are exchanged, the number of D3-branes between them is changed so that the linking number associated with each 5-brane is unchanged [10]. Here the linking number for an NS5-brane is defined by ] is the number of all D5-branes to the right[left] of the NS5-brane and r D3 [l D3 ] is the number of D3-branes emanating from the NS5-brane to the right[left], and similarly the linking number for a D5-brane is defined by where R NS5 [L NS5 ] is the number of all NS5-branes to the right[left] of the D5-brane and r D3 [l D3 ] is the number of D3-branes emanating from the D5-brane to the right [left]. After the rightmost NS5-brane moves to the leftmost, the brane configuration becomes Subsequently we move the NS5-brane at the 1st slot to the rightmost, which results in the following brane configuration The theory realized on the D3-branes of this brane configuration is Next let us include the FI terms for this duality argument. Following [49] we realize the FI term for a gauge group realized on a stack of D3-branes ending with two 5-branes by the difference of the positions of the 5-branes in a transverse direction. Then ζ 1 , ζ 2 are realized by the difference of the positions of the NS5 1 -brane and NS5 2 -brane, NS5 2 -brane and NS5 3 -brane, respectively. After the above interchange of 5-branes, the two FI terms are clearly exchanged with extra minus, which is unimportant in the partition function. This completes our argument of level/rank duality from the brane configuration, up to the overall factor e πi(−ζ 2 1 +ζ 2 2 ) (2 sin π((k+ M)/2 + i(ζ 1 + ζ 2 ))) k−2N −M in (4.5).

Free energy in the 't Hooft expansion
In this section we consider the 1/N expansion of the free energy of the N = 4 U(N) k × U(N + M) −k theory. To obtain a non-trivial result, we have to take the 't Hooft limit N, M, k → ∞ with λ 1 = N/k and λ 2 = (N + M)/k kept finite, where λ 1 and λ 2 are bounded as 0 ≤ λ i ≤ 1 due to the s-rule. As the 't Hooft expansion for the pure Chern-Simons free energy was already studied in full detail (see e.g. [50,51]), let us focus on Z (mat) ζ 1 +ζ 2 (k, N, N + M) (3.15), which is the remaining part in the factorization structure (3.12) and is also of our interest in the vector model limit M, k → ∞ with M/k and N kept finite.
For simplicity below we consider the massless case ζ 1 = ζ 2 = 0 and neglect the phases, where the summations coming from the products can be rearranged as (see (3.16)) Then the large k expansion of F (mat) (k, N, N + M) can be evaluated in the 't Hooft limit by using the standard technique employed for the pure Chern-Simons partition function [50,51], which will take the following form of genus expansion with each coefficients F (mat) g (λ 1 , λ 2 ) given in (5.5) and (5.9) below.
Below we first consider the planar limit, F (mat) 0 (λ 1 , λ 2 ). Subsequently we determine the all order correction in 1/k in an explicit form, where we also assume M ∈ 2N for a technical reason. Finally we further comment on the vector model limit, which corresponds to the limit λ 1 → 0 of the planar free energy in our notation.

Planar approximation
In the strict limit of k, N, M → ∞ we can approximate the summation by the integration over a segment: Here the integrations can be performed explicitly by expanding log 2 cos z formally as , (5.4) which are concisely expressed with the trilogarithm. By using the notation of (5.2) we obtain Note that the leading free energy is an odd function in terms of λ 1 .
In general, the replacement K α=1 f (α/K) corrections in the integrand. Hence the evaluation above is valid only for the leading part in the large k, N, M limit.

All order 't Hooft expansion
Let us compute the 1/k corrections by reorganizing the summations in (5.1) into an appropriate infinite series. Changing the order of summation, we obtain with B ℓ the Bernoulli numbers defined by x/(e x − 1) = ∞ ℓ=0 B ℓ /ℓ!x ℓ . Thus we can replace the original finite summation with an infinite series, where each term in summation has definite powers of k, N, M. Since N = kλ 1 , M = k(λ 2 − λ 1 ), this directly provides the 1/k expansion. Collecting the coefficients together, we finally obtain where g runs all positive integers. Note that although in the beginning there exist contributions with g half integer, they totally vanish due to the fact that B 2n+1 = 0 for n ≥ 1. The planar part F (mat) 0 (λ 1 , λ 2 ) is consistent with the result (5.5) which we have obtained by replacing the summations with one dimensional integrations.

Comments on higher Spin/Supergravity crossover
Now we comment on the vector model limit, where k, M → ∞ with N and λ 2 kept finite. This limit can be achieved by the strict 't Hooft limit with λ 1 → 0.
The large M scaling behavior of the free energy in the vector model limit can be deduced from the planar result in the following way. Suppose the planar part of the free energy is expanded in terms of λ 1 around λ 1 ∼ 0 as Then the free energy scales in the vector model limit as This is because by using λ 1 = N/k, this can be written as In our case, from (5.5), the small λ 1 expansion of F (mat) (k, N, N + M) starts with Hence F (mat) (k, N, N + M) is linear in M in the vector model limit, which is the expected result.
It is expected that the N = 4 Chern-Simons matter theories exhibit higher spin symmetry in the vector model limit [52,53] (see also [54,55]), where higher spin currents transform as the adjoint representation of U(N). In fact, there is a conjecture that the minimal N = 4 theory investigated in this paper corresponds to a parity-violated supersymmetric higher spin Vasiliev theory on AdS 4 [52]. Here the higher spin fields have the U(N) color index, where the gauge interaction is governed by the bulk 't Hooft coupling (N + M)/N. It was argued that at the weak 't Hooft coupling region (N ∼ ∞) the theory is well described by a certain colored supergravity, while at the strong bulk 't Hooft coupling (N ∼ 1) the bulk color interaction is confined so that the theory is described by the higher spin Vasiliev theory. In the case of ABJ theory [56] this transition was shown to occur in a smooth fashion [57]. We have showed that this is also the case in the minimal Gaiotto-Witten theory.

Discussion
In this paper we have computed the partition function of the N = 4 U(N) k ×U(N +M) −k linear quiver superconformal Chern-Simons theory on three sphere. By performing the 2N + M integrations in the localization formula (3.1) explicitly, we have obtained the closed form expression for the partition function (3.12). Interestingly, we have found that the factor Z (mat) ζ 1 +ζ 2 (k, N, N +M) in (3.12), which is a priori mere the remainder after the two pure Chern-Simons partition function, takes completely the same form as the one-loop determinant of the bifundamental hypermultiplet in the localization formula (3.1). This motivates us to refer to our result as "complete factorization". The closed form expression has also enabled us to show the invariance of the partition function under the level/rank duality, which was also confirmed from the Hanany-Witten transition in the type IIB brane configuration, and to compute the all order 't Hooft explanation of the free energy.
To understand the implication of "factorization", it is important to assign more concrete interpretation to Z (mat) . One tempting candidate may be to regard the whole partition function as the pure Chern-Simons partition function for the supergroup U(N|N + M). In fact, the partition function of pure Chern-Simons theory appears as a modular transformation matrix of the characters in the corresponding Affine Lie algebra associated with the gauge group [58]. Our simple result of the matter partition function (3.12) may have such a purely algebraic and representational origin. It would be interesting to elucidate the origin by generalizing modular transformation matrices of general Affine Lie algebras classified in [59] including matter contribution.
Although our result (3.12)-(3.15) may be valid only in the region N + M/2 − 1 < |k|/2, the result is meaningful beyond this region after the analytic continuation in k, which is indeed necessary to argue the level-rank duality. Once we define the partition function by the analytic continuation, we find that the partition function has poles at (3.17). This singular behavior is opposite from what happens in the violation of the bounds for the s-rule: the vanishing of the partition function which one might associate with the supersymmetry breaking in analogy of the argument with the Witten index on S 2 × S 1 . It will be interesting to provide a physical interpretation for these poles. 7 An analogous divergence is indeed known for the N = 4 U(N) gauge theory (without Chern-Simons term) coupled with N f matter multiplets, where the S 3 partition function diverges for N f /2 < N c ≤ N f [33,60]. In this case the theory is "bad" for N f /2 < N c ≤ N f and can be dual to a "good" theory with decoupled massless hypermultiplets. The contribution from the decoupled sector appears as a prefactor in the duality relation, which explains the divergence in the bad theory. The divergence in our theory as well as the relation under the level/rank exchange (4.5) might be interpreted in the same manner. To clarify this idea it would be important to identify the unitarity violating monopole operators or show their existence from e.g. the computation of the superconformal index of the U(N) k × U(N + M) −k linear quiver theory.
It is also interesting to compare our results with the scattering amplitudes in quantum field theories, where the analytic continuation is necessary as well to see their crossing symmetry. Indeed, a similar kind of prefactor as (4.5) appears in the crossing symmetry relation of scattering amplitude in Chern-Simons vector theory [61,62,63]. That prefactor is important not only for crossing symmetry but also for unitarity in scattering amplitudes. After noticing the similarity, it would be interesting to interpret the poles (3.17) in analogy to the bound states in scattering amplitude.
In [64] the free energy of the pure Chern-Simons theory was found to coincide with the topological string on the resolved conifold in the 't Hooft expansion. It would be fascinating to provide a topological string interpretation to (5.9) as well. The obstruction of this program will be in the fermionic nature of the matter partition function (3.12) in contrast to the bosonic nature of that of pure Chern-Simons theory. This difference is usually crucial in the free energy, which can be seen for example in cancellation of the contribution of high energy modes. Reproducing the matter partition function from the topological string may require some new topological object.
The minimal Gaiotto-Witten theory we have considered can be obtained from the ABJM theory by removing one bifundamental hypermultiplet and cutting open the circular quiver. Indeed, our computation in §3 for M = 0 can be reorganized as the "open version" of the Fermi gas formalism [36], where we compute the matrix elements over only N distinctive oneparticle states instead of the trace over whole Hilbert space of N fermions. Once we notice this viewpoint we can perform similar computation also for the linear quivers arising from the brane setup (2.1) with additional (1, k ′ )5-branes as well, where k ′ is not necessarily equal to k. We shall report the detailed analysis for these generalizations, together with the rank deformations and the addition of fundamental matter fields, in a separated work [65].
We hope to report any progress on these issues in near future.

A Determinant Formulas
• Cauchy-Vandermonde determinant formula: • Convolution of determinants (formula (A.1) in [66]) B Comment on computation of M 1,a,j In this section we revisit the computation of the matrix element M 1,a,j (3.9). In the computation in §3.1 we have formally applied the following integration formula to the µ + -integration to produce the delta function Though the choice of the integration contour s ∈ (−∞, ∞) together with the reality of t is crucial in this formula, these assumptions are not satisfied in the original integration (3.6).
Below we explain the computation in more detail. After a trivial change of the integration variables (µ, ν) → (x, y) = (µ + ν, µ − ν), (3.6) is written as Starting from (B.2) we shall follow the following two different ways: deform of integration contour so that the formula (B.1) is applicable ( §B.1) and change the order of integration so that the computation does not involve delta functions ( §B.2). Although giving two explanations would be redundant, we hope they are helpful in future for a more rigorous justification or generalizations of the model. Also note that both of the two computations assume |Re[B]| < 1/2 as a regularization, whose justification is a subtle issue. 8

B.1 Deformation of integration contour
In §3.1 we have naively applied the integration formula (B.1) to the x-integration in the expression (B.2). Here let us look this step carefully. In order to apply the formula (B.1), we have to deform the contour of y-integration so that ky/4π − iA ∈ R. For this purpose we denote y as y = z + 4πiA/k and decompose the z-integration as

B.2 Computation without delta function
We can also compute M 1,a,j , without explicit appearance of a delta function, with the help of the following Fourier transformation formula where z = kx/2 + πiB. Here we have deformed the integral contour so as to pass through the origin without picking up any pole due to the assumption |Re[B]| < 1/2. Now we encounter the same constraint |Re[A]| < |k|/4, or (B.7), which is required here for the convergence of z-integration. When the condition (B.7) is satisfied we can apply the formula again (B.8) to obtain the same expression for M 1,a,j as (B.6).