Lattice Quantum Chromodynamics for Quantum Simulations

Resumo Técnico: Cromodinâmica Quântica em Rede para Simulações Quânticas

Enunciado do Problema
A Cromodinâmica Quântica (QCD) é a teoria fundamental da interação forte, porém seu regime não perturbativo — essencial para compreender a estrutura de hádrons e a matéria de alta densidade — permanece computacionalmente desafiador. As abordagens padrão de QCD em rede dependem de integrais de caminho euclidianas e amostragem por importância, que falham na presença de um "problema de sinal", como em dinâmicas de tempo real ou quando um ângulo $\theta$ é incluído. Embora as formulações hamiltonianas ofereçam um caminho para simulações em tempo real, o crescimento exponencial do espaço de Hilbert com o volume da rede torna a simulação clássica intratável para sistemas grandes. Consequentemente, há uma necessidade crítica de frameworks que mapeiem teorias de gauge em rede para computadores quânticos, abordando especificamente a codificação dos graus de liberdade de gauge e fermiônicos, a construção de estados físicos (invariantes de gauge) e o desenvolvimento de algoritmos eficientes de evolução temporal.

Metodologia
Os autores desenvolvem um framework abrangente para simulações quânticas da teoria de gauge $SU(N_c)$ com quarks, focando em $N_c=3$ (QCD) em duas e três dimensões espaciais. A metodologia é construída sobre os seguintes pilares:

1. Formulação Hamiltoniana e Discretização:
   O trabalho utiliza a formulação hamiltoniana de Kogut-Susskind. Os autores consideram tanto fermiões staggered quanto de Wilson para abordar o problema do duplo de fermiões (fermion doubling). A teoria é discretizada em uma rede espacial cúbica com tempo contínuo. O Hamiltoniano é decomposto em termos elétricos ($H_E$), magnéticos ($H_B$), $\theta$ ($H_\theta$), cinéticos ($H_{kin}$) e de massa ($H_{mass}$). O termo $\theta$ é incluído como um acoplamento à densidade de carga topológica, introduzindo coeficientes complexos no Hamiltoniano.

2. Base de Representação e Espaço de Hilbert Físico:
   Uma contribuição central é a formulação do espaço de Hilbert físico em uma base de representação. Em vez de mapear elementos de grupo diretamente para qubits, os graus de liberdade de gauge nos elos (links) e os graus de liberdade fermiônicos nos sítios são codificados usando representações irredutíveis (irreps) de $SU(N_c)$.

   • Espaço de Hilbert do Elo (Link): Os elos são descritos por estados $|(R, r_R, r_L)\rangle$, onde $R$ é uma irrep e $r_{R,L}$ são vetores de base para as metades esquerda e direita do elo. Esta base diagonaliza os operadores de Casimir quadráticos.
   • Espaço de Hilbert do Sítio: Os fermiões são tratados via uma base de número de ocupação, que é então transformada em uma base de representação usando coeficientes de Clebsch-Gordan (CGCs). Isso permite que estados de partícula e antipartícula sejam classificados por suas propriedades de transformação sob $SU(N_c)$.
   • Invariância de Gauge: Estados físicos são construídos como singletos (representações triviais) formados pelo produto tensorial de todas as irreps que se encontram em um sítio da rede (elos e fermiões). Os autores definem explicitamente a convenção de "ordem do sítio" necessária para construir de forma inequívoca esses estados singletos, garantindo que a restrição da lei de Gauss seja satisfeita.

3. Truncamento e Codificação:
   Dado que o espaço de Hilbert de $SU(N_c)$ é de dimensão infinita, os autores implementam esquemas de truncamento. Estes incluem cortes no peso máximo das irreps ($n_1$) ou no autovalor do Casimir quadrático ($C_2$) para elos de gauge, e limites nos números de ocupação de cor para os fermiões. Os estados físicos são codificados em qubits usando um esquema de codificação local, onde os qubits representam irreps específicas e índices de multiplicidade de singlete. Este método é observado como sendo mais eficiente em termos de qubits do que codificações diretas para grandes espaços de Hilbert.

4. Algoritmos de Simulação Quântica:
   Os autores propõem um algoritmo de Trotterização para a evolução temporal. O operador de evolução temporal $e^{-iHt}$ é aproximado como um produto de unitários correspondentes aos termos decompostos do Hamiltoniano.

   • Construção de Portas: A evolução de operadores de plaqueta e de elo é mapeada para rotações de Givens (unitários de dois níveis) e deslocamentos de fase. Os autores detalham a construção de circuitos quânticos para estas operações, incluindo o tratamento de cordas de Jordan-Wigner (inversões de fase) necessários para as relações de anticomutação fermiônica.
   • Elementos de Matriz: O artigo fornece fórmulas explícitas para os elementos de matriz de todos os operadores do Hamiltoniano, expressos como produtos de "fatores de sítio" envolvendo CGCs. Estas fórmulas são essenciais para calcular os ângulos de rotação e as fases necessárias para as portas quânticas.

