Simulation of additive binding energies in asphalt using quantum-selected configuration interaction (QSCI)

Resumo Técnico: Simulação de Energias de Ligação Aditivas em Asfalto Usando Interação de Configuração Selecionada por Quantum (QSCI)

Declaração do Problema
A durabilidade de longo prazo e o envelhecimento oxidativo do ligante asfáltico, um componente crítico da infraestrutura rodoviária, são governados por interações moleculares complexas. Compreender como os aditivos interferem nas vias de envelhecimento exige o cálculo preciso das energias de ligação entre aditivos e componentes-chave do asfalto. No entanto, prever essas energias é desafiador devido à necessidade de contabilizar efeitos de correlação forte e interações de dispersão, que são frequentemente mal descritas por aproximações de campo médio padrão. Os autores visam demonstrar um fluxo de trabalho "centrado em quantum" capaz de calcular essas energias de ligação para um problema de materiais industrialmente relevante, focando especificamente na interação entre componentes do asfalto representados por um complexo piridina-fenol.

Metodologia
O estudo introduz o QuantumPave, um fluxo de trabalho híbrido quântico-clássico projetado para calcular energias de ligação de aditivos. A metodologia integra três abordagens complementares:

1. Seleção do Sistema: Os autores selecionaram um complexo piridina-fenol de 24 átomos ($C_{11}H_{11}NO$) como modelo representativo. Este sistema apresenta uma forte ligação de hidrogênio O-H$\cdots$N e sistemas eletrônicos $\pi$-conjugados, espelhando grupos funcionais identificados em estudos de dinâmica molecular do envelhecimento do asfalto.
2. Otimização Geométrica: Potenciais interatômicos de aprendizado de máquina (ORB v3) foram utilizados para otimização geométrica rápida, proporcionando um equilíbrio entre eficiência computacional e precisão química.
3. Cálculo de Energia de Ligação: Três métodos foram empregados para calcular energias de ligação:
   • Aprendizado de Máquina (ORB v3): Utilizado para triagem rápida e estimativas iniciais.
   • Teoria do Funcional da Densidade (DFT): Cálculos B3LYP-D3/def2-SVP serviram como referência clássica, incorporando correções de dispersão para empilhamento $\pi$-$\pi$ e ligação de hidrogênio.
   • Interação de Configuração Selecionada por Quantum (QSCI): Também conhecida como Diagonalização Quântica Baseada em Amostragem (SQD), este é o método central inspirado em quantum. Emprega uma aproximação de espaço ativo (10 elétrons em 10 orbitais, mapeados para 20 qubits) para reduzir a complexidade enquanto visa interações dominadas por valência.
4. Execução Quântica: O algoritmo QSCI foi executado no processador quântico supercondutor de 54 qubits IQM Emerald. O fluxo de trabalho envolve:
   • Amostragem: O processador quântico amostra configurações eletrônicas dominantes (strings de bits) a partir do estado preparado.
   • Mitigação de Erro: Em vez de extrapolação de ruído zero, o método utiliza recuperação de configuração autoconsistente. Este processo iterativo repara configurações inválidas (aquelas que violam a simetria do número de partículas) usando ocupações orbitais médias das amostras retidas.
   • Diagonalização: Os recursos clássicos de Computação de Alto Desempenho (HPC) realizam a diagonalização dentro do subespaço selecionado.
   • Correção de Dispersão: Correções empíricas D3BJ foram adicionadas aos resultados do QSCI para contabilizar interações de van der Waals de longo alcance não capturadas pela estrutura eletrônica do espaço ativo isoladamente.

Principais Resultados

• Precisão do Hardware: No processador IQM Emerald, o método SQD reproduziu exatamente a referência clássica de Interação de Configuração de Espaço Ativo Completo (CASCI) de espaço ativo. Ambos produziram uma energia de ligação de -3,52 kcal/mol (-0,153 eV).
• Robustez ao Ruído: O ruído do dispositivo, em vez de degradar o resultado, alargou a distribuição de amostragem para abranger o espaço ativo. Consequentemente, a diagonalização clássica recuperou a energia exata do espaço ativo sem exigir extrapolação de ruído zero.
• Comparação com o Experimento: O valor calculado do espaço ativo (-3,52 kcal/mol) subestima a entalpia calorimétrica experimental de aproximadamente -6,25 kcal/mol. Os autores atribuem essa discrepância à natureza compacta do espaço ativo, que captura correlações internas de valência, mas omite contribuições do espaço orbital mais amplo envolvido na ligação de hidrogênio eletrostática.
• Comparação com DFT: O método B3LYP-D3 superestimou a energia de ligação (-9,08 kcal/mol), consistente com tendências conhecidas de funcionais híbridos de superestimar ligações de hidrogênio moderadas.

Principais Contribuições

• Prova de Conceito para Supercomputação Centrada em Quantum: O artigo demonstra uma instanciação totalmente reproduzível e em pequena escala de supercomputação centrada em quantum, onde um QPU amostra configurações e o HPC clássico realiza a diagonalização.
• Validação de QSCI/SQD em Hardware: Fornece evidência experimental de que o QSCI pode alcançar precisão numérica contra seu limite clássico de espaço ativo em hardware quântico de escala intermediária ruidosa (NISQ) atual.
• Fluxo de Trabalho Tolerante a Ruído: O estudo destaca um fluxo de trabalho onde o ruído do dispositivo auxilia a convergência ao alargar a distribuição amostrada, eliminando a necessidade de técnicas complexas de extrapolação de ruído zero neste regime específico.
• Aplicação à Ciência dos Materiais: Aplica com sucesso esses algoritmos quânticos a um problema relevante para a indústria do asfalto, preenchendo a lacuna entre demonstrações de computação quântica e problemas industriais de materiais.

Significado e Alegações
O artigo alega que energias de ligação quimicamente significativas para problemas de materiais industrialmente relevantes são alcançáveis em hardware quântico atual dentro de um fluxo de trabalho centrado em quantum. O significado reside em demonstrar que:

1. Viabilidade: Processadores quânticos reais podem executar cálculos de estrutura eletrônica correlacionada que correspondem a referências clássicas de espaço ativo.
2. Praticidade: A abordagem é robusta contra ruído, sugerindo que estratégias de mitigação de erro, como extrapolação de ruído zero, podem não ser sempre necessárias se a estratégia de amostragem e o pós-processamento clássico forem otimizados.
3. Escalabilidade: Embora o espaço ativo atual seja limitado, a metodologia oferece um caminho sistemático para melhoria. Trabalhos futuros podem utilizar formulações eficientes em qubits e esquemas de incorporação (por exemplo, incorporação de matriz de densidade) para abordar ambientes de fase condensada maiores e mais realistas, como a matriz asfáltica.

Os autores concluem que a limitação primária é o tamanho do espaço ativo, e não o hardware em si, e que a abordagem fornece uma base promissora para aplicações de química computacional de próxima geração em ciência dos materiais e catálise.