Monopole operators from the $4-\epsilon$ expansion

Three-dimensional quantum electrodynamics with $N$ charged fermions contains monopole operators that have been studied perturbatively at large $N$. Here, we initiate the study of these monopole operators in the $4-\epsilon$ expansion by generalizing them to codimension-3 defect operators in $d = 4-\epsilon$ spacetime dimensions. Assuming the infrared dynamics is described by an interacting CFT, we define the"conformal weight"of these operators in terms of the free energy density on $S^2 \times \mathbb{H}^{2-\epsilon}$ in the presence of magnetic flux through the $S^2$, and calculate this quantity to next-to-leading order in $\epsilon$. Extrapolating the conformal weight to $\epsilon = 1$ gives an estimate of the scaling dimension of the monopole operators in $d=3$ that does not rely on the $1/N$ expansion. We also perform the computation of the conformal weight in the large $N$ expansion for any $d$ and find agreement between the large $N$ and the small $\epsilon$ expansions in their overlapping regime of validity.


.1 Motivation and setup
A fascinating class of relatively simple, yet sufficiently non-trivial, quantum field theories (QFTs) in 2 + 1 space-time dimensions can be obtained by coupling a U (1) gauge theory to charged matter. These theories share certain features with 3+1 dimensional QCD, such as asymptotic freedom and flowing to strong coupling in the infrared, but are easier to study because their gauge group is abelian. In addition to the local gauge-invariant operators that can be constructed as polynomials in the fundamental fields appearing in the Lagrangian and their derivatives, such QFTs contain also monopole operators [1]. These are local, gauge invariant operators distinguished by the fact that they carry non-zero charge under a topological global symmetry group U (1) top , whose conserved current is Here, F νρ = ∂ ν A ρ − ∂ ρ A ν is the field strength of the Abelian gauge field A µ . The current (1.1) is conserved due to the Bianchi identity satisfied by the gauge field strength. The existence of operators charged under U (1) top is intimately tied to the non-trivial topology of the gauge group, in particular to its non-vanishing fundamental group π 1 (U (1)) ∼ = Z. The Dirac quantization condition implies that the conserved charge q = d 2 x j 0 (1.2) satisfies q ∈ Z/2, or equivalently, that the integrated magnetic flux through a small two sphere surrounding a monopole operator is 4πq ∈ 2πZ. Monopole operators can play an important role in the dynamics of these theories-see, for instance, [2][3][4][5][6][7][8][9][10][11][12][13][14][15] for examples arising in various condensed matter systems.
Monopole operators are remarkably difficult to study. Even when the QFT of interest has a weakly coupled description, which is the case when the number of charged matter fields is very large, monopole operators cannot be studied using traditional techniques. Progress can be made in the deep infrared (IR), if one assumes that the IR limit of the renormalization group (RG) flow is described by a conformal field theory (CFT) with unbroken U (1) top sym-metry. One can then use the state-operator correspondence to map a monopole operator of topological charge q inserted at x = 0 to the ground state on S 2 × R in the presence of magnetic flux equal to 4πq though the S 2 [16]. The scaling dimension of the monopole operator is determined by the ground state energy (or free energy) on S 2 × R. 1 At weak coupling, this ground state energy can be computed by performing a saddle point approximation around the configuration where the magnetic flux is uniformly distributed on the S 2 .
Monopole operators in conformal field theories have been studied using the 1/N expansion [16,18,19]. BPS monopole operators, which can exist in theories with N ≥ 2 supersymmetry, have been studied using supersymmetric localization [20,21]. In this paper, we study monopole operators using the 4 − expansion. 23 In this approach, one formally continues the CFT to d = 4 − dimensions and takes to be small. The 4 − expansion is a well-established method for studying more conventional operators [27]. It has recently been used to provide an approximation for the S 3 free energy of various 3d CFTs with vector-like matter [28][29][30][31][32][33]. To the best of our knowledge, the only other example of a disorder operator that has been studied in the 4 − expansion is the twist defect of the 3d Ising model [34].
The generalization of monopole operators to d = 4− dimensions requires a little thought.
In going away from d = 3, it is natural to retain the notion of a conserved U (1) top quantum number and to think of a charge q monopole operator as a codimension-3 defect operator that creates magnetic flux 4πq through the S 2 that surrounds the codimension-3 defect. In d = 4 the conserved current in (1.1) becomes a conserved 2-form equal to the Hodge dual of the gauge field strength-see [35] for a description of such higher-form symmetries. The generalized monopole operator is a 't Hooft line operator that carries non-vanishing charge under a global symmetry associated to a conserved current two-form.
Having defined these codimension-3 operators in d = 4 − dimensions, we would also like to define a "conformal weight" that characterizes the transformation properties of these operators under dilatations and agrees with the scaling dimension of the monopole operators when = 1. As in [36], we consider a planar defect operator extending along a flat R 1− that passes through the origin. Using cylindrical coordinates around it, we can make a conformal 1 The scaling dimensions of monopole operators were first computed in [17] in the CP N b −1 model without making use of the state-operator correspondence. The computation in [17] was only at leading order in 1/N b . 2 It is also interesting to approach the problem of monopole scaling dimensions using numerical methods such as Quantum Monte Carlo [22,23] and the conformal bootstrap [?, 24]. 3 Besides the 1/N and expansions, monopole operators can be studied in the recently proposed large charge expansion [26], where it was observed that a subleading term in the large charge limit of operator dimensions can be calculated using effective field theory considerations. It would be very interesting to understand if it is possible to calculate the coefficients of the other terms in the large charge expansion using this method. mapping ds 2 = ds 2 Note that the defect operator sitting at r = 0 is mapped to the conformal boundary of H 2− . This configuration has magnetic flux 4πq through S 2 , and some free energy density F (q) / Vol(H 2− ) in the presence of this flux. We normalize this quantity by considering the difference in the free energy density with the curvature radii set to 1 between the configuration with magnetic flux 4πq threading S 2 and the vacuum. In four dimensions, this difference gives the scaling weight h q of the 't Hooft line operator [36], while in three dimensions it gives the scaling dimension ∆ q of the monopole operator. 4 In the rest of this paper, we restrict our attention to quantum electrodynamics with N flavors of two-component complex fermions of unit charge under the U (1) gauge group. 5 We choose N to be even in order to avoid a parity anomaly in d = 3 [38][39][40]. The bare Lagrangian of this theory on a d-dimensional curved manifold is where e 0 is the bare gauge coupling, and we grouped our N two-component fermions into N/2 four-component fermions.

