Analytical model for the photomultiplier single photoelectron response including the electron back-scattering contribution

Cet article propose un modèle analytique décrivant la réponse d'un électron photounique dans les photomultiplicateurs, qui intègre physiquement le phénomène de rétrodiffusion des électrons au premier dynode et a été validé expérimentalement sur deux dispositifs.

L'essentiel

Titre : Le Grand Voyage des Électrons : Comment on a appris à compter les photons avec précision

Imaginez que vous êtes dans une salle de concert immense et sombre. Soudain, un seul photomètre (un flash de lumière) éclaire la scène. Votre travail est de compter exactement combien de photons (les particules de lumière) ont touché le mur. Pour cela, vous utilisez un détecteur spécial appelé Photomultiplicateur.

Ce détecteur est comme une chaîne de montagnes russes géante. Quand un photon arrive, il libère un petit électron (le passager). Cet électron doit descendre la chaîne, en sautant de plateforme en plateforme (les "dynodes"), et à chaque saut, il enclenche une réaction en chaîne : un électron devient deux, deux deviennent quatre, quatre deviennent huit, etc. À la fin, un seul électron est devenu des millions, créant un signal électrique que nous pouvons mesurer.

Le Problème : Le Bruit de Fond et les Électrons "Timides"

Jusqu'à présent, les physiciens avaient un modèle pour décrire ce signal. Ils disaient : "Voici le pic principal (les millions d'électrons), et le reste, c'est juste du bruit." C'était un peu comme dire : "La montagne est là, et tout ce qui est en dessous, c'est de la poussière."

Mais en réalité, ce n'est pas si simple. Certains électrons, au premier saut (la première plateforme), ne font pas le grand saut complet. Ils ricochent (comme une balle de tennis qui rebondit sur un mur) ou perdent un peu d'énergie. Ils ne génèrent pas des millions d'électrons, mais seulement quelques milliers.

Dans le passé, on traitait ces signaux "mi-pleins" comme du bruit aléatoire. Mais ce papier dit : "Non ! Ce n'est pas du bruit, c'est une physique précise !"

La Nouvelle Découverte : La Carte au Trésor

Les auteurs de ce papier ont créé une nouvelle recette mathématique (un modèle analytique) pour décrire exactement ce qui se passe.

1. Les Électrons "Gagnants" (Fully Amplified) : Ce sont les passagers qui font le grand saut parfait. Ils arrivent au bas de la montagne avec un gros sac de millions d'électrons. C'est le pic principal.
2. Les Électrons "Ricochets" (Partially Amplified) : C'est ici que la magie opère. Certains électrons heurtent la première plateforme et rebondissent (back-scattering). Ils perdent un peu d'énergie. Au lieu de générer 15 millions d'électrons, ils n'en génèrent que 10, ou 5, ou 2. Ils ne sont pas "bruits", ils sont juste des voyageurs fatigués.
   • L'analogie : Imaginez que vous lancez une balle dans un escalier. Parfois, elle tombe directement sur la dernière marche (signal fort). Parfois, elle heurte la première marche, rebondit, perd de la vitesse, et finit par s'arrêter au milieu de l'escalier (signal faible). Le nouveau modèle permet de prédire exactement où elle va s'arrêter.
3. Les Électrons "Perdus" ou "Pré-pulses" : Il y a aussi des cas où l'électron arrive trop tôt ou trop tard, ou où un photon touche directement la première plateforme sans passer par la photocathode. Le modèle décrit aussi ces cas bizarres.

Pourquoi est-ce important ?

Avant, pour analyser ces signaux "mi-pleins", les scientifiques devaient inventer des formules compliquées et arbitraires (des "fonctions ad hoc") pour coller les courbes. C'était comme essayer de deviner la forme d'un nuage en regardant juste une ombre.

Maintenant, avec ce nouveau modèle :

• C'est plus précis : On utilise les vraies lois de la physique (la façon dont les électrons rebondissent) au lieu de deviner.
• C'est universel : Ils ont testé leur recette sur deux types de détecteurs différents (des modèles Hamamatsu) et ça a marché parfaitement dans les deux cas.
• C'est utile pour l'avenir : Des expériences géantes comme XENONnT (qui cherche la matière noire) ou LZ utilisent ces détecteurs. Pour trouver la matière noire, il faut être capable de distinguer un signal réel d'un faux signal avec une précision chirurgicale. Ce nouveau modèle aide à ne pas rater le moindre indice.

En Résumé

Imaginez que vous essayez de compter des gouttes de pluie tombant sur un toit.

• L'ancienne méthode : "On entend le bruit des grosses gouttes. Le reste, c'est juste le vent."
• La nouvelle méthode (ce papier) : "Attendez ! Les petites gouttes qui glissent sur les tuiles sans faire de bruit ont une trajectoire précise. Si on comprend comment elles glissent, on peut compter toutes les gouttes, même les plus petites, et reconstruire exactement la tempête."

Les auteurs ont donc écrit une "carte" mathématique qui explique pourquoi certains signaux sont plus faibles que d'autres, en se basant sur le fait que les électrons aiment parfois faire des ricochets. C'est une avancée majeure pour rendre nos yeux électroniques encore plus précis pour observer l'univers.

Résumé technique

1. Problématique

Les photomultiplicateurs (PMT) restent des capteurs de choix pour la détection de très faibles niveaux de lumière sur de grandes surfaces. Une caractéristique fondamentale pour leur utilisation est la réponse à un seul photoélectron (SPE). Cependant, la modélisation de cette réponse, en particulier dans la région spectrale située entre le pic de bruit électronique (pedestal) et le pic principal des photoélectrons entièrement amplifiés, pose problème.

