text:
# Konishi anomaly approach to gravitational 𝐹-terms

## 1 Introduction

Recently much progress has been made in the computation and understanding of F-terms which describe the coupling of $`N=1`$ gauge theories to $`N=1`$ supergravity . In particular the original conjecture of Djikgraaf and Vafa , relating these $`F`$-terms to non-planar corrections to the free energy of a related matrix model has been proved by  using diagrammatic techniques, extending to the gravitational case the super-Feynman diagram techniques of  Crucial ingredient in the proof was the modification of the chiral ring relations due to the coupling of the gauge theory to supergravity. In particular, if one restricts to the first non-trivial gravitational $`F`$-term contribution, corresponding to the genus one correction in the related matrix model, one needs to take into account just the modification which follows from standard $`N=1`$ supergravity tensor calculus. The purpose of the present work is rather to study the problem from the viewpoint of generalized Konishi anomaly relations in the chiral ring, by extending to the case of $`N=1`$ gauge theories coupled to $`N=1`$ supergravity the strategy of  The key point in our analysis will be,together with the modification of the chiral ring mentioned above, the observation that in the presence of a non-trivial supergravity background the usual factorization property of chiral correlators does not hold. In particular connected two point functions are generically non-vanishing, much like in matrix models, where connected correlators receive a subleading $`1/N^2`$ contribution, in the $`1/N`$ expansion

The paper is organized as follows: in section 2 we study in detail the deformation of the chiral ring due to the coupling to supergravity and some of the resulting relations that we will need subsequently. In section 3 we estimate the chiral, connected two- and three-point functions that will enter later in the equations obtained using the generalized Konishi anomaly. Section 4 contains the main results: we derive there the generalized loop equations which enables us to solve for the relevant connected two-point functions. We show the uniqueness of the solution and use it to compute the first non-trivial gravitational correction to the effective superpotential.. We finally show that this agrees with the genus one correction to the matrix model free energy. In section 5 we make our conclusions and mention some open problems.

## 2 The chiral ring

A basic ingredient in deriving the effective $`F`$-terms in $`N=1`$ gauge theories after integrating out the adjoint matter is the chiral ring . All quantities in the chiral ring are defined modulo $`\overline{D}_{\dot{\alpha }}`$ exact terms, here $`\overline{D}_{\dot{\alpha }}`$ refers to the differential operator conjugate to the supersymmetric current $`\overline{Q}_{\dot{\alpha }}`$. It is sufficient to restrict the quantities to be in the chiral ring as $`\overline{D}_{\dot{\alpha }}`$ exact terms do not contribute to $`F`$ terms. For the $`N=1`$ gauge theory with adjoint matter on $`R^4`$ the chiral ring relations are given by

$$[W_\alpha ,\mathrm{\Phi }]=0\text{mod }\overline{D},\{W_\alpha ,W_\beta \}=0\text{mod }\overline{D},$$
(1)
where $`W_\alpha `$ is the $`N=1`$ gauge multiplet and $`\mathrm{\Phi }`$ is the chiral multiplet in the adjoint representation of the gauge group.

The chiral ring relations given above are modified in the presence of a curved background . We review here the derivation of the modified chiral ring in presence of background $`N=1`$ gravity. Consider the following $`\overline{D}`$ exact quantity

$$\{\overline{D}^{\dot{\alpha }},[D_{\alpha \dot{\alpha }},W_\beta ]\}=\{[\overline{D}^{\dot{\alpha }},D_{\alpha \dot{\alpha }}],W_\beta \}$$
(2)
where $`D_{\alpha \dot{\alpha }}`$ is the full covariant derivative containing the gauge field and the spin connection. We used the Jacobi identity and the fact that $`W_\beta `$ is chiral to obtaining the second term in the above equality. From the Bianchi identity  for covariant derivatives one has

$$[\overline{D}^{\dot{\alpha }},D_{\alpha \dot{\alpha }}]=4iW_\alpha 8iG_{\alpha \beta \gamma }M^{\beta \gamma }$$
(3)
where $`G_{\alpha \beta \gamma }`$ is the $`N=1`$ Weyl multiplet and $`M^{\alpha \beta }`$ refers to the Lorentz generator on spinors, its action on a spinor is given by

$$[M^{\alpha \beta },\psi _\gamma ]=\frac{1}{2}(\delta _\gamma ^\alpha \psi ^\beta +\delta _\gamma ^\beta \psi ^\alpha )$$
(4)
Substituting the Bianchi identity (3) in (2) and using the action of the Lorentz generator we find the deformed chiral ring given below

$`\{W_\alpha ,W_\beta \}=2G_{\alpha \beta \gamma }W^\gamma \text{ mod }\overline{D},\{W^\alpha ,\mathrm{\Phi }\}=0\text{ mod }\overline{D}.`$ (5)
The second equation is obtained by replacing the $`W_\beta `$ in (2) by $`\mathrm{\Phi }`$. For the conventional $`N=1`$ supergravity theory, in the first equation above, only the $`SU(N)`$ part of the gauge field $`W^\gamma `$ appears in the right hand side. In , the modification of the ring involved the $`U(1)`$ part of the gauge field as well, which corresponds to a non-standard $`N=1`$ supergravity theory relevant to D-brane gauge theories. In the present paper we will always be restricting ourselves to the standard $`N=1`$ supergravity. In order to avoid explicitly including the $`SU(N)`$ projectors in all the formulae below, we shall always take gauge field backgrounds to be in the $`SU(N)`$ part of $`U(N)`$.

Using the deformed chiral ring we can derive many identities valid in the ring which are used crucially in the next sections 111The conventions followed in this paper are given in Appendix A.. From the definition of $`W^2`$ and (5) we have

$`W_\alpha W_\beta `$ $`=`$ $`{\displaystyle \frac{ϵ_{\alpha \beta }}{2}}W^2+G_{\alpha \beta \gamma }W^\gamma ,`$ (6)
$`W_\alpha W_\beta W_\gamma `$ $`=`$ $`{\displaystyle \frac{1}{2}}ϵ_{\alpha \beta }W^2W_\gamma +{\displaystyle \frac{1}{2}}G_{\alpha \beta \gamma }W^2+G_{\alpha \beta \delta }G_{\gamma \sigma }^\delta W^\sigma ,`$
using the above identities we are led to the following relations

$`W_\alpha W^2`$ $`=`$ $`{\displaystyle \frac{1}{2}}W^2W_\beta {\displaystyle \frac{1}{2}}G^2W_\alpha ,`$ (7)
$`W^2W_\alpha `$ $`=`$ $`{\displaystyle \frac{1}{2}}W_\alpha W^2{\displaystyle \frac{1}{2}}G^2W_\alpha ,`$
$`[W^2,W_\alpha ]=0,`$ $`W^2W_\alpha ={\displaystyle \frac{1}{3}}G^2W_\alpha ,W^2W^2={\displaystyle \frac{1}{3}}G^2W^2,`$
$`(G^2)^2`$ $`=`$ $`G^{\alpha \beta \gamma }G_{\alpha \delta \sigma }G_{\beta \rho }^\sigma G_\gamma ^{\delta \rho }.`$
These identities imply that the gauge invariant combination of certain chiral operators vanish in the chiral ring. The following chiral operator vanishes in the chiral ring.

$$G_{\alpha \beta \gamma }\mathrm{Tr}(\mathrm{W}^\gamma \mathrm{\Phi }\mathrm{\Phi }\mathrm{})=0.\text{mod }\overline{\mathrm{D}}$$
(8)
It is clear that if there are no $`\mathrm{\Phi }`$’s in the trace, the above equation is true for the gauge group $`SU(N)`$. To proof the above identity for arbitrary number of $`\mathrm{\Phi }`$ we use the following equation

$`\mathrm{Tr}(\mathrm{W}_\alpha \mathrm{W}_\beta \mathrm{\Phi }\mathrm{\Phi }\mathrm{})`$ $`=`$ $`\mathrm{Tr}(\mathrm{W}_\beta \mathrm{W}_\alpha \mathrm{\Phi }\mathrm{\Phi }\mathrm{}),`$
$`=`$ $`{\displaystyle \frac{1}{2}}ϵ_{\alpha \beta }\mathrm{Tr}(\mathrm{W}^2\mathrm{\Phi }\mathrm{\Phi }\mathrm{}).`$
To obtain the first equation above we have used the cyclic property of trace and (5). Now multiplying the first equation of (6) with arbitrary number of $`\mathrm{\Phi }`$’s and using (2) we obtain (8). Multiplying (8) by $`G^{\alpha \beta }`$ and using (96) we obtain

$$G^2\mathrm{Tr}(\mathrm{W}_\alpha \mathrm{\Phi }\mathrm{\Phi }\mathrm{})=0,\text{mod }\overline{\mathrm{D}}.$$
(10)
Another important identity in the chiral ring is

$$G^4=(G^2)^2=0,\text{mod }\overline{D}$$
(11)
The proof goes along the same lines as the derivation of the deformed chiral ring. Consider the following $`\overline{D}`$ exact quantity

