Correlation-consistent Gaussian basis sets for copper solids from material-constrained atomic optimization

Resumo Técnico: Conjuntos de Bases Gaussianas Consistentes com a Correlação para Sólidos de Cobre a partir da Otimização Atômica com Restrição de Material

Enunciado do Problema
Conjuntos de bases Gaussianas consistentes com a correlação (por exemplo, a série cc-pVXZ de Dunning) são o padrão para a química quântica molecular sistemática, permitindo a extrapolação confiável para o limite de base completa (CBS). No entanto, sua aplicação direta a sólidos periódicos, particularmente sistemas metálicos ou que contêm metais, é severamente limitada pela dependência linear. Este problema surge porque a otimização atômica frequentemente gera funções gaussianas primitivas difusas (com expoentes muito pequenos, $\alpha \ll 0,1$ u.a.) que se tornam quase linearmente dependentes quando os átomos são aproximados em uma rede sólida. Isso leva a números de condição ($\kappa$) elevados na matriz de sobreposição, causando instabilidades numéricas, como convergência difícil do campo autoconsistente (SCF) e superfícies de energia potencial não suaves.

Remédios existentes, como a truncagem de expoentes pequenos ou a realização de otimizações específicas para o sistema, frequentemente comprometem a transferibilidade dos conjuntos de bases ou destroem a estrutura consistente com a correlação necessária para cálculos de função de onda correlacionada de alta precisão. Embora alguns esforços recentes (por exemplo, GTH-cc-pVXZ) tenham preservado a consistência com a correlação ao limitar o tamanho do conjunto de valência, a natureza discreta de tais protocolos de seleção levanta questões sobre o equilíbrio ideal entre precisão e estabilidade.

Metodologia: Otimização Atômica com Restrição de Material (MCAO)
Para abordar esses desafios, os autores introduzem a Otimização Atômica com Restrição de Material (MCAO). Este framework aumenta a otimização padrão de conjuntos de bases atômicas com uma restrição de estabilidade dependente do material. O núcleo da MCAO é uma função de perda projetada para minimizar a energia atômica enquanto penaliza grandes números de condição em sólidos representativos:

$$L(\alpha; \{M\}, \kappa_0, \varepsilon_0) = E_{\text{átomo}}(\alpha) + \frac{\varepsilon_0}{\kappa_0} \times \max_{M \in \{M\}} \kappa(\alpha; M)$$

Onde:

• $E_{\text{átomo}}(\alpha)$ é a energia atômica padrão (HF para valência, MP2 para polarização).
• $\kappa(\alpha; M)$ é o número de condição da matriz de sobreposição para um sólido de referência $M$ (neste trabalho, cobre (Cu) de estrutura cúbica de faces centradas).
• $\kappa_0$ e $\varepsilon_0$ são parâmetros selecionados pelo usuário que controlam o número de condição alvo e a escala de energia da penalidade, respectivamente.

Ao contrário da truncagem de expoentes, a MCAO determina a compacidade apropriada através de otimização contínua em vez de cortes rígidos. O procedimento preserva o protocolo padrão de otimização atômica tanto quanto possível, garantindo que os conjuntos de bases resultantes mantenham a estrutura consistente com a correlação essencial para cálculos correlacionados.

Contribuições Principais e Implementação
Como prova de conceito, os autores aplicaram a MCAO para gerar conjuntos de bases consistentes com a correlação do tipo Dunning para o Cobre (Cu), denominados MCAO-cc-pVXZ ($X = \text{D, T, Q}$). O estudo abrange uma gama abrangente de tratamentos de potencial nuclear e relativísticos:

1. Todo-elétron (AE): Não-relativístico (NR) e escalar-relativístico (SFX2C-1e).
2. Potenciais de Core Efetivos (ECP): ccECPs duros (hard) e suaves (soft).
3. Pseudopotenciais: Pseudopotenciais GTH-PBE de núcleo pequeno (small-core, sc) e núcleo grande (large-core, lc).

A otimização visa a configuração [Ar]$3d^{10}4s^1$. O conjunto de valência é otimizado no nível HF, e o conjunto de polarização no nível MP2, correlacionando os 11 elétrons de valência. O esquema de contração utiliza orbitais naturais atômicos (ANOs) de HF com média esférica, com descontração parcial para formar conjuntos DZ, TZ e QZ.