Les modèles existants utilisent souvent des fonctions ad hoc (comme des gaussiennes modifiées exponentiellement) pour décrire cette région intermédiaire, l'attribuant généralement au « bruit ». Or, cette région contient des signaux physiques réels, notamment dus au rétrodiffusion (back-scattering) des photoélectrons primaires sur la première dynode, qui ne déposent qu'une fraction de leur énergie et génèrent donc un gain réduit. L'absence de modèle physique rigoureux pour cette composante limite la précision des calibrations et des simulations.

2. Méthodologie

Les auteurs proposent un modèle analytique complet de la réponse SPE, basé sur la description physique du processus de rétrodiffusion des électrons sur la première dynode, tel que décrit dans le « Photomultiplier Handbook » d'A. G. Wright.

Le modèle décompose le processus d'amplification en deux étapes distinctes :

1. La première dynode : Elle détermine la forme globale du spectre. Les photoélectrons primaires sont divisés en deux catégories :
   • Photoélectrons entièrement amplifiés (FA) : L'électron dépose toute son énergie, générant un nombre moyen $G_1$ d'électrons secondaires suivant une statistique de Poisson.
   • Photoélectrons partiellement amplifiés (PA) : Une fraction $\eta$ (environ 30 %) subit une rétrodiffusion inélastique. L'électron perd une partie de son énergie ($E_p - E'_p$) dans la dynode, générant moins d'électrons secondaires. La distribution d'énergie des électrons rétrodiffusés est modélisée comme uniforme, conduisant à une distribution de probabilité discrète spécifique pour le nombre d'électrons secondaires.
2. La chaîne des dynodes suivantes ($N-1$ dynodes) : Elles lissent la distribution discrète de la première dynode en une distribution continue (approximée par une gaussienne) et contribuent aux signaux de très faible charge.

Le modèle analytique final ($\mathcal{P}_{SPE}$) est une somme pondérée de quatre composantes :

• Un pic entièrement amplifié décrit par une distribution de Poisson continue mise à l'échelle.
• Une composante partiellement amplifiée décrite par une fonction en forme de « boîte arrondie » (combinant des fonctions d'erreur erf), résultant de la convolution de la distribution discrète des électrons PA avec le bruit électronique et l'amplification ultérieure.
• Des pré-impulsions (pre-pulses) dues à des photons atteignant directement la première dynode, modélisées par une distribution de Poisson à gain réduit.
• Des signaux de très faible charge (exponentiels) dus à des trajectoires erronées ou d'autres effets de multiplication.

3. Contributions Clés

• Modélisation physique de la rétrodiffusion : Contrairement aux approches empiriques, ce modèle dérive analytiquement la forme de la distribution des électrons partiellement amplifiés à partir des principes physiques de la rétrodiffusion inélastique sur la première dynode.
• Fonction analytique complète : Les auteurs fournissent des expressions mathématiques fermées pour toutes les composantes du spectre SPE (pics FA, PA, pré-impulsions, bruit), ne dépendant que d'un petit nombre de paramètres intrinsèques au PMT (gain $G_1$, fraction de rétrodiffusion $\eta$, résolution de la chaîne $R$, etc.).
• Validation expérimentale rigoureuse : Le modèle a été testé sur deux types de photomultiplicateurs différents (Hamamatsu R5912-100 de 8 pouces et Hamamatsu 6233 de 3 pouces) utilisant une chaîne de lecture à très faible bruit (LNA + numériseur rapide), permettant de résoudre des structures spectrales à très faible charge souvent masquées par le bruit électronique.

4. Résultats

Les expériences ont été menées à l'INFN Laboratori Nazionali del Gran Sasso avec des impulsions laser sub-nanoseconde.

• Adéquation du modèle : Le modèle ajuste parfaitement les données expérimentales pour les deux types de PMT, avec des valeurs de $\chi^2/\text{ndf}$ proches de 1 (ex: 1.04 pour le R5912-100).
• Paramètres extraits :
  • Pour le R5912-100, la fraction d'électrons rétrodiffusés ($\eta$) est d'environ 27 %. La composante exponentielle à faible charge est significative (~4 %), tandis que les pré-impulsions sont négligeables.
  • Pour le 6233, la fraction de rétrodiffusion est plus faible (17 %), mais une contribution significative de pré-impulsions (7 %) est observée, ce qui est cohérent avec la géométrie du tube.
• Dépendance aux conditions : Les paramètres de gain ($G_1$, $f$) dépendent linéairement de la haute tension (HV), tandis que les fractions d'événements ($\eta$, $A_{Exp}$) sont indépendantes de la HV mais dépendent de la position d'illumination sur la photocathode.
• Résolution spectrale : Le modèle permet de distinguer clairement les structures temporelles associées aux différents types d'événements (ex: les pré-impulsions arrivent ~10 ns plus tôt que le signal principal).

5. Signification

Ce travail offre une description physiquement motivée et robuste de la réponse SPE des photomultiplicateurs, éliminant le besoin de paramétrisations empiriques arbitraires pour la région intermédiaire du spectre.

• Précision accrue : En modélisant correctement la queue de distribution à faible charge, ce modèle permet une détermination plus précise du gain effectif, de la résolution et de l'acceptance des détecteurs.
• Applications : Le modèle est particulièrement pertinent pour les expériences de physique des particules et d'astroparticules (comme XENONnT ou LZ, mentionnées dans l'article) qui utilisent des PMT pour la détection de scintillations faibles et nécessitent une caractérisation photonique de haute précision pour la simulation et l'analyse des données.
• Généralité : Le modèle, avec son nombre limité de paramètres intrinsèques, est applicable à différents types de PMT et conditions opératoires, facilitant son intégration dans les cadres de calibration et de simulation.