Tailored coupled cluster method with sample-based quantum diagonalization: Application to titanium-based metallocene catalytic reactions for 1-hexene production

Resumo Técnico: Método de Coupled Cluster Customizado com Diagonalização Quântica Baseada em Amostragem

Definição do Problema
A previsão precisa da seletividade em reações catalíticas de metalocenos à base de titânio para a produção de 1-hexeno requer o cálculo de diferenças de energia livre entre estados de transição concorrentes (inserção de etileno vs. transferência de beta-hidreto) com uma precisão dentro de 1 kcal/mol. Os métodos atuais de estrutura eletrônica enfrentam um dilema: a Teoria do Funcional da Densidade (DFT) é computacionalmente eficiente, mas frequentemente carece da precisão necessária devido à dependência do funcional, enquanto métodos de função de onda de alto nível (ex: CCSD(T)) são computacionalmente proibitivos para sistemas de metais de transição realistas. Além disso, abordagens de computação quântica como a Diagonalização Quântica Baseada em Amostragem (SQD) são atualmente limitadas por restrições de hardware a espaços ativos que são pequenos demais para capturar toda a correlação dinâmica necessária para resultados quimicamente precisos em sistemas catalíticos. Consequentemente, tratar tanto a correlação estática (dentro do espaço ativo) quanto a correlação dinâmica (de orbitais externos) simultaneamente permanece um desafio significativo.

Metodologia
Este estudo propõe uma estrutura híbrida quântico-clássica que integra a Diagonalização Quântica Baseada em Amostragem (SQD) com a teoria de Tailored Coupled Cluster (TCC) para abordar as limitações do tamanho do espaço ativo e do tratamento da correlação.

1. Sistema Alvo: Um modelo molecular mínimo de um catalisador de titânio ansa-ciclopentadienil-areno catiônico foi utilizado para investigar as energias relativas de dois estados de transição, T3 (inserção de etileno) e T4 (transferência de beta-hidreto). Para garantir um arcabouço orbital consistente, o estado T3 foi aumentado com uma molécula de etileno explícita (T3(+C2H4)), igualando a contagem de átomos e elétrons de T4.
2. SQD para Correlação Estática: A SQD foi empregada para tratar a estrutura eletrônica do espaço ativo. Utilizando o ansatz LUCJ (Local-Unitary Cluster Jastrow) no processador quântico IBM Heron, a SQD amostra estocasticamente bitstrings para construir um Hamiltoniano projetado dentro de um subespaço. Isso permite o tratamento de espaços ativos (até 30 orbitais espaciais) que são combinatorialmente muito grandes para a Função de Configuração de Espaço Ativo Completo (CASCI) clássica. O método utiliza recuperação de configuração autoconsistente para mitigar o ruído do hardware e garantir números de elétrons e spin corretos.
3. TCC para Correlação Dinâmica: As funções de onda derivadas da SQD serviram como estados de referência para cálculos de Tailored Coupled Cluster (TCC). Nesta estrutura, as amplitudes de cluster para o espaço ativo ($\hat{T}_{AS}$) são "customizadas" (fixadas) com base nos resultados da SQD, capturando a correlação estática e excitações de ordem superior (triplos, quadruplos) dentro do espaço ativo. As amplitudes restantes para o espaço externo ($\hat{T}_{Ext}$) são otimizadas via equações padrão de Coupled Cluster Singles and Doubles (CCSD) para capturar a correlação dinâmica.
4. Extensões: O estudo aplicou adicionalmente uma extensão de triplos perturbativos, TCC(T), para melhorar a precisão. Um esquema de correção de energia foi utilizado para garantir que a energia do espaço ativo reflita o nível da SQD em vez da aproximação de nível inferior do CCSD.
5. Configuração Computacional: Todos os cálculos de ponto único pós-DFT foram realizados na fase gasosa para isolar os efeitos de correlação eletrônica, embora as otimizações de geometria tenham incluído efeitos de solvente (tolueno) via Modelo de Continuidade Polarizável (PCM).

Principais Resultados

• Problemas de Convergência com SQD Isolada: Quando a SQD foi aplicada isoladamente, as energias relativas dos estados de transição não puderam ser convergidas dentro dos espaços ativos atualmente viáveis. Os resultados foram sensíveis ao tamanho do espaço ativo, indicando que a correlação estática sozinha é insuficiente para este sistema.
• Sucesso do SQD-TCC: A abordagem híbrida SQD-TCC yielded resultados razoavelmente convergidos para as energias relativas. Ao incorporar a correlação dinâmica via estrutura TCC, o método equilibrou com sucesso o tratamento da correlação eletrônica.
• Impacto da Correlação Estática: Uma discrepância significativa (excedendo 1 kcal/mol) foi observada entre os resultados do CCSD(T) padrão e os do SQD-TCC(T). Isso destaca que os efeitos de correlação estática são críticos neste sistema de titânio e que negligenciá-los (como no CCSD(T) de referência única padrão) leva a erros significativos nas previsões de seletividade.
• Escalabilidade do Espaço Ativo: O estudo demonstrou que a SQD permite o tratamento de espaços ativos (ex: 30 orbitais) que são praticamente inacessíveis aos métodos CASCI clássicos dados os recursos computacionais disponíveis (até 192 núcleos de CPU e 4 TiB de memória).

Significância e Alegações
O artigo afirma que a integração de funções de onda multirreferência derivadas de sistemas quânticos com a teoria de coupled cluster clássica fornece um caminho viável para simulações quimicamente precisas de sistemas catalíticos complexos. Especificamente, o trabalho demonstra que:

1. Necessidade Híbrida: Tanto as correlações eletrônicas estáticas quanto as dinâmicas são essenciais para descrições confiáveis da química subjacente na catálise de metalocenos de titânio.
2. Avanço Metodológico: A estrutura SQD-TCC preenche com sucesso a lacuna entre as capacidades atuais de computação quântica (espaços ativos limitados) e os requisitos para simulações catalíticas realistas (grandes espaços ativos + correlação dinâmica).
3. Previsão de Seletividade: O tratamento da correlação eletrônica influencia significativamente a seletividade prevista da reação, com o método híbrido oferecendo uma abordagem mais robusta do que os métodos clássicos de alto nível para sistemas com forte correlação estática.

Os autores concluem que, embora o presente estudo foque em superfícies de energia potencial na fase gasosa, a metodologia estabelece uma base para trabalhos futuros que incorporem contribuições de energia livre e efeitos de solvente ao nível da função de onda correlacionada.