Sensitivity study of anomalous HZZ couplings at a future Higgs factory

We study the sensitivity of constraining the model independent HZZ coupling based on the effective theory up to dimension-6 operators at a future Higgs factory. Using the current conceptual design parameters of the Circular Electron Positron Collider, we give the experimental limits for the model independent operators given by the total Higgsstrahlung cross-section and the angular distribution of Z boson decays. In particular, we give the very small sensitivity limit for the CP violation parameter , which will be a clear window to test the Standard Model and look for new physics signals.


Introduction
Since the Large Hadron Collider (LHC) announced the discovery of a scalar-like resonance [1,2], many subsequent measurements confirmed that it is the Higgs particle, the last brick of the Standard Model (SM) [3]. Among these measurements, the generic Higgs coupling to vector gauge bosons presented the largest (~ 7%) deviation from the SM prediction, but its experimental error is so big (±56%) that the good agreement with SM still stands. Both the unexplained phenomena, like the dark matter, and the theoretical issues like the hierarchy (naturalness) problem, still need extensions of SM. The interaction between the Higgs scalar and vector gauge bosons is a key ingredient for the underlying nature of spontaneous breaking of electroweak gauge symmetry. In addition to the suggested experiments of kinematic distributions at LHC Run2 [4,5] and high-luminosity run [6], future Higgs factories are under considerations, such as the International Linear Collider (ILC) in Japan [7], the Circular Electron Positron Collider (CEPC) in China [8,9], and Future Circular Collider (FCC-ee) in Europe [10]. Along these lines, many pre-analyses have been made to unveil the nature of the gauge boson and Higgs couplings [11,12]. Future lepton colliders will accumulate events with full kinematics and less background at high luminosity, which will enable precision tests in the Higgs sector. Hence, it becomes indispensable to unfold and use the event shapes in the experimental data to explore the subtleties of the Higgs properties. Theoretically, the effective field theory (EFT) provides a natural framework. A comprehensive study usually lists the relevant higher dimension operators and constrains their Wilson coefficients separately [13][14][15]. Since the EFT operators are defined by the interaction bases before gauge symmetry breaking, all effective HVV couplings or become unavoidable after symmetry breaking when compared to experiments at the electroweak scale. An alternative approach is to start from phenomenological parameters defined by mass eigenstates as a convenient bridge between EFT and experiments, but at a price of losing the direct map to the underlying operators [16][17][18][19][20]. Many works have been done theoretically for the LHC experiments to demonstrate how these phenomenological parameters can be constrained using the distribution of the polarization angle [21][22][23].
In lepton colliders, the Higgs will be produced via strahlung [24,25] in association with the Z gauge boson (mainly via off-shell Z), so that at a given luminosity the

HWW, Hγγ
HZγ events via HZZ vertex will undoubtedly play a more important role due to larger statistics and lower experimental uncertainties than or . In order to connect the HZZ event shape to the phenomenological parameters, we followed the strongly interacting light Higgs (SILH) scenario [16,26,27], without any modifications in the underlying theory. Further developments were described in [28], where only four phenomenological parameters are used as HZZ couplings. Since the present designs of the circular electron positron colliders and of the ILC are based on a center-of-mass energy around 240 GeV, it is necessary to perform the detailed sensitivity study of the Higgs production for a specific detector design.

HZZ HZγ HWW
In this work, we investigate the polarization angle of the Z decays associated with Higgsstrahlung, as well as the azimuthal angle of the Z decays. Our work is focused on the physics parameters that are extracted from the angular distributions, and on the sensitivity study, which relies on the characteristics of detector design and on the systematical errors of the experiment. A comprehensive study would certainly be helpful at this stage, although the , and events will have to be analyzed separately by the experiments.

HZZ
The paper is organized as follows: In Section 2, analytical formulas for the angular distribution and CP-violation terms for the Z to 2 lepton decays affected by the new coupling are presented. Section 3 gives the numerical limit of the sensitivity of new physics parameters at CEPC. A summary is given in Section 4.

HZZ
The generic effective Hamiltonian in the sector is written as [27] L HZZ = − 1 4 where and . The effective Feynman rule can be derived from Eq (1) as (2) In this parameterization, is the HZZ coupling in the Standard Model. Taking the convention of [28], is a small number in units of , while are small numbers in units of , so that the interaction is consistent with the dimension of mass. The new type of couplings should be smaller than one in SM, since most experimental data are consistent with κ Z SM. The number of free parameters in new physics is then reduced from 12 to only 4 [28], while keeping a sufficiently general structure for the interaction between the Higgs and vector bosons. In contrast to the parametrization [29], which has only one parameter, these four parameters are effective for revealing the details of potential new physics.

