⚛️ quantum physics

Quantum-Computing Self-Consistent Kohn-Sham DFT: Plane-Wave-Orthonormalized Orbitals with a Dual-Basis Quantum Eigensolver

Este artigo apresenta um método de DFT Kohn-Sham de todos os elétrons e de primeira quantização para computadores quânticos que utiliza uma base de orbitais ortonormalizados por ondas planas e um resolvedor de autovalores quântico de base dupla para alcançar cálculos de campo autoconsistente de alta precisão, com resultados validados através de emulação clássica e demonstrações de hardware em um dispositivo IonQ Forte-1.

Autores originais: Lazaro Calderin

Publicado 2026-09-29

📖 6 min de leitura🧠 Leitura aprofundada

CC BY 4.0

Autores originais: Lazaro Calderin

Artigo original sob licença CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Esta é uma explicação gerada por IA do artigo abaixo. Não foi escrita nem endossada pelos autores. Para precisão técnica, consulte o artigo original. Ler aviso legal completo

Para compreender o trabalho aqui descrito, deve-se primeiro entender o desafio fundamental de simular a matéria. No cerne da química e da ciência dos materiais reside a necessidade de calcular como os elétrons se organizam em torno dos núcleos atômicos. Esse arranjo determina tudo, desde a cor de uma pedra preciosa até a resistência de uma viga de aço. Por décadas, os cientistas confiaram em uma poderosa estrutura matemática chamada teoria do funcional da densidade para resolver esse problema em computadores convencionais. Ela funciona tratando a complexa dança dos elétrons como uma nuvem de densidade, permitindo que pesquisadores predigam a energia e o comportamento das moléculas com alta precisão. No entanto, à medida que os sistemas crescem em tamanho ou complexidade, os cálculos tornam-se tão exigentes que até mesmo os supercomputadores mais poderosos do mundo lutam para acompanhar. Essa limitação desencadeou a busca por um novo tipo de máquina: o computador quântico. Ao contrário das máquinas tradicionais que processam informações em bits binários, os computadores quânticos usam bits quânticos, ou qubits, que podem existir em múltiplos estados simultaneamente, oferecendo teoricamente uma maneira de simular sistemas quânticos, como moléculas, de forma direta e eficiente.

Os pesquisadores deste estudo, trabalhando na IonQ, deram um passo significativo para tornar essa visão realidade ao demonstrar uma simulação completa e autoconsistente de uma molécula em um computador quântico. Eles focaram em uma cadeia de oito átomos de hidrogênio, um sistema pequeno o suficiente para ser manejável, mas complexo o suficiente para testar os limites da tecnologia atual. O objetivo deles não era apenas executar um único cálculo, mas construir um ciclo completo onde o computador resolve o arranjo eletrônico, usa esse resultado para atualizar as forças que atuam sobre os elétrons e repete o processo até que o sistema se estabilize em um estado natural estável. Esse processo, conhecido como método do campo autoconsistente, é a maneira padrão pela qual os químicos resolvem esses problemas, mas realizá-lo inteiramente em um processador quântico nunca havia sido alcançado antes. A equipe executou com sucesso todo esse ciclo, provando que um computador quântico pode refinar iterativamente sua própria compreensão da energia de uma molécula até atingir uma solução precisa.

Um grande obstáculo nessas simulações é o número absoluto de funções matemáticas necessárias para descrever os elétrons, particularmente perto dos núcleos atômicos, onde seu comportamento muda rapidamente. Métodos tradicionais usando um único tipo de função matemática, conhecidos como ondas planas, exigiriam dezenas de milhares dessas funções para obter uma resposta precisa até para uma molécula pequena. Isso demandaria um computador quântico com muito mais qubits do que os existentes atualmente. Para superar isso, os pesquisadores desenvolveram uma abordagem híbrida inteligente. Eles combinaram as ondas planas padrão com um segundo conjunto de funções, chamadas orbitais atômicos, que são especificamente projetadas para descrever os elétrons próximos ao núcleo. Ao ortonormalizar esses dois conjuntos entre si — garantindo que sejam matematicamente independentes e não se sobreponham de forma confusa — eles criaram uma base mista que captura o detalhe necessário com muito menos funções no total. Isso permitiu que descrevessem o mesmo sistema físico com uma fração dos recursos computacionais, tornando a simulação viável na tecnologia atual.

