Realization of all-to-all fermion propagator for the first principle high accuracy strong interaction prediction

Problema

A Cromodinâmica Quântica em Rede (LQCD) é uma ferramenta primária para calcular a estrutura de hádrons a partir de primeiros princípios, como o acoplamento axial vetorial do núcleon ($g_A$). No entanto, alcançar precisão sub-percentual exige enfrentar dois grandes desafios:

1. Custo Computacional de Propagadores All-to-All: O cálculo de funções de correlação de $N$-pontos (onde $N \ge 4$), que são essenciais para correções de QED, quebra de simetria de isospin (ISB) e loops de quarks desconectados, requer propagadores de férmions all-to-all. O cálculo exato escala como $O(V^2)$ com o volume da rede $V$, tornando-se proibitivamente caro.
2. Contaminação por Estados Excitados (ESC): A extração de elementos de matriz do estado fundamental é complicada pela contaminação de estados excitados (ex: estados do tipo Roper, estados multi-partículas como $\pi N$). Embora o método de destilação (projeção sobre modos de baixo momento) suprima efetivamente a ESC para espectros de hádrons, ele descarta modos de alto momento essenciais para inserções de correntes locais, limitando sua utilidade para $N$-pontos e elementos de matriz.

Técnicas existentes, como fontes sequenciais, abordagens estocásticas com substituição de modo baixo (LMS) e esparsificação de campo, oferecem soluções parciais, mas lutam para lidar simultaneamente com $N$-pontos de alta precisão, diagramas desconectados e supressão de ESC sem incorrer em custos computacionais massivos.

Metodologia: O Algoritmo de "Blending" (Mistura)

Os autores propõem um novo método estocástico chamado algoritmo de "blending" para aproximar propagadores de férmions all-to-all de forma eficiente e sem viés.

• Construção do Espaço de Blending: O método constrói um "espaço de blending" para uma fatia de tempo fixa combinando dois subespaços:
  1. Modos de Baixa Frequência (LFM): O subespaço $L_1$ gerado pelos $N_e$ autovetores de menor valor do Laplaciano espacial covariante de gauge (semelhante à destilação). Isso captura a física de longa distância e permite a construção de estados de hádrons espacialmente estendidos.
  2. Modos de Alta Frequência (SHFM): O complemento ortogonal $L_2$ é aproximado estocasticamente usando $N_{st}$ vetores aleatórios ($\eta_j$) amostrados uniformemente de $L_2$.
• Estimativa Sem Viés: O operador identidade no espaço total é aproximado por uma soma ponderada de projetores sobre estes vetores de base.
  • Os projetores de modo baixo retêm seu peso exato (1).
  • Os projetores estocásticos de modo alto recebem um peso de $\dim(L_2)/N_{st}$ para garantir uma estimativa sem viés do operador identidade no subespaço complementar.
  • Para operadores envolvendo produtos de campos (ex: em funções de $N$-pontos), um tensor de peso $\Omega^{(2)}_{ij}$ é aplicado para corrigir a contagem dupla e garantir que o estimador permaneça sem viés, particularmente para termos estocásticos diagonais.
• Compressão: O propagador de quark all-to-all $S$ é comprimido projetando-o neste espaço de blending. Isso reduz os requisitos de armazenamento de $O((4N_c N_L^3 N_T)^2)$ para $O((4N_T(N_e + N_{st}))^2)$, permitindo o cálculo flexível de funções de correlação com arbitrárias combinações de operadores de interpolação e inserções de corrente.
• Controle de Estado Excitado: O método facilita o uso de análise de Problema de Autovalor Generalizado (GEVP) com uma base enriquecida de interpoladores (incluindo operadores "envolvidos com a corrente", como estados de espalhamento $N f_1$) sem custos adicionais de inversão, pois os propagadores comprimidos podem ser reutilizados para múltiplas combinações de operadores.

Principais Contribções e Resultados

Os autores demonstram a eficiência e a precisão do método de blending através de aplicações em dois conjuntos de CLQCD: C24P29 ($m_\pi \approx 292$ MeV, $a \approx 0.105$ fm) e F64P13 ($m_\pi \approx 134$ MeV, $a \approx 0.077$ fm).

