Ortho-positronium observation in the Double Chooz experiment

Abstract

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

Abstract

The Double Chooz experiment measures the neutrino mixing angle [theta] ^sub 13^ by detecting reactor ...... via inverse beta decay. The positron-neutron space and time coincidence allows for a sizable background rejection, nonetheless liquid scintillator detectors would profit from a positron/electron discrimination, if feasible in large detector, to suppress the remaining background. Standard particle identification, based on particle dependent time profile of photon emission in liquid scintillator, can not be used given the identical mass of the two particles. However, the positron annihilation is sometimes delayed by the ortho-positronium (o-Ps) metastable state formation, which induces a pulse shape distortion that could be used for positron identification. In this paper we report on the first observation of positronium formation in a large liquid scintillator detector based on pulse shape analysis of single events. The o-Ps formation fraction and its lifetime were measured, finding the values of 44 % ±12 % (sys.) ±5 % (stat.) and 3.68 ns ±0.17 ns (sys.) ±0.15 ns (stat.) respectively, in agreement with the results obtained with a dedicated positron annihilation lifetime spectroscopy setup.

Published for SISSA by Springer

Received: July 28, 2014

Revised: September 12, 2014

Accepted: September 14, 2014 Published: October 6, 2014

Y. Abe,aa J.C. dos Anjos,e J.C. Barriere,n E. Baussan,v I. Bekman,a M. Bergevin,i T.J.C. Bezerra,y L. Bezrukov,m E. Blucher,f C. Buck,s J. Busenitz,b A. Cabrera,dE. Caden,j L. Camilleri,h R. Carr,h M. Cerrada,g P.-J. Chang,o E. Chauveau,yP. Chimenti,ae A.P. Collin,s E. Conover,f J.M. Conrad,r J.I. Crespo-Anadn,gK. Crum,f A.S. Cucoanes,w E. Damon,j J.V. Dawson,d J. Dhooghe,i D. Dietrich,adZ. Djurcic,c M. Dracos,v M. Elnimr,b A. Etenko,q M. Fallot,w F. von Feilitzsch,acJ. Felde,i,1 S.M. Fernandes,b V. Fischer,n D. Franco,d M. Franke,ac H. Furuta,yI. Gil-Botella,g L. Giot,w M. Gger-Ne ,ac L.F.G. Gonzalez,af L. Goodenough,c M.C. Goodman,c C. Grant,i N. Haag,ac T. Hara,p J. Haser,s M. Hofmann,ac G.A. Horton-Smith,o A. Hourlier,d M. Ishitsuka,aa J. Jochum,ad C. Jollet,v,2F. Kaether,s L.N. Kalousis,ag Y. Kamyshkov,x D.M. Kaplan,l T. Kawasaki,tE. Kemp,af H. de Kerret,d D. Kryn,d M. Kuze,aa T. Lachenmaier,ad C.E. Lane,jT. Lasserre,n,d A. Letourneau,n D. Lhuillier,n H.P. Lima Jr,e M. Lindner,sJ.M. Lpez-Castao,g J.M. LoSecco,u B. Lubsandorzhiev,m S. Lucht,a J. Maeda,ab,3C. Mariani,ag J. Maricic,j,4 J. Martino,w T. Matsubara,ab G. Mention,nA. Meregaglia,v T. Miletic,j R. Milincic,j,4 A. Minotti,v Y. Nagasaka,k Y. Nikitenko,mP. Novella,d L. Oberauer,ac M. Obolensky,d A. Onillon,w A. Osborn,x C. Palomares,gI.M. Pepe,e S. Perasso,d P. Pfahler,ac A. Porta,w G. Pronost,w J. Reichenbacher,bB. Reinhold,s,4 M. Rhling,ad R. Roncin,d S. Roth,a B. Rybolt,x Y. Sakamoto,zR. Santorelli,g A.C. Schilithz,e S. Schnert,ac S. Schoppmann,a M.H. Shaevitz,hR. Sharankova,aa S. Shimojima,ab D. Shrestha,o V. Sibille,n V. Sinev,mM. Skorokhvatov,q E. Smith,j J. Spitz,r A. Stahl,a I. Stancu,b L.F.F. Stokes,adM. Strait,f A. Stken,a F. Suekane,y S. Sukhotin,q T. Sumiyoshi,ab Y. Sun,b,4R. Svoboda,i K. Terao,r A. Tonazzo,d H.H. Trinh Thi,ac G. Valdiviesso,e

1Now at Department of Physics, University of Maryland, College Park, Maryland 20742, U.S.A.

2Corresponding author.

3Now at Department of Physics, Kobe University, Kobe, 657-8501, Japan.

4Now at Department of Physics & Astronomy, University of Hawaii at Manoa, Honolulu, Hawaii 96822, U.S.A.

Open Access, c

[circlecopyrt] The Authors.

Article funded by SCOAP3. doi:http://dx.doi.org/10.1007/JHEP10(2014)032

Web End =10.1007/JHEP10(2014)032

Ortho-positronium observation in the Double Chooz experiment

JHEP10(2014)032

The Double Chooz collaboration

