@@ -1072,8 +1072,7 @@ \subsubsection{Improving the script}
10721072coordinated atoms, and the second calculates the average over all type 1
10731073atoms. Add the following lines to \flecmd {improved.md.lmp}:
10741074\ begin{lstlisting}
1075- compute coor12 grp_t1 coord/atom cutoff 2 &
1076- group grp_t2
1075+ compute coor12 grp_t1 coord/atom cutoff 2 group grp_t2
10771076compute sumcoor12 grp_t1 reduce ave c_coor12
10781077\end {lstlisting }
10791078The \lmpcmd {compute reduce ave} command is used to average the per-atom
@@ -1128,7 +1127,8 @@ \subsubsection{Improving the script}
11281127\includegraphics [width=\linewidth ]{LJ-mixing}
11291128\caption {a)~Evolution of the numbers $ N_\text {1, in}$ and $ N_\text {2, in}$ of atoms
11301129of types 1 and 2, respectively, within the \lmpcmd {cyl\_ in} region as functions
1131- of time $ t$ . b)~Evolution of the coordination number $ C_{1-2}$ between atoms of types 1 and 2.}
1130+ of time $ t$ . b)~Evolution of the coordination number $ C_{1-2}$ (compute \lmpcmd {sumcoor12})
1131+ between atoms of types 1 and 2.}
11321132\label {fig:mixing }
11331133\end {figure }
11341134
@@ -1143,14 +1143,14 @@ \subsubsection{Improving the script}
11431143force applied to the atoms is equal to zero.
11441144
11451145\begin {note }
1146- The steps to ensure no initial linear momentum and no net random
1146+ The steps to ensure no initial linear momentum and no net total
11471147force are important to prevent the system from starting to drift or move as a
11481148whole. For a bulk system with periodic boundary conditions, it is
1149- expected to remain in place. Thus , when computing the temperature from the
1149+ expected to remain in place. Accordingly , when computing the temperature from the
11501150kinetic energy, we use $ 3 N-3 $ degrees of freedom since there is no
11511151global translation. In a drifting system, some of the kinetic energy is
11521152due to the drift, which means the system itself cools down. In
1153- extreme cases, it can freeze and drift very fast . This phenomenon is
1153+ extreme cases, the system can freeze while its center of mass drifts very quickly . This phenomenon is
11541154sometimes referred to as the `` flying ice cube syndrome'' ~\cite {wong2016good }.
11551155\end {note }
11561156
@@ -1175,7 +1175,7 @@ \subsubsection{Improving the script}
11751175\caption {Snapshot of the binary mixture simulated
11761176 during \hyperref [lennard-jones-label]{Tutorial 1} with atoms of type 1
11771177 colored according to their computed $ 1 -2 $ coordination
1178- number \lmpcmd {c \_ coor12 }, ranging from turquoise, \lmpcmd {c\_ coor12 = 0},
1178+ number from the compute \lmpcmd {coor12 }, ranging from turquoise,\lmpcmd {c\_ coor12 = 0},
11791179 to yellow, \lmpcmd {c\_ coor12 = 1}, and red, \lmpcmd {c\_ coor12 = 2}.}
11801180\label {fig:coords-viz }
11811181\end {figure }
@@ -1186,15 +1186,16 @@ \subsubsection{Improving the script}
11861186may lead to additional insights into how different systems are configured
11871187and how various features function:
11881188\begin {itemize }
1189- \item Use Nos\' e--Hoover thermostat (\lmpcmd {fix nvt}) instead of Langevin
1189+ \item Use the Nos\' e--Hoover thermostat (\lmpcmd {fix nvt}) instead of the Langevin thermostat
11901190 (\lmpcmd {fix nve} + \lmpcmd {fix langevin}).
1191- \item Omit the energy minimization step for both, Nos\' e--Hoover and Langevin.
1191+ \item Omit the energy minimization step before starting the MD simulation using either
1192+ the Nos\' e--Hoover or the Langevin thermostat.
11921193\item Apply the thermostat to only one type of atoms and observe the
11931194 temperature for each type separately.
11941195\item Append an NVE run (i.e.~without any thermostat) and observe the energy levels.
11951196\end {itemize }
11961197
1197- Another useful experiment to try is to color the atoms in the \guicmd {Slide Show}
1198+ Another useful experiment is coloring the atoms in the \guicmd {Slide Show}
11981199according to an observable, such as their respective coordination
11991200numbers. To do this, replace the
12001201\lmpcmd {dump} and \lmpcmd {dump\_ modify} commands with the following lines:
@@ -1205,7 +1206,7 @@ \subsubsection{Improving the script}
12051206dump_modify viz adiam 1 1 adiam 2 3 backcolor white &
12061207 amap -1 2 ca 0.0 4 min royalblue 0 turquoise 1 yellow max red
12071208\end {lstlisting }
1208- Run LAMMPS again. Now, the atoms of type 1 are colored based on the value
1209+ Run LAMMPS again. Atoms of type 1 are now colored based on the value
12091210of \lmpcmd {c\_ coor12}, which is mapped continuously from turquoise to yellow
12101211and red for atoms with the highest coordination (Fig.~\ref {fig:coords-viz }).
12111212In the definition of the variable \lmpcmd {v\_ coor12}, atoms of type 2 are
0 commit comments