Abstract

We present a study to examine the sensitivity of a future collider to the anomalous top flavour-changing neutral current (FCNC) to the gluon. To separate signal from background a multivariate analysis is performed on top quark pair and background events, where one top quark is considered to follow the dominant standard model (SM) decay, , and the other top decays through FCNC, , where is a - or a -quark. The analysis of fully hadronic FCNC decay of the pair is also presented. The 95% confidence level limits on the top quark anomalous couplings are obtained for different values of the center-of-mass energies and integrated luminosities.

1. Introduction

The top quark, which is the heaviest known elementary particle up to now, plays a special role in search for new physics beyond the standard model (SM) in particular through precise measurement of its couplings with other particles. The large mass of the top quark, [1], that is close to the scale of electroweak symmetry breaking and its interactions with other particles such as the Higgs boson make it an excellent object to investigate the validity of the SM. The anomalous interactions of the top quark can occur in various flavour-changing neutral current (FCNC) processes like , where , , , or Higgs. In [2], the anomalous and have been probed at a future electron-positron collider. In the present study, we focus on the top quark FCNC interactions involving the top quark, a light quark , (- or -quark), and a gluon. In the SM, the FCNC transition of () is forbidden at tree level due to the Glashow-Iliopoulos-Maiani (GIM) mechanism [3] and only can proceed through the loop corrections. In the SM framework, the loop-level branching ratio for () is of the order of 10−12 [4, 5]. Clearly, a lot of data is needed to enable us to observe such a decay process and measure this small branching ratio. Various models beyond the SM could lead to a very large increase in FCNC processes involving the top quark. Thus, any evidence of such processes will indicate the existence of new physics. In models beyond SM such as MSSM, Technicolor, extra dimensions models higher branching ratios up to 10−3–10−5 are predicted [610] which can be tested by present high energy experiments. There are several phenomenological studies in search for the anomalous couplings at the Tevatron and LHC and other experiments through different channels [1114]. At present the best and up-to-date experimental limits on the branching fractions come from the direct top production process at the Large Hadron Collider (LHC) by the ATLAS Collaboration, and at a center-of-mass energy of  TeV corresponding to an integrated luminosity of [15, 16].

It is expected that the future scale linear colliders such as Compact Linear Collider (CLIC) or International Linear Collider (ILC) would complete the LHC probes and even in some processes can improve the measurements and limits. The high luminosity and clean experimental environments of the scale collider make it an excellent precision machine for the investigation of the top quark properties. It also provides us with an important opportunity for precise measurements of the FCNC couplings in top quark sector [2, 17]. For example, in [2, 17] it has been shown that the branching ratios of the top quark decay into a photon and a boson can be measured up to the order of 10−6 at a linear electron-positron collider.

The collider, CLIC, is designed to operate with the center-of-mass energies of = 0.5, 1.5, and corresponding to total luminosity of = 2.3, 3.2, and 5.9 10−34 cm−2, respectively [1821]. The design of ILC is to work at the center-of-mass energies from to with the option of upgrading to . The plan for the ILC instantaneous luminosities is to reach 10−33-10−34 cm−2 [22, 23]. One of the main differences between the ILC and CLIC is the difference in the luminosity spectrum (LS) of these machines. The ILC luminosity spectrum has a narrower peak of luminosity. This leads to an increase of the total luminosity and consequently reducing the statistical uncertainty in the measurements.

In this work, we study the sensitivity of a future electron-positron collider (CLIC or ILC) to the anomalous top flavour-changing neutral current (FCNC) to the gluon, . To separate signal from backgrounds, a multivariate technique is used. We consider (semileptonic) and (full-hadronic) separately to search for the anomalous FCNC interactions in vertex. The analysis can also be done in the full-hadronic case with one of the top quarks decays into and the other top decays through anomalous couplings. Because of the large background contribution the results would be better than the semileptonic case; therefore we do not perform the analysis for this decay mode. We consider the center-of-mass energies of 0.5, 1, and , and for these energies we analyse two cases, semileptonic and fully hadronic decays of top quark.