N. Vassilopoulos,v C. Veyssiere,n M. Vivier,n S. Wagner,s N. Walsh,i H. Watanabe,sC. Wiebusch,a L. Winslow,r M. Wurm,ad,5 G. Yang,c,l F. Yermiaw and V. Zimmerac

aIII. Physikalisches Institut, RWTH Aachen University, 52056 Aachen, Germany

bDepartment of Physics and Astronomy, University of Alabama, Tuscaloosa, Alabama 35487, U.S.A.

cArgonne National Laboratory, Argonne, Illinois 60439, U.S.A.

dAstroParticule et Cosmologie, Universit Paris Diderot, CNRS/IN2P3, CEA/IRFU, Observatoire de Paris, Sorbonne Paris Cit, 75205 Paris Cedex 13, France

eCentro Brasileiro de Pesquisas Fsicas, Rio de Janeiro, RJ, 22290-180, Brazil

f The Enrico Fermi Institute, The University of Chicago, Chicago, Illinois 60637, U.S.A.

gCentro de Investigaciones Energticas, Medioambientales y Tecnolgicas, CIEMAT, 28040, Madrid, Spain

hColumbia University; New York, New York 10027, U.S.A.

iUniversity of California, Davis, California 95616, U.S.A.

jDepartment of Physics, Drexel University, Philadelphia, Pennsylvania 19104, U.S.A.

kHiroshima Institute of Technology, Hiroshima, 731-5193, Japan

lDepartment of Physics, Illinois Institute of Technology, Chicago, Illinois 60616, U.S.A.

mInstitute of Nuclear Research of the Russian Academy of Sciences, Moscow 117312, Russia

nCommissariat lEnergie Atomique et aux Energies Alternatives, Centre de Saclay, IRFU, 91191 Gif-sur-Yvette, France

oDepartment of Physics, Kansas State University, Manhattan, Kansas 66506, U.S.A.

pDepartment of Physics, Kobe University, Kobe, 657-8501, Japan

qNRC Kurchatov Institute, Moscow 123182, Russia

rMassachusetts Institute of Technology, Cambridge, Massachusetts 02139, U.S.A.

sMax-Planck-Institut fr Kernphysik, 69117 Heidelberg, Germany

tDepartment of Physics, Niigata University, Niigata, 950-2181, Japan

uUniversity of Notre Dame, Notre Dame, Indiana 46556, U.S.A.

vIPHC, Universit de Strasbourg, CNRS/IN2P3, 67037 Strasbourg, France

wSUBATECH, CNRS/IN2P3, Universit de Nantes, Ecole des Mines de Nantes,44307 Nantes, France

xDepartment of Physics and Astronomy, University of Tennessee,Knoxville, Tennessee 37996, U.S.A.

yResearch Center for Neutrino Science, Tohoku University, Sendai 980-8578, Japan

zTohoku Gakuin University, Sendai, 981-3193, Japan

aaDepartment of Physics, Tokyo Institute of Technology, Tokyo, 152-8551, Japan

abDepartment of Physics, Tokyo Metropolitan University, Tokyo, 192-0397, Japan

acPhysik Department, Technische Universitat Mnchen, 85748 Garching, Germany

adKepler Center for Astro and Particle Physics, Universitat Tbingen, 72076 Tbingen, Germany

aeUniversidade Federal do ABC, UFABC, Santo Andr, SP, 09210-580, Brazil

af Universidade Estadual de Campinas-UNICAMP, Campinas, SP, 13083-970, Brazil

agCenter for Neutrino Physics, Virginia Tech, Blacksburg, Virginia 24061, U.S.A.E-mail: mailto:[email protected]

Web End [email protected] , mailto:[email protected]

Web End [email protected]

5Now at Institut fr Physik and Excellence Cluster PRISMA, Johannes Gutenberg-Universitat Mainz, 55128 Mainz, Germany.

JHEP10(2014)032

Abstract: The Double Chooz experiment measures the neutrino mixing angle 13 by de

tecting reactor

e via inverse beta decay. The positron-neutron space and time coincidence allows for a sizable background rejection, nonetheless liquid scintillator detectors would prot from a positron/electron discrimination, if feasible in large detector, to suppress the remaining background. Standard particle identication, based on particle dependent time prole of photon emission in liquid scintillator, can not be used given the identical mass of the two particles. However, the positron annihilation is sometimes delayed by the orthopositronium (o-Ps) metastable state formation, which induces a pulse shape distortion that could be used for positron identication. In this paper we report on the rst observation of positronium formation in a large liquid scintillator detector based on pulse shape analysis of single events. The o-Ps formation fraction and its lifetime were measured, nding the values of 44 % 12 % (sys.) 5 % (stat.) and 3.68 ns 0.17 ns (sys.) 0.15 ns (stat.) respectively, in agreement with the results obtained with a dedicated positron annihilation lifetime spectroscopy setup.

Keywords: Neutrino Detectors and Telescopes

ArXiv ePrint: 1407.6913

JHEP10(2014)032

Contents

1 Introduction 1

2 Ortho-positronium properties in Double Chooz scintillators 2

3 Pulse shape reconstruction in Double Chooz 4

4 Ortho-positronium tagging algorithm 5

