Monte-Carlo-Methode: Was ist damit gemeint?

Verfahren, das Lösungen durch wiederholte stochastische Versuche mit Zufallszahlen ermittelt. In der Werkstoffwissenschaft gehören dazu Gleichgewichtsberechnungen mit dem Metropolis-Algorithmus und die kinetische Monte-Carlo-Methode zur Verfolgung der Diffusion.

Japanische Bezeichnung: モンテカルロ法

„Monte-Carlo-Methode“ kommt in 2 Übungsfragen aus 1 Thema vor, darunter Computergestützte Materialwissenschaft.

Fachgebiet: Computergestützte Materialwissenschaft

Selbst ausprobieren: Fragen zu Monte-Carlo-Methode

Bei der Monte-Carlo-Methode nach Metropolis erhöht eine Probeverschiebung eines Atoms die Energie des Systems um 0,10 eV. Welcher Wert kommt bei einer Temperatur von 1000 K der Wahrscheinlichkeit am nächsten, diesen Schritt anzunehmen? k = 8,617×10⁻⁵ eV/K.

  1. Aca. 0,69
  2. B0 (energieerhöhende Schritte werden nicht angenommen)
  3. Cca. 0,31
Lösung und Erläuterung anzeigen

Lösung: C. ca. 0,31

Beim Metropolis-Verfahren werden Schritte, die die Energie senken, immer angenommen, solche, die sie erhöhen, mit der Wahrscheinlichkeit exp(−ΔE/kT). Mit kT = 8,617×10⁻⁵ × 1000 ≈ 0,0862 eV ist ΔE/kT ≈ 1,16 und die Wahrscheinlichkeit exp(−1,16) ≈ 0,31. Praktisch wird eine Zufallszahl zwischen 0 und 1 erzeugt und der Schritt angenommen, wenn sie kleiner als diese Wahrscheinlichkeit ist. Da auch energieerhöhende Schritte mit einer gewissen Wahrscheinlichkeit angenommen werden, besucht das System Zustände gemäß der Boltzmann-Verteilung, und im Mittel ergeben sich die Gleichgewichtseigenschaften bei endlicher Temperatur. Würden energieerhöhende Schritte nie angenommen, bliebe das System in einem lokalen Energieminimum gefangen, was einer Optimierung bei 0 K entspräche. Ca. 0,69 verwechselt die Ablehnungswahrscheinlichkeit 1 − 0,31 mit der Annahmewahrscheinlichkeit.

Welche Aussage über die kinetische Monte-Carlo-Methode (kMC) trifft zu?

  1. AEreignissen wie Atomsprüngen werden Raten ν exp(−E_a/kT) zugeordnet, Ereignisse mit ratenproportionaler Wahrscheinlichkeit gewählt und die Zeit fortgeschrieben, sodass Diffusion und Ausscheidung über lange Zeiten verfolgbar sind
  2. BSie integriert die Bewegungsgleichungen der Atome in Femtosekundenschritten, verfolgt die thermischen Schwingungen genau und behandelt bei gleichem Aufwand längere Diffusionszeiten als die Molekulardynamik
  3. CSie kennt keine Zeit, nimmt energieerhöhende Versuche mit der Wahrscheinlichkeit exp(−ΔE/kT) an und dient nur zur Bestimmung mittlerer Eigenschaften im Gleichgewicht
Lösung und Erläuterung anzeigen

Lösung: A. Ereignissen wie Atomsprüngen werden Raten ν exp(−E_a/kT) zugeordnet, Ereignisse mit ratenproportionaler Wahrscheinlichkeit gewählt und die Zeit fortgeschrieben, sodass Diffusion und Ausscheidung über lange Zeiten verfolgbar sind

Bei der Diffusion im Festkörper führen Atome sehr viele Schwingungen aus, bevor sie selten in eine benachbarte Leerstelle springen. Die Molekulardynamik verfolgt alle diese Schwingungen, sodass Vorgänge über mehr als Mikrosekunden schwer zugänglich sind. Die kMC lässt die Schwingungen weg: Sie stellt eine Liste von Ereignissen wie Leerstellensprüngen mit ihren Aktivierungsenergien auf, wählt das nächste Ereignis mit einer zur Rate ν exp(−E_a/kT) proportionalen Wahrscheinlichkeit und schreitet in der Zeit gemäß dem Kehrwert der Summe aller Raten voran. Als Aktivierungsenergien werden häufig DFT-Werte verwendet. So lassen sich das Wachstum von Ausscheidungen beim Auslagern oder die Entwicklung von Strahlenschäden in realer Zeit verfolgen. Die Bewegungsgleichungen zu integrieren, kennzeichnet die Molekulardynamik, und ohne Zeitbegriff Gleichgewichtsmittelwerte zu bestimmen, das Metropolis-Verfahren.

Fragen, in denen „Monte-Carlo-Methode“ vorkommt (2)

  1. Bei der Monte-Carlo-Methode nach Metropolis erhöht eine Probeverschiebung eines Atoms die Energie des Systems um 0,10 eV. Welcher Wert kommt bei einer Temperatur von 1000 K der Wahrscheinlichkeit am nächsten, diesen Schritt anzunehmen? k = 8,617×10⁻⁵ eV/K.
  2. Welche Aussage über die kinetische Monte-Carlo-Methode (kMC) trifft zu?

Verwandte Begriffe

Leerstelle Molekulardynamik Dichtefunktionaltheorie CALPHAD-Methode Finite-Elemente-Methode (FEM) Interatomares Potential Materialinformatik Phasenfeldmethode

Zurück zum Glossar