@@ -3229,12 +3229,13 @@ \subsubsection{Prepare and relax}
32293229
32303230run 5000
32313231
3232- write_data relax.data
3232+ write_data relax.data nofix
32333233\end {lstlisting }
3234+ {\color {blue}The \lmpcmd {write\_ data} command is used with the \lmpcmd {nofix} keyword to
3235+ print a data file without extra sections from the \lmpcmd {reaxff/species} command.}
32343236Run the \flecmd {relax.lmp} file using LAMMPS. As seen from \flecmd {relax.species},
32353237only one species is detected, called \lmpcmd {O384Si192}, representing the entire system.
32363238
3237-
32383239\begin {figure }
32393240\centering
32403241\includegraphics [width=\linewidth ]{SIO-charge}\\ [-2ex]
@@ -3338,7 +3339,7 @@ \subsubsection{Deform the structure}
33383339
33393340run 25000
33403341
3341- write_data deform.data
3342+ write_data deform.data nofix
33423343\end {lstlisting }
33433344{\color {blue}The \lmpcmd {fix deform} command applies a continuous deformation
33443345by elongating the simulation box along the x-axis at a constant engineering
@@ -4028,7 +4029,7 @@ \subsubsection{Method 1: Free sampling}
40284029
40294030To begin this tutorial, select \guicmd {Start Tutorial 7} from the
40304031\guicmd {Tutorials} menu of \lammpsgui {} and follow the instructions.
4031- The editor should display the following content corresponding to \flecmd {free-energy .lmp}:
4032+ The editor should display the following content corresponding to \flecmd {free-sampling .lmp}:
40324033\ begin{lstlisting}
40334034variable sigma equal 3.405
40344035variable epsilon equal 0.238
@@ -4064,7 +4065,7 @@ \subsubsection{Method 1: Free sampling}
40644065\paragraph {System creation and settings }
40654066
40664067Let us define the simulation box and randomly add atoms by addying the
4067- following lines to \flecmd {free-energy .lmp}:
4068+ following lines to \flecmd {free-sampling .lmp}:
40684069\ begin{lstlisting}
40694070region myreg block -50 50 -15 15 -50 50
40704071create_box 1 myreg
@@ -4112,7 +4113,7 @@ \subsubsection{Method 1: Free sampling}
41124113
41134114We impose the force $ F(x)$ to the atoms in the simulation
41144115using the \lmpcmd {fix addforce} command. Add the following
4115- lines to \flecmd {free-energy .lmp}:
4116+ lines to \flecmd {free-sampling .lmp}:
41164117% SG: apparently variable can't be cutted in 2 lines, either warn the reader, or find
41174118% a way to write these commands on the entire width of the page, or divide the variables
41184119% into several variables
@@ -4164,7 +4165,7 @@ \subsubsection{Method 1: Free sampling}
41644165thermo_style custom step temp etotal v_n_center
41654166thermo 10000
41664167
4167- dump viz all image 50000 myimage-*.ppm type type &
4168+ dump viz all image 5000 myimage-*.ppm type type &
41684169 shiny 0.1 box yes 0.01 view 180 90 zoom 6 &
41694170 size 1600 500 fsaa yes
41704171dump_modify viz backcolor white acolor 1 cyan &
@@ -4187,7 +4188,7 @@ \subsubsection{Method 1: Free sampling}
41874188
41884189Once the system is equilibrated, we will record the density profile of
41894190the atoms along the $ x$ -axis using the \lmpcmd {ave/chunk} command.
4190- Add the following line to \flecmd {free-energy .lmp}:
4191+ Add the following line to \flecmd {free-sampling .lmp}:
41914192\ begin{lstlisting}
41924193reset_timestep 0
41934194
@@ -4230,7 +4231,7 @@ \subsubsection{Method 1: Free sampling}
42304231and compare it
42314232with the imposed potential $ U$ from Eq.~\eqref {eq:U } (Fig.~\ref {fig:US-density }\, b).
42324233The reference density, $ \rho _\text {bulk} = 0.0009 ~\text {\AA {}}^{-3}$ ,
4233- was estimated by measuring the density of the reservoir from the raw density
4234+ was estimated by measuring the density of the reservoir from the density
42344235profiles. The agreement between the MD results and the imposed energy profile
42354236is excellent, despite some noise in the central part, where fewer data points
42364237are available due to the repulsive potential.
0 commit comments