format: text

text:
# S-duality of the Leigh-Strassler Deformation via Matrix Models

## 1 Introduction

It has been known for some time that the $`𝒩=4`$ supersymmetric Yang-Mills theory in four dimensions possesses exactly marginal deformations. Leigh and Strassler  identified two $`𝒩=1`$ SUSY-preserving exactly marginal directions other than the $`𝒩=4`$ gauge coupling itself, giving rise to a three (complex) dimensional renormalization-group-fixed manifold of $`𝒩=1`$ superconformal theories (SCFTs) of which the $`𝒩=4`$ fixed line is a subset. Given that strong and weak coupling on the $`𝒩=4`$ fixed line are related by the action of $`\text{SL}(2,)`$ on $`\tau `$, the gauge coupling of the $`𝒩=4`$ theory, it is then natural to ask whether this duality extends in some non-trivial way over the entire $`𝒩=1`$ fixed manifold.

In this article we answer this question for a particular Leigh-Strassler deformation of the $`𝒩=4`$ theory: the so-called “$`q`$-deformation”  111One such deformation has been analysed recently using matrix models in . . By studying a mass deformation of this Leigh-Strassler SCFT using a recently proposed connection between $`𝒩=1`$ theories and matrix models , we demonstrate that Montonen-Olive electric-magnetic duality, extended to the full $`\text{SL}(2,)`$ modular group, extends non-trivially to this space of $`𝒩=1`$ superconformal theories accompanied by a well-defined action on the marginal parameters. Not only will we uncover the duality group action in these theories, but we will also provide rather powerful checks on the applicability of the matrix model approach towards solving these $`𝒩=1`$ SUSY theories. (Recent checks of the matrix model proposal for other $`𝒩=1`$ models have been performed in ). In particular we will test the superpotential and eigenvalue spectrum obtained from the matrix model against classical field theory results and find nontrivial agreement between the two. Furthermore, rather remarkably the matrix model results for the vacuum expectation values of the effective superpotential can be used to directly infer an exact, dynamical quantum superpotential for the mass deformation of the Leigh-Strassler SCFT compactified on $`^3\times S^1`$. This superpotential written in terms of the effective fields of the 3-dimensional theory encodes the entire vacuum structure of the four-dimensional theory and satisfies extremely non-trivial checks to be discussed below. It turns out to be a natural deformation of the elliptic superpotential of  for the mass-deformed $`𝒩=4`$, or $`𝒩=1^{}`$, theory.

Relevant deformations of SCFTs can tell us a great deal about the CFTs themselves; in particular, as is well-known in the context of $`𝒩=4`$ SUSY Yang-Mills and the so-called $`𝒩=2`$ elliptic quiver models , the duality properties of the parent theories are inherited by the perturbed theories and manifest themselves as modular properties of various exactly calculable observables in the perturbed theories. In fact, what happens in these known examples is that the holomorphic observables in one vacuum get mapped into the corresponding observables in a different vacuum (realized in a different phase) under the action of $`\text{SL}(2,)`$ with well-defined modular weights. The holomorphic quantities in question can often be computed exactly and reflect the duality symmetries of the parent theory in this simple fashion.

The lesson to be drawn from this limited set of examples is that the consequences of dualities in these SCFTs are readily visible and more importantly, are readily calculable for certain relevant deformations of these SCFTs. With this in mind, rather than focusing on the SCFT itself we investigate a special mass perturbation of the Leigh-Strassler fixed points. This theory is the $`\text{SU}(N)`$ gauge theory with $`𝒩=1`$ SUSY and three adjoint-valued chiral multiplets $`\widehat{\mathrm{\Phi }}^+`$, $`\widehat{\mathrm{\Phi }}^{}`$ and $`\widehat{\mathrm{\Phi }}`$; the same matter content as the $`𝒩=4`$ theory, but with the following classical superpotential222We choose a normalization in which the kinetic terms of the chiral multiplets do not have factors of $`1/g_{YM}^2`$ in front. In the following, hatted quantities refer to $`\text{SU}(N)`$ to distinguish them from $`\text{U}(N)`$.

$$W_{\mathrm{cl}}=\mathrm{Tr}\left(i\lambda \widehat{\mathrm{\Phi }}[\widehat{\mathrm{\Phi }}^+,\widehat{\mathrm{\Phi }}^{}]_\beta +M\widehat{\mathrm{\Phi }}^+\widehat{\mathrm{\Phi }}^{}+\mu \widehat{\mathrm{\Phi }}^2\right),$$
(1)
where we have defined the $`q`$-commutator

$$[\widehat{\mathrm{\Phi }}^+,\widehat{\mathrm{\Phi }}^{}]_\beta \widehat{\mathrm{\Phi }}^+\widehat{\mathrm{\Phi }}^{}e^{i\beta /2}\widehat{\mathrm{\Phi }}^{}\widehat{\mathrm{\Phi }}^+e^{i\beta /2}$$
(2)
and where $`\lambda `$ and $`\beta `$ are complex bare couplings. We also introduce the complex bare gauge coupling of this theory $`\tau 4\pi i/g_{YM}^2+\theta /2\pi `$. This superpotential represents a deformation away from the $`𝒩=4`$ point which is at $`\lambda =1`$, $`\beta =0`$ and $`M=\mu =0`$. It is worth noting that there are alternative relevant deformations we could have considered involving the operators $`\mathrm{Tr}\widehat{\mathrm{\Phi }}^{+2}`$ and $`\mathrm{Tr}\widehat{\mathrm{\Phi }}^2`$. These are equivalent to (1) as a deformation of the $`𝒩=4`$ theory, i.e. when $`\beta =0`$, but differ once the Leigh-Strassler marginal deformation is present. In particular, the resulting theories differ in the IR. However, as they only differ by strictly relevant operators they both flow to the same fixed point in the UV: namely the Leigh-Strassler SCFT. The reason for choosing to study the specific relevant deformation (1) is simply that the resulting matrix model turns out to be exactly soluble .

Aside from masses for the chiral multiplets $`\widehat{\mathrm{\Phi }}^\pm `$ and $`\widehat{\mathrm{\Phi }}`$, our theory also involves two trilinear deformations, $`𝒪_1=\mathrm{Tr}\widehat{\mathrm{\Phi }}[\widehat{\mathrm{\Phi }}^+,\widehat{\mathrm{\Phi }}^{}]`$ and $`𝒪_2=\mathrm{Tr}\widehat{\mathrm{\Phi }}\{\widehat{\mathrm{\Phi }}^+,\widehat{\mathrm{\Phi }}^{}\}`$. Both the operators $`𝒪_1`$ and $`𝒪_2`$ are of course, marginal by power counting but only one of them, or more precisely, only one linear combination of these operators is an exactly marginal deformation of the $`𝒩=4`$ theory. Adding $`𝒪_1`$ to the $`𝒩=4`$ Lagrangian only changes the coefficient in front of the $`𝒩=4`$ superpotential, but it is an irrelevant operator even at one-loop at the $`𝒩=4`$ fixed point . The operator $`𝒪_1`$ is actually a descendant and not a chiral primary field in the $`𝒩=4`$ theory, hence its dimensions are not protected.333In the context of the AdS/CFT correspondence it is expected that these operators will get mass dimensions $`(g_{YM}^2N)^{1/4}`$  in the strongly coupled $`𝒩=4`$ theory in the large-$`N`$ limit. The second operator $`𝒪_2`$ on the other hand is known to be exactly marginal at one-loop at the $`𝒩=4`$ point . More generally, away from the $`𝒩=4`$ line one should expect, for fixed $`g_{YM}`$, a particular linear combination of $`𝒪_1`$ and $`𝒪_2`$ to be exactly marginal. We have parametrized this particular linear combination via two complex numbers $`\lambda `$ and $`\beta `$. In principle, $`\lambda `$ should be determined as a function of $`\beta `$ and $`g_{YM}`$ on the fixed manifold. Fortunately the specifics of this will not be important for us. Note that in the $`\text{SU}(2)`$ theory the operator $`𝒪_2`$ vanishes identically and the above superpotential does not yield a marginal deformation. We will assume that $`N>2`$ throughout this paper.

The theory with masses for all the fields, as in Eq. (1), has a number of vacua, many of which are massive. The canonical examples are the Higgs and confining vacua which will play an important rôle in our discussion. The confining vacua correspond to the trivial classical solution of the $`F`$-term conditions with $`\widehat{\mathrm{\Phi }}=\widehat{\mathrm{\Phi }}^\pm =0`$ and preserves the full $`\text{SU}(N)`$ gauge symmetry classically. At low energies the only classically massless fields comprise of the $`𝒩=1`$ gauge multiplet which confines and generates a mass gap. There is also a Higgs vacuum where at the classical level the gauge symmetry is completely broken. In addition, there are other massive vacua which are visible classically as solutions that leave a non-abelian gauge subgroup unbroken. Their classification is similar to that of the $`𝒩=1^{}`$ theory . What is interesting about the Leigh-Strassler deformation is that we shall find other vacua which are not present in the $`𝒩=1^{}`$ theory itself. In tandem with this, for special values of the deformation parameter $`\beta `$, such that $`e^{i\beta }`$ is a root of unity, vacua can disappear as a result of the condensates diverging.

In this paper, we compute the values of the condensates $`\mathrm{Tr}\widehat{\mathrm{\Phi }}^k`$ and the effective superpotential in all the massive vacua of the $`\text{SU}(N)`$ theory and we find that Montonen-Olive duality indeed relates the holomorphic condensates in different vacua with a special action on the couplings that we discuss below. This explicit computation is made possible by the recent proposal of Dijkgraaf and Vafa  wherein the effective superpotential of the gauge theory is computed by the genus zero free energy of a holomorphic three-matrix integral . The matrix model however, computes the superpotentials for the $`\text{U}(N)`$ gauge theory. One of our important conclusions is that the $`\text{SU}(N)`$ superpotential differs non-trivially from its $`\text{U}(N)`$ counterpart. (See also the work of  where a similar issue is discussed). We show however that the former can be unambiguously extracted from the latter by a process of “integrating in” of additional fields present in the $`\text{U}(N)`$ gauge theory. This difference turns out to be crucial as it is the $`\text{SU}(N)`$ results that clearly exhibit Montonen-Olive duality while the $`\text{U}(N)`$ results do not. For the $`\text{SU}(N)`$ theory we find in a generic $`(p,k)`$ massive vacuum (up to inconsequential vacuum-independent additive constants)

$$W_{\mathrm{eff}}^{\text{SU}(N)}=\frac{pN\mu M^2}{2\lambda ^2\mathrm{sin}\beta }\frac{\theta _1^{}(p\beta /2|\stackrel{~}{\tau })}{\theta _1(p\beta /2|\stackrel{~}{\tau })}.$$
(3)
with

$$\stackrel{~}{\tau }=\frac{p\tau _R+k}{q};\tau _R\tau \frac{iN}{\pi }\mathrm{ln}\lambda ;pq=N;k=0,1,\mathrm{},q1.$$
(4)
The main consequence of this result is that the values of the effective superpotential (and indeed all the eigenvalues of $`\widehat{\mathrm{\Phi }}`$) in different massive vacua of the theory with deformation parameter $`\beta `$ are mapped into one another by the action of the $`\text{SL}(2,)`$ transformations on the couplings:

$$\tau _R\frac{a\tau _R+b}{c\tau _R+d};\beta \frac{\beta }{c\tau _R+d};\lambda ^2\mathrm{sin}\beta \frac{\lambda ^2\mathrm{sin}\beta }{c\tau _R+d},$$
(5)
$`adbc=1`$; $`a,b,c,d`$. In particular, duality of the underlying SCFT is actually realized via modular transformations on a particular combination of the bare couplings $`\tau _R`$ rather than the gauge coupling $`\tau `$. In addition, $`\text{SL}(2,)`$ permutes vacua of a theory with deformation parameter $`\beta `$ provided $`\beta `$ transforms with modular weight $`1`$ as above along with a specific action on $`\lambda `$ deduced from above.

The vacuum structure and modular transformations described above can also be understood via an exact elliptic superpotential for the theory compactified on $`^3\times S^1`$. The superpotential is a function of $`N`$ chiral superfields $`X_a`$ ($`a=1,2,\mathrm{}N`$, $`_{a=1}^NX_a=0`$) which parameterize the Coulomb branch of the compactified theory. These are a complex combination of the Wilson lines and dual photons of the compactified theory. The superpotential

