Low-Depth Initial-State Preparation for Ground-State Energy Estimation of Two-Dimensional Strongly Correlated Systems

Resumo Técnico: Preparação de Estado Inicial de Baixa Profundidade para Estimativa de Energia do Estado Fundamental de Sistemas Fortemente Correlacionados Bidimensionais

Definição do Problema
Algoritmos quânticos para estimativa da energia do estado fundamental, como a Estimativa de Fase Quântica (QPE), requerem um estado inicial com fidelidade não desprezível ao verdadeiro estado fundamental. A probabilidade de extrair com sucesso a energia do estado fundamental escala com essa fidelidade ($F_0$); se $F_0$ decair exponencialmente com o tamanho do sistema, o algoritmo torna-se intratável devido ao número de repetições exigido ($1/F_0$). Embora os Solucionadores de Autovalores Variacionais (VQE) ofereçam uma abordagem híbrida, eles enfrentam desafios em relação aos custos de avaliação de gradiente e à simulabilidade clássica de circuitos rasos. Além disso, compilar aproximações clássicas precisas (ex: Estados de Produto de Matrizes) em circuitos quânticos pode resultar em uma profundidade de circuito excessiva, o que é proibitivo para computadores quânticos tolerantes a falhas precoces (FTQC), onde recursos como qubits lógicos e estados mágicos são limitados. O desafio central é projetar circuitos quânticos rasos e eficientes em termos de recursos que preparem estados iniciais de alta fidelidade para sistemas fortemente correlacionados, especificamente o modelo de Hubbard bidimensional (2D), sem depender de loops de otimização profundos ou intratabilidade clássica.

Metodologia
Os autores propõem uma estratégia para construir circuitos quânticos rasos aproveitando a conexão de acoplamento forte entre o modelo de Heisenberg 2D e o modelo de Hubbard em meio preenchimento via transformação de Schrieffer–Wolff (SW). A metodologia procede em três estáções principais:

1. Construção do Circuito de Heisenberg:

   • Ansatz: O circuito visa o subespaço $S=0$ (singlete) do modelo de Heisenberg antiferromagnético. O estado inicial é um produto de dímeros de singlete ($| \Phi \rangle$) definidos em um emparelhamento perfeito da rede.
   • Dinâmica: O estado é evoluído usando camadas de portas SWAP exponencializadas ($e$SWAP). Essas portas respeitam a simetria SU(2), garantindo que o estado permaneça no subespaço de singlete.
   • Estrutura: As ligações da rede são particionadas em quatro subconjuntos disjuntos ($B_{x0}, B_{x1}, B_{y0}, B_{y1}$). O circuito aplica portas $e$SWAP em uma sequência específica ($B_{y1} \to B_{y0} \to B_{x1} \to B_{x0}$) através de múltiplos blocos ($N_b$). Esta estrutura mimetiza a física de ligação de valência ressonante (RVB) do estado fundamental de Heisenberg.
   • Mapeamento Fermiônico: Os estados de spin são mapeados para o subespaço fermiônico de ocupação única usando um embedding unitário local ($V$) que introduz qubits ancilares.

2. Extensão de Hubbard via Portas SW:

   • Para capturar flutuações de carga (dublons e holons) presentes no modelo de Hubbard em $U/t$ finito, o circuito de Heisenberg é "vestido" com portas de flutuação de carga.
   • Essas portas são derivadas diretamente do gerador de primeira ordem de SW ($\Omega$), aproximado via decomposição de Suzuki–Trotter em unitárias de dois sítios ($W_{\langle i,j \rangle}$).
   • O ansatz variacional final combina a evolução de Heisenberg ($U_B$) e o bloco de flutuação de carga ($W_B$), parametrizados pelos ângulos de rotação $\theta$ e $\phi$.

