\mathcal{O}(\alpha^{2}) ISR effects with a full electroweak one-loop correction for a top pair-production at the ILC

O(α2) ISR effects with a full electroweak one-loop correction for a top pair-production at the ILC

Junpei Fujimoto High Energy Accelerator Research Organization (KEK), Tsukuba, Ibaraki 305-0801, Japan. The Graduate University for Advanced Studies (SOKENDAI), Hayama, Kanagawa 240-0193, Japan.    Yoshimasa Kurihara111Corresponding author. High Energy Accelerator Research Organization (KEK), Tsukuba, Ibaraki 305-0801, Japan. The Graduate University for Advanced Studies (SOKENDAI), Hayama, Kanagawa 240-0193, Japan.    Nhi M. U. Quach High Energy Accelerator Research Organization (KEK), Tsukuba, Ibaraki 305-0801, Japan. The Graduate University for Advanced Studies (SOKENDAI), Hayama, Kanagawa 240-0193, Japan.
Abstract

Precise predictions for an cross section are presented at an energy region from 400 GeV to 800 GeV. Cross sections are estimated including the beam-polarization effects with full , and also with effects of the initial-state photon emission. A radiator technique is used for the initial-state photon emission up to two-loop order. A weak correction is defined as the full electroweak corrections without the initial-state photonic corrections. As a result, it is obtained that the total cross section of a top quark pair-production receives the weak corrections of over the trivial initial state corrections at a centre of mass energy of 500 GeV. Among the initial state contributions, a contribution from two-loop diagrams gives less than correction over the one-loop ones at the center of mass energies of from GeV to GeV. In addition, an effect of a running coupling constant is also discussed.

- Experiments

1 Introduction

The standard theory of particle physics are established finally by a discovery of the Higgs bosonAad:2012tfa (); Chatrchyan201230 () in 2012. A current target of particle physics is searching for a more fundamental theory beyond the standard model (BSM). A keystone along this direction must be the Higgs boson and top quark. Since the top quark is the heaviest fermion with a mass at the electroweak symmetry-breaking scale, it is naturally expected to have a special role in the BSM.

The international linear colliderBehnke:2013xla () (ILC), which is an electron-positron colliding experiment with centre of mass (CM) energies above 250 GeV, is proposed and intensively discussed as a future project of high-energy physics. One of main goals of ILC experiments is a precise measurement of top quark properties. Detailed Monte Carlo studies have shown that the ILC would be able to measure most of the standard model parameters to within sub-percent levelsBaer:2013cma (). Because of the improvement of the experimental accuracy of the ILC, theoretical predictions are required be given with new level of precision. In particular, a radiative correction due to the electroweak interaction (including spin polarizations) is mandatory for such requirements. Before the discovery of the top quark, a full electroweak radiative correction was conducted for an process at a lower energyFujimoto:1987hu (), and was then obtained independently for higher energiesFleischer:2002rn (); Fleischer2003 (). The same correction fo a process has also been reportedKhiem2013 (). Recently, full electroweak radiative corrections for the process using a narrow-width approximation for the top quarks including the spin-polarization effects are reportedQuach:2017ijt () by authors including those of a present report. A possible search of the minimum SUSY particles trough loop corrections of the top-quark pair-production at the ILC is reported in doi:10.1093/ptep/ptx048 ().

Among several sources of the radiative correction, it is known that the initial state photonic correction gives the largest contribution in general, and thus it is important for a precise estimation of production cross sections. In this report, the precise estimation of the effect due to the initial-state photon-radiations is discussed in detail

2 Calculation Method

2.1 GRACE system

