X(16.7) production in electron–positron collision

The anomaly found in the excited 8Be\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^8\text {Be}$$\end{document} nuclear transition to its ground state is attributed to a spin-1 gauge boson X(16.7). To hunt for this boson, we propose two traps: e+e-→Xγ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e^+e^-\rightarrow X\gamma $$\end{document} and J/ψ→Xγ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J/\psi \rightarrow X\gamma $$\end{document}, both following with X→e+e-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$X\rightarrow e^+e^-$$\end{document} decay. We adopt the “vector minus axial-vector” interaction hypothesis. Analysis on the X(16.7) decay length, production rates, differential distribution with respect to the e+e-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e^+e^-$$\end{document} invariant-mass spectrum, and the signal-to-noise ratios (SNR) after the smearing at BESIII detector are discussed in detail. Given the coupling strength of X to vector/axial-vector currents gfv/a∼10-3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$g_f^{v/a}\sim 10^{-3}$$\end{document} at BESIII: (1) there would be about 6000 X measurable events per year in electron–positron collision, yet with a large background after smearing; (2) while in J/ψ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J/\psi $$\end{document} decays, we find that the axial-vector current may come into play; though merely 52 events may appear, the SNR are inspiring even after smearing.


Introduction
So far, the Standard Model (SM) has achieved a great success in experiencing numerous experimental tests [1]. And the search of new physics beyond SM is now a major activity for both experimental and theoretical physicists. Gauge bosons, γ, Z 0 , W ± and gluons, as the messengers of the SM interaction forces, play crucial role in revealing the nature of interactions. A new gauge boson (or the fifth force), then becomes the prey which we all are hunting for.
Finding out whether there exists a new type of force beyond SM is a very tempting task. The ATLAS and CMS Collaborations at LHC searched for the TeV Z boson using the Run 2 data [2][3][4][5][6]. For B-factories, the BaBar Collaboration explored the 0.02-10.2 GeV region using the e + e − a e-mail: jiangjun13b@mails.ucas.ac.cn b e-mail: chenglogbin10@mails.ucas.ac.cn c e-mail: alvan@ucas.ac.cn d e-mail: qiaocf@ucas.ac.cn invariant-mass spectrum and upper limits on mixing strength between the dark photon (U) and SM photon (γ ) were placed [7]. At BESIII detector, the 1.5-3.4 GeV mass range was explored and limits on the U−γ mixing strength were also set [8]. Similarly in the KLOE-2 experiment, they searched for the dark photon in both e + e − → Uγ [9] and φ → Uη processes [10]. The NA48/2 experiment at CERN also searched for the dark photon in π 0 decays [11]. While in HADES experiment, the dark photon search was carried out using e + e − spectrum in the p-p, N-b reactions, as well as the Ar + KCl reaction [12]. These new and dark bosons, Z and U, were also probed with the precise electroweak data [13], in neutrino-electron scattering experiments [14], planned and future experiments [15][16][17][18], etc., and have attracted the bright physicists [19][20][21].
In 2016, an extraordinary experimental phenomenon was observed in the isoscalar 8 Be nuclear transition, 8 Be * → 8 Be [22], and new measurements are presented three times recently [23][24][25]. A significant enhancement relative to the internal pair creation was observed at large angles in the invariant-mass distribution of electron-positron pairs production. However, no anomaly is seen in the isovector 8 Be * → 8 Be transition. This observation is hard to be understood within the regime of conventional theory, but could be attributed to a neutral isoscalar particle X with the mass of 16.7 MeV and the saturating decay X → e + e − beyond SM.
At the beginning, this anomaly has been interpreted as a new vector boson which mediates a weak fifth force beyond SM [26], and a realistic model for the fifth force is also proposed [27]. Other possible explanations on this new boson, light pseudoscalar boson, protophobic vector boson, axialvector particle, etc., are also widely studied [28][29][30][31][32][33][34][35]. Among them, the axial-vector particle proposal is suggested in Ref.s [26,29] but not pursued systematically. In Ref. [30], the production of a hidden vector boson with axial-vector couplings to leptons and light quarks in the isoscalar 8 Be * → 8 Be nuclear transition is investigated. Note, along with the axialvector couplings to the Standard Model fermions, we need to devote some effort to obtain a UV-complete anomaly-free theory [36,37]. Besides, models including this X boson as a mediator to the dark sector or giving constrains on dark matters are discussed [38][39][40][41]. This boson is also introduced to account for some other anomaly observations [42,43] or as the light gauge boson mediator in the rare kaon and pion decays [44,45]. Meanwhile, suggestions for future experimental research are proposed [46][47][48][49][50].
Scientifically, to investigate further in experiments and get more knowledge of this-yet-not-independently-verified particle is currently the most important task among all the studies. In this work, we estimate the production of this X boson associated with a photon in electron-positron collision, as well as in the J/ψ decays. The complete "vector minus axialvector (V-A)" interaction is considered seriously, and we find that the axial-vector current contribution plays an important role, especially in J/ψ decays. Various analysis on its production are presented in detail, particularly the invariantmass spectrum of the electron-positron pairs in final states, and the signal-to-noise ratios after the "smearing".

