Forward-backward asymmetry in the gauge-Higgs unification at the International Linear Collider

TL;DR

This paper investigates how the gauge-Higgs unification model affects fermion pair production asymmetries at the ILC, highlighting observable deviations from the Standard Model predictions due to Kaluza-Klein modes.

Contribution

It provides the first detailed analysis of forward-backward asymmetry deviations in the gauge-Higgs unification model at the ILC with polarized beams.

Findings

Deviations in asymmetries are significant and detectable.

Differences depend on quark bulk mass parameters.

Predictions for two quark mass parameter cases.

Abstract

The signals of the SO(5) \times U(1) gauge-Higgs unification model at the International Linear Collider are studied. In this model, Kaluza-Klein modes of the neutral gauge bosons affect fermion pair productions. The deviations of the forward-backward asymmetries of the $e^+e^-\to \bar{b}b$, $\bar{t}t$ processes from the standard model predictions are clearly seen by using polarised beams. The deviations of these values are predicted for two cases, the bulk mass parameters of quarks are positive and negative case.

Figures

Tables

Table 1. Table 1 : Input parameters. Masses of Z 𝑍 Z boson, the Higgs boson, leptons and quarks in the unit of GeV at the Z 𝑍 Z mass scale.

${\alpha }_{\text{EM}}^{-1}$	${\sin }^2{\theta }_W$
$127.96$	$0.23122$
$m_Z$	$m_H$	$m_{\nu }$	$m_e$	$m_{\mu }$	$m_{\tau }$
$91.188$	$125$	${10}^{-12}$	$0.48657\times {10}^{-3}$	$102.72\times {10}^{-3}$	$1.7462$
$m_u$	$m_d$	$m_s$	$m_c$	$m_b$	$m_t$
$1.27\times {10}^{-3}$	$2.90\times {10}^{-3}$	$0.055$	$0.619$	$2.89$	$171.7$

Table 2. Table 2 : The warp factor, the bulk mass parameters of the fermion and AdS curvature with each values of θ H subscript 𝜃 𝐻 \theta_{H} are listed. The resultant W 𝑊 W -boson mass, Higgs boson mass and KK scale given from the model parameters are also summarized.

${\theta }_H$	$e^{kL}$	$\left|c_F\right|$	$k$ (GeV)	$m_W$ (GeV)	$m_H$ (GeV)	$m_{\text{KK}}$ (GeV)
$0.10$	$2.90\times {10}^4$	$0.29617$	$7.4431\times {10}^7$	$79.957$	$125.1$	$8063$
$0.09$	$1.70\times {10}^4$	$0.27670$	$4.7190\times {10}^7$	$79.958$	$125.1$	$8721$
$0.08$	$1.01\times {10}^4$	$0.25356$	$3.0679\times {10}^7$	$79.951$	$125.4$	$9544$
${\theta }_H$	$\left|c_e\right|$	$\left|c_{\mu }\right|$	$\left|c_{\tau }\right|$	$\left|c_u\right|$	$\left|c_c\right|$	$\left|c_t\right|$
$0.10$	$1.8734$	$1.3139$	$1.0060$	$1.6796$	$1.1200$	$0.16116$
$0.09$	$1.9504$	$1.3599$	$1.0348$	$1.7459$	$1.1552$	$0.11646$
$0.08$	$2.0342$	$1.4100$	$1.0663$	$1.8180$	$1.1936$	$0.0089140$

Table 3. Table 3 : Couplings of neutral vector bosons ( Z ′ superscript 𝑍 ′ Z^{\prime} bosons) to fermions in unit of g w = e / sin ⁡ θ W subscript 𝑔 𝑤 𝑒 subscript 𝜃 𝑊 g_{w}=e/\sin\theta_{W} for θ H = 0.10 subscript 𝜃 𝐻 0.10 \theta_{H}=0.10 , c l subscript 𝑐 𝑙 c_{l} , c q > 0 subscript 𝑐 𝑞 0 c_{q}>0 .

$f$	$g_{Zf}^L$	$g_{Zf}^R$	$g_{Z^{\left(1\right)}f}^L$	$g_{Z^{\left(1\right)}f}^R$	$g_{Z_R^{\left(1\right)}}^L$	$g_{Z_R^{\left(1\right)}f}^R$	$g_{{\gamma }^{\left(1\right)}f}^L$	$g_{{\gamma }^{\left(1\right)}f}^R$
${\nu }_e$	$+0.57037$	$0$	$-0.20943$	$0$	$0$	$0$	$0$	$0$
${\nu }_{\mu }$	$+0.57037$	$0$	$-0.20943$	$0$	$0$	$0$	$0$	$0$
${\nu }_{\tau }$	$+0.57037$	$0$	$-0.20928$	$0$	$0$	$0$	$0$	$0$
$e$	$-0.30661$	$+0.26384$	$+0.11258$	$+1.04332$	$0$	$-1.4357$	$+0.17720$	$-1.8962$
$\mu$	$-0.30661$	$+0.26384$	$+0.11258$	$+0.97948$	$0$	$-1.3582$	$+0.17720$	$-1.7801$
$\tau$	$-0.30661$	$+0.26383$	$+0.11250$	$+0.92684$	$0$	$-1.2940$	$+0.17708$	$-1.6844$
$u$	$+0.39453$	$-0.17589$	$-0.14486$	$-0.68311$	$0$	$+0.94208$	$-0.11813$	$+1.2415$
$c$	$+0.39453$	$-0.17589$	$-0.14485$	$-0.63219$	$0$	$+0.88013$	$-0.11812$	$+1.1489$
$t$	$+0.39353$	$-0.17694$	$+0.57109$	$-0.42117$	$+1.1369$	$+0.62142$	$+0.46722$	$+0.76730$
$d$	$-0.48245$	$+0.087946$	$+0.17715$	$+0.34156$	$0$	$-0.47104$	$+0.059066$	$-0.62077$
$s$	$-0.48245$	$+0.087945$	$+0.17713$	$+0.31609$	$0$	$-0.44006$	$+0.059060$	$-0.57445$
$b$	$-0.48252$	$+0.087939$	$-0.70659$	$+0.21112$	$+1.1347$	$-0.31045$	$-0.23377$	$-0.38353$

Table 4. Table 6 : Z ′ superscript 𝑍 ′ Z^{\prime} couplings of fermions for θ H = 0.10 subscript 𝜃 𝐻 0.10 \theta_{H}=0.10 , c l subscript 𝑐 𝑙 c_{l} , c q < 0 subscript 𝑐 𝑞 0 c_{q}<0 . The same unit as in Table  5 .

