Identification of extra neutral gauge bosons at the International Linear Collider

Identification of extra neutral gauge bosons at the International Linear Collider

P. Osland,111E-mail: A. A. Pankov,222E-mail: and A. V. Tsytrinov333E-mail:

Heavy neutral gauge bosons, s, are predicted by many theoretical schemes of physics beyond the Standard Model, and intensive searches for their signatures will be performed at present and future high energy colliders. It is quite possible that s are heavy enough to lie beyond the discovery reach expected at the CERN Large Hadron Collider LHC, in which case only indirect signatures of exchanges may occur at future colliders, through deviations of the measured cross sections from the Standard Model predictions. We here discuss in this context the foreseeable sensitivity to s of fermion-pair production cross sections at an linear collider, especially as regards the potential of distinguishing different models once such deviations are observed. Specifically, we assess the discovery and identification reaches on gauge bosons pertinent to the , LR, ALR and SSM classes of models, that should be attained at the planned International Linear Collider (ILC). With the high experimental accuracies expected at the ILC, the discovery and the identification reaches on the models under consideration could be increased substantially. In particular, the identification among the different models could be achieved for values of masses in the discovery (but beyond the identification) reach of the LHC. An important role in enhancing such reaches is played by the electron (and possibly the positron) longitudinally polarized beams. Also, although the purely leptonic processes are experimentally cleaner, the measurements of - and -quark pair production cross sections are found to carry important, and complementary, information on these searches.

Department of Physics and Technology, University of Bergen, Postboks 7803, N-5020 Bergen, Norway

The Abdus Salam ICTP Affiliated Centre, Technical University of Gomel, 246746 Gomel, Belarus

1 Introduction

Electroweak theories beyond the Standard Model (SM) based on spontaneously broken extended gauge symmetries naturally envisage the existence of heavy, neutral, vector bosons . The variety of the proposed models is somewhat broad, and for definiteness in the sequel we shall focus on the so-called , , and models. Particular attention has recently been devoted to the phenomenological properties and the search reaches on such scenarios, and in some sense we may consider these models as representative of this New Physics (NP) sector [1, 2].

A typical manifestation of the production of such states is represented by (narrow) peaks observed in the cross sections for processes among SM particles at high energy accelerators, for example in the invariant mass distributions for Drell-Yan dilepton pair production at the Fermilab Tevatron or at the CERN LHC hadronic colliders. Current experimental search limits on at 95% C.L., from Drell-Yan cross sections at the Tevatron, generally range in the interval 0.8–1 TeV, depending on the particular model being tested [3]. Even higher 95% C.L. limits, of the order of 1.14–1.4 TeV are obtained for the , , and models, from electroweak high precision data [2].

Clearly, the eventual discovery of a peak should be supplemented by the verification of the spin-1 of the assumed underlying , vs. the alternative spin-2 and spin-0 hypotheses corresponding, e.g., to exchanges of a Randall-Sundrum graviton resonance [4] or a sneutrino [5]. This kind of analysis relies on appropriate angular differential distributions and/or angular asymmetries. Finally, once the spin-1 has been established, the particular scenario pertinent to the observed signal should be identified, see, e.g., Refs. [6, 7, 8, 9, 10, 11, 12, 13, 14, 15]. From studies of Drell-Yan processes at the LHC with a time-integrated luminosity of 100 fb, it turns out that one can expect, at the 5- level, discovery limits on of the order of 4–4.5 TeV, spin-1 identification up to –3 TeV and potential of distinction among the individual models up to  TeV (95% C.L.).

For masses above the direct search limits mentioned above, and LHC luminosity at the design value, access to manifestations may be provided by indirect, virtual exchange effects causing deviations of cross sections from the SM predictions, if is not excessively heavy. However, at the LHC, model identification from Drell-Yan dilepton mass distributions and forward-backward asymmetries may be problematic due to limited statistics [16].

An alternative resource for the observation of virtual heavy gauge boson exchanges should be represented by the next generation International Linear Collider (ILC), with center of mass energy –1 TeV and typical time-integrated luminosities –1 ab [17, 18], and the really high precision measurements that will be possible there. Indeed, the baseline configuration envisages a very high electron beam polarization (larger than 80%). Also positron beam polarization, around 30%, might be initially obtainable and perhaps already available for physics. This polarization could be raised to about 60% or higher in the ultimate upgrade of the machine. The polarization option might represent an asset in order to enhance the discovery reaches and identification sensitivities on NP models of any kind [19], therefore also on exchanges in interactions of SM particles. Previous analyses, based on various final state channels and possible experimental observables, show that sensitivities to quite high masses could in principle be attained at the ILC (qualitatively, of the order of for the highest planned luminosity, see, e.g., [1, 20, 21, 22, 23, 24], and references therein). The ILC parameters have recently been fixed in the Reference Design Report [17], so that it should be interesting to reconsider the identification of models in the light of the numbers reported there.

