(2+1)D quantum electrodynamics at finite density on a quantum computer

Problema

A simulação de teorias de gauge em rede, particularmente a Cromodinâmica Quântica (QCD), em densidade finita permanece um desafio significativo para os métodos computacionais clássicos devido ao problema do sinal e à complexidade da dinâmica em tempo real. Embora modelos de dimensões inferiores como o modelo de Schwinger (1+1)D tenham sido simulados com sucesso em hardware quântico, estender esses estudos para a Eletrodinâmica Quântica (QED3) em (2+1)D com múltiplas sabores de férmions em densidade finita ainda não foi demonstrado em dispositivos quânticos. A QED3 é fisicamente relevante como uma descrição efetiva para líquidos de spins quânticos e supercondutores de cuprato, embora sua estrutura de fase em densidade finita não seja totalmente compreendida analiticamente.

Metodologia

Os autores utilizam o Algoritmo Variacional de Eigensolver Quântico (VQE) para simular um sistema QED3 em (2+1)D com duas sabores de férmions em uma pequena rede.

• Formulação do Modelo: O estudo utiliza um Hamiltoniano de rede com férmions escalonados e links de gauge compactos. O Hamiltoniano total inclui termos de campo elétrico, plaqueta magnética, massa, salto cinético e potencial químico. Os operadores de campo elétrico são truncados para um intervalo finito ($l=1$).
• Invariância de Gauge: Para impor a lei de Gauss, os autores utilizam um ansatz invariante de gauge. Para a rede específica de $2 \times 2$ com condições de contorno abertas (OBC), a lei de Gauss permite que três dos quatro links de gauge sejam expressos em termos do link de gauge dinâmico restante. Isso reduz significativamente o espaço de Hilbert do campo de gauge, exigindo a simulação de apenas um link de gauge dinâmico, enquanto retém a conectividade 2D completa do setor fermiônico.
• Ansatz do Circuito Quântico: O ansatz variacional consiste em três componentes:
  1. Codificação de Férmions: Utiliza o mapeamento de Jordan-Wigner para converter os graus de liberdade fermiônicos (duas sabores por sítio) em qubits.
  2. Codificação de Gauge: Emprega uma codificação de código Gray para os estados de campo elétrico truncados em dois qubits.
  3. Emaranhamento: Inclui portas controladas $U^{xy}$ para acoplar os setores de matéria e gauge enquanto preserva a invariância de gauge, juntamente com camadas de mistura entre sabores e portas de correção de fase ($R_Z$).
• Otimização: O circuito contém 24 parâmetros variacionais. Simulações clássicas são realizadas usando o otimizador COBYLA para encontrar os parâmetros ótimos, que são então implantados em hardware quântico. Os autores observam que otimizadores mais avançados como NFT foram considerados inadequados devido à estrutura específica de seus geradores de portas.
• Execução de Hardware: As execuções de inferência foram conduzidas no dispositivo IBM ibm_marrakesh. Os autores aplicaram Extrapolação de Zero-Ruído (ZNE) usando dobramento de porta (gate folding) para mitigar o ruído do hardware.

Principais Contribuições

1. Primeira Prova de Conceito: Este trabalho apresenta a primeira demonstração de QED3 com múltiplas sabores de férmions em densidade finita em hardware quântico.
2. Desenvolvimento de Protocolo: Os autores generalizam uma construção anterior de QED3 para incluir duas sabores de férmions e densidade finita, codificando explicitamente a lei de Gauss na estrutura do circuito para reduzir os requisitos de qubits.
3. Análise da Estrutura de Fase: O artigo deriva uma fórmula para os potenciais químicos críticos onde ocorrem cruzamentos de níveis induzidos por densidade em (2+1)D, analogamente a estudos anteriores em (1+1)D.
4. Análise de Ruído: Os autores fornecem uma análise quantitativa da sensibilidade do operador, demonstrando que o operador de número é significativamente mais robusto contra o ruído do hardware do que o Hamiltoniano devido às diferenças em suas normas de operador (números de condição).

Resultados

• Simulações Clássicas: Simulações clássicas sem ruído em uma rede $2 \times 2$ com acoplamento de gauge $g=5$ identificaram com sucesso cruzamentos de níveis induzidos por densidade. A diferença no número de partículas entre as duas sabores ($\langle \hat{N}_1 - \hat{N}_2 \rangle$) exibiu comportamento em degraus, indicando transições entre diferentes setores de número de partículas.
• Região Intermediária: Uma pequena fase intermediária foi observada perto da região de transição. Os autores descobriram que essa região torna-se mais estreita e difícil de resolver conforme a força de acoplamento $g$ aumenta, enquanto a entropia de emaranhamento entre o link de gauge e o sistema fermiônico é máxima nesta região intermediária.
• Inferência de Hardware: Quando implantada no processador quântico IBM, a simulação capturou com sucesso a tendência geral dos cruzamentos de níveis. O ponto crítico foi estimado em $\mu \approx 13 \pm 0.7$.
  • Os resultados do operador de número ($\hat{N}_1, \hat{N}_2$) permaneceram próximos dos valores esperados, apesar do ruído.
  • Os valores esperados de energia mostraram desvios maiores dos resultados de diagonalização exata (ED), consistentes com o maior número de condição do Hamiltoniano.
  • Devido às limitações de profundidade de circuito, o hardware não conseguiu resolver a estreita região intermediária, resultando em um ponto de transição médio em vez de dois pontos críticos distintos.

Significância e Alegações

O artigo afirma ser um "primeiro passo" em direção à evolução em tempo real de sistemas de densidade finita em (2+1)D. Os autores enfatizam que, embora as limitações atuais do hardware (ruído e profundidade de circuito) impeçam a resolução de estruturas de fase de grão fino em redes maiores, o estudo valida a viabilidade do protocolo VQE invariante de gauge proposto.

Os autores afirmam explicitamente que não existem resultados analíticos ou numéricos para os pontos de cruzamento de nível em QED (2+1)D para o regime de parâmetros estudado, tornando a identificação dessas transições em hardware quântico uma contribuição novel. Eles concluem que, embora o sistema $2 \times 2$ atual seja pequeno, o método é escalável em princípio usando ansatz adaptativos avançados (como SC-ADAPT-VQE) para sistemas maiores, embora isso permaneça um assunto para investigações futuras. O trabalho serve como um benchmark para estudos futuros que visam explorar o diagrama de fase de QED3 e sistemas relacionados de matéria condensada em computadores quânticos.