A feasibility study of ortho-positronium decays measurement with the J-PET scanner based on plastic scintillators

Abstract

(ProQuest: ... denotes formulae and/or non-USASCII text omitted; see image)

We present a study of the application of the Jagiellonian positron emission tomograph (J-PET) for the registration of gamma quanta from decays of ortho-positronium (o-Ps). The J-PET is the first positron emission tomography scanner based on organic scintillators in contrast to all current PET scanners based on inorganic crystals. Monte Carlo simulations show that the J-PET as an axially symmetric and high acceptance scanner can be used as a multi-purpose detector well suited to pursue research including e.g. tests of discrete symmetries in decays of ortho-positronium in addition to the medical imaging. The gamma quanta originating from o-Ps decay interact in the plastic scintillators predominantly via the Compton effect, making the direct measurement of their energy impossible. Nevertheless, it is shown in this paper that the J-PET scanner will enable studies of the ...... decays with angular and energy resolution equal to ...... and ......, respectively. An order of magnitude shorter decay time of signals from plastic scintillators with respect to the inorganic crystals results not only in better timing properties crucial for the reduction of physical and instrumental background, but also suppresses significantly the pile-ups, thus enabling compensation of the lower efficiency of the plastic scintillators by performing measurements with higher positron source activities.

Eur. Phys. J. C (2016) 76:445DOI 10.1140/epjc/s10052-016-4294-3

Regular Article - Experimental Physics

http://crossmark.crossref.org/dialog/?doi=10.1140/epjc/s10052-016-4294-3&domain=pdf

Web End = http://crossmark.crossref.org/dialog/?doi=10.1140/epjc/s10052-016-4294-3&domain=pdf

Web End = A feasibility study of ortho-positronium decays measurement with the J-PET scanner based on plastic scintillators

D. Kamiska1, A. Gajos1, E. Czerwiski1,a, D. Alfs1, T. Bednarski1, P. Biaas1, C. Curceanu2, K. Dulski1, B. Gowacz1,N. Gupta-Sharma1, M. Gorgol4, B. C. Hiesmayr3, B. Jasiska4, G. Korcyl1, P. Kowalski5, W. Krzemie6,N. Krawczyk1, E. Kubicz1, M. Mohammed1, Sz. Niedwiecki1, M. Pawlik-Niedwiecka1, L. Raczyski5,Z. Rudy1, M. Silarski2, A. Wieczorek1, W. Wilicki5, B. Zgardziska4, M. Zieliski1 , P. Moskal1

1 Faculty of Physics, Astronomy and Applied Computer Science, Jagiellonian University, S. ojasiewicza 11, 30-348 Krakw, Poland

2 INFN, Laboratori Nazionali di Frascati, CP 13, Via E. Fermi 40, 00044 Frascati, Italy

3 Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090 Vienna, Austria

4 Department of Nuclear Methods, Institute of Physics, Maria Curie-Sklodowska University, Pl. M. Curie-Sklodowskiej 1, 20-031 Lublin, Poland

5wierk Computing Centre, National Centre for Nuclear Research, 05-400 Otwock-wierk, Poland

6 High Energy Department, National Centre for Nuclear Research, 05-400 Otwock-wierk, Poland

Received: 4 June 2016 / Accepted: 1 August 2016 The Author(s) 2016. This article is published with open access at Springerlink.com

Abstract We present a study of the application of the Jagiellonian positron emission tomograph (J-PET) for the registration of gamma quanta from decays of orthopositronium (o-Ps). The J-PET is the rst positron emission tomography scanner based on organic scintillators in contrast to all current PET scanners based on inorganic crystals. Monte Carlo simulations show that the J-PET as an axially symmetric and high acceptance scanner can be used as a multi-purpose detector well suited to pursue research including e.g. tests of discrete symmetries in decays of orthopositronium in addition to the medical imaging. The gamma quanta originating from o-Ps decay interact in the plastic scintillators predominantly via the Compton effect, making the direct measurement of their energy impossible. Nevertheless, it is shown in this paper that the J-PET scanner will enable studies of the o-Ps 3 decays with angular and

energy resolution equal to () 0.4 and (E) 4.1 keV,

respectively. An order of magnitude shorter decay time of signals from plastic scintillators with respect to the inorganic crystals results not only in better timing properties crucial for the reduction of physical and instrumental background, but also suppresses signicantly the pile-ups, thus enabling compensation of the lower efciency of the plastic scintillators by performing measurements with higher positron source activities.

a e-mail: mailto:[email protected]

Web End [email protected]

1 Introduction

The positron emission tomography (PET) is based on registration of two gamma quanta originating from a positron annihilation in matter. However, the e+e 2 process

is not the only possible route of positron annihilation. Electron and positron may annihilate also to a larger number of gamma quanta with lower probability, or form a bound state called positronium. In the ground state with angular momentum equal to zero positronium may be formed in the triplet state (with spin S = 1) referred to as ortho-positronium (o-Ps), or singlet state (S = 0) referred to as para-positronium (p-Ps). Positronium, being a bound-state built from electron and anti-electron bound by the central potential, is an eigenstate of both charge (C) and spatial parity (P) operators, as well as of their combination (CP). Therefore, it is well suited for the studies of these discrete symmetries in the leptonic sector. These symmetries may be studied by the measurement of the expectation values of various operators (odd with respect to the studied symmetry) constructed from the momenta of photons and the spin of the ortho-positronium [1]. Such studies are limited by the photonphoton interaction, however it was estimated that the vacuum polarisation effects may mimic the CP and CPT symmetries violation only at the level of 109 [2], which is still by six orders of magnitude less than the presently best known experimental limits for

CP and CPT violations in the positronium decays which are at the level of 0.3 % [3,4]. Ortho-positronium is symmetric in space and spin and, therefore, as a system built from fermions it must be charge symmetry odd. Para-positronium,

123

445 Page 2 of 14 Eur. Phys. J. C (2016) 76:445

y

x

z