For precise cross-section calculations of the target process, a GRACE-Loop system is used in this study. The GRACE system is an automatic system to calculate cross sections of scattering processes at one-loop level for the standard theoryBelanger2006117 () and the minimum SUSY modelPhysRevD.75.113002 (). The GRACE system has treated electroweak processes with two, three or four particles in the final state are calculatedBelanger2003152 (); Belanger2003163 (); Belanger2003353 (); Kato:2005iw () at the one-loop order. The renormalization of the electroweak interaction is carried out using on-shell schemedoi:10.1143/PTPS.73.1 (); doi:10.1143/PTPS.100.1 (). Infrared divergences are regulated using fictitious photon-massdoi:10.1143/PTPS.100.1 (). The symbolic manipulation package FORMVermaseren:2000nd () is used to handle all Dirac and tensor algebra in -dimensions. For loop integrations, all tensor one-loop integrals are reduced to scalar integrals using our own formalismBelanger2006117 (), then performed integrations using packages FFvanOldenborgh:1990yc () or LoopToolsHahn:1998yk (). Phase-space integrations are done using an adaptive Monte Carlo integration package BASESKawabata1986127 (); KAWABATA1995309 (). For numerical calculations, we use a quartic precision for floating variables.

While using -gauge in linear gauge-fixing terms in the GRACE system,, the non-linear gauge fixing LagrangianBoudjema:1995cb (); Belanger2006117 () is also employing for the sake of the system checking. Before calculating cross sections, we performed numerical tests to confirm that the amplitudes are independent of all redundant parameters around digits at several randomly chosen phase points. In addition to above checks, soft-photon cut-off independence was examined: cross sections at the one-loop level, results must be independent from a head-photon cut-off parameter . We confirmed that, while varying a parameter from GeV to GeV, the results of numerical phase-space integrations are consistent each other within the statistical errors of numerical integrations, that is typically around order.

The effect of the initial photon emission can be factorized when a total energy of emitted photons are small enough compared with a beam energy or a small angle (co-linear) emission . The calculations under a such approximation is referred as to the “soft-colinear photon approximation(SPA)”. Under the SPA, the corrected cross sections with the initial state photon radiation(ISR), , can be obtained from the tree cross sections using a structure function as follows:

 σISR = ∫10dx H(x,s)σTree(s(1−x)), (1)

where is the CM energy square and is a energy fraction of an emitted photon. The structure function can be calculated using the perturbative method with the SPA. Concrete formulae of the structure function are calculated up to two loop orderdoi:10.1143/PTPS.100.1 (). A further improvement of the cross section estimation is possible using the “exponentiation method”. Initial state photon emissions under the SPA, a probability to emit each photon should be independent each other. Thus, the probability to emitt any number of photons can be calculated as;

 p = ∞∑k=0pk=∞∑k=01k!(p1)k, (2)

where is a probability to emit photons. A factor is necessary due to identical particles (photons) appearing in the final state. This is nothing more than the Taylor expansion of the exponential function. Therefore, the effect of the multiple photon emissions can be estimated by putting the one photon emission probability in an argument of the exponential function. This technique is referred as to the exponentiation method. When the exponentiation method is applied to the cross section calculations at loop level, the corrected cross sections can not be expressed simply like the formula (1), because the same loop corrections are included in both of the structure function and loop amplitudes. To avoid a double counting of the same corrections both in the structure function and loop amplitudes, one have to rearrange terms of corrections.

The total cross section at one-loop (fixed) order without the exponentiation, which is denoted as , can be expressed as;

 σNLO;fixed = σLoop+σSoft+σHard+σTree, (3)

where , and are the cross sections from the loop diagram, and soft and hard real-emission corrections, respectively. A photon whose energy is greater (less) than the threshold energy is defined as a hard (soft) photon, respectively. The soft-photon cross section can be expressed as , where a factorized function is obtained from the real-radiation diagrams using the SPA. The SPA consists of three parts, the initial-state, final-state radiations and their interference terms. For the initial state radiation, it can be written as;

 δSPA = απ(2(L−1)logmeλ+12L2−2l(L−1)), (4)

