@@ -905,13 +905,6 @@ \subsubsection{Improving the script}
905905\textit {initial.lmp } file \emph {without } the \textit {create\_ atoms }
906906commands in the \textit {System definition } part.
907907
908- \begin {figure }
909- \centering
910- \includegraphics [width=0.55\linewidth ]{LJ-cylinder}
911- \caption {Visualization of the improved system after minimization}
912- \label {fig:improved-min }
913- \end {figure }
914-
915908We want to create the the atoms of types 1 and 2 in two separate
916909regions. Thus we need to add two \textit {region } commands and then
917910re-add \textit {create\_ atoms } commands, but using the new regions
@@ -986,6 +979,13 @@ \subsubsection{Improving the script}
986979
987980\paragraph {Restarting from a saved configuration }
988981
982+ \begin {figure }
983+ \centering
984+ \includegraphics [width=0.55\linewidth ]{LJ-cylinder}
985+ \caption {Visualization of the improved system after minimization}
986+ \label {fig:improved-min }
987+ \end {figure }
988+
989989To continue a simulation from the saved configuration, open the
990990\textit {improved.md.lmp } which was downloaded during the tutorial setup.
991991This contains the \textit {Initialization } part from \textit {initial.lmp }
@@ -1087,6 +1087,16 @@ \subsubsection{Improving the script}
10871087first compute is referenced by the second. And we reference the second
10881088in a \textit {thermo\_ style custom } command.
10891089
1090+ \begin {figure }
1091+ \centering
1092+ \includegraphics [width=\linewidth ]{LJ-evolution}
1093+ \caption {Evolution of the system during mixing. The three snapshots show
1094+ respectively the system at $ t=0 $ (left panel), $ t=75 $ (middle panel), and $ t=1500 $
1095+ (right panel).}
1096+ \label {fig:evolution-population }
1097+ \end {figure }
1098+
1099+ Finally, let us complete the script by adding the following lines into
10901100Add the following lines into \textit {improved.md.lmp }. Note there is no
10911101need for a \textit {Settings } section, because the settings are taken
10921102from the \textit {.data } file.
@@ -1109,16 +1119,6 @@ \subsubsection{Improving the script}
11091119\textit {sumcoor12 }, were added to the list of information printed during
11101120the simulation.
11111121
1112- \begin {figure }
1113- \centering
1114- \includegraphics [width=\linewidth ]{LJ-evolution}
1115- \caption {Evolution of the system during mixing. The three snapshots show
1116- respectively the system at $ t=0 $ (left panel), $ t=75 $ (middle panel), and $ t=1500 $
1117- (right panel).}
1118- \label {fig:evolution-population }
1119- \end {figure }
1120-
1121- Finally, let us complete the script by adding the following lines into
11221122\textit {improved.md.lmp }:
11231123{\normalsize
11241124\begin {verbatim }
@@ -1187,6 +1187,20 @@ \subsubsection{Improving the script}
11871187\label {fig:coords-viz }
11881188\end {figure }
11891189
1190+ \paragraph {Experiments }
1191+
1192+ Here are some suggestions for further experiments with this system that
1193+ may lead to additional insights into how different systems are set up
1194+ and how different features work:
1195+ \begin {itemize }
1196+ \item Use Nos\' e--Hoover thermostat (\textit {fix nvt }) instead of Langevin
1197+ (\textit {fix nve } + \textit {fix langevin })
1198+ \item Skip the \textit {minimize } step for both, Nos\' e--Hoover and Langevin
1199+ \item Apply the thermostat to only one type of atoms and observe the
1200+ temperature for each type separately
1201+ \item append an NVE run (i.e.~without any thermostat) and observe the energies
1202+ \end {itemize }
1203+
11901204It is also possible to color the atoms in the \textit {Slide Show }
11911205according to their coordination number. For that purpose, replace the
11921206\textit {dump } and \textit {dump\_ modify } commands with the following
@@ -1207,18 +1221,6 @@ \subsubsection{Improving the script}
12071221one type 2 neighbor in yellow, and atoms of type 1 with two type 2
12081222neighbors in red. Figure \ref {fig:coords-viz } shows an example.
12091223
1210- \paragraph {Experiments }
1211-
1212- Here are some suggestions for further experiments with this system:
1213- \begin {itemize }
1214- \item Use Nos\' e--Hoover thermostat (\textit {fix nvt }) instead of Langevin
1215- (\textit {fix nve } + \textit {fix langevin })
1216- \item Skip the \textit {minimize } step for both, Nos\' e--Hoover and Langevin
1217- \item Apply the thermostat to only one type of atoms and observe the
1218- temperature for each type separately
1219- \item append an NVE run (i.e.~without any thermostat) and observe the energies
1220- \end {itemize }
1221-
12221224\subsection {Tutorial 2: Pulling on a carbon nanotube }
12231225\label {carbon-nanotube-label }
12241226
0 commit comments