We will here focus on the fermion-antifermion production reactions at the polarized ILC:


As basic experimental observables for the analysis, as an alternative to integrated observables like the total cross sections and/or angular-integrated asymmetries, we here choose the differential angular distributions for the above processes, that allow to exploit the information contained in the different portions of the final state phase space by a binned analysis. Particular emphasis will be given to the comparison between the cases of unpolarized and polarized initial beams, as regards the expected potential of ILC in identifying the models of interest here, for values of the order of and beyond the limits accessible at the LHC.

Indeed, concerning the mass, there are two scenarios. The first one is represented by the interval in between the expected identification and discovery limits at the LHC: here, we can assume the to have already been discovered at some (but the model not identified), so that the model identification (or equivalently the determination of the coupling constants) could be performed at the ILC, based on the deviations of cross sections from the SM predictions for the determined mass. For earlier attempts along this line see, e.g., Ref. [25]. The second mass range is above the LHC discovery limit and, here, with unknown, both discovery and identification reaches should be assessed for the ILC.

In the following, in Sec. 2 we give a brief introduction to the different models considered in the analysis, and give the corresponding leading order expressions of the polarized differential cross sections for processes (1), mostly in order to establish the notations. In Secs. 3 and 4 we present the results of our analysis for the discovery and identification reaches on the individual models at the ILC; and finally, Sec. 5 contains some concluding remarks.

2 Polarized observables and models

The analysis at the ILC is somewhat different from the corresponding studies of Drell-Yan processes at the LHC. Deviations of the various observables from SM predictions, such as cross sections and asymmetries, due to the interference of the SM amplitude with -channel exchanges of the , graviton resonance or sneutrino , might be observed at the ILC. However, in the latter case, there is no interference of with the SM exchanges [5]. Conversely, in the case of the spin-2 KK graviton exchange, the interference with the SM exchanges vanishes when one integrates over the full angular range [26], whereas for differential observables such interference survives. Nevertheless, it turns out [27] that the sensitivity at the ILC with TeV and fb to a KK graviton resonance in processes (1) is of the order 0.8 TeV (1.9 TeV) for the graviton coupling constant (), which is well within the expectations for discovery and identification at the LHC. Accordingly, the KK excitation would have been either discovered or excluded by the time the ILC will be operating. Therefore, from the considerations above, one can conclude that for , as will be the case for the ILC, only the interference of the amplitude with the SM one could be visible at the ILC in processes (1). Accordingly, it might not be so indispensible to perform angular analyses such as those foreseen for Drell-Yan processes at the LHC, in order to differentiate the spin of the exchanged intermediate heavy quantum states, because only should be able to lead to appreciable interference effects.

The polarized differential cross section for the Bhabha process , where and can be exchanged also in the -channel, can be written at leading order as (see, e.g., Refs. [28, 29]):


with the decomposition


In Eqs. (2) and (3), the subscripts and denote helicity cross sections with SM and exchanges in the corresponding channels, and the subscripts , denote the respective helicities, and denote the degrees of longitudinal polarization of the and beams, respectively.111In the recent review [19], the opposite sign convention for positron polarization was adopted. In terms of helicity amplitudes:


According to the previous considerations the amplitudes , with and , are given by the sum of the SM exchanges plus deviations representing the effect induced by a boson:


Here, (we are neglecting fermion masses), and with the electroweak mixing angle, whereas and are characteristic of the particular model. In the annihilation channel, below the mass, the interference with the SM will be destructive in the LL and RR cross sections, whereas it can be of either sign in the LR and RL cross sections.

The polarized differential cross section for the leptonic channels with can be obtained directly from Eq. (2), basically by dropping the -channel contributions. The same is true, after some obvious substitutions, for the annihilations into and final states, in which case also the color () and QCD correction factors, , must be taken into account. The -channel helicity amplitudes for the process (1) with can be written as [28]:


where in the latter expression .

As anticipated, the models that will be considered in our analysis are the following:

  • The scenarios originating from the exceptional group spontaneous breaking are defined in terms of a mixing angle . The specific values , and , correspond to different breaking patterns and define the popular scenarios , and , respectively.

  • The left-right models, originating from the breaking of an grand-unification symmetry, and where the corresponding couples to a combination of right-handed and neutral currents ( and denote baryon and lepton currents), specified by a real parameter bounded by . The particular value corresponds to a pure L-R symmetric model (LRS).

  • The predicted by the ‘alternative’ left-right scenario.

  • The so-called sequential , where the couplings to fermions are the same as those of the SM .