5 Results on the ortho-positronium properties 8

6 Possible physics impact 11

7 Conclusions 12

1 Introduction

Recent results of anti-neutrino reactor experiments, i.e. Double Chooz [13], Daya Bay [4] and RENO [5], clearly proved the non-zero value of the 13 mixing angle. Anti-neutrino detection, based on the inverse beta decay (IBD) process (i.e.

e + p ! e+ + n), allowed a

The space and time correlation between the prompt signal given by the positron, and the delayed one given by the neutron absorption on hydrogen or gadolinium, allows for a clear signal signature and strong background reduction. Nonetheless correlated events due to fast neutrons or cosmogenic generated radio-nuclides such as 9Li or 8He still remain as a background.

A technique used in liquid scintillator detectors to extract the signal from background is pulse shape discrimination (PSD) (see ref. [6] and references therein). Di erent particles have a di erent energy loss while crossing the scintillator media, resulting in a particle-dependent time prole of photon emission. This technique is quite e ective in separating light particles (e.g. electrons or positrons) from heavy ones (e.g. protons or alphas), however it is not adequate to distinguish between particles with similar energy losses such as positrons and electrons.

An alternative PSD, based on the observation of the ortho-positronium (o-Ps) was proposed [7] for +/ discrimination. The positron emitted in the IBD process annihilates with an electron in matter, sometimes forming a positronium metastable state, which leads to a delayed annihilation. o-Ps is the positronium triplet state which decays into three s with a lifetime in vacuum of 142 ns. In matter, however, the o-Ps lifetime can be quenched by several factors, such as chemical reactions (oxidation or compound formation), magnetic e ects (spin-ip), or by positronium interactions with the surrounding electrons (pick-o ),

1

JHEP10(2014)032

measurement of sin2(213) of 0.1.

yielding a two- decay [8]. Positronium can also be formed as a singlet state, which decays into two s, called para-positronium (p-Ps), however its lifetime of 125 ps is too short to exploit it for particle discrimination. The o-Ps lifetime has been measured in the most commonly used solvents for organic liquid scintillators in neutrino physics [7, 9], and as a function of the dopants concentration [10]. Although the o-Ps lifetime is typically shortened down to a few nanoseconds, the distortion in the photon emission time distribution is still observable. The Borexino collaboration [11] has exploited the signature provided by the o-Ps induced pulse shape distortion to statistically identify and reject cosmogenic 11C +

decays, the dominant background in the solar pep neutrino rate measurement.

In this paper we report the observation of o-Ps in Double Chooz, which is performed, for the rst time in a large liquid scintillator detector, on pulse shape analysis of single events. The identication of the o-Ps could be used in the future as additional handle in the signal selection, reducing the cosmogenic background due to 9Li or uncorrelated background due to accidentals (+ + n vs. + n chain). Our results fully demonstrate the proof of principle and the capability of this technique, which could be a key point for the development of scintillators for future neutrino experiments.

The analysis is performed on the data set published in ref. [12], corresponding to a live-time of about 228 days.

2 Ortho-positronium properties in Double Chooz scintillators

For the Double Chooz detector [12] two di erent scintillators are used: one for the Target and one for the Gamma Catcher. The liquid scintillator used for the Target is a mixture of n-dodecane, PXE, PPO, bis-MSB and 1 g gadolinium/l as a beta-diketonate complex. The Target volume is surrounded by the Gamma Catcher scintillator, which is similar to the Target (mineral oil, n-dodecane, PXE, PPO, bis-MSB) but Gd-free [13, 14]. The light yield of the Gamma Catcher was chosen to provide identical photoelectron yield across these two layers.

To measure the formation fraction of o-Ps and its lifetime in samples of the Double Chooz scintillators, a standard PALS (Positron Annihilation Lifetime Spectroscopy) system made of two plastic scintillator (BaF2) detectors has been used. The apparatus, at IPHC Strasbourg, is very similar to the one described in detail in ref. [10]. The positron source is a 1 MBq 22Na source, inserted between two 12.7 m titanium layers, immersed in the liquid scintillator. The 1.27 MeV -ray emitted in association with the positron is detected by one detector (lower threshold at 950 keV) and used as trigger, whereas the second detector (400700 keV energy range) is dedicated to the measurement of one -ray of 511 keV coming from the positron annihilation or one from the o-Ps decay; the time between the two signals is measured to reconstruct the o-Ps lifetime. The overall time resolution of the apparatus is 180 ps.

The measured time distribution is tted with a combination of three exponentials (Ai et/i) and a constant C, all convoluted with a Gaussian spread to model the detector

2

JHEP10(2014)032

Events / (0.05 ns)

A

= 0.4126

0.0033

A

1 = 0.2328

0.0031

A

2 = 0.32731

3

0.00095

t

= 0.1477

0.0014

t

1 = 0.4394

0.0052

t

2 = 3.414

0.012

3

10

10

10

JHEP10(2014)032

5 10 15 20 25

Time (ns)

Events / (0.05 ns)

A

= 0.3765

0.0032

A

1 = 0.2832

0.0031

A

2 = 0.31327

3

0.00098

t