Traps arranged for X (16.7)
To be a new member of the Particle Zoo, X (16.7) needs to be tagged. Since the proposals of scalar bosons are excluded [26], we here focus on the spin-1 gauge boson hypothesis. And the Lagrangian added to the Standard Model one can be formulated as Here, e is the electron charge and g v/a f denote the coupling strength of X to vector/axial-vector currents, which means X boson can either be a massive γ -like particle or a Z 0 -like one. For the g v e , it has been constrained by experimental data to the region of 2 × 10 −4 ≤ |g v e | ≤ 1.4 × 10 −3 [26]. In this paper, we will take the "vector minus axial-vector (V-A)" interaction vertex as a general situation, one can obtain the vector one readily by taking the coupling parameter g a f = 0.
Since the mass of X boson m X is much smaller than the energy of usual electron-positron colliders, its production always associates with another gauge boson radiation, i.e. the photon γ , as shown in Fig. 1.
The differential cross-section of the process can be readily obtained, In which θ is the emitting angle of photon with respect to the e + e − beam axis. Taking the fine structure constant α = 1/137 and center-of-mass energy (CMS) √ s = 3.7 GeV, we can obtain the differential cross-section as being displayed in Fig. 2 (left). One may notice from the figure that the main contribution comes from the region where | cos θ | is large. Taking the high-energy limit ( √ s m X , m e ), the differential cross-section turns to which agrees with the result for e + e − → 2γ by taking g v e = 1, g a e = 0 [51]. Integrating over cos θ of the differential cross-section Eq.
(2), one can get the colliding energy dependence of the crosssection, which is presented in Fig. 2 (right). Where two CMS energies ( √ s = 3.7, 10.6 GeV) are highlighted. The first one is the typical CMS energy of BESIII detector, which has the luminosity of 10 33 cm −2 s −1 at √ s = 3.7 GeV. While the second one is that of B-factories, i.e. the BaBar and Belle Collaborations. We find that the cross-section drops by about one order of magnitude when CMS energy increases from 3.7 to 10.6 GeV. Another important conclusion drawn from Fig. 2 is that the contribution of axial-vector current ("A only") and the vector current part ("V only") have almost the same share when taking g v e = g a e . However, the "A only" contribution has not been taken into account in the literature seriously.
At the CMS energy of √ s = 3.7 GeV, and taking g v/a e = 10 −4 -10 −3 , we find the cross-section of e + e − → X + γ lies The legend "V-A" refers to the complete "vector minus axial-vector" contributions, while "V only" stands for the vector case only in the region of 0.0064-0.64 pb, which happens to be two times the amount of the "V only" result. Since the BESIII detector can cover about 93% of the 4π solid angle, and with the luminosity of L-10 4 pb −1 year −1 , there will be about 60-6000 X bosons to be produced per year and more in yet collected data. 1 At Babar, with the copious 514 fb −1 data [7] and about 90% solid angle coverage [53], the X boson events would be around 4.1 × (10 2 -10 4 ).
Experimentally, the X boson would be reconstructed by its decay products. A priori the X boson can decay into e + e − , νν, 3γ or unknown particles, and νν and 3γ decay modes are highly suppressed [26]. So we assume that the X boson decays to electron-positron pairs in saturation. And the decay width reads which is consistent with the "V only" result by taking g a e = 0 [54]. Since m 2 X m 2 e , the decay width is about two times the amount of that in "V only" situation when taking g v e = g a e . Given the order of magnitude of g v/a e as 10 −4 -10 −3 , the X boson decay width varies as 8.1 × (10 −4 -10 −2 ) eV, which corresponds to the lifetime τ = 8.1 × (10 −15 -10 −13 ) s at 1 Roughly, a collider operates only about 10 7 s in 1 year [52]. So it is customary to estimate that 10 33 cm −2 s −1 10 4 pb −1 year −1 . Moreover, when evaluating the events, we simply multiply the cross-section with the integrated luminosity and the percentage of the solid angle coverage. the X boson rest frame. While in the experiment frame, the velocity of X boson and the energy of the emitting photon s respectively. After performing the Lorentz boost, the lifetime of X boson could increase by about two orders of magnitude, i.e., 9.0×(10 −13 -10 −11 ) s. Hence, the decay length in the experiment frame would be where g v e = g a e is adopted and symbolized as g e ,h and c are the reduced Planck constant and velocity of light respectively. Then numerically we have 0.27 mm < L < 27 mm at √ s = 3.7 GeV, which is measurable at BESIII. At the B-factories, whose CMS energy is about 10.6 GeV, the decay length can reach 0.77-77 mm. While at CLOE with √ s = 1.0195 GeV, the decay length of X boson is about 10% of that at Bfactories. Evidently, it is more attainable to measure the decay length of X boson at BaBar/Belle than at BESIII or CLOE. The measurement of the decay length is meaningful not only for the aim of disentangle signals from the background, but also for the determination of the coupling strength g v/a e . It's time to set the trap with the X → e + e − decay mode. And the dominant background originates from the virtual photon propagated scattering process, e + e − γ * − → e + e − γ . We consider Feynman diagrams that contain both X and γ * as inner propagators, which have been displayed in Fig. 3. Fig. 3 The Feynman diagrmas for the signal/background processes e + e − X/γ * − −− → e + e − γ Fig. 4 The differential cross-section with respect to the invariantmass distribution of e + e − pairs for the e + e − X/γ * − −− → e + e − γ processes, where we have adopted the good tracks conditions for both electron/positron and photon. Left: "Total" means it contains the con-tribution of cross-terms between the X -propagated Feynman diagrams and the γ * -propagated ones. Right: the contribution of cross-terms is excluded Here, we would evaluate the differential distribution with respect to the invariant-mass spectrum of electron-positron pairs in final states under the BESIII experiment conditions as a fine example. At BESIII, the good charged tracks are constrained in the region of | cos α| < 0.93, while the photon selection condition is | cos β| < 0.8 with the energy E > 25 MeV for the barrel [55], 2 here α/β are the polar angles with respect to the e + e − beam axis. In Fig. 4 we show the invariant-mass ( √ s ee or M ee ) distribution of e + e − pairs in final states for the differential cross-section, where the constrains on α/β are adopted. Given g v/a e = 10 −3 , we present the contribution of complete vector minus axialvector current ("V-A") and the vector current ("V only") one separately. Here, the legend "Background" refers to the contribution of e + e − γ * − → e + e − γ process; "Signal" stands for the contribution of e + e − X − → e + e − γ only; and the "Total" contains both the previous two parts and those from the crossterms between them. By comparing the "Total" and "Sig-nal+Background" figures in the left and right diagrams of Fig. 4, one can easily obtain the running line-shape of the cross-terms. From Fig. 4, we can see that the X boson acts as a rather sharp peak above the background around the mass of 16.7 MeV, which tells that a precise measurement on the invariantmass spectrum of e + e − pairs in final states will greatly suppress the background. Then, one may wonder if the energy resolution (δ E ) of the BESIII detector is sufficient to undertake this invariant-mass analysis, i.e., the capability of identifying the signals over the background after "smearing". In the following, we take the sample of 6000 X events as an example to illustrate the smearing and estimate the signalto-noise ratios after smearing at BESIII.
Theoretically, these 6000 X boson events would be reconstructed precisely at M ee = 16.7 MeV in the e + e − invariantmass spectrum, as is shown in Fig. 5a. However, experimentally the detector has the restricted energy resolution. So, some X events would not be found at M ee = 16.7 MeV, but The X resonance bump of "Signal V −A " before and after smearing simulated with 6000 X events. The b-d diagrams correspond to the energy resolution of δ E = 2, 4, 5 MeV after smearing respectively positioned at somewhere deviated from it. And the deviation amounts comply with a Gaussian distribution with the expectation of zero and the standard deviation of the energy resolution (δ E ). In Fig. 5b-d, we present how the 6000 X events are distributed in a wider invariant-mass range around M ee = 16.7 MeV in the spectrum after smearing, where δ E = 2, 4, 5 MeV accordingly. As is known, when the resonance peak is smeared to a wider M ee range, more background events would be included. To make it clear, we have considered four ranges of the e + e − invariant-mass spectrum M ee around the X resonance, and the corresponding signal-to-noise ratios (SNR) are presented in Table 1. Note, for the signal events (S), we adopt the rough estimate of X boson events produced in e + e − → X γ , since the following decay X → e + e − is in saturation; while the background events (B) are obtained within the stringent constrains of the good tracks conditions. In Table 1, we find that the SNR is quite encouraging in the first M ee range (near the X boson peak). And the last three mass ranges are roughly corresponding to those with the energy resolution of δ E = 2, 4, 5 MeV (Fig. 5b-d) respectively, where the background is much more noisy. With the decay length of X boson to further increase the SNR, we believe that these results are inspiring in the future X boson hunting if g v/a e -10 −3 . In fact, the BaBar, BESIII and KLOE experiments had ever searched for the dark photon in e + e − collision [7][8][9]. However, the mass range in BaBar experiment was in between 20 MeV-10.2 GeV, and BESIII searched for the dark photon in the 1.5-3.4 GeV mass region, both overshot the 16.7 MeV X boson. The CLOE experiment explored the 5-520 MeV range, yet did not find clear signatures of the dark photon, whereas it constrained the parameter g v e to be |g v e | ≤ 2×10 −3 , which is in accordance with what we employ in this work.
Note that, other proposals of searching for an extra U(1) gauge boson U with its decays into e + e − or invisible at electron-positron colliders are also provided in Ref. [47,49], as well as the background analysis [50], yet with different methods/focus in comparison with the present work.
Up to date the J/ψ events collected at BESIII are (1310.6 ± 7.0) × 10 6 [56], which is quite abundant. Hence in this subsection, we will consider the X (16.7) boson production in J/ψ decays seriously. Table 1 The signal-to-noise ratios (SNR) around M ee = 16.7 MeV, where we adopt four invariant-mass ranges M ee for the background process. The last three ranges roughly correspond to the energy res-olution of δ E = 2, 4, 5 MeV (Fig. 5b-d)  In this subsection, we continue to take the Lagrangian terms as shown in Eq. 1, i.e. the interaction vertex is still the "vector minus axial-vector (V-A)" type. The Feynman diagrams can be easily obtained from Fig. 6c, d, where the X boson decays into the e + e − pairs. Then the decay width reads where the squared wave function at the origin 2 = (2m c ) 2 (J/ψ→e + e − ) 16π e 2 c α 2 [57], with e c = 2/3 and (J/ψ → e + e − ) = 5.55 × 10 −6 GeV [1]. Obviously, the contribution of vector current ("V only") vanishes. 3 Taking the total decay width of J/ψ as 92.9 keV, then we obtain the branching fraction, 0.0398×(g a c ) 2 . Given the 1.3×10 9 J/ψ events and g a c = 10 −4 -10 −3 , the expected events in J/ψ → X γ would be 1(0.52)-52 (values in brackets are the estimate before rounding up).
The X boson signals will also be captured through X → e + e − decay in experiments. Feynman diagrams of the related processes J/ψ X/γ * − −− → e + e − γ are presented in Fig. 6. One may notice that if X boson was a massive γ -like particle, there would be no contributions from the last two figures of Fig.  6c, d. And we also find that the contribution of axial-vector current ("A only") from Fig. 6a, b vanishes at tree amplitudes level. According to the numerical evaluation (assuming g v/a c = g v/a e = 10 −3 ), it is found that contribution of vector current ("V only") of Fig. 6a, b is suppressed by a factor of 10 −12 in comparison with that of the background process J/ψ γ * − → e + e − γ . So it can be concluded that the main contributions to the decay width of J/ψ X − → e + e − γ come from the axial-vector current interaction in Fig. 6c, d. The invariant-mass distribution of electron-positron pairs in final states of processes J/ψ X/γ * − −− → e + e − γ has been presented in Fig. 7 (left), where we adopt no cuts on emitting polar angles of the final particles. Obviously, the "Signal V −A " line shape overlaps with the "Total" one, which implies that the contribution of cross-terms between the Xpropagated Feynman diagrams and the γ * -propagated ones is negligible. And we find that around the X resonance, the differential decay widths of the signals are much larger than the ones of background, i.e., the signal to noise ratio is dramatic before smearing. Here, we evaluate the contribution from the axial-vector current purposely (the "Signal A " dashed line), i.e., the contribution of the axial-vector current in Fig. 6c, d.
And it is about one half of the contribution of "Signal V −A ", out of our expectation discussed in above paragraph. We find that the other half also comes from the Fig. 6c, d, which is the contribution of the cross-term between the axial-vector term ig a c γ μ γ 5 in X − cc vertex and the vector term ig v e γ μ in X − e + e − vertex, which has been excluded in the "Signal A " case.
In Fig. 7 (right), we present the distribution of the estimated 52 X boson events in the M ee spectrum after the smearing with the energy resolution of δ E = 4 MeV. In Table 2, we have considered several M ee ranges around X resonance, and the related signal-to-noise ratios (SNR) are presented. Note that we adopt no cuts on the emitting polar angles of the final particles here. And the third range is roughly corresponding to the energy resolution of δ E = 4 MeV (Fig. 7, right). Given the coupling strength ∼ 10 −3 , although we have only 52 X events, the amount of the background events are even fewer comparatively, which implies that there might be some good news of trapping the X boson in J/ψ → e + e − γ process.
Note that, in contrast with the vector current coupling, such X boson accounting for the isoscalar 8   transition anomaly with its couplings to quarks and leptons through an axial-vector current provides a natural suppression of the isovector 8 Be * → 8 Be transition [30]. Moreover, a theory with axial-vector couplings motivated by several MeV-scale anomalies can be UV-completed consistent with the Standard Model gauge invariance, see Refs. [36,37] for details. In Ref. [30], according to their simple UV-complete model, the axial-vector couplings to u, d quarks were constrained to be at the order of 10 −5 -10 −4 , but this still leaves room to have a large enough coupling to explain the 8 Be anomaly, say if the assumptions about the couplings they adopted were relaxed as mentioned there, or in other consistent UV models. In the minimal flavor violating limit, one may assume the charm axial-vector coupling to be same as the up quark one. In this case, and consider the X − cc axial-vector coupling to be at the order of 10 −4 , the expected X (16.7) events in J/ψ decays would be about 1 as estimated above. This implies that the J/ψ search for the X (16.7) boson at BESIII will not reach the sensitivity below 10 −4 that Ref. [30] suggests, with the 1.3 × 10 9 J/ψ events collected to date.