HZ
We will focus only on the real Z that is produced in association with the Higgs boson. It decays into a pair of leptons, either or , since they are the particles with the highest detection efficiency and carry the polarization message of a Z boson. Even in lepton colliders it will be hard to tag the electric charge of jets, and we have to choose between the electron or muon as the spin analyzer at a price of reduced statistics. The kinematics of this process is illustrated in Fig. 1. Obviously, the above mentioned new physics coupling of beyond SM may make the cross-section different from SM. Furthermore, the complicated new physics structure in Eq.
(2) may also change the polarization fraction of the Z boson, making the angular distribution of the final lepton pairs different from SM. γ *

ZH
In the Standard Model, the off-shell photon can also contribute to the strahlung production via the for a given design of a lepton collider. In this regard, we do not need to worry about the constraint of the coupling on the electric dipole moment (EDM) of electron [30,31].
The momenta and helicities of the incoming (anti-) electron and outgoing bosons are defined as: where and . The invariant amplitude for the Higgs production is where is the propagator of the virtual Z boson in unitary gauge, and the polarization vector of the real Z is where is the energy of the Z boson.
In the rest frame of the real Z, the decay (helicity) amplitude is written as [32] where is the helicity of the spin analyzer in the Z decay, and is the usual representation of the rotation group. There is also a Breit-Wigner factor, but it is left out as it is an overall factor. The scatting angle , polarization angle and azimuthal angle are defined in Fig. 1.

Total cross-section for Higgsstrahlung
The differential cross-section for Higgs production in the Born approximation reads The kinetic factor K reads: , with s the center-of-mass energy squared, and the momentum of the lepton which is the Z spin analyzer. The definition of spin matrix respects the fact that no beam polarization is expected from the incoming leptons in CEPC. In the case of longitudinal beam polarizations in a linear collider like ILC, a detailed study was published in [17,32].
After integration in the phase space, the total crosssection is: where The SM coupling for the leptonic final states of Z will be defined by the experiment. Since new physics couplings are a small perturbation of the SM couplings, we only keep the leading order linear terms. It is interesting that the anomalous couplings appear as a combination: This relation further reduces the number of free parameters to three, , and . We would also like to point out that this relation takes place at the level of the amplitude of ZH associated production, so it can be regarded as a new parametrization for analyzing the Higgsstrahlung channel. To isolate the contribution, one has to invest-igate the decay channel Higgs into Z pair, whose yield seems relatively small, and is an independent issue beyond the scope of this work.

Polarization in Z boson decays
Although there are only three effective couplings left, one can not distinguish their contribution by the total cross-section measurement only. Different kinds of new physics will give more information on the angular distributions of the decay products of the Z boson, which characterize its polarization fractions. The polar angle distribution of the outgoing lepton is derived as where is the normalized partial width of the Z boson in the helicity state, defined as θ The fraction of each spin polarization, characterized by the distribution of the polarization angle , is obtained by integrating out the scatting angle . In principle, the final lepton coupling may receive extra contributions if anomalous interaction is included. Since CEPC proposes a better option, a Z-factory run to explore this possibility, we keep in this work as in SM . It is interesting to note that the fraction of transverse polarization can be increased if the integration in scattering angle is performed in a reduced interval, for example, in a forward region defined by , Γ ± (ϑ) It is obvious that the contribution from is enhanced by a factor of 3.3 in the forward region. In the experiments, this polarization distribution, together with the total cross-section, will be used to fit the parameters and .

Azimuthal angle distribution for CP violation
Up to now, all the analyses are independent of the CP violation term in the effective Hamiltonian of Eq. (1), which characterizes the CP violating effects in new physics beyond the Standard Model. We need to study the azimuthal angle dependence of the Z boson decays in order to study CP violation effects: sin 2φ cos 2φ Here, the first two terms are background from SM, and CP violation appears in the third term, with the dependence of the signal against the background of .
There is no term in the above equation. However, it can be recovered by breaking the symmetry in decay angle integration, or , at a price of background in SM dσ σdφ sin φ sin 2φ One can see from Fig. 2 that the distribution with reveals the CP violation as breaking of the height equality of the two peaks, while the term makes a phase shift with respect to CP conserving SM background.