1. Validação da Ausência de Viés:

   • O método foi testado usando o elemento de matriz da corrente vetorial conservada do píon ($\pi^+$). Os resultados mostraram que o uso apenas de modos baixos ($N_{st}=0$) produziu um resultado enviesado ($\approx 0.64$), enquanto a inclusão de modos estocásticos ($N_{st} \ge 50$) restaurou o resultado para a expectativa teórica de 1.000(29).
   • Verificações de consistência foram realizadas para o fator de forma elétrico do píon ($f_{\pi\pi}(Q^2)$) e raio de carga usando tanto funções de correlação de 3-pontos quanto de 4-pontos. Os resultados foram consistentes dentro das incertezas, validando o cálculo da função de 4-pontos.

2. Acoplamento Axial do Nucleon ($g_A$):

   • Isovetor ($g_A^{u-d}$): Usando um ajuste conjunto de 3 estados através de vários números de modos baixos ($\bar{N}_e$) e separações fonte-sumidouro, os autores obtiveram $g_A^{u-d} = 1.2337(84)$ na massa física do píon. Este resultado elimina a ambiguidade associada à seleção de um $\bar{N}_e$ específico e mostra excelente concordância com outros estudos de alta precisão.
   • Decomposição de Singlet e de Sabor: O método permitiu o cálculo de loops de quarks desconectados para a carga axial de singlet ($g_A^{u+d+s}$) e contribuições individuais de sabor.
     • Usando uma análise GEVP com uma base incluindo interpoladores $N$, $N f_1$ e $N f'_1$, eles suprimiram efetivamente a contaminação por estados excitados.
     • Os resultados finais no ponto físico ($a=0.077$ fm) são:
       • $g_A^u = 0.8408(86)$
       • $g_A^d = -0.3929(86)$
       • $g_A^s = -0.0381(57)$
       • $g_A^{u+d+s} = 0.410(20)$

3. Eficiência Computacional:

   • O método de blending reduz significativamente o número de inversões do operador de Dirac necessárias em comparação com métodos tradicionais de fonte sequencial ou estocástica.
   • Para o $g_A$ isovetor, o custo foi reduzido em aproximadamente uma ordem de magnitude comparado a estudos recentes de alta precisão (ex: ETM, RQCD, PNDME) visando precisão semelhante.
   • Para o $g_A$ de singlet (envolvendo loops desconectados), o método alcançou uma precisão de 4% com 0,70 milhões de inversões, comparado a 13,7–34,1 milhões de inversões em trabalhos anteriores comparáveis. Esta eficiência advém da capacidade de reutilizar o propagador comprimido para múltiplos interpoladores e inserções de corrente sem novas inversões.

Significância e Alegações

O artigo afirma que o método de blending fornece um arcabouço robusto para cálculos de alta precisão de estrutura de hádrons em LQCD, abordando especificamente os duplos desafios do custo computacional e da contaminação por estados excitados.

• Controle Sistemático: Ao permitir a construção de funções de correlação com múltiplos operadores de interpolação compartilhando os mesmos números quânticos, o método facilita a extração robusta de elementos de matriz de estado fundamental via análise variacional (GEVP), suprimindo efetivamente a ESC.
• Escalabilidade: Na massa física do píon, o método torna previsíveis as predições de precisão sub-percentual para elementos de matriz do núcleon com recursos computacionais significativamente reduzidos.
• Versatilidade: A abordagem é aplicável a várias funções de $N$-pontos, incluindo diagramas desconectados, correções de QED e efeitos de ISB, que são cruciais para testar o Modelo Padrão.
• Limitações: Os autores observam que o custo computacional permanece proibitivo para férmions Overlap devido aos seus maiores custos de inversão, sugerindo a necessidade de otimizações adicionais para uma aplicabilidade mais ampla.

Em resumo, o método de blending representa um passo significativo para tornar viáveis cálculos de primeira princípio de observáveis hadrônicos complexos de alta precisão, permitindo futuros estudos de espalhamento de múltiplos corpos e física de partons com incertezas estatísticas e sistemáticas reduzidas.