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

Technische Samenvatting: Multidimensionale Afgeleide-vrije Optimalisatie voor Hartree-Fock-Roothaan Energie-minimalisatie

Probleemstelling
Deze studie behandelt de uitdaging van het minimaliseren van Hartree-Fock-Roothaan (HFR) energiefuncties voor atomaire systemen die gebruikmaken van niet-gehele Slater-type orbitalen (n*-STOs). In deze berekeningen is de objectieve functie afhankelijk van niet-lineaire orbitale parameters ($\zeta$) en hoofdkwantumgetallen van niet-gehele orde ($n^*$). Een kritieke moeilijkheid ontstaat omdat analytische afgeleiden van de energiefunctie met betrekking tot deze parameters niet gemakkelijk beschikbaar zijn. Dit komt door de niet-analytische aard van machtsreeksontwikkelingen voor functies met reële exponenten en de computationele complexiteit van het evalueren van atomaire integralen over n*-STOs, die vaak gespecialiseerde numerieke methoden (bijv. hyperradiale functies) vereisen in plaats van gesloten vorm-uitdrukkingen. Bijgevolg zijn traditionele gradiënt-gebaseerde optimalisatietechnieken niet toepasbaar, wat het gebruik van afgeleide-vrije optimalisatiemethoden (DFO) noodzakelijk maakt.

Methodologie
De auteurs hebben vier verschillende DFO-algoritmen geïmplementeerd en vergeleken onder identieke numerieke condities met behulp van de programmeertaal Mathematica. De implementaties waren expliciet, waarbij ingebouwde optimalisatiecommando's werden vermeden om algoritmische transparantie te waarborgen. De geëvalueerde methoden waren:

1. Powell's Conjugate Direction (Powell CD): Een directe zoekmethode die gebruikmaakt van een set geconjugeerde richtingen die via een rotatiemechanisme worden bijgewerkt, waarbij gebruik wordt gemaakt van lijnzoekopdrachten (Golden Section en Brent's methode) binnen box-restricties.
2. Nelder-Mead Simplex (NM Simplex): Een directe zoekmethode die een simplex evolueert via reflectie-, expansie-, contractie- en krimpoperaties, voorzien van een automatisch herstartmechanisme om voortijdige convergentie te voorkomen.
3. Coordinate-based Pattern Search (PS-C): Een directe zoekmethode die pollt langs coördinatierichtingen met adaptief aangepaste stapgroottes, gebruikmakend van een opportunistische strategie en acceleratiestappen.
4. Model-Based Radial Basis Function (MB-RBF): Een surrogaatmodel-benadering die een interpolant construeert met behulp van radiale basisfuncties (RBF) binnen een trust-region kader. Deze methode maakte gebruik van adaptieve diagonale regularisatie om numerieke stabiliteit te garanderen zonder polynomiale staarten.

De algoritmen werden eerst getest op de Powell Singular Function (PSF) in 4 en 8 dimensies, een testgeval dat bekend staat om zijn Hessiaanse singulariteit bij het globale minimum. Vervolgens werden ze toegepast op HFR energie-minimalisatie voor gesloten-schil atomaire configuraties, specifiek de He-achtige en Be-achtige isoelektrische series. De optimalisatievectoren omvatten niet-gehele hoofdkwantumgetallen en orbitale exponenten, waarbij de dimensies van het probleem varieerden van 2 tot 8 niet-lineaire parameters.

Belangrijkste Resultaten

• Benchmarking (PSF): In de 4-dimensionale PSF-test convergeerden Powell CD en NM Simplex naar kleine residuwaarden met een matige computationele kostenpost. PS-C vertoonde tragere convergentie en een lagere nauwkeurigheid. MB-RBF bereikte een intermediaire nauwkeurigheid, maar tegen een aanzienlijk hogere computationele kosten door de constructie van het surrogaatmodel. In het 8-dimensionale geval nam de computationele kosten (functie-evaluaties en CPU-tijd) voor PS-C en MB-RBF substantieel toe.
• Atomaire Berekeningen (Minimale Basis): Voor He-achtige en Be-achtige ionen met minimale basissets vereiste het NM Simplex-algoritme consequent minder functie-evaluaties en minder CPU-tijd dan Powell CD en PS-C. Hoewel MB-RBF aanvaardbare energieën produceerde, werd de computationele overhead als te hoog beschouwd voor praktische herhaalde atomaire berekeningen.
• Uitgebreide Basis en Hoge Dimensies: Naarmate de dimensionaliteit toenam (bijv. Be-atoom met 8 orbitale parameters), werd de PS-C methode de primaire bron van computationele kosten, wat leidde ertoe dat de auteurs deze uitsloten van verdere vergelijkingen in hoge dimensies. Powell CD bleef effectief voor lager-dimensionale problemen, maar vertoonde een toenemende gevoeligheid voor tolerantie-instellingen naarmate de parameterruimte groeide. NM Simplex behield stabiliteit en leverde geconvergeerde grondtoestandsenergieën op, hoewel het in specifieke hoog-precieze uitgebreide basisset-tests iets meer CPU-tijd vereiste dan Powell CD.
• Nauwkeurigheid: Het NM Simplex-algoritme optimaliseerde succesvol de 8-parameter configuratie voor het Be-atoom, wat een grondtoestandsenergie opleverde ($E = -14.573023164$ a.u.) die overeenkwam met de geheroptimaliseerde analytische Hartree-Fock benchmarks gerapporteerd in de literatuur.

Claims en Betekenis
Het artikel beweert de eerste systematische onderzoek te presenteren naar afgeleide-vrije optimalisatiemethoden specificaal voor de minimalisatie van Hartree-Fock-Roothaan energiefuncties met niet-gehele Slater-type orbitalen. De studie benadrukt dat:

• Het Nelder-Mead Simplex algoritme de meest consistente en betrouwbare prestaties levert onder de geteste DFO-methoden voor deze specifieke klasse problemen, waarbij een balans wordt gevonden tussen convergentiesnelheid en computationele kosten.
• Hoewel Powell's CD methode effectief is voor lager-dimensionale problemen, verslechtert de prestatie bij een toenemende omvang van de parameterruimte en de gevoeligheid voor stopcriteria.
• Model-gebaseerde benaderingen (MB-RBF) en varianten van pattern search (PS-C) brengen prohibitieve computationele kosten met zich mee, wat hun praktische bruikbaarheid in deze context beperkt.
• De ontwikkelde algoritmen zijn overdraagbaar naar andere programmeeromgevingen en zijn geschikt voor grotere atomaire systemen en potentieel moleculaire berekeningen met exponentiële type basisfuncties.
• De studie biedt een benchmark voor toekomstige ontwikkelingen in variationele optimalisatiestrategieën voor niet-gehele Slater-type orbitalen, een gat dat voorheen onbehandeld bleef in de literatuur.

Het werk demonstreert dat een effectieve behandeling van n*-STOs mogelijk is door directe minimalisatie van de energiefunctie met behulp van robuuste DFO-strategieën, wat een pad biedt om de Hartree-Fock limiet te benaderen met minder basisfuncties door niet-lineaire parameters te optimaliseren.