Orbital-optimized spin-adapted multistate contracted VQE for excited states and properties on quantum hardware

Este artigo introduz o método VQE contratado de múltiplos estados com otimização orbital (oo-MC-VQE), que utiliza operadores adaptados ao spin para computar eficientemente estados fundamentais e excitados, juntamente com suas propriedades, em hardware quântico, enquanto equilibra precisão e complexidade de circuito através de uma escala de parâmetros linear com o número de estados.

O essencial

Imagine que você está tentando afinar uma orquestra massiva e complexa para tocar uma sinfonia perfeita. No mundo da química, a "orquestra" é uma molécula, e a "música" é a forma como seus elétrons se moveem e interagem. Para entender como uma molécula absorve luz (o que nos dá cores e impulsiona coisas como a fotossíntese), os cientistas precisam calcular as notas exatas que esses elétrons tocam.

Por muito tempo, fazer isso para moléculas com muitos elétrons foi como tentar resolver um quebra-cabeça que se torna exponencialmente mais difícil quanto mais peças você adiciona. Computadores clássicos (os que usamos hoje) eventualmente atingem um limite e não conseguem resolver esses quebra-cabeças para moléculas complexas.

Este artigo apresenta uma nova maneira de resolver esses quebra-cabeças usando computadores quânticos, que são máquinas especiais projetadas para lidar com esse tipo de complexidade naturalmente. Aqui está uma divisão simples do que os autores fizeram e descobriram:

1. O Problema: Afinando Muitas Notas ao Mesmo Tempo

Normalmente, os cientistas tentam afinar a orquestra para tocar apenas uma nota (o estado fundamental) perfeitamente. Mas, para entender como uma molécula reage à luz, precisamos conhecer muitas notas diferentes (estados excitados) ao mesmo tempo.

• O Desafio: Se você tentar afinar a orquestra para 10 músicas diferentes simultaneamente, as instruções (o circuito do computador) tornam-se incrivelmente longas e complicadas. Se as instruções forem muito longas, o computador quântico fica confuso devido ao "ruído" (estática ou erros), e a música desmorona.
• A Troca (Trade-off): Você precisa de um circuito complexo para obter uma resposta precisa, mas um circuito complexo tem mais chances de falhar em máquinas ruidosas atuais.

2. A Solução: Um Maestro Inteligente e Simétrico

Os autores desenvolveram um novo método chamado oo-MC-VQE. Pense nisso como um "maestro inteligente" para a orquestra quântica.

• Adaptado ao Spin (Spin-Adapted): Na química quântica, os elétrons têm uma propriedade chamada "spin" (como piões girando). Os autores construíram seu método de modo que o maestro sempre mantenha os piões girando da maneira simétrica correta. Isso evita que a música fique "fora de tom" devido a erros de simetria.
• Otimização de Orbitais (Orbital Optimized): Eles também permitiram que o maestro rearranjasse o mapa de assentos dos músicos (os orbitais) para fazer a música soar melhor antes mesmo de iniciar o ajuste complexo.
• Contração de Múltiplos Estados (Multistate Contracted): Em vez de tentar afinar 10 músicas com 10 manuais de instruções separados e massivos, eles encontraram uma maneira de usar um conjunto de instruções compartilhado e eficiente que funciona para todas as músicas ao mesmo tempo.

3. A Descoberta: Crescimento Linear

Uma das maiores questões era: Se eu quiser calcular 10 estados em vez de 1, precisarei de 10 vezes mais poder computacional?

• A Descoberta: Os autores descobriram que a resposta é surpreendentemente simples. A quantidade de "esforço" computacional (parâmetros do circuito) necessário cresce de forma linear. Se você dobrar o número de estados que deseja calcular, você aproximadamente apenas dobra o comprimento do manual de instruções. Isso é uma ótima notícia porque significa que o método é escalável.

4. O Teste no Mundo Real: Tocando em um Palco Ruidoso

Os autores não apenas simularam isso em um computador perfeito; eles realmente executaram seu método em hardware quântico real (computadores quânticos da IBM).

