Stochastisches Tunneln

Stochastisches Tunneln (STUN) ist eine Methode zur globalen Optimierung, in der die zu minimierende Funktion mit der Monte-Carlo-Methode abgetastet wird.

Prinzip

Schematische eindimensionale Testfunktion (Schwarz) und effektives STUN-Potential (Rot und Blau), wobei das mit Pfeilen markierte Minimum das beste bisher gefundene Minimum darstellt. Alle Potentialtöpfe, die über dem besten gefundenen Minimum liegen, sind unterdrückt<ref>K. Hamacher and W. Wenzel: *The Scaling Behaviour of Stochastic Minimization Algorithms in a Perfect Funnel Landscape*. Phys. Rev. E 59(1):938-941, 1999</ref>. Wenn der dynamische Prozess dem Potentialtopf um das aktuelle Minimum entkommen kann, wird er nicht von anderen, höheren lokalen Minima aufgefangen. Potentialtöpfe mit niedrigeren Minima sind verstärkt, was den dynamischen Prozess beschleunigt.

Optimierungsalgorithmen, die auf der Monte-Carlo-Methode beruhen, tasten die untersuchte Funktion ab, indem sie zufällig von der aktuellen Lösung zu einer anderen Lösung springen, wobei die Funktionswerte sich um <math>\Delta E</math> unterscheiden. Als Wahrscheinlichkeit für solch einen Versuchssprung wird meist das als <math> \min\left(1;\exp\left(-\beta\cdot\Delta E\right)\right) </math> definierte Metropolis-Kriterium gewählt, wobei der Parameter <math>\beta</math> passend gewählt wird.

Das Prinzip von STUN ist, die langsame Dynamik von ungünstig geformten Energiefunktionen, die zum Beispiel in Spingläsern angetroffen werden<ref>W. Wenzel and K. Hamacher: A Stochastic tunneling approach for global minimization. Phys. Rev. Lett. 82(15):3003-3007, 1999</ref>, zu umgehen, indem solche Barrieren durchtunnelt werden. Dieses Ziel wird durch die Monte-Carlo-Abtastung<ref>Metropolis, M. ; Rosenbluth, A. ; Rosenbluth, M. ; Teller, A. ; Teller, E.: Equation of state calculations by fast computing machines. J. Chem. Phys 21 (1953), S. 1087–1092</ref> der transformierten Funktion erreicht, die dieser langsamen Dynamik nicht unterliegt. Die Transformation ist in der "Standardform" definiert durch <math>f_{STUN}:=1-\exp\left( -\gamma\cdot\left( f(x)-f_o\right) \right)</math>, wobei <math>f_o</math> der bislang niedrigste gefundene Funktionswert ist. Diese Transformation bewahrt die geometrischen Orte der Minima. Der Effekt solch einer Transformation ist in der Abbildung dargestellt.

Eine adaptive Variante des Verfahrens gestattet es, die Parameter des Verfahrens selber "online" zu schätzen und dadurch die Effizienz des Algorithmus zu verbessern<ref>K. Hamacher: Adaptation in Stochastic Tunneling Global Optimization of Complex Potential Energy Landscapes. Europhys. Lett. 74(6):944, 2006</ref>.

Andere Ansätze

Einzelnachweise