Summary and prospect
In summary, inspired by the 6.8σ anomaly in 8 Be nuclear transition experiment and the passion of searching for the new gauge boson, we investigate the possibility of hunt-ing for this yet-not-verified gauge boson X (16.7) in both e + e − collision and J/ψ decays. In our model, the complete "vector minus axial-vector" interaction vertex is taken into account. We have set two traps for the X (16.7) hunting: e + e − → X γ and J/ψ → X γ , both following with the saturated X → e + e − decay. Analytical results of differential distribution in e + e − → X γ and the decay widths of both J/ψ → X γ and X → e + e − are presented and discussed. Phenomenologically, we evaluate the decay length of X boson in the experiment frame, the production rates, the cross-section/decay-width differential distribution with respect to the e + e − invariant-mass spectrum, and particularly the signal-to-noise ratios, which may be increased by the X boson decay length, before and after the smearing for both processes.
For X boson production in e + e − collision at √ s = 3.7 GeV, given the coupling strength g v/a e = 10 −4 ∼ 10 −3 , it is found that the decay length is 0.27 mm < L < 27 mm, the expected X boson events are 60 ∼ 6000 per year, and the signal-to-noise ratios (SNR) decreases significantly after smearing. We also find that contributions of vector and axialvector currents are of the equal importance. While in J/ψ decays, we find that the axial-vector current dominates the X (16.7) production associated with a photon. Given copious 1.3 × 10 9 J/ψ events at BESIII and g a c = 10 −3 , there would be only 52 X events within the data, while the SNR are inspiring after smearing. That is, though e + e − → X γ process may Table 2 The signal-to-noise ratios (SNR) around the X resonance mass (16.7 MeV), where we adopt four invariant-mass ranges for the background process. Note, signal events would be suppressed by two orders of magnitude when taking g a c = 10 −4 yield more signals, it has a relatively larger background in regard to J/ψ → X γ process. In all, by measuring the final state e + e − invariant-mass spectrum, we find it is possible for BESIII experiment to perform a decisive measurement on the X (16.7) boson.