PROCEEDINGS OF THE 31st ICRC, ŁOD´ Z´ 2009 1 Search for GUT monopoles at Super-Kamiokande Koh Ueno∗ for the Super-Kamiokande collaboration ∗Kamioka Observatory, ICRR, University of Tokyo, Higashi-Mozumi, Kamioka, Hida, Gifu 506-1205, Japan Abstract. GUT monopoles captured by the Sun’s process above subsequently decay into positive pions, 0 + + gravitation are expected to catalyze proton decays (ρ ,ω,η,K , · · ·) → π , νe, νµ and ν¯µ are produced via the Callan-Rubakov process. In this scenario, pro- as described below, tons, which initially decay into pions, will ultimately + + π → µ + νµ (2) produce νe, νµ and ν¯µ. After going through neutrino + + oscillation, all neutrino species appear when they µ → e + νe +¯νµ (3) arrive at the Earth, and can be detected by a 50,000 ton water Cherenkov detector, Super-Kamiokande (SK). A search for neutrinos produced through such a process has been carried out with SK and is expected 0.5 to give the cross section and velocity dependent limit 0.45 −23 −2 −1 on the monopole flux, F < 1 × 10 cm s 0.4 sr−1 for β = 10−3 and σ = 1 mb at 90%C.L. This is one order of magnitude more stringent a limit 0.35 than the previous result, which was obtained by the 0.3 Kamiokande experiment. 0.25 Keywords: GUT monopole, neutrino, water 0.2 ν ν± e h Cherenkov detector Arbitrary Unit 0.15 ν h I. INTRODUCTION 0.1 Grand Unified Theories (GUT) predict superheavy 0.05 ν± e monopoles (GUT monopoles) produced in the very 0 early universe. GUT monopoles are predicted to have 0 10 20 30 40 50 60 70 appeared as topological defects at the phase transition of Eν (MeV) vacuum, where GUT gauge group spontaneously broke to leave the U(1) of electromagnetism. If we take as Fig. 1. Expected neutrino spectra at SK detector. Neutrino oscillation an example the temperature of the phase transition at effect is taken into account. µ, τ components are grouped together as h 15 since these two flavor neutrinos have the same cross section in electron which the monopoles were formed equal to 10 GeV scattering. and assume the average production rate of about one monopole per horizon at that time [1], the density of the monopoles today exceeds critical density of the universe After going through neutrino oscillation, all neutrino by more than 14 orders of magnitude. Even though the species appear when they arrive at the Earth (Fig. 1), inflationary universe scenario [2], [3] overcomes this and such low energy neutrino events can be detected by problem, the monopole flux in the universe depends on a water Cherenkov detector. some parameters such as monopole mass or the reheating II. DETECTOR AND EXPECTED SIGNAL temperature, and therefore the uncertainty remains. In fact, due to the wide variety of elementary particle We searched for monopole-induced neutrinos using models, several models coexist with the Parker limit (∼ Super-Kamiokande (SK) [10], a large water Cherenkov 10−15 cm−2 s−1 sr−1) [4], [5], [6], and a flux in that detector located at Kamioka mine in Japan. SK is a range can be relatively easily detected by underground high-performance neutrino detector consisting of 11,146 experiments. Arafune et al. [7] pointed out copious low 20-inch PMTs and 50,000 tons of pure water. The energy neutrinos might be emitted when monopoles fiducial volume of this search, which amounts to 22,500 accumulating inside the Sun catalyze proton decay; tons, is defined to be more than 2 meters from the walls of the inner detector as there are gamma ray p → (ρ0,ω,η,K+, · · ·)+ e+(or µ+) (1) backgrounds near the wall. Monopole-induced neutrinos along their paths with cross sections typical of strong includes all six types of neutrinos, and electron scatter- − − interactions via so called Callan-Rubakov process [8], ing (νx(¯νx)+e → νx(¯νx)+e ) in the fiducial volume [9]. But, contrary to spontaneous proton decay, direct is used for the monopole-induced neutrino search. Fig. 2 neutrino production process , such as p → ν¯ + π+, shows the recoil electron spectra in SK, where the values is forbidden [8]. When decay mesons produced by the of the reference [11] are used as neutrino oscillation 2 UENO et al. GUT MONOPOLE SEARCH AT SUPER-KAMIOKANDE parameters. Matter effects on oscillation in the Sun are this cut. The source of muons is from the interaction also taken into account. of atmospheric muon neutrinos. Electrons with E > 18 ◦ MeV have a Cherenkov angle of θC ∼ 42 . Muons that III. DATA REDUCTION still remain in the monopole-induced neutrino search In this analysis, 1496 days’ data of SK-I taken from have pµ less than ∼ 350MeV, which corresponds to θC ◦ ◦ April 1996 to July 2001 is used. The backgrounds to ∼ 38 ; thus, the events with θC < 38 were removed the monopole-induced neutrino search are mainly atmo- from the data. The Cherenkov angle cut is also used to spheric neutrinos and muon-induced spallation products. remove events with θC > 50 degrees, which eliminates To minimize the contribution of the latter backgrounds, events without a clear Cherenkov ring pattern, such as the lower energy threshold in this analysis is 18 MeV multiple γ rays emitted by a nuclear de-excitation. in total electron energy. The higher energy threshold is Some events originating from outside of the fidu- set to be 55 MeV, which is the end-point of the recoil cial volume have the possibility of being reconstructed electron spectrum (Fig. 2). within the fiducial volume of SK. γ-rays from the mate- rials of the detector structure and surrounding rock are the source of these backgrounds. To remove such events, 0.1 the reconstructed direction is projected backwards from 0.09 the vertex position until it reaches the inner detector ν ν± ν ν± wall, and events with this distance between the vertex 0.08 e+ e+ h+ h and the wall greater than 450 cm are removed from the 0.07 data. 0.06 Finally, charged current interactions of atmospheric νµ ν 0.05 e produce muons and decay electrons. Some low energy muons ( 160 MeV) or decay electrons survive in this 0.04 ∼ Arbitrary Unit analysis but are tagged by preceding muons. Some frac- 0.03 tion of those muons which passed the Cherenkov angle ν 0.02 h cut can be identified by the following decay electron. ν± 0.01 h In order to remove those backgrounds we eliminate ν± e 0 the events which have time-correlated events before or 0 10 20 30 40 50 60 70 after the candidate events. However, the decay electrons Ee (MeV) whose mother muons are invisible (below Cherenkov threshold), Michel electrons, cannot be dropped by this Fig. 2. Recoil electron spectrum in the SK described analytically cut; how to handle these events is described later. without energy resolution included. How to read this figure is the same of Fig. 1 IV. ANALYSIS AND RESULTS The energy spectrum after each selection criterion is shown in Fig. 3. After all cuts are applied, 163 Besides setting the fiducial volume and energy crite- monopole-induced neutrino candidates remain in the ria, some other background reductions are additionally energy range from 18 MeV to 55 MeV. Most these candi- applied to the data. After a “standard cuts” (total charge dates are due to three kinds of irreducible backgrounds. cut, outer detector trigger cut, flasher event cut, etc...), The first background is spallation events from cosmic the data are subjected to a normal spallation cut. Cosmic ray muons which still remain past the two spallation cuts ray muons can spall oxygen and induce unstable nuclei described above. called spallation products (µ + 16O → µ + X). It The other two backgrounds are ν and ν com- is one of the most abundant backgrounds in < ∼20 e µ ponents of atmospheric neutrinos. The “ν component” MeV region, and as described above, the ability to e means atmospheric ν and ν¯ , while the “ν component” remove the background determines the lower threshold e e µ is atmospheric ν and ν¯ which produce muons below of the monopole-induced neutrino search. The spallation µ µ Cherenkov threshold whose decay electrons are then background is reduced by a likelihood method that uses observed in the detector (invisible muon). The spectrum timing and track information of the muons preceding of these “invisible muons” is a Michel spectrum. the candidate events. The same algorithm is used in the To subtract these backgrounds, angular distribution standard solar neutrino analysis of SK. In addition to the with respect to the direction of the Sun was used. As normal spallation cut, a tighter criterion is then applied shown in Fig. 4, the data are divided into 20 angle bins in order to enhance the rejection efficiency of spallation and the following χ2 function is minimized with respect background. The spallation products with the shortest to α for each β. half lives, such as 11Li and 12N, tend to create spallation events with the highest energies. So, we remove events 20 2 {Nreal(i) − αNflat(i) − βNMC(i)} which occur within 0.15 sec from cosmic ray muons. χ2(α, β) ≡ Next, a Cherenkov angle cut is applied to the re- σ2 i=1 i maining data set. The remaining muons are removed by X (4) PROCEEDINGS OF THE 31st ICRC, ŁOD´ Z´ 2009 3 10 3 2.25 after 1st reduction after spallation cut 2 after ultimate spallation cut after goodness cut 1.75 after Cherenkov angle cut 10 2 after effective wall cut 1.5 after subevent cut 1.25 1 10 0.75 0.5 Events / bin 1496 days 22.5 kton 0.25 Events / MeV 1496 days 22.5 kton 1 0 20 25 30 35 40 45 50 55 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 θ Energy (MeV) cos Sun Fig.