Amplitude analysis and branching fraction measurement of the decay $D^0 \to K^+K^-\pi^0\pi^0$

Usando dados do detector BESIII, este estudo realiza pela primeira vez uma análise de amplitude do decaimento $D^0 \to K^+K^-\pi^0\pi^0$, determinando sua fração de decaimento absoluta, identificando o processo intermediário dominante como $D^0 \to K^{*}(892)^+K^{*}(892)^-$ (que é dominado por onda S) e medindo sua fração de polarização longitudinal.

O essencial

Imagine que o universo subatômico é uma enorme sala de dança cheia de partículas que nascem, dançam e se transformam em outras coisas. Os cientistas do experimento BESIII (localizado na China) são como fotógrafos e coreógrafos que observam essa sala de dança para entender as regras do movimento.

Este artigo é o relato de uma "sessão de fotografia" muito específica, onde eles estudaram uma dança particular chamada $D^0 \to K^+K^-\pi^0\pi^0$.

Aqui está a explicação simplificada, passo a passo:

1. O Cenário: A Fábrica de Dança

Os cientistas usaram um acelerador de partículas chamado BEPCII, que funciona como uma pista de dança onde elétrons e pósitrons (partículas de luz e antiluz) colidem. Quando eles colidem com a energia certa, eles criam pares de partículas chamadas D0 e anti-D0.

Pense no D0 como um dançarino que, ao final de sua vida, decide se transformar em quatro outros dançarinos menores: dois kaons ($K^+$ e $K^-$) e dois píons neutros ($\pi^0$). O objetivo do estudo foi entender como essa transformação acontece.

2. O Mistério: A Coreografia Oculta

Quando o D0 se transforma, ele não faz isso de uma vez só. Geralmente, ele passa por "intermediários", como se fosse uma coreografia com passos intermediários.

• O problema: É difícil saber quais passos foram dados. Será que o D0 virou dois pares de parceiros que depois se separaram? Ou será que ele virou uma partícula estranha que depois explodiu em quatro?
• A solução (Análise de Amplitude): Os cientistas usaram uma técnica matemática avançada (chamada "análise de amplitude") para reconstruir a coreografia. Eles olharam para milhões de eventos e tentaram descobrir quais "danças intermediárias" (ressonâncias) estavam acontecendo.

3. As Descobertas Principais

A. O Passo Dominante: O "Casal de K*"

A descoberta mais importante foi que a maioria das vezes (cerca de 40% das vezes), o D0 faz um movimento específico: ele se transforma temporariamente em dois parceiros chamados $K^*(892)$ (um positivo e um negativo), que depois se separam nos quatro dançarinos finais.

• Analogia: Imagine que o D0 é um casal que se separa em dois casais intermediários, que então se separam em quatro pessoas sozinhas. Esse foi o movimento mais comum.

B. A Medida da "Polarização" (Quem está no centro?)

Quando duas partículas vetoriais (como os $K^*$) dançam juntas, elas podem girar de diferentes maneiras.

• Imagine dois patinadores girando. Eles podem girar de lado (transversal) ou girar em pé, um em cima do outro (longitudinal).
• Os cientistas mediram que, neste caso, a dança é 47% "vertical" (longitudinal).
• Por que isso importa? Teorias antigas diziam que a dança deveria ser majoritariamente "horizontal" (transversal). O fato de ser quase metade vertical é uma surpresa! Isso sugere que as regras da física que usamos para prever essas danças precisam ser ajustadas ou que há forças ocultas (interações finais) influenciando o movimento.

C. A Frequência da Dança (Razão de Branching)

Os cientistas também contaram quantas vezes essa dança acontece em relação a todas as outras que o D0 pode fazer.

• Eles descobriram que essa dança específica acontece em 0,073% dos casos.
• É uma dança rara, mas não impossível. A precisão dessa medida é muito maior do que em estudos anteriores (o dobro de precisão).

4. O Método: Como eles fizeram isso?

1. Coleta de Dados: Eles usaram dados de 20,3 "femtobarns" (uma unidade de quantidade de dados) coletados pelo detector BESIII. É como ter uma câmera de alta velocidade filmando milhões de colisões.
2. Filtragem: Eles usaram um "filtro" (seleção de eventos) para pegar apenas os casos onde o D0 realmente fez a dança correta, descartando o "ruído" (outros eventos que parecem semelhantes, mas não são).
3. Simulação: Eles criaram milhões de "danças virtuais" no computador para comparar com a realidade e entender onde o detector pode ter cometido erros.
4. Ajuste Fino: Eles ajustaram um modelo matemático complexo até que ele se encaixasse perfeitamente nos dados reais, revelando os passos secretos da coreografia.

