Resumo Técnico: Um Algoritmo Quântico sem Penalidade para Encontrar Autoestados de Energia

Declaração do Problema

Encontrar os autovalores e autoestados de Hamiltonianos de muitos corpos é um desafio fundamental na física e na ciência computacional. Métodos clássicos enfrentam escalamento exponencial ($O(2^{2n})$ de armazenamento e $O(2^{3n})$ de tempo) à medida que o número de qubits $n$ aumenta. Embora algoritmos quânticos ofereçam um caminho a seguir, as abordagens existentes possuem limitações significativas:

• Eigensolver Quântico Variacional (VQE): Uma abordagem híbrida quântico-clássica que frequentemente sofre de "platôs áridos" (gradientes que desaparecem) durante a otimização. Encontrar estados excitados tipicamente requer termos de penalidade para suprimir contribuições do estado fundamental, aumentando a complexidade e o potencial de erro.
• Recozimento Quântico: Requer tempos de evolução longos para garantir adiabaticidade e enfrenta dificuldades com pequenos gaps de energia entre o estado fundamental e os estados excitados.
• Evolução no Tempo Imaginário (ITE) & QITE: Embora promissoras, a ITE padrão enfrenta problemas de escalabilidade na realização de operadores não unitários. A QITE frequentemente depende de otimização variacional ou termos de penalidade para estados excitados, reintroduzindo os desafios dos platôs áridos e exigindo ansatzes complexos.
• Outros Métodos: Técnicas como Diagonalização por Filtro Quântico (QFD) e Estimativa de Fase Quântica Filtrada (FQPE) dependem fortemente de estados de referência específicos ou funções de filtro otimizadas.

Os autores identificam a necessidade de um algoritmo totalmente quântico que possa encontrar tanto estados fundamentais quanto excitados sem funções de penalidade, etapas variacionais ou loops híbridos clássico-quânticos, evitando assim os platôs áridos.

Metodologia

O algoritmo proposto é um método totalmente quântico e sem penalidade baseado em amostragem estocástica de uma matriz de densidade de estado misto combinada com simulação baseada em estado.

1. Decomposição do Hamiltoniano e Mapeamento da Matriz de Densidade

O algoritmo assume que o Hamiltoniano alvo $H$ pode ser decomposto em dois componentes solucionáveis, $H = h_1 + h_2$, onde os autovalores ($\lambda^j_i$) e autoestados ($|\psi^j_i\rangle$) de $h_1$ e $h_2$ são conhecidos e fáceis de preparar.

• Os autores deslocam e reescalam os autovalores de $h_1$ e $h_2$ para criar uma distribuição de probabilidade válida.
• Eles constroem uma matriz de densidade de estado misto $\rho$ equivalente ao Hamiltoniano deslocado:
  $$\rho = \frac{1}{C} (H + (c_1 + c_2)I) = \sum_{i,j} p^j_i |\psi^j_i\rangle\langle\psi^j_i|$$
  onde $p^j_i$ são probabilidades derivadas dos autovalores deslocados. Os autoestados de $\rho$ são idênticos aos de $H$.

2. Evolução Estocástica no Tempo Imaginário (ITE) para o Estado Fundamental

Para encontrar o estado fundamental, o algoritmo realiza ITE em $\rho$ usando amostragem estocástica:

• Amostragem: Em cada passo de tempo $\eta$, um projetor $P^j_i = |\psi^j_i\rangle\langle\psi^j_i|$ é amostrado da distribuição $\{p^j_i\}$.
• Operador de Evolução: Em vez de implementar o operador não unitário completo $e^{-\eta \rho}$, o algoritmo implementa o operador $I - \eta P^j_i$ para o projetor amostrado.
• Implementação do Circuito: Isso é alcançado por meio de um circuito unitário com pós-seleção envolvendo um qubit auxiliar. A unitária $U$ é uma porta SU(2) controlada por $n$ qubits. Se o sistema estiver no estado correspondente ao projetor amostrado, o ancila sofre uma rotação; caso contrário, permanece inalterado. A pós-seleção do ancila no estado $|0\rangle$ produz a operação desejada $I - \eta P^j_i$.
• Convergência: Repetir esse processo leva o sistema ao estado fundamental de $\rho$ (e, portanto, de $H$) sem otimização variacional.

3. Simulação Baseada em Estado para Estados Excitados