Summary of results
Before delving into the details of our computation, let us summarize our results. 4 Note that in three dimensions (1.3) describes S 2 × R, hence we recover the conventional state-operator correspondence. 5 The theory of N b complex scalars coupled to a U (1) gauge theory only has a non-trivial fixed point in the 4 − expansion for N b ≥ 183 [37], which makes the 4 − expansion unsuitable for studying such theories at small N b .

Results from the 4 − expansion
The result of the 4 − expansion is reached in two steps. The first step is to expand the free energy at small e 2 0 . The expansion takes the form: where, as we will explain, these terms represent contributions from the classical action, the one-loop determinant of the fermions, the two-loop current-current vacuum diagram, etc. All these terms are explicit functions of , and we obtain them by first expanding the relevant diagrams in and then using zeta-function regularization to regularize any divergences. We only obtain these terms up to some small order in the small expansion. These terms also depend on the ultraviolet (UV) cutoff scale µ that appears implicitly in the definition of the theory and on the curvature radii of S 2 and H 2− , which are both taken to be equal to R.
The dependence on µ and R is only through the dimensionless combination µR.
The second step in our calculation is to find the value of e 0 corresponding to an RG fixed point and plug it in (1.6). This critical value is customarily found from the vanishing of the beta function. However, beta functions are often scheme dependent, and one would have to make sure that the scheme used for obtaining the beta function matches the regularization procedure used in other parts of the computation. To circumvent this potential annoyance, we find the critical value of the bare coupling e 0 by requiring that ∆F (q) is µ-independent.
From this requirement, we obtain: where γ is the Euler-Mascheroni constant. By plugging this into (1.6) we turn the double expansion in e 0 and into an expansion in only : We find that the first two coefficients are and f(q) is given in Table 1 for several values of q (a general formula is given by (B.4) in Appendix B). In obtaining (1.7), we needed to determine the divergent piece of f Our setup contains two parameters that can be used to gain some insight into the resummation of the 4 − expansion, namely the number of matter fields N and the monopole charge q, which we now discuss.

