X-Git-Url: https://hackdaworld.org/gitweb/?p=lectures%2Flatex.git;a=blobdiff_plain;f=posic%2Ftalks%2Fmd_simulation_von_silizium.tex;fp=posic%2Ftalks%2Fmd_simulation_von_silizium.tex;h=61eadb65e98302a4d4316b38b7ab5a3256c07578;hp=3282577aae68f5b969957cbec4adbdbc2c209c53;hb=ab8d9f996a4bbe6444f7d6a76eec05209f10076f;hpb=a8210d608f38e87d1ec09cf134207dd4af596cda
diff --git a/posic/talks/md_simulation_von_silizium.tex b/posic/talks/md_simulation_von_silizium.tex
index 3282577..61eadb6 100644
--- a/posic/talks/md_simulation_von_silizium.tex
+++ b/posic/talks/md_simulation_von_silizium.tex
@@ -139,7 +139,7 @@ Ergodenhypothese: Gleichheit der zwei Mittelwerte
\item System Hamilton-Funktion $\mathcal{H}({\bf q},{\bf p})$
\item Hamilton'sche Bewegungsgleichungen:\\
\[
- \dot{p}_i = \frac{\partial \mathcal{H}}{\partial q_i},
+ \dot{p}_i = - \frac{\partial \mathcal{H}}{\partial q_i},
\qquad
\dot{q}_i = \frac{\partial \mathcal{H}}{\partial p_i}
\]\\
@@ -318,9 +318,11 @@ Kraft: ${\bf F}_i = - \nabla_{{\bf r}_i} \mathcal{V}$
{\large\bf
Wahl/Kontrolle des Ensembles
}\\
+\begin{picture}(350,10)
+\end{picture}
Erinnerung:
\begin{itemize}
- \item Stichproben aus Zust"anden im Phasenraum
+ \item Stichproben aus Zust"anden im Phasenraum, $_{ens} = _t$
\item Bewegungsgleichung als Propagationsvorschrift $\Rightarrow$ Gesamtenergie erhalten
\item Au"serdem konstant: $N$ und $V$
\end{itemize}
@@ -328,52 +330,139 @@ $\Rightarrow$ Simulation eines NVE-Ensembles
\[
\rho_{ens}=\delta(H(t)-E)
\]
-F"ur ander Ensembles:
+\begin{picture}(350,10)
+\end{picture}
+F"ur andere Ensembles:
\begin{itemize}
- \item Anpassung der Bewegungsgleichungen
- \item Tricks
+ \item Anpassung der Bewegungsgleichungen f"ur eine Sequenz von Konfigurationen
+ im gew"unschten Ensemble
+\end{itemize}
+\begin{center}
+{\color{red} oder}
+\end{center}
+\begin{itemize}
+ \item Tricks zur Kontrolle von $T$ und $p$
+ $\Rightarrow$ $NVE \rightarrow NVT,NpT$\\
+ Anmerkung: $T$ und $p$ fluktuieren,
+ Mittelwerte entsprechen den gew"unschten Werten
\end{itemize}
\end{slide}
\begin{slide}
{\large\bf
kanonisches Ensemble (NVT)
-}
+}\\
+\begin{picture}(350,10)
+\end{picture}
+Trick: {\em temperature scaling}
+\begin{itemize}
+ \item forcieren der gew"unschten Temperatur in jedem Schritt
+ \item $E_{kin} = 3/2 Nk_BT$
+ \item eigentlich {\em velocity scaling}
+ \item Berendsen Thermostat:
+ \[
+ \lambda = \sqrt{1+\frac{\delta t}{\tau_T}\Big(\frac{T_{ref}}{T}-1\Big)}
+ \]
+ \begin{center}
+ $\tau_T>100\times\delta t \Rightarrow$ reale thermische Fluktuationen\\
+ {\scriptsize Berendsen et al. J. Chem. Phys. 81 (1984) 3684.}
+ \end{center}
+\end{itemize}
+Andersen:
+\begin{itemize}
+ \item Zuf"alliges "Andern der Geschwindigkeit eines Atoms entsprechend
+ der Temperatur
+ \item Physikalische Interpretation: Kopplung an W"armebad
+ \item {\color{green} n"utzlich zum Berechnen thermodynamischer Gr"o"sen}
+ \item {\color{red} nicht geeignet zur Beschreibung atomistischer Prozesse}\\
+ (unphysikalische St"orung der Bewegung des einzelnen Atoms)
+\end{itemize}
\end{slide}
\begin{slide}
{\large\bf
isothermales isobares Ensemble (NpT)
-}
+}\\
+\begin{picture}(350,10)
+\end{picture}
+Trick: {\em pressure scaling}
+\begin{itemize}
+ \item analog zum {\em temperature scaling}
+ \item $p = - \frac{\partial \mathcal{V}}{V}$ (Alternative sp"ater)
+ \item eigentlich {\em volume scaling}
+ \item Berendsen Barostat:
+ \[
+ \mu = \Big[1-\frac{\delta t}{\tau_p}\beta (p_0-p)\Big]^{1/3}
+ \]
+\end{itemize}
+\begin{picture}(350,10)
+\end{picture}
+Andersen:
+\begin{itemize}
+ \item modifizierte Bewegungsgleichung
+ (neue Variable $Q$, ${\bf \rho}_i = {\bf r}_i/V^{1/3}$)
+ \[
+ \mathcal{L}(\rho^N,\dot{\rho}^N,Q,\dot{Q})
+ =\frac{1}{2}mQ^{2/3}\sum_i \dot{\rho}_i^2 -
+ \sum_{i