SISSA 40/2014/FISI
Bounds on OPE Coefficients in 4D Conformal Field Theories
Francesco Caracciolo, Alejandro Castedo Echeverri,
Benedict von Harling and Marco Serone
SISSA and INFN, Via Bonomea 265, I34136 Trieste, Italy
ICTP, Strada Costiera 11, I34151 Trieste, Italy
Abstract
We numerically study the crossing symmetry constraints in 4D CFTs, using previously introduced algorithms based on semidefinite programming. We study bounds on OPE coefficients of tensor operators as a function of their scaling dimension and extend previous studies of bounds on OPE coefficients of conserved vector currents to the product groups SOSO. We also analyze the bounds on the OPE coefficients of the conserved vector currents associated with the groups SO, SU and SOSO under the assumption that in the singlet channel no scalar operator has dimension less than four, namely that the CFT has no relevant deformations. This is motivated by applications in the context of composite Higgs models, where the strongly coupled sector is assumed to be a spontaneously broken CFT with a global symmetry.
To the memory of Francesco who tragically passed away during the completion of this project.
Contents
1 Introduction
There has recently been a renewed interest in studying general properties of fourdimensional (4D) Conformal Field Theories (CFTs) after the seminal paper [2] revived the bootstrap program advocated in the early 70s [3, 4]. Imposing the associativity of the Operator Product Expansion (OPE) and unitarity, ref.[2] has shown how one can set bounds on scalar operator dimensions in 4D CFTs. Although these constraints are based on numerical methods, they come from first principles, with no further assumptions. Since then, various generalizations of this result have been developed in order to improve the above bounds, to put bounds on OPE coefficients and on CFT data in presence of a global symmetry [5, 6, 7, 8, 9, 10, 11]. In particular, bounds were derived on the OPE coefficients associated with the energy momentum tensor (the central charge) and the conserved vector current of a global symmetry [7, 8, 9, 10, 11]. Superconformal field theories and CFTs in have also been considered [7, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22]. In addition to the above mostly numerical results, considerable progress has also been made on more analytic aspects of CFTs, see e.g. refs.[23, 24, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 25, 26].
The aim of this paper is to numerically study the bounds on the coefficient of the twopoint function between two conserved currents associated with a global symmetry of a CFT. Our main motivation comes from theoretical considerations in the context of composite Higgs models, in which the CFT is the hidden sector which gives rise to the Higgs, and a subgroup of the global symmetry of the CFT is weakly gauged in order to get the Standard Model gauge interactions. These composite Higgs models models are related, through the AdS/CFT correspondence, to RandallSundrum theories [37] with matter in the bulk, which are a promising solution to the gauge hierarchy problem. Particularly interesting are the models where the Higgs is a pseudo NambuGoldstone Boson (pNGB) of an approximate spontaneously broken global symmetry of the CFT, which correspond to gaugeHiggs unification models in 5D warped theories. Neither the UV completion of the 5D models nor the explicit form of the 4D CFT is known so far. Calculability of the dual 5D models would require that the CFT is in some large N limit, but this is not a necessary requirement. On the contrary, various phenomenological bounds tend to favour models at small N, so we will not assume the existence of a large N limit in the CFT. Constructing such a CFT is not a trivial task, so we look for possible consistency relations. When the global symmetry of the CFT is gauged, the coefficient of the currentcurrent twopoint function governs the leading contribution of the CFT to the oneloop evolution of the corresponding gauge coupling. This contribution should not lead to Landau poles for the SM gauge couplings. We also require that the CFT has no relevant deformations, in order not to reintroduce the hierarchy problem. This leads to the constraint that the dimension of the lowestlying scalar singlet operator should be . All our considerations apply independently of the pNGB nature or not of the Higgs.
Motivated by the above considerations, we extend the analysis of ref.[11], where lower bounds on have been set starting from crossing constraints imposed on a fourpoint function of scalar operators in the fundamental representation of SO and SU, in two ways. First, we see how the bounds found in ref.[11] are modified when the lowestlying singlet scalar operator is assumed to have a scaling dimension , where we choose for concreteness. Second, we extend the analysis to nonsimple groups of the form SOSO. We study nonsimple groups because they easily allow to generalize the bounds for the groups SO and SU, which are obtained by considering a single field in the fundamental representation of the group, to multiple fields. Analogous to what was found in ref.[11] for singlet operators, the lower bounds on vector currents for SU groups that we find are, within the numerical precision, identical to those obtained for SO. Hence we only report lower bounds for SO and SOSO global symmetries. We have derived the bootstrap equations also for groups of the form SOSU, but no bounds are reported for this case, since we have numerical evidence that the lower bounds for SOSU are essentially identical to those obtained for SOSO, similarly to the above equality between SO and SU bounds. In addition to that, we study the constraints on the OPE coefficients of spin and tensors coming from two identical scalar operators , as a function of the scaling dimension of the tensors, in the general case in which no global symmetry is assumed. In analogy to the vectorcurrent case, we analyze how these bounds change when one assumes a lower bound on the scale dimension of the scalar operators appearing in the OPE.
All our numerical results are based on semidefinite programming methods, as introduced in ref.[11] in the context of the bootstrap approach, with a few technical modifications which are discussed in subsection 3.2 and in appendix A.
The structure of the paper is as follows. In section 2 we describe the phenomenological motivations behind our work. In section 3 we briefly review the basic properties of the crossing constraints coming from fourpoint functions of identical scalars and review how bounds on OPE coefficients are numerically obtained. In section 4 we report our results for the OPE coefficients of tensor and operators. Section 5 contains the most important results of the paper. We report here the lower bounds on associated with SO (or SU) vector currents, when the global symmetry of the CFT is SO (or SU) and . ^{1}^{1}1Results for SO groups with odd are analogous to those for SO and do not need any special treatment. In section 6 we conclude. Two appendices complete the paper. In appendix A we discuss various technical details about our implementation of the bootstrap equations in the semidefinite programming method, while in appendix B we report the crossing equations for SOSO and SOSU.
2 Motivation for a Gap in the Scalar Operator Dimension
The motivation to consider CFTs with a gap in the scaling dimension of scalar gaugesinglet operators comes from applications in the context of physics beyond the Standard Model (SM) that addresses the gauge hierarchy problem. The latter can be formulated from a CFT point of view, see e.g. ref.[2]. Neglecting the cosmological constant, the SM can be seen as an approximate CFT with one relevant deformation of classical mass dimension , corresponding to the Higgs mass term . Relevant deformations grow in going from the UV towards the IR. If we assume that the Higgs mass term is generated at some high scale , we would expect from naturalness that the Higgs masssquared term is of order in the IR. There are essentially two ways to solve this hierarchy problem: i) invoke additional symmetries that keep the relevant deformation small in the IR (e.g. supersymmetry); ii) assume that the Higgs is a composite field of a strongly interacting sector, in which case the operator can have a large anomalous dimension that makes it effectively marginal or irrelevant.
A model along the lines of ii), conformal technicolor [38], where the strongly coupled sector is assumed to be a CFT in the UV, was in fact the motivation for the pioneering work [2]. Conformal technicolor is an interesting attempt to solve one of the longstanding problems of standard technicolor theories: how to reconcile the top mass with Flavour Changing Neutral Current (FCNC) bounds. In order to get a sizable top mass and at the same time avoid dangerous FCNCs, one has to demand that the scale dimension of the Higgs field is as close to one as possible. In order not to reintroduce the hierarchy problem, however, one has to keep at the same time. The analyses in refs.[2, 11] have shown that generally these two conditions are in tension and that one needs in order to have .
An alternative, phenomenologically more promising, solution is to rely on a different mechanism to generate SM fermion masses: partial compositeness [39]. To this end, one assumes that the SM fermions mix with fermion resonances of the strongly coupled sector. Due to this mixing, SM vectors and fermions become partially composite. In particular, the lighter the SM fermions are, the weaker is the mixing. This simple, yet remarkable, observation allows to significantly alleviate most flavor bounds. The Yukawa couplings are effective couplings that arise from the mixing terms once the strongly coupled states are integrated out.
This idea is particularly appealing when one assumes that the strongly coupled sector is an approximate CFT spontaneously broken at some scale . In this case, the hierarchy of the SM Yukawa couplings is naturally obtained by assigning different scale dimensions to the fermion operators mixing with the different SM fermions [40]. In particular, there is no longer the need to keep close to one since the effective size of the SM Yukawa couplings is governed by . One assumes that for all SM fermions except the top, so that the mixing terms are irrelevant deformations of the CFT and naturally give rise to suppressed Yukawa couplings in the IR. For the top, on the other hand, one assumes that , corresponding to a nearly marginal deformation of the CFT.
One might wonder whether CFTs with all the necessary requirements to give rise to theoretically and phenomenologically viable composite Higgs models exist at all. A possible issue might arise in weakly gauging the SM subgroup of the global symmetry of the CFT. Since partial compositeness requires a fermion operator in the CFT for each SM fermion, dangerous Landau poles can potentially appear in the theory. Indeed, it has recently been shown that Landau poles represent the main obstruction in obtaining UV completions of composite Higgs models with a pNGB Higgs, based on supersymmetry [41]. It is then of primary importance to try to understand if and at what scale Landau poles will arise. In theories with a pNGB Higgs, the relevant deformation can naturally be small, since it is protected by a shift symmetry. Moreover, it is not defined in the UV, where the global symmetry is restored. Nevertheless, in order not to introduce other possible finetunings, one should demand that any scalar operator which is not protected by any symmetry, namely which is neutral under all possible global symmetries of the CFT, should be marginal or irrelevant.
Summarizing, we can identify four properties that a CFT needs to have for a theoretically and phenomenologically viable composite Higgs model with partial compositeness:

