A Radioatividade Emitida Por Determinadas Amostras De Substâncias Provém

Como analisar amostras radioativas sem perder horas em medições que não fazem sentido

A análise de radioatividade em amostras de substâncias exige preparação, mas a maioria das pessoas começa pelo instrumento errado ou interpreta os dados de forma ingênua. Vou explicar o processo do jeito que ele funciona na prática, com os detalhes que aparecem só depois que você erra algumas vezes. Isso é o básico, mas é onde muitos cometem o primeiro erro: assumem que a taxa de contagem do seu detector representa diretamente a atividade da amostra. Não representa. O que o detector registra é uma fração minúscula da radiação emitida, determinada pela geometria, pela eficiência do cristal, pela aut absorção dentro da própria amostra e por fatores que muitas vezes nem estão no manual do equipamento. O que você precisa fazer antes de ligar qualquer coisa é montar uma tabela de calibração usando fontes padrões conhecidas. Se você pula essa etapa, cada leitura que fizer terá um erro sistemático que pode variar de 20% a 400%, dependendo da energia da radiação e da matriz da sua amostra. Fontes padrão de Cs-137, Co-60 e Am-241 são o mínimo viável. Use as energias características de cada uma para construir uma curva de eficiência versus energia para o seu detector específico.

Aqui vai uma informação que poucos mencionam: a eficiência do seu detector não cai de forma linear com a energia. Em detectores de centelhamento, por exemplo, a queda é mais acentuada em energias abaixo de 50 keV do que em energias acima de 500 keV. Já em detectores de germânio hiperpuro, o comportamento é diferente porque a probabilidade de interação por efeito fotoelétrico escala com Z elevado a potências altas. Se você está medindo amostras com emissão mista, como urânio natural com seus filhotos, a correção de eficiência mal feita distorce completamente o espectro resultante. Outro detalhe que parece óbvio mas é ignorado com frequência: o tempo de contagem. A estatística de decaimento segue uma distribuição de Poisson, o que significa que a incerteza relativa é proporcional a 1 dividido pela raiz quadrada do número de eventos. Para obter uma incerteza de 1% em uma amostra de baixa atividade, você precisa de 10 mil contagens. Se a taxa de contagem for de 10 cps, isso dá 1000 segundos, pouco mais de 16 minutos, sem contar o fundo. Se for 1 cps, são 166 minutos. Planeje seu tempo de medição com antecedência.

Preparação da amostra e montagem experimental

A forma como você prepara a amostra pode determinar se seus dados são úteis ou se precisam ser descartados. Amostras sólidas precisam ter geometria reprodutível. Espessura uniforme, área de deposição conhecida, posição fixa em relação ao detector. Qualquer variação altera a distância efetiva e a aut absorção, e esses dois fatores não são triviais de corrigir depois. Para amostras líquidas, use vials de polietileno ou plástico fino, preferencialmente da mesma posição e volume. A água absorve radiação beta com significativa intensidade. Se você está medindo puramente emissões beta, a espessura do líquido pode bloquear uma fração considerável das partículas de menor energia. Uma camada de 5 mm de água já reduz significativamente o fluxo de partículas beta abaixo de 100 keV. Medições de fundo são obrigatórias. Sem elas, não há como subtrair corretamente a contribuição ambiental. O fundo varia conforme a localização do laboratório, o nível do solo, a presença de materiais de construção próximos e até a umidade. Um fundo típico em um laboratório urbano pode variar entre 0.2 e 1.5 cps dependendo da região energética do espectro. Meça o fundo durante pelo menos 30 minutos, idealmente na mesma configuração geométrica da amostra, e subtraia o resultado bruto. A incerteza na subtração é a raiz quadrada da soma das incertezas estatísticas de ambas as medições.

Processo de medição e análise espectral

Coloque a amostra na posição definida, registre o tempo exato de início e fim da contagem, e anote também a data. A correção de decaimento é irrelevante para isótopos de meia-vida longa, mas para isótopos como o I-131 com meia-vida de 8 dias ou o Ra-226 cujos filhotos têm meias-vidas curtas, a data importa sim se você quer comparar resultados ao longo de semanas ou meses. Após a aquisição dos dados, o passo seguinte é a análise espectral. No caso de detectores de germânio hiperpuro, você precisará de software de desconvolução espectral. O Genie 2000 da ORTEC, o GammaVision da Canberra e o GammaW são opções comuns. Configure o energy calibration com as linhas conhecidas das fontes padrão, verifique o resolução em FWHM no pico do Cs-137 a 662 keV e confirme se está dentro da especificação do fabricante, tipicamente entre 1.8 e 2.1 keV para um detector de boa qualidade. Defina a região de interesse (ROI) ao redor de cada pico que pretende quantificar. A largura da ROI deve abranger aproximadamente 2 a 3 vezes o FWHM do pico para capturar a maior parte da área, sem incluir ruído excessivo naslaterais. A área do pico é o número total de contagens dentro da ROI, menos a linha de base. A subtração de linha de base é um ponto crítico. Use um método consistente, seja linear ou polinomial, e documente qual escolheu porque diferentes métodos podem dar resultados ligeiramente distintos.

