Feynman's clock and hierarchy-informed sampling for quantum error mitigation

Resumo Técnico: O Relógio de Feynman e Amostragem Informada pela Hierarquia para Mitigação de Erros Quânticos

Definição do Problema
Algoritmos quânticos de curto prazo enfrentam limitações significativas devido à natureza ruidosa (NISQ - Noisy Intermediate-Scale Quantum) do hardware atual, o que prejudica a execução física de circuitos quânticos. Embora a correção de erros exija dispositivos de grande escala ainda não disponíveis, esquemas eficientes de mitigação de erro quântico (QEM) são necessários para aplicações de curto prazo. Os métodos existentes são fundamentalmente limitados, mas permanecem como o principal meio de reduzir o ruído quântico. Este artigo aborda a necessidade de uma técnica geral de QEM capaz de mitigar o ruído proveniente da execução de circuitos quânticos arbitrários, estendendo-se além das aplicações anteriores limitadas a simulações específicas de cadeias de spin.

Metodologia
A técnica proposta, uma extensão do esquema de Mitigação por Amostragem Informada por BBGKY (BBGKY-ISM), opera através de um processo de três estágios:

1. Mapeamento de Circuito-para-Hamiltoniano via Relógio de Feynman:
   Os autores mapeiam um circuito arbitrário de $N_Q$ qubits $U$, composto por $N_G$ portas, na evolução temporal de um sistema quântico governado por um Hamiltoniano independente do tempo. Isso utiliza a construção de computador quântico de Feynman, introduzindo um registro de relógio $H_C$ com $N_C = \lceil \log_2(N_G + 1) \rceil$ qubits auxiliares. O sistema evolui sob o Hamiltoniano do relógio de Feynman $H$, onde os estados do relógio rastreiam o progresso da computação do circuito. A dinâmica completa é descrita por $|\psi(t)\rangle = e^{-itH}|\psi(0)\rangle$.

2. Reconstrução via Testes de Hadamard:
   Em vez de simular a evolução temporal via Trotterização (o que amplificaria o ruído), o método reconstrói os valores esperados temporais de interesse usando testes de Hadamard em execuções parciais do circuito original.

• As medições Z no registro de dados são reinterpretadas como valores esperados dependentes do tempo de operadores Pauli-Z únicos, projetados sobre o estado final do relógio $|N_G\rangle_C$.
• Essas quantidades são decompostas em um conjunto de "quantidades de interesse" envolvendo strings de Pauli.
• As partes real e imaginária dos elementos de matriz necessários são medidas usando testes de Hadamard, evitando a profundidade de circuito elevada associada aos fatias de Trotter.
• Os valores finais de medição Z são extraídos somando essas quantidades de evolução temporal reconstruídas, ponderadas pela amplitude de probabilidade do estado do relógio, efetivamente removendo a dependência temporal.

3. Mitigação por Amostragem Informada por BBGKY:
   As estimativas ruidosas obtidas através dos testes de Hadamard são mitigadas usando o esquema BBGKY-ISM.

• Distribuição Informada pela Física: O método amostra candidatos de mitigação de uma distribuição de probabilidade codificada com a dinâmica sem ruído. Esta distribuição é definida por uma ação $S(\vec{x})$ que é minimizada quando a configuração $\vec{x}$ aproxima-se da evolução ideal e sem ruído.
• Restrições Hierárquicas: As restrições são derivadas de uma hierarquia de equações de movimento do tipo BBGKY para correladores de strings de Pauli. A hierarquia conecta as quantidades de interesse a correladores de ordem superior.
• Recozimento Simulado (Simulated Annealing): Um procedimento de Cadeia de Markov Monte Carlo (Recozimento Simulado) é usado para amostrar configurações que minimizam a ação. Um parâmetro $r$ (raio) determina quantos níveis da hierarquia são incluídos nas restrições. Um $r$ maior codifica uma porção maior da dinâmica ideal, tornando as restrições sobre a amostragem mais rigorosas.

Principais Contribuições

• Generalização do BBGKY-ISM: O artigo estende o esquema BBGKY-ISM de simulações específicas de cadeias de spin para circuitos quânticos arbitrários ao utilizar a construção do relógio de Feynman.
• Escalonamento de Sobrecarga Eficiente: Os autores demonstram que tanto as sobrecargas clássicas quanto as quânticas do método proposto escalam polinomialmente com o tamanho do circuito ($N_G$), o número de qubits ($N_Q$) e o número de passos de tempo ($N_T$). Especificamente, o número de testes de Hadamard necessários escala como $N_Q \cdot \text{poly}(N_G)$.
• Natureza de Pós-processamento: A técnica é inteiramente um método de pós-processamento, permitindo que seja acoplada a esquemas de mitigação existentes ou aplicada após a aquisição de dados sem modificar a execução do hardware quântico.
• Redução de Erro Controlável: O método introduz um parâmetro ajustável $r$ (o raio da hierarquia) que permite aos usuários controlar a intensidade da mitigação.

Resultos
O método foi testado numericamente em circuitos de preparação de estado Bell ajustáveis ($N_Q=2, N_G=2$) sob um modelo de ruído que mimetiza o processador IBM Fez.

• Redução Sistemática: Os resultados mostram uma redução sistemática do erro quântico ($\Delta' < \Delta$) em todos os parâmetros de circuito e raios de hierarquia testados.
• Ajustabilidade: A capacidade de redução de erro aumenta conforme o raio da hierarquia $r$ aumenta. Isso confirma que codificar uma porção maior da hierarquia BBGKY restringe a amostragem de forma mais eficaz em direção à dinâmica sem ruído.
• Autoconsistência: Em um raio específico ($r=4$ no caso testado), o sistema de equações torna-se completamente autoconsistente ($z=1$). Neste ponto, o método reduz-se a um cálculo numérico puramente clássico que recupera a dinâmica exata sem ruído, efetivamente ignorando a contribuição do ruído quântico.

Significância e Alegações
O artigo alega fornecer uma nova técnica de mitigação de erro quântico de propósito geral que é tanto eficiente (escalonamento polinomial) quanto controlável. Ao mapear circuitos para dinâmicas hamiltonianas e alavancar a amostragem informada pela física, o método oferece uma maneira sistemática de reduzir o ruído sem exigir recursos quânticos adicionais além da execução inicial do circuito e do pós-processamento padrão. Os autores enfatizam que a eficácia da técnica é ajustável via raio da hierarquia, tornando-a adaptável a diferentes níveis de ruído e restrições computacionais. O trabalho sugere potenciais aplicações futuras para algoritmos mais complexos (ex: Grover, Shor) e ansatzes variacionais mais profundos, observando sua compatibilidade como uma camada de pós-processamento com outras estratégias de mitigação.