$$W_{^3\times S^1}=\frac{M^2\mu }{2\lambda ^2\mathrm{sin}\beta }\left(\frac{\omega _1}{\pi }\right)\underset{ab}{}\left(\zeta (X_aX_b+\beta \omega _1/\pi )\zeta (X_aX_b\beta \omega _1/\pi )\right)$$
(6)
reproduces the vacua described in Eqs. (1.3) and (1.4), taking the value (1.3) described above. Here, $`\zeta (z)`$ is the Weierstrass-$`\zeta `$ function for the torus with half-periods $`\omega _1=i\pi `$ and $`\omega _2=i\pi \tau _R`$, i.e. with complex structure $`\tau _R`$. This superpotential also predicts new $`\text{SL}(2,)`$-invariant vacua, not present in the $`𝒩=1^{}`$ theory, whose existence is confirmed by the classical analysis of the four-dimensional field theory.

We also remark that the properties of $`\beta `$ under modular transformations are similar to those of individual gauge couplings in elliptic quiver theories . In addition, at non-zero $`\beta `$ the matrix model solution naturally involves a torus with two marked points with $`\beta `$ parameterizing the separation. These suggest a deeper connection between the Leigh-Strassler and quiver SCFTs.

The layout of this paper is as follows: in Section 2 we show how one can relate the superpotentials for the $`\text{U}(N)`$ and $`\text{SU}(N)`$ theories by a process of integrating-in the trace part of the fields. This is very important because the matrix model approach yields the superpotential in the $`\text{U}(N)`$ theory, whereas, we are primarily interested in the $`\text{SU}(N)`$ theory. In Section 3, following Dijkgraaf and Vafa, we show how a matrix model can be used to find the superpotential in the confining vacua of the $`\text{U}(N)`$ theory and therefore by implication in $`\text{SU}(N)`$ theory. In particular, we find the spectrum of eigenvalues of the adjoint-valued field $`\widehat{\mathrm{\Phi }}`$ from which all the holomorphic condensates can be calculated. Section 4 is devoted to an analysis of the vacuum structure of the mass deformed Leigh-Strassler theory from the point-of-view of the tree-level superpotential. In particular, this leads to expressions for the classical limit of the condensates in each of the massive vacua. In addition, for the case of gauge group $`\text{SU}(3)`$ we find new vacua that are not present in the $`𝒩=1^{}`$ theory. We also show how some of the massive vacua can disappear when the coupling $`\beta `$ takes particular values. In Section 5, we return to the matrix model and consider multi-cut solutions that describe all the massive vacua. We then go on to how the structure of vacua can be used to deduce the action of the $`\text{SL}(2,)`$ duality group on the Leigh-Strassler theory itself. The final Section briefly reports on how the results from the matrix model can be used to deduce the exact elliptic superpotential (1.6) for the theory compactified on a circle to three dimensions, generalizing the one for the $`𝒩=1^{}`$ theory constructed in .

## 2 Relation between the $`\text{U}(N)`$ and $`\text{SU}(N)`$ theories

As we have explained one of our aims is to compute the quantum superpotential and a certain set of holomorphic observables in the massive vacua of the mass deformed $`\text{SU}(N)`$ Leigh-Strassler SCFT with $`𝒩=1`$ SUSY. For this purpose, we would like to employ matrix model techniques that have been proposed recently by Dijkgraaf and Vafa . However, already there is a subtlety. A direct application of the matrix model approach will solve the $`\text{U}(N)`$ theory since the DV proposal relates the effective superpotentials for $`\text{U}(N)`$ $`𝒩=1`$ gauge theories to the planar diagram expansion of corresponding matrix models. So the first question that we must address is how are the $`\text{SU}(N)`$ results related to those of the $`\text{U}(N)`$ gauge theory? We will now show that there is a very specific relation between the effective superpotential of the $`\text{SU}(N)`$ theory with classical superpotential Eq. (1), and the effective superpotential of the $`\text{U}(N)`$ theory with the same classical superpotential. This relation will eventually allow us to extract the $`\text{SU}(N)`$ results from the matrix model for the $`\text{U}(N)`$ gauge theory.

### 2.1 The $`\text{U}(N)`$ gauge theory

The Leigh-Strassler theory with $`\text{SU}(N)`$ gauge group differs non-trivially from its $`\text{U}(N)`$ counterpart. The $`\text{U}(N)`$ theory contains additional neutral chiral multiplets that couple to the chiral superfields transforming in the adjoint representation of the $`\text{SU}(N)\text{U}(N)`$. These interactions can modify the superpotential and other holomorphic observables of the $`\text{SU}(N)\text{U}(N)`$ theory. The way this happens is clearly seen from the point-of-view of the $`\text{U}(N)`$ theory. Let us begin by considering the $`\text{U}(N)`$ version of the theory, namely, the $`𝒩=1`$ SUSY gauge theory with tree level superpotential

$$W_{\mathrm{cl}}^{\text{U}(N)}=\mathrm{Tr}\left(i\lambda \mathrm{\Phi }[\mathrm{\Phi }^+,\mathrm{\Phi }^{}]_\beta +M\mathrm{\Phi }^+\mathrm{\Phi }^{}+\mu \mathrm{\Phi }^2\right),$$
(7)
where $`\mathrm{\Phi }^\pm `$ and $`\mathrm{\Phi }`$ are the fields in the adjoint of $`\text{U}(N)`$. Now the fields $`\mathrm{\Phi }^\pm `$ and $`\mathrm{\Phi }`$ can be naturally split into their traceless and trace parts:

$$\mathrm{\Phi }\widehat{\mathrm{\Phi }}+a;\mathrm{\Phi }^\pm \widehat{\mathrm{\Phi }}^\pm +a^\pm ;\mathrm{Tr}\widehat{\mathrm{\Phi }}0;\mathrm{Tr}\widehat{\mathrm{\Phi }}^\pm 0.$$
(8)
Here, $`\widehat{\mathrm{\Phi }}^\pm `$ and $`\widehat{\mathrm{\Phi }}`$ are traceless and so transform in the adjoint representation of $`\text{SU}(N)\text{U}(N)`$ while $`a^\pm =\mathrm{Tr}\mathrm{\Phi }^\pm /N`$ and $`a=\mathrm{Tr}\mathrm{\Phi }/N`$ are neutral. Rewritten in terms of these variables the tree level superpotential (7) for the $`\text{U}(N)`$ theory is

$$W_{\mathrm{cl}}^{\text{U}(N)}=N\left(M2\lambda a\mathrm{sin}\frac{\beta }{2}\right)a^+a^{}+N\mu a^2+W_{\mathrm{cl}}^{\text{SU}(N)}(a,a^\pm )$$
(9)
where

$$\begin{array}{cc}\hfill W_{\mathrm{cl}}^{\text{SU}(N)}(a,a^\pm )=& \mathrm{Tr}(i\lambda \widehat{\mathrm{\Phi }}[\widehat{\mathrm{\Phi }}^+,\widehat{\mathrm{\Phi }}^{}]_\beta +(M2\lambda a\mathrm{sin}\frac{\beta }{2})\widehat{\mathrm{\Phi }}^+\widehat{\mathrm{\Phi }}^{}+\mu \widehat{\mathrm{\Phi }}^2\hfill \\ & 2\lambda \mathrm{sin}\frac{\beta }{2}a^{}\widehat{\mathrm{\Phi }}\widehat{\mathrm{\Phi }}^+2\lambda \mathrm{sin}\frac{\beta }{2}a^+\widehat{\mathrm{\Phi }}\widehat{\mathrm{\Phi }}^{}).\hfill \end{array}$$
(10)
The main point here is that the neutral trace fields $`a`$ and $`a^\pm `$ have the effect of modifying the couplings of the $`\text{SU}(N)`$ fields. For example, the mass $`M`$ has been renormalized to $`M2\lambda a\mathrm{sin}\frac{\beta }{2}`$. In addition, there are new bilinears in $`\widehat{\mathrm{\Phi }}^\pm `$ and $`\widehat{\mathrm{\Phi }}`$ whose couplings444Strictly speaking, of course, these “couplings” are actually chiral superfields. But from the point of view the $`\text{SU}(N)`$ sub-sector of the $`\text{U}(N)`$ theory, these neutral chiral superfields do appear like couplings that have been elevated to chiral superfields otherwise known as “spurions”. depend on $`a,a^\pm `$. (Notice that when $`\beta =0`$, that is for the $`𝒩=1^{}`$ theory, these modifications disappear and there is no real difference between $`\text{SU}(N)`$ and $`\text{U}(N)`$.) In fact, we show below that because of the symmetries of the theory the VEVs of $`a^+`$ and $`a^{}`$ are forced to be zero self-consistently in the full theory. Thus the only effect of the neutral fields is to modify the mass term $`MM^{}M2a\lambda \mathrm{sin}\frac{\beta }{2}`$. One might then suspect that answers for pure $`\text{SU}(N)`$ theory of Eq. (1) may be obtained by simply rescaling the $`\text{U}(N)`$ results by appropriate powers of $`M^{}=M2a\lambda \mathrm{sin}\frac{\beta }{2}`$. We will now see that this is almost correct.

First of all, the fields $`a,a^\pm `$ are blind to all gauge interactions while $`\widehat{\mathrm{\Phi }},\widehat{\mathrm{\Phi }}^\pm `$ experience only $`\text{SU}(N)`$ gauge interactions, the remaining $`\text{U}(1)`$ gauge multiplet being decoupled from everything else. Hence in a vacuum where the $`\text{SU}(N)`$ fields are rendered massive and the gauge interactions generate an effective superpotential for that sector of the theory, we may readily write the effective superpotential for the $`\text{U}(N)`$ gauge theory as

$$W_{\mathrm{eff}}^{\text{U}(N)}=N\left(M2\lambda a\mathrm{sin}\frac{\beta }{2}\right)a^+a^{}+N\mu a^2+W_{\mathrm{eff}}^{\text{SU}(N)}(a,a^\pm ).$$
(11)
Note that the chiral fields $`a,a^\pm `$ appear as parameters or couplings for the $`\text{SU}(N)`$ sub-sector, however they are actually dynamical variables in the full theory and their values must be determined by extremizing the full effective superpotential. This procedure can actually be implemented formally by first noting that the $`\text{SU}(N)`$ sub-sector, from Eq. (10), has certain abelian symmetries. Let us define for the sake of convenience $`M^{}M2\lambda a\mathrm{sin}\frac{\beta }{2}`$. It is then sufficient to consider the following discrete symmetries: (i) $`(\widehat{\mathrm{\Phi }}^\pm ,a^\pm )`$ $``$ $`(\widehat{\mathrm{\Phi }}^\pm ,a^\pm )`$; (ii) $`(\widehat{\mathrm{\Phi }}^+,\widehat{\mathrm{\Phi }},a^+,M^{},\mu )`$ $``$ $`(i\widehat{\mathrm{\Phi }}^+,i\widehat{\mathrm{\Phi }},ia^+,iM^{},\mu )`$; and (iii) $`(\widehat{\mathrm{\Phi }}^{},\widehat{\mathrm{\Phi }},a^{},M^{},\mu )`$ $``$ $`(i\widehat{\mathrm{\Phi }}^{},i\widehat{\mathrm{\Phi }},ia^{},iM^{},\mu )`$. The only possible form for the $`\text{SU}(N)`$ effective superpotential, consistent with these symmetries, analytic in the parameters and which has mass dimension three, is

$$W_{\mathrm{eff}}^{\text{SU}(N)}=\mu M^2F(a^+a^{}/\mu M^{}),$$
(12)
for some unknown function $`F`$. Plugging this back into Eq. (11) for the effective $`\text{U}(N)`$ superpotential, and imposing the $`F`$-term conditions by extremizing with respect $`a,a^\pm `$, the only solution that generates a nontrivial effective superpotential is the one where

$$a^\pm =0;\mathrm{and}\frac{NM\mu }{2\lambda \mathrm{sin}\frac{\beta }{2}}a=W_{\mathrm{eff}}^{\text{U}(N)}.$$
(13)
This tells us that in a vacuum of the $`\text{U}(N)`$ gauge theory, the trace of the adjoint scalar $`\mathrm{\Phi }`$ must be related to the value of the effective superpotential of the $`\text{U}(N)`$ theory in that vacuum, precisely according to the above equation. This already constitutes a non-trivial prediction for the Dijkgraaf-Vafa approach for solving the $`\text{U}(N)`$ gauge theory, one that our results from the matrix model must satisfy.

### 2.2 The relation between the $`\text{U}(N)`$ and $`\text{SU}(N)`$ superpotentials

The symmetry arguments above also imply that the effective superpotential of the $`\text{SU}(N)`$ theory with a classical action as in Eq. (1) must be

$$W_{\mathrm{eff}}^{\text{SU}(N)}=\mu M^2F(0).$$
(14)
But we can easily determine the unknown function $`F`$ in terms of the $`\text{U}(N)`$ superpotential using Eqs. (11), (12) and (13) and we find:

