Optimized filtering for pulse-shape based pile-up rejection applied to $0\nu\beta\beta$ search with $^{100}$Mo

Este artigo apresenta um algoritmo de filtragem digital otimizado para detectores criogênicos que reduz significativamente o fundo induzido por pile-up em buscas de dupla beta sem neutrinos em 31% em comparação com métodos anteriores, mantendo uma alta eficiência de sinal.

O essencial

O Panorama Geral: Ouvindo um Sussurro em uma Tempestade

Imagine que você está tentando ouvir um sussurro específico (um evento nuclear raro chamado decaimento duplo beta sem neutrinos) em uma sala muito barulhenta. O problema é que a sala está cheia de pessoas conversando constantemente (ruído de fundo).

Pior ainda, às vezes duas pessoas falam quase ao mesmo tempo. As vozes delas se sobrepõem, criando um som que parece exatamente com o sussurro especial que você está procurando. Na física, essa sobreposição é chamada de "pile-up" (empilhamento).

Se você não conseguir distinguir o sussurro real do falatório sobreposto, terá alarmes falsos e seu experimento falhará. Este artigo apresenta uma nova maneira mais inteligente de filtrar esse falatório sobreposto para que os cientistas possam ouvir o sinal real.

O Problema: A Confusão do "Duplo Aplauso"

Os cientistas estão usando um tipo especial de detector feito de cristais (especificamente Molibdênio-100). Esses cristais são como microfones super sensíveis que aquecem ligeiramente quando uma partícula os atinge.

• O Sinal Real: Uma única partícula atinge o cristal, criando um "pulso" de calor limpo e suave.
• O Pile-up: Duas partículas atingem o cristal quase instantaneamente (em menos de um milissegundo). O detector vê um grande "pulso" bagunçado que parece um único evento, mas que é, na verdade, dois eventos esmagados um contra o outro.

Como o "pulso bagunçado" se parece tanto com o "sinal real", os métodos antigos de filtragem não eram bons o suficiente. Eles estavam descartando muitos sinais reais apenas para capturar alguns sinais falsos.

A Solução: Um "Engenheiro de Som" Inteligente

Os autores criaram um novo filtro digital (uma ferramenta matemática) que atua como um engenferieiro de som super inteligente. Em vez de apenas ouvir o volume, este engenheiro analisa a forma e a textura da onda sonora.

Veja como o novo método deles funciona, usando uma analogia:

O Truque dos "Dois Filtros"
Imagine que você está tentando distinguir entre uma única batida de tambor e duas batidas de tambor que aconteceram tão rápido que pareceram uma só.

1. O Filtro A é sintonizado para ouvir a "suavidade" do som. Ele adora a batida de tambor única e limpa.
2. O Filtro B é sintonizado para ouvir a "aspereza" ou as bordas afiadas que ocorrem quando duas batidas se sobrepõem.

O novo algoritmo não olha apenas para o resultado do Filtro A ou do Filtro B isoladamente. Em vez disso, ele pega a razão (a comparação) dos dois.

• Se o som for uma batida única, a razão permanece estável.
• Se o som for uma sobreposição bagunçada, a razão muda dramaticamente.

Ao comparar constantemente essas duas perspectivas, o algoritmo consegue detectar os eventos sobrepostos "falsos" com muito mais precisão do que antes.

Como Eles Construíram Isso

Os cientistas não apenas adivinharam as configurações para este filtro. Eles usaram um computador para "treiná-lo", de forma semelhante a como você treinaria um cachorro.

• Eles alimentaram o computador com milhares de exemplos de sinais reais e sinais sobrepostos falsos.
• O computador ajustou seus botões internos (pesos matemáticos) milhões de vezes por segundo para encontrar a configuração perfeita que captura o máximo de falsos enquanto deixa passar o máximo de sinais reais.
• Eles usaram uma técnica poderosa chamada descida do gradiente (um método de descer lentamente uma colina para encontrar o ponto mais baixo) para encontrar as melhores configurações de filtro possíveis.

Os Resultados: Um Sinal Mais Claro

Quando testaram este novo filtro nos dados do experimento CROSS (uma fase de teste para o futuro experimento CUPID):

