The leading Lyapunov exponent in the glasma

Resumo Técnico: O Expoente de Lyapunov Líder no Glasma

Enunciado do Problema
O estágio inicial de colisões de íons pesados é dominado por um estado de glúons altamente sobreocupado e longe do equilíbrio, conhecido como glasma. Embora a evolução deste sistema em direção a um plasma de quarks e glúons termalizado seja um tópico central na física nuclear de alta energia, os mecanismos que impulsionam a termalização e a produção de entropia permanecem mal compreendidos. Especificamente, a relação entre a chaoticidade clássica dos campos de gauge iniciais e a escala de tempo de termalização requer esclarecimento. A chaoticidade clássica é caracterizada pelo crescimento exponencial de pequenas perturbações, quantificada por expoentes de Lyapunov. Embora tal comportamento tenha sido estudado em sistemas de Yang-Mills clássicos (CYM) em equilíbrio térmico, sua aplicabilidade ao glasma não-equilibrado e invariante de boost é distinta devido à sobreocupação dos modos de momento dominantes, o que torna o sistema tratável via teoria de campo clássica.

Metodologia
Os autores investigam a dinâmica caótica do glasma invariante de boost em 2+1 dimensões usando o modelo de McLerran-Venugopalan (MV) dentro do framework de campos de gauge clássicos. O estudo emprega uma abordagem numérica "pedestre" e robusta:

1. Condições Iniciais: O sistema é inicializado usando o modelo MV, onde as cargas de cor são distribuídas como um ensemble Gaussiano em uma rede transversal. Os campos de gauge iniciais são derivados resolvendo as equações de Yang-Mills no gauge de luz e transformando para o gauge temporal ($A_\tau = 0$) dentro do cone de luz futuro.
2. Estratégia de Perturbação: Para sondar a chaoticidade, os autores inicializam duas configurações de campo quase idênticas em $\tau=0$. Uma configuração é perturbada por flutuações infinitesimalmente pequenas no campo elétrico longitudinal ($E_\eta$) ou no campo magnético longitudinal ($B_\eta$).
3. Filtragem de Momento: Para testar a dependência de escala do crescimento caótico, as perturbações são inicializadas usando três filtros distintos de espaço de momento:
   • Ruído Branco: Ruído não correlacionado através de todos os modos de momento.
   • Filtro Gaussiano: Suprime modos de alto momento (UV).
   • Filtro de Lei de Potência: Permite uma cauda UV mais pesada.
   • Filtro de Casca (Band-Pass): Isola cascas de momento específicas.
4. Implementação Numérica: A evolução é simulada usando um solver de teoria de rede de gauge em tempo real com um algoritmo leapfrog. O estudo utiliza o grupo de cor $SU(2)$ em redes de tamanho $N_\perp = 256$ (com verificações para $N_\perp = 128$) para garantir estatísticas suficientes ($N_{events} = 5000$) e para verificar a independência de parâmetros de discretização de rede ($g^2\mu a_\perp$ e $g^2\mu L_\perp$).
5. Extração do Expoente de Lyapunov: O crescimento da norma da perturbação, $\langle \text{Tr}(\delta F_\eta)^2 \rangle$, é monitorado. Os autores ajustam o comportamento de tempo tardio à forma funcional $\exp(\lambda \sqrt{g^2\mu \tau})$, onde $\lambda$ é o expoente de Lyapunov adimensional. Uma análise de regressão de múltiplos intervalos sistemáticos é empregada para identificar a janela de ajuste ideal que evita transientes iniciais e saturação não-linear.

Resultados Principais

• Crescimento Exponencial com a Raiz Quadrada do Tempo: As perturbações exibem crescimento exponencial dependente da raiz quadrada do tempo próprio, $\tau$. Isso é uma consequência direta da expansão invariante de boost do sistema. A dependência temporal é bem descrita por $\delta F(\tau) \sim \exp(\lambda \sqrt{g^2\mu \tau})$.
• Expoente de Lyapunov Líder: O expoente de Lyapunov líder extraído é encontrado como $\lambda \approx 0,39 \pm 0,02$. A incerteza é dominada por variações sistemáticas nas condições iniciais e parâmetros de rede, em vez de erros estatísticos.
• Independência de Escala: A taxa de crescimento $\lambda$ é notavelmente insensível à estrutura de momento da perturbação inicial. Quer seja inicializada com ruído branco, filtro Gaussiano, de lei de potência ou filtros de casca de momento específicos (abrangendo escalas infravermelhas a ultravioletas), o expoente extraído permanece consistente dentro da incerteza declarada. Isso sugere que a instabilidade caótica é uma propriedade coletiva e invariante de escala da dinâmica não-linear de Yang-Mills, em vez de ser impulsionada por modos de momento específicos.
• Acoplamento de Estados de Polarização: O modo instável acopla-se tanto aos componentes elétricos quanto magnéticos longitudais. Perturbar um componente (ex: $E_\eta$) leva ao crescimento exponencial de ambos $E_\eta$ e $B_\eta$ com o mesmo expoente de Lyapunov, indicando que a dinâmica caótica mistura os estados de polarização que são independentes em uma aproximação linearizada.
• Robustez à Discretização: O valor de $\lambda$ mostra-se independente do espaçamento da rede (cutoff UV) e do volume físico (cutoff IR) dentro dos intervalos testados, confirmando que o resultado é uma propriedade física genuína do glasma e não um artefato numérico.

Significância e Alegações
O artigo afirma demonstrar a existência de modos caóticos de crescimento exponencial no glasma invariante de boost e fornecer um valor numérico para o expoente de Lyapunov líder. Os autores interpretam essa taxa de crescimento como uma escala de tempo para a produção de entropia e termalização nos estágios iniciais de colisões de íons pesados.

A principal significância do trabalho reside em estabelecer que a taxa de crescimento caótico é uma propriedade universal do glasma sobreocupado, independente dos detalhes específicos de como as perturbações são inicializadas (amplitude, escala de momento ou componente de campo). Os autores notam modestamente que, embora a taxa de crescimento seja robusta, a amplitude do modo instável depende da estrutura de momento inicial, com maiores sobreposições para graus de liberdade infravermelhos.

O artigo conclui que a taxa de crescimento $\gamma(\tau) = \frac{d}{d\tau} \ln \delta F \approx (0,1 \pm 0,005) \sqrt{g^2\mu/\tau}$ é parametricamente semelhante às escalas de massa de plásmon ou Debye no glasma em expansão. Os autores identificam direções de pesquisas futuras, incluindo a extensão do estudo para $SU(3)$, a análise do espectro completo de expoentes de Lyapunov e a investigação da estrutura de momento do modo instável via análise de Fourier com fixação de gauge, mas não alegam ter resolvido o mecanismo completo de termalização ou a conexão com a instabilidade de Weibel neste trabalho.