$$W_{\mathrm{eff}}^{\text{SU}(N)}=\frac{W_{\mathrm{eff}}^{\text{U}(N)}}{1\frac{4\lambda ^2\mathrm{sin}^2\beta /2}{N\mu M^2}W_{\mathrm{eff}}^{\text{U}(N)}}.$$
(15)
Anticipating the matrix model results of the following section, we will write this in a slightly different way that may turn out to be more illuminating:

$$W_{\mathrm{eff}}^{\text{SU}(N)}=\frac{M^2}{M^2(a)}\left[W_{\mathrm{eff}}^{\text{U}(N)}+\frac{N\mu M^2}{4\lambda ^2\mathrm{sin}^2\frac{\beta }{2}}\right]\frac{N\mu M^2}{4\lambda ^2\mathrm{sin}^2\frac{\beta }{2}}.$$
(16)
In this form, it is apparent that the effective superpotentials of the $`\text{U}(N)`$ and $`\text{SU}(N)`$ gauge theories are indeed related by the replacement $`M^{}(a)M`$ after subtracting off certain additive constants. This form of the relation will turn out to be quite suggestive and useful when we discuss the matrix model results. Importantly, this simple relation tells us how to extract the $`\text{SU}(N)`$ answer from the $`\text{U}(N)`$ result which the matrix model naturally computes.

## 3 The $`\text{U}(N)`$ theory from the matrix model

According to the Dijkgraaf-Vafa proposal , in any given vacuum, the effective superpotential of the $`\text{U}(N)`$ $`𝒩=1`$ gauge theory with classical superpotential Eq. (1) is computed in terms of the planar diagram expansion of the three-matrix model partition function expanded around that vacuum

$$Z=[d\mathrm{\Phi }^+][d\mathrm{\Phi }^{}][d\mathrm{\Phi }]\mathrm{exp}g_s^1\mathrm{Tr}\left(i\lambda \mathrm{\Phi }[\mathrm{\Phi }^+,\mathrm{\Phi }^{}]_\beta +M\mathrm{\Phi }^+\mathrm{\Phi }^{}+\mu \mathrm{\Phi }^2\right).$$
(17)
We use the same notation for the matrix fields and the associated superfields in the $`\text{U}(N)`$ theory. In the matrix model, unlike the field theory, one takes $`\mathrm{\Phi }^+=(\mathrm{\Phi }^{})^{}`$ with the fluctuations of $`\mathrm{\Phi }`$ around the saddle point to be Hermitian. This matrix model has been actually solved in a different context  in the large-$`N`$ limit. First one integrates out $`\mathrm{\Phi }^\pm `$ and performs a field rescaling $`\mathrm{\Phi }\mathrm{\Phi }/\lambda `$ to get

$$Z=\lambda ^{N^2}[d\mathrm{\Phi }]\frac{\mathrm{exp}g_s^1\mu \mathrm{Tr}\mathrm{\Phi }^2/\lambda ^2}{|\mathrm{det}(M\mathrm{𝟏}\mathrm{𝟏}ie^{i\beta /2}\mathrm{\Phi }\mathrm{𝟏}+ie^{i\beta /2}\mathrm{𝟏}\mathrm{\Phi })|}.$$
(18)
Now we can follow  to obtain the saddle-point equation in the large-$`N`$ limit. For completeness we will now follow the steps required.

Let $`\{\varphi _i\}`$ denote the eigenvalues of the matrix $`\mathrm{\Phi }`$. Changing the integration variables in the matrix integral by going to the eigenvalue basis introduces a Jacobian: the famous Van der Monde determinant which leads to a repulsive force between the eigenvalues. The second step is a variable change that will eventually yield a simplified form for the saddle-point equation:

$$\varphi _i=Me^{\delta _i}+\frac{M}{2\mathrm{sin}\frac{\beta }{2}}.$$
(19)
In terms of these variables the classical potential for the matrix model eigenvalues takes the following form

$$\frac{\mu }{\lambda ^2}\mathrm{Tr}\mathrm{\Phi }^2=\underset{i}{}V(\delta _i)+\frac{\mu NM^2}{4\lambda ^2\mathrm{sin}^2\frac{\beta }{2}},$$
(20)
where we have defined

$$V(\delta )\frac{\mu M^2}{\lambda ^2}\left(e^{2\delta }\frac{e^\delta }{\mathrm{sin}\frac{\beta }{2}}\right).$$
(21)
In addition to this classical potential, the eigenvalues $`\delta _i`$ also experience pairwise effective interactions induced both by the Van der Monde determinant and the determinant resulting from integrating out $`\mathrm{\Phi }^\pm `$ in Eq. (18). The eigenvalues $`\delta _i`$ are naturally defined on the complex-$`z`$ plane with the identification $`zz+2\pi i`$, i.e. a cylinder.

### 3.1 Solving the large-$`N`$ matrix model for the confining vacuum

In the large-$`N`$ limit, the eigenvalues form a continuum and condense onto cuts in the complex plane. On can think of these cuts as arising from a quantum smearing-out of the classical eigenvalues of $`\mathrm{\Phi }`$. For the confining vacuum all the classical eigenvalues are degenerate, since $`\mathrm{\Phi }\mathrm{𝟏}`$, and so we expect a solution in the matrix model involving a single cut. Multi-cut solutions will be discussed later. The extent of the cut and the matrix model density of eigenvalues $`\rho (\delta )`$ can be determined from the saddle-point equation in terms of the parameters of the classical potential and the ’t Hooft coupling of the matrix model $`S=g_sN`$. The saddle-point equation is most conveniently formulated after defining the resolvent function:

$$\omega (z)=\frac{1}{2}_a^b𝑑\delta \frac{\rho (\delta )}{\mathrm{tanh}\frac{z\delta }{2}},\delta [a,b],_a^b\rho (\delta )𝑑\delta =1.$$
(22)
This function is analytic in $`z`$ and its only singularity is along a branch cut extending between $`[a,b]`$. The matrix model spectral density $`\rho (\delta )`$ is equal to the discontinuity across the cut

$$\omega (\delta +iϵ)\omega (\delta iϵ)=2\pi i\rho (\delta );\delta [a,b].$$
(23)
The saddle-point equation expresses the zero force condition on a test eigenvalue in the presence of the large-$`N`$ distribution of eigenvalues along the cut. For the matrix model that we are studying, this condition is best written in terms of the function $`G(z)`$ defined as

$$G(z)=U(z)+iS\left(\omega (z+i\frac{\beta }{2})\omega (zi\frac{\beta }{2})\right)$$
(24)
and $`\text{U}(z)`$ is a polynomial in $`e^z`$ such that

$$V^{}(z)=i\left(U(z+i\frac{\beta }{2})U(zi\frac{\beta }{2})\right).$$
(25)
From this definition, one may easily deduce that $`G(z)`$ is an analytic function on the cylinder and that in the interval $`|\mathrm{Im}z|\pi `$ it has two cuts $`[a+i\frac{\beta }{2},b+i\frac{\beta }{2}]`$ and $`[ai\frac{\beta }{2},bi\frac{\beta }{2}]`$. This is illustrated in Figure 1.

In terms of $`G(z)`$, the matrix model saddle-point equation is

$$G(\delta +i\frac{\beta }{2}\pm iϵ)=G(\delta i\frac{\beta }{2}iϵ);\delta [a,b].$$
(26)
The saddle-point equations actually provide a gluing condition on the cylinder with two branch cuts. In particular, this defines a torus with two marked points. $`G(z)`$ is then uniquely specified by gluing condition Eq. (26) and asymptotic behaviour at large $`z`$

$$\underset{z\mathrm{}}{lim}G(z)\frac{\mu M^2}{\lambda ^2}\left[\frac{1}{\mathrm{sin}\beta }e^{2z}\frac{1}{2\mathrm{sin}^2\frac{\beta }{2}}e^z+𝒪(e^z)\right],$$
(27)
which is a consequence of Eqs. (24) and (25). With this data we can determine $`G(z)`$ in a suitable parametrization so that we can then implement the DV proposal and obtain the $`\text{U}(N)`$ gauge theory superpotential.

### 3.2 The $`\stackrel{~}{\tau }`$-torus and the elliptic parametrization

As we stated earlier, the two-cut cylinder along with the asymptotic properties of $`G(z)`$ in Eq. (27) uniquely specifies a torus $`E_{\stackrel{~}{\tau }}`$ with complex structure parameter $`\stackrel{~}{\tau }`$. As shown in Figure 2, the contour $`𝒞_A`$ enclosing one of the cuts $`[a+i\frac{\beta }{2},b+i\frac{\beta }{2}]`$ anticlockwise maps to the $`A`$-cycle of the torus while the contour $`𝒞_B`$ joining the two cuts maps to the $`B`$-cycle. We will think of this torus as a quotient of the complex $`u`$-plane ($`u`$ being an auxiliary variable) by a lattice with periods $`2\stackrel{~}{\omega }_1`$ and $`2\stackrel{~}{\omega }_2`$ and with complex structure $`\stackrel{~}{\tau }=\stackrel{~}{\omega }_2/\stackrel{~}{\omega }_1`$.

¿From the field theory view-point, the appearance of this torus is not entirely obvious except for the fact that the gauge theory we are studying is connected to the $`𝒩=4`$ theory. In the $`𝒩=1^{}`$ case discussed in , one might have expected the torus to appear due to the relation with the Seiberg-Witten curve of the associated $`𝒩=2^{}`$ theory. For the Leigh-Strassler SCFT the fact that the torus structure is also present is one of the most important conclusions of this paper: as we shall see it leads to an understanding of how $`\text{SL}(2,)`$ acts on the underlying SCFT.

The map $`z(u)`$ from the $`u`$-plane to the $`z`$-plane is specified by the requirements that going around $`𝒞_A`$ (the $`A`$-cycle of the torus) returns $`z`$ to its original value, while traversing the contour $`𝒞_B`$ (the $`B`$-cycle of the torus) causes $`z`$ to jump by an amount $`i\beta `$, which is the distance between the two cuts in the $`z`$-plane. Both these operations leave $`G`$ unchanged which is therefore an elliptic function on the $`u`$-plane. Thus

$`A\mathrm{cycle}:z(u+2\stackrel{~}{\omega }_1)=z(u);G(z(u+2\stackrel{~}{\omega }_1))=G(z(u)),`$ (28a)
$`B\mathrm{cycle}:z(u+2\stackrel{~}{\omega }_2)=z(u)+i\beta ;G(z(u+2\stackrel{~}{\omega }_2))=G(z(u)).`$ (28b)
This determines the following unique map $`z(u)`$ from the $`u`$-plane to the two-cut cylinder:

$$\mathrm{exp}z(u)=\frac{1}{2\mathrm{tan}\frac{\beta }{2}}\frac{\theta _1^{}(0|\stackrel{~}{\tau })\theta _1(\pi u/2\stackrel{~}{\omega }_1\beta /4|\stackrel{~}{\tau })}{\theta _1^{}(\beta /2|\stackrel{~}{\tau })\theta _1(\pi u/2\stackrel{~}{\omega }_1+\beta /4|\stackrel{~}{\tau })}$$
(29)
The fact that $`G(z(u))`$ is an elliptic function on the $`u`$-plane along with its asymptotic behaviour at large $`z`$ (27) (corresponding to a pole in the $`u`$-plane) leads uniquely to

$$G(z(u))=\frac{\mu M^2\stackrel{~}{\omega }_1^2\mathrm{cos}\frac{\beta }{2}}{4\pi ^2\lambda ^2\mathrm{sin}^3\frac{\beta }{2}}\frac{\theta _1(\beta /2|\stackrel{~}{\tau })^2}{\theta _1^{}(\beta /2|\stackrel{~}{\tau })^2}\left[\mathrm{}(u+\stackrel{~}{\omega }_1\beta /2\pi )\mathrm{}(\stackrel{~}{\omega }_1\beta /\pi )\right].$$
(30)
Here, $`\mathrm{}(u)`$ is the Weierstrass function, an elliptic function defined on $`E_{\stackrel{~}{\tau }}`$, and $`\theta _1(u|\tau )`$ is the Jacobi theta function of the first kind (see  and Appendix A for more details). As expected, $`G(z(u))`$ is an elliptic function of $`u`$ with a second order pole in the $`u`$-plane; $`z(u)`$ on the other hand is only quasi-elliptic. For a different choice of the deformation potential the asymptotic behaviour of $`G(z)`$ will change and $`G(z(u))`$ will end up being an elliptic function of a different order, but significantly $`z(u)`$ will remain unchanged. There are two special points on the torus $`E_{\stackrel{~}{\tau }}`$ at $`u=\pm \beta \stackrel{~}{\omega }_1/2\pi `$ which map to the points $`z=\mathrm{}`$. Having determined $`G(z)`$ in the elliptic parametrization we can now implement the DV proposal to compute the superpotential in the confining vacuum of the $`\text{U}(N)`$ gauge theory.