= 0.1360

0.0019

t

1 = 0.4273

0.0040

t

2 = 3.446

0.014

3

10

10

10

5 10 15 20 25

Time (ns)

Figure 1. Time distribution (black dots) and corresponding t (red line) for the Target (top) and

Gamma Catcher (bottom) scintillator.

resolution:

F (t) =

3

Xi=1Ai et/i + C. (2.1)

The two short exponentials correspond to the positron annihilation and p-Ps formation and decay in the source structure and in the liquid, whereas the long one corresponds to the o-Ps lifetime. The constant function is used to take into account pile-up and random noise events. The details of the analysis can be found in ref. [10]. The obtained time distributions and the corresponding ts are shown in gure 1 for Target and Gamma Catcher scintillator. The lifetime can be taken immediately from the t parameter 3 whereas the o-Ps formation fraction has to be computed renormalizing correctly the number of events, neglecting the events in the titanium source support, and taking into account the di erent e ciencies for the three and two- decay modes.

Once the systematics are considered (they account typically for about 1.7 % on the formation fraction and for about 1 % on the lifetime [10]) we obtain for the Gamma Catcher

3

Scintillator o-Ps formation fraction o-Ps lifetime

Target 47.6 1.3 % 3.42 0.03 ns Gamma Catcher 45.6 1.3 % 3.45 0.03 ns

Table 1. o-Ps formation fraction and lifetime in Double Chooz scintillator of Target and Gamma Catcher measured with a dedicated PALS setup [10].

scintillator an o-Ps fraction formation of 45.61.3 % and lifetime of 3.450.03 ns. Slightly di erent values were found for the Target: a formation fraction of 47.61.3 % and a lifetime of 3.42 0.03 ns. All the results are summarized in table 1.

Although Target and Gamma Catcher yield values compatible within the errors for both o-Ps formation and lifetime, high precision measurements on di erent liquid scintillators showed that o-Ps formation has a clear dependence on the scintillator loading whereas the lifetime is much more stable [10].

3 Pulse shape reconstruction in Double Chooz

In Double Chooz, the scintillation signal is recorded by 390 10-inch PMTs, installed on the inner wall of the stainless steel bu er tank (see ref. [12] for details). At each trigger a waveform of 256 ns is recorded for each PMT using ash-ADC (FADC) [15]. The PMT baseline is computed using the rst thousands events of the run, selecting the ones in which the PMT itself did not record any signal above threshold (the threshold was set to 2 FADC counts which corresponds to approximately 0.3 photo-electrons). When a signal is observed on a PMT, a linear t is performed on its rising edge (see gure 2). The intercept of the t line with the baseline provides the pulse starting time and represents one entry of the time prole distribution. The obtained event time prole is the distribution of the arrival time of the pulses recorded by all PMTs.

To correctly build the time prole distribution, the time of ight between the reconstructed vertex position and the PMT is subtracted for each pulse. In addition, the calibration of time o set for each channel, as measured with a laser and monitored with a LED system, is accounted for.

Note that the same PMT can of course record more than one pulse in the same event. The minimum time between two pulses on a single PMT was set to 25 ns (i.e. all that happens within this time window is merged into a single pulse), however tests with a narrower time window set to 10 ns showed no signicant inuence on the result on this parameter.

Once all the time pulses in one event are computed, they are sorted and shifted in order to have the rst one equal to zero. This shift is needed to correctly compare time prole of di erent events.

This procedure was carried out on 60Co, 137Cs and 68Ge calibration runs, with the sources located at the center of the detector, and the resulting time prole of all the events are shown in gure 3. The small di erence between the obtained time proles can be interpreted as a slight energy dependence. This was conrmed looking at the changes of

4

JHEP10(2014)032

FADC count

212

0 50 100 150 200 250

211

210

209

208

207

206

205

Time (ns)

Figure 2. Example of the time determination of the pulse for one PMT. The green solid line represents the pedestal value (FADC=210), the blue dashed line shows the t of the pulse edge. The pulse starting time (84 ns) is represented by a red star.

Arbitrary scale

0.045

JHEP10(2014)032

0.035

60 Co source

68 Ge source

137 Cs source

0.04

0.03

0.025

0.02

0.015

0.01

0.005

0 0 5 10 15 20 25 30 35 40

Time (ns)

Figure 3. Time prole distribution for 60Co (thick black line), 137Cs (dashed blue line) and 68Ge (thin red line) events in the center of the Target. The distributions are normalized to one.

the reconstructed time prole for neutrino candidates when computed in di erent ranges of visible energy.

4 Ortho-positronium tagging algorithm

All the prompt signals of IBD events are composed of an e+ ionization signal, followed by the two 511 keV- rays emission. If o-Ps is not formed, the time between these two processes is too short to distinguish them, however, in case of o-Ps formation the delay between the two signals is no longer negligible since o-Ps has a lifetime of about 3.4 ns in the Double Chooz liquid scintillators.

5

The discrimination between ionization and annihilation signals becomes in principle possible, but it is not trivial given the scintillator fast decay time of 2.6 ns in the Target and5.4 ns in the Gamma Catcher (GC) [13]. In addition, faster electronics than used which has a time sample width of 2 ns, would be desirable to provide a better time resolution.