$f$	$g_{Zf}^L$	$g_{Zf}^R$	$g_{Z^{\left(1\right)}f}^L$	$g_{Z^{\left(1\right)}f}^R$	$g_{Z_R^{\left(1\right)}}^L$	$g_{Z_R^{\left(1\right)}f}^R$	$g_{{\gamma }^{\left(1\right)}f}^L$	$g_{{\gamma }^{\left(1\right)}f}^R$
${\nu }_e$	$+0.57054$	$0$	$+2.2561$	$0$	$-3.1047$	$0$	$0$	$0$
${\nu }_{\mu }$	$+0.57053$	$0$	$+2.1181$	$0$	$-2.9371$	$0$	$0$	$0$
${\nu }_{\tau }$	$+0.57052$	$0$	$+2.0042$	$0$	$-2.7982$	$0$	$0$	$0$
$e$	$-0.30384$	$+0.26662$	$-1.2015$	$-0.09790$	$-3.1130$	$0$	$-1.8962$	$+0.17720$
$\mu$	$-0.30384$	$+0.26662$	$-1.1280$	$-0.09790$	$-2.9450$	$0$	$-1.7801$	$+0.17720$
$\tau$	$-0.30383$	$+0.26662$	$-1.0674$	$-0.09783$	$-2.8057$	$0$	$-1.6844$	$+0.17708$
$u$	$+0.39419$	$-0.17630$	$+1.5309$	$+0.06473$	$-1.3595$	$0$	$+1.2415$	$-0.11813$
$c$	$+0.39180$	$-0.17868$	$+1.4082$	$+0.06560$	$+2.3878$	$0$	$+1.1489$	$-0.11812$
$t$	$+0.39281$	$-0.17766$	$+0.9400$	$-0.25520$	$+1.7082$	$+0.41431$	$+0.7673$	$+0.46722$
$d$	$-0.48019$	$+0.09032$	$-1.8649$	$-0.03316$	$-1.3614$	$0$	$-0.62077$	$+0.05907$
$s$	$-0.48256$	$+0.08794$	$-1.7344$	$-0.03229$	$+2.3805$	$0$	$-0.57445$	$+0.05906$
$b$	$-0.48255$	$+0.08793$	$-1.1585$	$+0.12877$	$+1.7027$	$-0.20689$	$-0.38353$	$-0.23377$

Table 5. Table 9 : Masses of Z ( 1 ) superscript 𝑍 1 Z^{(1)} , Z R ( 1 ) superscript subscript 𝑍 𝑅 1 Z_{R}^{(1)} and γ ( 1 ) superscript 𝛾 1 \gamma^{(1)} and total decay width of γ ( 1 ) superscript 𝛾 1 \gamma^{(1)} in the unit of GeV. Γ γ ( 1 ) subscript Γ superscript 𝛾 1 \Gamma_{\gamma^{(1)}} is independent of the sign of the bulk fermion parameters. Γ Z ( 1 ) / Z R ( 1 ) ​ ( ± , ± ) subscript Γ superscript 𝑍 1 superscript subscript 𝑍 𝑅 1 plus-or-minus plus-or-minus \Gamma_{Z^{(1)}/Z_{R}^{(1)}}(\pm,\pm) represent that left and right sign is sign of c l subscript 𝑐 𝑙 c_{l} and c q subscript 𝑐 𝑞 c_{q} .

${\theta }_H$	$m_{Z^{\left(1\right)}}$	$m_{Z_R^{\left(1\right)}}$	$m_{{\gamma }^{\left(1\right)}}$	${\Gamma }_{{\gamma }^{\left(1\right)}}$
$0.10$	$6585$	$6172$	$6588$	$905$
$0.09$	$7149$	$6676$	$7152$	$940$
$0.08$	$7855$	$7305$	$7858$	$986$
${\theta }_H$	${\Gamma }_{Z^{\left(1\right)}}(+,+)$	${\Gamma }_{Z^{\left(1\right)}}(+,-)$	${\Gamma }_{Z^{\left(1\right)}}(-,+)$	${\Gamma }_{Z^{\left(1\right)}}(-,-)$
$0.10$	$429$	$1632$	$959$	$2162$
$0.09$	$463$	$1674$	$1014$	$2225$
$0.08$	$534$	$1705$	$1112$	$2283$
${\theta }_H$	${\Gamma }_{Z_R^{\left(1\right)}}(+,+)$	${\Gamma }_{Z_R^{\left(1\right)}}(+,-)$	${\Gamma }_{Z_R^{\left(1\right)}}(-,+)$	${\Gamma }_{Z_R^{\left(1\right)}}(-,-)$
$0.10$	$784$	$2437$	$2398$	$4051$
$0.09$	$856$	$2480$	$2529$	$4153$
$0.08$	$1005$	$2485$	$2758$	$4238$

Table 6. Table 10 : Deviations of the cross sections at s = 250 𝑠 250 \sqrt{s}=250 GeV with P e − = + 0.8 subscript 𝑃 superscript 𝑒 0.8 P_{e^{-}}=+0.8 and P e + = − 0.3 subscript 𝑃 superscript 𝑒 0.3 P_{e^{+}}=-0.3 and 250 fb -1 luminosity for e + ​ e − → q ¯ ​ q ​ ( q ≠ b , t ) → superscript 𝑒 superscript 𝑒 ¯ 𝑞 𝑞 𝑞 𝑏 𝑡 e^{+}e^{-}\to\bar{q}q\leavevmode\nobreak\ (q\neq b,t) and e + ​ e − → b ¯ ​ b → superscript 𝑒 superscript 𝑒 ¯ 𝑏 𝑏 e^{+}e^{-}\to\bar{b}b processes. The statistical uncertainties calculated by the SM prediction are 0.074 % for q ¯ ​ q ¯ 𝑞 𝑞 \bar{q}q and 0.20 % for b ¯ ​ b ¯ 𝑏 𝑏 \bar{b}b .

${\theta }_H$	${\sigma }^{c_q>0}/{\sigma }^{\text{SM}}\left(\overline{q}q\right)-1$
$0.10$	$-4.56\%$ ($-61.5\sigma$)	$+0.16\%$ ($+2.21\sigma$)
$0.09$	$-3.70\%$ ($-49.9\sigma$)	$+0.14\%$ ($+1.90\sigma$)
$0.08$	$-2.98\%$ ($-40.0\sigma$)	$+0.11\%$ ($+1.53\sigma$)
${\theta }_H$	${\sigma }^{c_q>0}/{\sigma }^{\text{SM}}\left(\overline{b}b\right)-1$
$0.10$	$-4.18\%$ ($-21.1\sigma$)	$-3.96\%$ ($-19.9\sigma$)
$0.09$	$-3.41\%$ ($-17.2\sigma$)	$-3.29\%$ ($-16.6\sigma$)
$0.08$	$-2.84\%$ ($-14.3\sigma$)	$-2.65\%$ ($-13.3\sigma$)

Table 7. Table 11 : Deviation of A FB ​ ( c ¯ ​ c ) subscript 𝐴 FB ¯ 𝑐 𝑐 A_{\text{FB}}(\bar{c}c) and A FB ​ ( b ¯ ​ b ) subscript 𝐴 FB ¯ 𝑏 𝑏 A_{\text{FB}}(\bar{b}b) at s = 250 𝑠 250 \sqrt{s}=250 GeV with P eff = + 0.887 subscript 𝑃 eff 0.887 P_{\text{eff}}=+0.887 and 250 fb -1 luminosity. The statistical uncertainties σ 𝜎 \sigma calculated by the event number of the SM prediction are 0.16 % and 0.70 %, respectively.