Detailed descriptions of these models, as well as the specific references, can be found, e. g., in Ref. [1]. All numerical values of the couplings needed in Eq. (5) are collected, for example, in Table 1 of Ref. [14].

3 Discovery of

In the absence of available data, the assessment of the expected ‘discovery reaches’ on the various s needs the definition of a ‘distance’ between the NP model predictions and those of the SM for the basic observables that will be measured. The former predictions parametrically depend on the mass and its corresponding coupling constants, while the latter ones are calculated using the parameters known from the SM fits. Such a comparison can be performed by a standard -like procedure. As anticipated in Sec. 1, we divide the full angular range into bins and identify the basic observables with the polarized differential angular distributions for processes (1), , in each bin. Correspondingly, the relevant can symbolically be defined as:


Figure 1: Discovery reaches on models obtained from combined analysis of the unpolarized and polarized processes (1) (95% C.L.) at the ILC with  TeV (1 TeV) and fb (1000 fb), compared to the results expected from Drell-Yan processes at the LHC at the 5- level [14]. Three options of polarization are considered at the ILC: unpolarized beams, ; polarized electron beam, ; both beams polarized, and .

Notice that not only the different beam longitudinal polarizations, but eventually also the various processes ‘’ in Eq. (1) are combined in the definition (7). Here, one assumes to have produced a set of ‘data’, for example by using the dynamics specified by a given model, and in the denominator denotes the corresponding ‘experimental’ uncertainty on , combining statistical and, if possible, systematical ones. According to the previous considerations, the , besides the number of degrees of freedom, is basically a function of the chosen model parameters. In particular, if the coupling constants are fixed at specific values, it will depend solely on the mass, and we vary this parameter. The discovery sensitivity to the under consideration can in this case be identified as the limiting value of for which the value of has the probability needed for exclusion of the SM at a desired confidence level (in what follows, we shall impose 95% C.L.). In the cases where - or -dependent couplings are considered, SM exclusion regions can be defined analogously.

To derive the expected ‘discovery’ limits on models at the ILC, for the ‘annihilation’ channels in Eq. (1), with , we restrict ourselves to combining in Eq. (7) the and beam polarization configurations, that are the predominant ones. For the Bhabha process, , we combine in (7) the cross sections with all four possible polarization configurations, i.e., , ), , . Numerically, following the ILC Design Report [17], we take for the electron beam . For the positron beam, is discussed as possibly available ‘free of charge’ already in the ILC initial running conditions. However, such a small positron polarization will turn out not to affect our evaluated discovery and identification reaches on s considerably. We shall therefore present numerical results for two cases, unpolarized positron beam , and representing the ‘ultimate’ upgrade.

Regarding the ILC energy and the time-integrated luminosity (which, for simplicity, we assume to be equally distributed among the different polarization configurations defined above), still according to Ref. [17], we will give explicit numerical results for c.m. energy TeV with time-integrated luminosity , and for the ‘ultimate’ upgrade values TeV with . The assumed final state identification efficiencies governing, together with the luminosity, the expected statistical uncertainties, are: 100% for pairs; 95% for events (); 35% and 60% for and , respectively [17, 18].

As for the major systematic uncertainties, they originate from errors on beam polarizations, on the time-integrated luminosity, and the final-state reconstruction and energy efficiencies. For the longitudinal polarizations, we adopt the values %, rather ambitious, especially as far as is concerned, but strictly needed for conducting the planned measurements at the permille level, see, e.g., Refs.  [30, 31, 32].222For simplicity we here assume equal precisions on beam polarizations. Clearly, much larger systematic errors on could partially spoil the advantages expected from the availability of also positron beam polarization. As regards the other systematic uncertainties mentioned above, we assume for the combination the (perhaps conservative) lumpsum value of 0.5%. The systematic uncertainties are included using the covariance matrix approach [33, 34, 35].

Figure 2: Relative deviation of the unpolarized differential cross sections from the ALR-model prediction for the process as a function of for the SM and considered models with at the ILC with and . The error bars are statistical uncertainties at the 1- level.

Concerning the theoretical inputs, for the SM amplitudes we use the effective Born approximation [36] vertices, with and . The numerically dominant QED corrections are generated by initial-state radiation, for both Bhabha scattering and the annihilation processes in (1). They are accounted for by a structure function approach including both hard and soft photon emission [37], and by a flux factor method [38], respectively. Effects of radiative flux return to the -channel exchange are minimized by the cut on the radiated photon energy, with . In this way, only interactions that occur close to the nominal collider energy are included in the analysis and, accordingly, the sensitivity to the manifestations of the searched-for nonstandard physics can be optimized. By numerical studies based on the ZFITTER code [39], other QED effects such as final-state and initial-final state emission are found, in the processes () and (), to be numerically unimportant for the chosen kinematical cuts. Finally, correlations between the different polarized cross sections (but not between the individual angular bins) are taken into account in the derivation of the numerical results, that we present in Fig. 1. The figure includes a comparison with the discovery potential of the LHC with luminosity 100 fb, from the Drell-Yan processes () (at the 5- level). These values provide a representative overview of the sensitivities of the reach in on the planned energy and luminosity, as well as on beam polarization.