The presented paper is organised as follows. In Section 2, we introduce the theoretical formalism which describes the FCNC processes. Section 3 provides a full detailed description of the semileptonic channel in search for the FCNC. The event selection and the methods of event classification into signal- and background-like events using a multivariate analysis are also discussed in this section. Our fully hadronic analysis is presented in Section 4. The results of the investigated FCNC processes, including expected sensitivities on the anomalous couplings and corresponding branching fractions, are given in Section 5. Discussions on some detector effects and systematic uncertainties are also presented in Section 5. Section 6 contains a summary and conclusions of the analysis.

2. Theoretical Formalism

In this section, we give a brief overview of the theoretical framework for top FCNC which this analysis is based on. In this work, to describe the FCNC couplings amongst the top quark, a light quark, and a gluon () an effective Lagrangian approach is used. The FCNC anomalous interaction in the vertex of can be written as follows [11, 2427]: where the with are dimensionless real parameters that present the strength of the anomalous couplings and strong coupling constant is denoted by . In (1), where are the Gell-Mann matrices, is the new physics scale, is the gluon field tensor, and . In the effective Lagrangian with operators perform the left- (right-) handed projection and are chiral parameters normalized to . In Figure 1, we show the top pair production cross section times branching ratios of one top decays anomalously into and another one decays leptonically (electron and muon) as well as the top pair production cross section times the branching ratios of both tops decay anomalously into . It is presented for different center-of-mass energies, , 1, and , versus the anomalous coupling . As it can be seen the and the for the center-of-mass energy of . In order to calculate the cross section and simulate the events for the analysis, the FCNC effective Lagrangian has been implemented in the FeynRules package [28, 29]; then the model has been imported to a Universal FeynRules Output (UFO) module [30] and finally inserted to the MadGraph 5 [31]. The values of the cross sections are found to be in agreement with CompHEP package [32, 33].

3. Semileptonic Channel

This section presents the analysis of our signal and the corresponding backgrounds of semileptonic channel of events at the collider. In this channel, one of the top quarks decays through SM decay mode of and the other one is considered to decay through FCNC into , where is - or -quark. The hadronic final states of boson have larger background contribution which would not lead to better sensitivity with respect to the semileptonic channel. Therefore, the leptonic decay modes of the boson that provide cleaner signature are considered. The final state signal topology consists of an energetic lepton, neutrino (appears as missing momentum), and three hadronic jets. One of the jets is originated from a -quark. The representative Feynman diagram for the signal process is depicted in Figure 2.

Based on the expected signature of the signal events, the background topology is therefore given by . In order to investigate the possibility of separating signal from the background events, we use Monte Carlo (MC) simulation. The MC generation of the signal sample, , is done with CompHEP [32, 33] and the complete set of backgrounds including the SM process is generated using [31]. The symbol represents any jet that originates from quarks and gluon. To account for the resolution of detectors, we apply energy smearing effects to the final state particles according to the following relations [19, 34]: where and represent the energy of the jets and leptons, respectively. The energies are in GeV and the terms are added in quadrature. The jet energies are smeared according to a Gaussian distribution. We smear the energies of muons similar to the electrons for simplicity. Notice that better resolutions for leptons and jet lead to improvement of the results. We apply the detector acceptance cuts on the transverse momenta of leptons (jets),  GeV, and pseudorapidities, . In order to have well-isolated objects, it is required that the distances in () space between each of the two objects satisfy . It is assumed that the presence of a high- electron or muon plays would be sufficient for triggering the signal events. Now, the signal events are reconstructed as follows. A full reconstruction of the boson four-momentum () is needed to be able to reconstruct the semileptonic decaying top that is the combination of the reconstructed and -jet, . It should yield a distribution consistent with the top quark invariant mass, .