${\theta }_H$	$A_{\text{FB}}^{c_q>0}/A_{\text{FB}}^{\text{SM}}\left(\overline{c}c\right)-1$
$0.10$	$-0.38\%\left(-2.32\sigma \right)$	$+1.68\%\left(+10.1\sigma \right)$
$0.09$	$-0.31\%\left(-1.88\sigma \right)$	$+1.25\%\left(+7.58\sigma \right)$
$0.08$	$-0.25\%\left(-1.49\sigma \right)$	$+1.03\%\left(+6.21\sigma \right)$
${\theta }_H$	$A_{\text{FB}}^{c_q>0}/A_{\text{FB}}^{\text{SM}}\left(\overline{b}b\right)-1$
$0.10$	$+4.24\%\left(+6.04\sigma \right)$	$+7.33\%\left(+10.4\sigma \right)$
$0.09$	$+4.07\%\left(+5.81\sigma \right)$	$+5.69\%\left(+8.10\sigma \right)$
$0.08$	$+4.15\%\left(+5.92\sigma \right)$	$+3.83\%\left(+5.45\sigma \right)$

Table 8. Table 12 : Deviations of A FB ​ ( c ¯ ​ c ) subscript 𝐴 FB ¯ 𝑐 𝑐 A_{\text{FB}}(\bar{c}c) , A FB ​ ( b ¯ ​ b ) subscript 𝐴 FB ¯ 𝑏 𝑏 A_{\text{FB}}(\bar{b}b) and A FB ​ ( t ¯ ​ t ) subscript 𝐴 FB ¯ 𝑡 𝑡 A_{\text{FB}}(\bar{t}t) at s = 500 𝑠 500 \sqrt{s}=500  GeV with P eff = + 0.887 subscript 𝑃 eff 0.887 P_{\text{eff}}=+0.887 and 500 fb -1 luminosity. The statistical uncertainties calculated by the event number of the SM prediction are 0.26 %, 0.79 % and 0.53 %, respectively.

${\theta }_H$	$A_{\text{FB}}^{c_q>0}/A_{\text{FB}}^{\text{SM}}\left(\overline{c}c\right)-1$
$0.10$	$-3.25\%\left(-12.3\sigma \right)$	$+7.49\%\left(+28.3\sigma \right)$
$0.09$	$-2.58\%\left(-9.77\sigma \right)$	$+6.50\%\left(+24.6\sigma \right)$
$0.08$	$-2.01\%\left(-7.59\sigma \right)$	$+5.53\%\left(+20.9\sigma \right)$
${\theta }_H$	$A_{\text{FB}}^{c_q>0}/A_{\text{FB}}^{\text{SM}}\left(\overline{b}b\right)-1$
$0.10$	$+15.4\%\left(+19.4\sigma \right)$	$+23.1\%\left(+29.1\sigma \right)$
$0.09$	$+14.3\%\left(+18.1\sigma \right)$	$+18.4\%\left(+23.3\sigma \right)$
$0.08$	$+14.1\%\left(+17.8\sigma \right)$	$+12.9\%\left(+16.3\sigma \right)$
${\theta }_H$	$A_{\text{FB}}^{c_q>0}/A_{\text{FB}}^{\text{SM}}\left(\overline{t}t\right)-1$
$0.10$	$+5.36\%\left(+9.96\sigma \right)$	$+9.25\%\left(+17.2\sigma \right)$
$0.09$	$+5.03\%\left(+9.35\sigma \right)$	$+7.23\%\left(+13.4\sigma \right)$
$0.08$	$+5.14\%\left(+9.55\sigma \right)$	$+4.91\%\left(+9.13\sigma \right)$

Table 9. Table 13 : Deviation of A FB ​ ( c ¯ ​ c ) subscript 𝐴 FB ¯ 𝑐 𝑐 A_{\text{FB}}(\bar{c}c) , A FB ​ ( b ¯ ​ b ) subscript 𝐴 FB ¯ 𝑏 𝑏 A_{\text{FB}}(\bar{b}b) and A FB ​ ( t ¯ ​ t ) subscript 𝐴 FB ¯ 𝑡 𝑡 A_{\text{FB}}(\bar{t}t) at s = 1 𝑠 1 \sqrt{s}=1  TeV with P eff = + 0.887 subscript 𝑃 eff 0.887 P_{\text{eff}}=+0.887 and 1000 fb -1 luminosity. The statistical uncertainties calculated by the event number of the SM prediction are 0.39 %, 1.07 % and 0.45 %, respectively.

${\theta }_H$	$A_{\text{FB}}^{c_q>0}/A_{\text{FB}}^{\text{SM}}\left(\overline{c}c\right)-1$
$0.10$	$-25.9\%\left(-67.3\sigma \right)$	$-2.46\%\left(-6.37\sigma \right)$
$0.09$	$-18.3\%\left(-47.6\sigma \right)$	$+4.61\%\left(+12.0\sigma \right)$
$0.08$	$-12.8\%\left(-33.3\sigma \right)$	$+9.02\%\left(+23.4\sigma \right)$
${\theta }_H$	$A_{\text{FB}}^{c_q>0}/A_{\text{FB}}^{\text{SM}}\left(\overline{b}b\right)-1$
$0.10$	$+46.1\%\left(+43.1\sigma \right)$	$+21.1\%\left(+19.7\sigma \right)$
$0.09$	$+45.7\%\left(+42.8\sigma \right)$	$+36.8\%\left(+34.4\sigma \right)$
$0.08$	$+44.9\%\left(+42.0\sigma \right)$	$+38.8\%\left(+36.3\sigma \right)$
${\theta }_H$	$A_{\text{FB}}^{c_q>0}/A_{\text{FB}}^{\text{SM}}\left(\overline{t}t\right)-1$
$0.10$	$+11.3\%\left(+24.3\sigma \right)$	$+9.36\%\left(+20.6\sigma \right)$
$0.09$	$+11.6\%\left(+25.7\sigma \right)$	$+12.3\%\left(+27.1\sigma \right)$
$0.08$	$+12.2\%\left(+26.9\sigma \right)$	$+11.3\%\left(+25.0\sigma \right)$

Equations

$$d{s}^2=g_{MN}d{x}^{M}d{x}^N=e^{-2\sigma (y)}{\eta }_{\mu \nu }d{x}^{\mu }d{x}^{\nu }+d{y}^2,$$

$$S_{\text{bulk}}$$

$$+\frac{1}{2{\xi }_{(A)}}{\left(f_{\text{gf}}^{(A)}\right)}^2+\frac{1}{2{\xi }_{(B)}}{\left(f_{\text{gf}}^{(B)}\right)}^2+\frac{1}{2{\xi }_{(G)}}{\left(f_{\text{gf}}^{(G)}\right)}^2+{\mathcal{L}}_{GH}^{(A)}+{\mathcal{L}}_{GH}^{(B)}+{\mathcal{L}}_{GH}^{(G)}$$

