Highly improved staggered quarks on anisotropic lattices

Resumo Técnico: Quarks Staggered Altamente Melhorados em Lattices Anisotrópicas

Problema
A motivação primária deste trabalho é o desenvolvimento de um formalismo para gerar ensembles de campos de gauge otimizados para reconstrução espectral em QCD de rede (lattice QCD), especificamente para o estudo dos padrões de fusão de quarkonia pesados e estados de sabor aberto pesado no plasma de quarks-glúons. A reconstrução espectral a partir de funções de correlação euclidianas é um problema mal posto que requer ensembles com estatísticas muito altas e extensões temporais grandes ($N_\tau$). Enquanto ensembles isotrópicos com grande $N_\tau$ são computacionalmente proibitivos, lattices anisotrópicas (onde o espaçamento de rede espacial $a_\sigma$ é significativamente maior que o espaçamento temporal $a_\tau$, definido pela anisotropia renormalizada $\xi = a_\sigma/a_\tau \sim 6-8$) oferecem uma solução.

Simulações anisotrópicas de grande escala existentes tipicamente empregam férmions de Wilson, que são computacionalmente caros e limitam a precisão estatística. Por outro vez, embora os férmions staggered sejam eficientes para QCD de temperatura finita, tentativas anteriores usando férmions staggered anisotrópicos naïve (não melhorados) foram abandonadas devido a grandes efeitos de discretização, particularmente a quebra de sabor (taste-breaking). Este artigo aborda a lacuna desenvolvendo e ajustando o parâmetro da ação aHISQ (Highly Improved Staggered Quark) anisotrópica, que visa combinar a eficiência computacional dos férmions staggered com os erros de discretização reduzidos de ações melhoradas.

Metodologia
Os autores apresentam um estudo abrangente do ajuste dos parâmetros tanto para o setor de gauge quanto para o de férmions da ação aHISQ em ensembles de gauge puro (quenched).

1. Ajuste de Anisotropia de Gauge:

   • O estudo utiliza a ação de gauge Symanzik-melhorada ao nível da árvore com um parâmetro de anisotropia de gauge bruto $\xi_0$.
   • O ajuste é realizado utilizando o método de fluxo de gradiente (gradient flow). Os autores comparam múltiplos esquemas de fluxo (Wilson, Symanzik, Zeuthen) e discretizações de observáveis (Wilson, Symanzik, Clover, Clover Melhorado).
   • A anisotropia renormalizada $\xi$ é determinada impondo a condição $w_{0,\sigma}/w_{0,\tau} = 1$, onde $w_0$ é a escala de fluxo de gradiente.
   • Os autores identificam o esquema "Wilson flow com observável Clover" (WC) como ideal devido ao seu comportamento monotônico e estabilidade ao atingir o limite contínuo.

2. Ajuste de Anisotropia de Férmion:

   • O operador de Dirac aHISQ é construído introduzindo a anisotropia de férmion bruta $\xi^f_0$ no nível mais externo do smearing (termos Fat7 e Naik) para preservar a supressão do acoplamento com glúons nos cantos da zona de Brillouin.
   • O ajuste envolve o mapeamento dos parâmetros brutos $(a_\sigma m_s, \xi^f_0)$ para a anisotropia renormalizada $\xi^f$ e a massa física do quark strange.
   • Dois métodos para determinar $\xi^f$ são discutidos: a razão das energias do estado fundamental nas direções temporal vs. espacial e o método da relação de dispersão. Os autores preferem o método da relação de dispersão, ajustando a relação energia-momento $a_\tau^2 E^2(p) = a_\tau^2 E^2(0) + a_\sigma^2 p^2 / (\xi^f)^2$.
   • Um algoritmo completo para a força de férmion aHISQ é derivado e implementado para o algoritmo de Monte Carlo de Hibridismo Racional (RHMC), utilizando uma regra de cadeia recursiva para lidar com os fatores de anisotropia.