O núcleo de seu método envolve uma estratégia de base dupla que atua como um filtro inteligente para a memória do computador quântico. Em vez de tentar resolver todo o problema complexo de uma vez em um grande registro quântico, o sistema primeiro realiza alguns passos iniciais usando o conjunto de funções completo e detalhado. A partir desses passos, ele identifica os estados mais importantes — aqueles que realmente importam para a energia da molécula — e os comprime em um conjunto muito menor e contraído. Para a vasta maioria do cálculo, o computador quântico trabalha apenas com este conjunto de funções reduzido e comprimido, que requer apenas três qubits para representação. Ele verifica periodicamente se essa visão comprimida ainda é precisa; se o erro crescer demais, ele retorna brevemente ao conjunto completo para atualizar a visão comprimida e, então, retorna ao cálculo eficiente de pequena escala. Essa técnica permitiu que eles simulassem um sistema com mais de 32.000 ondas planas usando apenas um registro de três qubits para o trabalho pesado, um feito que seria impossível sem essa compressão.

A equipe testou seu método de duas maneiras: simulando o comportamento do computador quântico em um supercomputador clássico e realizando experimentos reais em um processador quântico real. Nas simulações, mostraram que seu método podia reproduzir os resultados de cálculos clássicos com incrível precisão, diferindo por menos de um centésimo de um milihartree, uma unidade padrão de energia na química. Eles também demonstraram que o método funciona mesmo quando o computador quântico está sujeito ao "ruído" e aos erros inerentes ao hardware atual. Quando executaram a simulação com ruído de medição simulado, o sistema ainda conseguiu convergir para uma solução estável, pairando muito próximo do valor de energia correto. Isso provou que o algoritmo é robusto o suficiente para lidar com as imperfeições de dispositivos quânticos do mundo real.

Finalmente, os pesquisadores levaram sua simulação para o mundo real, executando as etapas finais do cálculo em um processador quântico IonQ Forte-1. Eles prepararam os circuitos quânticos que representam os estados eletrônicos convergidos e mediram suas energias diretamente no dispositivo. No pequeno registro de três qubits comprimido, o hardware reproduziu com sucesso a diferença de energia entre o método híbrido e um método padrão baseado apenas em átomos dentro da margem de erro experimental. Esta foi uma validação crítica, mostrando que o framework teórico se sustenta quando executado em hardware real. No entanto, quando tentaram medir o sistema completo, não comprimido, no registro maior, o ruído do hardware foi grande demais para obter um resultado preciso, mesmo após a aplicação de técnicas de correção de erro. Isso destacou um achado fundamental: embora o potencial total da computação quântica para a química seja imenso, os dispositivos atuais são mais adequados para esse tipo de cálculo quando o problema é inteligentemente comprimido e as partes mais críticas são medidas na menor escala possível.

O trabalho não afirma ter resolvido o problema da simulação química ou ter superado os computadores clássicos em velocidade. De fato, os pesquisadores observam que sua implementação atual é mais lenta que os métodos clássicos porque o hardware quântico ainda está em seus estágios iniciais. Na verdade, a conquista reside em provar o conceito: um computador quântico agora pode participar do ciclo iterativo completo de uma simulação química moderna, refinando suas próprias respostas até que sejam consistentes. Ao combinar uma base matemática inteligente com uma estratégia de compressão de base dupla, a equipe mostrou um caminho viável para executar cálculos de estrutura eletrônica autoconsistentes complexos em hardware quântico. Isso abre as portas para futuros experimentos onde computadores quânticos possam enfrentar problemas que estão atualmente fora de alcance, desde que os algoritmos continuem sendo projetados com as forças e limitações específicas do hardware em mente.