4 Distinction of models

Basically, in the previous subsection we have assessed the extent to which models can give values of differential cross sections that can exclude the SM hypothesis to a prescribed C.L. Such ‘discovery reaches’ are represented by upper limits on masses, for which the observable deviations between the corresponding models and SM predictions are sufficiently large compared to the foreseeable experimental uncertainties on the cross sections at the ILC.

However, since different models can give rise to similar deviations, we would like to determine the ILC potential of identifying, among the various competing possibilities, the source of a deviation, should it be effectively observed. These ID-limits should obviously be expected to lie below the corresponding ILC discovery reaches and, for an approximate but relatively simple assessment, we adapt the naive -like procedure applied in the previous subsection.

To this purpose, we start by defining a ‘distance’ between pairs of models, and with denoting any of the SSM, SM, ALR, LRS, , , , but . We assume for example model to be the ‘true’ model, namely, we consider ‘data’ sets obtained from the dynamics , with corresponding ‘experimental’ uncertainties, compatible with the expected ‘true’ experimental data. The assessment of its distinguishability from a model, that we call ‘tested’ model, can be performed by a comparison analogous to (7), with the defined as:


As an illustration, the angular behavior of the deviations in the numerator of Eq. (8) for the unpolarized annihilation is depicted in Fig. 2, for the case where the ‘true’ model is ALR, with TeV for all models, at the ILC with TeV and (actually, in this figure, is the relative deviation, ).

Basically, considering that the ILC will start when the LHC will already be operating at the design energy and luminosity, as anticipated previously, we can envisage two cases requiring somewhat different strategies.

4.1 mass known

In the first case we assume that the mass is already measured at the LHC, but perhaps not ‘identified’ there, and the value is within the ILC discovery reaches for both models and . In this case one should set in Eq. (8) and, accordingly, the becomes a function of only the and coupling constants. If both the and couplings are fixed numerically, like in the example of Fig. 2, distinguishability can be assessed by varying , up to the point where the reaches the critical value suitable for exclusion of the ‘tested’ model by the ‘true’ model at the desired confidence level.

If the above mentioned couplings are, instead, the - or - dependent ones, ‘confusion’ domains between ‘true’ and ‘tested’ models can analogously be determined by means of Eq. (8) in the model parameter plane for fixed values of . By definition, in these ‘confusion’ domains, that depend on the actual value assumed for , the cross sections corresponding to definite values of and cannot be distinguished from each other at the desired confidence level. Correspondingly, the ‘complementary’ regions in the above mentioned parameter plane can define the ‘resolution’ domain of the ‘tested’ model by the the ‘true’ model hypothesis, and determine in this way the identification limit on the latter.

Figure 3: Regions of resolution (red hatched) of vs. left-right LR models in the model parameter plane () for obtained with Eq. (8) from the processes with different final states , () , , at , . The role of polarization is demonstrated. Using also a polarized electron (and positron) beam, the resolution region would enlarge to include also the blue hatched (plus the shaded, unhatched) regions, but not the white ones. Those would remain as ‘confusion’ regions.

Figure 4: Similar to Fig. 3 but for .

Figure 5: Similar to Fig. 3 but for .

As an illustrative example of application of the ID-criteria exposed above, we evaluate the ‘resolution’ regions between (‘true’) and left-right LR (‘tested’) models in the plane () for different values of . Figures 3-6 show the regions of ‘resolution’ obtained from the processes , (), and , for 4.5 TeV, 3.5 TeV and 2.5 TeV at TeV and 500 fb, and for different values of beam polarization. Notice that, in these figures, the horizontal axis includes also the values of specific of the , and models, while the vertical axis includes the value of representative of the LRS model.333Actually, we should recall that the point corresponds to the model which, indeed, exists in both classes of models. Therefore, since one would be testing the model against itself, that point must not be considered in an analysis based on Eq. (8).

Figures 35 clearly demonstrate the complementary roles of the processes with different final states, in particular, that the process can potentially be the most efficient one in distinguishing and left-right models from each other (it provides the largest resolution domains). Conversely, the purely leptonic processes, in (1), turn out to determine much less extended ‘resolution’ areas, in particular they cannot discriminate the LR model from the models at and (for any ).

Figure 6: Left: Resolution domains at , exploiting all final states. Right: Corresponding resolution domains for Drell-Yan production at the LHC, for the three values of .

