Arxiv:1910.10459V1 [Nucl-Ex]

Home , Axion, X17 particle

arXiv:1910.10459v1 [nucl-ex] 23 Oct 2019 onqak,witna utpe ftepsto charge 2 positron relation the the of satisfy multiples e, as written quarks, down hog ih suoclrpril 9.Gvnthe Given [9]. particle results pseudoscalar experimental light the a of through explanation possible other [8]. model doublet Higgs two includ- extended performed an been have ing models such of studies many erhsfor searches epoohbc .. h charges the i.e., protophobic, be the of couplings the require o stu rdcdi h ea fa xie tt to state excited an bo- of decay X17 state, the ground The in the cou- produced particles. thus some is Model(SM) with son force Standard fundamental to fifth pling a mediate may hog h X17 the through iha1. e,J MeV, 16.7 a with h nml,adfidterqie omfco ob un- be to factor of the form origin for required possible realistic the a find and as anomaly, factor the form explored transition They nuclear physics. the nuclear within anomaly the plain 3]. [2, community physics the hadron in interest and experimental particle and theoretical enhanced an eaino h oo otedr atrpolmand problem (g-2) the matter explain dark might the it to that fact boson the X the of relation n ea fa nemdaepril 1 ihmass with X17 particle intermediate 18.15 creation an the the m of of as interpreted decay correlation was and angular This the [1]. transition observed in MeV was peak-like angles creation Significant large pair at internal [1]. the of observed enhancement was correlation lar J eain o h 76MV n 81 e,J MeV, 18.15 and MeV, 17.6 the for relations π X ttesm ie lwne n oet aean- made Moretti and Ellwanger time, same the At osrit nsc e atce oal from notably particle, new a such on Constraints h aawr xlie yFn n owres[,6] [5, co-workers and Feng by explained were data The hn n ilr[]ivsiae h osblt oex- to possibility the investigated [4] Miller and Zhang eety emaue lcrnpsto nua cor- angular electron-positron measured we Recently, c 0 = 2 =16.70 + 1tastosin transitions M1 , e vdnespotn h xsec ftehptei X17 hypothetic the of existence the supporting evidence New ..Krasznahorkay, A.J. π ± falgtpril ihms of mass with particle light a of 115 ltn h 10 e 0 MeV 21.01 the ulating et 1,020 21) o eciigteaoayobserv whi anomaly particle, the X17 describing same for the (2016)] likely 052501 is 116, it Lett. mass, the to According ASnmes 32.a 32.n 14.70.Pw 23.20.En, 23.20.Ra, numbers: PACS 0 0 eosre lcrnpsto ar rmteelectro-magn the from pairs electron-positron observed We → . ◦ 8 5sa ) 35(stat enucleus. Be → ih7.2 with Z π 8 nttt fNcerRsac,(tmi,PO o 1 H-400 51, Box P.O. (Atomki), Research, Nuclear of Institute Be nttt fNcerRsac Aok) ..Bx5,H-4001 51, Box P.O. (Atomki), Research Nuclear of Institute ′ 1 = e + + ǫ ∗ e u γ + − + → σ yteN4/ xeiet[7], experiment NA48/2 the by etrguebsnX7 which X17, boson gauge vector ǫ ± process. 8 infiac,adcudb ecie yasmn h creatio the assuming by described be could and significance, d e+X7 n hndecays then and X17, + Be Z 0 nvriyo ercn 00Dbee,P o 0,Hungary 105, Box PO Debrecen, 4010 Debrecen, of University ≤ . ′ 8 (y)MV h possible The MeV. 5(sys) ou n onqak to quarks down and up to eadaoaosangu- anomalous and Be 10 ∗ − .Calo,L sg,J uya,M oza .Siam,adJ Tim´ J. and Szihalmi, B. Koszta, Gulyás, M. J. Csige, Csatlós, L. M. eǫ 3 − u 5 ] Subsequently, 6]. [5, tt in state and µ ..Firak, D.S. uze triggered puzzle, m eǫ EN eea wteln and Switzerland Geneva, CERN, X d c 4 2 e ekwsosre ntheir in observed was peak A He. π fu and up of =16.84 1 = .Krasznahorkay A. + .Ng,adNJ Sas N.J. and Nagy, A. ´ ± → 0 . 16( stat ihete nelre ig,o ag etr rboth. the or investigated sector, gauge or Higgs, enlarged frameworks that an BSM either scenarios include with They (BSM) these. Model accommodate Standard can describe and the [12] possible Beyond showed are solutions the also 1 and low They spin-0 of both variety how constraints. a evade experimental to as scale so couplings with axial breaking, symmetry significant U(1)0 the from stemming state, ehv ofimdtesga fteasmdX7particle X17 assumed the of signal the confirmed have We nml a edsrbdwt eylgtZ light very a with described be can anomaly etdt ea rdmnnl lob h msinof emission the by GeV ex- also the are predominantly (ALPs) e at decay particles like problem to axion CP pected or strong axions Such the pos- scale. the solving reopens of This sibility a emission. and axion nuclear via including hadronically, signatures, de-excitations coupled axion scale low-energy new of predicts GeV variety model the This at compat- particular, remains particles constraints. In that existing model with window. axion ible variant mass a MeV) found O(10 QCD they the on and Alves constraints in problem, experimental axions CP revisited strong of [10] the excluded significance Weiner to long the solution Considering axion been the have searches. to elec- laboratory the believed by below are or but at scale axion, constants ef- troweak invisible decay current ultra-light, with major axions an are visible for There searching in problem. forts CP strong the results. to our MeV with 18.15 agreement the nice to in compared is They which transition one, momentum. MeV case orbital 17.6 in ratio 1 the branching of = smaller times L ten with about predicted emitted was it if unu-ubr fthe of quantum-numbers oo ol nedb a be indeed could boson + nprle oteercn hoeia tde,w re- we studies, theoretical recent these to parallel In el oeadc-okr 1]soe htthe that showed [11] co-workers and Rose Delle h C xo soeo h otcmeln solutions compelling most the of one is axion QCD The e ) − ± pairs. 0 . din ed 20( tclyfridnM rniindepop- transition M0 forbidden etically hw eetysgetd[hs Rev. [Phys. suggested recently we ch syst 8 Be. e n Γ and MeV ) 8 ercn Hungary Debrecen, 1 e eaoaywt nipoe setup. improved an with anomaly Be + ercn Hungary Debrecen, e − n usqetdecay subsequent and n nua orltosat correlations angular J 8 X π Be 3 = 0 = ∗ . and 9 − × suoclrparticle, pseudoscalar particle 10 8 esae,teX17 the states, Be − 5 eV. ar 0 bosonic 2

2 and constrained its mass (mXc = 17.01(16) MeV) and for both cases (in-beam and off-beam). In this energy branching ratio compared to the γ-decay (Bx = 6(1) × range, no in-beam counts were expected. The elimina- 10−6) [13, 14]. We also re-investigated the e+e− pair tion of cosmic-ray background was then performed until correlation in the 17.6 MeV transition of 8Be, in which a all events disappeared with this high-energy gate. much smaller deviation was observed [15]. In order to reduce the cosmic-ray background, an ac- In order to confirm the existence of the X17 particle tive shield was also installed above the e+e− spectrom- we have conducted a search for its creation and decay in eter. It consisted of 12 pieces of 1.0 cm thick, 4.5 cm − the 21.01 MeV 0 → 0+ transition of 4He. Emission of a wide and 100 cm long plastic scintillators. Half of the 2 + mXc = 17 MeV vector boson (Jπ=1 ) or pseudoscalar cosmic-ray yield could be suppressed in this way. − particle (Jπ=0 ) is allowed in this transition with orbital In the original total energy spectrum of the e+e− pairs angular momentum 1 or 0, respectively. In this Letter we determined using all combinations of the telescopes we report on anomalous creation and angular correlation of have got a very large background from external pairs electron-positron pairs in this transition, which is in good created by γ-rays coming from the direct proton cap- agreement with the scenario of its decay by the assumed ture, which has a cross section of about 2 µb/sr at Ep= X17 particle. 900 keV [18]. In order to reduce the background from 3 4 The H(p,γ) He reaction at Ep=900 keV bombarding the external pair creation, we have created two total- energy was used to populate the wide (Γ = 0.84 MeV) energy spectra. One was constructed from e+e− pairs, − 4 0 second excited state in He [16], located at Ex= 21.01 which were detected by telescope pairs with relative an- MeV. This bombarding energy is below the threshold of gles of 120◦, while the other from e+e− pairs which were 4 ◦ the (p,n) reaction (Ethr=1.018 MeV) and excites the He detected by telescope pairs with relative angles of 60 . + − nucleus to Ex=20.49 MeV, which is below the centroid of Since the e e pairs from the X17 boson are expected − the wide 0 state. A proton beam with a typical current at around 115◦, the first spectrum is expected to contain 3 3 of 1.0 µA was impinged on a H target. The H was the majority of such events, while the second is expected 2 absorbed in a 3 mg/cm thick Ti layer evaporated onto to be mainly background. To enhance the X17 boson 3 a 0.4 mm thick Mo disc. The density of the H atoms events, we subtracted the second spectrum after appro- 20 2 was 2.66 × 10 atoms/cm . The disk was cooled down priate normalization from the first spectrum. 3 to liquid N2 temperature to prevent H evaporation. Fig. 1 shows the resulting spectrum of the e+e− pairs. The investigated 0− state overlaps with the first ex- 4 π + cited state in He (J =0 ,Ex=20.21 MeV, Γ=0.50 MeV), which was also excited but but give only a man- agable background to the e+e− spectra. The experiment was performed at the 5 MV Van de 160 Graaff accelerator in Debrecen. Compared to our pre- Counts vious experiment [1, 17], we increased the number of 140 telescopes (from 5 to 6) and we replaced the gas-filled MWPC detectors to a double-sided silicon strip detector 120 (DSSD) array. 100 The e+e− pairs were detected by six plastic scintillator + DSSD detector telescopes placed perpendicularly to 80 the beam direction at azimuthal angles of 0◦, 60◦, 120◦, 180◦, 240◦ and 300◦. The sizes of the scintillators are 60 × × 3 82 86 80 mm each. The positions of the hits were 40 registered by the DSSDs having strip widths of 3 mm Signal region and a thickness of 500 µm. The telescope detectors were 20 placed around a vacuum chamber made of a carbon fiber 0 tube with a wall thickness of 1 mm. 18 18.5 19 19.5 20 20.5 21 21.5 22 22.5 γ rays were also detected for monitoring. A ǫrel=50% E(sum) MeV HPGe detector was used at 25 cm from the target. In the investigation of such rare processes, the cosmic FIG. 1. Background subtracted total energy spectrum of the − ray background needs to be taken into account. The e+e pairs. background was measured for two weeks before and after the experiment, and was subtracted out by using the The spectrum (black dots with error bars) is originated same gates and conditions as for the in-beam data. In from the 20.21 MeV E0 transition and from the low- order to determine the normalization factor for the sub- energy part (Ex=20.49 MeV) of the broad 21.01 MeV traction, a gate of 25 MeV ≤E(sum)≤ 50 MeV was used electro-magnetically forbidden 0− → 0+ M0 transition to determine the angular correlations of the cosmic rays in 4He. 3

The efficiency calibration of the telescopes was per- sharp bump at around 115◦. This feature is similar to formed by using the same dataset but with uncorrelated the anomaly observed in 8Be, and seems to be in agree- pairs from consecutive events. Accordingly, an energy- ment with the X17 boson decay scenario. The green full independent efficiency curve could be extracted. The curve shows the simulated angular correlation including energy dependence of the efficiency calibration was sim- the decay of the expected X17 particle, which was fitted ulated by the GEANT3 code (for the same e+e− sum- to the data. In order to derive the exact value for the energy gate as we used in the experimental data reduc- mass of the decaying particle from the present data, we tion) and taken into account as a minor correction on the carried out a fitting procedure for both the mass value experimentally determined efficiency curve. and the height of the observed peak. Fig. 2 shows our experimental results (red asterisks The fit was performed with RooFit [19] by describing with error bars) for the angular correlation of e+e− the e+e− angular correlation distribution with the fol- pairs gated by the total energy of the signal region (19.5 lowing probability density function (PDF): MeV≤ Etot ≤22.0 MeV), using the asymmetry parame- + − ter (−0.5 ≤ y ≤ 0.5) as defined in Ref.[1] and corrected P DF (e e )= NBg ∗ P DF (exp)+ NSig ∗ P DF (sig) , for the relative efficiency of the spectrometer. Black stars (1) with error bars show the angular correlation of e+e− pairs where P DF (exp) was determined experimentally for for the background region (5 MeV≤ Etot ≤19 MeV and the background region, P DF (sig) was simulated by −0.5 ≤ y ≤ 0.5). GEANT4 for the two-body decay of the X particle as a function of its mass, and NBg and NSig are the fitted number of background and signal events, respectively. The signal PDF was constructed as a 2-dimensional model function of the e+e− opening angle and the mass of the simulated particle. To construct the mass depen- − 3H(p,e+e-)4He dence, the PDF linearly interpolates the e+e opening angle distributions simulated for discrete particle masses. E = 900 keV Using the composite PDF described in Equation 1 we

PCC (relative) p -1 first performed a list of fits by fixing the simulated parti- 10 cle mass in the signal PDF to a certain value, and letting RooFit estimate the best values for NSig andNBg. Let- ting the particle mass lose in the fit, the best fitted mass is calculated for the best fit and shown also in Fig. 2. in green. The significance of the peak observed in the e+e− angular correlations was found to be 7.2σ. The mass of 2 the particle derived from the fit is: mXc =16.84±0.16 MeV. -2 − 10 The partial width of the boson-decay ΓX from the 0 - state to the ground state is estimated as follows:

σ(X17) σ(0+) 40 50 60 70 80 90 100 110 120 130 ΓX/ΓE0 = − . (2) Θ (degree) σ(E0) exp. σ(0 )th

The σ(X17)/σ(E0)= 0.20 was obtained from the fit of the + − − simulated e e angular correlations to our data. The FIG. 2. Angular correlations for the e+e pairs measured in 3 4 exciataion energy of the nucleus in the case of Ep=900 the H(p,γ) He reaction at the Ep=900 keV. keV is 20.49 MeV. At that energy the contribution of the According to our simulations such background is orig- − 0+ and 0 resonances are equal, so the ratio of the cross inated from external pair creation on the target backing sections in the resonant proton capture reaction could be and on the surrounding materials (black histogram) and + − π + + calculated as follows: from the IPC e e pairs crated in the J =0 → 0gs. E0 transition (magenta histogram). The sum of that two σ(0+) Γ (0+) = tot =0.59. (3) components fitted to the experimental data is shown as a − − σ(0 ) th Γtot(0 ) blue histogram. The data measured for the background −4 were fitted by a 4-th order exponential polynomial, and Since ΓE0=(3.3 ± 1) × 10 is known [20], than −4 −5 the result is shown in a blue full curve. This blue curve ΓX=0.2 × 0.59 × 3.3 × 10 =3.9 × 10 eV. 8 −6 was rescaled to fit the background of the angular correla- In the case of Be it was ΓX=Γγ ×BX =1.9×6×10 tion shown in red in the range of 40◦ ≤ θ ≤ 90◦. The ob- eV = 1.2×10−5 eV. This value is indeed expected to be tained experimental angular correlation exhibits a rather smaller due to the phase space correction factor. 4

The signiﬁcance of the peak observed in the e+e− in- 100 variant mass distribution was found to be 7.1σ. The mass 2 3 + - 4 of the particle derived from the ﬁt is: mXc =17.00±0.13 90 H(p,e e ) He

Counts MeV. This value agrees within the erroor bar with the one we derived from the fit of the angular correlation. 80 Ep= 900 keV The systematic uncertainties was estimated by taking 70 into account the uncertainty of the target position along the beam line, which was estimated to be ± 2 mm, which 60 2 may cause mXc ± 0.06 MeV uncertainty. The uncer- 50 tainty of the place of the beam spot perpendicular to the ± 40 beam axis was estimated to be in worst case also 2 mm, which may cause a shift in the invariant mass of 2 2 30 mXc ± 0.15 MeV/c . The whole systematic error was 2 conservatively estimated as: mXc ±0.20 MeV. 20 The obtained mass value agrees very well with that 10 of derived for the X17 boson from the 8Be experiments. This is remarkable taking into account the fact that in 0 10 11 12 13 14 15 16 17 18 19 4 2 the present He transition the anomalous bump in the Invariant mass (MeV/c ) angular correlation spectrum appears at a quite different angle than it appears in the 8Be experiments due to the FIG. 3. Invariant mass distribution derived for the 20.49 MeV 4 different energies of the two excited states. The good transition in He. agreement between the two masses leads to the scenario of decaying both studied excited states by the same X17 The invariant mass distribution was also calculated particle. This strengthens the validity of the X17 bo- from the measured energies and angles of the same son hypothesis. It is also worth mentioning that strictly dataset: speaking it cannot be proved that in the 4He case the − + 2 anomalous decay belongs to the 21.01 MeV 0 → 0 2 2 2 (1 + y ) + mXc = 1 − y E sin(θ/2)+2m 1+ cos(θ) , transition. The wide 20.21 MeV 0 first excited state e − 2 − p (1 y ) overlaps with the 21.01 MeV 0 state, and they both were populated in the experiment. However, the anoma- where E = E + + E − and y = (E + − E − )/(E + + + e e e e e lous decay of the 0 state would result a different new E − ). The result is shown in Fig. 3 for the signal − e particle than the decay of the 0 state or the decay of (19.5 MeV≤ E ≤22.0 MeV, in red) and background + 8 tot the 1 state in the Be. Assuming two new particles with (5 MeV≤ E ≤19 MeV, in black) regions. tot the same mass is a less probable scenario than assuming The observed local p0 probability as a function of mX, only one X17 particle, which explains both anomalies. associated to the invariant mass distribution is shown in We are expecting independent (particle physics) ex- Fig. 4. It is the probability that the observed excess perimental results to come in the coming years. In the is due to a statistical oscillation of the background, as following we cite a few of them. defined and used in high energy physics [21]. Recently, the NA64 experiment [22] at CERN presented the first direct search with a 100 GeV/c e− beam 2 for this hypothetical mXc =16.7 MeV boson and ex- 0 1 cluded part of its allowed parameter space, but left the − 10 1 2σ × −4 ≤ ≤ × −3 −2 still unexplored region 4.2 10 ǫe 1.4 10 as Local p 10 σ − 3 10 3 quite an exciting prospect for further research. Experi- 10−4 4σ ment will be continued [23, 24]. − 10 5 − The goal of ForwArd Search ExpeRiment (FASER) 10 6 5σ 10−7 [25] at the LHC is to discover light, weakly interacting − 3 10 8 particles with a small (1 m ) detector placed in the far- − 6σ 10 9 − forward region of ATLAS. In particular, Ariga and his 10 10 −11 coauthors [26–30] considered the discovery prospects for 10 7σ 10−12 ALPs. The project has already been approved, and the − 10 13 experiment will start in 2023. 10 12 14 16 18 2 mXc (MeV) Jiang, Yang and Qiao [31] presented a comprehensive investigation on the possibility of search for the X boson + − FIG. 4. The observed local p0 as a function of the hypothesized directly in e −e collisions, and through the decay of the − X17 boson mass for the X17 → e+e channel. created J/ψ particles at the BESIII experiment for both 5 spin-0 and spin-1 hypotheses. They suggest that Z0-like [4] Xilin Zhang, Gerald A Miller Phys. Lett. B773 159 boson signal might be found or excluded in the present (2017). run of BESIII. The BESIII experiment has accumulated [5] J. Feng et al., Phys. Rev. Lett. 2016 117, 071803 (2016). D 95 the largest J/ψ dataset ( 1010 J/ψ events) worldwide. [6] J. Feng et al. Phys. Rev. , 035017 (2017). [7] J. Batley et al. (NA48/2 Collaboration), Phys. Lett. B They found that this is an ideal channel to test the spin 3 746, 178 (2015). of the particle. They are expecting ≈ 10 scalar/Z0-like [8] Luigi Delle Rose et al., Phys. Rev. D 96, 115024 (2017) X bosons when setting the reduced Yukawa coupling pa- and references therein. − rameters to 10 3, which is within the analysis sensitivity [9] U. Ellwanger and S. Moretti, JHEP 11 39 (2016). of BESIII. [10] D.S.M. Alves, and N. J. Weiner, High Energy Phys. 92, Nardi and coauthors [32] suggested the resonant pro- (2018). [11] Lugi Delle Rose et al., Phys. Rev. D 99 055022 (2019). duction of X17 in positron beam dump experiments. [12] Lugi Delle Rose et al., Frontiers in Physics 7 73 (2019). They explored the foreseeable sensitivity of the Frascati [13] A.J. Krasznahorkay et al., J. Phys.: Conf. Series 1056, PADME experiment in searching with this technique for 012028 (2017). 8 the X17 boson invoked to explain the Be anomaly in nu- [14] A.J. Krasznahorkay et al., Acta Phys. Pol. B 50, 675 clear transitions. PADME already took some test data (2019). and is running until the end of 2019 [33–37]. After that, [15] A.J. Krasznahorkay et al., Proceedings of Science, the experimental setup will be moved to Cornell and/or (Bormio2017) 036 (2017). [16] D.R. Tilley, H.R. Weller, G.M. Hale, Nucl. Phys. A541 JLAB to get higher intensity positron beams. 2 1 (1992). DarkLight will search for 10 - 100 MeV/c dark pho- [17] J. Gulyás et al., Nucl. Instr. and Meth. in Phys. Res. A 8 tons [38]. The sensitivity is projected to reach the Be 808, 21 (2016). anomaly region. The first beam was already used in sum- [18] K.I. Hahn, C.R. Brune, R.W. Kavanagh, Phys. Rev C mer 2016. Currently, they are doing proof-of-principle 51, 1624 (1995). measurements [39]. [19] W. Verkerke and D. P. Kirkby, “The RooFit toolkit for + − data modeling,” eConf C 0303241 (2003) MOLT007 In summary, we have observed e e pairs from an [physics/0306116]. electro-magnetically forbidden M0 transition depopulat- [20] Th. Walcher, Phys. Lett. 31B 442 (1970). − 4 ing the 21.01 MeV 0 state in He. The energy sum of [21] E. Gross, CERN Yellow Reports: School Proceedings, the pairs corresponds to the energy of the transition. The Vol. 4/2017, CERN-2017-008-SP (CERN, Geneva, 2017). measured e+e− angular correlation for the pairs shows [22] D. Banerjee, et al., on behalf of the NA64 collaboration a peak at 115◦, supporting the creation and decay of Phys. Rev. Lett., 120, 231802 (2018). 2 [23] M. Kisarov., EPJ Web of Conferences 212, 06005 (2019). the X17 particle with mass of mXc =16.84±0.16(stat) ± [24] D Banerjee, on behalf of the NA64 collaboration, 0.20(syst) MeV. This mass agrees nicely with the value 2 arXiv:1909.04363, 2019. of mXc =17.01 ±0.16 MeV we previously derived in the [25] J. L. Feng, et al., Phys. Rev. D 98, 055021 (2018). 8 Be experiment [1, 13, 14]. The partial width of the X17 [26] A. Ariga et al., for the PHASER collaboration, arXiv: −5 particle decay is esimated to be: ΓX=3.9 × 10 eV. We 1811.10243 (2018). are expecting more, independent experimental results to [27] A. Ariga et al., for the PHASER collaboration, arXiv: come for the X17 particle in the coming years. 1901.04468 (2019). [28] A. Ariga et al., for the PHASER collaboration, We wish to thank D. Horváth for the critical reading arXiv:1811.12522 (2019). the manuscript and for the many useful discussions. We [29] A. Ariga et al., for the PHASER collaboration, wish to thank also for Z. Pintye for the mechanical de- arXiv:1812.09139 (2019). sign of the experiment. This work has been supported by [30] A. Ariga et al., for the PHASER collaboration, Phys. the Hungarian NKFI Foundation No. K124810, by the Rev. D99, 095011 (2019). GINOP-2.3.3-15-2016-00034 and by the János Bolyai Re- [31] Jun Jiang, Hao Yang, Cong-Feng Qiao, Eur. Phys. J. search Fellowship of the Hungarian Academy of Sciences C79, 404 (2019). [32] E. Nardi etal., Phys. Rev. D97,095004 (2018). (L. Csige). [33] C. Taruggi, for the PADME collaboration, Frascati Physics Series Vol.67, 17, (2018). [34] C.D.R. Carvajal, Frascati Physics Series Vol. 67 (2018). [35] G. Piperno for the PADME collaboration, Journal of Physics: Conf. Series1162 (2019) 012031. [36] V. Kozhuharov, EPJ Web of Conferences 212, 06001 ∗ [email protected] (2019). [1] A.J. Krasznahorkay et al., Phys. Rev. Lett. 116, 042501 [37] S. Spagnolo on behalf of the PADME Collaboration, (2016). Journal of Physics: Conf. Series1137 (2018) 012043. [2] A. Datta, J. Feng, S. Kamali, J. Kumar, [38] R.Corliss, DarkLight Collaboration, Nucl. Inst. Meth. arXiv:1908.08625 865, 125 (2017). [3] http://inspirehep.net/search?ln=en&p=refersto%3- [39] Y. Wang et al., (DarkLight collaboration) Nuclear Inst. Arecid%3A1358248 and Methods in Physics Research, A935, 1 (2019).