$$\displaystyle+\sum_{g=1}^{3}\sum_{a=1}^{4}\bar{\Psi}^{g}_{a}\mathcal{D}(c^{g}_{a})\Psi^{g}_{a}+\sum_{i=1}^{N_{F}}\bar{\Psi}_{F_{i}}\mathcal{D}(c_{F_{i}})\Psi_{F_{i}}\biggr{\}},$$

$$\mathcal{D}(c_a)$$

$$\displaystyle\hskip 60.00009pt-ig_{A}A_{M}-ig_{B}Q_{X,a}B_{M}-ig_{C}Q_{C,a}G_{M}\bigg{)},$$

$$A_M=\sum _{a_L=1}^3{A}_M^{a_L}{T}^{a_L}+\sum _{a_R=1}^3{A}_M^{a_R}{T}^{a_R}+\sum _{\hat{a}=1}^4{A}_M^{\hat{a}}{T}^{\hat{a}},$$

$${\Gamma }^{\mu }={\gamma }^{\mu }=\left(\begin{array}{cc}& {\sigma }^{\mu }\\ {\bar{\sigma }}^{\mu }& \end{array}\right)\text{leavevmode},\text{leavevmode}\text{leavevmode}{\Gamma }^5={\gamma }^5=\left(\begin{array}{cc}1& \\ & -1\end{array}\right)\text{leavevmode},$$

$${\sigma }^{\mu }=(1,\vec{\sigma })\text{leavevmode},\text{leavevmode}\text{leavevmode}{\bar{\sigma }}^{\mu }=(-1,\vec{\sigma })\text{leavevmode}.$$

$${\Psi }_1$$

$${\Psi }_2$$

$${\Psi }_3$$

$${\Psi }_4$$

$$\left(\begin{array}{c}{A}_{\mu }\\ A_y\end{array}\right)(x,y_j-y)$$

$$\left(\begin{array}{c}{B}_{\mu }\\ B_y\end{array}\right)(x,y_j-y)$$

$$\left(\begin{array}{c}{G}_{\mu }\\ G_y\end{array}\right)(x,y_j-y)$$

$${\Psi }_a^g(x,y_j-y)$$

$${\Psi }_{F_i}(x,y_j-y)$$

$$P_{\text{ve}}=\text{diag}(-1,-1,-1,-1,$$

$$S_{\text{brane}}$$

$$\displaystyle+\sum_{\alpha=1}^{3}\big{(}\hat{\chi}^{q\dagger}_{\alpha R}\,i\bar{\sigma}^{\mu}D_{\mu}\hat{\chi}^{q}_{\alpha R}+\hat{\chi}^{l\dagger}_{\alpha R}i\bar{\sigma}^{\mu}D_{\mu}\hat{\chi}^{l}_{\alpha R}\big{)}$$

$$\displaystyle-i\Bigl{[}\kappa_{1}^{q}\hat{\chi}_{1R}^{q\dagger}\hat{\Psi}_{1L}\tilde{\hat{\Phi}}+\kappa_{2}^{q}\hat{\chi}_{2R}^{q\dagger}\hat{\Psi}_{2L}\tilde{\hat{\Phi}}+\kappa_{3}^{q}\hat{\chi}_{3R}^{q\dagger}\hat{\Psi}_{2L}\hat{\Phi}+\tilde{\kappa}^{q}\hat{\chi}_{2R}^{q\dagger}\hat{\Psi}_{1L}\hat{\Phi}-(\text{h.c.})\Bigr{]}$$

$$\displaystyle-i\Bigl{[}\kappa_{1}^{l}\hat{\chi}_{1R}^{l\dagger}\hat{\Psi}_{3L}\hat{\Phi}+\kappa_{2}^{l}\hat{\chi}_{2R}^{l\dagger}\hat{\Psi}_{4L}\tilde{\hat{\Phi}}+\kappa_{3}^{l}\hat{\chi}_{3R}^{l\dagger}\hat{\Psi}_{4L}\hat{\Phi}+\tilde{\kappa}^{l}\hat{\chi}_{3R}^{l\dagger}\check{\Psi}_{3L}\tilde{\hat{\Phi}}-(\text{h.c.})\Bigr{]}\bigg{\}},$$

$$\displaystyle D_{\mu}\hat{\Phi}=\Big{(}\partial_{\mu}-ig_{A}\sum^{3}_{a_{R}=1}A^{a_{R}}_{\mu}T^{a_{R}}-ig_{B}Q_{X}B_{\mu}\Big{)}\hat{\Phi}\leavevmode\nobreak\ ,$$

$$\displaystyle D_{\mu}\hat{\chi}_{\alpha R}=\Big{(}\partial_{\mu}-ig_{A}\sum^{3}_{a_{L}=1}A^{a_{L}}_{\mu}T^{a_{L}}-ig_{B}Q_{X,\alpha}B_{\mu}-ig_{C}Q_{C,\alpha}G_{\mu}\Big{)}\hat{\chi}_{\alpha R},$$

$$\tilde{\hat{\Phi }}=i{\sigma }_2{\hat{\Phi }}^{\ast },$$

$${\hat{\chi }}_{1R}^q$$

$${\hat{\chi }}_{1R}^{\ell }$$

$$⟨\Phi ⟩=\left(\begin{array}{c}0\\ w\end{array}\right).$$

$${\Psi }_{L,R}(x,y)=\frac{e^{\frac{3}{2}ky}}{\sqrt{L}}\sum _{n=0}^{\infty }{\psi }_{L,R}^{(n)}(x)\frac{f_{L,R}^{(n)}(y)}{\sqrt{N^{(n)}}},$$

$$N^{(n)}={\int }_0^L\frac{dy}{L}{\left\{f_{L,R}^{(n)}(y)\right\}}^2,$$

$$f_{L,R}^{(0)}(y)$$

$$f_{L,R}^{(n)}(y)$$

$$(t_L:B_L:U_L)=\left(\frac{1+c_H}{2}:\frac{1-c_H}{2}:-\frac{\tilde{\kappa }}{{\kappa }_2}\right).$$

$$t_L:\left(+\frac{1}{2},-\frac{1}{2}\right),B_L:\left(-\frac{1}{2},+\frac{1}{2}\right),U_L:\left(+\frac{1}{2},+\frac{1}{2}\right)$$

$${\left(\frac{1+c_H}{2}\right)}^2+{\left(\frac{1-c_H}{2}\right)}^2+{\left(\frac{m_b}{m_t}\right)}^2,$$

$${\left(\frac{1+c_H}{2}\right)}^2-{\left(\frac{1-c_H}{2}\right)}^2+{\left(\frac{m_b}{m_t}\right)}^2,$$

$$-{\left(\frac{1+c_H}{2}\right)}^2+{\left(\frac{1-c_H}{2}\right)}^2+{\left(\frac{m_b}{m_t}\right)}^2.$$

$$\frac{d\sigma }{d\cos \theta }$$

Full text

1 Introduction