• O Objetivo: Eles queriam manter 90% dos sinais reais (eficiência).
• A Vitória: Com o novo filtro, eles foram capazes de rejeitar 31% a mais dos sinais sobrepostos falsos em comparação ao método antigo.

Pense nisso da seguinte forma: Se o método antigo capturasse 100 sinais falsos, o novo método captura 131, sem perder nenhum dos sussurros reais. Este é um grande avanço para um experimento onde cada alarme falso conta.

Por Que Isso Importa

O artigo foca especificamente em detectores criogênicos (detectores que devem ser mantidos próximos ao zero absoluto) que buscam o decaimento duplo beta sem neutrinos.

• Os autores afirmam que este método é amplamente aplicável a outras tecnologias de detectores, mas seus resultados e números específicos referem-se estritamente à melhoria da busca por este decaimento nuclear específico usando cristais de Molibdênio.
• Eles não alegaram que isso funciona para imagens médicas ou outros campos neste artigo; o foco é inteiramente em tornar os experimentos de física mais sensíveis.

Resumo

O artigo introduz uma nova "orelha" matemática para os cientistas. Ao comparar duas maneiras diferentes de ouvir um sinal, esta nova ferramenta consegue detectar quando dois eventos se fundiram acidentalmente em um só. Isso permite que eles ignorem o ruído e foem nos sinais raros e preciosos que podem desbloquear segredos sobre o universo, especificamente em relação à natureza dos neutrinos.

Resumo técnico

Resumo Técnico: Filtragem Otimizada para Rejeição de Pile-up Baseada na Forma do Pulso em Pesquisas de $0\nu\beta\beta$ com $^{100}$Mo

Definição do Problema
Em buscas de eventos raros, como a busca pela dupla desintegração beta sem neutrinos ($0\nu\beta\beta$) utilizando calorímetros criogênicos, os eventos de pile-up — onde sinais distintos se sobrepõem temporalmente — representam um desafio crítico. Para detectores que utilizam $^{100}$Mo, como o experimento CUPID, a fonte dominante de fundo espera-se que surja do pile-up de eventos de dupla desintegração beta com dois neutrinos ($2\nu\beta\beta$). Esses eventos podem mimetizar o sinal monoenergético de $0\nu\beta\beta$ em 3034 keV. Melhorias de hardware (por exemplo, sensores mais rápidos como TESs e MMCs) visam reduzir a duração dos pulsos, mas o pile-up residual permanece inevitável em experimentos de grande escala. As técnicas de análise existentes, incluindo a filtragem ótima padrão e a rejeição baseada em $\chi^2$, perdem poder de discriminação quando a separação temporal entre os pulsos é menor que o tempo de subida do sinal (especificamente $\Delta t < 0,8$ ms para os detectores estudados). Isso torna necessária a utilização de métodos avançados de discriminação baseados em software para mitigar o fundo sem sacrificar a eficiência do sinal.

Metodologia
Os autores propõem um arcabouço para derivar um filtro digital otimizado para discriminação de pile-up baseado em uma descrição do sinal no domínio da frequência. O método assume:

1. Uma resposta ao impulso $s(t)$ do detector totalmente conhecida e independente da energia.
2. Ruído gaussiano estacionário de média zero, caracterizado por uma densidade espectral de potência (PSD) $P_N(\omega)$.
3. Eventos de pile-up modelados como a soma de dois pulsos com razão de amplitude relativa $r$ e separação temporal $\Delta t$.

O núcleo do método envolve a construção de um parâmetro discriminante $Y$ baseado na razão de dois estimadores de amplitude, $X[m_1]$ e $X[m_2]$, derivados do sinal medido $\hat{x}(\omega)$:
$$Y = \frac{X[m_1]}{X[m_2]}$$
onde $X[m_i]$ é o máximo da parte real da transformada de Fourier inversa do sinal filtrado por uma função de transferência $H(\omega)$ e uma função de ponderação $m_i(\omega)$. O filtro ótimo padrão $H(\omega)$ maximiza a relação sinal-ruído (SNR), enquanto as funções de ponderação $m_1(\omega)$ e $m_2(\omega)$ são as variáveis a serem otimizadas.