Fig. 1 Photo of the Jagiellonian positron emission tomograph (J-PET). The J-PET detector is made of three cylindrical layers of EJ-230 plastic scintillator strips (black) with dimension of 7 19 500 mm3 and

Hamamatsu R9800 vacuum tube photomultipliers (grey). The signals from photomultipliers are probed in the voltage domain at four thresholds with the timing accuracy of 30 ps [10] and the data acquisition is working in the trigger-less mode [11,12]

in turn, as anti-symmetric in spin and symmetric in space, must be charge symmetry even. C symmetry conservation implies that the ortho-positronium annihilate into odd number of gamma quanta, 3 being the most probable, with lifetime 142 ns and para-positronium decays into even number of gamma quanta with lifetime 125 ps [58]. Such a huge difference in the life-times enables an efcient experimental disentangling of o-Ps from p-Ps decays.

With the recently constructed J-PET detector (see Fig. 1) we intend to study the o-Ps 3 process in order to exam

ine discrete symmetries and to test new medical imaging techniques based on the detection of three photons [9].

In the ortho-positronium decay the additional information carried by the 3rd allows more precise annihilation point reconstruction. Schematic view of p-Ps and o-Ps annihilation is shown in Fig. 2.

Moreover, the observed yield of three gamma annihilation depends on materials properties (see Sect. 2.1), therefore it may allow to gain some information not only about location but also about properties of tumors [13]. In fundamental physics, studies of the three gamma annihilation allows not only to test the discrete symmetry violation [14] but also enables searches of physics beyond the Standard Model: extra dimensions [15], dark matter [16] and a new light vector gauge boson [17]. Since a detailed physics program of J-PET and its motivation is described elsewhere in a dedicated article [1], here as an example we would like only to discuss briey experimental approach to determining the expectation value of the odd operator for the CPT symmetry, whose violation has not been observed so far. As it was recently shown [18] the J-PET detector allows for a spin direction (S) determination of o-Ps created in cylindri-

z

y

x

Fig. 2 Schematic view of a single layer of the J-PET detector with two (up) or three (down) gamma quanta annihilation. In presently built geometry the rst layer consists of 48 plastic scintillators (green bars). In this pictorial representation, for clarity, a smaller number of strips is shown. Solid dark blue lines indicate annihilation quanta and dashed brown line indicates de-excitation gamma quantum e.g. from the

22Na 22Ne +e+ + 22Ne+ +e+ + decay chain. Due to the

momentum conservation annihilation quanta are moving along the same line in the case of e+e 2 , while in the case of the e+e 3

they are included in a single plane. The de-excitation photon (dashed line) is not correlated with the annihilation photons and is isotropically distributed with respect to the annihilation plane-of-response. Due to the fact that annihilation and de-excitation occur in a good approximation at the same place the photons from the e+e 2 form a plane

with the de-excitation photon

cal target. Additionally as it is described in Sect. 5 the J-PET detector enables determination of the momentum vectors of gamma quanta originating from the o-Ps 3 process.

These properties allow for construction of the following operator odd under CPT transformation: S (k1 k2), where k1

and k2 denote momenta of the most and second most energetic quanta, respectively. The non-zero expectation value (indicating violation of CPT symmetry) would manifest itself as an asymmetry between numbers of events with spin direction pointing to opposite sides of the decay plane (k1 k2).

In this paper we focus on the feasibility study of the detection of o-Ps annihilation.

The registration of three-gamma annihilation and conducting of the above mentioned research is possible by the J-PET detector whose novelty lies in application of plastic scintillators instead of crystals [19]. This solution allows to sample fast signals (5 ns) [10,2023] and build more extended geometries, in comparison to commercially used PET detectors [20]. In this paper we study the feasibility of the three

123

Eur. Phys. J. C (2016) 76:445 Page 3 of 14 445

Table 1 Summary of major physical characteristics of beta-plus isotopes useful for PET imaging and positron annihilation lifetime spectroscopy (PALS) investigations. For isotopes that decay into excited

states the properties of emitted gamma quanta are denoted. Data were adapted from [27]

Isotope Half-life + decay E (MeV) Emaxe+ (MeV) Excited nuclei lifetime

Isotopes for PALS and PET imaging22Na 2.6 (years) 22Na 22 Ne + e+ + e + 1.27 0.546 3.63 (ps)

68Ga 67.8 (min) 68Ga 68 Zn + e+ + e + 1.08 0.822 1.57 (ps)

44Sc 4.0 (h) 44Sc 44 Ca + e+ + e + 1.16 1.474 2.61 (ps)

Isotopes for PET imaging68Ga 67.8 (min) 68Ga 68 Zn + e+ + e 1.899

11C 20.4 (min) 11C 11 B + e+ + e 0.961

13N 10.0 (min) 13N 13 C + e+ + e 1.198

15O 2.0 (min) 15O 15 N + e+ + e 1.735

18F 1.8 (h) 18F 18 O + e+ + e 0.634

gamma annihilation measurements using the J-PET detector. To this end we have developed Monte Carlo simulations accounting for:

(i) positron emission and thermalisation in the target material,

(ii) angular and energy distributions of gamma quanta originating from ortho-positronium annihilation,

(iii) Compton interactions of emitted gamma quanta in the detector built from plastic scintillators,

(iv) determination of gamma quanta hit-position and hit-time in the detector with experimentally determined resolutions,

(v) multiple scattering and accidental coincidences,(vi) reconstruction of registered gamma quanta four-momenta,

and used four possible geometrical congurations of the JPET detector.

Section 2 gives a general introduction of positron emission and interaction with matter together with the formation of positronium and the description of ortho-positronium annihilation into three gamma quanta. Possible detector geometries are summarized in Sect. 3. Properties of J-PET detector, comparison between simulated and experimental spectra and the method of background rejection are presented in Sect. 4.Section 5 contains the detector efciency estimation as well as the energy and angular resolutions.

2 Performance assessment: Monte Carlo simulations

The following paragraphs contain the description of Monte Carlo simulations of positrons emitted from + source (22Na) that bind with electron and form positronium. Sim-