After the discovery of the Higgs boson, search for new physics in the electroweak sector is one of the most important topics of the particle physics. For this purpose, high-precision measurements of the electroweak sector are necessary. The International Linear Collider (ILC) has the capabilities for the high-precision measurements and enables us to test the standard model (SM) [1, 2, 3, 4]. New physics at the ILC are predicted by many alternative models such as the Higgs portal dark matter models [5, 6, 7], two-Higgs doublet models [8, 9, 10], Georgi-Machacek model [11], supersymmetric models [12, 13, 14, 15], littlest Higgs model [16, 17], universal extra dimensional model [18], warped extra dimensional models [19, 20, 21], composite Higgs models [22, 23, 24] and other models [25].

The gauge-Higgs unification (GHU) models [26, 27, 28, 29, 30, 31, 44, 32, 45, 46, 35, 47, 36, 38, 41, 42, 43, 33, 37, 48, 49, 39, 40, 34, 50, 51, 52, 53] are alternative models of the SM, in which the Higgs boson appears as an extra-dimensional component of the higher-dimensional gauge boson. Hence the Higgs sector is governed by the gauge principle and the Higgs boson mass is protected against the radiative corrections in the GHU models. The Higgs boson is massless at the tree-level and acquires the finite mass by the radiative corrections [26, 27, 28]. Note that also the Yukawa interactions in the GHU models are the gauge interactions of the extra-dimensional component of the gauge boson and the fermions. The phenomenologically most well-studied GHU model is the $SO(5)\times U(1)$ GHU models [31, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53], which are defined on the warped metric [60]. In the $SO(5)\times U(1)$ GHU models, the Higgs boson mass is protected by the custodial symmetry and the Higgs boson potential is calculable. One of the features of the warped extra dimensional models is that the first Kaluza-Klein (KK) excited states of the gauge bosons have the large asymmetry between the couplings with the left-handed and the right-handed fermions. The $SO(5)\times U(1)$ GHU models have the same feature.

There are several kinds of the $SO(5)\times U(1)$ GHU models depending on symmetry breaking patterns and the embedding patterns of the quarks and leptons. In the $SO(5)\times U(1)$ GHU model discussed here, the quarks and leptons are embedded in the 5 representation of $SO(5)$. Then the Yukawa, $HWW$ and $HZZ$ couplings are suppressed by $\cos\theta_{H}$ from the SM values, where $\theta_{H}$ is the Wilson line phase. The KK excited states of the neutral gauge bosons appear as the so-called $Z^{\prime}$ bosons. Besides the KK excited states of the photon $\gamma^{(n)}$ and that of the $Z$ boson $Z^{(n)}$, that of the $SU(2)_{R}$ gauge boson $Z_{R}^{(n)}$, which does not have a zero mode, exist as the neutral gauge bosons. Therefore the $\gamma^{(n)}$, $Z^{(n)}$ and $Z_{R}^{(n)}$ appear as the $Z^{\prime}$ bosons. From the result at the Large Hadron Collider, the parameter region is constrained to be $\theta_{H}\lesssim 0.1$, where $m_{\text{KK}}\gtrsim 8$ TeV [48]. For this region, the contributions of the KK excited states to the decay widths, $\Gamma(H\to\gamma\gamma)$ and $\Gamma(H\to Z\gamma)$ are less than 0.2%, thus negligible [44, 47]. The decay widths $\Gamma(H\to WW)$, $\Gamma(H\to ZZ)$, $\Gamma(H\to\bar{q}q)$, $\Gamma(H\to l^{+}l^{-})$, $\Gamma(H\to\gamma\gamma)$ and $\Gamma(H\to Z\gamma)$ are approximately suppressed by the common factor $\cos^{2}\theta_{H}$ at the leading order. Hence the branching ratio of the Higgs boson in this model is almost equivalent to that in the SM. At the ILC 500 GeV, the branching ratios of the processes $H\to WW$, $H\to ZZ$, $H\to\bar{b}b$, $H\to\bar{c}c$ and $H\to\tau^{+}\tau^{-}$ are measured at the $O(1)\%$ accuracy [2]. The branching ratios of the Higgs boson are consistent between the GHU model and the SM at the ILC.

At the ILC, it is possible to measure effects of new physics on the the cross sections and the forward-backward asymmetries of the $e^{+}e^{-}\to f\bar{f}$ processes [16, 17, 19, 20, 54, 55, 56, 57]. The forward-backward asymmetries are important observables for indirect search of physics beyond the SM as shown in the studies of that of top quark production at the Tevatron and Large Hadron Collider [58, 59]. In this paper, the effects of $Z^{\prime}$ bosons in the GHU model on these values are shown. The same topic is studied previously [49] and clear signals are predicted although $Z^{\prime}$ masses are much larger than the centre-of-mass of the ILC. However in the previous study, the bulk mass parameters are assumed to be positive although negative values are also allowed. Thus the negative region of the bulk mass parameters is also checked in this study. The following calculations are done at the tree-level without the one-loop effective potential of the Higgs boson, so the corrections from the strong interaction are not included.

This paper is constructed as follows. In Section 2, the model is shortly introduced. In Section 3, the parameters and the couplings and decay widths of the $Z^{\prime}$ bosons are shown. The bulk mass parameter dependence of the fermion mode functions are also shortly reviewed. In Section 4, the cross sections and the forward-backward asymmetries at the ILC are shown. In section 5, the results are summarised.

2 Model

There are several types of the $SO(5)\times U(1)$ GHU models depending on the symmetry breaking patterns and the embedding patterns of the quarks and leptons. The $SO(5)\times U(1)$ GHU model discussed in this paper is constructed as follows. The model is defined on the warped spacetime [60]. The metric is

[TABLE]

($0\leq|y|\leq+L$) where $k$ is the AdS curvature. The action has the $SO(5)\times U(1)_{X}$ local symmetry. The bulk action is given by

[TABLE]

where $A_{M},B_{M}$ and $G_{M}$ are the $SO(5)$, $U(1)_{X}$ and $SU(3)_{C}$ gauge fields, $F_{MN}^{(A)}=\partial_{M}A_{N}-\partial_{N}A_{M}-ig_{A}[A_{M},A_{N}]$, $F_{MN}^{(B)}=\partial_{M}B_{N}-\partial_{N}B_{M}$ and $F_{MN}^{(G)}=\partial_{M}G_{N}-\partial_{N}G_{M}-ig_{C}[G_{M},G_{N}]$, $g_{A}$, $g_{B}$ and $g_{C}$ are the five-dimensional gauge couplings of $SO(5)$, $U(1)_{X}$ and $SU(3)_{C}$. $f_{\text{gf}}^{(A)}$, $f_{\text{gf}}^{(B)}$ and $f_{\text{gf}}^{(G)}$ are gauge-fixing functions and $\xi_{(A)}$, $\xi_{(B)}$ and $\xi_{(G)}$ are gauge parameters. $\mathcal{L}_{GH}^{(A)}$, $\mathcal{L}_{GH}^{(B)}$ and $\mathcal{L}_{GH}^{(G)}$ denote ghost Lagrangians, respectively. $\Psi_{a}^{g}$ ($a=1,2,3,4$ and $g=1,2,3$) are the four $SO(5)$-vector (5 representation) fermions for each generation and $\Psi_{F_{i}}$ ($i=1,\cdots,N_{F}$) are the $N_{F}$ number of $SO(5)$-spinor (4 representation) fermions, which exist in the bulk space. The colour indices are not shown. The $SO(5)$ gauge fields $A_{M}$ are decomposed as