The idea is to look for a double peak pattern in the time prole: the rst one is due to positron ionization, and the second one to the emission of the two 511 keV- rays. Such a pattern can be better observed for events in which the o-Ps has a long lifetime (large time interval between the peaks on the time prole distribution) and for low energy events (second peak less hidden by the tail of the rst signal).

Based on this idea a specic algorithm was developed: a t function was built combining two reference time proles separated by a delay corresponding to the time taken by o-Ps to decay. The parameters considered and their allowed range are quoted here below.

t: time interval between the two reference time proles. It is allowed to vary between 0 (no o-Ps formation observed) and the shortest time for which the number of pulses in the following 50 ns is lower than 70 % of the number of pulses corresponding to the 1.022 MeV signal.

1, 2: normalization of the two reference time proles. The normalization evaluation is based on the visible energy: the second peak corresponds to a signal of 1.022 MeV and the rst one to the remaining energy. The relative normalization is computed assuming that the number of pulses is proportional to the visible energy. This is not completely correct, in particular as the energy increases, since a pulse can correspond to more than a single photo-electron. To overcome this issue, each normalization was allowed to vary in a range of = 0.6 , and the uncertainty of 0.2 was used in the 2 computation. The number of multiple photo-electrons pulses could be estimated from the visible energy reducing the error on the normalization, however this would rely on the vertex position introducing additional systematics related to the vertex reconstruction.

: shift of the rst reference time prole. Since not all the event time proles start exactly at zero, a possible shift is taken into account in a range between 10 ns and the time of the rst pulse in the event time prole.

If fTP is the function representing the reference time prole, each event time prole is therefore tted using the following function ft:

ft(t) = (1 + 1)fTP(t ) + (2 + 2)fTP(t t). (4.1)

The number of pulses per event is of the order of 400, therefore the number of entries per bin (300 bins of 0.5 ns) is rather small and it is Poisson distributed. The following 2 denition [16] was therefore used:

2 = 2

6

JHEP10(2014)032

N

Xi=1

i ni + ni lnni i

[bracketrightbigg] +

2

Xj=1 2j 2j

(4.2)

Entries/ns

18

16

Prompt signal fit Delayed signal fit Global fitt = 0.2 ns

D

Fitted

14

12

10

8

JHEP10(2014)032

6

4

2

0 0 10 20 30 40 50 60

Time (ns)

Entries/ns

12

Prompt signal fit Delayed signal fit Global fitt = 8.4 ns

D

Fitted

14

10

8

6

4

2

0 0 10 20 30 40 50 60

Time (ns)

Entries/ns

9

Prompt signal fit Delayed signal fit Global fitt = 16.0 ns

D

Fitted

8

7

6

5

4

3

2

1

0 0 10 20 30 40 50 60

Time (ns)

Figure 4. Three examples of o-Ps t. The dashed blue line represents the t of the rst time prole, the thin red line the t of the second one and the thick black line is the total t.

7

where N is the number of bins (i.e. 300), ni and i are the number of pulses observed and expected in the ith bin respectively, and only bins for which ni > 0 are considered. The term j with j = 1, 2 represents the normalization variations of the two reference time proles with respect to the computed ones, and j the 20 % normalization error. The 2 is computed using the MINUIT minimization embedded in the ROOT package [17].

To discard from the analysis events for which the t did not converge properly, an upper limit of 2 was applied on the reduced 2. For the same reason, an additional constraint was applied on the ratio between the integral of the tted function and the number of pulses, in the rst and second peak region separately. Events with an integral ratio larger than 2 in either the rst or the second peak region (i.e. the interval between 3 ns before the peak and 10 ns after it) are discarded. The applied cuts removed 0.2 % of the signal and resulted in no bias in the event selection.

As examples, the result of the t obtained for events with t of 0.2 ns, 8.4 ns and16.0 ns are shown in gure 4.

5 Results on the ortho-positronium properties

To unambiguously demonstrate the capability of Double Chooz to observe the o-Ps formation, the algorithm is applied to a pure 60Co sample and compared to the results obtained for the neutrino candidates. The 60Co reference time prole obtained locating the radioactive source at the center of the detector was used for this analysis. Systematic uncertainty related to the choice of reference time prole will be evaluated later on in this section.

Since an increase of the error of the reconstructed vertex is known to rapidly degrade the time prole, a maximal distance of 20 cm between the position where the source was deployed and the reconstructed vertex was required for a clean sample selection. Such a selection could only be applied on the source calibration data used to build the reference time prole.

As far as the neutrino candidate selection is concerned, an additional requirement on the energy with respect to the selection cut of ref. [12] is applied: only events with a visible energy between 1.2 and 3 MeV are analyzed. Below 1.2 MeV the rst peak energy is too small (i.e. below 200 keV) for a correct t convergence, whereas above 3 MeV the second peak is typically hidden by the rst peaks tail and therefore di cult to identify.