ulation takes into account the effects of nite positronium range and non-zero residual momentum of the annihilation positron-electron pair. Special emphasis is put on a proper description of available phase-space of photons from orthopositronium annihilation and their further detection in the J-PET detector that consists of plastic scintillators.

2.1 Positron source and positronium formation

Table 1 summarizes the important characteristics of the isotopes used for different types of imaging as well as in laboratory studies. Those isotopes decay through + transitions emitting a positron that travels through matter, scatters and slows down reaching thermal energies. Then it undergoes free annihilation or forms a positronium [24]. In water at 20 C the positron has about 64 % chance of undergoing free annihilation [25]. The positronium is produced mostly in the ground state forming para-positronium (1S0, p-Ps) or orthopositronium (3S0, o-Ps) with probability of 25 and 75 %, respectively. The annihilation of those states is leading predominantly to an emission of two or three gamma quanta for p-Ps or o-Ps states, respectively. However, the interactions with matter can lead to inversion of the ortho-positronium spin or to the pick-off processes and, as a result, can affect the relative ratio of 3 /2 annihilation. The effective yield of annihilation into 3 in most of non-metallic substances is of the order of 1 %, although in some cases, as for example ne powders of alkaline oxides, it can reach even 29 % as recently shown for the amberlite porous polymer XAD-4 (CAS 37380-42-0) [26].

Some of the + emitters, e.g. 22Na or 44Sc, decay to daughter nucleus in excited states and emit prompt gamma with a well dened energy. In plastic scintillators gamma quanta interact mostly via the Compton scattering. Figure 3 shows the energy loss spectrum expected for the gamma

123

445 Page 4 of 14 Eur. Phys. J. C (2016) 76:445

E

+

= 511 keV (e

e

Events/2 keV

70

3

10

)

-

60

Sc)

Na)

50

40

30

20

10

0 0 200 400 600 800 1000 1200

Edep [keV]

Fig. 3 Simulated spectra of deposited energy in plastic scintillators for gamma quanta from e+e 2 annihilation and for de-excitation

gamma quanta originating from isotopes indicated in the legend. The spectra were simulated including the energy resolution of the J-PET detector [20] and were normalized to the same number of events

quanta from the e+e 2 annihilation compared to the

spectra expected from the de-excitation quanta from 22Na

and 44Sc isotopes.

The results were obtained taking into account the experimental energy resolution of the J-PET detector [28]. The identication of de-excitation and annihilation photons is based on the energy loss and angular correlations. Using the energy loss criterion (e.g. Edep > 0.370 MeV) we can uniquely identify de-excitation quantum from the 44Sc and

22Na decays with a selection efciency of 0.66 and 0.70, respectively. The second selection method is, however, much more efcient. It will be based on the relation between the relative angles of the photons directions. The trilateration method allows reconstruction of an emission point [18] and

the relative angles between the gamma quanta. After assigning the numbers to the photons such that the relative angles are arranged in the ascending order (12 < 23 < 31), in the case of the 2 annihilation (Fig. 2, left) the largest angle 31 will be equal to 180 degrees and will correspond to the photons from the e+e 2 process. Therefore, the de-

excitation gamma quantum can be identied as photon number 2. This second selection method is independent of the energy loss criteria, and due to the high angular resolution of the J-PET tomograph (see Sect. 3), it will allow an identication with close to 100 % selection efciency. In case of the 3 annihilation, photons originating from o-Ps 3 pro

cess are emitted in a single plane. The gammas directions are not correlated with the de-excitation photon (see Fig. 2), so probability of their miss-identication with the de-excitation gamma quanta is at the level of few percent only, and for long lifetime of o-Ps (larger than few ns) it is negligible due to the large difference between the hit-times of annihilation and de-excitation photons which may be used as additional third criterion.

In further considerations we will focus on sodium isotope, which is commonly used as a source of positrons for various experiments and tests of detectors. Pictorial representation of the studied o-Ps 3 process is shown in Fig. 4.

In the conducted simulations we took into account the description of positron properties after thermalisation. Its energy was simulated according to the distribution presented in Fig. 5 [29].

The distribution of the initial positron kinetic energy depends only on thermalisation processes. This distribution is taken into account in the transformation of gamma quanta four-momenta from the rest frame of ortho-positronium to the laboratory frame. In addition, the small distance traveled by positron in matter was taken into account. Positron range depends on material properties and can be generated from proles known in the literature [30] provided by many simulation packages, such as GATE [31] or PeneloPET [32]. In this work the positron range distribution obtained by PeneloPET was adopted. Abovementioned effects introduce

3.7 ps

time

e+

e+ thermalization

oPs formation

oPs lifetime

Fig. 4 Scheme of sodium decay and formation of ortho-positronium

123

Eur. Phys. J. C (2016) 76:445 Page 5 of 14 445

10

3

3

10

Probability

22

Events/5 keV

25

20

18

20

16

14

15

12

10

10

8

6

5

4

2

0 0 10 20 30 40

50 60 70 80 90 100

0 0 100 200

300 400 500

E

e

[eV]

E

Fig. 5 Simulated probability density function of positronium formation as a function of positron energy after thermalisation in the water. The distribution is adapted from reference [29]

additional smearing of o-Ps annihilation position (see Fig. 6) and are included into performed simulations.

2.2 o-Ps 3 process

Positronium is the lightest purely leptonic system, and it can annihilate only into gamma quanta. Those photons are coplanar in the Center of Mass (CM) frame due to the momentum conservation. The cross-section for annihilation with formation of photons having frequencies i can be expressed as [33]:

[keV]

i

Fig. 7 Energy spectrum of photons originating from three-photon annihilation of an electron and a positron

3 =

4e6 vm2e [integraldisplay]

[integraldisplay]

me

0 me1

me

(1 + 2 me)2 2122

d1d2

=

4e6 vm2e

2 9

3 (1)