$$\{\overline{D}^{\dot{\alpha }},[D_{\alpha \dot{\alpha }},G_{\beta \gamma \delta }]\}=\{[\overline{D}^{\dot{\alpha }},D_{\alpha \dot{\alpha }}],G_{\beta \gamma \delta }\}$$
(12)
As $`G_{\beta \gamma \delta }`$ is uncharged with respect to gauge field the covariant derivative $`D_{\alpha ,\dot{\alpha }}`$ contains only the spin connection. The Bianchi identity for covariant derivatives now implies

$$[\overline{D}^{\dot{\alpha }},D_{\alpha \dot{\alpha }}]=8iG_{\alpha \beta \gamma }M^{\beta \gamma }$$
(13)
Substituting the above equation in (12) we obtain the following equation in the chiral ring

$$G_{\alpha \beta \sigma }G_{\gamma \delta }^\sigma +G_{\alpha \gamma \sigma }G_{\beta \delta }^\sigma +G_{\alpha \delta \sigma }G_{\beta \gamma }^\sigma =0$$
(14)
Multiplying this equation by $`G^{\alpha \gamma \rho }G_\rho ^{\beta \delta }`$ so that all the free indices are contracted and using the last equation in (7) and (96) we obtain (11). As a result the gravitational corrections to the $`F`$-terms truncate at order $`G^2`$.

## 3 Estimates of connected part of correlators

In the absence of gravity correlators of gauge invariant operators factorize in the chiral ring . This factorization enabled one to show that the loop equations satisfied by the resolvent on the gauge theory agreed with the loop equations of the matrix model in the large $`N`$ limit. On the matrix model side factorization in the loop equations was due to the the large $`N`$ limit. However the correspondence of the gauge theory with the matrix model proposed by Dijkgraaf and Vafa in  goes beyond the large $`N`$ approximation. If the $`N=1`$ gauge theory is placed in a background gravitational field then the gravitational corrections to the $`F`$ terms of the gauge theory is of the form

$$\mathrm{\Gamma }_1=d^4xd^2\theta G^2F_1(S),$$
(15)
here $`S`$ refers to the gaugino condensate. The Dijkgraaf and Vafa proposal states that the one can calculate $`F_1`$ from the genus one diagrams of the corresponding matrix model. It has been shown in  that the modification of the chiral ring in the presence of gravity allows the computation of the gravitational corrections to the $`F`$-term. The gravitational corrections enter at genus one on the gauge theory side and they reduce to the genus one diagrams of the corresponding matrix model. This implies that the loop equations of the gauge theory would not factorize in presence of gravity, as the loop equations of the matrix model do not factorize at genus one. Therefore a priori we expect that gauge invariant operators do not factorize in the presence of gravity. In this section we estimate the gravitational corrections to the connected parts of various correlators that can appear in the loop, equations using the deformation of the chiral ring in presence of gravity (5).

To estimate the gravitational corrections to the $`F`$-term obtained by integrating out the chiral multiplet we require the two point function the chiral scalars $`\mathrm{\Phi }\mathrm{\Phi }`$. This is obtained after integrating out the antichiral scalar and it is given by

$$\mathrm{\Phi }(x,\theta )\mathrm{\Phi }(x^{},\theta ^{})=\frac{\overline{m}}{\mathrm{}+\frac{i}{2}D^\alpha W_\alpha +iW^\alpha D_\alpha +m\overline{m}}\delta ^4(xx^{})\delta ^2(\theta \theta ^{})$$
(16)
This propagator in the presence of a gravitational background was derived in  and in writing the above equation we have assumed that the gravitational background is on shell which allows one to set the other terms which occur in the propagator to zero, $`\mathrm{}`$ is (16) stands of the full gauge and gravitational covariant derivative. The action of $`W_\alpha `$ is dictated by the representation of $`\mathrm{\Phi }`$, in this paper we will restrict ourselves to the adjoint action. The delta functions in (16) refer to the full covariant delta function in curved superspace. In order to obtain the leading estimates for the connected component of various correlators due to the modification of the chiral ring it is sufficient to use the free d’ Alembertian operator and a constant gaugino background. We argue this as follows, we can expand the propagator of (16) in a weak field as

$`{\displaystyle \frac{\overline{m}}{\mathrm{}+iD^\alpha W_\alpha +m\overline{m}}}`$ $`=`$ $`{\displaystyle \frac{\overline{m}}{\mathrm{}_0+iW^\alpha D_\alpha +m\overline{m}}}`$
$`+`$ $`{\displaystyle \frac{\overline{m}}{\mathrm{}_0+iW^\alpha D_\alpha +m\overline{m}}}(\mathrm{}\mathrm{}_0){\displaystyle \frac{1}{\mathrm{}_0+iW^\alpha D_\alpha +m\overline{m}}}+\mathrm{}`$
here $`\mathrm{}_0`$ refers to the free d’Alembertian operator. We have also dropped the terms $`D^\alpha W_\alpha `$ in (16) as we have considered a covaraintly constant gaugino background. From the expansion we see that the corrections in using the free d’ Alembertian operator in the propagator always occur with the factor $`(\mathrm{}\mathrm{}_0)`$, which is proportional to the gravitational background and therefore subleading to the estimate obtained using the free d’ Alembertian operator in the first term. One has to make a similar expansion for the covariant delta function in (16) and again one can see corrections in using the flat space delta function are subleading. However we will see later that if there is no gravitational contribution form the deformed chiral ring the leading estimate for the connect part of correlators arise from corrections in the propagator due to the presence of the full covariant $`\mathrm{}`$ and the covariant delta function. For the free d’Alembertian operator in the propagator it is possible to go over to momentum space and to write a Schwinger parameterization of the propagator as follows.

$$\mathrm{\Phi }(x,\theta )\mathrm{\Phi }(x^{},\theta )=𝑑sd^4pd^2\pi e^{ip(xx^{})}e^{\frac{s}{\overline{m}}(p^2+W^\alpha \pi _\alpha +m\overline{m})}$$
(18)
where $`\pi ^\alpha =iD^\alpha `$. In the above equation we have restricted to the superspace variable $`\theta `$ to be the same at $`x`$ and $`x^{}`$ as we will be interested in correlators at the same point in the superspace variable $`\theta `$. We will now use this propagator and the modified chiral ring (5) to make estimates for the connected part of various correlators. The modified ring allows more than two insertions of $`W_\alpha `$ in a given index loop, using the identities in (7) such contributions can be converted to gravitational corrections. At this point one might wonder if contributions to the connected diagrams of gauge invariant operators in presence of a gravitational background are in contradiction with the result found in . There it was found that on an arbitrary Kähler manifold gauge invariant operators of the $`N=1`$ theory factorize. The background considered in  was entirely bosonic, we find the estimates of contribution to the connected diagram indeed vanish for a purely bosonic background, thus there is no contradiction with . We will indicate this as we evaluate the estimates of various correlators.

The various operators involved in the correlators of interest are

$`(z)_{ij}={\displaystyle \frac{1}{32\pi ^2}}\left({\displaystyle \frac{W^2}{z\mathrm{\Phi }}}\right)_{ij},`$ $`R(z)=\mathrm{Tr}(z)`$
$`\rho _{\alpha ij}(z)={\displaystyle \frac{1}{4\pi }}\left({\displaystyle \frac{W_\alpha }{z\mathrm{\Phi }}}\right)_{ij},`$ $`w_\alpha (z)=\mathrm{Tr}\rho _\alpha (z)`$
$`𝒯(z)_{ij}=\left({\displaystyle \frac{1}{z\mathrm{\Phi }}}\right)_{ij},`$ $`T(z)=\mathrm{Tr}𝒯(z)`$ (19)
here we have defined separate symbols for the matrix elements and the trace for later convenience, the gauge invariant operators we will consider are the ones with the trace in the above equation. The contour integrals of $`R`$, $`w_\alpha `$ and $`T`$ around $`i`$-th branch cut define the gaugino bilinear $`S_i`$, the $`U(1)`$ gauge field $`w_{\alpha i}`$ and $`N_i`$ respectively in $`U(N_i)`$ subgroup of $`U(N)`$ as in . The fact that we are here restricting the background gauge field to be in $`SU(N)`$ rather than $`U(N)`$ implies that $`_iw_{\alpha i}=0`$. The chiral ring relation $`G^2w_\alpha (z)=0`$ implies that $`G^2w_{\alpha i}=0`$ for all $`i`$.

We first consider estimates of the connected part of two point function, we will discuss in detail the estimate for the follow correlator

