Production of double heavy quarkonia at super Z factory
Authors: Liao Qi-Li , Jiang Jun , Zhao Yu-Han 
Published: January 2023

Within the color singlet model, we calculate the exclusive production of double charmonia, double bottomonia, and double $$B_c$$

B
c

mesons at future super Z factory. The two heavy quarkonia or $$B_c$$

B
c

’s are either two S-wave Fock states ($$^1S_0, ~^3S_1$$

1

S
0

,

3

S
1

), or one S-wave and one P-wave states ($$^1P_1,~^3P_J~(J=0,1,2)$$

1

P
1

,

3

P
J

(
J
=
0
,
1
,
2
)

). The top three $$Z^0$$

Z
0

propagated channels in cross sections for double charmonia are $$\sigma {(J/\psi +h_c)}_{Z^0}=(454.0^{+39.1}_{-59.2})\times 10^{-5} fb, \sigma {(\eta _c+\chi _{c2})}_{Z^0}=(284.2^{+22.2}_{-35.3})\times 10^{-5} fb$$

σ

(
J
/
ψ
+

h
c

)

Z
0

=

(
454
.

0

-
59.2

+
39.1

)

×

10

-
5

f
b
,
σ

(

η
c

+

χ

c
2

)

Z
0

=

(
284
.

2

-
35.3

+
22.2

)

×

10

-
5

f
b

, and $$\sigma {(\eta _c+\chi _{c0})}_{Z^0}=(137.0^{+9.9}_{-16.4})\times 10^{-5} fb$$

σ

(

η
c

+

χ

c
0

)

Z
0

=

(
137
.

0

-
16.4

+
9.9

)

×

10

-
5

f
b

. For double bottomonia, they are $$\sigma {(\eta _b+\Upsilon )}_{Z^0}=(428.6^{+67.9}_{-60.9})\times 10^{-4} fb, \sigma {(\Upsilon +\chi _{b2})}_{Z^0}=(230.6^{+23.6}_{-22.2})\times 10^{-4} fb$$

σ

(

η
b

+
Υ
)

Z
0

=

(
428
.

6

-
60.9

+
67.9

)

×

10

-
4

f
b
,
σ

(
Υ
+

χ

b
2

)

Z
0

=

(
230
.

6

-
22.2

+
23.6

)

×

10

-
4

f
b

, and $$\sigma {(\Upsilon +\Upsilon )}_{Z^0}=(119.9^{+18.4}_{-16.6})\times 10^{-4} fb$$

σ

(
Υ
+
Υ
)

Z
0

=

(
119
.

9

-
16.6

+
18.4

)

×

10

-
4

f
b

. For double $$B_c$$

B
c

mesons, they are $$\sigma {(B_c^{*+}+B_c^{*-})}_{Z^0}=(1150^{+87}_{-84})\times 10^{-3} fb$$

σ

(

B
c

∗
+

+

B
c

∗
-

)

Z
0

=

(

1150

-
84

+
87

)

×

10

-
3

f
b

, $$\sigma {(B_c^{*+}+\chi _{bc2}^-)}_{Z^0}=(1133^{-16}_{+22})\times 10^{-3} fb$$

σ

(

B
c

∗
+

+

χ

b
c
2

-

)

Z
0

=

(

1133

+
22

-
16

)

×

10

-
3

f
b

, and $$\sigma {(\eta _{bc}^++B_c^{*-})}_{Z^0}=(634.6^{+47.8}_{-45.8})\times 10^{-3} fb$$

σ

(

η

bc

+

+

B
c

∗
-

)

Z
0

=

(
634
.

6

-
45.8

+
47.8

)

×

10

-
3

f
b

. Here the uncertainties come from the varying masses of constituent heavy quarks, which bring up to 20% corrections. The cross sections of double $$B_c$$

B
c

mesons are roughly one order of magnitude larger than those of the double bottomonia, and two orders of magnitude larger than those of the double charmonia. To make it helpful for experimental study, we present the total cross sections $$\sigma $$
σ
as functions of CM energy $$\sqrt{s}$$

s

, $$\sigma $$
σ
as functions of the renormalization scale $$\mu $$
μ
, the angle distributions $$d\sigma /dcos\theta $$

d
σ
/
d
c
o
s
θ

, and the $$p_T$$

p
T

distributions $$d\sigma /dp_{t}$$

d
σ
/
d

p
t

. We also find that the initial state radiation can bring about 30–40% suppresions when 1%$$m_Z$$

m
Z

energy is losing, and cross sections can increase by about 2–3 times or decrease by an order of magnitude when adopting different potential models which becomes the major source of uncertainty. The numerical results show that it might be not optimistic for the experimental observation, but it is still far from excluded at the FCC-ee and also the CEPC running in the Z factory mode.