O objetivo da otimização é minimizar a taxa global de identificação incorreta de pile-up $\langle M \rangle$, definida como a probabilidade de um evento de pile-up ser classificado incorretamente como um pulso único, sujeita a uma eficiência de aceitação de pulso único fixa $\epsilon$. Os autores utilizam PyTorch e diferenciação automática para realizar o gradiente descendente nas funções de ponderação $m_1(\omega)$ e $m_2(\omega)$, tratando-as como parâmetros treináveis discretizados sobre bins de frequência. A função de custo é avaliada sobre uma grade de separações de tempo e razões de amplitude, ponderadas pelas distribuições esperadas de eventos $2\nu\beta\beta$.

Principais Contribuições

• Arcabouço de Otimização: Desenvolvimento de um arcabouço geral e quantitativo para modelar, identificar e rejeitar eventos de pile-up em calorímetros criogênicos utilizando filtros otimizados no domínio da frequência.
• Discriminante Baseado em Razão: Implementação de um parâmetro discriminante baseado em razão, que aproveita a correlação entre dois observáveis filtrados para reduzir a sensibilidade a flutuações de amplitude e melhorar o poder de discriminação em comparação com métodos de amplitude única.
• Implementação Computacional: Uma implementação modular em Python utilizando PyTorch que permite a filtragem eficiente via aceleração por GPU, capaz de lidar com grandes grades de parâmetros.
• Validação Sistemática: Um estudo abrangente das incertezas sistemáticas decorrentes da seleção do template do pulso, tamanho da janela de processamento, posição do pulso e discretização da grade.

Resultos
O método foi validado utilizando dados do experimento CROSS, que operou uma torre de cristais de $Li_2^{100}MoO_4$ centelhantes acoplados a detectores de luz (LDs) com amplificação Neganov–Trofimov–Luke (NTL).

• Melhoria de Desempenho: Comparado a um método de referência (apresentado anteriormente em Eur. Phys. J. C 83(5), 373 (2023)), o novo filtro otimizado reduz o índice de fundo induzido por pile-up em uma média de 31% com 90% de eficiência de pulso único.
• Índice de Fundo: O índice de fundo residual médio extrapolado usando o novo método é de $5,9 \times 10^{-5}$ cts/keV/kg/ano (mediana: $5,7 \times 10^{-5}$), aproximando-se da meta de design do CUPID de $5,0 \times 10^{-5}$ cts/keV/kg/ano.
• Características do Filtro: Os filtros otimizados $m_1(\omega)$ e $m_2(\omega)$ exibem uma ponderação de frequência complementar. $m_1(\omega)$ mostra uma resposta crescente na faixa de 1–3 kHz para amplificar características de alta frequência essenciais para resolver o pile-up, enquanto $m_2(\omega)$ minimiza o ganho nesta região para enfatizar a forma geral do sinal.
• Sistemáticas: O estudo identificou que a escolha do template do pulso impacta significativamente os resultados (variação de até 13% no BI). Os autores recomendam o uso de um template de "pulso médio" como uma opção conservadora para evitar artefatos de impressões de ruído de alta frequência em pulsos únicos de alta energia.
• Dependência do Detector: Um estudo preliminar sobre pontos de operação de detectores (corrente de polarização) sugere que a otimização individual do ponto de operação pode melhorar ainda mais a rejeição de pile-up, com variações na rejeição de fundo atingindo um fator de dois para diferentes detectores.

Significância e Alegações
O artigo afirma que o proposto discriminante baseado em razão captura as características dominantes de pile-up relevantes para calorímetros criogênicos com respostas térmicas lentas. Ele fornece uma ferramenta robusta e computacionalmente eficiente para otimizar técnicas de análise na ausência de mudanças de hardware. Os autores afirmam que, embora o método dependa de suposições de ruído e formas de pulso estacionárias, ele demonstra redução de fundo significativa e robusta através de diversos detectores.

Este trabalho posiciona esta técnica como um "controle de software" complementar aos desenvolvimentos de hardware (como sensores mais rápidos ou melhor cobertura de amplificação NTL). Os autores observam que a extrapolação dos resultados para o experimento CUPID completo depende do ambiente de ruído, particularmente na ordem de O(1) kHz, e que futuros upgrades à infraestrutura criogênica do CUPID e ao desempenho dos detectores de luz devem render novas melhorias. O arcabouço é apresentado como amplamente aplicável a diversas tecnologias de detectores e contextos experimentais além da pesquisa específica de $^{100}$Mo.