On the generalized eigenvalue problem in subspace-based excited state methods for quantum computers

O artigo demonstra que métodos de estado excitado baseados em expansão de subespaço, como QSE e qEOM, tornam-se instáveis ou insolúveis em dispositivos quânticos devido à amplificação de erros estatísticos em matrizes de sobreposição mal condicionadas, sugerindo que abordagens como q-sc-EOM, que utilizam equações de autovalor padrão, são alternativas mais estáveis e adequadas para cálculos de química quântica.

O essencial

Imagine que você é um chef de cozinha tentando descobrir a receita perfeita para um prato novo (o estado excitado de uma molécula). Você tem uma cozinha super avançada (o computador quântico), mas ela tem um pequeno defeito: às vezes, quando você pede para medir o peso de um ingrediente, a balança treme um pouco e dá um número levemente errado. Isso é o "ruído de amostragem" mencionado no artigo.

O problema que os autores deste estudo estão resolvendo é: como encontrar a receita certa mesmo quando a balança não é 100% precisa?

Eles compararam duas formas diferentes de tentar descobrir essa receita: o método antigo (chamado QSE) e um método mais novo e robusto (chamado q-sc-EOM).

Aqui está a explicação do que eles descobriram, usando analogias do dia a dia:

1. O Problema da "Balança Trêmula" (O Método QSE)

O método QSE funciona como se você estivesse tentando resolver um quebra-cabeça onde as peças precisam se encaixar perfeitamente em uma moldura.

• A Analogia: Imagine que você tem uma moldura (a "matriz de sobreposição") que deve ser perfeitamente plana para você colocar as peças do quebra-cabeça (os dados da molécula).
• O que acontece: Quando a balança treme (erro de medição), a moldura fica levemente torta.
  • Se a moldura estiver quase reta (condição de número baixa), você ainda consegue montar o quebra-cabeça, mesmo com um pouco de tremedeira.
  • Mas, se a moldura já estiver um pouco torta (condição de número média), qualquer tremedeira na balança faz a moldura ficar extremamente torta. O erro se amplifica! Você tenta colocar as peças, mas elas não encaixam, ou você acaba com um prato que não faz sentido.
  • O Pior Cenário: Se a moldura estiver muito torta (condição de número alta), o quebra-cabeça se torna impossível de montar. A matemática "quebra" e você não consegue nenhuma receita. Para consertar isso, você é obrigado a jogar fora algumas peças do quebra-cabeça (uma técnica chamada "limiar" ou thresholding). O resultado? Você consegue montar algo, mas falta metade das peças. No mundo da química, isso significa que você descobre alguns estados de energia, mas esquece completamente de outros estados importantes. É como se você descobrisse a receita do bolo, mas esquecesse que precisava de chocolate.

2. A Solução Inteligente (O Método q-sc-EOM)

O método q-sc-EOM é como se você tivesse uma moldura feita de um material indestrutível e perfeitamente plano.

• A Analogia: Neste método, a "moldura" é construída de tal forma que ela é sempre perfeita (matematicamente, é uma "matriz identidade"). Não importa o quanto a balança treme ao pesar os ingredientes, a moldura nunca fica torta.
• O Resultado: Como a moldura é estável, você não precisa fazer aquele passo perigoso de tentar "endireitar" a moldura (inverter a matriz) antes de montar o quebra-cabeça.
• Vantagem: Mesmo com a balança tremendo, você consegue montar o quebra-cabeça completo. Você encontra todos os estados de energia, sem perder nenhum, e com muito menos erros do que o método antigo.

3. O Que Significa Isso para o Futuro?

Os autores testaram isso em moléculas reais (como um grupo de 4 átomos de hidrogênio e amônia) e provaram matematicamente que:

