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

Resumo Técnico: Uma Função Gaussiana Modificada de Quatro Parâmetros para Linhas Alfa Assimétricas

Definição do Problema
Na espectrometria alfa, as linhas espectrais frequentemente exibem formas assimétricas em vez de distribuições gaussianas ideais. Essa assimetria surge primariamente da degradação de energia à medida que as partículas alfa atravessam a camada morta (janela de entrada) do detector de semicondutor e o próprio material da fonte. A extensão dessa degradação depende do ângulo de emissão aleatório da partícula em relação ao detector. Consequentemente, o lado de baixa energia do pico torna-se mais largo, aproximando-se de uma queda exponencial, enquanto o lado de alta energia permanece gaussiano, modificado apenas por flutuações estatísticas.

Historicamente, não houve uma fórmula matemática única e geralmente aceita para descrever essa forma. Soluções práticas normalmente dependem de conjuntos de funções de múltiplos parâmetros aplicados a seções específicas da linha espectral. Embora vários pacotes de software (por exemplo, ALFUN-A, MicroSAMPO, ALFA) existam para abordar esses desafios, a função analítica específica utilizada no amplamente adotado código de computador ALF não havia sido detalhada na literatura antes deste trabalho.

Metodologia
O artigo apresenta uma função analítica original projetada para descrever a forma assimétrica do pico alfa usando apenas quatro parâmetros livres. Esta função foi implementada no código de computador ALF, uma ferramenta baseada em Fortran desenvolvida por volta de 1993 (modificada de um código anterior de raios gama, P.I.M.P.) para avaliação de espectros alfa.

O núcleo da metodologia é a seguinte função gaussiana modificada, $N(x)$, onde $x$ representa o número do canal:

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

Parâmetros:

• $\alpha$: Amplitude (contagens no máximo).
• $M$: Posição do máximo do pico.
• $\beta$ e $\gamma$: Parâmetros de forma (ambos $>0$) específicos da geometria da fonte e do detector.

Comportamento da Função:

• Lado de alta energia ($x > M$): A função converge para uma distribuição gaussiana padrão.
• Lado de baixa energia ($x < M$): A função torna-se mais complexa, transicionando para uma função exponencial ordinária em distâncias significativas do máximo do pico ($x \ll M$).
• Papel dos Parâmetros: $\gamma$ governa a forma do lado direito do pico, enquanto $\beta$ e $\gamma$ governam conjuntamente o lado esquerdo.

Procedimento de Ajuste (Deconvolução):
O código ALF utiliza um procedimento de ajuste de função não linear para deconvolver até 10 picos alfa sobrepostos. O processo é executado em duas etapas distintas:

1. Determinação da Forma Inicial: Uma linha espectral única e bem definida (ou um fragmento contendo o máximo) é ajustada para determinar os quatro parâmetros ($\alpha, M, \beta, \gamma$). Isso utiliza um ajuste de mínimos quadrados não linear padrão via 10 iterações de uma aproximação linear (expansão da série de Taylor).
2. Deconvolução de Picos: Uma vez que $\beta$ e $\gamma$ são determinados, eles são tratados como constantes fixas para todo o espectro, pois dependem da qualidade da fonte e da geometria fonte-detector, e não de picos individuais. A segunda etapa envolve o ajuste das amplitudes ($\alpha_i$) e das posições ($M_i$) para todos os picos sobrepostos simultaneamente. Isso é alcançado expandindo a soma das funções em uma série e resolvendo para os incrementos do vetor de variáveis usando um algoritmo de Cracoviana.

O algoritmo de Cracoviana é destacado como uma vantagem específica, pois calcula as incertezas dos parâmetros ajustados dentro do mesmo procedimento usado para determinar os valores dos parâmetros.

Principais Contribuições e Resultados

• Simplicidade Analítica: O artigo introduz uma função compacta de quatro parâmetros que modela efetivamente a assimetria complexa das linhas alfa sem exigir o ajuste complexo de múltiplos segmentos.
• Validação e Longevidade: A função tem sido de uso rotineiro desde 1993. Nos últimos 30 anos, foi aplicada a mais de 15.000 espectros. A validação foi realizada comparando as áreas de pico calculadas com a soma das contagens, com todos os testes passando positivamente.
• Capacidade de Deconvolução: O método separa com sucesso picos parcialmente sobrepostos, conforme demonstrado com amostras de frações de Th de ossos animais.
• Utilidade Inesperada: Os autores observam que a função provou ser surpreendentemente útil para analisar espectros de radiação beta (especificamente $^{90}\text{Sr}$/$^{90}\text{Y}$) obtidos com um espectrômetro de cintilação líquida usando uma escala de canais logarítmica, alcançando resultados dentro de uma margem de erro de 3%.

Significância
O artigo reivindica significância principalmente através da validação bem-sucedida e de longo prazo desta função de quatro parâmetros específica dentro do software ALF. Ao formalizar a base analítica anteriormente não documentada do ALF, este trabalho fornece uma ferramenta matematicamente definida e robusta para a deconvolução de picos alfa sobrepostos. O método oferece um equilíbrio entre rigor analítico e eficiência computacional, permitindo a determinação precisa das áreas de pico em espectros complexos onde o ajuste gaussiano tradicional falha. O uso do algoritmo de Cracoviana garante ainda que as incertezas dos parâmetros sejam calculadas rigorosamente juntamente com o ajuste.