• A Configuração: Eles testaram duas pequenas moléculas: Formaldeído (um composto químico comum) e um cátion de tri-hidrogênio ($H_3^+$).
• A Questão do Ruído: Computadores quânticos reais são como um palco com uma multidão barulhenta e luzes piscando. Os resultados ficaram bagunçados sem ajuda.
• A Correção: Eles usaram técnicas de "mitigação de erro". Imagine isso como um engenheiro de som usando software para filtrar o ruído da multidão e as luzes piscando após a performance.
• O Resultado:
  • Para o Formaldeído, o método funcionou muito bem. Mesmo com o ruído, eles puderam ver claramente os "picos" no espectro de absorção (as cores que a molécula absorve).
  • Para o $H_3^+$, o ruído foi um problema maior, deslocando os resultados significativamente. Os autores observaram que isso ocorreu porque a matemática deste modelo específico é mais sensível ao ruído (como um instrumento delicado que desafina facilmente).
  • Lição Principal: Embora os números não tenham sido perfeitos nas máquinas reais, a forma dos resultados estava correta. Eles ainda puderam ver as principais características do comportamento da molécula.

Resumo

O artigo mostra que, ao usar uma abordagem inteligente e simétrica, os cientistas podem calcular o comportamento de elétrons excitados em moléculas usando computadores quânticos atuais e imperfeitos. Eles provaram que calcular múltiplos estados não exige uma quantidade impossível de recursos e que, com alguns truques de "cancelamento de ruído", é possível obter insights químicos úteis de dispositivos quânticos reais hoje.

O que eles NÃO alegaram:
O artigo não afirma que este método possa imediatamente projetar novos painéis solares, curar doenças ou criar novos materiais. Ele foca estritamente em provar que o método funciona para calcular espectros em hardware quântico. Quaisquer aplicações futuras são implícitas pelo campo de estudo geral, mas não são alegações específicas deste estudo.

Resumo técnico

Resumo Técnico: VQE de múltiplos estados contratado, com adaptação de spin e otimização de orbitais

Definição do Problema
A computação robusta e equilibrada de estados eletronicamente excitados e suas propriedades de transição continua sendo um desafio significativo, particularmente em regimes que envolvem reorganização eletrônica forte, quase-degenerescências ou caráter multiconfiguracional. Embora os métodos de resposta linear (LR) sejam bem-sucedidos em muitos contextos, eles tratam os estados excitados como perturbações de uma referência de estado fundamental, o que leva a limitações conceituais próximas a interseções cônicas e em sistemas fortemente correlacionados. Além disso, as abordagens atuais de algoritmo de eigensolver variacional quântico (VQE) para estados excitados, como o VQE de busca de subespaço (SS-VQE) e o VQE contratado de múltiplos estados (MC-VQE), frequentemente enfrentam um compromisso entre a precisão na descrição de um número crescente de estados eletrônicos e a expressividade (complexidade do circuito) do ansatz subjacente. Adicionalmente, o hardware quântico atual é limitado por ruído e decoerência, o que restringe a profundidade do circuito e exige o desenvolvimento de métodos que sejam tanto eficientes quanto compatíveis com técnicas de mitigação de erro.

Metodologia
Os autores introduzem o método VQE contratado de múltiplos estados com otimização de orbitais (oo-MC-VQE), que incorpora operadores adaptados ao spin para computar estados fundamentais e excitados, bem como propriedades específicas de estado e de transição.

• Adaptação de Spin: Para garantir que a simetria de spin das referências seja conservada durante toda a otimização do VQE, o método utiliza um pool de operadores adaptados ao spin. Especificamente, os autores empregam operadores de excitação simples e de excitação dupla de par adaptados ao spin singlete. Isso evita a necessidade de termos de penalidade de spin e garante a simetria de spin correta de todos os estados incluídos sem introduzir parâmetros redundantes.
• Estrutura Contratada de Múltiplos Estados: O método otimiza múltiplos estados eletrônicos simultaneamente usando um funcional de energia média de estado. Cada estado é construído usando um circuito de preparação de estado fixo e livre de parâmetros (inicializando bases ortogonais) seguido por um circuito correlacionador parametrizado (usando um ansatz de Produto de Unidade, ou UPS, que introduz correlação eletrônica). A otimização visa a energia média de estado, seguida por uma diagonalização clássica da matriz hamiltoniana de múltiplos estados para resolver as energias e funções de onda individuais.
• Otimização de Orbitais: A estrutura é estendida para incluir a otimização de orbitais (oo-MC-VQE). Rotações orbitais são tratadas como uma função de custo clássica, permitindo a otimização simultânea dos parâmetros do circuito ($\theta$) e dos parâmetros de rotação orbital ($\kappa$). Os autores demonstram que, para o ansatz específico utilizado, os parâmetros de circuito de excitação simples adjacentes à parametrização orbital são redundantes e podem ser removidos, reduzindo a complexidade do circuito.
• Avaliação de Elementos de Matriz: Para avaliar elementos de matriz entre diferentes estados (necessários para o hamiltoniano de múltiplos estados e propriedades de transição), os autores empregam uma técnica de decomposição de superposição. Esta abordagem evita circuitos de preparação de estado profundos para estados multi-determinantes ao avaliar superposições mais simples de dois determinantes, reduzindo assim o número de portas CNOT ao custo de um aumento nas medições de valor esperado.
• Mitigação de Erro: A implementação integra o Mansatz0, uma técnica de mitigação de erro de leitura e de porta baseada em ansatz. Para cálculos de múltiplos estados envolvendo superposições de estados de referência, os autores propõem uma extensão (M0+) que caracteriza erros para todos os estados de referência multi-determinantes únicos, embora notem a escala exponencial desta abordagem.