### 3.3 Superpotential of the $`\text{U}(N)`$ theory via Dijkgraaf-Vafa

According the Dijkgraaf-Vafa proposal, the gluino condensate of the gauge theory gets identified with the ’t Hooft coupling $`S`$ of the matrix model. From Eqs. (22) and (24), the integral of $`G(z)`$ around one of its cuts, say $`𝒞_A`$, is equal to $`2\pi S`$, so we can directly compute the gluino condensate of the $`\text{U}(N)`$ gauge theory in the elliptic parametrization:

$$\mathrm{\Pi }_A=2\pi iS=i_AG(z(u))\frac{dz(u)}{du}𝑑u.$$
(31)
The second ingredient required to determine the QFT superpotential is the variation in the genus zero free energy of the matrix model in transporting a test eigenvalue from infinity to the endpoint of the cut. This is obtained by integrating the force on a test eigenvalue from infinity to the endpoint of the cut and can be expressed as $`\mathrm{\Pi }_B=i_BG(z)𝑑z`$. There is a slight subtlety here as there is a “zero-point” energy associated with additive constants in the matrix model potential. First note that there is a multiplicative factor of $`\lambda ^{N^2}`$ in the partition function (18) and an additive constant in the classical scalar potential (20). Both these appear as $`S`$-dependent contributions to the matrix model genus zero free energy defined as $`_0/g_s^2=N^2_0/S^2`$, so that

$$\frac{_0}{S}=\mathrm{\Pi }_B2S\mathrm{ln}\lambda +\frac{\mu M^2}{4\lambda ^2\mathrm{sin}^2\frac{\beta }{2}};\mathrm{\Pi }_B=i_BG(z(u))\frac{dz(u)}{du}𝑑u.$$
(32)
The effective superpotential of the mass deformation of the $`\text{U}(N)`$ Leigh-Strassler theory in its confining vacuum is obtained by extremizing the following expression with respect to $`S`$:

$$W_{\mathrm{eff}}^{\text{U}(N)}=N\frac{_0}{S}2\pi i\tau S;\frac{}{S}W_{\mathrm{eff}}^{\text{U}(N)}=0.$$
(33)
where $`\tau `$ is the bare coupling of the theory. The $`S`$-independent constant in Eq. (32) appears as a vacuum-independent additive constant in the superpotential. On the other hand, the term linear in $`S`$ in Eq. (32) actually has the effect of renormalizing the bare coupling $`\tau `$ so that

$$\tau \tau _R=\tau \frac{iN}{\pi }\mathrm{ln}\lambda .$$
(34)
Both $`\mathrm{\Pi }_A`$ and $`\mathrm{\Pi }_B`$ are easy to evaluate in our elliptic parametrization:

$$\mathrm{\Pi }_A=\frac{dh(\stackrel{~}{\tau })}{d\stackrel{~}{\tau }};\mathrm{\Pi }_B=\stackrel{~}{\tau }\frac{dh(\stackrel{~}{\tau })}{d\stackrel{~}{\tau }}h(\stackrel{~}{\tau }),$$
(35)
where

$$h(\stackrel{~}{\tau })=\frac{\mu M^2\mathrm{cos}\frac{\beta }{2}}{4\lambda ^2\mathrm{sin}^3\frac{\beta }{2}}\frac{\theta _1(\beta /2|\stackrel{~}{\tau })}{\theta _1^{}(\beta /2|\stackrel{~}{\tau })}.$$
(36)
It follows straightforwardly that

$$\frac{}{S}W_{eff}^{\text{U}(N)}=0\stackrel{~}{\tau }=\frac{\tau _R}{N}$$
(37)
so that

$$W_{\mathrm{eff}}^{\text{U}(N)}=Nh(\tau _R/N)+\frac{N\mu M^2}{4\lambda ^2\mathrm{sin}^2\frac{\beta }{2}}.$$
(38)
Finally we have, therefore, the quantum superpotential of the $`\text{U}(N)`$ gauge theory in its confining vacuum,

$$W_{\mathrm{eff}}^{\text{U}(N)}=\frac{N\mu M^2\mathrm{cos}\frac{\beta }{2}}{4\lambda ^2\mathrm{sin}^3\frac{\beta }{2}}\frac{\theta _1(\beta /2|\tau _R/N)}{\theta _1^{}(\beta /2|\tau _R/N)}+\frac{N\mu M^2}{4\lambda ^2\mathrm{sin}^2\frac{\beta }{2}},$$
(39)
with $`\tau _R`$ defined in (34). This result satisfies the first obvious check: in the $`\beta 0`$ limit with $`\lambda =1`$ it reproduces the superpotential of the $`𝒩=1^{}`$ theory. The superpotential in the other $`N1`$ confining vacua are obtained by repeated $`T`$-dualities $`\tau _R\tau _R+1`$. Note that the $`T`$-dualities in question naturally act on $`\tau _R`$ rather than the bare coupling $`\tau `$. Similarly the result for the superpotential has a modular parameter $`\tau _R/N`$ and so we expect that any modular properties we uncover should naturally be a consequence of the modular group acting on $`\tau _R`$ the “renormalized” coupling.

### 3.4 Eigenvalues of $`\mathrm{\Phi }`$ in the $`\text{U}(N)`$ theory and a test

One of the main results of our earlier work  in the context of $`𝒩=1^{}`$ theory was that the map $`z(u)`$ from the $`u`$-plane to the $`z`$-plane, when evaluated along the $`A`$-cycle of the $`\stackrel{~}{\tau }`$-torus (which maps to the upper cut on the $`z`$-plane), precisely yields the eigenvalues of the adjoint scalar $`\mathrm{\Phi }`$ of the gauge theory in its confining vacuum, distributed uniformly along this cycle.555In contrast the matrix model eigenvalues are not distributed uniformly along this cycle signalling the crucial difference between the spectral densities of the two systems. Nevertheless, the matrix model provides a means to extract the field theory distributions in a simple way. We were able to prove this assertion for a class of deformations of the $`𝒩=4`$ theory involving masses for two of the adjoint chiral superfields and arbitrary polynomials in the third adjoint field $`\mathrm{\Phi }`$. We will not pursue this generalized analysis in the present case, but we will present strong evidence that the same relation with the field theory eigenvalues holds under the Leigh-Strassler deformation.

For the Leigh-Strassler theory, we first recall that the $`\delta _i`$ are valued in the $`z`$-plane and these are in turn related to the eigenvalues of the $`\mathrm{\Phi }`$ matrix through Eq. (19). Taking into account the fact that the function $`G(z)`$ is defined in Eq. (24) in terms of the resolvent functions $`\omega (z\pm i\frac{\beta }{2})`$ where the positions of the cuts of the latter have been displaced by $`\pm i\frac{\beta }{2}`$, we conjecture that the function $`M\mathrm{exp}(z(u)i\frac{\beta }{2})/\lambda `$ evaluated along the $`A`$-cycle gives the eigenvalues of $`\mathrm{\Phi }`$ in the $`\text{U}(N)`$ theory:

$$\varphi (x)=\frac{M}{2\lambda \mathrm{sin}\frac{\beta }{2}}\left[\mathrm{cos}\frac{\beta }{2}\frac{\theta _1^{}(0|\tau _R/N)}{\theta _1^{}(\beta /2|\tau _R/N)}\frac{\theta _3(\pi x\beta /4|\tau _R/N)}{\theta _3(\pi x+\beta /4|\tau _R/N)}+1\right];x[\frac{1}{2},\frac{1}{2}].$$
(40)
Note that the factor $`\lambda `$ takes into account the field rescaling prior to Eq. (19) and we also incorporate the additive shift of $`M/(2\lambda \mathrm{sin}\frac{\beta }{2})`$ involved in the variable change in Eq. (19). This is of course a continuous distribution of eigenvalues appropriate for the $`\text{U}(N)`$ theory in the large-$`N`$ limit. But the finite $`N`$ answer is obtained by simply replacing the continuous index $`x`$ by a discrete index. For this to be legitimate, it is essential that the field theory eigenvalues be uniformly distributed in the interval $`x[\frac{1}{2},\frac{1}{2}]`$ as we have indeed pointed out earlier. The difference between finite-$`N`$ and large-$`N`$ manifests itself in gauge-invariant observables as vacuum-independent operator mixing ambiguities that vanish in the large-$`N`$ limit .

Our results for the eigenvalues in the confining vacuum of the $`\text{U}(N)`$ theory satisfy several checks. Firstly in the classical limit $`\tau _Ri\mathrm{}`$, using the properties of the Jacobi theta functions (see Appendix A), Eq. (40) vanishes as expected, since the confining vacuum should correspond to the trivial classical solution. Secondly, in the $`\beta 0`$ limit with $`\lambda =1`$ the eigenvalues of $`𝒩=1^{}`$ theory in its confining vacuum reproduce the expressions in . Two further nontrivial tests involve the computation of $`\mathrm{Tr}\mathrm{\Phi }`$ and $`\mathrm{Tr}\mathrm{\Phi }^2`$ in the $`\text{U}(N)`$ theory using the eigenvalues above. With a large-$`N`$ distribution, the trace reduces to an integral and we find

$`\mathrm{Tr}\mathrm{\Phi }=N{\displaystyle _{1/2}^{1/2}}𝑑x\varphi (x)={\displaystyle \frac{2\lambda \mathrm{sin}\frac{\beta }{2}}{\mu M}}W_{\mathrm{eff}}^{\text{U}(N)},`$ (41a)
$`\mathrm{Tr}\mathrm{\Phi }^2=N{\displaystyle _{1/2}^{1/2}}𝑑x\varphi ^2(x)={\displaystyle \frac{1}{\mu }}W_{\mathrm{eff}}^{\text{U}(N)},`$ (41b)
where $`W_{\mathrm{eff}}^{\text{U}(N)}`$ is the expression obtained in Eq. (39). We show how to compute the associated integrals explicitly in Appendix B. The first (41a) demonstrates that the trace of $`\mathrm{\Phi }`$ calculated using the matrix model results satisfies precisely the relation Eq. (13) predicted on general grounds from field theory. The second integral is also exactly in accord with field theory expectations since from our general arguments in Section 2, we have $`W_{\mathrm{eff}}^{\text{U}(N)}\mu `$ and therefore differentiating with respect to $`\mu `$ gives $`\mathrm{Tr}\mathrm{\Phi }^2=W_{\mathrm{eff}}^{\text{U}(N)}/\mu `$ precisely (41b).

In summary, the eigenvalues of the adjoint scalar $`\mathrm{\Phi }`$ in the confining vacuum of the mass deformed $`\text{U}(N)`$ Leigh-Strassler theory at large $`N`$ are given by Eq. (40). The eigenvalues of the $`\text{U}(N)`$ theory at finite $`N`$ are obtained by substituting the continuous index $`x`$ with an appropriate discrete index.

### 3.5 The $`\text{SU}(N)`$ superpotential

We argued in Section 2 that the superpotentials of the $`\text{U}(N)`$ and $`\text{SU}(N)`$ theories differ non-trivially, but are nevertheless simply related by Eq. (16). The main point there was that the $`\text{U}(N)`$ theory allows a non-zero $`\mathrm{Tr}\mathrm{\Phi }Na0`$ to develop in the confining vacuum of the quantum theory. This leads to what is effectively a renormalization of the mass parameter $`MM^{}=M2a\lambda \mathrm{sin}\frac{\beta }{2}`$. According to the result of our field theory analysis of Section 2, the $`\text{SU}(N)`$ superpotential may be recovered by a rescaling of the $`\text{U}(N)`$ result by $`M^2/M^2`$, after subtracting off a vacuum independent constant. It is quite interesting to note that precisely the same vacuum independent constant makes an appearance in the $`\text{U}(N)`$ superpotential (39) computed in the Dijkgraaf-Vafa approach.

Using our results for the $`\text{U}(N)`$ superpotential from the previous section and Eq. (16), we thus find that in the confining vacuum the $`\text{SU}(N)`$ theory develops the superpotential,

$$W_{\mathrm{eff}}^{\text{SU}(N)}=\frac{N\mu M^2}{2\lambda ^2\mathrm{sin}\beta }\frac{\theta _1^{}(\beta /2|\tau _R/N)}{\theta _1(\beta /2|\tau _R/N)}\frac{N\mu M^2}{4\lambda ^2\mathrm{sin}\frac{\beta }{2}},$$
(42)
where the renormalized coupling is defined in (34). This expression for the $`\text{SU}(N)`$ superpotential has the correct $`\beta 0`$ limit in that it reproduces the $`𝒩=1^{}`$ superpotential exactly.

### 3.6 $`\text{SU}(N)`$ eigenvalues and condensates