Resumo Técnico: DFT de Kohn–Sham de Autoconsistência para Computação Quântica

Enunciado do Problema
A teoria do funcional da densidade de Kohn–Sham (KS-DFT) é o padrão para simulações de estrutura eletrônica em computadores clássicos, mas sua aplicação à computação quântica enfrenta obstáculos significativos. As abordagens quânticas existentes tipicamente mapeiam Hamiltonianos de um único partícula de base fixa para qubits usando segunda quantização, frequentemente dependendo de solucionadores de autovalores variacionais (VQE) ou estimativa de fase quântica (QPE). Embora existam algoritmos de ondas planas (PW) de primeira quantização tolerantes a falhas, eles enfrentam dificuldades com o alto custo computacional de descrever sistemas de todos os elétrons, particularmente a convergência lenta das bases de ondas planas (PW) próximas aos núcleos. Além disso, demonstrações anteriores de campo autoconsistente (SCF) em hardware quântico foram limitadas a bases de orbitais atômicos (AO) pequenas e localizadas ou modelos de rede específicos, carecendo de uma formulação de potencial total e de todos os elétrons que possa convergir sistematicamente para tolerâncias de energia arbitrárias.

Metodologia
O autor apresenta um método de KS-DFT de potencial total, de todos os elétrons e de primeira quantização projetado para computadores quânticos. O núcleo da abordagem é um Solucionador de Autovalores de Circuito Quântico (QC eigensolver) que opera dentro de um loop autoconsistente. Componentes metodológicos principais incluem:

Bases de Funções:
- Onda Plana Pura (PW): Utiliza uma expansão sistemática onde o tamanho da base é aumentado até que a convergência de energia seja alcançada.
- Orbital Ortonormalizado por Onda Plana (PWOO): Uma base híbrida que combina um bloco retido de ondas planas de baixa energia com orbitais atômicos (AOs) que são ortonormalizados contra as ondas planas e entre si. Isso permite a representação de singularidades nucleares com AOs, reduzindo significativamente o número de ondas planas necessárias para uma determinada precisão de energia em comparação com a PW pura.
Contração de Base Dupla: Para mitigar o escalonamento exponencial dos requisitos de qubits, o método emprega uma estratégia de base dupla. Após alguns passos iniciais de "bootstrap" na base completa (PW ou PWOO), as iterações SCF são realizadas em uma base contraída de $K$ funções (exigindo apenas $\lceil \log_2 K \rceil$ qubits). A base contraída é atualizada apenas quando a norma do resíduo dos autovetores fora do subespaço excede um limiar específico. Isso mantém a largura do registro de qubits fixa enquanto permite que a base completa cresça.
Algoritmo do Solucionador de Autovalores: O QC eigensolver utiliza um ansatz de árvore de amplitude de Möttönen para codificar orbitais. A otimização é realizada via Rotosolve (descida de coordenadas), utilizando atualizações paralelas em bloco para eficiência. Soluções de múltiplos orbitais são obtidas via deflação por penalidade e ortonormalização de Gram–Schmidt modificada (MGSO). O algoritmo inclui um passo de Rayleigh–Ritz combinando candidatos de bloco e direções de resíduo computadas classicamente (em simulações de vetor de estado) para acelerar a convergência.
Loop SCF: O loop monta o Hamiltoniano de Kohn–Sham $H_{KS}[\rho]$ a partir da densidade gerada pelo solucionador de autovalores quântico, mistura a densidade (usando mistura linear ou mistura de Anderson) e itera até que a autoconsistência seja alcançada.

Resultos Principais
O artigo relata resultados para uma cadeia linear isolada de $H_8$ (espaçamento de 1.4 bohr) usando o funcional de densidade local SVWN.