A global symmetry .

No scalar operator with dimension which is neutral under .

No Landau poles for the SM gauge couplings below the scale when we gauge .^{2}^{2}2Ideally, we might want to have .

The presence of fermion operators with in some representation of , such that some of its components can mix with each of the SM fermion fields. At least one fermion operator should have dimension .^{3}^{3}3The righthanded top, in principle, might be directly identified with a field of the CFT.
Of course, these are only necessary but not sufficient conditions to get a viable CFT. In particular, one might want to address the mechanism which gives rise to the spontaneous breaking of the conformal symmetry as well as of the global symmetry in CFTs with a pNGB Higgs.
The consistency of a CFT which fulfils the above four conditions can be checked using crossing symmetry of fourpoint functions of the CFT. The first and second condition can be imposed by hand, assuming the existence of the global symmetry and that the lowestdimensional scalar operator in the singlet channel has dimension . One can extract information on the third condition by analyzing the bounds on the coefficients of currentcurrent twopoint functions. Finally, the fourth condition can again be implemented by assumption. The ideal configuration would be to analyze fourpoint functions involving fermion operators, which by assumption should appear in the CFT, and to extract any possible information from these correlators. Although this is in principle possible to do, correlation functions involving fermions in a nonsupersymmetric setting have not been worked out so far. Postponing to a future project the analysis of such correlation functions, in this paper we start to address these issues by replacing fermions with scalars with dimension in the third requirement.
Let us estimate how severe the Landau pole problem can be in the simplest composite Higgs model where the Higgs is the pNGB associated with the SO(5)SO symmetry breaking pattern. Let us consider the SU coupling , because it runs fastest and possibly leads to the lowestlying Landau pole, and let us denote by
(2.1) 
the CFT contribution to its oneloop function. Assuming that the only nonSM fields which are charged under SU arise from the CFT, a Landau pole develops at around
(2.2) 
for , where and is the scale where the CFT breaks spontaneously. Composite fermions coming from the CFT and mixing with SM fermions must be color triplets and in representations of SO that give rise to electroweak SU doublets and singlets. If we assume them to be in the fundamental representation of SO, the fermion components in a given can mix with both the lefthanded and righthanded components of a quark field. We then need s, one for each quark field, for a total of SU triplet Dirac fermions. In order to have an idea of the scales which are involved, it is useful to consider the (unrealistic) limit of a free CFT. In this case, we get
(2.3) 
corresponding to
(2.4) 
for TeV. It is clearly very important to set lower bounds on in a generic CFT, given the exponential sensitivity of on this quantity.
In the following, we will analyze bounds on the coefficients for SO (or, equivalently, SU) currents obtained from fourpoint functions of scalar operators in the fundamental representation of the group in presence of a gap in the operator dimension in the scalar gaugesinglet channel. In order to mimic the presence of more than one field multiplet, we will also consider fields in the bifundamental representation of the product group SOSO.
3 Review of the Bootstrap Program
In this section, we briefly review the equations that one obtains by imposing crossing symmetry on fourpoint functions of scalar operators in a unitary CFT, and how these are numerically handled to get bounds on the CFT data. We refer the reader to ref.[2] and references therein for more details and background material.
3.1 Bootstrap Equations
The fourpoint function of identical real scalar operators with scale dimension in a 4D CFT can be written as
(3.1) 
where and
(3.2) 
are two conformally invariant variables. All the dynamics is encoded in the function . Using the OPE between and (the channel), in the region this function is found to be
(3.3) 
where the sum is over all primary (traceless symmetric) operators of dimension and (even) spin that appear in the OPE and is the coefficient of the threepoint function . The in eq. (3.3) results from the contribution of the identity operator which is present in the OPE of two identical operators. The threepoint function in this case simplifies to the twopoint function which is normalized to unity. For each primary operator , the function , which is called a conformal block, takes into account the contribution of all descendants of . In 4D, the explicit form of is known and reads [42, 43]
(3.4) 
where and .^{4}^{4}4Here we have used the normalization of the conformal blocks introduced in ref.[9] which differs by a factor from the one of refs.[42, 43]. Alternatively, we can obtain an expression for using the OPE between and (the channel). Demanding that the channel and channel results for the fourpoint function agree gives the crossing symmetry constraint (or bootstrap equation)
(3.5) 
with
(3.6) 
and . When the CFT has a global symmetry , the above analysis can be generalized using scalar fields (real or complex) in some representation of [9]. The symmetry implies that all the field components of the multiplet must have the same dimension . Moreover, it allows to easily classify the operators appearing in the OPE in terms of the irreducible representations appearing in the product . A similar analysis applies for complex fields in the OPE. It is useful to introduce another function, similar to the of eq. (3.6):
(3.7) 
In presence of a global symmetry , eq. (3.5) generalizes to a system of equations of the form
(3.8) 
Here, runs over all possible irreducible representations that can appear in the  and channel decomposition, and are numerical factors that depend on and is a shorthand notation for the threepoint function coefficient. Furthermore, and if the singlet representation appears in the lefthand side of eq. (3.8), and otherwise. The explicit form of eq. (3.8) for the cases of interest will be given in section 5 and appendix B.
3.2 Bounds on OPE Coefficients and Numerical Implementation
The bootstrap equation (3.5) has originally been used to set bounds on the scalar operator dimensions that can appear in a CFT. Shortly after that, ref.[6] has shown how to obtain bounds on the OPE coefficient of an operator appearing in the OPE. Let us assume that a linear functional can be found, such that
(3.9) 
Applying such a functional to eq. (3.5) gives
(3.10) 
The optimal bound is obtained by minimizing among all the functionals which satisfy eq. (3.9). One can use the functional also to rule out the existence of certain CFTs. For instance, if under a certain assumption on the CFT data one finds a functional and an operator for which , then that CFT is ruled out.
The above procedure is easily generalized in presence of global symmetries. Let us assume that we want to bound the OPE coefficient of an operator with dimension and spin in the representation . We look for a set of linear functionals () such that
(3.11) 
Applying such a functional to eq. (3.8) gives
(3.12) 
In our paper, we will mainly be interested in the OPE coefficient associated with a conserved vector current of a global symmetry, which has and . We shall denote this coefficient by . As we will discuss in section 5 (see eq. (5.9)), upper bounds on turn into lower bounds on the coefficient introduced in eq. (2.1).
Following ref.[2], we consider functionals that act as linear combinations of derivatives on a generic function ,
(3.13) 
where are real coefficients. Due to the symmetries of the conformal blocks and , the sum can be restricted to and even values of when acts on , and and odd values of when it acts on .
We numerically search for functionals which satisfy eqs. (3.9) and (3.11) by following the method developed in refs.[7, 11]. We refer the reader to these references for further details. For this method, one approximates the derivatives of the conformal blocks and in eq. (3.13) with polynomials , where and is the unitarity bound on the scaling dimension for an operator of spin ( for ). The requirements in eq. (3.9) or eq. (3.11) imply that the linear combination of polynomials of the form must be positivesemidefinite on the positive real axis, for any value of . There are two great virtues in setting up the problem in this way. Firstly, there is no need to discretize the dimension and to put a cutoff value , like in the linear programming methods used in ref.[2]. In particular, we can probe all continuously up to infinity. Secondly, one can exploit numerical packages that allow to handle very large systems of equations quite efficiently. A key variable in the numerical algorithm is the coefficient entering in eq. (3.13). The larger , the larger is the space of possible viable functionals, and hence the stronger are the bounds. Of course, the larger , the more timeconsuming is the numerical evaluation. For our computations, we have chosen , depending on the complication of the problem.
The above algorithm, however, still requires to truncate the system at a given maximal spin . This is in principle a serious problem, because one might have
(3.14) 
If is chosen sufficiently large, or more, we do not expect possible violations in the semidefinite positiveness of of the form (3.14) to be important for the numerical value of the bound. Indeed, large spin implies large dimensions according to the unitarity bound, and the contribution to the fourpoint function of operators with large is exponentially suppressed in [30]. Nevertheless, it would be more reassuring to have more control on such effects. For parametrically large , the conformal blocks and and their derivatives allow for simple analytic expressions. For large , the terms involving the highest derivatives dominate. Using these analytic expressions, we can find the value , which depends on , for which the contribution of the large conformal blocks is largest. We find, for (see appendix A for details)
(3.15) 
where . For , eq. (3.15) gives . Ideally, one would include all spins from up to . This is computationally quite demanding. Fortunately, we have found that it is sufficient to take to get numerically stable bounds. Changing to or does not significantly alter the bounds. Nevertheless, in order to have more control on the higher states, we have included two other states in the constraints, at and at an intermediate value .^{5}^{5}5More precisely, we include spins for calculations with , for and for . We have numerically tested that this implementation works better than including states at very large values of , such as as done in e.g. ref.[11]. We can always check the positivity of a posteriori. We have found that by imposing constraints at the functional often becomes negative for values whereas for our implementation remains positive for most of the that we have checked. In practice, however, we have not detected deviations in the results among the two different implementations, confirming that values of are numerically negligible.
4 Bounds on OPE Coefficients for Tensor Operators
In this section, we report our results for the upper bounds on the threepoint function coefficient appearing in the OPE of two identical scalar operators of scaling dimension . The operator is a traceless symmetric tensor of even spin . The coefficient is normalized such that its freetheory value is
(4.1) 
We do not report the results for the case, which were first derived in ref.[6] and subsequently improved in ref.[11]. Our results agree with fig. 10 of ref.[11]. These bounds change if we assume that the first scalar operator which appears in the OPE has a dimension , where is the unitarity bound on . As expected, the upper bounds do not significantly change when is close to 1, since by continuity the theory is close to the free theory, where the only scalar operator arises exactly at . For values of not too close to 1, on the other hand, the bound is significantly improved and becomes more stringent as increases. In fig. 1, we report the bounds for and .
Analogously, one can study the upper bounds on for generic tensor operators with spin and dimension . Upper bounds on the central charge associated with the energymomentum tensor (the lowestdimensional operator in the sector, with ), have been extensively analyzed in refs.[7, 8, 11], with and without the assumption of a lower bound on the dimension of the lowestlying scalar operator appearing in the OPE. In fig. 2 (a), we report the upper bounds on the coupling between two scalar operators of dimension and a tensor operator with spin and dimension for different and as a function of . As can be seen, the larger is, the less stringent is the upper bound, in agreement with the naive expectation for which can be seen as a measure (for not too far from 1) of how strongly coupled the CFT is. Like for scalar operators, the bounds change if we make some assumptions on the CFT spectrum. As for the scalar case, the upper bounds do not significantly change when is very close to 1, but for values of not too close to 1, they become more stringent as increases. For illustration, in fig. 2 (b), we report the upper bounds on as a function of for , assuming that the lowest scalar operator appearing in the OPE has a dimension , and .^{6}^{6}6The value is roughly the minimal one compatible with the assumption , see e.g. fig. 2 of ref.[11].
Similarly, one can analyze tensor operators at higher . In figs. 3 (a) and (b), we report the same as above for operators. As expected, the absolute scale of becomes lower and lower as increases, with the allowed values of quickly decreasing as becomes larger. Notice that the maximal allowed value of both and is centered at values of that increase as is increased.
5 Bounds on CurrentCurrent TwoPoint Functions
At leading order, the CFT contribution to the oneloop beta function of a gauge field , external to the CFT, is governed by the coefficient of the two pointfunction of the corresponding current. Denoting by
(5.1) 
the total Lagrangian after the gauging, we can consider the effective action defined as (in euclidean signature)
(5.2) 
where the functional integration is over all the CFT states and we have omitted color indices. In general
(5.3) 
where and is the CFT contribution to the wave function renormalization of the gauge field, which in turns gives us the oneloop contribution of the CFT to the RG running of :
(5.4) 
By taking two functional derivatives with respect to and in eq. (5.2), we readily get
(5.5) 
where the subscript in the correlator specifies that the twopoint function is computed in the unperturbed CFT setting . The normalization of the current is uniquely fixed by Ward identities. Following the notation of ref.[11], we parametrize the twopoint function in configuration space as follows:^{7}^{7}7Notice that the definition of here is not identical to that of ref.[11] which tacitly applies to CFTs with one charged multiplet only. In general, where runs over all the charged fields of the CFT in the representations and .
(5.6) 
The “vector central charge” is roughly a measure of how many charged degrees of freedom are present in the CFT, similar to the standard central charge being a measure of the total number of degrees of freedom of the CFT. Modulo irrelevant contact terms, the momentum space correlation function reads
(5.7) 
and hence
(5.8) 
We extract by rescaling the vector current so that it appears as the coefficient of the threepoint function :
(5.9) 
Upper bounds on turn into lower bounds on . The constant factor is easily found by matching the result with the freetheory case, in which both and are calculable. In what follows, we will analyze the lower bounds on for different vector currents that come from the crossing symmetry constraints applied to fourpoint functions of scalars.
5.1 So Global Symmetry
We consider a fourpoint function of real scalars that are taken to be the components of a single field in the fundamental representation of SO with dimension . The crossing symmetry relations have been derived in ref.[9]. We report them here for completeness:
(5.10) 
Here , and refer respectively to the singlet, rank2 symmetric and antisymmetric (adjoint) representations of the operators which define the different conformal blocks. For the superscript , only even spins are included in the sum whereas for only odd spins are summed over. For simplicity, we have omitted the labels and the arguments of the conformal blocks and . Bounds on (as defined in eq. (5.6)) in this setup have already been found in ref.[11]. In this subsection we will see how these bounds change if assumptions on the dimensionality of the lowestlying scalar operator in the singlet channel are made.
First of all, let us consider the free theory of a real scalar in the fundamental representation of SO in order to fix the constant in eq. (5.9). The freetheory values of the OPE coefficients in the three different channels read
(5.11) 
where is given in eq. (4.1). We in particular get . Matching eq. (5.8) with the oneloop contribution to the function of a scalar in an SO gauge theory gives
(5.12) 
where we have taken (cf. footnote 7). From this it follows that in eq. (5.9).
In fig. 4, we report our results in terms of lower bounds on . We have considered the five different values and report the lower bounds on for the case where no assumption on the spectrum is made (the lines starting from ) and the case where the lowestlying scalar operator in the singlet channel is assumed to have dimension (the other lines). The former bounds agree with previous results (e.g. compare with fig. 18 of ref.[11]). Although it is not clearly visible from the figure, we have checked that all the bounds consistently tend to the freetheory value for . The latter bounds start from a given that depends on . This is of course expected, given the known results for the upper bound on the dimension of the lowestlying scalar singlet operator at a given : CFTs at are excluded under the assumption of a gap in the scalar singlet sector. The values of that we find agree with the values given in the literature (compare e.g. with the dimensions for which in fig. 4 of ref.[11]). The lower bounds on become significantly more stringent when we impose that . They also decrease less rapidly when increases compared to the unconstrained case.
In order to show how the assumption on affects the lower bounds on , in fig. 5, we fix and consider the three cases , and . As expected, the lower bound consistently becomes more severe as we increase . As before, the bounds start at certain dimensions which agree with previous results (compare e.g. with the dimensions for which in fig. 4 of ref.[11]).
5.2 Su Global Symmetry
We consider a fourpoint function of complex scalars that are taken to be the components of a field in the fundamental representation of SU with dimension . The crossing symmetry relations have been derived in ref.[9]. We report them here for completeness:
(5.13) 
Here , , and refer respectively to the singlet, adjoint, rank2 symmetric and rank2 antisymmetric representations of the operators which define the different conformal blocks. For the superscript , even spins are included in the sum, and for , odd spins are summed over. We consider here the lower bounds on (as defined in eq. (5.6)) associated with the adjoint current.
As in subsection 5.1, we start by looking at the free theory of a complex scalar in the fundamental representation of SU in order to fix the constant in eq. (5.9). The freetheory values of the OPE coefficients in the four different channels read
(5.14) 
where is given in eq. (4.1). We in particular get . Matching eq. (5.8) with the oneloop contribution to the function of a complex scalar in an SU gauge theory gives
(5.15) 
where we have taken (cf. footnote 7). From this it follows that in eq. (5.9) as for SO.
The six crossing symmetry equations (5.13) should reduce to the three equations (5.10) when the group SU is embedded in an underlying SO group. The decomposition of the singlet, adjoint and rank2 symmetric representations of SO in terms of SU representations reads
(5.16) 
If , for each primary operator in the () representation of SU, there is a corresponding operator in the and () representation as follows from eq. (5.16). The OPE coefficients of these operators are related by the underlying SO symmetry, , . It is straightforward to check with these identifications that eqs. (5.13) reduce to eqs. (5.10).
As we have already mentioned, the numerical results for the lower bounds on for SU are identical to those for SO, see fig. 4. This suggests that, given a set of three functionals that satisfy eq. (3.11) with , and as given by eq. (5.10), one should be able to construct a set of six functionals as linear combinations of the such that these functionals satisfy eq. (3.11) with , and as given by eq. (5.13). It would be interesting to find such a mapping and hence to understand in more analytical terms why the bounds on for SO( and SU() are equal.
5.3 Global Symmetries
The lower bounds on found in subsections 5.1 and 5.2 apply to CFTs in presence of at least one scalar field in the fundamental representation of , where or SU. Of course, the CFT can contain additional charged fields, for example a number of scalars in the fundamental representation of , with dimensions . In the freetheory limit of decoupled scalars (real for SO, complex for SU) we would simply have
(5.17) 
The larger is, the more constraining (and interesting) the lower bounds are. One cannot naively rescale the results of fig. 4 by a factor of in order to match the new free theory limit, however, because the interactions among the scalars will not be taken into account in this way. A more constraining bound could likely be obtained by studying the coupled set of fourpoint functions involving all scalars. This is in general not straightforward to do, since the crossing symmetry constraints are significantly more involved in presence of fields with different scaling dimensions.^{8}^{8}8See ref.[44] for very recent progress in this direction. A simple way to mimic the presence of more fields charged under a given group, though at the cost of assuming identical scaling dimensions , is obtained by introducing a further global symmetry group and assuming that the fields transform under some representation of . This is the main motivation for us to consider global symmetries which are direct products of two simple groups: it is a way to obtain lower bounds on in presence of more than one field charged under . More specifically, in the following we will consider fields in the fundamental representation of .
5.3.1 SoSo
Consider a CFT with global symmetry SOSO and one real scalar in the bifundamental representation of SOSO, and being SO and SO indices, respectively. In complete analogy to the SO case discussed in ref.[9], we can impose crossing symmetry in the  and channel on the fourpoint function in order to obtain the bootstrap equations. The operators appearing in the OPE transform under SOSO according to the decomposition of , where N and M denote the fundamental representations of respectively SO and SO. This gives 9 different representations, consisting of pairs , where refer to the singlet (), symmetric () and antisymmetric () representations of respectively SO and SO. Correspondingly, we get a total of equations. We report them in eq. (B.2) in appendix B.1.
The SO conserved current that we analyze is in the representation and is the lowestdimensional operator appearing in the functions and defined in eq. (B.1). In fig. 6, we show the lower bounds on for the three cases . While for the case of SO or SU considered before the lower bound first becomes significantly more stringent with growing and only from a certain onwards becomes less stringent, here a slight increase in the bound arises only for very close to 1 after which the bound decreases with . The lines starting from correspond to the case where no assumption on the spectrum has been made, while for the other lines we have assumed that the lowestlying scalar operator in the SS channel has dimension . The lower bounds for the latter case are stronger than for the former, but the difference is less substantial than for the groups SO or SU. The lower bound on the dimension of above which the lowestlying operator in the channel can have a dimension is also weaker than what was found for SO or SU. This is expected, since this bound becomes the weaker the larger the group is. The correct freetheory limit is obtained in all three cases. The shape of the lower bound on with no assumption on the spectrum in fig. 6 resembles the bound found in ref.[11] for SU singlet currents (see e.g. their fig. 19). From the SO point of view, the SO current is in fact a collection of singlet currents. On the other hand, for , we find that the lower bound on for SO currents shows the characteristic bump of single SO or SU currents, well above the freetheory value, as in fig. 4 (a). For illustration, we show the bound on obtained for and in fig. 7. It would be interesting to further explore these bounds and to understand the origin of their different behaviours in the regimes and .
As a final application of our results, in fig. 8, we report the lower bounds on for the group SOSO. We choose SO because the contribution of 120 free complex scalar triplets to the currentcurrent twopoint function gives . This in turn is the same value found in eq. (2.3) for the number of free fermion triplets which are needed to give mass to all the SM quarks in the SO(5)SO pNGB composite Higgs model mentioned in section 2. We consider SO(6)SO and not SUSO because the latter case is computationally very demanding (incidentally, in one of the models presented in ref.[41], SU was actually embedded in an underlying SO flavour global symmetry). Anyhow, given the equivalence between the SO and SU lower bounds on , we believe that these results would also hold for the SU(3)SO case.
As we see in fig. 8, assuming the absence of a relevant scalar singlet operator in the CFT does not significantly change the bounds. Furthermore, the most dangerous region regarding Landau poles, which is the region close to , is not consistent with the assumption of absence of relevant deformations. If we demand that no subPlanckian Landau pole arises, then we need , while for , remains asymptotically free.
5.3.2 SoSu
Consider a CFT with global symmetry SOSU and one complex scalar in the bifundamental representation of SOSU, and being SO and SU indices, respectively. We impose crossing symmetry in the  and channel on the fourpoint function and the fourpoint function with . The operators appearing in the OPE transform under SOSU in representations , where refer to the singlet (), symmetric () and antisymmetric () representations of SO and refer to the symmetric () and antisymmetric () representations of SU. This gives 6 different representations, with even and/or odd spin operators, depending on the representation. The operators appearing in the OPE transform in representations , with as before, whereas refer to the singlet () and adjoint () representations of SU. Taking into account that SU singlet and adjoint operators appear with both even and odd spins, we get 12 different conformal blocks, for a total of 18 bootstrap equations. We report them in eqs. (B.5) – (B.7) in appendix B.2. The SO conserved current that we are interested in transforms under the representation and is the lowestdimensional operator appearing in the functions and defined in eq. (B.4).
The large number of bootstrap equations makes the numerical analysis of the SOSU case computationally very demanding. We have however been able to compute the lower bounds on at level with , using and and either assuming that no scalar operator has dimension or without any assumption. The results are essentially indistinguishable from those obtained for when considering SOSO, though the accuracy obtained is not enough to claim that they are identical. Given also the observed equivalence of the bounds for SU and SO, we take this result as enough evidence to conjecture that the bounds for SOSU are identical to those for SOSO, for any and .
6 Conclusions
In this paper, we have numerically studied bounds on various OPE coefficients in 4D CFTs using the bootstrap approach applied to fourpoint functions of scalar operators. We have first studied bounds on OPE coefficients of symmetric traceless tensor operators with spins and as a function of their scaling dimension. Furthermore, we have analyzed how an assumption on the dimension of the lowestlying scalar operator affects such bounds.
We have then considered 4D CFTs with a global symmetry . When this group, or a subgroup of it, is gauged by weakly coupling external gauge fields to the CFT, the coefficient which enters in the twopoint function of the associated conserved vector currents governs the leading CFT contribution to the oneloop function of the corresponding gauge coupling. In particular, if this contribution is too large, it gives rise to unwanted subPlanckian Landau poles. Motivated by physics beyond the Standard Model, where , we have numerically studied the lower bounds on the coefficient using techniques developed in ref.[11]. Possible hierarchy problems are avoided by demanding that all scalar operators in the spectrum which are singlets under the global symmetry have dimensions . More specifically, we have considered lower bounds on extracted from crossing symmetry constraints on fourpoint functions of scalar operators in the fundamental representation of SO, or the bifundamental representation of SO