[TABLE]

where $T^{a_{L},a_{R}}(a_{L},a_{R}=1,2,3)$ and $T^{\hat{a}}(\hat{a}=1,2,3,4)$ are the generators of $SO(4)\simeq SU(2)_{L}\times SU(2)_{R}$ and $SO(5)/SO(4)$, respectively. For the fermion Lagrangian, $\Gamma^{M}$ denotes 5D gamma matrices which is defined by $\{\Gamma^{M},\Gamma^{N}\}=2\eta^{MN}$ ($\eta^{55}=+1$),

[TABLE]

and $\bar{\Psi}\equiv i\Psi^{\dagger}\Gamma^{0}$. $e_{A}{}^{M}$ is an inverse fielbein and $\Omega_{MBC}$ is the spin connection. $\epsilon(y)\equiv\sigma^{\prime}/k$ is the sign function and $c_{a}$ are the bulk mass parameters. The bulk mass parameters are set as $c_{1}^{g}=c_{2}^{g}$ and $c_{3}^{g}=c_{4}^{g}$.

The $SO(5)$-vector fermions are decomposed to the $SU(2)_{L}\times SU(2)_{R}$ bidoublet and singlet. The multiplet of the third generation are denoted as

[TABLE]

The first and second generation quarks and leptons are abbreviated.

The orbifold boundary conditions at $y_{0}=0$ and $y_{1}=L$ are given by

[TABLE]

By these boundary conditions, the $SO(5)\times U(1)_{X}$ symmetry is broken to $SO(4)\times U(1)_{X}\simeq SU(2)_{L}\times SU(2)_{R}\times U(1)_{X}$. For quarks and leptons, the bidoublets have the left-handed zero modes and the singlets have the right-handed zero modes.

The brane action is given by

[TABLE]

where

[TABLE]

are the right-handed fermions localised on $y=0$. $\kappa^{q}_{1,2,3}$, $\kappa^{\ell}_{1,2,3}$, $\tilde{\kappa}^{q}$ and $\tilde{\kappa}^{\ell}$ are the brane Yukawa couplings. The brane scaler breaks the $SU(2)_{R}\times U(1)_{X}$ symmetry to the $U(1)_{Y}$ symmetry spontaneously by acquiring the vacuum expectation value,

[TABLE]

The details of the brane interactions are shown in Ref. [41, 42]. The formulae of the KK expansion, the mass spectrum of the KK modes and the couplings are shown in Ref [48].

3 Parameters, couplings and decay widths

3.1 Bulk mass parameter dependence

I would like to emphasize the parameter dependence of fermion mode functions and fermion couplings with the KK gauge boson in this subsection. Generally, fermion with the bulk mass parameter $c$ on the warped metric is expanded as [61]

[TABLE]

where $N^{(n)}$ is the normalisation factor defined as

[TABLE]

and $f_{L,R}^{(n)}(y)$ is expressed as

[TABLE]

where $a_{n}$ and $b_{n}$ are constants determined by boundary conditions. The upper and lower sign of $\pm$ and $\mp$ correspond to the left- and right-handed fermions, respectively. Therefore the left-handed zero-mode is localised towards $y=0$ and $y=L$ for $\frac{1}{2}<c$ and $0<c<\frac{1}{2}$. In contrast, the right-handed zero-mode is localised towards $y=L$ for $0<c$. These behaviour of the left- and right-handed fermions are reversed for $c<0$.

From (13) and (14), it is straightforwardly derived that the product of the left- and right-handed fermion mode functions with the same KK number and the bulk mass, $f_{L}^{(n)}(y)\times f_{R}^{(n)}(y)$ is invariant under changing the sign of the bulk mass parameter $c\to-c$. Consequently, the Yukawa couplings and the fermion masses obtained from the Higgs boson vacuum expectation value are independent of the sign of the bulk mass parameters in the case that the left- and right- handed quarks have the same bulk mass. In this model, the above arguments are applicable.

Considering gauge boson, mode function of zero-mode is independent of $y$-coordinate and 1st KK gauge bosons have peaks near $y=L$ [62]. Therefore the right-handed fermions with $0<c$, the left-handed fermions with $c<0$ and the left-handed fermions with $-\frac{1}{2}<c<\frac{1}{2}$ couple to $Z^{\prime}$ bosons rather largely. As shown in the next subsection, the bulk mass parameter of the third generation quarks is $|c_{t}|<\frac{1}{2}$. Therefore the left-handed third generation quarks couplings with $Z^{\prime}$ bosons are large.

3.2 Parameters

The free parameters of this model are the warp factor $e^{kL}$ and $N_{F}$ which is $\Psi_{F}$’s degrees of the freedom, so once the $e^{kL}$ and $N_{F}$ are set, $\theta_{H}$ is determined. The physics of the quarks and leptons are almost independent of $N_{F}$ and determined by $\theta_{H}$ [44, 45]. The input parameters and the model parameters to realise the input parameters at the tree level are listed in Table 1 and Table 2 respectively, where $c_{e}$, $c_{\mu}$, $c_{\tau}$, $c_{u}$, $c_{c}$, $c_{t}$ and $c_{F}$ are the bulk mass parameters of leptons and quarks for the each generations and $\Psi_{F}$. As explained in previous subsection, the sign of the bulk mass parameter $c$ doesn’t affect the fermion mass in this model. Therefore only the absolute values of $c$ is determined from the input values. In the following, the bulk mass parameters of leptons and quarks are abbreviated as $c_{l}\equiv(c_{e},c_{\mu},c_{\tau})$ $c_{q}\equiv(c_{u},c_{c},c_{t})$. The resultant $W$-boson mass at the tree-level calculated by the boundary condition is $m_{W}^{\text{tree}}=80.0$ GeV. To realise the input parameters, the parameter region of $\theta_{H}$ is found to be $0.08\leq\theta_{H}\leq 0.10$. The lower limit of $\theta_{H}$ becomes slightly smaller for $N_{F}>4$. In $N_{F}=8$ case, lower limit of $\theta_{H}=0.078$.

3.3 Couplings and decay widths

The fermion couplings to $Z^{\prime}$s are listed in the Table 5, 5, 5, 8, 8 and 8. The 1st KK photon couplings to the left- and right-handed fermions for $c>0$ is equal to the left- and right-handed fermions for $c<0$ within 5 digits, respectively. The $Z_{R}^{(1)}$ coupling with ($u_{L}$, $d_{L}$) for $c<0$ is negative, although the couplings with ($c_{L}$, $s_{L}$) and ($t_{L}$, $b_{L}$) for $c<0$ are positive. This behaviour comes from the mass ratio of up-type and down-type quarks. In this model, left-handed top-quark is the mixing state of the ($t_{L}$, $B_{L}$, $U_{L}$) in (6) and their mixing ratio is

[TABLE]

where $\tilde{\kappa}/\kappa_{2}\simeq m_{b}/m_{t}$ in good accuracy. Their $SU(2)_{L}\times SU(2)_{R}$ isospins are

