Loading [MathJax]/jax/output/HTML-CSS/jax.js

Categorization of two-loop Feynman diagrams in the O(α2) correction to e+eZH

  • The e+eZH process is the dominant process for the Higgs boson production at the future Higgs factory. In order to match the analysis on the Higgs properties with highly precise experiment data, it will be crucial to include the theoretical prediction to the full next-to-next-to-leading order electroweak effect in the production rate σ(e+eZH). In this inspiring work, we categorize the two-loop Feynman diagrams of the O(α2) correction to e+eZH into 6 categories according to relevant topological structures. Although 25377 diagrams contribute to the O(α2) correction in total, the number of the most challenging diagrams with seven denominators is 2250, which contain only 312 non-planar diagrams with 155 independent types. This categorization could be a valuable reference for the complete calculation in future.
  • The discovery of the Higgs boson in 2012 at the Large Hadron Collider (LHC) [1,2] has ushered in a new era in particle physics. In particular, the Higgs boson is regarded as the key in solving some challenging problems such as the problem of hierarchy, the origin of the neutrino mass, and the dark matter problem. Consequently, the precise measurements of the Higgs boson properties have become top priorities for both experimental and theoretical particle physics.

    Although the LHC can produce a lot of Higgs bosons, the enormously complicated QCD backgrounds make sufficiently precise measurements hard to achieve. To precisely measure the properties of the Higgs boson, the next generation e+e colliders have been proposed owing to the Higgs factories aiming at much more accurate measurements. Compared to the LHC, the e+e colliders will have cleaner experiment conditions and higher luminosity. The candidates of the next generation e+e colliders include the Circular Electron Positron Collider (CEPC) [3,4], International Linear Collider (ILC) [5-7], and Future Circular Collider (FCC-ee) [8-10]. All of these are designed to operate at the center-of-mass energy s240250 GeV. In this energy range, the processes to produce Higgs bosons are e+eZH (Higgsstrahlung), e+eνeˉνeH (W fusion), and e+ee+eH (Z fusion). The dominant Higgs production process is the Higgsstrahlung. The recoil mass method can be applied to identify the Higgs boson candidates [3,9,11-13]. Subsequently, the Higgs boson production can be disentangled in a model-independent way.

    At the CEPC, over one million Higgs bosons can be produced in total with an expected integrated luminosity of 5.6 ab1 [14]. With these sizable events, many important properties of the Higgs boson can be measured highly accurately. For instance, the cross section σ(e+eZH) can be measured to the extremely high precision of 0.51% [14]. Since the Higgs boson candidate events can be identified by the recoil mass method, the measurements of the HZZ coupling mainly depend on the precise measurement of σ(e+eZH). Consequently, it is expected that the experiment error in the HZZ coupling may be 0.25% at the CEPC, which is much better than at the HL-LHC [15,16].

    Furthermore, such precise measurements can give the CEPC an unprecedented reach into new physics scenarios that are difficult to probe at the LHC [17]. In the natural supersymmetry (SUSY), typically, the dominant effect on Higgs precision from colored top partners may have blind spots [18,19]. The blind spots can be filled in by the measurement of σ(e+eZH), which is sensitive to loop-level corrections to the tree-level HZZ coupling [17]. When the δσZH approaches 0.2%, further constraints for m˜t1 and m˜t2 can be observed [20]. The HZZ coupling plays an important role in the study of Electroweak Phase Transition (EWPT). In the real scalar singlet model [21-23], the first order phase transition tends to predict a large suppression of the HZZ coupling, ranging from 1% to as much as 30% [24]. With the expected sensitivity of δgHZZ at the CEPC, models with a strong first order phase transition can be tested.

    In addition to the improvement of the experiment accuracy, a more precise theoretical prediction for σ(e+eZH) is also demanded to match the precision of the experiment measurements. The next-to-leading-order (NLO) electroweak (EW) corrections to σ(e+eZH) have been investigated two decades ago [25-27]. The next-to-next-to-leading-order (NNLO) EW-QCD corrections have also been calculated in recent years [28-30]. The results show that the NNLO EW-QCD corrections increase the cross section by more than one percent, which is larger than the expected experiment accuracies of the CEPC. Moreover, it indicates that the NNLO EW corrections can be significant. It is necessary to emphasize that the corrections depend crucially on the renormalization schemes. In α(0) scheme, the NNLO EW-QCD corrections are approximately 1.1% of the leading-order (LO) cross section. However, in the Gμ scheme, the NNLO EW-QCD corrections amount to only 0.3% of the LO cross section [29], and the sensitivity to the different scheme is reduced by NNLO EW-QCD corrections compared to NLO EW corrections. Consequently, the EW-QCD σNNLO ranges from 231 fb to 233 fb.

    Therefore, the missing two-loop corrections to σ(e+eZH) can lead to an intrinsic uncertainty of O(1%) [31], which is still larger than the experiment accuracy. Since σ(e+eZH) is proportional to the square of the HZZ coupling, the theoretical uncertainties also have a significant impact on the extraction of the HZZ coupling. Therefore, the accuracy of the HZZ coupling (0.25%) in the CEPC may not be achieved due to large theoretical uncertainties. Recently, some interesting calculations toward NNLO EW correction have been made, such as in [32]. We are convinced that, if the full NNLO EW corrections to σ(e+eZH) can be calculated, the scheme dependence can be further reduced to stabilize the theoretical prediction.

    Due to the complicacy of the EW interaction, there are more than 20 thousand Feynman diagrams that contribute to the O(α2) correction of e+eZH. The complete calculation of these Feynman diagrams are huge projects. Therefore, in this paper, we focus on the categorization of these Feynman diagrams. This categorization could be helpful for future calculations and analyses. In the next section, we categorize the Feynman diagrams into six categories and numberous subcategories. Finally, the conclusions are presented.

    In the Feynman gauge, we obtained 25377 diagrams contributing to the O(α2) corrections of e+eZH by using FeAmGen which interfaced to Qgraf [33]. We have chosen the Yukawa couplings of light fermions (all fermions except the top quark) to be zero. FeynArts [34] is used to check the correctness to this procedure.

    To categorize the diagrams, we put the diagrams which can be factorized into two one-loop diagrams in the first category C1. Then, according to the number of denominators in each diagram, we categorize the remaining non-factorizable two-loop diagrams into five categories C2,,C6. Furthermore, according to the topologies of loop structures, Ci can be categorized into several subcategories {Ci,j}. Since the light quarks are regarded as massless except for the top quark, we use Ci,j,a (Ci,j,b) to denote the diagrams without/with the top quark. Since some amplitudes can be obtained by replacing coupling factors or masses from other diagrams, Ci,j can be reduced to subset Cindi,j , which only includes the "independent" diagrams. Due to the color structure and conservation laws, 153 diagrams in total have amplitudes equal to zero.

    In this paper we use the Nickel index [35-37] to describe the topologies of loop structures. For the reader's convenience, we briefly explain the Nickel notation and the Nickel index. The Nickel notation is a labeling algorithm to describe connected undirected graphs with "simple" edges and vertices such as the topological structures of the Feynman diagrams. First, one should consider a connected graph with n vertices and label these n vertices by the integers 0 through n1 at random. Therefore, the sequence can be constructed according to [35,36]:

    vertices connected to vertex 0 |vertices connected to 1 excluding 0 |  |vertices connected to vertex i excluding 0 through i1 | |.

    (1)

    For instance, Fig. 1(a) can be represented by 12|223|3|. Otherwise, we can use the label "e" to describe the external lines in the diagrams. Moreover, the Nickel notation of Fig. 1(b) is ee11|ee|. With different labeling strategies, one diagram can be represented as different Nickel notations, which describe the same diagram up to a topological homeomorphism. For simplicity, the Nickel index algorithm can be used to find the "minimal" Nickel notation, which is called the Nickel index. Consequently, the diagram and its Nickel index are in a one-to-one correspondence. For instance, the Nickel index of Fig. 1(a) is 1123|23|||. The package GraphState [35] is a useful tool for constructing the Nickel index. The details of the Nickel index algorithm can be found in Ref. [35].

    Figure 1

    Figure 1.  Nickel notation and Nickel index.

    In this paper, the topological structures of the diagrams with one vertex connecting to one or two external legs are regarded as equivalent. For instance, the topological structures of two diagrams in Fig. 2 can be regarded as equivalent.

    Figure 2

    Figure 2.  Example of equivalent topological structures.

    The category C1 includes 7908 Feynman diagrams that can be factorized into two one-loop diagrams. Therefore, the calculations of diagrams in C1 can be regarded as one-loop level calculations. According to the topologies of loop structures in C1, they are categorized into 3 subcategories.

    The subcategory C1,1 includes 2117 diagrams which contain at least one one-loop vacuum bubble diagram. Furthermore, C1,1,a includes 2055 diagrams, and C1,1,b includes 62 diagrams. In C1,1,b, there are 14 diagrams whose amplitudes are equal to zero. Cind1,1,a has 449 independent diagrams, and Cind1,1,b has 36 independent diagrams. We choose diagram #47 as the representative of C1,1,a and diagram #4418 as the representative of C1,1,b.

    Figure 3

    Figure 3.  Diagram #47 (representative of C1,1,a).

    Figure 4

    Figure 4.  Diagram #4418 (representative of C1,1,b).

    The subcategory C1,2 includes 5513 diagrams that contain self-energy corrections. The diagrams C1,2 do not contain vacuum bubble diagrams. C1,2,a includes 4775 diagrams, and C1,2,b includes 738 diagrams. In C1,2,b, there are 131 diagrams whose amplitudes are equal to zero. Cind1,2,a has 740 independent diagrams, and Cind1,2,b has 278 independent diagrams. We choose diagram #36 as the representative of C1,2,a and diagram #1035 as the representative of C1,2,b.

    Figure 5

    Figure 5.  Diagram #36 (representative of C1,2,a).

    Figure 6

    Figure 6.  Diagram #1035 (representative of C1,2,b).

    The subcategory C1,3 includes 278 diagrams, which contain two vertex corrections. C1,3,a includes 260 diagrams, and C1,3,b includes 18 diagrams. Cind1,3,a has 68 independent diagrams, and Cind1,3,b has 14 independent diagrams. We choose diagram #6983 as the representative of C1,3,a and diagram #23660 as the representative of C1,3,b.

    Figure 7

    Figure 7.  Diagram #6983 (representative of C1,3,a).

    Figure 8

    Figure 8.  Diagram #23660 (representative of C1,3,b).

    The category C2 includes non-factorizable two-loop Feynman diagrams with three denominators. We found that all diagrams in C2 are two-loop self-energy diagrams. C2 includes 18 diagrams, none of which contains a top quark. Cind2 has 8 independent diagrams. We choose diagram #519 as the representative of C2.

    Figure 9

    Figure 9.  Diagram #519 (representative of C2).

    The category C3 includes 593 non-factorizable two-loop Feynman diagrams with four denominators. According to the topologies of loop structures in C3, we categorize them into 3 subcategories.

    The subcategory C3,1 includes 142 diagrams that can be separated into the tree-level diagrams and the two-loop vacuum bubble diagrams. The topology of their loop structures can be noted as 112|2|| in the Nickel index. The calculation of the two-loop vacuum bubble diagram has been well studied [38]. Cind3,1 has 51 independent diagrams. We choose diagram #3961 as the representative of C3,1.

    Figure 10

    Figure 10.  Diagram #3961 (representative of C3,1).

    The subcategory C3,2 includes 337 two-loop self-energy diagrams, none of which contains the top quark. Cind3,2 has 93 independent diagrams. We choose diagram #1 as the representative of C3,2.

    Figure 11

    Figure 11.  Diagram #1 (representative of C3,2).

    The subcategory C3,3 includes 114 two-loop vertex correction diagrams, none of which contains the top quark. The topology of their loop structures can be noted as e112|e2|e| in the Nickel index. The denominators of diagrams in C3,3 only depend on two external momenta. Cind3,3 has 24 independent diagrams. We choose diagram #191 as the representative of C3,3.

    Figure 12

    Figure 12.  Diagram #191 (representative of C3,3).

    The category C4 includes 4773 non-factorizable two-loop Feynman diagrams with five denominators. According to the topologies of loop structures in C4, we categorize them into 3 subcategories.

    The subcategory C4,1 includes 3266 two-loop self-energy diagrams, some of which contain the top quark. C4,1,a includes 2565 diagrams and C4,1,b includes 701 diagrams. Cind4,1,a has 753 independent diagrams, and Cind4,1,b has 249 independent diagrams. We choose diagram #603 as the representative of C4,1,a and diagram #611 as the representative of C4,1,b.

    Figure 13

    Figure 13.  Diagram #603 (representative of C4,1,a).

    Figure 14

    Figure 14.  Diagram #611 (representative of C4,1,b).

    The subcategory C4,2 includes 637 two-loop vertex correction diagrams, none of which contains the top quark. The topology of their loop structures can be noted as e12|e23|3|e| in the Nickel index. The denominators of diagrams in C4,2 only depend on two external momenta. Cind4,2 has 140 independent diagrams. We choose diagram #2676 as the representative of C4,2.

    Figure 15

    Figure 15.  Diagram #2676 (representative of C4,2).

    The subcategory C4,3 includes 870 two-loop vertex correction diagrams, none of which contains the top quark. The topology of their loop structures can be noted as e112|3|e3|e| in the Nickel index. The denominators of diagrams in C4,3 only depend on two external momenta. Cind4,3 has 278 independent diagrams. We choose diagram #3063 as the representative of C4,3.

    Figure 16

    Figure 16.  Diagram #3063 (representative of C4,3).

    The category C5 includes 9835 non-factorizable two-loop Feynman diagrams with six denominators. According to the topologies of loop structures in C5, we categorize them into six subcategories.

    The subcategory C5,1 includes two-loop planar triangle diagrams. The topology of their loop structures can be noted as e12|e3|34|4|e| in the Nickel index. The denominators of diagrams in C5,1 only depend on two external momenta. C5,1 includes 4897 diagrams, some of which contain the top quark. Then, C5,1,a includes 3966 diagrams, and C5,1,b includes 931 diagrams. Cind5,1,a has 1039 independent diagrams, and Cind5,1,b has 397 independent diagrams. We choose diagram #1325 as the representative of C5,1,a and diagram #16206 as the representative of C5,1,b.

    Figure 17

    Figure 17.  Diagram #1325 (representative of C5,1,a).

    Figure 18

    Figure 18.  Diagram #16206 (representative of C5,1,b).

    The subcategory C5,2 includes 184 two-loop planar diagrams, none of which contains the top quark. The topology of their loop structures can be noted as e12|e23|4|e4|e| in the Nickel index. Cind5,2 has 90 independent diagrams. We choose diagram #3613 as the representative of C5,2.

    Figure 19

    Figure 19.  Diagram #3613 (representative of C5,2).

    The subcategory C5,3 includes two-loop planar diagrams. The topology of their loop structures can be noted as e12|e3|e4|44|| in the Nickel index. The denominators of diagrams in C5,3 only depend on two external momenta. C5,3 includes 4067 diagrams, some of which contains the top quark. C5,3,a includes 3260 diagrams, and C5,3,b includes 807 diagrams. In C5,3,b, there are 131 diagrams whose amplitudes are equal to zero. Cind5,3,a has 1077 independent diagrams, and Cind5,3,b has 264 independent diagrams. We choose diagram #14794 as the representative of C5,3,a and diagram #14812 as the representative of C5,3,b.

    Figure 20

    Figure 20.  Diagram #14794 (representative of C5,3,a).

    Figure 21

    Figure 21.  Diagram #14812 (representative of C5,3,b).

    The subcategory C5,4 includes two-loop planar diagrams. The topology of their loop structures can be noted as e112|3|e4|e4|e| in the Nickel index. C5,4 includes 116 Feynman diagrams, none of which contains the top quark. Cind5,4 has 70 independent diagrams. We choose diagram #3845 as a representative of C5,4.

    Figure 22

    Figure 22.  Diagram #3845 (representative of C5,4).

    The subcategory C5,5 includes two-loop non-planar triangle diagrams. The topology of their loop structures can be noted as e12|34|34|e|e| in the Nickel index. The denominators of diagrams in C5,5 only depend on two external momenta. C5,5 includes 560 diagrams, some of which contain the top quark. C5,5,a includes 442 diagrams and C5,5,b includes 118 diagrams. Cind5,5,a has 140 independent diagrams and Cind5,5,b has 54 independent diagrams. We choose diagram #1267 as the representative of C5,5,a and diagram #11100 as the representative of C5,5,b.

    Figure 23

    Figure 23.  Diagram #1267 (representative of C5,5,a).

    Figure 24

    Figure 24.  Diagram #11100 (representative of C5,5,b).

    The subcategory C5,6 includes two-loop non-planar diagrams. The topology of their loop structures can be noted as e12|e34|34|e|e| in the Nickel index. C5,6 includes 11 diagrams, none of which contains the top quark. Cind5,6 has 8 independent diagrams. We choose diagram #3602 as the representative of C5,6.

    Figure 25

    Figure 25.  Diagram #3602 (representative of C5,6).

    The category C6 includes 2250 non-factorizable two-loop Feynman diagrams with seven denominators. According to the topologies of loop structures in C6, we categorize them into 4 subcategories.

    The subcategory C6,1 includes 446 two-loop planar double-box diagrams. The topology of their loop structures can be noted as e12|e3|34|5|e5|e| in the Nickel index. C6,1,a includes 424 diagrams, and C6,1,b includes 22 diagrams. Cind6,1,a has 194 independent diagrams, and Cind6,1,b has 18 independent diagrams. We choose diagram #23202 as the representative of C6,1,a and diagram #23228 as the representative of C6,1,b.

    Figure 26

    Figure 26.  Diagram #23202 (representative of C6,1,a).

    Figure 27

    Figure 27.  Diagram #23228 (representative of C6,1,b).

    The subcategory C6,2 includes 688 two-loop planar diagrams. The topology of their loop structures can be noted as e1|22|3|e4|e5|e6|| in the Nickel index. C6,2,a includes 580 diagrams and C6,2,b includes 108 diagrams. In C6,2,b, there are 4 diagrams whose amplitudes are equal to zero. Cind6,2,a has 299 independent diagrams, and Cind6,2,b has 48 independent diagrams. We choose diagram #24690 as the representative of C6,2,a and diagram #24708 as the representative of C6,2,b.

    Figure 28

    Figure 28.  Diagram #24690 (representative of C6,2,a).

    Figure 29

    Figure 29.  Diagram #24708 (representative of C6,2,b).

    The subcategory C6,3 includes 804 two-loop planar diagrams. The topology of their loop structures can be represented as e12|e3|e4|45|5|e| in the Nickel index. C6,3,a includes 733 diagrams, and C6,2,b includes 71 diagrams. Cind6,3,a has 302 independent diagrams, and Cind6,3,b has 42 independent diagrams. We choose diagram #23886 as the representative of C6,3,a and diagram #23907 as the representative of C6,3,b.

    Figure 30

    Figure 30.  Diagram #23886 (representative of C6,3,a).

    Figure 31

    Figure 31.  Diagram #23907 (representative of C6,3,b).

    The subcategoryC6,4 is the most challenging subcategory which includes 312 two-loop non-planar double-box diagrams. The topology of their loop structures can be noted as e12|e3|45|45|e|e| in the Nickel index. C6,4,a includes 301 diagrams, and C6,4,b includes 11 diagrams. Cind6,4,a has 146 independent diagrams, and Cind6,4,b has 9 independent diagrams. We choose diagram #22890 as the representative of C6,4,a and diagram #22909 as the representative of C6,4,b.

    Figure 32

    Figure 32.  Diagram #22890 (representative of C6,4,a).

    Figure 33

    Figure 33.  Diagram #22909 (representative of C6,4,b).

    Finally, the information for all subcategories has been summarized in Table 1.

    Table 1

    Table 1.  Summary table for all subcategories.
    subcategory name number of diagrams number of denominators non-planar diagrams contains the top quark number of independent diagrams
    C1,1 2117 No Yes 485
    C1,2 5513 No Yes 1018
    C1,3 278 No Yes 82
    C2 18 3 No No 8
    C3,1 142 4 No No 51
    C3,2 337 4 No No 93
    C3,3 114 4 No No 24
    C4,1 3266 5 No Yes 1002
    C4,2 637 5 No No 140
    C4,3 870 5 No No 278
    C5,1 4897 6 No Yes 1436
    C5,2 184 6 No No 90
    C5,3 4067 6 No Yes 1341
    C5,4 116 6 No No 70
    C5,5 560 6 Yes Yes 194
    C5,6 11 6 Yes No 8
    C6,1 446 7 No Yes 212
    C6,2 688 7 No Yes 347
    C6,3 804 7 No Yes 344
    C6,4 312 7 Yes Yes 155
    DownLoad: CSV
    Show Table

    In this paper, we categorize the two-loop Feynman diagrams contributing to the O(α2) corrections in the Higgsstrahlung e+eZH into 6 categories and numerous subcategories. The most challenging subcategory is C6,4, which includes 312 two-loop non-planar double-box diagrams. There are only 155 independent diagrams in C6,4. We hope that the calculations of these Feynman diagrams can be conveniently organized with the help of this categorization.

    The authors would like to thank Ayres Freitas, Hao Liang, and Tao Liu for helpful discussions.

    [1] G. Aad et al., Phys. Lett. B 716, 1-29 (2012 doi: 10.1016/j.physletb.2012.08.020
    [2] S. Chatrchyan et al., Phys. Lett. B 716, 30-61 (2012 doi: 10.1016/j.physletb.2012.08.021
    [3] The CEPC Study Group, CEPC Conceptual Design Report: Volume 2 - Physics & Detector, arXiv: 1811.10545
    [4] The CEPC Study Group, CEPC Conceptual Design Report: Volume 1 - Accelerator, arXiv: 1809.00285
    [5] H. Baer et al., arXiv: 1306.6352
    [6] T. Behnke et al., arXiv: 1306.6327
    [7] P. Bambade et al., arXiv: 1903.01629
    [8] M. Bicer et al., JHEP 01, 164 (2014
    [9] A. Abada et al., Eur. Phys. ST J. 228(2), 261 (2019 doi: 10.1140/epjst/e2019-900045-4
    [10] A. Abada et al., Eur. Phys. J. C 79(6), 474 (2019 doi: 10.1140/epjc/s10052-019-6904-3
    [11] B. Ioffe and V. A. Khoze, Sov. J. Part. Nucl. 9, 50 (1978
    [12] M. McCullough, Phys. Rev. D, 90(1): 015001 (2014), [Erratum: Phys. Rev. D 92, 039903 (2015)]
    [13] J. Yan, S. Watanuki, K. Fujii et al., Phys. Rev. D 94(11), 113002 (2016 doi: 10.1103/PhysRevD.94.113002
    [14] F. An et al., Chin. Phys. C 43(4), 043002 (2019 doi: 10.1088/1674-1137/43/4/043002
    [15] ATLAS Collaboration, Projections for measurements of Higgs boson signal strengths and coupling parameters with the ATLAS detector at a HL-LHC, ATL-PHYS-PUB-2014-016 (2014), http://cds.cern.ch/record/1956710
    [16] Projected performance of an upgraded cms detector at the lhc and hl-lhc: Contribution to the snowmass process, in Community Summer Study 2013, : Snowmass on the Mississippi, arXiv: 1307.7135
    [17] N. Craig, M. Farina, M. McCullough et al., JHEP 03, 146 (2015
    [18] U. Ellwanger, C. Hugonie, and A. M. Teixeira, Phys. Rept. 496, 1-77 (2010 doi: 10.1016/j.physrep.2010.07.001
    [19] J. Fan, M. Reece, and L.-T. Wang, JHEP 08, 152 (2015
    [20] R. Essig, P. Meade, H. Ramani et al., JHEP 09, 085 (2017
    [21] J. Choi and R. Volkas, Phys. Lett. B 317, 385-391 (1993 doi: 10.1016/0370-2693(93)91013-D
    [22] J. McDonald, Phys. Lett. B 323, 339-346 (1994 doi: 10.1016/0370-2693(94)91229-7
    [23] S. Profumo, M. J. Ramsey-Musolf, and G. Shaughnessy, JHEP 08, 010 (2007
    [24] P. Huang, A. J. Long, and L.-T. Wang, Phys. Rev. D 94(7), 075008 (2016 doi: 10.1103/PhysRevD.94.075008
    [25] J. Fleischer and F. Jegerlehner, Nucl. Phys. B 216, 469-492 (1983 doi: 10.1016/0550-3213(83)90296-1
    [26] B. A. Kniehl, Z. Phys. C 55, 605-618 (1992 doi: 10.1007/BF01561297
    [27] A. Denner, J. Kublbeck, R. Mertig et al., Z. Phys. C 56, 261-272 (1992 doi: 10.1007/BF01555523
    [28] Y. Gong, Z. Li, X. Xu et al., Phys. Rev. D 95(9), 093003 (2017 doi: 10.1103/PhysRevD.95.093003
    [29] Q.-F. Sun, F. Feng, Y. Jia et al., Phys. Rev. D 96(5), 051301 (2017 doi: 10.1103/PhysRevD.96.051301
    [30] W. Chen, F. Feng, Y. Jia et al., Chin. Phys. C 43(1), 013108 (2019 doi: 10.1088/1674-1137/43/1/013108
    [31] A. Freitas et al., arXiv: 1906.05379
    [32] Q. Song and A. Freitas, arXiv: 2101.00308
    [33] P. Nogueira, J. Comput. Phys. 105, 279-289 (1993 doi: 10.1006/jcph.1993.1074
    [34] T. Hahn, Comput. Phys. Commun. 140, 418-431 (2001 doi: 10.1016/S0010-4655(01)00290-9
    [35] D. Batkovich, Y. Kirienko, M. Kompaniets et al., arXiv: 1409.8227
    [36] B. Nickel, D. Meiron, and G. B. Jr., Compilation of 2-pt and 4-pt graphs for continuous spin model, University of Guelph report
    [37] J. F. Nagle, J.Math.Phys. 7, 1588-1592 (1966 doi: 10.1063/1.1705069
    [38] A. I. Davydychev and J. Tausk, Nucl. Phys. B 397, 123-142 (1993 doi: 10.1016/0550-3213(93)90338-P
  • [1] G. Aad et al., Phys. Lett. B 716, 1-29 (2012 doi: 10.1016/j.physletb.2012.08.020
    [2] S. Chatrchyan et al., Phys. Lett. B 716, 30-61 (2012 doi: 10.1016/j.physletb.2012.08.021
    [3] The CEPC Study Group, CEPC Conceptual Design Report: Volume 2 - Physics & Detector, arXiv: 1811.10545
    [4] The CEPC Study Group, CEPC Conceptual Design Report: Volume 1 - Accelerator, arXiv: 1809.00285
    [5] H. Baer et al., arXiv: 1306.6352
    [6] T. Behnke et al., arXiv: 1306.6327
    [7] P. Bambade et al., arXiv: 1903.01629
    [8] M. Bicer et al., JHEP 01, 164 (2014
    [9] A. Abada et al., Eur. Phys. ST J. 228(2), 261 (2019 doi: 10.1140/epjst/e2019-900045-4
    [10] A. Abada et al., Eur. Phys. J. C 79(6), 474 (2019 doi: 10.1140/epjc/s10052-019-6904-3
    [11] B. Ioffe and V. A. Khoze, Sov. J. Part. Nucl. 9, 50 (1978
    [12] M. McCullough, Phys. Rev. D, 90(1): 015001 (2014), [Erratum: Phys. Rev. D 92, 039903 (2015)]
    [13] J. Yan, S. Watanuki, K. Fujii et al., Phys. Rev. D 94(11), 113002 (2016 doi: 10.1103/PhysRevD.94.113002
    [14] F. An et al., Chin. Phys. C 43(4), 043002 (2019 doi: 10.1088/1674-1137/43/4/043002
    [15] ATLAS Collaboration, Projections for measurements of Higgs boson signal strengths and coupling parameters with the ATLAS detector at a HL-LHC, ATL-PHYS-PUB-2014-016 (2014), http://cds.cern.ch/record/1956710
    [16] Projected performance of an upgraded cms detector at the lhc and hl-lhc: Contribution to the snowmass process, in Community Summer Study 2013, : Snowmass on the Mississippi, arXiv: 1307.7135
    [17] N. Craig, M. Farina, M. McCullough et al., JHEP 03, 146 (2015
    [18] U. Ellwanger, C. Hugonie, and A. M. Teixeira, Phys. Rept. 496, 1-77 (2010 doi: 10.1016/j.physrep.2010.07.001
    [19] J. Fan, M. Reece, and L.-T. Wang, JHEP 08, 152 (2015
    [20] R. Essig, P. Meade, H. Ramani et al., JHEP 09, 085 (2017
    [21] J. Choi and R. Volkas, Phys. Lett. B 317, 385-391 (1993 doi: 10.1016/0370-2693(93)91013-D
    [22] J. McDonald, Phys. Lett. B 323, 339-346 (1994 doi: 10.1016/0370-2693(94)91229-7
    [23] S. Profumo, M. J. Ramsey-Musolf, and G. Shaughnessy, JHEP 08, 010 (2007
    [24] P. Huang, A. J. Long, and L.-T. Wang, Phys. Rev. D 94(7), 075008 (2016 doi: 10.1103/PhysRevD.94.075008
    [25] J. Fleischer and F. Jegerlehner, Nucl. Phys. B 216, 469-492 (1983 doi: 10.1016/0550-3213(83)90296-1
    [26] B. A. Kniehl, Z. Phys. C 55, 605-618 (1992 doi: 10.1007/BF01561297
    [27] A. Denner, J. Kublbeck, R. Mertig et al., Z. Phys. C 56, 261-272 (1992 doi: 10.1007/BF01555523
    [28] Y. Gong, Z. Li, X. Xu et al., Phys. Rev. D 95(9), 093003 (2017 doi: 10.1103/PhysRevD.95.093003
    [29] Q.-F. Sun, F. Feng, Y. Jia et al., Phys. Rev. D 96(5), 051301 (2017 doi: 10.1103/PhysRevD.96.051301
    [30] W. Chen, F. Feng, Y. Jia et al., Chin. Phys. C 43(1), 013108 (2019 doi: 10.1088/1674-1137/43/1/013108
    [31] A. Freitas et al., arXiv: 1906.05379
    [32] Q. Song and A. Freitas, arXiv: 2101.00308
    [33] P. Nogueira, J. Comput. Phys. 105, 279-289 (1993 doi: 10.1006/jcph.1993.1074
    [34] T. Hahn, Comput. Phys. Commun. 140, 418-431 (2001 doi: 10.1016/S0010-4655(01)00290-9
    [35] D. Batkovich, Y. Kirienko, M. Kompaniets et al., arXiv: 1409.8227
    [36] B. Nickel, D. Meiron, and G. B. Jr., Compilation of 2-pt and 4-pt graphs for continuous spin model, University of Guelph report
    [37] J. F. Nagle, J.Math.Phys. 7, 1588-1592 (1966 doi: 10.1063/1.1705069
    [38] A. I. Davydychev and J. Tausk, Nucl. Phys. B 397, 123-142 (1993 doi: 10.1016/0550-3213(93)90338-P
  • 加载中

Cited by

1. Wu, Q.-F., Li, Z. FEAMGEN.JL: A Julia program for Feynman Amplitude Generation[J]. Computer Physics Communications, 2024. doi: 10.1016/j.cpc.2024.109230
2. Wang, X.-R., Yan, B. Probing the Hgg coupling through the jet charge correlation in Higgs boson decay[J]. Physical Review D, 2023, 108(5): 056010. doi: 10.1103/PhysRevD.108.056010
3. Liu, X., Zhang, J. Quantum chromodynamics studies toward next generation colliders | [面向下一代对撞机的量子色动力学研究][J]. Scientia Sinica: Physica, Mechanica et Astronomica, 2022, 52(7): 270003. doi: 10.1360/SSPMA-2021-0369

Figures(33) / Tables(1)

Get Citation
Zhao Li, Yefan Wang and Quan-feng Wu. Categorization of two-loop Feynman diagrams in the O(α2) correction to e+eZH[J]. Chinese Physics C. doi: 10.1088/1674-1137/abe84d
Zhao Li, Yefan Wang and Quan-feng Wu. Categorization of two-loop Feynman diagrams in the O(α2) correction to e+eZH[J]. Chinese Physics C.  doi: 10.1088/1674-1137/abe84d shu
Milestone
Received: 2002-12-25
Article Metric

Article Views(2606)
PDF Downloads(58)
Cited by(3)
Policy on re-use
To reuse of Open Access content published by CPC, for content published under the terms of the Creative Commons Attribution 3.0 license (“CC CY”), the users don’t need to request permission to copy, distribute and display the final published version of the article and to create derivative works, subject to appropriate attribution.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

Categorization of two-loop Feynman diagrams in the O(α2) correction to e+eZH

    Corresponding author: Zhao Li, zhaoli@ihep.ac.cn
    Corresponding author: Yefan Wang, wangyefan@ihep.ac.cn
    Corresponding author: Quan-feng Wu, wuquanfeng@ihep.ac.cn
  • 1. Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049, China
  • 2. School of Physics Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
  • 3. Center of High Energy Physics, Peking University, Beijing 100871, China

Abstract: The e+eZH process is the dominant process for the Higgs boson production at the future Higgs factory. In order to match the analysis on the Higgs properties with highly precise experiment data, it will be crucial to include the theoretical prediction to the full next-to-next-to-leading order electroweak effect in the production rate σ(e+eZH). In this inspiring work, we categorize the two-loop Feynman diagrams of the O(α2) correction to e+eZH into 6 categories according to relevant topological structures. Although 25377 diagrams contribute to the O(α2) correction in total, the number of the most challenging diagrams with seven denominators is 2250, which contain only 312 non-planar diagrams with 155 independent types. This categorization could be a valuable reference for the complete calculation in future.

    HTML

    I.   INTRODUCTION
    • The discovery of the Higgs boson in 2012 at the Large Hadron Collider (LHC) [1,2] has ushered in a new era in particle physics. In particular, the Higgs boson is regarded as the key in solving some challenging problems such as the problem of hierarchy, the origin of the neutrino mass, and the dark matter problem. Consequently, the precise measurements of the Higgs boson properties have become top priorities for both experimental and theoretical particle physics.

      Although the LHC can produce a lot of Higgs bosons, the enormously complicated QCD backgrounds make sufficiently precise measurements hard to achieve. To precisely measure the properties of the Higgs boson, the next generation e+e colliders have been proposed owing to the Higgs factories aiming at much more accurate measurements. Compared to the LHC, the e+e colliders will have cleaner experiment conditions and higher luminosity. The candidates of the next generation e+e colliders include the Circular Electron Positron Collider (CEPC) [3,4], International Linear Collider (ILC) [5-7], and Future Circular Collider (FCC-ee) [8-10]. All of these are designed to operate at the center-of-mass energy s240250 GeV. In this energy range, the processes to produce Higgs bosons are e+eZH (Higgsstrahlung), e+eνeˉνeH (W fusion), and e+ee+eH (Z fusion). The dominant Higgs production process is the Higgsstrahlung. The recoil mass method can be applied to identify the Higgs boson candidates [3,9,11-13]. Subsequently, the Higgs boson production can be disentangled in a model-independent way.

      At the CEPC, over one million Higgs bosons can be produced in total with an expected integrated luminosity of 5.6 ab1 [14]. With these sizable events, many important properties of the Higgs boson can be measured highly accurately. For instance, the cross section σ(e+eZH) can be measured to the extremely high precision of 0.51% [14]. Since the Higgs boson candidate events can be identified by the recoil mass method, the measurements of the HZZ coupling mainly depend on the precise measurement of σ(e+eZH). Consequently, it is expected that the experiment error in the HZZ coupling may be 0.25% at the CEPC, which is much better than at the HL-LHC [15,16].

      Furthermore, such precise measurements can give the CEPC an unprecedented reach into new physics scenarios that are difficult to probe at the LHC [17]. In the natural supersymmetry (SUSY), typically, the dominant effect on Higgs precision from colored top partners may have blind spots [18,19]. The blind spots can be filled in by the measurement of σ(e+eZH), which is sensitive to loop-level corrections to the tree-level HZZ coupling [17]. When the δσZH approaches 0.2%, further constraints for m˜t1 and m˜t2 can be observed [20]. The HZZ coupling plays an important role in the study of Electroweak Phase Transition (EWPT). In the real scalar singlet model [21-23], the first order phase transition tends to predict a large suppression of the HZZ coupling, ranging from 1% to as much as 30% [24]. With the expected sensitivity of δgHZZ at the CEPC, models with a strong first order phase transition can be tested.

      In addition to the improvement of the experiment accuracy, a more precise theoretical prediction for σ(e+eZH) is also demanded to match the precision of the experiment measurements. The next-to-leading-order (NLO) electroweak (EW) corrections to σ(e+eZH) have been investigated two decades ago [25-27]. The next-to-next-to-leading-order (NNLO) EW-QCD corrections have also been calculated in recent years [28-30]. The results show that the NNLO EW-QCD corrections increase the cross section by more than one percent, which is larger than the expected experiment accuracies of the CEPC. Moreover, it indicates that the NNLO EW corrections can be significant. It is necessary to emphasize that the corrections depend crucially on the renormalization schemes. In α(0) scheme, the NNLO EW-QCD corrections are approximately 1.1% of the leading-order (LO) cross section. However, in the Gμ scheme, the NNLO EW-QCD corrections amount to only 0.3% of the LO cross section [29], and the sensitivity to the different scheme is reduced by NNLO EW-QCD corrections compared to NLO EW corrections. Consequently, the EW-QCD σNNLO ranges from 231 fb to 233 fb.

      Therefore, the missing two-loop corrections to σ(e+eZH) can lead to an intrinsic uncertainty of O(1%) [31], which is still larger than the experiment accuracy. Since σ(e+eZH) is proportional to the square of the HZZ coupling, the theoretical uncertainties also have a significant impact on the extraction of the HZZ coupling. Therefore, the accuracy of the HZZ coupling (0.25%) in the CEPC may not be achieved due to large theoretical uncertainties. Recently, some interesting calculations toward NNLO EW correction have been made, such as in [32]. We are convinced that, if the full NNLO EW corrections to σ(e+eZH) can be calculated, the scheme dependence can be further reduced to stabilize the theoretical prediction.

      Due to the complicacy of the EW interaction, there are more than 20 thousand Feynman diagrams that contribute to the O(α2) correction of e+eZH. The complete calculation of these Feynman diagrams are huge projects. Therefore, in this paper, we focus on the categorization of these Feynman diagrams. This categorization could be helpful for future calculations and analyses. In the next section, we categorize the Feynman diagrams into six categories and numberous subcategories. Finally, the conclusions are presented.

    II.   CATEGORIZATION
    • In the Feynman gauge, we obtained 25377 diagrams contributing to the O(α2) corrections of e+eZH by using FeAmGen which interfaced to Qgraf [33]. We have chosen the Yukawa couplings of light fermions (all fermions except the top quark) to be zero. FeynArts [34] is used to check the correctness to this procedure.

      To categorize the diagrams, we put the diagrams which can be factorized into two one-loop diagrams in the first category C1. Then, according to the number of denominators in each diagram, we categorize the remaining non-factorizable two-loop diagrams into five categories C2,,C6. Furthermore, according to the topologies of loop structures, Ci can be categorized into several subcategories {Ci,j}. Since the light quarks are regarded as massless except for the top quark, we use Ci,j,a (Ci,j,b) to denote the diagrams without/with the top quark. Since some amplitudes can be obtained by replacing coupling factors or masses from other diagrams, Ci,j can be reduced to subset Cindi,j , which only includes the "independent" diagrams. Due to the color structure and conservation laws, 153 diagrams in total have amplitudes equal to zero.

      In this paper we use the Nickel index [35-37] to describe the topologies of loop structures. For the reader's convenience, we briefly explain the Nickel notation and the Nickel index. The Nickel notation is a labeling algorithm to describe connected undirected graphs with "simple" edges and vertices such as the topological structures of the Feynman diagrams. First, one should consider a connected graph with n vertices and label these n vertices by the integers 0 through n1 at random. Therefore, the sequence can be constructed according to [35,36]:

      vertices connected to vertex 0 |vertices connected to 1 excluding 0 |  |vertices connected to vertex i excluding 0 through i1 | |.

      (1)

      For instance, Fig. 1(a) can be represented by 12|223|3|. Otherwise, we can use the label "e" to describe the external lines in the diagrams. Moreover, the Nickel notation of Fig. 1(b) is ee11|ee|. With different labeling strategies, one diagram can be represented as different Nickel notations, which describe the same diagram up to a topological homeomorphism. For simplicity, the Nickel index algorithm can be used to find the "minimal" Nickel notation, which is called the Nickel index. Consequently, the diagram and its Nickel index are in a one-to-one correspondence. For instance, the Nickel index of Fig. 1(a) is 1123|23|||. The package GraphState [35] is a useful tool for constructing the Nickel index. The details of the Nickel index algorithm can be found in Ref. [35].

      Figure 1.  Nickel notation and Nickel index.

      In this paper, the topological structures of the diagrams with one vertex connecting to one or two external legs are regarded as equivalent. For instance, the topological structures of two diagrams in Fig. 2 can be regarded as equivalent.

      Figure 2.  Example of equivalent topological structures.

    • A.   Category C1

    • The category C1 includes 7908 Feynman diagrams that can be factorized into two one-loop diagrams. Therefore, the calculations of diagrams in C1 can be regarded as one-loop level calculations. According to the topologies of loop structures in C1, they are categorized into 3 subcategories.

      The subcategory C1,1 includes 2117 diagrams which contain at least one one-loop vacuum bubble diagram. Furthermore, C1,1,a includes 2055 diagrams, and C1,1,b includes 62 diagrams. In C1,1,b, there are 14 diagrams whose amplitudes are equal to zero. Cind1,1,a has 449 independent diagrams, and Cind1,1,b has 36 independent diagrams. We choose diagram #47 as the representative of C1,1,a and diagram #4418 as the representative of C1,1,b.

      Figure 3.  Diagram #47 (representative of C1,1,a).

      Figure 4.  Diagram #4418 (representative of C1,1,b).

      The subcategory C1,2 includes 5513 diagrams that contain self-energy corrections. The diagrams C1,2 do not contain vacuum bubble diagrams. C1,2,a includes 4775 diagrams, and C1,2,b includes 738 diagrams. In C1,2,b, there are 131 diagrams whose amplitudes are equal to zero. Cind1,2,a has 740 independent diagrams, and Cind1,2,b has 278 independent diagrams. We choose diagram #36 as the representative of C1,2,a and diagram #1035 as the representative of C1,2,b.

      Figure 5.  Diagram #36 (representative of C1,2,a).

      Figure 6.  Diagram #1035 (representative of C1,2,b).

      The subcategory C1,3 includes 278 diagrams, which contain two vertex corrections. C1,3,a includes 260 diagrams, and C1,3,b includes 18 diagrams. Cind1,3,a has 68 independent diagrams, and Cind1,3,b has 14 independent diagrams. We choose diagram #6983 as the representative of C1,3,a and diagram #23660 as the representative of C1,3,b.

      Figure 7.  Diagram #6983 (representative of C1,3,a).

      Figure 8.  Diagram #23660 (representative of C1,3,b).

    • B.   Category C2

    • The category C2 includes non-factorizable two-loop Feynman diagrams with three denominators. We found that all diagrams in C2 are two-loop self-energy diagrams. C2 includes 18 diagrams, none of which contains a top quark. Cind2 has 8 independent diagrams. We choose diagram #519 as the representative of C2.

      Figure 9.  Diagram #519 (representative of C2).

    • C.   Category C3

    • The category C3 includes 593 non-factorizable two-loop Feynman diagrams with four denominators. According to the topologies of loop structures in C3, we categorize them into 3 subcategories.

      The subcategory C3,1 includes 142 diagrams that can be separated into the tree-level diagrams and the two-loop vacuum bubble diagrams. The topology of their loop structures can be noted as 112|2|| in the Nickel index. The calculation of the two-loop vacuum bubble diagram has been well studied [38]. Cind3,1 has 51 independent diagrams. We choose diagram #3961 as the representative of C3,1.

      Figure 10.  Diagram #3961 (representative of C3,1).

      The subcategory C3,2 includes 337 two-loop self-energy diagrams, none of which contains the top quark. Cind3,2 has 93 independent diagrams. We choose diagram #1 as the representative of C3,2.

      Figure 11.  Diagram #1 (representative of C3,2).

      The subcategory C3,3 includes 114 two-loop vertex correction diagrams, none of which contains the top quark. The topology of their loop structures can be noted as e112|e2|e| in the Nickel index. The denominators of diagrams in C3,3 only depend on two external momenta. Cind3,3 has 24 independent diagrams. We choose diagram #191 as the representative of C3,3.

      Figure 12.  Diagram #191 (representative of C3,3).

    • D.   Category C4

    • The category C4 includes 4773 non-factorizable two-loop Feynman diagrams with five denominators. According to the topologies of loop structures in C4, we categorize them into 3 subcategories.

      The subcategory C4,1 includes 3266 two-loop self-energy diagrams, some of which contain the top quark. C4,1,a includes 2565 diagrams and C4,1,b includes 701 diagrams. Cind4,1,a has 753 independent diagrams, and Cind4,1,b has 249 independent diagrams. We choose diagram #603 as the representative of C4,1,a and diagram #611 as the representative of C4,1,b.

      Figure 13.  Diagram #603 (representative of C4,1,a).

      Figure 14.  Diagram #611 (representative of C4,1,b).

      The subcategory C4,2 includes 637 two-loop vertex correction diagrams, none of which contains the top quark. The topology of their loop structures can be noted as e12|e23|3|e| in the Nickel index. The denominators of diagrams in C4,2 only depend on two external momenta. Cind4,2 has 140 independent diagrams. We choose diagram #2676 as the representative of C4,2.

      Figure 15.  Diagram #2676 (representative of C4,2).

      The subcategory C4,3 includes 870 two-loop vertex correction diagrams, none of which contains the top quark. The topology of their loop structures can be noted as e112|3|e3|e| in the Nickel index. The denominators of diagrams in C4,3 only depend on two external momenta. Cind4,3 has 278 independent diagrams. We choose diagram #3063 as the representative of C4,3.

      Figure 16.  Diagram #3063 (representative of C4,3).

    • E.   Category C5

    • The category C5 includes 9835 non-factorizable two-loop Feynman diagrams with six denominators. According to the topologies of loop structures in C5, we categorize them into six subcategories.

      The subcategory C5,1 includes two-loop planar triangle diagrams. The topology of their loop structures can be noted as e12|e3|34|4|e| in the Nickel index. The denominators of diagrams in C5,1 only depend on two external momenta. C5,1 includes 4897 diagrams, some of which contain the top quark. Then, C5,1,a includes 3966 diagrams, and C5,1,b includes 931 diagrams. Cind5,1,a has 1039 independent diagrams, and Cind5,1,b has 397 independent diagrams. We choose diagram #1325 as the representative of C5,1,a and diagram #16206 as the representative of C5,1,b.

      Figure 17.  Diagram #1325 (representative of C5,1,a).

      Figure 18.  Diagram #16206 (representative of C5,1,b).

      The subcategory C5,2 includes 184 two-loop planar diagrams, none of which contains the top quark. The topology of their loop structures can be noted as e12|e23|4|e4|e| in the Nickel index. Cind5,2 has 90 independent diagrams. We choose diagram #3613 as the representative of C5,2.

      Figure 19.  Diagram #3613 (representative of C5,2).

      The subcategory C5,3 includes two-loop planar diagrams. The topology of their loop structures can be noted as e12|e3|e4|44|| in the Nickel index. The denominators of diagrams in C5,3 only depend on two external momenta. C5,3 includes 4067 diagrams, some of which contains the top quark. C5,3,a includes 3260 diagrams, and C5,3,b includes 807 diagrams. In C5,3,b, there are 131 diagrams whose amplitudes are equal to zero. Cind5,3,a has 1077 independent diagrams, and Cind5,3,b has 264 independent diagrams. We choose diagram #14794 as the representative of C5,3,a and diagram #14812 as the representative of C5,3,b.

      Figure 20.  Diagram #14794 (representative of C5,3,a).

      Figure 21.  Diagram #14812 (representative of C5,3,b).

      The subcategory C5,4 includes two-loop planar diagrams. The topology of their loop structures can be noted as e112|3|e4|e4|e| in the Nickel index. C5,4 includes 116 Feynman diagrams, none of which contains the top quark. Cind5,4 has 70 independent diagrams. We choose diagram #3845 as a representative of C5,4.

      Figure 22.  Diagram #3845 (representative of C5,4).

      The subcategory C5,5 includes two-loop non-planar triangle diagrams. The topology of their loop structures can be noted as e12|34|34|e|e| in the Nickel index. The denominators of diagrams in C5,5 only depend on two external momenta. C5,5 includes 560 diagrams, some of which contain the top quark. C5,5,a includes 442 diagrams and C5,5,b includes 118 diagrams. Cind5,5,a has 140 independent diagrams and Cind5,5,b has 54 independent diagrams. We choose diagram #1267 as the representative of C5,5,a and diagram #11100 as the representative of C5,5,b.

      Figure 23.  Diagram #1267 (representative of C5,5,a).

      Figure 24.  Diagram #11100 (representative of C5,5,b).

      The subcategory C5,6 includes two-loop non-planar diagrams. The topology of their loop structures can be noted as e12|e34|34|e|e| in the Nickel index. C5,6 includes 11 diagrams, none of which contains the top quark. Cind5,6 has 8 independent diagrams. We choose diagram #3602 as the representative of C5,6.

      Figure 25.  Diagram #3602 (representative of C5,6).

    • F.   Category C6

    • The category C6 includes 2250 non-factorizable two-loop Feynman diagrams with seven denominators. According to the topologies of loop structures in C6, we categorize them into 4 subcategories.

      The subcategory C6,1 includes 446 two-loop planar double-box diagrams. The topology of their loop structures can be noted as e12|e3|34|5|e5|e| in the Nickel index. C6,1,a includes 424 diagrams, and C6,1,b includes 22 diagrams. Cind6,1,a has 194 independent diagrams, and Cind6,1,b has 18 independent diagrams. We choose diagram #23202 as the representative of C6,1,a and diagram #23228 as the representative of C6,1,b.

      Figure 26.  Diagram #23202 (representative of C6,1,a).

      Figure 27.  Diagram #23228 (representative of C6,1,b).

      The subcategory C6,2 includes 688 two-loop planar diagrams. The topology of their loop structures can be noted as e1|22|3|e4|e5|e6|| in the Nickel index. C6,2,a includes 580 diagrams and C6,2,b includes 108 diagrams. In C6,2,b, there are 4 diagrams whose amplitudes are equal to zero. Cind6,2,a has 299 independent diagrams, and Cind6,2,b has 48 independent diagrams. We choose diagram #24690 as the representative of C6,2,a and diagram #24708 as the representative of C6,2,b.

      Figure 28.  Diagram #24690 (representative of C6,2,a).

      Figure 29.  Diagram #24708 (representative of C6,2,b).

      The subcategory C6,3 includes 804 two-loop planar diagrams. The topology of their loop structures can be represented as e12|e3|e4|45|5|e| in the Nickel index. C6,3,a includes 733 diagrams, and C6,2,b includes 71 diagrams. Cind6,3,a has 302 independent diagrams, and Cind6,3,b has 42 independent diagrams. We choose diagram #23886 as the representative of C6,3,a and diagram #23907 as the representative of C6,3,b.

      Figure 30.  Diagram #23886 (representative of C6,3,a).

      Figure 31.  Diagram #23907 (representative of C6,3,b).

      The subcategoryC6,4 is the most challenging subcategory which includes 312 two-loop non-planar double-box diagrams. The topology of their loop structures can be noted as e12|e3|45|45|e|e| in the Nickel index. C6,4,a includes 301 diagrams, and C6,4,b includes 11 diagrams. Cind6,4,a has 146 independent diagrams, and Cind6,4,b has 9 independent diagrams. We choose diagram #22890 as the representative of C6,4,a and diagram #22909 as the representative of C6,4,b.

      Figure 32.  Diagram #22890 (representative of C6,4,a).

      Figure 33.  Diagram #22909 (representative of C6,4,b).

      Finally, the information for all subcategories has been summarized in Table 1.

      subcategory name number of diagrams number of denominators non-planar diagrams contains the top quark number of independent diagrams
      C1,1 2117 No Yes 485
      C1,2 5513 No Yes 1018
      C1,3 278 No Yes 82
      C2 18 3 No No 8
      C3,1 142 4 No No 51
      C3,2 337 4 No No 93
      C3,3 114 4 No No 24
      C4,1 3266 5 No Yes 1002
      C4,2 637 5 No No 140
      C4,3 870 5 No No 278
      C5,1 4897 6 No Yes 1436
      C5,2 184 6 No No 90
      C5,3 4067 6 No Yes 1341
      C5,4 116 6 No No 70
      C5,5 560 6 Yes Yes 194
      C5,6 11 6 Yes No 8
      C6,1 446 7 No Yes 212
      C6,2 688 7 No Yes 347
      C6,3 804 7 No Yes 344
      C6,4 312 7 Yes Yes 155

      Table 1.  Summary table for all subcategories.

    III.   CONCLUSION
    • In this paper, we categorize the two-loop Feynman diagrams contributing to the O(α2) corrections in the Higgsstrahlung e+eZH into 6 categories and numerous subcategories. The most challenging subcategory is C6,4, which includes 312 two-loop non-planar double-box diagrams. There are only 155 independent diagrams in C6,4. We hope that the calculations of these Feynman diagrams can be conveniently organized with the help of this categorization.

    ACKNOWLEDGMENTS
    • The authors would like to thank Ayres Freitas, Hao Liang, and Tao Liu for helpful discussions.

Reference (38)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return