Quantum speed-up for solving the one-dimensional Hubbard model using quantum annealing

Resumo Técnico: Aceleração Quântica para Resolver o Modelo de Hubbard Unidimensional Usando Recozimento Quântico

Enunciação do Problema
O modelo de Hubbard unidimensional é uma estrutura fundamental na física da matéria condensada para compreender sistemas de elétrons correlacionados. Embora o modelo seja "integrável"—ou seja, existam resultados analíticos exatos para seu estado fundamental no limite termodinâmico por meio das equações de Bethe-ansatz—resolver essas equações para sistemas finitos com fronteiras abertas apresenta desafios computacionais significativos. Especificamente, embora a energia do estado fundamental possa ser encontrada em tempo polinomial resolvendo as equações de Bethe-ansatz não lineares acopladas, determinar a completa função de onda do estado fundamental (os coeficientes) exige calcular um número exponencialmente grande de termos, demandando recursos computacionais exponenciais. Os autores investigam se algoritmos quânticos, especificamente o recozimento quântico, podem oferecer uma aceleração sobre esses métodos clássicos para encontrar o estado fundamental de sistemas de Hubbard unidimensionais finitos e semi-preenchidos.

Metodologia
Os autores implementam uma simulação baseada em portas do recozimento quântico em um simulador de computador quântico universal (JUQCS) para resolver o Hamiltoniano de Hubbard. A metodologia envolve três etapas principais:

1. Mapeamento do Hamiltoniano: O Hamiltoniano de Hubbard fermiónico é mapeado para um Hamiltoniano de qubits usando a transformação de Jordan-Wigner. Isso converte operadores de criação e aniquilação de férmions em matrizes de Pauli atuando sobre $2L$ qubits (onde $L$ é o número de sítios da rede), preservando as relações de anticomutação.
2. Preparação do Estado Inicial: O processo de recozimento começa com o estado fundamental de um Hamiltoniano de salto não interativo ($H_I$). Os autores utilizam um algoritmo baseado em rotações de Givens para preparar esse estado inicial, que corresponde ao preenchimento dos estados de momento de menor energia para um número dado de elétrons com spin para cima e spin para baixo.
3. Dinâmica de Recozimento: O sistema evolui sob um Hamiltoniano dependente do tempo $H(s) = (1-s)H_I + sH_H$, onde $s$ varia de 0 a 1. A evolução temporal é simulada usando a decomposição de Suzuki-Trotter (especificamente uma aproximação de segunda ordem) para dividir a evolução contínua em passos de tempo discretos. Os operadores unitários para cada passo são decompostos em portas universais de 1 qubit e 2 qubits (portas Hadamard, $R_Z$, $R_{ZZ}$ e $X$).
4. Análise de Escala: Simulações foram realizadas para tamanhos de sistema até $L=20$ (40 qubits). A métrica de desempenho é a energia residual $\Delta E = \langle \psi(1) | H_H | \psi(1) \rangle - E_0$, onde $E_0$ é a energia exata do estado fundamental calculada via Bethe-ansatz. Os autores analisam como o tempo total de recozimento necessário ($T_A$) escala com o tamanho do sistema ($L$) para atingir uma precisão específica.

Principais Contribuições e Resultados

• Implementação Baseada em Portas: O artigo fornece uma construção concreta para simular o recozimento quântico para o modelo de Hubbard em um computador quântico baseado em portas, detalhando a decomposição do circuito e os requisitos de recursos (contagens de portas).
• Escala do Tempo de Recozimento: Para sistemas semi-preenchidos ($2N_\downarrow = N = L$), os autores encontram que a energia residual escala como $\Delta E \propto L^{1.234} / T_A^2$ para cronogramas de recozimento lineares. Para atingir uma precisão fixa $\Delta E^*$, o tempo de recozimento necessário escala como $T_A \propto L^{0.62}$ para tamanhos de sistema abaixo de um certo limiar $L^*$.
• Regime de Acoplamento Forte: Simulações para fortes intensidades de interação ($U/t_H = 8, 16$) confirmam que o início do regime de escala permanece benigno, com o tempo de recozimento limitado por uma função linear do tamanho do sistema ($T_A \propto L^{0.995}$) para $L$ arbitrariamente grande.
• Comparação com Métodos Clássicos: Os autores comparam suas descobertas com algoritmos clássicos de busca de raízes para as equações de Bethe-ansatz. Embora encontrar a energia via Bethe-ansatz seja polinomial ($T_{BA} \propto L^{3.3}$), obter os coeficientes completos do estado fundamental é exponencialmente caro. Em contraste, a abordagem de recozimento quântico produz o estado fundamental com uma escala de tempo e recursos sub-linear a linear em relação ao tamanho do sistema.

Significado e Alegações
O artigo alega que o recozimento quântico fornece uma aceleração quântica substancial sobre algoritmos clássicos para encontrar o estado fundamental do modelo de Hubbard unidimensional.

• Aceleração Exponencial para Preparação de Estado: Os autores argumentam que, como os métodos clássicos exigem recursos exponenciais para determinar os coeficientes do estado fundamental (mesmo que a energia seja encontrada polinomialmente), a escala polinomial (especificamente sub-linear a linear) do tempo de recozimento quântico representa uma aceleração quântica exponencial no contexto da preparação do estado fundamental.
• Validação da Computação Adiabática: Os resultados suportam a utilidade da computação quântica adiabática para problemas quânticos de muitos corpos integráveis. A escala polinomial observada é atribuída à ausência de transições de fase no modelo de Hubbard unidimensional em intensidade de interação não nula, o que impede que o gap de energia feche exponencialmente.
• Perspectiva Futura: Embora o estudo foque no caso integrável 1D, os autores sugerem que esses resultados encorajam investigações adicionais em sistemas não integráveis (como o modelo de Hubbard 2D). Eles observam que, embora transições de fase em dimensões superiores possam alterar o comportamento do gap, a vantagem demonstrada no caso 1D reforça o potencial do recozimento quântico para resolver problemas complexos de muitos corpos.

Os autores permanecem modestos quanto à aplicação prática imediata, reconhecendo que os cálculos atuais com milhares de portas estão além das capacidades do hardware de curto prazo, mas que a escala de recursos augura bem para futuras arquiteturas tolerantes a falhas.