Counterdiabatic ADAPT-VQE for molecular simulation

Resumo Técnico: ADAPT-VQE Contra-diabático para Simulação Molecular

Definição do Problema
Os Solucionadores Variacionais de Autovalores Quânticos (VQE) são algoritmos proeminentes para estimar estados fundamentais moleculares em dispositivos de Escala Intermediária com Ruído (NISQ). No entanto, os ansatzes fixos padrão, como aqueles baseados em Excitação de Cluster Unitário (UCCSD) com excitações de singlete e dupla, frequentemente sofrem do problema do "platô estéril" (barren plateau), onde a variância do gradiente da função de custo desaparece exponencialmente com o tamanho do sistema, dificultando o treinamento. Embora o ADAPT-VQE (VQE Adaptativo) mitigue isso construindo dinamicamente o ansatz a partir de um pool de operadores, a escolha do pool permanece crítica. Pools padrão baseados em excitações fermiônicas mapeadas para qubits podem resultar em circuitos profundos. Por outro lado, a otimização quântica contra-diabática digitalizada (DCQO) oferece vantagens de desempenho sobre a evolução adiabática padrão ao suprimir transições não-adiabáticas, mas frequentemente requer circuitos profundos quando aproximações de ordem superior do potencial de gauge adiabático (AGP) são utilizadas, uma vez que o número de termos cresce rapidamente com a ordem da aproximação.

Metodologia
Os autores propõem um algoritmo híbrido quântico-clássico, denominado CD-ADAPT, que integra a construção de ansatz adaptativo do ADAPT-VQE com os princípios de seleção de operadores da condução contra-diabática. A metodologia procede da seguinte forma:

1. Construção do Hamiltoniano: O Hamiltoniano molecular ($\hat{H}_f$) é mapeado para uma representação de qubits usando transformações padrão (ex: Jordan-Wigner). Um Hamiltoniano adiabático é construído interpolando entre um Hamiltoniano inicial ($\hat{H}_i$) e $\hat{H}_f$. O Hamiltoniano inicial é projetado de modo que seu estado fundamental corresponda ao estado de Hartree-Fock, utilizando operadores Pauli-Z locais.
2. Condução Contra-diabática e AGP: Para acelerar a evolução adiabática, um termo contra-diabático é adicionado, governado por um potencial de gauge adiabático (AGP), $\hat{A}_\lambda$. Como o AGP exato é computacionalmente intratável, os autores empregam uma aproximação baseada em comutadores aninhados de ordem $l$:
   $$\hat{A}^{(l)}_\lambda = i \sum_{k=1}^l \alpha_k(t) \hat{O}_{2k-1}(t)$$
   onde $\hat{O}_k$ envolve comutadores aninhados do Hamiltoniano adiabático e sua derivada em relação à função de cronograma (schedule function).
3. Geração do Pool de Operadores: Em vez de usar o operador de evolução temporal completo do protocolo contra-diabático (o que fixaria a profundidade do circuito), os autores extraem os termos geradores (strings de Pauli) a partir do AGP aproximado. Estes termos formam um novo pool de operadores $\{ \hat{G}_j \}$. O tamanho do pool depende da ordem de aproximação $l$ e, para $l \geq 2$, pode ser modulado pela avaliação dos coeficientes dependentes do tempo em tempos intermediários específicos $t'$.
4. Seleção Adaptativa: O algoritmo CD-ADAPT inicia com o estado de Hartree-Fock. Ele calcula iterativamente o gradiente de energia para cada operador no pool derivado do contra-diabático. O operador que produz a maior magnitude de gradiente é anexado ao ansatz, e uma otimização VQE é realizada para atualizar os parâmetros variacionais. Este processo se repete até que a norma do gradiente caia abaixo de um limiar $\epsilon$.

Principais Contribuições

