2 \documentclass[landscape,semhelv]{seminar}
5 \usepackage[german]{babel}
6 \usepackage[latin1]{inputenc}
7 \usepackage[T1]{fontenc}
11 \usepackage{calc} % Simple computations with LaTeX variables
12 \usepackage[hang]{caption2} % Improved captions
13 \usepackage{fancybox} % To have several backgrounds
15 \usepackage{fancyhdr} % Headers and footers definitions
16 \usepackage{fancyvrb} % Fancy verbatim environments
17 \usepackage{pstricks} % PSTricks with the standard color package
20 \graphicspath{{../img/}}
23 \usepackage{semlayer} % Seminar overlays
24 \usepackage{slidesec} % Seminar sections and list of slides
26 \input{seminar.bug} % Official bugs corrections
27 \input{seminar.bg2} % Unofficial bugs corrections
33 \extraslideheight{10in}
39 \def\slideleftmargin{5.1cm}
40 \def\slidetopmargin{-0.6cm}
47 \begin{minipage}{2.0cm}
49 \includegraphics[height=1.6cm]{uni-logo.eps}
52 \begin{minipage}{8.0cm}
55 Molekulardynamische Simulation\\
60 \begin{minipage}{2.3cm}
62 \includegraphics[height=1.5cm]{Lehrstuhl-Logo.eps}
66 \begin{picture}(350,40)
71 \begin{picture}(350,40)
85 \begin{picture}(350,10)
88 {\color{blue} Molekulardynamische Simulation} von
89 {\color{green} Silizium}
91 \begin{picture}(350,10)
94 \item {\color{blue} Grundlagen} der {\color{blue}MD-Simulation}
95 \item MD in der Praxis - die Numerik im Detail
96 \item {\color{green} Potential} f"ur {\color{green} kovalent gebundene Materialien}
97 \item Zusammengefasst: Die {\em moldyn} Bibliothek
103 Atomistische Simulationen - Unterschied MC / MD
105 \begin{picture}(350,10)
108 Methoden zur Beschreibung eines Vielteilchen-Systems\\
110 \begin{picture}(350,10)
114 \item stochastische Simulation
115 \item Abrastern des Phasenraums durch Metropolis Algorithmus
116 \item Systemgr"o"se: Mittelwert aus allen Monte Carlo Schritten
118 \begin{picture}(350,10)
122 \item deterministische Simulation
123 \item N"achster Phasenraumpunkt bestimmt durch Bewegungsgleichung
124 \item Systemgr"o"se: Mittelwert aus allen Zeitschritten
126 \begin{picture}(350,10)
128 Ergodenhypothese: Gleichheit der zwei Mittelwerte
133 Prinzip der MD-Simulation
136 \item System von $N$ Teilchen (Molek"ulen)
137 \item zeitliche Entwicklung von Orten und Geschwindigkeiten
138 $\{{\bf q}_i,{\bf p}_i\}$
139 \item System Hamilton-Funktion $\mathcal{H}({\bf q},{\bf p})$
140 \item Hamilton'sche Bewegungsgleichungen:\\
142 \dot{p}_i = \frac{\partial \mathcal{H}}{\partial q_i},
144 \dot{q}_i = \frac{\partial \mathcal{H}}{\partial p_i}
146 Propagationsvorschrift im $6N$-dimensionalen Phasenraum
148 $\Rightarrow$ mikroskopische Beschreibung des Systems\\
149 $\Rightarrow$ observable Gr"o"sen durch zeitlichen Mittelwert\\
150 \begin{picture}(350,10)
153 Notwendige Bestandteile der MD-Simulation
156 \item Methode zum Integrieren der Bewegungsgleichungen (Integrator)
157 \item Modell zur Wechselwirkung (analytische Potentialfunktion)
158 \item Zusatz zur \dq Kontrolle\dq{} des gew"unschten Ensembles
164 Integration der Bewegungsgleichungen
167 \item Keine analytische L"osung f"ur $N>3$ $\Rightarrow$ numerische Integration
168 \item $3N$ DGLs zweiter Ordnung {\color{blue}oder} $6N$ DGLs erster Ordnung\\
170 m_i \ddot{{\bf r}}_i = {\bf f}_i \qquad
171 \textrm{{\color{blue}oder}} \qquad
172 m_i \dot{{\bf r}}_i = {\bf p}_i, \quad \dot{{\bf p}}_i = {\bf f}_i
175 (${\bf f}_i = - \nabla_{{\bf r}_i} \mathcal{V}$, kartesische Koordinaten)
177 \item Prinzip der Finite-Differenzen-Methode\\
179 \Gamma(t) \rightarrow \Gamma(t+\delta t)
182 (mit Anfangsbedingungen ${\bf r}(0)$, ${\bf p}(0)$)
185 ${\bf v}(t+\delta t) = {\bf v}(t) + \delta t {\bf a}(t)$\\
186 ${\bf r}(t+\delta t) = {\bf r}(t) + \delta t {\bf v}(t)$\\
189 Anforderungen an den Integrator
192 \item fehlerfreie Reproduktion der \dq echten\dq{} Trajektorie
193 \item Erhaltung der Energie, Reversibel in der Zeit
194 \item schnell \& nur eine Kraftberechnug pro Zeitschritt $\delta t$
200 {\em Predictor-Corrector} Algorithmus
203 \item Vorhersage der Orte, Geschwindigkeiten, Beschleunigungen etc ...
205 {\bf r}^p(t + \delta t) &=& {\bf r}(t) + \delta t {\bf v}(t) +
206 \frac{1}{2} \delta t^2 {\bf a}(t) +
207 \frac{1}{6} \delta t^3 {\bf b}(t) + \ldots
209 {\bf v}^p(t + \delta t) &=& {\bf v}(t) + \delta t {\bf a}(t) +
210 \frac{1}{2} \delta t^2 {\bf b}(t) + \ldots
212 {\bf a}^p(t + \delta t) &=&{\bf a}(t) + \delta t {\bf b}(t) + \ldots
214 {\bf b}^p(t + \delta t) &=&{\bf b}(t) + \ldots
217 \item Brechnung der tats"achlichen Kraft/Beschleunigung ${\bf a}^c$
218 f"ur die vorhergesagten Orte ${\bf r}^p$ \\
219 $\Rightarrow$ Korrekturfaktor:
220 $\Delta {\bf a}(t + \delta t) =
221 {\bf a}^c(t + \delta t) - {\bf a}^p(t + \delta t)$
224 {\bf r}^c(t + \delta t) &=& {\bf r}^p(t + \delta t) +
225 c_0 \Delta {\bf a}(t + \delta t) \nonumber \\
226 {\bf v}^c(t + \delta t) &=& {\bf v}^p(t + \delta t) +
227 c_1 \Delta {\bf a}(t + \delta t) \nonumber \\
228 {\bf a}^c(t + \delta t) &=& {\bf a}^p(t + \delta t) +
229 c_2 \Delta {\bf a}(t + \delta t) \nonumber \\
230 {\bf b}^c(t + \delta t) &=& {\bf b}^p(t + \delta t) +
231 c_3 \Delta {\bf a}(t + \delta t) \nonumber
233 \item Optional: Iteration des Korrekturschrittes
237 The numerical integration of ordinary differential equations of various orders.
240 Numerical initial value problems in ordinary differential equations.
249 Aus formaler L"osung der Liouville-Gleichung f"ur Ensemble Zeitentwicklung:
251 {\bf r}(t+\delta t) &=& {\bf r}(t) + \delta t {\bf v}(t) +
252 \frac{1}{2} \delta t^2 {\bf a}(t) \nonumber \\
253 {\bf v}(t+\delta t) &=& {\bf v}(t) + \frac{1}{2} \delta t (
254 {\bf a}(t) + {\bf a}(t+\delta t)) \nonumber
258 \item Berechnung der neuen Ortskoordinaten ${\bf r}(t+\delta t)$
259 \item Erste Berechnung der Geschwindigkeiten
261 {\bf v}(t+\delta t/2) = {\bf v}(t) + \frac{1}{2} \delta t {\bf a}(t)
263 \item Berechnung der Kr"afte f"ur die Orte ${\bf r}(t+\delta t)$
264 $\Rightarrow {\bf a}(t+\delta t)$
265 \item Update der Geschwindigkeiten
267 {\bf v}(t+\delta t) = {\bf v}(t+\delta t/2) +
268 \frac{1}{2} \delta t {\bf a}(t+\delta t)
273 \item entspricht {\em GEAR-3} mit Ortskorrekturfaktor $c_0=0$
274 \item einfach, schnell, wenig Speicheraufwand $(9N)$
275 \item verh"altnism"a"sig pr"azise
281 Modell zur Wechselwirkung - Das Potential
283 Klassisches Potential:
285 {\mathcal V} = \sum_i {\mathcal V}_1({\bf r}_i) +
286 \sum_{i,j} {\mathcal V}_2({\bf r}_i,{\bf r}_j) +
287 \sum_{i,j,k} {\mathcal V}_3({\bf r}_i,{\bf r}_j,{\bf r}_k) +
290 \begin{minipage}{8.3cm}
292 \item ${\mathcal V}_1$: Eink"orperpotential (Gravitation, elektrisches Feld)
293 \item ${\mathcal V}_2$: Paarpotential
294 (nur abh"angig vom Abstand $r_{ij}$)
295 \item ${\mathcal V}_3$: Dreik"orperpotential
296 (oft ${\mathcal V}_3(r_{ij},r_{ik},\theta_{ijk})$)
297 \item oft umgebungsabh"angiger Term in ${\mathcal V}_2$ eingebaut
298 \item Terme h"oherer Ordnung vermutlich klein (Biologie, Chemie)
299 \item nur Paarpotential\\
300 $\Rightarrow$ hcp im Grundzustand\\
301 $\Rightarrow$ ungen"ugend f"ur kovalent gebundene Materialien
304 \begin{minipage}{4cm}
305 \includegraphics[width=4.3cm]{tersoff_angle.eps}
307 Cut-Off Radius $r_c$:
309 \item Atome $i$, $j$ mit $r_{ij} > r_c$ wechselwirken nicht
310 \item Korrektur f"ur Potentiale die erst im Unendlichen verschwinden
312 \begin{picture}(350,5)
314 Kraft: ${\bf F}_i = - \nabla_{{\bf r}_i} \mathcal{V}$
319 Wahl/Kontrolle des Ensembles
328 kanonisches Ensemble (NVT)
334 isothermales isobares Ensemble (NpT)
340 Die Simulationszelle \& Randbedingungen
346 Trick: Nachbarlisten \& Zell-Methode
352 Thermodynamische Gr"o"sen
376 Albe Reparametrisierung