$$R(z,x,\theta ),R(w,y,\theta )_c,$$
(20)
where the subscript stands for the connected part and out line the derivation of the estimates for the other two point functions. The various contribution to this correlator in (20) can be found by expanding in $`z`$ and $`w`$, by definition of the connected correlator the expansion starts of with the power $`1/z^2w^2`$. Let us focus on a Feynman diagram consisting of $`l`$ loops, there will be $`l+1`$ bosonic and fermionic momentum integrations in this diagram. The extra momentum integral comes from the final Fourier transform to convert to the position space representation for the above correlator. We can organize this diagram into index loops due to the adjoint action of $`W_\alpha `$. The fermionic momentum integral forces us to bring down $`2(l+1)`$ powers of $`W`$ from the propagator in (18). To obtain the leading gravitational correction we would like the number of index loops to be as large as possible so that we can avoid having more than two $`W`$’s in a given index loop. The number of index loops $`h`$ and the number of loops are related by $`l=h1+2g`$, where $`g`$ is the genus of the diagram. For a given number of Feynman loops the number of index loops is largest for genus zero, thus the leading estimate to the connected graph arises from the planar diagram. For a planar diagram we need to saturate the fermionic momentum integrals by bringing down $`2h`$ $`W`$’s. This can be done by inserting $`W^2`$ in $`h`$ index loops. We still have two more external $`W^2`$ in (20). This can at best be inserted in two different index loops. Thus we have two index loops with $`(W^2)^2`$ insertions. Using the identities in (7) we see that each of them reduces to $`G^2W^2`$. Thus there is a term proportional to $`G^4W^2`$ on one of the index loops. Note that if one had a purely bosonic background both $`G`$ and $`W`$ would start at $`\theta `$ in the superspace expansion, thus $`G^4W^2`$ would vanish in agreement with . In fact $`G^4W^2`$, as is trivial in the chiral ring by (11), the leading estimate for the correlator in (20) vanishes. In the next section it is shown that the two point function in (20) in fact vanishes. Next we consider the following two point function

$$R(z,x,\theta )w_\alpha (w,x^{},\theta )_c,$$
(21)
we can apply the same counting again, finally in the planar diagram we will be left with at best with one index loop with $`(W^2)^2`$ and one with $`W^2W_\alpha `$ insertions. This reduces to a $`G^2W^2`$ insertion and a $`G^2W_\alpha `$ insertions, which implies that the leading gravitational contribution to (21) is proportional to $`G^4`$. Now we have seen in (10) that $`G^2W_\alpha `$ is zero in the chiral ring, thus this leading estimate in fact vanishes in the chiral ring. Similarly, consider the correlator

$$R(z,x,\theta )T(w,x^{},\theta )_c.$$
(22)
Here we will be left with $`(W^2)^2`$ in a single index loop, which reduces to $`G^2W^2`$. Thus the above correlator is proportional to $`G^2W^2`$. For a purely bosonic background we see that this contribution again vanishes, consistently with . For the case of

$$w^\alpha (z,x,\theta )w_\alpha (w,x^{},\theta )_c,$$
(23)
we will be left with either $`W^\alpha W^2`$ insertion in two different index loops or a $`(W^2)^2`$ insertion in a single index loop. The former case vanishes in the chiral ring, but the latter case survives, with a contribution proportional to $`G^2W^2`$. For the following two point function

$$w_\alpha (z,x,\theta )T(w,x^{},\theta )_c,$$
(24)
there is a $`W^\alpha W^2`$ insertion in a single index loop, which is proportional to $`G^2W_\alpha `$. Note that this leading contribution vanishes in the chiral ring due to (10) and also for a purely bosonic background. Finally, we have the two point function

$$T(z,x,\theta )T(w,x^{},\theta )_c$$
(25)
For this case all the $`h`$ index loops are saturated with one $`W^2`$ and there are extra insertions of $`W^2`$ for any of the $`h`$ index loops, as there is no external $`W`$. Thus we have no contribution for the connected part of this correlator from the modified ring. However we will see later that there is a direct gravitational contribution to the above correlator. This can be seen roughly as follows: the d’ Alembertian in (16) carries the covariant derivatives which can possibly contribute to the connected two point function, as seen in the expansion of the propagator in (3). This fact can be further justified by the evaluation of the 1-loop effective action obtained by integrating out the chiral multiplet in the absence of the gauge field background, which gives a term proportional to $`G^2\mathrm{ln}(m)`$ . Therefore we expect the leading term in the correlator in (25) to be proportional to $`G^2`$ and this will be shown explicitly in the next section. Again we see that for a purely bosonic background $`G^2`$ is proportional to $`\theta ^2`$, which implies that the lowest components of the superfields in (25) factorize consistently with  222Note that the fact that the lowest component factorize cannot be used to promote it to a superfield equation as the is no $`Q_\alpha `$ which preserves the background.. We summarize the estimates of the various connected two point correlators in the following table for future reference

| Correlator | Estimate |
| --- | --- |
| $`RR_c`$ | $`G^4=0`$ |
| $`Rw_\alpha _c`$ | $`G^4=0`$ |
| $`RT_c`$ | $`G^2S^h`$ |
| $`w^\alpha w_\alpha _c`$ | $`G^2S^h`$ |
| $`w_\alpha T_c`$ | $`G^2S^{h1}w_{\alpha i}=0`$ |
| $`TT_c`$ | $`G^2S^{h1}`$ |

Table 1.

where $`S`$ represents schematically any of the $`S_i`$’s, and we have used the chiral ring relations $`G^4=0`$ and $`G^2w_{\alpha i}=0`$.

Now we provide estimates for the fully connected part of various three point functions. All the fully connected part of the three point functions are proportional to at least $`G^4`$, therefore using (11) they all vanish in the chiral ring. We discuss the method of arriving at the estimates for one case in detail and just outline the results for the others. Consider the following fully connected three point function.

$$R(z,x,\theta )R(w,x^{},\theta )R(u,x^{\prime \prime },\theta )_c.$$
(26)
By the definition of the full connected three point function, the first possibly non zero term of the expansion in $`z,w,u`$ starts at $`1/(zwu)^2`$. Consider a contribution to any of the correlators appearing in this expansion. A Feynman diagram consisting of $`l`$ loops will now have $`l+2`$ bosonic and fermionic momentum integrations. This is because a three point function in momentum space will in general have two external independent momenta, and then converting that to position space will involve these additional momentum integrals. As we have argued earlier for the case of the two point function, the leading contribution will be from the genus zero graphs. For a planar graph then there are $`2(h+1)`$ fermionic momentum integrals to be done. Therefore in addition to inserting $`h`$ index loops by $`W^2`$, at best four different index loops will have $`(W^2)^2`$ insertions (we assume $`l`$ is large enough). Using the identities in (7) we see that each $`(W^2)^2`$ insertion is proportional to $`G^2W^2`$. Thus the above three point function is proportional to $`G^8`$, but $`G`$ is fermionic and has $`4`$ independent components, thus $`G^8`$ vanishes due to Fermi statistics. Similar arguments show that the following correlators vanish

$`R(z,x,\theta )R(w,x^{},\theta )w_\alpha (u,x^{\prime \prime },\theta )_c=0,`$ (27)
$`R(z,x,\theta )R(w,x^{},\theta )T(u,x^{\prime \prime },\theta )_c=0,`$
$`R(z,x,\theta )w^\alpha (w,x^{},\theta )w_\alpha (u,x^{\prime \prime },\theta )_c=0,`$
$`R(z,x,\theta )w^\alpha (w,x^{},\theta )T(u,x^{\prime \prime },\theta )_c=0,`$
The first correlator in (27) is proportional to $`G^8`$ and the rest are proportional to $`G^6`$, thus they vanish due to Fermi statistics. Now consider the following three point function

$$w^\alpha (z,x,\theta )T(w,x^{},\theta )T(u,x^{\prime \prime },\theta )_c.$$

We have seen that in a planar graph the fermionic momentum integrations force one to insert at least one factor of $`W^2`$ in all of the $`h`$ index loops and there is at least one index loop with a $`(W^2)^2`$ insertion. Using the identities in (7) this can be manipulated to a gravitational contribution proportional to $`G^2W^2`$. For the above correlator there is one more index loop with an insertion of $`W^2W^\alpha `$ and again using the identities in (7), this term is proportional to $`G^2W^\alpha `$. Thus the leading gravitational contribution to the three point function in (3) is proportional to $`G^4W_\alpha `$ and thus it is zero in the chiral ring using (10). Next we consider the following three point function

$$R(z,x,\theta )T(w,x^{},\theta )T(u,x^{\prime \prime },\theta )_c,$$

As discussed above, since there are $`2(h+1)`$ momentum integration, all the $`h`$ index loops have at least a $`W^2`$ insertion with one having a $`(W^2)^2`$ insertions, the above correlator has an external $`W^2`$. Therefore the leading gravitational contribution arises with two different index loops each with a $`(W^2)^2`$ insertion and using the identities in (7) this can be manipulated in the chiral ring to give a factor of $`G^4`$. The following three point function is also proportional to $`G^4`$

$$W^\alpha (z,x,\theta )W_\alpha (w,x^{},\theta )T(u,x^{\prime \prime },\theta )_c,$$

Here again due to the momentum integrations there is already a factor of $`G^2`$, the two external $`W`$’s can be inserted in another index loop, but this loop already has an insertion of $`W^2`$, which gives rise to a contribution proportional to $`G^2`$. Thus the three point function in (3) is proportional to $`G^4`$. Therefore the correlators in (3) and (3) vanish in the chiral ring due to (11). Finally, let us consider the following three point function

$$T(z,x,\theta )T(w,x^{},\theta )T(u,x^{\prime \prime },\theta )_c.$$