3. Análise do Espectro:

   • O estudo mede os desdobramentos de massa de sabor de píon (pion taste mass splittings) para as ações naïve e aHISQ através de uma gama de anisotropias ($\xi = 1$ a $8$).
   • Os autores analisam a dependência dos desdobramentos de sabor em relação a $\xi$ e comparam o comportamento das duas ações.

Principais Contribuições e Resultados

• Ajuste de Gauge: Os autores mapearam com sucesso as Linhas de Anisotropia Renormalizada Constante (LCRA) para $\xi=2$ e $\xi=4$ no plano $(\beta, \xi_0)$. Eles reproduziram resultados existentes para $\xi=2$ e estenderam o ajuste para anisotropias mais altas ($\xi=8$) usando o esquema WC.
• Eficiência do Ajuste de Férmion: Um achado significativo é o comportamento da ação aHISQ comparada à ação naïve. Para férmions staggered naïve, a anisotropia de férmion bruta $\xi^f_0$ depende fortemente da massa do quark, exigindo um escaneamento bidimensional no espaço de parâmetros $(a_\sigma m_s, \xi^f_0)$. Em contraste, para aHISQ, a dependência de $\xi^f_0$ em relação à massa do quark é desprezível (plana) através de uma ampla gama de massas. Isso permite que o ajuste da aHISQ seja reduzido a dois escaneamentos unidimensionais, simplificando drasticamente o custo computacional para futuras simulações dinâmicas.
• Comportamento do Espectro de Sabor:
  • Ação Naïve: À medida que a anisotropia aumenta, sabores de píon mais pesados descem em massa, aproximando-se de um novo padrão de degenerescência esperado para anisotropia infinita.
  • Ação aHISQ: O comportamento é qualitativamente diferente. Enquanto alguns sabores (ex: $\xi_4\xi_5$) descem em massa, outros (ex: $\xi_i\xi_5$, $\xi_i\xi_4$, $\xi_4$) sobem ou permanecem relativamente constantes em baixa anisotropia antes de eventualmente seguirem o novo padrão.
• Modelagem Empírica: Os autores desenvolveram dois modelos para descrever os desdobramentos de sabor aHISQ:
  1. Inspirado em Teoria de Perturbação de Chiral (ChPT): Uma extensão da análise de Lee-Sharpe onde as constantes de energia baixa ($C_i, D_i$) tornam-se dependentes de $\xi$. Eles descobriram que, para aHISQ, os desdobramentos são dominados por $C_4(\xi)$ e $D_4(\xi)$, com outras constantes próximas de zero.
  2. Modelo de Perturbação de Quebra de Simetria: Um modelo empírico que trata a anisotropia como uma perturbação que redistribui a energia baseada no número de elos (links) espaciais ($l_\sigma$) e temporais ($l_\tau$) desbalanceados nos operadores interpoladores. Este modelo captura com sucesso o padrão equidistante observado dos desdobramentos e a dependência específica de $C_4$ e $D_4$ em relação a $\xi$.

Significância
O artigo afirma que o formalismo aHISQ desenvolvido e os procedimentos de ajuste fornecem uma base robusta para gerar ensembles anisotrópicas totalmente dinâmicas. A principal significância reside na eficiência demonstrada: a capacidade de ajustar a ação aHISQ com uma busca unidimensional simplificada para a anisotropia de férmion torna a geração de ensembles de alta estatística e grande $N_\tau$ viável. Essa capacidade é crucial para abordar o problema da reconstrução espectral no regime de plasma de quarks-glúons, onde as abordagens atuais de férmions de Wilson são limitadas pelo custo computacional e pela extensão temporal. O artigo valida o código reproduzindo resultados da literatura para gauge puro, espectros naïve quenched e espectros dinâmicos naïve quenched, e apresenta o primeiro ajuste detalhado e análise espectral para a ação aHISQ em lattices anisotrópicas.