• O método antigo (QSE) é muito sensível. Em situações onde a molécula é complexa, ele exige que você faça milhões de medições (milhões de "tiro" ou shots) apenas para tentar compensar o erro e chegar perto da verdade. E mesmo assim, se a molécula for muito difícil, você perde dados importantes.
• O método novo (q-sc-EOM) é muito mais eficiente. Ele precisa de muito menos medições para dar um resultado preciso e, o mais importante, não perde nenhum estado excitado.

Resumo em uma Frase

Enquanto o método antigo tenta adivinhar a receita tentando corrigir uma balança quebrada (o que faz os erros explodirem e esconderem ingredientes), o método novo usa uma balança que é matematicamente perfeita, garantindo que você descubra a receita completa sem perder nenhum detalhe, mesmo com a cozinha um pouco barulhenta.

Conclusão: Para os computadores quânticos atuais (que ainda têm um pouco de ruído), o método q-sc-EOM é a escolha mais segura e confiável para estudar estados excitados de moléculas, evitando que cientistas "esqueçam" partes cruciais da química das moléculas.

Resumo técnico

Aqui está um resumo técnico detalhado do artigo "On the generalized eigenvalue problem in subspace-based excited state methods for quantum computers", apresentado em português:

Título: Sobre o problema de autovalores generalizado em métodos de estados excitados baseados em subespaço para computadores quânticos

1. O Problema

O artigo aborda um desafio crítico na química quântica em computadores quânticos de pré-falta de tolerância a falhas: a instabilidade numérica de métodos baseados em subespaço para calcular estados excitados.

• Contexto: Algoritmos como Quantum Subspace Expansion (QSE) e Quantum Equation of Motion (qEOM) requerem a resolução de uma equação de autovalores generalizada ($HC = SCE$) em um computador clássico, onde os elementos das matrizes de Hamiltoniano ($H$) e de sobreposição ($S$) são medidos em um computador quântico.
• A Causa da Instabilidade: Em computadores quânticos, as medições estão sujeitas a erros de amostragem estatística (ruído de "shots"). A resolução da equação generalizada exige a inversão da matriz de sobreposição $S$. Quando o número de condição de $S$ é alto (o que ocorre frequentemente em geometrias moleculares específicas ou devido a ruído), a inversão torna-se extremamente sensível a erros.
• Consequência: Pequenos erros estatísticos nas medições quânticas são amplificados drasticamente durante a inversão de $S$, levando a autovalores (energias) imprecisos, variâncias elevadas ou até singularidades (matrizes não invertíveis), tornando o método QSE instável e, em casos extremos, insolúvel sem técnicas de estabilização que comprometem a completude espectral.

2. Metodologia

Os autores realizaram uma análise combinada teórica (analítica) e numérica para comparar dois métodos:

1. QSE (Quantum Subspace Expansion): Utiliza uma equação de autovalores generalizada ($HC = SCE$), onde $S$ é uma matriz de sobreposição não trivial que deve ser invertida.
2. q-sc-EOM (Quantum Self-consistent Equation of Motion): Utiliza uma formulação onde o conjunto de vetores de base é ortonormal por construção analítica. Isso resulta em uma matriz de sobreposição que é identicamente a matriz identidade ($S = I$), transformando o problema em uma equação de autovalores padrão ($HC = CE$), eliminando a necessidade de inverter $S$.

Configuração Experimental:

• Sistemas Estudados: Molécula de $H_4$ (com diferentes geometrias lineares e quadradas) e amônia ($NH_3$) plana.
• Simulação: Utilização de bases STO-3G e STO-6G, mapeamento Jordan-Wigner e ansatz UCCSD via VQE para o estado fundamental.
• Variáveis: Foram introduzidos ruídos de amostragem (shot noise) simulando um número finito de medições. Analisaram-se cenários de baixo, médio e alto número de condição da matriz $S$ no QSE.
• Ferramentas: Códigos desenvolvidos em PennyLane, PySCF e Qiskit.

3. Contribuições Principais

