🔬 physics

A four parameters modified Gaussian function describing shape of asymmetric alpha lines

Cet article présente une nouvelle fonction analytique à quatre paramètres pour modéliser avec précision les pics alpha asymétriques, laquelle a été intégrée au logiciel établi ALF pour permettre la déconvolution de jusqu'à dix pics chevauchants.

Auteurs originaux : Jerzy Wojciech Mietelski

Publié 2026-07-09

📖 4 min de lecture☕ Lecture pause café

CC BY 4.0

Auteurs originaux : Jerzy Wojciech Mietelski

Article original sous licence CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ Ceci est une explication générée par l'IA de l'article ci-dessous. Elle n'a pas été rédigée ni approuvée par les auteurs. Pour une précision technique, consultez l'article original. Lire la clause de non-responsabilité complète

Imaginez que vous écoutiez une chorale de particules alpha chantant leur chemin vers un détecteur. Dans un monde parfait, chaque chanteur frappe exactement la même note, créant une belle courbe en cloche symétrique — une forme gaussienne classique. Mais dans le monde réel, les choses sont désordonnées. Tandis que ces minuscules particules voyagent de leur source vers le détecteur, elles percutent des objets : une « couche morte » sur la fenêtre d'entrée du détecteur et le matériau de la source elle-même. Ces collisions leur volent un peu d'énergie, mais pas toujours la même quantité.

Voyez cela comme des coureurs dans une course où certains trébuchent sur une pierre instable et d'autres non. Les coureurs qui ne trébuchent pas arrivent exactement à l'heure, maintenant le côté haute énergie du pic net et symétrique. Mais les coureurs qui trébuchent arrivent en retard, étirant le côté basse énergie du pic en une longue traîne désordonnée. Cela crée une forme « asymétrique » qui ressemble à une cloche qui aurait été traînée dans la boue d'un côté.

Pendant des décennies, les scientifiques ont lutté pour décrire cette forme désordonnée avec des mathématiques. Ils devaient généralement utiliser une boîte à outils compliquée composée de nombreuses formules différentes pour assembler les différentes parties de la courbe. Mais Jerzy Wojciech Mietelski, un chercheur de l'Institut de Physique Nucléaire de Pologne, a décidé d'essayer une approche différente. Il a introduit une fonction mathématique unique et ingénieuse qui agit comme un « changeur de forme ».

Cette nouvelle fonction est une gaussienne modifiée, mais avec une particularité. Elle ne possède que quatre boutons de réglage (paramètres) à tourner :

L'amplitude : La hauteur du pic.
La position : L'endroit où le pic se situe sur le graphique.
Le paramètre de forme 1 : Contrôle le comportement de la courbe sur le côté droit (le côté gaussien, propre).
Le paramètre de forme 2 : Contrôle la traîne longue et désordonnée sur le côté gauche.

Voici le tour de magie : sur le côté droit du pic, la fonction se comporte comme une cloche parfaite standard. Mais sur le côté gauche, elle se transforme en une formule plus complexe qui crée naturellement cette longue traîne exponentielle causée par la perte d'énergie. C'est comme avoir un seul morceau d'argile qui peut être lisse et rond d'un côté, mais étiré et effilé de l'autre, sans avoir besoin de coller deux formes différentes ensemble.

Mietelski n'a pas seulement imaginé cela ; il l'a mis en pratique dans un code informatique appelé ALF, qui tourne dans son laboratoire depuis 1993. Le processus se déroule en deux étapes, comme l'accordage d'une radio. D'abord, l'ordinateur examine un pic unique et clair pour déterminer les réglages exacts des paramètres de forme (le « côté boueux » et le « côté propre »). Une fois qu'il sait à quoi ressemble la « boue » pour cette expérience spécifique, il verrouille ces réglages. Ensuite, il utilise cette forme fixe pour démêler un fouillis de pics qui se chevauchent — jusqu'à 10 à la fois — en ajustant simplement leur hauteur et leur position.

L'article rapporte que cette méthode a été testée sur 30 ans et a analysé avec succès plus de 15 000 spectres. Elle a même passé des « tests de compétence » rigoureux où le laboratoire devait prouver qu'il pouvait mesurer correctement par rapport à d'autres experts. Les auteurs notent que les résultats sont si fiables que l'aire sous la courbe (qui indique la quantité de matériau radioactif présent) correspond presque parfaitement au nombre total de particules comptées.

Étonnamment, ce « changeur de forme alpha » n'était pas réservé qu'aux particules alpha. L'article mentionne que l'équipe a également essayé d'utiliser ces mêmes mathématiques pour séparer les signaux de rayonnement bêta d'un spectromètre à scintillation liquide. Dans un test impliquant le 90Sr (Strontium-90) en équilibre avec le 90Y (Yttrium-90), la méthode a réussi à séparer les signaux avec une marge d'erreur de seulement 3 %.

Ainsi, bien que l'article ne prétende pas que c'est un remède miracle pour tous les problèmes de la physique, il présente un outil robuste à quatre paramètres qui a fait ses preuves pendant des décennies. Il transforme un problème asymétrique et désordonné en un puzzle soluble, permettant aux scientifiques de voir clairement à travers le bruit des signaux qui se chevauchent, qu'ils étudient des particules alpha provenant de lichens de l'Antarctique ou des particules bêta provenant d'ossements d'animaux.

Résumé technique : Une fonction gaussienne modifiée à quatre paramètres pour les raies alpha asymétriques