Also, as can be seen from these figures, the leptonic processes are found to provide ‘confusion’ domains (white) located in the ‘central’ part of the plane , around , whereas the processes into final states exhibit the opposite feature. Therefore, as shown in Fig. 6, the combination of all processes is expected to dramatically reduce the ‘confusion’ area in the above mentioned plane and to determine the largest possible domain in which the considered models can be mutually distinguished from one another. The substantial role of electron polarization and, to a somewhat lesser extent, of positron polarization in shrinking ‘confusion’ domains, leading to enlarged model ‘resolution’ domains, can also be seen in these panels. Combining all processes, and with both beams polarized, a ‘confusion’ turns out to persist only in the minute corners shown in Fig. 6, for , and nothing at the lower masses.


Figure 7: Left panel: E6 identification reaches on at 95% C.L. obtained from combination of all unpolarized processes at =0.5 TeV and =500 fb. The model is assumed to be ‘true’ while the others (SSM, LR, ALR, SSM and SM) are taken as tested models. The identification range is indicated as the shaded (yellow) area. Right panel: Similar, but for the polarized processes.


Figure 8: Similar to Fig. 7, for LR models.

It is interesting to compare these resolution regions with the corresponding ones resulting from the assumed discovery in the Drell-Yan process at the LHC, shown in the lower-right panel of Fig. 6. This figure shows that, at the LHC, for the discovery of a 4.5 TeV the corresponding resolution region is found to cover only a narrow strip, and . Even for , at the LHC, there are only modest hyperbola-like strips at where models can be distinguished.

Figure 9: Exclusion reaches on at 95% C.L. obtained from combination of all processes at =0.5 (1.0) TeV and =500 (1000) fb in the case when the ALR model (top panel) or SSM model (bottom) is assumed to be the ‘true’ models while the others are taken as tested models. The unpolarized and polarized cases are compared.

We now continue the above analysis in a somewhat different direction, namely, we wish to determine the limiting value of up to which a particular -dependent model, assumed to be true, can be identified at the ILC in the sense that all the other, potentially competing, models can be excluded. The results from this analysis are shown in Figs. 79.

Figure 7 exhibits the exclusion limits vs.  on the models LR, ALR, SSM and SM (recall that exclusion of the SM determines the discovery reaches), once a model is ‘true’ (all processes combined). For the LR ‘tested’ model, the corresponding curve in Fig. 7 is obtained, for each , by varying in the full allowed range, which gives LR exclusion limits , and choosing the minimum value of such s (in this way the whole class of LR models, as well as the LRS, are excluded). The solid line labelled ‘SM’ represents the discovery reach, i.e., the mass up to which the SM can be excluded. The overall identification range is shown as the shaded (yellow) region. One can see that, in this case, the identification of the class of models considered here is basically determined by the exclusion of the class of the LR models and, for the ‘central’ values of , by the SM (i.e., by the discovery reach). Consequently the ID-limit is, in a (somewhat broad) range around , essentially identical to the discovery limit, whereas it is substantially smaller in the two intervals close to . Figure 7 shows that, numerically, for the ID-limit is as large as , and for near 1 models become more and more difficult to distinguish from the competitor ones. The right panel of Fig. 7 shows the corresponding identification reaches for the polarized case, and , and the quantitative improvements that can be achieved in this case.

Similarly, the identification limits on LR models vs. the parameter can be read off from Fig. 8. The curve labelled as ‘’ is obtained by a procedure analogous to the curve ‘LR’ in Fig. 7, and the solid curve ‘SM’ represents the exclusion limits of the SM (hence the discovery reaches). In this case, the identification of the class of LR s turns out to be determined basically by the exclusion of the class of models, generally not so much below the discovery limit for all values of . On the other hand, the figure shows rather high identification limits, of the order of in the range, say, , whereas they substantially decrease for smaller .

2.5 3.5 4.5 2.5 3.5 4.5 2.5 3.5 4.5 2.5 3.5 4.5
0.5 TeV 49.6 225 785 51.0 241 944 2.3 9.4 27.1 2.1 8.5 24.0
1 TeV 9.1 41.6 125 9.4 42.3 128 0.4 2.0 5.9 0.4 1.8 5.3
Table 1: Required integrated luminosity, , at the two energies and 1 TeV and with polarized beams, required for model identification. Three mass values, , 3.5 TeV and 4.5 TeV, assumed determined at the LHC, are considered.

We can conclude, from Figs. 7 and 8, that the identification reach at the ILC, already at and =500 fb, exceeds the corresponding discovery reach at the LHC. In fact, the full integrated luminosity considered here might be not quite indispensible for this identification. In Table 1 we show the required integrated luminosity, at the two ILC energies of 0.5 and 1 TeV, for the identification of these different models, realized as a at 2.5, 3.5 or 4.5 TeV (within the discovery reach of the LHC).

