CHESS: CHEbyshev pSeudo-Spectral transport for Feynman integral differential equations

Resumo Técnico: Chess (CHEbyshev pSeudo-Spectral transport for Feynman integral differential equations)

Enunciado do Problema
A avaliação de integrais mestres de Feynman de múltiplos loops é um gargalo computacional na física de colisores de precisão. À medida que os cálculos avançam para ordens de loop e multiplicidades mais elevadas, as famílias de integrais envolvem inúmeras escalas cinemáticas, grandes conjuntos de integrais mestres e estruturas de singularidade complexas. Embora existam várias estratégias (ex: decomposição de setores, fluxo de massa auxiliar, expansões de série local), um subproblema numérico específico permanece desafiador: o transporte de alta precisão unidimensional de equações diferenciais fatorizadas em $\epsilon$ ao longo de um caminho prescrito no espaço cinemático.

Os métodos de propagação local padrão frequentemente exigem muitos segmentos de expansão ou passos pequenos para controlar a alta precisão e o comportamento nos pontos finais, particularmente quando a matriz de coeficientes puxada envolve combinações esparsas de matrizes constantes e pullbacks escalares de uma forma, ou quando o sistema inclui ingredientes algébricos ou elípticos. O desafio reside em resolver o sistema de equações diferenciais ordinárias (EDO) resultante de forma eficiente, mantendo a acurácia espectral e lidando com potenciais singularidades regulares espúrias nos pontos finais da integração.

Metodologia
O artigo introduz o Chess (CHEbyshev pSeudo-Spectral transport), um pacote da linguagem Wolfram projetado para resolver o problema de transporte unidimensional:
$$\frac{d}{dt}J(t, \epsilon) = \epsilon B(t)J(t, \epsilon)$$
onde $J$ é um vetor de integrais mestres (ou seus coeficientes de expansão em $\epsilon$) e $B(t)$ é a matriz de coeficientes puxada.

A metodologia central baseia-se na colocação espectral de Chebyshev–Lobatto. O método explora a analiticidade da conexão puxada no parâmetro do caminho complexificado. Se a matriz de coeficientes $B(t)$ for analítica em uma vizinhança complexa do intervalo de transporte, a solução converge geometricamente com o número de nós de colocação.

Principais características algorítmicas incluem:

1. Hierarquia $\epsilon$-Fatorizada: O sistema é expandido em uma janela Laurent finita. Devido à forma fatorizada em $\epsilon$, a hierarquia de coeficientes é triangular. Cada coeficiente $J^{(r)}$ é obtido a partir do anterior $J^{(r-1)}$ resolvendo um sistema de primeira ordem não homogêneo. Isso permite que o solver utilize um único nó escalar fatorizado atuando sobre muitos lados direitos, em vez de construir um sistema monolítico para todos os nós e integrais simultaneamente.
2. Modos de Entrada Duplos:
   • Modo de Dados Preparados: O usuário fornece matrizes esparsas constantes $A_\alpha$ e pullbacks escalares pré-computados $b_\alpha(t)$. O pacote monta $B(t) = \sum A_\alpha b_\alpha(t)$ nos nós de colocação. Isso separa as avaliações de escala dependentes do caminho dos dados de matriz constante, permitindo cache e montagem esparsa.
   • Modo de Matriz Direta: O usuário fornece um avaliador numérico para a matriz completa $B(t)$. Isso permite que o solver lide com sistemas contendo períodos elípticos ou outras funções especiais sem exigir uma decomposição logarítmica específica.
3. Regularização de Pontos Finais: Para lidar com singularidades regulares espúrias nos pontos finais (onde $B(t)$ pode divergir, mas a solução permanece finita), o Chess implementa uma regularização baseada em resíduos.
   • Para um ponto final à direita, a linha da equação diferencial é substituída por uma condição de fonte finita derivada do resíduo e da extração de parte finita da conexão.
   • Para um ponto final à esquerda, o valor de contorno conhecido é removido, e um "levantamento" (lifted) desconhecido (o quociente de diferença de um lado regular) é resolvido, substituindo a equação diferencial singular por uma condição de derivada finita.