Para encontrar estados excitados, o algoritmo deve suprimir as contribuições de estados de menor energia já encontrados (por exemplo, o estado fundamental $P_g$).

• Desafio: Implementar diretamente projetores para estados fundamentais altamente emaranhados é difícil.
• Solução: Os autores utilizam um protocolo de simulação baseada em estado. Isso envolve um sistema dual de $n$ qubits (ancila) preparado no estado fundamental conhecido $P_g$, um qubit de controle e o sistema principal.
• Mecanismo: Uma porta SWAP controlada é aplicada entre o sistema principal e o ancila, condicionada ao qubit de controle. Medir o qubit de controle na base $|+\rangle$ implementa efetivamente a evolução $e^{-\eta P_g}$ (aproximada como $I - \eta P_g$), projetando fora a componente do estado fundamental.
• Procedimento de Elevação: Para encontrar o primeiro estado excitado, a probabilidade de amostrar o projetor do estado fundamental é artificialmente "elevada" (aumentada) na distribuição da matriz de densidade. O algoritmo então realiza ITE nessa distribuição modificada. Quando o projetor do estado fundamental é amostrado, a simulação baseada em estado é usada para filtrá-lo. Esse processo pode ser iterado para encontrar estados excitados mais altos.

4. Variante para Hamiltonianos Gerais

Para Hamiltonianos que não são facilmente divididos em duas partes solucionáveis, os autores propõem uma variante onde $H$ é tratado como uma combinação linear de operadores de Pauli ($H = \sum \alpha_i O_i$). A amostragem estocástica é aplicada diretamente a esses operadores de Pauli, com implementações de circuito semelhantes para as etapas de evolução.

Resultados Chave

Os autores validaram o algoritmo usando simulações clássicas do Modelo de Ising com Campo Transverso (uma cadeia 1D com condições de contorno periódicas), que serve como benchmark porque pode ser dividido em componentes solucionáveis nas bases $Z$ e $X$.

• Estado Fundamental (sistema de 4 qubits): Alcançou uma fidelidade de 99,9% com o estado fundamental exato.
• Estados Excitados (sistema de 4 qubits):
  • 1º estado excitado: 98,6% de fidelidade.
  • 2º estado excitado: 96,8% de fidelidade.
  • 3º estado excitado: 96,6% de fidelidade.
• Teste de Escalabilidade (sistema de 10 qubits): O algoritmo encontrou com sucesso o estado fundamental (98,6% de fidelidade) e o primeiro estado excitado (97,8% de fidelidade). Também resolveu o terceiro autoestado mais baixo (90,2% de fidelidade) apesar da quase degenerescência entre os dois primeiros estados.
• Variante de Pauli: A variante de amostragem de operadores de Pauli também demonstrou altas fidelidades (97,1%–99,5%) para os quatro primeiros autoestados no caso de 4 qubits.
• Escalamento de Recursos: A análise teórica indica que o algoritmo escala com $O(n)$ portas de dois qubits por passo, uma melhoria significativa em relação aos métodos ITE baseados em Trotter, que escalonam como $O(n^2)$.

Significado e Alegações

O artigo alega apresentar um algoritmo totalmente quântico que aborda o problema dos autoestados sem as desvantagens dos métodos híbridos ou variacionais atuais:

1. Sem Penalidade: Encontra estados excitados sem adicionar termos de penalidade a uma função de perda, evitando a necessidade de ansatzes complexos e o risco de platôs áridos.
2. Sem Etapas Híbridas: Elimina a necessidade de loops de otimização clássica, tornando-o um procedimento puramente quântico.
3. Ampla Aplicabilidade: A suposição de que $H$ pode ser dividido em partes solucionáveis é argumentada como não restritiva, cobrindo muitos modelos físicos e até problemas QMA-completos.
4. Robustez: O método demonstra a capacidade de encontrar múltiplos autoestados de baixa energia sequencialmente, mesmo em sistemas com níveis de energia quase degenerados.

Os autores reconhecem que a taxa de sucesso geral depende da pós-seleção, o que pode ser um desafio em dispositivos Quânticos de Escala Intermediária e Ruidosa (NISQ) atuais devido a erros de porta e decoerência. No entanto, eles postulam que o método oferece um caminho concreto para usar circuitos quânticos para resolver problemas de autoestados intratáveis classicamente, servindo como uma adição fundamental à caixa de ferramentas da computação quântica.