where me is electron mass, v denotes electron-positron relative velocity, e is the elementary charge. In above formula the conservation of 4-momentum allows to eliminate one of the frequencies (3). Equation 1 results in the characteristic energy distribution of gamma quanta (see Figs. 7, 8).

Ps kinetic energy [eV]

100

90

80

2

10

70

60

50

10

40

30

20

10

1

x

0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7

d [cm]

Fig. 6 Left Scheme of the ortho-positronium annihilation into three gamma quanta i in the detector reference frame. Gamma quanta are not contained in a single plane due to non-zero kinetic energy of the ortho-positronium. In the experiment a plane of response can be determined from gamma quanta interaction position in the scintillators (green

bars). The distance d between plane of response and annihilation vertex gives information about annihilation position uncertainty. Right Distribution of distance d as a function of kinetic energy of ortho-positronium. Taking into account resolution of the J-PET annihilation point reconstruction [18], the uncertainty caused by o-Pss boost is negligible

123

445 Page 6 of 14 Eur. Phys. J. C (2016) 76:445

Fig. 8 Distribution of angles (left) and Dalitz plot (right) of o-Ps 3 annihilation. Boundaries are determined by kinematic constraints

Fig. 9 Transverse view of simulated geometries

3 Simulated geometries

A few possible detector geometries were simulated. They are referred to as J-PET, J-PET+1, J-PET+2 and J-PET-full:

J-PET corresponds to the already built detector [1] with 3 layers of scintillators (from in to out: 48 + 48 + 96 scintillators).

J-PET+1 the J-PET geometry extended by an additional layer lled by 96 scintillators.

J-PET+2 geometry assumes complete fulllment of all available layers in the J-PET detector (48 + 48 + 96 + 96 + 96).

J-PET-full detector with fully coverage of four plastic scintillator layers.

The details of the simulated geometries are presented in Fig. 9 and Table 2. Comparison of the results for all above mentioned options shows the accuracy of o-Ps 3 reg

istration achievable at current J-PET setup and upgrades planned in the next two years, as well as for the J-PET-full detector. In all the cases simulations assume the usage of EJ-230 plastic scintillator strips (dimensions 1.9 cm0.7 cm

50.0 cm), with the longest side of the scintillator arranged along the z axis.

123

Eur. Phys. J. C (2016) 76:445 Page 7 of 14 445

Table 2 Details of simulated layers of the J-PET geometry. J-PET detector has been already built [1]. The mechanical construction for the next phases J-PET+1 and J-PET+2 is also prepared and the hardware upgrade is planned within the next 2 years

Layer number Layer radius with respect to the center of scintillator (cm)

Number of scintillators in the layer

Angular displacement of ni scintillator

z1

z2

J-PET

1 42.50 48 ni 7.5

2 46.75 48 ni 7.5 + 3.75

3 57.50 96 ni 3.75 + 1.875

J-PET+2

1 42.50 48 ni 7.5

2 46.75 48 ni 7.5 + 3.75

3 50.90 96 ni 3.75

4 53.30 96 ni 3.75 + 1.875

5 57.50 96 ni 3.75 + 1.875

J-PET+1

1 42.50 48 ni 7.5

2 46.75 48 ni 7.5 + 3.75

3 53.30 96 ni 3.75 + 1.875

4 57.50 96 ni 3.75 + 1.875

J-PET-full

1 43.0 400 ni 0.9

2 45.0 437 ni 0.82

3 47.0 473 ni 0.76

4 49.0 508 ni 0.71

4 J-PET detector properties

The multipurpose detector (J-PET) constructed at the Jagiellonian University of which novelty lies in using large blocks of plastic scintillators instead of crystals as detectors of annihilation quanta, requires the usage of the time of signals, instead of their amplitude, and allows to obtain time resolution better than 100 ps [28].

4.1 Determination of hit and time position at J-PET

Reconstruction of time and gamma quanta hit position in ith plastic scintillator can be based on the time values (t Ai, t Bi) of scintillation light registration in photomultipliers located at the ends of single plastic scintillator strip. Then the distance ( zi) along the strip between its center and the hit position can be expressed as:

zi =

(t Ai t Bi) v

2 , (2)

Fig. 10 Registration of the signals arrival time on the two ends of a single scintillator (t A1, t B1 and t A2, t B2 for the rst and second strip, respectively) allows to determine the distance from the scintillators centers ( z1,2) and times (t1,2) when gamma quanta interacts with scintillators.

Then the line of response can be determined as well as the displacement of the annihilation position from its center ( LOR)

where v is the light velocity in the plastic scintillator. Based on this information, in case of two-gamma quanta annihilation, the line of response (LOR) and the annihilation position along it can be determined (see Fig. 10).

In case of three-gamma annihilation, the registered gamma quanta are coplanar (o-Ps kinetic energy can be neglected, see Sect. 2.1) and the registered hit-points form the plane-of-response (POR) (see Fig. 2, right panel). In this case the annihilation position can be determined using the novel reconstruction based on trilateration method (see Sect. 5.1.2). The obtained energy and time resolution of registered gamma quanta were experimentally determined and within the range of deposited energy (Edep) (200, 340) keV, are

equal to [28]:

(T 0hit) 80 ps, (3)

(E)

E =

0.44E [MeV]. (4)

For lower energies the time resolution can be expressed as a function of deposited energy (Edep):

(Thit(Edep)) =

(T 0hit)[ps]

Edep[keV]270. (5)

Considering the most challenging time reconstruction for gamma quanta with low energies (around 50 keV), one can see that the J-PET detector provides a precision on the level of two hundred picoseconds. In the commercial PET systems the events with an energy deposition lower than about 400 keV [36,37] are discarded.

123

445 Page 8 of 14 Eur. Phys. J. C (2016) 76:445

6

10

]

-1 0.4

1