As in the case for (25) we can not estimate the $`G`$ dependence of this correlator solely using the chiral ring. But from the fact that the correlator in (25) is proportional to $`G^2`$ and since there is at least one index loop with $`(W^2)^2`$ insertion, we can arrive at the conclusion that the above three point function will be proportional to $`G^4`$. To sum up we have examined all possibly non-vanishing full connected three point fully and found them to be least proportional to $`G^4`$, and therefore using (11) they vanish in the chiral ring.

## 4 Anomaly equations and matrix model loop equations

In this section we will use the generalized Konishi anomaly to extract the gravitational corrections to the effective superpotential. We recall that the anomaly equations obtained in  in the absence of gravitational fields are as follows

$`R(z)R(z)`$ $``$ $`\mathrm{tr}(V^{}(\mathrm{\Phi })(z))=0`$
$`2R(z)w_\alpha (z)`$ $``$ $`\mathrm{tr}(V^{}(\mathrm{\Phi })\rho _\alpha (z)=0`$
$`2R(z)T(z)`$ $``$ $`\mathrm{tr}(V^{}(\mathrm{\Phi })𝒯(z)+w^\alpha (z)w_\alpha (z)=0`$ (28)
Here $`V`$ denotes the classical superpotential of degree $`n+1`$. We have indicated above the full two point functions that include the disconnected and connected two point functions. The latter vanish in the absence of gravitational field but as we will show below do not vanish in the presence of the gravitational field. These equations were obtained by the generalized Konishi anomalies upon transforming the adjoint chiral field $`\mathrm{\Phi }`$ as $`\delta \mathrm{\Phi }_{ij}`$ equal to $`_{ij}(z)`$, $`\eta ^\alpha \rho _{\alpha ij}(z)`$ and $`𝒯_{ij}(z)`$ respectively, with $`\eta ^\alpha `$ being an arbitrary field independent spinor. In general for the infinitesimal transformation $`\delta \mathrm{\Phi }_{ij}=f_{ij}`$, the generalized Konishi anomaly is given by

$$\frac{\delta f_{ji}}{\delta \mathrm{\Phi }_k\mathrm{}}A_{ij,k\mathrm{}}$$
(29)
In the absence of gravitation

$$A_{ij,k\mathrm{}}=(W^2)_{kj}\delta _i\mathrm{}+\delta _{kj}(W^2)_i\mathrm{}2W_{kj}^\alpha W_{\alpha i\mathrm{}}$$
(30)
Using the above anomaly and the equation $`\{W_\alpha W_\beta \}=0`$ in the chiral ring, one obtains the equation (28).

In the presence of the gravitational field $`G_{\alpha \beta \gamma }`$, these equations are modified for two reasons: firstly there is a direct gravitational anomaly (ie. even in the absence of gauge fields, in other words when chiral multiplets couple only to gravitational fields) and secondly due to the modification of the ring (5).

Under the infinitesimal transformation $`\delta \mathrm{\Phi }_{ij}=f_{ij}`$ the pure gravitational contribution to the anomaly is

$$\frac{\delta f_{ji}}{\delta \mathrm{\Phi }_k\mathrm{}}\alpha G^2\delta _{kj}\delta _i\mathrm{}$$
(31)
where $`\alpha =\frac{1}{3}`$ as shown in the appendix B, and from now on $`G^2`$ will denote $`\frac{1}{32\pi ^2}G^{\alpha \beta \gamma }G_{\alpha \beta \gamma }`$. This changes

$$A_{ij,k\mathrm{}}A_{ij,k\mathrm{}}+\frac{1}{3}G^2\delta _{kj}\delta _i\mathrm{}$$
(32)
It is easy to see that the pure gravitational anomaly and the modification of the chiral ring (5) together with the consequent identities given in the equations (7),(10),(11), give rise to the following modification of the equations (28):

$`R(z)R(z)`$ $``$ $`tr(V^{}(\mathrm{\Phi })(z))=0`$
$`2R(z)w_\alpha (z)`$ $``$ $`tr(V^{}(\mathrm{\Phi })W_\alpha (z))=0`$
$`2R(z)T(z)`$ $``$ $`tr(V^{}(\mathrm{\Phi })𝒯(z))`$ (33)
$`+`$ $`w^\alpha (z)w_\alpha (z){\displaystyle \frac{1}{3}}G^2T(z)T(z)=0`$
Note that in the first two equations above the pure gravitational anomaly cancels with the contributions coming from the modification of the chiral ring via eq(7). This happens due to the remarkable fact that the pure gravitational anomaly (32) comes with a factor of $`\frac{1}{3}`$ !! We have also used the fact that $`G^2w_\alpha `$ vanishes in the chiral ring. In the last equation however there is no contribution due to the modification of the ring and hence the last two terms on the left hand side arise solely from the pure gravitational anomaly. Here again, a priori, the two point functions are the sum of connected and disconnected parts.

For later purposes let us write the first equation in (33) more explicitly

$$R(z)^2V^{}(z)R(z)\frac{1}{4}f(z)=R(z)R(z)_c$$
(34)
where $`f(z)`$ is a polynomial of degree $`n1`$ defined by