Énoncé du problème
En spectrométrie alpha, les raies spectrales présentent fréquemment des formes asymétriques plutôt que des distributions gaussiennes idéales. Cette asymétrie provient principalement de la dégradation de l'énergie lorsque les particules alpha traversent la couche morte (fenêtre d'entrée) du détecteur à semi-conducteur et le matériau de la source elle-même. L'étendue de cette dégradation dépend de l'angle d'émission aléatoire de la particule par rapport au détecteur. Par conséquent, le côté basse énergie du pic devient plus large, tendant vers une décroissance exponentielle, tandis que le côté haute énergie reste gaussien, modifié uniquement par les fluctuations statistiques.

Historiquement, il n'existait aucune formule mathématique unique et généralement acceptée pour décrire cette forme. Les solutions pratiques reposent généralement sur l'application de fonctions à ensembles de paramètres multiples à des sections spécifiques de la raie spectrale. Bien que divers progiciels (par exemple, ALFUN-A, MicroSAMPO, ALFA) existent pour relever ces défis, la fonction analytique spécifique utilisée par le code informatique ALF, largement adopté, n'avait pas été détaillée dans la littérature avant ce travail.

Méthodologie
L'article présente une fonction analytique originale conçue pour décrire la forme asymétrique du pic alpha en utilisant seulement quatre paramètres libres. Cette fonction a été implémentée dans le code informatique ALF, un outil basé sur Fortran développé vers 1993 (modifié à partir d'un ancien code pour les rayons gamma, P.I.M.P.) pour l'évaluation des spectres alpha.

Le cœur de la méthodologie est la fonction gaussienne modifiée suivante, $N(x)$ , où $x$ représente le numéro de canal :

$N(x) = \alpha \cdot \exp\left( \frac{-(x-M)^2}{\beta \{|x-M| - (x-M)\} + \gamma} \right)$

Paramètres :

$\alpha$ : Amplitude (comptages au maximum).
$M$ : Position du maximum du pic.
$\beta$ et $\gamma$ : Paramètres de forme (tous deux $>0$ ) spécifiques à la géométrie de la source et du détecteur.

Comportement de la fonction :

Côté haute énergie ( $x > M$ ) : La fonction converge vers une distribution gaussienne standard.
Côté basse énergie ( $x < M$ ) : La fonction devient plus complexe, transitant vers une fonction exponentielle ordinaire à des distances significatives du maximum du pic ( $x \ll M$ ).
Rôles des paramètres : $\gamma$ régit la forme du côté droit du pic, tandis que $\beta$ et $\gamma$ régissent conjointement le côté gauche.

Procédure d'ajustement (déconvolution) :
Le code ALF utilise une procédure d'ajustement de fonction non linéaire pour déconvoluer jusqu'à 10 pics alpha chevauchants. Le processus s'exécute en deux étapes distinctes :

Détermination de la forme initiale : Une raie spectrale unique et bien définie (ou un fragment contenant le maximum) est ajustée pour déterminer les quatre paramètres ( $\alpha, M, \beta, \gamma$ ). Cela utilise un ajustement standard par moindres carrés non linéaires via 10 itérations d'une approximation linéaire (développement de Taylor).
Déconvolution de pics : Une fois que $\beta$ et $\gamma$ sont déterminés, ils sont traités comme des constantes fixes pour l'ensemble du spectre, car ils dépendent de la qualité de la source et de la géométrie source-détecteur plutôt que des pics individuels. La deuxième étape consiste à ajuster les amplitudes ( $\alpha_i$ ) et les positions ( $M_i$ ) pour tous les pics chevauchants simultanément. Ceci est réalisé en développant la somme des fonctions en une série et en résolvant pour les incréments du vecteur de variables à l'aide d'un algorithme de Cracovienne.

L'algorithme de Cracovienne est mis en évidence comme un avantage spécifique, car il calcule les incertitudes des paramètres ajustés au sein de la même procédure utilisée pour déterminer les valeurs des paramètres.

Contributions clés et résultats

Simplicité analytique : L'article introduit une fonction compacte à quatre paramètres qui modélise efficacement l'asymétrie complexe des raies alpha sans nécessiter d'ajustements multi-segments complexes.
Validation et longévité : Cette fonction est utilisée de routine depuis 1993. Au cours des 30 dernières années, elle a été appliquée à plus de 15 000 spectres. La validation a été effectuée en comparant les aires de pics calculées à la somme des comptages, tous les tests ayant été concluants.
Capacité de déconvolution : La méthode sépare avec succès des pics partiellement chevauchants, comme démontré avec des échantillons de Th de restes osseux d'animaux.
Utilité inattendue : Les auteurs notent que la fonction s'est avérée étonnamment utile pour analyser les spectres de rayonnement bêta (spécifiquement $^{90}\text{Sr}$ / $^{90}\text{Y}$ ) obtenus avec un spectromètre à scintillation liquide utilisant une échelle de canaux logarithmique, atteignant des résultats avec une marge d'erreur de 3 %.

Signification
L'article revendique sa signification principalement par la validation réussie et à long terme de cette fonction spécifique à quatre paramètres au sein du logiciel ALF. En formalisant la base analytique précédemment non documentée de ALF, ce travail fournit un outil mathématiquement défini et robuste pour la déconvolution de pics alpha chevauchants. La méthode offre un équilibre entre rigueur analytique et efficacité de calcul, permettant la détermination précise des aires de pics dans des spectres complexes où l'ajustement gaussien traditionnel échoue. L'utilisation de l'algorithme de Cracovienne garantit en outre que les incertitudes des paramètres sont calculées rigoureusement parallèlement à l'ajustement.

Résumé technique : Une fonction gaussienne modifiée à quatre paramètres pour les raies alpha asymétriques

Articles similaires