Estimation of constraints at future Higgs factories
At future Higgs factories, millions of Higgs events are expected, which will give signals of new physics, or at least provide constraints for new physics, presented in the form of Eq. (2). For example, the circular electron positron collider may deliver a luminosity of at center-of-mass energy . In the conceptual design report [8], the exclusive channel is investigated with phase space cuts: • . Furthermore, CEPC simulations provided the expected performance for: • lepton identification efficiency: 85% • bottom jet tagging efficiency: 75%. We adopt tighter cuts on the phase space and particle tagging (identification), so that the background (mainly ZZ production) can be suppressed to such a level that it can be ignored. At least, its contamination can be well estimated and subtracted in future experiments.
Before Higgs factory data and details of systematical studies become available, we perform Monte Carlo simulations and compare the new physics contribution with SM. Based on the histograms of angular distributions in SM, a Pearson is defined for the event number in each bin based on the hypothesis of new physics. When the parameters , and are sufficiently small, the effects of new physics are hidden by the (mainly statistical) SM errors (reflected in ) , so that the limits of sensitivity can be estimated accordingly.
Instead of making additional assumptions or including more complicated procedures, we assume that the systematical errors are at the same level as the statistical ones, so that the sensitivity limits estimated in following subsections are conservative. Optimized and reliable sensitivity limits should be left to the actual data analysis using experimental input.

Limits from the total cross-section
Using the above mentioned cuts, we scanned the new physics parameters simultaneously. Their sensitivity limits occur when , where is the observed number of events (signal) in , as shown in Fig. 3. The new physics parameters inside the parallelogram will be difficult to distinguish from SM within experimental errors. It can be seen that CEPC can set limits at and .
One may also set lower limits and , by discarding systematical errors and relaxing event selection.
bb On the other hand, the reconstruction of the recoil Z boson leads to an inclusive analysis with Higgs decaying to anything other than the final states. In this case, a tighter limit can be set with about three times larger statistics. In such a reconstruction, it is possible to reduce the dependence on invisible Higgs decays, but the details are beyond the scope of this paper.

Limits from Z polarization ϑ
With the planned luminosity and the analysis outlined at the beginning of this section, we study the sensitivity to new physics using the polarization angle distribution, given in Eq. (12) and (15). The expected event number distribution with the polarization angle is shown in Fig. 4, with blue points for Eq. (12) and black points for Eq. (15). The polarization angle distribution will be different from the black plot in Fig. 4 if only the forward region of the decay angle is taken into account, as in Eq. (15). Since it comes with lower statistics (only half of the number of events), it is omitted in the present analysis until there is better input for experimental systematics.
We show the experimental limits from polarization angle distributions for new physics parameters and in Fig. 5. The parameter region inside the blue lines is indistinguishable from the Standard Model. We also plot the limits from the cross-section study in Fig. 3. Fig. 5 shows that the two limiting regions overlap, which means that the sensitivity limits are further restricted to the meshed region. The polarization angle distribution will anyway be helpful since it will constrain new physics from a different direction than the cross-section. It is also worth to point out that in Eq. (12) the polarization angle distribution is normalized by the cross-section. This means that there is less dependence on the uncertainties from Higgs production or decay, because the actual ana-   lysis is done solely by fitting the shape. As a result, the HZZ couplings are better determined.g

Limits for CP Violation parameter
Using Eq. (16), we show the expected event number distribution with the azimuthal angle in Fig. 6, where CP violation effect may appear. Without loss of generality, only one of the new physics parameters or is scanned at a time together with the CP violation parameter .
Again, the forward region defined in Eq. (17) is omit-g ted in the present analysis due to its lower statistics. After a careful study of the background, we derive the experimental limit for in correlation with , shown in Fig. 7. The correlation sensitivity of and is shown in Fig. 8.
These figures indicate that the experimental sensitivity limit for is .

Summary
e + e − → HZg −0.04 ∼ 0.01 We have studied the sensitivity to new physics in the process at a future Higgs factory. On the basis of measurements of the cross-sections and angular distributions, we set the experimental limits for the dimension-6 operators of the effective field theory in a model independent way. In particular, using the azimuthal angle distribution of the Z boson decays, we found that a future Higgs factory could set a stringent limit on the CP violation effective operators in new physics, with the sensitivity limit of . Our study showed that a future electron positron collider is an ideal machine for the search of new physics signals.