$$\mathrm{tr}(V^{}(\mathrm{\Phi })V^{}(z)(z)=\frac{1}{4}f(z)$$
(35)
and the subscript $`c`$ denotes the connected part of the correlation function. The latter, as we have seen in the last section, goes as $`G^4`$ and therefore trivial in the chiral ring. As a result the equation for $`R`$ is unmodified by the gravitational field. Since the finite polynomial $`f`$ is determined completely by the periods of $`R`$, i.e. the contour integrals around the various branch cuts $`C_i`$, $`i=1,\mathrm{},n`$,

$$\frac{1}{2\pi i}_{C_i}𝑑zR(z)=S_i$$
(36)
we conclude that $`R`$ does not receive any gravitational corrections.

The strategy now is to expand all of the quantities appearing above in a perturbation series in $`G^2`$. Of course, this series ends at order $`G^2`$ since $`G^4`$ is trivial. Thus for example we write

$$T=T^{(0)}+G^2T^{(1)}$$
(37)
and similarly for the connected parts of the 2-point functions appearing in the above equations. As discussed in the last section, the latter start from order $`G^2`$, and therefore the equations for $`T^{(0)}`$ and $`w_\alpha ^{(0)}`$ are the same as in . To go beyond this we need to solve for the connected two point functions.

### 4.1 Equations for the connected two point functions in the presence of gravitational fields

We will now derive equations for the connected two point functions that appear in eq.(33). Although, in eq.(33) the connected 2-point functions are of the form $`R(z)T(z)`$, i.e. both the operators are at same $`z`$, it turns out to be more convenient to consider the two operators at different points in the complex plane (say $`z`$ and $`w`$). The reason is that we can impose conditions on a connected 2-point function, like $`R(z)T(w)`$ that their integrals around various branch cuts in $`z`$ and $`w`$ plane vanish separately. As a result, we will be able to solve completely the corresponding generalized Konishi anomaly equations.

We illustrate the general method of obtaining the generalized Konishi anomaly equations for the connected 2-point functions in one example and then give the complete set of equations which can easily be derived following the methods given below.

Consider the infinitesimal transformation (local in superspace coordinates $`(x^\mu ,\theta ,\overline{\theta })`$)

$$\delta \mathrm{\Phi }_{ij}=_{ij}(z)T(w)$$
(38)
The Jacobian of this transformation has two parts

$$\frac{\delta (\delta \mathrm{\Phi }_{ji})}{\delta \mathrm{\Phi }_{kl}}=\frac{\delta R_{ji}(z)}{\delta \mathrm{\Phi }_{kl}}T(w)+\underset{m}{}R_{ji}(z)T_{mk}(w)T_{lm}(w)$$
(39)
The first term in the equation above together with the variation of the classical superpotential gives rise to

$`(R(z)R(z)`$ $``$ $`tr(V^{}(\mathrm{\Phi })(z)))T(w)`$ (40)
$`=(R(z)R(z)tr(V^{}(\mathrm{\Phi })(z)))T(w)+2R(z)R(z)T(w)_c`$
$`tr(V^{}(\mathrm{\Phi })(z)))T(w)_c+R(z)R(z)T(w)_c`$
where the subscript $`c`$ denotes completely connected 2 or 3 point functions as indicated. The first term on the right hand side vanishes by virtue of the first equation of (33).

The second term in the Jacobian when combined with the anomaly (30,32) gives rise to a single trace contribution

$$\frac{1}{3}G^2tr((z)𝒯(w)𝒯(w))=\frac{1}{3}G^2_w\frac{R(z)R(w)}{zw}$$
(41)
Combining eqs.(40),(41) and the first equation of (33), one obtains the following equation for the connected correlation functions:

$$(2R(z)I(z))R(z)T(w)_c+R(z)R(z)T(w)_c\frac{1}{3}G^2_w\frac{R(z)R(w)}{zw}=0$$
(42)
Here the integral operator $`I(z)`$ denotes the following:

$$I(z)A(z)=\frac{1}{2\pi i}_{C_z}𝑑y\frac{V^{}(y)A(y)}{yz}$$
(43)
with the contour $`C_z`$ encircling $`z`$ and $`\mathrm{}`$. It is clear that for $`A`$ equal to $``$, $`\rho _\alpha `$ or $`𝒯`$, the integral operator reduces to

$$I(z)A(z)=V^{}(\mathrm{\Phi })A(z)$$
(44)
as is familiar in the matrix model works.

Since the last term in the eq.(42) is of order $`G^2`$ and involves one point function of $`R`$ which is certainly not zero, we conclude that the sum of the remaining terms which involve connected correlations functions cannot all vanish at order $`G^2`$. This proves our basic assertion. In fact the connected 3-pt. function $`R(z)R(z)T(w)_c`$ vanishes, as argued in eq.(27) in the last section, so eq.(42) implies that the connected 2-pt function $`R(z)T(w)_c`$ does not vanish at order $`G^2`$.

In order to completely solve the relevant connected correlation functions, we need to consider all transformations of the form

$$\delta \mathrm{\Phi }_{ij}=A_{ij}(z)B(w)$$
(45)
where $`A`$ is $``$, $`\rho _\alpha `$ or $`𝒯`$ and $`B`$ is $`R`$, $`w_\beta `$ or $`T`$. The resulting generalized Konishi anomaly equations can be derived in the same way as above and can be summarized in the following matrix equation:

$`\left[\begin{array}{ccc}\widehat{M}(z)& 2T(z)& 0\\ 0& M(z)& 0\\ 0& 0& M(z)\end{array}\right]`$ $`\left[\begin{array}{ccc}T(z)T(w)_c& T(z)R(w)_c& T(z)w_\beta (w)_c\\ R(z)T(w)_c& R(z)R(w)_c& R(z)w_\beta (w)_c\\ w_\alpha (z)T(w)_c& w_\alpha (z)R(w)_c& w_\alpha (z)w_\beta (w)_c\end{array}\right]=`$ (52)
$`={\displaystyle \frac{1}{3}}G^2_w\left[\begin{array}{ccc}T(z,w)& R(z,w)& 0\\ R(z,w)& 0& 0\\ 0& 0& 5ϵ_{\alpha \beta }R(z,w)\end{array}\right]`$ (56)
Here we have used the chiral ring equations $`G^2w_\alpha =0`$ and $`(G^2)^2=0`$ as shown in the section 2. We have also dropped various connected 3-pt functions that vanish via eq.(27). $`M(z)`$ denotes the integral operator $`(2R(z)I(z))`$, $`\widehat{M}(z)`$ denotes $`M(z)\frac{2}{3}G^2T(z)`$ and finally $`R(z,w)`$ and $`T(z,w)`$ denote $`(R(z)R(w))/(zw)`$ and $`(T(z)T(w))/(zw)`$ respectively.

There are a few points to note about these equations:
1) Chiral ring equations are consistent with the above matrix equation. For example, if one takes the equation for $`w_\alpha (z)B(w)_c`$ and multiplies by $`G^2`$ one finds that the equation identically vanishes in the chiral ring.
2) We have used the estimates given in the last section only to drop all the completely connected 3-pt. functions as they were shown to go as $`G^4`$. The estimates for the connected 2-pt. functions given in the Table 1, says that $`R(z)R(w)_c`$, $`R(z)w_\alpha (w)_c`$, and $`T(z)w_\alpha (w)_c`$ all vanish in the chiral ring. We have not used these estimates in the above equation but we note that the (1,3), (2,2), (2,3), (3,1) and (3,2) matrix elements on the right hand side vanish. Thus the estimates given in the table are indeed consistent with the matrix equation above. In fact, in the next subsection, we will show that the solutions to this equation are unique thereby proving that these connected 2-pt. functions vanish. Similarly had we used the estimate for the 2-pt. function $`T(z)T(w)_c`$ which was shown in the last section to go as $`G^2`$, we could have replaced $`\widehat{M}`$ by $`M`$ in the above equation.
3) The integrability condition is satisfied: the above equation is of the form

$$(z)N(z,w)=_wK(z,w)$$
(57)
where $`(z)`$ is the first matrix operator appearing on the left hand side of eq(56), $`N(z,w)`$ and $`K(z,w)`$ satisfy $`N(z,w)=N^t(w,z)`$ and $`K(z,w)=K^t(w,z)`$. The non-trivial consistency condition then is

$$(_wK(z,w))^t(w)=(z)_zK(z,w)$$
(58)
The crucial identity needed for this is the one involving the integral operator $`I(z)`$ and is as follows:

$$(I(z)_zI(w)_w)\frac{A(z)A(w)}{zw}=(_z_w)\frac{I(z)A(z)I(w)A(w)}{zw}$$
(59)
for any function $`A`$ which is smooth at $`z`$ and $`w`$. This can be proved by using the definition of the contours involved in $`I(z)`$ and $`I(w)`$. It follows that the (1,2), (2,1) and (3,3) components of the integrability condition (58) implies the following equation:

$$G^2(R(z)^2I(z)R(z))=0$$
(60)
This equation is just $`G^2`$ times the first equation of (33) if one takes into account the fact that the connected part of the correlation function appearing in the latter already is of order $`(G^2)^2`$. The only other non-trivial part of the integrability condition is its (1,1) component:

$$G^2[(2R(z)I(z))T(z)\frac{1}{3}G^2T(z)^2]=0$$
(61)
which is just $`G^2`$ times the disconnected part of the third equation of (33) thereby proving the integrability condition for the (1,1) component. This is because all the connected parts appearing in that equation will be trivial when multiplied by $`G^2`$. Note also that the last term in eq.(61) could have been dropped as it is trivial. Its origin comes from the extra term in $`\widehat{M}`$ appearing in the (1,1) component of $``$ which as argued in the point 2) above could have been dropped.

### 4.2 Uniqueness of the solutions for the connected two point functions

Since the integrability conditions are satisfied, solution to eq.(56) exists. However the solution has a finite ambiguity which will be fixed by the physical requirement that the contour integrals around all the branch cuts of the connected two point functions in the the $`z`$ and $`w`$ planes must vanish separately. The reason for this is that the following operator equations hold:

$$\frac{1}{2\pi i}_{C_i}𝑑zR(z)=S_i,\frac{1}{2\pi i}_{C_i}𝑑wT(w)=N_i,\frac{1}{2\pi i}_{C_i}𝑑zw_\alpha (z)=w_{\alpha i}.$$
(62)
where $`S_i`$ is the chiral superfield whose lowest component is the gaugino bilinear in the $`i`$-th gauge group factor in the broken phase $`U(N)_{i=1}^nU(N_i)`$ and $`w_{\alpha i}`$ is the $`U(1)`$ chiral gauge superfield of the $`U(N_i)`$ subgroup. Since these fields are background fields, in the connected correlation functions the contour integrals around the branch cuts must vanish. Similarly since these background fields are independent of the gravitational fields, order $`G^2`$ corrections to the one point functions of $`R`$, $`w_\alpha `$ and $`T`$ must also have vanishing contour integrals around the branch cuts.

For later use, we can write a complete set of normalized differentials using eq. (34) upto order $`G^2`$ as

$$\omega _j=\frac{1}{4}\frac{dz}{[2R(z)V^{}(z)]}\frac{}{S_j}f(z),\frac{1}{2\pi i}_{C_i}\omega _j=\delta _{ij}$$
(63)
To illustrate the method, we will again focus on $`R(z)T(w)_c`$. The action of operator $`I(z)`$ is given as:

$$I(z)R(z)T(w)_c=V^{}(z)R(z)T(w)_c+\underset{k=0}{\overset{n1}{}}c_k(w)z^k$$
(64)
where the first term on the right hand side comes from the contour integral around $`z`$ and the second term from that around $`\mathrm{}`$, with $`n`$ being the order of $`V^{}(z)`$. Here we have used the fact that $`R(z)`$ asymptotically vanishes as $`1/z`$ 333 Actually the coefficient of $`1/z`$ is $`\mathrm{TrW}^2`$ and hence in the connected 2-pt function it vanishes. As a result the connected 2-pt function goes as $`1/z^2`$ (and for similar reasons $`1/w^2`$ asymptotically in $`w`$) which means that the sum over $`k`$ is between $`0`$ and $`n2`$. However in the above expression we have kept the sum upto $`n1`$ since, as it will turn out, the condition of vanishing contour integrals around all the branch cuts will in particular imply that $`c_{n1}=0`$. The coefficients $`c_k(w)`$ are arbitrary functions of $`w`$ which asymptotically vanish as $`1/w^2`$. Similarly we have

$$I(w)R(z)T(w)_c=V^{}(w)R(z)T(w)_c+\underset{k=0}{\overset{n1}{}}\stackrel{~}{c}_k(z)w^k$$
(65)
with $`\stackrel{~}{c}_k(z)`$ being arbitrary functions of $`z`$ that vanish asymptotically as $`1/z^2`$.

From (56), the two equations that this two point function satisfies are as follows:

