Search for relativistic magnetic monopoles with five years of the ANTARES detector data

Journal of High Energy Physics, Jul 2017

A search for magnetic monopoles using five years of data recorded with the ANTARES neutrino telescope from January 2008 to December 2012 with a total live time of 1121 days is presented. The analysis is carried out in the range β > 0.6 of magnetic monopole velocities using a strategy based on run-by-run Monte Carlo simulations. No signal above the background expectation from atmospheric muons and atmospheric neutrinos is observed, and upper limits are set on the magnetic monopole flux ranging from 5.7 × 10−16 to 1.5 × 10−18 cm−2·s−1·sr−1.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

https://link.springer.com/content/pdf/10.1007%2FJHEP07%282017%29054.pdf

Search for relativistic magnetic monopoles with five years of the ANTARES detector data

Revised: May Search for relativistic magnetic monopoles with ve years of the ANTARES detector data A. Albert 0 1 2 3 4 7 8 9 10 11 12 13 M. Andre 0 1 2 3 4 6 8 9 10 11 12 13 M. Anghinol 0 1 2 3 4 8 9 10 11 12 13 G. Anton 0 1 2 3 4 8 9 10 11 12 13 M. Ardid 0 1 2 3 4 5 8 9 10 11 12 13 J.-J. Aubert 0 1 2 3 4 8 9 10 11 12 13 T. Avgitas 0 1 2 3 4 8 9 10 11 12 13 g B. Baret 0 1 2 3 4 8 9 10 11 12 13 g J. Barrios-Mart 0 1 2 3 4 8 9 10 11 12 13 h S. Basa 0 1 2 3 4 8 9 10 11 12 13 i V. Bertin 0 1 2 3 4 8 9 10 11 12 13 S. Biagi 0 1 2 3 4 8 9 10 11 12 13 R. Bormuth 0 1 2 3 4 8 9 10 11 12 13 l S. Bourret 0 1 2 3 4 8 9 10 11 12 13 g M.C. Bouwhuis 0 1 2 3 4 8 9 10 11 12 13 k R. Bruijn 0 1 2 3 4 8 9 10 11 12 13 m J. Brunner 0 1 2 3 4 8 9 10 11 12 13 J. Busto 0 1 2 3 4 8 9 10 11 12 13 0 A. Marinelli 1 A. Deschamps 2 J.A.B. Coelho 3 C. Pellegrino 4 P. Migliozzi 5 I. Salvadori 6 Technical University of Catalonia, Laboratory of Applied Bioacoustics 7 G. Riccobene 8 Rambla Exposicio , 08800 Vilanova i la Geltru, Barcelona , Spain 9 I. Kreykenbohm 10 J. Hofestadt 11 D.F.E. Samtleben 12 H. van Haren 13 A. Enzenhofer aGRPHE, Universite de Haute Alsace, Institut universitaire de technologie de Colmar, 34 rue du Grillenbreit BP 50568, 68008 Colmar, France - dFriedrich-Alexander-Universitat Erlangen-Nurnberg, Erlangen Centre for Astroparticle Physics, Erwin-Rommel-Str. 1, 91058 Erlangen, Germany eInstitut d'Investigacio per a la Gestio Integrada de les Zones Costaneres (IGIC), Universitat Politecnica de Valencia, C/ Paranimf 1, 46730 Gandia, Spain f Aix Marseille Univ, CNRS/IN2P3, CPPM, Marseille, France hIFIC, Instituto de F sica Corpuscular (CSIC, Universitat de Valencia) c/Catedratico Jose Beltran, 2 E-46980 Paterna, Valencia, Spain qGran Sasso Science Institute, Viale Francesco Crispi 7, 00167 L'Aquila, Italy Emil-Fischer Str. 31, 97074 Wurzburg, Germany yDipartimento di Fisica e Astronomia dell'Universita, Viale Berti Pichat 6/2, 40127 Bologna, Italy zLaboratoire de Physique Corpusculaire, Clermont Universite, Universite Blaise Pascal, CNRS/IN2P3, BP 10448, F-63000 Clermont-Ferrand, France aaINFN, Sezione di Catania, Viale Andrea Doria 6, 95125 Catania, Italy abLSIS, Aix Marseille Universite CNRS ENSAM LSIS UMR 7296 13397 Marseille, France acInstitut Universitaire de France, 75005 Paris, France adRoyal Netherlands Institute for Sea Research (NIOZ), Landsdiep 4, 1797 SZ 't Horntje (Texel), The Netherlands aeDr. Remeis-Sternwarte and ECAP, Universitat Erlangen-Nurnberg, Sternwartstr. 7, 96049 Bamberg, Germany af Moscow State University, Skobeltsyn Institute of Nuclear Physics, Leninskie gory, 119991 Moscow, Russia agMediterranean Institute of Oceanography (MIO), Aix-Marseille University, 13288, Marseille, Cedex 9, France ahDipartimento di Fisica ed Astronomia dell'Universita, Viale Andrea Doria 6, 95125 Catania, Italy aiDirection des Sciences de la Matiere, Institut de recherche sur les lois fondamentales de l'Univers, Service de Physique des Particules, CEA Saclay, 91191 Gif-sur-Yvette Cedex, France ajINFN, Sezione di Pisa, Largo B. Pontecorvo 3, 56127 Pisa, Italy akDipartimento di Fisica dell'Universita, Largo B. Pontecorvo 3, 56127 Pisa, Italy alINFN, Sezione di Napoli, Via Cintia 80126 Napoli, Italy amDipartimento di Fisica dell'Universita Federico II di Napoli, Via Cintia 80126, Napoli, Italy anUniversite de Strasbourg, CNRS, IPHC UMR 7178, F-67000 Strasbourg, France aoUniversity Mohammed V in Rabat, Faculty of Sciences, 4 av. Ibn Battouta, B.P. 1014, R.P. 10000 Rabat, Morocco apDepartamento de F sica Teorica y del Cosmos & C.A.F.P.E., Universidad de Granada, Av. Severo Ochoa s/n, 18071 Granada, Spain aqUniversite de Toulon CNRS LSIS UMR 7296, 83957 La Garde, France arUniversite du Sud Toulon-Var, CNRS-INSU/IRD UM 110, 83957, La Garde Cedex, France ANTARES neutrino telescope from January 2008 to December 2012 with a total live time of 1121 days is presented. The analysis is carried out in the range > 0:6 of magnetic monopole velocities using a strategy based on run-by-run Monte Carlo simulations. No signal above the background expectation from atmospheric muons and atmospheric neutrinos is observed, and upper limits are set on the magnetic monopole ux ranging from 5:7 10 16 to 1:5 10 18 cm 2 s 1 sr 1. ArXiv ePrint: 1703.00424 1 Introduction 2 The ANTARES telescope 3 Detection of magnetic monopoles 4 Monte Carlo simulation 4.1 4.2 Magnetic monopole simulation Background simulation 5 Trigger and reconstruction 6 Event selection 7 Optimization of cuts 8 Results and discussion 9 Conclusion 1 Introduction { 1 { The concept of a particle with a magnetic charge (the magnetic monopole, MM in the following) was introduced by P.A.M. Dirac in 1931 [1] to explain the quantization of the elementary electric charge, e. The Dirac basic relation between e and the magnetic charge g is eg c = n~ 2 ! g = k gD = k e 2 ; where gD is the unit Dirac charge, k is an integer and ' 1=137 is the ne structure constant. The existence of magnetic charges and currents would symmetrize the Maxwell's equations. However, the symmetry would not be perfect, as gD is numerically much larger than e. In 1974, G. 't Hooft [2] and A.M. Polyakov [3] showed that the electric charge is naturally quantized in Grand Uni cation Theories (GUTs). MMs appear at the phase transition corresponding to the spontaneous breaking of the uni ed group into subgroups, one of which is U(1), describing electromagnetism. While there is no indication of the mass of the Dirac's magnetic monopole, in the context of GUTs the MM mass M is related to the mass of the X-boson carrier of the uni ed interaction (mX 1015 GeV/c2), yielding M & mX = ' 1017 GeV/c2. MMs with masses 1012 GeV/c2 (called intermediate-mass MMs) are predicted by theories with an intermediate energy scale between the GUT and the electroweak scales and would appear in the early Universe at a considerably later time than the GUT epoch [4]. More recently, it has been proposed [5] that solutions yielding MMs could arise within the electroweak theory itself. This Cho-Maison, or electroweak MM, would be expected to have a mass of the order of several TeV. Guided mainly by Dirac's argument and their predicted existence from spontaneous symmetry breaking mechanisms, searches have been routinely made for MMs produced at accelerators, in cosmic rays, and bound in matter [6, 7]. Eq. (1.1) de nes most of the MM properties, as they are assumed as point-like particles, of magnetic charge equal g, with unknown mass and with unknown relic cosmic abundance. To date, there are no con rmed observations of exotic particles possessing magnetic charge. MMs at the electroweak scale with M < 10 TeV are very good candidates for searches at the CERN Large Hadron Collider (LHC). The ATLAS collaboration [8] searched for MMs as highly ionizing particles produced in proton-proton collisions, leading to new cross section upper limits for spin 1/2 and spin 0 particles. MoEDAL is a dedicated experiment searching for MMs produced in high-energy collisions at the LHC using stacks of nucleartrack detectors and a trapping detector. Recently, limits on MM production cross sections have been reported both for the 8 TeV and 13 TeV LHC runs [9, 10]. GUT MMs are very massive and composite objects, well beyond the reach of any existing or foreseen accelerator. They could have been produced in a phase transition in the early Universe [ 11 ], and appeared as topological defects, about one pole for each causal domain. This would lead to a present-day overabundance [12]: the reduction of the number of MMs in the Universe was one of the motivating factors for cosmological in ation in Guth's original work [13]. As the Universe expanded and cooled down, the energy of MMs decreased: they would have reached a speed = v=c 10 10 during the epoch of galaxy formation (v is the MM speed and c is the speed of light in vacuum). After the gravitationally-driven galaxy formation epoch, galactic magnetic elds developed through the dynamo mechanism. Then, MMs were re-accelerated by these magnetic elds, yielding an isotropic intergalactic ux of relatively high-energy MMs. A magnetic eld B acting over a length ` increases the MM kinetic energy by a quantity gB`. The nal speed depends on galactic magnetic eld strength, on the coherent length ` and on MM mass and magnetic charge. For the typical values in our Galaxy, i.e. B are relativistic up to M 3 10 6 G and ` 300 pc = 1021 cm, MMs of g = gD 1011 GeV/c2. Then, their velocity decreases to reach the value ' 10 3 for M & 1017 GeV/c2. In models in which the cosmic magnetic eld, instead of being uniformly distributed, is strongly correlated with the large scale structure of the universe, MMs are relativistic up to 2 1013 GeV/c2 for g = gD [14]. The above MM acceleration process drains energy from the galactic magnetic eld. An upper bound on the ux of MMs in the galaxy (called the Parker bound [15]) has been obtained by requiring the rate of this energy loss to be small compared to the time scale on which the galactic eld can be regenerated. With reasonable choices for the astrophysical { 2 { HJEP07(21)54 parameters [6], the Parker bound corresponds to M . ( 10 15 description of the techniques used for the search of these particles is in [7], and a complete list of the results in [6]. Several searches were carried out also using neutrino telescopes. The ANTARES neutrino telescope [16] was completed in 2008 and the collected data can be used to search for MMs with energies high enough to yield light emission. The results of the analysis published in [17] using a data set of 116 days live time, lead to upper limits on the ux in the range between 1:3 10 17 and 5:7 10 16 cm 2 s 1 sr 1 for MMs with > 0:6. The IceCube collaboration has set upper limits on the ux for relativistic MMs ranging from of 1121 days collected from 2008 to 2012, increasing by a factor of 10 the live time of the previous published result. This analysis is based on a new selection of cuts, yielding a better separation of the MM signal from the background of atmospheric muons and neutrinos. Further, it relies on a new simulation strategy that reproduces each data run individually, allowing for an accurate reproduction of the data taking conditions. The paper is organized as follows: a brief description of the ANTARES telescope and the MM expected signatures are given in sections 2 and 3, respectively. The simulation and reconstruction algorithms are described in sections 4 and 5. The MM-sensitive observables, the selection strategy and the upper limit calculation are discussed in sections 6 and 7. Finally, the results are presented and discussed in section 8. 2 The ANTARES telescope The ANTARES detector [16] is an undersea neutrino telescope anchored 2475 m below the surface of the Mediterranean Sea and 40 km o shore from Toulon (France). It consists of 12 detection lines with 25 storeys per line and 3 optical modules (OMs) with 10-inch photomultipliers (PMTs) per storey. The detection lines are 450 m long and spaced 60 75 m apart horizontally. The main channel for neutrino detection is via the muons produced from high-energy muon neutrinos interacting inside, or in the vicinity of the detector. These muons move at relativistic velocities and induce the emission of Cherenkov light along their paths, detected by the optical modules. PMT signals corresponding to a charge above a threshold of 0.3 photo-electrons are integrated with a time window of 40 ns, digitised and denoted as hits. The readout of OMs is performed in the storey's Local Control Module, which collects the data in packages of 104 ms. These packages are sent to an on-shore farm of computers for further data processing and ltering. Each detector storey has one local clock that is synchronized to the on-shore master clock [19]. Furthermore, at the computer { 3 { farm a system of triggers is applied on the data (see section 5), selecting signatures which may correspond to the passage of relativistic particles. 3 Detection of magnetic monopoles The signature of a MM in a neutrino telescope like ANTARES is similar to that of a highly energetic muon. Thus, as in the case of electrically-charged particles, magneticallycharged particles induce the polarization of the dielectric medium. Coherent light emission (Cherenkov e ect) is induced by the restoring medium if the particle travels with a speed above the Cherenkov threshold th = 1=n, where n is the refractive index of the medium [20]. In water the threshold is th 0:74. The number of photons emitted from a MM with magnetic charge g in a small interval of path length, dx, and in the range d of wavelength, for th can be expressed as d2n d dx = 2 2 ng 2 e 1 1 where n is the number of photons emitted and is their wavelength; the remaining quantities are already de ned in eq. (1.1). For a given velocity, the Cherenkov radiation yield by a MM is a factor Ze ng 2 larger than that from a particle with electric charge Ze. Thus, for the refractive index of sea water, fast MMs with g = gD are expected to emit about 8550 times more Cherenkov photons than relativistic muons. In addition to the direct Cherenkov radiation, MMs can knock o atomic electrons ( -ray electrons) that can have velocities above the Cherenkov threshold, contributing to the total light yield. The production of -electrons is described by the di erential crosssection of Kasama, Yang and Goldhaber (KYG) [21] or by the more conservative (in terms of photon yield) Mott cross section [22]. The contributions to the light yield from these mechanisms are shown in gure 1. In both cases, some commonly accepted assumptions for the quantum-mechanical aspects of the interaction between a MM and an electron are used that must be implemented in the simulations. In this work, the Mott cross section is used, starting for the minimum velocity of = 0:5945: this allows a simpler application in the Monte Carlo simulation of the spectrum of the produced -ray electrons, yielding a safer estimate of the light yield. Contributions from radio-luminescence of water, pair production, Bremsstrahlung and photo-nuclear reactions induced by relativistic MMs are negligible compared to the direct and indirect Cherenkov light presented in gure 1, and are not taken into account in this analysis. In neutrino telescopes, the background of atmospheric muons dominates the solid angle region corresponding to down-going events. In particular, muons in bundle can easily be misidenti ed with the passage of a relativistic highly ionizing particle. On the opposite, the solid angle region corresponding to up-going events is almost background free, apart from the events induced by atmospheric neutrinos and the surviving down-going atmospheric muons misreconstructed as up-going. Due to the energy spectrum of atmospheric muon neutrinos, they usually induce minimum ionizing muons that can be easily distinguished from fast MMs. In order to suppress the irreducible background of atmospheric muons, only up-going MMs were considered. { 4 { that are directly produced per centimeter path length by a MM with g = gD, as a function of its velocity ( ). The number of photons produced by -rays with Mott cross section model [22] and KYG cross section model [21] and by a minimum ionizing muon are also shown. The request of up-going MMs reduces the range of masses M that can be observed in a neutrino telescope. The stopping-power de ned by S.P. Ahlen [23] has been used to estimate the absorption and energy loss of a MM when crossing the Earth. This work has established for MMs the equivalent of the Bethe-Bloch formula that describes the energy loss in the passage of a heavy electric charge by ionization and excitation in a non-conductive medium. Thus, the stopping-power of a MM crossing the Earth could be estimated using the simpli ed density pro le established by Derkaoui et al. [24]. Despite the high energy loss, MMs would remain relativistic and detectable as up-going events if M & 1010 GeV/c2 (see for instance gure 3 of [7]). As discussed in section 1, the MM speed depends on the characteristic of the galactic magnetic elds and on the mass M . Within reasonable astrophysical considerations, only MMs with a mass M . 1014 GeV/c2 can be expected in neutrino telescopes as an up-going event with a speed exceeding the Cherenkov threshold. Thus, the limits presented in this paper hold for MM in the mass range 1010 GeV/c2 . M . 1014 GeV/c2. 4 Monte Carlo simulation background events are discussed. 4.1 Magnetic monopole simulation In this section, the simulation of the MM signal and the atmospheric (neutrino and muon) Up-going MMs with one unit of Dirac charge, g = gD, have been simulated using nine equal width ranges of velocity in the region = [0:5945; 0:9950]. The nine intervals of the velocity are de ned in the rst column of table 1. { 5 { HJEP07(21)54 MMs have been simulated using a Monte Carlo program based on GEANT3 [25]. The simulation is independent of the MM mass and the incoming direction of MMs was distributed isotropically over the lower hemisphere. The propagation and detection of emitted photons is processed inside a virtual cylindrical surface surrounding the instrumented volume around the detector. A radius of 480 m is chosen to take into account the large amount of light emitted by MMs. The main source of background comes from up-going muons induced by atmospheric neutrinos and down-going atmospheric muons wrongly reconstructed as up-going tracks. The simulation of atmospheric muons is carried out using the generator MUPAGE [26] based on the parametrisation of the angle and energy distributions of muons under-water as a function of the muon bundle multiplicity [27]. MUPAGE produces muon events on the surface of the virtual cylinder. Up-going atmospheric neutrinos from the decay of pions and kaons are simulated using the package GENHEN [28, 29] assuming the model from the Bartol group [30, 31] which does not include the decay of charmed particles. The analysis presented in this paper is based on a run-by-run Monte Carlo simulation [32], which takes into consideration the real data taking conditions of the detector (e.g. sea water conditions, bioluminescence variability, detector status). 5 Trigger and reconstruction The applied triggers are based on local coincidences de ned as the occurrence of either two hits on two separate optical modules of a single storey within 20 ns, or one single hit of large amplitude, typically more than 3 photo-electrons. The trigger used for this analysis is de ned as a combination of two local coincidences in adjacent or next-to-adjacent storeys within 100 ns or 200 ns, respectively. In this analysis, only events passing such a trigger, well suited for MMs, are considered. The event reconstruction has been done with a slightly modi ed version of the algorithm described in [33]. By default, it assumes that particles travel at the speed of light. In order to improve the sensitivity for MMs travelling with lower velocities, the algorithm was modi ed such as to leave the reconstructed velocity of the particle t as a free parameter to be derived by the track t. The algorithm performs two independent ts: a track t and a bright-point t. The former reconstructs particles crossing the detector, while the latter reconstructs showering events, as those induced by e charged current interactions. Both ts minimize the same 2 quality function, thus, two parameters de ning the quality of these reconstructions are introduced, t 2 for the track t, and b 2 for the bright-point t. Some basic quality cuts have been applied to the data to ensure good data taking conditions [34]. To avoid any experimental bias, the search strategy is based on a blind analysis. The selection cuts applied on the analysis are established on Monte Carlo simula{ 6 { with error bars). For comparison, the distributions of the reconstructed t for MMs simulated in the velocity ranges [0:7280; 0:7725] (magenta histogram) and [0:7725; 0:8170] (green histogram) are also shown. All distributions correspond to events reconstructed as up-going. tions and using a test data sample of about 10% of the total data set, equivalent to 109 days out of the total 1121 days of live time. These runs are not used later for setting the limits. In the following comparisons between the test data sample and simulation, the full collection of Monte Carlo runs is used, and the 10% of test data is scaled to the total live time. Figure 2 shows the distribution of the reconstructed velocity t for MM events, atmospheric muons and neutrinos and compared to the test data sample. The neutrino distribution represents electron neutrinos and muon neutrinos for both neutral and charged currents. 6 Event selection In order to remove the bulk of down-going events, only up-going events with reconstructed zenith angles 90 are selected ( gure 3). Thus, the comparison shows a good agreement between the test data sample and simulation. The systematic uncertainties a ecting the predictions of atmospheric neutrino and atmospheric muon uxes are discussed in section 8. Accordingly, the event distributions of these two channels shown in this paper su er from an overall normalization uncertainty of about 30% and 35%, respectively. Additional cuts on the track t quality parameter are implemented to remove misreconstructed atmospheric muon tracks. In particular, the requirement t 2 b 2 is applied to favour events reconstructed as a track rather than those reconstructed as a bright point. The further event selections were optimized for di erent MM velocities. A di erent event selection was performed for each of the nine bins of reported in the rst column of table 1. The modi ed reconstruction algorithm which treats t as a free parameter was used only in the regions of low velocities between = 0:5945 and = 0:8170 ( ve bins). Thus, MMs with these velocities could be distinguished from particles traveling with the speed of light ( t = 1). For each of the ve low beta bins, only events reconstructed with t { 7 { band of 35% ( lled in gray), atmospheric neutrinos (blue histogram) and data (points with error bars). For comparison, the distributions of the reconstructed zenith angle for MMs simulated in the velocity ranges [0:7280; 0:7725] (magenta histogram) and [0:7725; 0:8170] (green histograms) are also shown. The peak at zenith = 0 represents wrongly reconstructed events. in the range of simulated were used in the nal selection. For example, at the range = [0:5945; 0:6390], only events with reconstructed velocity t = [0:5945; 0:6390] were selected. In the high velocity interval ranging from = 0:8170 to = 0:9950 (four bins), the t is not a discriminant variable anymore. However, MMs emit a large amount of light compared to that emitted from other particles, which allows them to be distinguished. In the used reconstruction algorithm, the hits from the optical modules belonging to the same storey are summed together to form a track hit. The coordinates of its position are coincident with the center of the storey, the time is equal to the time of the rst hit and the charge equal to the sum of the hits charges. For all velocity bins, the number of storeys with selected track hits Nhit, is used as a powerful discriminant variable since it refers to the amount of light emitted in the event (see gure 4). A second discriminative variable is introduced to further reduce the background, in particular for the velocities below the threshold for direct Cherenkov radiation where the light emission is lower. This variable, named , is de ned from a combination of the track t quality parameter t 2 and Nhit following [33]: = where Ndf is the number of free parameters in the reconstruction algorithm. It is equal to 6 when Example of t is included in the reconstruction, and 5 when the velocity is not reconstructed. distribution is shown at gure 5. This parameter has the advantage of including the track t quality parameter balanced with the brightness of the events, avoiding that bright events get cut by the condition applied on the t 2 variable. { 8 { 0 50 100 150 200 250 35% ( lled in gray), atmospheric neutrinos (blue histogram) and data (points with error bars). For comparison, the distributions of Nhit for MMs simulated in the velocity ranges [0:8170; 0:8615] (magenta histogram) and [0:9505; 0:9950] (green histogram) are also shown. At high velocities, Nhit provides a good discrimination for MM signals after applying the cuts zenith 105 α 1 . /0104 10­20 variable for atmospheric muons (red histogram) with an uncertainty band of 35% ( lled in gray), atmospheric neutrinos (blue histogram) and data (points with error bars). For comparison, the distribution of the variable for MMs simulated in the velocity range [0:7725; 0:8170] (magenta histogram) is also shown. Only events with reconstructed velocity t = [0:7725; 0:8170] were selected, and the cuts zenith b 2 have been applied. 7 Optimization of cuts The following step to suppress the atmospheric background is to use speci c cuts on the Nhit and parameters in order to maximize the signal-to-noise ratio. In gure 6, the event distribution of as a function of Nhit is shown for one range of MM velocity. This distribution indicates that a good separation of MM signal from background is achievable. { 9 { and Nhit, for atmospheric muons, atmospheric neutrinos, and MMs simulated in the velocity range [0:7280; 0:7725]. The cuts zenith t = [0:7280; 0:7725]. The vertical and horizontal lines show the cuts applied after optimization. No neutrinos survived at this range of . The horizontal and vertical lines show the e ect of the cuts. The signal region corresponds to the left upper quadrant. The 90% con dence level interval 90(nb; nobs), where nb is the number of background events is the 90% con dence interval de ned by the Feldman-Cousins approach [35]. It depends on the number of observed events nobs which is not known at this point because of the blind approach. Instead, the average con dence interval 90(nb) is calculated, from which the sensitivity of the analysis can be derived, by assuming a Poissonian probability distribution for the number of observed events nobs. The selection cuts are optimized by minimizing the so-called Model Rejection Factor (MRF) [36]: MRF = 90(nb) ; nMM (7.1) where nMM is the number of signal events remaining after the cuts, assuming an isotropic MM ux with 0MM = 1:7 10 13 cm 2 s 1 sr 1. In addition to the speci c values of the cuts, nMM depends on the detector acceptance Se (cm2 sr) and on the time period over which data was collected T (s). In order to compensate for the lack of statistics in the remaining sample of atmospheric muon background, an extrapolation has been performed in the region of interest for the signal. An example of extrapolation performed is shown in gure 7. After tting the Nhit distribution for muons with a Landau type function (red), the latter is extrapolated to the region of interest (pink), then the number of muons remaining after the nal cut on Nhit is given by the sum of the events from the muon histogram (blue) and the extrapolation (pink). Columns 3 and 4 of table 1 shows the background expectation, dominated by f o vee105 rbe104 m u N103 102 10 fit extrapolation N_!{hi Nhit cut Events remaining after this cut 0 50 100 150 200 250 300 Nhit The contribution of the extrapolation in the total number of events was taken into account in the optimization and the extrapolation uncertainties were computed. For this bin = [0:8170; 0:8615], 1.4 events are found after the cut Nhit > 91. atmospheric muons, for each bin of . After the optimization procedure and the estimation of the background, the 90% con dence level upper limit on the MM ux is obtained from the values of the cuts yielding the minimum value of the Model Rejection Factor MRF: 90% = 0 MM MRF: (7.2) 8 Results and discussion The unblinding was performed on the total set of data collected by the ANTARES telescope during ve years, which corresponds to 1012 active days live time after subtracting the 10% burn sample. No signi cant excess of data events is observed over the expected background, and the upper limits on ux have been found using eq. (7.2). Table 1 summarizes, for each of the nine bins of , the selection cuts, the number of expected background and observed events, and the 90% C.L. upper limits on the MM ux. The computation of the 90% C.L upper limits through eq. (7.2) includes the statistical uncertainties on the expected atmospheric muon events in column 3 of table 1. These uncertainties are dominant over the uncertainties related to the detector response. The e ects on the muon and neutrino rates due to the detector uncertainties are widely discussed elsewhere, particularly in [34, 37{39]. For the atmospheric neutrinos, the systematic uncertainties as a function of the energy are detailed in [38]. As shown in table 1, the contribution of atmospheric neutrinos is almost negligible with respect to atmospheric muons and the e ects of these uncertainties have been ignored. Concerning atmospheric muons, the dominant detector e ects are connected to the angular acceptance of the optical module [40] and to the absorption and scattering lengths in water [41]. The maximum 15% uncertainty on the optical module acceptance and the 10% on the light absorption length in water over the whole wavelength spectrum yields an overall +3350%% e ect on the expected range Selection cuts Number of Number of Number of atm. muons atm. neutrinos obs. events 90% C.L. (cm 2 s 1 sr 1) 5:9 10 16 of data taking). The selection cuts, the number of expected (muons and neutrinos) background and observed events and the upper limits on the ux are presented for each range of velocity ( ). The table was divided into two parts to distinguish the rst ve bins where t was assumed as a free parameter from the four bins where muon rate [37]. However, as already stated, in this case the dominant e ect (in most cases, with e ects larger than 50% on the number of events) is due to the lack in the statistics of the surviving muons and to the procedure for the background extrapolation, as described in gure 7. The values reported in column 3 represent the overall uncertainties on the surviving muon background in each bin. The e ect of a third uncertainty, due to the use of the Mott cross-section instead of the KYG (as discussed in section 3) has not been considered. In this case, a more conservative choice in terms of photon yield has been made. The outcome is to neglect a possible larger photon yield, that has the e ect of decreasing the detection thresholds towards smaller values of in gure 1. In the rst ve bins, the reconstructed velocity t was restricted to be compatible with the range of the MM velocity. Therefore, the event samples in these ranges are exclusive and must be added. As shown in table 1, the sum of background events in the rst ve ranges adds up to 5.4 events whereas only one event has been observed. This indicates a rather conservative method of extrapolating the atmospheric muon sample into the region de ned by the nal cuts. For the last four bins, t = 1 and cuts on and Nhit are tightened from bin to bin, that means bin 7 is a subset of bin 6 and so on. Thus, the total background is given here by bin 6 already. In gure 8 the ANTARES upper limits as a function of are presented, together with other experimental results from IceCube [18], MACRO [42] and Baikal [43], as well as the previous result from ANTARES [17] and the theoretical Parker bound [15]. The MACRO experiment was sensitive also to down-going candidates, surviving the 3000 meters of water equivalent of the Gran Sasso mountain overburden. Thus, their limit holds for MMs of lower mass (starting from 106 GeV/c2). For MMs that have to cross the Earth, as in ANTARES 2008 IceCube 86 IceCube 40 MACRO Baikal Parker bound HJEP07(21)54 r s . 10­18 10­109.55 0.6 active days live time (solid red line), compared to the upper limits obtained by other experiments [18, 42, 43], as well as the previous analysis of ANTARES (dashed red line) [17] and the theoretical Parker bound [15]. In [18] a more optimistic model for -rays production of MMs is used, making a direct comparison di cult. the case of the present paper, the limit is valid for M > 1010 GeV/c2. After applying the nal cuts to the unblinded data, two events have been observed. There is one event with Nhit = 93, = 0.5 and zenith = 27.4 which passes the cuts optimized of two bins of . It is identi ed as a bright well-reconstructed neutrino event regarding its physical properties, compatible with the total background observed at this range of high velocities. The second event with 0:728 is consistent with a down-going (zenith = 108.1 ) atmospheric muon yielding a bright shower. 9 Conclusion A search for relativistic MMs with the ANTARES neutrino telescope has been performed, using data collected during ve years (from 2008 to 2012) and corresponding to a total live time of 1012 days. No signal has been observed above the atmospheric background expectation and new upper limits on the MM ux have been set. Above the threshold for direct Cherenkov radiation 0:74, the limits found are better than those of other neutrino experiments. Below Cherenkov threshold, direct comparison is not straightforward due to the model of cross section used. Neutrino telescopes are well suited for the search for MMs. The future detector KM3NeT [44] will improve the sensitivity to the detection of MMs due to its large volume and high detection performance. Acknowledgments The authors acknowledge the nancial support of the funding agencies: Centre National de la Recherche Scienti que (CNRS), Commissariat a l'energie atomique et aux energies alternatives (CEA), Commission Europeenne (FEDER fund and Marie Curie Program), Institut Universitaire de France (IUF), IdEx program and UnivEarthS Labex program at Sorbonne Paris Cite (ANR-10-LABX-0023 and ANR-11-IDEX-0005-02), Labex OCEVU (ANR-11LABX-0060) and the A*MIDEX project (ANR-11-IDEX-0001-02), Region ^Ile-de-France (DIM-ACAV), Region Alsace (contrat CPER), Region Provence-Alpes-C^ote d'Azur, Departement du Var and Ville de La Seyne-sur-Mer, France; Bundesministerium fur Bildung und Forschung (BMBF), Germany; Istituto Nazionale di Fisica Nucleare (INFN), Italy; Stichting voor Fundamenteel Onderzoek der Materie (FOM), Nederlandse organisatie voor Wetenschappelijk Onderzoek (NWO), the Netherlands; Council of the President of the Russian Federation for young scientists and leading scienti c schools supporting grants, Russia; National Authority for Scienti c Research (ANCS), Romania; Ministerio de Econom a y Competitividad (MINECO): Plan Estatal de Investigacion (refs. FPA2015-65150-C3-1-P, -2-P and -3-P, (MINECO/FEDER)), Severo Ochoa Centre of Excellence and MultiDark Consolider (MINECO), and Prometeo and Grisol a programs (Generalitat Valenciana), Spain; Ministry of Higher Education, Scienti c Research and Professional Training, Morocco. We also acknowledge the technical support of Ifremer, AIM and Foselev Marine for the sea operation and the CC-IN2P3 for the computing facilities. 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. 133 (1931) 60 [INSPIRE]. [1] P.A.M. Dirac, Quantized Singularities in the Electromagnetic Field, Proc. Roy. Soc. Lond. A [2] G. 't Hooft, Magnetic Monopoles in Uni ed Gauge Theories, Nucl. Phys. B 79 (1974) 276 [3] A.M. Polyakov, Particle Spectrum in the Quantum Field Theory, JETP Lett. 20 (1974) 194 360 [hep-th/9601028] [INSPIRE]. Chin. Phys. C 40 (2016) 100001 [INSPIRE]. Sci. 65 (2015) 279 [arXiv:1510.07125] [INSPIRE]. [4] G. Lazarides, C. Panagiotakopoulos and Q. Sha , Magnetic Monopoles From Superstring Models, Phys. Rev. Lett. 58 (1987) 1707 [INSPIRE]. [5] Y.M. Cho and D. Maison, Monopoles in Weinberg-Salam model, Phys. Lett. B 391 (1997) [6] Particle Data Group collaboration, C. Patrignani et al., Review of Particle Physics, [7] L. Patrizii and M. Spurio, Status of Searches for Magnetic Monopoles, Ann. Rev. Nucl. Part. [8] ATLAS collaboration, Search for magnetic monopoles and stable particles with high electric charges in 8 TeV pp collisions with the ATLAS detector, Phys. Rev. D 93 (2016) 052009 [arXiv:1509.08059] [INSPIRE]. [9] MoEDAL collaboration, B. Acharya et al., Search for magnetic monopoles with the MoEDAL prototype trapping detector in 8 TeV proton-proton collisions at the LHC, JHEP 08 (2016) 067 [arXiv:1604.06645] [INSPIRE]. [10] MoEDAL collaboration, B. Acharya et al., Search for Magnetic Monopoles with the MoEDAL Forward Trapping Detector in 13 TeV Proton-Proton Collisions at the LHC, Phys. Rev. Lett. 118 (2017) 061801 [arXiv:1611.06817] [INSPIRE]. HJEP07(21)54 [INSPIRE]. (1979) 1365 [INSPIRE]. [12] J. Preskill, Cosmological Production of Superheavy Magnetic Monopoles, Phys. Rev. Lett. 43 [13] A.H. Guth, The In ationary Universe: A Possible Solution to the Horizon and Flatness Problems, Phys. Rev. D 23 (1981) 347 [INSPIRE]. [14] D. Ryu, H. Kang and P.L. Biermann, Cosmic magnetic elds in large scale laments and sheets, Astron. Astrophys. 335 (1998) 19 [astro-ph/9803275] [INSPIRE]. [16] ANTARES collaboration, M. Ageron et al., ANTARES: the rst undersea neutrino telescope, Nucl. Instrum. Meth. A 656 (2011) 11 [arXiv:1104.1607] [INSPIRE]. [17] ANTARES collaboration, S. Adrian-Martinez et al., Search for Relativistic Magnetic Monopoles with the ANTARES Neutrino Telescope, Astropart. Phys. 35 (2012) 634 [arXiv:1110.2656] [INSPIRE]. [18] IceCube collaboration, M.G. Aartsen et al., Searches for Relativistic Magnetic Monopoles in IceCube, Eur. Phys. J. C 76 (2016) 133 [arXiv:1511.01350] [INSPIRE]. [19] ANTARES collaboration, J.A. Aguilar et al., The data acquisition system for the ANTARES Neutrino Telescope, Nucl. Instrum. Meth. A 570 (2007) 107 [astro-ph/0610029] [INSPIRE]. [20] D.R. Tompkins, Total energy loss and Cerenkov emission from monopoles, Phys. Rev. 138 (1965) B248. 2935 [INSPIRE]. [INSPIRE]. [21] Y. Kazama, C.N. Yang and A.S. Goldhaber, Scattering of a Dirac Particle with Charge Ze by a Fixed Magnetic Monopole, Phys. Rev. D 15 (1977) 2287 [INSPIRE]. [22] S.P. Ahlen, Monopole Track Characteristics in Plastic Detectors, Phys. Rev. D 14 (1976) [23] S.P. Ahlen, Stopping Power Formula for Magnetic Monopoles, Phys. Rev. D 17 (1978) 229 Phys. 9 (1998) 173 [INSPIRE]. CERN Program Library Long Writeup W5013 (1993). 915 [arXiv:0802.0562] [INSPIRE]. [24] J. Derkaoui et al., Energy losses of magnetic monopoles and of dyons in the earth, Astropart. [25] CERN Application Software Group, GEANT 3.21 Detector Description and Simulation Tool, [26] G. Carminati, A. Margiotta and M. Spurio, Atmospheric MUons from PArametric formulas: A fast GEnerator for neutrino telescopes (MUPAGE), Comput. Phys. Commun. 179 (2008) reconstruction and its application to the ANTARES neutrino telescope, Astropart. Phys. 34 (2011) 652 [arXiv:1105.4116] [INSPIRE]. neutrino sources close to the Galactic Centre using the ANTARES neutrino Telescope, [astro-ph/0209350] [INSPIRE]. [arXiv:1007.1777] [INSPIRE]. muons measured with the 5-line ANTARES detector, Astropart. Phys. 34 (2010) 179 energy spectrum from 100 GeV to 200 TeV with the ANTARES telescope, Eur. Phys. J. C 73 L14 [arXiv:1108.0292] [INSPIRE]. [11] T.W.B. Kibble , Topology of Cosmic Domains and Strings, J. Phys. A 9 ( 1976 ) 1387 [15] E.N. Parker , The Origin of Magnetic Fields, Astrophys. J 160 ( 1970 ) 383 . Amsterdam ( 2003 ), http://www.vlvnt.nl/proceedings.pdf. detectors , Nucl. Instrum. Meth. A 725 ( 2013 ) 98 [INSPIRE]. Phys. Rev . D 53 ( 1996 ) 1314 [ hep -ph/9509423] [INSPIRE]. [27] Y. Becherini , A. Margiotta , M. Sioli and M. Spurio , A parameterisation of single and multiple muons in the deep water or ice, Astropart . Phys. 25 ( 2006 ) 1 [ hep -ph/0507228] [INSPIRE]. [28] J. Brunner , ANTARES simulation tools , in proceedings of The VLVnT workshop , [29] ANTARES collaboration, A. Margiotta, Common simulation tools for large volume neutrino [30] V. Agrawal , T.K. Gaisser , P. Lipari and T. Stanev , Atmospheric neutrino ux above 1-GeV , [31] G.D. Barr , T.K. Gaisser , S. Robbins and T. Stanev , Uncertainties in Atmospheric Neutrino Fluxes, Phys. Rev. D 74 ( 2006 ) 094009 [ astro -ph/0611266] [INSPIRE].


This is a preview of a remote PDF: https://link.springer.com/content/pdf/10.1007%2FJHEP07%282017%29054.pdf

A. Albert, M. André, M. Anghinolfi, G. Anton, M. Ardid, J.-J. Aubert, T. Avgitas, B. Baret, J. Barrios-Martí, S. Basa, V. Bertin, S. Biagi, R. Bormuth, S. Bourret, M. C. Bouwhuis, R. Bruijn, J. Brunner, J. Busto, A. Capone, L. Caramete, J. Carr, S. Celli, T. Chiarusi, M. Circella, J. A. B. Coelho, A. Coleiro, R. Coniglione, H. Costantini, P. Coyle, A. Creusot, A. Deschamps, G. De Bonis, C. Distefano, I. Di Palma, A. Domi, C. Donzaud, D. Dornic, D. Drouhin, T. Eberl, I. El Bojaddaini, D. Elsässer, A. Enzenhöfer, I. Felis, L. A. Fusco, S. Galatà, P. Gay, V. Giordano, H. Glotin, T. Grégoire, R. Gracia Ruiz, K. Graf, S. Hallmann, H. van Haren, A. J. Heijboer, Y. Hello, J. J. Hernández-Rey, J. Hößl, J. Hofestädt, C. Hugon, G. Illuminati, C. W James, M. de Jong, M. Jongen, M. Kadler, O. Kalekin, U. Katz, D. Kießling, A. Kouchner, M. Kreter, I. Kreykenbohm, V. Kulikovskiy, C. Lachaud, R. Lahmann, D. Lefèvre, E. Leonora, M. Lotze, S. Loucatos, M. Marcelin, A. Margiotta, A. Marinelli, J. A. Martínez-Mora, R. Mele, K. Melis, T. Michael, P. Migliozzi, A. Moussa, S. Navas, E. Nezri, M. Organokov, G. E. Păvălas, C. Pellegrino, C. Perrina, P. Piattelli, V. Popa, T. Pradier, L. Quinn, C. Racca, G. Riccobene, A. Sánchez-Losa, M. Saldaña, I. Salvadori, D. F. E. Samtleben, M. Sanguineti, P. Sapienza, F. Schüssler, C. Sieger, M. Spurio, Th. Stolarczyk, M. Taiuti, Y. Tayalati, A. Trovato, D. Turpin, C. Tönnis, B. Vallage, V. Van Elewyck, F. Versari, D. Vivolo, A. Vizzoca, J. Wilms, J. D. Zornoza, J. Zúñiga, The ANTARES collaboration. Search for relativistic magnetic monopoles with five years of the ANTARES detector data, Journal of High Energy Physics, 2017, 54, DOI: 10.1007/JHEP07(2017)054