[cm

0.35

Events/5 keV

0.8

0.3

0.25

0.6

0.2

0.15

0.4

0.1

0.2

0.05

0 0 200 400 600 800 1000 1200 1400

0 0 50 100 150 200 250 300 350 400

E

[keV]

E

[keV]

dep

Fig. 11 Left dependency of attenuation coefcient on incident gamma quanta energy. Data taken from [34]. Right distribution of energy deposited by gamma quanta in plastic scintillators originating from o-Ps 3 annihilations. The shown spectrum is a convolution of

the energy distribution of gamma quanta from the o-Ps 3 decay

(Fig. 7) and the KleinNishina distribution of kinetic energy of electrons acquired via Compton scattering [35]. Spectrum includes the absorption dependence on the energy (left panel) and the detector energy resolution

4.2 Spectra of deposited energy

The probability of incident gamma quanta registration is a function of the attenuation coefcient and distance that gamma quantum travels through the material. In the simulations the attenuation coefcient was parametrized as a function of incident gamma quanta energy (see Fig. 11, left panel).

Gamma quanta interact with plastic scintillators mainly via Compton effect and the characteristic spectra of deposited energy are described by Klein-Nishina formula [35,38]. The distribution for 511 keV incident gamma quantum is shown in Fig. 12. Energy of single gamma quanta from orthopositronium annihilation is within [0, 511] keV energy range

and the spectrum of deposited energy via Compton effect for the corresponding energy range is presented in Fig. 11 (right panel).

4.3 Background rejection

Direct annihilation of positron with electron, as well as intrinsic annihilation of para-positronium, are both characterized by short times of 400 and 125 ps, respectively. For com

parison, an ortho-positronium lifetime in vacuum amounts to about 142 ns [68]. Therefore, events corresponding to direct annihilation and decay of para-positronium can be reduced to a negligible level by requiring the time difference between de-excitation photon and annihilation photons detection to be larger than e.g. 20 ns. However, such lifetime criterion cannot discriminate pick-off and conversion processes of o-Ps which may lead to the annihilation into 2 quanta.

Annihilation into 2 may mimic a registration of 3 annihilation due to the secondary scatterings in the detector. Such scattering is shown pictorially in Fig. 13. For the reduction of this background the following complementary methods can be considered, based on information of:

relation between position of the individual detectors and the time difference between registered hits,

angular correlation of relative angles between the gamma quanta propagation directions,

Fig. 12 Spectra of simulated (red, dashed line) and measured energy (solid, black line) deposition by 511 keV gamma quanta in J-PET detector. The simulated spectrum was normalized to the experimental one, and simulations were performed taking into account the energy resolution (Eq. 4). The left part of the experimental spectrum was cut due to the triggering threshold applied in the experiment

123

Eur. Phys. J. C (2016) 76:445 Page 9 of 14 445

12

y

y

y

x e+e- 2 single scattered

x e+e- 2 double scattered

x

o-Ps 3

Fig. 13 Pictorial illustration of the possible response of the detector to o-Ps 3 and e+e annihilation into 2 . Arranged circularly squares

represents scintillator stripspurple and green colors indicate strips where the gamma quanta were or were not registered, respectively. The

arrows represents gamma quanta occurring in the events, while dotted lines indicate naively reconstructed gamma quanta. Examples of primary and secondary scatterings are depicted

Fig. 14 Distribution of o-Ps 3 (green) and scattered events

(brown) as a function of 12 vs 23 angles. Events, where one of the gamma from e+e 2 annihilation is registered in the detector

while the other is scattered and cause signals in two detectors, lies on the diagonal of the plot. Events where one gamma is missing detection, and the other undergoes two scatterings are localized below the diagonal line. Example of analysis cut, rejecting 3 % of signal and reducing background by factor 104, is shown as a dashed purple line. Distribution includes the angular resolution of the J-PET detector

the distance between the origin of the annihilation (position of the annihilation chamber) and the decay plane.

In Fig. 14 we show as an example spectra the 23 vs 12 distribution, where i j are the ordered opening angles (12 < 23 < 13) between registered gammas. For the o-Ps 3 process, due to the momentum conservation,

23 > 180 12 and therefore events corresponding to the

o-Ps 3 decay will lie above the diagonal, as shown

in green colour in Fig. 14. Background events will correspond to points at the diagonal (23 = 180 12) and below

diagonal (23 < 180 12) as can be inferred from the

middle and left panel of Fig. 13. Therefore, one of the possible selection cuts can be applied on ordered opening angles

(12 < 23 < 13) between registered gammas, and is resulting in a decrease of background by a factor 104 while, rejecting only 3 % of signal events (see Fig. 14). Combining aforementioned criterion with requirement that registered time difference ( t) as a function of detector number ( I D) is small ( t < 0.3 ns), allows for total reduction of the instrumental background by a factor of 109. However, we have to take into account that the remaining background is caused not only by misidentied 2 events, but also by true annihilations into 3 which may originate from the interaction of the positronium with surrounding electrons and hence will constitute a background for studies of discrete symmetries. Interaction of ortho-positronium with matter is classied into: pick-off annihilations and ortho-para spin conversion. Contribution from these processes depends on the used target material, e.g. in aerogel IC3100 and amberlite porous polymer XAD-4 about 7 and 36 % of ortho-positronium undergo through it, respectively [26]. The events originating from the true of o-Ps 3 annihilation process (NoPs) can be misidenti

ed with the events from the following processes: pick-off process with direct annihilation to 3 (N3 pick-of f ); pick-off process with annihilation to 2 misidentied as 3 due to secondary scatterings (N2 pick-of f ); conversion of orthopositronium to para-positronium with subsequent C symmetry violating decay to 3 (N3 conv); conversion of orthopositronium to para-positronium with subsequent annihilation to 2 misidentied as 3 due to the secondary scatterings (N2 conv).

The conservative upper limit of these background contributions may be estimated as:

N2 conv/NoPs < N2 pick-of f /NoPs

< N3 conv/NoPs < N3 pick-of f /NoPs, (6)

where:

123

445 Page 10 of 14 Eur. Phys. J. C (2016) 76:445

N3 pick-of f /NoPs < (1

matter

vacuum )/370 2

104(IC3100) < 103(XAD-4); N2 pick-of f /NoPs < 0.07 109(IC3100) < 0.36

109(XAD-4); N3 conv/NoPs < 0.072.8106(IC3100) < 0.36

2.8 106(XAD-4);

N2 conv/NoPs < 0.07 109(IC3100) < 0.36

109(XAD-4).

In the above estimations the factor 109 denotes the reduction power of the 2 events and 2.8 106 stands for the

upper limit of the C symmetry violation via the p-Ps 3

process [39]. The precise control of these contributions will be provided by the measurement of the true 2 events with high statistics.

5 J-PET performance in o-Ps 3 decay

measurements

In order to determine the angular and energy resolution we have performed simulations of point-like 22Na source surrounded by water and localized in the geometrical center of the J-PET detector. The conducted simulations accounted for positron emission and thermalisation in the target material, angular and energy distributions of gamma quanta originating from ortho-positronium annihilation and Compton interactions of emitted gamma quanta in the J-PET detector. Details were presented in the Sect. 2. In the next step, based on the simulated data, we reconstructed hit-time and hit-position of the registered gamma quantum interaction in the detector, taking into account the experimentally determined resolutions. Based on obtained informations the reconstruction of angles between gamma quanta and of their energies is performed, as described in the next paragraph.

5.1 Angular and energy resolution

Incident gamma quantum transmits energy as well as momentum to an electron in the plastic scintillator via Compton effect. Due to that, registered signals at the end of the scintillator strips cannot give information about the energy of the incident gamma quantum on the event-by-event basis. However, registration of three gamma quanta hit-position from o-Ps 3 annihilation allows reconstruction of their ener

gies based on the energy and momentum conservation.

In CM frame, energies of three gamma quanta from an ortho-positronium annihilation, can be expressed as a functions of angles (12, 23, 13) between momentum vectors (see also Fig. 8, right panel), as follows:

E1 = 2me

E2 = 2me

cos 12 cos 13cos 23

(1+cos 12)(1+cos 12cos 13 cos 23)

,

E3 = 2me

1+cos 121+cos 12cos 13cos 23

. (7)

The measured positions of gamma interaction in the detector, together with known or reconstructed position of annihilation, allow for Ei determination. The determination of angles requires reconstruction of interaction points and annihilation position. As regards annihilation position we may distinguish two cases, discussed in the next paragraphs.

5.1.1 Point-like positronium source

In some cases of discrete symmetries studies positronium will be produced in the well localized material surrounding the point-like positron source [1]. Assuming that + emitter position corresponds to the ortho-positronium annihilation point, the angles (12, 13 and 23) between gamma quanta can be determined from registered gamma quanta interaction points (rhit) in the detector. Coordinates x and y are determined as the centre of the scintillator strip, and therefore the precision of their determination correspond to the geometrical cross section of the scintillator strip. The z coordinate is determined from signals arrival time to photomultipliers at the ends of scintillator strip, and its uncertainty is equal to about (z) = 0.94 cm [22,23]. Uncertainty of (rhit)

determination gives the main contribution to estimation of angular and energy resolutions. The second order effect is an uncertainty originating from non zero boost and distance traveled by positron in matter.

5.1.2 Spatially extended positronium source

The angles (12, 23, 13) and hence a full kinematics of o-Ps 3 decay can be also reconstructed in the case of

the extended positronium target. For example a target of a cylindrical shape with the diameter of 20 cm was proposed for the production of a linearly polarized positronium [1]. Polarisation can be determined provided that positron emission and positronium formation (approximately the same as annihilation) position are known.

A new reconstruction algorithm that allows reconstruction of ortho-positronium annihilation position for an event by event basis was recently reported [9,18]. The method based on trilateration allows for a simultaneous reconstruction of both location and time of the annihilation based on time and interaction position of gamma quanta in the J-PET detector. The reconstruction performance strongly depends on detector time resolution ((Thit)). Using aforementioned reconstruction algorithm, current J-PET spatial resolution for annihilation point reconstruction is at the level of 1.5 cm along the main detector axis and 2 cm in the transverse plane [18].

cos 13+cos 12 cos 23

(1+cos 12)(1 + cos 12cos 13 cos 23)

,

123

Eur. Phys. J. C (2016) 76:445 Page 11 of 14 445

/ ndf

178.4 / 197 N 123.8

1504

0.01346

0.03927

0.0574

0.8453

193.9

2643

0.00411

0.0102

0.4703

120.0

829.7

0.01683

1.468

0.041

/ ndf

160 / 111 N 9.164e+02

3

10

1.973e+04

0.020

0.116

0.053

Events/0.5 keV

]