Finally, in Fig. 9 we summarize the information, of a similar kind as represented in Figs. 7 and 8, relevant to the cases where the ALR model or the SSM model is assumed ‘true’ (upper and lower panels, respectively). As usual, the figure shows the limiting value of at which the other models can be excluded, the SM being one of them. The two energies 0.5 TeV and 1 TeV are considered, the different processes (1) are combined in the , and the unpolarized and polarized cases (, ) are compared. The entries for the LR and models in this figure refer to the worst case, i.e., similar to the procedures adopted for the ‘LR’ and ‘’ curves in Figs. 7 and 8, adopting the lowest value of as and are varied, in order to represent a whole class of models. Specifically, in Fig. 9, one can easily read off the identification reaches for ALR and SSM models at 0.5 TeV: and are both 8–9 TeV if polarization is available. In a sense, the ALR and SSM models are the most ‘orthogonal’ ones since, if either of them is assumed ‘true’, the other one can be excluded up to a really high value of .

4.2 mass not known

Figure 10: Regions of confusion (yellow shaded areas) between a ‘true’ model and the ‘tested’ models in the mass plane () at 95% C.L. obtained from combination of all polarized processes at =0.5 TeV and =500 fb. The dashed lines indicate the identification reach.

The second kind of situation is met in the case where the mass cannot be known a priori, e.g., the is too heavy to be discovered at the LHC [say, 4–5 TeV], but deviations from the SM predictions can still be observed at the ILC. Actually, models with different masses and coupling constants can in principle be the source of a deviation from the SM predictions observed at the ILC. With the coupling constants held fixed numerically at the theoretical values pertinent to the and models under consideration, the of Eq. (8) becomes a function of the two masses, and , both assumed to lie in the respective ILC discovery ranges. In this case, one can derive a contour in the two-dimensional plane where to each value of is associated a value such that, for all , the value of in (8) is consistent with ‘confusion’ of and at the desired confidence level. The region encircled by such a contour will be the ‘confusion’ (or ‘no distinction’) domain between the ‘true’ model and the ‘tested’ model and correspondingly, in the complementary domain the hypothesis could be excluded if is assumed to be ‘true’. We refer to this latter, complementary, region as ‘resolution’ region.

One can iterate this procedure and generate pairwise ‘confusion’ and ‘exclusion’ regions in the two-dimensional planes of parameters for all models . As will be illustrated graphically in the remaining part of the paper, a common feature of such ‘exclusion’ regions is that the relevant contours admit, for each (and obviously fixed ) a minimum value such that, for any value , the ‘tested’ model can be excluded regardless of . We finally assume, as identification limit on the model at ILC, the smallest of the values for , for which all tested models will be excluded by the hypothesis of being ‘true’. Of course, such ID-value of should be smaller (or at most equal), than the ILC discovery reach on model . This procedure can finally be iterated, in turn, to all the different models and the assessment of corresponding ID-reaches. This naive procedure can also be extended in a straightforward way to estimating exclusion ranges—and corresponding identification limits—in the cases where - and/or -dependent models are considered in Eq. (8).

Examples of pairwise ‘confusion’ regions and corresponding contours, relevant to the models chosen in Fig. 1, are shown in Fig. 10. In this figure, the various steps of the procedure outlined above, as well as the final derivation of the ID-limits, can easily be followed. As an example of how to read this figure, consider the hypothesis that the model is ‘true’ (lower left panel), with . Then, if instead the or model should be true, the mass would have to exceed 4.2 or 6.3 TeV, respectively.

Figure 11: Comparison of the -model distinction bounds on obtained from combined analysis of the unpolarized and polarized processes (1) at the ILC with  TeV (1 TeV) and fb (1000 fb), compared to the results expected from Drell-Yan processes at the LHC at 95% C.L. [14]. Two options of polarization are considered: unpolarized beams and both beams are polarized, and .

Finally, Fig. 11 shows the comparison of identification reaches or distinction bounds on the -models considered in Fig. 1, together with the corresponding bounds on obtained from the process at the LHC with c.m. energy 14 TeV and time-integrated luminosity 100 . We assume, for the ILC, the same c.m. energy, luminosty and beam polarization as in Fig. 1. The figure speaks for itself, and in particular clearly exhibits the roles of the ILC parameters. In summary, one might be able to distinguish among the considered models at 95% C.L. up to  TeV (4.0 TeV) for unpolarized (polarized) beams at the ILC (0.5 TeV) and 5.3 TeV (7.0 TeV) at the ILC (1 TeV), respectively. In particular, the figure explicitly manifests the substantial role of electron beam polarization in sharpening the identification reaches. Positron polarization can also give a considerable enhancement in this regard (if measurable with the same high accuracy as for electron polarization), although to a more limited extent in some cases.

