spelling, formatting issues

Author: jrgissing

lammps-tutorials.tex

@@ -341,7 +341,7 @@ \subsection{Background knowledge}
   with Applications to Soft Matter} by Jean-Pierre Hansen and Ian Ranald
 McDonald~\cite{hansen2013theory}.  For more resources, the SklogWiki
 platform provies a wide range of information information on statistical mechanics
-and molecular simulations~\cite{sklogwiki_main_page}. 
+and molecular simulations~\cite{sklogwiki_main_page}.
 
 \subsection{Software/system requirements}
 
@@ -498,7 +498,7 @@ \subsubsection{My first input}
 a popular choice for simulations that explore general statistical mechanical
 principles, as it emphasizes relative differences between parameters rather than
 representing any specific material.  The second line, \lmpcmd{dimension 3}, specifies that the simulation is conducted
-in 3D space, as opposed to 2D, where atoms are confined to move only in the 
+in 3D space, as opposed to 2D, where atoms are confined to move only in the
 xy-plane.  The third line, \lmpcmd{atom\_style atomic}, designates
 the atomic style for representing simple, individual particles.
 In this style, each particle is treated as a point with a mass, making it the
@@ -681,7 +681,7 @@ \subsubsection{My first input}
 algorithm~\cite{hestenes1952methods}.  The simulation will stop as soon
 as the minimizer algorithm cannot find a way to lower the potential
 energy. % SG: I don't think that its true, its rather the algorithm
-% will stop when one of the four criteria is met. Axel, what do you think? 
+% will stop when one of the four criteria is met. Axel, what do you think?
 % I propose to replace by "when specific convergence criteria are met"
 Note that, except for trivial systems, minimization algorithms will find a
 local minimum rather than the global minimum.
@@ -1135,7 +1135,7 @@ \subsubsection{Improving the script}
 \caption{a)~Evolution of the numbers $N_\text{1, in}$ and $N_\text{2, in}$ of atoms
 of types 1 and 2, respectively, within the \lmpcmd{cyl\_in} region as functions
 of time $t$.  b)~Evolution of the coordination number $C_{1-2}$ (compute \lmpcmd{sumcoor12})
-between atoms of types 1 and 2.} 
+between atoms of types 1 and 2.}
 \label{fig:mixing}
 \end{figure}
 
@@ -1699,7 +1699,7 @@ \subsubsection{Breakable bonds}
 bonds break, the energy relaxes abruptly, as can be seen near $t=32~\text{ps}$ in Fig.~\ref{fig:CNT-breakable-energy-stress}\,a.
 Using a similar script as previously,
 i.e.,~\href{\filepath tutorial2/unbreakable-yaml-reader.py}{\dwlcmd{unbreakable-yaml-reader.py}},
-import the data into Python and generate the stress-strain curve (Fig.~\ref{fig:CNT-breakable-energy-stress}\,b).  The 
+import the data into Python and generate the stress-strain curve (Fig.~\ref{fig:CNT-breakable-energy-stress}\,b).  The
 stress-strain curve reveals a linear (elastic) regime where $F_\text{cnt} \propto \Delta L_\text{cnt}$
 for $\Delta L_\text{cnt} < 5\,\%$, and a non-linear (plastic) regime
 for $5\,\% < \Delta L_\text{cnt} < 25\,\%$.
@@ -2111,7 +2111,7 @@ \subsubsection{Stretching the PEG molecule}
 dump mydmp all local 100 pull.dat index c_dphi c_prop
 \end{lstlisting}
 By contrast with the radius of gyration (compute \lmpcmd{rgyr}), the dihedral angle  % SG: why is c_prop not printed?
-$\phi$ (compute \lmpcmd{dphi}) is returned as a vector by the \lmpcmd{compute dihedral/local} 
+$\phi$ (compute \lmpcmd{dphi}) is returned as a vector by the \lmpcmd{compute dihedral/local}
 command and must be written to a file using the \lmpcmd{dump local} command.
 
 Finally, let us simulate 15 picoseconds without any external force:
@@ -2750,7 +2750,7 @@ \subsubsection{Imposed shearing}
 the average force on each wall is given by \lmpcmd{f\_mysf1[1]} and \lmpcmd{f\_mysf2[1]}
 and is approximately $2.7\,\mathrm{kcal/mol/\AA}$ in magnitude.  Using a surface area
 for the walls of $A = 6 \cdot 10^{-18}\,\text{m}^2$, one obtains an estimate for