4. Paralelismo: O paralelismo é aplicado apenas a tarefas independentes, como a avaliação de pullbacks escalares nos nós ou a extração de dados de ponto final. A recursão sequencial em $\epsilon$ é preservada, pois a camada $r$ depende da camada $r-1$.

Contribuições Principais

• Eficiência Algorítmica: Ao utilizar a estrutura triangular de sistemas fatorizados em $\epsilon$ e a montagem de matrizes esparsas, o Chess evita o overhead de memória e computação de solvers monolíticos. Ele alcança uma convergência rápida de nós, frequentemente atingindo os limites de precisão da máquina com relativamente poucos nós de colocação.
• Flexibilidade: O pacote não é restrito a sistemas logarítmicos. A interface de matriz direta permite transportar sistemas que contêm ingredientes elípticos ou modulares, desde que um avaliador numérico confiável para $B(t)$ esteja disponível.
• Robustez: A implementação inclui tratamento específico para pontos finais de singularidade regular espúria, um problema comum em transportes cinemáticos que pode desestabilizar métodos de colocação padrão.
• Implementação: O pacote é escrito na linguagem Wolfram, oferecendo uma interface de alto nível para usuários que já realizaram a redução de integrais, construção de base e determinação de constantes de contorno.

Resultados e Benchmarks
O artigo valida o Chess através de três grupos de benchmarks:

1. Duplo Pentágono de Seis Pontos de Dois Loops: Transportado na região física e comparado contra dados independentes de AMFlow. Com 24 nós, a discrepância foi de $1.37 \times 10^{-31}$, diminuindo para $10^{-101}$ com 96 nós, saturando na precisão dos dados de referência.
2. Famílias de Cinco Pontos de Três Loops: Quatro famílias planas sem massa (PBB, BPB, BHB, PBP) foram testadas com regularização de ponto final à esquerda. O método demonstrou refinamento estável de nós com diferenças em torno de $10^{-70}$ (aprox. 70 dígitos). Uma comparação de protocolo fixo com o código de série local DiffExp mostrou o Chess sendo 8,8 a 17,9 vezes mais rápido em tempo de execução (wall time). Para o maior sistema (734 integrais mestres), o Chess também demonstrou menor uso de memória da árvore de processos em comparação ao DiffExp.
3. Triângulo Não-Planar (Caso Elíptico): Uma família de triângulo não-planar de dois loops com períodos elípticos foi transportada usando o modo de matriz direta. O método alcançou convergência geométrica até o limite de precisão de 200 dígitos da referência independente AMFlow, atingindo erros abaixo de $10^{-200}$ com 512 nós.

Significância e Escopo
O artigo posiciona o Chess como uma ferramenta especializada de alto desempenho dentro do ecossistema mais amplo de computação de integrais de Feynman. Ele afirma explicitamente que o pacote não realiza redução de integrais, constrói bases canônicas, determina constantes de contorno ou seleciona caminhos analíticos. Essas etapas preparatórias são assumidas como concluídas pelo usuário.

A significância do Chess reside em sua capacidade de realizar o passo final de transporte, que é computacionalmente intensivo, com alta precisão e eficiência. Ele aproveita as propriedades analíticas do parâmetro do caminho para alcançar convergência espectral, tornando-o particularmente eficaz para sistemas de múltiplas escalas e grandes, onde métodos tradicionais de série local podem exigir segmentação excessiva. Os autores enfatizam que o sucesso do método depende da analiticidade das funções de coeficiente fornecidas; não é um algoritmo automático de planejamento de caminho para atravessar singularidades físicas. O pacote serve como uma camada de transporte robusta para usuários que já estabeleceram o framework analítico necessário.