where , and are the fine structure constant, electron mass and fictitious photon mass, respectively. Here two large log-factors appear as and , where is a beam energy and . The threshold energy is included in both and , and then the total cross section must be independent of the value of after summing up all contributions. At the same time, final results are also independent from the photon mass due to the cancellation among contributions in and .

The cross section with fixed order correction (3) can be improved using the exponentiation method. To avoid a double counting of the terms appearing in both loop corrections and the structure function, terms must be re-arranged as follows:

 σNLO;ISR = (σLoop−σTreeδISL)+˜σSoft+σHard+˜σISR, (5)

where is a correction factor from the initial state photon-loop diagrams, that can be given as;

 δISL = 2απ(−(L−1)logmeλ−14L2+34L+π23−1). (6)

This term must be subtracted from because the same contribution is also included in the structure function . A term includes only the final-state radiation and interference terms between the initial and final state radiations. Instead, gives the improved cross section including the initial-state radiation using the radiator methoddoi:10.1143/PTPS.100.1 (). The total cross section can be calculated using the radiator function as;

 σISR = ∫k2c/s0dx1∫1−x10dx2 D(x1,s)D(x2,s)σTree(sx1x2)). (7)

The radiator function , which is corresponding to square root of the structure function, gives a probability to emit a photon with energy fraction of at the CM energy square . In this method, electron and positron can emit different energies, and thus finite boost of the CM system can be treated. The radiator function can be obtained asdoi:10.1143/PTPS.100.1 ().

 D(1−x,s)2 = H(x,s)=Δ2βxβ−1−Δ1β(1−x2) (8) +β28[−4(2−x)logx−1+3(1−x)2xlog(1−x)−2x],

where

 β = 2απ(logsm2e−1), Δ2=1+δ1+δ2, Δ1=1+δ1 δ1=απ(32L+π23−2), δ2=(αLπ)2(−118L+11972−π23).

In this case, the photon mass is cancelled between and , and the threshold energy is cancelled between and . This result is obtained based on perturbative calculations for initial-state photon emission diagrams up to two-loop order doi:10.1143/PTPS.100.1 (). Terms with in (8) are obtained from two-loop diagrams.

3 Results and discussions

3.1 Input parameters

The input parameters used in this report are listed in Table 1. The masses of the light quarks (i.e., other than the top quark) and boson are chosen to be consistent with low-energy experimentsKhiem2015192 (). Other particle masses are taken from recent measurementsOlive:2016xmw (). The weak mixing-angle is given using the on-shell condition because of our renormalization scheme. The fine-structure constant is taken from the low-energy limit of Thomson scattering, again because of our renormalization scheme.

Because the initial state photonic corrections are independent of the beam polarization, all cross section calculations in this report are performed with left (right) polarization of an electron (positron), respectively.

3.2.1 total cross sections

At first, the fixed order correction without using an exponentiation method is investigated. Total cross sections obtained at leading (tree) and next-to-leading order (NLO) calculations are shown in Figure 1. The NLO calculations near the top-quark production threshold (around the CM energy of GeV) shows negative corrections about . This effect is called “coulomb correction”, and it can be large because produced particles are moving slowly and have enough time-duration to interact with each other. The corrections become very small around the CM energy of GeV and increase at the high energy region as to at the CM energy of GeV. Among several types of radiative corrections, e.g., the initial and final state photon radiations, the vertex and box correction, and so on, the initial-state photonic correction give the largest contribution at the high energy region. As shown in Figure 1, the cross sections at tree level including the ISR correction (a dotted line in the figure) are almost the same as the full order electroweak correction (a dashed line in the figure) at CM energies above GeV. This means that the main contribution of the higher order corrections is caused by the initial state photonic corrections. On the other hand, other corrections from the loop diagrams also give a large correction near the threshold region. Around the CM energy GeV, these effects are cancelled accidentally and give a small correction on the total cross section.