Resultados
Os conjuntos de bases MCAO-cc-pVXZ demonstram um equilíbrio bem-sucedido entre estabilidade numérica e precisão:

• Estabilidade Numérica: Os números de condição ($\kappa$) para o Cu fcc são reduzidos significativamente em comparação com conjuntos otimizados atomicamente padrão (por exemplo, cc-pVXZ, def2-XZVP), trazendo-os abaixo do limiar de instabilidade de $10^9$. Os conjuntos de bases permanecem estáveis para sólidos e superfícies de Cu sem a necessidade de truncagem ad hoc.
• Precisão Molecular: Para o dímero de Cu, os conjuntos MCAO reproduzem as energias de ligação dentro de 0,02 eV/átomo dos conjuntos cc-pVXZ padrão em todos os níveis zeta e no limite de CBS. Em contraste, conjuntos de bases gerados via simples truncagem de expoentes (por exemplo, def2-TZVP modificado) mostraram sobreestimar as energias de ligação, indicando uma descrição menos equilibrada.
• Propriedades de Estado Sólido: Para o Cu bulk, os conjuntos MCAO convergem rapidamente as constantes de rede ($a_0$) e os módulos de compressão ($B_0$) para os valores de referência de ondas planas (PW) (dentro de 0,01 Å e 1 GPa, respectivamente). As estruturas de bandas, incluindo as bandas de condução de baixa energia críticas para cálculos correlacionados, são reproduzidas com precisão.
• Avaliação de Erro no Cu Bulk: Usando cálculos de Aproximação de Campo de Fase Aleatória extrapolados para CBS (RPA@PBE), os autores avaliaram os erros provenientes de diferentes aproximações. O pseudopotencial GTH-PBE-sc e os ccECPs suaves concordaram bem com a referência de todo-elétron SFX2C-1e. No entanto, o pseudopotencial GTH-PBE-lc sobreestimou substancialmente a força de interação (subestimando $a_0$ em ~0,1 Å e sobreestimando $B_0$ em ~70 GPa), destacando os desafios na parametrização de pseudopotenciais para métodos correlacionados.
• Adsorção de CO em Cu(111): Os conjuntos de bases foram aplicados ao "problema da adsorção de CO", onde o DFT semilocal prevê incorretamente uma preferência de sítio hollow sobre o sítio top observado experimentalmente. Cálculos RPA@PBE extrapolados para CBS confirmaram a preferência pelo sítio top. O estudo forneceu benchmarks de todo-elétron escalar-relativísticos de alta qualidade, mostrando que a maioria dos pseudopotenciais e ECPs reproduz a energia de preferência de sítio dentro de 0,02 eV da referência de todo-elétron, enquanto o GTH-PBE-lc falhou significativamente.

Significância e Alegações
O artigo afirma que a MCAO fornece uma rota robusta e sistemática para gerar conjuntos de bases Gaussianas para sólidos que mantêm a precisão e a estrutura consistente com a correlação da otimização atômica, enquanto controla a dependência linear. Os autores argumentam que:

1. Estabilidade sem Sacrifício: É possível alcançar estabilidade numérica em sólidos periódicos sem abandonar o framework consistente com a correlação ou recorrer a otimizações específicas de sistema que comprometem a transferibilidade.
2. Capacidade de Benchmarking: A extrapolação CBS confiável possibilitada por estes conjuntos de bases permite o isolamento controlado e a avaliação de erros decorrentes de pseudopotenciais, ECPs e aproximações relativísticas em cálculos de função de onda correlacionada.
3. Resolução do Problema do CO: O trabalho fornece benchmarks de todo-elétron escalar-relativístico de alta qualidade para a adsorção de CO em Cu(111), apoiando a validade do RPA@PBE para a química de superfícies de metais de transição.

Os autores reconhecem limitações, observando que a formulação atual depende do número de condição (uma métrica única) e otimiza apenas expoentes primitivos. Eles sugerem que trabalhos futuros poderiam estender a otimização para coeficientes de contração e refinar a regularização para lidar com múltiplos modos mal condicionados. No entanto, eles postulam que a estratégia MCAO oferece uma abordagem fundamental para a futura otimização de conjuntos de bases para estado sólido.