Large N limit
The free energy on S 2 × H 2− can be evaluated in the large N expansion for any fixed d = 4 − , with not necessarily small. In d = 3, this approach was used for theories with fermionic matter in [18,19]. The structure of the expansion is: From the consistency of (1.10) with (1.8) in the overlapping region of large N and small we expect that (1.11) In Section 3, we compute explicitly g (q) 0 ( ), and by expanding it at small we do indeed recover (1.11). In other words, the N → ∞ and → 0 limits commute. The fact that these limits commute may seem rather trivial, but the derivation we present uncovers an interesting subtlety. Indeed, if > 0, the Maxwell term in the action (1.5) is irrelevant and has to be dropped, hence the entire contribution to g Nevertheless, expanding the regularized fermion determinant at small , we do reproduce the classical contribution of the Maxwell term. We find this agreement quite remarkable! Since the large N and small expansions commute, it makes sense to consider the double expansion in 1/N and . The non-zero terms in this expansion are shown in Table 2. As   can be seen from this table, at fixed order k the 1/N expansion is finite; it starts with a term of order N and ends with a term of order 1/N k . (This pattern is already visible in the expansions given above in (1.8) and (1.10).) The fact that at fixed , the 1/N expansion terminates is a property of the fermionic theory-such a feature would not arise in similar computations in theories with bosonic matter.
The difference between the large N and small expansions is in how precisely we resum the terms in Table 2. In the large N expansion at fixed , one essentially resums each row of Table 2 first and then extrapolates the sum of the first few rows to finite N . (See [18,19] for results in d = 3 from the large N expansion.) In the small expansion at fixed N , one resums each column of Table 2 first and then extrapolates the sum of these columns to finite . It is of course possible to contemplate other ways of resumming the terms in Table 2-for instance, first resumming each NE-SW diagonal or, alternatively, first resumming each NW-SE diagonal. Since the full resummation is believed to be only an asymptotic series, it could happen, in principle, that some ways of resumming the entries of Table 2 would provide better approximations than others. 6

Large q limit
Our setup has another parameter, q, which takes only discrete values but which can be taken to be large. In the limit q → ∞ we can derive some analytic formulas. In [41] it was pointed out that q → ∞ is a flat space limit: it corresponds to a uniform large magnetic field on S 2 , which makes particles move on highly localized Landau levels that are ignorant about the curvature of S 2 . 7 Then from flat space dimensional analysis we conclude that Ref. [26] gave an effective field theory argument to the same effect in d = 3.
Let us see how our results match these expectations. Although we only know f(q) in terms of a sum and integral, we can derive its large q behavior, from which we get: (1.13) 6 It is worth noting that in [28] it was noticed that in the same theory as the one considered here, the expansion and the large N expansion resummations of the free energy on S 4− give similar results. In the double expansion in and 1/N , the S 4− free energy has a similar structure as the free energy on S 2 × H 2− in the presence of monopole flux: in particular, for fixed order in , the 1/N expansion contains only a finite number of terms. So based on this evidence, one may conjecture that the resummations of the free energy on S 2 × H 2− using the expansion and the 1/N expansion would give similar results. 7 The particles are free to move in the rest of the 1 − spatial directions, but this doesn't change our conclusions.
where γ is the Euler-Mascheroni constant and A is the Glaisher constant. From this expression we readily see that q 2 1 − 2 log q + . . . resums to (1.12), and we obtain an expansion for the -dependent prefactor. (Of course, this analysis does not say anything about what the correct resummation for the prefactor might be.) The large-q expansion thus provides another clue as to how to resum the 4 − expansion.
It is worth mentioning that the first few term of the large q expansion written down in (C.4) approximates the f(q) values in Table 1 to 4 digits precision even at q = 1/2, so in practice (C.4) is valid for all allowed values of q. We do not currently have an explanation for this observation.

Organization of the paper
The rest of this paper is organized as follows. Section 2 represents the main part of this paper where we compute the free energy density in the presence of a monopole background on S 2 × H 2− . In Section 3 we perform the leading large N expansion computation for arbitrary and compare it to the computation of Section 2. In Section 4 we remark on how the defect operators transform under the global symmetries of the theory. Several technical details of our computation are included in the Appendices.