3. Estimativa de Fidelidade e Benchmarking:

   • Transferência de Parâmetros: Parâmetros otimizados em uma rede pequena de $4 \times 4$ são transferidos para redes maiores ($6 \times 6$ a $10 \times 10$) sem nova otimização.
   • Simulação: Os estados do circuito são simulados usando Estados de Pares Entrançados Projetados (PEPS) com atualizações simples guiadas por propagação de crença (BP).
   • Limites de Fidelidade: Como o estado fundamental exato é desconhecido para sistemas grandes, os autores estimam um limite inferior na fidelidade do estado fundamental ($F_0$). Isso é alcançado por:
     • Calculando a fidelidade ($F_{ref}$) entre o estado do circuito e um estado de referência de Monte Carlo Variacional (VMC).
     • Limitando a fidelidade do estado VMC em relação ao estado fundamental exato ($F_{ref,0}$), usando a energia de referência, a energia do estado fundamental exata (de Expansão de Séries Estocásticas) e o gap de excitação de singlete ($\Delta_0$).
     • Aplicando a desigualdade triangular para ângulos de Fubini–Study para derivar um limite inferior rigoroso para $F_0$.

Principais Contribuições e Resultados

• Transferência de Parâmetros Eficaz: Para o modelo de Heisenberg, parâmetros otimizados em uma rede $4 \times 4$ produzem limites inferiores de fidelidade robustos para redes de até $10 \times 10$ sem reotimização. Especificamente, para $N_b=3$, o limite inferior da fidelidade do estado fundamental estimado permanece positivo ($\approx 0,028$) mesmo para a rede $10 \times 10$, sugerindo que o estado inicial é adequado para amplificação de amplitude.
• Melhoria da Fidelidade de Hubbard: Cálculos exatos em uma rede $4 \times 4$ demonstram que incorporar as portas de flutuação de carga derivadas de SW melhora significativamente a fidelidade para o estado fundamental de Hubbard. Para $N_b=4$, a fidelidade aumenta de $\approx 0,44$ (apenas Heisenberg) para $\approx 0,89$ (com ângulos de carga otimizados).
• Eficiência de Recursos: Os autores estimam o custo de preparação de estado para uma rede $10 \times 10$. Assumindo síntese Clifford+$T$ com um erro por rotação de $10^{-5}$, o custo total é da ordem de $10^5$ portas $T$. Isso coloca o custo de preparação bem dentro do regime de "megaquops" para FTQCs iniciais e significativamente abaixo dos custos associados aos algoritmos subsequentes de estimativa de energia do estado fundamental.
• Escalabilidade: A profundidade do circuito escala linearmente com o número de ligações da rede ($|B|$), e a estratégia de compartilhamento de parâmetros permite a construção de circuitos para sistemas grandes baseados na otimização de sistemas pequenos.

Significância e Alegações
O artigo afirma oferecer uma rota viável para a preparação de estado inicial de baixo custo para a estimativa da energia do estado fundamental em sistemas fortemente correlacionados. Ao utilizar a conexão física entre os modelos de Heisenberg e Hubbard, os autores demonstram que circuitos rasos podem alcançar fidelidade de estado fundamental não desprezível em redes além do alcance da diagonalização clássica exata. Os resultados sugerem que esses circuitos são candidatos atraentes para o regime de FTQC inicial, onde minimizar a profundidade do circuito e a contagem de portas $T$ é crítico. O trabalho apoia a hipótese de que simulações de redes pequenas podem efetivamente guiar o design de circuitos quânticos para sistemas maiores, classicamente intratáveis, desde que os parâmetros sejam transferíveis e a profundidade do circuito seja mantida rasa. Os autores observam explicitamente que, embora seus limites sejam rigorosos, as fidelidades reais podem ser maiores, e trabalhos futuros poderiam envolver ansatzes variacionais de Hamiltoniana mais profundos ou técnicas de simulação mais precisas (ex: PEPS de atualização total) para refinar ainda mais essas estimativas.