As can be seen in gure 5 (red squares), the t distribution for the 60Co events is, as expected, peaked at zero. Nonetheless, a smearing can be seen resulting in t values up to 7 ns. This smearing shows the present limitation of the developed algorithm, which sometimes tends to increase the t parameter reducing the value of the shift of the rst reference time prole. Indeed the uctuations of the tail of a time prole can sometimes be wrongly identied as a second signal.

In the same gure (gure 5) the distribution obtained for the neutrino selection of the gadolinium analysis of ref. [12] is shown (black dots) for a direct comparison between the two. A clear excess of events at large t is present in the neutrino sample with respect to the cobalt one. This shows indeed the Double Chooz capability to observe o-Ps formation.

8

JHEP10(2014)032

10

10

10

60 Co source

Target candidates

n

JHEP10(2014)032

0 5 10 15 20 25 30

t (ns)

Figure 5. Distribution of the t value determined by the t for the cobalt sample (red squares), and for the neutrino sample (black dots), normalized to one.

To estimate the lifetime and the o-Ps production fraction, the obtained t distribution is tted with an exponential function with a lower bound of 5 ns, since this corresponds to the end of the smearing of the zero t events as observed in the 60Co sample study. A higher value for the t lower bound could be chosen granting a higher purity o-Ps sample, however the reduced statistics would lead to a larger error. Tests were made using values between 4 and 8 ns for the lower boundaries and the results were used to compute systematics related errors. Moreover, the algorithm outcome has a dependence on the choice of the reference time prole and this has also to be evaluated as systematic uncertainties. Additional parameters having an impact on the nal results are the vertex distance cut applied when building the reference time prole for a given calibration source, and its position inside the detector. To study all these aspects the analysis of the neutrino candidates was carried out under di erent conditions, namely:

1. Using as reference time prole the one obtained with a 60Co and a 137Cs source.

2. Modifying the distance cut between the reconstructed vertex and the true source position between 10 and 30 cm when building the reference time prole.

3. Selecting di erent positions inside the detector of the same source to build the reference time prole.

4. Changing the lower bound of time interval of the t between 4 and 8 ns.

Whereas for items 2 to 4 it is reasonable to assume that the best estimate is the mean of the results and to take as systematic uncertainty the spread, i.e. the r.m.s., of the di erent measurements, for item number 1 the situation is di erent. The best estimate corresponds to the measurement obtained with the reference time prole which has the energy that best matches the one of the neutrino candidate events.

9

Fitted

D

Entries/ns

n

Exponential fit

0.19 ns (stat.)

Target candidates

oPs

= 3.59 ns

t

oPs fraction = 53%

7% (stat.)

10

10

10

JHEP10(2014)032

1

0 5 10 15 20 25 30

Fitted

D

t (ns)

Entries/ns

10

Target candidates

n Exponential fit

0.24 ns (stat.)

oPs

t

= 3.77 ns

oPs fraction = 34.3%

6% (stat.)

10

10

1

0 5 10 15 20 25 30

Fitted

D

t (ns)

Figure 6. Fit of the t distribution, for neutrino candidates with energy between 1.2 and 3 MeV, with an exponential function for t greater than 5 ns. The 60Co (top) and the 137Cs (bottom) reference time proles are used in the oPs tagging algorithm. The errors quoted on the gures are statistical only.

For the neutrino candidates sample, the mean visible energy in the selected region from 1.2 to 3 MeV, is 2.237 MeV, which is more similar to the high energy bound given by the 60Co (2.5 MeV) rather than the low energy one of 137Cs (0.66 MeV). However, if o-Ps is formed we have two separate signals of 1.215 MeV (mean positron prompt energy) and1.022 MeV (o-Ps decay) respectively, therefore the best reference time prole would be in between the two energy regimes.

The systematics due to the choice of the source for the reference time prole (Ref)

is therefore evaluated as half the di erence between the low energy (i.e. 137Cs) and high

energy (i.e. 60Co) regimes:

Ref = 1

2(VCo

VCs) (5.1)

where VX is the value (i.e. either the o-Ps lifetime or its formation fraction) measured for the reference time prole obtained with the source X.

10

Error o-Ps formation o-Ps lifetime type fraction error [%] error [ns] Source element 9 0.09 Cut on the vertex distance 1.25 0.019

Source position 5 0.055 Fit interval 7 0.14 Total systematics 12 0.17

Statistics 5 0.15

Table 2. Evaluation of the di erent statistical and systematic uncertainties.

The ts obtained using the 60Co and 137Cs reference time proles are shown in gure 6, together with the results and the statistical errors, whereas a summary of all the di erent contributions to the systematic error is given in table 2. Note that in the fraction formation evaluation, the presence of background events in the neutrino candidate sample is accounted for. The expected background, in the 1.2 to 3 MeV energy region, is 0.5 events per day, whereas the neutrino candidates rate in the same region, after the applied analysis cuts reported at the end of section 4, is 13.3 events per day. The correction due to the presence of background corresponds therefore to about 3.8 %.

Adding quadratically the four systematic contributions, and taking the mean between the o-Ps properties measured using 137Cs and 60Co, it can be stated that the formation

fraction and the lifetime in the Target scintillator were observed to be 44 % 12 % (sys.) 5 % (stat.) and 3.68 ns 0.17 ns (sys.) 0.15 ns (stat.) respectively, in good agreement with the expectations.