Desempenho de Onda Plana Pura:
- Uma emulação de vetor de estado com $n_q = 15$ qubits (32.768 ondas planas) alcançou concordância com cálculos de PW clássicos dentro de 0,015 mHa.
- A contração de base dupla permitiu que esses cálculos de base grande rodassem em um registro contraído de 3 qubits para a maioria dos passos de SCF.
- A norma de Pauli $L_1$ do Hamiltoniano contraído permaneceu baixa (~0,6 Ha), enquanto a norma do Hamiltoniano completo cresceu para 233 Ha, reduzindo drasticamente o limite de disparos (shot bound) para a estimativa de energia.
Eficiência e Precisão da PWOO:
- A PWOO alcançou energias abaixo da referência STO-3G usando apenas 16 funções de base, enquanto a PW pura exigiu mais de 4.000 ondas planas para precisão comparável.
- Cálculos usando bases AO cc-pVDZ e cc-pVTZ satisfizeram o limite variacional de base aninhada (energia PWOO < energia apenas de AO), confirmando a montagem correta dos blocos de operadores.
- As energias de SCF quânticas usando a PWOO coincidiram com as energias de PWOO clássicas dentro de 0,44 µHa.
Simulações Baseadas em Disparos (Shot-Based):
- Execuções de ruído de disparo emuladas na base PWOO STO-3G ( $n_q=4$ ) demonstraram convergência para o ponto fixo sem ruído ao usar um parâmetro de mistura reduzido ( $\alpha=0.2$ ) e altos orçamentos de disparos ( $10^{13}$ disparos por lote).
- A mudança do projetor flutuou próximo ao limiar de convergência, com energias finais dentro de 0,05–0,31 mHa do resultado sem ruído.
Demonstração em Hardware (IonQ Forte-1):
- Um experimento em dispositivo real mediu os estados PWOO STO-3G convergidos em um registro contraído de 3 qubits.
- A diferença de energia medida entre os estados PWOO e apenas de AO foi de $-8,53 \pm 1,56$ mHa, comparado ao valor exato de $-7,06$ mHa. Este resultado situa-se dentro de um desvio padrão do valor exato.
- A leitura do registro completo (exigindo mais qubits e circuitos mais profundos) sofreu um viés significativo de hardware que não pôde ser totalmente corrigido por calibração de depolarização simples, destacando a vantagem da abordagem contraída para dispositivos de curto prazo.

Significância e Alegações
O artigo apresenta uma formulação de primeira quantização, de todos os elétrons e de potencial total de KS-DFT para computação quântica que utiliza uma base sistematicamente aperfeiçoável (PW e PWOO). Embora o autor reconheça uma demonstração molecular simulada anterior de Senjean et al. [21] (que usou uma base AO STO-3G e VQE de busca de subespaço), este trabalho distingue-se por implementar um loop autoconsistente que suporta cálculos de todos os elétrons e de potencial total com bases de ondas planas e PWOO sistematicamente convergentes.

Escalabilidade: O método demonstra que problemas de estrutura eletrônica de base grande (ex: 32.768 ondas planas) podem ser resolvidos em um computador quântico usando um registro de qubits pequeno e fixo (3 qubits) via contração de base dupla, evitando a necessidade de que a contagem de qubits escale com o tamanho da base.
Eficiência: A base PWOO reduz significamente o tamanho da base necessário para precisão química em comparação com as ondas planas puras, tornando-a mais adequada para o hardware quântico atual ou de curto prazo.
Validação: O limite variacional de base aninhada serve como uma verificação de consistência rigorosa para a montagem do operador na base mista.
Viabilidade de Hardware: A reprodução bem-sucedida de diferenças de energia no IonQ Forte-1 usando registros contraídos sugere que esta abordagem é viável para demonstrar vantagem quântica em estrutura eletrônica em dispositivos NISQ atuais, desde que o problema seja mapeado para um subespaço contraído suficientemente pequeno.

O autor afirma explicitamente que não reivindica uma vantagem de tempo de execução (wall-time) sobre o DFT clássico; o foco é a viabilidade algorítmica de realizar cálculos de estrutura eletrônica quântica autoconsistentes com recursos fixos de qubits. Trabalhos futuros estão planejados para abordar polarização de spin, aproximações de gradiente generalizado e funcionais híbridos.

Mais como este