Clearly, our analysis is greatly simplified by the fact that the vector and axial vector couplings of the considered s are fixed theoretically. If we wanted to determine them in general, namely, with both masses and coupling constants a priori free variables, the analysis should be five-dimensional with, in addition, the limitation that for (contact-interaction regime), could not be simultaneously extracted. In principle, data at different collider energies could be utilized in this regard, for masses not too far from [40].

5 Concluding remarks

We have explored in some detail how the discovery reach at the ILC depends on the c.m. energy, on the available polarization, as well as on the model actually realized in Nature. The lower part of this range, up to  TeV, will also be covered by the LHC, but the identification reach at the LHC is only up to TeV.

In this LHC discovery range, the cleaner ILC environment, together with the availability of beam polarization, allow for an identification of the particular version realized. Actually, this ILC identification range extends considerably beyond the LHC discovery range. Specifically, the ILC with polarized beams at TeV and 1 TeV allows to identify all considered bosons if . This represents a substantial extension of the the LHC reach.


It is a great pleasure to thank Nello Paver for his contributions to this study, in particular his critical and constructive comments. This research has been partially supported by the Abdus Salam ICTP and the Belarusian Republican Foundation for Fundamental Research. The work of PO has been supported by the Research Council of Norway.


  • [1] For reviews and original references see, e.g., P. Langacker, arXiv:0801.1345 [hep-ph]; T. G. Rizzo, arXiv:hep-ph/0610104;
    A. Leike, Phys. Rept. 317, 143 (1999) [arXiv:hep-ph/9805494];
    J. L. Hewett and T. G. Rizzo, Phys. Rept. 183, 193 (1989).
  • [2] For recent constraints from electroweak data, see J. Erler, P. Langacker, S. Munir and E. R. Pena, JHEP 0908, 017 (2009) [arXiv:0906.2435 [hep-ph]].
  • [3] T. Aaltonen et al. [CDF Collaboration], Phys. Rev. Lett. 99, 171802 (2007) [arXiv:0707.2524 [hep-ex]]; Phys. Rev. Lett. 102, 091805 (2009) [arXiv:0811.0053 [hep-ex]];
    R. J. Hooper [D0 Collaboration], Int. J. Mod. Phys. A 20, 3277 (2005).
  • [4] L. Randall and R. Sundrum, Phys. Rev. Lett. 83 (1999) 3370 [arXiv:hep-ph/9905221]; Phys. Rev. Lett. 83 (1999) 4690 [arXiv:hep-th/9906064].
  • [5] J. Kalinowski, R. Ruckl, H. Spiesberger and P. M. Zerwas, Phys. Lett. B 406 (1997) 314 [arXiv:hep-ph/9703436]; Phys. Lett. B 414 (1997) 297 [arXiv:hep-ph/9708272];
    T. G. Rizzo, Phys. Rev. D 59 (1999) 113004 [arXiv:hep-ph/9811440].
  • [6] B. C. Allanach, K. Odagiri, M. A. Parker and B. R. Webber, JHEP 0009, 019 (2000) [arXiv:hep-ph/0006114];
    B. C. Allanach, K. Odagiri, M. J. Palmer, M. A. Parker, A. Sabetfakhri and B. R. Webber, JHEP 0212, 039 (2002) [arXiv:hep-ph/0211205].
  • [7] M. Dittmar, A. S. Nicollerat and A. Djouadi, Phys. Lett. B 583, 111 (2004) [arXiv:hep-ph/0307020].
  • [8] R. Cousins, J. Mumford, J. Tucker and V. Valuev, JHEP 0511, 046 (2005).
  • [9] S. Godfrey, P. Kalyniak and A. Tomkins, arXiv:hep-ph/0511335; S. Godfrey and T. A. W. Martin, Phys. Rev. Lett. 101, 151803 (2008) [arXiv:0807.1080 [hep-ph]]; R. Diener, S. Godfrey and T. A. W. Martin, arXiv:0910.1334 [hep-ph].
  • [10] D. Feldman, Z. Liu and P. Nath, JHEP 0611, 007 (2006) [arXiv:hep-ph/0606294].
  • [11] F. Petriello and S. Quackenbush, Phys. Rev. D 77, 115004 (2008) [arXiv:0801.4389 [hep-ph]]; Y. Li, F. Petriello and S. Quackenbush, Phys. Rev. D 80, 055018 (2009) [arXiv:0906.4132 [hep-ph]].
  • [12] P. Osland, A. A. Pankov, N. Paver and A. V. Tsytrinov, Phys. Rev. D 78, 035008 (2008) [arXiv:hep-ph/0805.2734], and arXiv:0902.1593 [hep-ph].
  • [13] H. Murayama and V. Rentala, arXiv:0904.4561 [hep-ph].
  • [14] P. Osland, A. A. Pankov, A. V. Tsytrinov and N. Paver, Phys. Rev. D 79, 115021 (2009) [arXiv:0904.4857 [hep-ph]].
  • [15] E. Salvioni, G. Villadoro and F. Zwirner, arXiv:0909.1320 [hep-ph].
  • [16] T. G. Rizzo, JHEP 0908, 082 (2009) [arXiv:0904.2534 [hep-ph]].
  • [17] J. Brau et al. [ILC Collaboration], “ILC Reference Design Report Volume 1 - Executive Summary,” arXiv:0712.1950 [physics.acc-ph].
  • [18] G. Aarons et al. [ILC Collaboration], “International Linear Collider Reference Design Report Volume 2: PHYSICS AT THE ILC,” arXiv:0709.1893 [hep-ph].
  • [19] G. A. Moortgat-Pick et al., Phys. Rept. 460, 131 (2008) [arXiv:hep-ph/0507011].
  • [20] M. Cvetic and S. Godfrey, “Discovery and identification of extra gauge bosons,” arXiv:hep-ph/9504216;
  • [21] S. Riemann, LC report LC-TH-2001-007.
  • [22] A. Djouadi, A. Leike, T. Riemann, D. Schaile and C. Verzegnassi, Z. Phys. C 56, 289 (1992).
  • [23] A. V. Gulov and V. V. Skalozub, Phys. Rev. D 70, 115010 (2004) [arXiv:hep-ph/0408076].
  • [24] A. A. Pankov, N. Paver and A. V. Tsytrinov, Phys. Rev. D 73, 115005 (2006) [arXiv:hep-ph/0512131].
  • [25] G. Weiglein et al. [LHC/LC Study Group], Phys. Rept. 426, 47 (2006) [arXiv:hep-ph/0410364].
  • [26] J. L. Hewett, Phys. Rev. Lett. 82 (1999) 4765 [arXiv:hep-ph/9811356].
  • [27] H. Davoudiasl, J. L. Hewett and T. G. Rizzo, Phys. Rev. Lett. 84 (2000) 2080 [arXiv:hep-ph/9909255].
  • [28] B. Schrempp, F. Schrempp, N. Wermes and D. Zeppenfeld, Nucl. Phys. B 296, 1 (1988).
  • [29] A. A. Pankov and N. Paver, Eur. Phys. J. C 29, 313 (2003) [arXiv:hep-ph/0209058].
  • [30] S. Riemann, A. Schalicke and A. Ushakov, DESY 09-038, arXiv:0903.2366 [physics.ins-det].
  • [31] B. Aurand et al., Report DESY-09-042, arXiv:0903.2959 [physics.acc-ph].
  • [32] S. Boogert et al., ILC-NOTE-2009-049, arXiv:0904.0122 [physics.ins-det].
  • [33] W.T. Eadie, D. Drijard, F.E. James, M. Roos, B. Sadoulet, Statistical methods in experimental physics (American Elsevier, 1971).
  • [34] F. Cuypers and P. Gambino, Phys. Lett. B 388 (1996) 211 [hep-ph/9606391];
    F. Cuypers, hep-ph/9611336.
  • [35] A. A. Babich, P. Osland, A. A. Pankov and N. Paver, Phys. Lett. B 518, 128 (2001) [arXiv:hep-ph/0107159].
  • [36] M. Consoli, W. Hollik and F. Jegerlehner, CERN-TH-5527-89, Presented at Workshop on Z Physics at LEP;
    G. Altarelli, R. Casalbuoni, D. Dominici, F. Feruglio and R. Gatto, Nucl. Phys. B 342, 15 (1990).
  • [37] For reviews see, e.g., O. Nicrosini and L. Trentadue, in Radiative Corrections for Collisions, ed. J. H. Kühn 25 (Springer, Berlin, 1989), p. 25; in QED Structure Functions, Ann Arbor, MI, 1989, ed. G. Bonvicini, AIP Conf. Proc. No. 201 (AIP, New York, 1990), p. 12.
  • [38] For a review see, e.g., W. Benakker and F. A Berends: Proc. of the Workshop on Physics at LEP2, CERN 96-01, vol. 1, p. 79 and references therein.
  • [39] D. Y. Bardin, P. Christova, M. Jack, L. Kalinovskaya, A. Olchevski, S. Riemann and T. Riemann, Comput. Phys. Commun. 133, 229 (2001) [arXiv:hep-ph/9908433].
  • [40] T. G. Rizzo, Phys. Rev. D 55, 5483 (1997) [arXiv:hep-ph/9612304].
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description