Conventions
We will investigate U (1) gauge theories with N/2 four-component Dirac fermions of unit charge. In d = 4, these theories have an SU (N/2) L × SU (N/2) R flavor symmetry. In d = 3, the symmetry is enhanced to SU (N ). The defect operators are expected to transform in representations of this symmetry. 8 The action on a manifold of dimension d = 4 − is where ∇ µ is the covariant derivative compatible with the background metric g µν , ψ α are the N/2 Dirac fermions 9 , and the gauge field has been written as the sum of a background A µ , with field strength F µν , and a small fluctuation a µ , with field strength f µν , around this 8 We discuss global symmetries more extensively in Section 4. 9 Recall that N has to be an even number in order to avoid a parity anomaly in d = 3.
background. As usual for a Wilson-Fisher type expansion, we have included all terms that are marginal in 4d, the dimension that we are perturbing away from.
We are interested in studying this theory on S 2 × H 2− in the background of magnetic flux 4πq through S 2 . We may write the metric as where R S and R H denote the curvature radii of S 2 and H 2− , respectively. We are eventually interested in taking these curvature radii equal, R H = R S , because that is the case that is related by a conformal transformation to R 4− , but for now we will keep them distinct.
In order to have magnetic flux 4πq through S 2 , we choose a background gauge field given by: The expression for A is well-defined everywhere away from θ = π. The singularity at θ = π is not physically observable provided that the Dirac quantization condition q ∈ Z/2 is obeyed.
Since we will be working with spinors on a curved manifold, we must specify our conventions for the frame and gamma matrices. When applying the expansion to spinors, which must be defined in integer dimensions, it is standard practice to use the spinor conventions of the = 0 dimension, i.e. 4d in our case. In 4d, we can write the metric (2.2) explicitly in geodesic polar coordinates as We then define our frame as We define the gamma matrices on the 4d space S 2 × H 2 as a tensor product of the standard 2d gamma matrices (Pauli matrices σ 1 and σ 2 ) on the separate spaces S 2 and H 2 :

First step: Setup of the loop expansion
In the following, we use the 4 − expansion to compute the free energy density The 4 − expansion is a procedure in which we expand every quantity in first, and then regularize divergences and extract the renormalized quantities. We will see in Section 3 that this order of steps is rather crucial.
The free energy on S 2 × H 2− can be written as 10 Because the free energy is dimensionless and we want to obtain a dimensionless number for the density, by Vol(M ), with M = H 2− or S 2 , we always mean the volume of the unit-radius M . 11 More precisely, we are only considering gauge field fluctuations that alone have vanishing total flux. Then a is a well defined 1-form, hence f = da is a total derivative and integrates to zero. (Recall that F is proportional to the volume form on S 2 .) where the linear in a µ term is absent because A is a saddle point, and denotes, up to an overall sign, the current-current two-point function in the theory of a free Dirac fermion in the monopole flux background on S 2 × H 2− . Plugging this expansion back into (2.8), and integrating over the gauge field fluctuations, we find that we can write the free energy as: The expression (1.6) can then be calculated from (2.12).
Note that in (2.12) we absorbed the q-independent functional determinant coming from the term f 2 /4e 2 0 in the path integral measure for a µ , so that it does not show up in (2.12). In the following we will explicitly calculate the first two terms and draw some conclusions about the third term. We emphasize again that these terms should be expanded in first, then regularized and renormalized.

Classical contribution
The classical contribution of the monopole background to the free energy is at all orders in given by (2.13) In the case of H 2− , the volume Vol(H 2− ) is divergent, but the free energy density in (2.13) is finite. The resulting free energy density is dimensionless because [e 2 0 ] = in 4 − dimensions.
For future reference, when R S = R H = R, the quantity f (2.14) The free energy F (q) /Vol(H 2− ) interpolates between the scaling dimension ∆ q of monopole operators in 3d and the scaling weight h q of line operators in 4d as defined in [36]. Indeed, it can easily be checked that the definition in [36] agrees with (2.13) when = 0 as follows.
From (2.13), we have Quite generally, the same quantity can be computed from the definition of T µν as where δg µν is the change in the inverse metric when R S and R H change by δR S and δR H , respectively. From [36] we have that the scaling weight of monopole operators is defined through where the curvature radii are set to R 2 S = R 2 H = R 2 , and it equals h q = q 2 2e 2 0 in pure Maxwell theory. Plugging this information into (2.16) and setting the curvature radii equal, one immediately reproduces (2.15). 12