• Estrutura Híbrida: O artigo introduz um novo método que utiliza termos de condução contra-diabática para definir o pool de operadores para um algoritmo ADAPT-VQE. Isso permite a inclusão de aproximações de AGP de ordem superior ($l > 1$) sem o aumento exponencial na profundidade do circuito tipicamente associado ao DCQO, uma vez que a seleção adaptativa filtra os termos menos relevantes.
• Modulação de Pool Dependente do Tempo: Os autores propõem uma técnica para gerar pools de operadores de tamanhos intermediários entre as ordens de aproximação $l=1$ e $l=2$ através da avaliação dos coeficientes dependentes do tempo do AGP em tempos específicos $t'$. Isso é alcançado usando o algoritmo de Decomposição de Pauli Tensorizada (TPD) para extrair strings de Pauli.
• Hamiltoniano Inicial Local: O trabalho define um Hamiltoniano inicial específico composto por operadores Pauli-Z que codifica o estado de Hartree-Fock como seu estado fundamental, facilitando o protocolo adiabático.

Resultados
O algoritmo foi simulado numericamente para três sistemas moleculares: Hidreto de Lítio (LiH), Fluoreto de Hidrogênio (HF) e Hidreto de Berílio de Linha (BeH$_2$), utilizando modelos de 10 qubits com transformações de espaço ativo.

• Precisão: O CD-ADAPT superou significativamente tanto o DCQO padrão quanto o ADAPT-VQE fermiônico.
  • Comparado ao DCQO (com 2 passos de Trotter), o CD-ADAPT reduziu o erro da energia do estado fundamental em aproximadamente três a cinco ordens de magnitude.
  • Comparado ao ADAPT-VQE fermiônico, o CD-ADAPT alcançou uma melhoria de cerca de uma a três ordens de magnitude em precisão.
  • Para aproximações $l=2$, os erros foram reduzidos para a faixa de $10^{-7}$ a $10^{-8}$ Hartree, bem abaixo da precisão química ($10^{-3}$ Ha).
• Eficiência de Circuito: Apesar de utilizar um pool de operadores potenciais maior, a seleção adaptativa resultou em circuitos mais compactos. O CD-ADAPT exigiu significativamente menos portas CNOT (aproximadamente 100–500 menos) do que o ADAPT-VQE fermiônico para as mesmas moléculas e limiares de convergência, mantendo uma precisão superior.
• Dependência de Ordem: O aumento da ordem de aproximação do AGP de $l=1$ para $l=2$ melhorou consistentemente a precisão. A escolha do tempo de avaliação $t'$ dentro do regime $l=2$ influenciou o desempenho, sendo que $t'=0.75$ geralmente produziu erros menores do que $t'=0.25$ para LiH e HF.

Significância e Alegações
Os autores alegam que combinar a condução contra-diabática com a construção de ansatz adaptativo do ADAPT-VQE fornece uma rota mais eficaz para estados fundamentais moleculares do que depender apenas de termos de interação motivados fisicamente, como excitações de singlete e dupla. A significância deste trabalho reside em demonstrar que a dinâmica adiabática acelerada pode guiar a seleção de operadores do ansatz, levando a:

1. Precisão Superior: Alcançando erros várias ordens de magnitude menores do que os benchmarks padrão.
2. Redução da Profundidade do Circuito: Gerando circuitos mais rasos comparados tanto ao DCQO padrão quanto ao ADAPT-VQE fermiônico, tornando a abordagem mais viável para dispositivos NISQ e de tolerância a falhas iniciais.
3. Escalabilidade de Aproximação: Permitindo o uso de aproximações de AGP de ordem superior ($l>1$) sem os custos proibitivos de profundidade de circuito normalmente associados a protocolos contra-diabáticos digitalizados.

O artigo conclui que esta abordagem híbrida é um candidato adequado para simulações moleculares em hardware quântico atual e futuro próximo, oferecendo uma alternativa robusta aos métodos variacionais e adiabáticos existentes.