Molekulardynamik: Was ist damit gemeint?

Verfahren, das die Newtonschen Bewegungsgleichungen vieler Atome mit Kräften aus interatomaren Potentialen numerisch integriert und so die Bewegung der Atome im zeitlichen Verlauf berechnet.

Japanische Bezeichnung: 分子動力学法

„Molekulardynamik“ kommt in 8 Übungsfragen aus 1 Thema vor, darunter Computergestützte Materialwissenschaft.

Fachgebiet: Computergestützte Materialwissenschaft

Selbst ausprobieren: Fragen zu Molekulardynamik

In der Molekulardynamik sollen die Bewegungsgleichungen mit einem Zeitschritt von 2 fs integriert werden, um die Bewegung der Atome über 1 ns zu berechnen. Wie viele Zeitschritte sind nötig?

  1. A5×10²
  2. B5×10⁸
  3. C5×10⁵
Lösung und Erläuterung anzeigen

Lösung: C. 5×10⁵

Mit 1 ns = 10⁻⁹ s und 2 fs = 2×10⁻¹⁵ s ist die Zahl der Schritte 10⁻⁹/(2×10⁻¹⁵) = 5×10⁵. Der Zeitschritt muss so klein sein, dass die Periode der Atomschwingungen (etwa 100 fs) in einige Dutzend Schritte unterteilt wird; in Systemen mit leichtem Wasserstoff wählt man 1 fs oder weniger. Ist der Schritt zu groß, bleibt die Energie wegen der Fehler der numerischen Integration (z. B. mit dem Verlet-Algorithmus) nicht erhalten, und die Rechnung bricht zusammen. Der Wert 5×10² verwechselt fs mit 10⁻¹² s (ps), der Wert 5×10⁸ verwechselt ns mit 10⁻⁶ s (µs). Da in jedem Schritt die Kräfte auf alle Atome berechnet werden, liegt die mit Molekulardynamik erreichbare Zeit als Richtwert bei Nanosekunden bis Mikrosekunden.

Welche Aussage über die in der Molekulardynamik verwendeten interatomaren Potentiale trifft zu?

  1. ADas Lennard-Jones-Potential ist ein Paarpotential, das nur vom Abstand zweier Atome abhängt, und beschreibt auch die metallische Bindung und die Richtungsabhängigkeit kovalenter Bindungen gut
  2. BDas Potential der Embedded-Atom-Methode (EAM) fügt eine Einbettungsenergie hinzu, die von der Elektronendichte der umgebenden Atome abhängt, und beschreibt so den Vielteilchencharakter der metallischen Bindung
  3. CFür kovalente Kristalle wie Silizium ist ein Paarpotential, das nicht von den Bindungswinkeln abhängt, am besten geeignet
Lösung und Erläuterung anzeigen

Lösung: B. Das Potential der Embedded-Atom-Methode (EAM) fügt eine Einbettungsenergie hinzu, die von der Elektronendichte der umgebenden Atome abhängt, und beschreibt so den Vielteilchencharakter der metallischen Bindung

Das Lennard-Jones-Potential ist ein Paarpotential, das nur vom Abstand zweier Atome abhängt; es beschreibt schwache, ungerichtete Bindungen wie bei Edelgasen gut, nicht aber die Eigenschaften metallischer oder kovalenter Bindungen. In Metallen spielt ein Vielteilcheneffekt eine große Rolle: Die Bindungsstärke eines Atoms hängt von der Zahl seiner Nachbarn ab. Die EAM berücksichtigt dies als Funktion der von den umgebenden Atomen erzeugten Elektronendichte. Mit reinen Paarpotentialen gilt für spannungsfreie kfz- und krz-Kristalle, in denen jedes Atom ein Symmetriezentrum ist, für die elastischen Konstanten stets C₁₂ = C₄₄ (Cauchy-Relation), was viele Metalle nicht erfüllen. Bei kovalenten Bindungen sind die Bindungswinkel wichtig, daher werden Potentiale mit Winkeltermen wie das Stillinger-Weber- oder das Tersoff-Potential verwendet. Verbreitet ist auch, ein interatomares Potential durch maschinelles Lernen an DFT-Ergebnisse anzupassen (Machine-Learning-Potentiale).

„Molekulardynamik“ in den Erläuterungen

Fragen, in denen „Molekulardynamik“ vorkommt (8)

  1. In der Molekulardynamik sollen die Bewegungsgleichungen mit einem Zeitschritt von 2 fs integriert werden, um die Bewegung der Atome über 1 ns zu berechnen. Wie viele Zeitschritte sind nötig?
  2. Welche Aussage über die in der Molekulardynamik verwendeten interatomaren Potentiale trifft zu?
  3. In einer Molekulardynamik-Rechnung eines einatomigen Systems betrug die mittlere kinetische Energie pro Atom 0,0517 eV. Welcher Wert kommt der nach dem Gleichverteilungssatz bestimmten Temperatur des Systems am nächsten? Die Boltzmann-Konstante sei k = 8,617×10⁻⁵ eV/K.
  4. Welche Zuordnung eines Rechenverfahrens der Werkstoffsimulation zu den räumlichen und zeitlichen Größenordnungen, die es üblicherweise behandeln kann, ist am zutreffendsten?
  5. Welche Aussage über die Grundideen der Dichtefunktionaltheorie (DFT) trifft zu?
  6. Welche Tendenz zeigen die Ergebnisse von Standard-DFT-Rechnungen mit LDA- oder GGA-Austausch-Korrelations-Funktionalen?
  7. Welche Aussage über die kinetische Monte-Carlo-Methode (kMC) trifft zu?
  8. Welche Aussage über die Phasenfeldmethode trifft zu?

Verwandte Begriffe

Dichtefunktionaltheorie Interatomares Potential Phasenfeldmethode Kubisch-flächenzentriertes Gitter (kfz) Kubisch-raumzentriertes Gitter (krz) Finite-Elemente-Methode (FEM) Phonon Monte-Carlo-Methode Leerstelle CALPHAD-Methode

Zurück zum Glossar