One loop contribution at O( 0 )
The next term in (2.12) is the fermionic 1-loop determinant. In this section we calculate it to leading order in , namely 0 , while in Section 2.5 we extract some information at nextto-leading order. Hence, corresponding to setting = 0, we want to calculate the following 12 In terms of the free energy F (q) (R S , R H ) the scaling weight h q is given by: 13 determinant in four dimensions: where the trace is to be computed on the space S 2 × H 2 . We will use zeta function regularization to regulate the divergences in this trace.
We will find that the functional determinant depends on the logarithm of an arbitrary scale µ; the coefficient of log µ will be found to be the trace anomaly. Note, however, that is a well-defined quantity in a CFT, so it should be independent of µ.
In Section 2.7 we will indeed find that the residual log µ dependence, which remains after the subtraction of the vacuum energy F (0) , gets absorbed into the expression for the bare coupling at the conformal fixed point order by order in . With this preview, let us turn to a brief review of the expected trace anomaly.
Expectation from trace anomaly. The trace anomaly for N Weyl fermions is [42]: where β(1/e 2 ) is the beta function of 1/e 2 . 13 If we integrate the trace anomaly over S 2 × H 2 , we get the response of the free fermions to the rescaling of µ. 14 Then a standard argument 13 The more conventional form of the beta function is 14 On S 2 × H 2 the curvature invariants are calculated from the metric (2.4) to be in our monopole background.
implies that the free energy contains a logarithmic term in µ: where, in order to avoid clutter, we introduced the notation We will see that our explicit computation of F  On S 2 , we can find an orthonormal basis of eigenspinors of (i / ∇ + / A) obeying  2 (λ, n), which is in fact the same for each mode, so that the combined density of states is where we assumed that the n = 0 modes satisfy ψ (I) λ0 (0) = 1. We will not need an expression for µ 2 (λ, n) for n > 0.
where x is any point of our choosing. 15 We can take x to correspond to θ = η = 0 in the coordinates given in (2.4). With this choice, we can insert a complete set of modes in (2.33) and find that only the modes with m = n = 0 contribute. The density of these modes is: and hence the free energy is given by: Here, µ is the UV cutoff scale we introduced to make the formulas dimensionally correct, and the argument of the logarithm is just the (absolute value of the) determinant of (2.31).
Although λ ∈ R, we only integrate over R + , which cancels a 1/2 factor coming from the fact that the determinant of (2.31) corresponds to the contribution of two modes.
Regularization of the functional determinant. Both the sum and the integral in (2.36) are divergent. To regularize them, it is convenient to introduce a slight variation of the usual zeta function regularization procedure [43] and write − log A, for some quantity A, as where D is a linear functional defined by where any n > m in (2.37) is a good choice. In fact, we will only act with D on Laurent series, so (2.39) can be regarded as the definition of the action of D.
We will thus rewrite (2.36) as Let us focus on I (q) first. After performing the λ integral, we have: where we have separated out the = q term from the rest. Next, we add and subtract The full expression when r = 1 is given in Appendix B.
The λ integrals in II (q) are convergent, so we need only regularize the divergent sum over , which we do as in I (q) . The answer is again too lengthy to be reproduced here, but it takes the form II (q) = − N (35 + 93r 2 ) 360π log(µR H ) + terms independent of µ . (2.44) In Appendix B, we give the full expression for II (q) when r = 1.
Quite nicely, comparing the sum of (2.44) and (2.43) to (2.23) we find the expected log µ dependence for the functional determinant to this order. So far, we have kept the radii of S 2 and H 2 different in order to provide a detailed check of our results using the trace anomaly (2.23). From now on we will set R S = R H = 1, hence r = 1. 16 Using the expressions in Appendix B, we can write explicitly the quantity f (q) 1-loop appearing in (1.6) as where the second line can be taken as the definition of f(q). The reason for the slightly awkward definition will become clear in eq. (2.62) below. The quantity f(q) is evaluated numerically in Table 1 for a few values of q, and its large q behavior is derived in Appendix C.