A atividade específica se calcula dividindo a taxa líquida de contagem (cps menos fundo) pela eficiência de detecção na energia correspondente e pelo fator de rendimento de ramo (branching ratio). Para o Cs-137, o branching ratio do gama de 662 keV é 85.1%. Isso significa que apenas 85.1% dos decaimentos produzem esse fóton. Se você ignorar esse fator, superestima a atividade em cerca de 18%. Para o Pb-214 a 352 keV, o branching ratio é 35.9%. A diferença é enorme e afeta diretamente o resultado final.

Problema real que encontrei e como resolvi

Num projeto de caracterização de sedimentos ribeirinhos, comecei a obter valores inconsistentes entre amostras do mesmo local quando comparadas em dias diferentes. As contagens pareciam estáveis, mas a atividade calculada variava de forma inexplicável. O problema era um efeito de deadtime que eu estava ignorando. O detector estava operando próximo de 15% de deadtime sem eu perceber, porque o software exibia apenas o tempo corrido, não o tempo vivo. Após a correção adequada pelo modelo deadtime do sistema, os valores estabilizaram e as variações diárias caíram para dentro da incerteza estatística esperada. A lição foi registrar sempre o parâmetro de deadtime reportado pelo sistema e garantir que ele esteja abaixo de 10% para medições quantitativas confiáveis. Acima disso, a correção introduz mais incerteza do que resolve. Este método tem limitações sérias que você precisa aceitar antes de começar. Detecção de alfa por si só exige preparações especiais, como deposição em disco plano e uso de detector de barreira superficial. A eficiência geométrica para partículas alfa é muito menor que para gama, e a aut absorção em amostras mesmo finas de 0.5 mm pode eliminar grande parte do sinal se a densidade do material for alta. Para medições de trítio ou C-14, a cintilação líquida é praticamente obrigatória, pois a energia máxima do beta é tão baixa que a radiação não consegue sair da amostra sólida nem atravessar a parede de um vial comum até o cristal. Outro ponto: a determinação de U-238 e Th-232 por espectrometria gama indireta depende dos filhotos em equilíbrio. Se a amostra foi quimicamente processada ou contaminou-se com sais lixiviados, o equilíbrio secular pode estar quebrado. Medir apenas o Pb-214 para inferir U-238 ou o Tl-208 para inferir Th-232 em amostras não equilibradasm gera erros sistemáticos grandes. Nesse caso, a solução é usar espectrometria alfa após separação química por iontropeografia ou extração com solvente, ou recorrer a técnicas como ICP-MS, que medem o isótopo diretamente sem depender do equilíbrio radioativo.

Um terceiro problema prático: a interferência espectral. Picos sobrepostos são frequentes. O Bi-214 tem uma linha a 609 keV e outra a 1120 keV. O K-40 emite a 1460 keV. Se sua amostra contém tanto K-40 quanto uma trilha de radiação próxima, a área do pico fica comprometida. O software tenta resolver com ajustes de função de picos, mas quando a sobreposição é maior que 50% da largura do pico, a incerteza na área aumenta rapidamente. Nesses casos, use outra linha energética do mesmo nuclídeo para confirmação cruzada ou isole quimicamente o elemento de interesse antes da medição.

Cálculo de limite de detecção

Antes mesmo de medir sua primeira amostra desconhecida, calcule o limite mínimo de detecção (MDA) para cada nuclidio que pretende detectar. A fórmula de Currie, amplamente adotada na área, define o MDA como a atividade mínima que pode ser distinguida do fundo com um nível de confiança estatística razoável. Ele depende da taxa de fundo, do tempo de contagem, da eficiência e do branching ratio. Se o MDA da sua configuração é 5 Bq/kg para o Pb-214, qualquer valor abaixo disso nos seus resultados não é confiável, mesmo que o pico apareça no espectro. Relatórios formais de dosimetria ou caracterização ambiental exigem a apresentação dos MDA calculados junto com os resultados, mesmo quando as amostras estão abaixo do limite. Isso evita interpretações equivocadas. Relatar zero ou menos que LDM sem explicitar o valor do limite é impróprio e frequentemente rejeitado em revisões técnicas.

Referências para software e fontes

Para quem precisa de ferramentas de análise, além dos pacotes comerciais citados, existem alternativas gratuitas ou de código aberto que merecem menção. O GammaW é gratuito e bastante capaz para análises rotineiras. O QGDA da Quantus é uma opção sólida sem custo para instituições acadêmicas. Para planejamento de medições e cálculo de MDA, a calculadora disponível no site do NIST oferece uma abordagem prática. Para normas técnicas, a ISO 11929-1 estabelece os critérios para a declaração de valores de atividades e a forma correta de comunicar incertezas e limites de decisão, e a ASTM E244 cobre procedimentos específicos para espectrometria gama de alto resolução com germânio. O acompanhamento dos dados e a documentação do processo completo, desde a preparação da amostra até o cálculo final da atividade, são tão importantes quanto a qualidade do hardware utilizado. Cada decisão tomada durante a medição impacta a incerteza final, e registros consistentes permitem que outro técnico reprova o procedimento e obtenha resultados compatíveis. Quando isso não acontece, os dados ficam presos ao conhecimento de uma única pessoa e perdem validade ao longo do tempo.