$`(2R(z)I(z))R(z)T(w)_c`$ $`=`$ $`{\displaystyle \frac{1}{3}}G^2_wR(z,w)`$
$`(2R(w)I(w))R(z)T(w)_c`$ $`=`$ $`{\displaystyle \frac{1}{3}}G^2_zR(z,w)`$ (66)
The solutions to these two equations are

$`R(z)T(w)_c`$ $`=`$ $`{\displaystyle \frac{1}{2R(z)V^{}(z)}}[{\displaystyle \frac{1}{3}}G^2_wR(z,w)+{\displaystyle \underset{k=0}{\overset{n1}{}}}c_k(w)z^k]`$ (67)
$`=`$ $`{\displaystyle \frac{1}{2R(w)V^{}(w)}}[{\displaystyle \frac{1}{3}}G^2_zR(z,w)+{\displaystyle \underset{k=0}{\overset{n1}{}}}\stackrel{~}{c}_k(z)w^k]`$
Equating the two right hand sides , we see that $`c_k`$ and $`\stackrel{~}{c}_k`$ are not arbitrary functions of the respective arguments but are fixed upto a finite polynomial ambiguity in $`z`$ as well as $`w`$. They must be of the form

$`{\displaystyle \underset{k=0}{\overset{n1}{}}}c_k(w)z^k`$ $`=`$ $`{\displaystyle \frac{1}{(2R(w)V^{}(w))}}{\displaystyle \frac{G^2}{3}}[R(w){\displaystyle \frac{(V^{}(w)V^{}(z)+(zw)V^{\prime \prime }(w))}{(zw)^2}}`$
$`+{\displaystyle \frac{1}{4}}{\displaystyle \frac{(f(w)f(z)+(zw)f^{}(w))}{(zw)^2}}+{\displaystyle \underset{k,\mathrm{}=0}{\overset{n1}{}}}c_k\mathrm{}z^kw^{\mathrm{}}]`$
$`\stackrel{~}{c}_k(z)`$ $`=`$ $`c_k(z)+{\displaystyle \underset{\mathrm{}}{}}(c_k\mathrm{}c_\mathrm{}k)z^{\mathrm{}}`$ (68)
where $`c_k\mathrm{}`$ are arbitrary coefficients to be determined later. In deriving the above equation we have used the equation (34) with the right hand side set to zero. This is because the correction coming from the right hand side is of order $`G^4`$ and hence trivial.

Substituting this expression in eq.(67) and repeatedly using eq.(34), we obtain after some algebra:

$`R(z)T(w)_c`$ $`=`$ $`{\displaystyle \frac{G^2}{3}}{\displaystyle \frac{1}{[2R(z)V^{}(z)][2R(w)V^{}(w)]}}[{\displaystyle \underset{k,\mathrm{}=0}{\overset{n1}{}}}c_k\mathrm{}z^kw^{\mathrm{}}+`$ (69)
$`{\displaystyle \frac{[R(z)R(w)V^{}(z)R(w)\frac{1}{4}f(z)+(zw)]}{(zw)^2}}]`$
Note that the second term in the bracket is symmetric in $`zw`$.

As mentioned earlier, the connected two point function must obey the following conditions:

$$_{C_i}𝑑zR(z)T(w)_c=_{C_i}𝑑wR(z)T(w)_c=0$$
(70)
for all $`i`$, where $`C_i`$ is the contour around the $`i`$-th branch cut.

Let us first consider contour integrals around the branch cuts in $`w`$-plane. To this end we can use the first equation in (67). The first term on the right hand side is a total derivative in $`w`$ and therefore its contribution to the contour integral vanishes. Thus we arrive at the condition:

$$_{C_i}c_k(w)=0$$
(71)
Using the expression (68) for $`c_k`$, and the fact that $`w^{\mathrm{}}/(2R(w)+V^{}(w))`$ for $`\mathrm{}=0,\mathrm{},n1`$ form a complete basis of holomorphic 1-forms in the present case, these equations determine $`c_k\mathrm{}`$ completely. Very explicitly if

$`{\displaystyle \underset{k=0}{\overset{n1}{}}}t_k^{(i)}z^k`$ $`=`$ $`{\displaystyle _{C_i}}dw{\displaystyle \frac{1}{(2R(w)V^{}(w))}}[R(w){\displaystyle \frac{(V^{}(w)V^{}(z)+(zw)V^{\prime \prime }(w))}{(zw)^2}}`$ (72)
$`+{\displaystyle \frac{1}{4}}{\displaystyle \frac{(f(w)f(z)+(zw)f^{}(w))}{(zw)^2}}]`$
then using the basis of normalized differentials eq.(63),

$$\underset{k,\mathrm{}=0}{\overset{n1}{}}c_k\mathrm{}z^kw^{\mathrm{}}=\frac{1}{4}\underset{i=1}{\overset{n}{}}\underset{k=0}{\overset{n1}{}}t_k^{(i)}z^k\frac{}{S_i}f(w)$$
(73)
We will now show that $`c_k\mathrm{}`$ are symmetric in $`k`$ and $`\mathrm{}`$ exchange. Let us define a matrix $`G_{ij}`$ by the following equation

$$G_{ij}=\frac{1}{2\pi i}_{C_j}𝑑z\frac{1}{2R(z)V^{}(z)}\underset{k=0}{\overset{n1}{}}t_k^{(i)}z^k$$
(74)
We first simplify eq.(72) by using (34) as

$$\underset{k=0}{\overset{n1}{}}t_k^{(i)}z^k=_{C_i}𝑑w\frac{1}{(2R(w)V^{}(w))}\frac{2V^{}(z)V^{}(w)+f(z)+f(w)}{4(zw)^2}$$
(75)
where we have omitted a total derivative term with respect to $`w`$ since it does not contribute to the contour integral. Substituting this in eq.(74) and noting that the residue of the first order pole $`1/(zw)`$ vanishes due to eq.(34), we find that $`G_{ij}`$ is symmetric. Eq.(74) can be solved explicitly for $`t_k^{(i)}`$ as

$$\underset{k=0}{\overset{n1}{}}t_k^{(i)}z^k=\frac{1}{4}\underset{j=1}{\overset{n}{}}G_{ij}\frac{}{S_j}f(z)$$
(76)
Plugging this equation in eq.(73) and using the fact that $`G_{ij}`$ is symmetric, we find that the coefficients $`c_k\mathrm{}`$ are symmetric in $`k`$ and $`\mathrm{}`$ exchange as claimed above. .

Finally note that the symmetry of $`c_k\mathrm{}`$ implies that $`R(z)T(w)_c`$ is symmetric in $`z`$ and $`w`$.This will be crucial in the following. In particular this also implies that the contour integrals around the branch cuts in the $`z`$-plane vanish.

To summarize this subsection, although we have discussed in detail the example of $`R(z)T(w)_c`$, it is easily seen that the eq.(56) and the conditions like eq.(70) fix the solutions for all the connected two point functions uniquely.

### 4.3 Solutions for the connected two functions

Note that the right hand side of equation (56) vanishes for all components except $`(1,1),(1,2),(2,1)`$ and (3,3) (the last being proportional to $`ϵ_{\alpha \beta }`$). The uniqueness of the solution then implies

$$R(z)R(w)_c=w_\alpha (z)T(w)_c=w_\alpha (z)R(w)_c=w_{(\alpha }(z)w_{\beta )}(w)_c=0$$
(77)
This is in accordance with the estimates given in the Table 1. We have already obtained the solution for $`R(z)T(w)_c`$ in equations (69,72,73). Let us denote this solution as $`G^2H(z,w)`$.

The remaining two equations are:

$`M(z)w^\alpha (z)w_\alpha (w)_c`$ $`=`$ $`{\displaystyle \frac{10}{3}}G^2_wR(z,w)`$
$`M(z)T(z)T(w)_c`$ $`+`$ $`2T(z)R(z)T(w)_c={\displaystyle \frac{1}{3}}G^2_wT(z,w)`$ (78)
Comparing the first equation with that of $`R(z)T(w)_c`$ namely eq.(66) we conclude that

$$w^\alpha (z)w_\alpha (w)_c=10G^2H(z,w)$$
(79)
Finally taking the derivative of eq.(66) with respect to $`N_i\frac{}{S_i}`$, and using the fact that to the leading order (i.e. order(1)) $`T(z)=(N_i\frac{}{S_i}+\frac{1}{2}w_i^\alpha w_{\alpha j}\frac{^2}{S_iS_j})R(z)`$, we obtain the following equation:

$$M(z)N_i\frac{}{S_i}R(z)T(w)_c+2T(z)R(z)T(w)_c=\frac{1}{3}G^2_wT(z,w)$$
(80)
where we have used the fact that $`G^2w_{\alpha i}`$ is trivial. This equation is the same as the second equation of (78). Uniqueness of the solution then implies that upto order $`G^2`$,