The eigenvalues $`\widehat{\varphi }_i`$ of the adjoint scalar $`\widehat{\mathrm{\Phi }}`$ in the confining vacuum of the $`\text{SU}(N)`$ theory are related to the $`\text{U}(N)`$ eigenvalues by a simple extension of our arguments in Section 2 as

$$\widehat{\varphi }(x)=\frac{M}{M^{}}\left(\varphi (x)\frac{M}{2\lambda \mathrm{sin}\frac{\beta }{2}}\right)+\frac{M}{2\lambda \mathrm{sin}\frac{\beta }{2}}.$$
(43)
This can be understood simply as a rescaling of the $`\text{U}(N)`$ spectrum of eigenvalues by $`M/M^{}(a)`$ after subtracting out a constant and is in line with the arguments of Section 2 for extracting the $`\text{SU}(N)`$ superpotential.

Evaluating this explicitly yields

$$\widehat{\varphi }(x)=\frac{M}{2\lambda }\left[\frac{\theta _1^{}(0|\tau _R/N)}{\theta _1(\beta /2|\tau _R/N)}\frac{\theta _3(\pi x\beta /4|\tau _R/N)}{\theta _3(\pi x+\beta /4|\tau _R/N)}+\frac{1}{\mathrm{sin}\frac{\beta }{2}}\right];x[\frac{1}{2},\frac{1}{2}].$$
(44)
We show in Appendix B that with these eigenvalues $`\mathrm{Tr}\widehat{\mathrm{\Phi }}=0`$ and $`\mathrm{Tr}\widehat{\mathrm{\Phi }}^2=W_{\mathrm{eff}}^{\text{SU}(N)}/\mu `$. Once again we remark that with a continuous label $`x`$ these are really the eigenvalues of the large-$`N`$ gauge theory but finite-$`N`$ results are obtained by simply replacing $`x`$ with a discrete index.

Now that we have the complete distribution of eigenvalues of $`\widehat{\mathrm{\Phi }}`$ in the confining vacuum of the $`\text{SU}(N)`$ theory, every condensate of the form $`\mathrm{Tr}\widehat{\mathrm{\Phi }}^k`$ is automatically computed via the integrals

$$\mathrm{Tr}\widehat{\mathrm{\Phi }}^k=N_{1/2}^{1/2}𝑑x\widehat{\varphi }^k(x).$$
(45)
We will not undertake the exercise of computing these but, if desired, they can be computed following the methods in Appendix B.

### 3.7 Summary

In this section we have computed the value of the effective superpotential and the spectrum of eigenvalues of $`\widehat{\mathrm{\Phi }}`$ in the confining vacuum of the $`\text{SU}(N)`$ theory. We did this by obtaining the $`\text{U}(N)`$ superpotential via a matrix model and subsequently extracting the $`\text{SU}(N)`$ results by a process of integrating-in666We remark that one might have expected to obtain the $`\text{SU}(N)`$ superpotential and eigenvalues directly from the matrix model by imposing a constraint on the matrices via a Lagrange multiplier. However this turns out not to work in any simple way, the primary reason being, as pointed out in , that the matrix model eigenvalue spectrum is completely different from the field theory spectrum and tracelessness of the latter translates into a rather complicated constraint on the former. discussed in Section 2. The results of the computation have satisfied several general field theory checks. In addition, the emergence of a quasi-modular function ($`\theta _1^{}/\theta _1`$) as the superpotential of the $`\text{SU}(N)`$ theory is strongly suggestive of Montonen-Olive duality being realized in this theory. (In contrast, the $`\text{U}(N)`$ superpotential being basically the reciprocal of the $`\text{SU}(N)`$ result, transforms in a complicated way under the modular group.) In order to establish electric-magnetic duality we need to be able to compute the corresponding quantities in other massive vacua of the theory (most notably the Higgs vacuum which should be $`S`$-dual to the confining phase). We will address this in Sections 4 and 5 which will be devoted to understanding the classical and quantum vacuum structure respectively.

## 4 Vacuum structure of deformed $`\text{SU}(N)`$ Leigh-Strassler theory

In this section we exhibit the classical vacuum structure of the $`\text{SU}(N)`$ theory. Importantly these classical configurations will serve as the starting point for analysing other saddle points of the matrix model and thence obtaining the effective quantum superpotential in all the massive vacua. Interestingly the classical solutions (the eigenvalues) for the Higgs vacuum that we find here will be shown to be related by $`S`$-duality to the quantum eigenvalue spectrum of the confining vacuum that we have already obtained in the previous section.

The vacuum structure of the mass deformed $`\text{SU}(N)`$ Leigh-Strassler theory turns out to be subtly different from that of the $`𝒩=1^{}`$ theory: there are new vacua and old vacua can disappear. In order to describe the vacuum structure, we start with the tree-level superpotential (1). At the classical level vacua are determined by solving the $`F`$\- and $`D`$-flatness conditions modulo gauge transformations. As usual, if one is interested in holomorphic information it is sufficient to solve the $`F`$-flatness conditions modulo complex gauge transformations. These conditions follow from extremizing the tree-level superpotential (1):777Note that the $`F`$-term equations are rather similar to the “$`q`$-deformed” $`\text{SU}(2)`$ algebra. They are similar to, but not the same as the equations encountered in . The authors of  considered a slightly different mass deformation and a $`\text{U}(N)`$ gauge group.

$$[\widehat{\mathrm{\Phi }}^+,\widehat{\mathrm{\Phi }}^{}]_\beta =i\frac{2\mu }{\lambda }\widehat{\mathrm{\Phi }},[\widehat{\mathrm{\Phi }},\widehat{\mathrm{\Phi }}^+]=i\frac{M}{\lambda }\widehat{\mathrm{\Phi }}^+,[\widehat{\mathrm{\Phi }}^{},\widehat{\mathrm{\Phi }}]=i\frac{M}{\lambda }\widehat{\mathrm{\Phi }}^{}.$$
(46)
In the above, $``$ is the projector onto the traceless part which must be imposed since we are dealing with the $`\text{SU}(N)`$, rather than the $`\text{U}(N)`$, theory. Complex gauge transformation can be used to diagonalize $`\widehat{\mathrm{\Phi }}`$.

In the limit $`\beta =0`$, the theory becomes the basic $`𝒩=1^{}`$ theory and the conditions (46) reduce to the Lie algebra of $`\text{SU}(2)`$. Let us briefly recall the solutions in that case. Each solution is associated to a representation of $`\text{SU}(2)`$ of dimension $`N`$ where $`\widehat{\mathrm{\Phi }}`$ and $`\widehat{\mathrm{\Phi }}^\pm `$ are, up to a trivial re-scaling, $`J_3`$ and $`J_\pm `$, respectively. Each representation of $`\text{SU}(2)`$ of dimension $`N`$ is associated to a partition $`Nn_1+\mathrm{}+n_q`$, where each integer $`n_i`$ is the dimension of each irreducible component. In general the unbroken gauge group has abelian factors and the phase is Coulomb, or massless. The massive vacua are associated to the the equi-partitions $`N=qp`$ for which the unbroken gauge group is $`\text{SU}(q)`$. In such a vacuum, the low energy description is pure $`𝒩=1`$ with gauge group $`\text{SU}(q)`$ which has $`q`$ inequivalent quantum vacua. Hence, there are $`_{q|N}q`$ vacua with a mass gap (the sum being over the positive integer divisors of $`N`$).

Now we consider the case when $`\beta 0`$. Each of the vacua of the $`𝒩=1^{}`$ theory remains a vacuum of the deformed theory, at least for generic values of $`\beta `$. These deformed solutions have the following form. Firstly, the block form of the solutions associated to the partition $`Nn_1+\mathrm{}+n_q`$ is preserved. Let us concentrate on one of the blocks of dimension $`n`$. For $`\beta =0`$, $`\widehat{\mathrm{\Phi }}`$ and $`\widehat{\mathrm{\Phi }}^\pm `$ are equal to $`J_3`$ and $`J_\pm `$, up to trivial re-scaling, of the irreducible $`\text{SU}(2)`$ representation of dimension $`n`$:

$$\beta =0:\widehat{\mathrm{\Phi }}=\frac{M}{\lambda }J_3;\mathrm{\Phi }^\pm =\frac{2\mu M}{\lambda }J_\pm .$$
(47)
When $`\beta 0`$, the non-zero elements of $`\widehat{\mathrm{\Phi }}^\pm `$ have the same pattern as $`J_\pm `$ (so only the elements just above and below the diagonal, respectively are non-zero) but the numerical values are changed. The elements of diagonal matrix $`\widehat{\mathrm{\Phi }}`$ in the block in question are found to be

$$\widehat{\varphi }_j=\frac{M}{2\lambda }\left(\frac{1}{\mathrm{sin}\frac{\beta }{2}}\frac{n}{\mathrm{sin}\frac{n\beta }{2}}e^{i\beta (n+12j)/2}\right),j=1,\mathrm{},n.$$
(48)
Notice that when $`\beta 0`$, $`\widehat{\mathrm{\Phi }}`$, in this block, becomes equal to $`(M/\lambda )J_3`$ as expected.

For generic $`\beta `$, these deformed vacua are massive/massless according to their type when $`\beta =0`$. However, for specific values of $`\beta `$ vacua can migrate off to infinity in field space and hence cease to be genuine vacua of the field theory. This occurs when the eigenvalues of $`\widehat{\mathrm{\Phi }}`$ diverge. It is easy to see that this will occur in a block of dimension $`n`$ when the deformation parameter

$$\beta =\frac{2\pi k}{n},k.$$
(49)
In one of the $`N/p`$ massive vacua with unbroken $`\text{SU}(N/p)`$ symmetry the eigenvalues of $`\widehat{\mathrm{\Phi }}`$ are given by (48) with $`n=p`$ and where each eigenvalue has degeneracy $`N/p`$. From this spectrum, we can calculate the classical values for the condensates; for instance, for the quadratic condensate:

$$\mathrm{Tr}\widehat{\mathrm{\Phi }}^2=\frac{M^2N}{4\lambda ^2\mathrm{sin}^2\frac{\beta }{2}}\left(p\frac{\mathrm{tan}\frac{\beta }{2}}{\mathrm{tan}\frac{p\beta }{2}}1\right).$$
(50)
Notice that this condensate diverges when $`\beta =2\pi k/p`$, $`k`$ in accord with (49).

### 4.1 New classical branches

As well as the original massive vacua, there are also new vacua which are not present when $`\beta =0`$. We illustrate this is in the simplest case where the gauge group is $`\text{SU}(3)`$. In this case there is a new vacuum, more precisely a moduli space of vacua, in which $`\widehat{\mathrm{\Phi }}`$ and $`\widehat{\mathrm{\Phi }}^\pm `$ are all, up to complex gauge transformations, diagonal:

$$\begin{array}{cc}\hfill \widehat{\mathrm{\Phi }}& =\frac{M}{2\lambda \mathrm{sin}\frac{\beta }{2}}\mathrm{diag}(x_1,x_2,x_3),\hfill \\ \hfill \widehat{\mathrm{\Phi }}^+& =\frac{\sqrt{\mu M/2}}{\lambda \mathrm{sin}\frac{\beta }{2}}\rho \mathrm{diag}((1x_1)^1,(1x_2)^1,(1x_3)^1),\hfill \\ \hfill \widehat{\mathrm{\Phi }}^{}& =\frac{\sqrt{\mu M/2}}{\lambda \mathrm{sin}\frac{\beta }{2}}\xi \mathrm{diag}((1x_1)^1,(1x_2)^1,(1x_3)^1).\hfill \end{array}$$
(51)
Here, $`\rho `$ and $`\xi `$ are two complex moduli of the solution and $`x_i`$, $`i=1,2,3`$ are the three distinct roots of the cubic

$$x^33x+2\rho \xi =0.$$
(52)
Note that the tracelessness of $`\widehat{\mathrm{\Phi }}`$ and $`\widehat{\mathrm{\Phi }}^\pm `$ is guaranteed by the form of the cubic (52). Clearly these solutions diverge to infinity when $`\beta 0`$. However, the moduli space of the solutions also has additional branches which occur whenever

$$e^{i\beta /2}x_ie^{i\beta /2}x_j=2i\mathrm{sin}\frac{\beta }{2}.$$
(53)
In this case, an off-diagonal component of either of $`(\widehat{\mathrm{\Phi }}^\pm )_{ij}`$ can become non-zero. The quadratic condensate in this space of vacua is

$$\text{SU}(3):\mathrm{Tr}\widehat{\mathrm{\Phi }}^2=\frac{3M^2}{2\lambda ^2\mathrm{sin}^2\frac{\beta }{2}}.$$
(54)
Clearly this diverges when $`\beta =0`$. We will see that these new vacua are important for checking a conjecture that we make later.