After subtracting a trivial ISR correction from the total corrections, one can discuss the “pure” weak correction in the full electroweak radiative corrections. While the NLO correction degree is defined as

 δNLO = σNLO;fixed−σTreeσTree, (9)

the weak correction degree is defined as

 δweak = σNLO;ISR−σISRσISR. (10)

In the definition of , the trivial initial-state photonic corrections are subtracted from the full electroweak radiative corrections, and thus, shows a fraction of the mainly weak-correction in the full electroweak radiative corrections. The behaviors of and are shown in Figure 2 with respect to the CM energies. One can see that the weak correction becomes smaller and smaller at a high energy region, and arrives at almost zero at the CM energy of GeV. On the other hand, at the CM energy of GeV, the pure-weak corrections gives correction over the trivial ISR corrections.

3.2.2 angular distribution

As mentioned above, the radiative correction does not change the total cross section accidentally at the CM energy GeV. While the total cross section stays the same after the radiative correction, an angular distribution is not the case. In general, the real-photon emission affects the angular distribution such that a steep peak being mild. In reality, the ISR correction (a dotted line with rectangle points) makes a forward peak of a top quark production at a tree level (a solid line with triangle points) smaller as shown in Figure 3. In addition to that, the weak correction raises the backward scattering as shown by a dashed line with circle points in Figure 3. As the result, the total cross section does not change so much.

3.3 Photonic correction at two-loop order

The structure function given in (8) are including two-loop effects. All of above results are obtained using full formula (8). As mentioned in previous subsection, a main contribution of the radiative corrections comes from the initial state photonic-correction at the high energy region. Therefore, the ISR correction is important among many terms of the radiative corrections around these energies. If the two-loop contribution has a significant fraction in the full correction, even higher-loop corrections must be considered for the future experiments. The fraction of two-loop contribution over the one-loop one is defined as

 δ2−loop = σISR−σ(1)ISRσ(1)ISR, (11)

where shows the ISR corrected cross sections using the structure function (8) with omitting and terms. Numerical results are shown in Figure 4. The two-loop contribution is smaller than in the energy region between GeV to GeV as shown in Figure 4.

3.4 Running coupling

Yet another improvement of the cross section estimation is known as a running coupling method. This method is also taking a higher order effect of vacuum polarization diagrams into account as an effective coupling constant. After summing up contributions from fermion-loops on gauge-boson propagators, those form solely a gauge-invariant subset, an electro-weak coupling effectively varies according to the four-momentum square of propagators. The effective coupling at the energy scale can be written as;

 α(|q2|) = α(μ2)1−α(μ2)3πlog(|q2|μ2), (12)

where is the four-momentum square of a propagator of the target process. At the same time, the weak mixing angle is obtained from the Fermi weak-coupling constant asAltarelli:473529 ();

 sin2θW = πα(|q2|)√2GFm2W. (13)

The improved cross section can be written as;

 σimp = ∫10dx1∫1−x10dx2 D(x1,s)D(x2,s)σTree(α(sx1x2);sx1x2), (14)

where is the tree cross section at the CM energy-square with the coupling constant . This method is referred to as the improved Born approximation. Numerical results are summarized in Figure 5. Here the measured value of Altarelli:473529 () is used. Above the CM energy of GeV, the improved Born method gives approximated values better than with respect to the NLO cross sections with the ISR.

4 Summary

We calculated the precise cross sections of an process at an energy region from 400 GeV to 800 GeV. Especially, the initial-state photon emissions are discussed in details. An exponentiation technique is applied for the initial-state photon emission up to two-loop order. We found that the total cross section of a top quark pair-production at a center of mass energy of 500 GeV receives the weak corrections of over the trivial ISR corrections. Among the ISR contributions, two-loop diagrams gives less than correction with respect to the one-loop ones at the CM energies from GeV to GeV.

The improved Born approximation gives cross sections better than compared with the NLO cross sections with the ISR.

References

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