[TABLE]

Therefore up-type quark coupling with the $U(1)_{X}$ gauge boson is proportional to

[TABLE]

up-type quark coupling with the $SU(2)_{L}$ neutral gauge boson is proportional to

[TABLE]

and up-type quark coupling with the $SU(2)_{R}$ neutral gauge boson is proportional to

[TABLE]

For small $\theta_{H}$, the sign of the last one depends on the ratio of the up- and down-type quarks. Thus the $Z_{R}^{(1)}$ couplings with the up-quark have different sign with that of the charm- and top-quarks The behaviour of the $Z_{R}^{(1)}$ couplings with the down-type quarks are derived by the same reason.

The $Z^{\prime}$ masses obtained by the boundary conditions and the decay widths calculated from the couplings shown in Table 5, 5, 5, 8, 8 and 8 are summarized in Table 9. The 1st KK photon decay width is independent of the sign of the bulk fermion parameters.

4 Cross section and forward-backward asymmetry

The parameters are constrained by the experimental results of the forward-backward asymmetry at the $Z$-pole. For $c_{l}>0$, the deviations of the $Z$-boson couplings are $O(0.01)\%$. In contrast, for $c_{l}<0$, their deviations are $O(0.1)\%$. Thus the forward-backward asymmetry of $e^{+}e^{-}\to\mu^{+}\mu^{-}$ process at the $Z$-pole in the GHU model deviate nearly 10 % from the observed value for $c_{l}<0$. Consequently, the value $\sin^{2}\theta_{W}=0.23122$ is not valid and the value of $\theta_{W}$ which consistently explain the experimental results at the $Z$-pole must be searched. In this paper, $c_{l}<0$ case is not considered further.

The longitudinal polarisation $P_{e^{\pm}}$ ($-1\leq P_{e^{\pm}}\leq 1$) is introduced, where the electron and positron is purely right-handed when $P_{e^{\pm}}=1$. The cross section of $e^{-}e^{+}\to Z^{\prime}\to f\bar{f}$ at the centre-of-mass frame is given by

[TABLE]

where $\sigma_{LR}$ ($\sigma_{RL}$) is $e_{L}^{-}e_{R}^{+}(e_{R}^{-}e_{L}^{+})\to f\bar{f}$ cross section. The formula (20) is rewritten by using $P_{\text{eff}}=(P_{e^{-}}-P_{e^{+}})/(1-P_{e^{-}}P_{e^{+}})$ as $\sigma(P_{\text{eff}},0)=\sigma(P_{e^{-}},P_{e^{+}})/(1-P_{e^{-}}P_{e^{+}})$, then the ratio of $\sigma$ is parametrised by one polarisation parameter $P_{\text{eff}}$. Typical values of polarisation parameters are $(P_{e^{-}},P_{e^{+}})=(\pm 0.8,\mp 0.3)$ ($P_{\text{eff}}=\pm 0.887$).

Considering the $e^{+}e^{-}\to\mu^{+}\mu^{-}$ process, the difference of the cross sections between $c_{q}>0$ case ($\sigma^{c_{q}>0}$) and $c_{q}<0$ case ($\sigma^{c_{q}<0}$) arises from only the $Z^{\prime}$ decay widths. Consequently the deviation of $\sigma^{c_{q}<0}(\mu^{+}\mu^{-})$ from $\sigma^{c_{q}>0}(\mu^{+}\mu^{-})$ is small. As shown in  [49], at $\sqrt{s}=250$ GeV with 250 $\text{fb}^{-1}$ unpolarised beam, $4.66\times 10^{5}$ events are expected in the SM. Therefore the statistical uncertainty is 0.15 %. The difference of the cross sections of the two cases over the SM value, $(\sigma^{c_{q}>0}-\sigma^{c_{q}<0})/\sigma^{\text{SM}}(\mu^{+}\mu^{-})$ is less than 0.11% at $\sqrt{s}=250$ GeV with unpolarised beam. For the forward-backward asymmetry, $(A_{\text{FB}}^{c_{q}>0}-A_{\text{FB}}^{c_{q}<0})/A_{\text{FB}}^{\text{SM}}(\mu^{+}\mu^{-})$ is less than 0.04 % at $\sqrt{s}=250$ GeV with unpolarised beam. Thus the two cases are difficult to distinguish at the $e^{+}e^{-}\to\mu^{+}\mu^{-}$ process. The detailed analysis of the $e^{+}e^{-}\to\mu^{+}\mu^{-}$ process in the GHU model is shown in [49].

For the $e^{+}e^{-}\to\bar{q}q\leavevmode\nobreak\ (q\neq b,t)$ process, the cross section in the SM is $\sigma^{\text{SM}}(\bar{q}q)$= 9.75 pb and 7.26 pb at $\sqrt{s}=250$ GeV with unpolarised and polarised ($P_{e^{-}}=+0.8$ and $P_{e^{+}}=-0.3$) beams, respectively. For $\sqrt{s}=O(100)$ GeV region, the cross section in the GHU model with $c_{q}>0$ are smaller than that in the SM. In contrast, that in the GHU model with $c_{q}<0$ are larger than that in the SM as shown in figure 1. In figure 2, the ratio of cross sections in the GHU model to that in the SM with polarised beams are plotted. At $\sqrt{s}=250$ GeV with $P_{\text{eff}}=+0.887$ and 250 fb*-1* luminosity, the event number of the SM is $1.814\times 10^{6}$ and the statistical uncertainty $\sigma$ over the event number is $0.074\%$. $\sigma(\bar{q}q)$ in the GHU model deviates from the that in the SM larger than 3 % for $c_{q}>0$ at $\sqrt{s}=250$ GeV with $P_{\text{eff}}=+0.887$ and 250 fb*-1* luminosity. The deviations at $\sqrt{s}=250$ GeV are summarised in Table 10.

For the $\bar{c}c$ final state, the forward-backward asymmetry is also measured. The $A_{\text{FB}}^{\text{GHU}}(\bar{c}c)$ decrease from the $A_{\text{FB}}^{\text{SM}}(\bar{c}c)$ for $c_{q}>0$. In contrast, $A_{\text{FB}}^{\text{GHU}}(\bar{c}c)$ increase from the $A_{\text{FB}}^{\text{SM}}(\bar{c}c)$ for $c_{q}<0$. At $\sqrt{s}=250$ GeV with $P_{\text{eff}}=+0.887$ polarised beam, $A_{\text{FB}}^{\text{SM}}(\bar{c}c)=0.700$ and with 250 $\text{fb}^{-1}$ luminosity, the statistical uncertainty is 0.16%. The deviations of the $A_{\text{FB}}^{\text{GHU}}(\bar{c}c)$ from $A_{\text{FB}}^{\text{SM}}(\bar{c}c)$ are $-0.38\%$, $-0.31\%$ and $-0.25\%$ for $c_{q}>0$ and $\theta_{H}=0.10$, $0.09$ and $0.08$ at $\sqrt{s}=250$ GeV with $P_{\text{eff}}=+0.887$. For $c_{q}<0$ and $\theta_{H}=0.10$, $0.09$ and $0.08$ at $\sqrt{s}=250$ GeV with $P_{\text{eff}}=+0.887$, the deviations are $+1.68\%$, $+1.25\%$ and $+1.03\%$. The polarisation dependence of the ratio $A_{\text{FB}}^{\text{GHU}}/A_{\text{FB}}^{\text{SM}}(\bar{c}c)$ are plotted in figure 3 and the deviations at $\sqrt{s}=250$ GeV, 500 GeV and 1 TeV are summarised in Table 11, Table 12 and Table 13, respectively.