## 5 Electric-magnetic duality in the Leigh-Strassler Deformation

The question we now address is what becomes of the $`S`$-duality of the $`𝒩=4`$ theory under the Leigh-Strassler deformation? Specifically we are considering a deformation of the $`𝒩=4`$ Lagrangian of the form;

$$=_{𝒩=4}+d^2\theta \mathrm{\Delta }W$$
(55)
where the deforming superpotential can be written as

$$\mathrm{\Delta }W=i\left(\lambda \mathrm{cos}\frac{\beta }{2}1\right)𝒪_1\lambda \mathrm{sin}\frac{\beta }{2}𝒪_2+M\mathrm{Tr}\widehat{\mathrm{\Phi }}^+\widehat{\mathrm{\Phi }}^{}+\mu \mathrm{Tr}\widehat{\mathrm{\Phi }}^2$$
(56)
where, as in Section 1, $`𝒪_1=\mathrm{Tr}\widehat{\mathrm{\Phi }}[\widehat{\mathrm{\Phi }}^+,\widehat{\mathrm{\Phi }}^{}]`$ and $`𝒪_2=\mathrm{Tr}\widehat{\mathrm{\Phi }}\{\widehat{\mathrm{\Phi }}^+,\widehat{\mathrm{\Phi }}^{}\}`$. At the $`𝒩=4`$ point, the operators appearing in the deforming superpotential transform under $`S`$-duality in a definite way. Adding these operators to the superpotential therefore breaks modular invariance. However, as is standard when dealing with broken symmetries, we should be able to restore invariance by allowing the corresponding couplings to transform in an appropriate way. This strategy was implemented in detail for the $`𝒩=1^{}`$ case $`\beta =0`$ in . We will now review this briefly and extend this discussion to the marginal deformation.

Using the conventions described in ,888See the final paragraph of Section 3 of this reference. the superspace measure $`d^2\theta `$ has holomorphic modular weight $`2`$. Hence for modular invariance of the Lagrangian, the superpotential must have weight $`+2`$. In the $`𝒩=4`$ theory, the mass operators $`\mathrm{Tr}\widehat{\mathrm{\Phi }}^+\widehat{\mathrm{\Phi }}^{}`$ and $`\mathrm{Tr}\widehat{\mathrm{\Phi }}^2`$ are both part of the multiplet of chiral primary operators which transform as second rank symmetric traceless tensors under the $`\text{SO}(6)`$ R-symmetry. Their modular weights can be deduced via the AdS/CFT correspondence by comparison with the dual fields in IIB supergravity . In fact they have holomorphic modular weight $`+2`$ and the Lagrangian will be invariant if we assign the mass parameters $`M`$ and $`\mu `$ weight zero. Although there is no explicit breaking of modular invariance, it is important to remember that the theory has various vacuum states in which these mass operators acquire expectation values, thereby ‘spontaneously’ breaking $`S`$-duality. In other words, modular transformations map the observables in one vacuum to those in another. Indeed, this property is manifest in the exact results of .

Similar arguments should apply to the Leigh-Strassler theory at least near the $`𝒩=4`$ point where the modular properties of the deforming operators are known. If we set $`\lambda =1`$ and work to linear order in $`\beta `$, the deformation only contains the cubic operator $`𝒪_2`$. This operator is also part of an $`𝒩=4`$ chiral primary multiplet consisting, in this case, of third rank symmetric traceless tensors of $`SO(6)`$. In the conventions of  it has modular weight $`+3`$. To restore invariance under $`S`$-duality we should therefore assign the parameter $`\beta `$ holomorphic modular weight $`1`$. In Section 5.2, we will find the appropriate generalization of this transformation property throughout the parameter space, by using the matrix model to find the exact superpotential in each of the massive vacua. In Section 5.1 we see what can be deduced by looking at the Higgs and confining vacua alone.

### 5.1 A first look at $`S`$-duality

Let us concentrate on the Higgs and confining vacua. In the $`𝒩=1^{}`$ limit, it is well known that each observable in the Higgs vacuum is simply related by $`S`$-duality, $`\tau 1/\tau `$, to the corresponding observable in the confining vacuum. This is an explicit realization of the Montonen-Olive electric-magnetic duality of the $`𝒩=4`$ theory. The question is whether a similar duality persists under the Leigh-Strassler deformation.

If $`S`$-duality is realized in this theory it should exchange the Higgs and confining vacua. Since we already have the eigenvalues in the confining vacuum of the $`\text{SU}(N)`$ theory a simple and powerful test emerges: use $`S`$-duality to deduce the eigenvalues in the Higgs vacuum and then take the classical limit and compare with the spectrum (48) ($`n=N`$) deduced from the $`F`$-flatness conditions. In the next section, we shall use the matrix model to solve for the superpotential in each of the massive vacua and this uncovers the action of the whole of the modular group.

The eigenvalues of $`\widehat{\mathrm{\Phi }}`$ in the confining vacuum are given in (44). (Recall that the eigenvalues in this vacuum have the property that they vanish in the classical limit.) The first thing to notice is that the modular parameter of the elliptic function parametrizing the eigenvalue spectrum is $`\tau _R/N`$ with $`\tau _R=\tau (iN/\pi )\mathrm{ln}\lambda `$, where $`\tau `$ is the gauge coupling. Hence $`S`$-duality, if realized, must act on the renormalized coupling $`\tau _R`$ as $`\tau _R1/\tau _R`$ rather than $`\tau `$ itself. Under this action, the modular parameter of the elliptic functions changes as $`\tau _R/N1/(\tau _RN)`$. But now we want to take the classical limit $`\tau _Ri\mathrm{}`$ and this can be done easily after performing a modular transformation $`1/(\tau _RN)\tau _RN`$ on the (quasi-)elliptic functions. Using the modular transformation properties of the theta functions (see Appendix A) we find

$$\widehat{\varphi }(x)\tau _R\frac{NM}{2\lambda \mathrm{sin}\frac{N\tau _R\beta }{2}}e^{iN\tau _R\beta x}+\frac{M}{2\lambda \mathrm{sin}\frac{\beta }{2}};x[\frac{1}{2},\frac{1}{2}].$$
(57)
Apart form the additive constant, this coincides precisely with the expression Eq. (48) ($`n=N`$) for the eigenvalues in the Higgs vacuum of the deformed $`\text{SU}(N)`$ Leigh-Strassler theory with $`\beta `$ replaced by $`\tau _R\beta `$. Equivalently, the Higgs and confining vacua of the theory with deformation parameter $`\beta `$ are simply exchanged under the combined operation $`\tau _R1/\tau _R`$ and $`\beta \beta /\tau _R`$

The overall factor of $`\tau _R`$ in the above expression is expected because the eigenvalues have modular weight one. The fact that the additive constants are different is to be expected: even in the context of $`𝒩=1^{}`$ theory the matrix model approach only yields results that are $`S`$-duality covariant up to a vacuum independent shift . This additive shift will be discussed in more detail in the next section. We should also point out that after $`S`$-duality on the spectrum of eigenvalues in the confining vacuum we have taken a classical limit above. In this limit, the exact transformation properties of the couplings $`\lambda ,\beta `$ are not visible. These will become clear in the next section.

The results of this simple analysis are that there is indeed a form of $`S`$-duality that relates the Higgs and confining vacua of the theory with deformation parameter $`\beta `$. In particular this duality acts on the bare couplings of the Leigh-Strassler theory as

$$\tau _R1/\tau _R,\beta \beta /\tau _R.$$
(58)
This is consistent with $`\beta `$ having modular weight -1 as expected. Recall that the bare theory has 3 couplings $`\tau `$, $`\lambda `$ and $`\beta `$ with one non-trivial relation, so that (58) determines the action of $`S`$-duality on the bare theory up to the unknown relation. In the next section we expand the discussion of duality to include the whole of the modular group $`\text{SL}(2,)`$ and deduce the action on all three of the bare couplings.

### 5.2 Multi-cut solutions of the matrix model and $`\text{SL}(2,)`$

Now that we have described the structure of massive vacua of the Leigh-Strassler deformed theory, we can return to the matrix model and ask whether we can solve for the exact superpotential and condensates in each of these vacua. This will turn out to be very important because it will prove that there is an action of electric-magnetic duality on the Leigh-Strassler deformation of the $`𝒩=4`$ theory. In fact this extends to the action of the full $`\text{SL}(2,)`$ modular group.

The way to use the matrix model to solve for the theory in an arbitrary massive vacuum is described for the $`𝒩=1^{}`$ theory in . One finds that the same techniques can be used in the deformed matrix model described in Section 3. In order to describe the massive vacua with unbroken $`\text{SU}(N/p)`$ symmetry, one needs a multi-cut solution of the matrix model. Of course the matrix model computes the superpotential in the $`\text{U}(N)`$ theory and we then use (15) to deduce it in the $`\text{SU}(N)`$ theory.

We start with the observation that the classical eigenvalues in the massive vacua with unbroken $`\text{SU}(N/p)`$ are given by (48) with $`n=p`$, each with degeneracy $`N/p`$. The only thing that will be important for us will be the range of the phases of these classical eigenvalues. The classical eigenvalues of the $`\text{U}(N)`$ gauge theory also have exactly the same phases in the corresponding massive vacua. As we have seen explicitly in the case of the confining vacuum, the eigenvalues of the field theory are essentially proportional to $`\mathrm{exp}[z]`$. This directs us to consider a $`p`$-cut solution of the matrix model with the cuts in the $`z`$-plane situated at

$$[a+i\beta (p+12j)/2,b+i\beta (p+12j)/2];j=1,\mathrm{},p$$
(59)
where these positions are dictated precisely by the phases of the classical eigenvalues. The classical eigenvalues are clustered in groups of $`N/p`$ at each of these points, but the pairwise effective interactions induced in the matrix model will cause each cluster to spread out into a cut.

Following , we then make an ansatz where the cuts have the same length and the same filling fractions. It can then be shown that $`G(z)`$ defined in (24) again only has two cuts, as in the confining vacuum, but where the cuts are shifted to $`\pm [a+i\beta p/2,b+i\beta p/2]`$. This modifies the map $`z(u)`$ in (29) by a shift $`\beta p\beta `$:

$$\mathrm{exp}z(u)=\frac{1}{2\mathrm{tan}\frac{p\beta }{2}}\frac{\theta _1^{}(0|\stackrel{~}{\tau })\theta _1(\pi u/2\stackrel{~}{\omega }_1p\beta /4|\stackrel{~}{\tau })}{\theta _1^{}(p\beta /2|\stackrel{~}{\tau })\theta _1(\pi u/2\stackrel{~}{\omega }_1+p\beta /4|\stackrel{~}{\tau })}.$$
(60)
Similarly

$$G(z(u))=\frac{\mu M^2\stackrel{~}{\omega }_1^2\mathrm{cos}\frac{p\beta }{2}}{4\pi ^2\lambda ^2\mathrm{sin}^3\frac{p\beta }{2}}\frac{\theta _1(p\beta /2|\stackrel{~}{\tau })^2}{\theta _1^{}(p\beta /2|\stackrel{~}{\tau })^2}\left[\mathrm{}(u+\stackrel{~}{\omega }_1p\beta /2\pi )\mathrm{}(\stackrel{~}{\omega }_1p\beta /\pi )\right].$$
(61)
When the Dijkgraaf-Vafa recipe is followed one finds that for each cut there is an associated filling fraction of eigenvalues $`N_j/N`$ which translates to an effective ’t Hooft coupling $`S_j=g_sN_j`$ with $`S_j=S`$. For the massive vacua we note that the filling fractions in each of the $`p`$ cuts are the same and equal to $`N/p=q`$

$$2\pi iS_j=\frac{dh(\stackrel{~}{\tau })}{d\stackrel{~}{\tau }};j=1,\mathrm{},p,$$
(62)
where $`h(\stackrel{~}{\tau })`$ is the same function that we encountered in (36) with $`\beta p\beta `$:

$$h(\stackrel{~}{\tau })=\frac{\mu M^2\mathrm{cos}\frac{p\beta }{2}}{4\lambda ^2\mathrm{sin}^3\frac{p\beta }{2}}\frac{\theta _1(p\beta /2|\stackrel{~}{\tau })}{\theta _1^{}(p\beta /2|\stackrel{~}{\tau })}.$$
(63)
In addition one also needs the variation in the genus zero free energy upon transporting an eigenvalue from any given cut to infinity. After minor modifications this is also given by the same expressions that were encountered in the one-cut solution for the confining vacuum:

