Phenomenological Modeling of the $^{163}$Ho Calorimetric Electron Capture Spectrum from the HOLMES Experiment

Cet article présente une analyse phénoménologique complète du spectre de capture électronique calorimétrique de $^{163}$Ho mesuré par l'expérience HOLMES, modélisant le spectre comme une somme de résonances de Breit-Wigner et de continua de secousse pour fournir une description précise de la région de l'extrémité cruciale pour la détermination de la masse du neutrino.

L'essentiel

🕵️‍♂️ L'Enquête : Chasser le fantôme de la masse du neutrino

Imaginez que vous essayez de peser un fantôme. C'est exactement ce que fait l'expérience HOLMES. Les physiciens cherchent à mesurer la masse du neutrino, une particule fantôme qui traverse tout (y compris votre corps) sans presque rien laisser derrière elle.

Pour y parvenir, ils utilisent un atome spécial : le Holmium-163. C'est un peu comme un piège à souris atomique. Quand cet atome se transforme (il "capture" un électron), il libère de l'énergie. Si le neutrino a une masse, il emporte un peu de cette énergie avec lui, comme un voleur qui s'échapperait avec un bout du butin.

Le but du papier est de comprendre parfaitement toute la scène du crime (le spectre d'énergie) pour savoir exactement combien le "voleur" (le neutrino) a pu prendre.

---

🎻 Le Concert des Atomes : Une symphonie complexe

Jusqu'à présent, les scientifiques pensaient que l'atome de Holmium, en se transformant, jouait une mélodie simple : une série de notes claires et distinctes (comme des cloches qui sonnent). C'était ce qu'on appelait l'approximation "trou unique" (un seul électron manque).

Mais le papier de HOLMES nous dit : "Attendez, c'est beaucoup plus compliqué !"

En réalité, quand l'atome change, c'est comme si un chef d'orchestre avait donné un coup de sifflet soudain. Non seulement la note principale sonne, mais tout l'orchestre réagit :

1. Les résonances (Breit-Wigner) : Ce sont les notes principales, les cloches qui sonnent fort.
2. Les "Shake-up" (Secousses) : Imaginez que le tremblement de la terre fait tomber des meubles sur le sol. Ce sont des électrons qui sont "secoués" et passent d'une place à une autre, créant des notes supplémentaires.
3. Les "Shake-off" (Éjections) : C'est encore plus violent. Certains électrons sont tellement secoués qu'ils sont éjectés de l'atome, créant un bruit de fond continu, comme une pluie fine qui tombe sur les notes principales.

L'analogie du gâteau :
Imaginez que vous essayez de goûter un gâteau (l'énergie totale).

• L'ancienne théorie disait : "Le gâteau est juste une couche de chocolat."
• La nouvelle théorie (celle de HOLMES) dit : "Non, le gâteau a une couche de chocolat, mais aussi des éclats de noix, du sucre cristallisé, et une crème fouettée qui coule sur les côtés. Si vous ne comprenez pas la recette exacte, vous ne saurez jamais si le gâteau est vraiment sucré ou s'il vous manque un ingrédient (la masse du neutrino)."

---

🔍 La Méthode : Déballer le cadeau (Déconvolution)

Le problème, c'est que les détecteurs de HOLMES sont comme des oreilles un peu bouchées. Ils entendent le concert, mais les notes se mélangent un peu à cause du bruit de fond et de la qualité de l'oreille.

Les auteurs ont utilisé une technique mathématique puissante appelée "déconvolution" (ou "déballage").

• Avant : Ils avaient un signal flou, un gros tas de données où tout se mélangeait.
• L'astuce : Ils ont utilisé un algorithme (une sorte de détective mathématique) pour "nettoyer" le signal. Ils ont dit : "Si je sais comment mon détecteur déforme le son, je peux inverser le processus pour retrouver le son original."

Résultat : Ils ont pu voir clairement les 25 composantes différentes de ce concert atomique (les cloches, les chocs, les éjections) et mesurer leur hauteur et leur volume avec une précision incroyable.

---

🤝 Le Duel : Théorie vs Réalité

Les auteurs ont comparé leur "carte du concert" (leurs données réelles) avec les meilleures prédictions des ordinateurs (les calculs ab initio).

• Le verdict : C'est un match serré ! Les ordinateurs ont fait du bon travail pour prédire la forme générale, mais ils ont raté certains détails fins. C'est comme si un compositeur avait écrit une partition magnifique, mais qu'il avait oublié de noter quelques trilles ou des variations subtiles que les musiciens réels jouent spontanément.
• L'importance : Cette comparaison est cruciale. Elle montre que nous ne pouvons pas nous fier aveuglément aux théories actuelles pour mesurer la masse du neutrino. Il faut utiliser les données réelles pour corriger les théories.

---

🎯 Pourquoi c'est important pour nous ?

Pourquoi se soucier de toutes ces notes et de ces secousses d'électrons ?

1. La précision du vol : Pour savoir si le neutrino a volé de l'énergie (et donc a une masse), il faut connaître la valeur exacte de l'énergie totale disponible. Si on ne comprend pas les "bruits de fond" (les secousses d'électrons), on risque de confondre un bruit avec le vol du neutrino.
2. L'avenir : Ce papier fournit la "recette de cuisine" parfaite pour les futures expériences. Grâce à ce modèle, les prochains détecteurs pourront simuler exactement ce à quoi ils s'attendent, éliminer les erreurs et enfin peser ce fantôme insaisissable.

En résumé :
Ce papier est comme un manuel d'instructions ultra-détaillé pour comprendre le bruit d'une explosion atomique. En décortiquant chaque petit bruit et chaque résonance, les scientifiques de HOLMES nous donnent les outils nécessaires pour enfin répondre à l'une des plus grandes questions de la physique : Quelle est la masse du neutrino ?

Résumé technique

1. Problématique et Contexte

La mesure de la masse absolue du neutrino électronique ($m_\nu$) est un enjeu majeur en physique des particules et en cosmologie. L'expérience HOLMES vise à déterminer cette masse en étudiant le spectre de désintégration par capture électronique (CE) de l'isotope Holmium-163 ($^{163}\text{Ho}$) à l'aide de microcalorimètres.

Le défi principal réside dans la modélisation précise du spectre de désexcitation calorimétrique ($E_c$). Le modèle théorique standard, basé sur l'approximation d'un seul trou (Single Hole Approximation - SHA) et la règle d'or de Fermi, décrit le spectre comme une somme de résonances de Breit-Wigner correspondant aux trous dans les couches électroniques (M, N, O, etc.). Cependant, ce modèle simple s'avère insuffisant pour expliquer les données expérimentales à haute statistique, notamment :

• La présence de structures complexes et de queues spectrales non prédites.
• Les effets de shake-up (excitation d'un électron vers une orbitale liée) et de shake-off (ionisation d'un électron vers le continuum) qui créent des états à deux trous.
• La nécessité d'une description analytique précise de la région de l'extrémité (endpoint) du spectre, cruciale pour l'extraction de la masse du neutrino, car les effets de fond (piles-up, queues de basses énergies) peuvent masquer le signal.

2. Méthodologie

A. Expérience et Acquisition de Données

L'expérience HOLMES utilise des microcalorimètres à capteur à transition (TES) en couches minces Mo/Cu, couplés à des absorbeurs en or contenant l'isotope $^{163}\text{Ho}$ implanté par ionisation.

• Données : L'analyse porte sur des données à haute statistique collectées lors de trois campagnes de mesure.
• Calibration : Une procédure itérative basée sur l'apprentissage bayésien a été utilisée pour calibrer les pics principaux (M1, M2, N1, etc.) en utilisant les rayons X de fluorescence (Al, Cl) et en corrigeant les non-linéarités des détecteurs.
• Déroulement (Unfolding) : Pour éviter la distorsion des caractéristiques spectrales intrinsèques due à la résolution énergétique variable des ~1000 spectres individuels, une méthode de déroulement (unfolding) itérative basée sur le théorème de Bayes (via le framework PyUnfold) a été appliquée. Cela permet de reconstruire la distribution d'énergie vraie avant la sommation des spectres.

B. Modélisation Phénoménologique

Les auteurs proposent un modèle analytique complet pour ajuster le spectre déroulé :

1. Résonances (Breit-Wigner) : Les pics discrets sont modélisés par des fonctions de Breit-Wigner asymétriques. Une paramétrisation spécifique (Éq. 3.1-3.2) est utilisée pour capturer l'asymétrie forte observée, incluant un terme de suppression de la queue gauche ($p$) et des largeurs différentes à gauche et à droite ($\Gamma_L, \Gamma_R$).
2. Continus (Shake-off) : Les contributions continues dues aux processus de shake-off sont modélisées par des distributions analytiques dérivées des travaux de Levinger (adaptés pour les atomes de haute Z), convoluées avec la largeur naturelle des états à deux trous.
3. Ajustement Global : L'ajustement est réalisé en deux étapes :
   • Une optimisation initiale avec Minuit2 pour identifier les composants.
   • Une estimation bayésienne simultanée de tous les paramètres (103 paramètres libres pour 25 composants) utilisant l'algorithme Hamiltonian Markov Chain Monte Carlo (HMC) implémenté dans STAN.
   • Le modèle inclut le facteur d'espace des phases $(Q - E_c)^2$ et un terme de fond constant.

3. Résultats Clés

• Décomposition Spectrale : Le spectre déroulé est parfaitement reproduit par une combinaison de 19 pics Breit-Wigner (symétriques et asymétriques) et de 6 spectres de shake-off.
• Identification des Composantes : Les auteurs attribuent les composantes observées à des excitations à un trou (single-hole) et à deux trous (double-hole).
  • Le modèle identifie des excitations dans des couches non accessibles directement par la capture électronique dans le modèle SHA (ex: M3, M4, M5, N3, N4, N5), expliquées par des processus de relaxation inter-cœur (Coulomb scattering).
  • Des structures complexes sont résolues dans les groupes O, N et M, notamment des multiplets et des satellites de shake-up/shake-off.
• Comparaison avec la Théorie Ab Initio :
  • Une comparaison est faite avec des calculs ab initio récents (Brass et al., utilisant la formalisme de Kubo et les fonctions de Green).
  • Bien que les calculs ab initio reproduisent qualitativement la forme générale et l'asymétrie, ils ne correspondent pas parfaitement aux énergies de résonance et aux largeurs sans ajustements phénoménologiques.
  • Le modèle phénoménologique de HOLMES confirme la prédiction selon laquelle la région de l'extrémité est exempte de pics discrets, ne contenant que des queues de résonances et des contributions de shake-off.
• Paramètres de Capture : Les probabilités de capture pour les différentes couches (M1, M2, N1, etc.) sont recalculées en sommant les intensités relatives de toutes les excitations impliquant un trou de cœur donné. Les résultats sont cohérents avec les attentes théoriques mais affinés par les données expérimentales.

4. Contributions Majeures

1. Première analyse à haute statistique : Fourniture du premier spectre calorimétrique complet de la CE du $^{163}\text{Ho}$ avec une résolution et une statistique suffisantes pour révéler des structures fines.
2. Modèle analytique robuste : Développement d'un modèle paramétrique complet (pics + continu) capable de décrire l'ensemble du spectre de 30 eV jusqu'à l'extrémité $Q$. Ce modèle est essentiel pour les analyses futures de masse de neutrino.
3. Outil pour les études de sensibilité : La décomposition permet la génération de spectres "toy" Monte Carlo réalistes pour évaluer la sensibilité des futures expériences (ECHo, HOLMES Phase 2) et pour étudier les incertitudes systématiques liées aux effets du solide et du détecteur.
4. Validation de l'extrémité : Confirmation que la région de l'extrémité est lisse, ce qui valide l'approche de l'extraction de la masse du neutrino sans biais de structures non modélisées.

5. Signification et Impact

Ce travail établit une fondation robuste pour les futures expériences de mesure de la masse du neutrino utilisant la capture électronique de l'holmium.

• Précision de l'extrémité : Le modèle permet une description précise de la région critique près de $Q$, incluant les effets de fond (piles-up, queues de basses énergies), réduisant ainsi les incertitudes systématiques sur la détermination de $m_\nu$.
• Compréhension Physique : Il apporte des éclaircissements sur les processus atomiques complexes (shake-up, shake-off, relaxation inter-cœur) qui affectent le spectre, comblant le fossé entre les modèles théoriques simplifiés et la réalité expérimentale.
• Futur : Le modèle sera utilisé pour l'analyse des données en cours et pour la conception optimisée des détecteurs des prochaines générations d'expériences de neutrino.

En résumé, cette étude transforme la compréhension du spectre de désintégration du $^{163}\text{Ho}$ d'une approximation théorique simplifiée vers un modèle phénoménologique précis et validé expérimentalement, indispensable pour atteindre la précision requise dans la mesure de la masse du neutrino.