Multidimensional derivative-free optimization. A case study on minimization of Hartree-Fock-Roothaan energy functionals

Resumo Técnico: Otimização Multidimensional Livre de Derivadas para a Minimização da Energia de Hartree-Fock-Roothaan

Definição do Problema
Este estudo aborda o desafio de minimizar os funcionais de energia de Hartree-Fock-Roothaan (HFR) para sistemas atômicos empregando orbitais do tipo Slater não inteiros (n*-STOs). Nestes cálculos, a função objetivo depende de parâmetros orbitais não lineares ($\zeta$) e números quânticos principais de ordem não inteira ($n^*$). Uma dificuldade crítica surge porque as derivadas analíticas do funcional de energia em relação a esses parâmetros não estão prontamente disponíveis. Isso se deve à não analiticidade das expansões de séries de potências para funções com expoentes reais e à complexidade computacional da avaliação de integrais atômicos sobre n*-STOs, que frequentemente requerem métodos numéricos especializados (por exemplo, funções hiper-radiais) em vez de expressões em forma fechada. Consequentemente, técnicas de otimização tradicionais baseadas em gradiente são inaplicáveis, necessitando do uso de métodos de otimização livre de derivadas (DFO).

Metodologia
Os autores implementaram e compararam quatro algoritmos distintos de DFO sob condições numéricas idênticas utilizando a linguagem de programação Mathematica. As implementações foram explícitas, evitando comandos de otimização integrados para garantir a transparência algorítmica. Os métodos avaliados foram:

1. Direção Conjugada de Powell (Powell CD): Um método de busca direta utilizando um conjunto de direções conjugadas atualizadas via um mecanismo de rotação, empregando buscas de linha (Seção Áurea e método de Brent) dentro de restrições de caixa.
2. Símplex de Nelder-Mead (NM Simplex): Um método de busca direta que evolui um símplex através de operações de reflexão, expansão, contração e encolhimento, apresentando um mecanismo de reinicialização automática para evitar a convergência prematura.
3. Busca de Padrão Baseada em Coordenadas (PS-C): Um método de busca direta que realiza sondagens ao longo de direções de coordenadas com tamanhos de passo ajustados adaptativamente, utilizando uma estratégia oportunista e etapas de aceleração.
4. Função de Base Radial Baseada em Modelo (MB-RBF): Uma abordagem de modelo substituto construindo um interpolante usando funções de base radial (RBF) dentro de um arcabouço de região de confiança. Este método utilizou regularização diagonal adaptativa para garantir estabilidade numérica sem caudas polinomiais.

Os algoritmos foram primeiramente testados no Powell Singular Function (PSF) em 4 e 8 dimensões, um caso de teste conhecido por sua singularidade de Hessiana no mínimo global. Subsequentemente, foram aplicados à minimização da energia HFR para configurações atômicas de camada fechada, especificamente para as séries isoeletrônicas de He-like e Be-like. Os vetores de otimização incluíram números quânticos principais não inteiros e expoentes orbitais, com dimensões de problema variando de 2 a 8 parâmetros não lineares.

Resultos Principais

• Benchmarking (PSF): No teste PSF de 4 dimensões, o Powell CD e o NM Simplex convergiram para valores residuais pequenos com custo computacional moderado. O PS-C exibiu convergência mais lenta e precisão reduzida. O MB-RBF alcançou precisão intermediária, mas com um custo computacional significativamente maior devido à construção do modelo substituto. No caso de 8 dimensões, o custo computacional (avaliações de função e tempo de CPU) para o PS-C e MB-RBF aumentou substancialmente.
• Cálculos Atômicos (Base Mínima): Para íons He-like e Be-like usando conjuntos de bases mínimas, o algoritmo NM Simplex consistentemente exigiu menos avaliações de função e menos tempo de CPU do que o Powell CD e o PS-C. Embora o MB-RBF tenha produzido energias aceitáveis, seu overhead computacional foi considerado muito alto para cálculos atômicos repetidos na prática.
• Base Estendida e Altas Dimensões: À medida que a dimensionalidade aumentava (por exemplo, átomo de Be com 8 parâmetros orbitais), o método PS-C tornou-se a principal fonte de custo computacional, levando os autores a excluí-lo de comparações subsequentes de alta dimensão. O Powell CD permaneceu eficaz para problemas de menor dimensão, mas mostrou crescente sensibilidade às configurações de tolerância conforme o espaço de parâmetros crescia. O NM Simplex manteve a estabilidade e forneceu energias de estado fundamental convergidas, embora tenha exigido um pouco mais de tempo de CPU que o Powell CD em testes específicos de conjunto de base estendida de alta precisão.
• Precisão: O algoritmo NM Simplex otimizou com sucesso a configuração de 8 parâmetros para o átomo de Be, produzindo uma energia de estado fundamental ($E = -14.573023164$ u.a.) que concordou com os benchmarks de Hartree-Fock analíticos reotimizados relatados na literatura.

Alegações e Significância
O artigo afirma apresentar a primeira investigação sistemática de métodos de otimização livre de derivadas especificamente para a minimização de funcionais de energia de Hartree-Fock-Roothaan envolvendo orbitais do tipo Slater não inteiros. O estudo destaca que:

• O algoritmo Nelder-Mead Simplex entrega o desempenho mais consistente e confiável entre os métodos DFO testados para esta classe específica de problemas, equilibrando velocidade de convergência e custo computacional.
• Embora o método de Direção Conjugada de Powell seja eficaz para problemas de menor dimensão, seu desempenho degrada-se com o aumento do tamanho do espaço de parâmetros e da sensibilidade aos critérios de parada.
• Abordagens baseadas em modelos (MB-RF) e variantes de busca de padrão (PS-C) incorrem em custos computacionais proibitivos para estes cálculos atômicos específicos, limitando sua utilidade prática neste contexto.
• Os algoritmos desenvolvidos são transferíveis para outros ambientes de programação e são adequados para sistemas atômicos maiores e potencialmente para cálculos moleculares envolvendo funções de base do tipo exponencial.
• O estudo fornece um benchmark para futuros desenvolvimentos em estratégias de otimização variacional para orbitais do tipo Slater não inteiros, uma lacuna anteriormente não abordada na literatura.

O trabalho demonstra que o tratamento eficaz de n*-STOs é possível através da minimização direta do funcional de energia usando estratégias DFO robustas, oferecendo um caminho para se aproximar do limite de Hartree-Fock com menos funções de base através da otimização de parâmetros não lineares.