$$\underset{j=1}{\overset{p}{}}\frac{_0}{S_j}=\stackrel{~}{\tau }\frac{dh(\stackrel{~}{\tau })}{d\stackrel{~}{\tau }}h(\stackrel{~}{\tau })2\mathrm{ln}\lambda \underset{j=1}{\overset{p}{}}S_j+p\frac{\mu M^2}{4\lambda ^2\mathrm{sin}^2\frac{\beta }{2}}.$$
(64)
Then the superpotential for the $`\text{U}(N)`$ theory at any of its massive vacua is

$$W_{\mathrm{eff}}^{\text{U}(N)}=\frac{N}{p}\underset{j=1}{\overset{p}{}}\frac{_0}{S_j}2\pi i\tau \underset{j=1}{\overset{p}{}}S_j.$$
(65)
Extremizing with respect to $`\stackrel{~}{\tau }`$ gives

$$\stackrel{~}{\tau }=\frac{p^2}{N}\tau _R$$
(66)
and therefore we find that the value of the quantum effective superpotential for the $`\text{U}(N)`$ theory is

$$W_{\mathrm{eff}}^{\text{U}(N)}=\frac{N\mu M^2\mathrm{cos}\frac{p\beta }{2}}{4p\lambda ^2\mathrm{sin}^3\frac{p\beta }{2}}\frac{\theta _1(p\beta /2|p^2\tau _R/N)}{\theta _1^{}(p\beta /2|p^2\tau _R/N)}+\frac{N\mu M^2}{4\lambda ^2\mathrm{sin}^2\frac{\beta }{2}}.$$
(67)
In all of the above we should remember that $`p`$ is a divisor of $`N`$.

The result for $`\text{SU}(N)`$ follows by the application of Eq. (15):

$$W_{\mathrm{eff}}^{\text{SU}(N)}=\frac{pN\mu M^2}{2\lambda ^2\mathrm{sin}\beta }\frac{\theta _1^{}(p\beta /2|p^2\tau _R/N)}{\theta _1(p\beta /2|p^2\tau _R/N)}\frac{N\mu M^2}{4\lambda ^2\mathrm{sin}^2\frac{\beta }{2}}.$$
(68)
This gives the result in only one of $`q=N/p`$ vacua of this type. The results in the remaining vacua follow by the replacement

$$\tau _R\tau _R+k/p;k=0,1,\mathrm{}q1.$$
(69)
The Higgs vacuum corresponds to $`p=N`$ whilst the confining vacuum considered previously has $`p=1`$ and $`k=0`$.

One interesting and immediate check on the form of the superpotential follows from the fact that it is manifestly divergent when the deformation parameter takes the special values

$$\beta =\frac{2\pi k}{p},k.$$
(70)
But this matches precisely the values of $`\beta `$ in Eq. (49) for which the classical eigenvalues of $`\widehat{\mathrm{\Phi }}`$ for the massive vacua diverge (since they consist of $`q`$ blocks of size $`n=p`$). Another check follows by taking the classical limit of the superpotential (68). Using (90c), one finds

$$W_{\mathrm{eff}}^{\text{SU}(N)}\stackrel{\tau _Ri\mathrm{}}{}\mathrm{Tr}\widehat{\mathrm{\Phi }}^2=\frac{N\mu M^2}{4\lambda ^2\mathrm{sin}^2\frac{\beta }{2}}\left(p\frac{\mathrm{tan}\frac{\beta }{2}}{\mathrm{tan}\frac{p\beta }{2}}1\right),$$
(71)
which matches the values of $`\mu \mathrm{Tr}\widehat{\mathrm{\Phi }}^2`$ obtained in (50) from the $`F`$-flatness condition precisely.

Now that we have the superpotential in each of the massive vacua, we can now investigate how the full $`\text{SL}(2,)`$ duality group acts on the theory. In the $`𝒩=1^{}`$ theory described in , we established that the matrix model yields a result for the superpotential which is only duality covariant when shifted by a vacuum independent constant. This kind of additive ambiguity is ubiquitous in this subject: see . Naturally, the same kind of ambiguity will be present in the Leigh-Strassler deformation where we claim that the appropriate definition of a duality covariant superpotential is

$$\stackrel{~}{W}_{\mathrm{eff}}^{\text{SU}(N)}=W_{\mathrm{eff}}^{\text{SU}(N)}+C(\tau _R,\beta ,\lambda ),$$
(72)
where the vacuum-independent constant is given by either of the expressions

$$\begin{array}{cc}\hfill C(\tau _R,\beta ,\lambda )& =\frac{N\mu M^2}{4\lambda ^2\mathrm{sin}^2\frac{\beta }{2}}\frac{N\mu M^2}{2\lambda ^2\mathrm{sin}\beta }\frac{\theta _1^{}(\beta /2|\tau _R)}{\theta _1(\beta /2|\tau _R)}+\frac{N(N1)\mu M^2\beta }{12\lambda ^2\mathrm{sin}\beta }E_2(\tau _R)\hfill \\ & =\frac{N\mu M^2}{4\lambda ^2\mathrm{sin}^2\frac{\beta }{2}}\frac{N\mu M^2\omega _1}{\pi \lambda ^2\mathrm{sin}\beta }\left(\zeta (\omega _1\beta /\pi )\frac{N\beta }{\pi }\zeta (\omega _1)\right).\hfill \end{array}$$
(73)
Here, the Weierstrass $`\zeta `$-function is defined relative to a torus with periods $`\omega _1`$ and $`\omega _2`$ and $`\tau _R=\omega _2/\omega _1`$ and $`E_2(\tau _R)`$ is the $`2^{\mathrm{nd}}`$ Eisenstein series (see the Appendix of  for definitions and references). One property that the $`C`$ has is that it reduces to the correct form needed in the $`𝒩=1^{}`$ theory in the limit $`\beta 0`$.

Our claim is that the duality-covariant definition of the superpotential in the each of the massive vacua labelled by $`(p,k)`$, $`pq=N`$ and $`k=0,\mathrm{},q`$, is then

$$\stackrel{~}{W}_{\mathrm{eff}}^{\text{SU}(N)}=\frac{N\mu M^2}{2\lambda ^2\mathrm{sin}\beta }\left(p\frac{\theta _1^{}(p\beta /2|\stackrel{~}{\tau })}{\theta _1(p\beta /2|\stackrel{~}{\tau })}\frac{\theta _1^{}(\beta /2|\tau _R)}{\theta _1(\beta /2|\tau _R)}+\frac{\beta }{6}(N1)E_2(\tau _R)\right),$$
(74)
where

$$\stackrel{~}{\tau }=(p\tau _R+k)/q.$$
(75)
Using the modular properties of the theta functions and the $`2^{\mathrm{nd}}`$ Eisenstein series, (89a) and (92), one can easily show, for example, that

$$\stackrel{~}{W}_{\mathrm{eff}}(\beta /\tau _R,1/\tau _R,\lambda )|_{\mathrm{Conf}}=\tau _R\frac{\mathrm{sin}\beta }{\mathrm{sin}(\beta /\tau _R)}\stackrel{~}{W}_{\mathrm{eff}}(\beta ,\tau _R,\lambda )|_{\mathrm{Higgs}}.$$
(76)
Notice that the coupling $`\beta `$ transforms with unit weight under the modular group. Note also the rather unconventional transformation of the superpotential itself: this we will have to cure by an appropriate definition of the transformation of the coupling $`\lambda `$.

A consideration of the other massive vacua show that they all lie on a single orbit of the $`\text{SL}(2,)`$ duality group. For the element

$$\sigma =\left(\begin{array}{cc}a& b\\ c& d\end{array}\right)\text{SL}(2,),$$
(77)
which relates two vacua $`A=\sigma (B)`$, we have999The actual action of $`\sigma `$ on a vacuum with labels $`(p,k)`$ is not difficult to elucidate.

$$\stackrel{~}{W}_{\mathrm{eff}}(\frac{\beta }{c\tau _R+d},\frac{a\tau _R+b}{c\tau _R+d},\lambda )|_A=(c\tau _R+d)\frac{\mathrm{sin}\beta }{\mathrm{sin}[\beta /(c\tau _R+d)]}\stackrel{~}{W}_{\mathrm{eff}}(\beta ,\tau _R,\lambda )|_B.$$
(78)
From this we can deduce the action of $`\sigma \text{SL}(2,)`$ on the bare couplings of the theory; clearly

$$\sigma (\tau _R)=\frac{a\tau _R+b}{c\tau _R+d},\sigma (\beta )=\frac{\beta }{c\tau _R+d},$$
(79)
confirming that $`\beta `$ transforms with modular weight -1. However, we can say more: in order that the superpotential transforms with a definite modular weight, so that it is consistent with the group action, we can deduce the transformation of the coupling $`\lambda `$. Recall, that on the marginal surface of the Leigh-Strassler deformation $`\lambda `$ is a function—albeit unknown—of the other two couplings and so must transform. In addition, since in the limit $`\beta 0`$, the superpotential transforms with modular weight 2, we expect this to be preserved for $`\beta 0`$. This requires the combination $`\lambda ^2\mathrm{sin}\beta `$ that appears as an overall factor in the superpotential must transform with unit modular weight:

$$\sigma (\lambda ^2\mathrm{sin}\beta )=\frac{\lambda ^2\mathrm{sin}\beta }{c\tau _R+d}.$$
(80)
## 6 An exact elliptic superpotential

Another method that has been used to calculate the exact values of the superpotential in all the vacua of the $`𝒩=1^{}`$ theory is to compactify it on a circle of finite radius . The effective superpotential is then a function of the dual photons and Wilson lines of the abelian subgroup $`\text{U}(1)^{N1}\text{SU}(N)`$. These comprise $`N1`$ complex scalar fields $`X_a`$, $`a=\mathrm{},N`$ (with $`_{a=1}^NX_a=0`$) which naturally live on a torus of complex structure $`\tau `$ because of the periodicity of each dual photon and Wilson line. The superpotential describing the $`𝒩=1^{}`$ deformation is therefore constrained to be an elliptic function of the complex scalars $`X_a`$. In  the exact answer was found to be

$$W_{^3\times S^1}\underset{ab}{}\mathrm{}(X_aX_b),$$
(81)
where $`\mathrm{}(z)`$ is the Weierstrass function defined on a torus with half periods $`\omega _{1,2}`$ with $`\tau =\omega _2/\omega _1`$. What is particularly useful about this superpotential is that it is independent of the compactification radius, and therefore, yields results that are valid in the four-dimensional limit, and it encodes all the vacua of the $`𝒩=1^{}`$ theory: both massive and massless.

The question is whether there is an elliptic superpotential that describes the Leigh-Strassler deformation? Similar arguments to those used in  suggest that there should be, but are not powerful enough to determine the exact answer. However, there is an obvious candidate for appropriate marginal deformation of (81), which turns out to be correct. Specifically, we can replace the Weierstrass function $`\mathrm{}(z)`$, by its unique deformation in the space of even elliptic functions.101010The function must be even to respect the Weyl group of $`\text{SU}(N)`$ In particular, the Weierstrass function has one double pole at the origin and we can consider the elliptic function of order two obtained by splitting the double pole into two single poles. Remarkably, a superpotential of precisely this form reproduces our results from the matrix model for the superpotential in each of the massive vacua. This superpotential can be written as,

$$W_{^3\times S^1}=\frac{M^2\mu }{2\lambda ^2\mathrm{sin}\beta }\left(\frac{\omega _1}{\pi }\right)\underset{ab}{}\left(\zeta (X_aX_b+\beta \omega _1/\pi )\zeta (X_aX_b\beta \omega _1/\pi )\right).$$
(82)
The Weierstrass-zeta function $`\zeta (u)`$ can be defined via $`\mathrm{}(u)=\zeta ^{}(u)`$ where $`\mathrm{}(u)`$ is the Weierstrass function for the torus with complex structure $`\tau _R`$. Then in the limit $`\beta =0`$ this reproduces (81). What is remarkable is that the massive vacua of the $`𝒩=1^{}`$ theory for which the $`X_a`$ lie on sub-lattices are also critical points of (82) since this only relies on the fact that the two-body interaction is an even elliptic function. For instance, for the $`pq=N`$ massive vacuum with $`k=0`$ we have

$$X_a\{\frac{2\omega _1r}{p}+\frac{2\omega _2s}{q},r=0,\mathrm{},p1,s=0,\mathrm{},q1\}.$$
(83)
Furthermore, the values of the superpotential are precisely given by (74).

A very strong check on the form of the superpotential we have conjectured in (82) is possible because of the existence of new vacua which are not present in the limit $`\beta =0`$. For example, it is easy to show that $`X_a=0`$, $`a`$, is a critical point. Actually this is an $`\text{SL}(2,)`$-invariant vacuum. Now we can take the classical limit, $`\tau i\mathrm{}`$, after having shifted by minus the vacuum-independent shift (73), and compare with the analysis of the vacuum structure in Section 4. In the case of $`\text{SU}(3)`$, the classical limit of the superpotential $`W_{^3\times S^1}C`$ in the vacuum with $`X_a=0`$ precisely equals the value calculated from the new vacuum solution of the $`F`$-flatness conditions (54). We expect this to generalize to $`\text{SU}(N)`$, however, the analysis remains to be done.