Due to undetected neutrino which leaves no track in the detector, we have difficulties in reconstruction of the boson. The transverse components of the neutrino momentum () can be identified by the missing transverse momentum of the events. The longitudinal component of the momentum of the neutrino, , can be found by solving the following quadratic equation: where we put a mass constraint on , . Solving the above quadratic equation allows us to obtain the longitudinal component of the neutrino momentum. This equation has up to two real solutions. In the case of having two real solutions, the one with minimum absolute value is taken. For the events with complex solution only the real part of the solution is considered as the -component of the neutrino momentum.

We assume a -tagging efficiency of about 60% for -jets, 5% for -quarks, and 1% for lighter quarks to be mistagged as -quark jets [35, 36]. In order to reconstruct the invariant mass of the semileptonic decaying top quark, , we require a completely reconstructed boson and -tagged jet. The anomalously decaying top quark will be reconstructed by combination of the two other remaining jets, which are not tagged as -jets. The reconstructed invariant masses of both top quarks should have mass closest to the physical top quark mass, . In some events, there can be more than one -tagged jet. To make the correct combination of the jets, the event reconstruction is completed by minimizing defined as where is the reconstructed mass of the anomalously decaying top quark and is the reconstructed mass of the top quark decaying through SM. is the angular distance between lepton and jets. It is expected that the jet originating from the semileptonic top quark decays to be too close to the charged lepton. Various combinations of are made and the one with the minimum is chosen. The minimum value of implies that the reconstructed particles fit the requirement of coming from FCNC or SM top quark decay. The reconstructed top quark mass distributions, and , for signal and corresponding background at center-of-mass energy of are shown in Figure 3. In Table 1 (left side), we show the number of signal and background events before and after the kinematical cuts for an integrated luminosity of . In this table, the numbers are presented after including the -tagging efficiency.

Certainly, a detailed background study is essential in order to separate the signal from the background events. We use TMVA [37, 38] as a toolkit for multivariate analysis to separate signal from background. Indeed, a multivariate analysis technique is necessary because a single variable does not have sufficient discrimination power to separate signal from background events. Among the multivariate analysis techniques that are usually used to separate the signal events from the backgrounds, the boosted decision tree (BDT) method is chosen [3941]. We choose the variables which have the most possible separation power between the signals from background events for the BDT input. The kinematical variables are selected as input to the BDT as follow: the reconstructed top quarks masses, the transverse momenta of the -jet and charged lepton difference of the azimuthal angle between the reconstructed boson and the -jet , the pseudorapidity distribution of -jet , the invariant mass of the reconstructed boson , and finally the angular separation between charged lepton and -jet .

The variable list can be extended and more variables could be given in the BDT input for better discrimination between signal and backgrounds. The kinematic distributions of some used variables which have the most discrimination power are presented in Figure 4 before applying the acceptance cuts. These distributions are normalized to unity.

These variables are given to the BDT and the multivariate analysis is performed to achieve the best separation between signal and backgrounds and enhance signal significance. The test and training processes are done using a mixture of 50% of signal and 50% of background events. Due to the sensitivity of the BDT classifier to the statistical fluctuation of the training data sample, we use Adaptive Boosting algorithm to increase the performance. In order to avoid overtraining and to improve the quality of the analysis, the BDT built-in options such as Cost Complexity pruning methods are implemented during the training process. The goal is to find the best cuts that enhance the signal and reduce the background. Obtaining best cuts is generally done by finding the maximum value of the statistical significance, , where is the number of signal events and is the number of background events. By choosing the optimum cuts on the BDT output spectrum, we determine the number of selected signal and background events that provide the best signal significance. The results will be discussed in Section 5.

4. Fully Hadronic Decays of the

In pervious section, we discussed the signal and background where one of the top quarks decays through SM decay mode and the other one is considered to decay through FCNC into . In this section, we consider FCNC decay of both top quarks where the final state consists four jets at the center-of-mass energies of , 1, and . The representative Feynman diagram for the signal process in the full-hadronic channel is depicted in Figure 5.