One loop contribution at O( )
The computation of the previous section can be easily extended to higher orders in . Indeed, the functional determinant can be formally written down in any fractional dimension; the only difference from the = 0 case is that now the density of states for spinors on H 2− and n = 0 is [44]: (2.46) 16 The common radius R S = R H = R can always be easily reintroduced using dimensional analysis. By not dropping the r-dependence, one could have also extracted some information on integrated n-point functions of the stress tensor in the presence of the defect.
which is a simple generalization of (2.29). 17 Then the functional determinant (2.36) for finite is where µ is the UV cutoff.
Expanding the H 2− density of states (2.46) to linear order in , we get where ψ(z) is the digamma function. We see that the second term in the parenthesis gives a contribution to the free energy proportional to the leading term: The contribution from the first term in the parenthesis in (2.48) requires a more detailed analysis. Using experience from Section 2.4, to determine the q-dependent terms that are proportional to log µ, we only need to consider the asymptotic behavior (to order 1/λ 3 ) of the spectral density. In analogy with term I in (2.41) we write: where the first term is the contribution of = q. As in Section 2.4 we expand for large and use zeta function regularization for the resulting terms. The remaining finite sum does not give terms proportional to log µ, so we will not consider it here. At the end of the day, we 17 It is easy to check that this is the correct expression, as for = 0 it gives back (2.35), whereas for = 1 it gives µ(λ) = 1 π , which is the appropriate result for three dimensions. In three dimensions we usually write x, µ have where the linear in q term from the first term in (2.50) cancels.
Adding up (2.50) and (2.49) and copying over the 0 result in (2.45), we obtain an improved version of (2.45): (2.52) In principle, one can straightforwardly go to higher orders in and also compute the terms independent of µ.

Two loop contribution at O( 0 )
Having determined the first two coefficients in (1.8) in (2.14) and (2.52) respectively, we can now move on to the two-loop term. This is the contribution of the current-current correlator to the free energy appearing in the last line of (2.12); it can be thought of as the vacuum diagram shown in Figure 1. We will determine its singular piece which is proportional to log µ. In Appendix D we show that the short distance behavior of this diagram is determined by the flat space limit where i, j are frame indices, p is the flat space momentum, and we give formulas for The p i p j term doesn't contribute because K ij is conserved. The first two terms in (D.13) don't contribute if we use symmetric integration over the momenta. We're left with the last term: N e 2 0 q 2 32π 3 log µ . (2.54) We leave the computation of the finite part of D (q) to future work.
From (2.54), we deduce that the coefficient f By tuning the bare coupling e 0 to the weakly coupled conformal fixed point, we should be able to get rid of the µ-dependence of (2.56), so that we get an unambiguous answer for the free energy.
We will go through the procedure of determining e 0 at the fixed point in some detail to highlight the nontrivial cancellations that occur in the process. The reason why it made sense to further expand each order in the loop expansion in is because the theory is weakly coupled, e 2 0 ∼ . Hence we look for e 0 in the following form: Requiring that the log µ dependence cancels at O( 0 ) fixes With the value of α given in (2.58) this term also cancels. One of the main motivation for the calculations presented in Section 2.5 was to see this cancellation. From matching the log µ terms at O( ) we get that We conclude that where in the second line we resummed part of c 1 into a prefactor familiar from dimensional regularization. 18 Using the fixed point value of e 0 , (2.56) becomes: This concludes the calculation in the expansion. Before turning our attention to two consistency checks in the large q limit and in the large N limit of (2.62), respectively, we digress to derive the fixed point coupling constant (2.60) from another point of view. 18 In dimensional regularization logarithmic divergent integrals give 1/ poles that are accompanied by the prefactor in (2.60): log. div. int. = 1 + 1 2 log µ 2 4πe −γ + finite . (2.61)