Events/0.1 [

2.615N 7.346e+02

3

4.5

3

10

0.2905

0.169

0.0511

4

5.574

3.622e+02

2.5

3.5

1.025

0.43

0.156

2

3

15.77

2.5

1.5

2

1

1.5

1

0.5

0.5

0

0 5 0 5

20 0 20

rec

-

[

true

]

E

rec

-E

[keV]

true

Fig. 15 Resulting angular (left) and energy (right) resolution spectra for point-like positronium source with known location and assumed detector resolution (T 0hit) = 80 ps

9

60

]

) [ (

8

(E) [keV]

50

7

6

40

5

30

4

3

20

2

10

1

0 0 20 40 60 80 100 120 140 160 180 200

0 0 20 40 60 80 100 120 140 160 180 200

(T

0 ) [ps]

hit

(T

0 ) [ps]

hit

Fig. 16 Angular (left) and energy (right) resolution for the registration of the gamma quanta originating from ortho-positronium annihilation as a function of detector time resolution for point-like (blue box) and extended (black triangle) positronium source

5.1.3 Performance studies

The angular and energy resolutions for the registration of the gamma quanta from the o-Ps 3 decay are estab

lished from simulations i.e. the distributions of the differences between generated and reconstructed values of angles and energies. Figure 15 show results obtained under assumption that the hit-time resolution is given by Eqs. 4 and 5.In order to determine the angular and energy resolution the triple Gaussian model, which effectively describes obtained distributions, was applied:

f (x) =

3

[summationdisplay]

i=1

[bracketleftBigg][parenleftBigg]

Ni

3j=1 Nj

[parenrightBigg] i[bracketrightBigg]2. (9)

Since the angular and energy resolution strongly depend on hit-time resolution registered in the J-PET detector, the studies of resolution were made for (T 0hit) in the range from 0 ps to 190 ps. Comparison between obtained resolutions for the point-like and extended positronium source is shown in Fig. 16. In both cases energy and angular resolutions are improving with decreasing (T 0hit), and for presently

Ni 2i

e

1

2

[parenleftBig]

xi

i

2, (8)

where Ni, i and i were varied in the t. The total uncertainty was obtained as a standard deviation of the total distribution equivalent to:

total = [radicaltp]

[radicalvertex]

[radicalvertex]

[radicalbt]

3

[summationdisplay]

i=1

123

445 Page 12 of 14 Eur. Phys. J. C (2016) 76:445

exp> E

i

1 Ee

th

exp > E

i

2 Ee

th

exp > E

i

3 Ee

th

exp> E

i

1 Ee

det

det

th

exp > E

i

2 Ee

th

10

1

10

1

exp > E

i

3 Ee

th

10

1

10

1

10

2

10

2

10

3

10

3

10

4

10

4

10

5

6

10

5

6

7

7

10 0 50 100 150 200 250 300 350 400

10 0 50 100 150 200 250 300 350 400

(a) (b)

E

[keV]

th

E

[keV]

th

exp> E

i

1 Ee

th

exp > E

i

2 Ee

th

exp > E

i

3 Ee

th

exp> E

i

1 Ee

det

det

th

exp > E

i

2 Ee

th

10

1

10

1

3 Ee

exp > E

th

10

1

10

1

10

2

10

2

10

3

10

3

10

4

10

4

10

5

6

10

5

6

7

7

10 0 50 100 150 200 250 300 350 400

10 0 50 100 150 200 250 300 350 400

E

[keV]

th

E

[keV]

th

(c) (d)

Fig. 17 oPs 3 registration efciency (determined taking into

account geometrical acceptance, probability of gamma quanta registration in the plastic scintillator and J-PET detector resolution) as a function of applied threshold for different types of simulated geometries.

The shown dotted, dashed and solid lines indicate efciency assuming that at least one, two or three photons deposited energy above the threshold, respectively

achieved time resolution of (T 0hit), and well a localized point-like positronium source, they amount to () = 0.4

and (Ehit) = 4.1 keV, respectively. In case of the extended

positronium source, when the reconstruction of the annihilation point is needed both resolutions increases to () = 4.2

and (Ehit) = 30 keV, respectively.

5.2 J-PET efciency studies with Monte Carlo simulations

The rate of registered o-Ps 3 events in general can be

expressed by the formula:

RoPs3 = A foPs3 det(th) ana, (10)

where A is the total annihilation rate (fast timing of applied plastic scintillators allows for usage of the 10 MBq positron source), foPs3 is the fraction of annihilations via o-Ps

3 process in the target material, det(th) is the detector efciency as a function of applied detection threshold while ana denotes selection efciency used to discriminate between 3 and 2 events.

The det efciency of the o-Ps 3 reconstruction will

depend on the energy deposition threshold used in the analysis (see Fig. 12). The hardware threshold at the order of 10 keV [28] will be set to discriminate the experimental noise and later on we will apply further selection threshold based on the measured energy deposition. The probability of registration of 1, 2 or 3 gamma quanta originating from o-Ps 3

123

Eur. Phys. J. C (2016) 76:445 Page 13 of 14 445

Table 3 Expected rate of registered signal events in different geometries and target materials assuming 106 annihilations per second and requiring energy deposition above 50 keV for all three gamma quanta from o-Ps 3 decay

Target material Rate of registered o-Ps 3 events (s1)J-PET J-PET+1 J-PET+2 J-PET-full

IC3100 15 70 130 10600

XAD-4 25 115 230 18300

annihilation ( det) as a function of applied selection threshold in different geometries is shown in Fig. 17. Efciency det contains contribution from geometrical acceptance, probabilities of gamma quanta interaction in applied plastic scintillators and it was determined taking into account the J-PET detector resolution. In our evaluation we assume conservatively that the event selection threshold will be set to 50 keV. A fraction of annihilations via o-Ps 3 process is esti

mated taking into account only longest lived component in two selected materials IC3100 ( foPs3 = 16.6 %) and

XAD-4 ( foPs3 = 28.6 %) [26]. The expected rate of reg

istered signal events is shown in Table 3. Using in the experiment amberlite porous polymer XAD-4 instead of aerogel IC3100 as target material, allows to collect the required statistics almost twice faster, however, resulting with higher systematic uncertainties due to the interaction of positronium with the target material, as discussed in Sect. 4.3.

6 Conclusions

We presented results of Monte Carlo simulations showing that the Jagiellonian-PET multipurpose detector constructed at the Jagiellonian University allows exclusive registration of the decays of ortho-positronium into three photons (o-Ps

3 ) providing angular and energy resolution of ()

0.4 and (E) 4.1 keV, respectively.

The achieved results indicate that the J-PET detector gives a realistic chance to improve the best present limits established for the CP and CPT symmetry violations in the decays of positronium [3,4] by more than an order of magnitude. This can be achieved by (1) collecting at least two orders of magnitude higher statistics, due to the possibility of using a + source with higher rate (10 MBq at J-PET vs 0.37 MBq at Gammasphere [3] or 1 MBq at Tokyo University experiment [4]), (2) the enhanced fraction of 3 events by the use of the amberlite polymer XAD-4, (3) a measurements with a few times improved angular resolution and (4) about two times higher degree of o-Ps polarization, as shown recently in reference [18]. The limitation on the source activity can be overcome by the J-PET due to the application of plastic scintillators that are characterized by about two orders of

magnitude shorter duration of signals, thus decreasing significantly the pile-ups problems with respect to the crystal based detector systems. In addition, the improved angular resolution combined with the superior timing of the J-PET detector (by more than order of magnitude improved with respect to the crystal detectors) and with the possibility of the triggerless registrations [11,12] of all kind of events with no hardware coincidence window allow suppression and monitoring of the background, due to misidentication of 2 events and possible contribution from 3 pick-off annihilations.

Acknowledgments We acknowledge valuable discussions with Dr.J. Wawryszczuk and technical and administrative support by A. Heczko, M. Kajetanowicz, W. Migda, and the nancial support by the Polish National Center for Research and Development through Grants INNOTECH-K1/IN1/64/159174/NCBR/12 and LIDER-274/L-6/14/NCBR/2015, the Foundation for Polish Science through MPD program and the EU, MSHE Grant No. POIG .02.03.00-161 00-013/09, Marian Smoluchowski Krakw Research Consortium MatterEnergy Future, and the Polish Ministry of Science and Higher Education through Grant 7150/E-338/M/2015. BCH gratefully acknowledges the Austrian Science Fund FWF-23627.

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecomm http://creativecommons.org/licenses/by/4.0/

Web End =ons.org/licenses/by/4.0/ ), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Funded by SCOAP3.

References

1. P. Moskal et al., Acta Phys. Polon. B 47, 509 (2016)2. W. Bernreuther, U. Low, J.P. Ma, O. Nachtmann, Z. Phys. C 41, 143 (1988)

3. P.A. Vetter, S.J. Freedman, Phys. Rev. Lett. B 91, 263401 (2003)4. T. Yamazaki, T. Namba, S. Asai, T. Kobayashi, Phys. Rev. Lett. 104, 083401 (2010)

5. M.D. Harpen, Med. Phys. 31, 57 (2004)6. A.H. Al-Ramadhan, D.W. Gidley, Phys. Rev. Lett. 72, 1632 (1994)7. R.S. Vallery, P.W. Zitzewitz, D.W. Gidley, Phys. Rev. Lett. 90, 203402 (2003)

8. O. Jinnouchi, S. Asai, T. Kobayashi, Phys. Lett. B 572, 117 (2003)

9. A. Gajos, E. Czerwiski, D. Kamiska, P. Moskal, Multi-tracer morphometric image reconstruction method and system using fast analytical algorithm for calculating time and position of positron-electron annihilation into three gamma quanta. PCT/PL2015/050038

10. M. Paka et al., Bioalgorithms Med. Syst. 10(1), 41 (2014)11. G. Korcyl et al., Bioalgorithms Med. Syst. 10(1), 37 (2014)12. G. Korcyl et al., Acta Phys. Polon. B 47, 491 (2016)13. K. Kacperski, N.M. Spyrou, F.A. Smith, IEEE Trans. Med. Im. 23, 525 (2004)

14. D. Kamiska et al., Nukleonika 60(4), 729 (2015)15. S.L. Dubovsky, V.A. Rubakov, P.G. Tinyakov, Phys. Rev. D 62, 105011 (2000)

16. P. Crivelli, A. Belov, U. Gendotti, S. Gninenko, A. Rubbia, JINST 5, P08001 (2010)

17. S.N. Gninenko, N.V. Krasnikov, A. Rubbia, Mod. Phys. Lett. A 17, 1713 (2002)

123

445 Page 14 of 14 Eur. Phys. J. C (2016) 76:445

18. A. Gajos et al., Nucl. Instrum. Methods A 819, 54 (2016)19. P. Moskal, Strip device and the method for the determination of the place and response time of the gamma quanta and the application of the device for the positron emission thomography. 2011008119, US2012112079, PL388555, JP2012533734, EP2454612

20. P. Moskal et al., Phys. Med. Biol. 61, 2025 (2016)21. P. Moskal et al., Nucl. Instrum. Methods A 775, 54 (2015)22. L. Raczyski et al., Nucl. Instrum. Methods A 786, 105 (2015)23. L. Raczyski et al., Nucl. Instrum. Methods A 764, 186 (2014)24. Y.C. Jean, P.E. Mallon, D.M. Schrader, Positron and Positronium Chemistry (World Scientic Publishing, Singapore, 2003)

25. P. Colombino, B. Fiscella, Nuovo Cimento 3, 1 (1971)26. B. Jasiska et al., Acta Phys. Polon. B 47, 453 (2016)27. National Nuclear Data Center. http://www.nndc.bnl.gov/

Web End =http://www.nndc.bnl.gov/ . Access:20.01.2016r28. P. Moskal et al., Nucl. Instrum. Methods A 764, 317 (2014)29. Ch. Champion, Braz. Arch. Biol. Technol. 48, 191 (2005)30. J. Cal-Gonzalez, J.L. Herraiz, S. Espana, M. Desco, J.M. Udias, Phys. Med. Biol. 58, 5127 (2013)

31. S. Jan, G. Santin, D. Strul, S. Staelens, K. Assi, D. Autret, S. Avner, Phys. Med. Biol. 49, 4543 (2004)

32. J. Espaa, L. Herraiz, E. Vicente, J.J. Vaquero, M. Desco, J.M. Udas, Phys. Med. Biol. 54, 1723 (2009)

33. V.B. Berestetskii, E.M. Lifshitz, L.P. Pitaevskii, Relativistic Quantum Theory (Headington Hill Hall, Pergamon Press, Oxford, 1971)

34. Saint-Gobain Crystals, Organic Scintillation Materials [Brochure] (Saint-Gobain Crystals, Hiram, 2013)

35. C. Leroy, P. Rancoita, Principles of radiation interaction in matter and detection (World Scientic, Singapore, 2009)

36. V. Bettinardi et al., Med. Phys. 38(10), 5394 (2011)37. S. Surti et al., J. Nucl. Med. 48, 471 (2007)38. O. Klein, T. Nishina, Z. Phys. 52, 853 (1929)39. A.P. Mills, S. Berko, Phys. Rev. Lett. 18, 420 (1967)

123