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 Intgrieren 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) +
291 \item ${\mathcal V}_1$: Eink"orperpotential (Gravitation, elektrisches Feld)
292 \item ${\mathcal V}_2$: Paarpotential
293 (nur abh"angig vom Abstand ${\bf r}_{ij}$)
294 \item ${\mathcal V}_3$: Dreik"orperpotential
300 Wahl/Kontrolle des Ensembles
306 kanonisches Ensemble (NVT)
312 isothermales isobares Ensemble (NpT)
318 Die Simulationszelle \& Randbedingungen
324 Trick: Nachbarlisten \& Zell-Methode
330 Thermodynamische Gr"o"sen
354 Albe Reparametrisierung