Volume 236, number 1
PHYSICS LETTERS B
8 February 1990
SEARCH FOR NEUTRAL BOSON IN ORTHOPOSITRONIUM DECAY M. T S U C H I A K I , S. ORITO, T. Y O S H I D A and M. M I N O W A Department of Physics, Faculty of Science, Universityof Tokyo, Tokyo 113, Japan Received 20 October 1989
We have searched for the exotic decay of orthopositronium into a photon and a neutral boson qb, without any assumptions on its lifetime. Our negative result rules out unambiguously the possibility of this exotic decay being the origin of the reported discrepancy of the orthopositronium decay rate.
Recently a precise measurement of the decay rate of orthopositronium (o-Ps) was reported. The measured rate is larger than the theoretical prediction by ten standard deviations [ 1 ]. As a possibility to explain this discrepancy one can speculate on an exotic decay of o-Ps into a photon and a hypothetical neutral boson do. Samuel [2] claimed that the mass of such a particle must be below 5.7 keV for the case of a pseudoscalar in order to be compatible with the electron g - 2 constraint and with negative results of beam d u m p experiments [ 3 ]. However, the g - 2 constraint is not valid in the case cancellations exist among new bosons of different spins and couplings [4]. Samuel's conclusion is also based on the assumption that do decays into two photons only via the electron triangle diagram. Since such bosons would be extraordinary objects not predicted by the standard theory, we should not a priori exclude dynamics which might enhance the two photon decay rate of do relative to the electron triangle diagram contribution. In searches of such an object we should rather regard the coupling to two photons as a free phenomenological parameter. It should also be noted that beam d u m p experiments are only valid for particles which pass through the shield and decay within the decay region. Since the production mechanism and interactions of qb in the shield are very complicated, assumptions on the interrelations between dOcouplings to the photon, the electron and the nucleons are needed to derive the limits [ 3 ]. Therefore Samuel's argument does not completely rule out, in a model independent way, the
possible existence of do with a mass above 5.7 keV. Carlson et al. obtained astrophysical constraints from the evolution of stars [5]. However, these constraints apply to particles with mass below 300 keV. Other experiments relevant to do are searches for the decay of o-Ps into a photon plus an invisible particle (single photon searches) [ 6,7 ]; especially ref. [ 7 ] rules out Samuel's hypothesis in the low mass region. However, if the particle has a strong coupling to photons, therefore a short life, these searches are not reliable especially in the high mass region. In this situation, a search which is independent of the do lifetime is very desirable, We report in this paper on a new search for the 7 + do decay of o-Ps in the high do mass region above 300 keV. As we searched for peaks in the inclusive 7 spectrum, we are sensitive to dOeven with a very short lifetime and free from any assumptions on the decay of do, i.e. the coupling to photons, etc. The only relevant free parameter is the do coupling to the electron. The set-up of the experiment is schematically shown in fig. la. The positron source was 22Na of 5 gCurie, sandwiched between two plastic scintillators ( N E 1 0 2 A ) of 100 Jam thickness each. These scintillators were further sandwiched between two blocks of silica aerogel ( p = 0 . 2 3 g / c m 3) of I cm thickness each. About 60% of the emitted positrons pass through the scintillators and stop in the silica aerogel, forming positronia. The scintillation light signalling the emission of a positron was detected by two photomultipliers ( H a m a m a t s u R647). To avoid the quenching annihilation induced by oxygen in the air
0370-2693/90/$ 03.50 © Elsevier Science Publishers B.V. ( North-Holland )
81
Volume 236, n u m b e r 1
P H Y S I C S LETTERS B eeNo +scintiRotors
N2,
5
I r°geLL (
5cm
>
a
b >0 Z Ill
LL LU
10
TM
w~ I LU rd
[ I00
I
i i [ 11 i l l 500
1000
ENERGY(KEV) Fig. 1. (a) Experimental set-up. (b) Relative photopeak efficiency as a function of energy. The points represent the measured values and the line was obtained by fitting these values.
[ 8 ], the source assembly was placed in nitrogen gas of 1 atm in an aluminum can of 0.5 mm thickness. The photons from the positronium decays were detected by a pure coaxial germanium detector of dimensions 44 mm in diameter and 42 mm in length. The energy calibration, the measurements of the energy resolution as well as the relative photopeak efficiency, ~(k), of the germanium detector were performed by using line 7 peaks from a weak '52Eu source placed near to the target. The energy resolution obtained was 1.1 keV and 1.3 keV FWHM at 100 keV and 500 keV respectively. The relative photopeak efficiency obtained, E(k) as a function of the energy k, is shown in fig. lb. The points represent the measured efficiencies and the efficiency curve was determined by fitting these values to the spline function of degree 2 in the ln(E)-ln(k) plot. One output from the preamplifier of the germa82
8 February 1990
nium detector was fed into a shaping amplifier (Ortec-673 ). The gated integrator output of the shaping amplifier was sent to a PHA for the energy measurement. The data stored in the PHA were read out by a computer for further analysis. The second output from the preamplifier was fed into a timing filter amplifier (Ortec-579) and was then discriminated by a constant fraction discriminator (Ortec-583) for the timing measurement. The decay time distribution of the positronium in the silica aerogel was measured by using the scintillator signal as the start and the germanium signal as the stop of a TAC. Fig. 2a shows the TAC spectrum. A sharp peak of the prompt annihilation is followed by the exponential decay curve of the o-Ps and subsequently by the constant accidental background. The lifetime of o-Ps in the silica aerogel was measured to be 127 + 2 ns by fitting the spectrum between 100 ns and 700 ns after the prompt peak with an exponential decay curve plus a constant. In order to obtain a pure sample of o-Ps decay, a delayed coincidence was required which was issued if the germanium signal arrived between 50 and 230 ns after the prompt annihilation. 3.3 × 108 events were collected in four runs for a total of 1.03× 106 s. The stability check and energy calibrations were performed by utilizing the 511 and 1274 keV peaks during data taking. The resulting energy spectrum is shown in fig. 2b. Narrow peak structures were searched for in this raw data spectrum by checking if the data were consistent with an overall smooth curve plus a sharp peak. The shape of the resonance peak was assumed to be gaussian with an FWHM of Ak. The expected energy spread Ak of the peak was mainly determined by two factors, i.e. the thermal motion of positronium in silica aerogel and the resolution of the germanium detector. The thermal kinetic energies of o-Ps in a silica aerogel are known to decrease rapidly to 0.03 eV within 60 ns after the formation [9]. We assume in this letter a value of 1 eV as a conservative averaged kinetic energy [8]. The energy spread due to this thermal motion has been further folded by the energy resolution of the germanium detector, resulting in an energy spread Ak of 1.2 keV and 1.6 keV FWHM at 100 keV and 500 keV, respectively. The smooth curve was obtained by fitting the data while excluding the region within twice the FWHM of the assumed peak position k o. Normally, the fit-
Volume 236, number 1
PHYSICS LETTERS B
8 February 1990
i0 ~
a
C 2.0
iO4 t6
co t-Z
io 3
o (_)
r%,~
0.8
10Z
...... ,'~,~,~,,i,17i~,~it~' 0.4
0 cn I
~0
100
I
200
I
300 400
I
500
I
u
600
700
I
O.O 0
100
I
I
I
200
300
400
500
ENERGY(KEV)
TIME(ns) 10 ="
I
511Kev
b
10 6
1274 backsc~tler ~'
(,f)
key
Complon edge
I.ZlO ~ 0 C~
0.0
0.2
0.4
0.6
0.8
1.0
1.2
ENERGY(MEV)
ting was done by a polynomial curve of order up to 6. The fitting was done by a spline of degree 3 when the assumed peak position kp is in the energy region between 150 keV and 200 keV and between 320 keV and 370 keV corresponding to the backscatter peak and Compton edge of 511 keV, respectively. Then a resonance at the assumed position kp was superimposed on this smooth curve. The best fit to the data was searched for by varying the strength of the resonance peak.
Fig. 2. (a) A time difference spectrum between the scintillator and germanium pulses. The timing of the delayed coincidence are indicated. (b) The 7 ray energy spectrum taken with delayed coincidence. The vertical scale is counts per 0.336 KeV. The expected 37 spectrum of o-Ps decays (photopeak counts) is also indicated. (c) The original energy spectrum ofo-Ps deduced from the raw data in the valley region (points). The theoretical spectrum is also indicated (line).
The energy region between 100 keV and 500 keV was scanned by moving the resonance position kp by steps of 1 keV. No statistically significant peak was found in this energy region except for around 170 keV, 340 keV and 252 keV, corresponding to the backscattering, the Compton edge of 511 keV 7 and the double escape peak of 1274 keV Y- The upper limits on the number of peak events np (kp) thus obtained from the fitting procedure were converted into the limits on the branching ratio BR (kp) using the formula 83
Volume 236, number 1
np(kp) = BR (kp) X E(kp) X N o . ,
PHYSICS LETTERS B
(1)
where No,, is the total number of o-Ps decays occurring within the timing defined by delayed coincidence. No, can be determined from the 3y spectrum ofo-Ps decays which is deduced from the raw data of fig. 2b in the energy region between 400 keV and 500 keV (valley region). In the valley region the Compton contribution of 511 keV is mainly due to the multiple scattering and rather small. For this purpose the random accidental contribution was first subtracted utilizing a spectrum taken with an accidental coincidence which is defined as 1000 to 1180 ns after the prompt peak. The Compton contribution in the subtracted spectrum was further removed in the valley region by the unfolding process. The response functions of the germanium detector to the various y rays in the valley region were deduced from the interpolation of the measured spectra of monochromatic y ray ofSSSr(514 keV) and ~3Sn(392 keV). After this unfolding process the 3y spectrum n(k) (photopeak counts) of o-Ps decay was reproduced in this energy region. The original energy spectrum could be obtained by dividing n (k) by the energy dependent efficiency, e(k). Fig. 2c shows the original energy spectrum, n(k)/e(k), together with the theoretical spectrum [ 10]. The agreement is excellent. No~, can then be obtained from the normalization factor. It should be noted that the final BR(kp) thus determined depends only on the ratio of the efficiencies, not on their absolute values. Fig. 3a shows the obtained upper limits (95% CL) on the branching ratio as a function of m . [6 ]. The top line indicates the branching ratio which corresponds to the reported discrepancy of o-Ps lifetime [ 1 ]. Our negative result rules out unambiguously the possibility of this exotic decay being the origin of the reported discrepancy of o-Ps decay rate, in the mass region 300 keV~< m¢~< 900 keV. It should be noted that the region m . ~<300 keV is disfavored by astrophysical constraints [ 5 ]. When we consider the spin 0 boson (scalar or pseudo-scalar), the branching ratio can be translated into the coupling constant 6% =g2e/47C [ 11 ]~1. The resulting upper bounds on &~e are of order 10 - 9 in this mass region. ~ The formula relating the branching ratio to the coupling constant can be applied to the case of a scalar, although the formula was derived for the case o f a pseudoscalar.
84
8 February 1990
o-Ps d e c a ~ _ s . c r e p a n c y
10-3 Complon edge
~ °o~ -I--
doubte escape peak
///i Ii/i/~/,,/~'I I////I
L.
upper tlmlfs//
E~
~"
"~
.~ lO-°
.~~7
~
/1~//~.~
b a c k s c o t t e r peak
r~J L. ..Q
10-~
400
0
~bo
~do
~o
rn~(KEV) X 10-9 10
bS ,
(m~:400
kev)
6
4
///
;/////// ////)~/.j /, O,e!br,O,ck///////// ,// newly excl.uded region
2
0
........
0
~
~.............................s.iff..!L P~°~°~ .,............ 2
4
6
'~r'ee
8
10 -9 X 10
Fig. 3. (a) The resulting upper limits at 95% confidence level on the branching ratio of 7 + ~ decay. The top line indicates the branching ratio which corresponds to the reported discrepancy of the o-Ps lifetime. (b) Our limits in the &~e-&r~ plane for m . = 4 0 0 keV. The region excluded by single photon experiments is below the bottom line for the pseudoscalar particle. The upper bounds for &rt derived from the Delbdick scattering experiment and that for &ee(g-- 2 ( 1 ) ) and &c~&ry(g-- 2 (2) ) from g - - 2 constraints are also indicated for the case of the pseudoscalar particle.
In order to clarify the generality of the experiment, we have expressed in fig. 3b our result in the &ee--&~ plane for m ¢ = 4 0 0 keV as an example. 6~w-f2/47r is the coupling constant of q) to photons defined by the interaction lagrangian of•yy= (g4mo)F.~,F"."~ for
Volume 236, number 1
PHYSICS LETTERS B
the pseudoscalar and ~ = (f/2m.)F, uF"J'Ofor the scalar. The region excluded by single photon searches is below the bottom line in the case of the pseudoscalar particle [6,12]. Also indicated in this figure are the regions excluded by g - 2 constraints and the DelbriJck scattering experiment in the case of the pseudoscalar particle as well as the region newly excluded by this experiment [ 2,4,12 ]. In conclusion, our result together with the existing constraints completely rule out the possibility of the exotic decay 7 + 0 being the origin of the observed discrepancy of the o-Ps decay rate in the mass region below 900 keV [5,7]. We are deeply indebted to Professor T. Hyodo for many helpful suggestions. Sincere thanks are due to the staffs of the Radioisotope Center of the University of Tokyo for providing space in the early stage of this experiment.
8 February 1990
References [ 1] c.I. Westbrook, D.W. Gidley, R.S. Conti and A. Rich, Phys. Rev. Len. 58 (1987) 1328. [2] M.A. Samuel, Mod. Phys. Lett. 3 (1988) 1117. [3] J.D. Bjorken et al., Phys. Rev. D 38 (1988) 3375; D.J. Bechis et al., Phys. Rev. Lett. 42 ( 1979 ) 1511. [4] A Sch/ifer, J. Reinhardt, W. Greiner and B. Mtiller, Mod. Phys. Lett. 1 (1986) 1. [ 5 ] E.D. Carlsonand P. Salati, preprint UCB PTII-88/28, LBL26215. [6] G. Carboni and W. Dahme, Phys. Lett. B 123 (1983) 349; U. Amaldi, G. Carboni, B. Jonson and J. Thun, Phys. Len B 153 (1985) 444; V. Metag et al., Nucl. Phys. A409 (1983) 331c. [7] S. Orito, K. Yoshimura, T. Haga, M. Minowa and M. Tsuchiaki, Phys. Rev. Lett. 63 (1989) 597. [8 ] M. Kakimoto, DoctoralThesis, Universityof Tokyo ( 1987). [9IT. Chang, M. Xu and Z. Zeng, Phys. Lett. A 126 (1987) 189. [ 10] G.S. Adkins, Ann. Phys. 146 (1983) 78. [ 11 ] J. Cleymansand P.S. Ray, Lett. Nuovo Cimento 37 (1983) 569. [12] M. Suzuki, Phys. Lett. B 175 (1986) 364.
85