Processing math: 40%

Precision Higgs physics at the CEPC

Cited by

1. Zhou, C., Zhu, Y.-F., Guo, Q.-Y. et al. Higgs physics research: yesterday, today, and tomorrow | [希格斯物理研究: 昨天、今天、明天][J]. Wuli Xuebao/Acta Physica Sinica, 2024, 73(24): 241401. doi: 10.7498/aps.73.20241207
2. Liang, H., Li, L., Zhu, Y. et al. Measurement of CKM element |Vcb| from W boson decays at the future Higgs factories[J]. Journal of High Energy Physics, 2024, 2024(12): 71. doi: 10.1007/JHEP12(2024)071
3. Li, L.K., Dong, M.Y., Zhang, H.Y. et al. Study of double-sided silicon pixel ladders with low material budget[J]. Nuclear Instruments and Methods in Physics Research, Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 2024. doi: 10.1016/j.nima.2024.169846
4. Ouazghour, B.A., Arhrib, A., Cheung, K. et al. Associated charged Higgs boson production within the 2HDM: E-e+ versus μ-μ+ colliders[J]. Physical Review D, 2024, 110(9): 095026. doi: 10.1103/PhysRevD.110.095026
5. Su, Y.-H., Cai, C., Zeng, Y.-P. et al. Complex scalar dark matter in a new gauged U(1) symmetry with kinetic and direct mixings[J]. Physical Review D, 2024, 110(9): 095014. doi: 10.1103/PhysRevD.110.095014
6. Coloretti, G., Crivellin, A., Mellado, B. Combined explanation of LHC multilepton, diphoton, and top-quark excesses[J]. Physical Review D, 2024, 110(7): 073001. doi: 10.1103/PhysRevD.110.073001
7. Hashemi, M., Roushandel, L. Charged Higgs decay to W±H at a high energy lepton collider[J]. European Physical Journal C, 2024, 84(9): 966. doi: 10.1140/epjc/s10052-024-13316-1
8. Celada, E., Giani, T., ter Hoeve, J. et al. Mapping the SMEFT at high-energy colliders: from LEP and the (HL-)LHC to the FCC-ee[J]. Journal of High Energy Physics, 2024, 2024(9): 91. doi: 10.1007/JHEP09(2024)091
9. Fu, B., Ghoshal, A., King, S.F. et al. Type-I two-Higgs-doublet model and gravitational waves from domain walls bounded by strings[J]. Journal of High Energy Physics, 2024, 2024(8): 237. doi: 10.1007/JHEP08(2024)237
10. Liu, M.-Y., Li, W.-D., Huang, X.-T. et al. Simulation and reconstruction of particle trajectories in the CEPC drift chamber[J]. Nuclear Science and Techniques, 2024, 35(8): 128. doi: 10.1007/s41365-024-01497-z
11. Yang, T.-Z., Zhang, X. Three-point energy correlators in hadronic Higgs boson decays[J]. Physical Review D, 2024, 109(11): 114036. doi: 10.1103/PhysRevD.109.114036
12. Ait Ouazghour, B., Arhrib, A., Cheung, K. et al. Comparison between μ-μ+ and e-e+ colliders for charged Higgs production in the 2HDM[J]. Physical Review D, 2024, 109(11): 115009. doi: 10.1103/PhysRevD.109.115009
13. Chun, E.J., Mandal, S., Padhan, R. Collider imprints of a right handed neutrino magnetic moment operator[J]. Physical Review D, 2024, 109(11): 115002. doi: 10.1103/PhysRevD.109.115002
14. Liang, H., Zhu, Y., Wang, Y. et al. Jet-Origin Identification and Its Application at an Electron-Positron Higgs Factory[J]. Physical Review Letters, 2024, 132(22): 221802. doi: 10.1103/PhysRevLett.132.221802
15. Crivellin, A., Mellado, B. Anomalies in particle physics and their implications for physics beyond the standard model[J]. Nature Reviews Physics, 2024, 6(5): 294-309. doi: 10.1038/s42254-024-00703-6
16. Liu, D., Wang, L.-T., Xie, K.-P. Composite resonances at a 10 TeV muon collider[J]. Journal of High Energy Physics, 2024, 2024(4): 84. doi: 10.1007/JHEP04(2024)084
17. Li, P., Liu, Z., Lyu, K.-F. Higgs boson width and couplings at high energy muon colliders with forward muon detection[J]. Physical Review D, 2024, 109(7): 073009. doi: 10.1103/PhysRevD.109.073009
18. Wang, J., Wang, Y., Zhang, D.-J. Analytic decay width of the Higgs boson to massive bottom quarks at next-to-next-to-leading order in QCD[J]. Journal of High Energy Physics, 2024, 2024(3): 68. doi: 10.1007/JHEP03(2024)068
19. Hong, T.T., Le, V.K., Phuong, L.T.T. et al. Decays of Standard Model–Like Higgs Boson h → γγ, Zγ in a Minimal Left-Right Symmetric Model[J]. Progress of Theoretical and Experimental Physics, 2024, 2024(3): 033B04. doi: 10.1093/ptep/ptae029
20. Dev, P.S.B., Mohapatra, R.N., Zhang, Y. Explanation of the 95 GeV γγ and bb¯ excesses in the minimal left-right symmetric model[J]. Physics Letters, Section B: Nuclear, Elementary Particle and High-Energy Physics, 2024. doi: 10.1016/j.physletb.2024.138481
21. Hu, P., Wang, Y., Du, D. et al. GSHCAL at future e+e− Higgs factories[J]. Nuclear Instruments and Methods in Physics Research, Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 2024. doi: 10.1016/j.nima.2023.168944
22. Fadol, A., Sha, Q., Fang, Y. et al. Application of quantum machine learning in a Higgs physics study at the CEPC[J]. International Journal of Modern Physics A, 2024, 39(1): 2450007. doi: 10.1142/S0217751X24500076
23. Shi, Y., Xia, S., Ji, Q. et al. The optimal design of structural stiffness for the CEPC vertex detector supports based on a multi-objective genetic algorithm[J]. Radiation Detection Technology and Methods, 2024. doi: 10.1007/s41605-024-00498-1
24. Belyaev, A., Deandrea, A., Moretti, S. et al. A fermionic portal to a non-abelian dark sector[J]. Frontiers in Physics, 2024. doi: 10.3389/fphy.2024.1339886
25. Forslund, M., Meade, P. Precision Higgs width and couplings with a high energy muon collider[J]. Journal of High Energy Physics, 2024, 2024(1): 182. doi: 10.1007/JHEP01(2024)182
26. Mosala, K., Sharma, P., Kumar, M. et al. Axion-like particles at future e-p collider[J]. European Physical Journal C, 2024, 84(1): 44. doi: 10.1140/epjc/s10052-024-12401-9
27. Belyaev, A., Deandrea, A., Moretti, S. et al. Fermionic portal to vector dark matter from a new gauge sector[J]. Physical Review D, 2023, 108(9): 095001. doi: 10.1103/PhysRevD.108.095001
28. Ashanujjaman, S., Banik, S., Coloretti, G. et al. SU (2)L triplet scalar as the origin of the 95 GeV excess[J]. Physical Review D, 2023, 108(9): L091704. doi: 10.1103/PhysRevD.108.L091704
29. Sun, K., Wang, Y., Liu, J. et al. A conceptual design of TOF based on MRPC technology for the future electron–positron Higgs factory[J]. Nuclear Instruments and Methods in Physics Research, Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 2023. doi: 10.1016/j.nima.2023.168656
30. Wang, X.-R., Yan, B. Probing the Hgg coupling through the jet charge correlation in Higgs boson decay[J]. Physical Review D, 2023, 108(5): 056010. doi: 10.1103/PhysRevD.108.056010
31. Hashemi, M., Molanaei, M. Heavy neutral 2HDM Higgs boson pair production at CLIC energies[J]. Physical Review D, 2023, 108(3): 035012. doi: 10.1103/PhysRevD.108.035012
32. Barman, R.K., Bélanger, G., Bhattacherjee, B. et al. Is Light Neutralino Thermal Dark Matter in the Phenomenological Minimal Supersymmetric Standard Model Ruled Out[J]. Physical Review Letters, 2023, 131(1): 011802. doi: 10.1103/PhysRevLett.131.011802
33. Ovchynnikov, M., Zhu, J.-Y. Search for the dipole portal of heavy neutral leptons at future colliders[J]. Journal of High Energy Physics, 2023, 2023(7): 39. doi: 10.1007/JHEP07(2023)039
34. Fang, Y., Xin, S. Physics Highlights at CEPC[J]. Nuclear and Particle Physics Proceedings, 2023. doi: 10.1016/j.nuclphysbps.2023.04.011
35. Mou, J.-L., Lü, J.-G., Sun, X.-L. et al. Time of flight detector for charged particle identification based on circular electron-positron collider | [环形正负电子对撞机带电粒子鉴别的飞行时间探测器][J]. Wuli Xuebao/Acta Physica Sinica, 2023, 72(12): 122901. doi: 10.7498/aps.72.20222271
36. Wu, S.-Z., Wu, P., Li, Y.-W. Production of the triply heavy ωccc and ωbbb baryons at e+e- colliders[J]. Physical Review D, 2023, 107(11): 116022. doi: 10.1103/PhysRevD.107.116022
37. Rao, K., Rindani, S.D., Sarmah, P. et al. Polarized Z cross sections in Higgsstrahlung for the determination of anomalous ZZH couplings[J]. Physics Letters, Section B: Nuclear, Elementary Particle and High-Energy Physics, 2023. doi: 10.1016/j.physletb.2023.137847
38. Hashino, K., Ueda, D. SMEFT effects on the gravitational wave spectrum from an electroweak phase transition[J]. Physical Review D, 2023, 107(9): 095022. doi: 10.1103/PhysRevD.107.095022
39. Li, P., Liu, Z., Lyu, K.-F. Heavy neutral leptons at muon colliders[J]. Journal of High Energy Physics, 2023, 2023(3): 231. doi: 10.1007/JHEP03(2023)231
40. Guo, F., Fang, Y., Li, G. et al. Expected measurement precision of the branching ratio of the Higgs boson decaying to the di-photon at the CEPC[J]. Chinese Physics C, 2023, 47(4): 043002. doi: 10.1088/1674-1137/acaa22
41. Antonov, E.S., Drutskoy, A.G. Measurement of Br(H →Zγ) at the 250 GeV ILC[J]. JETP Letters, 2023, 117(3): 177-183. doi: 10.1134/S002136402260330X
42. ZhiFeng, L., YanQing, M. Feynman integral calculation promoting the era of high-precision particle physics | [费曼积分计算助力高精度粒子物理时代][J]. Scientia Sinica: Physica, Mechanica et Astronomica, 2023, 53(10): 100006. doi: 10.1360/SSPMA-2022-0454
43. Wang, Y.. Circular Colliders in China[J]. Springer Proceedings in Physics, 2023. doi: 10.1007/978-3-031-23042-4_21
44. Li, G., Liao, L., Lou, X. et al. Classify the Higgs decays with the PFN and ParticleNet at electron-positron colliders*[J]. Chinese Physics C, 2022, 46(11): 113001. doi: 10.1088/1674-1137/ac7f21
45. Zhu, Y., Cui, H., Ruan, M. The Higgs→ bb¯ , cc¯ , gg measurement at CEPC[J]. Journal of High Energy Physics, 2022, 2022(11): 100. doi: 10.1007/JHEP11(2022)100
46. De Blas, J., Gu, J., Liu, Z. Higgs boson precision measurements at a 125 GeV muon collider[J]. Physical Review D, 2022, 106(7): 073007. doi: 10.1103/PhysRevD.106.073007
47. Li, K., Liu, Y., Yan, Q.-S. et al. Shifts in the BCFW method for QED *[J]. Chinese Physics C, 2022, 46(9): 093111. doi: 10.1088/1674-1137/ac745a
48. Liu, Q., Ran, K., Huang, Y. et al. Expected H → μ+μ-measurement precision with e+e-→ Z(qq)H production[J]. Chinese Physics C, 2022, 46(9): 093001. doi: 10.1088/1674-1137/ac74a9
49. Forslund, M., Meade, P. High precision higgs from high energy muon colliders[J]. Journal of High Energy Physics, 2022, 2022(8): 185. doi: 10.1007/JHEP08(2022)185
50. Algueró, M., Matias, J., Crivellin, A. et al. Unified explanation of the anomalies in semileptonic B decays and the W mass[J]. Physical Review D, 2022, 106(3): 033005. doi: 10.1103/PhysRevD.106.033005
51. Al Ali, H., Arkani-Hamed, N., Banta, I. et al. The muon Smasher's guide[J]. Reports on Progress in Physics, 2022, 85(8): 084201. doi: 10.1088/1361-6633/ac6678
52. Hashemi, M., Nowbakht Ghalati, N. Parameter dependence and analysis of the 2HDM neutral Higgs boson pair production and decay at future lepton colliders[J]. Physics Letters, Section B: Nuclear, Elementary Particle and High-Energy Physics, 2022. doi: 10.1016/j.physletb.2022.137153
53. Bedeschi, F., Gouskos, L., Selvaggi, M. Jet flavour tagging for future colliders with fast simulation[J]. European Physical Journal C, 2022, 82(7): 646. doi: 10.1140/epjc/s10052-022-10609-1
54. Ouhammou, M., Ouali, M., Taj, S. et al. Laser-assisted neutral Higgs-boson pair production in Inert Higgs Doublet Model (IHDM)[J]. Chinese Journal of Physics, 2022. doi: 10.1016/j.cjph.2021.09.012
55. Alasfar, L., de Blas, J., Gröber, R. Higgs probes of top quark contact interactions and their interplay with the Higgs self-coupling[J]. Journal of High Energy Physics, 2022, 2022(5): 111. doi: 10.1007/JHEP05(2022)111
56. Yang, J.-C., Guo, Y.-C., Cai, L.-H. Using a nested anomaly detection machine learning algorithm to study the neutral triple gauge couplings at an e+e− collider[J]. Nuclear Physics B, 2022. doi: 10.1016/j.nuclphysb.2022.115735
57. Rao, K., Rindani, S.D., Sarmah, P. et al. Use of Z polarization in e+e− → ZH to measure the triple-Higgs coupling[J]. Nuclear Physics B, 2022. doi: 10.1016/j.nuclphysb.2021.115649
58. Yue, C.-X., Zhang, H.-Y., Wang, H. Production of axion-like particles via vector boson fusion at future electron-positron colliders[J]. European Physical Journal C, 2022, 82(1): 88. doi: 10.1140/epjc/s10052-022-10007-7
59. Giffin, P., Lewis, I.M., Zheng, Y.-J. Higgs production in association with a dark-Z at future electron positron colliders[J]. Journal of Physics G: Nuclear and Particle Physics, 2022, 49(1): 015003. doi: 10.1088/1361-6471/ac38c1
60. Manzari, C.A.. Flavour and collider signals from a singly charged scalar[J]. Nuovo Cimento della Societa Italiana di Fisica C, 2021, 44(6): 183. doi: 10.1393/ncc/i2021-21183-5
61. Azzi, P., Perez, E. Exploring requirements and detector solutions for FCC-ee[J]. European Physical Journal Plus, 2021, 136(11): 1195. doi: 10.1140/epjp/s13360-021-02141-0
62. Fu, Q., Yang, J.-C., Yue, C.-X. et al. The study of neutral triple gauge couplings in the process e+e− → Zγ including unitarity bounds[J]. Nuclear Physics B, 2021. doi: 10.1016/j.nuclphysb.2021.115543
63. Chen, J., Lu, C.-T., Wu, Y. Measuring Higgs boson self-couplings with 2 → 3 VBS processes[J]. Journal of High Energy Physics, 2021, 2021(10): 99. doi: 10.1007/JHEP10(2021)099
64. Kiuchi, R., Gu, Y., Zhong, M. et al. Physics potential for the H → ZZ ∗ decay at the CEPC[J]. European Physical Journal C, 2021, 81(10): 879. doi: 10.1140/epjc/s10052-021-09653-0
65. Wang, Y., Li, C.S., Huang, F.P. Complementary probe of dark matter blind spots by lepton colliders and gravitational waves[J]. Physical Review D, 2021, 104(5): 053004. doi: 10.1103/PhysRevD.104.053004
66. Crivellin, A., Hoferichter, M., Manzari, C.A. Fermi Constant from Muon Decay Versus Electroweak Fits and Cabibbo-Kobayashi-Maskawa Unitarity[J]. Physical Review Letters, 2021, 127(7): 071801. doi: 10.1103/PhysRevLett.127.071801
67. Hu, Y.L., Sun, C.L., Shen, X.M. et al. Hadronic decays of Higgs boson at NNLO matched with parton shower[J]. Journal of High Energy Physics, 2021, 2021(8): 122. doi: 10.1007/JHEP08(2021)122
68. Crivellin, A., Müller, D., Saturnino, F. Correlating h →μ+μ- To the Anomalous Magnetic Moment of the Muon via Leptoquarks[J]. Physical Review Letters, 2021, 127(2): 021801. doi: 10.1103/PhysRevLett.127.021801
69. Crivellin, A., Manzari, C.A., Algueró, M. et al. Combined Explanation of the Z →b b ¯ Forward-Backward Asymmetry, the Cabibbo Angle Anomaly, and τ →μνν and b →sℓ+ℓ- Data[J]. Physical Review Letters, 2021, 127(1): 011801. doi: 10.1103/PhysRevLett.127.011801
70. Arbuzov, A., Bondarenko, S., Kalinovskaya, L. et al. Electroweak effects in e+e− → ZH process[J]. Symmetry, 2021, 13(7): 1256. doi: 10.3390/sym13071256
71. Crivellin, A., Hoferichter, M. Consequences of chirally enhanced explanations of (g − 2) μ for h → μμ and Z → μμ[J]. Journal of High Energy Physics, 2021, 2021(7): 135. doi: 10.1007/JHEP07(2021)135
72. Liu, J., Wang, X.-P., Xie, K.-P. Searching for lepton portal dark matter with colliders and gravitational waves[J]. Journal of High Energy Physics, 2021, 2021(6): 149. doi: 10.1007/JHEP06(2021)149
73. Gray, H.M.. Future colliders for the high-energy frontier[J]. Reviews in Physics, 2021. doi: 10.1016/j.revip.2021.100053
74. Ge, S.-F., Li, G., Pasquini, P. et al. CP -violating Higgs boson ditau decays: Baryogenesis and Higgs factories[J]. Physical Review D, 2021, 103(9): 095027. doi: 10.1103/PhysRevD.103.095027
75. Abouabid, H., Arhrib, A., Benbrik, R. et al. One-loop radiative corrections to e+e− → Zh0 /H0A0 in the Inert Higgs Doublet Model[J]. Journal of High Energy Physics, 2021, 2021(5): 100. doi: 10.1007/JHEP05(2021)100
76. Li, Z., Wang, Y., Wu, Q.-F. Categorization of two-loop Feynman diagrams in the O(a2) correction to e+e-→ ZH[J]. Chinese Physics C, 2021, 45(5): 053102. doi: 10.1088/1674-1137/abe84d
77. Pandurovic, M.. Measurement of Higgs to WW in the all-jet final state at CEPC vs=250 GeV[J]. Proceedings of Science, 2021.
78. Crivellin, A., Kirk, F., Manzari, C.A. et al. Searching for lepton flavor universality violation and collider signals from a singly charged scalar singlet[J]. Physical Review D, 2021, 103(7): 073002. doi: 10.1103/PhysRevD.103.073002
79. Chen, C.-H., Nomura, T. Electron and muon g − 2, radiative neutrino mass, and ℓ′ → ℓγ in a U(1)e−μ model[J]. Nuclear Physics B, 2021. doi: 10.1016/j.nuclphysb.2021.115314
80. Yang, B., Wang, M., Bi, H. et al. Single production of vectorlike T quark decaying into Wb at the LHC and the future pp colliders[J]. Physical Review D, 2021, 103(3): 036006. doi: 10.1103/PhysRevD.103.036006
81. Crivellin, A., Greub, C., Müller, D. et al. Scalar leptoquarks in leptonic processes[J]. Journal of High Energy Physics, 2021, 2021(2): 182. doi: 10.1007/JHEP02(2021)182
82. Gao, J., Shtabovenko, V., Yang, T.-Z. Energy-energy correlation in hadronic Higgs decays: analytic results and phenomenology at NLO[J]. Journal of High Energy Physics, 2021, 2021(2): 210. doi: 10.1007/JHEP02(2021)210
83. Ma, S., Wang, K., Zhu, J. Higgs decay to light (pseudo)scalars in the semi-constrained NMSSM[J]. Chinese Physics C, 2021, 45(2): 023113. doi: 10.1088/1674-1137/abce4f
84. Bi, Q., Chai, K., Gao, J. et al. Investigating bottom-quark Yukawa interaction at Higgs factory[J]. Chinese Physics C, 2021, 45(2): 023105. doi: 10.1088/1674-1137/abcd2c
85. Ellis, J., He, H.-J., Xiao, R.-Q. Probing new physics in dimension-8 neutral gauge couplings at e + e − colliders[J]. Science China: Physics, Mechanics and Astronomy, 2021, 64(2): 221062. doi: 10.1007/s11433-020-1617-3
86. Han, T., Liu, D., Low, I. et al. Electroweak couplings of the Higgs boson at a multi-TeV muon collider[J]. Physical Review D, 2021, 103(1): 013002. doi: 10.1103/PhysRevD.103.013002
87. Potrebko, A., Polaka, I. Application of Gene Expression Programming in Improving the Event Selection of the Semi-leptonic Top Quark Pair Process[J]. 2021. doi: 10.1109/ITMS52826.2021.9615317
88. Han, T., Li, S., Su, S. et al. Comparative studies of 2HDMs under the Higgs boson precision measurements[J]. Journal of High Energy Physics, 2021, 2021(1): 45. doi: 10.1007/JHEP01(2021)045
89. Tan, Y., Shi, X., Kiuchi, R. et al. Search for invisible decays of the Higgs boson produced at the CEPC[J]. Chinese Physics C, 2020, 44(12): 123001. doi: 10.1088/1674-1137/abb4d8
90. Banerjee, A., Bhattacharyya, G. Probing the Higgs boson through Yukawa force[J]. Nuclear Physics B, 2020. doi: 10.1016/j.nuclphysb.2020.115261
91. Crivellin, A., Müller, D., Saturnino, F. Leptoquarks in oblique corrections and Higgs signal strength: status and prospects[J]. Journal of High Energy Physics, 2020, 2020(11): 94. doi: 10.1007/JHEP11(2020)094
92. Li, L., Li, Y.-Y., Liu, T. et al. Learning physics at future e − e + colliders with machine[J]. Journal of High Energy Physics, 2020, 2020(10): 18. doi: 10.1007/JHEP10(2020)018
93. Draper, P., Kozaczuk, J., Thomas, S. Precision inclusive Higgs physics at e + e − colliders with tracking detectors and without calorimetry[J]. Journal of High Energy Physics, 2020, 2020(9): 174. doi: 10.1007/JHEP09(2020)174
94. Carena, M., Liu, Z., Wang, Y. Electroweak phase transition with spontaneous Z 2-breaking[J]. Journal of High Energy Physics, 2020, 2020(8): 107. doi: 10.1007/JHEP08(2020)107
95. Wang, Z.S., Wang, K. Long-lived light neutralinos at future Z factories[J]. Physical Review D, 2020, 101(11): 115018. doi: 10.1103/PhysRevD.101.115018
96. Ellis, J., Ge, S.-F., He, H.-J. et al. Probing the scale of new physics in the ZZγ coupling at e+e-colliders[J]. Chinese Physics C, 2020, 44(6): 063106. doi: 10.1088/1674-1137/44/6/063106
97. Niu, Y., Zhao, S., Zhang, Y. et al. Design of Sc-ECAL prototype for CEPC and performance of first two layers[J]. Journal of Instrumentation, 2020, 15(5): C05036. doi: 10.1088/1748-0221/15/05/C05036
98. Wang, Z.S., Wang, K. Physics with far detectors at future lepton colliders[J]. Physical Review D, 2020, 101(7): 075046. doi: 10.1103/PhysRevD.101.075046
99. Rindani, S.D., Singh, B. Indirect measurement of triple-Higgs coupling at an electron-positron collider with polarized beams[J]. International Journal of Modern Physics A, 2020, 35(4): 2050011. doi: 10.1142/S0217751X20500116
100. Chiu, W.H., Liu, Z., Wang, L.-T. Probing flavor nonuniversal theories through Higgs physics at the LHC and future colliders[J]. Physical Review D, 2020, 101(3): 035045. doi: 10.1103/PhysRevD.101.035045
101. Arco, F., Heinemeyer, S., Schappacher, C. Higgs-photon production at future e+e− colliders in the complex MSSM[J]. 2020.
102. de Blas, J., Cepeda, M., D’Hondt, J. et al. Higgs Boson studies at future particle colliders[J]. Journal of High Energy Physics, 2020, 2020(1): 139. doi: 10.1007/JHEP01(2020)139
103. Bai, Y., Chen, C.-H., Fang, Y.-Q. et al. Measurements of decay branching fractions of in associated production at the CEPC[J]. Chinese Physics C, 2020, 44(1): 013001. doi: 10.1088/1674-1137/44/1/013001
104. de Blas, J., Durieux, G., Grojean, C. et al. On the future of Higgs, electroweak and diboson measurements at lepton colliders[J]. Journal of High Energy Physics, 2019, 2019(12): 117. doi: 10.1007/JHEP12(2019)117
105. Li, R., Wang, B.-W., Wang, K. et al. Probing the charm Yukawa coupling at future e-p and e+e- colliders[J]. Physical Review D, 2019, 100(5): 053008. doi: 10.1103/PhysRevD.100.053008
106. Gao, J.. Higgs boson decay into four bottom quarks in the SM and beyond[J]. Journal of High Energy Physics, 2019, 2019(8): 174. doi: 10.1007/JHEP08(2019)174
107. Lee, S.J., Park, M., Qian, Z. Probing unitarity violation in the tail of the off-shell Higgs boson in VLVL mode[J]. Physical Review D, 2019, 100(1): 011702. doi: 10.1103/PhysRevD.100.011702
108. Luo, M.-X., Shtabovenko, V., Yang, T.-Z. et al. Analytic next-to-leading order calculation of energy-energy correlation in gluon-initiated Higgs decays[J]. Journal of High Energy Physics, 2019, 2019(6): 37. doi: 10.1007/JHEP06(2019)037

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(2565)
PDF Downloads(90)
Cited by(108)
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

      \delta = c\frac{{{v^2}}}{{M_{{\rm{NP}}}^2}},

      (1)

      where v and M_{\rm NP} 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 {\mathcal{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 \sqrt{s}\sim 240 GeV, near the maximum of the Higgs boson production cross section through the e^+e^-\to ZH 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.6\,{\rm{ab}}^{-1}, 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/C\!P 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 \sqrt{s} = 240 GeV and as a Z factory at \sqrt{s} = 91.2 GeV. It will also perform WW threshold scans around \sqrt{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
      \sqrt{s}/GeV91.2160240
      run time/y217
      instantaneous luminosity/(10^{34}\, {\rm{cm}}^{-2}\, {\rm{s}}^{-1})16–32103
      integrated luminosity/({\rm{ab}}^{-1})8–162.65.6
      Higgs boson yield106
      W boson yield107108
      Z boson yield10^{11}10^{12}108108

      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 \sqrt{s} = 160 GeV, 10^{7} W bosons will be produced in one year. Similarly running at the Z pole around \sqrt{s} = 91.2 GeV (the Z factory), CEPC will produce 10^{11}10^{12} Z bosons. These large samples will enable high precision measurements of the electroweak observables such as A_{FB}^{b}, R_b, 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 \pi^\pm and K^\pm 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 \tau 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,24X_0,5\!\times\! 5 mm2cell with 30 layers
       HCALFe/RPC,6\lambda,10\!\times\! 10 mm2 cell with 40 layers
      performance
       track momentum resolution\Delta(1/p_T)\sim 2\times 10^{-5}(1/GeV)
       impact parameter resolution5\,{\rm \mu m} \oplus 10\,{\rm \mu m} /[(p/{\rm GeV})\, (\sin\theta)^{3/2}]
       ECAL energy resolution\Delta E / E \sim 16%/\sqrt{E/{\rm GeV}} \oplus 1%
       HCAL energy resolution\Delta E / E \sim 60%/\sqrt{E/{\rm GeV}} \oplus 1%

      Table 2.  Basic parameters and performance of the CEPC-v1 detector. The radiation length (X_0) and the nuclear interaction length (\lambda) are measured at the normal incidence. The cell sizes are for transverse readout sensors and the layer numbers are for longitudinal active readouts. The \theta 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^+e^- \to ZH\to q\bar{q}\, b\bar{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^+e^-\to ZH\to q\bar{q}\ b\bar{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 (\ell) 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^+e^-\to ZH\to \ell^+\ell^-\, 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 \tau 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\to \mu^+\mu^- decays as shown in Fig. 4(a).

      Figure 4.  Simulated invariant mass distributions of (a) muon pairs from H\to \mu^+\mu^- and (b) photon pairs from H\to\gamma\gamma, both from the e^+e^-\to ZH process with the Z\to \nu\bar{\nu} decay. The M_{ \mu^+\mu^-} 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_{\gamma\gamma} 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\to\gamma\gamma and H\to Z\gamma decays. They are also important for the reconstruction and measurements of \tau leptons and jets. The H\to\gamma\gamma 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\to\gamma\gamma decays.

    • 2.3.2.   Jets
    • Approximately 70% of Higgs bosons decay directly into jets (b\bar{b}, c\bar{c}, gg) and an additional 22% decay indirectly into final states with jets through the H\to WW^*, 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 W\to q\bar{q}, Z\to q\bar{q} and H\to b\bar{b}/ c\bar{c}/ gg decays from the ZZ\to \nu\bar{\nu}\, q\bar{q}, WW\to\ell\nu\, q\bar{q} and ZH\to \nu\bar{\nu} ( b\bar{b}/c\bar{c}/gg) processes, respectively. Compared with W\to q\bar{q}, the Z\to q\bar{q} and H \to b\bar b/c\bar 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\sigma or better of the the hadronically decaying W and Z bosons.

      Figure 5.  (color online) Distributions of the reconstructed dijet invariant mass for the W\to q\bar{q}, Z\to q\bar{q} and H \to b\bar b/c\bar c/gg decays from, respectively, the WW\to \ell\nu q\bar{q}, ZZ\to \nu\bar{\nu} q\bar{q} and ZH\to \nu\bar{\nu}( b\bar{b}/c\bar{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 Z\to q\bar{q} decays from the Z factory operation is shown in Fig. 6. For an inclusive Z\to q\bar{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 Z\to q\bar{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\times 5 mm2 to 10\times 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 \sqrt{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 \sqrt{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 \sqrt{s}\sim 240-250 GeV are e^+e^-\to ZH (ZH associate production or Higgsstrahlung), e^+e^-\to \nu_e\bar{\nu}_e H (W fusion) and e^+e^-\to e^+e^- H (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^+e^-\to ZH, (b) e^+e^-\to \nu_e\bar{\nu}_e H and (c) e^+e^-\to e^+e^- H.

      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^+e^-\to ZH process reaches its maximum at \sqrt{s}\sim 250 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 \ln^2(s/M_V^2). 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^+e^-\to ZH and e^+e^-\to ( e^+e^-/ \nu\bar{\nu}) H as functions of \sqrt{s} for a 125 GeV SM Higgs boson. The vertical indicates \sqrt{s} = 250 GeV, the energy assumed for most of the studies summarized in this paper.

      decay modebranching ratiorelative uncertainty
      H\to b\bar{b}57.7%+3.2%, −3.3%
      H\to c\bar{c}2.91%+12%, −12%
      H\to \tau^+\tau^-6.32%+5.7%, −5.7%
      H\to \mu^+\mu^-2.19×10−4+6.0%, −5.9%
      H\to WW^*21.5%+4.3%, −4.2%
      H\to ZZ^*2.64%+4.3%, −4.2%
      H\to \gamma\gamma2.28×10−3+5.0%, −4.9%
      H\to Z\gamma1.53×10−3+9.0%, −8.8%
      H\to gg8.57%+10%, −10%
      \Gamma_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 \sqrt{s} = 250 GeV are listed in Table 4. Because of the interference effects between e^+e^-\to ZH and e^+e^-\to \nu_e\bar{\nu}_e H for the Z\to\nu_e\bar{\nu}_e decay and between e^+e^-\to ZH and e^+e^-\to e^+e^- H for the Z\to e^+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^+e^-\to ZH cross section shown is from Fig. 7(a) only whereas the e^+e^-\to \nu_e\bar{\nu}_e H and e^+e^-\to e^+e^- H cross sections include contributions from their interferences with the e^+e^-\to ZH process.

      processcross sectionevents in 5.6 ab−1
      Higgs boson production, cross section in fb
      e^+e^-\to ZH204.71.15\times{10^{6}}
      e^+e^-\to \nu_e\bar{\nu}_e H6.853.84\times{10^{4}}
      e^+e^-\to e^+e^- H0.633.53\times{10^{3}}
      total212.11.19\times{10^{6}}
      background processes, cross section in pb
      e^+e^-\to e^+e^-\,(\gamma)(Bhabha)8504.5\times{10^{9}}
      e^+e^-\to q\bar{q}\,(\gamma)50.22.8\times{10^{8}}
      e^+e^-\to \mu^+\mu^-\,(\gamma)[or\tau^+\tau^-\,(\gamma)]4.402.5\times{10^{7}}
      e^+e^-\to WW15.48.6\times{10^{7}}
      e^+e^-\to ZZ1.035.8\times{10^{6}}
      e^+e^-\to e^+e^- Z4.732.7\times{10^{7}}
      e^+e^-\to e^+\nu W^-/e^-\bar{\nu}W^+5.142.9\times{10^{7}}

      Table 4.  Cross sections of the Higgs boson production and other SM processes at \sqrt{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\theta|<0.99. Photons, if any, are required to have E_\gamma>0.1 GeV and |\cos\theta_{e^\pm\gamma}|<0.99.

      The CEPC as a Higgs boson factory is designed to deliver a combined integrated luminosity of 5.6\,{\rm{ab}}^{-1} to two detectors in 7 years. Over 10^6 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^+e^-\to ZH 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^+e^-\to ZH production cross section, \sigma(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^+e^-\to e^+e^- (Bhabha scattering), e^+e^-\to Z\gamma (initial-state radiation return), e^+e^-\to WW/ZZ (diboson) as well as the single boson production of e^+e^-\to e^+e^- Z and e^+e^-\to e^+\nu W^-/e^-\bar{\nu}W^+. Their cross sections and expected numbers of events for an integrated luminosity of 5.6\,{\rm{ab}}^{-1} at \sqrt{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^+e^-\to e^+\nu_e W^-\to e^+\nu_e e^-\bar{\nu}_e and e^+e^-\to e^+e^- Z\to e^+e^-\nu_e\bar{\nu}_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 \sqrt{s} obtained from the Whizard program [34]. The calculations include initial-state radiation (ISR). The W and Z fusion processes refer to e^+e^- \to \nu\bar{\nu} H and e^+e^-\to e^+e^- H production, respectively. Their numerical values at \sqrt{s} = 250 GeV can be found in Table 4.

      Along with 1.2\times 10^{6} Higgs boson events, 5.8\times{10^6} ZZ, 8.6\times{10^7} WW and 2.8\times{10^8} q\bar{q}(\gamma) 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 \sqrt{s} and zero total momentum:

      M_{{\rm{recoil}}}^2 = {(\sqrt s - {E_{ff}})^2} - p_{ff}^2 = s - 2{E_{ff}}\sqrt s + m_{ff}^2.

      (2)

      Here, E_{ff}, p_{ff} and m_{ff} are, respectively, the total energy, momentum and invariant mass of the fermion pair. The M_{\rm{recoil}} distribution should show a peak at the Higgs boson mass m_H for e^+e^-\to ZH and e^+e^-\to e^+e^- H 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 M_{\rm{recoil}} 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\to \ell^+\ell^-\, (\ell = e,\mu) decays. The height of the resonance is proportional to the Higgs boson production cross section \sigma(ZH). By fitting the M_{\rm{recoil}} spectrum, the e^+e^-\to ZH event yield, and therefore \sigma(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^+e^-\to ZH candidates. The recoil mass spectrum has been investigated for both leptonic and hadronic Z boson decays as presented below.

    • 4.1.   Z\to \ell^+\ell^-

    • The leptonic Z boson decay is ideal for studying the recoil mass spectrum of the e^+e^-\to ZX events. The decay is easily identifiable and the lepton momenta can be precisely measured. Figure 10 shows the reconstructed recoil mass spectra of e^+e^-\to ZX candidates for the Z\to \mu^+\mu^- and Z\to e^+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^+e^-\to ZH 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\to \mu^+\mu^- 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^\circ. 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\gamma 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^+e^-\to ZX candidates for (a) Z\to \mu^+\mu^- and (b) Z\to e^+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.6\,{\rm{ab}}^{-1}, 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\to \mu^+\mu^- decay, the analysis of the Z\to e^+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 Z\to e^+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^{+}e^{-}\to e^+e^-(\gamma), e^+\nu W^-\,(e^-\bar{\nu}W^+),\, e^+e^- Z 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 \sigma(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 Z\to e^+e^- decay. These improvements are not implemented in the current study.

    • 4.2.   Z\to q\bar{q}

    • The recoil mass technique can also be applied to the hadronic Z boson decays (Z\to q\bar{q}) of the e^+e^-\to ZX candidates. This analysis benefits from a larger Z\to q\bar{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 Z\to q\bar{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\gamma 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 \sigma(ZH) and m_H

    • Both the inclusive e^+e^-\to ZH production cross section \sigma(ZH) and the Higgs boson mass m_H can be extracted from fits to the recoil mass distributions of the e^+e^-\to Z+X \to\ell^+\ell^-/ q\bar{q} + X candidates. For the leptonic Z\to \ell^+\ell^- decays, the recoil mass distribution of the signal process e^+e^-\to ZH (and e^+e^-\to e^+e^- H in case of the Z\to e^+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\to \mu^+\mu^- and Z\to e^+e^- decay modes, respectively. After combining all channels, an uncertainty of 5.9 MeV can be achieved.

      The process e^+e^-\to Z+X\to q\bar{q}+X contributes little to the precision of the m_H measurement due to the poor Z\to q\bar{q} mass resolution, but dominates the sensitivity of the e^+e^-\to ZH cross section measurement because of the large statistics. A simple event counting analysis shows that the expected relative precision on \sigma(ZH) is 0.61%. In comparison, the corresponding relative precision from the Z\to e^+e^- and Z\to \mu^+\mu^- 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 m_H and \sigma(ZH) from a CEPC dataset of 5.6 ab−1.

      Z decay mode\Delta m_H/MeV\Delta\sigma(ZH)/\sigma(ZH)
      e^+e^-141.4%
      \mu^+\mu^-6.50.9%
      q\bar{q}0.6%
      combination5.90.5%

      Table 5.  Estimated measurement precision for the Higgs boson mass m_H and the e^+e^-\to ZH production cross section \sigma(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^+e^-\to ZH 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 H \to b\bar b/c\bar c/gg, H\to WW^*, H\to ZZ^*, H\to Z\gamma, H\to \tau^+\tau^-, H\to \mu^+\mu^- and Higgs boson to invisible particles (H\to {\rm{inv}}). The large number of the decay modes of the H, W and Z boson as well as the \tau-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^-\to \nu_e\bar{\nu}_e H and e^+e^-\to e^+e^- H processes are included whenever applicable. The production cross sections of the individual decay modes, \sigma(ZH)\times {\rm{BR}}, are extracted. These measurements combined with the inclusive \sigma(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 \sqrt{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, \sigma(ZH)\times {\rm 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.   H \to b\bar b/c\bar 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 H\to b\bar{b} decay has recently been observed at the LHC [40,41], the H\to c\bar{c} and H\to gg 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 H\to c\bar{c} decay is likely the only process for studying Higgs boson coupling to the second-generation quarks at collider experiments. The identification of H \to b\bar b/c\bar 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^+e^-\to ZH production with the leptonic decays of the Z bosons. The contribution from the Z-fusion process of e^+e^-\to e^+e^- H is included in the e^+e^-\to ZH\to e^+e^- H study. The analysis is based on full simulation for the Higgs boson signal samples and fast simulation for the \ell^+\ell^- q\bar{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 \mu^+\mu^- 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, \vert \cos\theta_{H}\vert < 0.8, is applied.

      In order to identify the flavors of the two jets of the Higgs boson candidate, variables L_B and L_C are constructed using information such as those from LCFIPlus jet flavor tagging algorithm. The values of L_{B} (L_{C}) 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 M_{\rm recoil}, L_B and L_C distributions of candidate events is used to extract the individual signal yields of the H\to b\bar{b}, H\to c\bar{c} and H\to gg decay modes. The total probability density function (PDF) is the sum of signal and background components. For signals, their M_{\rm recoil} 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^+e^-\to ZZ\to \ell\ell q\bar{q} production. The background from other Higgs boson decay channels has the same M_{\rm recoil} PDF as the signals. The M_{\rm recoil} distribution of the combinatorial background is modeled by a second order polynomial. The PDFs of the signal L_B and L_C distributions are described by two dimensional histograms, taken from the MC simulated events. The L_B and L_C 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 \sigma(ZH)\times {\rm{BR}}( H \to b\bar b/c\bar c/gg) are listed in Table 6.

      Z decay modeH \to b\bar{b}H \to c\bar{c}H \to gg
      Z\to e^+e^-1.3%12.8%6.8%
      Z\to \mu^+\mu^-1.0%9.4%4.9%
      Z\to q\bar{q}0.5%10.6%3.5%
      Z\to \nu\bar{\nu}0.4%3.7%1.4%
      combination0.3%3.1%1.2%

      Table 6.  Expected relative precision on \sigma(ZH)\times{\rm BR} for the H\to b\bar{b},\,c\bar{c} and gg decays from a CEPC dataset of 5.6\,{\rm{ab}}^{-1}.

      Figure 11.  (color online) ZH production with H \to b\bar b/c\bar c/gg: the recoil mass distributions of (a) Z\to e^+e^- and (b) Z\to \mu^+\mu^-; the dijet mass distributions of Higgs boson candidates for (c) Z\to q\bar{q} and (d)Z\to \nu\bar{\nu}. The markers and their uncertainties represent expectations from a CEPC dataset of 5.6\,{\rm{ab}}^{-1} 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\to \nu\bar{\nu} and Z\to q\bar{q} decays. For the Z\to q\bar{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 H\to c\bar{c} and H\to gg decays. As for the Z\to \nu\bar{\nu} 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\to \ell^+\ell^- channel is subsequently performed to statistically separate the H\to b\bar{b},\, c\bar{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 Z\to q\bar{q} and Z\to \nu\bar{\nu}.

      Combining all Z boson decay modes studied, a relative statistical precision for \sigma(ZH)\times\mathrm{BR} of 0.3%, 3.3% and 1.3% can be achieved for the H\to b\bar{b},\, c\bar{c} and gg decays, respectively.

    • 5.2.   H\to WW^*

    • For a 125 GeV SM Higgs boson, the H\to WW^* decay has the second largest branching ratio of 21.5% [33]. The sensitivity of the \sigma(ZH)\times{\rm BR}( H\to WW^*) measurement is estimated by combining results from the studies of a few selected final states (Table 7) of the H\to WW^* decay of ZH production. SM diboson production is the main background source in all cases.

      ZH final stateprecision
      Z\to e^+e^-H\to WW^*\to \ell\nu\ell'\nu,\, \ell\nu q\bar{q}2.6%
      Z\to \mu^+\mu^-H\to WW^*\to \ell\nu\ell'\nu,\, \ell\nu q\bar{q}2.4%
      Z\to \nu\bar{\nu}H\to WW^*\to \ell\nu q\bar{q}, q\bar{q} q\bar{q}1.5%
      Z\to q\bar{q}H\to WW^*\to q\bar{q} q\bar{q}1.7%
      combination0.9%

      Table 7.  Expected relative precision on the \sigma(ZH)\times{\rm BR}( H\to WW^*) measurement from a CEPC dataset of 5.6\,{\rm{ab}}^{-1}.

      For Z\to \ell^+\ell^-, the H\to WW^* decay final states studied are \ell\nu\ell'\nu and \ell\nu q\bar{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\to \nu\bar{\nu}, the \ell\nu q\bar{q} and q\bar{q} q\bar{q} final states are considered for the H\to WW^* 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 \ell\nu q\bar{q} (q\bar{q} q\bar{q}) final state. The total visible mass of the event must be in the range of 100–150 GeV for both \ell\nu q\bar{q} and q\bar{q} q\bar{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 Z\to q\bar{q}, the H\to WW^*\to q\bar{q} q\bar{q} decay is studied. Candidate events are reconstructed into 6 jets. Jets from Z\to q\bar{q}, W\to q\bar{q} and H\to WW^*\to q\bar{q} q\bar{q} decays are selected by minimizing the \chi^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\to \nu\bar{\nu} and H\to WW^*\to q\bar{q} q\bar{q} final state.

      Figure 12.  (color online) ZH production with Z\to \nu\bar{\nu} and H\to WW^*\to q\bar{q} q\bar{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.6\,{\rm{ab}}^{-1}, 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 \sigma(ZH)\times\mathrm{BR}( H\to WW^*) 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 H\to WW^* decay remain to be explored. Including these missing final states will no doubt improve the precision.

    • 5.3.   H\to ZZ^*

    • The H\to ZZ^* decay has a branching ratio 2.64% [33] for a 125 GeV Higgs boson in the SM. Events from e^+e^-\to ZH production with the H\to ZZ^* 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^+e^-\to ZH\to ZZZ^* process has a very rich variety of topologies.

      Studies are performed for a few selected ZH final states: Z\to \mu^+\mu^- and H\to ZZ^*\to \nu\bar{\nu} q\bar{q}; Z\to \nu\bar{\nu} and H\to ZZ^*\to \ell^+\ell^- q\bar{q}. The W and Z boson fusion processes, e^+e^-\to e^+e^- H and e^+e^-\to \nu\bar{\nu} H, are included in the Z( e^+e^-)H and Z( \nu\bar{\nu})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\to \mu^+\mu^- and H\to ZZ^*\to \nu\bar{\nu} q\bar{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 Z^*\to q\bar{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\to \mu^+\mu^- after the selection. The background is negligible in this final state.

      Figure 13.  (color online) ZH production with H\to ZZ^*: a) the recoil mass distribution of the \mu^+\mu^- system for Z\to \mu^+\mu^-, H\to ZZ^*\to \nu\bar{\nu} q\bar{q}; b) the invariant mass distribution of the \mu^+\mu^- q\bar{q} system for Z\to \nu\bar{\nu},\, H\to ZZ^*\to \mu^+\mu^- q\bar{q}. The markers and their uncertainties represent expectations from a CEPC dataset of 5.6\,{\rm{ab}}^{-1}, 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\to \nu\bar{\nu} and H\to ZZ^*\to\ell^+\ell^- q\bar{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 \sigma(ZH)\times{\rm BR}( H\to ZZ^*) 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 Z\to q\bar{q} and H\to ZZ^*\to q\bar{q} q\bar{q} which accounts for a third of all ZH\to ZZZ^* decay is not studied. Moreover, there are further potential improvements by using multivariate techniques.

      ZH final stateprecision
      Z\to \mu^+\mu^-H\to ZZ^*\to \nu\bar{\nu} q\bar{q}7.2%
      Z\to \nu\bar{\nu}H\to ZZ^*\to \ell^+\ell^- q\bar{q}7.9%
      combination4.9%

      Table 8.  Expected relative precision for the \sigma(ZH)\times \mathrm{BR}( H\to ZZ^*) measurement with an integrated luminosity 5.6\;{\rm{ab}}^{-1}.

    • 5.4.   H\to\gamma\gamma

    • 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\to\gamma\gamma and four different Z boson decay modes: Z\to \mu^+\mu^-,\, \tau^+\tau^-,\, \nu\bar{\nu} and q\bar{q}. The Z\to e^+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 Z\to q\bar{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\theta|<0.9. The photon pair with the highest invariant mass is retained as the H\to\gamma\gamma candidate and its recoil mass must be consistent with the Z boson mass. For the Z\to \mu^+\mu^- and Z\to \tau^+\tau^- decays, a minimal angle of 8^\circ 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^-\to (Z/\gamma^*)\gamma\gamma process where the \gamma'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\to \nu\bar{\nu} decay mode is shown in Fig. 14. A relative precision of 6.2% on \sigma(ZH)\times{\rm BR}(H\to\gamma\gamma) can be achieved.

      Figure 14.  (color online) ZH production with H\to \gamma\gamma: the diphoton invariant mass distribution for the Z\to \nu\bar{\nu} decay. The markers and their uncertainties represent expectations from a CEPC dataset of 5.6\,{\rm{ab}}^{-1}, whereas the solid blue curve is the fit result. The dashed curves are the signal and background components.

    • 5.5.   H\to Z\gamma

    • Similar to the H\to \gamma\gamma decay, the H\to Z\gamma decay in the SM is mediated by W-boson and top-quark loops and has a branching ratio of 0.154%. The H\to Z\gamma analysis targets the signal process of ZH\to ZZ\gamma\to \nu\bar{\nu} q\bar{q}\gamma, 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 \pm 12 GeV and \pm 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\bar{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, \Delta M = M_{ q\bar{q}\gamma}-M_{ q\bar{q}} and \Delta M = M_{ \nu\bar{\nu}\gamma}-M_{ \nu\bar{\nu}}, are calculated. Here the energy and momentum of the \nu\bar{\nu} 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 M_H-M_Z\sim 35 GeV whereas the other should be part of the continuum background. Figure 15 shows the \Delta M distribution expected from an integrated luminosity of 5.6~{\rm{ab}}^{-1}. 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 \sigma(ZH)\times {\rm BR}( H\to Z\gamma).

      Figure 15.  (color online) The distribution of the mass difference \Delta M (M_{ q\bar{q}\gamma}-M_{ q\bar{q}} and M_{ \nu\bar{\nu}\gamma}-M_{ \nu\bar{\nu}}) of the selected e^+e^-\to ZH\to ZZ\gamma\to \nu\bar{\nu} q\bar{q}\gamma candidates expected in a dataset with an integrated luminosity of 5.6\;{\rm{ab}}^{-1}. 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 ZH\to ZZ\gamma\to q\bar{q}\, q\bar{q}\gamma should further improve the precision on the \sigma(ZH)\times{\rm BR}( H\to Z\gamma) measurement.

    • 5.6.   H\to \tau^+\tau^-

    • The H\to \tau^+\tau^- decay has a branching ratio of 6.32% [33] at m_H = 125 GeV in the SM. The \tau-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^+e^-\to ZH production with H\to \tau^+\tau^- and Z\to \mu^+\mu^-, \nu\bar{\nu} and q\bar{q} decays. For Z\to \mu^+\mu^-, 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\to \nu\bar{\nu}, 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\to \mu^+\mu^- and contributes to approximately 40% of the total background for Z\to \nu\bar{\nu}. The rest of the background in the Z\to \nu\bar{\nu} channel comes from diboson production. For Z\to q\bar{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 ZH\to q\bar{q}\tau^+\tau^- 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 \sigma(ZH)\times {\rm BR}( H\to \tau^+\tau^-) expected from a CEPC dataset of 5.6\,{\rm{ab}}^{-1} for the three Z boson decay modes studied. The precision from the Z\to e^+e^- decay mode extrapolated from the Z\to \mu^+\mu^- study is also included. The e^+e^-\to e^+e^- H 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 \sigma(ZH)\times{\rm BR}( H\to \tau^+\tau^-).

      ZH final statePrecision
      Z\to \mu^+\mu^-H\to \tau^+\tau^-2.6%
      Z\to e^+e^-H\to \tau^+\tau^-2.7%
      Z\to \nu\bar{\nu}H\to \tau^+\tau^-2.5%
      Z\to q\bar{q}H\to \tau^+\tau^-0.9%
      combination0.8%

      Table 9.  Expected relative precision for the \sigma(ZH)\times \mathrm{BR}( H\to \tau^+\tau^-) measurement from a CEPC dataset of 5.6\,{\rm{ab}}^{-1}.

      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\to \mu^+\mu^- and (b) Z\to \nu\bar{\nu}. Here D_0 and Z_0 are the transverse and longitudinal impact parameters, respectively. The markers and their uncertainties represent expectations from a CEPC dataset of 5.6\,{\rm{ab}}^{-1}, 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\to \mu^+\mu^-

    • The dimuon decay of the Higgs boson, H\to \mu^+\mu^-, 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\times 10^{-4} [33] for m_{H} = 125 GeV.

      To estimate CEPC's sensitivity for the H\to \mu^+\mu^- decay, studies are performed for the ZH production with the Z decay modes: Z\to \ell^+\ell^- , Z\to \nu\bar{\nu} , and Z\to q\bar{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\to \ell^+\ell^- , candidate events must have another lepton pair with its mass consistent with m_Z . In the case of Z\to \mu^+\mu^- , the muon pairs of the Z\to \mu^+\mu^- and H\to \mu^+\mu^- decays are selected by minimizing a \chi^2 based on their mass differences with m_Z and m_H . For the Z\to \nu\bar{\nu} decay, a requirement on the missing energy is applied. For the Z\to q\bar{q} decay, candidate events must have two jets with their mass consistent with m_Z . 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\to \mu^+\mu^- muon pair, are exploited. Simple criteria on these variables are applied for the Z\to \ell^+\ell^- and Z\to \nu\bar{\nu} decay mode whereas a BDT is used for the Z\to q\bar{q} decay.

      In all analyses, the signal is extracted through unbinned likelihood fits to the M_{ \mu^+\mu^- } 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 \sigma(ZH)\times {\rm BR}( H\to \mu^+\mu^- ) measurement is estimated to be about 16% for data corresponding to an integrated luminosity of 5.6 \;{\rm ab}^{-1} .

      Figure 17.  (color online) ZH production with the H\to \mu^+\mu^- decay: dimuon invariant mass distribution of the selected H\to \mu^+\mu^- candidates expected from an integrated luminosity of 5.6 {\rm ab}^{-1} at the CEPC. The distribution combines contributions from Z\to \ell^+\ell^- , Z\to \nu\bar{\nu} , and Z\to q\bar{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: H\to {\rm inv}

    • In the SM, the Higgs boson can decay invisibly via H\to ZZ^* \to \nu\bar{\nu} \nu\bar{\nu} . For a Higgs boson mass of 125 GeV, this decay has a branching ratio of 1.06\times 10^{-3} . 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 {\rm BR}( H\to {\rm inv} ) measurement is studied for the Z\to \ell^+\ell^- and Z\to q\bar{q} decay modes. The H\to ZZ^* \to \nu\bar{\nu} \nu\bar{\nu} decay is used to model the H\to {\rm inv} 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\to \ell^+\ell^- 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 Z\to q\bar{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 \sigma(ZH)\times{\rm BR}( H\to {\rm inv} ) and the 95% confidence-level (CL) upper limit on {\rm BR}( H\to {\rm inv} ) from a CEPC dataset of 5.6\; {\rm ab}^{-1} . Subtracting the SM H\to ZZ^* \to \nu\bar{\nu} \nu\bar{\nu} contribution, a 95% CL upper limit of 0.30% on {\rm BR_{inv}^{BSM}} , the BSM contribution to the H\to {\rm inv} decay can be obtained.

      ZH final state studiedrelative precision on \sigma\times {\rm BR} upper limit on {\rm BR} ( H\to {\rm inv} )
      Z\to e^+e^- H\to {\rm inv} 339%0.82%
      Z\to \mu^+\mu^- H\to {\rm inv} 232%0.60%
      Z\to q\bar{q} H\to {\rm inv} 217%0.57%
      combination143%0.41%

      Table 10.  Expected relative precision on \sigma(ZH)\times {\rm BR}( H\to {\rm inv} ) and 95% CL upper limit on {\rm BR}( H\to {\rm inv} ) from a CEPC dataset of 5.6 \;{\rm ab}^{-1} .

    • 5.9.   Measurement of \sigma( e^+e^-\to \nu_e\bar{\nu}_e H )\times{\rm BR}( H\to b\bar{b} )

    • The W-fusion e^+e^-\to \nu_e\bar{\nu}_e H process has a cross section of 3.3% of that of the ZH process at \sqrt{s}=250 GeV. The product of its cross section and {\rm BR}( H\to b\bar{b} ) , \sigma( \nu\bar{\nu} H)\times\mathrm{BR}( H\to b\bar{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^- \to \nu\bar{\nu} H\to \nu\bar{\nu} b\bar{b} process has the same final state as the e^+e^- \to ZH\to \nu\bar{\nu} b\bar{b} process, but has a rate that is approximately one sixth of e^+e^- \to ZH\to \nu\bar{\nu} b\bar{b} at \sqrt{s}=250 GeV. The main non-Higgs boson background is the SM ZZ production.

      The Z( \nu\bar{\nu} )H background is irreducible and can also interfere with \nu\bar{\nu} H in the case of Z\to\nu_e\bar{\nu}_e . However, the interference effect is not considered in the current study. The \nu\bar{\nu} H and Z( \nu\bar{\nu} )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( \nu\bar{\nu} )H and show a continuum spectrum for \nu\bar{\nu} H . Furthermore, Higgs bosons are produced with different polar angular distributions, see Fig. 18(a).

      Figure 18.  (color online) Distributions of the b\bar{b} system of the e^+e^- \to \nu\bar{\nu} b\bar{b} events: (a) cosine of the polar angle \theta before the event selection and (b) the recoil mass after the event selection. Contributions from e^+e^-\to \nu_e\bar{\nu}_e H , ZH and other SM processes are shown. The \cos\theta 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 L_B described in Section 5.1. To separate \nu\bar{\nu} H and Z( \nu\bar{\nu} )H contributions, a 2-dimensional simultaneous fit in the plane of the recoil mass and polar angle of the b\bar{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\bar{b} system after the kinematic fit. A fit to the M_{ b\bar{b} }-\cos\theta distribution with both rates of \nu\bar{\nu} H and Z( \nu\bar{\nu} ) H processes as free parameters leads to relative precision of 2.9% for \sigma( \nu\bar{\nu} H)\times {\rm BR}( H\to b\bar{b} ) and 0.30% for \sigma(ZH)\times {\rm BR}( H\to b\bar{b} ) . The latter is consistent with the study of the H \to b\bar b/c\bar c/gg decay described in Section 5.1. Fixing the Z( \nu\bar{\nu} )H( b\bar{b} ) contribution to its SM expectation yields a relative precision of 2.6% on \sigma( e^+e^-\to \nu_e\bar{\nu}_e H )\times {\rm BR}( H\to b\bar{b} ) .

    6.   Combinations of individual measurements

      6.1.   Combined measurements of \sigma\times\mathrm{BR} and BR

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

      L({\rm{BR}},\theta ) = {\rm{Poisson}}\left. {[{N^{{\rm{obs}}}}|{N^{{\rm{exp}}}}({\rm{BR}},\theta )} \right] \cdot G(\theta ),

      (3)

      where \mathrm{BR} is the parameter of interest and \theta represents nuisance parameters associated with systematic uncertainties. The number of observed events is denoted by N^{\rm obs} , N^{\rm exp}(\mathrm{BR},\theta) is the expected number of events, and G(\theta) 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^+e^-\to ZH cross section \sigma(ZH) measured from the recoil method, Higgs boson branching ratio \mathrm{BR} , the event selection efficiency \epsilon . The number of the expected background events, N^b , is estimated using Monte Carlo samples. Thus:

      {N^{{\rm{exp}}}}({\rm{BR}},\theta ) = {\rm{Lumi}}({\theta ^{{\rm{lumi}}}}) \times {\sigma _{ZH}}({\theta ^\sigma }) \times {\rm{BR}} \times \epsilon ({\theta ^\epsilon }) + {N^b}({\theta ^b}),

      (4)

      where \theta^X\ (X={\rm lumi}, \sigma , \epsilon and b ) are the nuisance parameters of their corresponding parameters or measurements. Even with 10^6 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 H \to b\bar b/c\bar c/gg analysis suffers from contamination from the H\to WW^*,ZZ^*\to q\bar{q} q\bar{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 H\to WW^* and H\to ZZ^* 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 \sigma(ZH)\times {\rm BR}(H\to ZZ^*) 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\bar{b} , c\bar{c} , gg , WW^* , ZZ^* and \tau^+\tau^- , percent level precision is expected. The best achievable statistical uncertainties for a dataset of 5.6~{\rm ab}^{-1} are 0.26% for \sigma(e^+e^-\to ZH)\times \mathrm{BR}( H\to b\bar{b} ) and 0.5% for \sigma(e^+e^-\to ZH) . 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
      m_H 5.9 MeV5.9 MeV
      \Gamma_H 2.7%2.8%
      \sigma(ZH) 0.5%0.5%
      \sigma( \nu\bar{\nu} H) 3.0%3.2%
      decay mode \sigma\times\mathrm{BR} {\rm BR} \sigma\times\mathrm{BR} {\rm BR}
      H\to b\bar{b} 0.26%0.56%0.27%0.56%
      H\to c\bar{c} 3.1%3.1%3.3%3.3%
      H\to gg 1.2%1.3%1.3%1.4%
      H\to WW^* 0.9%1.1%1.0%1.1%
      H\to ZZ^* 4.9%5.0%5.1%5.1%
      H\to \gamma\gamma 6.2%6.2%6.8%6.9%
      H\to Z\gamma 13%13%16%16%
      H\to \tau^+\tau^- 0.8%0.9%0.8%1.0%
      H\to \mu^+\mu^- 16%16%17%17%
      {\rm BR_{inv}^{BSM}} - <0.28% - <0.30%

      Table 11.  Estimated precision of Higgs boson property measurements for the CEPC-v1 detector concept operating at \sqrt{s}=250 GeV. All precision are relative except for m_H and {\rm BR_{inv}^{BSM}} for which \Delta m_H and 95% CL upper limit are quoted respectively. The extrapolated precision for the CEPC-v4 concept operating at \sqrt{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 ZH\to ZWW^* and ZH\to ZZZ^* 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 \sqrt{s}=250 GeV. However, given the relative small differences in the performance of the two detector concepts and in \sqrt{s} , the results for CEPC-v4 operating at \sqrt{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^+e^-\to ZH and e^+e^-\to \nu_e\bar{\nu}_e H 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\to \mu^+\mu^- 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 \sqrt{s}=240 GeV are used.

    • 6.3.   Measurement of Higgs boson width

    • The Higgs boson width ( \Gamma_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^+e^-\to ZH cross section \sigma(ZH) can be measured from the recoil mass distribution, independent of the Higgs boson decays.

      Measurements of \sigma(ZH) and \mathrm{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:

      {\Gamma _H} = \frac{{\Gamma (H \to Z{Z^*})}}{{{\rm{BR}}(H \to Z{Z^*})}} \propto \frac{{\sigma (ZH)}}{{{\rm{BR}}(H \to Z{Z^*})}},

      (5)

      where \Gamma( H\to ZZ^* ) is the partial width of the H\to ZZ^* decay. Because of the small expected {\rm BR}( H\to ZZ^* ) value for a 125 GeV Higgs boson (2.64% in the SM), the precision of \Gamma_H is limited by the H\to ZZ^* analysis statistics. It can be improved including the decay final states with larger branching ratios, e.g. the H\to b\bar{b} decay:

      {\Gamma _H} = \frac{{\Gamma (H \to b\bar b)}}{{{\rm{BR}}(H \to b\bar b)}},

      (6)

      where the partial width \Gamma( H\to b\bar{b} ) can be independently extracted from the cross section of the W fusion process e^+e^- \to \nu\bar{\nu} H\to \nu\bar{\nu} b\bar{b} :

      \sigma (\nu \bar \nu H \to \nu \bar \nu b\bar b) \propto \Gamma (H \to W{W^*}) \cdot {\rm{BR}}(H \to b\bar b)

      (7)

      = \Gamma (H \to b\bar b) \cdot {\rm{BR}}(H \to W{W^*}).

      (8)

      Thus, the Higgs boson total width is:

      {\Gamma _H} = \frac{{\Gamma (H \to b\bar b)}}{{{\rm{BR}}(H \to b\bar b)}} \propto \frac{{\sigma ({e^ + }{e^ - } \to {\nu _e}{{\bar \nu }_e}H)}}{{{\rm{BR}}(H \to W{W^*})}},

      (9)

      where {\rm BR}( H\to b\bar{b} ) and {\rm BR}( H\to WW^* ) are measured from the e^+e^-\to ZH process. The limitation of this method is the precision of the \sigma(e^+e^-\to \nu\bar{\nu} H \to \nu\bar{\nu} b\bar{b} ) measurement.

      The expected precision on \Gamma_H is 5.1% from the measurements of \sigma(ZH) and {\rm BR}( H\to ZZ^* ) and is 3.5% from the measurements of \sigma( \nu\bar{\nu} H\to \nu\bar{\nu} b\bar{b} ) , {\rm BR}( H\to b\bar{b} ) and {\rm BR}( H\to WW^* ) . The quoted precision is dominated by the {\rm BR}( H\to ZZ^* ) measurement for the former case and the \sigma( \nu\bar{\nu} H\to \nu\bar{\nu} b\bar{b} ) measurement for the latter case. The combined \Gamma_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 \kappa -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 \kappa -framework

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

      {\kappa _f} = \frac{{g(Hff)}}{{{g_{{\rm{SM}}}}(Hff)}},\;\;{\kappa _V} = \frac{{g(HVV)}}{{{g_{{\rm{SM}}}}(HVV)}},

      (10)

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

      \begin{array}{l} \sigma (ZH) = \kappa _Z^2 \cdot {\sigma _{{\rm{SM}}}}(ZH)\\ \sigma (ZH) \times {\rm{BR}}(H \to ff) = \dfrac{{\kappa _Z^2\kappa _f^2}}{{\kappa _\Gamma ^2}} \cdot {\sigma _{{\rm{SM}}}}(ZH) \times {\rm{B}}{{\rm{R}}_{{\rm{SM}}}}(H \to ff) \end{array}

      (11)

      Here \kappa_\Gamma^2(\equiv \Gamma_H/\Gamma_H^{\rm SM}) 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\gamma\gamma and HZ\gamma in the SM. In the absence of new physics, these couplings, often referred to as the effective couplings, can be expressed using the \kappa parameters, described previously. However, new physics states in the loops can alter these couplings. For this reason, three additional \kappa parameters: \kappa_g , \kappa_\gamma and \kappa_{Z\gamma} are introduced to parametrize the potential deviations from the SM for the three effective Hgg, H\gamma\gamma and HZ\gamma 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 {\rm BR_{inv}^{BSM}} .

      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 \Gamma_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 ( > \mathcal{O}(1) ) is unlikely. For smaller deviations, the Higgs phenomenology is not sensitive to the deviations of \kappa_e , \kappa_u , \kappa_d and \kappa_s as the Higgs boson couplings to these particles are negligible compared with the couplings to other particles [78]. Therefore, these \kappa 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\gamma\gamma and HZ\gamma amplitudes. However, this effect is parametrized by \kappa_g , \kappa_\gamma and \kappa_{Z\gamma} already. Therefore, \kappa_t is not considered as an independent parameter. For simplicity, previous studies often do not include \kappa_{Z\gamma} in the fit. We will follow this approach here. This leaves the following set of 10 independent parameters:

      {\kappa _b},\;{\kappa _c},\;{\kappa _g},\;{\kappa _W},\;{\kappa _\tau },\;{\kappa _Z},\;{\kappa _\gamma },\;{\kappa _\mu },\;{\rm{BR}}_{{\rm{inv}}}^{{\rm{BSM}}},\;{\Gamma _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 H\to ZZ^* \to \nu\bar{\nu} \nu\bar{\nu} ) [47, 79]:

      {\kappa _b},\;{\kappa _c},\;{\kappa _g},\;{\kappa _W},\;{\kappa _Z},\;{\kappa _\gamma },\;{\kappa _\tau } = {\kappa _\mu }.

      (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 \kappa_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 \kappa -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 \kappa modifications violate gauge invariance. Higher order corrections in the \kappa -framework cannot be easily defined. Moreover, the \kappa_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 \kappa_Z and \kappa_W , and further expands them into three different Lorentz structures. Moreover, some of these higher dimensional HVV couplings are also connected with \kappa_\gamma and anomalous trilinear gauge couplings. The current EFT analysis does not include any new light degrees of freedom, in contrast to the \kappa -framework with independent parameters {\rm BR_{inv}^{BSM}} and \Gamma_H . Overall, \kappa -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\bar{t} H 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 \kappa -framework, the LHC, with its large statistics, improves the precision of rare processes such as H\to \gamma\gamma . 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 pp \to H\to \gamma\gamma and pp \to H\to ZZ^* at the LHC and the measurement of the HZZ coupling at the CEPC can significantly improve the \kappa_\gamma 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.6 \;{\rm ab}^{-1} 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 \sqrt{s}=14 TeV and accumulate an integrated luminosity of 3000 {\rm fb}^{-1} . 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 {\rm BR_{inv}^{BSM}}
      Quantity10-parameter fit7-parameter fit
      CEPCCEPC+HL-LHCCEPCCEPC+HL-LHC
      \kappa_b 1.3%1.0%1.2%0.9%
      \kappa_c 2.2%1.9%2.1%1.9%
      \kappa_g 1.5%1.2%1.5%1.1%
      \kappa_W 1.4%1.1%1.3%1.0%
      \kappa_\tau 1.5%1.2%1.3%1.1%
      \kappa_Z 0.25%0.25%0.13%0.12%
      \kappa_\gamma 3.7%1.6%3.7%1.6%
      \kappa_\mu 8.7%5.0%
      {\rm BR_{inv}^{BSM}} <0.30%<0.30%
      \Gamma_H 2.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 {\rm BR_{inv}^{BSM}} 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 \kappa_Z is more than a factor of 10 better. The CEPC can also improve significantly the precision on a set of \kappa parameters that are affected by large backgrounds at the LHC, such as \kappa_b , \kappa_c , and \kappa_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 \kappa_\gamma , 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.6 \;{\rm ab}^{-1} 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 {\rm fb}^{-1} 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~{\rm ab}^{-1} 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 \kappa_Z . It can also set a much stronger constraint on {\rm BR_{inv}^{BSM}} .

      Figure 20.  (color online) The 10 parameter fit results for the CEPC at 240 GeV with an integrated luminosity of 5.6 {\rm ab}^{-1} (light red bars) and for the combination with the HL-LHC inputs (dark red bars). All the numbers are relative precision except for {\rm BR_{inv}^{BSM}} 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:

      {{\cal L}_{{\rm{EFT}}}} = {{\cal L}_{{\rm{SM}}}} + \sum\limits_i {\frac{{c_i^{(6)}}}{{{\Lambda ^2}}}} {\cal O}_i^{(6)} + \sum\limits_j {\frac{{c_j^{(8)}}}{{{\Lambda ^4}}}} {\cal O}_j^{(8)} + \cdots

      (14)

      where \Lambda 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, c_i , 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 \kappa -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 ( \sim 10^{11} 10^{12} ) 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\bar{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^+e^- \to WW 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 \delta g_{1, Z} , \delta \kappa_\gamma and \lambda_Z [100, 101]. Among them, \delta g_{1,Z} and \delta \kappa_\gamma are generated by operators that also contribute to the Higgs boson processes. At 240 GeV, the e^+e^- \to WW process cross section is almost two orders of magnitude larger than that of the e^+e^-\to ZH 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 \chi^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.6 {\rm ab}^{-1} )
      uncertaintycorrelation matrix
      \delta g_{1,Z} \delta \kappa_\gamma \lambda_Z
      \delta g_{1,Z} 1.2\times10^{-3} 10.08−0.90
      \delta \kappa_\gamma 0.9\times10^{-3} 1−0.42
      \lambda_Z 1.3\times10^{-3} 1

      Table 13.  The estimated constraints on aTGCs from the measurements of the diboson process ( e^+e^- \to WW ) in the semi-leptonic channel at the CEPC 240 GeV with 5.6 {\rm ab}^{-1} 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, \mathcal{O}_{3W} corresponds to the aTGC parameter \lambda_Z , \mathcal{O}_{HW} and \mathcal{O}_{HB} generate the aTGC parameters \delta g_{1,Z} and \delta \kappa_\gamma 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\gamma\gamma and HZ\gamma vertices normalized to the SM one-loop contributions, and denoted as \bar{c}_{gg} , \bar{c}_{\gamma\gamma} and \bar{c}_{Z\gamma} (as opposed to c_{gg} , c_{\gamma\gamma} and c_{Z\gamma} in Ref. [103]). The parameter \bar{c}^{\rm eff}_{gg} is further defined to absorb all contributions to the Hgg vertex. With these redefinitions, the set of twelve parameters is given by

      \mathcal{O}_H = \frac{1}{2} (\partial_\mu |H^2| )^2 \mathcal{O}_{GG} = g_s^2 |H|^2 G^A_{\mu\nu} G^{A,\mu\nu}
      \mathcal{O}_{WW} = g^2 |H|^2 W^a_{\mu\nu} W^{a,\mu\nu} \mathcal{O}_{y_u} = y_u |H|^2 \bar{Q}_L \tilde{H} u_R + {\rm h.c.} (u \to t, c)
      \mathcal{O}_{BB} = g'^2 |H|^2 B_{\mu\nu} B^{\mu\nu} \mathcal{O}_{y_d} = y_d |H|^2 \bar{Q}_L H d_R + {\rm h.c.} (d \to b)
      \mathcal{O}_{HW} = ig(D^\mu H)^\dagger \sigma^a (D^\nu H) W^a_{\mu\nu} \mathcal{O}_{y_e} = y_e |H|^2 \bar{L}_L H e_R + {\rm h.c.} (e \to \tau, \mu)
      \mathcal{O}_{HB} = ig'(D^\mu H)^\dagger (D^\nu H) B_{\mu\nu} \mathcal{O}_{3W} = \frac{1}{3!} g \epsilon_{abc} W^{a \nu}_\mu W^b_{\nu \rho} W^{c \rho\mu}

      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. G^A_{\mu\nu} , W^a_{\mu\nu} and B_{\mu\nu} are the field strength tensors for the SM SU(3)_c , SU(2)_L and U(1)_Y gauge fields, respectively. For \mathcal{O}_{y_u} , \mathcal{O}_{y_d} and \mathcal{O}_{y_e} , only the contributions to the diagonal elements of the Yukawa matrices that corresponds to the top, charm, bottom, tau, and muon couplings are considered.

      \delta {c_Z},{c_{ZZ}},{c_{Z\square}},{\bar c_{\gamma \gamma }},{\bar c_{Z\gamma }},\bar c_{gg}^{{\rm{eff}}},\delta {y_t},\delta {y_c},\delta {y_b},\delta {y_\tau },\delta {y_\mu },{\lambda _Z}.

      (15)

      These parameters can be conveniently interpreted as the precision of the Higgs boson couplings analogous to those in the \kappa -framework. In particular, \delta c_Z , \bar{c}_{\gamma\gamma} , \bar{c}_{Z\gamma} , \bar{c}^{\rm eff}_{gg} and \delta y_{t, c, b, \tau, \mu} modifies the sizes of the SM Higgs boson couplings to ZZ, \gamma\gamma , Z\gamma , gg and fermions, respectively. c_{ZZ} and c_{Z\square} parametrize the anomalous HZZ couplings:

      {\cal L} = \frac{H}{v}[{c_{ZZ}}\frac{{{g^2} + {{g'}^2}}}{4}{Z_{\mu \nu }}{Z^{\mu \nu }} + {c_{Z\square}}{g^2}{Z_\mu }{\partial _\nu }{Z^{\mu \nu }}] + ...,

      (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, \lambda_Z is kept in Eq. 15, while \delta g_{1,Z} and \delta \kappa_\gamma are written in terms of the linear combinations of c_{ZZ} , c_{Z\square} , \bar{c}_{\gamma\gamma} and \bar{c}_{Z\gamma} . 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^+e^-\to ZH , Z\to \ell^+\ell^-, H\to b\bar{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 \chi^2 is obtained by summing over the \chi^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 \chi^2 can be written as

      {\chi ^2} = \sum\limits_{ij} {{{(c - {c_0})}_i}} \sigma _{ij}^{ - 2}{(c - {c_0})_j},{\rm{where}}\;\sigma _{ij}^{ - 2} \equiv {\left( {\delta {c_i}{\rho _{ij}}\delta {c_j}} \right)^{ - 1}},

      (17)

      where c_i 's are the EFT parameters, c_0 's are the corresponding central values which are zero by construction, as the measurements are assumed to be SM-like. The one-sigma uncertainties \delta c_i and the correlation matrix \rho can be obtained from \sigma^{-2}_{ij} = {\partial^2 \chi^2} \big/ {\partial c_i \partial c_j} .

      For comparison, the sensitivities of the LHC 14 TeV with total luminosities of 300~{\rm fb}^{-1} and 3000~{\rm fb}^{-1} are also considered. These are combined with the diboson ( e^+e^- \to WW ) 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~{\rm fb}^{-1} and 3000~{\rm fb}^{-1} ) combined with the LEP diboson measurements are shown in comparison. The results of the combination of the CEPC with the HL-LHC ( 3000\; {\rm fb}^{-1} ) 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 300 \;{\rm fb}^{-1} ( 3000 \;{\rm fb}^{-1} ) 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 300\; {\rm fb}^{-1} (light gray bar) and 3000 \;{\rm fb}^{-1} (dark gray bar) combined with the LEP diboson ( e^+e^- \to WW ) measurement. The second column shows the results from the CEPC with 5.6\; {\rm ab}^{-1} 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, \delta y_c is fixed to zero.

      Figure 22.  (color online) The 95%, CL sensitivity to \Lambda/\sqrt{|c_i|} for the operators in the basis defined in Table 14. The first two columns show the results from the LHC Higgs boson measurements with 300\;{\rm fb}^{-1} and 3000\;{\rm fb}^{-1} combined with the LEP diboson (e^+e^- \to WW) measurement. The last two columns show the results from the CEPC alone and the combination of the CEPC and the HL-LHC (3000\;{\rm fb}^{-1}). 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, \delta y_c is fixed to zero.

      In Fig. 22, the results are presented in terms of the sensitivity to \Lambda/\sqrt{|c_i|} at 95% CL for each operator as defined in Eq. (14), where \Lambda is the scale of new physics and c_i is the corresponding Wilson coefficient. Four columns are shown separately for the LHC 300\; {\rm fb}^{-1} , the HL-LHC 3000\; {\rm fb}^{-1} , 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 \bar{c}_{\gamma\gamma} , \bar{c}_{Z\gamma} and \delta y_\mu , the clean signal and small branching ratios of the corresponding channels ( H\to \gamma\gamma/Z\gamma/\mu\mu ) 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 \delta y_t 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 \bar{c}^{\rm eff}_{gg} absorbs both contributions and reflects the overall precision of the Hgg coupling. The combination with the LHC t\bar{t}H measurements can resolve this flat direction. The CEPC measurements, in turn, can improve the constraint on \delta y_t set by the LHC by providing much better constraints on the other parameters that contribute to the t\bar{t}H 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, \delta y_c is fixed to be zero for the LHC-only fits, as treating \delta y_c 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 \delta y_c to about 2%.

      Regarding the sensitivity to \Lambda/\sqrt{|c_i|} 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 ( c_i \sim 1 ), 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., O_{WW} and O_{BB} from the measurement of H\to \gamma\gamma ), 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 O_{GG} and O_{y_t} 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\bar{t}H 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 \Lambda 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.6\; {\rm ab}^{-1} at 240 GeV with unpolarized beams, both without and with the combination with the HL-LHC ( 3000\; {\rm fb}^{-1} ) Higgs boson measurements. With both the one-sigma bounds and the correlation matrix, the corresponding \chi^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
      \delta c_Z c_{ZZ} c_{Z\square} \bar{c}_{\gamma\gamma} \bar{c}_{Z\gamma} \bar{c}^{\rm eff}_{gg} \delta y_t \delta y_c \delta y_b \delta y_\tau \delta y_\mu \lambda_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
      c_i/\Lambda^2 [{\rm TeV}^{-2}] of dimension-six operators
      c_H c_{WW} c_{BB} c_{HW} c_{HB} c_{GG} c_{y_t} c_{y_c} c_{y_b} c_{y_\tau} c_{y_\mu} c_{3W}
      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.6 {\rm ab}^{-1} ) 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\bar{t}H measurements, \delta y_t can not be constrained in a global fit, thus c_{GG} and c_{y_t} 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 \mathcal{O}(10) level, while the HL-LHC is expected to improve the precision to the level of \mathcal{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 ( \gtrsim 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^+e^-\to ZH 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 \delta \kappa_\lambda \equiv \kappa_\lambda -1 , where \kappa_\lambda is the ratio of the Higgs boson trilinear coupling to its SM value,

      {\kappa _\lambda } \equiv \frac{{{\lambda _3}}}{{\lambda _3^{{\rm{sm}}}}},\quad \lambda _3^{{\rm{sm}}} = \frac{{m_H^2}}{{2{v^2}}}.

      (18)

      The global fit is performed simultaneously with \delta \kappa_\lambda 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 \delta \kappa_\lambda is around \pm 3 , significantly worse than the 35% in the \delta \kappa_\lambda -only fit. This is a clear indication that it is difficult to resolve the effects of \delta \kappa_\lambda from other Higgs boson couplings. For the HL-LHC, the reach on \delta \kappa_\lambda 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 \delta \kappa_\lambda 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 \chi^2 as a function of \delta \kappa_\lambda for the two colliders.

      bounds on \delta\kappa_\lambda \Delta\chi^2=1 \Delta\chi^2=4
      CEPC 240 GeV ( 5.6~{\rm ab}^{-1} )[−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 \Delta\chi^2=1 (one-sigma) and \Delta\chi^2=4 (two-sigma) bounds of \delta\kappa_\lambda 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 \delta\kappa_\lambda 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 \mathcal{O}_{y_t} and \mathcal{O}_{y_b} ,

      \Delta {y_t} = y_t^{{\rm{SM}}}\left( {\Re [{C_{{y_t}}}]\frac{{{v^3}}}{{2{m_t}{\Lambda ^2}}} + i\Im [{C_{{y_t}}}]\frac{{{v^3}}}{{2{m_t}{\Lambda ^2}}}} \right),

      (19)

      \Delta {y_b}{\rm{ }} = y_t^{{\rm{SM}}}\left( {\Re [{C_{{y_b}}}]\frac{{{v^3}}}{{2{m_b}{\Lambda ^2}}} + i\Im [{C_{{y_b}}}]\frac{{{v^3}}}{{2{m_b}{\Lambda ^2}}}} \right).

      (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\to \gamma\gamma and H\to gg for \mathcal{O}_{y_t} , and H\to gg and H\to b\bar b for \mathcal{O}_{y_b} . Given that H\to gg measurements are sensitive to both operators, a joint analysis of \mathcal{O}_{y_t} and \mathcal{O}_{y_b} 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 {\cal O}_{tH} and {\cal O}_{bH} 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\bar tH measurements, red bands are constraints from the CEPC H\to gg measurements and blue bands are constraints from the CEPC H\to \gamma\gamma measurements. Similarly, in the right panel, the cyan bands are constraints from H\to b\bar b and the red bands are constraints from H\to gg at the CEPC.

      Figure 24.  (color online) Results for analysis on C_{y_t} and C_{y_b} 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 \mathcal{O}_{y_t} and \mathcal{O}_{y_b} 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 \pm 0.16\pi . This constraint is driven by the H\to \gamma\gamma measurement, where a sizable phase shift will enlarge the H\to\gamma\gamma 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 H\to gg measurement which is dominated by the top-quark loop contribution. Note that constraints from the H\to gg 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 \pm 0.47\pi 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^-(\to Z^*)\to ZH\to \mu^+\mu^- b\bar{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, f_{a2} , and the fraction of a CP-odd contribution due to new physics, f_{a3} . The following two types of observables can be used to measure these anomalous couplings of the Higgs bosons.

      1) The dependence of the e^+e^- \to Z^*\to ZH cross section on \sqrt{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 \sqrt{s} is assumed for the CEPC operating as a Higgs boson factory.

      2) Angular distributions, \cos\theta_1 or \cos\theta_2 and \Phi 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\theta_1 or \cos\theta_2 and non-trivial phase in the \Phi distributions can lead to an unambiguous interpretation of CP violation.

      Figure 25.  (color online) The Higgs boson production and decay angles for the e^+e^- \to Z^* \to ZH \to \mu^+\mu^- b\bar{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^+e^-\to ZZ\to \ell^+\ell^- b\bar{b} process in MadGraph [134].

      Several thousand statistically independent pseudo-experiments are generated and fitted to estimate the sensitivity to f_{a2} and f_{a3} , defined as the smallest values that can be measured with 3\sigma 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\sigma discovery is estimated to be 0.018 for f_{a2} and 0.007 for f_{a3} . Figure 26(a,b) show the distributions of the fitted values of f_{a2} and f_{a3} from the pseudo-experiments expected for f_{a2}=0.018 and f_{a3}=0.008 , respectively. A simultaneous fit of f_{a2} and f_{a3} is also performed with the 68% and 95% CL contours shown in Fig. 26(c).

      Figure 26.  (color online) Distribution of fitted values of f_{a2} and f_{a3} 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 f_{a2} and f_{a3} , with 68% and 95% confidence level contours shown.

      The sensitivities for f_{a2} and f_{a3} are then converted to the corresponding parameters defined for the on-shell H\to ZZ^* decays, f_{a2}^{\rm dec} and f_{a3}^{\rm dec} , in order to compare with the sensitivities from the LHC experiments as described in Ref. [133]. The corresponding sensitivities of f_{a2}^{\rm dec} and f_{a3}^{\rm dec} are 2\times 10^{-4} and 1.3\times 10^{-4} , respectively. The much smaller values in the f_{a2,a3}^{\rm dec} are due to the much larger m_{Z^*}^2 in the e^+e^- \to Z^*\to ZH 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 f_{a2} and f_{a3} at the CEPC are better by a factor of 300 and 3. Further improvements can be achieved by exploring kinematics in the H\to b\bar{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 \sqrt{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, H\to X_1 X_2 , where the BSM particles X_i 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^+e^-\to ZH 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.6 {\rm ab}^{-1} 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
      {\not\!\! E_{\rm{T}}} 0.230.0560.0030
      (b\bar b)+{\not \!\!E_{\rm{T}}} [0.2]1×10−4
      (jj)+{\not\!\! E_{\rm{T}}} 4×10−4
      (\tau^+\tau^-)+{\not\!\! E_{\rm{T}}} [1]8×10−5
      b\bar b+{\not\!\! E_{\rm{T}}} [0.2]2×10−4
      jj+{\not\!\! E_{\rm{T}}} 5×10−4
      \tau^+\tau^-+{\not\!\! E_{\rm{T}}} 8×10−5
      (b\bar b)(b\bar b) 1.7 (0.2)6×10−4
      (c\bar c)(c\bar c) (0.2)8×10−4
      (jj)(jj) [0.1]2×10−3
      (b\bar b)(\tau^+\tau^-) 0.1 [0.15]4×10−4
      (\tau^+\tau^-)(\tau^+\tau^-) 1.2 [0.2~0.4]2×10−4
      (jj)(\gamma\gamma) [0.01]1×10−4
      (\gamma\gamma)(\gamma\gamma) 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.6 {\rm ab}^{-1} , based on the results from Ref. [136]. The second column shows the current LHC results and the projections for 100 {\rm fb}^{-1} (in parentheses) and 300 {\rm fb}^{-1} (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\bar{b} )+{\not \!\!E_{\rm{T}}} stands for X+{\not\!\! E_{\rm{T}}}\to ( b\bar{b} )+{\not\!\! E_{\rm{T}}} 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\bar b+{\not E_{\rm{T}}} , jj+{\not E_{\rm{T}}} and \tau^+\tau^-+{\not E_{\rm{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) = - \frac{1}{2}{\mu ^2}|H{|^2} + \frac{\lambda }{4}|H{|^4},

      (21)

      with the vacuum expectation value (VEV) of the Higgs field spontaneously breaking the SU(2)_{\rm L} \times U(1)_{\rm Y} gauge symmetry down to U(1)_{\rm 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), \mu^2 and \lambda , 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 \Lambda_{\rm weak} \sim 10^2 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, M_{\rm Planck} = 10^{19} 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

      \Lambda _{{\rm{weak}}}^2 \sim {c_1}M_1^2 + {c_2}M_2^2 + ...,

      (22)

      where M_i \sim M_{\rm NP} are the scale of new physics. They are typically the masses of the new physics particles. The c_i 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 M_{\rm NP} \gg \Lambda_{\rm EW} , with the level of tuning proportional to M_{\rm NP}^2 . 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

      \delta {\kappa _W},\;\delta {\kappa _Z} \sim {\cal O}\left( {\frac{{{v^2}}}{{{f^2}}}} \right).

      (23)

      Therefore, the measurement of \kappa_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 \kappa_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_\rho \sim g_\rho f , where g_\rho 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_\rho , coupling parameter g_{\rho L} and mixing parameter \xi \equiv v^2/f^2.

      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, \tilde{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 \gamma \gamma couplings since they are also generated at the one-loop level in the SM. The dominant effect is on the Hgg coupling,

      {\kappa _g} - 1 \simeq \frac{{m_t^2}}{{4m_{\tilde t}^2}}.

      (24)

      The measurement of \kappa_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 \phi_t with its only interaction to the SM fields given by H^\dagger H \phi_t^\dagger \phi_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 \kappa_Z and provides the best constraint on the mass of the top-quark partner, m_\phi . As we can see from the left panel of Fig. 31, the CEPC will be able to probe m_\phi 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 \sigma_{ZH} at the Higgs factory in the scalar singlet top-quark partner model with the H^\dagger H \phi^\dagger_t \phi_t interaction, reproduced from Ref. [144]. (b) Projected constraints in the folded stop mass plane from the H\gamma\gamma 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 g_{111} 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 v^2/M_{\rm NP}^2 . 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

      \begin{split} V(H,S) =& \displaystyle\frac{1}{2}{\mu ^2}|H{|^2} + \frac{\lambda }{4}|H{|^4} + m_S^2{S^2} \\ &+\tilde aS|H{|^2} + \tilde \kappa {S^2}|H{|^2} + \tilde b{S^3} + {\lambda _S}{S^4}. \end{split}

      (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 v^2/m_S^2 . At the same time, it will also generate the operator |H^\dagger \partial H|^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 \sim v^2 / m_{S}^2 . We thus expect \kappa_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 \sigma(ZH) at CEPC, the sensitivities of the current and future electroweak precision tests are also presented [152]. The \sigma(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 \sigma(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, |H^\dagger \partial H|^2 .

      A more restricted scenario, in which a discrete Z_2 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 \sigma(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 \gamma \gamma couplings. One such example is shown in Fig. 37(c), with a 6% deviation in the H\gamma\gamma 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 \kappa_\gamma 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, \sigma(ZH) and H\gamma\gamma couplings in models with Z_2 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 \xi = v(T_c)/T_c , where v(T_c) is the Higgs VEV at temperature T_c . 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, H^\dagger H {\mathcal{O}}_{\rm NP} , where {\mathcal{O}}_{\rm NP} is an operator composed out of new physics fields. Since H^\dagger H 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 m_H /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