6 Possible physics impact

The Borexino experiment has already proven the possibility to use o-Ps observation on a statistical way for e+/e discrimination reducing the cosmogenic 11C background in the solar pep neutron observation [11]. The demonstrated capability of Double Chooz to observe o-Ps formation on pulse shape analysis of single events is an important step forward in this technique, that could be exploited for the neutrino event selection and background reduction in liquid scintillation detectors.

In experimental projects aiming at the electron antineutrino detection, the e+/e discrimination could be directly used to reduce the background. This is not the case of IBD based experiments. Since o-Ps is formed only in about half of the events, a selection cut based on o-Ps observation is not conceivable for antineutrino events selection in IBD based experiments, even assuming a perfect detection e ciency. In this case, however, the o-Ps detection could be used to select a clean antineutrino subsample to validate experimental results. For example, one of the largest background contribution in IBD experiments, comes from the cosmogenic 9Li events, for which the prompt signal is an electron resulting therefore in no o-Ps formation observation.

11

JHEP10(2014)032

In order to prot the most from the o-Ps signature, several improvements would be desirable, in particular the liquid scintillator should be optimized to have a fast scintillating time and the longest possible o-Ps lifetime. Another critical parameter is the detector light yield: a large light yield would provide an easier o-Ps observation over a broader energy range. In addition, further improvements could come from the detector electronics, since a fast (sub-nanosecond) readout electronics would allow for the observation of small pulse shape distortions.

7 Conclusions

o-Ps formation could be exploited in anti-neutrino detector for additional background rejection. Although its detection is quite challenging in large liquid scintillator detectors due to its short lifetime of about 3 ns, this process has already been used in neutrino physics for a e+/e discrimination [11] on a statistical basis.

Double Chooz has demonstrated for the rst time the possibility to tag such a process on pulse shape analysis of single events. Relying on a selection algorithm based on the pulse shape distortion of the event, o-Ps formation was observed and its lifetime and formation fraction were measured, nding a good agreement with the values obtained in laboratory measurements.

The obtained result is so far energy dependent and it can not be used directly for a background reduction in Double Chooz, however it is now possible to assign a probability of each event of being an o-Ps decay for dedicated studies on pure samples.

Considering that the Double Chooz detector was not conceived for such a measurement (both in the choice of electronics and scintillator), this result is also an excellent starting point for future projects aiming at the liquid scintillation technology for anti-neutrino detection.

Acknowledgments

We thank the French electricity company EDF; the European fund FEDER; the Rgion de Champagne Ardenne; the Dpartement des Ardennes; and the Communaut des Communes Ardennes Rives de Meuse. We acknowledge the support of the ANR NuToPs project (grant 2011-JS04-009-01), the CEA, CNRS/IN2P3, the computer center CCIN2P3, and the UnivEarthS LabEx program of Sorbonne Paris Cit (ANR-10-LABX-0023 and ANR-11-IDEX-0005-02) in France; the Ministry of Education, Culture, Sports, Science and Technology of Japan (MEXT) and the Japan Society for the Promotion of Science (JSPS); the Department of Energy and the National Science Foundation of the United States; the Ministerio de Ciencia e Innovacin (MICINN) of Spain; the Max Planck Gesellschaft, and the Deutsche Forschungsgemeinschaft DFG (SBH WI 2152), the Transregional Collaborative Research Center TR27, the excellence cluster Origin and Structure of the Universe, and the Maier-Leibnitz-Laboratorium Garching in Germany; the Russian Academy of Science, the Kurchatov Institute and RFBR (the Russian Foundation for Basic Research); the Brazilian Ministry of Science, Technology and Innovation (MCTI), the Financiadora

12

JHEP10(2014)032

de Estudos e Projetos (FINEP), the Conselho Nacional de Desenvolvimento Cientco e Tecnolgico (CNPq), the So Paulo Research Foundation (FAPESP), and the Brazilian Network for High Energy Physics (RENAFAE) in Brazil.

Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.

References

[1] Double CHOOZ collaboration, Y. Abe et al., Improved measurements of the neutrino mixing angle 13 with the Double CHOOZ detector, arXiv:1406.7763 [http://inspirehep.net/search?p=find+EPRINT+arXiv:1406.7763

Web End =INSPIRE ].

[2] Double CHOOZ collaboration, Y. Abe et al., First measurement of 13 from delayed neutron capture on hydrogen in the Double CHOOZ experiment, http://dx.doi.org/10.1016/j.physletb.2013.04.050

Web End =Phys. Lett. B 723 (2013) 66 [arXiv:1301.2948] [http://inspirehep.net/search?p=find+J+Phys.Lett.,B723,66

Web End =INSPIRE ].

[3] Double CHOOZ collaboration, Y. Abe et al., Background-independent measurement of 13 in Double CHOOZ, http://dx.doi.org/10.1016/j.physletb.2014.04.045

Web End =Phys. Lett. B 735 (2014) 51 [arXiv:1401.5981] [http://inspirehep.net/search?p=find+EPRINT+arXiv:1401.5981

Web End =INSPIRE ].

