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

Precision Higgs physics at the CEPC

  • The discovery of the Higgs boson with its mass around 125 GeV by the ATLAS and CMS Collaborations marked the beginning of a new era in high energy physics. The Higgs boson will be the subject of extensive studies of the ongoing LHC program. At the same time, lepton collider based Higgs factories have been proposed as a possible next step beyond the LHC, with its main goal to precisely measure the properties of the Higgs boson and probe potential new physics associated with the Higgs boson. The Circular Electron Positron Collider (CEPC) is one of such proposed Higgs factories. The CEPC is an e+e circular collider proposed by and to be hosted in China. Located in a tunnel of approximately 100 km in circumference, it will operate at a center-of-mass energy of 240 GeV as the Higgs factory. In this paper, we present the first estimates on the precision of the Higgs boson property measurements achievable at the CEPC and discuss implications of these measurements.
  • [1] ATLAS Collaboration, Phys. Lett. B, 716: 1-29 (2012), arXiv:1207.7214[hep-ex
    [2] CMS Collaboration, Phys. Lett. B, 716: 30-61 (2012), arXiv:1207.7235[hep-ex
    [3] ATLAS Collaboration, Phys. Lett. B, 726: 88-119 (2013), arXiv:1307.1427[hep-ex
    [4] ATLAS Collaboration, Phys. Lett. B, 726: 120-144 (2013), arXiv:1307.1432[hep-ex
    [5] CMS Collaboration, JHEP, 1306: 081 (2013), arXiv:1303.4571[hep-ex
    [6] CMS Collaboration, Nature Phys., 10: (2014), arXiv:1401.6527[hep-ex
    [7] CMS Collaboration, Phys. Rev. D, 92: 012004 (2015), arXiv:1411.3441[hep-ex
    [8] ATLAS and CMS Collaborations, JHEP, 08: 045 (2016), arXiv:1606.02266[hep-ex
    [9] ATLAS and CMS Collaborations, Phys. Rev. Lett., 114: 191803 (2015), arXiv:1503.07589[hep-ex doi: 10.1103/PhysRevLett.114.191803
    [10] 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
    [11] CMS Collaboration, Projected Performance of an Upgraded CMS Detector at the LHC and HL-LHC: Contribution to the Snowmass Process, arXiv:1307.7135[hep-ex]. http://www.slac.stanford.edu/econf/C1307292/docs/submittedArxivFiles/1307.7135.pdf
    [12] CEPC Study Group, CEPC Conceptual Design Report, arXiv:1809.00285[physics.acc-ph]
    [13] ILD Concept Group-Linear Collider Collaboration, The International Large Detector: Letter of Intent, arXiv:1006.3396[hep-ex]
    [14] T. Behnke, J. E. Brau, P. N. Burrows, J. Fuster, M. Peskin, et al, The International Linear Collider Technical Design Report-Volume 4: Detectors, arXiv:1306.6329[physics.ins-det]
    [15] CEPC-SPPC Study Group, CEPC-SPPC Preliminary Conceptual Design Report. 1. Physics and Detector (2015), http://cepc.ihep.ac.cn/preCDR/volume.html
    [16] P. Janot, Particle Flow Event Reconstruction from LEP to LHC, Presented at Excellence in Detectors and Instrumentation Technologies workshop, CERN (2011), https://indico.cern.ch/event/96989/contribution/15
    [17] M. Minard, Jet energy measurement with the ALEPH detector at LEP2, Presented at CALOR2002 Conference, Pasadena, California, USA (2002), http://inspirehep.net/record/608013
    [18] M. Thomson, Nucl. Instrum. Meth. A, 611: 25 (2009), arXiv:0907.3577[physics.ins-det
    [19] M. Ruan and H. Videau, Arbor, a new approach of the Particle Flow Algorithm, arXiv:1403.4784[physics.ins-det]
    [20] CMS Collaboration, Particle-Flow Event Reconstruction in CMS and Performance for Jets, Taus, and MET, CMS-PAS-PFT-09-001 (2009), http://cds.cern.ch/record/1194487
    [21] F. Beaudette, The CMS Particle Flow Algorithm, pp. 295-304. Proceedings, International Conference on Calorimetry for the High Energy Frontier (CHEF 2013): Paris, France, April 22-25, 2013, arXiv:1401.8155[hep-ex]
    [22] P. Mora de Freitas and H. Videau, Detector simulation with MOKKA/GEANT4: Present and future, Presented at the International Workshop on physics and experiments with future electron-positron linear colliders, Jeju Island, Korea (2002), http://inspirehep.net/record/609687
    [23] S. Agostinelli et al, Nucl. Instrum. Meth. A, 506: 250 (2003)
    [24] M. Ruan et al, Eur. Phys. J. C, 78: 426 (2018)
    [25] H. Zhao, Y.-F. Zhu, C.-D. Fu, D. Yu, and M.-Q. Ruan, The Higgs Signatures at the CEPC CDR Baseline, arXiv:1806.04992[hep-ex]
    [26] D. Yu, M. Ruan, V. Boudry, and H. Videau, Eur. Phys. J. C, 77: 591 (2017), arXiv:1701.07542[physics.ins-det
    [27] S. Catani, Y. L. Dokshitzer, M. Olsson, G. Turnock, and B. Webber, Phys. Lett. B, 269: 432 (1991)
    [28] ATLAS Collaboration, Eur. Phys. J. C, 73: 2306 (2013), arXiv:1210.6210[hep-ex
    [29] CMS Collaboration, JINST, 12: P02014 (2017), arXiv:1607.03663[hep-ex doi: 10.1088/1748-0221/12/02/P02014
    [30] T. Tanabe and T. Suehara, LCFIPlus, Presented at ILD workshop at Kyushu University (2012). https://agenda.linearcollider.org/event/5496/session/1/contribution/16.
    [31] LHC Higgs Cross Section Working Group, Handbook of LHC Higgs Cross Sections: 1. Inclusive Observables, arXiv: 1101.0593[hep-ph]
    [32] LHC Higgs Cross Section Working Group, Handbook of LHC Higgs Cross Sections: 2. Differential Distributions, arXiv: 1201.3084[hep-ph]
    [33] LHC Higgs Cross Section Working Group, Handbook of LHC Higgs Cross Sections: 3. Higgs Properties, arXiv: 1307.1347[hep-ph]
    [34] W. Kilian, T. Ohl, and J. Reuter, Eur. Phys. J. C, 71: 1742 (2011), arXiv:0708.4233[hep-ph
    [35] C. M. Carloni Calame, G. Montagna, O. Nicrosini, and F. Piccinini, Nucl. Phys. Proc. Suppl., 131: 48-55 (2004), arXiv:hep-ph/0312014 doi: 10.1016/j.nuclphysbps.2004.02.008
    [36] S. Aplin, J. Engels, and F. Gaede, A production system for massive data processing in ILCSoft,. http://inspirehep.net/record/889841
    [37] D. Asner, T. Barklow, C. Calancha, K. Fujii, N. Graf, et al., ILC Higgs White Paper, arXiv:1310.0763[hep-ph].
    [38] Y. Haddad, Feasibility of a minimum bias analysis of e+e → ZH → qq+X at a 250 GeV ILC, 010201-44 Chinese Physics C, to be publishing, arXiv:1404.3164[hep-ph]
    [39] M. Oreglia, A Study of the Reactions ψ’→ γγψ, SLAC-R-0236(1980), http://www.slac.stanford.edu/cgi-wrap/getdoc/slac-r-23\6.pdf
    [40] ATLAS Collaboration, Submitted to: Phys. Lett. B, (2018), arXiv:1808.08238[hep-ex
    [41] CMS Collaboration, Submitted to: Phys. Rev. Lett., (2018), arXiv:1808.08242[hep-ex
    [42] R. E. Shrock and M. Suzuki, Phys. Lett. B, 110: 250 (1982)
    [43] K. Griest and H. E. Haber, Phys. Rev. D, 37: 719 (1988)
    [44] C. Englert, T. Plehn, D. Zerwas, and P. M. Zerwas, Phys. Lett. B, 703: 298-305 (2011), arXiv:1106.3097[hep-ph
    [45] C. Bonilla, J. W. F. Valle, and J. C. Romo, Phys. Rev. D, 91: 113015 (2015), arXiv:1502.01649[hep-ph
    [46] ALEPH, Phys. Rept., 532: 119 (2013), arXiv:1302.3415[hep-ex doi: 10.1016/j.physrep.2013.07.004
    [47] S. Dawson, A. Gritsan, H. Logan, J. Qian, C. Tully et al, Working Group Report: Higgs Boson, arXiv:1310.8361[hep-ex]
    [48] LHC Higgs Cross Section Working Group, Handbook of LHC Higgs Cross Sections: 4. Deciphering the Nature of the Higgs Sector, arXiv:1610.07922[hep-ph]
    [49] T. Han, Z. Liu, and J. Sayre, Phys. Rev. D, 89: 113006 (2014), arXiv:1311.7155[hep-ph
    [50] M. E. Peskin, Estimation of LHC and ILC Capabilities for Precision Higgs Boson Coupling Measurements, arXiv:1312.4974[hep-ph]
    [51] P. P. Giardino, K. Kannike, I. Masina, M. Raidal, and A. Strumia, JHEP, 05: 046 (2014), arXiv:1303.3570[hep-ph
    [52] G. Belanger, B. Dumont, U. Ellwanger, J. F. Gunion, and S. Kraml, Phys. Rev. D, 88: 075008 (2013), arXiv:1306.2941[hep-ph
    [53] P. Bechtle, S. Heinemeyer, O. Stal, T. Stefaniak, and G. Weiglein, JHEP, 11: 039 (2014), arXiv:1403.1582[hep-ph
    [54] K. Cheung, J. S. Lee, and P.-Y. Tseng, Phys. Rev. D, 90: 095009 (2014), arXiv:1407.8236[hep-ph
    [55] S. Fichet and G. Moreau, Nucl. Phys. B, 905: 391-446 (2016), arXiv:1509.00472[hep-ph
    [56] R. Lafaye, T. Plehn, M. Rauch, and D. Zerwas, Phys. Rev. D, 96(no. 7): 075044 (2017), arXiv:1706.02174[hep-ph
    [57] V. Barger, T. Han, P. Langacker, B. McElrath, and P. Zerwas, Phys. Rev. D, 67: 115001 (2003), arXiv:hep-ph/0301097
    [58] T. Corbett, O. J. P. Eboli, J. Gonzalez-Fraile, and M. C. Gonzalez-Garcia, Phys. Rev. D, 87: 015022 (2013), arXiv:1211.4580[hep-ph
    [59] J. Elias-Miro, J. R. Espinosa, E. Masso, and A. Pomarol, JHEP, 11: 066 (2013), arXiv:1308.1879[hep-ph
    [60] A. Pomarol and F. Riva, JHEP, 01: 151 (2014), arXiv:1308.2803[hep-ph
    [61] G. Amar, S. Banerjee, S. von Buddenbrock, A. S. Cornell, T. Mandal, B. Mellado, and B. Mukhopadhyaya, JHEP, 02: 128 (2015), arXiv:1405.3957[hep-ph
    [62] J. Ellis, V. Sanz, and T. You, JHEP, 03: 157 (2015), arXiv:1410.7703[hep-ph
    [63] A. Falkowski, Pramana, 87: 39 (2016), arXiv:1505.00046[hep-ph doi: 10.1007/s12043-016-1251-5
    [64] A. Butter, O. J. P. Éboli, J. Gonzalez-Fraile, M. C. Gonzalez-Garcia, T. Plehn, and M. Rauch, JHEP, 07: 152 (2016), arXiv:1604.03105[hep-ph
    [65] N. Craig, M. Farina, M. McCullough, and M. Perelstein, JHEP, 03: 146 (2015), arXiv:1411.0676[hep-ph
    [66] M. Beneke, D. Boito, and Y.-M. Wang, JHEP, 11: 028 (2014), arXiv:1406.1361[hep-ph
    [67] N. Craig, J. Gu, Z. Liu, and K. Wang, JHEP, 03: 050 (2016), arXiv:1512.06877[hep-ph
    [68] J. Ellis and T. You, JHEP, 03: 089 (2016), arXiv:1510.04561[hep-ph
    [69] S.-F. Ge, H.-J. He, and R.-Q. Xiao, JHEP, 10: 007 (2016), arXiv:1603.03385[hep-ph
    [70] J. Ellis, P. Roloff, V. Sanz, and T. You, Dimension-6 Operator Analysis of the CLIC Sensitivity to New Physics, arXiv:1701.04804[hep-ph]
    [71] G. Durieux, C. Grojean, J. Gu, and K. Wang, JHEP, 09: 014 (2017), arXiv:1704.02333[hep-ph
    [72] T. Barklow, K. Fujii, S. Jung, R. Karl, J. List, T. Ogawa, M. E. Peskin, and J. Tian, Phys. Rev. D, 97(no. 5): 053003 (2018), arXiv:1708.08912[hep-ph
    [73] T. Barklow, K. Fujii, S. Jung, M. E. Peskin, and J. Tian, Phys. Rev. D, 97(no. 5): 053004 (2018), arXiv:1708.09079[hep-ph
    [74] S. Di Vita, G. Durieux, C. Grojean, J. Gu, Z. Liu, G. Panico, M. Riembau, and T. Vantalon, JHEP, 02: 178 (2018), arXiv:1711.03978[hep-ph
    [75] W. H. Chiu, S. C. Leung, T. Liu, K.-F. Lyu, and L.-T. Wang, JHEP, 05: 081 (2018), arXiv:1711.04046[hep-ph
    [76] J. Ellis, C. W. Murphy, V. Sanz, and T. You, JHEP, 06: 146 (2018), arXiv:1803.03252[hep-ph
    [77] G. Durieux, J. Gu, E. Vryonidou, and C. Zhang, Chin. Phys. C, 42(no. 12): 123107 (2018), arXiv:1809.03520[hep-ph
    [78] J. Gao, JHEP, 01: 038 (2018), arXiv:1608.01746[hep-ph
    [79] LHC Higgs Cross Section Working Group, LHC HXSWG interim recommendations to explore the coupling structure of a Higgs-like particle, Chinese Physics C, to be publishing, arXiv:1209.0040[hep-ph]
    [80] CMS Collaboration, Phys. Rev. Lett., 120: 231801 (2018), arXiv:1804.02610[hep-ex doi: 10.1103/PhysRevLett.120.231801
    [81] ATLAS Collaboration, Phys. Lett. B, 784: 173-191 (2018), arXiv:1806.00425[hep-ex
    [82] A. Banfi, A. Martin, and V. Sanz, JHEP, 1408: 053 (2014), arXiv:1308.4771[hep-ph
    [83] A. Azatov and A. Paul, JHEP, 1401: 014 (2014), arXiv:1309.5273[hep-ph
    [84] C. Grojean, E. Salvioni, M. Schlaffer, and A. Weiler, JHEP, 1405: 022 (2014), arXiv:1312.3317[hep-ph
    [85] M. Buschmann, C. Englert, D. Goncalves, T. Plehn, and M. Spannowsky, Phys. Rev. D, 90: 013010 (2014), arXiv:1405.7651[hep-ph
    [86] J. Ellis, V. Sanz, and T. You, JHEP, 1407: 036 (2014), arXiv:1404.3667[hep-ph
    [87] M. Klute, R. Lafaye, T. Plehn, M. Rauch, and D. Zerwas, Europhys. Lett., 101: 51001 (2013), arXiv:1301.1322[hep-ph doi: 10.1209/0295-5075/101/51001
    [88] K. Fujii et al, Physics Case for the 250 GeV Stage of the International Linear Collider, arXiv:1710.07621[hep-ex]
    [89] Q.-F. Sun, F. Feng, Y. Jia, and W.-L. Sang, Phys. Rev. D, 96(no. 5): 051301 (2017), arXiv:1609.03995[hep-ph
    [90] Y. Gong, Z. Li, X. Xu, L. L. Yang, and X. Zhao, Phys. Rev. D, 95(no. 9): 093003 (2017), arXiv:1609.03955[hep-ph
    [91] G. P. Lepage, P. B. Mackenzie, and M. E. Peskin, Expected Precision of Higgs Boson Partial Widths within the Standard Model, arXiv:1404.0319[hep-ph]
    [92] A. Denner, S. Heinemeyer, I. Puljak, D. Rebuzzi, and M. Spira, Eur. Phys. J. C, 71: 1753 (2011), arXiv:1107.5909[hep-ph
    [93] L. G. Almeida, S. J. Lee, S. Pokorski, and J. D. Wells, Phys. Rev. D, 89: 033006 (2014), arXiv:1311.6721[hep-ph
    [94] C. Bernaciak, T. Plehn, P. Schichtel, and J. Tattersall, Phys. Rev. D, 91: 035024 (2015), arXiv:1411.7699[hep-ph
    [95] S. M. Barr and A. Zee, Phys. Rev. Lett., 65: 21-24 (1990) doi: 10.1103/PhysRevLett.65.21
    [96] J. Fan and M. Reece, JHEP, 06: 004 (2013), arXiv:1301.2597[hep-ph
    [97] ACME Collaboration, Science, 343: 269-272 (2014), arXiv:1310.7534[physics.atom-ph doi: 10.1126/science.1248213
    [98] Y. T. Chien, V. Cirigliano, W. Dekens, J. de Vries, and E. Mereghetti, JHEP, 02: 011 (2016), arXiv:1510.00725[hep-ph
    [99] R. Harnik, A. Martin, T. Okui, R. Primulando, and F. Yu, Phys. Rev. D, 88: 076009 (2013), arXiv:1308.1094[hep-ph
    [100] K. Hagiwara, S. Ishihara, R. Szalapski, and D. Zeppenfeld, Phys. Rev. D, 48: 2182-2203 (1993)
    [101] G. Gounaris et al, Triple gauge boson couplings, arXiv:hep-ph/9601233[hep-ph]
    [102] L. Bian, J. Shu, and Y. Zhang, JHEP, 09: 206 (2015), arXiv:1507.02238[hep-ph
    [103] A. Falkowski, Higgs Basis: Proposal for an EFT basis choice for LHC HXSWG, LHCHXSWG-INT-2015-001(March, 2015), https://cds.cern.ch/record/2001958
    [104] ATLAS Collaboration, Projections for measurements of Higgs boson cross sections, branching ratios and coupling parameters with the ATLAS detector at a HL-LHC, ATL-PHYS-PUB-2013-014(2013). https://cds.cern.ch/record/1611186
    [105] ATLAS Collaboration, HL-LHC projections for signal and background yield measurements of the H → γγ when the Higgs boson is produced in association with t quarks, W or Z bosons, ATL-PHYS-PUB-2014-012(2014). https://cds.cern.ch/record/1741011
    [106] ATLAS Collaboration, Update of the prospects for the H → Zγ search at the High-Luminosity LHC, ATL-PHYS-PUB-2014-006(2014). https://cds.cern.ch/record/1703276
    [107] ATLAS Collaboration, Prospects for the study of the Higgs boson in the VH(bb) channel at HL-LHC, ATL-PHYS-PUB-2014-011(2014). https://cds.cern.ch/record/1740962
    [108] ATLAS Collaboration, Studies of the VBF H → τlτhad analysis at High Luminosity LHC conditions, ATL-PHYS-PUB-2014-018(2014). https://cds.cern.ch/record/1956732
    [109] A. Falkowski, M. Gonzalez-Alonso, A. Greljo, and D. Marzocca, Phys. Rev. Lett., 116: 011801 (2016), arXiv:1508.00581[hep-ph doi: 10.1103/PhysRevLett.116.011801
    [110] R. Contino, A. Falkowski, F. Goertz, C. Grojean, and F. Riva, JHEP, 07: 144 (2016), arXiv:1604.06444[hep-ph
    [111] A. Falkowski, M. Gonzalez-Alonso, A. Greljo, D. Marzocca, and M. Son, JHEP, 02: 115 (2017), arXiv:1609.06312[hep-ph
    [112] Z. Zhang, Phys. Rev. Lett., 118: 011803 (2017), arXiv:1610.01618[hep-ph doi: 10.1103/PhysRevLett.118.011803
    [113] ATLAS Collaboration, Search for the Standard Model Higgs and Z Boson decays to J/ψ γ: HL-LHC projections, ATL-PHYS-PUB-2015-043(2015). http://cds.cern.ch/record/2054550
    [114] G. T. Bodwin, F. Petriello, S. Stoynev, and M. Velasco, Phys. Rev. D, 88: 053003 (2013), arXiv:1306.5770[hep-ph
    [115] G. Perez, Y. Soreq, E. Stamou, and K. Tobioka, Phys. Rev. D, 92: 033016 (2015), arXiv:1503.00290[hep-ph
    [116] I. Brivio, F. Goertz, and G. Isidori, Phys. Rev. Lett., 115: 211801 (2015), arXiv:1507.02916[hep-ph doi: 10.1103/PhysRevLett.115.211801
    [117] F. Bishara, U. Haisch, P. F. Monni, and E. Re, Constraining Light-Quark Yukawa Couplings from Higgs Distributions, arXiv:1606.09253[hep-ph]
    [118] L. M. Carpenter, T. Han, K. Hendricks, Z. Qian, and N. Zhou, Phys. Rev. D, 95: 053003 (2017), arXiv:1611.05463[hep-ph
    [119] ATLAS Collaboration, Study of the double Higgs production channel H(→ bb)H(→ γγ) with the ATLAS experiment at the HL-LHC, ATL-PHYS-PUB-2017-001, CERN, Geneva, 010201-46 Chinese Physics C, 010201 Jan, 2017. https://cds.cern.ch/record/2243387
    [120] R. Contino et al, CERN Yellow Report,255-440 (2017), arXiv:1606.09408[hep-ph
    [121] M. McCullough, Phys. Rev. D, 90: 015001 (2014), arXiv:1312.3322[hep-ph
    [122] S. Di Vita, C. Grojean, G. Panico, M. Riembau, and T. Vantalon, JHEP, 09: 069 (2017), arXiv:1704.01953[hep-ph
    [123] A. Azatov, R. Contino, G. Panico, and M. Son, Phys. Rev. D, 92: 035001 (2015), arXiv:1502.00539[hep-ph
    [124] Z. Liu, I. Low, and L.-T. Wang, Higgs-Top Interactions at Future Circular e+e Colliders, arXiv:2018.nnnn
    [125] G. Li, H.-R. Wang, and S.-h. Zhu, Probing CP-violating htt coupling in e+e → hγ, arXiv:1506.06453[hep-ph]
    [126] E. Vryonidou and C. Zhang, Dimension-six electroweak top-loop effects in Higgs production and decay, arXiv:1804.09766[hep-ph]
    [127] G. Durieux, M. Perello, M. Vos, and C. Zhang, JHEP, 10: 168 (2018), arXiv:1807.02121[hep-ph
    [128] S. Boselli, R. Hunter, and A. Mitov, Prospects for the determination of the top-quark Yukawa coupling at future e+e colliders, arXiv:1805.12027[hep-ph]
    [129] J. A. Aguilar-Saavedra, Nucl. Phys. B, 812: 181-204 (2009), arXiv:0811.3842[hep-ph
    [130] J. A. Aguilar-Saavedra, Nucl. Phys. B, 821: 215-227 (2009), arXiv:0904.2387[hep-ph
    [131] Y. Gao, A. V. Gritsan, Z. Guo, K. Melnikov, M. Schulze, and N. V. Tran, Phys. Rev. D, 81: 075022 (2010), arXiv:1001.3396[hep-ph
    [132] S. Bolognesi, Y. Gao, A. V. Gritsan, K. Melnikov, M. Schulze, N. V. Tran, and A. Whitbeck, Phys. Rev. D, 86: 095031 (2012), arXiv:1208.4018[hep-ph
    [133] I. Anderson, et al, Phys. Rev. D, 89: 035007 (2014), arXiv:1309.4819[hep-ph
    [134] J. Alwall, P. Demin, S. de Visscher, R. Frederix, M. Herquet, F. Maltoni, T. Plehn, D. L. Rainwater, and T. Stelzer, JHEP, 09: 028 (2007), arXiv:0706.2334[hep-ph
    [135] D. Curtin, R. Essig, S. Gori, P. Jaiswal, A. Katz, et al, Phys.Rev. D, 90: 075004 (2014), arXiv:1312.4992[hep-ph
    [136] Z. Liu, L.-T. Wang, and H. Zhang, Exotic decays of the 125 GeV Higgs boson at future e+e lepton colliders, arXiv:1612.09284[hep-ph]
    [137] J. Liu, Z. Liu, and L.-T. Wang, Long-lived particles at the LHC: catching them in time, arXiv:1805.05957[hep-ph]
    [138] A. Thamm, R. Torre, and A. Wulzer, JHEP, 07: 100 (2015), arXiv:1502.01701[hep-ph
    [139] J. Gu, H. Li, Z. Liu, S. Su, and W. Su, JHEP, 12: 153 (2017), arXiv:1709.06103[hep-ph
    [140] J. Fan, M. Reece, and L.-T. Wang, Precision Natural SUSY at CEPC, FCC-ee, and ILC, arXiv:1412.3107[hep-ph]
    [141] R. Essig, P. Meade, H. Ramani, and Y.-M. Zhong, JHEP, 09: 085 (2017), arXiv:1707.03399[hep-ph
    [142] Z. Chacko, H.-S. Goh, and R. Harnik, Phys. Rev. Lett., 96: 231802 (2006), arXiv:hep-ph/0506256 doi: 10.1103/PhysRevLett.96.231802
    [143] G. Burdman, Z. Chacko, H.-S. Goh, and R. Harnik, JHEP, 02: 009 (2007), arXiv:hep-ph/0609152
    [144] N. Craig, C. Englert, and M. McCullough, Phys. Rev. Lett., 111: 121803 (2013), arXiv:1305.5251[hep-ph doi: 10.1103/PhysRevLett.111.121803
    [145] K. Kajantie, M. Laine, K. Rummukainen, and M. E. Shaposhnikov, Phys. Rev. Lett., 77: 2887-2890 (1996), arXiv:hep-ph/9605288 doi: 10.1103/PhysRevLett.77.2887
    [146] S. Profumo, M. J. Ramsey-Musolf, C. L. Wainwright, and P. Winslow, Phys. Rev. D, 91: 035018 (2015), arXiv:1407.5342[hep-ph
    [147] A. Katz and M. Perelstein, JHEP, 07: 108 (2014), arXiv:1401.1827[hep-ph
    [148] A. Noble and M. Perelstein, Phys. Rev. D, 78: 063518 (2008), arXiv:0711.3018[hep-ph
    [149] B. Henning, X. Lu, and H. Murayama, What do precision Higgs measurements buy us?, arXiv:1404.1058[hep-ph]
    [150] S. Profumo, M. J. Ramsey-Musolf, and G. Shaughnessy, JHEP, 08: 010 (2007), arXiv:0705.2425[hep-ph
    [151] D. Curtin, P. Meade, and C.-T. Yu, JHEP, 11: 127 (2014), arXiv:1409.0005[hep-ph
    [152] Q.-H. Cao, F. P. Huang, K.-P. Xie, and X. Zhang, Chin. Phys. C, 42: 023103 (2018), arXiv:1708.04737[hep-ph
    [153] P. Huang, A. J. Long, and L.-T. Wang, Phys. Rev. D, 94: 075008 (2016), arXiv:1608.06619[hep-ph
  • [1] ATLAS Collaboration, Phys. Lett. B, 716: 1-29 (2012), arXiv:1207.7214[hep-ex
    [2] CMS Collaboration, Phys. Lett. B, 716: 30-61 (2012), arXiv:1207.7235[hep-ex
    [3] ATLAS Collaboration, Phys. Lett. B, 726: 88-119 (2013), arXiv:1307.1427[hep-ex
    [4] ATLAS Collaboration, Phys. Lett. B, 726: 120-144 (2013), arXiv:1307.1432[hep-ex
    [5] CMS Collaboration, JHEP, 1306: 081 (2013), arXiv:1303.4571[hep-ex
    [6] CMS Collaboration, Nature Phys., 10: (2014), arXiv:1401.6527[hep-ex
    [7] CMS Collaboration, Phys. Rev. D, 92: 012004 (2015), arXiv:1411.3441[hep-ex
    [8] ATLAS and CMS Collaborations, JHEP, 08: 045 (2016), arXiv:1606.02266[hep-ex
    [9] ATLAS and CMS Collaborations, Phys. Rev. Lett., 114: 191803 (2015), arXiv:1503.07589[hep-ex doi: 10.1103/PhysRevLett.114.191803
    [10] 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
    [11] CMS Collaboration, Projected Performance of an Upgraded CMS Detector at the LHC and HL-LHC: Contribution to the Snowmass Process, arXiv:1307.7135[hep-ex]. http://www.slac.stanford.edu/econf/C1307292/docs/submittedArxivFiles/1307.7135.pdf
    [12] CEPC Study Group, CEPC Conceptual Design Report, arXiv:1809.00285[physics.acc-ph]
    [13] ILD Concept Group-Linear Collider Collaboration, The International Large Detector: Letter of Intent, arXiv:1006.3396[hep-ex]
    [14] T. Behnke, J. E. Brau, P. N. Burrows, J. Fuster, M. Peskin, et al, The International Linear Collider Technical Design Report-Volume 4: Detectors, arXiv:1306.6329[physics.ins-det]
    [15] CEPC-SPPC Study Group, CEPC-SPPC Preliminary Conceptual Design Report. 1. Physics and Detector (2015), http://cepc.ihep.ac.cn/preCDR/volume.html
    [16] P. Janot, Particle Flow Event Reconstruction from LEP to LHC, Presented at Excellence in Detectors and Instrumentation Technologies workshop, CERN (2011), https://indico.cern.ch/event/96989/contribution/15
    [17] M. Minard, Jet energy measurement with the ALEPH detector at LEP2, Presented at CALOR2002 Conference, Pasadena, California, USA (2002), http://inspirehep.net/record/608013
    [18] M. Thomson, Nucl. Instrum. Meth. A, 611: 25 (2009), arXiv:0907.3577[physics.ins-det
    [19] M. Ruan and H. Videau, Arbor, a new approach of the Particle Flow Algorithm, arXiv:1403.4784[physics.ins-det]
    [20] CMS Collaboration, Particle-Flow Event Reconstruction in CMS and Performance for Jets, Taus, and MET, CMS-PAS-PFT-09-001 (2009), http://cds.cern.ch/record/1194487
    [21] F. Beaudette, The CMS Particle Flow Algorithm, pp. 295-304. Proceedings, International Conference on Calorimetry for the High Energy Frontier (CHEF 2013): Paris, France, April 22-25, 2013, arXiv:1401.8155[hep-ex]
    [22] P. Mora de Freitas and H. Videau, Detector simulation with MOKKA/GEANT4: Present and future, Presented at the International Workshop on physics and experiments with future electron-positron linear colliders, Jeju Island, Korea (2002), http://inspirehep.net/record/609687
    [23] S. Agostinelli et al, Nucl. Instrum. Meth. A, 506: 250 (2003)
    [24] M. Ruan et al, Eur. Phys. J. C, 78: 426 (2018)
    [25] H. Zhao, Y.-F. Zhu, C.-D. Fu, D. Yu, and M.-Q. Ruan, The Higgs Signatures at the CEPC CDR Baseline, arXiv:1806.04992[hep-ex]
    [26] D. Yu, M. Ruan, V. Boudry, and H. Videau, Eur. Phys. J. C, 77: 591 (2017), arXiv:1701.07542[physics.ins-det
    [27] S. Catani, Y. L. Dokshitzer, M. Olsson, G. Turnock, and B. Webber, Phys. Lett. B, 269: 432 (1991)
    [28] ATLAS Collaboration, Eur. Phys. J. C, 73: 2306 (2013), arXiv:1210.6210[hep-ex
    [29] CMS Collaboration, JINST, 12: P02014 (2017), arXiv:1607.03663[hep-ex doi: 10.1088/1748-0221/12/02/P02014
    [30] T. Tanabe and T. Suehara, LCFIPlus, Presented at ILD workshop at Kyushu University (2012). https://agenda.linearcollider.org/event/5496/session/1/contribution/16.
    [31] LHC Higgs Cross Section Working Group, Handbook of LHC Higgs Cross Sections: 1. Inclusive Observables, arXiv: 1101.0593[hep-ph]
    [32] LHC Higgs Cross Section Working Group, Handbook of LHC Higgs Cross Sections: 2. Differential Distributions, arXiv: 1201.3084[hep-ph]
    [33] LHC Higgs Cross Section Working Group, Handbook of LHC Higgs Cross Sections: 3. Higgs Properties, arXiv: 1307.1347[hep-ph]
    [34] W. Kilian, T. Ohl, and J. Reuter, Eur. Phys. J. C, 71: 1742 (2011), arXiv:0708.4233[hep-ph
    [35] C. M. Carloni Calame, G. Montagna, O. Nicrosini, and F. Piccinini, Nucl. Phys. Proc. Suppl., 131: 48-55 (2004), arXiv:hep-ph/0312014 doi: 10.1016/j.nuclphysbps.2004.02.008
    [36] S. Aplin, J. Engels, and F. Gaede, A production system for massive data processing in ILCSoft,. http://inspirehep.net/record/889841
    [37] D. Asner, T. Barklow, C. Calancha, K. Fujii, N. Graf, et al., ILC Higgs White Paper, arXiv:1310.0763[hep-ph].
    [38] Y. Haddad, Feasibility of a minimum bias analysis of e+e → ZH → qq+X at a 250 GeV ILC, 010201-44 Chinese Physics C, to be publishing, arXiv:1404.3164[hep-ph]
    [39] M. Oreglia, A Study of the Reactions ψ’→ γγψ, SLAC-R-0236(1980), http://www.slac.stanford.edu/cgi-wrap/getdoc/slac-r-23\6.pdf
    [40] ATLAS Collaboration, Submitted to: Phys. Lett. B, (2018), arXiv:1808.08238[hep-ex
    [41] CMS Collaboration, Submitted to: Phys. Rev. Lett., (2018), arXiv:1808.08242[hep-ex
    [42] R. E. Shrock and M. Suzuki, Phys. Lett. B, 110: 250 (1982)
    [43] K. Griest and H. E. Haber, Phys. Rev. D, 37: 719 (1988)
    [44] C. Englert, T. Plehn, D. Zerwas, and P. M. Zerwas, Phys. Lett. B, 703: 298-305 (2011), arXiv:1106.3097[hep-ph
    [45] C. Bonilla, J. W. F. Valle, and J. C. Romo, Phys. Rev. D, 91: 113015 (2015), arXiv:1502.01649[hep-ph
    [46] ALEPH, Phys. Rept., 532: 119 (2013), arXiv:1302.3415[hep-ex doi: 10.1016/j.physrep.2013.07.004
    [47] S. Dawson, A. Gritsan, H. Logan, J. Qian, C. Tully et al, Working Group Report: Higgs Boson, arXiv:1310.8361[hep-ex]
    [48] LHC Higgs Cross Section Working Group, Handbook of LHC Higgs Cross Sections: 4. Deciphering the Nature of the Higgs Sector, arXiv:1610.07922[hep-ph]
    [49] T. Han, Z. Liu, and J. Sayre, Phys. Rev. D, 89: 113006 (2014), arXiv:1311.7155[hep-ph
    [50] M. E. Peskin, Estimation of LHC and ILC Capabilities for Precision Higgs Boson Coupling Measurements, arXiv:1312.4974[hep-ph]
    [51] P. P. Giardino, K. Kannike, I. Masina, M. Raidal, and A. Strumia, JHEP, 05: 046 (2014), arXiv:1303.3570[hep-ph
    [52] G. Belanger, B. Dumont, U. Ellwanger, J. F. Gunion, and S. Kraml, Phys. Rev. D, 88: 075008 (2013), arXiv:1306.2941[hep-ph
    [53] P. Bechtle, S. Heinemeyer, O. Stal, T. Stefaniak, and G. Weiglein, JHEP, 11: 039 (2014), arXiv:1403.1582[hep-ph
    [54] K. Cheung, J. S. Lee, and P.-Y. Tseng, Phys. Rev. D, 90: 095009 (2014), arXiv:1407.8236[hep-ph
    [55] S. Fichet and G. Moreau, Nucl. Phys. B, 905: 391-446 (2016), arXiv:1509.00472[hep-ph
    [56] R. Lafaye, T. Plehn, M. Rauch, and D. Zerwas, Phys. Rev. D, 96(no. 7): 075044 (2017), arXiv:1706.02174[hep-ph
    [57] V. Barger, T. Han, P. Langacker, B. McElrath, and P. Zerwas, Phys. Rev. D, 67: 115001 (2003), arXiv:hep-ph/0301097
    [58] T. Corbett, O. J. P. Eboli, J. Gonzalez-Fraile, and M. C. Gonzalez-Garcia, Phys. Rev. D, 87: 015022 (2013), arXiv:1211.4580[hep-ph
    [59] J. Elias-Miro, J. R. Espinosa, E. Masso, and A. Pomarol, JHEP, 11: 066 (2013), arXiv:1308.1879[hep-ph
    [60] A. Pomarol and F. Riva, JHEP, 01: 151 (2014), arXiv:1308.2803[hep-ph
    [61] G. Amar, S. Banerjee, S. von Buddenbrock, A. S. Cornell, T. Mandal, B. Mellado, and B. Mukhopadhyaya, JHEP, 02: 128 (2015), arXiv:1405.3957[hep-ph
    [62] J. Ellis, V. Sanz, and T. You, JHEP, 03: 157 (2015), arXiv:1410.7703[hep-ph
    [63] A. Falkowski, Pramana, 87: 39 (2016), arXiv:1505.00046[hep-ph doi: 10.1007/s12043-016-1251-5
    [64] A. Butter, O. J. P. Éboli, J. Gonzalez-Fraile, M. C. Gonzalez-Garcia, T. Plehn, and M. Rauch, JHEP, 07: 152 (2016), arXiv:1604.03105[hep-ph
    [65] N. Craig, M. Farina, M. McCullough, and M. Perelstein, JHEP, 03: 146 (2015), arXiv:1411.0676[hep-ph
    [66] M. Beneke, D. Boito, and Y.-M. Wang, JHEP, 11: 028 (2014), arXiv:1406.1361[hep-ph
    [67] N. Craig, J. Gu, Z. Liu, and K. Wang, JHEP, 03: 050 (2016), arXiv:1512.06877[hep-ph
    [68] J. Ellis and T. You, JHEP, 03: 089 (2016), arXiv:1510.04561[hep-ph
    [69] S.-F. Ge, H.-J. He, and R.-Q. Xiao, JHEP, 10: 007 (2016), arXiv:1603.03385[hep-ph
    [70] J. Ellis, P. Roloff, V. Sanz, and T. You, Dimension-6 Operator Analysis of the CLIC Sensitivity to New Physics, arXiv:1701.04804[hep-ph]
    [71] G. Durieux, C. Grojean, J. Gu, and K. Wang, JHEP, 09: 014 (2017), arXiv:1704.02333[hep-ph
    [72] T. Barklow, K. Fujii, S. Jung, R. Karl, J. List, T. Ogawa, M. E. Peskin, and J. Tian, Phys. Rev. D, 97(no. 5): 053003 (2018), arXiv:1708.08912[hep-ph
    [73] T. Barklow, K. Fujii, S. Jung, M. E. Peskin, and J. Tian, Phys. Rev. D, 97(no. 5): 053004 (2018), arXiv:1708.09079[hep-ph
    [74] S. Di Vita, G. Durieux, C. Grojean, J. Gu, Z. Liu, G. Panico, M. Riembau, and T. Vantalon, JHEP, 02: 178 (2018), arXiv:1711.03978[hep-ph
    [75] W. H. Chiu, S. C. Leung, T. Liu, K.-F. Lyu, and L.-T. Wang, JHEP, 05: 081 (2018), arXiv:1711.04046[hep-ph
    [76] J. Ellis, C. W. Murphy, V. Sanz, and T. You, JHEP, 06: 146 (2018), arXiv:1803.03252[hep-ph
    [77] G. Durieux, J. Gu, E. Vryonidou, and C. Zhang, Chin. Phys. C, 42(no. 12): 123107 (2018), arXiv:1809.03520[hep-ph
    [78] J. Gao, JHEP, 01: 038 (2018), arXiv:1608.01746[hep-ph
    [79] LHC Higgs Cross Section Working Group, LHC HXSWG interim recommendations to explore the coupling structure of a Higgs-like particle, Chinese Physics C, to be publishing, arXiv:1209.0040[hep-ph]
    [80] CMS Collaboration, Phys. Rev. Lett., 120: 231801 (2018), arXiv:1804.02610[hep-ex doi: 10.1103/PhysRevLett.120.231801
    [81] ATLAS Collaboration, Phys. Lett. B, 784: 173-191 (2018), arXiv:1806.00425[hep-ex
    [82] A. Banfi, A. Martin, and V. Sanz, JHEP, 1408: 053 (2014), arXiv:1308.4771[hep-ph
    [83] A. Azatov and A. Paul, JHEP, 1401: 014 (2014), arXiv:1309.5273[hep-ph
    [84] C. Grojean, E. Salvioni, M. Schlaffer, and A. Weiler, JHEP, 1405: 022 (2014), arXiv:1312.3317[hep-ph
    [85] M. Buschmann, C. Englert, D. Goncalves, T. Plehn, and M. Spannowsky, Phys. Rev. D, 90: 013010 (2014), arXiv:1405.7651[hep-ph
    [86] J. Ellis, V. Sanz, and T. You, JHEP, 1407: 036 (2014), arXiv:1404.3667[hep-ph
    [87] M. Klute, R. Lafaye, T. Plehn, M. Rauch, and D. Zerwas, Europhys. Lett., 101: 51001 (2013), arXiv:1301.1322[hep-ph doi: 10.1209/0295-5075/101/51001
    [88] K. Fujii et al, Physics Case for the 250 GeV Stage of the International Linear Collider, arXiv:1710.07621[hep-ex]
    [89] Q.-F. Sun, F. Feng, Y. Jia, and W.-L. Sang, Phys. Rev. D, 96(no. 5): 051301 (2017), arXiv:1609.03995[hep-ph
    [90] Y. Gong, Z. Li, X. Xu, L. L. Yang, and X. Zhao, Phys. Rev. D, 95(no. 9): 093003 (2017), arXiv:1609.03955[hep-ph
    [91] G. P. Lepage, P. B. Mackenzie, and M. E. Peskin, Expected Precision of Higgs Boson Partial Widths within the Standard Model, arXiv:1404.0319[hep-ph]
    [92] A. Denner, S. Heinemeyer, I. Puljak, D. Rebuzzi, and M. Spira, Eur. Phys. J. C, 71: 1753 (2011), arXiv:1107.5909[hep-ph
    [93] L. G. Almeida, S. J. Lee, S. Pokorski, and J. D. Wells, Phys. Rev. D, 89: 033006 (2014), arXiv:1311.6721[hep-ph
    [94] C. Bernaciak, T. Plehn, P. Schichtel, and J. Tattersall, Phys. Rev. D, 91: 035024 (2015), arXiv:1411.7699[hep-ph
    [95] S. M. Barr and A. Zee, Phys. Rev. Lett., 65: 21-24 (1990) doi: 10.1103/PhysRevLett.65.21
    [96] J. Fan and M. Reece, JHEP, 06: 004 (2013), arXiv:1301.2597[hep-ph
    [97] ACME Collaboration, Science, 343: 269-272 (2014), arXiv:1310.7534[physics.atom-ph doi: 10.1126/science.1248213
    [98] Y. T. Chien, V. Cirigliano, W. Dekens, J. de Vries, and E. Mereghetti, JHEP, 02: 011 (2016), arXiv:1510.00725[hep-ph
    [99] R. Harnik, A. Martin, T. Okui, R. Primulando, and F. Yu, Phys. Rev. D, 88: 076009 (2013), arXiv:1308.1094[hep-ph
    [100] K. Hagiwara, S. Ishihara, R. Szalapski, and D. Zeppenfeld, Phys. Rev. D, 48: 2182-2203 (1993)
    [101] G. Gounaris et al, Triple gauge boson couplings, arXiv:hep-ph/9601233[hep-ph]
    [102] L. Bian, J. Shu, and Y. Zhang, JHEP, 09: 206 (2015), arXiv:1507.02238[hep-ph
    [103] A. Falkowski, Higgs Basis: Proposal for an EFT basis choice for LHC HXSWG, LHCHXSWG-INT-2015-001(March, 2015), https://cds.cern.ch/record/2001958
    [104] ATLAS Collaboration, Projections for measurements of Higgs boson cross sections, branching ratios and coupling parameters with the ATLAS detector at a HL-LHC, ATL-PHYS-PUB-2013-014(2013). https://cds.cern.ch/record/1611186
    [105] ATLAS Collaboration, HL-LHC projections for signal and background yield measurements of the H → γγ when the Higgs boson is produced in association with t quarks, W or Z bosons, ATL-PHYS-PUB-2014-012(2014). https://cds.cern.ch/record/1741011
    [106] ATLAS Collaboration, Update of the prospects for the H → Zγ search at the High-Luminosity LHC, ATL-PHYS-PUB-2014-006(2014). https://cds.cern.ch/record/1703276
    [107] ATLAS Collaboration, Prospects for the study of the Higgs boson in the VH(bb) channel at HL-LHC, ATL-PHYS-PUB-2014-011(2014). https://cds.cern.ch/record/1740962
    [108] ATLAS Collaboration, Studies of the VBF H → τlτhad analysis at High Luminosity LHC conditions, ATL-PHYS-PUB-2014-018(2014). https://cds.cern.ch/record/1956732
    [109] A. Falkowski, M. Gonzalez-Alonso, A. Greljo, and D. Marzocca, Phys. Rev. Lett., 116: 011801 (2016), arXiv:1508.00581[hep-ph doi: 10.1103/PhysRevLett.116.011801
    [110] R. Contino, A. Falkowski, F. Goertz, C. Grojean, and F. Riva, JHEP, 07: 144 (2016), arXiv:1604.06444[hep-ph
    [111] A. Falkowski, M. Gonzalez-Alonso, A. Greljo, D. Marzocca, and M. Son, JHEP, 02: 115 (2017), arXiv:1609.06312[hep-ph
    [112] Z. Zhang, Phys. Rev. Lett., 118: 011803 (2017), arXiv:1610.01618[hep-ph doi: 10.1103/PhysRevLett.118.011803
    [113] ATLAS Collaboration, Search for the Standard Model Higgs and Z Boson decays to J/ψ γ: HL-LHC projections, ATL-PHYS-PUB-2015-043(2015). http://cds.cern.ch/record/2054550
    [114] G. T. Bodwin, F. Petriello, S. Stoynev, and M. Velasco, Phys. Rev. D, 88: 053003 (2013), arXiv:1306.5770[hep-ph
    [115] G. Perez, Y. Soreq, E. Stamou, and K. Tobioka, Phys. Rev. D, 92: 033016 (2015), arXiv:1503.00290[hep-ph
    [116] I. Brivio, F. Goertz, and G. Isidori, Phys. Rev. Lett., 115: 211801 (2015), arXiv:1507.02916[hep-ph doi: 10.1103/PhysRevLett.115.211801
    [117] F. Bishara, U. Haisch, P. F. Monni, and E. Re, Constraining Light-Quark Yukawa Couplings from Higgs Distributions, arXiv:1606.09253[hep-ph]
    [118] L. M. Carpenter, T. Han, K. Hendricks, Z. Qian, and N. Zhou, Phys. Rev. D, 95: 053003 (2017), arXiv:1611.05463[hep-ph
    [119] ATLAS Collaboration, Study of the double Higgs production channel H(→ bb)H(→ γγ) with the ATLAS experiment at the HL-LHC, ATL-PHYS-PUB-2017-001, CERN, Geneva, 010201-46 Chinese Physics C, 010201 Jan, 2017. https://cds.cern.ch/record/2243387
    [120] R. Contino et al, CERN Yellow Report,255-440 (2017), arXiv:1606.09408[hep-ph
    [121] M. McCullough, Phys. Rev. D, 90: 015001 (2014), arXiv:1312.3322[hep-ph
    [122] S. Di Vita, C. Grojean, G. Panico, M. Riembau, and T. Vantalon, JHEP, 09: 069 (2017), arXiv:1704.01953[hep-ph
    [123] A. Azatov, R. Contino, G. Panico, and M. Son, Phys. Rev. D, 92: 035001 (2015), arXiv:1502.00539[hep-ph
    [124] Z. Liu, I. Low, and L.-T. Wang, Higgs-Top Interactions at Future Circular e+e Colliders, arXiv:2018.nnnn
    [125] G. Li, H.-R. Wang, and S.-h. Zhu, Probing CP-violating htt coupling in e+e → hγ, arXiv:1506.06453[hep-ph]
    [126] E. Vryonidou and C. Zhang, Dimension-six electroweak top-loop effects in Higgs production and decay, arXiv:1804.09766[hep-ph]
    [127] G. Durieux, M. Perello, M. Vos, and C. Zhang, JHEP, 10: 168 (2018), arXiv:1807.02121[hep-ph
    [128] S. Boselli, R. Hunter, and A. Mitov, Prospects for the determination of the top-quark Yukawa coupling at future e+e colliders, arXiv:1805.12027[hep-ph]
    [129] J. A. Aguilar-Saavedra, Nucl. Phys. B, 812: 181-204 (2009), arXiv:0811.3842[hep-ph
    [130] J. A. Aguilar-Saavedra, Nucl. Phys. B, 821: 215-227 (2009), arXiv:0904.2387[hep-ph
    [131] Y. Gao, A. V. Gritsan, Z. Guo, K. Melnikov, M. Schulze, and N. V. Tran, Phys. Rev. D, 81: 075022 (2010), arXiv:1001.3396[hep-ph
    [132] S. Bolognesi, Y. Gao, A. V. Gritsan, K. Melnikov, M. Schulze, N. V. Tran, and A. Whitbeck, Phys. Rev. D, 86: 095031 (2012), arXiv:1208.4018[hep-ph
    [133] I. Anderson, et al, Phys. Rev. D, 89: 035007 (2014), arXiv:1309.4819[hep-ph
    [134] J. Alwall, P. Demin, S. de Visscher, R. Frederix, M. Herquet, F. Maltoni, T. Plehn, D. L. Rainwater, and T. Stelzer, JHEP, 09: 028 (2007), arXiv:0706.2334[hep-ph
    [135] D. Curtin, R. Essig, S. Gori, P. Jaiswal, A. Katz, et al, Phys.Rev. D, 90: 075004 (2014), arXiv:1312.4992[hep-ph
    [136] Z. Liu, L.-T. Wang, and H. Zhang, Exotic decays of the 125 GeV Higgs boson at future e+e lepton colliders, arXiv:1612.09284[hep-ph]
    [137] J. Liu, Z. Liu, and L.-T. Wang, Long-lived particles at the LHC: catching them in time, arXiv:1805.05957[hep-ph]
    [138] A. Thamm, R. Torre, and A. Wulzer, JHEP, 07: 100 (2015), arXiv:1502.01701[hep-ph
    [139] J. Gu, H. Li, Z. Liu, S. Su, and W. Su, JHEP, 12: 153 (2017), arXiv:1709.06103[hep-ph
    [140] J. Fan, M. Reece, and L.-T. Wang, Precision Natural SUSY at CEPC, FCC-ee, and ILC, arXiv:1412.3107[hep-ph]
    [141] R. Essig, P. Meade, H. Ramani, and Y.-M. Zhong, JHEP, 09: 085 (2017), arXiv:1707.03399[hep-ph
    [142] Z. Chacko, H.-S. Goh, and R. Harnik, Phys. Rev. Lett., 96: 231802 (2006), arXiv:hep-ph/0506256 doi: 10.1103/PhysRevLett.96.231802
    [143] G. Burdman, Z. Chacko, H.-S. Goh, and R. Harnik, JHEP, 02: 009 (2007), arXiv:hep-ph/0609152
    [144] N. Craig, C. Englert, and M. McCullough, Phys. Rev. Lett., 111: 121803 (2013), arXiv:1305.5251[hep-ph doi: 10.1103/PhysRevLett.111.121803
    [145] K. Kajantie, M. Laine, K. Rummukainen, and M. E. Shaposhnikov, Phys. Rev. Lett., 77: 2887-2890 (1996), arXiv:hep-ph/9605288 doi: 10.1103/PhysRevLett.77.2887
    [146] S. Profumo, M. J. Ramsey-Musolf, C. L. Wainwright, and P. Winslow, Phys. Rev. D, 91: 035018 (2015), arXiv:1407.5342[hep-ph
    [147] A. Katz and M. Perelstein, JHEP, 07: 108 (2014), arXiv:1401.1827[hep-ph
    [148] A. Noble and M. Perelstein, Phys. Rev. D, 78: 063518 (2008), arXiv:0711.3018[hep-ph
    [149] B. Henning, X. Lu, and H. Murayama, What do precision Higgs measurements buy us?, arXiv:1404.1058[hep-ph]
    [150] S. Profumo, M. J. Ramsey-Musolf, and G. Shaughnessy, JHEP, 08: 010 (2007), arXiv:0705.2425[hep-ph
    [151] D. Curtin, P. Meade, and C.-T. Yu, JHEP, 11: 127 (2014), arXiv:1409.0005[hep-ph
    [152] Q.-H. Cao, F. P. Huang, K.-P. Xie, and X. Zhang, Chin. Phys. C, 42: 023103 (2018), arXiv:1708.04737[hep-ph
    [153] P. Huang, A. J. Long, and L.-T. Wang, Phys. Rev. D, 94: 075008 (2016), arXiv:1608.06619[hep-ph
  • 加载中

Figures(37) / Tables(17)

Get Citation
Fenfen An, Yu Bai, Chunhui Chen, Xin Chen, Zhenxing Chen, Joao Guimaraes da Costa, Zhenwei Cui, Yaquan Fang, Chengdong Fu, Jun Gao, Yanyan Gao, Yuanning Gao, Shao-Feng Ge, Jiayin Gu, Fangyi Guo, Jun Guo, Tao Han, Shuang Han, Hong-Jian He, Xianke He, Xiao-Gang He, Jifeng Hu, Shih-Chieh Hsu, Shan Jin, Maoqiang Jing, Susmita Jyotishmati, Chia-Ming Kuo, Pei-Zhu Lai, Boyang Li, Congqiao Li, Gang Li, Haifeng Li, Liang Li, Shu Li, Tong Li, Qiang Li, Hao Liang, Zhijun Liang, Libo Liao, Bo Liu, Jianbei Liu, Tao Liu, Zhen Liu, Xinchou Lou, Lianliang Ma, Bruce Mellado, Xin Mo, Mila Pandurovic, Jianming Qian, Zhuoni Qian, Nikolaos Rompotis, Manqi Ruan, Alex Schuy, Lian-You Shan, Jingyuan Shi, Xin Shi, Shufang Su, Dayong Wang, Jin Wang, Lian-Tao Wang, Yifang Wang, Yuqian Wei, Yue Xu, Haijun Yang, Ying Yang, Weiming Yao, Dan Yu, Kaili Zhang, Zhaoru Zhang, Mingrui Zhao, Xianghu Zhao and Ning Zhou. Precision Higgs physics at the CEPC[J]. Chinese Physics C. doi: 10.1088/1674-1137/43/4/043002
Fenfen An, Yu Bai, Chunhui Chen, Xin Chen, Zhenxing Chen, Joao Guimaraes da Costa, Zhenwei Cui, Yaquan Fang, Chengdong Fu, Jun Gao, Yanyan Gao, Yuanning Gao, Shao-Feng Ge, Jiayin Gu, Fangyi Guo, Jun Guo, Tao Han, Shuang Han, Hong-Jian He, Xianke He, Xiao-Gang He, Jifeng Hu, Shih-Chieh Hsu, Shan Jin, Maoqiang Jing, Susmita Jyotishmati, Chia-Ming Kuo, Pei-Zhu Lai, Boyang Li, Congqiao Li, Gang Li, Haifeng Li, Liang Li, Shu Li, Tong Li, Qiang Li, Hao Liang, Zhijun Liang, Libo Liao, Bo Liu, Jianbei Liu, Tao Liu, Zhen Liu, Xinchou Lou, Lianliang Ma, Bruce Mellado, Xin Mo, Mila Pandurovic, Jianming Qian, Zhuoni Qian, Nikolaos Rompotis, Manqi Ruan, Alex Schuy, Lian-You Shan, Jingyuan Shi, Xin Shi, Shufang Su, Dayong Wang, Jin Wang, Lian-Tao Wang, Yifang Wang, Yuqian Wei, Yue Xu, Haijun Yang, Ying Yang, Weiming Yao, Dan Yu, Kaili Zhang, Zhaoru Zhang, Mingrui Zhao, Xianghu Zhao and Ning Zhou. Precision Higgs physics at the CEPC[J]. Chinese Physics C.  doi: 10.1088/1674-1137/43/4/043002 shu
Milestone
Received: 2018-11-09
Revised: 2019-01-21
Article Metric

Article Views(2765)
PDF Downloads(91)
Cited by(0)
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:

Precision Higgs physics at the CEPC

    Corresponding author: Yaquan Fang, fangyq@ihep.ac.cn
  • 1. Department of Modern Physics, University of Science and Technology of China, Anhui 230026, China
  • 2. China Institute of Atomic Energy, Beijing 102413, China
  • 3. School of Physics, Peking University, Beijing 100871, China
  • 4. Institute of High Energy Physics, Beijing 100049, China
  • 5. Department of Engineering Physics, Physics Department, Tsinghua University, Beijing 100084, China
  • 6. University of Chinese Academy of Science (UCAS), Beijing 100049, China
  • 7. School of Nuclear Science and Technology, University of South China, Hengyang 421001, China
  • 8. Department of Physics, Nanjing University, Nanjing 210093, China
  • 9. Department of Physics, Southeast University, Nanjing 210096, China
  • 10. School of Physics and Astronomy, Shanghai Jiao Tong University, KLPPAC-MoE, SKLPPC, Shanghai 200240, China
  • 11. Tsung-Dao Lee Institute, Shanghai 200240, China
  • 12. Institute of Frontier and Interdisciplinary Science and Key Laboratory of Particle Physics and Particle Irradiation (MOE), Shandong University, Qingdao 266237, China
  • 13. PRISMA Cluster of Excellence & Mainz Institute of Theoretical Physics, Johannes Gutenberg-Universität Mainz, Mainz 55128, Germany
  • 14. Department of Physics, Hong Kong University of Science and Technology, Hong Kong
  • 15. Kavli IPMU (WPI), UTIAS, The University of Tokyo, Kashiwa, Chiba 277-8583, Japan
  • 16. Vinca Institute of Nuclear Sciences, University of Belgrade, Belgrade 11000, Serbia
  • 17. School of Physics and Institute for Collider Particle Physics, University of the Witwatersrand, Johannesburg 2050, South Africa
  • 18. iThemba LABS, National Research Foundation, PO Box 722, Somerset West 7129, South Africa
  • 19. Center for Theoretical Physics of the Universe, Institute of Basic Science, Daejeon 34126, South Korea
  • 20. Department of Physics, National Taiwan University, Taipei 10617
  • 21. Department of Physics and Center for High Energy and High Field Physics, National Central University, Taoyuan City 32001
  • 22. Department of Physics, University of Liverpool, Liverpool L69 7ZX, United Kingdom
  • 23. Department of Physics and Astronomy, Iowa State University, Ames 50011-3160, USA
  • 24. Department of Physics, University of Michigan, Ann Arbor, Michigan 48109, USA
  • 25. Department of Physics, University of Arizona, Arizona 85721, USA
  • 26. Theoretical Physics Department, Fermi National Accelerator Laboratory, Batavia 60510, USA
  • 27. Department of Physics, University of Chicago, Chicago 60637, USA
  • 28. Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
  • 29. Department of Physics, University of California, Berkeley, California 94720, USA
  • 30. Maryland Center for Fundamental Physics, Department of Physics, University of Maryland, College Park, Maryland 20742, USA
  • 31. Department of Physics & Astronomy, University of Pittsburgh, Pittsburgh 15260, USA
  • 32. Department of Physics, University of Washington, Seattle 98195-1560, USA
  • 33. Department of Physics, University of Texas at Dallas, Texas 75080-3021, USA
  • 34. Physical Science Laboratory, Huairou National Comprehensive Science Center, Beijing, 101400, China

Abstract: The discovery of the Higgs boson with its mass around 125 GeV by the ATLAS and CMS Collaborations marked the beginning of a new era in high energy physics. The Higgs boson will be the subject of extensive studies of the ongoing LHC program. At the same time, lepton collider based Higgs factories have been proposed as a possible next step beyond the LHC, with its main goal to precisely measure the properties of the Higgs boson and probe potential new physics associated with the Higgs boson. The Circular Electron Positron Collider (CEPC) is one of such proposed Higgs factories. The CEPC is an e+e circular collider proposed by and to be hosted in China. Located in a tunnel of approximately 100 km in circumference, it will operate at a center-of-mass energy of 240 GeV as the Higgs factory. In this paper, we present the first estimates on the precision of the Higgs boson property measurements achievable at the CEPC and discuss implications of these measurements.

    HTML

    1.   Introduction
    • The historic discovery of a Higgs boson in 2012 by the ATLAS and CMS collaborations [1,2] at the Large Hadron Collider (LHC) has opened a new era in particle physics. Subsequent measurements of the properties of the new particle have indicated compatibility with the Standard Model (SM) Higgs boson [3-9]. While the SM has been remarkably successful in describing experimental phenomena, it is important to recognize that it is not a complete theory. In particular, it does not predict the parameters in the Higgs potential, such as the Higgs boson mass. The vast difference between the Planck scale and the weak scale remains a major mystery. There is not a complete understanding of the nature of electroweak phase transition. The discovery of a spin zero Higgs boson, the first elementary particle of its kind, only sharpens these questions. It is clear that any attempt of addressing these questions will involve new physics beyond the SM (BSM). Therefore, the Higgs boson discovery marks the beginning of a new era of theoretical and experimental explorations.

      A physics program of the precision measurements of the Higgs boson properties will be a critical component of any road map for high energy physics in the coming decades. Potential new physics beyond the SM could lead to observable deviations in the Higgs boson couplings from the SM expectations. Typically, such deviations can be parametrized as

      δ=cv2M2NP,

      (1)

      where v and MNP are the vacuum expectation value of the Higgs field and the typical mass scale of new physics, respectively. The size of the proportionality constant c depends on the model, but it should not be much larger than O(1). The high-luminosity LHC (HL-LHC) will measure the Higgs boson couplings to about 5% [10,11]. At the same time, the LHC will directly search for new physics from a few hundreds of GeV to at least one TeV. Eq. (1) implies that probing new physics significantly beyond the LHC reach would require the measurements of the Higgs boson couplings at least at percent level accuracy. To achieve such precision will need new facilities, a lepton collider operating as a Higgs factory is a natural next step.

      The Circular Electron-Positron Collider (CEPC) [12], proposed by the Chinese particle physics community, is one of such possible facilities. The CEPC will be placed in a tunnel with a circumference of approximately 100 km and will operate at a center-of-mass energy of s240 GeV, near the maximum of the Higgs boson production cross section through the e+eZH process. At the CEPC, in contrast to the LHC, Higgs boson candidates can be identified through a technique known as the recoil mass method without tagging its decays. Therefore, the Higgs boson production can be disentangled from its decay in a model independent way. Moreover, the cleaner environment at a lepton collider allows much better exclusive measurements of Higgs boson decay channels. All of these give the CEPC an impressive reach in probing Higgs boson properties. With the expected integrated luminosity of 5.6ab1, over one million Higgs bosons will be produced. With this sample, the CEPC will be able to measure the Higgs boson coupling to the Z boson with an accuracy of 0.25%, more than a factor of 10 better than the HL-LHC [10,11]. Such a precise measurement gives the CEPC unprecedented reach into interesting new physics scenarios which are difficult to probe at the LHC. The CEPC also has strong capability in detecting Higgs boson invisible decay. It is sensitive to the invisible decay branching ratio down to 0.3%. In addition, it is expected to have good sensitivities to exotic decay channels which are swamped by backgrounds at the LHC. It is also important to stress that an e+e Higgs factory can perform model independent measurement of the Higgs boson width. This unique feature in turn allows for the model independent determination of the Higgs boson couplings.

      This paper documents the first studies of a precision Higgs boson physics program at the CEPC. It is organized as follows: Section 2 briefly summarizes the collider and detector performance parameters assumed for the studies. Section 3 gives an overview of relevant e+e collision processes and Monte Carlo simulations. Sections 4 and 5 describe inclusive and exclusive Higgs boson measurements. Section 6 discusses the combined analysis to extract Higgs boson production and decay properties. Section 7 interprets the results in the coupling and effective theory frameworks. Section 8 estimates the reaches in the test of Higgs boson spin/CP properties and in constraining the exotic decays of the Higgs boson based on previously published phenomenological studies. Finally the implications of all these measurements are discussed in Section 9.

    2.   CEPC detector concept

      2.1.   The CEPC operating scenarios

    • The CEPC is designed to operate as a Higgs factory at s=240 GeV and as a Z factory at s=91.2 GeV. It will also perform WW threshold scans around s=160 GeV. Table 1 shows potential CEPC operating scenarios and the expected numbers of H, W and Z bosons produced in these scenarios.

      operation modeZ factoryWW thresholdHiggs factory
      s/GeV91.2160240
      run time/y217
      instantaneous luminosity/(1034cm2s1)16–32103
      integrated luminosity/(ab1)8–162.65.6
      Higgs boson yield106
      W boson yield107108
      Z boson yield10111012108108

      Table 1.  CEPC operating scenarios and the numbers of Higgs, W and Z bosons produced. The integrated luminosity and the event yields assume two interaction points. The ranges of luminosities and the Z yield of the Z factory operation correspond to detector solenoid field of 3 and 2 Tesla.

      The CEPC operation as a Higgs factory will run for 7 years and produce a total of 1 million Higgs bosons with two interaction points. Meanwhile, approximately 100 million W bosons and 1 billion Z bosons will also be produced in this operation. These large samples of W and Z bosons will allow for in-situ detector characterization as well as for the precise measurements of electroweak parameters.

      Running at the WW threshold around s=160 GeV, 107 W bosons will be produced in one year. Similarly running at the Z pole around s=91.2 GeV (the Z factory), CEPC will produce 10111012 Z bosons. These large samples will enable high precision measurements of the electroweak observables such as AbFB, Rb, the Z boson line-shape parameters, the mass and width of the W boson. An order of magnitude or more improvement in the precision of these observables are foreseen.

    • 2.2.   Conceptual detector design

    • The primary physics objective of the CEPC is the precise determination of the Higgs boson properties. Therefore CEPC detectors must be able to reconstruct and identify all key physics objects that the Higgs bosons are produced with or decay into with high efficiency, purity and accuracy. These objects include charged leptons, photons, jets, missing energy and missing momentum. Moreover, the flavor tagging of jets, such as those from b, c and light quarks or gluons, are crucial for identifying the hadronic decays of the Higgs bosons. The detector requirements for the electroweak and flavor physics are similar. One notable additional requirement is the identification of charged particles such as π± and K± for the flavor physics program.

      Using the International Large Detector (ILD) [13,14] as a reference, a particle flow oriented conceptual detector design, CEPC-v1 (see Fig. 1), has been developed for the CEPC. A detailed description of the CEPC-v1 detector can be found in Ref. [15]. Originally developed for the LEP experiments [16,17], particle flow is a proven concept for event reconstruction [18-21], based on the principle of reconstructing all visible final-state particles in the most sensitive detector subsystem. Specifically, a particle-flow algorithm reconstructs charged particles in the tracking system, measures photons in the electromagnetic calorimeter and neutral hadrons in both electromagnetic and hadronic calorimeters. Physics objects are then identified or reconstructed from this list of final state particles. Particle flow reconstruction provides a coherent interpretation of an entire physics event and, therefore, is particularly well suited for the identification of composite physics objects such as the τ leptons and jets.

      Figure 1.  (color online) Conceptual CEPC detector, CEPC-v1, implemented in MOKKA [22] and GEANT 4 [23]. It is comprised of a silicon vertexing and tracking system of both pixel and strips geometry, a Time-Projection-Chamber tracker, a high granularity calorimeter system, a solenoid of 3.5 Tesla magnetic field, and a muon detector embedded in a magnetic field return yoke.

      The particle-flow algorithm requires good spatial separations of calorimeter showers induced by different final state particles for their reconstruction. It is imperative to minimize the amount of material before the calorimeter to reduce the uncertainty induced by the nuclear interactions and Bremsstrahlung radiations. Therefore, a high granularity calorimeter system and low material tracking system are implemented in the CEPC-v1 detector concept. The tracking system consists of silicon vertexing and tracking detectors as well as a Time Projection Chamber (TPC). The calorimetry system is based on the sampling technology with absorber/active-medium combination of Tungsten-Silicon for the electromagnetic calorimeter (ECAL) and Iron-Resistive Plate Chamber (RPC) for the hadronic calorimeter (HCAL). The calorimeters are segmented at about 1 channel/cm3, three orders of magnitude finer than those of the LHC detectors. Both the tracking and the calorimeter system are housed inside a solenoid of 3.5 Tesla magnetic field. The CEPC-v1 detector has a sophisticated machine-detector interface with an 1.5-meter L (the distance between the interaction point and the final focusing quadrupole magnet) to accommodate the high design luminosity. Table 2 shows the geometric parameters and the benchmark detector subsystem performance of the CEPC-v1 detector. A schematic of the detector is shown in Fig. 2.

      Figure 2.  (color online) The layout of one quarter of the CEPC-v1 detector concept.

      tracking system
       vertex detector6 pixel layers
       Silicon tracker3 barrel layers, 6 forward disks on each side
       time projection chamber220 radial readouts
      calorimetry
       ECALW/Si,24X0,5×5 mm2cell with 30 layers
       HCALFe/RPC,6λ,10×10 mm2 cell with 40 layers
      performance
       track momentum resolutionΔ(1/pT)2×105(1/GeV)
       impact parameter resolution5μm10μm/[(p/GeV)(sinθ)3/2]
       ECAL energy resolutionΔE/E16%/E/GeV1%
       HCAL energy resolutionΔE/E60%/E/GeV1%

      Table 2.  Basic parameters and performance of the CEPC-v1 detector. The radiation length (X0) and the nuclear interaction length (λ) are measured at the normal incidence. The cell sizes are for transverse readout sensors and the layer numbers are for longitudinal active readouts. The θ is the track polar angle.

    • 2.3.   Object reconstruction and identification

    • A dedicated particle flow reconstruction toolkit, ARBOR [19], has been developed for the CEPC-v1 detector. Inspired by the tree structure of particle showers, ARBOR attempts to reconstruct every visible final state particle. Figure 3 illustrates a simulated e+eZHqˉqbˉb event as reconstructed by the ARBOR algorithm. The algorithm's performance for leptons, photons and jets are briefly summarized here. More details can be found in Refs. [24,25].

      Figure 3.  (color online) A simulated e+eZHqˉq bˉb event reconstructed with the ARBOR algorithm. Different types of reconstructed final state particles are represented in different colors.

    • 2.3.1.   Leptons and Photons
    • Leptons () are fundamental for the measurements of the Higgs boson properties at the CEPC. About 7% of the Higgs bosons are produced in association with a pair of leptons through the e+eZH+H process. These events allow for the identifications of Higgs bosons using the recoil mass information and therefore enable the measurement of the ZH production cross section and the Higgs boson mass. Moreover, a significant fraction of Higgs bosons decay into final states with leptons indirectly through the leptonic decays of the W or Z bosons as well as the τ leptons. These leptons serve as signatures for identifying different Higgs boson decay modes.

      A lepton identification algorithm, LICH [26], has been developed and integrated into ARBOR. Efficiencies close to 99.9% for identifying electrons and muons with energies above 2 GeV have been achieved while the mis-identification probabilities from hadrons are limited to be less than 1%. The CEPC-v1 tracking system provides an excellent momentum resolution that is about ten times better than those of the LEP and LHC detectors. The good resolution is illustrated in the narrow invariant mass distribution of the muon pairs from the Hμ+μ decays as shown in Fig. 4(a).

      Figure 4.  Simulated invariant mass distributions of (a) muon pairs from Hμ+μ and (b) photon pairs from Hγγ, both from the e+eZH process with the Zνˉν decay. The Mμ+μ distribution is fit with a Gaussian core plus a small low-mass tail from the Bremsstrahlung radiation. The Gaussian has a width of 0.2 GeV, corresponding to a relative mass resolution of 0.16%. The Mγγ distribution is described well by a Crystal Ball function with a width of 3.1 GeV, corresponding to a relative mass resolution of 2.5%.

      Photons are essential for the studies of Hγγ and HZγ decays. They are also important for the reconstruction and measurements of τ leptons and jets. The Hγγ decay is an ideal process to characterize the photon performance of the CEPC-v1. Figure 4(b) shows the invariant mass distribution of the photon pairs from the Hγγ decays.

    • 2.3.2.   Jets
    • Approximately 70% of Higgs bosons decay directly into jets (bˉb,cˉc,gg) and an additional 22% decay indirectly into final states with jets through the HWW,ZZ cascades. Therefore, efficient jet reconstruction and precise measurements of their momenta are pre-requisite for a precision Higgs physics program. In ARBOR, jets are reconstructed using the Durham algorithm [27]. As a demonstration of the CEPC-v1 jet performance, Fig. 5 shows the reconstructed dijet invariant mass distributions of the Wqˉq, Zqˉq and Hbˉb/cˉc/gg decays from the ZZνˉνqˉq, WWνqˉq and ZHνˉν(bˉb/cˉc/gg) processes, respectively. Compared with Wqˉq, the Zqˉq and Hbˉb/cˉc/gg distributions have long low-mass tails, resulting from the heavy-flavor jets in these decays. The jet energy resolution is expected to be between 3–5% for the jet energy range relevant at the CEPC. This resolution is approximately 2–4 times better than those of the LHC experiments [28,29]. The dijet mass resolution for the W and Z bosons is approximately 4.4%, which allows for an average separation of 2σ or better of the the hadronically decaying W and Z bosons.

      Figure 5.  (color online) Distributions of the reconstructed dijet invariant mass for the Wqˉq, Zqˉq and Hbˉb/cˉc/gg decays from, respectively, the WWνqˉq, ZZνˉνqˉq and ZHνˉν(bˉb/cˉc/gg) processes. All distributions are normalized to unit area.

      Jets originating from heavy flavors (b- or c-quarks) are tagged using the LCFIPlus algorithm [30]. The algorithm combines information from the secondary vertex, jet mass, number of leptons etc to construct b-jet and c-jet discriminating variables. The tagging performance characterized using the Zqˉq decays from the Z factory operation is shown in Fig. 6. For an inclusive Zqˉq sample, b-jets can be tagged with an efficiency of 80% and a purity of 90% while the corresponding efficiency and purity for tagging c-jets are 60% and 60%, respectively.

      Figure 6.  (color online) Efficiency for tagging b-jets vs rejection for light-jet background (blue) and c-jet background (red), determined from an inclusive Zqˉq sample from the Z factory operation.

    • 2.4.   Ongoing optimization

    • The CEPC-v1 detector design is used as the reference detector for the studies summarized in this paper. A series of optimizations have been performed meanwhile, aiming to reduce power consumption and construction cost and to improve the machine-detector interface while minimizing the impact on Higgs boson physics. An updated detector concept, CEPC-v4, has thus been developed. The CEPC-v4 has a smaller solenoidal field of 3 Tesla and a reduced calorimeter dimensions along with fewer readout channels. In particular, the ECAL readout senor size is changed from 5×5 mm2 to 10×10 mm2. A new Time-of-Flight measurement capability is added to improve the flavor physics potential.

      The weaker magnetic field degrades momentum resolution for charged particles by 14%, which translates directly into a degraded muon momentum resolution. The impact on other physics objects such as electrons, photons and jets are estimated to be small as the track momentum resolution is not a dominant factor for their performance. In parallel with the detector optimization, the accelerator design has chosen 240 GeV as the nominal center-of-mass energy for the Higgs factory. However, the simulation of CEPC-v1 assumes s=250 GeV. The estimated precision of Higgs boson property measurements for CEPC-v1 operating at 250 GeV are therefore extrapolated to obtain those for CEPC-v4 at s=240 GeV, as discussed Section 6.2.

    3.   Theory and Monte Carlo samples

      3.1.   Higgs boson production and decay

    • Production processes for a 125 GeV SM Higgs boson at the CEPC operating at s240250 GeV are e+eZH (ZH associate production or Higgsstrahlung), e+eνeˉνeH (W fusion) and e+ee+eH (Z fusion) as illustrated in Fig. 7. In the following, the W and Z fusion processes are collectively referred to as vector-boson fusion (VBF) production.

      Figure 7.  Feynman diagrams of the Higgs boson production processes at the CEPC: (a) e+eZH, (b) e+eνeˉνeH and (c) e+ee+eH.

      The SM Higgs boson production cross sections as functions of center-of-mass energy are shown in Fig. 8, assuming that the mass of the Higgs boson is 125 GeV. Similarly, the Higgs boson decay branching ratios and total width are shown in Table 3. As an s-channel process, the cross section of the e+eZH process reaches its maximum at s250 GeV, and then decreases asymptotically as 1/s. The VBF process proceeds through t-channel exchange of vector bosons and its cross section increases logarithmically as ln2(s/M2V). Because of the small neutral-current Zee coupling, the VBF cross section is dominated by the W fusion process.

      Figure 8.  (color online) Production cross sections of e+eZH and e+e(e+e/νˉν)H as functions of s for a 125 GeV SM Higgs boson. The vertical indicates s=250 GeV, the energy assumed for most of the studies summarized in this paper.

      decay modebranching ratiorelative uncertainty
      Hbˉb57.7%+3.2%, −3.3%
      Hcˉc2.91%+12%, −12%
      Hτ+τ6.32%+5.7%, −5.7%
      Hμ+μ2.19×10−4+6.0%, −5.9%
      HWW21.5%+4.3%, −4.2%
      HZZ2.64%+4.3%, −4.2%
      Hγγ2.28×10−3+5.0%, −4.9%
      HZγ1.53×10−3+9.0%, −8.8%
      Hgg8.57%+10%, −10%
      ΓH4.07 MeV+4.0%, −4.0%

      Table 3.  Standard model predictions of the decay branching ratios and total width of a 125 GeV Higgs boson [31-33]. The quoted uncertainties include contributions from both theoretical and parametric sources.

      Numerical values of these cross sections at s=250 GeV are listed in Table 4. Because of the interference effects between e+eZH and e+eνeˉνeH for the Zνeˉνe decay and between e+eZH and e+ee+eH for the Ze+e decay, the cross sections of these processes cannot be separated. The breakdowns in Fig. 8 and Table 4 are for illustration only. The e+eZH cross section shown is from Fig. 7(a) only whereas the e+eνeˉνeH and e+ee+eH cross sections include contributions from their interferences with the e+eZH process.

      processcross sectionevents in 5.6 ab−1
      Higgs boson production, cross section in fb
      e+eZH204.71.15×106
      e+eνeˉνeH6.853.84×104
      e+ee+eH0.633.53×103
      total212.11.19×106
      background processes, cross section in pb
      e+ee+e(γ)(Bhabha)8504.5×109
      e+eqˉq(γ)50.22.8×108
      e+eμ+μ(γ)[orτ+τ(γ)]4.402.5×107
      e+eWW15.48.6×107
      e+eZZ1.035.8×106
      e+ee+eZ4.732.7×107
      e+ee+νW/eˉνW+5.142.9×107

      Table 4.  Cross sections of the Higgs boson production and other SM processes at s=250 GeV and numbers of events expected in 5.6 ab-1. Note that there are interferences between the same final states from different processes after the W or Z boson decays, see text. With the exception of the Bhabha process, the cross sections are calculated using the Whizard program [34]. The Bhabha cross section is calculated using the BABAYAGA event generator [35] requiring final-state particles to have |cosθ|<0.99. Photons, if any, are required to have Eγ>0.1 GeV and |cosθe±γ|<0.99.

      The CEPC as a Higgs boson factory is designed to deliver a combined integrated luminosity of 5.6ab1 to two detectors in 7 years. Over 106 Higgs boson events will be produced during this period. The large statistics, well-defined event kinematics and clean collision environment will enable the CEPC to measure the Higgs boson production cross sections as well as its properties (mass, decay width and branching ratios, etc.) with precision far beyond those achievable at the LHC. In contrast to hadron collisions, e+e collisions are unaffected by underlying events and pile-up effects. Theoretical calculations are less dependent on higher order QCD radiative corrections. Therefore, more precise tests of theoretical predictions can be performed at the CEPC. The tagging of e+eZH events using the recoil mass method (see Section 4), independent of the Higgs boson decay, is unique to lepton colliders. It provides a powerful tool to perform model-independent measurements of the inclusive e+eZH production cross section, σ(ZH), and of the Higgs boson decay branching ratios. Combinations of these measurements will allow for the determination of the total Higgs boson decay width and the extraction of the Higgs boson couplings to fermions and vector bosons. These measurements will provide sensitive probes to potential new physics beyond the SM.

    • 3.2.   Background processes

    • Apart from the Higgs boson production, other SM processes include e+ee+e (Bhabha scattering), e+eZγ (initial-state radiation return), e+eWW/ZZ (diboson) as well as the single boson production of e+ee+eZ and e+ee+νW/eˉνW+. Their cross sections and expected numbers of events for an integrated luminosity of 5.6ab1 at s=250 GeV are shown in Table 4 as well. The energy dependence of the cross sections for these and the Higgs boson production processes are shown in Fig. 9. Note that many of these processes can lead to identical final states and thus can interfere. For example, e+ee+νeWe+νeeˉνe and e+ee+eZe+eνeˉνe have the same final state after the decays of the W or Z bosons. Unless otherwise noted, these processes are simulated together to take into account interference effects for the studies presented in this paper. The breakdowns shown in the table and figure assume stable W and Z bosons, and thus are, therefore, for illustration only.

      Figure 9.  (color online) Cross sections of main SM processes of e+e collisions as functions of center-of-mass energy s obtained from the Whizard program [34]. The calculations include initial-state radiation (ISR). The W and Z fusion processes refer to e+eνˉνH and e+ee+eH production, respectively. Their numerical values at s=250 GeV can be found in Table 4.

      Along with 1.2×106 Higgs boson events, 5.8×106 ZZ, 8.6×107 WW and 2.8×108 qˉq(γ) events will be produced. Though these events are backgrounds to Higgs boson events, they are important for the calibration and characterization of the detector performance and for the measurements of electroweak parameters.

    • 3.3.   Event generation and simulation

    • The following software tools have been used to generate events, simulate detector responses and reconstruct simulated events. A full set of SM samples, including both the Higgs boson signal and SM background events, are generated with WHIZARD [34]. The generated events are then processed with MokkaC [22], the official CEPC simulation software based on the framework used for ILC studies [36]. Limited by computing resources, background samples are often pre-selected with loose generator-level requirements or processed with fast simulation tools.

      All Higgs boson signal samples and part of the leading background samples are processed with Geant4 [23] based full detector simulation and reconstruction. The rest of backgrounds are simulated with a dedicated fast simulation tool, where the detector acceptance, efficiency, intrinsic resolution for different physics objects are parametrized. Samples simulated for ILC studies [37] are used for cross checks of some studies.

    4.   Higgs boson tagging using recoil mass
    • Perhaps the most striking difference between hadron-hadron and e+e collisions is that electrons are fundamental particles whereas hadrons are composite. Consequently the energy of e+e collisions is known. Therefore through conservation laws, the energy and momentum of a Higgs boson can be inferred from other particles in an event without examining the Higgs boson itself. For a Higgsstrahlung event where the Z boson decays to a pair of visible fermions (ff), the mass of the system recoiling against the Z boson, commonly known as the recoil mass, can be calculated assuming the event has a total energy s and zero total momentum:

      M2recoil=(sEff)2p2ff=s2Effs+m2ff.

      (2)

      Here, Eff, pff and mff are, respectively, the total energy, momentum and invariant mass of the fermion pair. The Mrecoil distribution should show a peak at the Higgs boson mass mH for e+eZH and e+ee+eH processes, and is expected to be smooth without a resonance structure for other processes in the mass region around 125 GeV.

      Two important measurements of the Higgs boson can be performed from the Mrecoil mass spectrum. The Higgs boson mass can be measured from the peak position of the resonance. The width of the resonance is dominated by the beam energy spread (including ISR effects) and energy/momentum resolution of the detector if the Higgs boson width is only 4.07 MeV as predicted in the SM. The best precision of the mass measurement can be achieved from the leptonic Z+(=e,μ) decays. The height of the resonance is proportional to the Higgs boson production cross section σ(ZH). By fitting the Mrecoil spectrum, the e+eZH event yield, and therefore σ(ZH), can be extracted, independent of the Higgs boson decays. The Higgs boson branching ratios can then be determined by studying Higgs boson decays in selected e+eZH candidates. The recoil mass spectrum has been investigated for both leptonic and hadronic Z boson decays as presented below.

    • 4.1.   Z+

    • The leptonic Z boson decay is ideal for studying the recoil mass spectrum of the e+eZX events. The decay is easily identifiable and the lepton momenta can be precisely measured. Figure 10 shows the reconstructed recoil mass spectra of e+eZX candidates for the Zμ+μ and Ze+e decay modes. The analyses are based on the full detector simulation for the signal events and on the fast detector simulation for background events. They are performed with event selections entirely based on the information of the two leptons, independent of the final states of Higgs boson decays. This approach is essential for the measurement of the inclusive e+eZH production cross section and the model-independent determination of the Higgs boson branching ratios. The SM processes with at least 2 leptons in their final states are considered as backgrounds.

      The event selection of the Zμ+μ decay mode starts with the requirement of a pair of identified muons with opposite charges. Events must have the dimuon invariant mass in the range of 80–100 GeV and the recoil mass between 120 GeV and 140 GeV. The muon pair is required to have its transverse momentum larger than 20 GeV, and its opening angle smaller than 175. A Boosted Decision Tree (BDT) technique is employed to enhance the separation between signal and background events. The BDT is trained using the invariant mass, transverse momentum, polar angle and acollinearity of the dimuon system. Leading background contributions after the selection are from ZZ, WW and Zγ events. As shown in Fig. 10(a), the analysis has a good signal-to-background ratio. The long high-mass tail is largely due to the initial-state radiation.

      Figure 10.  (color online) The inclusive recoil mass spectra of e+eZX candidates for (a) Zμ+μ and (b) Ze+e. No attempt to identify X is made. The markers and their uncertainties (too small to be visible) represent expectations from a CEPC dataset of 5.6ab1, whereas the solid blue curves are the fit results. The dashed curves are the signal and background components.

      Compared to the analysis of the Zμ+μ decay, the analysis of the Ze+e decay suffers from additional and large background contributions from Bhabha scattering and single boson production. A cut based event selection is performed for the Ze+e decay. The electron-positron pair is required to have its invariant mass in the range of 86.2–96.2 GeV and its recoil mass between 120 GeV and 150 GeV. Additional selections based on the kinematic variables of the electron-positron system, the polar angles and the energies of the selected electron and positron, are applied. Events from e+ee+e(γ),e+νW(eˉνW+),e+eZ production are the dominant backgrounds after the selection. The recoil mass distribution of the selected events is shown in Fig. 10(b).

      While event selections independent of the Higgs boson decays are essential for the model-independent measurement of σ(ZH), additional selection criteria using the Higgs boson decay information can, however, be applied to improve the Higgs boson mass measurement. This will be particularly effective in suppressing the large backgrounds from Bhabha scattering and single W or Z boson production for the analysis of the Ze+e decay. These improvements are not implemented in the current study.

    • 4.2.   Zqˉq

    • The recoil mass technique can also be applied to the hadronic Z boson decays (Zqˉq) of the e+eZX candidates. This analysis benefits from a larger Zqˉq decay branching ratio, but suffers from the fact that jet energy resolution is worse than the track momentum and electromagnetic energy resolutions. In addition, ambiguity in selecting jets from the Zqˉq decay, particularly in events with hadronic decays of the Higgs boson, can degrade the analysis performance and also introduce model-dependence to the analysis. Therefore, the measurement is highly dependent on the performance of the particle-flow reconstruction and the jet clustering algorithm.

      Following the same approach as the ILC study [38], an analysis based on the fast simulation has been performed. After the event selection, main backgrounds arise from Zγ and WW production. Compared with the leptonic decays, the signal-to-background ratio is considerably worse and the recoil mass resolution is significantly poorer.

    • 4.3.   Measurements of σ(ZH) and mH

    • Both the inclusive e+eZH production cross section σ(ZH) and the Higgs boson mass mH can be extracted from fits to the recoil mass distributions of the e+eZ+X+/qˉq+X candidates. For the leptonic Z+ decays, the recoil mass distribution of the signal process e+eZH (and e+ee+eH in case of the Ze+e decay) is modeled with a Crystal Ball function [39] whereas the total background is modeled with a polynomial function. As noted above, the recoil mass distribution is insensitive to the intrinsic Higgs boson width should it be as small as predicted by the SM. The Higgs boson mass can be determined with precision of 6.5 MeV and 14 MeV from the Zμ+μ and Ze+e decay modes, respectively. After combining all channels, an uncertainty of 5.9 MeV can be achieved.

      The process e+eZ+Xqˉq+X contributes little to the precision of the mH measurement due to the poor Zqˉq mass resolution, but dominates the sensitivity of the e+eZH cross section measurement because of the large statistics. A simple event counting analysis shows that the expected relative precision on σ(ZH) is 0.61%. In comparison, the corresponding relative precision from the Ze+e and Zμ+μ decays is estimated to be 1.4% and 0.9%, respectively. The combined relative precision of the three measurements is 0.5%. Table 5 summarizes the expected precision on mH and σ(ZH) from a CEPC dataset of 5.6 ab−1.

      Z decay modeΔmH/MeVΔσ(ZH)/σ(ZH)
      e+e141.4%
      μ+μ6.50.9%
      qˉq0.6%
      combination5.90.5%

      Table 5.  Estimated measurement precision for the Higgs boson mass mH and the e+eZH production cross section σ(ZH) from a CEPC dataset of 5.6 ab-1.

    5.   Analyses of individual decay modes
    • Different decay modes of the Higgs boson can be identified through their unique signatures, leading to the measurements of production rates for these decays. For the e+eZH production process in particular, the candidate events can be tagged from the visible decays of the Z bosons, the Higgs boson decays can then be probed by studying the rest of the events. Simulation studies of the CEPC baseline conceptual detector have been performed for the Higgs boson decay modes of Hbˉb/cˉc/gg, HWW, HZZ, HZγ, Hτ+τ, Hμ+μ and Higgs boson to invisible particles (Hinv). The large number of the decay modes of the H, W and Z boson as well as the τ-lepton leads to a very rich variety of event topologies. This complexity makes it impractical to investigate the full list of final states stemming from the Higgs boson decays. Instead, a limited number of final states of individual Higgs boson decay modes have been considered. For some decay modes, the chosen final states may not be the most sensitive ones, but are nevertheless representatives of the decay mode. In most cases, the dominant backgrounds come from the SM diboson production and the single Z production with the initial and final state radiation.

      The studies are optimized for the dominant ZH process, however, the e+eνeˉνeH and e+ee+eH processes are included whenever applicable. The production cross sections of the individual decay modes, σ(ZH)×BR, are extracted. These measurements combined with the inclusive σ(ZH) measurement discussed in Section 4, will allow the determination of the Higgs boson decay branching ratios in a model-independent way.

      In this section, the results of the current CEPC simulation studies of different Higgs boson decay modes are summarized. The studies are based on the CEPC-v1 detector concept and e+e collisions at s=250 GeV. The expected relative precision from a CEPC dataset of 5.6 ab−1 on the product of the ZH cross section and the Higgs boson decay branching ratio, σ(ZH)×BR, are presented. Detailed discussions of individual analyses are beyond the scope of this paper and therefore only their main features are presented. For the study of a specific Higgs boson decay mode, the other decay modes of the Higgs boson often contribute as well. Those contributions are fixed to their SM expectations and are included as backgrounds unless otherwise noted. However, for the combination of all the decay modes under study, they are allowed to vary within the constraints of the measurements of those decays, see Section 6.

      In addition to the invariant and recoil mass, two other mass observables, visible mass and missing mass, are often used in analyses described below. They are defined, respectively, as the invariant mass and recoil mass of all visible experimental objects such as charged leptons, photons and jets, i.e. practically all particles other than neutrinos.

    • 5.1.   Hbˉb/cˉc/gg

    • For a SM Higgs boson with a mass of 125 GeV, nearly 70% of all Higgs bosons decay into a pair of jets: b-quarks (57.7%), c-quarks (2.9%) and gluons (8.6%). While the Hbˉb decay has recently been observed at the LHC [40,41], the Hcˉc and Hgg decays are difficult, if not impossible, to be identified there due to large backgrounds. In comparison, all these three decays can be isolated and studied at the CEPC. The Hcˉc decay is likely the only process for studying Higgs boson coupling to the second-generation quarks at collider experiments. The identification of Hbˉb/cˉc/gg decays poses critical challenges to the CEPC detector performance, particularly the ability to tag b- and c-quark jets against light-flavored jets (u,d,s,g). Thus they are good benchmarks for the design and optimization of the jet flavor tagging performance of the CEPC detector.

      Studies are performed in detail for e+eZH production with the leptonic decays of the Z bosons. The contribution from the Z-fusion process of e+ee+eH is included in the e+eZHe+eH study. The analysis is based on full simulation for the Higgs boson signal samples and fast simulation for the +qˉq background samples. After selecting two leading leptons with opposite charge, the rest of the reconstructed particles are clustered into two jets to form a hadronically decaying Higgs boson candidate, whose invariant mass is required to be between 75 GeV and 150 GeV. The dilepton invariant mass is required to be within 70–110 GeV for the e+e channel and 81–101 GeV for the μ+μ channel. Moreover, the dilepton system must have its transverse momentum in the range 10–90 GeV and its recoil mass between 120 GeV and 150 GeV. In addition, a requirement on the polar angle of the Higgs boson candidate, |cosθH|<0.8, is applied.

      In order to identify the flavors of the two jets of the Higgs boson candidate, variables LB and LC are constructed using information such as those from LCFIPlus jet flavor tagging algorithm. The values of LB (LC) are close to one if both jets are originated from b(c) quarks and are close to zero if both have light-quark or gluon origins. An unbinned maximum likelihood fit to the Mrecoil, LB and LC distributions of candidate events is used to extract the individual signal yields of the Hbˉb, Hcˉc and Hgg decay modes. The total probability density function (PDF) is the sum of signal and background components. For signals, their Mrecoil PDFs are modeled by Crystal Ball functions [39] with small exponential tails. The background PDF is taken as a sum of two components: a background from Higgs boson decays to other final states such as WW and ZZ, and a combinatorial background from other sources, dominated by the e+eZZqˉq production. The background from other Higgs boson decay channels has the same Mrecoil PDF as the signals. The Mrecoil distribution of the combinatorial background is modeled by a second order polynomial. The PDFs of the signal LB and LC distributions are described by two dimensional histograms, taken from the MC simulated events. The LB and LC distributions of both background components are modeled by 2-dimensional histogram PDFs based on the MC simulation. The dilepton recoil mass distributions of the simulated data and the fit results are shown in Fig. 11(a,b). The estimated relative statistical precision of the measurements of σ(ZH)×BR(Hbˉb/cˉc/gg) are listed in Table 6.

      Z decay modeHbˉbHcˉcHgg
      Ze+e1.3%12.8%6.8%
      Zμ+μ1.0%9.4%4.9%
      Zqˉq0.5%10.6%3.5%
      Zνˉν0.4%3.7%1.4%
      combination0.3%3.1%1.2%

      Table 6.  Expected relative precision on σ(ZH)×BR for the Hbˉb,cˉc and gg decays from a CEPC dataset of 5.6ab1.

      Figure 11.  (color online) ZH production with Hbˉb/cˉc/gg: the recoil mass distributions of (a) Ze+e and (b) Zμ+μ; the dijet mass distributions of Higgs boson candidates for (c) Zqˉq and (d)Zνˉν. The markers and their uncertainties represent expectations from a CEPC dataset of 5.6ab1 whereas the solid blue curves are the fit results. The dashed curves are the signal and background components. Contributions from other decays of the Higgs boson are included in the background.

      Table 6 also includes the results of the Zνˉν and Zqˉq decays. For the Zqˉq final state, events are clustered into four jets and the mass information of jet pairs are used to select the Higgs and Z boson candidates. In addition to ZZ, WW is also a major background for this analysis, particularly for the Hcˉc and Hgg decays. As for the Zνˉν final state, events are clustered into two jets are to form the Higgs boson candidate, the invisibly decaying Z boson is inferred from the missing mass of the event. Fits similar to the one used in the analysis of the Z+ channel is subsequently performed to statistically separate the Hbˉb,cˉc and gg decay components. The simulated data and the fitted dijet mass distributions of the Higgs boson candidates are shown in Fig. 11(c,d) for Zqˉq and Zνˉν.

      Combining all Z boson decay modes studied, a relative statistical precision for σ(ZH)×BR of 0.3%, 3.3% and 1.3% can be achieved for the Hbˉb,cˉc and gg decays, respectively.

    • 5.2.   HWW

    • For a 125 GeV SM Higgs boson, the HWW decay has the second largest branching ratio of 21.5% [33]. The sensitivity of the σ(ZH)×BR(HWW) measurement is estimated by combining results from the studies of a few selected final states (Table 7) of the HWW decay of ZH production. SM diboson production is the main background source in all cases.

      ZH final stateprecision
      Ze+eHWWνν,νqˉq2.6%
      Zμ+μHWWνν,νqˉq2.4%
      ZνˉνHWWνqˉq,qˉqqˉq1.5%
      ZqˉqHWWqˉqqˉq1.7%
      combination0.9%

      Table 7.  Expected relative precision on the σ(ZH)×BR(HWW) measurement from a CEPC dataset of 5.6ab1.

      For Z+, the HWW decay final states studied are νν and νqˉq. The ZH candidate events are selected by requiring the dilepton invariant mass in the range of 80–100 GeV and their recoil mass in 120–150 GeV. For Zνˉν, the νqˉq and qˉqqˉq final states are considered for the HWW decay. The presence of neutrinos in the event results in large missing mass, which is required to be in the range of 75–140 (75–150) GeV for the νqˉq (qˉqqˉq) final state. The total visible mass of the event must be in the range of 100–150 GeV for both νqˉq and qˉqqˉq final states. In addition, the total transverse momentum of the visible particles must be in the range of 20–80 GeV. Additional requirements are applied to improve the signal-background separations. For Zqˉq, the HWWqˉqqˉq decay is studied. Candidate events are reconstructed into 6 jets. Jets from Zqˉq, Wqˉq and HWWqˉqqˉq decays are selected by minimizing the χ2 of their mass differences to the masses of Z, W and H boson. Figure 12 shows the visible and missing mass distributions after the selection of the Zνˉν and HWWqˉqqˉq final state.

      Figure 12.  (color online) ZH production with Zνˉν and HWWqˉqqˉq: distributions of (a) the visible mass and (b) the missing mass of selected events. The markers and their uncertainties represent the expected number of events in a CEPC dataset of 5.6ab1, whereas the solid blue curves are the fit results. The dashed curves are the signal and background components. Contributions from other decays of the Higgs boson are included in the background.

      The relative precision on σ(ZH)×BR(HWW) from the decay final states studied is summarized in Table 7.

      The combination of these decay final states leads to a precision of 0.9%. This is likely a conservative estimate as many of the final states of the HWW decay remain to be explored. Including these missing final states will no doubt improve the precision.

    • 5.3.   HZZ

    • The HZZ decay has a branching ratio 2.64% [33] for a 125 GeV Higgs boson in the SM. Events from e+eZH production with the HZZ decay have three Z bosons in their final states with one of them being off-shell. Z bosons can decay to all lepton and quark flavors, with the exception of the top quark. Consequently, the e+eZHZZZ process has a very rich variety of topologies.

      Studies are performed for a few selected ZH final states: Zμ+μ and HZZνˉνqˉq; Zνˉν and HZZ+qˉq. The W and Z boson fusion processes, e+ee+eH and e+eνˉνH, are included in the Z(e+e)H and Z(νˉν)H studies assuming their SM values for the production rates. For the final states studied, the SM ZZ production is the main background. For Zμ+μ and HZZνˉνqˉq, the muon pairs must have their invariant masses between 80–100 GeV, recoil masses between 120–160 GeV and transverse momenta larger than 10 GeV. The jet pairs of the Zqˉq decay candidates are required to have their invariant masses in the range of 10–38 GeV. Figure 13(a) shows the recoil mass distribution of Zμ+μ after the selection. The background is negligible in this final state.

      Figure 13.  (color online) ZH production with HZZ: a) the recoil mass distribution of the μ+μ system for Zμ+μ,HZZνˉνqˉq; b) the invariant mass distribution of the μ+μqˉq system for Zνˉν,HZZμ+μqˉq. The markers and their uncertainties represent expectations from a CEPC dataset of 5.6ab1, whereas the solid blue curves are the fit results. The dashed curves are the signal and background components. Contributions from other decays of the Higgs boson are included in the background.

      The candidates of Zνˉν and HZZ+qˉq are selected by requiring a same-flavor lepton pair and two jets. The total visible energy must be smaller than 180 GeV and the missing mass in the range of 58–138 GeV. Additional requirements are applied on the mass and transverse momenta of the lepton and jet pairs. After the selection, the background is about an order of magnitude smaller than the signal as shown in Fig. 13(b).

      Table 8 summarizes the expected precision on σ(ZH)×BR(HZZ) from the final states considered. The combination of these final states results in a precision of about 4.9%. The sensitivity can be significantly improved considering that many final states are not included in the current study. In particular, the final state of Zqˉq and HZZqˉqqˉq which accounts for a third of all ZHZZZ decay is not studied. Moreover, there are further potential improvements by using multivariate techniques.

      ZH final stateprecision
      Zμ+μHZZνˉνqˉq7.2%
      ZνˉνHZZ+qˉq7.9%
      combination4.9%

      Table 8.  Expected relative precision for the σ(ZH)×BR(HZZ) measurement with an integrated luminosity 5.6ab1.

    • 5.4.   Hγγ

    • The diphoton decay of a 125 GeV Higgs boson has a small branching ratio of 0.23% in the SM due to its origin involving massive W boson and top quark in loops. However, photons can be identified and measured well, thus the decay can be fully reconstructed with a good precision. The decay also serves as a good benchmark for the performance of the electromagnetic calorimeter.

      Studies are performed for the ZH production with Hγγ and four different Z boson decay modes: Zμ+μ,τ+τ,νˉν and qˉq. The Ze+e decay is not considered because of the expected large background from the Bhabha process. The studies are based on the full detector simulation for the Zqˉq decay channel and the fast simulation for the others. Photon candidates are required to have energies greater than 25 GeV and polar angles of |cosθ|<0.9. The photon pair with the highest invariant mass is retained as the Hγγ candidate and its recoil mass must be consistent with the Z boson mass. For the Zμ+μ and Zτ+τ decays, a minimal angle of 8 between any selected photon and lepton is required to suppress backgrounds from final state radiations. After the selection, the main SM background is the e+e(Z/γ)γγ process where the γ's arise from the initial and final state radiation.

      The diphoton mass is used as the final discriminant for the separation of signal and backgrounds. The distribution for the Zνˉν decay mode is shown in Fig. 14. A relative precision of 6.2% on σ(ZH)×BR(Hγγ) can be achieved.

      Figure 14.  (color online) ZH production with Hγγ: the diphoton invariant mass distribution for the Zνˉν decay. The markers and their uncertainties represent expectations from a CEPC dataset of 5.6ab1, whereas the solid blue curve is the fit result. The dashed curves are the signal and background components.

    • 5.5.   HZγ

    • Similar to the Hγγ decay, the HZγ decay in the SM is mediated by W-boson and top-quark loops and has a branching ratio of 0.154%. The HZγ analysis targets the signal process of ZHZZγνˉνqˉqγ, in which one of the Z bosons decays into a pair of quarks and the other decays into a pair of neutrinos.

      The candidate events are selected by requiring exactly one photon with transverse energy between 20–50 GeV and at least two jets, each with transverse energy greater than 10 GeV. The dijet invariant mass and the event missing mass must be within windows of ±12 GeV and ±15 GeV of the Z boson mass, respectively. Additional requirements are applied on the numbers of tracks and calorimeter clusters as well as on the transverse and longitudinal momenta of the Z boson candidates. The backgrounds are dominated by the processes of single boson, diboson, qˉq, and Bhabha production.

      After the event selection, the photon is paired with each of the two Z boson candidates to form Higgs boson candidates and the mass differences, ΔM=MqˉqγMqˉq and ΔM=MνˉνγMνˉν, are calculated. Here the energy and momentum of the νˉν system are taken to be the missing energy and momentum of the event. For signal events, one of the mass differences is expected to populate around MHMZ35 GeV whereas the other should be part of the continuum background. Figure 15 shows the ΔM distribution expected from an integrated luminosity of 5.6 ab1. Modeling the signal distribution of the correct pairing with a Gaussian and the background (including wrong-pairing contribution of signal events) with a polynomial, a likelihood fit results in a relative precision of 13% on σ(ZH)×BR(HZγ).

      Figure 15.  (color online) The distribution of the mass difference ΔM(MqˉqγMqˉq and MνˉνγMνˉν) of the selected e+eZHZZγνˉνqˉqγ candidates expected in a dataset with an integrated luminosity of 5.6ab1. The signal distribution shown is for the correct pairings of the Higgs boson decays.

      This analysis can be improved with further optimizations and the use of multivariate techniques. Other decay modes such as ZHZZγqˉqqˉqγ should further improve the precision on the σ(ZH)×BR(HZγ) measurement.

    • 5.6.   Hτ+τ

    • The Hτ+τ decay has a branching ratio of 6.32% [33] at mH=125 GeV in the SM. The τ-lepton is short-lived and decays to one or three charged pions along with a number of neutral pions. The charged and neutral pions, as well as the two photons from the decay of the latter, can be well resolved and measured by the CEPC detector.

      Simulation studies are performed for e+eZH production with Hτ+τ and Zμ+μ,νˉν and qˉq decays. For Zμ+μ, candidates are first required to have a pair of oppositely charged muons with their invariant mass between 40–180 GeV and their recoil mass between 110–180 GeV. For Zνˉν, candidates are preselected by requiring a missing mass in the range of 65–225 GeV, a visible mass greater than 50 GeV and an event visible transverse momentum between 10–100 GeV. For both decays, a BDT selection is applied after the preselection to identify ditau candidates. The BDT utilizes information such as numbers of tracks and photons and the angles between them. After these selections, the ZH production with the non-tau decays of the Higgs boson is the dominant (>95%) background for Zμ+μ and contributes to approximately 40% of the total background for Zνˉν. The rest of the background in the Zνˉν channel comes from diboson production. For Zqˉq, candidates are required to have a pair of tau candidates with their invariant mass between 20–120 GeV, a pair of jets with their mass between 70–110 GeV and their recoil mass between 100–170 GeV. The main background is again from ZH production originating from the decay modes other than the intended ZHqˉqτ+τ decay. The rest of the background is primarily from ZZ production.

      The final signal yields are extracted from fits to the distributions of variables based on the impact parameters of the leading tracks of the two tau candidates as shown in Fig. 16. Table 9 summarizes the estimated precision on σ(ZH)×BR(Hτ+τ) expected from a CEPC dataset of 5.6ab1 for the three Z boson decay modes studied. The precision from the Ze+e decay mode extrapolated from the Zμ+μ study is also included. The e+ee+eH contribution from the Z fusion process is fixed to its SM value in the extrapolation. In combination, the relative precision of 0.8% is expected for σ(ZH)×BR(Hτ+τ).

      ZH final statePrecision
      Zμ+μHτ+τ2.6%
      Ze+eHτ+τ2.7%
      ZνˉνHτ+τ2.5%
      ZqˉqHτ+τ0.9%
      combination0.8%

      Table 9.  Expected relative precision for the σ(ZH)×BR(Hτ+τ) measurement from a CEPC dataset of 5.6ab1.

      Figure 16.  (color online) Distributions of the impact parameter variable of the leading tracks from the two tau candidates in the Z decay mode: (a) Zμ+μ and (b) Zνˉν. Here D0 and Z0 are the transverse and longitudinal impact parameters, respectively. The markers and their uncertainties represent expectations from a CEPC dataset of 5.6ab1, whereas the solid blue curves are the fit results. The dashed curves are the signal and background components. Contributions from other decays of the Higgs boson are included in the background.

    • 5.7.   Hμ+μ

    • The dimuon decay of the Higgs boson, Hμ+μ, is sensitive to the Higgs boson coupling to the second-generation fermions with a clean final-state signature. In the SM, the branching ratio of the decay is 2.18×104 [33] for mH=125 GeV.

      To estimate CEPC's sensitivity for the Hμ+μ decay, studies are performed for the ZH production with the Z decay modes: Z+, Zνˉν, and Zqˉq. In all cases, the SM production of ZZ is the dominant background source. Candidate events are selected by requiring a pair of muons with its mass between 120–130 GeV and their recoiling mass consistent with the Z boson mass (in the approximate range of 90–93 GeV, depending on the decay mode). Additional requirements are applied to identify specific Z boson decay modes. For Z+, candidate events must have another lepton pair with its mass consistent with mZ. In the case of Zμ+μ, the muon pairs of the Zμ+μ and Hμ+μ decays are selected by minimizing a χ2 based on their mass differences with mZ and mH. For the Zνˉν decay, a requirement on the missing energy is applied. For the Zqˉq decay, candidate events must have two jets with their mass consistent with mZ. To further reduce the ZZ background, differences between the signal and background in kinematic variables, such as the polar angle, transverse momentum and energy of the candidate Hμ+μ muon pair, are exploited. Simple criteria on these variables are applied for the Z+ and Zνˉν decay mode whereas a BDT is used for the Zqˉq decay.

      In all analyses, the signal is extracted through unbinned likelihood fits to the Mμ+μ distributions in the range of 120–130 GeV with a signal-plus-background model. Analytical functions are used model both the signal and background distributions. The signal model is a Crystal Ball function while the background model is described by a second-order Chebyshev polynomial. The dimuon mass distribution combining all Z boson decay modes studied is shown in Fig. 17 with the result of the signal-plus-background fit overlaid. The combined relative precision on the σ(ZH)×BR(Hμ+μ) measurement is estimated to be about 16% for data corresponding to an integrated luminosity of 5.6ab1.

      Figure 17.  (color online) ZH production with the Hμ+μ decay: dimuon invariant mass distribution of the selected Hμ+μ candidates expected from an integrated luminosity of 5.6ab1 at the CEPC. The distribution combines contributions from Z+, Zνˉν, and Zqˉq decays. The markers and their uncertainties represent expectations whereas the solid curve is the fit result. The dashed curves are the signal and background components.

    • 5.8.   The invisible decay of the Higgs boson: Hinv

    • In the SM, the Higgs boson can decay invisibly via HZZνˉννˉν. For a Higgs boson mass of 125 GeV, this decay has a branching ratio of 1.06×103. In many extensions to the SM, the Higgs boson can decay directly to invisible particles [42-45]. In this case, the branching ratio can be significantly enhanced.

      The sensitivity of the BR(Hinv) measurement is studied for the Z+ and Zqˉq decay modes. The HZZνˉννˉν decay is used to model the Hinv decay both in the context of the SM and its extensions. This is made possible by the fact that the Higgs boson is narrow scalar so that its production and the decay can be treated separately. The main background is SM ZZ production with one of the Z bosons decay invisibly and the other decays visibly. Candidate events in the Z+ decay mode are selected by requiring a pair of lepton with its mass between 70–100 GeV and event visible energy in the range 90–120 GeV. Similarly, candidate events in Zqˉq are selected by requiring two jets with its mass between 80–105 GeV and event visible energy in the range 90–130 GeV. Additional selections, including using a BDT to exploit the kinematic differences between signal and background events, are also implemented.

      Table 10 summarizes the expected precision on the measurement of σ(ZH)×BR(Hinv) and the 95% confidence-level (CL) upper limit on BR(Hinv) from a CEPC dataset of 5.6ab1. Subtracting the SM HZZνˉννˉν contribution, a 95% CL upper limit of 0.30% on BRBSMinv, the BSM contribution to the Hinv decay can be obtained.

      ZH final state studiedrelative precision on σ×BRupper limit on BR(Hinv)
      Ze+eHinv339%0.82%
      Zμ+μHinv232%0.60%
      ZqˉqHinv217%0.57%
      combination143%0.41%

      Table 10.  Expected relative precision on σ(ZH)×BR(Hinv) and 95% CL upper limit on BR(Hinv) from a CEPC dataset of 5.6ab1.

    • 5.9.   Measurement of σ(e+eνeˉνeH)×BR(Hbˉb)

    • The W-fusion e+eνeˉνeH process has a cross section of 3.3% of that of the ZH process at s=250 GeV. The product of its cross section and BR(Hbˉb), σ(νˉνH)×BR(Hbˉb), is a key input quantity to one of the two model-independent methods for determining the Higgs boson width at the CEPC, see Section 6. The e+eνˉνHνˉνbˉb process has the same final state as the e+eZHνˉνbˉb process, but has a rate that is approximately one sixth of e+eZHνˉνbˉb at s=250 GeV. The main non-Higgs boson background is the SM ZZ production.

      The Z(νˉν)H background is irreducible and can also interfere with νˉνH in the case of Zνeˉνe. However, the interference effect is not considered in the current study. The νˉνH and Z(νˉν)H contributions can be separated through the exploration of their kinematic differences. While the invariant mass distributions of the two b-quark jets are expected to be indistinguishable, the recoil mass distribution should exhibit a resonance structure at the Z boson mass for Z(νˉν)H and show a continuum spectrum for νˉνH. Furthermore, Higgs bosons are produced with different polar angular distributions, see Fig. 18(a).

      Figure 18.  (color online) Distributions of the bˉb system of the e+eνˉνbˉb events: (a) cosine of the polar angle θ before the event selection and (b) the recoil mass after the event selection. Contributions from e+eνeˉνeH, ZH and other SM processes are shown. The cosθ distributions are normalized to unity and therefore only shapes are compared.

      Candidate events are selected by requiring their visible energies between 105 GeV and 155 GeV, visible masses within 100–135 GeV, and missing masses in the range of 65–135 GeV. The two b-quark jets are identified using the variable LB described in Section 5.1. To separate νˉνH and Z(νˉν)H contributions, a 2-dimensional simultaneous fit in the plane of the recoil mass and polar angle of the bˉb system is performed. The recoil mass resolution is improved through a kinematic fit by constraining the invariant mass of the two b-jets within its resolution to that of the Higgs boson mass. Figure 18(b) shows the recoil mass distribution of the bˉb system after the kinematic fit. A fit to the Mbˉbcosθ distribution with both rates of νˉνH and Z(νˉν)H processes as free parameters leads to relative precision of 2.9% for σ(νˉνH)×BR(Hbˉb) and 0.30% for σ(ZH)×BR(Hbˉb). The latter is consistent with the study of the Hbˉb/cˉc/gg decay described in Section 5.1. Fixing the Z(νˉν)H(bˉb) contribution to its SM expectation yields a relative precision of 2.6% on σ(e+eνeˉνeH)×BR(Hbˉb).

    6.   Combinations of individual measurements

      6.1.   Combined measurements of σ×BR and BR

    • With the measurements of the inclusive cross section σ(ZH) and the cross section times the branching ratio σ(ZH)×BR for the individual Higgs boson decay modes, the branching ratio BR can be extracted. Most of the systematic uncertainties associated with the measurement of σ(ZH) cancel in this procedure. A maximum likelihood fit is used to estimate the precision on the BRs. For a given Higgs boson decay mode, the likelihood has the form:

      L(BR,θ)=Poisson[Nobs|Nexp(BR,θ)]G(θ),

      (3)

      where BR is the parameter of interest and θ represents nuisance parameters associated with systematic uncertainties. The number of observed events is denoted by Nobs, Nexp(BR,θ) is the expected number of events, and G(θ) is a set of constraints on the nuisance parameters due to the systematic uncertainties. The number of expected events is the sum of signal and background events. The number of signal events is calculated from the integrated luminosity, the e+eZH cross section σ(ZH) measured from the recoil method, Higgs boson branching ratio BR, the event selection efficiency ϵ. The number of the expected background events, Nb, is estimated using Monte Carlo samples. Thus:

      Nexp(BR,θ)=Lumi(θlumi)×σZH(θσ)×BR×ϵ(θϵ)+Nb(θb),

      (4)

      where θX (X=lumi,σ, ϵ and b) are the nuisance parameters of their corresponding parameters or measurements. Even with 106 Higgs boson events, statistical uncertainties are expected to be dominant and thus systematic uncertainties are not taken into account for the current studies. The nuisance parameters are fixed to their nominal values.

      For the individual analyses discussed in Section 5, contamination from Higgs boson production or decays other than the one under study are fixed to their SM values for simplicity. In the combination, however, these constraints are removed and the contamination are constrained only by the analyses targeted for their measurements. For example, the Hbˉb/cˉc/gg analysis suffers from contamination from the HWW,ZZqˉqqˉq decays. For the analysis discussed in Section 5.1, these contaminations are estimated from SM. In the combination fit, they are constrained by the HWW and HZZ analyses described in Sections 5.2 and 5.3, respectively. Taking into account these across-channel contaminations properly generally leads to small improvements in precision. For example, the precision on σ(ZH)×BR(HZZ) is improved from 5.3% of the standalone analysis to 4.9% after the combination.

      Table 11 summarizes the estimated precision of Higgs boson property measurements, combining all studies described in this paper. For the leading Higgs boson decay modes, namely bˉb, cˉc, gg, WW, ZZ and τ+τ, percent level precision is expected. The best achievable statistical uncertainties for a dataset of 5.6 ab1 are 0.26% for σ(e+eZH)×BR(Hbˉb) and 0.5% for σ(e+eZH). Even for these measurements, statistics is likely to be the dominant uncertainty source. Systematic uncertainties due to the acceptance of the detector, the efficiency of the object reconstruction/identification, the luminosity and the beam energy determination are expected to be small. The integrated luminosity can be measured with a 0.1% precision, a benchmark already achieved at the LEP [46], and can be potentially improved in the future. The center-of-mass energy will be known better than 1 MeV, resulting negligible uncertainties on the theoretical cross section predictions and experimental recoil mass measurements.

      propertyestimated Precision
      CEPC-v1CEPC-v4
      mH5.9 MeV5.9 MeV
      ΓH2.7%2.8%
      σ(ZH)0.5%0.5%
      σ(νˉνH)3.0%3.2%
      decay modeσ×BRBRσ×BRBR
      Hbˉb0.26%0.56%0.27%0.56%
      Hcˉc3.1%3.1%3.3%3.3%
      Hgg1.2%1.3%1.3%1.4%
      HWW0.9%1.1%1.0%1.1%
      HZZ4.9%5.0%5.1%5.1%
      Hγγ6.2%6.2%6.8%6.9%
      HZγ13%13%16%16%
      Hτ+τ0.8%0.9%0.8%1.0%
      Hμ+μ16%16%17%17%
      BRBSMinv<0.28%<0.30%

      Table 11.  Estimated precision of Higgs boson property measurements for the CEPC-v1 detector concept operating at s=250 GeV. All precision are relative except for mH and BRBSMinv for which ΔmH and 95% CL upper limit are quoted respectively. The extrapolated precision for the CEPC-v4 concept operating at s=240 GeV are included for comparisons, see Section 6.2.

      The estimated precision is expected to improve as more final states are explored and analyses are improved. This is particularly true for ZHZWW and ZHZZZ with complex final states. Therefore, Table 11 represents conservative estimates for many Higgs boson observables.

    • 6.2.   Extrapolation to CEPC-v4

    • As discussed in Section 2.4, the CEPC conceptual detector design has evolved from CEPC-v1 to CEPC-v4 with the main change being the reduction of the solenoidal field from 3.5 Tesla to 3.0 Tesla. In the meantime, the nominal CEPC center-of-mass energy for the Higgs boson factory has been changed from 250 GeV to 240 GeV. The results presented above are based on CEPC-v1 operating at s=250 GeV. However, given the relative small differences in the performance of the two detector concepts and in s, the results for CEPC-v4 operating at s=240 GeV can be estimated through extrapolation taking into account changes in signal and background cross sections as well as track momentum resolution. From 250 GeV to 240 GeV, the e+eZH and e+eνeˉνeH cross sections are reduced, respectively, by approximate 5% and 10% while cross sections for background processes are increased by up to 10%. The change in magnetic field affects the Hμ+μ analysis the most whereas its effect on other analyses are negligible. The extrapolated results for CEPC-v4 at 240 GeV are included in Table 11. In most cases, small relative degradations of a few percent are expected. For the following analyses, the extrapolated results for CEPC-v4 at s=240 GeV are used.

    • 6.3.   Measurement of Higgs boson width

    • The Higgs boson width (ΓH) is of special interest as it is sensitive to BSM physics in Higgs boson decays that are not directly detectable or searched for. However, the 4.07 MeV width predicted by the SM is too small to be measured with a reasonable precision from the distributions of either the invariant mass of the Higgs boson decay products or the recoil mass of the system produced in association with the Higgs boson. In a procedure that is unique to lepton colliders, the width can be determined from the measurements of Higgs boson production cross sections and its decay branching ratios. This is because the inclusive e+eZH cross section σ(ZH) can be measured from the recoil mass distribution, independent of the Higgs boson decays.

      Measurements of σ(ZH) and BR's have been discussed in Sections 4 and 5. By combining these measurements, the Higgs boson width can be calculated in a model-independent way:

      ΓH=Γ(HZZ)BR(HZZ)σ(ZH)BR(HZZ),

      (5)

      where Γ(HZZ) is the partial width of the HZZ decay. Because of the small expected BR(HZZ) value for a 125 GeV Higgs boson (2.64% in the SM), the precision of ΓH is limited by the HZZ analysis statistics. It can be improved including the decay final states with larger branching ratios, e.g. the Hbˉb decay:

      ΓH=Γ(Hbˉb)BR(Hbˉb),

      (6)

      where the partial width Γ(Hbˉb) can be independently extracted from the cross section of the W fusion process e+eνˉνHνˉνbˉb:

      σ(νˉνHνˉνbˉb)Γ(HWW)BR(Hbˉb)

      (7)

      =Γ(Hbˉb)BR(HWW).

      (8)

      Thus, the Higgs boson total width is:

      ΓH=Γ(Hbˉb)BR(Hbˉb)σ(e+eνeˉνeH)BR(HWW),

      (9)

      where BR(Hbˉb) and BR(HWW) are measured from the e+eZH process. The limitation of this method is the precision of the σ(e+eνˉνHνˉνbˉb) measurement.

      The expected precision on ΓH is 5.1% from the measurements of σ(ZH) and BR(HZZ) and is 3.5% from the measurements of σ(νˉνHνˉνbˉb), BR(Hbˉb) and BR(HWW). The quoted precision is dominated by the BR(HZZ) measurement for the former case and the σ(νˉνHνˉνbˉb) measurement for the latter case. The combined ΓH precision of the two measurements is 2.8%, taking into account the correlations between the two measurements.

    7.   Higgs boson coupling measurements
    • To understand the implications of the estimated CEPC precision shown in Table 11 on possible new physics models, the results need to be interpreted in terms of constraints on the parameters in the Lagrangian. This is often referred to as the “Higgs boson coupling measurements”, even though the term can be misleading as discussed below.

      There is no unique way to present the achievable precision on the couplings. Before going into the discussion of the CEPC results, we briefly comment on the choices made here. The goal of the theory interpretation here is to obtain a broad idea of the CEPC sensitivity to the Higgs couplings. The interpretation should be simple with intuitive connections between the models and the experimental observables. Ideally, it should have as little model assumptions as possible. Furthermore, it would be convenient if the results can be interfaced directly with the higher order theoretical calculations, renormalization group equation evolutions, etc. Unfortunately, it is impossible to achieve all of these goals simultaneously.

      Two popular frameworks are, instead, chosen for the interpretation of the CEPC results: the so-called κ-framework [47-56] and the effective field theory (EFT) frameworks [57-77]. As discussed in more detail later, none of these is perfect. But neither of these is wrong as long as one is careful not to over interpret the results. Another important aspect of making projections on the physics potential of a future experiment is that they need to be compared with other experiments. The choices made here follows the most commonly used approaches to facilitate such comparisons. In the later part of this section, Higgs physics potential beyond coupling determination is also discussed.

    • 7.1.   Coupling Fits in the κ-framework

    • The Standard Model makes specific predictions for the Higgs boson couplings to the SM fermions, gSM(Hff), and to the SM gauge bosons, gSM(HVV). In the κ-framework, the potential deviations from the SM are parametrized using the κ parameters defined as:

      κf=g(Hff)gSM(Hff),κV=g(HVV)gSM(HVV),

      (10)

      with κi=1 being the SM prediction. The rates of the Higgs boson production and decays are modified accordingly. For example,

      σ(ZH)=κ2ZσSM(ZH)σ(ZH)×BR(Hff)=κ2Zκ2fκ2ΓσSM(ZH)×BRSM(Hff)

      (11)

      Here κ2Γ(ΓH/ΓSMH) parametrizes the change in the Higgs boson width due to both the coupling modifications and the presence of BSM decays.

      Apart from the tree-level couplings, there are also loop-level couplings of Hgg, Hγγ and HZγ in the SM. In the absence of new physics, these couplings, often referred to as the effective couplings, can be expressed using the κ parameters, described previously. However, new physics states in the loops can alter these couplings. For this reason, three additional κ parameters: κg, κγ and κZγ are introduced to parametrize the potential deviations from the SM for the three effective Hgg, Hγγ and HZγ couplings, respectively.

      It is possible that the Higgs boson can decay directly into new particles or have BSM decays to SM particles. In this case, two types of new decay channels should be distinguished:

      1. Invisible decay. This is a specific channel in which Higgs boson decay into new physics particles that are “invisible” in the detector. Such decays can be specifically searched for. If detected, its rate can be measured. The CEPC sensitivity to this decay channel is quantified by the upper limit on BRBSMinv.

      2. Exotic decays. These include all the other new physics channels. Whether they can be observed, and, if so, to what precision, depends sensitively on the final states. In one extreme, the final states can be very distinct, and the rate can be well measured. In the another extreme, they can be completely swamped by the background. Without the knowledge of the final states and the corresponding expected CEPC sensitivity, the exotic decays are accounted for by treating the Higgs boson width ΓH as an independent free parameter in the interpretation.

      In general, possible deviations of all SM Higgs boson couplings should be considered. However, in the absence of obvious light new physics states with large couplings to the Higgs boson or to other SM particles, a very large deviation (>O(1)) is unlikely. For smaller deviations, the Higgs phenomenology is not sensitive to the deviations of κe, κu, κd and κs as the Higgs boson couplings to these particles are negligible compared with the couplings to other particles [78]. Therefore, these κ parameters are set to unities.

      The CEPC will not be able to directly measure the Higgs boson coupling to top quarks. A deviation of this coupling from its SM value does enter the Hgg, Hγγ and HZγ amplitudes. However, this effect is parametrized by κg, κγ and κZγ already. Therefore, κt is not considered as an independent parameter. For simplicity, previous studies often do not include κZγ in the fit. We will follow this approach here. This leaves the following set of 10 independent parameters:

      κb,κc,κg,κW,κτ,κZ,κγ,κμ,BRBSMinv,ΓH.

      (12)

      Additional assumptions can be made to reduce the number of parameters [33, 79]. For example, it can be reduced to a 7-parameter set, by assuming lepton universality, and the absence of exotic and invisible decays (excluding HZZνˉννˉν) [47, 79]:

      κb,κc,κg,κW,κZ,κγ,κτ=κμ.

      (13)

      This is useful for studies at hadron colliders as the Higgs boson total width cannot be measured with good precision. The interpretation of the CEPC results is also performed using this reduced set to allow for direct comparisons with the expected HL-LHC sensitivity.

      The κi parameters give a simple and intuitive parametrization of the potential deviations. It has a direct connection with the observables shown in Table 11 and does cover many possible modifications of the couplings. However, the κ-framework has its limitations as well. Strictly speaking, it should not be understood as the modification of the SM renormalizable couplings by a multiplicative factor. For instance, some of such κ modifications violate gauge invariance. Higher order corrections in the κ-framework cannot be easily defined. Moreover, the κi parameters do not include all possible effects of new physics either. For example, apart from the overall size, potential new physics can also introduce form factors which can change the kinematics of particles that couple to a particular vertex. Manifestations of this effect can be seen in the EFT analysis. It is useful to compare with the EFT analysis discussed in the next subsection. The EFT relates κZ and κW, and further expands them into three different Lorentz structures. Moreover, some of these higher dimensional HVV couplings are also connected with κγ and anomalous trilinear gauge couplings. The current EFT analysis does not include any new light degrees of freedom, in contrast to the κ-framework with independent parameters BRBSMinv and ΓH. Overall, κ-framework does capture the big picture of the CEPC capability in precision Higgs boson measurements. It is useful as long as its limitations are understood.

      The LHC and especially the HL-LHC will provide valuable and complementary information about the Higgs boson properties. For example, the LHC is capable of directly measuring the tˉtH process [80, 81]. It can also use differential cross sections to differentiate contributions between the top-quark and other heavy particle states in the loop of the Hgg vertex [82-85]. Moreover, it can separate contributions from different operators in the couplings between the Higgs and vector bosons [86]. For the purpose of the coupling fit in the κ-framework, the LHC, with its large statistics, improves the precision of rare processes such as Hγγ. Note that a large portion of the systematic uncertainties intrinsic to a hadron collider can be canceled by taking ratios of measured cross sections. For example, combining the ratio of the rates of ppHγγ and ppHZZ at the LHC and the measurement of the HZZ coupling at the CEPC can significantly improve the κγ precision. These are the most useful inputs from the LHC to combine with the CEPC. Similar studies of combination with the LHC for the ILC can be found in Refs. [49, 50, 72, 87, 88].

      The results of the 10-parameter and the 7-parameter fits for the CEPC with an integrated luminosity of 5.6ab1 are shown in Table 12. The combined precision with the HL-LHC estimates (using fit result number 15 of Ref. [10]) are also shown. The HL-LHC estimates used assume no theoretical uncertainties and thus represent the aggressive HL-LHC projection. It is assumed that the HL-LHC will operate at s=14 TeV and accumulate an integrated luminosity of 3000 fb1. For the 7-parameter fit, the Higgs boson width is a derived quantity, not an independent parameter. Its precision, derived from the precision of the fitted parameters, is 2.4% for the CEPC alone and 1.8% when combined with the HL-LHC projection.

      Relative coupling measurement precision and the 95% CL upper limit on BRBSMinv
      Quantity10-parameter fit7-parameter fit
      CEPCCEPC+HL-LHCCEPCCEPC+HL-LHC
      κb1.3%1.0%1.2%0.9%
      κc2.2%1.9%2.1%1.9%
      κg1.5%1.2%1.5%1.1%
      κW1.4%1.1%1.3%1.0%
      κτ1.5%1.2%1.3%1.1%
      κZ0.25%0.25%0.13%0.12%
      κγ3.7%1.6%3.7%1.6%
      κμ8.7%5.0%
      BRBSMinv<0.30%<0.30%
      ΓH2.8%2.3%

      Table 12.  Coupling measurement precision from the 10-parameter fit and 7-parameter fit described in the text for the CEPC, and corresponding results after combination with the HL-LHC. All the numbers refer to are relative precision except for BRBSMinv for which the 95% CL upper limit are quoted respectively. Some entries are left vacant for the 7-parameter fit as they are not dependent parameters under the fitting assumptions.

      The CEPC Higgs boson property measurements mark a giant step beyond the HL-LHC. First of all, in contrast to the LHC, a lepton collider Higgs factory is capable of measuring the Higgs boson width and the absolute coupling strengths to other particles. A comparison with the HL-LHC is only possible with model dependent assumptions. One of such comparisons is within the framework of the 7-parameter fit, shown in Fig. 19. Even with this set of restrictive assumptions, the advantage of the CEPC is still significant. The measurement of κZ is more than a factor of 10 better. The CEPC can also improve significantly the precision on a set of κ parameters that are affected by large backgrounds at the LHC, such as κb, κc, and κg. Note that this is in comparison with the HL-LHC projection with large systematic uncertainties. Such uncertainties are typically under much better control at lepton colliders. Within this 7-parameter set, the only coupling that the HL-LHC can give a competitive measurement is κγ, for which the CEPC sensitivity is statistically limited. This is also the most valuable input that the HL-LHC can give to the Higgs boson coupling measurements at the CEPC, which underlines the importance of combining the results from these two facilities.

      Figure 19.  (color online) The results of the 7-parameter fit and comparison with the HL-LHC [10]. The projections for the CEPC at 240 GeV with an integrated luminosity of 5.6ab1 are shown. The CEPC results without combination with the HL-LHC input are shown as light red bars. The LHC projections for an integrated luminosity of 300 fb1 are shown in light gray bars.

      The direct search for the Higgs boson decay to invisible particles from BSM physics is well motivated and closely connected to the dark sectors. The CEPC with an integrated luminosity of 5.6 ab1 has a sensitivity of 0.30% expressed in terms of the 95% CL upper limit on the decay branching ratio, as shown in Table 12. The HL-LHC, on the other hand, has a much lower sensitivity of 6%–17% [47] while optimistically may reach 2%–3.5% [94].

      As discussed above, one of the greatest advantages of a lepton collider Higgs factory is its capability to measure the Higgs boson width and couplings in a model-independent way. The projection of such a determination at the CEPC is shown in Fig. 20. For most of the measurements, an order of magnitude improvements over the HL-LHC are expected. The CEPC has a clear advantage in the measurement of κZ. It can also set a much stronger constraint on BRBSMinv.

      Figure 20.  (color online) The 10 parameter fit results for the CEPC at 240 GeV with an integrated luminosity of 5.6ab1 (light red bars) and for the combination with the HL-LHC inputs (dark red bars). All the numbers are relative precision except for BRBSMinv for which the 95% CL upper limit are quoted.

    • 7.2.   Effective-field-theory analysis

    • With the assumption that the scale of new physics is higher than the relevant energy directly accessible at the Higgs factory, the effect of new physics can be characterized within the EFT framework. In this framework, operators with dimension greater than four supplement the SM Lagrangian. Imposing baryon and lepton numbers conservation, all higher dimensional operators are of even dimension:

      LEFT=LSM+ic(6)iΛ2O(6)i+jc(8)jΛ4O(8)j+

      (14)

      where Λ is the new physics scale. The leading new physics effects at the electroweak scale would be from the dimension-six operators. To obtain robust constraints on the Wilson coefficients, ci, a global analysis is required, which includes contributions from all possible dimension-six operators. While a large number of dimension-six operators can be written down, only a subset of them contribute to the Higgs boson processes at the leading order. Among these operators, some are much better constrained by other measurements. It is thus reasonable to focus on the operator that primarily contribute to the Higgs boson processes and therefore reduce the parameter space by making appropriate assumptions, as done in the recent studies for future lepton colliders [68, 70-75]. Following these studies, the CP-violating operators as well as the ones that induce fermion dipole interactions are discarded. At the leading order, CP-violating operators do not have linear contributions to the rates of the Higgs boson processes. While they do contribute to angular observables at the leading order [66, 67], these operators are usually much better constrained by the Electric Dipole Moment (EDM) experiments [95-97], though some parameter space is still available for the CP-violating couplings of the Higgs boson to heavy flavor quarks and leptons [98, 99]. The interference between the fermion dipole interactions with SM terms are suppressed by the fermion masses. The corresponding operators also generate dipole moments, which are stringently constrained, especially for light fermions. For the operators that modify the Yukawa coupling matrices, only the five diagonal ones that correspond to the top, charm, bottom, tau, and muon Yukawa couplings are considered, which are relevant for the Higgs boson measurements at the CEPC.

      Before presenting the projections, some brief comments on the EFT framework are in order. In comparison with the κ-framework, a significant advantage of the EFT is that it gives physical parametrization of potential new physics effects. EFT operators can be used directly in computations. The EFT framework also allow for a natural inclusion of new observables, with possible correlations automatically taken into account. At the same time, the connections with experimental observables are less direct and intuitive. Sometimes, the EFT approach is referred to as model-independent. This is only accurate to a certain extent. It assumes that there are no new light degrees of freedom. In practice, assumptions are often made to simplify the set of EFT operators, as also done here.

      The electroweak precision observables are already tightly constrained by the LEP Z-pole and W mass measurements. The CEPC Z-pole run can further improve the constraints set by the LEP, thanks to the enormous amount (10111012) of Z bosons. The W mass can also be measured with a precision of a few MeVs at the CEPC even without a dedicated WW threshold run. Given that the expected precision of the Z-pole observables and the W mass are much higher than the ones of Higgs boson observables, it is assumed that the former ones are perfectly constrained, which significantly simplifies the analysis. In particular, in a convenient basis all the contact interaction terms of the form HVfˉf can be discarded since they also modify the fermion gauge couplings. Realistic Z-pole constraints have also been considered in recent studies [72, 73, 75], but certain assumptions (such as flavor-universality) and simplifications are made. Future studies with more general frameworks are desired to fully determine the impact of the Z-pole measurements on the Higgs boson analysis.

      The measurements of the triple gauge couplings (TGCs) from the diboson process e+eWW play an important role in the Higgs boson coupling analysis under the EFT framework. Focusing on CP-even dimension-six operators, the modifications to the triple gauge vertices from new physics can be parametrized by three anomalous TGC parameters (aTGCs), conventionally denoted as δg1,Z, δκγ and λZ [100, 101]. Among them, δg1,Z and δκγ are generated by operators that also contribute to the Higgs boson processes. At 240 GeV, the e+eWW process cross section is almost two orders of magnitude larger than that of the e+eZH process. The measurements of the diboson process thus provide strong constraints on the operators that generate the aTGCs. A dedicated study on the TGC measurements at the CEPC is not currently available. A simplified analysis is thus performed to estimate the aTGC sensitivity. The results are shown in Table 13. The analysis roughly follows the methods in Refs. [71, 102]. Only the WW events in the semi-leptonic (electron or muon) channel are used, which are easier to reconstruct and have a sizable branching ratio (29%). In particular, the production polar angle, as well as the two decay angles of the leptonically decaying W boson, can be fully reconstructed, which contain important information on the aTGCs. The two decay angles of the hadronically decaying W boson can only be reconstructed with a two-fold ambiguity. A χ2 fit of the three aTGC parameters to the binned distribution of all five angles is performed, from which the one-sigma interval for each of the three aTGCs as well as the correlations among them are extracted. A signal selection efficiency of 80% is assumed. The effects of systematic uncertainties and backgrounds are not considered, assuming they are under control after the selection cuts.

      CEPC 240 GeV (5.6ab1)
      uncertaintycorrelation matrix
      δg1,ZδκγλZ
      δg1,Z1.2×10310.08−0.90
      δκγ0.9×1031−0.42
      λZ1.3×1031

      Table 13.  The estimated constraints on aTGCs from the measurements of the diboson process (e+eWW) in the semi-leptonic channel at the CEPC 240 GeV with 5.6ab1 data and unpolarized beams. All angular distributions are used in the fit. Only the statistical uncertainties of the signal events are considered, assuming a selection efficiency of 80%.

      Under the assumptions specified above, the dimension-six operator contribution to the Higgs boson and diboson processes consists of a total of twelve degrees of freedom. While all non-redundant bases are equivalent, it is particularly convenient to choose a basis in which the twelve degrees of freedom can be mapped to exactly twelve operators, whereas the rest are removed by the assumptions. Two such bases are considered in this analysis. The first is defined by the set of dimension-six operators in Table 14. Among them, O3W corresponds to the aTGC parameter λZ, OHW and OHB generate the aTGC parameters δg1,Z and δκγ as well as Higgs boson anomalous couplings, while the rest operators can only be probed by the Higgs boson measurements at the leading order. The second basis is the so-called “Higgs basis,” proposed in Ref. [103]. In the Higgs basis, the parameters are defined in terms of the mass eigenstates after the electroweak symmetry breaking, and can be directly interpreted as the size of the Higgs boson couplings. Different from the original Higgs basis, this analysis follows Ref. [71], with the parameters associated with the Hgg, Hγγ and HZγ vertices normalized to the SM one-loop contributions, and denoted as ˉcgg, ˉcγγ and ˉcZγ (as opposed to cgg, cγγ and cZγ in Ref. [103]). The parameter ˉceffgg is further defined to absorb all contributions to the Hgg vertex. With these redefinitions, the set of twelve parameters is given by

      OH=12(μ|H2|)2OGG=g2s|H|2GAμνGA,μν
      OWW=g2|H|2WaμνWa,μνOyu=yu|H|2ˉQL˜HuR+h.c.(ut,c)
      OBB=g2|H|2BμνBμνOyd=yd|H|2ˉQLHdR+h.c.(db)
      OHW=ig(DμH)σa(DνH)WaμνOye=ye|H|2ˉLLHeR+h.c.(eτ,μ)
      OHB=ig(DμH)(DνH)BμνO3W=13!gϵabcWaνμWbνρWcρμ

      Table 14.  A complete set of CP-even dimension-six operators that contribute to the Higgs boson and TGC measurements, assuming there is no correction to the Z-pole observables and the W mass, and also no fermion dipole interaction. GAμν, Waμν and Bμν are the field strength tensors for the SM SU(3)c, SU(2)L and U(1)Y gauge fields, respectively. For Oyu, Oyd and Oye, only the contributions to the diagonal elements of the Yukawa matrices that corresponds to the top, charm, bottom, tau, and muon couplings are considered.

      δcZ,cZZ,cZ,ˉcγγ,ˉcZγ,ˉceffgg,δyt,δyc,δyb,δyτ,δyμ,λZ.

      (15)

      These parameters can be conveniently interpreted as the precision of the Higgs boson couplings analogous to those in the κ-framework. In particular, δcZ, ˉcγγ, ˉcZγ, ˉceffgg and δyt,c,b,τ,μ modifies the sizes of the SM Higgs boson couplings to ZZ, γγ, Zγ, gg and fermions, respectively. cZZ and cZ parametrize the anomalous HZZ couplings:

      L=Hv[cZZg2+g24ZμνZμν+cZg2ZμνZμν]+...,

      (16)

      which are not present in the SM at the leading order. The HWW couplings are written in terms of the parameters shown in Eq. 15 via gauge invariance and are not shown explicitly. For the three aTGC parameters, λZ is kept in Eq. 15, while δg1,Z and δκγ are written in terms of the linear combinations of cZZ, cZ, ˉcγγ and ˉcZγ. The exact definitions of the Higgs basis and the translation to the basis in Table 14 can be found in Ref. [71].

      The estimated precision of all the Higgs boson rate measurements in Section 6 (Table 11), along with their correlations, are included as inputs to the EFT global analysis. In addition, the angular observables of the e+eZH,Z+,Hbˉb channel are included, following the studies in Refs. [66, 67]. This channel is almost background-free after the selection, with a signal selection efficiency of about 40%. For the TGC measurements, the results in Table 13 are used as inputs. The global χ2 is obtained by summing over the χ2 of all the measurements. Due to the high precision of the measurements, it is shown that for all observables, keeping only the linear terms of all EFT parameters gives a very good approximation [71]. This greatly simplifies the fitting procedure, as the total χ2 can be written as

      χ2=ij(cc0)iσ2ij(cc0)j,whereσ2ij(δciρijδcj)1,

      (17)

      where ci's are the EFT parameters, c0's are the corresponding central values which are zero by construction, as the measurements are assumed to be SM-like. The one-sigma uncertainties δci and the correlation matrix ρ can be obtained from σ2ij=2χ2/cicj.

      For comparison, the sensitivities of the LHC 14 TeV with total luminosities of 300 fb1 and 3000 fb1 are also considered. These are combined with the diboson (e+eWW) measurements at the LEP as well as the LHC 8 TeV Higgs boson measurements. For the LHC 14 TeV Higgs boson measurements, the projections by the ATLAS collaboration [10] are used, while the composition of each channel is obtained from Refs. [104-108]. The constraints from the LHC 8 TeV Higgs boson measurements and the diboson measurements at the LEP are obtained directly from Ref. [109]. While the LHC diboson measurements can potentially improve the constraints on aTGCs set by the LEP [64], they are not included in this analysis due to the potential issues related to the validity of the EFT [110, 111] and the assumption that the TGCs dominated by the non-anomalous terms [112].

      The results of the 12-parameter fit at the CEPC are shown in Fig. 21 for the Higgs basis and Fig. 22 for the basis in Table 14. The results from the LHC Higgs boson measurements (both 300 fb1 and 3000 fb1) combined with the LEP diboson measurements are shown in comparison. The results of the combination of the CEPC with the HL-LHC (3000fb1) are also shown in addition to the ones from the CEPC alone. In Fig. 21, the results are shown in terms of the one-sigma precision of each parameter. The LHC results are shown with gray columns with 300fb1 (3000fb1) in light (dark) bars, while the CEPC ones are shown with the red columns, with the CEPC-alone (combination with the HL-LHC) results shown in light (dark) bars.

      Figure 21.  (color online) One-sigma precision of the twelve parameters in the Higgs basis. The first column shows the results from the LHC Higgs boson measurements with 300fb1 (light gray bar) and 3000fb1 (dark gray bar) combined with the LEP diboson (e+eWW) measurement. The second column shows the results from the CEPC with 5.6ab1 data collected at 240 GeV with unpolarized beam. The results from the CEPC alone are shown in light red bars, and the ones from a combination of the CEPC and the HL-LHC are shown in dark red bars. For the LHC fits, δyc is fixed to zero.

      Figure 22.  (color online) The 95%, CL sensitivity to Λ/|ci| for the operators in the basis defined in Table 14. The first two columns show the results from the LHC Higgs boson measurements with 300fb1 and 3000fb1 combined with the LEP diboson (e+eWW) measurement. The last two columns show the results from the CEPC alone and the combination of the CEPC and the HL-LHC (3000fb1). The results of the global fits are shown with dark colored bars. The results from individual fits (by switching on one operator at a time) are shown with light colored bars.For the LHC fits, δyc is fixed to zero.

      In Fig. 22, the results are presented in terms of the sensitivity to Λ/|ci| at 95% CL for each operator as defined in Eq. (14), where Λ is the scale of new physics and ci is the corresponding Wilson coefficient. Four columns are shown separately for the LHC 300fb1, the HL-LHC 3000fb1, the CEPC alone and the CEPC combined with the HL-LHC. The results of the global fits, i.e. simultaneous fits to the 12 parameters, are shown with dark colored bars. The results from individual fits are shown with light colored bars, which are obtained by switching on one operator at a time with the rest fixed to zero.

      It is transparent from Fig. 21 that the CEPC can measure the Higgs boson couplings with precision that is one order of magnitude better than the LHC [10, 11]. For the parameters ˉcγγ, ˉcZγ and δyμ, the clean signal and small branching ratios of the corresponding channels (Hγγ/Zγ/μμ) makes the HL-LHC precision comparable to the CEPC. The combination with the LHC measurements thus provides non-negligible improvements, especially for those parameters. It should be noted that, while δyt modifies the Hgg vertex via the top-quark loop contribution, the CEPC alone cannot discriminate it from the Hgg contact interaction obtained from integrating out a heavy new particle in the loop. The parameter ˉceffgg absorbs both contributions and reflects the overall precision of the Hgg coupling. The combination with the LHC tˉtH measurements can resolve this flat direction. The CEPC measurements, in turn, can improve the constraint on δyt set by the LHC by providing much better constraints on the other parameters that contribute to the tˉtH process. It should also be noted that the measurement of the charm Yukawa coupling is not reported in Ref. [10], while the projection of its constraint has a large variation among different studies and can be much larger than one [113-118]. Therefore, δyc is fixed to be zero for the LHC-only fits, as treating δyc as an unconstrained free parameter generates a flat direction in the fit which makes the overall sensitivity much worse. The CEPC, on the other hand, provides excellent measurements of the charm Yukawa coupling and can constrain δyc to about 2%.

      Regarding the sensitivity to Λ/|ci| in Fig. 22, it is also clear that the CEPC has a significantly better performance than the LHC. If the couplings are naïvely assumed to be of order one (ci1), the Higgs boson measurements at the CEPC would be sensitive to new physics scales at several TeV. While the individual sensitivity to some of the operators at the LHC can be comparable to the CEPC (e.g., OWW and OBB from the measurement of Hγγ), the CEPC sensitivity is much more robust under a global framework. This is due to its comprehensive measurements of both the inclusive ZH cross section and the exclusive rates of many Higgs boson decay channels. Operators OGG and Oyt both contribute to the Hgg vertex. While the CEPC can provide strong constraints on either of them if the other is set to zero, they can only be constrained in a global fit if the tˉtH measurements at the LHC are also included. It is also important to note that the validity of EFT can be a potential issue for the LHC measurements [110]. Depending on the size of the couplings, the inferred bounds on the new physics scale Λ can be comparable with or even smaller than the energy scale probed by the LHC. The CEPC has a smaller center of mass energy and much better precision, which ensures the validity of EFT for most new physics scenarios.

      In Table 15, the numerical results of the global fit are presented for the CEPC in terms of the one-sigma band of the 12 parameters and the correlations among them. The results assume an integrated luminosity of 5.6ab1 at 240 GeV with unpolarized beams, both without and with the combination with the HL-LHC (3000fb1) Higgs boson measurements. With both the one-sigma bounds and the correlation matrix, the corresponding χ2 can be reconstructed, which can be used to derive the constraints in any other EFT basis or any particular model that can be matched to the EFT. This offers a convenient way to study the sensitivity to new physics models, as detailed knowledge of the experimental measurements are not required.

      Higgs basis
      δcZcZZcZˉcγγˉcZγˉceffggδytδycδybδyτδyμλZ
      0.00540.00510.00320.0350.0800.00920.0180.00600.00770.0860.0012
      0.00480.00480.00300.0150.0680.00790.0500.0180.00550.00720.0500.0012
      ci/Λ2[TeV2]of dimension-six operators
      cHcWWcBBcHWcHBcGGcytcyccybcyτcyμc3W
      0.180.0400.0400.130.180.280.0770.111.40.19
      0.160.0350.0350.120.170.00180.820.280.0760.110.830.19

      Table 15.  The one-sigma uncertainties for the 12 parameters from the CEPC (240 GeV, 5.6ab1) in the Higgs basis and the basis of dimension-six operators. For both cases, the upper (lower) row correspond to results without (with) the combination of the HL-LHC Higgs boson measurements.. Note that, without the tˉtH measurements, δyt can not be constrained in a global fit, thus cGG and cyt can not be resolved.

      In the EFT framework, it is explicitly assumed that the Higgs boson width is the sum of all partial widths of its SM decay channels. This is because the EFT expansion in Eq. (14) relies on the assumption that the new physics scale is sufficiently high, while any potential Higgs boson exotic decay necessarily introduces light BSM particles, thus in direct conflict with this assumption. One can nevertheless treat the Higgs boson total width as a free parameter in the EFT global fit and obtain an indirect constraint of it, as done in Ref. [72]. With this treatment, the CEPC can constrain the Higgs boson width to a precision of 1.7% (1.6% if combined with the HL-LHC). This result is significantly better than the one from the 10-parameter coupling fit in Table 12 (3.4%/2.6%). The improvement is mainly because the HWW and HZZ couplings are treated as being independent in the 10-parameter coupling fit, while in the EFT framework they are related to each other under gauge invariance and custodial symmetry. It should also be noted that the Higgs boson width determined using Eqs. (5) and (9) explicitly assumes that the HWW and HZZ couplings are independent of the energy scale. Such an assumption is not valid in the EFT framework with the inclusion of the anomalous couplings.

    • 7.3.   The Higgs boson self-coupling

    • The Higgs boson self-coupling is a critical parameter governing the dynamics of the electroweak symmetry breaking. In the SM, the Higgs boson trilinear and quadrilinear couplings are fixed once the values of the electroweak vacuum expectation value and the Higgs boson mass are known. Any deviation from the SM prediction is thus clear evidence of new physics beyond the SM. The Higgs trilinear coupling is probed at the LHC by the measurement of the di-Higgs production. Current bounds on the Higgs trilinear coupling is at the O(10) level, while the HL-LHC is expected to improve the precision to the level of O(1) [119]. The prospects for extracting the Higgs boson quadrilinear coupling are much less promising, even for a 100 TeV hadron collider [120].

      To measure the di-Higgs production at a lepton collider, a sufficiently large center of mass energy (400 GeV) is required, which is likely to be achieved only at a linear collider. The CEPC, instead, can probe the Higgs boson trilinear coupling via its loop contributions to the single Higgs boson processes. This indirect approach, nevertheless, provides competitive sensitivity, since the loop suppression is compensated by the high precision of the Higgs boson measurements at the CEPC [121]. With a precision of 0.5% on the inclusive ZH cross section at 240 GeV, the Higgs boson trilinear coupling can be constrained to a precision of 35%, assuming all other Higgs boson couplings that contribute to e+eZH are SM-like. While this indirect bound is comparable to the direct ones at linear colliders, it relies on strong assumptions which are only applicable to some specific models.

      A more robust approach is to include all possible deviations on the Higgs boson couplings simultaneously and constrain the Higgs boson trilinear coupling in a global fit. The EFT framework presented in Section 7.2 is ideal for such an analysis. Under this framework, the one-loop contributions of the trilinear Higgs boson coupling to all the relevant Higgs boson production and decay processes are included, following Ref. [74]. The new physics effect is parametrized by the quantity δκλκλ1, where κλ is the ratio of the Higgs boson trilinear coupling to its SM value,

      κλλ3λsm3,λsm3=m2H2v2.

      (18)

      The global fit is performed simultaneously with δκλ and all the 12 EFT parameters defined in Section 7.2. The results are presented in Table 16. The results for the HL-LHC are also shown, which were obtained in Ref. [122] under the same global framework. For the CEPC 240 GeV, the one-sigma bound on δκλ is around ±3, significantly worse than the 35% in the δκλ-only fit. This is a clear indication that it is difficult to resolve the effects of δκλ from other Higgs boson couplings. For the HL-LHC, the reach on δκλ is still dominated by di-Higgs production. However, as a result of the destructive interferences among diagrams, di-Higgs production at the LHC cannot constrain δκλ very well on its positive side, even with the use of differential observables [123]. The combination of the HL-LHC and the CEPC 240 GeV thus provides a non-trivial improvement to the HL-LHC result alone, in particular for the two-sigma bound on the positive side, which is improved from +6.1 to +2.7. This is illustrated in Fig. 23, which plots the profiled χ2 as a function of δκλ for the two colliders.

      bounds on δκλΔχ2=1Δχ2=4
      CEPC 240 GeV (5.6 ab1)[−3.0, +3.1][−5.9, +6.2]
      HL-LHC[−0.9, +1.3][−1.7, +6.1]
      HL-LHC+CEPC 240 GeV[−0.8, +1.0][−1.5, +2.7]

      Table 16.  The Δχ2=1 (one-sigma) and Δχ2=4 (two-sigma) bounds of δκλ for various scenarios, obtained in a global fit by profiling over all other EFT parameters.

      Figure 23.  (color online) Chi-square as a function of δκλ after profiling over all other EFT parameters for the HL-LHC, the CEPC and their combination. The results for the HL-LHC are obtained from Ref. [122].

    • 7.4.   Higgs boson and top-quark couplings

    • Interactions of the Higgs boson with the top quark are widely viewed as a window to new physics beyond the SM. The CEPC potential on the interactions between the Higgs boson and the top quark can be evaluated [77, 124-128] by parametrizing these interactions in terms of dimension-six gauge-invariant operators [129, 130]. This EFT basis enlarges the Higgs basis EFT considered above. Moreover, the CP violation effects in the third generation Yukawa couplings are reflected in the imaginary parts of the Wilson coefficients of operators Oyt and Oyb,

      Δyt=ySMt([Cyt]v32mtΛ2+i[Cyt]v32mtΛ2),

      (19)

      Δyb=ySMt([Cyb]v32mbΛ2+i[Cyb]v32mbΛ2).

      (20)

      In this section, the effect of introducing CP phases in the Yukawa operators in Higgs boson physics is discussed. For more detailed discussion on a complete set of Higgs boson and Top quark operators, see Ref. [124]. The dominant sources of constraints are from Hγγ and Hgg for Oyt, and Hgg and Hbˉb for Oyb. Given that Hgg measurements are sensitive to both operators, a joint analysis of Oyt and Oyb will yield a significantly different result comparing to individual operator analysis. A joint analysis for these two operators in terms of Yukawa coupling strengths and the associated CP phases is performed at the CEPC. The important physics cases for such considerations are highlighted.

      Constraints on the top-quark and bottom-quark Yukawa couplings, including their CP phases, are presented, respectively, in the left and right panels of Fig. 24, respectively. The 68% and 95% CL exclusion bands are shown in dashed and solid lines. The limits for the CEPC are shown in bright black and magenta lines for individual operator analysis and the bright green and yellow shaded regions representing the allowed parameter space at 68% and 95% CL, respectively. The dimmed thick black curves represent the results after turning on both operators OtH and ObH at the same time, using a profile-likelihood method profiling over other parameters. Furthermore, in the left panel the cyan band represents constraints from the HL-LHC tˉtH measurements, red bands are constraints from the CEPC Hgg measurements and blue bands are constraints from the CEPC Hγγ measurements. Similarly, in the right panel, the cyan bands are constraints from Hbˉb and the red bands are constraints from Hgg at the CEPC.

      Figure 24.  (color online) Results for analysis on Cyt and Cyb in the projected allowed regions for modification to the top-quark and bottom-quark Yukawa coupling magnitude and CP phase at 68% and 95% CL. The combined results for the CEPC are shown in black curves. The source of individual constraints for the single operator analysis are labeled correspondingly. For a joint analysis of simultaneous appearance of both Oyt and Oyb operators, the results for the CEPC are shown in the enlarged yellow (95% CL) and green regions (68% CL) with thick brown boundary lines.

      The left panel of Fig. 24 shows that the expected sensitivity on the modification in the magnitude of top-quark Yukawa coupling is around ±3% for the single operator analysis. This is relaxed to [-9.5%, +3%] assuming zero CP phase for the top-quark Yukawa coupling and allowing the bottom-quark Yukawa coupling and its phase to vary freely. The phase of the top-quark Yukawa coupling can be constrained to ±0.16π. This constraint is driven by the Hγγ measurement, where a sizable phase shift will enlarge the Hγγ decay rate via reducing the interference with the SM W boson loop. The constraint on the magnitude of the top-quark Yukawa coupling is driven by the Hgg measurement which is dominated by the top-quark loop contribution. Note that constraints from the Hgg measurement are not constant with respect to the Yukawa coupling magnitude. This is due to the different sizes of the top-quark loop contribution to Hgg through scalar and pseudoscalar couplings. Similarly, as shown in the right panel of Fig. 24 for the bottom-quark Yukawa coupling, the constraint for the magnitude is ±2.5%. For the CP phase, the constraint changes from ±0.47π to zero when the top-quark Yukawa coupling is left free.

    8.   Higgs boson CP test and exotic decays
    • In addition to the studies based on the simulation of the CEPC baseline conceptual detector, the sensitivity of tests on Higgs boson spin/CP properties and in constraining branching ratios of Higgs boson exotic decays are also estimated. These estimates are based on previously published phenomenological studies and are summarized in this section.

    • 8.1.   Tests of Higgs boson spin/CP property

    • The CP properties of the Higgs boson and, more generally, its anomalous couplings to gauge bosons in the presence of BSM physics, can be measured at the CEPC using the e+e(Z)ZHμ+μbˉb process. It is convenient to express the effects of the anomalous couplings in terms of the fractions of events from the anomalous contribution relative to the SM predictions. These fractions are invariant under the independent rescalings of all couplings, see Refs. [131-133].

      Two of the anomalous HZZ coupling measurements are of particular interest at the CEPC: the fraction of the high-order CP-even contribution due to either SM contribution or new physics, fa2, and the fraction of a CP-odd contribution due to new physics, fa3. The following two types of observables can be used to measure these anomalous couplings of the Higgs bosons.

      1) The dependence of the e+eZZH cross section on s is different for different CP property of the Higgs boson [133]. Therefore, measurements of the cross section at several different energies will yield useful information about anomalous HZZ couplings. However this has non-trivial implications to the accelerator design and is not included in this study as a single value of s is assumed for the CEPC operating as a Higgs boson factory.

      2) Angular distributions, cosθ1 or cosθ2 and Φ as defined in Fig. 25. These angles are also sensitive to interference between CP-even and CP-odd couplings. In particular forward-backward asymmetry with respect to cosθ1 or cosθ2 and non-trivial phase in the Φ distributions can lead to an unambiguous interpretation of CP violation.

      Figure 25.  (color online) The Higgs boson production and decay angles for the e+eZZHμ+μbˉb process [133].

      To estimate the sensitivity on the anomalous couplings, a maximum likelihood fit [133] is performed to quantify the compatibility of the observed angular distributions to the theory predictions, including both signal and background processes. In this likelihood fit, the signal probability density functions are taken from analytical predictions that are validated using a dedicated MC program, the JHU generator [131, 132], which incorporates all the anomalous couplings, spin correlations, the interference of all contributing amplitudes. The background probability density function is modeled using simulation based on e+eZZ+bˉb process in MadGraph [134].

      Several thousand statistically independent pseudo-experiments are generated and fitted to estimate the sensitivity to fa2 and fa3, defined as the smallest values that can be measured with 3σ away from 0. All other parameters in the fit, including the number of expected signal and background events, are fixed. The expected sensitivity of 3σ discovery is estimated to be 0.018 for fa2 and 0.007 for fa3. Figure 26(a,b) show the distributions of the fitted values of fa2 and fa3 from the pseudo-experiments expected for fa2=0.018 and fa3=0.008, respectively. A simultaneous fit of fa2 and fa3 is also performed with the 68% and 95% CL contours shown in Fig. 26(c).

      Figure 26.  (color online) Distribution of fitted values of fa2 and fa3 in a large number of pseudo-experiments. In the left and middle plots, only the parameter shown is floated. Other parameters are fixed to their SM values. Right plot: simultaneous fit of non-zero fa2 and fa3, with 68% and 95% confidence level contours shown.

      The sensitivities for fa2 and fa3 are then converted to the corresponding parameters defined for the on-shell HZZ decays, fdeca2 and fdeca3, in order to compare with the sensitivities from the LHC experiments as described in Ref. [133]. The corresponding sensitivities of fdeca2 and fdeca3 are 2×104 and 1.3×104, respectively. The much smaller values in the fdeca2,a3 are due to the much larger m2Z in the e+eZZH process compared to the value from the Higgs boson decays.

      Compared to the ultimate sensitivity of HL-LHC as shown in Ref. [133], the sensitivities in the fa2 and fa3 at the CEPC are better by a factor of 300 and 3. Further improvements can be achieved by exploring kinematics in the Hbˉb decays, including other Z decay final states, and combining with the overall cross-section dependence of the signal as obtained by a threshold scan in s.

    • 8.2.   Higgs boson exotic decays

    • The Higgs boson can be an important portal to new BSM physics. Such new physics could manifest itself through the exotic decays of the Higgs boson if some of the degrees of freedom are light. The Higgs boson BSM decays have a rich variety of possibilities. The two-body decays of the Higgs boson into BSM particles, HX1X2, where the BSM particles Xi are allowed to subsequently decay, are considered here. These decay modes are classified into four cases, schematically shown in Fig. 27. These processes are well-motivated by BSM models such as singlet extensions of the SM, two-Higgs-doublet-models, SUSY models, Higgs portals, gauge extensions of the SM, and so on [48, 135, 136]. In this study, only prompt decays of the BSM particles are considered. For the Higgs boson decaying into long-lived particles, novel search strategies have to be developed in the future, using also the latest advances in the detector development [137].

      Figure 27.  (color online) The topologies of exotic decays of the Higgs boson.

      For the CEPC running at 240 GeV, the most important Higgs boson production mechanism is the e+eZH production. The Z boson with visible decays enables the Higgs boson tagging using the “recoil mass” technique as described in Section 4. A cut around the peak of the recoil mass spectrum would remove the majority of the SM background. Further selection and tagging on the Higgs boson decay products can ensure that the major background would be from the SM decays of the Higgs bosons. The details of these analysis can be found in Ref. [136].

      The set of Higgs boson exotic decays with their projected LHC constraints and limits from the CEPC with an integrated luminosity of  5.6ab1 are summarized in Table 17. For the LHC constraints, both the current limits and projected limits on these exotic decay channels from various references are tabulated. The comparison are performed for particular benchmark points to demonstrate the qualitative difference between the (HL-)LHC and CEPC.

      decay mode95% C.L. limit on BR
      LHCHL-LHCCEPC
      ET0.230.0560.0030
      (bˉb)+ET[0.2]1×10−4
      (jj)+ET4×10−4
      (τ+τ)+ET[1]8×10−5
      bˉb+ET[0.2]2×10−4
      jj+ET5×10−4
      τ+τ+ET8×10−5
      (bˉb)(bˉb)1.7 (0.2)6×10−4
      (cˉc)(cˉc)(0.2)8×10−4
      (jj)(jj)[0.1]2×10−3
      (bˉb)(τ+τ)0.1 [0.15]4×10−4
      (τ+τ)(τ+τ)1.2 [0.2~0.4]2×10−4
      (jj)(γγ)[0.01]1×10−4
      (γγ)(γγ)7×10−34×10−48×10−5

      Table 17.  The current and projected limits on the Higgs boson exotic decay modes for the (HL-)LHC and the CEPC with an integrated luminosity of 5.6ab1, based on the results from Ref. [136]. The second column shows the current LHC results and the projections for 100 fb1 (in parentheses) and 300 fb1 (in square brackets). The available projections for the HL-LHC are listed in the third column. Pairs of objects in parentheses indicate that they are decay products of intermediate resonances, e.g. (bˉb)+ET stands for X+ET(bˉb)+ET where X is an intermediate resonance.

      A selection of results for channels, which are hard to be constrained at the LHC, is shown in Table 17 and Fig. 28. The red bars in the figure correspond to the results using the leptonic Z boson decays that are produced in association with the Higgs bosons. The hadronic decaying Z-boson provides around ten times more statistics and hence further inclusion will improve the results significantly. Based upon the study of the Higgs bosons decaying into WW, ZZ and invisible particles, hadronic decaying Z bosons are conservatively assumed to provide the same upper limits as the leptonic Z boson decays and, hence, improve the limits by around 40% when combined. This extrapolated results are shown in yellow bars.

      Figure 28.  (color online) The 95% C.L. upper limits on selected Higgs boson exotic decay branching ratios at the HL-LHC and the CEPC, based on Ref. [136]. The benchmark parameter choices are the same as in Table 17. The red bars correspond to the results using leptonically decaying spectator Z-boson alone. The yellow bars further include extrapolation with the inclusion of the hadronically decaying Z-bosons. Several vertical lines are drawn in this figure to divide different types of Higgs boson exotic decays.

      In comparison with the HL-LHC, the improvement on the Higgs boson exotic decay branching ratios is significant, varying from one to four orders of magnitude for the channels considered. For the Higgs boson exotic decays into hadronic final states plus missing energy, bˉb+T, jj+T and τ+τ+T, the CEPC improves the HL-LHC sensitivity by three to four orders of magnitude. These significant improvements benefit from the low QCD background and the Higgs boson tagging from the recoil mass reconstruction at the CEPC. Final states with leptons and photons have smaller QCD background at the LHC and therefore the improvements from the CEPC are limited for these final states.

    9.   Implications
    • In this section, we briefly discuss the most important physics implications of the Higgs boson measurements at the CEPC. The measurements of the Higgs boson properties are essential to the understanding of the nature of electroweak symmetry breaking, which remains to be a central and open question. In the SM, it is parametrized by the so-called “Mexican Hat” Higgs potential,

      V(H)=12μ2|H|2+λ4|H|4,

      (21)

      with the vacuum expectation value (VEV) of the Higgs field spontaneously breaking the SU(2)L×U(1)Y gauge symmetry down to U(1)em, and generating masses for the W and Z bosons. With the measurements of the Fermi constant (from muon decay) and the Higgs boson mass, the two parameters in Eq. (21), μ2 and λ, are determined to a very good precision, and thus the SM Higgs potential is fully determined. However, we would like to emphasize that this simplicity is somewhat misleading, as our knowledge of the electroweak symmetry breaking is far from complete. First of all, even though the values of these parameters can be fixed by the experimental measurement, the SM does not contain an explanation of their sizes, and in particular why the electroweak scale appears to be many orders of magnitude smaller than the Planck scale. Furthermore, the Mexican Hat potential as well as the SM itself are based on assumptions, which need to be explicitly tested by experiments before they are established to be correct. In this section, we will focus on the potential of using the precision measurements of Higgs boson properties at the CEPC to address these important questions.

    • 9.1.   Naturalness of the electroweak scale

    • An important question associated with the electroweak symmetry breaking is naturalness. It arises from the need to explain the presence of the weak scale Λweak102 GeV in terms of a more fundamental theory. New physics is necessarily involved in such a theory. The SM by itself cannot answer this question, however, there are many new physics models with the potential to provide an answer. However, a key question for any model of electroweak symmetry breaking, regardless of the model details, is what the scale of new physics is. For instance, if the new physics is the quantum gravity scale, MPlanck=1019 GeV, then an immediate question is how to explain the 17 orders of magnitude difference between it and the electroweak scale. This is often denoted as the naturalness/hierarchy/fine-tuning problem. More generally, the weak scale in any such model can be expressed using dimensional analysis as

      Λ2weakc1M21+c2M22+...,

      (22)

      where MiMNP are the scale of new physics. They are typically the masses of the new physics particles. The ci are numerical coefficients that depend on the details of the model. However, we do note expect them to be very different from order one. Therefore, a large and precise cancellation is needed if MNPΛEW, with the level of tuning proportional to M2NP. The discovery of the spin-zero Higgs boson deepens this mystery. While it is possible to generate a large cancellation by imposing symmetries instead of tuning – one well-known example is the chiral symmetry which protects the masses of the light fermions from receiving large quantum corrections – there is no obvious symmetry that protects the mass of the Higgs boson if it is an elementary scalar particle. To avoid an excessive amount of fine tuning in the theory, the new physics cannot be too heavy, and should preferably be below the TeV scale. This is the main argument for TeV new physics based on naturalness.

      Searching for new physics which leads to a natural electroweak symmetry breaking has been and will continue to be a main part of the physics program at the LHC. Looking for signals from the direct production of the new physics particles, the LHC will probe the new physics scale up to a few TeV. At the same time, as we will show below, the precision measurements at the CEPC can provide competitive sensitivity reach, and has the potential of probing significant higher new physics scales for many scenarios. In addition, the reach of the LHC searches has a strong dependence on the production and decay modes of the new physics particles. The measurements at the CEPC thus provides crucial complementary information and can probe scenarios that are difficult at the LHC. Indeed, the precision measurement of the Higgs boson couplings offers a very robust way of probing new physics related to electroweak symmetry breaking. Any such new physics would necessarily contain particles with sizable couplings to the Higgs boson, which leave their imprints in the Higgs boson couplings. Such a model independent handle is of crucial importance, given the possibility that the new physics could simply be missed by the LHC searches designed based on our wrong expectations.

      In the following, we demonstrate the sensitivity potential to new physics in several broad classes of models, which can address the naturalness of the electroweak symmetry breaking.

      One obvious idea is that the Higgs boson is a composite particle instead of an elementary one. After all, many composite light scalars already exist in nature, such as the QCD mesons. The composite Higgs boson can thus be regarded as a close analogy of the QCD mesons. A light Higgs boson can be naturally obtained if it is implemented as a pseudo-Nambu-Goldstone boson with new dynamics at scale f. Its physics can be described by a chiral Lagrangian similar to that of the low energy QCD. The explicit breaking comes from the couplings which are responsible for the SM fermion masses, and the SM gauge couplings. In this case, the Higgs boson would not unitarize the WW scattering amplitude completely, and its coupling to W and Z will be shifted approximately by

      δκW,δκZO(v2f2).

      (23)

      Therefore, the measurement of κZ provides a strong and robust constraint on f. Taking the results of the 10-parameter fit in Table 12, a precision of 0.21% on κZ implies that values of f below 2.7 TeV are excluded at 95% CL. For specific models, an even stronger bound on f, up to around 5 TeV, can be obtained by exploiting also its contributions to other Higgs boson couplings [139]. The masses of the composite resonances are given by mρgρf, where gρ is the coupling of the new strong interaction, with a size typically much larger than one. This indicates that the CEPC has the potential to probe composite resonance scales much above 10 TeV, which is far beyond the reach of the LHC direct searches. The Higgs boson measurements at the CEPC thus provides a strong and robust test of the idea of naturalness in the composite Higgs boson models. The detailed exclusion regions from the CEPC and the LHC are shown in Fig. 29, in terms of resonance mass mρ, coupling parameter gρL and mixing parameter ξv2/f2.

      Figure 29.  (color online) Limits on the composite Higgs boson model from both direct searches at the LHC and precision measurement at the CEPC. The figures are updated versions of the ones presented in Ref. [138].

      Due to the large Higgs boson coupling to the top quark, arguably the most important particle in addressing the naturalness problem is the top-quark partner. For example, in supersymmetric models (SUSY), the particle mainly responsible for stabilizing the electroweak scale is the scalar top, ˜t (stop). The presence of stop will modify the Higgs boson couplings via a loop contribution, which is most notable for the Hgg and Hγγ couplings since they are also generated at the one-loop level in the SM. The dominant effect is on the Hgg coupling,

      κg1m2t4m2˜t.

      (24)

      The measurement of κg at the CEPC, up to 1% accuracy, will allow us to probe stop mass up to 900 GeV [140, 141]. The situation is also very similar for non-SUSY models with fermionic top-quark partners, with the bounds on the top-quark partner mass being even stronger than the stop one [141]. The more detailed exclusion region in the top-quark partner parameter space is presented in Fig. 30 for both scenarios.

      Figure 30.  (color online) 95% CL Limits on the stop (a) and fermionic top-quark partner (b) from Higgs boson coupling measurements at various current and future collider scenarios, including the CEPC. This figure is reproduced from Ref. [141].

      This gives us another important handle to test the idea of naturalness. We note that, in favorable cases, the search of stop at the LHC run 2 can set a stronger limit on the stop mass. However, this limit depends strongly on the assumption of the mass spectrum of the other superpartners, as well as the relevant decay modes of the stop. As a result, there will still be significant gaps remaining in the parameter space after the upcoming runs of the LHC, and even very light stops cannot be completely excluded. On the other hand, the measurement of the Hgg coupling offers a complementary way of probing the stop that is independent of the decay modes of the stop.

      It is also possible that the top-quark partner does not have the same SM gauge quantum numbers as the top quark. A particularly interesting possibility is that the top-quark partner is a SM singlet. In such scenarios, it is very difficult to search for the top-quark partner at the LHC. It is nontrivial to construct models with SM-singlet top-quark partners that resolve the fine-tuning problem of the electroweak scale [142, 143]. Nevertheless, they offer an extreme example that new physics with a scale of a few hundred GeVs could still be alive after the current and future LHC runs. However, as mentioned earlier, any model that addresses the electroweak naturalness problem would inevitably contain sizable couplings to the Higgs boson. The Higgs boson coupling measurements at the CEPC thus offer an ideal way of testing this type of models, which is very important for making robust arguments on the naturalness problem. As an example, we consider a scalar top-quark partner ϕt with its only interaction to the SM fields given by HHϕtϕt [65, 144]. This interaction contributes to the Higgs propagator at one-loop level, and induces a universal shift to all Higgs boson couplings. The precise measurement of the inclusive ZH cross section imposes a strong constraint on κZ and provides the best constraint on the mass of the top-quark partner, mϕ. As we can see from the left panel of Fig. 31, the CEPC will be able to probe mϕ up to around 700 GeV, giving an non-trivial test of naturalness even in this very difficult scenario. A more concrete model is the so-called “folded SUSY” [143], in which the top-quark partners are scalars analogous to the stops in SUSY. The projected constraints in the folded stop mass plane is shown on in the right panel of Fig. 31, which are at least around 350 GeV for both stops.

      Figure 31.  (color online) (a) The fractional deviation of σZH at the Higgs factory in the scalar singlet top-quark partner model with the HHϕtϕt interaction, reproduced from Ref. [144]. (b) Projected constraints in the folded stop mass plane from the Hγγ coupling measurements at HL-LHC and CEPC, reproduced from Ref. [140]. The dot-dashed red contours indicates the fine-tuning in the Higgs boson mass from the quadratic sensitivity to stop soft terms.

    • 9.2.   Electroweak phase transition

    • The measurement of the properties of the Higgs boson at the LHC has been consistent with the SM so far. At the same time, the nature of the electroweak phase transition remains unknown. While we have a very good knowledge of the sizes of the electroweak VEV and the Higgs boson mass, they only allow to probe a small region of the Higgs potential near the minimum, whereas the global picture of the potential is largely undetermined. This is shown schematically in Fig. 32.

      Figure 32.  (color online) A schematic drawing illustrating the question of the nature of the electroweak phase transition. Left: Our current knowledge of the Higgs potential. Right: Based on our current knowledge, we could not distinguish the SM Mexican Hat potential from an alternative one with more wiggles.

      The remaining region of the Higgs potential is difficult to probe, even with an upgraded LHC. Meanwhile, it has important consequences on the early universe cosmology and the understanding of our observable world. For example, it is crucial in determining whether the electroweak phase transition is of first or second order. The nature of the electroweak phase transition can also be relevant for the matter anti-matter asymmetry in the Universe, as a large class of models of baryogenesis rely on a first order electroweak phase transition. The CEPC has the capability of probing many of these models and potentially revealing the nature of the electroweak phase transition and the origin of baryogenesis.

      It is well known that with a minimal Higgs potential and the SM Higgs sector, the electroweak phase transition is of second order [145]. New physics with sizable couplings to the Higgs boson are needed to make the phase transition a first order one. The measurement of the triple Higgs boson coupling offers an ideal testing ground for these new physics models. Being the third derivative, it carries more information about the global shape of the Higgs potential than the mass. It can also be determined to a reasonable precision at the future colliders, unlike the quartic Higgs boson coupling. Indeed, most models with first order electroweak phase transition predict a triple Higgs boson coupling with large deviations from the SM prediction. This is demonstrated with a simple example in Fig. 33, which shows the deviation in the triple Higgs boson coupling for a generic singlet model. For the model points that produces a first order phase transition, the value of triple Higgs boson coupling indeed covers a wide range and can be different from the SM prediction by up to 100%.

      Figure 33.  (color online) The deviation in the triple Higgs boson coupling in a generic singlet model that could produce first order electroweak phase transition, reproduced from Ref. [146]. Black dots are points where the phase transition is of first order. The parameter g111 is the triple Higgs boson coupling.

      The CEPC could probe the triple Higgs boson coupling via its loop contributions to single Higgs boson processes. As pointed out in Section 7.3, it will have a limited reach to the most general scenario in which all Higgs boson couplings are allowed to deviate from their SM values. An additional run at 350 GeV will help to improve the sensitivity, while a direct measurement using di-Higgs production would have to wait for a future proton-proton collider, or a lepton collider running at much higher energies. However, it should be noted that the model independent approach in Section 7.3 makes no assumption on any possible connection between the triple Higgs boson coupling and other couplings. In practice, to induce large deviation in triple Higgs boson coupling requires the new physics to be close to the weak scale, while the presence of such new physics will most likely induce deviations in other Higgs boson couplings as well, such as the couplings to the electroweak gauge bosons. Without some symmetry or fine tuning, both deviations are expected to come in at the order of v2/M2NP. Such deviations can be probed very well at lepton colliders.

      We will now demonstrate this in the context of models. Instead of a comprehensive survey, we will focus here on some of the simplest possibilities which are also difficult to probe. The minimal model that has been well studied in this class introduces an additional singlet scalar which couples to the Higgs boson [146-151]. The general potential of the Higgs boson and the new scalar S is

      V(H,S)=12μ2|H|2+λ4|H|4+m2SS2+˜aS|H|2+˜κS2|H|2+˜bS3+λSS4.

      (25)

      After integrating out the singlet, it will generate an |H|6 interaction (shown in panel (a) in Fig. 34), which, after electroweak symmetry breaking, leads to a modification of the triple Higgs boson coupling on the order of v2/m2S. At the same time, it will also generate the operator |HH|2. This leads to a wave function renormalization, which gives rises to universal shift of the Higgs boson couplings. In particular, the modification of the HZZ coupling is also of order v2/m2S. We thus expect κZ, which is constrained within 0.25% even with the inclusive ZH measurement alone, to provide the best constraining power on this model. This is explicitly verified with a scan in the model parameter space, shown in Fig.35. The model points with a first order phase transition are projected on the plane of the HZZ and triple Higgs boson couplings. Indeed, for model points with a large deviation in the triple Higgs boson coupling, a sizable deviation in the HZZ coupling is also present. In this model, constraining power of the HZZ coupling measurement at CEPC is almost the same as the triple Higgs boson coupling measurement at a future 100 TeV hadron collider. A more detailed view of the parameter space of the real singlet model is presented in Fig.36. In addition to the deviations in σ(ZH) at CEPC, the sensitivities of the current and future electroweak precision tests are also presented [152]. The σ(ZH) measurement, with a projected precision of 0.5%, indeed provides the best sensitivity in this scenario. We thus conclude that CEPC has an excellent coverage in the full model space that gives a first order electroweak phase transition.

      Figure 35.  (color online) The HZZ and HHH couplings in the real scalar singlet model of Eq. 25. The points in this figure represent models with a first order electroweak phase transition, and are obtained by scanning over the theory space. Points with a first order phase transition are shown in orange, points with a strongly first order phase transition are shown in blue, and points with a strongly first order phase transition that also produces detectable gravitational waves are shown in red. This figure is reproduced from Ref. [153].

      Figure 36.  (color online) The parameter space compatible with a strong first order phase transition (green region) and the deviations in σ(ZH) (dashed red contours) in the real singlet scalar model, reproduced from Ref. [152]. The solid blue region is excluded by current EW and Higgs boson data, and the region with dashed blue lines can be probed by the CEPC Z-pole run.

      Figure 34.  (a) Induced |H|6 couplings after integrating out the singlet. (b) Induced wave function renormalization of the Higgs, |HH|2.

      A more restricted scenario, in which a discrete Z2 symmetry is imposed on the singlet, has also been considered [147, 151]. It is significantly more difficult to achieve a first order electroweak phase transition in this scenario, since the singlet could only modify the Higgs potential at loop levels. To produce the same level of deviation in the Higgs potential, a much stronger coupling between the Higgs boson and the singlet is required, which often exceeds the limits imposed by the requirement of perturbativity. For the same reason, the expected loop induced deviation in the triple Higgs boson coupling is also generically smaller in this case, and is about 10%-15%, as shown in Fig. 37(a). Even in this difficult case, we see in Fig. 37(b) that the expected deviation of the cross section σ(ZH) is about 0.6%. Therefore, the CEPC will see the first evidence of new physics even in this very difficult case. In the more general classes of models, the new physics which modifies the Higgs boson coupling could carry other SM gauge quantum numbers, such as electric charge and/or color. In such cases, there will be significant modifications to the Hgg and Hγγ couplings. One such example is shown in Fig. 37(c), with a 6% deviation in the Hγγ coupling expected in order to obtain a first order phase transition. As shown in Table 12, the combination of CEPC and HL-LHC measurements could constrain κγ to a precision of 1.7%, and would test this scenario with a sensitivity of more than three standard deviations.

      Figure 37.  (color online) Deviations in the triple Higgs, σ(ZH) and Hγγ couplings in models with Z2 symmetry. In each plot, the dashed orange lines are contours of constant deviations in the corresponding quantity, the solid black lines are contours of constant electroweak phase transition strength parameter ξ=v(Tc)/Tc, where v(Tc) is the Higgs VEV at temperature Tc. The shaded region is excluded for producing a color-breaking vacuum.

      In general, the newly discovered Higgs particle could serve as a gateway to new physics. One generic form of the Higgs boson coupling to new physics is the so-called Higgs portal, HHONP, where ONP is an operator composed out of new physics fields. Since HH is the lowest dimensional operator that is consistent with all the symmetries in the SM, it is easy to construct scenarios in which such Higgs portal couplings are the most relevant ones for the low energy phenomenology of new physics. The singlet extended Higgs sector and the scalar top-quark partner, discussed earlier, are special examples of this scenario. In general, the Higgs portal interactions will shift the Higgs boson couplings, and can be thoroughly tested at the CEPC. Moreover, if the new physics is lighter than mH/2, the Higgs portal coupling will lead to new Higgs boson decay channels. We have already seen in Section 8.2 that the CEPC has an excellent capability of probing such exotic decays, and could cover a vast range of decay signals.

    10.   Conclusion
    • The Higgs boson is responsible for the electroweak symmetry breaking. It is the only fundamental scalar particle in the Standard Model observed so far. The discovery of such a particle at the LHC is a major breakthrough on both theoretical and experimental fronts. However, the Standard Model is likely only an effective theory at the electroweak scale. To explore potential new physics at the electroweak scale and beyond, complementary approaches of direct searches at the energy frontier as well as precision measurements will be needed. The current LHC and the planned HL-LHC have the potential to significantly extend its new physics reach and to measure many of the Higgs boson couplings with precision of a few percent.

      However, many new physics models predict Higgs boson coupling deviations at the sub-percent level, beyond those achievable at the LHC. The CEPC complements the LHC and will be able to study the properties of the Higgs boson in great detail with unprecedented precision. Therefore it is capable of unveiling the true nature of this particle. At the CEPC, most Higgs boson couplings can be measured with precision at a sub-percent level. More importantly, the CEPC will able to measure many of the key Higgs boson properties such as the total width and decay branching ratios in a model-independent way, greatly enhancing the coverage of new physics searches. Furthermore, the clean event environment of the CEPC will allow the detailed study of known decay modes and the identification of potential unknown decay modes that are impractical to test at the LHC.

      This paper provides a snapshot of the current studies, many of which are still ongoing. More analyses are needed to fully understand the physics potential of the CEPC. Nevertheless, the results presented here have already built a strong case for the CEPC as a Higgs factory. The CEPC has the potential to characterize the Higgs boson in the same way LEP did with the Z boson, and potentially shed light on new physics.

Reference (153)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return