For the $e^{+}e^{-}\to\bar{b}b$ process, the cross section in the SM is $\sigma^{\text{SM}}(\bar{b}b)$= 1.77 pb and 1.02 pb at $\sqrt{s}=250$ GeV with unpolarised and polarised ($P_{e^{-}}=+0.8$ and $P_{e^{-}}=-0.3$) beams, respectively. The statistical uncertainty at $\sqrt{s}=250$ GeV and 250 fb*-1* luminosity with $P_{e^{-}}=+0.8$ and $P_{e^{+}}=-0.3$ beam are 0.20 %. For this process, the cross sections in the GHU model with $c_{q}>0$ and $c_{q}<0$ cases both decrease from that in the SM. The deviations at $\sqrt{s}=250$ GeV are summarised in Table 10. The deference between the $c_{q}>0$ and $c_{q}<0$ cases more obviously appear at the forward-backward asymmetry. In the SM, $A_{\text{FB}}^{\text{SM}}(\bar{b}b)=0.618$ and $0.366$ at $\sqrt{s}=250$ GeV with unpolarised and $P_{\text{eff}}=+0.887$ beams, respectively. In the GHU model the $A_{\text{FB}}(\bar{b}b)$ increase from the SM value, $A_{\text{FB}}^{c_{q}>0}(\bar{b}b)$ increase $4.24\%$, $4.07\%$, $4.15\%$ and $A_{\text{FB}}^{c_{q}<0}(\bar{b}b)$ increase $7.33\%$, $5.69\%$, $3.83\%$ at $\sqrt{s}=250$ GeV with $P_{\text{eff}}=+0.887$ for $\theta_{H}=0.10$, $0.09$ and $0.08$. In figure 4, the ratio of the $A_{\text{FB}}(\bar{b}b)$ in the GHU model to that in the SM with polarised beams are plotted. The $c_{q}<0$ case predicts larger deviation than the $c_{q}>0$ case for $\theta_{H}=0.10$ and $0.09$. At $\sqrt{s}=250$ GeV with $P_{\text{eff}}=+0.887$ and 250 fb*-1* luminosity, the statistical uncertainty of the $A_{\text{FB}}(\bar{b}b)$ in the SM is 0.70%. $A_{\text{FB}}(\bar{b}b)$ in the GHU model deviates from the that in the SM larger than $5.4\sigma$. The deviations at $\sqrt{s}=250$ GeV, 500 GeV and 1 TeV are summarised in Table 11, Table 12 and Table 13, respectively.

The forward-backward asymmetry of the $e^{+}e^{-}\to\bar{t}t$ process is measured at $\sqrt{s}=500$ GeV. The polarisation dependence of $A_{\text{FB}}^{\text{GHU}}/A_{\text{FB}}^{\text{SM}}(\bar{t}t)$ is qualitatively similar to that of $A_{\text{FB}}^{\text{GHU}}/A_{\text{FB}}^{\text{SM}}(\bar{b}b)$, which also increase from that in the SM. $A_{\text{FB}}^{c_{q}>0}(\bar{t}t)/A_{\text{FB}}^{\text{SM}}(\bar{t}t)-1=5.36\%$, $5.03\%$, $5.14\%$ and $A_{\text{FB}}^{c_{q}<0}(\bar{t}t)/A_{\text{FB}}^{\text{SM}}(\bar{t}t)-1=9.25\%$, $7.23\%$, $4.91\%$ at $\sqrt{s}=500$ GeV with $P_{\text{eff}}=+0.887$ for $\theta_{H}=0.10$, $0.09$ and $0.08$. At $\sqrt{s}=500$ GeV with $P_{\text{eff}}=+0.887$ and 500 fb*-1* luminosity, $\sigma^{\text{SM}}(\bar{t}t)=479$ $\text{fb}^{-1}$, $A_{\text{FB}}^{\text{SM}}(\bar{t}t)=0.463$ and the uncertainty of the $A_{\text{FB}}^{\text{SM}}(\bar{t}t)$ is 0.538 %. The deviations of the forward-backward asymmetries of the $e^{+}e^{-}\to\bar{c}c,\ \bar{b}b,\ \bar{t}t$ processes at $\sqrt{s}=500$ GeV and 1 TeV are summarised in Table 12 and Table 13.

5 Summary

In the above calculations, the quark bulk mass parameters $(c_{u},c_{c},c_{t})$ are assumed to be all positive or all negative. It is also allowed to be only one of them is positive or negative. In the case, the $Z^{\prime}$ decay widths change from the values shown in Table 9, therefore the cross sections and the forward-backward asymmetries slightly change from the results in this paper. Neglecting the difference arising from the $Z^{\prime}$ decay widths, the sign of the $c_{c}$ is determined by measuring the $A_{\text{FB}}(\bar{c}c)$ and the sign of the $c_{u}$ is determined by the $A_{\text{FB}}(\bar{c}c)$ and $\sigma(\bar{q}q)$ at $\sqrt{s}=250$ GeV. It is difficult to determine the sign of the $c_{t}$ by measuring $A_{\text{FB}}(\bar{b}b)$ and $A_{\text{FB}}(\bar{t}t)$. At the ILC 250 GeV, the $c_{q}<0$ case predicts $4\sigma$ larger deviation of the $A_{\text{FB}}(\bar{b}b)$ than the $c_{q}<0$ case for $\theta_{H}=0.10$, and at the ILC 500 GeV $5\sigma$ larger deviation for $\theta_{H}=0.09$. For $\theta_{H}=0.08$, to clearly determine the sign of $c_{t}$ by observing the $A_{\text{FB}}(\bar{t}t)$, higher energy and luminosity such as the ILC 1 TeV are necessary.

In this paper, the forward-backward asymmetry of the $e^{+}e^{-}\to\bar{q}q$ processes in the GHU model are studied for two cases where all of the quark bulk mass parameters are positive or negative. The GHU model predicts large deviations at the $\sqrt{s}=250$ GeV with polarised beams. Therefore the GHU model is testable at the ILC 250 GeV. The signs of the bulk mass parameters are distinguished at the ILC 500 GeV or ILC 1 TeV. For the case where the lepton bulk mass parameters are negative, detail is going to be analysed in near future.

Acknowledgements

I thank Yutaka Hosotani, Hisaki Hatanaka, Yuta Orikasa and Naoki Yamatsu for many advices and important discussions. This work is supported by the National Natural Science Foundation of China (Grant Nos. 11775092, 11675061, 11521064 and 11435003), and the International Postdoctoral Exchange Fellowship Program (IPEFP).