Lin Fei Thesis
Lin Fei Thesis
Lin Fei Thesis
Lin Fei
A Dissertation
Presented to the Faculty
of Princeton University
in Candidacy for the Degree
of Doctor of Philosophy
September 2017
c Copyright by Lin Fei, 2017.
scalar O(N ) theory in d > 4. We study this theory to three-loop order and show
that various operator dimensions are consistent with large- N results. This theory
possesses an IR stable xed point at real couplings for N > 1038, suggesting the
CFT in d = 5. For the special case of Sp(2), the IR xed point possesses an enhanced
symmetry given by the supergroup OSp(1|2). We also observe that various operator
dimensions of the Sp(2) theory match those from the 0-state Potts model. We provide
a graph theoretic proof showing that the zero, two, and three-point functions in the
Sp(2) model and the 0-state Potts model indeed match to all orders in perturbation
theory, strongly suggesting their equivalence. We then study two fermionic theories
tain operator dimensions, passing all checks against large- N results. We use two
sided Padé approximations with our -expansion results to obtain estimates of vari-
ous quantities in the physical dimension d = 3. Finally, we provide evidence that the
iii
Acknowledgements
I would like to express my sincerest gratitude to my advisor Simone Giombi, and Igor
Klebanov for their guidance, insights, and encouragement throughout the years. I
would also like to thank my colleague Grigory Tarnopolsky for all the help he has
low graduate students, especially Mykola Dedushenko, Kenan Diab, Po-Shen Hsin,
family.
research.
The National Science Foundation for providing funding for my work as well as
iv
To Helen and Neil
v
Contents
Abstract . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . iii
Acknowledgements . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . iv
1 Introduction 1
vi
3 Three-Loop Analysis of the Critical O(N ) Models in 6− Dimensions 53
3.6.1 Counterterms . . . . . . . . . . . . . . . . . . . . . . . . . . . 79
5 Equivalence between the Sp(2) model and the 0-state Potts model 95
5.1 Preliminaries . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
vii
6 Yukawa CFTs and Emergent Supersymmetry 111
ries 169
Bibliography 176
viii
Chapter 1
Introduction
string theory to condensed matter physics. They are Quantum Field Theories (QFT)
the critical points of statistical systems, the correlation length scale ξ diverges. In the
continuum theory that describes the physics across distances that are large compared
to the lattice spacing, the correlation functions are insensitive to the microscopic de-
tails of the model and become scale-invariant. Thus, CFTs play a primary role in the
of seemingly unrelated statistical systems share some universal large scale proper-
ties despite their microscopic dierences. In order to understand the origin of these
are summing over a large number of random variables with arbitrary distribution
(provided that large deviations with respect to the mean decrease fast enough), the
1
asymptotic behavior is always a Gaussian, regardless of the distribution of the orig-
inal random variables. One can understand this using a renormalization strategy by
averaging over two identical random variables with distribution ρ results in a random
variable with a transformed distribution T ρ. However, the limiting distribution ρ∗
is necessarily a xed point of this operation T ρ∗ = ρ∗ . One nds that the Gaussian
the initial distribution is, after repeated application of the averaging operation, all of
described by the Hamiltonian H, and a microscopic scale 1/Λ describing, for exam-
ple, the lattice spacing for a lattice model. Now we recursively integrate out short
d
λ Hλ = T [Hλ ] (1.1)
dλ
T [H∗ ] = 0 (1.2)
If the above equations have an attractive xed point, then all the Hamiltonians
within the basin of attraction will eventually ow to the same xed point - they
belong to the same universality class. The Gaussian xed point is again a stable xed
point in QFTs, but they describe free theories without interactions. However, there
2
The linear partial dierential equations for the correlators of the renormalized
using these RG equations, one can compute two functions: the β -function which
describe the rate at which the couplings in H ows, and the anomalous dimension η,
which describes the correction to the classical scaling dimensions of various operators
when interactions are turned on. All of these quantities are computed order by order in
perturbation theory, with higher order corrections coming from higher loop Feynman
diagrams.
Ising model. Near the critical point, it is described by a Euclidean QFT of a real
Z
3 1 λ
S= dx (∂φ)2 + φ4 . (1.3)
2 4!
The coupling constant λ is the strength of the interaction, which will change with scale
due to renormalization. It turns out, this model is weakly coupled at short distances,
but it becomes strongly coupled at long distances. It is the long distance regime that
is needed for describing the critical behavior, but perturbation theory in λ cannot be
various operators. However, theorists have invented ingenious expansion schemes that
d = 4− where 1. Then the IR stable xed point of the Renormalization Group
3
may be developed [1]. The coecients of the rst few terms fall o rapidly, so that
setting =1 provides a rather precise approximation that is in good agreement with
Another important idea has been the large- N expansion in the O(N ) symmetric
λ
QFT of N real scalar elds φi , i = 1, . . . , N , with interaction
4
(φi φi )2 . Using a gen-
eralized Hubbard-Stratonovich transformation with an auxiliary eld σ, it is possible
to develop expansion in powers of 1/N (for a comprehensive review, see [3]). The
large N expansion may be developed for a range of d [4, 5, 6, 7, 8, 9, 10, 11, 12] and
compared with the regimes where other perturbative expansions such as dimensional
continuation.
any N [2], and the results can be matched with the large N techniques. Also, for
d = 2 + the large N results match with the perturbative UV xed point of the O(N )
Non-linear Sigma Model (NL σ M). Using the rst few terms in the large N expansion
directly in d=3 provides another approach to estimating the scaling dimensions for
low values of N. Thus, a combination of the large N and expansions provides good
approximations for the critical behavior in the entire range 2 < d < 4. One should
note, however, that both expansions are not convergent but rather provide asymptotic
series. There are continued eorts towards obtaining a more rigorous approach to the
O(N ) symmetric CFT's using conformal bootstrap ideas [13, 14, 15, 16], and recently
it has led to more precise numerical calculations of the operator scaling dimensions
in three-dimensional CFT's [17, 18]. The bootstrap approach has been successfully
applied not only in the range 2<d<4 [19], but also for 4<d<6 [20].
4
1.2 Review of Wilson-Fisher -expansion
Let us consider the Euclidean eld theory of N real massless scalar elds with an
Z
d 1 λ
S= d x (∂φi )2 + (φi φi )2 . (1.4)
2 4
As follows from dimensional analysis, the interaction term is relevant for d<4 and
irrelevant for d > 4. Hence, for 2 < d < 4, it is expected that the quartic interaction
generates a ow from the free UV xed point to an interacting IR xed point. In
d = 4− this xed point can be studied perturbatively in the framework of the
Wilson-Fisher -expansion [1, 2]. Indeed, the one-loop beta function for the theory in
d=4− reads
λ2
βλ = −λ + (N + 8) , (1.5)
8π 2
8π 2
λ∗ = . (1.6)
N +8
Higher order corrections in will change the value of the critical coupling, but not
its existence, at least in perturbation theory. The anomalous dimensions of the fun-
damental eld φi and the composite φi φi at the xed point can be computed to be,
to leading order in
N +2 2 N +2
γφ = 2
+ O(3 ) , γφ2 = + O(2 ) , (1.7)
4(N + 8) N +8
5
corresponding to the scaling dimensions
d N +2 2
∆φ = − 1 + γφ = 1 − + 2
+ O(3 ) , (1.8)
2 2 4(N + 8)
6
∆φ2 = d − 2 + γφ2 = 2 − + O(2 ) . (1.9)
N +8
large N , a result that follows from the large N analysis reviewed below. Higher order
corrections in for general N may be derived by higher loop calculations in the theory
For d > 4, the interaction is irrelevant and so the IR xed point is the free theory;
however, one may ask about the existence of interacting UV xed points. Working
in d = 4+ for small , one indeed nds a perturbative UV xed point at (see, for
example [23])
8π 2
λ∗ = − . (1.10)
N +8
The anomalous dimensions at this critical point are given by the same expressions
(1.7) with → −. Note that, because γφ starts at order 2 , the dimension of φ stays
above the unitarity bound for all N, at least for suciently small . However, since
the xed point requires a negative coupling, one may worry about its stability, and it
dimension d is the large N expansion. The standard technique to study the theory
Z
d 1 1 σ2
S= d x (∂φi )2 + σφi φi − . (1.11)
2 2 4λ
6
Integrating out σ via its equation of motion σ = λφi φi , one gets back to the original
example of the double trace deformations studied in [24]. One can then show that at
1
be dropped , and the eld σ plays the role of the composite operator φi φi . One may
Z
d 1 1
Scrit = d x (∂φi )2 + √ σφi φi (1.12)
2 2 N
√
where we have rescaled σ by a factor of N for reasons that will become clear
momentarily. The 1/N perturbation theory can be developed by integrating out the
Z (1.13)
1
dd xdd y σ(x)σ(y) hφi φi (x)φj φj (y)i0 +O(σ 3 )
R
= Dσ e 8N
Z
i i j j 2 dd p eip(x−y)
hφ φ (x)φ φ (y)i0 = 2N [G(x − y)] , G(x − y) = . (1.14)
(2π)d p2
Z
2 dd p ip(x−y)
[G(x − y)] = e G̃(p)
(2π)d
Z
dd q 1 (p2 )d/2−2
G̃(p) = = − d−3 (1.15)
(2π)d q 2 (p − q)2 2d (4π) 2 Γ( d−1 ) sin( πd ) 2 2
1 This applies formally to both IR and UV xed points.
7
and so from (1.13) one nds the two-point function of σ in momentum space
d+1 d−3 d−1 πd d d
hσ(p)σ(−p)i = 2 (4π) 2 Γ sin( ) (p2 )2− 2 ≡ C̃σ (p2 )2− 2 . (1.16)
2 2
d−1
2d+2 Γ 2
sin( πd ) 1 Cσ
hσ(x)σ(y)i = 3 2 ≡ . (1.17)
π 2 Γ d2 −2 |x − y|4 |x − y|4
∆ = 2. Note also that the coecient Cσ is positive in the range 2 < d < 6.
The large N perturbation theory can then be developed using the propagator
(1.16)-(1.17) for σ, the canonical propagator for φ and the interaction term σφi φi in
(1.12). For instance, the 1/N term in the anomalous dimension of φi can be computed
from the one-loop correction to the φ-propagator
Z
1 dd q 1 C̃σ
, (1.18)
N (2π) (p − q) (q 2 ) d2 −2+δ
d 2
2
regulator . Doing the momentum integral by using (A.1), one obtains the result
C̃σ (d − 4)
(p2 )1−δ Γ(δ − 1) (1.19)
N (4π) d2 dΓ( d )
2
where we have set δ = 0 in the irrelevant factors. The 1/δ pole corresponds to a
d 1 1
∆φ = − 1 + η1 + 2 η2 + . . . (1.20)
2 N N
2 One may perform the calculation using a momentum cuto, yielding the same nal result.
8
η1 Cσ
0.3
20
0.2
15
0.1
10
0.0 d
2 3 4 5 6 5
- 0.1 0 d
2 3 4 5 6
Figure 1.1: The 1/N anomalous dimension of φi and the coecient of the two-point
d−1
πd
C̃σ (d − 4) 2d−3 (d − 4)Γ 2
sin 2
η1 = d = 3 d
. (1.21)
(4π) 2 dΓ( d2 ) π2Γ 2
+1
Setting d = 4− and expanding for small , it is straightforward to check that this
agrees with the 1/N term of (1.8). The leading anomalous dimension of σ also takes
1 4(d − 1)(d − 2) 1
∆σ = 2 + η1 + O( 2 ) . (1.22)
N d−4 N
while the focus of most of the exisiting literature is on the range 2 < d < 4, we see
no obvious problems with unitarity in continuing the large N critical O(N ) theory to
the range 4 < d < 6.3 A plot of η1 and of the σ two-point function coecient Cσ is
given in Figure 1.1, showing that they are both positive for 2 < d < 6.
3 Recall that the unitarity bound for a scalar operator is ∆ ≥ d/2 − 1. One can also check that
the order 1/N term in the anomalous dimension of the O(N ) invariant higher spin currents, given
in [10], is positive in the range 2 < d < 6, consistently with the unitarity bound ∆s ≥ s + d − 2 for
spin s operators (s > 1/2).
9
1.4 The Gross-Neveu Model and its UV completion
In order to motivate our proposed cubic model as the UV completion of φ4 theory, we
rst review the relation between the Gross-Neveu (GN) model and its UV completion,
The same and 1/N -expansion schemes can be applied to fermionic theories as
g
LGN = ψ̄j 6 ∂ψ j + (ψ̄j ψ j )2 . (1.23)
2
N −2 2
β = g − g + ... , (1.24)
2π
1+ N −1 2
∆ψ = + g . (1.25)
2 8π 2
2π
Therefore this theory possesses a non-trivial UV xed point given by g∗ = N −2
in
d = 2 + . And the anomalous dimension at that xed point is:
1 1 N −1
∆ψ = + + + ... . (1.26)
2 2 4(N − 2)2
sions greater than two, and therefore in the UV, the large momentum behavior of the
to describe this xed point through some other means, such as the IR xed point of
10
some other theory. Indeed such a theory exists, and it is the Gross-Neveu-Yukawa
model, although it requires a bit more care. Introducing the auxiliary eld σ, and
dropping the quadratic term in the critical limit, we have the action
Z
1
Scrit ferm = d / i
d x − ψ̄0i ∂ ψ0 + √ i
σ0 ψ̄0i ψ0 , (1.27)
N
ip/
hψ0i (p)ψ̄0j (−p)i0 = δji . (1.28)
p2
The σ eective propagator obtained after integrating over the fundamental elds ψ0i
reads
d
hσ0 (p)σ0 (−p)i0 = C̃σ0 /(p2 ) 2 −1+∆ , (1.29)
where
d−3
d − 1 πd
C̃σ0 ≡ −2d+1 (4π) 2 Γ sin (1.30)
2 2
and we have introduced the regulator ∆. Note that the power of p2 in the propagator
d d
is
2
−1+∆ instead of
2
−2+∆ found in the scalar case. In order to cancel the
1/2
ψ = Zψ ψ0 , σ = Zσ1/2 σ0 , (1.31)
11
where
1 Zψ1 1 Zσ1
Zψ = 1 + + O(1/N 2 ), Zσ = 1 + + O(1/N 2 ) . (1.32)
N ∆ N ∆
ip/ C̃σ
hψ i (p)ψ̄j (−p)i = δji C̃ψ , hσ(p)σ(−p)i = , (1.33)
d
−∆ψ + 12 d
(p2 ) 2 (p2 ) 2 −∆σ
− + η GN ,
d 1
∆ψ = (1.34)
2 2
We can expand this expression in d=2+ and see agreement with our -expansion
result.
GNY model [26, 27] contains a scalar eld with cubic interaction with the fermions
1 1
LGNY = (∂µ σ)2 + ψ̄j 6 ∂ψ j + g1 σ ψ̄j ψ j + g2 σ 4 . (1.36)
2 24
12
In d = 4 − , the β -functions to one-loop are given by[3]:
1 2 2 4
βg 2 = −g2 + 3g2 + 2N g1 g2 − 12N g1 ,
(4π)2
N +6 3
βg 1 = − g1 + g , (1.37)
2 2(4π)2 1
where N = Nf tr1 = 4Nf . The anomalous dimension for the ψ eld is given by:
3− 1
∆ψ = + 2
g12 . (1.38)
2 2(4π)
The model possesses an IR stable xed point at the critical couplings gi∗ given by
(g1∗ )2 1
2
= ,
(4π) N +6
√
g2∗ −N + 6 + N 2 + 132N + 36
= . (1.39)
(4π)2 6(N + 6)
3 N +5
∆ψ = − . (1.40)
2 2(N + 6)
Which again matches against the large N result expanded near d = 4−. Therefore
the UV xed point of the GN model is described by the IR xed point of the GNY
Similarly, in this thesis, we are trying to show that the O(N ) symmetric ((φi )2 )2
theory with a UV xed point in d = 4+ has a UV completion given by a cubic
1 1 g1 g2
L = (∂µ φi )2 + (∂µ σ)2 + σφi φi + σ 3 . (1.41)
2 2 2 6
13
We will demonstrate this by comparing against large N results of various operator
dimensions.
work published in [28, 29, 30, 31] respectively, with co-authors Simone Giombi, Igor
in 6− dimensions shows that this theory has an IR stable xed point at real values
of the coupling constants for N > 1038. We show that the 1/N expansions of various
operator scaling dimensions match the known results for the critical O(N ) theory
continued to d = 6 − . These results suggest that, for suciently large N, there are
5-dimensional unitary O(N ) symmetric interacting CFT's; they should be dual to the
Vasiliev higher-spin theory in AdS6 with alternate boundary conditions for the bulk
scalar. Using these CFT's we provide a new test of the 5-dimensional F -theorem,
and also nd a new counterexample for the CT theorem.
In chapter 3, we continue the study of the same theory to three loop order. We
calculate the 3-loop beta functions for the two couplings and use them to determine
also use the beta functions to determine the corrections to the critical value of N
below which there is no xed point at real couplings. We also study the theory with
IR stable xed point at imaginary couplings which can be reached by ow from a
14
scaling dimensions at the IR xed point of the N =1 theory and suggest that, upon
Borrowing our three loop calculations in 3, we develop expansions for several op-
erator dimensions and for the sphere free energy F. The conjectured F -theorem is
obeyed in spite of the non-unitarity of the theory. Our results point to the exis-
operator dimensions are real. For N =2 we show that the IR xed point possesses
an enhanced global symmetry given by the supergroup OSp(1|2). This suggests the
the 6 − expansions of the scaling dimensions and sphere free energy in our OSp(1|2)
model are the same as in the q→0 limit of the q -state Potts model.
order, the anomalous dimensions and beta functions of these two theories are equal
up to a sign alternating with successive loop order. Using a graph theoretic approach,
we prove this equality as well as the alternating sign. This argument holds for all
loop orders in perturbation theory. This strongly suggests that these two models are
the sphere free energy and certain operator scaling dimensions using dimensional
continuation. We provide new evidence that the 4− expansion of the N =1 Gross-
15
appear to be in good agreement with the available results obtained using the numer-
ical conformal bootstrap. We apply a similar approach to calculate the sphere free
one 2-component Dirac fermion, this theory has an emergent supersymmetry with 4
16
Chapter 2
Giombi, and Igor Klebanov. We also thank David Gross, Daniel Harlow, Igor
Herbut, Juan Maldacena, Giorgio Parisi, Silviu Pufu, Leonardo Rastelli, Dam Son,
We would like to extend the results of the -expansion and large- N expansion to
interacting O(N ) models in d > 4. At rst glance, such extensions seem impossible:
the (φi φi )2 interaction is irrelevant at the Gaussian xed point, since it has scaling
the free eld theory. However, at least for large N, the theory possesses a UV stable
transformation. At the interacting xed point, the scaling dimension of the operator
φi φi is 2 + O(1/N ) for any d. For d > 6, this dimension is below the unitarity bound
d/2 − 1. It follows that the interesting range, where the UV xed point may be
17
unitary, is [32, 33, 34, 35]
4<d<6. (2.1)
We will examine the structure of the O(N ) symmetric scalar eld theory in this range
from various points of view.
We rst note that the coecients in the 1/N expansions derived for various quan-
tities in [4, 5, 6, 7, 8, 9, 10, 11, 12] may be continued to the range (2.1) without any
obvious diculty. Thus, the UV xed points make sense at least in the 1/N expan-
sion. However, one should be concerned about the stability of the xed points in this
2
coupled, it occurs for the negative quartic coupling λ∗ = − N8π+8 + O(2 ) [36, 23].
Thus, it seems that the UV xed point theory is unlikely to be completely stable,
point theory, it would be helpful to describe it via RG ow from another theory. In
such a UV complete" description, the O(N ) symmetric theory we are after should
1 1 g1 g2
L = (∂µ φi )2 + (∂µ σ)2 + σφi φi + σ 3 . (2.2)
2 2 2 6
The cubic interaction terms are relevant for d < 6, so that the theory ows from the
weakly coupled for d = 6 − . The idea to study a cubic scalar theory in d=6− is
not new. Michael Fisher has explored such an -expansion in the theory of a single
scalar eld as a possible description of the Yang-Lee edge singularity in the Ising
model [37]. In that case, which corresponds to the N = 0 version of (2.2), the IR
18
xed point is at an imaginary value of g2 . 1 This is related to the lack of unitarity of
the xed point theory. Using the one-loop beta functions for g1 and g2 , we will show
in Section 2.3 that for large N the IR stable xed point instead occurs for real values
of the couplings, thus removing conict with unitarity. Remarkably, such unitary
xed points exist only for N > 1038. As we show in Section 2.4, for N ≤ 1038 the
IR stable xed point becomes complex, and the theory is no longer unitary. Thus,
the breakdown of the large N expansion in d = 6− occurs at a very large value,
Ncrit = 1038. We will provide evidence, however, that for the physically interesting
Besides its intrinsic interest, the O(N ) invariant scalar CFT in d=5 has inter-
esting applications to higher spin AdS/CFT dualities. There exists a class of Vasiliev
theories in AdSd+1 [41, 42, 43, 44, 45, 46] that is naturally conjectured to be dual to
the O(N ) singlet sector of the d-dimensional CFT of N free scalars [47, 23]. In order to
λ
extend the duality to interacting CFT's, one adds the O(N ) invariant term
4
(φi φi )2 .
For d=3 this leads to the well-known Wilson-Fisher xed points [1, 2]. In the dual
Vasiliev theory in AdS6 [23]. One can adopt the ∆+ = 3 boundary conditions on the
bulk scalar, which are necessary for the duality to the free O(N ) theory. Alterna-
tively, the ∆− = 2 + O(1/N ) are allowed as well [48, 49]. This suggests that the dual
interacting O(N ) CFT should exist in d = 5, at least for large N [36, 23]. Our RG
calculations lend further support to the existence of this interacting d=5 CFT.
In Sections 2.3 and 2.5, using one-loop calculations for the theory (2.2) in d = 6−,
we nd some IR operator dimensions to order , while keeping track of the dependence
on 1/N to any desired order. We will then match our results with the 1/N expansions
1 Renormalization group calculations for similar cubic theories in d = 6− were carried out in
[38, 39, 40].
19
3dFCriticalFO(N)FCFT 5dFCriticalFO(N)FCFT
2 4 6
Figure 2.1: Interacting unitary O(N ) symmetric scalar CFT's exist for dimensions
2 < d < 6, with d = 4 excluded. In 6− and 4− dimensions they may be described as
weakly coupled IR xed points of the cubic and quartic scalar theories, respectively.
In 4+ and 2+ dimensions they are weakly coupled UV xed points of the quartic
derived for the (φi φi )2 theory in [4, 5, 6, 7, 8, 9, 10, 11, 12], evaluating them in d = 6−.
The perfect match of the coecients in these two 1/N expansions provides convincing
evidence that the IR xed point of the cubic O(N ) theory (2.2) indeed describes the
same physics as the UV xed point of the (φi φi )2 theory. Our results thus provide
evidence that, at least for large N, the interacting unitary O(N ) symmetric scalar
CFT's exist not only for 2 < d < 4, but also for 4 < d < 6. In Fig. 1 we sketch the
entire available range 2 < d < 6, pointing out the various perturbative descriptions
In Section 2.6 we discuss the large N results for CT , the coecient of the two-
for our proposal that the IR xed point of (2.2) describes the O(N ) symmetric CFT.
We show that the RG ow from this interacting CFT to the free theory of N scalars
provides a counter example to the conjectured CT theorem. On the other hand, the
20
Our discussion of the (φi φi )2 scalar theory in the range 4<d<6 is analogous to
the much earlier results [26, 27, 3] about the Gross-Neveu model [25] in the range 2<
d < 4. The latter is a U (Ñ ) invariant theory of Ñ Dirac fermions with an irrelevant
not hard to show that this model has a UV xed point, at least for large Ñ [3].
This CFT was conjectured [51, 52] to be dual to the type B Vasiliev theory in AdS4
this CFT is via the Gross-Neveu-Yukawa (GNY) model [26, 27, 3], which is a theory
hard to show that the one-loop beta functions have an IR stable xed point for any
positive Ñ . The operator dimensions at this xed point match the large Ñ treatment
of the Gross-Neveu model [53]. These results will be reviewed in Section 2.7, where
we also discuss tests of the 3-d F-theorem [54, 50] provided by the GNY model.
compute the anomalous dimension of φi to order 1/N 3 and that of σ to order 1/N 2 .
These results have been successfully matched to all available orders in the d=2+
and d = 4− expansions, providing a strong test of their correctness. The explicit
21
form of η2 in general dimensions, dened as in (1.20), reads [4, 5]:
2
η2 = 2η12 (f1 + f2 + f3 ) ;
µ2 + µ − 1 µ 0 µ(3 − µ)
f1 = v 0 (µ) + , f2 = v (µ) + ,
2µ(µ − 1) 2−µ 2(2 − µ)2
µ(2µ − 3) 0 2µ(µ − 1)
f3 = v (µ) + ;
2−µ 2−µ
v 0 (µ) = ψ(2 − µ) + ψ(2µ − 2) − ψ(µ − 2) − ψ(2) ,
Γ0 (x) d
ψ(x) = , µ= . (2.3)
Γ(x) 2
The expressions for the dimension of σ at order 1/N 2 [5] and for the coecient η3
[6] in arbitrary dimensions are lengthy and we do not report them explicitly here.
3
However, since they are useful to test our general picture, we write below the explicit
1 44 1936 11 835 16352
∆φ = 2 − + + 2+ + ... − + + + . . . 2 + O(3 )(2.4)
2 N N N3 12N 6N 2 N3
and
40 6800 104 34190
∆σ = 2 + + + ... − + + . . . 2 + O(3 ) . (2.5)
N N2 3N 3N 2
In the next section, we will show that the order terms precisely match the one-loop
2 Ref. [5] contains an apparent typo: in the denition of the function v 0 (µ) given in Eq. (21),
α=µ−1 should be replaced by α = µ − 2. The correct formula for η2 may be found in the earlier
paper [4].
3 Note that [5] derives a result for the critical exponentν , which is related to the dimension of σ
1
by ∆σ = d − ν . We also note that a misprint in eq. (22) of [6] has been later corrected in eq. (11)
of [53].
22
Using the results in [4, 5, 6, 53] we may also compute the dimension of φi and
non-trivial self-energy integral that was not evaluated for general dimension in [6].
An explicit derivation of this integral in general d was later obtained in [55]. Using
3 32 1427456
∆φ = + 2
−
2 15π N
3375π 4 N 2
275255197696 89735168 32768 ln 4 229376ζ(3) 1
+ 6
− 4
+ 4
− 6
+ ...
759375π 2025π 9π 3π N3
3 0.216152 4.342 121.673
= + − − + ... (2.6)
2 N N2 N3
and
We note that the coecients of the 1/N expansion are considerably larger than
in the d = 3 case.
4 Assuming the result (3.54) to order 1/N 3 , one nds that the
dimension of φi goes below unitarity at Ncrit = 35. This is much lower than the value
Ncrit = 1038 that we will nd in d = 6−. This is just a rough estimate, since the 1/N
expansion is only asymptotic and should be analyzed with care. Note for instance
that in the d = 3 case, a similar estimate to order 1/N 3 would suggest a critical value
model, where it is known [56, 17, 18, 57] that ∆φ ≈ 0.518 > 1/2. Nevertheless, the
reduction from a very large value Ncrit = 1038 in d = 6− to a much smaller critical
to occur in the Abelian Higgs model containing Nf complex scalars. A xed point in
23
d=4− exists only for Nf ≥ 183 [58, 3], while non-perturbative studies directly in
d=3 suggest a much lower critical value of Nf . 5 Some evidence for the reduction
corrections in [59]. It would be nice to study such corrections for the O(N ) theory
in 6− dimensions. It would also be very interesting to explore the 5d xed point at
nite N by numerical bootstrap methods similar to what has been done in d=3 in
[18]. For the non-unitary theory with N =0 such bootstrap studies were carried out
The large N critical O(N ) theory in general d was further studied in a series of
works by Lang and Ruhl [7, 8, 9, 10] and Petkou [11, 12]. Using conformal symmetry
and self-consistency of the OPE expansion, various results about the operator spec-
trum of the critical theory were derived. As an example of interest to us, [10] derived
an explicit formula for the anomalous dimension of the operator σk (the k -th power
gives
∆(σ k ) = 2k + (130k − 90k 2 ) + O(2 ) . (2.10)
N
For k = 2, 3, we will be able to match this result with the one-loop operator mixing
24
In [11, 12], explicit results for the 3-point function coecients gφφσ and gσ 3 were
gφφσ
hφi (x1 )φj (x2 )σ(x3 )i = 2∆ −∆
δ ij .
|x12 | φ σ |x23 |∆σ |x13 |∆σ
(2.11)
gσ3
hσ(x1 )σ(x2 )σ(x3 )i =
(|x12 ||x23 ||x13 |)∆σ
The coecient gφφσ was given in [11, 12] to order 1/N 2 and arbitrary d. Expanding
2 6 44 1
gφφσ = 1+ + O( 2 ) . (2.12)
N N N
This result indeed matches the value of the coupling g12 in (2.2) at the IR xed point,
as will be shown in the next section. To leading order in 1/N , the 3-point function
g2∗
which, as we will see, is precisely consistent with the ratio
g1∗
= 6 + O(1/N ) of the
sions of the fundamental elds φi and σ at the xed point. The theory (2.2) has an
25
d = 6, and have dimension 2
in d = 6 − . The Feynman rules of the theory are shown
in Figure 2.2.
δij
i j 1
φ p2 σ
p2
It is not hard to compute the one-loop beta functions β1 , β2 for the couplings
that Gm,n denote the Green's function with mφ elds and nσ elds. The one-loop
Z
2 dd k 1 1 2 p2 g12 Γ(3 − d/2)
D1 = (−g1 ) = −g1 I1 = − . (2.14)
(2π)d (p + k)2 k 2 (4π)3 6 (M 2 )3−d/2
Here we have used the renormalization condition p2 = M 2 . This has a 1/ pole in
For δσ , we have two one-loop diagrams (2 and 3). However, other than having
dierent coupling constant factors, the integrals are identical to diagram 1. Note that
7 We will state approximate expressions for the one-loop integrals I and I that are sucient for
1 2
extracting the log M 2 terms in d = 6 − . The more precise expressions are given in Appendix A.
26
G2,0 = + 1 + + ···
−p2 δφ
G0,2 = + + + ···
3
−p2 δσ
G2,1 = + + + ···
5 −δg1
G0,3 = + + + ···
7 −δg2
these two diagrams have a symmetry factor of 2, and there is a factor of N associated
1 2 N 2 N g12 + g22
D2 + D3 = (−g2 ) I1 + (−g1 ) I1 = − I1 . (2.16)
2 2 2
27
Now let us compute corrections to the 3-point functions. Diagram 4 gives:
Z
2 dd k 1 1 1
D4 = (−g1 ) (−g2 )
(2π) (k − p) (k + q) k 2
d 2 2
Z
dd k
3 1 1 1
D5 = (−g1 )
(2π) (k − p) (k + q) k 2
d 2 2
Finally, we calculate δg2 . The diagrams have the same topology, and just dier in
the coupling factors. Also, in diagram 6, we have a factor of N from the φ loop. So
we nd
D6 + D7 = N (−g1 )3 I2 + (−g2 )3 I2
−(N g13 + g23 ) Γ(3 − d/2)
= −(N g13 + g23 )I2 = . (2.21)
2(4π)3 (M 2 )3−d/2
28
The Callan-Symanzik equation for the Green's function Gm,n is:
∂ ∂ ∂
(M + β1 + β2 + mγφ + nγσ )Gm,n = 0 (2.23)
∂M ∂g1 ∂g2
Let's rst apply it to G2,0 . Then, to leading order in perturbation theory, the Callan-
1 ∂ 1
− M δ φ + 2γφ = 0, (2.24)
p2 ∂M p2
1 ∂ 1 g12
γφ = M δφ = . (2.25)
2 ∂M (4π)3 6
Note that γ1 > 0 as long as the coupling constant g1 is real. Analogously, we obtain
1 ∂ 1 N g12 + g22
γσ = M δσ = . (2.26)
2 ∂M (4π)3 12
∂ 1
β1 = − g1 + M (−δg1 + g1 (2δφ + δσ )) , (2.27)
2 ∂M 2
where the rst term accounts for the bare dimension of g1 in d = 6 − . After some
29
Finally, let us compute β2 . The Callan-Symanzik equation gives:
∂ 1
β2 = − g2 + M (−δg2 + g2 (3δσ )) , (2.29)
2 ∂M 2
and so
As a check, let us note that for N =0 the beta function for g2 in d=6 reduces to
3g23
β2 = − , (2.31)
4(4π)3
which is the correct result for the single scalar cubic eld theory in d = 6.
The single scalar cubic eld theory in d = 6− has no xed points at real
coupling, due to the negative sign of the beta function (2.31). It has a xed point at
to the Yang-Lee edge singularity [37]. However, as we now show, for suciently large
N, our model has a stable interacting IR xed point. Note that for large N the beta
functions simplify to
N g13
β1 = − g1 + (2.32)
2 12(4π)3
−4N g13 + N g12 g2
β2 = − g2 + . (2.33)
2 4(4π)3
r
6(4π)3
g1∗ = (2.34)
N
g2∗ = 6g1∗ . (2.35)
30
It is straightforward to compute the subleading corrections at large N by solving the
exact beta function equations (4.5), (3.15) in powers of 1/N . This yields
r
6(4π)3 22 726 326180
g1∗= 1+ + 2 − + ... (2.36)
N N N N3
r
∗ 6(4π)3 162 68766 41224420
g2 = 6 1+ + + + ... . (2.37)
N N N2 N3
Note that the coecients in this expansion appear to increase quite rapidly. This
suggests that the large N expansion may break down at some nite N. Indeed, we
will see in Section 2.4 that this large N IR xed point disappears at N ≤ 1038 (the
coupling constants go o to the complex plane). For all values of N ≥ 1039, the
xed point has real couplings and is IR stable, see Section 2.4. At large N, the IR
∂βi
stability of the xed point can be seen from the fact that the matrix evaluated
∂gj
at the xed point has two positive eigenvalues. These eigenvalues are in fact related
to the dimensions of the two eigenstates coming from the operator mixing of σ3 and
We can now use the values of the xed point couplings to compute the dimensions
of the elementary elds φi and σ in the IR. From (2.25) and (2.26) we obtain
d 1 (g1∗ )2
∆φ = − 1 + γφ = 2 − +
2 2 (4π)3 6
44 1936
=2− + + 2 + + ... (2.38)
2 N N N3
and
d 1 N (g1∗ )2 + (g2∗ )2
∆σ = − 1 + γσ = 2 − +
2 2 (4π)3 12
40 6800
=2+ + + ... (2.39)
N N2
31
Note that both dimensions are above the unitarity bound in d = 6 − , namely
d
∆ > 2
− 1, since the anomalous dimensions are positive. Moreover, note that the
order term in the dimension of σ cancels out and we nd ∆σ = 2 + O(1/N ). This is
in perfect agreement with the large N description of the critical O(N ) CFT. In that
goes from ∆ = d−2 at the free xed point to d − ∆ + O(1/N ) = 2 + O(1/N ) at the
interacting xed point. Furthermore, comparing (2.38), (2.39) with (2.4), (2.5), we
see that the coecients of the 1/N expansion precisely match the available results
for the large N critical O(N ) theory [5, 4, 6] expanded at d = 6 − . This is a strong
test that the IR xed point of the cubic theory in d = 6− is identical at large N to
We may also match the values of the xed point couplings (2.36),(2.37) with the
large N results (2.12), (2.13) for the 3-point functions coecients in the critical O(N )
theory. At leading order in , the 3-point functions in our cubic model simply come
from a tree level calculation, and it is straightforward to verify the agreement of (2.36)
with (2.12),
9 and that the ratio g1∗ /g2∗ = 6 at leading order at large N agrees with
(2.13).
r
6(4π)3
g1 ≡ x, (2.40)
r N
6(4π)3
g2 ≡ y. (2.41)
N
9 An overall normalization factor comes from the normalization of the massless scalar propagators
in d = 6.
32
After this, the vanishing of the β -functions (4.5), (3.15) gives
N y = −12N x3 + 3N x2 y − 9y 3 . (2.43)
them is the trivial solution (0, 0). Two of them are purely imaginary; they occur at
q
N
(0, ± 9
i). There are no simple expressions for the remaining six solutions, but it is
straightforward to study them numerically.
Two of them are real solutions in the second and fourth quadrant: (−x1 , y1 ) and
(x1 , −y1 ), with x1 , y1 > 0. They exist for any N, but they are not IR stable (they are
saddles, i.e. there is one positive and one negative direction in the coupling space).
They can be seen in both the bottom left and right graphs in Figure 2.4, corresponding
The behavior of the remaining four solutions changes depending on the value of N.
For N ≤ 1038, we nd that all four solutions are complex. The bottom right graph
of Figure 2.4 shows that for N = 500 there are indeed only the two real solutions
For N ≥ 1039, all four of these solutions become real, and they lie in the rst and
third quadrants: they have the form (x3 , y3 ), (x4 , y4 ), (−x3 , −y3 ), (−x4 , −y4 ), with
the bottom left graph of Figure 2.4 we plot the zeroes of β1 , β2 for N = 2000, with
regions of their signs labeled. Combining these, we can get the ow directions in each
of these regions in the bottom left graph. We nd that for all N ≥ 1039, we have two
stable IR xed points that are symmetric with respect to the origin, and are labelled
as red dots in the gure. These correspond to the large N solution (2.36)-(2.37), and
its equivalent solution with opposite signs of the couplings. For very large N, we see
33
that these solutions satisfy g2∗ = 6g1∗ as in (2.35). The origin is a UV stable xed
point, and from the direction of the renormalization group ow we can see that all
other xed points are saddle points: they have one stable direction and one unstable
direction.
It is interesting to treat N as a continuous parameter and solve for the value Ncrit
where the real IR stable xed points disappear. First, we notice that we can factor
y
x
times (2.42), we will cancel the Ny term, making the second beta function equation
x = α, y = β + 6α . (2.46)
(N − 44)α2 + β 2 = N (2.47)
α3 α2 α
840 3
+ (464 − N ) 2
+ 84 + 5 = 0 . (2.48)
β β β
Ncrit occurs when the curves dened by the two equations are tangent to each other.
Notice that we have eectively decoupled the two equations, since if we write z = α/β ,
we get
α2 (N − 44 + z −2 ) = N (2.49)
34
We can just rst solve the second equation, then easily solve the rst. To determine
the critical N, we want the second equation to have exactly one real root, thus we
For completeness, we now discuss the large N behavior of the four xed points
which are not IR stable (two of them are present for any N, and two of them only
√
O( N ). Then, at the leading order in N , we get:
N x = N x3 + xy 2 , N y = 3N x2 y − 9y 3 . (2.55)
q q √ q q √ q q √
This can be solved to get four solutions: ( 6 , 6 N ), (− 6 , 6 N ), ( 6 , − 16 N ),
5 1 5 1 5
q q √
5 1
(− 6
,− 6
N ). They correspond to the four real xed points at large N that are
35
Figure 2.4: (a) zeros of β1 at N = 2000. (b) zeros of β2 N = 2000. (c)
at RG ow
saddle points. Using the equation from β1 , we can also check that:
N x2 + y 2
∆σ = 2 − +
2 2 N
√
5 1
= 2 ± √ + O( ) . (2.56)
N N
36
2.5 Operator mixing and anomalous dimensions
1 1
= hφ(p)φ(q)φ2 (p + q)i = p2 q 2
1 1
= hσ(p)σ(q)σ 2 (p + q)i = p2 q 2
Let OiM denote the operator renormalized at scale M, and Oi0 denote the bare
operator. We are looking for expressions similar to (18.53) of [62], i.e. of the form:
The δ ij counterterms in the above equations are obtained by extracting the loga-
rithmic divergence from the diagrams shown in Figure 2.6, while the δφ O10 and δσ O20
terms correspond to external leg corrections, which are omitted in Figure 2.6.
37
= + + + +
δ 11 δ 12
= + + + +
δ 21 δ 22
Each of the terms δ ij is given by canceling the divergent pieces of two of the above
diagrams, as shown. For example, let us compute the rst diagram of the δ 11 counter
Z
2 dd k 1 1 1
D1 = (−g1 )
(2π) (k − p) (k + q) k 2
d 2 2
Next, we compute the second diagram of the δ 11 counter term. Notice that there is a
Z
N 2 dd k 1 1 1
D2 = (−g1 )
2 (2π) (p + k) k 2 p2
d 2
11 g12 N g12 Γ(3 − d/2)
δ = − + . (2.61)
2(4π)3 12(4π)3 (M 2 )3−d/2
Similarly, we can calculate the other three δ ij . The integrals are the same, only
the factors of coupling constants and N (due to closed φ loops) are dierent. With
√
our normalization convention of O1 , we need to multiply δ 12 by a factor of N, and
38
√
divide δ 21 by a factor of N. Then the matrix is symmetric, δ 21 = δ 12 , and we nd
√ √ !
12 N g12 N g1 g2 Γ(3 − d/2)
δ = − + , (2.62)
2(4π)3 1(4π)3 (M 2 )3−d/2
22 g22 g22 Γ(3 − d/2)
δ = − + . (2.63)
2(4π)3 12(4π)3 (M 2 )3−d/2
√ 2 √
1 Γ(3 − d/2) + −6g12 N g12 −6 N g1 + N g1 g2
δ ij = √ √ . (2.64)
12(4π)3 (M 2 )3−d/2 −6 N g12 + N g1 g2 2 2
−6g2 + g2
∂
γ ij = M (−δ ij + δzij ) , (2.65)
∂M
0 δσ
√ 2 √
−1 − 4g12 N g12 6 N g 1 − N g 1 g2
γ ij = 3 √ √ . (2.67)
6(4π) 6 N g12 − N g1 g2 4g22 − N g12
The eigenvalues of this matrix will give the dimensions of the two eigenstates arising
39
After plugging in the coupling constants at the IR xed point from (2.36) and
(2.37), and keeping its entries to order 1/N , the matrix elements of γ ij become:
√ 2 √
−1 4g12
− N g12 6 N g1 − N g1 g2
γ ij = √ √ (2.68)
6(4π)3 6 N g12 − N g1 g2 4g22 − N g12
40 840
1+ N N 3/2
= . (2.69)
840 100
N 3/2
1− N
To order 1/N , the o-diagonal terms do not aect the eigenvalues, and we get the
scaling dimensions
100 1
∆− = d − 2 + γ− = 4 − +O
N N2
(2.70)
40 1
∆+ = d − 2 + γ+ = 4 + +O
N N2
√
+ N 1
O = O + O2 (2.71)
6
−6
O− = √ O1 + O2 . (2.72)
N
Satisfyingly, we see that to this order the dimension ∆− precisely matches the dimen-
sion of the only primary eld of dimension near 4 in the large N UV xed point of
the quartic theory. It was determined using large N methods in [10], and corresponds
Comparing (4.26) with (2.39), we see that to this order ∆+ = 2 + ∆σ , so that the
40
We can actually show that ∆+ = 2 + ∆σ to all orders in 1/N (and for all xed
points) just using the β -function equations. To start, we redene our variables as
(2.40), (2.41). With this denition, the condition satised by x and y at the IR xed
8 2 12 1
1 = x2 − x − xy + y 2 (2.73)
N N N
x3 9
1 = −12 + 3x2 − y 2 . (2.74)
y N
4 2 √1 (xy
x2 − N
x N
− 6x2 )
γ ij = , (2.75)
√1 (xy − 6x )2 2 4 2
x − Ny
N
4 2
p
(2x2 − N
x − N4 y 2 ) ± |γ|
λ± = , (2.76)
2
144 4 48 3 4 16 32 16
|γ| = x − x y + x2 y 2 + 2 x 4 − 2 x2 y 2 + 2 y 4 . (2.77)
N N N N N N
Also, from the rst line of (2.39), we see that in terms of x and y the dimension ∆σ
becomes
2
∆σ = 2 − + x2 + y , (2.78)
2 2 2N
∆σ + 2 = 4 − + λ+ (2.79)
41
which, after some algebra, is seen to imply
4 2 5 p
1 − x2 + x + y 2 = |γ| . (2.80)
N N
4 2 12 6 p
− x − xy + y 2 = |γ| . (2.81)
N N N
Squaring both sides, and plugging in the expression for |γ|, we get that the following
1
36x4 − 12x3 y + x2 y 2 = 24x3 y + 32x2 y 2 − 36xy 3 + 5y 4 . (2.82)
N
To prove this equality, we again go back to (2.73) and (2.74). If we multiply both of
y2
these equations by 3xy − 2
and subtract the former from the latter, we get exactly
1 σ(φi )2
O = √ , O2 = σ 3 , (2.83)
N
(∂φi )2
O3 = √ , O4 = (∂σ)2 , (2.84)
N
φi φi
O5 = √ , O6 = σσ. (2.85)
N
42
However, in d = 6 − , O1 and O2 have bare dimensions 3 d−2
2
= 6 − 32 , while O3 , O4 ,
O5 , O6 have bare dimensions 2 + 2 d−2
2
= 6 − . Therefore, in d = 6 − , O1 and O2
mix with each other, and O3 , O4 , O5 , O6 mix with themselves.
expressions for the β functions, (4.5), (3.15), we can very simply compute the mixed
3N g12 −24g12 −24g1 g2 +g22 −12g12 +2g1 g2 √
− 2 + 12(4π)3 12(4π)3
3N
M ij = . (2.86)
−12g12 +2g1 g2 √ N g12 −9g22
12(4π)3
3N − 2 + 4(4π)3
We plug in the xed point value of g1 and g2 from (2.36) and (2.37) to get:
√
√ 3
840
− 5040
N2
+ ... N N
+ ...
M ij = √ . (2.87)
840
√ 3 + ... 420 154560
N N
− N
− N2
+ ...
Notice that at this IR xed point, both eigenvalues of M ij are positive, which implies
3780 1
λ1 = − + O( ) (2.88)
N2 N3
420 155820 1
λ2 = − − 2
+ O( 3 ) . (2.89)
N N N
The relation between scaling dimensions and the eigenvalues of the matrix is given
by
10
∆ = d + λ. (2.90)
10 As a simple test of this formula, note that at the free UV xed point the eigenvalues of (2.86)
are λ1 = λ2 = − 2 , giving dimensions ∆1 = ∆2 = 3(d/2 − 1) as it should be.
43
Thus, the scaling dimensions of the mixture of the two operators are:
1
∆1 = d + λ1 = 6 + O( ) (2.91)
N2
420 1
∆2 = d + λ2 = 6 − + O( 2 ) . (2.92)
N N
check, we note that ∆2 agrees with the dimension of the k = 3 primary operator
one positive and one negative eigenvalues. If we use (2.40), (2.41) and plug them into
(2.86), we get
√
− 2 +
2N
(3N x2 2
− 24x − 24xy + y ) 2
2N
(−12x2 + 2xy) 3N
M ij = √ . (2.93)
2N
(−12x 2
+ 2xy) 3N − 2 + 2N
3
(N x2 − 9y 2 )
p p √
x = ± 5/6 , y=± 1/6 N . (2.94)
Plugging these in, we nd that the eigenvalues of M ij at these xed points are:
5
λ1 = , λ2 = − , (2.95)
3
which conrms our graphical analysis that all of these xed points are saddle points.
44
2.6 Comments on CT and 5-d F theorem
A quantity of interest in a CFT is the coecient CT of the stress-tensor 2-point
Iµν,ρσ (x12 )
hTµν (x1 )Tρσ (x2 )i = CT , (2.96)
x2d
12
where Iµν,ρσ (x12 ) is a tensor structure uniquely xed by conformal symmetry, see e.g.
[63]. For N free real scalar elds in dimension d with canonical normalization, one
Nd
has CT = (d−1)Sd2
, where Sd is the volume of the d-dimensional round sphere. In the
critical O(N ) theory, using large N methods one nds the result [11, 12]
Nd 1 1
CT = 1 + CT,1 + O( 2 )
(d − 1)Sd2 N N
2
(2.97)
2C(µ) µ + 3µ − 2 d
CT,1 = − + 2
η1 , µ=
µ+1 µ(µ − 1) 2
Γ0 (x)
where C(µ) = ψ(3 − µ) + ψ(2µ − 1) − ψ(1) − ψ(µ), and ψ(x) = Γ(x)
. It is interesting
in (1.21), is of order , but there is a pole in C(µ) from the term ψ(3 − µ) ∼ −2/.
Hence, one nds at d=6−
7
CT,1 = 1 − + O(2 ) (2.98)
4
and so
d 7
CT = N + 1 − + O(2 ) . (2.99)
(d − 1)Sd2 4
a nice check on our description of the critical O(N ) theory via the IR xed point of
(2.2). The O() correction may be calculated in this description as well, using the
1-loop xed point (2.36), (2.37), but we leave this for future work.
45
The appearance of an extra massless scalar eld as d → 6 is also suggested by
for σ (as well as the σ3 coupling) become classically marginal. Note also that the
negative sign of the order correction in (2.99) implies that CT decreases from the
UV xed point of N + 1 free elds to the interacting IR xed point, consistently with
the idea that CT may be a measure of degrees of freedom [12]. However, expanding
CT in d = 4 + , one nds
d 5 2
CT = N − + O(3 ) . (2.100)
(d − 1)Sd2 12
with the UV xed point of the −φ4 theory, while the IR xed point corresponds to
N free scalars. Then, (2.100) is seen to violate CTU V > CTIR . A plot of CT,1 in the
range 2<d<6 is given in Figure 2.7. Finally, let us quote the value of CT in the
5 1408
CTd=5 = N− . (2.101)
4S52 1575π 2
Again, we observe that this value is consistent with CTU V > CTIR for the ow from the
free UV theory of N + 1 massless scalars to the interacting xed point, but not for the
ow from the interacting xed point to the free IR theory of N massless scalars. Thus,
It was proposed in [50] (see also [65, 66]) that, for any odd dimensional Euclidean
d+1 d−1
F̃ = (−1) 2 F = (−1) 2 log ZS d , (2.102)
46
CT,1
1.0
0.5
d
3 4 5 6
- 0.5
- 1.0
Figure 2.7: The O(N 0 ) CT,1 in the coecient CT of the stress tensor two-
correction
point function in the critical O(N ) theory for 2 < d < 6; see (2.97). Note that CT,1
crit
is negative for 2 < d . 5.22, so that in this range of dimensions CT < N CTfree sc.
for large N . Therefore, the CT theorem is respected in 2 < d < 4, but violated for
4 < d . 5.22.
Here F = − log ZS d is the free energy of the CFT on the round sphere, which is a
nite number after the power law UV divergences are regulated away (for instance,
by ζ -function or dimensional regularization). Note that when we put the CFT on the
sphere, we have to add the conformal coupling of the scalar elds to the Sd curvature.
This eectively renders the vacuum metastable, both in the case of the −φ4 theory
and in the cubic theory. Similarly, the CFT is metastable on R × S d−1 , which is
For d=3 (where F̃ = F ), the F-theorem (2.103) was proved in [67]. However,
higher dimensional versions are less well established. Using the results of this chapter
we can provide a simple new test of the 5d version of the F-theorem (in this case
the IR xed point of the cubic theory (2.2) or the UV xed point of the quartic scalar
47
theory. This implies that F̃ should satisfy
where F̃free sc. is minus the free energy of a 5d free conformal scalar [50]
leading O(N ) term is the same as in the free theory, while the O(N 0 ) term arises from
the determinant of the non-local kinetic operator for the auxiliary eld [24, 68, 50].
3ζ(5) + π 2 ζ(3) 1
F̃crit. = N F̃free sc. + +O . (2.106)
96π 2 N
The same answer may be obtained by computing the ratio of determinants for the
bulk scalar eld in the AdS6 Vasiliev theory with alternate boundary conditions [23].
We see that the left inequality in (2.104) is satised, since the O(N 0 ) correction
in (2.106) is positive (this check was already made in [23]). More non-trivially, we
3ζ(5)+π 2 ζ(3)
observe that the right inequality also holds, because
96π 2
' 0.001601 is smaller
than the value of F̃ for one free scalar, eq. (2.105).
48
2.7 Gross-Neveu-Yukawa model and a test of 3-d F-
theorem
The action of the Gross-Neveu (GN) model [25] is given by
Z
d i/ i 1 i i 2
S(ψ̄, ψ) = − d x ψ̄ ∂ ψ + g(ψ̄ ψ ) . (2.107)
2
Here N = Ñ tr1, where tr1 is the trace of the identity in the Dirac matrix space, and
the actual number of fermion components, and it is the natural parameter to write
down the 1/N expansion (factors of tr1 never appear in the expansion coecients if
The beta functions and anomalous dimensions of this model in d = 2+ can be
g2 g3
β(g) = g − (N − 2) + (N − 2) 2 + O(g 4 ) (2.108)
2π 4π
N − 1 2 (N − 1)(N − 2) 3
ηψ (g) = 2
g − 3
g + O(g 4 ) (2.109)
8π 32π
N −1 N −1 2
ηM (g) = g− g + O(g 3 ), (2.110)
2π 8π 2
can solve the beta function for the critical value of g at the xed point:
∗ 2π
g = 1− + O(3 ) . (2.111)
N −2 N −2
49
Plugging this value of g∗, we can nd the dimensions of the fermion eld and the
composite eld:
1+ N −1 2
∆ψ = d − 1 + ηψ (g ∗ ) = + 2
+ O(3 ) , (2.112)
2 4(N − 2)
∆σ = d − 1 − ηM (g ∗ ) = 1 − + O(2 ) . (2.113)
N −2
The UV xed point of the GN model in 2 < d < 4 dimensions is related to the IR
xed point of the Gross-Neveu-Yukawa (GNY) model, which has the following action
[26, 27, 3]
Z
d i 1 2
g2 4 i
S(ψ̄, ψ, σ) = d x −ψ̄ ∂/ + g1 σ ψ + (∂µ σ) + σ . (2.114)
2 24
In d = 4 − , the one-loop beta functions of the GNY model are given by:
N +6 4
βg12 = −g12 + g
2 1
(2.115)
16π
1 3 2 2 4
βg2 = −g2 + 2 g + N g2 g1 − 6N g1 . (2.116)
8π 2 2
For small , one nds that there is a stable IR xed point for any positive N:
16π 2
(g1∗ )2 = (2.117)
N +6
g2∗ = 16π 2 R, (2.118)
where
24N
R= √ . (2.119)
(N + 6)[(N − 6) + N 2 + 132N + 36]
1 2 N 2
γψ = g , γσ = g . (2.120)
32π 2 1 32π 2 1
50
and plugging in the value of the xed point couplings, one obtains
d−1 3 N +5
∆ψ = + γψ = − (2.121)
2 2 2(N + 6)
d−2 3
∆σ = + γσ = 1 − . (2.122)
2 N +6
arbitrary dimension following similar lines as in the scalar case reviewed in Section 2:
one introduces an auxiliary eld σ to simplify the quartic interaction, and develops
the 1/N perturbation theory with the eective non-local propagator for σ and the
σ ψ̄ψ interaction. The anomalous dimensions of ψ and σ in the critical theory have
been calculated respectively to order 1/N 3 and 1/N 2 , and in arbitrary dimension d
[53, 69]. To leading order in 1/N , the explicit expressions are given by
d−1 Γ (d − 1) sin( πd
2
) 1 1
∆ψ = − d
d
+O (2.123)
2 πΓ 2 − 1 Γ 2 + 1 N N2
4(d − 1)Γ (d − 1) sin( πd
2
) 1 1
∆σ = 1+ d
d
+O . (2.124)
π(d − 2)Γ 2 − 1 Γ 2 + 1 N N2
with (2.113) and (2.122) respectively. Note that both ∆ψ and ∆σ go below their
respective unitarity bounds for d > 4; so, unlike in the scalar case, there is no unitary
critical fermion theory above dimension four. One may check that the available
subleading large N results [53, 69] also precisely agree with the 1/N expansion of
(2.113) and (2.122), giving strong support to the fact that the 3d critical fermionic
CFT may be viewed as either the IR xed point of the GNY model, or the UV xed
We may use the two alternative descriptions of the critical fermion theory to
provide a simple test of the 3d F-theorem, similar to the tests carried out in [50]. In
the UV, the GNY model (which has relevant interactions in d = 3) is a free theory
51
of N fermions and one conformal scalar, while the GN model is free in the IR. Thus,
N Ffree fermi < Fcrit. fermi < N Ffree fermi + Ffree sc. . (2.125)
In d = 3, one nds
log 2 3ζ(3)
Ffree sc. = − ' 0.0638 . (2.126)
8 16π 2
The value of F in the critical theory may be computed perturbatively in the 1/N
expansion, and one obtains the result [68, 50]
ζ(3) 1
Fcrit. fermi = N Ffree fermi + +O . (2.127)
8π 2 N
ζ(3)
Because
8π 2
' 0.0152 < Ffree sc. , we see that indeed (2.125) is satised in both
directions.
52
Chapter 3
Igor Klebanov, and Grigory Tarnopolsky. We also thank J. Gracey and I. Herbut for
We would like to extend the one-loop analysis was carried out in 2 to three loops.
We have the cubic O(N ) symmetric theory of N + 1 scalar elds σ and φi in 6−
dimensions with the Lagrangian:
1 1 1 1
L = (∂µ φi )2 + (∂µ σ)2 + g1 σ(φi )2 + g2 σ 3 , (3.1)
2 2 2 6
The one-loop beta functions showed that for N > Ncrit there exists an IR stable
xed point with real values of the two couplings. We argued that this IR xed point
of the cubic O(N ) theory is equivalent to the perturbatively unitary UV xed point
of the O(N ) model with interaction (φi φi )2 , which exists for large N in 4<d<6
[32, 33, 34, 35]. The 1/N expansions of various operator scaling dimensions we found
53
in chapter 2 agree with the corresponding results [4, 5, 6, 7, 8, 9, 10, 11, 12] in the
A surprising result of chapter 2 was that the one-loop value of Ncrit is very large:
Ncrit , in section 3.2 we calculate the three-loop β functions, following the earlier work
N:
Ncrit = 1038.266 − 609.840 − 364.1732 + O(3 ). (3.2)
Neglecting further corrections, this gives Ncrit ( = 1) ≈ 64, but higher orders in can
obviously change this value signicantly. It is our hope that a conformal bootstrap
approach [13, 14, 15, 16], perhaps along the lines of [61], can help determine Ncrit
more precisely in d = 5. The bootstrap approach may also be applied in non-integer
dimensions close to 6, but one should keep in mind that such theories are not strictly
unitary [71].
2
known about the Abelian Higgs model in 4 − dimensions.3 For the model containing
Nf complex scalars, the one-loop critical value of Nf is found to be large, Nf,crit ≈ 183
[58]. However, the O() correction found from two-loop beta functions has a negative
coecient and almost exactly cancels the leading term when = 1, suggesting that
1 These papers considered cubic eld theories of q−1 scalar elds that were shown in [70] to
describe the q -state Potts model. These theories possess only discrete symmetries and generally
dier from the O(N ) symmetric theories that we study.
2 Another possible non-perturbative approach to the theory in 4<d<6 is the Exact Renormal-
ization Group [72]. This approach does not seem to indicate the presence of a UV xed point in the
theory of N scalar elds, but a search for an IR xed point in the theory of N +1 scalar elds has
not been carried out yet.
3 We are grateful to Igor Herbut for pointing this out to us.
54
Another interesting property of the theories (4.2) is the existence of the lower crit-
0 0
ical value Ncrit such that for N < Ncrit there is an IR stable xed point at imaginary
edge singularity in the two-dimensional Ising model (this is the (2, 5) minimal model
[73, 74] with central charge −22/5). From the three-loop β functions we nd the
expansion
0
Ncrit = 1.02145 + 0.03253 − 0.001632 + O(3 ) . (3.3)
0
The smallness of the coecients suggests that Ncrit > 1 for a range of dimensions
that it possesses an unstable xed point with g1∗ = g2∗ where the lagrangian splits
into that of two decoupled N = 0 theories. There is also an IR stable xed point
where g2∗ = 6g1∗ /5 + O(). A distinguishing feature of this non-unitary CFT is that
it has a discrete Z2 symmetry, and it would be interesting to search for it using the
In Section 3.4.1 we also discuss unstable unitary xed points that are present in
6− dimensions for all N . For N = 1 the xed point has g1∗ = −g2∗ ; it is Z3 symmetric
and describes the critical point of the 3-state Potts model in 6− dimensions [75].
4
Z
1 1 1 1
S= dd x (∂µ φi0 )2 + (∂µ σ0 )2 + g1,0 σ0 φi0 φi0 + g2,0 σ03 , (3.4)
2 2 2 6
4 We are grateful to Yu Nakayama for valuable discussions on this issue.
55
where φi0 and σ0 are bare elds and g1,0 and g2,0 are bare coupling constants.
5 As
1/2
σ0 = Zσ1/2 σ, φi0 = Zφ φi ,
(3.5)
g1,0 = µ 2 Zg1 Zσ−1/2 Zφ−1 g1 , g2,0 = µ 2 Zg2 Zσ−3/2 g2 .
Z
1 1 g1 g2
S= dd x (∂µ φi )2 + (∂µ σ)2 + σφi φi + σ 3 (3.7)
2 2 2 6
δφ δσ δg1 i i δg2 3
+ (∂µ φi )2 + (∂µ σ)2 + σφ φ + σ .
2 2 2 6
[76] in d = 6 − and employ the minimal subtraction scheme [77]. In this scheme, the
∞
X ∞
X ∞
X ∞
X
an (g1 , g2 ) bn (g1 , g2 ) z φ (g1 , g2 )
n z σ (g1 , g2 )
n
δg1 = , δg2 = , δφ = , δσ = .
n=1
n n=1
n n=1
n n=1
n
(3.8)
The anomalous dimensions and β -functions are determined by the coecients of the
simple 1/ poles in the counterterms [77]. Specically, in our case we have that the
5 We do not include mass terms as we are ultimately interested in the conformal theory. In the
dimensional regularization that we will be using, mass terms are not generated if we set them to
zero from the start.
56
anomalous dimensions are given by
1 ∂ ∂ φ
γφ = − g1 + g2 z , (3.9)
4 ∂g1 ∂g2 1
1 ∂ ∂ σ
γσ = − g1 + g2 z (3.10)
4 ∂g1 ∂g2 1
1 ∂ ∂ 1
β1 (g1 , g2 ) = − g1 + g1 + g2 − 1 (a1 − g1 (2z1φ + z1σ )),
2 2 ∂g1 ∂g2 2
1 ∂ ∂ 3
β2 (g1 , g2 ) = − g2 + g1 + g2 − 1 (b1 − g2 z1σ ). (3.11)
2 2 ∂g1 ∂g2 2
have to calculate the coecients of the 1/-divergencies in the loop diagrams, from
which we can read o the residues a1 (g1 , g2 ), b1 (g2 , g2 ), z1φ (g1 , g2 ), z1σ (g1 , g2 ).
Working in perturbation theory, we will denote by a1i the term of i-th order in
the coupling constants, and similarly for the other residue functions. Then, using the
results for the Feynman diagrams collected in the Appendix, we nd the anomalous
dimensions
1 φ φ 3 φ
γφ = − z12 − z14 − z16
2 2
g 2 2
g1
= 1 3− g 2
1 (11N − 26) − 48g1 g2 + 11g2
2
6(4π) 432(4π)6
g12
− g 4 (N (13N − 232) + 5184ζ(3) − 9064) + g13 g2 6(441N − 544)
31104(4π)9 1
− 2g12 g22 (193N − 2592ζ(3) + 5881) + 942g1 g23 + 327g24 ,
(3.12)
57
1 σ σ 3 σ
γσ = − z12 − z14 − z16
2 2
N g12 + g22 1 4 3 2 2 4
= + 2N g 1 + 48N g 1 g 2 − 11N g 1 g2 + 13g2
12(4π)3 432(4π)6
1
+ 96N (12N + 11)g15 g2 − 1560N g13 g23 + 952N g12 g24 (3.13)
62208(4π)9
− 2N g16 (1381N − 2592ζ(3) + 4280) + g26 (2592ζ(3) − 5195)
+ 3N g14 g22 (N + 4320ζ(3) − 8882)
1 φ σ 1 φ σ 1 φ σ
β1 = − g1 + (a13 − g1 (2z12 + z12 )) + 2(a15 − g1 (2z14 + z14 )) + 3(a17 − g1 (2z16 + z16 ))
2 2 2 2
1
= − g1 + 3
g1 (N − 8)g12 − 12g1 g2 + g22
2 12(4π)
1
− 6
g1 (536 + 86N )g14 + 12(30 − 11N )g13 g2 + (628 + 11N )g12 g22 + 24g1 g23 − 13g24
432(4π)
1
+ g1 g26 (5195 − 2592ζ(3)) + 12g1 g25 (−2801 + 2592ζ(3))
62208(4π)9
− 8g12 g24 (1245 + 119N + 7776ζ(3))
58
3 σ 3 σ 3 σ
β2 = − g2 + (b13 − g2 z12 ) + 2(b15 − g2 z14 ) + 3(b17 − g2 z16 )
2 2 2 2
1
= − g2 + 3
−4N g13 + N g12 g2 − 3g23
2 4(4π)
1
+ 6
−24N g15 − 322N g14 g2 − 60N g13 g22 + 31N g12 g23 − 125g25
144(4π)
1
+ − 48N (713 + 577N )g17 + 6272N g12 g25 + 48N g13 g24 (181 + 432ζ(3))
20736(4π)9
− 5g27 (6617 + 2592ζ(3)) − 24N g15 g22 (1054 + 471N + 2592ζ(3))
+ 2N g16 g2 (19237N − 8(3713 + 324ζ(3))) + 3N g14 g23 (263N − 6(7105 + 2448ζ(3)) .
(3.15)
In the case N =0 (the single scalar cubic theory), our results are in agreement with
r r
6(4π)3 6(4π)3
g1 ≡ x, g2 ≡ y. (3.16)
N N
59
In terms of the new variables x and y, the condition that both β -functions be zero
reads
1
0 = x(−8x2 + N (x2 − 1) − 12xy + y 2 )
2
1
− x (536 + 86N )x4 + 12(30 − 11N )x3 y + (628 + 11N )x2 y 2 + 24xy 3 − 13y 4
12N
1
− x 12xy 5 (2801 − 2592ζ(3)) + y 6 (−5195 + 2592ζ(3))
288N 2
+ 48x3 y 3 (2037 − 95N + 2592ζ(3))
and
1
0=− 9y 3 + N (12x3 + y − 3x2 y)
2
1
− 125y 5 + N x2 (24x3 + 322x2 y + 60xy 2 − 31y 3 )
4N
1
− N 2 x4 (27696x3 − 38474x2 y + 11304xy 2 − 789y 3 ) + 5y 7 (6617 + 2592ζ(3))
96N 2
+ 34224N x7 − 6272N x2 y 5 + 16N x6 y(3713 + 324ζ(3)) − 48N x3 y 4 (181 + 432ζ(3))
+ 48N x5 y 2 (527 + 1296ζ(3)) + 18N x4 y 3 (7105 + 2448ζ(3)) 2 . (3.18)
These equations can be solved order by order in the expansion. Using also the 1/N
expansion, we nd the xed point values
60
162 68766 41224420 28762554870
y∗ = 6(1 + + + + + ...
N N2 N3 N4
215 86335 75722265 69633402510
+ − − − − + ...
2N N2 N3 N4
2781 270911 − 157140ζ(3)
+ + + ... 2 . (3.20)
48N 6N 2
This large N solution corresponds to an IR stable xed point and generalizes the
previous one-loop result. This xed point exists and is stable to all orders in the 1/N
expansion.
If results beyond the 1/N expansion are desired, one can determine the the
expansions of x∗ , y∗ for nite N as follows. Plugging the expansions
into (3.17)-(3.18), the leading order terms are found to be [28, 78]
s s
N N
x0 (N ) = z(N ) , y0 (N ) = (1+6z(N )) ,
(N − 44)z(N )2 + 1 (N − 44)z(N )2 + 1
(3.22)
with large N behavior z(N ) = 840N + O(N 0 )6 . This solution is real only if N >
1038.27, as can be seen from the discriminant of the above cubic equation. Once the
x0 (N ), y0 (N ) are known, one can then determine the higher order terms in (3.21) by
61
For N 1038 the nite N exact results are close to (3.19)-(3.20), but for N ∼
1038 they deviate somewhat, indicating that, close to the critical N, the large N
expansion is not a good approximation (see also Figure 3.1 below).
read
x2 x2
γφ = − 2
(26 − 11N )x2 + 48xy − 11y 2 2
N 12N
x2
+ 3
6(544 − 441N )x3 y − 942xy 3 − 327y 4
144N
+ x4 (9064 + (232 − 13N )N − 5184ζ(3))
+ 2x2 y 2 (5881 + 193N − 2592ζ(3))} 3 , (3.24)
N x2 + y 2 1 4 2 2 2
2
γσ = − 13y + N x (2x + 48xy − 11y )
2N 12N 2
1
+ N 2 x4 (2762x2 − 1152xy − 3y 2 )
288N 3
+ 2N x2 −528x3 y + 780xy 3 − 476y 4 + 3x2 y 2 (4441 − 2160ζ(3)) + 8x4 (535 − 324ζ(3))
+ y (5195 − 2592ζ(3)) 3 .
6
(3.25)
62
Plugging the xed point values (3.19)-(3.20) into these expressions, we get the con-
d
∆φ = − 1 + γφ
2
1 44 1936 11 835 16352
=2− + + 2+ + ... + − − − + ... 2
2 N N N3 12N 6N 2 N3
13 6865 54367/2 − 3672ζ(3)
+ − + + + ... 3 , (3.26)
144N 72N 2 N3
d
∆σ = − 1 + γσ
2
40 6800 104 34190
=2+ + 2
+ ... + − − 2
+ ... 2
N N 3N 3N
22 47695/18 − 2808ζ(3)
+ − + 2
+ ... 3 . (3.27)
9N N
One can verify that these results are in precise agreement with the large N calculation
of [4, 5, 6] for the critical O(N ) model in general d, analytically continued to d = 6−.
This provides a strong check on our calculations and on our interpretation of the IR
The 1/N expansions are expected to work well for N 1038. For any N larger
than the critical value, the expansions of the scaling dimensions may be determined
using (3.25) and the exact analytic solutions for the xed point location (x∗ , y∗ ). For
example, in Figure 3.1 we plot the coecient of the O(3 ) term in ∆σ as a function
chapter, we carry out an additional check, comparing with the O(1/N 2 ) correction
63
gs at OHe3L
0.008
0.006
exact result
0.004 up to 1 N 8
up to 1 N 4
0.002
0.000 N
1200 1400 1600 1800 2000
-0.002
Figure 3.1: The O(3 ) in ∆σ as a function of N for N ≥ 1039. The 1/N expansion
found in [79], but still working to the one-loop order in (it should be straightforward
to generalize the mixing calculation to higher loops, but we will not do it here).
λ
In the quartic O(N ) theory with interaction
4
(φi φi )2 , the derivative of the beta
function at the xed point coupling
ω1 ω2
ω = β 0 (λ∗ ) = 4 − d + + 2 + ... (3.28)
N N
∆φ4 = d + ω . (3.29)
64
The coecient ω2 has a more complicated structure for general d which was rst
2259
ω2 = 102 + − + 120ζ(3) 3 + . . . , d=4− (3.32)
2
237476 2 92480
ω2 = −49760 + + − + 32616ζ(3) 3 + . . . , d=6− (3.33)
3 9
In d = 4 − , one can check that the above results correctly reproduce the derivative
N + 8 2 3(3N + 14) 3
β = −λ + λ − λ
8π 2 64π 4
33N 2 + 480N ζ(3) + 922N + 2112ζ(3) + 2960 4
+ λ + O(λ5 ) (3.34)
4096π 6
matched to the primary operator arising from the mixing of the σ2 and φi φi operators
in our cubic theory. In chapter 2, the mixing matrix of σ2 and φi φi to one-loop order
was found to be
√ 2 √
−1 4g12 − N g12 6 N g1 − N g1 g2
γ ij = √ √ . (3.36)
6(4π)3 6 N g12 − N g1 g2 4g22 − N g12
65
Computing the eigenvalues γ± of this matrix, and inserting the values of one-loop
r
6(4π)3 22 726 326180
g1∗ = 1+ + 2 − + ... , (3.37)
N N N N3
r
6(4π)3 162 68766 41224420
g2∗ =6 1+ + + + ... (3.38)
N N N2 N3
100 49760 27470080
∆− = d − 2 + γ− = 4 + − − − + . . . + O(2 ), (3.39)
N N2 N3
40 6800 2637760
∆+ = d − 2 + γ+ = 4 + + 2
+ 3
+ . . . + O(2 ). (3.40)
N N N
∆− is a primary, and comparing with (3.30), (3.33), we see that its dimension precisely
agrees with the results of [79] to order 1/N 2 . The higher order terms in can be
We now calculate the mixed anomalous dimensions of the nearly marginal oper-
ators O1 = σφφ and O2 = σ 3 . Using the beta functions written in equations (4.5)-
(3.15), we can determine the anomalous dimensions of the nearly marginal operators
∂βi
Mij = . (3.41)
∂gj
Strictly speaking, this matrix is not exactly equal to the anomalous dimension mix-
√
dividing and multiplying the o-diagonal elements by a factor 3N + O(), which
corresponds to an appropriate rescaling of the couplings. This clearly does not change
66
the eigenvalues of the matrix, and hence we can directly compute the eigenvalues λ±
of (3.41), and obtain the dimensions of the eigenstate operators as
∆± = d + λ± . (3.42)
Plugging in the xed point values x∗ and y∗ from equations (3.19)-(3.20), we nd that
155 2 1777 3 1 1
∆+ = 6 + − + ... + O( 2 ) ,
3 36 N N
(3.43)
1051 3 1 1
∆− = 6 + −420 + 4992 − + ... + O( 2 ).
12 N N
The dimension of the σk operator in the quartic O(N ) model is known to order 1/N
as function of d [10], and may be written as
k k(d − 2) (k − 1)d2 − d(3k − 1) + 4 Γ (d) 1
∆(σ ) = 2k + 3 + O( 2 ) . (3.44)
N dΓ 2 − d2 Γ d2 N
Our result for ∆− agrees with the expansion of this formula for k=3 in d = 6 − .
values of the couplings. This can be dened as the value of N (formally viewed as a
continuous parameter) at which two real solutions of the β -function equations merge,
and subsequently go o to the complex plane. Geometrically, this means that the
curves on the (g1 , g2 ) plane dened by the zeroes of β1 and β2 are barely touching, i.e.
they are tangent to each other. Therefore the critical N , as well as the corresponding
67
critical value of the couplings, can be determined by solving the system of equations
β1 = 0 , β2 = 0,
(3.45)
∂β1 /∂g1 ∂β2 /∂g1
= .
∂β1 /∂g2 ∂β2 /∂g2
Note that the condition in the second line is equivalent to requiring that the de-
∂βi
Mij = ∂gj
, vanishes. This means that one of the two eigenstates becomes marginal.
expansion
x = x0 + x1 + x2 2 + O(3 ),
y = y0 + y1 + y2 2 + O(3 ), (3.46)
N = N0 + N1 + N2 2 + O(3 )
and plugging this into (3.45), we can solve for the undetermined coecients uniquely.
The above system of equations can be solved analytically, as was done in chapter 2.
We nd that, up to the signs of x0 and y0 , there are three inequivalent solutions
68
y
40
20
IR
0
IR UV
-20
-40
x
-3 -2 -1 0 1 2 3
Figure 3.2: The zeroes of the one-loop β functions and the RG ow directions for
N = 2000. The red dots correspond to the stable IR xed points, while the black
dots are unstable xed points. As N → Ncrit , the red dot merges with the nearby
black dot, and the two xed points move into the complex plane.
The rst of these solutions, with Ncrit = 1038.26605 + O(), is of most interest to
us because it is related to the large N limit of the theory. For N > Ncrit , we nd
a stable IR xed point at real couplings g1 and g2 . 7 This xed point is shown with
the red dot in Figure 3.2 (there is a second stable IR xed point obtained by the
also a nearby unstable xed point, shown with a black dot, which has one stable and
one unstable direction. As N approaches Ncrit from above, the nearby unstable xed
7 It is stable with respect to ows of the nearly marginal couplings g and g . As usual, there are
1 2
some O(N ) invariant relevant operators that render this xed point not perfectly stable.
69
point approaches the IR stable xed point, and they merge at Ncrit . At N < Ncrit ,
both xed points disappear into the complex plane. As discussed in [80], this is a
rather generic behavior at the lower edge of the conformal window: the conformality
is lost through the annihilation of a UV xed point and an IR xed point. In [80] this
was argued to happen at the lower (strongly coupled) edge of the conformal window
that the same type of behavior occurs at the lower edge of the conformal window of
Let us identify the operator that causes the ow between the unstable xed point
and the IR stable xed point of our primary interest. It is one of the two nearly
marginal operators cubic in the elds that were studied in section 3.2. By studying
N = Ncrit and causes the ow between the IR xed point and the nearby UV xed
point for N slightly above Ncrit . In bootstrap studies of the quartic O(N ) model
this operator was denoted by σ3 [10], i.e. it can be thought of as the triple-trace
operator" (φi φi )3 . The theory at the unstable xed point has an unconventional large
√
N behavior where x ∼ O(1) and y ∼ O( N ), so that corrections to scaling dimension
proceed in powers of N −1/2 as in chapter 2.
Let us now go back to nding the higher order corrections to Ncrit given by (3.48)
(the higher order corrections to the other critical values (3.49)-(3.50) will be discussed
in the next section). Once we have solved the leading order system (3.47), we can plug
the solution into (3.46) and expand (3.45) up to order 2 . From this we obtain simple
nd
70
Thus, to three-loop order, we conclude that
We have also checked these expansion coecients via a direct high-precision numerical
calculation of Ncrit for very small values of . The large and negative coecients
lower than the zeroth order value (this is analogous to the result [59] for the Abelian
Higgs model). If we just use the rst three terms and plug in = 1, we get:
For N < Ncrit the anomalous dimensions, such as γφ , are no longer positive (in
fact, they become complex). This loss of positivity of γφ can also be seen as N is
reduced in the quartic O(N ) model. For example, using the 1/N expansion of γφ in
d=5 [6]
32 1427456
γφ = 2
−
15π N
3375π 4 N 2
275255197696 89735168 32768 ln 4 229376ζ(3) 1
+ 6
− 4
+ 4
− 6
+ ...
759375π 2025π 9π 3π N3
3 0.216152 4.342 121.673
= + − − + ... (3.54)
2 N N2 N3
we nd that it stops being positive for N < 35. This critical value is not too far from
(3.53).
71
It is also instructive to study the theory using the 4 − expansion. The anomalous
dimensions of φi is [2]
N +2 2 N +2 2
3
γφ = + −N + 56N + 272
4(N + 8)2 16(N + 8)4
N +2 4 3 2
4
+ −5N − 230N + 1124N + 17920N + 46144 − 384ζ(3)(5N + 22)(N + 8)
64(N + 8)6
+ O(5 ) (3.55)
For positive this expansion gives accurate information about the Wilson-Fisher IR
xed points [1]. For negative there exist formal UV xed points at negative quartic
coupling where we can apply this formula as well. In that case, γφ becomes negative
for suciently large || and N < Ncrit , indicating that the operator φi violates the
d=5 are even lower than the three-loop estimate (3.53). It seems safe to conclude
that the true value is much lower than the one-loop estimate of 1038. To determine
Ncrit in d=5 more precisely, one needs a non-perturbative approach to the d=5
theory, perhaps along the lines of the conformal bootstrap calculation in [61].
xed points that are located in the upper left and lower right corners of Figure 3.2
exist for all positive N, and we would like to nd their interpretation.
The xed point with N = 1 has a particularly simple property that g1∗ = −g2∗ .
This property of the solution holds for the three-loop β functions, and we believe that
it is exact. Using this, we note that the action at the xed point is proportional to
(σ + iφ)3 + (σ − iφ)3 . Therefore, the theory at this xed point enjoys a Z3 symmetry
72
acting by the phase rotation on the complex combination σ + iφ. This cubic classical
action appears in the Ginzburg-Landau theory for the 3-state Potts model (see, for
example, [75]).
8 Therefore, we expect the Z3 symmetric xed point to describe the
3-state Potts model in d = 6 − . The dimensions of operators at this xed point are
1 2 443 3
∆φ = ∆σ = 2 − + 2 + + O(4 ) . (3.56)
3 3 54
This is in agreement with the result of [39]. By calculating the eigenvalues λ± of the
∂βi
matrix Mij = ∂gj
, we also nd the dimensions (3.42) of the two cubic operators to
order 3 :
14 158 2 17380
∆− = 6 − − − + 16ζ(3) 3 = 6 − 4.66667 − 17.55562 − 233.8013 ,
3 9 81
83 2 38183
∆+ = 6 − − + 4ζ(3) 3 = 6 − 4.611112 − 63.73263 . (3.57)
18 648
the Z3 symmetry and is slightly irrelevant for small . The dimension ∆− corresponds
to the relevant operator σ(σ 2 +φ2 ) which breaks the Z3 . Thus, the Z3 symmetry helps
stabilize the xed point at small .
Unfortunately, the 6− expansions (3.57) have growing coecients, and it is
not clear for what range of the xed point exists. Thus, one may not be able to
interpolate smoothly from the Z3 symmetric xed point in d=6− to d=2 where
the 3-state Potts model is described by the unitary (5, 6) minimal model [73].
For all N ≥2 we nd unstable xed points with O(N ) symmetry. These xed
∂βi
the matrix
∂gj
. Also, they exhibit an unconventional large N behavior involving half-
integer powers of N, similarly to the unstable xed points that appear for N > Ncrit
8 We are grateful to Yu Nakayama for pointing this out to us.
73
and are shown by the black dots in the upper right and lower left corners of Figure
00 0
appear for N > Ncrit , there exist non-unitary xed points for Ncrit < N < Ncrit . The
0 00
leading values of Ncrit and Ncrit are given in (3.49) and (3.50), respectively. Using the
method developed above for nding the higher order in corrections to Ncrit we get
0
Ncrit = 1.02145 + 0.03253 − 0.001632
x0 = i 0.23185 + 0.08887 − 0.039562 , y 0 = i 0.25582 + 0.11373 − 0.042762
(3.58)
and
00
Ncrit = −0.08750 + 0.34726 − 0.882742
Unfortunately, the latter expansion has growing coecients, and we cannot extract
any useful information from it. On the other hand, the higher order corrections to
0 0
Ncrit are very small, which suggests that Ncrit >1 for range of dimensions below 6.
The theory with N = 0, which contains only the eld σ, was originally studied by
Michael Fisher as an approach to the Yang-Lee edge singularity in the Ising model [37].
Since the coupling is imaginary, it describes a non-unitary theory where some operator
dimensions (e.g. σ ) are below the unitarity bounds. In d = 2, this CFT corresponds to
the (2, 5) minimal model [74], which has c = −22/5. A conformal bootstrap approach
to this model [60] has produced good results for a range of dimensions below 6.
74
The N =1 theory, which has two elds and two coupling constants, has a more
intricate structure. This theory is distinguished from the N =0 case by the presence
σ1 = σ + φ , σ2 = σ − φ , (3.60)
we note that for g1∗ = g2∗ the interactions of the N = 1 model decouple as ∼ σ13 +σ23 , i.e.
at this xed point the theory is a sum of two Fisher's N =0 theories. However, one
of the ow directions at this xed point is unstable, since the corresponding operator
has ∆O = 6 − 10/9 + O(2 ) and is relevant (this value of the dimension corresponds
∂βi
to the negative eigenvalue of the matrix Mij = ∂gj
at the g1∗ = g2∗ xed point). This
dimension has a simple explanation as follows. The ow away from the decoupled
xed point is generated by the operator O = σ1 σ22 + σ2 σ12 . This is allowed by the
∆O = ∆N
σ
=0
+ ∆N
σ2
=0
= 2 + 2∆N
σ
=0
, (3.61)
5 43 2 8ζ(3) 8375
∆σ = 2− − + − 3 = 2−0.555556−0.02949252 +0.0218453
9 1458 243 472392
(3.62)
Substituting this into (3.61) we nd the dimension of the relevant operator O, which
∂βi
Mij = ∂gj
at the g1∗ = g2∗ xed point. Using the expansion (3.62), we nd that O
continues to be relevant as is increased. For = 4, i.e. d = 2, we know the exact
75
result in the (2, 5) minimal model that ∆N
σ
=0
= −2/5, which implies ∆O = 6/5. This
strongly suggests that O is relevant, and the decoupled xed point is unstable, for
the entire range 2 ≤ d < 6. To describe this CFT in d = 2 more precisely, we note the
existence of the modular invariant minimal model M (3, 10), which is closely related
The ow away from the unstable xed point with g1∗ = g2∗ can lead the N =1
theory to the IR stable xed point where g2∗ = 6g1∗ /5 + O(). Using our results we
can deduce the expansion of various operator dimensions at this xed point. For
example,
∂βi
By calculating the eigenvalues λ± of the matrix Mij = ∂gj
, we nd the dimensions of
As is increased, these expansions suggest that the two operators become more
while the N =0 theory has no symmetries at all. As we have noted, in d=2 the
CFT can be obtained via deforming the (3, 10) minimal model by a Virasoro primary
eld of dimension 6/5 (this is the highest dimension relevant operator in that minimal
76
model). After analyzing the spectra of several candidate minimal models, we suggest
that the end point of this RG ow is described by the (3, 8) minimal model with
c = −21/4.9 Let us note that M (2, 5) and M (3, 8) are members of the series of
6− dimensions, we can tentatively identify the leading Z2 odd operator as φ and
the leading Z2 even one as σ. Obviously, further work is needed to check if the stable
α
α β
= −dαβγ = δαβ p12
β γ
α
α β
= −(δg)αβγ = −p2 (δz)αβ
β γ
9 Note that this value is greater than the central charge of the UV theory M (3, 10), which is equal
to −44/5. For ows between non-unitary theories the Zamolodchikov c-theorem does not hold, and
it is possible that cU V < cIR .
77
where we introduced symmetric tensor coupling dαβγ and counterterms (δg)αβγ , (δz)αβ
with α, β, γ = 0, 1, .., N as
where i = 1, ..., N . The general form of a Feynman diagram in our theory could be
schematically represented as
The Tensor structure factors are products of the tensors dαβγ and (δg)αβγ , (δz)αβ ,
with summation over the dummy indices. Their values for dierent diagrams are
symmetry factors and are the same as in the usual ϕ3 -theory; their values are listed
10 To nd the Tensor structure factor we used the fact that it is a polynomial in N , so we
calculated sums of products of dαβγ , (δg)αβγ , (δz)αβ explicitly for N = 1, 2, 3, 4, ..., using Wolfram
Mathematica. Having answers for N = 1, 2, 3, 4, .. it's possible to restore the general N form.
78
3.6.1 Counterterms
φ g12 2 2
z14 = g1 (11N − 26) − 48g 1 g2 + 11g 2 ,
432(4π)6
σ 1 4 3 2 2 4
z14 =− 2N g 1 + 48N g 1 g 2 − 11N g 1 g 2 + 13g 2 ,
432(4π)6
1 2 3 2 2 3
a15 =− g 1 g 1 (11N + 98) − 2g 1 g2 (7N − 38) + 101g 1 g 2 + 4g2 ,
144(4π)6
1 5 4 3 2 2 3 5
b15 =− 4N g 1 + 54N g 1 g 2 + 18N g 1 g 2 − 7N g 1 g2 + 23g 2 ,
48(4π)6
(3.68)
φ g12
z16 = g 4 (N (13N − 232) + 5184ζ(3) − 9064) + g13 g2 6(441N − 544)
46656(4π)9 1
− 2g12 g22 (193N − 2592ζ(3) + 5881) + 942g1 g23 + 327g24 ,
σ 1
z16 =− 9
2N g16 (1381N − 2592ζ(3) + 4280) − 96N (12N + 11)g15 g2
93312(4π)
− 3N g14 g22 (N + 4320ζ(3) − 8882) + 1560N g13 g23 − 952N g12 g24 − g26 (2592ζ(3) − 5195) ,
g12
a17 = 9
− g15 (N (531N + 10368ζ(3) − 2600) + 23968)
15552(4π)
+ g14 g2 (99N 2 + 2592(5N − 6)ζ(3) − 9422N − 2588) + 2g13 g22 (1075N + 2592ζ(3) − 16897)
+ 2g12 g23 (125N − 5184ζ(3) − 3917) − g1 g24 (5184ζ(3) + 721) + g25 (2592ζ(3) − 2801) ,
1
b17 = − 2g17 N (577N + 713) − 48g16 g2 N (31N − 59)
2592(4π)9
+ g15 g22 N (423N + 2592ζ(3) + 1010) − g14 g23 N (33N − 1296ζ(3) − 6439)
− 27g13 g24 N (32ζ(3) + 11) − 301N g12 g25 + g27 (432ζ(3) + 1595) . (3.69)
79
Chapter 4
Igor Klebanov, and Grigory Tarnopolsky. We also thank D. Anninos, D. Harlow and
The O(N ) invariant theories of N massless scalar elds φi , which interact via
λ
the potential
4
(φi φi )2 possess interacting IR xed points in dimensions 2 < d <
4 for any positive N [1]. These theories contain O(N ) singlet current operators
with all even spin, and when N is large the current anomalous dimensions are ∼
1/N [2]. Since all the higher spin currents are nearly conserved, the O(N ) models
possess a weakly broken higher spin symmetry. In the Anti-de Sitter/Conformal Field
Theory (AdS/CFT) correspondence [84, 85, 86], each spin s conserved current in a d
dimensional CFT is mapped to a massless spin s gauged eld in d+1 dimensional
AdS space. For these reasons, it was conjectured [47] that the singlet sector of the
critical O(N ) models in d=3 is dual to the interacting higher spin theory in d=4
80
containing massless gauge elds of all even positive spin [41, 42, 43, 46]. This minimal
Vasiliev theory also contains a scalar eld with m2 = −2/`2AdS , and the two admissible
boundary conditions on this eld [48, 49] distinguish the interacting O(N ) model from
the free one (in the latter all the higher spin currents are conserved exactly). A review
A remarkable feature of the Vasiliev theories [41, 42, 43, 46] is that they are
consistent not only in Anti-de Sitter, but also in de Sitter space. On general grounds
one expects a CFT dual to quantum gravity in dS4 to be a non-unitary theory dened
in three dimensional Euclidean space [88]. In [89] it was proposed that the CFT dual
to the minimal higher spin theory in dS4 is the theory of an even number N of
Z
1 λ
S= 3
d x Ωij ∂µ χi ∂ µ χj + (Ωij χi χj )2 . (4.1)
2 4
This theory possesses Sp(N ) symmetry, and Ωij is the invariant symplectic matrix.
This model was originally introduced and studied in [90, 91], where it was shown to
possess an IR xed point in 4− dimensions. The beta function of this model is
related to that of the O(N ) model via the replacement N → −N . According to the
proposal of [89], the free UV xed point of (4.1) is dual to the minimal higher spin
theory in dS4 with the Neumann future boundary conditions on the m2 = 2/`2dS bulk
scalar eld, and its interacting IR xed point to the same higher spin theory but with
the Dirichlet boundary conditions on the scalar eld. In the latter case, the higher
spin symmetry is slightly broken at large N, and the de Sitter higher spin gauge
1 A potential diculty with this picture is that unitarity in dSd+1 space requires that a massive
eld of spin s>1 should satisfy m2 > (s − 1)(s + d − 3)/`2dS [92, 93, 94]. In other words, there is a
nite gap between massive elds and massless ones (the latter are dual to exactly conserved currents
in the CFT). However, since the masses are generated by quantum eects and are parametrically
small at large N, this is perhaps not fatal for bulk unitarity. It would be interesting to clarify this
further.
81
of the de Sitter boundary conditions from the point of view of the wave function of
of the higher spin AdS/CFT. It was observed long ago [32, 33, 34] that in d>4 the
quartic O(N ) models possess UV xed points which can be studied in the large N
expansion. The UV completion of the O(N ) scalar theory in 4 < d < 6 was proposed
in [28]; it is the cubic O(N ) symmetric theory of N +1 scalar elds σ and φi :
Z
1 1 1 1
S= d x (∂µ φi )2 + (∂µ σ)2 + g1 σ(φi )2 + g2 σ 3
d
. (4.2)
2 2 2 6
For suciently large N, this theory has an IR stable xed point in 6− dimensions
with real values of g1 and g2 . The beta functions and anomalous dimensions were
1/N expansion of the quartic O(N ) model at its UV xed point [4, 6, 7, 10, 11] when
the quartic O(N ) model is continued to 6− dimensions. The conformal bootstrap
approach to the higher dimensional O(N ) model was explored in [61, 20, 98].
To extend the idea of [89] to dSd+1 with d > 4, we may consider non-unitary
CFTs (4.1) which in d>4 possess UV xed points for large N. The 1/N expansion
Stratonovich transformation, and one nds that it is related to the 1/N expansion in
the O(N ) models via the replacement N → −N . In d = 5 this interacting xed point
should be dual to the higher spin theory in dS6 with Neumann boundary conditions
82
2 + O(1/N ) on the CFT side).3 In search of the UV completion of these quartic CFTs
in 4 < d < 6, we introduce the cubic theory of one commuting real scalar eld σ and
Z
1 1 1 1
S= d x Ωij ∂µ χi ∂ µ χj + (∂µ σ)2 + g1 Ωij χi χj σ + g2 σ 3
d
. (4.3)
2 2 2 6
Alternatively, we may combine the elds χi into N/2 complex anti-commuting scalars
θα , α = 1, . . . , N/2 [90, 91]; then the action assumes the form
Z
d α µ α 1 2 α α 1 3
S= d x ∂µ θ ∂ θ̄ + (∂µ σ) + g1 σθ θ̄ + g2 σ . (4.4)
2 6
We study the beta functions for this theory in d = 6− and show that they are
related to the beta functions of the theory (4.2) via the replacement N → −N . For
all N there exists an IR xed point of the theory (4.3) with imaginary values of g1
and g2 . This is similar to the IR xed point of the single scalar cubic eld theory
(corresponding to the N =0 case of our models), which was used by Fisher [37] as
an approach to the Lee-Yang edge singularity. The fact that the couplings are purely
imaginary makes the integrand of the path integral oscillate rapidly at large σ; this
should be contrasted with real couplings giving a potential unbounded from below.
Our results allow us to study the 6 − expansion of the theory (4.3) with arbtrary
N, and we observe that at nite N there are qualitative dierences between the
Sp(N ) and O(N ) models which are not seen in the 1/N expansion. In fact, for the
Sp(N ) model there is no analogue of the lower bound Ncrit that was found in the
O(N ) case [28]. For the lowest value, N = 2, we observe some special phenomena.
formulate the quartic interaction (4.1); thus, the cubic lagrangian (4.20) seems to be
3 Thed = 5 free Sp(N ) model corresponds to Dirichlet scalar boundary conditions in the dS6 dual.
It should be possible to extend this free Sp(N ) model/dS higher spin duality to general dimensions,
since the Vasiliev equations in (A)dSd+1 are known for all d [46].
83
the only possible description of the interacting theory with global Sp(2) symmetry.
implies that the scaling dimensions of σ and χi are equal, and we check this to order
dimensions in our OSp(1|2) symmetric theory are the same as in the q→0 limit of
the q -state Potts model.5 This provides strong evidence that the OSp(1|2) symmetric
IR xed point of the cubic theory (4.20) describes the second order transitions in the
Using the results of [78], we also compute perturbatively the sphere free energies of
the models (4.3). The sphere free energy for the OSp(1|2) symmetric model is found to
be the same as for the q = 0 Potts model. In terms of the quantity F̃ = sin( πd
2
) log ZS d ,
which was introduced in [78] as a natural way to generalize the F -theorem [66, 54, 50]
to continuous dimensions, we nd that the RG ow in the cubic Sp(N ) models in
d = 6− satises F̃UV > F̃IR for all N ≥ 2. We show that the same result holds in
CFTs the inequality F̃UV > F̃IR is not always satised. It would be interesting to
to the fact that these models are presumably dual to unitary higher spin gravity
4 We are grateful to Giorgio Parisi for informing us about this and for important discussions.
5 We thank Sergio Caracciolo for suggesting this comparison to us and for informing us about the
paper [100].
6 The models (4.1) and (4.3) also satisfy the F -theorem to leading order in the large N expansion
(for all d), since the leading order correction to F̃UV − F̃IR is of order N 0 , and was shown to satisfy
the F -theorem in the corresponding unitary O(N ) models [50, 23, 78].
84
4.2 The IR xed points of the cubic Sp(N ) theory
The beta functions and anomalous dimensions for the Sp(N ) symmetric model (4.3)
O(N ) model (4.2), which were computed in [28, 29] to three-loop order. Indeed,
writing the action in the complex basis (4.4), we see that the Feynman rules and
propagators are identical to those of the O(N ) theory written in the U (N/2) basis,
the only dierence being that the N/2 complex scalars are anticommuting. Hence, for
each closed loop of the θα we get an extra minus sign, thus explaining the replacement
N → −N .
Using the results in [28, 29], the beta functions for the Sp(N ) model are then
found to be:
1 2 2
β1 = − g1 − g1 (N + 8)g 1 + 12g 1 g2 − g 2
2 12(4π)3
1
− 6
g1 (536 − 86N )g14 + 12(30 + 11N )g13 g2
432(4π)
+ (628 − 11N )g12 g22 + 24g1 g23 − 13g24 + . . . ,
1
β2 = − g2 + 3
4N g13 − N g12 g2 − 3g23
2 4(4π)
1
+ 6
24N g15 + 322N g14 g2 + 60N g13 g22 − 31N g12 g23 − 125g25 + . . . . (4.5)
144(4π)
We have omitted the explicit three-loop terms, which can be obtained from [28, 29].
Similarly, the anomalous dimensions of the elds χi and σ take the form
g12 g12 2 2
γχ = + g1 (11N + 26) + 48g 1 g 2 − 11g2 + ...,
6(4π)3 432(4π)6
(4.6)
N g12 − g22 1 4 3 2 2 4
γσ = − − 2N g 1 + 48N g 1 g 2 − 11N g 1 g2 − 13g 2 + ....
12(4π)3 432(4π)6
With the beta functions at hand, we can look for non-trivial xed points of the RG
ow satisfying β1 (g1∗ , g2∗ ) = 0, β2 (g1∗ , g2∗ ) = 0. For all positive N , we nd two physically
85
equivalent xed points with purely imaginary coupling constants, hence all operator
dimensions remain real. These xed points are IR stable for all N (the stability
∂βi
matrix Mij = ∂gj
has positive eigenvalues). Note that this is dierent from the O(N )
versions of these models [28, 29], where one nds a critical N , whose one-loop value is
' 1038, below which the IR stable xed points with real coupling constants disappear.
However, to all orders in the 1/N expansion, the xed point couplings and conformal
dimensions in the Sp(N ) models are related to the ones in the O(N ) models by the
replacement N → −N .
Figure 4.1 shows the RG ow directions for N = 2. The arrows indicate how
the coupling constants ow towards the IR. The two IR xed points are physically
equivalent because they are related by gi → −gi . At higher values of N , the qualitative
behavior of the RG ows and xed points remain the same. We still have a UV
A special structure emerges for N = 2. In this case we nd the xed point solution
r
(4π)3 67
g2∗ = 2g1∗ , g1∗ =i 1+ 2
+ O( ) , (4.7)
5 180
and the conformal dimensions of the fundamental elds are equal (the three-loop term
8 7 2 269 − 702ζ(3) 3
∆σ = ∆χ = 2 − − − + ... . (4.8)
15 450 33750
N. For instance, we can explicitly check for which N the dimensions of σ and χ are
equal. A direct calculation using the beta functions and anomalous dimensions up
to three-loops shows that this only happens for N = 2 and N = −1. The latter
86
6
y
IR
UV
IR
-2
-4
-6
x
-3 -2 -1 0 1 2 3
Figure 4.1: The zeroes of the one-loop β functions and the RG ow directions for
q q
3 (4π)3
the OSp(1|2) model. The coordinates are dened via g1 = i (4π)
5
x , g2 = i
5
y,
and the red dots correspond to the stable IR xed points.
case corresponds to the 3-state Potts model xed point of the the theory with two
from the results in [28, 29]. Let us quote the explicit result for the quadratic operators
arising from the mixture of σ2 and Ωij χi χj . These operators have the same classical
7 For N = −1, there is also a (non-unitary) solution with g1∗ = g2∗ which corresponds to two
decoupled Fisher models [37], and hence the dimensions of the fundamental elds are trivially equal.
87
we nd the mixing matrix
g 2 (N +4) g12 (g12 (21N −134)−6g1 g2 (2N +5)+5g22 ) N g1 (6g1 −g2 ) N g12 (18g12 +161g1 g2 +15g22 )
− 16(4π)3 + 108(4π)6 6(4π)3
+ 108(4π)6
g1 (6g13 (11N +20)−g12 g2 (11N −324)+12g1 g22 −13g23 )
g1 (g2 −6g1 ) N g12 +4g22 N g12 (160g12 +36g1 g2 −41g22 )−181g24
6(4π)3
− 216(4π)6
− 6(4π) 3 + 216(4π)6
(4.9)
1
where index ` ' corresponds to the operator Ωij χi χj , 2
and index ` ' corresponds to σ2.
The diagrams contributing to the calculation are listed in Fig. 4.2.
Figure 4.2: Diagrams contributing to the mixing of σ2 and Ωij χi χj operators to two-
loop order. Notice that in these diagrams, we have not distinguished the σ eld and
χi elds. Therefore, each one of these diagrams represent several diagrams where the
elds are dierent. The solid dots represent one-loop counter terms for either the
propagator or the vertex, and the crossed dots represent one-loop counter terms for
The two eigenvectors of γ ab give the two linear combinations of σ2 and Ωij χi χj ,
and the eigenvalues γ± give their anomalous dimensions, so that ∆± = d − 2 + γ± .
We nd that one of these combinations is a conformal primary, and the other one
−40 + 104
3
2 6800 − 34190
3
2
∆− = ∆σ + 2 = 4 + + + ... ,
N N2 (4.10)
100 − 395
3
2 −49760 + 237476
3
2
∆+ = 4 + + + ... .
N N2
the available large N results for the critical exponent ω in the quartic O(N ) theory
8
[10, 79].
8 This critical exponent is related to the derivative of the beta function in the O(N ) theory with
quartic interaction; see [29] for more details on the comparison to the large N results.
88
Figure 4.3: Dierent Padé approximations of ∆χ for the N =2 model.
A0 + A1 d + A2 d2 + . . . + Am dm
Padém,n [d] = . (4.11)
1 + B1 d + B2 d2 + . . . + Bn dn
For N = 2, by demanding that the expansion have the same behavior as (4.8), we
can x four coecients in the Padé approximant. Here we use Padé1,2 , and Padé2,1 ,
3
As expected, this is below the unitarity bound in d = 5, which is ∆= 2
. The plots
of dierent Padé approximants are shown in Fig. 4.3. For the next primary operator,
89
whose scaling dimension is given in (4.30), using the Padé2,1 approximant we estimate
cubic theory (4.3) is equivalent to the UV xed point of the quartic theory (4.1).
9 If
so, then not only do we know the 6 − expansion of ∆χ to O(3 ), we can also use the
(4.13)
Together, these expansions allow us to determine ten coecients of the Padé ap-
proximant. Disregarding the approximants that have a pole, we obtain the estimate
∆N
χ
=4
≈ 1.478 in d = 5.
d = 6 − . The leading order term in the expansion for the corresponding O(N )
models was computed in [78]. Sending N → −N , we nd the following result for
F̃ = sin( πd
2
) log ZS d in the cubic Sp(N ) models
π (g2∗ )2 − 3N (g1∗ )2
F̃IR = F̃UV − 3
+ O(3 ), (4.14)
17280 (4π)
where F̃UV = (1 − N )F̃s , and F̃s the value corresponding to a free conformal scalar.
Plugging in the explicit solutions for the xed point couplings g1∗ , g2∗ , it is straight-
9 This equivalence holds for Sp(N ) models with suciently large N , but it is not completely clear
if it applies for N = 4.
90
forward to verify that, for all N ≥ 2, we have
π 2
F̃IR = −F̃s − + O(3 ) . (4.16)
43200
Note that for the Fisher model [37] of the Lee-Yang edge singularity, corresponding to
N = 0 and imaginary g2∗ , the inequality (4.15) does not hold. For general non-unitary
theories the inequality does not have to hold; remarkably, it does hold for the Sp(N )
models with positive N.
Similarly, using the results of [78] for the quartic O(N ) theory in d = 4 − , we can
compute the sphere free energy at the IR xed point of the model (4.1) in d = 4 − .
We nd
π N (N − 2) 3
F̃quartic = F̃free − + O(4 ) , (4.17)
576 (N − 8)2
where F̃free = −N F̃s . We observe that F̃free > F̃quartic is satised for all N > 2. 10
The case N =8 is special and needs to be treated separately. Here the one-loop
15 3
βλ = −λ + λ + O(λ4 ). (4.18)
32π 4
32π 4
Therefore, at the IR xed point λ2∗ = 15
+ O(3/2 ), and we get
π 2
F̃quartic = F̃free − + .... (4.19)
360
10 For <0 the interacting theory has a UV xed point. So, in this case F̃free < F̃quartic , again in
agreement with the conjectured F̃ theorem.
91
4.3 Symmetry enhancement for N = 2
Let us write the cubic Sp(2) model in terms of a real scalar σ and a single complex
anti-commuting fermion θ:
Z
d µ 1 2 1 3
S= d x ∂µ θ∂ θ̄ + (∂µ σ) + g1 σθθ̄ + g2 σ . (4.20)
2 6
For g2 = 2g1 , i.e. the xed point relation (4.7), this action possesses a fermionic
seen explicitly in eq. (4.8). This complex fermionic symmetry enhances the Sp(2)
to OSp(1|2), which is the smallest supergroup. The full set of supergroup generators
+ 1 ∂ ∂ − 1 ∂ ∂
Q = σ +θ , Q = σ − θ̄ ,
2 ∂ θ̄ ∂σ 2 ∂θ ∂σ
(4.22)
∂ ∂ 1 ∂ ∂
J+ = θ , J − = θ̄ , 3
J = θ − θ̄ ,
∂ θ̄ ∂θ 2 ∂θ ∂ θ̄
and it is not hard to check that they satisfy the algebra of OSp(1|2):
3 ± + −
J , J = ±J ± , J , J = 2J 3 (4.23)
3 ± 1 ± ∓
J , Q = ± Q± , J , Q = −Q± (4.24)
2
1 1
{Q± , Q± } = ± J ± , {Q+ , Q− } = J 3 (4.25)
2 2
example, let us study the mixing of the Sp(2) singlet operators σ2 and θθ̄. Setting
92
N =2 in (4.9), we nd the scaling dimensions of the two eigenstates:
2 1 8 7 2
∆+ = 4 − + 2 + O(3 ) , ∆− = 4 − − + O(3 ). (4.26)
3 30 15 450
O+ = σ 2 + 2θθ̄, (4.27)
which is invariant under all the OSp(1|2) generators. The second corresponds to
6. In [99] it was argued that the q→0 limit of the q -state Potts model is described
by the OSp(1|2) sigma model [101, 102]. There is evidence that the upper critical
dimension of the spanning forest model is 6 [100], and the 6− expansions of the
1 7 2 269 − 702ζ(3) 3
η=− − − + O(4 ) , (4.28)
15 225 16875
1 1 173 − 864ζ(3) 3
ν −1 = 2 − − 2 − + O(4 ) . (4.29)
3 30 27000
Remarkably, the operator dimensions we have calculated (4.8) and (4.26) agree with
d η
∆χ = −1+ , ∆+ = d − ν −1 . (4.30)
2 2
93
Similarly, we can match the 6− expansions of the sphere free energies. For the
π(q − 1)(q − 2) 2
F̃q = (q − 1)F̃s + + O(3 ) , (4.31)
8640(3q − 10)
and for q = 0 this matches the F̃ of our OSp(1|2) model, (4.16). These results provide
strong evidence that the OSp(1|2) symmetric IR xed point of the cubic theory (4.20)
describes the q→0 limit of the q -state Potts model.
the q → 0 limit of the ferromagnetic q -state Potts model, indicate the existence
found in [100] are in good agreement with the Pade extrapolations of 6 − expansions
exhibited in section 2.1. This provides additional evidence for the existence of the
critical theories with OSp(1|2) symmetry. It would be of further interest to nd the
critical statistical models that are described by the Sp(N ) invariant theories with
N > 2.
94
Chapter 5
5.1 Preliminaries
In this chapter, we will show that the 0-state Potts model and the Sp(2) models
are indeed equivalent at the level of tensor structures. The N + 1-state Potts model
1 ijk i j k
gd φ φ φ . (5.1)
6
Where,
N
X +1
ijk
d = eiα ejα ekα . (5.2)
α=1
95
Where the eiα satises:
N
X +1
eiα = 0 , (5.3)
α=1
N
X +1
eiα ejα = (N + 1)δ ij . (5.4)
α=1
XN
eiα eiβ = (N + 1)δαβ − 1 . (5.5)
i=1
the O(N ) cubic model have been computed in [103]. For the O(N ) cubic model, the
Sp(2) model corresponds to setting N = −2, and g2 = 2g1 = 2g , the beta function
g 5 3 335 5 144889 53ζ(3)
βSp2 = − + g − g + + g7
2 4 72 5184 4
7502393 53π 4 5297ζ(3) 125ζ(5)
+ − + − − g9 , (5.7)
31104 144 18 3
1 2 25 4 14213 5ζ(3)
γφSp2 = − g + g + − + g6
6 108 7776 6
4
78667 8π 209ζ(3) 160ζ(5)
+ + + − g8 . (5.8)
5184 135 108 9
The 0-state Potts model corresponds to N = −1, but the coupling and other
96
get that the beta function and anomalous dimensions are:
g 5 2 3 335
βP otts = − − + O(∆) ∆ g + − + O(∆) ∆4 g 5
2 4 72
144889 53ζ(3)
− + + O(∆) ∆6 g 7
5184 4
7502393 53π 4 5297ζ(3) 125ζ(5)
+ − + − − + O(∆) ∆8 g 9 , (5.9)
31104 144 18 3
1 2 2 25 4 4 14213 5ζ(3)
γφP otts = + O(∆) ∆ g + + O(∆) ∆ g + − + O(∆) ∆6 g 6
6 108 7776 6
78667 8π 4 209ζ(3) 160ζ(5)
+ + + − + O(∆) ∆8 g 8 . (5.10)
5184 135 108 9
Notice that the terms are very similar to the Sp(2) model in the ∆ → 0 limit, if we
rescale g → g/∆. In fact, the only dierence is that the terms have alternating signs in
successive loop orders. Namely, one loop contributions have opposite signs, two loop
contributions have the same sign, three loop contributions have the opposite sign, etc.
one could easily verify. Therefore, they will give exactly the same answer for physical
This strongly suggests that the two models might be equivalent. We will show that
the tensor structures are indeed equal for general zero, two and three-point functions,
97
Figure 5.1: Diagram for what Gracey calls T5
X
T5 dijk = dil1 l2 djl3 l4 dkl5 l6 dl1 l3 l5 dl2 l4 l6 (5.11)
l1 ,...,l6
X X
= (eiα1 elα11 elα21 )(ejα2 elα32 elα42 ) . . . (elα25 elα45 elα65 ) . (5.12)
α1 ,...,α5 l1 ,...,l6
Now we perform the sum over the l's, using 5.5, for example:
X
elα11 elα14 = (N + 1)δα1 α4 − 1 . (5.13)
l1
X
T5 dijk = eiα1 ejα2 ekα3 [(N + 1)δα1 α4 − 1][(N + 1)δα1 α5 − 1][(N + 1)δα2 α4 − 1]
α1 ,...,α5
(5.14)
If we multiply the 6 factors out, we will get 64 terms, carrying various powers of
(N + 1) and −1. Each of these terms can be thought of as a coloring of the original
6 propagators in the graph. Each solid line corresponds to a factor of (N + 1), while
98
each dashed line corresponds to a factor −1. We will have one diagram where all 6
lines are solid, 6 diagrams where 5 lines are solid..., until one diagram where all lines
are dashed. Notice also that from 5.3, a three-point diagram will only be non-zero
if all its vertices are connected by solid lines, since an isolated vertex will give zero.
Also, when there are free clusters of isolated solid lines in the middle of the diagram,
each would contribute another factor of (N + 1) from having a free index in the sum.
Likewise, a vertex with three dashed lines can be considered as a trivial cluster, and
also contribute a factor of (N + 1) due to having a free index under the sum.
longer be connected.
This has a natural interpretation. It can be easily argued using induction [38],
that all terms contributing to the lowest power of (N + 1) have no solid line cycles,
Thus, in order to compute the tensor structures in the 0-state Potts model, we
need to enumerate the forests, while remembering that each one will carry +1 or −1
depending on how many dashed lines there are.
99
Our proof will be roughly as follows. The ma jority of the work will be to rst
show that all two-point tensor structures are the same in the Potts model and the
Sp(2) model. Then we can use this to deduce that it holds for three-point functions
only be non-zero if there's a solid line going through the middle and connecting the
two outer points. There must be one and only one such line for each forest, otherwise
there will be a cycle. For each forest, we call this line the principal line. There can
be trees that are connected to the principal line, which we call rooted trees, and
there can be trees not connected to the principal line, which we call isolated trees.
In the Sp(N ) model, we only consider the two-point function of θθ̄. For the same
topology, the tensor value for σσ must be the same due to Osp(1|2) invariance that
was proved in chapter 4. To motivate the analogy to the Potts model case, we will
denote lines of θ by solids, and σ by dashes. It can be easily reasoned that because all
connecting the two outer points. If this is not the case, then a solid line must either
isolated trees in a valid Feynman diagram. However, there can be other cycles not
connected to the principal line. Further. the cycles do not share any edges.
5.2 shows an example of possible structures of a given topology in the Potts case
100
Figure 5.2: Example illustrating various terminologies.
In the proof we will also need the following fact. We start with a given topology
and a principal line. Suppose there are P vertices on the principal line. V vertices
outside of the principal line, and E edges not on the principal line. We have the
following relation:
2E = 3V + P . (5.16)
This can be reasoned as follows. Suppose we rst add in the vertices without any
edges (except those on the principal line). Since the interaction is cubic, every single
vertex must have 3 edges. Each of the P vertex on the principal line already has
2 edges, so they need 1 more edge. Each of the additional V vertices not on the
principal line needs 3 edges each. Whenever we add an edge not on the principal line,
Therefore, we have:
101
We will prove a stronger statement than the original claim. We will show that for
a given topology and an arbitrary principal line, the tensor factors calculated from the
0-state Potts model agree with the Sp(2) model, with a factor of (−1)loop order . Since
this equivalence holds for any principal line, the equivalence extends to all tensor
Given a topology and a xed principal line, with V vertices and E edges outside of
the principle line, where each edge can be either dashed or solid, we dene three sets
of graphs.
1. The set of cycle graphs, {GC }, consisting of all graphs on the topology
where there can be zero or more solid line cycles. These cycles cannot intersect
with the principal line, or themselves. These clearly have a 1-to-1 correspondence
with Feynman diagrams in the Sp(2) model. Each vertex is either dash-dash-dash or
contained in n solid line cycles (hence V − L vertices not contained in them) we dene
its value to be V (gc ) = 2V −L × (−2)n . This corresponds to the value of the Feynman
diagram in the Sp(2) model. Since there will be V −L dash-dash-dash vertices, each
GC equals to the sum of all Feynman diagrams in the Sp(2) model with the given
2. The set of forest graphs, {GF }, consisting of all graphs on the topology where
there are no solid lines cycles, including those formed with the principal line. These
are just the set of all forest diagrams on the given topology, which have a 1-to-1
correspondence with the diagrammatic method of calculating the tensor factor of the
102
0-state Potts model. for each individual graph gF ∈ GF , suppose it has k solid edges
outside of the principal line (hence it has E −k dashed lines) we dene its value to be
(−1)E−k . Therefore the sum of the value of all graphs in GF gives exactly the tensor
factor in the 0-state Potts model with the given topology and principal line.
Therefore our goal is to show the sum of value of all graphs in GC is equal to
GF . To do so we need to dene the following quantity which, unlike the other two
either theory.
3. The set of arrow graphs, {GA }, consists of all graphs on the topology where
we draw arrows on each edge with the following rules. Imagine initially all edges not
on the principal line are dashes. For each of the V vertex not on the principle line, we
make one of four choices: to color one of the three lines connected to this vertex to be
solids (each of these choice carry a factor of 1, we say that such a vertex is owing),
or not coloring any of the lines and leaving them as dashed (this carries a factor of
−1, we call such a vertex empty). The value of an individual graph gA ∈ GA is just
the product of all the value on its vertices, which is V (gA ) = (−1)# of empty vertices .
Let's look at the properties of {GA }. We can make each of the 4 choices indepen-
dently for each vertex, getting a total of 4V possibilities. If we sum the value of all
such 4V diagrams, we will get (3 − 1)V = 2V , since the choice made on each vertex is
independent.
Notice that we allow a line to be colored twice by its two adjacent vertices. We
will keep track of which vertex colored which line with arrows. If a line connecting
vertices A and B is colored by A, then an arrow is drawn on the line pointing away
from A, likewise for B. If a line is colored by both vertices, we draw a double arrow.
This is illustrated in 5.3. Since each vertex can color up to 1 of its adjacent lines,
there can never be more than one arrow coming out of any vertex.
103
Figure 5.3: Line 1 is colored by vertex A, line 2 is double-colored by vertices B and
1. It is impossible to create a cycle with part of the principal line. Since in order
to form a cycle with the principle line, you would need to color (k + 1) lines with only
k vertices, which is impossible. See 5.4 (a).
must be a simple cycle. For example, if you want to form a double cycle, with k
vertices, you need at least (k+1) colored lines. See 5.4 (b). Thus, all cycles created
must be disjoint.
Our main proof proceeds as follows. First we will show that the sum of all diagrams
in GC corresponds to the sum of all diagrams in GA that does not result in any cycles,
104
Then we will show that sum of all the forest diagrams in GA equals to the sum of
diagrams in GF times the factor (−1)E−V , which we've shown to be (−1)number of loops .
This explains the observation we made in the beginning of this chapter that the
dashes, there is only one such diagram. That means all of the V vertices outside the
principal line are dash-dash-dash, with each contributing a factor of 2. This diagram
has value 2V . On the other hand, as we've already noted, the sum of all diagrams in
GA is exactly 2V .
Conclusion 1: the only diagram with no cycle in GC is equal to the sum
of all diagrams in GA
Now we look at a diagram in GC where there is exactly one solid line cycle.
Suppose the cycle has length L. This diagram has V −L triple dash vertices, each
of 1. Finally, it has a solid line cycle, which gives an additional factor of N = −2. In
care about the choice of other vertices (in particular, they can form additional cycles).
a loop must all be owing. Further they either all ow clockwise or counterclockwise,
giving 2 possibilities, and each owing vertex gives a factor of 1. In total the sum of
105
Conclusion 2: the sum of diagrams with exactly one cycle in GC is equal
to the negative of the sum of all diagrams with at least one cycle in GA
Now we look at a diagram in GC with exactly 2 solid line cycles. As we argued
before the two cycles must be disjoint. Suppose the length of the two cycles are L1 and
L2 . Then the value of this diagram is the product of: 2V −L1 −L2 (for the triple-dash
of the choices of the other V − L1 − L2 vertices. However, each of these cycles have 2
factor of 4.
Conclusion 3: the sum of diagrams with exactly two cycles in GC is
X
GC = GC,0 + GC,1 + GC2 + . . . . (5.18)
all diagrams
Then let GA,n+ be the sum of all diagrams in GA with at least n cycles, we have:
X
GC = GC,0 + GC,1 + GC,2 + GC,3 . . . , (5.19)
all diagrams
= GA,0+ − GA,1+ + GA,2+ − GA,3+ + . . . . (5.20)
In other words, we have the sum of: all diagrams in GA , minus all diagrams with
at least one cycle in GA , plus all diagrams with at least two cycles in GA , etc. By
the principle of inclusion and exclusion, the sum is exactly equal to the sum of all
106
Figure 5.5: There is only one possible coloring for a rooted tree: every arrow points
Therefore, we nally conclude that, the sum of all diagrams in GC is equal to the
the topology. Therefore if the sum of all diagrams in GF is also equal to the sum of
Case 1: a single rooted tree If the forest diagram contains a single rooted
tree with k vertices and nothing else, there are k vertices which need to color exactly
k lines. All k vertices must be owing and the direction of the arrow is xed, as
illustrated in 5.5. Thus, there is exactly one conguration which colors each rooted
tree. Suppose there are V vertices and E edges in total. in GA , k of the vertices are
107
Figure 5.6: There are (k+1) choices for the empty vertex, and k choice for the double-
Case 2: a single isolated tree For a single isolated tree in GA , there are k+1
vertices which need to color k lines. This means that either one of the vertices doesn't
color any lines, or one of the lines is double colored. There are k + 1 choices for where
the empty vertex is, and k choices for where the double-colored line is. However,
there is a sign dierence due to the fact that if a vertex doesn't color any lines, it carries
a factor −1. Thus, the total factor is (−1)V −k (k + 1) + (−1)V −k−1 k = (−1)V −k . This
is illustrated in 5.6 Notice that the same graph in GF will carry a factor of (−1)E−k .
Again, in order to match the Potts model, we need an overall factor of (−1)E−V . This
Case 3: multiple rooted tree and isolated trees This is a simple generaliza-
tion of the previous two cases, indeed the GA value is again consistent with the GF
value, up to the same overall factor.
Combining the results from the previous section, we see that, for a given topology
108
X X
GF = (−1)#of loops GA (5.21)
all diagrams
This proves our claim, that for a given two-point function topology and principal
line, the value of the tensor structure in the Sp(2) model is equal to that in the 0-
state Potts model. If we sum over all possible principal lines, the equality still holds.
Therefore all two-point function tensor structures in the Sp(2) model are equal to
point function T3 dijk in the Potts model. We would like to show that it's equal to a
three-point function in the Sp(2) model. Again, there are two choices of three-point
functions in the Sp(2) case: σ 3 or σθθ̄, we can choose either one due to supersymmetry
(except σ3 is bigger by a factor of 2), so for convenience, we choose the latter.
We can convert the 3-point function to a 2-point function by merging two of its
lines into one (In the Sp(2) case, we must merge the σ line with one of the two θ
lines). This eectively creates a new two-point function tensor structure T2new related
to T3 via:
X
T2new = T3 dijk dijl = T3 T2 δ il , (5.23)
j,k
true.
109
The generalization to zero point functions is also simple. Notice that relation 5.20
still holds when we ignore the existence of the principal line. The consideration for
diagrams in GF is also simpler because there will only be isolated trees, and no rooted
Notice that the equality between all zero point functions in the two models imply
110
Chapter 6
Supersymmetry
Igor Klebanov, and Grigory Tarnopolsky. We also thank David Gross, David Poland,
teractions date back to Fermi's theory of beta decay. The rst application to strong
interactions was the seminal Nambu and Jona-Lasinio (NJL) model. In the original
paper [104] they considered the model in 3 + 1 dimensions with a single 4-component
Dirac fermion and Lagrangian
g
LNJL = ψ̄6 ∂ψ + (ψ̄ψ)2 − (ψ̄γ5 ψ)2 . (6.1)
2
In addition to the U (1) symmetry ψ → eiβ ψ , this Lagrangian possesses the U (1)
chiral symmetry under ψ → eiαγ5 ψ . Using the gap equation it was shown that the
chiral U (1) can be broken spontaneously, giving rise to the massless Nambu-Goldstone
111
boson. This was the discovery of the crucial role of chiral symmetry breaking in the
g
LNJL = ψ̄j 6 ∂ψ j + (ψ̄j ψ j )2 − (ψ̄j γ5 ψ j )2 , (6.2)
2
When considered in 2<d<4 this model gives rise to an interacting conformal eld
theory which describes the second-order phase transition separating the phases where
In addition to studying the NJL model with the U (1) chiral symmetry, we will
present new results for the Gross-Neveu (GN) model [25], which has a simpler quartic
interaction
g
LGN = ψ̄j 6 ∂ψ j + (ψ̄j ψ j )2 . (6.3)
2
Instead of the continuous chiral symmetry it possesses the discrete chiral symmetry
theory which describes the second-order phase transition where the discrete chiral
symmetry is broken.
While they are renormalizable in the sense of the 1/N expansion [107], at nite N it
1 It is also often called the chiral Gross-Neveu model [25]; in d = 2 it is equivalent to the SU (2Nf )
Thirring model [105, 106].
112
is important to know the UV completion of these theories. In [26, 27] it was suggested
that the UV completion in 2 < d < 4 is provided by the appropriate Yukawa theories.
The UV completion of the GN model (6.3) contains a real scalar eld σ and Nf
4-component Dirac fermions ψj , j = 1, . . . Nf :
1 1
LGNY = (∂µ σ)2 + ψ̄j 6 ∂ψ j + g1 σ ψ̄j ψ j + g2 σ 4 . (6.4)
2 24
section 6.2.
The UV completion of the U (1) symmetric NJL model (6.2), which contains a
1 1 1
LNJLY = (∂µ φ1 )2 + (∂µ φ2 )2 + ψ̄j ∂/ψ j + g1 ψ̄j (φ1 + iγ5 φ2 )ψ j + g2 (φφ̄)2 . (6.5)
2 2 24
This theory, which we will call the Nambu-Jona-Lasinio-Yukawa model (NJLY), will
There is a large body of literature on the GN and NJL CFTs in d=3 and their
applications; see, for example, [109, 26, 27, 110, 111, 112, 113, 114, 115, 3, 116,
117, 118, 119]. We will carry out further studies of these CFTs using the 4− and
dimensions of some low-lying operators, we will calculate the sphere free energy F.
The latter determines the universal entanglement entropy across a circle [66], and is
2 In d = 3, one may express the Lagrangian (6.2) in terms of 2N 2-component Dirac spinors
f
χi , χi+Nf by writing ψ i = (χi , χi+Nf ), i = 1, . . . , Nf . See for instance [108] for the explicit relation
i±
between 4-component and 2-component notations in 3d. The 2-component spinors χ
√ = (χi ±
χi+Nf )/ 2 have charge ±1 under the U (1) symmetry (6.6).
113
the quantity that enters the F -theorem [120, 54, 50, 67]. We will also discuss CT ,
the normalization of the correlation function of two stress-energy tensors. For the
GN model, the 1/N and expansions of CT were studied in [119]; in this chapter we
extend these results to the NJL model and also provide the numerical estimates in
(which is equivalent to one Dirac fermion in d = 3), the NJLY model is expected to
d<4 this theory denes a CFT with emergent supersymmetry" [121, 122], in the
sense that the RG ow drives the interactions to a supersymmetric IR stable xed
point, where the global U (1) symmetry becomes the U (1)R symmetry (see gure 6.1).
We will provide additional evidence for this using the 4− expansion of the NJLY
model with Nf = 1/2 to two loops, both for certain scaling dimensions and for the
Even more intriguingly, when the GNY model is continued to Nf = 1/4, which
to a CFT with 2 supercharges [121, 124, 125, 118, 126]. We will show that the O(2 )
corrections to scaling dimensions of operators continue to respect this emergent super-
symmetry.
3 This provides new support for the existence of an N = 1 supersymmetric
CFT in d = 3. Our Padé extrapolations of operator dimensions including the 2 cor-
rections are in good agreement with the conformal bootstrap approach to the N =1
supersymmetric CFT in d = 3 [118]. We also estimate CT and F for this theory.
3 This is reminiscent of the symmetry enhancement from Sp(2) to the supergroup OSp(1|2) at
the IR stable xed point of a cubic theory in 6− dimensions in chapter 4.
114
GNY model with N=1
8
●
6 =1 SUSY
◆
4 Ising
g2
2
0 ■
UV
-2 ?
▼
-4
-1 0 1 2 3 4 5
2
(g1 )
0 ■
UV
-5 ?
▼
-1 0 1 2 3 4 5
2
(g1 )
Figure 6.1: RG ow and xed point structure for the GNY model with N =1 (one
Majorana fermion in d = 3) and the NJLY model with N =2 (one Dirac fermion in
d = 3), obtained from the one-loop β -functions (6.7) and (6.68) in d = 4 − . The
respectively. The red triangles denote unstable xed points with negative quartic
potential which can be seen in the one-loop analysis in d = 4 − ; their fate in d=3
is unclear.
115
6.2 The Gross-Neveu-Yukawa Model
The β -functions for the GNY model with action (6.4) in d = 4 − , up to two-loop
1 2 2 4
1 6 4 2 2 17g23
βg2 = −g2 + 3g2 + 2N g1 g2 − 12N g1 + 96N g1 + 7N g1 g2 − 3N g1 g2 − ,
(4π)2 (4π)4 3
N +6 3 1 3 5 3 g1 g22
βg1 = − g1 + g + − (4N + 3)g1 − 2g1 g2 + , (6.7)
2 2(4π)2 1 (4π)4 4 12
where N = Nf tr1 = 4Nf . The model possesses an IR stable xed point at the critical
√
(g1∗ )2 1 (N + 66) N 2 + 132N + 36 − N 2 + 516N + 882 2
2
= + 3
+ O(3 ) ,
(4π) N +6 108(N + 6)
∗
√
g2 2
−N + 6 + N + 132N + 36 1
2
= + √
(4π) 6(N + 6) 54(N + 6)3 N 2 + 132N + 36
× 3N 4 + 155N 3 + 2745N 2 − 2538N + 7344
√
3 2
− (3N − 43N − 1545N − 1224) N + 132N + 36 2 + O(3 ) .
2
(6.8)
16π 2
Of course, there is also a xed point g1∗ = 0, g2∗ = g2Ising = 3
+ O(2 ) which
and Nf free fermions. By looking at the derivative of the beta functions at the xed
points, one can verify that (6.8) is attractive for all Nf , so one can ow to it from the
Ising" xed point along a relevant direction. Let us mention that there is formally
√
also a third xed point obtained from (6.8) by changing the sign of N 2 + 132N + 36.
This xed point is unstable in d = 4− due to the negative quartic coupling, g2∗ < 0,
but its dimensional continuation may produce a CFT in d = 3.
4 We have reproduced them using the general two-loop results [127, 128, 129, 130, 131] for the
Yukawa theories, which are reviewed in the Appendix.
116
The scaling dimensions of σ and ψ are found to be [114]
1 2 1 1 2 5 4
∆σ = 1 − + Ng + g − Ng ,
2 2(4π)2 1 (4π)4 12 2 4 1
3− 1 1
∆ψ = + 2
g12 − 4
(3N + 1)g14 . (6.9)
2 2(4π) 8(4π)
√
3 52N 2 − 57N + 36 + (11N + 6) N 2 + 132N + 36 2
∆σ = 1 − + + O(3 ) ,
N +6 36(N + 6)3
√
3 N +5 −82N 2 + 3N + 720 + (N + 66) N 2 + 132N + 36 2
∆ψ = − + + O(3 ) .
2 2(N + 6) 216(N + 6)3
(6.10)
These dimensions agree with [112] after correcting some typos in eq. (11) of that
paper (in particular, the coecient 33 should be changed to 3). Our O() term in ∆ψ
corrects a typo in eq. (4.44) of [3].
16π 2
Setting g1∗ = 0, g2∗ = 3
+ O(2 ) in (6.9), we may also recover the result at the
2
Ising xed point ∆Ising
σ = 1− 2
+ 108
+ O(3 ). One can then see that the Yukawa
operator σ ψ̄ψ is relevant at this decoupled xed point, and can trigger a ow to the
g2 1
γσ2 = − (g 2 + N g2 g12 − 2N g14 ) + 2γσ . (6.11)
(4π) 2 (4π)4 2
5 Throughout the chapter, one can obtain the corresponding expansions at the unstable xed
point by changing the sign of the square root.
117
At the xed point (6.8) we nd
√
N 2 + 132N + 36 − N − 30 1
∆σ2 = d − 2 + γσ2 = 2 + + √
6(N + 6) 54(N + 6) N 2 + 132N + 36
3
√
× (3N 3 + 109N 2 + 510N + 684) N 2 + 132N + 36
− 3N 4 − 658N 3 − 333N 2 − 15174N + 4104 2 + O(3 ) . (6.12)
We have also calculated the one-loop anomalous dimensions of the operators ψ̄ψ and
σ3:
(N + 2)g12 3g2
∆ψ̄ψ = + 2∆ψ , ∆σ3 = + 3∆σ . (6.13)
(4π)2 (4π)2
At higher orders these operator will mix, and one has to nd the eigenvalues of their
√
3 N 2 + 132N + 36 − N − 12
∆ψ̄ψ =3− + O(2 ) , ∆σ3 =3+ + O(2 ) .
N +6 2(N + 6)
(6.14)
model, which was discussed in [119]; the diagrams contributing to the term ∼ g12 are
5N (g1∗ )2 d N 1 5N
CT = N CT,f + CT,s − 2
= 2 + − , (6.15)
12(4π) Sd 2 d − 1 12(N + 6)
118
where Sd = 2π d/2 /Γ(d/2), and we used the values of CT for free scalar and fermion
theories:
d d
CT,s = , CT,f = . (6.16)
(d − 1)Sd2 2Sd2
the relevant curvature couplings that are marginal in d=4− [132, 133]
Z
√ 1 d−2 g2,0 4
SGNY = dd x g ψ̄i 6 ∂ψ i + (∂µ σ)2 + Rσ 2 + g1,0 σ ψ̄i ψ i + σ
2 8(d − 1) 24
(6.17)
η0
+ Rσ 2 + a0 W 2 + b0 E + c0 R2 ,
2
where R is the scalar curvature, W 2 is the square of the Weyl tensor, and E the Euler
density
be xed order by order in perturbation theory. On a sphere, the Weyl square term
drops out, and to the order we work below we will only need the renormalization of
the Euler coupling b0 (the R2 term and the renormalization of conformal coupling
are expected to play a role at higher orders [132, 133, 130]). The corresponding beta
function can be extracted from the results of [130, 134], and we nd
11N/4 + 1 1 1 9 6 3 2 6 1 4 1 2 2
βb = b − − N g + N g1 + N g1 g2 − N g1 g2 + . . . ,
360(4π)2 (4π)8 36 32 1 8 4 96
(6.19)
119
where b is the renormalized coupling, and its the relation to the bare one b0 = µ− (b −
11N/4+1
360(4π)2
+. . .) can be inferred from the above beta function. The coupling independent
term is related to the a-anomaly of the free fermions and scalar.
The calculation of the sphere free energy now proceeds as in [135, 136]. Keeping
Z
1 2 √ √
F = N Ff + Fs − g1,0 dd xdd y gx gy hσ ψ̄ψ(x)σ ψ̄ψ(y)i0
2
Z
1 2 √ √
− 2
g2,0 dd xdd y gx gy hσ 4 (x)σ 4 (y)i0
2!(4!)
Z
1 4 √ √ √ √
− g1,0 dd xdd ydd zdd w gx gy gz gw hσ ψ̄ψ(x)σ ψ̄ψ(y)σ ψ̄ψ(z)σ ψ̄ψ(w)i0 + δFb ,
4!
(6.20)
R √
where δFb = b0 dd x gE is the contribution of the curvature term, and Ff , Fs are
the sphere free energies of free fermion and scalar, which can be found in [78]. Starting
1 γ µ (x − y)µ
hσ(x)σ(y)i = Cφ , hψi (x)ψ̄ j (y)i = δij Cψ , (6.21)
|x − y|d−2 |x − y|d
d d
where Cφ = Γ( d2 − 1)/(4π 2 ) and Cψ = Γ( d2 )/(2π 2 ), and then Weyl transforming to
Z Z √ √
d √ √
d 2 d d
gx gy 2 3d
d xd y gx gy hσ ψ̄ψ(x)σ ψ̄ψ(y)i0 = N Cφ Cψ d xd y = N C φ C ψ I2 ( − 2) ,
s(x, y)3d−4 2
Z Z √ √
d d √ √ 4 4 4 d d
gx gy
d xd y gx gy hσ (x)σ (y)i0 = 4!Cφ d xd y 4(d−2)
= 4!Cφ4 I2 (2d − 4) .
s(x, y)
(6.22)
120
Here I2 (∆) denotes the integrated 2-point function of an operator of dimension ∆,
which is given by [137, 138, 50]
Z √ √ 1−d d+ 21
dd xdd y gx gy 2(d−∆) 2 π Γ d2 − ∆
I2 (∆) = = (2R) . (6.23)
s(x, y)2∆ Γ d+1
2
Γ(d − ∆)
1
hσ ψ̄ψ(x)σ ψ̄ψ(y)σ ψ̄ψ(z)σ ψ̄ψ(w)i0 = 6N 2 Cφ2 Cψ4
s2d−2 2d−2 d−2 d−2
xy szw sxw syz
s2xw s2yz − s2xz s2yw + s2xy s2zw 2 1
− 3N Cφ2 Cψ4 × + ,
(sxy syz szw sxw )d (sxy szw )d−2 (sxz syw )d−2
(6.24)
where we used a shorthand notation for the chordal distance sxy ≡ s(x, y). The
integral of this 4-point function over the sphere cannot be calculated explicitly, but
one can nd it as a series in d = 4 − . For this we used the Mellin-Barnes approach,
which is described in [135, 136]. The result for the integral reads
Z
√ √ √ √
dd xdd ydd zdd w gx gy gz gw hσ ψ̄ψ(x)σ ψ̄ψ(y)σ ψ̄ψ(z)σ ψ̄ψ(w)i0 =
N (N + 6) N N − 6 + 6(N + 6)(3 + γ + log(4πR2 ))
=− − + O().
2(4π)4 12(4π)4
(6.25)
1 2 3 1 2 4
F = N Ff + Fs − g1,0 N Cφ Cψ2 I2 d − 2 − g C I2 (2d − 4)
2 2 2 · 4! 2,0 φ
1 4 N (N + 6) N N − 6 + 6(N + 6)(3 + γ + log(4πR2 ))
+ g1,0 + + O() + δFb .
4! 2(4π)4 12(4π)4
(6.26)
121
Now replacing the bare couplings with the renormalized ones
N + 6 g13 11N/4 + 1
−
g1,0 = µ 2 g1 + + ... , g2,0 = µ g2 + . . . , b0 = µ b− + ...
32π 2 360(4π)2
(6.27)
we nd that all pole cancels, and the free energy is a nite function of the renormalized
couplings g1 , g2 , b.6
As explained in [135], in order to calculate the free energy at the critical point we
should now tune all couplings, including b, to their xed point values. Using (6.8),
we get
N
F = N F f + Fs −
48(N + 6)
√
(N + 99N + 18) N 2 + 132N + 36 + (80N 2 + 2103N + 6381) N + 108 2
2
−
7776(N + 6)3
+ δFb + O(3 ) . (6.28)
From the curvature beta function (6.19), we nd at the critical point
11N/4 + 1 1 1 9 3 1 1 1
b∗ = + N (g1∗ )6 + N 2 (g1∗ )6 + N (g1∗ )4 g2∗ − N (g1∗ )2 (g2∗ )2
360(4π) (4π)8 36
2 32 8 4 96
(6.29)
R √
and, using (6.8) and dd x gE = 64π 2 +O(), we nd that the Euler term contributes
√ √
N (882 + 66 N 2 + 132N + 36 + N (516 − N + N 2 + 132N + 36)) 2
δFb = + O(3 ) .
15552(N + 6)3
(6.30)
6 Note that the coupling dependent part in the renormalization of b is necessary to cancel poles
coming from diagrams at the next order. However, we still have to carefully include the Euler term,
as in [135, 136], as it aects the free energy at the xed point to order 2 .
122
Substituting this into (6.28), and writing the result in terms of F̃ = − sin(πd/2)F ,
we nd
N π2 1
F̃ = N F̃f + F̃s − − 3
161N 3 + 3690N 2 + 11880N + 216
96(N + 6) 31104(N + 6)
√
+ N + 132N + 36 N + 132N + 36 π3 + O 4 .
2 2 (6.31)
6.2.2 2 + expansions
In this section we review the known results for operator dimensions at the UV xed
point of the Gross-Neveu model in d = 2 + , and then compute its sphere free energy
to order 3 .
The action for the Gross-Neveu model [25] in Euclidean space in terms of bare
Z Z
√ 1 √
SGN = − d i i 2
d x g ψ̄i 6 ∂ψ + g0 (ψ̄i ψ ) + b0 dd x gR , (6.32)
2
where i = 1, . . . 2Nf , and we have included the Euler term which is needed for the
be [139, 3, 140]
N − 2 2 N − 2 3 (N − 2)(N − 7) 4
β = g − g + 2
g + 3
g + O(g 5 ) . (6.33)
2π 4π 32π
Therefore, one can see that there is a perturbative UV xed point at a critical coupling
g∗ given by
7
2π 2π π(N + 1) 3
g∗ = + 2
2 + + O(4 ) . (6.34)
N −2 (N − 2) 2(N − 2)3
7 This corrects a typo in g∗ in [3] on page 59 (there it is denoted by uc ).
123
The scaling dimensions are found to be [140]
1 + N − 1 2 (N − 1)(N − 2) 3
∆ψ = + 2
g − 3
g + O(g 4 ) ,
2 8π 32π
N −1 N − 1 2 (N − 1)(2N − 3) 3
∆σ = 1 + − g+ g + g + O(g 4 ) , (6.35)
2π 8π 2 32π 3
1 1 N − 1 2 (N − 1)(N − 6) 3
∆ψ = + + − + O(4 ) , (6.36)
2 2 4(N − 2)2 8(N − 2)3
1 N − 1 2 N (N − 1) 3
∆σ = 1 − − + + O(4 ) , (6.37)
N −2 2(N − 2)2 4(N − 2)3
1 N −3 3
∆σ2 = d + β 0 (g∗ ) = 2 + 2 + 2
+ O(4 ) . (6.38)
N −2 2(N − 2)
Let us now turn to the calculation of the free energy on S 2+ . To order g4, we
have
Z
1 g0 2 1 g0 3 1 g0 4 √
F =N Ff − S2 − S3 − S4 + b0 dd x gR , (6.39)
2! 2 3! 2 4! 2
Z Y
n
dxi gxi hψ 4 (x1 )...ψ 4 (xn )iconn
√
Sn = 0 , (6.40)
i=1
with ψ 4 ≡ (ψ̄i ψ i )2 . Using the at space fermion propagator in (6.21), and then
124
S2 = 2N (N − 1)Cψ4 I2 (2d − 2) ,
where the integral I2 (∆) is given in (A.12), and I3 (∆) denotes the integrated 3-point
Z √ √ √ 3(1+d)
dd xdd ydd z gx gy gz 3(d−∆) 8π
2 Γ(d − 3∆
2
)
I3 (∆) = = R . (6.42)
[s(x, y)s(y, z)s(z, x)]∆ Γ(d)Γ( 1+d−∆2
)3
Z Y
4
√ 2(N 2 − 3N + 4) 2N
S4 = 24N (N − 1)Cψ8 dxi gxi +
i=1
(s13 s14 s23 s24 )2d−2 s2d−4 2d−4 d d d d
12 s34 s13 s14 s23 s24
d+4 d+4
2s413 s424 sd12 sd34 + 4s212 s234 sd+2
+ s12 s34 d+2 4 4 d d 2 2 d+2 d+2
13 s24 − 4s12 s34 s14 s23 − 4s14 s23 s12 s34
+
(s12 s213 s214 s223 s224 s34 )d
s2−2d 2d 2−2d 2d 2 2 2 2
12 s13 s34 s24 + s13 s24 + s14 s23 − s12 s34
2 2
− 4(N − 1) 3d−2 d d 3d−2
, (6.43)
s13 s14 s23 s24
where smn ≡ s(xm , xn ) is a chordal distance on a sphere. Using the methods described
4 7π 2
3
1−d d+1
8 4(2−d) 2 π 2 e 2 γ(2−d) 2 10
S4 =Cψ (2R) d+1 3 192N (N − 1) (N − 2) − + − 25 +
Γ( 2 )(π d/2 ) 2 12
3
11
− 2(N − 2) − + 1 + O() . (6.44)
5 5
N − 2 g2
− (N − 2)2 N − 2 3
g0 = µ g− + − g + ... , (6.45)
2π 4π 2 2 8π 2
125
we nd a surviving pole in F of order g 4 /. This pole can be cancelled provided we
N N (N − 1)(N − 2) g 4
b0 = µ b − + ,
48π 30(4π)5
N N (N − 1)(N − 2) 4
βb = −b + − g , (6.46)
48π 6(4π)5
where the coupling independent term is due to the trace anomaly of the free fermion
R √
eld, and we used dd x gR = 8π + O(). In order to obtain the correct expression
N (N − 1) 3
δFb = − 3
+ O(4 ) , (6.47)
48(N − 2)
N (N − 1) 2 N (N − 1)(N − 3) 3
F = N Ff + − + O(4 ) , (6.48)
24(N − 2)2 16(N − 2)3
This result agrees with the one obtained in [135] using a conformal perturbation
theory approach; this provides a non-trivial check on the procedure used here which
8 Note that even though we had to compute S4 to x the renormalization of b, we cannot obtain
F to order 4 . That would require xing βb to order g5 .
126
6.2.3 Large N expansions
In this section we test the 4− and 2+ expansions by comparing them with the
known 1/N expansions [109, 69, 141]. The general form of the large N expansions of
d−1 1 1
∆ψ = + γψ,1 + 2 γψ,2 + O(N −3 ) ,
2 N N
1 1
∆σ = 1 + γσ,1 + 2 γσ,2 + O(N −3 ) ,
N N
1 1
∆σ2 = 2 + γσ2 ,1 + 2 γσ2 ,2 + O(N −3 ) , (6.50)
N N
where
Γ(d − 1) 2 (d − 1)Ψ(d) 4d2 − 6d + 2
γψ,1 =− d , γψ,2 = 4γψ,1 + ,
Γ( 2 − 1)Γ(1 − d2 )Γ( d2 + 1)Γ( d2 ) d−2 (d − 2)2 d
(6.51)
d−1
γσ,1 = −4 γψ,1 ,
d −2
2 6d2 Θ(d) 16(d − 1)2 Ψ(d) 4(d − 1) (d4 − 2d3 − 12d2 + 20d − 8)
γσ,2 = −γψ,1 + −
d−2 (d − 2)2 (d − 2)3 d
(6.52)
and
127
In these equations we dened
d d
Ψ(d) = ψ(d − 1) − ψ(1) + ψ 2 − −ψ ,
2 2
0 d
Θ(d) = ψ − ψ 0 (1) ,
2
0 0 d 0 d
Φ(d) = ψ (d − 1) − ψ 2 − −ψ + ψ 0 (1) , (6.54)
2 2
For F̃ , theO(N 0 ) result can be obtained from the general formula [68, 78] for the
R 2
change in F under a double trace" deformation δS = g O∆ , where O∆ is a scalar
Z ∆− d2
1 d d
δ F̃ = du u sin(πu)Γ +u Γ −u . (6.58)
Γ (d + 1) 0 2 2
Z d
−1
1 2 d d
F̃ = N F̃f + du u sin(πu)Γ +u Γ − u + O(1/N ) . (6.59)
Γ (d + 1) 0 2 2
The 4− expansion of this agrees with the large N expansion of (6.31), i.e.
π2 π3 π2 π3 1 3π2 17π3 1
F̃ = N F̃f + F̃s − + + − − − + O(1/N 3 ),
96 192 16 32 N 8 32 N2
(6.60)
128
Similarly, the 2+ expansion of (6.59) agrees with the large N expansion of (6.49),
i.e.
π3 π4 1 π3 π4 1 π3 3π4
F̃ = N F̃f + − + − + 2 − + O(1/N 3 ) .
48 32 N 16 16 N 6 32
(6.61)
For CT , the relative O( N1 ) correction to the answer for free feermions is given in
[119]:
CT 1
CT = N CT,f 1 + + O(1/N 2 ) ,
N
Ψ(d) d−2
CT 1 = −4γψ,1 + . (6.62)
d + 2 (d − 1)d(d + 2)
Its expansion in d = 4− can be seen to match (6.15), in particular reproducing the
near d = 4 and d = 2 agrees with the available perturbative results. Clearly, the
degree" n+m of the approximant is bound by how many terms in the -expansion
are known. Such approximants may be derived for any nite N, and it is useful to
compare their large N behavior as a function of d to the 1/N -expansion results listed
in the previous section. When several approximants Pade[m,n] are possible for the
same quantity, we use such comparisons to large N results to choose the one which
129
For ∆ψ , ∆σ and ∆σ2 we know the 4− expansion to order 2 , and the 2+
expansion to order 3 . This allows to use Padé approximants with m + n = 6. For ∆ψ
and ∆σ , we nd that Pade[4,2] has no poles in 2<d<4 and for large N is in good
agreement with the results (6.51) and (6.52). For ∆σ2 , we perform the Padé on the
critical exponent ν −1 = d − ∆σ2 , and then translate to ∆σ2 at the end. In this case,
we nd that Pade[1,5] is the only approximant with no poles, and it matches well to
the large N result. The resulting estimates in d=3 for these scaling dimensions are
given in Table 6.1. In gure 6.2, we plot our 3d estimates for ∆ψ and ∆σ , compared to
the large N curve obtained from (6.57) by eliminating N to express ∆σ as a function
N 3 4 5 6 8 20 100
Table 6.1: Estimates of scaling dimensions and sphere free energy at the d = 3
interacting xed point of the GN model.
For the sphere free energy of the interacting CFT, F̃ , we nd it convenient to
which is essentially the interacting part of the free energy in the GN description, but
it includes the contribution of a free scalar from the GNY point of view. Using the
results (6.31) and (6.49), we can use Padé approximants with n + m = 8, and we
nd that Pade[4,4] has no poles and agrees well with the large N result (6.58). The
130
Δσ
1.0
0.9
◆
Large N
◆ SUSY
0.8 ◆
0.7
◆
◆
0.6
◆
N=1
xed point discussed in Section 6.4. The black dotted line is the SUSY relation
∆σ = ∆ψ − 1/2.
are given in Table 6.1 for a few values of N. In Figure 6.3, we also plot the result of
the constrained Padé approximants for F̃ − N F̃f as a function of 2<d<4 for a few
values of N, showing that they approach well the analytical large N formula (6.58).
˜ ˜
F - N Ff
0.030
0.025
0.020
0.015
N=4
N=6
0.010
N=20
N=100
0.005 Large N
Figure 6.3: Padé estimates of F̃ − N F̃f in 2<d<4 compared to the large N result
(6.58).
We can use our estimates to make some tests of the F -theorem. In the GNY
description, we can ow to the critical theory from the free UV xed point of N
fermions plus a scalar, while in the GN description one can ow to the critical theory
131
to the free fermions. Then, the F -theorem implies the inequalities
We veried that our estimates satisfy these inequalities for all values of N. As an
example, for N = 4 we get FGN /(N Ff + Fs ) ≈ 0.93 and FGN /(N Ff ) ≈ 1.06. In
the GNY description, we also see that we can ow to the critical GN point from the
decoupled product of the Ising CFT and N free fermions. This implies
Using the estimates for FIsing derived in [135], we have checked that this inequality
indeed holds. Using the Padé approximants as a function of 2 < d < 4, we can also
verify that both (6.64) and (6.65) are satised, in terms of F̃ , in the whole range of
cannot apply the strategy described above because the 2+ expansions (6.37) become
singular. Directly in d = 2, the GN model is equivalent to the Thirring model for a
single 2-component Dirac fermion which can be solved via bosonization and has a line
of xed points; the dimensions of ψ and σ = ψ̄ψ depend on the interaction strength.
Therefore, for these operators we can only perform the one-sided" Pade[1,1] on the
∆N
ψ
=2
≈ 1.076 , ∆N
σ
=2
≈ 0.656 , ∆N
σ2
=2
≈ 1.75 , (6.66)
132
and for the sphere free energy
1 N 3 1 8 3 1 2 5
β1 = − g1 + + 2 g1 + − g1 g2 + g1 g2 + (7 − 3N )g1 ,
2 (4π)2 2 (4π)4 3 9
1 10 2 2 4 1 20 3 6 4 10 2 2
β2 = −g2 + g + 2N g2 g1 − 12N g1 − g − 96N g1 − 2N g2 g1 + N g1 g2 .
(4π)2 3 2 (4π)4 3 2 3
(6.68)
The one-loop terms above agree with the Nb = 1 case of the results in [113]. Solving
√
(g1∗ )2 −N 2 + 448N − 1096 + (N + 76) N 2 + 152N + 16 2
2
= + 3
+ O(3 ) ,
(4π) N +4 100(N + 4)
∗
√
g2 2
3( N + 152N + 16 − N + 4) 9
2
= + √
(4π) 20(N + 4) 500(N + 4)3 N 2 + 152N + 16
× 3N 4 + 114N 3 + 764N 2 − 26192N + 1280
√
− (3N 3 − 114N 2 − 1932N − 320) N 2 + 152N + 16 2 + O(3 ) .
(6.69)
To get a positive solution for g2∗ we have picked the + sign for the square root.
The other choice of sign gives another xed point which is presumably unstable. In
addition, there is clearly a xed point with g1∗ = 0 and g2∗ corresponding to the O(2)
133
Wilson-Fisher xed point. The solution (6.69) yields an IR stable xed point for all
values of N.
The scaling dimensions of the elds are found to be
3− 1 2 1
∆ψ = + 2
g1 − 4
(3N + 2)g14 ,
2 (4π) 4(4π)
1 2 1 1 2 3 4
∆φ = 1 − + Ng + g − Ng . (6.70)
2 2(4π)2 1 (4π)4 9 2 2 1
√
3 N +2 −76N 2 + 98N − 1296 + (N + 76) N 2 + 152N + 16 2
∆ψ = − + 3
+ O(3 ) ,
2 2(N + 4) 100(N + 4)
2
√
2 56N − 498N + 16 + (19N + 4) N 2 + 152N + 16 2
∆φ = 1 − + + O(3 ) .
N +4 50(N + 4)3
(6.71)
The NJLY has two types of operators quadratic in the scalar elds: the U (1)
invariant operator φφ̄, and the charged operators φ2 and φ̄2 . The one-loop scaling
dimension of φφ̄ was determined in [113]. Using [131], we nd up to two loops
4 4 2 2
∆φφ̄ = d − 2 + g 2 − g + N g g
2 1 + 2γφ ,
3(4π)2 3(4π)4 2
2 2 4 2 2 4
∆φ2 =d−2+ g2 − g + N g g
2 1 − 6N g 1 + 2γφ . (6.72)
3(4π)2 3(4π)4 3 2
√
N 2 + 152N + 16 − N − 16 1
∆φφ̄ =2 + + √
5(N + 4) 250(N + 4)3 N 2 + 152N + 16
√
× (17N 3 + 104N 2 + 1252N + 1120) N 2 + 152N + 16
− 17N − 4646N + 2304N − 187712N + 4480 2 + O(3 )
4 3 2
(6.73)
134
and
√
N 2 + 152N + 16 − N − 36 1
∆φ2 =2 + + √
10(N + 4) 125(N + 4) N 2 + 152N + 16
3
√
× (3N 3 + 571N 2 + 688N + 240) N 2 + 152N + 16
− 3N 4 − 924N 3 + 7806N 2 − 47688N + 960 2 + O(3 ) . (6.74)
The 4− expansion of CT in the NJLY model proceeds similarly to that for the
GNY model presented in [119] and reviewed in section 6.2. In the NJLY model there
are two real scalar elds, and we need to replace Tσ by Tσ1 + Tσ2 . It is not hard
to check that each diagram contributing to the term ∼ g12 picks up a factor of 2
compared to the GNY model, since each of the internal scalar lines can be either φ1
or φ2 . Thus, we nd
5N (g1∗ )2 d N 2 5N
CT = N CT,f + 2CT,s − 2
= 2 + − . (6.75)
6(4π) Sd 2 d − 1 6(N + 4)
GNY model discussed earlier. The integrals are also nearly identical, except for
combinatorics factors due to the fact that we have two scalar elds, one of which has
iγ 5 coupling.
Z
1 2 √ √
F = N Ff + 2Fs − g1,0 dd xdd y gx gy hO1 (x)O1 (y)i0
2
Z
1 √ √
− 2
g2,0 dd xdd y gx gy hO2 (x)O2 i0
2
2!(4!)
Z
1 4 √ √ √ √
− g1,0 dd xdd ydd zdd w gx gy gz gw hO1 (x)O1 (y)O1 (z)O1 (w)i0 + δFb ,
4!
(6.76)
135
where we dened the operators O1 = ψ̄j (φ1 + iγ5 φ2 )ψ j and O2 = (φ21 + φ22 )2 , and
R
δFb = b0 dd xE is the Euler term. Evaluating the correlation functions above, we
nd
1 2 3 1 1 4
F = N Ff +2Fs − 2g1,0 N Cφ Cψ2 I2 d−2 + 2
2
64g2,0 Cφ4 I2 (2d−4)+ g1,0 I4 +δFb ,
2 2 2 · (4!) 4!
(6.77)
N + 4 g13
g1,0 = µ 2 g1 + + . . . , g2,0 = µ g2 + . . . (6.79)
32π 2
we nd that all poles cancel as they should. As explained in the GNY and GN
calculations, in order to obtain the correct expression for F at the IR xed point in
d = 4− we need to include the eect of the Euler term. Using an improved version
of the result from [130, 134], adapted to the presence of γ5 in the vertices, we get:
11N/4 + 2 1 1
βb = b − 2
− 8
N g12 (24g12 g2 − g22 + 9g14 (3N − 7)) . (6.80)
360(4π) (4π) 1296
From this we can solve for the xed point value b∗ with the couplings given in (6.69),
√ √
N (−1096 + 76 N 2 + 152N + 16 + N (448 − N + N 2 + 152N + 16)) 2
δFb = .
7200(N + 4)3
(6.81)
136
Putting this together with the integrals in (6.77), we nally nd in terms of F̃ :
N π2 1
F̃ = N F̃f + 2F̃s − − 3
149N 3 + 2372N 2 + 1312N + 64
48(N + 4) 14400(N + 4)
√
2
+ N + 152N + 16 N + 152N + 16 π3 + O(4 ) .
2 (6.82)
6.3.2 2 + expansions
In this section we consider the theory in 2 + dimensions with 4-fermion interactions
which respect the U (1) chiral symmetry. We begin with the action in Euclidean space
Z
1 1 1
SSPV = − d / i i 2 i 2 i 2
d x ψ̄i ∂ ψ + gS (ψ̄i ψ ) + gP (ψ̄i γ5 ψ ) + gV (ψ̄i γµ ψ ) , (6.83)
2 2 2
Gross-Neveu model [25] in d = 2. However, as we will see below, for our purposes it
is not consistent to set gV = 0 the corresponding operator respects the U (1) chiral
The one-loop beta-functions and anomalous dimension of the ψ eld were found
1 1
βS = gS − (N − 2)gS2 − gP gS − 2gV (gS + gP ) + . . . ,
π 2
1 1
βP = gP + (N − 2)gP2 − gP gS + 2gV (gS + gP ) + . . . , (6.84)
π 2
1
βV = gV + gP gS + . . . ,
π
137
and
1 1 N 2 2 2
1 2
∆ψ = + + g P + gS + 2gV − 4gV (gS − gP ) + (gS + gP ) .
2 2 4(2π)2 4(2π)2
(6.85)
We note that at the leading order in the 2+ expansion the evanescent operators do
not appear [106, 144, 145]. One of the UV xed points of (6.84) is
10
N ∗ 2π
gS∗ = −gP∗ = gV = , (6.86)
2 N
which corresponds to the SU (2Nf ) Thirring model [106, 146]. Indeed using the
1 i k 2 i k
relation for the SU (2Nf ) generators (T a )ij (T a )kl = (δ δ
2 l j
− δ δ ) and the Fierz
N j l
identity in d = 2 (γµ )αβ (γµ )γδ = δαδ δβγ − (γ5 )αδ (γ5 )βγ one nds
2
(ψ̄ψ)2 − (ψ̄γ5 ψ)2 + (ψ̄γµ ψ)2 = −2(ψ̄γµ T a ψ)2 . (6.87)
N
Z 1 2
SSU (2Nf )Thirring = − dd x ψ̄i ∂/ψ i + g (ψ̄i ψ i )2 − (ψ̄i γ5 ψ i )2 + (ψ̄i γµ ψ i )2 ,
2 N
(6.88)
N 2 1 1 N2 − 4 2
βg = g − g + O(g 3 ), ∆ψ = + + g + O(g 3 ) . (6.89)
2π 2 2 8π 2 N
It is plausible that this model is the continuation of the NJLY model (6.5) to d = 2+.
One nds that the UV xed point is g∗ = 2π/N , and the critical anomalous dimension
10 Equations (6.84) have four dierent non-trivial xed points. Two of them correspond to the GN
and SU (2Nf ) Thirring models.
138
reads
1 1 1 1 4 2
∆ψ = + + − 3 + O(3 ) . (6.90)
2 2 2 N N
We can also calculate the free energy of this model. To order g2, we have:
1 g0 2
F = N Ff − S2 + O(g03 ) , (6.92)
2! 2
Z
dx1 dx2 gx1 gx2 hO(x1 )O(x2 )iconn
√ √
S2 = 0 , (6.93)
where, under the Thirring description (6.87), O = −2(ψ̄γµ T a ψ)2 . Going through the
combinatorics, we nd
Evaluating this in d = 2+ and plugging in the xed point value (6.86) we nally
nd
(N 2 − 4)π3
F̃ = N F̃f + + O(4 ) . (6.95)
24N 2
139
6.3.3 Large N expansions
For the NJL model, the 1/N expansions of operator dimensions again assume the
2Γ(d − 1) 4 1
2
γψ,1 =− d , γψ,2 = 2γψ,1 Ψ(d) + + ,
Γ( 2 − 1)Γ(1 − d2 )Γ( d2 + 1)Γ( d2 ) d−2 d
(6.96)
4d2 (d2 − 5d + 7) 2
γφ,1 = −2γψ,1 , γφ,2 = −2γψ,2 + γψ,1 (6.97)
(d − 2)3
and
There is a considerable similarity between these results and the corresponding results
for the GN model. For example, γψ,1 in the NJL model is 2 times that in the GN
model. This is because the NJL lagrangian (6.2) contains two separate double-trace
operators, and hence there are two auxiliary elds, φ1 and φ2 , that couple to the
140
When (6.97) and (6.99) are expanded in d = 4− and d = 2+, they are consistent
with our -expansion results.
11 Setting d=3 in these equations gives [110, 147]
Turning to the sphere free energy, we can obtain the O(N 0 ) result for F̃ by simply
Z d
−1
2 2 d d
F̃ = N F̃f + du u sin(πu)Γ +u Γ − u + O(1/N ) . (6.101)
Γ (d + 1) 0 2 2
The 4− expansion of this agrees with the large N expansion of (6.82), i.e.
π2 π3 π2 π3 1 π2 17π3 1
F̃ = N F̃f + 2F̃s − + + − − − + O(1/N 3 ),
48 96 12 18 N 3 36 N2
(6.102)
and its 2+ expansion agrees with the large N expansion of (6.95), which yields
π3 π3
F̃ = N F̃f + − 2
+ O(1/N 4 ) . (6.103)
24 6N
For CT , the presense of the extra scalar eld compared to the GN case [119] again
poses no diculty. After some simple excercise commuting γ 5, we conclude that all
11 The form of (6.91) agrees with the large N result (6.99), but we haven't compared the coecient
of 2 .
141
due to the presence of two scalar elds. Hence, we nd for the NJL model
CT 1
2
CT = N CT,f 1 + + O(1/N ) ,
N
Ψ(d) d−2
CT 1 = −4γψ,1 + , (6.104)
d + 2 (d − 1)d(d + 2)
i.e. the correction CT 1 is just twice the corresponding term in the GN model (note
that γψ,1 in the NJL model is twice that of the GN model). This result is, in particular,
consistent with the fact that in the limit d → 4 we expect CT 1 to reproduce the
contribution of two free scalar elds. Its 4− expansion can be also seen to agree
with (6.75).
2+ expansions
12 derived in the previous sections to obtain rational approximants
of scaling dimensions and sphere free energy at the NJL xed point in 2 < d < 4,
for N >2 (the case N = 2, which displays the emergent supersymmetry, is treated
separately in Section 6.4). The results in d = 3 are given in Table 6.2, indicating which
approximants was chosen in each case. In gure 6.4, we also plot our 3d estimated for
∆ψ and ∆φ , compared to the large N results. From the expression for ∆φ in (6.100),
it appears that the expansion does not converge well already at O(1/N 2 ), therefore
we have just used (6.100) to order 1/N to produce the plot below.
For the sphere free energy F̃ , we again nd it convenient to perform the Padé
we plot the resulting Padé approximants as a function of 2<d<4 for a few values
of N, showing that they approach well the analytical large N formula (6.101).
142
N 4 6 8 10 12 20 100
Table 6.2: Estimates of scaling dimensions and sphere free energy at the d = 3
interacting xed point of the NJL model.
Δϕ
1.00
◆
◆
0.95
◆ Large N
◆
0.90 ◆
◆
0.85
◆
0.80
0.75
0.70
1.00 1.02 1.04 1.06 1.08
Δψ
Figure 6.4: Padé estimates in d=3 of ∆σ versus ∆ψ for N = 4, 6, 8, 10, 12, 20, 100,
compared to the large N results (6.100).
˜ ˜
F - N Ff
0.05
0.04
0.03
N=4
N=6
0.02
N=20
N=100
0.01 Large N
Figure 6.5: Padé estimates of F̃ − N F̃f in 2<d<4 compared to the large N result
(6.101).
We can again use our estimates to test the d = 3 F -theorem and its proposed
143
dimensions, the inequalities
hold, in accordance with the conjectured generalized F -theorem [78, 135]. Using the
results in [135] for the free energy of the O(2) Wilson-Fisher model, we have also
checked that
which is consistent with the fact that the theory can ow from the xed point con-
sisting of free fermions decoupled from the O(2) model, to the IR stable NJL xed
point.
in [78, 148, 149]. In d = 3, one nds the N =2 theory of a single chiral supereld (a
complex scalar and a 2-component Dirac fermion) with cubic superpotential, which
ows to a non-trivial CFT in the IR [150]. In d = 2, the model matches onto the (2, 2)
supersymmetric CFT with c = 1, which is the k=1 member of the superconformal
3k
discrete series with c= k+2
.
144
The component Lagrangian of the WZ model in 4d reads
1 1 1 λ λ2
LWZ = (∂µ φ1 )2 + (∂µ φ2 )2 + ψ̄ ∂/ψ + √ ψ̄(φ1 + iγ5 φ2 )ψ + (φ̄φ)2 , (6.107)
2 2 2 2 2 16
3 λ3 3 λ5 15 + 36ζ(3) λ7
βWZ = − λ + 2
− 4
+ 6
+ O(λ9 ) , (6.108)
2 2 (4π) 2 (4π) 8 (4π)
16π 2
from which one nds an IR xed point with λ2∗ = 3
+ O(2 ). The dimension of
d−1 d−1
∆φ = Rφ = . (6.109)
2 3
One also has the exact result ∆φ2 = ∆φ + 1, since the operator φ2 is obtained from φ
by acting twice with the supercharges (this is because, due to cubic superpotential,
one has the relation Φ2 = 0 in the chiral ring). It also follows from supersymmetry
[121, 123] that the dimension of the operator φφ̄ at the xed point is given by
1 1 + 12ζ(3) 3
∆φφ̄ = d − 2 + β 0 (λ∗ ) = 2 − 2 + + O(4 ) , (6.110)
3 18
which agrees with the explicit three-loop calculation of [123]. In performing Padé
extrapolation of this result to d = 3, we have found that Padé[1,2] and Padé[2,1] give
answers close to each other. Their average is ≈ 1.909, which is very close to the
value ≈ 1.91 reported using numerical bootstrap studies [148, 149]. We can also take
into account the fact that in d=2 the dimension of φφ̄ should approach 2. Since
145
ν −1 = d − ∆φφ̄ , and then return to ∆φφ̄ in d = 3, then we nd ≈ 1.94. This is
Nf = 1/2 and N =2 in the NJLY model [121, 122]. Note that the NJLY Lagrangian
(6.5), specialized to the case of a single Ma jorana fermion, coincides with the WZ
3
3g12 = g2 = λ2 . (6.111)
2
Indeed, setting N =2 in the result for the xed point couplings (6.8), we nd:
(g1∗ )2 1 1 2
2
= + + O(3 ) ,
(4π) 6 18
(6.112)
∗
g2 1 1 2
2
= + + O(3 ) .
(4π) 2 6
This is precisely consistent with the relation (6.111), and gives evidence of the emer-
gent supersymmetry in the N =2 NJLY model. Note that for this value of N, the
chiral U (1) symmetry of the NJLY model becomes the U (1) R-symmetry of the WZ
model.
14 Further evidence can be found by setting N =2 in the 4− expansions of
3
∆ψ = − , ∆φ = 1 − , ∆φ2 = 2 − , (6.113)
2 3 3 3
in agreement with the supersymmetry. In particular, the fact that the O(2 ) terms
vanish is consistent with the exact result (6.109), and we also see ∆φ2 = ∆φ + 1 as
146
discussed above. Furthermore, setting N = 2 in (6.73), we nd ∆φφ̄ = 2−2 /3+O(3 ),
in agreement with the result (6.110) in the WZ model.
It is also interesting to look at the sphere free energy. Setting N =2 in the 4−
expansion (6.82) of F̃ , we get
π2 π3
F̃N =2 = 2F̃s + 2F̃f − − + O(4 ) . (6.114)
144 162
This precisely agrees with the expansion of (5.23) in [78], which was derived using a
π 3
curvature term (6.81) contributes at O(3 ) order to F̃ : δ F̃b = − 1296 . This contribution
is crucial for agreement with [78]. Thus, (6.114) provides a nice perturbative test of
the exact formula for F̃ as a function of d proposed in [78] for the Wess-Zumino
model.
[135], the value of F for the O(2) Wilson-Fisher xed point in d=3 was estimated
1 3ζ(3)
to be FO(2) ≈ 0.124. Using also the value Ff = 8
log(2) + 16π 2
in d=3 [50], we see
that
[64, 154]:
√
CT 16(16π − 9 3)
= ≈ 0.7268 , (6.116)
CTUV 243π
where CTUV = 4CT,s is the value for the free UV theory of two scalars and one two-
component Dirac fermion. Let us compare this with a Padé extrapolation of the ratio
147
CT /CT,s using the boundary conditions
CT 1 in d = 2,
= (6.117)
CT,s
5 − 11
+ O( ) 2
in d = 4 − ,
6
where the 4− expansion was obtained by setting N =2 in (6.75). The Pade[1,1]
CT 49d − 76
= , (6.118)
CT,s d + 20
this result implies CTIR /CTUV ≈ 0.77. This is not far from the exact result (6.116),
demonstrating that the Padé approach works quite well. It would be useful to know
the next order in the 4− expansion, which may improve the agreement.
theory containing a 2-component Ma jorana fermion ψ [121, 124, 125, 118, 126]. This
theory must also contain a pseudoscalar operator σ , whose scaling dimension is related
to that of ψ by the supersymmetry,
1
∆σ = ∆ψ − . (6.119)
2
Some evidence for the existence of this N =1 supersymmetric CFT was found using
To describe the theory in d = 3, one can write down the Lagrangian [121, 124,
125, 118]
1 1 λ λ2
LN =1 = (∂µ σ)2 + ψ̄6 ∂ψ + σ ψ̄ψ + σ 4 . (6.120)
2 2 2 8
148
This model has N =1 supersymmetry in d = 3; the eld content can be packaged
in the real supereld Σ = σ + θ̄ψ + 12 θ̄θf , and the interactions follow from the cubic
superpotential W ∼ λΣ3 .15 It is natural to expect that this model ows to a non-
trivial N =1 SCFT in the IR. Note that the theory cannot be described as the UV
xed point of a lagrangian for a Ma jorana fermion with quartic interaction, because
The theory (6.120) is super-renormalizable in d = 3, and one may attempt its 4−
expansion [121]. To formulate a Yukawa theory in d = 4, one strictly speaking needs
√
simplication that occurs for this value is that N 2 + 132N + 36, which appears in
the 4− expansions (6.10), (6.12), equals 13 for N = 1. For this value of N, we nd
(g1∗ )2 1 3 2
2
= + + O(3 ) ,
(4π) 7 49
(6.121)
∗
g2 3 9 2
2
= + + O(3 ) .
(4π) 7 49
This is consistent with the exact relation 3g12 = g2 in the SUSY model. Indeed the
GNY Lagrangian (6.4), formally applied to the case of a single 2-component Majorana
fermion, coincides with (6.120) when 3g12 = g2 = 3λ2 . The result (6.121) gives a two-
loop evidence that the non-supersymmetric GNY model with N = 1 ows at low
us check that the supersymmetry relation (6.119) continues to hold at order 2 . Using
15 To obtain the component Lagrangian (6.120), one should eliminate the auxiliary eld f using
its equation of motion f ∼ λσ 2 .
149
(6.10) with N = 1, we indeed nd
3 1
∆σ = 1 − + 2 + O(3 ) , (6.122)
7 49
3 3 1
∆ψ = − + 2 + O(3 ) . (6.123)
2 7 49
1
∆σ2 = ∆σψ − 2
. Since σψ is a descendant, we also have ∆σψ = ∆ψ + 1. Thus, the
1
∆σ2 = ∆ψ + = ∆σ + 1 . (6.124)
2
3 1
∆σ2 = 2 − + 2 + O(3 ) , (6.125)
7 49
so that the supersymmetry relation (6.124) holds to order 2 . These non-trivial checks
provide strong evidence that the continuation of the GNY model to N =1 ows to a
8d − 11
∆σ = , (6.126)
25 − d
using the numerical bootstrap [118]. It is the value ∆σ ≈ 0.582 where the boundary
150
It is also important to know how the theory behaves when continued to d = 2.
It is plausible that the d = 2 theory has N = 1 superconformal symmetry, and
the obvious candidate is the tri-critical Ising model [124, 126], which is the simplest
supersymmetric minimal model [73, 155]. The Padé extrapolation (6.126) gives ∆σ ≈
0.217, which is quite close to the dimension 1/5 of the energy operator in the tri-
critical Ising model. This provides new evidence that the GNY model with N =1
extrapolates to the tri-critical Ising model in d = 2; in gure 6.6 we show how the
operator spectrum matches with the exact results in d = 2. Imposing the boundary
M4,5 N = 1 SCFT
∆
S S̄ = X( 3 , 3 ) ψ̄ψ, σ 3
2 2 3 ≈ 2.79
≈ 2.59 σψ, σ ψ̄
L−1 L̄−1 Φ
≈ 2.09 σ2
2
L−1 Ψ̄, L̄−1 Ψ
≈ 1.59 ψ, ψ̄
Φ̃( 3 , 3 )
5 5
≈ 1.09
1 σ
Ψ( 3 , 1 ) , Ψ̄( 1 , 3 )
5 10 10 5
≈ 0.59
Φ( 1 , 1 )
10 10
2 3 4 d
Figure 6.6: Qualitative picture of the interpolation of operator dimensions from the
From the β -functions (6.7) for the GNY model we may also deduce the dimensions
of two primary operators that are mixtures of σ4 and σ ψ̄ψ . They are determined by
151
the eigenvalues λ1 , λ2 of the matrix ∂βi /∂gj . At the xed point for N =1 we nd
3
∆1 = d + λ1 = 4 − 2 + O(3 ) ,
7
6 95
∆2 = d + λ2 = 4 + − 2 + O(3 ) , (6.127)
7 49
3
∆σ3 = ∆1 − 1 = 3 − 2 + O(3 ) . (6.128)
7
As a check, we may set N =1 in (6.14) and see that the order term in ∆σ3 indeed
vanishes. It is further possible to argue that this primary operator, when continued
to d = 2, matches onto the dimension 3 operator in the tri-critical Ising model, i.e.
the product of left and right supercurrents. This implies that ∆σ3 is not monotonic
setting =1 in (6.128) gives ∆σ3 ≈ 2.57, but this probably underestimates it. One
yield the estimate ∆σ3 ≈ 2.79 in d = 3. This slightly relevant operator is parity odd
Similarly, we can perform extrapolation of F̃ and check if the d=2 value is close
to πc/6 where c = 7/10 is the central charge of the tri-critical Ising model. Setting
N =1 in (6.31), we have
π 2 3 4
F̃ = F̃s + F̃f − + + O( ) . (6.129)
6 112 98
Performing a Padé approximation of the quantity f (d) = F̃ − F̃f , we nd that the
average of the standard Padé[2,1] and Padé[1,2] approximants yields F̃ /F̃s ≈ 0.68,
17 This is in line with a statement in [118] that a kink lies on the SUSY line when ∆σ 3 is slightly
smaller than 3.
152
which is quite close to c = 7/10. Therefore, in order to get a better estimate in d = 3,
it makes sense to impose it as an exact boundary condition in d = 2. Following this
F ≈ 0.158 . (6.130)
In the UV, we have the free CFT of one scalar and one Majorana fermion, which
1
has F UV = 4
log 2 ≈ 0.173, and therefore we nd F IR /F UV ≈ 0.91. This is a check
of the F -theorem for the ow from the free to the interacting N = 1 SCFT. It is
also interesting to compare the value of F at the SUSY xed point to the decoupled
Ising xed point in Figure 6.1. A plot comparing F̃ − F̃f with F̃Ising , which was
obtained in [135], is given in Figure 6.7. It shows that F̃ < F̃f + F̃Ising in the whole
range 2 < d < 4, in agreement with the generalized F -theorem [78, 135] and the
0.25
0.20
˜
0.15 F Ising
˜ ˜
F N=1 - F fer
0.10
0.05
Figure 6.7: Comparison of F̃ − F̃f at the SUSY xed point and F̃Ising (from [135]) in
2 ≤ d ≤ 4.
and we can use it to estimate the CT value at the N =1 SUSY xed point in d = 3.
We can perform a Padé approximation on the ratio CT /CT,s , where CT,s is the free
153
scalar value, using the boundary conditions
CT 7 in d = 2,
10
= (6.131)
CT,s
5 − 19
+ O(2 ) in d = 4 − ,
2 28
where we have imposed the d=2 value corresponding to the tri-critical Ising model.
CT 497d − 728
= , (6.132)
CT,s 62d + 256
which in d=3 gives CT /CT,s ≈ 1.73. This value lies well within the region allowed
by the bootstrap in g. 8 of [118] (note that ∆ψ ≈ 1.09). Comparing to the free UV
CFT of one scalar and one Ma jorana fermion, we have found CTIR /CTUV ≈ 0.86. It
would be useful to know the next order in the 4− expansion in order to obtain a
more precise estimate. Based on the comparison of the Padé for the N =2 model
with the exact result, we may expect the 2 correction to reduce CT somewhat in
d = 3.
Let us also include a brief discussion of the non-supersymmetric xed point of the
N =1 GNY model, marked by the red triangle in gure 6.1. While it has g2∗ < 0
in d = 4 − , it may become stable for suciently small d. Changing the sign of the
square root in (6.31), we nd the 4− expansion of the sphere free energy at this
xed point:
π 2 68753
F̃non−SUSY = F̃s + F̃f − + + O(4 ) . (6.133)
6 112 889056
charge 4/5 of the (5, 6) minimal model. This model includes primary elds of con-
formal weight 1/40 and 21/40, so that one may form a spin 1/2 eld with weights
154
(h, h̃) = (1/40, 21/40) that could correspond to ψ in the Yukawa theory (such a eld
is not present in the standard modular invariants that retain elds of integer spin
only).
18 In d=3 our extrapolation gives Fnon−SUSY ≈ 0.16. The latter quantity is
bigger than (6.130); therefore, the non-SUSY xed point, if it is stable, can ow to
supersymmetric xed point. Changing the sign of the square roots in (6.10) and
3 95 2
∆non−SUSY
σ =1− − + O(3 ) ,
7 6174
3 3 115 2
∆non−SUSY
ψ = − − + O(3 ) , (6.134)
2 7 37044
22 1117 2
∆non−SUSY
σ2 =2− + + O(3 )
21 9261
∆non−SUSY
σ ≈ 0.077 , ∆non−SUSY
ψ ≈ 0.63 , ∆σnon−SUSY
2 ≈ 0.297 . (6.135)
These numbers are not far from the corresponding operator dimensions in either the
(5, 6) or the (6, 7) minimal models. For example, in the (5, 6) interpretation the exact
scaling dimension of ψ in d = 2 should be 11/20, while in the (6, 7) it should be 19/28.
The Padé value we nd lies between the two, but the accuracy of the extrapolation
∆non−SUSY
ψ ≈ 1.068. Interestingly, these values are not far from the feature at
(∆σ , ∆ψ ) ≈ (0.565, 1.078), which lies below the supersymmetry line in gure 7 of
155
[118]. Changing the sign of the square root in (6.14), we also nd
13
∆σnon−SUSY
3 =3− + O(2 ) , (6.136)
7
which suggests that in d=3 the dimension of this parity-odd operator is less than 3.
156
Appendix A
calculation
Z
dd q 1
I(α, β) = d 2α 2β
Z 1(2π) q (p + q) Z
α−1 β−1 Γ (α + β) dd q 1
= dxx (1 − x)
Γ(α)Γ(β) (2π) [q + x(1 − x)p2 ]α+β
d 2
Z0 1 Z ∞ Z
xα−1 (1 − x)β−1 α+β−1 dd q −t(q2 +x(1−x)p2 )
= dx dtt e
0 Γ(α)Γ(β) 0 (2π)d
α+β− d
1 Γ d2 − α Γ d2 − β Γ α + β − d2 1 2
= d (A.1)
(4π) 2 Γ (α) Γ (β) Γ (d − α − β) p2
2 2− d
1 Γ 3 − d2 Γ d2 − 1 1 2
I1 ≡ I(1, 1) = d d
(A.2)
(4π) 2 (2 − 2 )Γ (d − 2) p2
p2 Γ(3 − d/2)
I1 → − . (A.3)
6(4π)3 (p2 )3−d/2
157
We will also need the following three-propagator integral
Z
dd k 1 1 1
I2 = (A.4)
(2π) (k − p) (k + q) k 2
d 2 2
Z 1 Z
dd k 1
=2 dxdydzδ(1 − x − y − z) (A.5)
0 (2π) [x(k − p) + y(k + q)2 + zk 2 ]3
d 2
Z 1 Z
dd k 1
=2 dxdydzδ(1 − x − y − z) (A.6)
0 (2π)d [k 2 − 2xkp + 2ykq + xp2 + yq 2 ]3
Dening l ≡ k − xp + yq , we get:
Z 1 Z
dd l 1
I2 = 2 dxdydzδ(1 − x − y − z) (A.7)
0 (2π) [l + ∆]3
d 2
with ∆ = x(1 − x)p2 + y(1 − y)q 2 + 2xypq . Evaluating the standard momentum
integral, we obtain
Z 1
1 Γ(3 − d/2)
I2 = dxdydzδ(1 − x − y − z) . (A.8)
0 (4π)d/2 ∆3−d/2
Z 1
1 Γ(3 − d/2) 1
I2 = dxdydzδ(1 − x − y − z) .
(4π)d/2 (M 2 )3−d/2 0 (x(1 − x) + y(1 − y) − xy)3−d/2
(A.9)
Z 1 1
2
dxdydzδ(1−x−y−z) 1 − log(x(1 − x) + y(1 − y) − xy) + O( ) = +O() .
0 2 2
(A.10)
Thus,
1 2 2
I2 = − log M + A + O() , (A.11)
2(4π)3
158
where A is an unimportant constant. For the purpose of extracting the log M 2 terms
1 Γ(3 − d/2)
I2 → . (A.12)
2(4π)3 (M 2 )3−d/2
159
Appendix B
The Feynman rules for our theory are depicted in Fig. B.1
α
α β
= −dαβγ = δαβ p12
β γ
α
α β
= −(δg)αβγ = −p2 (δz)αβ
β γ
where we introduced symmetric tensor coupling dαβγ and counterterms (δg)αβγ , (δz)αβ
with α, β, γ = 0, 1, .., N as
160
where i = 1, ..., N . The general form of a Feynman diagram in our theory could be
schematically represented as
The Tensor structure factors are products of the tensors dαβγ and (δg)αβγ , (δz)αβ ,
with summation over the dummy indices. Their values for dierent diagrams are
symmetry factors and are the same as in the usual ϕ3 -theory; their values are listed
1 To nd the Tensor structure factor we used the fact that it is a polynomial in N , so we
calculated sums of products of dαβγ , (δg)αβγ , (δz)αβ explicitly for N = 1, 2, 3, 4, ..., using Wolfram
Mathematica. Having answers for N = 1, 2, 3, 4, .. it's possible to restore the general N form.
161
B.1 Counterterms
φ g12 2 2
z14 = g1 (11N − 26) − 48g 1 g2 + 11g 2 ,
432(4π)6
σ 1 4 3 2 2 4
z14 =− 2N g 1 + 48N g 1 g 2 − 11N g 1 g 2 + 13g 2 ,
432(4π)6
1 2 3 2 2 3
a15 =− g 1 g 1 (11N + 98) − 2g 1 g2 (7N − 38) + 101g 1 g 2 + 4g2 ,
144(4π)6
1 5 4 3 2 2 3 5
b15 =− 4N g 1 + 54N g 1 g 2 + 18N g 1 g 2 − 7N g 1 g2 + 23g 2 ,
48(4π)6
(B.4)
φ g12
z16 = g 4 (N (13N − 232) + 5184ζ(3) − 9064) + g13 g2 6(441N − 544)
46656(4π)9 1
− 2g12 g22 (193N − 2592ζ(3) + 5881) + 942g1 g23 + 327g24 ,
σ 1
z16 =− 9
2N g16 (1381N − 2592ζ(3) + 4280) − 96N (12N + 11)g15 g2
93312(4π)
− 3N g14 g22 (N + 4320ζ(3) − 8882) + 1560N g13 g23 − 952N g12 g24 − g26 (2592ζ(3) − 5195) ,
g12
a17 = 9
− g15 (N (531N + 10368ζ(3) − 2600) + 23968)
15552(4π)
+ g14 g2 (99N 2 + 2592(5N − 6)ζ(3) − 9422N − 2588) + 2g13 g22 (1075N + 2592ζ(3) − 16897)
+ 2g12 g23 (125N − 5184ζ(3) − 3917) − g1 g24 (5184ζ(3) + 721) + g25 (2592ζ(3) − 2801) ,
1
b17 = − 2g17 N (577N + 713) − 48g16 g2 N (31N − 59)
2592(4π)9
+ g15 g22 N (423N + 2592ζ(3) + 1010) − g14 g23 N (33N − 1296ζ(3) − 6439)
− 27g13 g24 N (32ζ(3) + 11) − 301N g12 g25 + g27 (432ζ(3) + 1595) . (B.5)
162
Appendix C
Z
dd p 1 Ld (α, β)
I(α, β) = d 2α 2β
= 2 α+β−d/2 , (C.1)
(2π) p (p − k) (k )
where
1 Γ( d2 − α)Γ( d2 − β)Γ(α + β − d2 )
Ld (α, β) = . (C.2)
(4π)d/2 Γ(α)Γ(β)Γ(d − α − β)
For the more complicated integrals, we use the mathematica program FIRE [156],
which uses integration-by-parts (IBP) relations to turn them into simpler master
There are two categories of diagrams which show up quite frequently as subdia-
grams, the special KIT E diagrams and the ChT diagrams shown in Figure C.1.
163
1 1 α 1
SK(α) = α ChT (α, β) = 1
1 1 β 1
Figure C.1: The Special KIT E and ChT diagrams, the numbers labeling each prop-
integral:
Z
dd p dd q 1
SK(α) = . (C.3)
(2π) (2π) p (p + k) q (q + k)2 (p − q)2α
d d 2 2 2
Notice that the power of the middle propagator is arbitrary. Via the Gegenbauer
n=0
Γ(n + α + 1) (n + 1 − λ + α)
where λ = d/2 − 1. In the case of d=6− we have found that, for example:
d 1 1 −54 −71 + 24γ −14641 + 8520γ − 1440γ 2 + 120π 2
SK(2− ) = + + +. . .
2 (4π)d (k 2 )6−3d/2 2 1296 15520
(C.5)
The -expansion of the above result can also be veried indirectly with the mathe-
the integral:
Z
dd p dd q 1
ChT (α, β) = . (C.6)
(2π) (2π) p (p + k) q (q + k)2 (p − q)2
d d 2α 2β 2
164
In this diagram, one triangle of the KIT E diagram all have indices 1, and the other
two lines have arbitrary indices α and β. This diagram was evaluated in position
π d v(d − 2) 1
ChT (α, β) = d 2 d/2−3+α+β
Γ( 2 − 1) (x )
v(α)v(2 − α) v(β)v(2 − β) v(α + β − 1)v(3 − α − β)
× + + ,
(1 − β)(α + β − 2) (1 − α)(α + β − 2) (α − 1)(β − 1)
(C.7)
Γ(d/2−α)
where v(α) = Γ(α)
. For our purpose, we just need to fourier transform this
We also need variations of the SK and ChT diagrams, with a particular index
raised by 1, for example. However, we can use FIRE to relate them to the original
Z
dd pdd qdd r 1
LADDER = 3d
, (C.8)
(2π) p (p + k) (p − r) r (r + k)2 (r − q)2 q 2 (q + k)2
2 2 2 2
where the loop momenta are p, q , and r. The external momentum is k. Using FIRE,
165
r
p q
k p−r k MA MB MC
r−q
p+k q+k
r+k
MD ME
Figure C.2: The LADDER diagram can be reduced to ve master integrals
The diagrams corresponding to the master integrals are listed in Figure C.2. Among
these master integrals, only MD is non-primitive, the rest can be calculated easily.
However, if we integrate over the middle loop, we see that MD is in fact related to
k2 2 −115 + 36γ + log k 2
LADDER = − +
(4π)3d/2 93 1082
−4043 + 18(115 − 18γ)γ + 18π 2 − 18 log k 2 (−115 + 36γ + 18 log k 2 )
+ + ... .
1296
(C.11)
the same techniques used for the two-point functions, we impose that the momentum
running through the three points are p, −p, and 0, as a three-point function with
166
0 0 0
p p p p p p
Figure C.3: The three orientations of the same diagram topology correspond to dif-
ferent integrals.
However, since the momenta are asymmetric, it is necessary to consider all three
orientations of each topology of the diagram. Notice that the tensor factors men-
tioned in the previous section will also be dierent. As an illustration, let's denote the
into account that one of the external momenta is zero, they are each equivalent to a
two-point function as shown in Figure C.3. All lines have indices 1, except those with
I1 contains a subdiagram that is equivalent to ChT (1, 2), which can be evaluated
easily using our formula before; after that, the diagram is primitive. The other two
167
1 1 5 − 2γ − 2 log k 2
I1 = +
(4π)3d/2 63 82
173 + 18(γ − 5)γ − π 2 + 18 log k 2 (−5 + 2γ + log k 2 )
+ + ... (C.12)
96
1 1 5 − 2γ − 2 log k 2
I2 = +
(4π)3d/2 63 82
125 + 18(γ − 5)γ − π 2 + 18 log k 2 (−5 + 2γ + log k 2 )
+ + ... (C.13)
96
1 1 5 − 2γ − 2 log k 2
I3 = +
(4π)3d/2 63 82
125 + 18(γ − 5)γ − π 2 + 18 log k 2 (−5 + 2γ + log k 2 )
+ + ... . (C.14)
96
168
Appendix D
theories
In this appendix we list known general results for β- and γ -functions for general
1 1
L = (∂µ φi )2 + ψ̄ ∂/ψ + ψ̄Γi ψφi + gijkl φi φj φk φl , (D.1)
2 4!
component Dirac spinors, and the matrices Γi have the following general form
and act in the avor and spinor spaces and are not necessarily Hermitian. We also
the general Yukawa theoryfootnoteWe note that one can use results of [127, 128, 129]
169
for the four-component spinors, see sec. 4 in [128]. For the γ -functions the result
1 1 † 1 1 † † 3 † †
γψ = Γi Γi − Γi Γj Γj Γi + tr(Γi Γj )Γi Γj ,
(4π)2 2 (4π)4 8 8
1 1 1 3 1
† † † † †
γφ,ij = tr(Γi Γj ) + giklm gjklm − tr(Γj Γi Γk Γk ) − tr(Γj Γk Γi Γk ) .
2(4π)2 (4π)4 12 4 2
(D.3)
For the anomalous mixing matrix of the Nb (Nb + 1)/2 quadratic operators Oij = φi φj
with i6j we have [131]
M ij,kl + Mij,lk , k 6= l ,
γij,kl = γφ,mk δij,ml + γφ,ml δij,km + (i 6 j, k 6 l) ,
Mij,kk , k = l,
(D.4)
where δij,kl is the Kronecker delta (e.g. δ11,11 = 1, δ11,12 = 0, δ12,12 = 1, δ22,22 = 1,. . . )
and
1 1 † † †
Mij,kl = g ijkl − gikmn gjlmn + tr(Γl Γm )g ijkm − 2tr(Γi Γ k Γj Γl ) . (D.5)
(4π)2 (4π)4
For the β -functions we have (see (3.3) in [128] and (7.2) in [130])
1 1 2† † 1 †
1
βi = 2
2
(Γ Γi + Γi Γ ) + 2Γj Γi Γj + Γj tr(Γj Γi ) + 4
2Γk Γ†j Γi (Γ†k Γj − Γ†j Γk )
(4π) 2 2 (4π)
1
− Γj (Γ2 Γ†i + Γ†i Γ2† )Γj − (Γj Γ2 Γ†j Γi + Γi Γ†j Γ2† Γj ) − tr(Γ†i Γk )Γj Γ†k Γj
8
3 3 1
− tr(Γ†j Γk )(Γj Γ†k Γi + Γi Γ†k Γj ) − Γj tr (Γ2 Γ†j + Γ†j Γ2† )Γi + Γ†j Γk Γ†i Γk
8 8 2
1
− 2gijkl Γj Γ†k Γl + giklm gjklm Γj (D.6)
12
170
and (see (4.3) in [129])
1 1 X 1X † † 1 X †
βijkl = gijmn gmnkl − tr(Γi Γj Γk Γl ) + tr(Γa Γa )gijkl
(4π)2 8 perm 2 perm 2 a=i,j,k,l
1 1 X 1X 1X †
+ g ijkl gamnp gamnp − g ijmn gkmpr glnpr − Tr(Γn Γp )gijmn gklmp
(4π)4 12 a=ijkl
4 perm 8 perm
1X X 3 1
+ gijmn tr(Γk Γ†m Γl Γ†n ) − gijkl tr(Γa Γ†a Γm Γ†m ) + tr(Γa Γ†m Γa Γ†m )
2 perm a=ijkl
4 2
X
† † † † † † † † †
+ tr(Γm Γm Γi Γj Γk Γl ) + 2tr(Γm Γi Γm Γj Γk Γl ) + tr(Γi Γj Γm Γk Γl Γm ) ,
perm
(D.7)
P 2 †
where
perm denotes the sum over all permutation of the indices i, j, k, l and Γ = Γi Γi
P
and Γ2† = Γi Γ†i , also
a=ijkl f (a) ≡ f (i) + f (j) + f (k) + f (l). Traces tr(Γi . . . ) are
generalize the result of [130] for βb , when the matrices Γi are not Hermitian:
1 1 1
βb = − tr(Γ†i Γ2 Γi Γ2 ) + tr(Γ†i Γj Γ†i Γj Γ2 ) + tr(Γi Γ†j Γi Γ†k Γj Γ†k ) − tr(Γi Γ†j Γk Γ†i Γj Γ†k )
(4π)8 144 8
3 1
+ tr(Γ†i Γj )tr(Γ†i Γj Γ2 + Γ†i Γk Γ†j Γk ) + gijkl tr(Γi Γ†j Γk Γ†l ) − giklm gjklm tr(Γ†i Γj ) .
4 24
(D.8)
Now to use these formulas for the GNY model one simply takes
171
and nds the results (6.7), (6.9), (6.11) and (6.19). In the case of the NJL model, we
have
172
∂
Value of ∂p2
(sum of graphs)
1 1 Ng12 + g22
(a) − (4π)3 6ǫ ×
2g12
(
173
1− ×
2g13(2g13 + Ng12 g2 + g23 )
(
Figure D.1: Values of derivatives of two-point diagrams. The upper row in parenthesis is for hσσi and the lower is for hφφi.
Figure D.2: Values of three-point diagrams. The upper row in parenthesis is for hσσσi and the lower is for hσφφi.
1 1 Ng 3 + g23
(a) − (4π)3 ǫ × 2 1
g1 (g1 + g2 )
(
1 1 g2 (3Ng14 + g24 )
(d) − (4π)6 4ǫ ×
2g13(g12 + g22 )
(
1 1 3ǫ 3ǫ2 3 (N(N + 1)g17 + Ng16 g2 + Ng15 g22 + 2Ng14 g23 + Ng13 g24 + g27)
174
(e) − (4π)9 6ǫ3 1 − 4
+ 4 × 2
g1 ((2N + 3)g15 + (N + 7)g14g2 + (N + 6)g13g22 + 3g12 g23 + g25)
(
1 1 29ǫ 19ǫ2 3 (N(N + 2)g17 + Ng15 g22 + 2Ng14 g23 + Ng12 g25 + g27 )
(m) − (4π)9 108ǫ3 1 − 12
+ 9 × 2
g1 ((4N + 4)g15 + (5N + 2)g14 g2 + (N + 2)g13g22 + (2N + 1)g12g23 + g1 g24 + 2g25)
(
175
1 1 11ǫ 3 (N(N + 4)g16 g2 + Ng14 g23 + Ng12 g25 + g27 )
(p) − (4π)9 36ǫ2 1− 12 ×
2g13 ((N + 4)g14 + (2N + 3)g12 g22 + 2g24)
(
Figure D.3
Bibliography
[2] K. Wilson and J. B. Kogut, The Renormalization group and the epsilon ex-
[3] M. Moshe and J. Zinn-Justin, Quantum eld theory in the large N limit: A
the Critical Indices in the 1/ N Expansion, Theor.Math.Phys., vol. 46, pp. 104
113, 1981.
[7] K. Lang and W. Ruhl, Field algebra for critical O(N ) vector nonlinear sigma
models at 2 < d < 4, Z.Phys., vol. C50, pp. 285292, 1991.
176
[8] K. Lang and W. Ruhl, The Critical O(N ) sigma model at dimension 2 <
d < 4 and order 1/N 2 : Operator product expansions and renormalization,
[9] K. Lang and W. Ruhl, The Critical O(N) sigma model at dimensions 2 < d < 4:
A List of quasiprimary elds, Nucl.Phys., vol. B402, pp. 573603, 1993.
[10] K. Lang and W. Ruhl, The Critical O(N) sigma model at dimensions 2 < d < 4:
Fusion coecients and anomalous dimensions, Nucl.Phys., vol. B400, pp. 597
623, 1993.
[12] A. C. Petkou, C(T) and C(J) up to next-to-leading order in 1/N in the con-
formally invariant 0(N) vector model for 2 < d < 4, Phys.Lett., vol. B359,
and conformally covariant operator product expansion, Annals Phys., vol. 76,
177
[17] S. El-Showk, M. F. Paulos, D. Poland, S. Rychkov, D. Simmons-Dun,
et al., Solving the 3D Ising Model with the Conformal Bootstrap, Phys.Rev.,
Models, 2013.
[23] S. Giombi, I. R. Klebanov, and B. R. Safdi, Higher Spin AdSd+1 /CFTd at One
presence of double trace deformations, Nucl.Phys., vol. B656, pp. 2336, 2003.
lence of the top quark condensate and the elementary Higgs eld, Nucl.Phys.,
178
[27] J. Zinn-Justin, Four fermion interaction near four-dimensions, Nucl.Phys.,
[28] L. Fei, S. Giombi, and I. R. Klebanov, Critical O(N ) Models in 6− Dimen-
of the critical O(N) models in 6− dimensions, Phys. Rev., vol. D91, no. 4,
p. 045011, 2015.
in 6− dimensions and higher spin dS/CFT, JHEP, vol. 09, p. 076, 2015.
tum Field Theory and Statistical Mechanics Cargèse 1976, pp. 281305, Springer
US, 1977.
[34] X. Bekaert, E. Meunier, and S. Moroz, Towards a gravity dual of the unitary
2013.
179
[37] M. Fisher, Yang-Lee Edge Singularity and phi**3 Field Theory,
the Percolation Problem and the Yang-lee Edge Singularity, J.Phys., vol. A14,
p. 2391, 1981.
sitions in Potts φ3 eld theory with quadratic and trilinear symmetry breaking,
Phys.Rev., vol. B34, pp. 31653176, 1986.
[42] M. A. Vasiliev, Consistent equation for interacting gauge elds of all spins in
all spins in (3+1)-dimensions, Phys. Lett., vol. B285, pp. 225234, 1992.
[44] M. A. Vasiliev, Higher-spin gauge theories in four, three and two dimensions,
[45] M. A. Vasiliev, Higher spin gauge theories: Star-product and AdS space, 1999.
[46] M. Vasiliev, Nonlinear equations for symmetric massless higher spin elds in
[47] I. R. Klebanov and A. M. Polyakov, AdS dual of the critical O(N ) vector
180
[48] P. Breitenlohner and D. Z. Freedman, Stability in Gauged Extended Super-
[51] E. Sezgin and P. Sundell, Holography in 4D (super) higher spin theories and
a test via cubic scalar couplings, JHEP, vol. 07, p. 044, 2005.
2011.
[55] D. J. Broadhurst and A. Kotikov, Compact analytical form for nonzeta terms
[56] M. Hasenbusch, Finite size scaling study of lattice models in the three-
dimensional Ising universality class, Phys. Rev. B, vol. 82, p. 174433, Nov.
2010.
A. Vichi, Solving the 3d Ising Model with the Conformal Bootstrap II. c-
181
Minimization and Precise Critical Exponents, J. Stat. Phys., vol. 157, p. 869,
2014.
[58] B. Halperin, T. Lubensky, and S.-k. Ma, First order phase transitions in super-
conductors and smectic A liquid crystals, Phys.Rev.Lett., vol. 32, pp. 292295,
1974.
[59] I. Herbut et al., Herbut and tesanovic reply, Physical Review Letters, vol. 78,
p. 980, 1997.
[60] F. Gliozzi and A. Rago, Critical exponents of the 3d Ising and related models
[61] Y. Nakayama and T. Ohtsuki, Five dimensional O(N )-symmetric CFTs from
1995.
for general dimensions, Annals Phys., vol. 231, pp. 311362, 1994.
JHEP, vol. 1101, p. 125, 2011. 62 pages, 4 gures, few typo's corrected.
182
[68] D. E. Diaz and H. Dorn, Partition functions and double-trace deformations in
Lattice?, 2014.
[72] R. Percacci and G. P. Vacca, Are there scaling solutions in the O(N )-models
for large N in d>4 ?, 2014.
[74] J. L. Cardy, Conformal Invariance and the Yang-lee Edge Singularity in Two-
[75] D. Amit and D. Roginsky, Exactly Soluble Limit of φ3 Field Theory with
183
[79] D. J. Broadhurst, J. Gracey, and D. Kreimer, Beyond the triangle and unique-
[83] T. Quella, I. Runkel, and G. M. Watts, Reection and transmission for con-
from Non-Critical String Theory, Phys. Lett., vol. B428, pp. 105114, 1998.
[86] E. Witten, Anti-de Sitter Space and Holography, Adv. Theor. Math. Phys.,
[87] S. Giombi and X. Yin, The Higher Spin/Vector Model Duality, J.Phys.,
[88] A. Strominger, The dS / CFT correspondence, JHEP, vol. 0110, p. 034, 2001.
184
[90] A. LeClair, Quantum critical spin liquids, the 3D Ising model, and conformal
[92] A. Higuchi, Forbidden Mass Range for Spin-2 Field Theory in De Sitter Space-
[93] A. Higuchi, Symmetric Tensor Spherical Harmonics on the N Sphere and Their
Application to the De Sitter Group SO( N ,1), J.Math.Phys., vol. 28, p. 1553,
1987.
[98] J.-B. Bae and S.-J. Rey, Conformal Bootstrap Approach to O(N) Fixed Points
Fermionic eld theory for trees and forests, Phys.Rev.Lett., vol. 93, p. 080601,
2004.
185
[100] Y. Deng, T. M. Garoni, and A. D. Sokal, Ferromagnetic Phase Transition for
the Spanning-Forest Model ( q→0 Limit of the Potts Model) in Three or More
[102] H. Saleur and B. Wehefritz Kaufmann, Integrable quantum eld theories with
supergroup symmetries: The OSP (1/2) case, Nucl.Phys., vol. B663, p. 443,
2003.
Based on an Analogy with Superconductivity. 1., Phys. Rev., vol. 122, pp. 345
358, 1961.
[105] E. Witten, Chiral Symmetry, the 1/N Expansion, and the SU (N ) Thirring
[106] A. Bondi, G. Curci, G. Pauti, and P. Rossi, Metric and Central Charge in the
[108] K.-i. Kubota and H. Terao, Dynamical symmetry breaking in QED(3) from
the Wilson RG point of view, Prog. Theor. Phys., vol. 105, pp. 809825, 2001.
Neveu model, Int. J. Mod. Phys., vol. A6, pp. 395408, 1991. [Erratum: Int.
J. Mod. Phys.A6,2755(1991)].
186
[110] J. A. Gracey, The Nambu-Jona-Lasinio model at O(1/N 2 ), Phys. Lett.,
[111] S. Hands, A. Kocic, and J. B. Kogut, Four Fermi theories in fewer than four-
[112] B. Rosenstein, H.-L. Yu, and A. Kovner, Critical exponents of new universality
graphene and topological insulators, Phys. Rev., vol. B87, no. 4, p. 041401,
2013.
3-d Gross-Neveu and Higgs-Yukawa models, Nucl. Phys., vol. B415, pp. 781
187
[120] R. C. Myers and A. Sinha, Seeing a C-Theorem with Holography, Phys. Rev.,
symmetry at the Boundary of a Topological Phase, Science, vol. 344, no. 6181,
[126] H. Shimada and S. Hikami, Fractal dimensions of self-avoiding walks and Ising
[130] I. Jack and H. Osborn, Analogs for the c Theorem for Four-dimensional Renor-
malizable Field Theories, Nucl.Phys., vol. B343, pp. 647688, 1990.
188
[131] M. Pernici, M. Raciti, and F. Riva, Dimensional renormalization of Yukawa
theories via Wilsonian methods, Nucl. Phys., vol. B577, pp. 293324, 2000.
[133] S. Hathrell, Trace Anomalies and λφ4 Theory in Curved Space, Annals Phys.,
[134] I. Jack and H. Osborn, General Background Field Calculations With Fermion
[140] J. A. Gracey, Four loop MS-bar mass anomalous dimension in the Gross-Neveu
1994.
189
[142] J. A. Gracey, Computation of critical exponent eta at O(1/N**3) in the four
Fermi model in arbitrary dimensions, Int. J. Mod. Phys., vol. A9, pp. 727744,
1994.
index eta in order 1/n**3, Theor. Math. Phys., vol. 94, pp. 127136, 1993.
2-D four-fermion interactions, Theor. Math. Phys., vol. 107, pp. 441455, 1996.
symmetric model, Teor. Mat. Fiz., vol. 107N3, pp. 359371, 1996. [Theor.
Math. Phys.107,710(1996)].
O(1/N 2 ), Phys. Rev., vol. D50, pp. 28402847, 1994. [Erratum: Phys.
Rev.D59,109904(1999)].
Dimensional Supersymmetric Ising Model, Phys. Rev. Lett., vol. 115, no. 5,
p. 051601, 2015.
190
[149] N. Bobev, S. El-Showk, D. Mazac, and M. F. Paulos, Bootstrapping SCFTs
Strange Metals in One Spatial Dimension, Phys. Rev., vol. D86, p. 066003,
2012.
face states with emergent supersymmetry, Phys. Rev. Lett., vol. 116, no. 10,
p. 100402, 2016.
Dimensions and the Tricritical Ising Model, Phys. Lett., vol. B151, pp. 3743,
1985.
191