Principais Contribuições

1. oo-MC-VQE Adaptado ao Spin: O desenvolvimento de uma formulação adaptada ao spin para VQE de média de estado com otimização de orbitais, garantindo a conservação da simetria de spin correta em múltiplos estados.
2. Análise de Escala: Uma investigação explícita sobre o compromisso entre precisão e complexidade de circuito. Os autores descobrem que o número de parâmetros de circuito necessários para obter resultados precisos escala aproximadamente de forma linear com o número de estados incluídos na média de estado.
3. Implementação em Hardware Quântico: Uma implementação explícita de circuito quântico do método, incluindo a redução de portas CNOT via superposição de determinantes e a integração de técnicas de mitigação de erro.
4. Demonstração em Dispositivo Real: Execução do método em hardware quântico real (IBM ibm marrakesh e ibm aachen) para computar espectros de absorção para sistemas moleculares de referência (formaldeído e $H_3^+$).

Resultados

• Complexidade de Circuito: Simulações em LiH, $H_2O$ e $NH_3$ usando os ansatzes ADAPT e oo-ADAPT mostram que o número de parâmetros necessários aumenta aproximadamente de forma linear com o número de estados. Por exemplo, alcançar um erro abaixo de $10^{-4}$ Hartree para $H_2O$ exigiu 27, 32, 24 e 30 parâmetros por estado para 1, 2, 3 e 4 estados, respectivamente.
• Impacto da Otimização de Orbitais: Embora a otimização de orbitais melhore a convergência nas fases iniciais do procedimento ADAPT, ela tem um efeito negligenciável no número final de parâmetros de circuito para uma convergência mais rigorosa. Em alguns casos (ex: $NH_3$ com 3 ou 4 estados), dobrar excitações simples na otimização orbital teve um efeito deletério na contagem de parâmetros.
• Desempenho de Hardware Quântico:
  • Ruído Simulado: Em hardware emulado usando modelos de ruído da IBM, as técnicas de mitigação de erro (M0 e M0+) produziram resultados quase idênticos às simulações sem ruído. O sistema $H_3^+$ foi mais adversamente afetado pelo ruído do que o formaldeído, atribuído a um número de condição mais alto da matriz hamiltoniana de múltiplos estados.
  • Hardware Real: Cálculos em dispositivos reais mostraram que, embora os resultados mitigados pelo erro preservassem as tendências qualitativas e as principais características dos espectros de absorção, a precisão quantitativa era significativamente menor do que nos dispositivos simulados. A discrepância destaca as limitações dos modelos de ruído atuais em relação ao ruído dependente do tempo e ao crosstalk. No entanto, o método conseguiu extrair insights químicos significativos apesar do ruído do hardware.

Significância e Alegações
O artigo afirma que o método oo-MC-VQE fornece uma descrição equilibrada da estrutura eletrônica adequada para cenários desafiadores como interseções cônicas, abordando as limitações inerentes às abordagens de resposta linear. Os autores enfatizam que a escala linear de parâmetros de circuito com o número de estados sugere que o método é viável para visar múltiplos estados excitados em hardware quântico de curto prazo (near-term). Embora reconheçam que a precisão quantitativa permanece desafiadora nos dispositivos atuais, o trabalho demonstra que insights químicos significativos, como espectros de absorção, podem ser extraídos quando combinando o método com técnicas apropriadas de mitigação de erro. O estudo serve como uma prova de conceito para a execução de cálculos de múltiplos estados, adaptados ao spin e com otimização de orbitais em dispositivos quânticos reais.