Thèse présentée par. Détermination du flux de muons atmosphériques avec le télescope à neutrinos ANTARES - PDF Free Download\nDownload "Thèse présentée par. Détermination du flux de muons atmosphériques avec le télescope à neutrinos ANTARES"\n1 UNIVERSITE PARIS CEA/IRFU/SPP APC DIDEROT (Paris 7) Thèse présentée par Claire Picq Détermination du flux de muons atmosphériques avec le télescope à neutrinos ANTARES Pour l obtention du titre de Docteur de l Université Paris VII Spécialité : Champs, Particules, Matières Soutenue le 18 juin 2009, devant le jury composé de M. Etienne Parizot (Président) Mme. Annarita Margiotta (Rapporteur) M. Fabrice Feinstein (Rapporteur) M. Jorg R. Hörandel (Examinateur) M. Bertrand Vallage (Directeur) M. Antoine Kouchner (Directeur)\n3 3 A mes parents A Xavier\n5 Remerciements Tout d abord, je tiens à remercier mes deux directeurs de thèse, Bertrand et Antoine, pour avoir encadré cette thèse et m avoir fait confiance sur la façon de mener mes recherches. Je remercie Sotiris pour son aide aux moments critiques de ma thèse et ses bonnes idées qui m ont permis d avancer. Je remercie Philippe, mon "parrain" du CEA, pour son soutien lors mes coups de blues et ses précieux conseils dans l orientation de ma vie professionnelle. Je remercie Bruno Mansoulié et Ursula Bassler pour leur accueil au sein du Service de Physique des Particules au Commissariat à l Energie Atomique et Pierre Binétruy au sein du laboratoire AstroParticules et Cosmologie de l Université Denis Diderot. Je remercie les membres du jury : merci à Etienne Parizot pour avoir accepté d en prendre la présidence ; merci à mes rapporteurs, Annarita Margiotta qui a suivi mes travaux et m a apporté son aide tout au long de cette thèse en travaillant sur les simulations ANTARES des muons atmosphériques, et Fabrice Feinstein pour avoir envoyé ses commentaires ; enfin, merci à Jorg R. Hörandel de s être intéressé à mes travaux en participant à ce jury. Je remercie l équipe du CEA et tout particulièrement Niccolò avec qui j ai partagé mon bureau si agréablement, Gilles pour sa bonne humeur, Yvonne pour son aide au moment le plus crucial qu est le début d une thèse, ainsi que toutes les personnes de l équipe pour leur grande connaissance du détecteur. Je remercie l équipe de l APC : Corinne pour les études que nous avons faites sur l électronique, Bruny pour son soutien, Véronique pour sa sympathie. Je remercie la collaboration ANTARES pour son accueil et son écoute. Je remercie le Laboratoire d Informatique de Paris 6 et tout particulièrement Nicolas, Lom, Alexandre, Alban et Fabrice pour leur aide en C++ et pour avoir accueilli ma thèse sur leur serveur. Enfin je remercie énormément ma famille et tout particulièrement mes parents qui m ont soutenue tout au long de mes études et Xavier pour sa présence à mes côtés pendant cette thèse.\n7 RESUME Le télescope à neutrinos ANTARES est un détecteur sous-marin situé en mer Méditerranée. L univers étant transparent aux neutrinos, leur étude est un moyen unique d améliorer nos connaissances sur la nature et l origine des rayons cosmiques, ainsi que sur les mécanismes d émission des sources astrophysiques les plus puissantes de l univers. Les neutrinos offrent en outre la possibilité d ouvrir une nouvelle fenêtre d observation en énergie (> TeV) sur l univers. Cette thèse est consacrée à l étude du bruit de fond principal du détecteur, dû au passage de muons atmosphériques produits lors de la désintégration des rayons cosmiques de haute énergie avec les noyaux de l atmosphère. La première partie de cette thèse est dédiée à l étude du détecteur. Les différentes caractéristiques du détecteur, son étalonnage ainsi que les techniques de réglages de l électronique y sont détaillés. La seconde partie de cette thèse fait état des différents résultats obtenus sur les muons atmosphériques avec un détecteur à cinq lignes. Une présentation détaillée de la simulation utilisée y est présentée. La première difficulté de détection des muons atmosphériques est due à la géométrie du détecteur, la seconde au fait que les muons atmosphériques arrivent souvent en gerbes et que le nombre de muons dans une gerbe à 2500 m de profondeur est encore inconnu. Une première étude basée sur des simulations permet de discriminer entre les muons seuls et les gerbes de muons. Une seconde étude est dédiée à la mesure du flux de muons en fonction de la profondeur de détection. En prenant en compte les erreurs systématiques, la mesure est compatible avec les résultats d autres instruments. Mots-clés : télescope, ANTARES, rayons cosmiques, neutrinos, étalonnage, muons atmosphériques, flux de muons ABSTRACT The neutrino telescope ANTARES is a deep-sea detector located in the Mediterranean Sea. The universe is transparent to neutrinos, so their study provides a unique means of improving our knowledge of the nature of cosmic rays, their origins and their emission from the most powerful astrophysical sources in the cosmos. Neutrinos also offer the possibility of opening a new energy window (>TeV) for observation of the universe. This thesis is dedicated to the study of the main background noise of the detector, due to the passage of atmospheric muons produced by high energy cosmic rays interacting with atmospheric nuclei. The first part of this thesis focuses on the study of the detector. The different characteristics and the calibration of the detector as well as the techniques of monitoring the electronic are described. The second part of this thesis reports the various results obtained on the atmospheric muons with the five line detector. A detailed presentation of the simulations used is presented. The first difficulty of detecting atmospheric muons is due to the geometry of the detector. The second is due to the fact that the atmospheric muons often arrive in bundles and that the number of muons in these bundles is unknown at a depth of 2500 m. A first study based on simulations makes it possible to discriminate between the muons alone and the bundles of muons. A second study is dedicated to the measurement of the muon flux depending on the slant depth. The measurement is compatible with the results of other instruments when the systematic uncertainties are taken into account. Key-words : telescope, ANTARES, cosmic ray, neutrinos, calibration, atmosperic muons, muon flux\n9 Table des matières Remerciements 5 Introduction 13 I Rayons cosmiques 15 1 Les rayons cosmiques Historique Le spectre en énergie Composition Accélération des rayons cosmiques Origine des rayons cosmiques Propagation Et les neutrinos? Sources ponctuelles éventuellement productrices de neutrinos a Les supernovae et leurs restes b Les microquasars c Les noyaux actifs de galaxie d Les sursauts gamma e La matière noire f Etude des sources à l aide des gamma : l expérience H.E.S.S Méthodes de détection des rayons cosmiques L expérience BESS Les détecteurs au sol a L expérience KASCADE-grande b L expérience Auger Un détecteur souterrain : MACRO II Télescopes à neutrinos 37 1 Principe de détection Interaction du neutrino avec la matière Propagation dans la Terre Propagation des muons induits a Corrélation neutrino-muon b Interactions dominantes c Effet Tcherenkov Bruits de fond physiques Le détecteur ANTARES\n10 TABLE DES MATIÈRES 2.1 Description du détecteur a Les lignes b Modules optiques Caractéristiques optiques du milieu a Longueur d atténuation dans l eau b Le bruit de fond optique L électronique de lecture des photo-multiplicateurs Le système d acquisition des données Algorithmes de filtrage Reconstruction des traces de muons a Principe de reconstruction b Reconstruction des gerbes c Description des reconstructions Les autres télescopes à neutrinos Un télescope en eau douce : Baïkal Les télescopes sous la glace a AMANDA b IceCube Les télescopes en eau salée a Dumand b Nemo c Nestor d KM3Net III Etalonnage des lignes du détecteur Antares 71 1 Méthodes d étalonnage Etalonnage en temps Etalonnage du TVC Décalage temporel des étages a Horloge à Terre b Détermination des T 0 à Terre c Contrôle in situ Effet de Walk Non-linéarité différentielle des ADCs Le bit flip Etalonnage de la charge Le piédestal Le photo-électron Baisse du gain des modules optiques La diaphonie Détermination des seuils a Mesure des seuils b Nombre d événements par OM c La méthode des T : nouvel étalonnage des seuils Principe de la mesure 91\n11 TABLE DES MATIÈRES 11 IV Simulation des muons atmosphériques 93 1 Flux des rayons cosmiques interagissant dans la haute atmosphère Le modèle poly-gonato Le modèle NSU La paramétrisation d "HEMAS" Comparaison des différents modèles Les modèles d interactions hadroniques à haute énergie Propagation des muons jusqu au détecteur Simulation du détecteur Longueur d atténuation dans l eau Acceptance des modules optiques V Etude de la multiplicité La multiplicité dans la simulation Comparaison de la multiplicité à différents niveaux de la simulation Les muons et l algorithme de filtrage 3N a Composition et énergie des rayons cosmiques b L angle zénithal Comparaison de quelques variables entre les données et les simulations Qualité de la reconstruction : Λ Fraction de petits résidus Charge totale Longueur de la trace Tenseur d inertie Angle zénithal Etude classique par coupures Analyse multivariable Décorrélation des variables Variables utilisées Résultats VI Détermination du flux de muons Sélection des événements Nombre de lignes Evènements à une ligne a χ 2 de la trace reconstruite b χ 2 de la gerbe c Nombre de coups d Angle zénithal Evènements à plusieurs lignes a χ 2 de la trace reconstruite b χ 2 de la gerbe c Angles zénithal et azimuthal Efficacité et résolution angulaire Détermination du flux Comparaison entre les données et les simulations Multiplicité au niveau de la canette Calculs d erreurs statistiques sur les efficacités\n12 12 TABLE DES MATIÈRES 2.3.a Evènements non pondérés (MUPAGE) b Evènements pondérés (CORSIKA) Efficacité de déclenchement Surface effective des muons au niveau du déclencheur Efficacité de la reconstruction Erreurs statistiques sur le flux de muons Correction au flux de muons a Perte par temps mort du filtrage des évènements sur les données 5 lignes b Correction du centre de cellule : c Correction des seuils : Effets systématiques a Acceptance angulaire b Longueur d absorption c Etalonnage en charge, Réglage des seuils d Bio-salissure des photo-multiplicateurs e Simulation de la multiplicité f Reconstruction Flux de muons Vérification de la méthode Relation profondeur - intensité du flux Verticalisation Résultats Conclusion 167 A Différentes méthodes de discrimination utilisées dans TMVA Discriminant de Fisher Rapport de vraisemblance Réseau de neurones Les arbres de décision\n13 Introduction Depuis leur découverte au début du vingtième siècle, les rayons cosmiques n ont cessé d intriguer les physiciens. Leur origine et leur composition demeurent controversées. Le spectre en énergie des particules observées varie sur 32 ordres de grandeur et s étend sur une très grande gamme jusqu à des énergies macroscopiques qui frôlent les limites des modèles classiques d accélération. L étude des rayons cosmiques nécessite donc des approches variées. Aux plus hautes énergies, des détecteurs couvrant plusieurs milliers de km 2 au sol sont nécessaires, tandis que l étude détaillée de la composition aux énergies intermédiaires requiert d embarquer des détecteurs de plus en plus sophistiqués dans des ballons ou des satellites. Si la source des rayons cosmiques reste inconnue, c est que les rayons cosmiques, pour l essentiel électriquement chargés, sont défléchis par les champs magnétiques lors de leur propagation. À cet égard, l étude des neutrinos de haute énergie offre de nouvelles perspectives. En effet, en interagissant à la source avec la matière ou le rayonnement ambiant, les rayons cosmiques peuvent produire des neutrinos s échappant de la zone de production et se propagent en ligne droite. Le développement de cette nouvelle "astronomie neutrino" représente un défi technologique, conforté par la détection au MeV de neutrinos issus de l explosion de la supernova 1987A, en février A plus haute énergie (jusqu à 9 GeV) les détecteurs ANTARES (Astronomy with a Neutrino Telescope and Abyss environmental RESearch, dans l hémisphère nord) et IceCube (dans l hémisphère sud) scrutent désormais le ciel ouvrant ainsi une nouvelle fenêtre d observation sur l Univers. La détection des neutrinos est indirecte. Les neutrinos cosmiques de haute énergie qui interagissent lors de leur traversée de la Terre produisent des muons qui émettent un cône de lumière ténue -dite lumière Tcherenkov- lors de leur passage dans le milieu marin (ou la glace). Cette lumière est détectée par un réseau tridimensionnel de modules optiques (photomultiplicateurs enfermés dans une sphère en verre), ce qui permet de remonter aux caractéristiques du muon et donc du neutrino incident. Les muons descendants provenant des gerbes atmosphériques produites par l interaction des rayons cosmiques constituent une source intense de bruit de fond, cependant ils servent aussi comme source d étalonnage. Cette thèse est consacrée à l étude des muons descendants avec le télescope à neutrinos ANTARES. Le premier chapitre de la thèse décrit brièvement les enjeux de la physique des rayons cosmiques. Le spectre en énergie des rayons cosmiques, ainsi qu une sélection d expériences ayant participé à l élaboration de ce spectre sont décrits. Une attention particulière sera portée aux objets astrophysiques susceptibles de produire des rayons cosmiques (et des neutrinos) de haute énergie. Le deuxième chapitre expose le principe de détection des neutrinos et des muons avec les télescopes à neutrinos. Y sont aussi présentés les différents projets de détecteurs de neutrinos et notamment ANTARES. Les composants et l agencement du détecteur seront exposés, ainsi que les algorithmes de filtrage et de reconstruction des données. Ces connaissances sont nécessaires à la compréhension des études effectuées dans le cadre de cette thèse. L étalonnage des lignes du détecteur ANTARES fait l objet du troisième chapitre. Nous insisterons d une part sur les méthodes d étalonnage développées lors de l intégration du détecteur à terre, et sur 13\n14 14 INTRODUCTION les méthodes utilisées in situ d autre part. Une étude de stabilité de certains de ces étalonnages sera présentée. La deuxième partie de la thèse est consacrée à l étude des muons. Le principe de la mesure du flux de muons est tout d abord abordé. La simulation des muons atmosphériques est alors exposée (chapitre 4), depuis la simulation des rayons cosmiques en haut de l atmosphère jusqu à la propagation de la lumière Tcherenkov produite par des muons traversant le détecteur, en passant par la propagation des muons de chaque gerbe jusqu au détecteur. Une étude de la multiplicité de ces muons, dans le détecteur ANTARES, est proposée dans le cinquième chapitre de la thèse. Après une comparaison des variables entre simulations et données, une étude classique par coupures puis une analyse multi-variable sont développées afin de séparer les muons isolés des gerbes contenant plus d un muon. Le dernier chapitre est consacré à l étude du flux de muons atmosphériques. La mesure de la dépendance zénithale de ce flux est établie à partir de plusieurs algorithmes de reconstruction ainsi que de plusieurs simulations Monte Carlo. Les différents résultats obtenus sont comparés et transcrits en termes de relation intensité-profondeur pour comparaison avec d autres expériences.\n15 Chapitre I Rayons cosmiques 15\n16 16 CHAPITRE I. RAYONS COSMIQUES 1 Les rayons cosmiques 1.1 Historique Après la découverte de la radioactivité en 1896 par Henri Becquerel, la communauté scientifique pensait que l ionisation de l air contenu dans des électroscopes était due à des éléments radioactifs dans la Terre ou dans l atmosphère. De 1911 à 1913, plusieurs voyages en ballon permettent à Victor Hess de mesurer une augmentation du taux d ionisation dans l air avec l altitude, attribuable à une irradiation d origine extra-terrestre. Cette découverte lui vaudra le prix de Nobel de Physique en En 1926, Millikan baptise ce rayonnement : "rayons cosmiques". Millikan croit alors que les rayons cosmiques sont des photons de très haute énergie produisant des électrons par effet Compton, alors que Compton pense que les rayons cosmiques sont issus de particules primaires chargées. Les expériences de la décennie lui donneront raison . Les rayons cosmiques fournissent alors une source naturelle de particules énergétiques aux physiciens. Ainsi les muons seront découverts en 1938 par Carl David Anderson et les pions en 1947 par Cecil Powell . En 1938, Pierre Auger met en évidence l existence de grandes gerbes atmosphériques en détectant l arrivée en coïncidence de particules sur des détecteurs espacés de quelques dizaines de mètres au sol. Il en déduit que ces gerbes sont produites dans la haute atmosphère, initiant une cascade d interactions secondaires. En 1948, Gottlieb et Van Allen montrent que les rayons cosmiques sont essentiellement des protons, quelques noyaux d hélium et une petite proportion de noyaux plus lourds. Depuis, certaines observations ont mis en évidence l existence de particules dont l énergie dépasse 20 ev. L étude du rayonnement cosmique n est pas close. De nombreux points, comme l origine, la composition, les mécanismes d accélération et la propagation des rayons cosmiques sont encore incompris. L observation du ciel avec différents messagers (photons, neutrinos, particules chargées) permettra l amélioration de la connaissance dans ce domaine. 1.2 Le spectre en énergie Parmi toutes les caractéristiques décrivant les rayons cosmiques, la plus impressionnante est la régularité de son spectre d énergie en loi de puissance, dn de E α (avec un indice spectral α compris entre 2,7 et 3,0), sur près de 12 ordres de grandeur en énergie (figure I.1). Différentes régions du spectre présentent un intérêt particulier pour la physique et l astrophysique du rayonnement cosmique . Basse énergie et modulation solaire : Les rayons cosmiques d énergie comprise entre 0 MeV et quelques GeV sont de loin les plus nombreux. Ils participent le plus intensément à l équilibre du milieu interstellaire et conditionnent la nucléosynthèse des éléments légers dans l univers. Cette composante du spectre nous est en partie masquée par l effet répulsif du vent solaire, dont les champs magnétiques repoussent les particules peu énergétiques. Les rayons cosmiques à ces énergies sont largement influencés par l intensité de l activité magnétique solaire. On observe une modulation périodique de leur intensité en opposition de phase avec le cycle du Soleil de période 11 ans . Energies intermédiaires : Aux énergies légèrement supérieures, jusque vers 0 GeV, les sections efficaces d interaction des noyaux énergétiques avec le milieu ambiant sont suffisamment grandes pour affecter les compositions des rayons cosmiques. Des études sur les rapports d abondance des noyaux et sur la dépendance de ces rapports avec l énergie donnent des informations sur la propagation des rayons cosmiques dans la galaxie. Ces "rayons" témoignent des conditions dans lesquelles s effectue le transport des noyaux depuis leurs sources jusqu aux détecteurs. Lorsque l énergie augmente, on perd de\n17 1. LES RAYONS COSMIQUES 17 FIGURE I.1 Spectre d énergie global des rayons cosmiques, donnant leur flux différentiel en fonction de l énergie l information sur la propagation des rayons cosmiques mais on obtient un accès plus direct aux rayons cosmiques à la sortie des processus d accélération. Le domaine d énergies comprises entre 0 GeV et le PeV correspond à l énergie maximale qui peut être conférée aux protons par des restes de supernovae dans les modèles standards. Les résultats expérimentaux indiquent que les restes de supernovae seraient susceptibles d accélérer des protons jusqu au PeV (ou Z fois plus haut pour les noyaux plus lourds que l hydrogène). Ces résultats seraient plutôt en faveur d un reste de supernova unique, proche et récent, apportant sa contribution au niveau du genou. Le genou : Autour de ev se situe la zone du "genou" (rupture de pente dans le spectre). Dans cette gamme d énergie, les flux sont tellement faibles (une particule par m 2 et par an) que l on ne peut faire de mesure directe en détectant la particule primaire : on doit observer la gerbe produite par l interaction des rayons cosmiques avec l atmosphère (détection indirecte). Les hypothèses généralement avancées pour expliquer la rupture de pente sont les suivantes : la modification des conditions de propagation des rayons cosmiques (origine astrophysique) l ouverture d un nouveau canal d interaction (nouvelle physique au TeV) la présence d une source astrophysique dans notre proche voisinage qui ajouterait une contribution en cet endroit précis du spectre\n18 18 CHAPITRE I. RAYONS COSMIQUES L énergie du genou varie selon le noyau considéré. Une étude détaillée du spectre énergétique pour différents noyaux permettrait d identifier si l énergie du genou est proportionnelle à la charge Z ou à la masse A du noyau. D un côté si l énergie est proportionnelle à Z alors une origine astrophysique du genou est attendue, puisque les champs magnétiques influent de la même façon sur la trajectoire des particules de même rayon de Larmor (α1/z). D un autre côté si un nouveau canal d interaction s ouvre, on s attend à ce que que l énergie soit proportionnelle à A. 1 Autour de 17 ev un second "genou" pourrait avoir été observé dans le spectre. L indice spectral deviendrait 3,3. Ce changement de régime n est aujourd hui pas expliqué . La cheville : Comme expliqué auparavant, les rayons cosmiques de basse énergie sont d origine galactique. Aux alentours de ev une autre rupture de pente apparaît et l indice spectral devient 2,8. Le changement de pente peut s expliquer par le passage d une composante dominante à basse énergie à une composante, de pente plus faible, finissant par dominer à haute énergie. A haute énergie le rayon de giration des particules dans les champs magnétiques est plus grand que la taille de la galaxie. Les rayons cosmiques ne peuvent alors plus être confinés dans la galaxie. La cheville représente la transition entre une origine galactique du rayonnement cosmique, dont les noyaux lourds sont le dernier vestige, et une origine extragalactique à plus haute énergie dominée par les protons. Il existe un autre de modèle, expliquant cette rupture de pente, basé sur la production de paires (e + e ) des protons extragalactiques avec le fond diffus cosmologique.  Le domaine GZK : Les rayons cosmiques d énergies de l ordre de 20 ev ne sont plus confinés dans la galaxie. Les objets envisagés pour produire des énergies supérieures à 20 ev seraient des objets extragalactiques comme les noyaux actifs de galaxie ou des sursauts gamma. L interaction des particules d une énergie de 20 ev avec les photons du fond de rayonnement cosmologique dans le domaine des micro-ondes (Cosmological Microwave Background -CMB-) diminue l énergie des rayons cosmiques alors qu ils se propagent. Si les rayons cosmiques sont des protons, il devrait y avoir une photoproduction de pions, si les noyaux sont plus lourds, le processus dominant serait la photodésintégration. Il en résulte une limitation de la distance que peuvent parcourir les rayons cosmiques ultra-énergétiques : quelques dizaines de Mpc pour les énergies supérieures à 20 ev. Une réduction du flux des rayons cosmiques est attendue au dessus de quelques 19 ev, c est la coupure Greisen-Zatsepin-Kutzmin [8, 9]. L expérience AUGER a confirmé une coupure progressive au-delà de ev , comme le montre la figure I Composition Jusqu à quelques centaines de TeV, les rayons cosmiques sont composés d environ 1% d électrons, dont l abondance relative décroît avec l énergie, et de 99% de noyaux. L existence de plusieurs techniques d observation des rayons cosmiques a permis de s intéresser séparément aux différents types de particules. Les protons et les noyaux Les protons représentent 89% des noyaux, les noyaux d hélium %, et les noyaux plus lourds le pourcent restant. La figure I.3 représente les abondances relatives des différents noyaux dans les rayons cosmiques comparées à celles du système solaire à basse énergie (E<270 MeV). La fraction relative de ces éléments est à peu près semblable dans les rayons cosmiques et dans le système solaire à deux différences près  : une sous représentation des éléments H et He, peut-être due à la difficulté d ioniser ces éléments, rendant les atomes accélérés plus rares. 1. L énergie peut être proportionnelle à A dans certains modèles astrophysiques (cannonball) .\n19 1. LES RAYONS COSMIQUES 19 FIGURE I.2 En haut : flux en fonction de l énergie avec ses erreurs statistiques. En bas : comparaison d Auger et d Hires avec un flux d indice spectral 2,69 . FIGURE I.3 Abondance des éléments du rayonnement cosmique galactique (E<270 MeV) (trait pointillé) comparés à celles du système solaire (trait plein) \n20 20 CHAPITRE I. RAYONS COSMIQUES une surreprésentation des éléments Li, Be et B et des éléments au-dessous du pic du fer, provenant de la fragmentation du rayonnement primaire sur le milieu interstellaire. La nature des rayons cosmiques à haute énergie ne peut être déterminée que par des mesures indirectes dans lesquelles on détecte les particules secondaires de la gerbe. La nature des rayons cosmiques est indéterminée aujourd hui : il existe différents modèles et leur abondance relative à haute énergie est controversée. Les électrons Le spectre des électrons au-delà de l atmosphère suit également une loi de puissance, mais l indice spectral est de 3,3. Le flux d électrons/positons est toujours inférieur à celui des protons et le rapport e /p diminue avec l énergie . Étant chargés et légers, les électrons, de basse énergie (E < GeV), sont particulièrement sensibles aux champs créés par les vents solaires et aux champs magnétiques du milieu interstellaire et sont sujets à l émission synchrotron. On considère que les électrons détectés ont deux provenances : les électrons primaires accélérés depuis une "source". Dans ce cas, les électrons, étant soumis à l émission synchrotron dans le champ galactique, mais également à la diffusion Compton inverse avec les photons mous ambiants, perdent très rapidement leur énergie. des électrons secondaires créés par l interaction de noyaux lors de leur propagation dans la Galaxie. Ces derniers sont alors liés au flux de protons et de noyaux cosmiques. Les photons L astronomie gamma est un des domaines les plus développés dans le domaine des rayons cosmiques. Du fait de leur neutralité, les photons se propagent "en ligne droite" depuis leurs sources jusqu à nous. Leur trajectoire nous renseigne donc sur leur origine, ce qui permet de faire des spectres séparément pour les différents types de sources. Le spectre des sources est obtenu par combinaison de données de plusieurs détecteurs, permettant ainsi d étudier l émission des sources sur différents domaines énergétiques (des ondes radio aux émissions γ). Aux plus hautes énergies, la résolution angulaire est moindre mais celle du télescope HESS est suffisante pour effectuer une cartographie des sources (paragraphe 1.8.f). Les photons γ proviennent principalement des phénomènes astrophysiques les plus violents (figure I.4) tels que : les supernovae, les étoiles à neutrons, les trous noirs stellaires, les noyaux actifs de galaxie (AGN), les sursauts γ... L étude de leurs spectres a permis d observer, pour chaque type de source, des propriétés communes que l on retrouve systématiquement, aboutissant à une classification des sources. Cette classification est ensuite utilisée pour identifier les nouvelles sources observées. Il existe des sources sans contrepartie qui ne sont pas identifiées. 1.4 Accélération des rayons cosmiques Les sources ponctuelles connues ne sont pas suffisantes pour rendre compte du flux total des rayons cosmiques mesuré. Il est possible qu un type de source particulier domine tous les autres par sa puissance et produise l essentiel des rayons cosmiques de l univers. A des énergies intermédiaires, ce seraient les SNRs (restes de supernovae). A très haute énergie, d autres sources extragalactiques pourraient entrer en jeu comme les noyaux actifs de galaxie abritant un trou noir supermassif en leur centre ou encore les sources de sursauts gamma. Il existe différentes théories pour expliquer la production des rayons cosmiques de haute énergie, que l on peut classer en deux catégories : les mécanismes astrophysiques, appelés également "Bottom-Up", où l on part d une particule peu énergétique que l on accélère dans des objets astrophysiques tels que les noyaux actifs de galaxies, les sursauts gamma ou les étoiles à neutrons.\n21 1. LES RAYONS COSMIQUES 21 FIGURE I.4 Sources ponctuelles (E > 0 MeV) du troisième catalogue d EGRET en coordonnées galactiques. les mécanismes "Top-Down", dits exotiques, où l on part d une particule "super-massive" (M > 21 ev) qui se désintègre en particules secondaires connues ultra-relativistes. Les modèles de type "Bottom-Up" L accélération des particules dans les objets astrophysiques met en jeu des mécanismes de transfert d énergie de l objet astrophysique vers la particule. Seuls les champs magnétiques élevés permettent d atteindre de hautes énergies. Il existe deux modèles utilisant de grands champs électromagnétiques  : le premier accélère les particules chargées dans un seul gradient. Pour atteindre une énergie de 20 ev, le potentiel électrique doit être de l ordre de 20 Z Volts où Z est le nombre atomique de la particule accélérée. Actuellement seuls les pulsars (étoiles à neutrons en rotation rapide) peuvent présenter un tel champ. le second repose sur le mécanisme de Fermi. Enrico Fermi proposa en 1949  un mécanisme où les particules diffusées de manière stochastique sur les irrégularités magnétiques d un nuage de gaz en mouvement à une vitesse V par rapport à la source gagnent une énergie proportionnelle à V/C à chaque collision et statistiquement l énergie totale gagnée est proportionnelle à (V/c) 2, où c est la vitesse de la lumière. C est ce qu on appelle le mécanisme de Fermi du second ordre. L intérêt est que cela conduit naturellement à un spectre différentiel en loi de puissance. Le problème est qu étant donné les vitesses très faibles des nuages de gaz interstellaire et la faible densité de ceuxci, l efficacité du mécanisme est insuffisante. Quand on considère le même phénomène, mais dans des ondes de choc, on obtient un gain d énergie en V/c où V est maintenant la vitesse du front de l onde de choc. C est le mécanisme de Fermi du premier ordre. Celui-ci conduit naturellement à une loi de puissance avec un indice spectral de -2. Ces ondes de choc apparaissent naturellement lorsque la matière éjectée par une supernova rencontre le gaz environnant. Les modèles de type "Top-Down" Dans ce type de modèles, on suppose que les rayons cosmiques sont produits par la désintégration de particules super massives et métastables ayant des masses supérieures à 21 ev, dont la durée de vie est de l ordre de l âge de l univers . Les particules produites selon ces modèles sont des hadrons, des photons et des neutrinos. Un télescope à neutrinos serait donc\n22 22 CHAPITRE I. RAYONS COSMIQUES bien placé pour donner des limites sur ce genre de modèles. 1.5 Origine des rayons cosmiques La plupart des sources témoignent de spectres d énergie en loi de puissance avec des indices spectraux entre -1,8 et -2,4, indices similaires à celui du rayonnement cosmique à sa source. Pour qu une particule acquiert une grande énergie, elle doit être confinée dans une région d accélération . Or plus l énergie d une particule augmente, plus il est difficile pour les champs magnétiques qui les confinent de courber leur trajectoire. Une limite supérieure à l énergie maximale atteinte en un site quelconque d accélération peut être obtenue en exigeant que le rayon de Larmor des particules, r L = E/qBc, soit inférieur aux dimensions du site, R. Les quantités B et R sont obtenues par la relation : E max = ZecBR 21 ev On peut remarquer que si on définit le champ électrique E = Bc comme le champ effectif maximal associé à B, alors E max représente le travail de ce champ sur R. Cependant dans le cas d une accélération par onde de choc, l énergie maximale est donnée par E max = ZeBV c R où V c est la vitesse du choc. On voit alors que le champ électrique induit E = BV c est mis à profit pour accélérer les particules. Si l accélérateur présente un mouvement relativiste de facteur de Lorentz Γ, le champ électrique induit sera de la forme E = ΓBV c. Par conséquent, pour atteindre des énergies élevées, il faut se tourner vers les sources les plus vastes ou dans lesquelles on rencontre des différences de vitesse élevées ou de grands champs magnétiques. La figure I.5 représente différents types de sources énergétiques de l univers dans un diagramme (taille de l objet, champ magnétique de cet objet). A haute énergie, les pertes énergétiques ne sont pas prises en compte. Ces pertes peuvent être distinguées en deux parts : les pertes énergétiques par rayonnement synchrotron des électrons ou des protons et les pertes par photoproduction pγ π ou par production de paires pγ e e Propagation Une fois ces particules injectées dans le milieu interstellaire, elles vont se propager dans la galaxie. Les particules vont alors  : se disperser en raison des inhomogénéités du champ magnétique galactique. Elles perdent alors toute information de direction du rayonnement, interagir avec les nuages de gaz interstellaires en créant des particules secondaires par spallation, en perdant de l énergie par ionisation ou en étant éventuellement absorbées, être transportées par le vent galactique, se désintégrer, être réaccélérées de façon stochastique sur les inhomogénéités magnétiques (mécanisme de Fermi du second ordre). Ce modèle général aboutit à une équation non linéaire, ne présentant pas de solutions analytiques dans le cas général. Le modèle le plus simple pour prendre en compte ces différents éléments est alors le modèle de la boîte qui fuit (Leaking Box Model). Le modèle de la boîte qui fuit Dans ce modèle, les rayons cosmiques se propagent librement dans un volume avec une probabilité constante de s échapper de ce volume. Les convections ne sont pas pris en compte. Ce modèle, qui a le mérite d être très simple et de dépendre de peu de paramètres, rend compte assez fidèlement des flux de particules primaires et des particules secondaires qui leurs sont associées. Il\n23 1. LES RAYONS COSMIQUES 23 FIGURE I.5 Diagramme de Hillas pour les sites potentiellement capables d accélérer des rayons cosmiques aux énergies extrêmes. Si on fixe une charge et une énergie, les couples (taille caractéristique, intensité B du champ magnétique) possibles forment alors une droite sous cette représentation (formule 1.5). Les lignes rouges et vertes représentent l application de ce critère aux protons de 20 ev et 21 ev et à des noyaux de fer de 20 ev. ne permet pas l étude des primaires d origine exotique car ceux-ci sont très sensibles au profil de densité de matière, notamment dans la région centrale de la galaxie. 1.7 Et les neutrinos? Les γ de très haute énergie, de part leur propagation en ligne droite, ont permis d identifier récemment de nombreuses sources (voir paragraphe 1.8.f). Cependant les photons interagissent avec la matière interstellaire (nuages de poussière ou de gaz) et avec les rayonnements électromagnétiques par création de paires au dessus du TeV. L observation du ciel avec les photons est alors limitée en énergie et en distance. A contrario, le neutrino, qui n interagit que via l interaction faible, peut parcourir des distances cosmologiques sans altération (figure I.6). Comme nous l avons déjà observé précédemment, les protons issus du rayonnement cosmique sont déviés par les champs magnétiques galactiques et intergalactiques et ne permettent pas de déterminer facilement la source de ce rayonnement pour des énergies inférieures à 20 ev. On pourrait alors penser au neutron, mais sa courte durée de vie (866 s) l empêche de parcourir de longues distances. Comme le montre la figure I.7, l observation par neutrinos est possible sur plus de 20 ordres de grandeur. Les télescopes à neutrinos sont sensibles à partir de GeV pour l observation des neutrinos atmosphériques, et à partir du TeV pour les neutrinos issus de sources astrophysiques.\n24 24 CHAPITRE I. RAYONS COSMIQUES FIGURE I.6 Panorama des messagers permettant de sonder l Univers sur des grandes échelles en fonction de l énergie. Les neutrinos couvrent une large zone du spectre permettant les observations astronomiques. FIGURE I.7 Spectre des flux diffus de neutrinos issus de certaines catégories de sources. Dans les mécanismes Bottom-Up, les neutrinos sont issus de la collision de noyaux avec d autres noyaux ou photons de plus faible énergie. Ils ne peuvent être produits que par des processus hadroniques. Des mésons neutres et chargés (des pions principalement) sont alors produits et se désintègrent selon le mécanisme suivant : p + A p + γ } π 0 + π + + π +... γ + γ µ + + ν µ µ + ν µ e + + ν e + ν µ e + ν e + ν µ Ce mécanisme requiert en général des champs magnétiques intenses susceptibles de piéger les particules accélérées (électrons, protons, ions) dans des zones de confinement. Les électrons perdent leur énergie par rayonnement synchrotron et produisent ainsi des photons de basse énergie qui constituent alors une cible pour les protons de haute énergie. Les réactions photo-hadroniques nécessitent de fran-\n25 1. LES RAYONS COSMIQUES 25 chir un seuil en énergie. Les pions produits sont donc quasiment au repos dans le référentiel du plasma de particules qui peut-être, lui, ultra-relativiste. Le flux de neutrinos issus de ces désintégrations peut être établi à partir de la relation différentielle  : ( dn ν = N(E ν) de ν 1 Z NN A π,ν 1 + B π,ν cos θ Eν ǫ c π + 0,635 A K,ν 1 + B K,ν cos θ Eν ǫ c K où I représente le flux initial θ est l angle d incidence par rapport à la cible A, B et Z NN sont des constantes caractéristiques des interactions ǫ c N est l énergie critique pour les pions ou les kaons. Elle donne une estimation de l énergie au dessus de laquelle la probabilité d interaction prend le dessus sur celle de la désintégration. X prend en compte la désintégration des mésons de saveurs plus lourdes que les pions et les kaons. Le premier terme reflète la désintégration des pions en neutrinos, alors que le second terme représente la désintégration des kaons. Le rapport des flux de neutrinos Φ ν+ν et de photons Φ γ est variable suivant l indice spectral du primaire et le type de l interaction. Les simulations permettent d aboutir à la relation indépendante de l énergie : Φ ν+ν = AΦ γ où A 0, 5 dépend de l indice spectral . Etant donné la petitesse de la section efficace des neutrinos (de 11 à 9 barns pour des énergies comprises entre le TeV et le PeV), on comprend que l astronomie neutrino est un défi qui nécessite des détecteurs à l échelle de 1 km 3. ANTARES est le plus grand détecteur dans l hémispère Nord mais ses dimensions restent encore modestes (0,01 km 2 au sol). + X ) (I.1) (I.2) 1.8 Sources ponctuelles éventuellement productrices de neutrinos La détection de neutrinos provenant d une source ponctuelle permettrait de trouver les sources émettrices de rayons cosmiques. Les sources les plus diffuses émettraient des neutrinos comme conséquence de l interaction des protons du rayonnement cosmique avec un milieu étendu. Parmi les sources galactiques susceptibles d engendrer des neutrinos se trouvent des restes de supernovae (supernova remnants -SNRs-), qui pourraient être les sources d accélération des rayons cosmiques jusqu au genou, les microquasars et les pulsars, qui seraient d excellents candidats pour la production de rayons cosmiques d énergie inférieure à 0 TeV. Les sources extragalactiques pourraient être les noyaux actifs de galaxie (Active Galactic Nuclei -AGN-) ou les sursauts gamma (Gamma Ray Bursts -GRBs-).  1.8.a Les supernovae et leurs restes L observation des restes de supernova au TeV en neutrinos est la signature d une désintégration de π ±. On attend de 1 à événements dans un détecteur de taille kilométrique avec un seuil en énergie de 1 TeV, avec un bruit de fond provenant des neutrinos atmosphériques du même ordre de grandeur dans la ligne de visée.  Les explosions de supernova, qui devraient arriver tous les 30 ans environ dans notre galaxie, produisent des neutrinos à des énergies autour du MeV (ce qui a été observé lors de l explosion de SN1987A).\n26 26 CHAPITRE I. RAYONS COSMIQUES 1.8.b Les microquasars Les microquasars  sont des systèmes binaires émettant des jets colinéaires observés en rayons X. Ces émissions seraient dues à des instabilités dans le disque d accrétion du système lorsqu une grande quantité de matière s est accrétée autour de l étoile compagnon (figure I.8). Divers modèles prédisant des processus hadroniques produisant des neutrinos ont été proposés, dont : [20, 21] : des inhomogénéités dans le jet causent des chocs internes qui accélèrent les électrons et les protons leur donnant une distribution en énergie en loi de puissance. Les protons accélérés peuvent interagir aussi bien avec les photons X émis par le disque d accrétion qu avec les photons synchrotron produits dans le jet par les électrons accélérés  : dans des systèmes binaires massifs très compacts, des noyaux accélérés dans le jet perdent des neutrons dans les interactions avec le rayonnement stellaire ou du disque d accrétion. Ces neutrons se propagent vers le disque et l étoile massive, produisant des neutrinos par interaction dans  une anti-corrélation entre gammas et neutrinos de haute énergie est prédite pour le microquasar variable LS 5039 observé par HESS . La région serait optiquement mince pendant la période haute de gammas, tandis que la période basse serait due à une région absorbante pour les photons, servant de cible pour la production de neutrinos. 1.8.c Les noyaux actifs de galaxie Les noyaux actifs de galaxie (Active Galactic Nuclei -AGN-) sont des objets contenant un trou noir super-massif en rotation, situé au centre des galaxies. Le trou noir est entouré d un disque d accrétion visible en optique. Le disque est alimenté par un tore de poussière. Perpendiculairement au disque d accrétion, deux jets relativistes sont émis, visibles en radio (figure I.9). Ils sont des candidats pour l accélération des rayons cosmiques du fait de la grande quantité d énergie émise. Dans le disque d accrétion, des particules sont accélérées et peuvent alors interagir avec la matière ou avec les photons. De plus, les champs magnétiques dans les plasmas entourant le trou noir central piègent en partie les photons et les protons alors que les neutrinos s en échappent . 1.8.d Les sursauts gamma Les sursauts gamma (Gamma Ray Burst -GRB-) sont de courtes et très intenses émissions de photons γ. On classe généralement les sursauts suivant deux types : les longs et les courts (moins de 2 s). Le décalage vers le rouge élevé des galaxies hôtes des sursauts gamma et la distribution isotrope sur la carte du ciel de BATSE (figure I.) excluent une origine galactique des sursauts gamma. Trois phases d émission sont attendues : les heures précédant le GRB, la phase du sursaut, et l émission après le GRB. Quelques heures avant l émission principale, le sursaut est opaque aux photons et des neutrinos peuvent être produits par interaction des protons accélérés avec des rayons X. Lors du sursaut, des protons seraient accélérés par des interactions entre particules ultra-relativistes, ce qui induirait des productions de neutrinos. Après le sursaut principal, des neutrinos seraient produits en corrélation avec la production électromagnétique. Deux modèles tentant d expliquer les sursauts gamma impliquent une émission de neutrinos : le modèle dit de "la boule de feu" (fireball model ), le plus généralement utilisé, et celui plus controversé du "boulet de canon" (canon ball model ). Le premier se base l hypothèse d interactions entre les particules de plasmas relativistes, alors que le second suppose des corps denses relativistes, produit juste après une supernova, venant interagir avec la matière expulsée par l explosion de l étoile. Les neutrinos émis dans le second modèle sont de plus basse énergie que ceux émis dans le premier modèle.\n27 1. LES RAYONS COSMIQUES 27 FIGURE I.8 Modèle de microquasar FIGURE I.9 Classification des différents types de noyaux actifs de galaxie en fonction de leur intensité et de la ligne de visée de l observateur.\n28 28 CHAPITRE I. RAYONS COSMIQUES FIGURE I. Carte des sursauts gamma détectés par BATSE, en coordonnées galactiques 1.8.e La matière noire L annihilation de particules de matière noire dans des corps compacts émettrait des neutrinos . Lorsque des particules comme les neutralinos ou les WIMPs sont piégées dans un corps dense par la force gravitationnelle, ils peuvent s annihiler en des particules dont la désintégration donnerait des neutrinos de haute énergie. La détection indirecte de matière noire en direction du soleil prouverait l existence de ces particules. L annihilation de particules reliques des premiers instants de l Univers participe au flux diffus de neutrinos et un détecteur de taille kilométrique permettrait de contraindre les modèles concernant cette production de neutrinos. 1.8.f Etude des sources à l aide des gamma : l expérience H.E.S.S. Les différentes sources citées ainsi que des modèles de matière noire ont été étudiés par l expérience H.E.S.S. Nous présenterons ici les derniers résultats obtenus par H.E.S.S. en ce qui concerne l étude du centre de notre galaxie. Le détecteur H.E.S.S. est un détecteur au sol de lumière Tcherenkov. La lumière Tcherenkov est émise dans l atmosphère par les particules chargées dont la vitesse est supérieure à celle de la lumière dans l air. L indice de l air variant avec la densité de l air, les seuils d émission et l angle du cône de lumière varient avec l altitude. Au niveau de la mer, le seuil en énergie est de 22 MeV pour les électrons et 4,5 GeV pour les muons et l angle du cône de l ordre du degré. Une gerbe de particules émet de la lumière Tcherenkov et éclaire un détecteur au sol sur plusieurs dizaines de mètres 2. Le détecteur H.E.S.S. utilise 4 grand miroirs de 13 m de diamètre équipés au plan focal d une caméra de haute définition, constituée de 960 photomultiplicateurs. Les images obtenues permettent de remonter à la direction des photons. H.E.S.S. a fait de nombreuses cartes des différentes sources observées et notamment dans la direction du centre galactique (figure I.11). Les deux sources non identifiées d Egret apparaissent proches des émissions détectées par H.E.S.S. L émission le long du plan galactique, d environ 320 pc, est en très bonne corrélation avec les distributions de gaz moléculaires dans la région. Les corrélations observées entre les γ de H.E.S.S. et les nuages moléculaires suggèrent une origine hadronique du rayonnement γ observé.\n29 2. MÉTHODES DE DÉTECTION DES RAYONS COSMIQUES 29 FIGURE I.11 Carte en coordonnées galactiques de notre galaxie vue par H.E.S.S. En haut : vue globale du centre galactique avec le nombre de photons détectés. En bas : même carte après soustraction des deux sources dominantes, les contours blancs représentent la densité de gaz moléculaire. Le cercle jaune représente SNR G et l étoile noire la position du centre galactique Sgr A. Les régions (à 95% de niveau de confiance) des deux sources non identifiées d Egret sont indiquées en lignes vertes pointillées.  2 Méthodes de détection des rayons cosmiques Comme montré à la figure I.1, les mesures effectuées sur les rayons cosmiques dépendent fortement du flux : il est donc nécessaire d utiliser différents types d instruments (figure I.12). A basse énergie, le flux est suffisamment important pour permettre une observation directe des rayons cosmiques dans la haute atmosphère : on utilise des détecteurs placés dans des satellites ou dans des ballons (BESS 2.1). Lorsque l énergie augmente, le flux de particules diminue. Au delà de la taille des satellites et des détecteurs en ballon, des grands détecteurs au sol, KASCADE-Grande (paragraphe 2.2.a) et AUGER (paragraphe 2.2.b) observent les particules secondaires produites par la gerbe (détection indirecte des rayons cosmiques). Quelques expériences seront détaillées par la suite : l expérience BESS (Balloon-borne Experiment with a Superconducting Spectrometer) est une expérience en ballon dont le premier vol eu lieu en 1993 et le dernier en décembre BESS-Polar permet de mesurer une énergie comprise entre 0,2 et 20 GeV par nucléon. l expérience KASCADE-GRANDE (KArlsruhe Shower Core and Array DEtector and Grande array) étudie les gerbes de particules secondaires au niveau du sol. Les énergies couvertes par cette expérience sont comprises entre 14 et 18 ev. l expérience AUGER est aussi un détecteur au sol et fonctionne depuis Cette expérience se concentre sur les particules dont l énergie est supérieure à quelques 18 ev.\n30 30 CHAPITRE I. RAYONS COSMIQUES FIGURE I.12 Ordre de grandeur des intensités du rayonnement cosmique primaire (toutes particules chargées) à différentes énergies. Domaine respectif des différentes techniques de détection. l expérience souterraine MACRO (Monopole, Astrophysics, and Cosmic Ray Observatory) détecte des muons dans le domaine d énergie d ANTARES. 2.1 L expérience BESS La détection au dessus de l atmosphère permet l accès direct au rayon cosmique. On a donc une seule particule à identifier. Pour cela, on utilise les détecteurs de physique des particules pourvu qu ils soient suffisamment légers et petits pour être placés sur des satellites, des ballons atmosphériques ou des navettes spatiales. Une des techniques les plus légères, utilisée dans de nombreuses expériences en ballon, utilise un spectromètre magnétique (figures I.13 et I.14). Un premier trajectographe enregistre la direction et le temps d arrivée de la particule qui traverse ensuite un champ magnétique où elle est déviée. Un deuxième trajectographe permet alors de mesurer la rigidité R = pc/ze (où Z est la charge) et le temps de vol (qui donne β) de la particule dans le champ magnétique. Ce système est complété par un ou plusieurs détecteurs permettant de mesurer la charge du rayon cosmique : des compteurs proportionnels, des scintillateurs ou des chambres à ionisation, qui donnent la perte d énergie différentielle de/dx en fonction de β et de Z 2 . Le principal inconvénient de ce détecteur est qu il est limité à un intervalle d énergie très réduit. Une particule de trop faible énergie va avoir un rayon de giration très petit, et ne va pas traverser le deuxième trajectographe. La résolution du trajectographe est limitée par la taille du rayon de giration liée à l énergie du rayon cosmique. L étude des plus hautes énergies se fait soit en augmentant la résolution des trajectographes, soit en augmentant l intensité du champ magnétique. De plus, ce type de détecteurs n offre pas la même sensibilité entre les différents noyaux qui ont des charges différentes. Lors des mesures en vol de l expérience BESS, un champ magnétique de 1 Tesla engendré par un aimant supraconducteur a été utilisé. Le spectromètre magnétique de forme cylindrique était constitué de scintillateurs qui servaient de trajectographe et permettaient d identifier la nature de la particule (par la perte d énergie de/dx) . Un des résultats les plus importants est la mesure du flux de proton et d hélium à des énergies comprises entre 0,215 et 21,5 GeV par nucléon (figure I.15). Les données utilisées pour ces mesures ont été prises entre 1997 et 2002 avec le spectromètre BESS à bord d un ballon lancé à Lynn Lake au Canada .\n31 2. MÉTHODES DE DÉTECTION DES RAYONS COSMIQUES 31 FIGURE I.13 Schéma d un spectroscope à champ magnétique couplé à un système qui mesure la charge : compteurs proportionnels, scintillateurs, chambres à ionisation, etc. FIGURE I.14 Vues transversales des spectromètres BESS-98 et BESS-TeV. Les lignes épaisses représentent la longueur de la trace. FIGURE I.15 Flux de proton et d helium en haut de l atmosphère obtenus avec les données de BESS entre 1997 et Le flux d helium a été divisé par un facteur. \n32 32 CHAPITRE I. RAYONS COSMIQUES 2.2 Les détecteurs au sol Les gerbes susceptibles d engendrer des muons suffisamment énergétiques pour atteindre le détecteur ANTARES possèdent une énergie supérieure à 1 TeV. Pour les détectern on installe des détecteurs au sol. Lorsqu un rayon cosmique pénètre dans l atmosphère, il interagit avec le milieu et forme une gerbe de particules. Plus l énergie de la particule initiale est grande, plus la gerbe engendrée est pénétrante et large au niveau du sol. Lors du développement de la gerbe, le nombre de particules secondaires augmente, atteint un maximum, puis diminue à cause de l absorption atmosphérique. A partir d une dizaine de TeV pour les hadrons (0 TeV pour les γ), les gerbes deviennent suffisamment pénétrantes pour atteindre le sol avant d être totalement absorbées. 2.2.a L expérience KASCADE-grande KASCADE-Grande est une expérience multi-détecteur optimisée pour des énergies entre 14 et 18 ev. Elle est composée de plusieurs détecteurs : Le réseau de cuves KASCADE (figure I.17) mesure le nombre total d électrons et de muons d énergie supérieure à 230 MeV de la gerbe avec un réseau de 252 stations. Chaque station possède des détecteurs avec et sans bouclier sur une surface de m 2. Le traceur de muon mesure l angle d incidence, relativement à la direction de la gerbe, des muons d énergie supérieure à 800 MeV. Les muons d énergie supérieure à 2,4 GeV sont détectés par une chambre à fils proportionnelle. Le réseau Grande est constitué de scintillateurs plastiques distants de 130 m. Il permet des mesures de rayons cosmiques jusqu à 18 ev. Issu de cette expérience, un programme de simulation CORSIKA (COsmic Ray SImulations for KAscade) a été développé. Il autorise de choisir un flux de rayons cosmiques primaires et un modèle d interaction hadronique régissant le développement de la gerbe dans l atmosphère. Ce programme sera utilisé dans l analyse présentée dans cette thèse. Le spectre des particules obtenu par KASCADE en prenant pour un hypothèse un certain modèle d interaction hadronique est comparé aux résultats d autres expériences (figure I.16). On observe un bon accord sur le flux mesuré et la position du genou entre les expériences KASCADE, HEGRA, et EAS-TOP. 2.2.b L expérience Auger L expérience Auger est un détecteur hybride composé de deux parties : une constituée de quatre détecteurs de fluorescence, et l autre de 1600 cuves d eau (figure I.18). La technique "oeil de mouche" (ou détecteur de fluorescence) consiste à collecter la lumière de fluorescence qu émettent les atomes d azote de l air au passage d un rayon cosmique. On peut observer le développement de la gerbe en détectant les photons ultra violets émis. L installation de deux télescopes à quelques kilomètres l un de l autre permet une observation stéréoscopique. On estime la direction de la particule incidente ainsi que son énergie à partir de l image spatio-temporelle de la gerbe. La nature du rayon cosmique est identifiée par l altitude du maximum de développement dans l atmosphère. En effet, un noyau lourd atteint son maximum plus tôt et donc plus haut qu un noyau léger de même énergie et de même incidence. L inconvénient de cette méthode de détection est que la lumière émise est très faible. De plus cette observation est sensible aux variations atmosphériques et ne peut se faire que lors des nuits sans lune. Les cuves d eau permettent d étudier les particules de la gerbe au sol. Des photomultiplicateurs installés dans la cuve permettent la mesure de l émission Tcherenkov dans l eau des particules chargées.\n33 2. MÉTHODES DE DÉTECTION DES RAYONS COSMIQUES 33 FIGURE I.16 Spectre des particules obtenu par Kascade avec un modèle d interaction hadronique est comparé aux résultats d autres expériences.  Plus l énergie de la gerbe est grande, plus la surface au sol de la gerbe est grande pour une direction donnée. Le diamètre de la gerbe au sol peut être de plusieurs kilomètres. Les résultats de l expérience ont permis de confirmer une coupure du spectre des rayons cosmiques à ev (figure I.2). Les 27 événements de plus haute énergie (E>57 EeV) issus des données révèlent des corrélations avec la matière extra-galaxtique locale, mise en évidence par la distribution des AGNs au sein du catalogue Véron-Cetty & Véron. Ce résultat indique une origine extra-galactiques des sources de rayons cosmiques de haute énergie. FIGURE I.17 Réseau de l expérience Kascade FIGURE I.18 Schéma et principe de détection de l expérience Auger avec ses cuves et ses détecteurs de fluorescence\n34 34 CHAPITRE I. RAYONS COSMIQUES FIGURE I.19 Carte du ciel en coordonnées galactiques représentant les 27 événements d Auger de haute énergie (E>57 EeV) avec une résolution angulaire de 3,2 et 442 AGNs de redshift z 0,017. Le code de couleur bleue indique le temps d exposition du ciel par l observatoire Pierre Auger.  2.3 Un détecteur souterrain : MACRO Le détecteur MACRO (Monopole, Astrophysics, and Cosmic Ray Observatory) était un détecteur de muons souterrain au laboratoire du Gran Sasso (figure I.20). Les dimensions du détecteurs étaient 76,5 12 9,3 m 3, sa masse d environ 5300 t, et son acceptance à un flux isotrope de particules de 000 m 2. Le principal intérêt de MACRO était la redondance et la complémentarité des trois détecteurs qui le composaient : des scintillateurs liquides et des tubes à avalanches pour la détection des muons, des détecteurs de traces nucléaires pour la détection des particules exotiques. FIGURE I.20 Vue en coupe du détecteur montrant la structure du détecteur  A l aide des scintillateurs et des tubes à avalanche, MACRO a reconstruit les traces de muons en obtenant leur énergie et leur direction. Si le muon donnait un événement dans la partie supérieure du détecteur, le grand bras de levier donnait une grande précision sur la direction du muon. La prise de données de MACRO dura 5 ans et MACRO collecta plus de 60 millions de muons descendants de haute énergie. MACRO a mesuré le flux de muons dans la roche en fonction de la profondeur (figure I.21) entre 3200 et 7000 hg.cm 2. Les muons utilisés correspondent à des énergies comprises\n35 2. MÉTHODES DE DÉTECTION DES RAYONS COSMIQUES 35 entre 1 et 20 TeV au niveau de la surface et à des énergies du rayon cosmique primaire comprises entre 2 et 200 TeV, ce qui correspond à la gamme d énergie d ANTARES. La recherche de sources ponctuelles par MACRO avec 5, 8 années effectives de prises de données a permis de mettre une limite sur le flux de muons induits par des neutrinos de 1,7 15 µ cm 2 s 1 . FIGURE I.21 Flux vertical de muon obtenu par MACRO. Le trait plein représente un ajustement aux données de MACRO, le trait long pointillé aux données du FREJUS, et le trait pointillé l ajustement de Crouch .\n36 36 CHAPITRE I. RAYONS COSMIQUES\n37 Chapitre II Télescopes à neutrinos 37\n38 38 CHAPITRE II. TÉLESCOPES À NEUTRINOS 1 Principe de détection Le neutrino est sensible uniquement à l interaction faible : d un côté cela lui permet de parcourir des distances cosmologiques sans être absorbé par la matière interstellaire ou de s échapper plus facilement des sites de production des rayons cosmiques ; de l autre cela rend sa détection difficile. Les difficultés techniques associées à la méthode de détection constituent plusieurs défis à relever. La détection des neutrinos est indirecte. Le détecteur est une matrice tridimensionnelle de photomultiplicateurs qui, plongés dans un milieu transparent, captent la lumière Tcherenkov émise par les muons issus de l interaction profondément inélastique de neutrinos muoniques avec les noyaux environnants. Le nombre de photons Tcherenkov touchant les photo-multiplicateurs, ainsi que leurs temps d arrivée permettent une reconstruction précise de la trajectoire des muons. Pour diminuer l influence du bruit de fond que constituent les muons atmosphériques (descendants) le détecteur doit être placé sous un épais bouclier. En outre, le signal est recherché parmi les traces montantes (seuls les neutrinos peuvent traverser la Terre) ouvrant ainsi une fenêtre d observation instantanée sur la moitié de la voûte céleste. Compte tenu de la dimension du détecteur (environ 0,01 km 2 au sol et 450 m de hauteur pour ANTARES), imposée par la faiblesse des flux de neutrinos cosmiques et celle de leur section efficace, la seule possibilité est l enfouissement profond dans la glace ou dans la mer. Le site retenu pour le télescope ANTARES se trouve par 2475 m de fond, au large de Toulon (Var). La récupération des données, les opérations de contrôle, la mise en service ainsi que l alimentation électrique sont assurées par un câble électro-optique d une quarantaine de kilomètres, depuis la station côtière. Le télescope comporte 12 lignes de détection équipées de 75 photo-multiplicateurs chacune. La construction du télescope s est échelonnée sur plusieurs années. La prise de données a débuté dès la connexion de la première ligne en février Le télescope est complet depuis mai Interaction du neutrino avec la matière Les neutrinos interagissent avec la matière par interaction faible suivant deux voies possibles : par échange de courant neutre ou par échange de courant chargé. La section efficace neutrino-électron est petite devant la section efficace neutrino-nucléon parce que la masse de l électron est très petite devant la masse d un nucléon. Nous n étudierons donc que l interaction neutrino-nucléon à haute énergie avec production de muons issus de cette diffusion profondément inélastique. L interaction par courant neutre du neutrino sur un nucléon N se note : ν l (ν l ) + N Z0 ν l (ν l ) + X (II.1) Où X est une gerbe hadronique que l on peut détecter. L interaction par courant chargé produit un lepton chargé de même saveur que le neutrino initial : ν l (ν l ) + N W ± l ( l +) + X (II.2) Comme le muon parcourt des distances plus grandes, il a plus de chances d être détecté que les autres leptons. Il permet en outre de reconstruite avec davantage de précision la direction du neutrino incident. Il permet en outre de reconstruire avec davantage de précision la direction du neutrino incident (et des astronomiques). On se concentre donc sur le canal muonique. Pour un muon, le processus est illustré par le diagramme de Feynman (figure II.1). Dans le cadre de la théorie électrofaible et du modèle des quarks-partons de Björken, la section efficace de diffusion profondément inélastique d interaction d un neutrino muonique par échange de\n39 1. PRINCIPE DE DÉTECTION 39 FIGURE II.1 Diagramme de Feynman de l interaction d un neutrino muonique sur un neutron par courant chargé. courant chargé s exprime comme  : d 2 σ dxdy = 2G2 F ME ( ν M 2 ) 2 [ W π Q 2 + MW 2 xq ( x,q 2) + xq ( x,q 2) (1 y) 2] (II.3) avec M et M W les masses du nucléon et du boson intermediaire, N un nucléon isoscalaire (N n+p 2 ), et x = Q2 2mν, Q2 étant le carré de la quadri-impulsion du parton, y = ν E ν, ν étant l énergie hadronique (fraction d énergie du neutrino emportée par le boson) dans le référentiel du laboratoire donc ν = E ν E µ et G F est la constante de Fermi. Les fonctions de distribution des partons sont : q ( x,q 2) = u ( v x,q 2 ) ( + d v x,q 2 ) + u ( s x,q 2 ) ( + d s x,q 2 ) ( + s s x,q 2 ) ( + b s x,q 2 ) (II.4) ( 2 2 x,q 2 ) ( + d s x,q 2 ) q ( x,q 2) = u s 2 + c s ( x,q 2 ) + t s ( x,q 2 ) (II.5) où les indices v et s sont les contributions de valence et de la mer, et u, d, c, s, t, b les distributions des différents quarks dans un proton. 1.2 Propagation dans la Terre L augmentation de la section efficace d interaction du neutrino avec l énergie favorise la détection. Cependant cet effet est limité : la probabilité d interaction des neutrinos avec la Terre à des énergies de 0 TeV est proche de 1, la terre devient opaque. A 70 TeV, la longueur d absorption du neutrino est égale au diamètre de la Terre et les interactions des neutrinos ne sont plus détectées. Pour quantifier cet effet on prend en compte le profil de densité terrestre  illustré à la figure II.3. L épaisseur (en g.cm 2 ) de matière traversée par le neutrino est alors calculée en intégrant la densité des dix couches du profil de densité terrestre sur la hauteur de chacune de ces couches (figure II.4). Il est ensuite possible de calculer numériquement l absorption des neutrinos par la Terre en fonction de l énergie et de l angle zénithal (figures II.5 et II.6). L absorption par la Terre est négligeable pour les neutrinos horizontaux (angle zénithal proche de 90 ) car ils traversent peu de matière. La Terre étant opaque aux neutrinos de ultra hautes énergies, ces neutrinos sont visibles essentiellement à l horizontal. Par contre les neutrinos de plus basse énergie n interagissent pas avec la Terre et seront peu visibles à ces angles.\n40 40 CHAPITRE II. TÉLESCOPES À NEUTRINOS FIGURE II.2 Section efficace des interaction neutrino/nucleons à hautes énergies : la ligne en pointillés correspond aux échanges par courant neutre, la ligne fine aux échanges par courant chargé et l autre ligne est la section efficace totale . FIGURE II.3 Profil de densité terrestre utilisé dans . Dans ce modèle, la Terre est schématisée par couches concentriques dont la densité dépend de la distance radiale au centre de la terre. FIGURE II.4 Epaisseur à traverser pour atteindre le détecteur en fonction de l angle nadir.\n41 1. PRINCIPE DE DÉTECTION 41 Coefficient moyen d' absorption ds la terre cos(θ) FIGURE II.5 Effet moyen de l absorption de neutrinos d énergie de 1, et 0 TeV en fonction de leur angle zénithal. FIGURE II.6 Absorption moyennée pour les neutrinos montants (en noir) et pour les deux hémisphères (en rouge) en fonction de l énergie. 1.3 Propagation des muons induits 1.3.a Corrélation neutrino-muon La différence entre la direction du muon au niveau du détecteur et la direction initiale du neutrino est gouvernée par deux processus : la cinématique de la réaction et les diffusions coulombiennes multiples subies par le muon avant d atteindre le détecteur. Au point d interaction, la direction du muon est d autant plus alignée sur celle du neutrino que son énergie est élevée : cos θ labo µν = β + cos θcm µν 1 + β cos θ cm µν β 1 1 (II.6) avec β = E ν E ν + M E ν + 1 avec θµν cm l angle dans le centre de masse entre le neutrino initial et le muon produit. Pour une énergie donnée du neutrino, les muons de plus haute énergie (θµν cm petits) sont favorisés par la détection. L angle moyen entre le neutrino incident et le muon induit au vertex et au détecteur varie selon l énergie (figure II.7), cette différence est due aux différentes diffusions subies par le muon au cours de sa propagation jusqu au détecteur.\n42 42 CHAPITRE II. TÉLESCOPES À NEUTRINOS FIGURE II.7 Angle moyen entre le neutrino et le muon incident au niveau du vertex et du détecteur en fonction de l énergie du neutrino. La différence entre les deux courbes est due au fait que le muon subit des diffusions lors de sa propagation jusqu au détecteur. 1.3.b Interactions dominantes Au cours de son trajet du vertex jusqu au détecteur le muon diffuse et perd de l énergie. En dessous de 1 TeV, le muon est au minimum d ionisation et perd de l énergie de manière continue par ionisation et excitation du milieu. A haute énergie, il perd son énergie via des phénomènes "catastrophiques" tels que la création de paire, le bremsstrahlung et les interactions photonucléaires (figure II.8). Ionisation : Lorsque l énergie cinétique est trop faible pour initier les processus radiatifs, les muons interagissent par diffusion élastique avec les électrons du cortège, leur cédant une petite quantité d énergie, néanmoins suffisante pour les délier du noyau. La matière, le long du parcours du muon, se retrouve ionisée. La quantité d énergie transmise par ionisation est caractérisée par la formule de Bethe-Bloch . Création de paires : Elle consiste en une production d une paire électron - positron par échange d un photon virtuel avec un noyau. C est le processus dominant dans la perte d énergie du muon au-delà de 1 TeV. La paire produit des cascades électromagnétiques que les particules secondaires (électrons et positons) continuent à alimenter. Bremsstrahlung : C est le rayonnement de freinage que subit le muon lorsqu il passe à proximité d un centre diffuseur, le muon émet alors un photon de recul. Interactions nucléaires : Il s agit de la contribution la plus modeste. C est une interaction rare sur les noyaux du milieu via un photon virtuel, mais lorsqu elle se produit, la perte d énergie est importante. Si bien qu aux alentours de 1 TeV, elle contribue à hauteur d environ % à la perte d énergie totale dans la roche. Perte globale : La contribution globale de tous ces processus se résume par la relation de Barrett  : de α + βe (II.7) dx\n43 1. PRINCIPE DE DÉTECTION 43 FIGURE II.8 Contribution à la perte d énergies des différents types de processus radiatifs engendrés par le passage d un muon dans la roche (à gauche) et dans l eau (à droite). avec α 2,2 MeV.g 1.cm 2 (terme d ionisation déterminé par la formule de Bethe-Bloch) et β 4 6 g 1.cm 2 (terme de radiation : production de paires, bremsstrahlung et contributions photonucléaire) . A basse énergie (< 1 TeV), le terme dominant provient de l ionisation : le parcours du muon est proportionnel à son énergie. Au-delà de 1 TeV environ, le second terme domine l équation. La perte d énergie moyenne est proportionnelle à l énergie. 1.3.c Effet Tcherenkov Le passage de particules chargées à travers la matière donne lieu à une émission de radiation électromagnétique appelée effet Tcherenkov. Dans les milieux transparents, il apparaît sous la forme de lumière dont les propriétés sont liées à la vitesse de la particule émettrice et au milieu. Les physiciens P.A. Tcherenkov  et S. Ivanovich ont énoncé les principales propriétés de cette émission en 1934 : l émission se produit dans tous les milieux transparents sous la forme d un spectre continu ; la lumière est fortement polarisée et émise dans une direction déterminée par rapport à la trajectoire de la source du rayonnement. Lors de son déplacement dans un milieu transparent d indice optique n, une particule chargée agit sur les milieux di-électriques pour former une série de sources émettant à des fréquences différentes. Ces sources sont décalées dans le temps de t 1 = L/v : temps que met la particule se déplaçant à la vitesse v pour aller d un point à un autre de sa trajectoire distants de L (figure II.9). Pour un observateur qui regarde la particule incidente sous un angle θ, le retard entre ces deux sources est différent à cause du temps de propagation des photons, et il observe : t = t 1 d/(c/n), avec d la différence de marche entre deux sources. Comme d = L cos θ, alors ( 1 t = L v n cos θ ) c\n44 44 CHAPITRE II. TÉLESCOPES À NEUTRINOS FIGURE II.9 Sillage conique de lumière Tcherenkov émis par une particule de vitesse v > c/n. Lorsque la particule se déplace à une vitesse supérieure à la vitesse de la lumière dans le milieu, on peut trouver un angle d observation pour lesquels les ondes sont constructives : 1 v n cos θ = 0 c Toutes les sources sont synchrones : il y a émission de lumière Tcherenkov. Dans tous les autres cas, les ondes sont destructives, il n y aura pas de détection possible de lumière. La condition d émission s écrit donc : cos θ = 1 (II.8) βn où β est le facteur de Lorentz. Dans l eau l angle Tcherenkov vaut environ 43. Les équations de Maxwell décrivent ce phénomène et permettent de calculer le nombre de photons dn γ émis par unité de distance dx : dn 600 γ dx = 2π dλ sin2 θ c 330 λ 2 si on prend en compte des longueurs d ondes dans le visible comprises entre 330 et 600 nm (gamme sensible des photo-multiplicateurs), on obtient environ 300 photons par cm. Le principe de reconstruction des traces de muons repose entièrement sur ce processus. Il s agit de détecter la lumière Tcherenkov émise par le muon à l aide d un ensemble de photo-multiplicateurs dont le positionnement est connu. Le temps d arrivée d un photon sur le photo-multiplicateur est aussi connu. On peut ainsi reconstruire la trajectoire du muon : soit montante, le muon vient de l autre côté de la Terre, soit descendante, le muon a été produit du même côté de la Terre que le détecteur. 1.4 Bruits de fond physiques La détection des neutrinos cosmiques est affectée par deux types de bruits de fond : les muons induits par les neutrinos créés dans les gerbes atmosphériques. Ces muons seront reconstruits montants dans le détecteur. Ce bruit de fond est irréductible. La composante des neutrinos atmosphériques est la même pour les deux hémisphères (figure II.). La détection de sources ponctuelles s effectue en calculant les excès de neutrinos cosmiques par rapport aux neutrinos atmosphériques dans les cartes du ciel. Le spectre en énergie des neutrinos atmosphériques étant plus dur que celui des neutrinos cosmiques, une coupure en énergie permet l élimination des neutrinos atmosphériques lors de l étude du flux diffus de neutrinos.\n45 1. PRINCIPE DE DÉTECTION 45 les muons produits dans les gerbes atmosphériques. Ces muons reconstruits descendants dominent de plusieurs ordres de grandeur le flux de muons induits par les neutrinos, il est donc nécessaire de se restreindre aux traces montantes. Le détecteur est optimisé pour détecter des muons issus de neutrinos, il est donc peu efficace lorsqu il s agit de détecter les muons atmosphériques. Il est aussi possible de réduire cette composante en augmentant la profondeur du site comme le montre la figure II.11. Le flux de muons atmosphériques et sa dépendance en fonction de la profondeur, qui sont étudiés dans cette thèse, sont connus et pemettent d approfondir la compréhension de notre détecteur. FIGURE II. Dépendance angulaire des différentes composantes du bruit de fond pour le site d ANTARES. cos θ = 1 représente les muons descendants verticaux. FIGURE II.11 Evolution du bruit de fond suivant la profondeur et l angle. Le flux de muons issus des neutrinos atmosphériques est obtenu avec la paramétrisation de Volkova  alors que le flux de muons atmosphériques est obtenu avec la paramétrisation d Okada \n46 46 CHAPITRE II. TÉLESCOPES À NEUTRINOS 2 Le détecteur ANTARES La collaboration ANTARES est composée de 20 instituts et d environ 150 chercheurs en physique des particules et des astroparticules, astronomes, océanologues et ingénieurs. Elle regroupe des instituts français, italiens, hollandais, allemands, espagnols, russes et roumains (figure II.12). FIGURE II.12 Disposition géographique des 20 instituts de la collaboration ANTARES. 2.1 Description du détecteur Le site sélectionné pour l implantation du détecteur ANTARES se situe à environ 40 km au large de La Seyne-sur-Mer (dans le Var). Ses coordonnées géographiques sont N, 6 E et sa profondeur est de m II.13. Le télescope à neutrinos ANTARES , schématisé figure II.14, est constitué de 12 lignes de détection indépendantes espacées les unes des autres de 60 à 75 m (figure II.16). Onze des douze lignes se composent de cinq secteurs de cinq étages identiques, espacés verticalement de 14,5 m, qui regroupent trois modules optiques (Optical Module -OM-) contenant des photo-multiplicateurs de pouces. La dernière ligne comporte 4 secteurs identiques aux secteurs des autres lignes et un secteur consacré à la détection acoustique. Chaque ligne comporte donc 75 modules optiques et la dernière ligne 60, soit un total de 885 pour le détecteur complet. Les lignes sont reliées à une boîte de jonction (figure II.15), installée en décembre 2002, qui est reliée à la côte via un câble électro-optique de 44,2 km de long. Elle assure l alimentation des lignes au moyen de transformateurs et contient aussi le système de distribution de l horloge vers les différentes lignes. Elle regroupe les fibres venant des lignes vers le câble principal qui va à la terre. Les connexions sous-marines entre les lignes et la boîte de jonction sont réalisées grâce à un submersible. Chaque étage est contrôlé par un module local de contrôle (Local Control Module -LCM-) qui, en plus d assurer les fonctions de distribution de tension et de transmission du signal, contient des éléments\n47 2. LE DÉTECTEUR ANTARES 47 FIGURE II.13 Carte géographique du site sélectionné pour ANTARES. permettant de déterminer l orientation de l étage : un compas et un inclinomètre. Le deuxième LCM de chaque secteur, en partant du bas, assure le contrôle du secteur entier : c est le module de contrôle maître (Master Local Control Module -MLCM-). 2.1.a Les lignes La ligne est ancrée sur le sol avec un pied de ligne lesté. Ce pied de ligne contient le module de contrôle de la ligne (String Control Module -SCM-), le module d alimentation de la ligne (String Power Module -SPM-), un émetteur-récepteur acoustique et l électronique de contrôle. Chaque étage de la ligne contient trois modules optiques orientés à 45 vers le bas et séparés entre eux de 120 afin d homogénéiser l acceptance angulaire pour les muons montants. Le premier étage est placé à 0 m au dessus de l ancre afin de laisser suffisamment de milieu transparent aux muons montants pour émettre de la lumière Tcherenkov. Les 3 modules optiques sont contrôlés par un LCM placé dans un cylindre au centre de la structure (figure II.17). Chaque secteur contient un MLCM qui contient également un switch Ethernet pour transmettre les signaux du secteur à la boîte de jonction. Chaque secteur contient également des hydrophones placés sur les étages 1, 8, 14, 19 et 25 de la ligne permettant un positionnement spatial en temps réel de la ligne. Ces hydrophones sont couplés à des pyramides accoustiques (figure II.15) pour le positionnement du bas des ligne. De plus, chaque ligne contient quatre balises LED aux étages 2, 9, 15 et 21. Les LED illuminent localement le détecteur pour contrôler l étalonnage temporel des modules optiques des lignes. 2.1.b Modules optiques Les modules optiques [46, 47] sont composés d une sphère en verre, résistante aux hautes pressions, d un large photo-multiplicateur demi-sphérique placé dans cette sphère à l aide d un gel optique. Une cage en µ-métal est utilisée comme barrière contre le champ magnétique terrestre, afin d amoindrir la déviation des photo-électrons dans le photo-multiplicateur. L électronique de lecture du photomultiplicateur n est pas dans la sphère mais dans le LCM. Les composants à l intérieur du module optique sont alors réduits au minimum : la base active Crawford qui transforme du 48 V en une tension de\n48 48 CHAPITRE II. TÉLESCOPES À NEUTRINOS FIGURE II.14 Schéma du détecteur ANTARES. FIGURE II.15 Répartition des lignes au niveau du sol. Les traits de couleurs représentent la position exacte des différents câbles reliant les lignes à la boîte de jonction. Les "Pyr" sont l emplacement des pyramides accoustiques utilisées lors du positionnement accoustique.\n49 2. LE DÉTECTEUR ANTARES 49 FIGURE II.16 Répartition idéale des lignes selon un octogone. FIGURE II.17 Photo d un étage avant son immersion, on peut y voir les trois modules optiques et le module de contrôle local. plus de 00 V et une LED qui est utilisée pour l étalonnage interne du photo-multiplicateur. La sphère contenant le photo-multiplicateur doit satisfaire le cahier des charges suivants : résistance à des pressions entre 260 (fonctionnement normal) et 700 bars (tests). transparence aux photons de longueurs d onde comprises entre 400 (bleu) et 500 (vert) nm. transmissibilité supérieure à 95%. Un trou dans la sphère permet de connecter le photo-multiplicateur au reste du détecteur. Afin d éviter tout risque de corrosion, le pénétrateur est en titane pour effectuer les connexions : l alimentation et la transmission du signal reçu. indice optique proche de l eau dans le but de minimiser la réflexion. La sphère NAUTILUS réunit toutes ces caractéristiques. La sphère est séparée en deux parties : la partie opposée à la photo-cathode du photo-multiplicateur est peinte en noire et la partie inférieure est en verre transparent. L axe du photo-multiplicateur par rapport à la verticale du détecteur est incliné vers le bas d un angle de 45, les photons issus des muons descendants sont moins visibles par le photomultiplicateur de par l orientation de la sphère. Le contact optique entre la sphère et le phototube est assuré pas un gel de silicone transparent (gel de Wacker), avec peu de bulles d air dans le gel, les bulles étant sources de diffusion. Afin de minimiser la réflexion, l indice optique du gel (1,404) est compris entre celui de la sphère (1,47) et de l eau (1,35). Ce gel permet aussi d assurer le maintien mécanique du photo-multiplicateur : il devrait être suffisament dur pour maintenir le PM et suffisament élastique pour supporter la déformation de la sphère sous la pression de l eau (diminution de 1,25 mm sur le diamètre de la sphère à 250 bars). La stabilité de ces différentes propriétés a été étudiée pour tenir dix ans. Le photon arrache un électron à la photo-cathode qui, par un phénomène de cascade sur les dynodes, se multiplie pour donner en sortie un signal électrique mesurable. La cage en µ-métal, alliage de fer et de\n50 50 CHAPITRE II. TÉLESCOPES À NEUTRINOS FIGURE II.18 Schéma du photo-multiplicateur. nickel, réduit la largeur du signal en sortie du photo-multiplicateur en éliminant la déviation des électrons par le champ magnétique terrestre grâce à la grande perméabilité magnétique de la cage. Le photo-multiplicateur utilisé pour ANTARES est un PM Hamamatsu de ". La résolution temporelle, inférieure à 3 ns, est caracterisée par la dispersion du temps de transit des électrons entre différents points de la photocathode et de l anode. Le PM étant composé de 14 dynodes, les hautes tensions ont été réglées de façon à obtenir un facteur d amplification de 5. 7 correspondant à un signal du photo-électron de 45 mv. Un pénétrateur fixé dans la sphère et un câble de 1,9 m de long assurent les connexions électriques entre les modules optiques et le contrôle local des modules (LCM). Le pénétrateur possède un diamètre de 20 mm, le câble contient 12 fils pour : le signal d anode le signal de la dynode 12 le signal de la dynode 14 l alimentation 48 V la commande de hautes tensions le contrôle de la haute tension le pulser de la LED interne Des mesures en laboratoire ont été effectuées afin de caractériser le comportement des modules optiques. Les résultats montrent que la réponse est maximum quand l angle entre l axe du module optique et la direction du muon est à θ = 43 (figure II.19), qui est l angle Tcherenkov dans l eau. L étude de l acceptance angulaire est détaillée en IV Caractéristiques optiques du milieu 2.2.a Longueur d atténuation dans l eau Un des paramètres qu il faut connaître est la longueur d atténuation de la lumière Tcherenkov dans l eau. Pour une source lumineuse ponctuelle isotrope d intensité I 0, l intensité I mesurée par un photomultiplicateur avec une aire effective A à une distance R est définie par : I(R) = I R 0A 4πR 2e λ eff att avec 1 λ eff att = 1 λ abs + 1 λ eff diff et λ eff diff λ diff 1 cosθ (II.9)\n51 2. LE DÉTECTEUR ANTARES 51 FIGURE II.19 Nombre de photo-électrons détectés en fonction de l angle formé par la direction de visée du module optique et la direction incidente du µ, pour des trajectoires éloignées d une distance inférieure à 50 cm de la photocathode. où λ abs est la longueur d absorption et λ diff est la longeur de diffusion. Le nombre de photo-électrons détectés par les photo-multiplicateurs devrait suivre le même comportement après correction de l acceptance angulaire. Q(t) Q R 0 R 2e λ eff att (II.) Des mesures sous-marines ont été effectuées avec les LED installées sur les lignes en flashant les modules optiques . En mesurant le nombre de photons en fonction de la distance, on trouve alors une longueur d atténuation d environ 50 m compatible avec les mesures précédentes (IV.4.1) 2.2.b Le bruit de fond optique La connaissance du bruit de fond optique sur le site ANTARES a été déterminante pour le choix de la logique de l électronique embarquée et de déclenchement. Ce bruit de fond est décomposé en deux parties : la ligne de base du bruit de fond est principalement due à la décroissance de 40 K dans l eau de mer et à des organismes microscopiques (figures II.20 et II.21). les sursauts de bruit de fond seraient dus à des passages d animaux macroscopiques (figure II.22). La désintégration β du 40 K intervient à la hauteur d environ 30 khz pour des PM de pouces à cause de la lumière Tcherenkov émise par l électron produit via : 40 19K Ca + e + ν e (II.11) La désintégration du 40 K de la sphère intervient à hauteur de quelques khz et produit des événements à haute charge. La composante variable de la ligne de base est due à une présence continue d organismes microscopiques vivant au fond de la mer. Cette composante apporte en général une contribution de 30kHz. La ligne de base des modules optiques est d environ 60 khz, mais le taux moyen vu par les photo-multiplicateurs depuis leur mise à l eau est de 166 khz . La fraction du temps où le taux de bruit de fond optique est 1,2 fois supérieur à la ligne de base du bruit de fond optique définit la fraction de sursauts. Les sursauts de bruit de fond sont fortement corrélés avec le courant (figure II.23). Cet effet pourrait être dû à une agitation d organismes dans des\n52 52 CHAPITRE II. TÉLESCOPES À NEUTRINOS turbulences formés en aval des modules optiques en cas de courant élevé où à des variations infimes de température ( 0,1 C). Cette agitation déclencherait une production de lumière chez ces organismes vivants provoquant des pics de luminosité dans le détecteur. 2.3 L électronique de lecture des photo-multiplicateurs. L électronique frontale des photo-multiplicateurs a pour but la sélection des signaux analogiques qui en sont issus au dessus d un certain seuil, et la numérisation du temps auquel se produit le signal ainsi que de la charge recueillie. Contrairement aux mesures de temps "classiques" par convertisseur tempsnumérique (Time to Digital Converter -TDC-), la numérisation n est pas déclenchée par une source extérieure mais par le signal lui-même, et le temps absolu auquel se produit le signal doit être enregistré. Ceci est rendu possible par un système d horloge global de l expérience, séquencée à 20 MHz et distribuée à tous les étages du détecteur. Le comptage des alternances de l horloge par rapport au début d une période de prise de données constitue le "TimeStamp" en unités de 50 ns. L électronique frontale doit donc fournir à la fois ce temps grossier et un temps fin indiquant à quel moment le signal a passé le seuil dans ces 50 ns. Le seuil choisi pour la détection de photons Tcherenkov, dit seuil de niveau 0 ou L0, est de l ordre du tiers de photo-électron, et résulte en des taux de comptages de l ordre de quelques dizaines de khz à quelques MHz. En raison de ces larges fluctuations de taux de comptage dues au bruit de fond optique de bioluminescence, il n est pas possible d envisager une transmission synchrone des informations issues de l électronique frontale jusqu à un processeur de lecture, et celle-ci dut être munie d une mémoire tampon permettant d absorber les fluctuations avant lecture des données. Dans le cas du circuit utilisé, il s agit de l ASIC "ARS" pour "Analog Ring Sampler" [50, 51] : cette mémoire est analogique et les informations sont stockées sous forme de tensions dans un "pipeline" analogique capacitif avant conversion du temps et de la charge. En ce qui concerne la charge, le signal du PM est intégré sur une durée de l ordre de 40 ns ; pour le temps fin, le passage du seuil mémorise la valeur de la tension d une rampe synchrone de l horloge de l expérience. La conversion par une suite de comparaisons séquencée localement à 25 MHz résulte en un codage de la charge sur 8 bits (appelé " AVC " pour Amplitude to Voltage Converter) et du temps fin sur 8 bits (appelé "TVC" pour Time to Voltage Converter). Les temps d intégration et de conversion résultent en un temps mort de l ordre de 250 ns ; pour diminuer ce temps mort à chaque PM sont associés 2 puces ARS interconnectées par une logique de " token ring " (jeton). Le temps mort de chaque PM n est alors plus que de 40 ns entre les 2 premiers signaux consécutifs, puis 250 ns entre le deuxième et le troisième signal, et ainsi de suite. La majeure partie des signaux des PM sont des signaux simples, de charge 1 photo-électron. Ils sont donc traités par l ASIC en mode dit "SPE" pour "Single Photo Electron", et l information de chaque signal ne consiste qu en ses TimeStamp, AVC et TVC codés sur 6 octets. Néanmoins l ASIC a également été muni d une fonctionnalité de gabarit dit "PSD" pour Pulse Shape Discriminator. Celui-ci permet d étiqueter les signaux de grande charge ou de durée supérieure à celle attendue pour des signaux simples et dont l intégration serait incomplète en 40 ns (voir figure II.26). Dans ce cas, l ARS peut éventuellement traiter le signal en mode oscilloscope dit "WFA" pour Wave Form Anode, avec un échantillonnage allant jusqu à 1 GHz. Les 128 échantillons convertis donnent alors lieu à un événement codé sur 272 octets. La fréquence de tels signaux est de l ordre de 1 %, et l acquisition d ANTARES est dimensionnée pour en permettre l enregistrement. Néanmoins il a été décidé que le mode d enregistrement des données des photo-multiplicateurs serait exclusivement en mode " SPE " pour les prises de données de physique, l activation du gabarit n étant effective que durant certaines prises de données de l étalonnage.\n53 2. LE DÉTECTEUR ANTARES 53 FIGURE II.20 Ligne de base du bruit de fond optique du mois de juin FIGURE II.21 Ligne de base du bruit de fond optique du mois d août La ligne de base est instable ce qui est dû à une grande activité des organismes microscopiques. FIGURE II.22 Sursauts répétés de bruit de fond optique du mois de juin La fraction du temps où le taux de bruit de fond optique est 1,2 fois supérieur à la ligne de base du bruit de fond optique définit la fraction de sursauts. FIGURE II.23 Evolution de la proportion de sursauts en fonction du courant pour trois étages du détecteurs.\n54 54 CHAPITRE II. TÉLESCOPES À NEUTRINOS FIGURE II.24 Représentation schématique des composantes du système d acquisition présentes dans le conteneur d éléctronique de chaque étage. L ASIC est configurable et son "contrôle lent" permet de configurer 75 paramètres, en particulier le seuil L0 et les paramètres de numérisation. La gamme dynamique choisie est de 20 photo-électrons, ce qui pour un codage sur 8 bits et en tenant compte d une marge pour un piédestal de charge non nul, résulte en environ une dizaine de coups d ADC par photo-électron. De la même façon, les 50 ns de temps fin ne sont codés que sur 200 canaux d ADC, résultant en une résolution de 250 ps (la distribution du TVC est plate en regroupant 4 bins). Ces points seront développés dans le chapitre "étalonnage". 2.4 Le système d acquisition des données Le système d acquisition et de déclenchement d ANTARES  est basé sur un schéma classique pour les expériences de physique des hautes énergies. Il doit en principe accomplir les tâches d envoi des données correspondant à une condition de déclenchement locale (données dites de niveau 1), de rassembler les données des différents étages du détecteur appartenant à un même événement pour l envoyer à un système de déclenchement centralisé dit de niveau 2 dont le taux doit être de quelques Hz de façon à permettre l enregistrement sur disque des données dans la "région d intérêt" du déclenchement qui va constituer un événement. Néanmoins, la collaboration ANTARES a décidé d éviter un déclenchement local figé, et de fait toutes les données de niveau 0 (signaux de photo-multiplicateurs "L0") doivent être transférées au système de déclenchement centralisé à terre, système qui doit alors jouer les rôles de déclenchement de niveau 1, de niveau 2, voire de niveau 3 si on effectue une reconstruction en ligne. Ce concept est appelé "all data to shore" puisque toutes les données de niveau 0 sont acheminées à terre. En raison du large volume de données, la combinatoire de l algorithme de déclenchement ne peut être effectuée sur un simple processeur, et les données de chaque photo-multiplicateur sont regroupées en tranches temporelles appelées "TimeSlice" en fonction de leur temps d arrivée. Une tranche temporelle complète, comprenant toutes les données des 12 lignes, est alors acheminée vers un processeur d une ferme de PC à terre. L architecture qui en résulte se compose d une collection de 300 étages "envoyeurs" connectés à une collection d une cinquantaine de "cibles", où une cible est le processeur d un PC à terre, par l intermédiaire d un "aiguilleur ethernet" ("switch") totalement croisé. En d autres termes, il s agit d un réseau local Gigabit de quelques centaines de noeuds. Le dimensionnement du système de lecture est alors fonction du taux de données devant être en-\n55 2. LE DÉTECTEUR ANTARES 55 FIGURE II.25 Représentation schématique de l équipement (hardware) du système d acquisition. Une seule ligne de détection est représentée. A chaque étage, trois PMs et un module (storey module) pour l électronique associée sont représentés. Ce module contient une horloge locale et un processeur (CPU). Un module sur cinq (sector module) contient également une carte éthernet et un émetteur-récepteur DWDM. En bas de chaque ligne se trouve un module (string module) dans lequel on trouve également un démultiplexeur optique (DEMUX). Chaque ligne est connectée à la boîte de jonction principale par une câble interlink (ILC). Les données sont ensuite achminées à terre via le câble électro-optique reliant la boîte de jonction et la station d acquisistion à terre qui abritte la ferme de station de calculs et l horloge de référence principale (master clock).\n56 56 CHAPITRE II. TÉLESCOPES À NEUTRINOS FIGURE II.26 fonctionnement de la discrimination SPE ou WF. voyées à terre. Chaque étage est muni d une carte d acquisition spécifique comprenant un processeur, une interface ethernet 0 Mb optique, une RAM de 64 Moctets, un FPGA (Field Programmable Gate Array) et un bus reliant ces composants. Le FPGA est programmé comme machine esclave transférant les données des ARS à la RAM, et les formatant en tranches temporelles (figure II.24). Le processeur envoie de façon asynchrone les tranches temporelles de la RAM jusqu au PC associé à terre par le lien ethernet 0 Mb. Le séquençage de la lecture des ARS à 25 MHz limite le taux de données enregistrables de chaque PM à 1 MHz. La carte d acquisition elle-même ne peut transmettre en continu plus de 50 Mb/s, ce qui limite le taux moyen enregistrable de chaque PM à 350 khz. Ces valeurs sont toutefois bien au delà des taux nominaux qui sont de l ordre de 60 khz, et ne constituent pas une contrainte pour les prises de données. Afin de limiter le nombre de fibres optiques, 5 étages consécutifs d une ligne sont groupé en un secteur, dont l étage maître contient un concentrateur ethernet (" switch ") et une carte dite " DWDM " (Dense Wavelength Division Multiplexing) émettant cette fois en ethernet Gb avec un laser puissant pour aller jusqu à la terre sur une longueur d onde associée au numéro de secteur (figure II.25). Les 5 longueurs d ondes des 5 secteurs d une ligne sont multiplexées dans le module SCM de bas de ligne sur une fibre unique, puis démultiplexées à terre avant d être connectées à l aiguilleur ethernet. Du point de vue de l acquisition, la boite de jonction ne fait que regrouper les 24 fibres (12 aller et 12 retour) des différentes lignes dans le câble électro-optique qui relie la côte distante de 40 km. Par contre elle distribue les deux fibres (une nominale et une de rechange) du système d horloge sur les 16 connecteurs dont elle est munie. A terre, les processeurs de déclenchement sont des PC commerciaux (DELL PowerEdge 950) montés dans des racks situés dans une pièce spécifique consacrée à l Institut Michel Pacha. Ils exécutent différents algorithmes pour faire les coïncidences locales de niveau 1, puis recherchent des accumulations de ces coïncidences liées par relation de causalité, et compatibles avec la lumière Tcherenkov. Le taux de déclenchement de niveau 2 est de quelques Hz pour le détecteur à 5 lignes, et une dizaine de PC suffisent à traiter le volume de données. En cas de période agitée (augmentation de la bioluminescence), il peut y avoir engorgement du système. Dans ce cas la consommation de données se ralentit et les " envoyeurs " du détecteur ne peuvent vider les données aussi vite qu elles se produisent. Les RAM de la carte d acquisition se remplissent jusqu à un niveau limite où le processeur de la carte cesse l enregistrement des données des ARS. Dans ce cas, les tranches temporelles de données seront vides, mais étiquetées " XOFF " et constituent un temps mort macroscopique de l expérience. On limite la survenue de cet engorgement par un mécanisme de veto sur les taux élevés (" High rate Veto "), dans lequel le processeur de la carte\n57 2. LE DÉTECTEUR ANTARES 57 décide par lui-même l arrêt de l enregistrement des données dès qu un taux important est détecté durant une tranche temporelle. La valeur de ce taux est généralement 400 khz, et la condition est remplie pour tout sursaut de bioluminescence. L optimisation de l acquisition a conduit à adopter des tranches temporelles de ms, soit ns, permettant l envoi par TCP/IP de " fichiers " de taille convenable pour le système d exploitation VxWorks embarqué. 2.5 Algorithmes de filtrage Les données sont transférées à terre mais elles sont filtrées par différents types d algorithmes avant d être sauvegardées sur disques. Seules les données présentant les caractéristiques d un événement physique sont écrites sur disque. Les différents types d algorithme de filtrage sont : algorithme de filtrage 3N : le premier niveau de filtrage, L1, est composé de deux coups en coïncidence sur le même étage dans une fenêtre de 20 ns, ou de coups de grande charge (entre 2,5 et photo-électrons). Lorsque cinq L1 sont repérés en coïncidence dans une fenêtre de temps compatible avec la vitesse du muon, les conditions de l algorithme de filtrage sont satisfaites et les données sont sauvegardées. algorithme de filtrage T3 : l algorithme de filtrage T3 est défini par la coincidence de deux L1 dans une fenêtre de 80 ns entre deux étages adjacents ou une fenêtre de 160 ns entre deux étages distant d un étage. algorithme de filtrage du centre galactique : cet algorithme de filtrage demande 1 L1 et 4 L0 alignés selon une direction. Cet algorithme de filtrage suit le centre galactique. On observe une variation du nombre d événements acceptés par l algorithme galactique due à un effet de variation géométrique du détecteur à 5 lignes (figure II.27). algorithme de filtrage à "biais minimum" : afin d étudier le bruit de fond, on enregistre sans filtrage une tranche temporelle de 4 µs une fois par seconde. Prise de données spéciale 40 K : ces prises de données de calibration consistent à enregistrer des coïncidences de 2 L0 sur 2 modules optiques d un même étage dans une fenêtre de 50 ns. Elles servent à vérifier localement la stabilité de l étalonnage en temps, les gains et les seuils des PMs. En plus, des paires de L1 en coïncidence dans une fenêtre de temps compatible avec le temps de parcours d un muon permet une mesure simple du flux de muons. Les différents algorithmes de filtrage ont des caractéristiques différentes présentées à la figure II.27. Les données utilisées pour les études présentées aux chapitres suivants ont été obtenus à partir du filtrage 3N. Il n est pas possible de déterminer l efficacité de déclenchement à partir des données Minimum-bias car celles-ci sont trop peu nombreuses. On la déterminera donc par Monte Carlo en faisant subir les algorithmes de filtrage aux données simulées. Dans le cadre d analyse multi-messagers, ANTARES est connecté à un réseau d alertes (Global Coordinate Network -GCN). Lors des alertes SWIFT ou FERMI, la totalité des données sur une période de deux minutes centrée autour du GRB est sauvegardée sur disque. Nous nous servons de ces données dans le cadre de l étalonnage lorsque des données brutes sont nécessaires. 2.6 Reconstruction des traces de muons Il existe différentes façons d obtenir la direction du muon. La reconstruction de Aart  utilise l algorithme de filtrage 3N couplée à un maximum de vraisemblance. La reconstruction BBFit  utilise l algorithme de filtrage T3 et minimise les paramètres de la trace avec un χ 2.\n58 58 CHAPITRE II. TÉLESCOPES À NEUTRINOS FIGURE II.27 Taux des différents algorithmes de filtrage appliqués dans ANTARES. Données du 21 juillet a Principe de reconstruction La trajectoire du muon est caractérisée par sa direction d (d x,d y,d z ) et par sa position p (x 0,y 0,z 0 ) à un temps donné t 0, choisi arbitrairement. La direction s écrit en coordonnées sphériques d (sin θ cos φ,sin θ sin φ,cos θ) où θ et φ sont le zénith et l azimuth. Les paramètres à déterminer pour reconstruire une trace sont donc au nombre de cinq : x 0, y 0, z 0, θ, φ. Si le muon se propage en ligne droite à la vitesse c, le temps d arrivée du photon émis par lumière Tcherenkov direct est déterminé par : t théorique = t c ( l k tan θ C ) + 1 ( ) k v g sin θ C (II.12) avec v g la vitesse de groupe de la lumière à une longueur d onde de 460 nm [55, 56, 57], t 0 le temps de l événement, θ C l angle Tcherenkov. Les paramètres l et k sont définis figure II.28. Le second terme de cette équation est le temps nécessaire au muon pour atteindre le point d émission du photon Tcherenkov direct détecté. Le troisième terme est le temps mis par ce photon pour atteindre le photo-multiplicateur. FIGURE II.28 Illustration du calcul du temps d arrivée du photon sur le module optique. La formule II.12 est valable dans le cadre d une reconstrution tri-dimensionnelle. Si une seule ligne\n59 2. LE DÉTECTEUR ANTARES 59 est touchée, les faibles déviations de la ligne avec le courant et la faible distance entre deux modules optiques d un même étage impliquent que le problème se simplifie par symétrie autour de l axe de la ligne et la dégénerescence de l angle azimuthal s en suit. L équation II.12 s écrit alors : c(t j t 0 ) = (z j z 0 )cos θ + n 2 1 d 2 + (z j z 0 )sin 2 θ (II.13) 2.6.b Reconstruction des gerbes Lorsqu une gerbe électromagnétique ou hadronique est produite, la gerbe se développe sur une courte distance relativement aux distances caractéristiques du détecteur. La lumière est émise en un point (x 0,y 0,z 0 ) à un temps t 0 puis se propage à la vitesse v g de façon isotrope. L enveloppe sphérique de la lumière peut être reconstruite avec un ajustement à quatre paramètres (x 0,y 0,z 0,t 0 ), le temps t j d arrivée en (x j,y j,z j ) vérifiant alors : v g (t j t 0 ) = (x j x 0 ) 2 + (y j y 0 ) 2 + (z j z 0 ) 2 (II.14) Si une seule ligne est touchée, on reconstruit avec un ajustement à 3 paramètres (d,z 0,t 0 ), d étant la distance séparant le point de l émission de la lumière à la ligne, d où : v g (t j t 0 ) = d 2 + (z j z 0 ) 2 (II.15) 2.6.c Description des reconstructions Dans le cadre de l étude du chapitre VI, nous allons utiliser deux types de reconstruction pour les muons et un pour les gerbes. La première est basée sur une fonction de densité de probabilité et l autre sur une étude de χ 2 avec sélection de coups. La reconstruction de Aart, du renom de son développeur Aart Heijboer, utilise la géométrie réelle du détecteur (position des modules optiques) alors que la reconstruction BBfit utilise une géométrie idéalisée du détecteur (elle a été conçue pour une utilisation en ligne). Dans les deux cas, la reconstruction se fonde sur l hypothèse d un muon unique. La reconstruction de Aart est optimisée pour reconstruire des muons issus de neutrinos alors que BBfit devrait être valable pour la reconstruction des muons descendants et montants en étant moins efficace que Aart pour les muons montants. Reconstruction de Aart Elle s effectue en 4 étapes décrites ci-dessous. Tout d abord une pré-sélection des coups permet de ne garder que les coups qui proviennent d un muon : T (d/v g + 0ns) avec T l écart temporel entre l impulsion considérée et l impulsion de plus grande amplitude et d la distance entre les photo-multiplicateurs de ces deux impulsions. 1. Régression linéaire : On effectue une régression linéaire qui passe par la position des coups sur les modules optiques. Les coups à haute charge (> 3 photo-électrons) seront plus près de la trace du muon que les coups à basse charge, diminuant l erreur sur la régression linéaire. Ainsi cette régression ne prend en compte que les coups de grande charge et les coups en coïncidence. Une coïncidence est définie par la coïncidence de 2 impulsions (ou plus) sur un même étage en 25 ns. 2. Le M-estimateur : Afin d obtenir une estimation des paramètres de la trace, on minimise la fonction suivante : ( ) G = i α A ir 2 i 2 (1 α)f ang (a i ) (II.16)\n60 60 CHAPITRE II. TÉLESCOPES À NEUTRINOS ( p avec A i la charge du coup i en photo-électrons, r i = t i t th i, d) le résidu temporel en ns (t th i le temps attendu d arrivée du photon), a i le cosinus de l angle entre le photon et l axe du photo-multiplicateur et α = 0,05 est issu de la simulation. Pour calculer a i, on suppose que le photon est émis par rapport à la trace du muon avec un angle Tcherenkov. La fonction f ang décrit l acceptance angulaire du photo-multiplicateur. Le M-estimateur est un χ 2 auquel on associe des poids faibles pour les coups éloignés du préfit. 3. Maximum de vraisemblance : Pour chaque possibilité de combinaison de paramètres, il est possible de calculer une probabilité : la vraisemblance. On utilise une fonction de densité de probabilité obtenue à partir de la mesure du temps de propagation de la lumière Tcherenkov dans l eau. Comme les probabilités des coups sont non corrélées, la vraisemblance s exprime comme le produit des vraisemblances des coups individuels : P (événement trace) P ( coups p, ) d = i ( P t i t th i ) (II.17) avec t i le temps du coup et t th i le temps attendu d arrivée du photon. La combinaison de paramètres acceptée est celle qui maximise cette vraisemblance. La fonction de densité de probabilité des résidus temporels (figure II.29) a été paramétrisée par  : dp signal dr = Ae (r τ)2 2σ 2 si r < c 1 B(αr 3 + βr 2 + γr + 1) si c 1 < r < c 2 (II.18) C e r/η r+ρ si r > c 2 les variables c 1 et c 2 étant des variables libres de l ajustement. Le pic de la distribution est ajusté par une fonction gaussienne alors que la queue est ajustée par la fonction e r/η r+ρ. Les deux fonctions sont reliées par un polynôme du troisième degré dont les coefficients sont fixés pour obtenir la continuité de la fonction de densité de probabilité. Le pic à 25 ns de la figure II.29 représente le temps d intégration de l ARS et le passage d un ARS à un autre, cette valeur vaut désormais 40 ns dans la simulation. Cet effet n a pas encore été incorporé à la fonction de densité de probabilité. 4. les deux étapes 2 et 3 sont répétées 9 fois en effectuant des rotations et des translations à partir de la trace de départ. Le meilleur résultat est alors conservé. On effectue alors une dernière analyse de maximum de vraisemblance avec une fonction de densité de probabilité améliorée dans une fenêtre de résidus temporels de ±250ns. La figure II.30 présente l erreur de reconstruction à chaque étape. On passe de 1,1% des traces reconstruites avec une erreur inférieure à 1 à 38% après le M-estimateur, puis à 57% après l ajustement de l étape 3. Au final on obtient 59% des traces reconstruites avec une erreur inférieure à 1. Le second pic de α est un artefact de l espace des phases. La figure II.31 présente le nombre d événements restants après la coupure de qualité. Les muons reconstruits montants sont en jaune et la totalité des muons est représentée en vert. Il est alors possible de choisir une coupure qui minimise le nombre de muons montants et maximise le nombre de muons descendants. Reconstruction BBFit La reconstruction BBFit  commence par un filtrage identique à celui de l algorithme de filtrage T3. Les déviations de la ligne ou l orientation des photo-multiplicateurs ne sont pas prises en compte. Les positions des coups sont placées au centre de l étage. Dans cette reconstruction les événements sont séparés en quatre catégories : les muons reconstruits à une ligne, les muons à plusieurs lignes, les gerbes reconstruites à une ligne et les gerbes à plusieurs lignes.\n61 2. LE DÉTECTEUR ANTARES 61 FIGURE II.29 Fonction de densité de probabilité pour les résidus temporels pour différentes charges A en photo-électron et distribution de ces résidus. Dans l encart, la même distribution est représentée en échelle linéaire. \n62 62 CHAPITRE II. TÉLESCOPES À NEUTRINOS FIGURE II.30 Performance des différentes étapes de la reconstruction. FIGURE II.31 Nombre d événements restants en fonction de la coupure de qualité (histogramme intégral). Il reste les neutrinos atmosphériques et les muons atmosphériques passant le critère. L efficacité du signal est aussi indiquée. La première étape est une régression linéaire, dont les paramètres sont utilisés comme point de départ des ajustements suivants. La seconde étape est une minimisation du χ 2 : χ 2 = N coup i=1 [ 1 σ 2 (t fit t i ) 2 + q(q ] i,q 0,α)d(d fit,d 0 ) q d 1 (II.19) Le premier terme de l équation correspond aux résidus temporels. Les quantités q 0 et d 0 sont les valeurs de charge et de distance à saturation, α une fonction angulaire dépendante de la position de l étage, q la charge moyenne de l événement, d 1 = 50 m et d fit la distance par un ajustement en suivant les formules II.13, II.14 et II.15. Le second terme prend en compte la perte d énergie des muons par ionisation. Afin d éviter une coupure à haute énergie, la fonction de charge est définie par : q i q 0 q = (II.20) qi 2 + q2 0 avec q 0 = photo-électrons introduisant une saturation artificielle des ARS. La fonction de dépendance en distance est définie par : d = d d2 fit (II.21) avec d 0 = 5 m afin d éviter des biais de reconstruction. Le produit d 0.q 0 est constant à 50 m.photoélectrons, car un photo-multiplicateur mesure un signal de 50 photo-électrons issus de la lumière Tcherenkov émise à 1 mètre par un muon dont la direction est égale à 45 (confere acceptance). Le critère de qualité est le rapport du χ 2 (équation II.19) par le nombre de degrés de liberté. Cette variable permet de déterminer la qualité de l ajustement. Ce critère présenté à la figure II.32 permet de choisir facilement une valeur de discrimination entre les muons et les neutrinos.\n63 3. LES AUTRES TÉLESCOPES À NEUTRINOS 63 FIGURE II.32 Coupure de qualité de BBFit pour les muons et pour les neutrinos pour un détecteur à cinq lignes. 3 Les autres télescopes à neutrinos 3.1 Un télescope en eau douce : Baïkal Le premier détecteur à neutrino de haute énergie qui fut opérationnel est situé en Sibérie dans le lac Baïkal. Ce télescope est immergé à 10 m de profondeur, à 3,6 km de la côte. Les opérations de déploiement et de maintenance s effectuent en hiver, profitant de l épaisse couche de glace qui recouvre tout le lac. Après presque dix ans d études préliminaires et cinq ans de construction, le télescope Baïkal fonctionne depuis 1998 avec 192 modules constitués de photo-multiplicateurs de 15 pouces répartis sur 8 lignes (configuration NT-200 figure II.33). L ajout de trois lignes additionnelles (figure II.34) est un premier pas vers un détecteur de taille kilométrique dans ce lac . Au cours de son fonctionnement, les mesures effectuées ont permis d estimer le flux de neutrinos atmosphériques concordant avec les prédictions théoriques. De plus, elles ont abouti à la détermination d une limite sur le flux diffus de neutrinos cosmiques ainsi que sur le flux de neutrinos produit par annihilation de neutralinos au centre de la terre . 3.2 Les télescopes sous la glace 3.2.a AMANDA Une collaboration de physiciens américains, allemands, belges et suédois lancèrent le projet AMANDA à la fin des années 80 (Antartic Muon And Neutrino Detector Array). Ce télescope se situe au pôle Sud, il utilise l épaisse couche de glace de 3 km comme milieu transparent. L avantage de ce choix provient du faible bruit de fond optique, de l ordre de 1 khz par photo-multiplicateur de 8 pouces, et de la facilité de déploiement. Cependant l utilisation de la glace au lieu de l eau liquide, implique une forte diffusion de la lumière due à la présence dans la glace de micros bulles d air et de poussières. La glace est creusée\n64 64 CHAPITRE II. TÉLESCOPES À NEUTRINOS FIGURE NT200. II.33 Schéma du détecteur Baïkal FIGURE II.34 Schéma du détecteur Baïkal NT200+. à l aide de jets d eau chaude. Les modules optiques sont protégés de la pression exercée par la glace à l aide de sphères en verre. AMANDA-II est constitué de 676 modules optiques. Ce télescope est complété en surface par un détecteur de gerbes atmosphériques SPASE (South Pole Air Shower Experiment) constitué d un réseau de 30 stations séparées de 30 m et contenant chacune 4 scintillateurs de 0,2 m 2 et un télescope à imagerie gamma GASP (Gamma Astronomy South Pole). Ces deux détecteurs permettent de vérifier les mesures de flux de muons atmosphériques descendants et la résolution angulaire du télescope (de l ordre de 3 ). La collaboration a pu déterminer le flux de muons en fonction de l angle zénithal et en foncion de la profondeur (figure II.35) . Actuellement, la limite la plus contraignante sur le flux de neutrinos de haute énergie a été obtenue par cette collaboration. Toutefois la mauvaise résolution angulaire de ce détecteur et le faible nombre d événements attendus rendent la recherche de sources ponctuelles très difficile. C est pourquoi la collaboration a lancé la construction de la seconde génération de télescope, de taille kilométrique : IceCube. 3.2.b IceCube Comme expliqué ci-dessus, IceCube est la continuation d AMANDA-II à l échelle kilométrique (figure II.36). La première ligne fut installée en janvier 2005 et le détecteur est actuellement composé de 59 lignes. Le détecteur sera finalisé en 2011 et composé de 80 lignes de 60 modules optiques, pour un total de photo-multiplicateurs de pouces, figés dans la glace du Pôle Sud entre m et m de profondeur. À la surface, un détecteur, IceTop, sert de système d étalonnage pour IceCube. IceTop sera constitué de 80 stations placées près du haut des lignes IceCube, chaque station étant composée de deux cuves Tcherenkov. Les lignes et les stations sont déployées au même rythme, il y a\n65 3. LES AUTRES TÉLESCOPES À NEUTRINOS 65 FIGURE II.35 Flux de muons en fonction de l angle zénithal (A) et de la profondeur (B) mesuré par AMANDA-II . FIGURE II.36 Schéma d IceCube comprenant AMANDA-II en bleu foncé. IceTop est installé sur la glace.\n66 66 CHAPITRE II. TÉLESCOPES À NEUTRINOS donc actuellement 59 stations . IceTop étudie la distribution angulaire des gerbes et leur énergie afin de contraindre les modèles sur les flux des rayons cosmiques et leur propagation dans l atmosphère. La taille du détecteur IceCube permettra de détecter des traces de muons issus des ν µ caractérisées par de longues traces. La reconstruction permettra d obtenir une résolution angulaire de l ordre du degré et une résolution en énergie faible (facteur2 soit d(log )E = 0.3). des cascades reconstruites avec une bonne résolution sur l énergie (%) mais la directionnalité n est pas disponible. des ν τ par des événements "double bang", qui sont des événements contenus ou des "lollipop" semi-contenus. Lors de la prise de données du détecteur à 22 lignes (IC22), 8727 candidats neutrinos ont été sélectionnés en 276 jours. Des recherches de sources ponctuelles ont alors été effectuées  (figure II.37). L accumulation la plus élevée révèle une déviation de 2,2σ seulement par rapport au bruit de fond des neutrinos atmosphériques. IceCube est le plus gros détecteur de neutrinos au monde. Il observe le ciel de l hémisphère Nord. Une extension d IceCube est en cours de construction : IceCube Deep Core. Il s agit de rajouter 6 lignes dans le réseau des lignes IceCube afin d obtenir un réseau plus dense de lignes dans une sous-partie du détecteur. Abaissant ainsi le seuil en énergie, les études de matière noire (WIMPs) et d oscillations de neutrinos seront favorisées. Une étude est en cours pour détecter les neutrinos descendants (du ciel de l hémisphère Sud) en utilisant IceCube comme véto, ce qui permettrait de regarder l ensemble du ciel. 3.3 Les télescopes en eau salée 3.3.a Dumand Le premier essai de construction d un télescope à neutrinos sous-marin a été effectué par la collaboration Dumand (Deep Underwater Muon And Neutrino Detector) au large de la côte hawaienne . En 1987, la collaboration a déployé une ligne prototype attaché à un bateau et a effectué quelques mesures. A la suite de ce succès, un projet de construction d un réseau de 9 lignes ancrées au font de la mer à 4800 m de profondeur a été soumis. Seule la première ligne a été déployée. Elle a été opérationnelle pendant quelques heures. Une fuite s est alors produite dans l une des unités électriques. L expérience a finalement été définitivement arrêtée en 1996, par manque de financement, et après 18 années de travaux. 3.3.b Nemo La collaboration italienne Nemo (NEutrino Mediterranean Observatory), créée en 1999, projette la construction d un télescope de taille kilométrique au large de Capo Passero en Sicile, à 3350 m de profondeur. L architecture proposée est constituée de photo-multiplicateurs fixés sur des bras rigides et reliés par des câbles pour former des tours de 750 mètres. Chaque bras se compose de quatre modules optiques, deux orientés horizontalemen et deux orientés vers le bas. Ce télescope formé de 64 tours espacées les unes des autres de 200 mètres soit au total modules optiques, devrait obtenir une résolution angulaire meilleure que 0,3. Une mini-tour, composée de 4 étages, a été déployée à 2000 m de profondeur et exploitée entre le 18 décembre 2006 et le 18 Mai Chaque étage comporte une structure de 15 m comportant 2 modules optiques à chaque bout, un regardant horizontalement et l autre verticalement (figure II.38). Les étages sont séparés de 40 m, avec leur structure perpendiculaire les unes aux autres . Cette mini-tour constitue la phase 1 de Nemo, une deuxième phase est prévu avec le déploiement d une tour complète.\n67 3. LES AUTRES TÉLESCOPES À NEUTRINOS 67 FIGURE II.37 Recherches de sources ponctuelles dans IceCube avec le détecteur à 22 lignes et une prise de données de 250 jours. La résolution angulaire pour IC22 est de 1,5. La déviation la plus grande obtenue est une déviation de 2,2σ. FIGURE II.39 Schéma d un étage de la mini-tour Nemo. FIGURE II.38 Schéma de la mini tour du détecteur Nemo.\n68 68 CHAPITRE II. TÉLESCOPES À NEUTRINOS 3.3.c Nestor NESTOR (Neutrino Experimental Submarine Telescope with Oceanography Research) a pour but le déploiement d un détecteur de taille kilométrique au large de la baie de Navarino en Grèce, à près de mètres de profondeur. Ce détecteur est constitué d une tour semi-rigides autonomes, de 12 étages espacés de 30 m (figure II.40). Chaque étage comporte 6 bras de 16 m, disposés radialement, dont les extrémités sont équipées de photo-multiplicateurs, dont la moitié est orientée vers le haut et l autre vers le bas. Le détecteur final contiendra en tout modules optiques. Le câble principal de 28 km reliant le futur détecteur à la côte a été immergé en Puis en 2003, la collaboration a déployé un étage constitué de 12 photo-multiplicateurs de 15 pouces aux extrémités des bras d une étoile de 6 m de rayon, pour valider les choix technologiques. Les mesures effectuées avec une ligne test ont permis de mesurer le flux de muons atmosphériques à une profondeur de l ordre de 3500 m . FIGURE II.40 Représentation schématique de la tour Nestor 3.3.d KM3Net Si IceCube constituera le premier télescope de neutrinos de haute énergie de taille kilométrique au monde, il ne pourra, de par sa position au pôle Sud, observer l intégralité du ciel, en particulier le centre\n69 3. LES AUTRES TÉLESCOPES À NEUTRINOS 69 de notre galaxie. Il est donc impératif de construire un détecteur équivalent dans l hémisphère Nord dont la sensibilité aux sources sera meilleure que les détecteurs ANTARES ou Baïkal. Dans cette optique, des instituts de pays européens (Allemagne, Angleterre, Chypre, Espagne, France, Grèce, Irlande, Italie, Pays- Bas, Roumanie) se sont regroupés dans le consortium KM3NeT en Ce projet a pour ambition d unifier les efforts des différents instituts (généralement impliqués dans un des projets méditerranéens : ANTARES, NEMO, NESTOR) pour définir les caractéristiques d un futur télescope méditerranéen de taille kilométrique. Plusieurs configurations peuvent être envisagées (figure II.41) pour ce détecteur : une configuration homogène, une configuration en grappes ou une configuration en anneaux par exemple  ; ces différentes configurations ne sont pas efficaces pour les mêmes énergies et pour les même directions. Sont envisagés différents types de modules optiques (un seul gros ou plusieurs petits dans une même sphère), différents types d électronique, différentes formes de lignes, etc. Avec la configuration suivante : configuration homogène de taille lignes, distance de 95 m entre chaque ligne, 37 étages et une distance entre étages de 15,5m, un étage est composé d un module optique incluant 23 photo-multiplicateurs de 3. on peut alors faire une étude de sensibilité aux sources (figure II.42) et une étude sur le flux diffus de neutrinos muoniques (figure II.43). Dans ces études, le bruit de fond est composé de neutrinos atmosphériques uniquement. La sensibilité aux sources est alors e flux nécessaire pour qu au bout d un an de données on ait détecté à 90% de chances une source. Le détecteur de taille kilométrique en mer Méditerrannée sera deux fois plus sensible aux sources ponctuelles qu Icecube. La résolution angulaire de KM3Net sera de l ordre de 0,1 au-dessus de 30 TeV contre 0,6 au-dessus de 0 Tev dans IceCube ce qui rend KM3Net compétitif dans la recherche de sources ponctuelles.\n70 70 CHAPITRE II. TÉLESCOPES À NEUTRINOS FIGURE II.41 Vue du dessus (haut) et vue tridimensionnelle (bas) de plusieurs configurations possibles pour un détecteur de taille kilométrique. Les points sur les figures du haut représentent la position des lignes alors qu en bas ils représentent la position d un étage. FIGURE II.42 Sensibilité du détecteur KM3Net décrit comparé aux limites actuelles de MACRO and AMANDA-II et aux limites attendues d AN- TARES et IceCube. Les points indiquent les limites pour des sources spécifiques.  FIGURE II.43 Limite sur le flux diffux du détecteur KM3Net décrit comparé à différentes expériences pour différents modèles de flux. \n71 Chapitre III Etalonnage des lignes du détecteur Antares 71\n72 72 CHAPITRE III. ETALONNAGE DES LIGNES DU DÉTECTEUR ANTARES 1 Méthodes d étalonnage L étalonnage à terre des lignes ANTARES permet de vérifier le fonctionnement de chaque partie du détecteur et d obtenir une première détermination des paramètres d étalonnage de l ARS autant en charge qu en temps. Nous pouvons ainsi remonter des quantités numériques de l AVC, du TVC et du timestamp aux quantités physiques que sont la charge en nombre de photo-électrons et le temps en ns. Le détecteur est calibré d abord à terre en mettant les modules optiques dans le noir puis in situ à 2475 m de profondeur. La moitié des lignes du détecteur ANTARES a été intégrée et calibrée au CEA Saclay et l autre moitié au Centre de Physiques des Particules de Marseille (CPPM). Ici ne sera présentée que l intégration effectuée au CEA Saclay (figure III.1). Les méthodes d étalonnage développées au CEA sont similaires à celle du CPPM à cela près que l étalonnage se fait en chambre noire au CPPM alors qu elle se fait en boîtes noires au CEA. FIGURE III.1 A gauche, hall d intégration de Saclay. A droite, trois modules optiques dans une boîte noire. Le système de refroidissement des modules de contrôle local (Local Control Module -LCM-) permet d obtenir une température proche de celle in situ ( 28 C). Les valeurs de TVCs sont susceptibles de varier de façon importante avec la température, ce système permet de garder des valeurs constantes entre l étalonnage à terre et in situ. L intégration des lignes s effectue secteur par secteur (5 étages). Les tests des modules optiques du secteur sont effectués selon le schéma de la figure III.3 : Les trois modules optiques associés à leur module de contrôle local (Local Control Module - LCM-) sont placés dans une seule et même boîte noire. La partie sensible est orientée vers le bas et placée dans des pots de fleurs entourés de tissus noirs. Les LCMs sont réfrigérés pour éviter une surchauffe de l électronique qu ils contiennent. La figure III.2 montre la distribution de l horloge aux LCMs d un secteur. A chaque étage, un diviseur optique sépare le signal en 25% du signal vers le LCM de l étage et 75% du signal vers l étage suivant impliquant une perte d intensité à chaque étage. Au dernier étage, le diviseur optique est inversé et envoie 75% du signal vers l étage. Chaque secteur est relié au module de contrôle de la ligne (String Control Module -SCM-) par une fibre optique différente pour l horloge. Un générateur de 1 khz impose sa fréquence à un laser, lui-même relié par fibres optiques aux 15 modules optiques testés et à 1 module optique de référence. 15 fibres optiques sont reliées aux 15 modules optiques du secteur. Une seizième fibre relie le laser à un module optique d un LCM supplémentaire afin de contrôler la charge de la lumière émise par ce laser. Une photodiode déclenchée par le laser envoie un signal électrique utilisé comme référence pour l estimation du temps d émission de la lumière.\n73 1. MÉTHODES D ÉTALONNAGE 73 FIGURE III.2 Schéma du cablage optique de l horloge d un secteur. On distingue deux fibres, l une passe directement du bas d un secteur au secteur suivant, alors que la distribution dans un secteur passe par tous les diviseurs optiques induisant des retards.\n74 74 CHAPITRE III. ETALONNAGE DES LIGNES DU DÉTECTEUR ANTARES Une fois le secteur intégré et testé il est placé sur la palette de transport de la ligne qui sert à acheminer la ligne jusqu au site de préparation au déploiement à la Seyne-sur-mer (Fosslev Marine). FIGURE III.3 Dispositif d étalonnage d un secteur. Lors de l étalonnage à terre, différentes prises de données sont effectuées : Déclenchement par un ordre extérieur aléatoire (runs piédestaux pour la charge, paragraphe 3.1). Prises de données de bruit de fond pour la détermination de la charge du photo-électron (paragraphe 3.2). Prise de données avec une forte intensité laser pour l étalonnage en temps (paragraphe 2). In situ, les étalonnages en charge et en temps peuvent partiellement être choisis à l aide des données standard. Des prises de données spéciales sont réalisées de façon régulière : les prises de données avec les balises optiques (paragraphe 2.2.c). les prises de données 40 K, détaillées en II.2.5. Ces prises de données permettent de valider l étalonnage en temps sur un étage et peuvent révéler les variations de gain des modules optiques (paragraphe 3.3). Les programmes développés pour l étalonnage en laboratoire peuvent être utilisés in situ en utilisant les prises de données utilisées pour les études de physique.\n75 2. ETALONNAGE EN TEMPS 75 2 Etalonnage en temps L étalonnage en temps du détecteur consiste à déterminer les constantes des convertisseurs tempstension (Time to Voltage Converter -TVC-) de chacun des ARSs et les constantes de temps entre les différents étages (télé-propagation des signaux). Le TVC discrétise le temps de façon linéaire. L étalonnage consiste d une part à extraire la pente et l ordonnée à l origine de la fonction de transfert linéaire temps-tvc. Le démarrage du comptage de l horloge initialise le début de chaque prise de données. Or le début de la prise de données est différent selon les étages à cause des délais de propagation dans les fibres optiques. Donc la définition de cette origine est différente pour chaque étage selon : la longueur des fibres optiques caractérisée par les retards de l horloge. l électronique embarquée : les résidus temporels, T 0, mesurent la différence de temps entre l instant réel de mesure et l instant attendu en connaissant la longueur des fibres optiques. L étalonnage consistera donc, d autre part, à mesurer ces résidus temporels. 2.1 Etalonnage du TVC Un convertisseur analogique-numérique (Analog to Digital Converter -ADC-) code idéalement l information du temps sur 256 canaux et la détermination de la pente et de l ordonnée à l origine du TVC est alors essentielle. Les TVCs n étant pas tous identiques, il est nécessaire de connaître les minima et maxima de chacun des TVCs, la pente étant obtenue avec la formule : TV C pente = dynamique du TVC TV C max TV C min avec dynamique du TVC = 50 ns (III.1) Pour un signal de bruit de fond aléatoire, les distributions des TVCs devraient être plates. A cause d un effet de non linéarité différentielle (section 2.4), ces distributions (figure III.4) s avèrent chevelues. Par convention, le TVC minimum (maximum) est alors défini comme le premier (dernier) bin dont la valeur est supérieure à % de la valeur moyenne de la population. Ces valeurs peuvent être obtenues avec n importe quel type de prise de données, en ayant vérifié que les résultats étaient identiques Mean RMS FIGURE III.4 Distribution d une rampe d un Time to Voltage Converter. FIGURE III.5 Intervalle décrit par le TVC en ns : 50 TV C max TV C min\n76 76 CHAPITRE III. ETALONNAGE DES LIGNES DU DÉTECTEUR ANTARES Pour être acceptés, le minimum du TVC doit être compris entre et 80 et le maximum du TVC entre 180 et 240. En effet le TVC est défini sur 50 ns, dans l intervalle accepté, le TVC peut donc être défini sur un intervalle compris entre 0 canaux (180-80) et 230 canaux (240-), soit un facteur deux possible sur l intervalle des TVCs. Si les minima et maxima ne sont pas dans les intervalles définis, il est possible d ajuster l intervalle du TVC en utilisant les paramètres slow control de l ARS. Les résultats des TVCs minima et maxima obtenus à terre pour la ligne 7 et le module optique de référence sont présentés à la figure III.6. La distribution des pentes issues de la calibration des TVCs est représentée figure III.5. FIGURE III.6 Distribution des minima et maxima des TVCs de la ligne 7 lors de l étalonnage à terre. Le nombre d entrée correspond au nombre d ARS multiplié par 2 car chaque ARS utilise deux rampes de conversion TVC en interne. Une fois le TVC calibré, nous pouvons remonter au temps d arrivée d un coup sur le photo-multiplicateur : t = TS + TV C TV C min TV C max TV C min 50 (III.2) avec T S la valeur du timestamp (c est-à-dire le comptage de l horloge en unité de 50 ns), T V C la valeur du TVC, TV C max la valeur maximum du TVC et TV C min la valeur minimum du TVC. 2.2 Décalage temporel des étages Il existe une différence de temps entre les délais induits par le chemin optique connu et ce que l on mesure effectivement : c est le résidu temporel, appelé T 0. Cette différence est due à la réponse de l électronique et au temps de transit dans le photo-multiplicateur. Elle diffère peu d un module optique à un autre (quelques ns). 2.2.a Horloge à Terre Une horloge commune installée à terre et synchronisée à la ms près par GPS est envoyée à chaque étage du détecteur. Cette horloge joue différents rôles :\n77 2. ETALONNAGE EN TEMPS 77 distribution synchrone des ordres de début et de fin d acquisition aux cartes horloges de tous les LCMs qui synchronisent à leur tour l acquisition de chaque ARS. contrôle de la fréquence de remise à zéro du compteur de timestamp de chaque ARS (Reset TimeStamp, RTS). L horloge donne le début de la prise de données des runs. Entre le bas et le haut d une ligne, la prise de données ne sera pas démarrée aux mêmes moments pour tous les LCMs à cause du temps de propagation de l horloge. Ces temps de propagation sont mesurés par le système d horloge par aller-retour de la lumière entre le système à terre et chaque étage. Comme on peut le voir à la figure III.7, les retards de l horloge croissent linéairement dans chaque secteur. En revanche on remarque des marches entre les différents secteurs. Ces marches s expliquent à l aide du schéma III.2. La distribution de la fibre d horloge dans un secteur passe par les quatres diviseurs optiques du secteur, alors que celle qui monte au secteur suivant traverse directement. Au bout de cinq étages les retards des fibres passées par les diviseurs sont plus grands que les retards des fibres qui n y sont pas passées. 2.2.b Détermination des T 0 à Terre Pour mesurer les T 0, on utilise les coups liés à l émission du laser en faisant la différence entre l instant où la photodiode détecte un signal et celui de la réponse du photo-multiplicateur, en prenant en compte les délais d horloge et les délais de propagation dans les fibres optiques connectées au module optique. Par convention, on référence les T 0 au premier module optique de la ligne qui aura un T 0 nul. Ainsi : T 0 =< T PMi T Photodiode > < T PM0 T Photodiode > avec T PMi le temps enregistré par le module optique i, T PM0 le temps relatif au premier module optique de la ligne et T Photodiode le temps obtenu avec la photodiode. Les T 0 de la ligne 7 sont représentés à titre d exemple à la figure III.8. La figure III.9 montre que la distribution des T 0 obtenus à terre n est pas centrée en 0 ce qui reflète le choix du premier module optique de la ligne comme module optique de référence. Les données utilisées pour mesurer les T 0 sont prises avec des hautes intensités pour s affranchir de l effet de Walk (2.3). Il est aussi possible de mesurer les T 0 des balises LED ( T 0 60 ns qui correspond à la différence de temps de transit entre les photo-multiplicateurs de " et le petit photo-multiplicateur associé aux balises LED) afin de les prendre en compte dans la détermination in situ présentée au paragraphe suivant. 2.2.c Contrôle in situ Le système de balises optiques : Les LED Optical Beacons -LOB- (balises LEDs) permettent d illuminer les modules optiques du détecteur avec un temps d émission des photons connu à la ns près . Les balises optiques se situent sur les étages 2, 9, 15 et 21 de chaque ligne. Ces LOBs envoient des flashs de lumière bleue. Afin de contrôler le temps d émission de la lumière un petit photo-multiplicateur est placé dans ce système. On calcule alors les T 0 in situ comme la valeur moyenne des résidus ( t) obtenus par la formule suivante  : t = (t ARS T 0ARS ) (t LOB T 0LOB ) d(lob,ars) c eau (III.3) où les temps t ARS et t LOB sont corrigés par les T 0 correspondants mesurés à terre. Le terme d(lob,ars) est la distance entre la LOB et l ARS, et c eau la vitesse de la lumière dans l eau. A la figure III., la différence entre les t de deux modules optiques d un même étage est représentée. Si cette différence est supérieure à 3 ns, elle donne lieu à une correction des T 0 mesurés à Terre. Un autre système, consistant\n78 78 CHAPITRE III. ETALONNAGE DES LIGNES DU DÉTECTEUR ANTARES Graph phase de l horloge en ns numéro d étage de la ligne FIGURE III.7 Retards de l horloge pour la ligne 5 du détecteur ANTARES. FIGURE III.8 Il y a 6 ARSs par étage et 25 étages, on obtient donc 150 T 0 exprimés en ns pour la ligne 7. FIGURE III.9 Histogramme des T 0, cette distribution n est pas centrée en 0 reflétant le choix du module optique de référence.\n79 2. ETALONNAGE EN TEMPS 79 en un laser de même caractéristique que celui utilisé à l intégration et fixé sur les pieds des lignes 7 et 8, permet d illuminer simultanément tous les pieds de ligne du détecteur. Il permet donc de contrôler les dérives possibles des premiers étages des lignes. FIGURE III. Différences en temps (t ARS T 0ARS ) (t LOB T 0LOB ) pour deux modules optiques du troisième étage de la ligne 1 illuminés par la balise LED installée au deuxième étage de la ligne 2 . Etalonnage in situ avec la désintégration du 40 K de l eau : La désintégration du 40 K de l eau suit la réaction suivante : 40 K β 40 Ca + e + ν e (III.4) L électron issu de cette réaction rayonne de la lumière Tcherenkov détectable avec les modules optiques du détecteur au taux de 30 khz environ. Une fraction de ces désintégrations peut être vue simultanément par deux modules optiques, donnant lieu à une coïncidence dont le taux est de l ordre de 13 Hz. On calcule alors les différences en temps entre ces deux modules optiques de l étage (figure III.11). Des prises de données sont consacrées à cette étude (voir II.2.5). De plus, en regardant le taux de coïncidences entre deux modules optiques du même étage, on peut déterminer si le réglage des seuils des ARSs est correct ou non (voir 3.5). Si le réglage des seuils est correct, les taux de coïncidences entre deux modules optiques du même étage sont les mêmes sur les différents étages et les différentes lignes comme le montre la figure III.12 . 2.3 Effet de Walk Lorsqu un signal franchit le seuil L0, on mesure le temps de l impulsion au moment du passage du seuil. La grande majorité des impulsions a une charge d un photo-électron. Mais à grandes amplitudes, le temps du passage du seuil sera avancé d un certain temps noté t sur la figure III.13 . Les variations de l amplitude entraînent une variation du temps du passage du seuil (figure III.14). En ajustant la variation du temps en fonction de la charge, on peut corriger cet effet. On prend comme référence le temps de passage du seuil pour les grandes amplitudes et on corrige alors le temps de passage du seuil des coups de basse amplitude.\n80 80 CHAPITRE III. ETALONNAGE DES LIGNES DU DÉTECTEUR ANTARES FIGURE III.11 Différence en temps entre deux des trois modules optiques d un étage. Cette différence en temps permet également de modifier les T 0 déterminés à terre si nécessaire. FIGURE III.12 Taux des coïncidences entre deux modules optiques d un même étage pour les lignes 1 à 5. Les modules optiques pour lesquels les taux de coïncidence sont marginaux doivent être réglés (paragraphe 3.5) FIGURE III.13 Principe de l effet de walk. LCM_29 ARS_2 FIGURE III.14 Corrélation de la différence en temps de passage du seuil avec la charge de l impulsion, la ligne pleine étant la différence en temps théorique pour cet effet Hits rampe A du TVC de l ARS 2 du LCM 29 FIGURE III.15 Effet de non linéarité différentielle : Histogramme d un TVC.\n81 3. ETALONNAGE DE LA CHARGE Non-linéarité différentielle des ADCs La valeur du TVC est déterminée par un ADC qui effectue une comparaison par dichotomie. Celuici compare si l amplitude de la rampe est en-dessous ou au-dessus de 256/2 = 128 (une tension de référence). S il est au-dessus de 128, il compare avec /2 = 192 ; s il est en dessous de 128, il compare avec 128/2 = 64, etc. Le résultat de ce codage est le TVC. Il peut arriver que le comparateur soit trop lent à atteindre sa valeur de référence (chargement capacitif). Au moment t com de la comparaison, s il n a pas atteint cette valeur, il compare l amplitude du signal d entrée à une valeur inférieure à la valeur de référence. La hauteur du signal d entrée peut alors être surestimée. Lors des comparaisons suivantes, il trouvera la valeur du signal d entrée toujours inférieure aux valeurs de comparaisons et mettra tous les bits suivants à 0. Certains canaux de l ADC sont donc privilégiés de la même façon que si la largeur des canaux de l ADC n était pas fixe. On observe sur la distribution du TVC que certains canaux sont remplis et les canaux adjacents sont vides. Cet effet étant amplifié par de grandes vitesses de séquençage de l ADC, il suffirait de baisser la vitesse de séquençage pour limiter cet effet. Apparaîtrait alors un autre effet : le bit flip. 2.5 Le bit flip Le phénomène de bit flip est une particularité inexpliquée de l ARS, que nous allons décrire brièvement. Il correspond à un passage inopiné de 0 à 1 d un bit de timestamp de l ARS, les transitions de 1 à 0 arrivant beaucoup plus rarement. Les bits affectés par ce changement sont ceux de rang impair et les changements apportés sont alors +2 1, +2 3, +2 5,..., mais la transition n a lieu au rang n que si le bit de rang n + 1 est 1. La probabilité de transition augmente si en plus le bit de rang n 1 vaut 1. Plus le taux de données est élevé sur l ARS, plus la probabilité du bit flip augmente. La probabilité du bit flip dépend de l ARS mais elle reste faible ( 4 pour des fréquences de lecture de 25 MHz allant jusqu à quelques pourcents pour des fréquences de 16 MHz). Pour cette raison, on a préféré conserver une fréquence de lecture haute, au prix d une dégradation de la précision sur la mesure du temps (2.4). 3 Etalonnage de la charge L étalonnage en charge du détecteur consiste à déterminer la pente et l ordonnée à l origine du convertisseur amplitude-tension (Analog Voltage Converter -AVC-). Elle repose donc sur l hypothèse de linéarité de la fonction de transfert des AVC vérifiée en laboratoire. Pour cela, deux mesures sont nécessaires : l idéal serait d utiliser une source de lumière calibrée d intensité variable. En pratique on dispose du bruit de fond optique, qui fournit en moyenne un photo-électron. Des prises de données spéciales permettent par ailleurs d energistrer la réponse du détecteur à un signal nul (piédestal) en déclenchant aléatoirement l acquisition. On utilise donc les charges de 0 et de 1 photo-électron en identifiant la valeur correspondante de l AVC. 3.1 Le piédestal Afin de mesurer le piédestal, on envoie un signal sur l ARS à l aide d un ordre extérieur (à un temps aléatoire). Ce signal déclenche la conversion de l AVC à charge nulle. La valeur du seuil de l ARS est mise au maximum afin que seuls les rares signaux de très haute charge soient détectés en coïncidence : on pourra facilement les discriminer. La figure III.16 présente la distribution typique d un AVC lors d une prise de données de piédestal. Les événements mesurés aux plus grandes valeurs de l AVC proviennent des désintégrations de 40 K dans la sphère. La figure III.17 représente la figure III.16 intégrée utilisée pour la détermination du piédestal.\n82 82 CHAPITRE III. ETALONNAGE DES LIGNES DU DÉTECTEUR ANTARES A faibles valeurs de l AVC, on distingue le piédestal puis un plateau indiquant qu aucune charge n a été détectée à ces valeurs de l AVC. Le piédestal est alors ajusté à l aide de la fonction erf. Cette méthode permet de diminuer les effets de non-linéarité différentielle décrite en 2.4. La distribution des piédestaux pour la ligne 9 est alors représentée figure III canal AVC FIGURE III.16 Distribution d un AVC. A petite charge, la valeur moyenne du pic correspond à la valeur du piédestal. A grande charge, les signaux sont dus à la désintégration du 40 K dans la sphère. FIGURE III.17 Distribution intégrée de l AVC : à faibles valeurs de l AVC on distingue le piédestal ; aux hautes valeurs de l AVC, les grandes charges sont détectées. FIGURE III.18 Distribution des piédestaux de la ligne 9. On observe une large dispersion des valeurs de piédestal ce qui justifie les calibrations individuelles des ARSs.\n83 3. ETALONNAGE DE LA CHARGE 83 Stabilité du piédestal Les valeurs de piédestal sont stables dans le temps comme l attestent par exemple les mesures effectuées en février, en mai et en août La différence entre les piédestaux du mois de février et d août est présentée à la figure III.19 et la différence entre les piédestaux du mois de mai et d août est présentée à la figure III.20. Les valeurs moyennes des distributions sont 0, 18 ± 0, 40 et 0, 13 ± 0, 30. On peut donc déduire qu il n y a pas eu de variation significative des piédestaux en fonction du temps sur la durée de la prise de données 5 lignes. 3.2 Le photo-électron La plupart des signaux lumineux (environ 99%) reçus par le photo-multiplicateur sont issus de la bioluminescence ou de la désintégration de 40 K dans l eau de mer. Ces signaux fournissent en moyenne un photo-électron unique III.21. On utilise donc des données de bruit de fond pour mesurer la position du photo-électron en canaux d AVC ajustée par une gaussienne (on ajuste en réalité une fonction erf sur la distribution à l instar du piédestal). Avec les réglages de slow control de l ARS choisi, la distance entre le photo-électron et le piédestal se situe autour de canaux (figure III.23) ce qui assure à la fois une bonne précision sur la mesure du pic du photo-électron et une grande gamme dynamique. La gamme dynamique de la charge mesurée par l ARS est définie par : I = 255 C 0pe C 1pe C 0pe 20 pe (III.5) avec C 0pe le canal du piédestal et C 1pe le canal du photo-électron. Comme le piédestal est aux alentours du canal 50 et que C 1pe C 0pe =, la gamme dynamique de l AVC est d environ 20 pe (figure III.22). 3.3 Baisse du gain des modules optiques Cette méthode d étalonnage permet de surveiller l évolution du gain des PMs. Le piédestal étant stable, une variation de la position du pic du photo-électron indique une variation du gain de l électronique ou de celui du PM. Les mesures du taux de coïncidence (paragraphe II.2.5) 40 K favorisent la deuxième hypothèse. La figure III.24 montre une diminution du canal du photo-électron de 0,2-0,3 canal par mois. Les études à l aide des données 40 K sur les taux de coïncidences entre les modules optiques d un même étage montrent une diminution de ce taux de coïncidences (figure III.25). Cet effet pourrait être dû à un vieillissement prématuré des PM dû à un taux élevé de la bioluminescence de l ordre de 170 khz en moyenne. 3.4 La diaphonie Les tensions qui donneront les TVCs et AVCs sont temporairement stockées dans des mémoires capacitives dans le pipeline analogique avant d être numérisées par les ADCs. Lors de ce stockage, il se produit un couplage capacitif entre les voies temps et charges dont résulte une diaphonie TVC vers AVC. En pratique, les valeurs de l AVC sont influencées par celles du TVC comme le montre la figure III.26. La relation entre l AVC et le TVC est linéaire, on corrige donc la valeur de l AVC de la façon suivante : ( AV C corr = AV C XT slope TV C (TV C ) max + TV C min ) (III.6) 2 avec XT slope la pente présentée à la figure III.26. Cet effet doit y être corrigée avant d effectuer la correction de walk (paragraphe 2.3). Par convention, aucune correction n est appliquée au milieu de la gamme du TVC. Ainsi pour une distribution plate des TVCs (bruit de fond), la valeur moyenne de la\n84 84 CHAPITRE III. ETALONNAGE DES LIGNES DU DÉTECTEUR ANTARES hpedstabfevaout Mean RMS hpedstabmaiaout Mean RMS différence de valeur de piédestal entre février 2007 et août différence de valeur de piédestal entre mai 2007 et août 2007 FIGURE III.19 La différence de piédestal entre le mois de février et d août On n observe pas de différence entre les deux étalonnages. FIGURE III.20 La différence de piédestal entre le mois de mai et d août On n observe toujours pas de différence entre les deux étalonnages. LCM_1313 ARS_4 Entries Mean Entries 150 Mean x RMS Valeur du photoelectron en canal AVC Mean y 61 RMS x RMS y valeur de l AVC Numero des ARS FIGURE III.21 A gauche : distribution du bruit de fond, le pic de la gaussienne déterminant la valeur du photo-électron. A droite : les valeurs du photo-électron pour les différents ARSs de la ligne 9.\n85 3. ETALONNAGE DE LA CHARGE 85 FIGURE III.22 Gamme dynamique des ARSs. FIGURE III.23 Distance entre le photo-électron et le piédestal (canaux AVC). FIGURE III.24 Baisse du canal du photo-électron en fonction du temps. FIGURE III.25 Diminution du taux de coïncidence du 40 K en fonction du temps.\nQuelques aspects de l intégration du premier prototype de ligne d ANTARES\nQuelques aspects de l intégration du premier prototype de ligne d ANTARES Stage expérimental, MIP deuxième année Cédric Roux Laboratoire d accueil : Centre de Physique des Particules de Marseille 163,\nLe savoir-faire du Centre d Études Nucléaires de Bordeaux-Gradignan au service d une mission spatiale internationale\nLe savoir-faire du Centre d Études Nucléaires de Bordeaux-Gradignan au service d une mission spatiale internationale La naissance de la mission internationale GLAST Observer des étoiles et des planètes,\nUniversité Paris XI THESE DE DOCTORAT. pour obtenir le titre de DOCTEUR DE L UNIVERSITE Paris XI Spécialité : astroparticules.\nUniversité Paris XI THESE DE DOCTORAT pour obtenir le titre de DOCTEUR DE L UNIVERSITE Paris XI Spécialité : astroparticules Présentée par Damien DORNIC Institut de Physique Nucléaire d Orsay Développement