[4] Daya Bay collaboration, F.P. An et al., Spectral measurement of electron antineutrino oscillation amplitude and frequency at Daya Bay, http://dx.doi.org/10.1103/PhysRevLett.112.061801

Web End =Phys. Rev. Lett. 112 (2014) 061801 [arXiv:1310.6732] [http://inspirehep.net/search?p=find+J+Phys.Rev.Lett.,112,061801

Web End =INSPIRE ].

[5] RENO collaboration, J.K. Ahn et al., Observation of reactor electron antineutrino disappearance in the RENO experiment, http://dx.doi.org/10.1103/PhysRevLett.108.191802

Web End =Phys. Rev. Lett. 108 (2012) 191802 [arXiv:1204.0626] [http://inspirehep.net/search?p=find+J+Phys.Rev.Lett.,108,191802

Web End =INSPIRE ].

[6] G. Ranucci, A. Goretti and P. Lombardi, Pulse-shape discrimination of liquid scintillators, http://dx.doi.org/10.1016/S0168-9002(98)00456-2

Web End =Nucl. Instrum. Meth. A 412 (1998) 374 [http://inspirehep.net/search?p=find+J+Nucl.Instrum.Meth.,A412,374

Web End =INSPIRE ].

[7] D. Franco, G. Consolati and D. Trezzi, Positronium signature in organic liquid scintillators for neutrino experiments, http://dx.doi.org/10.1103/PhysRevC.83.015504

Web End =Phys. Rev. C 83 (2011) 015504 [arXiv:1011.5736] [http://inspirehep.net/search?p=find+J+Phys.Rev.,C83,015504

Web End =INSPIRE ].

[8] H.J. Ache, Positronium and muonium chemistry, http://dx.doi.org/10.1021/ba-1979-0175.fw001

Web End =Adv. Chem. 175 (1979) 1 , American Chemical Society, U.S.A. (1979).

[9] Y. Kino et al., Positron annihilation in liquid scintillator for electron antineutrino detection,J. Nucl. Radiochem. Sci. 1 (2000) 63.

[10] G. Consolati et al., Characterization of positronium properties in doped liquid scintillators, http://dx.doi.org/10.1103/PhysRevC.88.065502

Web End =Phys. Rev. C 88 (2013) 065502 [arXiv:1308.0493] [http://inspirehep.net/search?p=find+J+Phys.Rev.,C88,065502

Web End =INSPIRE ].

[11] Borexino collaboration, G. Bellini et al., First evidence of pep solar neutrinos by direct detection in Borexino, http://dx.doi.org/10.1103/PhysRevLett.108.051302

Web End =Phys. Rev. Lett. 108 (2012) 051302 [arXiv:1110.3230] [http://inspirehep.net/search?p=find+J+Phys.Rev.Lett.,108,051302

Web End =INSPIRE ].

[12] Double CHOOZ collaboration, Y. Abe et al., Reactor electron antineutrino disappearance in the Double CHOOZ experiment, http://dx.doi.org/10.1103/PhysRevD.86.052008

Web End =Phys. Rev. D 86 (2012) 052008 [arXiv:1207.6632] [http://inspirehep.net/search?p=find+J+Phys.Rev.,D86,052008

Web End =INSPIRE ].

[13] C. Aberle et al., Large scale Gd-beta-diketonate based organic liquid scintillator production for antineutrino detection, http://dx.doi.org/10.1088/1748-0221/7/06/P06008

Web End =2012 JINST 7 P06008 [arXiv:1112.5941] [http://inspirehep.net/search?p=find+J+JINST,7,P06008

Web End =INSPIRE ].

[14] C. Aberle, C. Buck, F.X. Hartmann and S. Schonert, Light yield and energy transfer in a new Gd-loaded liquid scintillator, http://dx.doi.org/10.1016/j.cplett.2011.09.067

Web End =Chem. Phys. Lett. 516 (2011) 257 [http://inspirehep.net/search?p=find+J+CHPLB,516,257

Web End =INSPIRE ].

13

JHEP10(2014)032

[15] Y. Abe et al., The waveform digitiser of the Double CHOOZ experiment: performance and quantisation e ects on photomultiplier tube signals, http://dx.doi.org/10.1088/1748-0221/8/08/P08015

Web End =2013 JINST 8 P08015 [arXiv:1307.4917] [http://inspirehep.net/search?p=find+J+JINST,8,P08015

Web End =INSPIRE ].

[16] Particle Data Group collaboration, J. Beringer et al., Review of particle physics (RPP), http://dx.doi.org/10.1103/PhysRevD.86.010001

Web End =Phys. Rev. D 86 (2012) 010001 [http://inspirehep.net/search?p=find+J+Phys.Rev.,D86,010001

Web End =INSPIRE ].

[17] R. Brun and F. Rademakers, ROOT: an object oriented data analysis framework, http://dx.doi.org/10.1016/S0168-9002(97)00048-X

Web End =Nucl. Instrum. Meth. A 389 (1997) 81 [http://inspirehep.net/search?p=find+J+Nucl.Instrum.Meth.,A389,81

Web End =INSPIRE ].

JHEP10(2014)032

14