5. Por que isso é importante?

A física de partículas é como tentar entender as regras de um jogo olhando apenas as jogadas, sem ver as regras escritas.

• Este estudo mostra que as "regras" que os teóricos escreveram (modelos matemáticos) não estão totalmente corretas para prever como essa dança acontece.
• Ao medir com precisão a "polarização" (como as partículas giram), eles estão testando se há alguma "nova física" ou se precisamos apenas refinar nossa compreensão de como as partículas interagem (interações finais).

Resumo em uma frase

Os cientistas do BESIII usaram um detector gigante para filmar a "dança" rara de uma partícula chamada D0 se transformando em quatro outras, descobrindo que a coreografia é mais complexa e "vertical" do que os teóricos esperavam, o que ajuda a refinar as leis da física que governam o universo subatômico.

Resumo técnico

Título: Análise de Amplitude e Medição da Fração de Ramificação do decaimento $D^0 \to K^+K^-\pi^0\pi^0$

Colaboração: BESIII
Fonte de Dados: Colisões $e^+e^-$ no detector BESIII, energia no centro de massa de 3,773 GeV (ressonância $\psi(3770)$), com luminosidade integrada de 20,3 fb$^{-1}$.

---

1. Problema e Contexto Físico

O estudo de decaimentos hadrônicos de mésons de charme é desafiador devido à massa do quark charm, que é intermediária: não é suficientemente pesada para permitir expansões de massa de quark pesado confiáveis, nem suficientemente leve para a aplicação direta da teoria de perturbação quiral.

• Objetivo Específico: O decaimento $D^0 \to K^+K^-\pi^0\pi^0$ é um decaimento suprimido por Cabibbo (SCS). A compreensão dos mecanismos de decaimento subjacentes, incluindo a dinâmica de interação forte e violação de CP, requer uma análise detalhada dos processos intermediários.
• Gap de Conhecimento: Até este trabalho, não existia uma análise de amplitude completa para este canal específico. Além disso, as previsões teóricas para a fração de polarização longitudinal em decaimentos $D \to VV$ (onde $V$ é um vetor) variam significativamente entre diferentes modelos (ex: Modelo de Quarks Constituintes Não-Relativísticos vs. Simetria SU(3) de Sabor).

2. Metodologia

Seleção de Eventos e Reconstrução

• Técnica de Dupla Marca (Double-Tag - DT): Utilizou-se a técnica de DT para isolar eventos de $D^0\bar{D}^0$. Um méson $D^0$ (lado da marca) é reconstruído em modos hadrônicos simples ($\bar{D}^0 \to K^+\pi^-$, $K^+\pi^-\pi^0$, $K^+\pi^-\pi^-\pi^+$), enquanto o outro ($D^0$) é reconstruído no modo sinal $K^+K^-\pi^0\pi^0$.
• Corte de Seleção: Foram aplicados cortes rigorosos em momento, identificação de partículas (PID), energia de fótons e massa invariante de $\pi^0$.
• Fita Cinemática: Um ajuste cinemático de restrição total foi realizado para melhorar a resolução de massa, restringindo os momentos das partículas finais aos do sistema $e^+e^-$ e as massas reconstruídas aos valores conhecidos do $D^0$ e $\pi^0$.
• Purificação: Após a seleção, 791 eventos foram retidos na região de sinal, com uma pureza de sinal estimada em $(94,2 \pm 0,2)\%$.

Análise de Amplitude

• Modelo de Isobaros: A análise utilizou um modelo de isobaros, onde a amplitude total é a soma coerente de amplitudes de processos intermediários ($M = \sum c_n A_n$).
• Função de Verossimilhança: Foi realizada uma maximização de verossimilhança não-binned, incorporando funções de densidade de probabilidade (PDF) para o sinal e o fundo. O fundo foi modelado usando amostras de Monte Carlo (MC) inclusivas e corrigido com técnicas de XGBoost para eficiência e forma.
• Componentes do Modelo:
  • Fatores de barreira de Blatt-Weisskopf.
  • Propagadores de Breit-Wigner relativísticos para ressonâncias (ex: $K^*(892)$, $f_1(1420)$, $\eta(1475)$).
  • Fatores de spin construídos via formalismo de tensores covariantes.