-the shear viscosity for the confined fluid of $\eta = 3.1\,\text{mPa.s}$ using Eq.~\eqref{eq:eta}. 
+the shear viscosity for the confined fluid of $\eta = 3.1\,\text{mPa.s}$ using Eq.~\eqref{eq:eta}.
 
 \begin{note}
 The viscosity calculated at such a high shear rate may differ from the expected
@@ -4021,7 +4021,7 @@ \subsubsection{Creating the system}
 The \textit{class2} bond, angle, dihedral, and improper styles are used as
 well, see the documentation for a description of their respective potentials.
 The \lmpcmd{mix sixthpower} imposes the following mixing rule for the calculation
-of the cross coefficients: 
+of the cross coefficients:
 \begin{eqnarray}
 \nonumber
 \sigma_{ij} & = & 2^{-1/6} (\sigma^6_i+\sigma_j^6)^{1/6}, ~ \text{and} \\
@@ -4055,7 +4055,7 @@ \subsubsection{Creating the system}
 by manually shrinking the simulation box at a constant rate.  The dimension parallel
 to the CNT axis is maintained fixed because the CNT is periodic in that direction.
 Add the following commands to \flecmd{mixing.lmp}:
-% SG: I removed the loop local, unless its important? But if it is, we have to 
+% SG: I removed the loop local, unless its important? But if it is, we have to
 % explain what it does and why it was chosen here.
 \begin{lstlisting}
 velocity all create 530 9845 dist gaussian rot yes
@@ -4071,7 +4071,7 @@ \subsubsection{Creating the system}
 run 9000
 \end{lstlisting}
 The \lmpcmd{fix halt} command is used to stop the box shrinkage once the
-target density is reached. 
+target density is reached.
 
 For the next stage of the simulation, we will use \lmpcmd{dump image} to
 output images every 200 steps:
@@ -4096,7 +4096,7 @@ \subsubsection{Creating the system}
 unfix myhal
 reset_timestep 0
 
-group CNT molecule 1 2 3
+group CNT molecule 1
 fix myrec CNT recenter NULL 0 0 units box shift all
 
 run 10000
@@ -4234,7 +4234,7 @@ \subsubsection{Simulating the reaction}
 \centering
 \includegraphics[width=0.7\linewidth]{REACT-final.png}
 \caption{Final configuration for \hyperref[bond-react-label]{Tutorial 8}.
-The atoms from the formed polymer named \lmpcmd{c1}, \lmpcmd{c2}, and 
+The atoms from the formed polymer named \lmpcmd{c1}, \lmpcmd{c2}, and
 \lmpcmd{c3} are colored in pink.}
 \label{fig:REACT-final}
 \end{figure}
@@ -4244,7 +4244,7 @@ \subsubsection{Simulating the reaction}
 to track which atoms are being stabilized and which atoms are undergoing
 dynamics with the system-wide time integrator (here, \lmpcmd{fix nvt}).
 When reaction stabilization is employed, there should not be a time integrator acting on
-the group \lmpcmd{all}.  Instead, the group of atoms not currently
+the group \mbox{\lmpcmd{all}.}  Instead, the group of atoms not currently
 undergoing stabilization is named by appending `\_REACT' to the user-provided prefix.
 \end{note}
 
@@ -4273,10 +4273,10 @@ \subsubsection{Simulating the reaction}
 \caption{a) Evolution of the system temperature, $T$,
 as a function of the time, $t$, during the polymerization step of
 \hyperref[bond-react-label]{Tutorial 8}.
-b) Evolution of the three reaction counts, corresponding respecively to
-the polymerization of two styrene monomers (Rnx~1), the  addition of a styrene
-monomer to the end of a growing polymer chain (Rnx~2), and to the linking
-of two polymer chains (Rnx~3).}
+b) Evolution of the three reaction counts, corresponding respectively to
+the polymerization of two styrene monomers (Rxn~1), the  addition of a styrene
+monomer to the end of a growing polymer chain (Rxn~2), and to the linking
+of two polymer chains (Rxn~3).}
 \label{fig:evolution-reacting}
 \end{figure}