Principais Contribuições

• Framework para $SU(3)$ com Ângulo $\theta$: O artigo apresenta a primeira formulação e simulação de QCD em rede com um ângulo $\theta$ em três dimensões espaciais usando uma abordagem hamiltoniana em um computador quântico.
• Formalismo de Base de Representação: Fornece uma construção rigorosa do espaço de Hilbert físico usando irreps e singletes CGC, oferecendo uma alternativa escalável às bases de elementos de grupo.
• Detalhes Algorítmicos: O trabalho detalha os circuitos quânticos específicos para a evolução Trotterizada, incluindo a decomposição de termos complexos do Hamiltoniano (como o termo $\theta$ e termos cinéticos de Wilson) em portas elementares.
• Cálculos de Elementos de Matriz: Os autores derivam e tabulam os elementos de matriz para operadores de plaqueta, $\theta$ e cinéticos fermiônicos, fornecendo os dados necessários para implementar os circuitos de simulação.

Resultados
Os autores realizaram simulações sem ruído em redes pequenas (até uma rede periódica de $2\times2$ e um cubo de fronteira aberta de $2\times2\times2$) usando até 32 qubits em simuladores clássicos (Qiskit Aer e CUDA-Q). As simulações demonstraram:

• Efeitos de Truncamento do Espaço de Hilbert: Análise de como diferentes truncamentos de gauge e fermiões afetam observáveis como o campo elétrico médio ($\langle \Pi^2 \rangle$) e o número de ocupação de cor ($\langle N_q \rangle$). Os resultados indicaram que o acoplamento forte suprime estados de alta população, tornando truncamentos modestos eficazes para fermiões staggered, enquanto os fermiões de Wilson mostraram maior sensibilidade ao truncamento.
• Dependência de Parâmetros: A evolução de observáveis foi estudada em função da massa fermiônica e do parâmetro de Wilson $r$. Constatou-se que massas mais pesadas suprimem oscilações para estados ricos em fermiões.
• Dinâmica do Ângulo $\theta$: Simulações variando o ângulo $\theta$ revelaram comportamento aperiódico no regime de acoplamento forte, com dinâmicas distintas para $\theta$ positivo e negativo.
• Estados Hadrônicos e Dinâmica de Cordas: Os autores demonstraram a preparação de estados de "píon carregado" e a evolução de estados de corda (tanto cordas de loop quanto de fermiões). Eles observaram a dinâmica de quebra de corda, onde tubos de fluxo decaem em pares férmion-antiférmion, e rastrearam a probabilidade de quebra de corda versus configurações não quebradas.
• Potencial Químico de Bárion: Usando diagonalização exata em um único cubo, os autores computaram o custo de energia para adicionar um bárion, servindo como uma prova de conceito para o estudo de matéria de alta densidade.

Significância e Alegações
O artigo posiciona-se como uma contribuição fundamental para o desenvolvimento de simulações quânticas práticas para física fundamental. Os autores alegam que seu trabalho fornece os "ingredientes" necessários para projetar simulações quânticas de QCD em rede na formulação hamiltoniana.

• Prova de Conceito: A principal significância reside em demonstrar que a complexa maquinaria da teoria de gauge em rede $SU(3)$, incluindo o ângulo $\theta$ e a matéria fermiônica, pode ser codificada e evoluída em hardware quântico (simulado aqui).
• Caminho para Escalabilidade: A codificação local e a base de representação são apresentadas como um caminho viável para simulações maiores, potencialmente exigindo menos qubits do que codificações globais.
• Perspectiva Futida: Os autores afirmam modestamente que os resultados atuais são exploratórios e ainda não fazem declarações físicas sérias sobre QCD em tempo real devido aos pequenos tamanhos de rede, regime de acoplamento forte e truncamentos agressivos utilizados. Eles identificam a necessidade de Hamiltonianos melhorados (contra-termos), melhores algoritmos de preparação de estado e hardware tolerante a falhas para alcançar o limite contínuo e simular regimes fisicamente relevantes. O trabalho é estruturado como um degrau para simular matéria densa (ex: estrelas de nêutrons) e a dinâmica de axion-QCD no futuro.