Constraining capability of Zγh production at the ILC Sher Alama∗, Subhasish Beherab†, Satendra Kumarc‡, Shibananda Sahoob§ a Center of High Energy Physics, University of Punjab, Pakistan, Japan-Pakistan Collaboration, 448-18 Ozone, Japan b Department of Physics, Indian Institute of Technology Guwahati, Guwahati 781 039, India and c Indian Institute of Technology Gandhinagar, Gandhinagar 382 355, India Higgs boson couplings with gauge bosons are probed through e−e+ → Zγh in an effective La- grangianframework. ForthisstudythebeampolarizationfacilityattheILCalongwiththetypical center-of-massenergyof500GeVisconsidered. ThereachoftheILCwithanintegratedluminosity of300fb−1 inthedeterminationofCP-conservingparametersisobtained. Sensitivityoftheprobe of each of these couplings in the presence of other couplings is investigated. The most influential coupling parameters are c¯ = −c¯ . Other parameters of significant effect are c¯ and c¯ . A W B HW HB 7 detailedstudyofthevariouskinematicdistributionsrepresentspossibilitiestodisentangletheeffect 1 of some of these couplings. 0 2 PACSnumbers: 12.15.-y,14.70.Fm,13.66.Fg Keywords: electronpositroncollisions,anomalouscouplings,effectiveHiggsLagrangian b e F 1 ] h p - p e h [ 2 v 0 5 2 8 0 . 1 0 7 1 : v i X r a ∗ [email protected] † [email protected] ‡ [email protected] § [email protected] 2 I. INTRODUCTION With the discovery of the new resonance of mass around 125 GeV by the ATLAS and CMS collaborations at the LHC in 2012, a change of paradigm has taken [1–13]. The discovery has unambiguously established the role of the Higgs mechanism in electroweak symmetry breaking (EWSB). All the properties of the new particle measured so far are consistent with that of the standard Higgs boson. Thus, one may be tempted to conclude that for all practical purposes, the newly found particle is like that of the Standard Model (SM) Higgs boson, and new physics effects are decoupled as far as the Higgs sector is concerned. At the same time, it is well known that there are difficulties associated with the Higgs sector of the SM that need to be addressed. The main difficulty is the hierarchy problem associated with the quadratically diverging quantum corrections to the mass of the Higgs boson when computed in the SM. There is no remedy to this difficulty within the SM, and for a Higgs boson of mass 125 GeV, the new physics effects should show up within the TeV range to cure this malady. Assuming that the new physics effects are expected to appear only indirectly in the Higgs sector, it is natural to consider these effects through effective couplings of the Higgs bosons, with itself as well as with the gauge bosons and heavy fermions. Precise measurement of these couplings is very essential to establish the true nature of the EWSB mechanism. While the LHC is capable of probing some of these couplings [14], especially the Higgs couplings with the gauge bosons and top quark, one may need to rely on a cleaner machine like the International Linear Collider (ILC) [15–18] for the required precision. Within an effective Lagrangian, the effect of new physics could be studied in the various couplings through the quantum corrections they acquire. Such an effective Lagrangian basicallyencodesthenewphysicseffectsinhigher-dimensionaloperatorswithanomalouscouplings. The study of the Higgs sector through an effective Lagrangian goes back to Refs.[19–31]. More recently, the Lagrangian including a complete set of dimension-six operators was studied by Refs. [32–35]. For some of the recent references discussing the constraints on the anomalous couplings within different approaches, please see Refs. [36–49]. Reference [47] studied the h+V, where V =Z,W associated production at the LHC and Tevatron to discuss the bounds obtainable from the global fit to the presently available data, whereas Ref. [48] has discussed the constraint on the parameters coming from the LHC results as well as other precision data from LEP, SLC, and Tevatron. Experimental studies on the Higgs couplings at the LHC are presented in, for example, Refs. [50, 51]. The measurement of trilinear Higgs couplings is best done through the process e+e− → Zhh [52–59, 62–64]. At the same time, this process also depends on the Higgs-gauge boson couplings, ZZh and ZZhh, which will affect the determination of the hhh coupling. Another process that could probe the hhh couplings is e+e− → ν ν¯ hh following the WW fusion [56–59], which is also e e affected by the WWh and WWhh couplings. In a recent study [60], we investigated the effect of the VVh coupling, where V = Z, W, in the extraction of the hhh coupling, and found that a precise knowledge of the WWh and ZZh couplings is necessary to derive information regarding the trilinear couplings. Apart from this, an investigation of Higgs-gauge boson couplings has been performedintherecentstudyofWWhproductionatILC[61]. TheanomalouscouplingslikehZγ and hγγ has been studied [65, 66] in the process of e−e+ → hγ. The Table I contains obtained limits on coupling parameters in previous studies. The process e+e− → Zγh is well suited to study the Higgs to neutral gauge boson couplings [52–59, 62–64]. This process is influenced by the trilinear couplings like Zγh, hγγ and ZZhγ which can contaminate the effects of ZZh couplings. In this paper we will focus our attention on this process in some detail within the framework of the effective Lagrangian. One goal of this study is to investigate CP-conserving Higgs sector through Higgs to gauge boson couplings and to understand the significance of other couplings in their measurement. 3 Coefficients Current Bound σ(γh) σ(Zh) c¯ [-0.042, 0.008] [-0.0050, 0.0033] [-1.8, 1.8]×10−4 HW c¯ [-0.053, 0.044] [-0.0033, 0.0050] [-1.8, 1.8]×10−4 HB c¯ [-4.0, 2.3]×10−4 [-0.0012, 0.0028] [-9.0,9.0]×10−4 γ TABLE I: Current 95% CL bounds (2nd column) and future Higgs factory 2σ exclusion sensitivities (3rd and 4th columns) on the coefficients of dim.-6 operators that contribute to e+e− →γh. The paper is presented in the following way. In Sec. II, the effective Lagrangian will be pre- sented, with the currently available constraint on the parameters. In Sec. III, the process under consideration will be presented, with details. In Sec. IV, the results will be summarized. II. GENERAL SETUP References [27–30, 34, 47, 71] present the most general effective Lagrangian with dimension-six operatorsinvolvingtheHiggsbosons. PartofthisLagrangianrelevanttotheprocesse+e− →Zγh considered in this paper is given by Lanom = c¯T (Φ†←D→µΦ)(Φ†←D→ Φ)+ c¯γ g(cid:48)2 Φ†ΦB Bµν Higgs 2v2 µ m2 µν W +c¯HWig (cid:0)DµΦ†σ DνΦ(cid:1)Wk + c¯HBig(cid:48) (cid:0)DµΦ†DνΦ(cid:1)B m2 k µν m2 µν W W + c¯W ig (cid:0)Φ†σ ←D→µΦ(cid:1)DνWk + c¯B ig(cid:48) (cid:0)Φ†←D→µΦ(cid:1)∂νB , (1) 2m2 k µν 2m2 µν W W ←→ where Φ†D Φ = Φ†(D Φ)−(D Φ†)Φ , D being the appropriate covariant derivative operator, µ µ µ µ and Φ, the usual Higgs doublet in the SM. Also, Wk and B are the field tensors corresponding µν µν to the SU(2) and U(1) of the SM gauge groups, respectively, with gauge couplings g and g(cid:48), in L Y that order. σ are the Pauli matrices, and λ is the usual (SM) quartic coupling constant of the k Higgs field. The above Lagrangian, leads to the following in the unitary gauge and mass basis [71] 1 Lanom = − g(1) Z Zµνh−g(2) Z ∂ Zµνh h,Z,γ 4 hZZ µν hZZ ν µ 1 1 + g(3) Z Zµh− g(1) Z Fµνh−g(2) Z ∂ Fµνh 2 hZZ µ 2 hγZ µν hγZ ν µ 1 1 − g(1) Z Fµνh2− g(2) Z ∂ Fµνh2 (2) 4 hhγZ µν 2 hhγZ ν µ Trilinear and quartic CP-conserving couplings g(1) = 2g (cid:2)c¯ sinθ2 −4c¯ sinθ4 +c2 c¯ (cid:3), hzz cosθW2 mW HB W γ W W HW g(2) = g (cid:2)(c¯ +c¯ )cosθ2 +(c¯ +c¯ )sinθ2 (cid:3), g(3) = gmZ [1−2c¯ ] hzz cosθW2 mW HW W W B HB W hzz cosθW T g(1) = gsinθW (cid:2)c¯ −c¯ +8c¯ sinθ2 (cid:3), g(2) = gsinθW [c¯ −c¯ −c¯ +c¯ ] hγz cosθWmW HW HB γ W hγz cosθWmW HW HB B W g(1) = g2sinθW (cid:2)c¯ −c¯ +8c¯ sinθ2 (cid:3), g(2) = g2sinθW [c¯ −c¯ −c¯ +c¯ ] hhγz 2cosθWm2W HW HB γ W hhγz 2cosθWm2W HW HB B W TABLE II: Physical couplings in Eqs. (2) are given in terms of the effective couplings in Eq. (1), where θ being the weak mixing angle. W 4 The physical couplings relevant to the process e+e− → Zγh, and associated with the La- grangian in Eqs. (2) expressed in terms of the effective couplings presented in Eq. (1) are listed in Table II. In total, there are five parameters which are relevant to the process considered, viz, c¯ ,c¯ ,c¯ ,c¯ ,c¯ ,c¯ . These anomalous coefficients c¯ ,c¯ ,c¯ ,c¯ are expected to be of the T γ B W HB HW T HW HB γ order (cid:18)g2 v2(cid:19) (cid:18)g2 M2 (cid:19) c¯ ∼O NP and c¯ ,c¯ ,c¯ ∼O NP W , (3) T M2 HW HB γ 16π2M2 where g denotes the generic coupling of the new physics, and M is the new physics scale. This NP indicates that these couplings can be significantly large for strongly coupled physics. In contrast the coefficients of the operators such as c¯ and c¯ are given by W B (cid:18)m2 (cid:19) c¯ ,c¯ ∼O W (4) B W M2 and therefore, expected to be relatively suppressed or enhanced according to the ratio g/g . NP Comingtotheexperimentalbounds,electroweakprecisiondataputthefollowingconstraints[32], c¯ (m )∈[−1.5,2.2]×10−3 and (c¯ (m )+c¯ (m ))∈[−1.4,1.9]×10−3 (5) T Z W Z B Z Thismeanswecansafelyignoretheeffectofc¯ inouranalysis. Ontheotherhand,c¯ andc¯ are T W B not independently constrained, leaving the possibility of having large values with a cancellation between them as per the above constraint. c¯ itself along with c¯ and c¯ are constrained W HW HB from LHC observations on the associated production of the Higgs along with W in Ref. [47]. Considering the Higgs-associated production along with W, ATLAS and CMS along with D0 put (cid:2) (cid:3) a limit of c¯ ∈ −0.03,0.01 , when all other parameters were set to zero. A global fit using W various information from ATLAS and CMS including signal-strength information constrains the region in the c¯ − c¯ plane, leading to a slightly more relaxed limit on c¯ and a limit of W HW W (cid:2) (cid:3) about c¯ ∈ −0.04,0.01 . The limit on c¯ estimated using a global fit in Ref. [47] is about HW HB c¯ ∈[−0.05,0.05], while a limit of about c¯ ∈[−0.04,0.03] with a one-parameter fit. HB γ The purpose of this study is to understand how to exploit a precision machine like the ILC to investigate a suitable process so as to derive information regarding these couplings. In the next section, we shall explain the process of interest in the present case and discuss the details to understand the influence of one or more of the couplings mentioned above. III. ANALYSES OF THE PROCESS CONSIDERED TheFeynmandiagramscorrespondingtotheprocesse−e+ →ZγhintheSMaregiveninFig.1. This process is basically influenced by Higgs to neutral gauge boson couplings like ZZh, Zγh, hγγ, and ZZγh apart from the fermionic couplings, which are taken to be the standard couplings in our study. FIG. 1: Feynman diagrams contributing to the process e−e+ →Zγh in the SM. The effective Lagrangian, Eq. (1), apart from allowing the existing Higgs and gauge boson cou- plingstobenonstandard,introducesnewcouplingswhichareabsentintheSM.Inaspecificmodel 5 such effects appear at higher orders with a new particle present in the loops. When the masses of such particles are taken to be large, the effect of such quantum correction can be considered in terms of changed couplings. Such effective couplings arising in the present analysis are presented in Table II. Our numerical analyses are carried out using madgraph [68, 69], with the effective Lagrangian implemented through feynrules [70, 71]. 10 SM Unpolarized 8 Pe- = -80%, Pe+ = +30% Pe- = -80%, Pe+ = +60% 6 b] σ [f 4 2 0 200 300 400 500 600 700 800 900 1000 √ √s [GeV] FIG. 2: The total cross section against s in the SM with and without polarized beams. As the first observable, we consider the cross section1. Figure 2 presents the total cross section against the center-of-mass energy for the Zγh production. The cross section peaks around the center-of-mass energy of 350 GeV with a value of about 5.2 fb, which slides down to about 3.8 fb at 500 GeV. In order to avoid any complications arising from the threshold effects, we perform our analysis for an ILC running at a center-of-mass energy of 500 GeV, sufficiently away from the thresholdvalue. Thisisoneoftheplannedcenter-of-massenergyoftheproposedILC.Asexpected, the polarization hugely improves the situation. The case of a typical polarization combination expected at the ILC, 80% left-polarized electron beam and 30% right-polarized positron beam, is considered [17], along with the case of an 80% left-polarized electron beam and a 60% right- polarized positron beams, which are expected in the upgraded version of the ILC. In Fig.3 the cross section against anomalous couplings parameters, c¯ or c¯ or c¯ or c¯ at fixed center- W HW HB γ of-mass energy of 500 GeV is considered along with the role of the polarized beams. In order to be consistent with the experimental constraint [Eq. (5], we choose c¯ = −c¯ throughout our B W analysis, showing the high sensitivity of the cross section on this parameter. Assuming that no other couplings affect the process, the single parameter reach corresponding to the 3σ limit with 300 fb−1 integrated luminosity is given in Table III, while two parameter reach can be seen in Table IV. The obtained limit in the case of unpolarized beam, which is improved with an 80% left-polarized electron beam and a 30% right-polarized positron beam. While the case with an 80% left-polarized electron beam and a 60% right-polarized positron beam does not change this limit significantly, the cross section is increased from about 5.4 fb to about 6.4 fb, enhancing the statistics. In our further analysis, we consider the baseline expectation of an 80% left-polarized electron beam and a 30% right-polarized positron beam. [1]ThemostgeneralformulaisavailableintheAppendix-A. 6 16 √s = 500 GeV Unpolarized 10 √s = 500 GeV unpolarized 14 cW = -cB Pe- = -80%, Pe+ = +30% 9 Pe- = -80%, Pe+ = +30% Pe- = +80%, Pe+ = -30% Pe- = +80%, Pe+ = -30% 12 Pe- = -80%, Pe+ = +60% 8 Pe- = -80%, Pe+ = +60% Pe- = +80%, Pe+ = -60% 7 Pe- = +80%, Pe+ = -60% 10 σ [fb] 8 σ [fb] 56 6 4 4 3 2 2 1 -0.03 -0.02 -0.01 0 0.01 -0.04 -0.03 -0.02 -0.01 0 0.01 cW cHW 10 18 √s = 500 GeV unpolarized √s = 500 GeV unpolarized 16 PPee-- == -+8800%%,, PPee++ == +-3300%% 9 PPee-- == -+8800%%,, PPee++ == +-3300%% 14 Pe- = -80%, Pe+ = +60% 8 Pe- = -80%, Pe+ = +60% Pe- = +80%, Pe+ = -60% Pe- = +80%, Pe+ = -60% 12 7 σ [fb] 10 σ [fb] 6 8 5 6 4 4 3 2 2 -0.05-0.04-0.03-0.02-0.01 0 0.01 0.02 0.03 0.04 0.05 -0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 cHB cγ x100 √ FIG.3: Thetotalcrosssectionagainstanomalouscouplingparametersat s=500GeV,wherethegray, red and blue bands correspond to 3σ deviations from the SM with unpolarized and polarized beams with an integrated luminosity of 300 fb−1, respectively.Note that here, as well as in all figures henceforth, we have removed the “bar” from the symbols denoting the CP-conserving parameters for convenience. Couplings Unpolarized e− =−80%,e+ =+30% e− =−80%,e+ =+60% [-0.022, -0.017] ———— ———— c¯ =−c¯ [-0.003, +0.002] [-0.001, +0.001] [-0.001, +0.0007] W B c¯ [-0.004, +0.003] [-0.0019, +0.0017] [-0.0018, +0.0016] HW c¯ [-0.035, +0.010] [-0.005, +0.005] [-0.004, +0.004] HB c¯ [-0.026, +0.014] [-0.011, +0.014] [-0.009, +0.012] γ TABLEIII:Showingthesingleparameterreachcorrespondingtothe3σlimitwithanintegratedluminosity of 300 fb−1 at center-of-mass energy of 500 GeV, assuming that no other couplings affect the process. 7 9 4 .55 ccWW==cc--HHccBBBBcc γ==γ== =--= 0000- 00....0000..000000111155 √s = 500 GeV 78 √Pse - == 5-0800 %G,e PVe+ = +30% ccWW==cc--HHccBBBBcc γ==γ== =--= 0000- 00....0000..000000111155 6 4 σ [fb] 3.5 σ [fb] 5 4 3 3 2.5 2 2 1 -0.04 -0.03 -0.02 -0.01 0 0.01 -0.04 -0.03 -0.02 -0.01 0 0.01 cHW cHW 8 12 7.5 ccHHWWc γ== =- 00-..000.000111 √s = 500 GeV 11 Pe- = -80%√s, P=e 5+ 0=0 + G3e0V% 6 .75 ccWW==--ccBBc =γ= -= 00 0..00.000111 1 90 ccHHWWc γ== =- 00-..000.000111 cγ = 0.01 6 8 cW=-cB =-0.001 σ [fb] 5.5 σ [fb] 7 cW=-cB = 0.001 5 6 5 4.5 4 4 3 3.5 -0.05-0.04-0.03-0.02-0.01 0 0.01 0.02 0.03 0.04 0.05 -0.05-0.04-0.03-0.02-0.01 0 0.01 0.02 0.03 0.04 0.05 cHB cHB 44 ..568 ccWWc==ccHc--HHHcWcWBBBB = ===== --- 000000......000000000000111551 √s = 500 GeV 67 ..755 √Pse - == 5-0800 %G,e PVe+ = +30%ccWW==cccc-H-HHHccWWBBBB ======--- 000000......000000000000111155 4.4 σ [fb] 4.2 σ [fb] 6 4 5.5 3.8 5 3.6 4.5 3.4 -0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 -0.04 -0.03 -0.02 -0.01 0 0.01 cγ x100 cγ x100 FIG.4: Crosssectionagainst,c¯ (top),c¯ (middle)andc¯ (bottom). Theblacksolidlinecorresponds HW HB γ √ tothecasewhenonlyc¯ orc¯ orc¯ ispresent. Thecenter-of-massenergyisassumedtobe s=500 HW HB γ GeV. In each case, all other parameters are set to zero. The gray band indicates the 3σ limit of the SM crosssectionwithunpolarized(left)andpolarizedbeams(rigth)withanintegratedluminosityof300fb−1. Coming to the Fig.3 & 4, here we present the single parameters analysis in the absence of other anomalous coupling parameters. Clearly, it can be seen that how the obtained limits are affected bythepresenceofotherparameters. Thegraybandindicatesthe3σ limit2 oftheSMcrosssection with unpolarized (left) and polarized beams (right) with an integrated luminosity of 300 fb−1. [2]TheformulacanbeseenintheAppendix-B. 8 Couplings No other couplings c¯ =−c¯ =−0.03 c¯ =−c¯ =+0.01 W B W B c¯ [-0.080, +0.000] [-0.065, +0.055] [-0.105, +0.000] HW c¯ [-0.100, +0.020] [-0.016, +0.010] [-0.011, +0.050] HB c¯ [-0.065, +0.020] [-0.175, +0.095] [-0.060, +0.035] γ c¯ [-0.035, +0.025] [-0.150, +0.032] [-0.020, +0.045] HB Couplings No other couplings c¯ =−0.04 c¯ =+0.01 HW HW c¯ [-0.070, +0.003] [-0.035, +0.025] [-0.082, +0.002] W c¯ [-0.170, +0.015] [-0.135, +0.030] [-0.200, +0.025] HB TABLE IV: Showing the two parameters reach correspond to the 3σ limit with an integrated luminosity of 300 fb−1 at center-of-mass energy of 500 GeV. Fig. 5 presents the correlations between various anomalous coupling paramters like c¯ −c¯ , HW HB c¯ −c¯ ,andc¯ −c¯ ,wheretheyellowandgraybandsshowthepresentlimitsderivedfromthe W HB γ HB LHCresultsontheassociatedproductionoftheHiggsbosonwiththeW boson[47]. Intheabsence ofanyotherparameter,theallowedregioninthec¯ −c¯ planeisrestrictedtoanarrowellipse HW HB (red). This ellipse is not affected much by the presence of c¯ if it is positive (green ellipse). On W the other hand, if c¯ is negative, within the present bounds, it can significantly affect the allowed W region (blue ellipse) in the c¯ −c¯ plane. Similarly, the allowed region in the c¯ −c¯ and HW HB W HB c¯ −c¯ planes are illustrating the the presence of c¯ and c¯ =−c¯ respectivelly. γ HB HW W B 0.05 √s = 500 GeV 0.04 SM (red) 0.15 0.02 ccHHWW ==- 00..0014 ((gbrlueee)n) 0.1 0 0 0.05 W -0.02 0 cH -0.05 cW cγ -0.04 -0.05 -0.1 SM (red) -0.06 -0.1 SM (red) ccWW ==--ccBB ==- 00..0013 ((gbrlueee)n) -0.08 -0.15 ccWW ==--ccBB ==- 00..0013 ((gbrlueee)n) -0.15 -0.1 -0.2 -0.2 -0.15 -0.1 -0.05 0 0.05 -0.2 -0.15 -0.1 -0.05 0 0.05 -0.2 -0.15 -0.1 -0.05 0 0.05 cHB cHB cHB FIG.5: Theellipsescorrespondtoregionsinthec¯ − c¯ (top),c¯ − c¯ (left,bottom)andc¯ − c¯ HB HW HB W HB γ (right, bottom) planes with the total cross section within the 3σ limit of the SM cross section (red), and cross sections with c¯ =−0.03 (blue) and c¯ =+0.01 (green). An integrated luminosity of 300 fb−1 is W W considered, and the center-of-mass energy is taken as 500 GeV. The yellow and gray bands correspond to the present limits of c¯ , c¯ and c¯ , respectively. W HW HB Presently, we would like to be content with the analysis at the production level, considering the limited scope of this work. However, the detector effects may lead to reduction in the sensitivities of the observables studied. It is thus needed to perform a full detector level simulation to estimate the realistic efficiencies of these observables. We have left this as a future work. As mentioned earlier, we shall focus on an ILC running at a center-of-mass energy of 500 GeV for our study. In order to understand the interplay of various CP-conserving couplings, we consider c¯ , c¯ , c¯ W HW HB and c¯ parameters related to anomalous couplings. γ TheeffectoftheanomalouscouplingsonthekinematicdistributionsarepresentedinFigs. 6, 7, 8, 9, 10 & 11. The CP-conserving couplings c¯ , c¯ , c¯ and c¯ parameter choices considered W HW HB γ for these numerical analyses are c¯ =−0.03,+0.01, c¯ =−0.04,+0.01, c¯ =−0.05,+0.05, c¯ =−0.04,+0.03 W HW HB γ 9 While for c¯ , the maximum allowed values as per the present bounds are used, in the case c¯ W HW orc¯ orc¯ ,itissomewhatarbitrarybutwithinthelimits. Whileconsideringbeampolarization, HB γ an 80% left-polarized electron beam and a 30% right-polarized positron beam are assumed, as is expected in the first phase of the ILC, according to the present baseline design. 0.065 0.065 √s = 500 GeV √s = 500 GeV 0.06 0.06 θos-eZ 0.055 θos-eZ 0 0.0.0555 ∆c 0.05 ∆c σσθ1/ (d/dcos) -eZ 00 .0.00.304545 ccWW==c-c-ccHccHHHBBWBBW ======--- 000000......000000551431 σσθ1/ (d/dcos) -eZ 00 0.0.00..03043545 ccWW==c-c-ccHccHHHBBWBBW ======--- 000000......000000551431 cγ x100= 0.03 cγ x100= 0.03 0.03 cγ x100=-0.04 0.025 cγ x100=-0.04 SM SM 0.025 0.02 -1 -0.5 0 0.5 1 -1 -0.5 0 0.5 1 cosθe-Z cosθe-Z FIG.6: Distributionofcosθ fordifferentanomalouscouplingswithunpolarized(left)andpolarizedwith Z P =−80%, P =+30% beams (right). A center-of-mass energy of 500 GeV is assumed. e− e+ 0.065 √s = 500 GeV 0.065 √s = 500 GeV 0.06 0.06 θos-eh 0.055 θos-eh 0.055 ∆c ∆c σσθ (d/dcos) -eh 0 0.00..040455 ccWW==c--ccHccHHBBWBW =====-- 00000.....0000051413 σσθ (d/dcos) -eh 0 0.00..040455 ccWW==c--ccHccHHBBWBW =====-- 00000.....0000054113 1/ cHB=-0.05 1/ cHB=-0.05 0.035 cγ x100= 0.03 0.035 cγ x100= 0.03 cγ x100=-0.04 cγ x100=-0.04 SM SM 0.03 0.03 -1 -0.5 0 0.5 1 -1 -0.5 0 0.5 1 cosθe-h cosθe-h FIG.7: Distributionofcosθ fordifferentanomalouscouplingswithunpolarized(left)andpolarizedwith h P =−80%, P =+30% beams (right). A center-of-mass energy of 500 GeV is assumed. e− e+ 10 0.056 √s = 500 GeV 0.056 √s = 500 GeV 0.054 0.054 θosZh 0.052 θosZh 0.052 σθ∆d/dcos) cZh 00 ..000.440685 ccWW==cc--HcHcBBWW ====-- 0000....00001413 σθ∆d/dcos) cZh 00 ..000.440685 ccWW==cc--HcHcBBWW ====-- 0000....00001413 σ ( cHB= 0.05 σ ( cHB= 0.05 1/ 0.044 cHB=-0.05 1/ 0.044 cHB=-0.05 cγ x100= 0.03 cγ x100= 0.03 0.042 cγ x100=-0.04 0.042 cγ x100=-0.04 SM SM 0.04 0.04 -1 -0.5 0 0.5 1 -1 -0.5 0 0.5 1 cosθZh cosθZh 0.14 0.18 √s = 500 GeV √s = 500 GeV 0.12 cW=-cB = 0.01 0.16 cW=-cB = 0.01 cW=-cB =-0.03 cW=-cB =-0.03 θ∆θdcos) cosγγhh 0 .00.81 ccγγ xxcc11ccHH0H0HWW0B0B======--- 000000S......000000M551443 θ∆θdcos) cosγγhh 000 .0..011.8124 ccγγ xxcc11ccHH0H0HWW0B0B======--- 000000S......000000M551443 σ/ 0.06 σ/ d d σ ( σ ( 0.06 1/ 0.04 1/ 0.04 0.02 0.02 -1 -0.5 0 0.5 1 -1 -0.5 0 0.5 1 cosθγh cosθγh 0.16 0.2 √s = 500 GeV √s = 500 GeV 0.14 cW=-cB = 0.01 0.18 cW=-cB = 0.01 σθ∆θ/dcos) cosγγZZ 00 .0.01.812 cccWγγ =xxcc1-1ccHcH0H0HBWW0B0B =======---- 0000000S.......0000000M5514433 σθ∆θ/dcos) cosγγZZ 0000 .0...0111.81246 cccWγγ =xxcc1-1ccHcH0H0HBWW0B0B =======---- 0000000S.......0000000M5514433 σ1/ (d 0.06 σ1/ (d 0.06 0.04 0.04 0.02 0.02 -1 -0.5 0 0.5 1 -1 -0.5 0 0.5 1 cosθZγ cosθZγ FIG.8: Thenormalizeddistributionsofcosθ (top),cosθ (middle)andcosθ fordifferentanomalous Zh γh Zγ couplings with unpolarized (left) and polarized with P = −80%, P = +30% beams (right). A e− e+ center-of-mass energy of 500 GeV is assumed. We first consider in Fig.6 the normalized cosθ distributions of the Higgs boson for the SM Z case, as well as different cases with anomalous couplings as indicated in the figure, while all other parameters are set to zero. The normalized distributions provide clear information on the shape of the distribution, bringing out the qualitative difference between the different cases considered. The shape of the distribution remains more or less the same as that of the SM case, except a small enhancement in the central regions when c¯ =0.05 (solid blue) in the case of unpolarized HB beams. The advantage of beam polarization is evident (figure on the right) when compared to the corresponding unpolarized (figure on the left) case. Here, the case of negative c¯ differs from the W other cases. This feature can be exploited to discriminate this case from others.