• Polarização: As amplitudes de helicidade ($H_+, H_0, H_-$) foram construídas independentemente para medir a fração de polarização longitudinal ($F_L$).

Medição da Fração de Ramificação (BF)

• A fração de ramificação absoluta foi calculada combinando o rendimento de eventos de dupla marca com as eficiências de detecção e as frações de ramificação conhecidas dos modos de marca.
• Correções sistemáticas foram aplicadas para eficiência de rastreamento, PID, reconstrução de $\pi^0$ e diferenças entre dados e simulação MC.

3. Principais Contribuições e Resultados

Fração de Ramificação Absoluta

A primeira medição precisa da fração de ramificação absoluta para este canal foi obtida:
$$\mathcal{B}(D^0 \to K^+K^-\pi^0\pi^0) = (0,73 \pm 0,03_{\text{stat}} \pm 0,01_{\text{syst}}) \times 10^{-3}$$
Este resultado é consistente com medições anteriores do BESIII, mas com uma precisão 2,5 vezes maior.

Análise de Amplitude e Processos Dominantes

A análise identificou os seguintes processos intermediários dominantes (com significância estatística $> 5\sigma$):

1. $D^0 \to K^*(892)^+ K^*(892)^-$: O processo dominante, contribuindo com cerca de 40% da fração de ajuste total (considerando interferências).
   • Polarização: A análise revelou que este decaimento é dominado pela onda S.
   • Fração de Polarização Longitudinal ($F_L$): Medida como $F_L = 0,468 \pm 0,046_{\text{stat}} \pm 0,011_{\text{syst}}$.
2. Outros Canais: Contribuições significativas foram observadas para:
   • $D^0 \to \eta(1475)\pi^0$ (com $\eta(1475) \to K^*\bar{K}$).
   • $D^0 \to \eta(1405)\pi^0$ (com $\eta(1405) \to a_0(980)\pi^0$).
   • $D^0 \to K_1(1400)^+ K^-$.
   • $D^0 \to f_1(1420)\pi^0$.

Medição de Polarização

A medição de $F_L \approx 0,47$ para $D^0 \to K^*(892)^+ K^*(892)^-$ é um resultado crucial.

• Comparação Teórica: O valor medido é significativamente maior (mais de 3 desvios padrão) do que as previsões de modelos de quarks constituintes não-relativísticos (NRCQM), que previam dominância de polarização transversal ($F_L \approx 0,31-0,32$).
• Implicação: Isso sugere que os efeitos de interação de estado final (FSI) ou dinâmicas não capturadas pelos modelos simples são importantes.

Frações de Ramificação de Processos Intermediários

Com base na análise de amplitude, foram calculadas as frações de ramificação absolutas para os canais intermediários, sendo o mais proeminente:
$$\mathcal{B}(D^0 \to K^*(892)^+ K^*(892)^-) = (2,61 \pm 0,31_{\text{stat}} \pm 0,15_{\text{syst}}) \times 10^{-3}$$
Este valor é significativamente menor que as previsões de modelos de fatorização e NRCQM, mas alinha-se bem com o modelo de simetria SU(3) de sabor.

4. Significância e Impacto

• Validação de Modelos Teóricos: Os resultados desafiam as previsões de modelos teóricos estabelecidos (NRCQM) sobre a polarização em decaimentos $D \to VV$, indicando a necessidade de refinamento nas teorias de dinâmica hadrônica e interações de estado final.
• Precisão Experimental: Fornece a primeira medição completa de amplitude para este canal, estabelecendo uma base sólida para futuros estudos de violação de CP em decaimentos de mésons D.
• Física de Axial-Vetores: A análise fornece informações únicas sobre a produção de mésons axiais ($K_1$, $f_1$, $\eta$) em decaimentos suprimidos por Cabibbo, ajudando a entender a mistura de estados de axial-vetor.
• Metodologia: O uso de técnicas avançadas de aprendizado de máquina (XGBoost) para modelagem de fundo e eficiência demonstra a evolução das metodologias de análise no experimento BESIII.

Em resumo, este trabalho estabelece novos padrões experimentais para o decaimento $D^0 \to K^+K^-\pi^0\pi^0$, revelando uma dinâmica de polarização inesperada e fornecendo dados cruciais para a compreensão da física de quarks pesados e interações fortes.