• Análise Teórica de Erros: Derivação de limites de perturbação que quantificam como o erro nos autovalores escala com o número de condição da matriz de sobreposição ($\kappa(S)$). Demonstraram que o erro escala com $1/\lambda_{min}(S)$, onde$\lambda_{min}(S)$é o menor autovalor de$S$.
• Diagnóstico Prático: Estabelecimento de diretrizes práticas ligando o número de condição, o ruído de shot e a confiabilidade dos estados excitados. Mostraram que, para manter a precisão, o número de shots necessário cresce quadraticamente à medida que $\lambda_{min}(S)$ diminui.
• Comparação de Estabilidade: Evidência numérica robusta de que o q-sc-EOM, ao evitar a inversão de $S$, é intrinsecamente mais robusto a erros de amostragem do que o QSE.
• Análise de Técnicas de Estabilização: Investigação do uso de thresholding (corte de autovalores pequenos de $S$) para resolver singularidades no QSE, demonstrando que, embora resolva a singularidade, isso leva à perda de estados excitados (incompletude espectral).

4. Resultados

• Cenário de Baixo Número de Condição ($\kappa(S) < 100$): O QSE e o q-sc-EOM apresentam desempenho comparável e convergem de forma semelhante com o aumento do número de shots. A amplificação de erro não é dominante.
• Cenário de Médio Número de Condição ($100 < \kappa(S) < 10^4$):
  • O QSE torna-se instável. Erros na média e na variância das energias aumentam drasticamente com o aumento do número de condição. Para atingir a mesma precisão do q-sc-EOM, o QSE exigiria uma ordem de magnitude mais de shots (ex: $10^5$vs$10^4$).
  • O q-sc-EOM mantém-se estável e preciso, pois não sofre com a amplificação de erro associada à inversão de $S$.
  • Verificação: Quando a matriz $S$ do QSE é artificialmente tornada exata (sem ruído), os erros do QSE diminuem drasticamente, confirmando que o ruído na inversão de $S$ é a fonte principal da instabilidade.
• Cenário de Alto Número de Condição ($\kappa(S) > 10^4$):
  • O QSE torna-se singular e insolúvel com ruído de shot padrão.
  • A aplicação de thresholding permite resolver a equação, mas resulta na perda de estados excitados (especialmente estados degenerados), comprometendo a completude espectral necessária para estudos químicos.
  • O q-sc-EOM continua a fornecer o conjunto completo de estados excitados com estabilidade.
• Caso $NH_3$: A tendência observada no $H_4$ foi confirmada na amônia, onde aumentos moderados no número de condição (devido ao estiramento da ligação) causaram grandes desvios no QSE, enquanto o q-sc-EOM permaneceu preciso.

5. Significância e Conclusão

O estudo conclui que métodos de estados excitados baseados em subespaço que dependem da resolução de equações de autovalores generalizadas (como QSE e qEOM padrão) são vulneráveis a erros de amostragem estatística inerentes aos computadores quânticos atuais. A instabilidade é exacerbada pelo número de condição da matriz de sobreposição.

Em contraste, o q-sc-EOM é identificado como um candidato superior para cálculos de química quântica em estados excitados em dispositivos de pré-falta de tolerância a falhas. Sua formulação com matriz de sobreposição identidade ($S=I$) elimina o mecanismo de amplificação de erro associado à inversão de matrizes, oferecendo:

1. Maior estabilidade numérica frente ao ruído.
2. Redução significativa nos requisitos de recursos (número de shots) para atingir uma dada precisão.
3. Completude espectral garantida, sem a necessidade de técnicas de corte (thresholding) que ocultam estados físicos.

O trabalho alerta que problemas semelhantes podem afetar outros algoritmos quânticos que utilizam equações generalizadas (como métodos de Krylov quântico) e sugere que a escolha de formulações que evitam a inversão de matrizes de sobreposição é crucial para a viabilidade prática de simulações químicas em hardware quântico atual.