$$T(z)T(w)_c=G^2N_i\frac{}{S_i}H(z,w)$$
(81)
We are now in a position to compute the gravitational corrections to the one point functions of $`R(z)`$, $`w_\alpha (z)`$ and $`T(z)`$ from eq.(33). Note that in this equation the two point functions contain both the disconnected and connected pieces. Since $`R(z)R(z)_c`$ and $`R(z)w_\alpha (z)_c`$ vanish, the first two equations do not contain any connected parts. Uniqueness then implies that one point function of $`R`$ and $`w_\alpha =w_\alpha ^i\frac{}{S_i}R`$ do not get any gravitational correction. The non-trivial equation is the third one. Using the results of this subsection we get

$$(M(z)\frac{1}{3}G^2T(z))T(z)+12G^2H(z,z)=0$$
(82)
Expanding $`T(z)=T^{(0)}+G^2T^{(1)}`$, with $`T^{(0)}=N_i\frac{}{S_i}R+\frac{1}{2}w_i^\alpha w_{\alpha j}\frac{}{S_i}\frac{}{S_j}R`$ we obtain

$$M(z)T^{(1)}(z)+[\frac{1}{3}(T^{(0)}(z))^2]+12H(z,z)=0$$
(83)
Here the term indicated in the square bracket goes as $`N^2`$ (note that only the $`N_i`$ dependent term in $`T^{(0)}`$ contribute since $`G^2w_\alpha `$ is trivial) and therefore represents genus 0 contribution. On the other hand the term proportional to $`H(z,z)`$ does not come with any factors of $`N_i`$ as is seen from the explicit solution given in (69), (72) and (73). This contribution therefore comes from genus 1. Writing $`T^{(1)}=T_0^{(1)}+T_1^{(1)}`$ where the subscript denotes the genus, we have the following solution to the above equation:

$`T_0^{(1)}(z)`$ $`=`$ $`{\displaystyle \frac{1}{6}}N_i{\displaystyle \frac{}{S_i}}T^{(0)}(z)`$
$`T_1^{(1)}(z)`$ $`=`$ $`{\displaystyle \frac{12}{(2R^{(0)}(z)V^{}(z))}}[H(z,z)+c^{(1)}(z)]`$ (84)
where $`c^{(1)}(z)`$ is a polynomial of degree $`n2`$ and is uniquely determined by the requirement that the contour integrals of $`T_1^{(1)}(z)`$ around every branch cut vanishes. In the next subsection we will show that this is exactly the answer the Matrix model provides.

Let us note that the genus 0 contribution $`T_0^{(1)}`$ above can be absorbed in $`T^{(0)}`$ by a field redefinition

$$S_iS_i+\frac{1}{6}G^2N_i$$
(85)
In particular this means that the contribution of $`T_0^{(1)}`$ to the effective superpotential can be absorbed by the above field redefinition into the original genus 0 effective superpotential in the absence of the gravitational field. This also implies that this term does not contribute to the superpotential when evaluated at the classical solution of $`S_i`$ in agreement with the statement made in .

### 4.4 Comparison with the Matrix model results

In the Matrix model a systematic approach to computing higher genus contributions has been developed in , however in the following we will rederive their results in a way parallel to the gauge theory discussion above. This will make the comparison between the two very transparent. Consider a hermitian matrix model with action given by $`S=\frac{\widehat{N}}{g_m}_k\frac{g_k}{k}\mathrm{tr}M^k\frac{\widehat{N}}{g_m}V`$, where $`M`$ is a hermitian $`\widehat{N}\times \widehat{N}`$ matrix. In Matrix model the resolvent $`\mathrm{\Omega }(z)\frac{g_m}{\widehat{N}}\mathrm{tr}\frac{1}{\mathrm{z}\mathrm{M}}`$ satisfies a loop equation similar to gauge theory $`R(z)`$.

$$\mathrm{\Omega }(z)^2I(z)\mathrm{\Omega }(z)+\mathrm{\Omega }(z)\mathrm{\Omega }(z)_c=0$$
(86)
Here $`I(z)`$ is the same integral operator as in the gauge theory discussion above. In the large $`\widehat{N}`$ limit, the two point function factorizes. However in the subleading order in $`1/\widehat{N}^2`$ the connected part of the two point function (in fact the planar connected graph) contributes and which in turn yields the genus 1 contribution to the resolvent via the above equation. By definition

$`\mathrm{\Omega }(z)\mathrm{\Omega }(w)_c`$ $`=`$ $`{\displaystyle \frac{1}{\widehat{N}^2}}{\displaystyle \underset{k=0}{\overset{\mathrm{}}{}}}{\displaystyle \frac{k}{z^{k+1}}}{\displaystyle \frac{}{g_k}}\mathrm{\Omega }(w)`$ (87)
$``$ $`{\displaystyle \frac{1}{\widehat{N}^2}}𝒪(z)\mathrm{\Omega }(w)`$
Let us expand the 1-point function as :

$$\mathrm{\Omega }(z)=\mathrm{\Omega }_{(0)}(z)+\frac{1}{\widehat{N}^2}\mathrm{\Omega }_{(1)}(z)+\mathrm{}$$
(88)
where dots represent terms of higher order in $`1/\widehat{N}^2`$. Inserting these expansions in the above equations we get:

$`\mathrm{\Omega }_{(0)}(z)^2I(z)\mathrm{\Omega }_{(0)}(z)`$ $`=`$ $`0`$
$`(2\mathrm{\Omega }_{(0)}(z)I(z))\mathrm{\Omega }_{(1)}(z)`$ $`+`$ $`𝒪(z)\mathrm{\Omega }_{(0)}(z)=0`$ (89)
Now we need to solve for $`𝒪(w)\mathrm{\Omega }_{(0)}(z)`$. This can be done by applying the differential operator $`𝒪(w)`$ on the first equation of (89). To this end we need the following identity:

$$𝒪(w)V^{}(y)=\underset{k=1}{\overset{\mathrm{}}{}}\frac{k}{w}(\frac{y}{w})^k=\frac{1}{(wy)^2}$$
(90)
which is valid for $`|w|>|y|`$. It follows that

$`{\displaystyle _{C_z}}𝑑y[𝒪(w)V^{}(y)]{\displaystyle \frac{\mathrm{\Omega }_{(0)}(y)}{yz}}`$ $`=`$ $`{\displaystyle _{C_z,|y|<|w|}}𝑑y{\displaystyle \frac{1}{(wy)^2}}{\displaystyle \frac{\mathrm{\Omega }_{(0)}(y)}{yz}}`$ (91)
$`=`$ $`_w{\displaystyle \frac{\mathrm{\Omega }_{(0)}(z)\mathrm{\Omega }_{(0)}(w)}{zw}}.`$
Using this, we obtain the following equation by applying $`𝒪(w)`$ on the first equation of (89)

$$(2\mathrm{\Omega }_{(0)}(z)I(z))𝒪(w)\mathrm{\Omega }_{(0)}(z)_w\frac{\mathrm{\Omega }_{(0)}(z)\mathrm{\Omega }_{(0)}(w)}{zw}.$$
(92)
Since $`\mathrm{\Omega }_{(0)}`$ of the matrix model is the same as the $`R^{(0)}`$ for the gauge theory, we see that $`𝒪(w)\mathrm{\Omega }_{(0)}(z)`$ satisfies the same equation (66) as $`3R(z)T(w)_c`$. We now impose the conditions

$$_{C_i}𝑑z𝒪(w)\mathrm{\Omega }_{(0)}(z)=_{C_i}𝑑w𝒪(w)\mathrm{\Omega }_{(0)}(w)=0$$
(93)
which are the analogues of the equations (70). It follows from the discussion of uniqueness that $`𝒪(w)\mathrm{\Omega }_{(0)}(z)`$ is equal to $`3H(z,w)`$. Note that as we showed in the last subsection, $`H(z,w)`$ is symmetric in $`z`$ and $`w`$. This is consistent with the fact that $`𝒪(w)\mathrm{\Omega }_{(0)}(z)`$ is symmetric in $`z`$ and $`w`$. Finally, substitution of $`𝒪(z)\mathrm{\Omega }_{(0)}(z)`$ in the second equation of (89), results in an equation for $`\mathrm{\Omega }_{(1)}(z)`$ which is identical to that for the genus 1 part of the gauge theory equation (83) for $`T_1^{(1)}`$. Using the fact that the integral of $`\mathrm{\Omega }_{(1)}(z)`$ around every branch cut is zero, we conclude, from the uniqueness of the solution, that

$$\mathrm{\Omega }_{(1)}(z)=\frac{1}{4}T_1^{(1)}(z),$$
(94)
with the right hand side being the genus 1 part of the solution given in (84). While in the matrix model the $`\frac{1}{\widehat{N}^2}`$ correction to the effective potential is obtained by integrating the asymptotic expansion of $`\mathrm{\Omega }_{(1)}(z)`$ with respect to the couplings $`g_k`$, the order $`G^2`$ correction to the effective superpotential in gauge theory is obtained by integrating the asymptotic expansion of $`T_1^{(1)}(z)`$ with respect to the coupling constants $`g_k`$ (we already argued in the last subsection that the genus 0 contribution coming from $`T_0^{(1)}`$ can be absorbed by a field redefinition of $`S_i`$). Eq.(94) implies therefore that the genus 1 contribution to the effective potential in matrix model is equal to the genus 1 contribution to the order $`G^2`$ term in the gauge theory effective superpotential. In fact, the relative coefficient $`4`$ in eq.(94) is exactly reproduced if one follows the numerical factors in the diagrammatic computations given in .