More conventional derivation of the fixed point coupling
The expression (2.60) for the fixed point coupling can also be understood from a more conventional point of view. Note that by computing the free energy in a flux background we have essentially performed the computation of the effective potential of the gauge field. Gauge invariance dictates that the effective potential can only depend on q. We have computed the one-loop contribution exactly and have determined the divergent piece of the two-loop term.
It is then standard procedure to read off the two-loop beta function from such a calculation.
Because our zeta function scheme is such that it doesn't produce any divergences at any order of perturbation theory (unlike dimensional regularization), we don't get 1/s poles at the end of the calculation, and in our scheme there is a very simple relation between the bare and renormalized coupling (

The large q limit
The first check of our final expression (2.62) comes from its large q behavior, which was also mentioned in the introduction. Following the techniques introduced in [41], we calculate the large q behavior of the functional determinant in Appendix C to order q −1 . The leading order result is: where A is the Glaisher constant. From (2.62), we then obtain the large q behavior of the free energy, which can be written as In Figure 2 we compare the analytic large q expansion to order O(q −1 ) determinant (C. 4) to the exact numerical expression. In this comparison we subtract the q = 0 part of the exact numerical free energy from the analytic large q expansion so that the comparison is accurate.
Note that the plot goes like q −1 and asymptotes to zero, as we would expect for a large q approximation to order O(q −1 ). The large q expansion seems to be extremely accurate even for small q. We see that even for the smallest possible value of q = 1/2, the error is only O(10 −4 ).

The large N expansion
As already mentioned, another perturbative approach to the theories of interest is the large N expansion. In non-supersymmetric theories, this method is pretty much the only way monopole operators in 3d gauge theories have been studied thus far. The large N expansion can be performed in any space-time dimension d = 4 − . As mentioned in the introduction, it has the structure given in (1.10), which we reproduce here for the reader's convenience: In this section, we will determine g (q) 0 ( ) and compare the → 0 limit of the result to the 4 − expansion result of the previous section.
An interesting subtlety is that, when > 0, the Maxwell term is irrelevant, and hence it shouldn't be included in the Lagrangian if we wish to describe the IR CFT. Consequently, in contrast to the 4 − expansion method, in the large N approach we do not have a tunable gauge coupling e 0 ; instead, the theory is automatically conformally invariant. It is thus quite non-trivial to see how the results from the 4 − expansion agree with those from the large N expansion in their overlapping regime of validity.
Even though there is no Maxwell term in the action, the gauge field nevertheless acquires a non-local kinetic term from quantum effects due to the massless matter fields; since at large N there are many matter fields, the gauge field fluctuations are suppressed and the theory is perturbative. The leading contribution to the free energy comes from the matter field fluctuations around the monopole background. At subleading order, the Gaussian gauge field fluctuations contribute to give: where / ∇ 4− is the background gauge covariant derivative and K is (up to contact terms) the current-current correlation function of the CFT (2.10), which becomes the kernel of the gauge field fluctuations. By comparing (3.1) with (3.2), we see that g (q) 0 ( ) is determined from the functional determinant, and g (q) 1 ( ) is calculable from K. We put tr log K in quotation marks in (3.2) to indicate that there are subtleties related to gauge invariance that will not be important here. Similar calculations were performed in three-dimensional theories with fermionic matter in [18,19] and in bosonic theories [41,46]. By contrast, in the large N method we first take s → 0 and then, if we wish to describe the theory close to four dimensions, we can take → 0. Indeed, expanding the regularized functional determinant at small and comparing to the relevant terms from the expansion method will allow us to check the results of the expansion method. In some sense, the large N results are interesting in that they provide an orthogonal viewpoint on our computation.

Functional determinant at finite
The generalization of the free energy (2.36) to finite is given in (2.47) in the case where the curvature radii are R S = R H = 1. For general R S and R H , it takes the form where the density of states µ 2− (λ, 0) was given in (2.46) and the degeneracy of the modes on S 2 was given in (2.27).
The method for computing (3.3) is as in Section 2.4. In particular, we first divide the integral into the large λ divergent part and the remaining convergent part, thus writing with Here, is just an approximation of µ 2− (λ) at large λ, and s is a regularization parameter. Then, we use zeta function regularization to regularize the divergent sums and integrals. (Note that now that is nonzero, term I (q) would have been infrared divergent had the lower limit of integration been 0; to avoid this problem, we divided the integration range into The final expressions are rather complicated, and we will not reproduce them here. We plot the final result in Figure 3 as a function of .
Let us now explain how the → 0 limit of this large N result is compatible with the small expansion previously obtained. This limit should be taken carefully, because the → 0 limit does not commute with the s → 0 limit that is implicit in the definition of the linear functional D. Actually, the two limits do not commute only for I (q) ; they commute for II (q) , and thus  Setting R S = R H = 1 from now on, we find 21 (3.9) In (2.45), we found Adding up (3.9)-(3.10), we obtain which at large N agrees precisely with the expression (2.62) we obtained in the expansion! It is quite remarkable how the fermion determinant computed in the order of limits s → 0 then → 0, reproduces, for instance, the free energy contribution from the Maxwell action evaluated on the classical background when the coupling constant e 0 is tuned to its critical value. 22 For > 0, there was no coupling constant that needed to be tuned to criticality, and yet the result is continuously connected to that in the 4 − expansion where such a coupling was explicit in the action. We therefore regard the computation in this section as a highly non-trivial check of the approach presented in the previous section.
In the Introduction we argued that from the consistency between the large N and expansions, the coefficients g 1 should follow from a similar computation of the gauge field fluctuations, where again the order in which we take → 0 and remove the regulator should make a difference. 21 More generally, we expect that the difference between I (q) − I (0) and the corresponding quantity computed by expanding in first and then regularizing resums to N q 2 (4πe −γ ) /2 6πµ . In (3.9) we see only the first two terms in the small expansion of this quantity. 22 It is worth noting that terms such as the log(4πe −γ ) appearing (3.9) were also encountered in [47] when comparing the results of zeta function and dimensional regularizations on product manifolds.

Remarks on the flavor symmetry transformation of defect operators
As discussed in detail in [48], It is thus natural to ask if the defect operators studied in this paper transform under the flavor symmetry.
In d = 3 it is known how the monopole operators transform under the flavor symmetry [16,19]. The Dirac operator on S 2 in the presence of monopole flux has 2 |q| zero modes. The eigenfunctions of the Dirac operator on hyperbolic space were computed in [49]. We

A.2 Spinor monopole harmonics
The spinor monopole spherical harmonics can be described using the usual scalar monopole spherical harmonics q Y m (θ, φ), as explained in [16,18]. A minor difficulty is that in these references, the spinor harmonics were written in a basis for the spin bundle which is the result of a conformal transformation from R 3 to R × S 2 . If τ parametrizes R, the frame used in [16,18] can be written as e 1 = cos θ cos φdθ − sin θ sin φdφ + sin θ cos φdτ , e 2 = cos θ sin φdθ + sin θ cos φdφ + sin θ sin φdτ , e 3 = − sin θdθ + cos θdτ . and then simply erase the R direction along with e 3 . To match our conventions, we should also take γ i = σ i , as was also done in [16,18].
The two frames (A.5) and (A.6) differ by a frame rotation. If we write In the conventions in the main text, the basis described in [18] is then (A.10) On the basis spinors, the gauge-covariant Dirac operator in a monopole background acts as such that the eigenvector/eigenvalue combinations are Term II is of the form where the bullet stands for the integrand. We will therefore list only the integrands, denoted by lower case Roman numerals, which all give finite results when integrated in (B.2): ii (q) = − 2ζ (−1, q + 1) + q 2 − λ 2 ψ(q + 1) + q log λ In terms of these quantities, we express f(q) as C O( ) functional determinant at large q In the limit q → ∞ we can evaluate the functional determinant in Section 2.4 exactly. We will evaluate this determinant to order O(q −1 ). In the following we set R S , R H , µ = 1, for simplicity.

(D.4)
One can check that i / ∇G = 0 at separated points. It is normalized so that it gives a deltafunction of unit strength.

D.2 Fermion propagator when q = 0
When q = 0, one can find a series expansion for the Green's function. To simplify the expressions, let's define as well as which is the quantity that appears in the denominator of the q = 0 Green's function (D.4).
Then expanding at small S, we have where c is an arbitrary constant and G analytic is an analytic piece. (The constant c multiplies analytic terms also.) This expression was obtained from requiring G q to satisfy i / ∇G = 0 at separated points.