Arxiv:2006.02969V3 [Hep-Ph] 17 Jan 2021 H Rnhn Rcin[]Frti Ea [6], Decay This Predicts for Model [5] Standard Fraction the Branching the 10:1

Home , Pionium

arXiv:2006.02969v3 [hep-ph] 17 Jan 2021 h rnhn rcin[]frti ea [6], decay this predicts for Model [5] Standard fraction The branching the 10:1. of ratio background hoeia rnhn rcin[]is [6] fraction branching theoretical eeao eerhCmlx twsdsge oobserve to the designed Ac- decay was at It the Proton conducted Japan Complex. the Research being celerator at is Facility Experimental [7] Hadron experiment KOTO The olc enough collect omdi e.[]sosta ftefu vnsi the to in equal branch- events mechanism a underlying four mean the the would of it if fraction real, that ing were shows region [9] signal KOTO Ref. per- analysis in The formed result. surprising this of interpretations vnswr xetd eyso,svrlppr ap- papers several soon, Very 0.10 mere peared, where expected. situation in were the events Col- events in four KOTO region of the signal presence the of the upper announced behalf September [8] the (on in laboration) Shinohara down However, Satoshi pushing fractions. 2019, pro- branching slowly experiments of expected, decay bounds kaon as rare ceeded both recently, Until tte6(5%cndnelevel. confidence 68(95)% the at hs xeiet,tepoim eepoue ythe by In produced were Synchrotron. pioniums Proton the CERN experiments, (DIRAC) the these Complex at Atomic [13] Energy Relativistic experiment High Dimeson studied of the intensively and Institute in [12] the Russia Serpukhov, at at 1993 Physics in discovered was ]daswt oiieycagdkaons charged NA62 [1– positively The Synchrotron with Proton deals neutrinos. Super two 4] CERN and the constrain pion at to experiment a and rare into aim very Model the decays main measuring Standard kaon precisely The the by theories test physics to running. new is currently both are of flight in decays The w motn xeiet netgtn h aekaon rare the investigating experiments important Two B π ( e.g. + B K B π ( K ( L K − es 91] ie tfidn e physics new finding at aimed [9–11], Refs. , K L tt a osiueantytcniee akrudt h v the to background considered yet not a constitute can state fectdpoim.T lutaeteiprac fti ba this unexpected of the importance assumptions, the some illustrate under To explain, pionium). excited of ato inusmyecp h ea einbfr decaying before region decay the escape may pioniums of part experiment. → + L tm inu uulydntda A as denoted (usually pionium atom, → esgetta h ea oe fkoswt inadapioni a and pion a with kaons of modes decay the that suggest We → → K π nttt fPyisadRsac etefrComputational for Centre Research and Physics of Institute 0 π + ν inu sasuc ffleeet nthe in events false of source a as Pionium π π 0 ν ¯ + 0 ν → ) ν ν ν ¯ KOTO ν ¯ ν ¯ (3 = ) fln-ie eta an.The kaons. neutral long-lived of (8 = ) π + ν 2 = ν ¯ ieinUiest nOaa 4 1Oaa zc Republic Czech Opava, 01 746 Opava, in University Silesian . . 4 4 vnst e inlto signal a get to events . ± 1 ± +2 − 0 1 1 . . . . 6) nttt fEprmna n ple Physics, Applied and Experimental of Institute 0(+4 1( 0) − × × 1 K zc ehia nvriyi Prague, in University Technical Czech . . 1) 7) 10 10 + 2 0Pau,CehRepublic Czech Prague, 00 128 × − − n isto aims and 11 11 10 . . − 9 ± ee Lichard Peter 0.02 2 (1) (2) (3) π ), and a esrdi e.[5 ihteresult the with [15] Ref. in measured was (A so-called [15] observed long-lived Collaboration DIRAC The nium. uhaln ieiei asdb h atta h decay states the that C-parity fact positive the the by caused to is modes lifetime long a Such et oti n,w eda es rd estimate crude a decay least at of experi- need fraction KOTO we branching the end, the in this of it To estimate will ment. we amount, tant o obde n h lw2p slow the and forbidden now sueta h rnhn ratio branching the that assume h akrudi the in background the pionium. any into decay kaon neutral two to decay a state, 1s the reaching After h A82Claoaini hi eia ae 1]de- ratio [14] branching that paper the seminal important termined their two- is in subsequent it Collaboration NA48/2 considerations, their the later and our decays For kaon decay. the in niums A 2 olw:A follows: n h esrdlftm s[13] is lifetime measured the and momenta. ap- relative with small constituents and target, energies same their equal the into proximately In broke-up target. pioniums a the on impinging beam proton studied nthe in K m 2 ′ L nte motn icvr ocrsteeie pio- exited the concerns discovery important Another nodrt hwta h inusmycnrbt to contribute may pioniums the that show to order In h ea ftecagdkosit xie pionium excited into kaons charged the of decay The h ea fpoimt w eta in sdominant is pions neutral two to pionium of decay The h A82Claoaina h ENSS[14] SPS CERN the at Collaboration NA48/2 The R 2 ′ π π π ea vnsta perdi h KOTO the in appeared that events decay + ihqatmnumbers quantum with ) = a o enrpre e.Nihrhsbe the been has Neither yet. reported been not has htcnb nepee stepouto fpio- of production the as interpreted be can that π Γ( 0 K Γ( π krud eso hti a even may it that show we ckground, notwo into K r aedecay rare ery 2 ′ 0 ± π K π + ± nain asdsrbto ntevcnt of vicinity the in distribution mass invariant hsc n aaProcessing, Data and Physics → ± → π → − K → τ τ A 2p R 1s π m( um ˜ π tm,wihaete2 tmcstates atomic 2p the are which atoms, 2 ± π π ± π → = 3 = 0 = 0 ± π γ π o to (or s 0 A π + Γ( → π Γ( + . . πν 2 π 15 K 45 0 π K π K − π K ) − − +0 easadfuda anomaly an found and decays − +1 L → ) 0 ν ¯ L 0 → 0 π tm ntefinal the in atom) π (1 = → . . . . 28 26 08 30 0 → πν 0 decays γ πν π → π × . × J 0 ν π ¯ 0 γ . ν ¯ stasto dominates. transition 1s C P 61 10 0 π 10 nfc,a fact, In . xeiet yimpor- by experiments ntecase the in A + − − ± 2 π K 1 = π 15 11 − 0 ) L ) . π s s 66) . −− . 0 → π 0 t lifetime Its . × π π and 10 0 0 A quickly s − 2 5 π γγ . We . are (5) (6) (4) π 0 2 has the same value as that for charged kaons (5). Then If we denote the length of the signal region as L, the we can write the branching-fraction estimate mean survival probability of pionium at the point where 0 0 + − it leaves the signal region is B(KL → π A2π)= R ×B(KL → π π π ) ≈ 2 × 10−6, (7) 1 L l S¯ = S(L − z)dz = [1 − exp{−L/l}] . (11) where we have also consulted Ref. [16]. L Z0 L We will also consider possibility that the pionium that Numerical values for two kinds of pionium are [18] appears in the kaon decays is not the ground state (1s) pionium, but its excited (2p) partner and that the corre- ¯ +0.16 −6 S1s =1.75−0.15 × 10 , sponding branching fraction are the same, ¯ +6.0 −3 S2p =2.5−1.7 × 10 . 0 ′ −6 B(KL → π A ) ≈ 2 × 10 . (8) 2π Multiplying these numbers by the assumed branching In what follows, we will pursue those two alternatives in fractions (7) and (8), we obtain the branching fractions 0 0 ′ parallel. of KL → π A2π and KL → π A2π events that look like 0 To simplify the reasoning, we will ignore the momen- the KL → π νν¯ events because pioniums left the KOTO tum spread of the KL beam in the KOTO experiment decay volume undecayed, and will use its peak value P = 1.4 GeV/c [7]. The +0.32 −12 length of the KOTO signal region is L=1.7 m. Another B1s =3.50−0.30 × 10 , (12) +12.0 −9 quantity that enters the game is the pionium mass. It B2p =5.0−3.4 × 10 . (13) is given by ma = 2mπ+ − b, where b is the binding energy. We will take its coulombic value, which can be cal- Branching fraction (12) is by 3 orders of magnitude culated from the hydrogen-atom-like formula. One ob- smaller than branching fraction (3), which characterizes tains b = 1.86 keV for the pionium ground state and the presence of unexpected events in the KOTO signal re- b =0.464 keV for the excited 2p state. gion. Therefore, we will not follow the 1s pionium option The laboratory energy of pionium with mass ma is any longer. uniformly distributed in an interval, the bounds of which A comparison of (13) with (3) suggests that the anoma- are given by the formula lous events in the KOTO experiment could be explained as undecayed 2p pioniums. But to make a realistic com- 1 ∗ ∗ parison, we must take into account that the experimental Ea± = (EEa ± Ppa) , 0 0 ′ M efficiencies for the KL → π νν¯ and KL → π A2π are dif- where M, E, and P are the mass, energy, and momentum ferent. The right quantity that should be compared with ∗ 2 2 ∗ (3) is of the KL, respectively, Ea = (M + ma)/(2M), pa = 2 2 (M − ma)/(2M). Numerically, Ea− = 0.497 GeV and ′ B2p = R ×B2p, (14) Ea+ =1.456 GeV. To simplify the consideration further [17], we will as- where R is the ratio of efficiencies, sume that all pioniums have the same laboratory momentum, given as the mean value, ǫπ0A′ R = 2π . (15) Ea+ ǫπ0νν¯ 1 2 2 1 pa = E − m dE = E − E a 2(E − E ) a+ a− ZEa− p a+ a− Recently, the authors of Ref. [9] considered a similar ratio, 2 Ea+ + pa+ × Ea+pa+ − Ea−pa− − ma log , (9) Ea− + pa− ǫπ0X0 R(mX0 )= , (16) ǫπ0νν 2 2 ¯ where pa± = Ea± − ma. Numerically, pa = q 0.880 GeV/c. where X0 is an invisible boson. They performed a Monte The probability that pionium travels the path s with- Carlo simulation taking into account the experimental out decaying is given by cuts, acceptances, and resolutions. The dependence of ratio (16) on the boson mass mX0 can be obtained from S(s) = exp{−s/l}, the curve in the left pane of Fig. 2 [9]. In principle, it should be possible to use their results and get ratio (15) where l is the mean decay length of pionium, 0 ′ by setting mX to mA2π . Unfortunately, the value for ′ paτ mA2π , which is close to 280 MeV, cannot be read from l = , (10) the curve. We will therefore make our own crude estimate ma of the efficiency ratio (15). τ being the pionium mean lifetime. For the two types The efficiencies are most influenced by the pT cut im- +0.27 of pioniums, we get l1s = 2.98−0.25 µm and l2p = posed in the KOTO experiment. To suppress the events +10.2 + − 0 0 4.2−2.8 mm. from the KL → π π π , a pT momentum of π greater 3 than 130 MeV/c is required over the majority of the sig- that its mass must be greater than 2mπ0 . The binding nal region; at the downstream edge, the cut is even higher energy b should thus be smaller than 2(mπ+ − mπ0 ) ≈ (up to 150 MeV/c). This cut also suppresses the events 9.19 MeV. It is possible [20] that the binding energy of 0 ′ with two-body final states containing, besides the π , a A2π may be around 9 MeV [21]. In that case, we have 0 massive particle X . It may be an invisible light boson, ma = 270.14 MeV, pTmax,A′ = 139.0 MeV/c, ǫπ0A′ = as in Ref. [9], or a pionium. In these events, the maximum (139 − 130)/139= 6.5%, and transverse momentum of π0 is equal to the momentum of the outgoing decay products in the K rest frame, L R(b =9 MeV)=0.29.

1 2 2 2 pT,max = λ(mK ,mπ,mX0 ) , (17) 2mK q The eﬃciency corrected branching fraction now comes where we use the usual notation out as

2 2 2 λ(x,y,z)= x + y + z − 2xy − 2xz − 2yz. ′ +3.5 −9 B2p(b = 9 MeV )=1.5−1.0 × 10 , 0 The maximum transverse momentum of π in the KL → π0νν¯ events is higher because the lowest invariant mass which is, within the errors, compatible with the branch- of the two-neutrino system is zero (neutrinos’ masses ne- ing fraction (3) characterizing the unexpected KOTO glected), namely, events. 2 2 The KOTO Collaboration continues in analyzing its mK − mπ pTmax,νν = = 230.5 MeV/c. older data and taking new ones [22]. One of the four 2m K events reported in [8] has been shown to be due to a mis- 0 Given a pT cut, we can restrict the masses of X s that take, so only three good candidate events remain. The can be detected as new background estimate is 1.05±0.28 events. The num- ber of expected Standard Model events is 0.04. 2 2 2 2 It is interesting that the KOTO experiment sees more mX0 < mK + mπ − 2mK mπ + pT,max . r q K± than expected. It may be worth remaining that the ± ∓ KL → K e ν¯(ν) decays [23] have not been observed yet. For p = 130 MeV/c, it implies m 0 < 281 MeV. It T,max X Nevertheless, the KOTO Collaboration considered them is obvious that ǫ 0 0 should go rapidly to zero with m 0 π X X [24] as a source of background and varied the branching approaching 281 MeV. Figure 2 in Ref. [9] confirms that. fraction in reasonable bounds. A negligible contribution In order to get estimates of the efficiencies, we first was found. randomly generate the particle momenta within the π0νν¯ system in the KL rest frame using the GENBOD program To conclude, we have suggested that events in which from the old CERN Program Library [19]. Counting the the kaon decays into pionium, which then leaves the de- events with the π0 transverse momentum greater than cay region without decaying, may contribute to the back- the KOTO cut of 130 MeV/c, we get an efficiency ǫπ0νν¯ ground in the K → πνν¯ experiments. In our opinion, of 22%. this source of background should be taken into account Next, we consider the 2p pionium with the Coulombic when analyzing the existing experiments or planning new binding energy, the mass of which is ma = 279.14 MeV. ones [25]. However, to get reliable estimates of this back- Using the formula (17) we get pTmax,π0 =132.0 MeV/c. ground, the discovery and measurement of the following 0 Assuming isotropic decay in the KL rest frame, we can decays will be very valuable: KL → π A2π (a neutral estimate the portion of events with pT,π0 greater than kaon partner of the charged kaon decay already observed 0 ′ ± ± ′ 130 MeV/c as ǫπ0A′ = (132 − 130)/132= 1.5%. It means [14]), KL → π A2π, and K → π A2π. an efficiency ratio (15) of As an illustration of the idea, we have shown that, under some assumptions, the production of 2p pioniums R(Coul.) = 0.069 ′ A2π in the KL decays and their subsequent escaping from and an efficiency corrected branching ratio the signal region may explain the production of unexpected events reported by the KOTO Collaboration [8]. ′ +0.83 −9 Our key assumptions, which may or may not be cor- B2p(Coul.) = 0.34−0.23 × 10 . rect, have been as follows: (i) main contribution comes A comparison with (3) shows that if the experimental from pioniums that are in the 2p state, which lives much cut on transverse momentum of π0 is taken into account, longer than the ground state; (ii) the branching ratio of 0 ′ 0 + − the 2p pioniums with the Coulombic binding energy can- KL → π A2π to KL → π π π is the same as that of + + + + + − not be a source of the unexpected events in the KOTO K → π A2π to K → π π π , which was determined ′ experiment. by the NA48/2 Collaboration [14]; (iii) the A2π binding Yet, no direct measurement of the 2p pionium mass energy is larger than the Coulombic one, which helps pi- or binding energy exists. From its decay chain, we know oniums to overcome the experimental pT cut. 4

ACKNOWLEDGMENTS This work was partly supported by the Ministry of Education, Youth and Sports of the Czech Re- I thank Teppei Kitahara for pointing me toward public Inter-Excellence Project No. LTI17018 and the importance of eﬃciencies and Augusto Ceccucci, European Regional Development Fund Project No. Matthew Moulson, and Taku Yamanaka for useful cor- CZ.02.1.01/0.0/0.0/16 019/0000766. respondence.

[1] E. Cortina Gil et al., J. Instrum. 12, P05025 (2017); S. Zhang, Phys. Rev. D 102, 055005 (2020). Martellotti (on behalf of the NA62 Collaboration), KnE [11] B. Dutta, S. Ghosh, and T. Li, Phys. Rev. D 102, 055017 Energ. Phys. 3, 372 (2018). (2020); X. Liu, Y. Li, T. Li, and B. Zhu, J. High Energy [2] E. Cortina Gil et al. (NA62 Collaboration), Phys. Lett. B Phys. 10 (2020) 197. 791, 156 (2019). [12] L.G. Afanasyev et al., Phys. Lett. B 308, 200 (1993). [3] E. Cortina Gil et al. (NA62 Collaboration), J. High En- [13] B. Adeva et al., Phys. Lett. B 619, 50 (2005); 704, 24 ergy Phys. 11 (2020) 042. (2011). [4] R. Marchevski (on behalf of the NA62 Collabora- [14] J.R. Batley et al. (NA48/2 Collaboration), Phys. Lett. B tion), Evidence for the decay K+ → π+νν¯ from the 633, 173 (2006). NA62 experiment at CERN, Proceedings of the ICHEP [15] B. Adeva et al., Phys. Lett. B 751, 12 (2015); 2020, Prague (Virtual Conference), July 28, 2020, Phys. Rev. Lett. 122, 082003 (2019). https://indico.cern.ch/event/868940/contributions/3815641/.[16] P.A. Zyla et al. (Particle Data Group), Prog. Theor. Exp. [5] We use the terminology of [16]: Branching fraction B is Phys. 2020, 083C01 (2020). the ratio of a partial decay rate to the total decay rate; [17] Alternatively, we also performed the evaluation by gener- branching ratio is the ratio of two partial decay rates or, ating the pionium energy randomly within given bounds. equivalently, of two branching fractions. The ﬁnal results diﬀered from those with mean pa by less [6] J. Brod, M. Gorbahn, and E. Stamou, Phys. Rev. D 83, than 10%. 034030 (2011); A.J. Buras, D. Buttazzo, J. Girrbach-Noe, [18] We did not perform a serious study of error propaga- and R. Knegjens, J. High Energy Phys 11 (2015) 33. tion in (11), but simply calculated the values of S¯ in the [7] E. Iwai et al. (KOTO Collaboration), Nucl. Phys. B, boundary points of error intervals of decay lengths. Proc. Suppl. 233, 279 (2012); J.K. Ahn et al. (J-PARC [19] F. James, Report No. CERN 68-15, 1968. KOTO Collaboration), Prog. Theor. Exp. Phys. 2017, [20] P. Lichard, Phys. Rev. D 102, 073004 (2020). 021C01 (2017); J.K. Ahn et al. (KOTO Collaboration), [21] J. Uretsky and J. Palfrey [Phys. Rev. 121, 1798 (1961)] Phys. Rev. Lett. 122, 021802 (2019). considered the pionium binding energies even higher than [8] S. Shinohara (on behalf of the KOTO Collaboration), J. 10 MeV. Phys. Conf. Ser. 1526, 012002 (2020). [22] N. Shimizu: Search for new physics via the 0 [9] T. Kitahara, T. Okui, G. Perez, Y. Soreq, and K. To- KL → π νν¯ decay at the J-PARC KOTO bioka, Phys. Rev. Lett. 124, 071801 (2020). experiment, Proceedings of the ICHEP 2020, [10] P.S.Bhupal Dev, R.N. Mohapatra, and Y. Zhang, Prague (Virtual Conference), July 28, 2020, Phys. Rev. D 101, 075014 (2020); D. Egana-Ugrinovic, https://indico.cern.ch/event/868940/contributions/3815582/. S. Homiller, and P. Meade, Phys. Rev. Lett. 124, 191801 [23] P. Lichard, Phys. Rev. D 55, 5385 (1997); P. Lichard and (2020); M. Fabbrichesi and E. Gabrielli, Eur. Phys. J. C J. Thompson, arXiv:hep-ph/0004224. 80, 532 (2020); Y. Liao, H.-L. Wang, C.-Y. Yao, and J. [24] T. Yamanaka (private communication). [25] M. Moulson, J. Phys. Conf. Ser. 1526, 012028 (2020).