## 5 Conclusions

In this paper we have analyzed the problem of computing the first non-trivial gravitational corrections to the effective superpotential, $`\mathrm{\Gamma }_1=d^4xd^2\theta G^2F_1(S)`$, resulting from integrating out an adjoint scalar superfield with polynomial tree level superpotential.and minimal coupling to $`N=1`$ supergravity. Whereas the problem has been earlier analyzed by  who showed, by using diagrammatic techniques, that the effective superpotential is in fact given by the genus one correction to free energy in the corresponding matrix model, we have here considered the issue from the point of view of “loop equations” arising from the generalized Konishi anomaly in presence of an $`N=1`$ gauge $`W_\alpha `$ and supergravity $`G_{\alpha \beta \gamma }`$ background superfields, whose lowest components are the gaugino and the gravitino field strength respectively. It should be possible to generalize these methods to the case of other classical groups, as well as different representations for the chiral fields.

From this point of view, the appearance of non-planar (genus one) diagrams is related to the lack of the usual factorization property of gauge invariant correlators in the presence of a non-trivial supergravity background, The fact that a non-trivial supergravity background deforms the chiral ring was stressed and used explicitly in the work by , where the failure of the factorization property was implicit. In this paper we have exploited both these facts and included also the supergravity contribution to the Konishi anomaly. First, we have estimated various two-and three- point connected correlators using the modified chiral ring relations, then we have derived loop equations which, due to the absence of factorization, involve connected two-point functions of “resolvents” $`R(z)`$, $`T(z)`$ and $`w_\alpha (z)`$, generalizing the equations obtained by  in two ways: first because we found additional contributions of order $`G^2`$, and second because we found it useful to consider two-point functions involving operators at two different points in the complex plane. We have then shown that these loop equations satisfy quite non-trivial consistency conditions and that there is a unique solution for the non-vanishing connected two-point functions. Using these results, we have then solved for the $`𝒪(G^2)`$, genus one, correction to the superpotential and found agreement with the matrix model result. We also find at $`𝒪(G^2)`$, a genus zero contribution to the superpotential, which can be removed by a field redefinition.

Concerning the issue of factorization (i.e. position independence) of chiral correlators in a supergravity background, we have observed in section 3 that our results do not contradict the known fact that this property holds in $`N=1`$ gauge theories on Kähler manifolds (with non-trivial twisting if the manifold is not Hyper-Kahler) , since we indeed find that for purely bosonic background connected correlators vanish.

Finally, it is a non-trivial open problem to extend the above results to higher genera. Indeed, whereas the chiral ring deformation that we have here employed follows directly from standard $`N=1`$ supergravity tensor calculus. it has been shown by , using diagrammatic techniques, that in order to capture these higher order corrections one has to modify the chiral ring relations rather drastically, i.e. one has to include a self-dual two form $`F_{\alpha \beta }`$, which is a remnant of the graviphoton field strength of the parent $`N=2`$ supergravity. The modified relation, $`\{W_\alpha ,W_\beta \}=2G_{\alpha \beta \gamma }W^\gamma +F_{\alpha \beta }`$,. does not have a conventional interpretation in $`N=1`$ supergravity. It would be interesting to see how this non-standard modification can be incorporated in the approach followed here, based on generalized Konishi anomaly relations.

###### Acknowledgments.

It is a pleasure to thank George Thompson for numerous valuable discussions and for collaboration in the early phase of this project. The work of the authors is supported in part by EEC contract EC HPRN-CT-2000-00148.

## Appendix A Conventions

In this paper we follow conventions of . Raising and lowering of spinor indices are done by the $`ϵ`$ tensor as follows

$`W^\alpha =ϵ^{\alpha \beta }W_\beta ,W_\alpha =ϵ_{\alpha \beta }W^\beta ,`$ (95)
$`ϵ^{\alpha \beta }ϵ_{\beta \alpha }=2,ϵ^{\alpha \beta }ϵ_{\beta \alpha ^{}}=\delta _\alpha ^{}^\alpha .`$
We also define the following products of the $`N=1`$ gauge multiplet and the $`N=1`$ Weyl multiplet

$`W^2`$ $`=`$ $`W^\alpha W_\alpha =W_\alpha W^\alpha ,`$ (96)
$`G^2=G^{\alpha \beta \gamma }G_{\alpha \beta \gamma },`$ $`G_{\delta \gamma \alpha }G_\beta ^{\delta \gamma }=G_{\delta \gamma \beta }G_\alpha ^{\delta \gamma }={\displaystyle \frac{ϵ_{\alpha \beta }}{2}}G^2.`$
## Appendix B The gravitational contribution to Konishi anomaly

In this appendix we fix the normalization constant of the pure gravitational contribution to the the Konishi anomaly which was used in (32). The Konishi anomaly equation including the pure gravitational contribution in superspace is given by

$$\overline{D}^2(\overline{\mathrm{\Phi }}e^V\mathrm{\Phi })=\frac{1}{32\pi ^2}\mathrm{Tr}_{\mathrm{Ad}}(\mathrm{W}^2)+\alpha \frac{1}{32\pi ^2}\mathrm{G}^2$$
(97)
Here $`\alpha `$ is the unknown normalization constant which we will fix below. The $`\theta ^2`$ component of the above equation together with its anti-holomophic counterpart should reduce to the familiar equation of the chiral anomaly including the gravitational contribution given below.

$$_\mu (\overline{\psi }\overline{\sigma }^\mu \psi )_\mu (\psi \sigma ^\mu \overline{\psi })=\frac{1}{32\pi ^2}\frac{i}{2}ϵ^{mnlk}\mathrm{Tr}_{\mathrm{Ad}}(F_{mn}F_{lk})+\frac{1}{32\pi ^2}\frac{1}{24}\frac{i}{2}ϵ^{mnlk}R_{mnst}R_{lk}^{st}$$
(98)
The coefficients in the above equation have been obtained from  Our strategy to fix the normalization constant $`\alpha `$ will be to extract the contribution of $`RR`$ from the superspace equation in (97) and require it to agree with the coefficient in (98). From  the lowest component of the $`N=1`$ Weyl multiplet starts off as the gravitino field strength and is given by

$$G_{\alpha \beta \gamma }=\frac{1}{12}(\sigma _{\alpha \beta }^{ab}\psi _{ab\gamma }+\sigma _{\beta \gamma }^{ab}\psi _{ab\alpha }+\sigma _{\gamma \alpha }^{ab}\psi _{ab\beta }),$$
(99)
where we have set the auxiliary field in the above formula to zero, as we are working on shell and $`a,b`$ refer to the local Lorentz indices. The gravitino field strength is defined by

$$\psi _{ab}^\alpha =\widehat{D}_a\psi _b^\alpha \widehat{D}_b\psi _a^\alpha ,\widehat{D}_a\psi _b^\alpha =_a\psi _b^\alpha +\psi _b^\beta \omega _{n\beta }^\alpha ,$$
(100)
where $`\omega _{n\beta }^\alpha `$ is the spin connection. As we are working with an on shell background, equations of motion for the gravitino imply $`\sigma _{\alpha \beta }^{ab}\psi _{ab\gamma }=\sigma _{\alpha \gamma }^{ab}\psi _{ab\beta }`$. Therefore on shell, the lowest component of the Weyl multiplet is given by

$$G_{\alpha \beta \gamma }=\frac{1}{4}\sigma _{\alpha \beta }^{ab}\psi _{ab\gamma }$$
(101)
The supersymmetric transformation on the gravitino field strength is given by

$$\delta \psi _{mn}^\alpha =\xi ^\beta R_{mnab}\sigma _\beta ^{ab\alpha }+\mathrm{}$$
(102)
The dots in the above equation all refer to terms that involve the fermions, which are not of interest for the present purpose of determining the coefficient of $`RR`$. Substituting this variation in (101) we obtain

$$\delta G_{\alpha \beta \gamma }=\frac{1}{4}\sigma _{\alpha \beta }^{ab}\sigma _{\delta \gamma }^{cd}\xi ^\delta R_{abcd}$$
(103)
The $`\theta ^2`$ component of $`G^2`$ contains the $`RR`$ term, which is imaginary and thus contributes to the anomaly. This is given by

$$G^2|_{\theta ^2}=\frac{i}{2}\frac{1}{16}ϵ^{a^{}b^{}ab}R_{a^{}b^{}cd}R_{ab}^{cd}$$
(104)
The total contribution to the anomalous current is obtained by subtracting this out with the $`G_{\dot{\alpha }\dot{\beta }\dot{\gamma }}`$, the anti-holomorphic contributions. Therefore the coefficient of the $`RR`$ term from $`W^2`$ is $`1/8`$. Comparing with the coefficient of $`RR`$ in (98) we see that $`\alpha =1/3`$ in (97) in order to reproduce the chiral anomaly, including the gravitational contribution.