The properties and construction of this elliptic superpotential will be described in a separate publication; however, we mention the fact that for $`𝒩=1^{}`$, the elliptic superpotential is, more or less, the (complexified) Hamiltonian of the elliptic Calogero-Moser integrable system which described $`N`$ particles interacting though pair-wise forces proportion to $`\mathrm{}(X_aX_b)`$. The intriguing question is whether this system can be deformed

$$\mathrm{}(X_aX_b)\zeta (X_aX_b+\beta \omega _1/\pi )\zeta (X_aX_b\beta \omega _1/\pi )$$
(84)
whilst preserving integrability.

Finally we comment on a remarkable correspondence between our results for the Leigh-Strassler deformation and the exact superpotential of a certain relevant deformation of the $`A_1`$-quiver theory studied in . Specifically this is the theory of an $`\text{SU}(N)\times \text{SU}(N)`$ gauge multiplet with gauge couplings $`g_1`$ and $`g_2`$, with two chiral multiplets $`A`$ and $`B`$ in the $`(N,\overline{N})`$ representation and two chiral multiplets $`\stackrel{~}{A}`$ and $`\stackrel{~}{B}`$ in the $`(\overline{N},N)`$. The theory also contains two chiral multiplets $`\widehat{\mathrm{\Phi }}_1`$ and $`\widehat{\mathrm{\Phi }}_2`$ in the $`(\mathrm{adj},0)`$ and $`(0,\mathrm{adj})`$ representations respectively with classical superpotential

$$W=\mathrm{Tr}\left\{(g_1\widehat{\mathrm{\Phi }}_1+M/2)(A\stackrel{~}{A}+B\stackrel{~}{B})+(g_2\widehat{\mathrm{\Phi }}_2+M/2)(\stackrel{~}{A}A+\stackrel{~}{B}B)\right\}.$$
(85)
This theory actually has $`𝒩=2`$ supersymmetry but we further deform it by adding the operator,

$$\mu \mathrm{Tr}\left(\widehat{\mathrm{\Phi }}_1^2\widehat{\mathrm{\Phi }}_2^2\right).$$
(86)
to the classical superpotential which preserves $`𝒩=1`$ supersymmetry. This theory is a relevant deformation of the $`𝒩=2`$ superconformal theory obtained by setting $`\mu =M=0`$. The exact vacuum structure and superpotential was determined in . The theory has $`N_{d|N}d`$ massive vacua. Remarkably, in $`_{d|N}d`$ of these vacua, we find exactly the same value of the superpotential as in each of the massive vacua of the $`\text{SU}(N)`$ Leigh-Strassler theory studied in this paper. Specifically, this agreement holds if we identify the parameter $`\beta `$ of the Leigh-Strassler deformation and the complexified coupling of one of the $`\text{SU}(N)`$ factors in the quiver theory $`\tau _1=4\pi i/g_1^2+\theta _1/2\pi `$ (or $`z=2i\pi \tau _1`$ to use the notation of ) :

$$z=\frac{N\omega _1\beta }{\pi }.$$
(87)
In the brane set-up for the quiver model $`z`$ is the relative separation of the two NS5-branes. We also identify the parameters $`M`$ and $`\mu `$ of the quiver superpotential with parameters of the same name in the deformed Leigh-Strassler superpotential (1).

This connection suggests that there should be an RG flow from the superconformal $`A_1`$-quiver to the Leigh-Strassler fixed line. At least in one special case, such a flow is known to exist. In , Klebanov and Witten considered the quiver theory in the case of equal gauge couplings $`g_1=g_2`$ and zero masses for the bi-fundamental multiplets, $`M=0`$. In this case one can integrate out the adjoint chiral fields $`\widehat{\mathrm{\Phi }}_1`$ and $`\widehat{\mathrm{\Phi }}_2`$ and induce a marginal quartic superpotential for the bi-fundamentals which corresponds to D3-branes placed at a resolved conifold singularity. Moving onto the Higgs branch where $`\text{SU}(N)\times \text{SU}(N)`$ is broken to the diagonal $`\text{SU}(N)`$, we find an effective theory of three chiral superfields transforming in the adjoint of the unbroken $`\text{SU}(N)`$ which coincides with $`𝒩=4`$ SUSY Yang-Mills. Repeating this exercise with unequal gauge couplings $`g_1g_2`$, we find that the Leigh-Strassler marginal deformation as well as other irrelevant operators are induced. However, it is not obvious whether this flow is related to the exact agreement we have found above.

Acknowledgements: TJH would like to thank Frank Ferrari for discussions. The authors would also like to thank Robbert Dijkgraaf and Cumrun Vafa for valuable discussions. We thank Robbert Dijkgraaf for his comments and observations on a preliminary version of this article. SPK acknowledges support from a PPARC Advanced Fellowship.

Appendix A: Some properties of elliptic functions

The (quasi-)elliptic functions that we need, $`\mathrm{}(z)`$, $`\zeta (z)`$ and $`\theta _i(z|\tau )`$, are all standard and their definitions and properties can be found in . In our conventions, the basic torus has periods $`2\omega _1`$ and $`2\omega _2`$ with complex structure $`\tau =\omega _2/\omega _1`$. We also define $`q=\mathrm{exp}(i\pi \tau )`$. Below we list some properties of these functions that will be useful.

An important relation between the Weierstrass and Jacobi theta functions is

$$\zeta (u)\frac{\zeta (\omega _1)}{\omega _1}u=\frac{\pi }{2\omega _1}\frac{\theta _1^{}(\pi u/(2\omega _1)|\tau )}{\theta _1(\pi u/(2\omega _1)|\tau )}.$$
(88)
The Jacobi theta functions have the following modular properties that will be useful

$`\theta _1(x|\tau )`$ $`=i(i\tau )^{{\scriptscriptstyle \frac{1}{2}}}\mathrm{exp}(ix^2/\pi \tau )\theta _1(x/\tau |1/\tau ),`$ (89a)
$`\theta _3(x|\tau )`$ $`=(i\tau )^{{\scriptscriptstyle \frac{1}{2}}}\mathrm{exp}(ix^2/\pi \tau )\theta _3(x/\tau |1/\tau ).`$ (89b)
In order to take the classical limit, we will need the expansion

$`\theta _1(x|\tau )`$ $`=2{\displaystyle \underset{n=0}{\overset{\mathrm{}}{}}}(1)^nq^{(n+1/2)^2}\mathrm{sin}(2n+1)x,`$ (90a)
$`\theta _3(x|\tau )`$ $`=1+2{\displaystyle \underset{n=0}{\overset{\mathrm{}}{}}}(1)^nq^{n^2}\mathrm{cos}2nx,`$ (90b)
$`{\displaystyle \frac{\theta _1^{}(x|\tau )}{\theta _1(x|\tau )}}`$ $`=\mathrm{cot}x+4{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{q^{2n}}{1q^{2n}}}\mathrm{sin}2nx.`$ (90c)
We also need the $`2^{\mathrm{nd}}`$ Eisenstein series

$$E_2(\tau )=124\underset{n=1}{}\sigma _1(n)q^{2n},$$
(91)
where $`\sigma _1(n)`$ is a sum over each positive integral divisor of $`n`$. It has the following non-trivial modular transformation

$$E_2(1/\tau )=\tau ^2E_2(\tau )+\frac{6\tau }{\pi i}.$$
(92)
Appendix B: Calculating condensates

In this appendix we show how the integrals of the type (41a) and (41b) can be computed to yield the condensates in the $`\text{U}(N)`$ and $`\text{SU}(N)`$ gauge theories. Let us first consider the function $`f(x)`$ in terms of which the eigenvalues have been defined

$$f(x)=\frac{\theta _4(\pi x\beta /4|\stackrel{~}{\tau })}{\theta _4(\pi x+\beta /4|\stackrel{~}{\tau })}\frac{\theta _3(\pi x\pi /2\beta /4|\stackrel{~}{\tau })}{\theta _3(\pi x\pi /2+\beta /4|\stackrel{~}{\tau })};x[0,1].$$
(93)
Clearly, the eigenvalues can be written either in terms of $`\theta _3`$ or $`\theta _4`$ since the only difference is in the range of the argument $`x`$. We also define the function $`g(x)`$

$$g(x)=\frac{\theta _1(\pi x\beta /4|\stackrel{~}{\tau })}{\theta _1(\pi x+\beta /4|\stackrel{~}{\tau })}.$$
(94)
The functions $`f`$ and $`g`$ have some nice properties

$$g(x+1)=g(x);g(x+\stackrel{~}{\tau }/2)=e^{i\beta /2}f(x);g(x\stackrel{~}{\tau }/2)=e^{i\beta /2}f(x),$$
(95)
which follow from standard identities. Note also that the function $`g(x)`$ has a simple pole in its fundamental period at $`x=1\beta /4\pi `$. Now it is a simple matter to convince oneself that because of the periodicity properties of $`g(x)`$ the contour integral around this pole and hence the residue there can also be written as

$$\mathrm{Res}g(x)|_{x=\beta /4\pi }=2i\mathrm{sin}\frac{\beta }{2}_0^1f(x)𝑑x.$$
(96)
Thus we see that

$$_0^1f(x)𝑑x=\frac{1}{\mathrm{sin}\frac{\beta }{2}}\frac{\theta _1(\beta /2|\stackrel{~}{\tau })}{\theta _1^{}(0|\stackrel{~}{\tau })}.$$
(97)
From this key result it follows using Eq. (40)

$$\mathrm{Tr}\mathrm{\Phi }=N_0^1𝑑x\varphi (x)=\frac{NM}{2\lambda \mathrm{sin}\frac{\beta }{2}}\left[\frac{1}{\mathrm{tan}\frac{\beta }{2}}\frac{\theta _1(\beta /2|\stackrel{~}{\tau })}{\theta _1^{}(\beta /2|\stackrel{~}{\tau })}+1\right]=\frac{2\lambda \mathrm{sin}\frac{\beta }{2}}{M\mu }W_{\mathrm{eff}}^{\text{U}(N)}.$$
(98)
(Note that because $`\varphi (x)`$ is a periodic function under $`xx+1`$ the actual limits of integration can always be changed from $`(1/2,1/2)`$ to $`(0,1)`$.) Similarly it also follows that in the $`\text{SU}(N)`$ theory

$$\mathrm{Tr}\widehat{\mathrm{\Phi }}=N_0^1\widehat{\varphi }(x)𝑑x=0.$$
(99)
To compute the quadratic condensates we need the integrals $`_0^1f^2(x)𝑑x`$. These can also be calculated using the same trick as above and we find that

$$_0^1f^2(x)𝑑x=\frac{1}{2i\mathrm{sin}\beta }\mathrm{Res}g^2(x)|_{x=\beta /4\pi }=\frac{2\theta _1^{}(\beta /2|\stackrel{~}{\tau })\theta _1(\beta /2|\stackrel{~}{\tau })}{\mathrm{sin}\beta \theta _1^{}(0|\stackrel{~}{\tau })^2}.$$
(100)
Using the results Eqs. (97) and (100) we can readily compute the quadratic condensates. For the $`\text{U}(N)`$ theory we find

$$\mathrm{Tr}\mathrm{\Phi }^2=N_0^1\varphi ^2(x)𝑑x=\frac{NM^2}{2\lambda ^2\mathrm{sin}\beta }\frac{\theta _1(\beta /2|\stackrel{~}{\tau })}{\theta _1^{}(\beta /2|\stackrel{~}{\tau })}+\frac{NM^2}{4\lambda ^2\mathrm{sin}^2\frac{\beta }{2}}=\frac{1}{\mu }W_{\mathrm{eff}}^{\text{U}(N)}.$$
(101)
For the $`\text{SU}(N)`$ theory we find similarly

$$\mathrm{Tr}\widehat{\mathrm{\Phi }}^2=N_0^1\widehat{\varphi }^2(x)𝑑x=\frac{NM^2}{2\lambda ^2\mathrm{sin}\beta }\frac{\theta _1^{}(\beta /2|\stackrel{~}{\tau })}{\theta _1(\beta /2|\stackrel{~}{\tau })}\frac{NM^2}{4\lambda ^2\mathrm{sin}^2\frac{\beta }{2}}=\frac{1}{\mu }W_{\mathrm{eff}}^{\text{SU}(N)}$$
(102)
Similarly higher condensates could be computed explicitly if required and will in general depend on the residue of $`g^n(x)`$ at the pole $`x=\beta /4\pi `$