It is worth mentioning that hadron colliders may not be a good area to study this fully hadronic process due to the extremely large QCD background contributions. The linear electron-positron colliders such as CLIC or ILC have a clean environment; consequently these fully hadronic final states can be probed at the CLIC or ILC easier than the hadron colliders.

The method of the channel is similar to semileptonic one that presented in the previous section. The MC generation of the signal sample is generated with CompHEP and the complete set of four-jet backgrounds is done using MadGraph 5.

Similar to the semileptonic case, to account for the resolution of the detectors a Gaussian energy smearing is performed on the final states jets. The jets are required to have transverse momenta greater than 30 GeV within the pseudorapidity acceptance range of . It is also required that . The number of events before and after the kinematical cuts is shown in the right side of Table 1 for an integrated luminosity of . The reconstructed top quark mass distributions for signal and the corresponding background at center-of-mass energy of are shown in Figure 6 for an integrated luminosity of . The number of signal events in these figures has been multiplied by a factor of 10.

Again we use the boosted decision tree (BDT) classifier of the TMVA package for discriminating signal from background events. For the BDT algorithm, the simulated events of the signal and background are split up into two similar samples for the training and test processes.

The input kinematical variables to the BDT are the reconstructed top quark masses , the transverse momentum of the highest jet, the corresponding pseudorapidity distribution of the highest jet, the angular separation between the two jets, and the scalar transverse energy, . The opening angles between the directions of the final state jets are correlated with the mentioned variables so we neglected them. For the fully hadronic top quark reconstruction, we take the pair of jets which have an invariant mass closest to the nominal top quark mass as well as having smaller angular distance . In the analysis, an angular resolution of around 100 mrad is assumed due to the expected high granularity design of the calorimeters of the future electron-positron collider [42].

In summary, in this section we concentrated on the channel of top pair production where both top quarks decay anomalously into two jets. After a rough detector simulation and applying the acceptance cuts, optimum kinematical variables are found and given to the BDT for discriminating between signal and backgrounds. In the next section, the limit on the branching fractions is given.

5. Results

With assuming of observation no signal events after performing the experiment or in other words if the number of observed events is equal to the number of expected background events, we proceed to set 95% CL upper limit on the signal cross section. Limits on the cross section of signal are calculated with a CLs approach [43]. The RooStats [44] program is used for statistical data analysis for the numerical evaluation of the CLs limits. The program returns the 95% CL upper limit on the signal cross section times branching ratios of the top quarks decays.

The sensitivity of the branching fractions as a function of the integrated luminosity for the future electron-positron collider at different center-of-mass energies is shown in Figures 7 and 8 for the semileptonic and fully hadronic analyses, respectively. As it can be seen from the figure, higher integrated luminosities lead to better bounds on the branching ratio up to around . The limits at the center-of-mass energy of are better than the ones at 1 and because of the larger cross sections at smaller energies. Comparing the semileptonic channel with the full-hadronic one, better sensitivity is achieved in the semileptonic channel. It is again due to the fact that in the full-hadronic channel statistics is poor with respect to the semileptonic one. The upper limits on the at 95% CL with are 0.00117 and 0.0236 for the semileptonic and full-hadronic channels, respectively. It is interesting to mention here that the dependence of the expected upper limit on the integrated luminosity becomes weaker at luminosities larger than .

Now, the sensitivity of the results on the detector performance is discussed. In this analysis, almost all subdetectors are involved to identify and reconstruct leptons, jets, -jets, and missing energy. Precise reconstruction of secondary vertex for an efficient -tagging is necessary in this analysis to suppress the backgrounds and obtain a pure signal sample. The variation of -tagging efficiency in this analysis by 10% leads to approximately 4% change in the expected upper limit on the branching ratio. The resolutions in measurement of the jet and lepton energies are less important than -jet identification. Varying the resolution in jet and lepton energy measurement by 10% and 5% (2) leads to changing the upper limit on the branching fraction by less than 1%.

