We present a phenomenological analysis of the recent Tevatron results on the forward-backward asymmetry and invariant-mass spectrum assuming a new contribution from an -channel gluon resonance with a mass in the range from to GeV. In contrast to most of the previous works, this analysis shows that for masses below TeV resonant New Physics could accommodate the experimental data. In general, we find that axial-like couplings are preferred for light and top quark couplings, and that only top quark couples strongly to New Physics. We find that composite model scenarios arise naturally from only phenomenological analyses of the experimental results. We show that our results are compatible with recent LHC limits in dijet and production, and find some tension for large resonance mass TeV. We indicate as best observables for discriminating a relatively light new gluon a better resolution in CDF forward-backward asymmetry, as well as the charge asymmetry and invariant-mass spectrum at the LHC.
Phenomenology of a light gluon resonance in top-physics at Tevatron and LHC
Ezequiel Álvarez111sequi@unsl.edu.ar, Leandro Da Rold222daroldl@cab.cnea.gov.ar, Juan Ignacio Sanchez Vietto333jisanchez@df.uba.ar, Alejandro Szynkman444szynkman@fisica.unlp.edu.ar
CONICET, INFAP and Departamento de Física, FCFMN, Universidad de San Luis
Av. Ejército de los Andes 950, 5700, San Luis, Argentina
[0.2cm] CONICET, Centro Atómico Bariloche and Instituto Balseiro
Av. Bustillo 9500, 8400, S. C. de Bariloche, Argentina
[0.2cm] CONICET, IFIBA and Departamento de Física, FCEyN, Universidad de Buenos Aires
Ciudad Universitaria, Pab.1, (1428) Buenos Aires, Argentina
[0.2cm] IFLP, CONICET - Dpto. de Física, Universidad Nacional de La Plata, C.C. 67, 1900 La Plata, Argentina
The Standard Model (SM) has accurately explained with great success the physics up to the GeV scale. However, there are several hints that indicate that New Physics (NP) could be expected at the TeV scale. Among the main motivations to look for NP at this scale, we can mention the quadratic divergences of the Higgs mass and the mechanism responsible for spontaneous symmetry breaking.
To explore the TeV-scale physics, physicists have built two colossal machines, Tevatron and LHC. Tevatron has not found signals of the Higgs boson, but it has studied as never before the top quark sector. The LHC has begun its search for the Higgs boson, which could take a few years. In the mean time, the LHC is becoming the largest top factory ever built.
Experiments in Tevatron and LHC have a special focus in top quark physics since, from a theoretical point of view, it is expected that the top quark could be a window to the NP . In fact, mainly because its extraordinary large mass and secondly because its poor exploration insofar, it is expected that the study of top quark production and the analysis of its properties, as couplings, mass and width, may throw some light on how is the NP. As a matter of fact, CDF and D0 have published several articles [2, 3, 4] in the last few years on pair-production, where in CDF results the forward-backward asymmetry () measurement did not agree with the SM NLO (next-to-leading order) expectation, whereas the cross section () and its invariant-mass () spectrum did agree with the SM NLO expectation. These results have raised the attention of the theoretical community which has essayed their hypothesis in several and varied publications . In a recent CDF article  on this issue, it was found that the total agrees with the SM expectation, but when studied as a function of the invariant mass of the top pair a large deviation () is found for large invariant mass. On the other hand, the cross section and invariant-mass spectrum still agree with the SM expectation. The results in that article, which have also motivated many theoretical publications , are the motivation for the present work.
CDF has reported in Ref.  that the forward-backward asymmetry measured in processes deviates from the SM prediction when the invariant mass of the pair is above GeV. This deviation, when analyzed from the not-unfolded results, consists in 4 bins in the range in which, although with little statistical significance, all the measurements are above the SM expected. When this data is unfolded and computed in only 2 bins for GeV, then a disagreement is found in the higher invariant mass bin.
Previous works  have proposed new particles in - and -channel, and non-resonant NP as effective operators to solve this anomaly. Moreover, in all the -channel proposals, the resonance mass is always pushed above TeV in order to avoid spoiling the invariant-mass spectrum111During the preparation of this article , a work analyzing the possibility of a light -channel resonance came out . In the region of the parameter space shared with the present article we reach similar conclusions.. In this work, motivated by a qualitative analysis of the not-unfolded data on the in Ref.  –which has much better resolution than the unfolded results–, we explore the possibility of having an -channel contribution of a new particle whose mass could be as light as GeV. In fact, the excess in the not-unfolded data in may be suggesting the well studied tail of a non-resonant contribution, as well as a light ( GeV) resonance peak. Henceforth, we have explored the phenomenology of an -channel colour-octet resonance with chiral couplings in a mass range from GeV to GeV. We have required for the light resonance to be weakly coupled to light quarks in order to avoid spoiling the invariant-mass spectrum and the colliders jet phenomenology. As it is shown in the article, these and other phenomenological requirements give rise naturally to composite model scenarios where only the top quark couples strongly to the NP. Moreover, we found that axial-like couplings for light and top quarks are preferred in all the studied mass range. Although we found that larger masses are preferred, these are more likely to be excluded by the dijet phenomenology.
It is worth stressing at this point that, although the motivation for this work came from a qualitative analysis of the above mentioned not-unfolded data, the quantitative analysis has been performed against the unfolded data.
This work is divided as follows. In the next section we present the effective model to be used in the analysis. In section 3 we perform an analytic study of the process , paying special attention to the resonant and interfering terms in order to determine which are the best regions of the parameter space to look for agreement with the CDF observables and -spectrum. In section 4 we present the results of Monte Carlo simulations, exploring the parameter space of the effective model which satisfies the CDF observables. Section 5 analyzes the favourable points in parameter space found in previous section within the existing constraints and perspectives for LHC. Section 6 contains an analysis and discussion of the favourable points within the framework of some popular models as warped extra-dimensions and composite models. Section 7 contains the conclusions and final ideas.
2 Effective model
We consider a model with a light gluon resonance transforming as an octet of SU(3) and mass TeV. Such a scenario is naturally obtained within the framework of composite models and extra dimensions, since they typically contain this class of resonances . The interactions with the SM-quarks are given by:
where are the chiral couplings between the gluon resonance and the SM quarks, in units of the strong coupling . For simplicity we will consider that the light quarks have the same chiral couplings:
and we will allow different couplings for the third generation:
In warped extra-dimensional and composite Higgs models is usually correlated with the Yukawa couplings, leading to . However the quark content of the proton is mostly and quarks, thus at leading order in our model and only enter in the width of the resonance. As we discuss next, also the recent experimental results from Tevatron prefer , leading to a resonance width mostly determined by the couplings . In this scenario top-pair creation in proton collisions are not sensitive to small differences between the couplings of the first and second generations, as long as they remain small compared with the couplings of the third generation, justifying our simplified assignment. The anomaly in measured at LEP and SLC can be solved if is partially composite, leading to larger couplings for than for the light generations. For this reason we consider as an independent coupling and allow it to be larger than Ã¢ [12, 13, 15, 16]. For , the resonance width is comparable to the resonance mass: , therefore we will consider an energy dependent width for the resonance.
3 Analytic study of the process
Within the framework of the effective model considered in the previous section, we perform an analytic parton-level analysis of both the forward-backward asymmetry and the cross section of pair production at Tevatron. The NP contribution to production in Tevatron comes from a gluon resonance exchange in the -channel of the process. Henceforth, in order to obtain a qualitative understanding of the asymmetry and the cross section at Tevatron, we study the -dependent tree-level parton formulae in the center of mass frame of the system. In the following section, we refine the analysis turning to a Monte-Carlo simulation of the process and include the SM-NLO result that gives a non-vanishing contribution to the SM expected value of .
When considering as the final state, we can write the differential angular cross section in the center of mass reference frame of the system as
where INT and NPS stand for SM-NP interference and new physics squared, respectively. At tree-level, including only LO QCD and the gluon resonance contributions, the three different terms read
The SM contribution corresponds to the -channel gluon exchange diagram, the NPS term is obtained by replacing the gluon by in that diagram, and INT stands for the gluon- interference. The angle is defined by the directions of motion of the top quark and the incoming light quark (up quark, for instance) in the center of mass frame of the top pair system. The coefficient denotes the vector(axial) coupling between a quark of the light generations and in units of , whereas is an analogous notation for the top couplings 222The vector and axial couplings are given in terms of the chiral couplings introduced in Section 2 as
and .. Finally, represents the mass of , is the energy in the center of mass frame of the system and is the velocity of the top in that frame. Since the narrow-width approximation is no longer valid for large couplings, the use of a modified Breit-Wigner function is needed for the propagator in Eqs. (6,7), where is an energy dependent width.
The forward and backward cross sections are defined as
and is given by
Note that just the constant and terms contribute to , whereas only the term gives a contribution to the numerator of .
If the NP effects in the cross section are small compared to the leading SM contribution, (where contains both INT and NPS contributions), we can approximate by
for each bin in the invariant-mass distribution. The numerator of includes INT+NPS terms given at tree-level by Eqs. (6) and (7), whereas the denominator is dominated by the SM tree-level contribution given by Eq. (5). We use the NLO result for in order to take into account the SM contribution to the asymmetry. We have checked that the condition is satisfied in our simulations.
The main point of this section is to gain a preliminar insight on which parameter space regions are favored by the experimental results. According to this, we analyze Eqs. (5-7) and look for general conditions to be fulfilled by the couplings. In particular, motivated by the not-unfolded data on (see Fig. 10 in Ref. ), we pursue to reproduce the peak around 650 GeV as a resonant effect arising from a with GeV and, at the same time, to produce small corrections to the invariant mass distribution in the cross section. As a first observation, we notice that the NPS contributions to dominate over the INT terms in the resonant region. Since a positive asymmetry is measured in that region, for typical values GeV, a first constraint is imposed on the couplings: . Subsequently, it is also required in order to assure a positive pre-peak since INT contributions may become competitive in the region where the pre-peak takes place. In addition, these two conditions imply that , giving positive and negative contributions to from INT terms in the regions with and , respectively. In order to avoid large distortions in the invariant mass profile and significant deviations from the dijet final state invariant mass distributions, we demand small axial and vector couplings for light quarks –we satisfy this condition by assuming small chiral couplings–. On the other side, we set large top chiral couplings for two reasons: (i) to make the effects detectable and (ii) to explain the shape by means of the presence of a relatively broad resonance –we will see below that a narrow resonance does not lead to an agreement with the experimental results. Besides, NPS contributions to are always positive for any choice of couplings as we can read from Eq. (7). Therefore, we expect the cross section to increase in the resonant region. Thus, we see that a delicate interplay between the size of width and its couplings to quarks for fixed values of dictates the experimentally allowed regions in the parameter space. It is worth to notice that a clear signal of the presence of a light resonance would be to detect an excess of events in the invariant mass distribution of and for values below the resonant peak and a defect above.
The previous discussion is quite general for a gluon resonance with TeV and 0.2-0.3. Now, with the aim to illustrate the way we use that argument as a guide to generate the simulations presented in the next section, we take a gluon resonance with GeV as a typical example and analyze the role played by the couplings in the search of suitable and distributions. We also set, in this example, the quarks Left couplings to zero to simplify the analysis. In Fig. 1, we show the analytical predictions for and at the parton-level (these results have to be convoluted with the PDF’s to obtain the physical values). For zero Right couplings, we obtain the SM cross section at tree-level (blue line) and a null contribution to the asymmetry. Now, to achieve the condition , we need and to have opposite signs –the other constraint, , is then directly fulfilled–. Besides, as a relative large asymmetry is required, we set first both couplings to have an absolute value equal to 1. This selection results in too large deviations from measurements and it is forbidden (green lines). We reduce then the value of the light quark right coupling, , to produce a lower impact on the observables, however, this choice leads to a rather narrow resonance that still gives a large deviation in the resonant region (brown lines). Finally, in red lines, we find a possible solution. We see that the inadequate effect produced by a narrow-width may be changed by making the resonance broader. We keep a small value for but increase up to 5, giving a width GeV, which is of the same order as the typical values for allowed parameter space points found in the simulations. Therefore, this analysis shows us the main ingredients to be taken into account and anticipates light-quark (top) chiral couplings of order . Besides, there could be new channels open, for example if new coloured fermions were present, the processes and could be open whenever the available phase space allows them . An alternative mechanism to explain the experimental data may occur in that case because would increase without the need of enlarging the couplings. In this work we assume that these channels are not open.
For large masses ( TeV), the INT term dominates in the whole range, always giving positive contributions to and . In order to have an agreement with , the couplings of the light quarks have to be 1 because of the lack of a significant NPS contribution in the resonant region. Finally, the dependence on is smoother for heavy gluons than in the case of light gluons.
We see from Eqs. (6) and (7) that the strength of the INT and NPS terms are controlled by the product of the light and top couplings. On the other hand, the width is given by a sum of light and top couplings squared (and a secondary dependence on the mass of the quarks). Therefore, when exploring the parameter space we have that if –as it is going to be the case– the light couplings are small and the top couplings large, then modifications to the width arise mainly from changes in the top couplings. On the other hand, variations in the strength of INT and NPS terms are equally sensitive to changes in both the light and top couplings. Therefore, in a first approximation, if we need to modify the strength of the INT and NPS contributions without producing large variations in the width, we should only change the light-quark couplings; if we need to modify the width, we should change the top couplings. This discussion will be helpful for the analysis performed in the next section.
4 Numerical scan and results
In this section we perform a numerical scan of the regions of parameter space pointed out by the previous paragraphs analytic study of the forward-backward asymmetry () and invariant-mass spectrum in the process. In this numerical scan we use Madgraph/Madevent  (MGME) 333We have slightly modified MGME to include a running width for the resonance. with CTEQ6L parton distribution functions  to simulate Tevatron process and compute and -spectrum at parton-level.
In order to test the different points in the parameter space, we compare our simulated results to the CDF two-bins unfolded result on this observable . On the other hand, given the lack of unfolded results in the recent experimental data on the -spectrum , we interpret the experimental data as that there is agreement with the SM expectation. Therefore, we simulate the -spectrum using our model and using the SM and compare both results. In this last comparison we use in both simulations the same values for the top mass, and the factorization and renormalization scale. (In any case, we have checked that changes in these values within their allowed limits do not modify our results within the statistical error.)
In order to compare the forward-backward asymmetry results, we perform a test on the two bins of the unfolded data for . To compare the invariant-mass spectrum we also perform a test, but as follows. We divide the domain in GeV bins up to GeV and from there up we divide in three bins in , GeV and greater than GeV, in order to assure at least 5 events per bin in the SM simulated -spectrum with (the collected luminosity in the CDF results). To avoid large fluctuations we always simulate events, and for the -spectrum we scale the results to the measured luminosity. In each test we require to state agreement of the NP model on the corresponding observable.
We perform a scanning of two different regions in parameter space. In the first case, in order to obtain clean results, we set all the Left couplings of the quarks to the new gluon to zero, and plot our results solely as a function of the light and top quarks Right couplings. The motivation for this choice in the top quarks is simplicity and the existence of precision measurements in -physics  which constrain the couplings of the doublet to new gluons. In this scenario we have that axial and vector couplings are equal; therefore the asymmetry at the resonance is positive (see Eq. 7). By setting opposite signs to the light and top Right couplings we get a positive pre-peak in the asymmetry, as sought. In the second case we perform a random scan of all couplings in order to see the regions which have better agreement with and the -spectrum .
In Figure 2 we plot the results for the scan where all the Left couplings are set to zero. We have performed a uniform density random scan in this region and use the black crosses to state no agreement neither in nor in -spectrum , the magenta squares to state only agreement, the blue triangles to state only -spectrum agreement and the red points to state agreement in both observables. After the first scan we have doubled the random density close to the red points region in order to increase the contrast therein.
These plots should be read as follows. First, observe that the SM, which is the point at the origin, does not agree with the forward-backward asymmetry. Then notice that as couplings increase (in absolute value) from zero, the only agreement is, as expected, with -spectrum (blue triangles). For larger couplings this agreement has to break down. On the other hand, from the top and left region of the plots we have –as expected from previous section– agreement with the forward-backward asymmetry (magenta squares). (In Figure 2a the top left corner begins with no agreement, since such a light resonance gives too large asymmetry for those couplings.) Again, this agreement has to break down as couplings decrease in absolute value. Henceforth, these two regions (blue triangles and magenta squares) could be joined mainly in two ways. Either through a region where both agreements ( and -spectrum) fail (black crosses), or either they overlap and there is an agreement in both observables (red points).
We see from the plots in Fig. 2 that, as it was the first motivation of the work, there is a strip of favourable (red) points for a light resonance beginning at GeV. Although in all cases the agreement region is a relatively thin strip, this agreement is slightly more convincing for larger . It is interesting to observe at this point that the favourable region that has raised from the phenomenological analysis is composite models friendly, with only the Right chirality of the top quark being composite in this case, since it has large couplings to New Physics. A possible way to broaden the red strip is to consider the presence of extra fermions that open new decay channels modifying the resonance width . However, for the case of light fermions, a deeper examination of the existing collider phenomenology should be carried on.
At this point, it is also interesting to understand the red points patterns in Fig. 2 from the previous section parton-level equations (6) and (7). The off-resonant case (large ) is easy to understand, since in this case the dominant contribution to -production comes from the INT term, Eq. (6). As previously discussed, the effect of the width in this term may be neglected in a first approximation and, therefore, the dependence on the couplings is just through the product. Therefore, regions where this product remains constant will have approximately the same outcome for and -spectrum . In fact, this is easily verified for GeV in Fig. 2d, where the red points strip is along a constant line.
On the other hand, in the resonant case (low ) it is the NPS term, Eq. 7, the one that rules, at least close to the resonance. In this case the width plays a major role and it has to be taken into account. We observe from Eq. 7 that, for the case of large and little , the NPS term close to the resonance goes essentially as . Therefore, an increase in will have the effect of diminish and broaden the peak in the -spectrum and . Since the resolution in is only in two large bins, then the change will be minor here. Henceforth, we may expect that if only is increased from a favourable (red) point in parameter space, then the agreement with -spectrum and may hold. On the other hand, if we only move we will have that the width remains unchanged and the NPS changes to disfavourable regions. In fact, if is decreased, then the peak in decreases but is not broaden, therefore spoiling the asymmetry. If is increased then the peak in -spectrum is also increased and spoils this observable. These observations are best seen in Figs. 2a and 2b for and GeV, where the red points lie along a vertical strip.
Finally, we have also performed a total random scan of the couplings. We present the results of this scan for masses of the resonance which could be resonant at Tevatron ( and GeV) in Fig. 3, and for larger masses ( and GeV) in Fig. 4. In all cases we have varied , , , and . The variation of , which is between and (in units of ) and is not plotted, modifies slightly the width of and henceforth the outcome of the simulation. At difference with the previous case, in each of these plots there are three couplings which change randomly and are not being plotted. Therefore, there is not a clear border between the different colour regions in this case. A general feature of all plots is that the favourable points show an axial-like preference for the couplings of light and top quarks. Moreover, we found that, as anticipated in the previous section, these axial couplings have in general opposite signs (this result cannot be seen from the plots). Also notice that, as expected, as increases also the couplings required to satisfy the CDF observables increase (in absolute value). And last, observe that the density of red points increases considerably as increases, which would mean that a large solution is favoured.
As predicted in the previous section, as the mass of increases, the NPS term becomes less important and therefore the large-width requirement to preserve the -spectrum agreement may be relaxed. This is translated in Figs. 3 and 4 in that magenta squares density decreases (where agrees, but -spectrum fails) and red points density increases ( and -spectrum agree), as increases from GeV to GeV.
5 LHC phenomenology
Our analysis shows that a light gluon resonance with GeV can accommodate the CDF data on production within . However, there are new data from LHC that show agreement with the SM predictions. The most important constraints for this type of gluon resonance searches arise from the production [21, 22] and dijet final states [23, 24] in LHC at 7 TeV. In this section we analyze those constraints for all the points of the parameter space that reproduce the Tevatron measurements, using the random scan of section 4. We also show the main predictions of our model for LHC.
Concerning top-pair production, Ref.  shows agreement of the cross section with the SM for an accumulated luminosity of 35 pb, whereas the recent analysis of  sets limits on the production of resonances with an accumulated luminosity of 200 pb. Ref.  excludes a gluon resonance with GeV for a specific choice of couplings: , and . Note that, as shown in the previous sections, the Tevatron data disfavours vector-like couplings as the ones chosen for the light quarks in , since they lead to too low . Since the Tevatron requires large couplings between the top-quark and the resonance, , the top-pair pair production is one of the better tests for a gluon resonance of this kind at LHC. We have simulated the production cross section of a gluon resonance decaying to at LO at parton-level, finding pb in all the cases. Since the analysis of  gives upper bounds for TeV only, we can not precisely check limits in the case of TeV. For a light resonance mass, TeV, the limits set by  are well beyond the results of our model, mainly because the Tevatron results can be obtained with small couplings for the light quarks, suppressing the production through at LHC. For TeV, the ATLAS sensitivity reported in  is near the region of the parameter space that can explain . We expect that with an accumulated luminosity greater than 1 fb, either an excess in can be measured or a large region of the parameter space with TeV can be excluded, mostly the region with larger top-couplings. Also a better resolution of the invariant mass distribution could be useful to test our predictions for .
In Ref.  the ATLAS collaboration has provided experimental results on the search for new resonances in dijet final states with 36 pb. A recent update with 163 pb provides stronger constraints on the search for new resonances and extends the analysis to larger masses . A light gluon resonance can give large contributions to jet pair creation, as shown in Fig. 2 of Ref.  for axial couplings of order 1. Since the Tevatron data selects the couplings of the light quarks, we analyze the predictions for dijet production at LHC, that have a strong dependence with the resonance mass. First notice that, since the events analyzed by ATLAS include -jets, the dijet production through is controlled by two terms, one including only light-jets proportional to and another including -jets proportional to . For GeV, we find agreement with the Tevatron data for regions of the parameter space with , whereas . These couplings lead to a suppression factor for light-jets, and a suppression factor for -jets, with a dispersion that depends on the precise value of . On the other hand, for TeV, larger are required to reproduce the Tevatron results on , , that is similar to the range for , but smaller than the range for . Therefore in this case, although there is a large suppression from , the large couplings induce a sizable contribution that can be tested with the present data. We have simulated at LO the dijet production in our model with MGME and Pythia  for hadronization and showering simulations. We have considered the following set of kinematic cuts, implemented in : , , GeV, GeV and GeV. Although we have not made a detailed simulation of the detectors, neither computed NLO corrections, we expect that our LO predictions can give the correct order of magnitude, mainly on the upper bound limits. Our results show that, for TeV, the dijet limits on the cross section are all beyond the present sensitivity. For TeV the ATLAS analysis put severe constraints on the parameter space. Taking the confidence level upper limit of , we obtain that a large region of the parameter space is ruled out in this case.
Let us comment also on the proportion of events with light- and -jets. For and 850 GeV we find that the ratio of events with -jets is between for the different points of the parameter space, whereas for GeV the ratio of -jets is and for GeV the ratio is . The dispersion in the frequency of -jets is due to the different values of and for each point of the parameter space with fixed mass. On the other hand, for TeV, we find points of the parameter space where the ratio of -jets is only , and we also find points where this ratio raises to , depending on the region of the parameter space. As expected, we also find that, for all the resonance masses, the points of the parameter space that have a lower ratio of events with -jets correspond to lower values of . For these reasons it will be very interesting to have experimental data on resonance searches with dijets and -tagging, that could help to better select the favourable region of the parameter space.
We summarize the predictions and constraints from resonances searches in and creation at LHC in Fig. 5a. The different colours codify the different values of the resonance mass: red for GeV, orange for GeV, brown for TeV, green for TeV, and blue for TeV. We also show with coloured lines the upper bounds for each mass in the contribution to dijets and top-pair production at confidence-level, taken from Fig. 2 of  and Fig. 6 of . There is no upper bound line on production for TeV because Ref.  only provides information up to TeV. Only the points on the lower-left corner of these lines are not ruled out. Note that although the limits from do not exclude any point, a mild increase of the accumulated luminosity may be enough to test a resonance with TeV. For TeV the dijet constraints are much above from the predictions of the model, whereas for TeV approximately of the points are excluded by dijet constraints. In this case only points with are allowed, therefore the non-observation of any resonance in dijet final states selects the smaller couplings.
Recently, the CMS Collaboration has reported the first measurement of the charge asymmetry in top quark pair production at LHC, with an integrated luminosity of 36 pb at 7 TeV. The measured asymmetry , being still dominated by the statistical uncertainties , is consistent with the SM prediction . Increasing the integrated luminosity to fb, CMS can reach the same sensitivity as Tevatron results, providing an independent measurement of top asymmetries that can complement the Tevatron results on . Moreover, as shown in Ref. , there is a correlation between both asymmetries, that depends on the details of NP and can thus discriminate between different models. The charge asymmetry measured by CMS is defined by:
with the number of events with positive or negative values of , being the pseudo-rapidity in the laboratory frame, .
We have computed the charge asymmetry for top-pair production at LHC at 7 TeV. The NP contribution to the total charge asymmetry in production, , is predicted to be positive and of order in our model, below the sensitivity of the present measurements. Although there is a correlation with , it is different from the one reported for the axial model in , essentially because the couplings are neither axial, nor universal, resulting in a more complicated pattern that is difficult to disentangle from the signals of other models by just analyzing these observables: and . In Fig. 5b we show our results for the contributions to and . The different colours codify the resonance mass as explained before.
In  the CMS Collaboration has also provided the unfolded results on the spectrum of , considering 6 bins, defined by the following limits for : . Thus we have computed the spectrum of , obtaining deviations of order compared with the SM prediction for these 6 bins. We obtain the largest deviations in the bin , with positive NP contributions of order . However, since the uncertainties are still large, all the points we have simulated are in agreement with the data. Either larger luminosity or a more refined resolution is needed to test our model with the spectrum.
There is also an interesting prediction for LHC concerning the invariant mass spectrum. For all the points that agree with the Tevatron data, our model results in an excess of events for and a lack of events for , compared with the SM. This can be seen from the parton-level analysis of section 3. Therefore, detailed data on the -spectrum could provide a strong test for this kind of gluon resonance.
6 Discussion and UV completions
We have considered an effective approach with just one gluon-resonance. There are many UV completions of the SM that contain a massive colour-octet vector as , as for example composite Higgs models and extra-dimensions. In these theories there are in general resonances associated to the other SM fields also, and even infinite towers of massive resonances. We have considered just the effects of the lightest resonance associated to the gluon for simplicity, and because in many cases it gives the leading contributions to the physical processes that we have studied in this work. Examples of UV completions of our model have been presented recently in Refs.  and , corresponding to two different models with a warped extra dimension, the main difference between them being the resonance scale, TeV and TeV, respectively. Our results are compatible with the analysis of those papers, but the region of the parameter space that we explore is much larger than the one allowed in those frameworks. Usually, in warped-extra dimensional models, the hierarchy in the fermionic spectrum is explained by the localization of the 0-mode Kaluza-Klein (KK) wave functions in the extra dimension. These wave functions also determine the strength of the couplings with the massive KK modes, as . In our model we have considered random couplings, without that kind of constraints. As explained in , in a large family of composite Higgs models and in simple extra dimensional models it is not possible to obtain large axial couplings, because is either positive or small and negative. For this reason our results go beyond this kind of models. Interestingly, the results on choose a region with large top-couplings, characteristic of models with a composite top.
We have shown that, to reproduce , the couplings of the light quarks can not be too small for TeV. For warped extra-dimensions, this implies that at least one of the chiralities of the light quarks has to be somewhat delocalized from the UV boundary, leading to sizable couplings not only with , but also with the whole set of resonances associated to the EW sector. The precision EW measurements put strong constraints on these couplings, that usually can only be satisfied by introducing new symmetries [28, 29] or more structure [30, 31]. In Refs. [12, 13] we have shown how, by properly choosing the quantum numbers of the fermions under the custodial symmetry and an extra discrete symmetry (called and ), one can protect the -couplings of the light fermions and reproduce the measured .444See also Refs. [15, 32, 16] for other solutions. The same mechanism can be implemented in the present set-up for TeV.
In this work we have studied the phenomenology at Tevatron and LHC of a TeV colour-octet vector resonance. We have focused our attention on the top-pair creation at Tevatron, making an analytical study of the new physics contribution to the cross section and forward-backward asymmetry as functions of for different configurations. This analysis shows the interplay between the interference and the new physics square terms depending on the size of the couplings and the resonance mass. It allows us to select the regions of the parameter space that can increase the asymmetry in the high invariant mass region without spoiling the cross section. We have also simulated the top-pair creation at Tevatron at LO and at partonic level with Madgraph/Madevent. The previous analysis also allows us to understand the results of our simulations. Finally, we have computed the contributions of our model to dijet final states and top-pair creation at LHC, contrasting our results with the limits on the production of resonances. We have also shown the LHC predictions of our model for the charge asymmetry in creation.
As a summary of our results for Tevatron, we have shown that an -channel gluon resonance with a mass TeV and chiral couplings with the SM quarks can reproduce the data on top-pair production with , alleviating the tension present in the high mass regime of , that exceeds the SM prediction by more than three standard deviations. Almost all the previous analysis have focused either on effective operators or heavy resonances to keep the agreement with the SM cross section: TeV (see footnote 1), whereas we have shown that resonant new physics can also do the job. We have determined the couplings that reproduce the data for different values of the resonance mass: and 2500 GeV, demanding small couplings for the light quarks to avoid spoiling the LHC jet phenomenology. We have found that in general axial-like couplings are preferred and that, since we have chosen small light-quark couplings, the top couplings are large , pointing towards top compositeness. The size of the couplings of the light-quarks increases as the resonance becomes heavier to keep a sizable contribution to the asymmetry, with for GeV and for TeV. We have also found that the density of points in the parameter space that can accommodate the data is much smaller for GeV than for TeV.
Concerning the LHC phenomenology, we have found that for and 1500 GeV, the limits on the search for new resonances are above the predictions of our model. However, for TeV, of the parameter space is ruled out by the limits on new resonances set by dijets, remaining only the points with light-quark couplings not larger than . Our results also show that a large region of the parameter space that lies below the present limits could be tested with a sensible rise of the accumulated luminosity and the center of mass energy. It will be interesting to have experimental data on dijet resonance searches with -tagging in order to discriminate regions in the parameter space. We have also computed the top charge asymmetry, obtaining new physics contributions of order , and corrections to the spectrum.
Note that, if the Tevatron data is explained by a gluon resonance, we have to demand GeV, with light-quark couplings increasing with the resonance mass. On the other hand, the limits from LHC resonance searches seem to avoid TeV if the couplings are large. Therefore, our analysis suggests that there is a preferred window for the resonance mass: TeV. Interestingly, a resonance like this can be tested at LHC in the next years.
Ref.  has provided a distribution of the events and of the asymmetry as functions of the invariant mass, with bins of 15 GeV for the cross section and 50-100 GeV for the asymmetry. However these distributions are not unfolded, making it difficult to test the models beyond the SM with this level of accuracy. Notice that the parton-level unfolded data presented for the asymmetry is divided in two bins only . A more refined distribution at the parton-level would be very useful to better test the new physics.
As a final consideration, note that for certain regions of the parameter space, our model can be embedded in a more complete theory solving the hierarchy problem, as composite Higgs models or extra dimensions. However there are some combinations of couplings that either can not be obtained in these models, or fail to reproduce the EW precision data when the full EW sector is considered. It would be very interesting to find a model able to cover those regions of the parameter space.
We thank José Santiago for a careful reading of the manuscript and suggestions. We also thank Carlos Wagner, Ricardo Piegaia, María Teresa Dova, Xavier Anduaga and Cedric Delaunay for useful conversations. This work has been partially supported by ANPCyT (Argentina) under grant # PICT-PRH 2009-0054 and by CONICET (Argentina) PIP-2011.
-  C. T. Hill, S. J. Parke, Phys. Rev. D49 (1994) 4454-4462. [arXiv:hep-ph/9312324 [hep-ph]]; C. T. Hill, E. H. Simmons, Phys. Rept. 381 (2003) 235-402. [hep-ph/0203079]; K. Agashe, R. Contino, R. Sundrum, Phys. Rev. Lett. 95 (2005) 171804. [hep-ph/0502222].
CDF Public Note: http://www-cdf.fnal.gov/physics/new/top/2010/tprop
-  T. Aaltonen et al. [ CDF Collaboration ], Phys. Rev. Lett. 101, 202001 (2008) [arXiv:0806.2472 [hep-ex]].
-  V. M. Abazov et al. [ D0 Collaboration ], Phys. Rev. Lett. 100, 142002 (2008) [arXiv:0712.0851 [hep-ex]].
-  G. Burdman, L. de Lima, R. D. Matheus, Phys. Rev. D83 (2011) 035012. [arXiv:1011.6380 [hep-ph]]; B. Xiao, Y. -k. Wang, S. -h. Zhu, [arXiv:1011.0152 [hep-ph]]; C. Degrande, J. -M. Gerard, C. Grojean, F. Maltoni, G. Servant, JHEP 1103 (2011) 125. [arXiv:1010.6304 [hep-ph]]; C. -H. Chen, G. Cvetic, C. S. Kim, PHLTA,B694,393-397. 2011 B694 (2011) 393-397. [arXiv:1009.4165 [hep-ph]]; J. A. Aguilar-Saavedra, Nucl. Phys. B843 (2011) 638-672. [arXiv:1008.3562 [hep-ph]]; M. Bauer, F. Goertz, U. Haisch, T. Pfoh, S. Westhoff, JHEP 1011 (2010) 039. [arXiv:1008.0742 [hep-ph]]; R. S. Chivukula, E. H. Simmons, C. -P. Yuan, Phys. Rev. D82 (2010) 094009. [arXiv:1007.0260 [hep-ph]]; B. Xiao, Y. -k. Wang, S. -h. Zhu, Phys. Rev. D82 (2010) 034026. [arXiv:1006.2510 [hep-ph]]; Q. -H. Cao, D. McKeen, J. L. Rosner, G. Shaughnessy, C. E. M. Wagner, Phys. Rev. D81 (2010) 114004; [arXiv:1003.3461 [hep-ph]]; J. Cao, Z. Heng, L. Wu, J. M. Yang, Phys. Rev. D81 (2010) 014016. [arXiv:0912.1447 [hep-ph]]; I. Dorsner, S. Fajfer, J. F. Kamenik, N. Kosnik, Phys. Rev. D81 (2010) 055009. [arXiv:0912.0972 [hep-ph]]; A. Arhrib, R. Benbrik, C. -H. Chen, Phys. Rev. D82 (2010) 034034. [arXiv:0911.4875 [hep-ph]]; J. Shu, T. M. P. Tait, K. Wang, Phys. Rev. D81 (2010) 034012. [arXiv:0911.3237 [hep-ph]]; P. H. Frampton, J. Shu, K. Wang, Phys. Lett. B683 (2010) 294-297. [arXiv:0911.2955 [hep-ph]]; K. Cheung, W. -Y. Keung, T. -C. Yuan, Phys. Lett. B682 (2009) 287-290. [arXiv:0908.2589 [hep-ph]]; S. Jung, H. Murayama, A. Pierce, J. D. Wells, Phys. Rev. D81 (2010) 015004. [arXiv:0907.4112 [hep-ph]]; P. Ferrario, G. Rodrigo, Phys. Rev. D80 (2009) 051701. [arXiv:0906.5541 [hep-ph]]; M. V. Martynov, A. D. Smirnov, Mod. Phys. Lett. A24 (2009) 1897-1905. [arXiv:0906.4525 [hep-ph]]; K. Kumar, T. M. P. Tait, R. Vega-Morales, JHEP 0905 (2009) 022. [arXiv:0901.3808 [hep-ph]]; and Ref. [13, 14].
-  T. Aaltonen et al. [ CDF Collaboration ], [arXiv:1101.0034 [hep-ex]].
-  E. Gabrielli, M. Raidal, [arXiv:1106.4553 [hep-ph]]; Y. Cui, Z. Han, M. D. Schwartz, [arXiv:1106.3086 [hep-ph]]; J. A. Aguilar-Saavedra, M. Perez-Victoria, [arXiv:1105.4606 [hep-ph]]; M. Duraisamy, A. Rashed and A. Datta, [arXiv:1106.5982 [hep-ph]]; D. Krohn, T. Liu, J. Shelton, L. -T. Wang, [arXiv:1105.3743 [hep-ph]]; K. S. Babu, M. Frank, S. K. Rai, [arXiv:1104.4782 [hep-ph]]; D. -W. Jung, P. Ko, J. S. Lee, [arXiv:1104.4443 [hep-ph]]; S. Jung, A. Pierce, J. D. Wells, [arXiv:1104.3139 [hep-ph]]; C. Degrande, J. -M. Gerard, C. Grojean, F. Maltoni, G. Servant, [arXiv:1104.1798 [hep-ph]]; X. -P. Wang, Y. -K. Wang, B. Xiao, J. Xu, S. -h. Zhu, PHRVA,D83,115010. 2011 D83 (2011) 115010. [arXiv:1104.1917 [hep-ph]]; A. E. Nelson, T. Okui, T. S. Roy, [arXiv:1104.2030 [hep-ph]]; J. A. Aguilar-Saavedra, M. Perez-Victoria, PHLTA,B701,93-100. 2011 B701 (2011) 93-100. [arXiv:1104.1385 [hep-ph]]; A. Rajaraman, Z. ’e. Surujon, T. M. P. Tait, [arXiv:1104.0947 [hep-ph]]; J. Shu, K. Wang, G. Zhu, [arXiv:1104.0083 [hep-ph]]; S. Jung, A. Pierce, J. D. Wells, [arXiv:1103.4835 [hep-ph]]; M. I. Gresham, I. -W. Kim, K. M. Zurek, [arXiv:1103.3501 [hep-ph]]; Z. Ligeti, M. Schmaltz, G. M. Tavares, JHEPA,1106,109. 2011 1106 (2011) 109. [arXiv:1103.2757 [hep-ph]]; C. Grojean, E. Salvioni, R. Torre, [arXiv:1103.2761 [hep-ph]]; J. A. Aguilar-Saavedra, M. Perez-Victoria, JHEP 1105 (2011) 034. [arXiv:1103.2765 [hep-ph]]; C. Delaunay, O. Gedalia, Y. Hochberg, G. Perez, Y. Soreq, [arXiv:1103.2297 [hep-ph]]; R. Foot, Phys. Rev. D83 (2011) 114013. [arXiv:1103.1940 [hep-ph]]; N. Craig, C. Kilic, M. J. Strassler, [arXiv:1103.2127 [hep-ph]]; E. R. Barreto, Y. A. Coutinho, J. Sa Borges, Phys. Rev. D83 (2011) 054006. [arXiv:1103.1266 [hep-ph]]; A. R. Zerwekh, [arXiv:1103.0956 [hep-ph]]; G. Isidori, J. F. Kamenik, Phys. Lett. B700 (2011) 145-149. [arXiv:1103.0016 [hep-ph]]; K. M. Patel, P. Sharma, JHEP 1104 (2011) 085. [arXiv:1102.4736 [hep-ph]]; B. Grinstein, A. L. Kagan, M. Trott, J. Zupan, PRLTA,107,012002. 2011 107 (2011) 012002. [arXiv:1102.3374 [hep-ph]]; V. Barger, W. -Y. Keung, C. -T. Yu, Phys. Lett. B698 (2011) 243-250. [arXiv:1102.0279 [hep-ph]]; M. I. Gresham, I. -W. Kim, K. M. Zurek, [arXiv:1102.0018 [hep-ph]]; Y. Bai, J. L. Hewett, J. Kaplan, T. G. Rizzo, JHEP 1103 (2011) 003. [arXiv:1101.5203 [hep-ph]]; K. Cheung, T. -C. Yuan, Phys. Rev. D83 (2011) 074006. [arXiv:1101.1445 [hep-ph]]; and Refs. [9, 16, 17].
-  E. Alvarez, L. Da Rold, J.I. Sanchez Vietto and A. Szynkman, “Weakly coupled resonant NP at the Tevatron t-tbar forward-backward asymmetry,” presented at PLANCK 2011 - From the Planck Scale to the ElectroWeak Scale, Lisbon, May 30th to June 3rd (2011).
-  R. Barcelo, A. Carmona, M. Masip, J. Santiago, [arXiv:1105.3333 [hep-ph]].
-  H. Davoudiasl, J. L. Hewett and T. G. Rizzo, Phys. Lett. B 473 (2000) 43 [arXiv:hep-ph/9911262]; A. Pomarol, Phys. Lett. B 486 (2000) 153 [arXiv:hep-ph/9911294].
-  R. Contino, T. Kramer, M. Son, R. Sundrum, JHEP 0705 (2007) 074. [hep-ph/0612180].
-  L. Da Rold, JHEP 1102 (2011) 034 [arXiv:1009.2392 [hep-ph]].
-  E. Alvarez, L. Da Rold and A. Szynkman, JHEP 1105 (2011) 070 [arXiv:1011.6557 [hep-ph]].
-  A. Djouadi, G. Moreau, F. Richard and R. K. Singh, Phys. Rev. D 82 (2010) 071702 [arXiv:0906.0604 [hep-ph]].
-  A. Djouadi, G. Moreau and F. Richard, Nucl. Phys. B 773 (2007) 43 [arXiv:hep-ph/0610173].
-  A. Djouadi, G. Moreau and F. Richard, arXiv:1105.3158 [hep-ph].
-  R. Barcelo, A. Carmona, M. Masip, J. Santiago, [arXiv:1106.4054 [hep-ph]].
-  J. Alwall, P. Demin, S. de Visscher, R. Frederix, M. Herquet, F. Maltoni, T. Plehn, D. L. Rainwater et al., JHEP 0709 (2007) 028, [arXiv:0706.2334 [hep-ph]].
-  J. Pumplin, D. R. Stump, J. Huston, H. L. Lai, P. M. Nadolsky and W. K. Tung, JHEP 0207 (2002) 012 [arXiv:hep-ph/0201195].
-  S. Chang, C. S. Kim, J. Song, JHEP 0702 (2007) 087. [hep-ph/0607313]; C. Csaki, D. Curtin, Phys. Rev. D80 (2009) 015027. [arXiv:0904.2137 [hep-ph]].
-  ATLAS-CONF-2011-023, the Single Lepton+Jets Channel without b-tagging” [ATLAS collaboration]; ATLAS-CONF-2011-034, cross section with ATLAS in pp collisions at sqrt(s) = 7 TeV in dilepton final states” [ATLAS collaboration]; ATLAS-CONF-2011-035, cross section with ATLAS in pp collisions at sqrt(s) = 7 TeV in the single-lepton channel using b-tagging” [ATLAS collaboration]; ATLAS-CONF-2011-040, production cross section using dilepton and single-lepton final states” [ATLAS collaboration]; ATLAS-CONF-2011-066, all-hadronic channel in ATLAS with = 7 TeV data” [ATLAS collaboration].
-  ATLAS-CONF-2011-087, “A Search for Resonances in the Lepton Plus Jets Channel using 200 pb of Collisions at =7 TeV”, [ATLAS collaboration].
-  G. Aad et al. [ ATLAS Collaboration ], New J. Phys. 13 (2011) 053044. [arXiv:1103.3864 [hep-ex]].
-  ATLAS-CONF-2011-081, “Update of the Search for New Physics in the Dijet Mass Distribution in 163 pb of pp Collisions at =7 TeV Measured with the ATLAS Detector,” [ATLAS collaboration].
-  T. Sjostrand, S. Mrenna and P. Z. Skands, JHEP 0605 (2006) 026 [arXiv:hep-ph/0603175].
-  CMS-PAS-TOP-10-010, “Measurement of the charge asymmetry in top quark pair production, with the CMS experiment”.
-  J. A. Aguilar-Saavedra and M. Perez-Victoria, arXiv:1105.4606 [hep-ph].
-  K. Agashe, A. Delgado, M. J. May, R. Sundrum, JHEP 0308 (2003) 050, [arXiv:hep-ph/0308036].
-  K. Agashe, R. Contino, L. Da Rold, A. Pomarol, Phys. Lett. B641 (2006) 62-66, [arXiv:hep-ph/0605341].
-  C. Delaunay, O. Gedalia, S. J. Lee, G. Perez, E. Ponton, [arXiv:1101.2902 [hep-ph]].
-  J. A. Cabrer, G. von Gersdorff, M. Quiros, Phys. Lett. B697 (2011) 208-214, [arXiv:1011.2205 [hep-ph]], JHEP 1105 (2011) 083. [arXiv:1103.1388 [hep-ph]].
-  C. Bouchart, G. Moreau, Nucl. Phys. B810 (2009) 66-96. [arXiv:0807.4461 [hep-ph]].