In this analysis, we have calculated the cross section for the energy at the vertex of electron-positron. Therefore, further effects such as the initial state radiation (ISR) and the luminosity spectrum (LS) of the collider need to be considered. Both the initial state radiation and the luminosity spectrum lead to reducing the cross section. We calculate the effect of ISR on the cross section at the center-of-mass energy of 500 GeV. The signal cross section decreases by around 2% which leads to losing the expected upper limit on the branching ratio from 0.00117 to 0.00118 with .

To have a more realistic analysis, the effects of the systematic uncertainties should be estimated. The uncertainties can arise from jet energy scale, lepton energy, lepton reconstruction identification efficiencies, -tagging efficiency, and uncertainties on the masses of the top quark and boson. We vary the -tagging efficiency by 5%. This leads to changing the expected upper limit by 1.5%. To estimate the uncertainty from jet energy scale, we vary the energy of each jet by 2% and recalculate the limit. It results in a change of 0.5% on the expected upper limit. The uncertainties on the top quark and boson masses are calculated as follows. We generate new signal samples with varied top mass (1 GeV) and boson mass (50 MeV) and redo the analysis. This leads to a change of 0.05% on the expected upper limit on the branching fraction of .

5.1. Comparison with the LHC Results

So far, we have examined the future collider potential to probe the anomalous in the decay of the top quark in top pair production. At hadron colliders, the anomalous couplings can be probed either in top production or top decay. The best limits have been obtained in the production processes. There are different production channels to search for the anomalous : (1) direct top quark production ( process), (2) single top quark production (), (3) double top pair () production, and (4) top plus vector boson production () [45]. Currently, the strongest experimental limits on the branching fractions come from the direct top production ( process) at the LHC by the ATLAS Collaboration, and at a center-of-mass energy of corresponding to an integrated luminosity of .

In [46], the anomalous couplings have been probed in top decay at the Tevatron. The obtained upper limit on the branching ratio is 5(2.7) 10−3 with 10 (30)  of data. These limits are weaker in comparison with the limits that can be obtained from the production processes. The future LHC bounds at center-of-energy using of data using various processes are compared with the ones obtained in this work which are compared in Table 2. As it can be seen, among the all processes the process provides the strongest limit (10−6). The limits that we have obtained in this study for collider are comparable with the ones that come out of the same-sign top () production at the LHC and the ones from top decays in top pair events at the Tevatron.

6. Summary and Conclusions

In this paper, we have studied the signals of top quark flavor-changing neutral current in the vertex of , where and , at a future electron-positron collider. This study has been done by looking at the top pair production and at three different center-of-mass energies of 0.5, 1, and in the top quarks decays. We have investigated two possible cases: first one is the case that one top quark decays anomalously to and another one follows SM decay to a boson and a -quark and boson decays leptonically . Second is that both top quarks decay anomalously through FCNC decay mode . Using the boosted decision tree (BDT) technique, we discriminate between signal and backgrounds. Then the CLs approach has been utilized to set upper limits on the branching ratio. The 95% CL upper limit on the branching ratio using of data at the center-of-mass energy of is 0.00117 (0.0236) in semileptonic (full-hadronic) channel. It is shown that the limit is improved with the integrated luminosity up to around and the dependence of the expected upper limit on the integrated luminosity becomes weaker at luminosities larger than .

We have found that sensitivity to the anomalous couplings decreases with increasing the center of energy of the collisions simply due to the decrease in the signal cross section with growing the center-of-mass energy. The expected bounds are comparable with the ones obtained from the double top production at the LHC () and from the anomalous top decay in top pair events at Tevatron.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors are grateful to R. Goldouzian for valuable discussions. They acknowledge financial support of the School of Particles and Accelerators, Institute for Research in Fundamental Sciences (IPM).