@@ -651,7 +651,7 @@ \subsubsection{My first input}
651651of the potential well that determines the interaction strength, and
652652$ \sigma _{ij}$ is the distance at which the potential energy equals zero.
653653The indices $ i$ and $ j$ refer to pairs of atoms {\color {blue} with the
654- respective atom type value }. The fourth line, \lmpcmd {pair\_ coeff 1 1
654+ corresponding atom types }. The fourth line, \lmpcmd {pair\_ coeff 1 1
655655 1.0 1.0}, specifies the Lennard-Jones coefficients for interactions
656656between pairs of atoms {\color {blue} that both have} atom type 1: the
657657energy parameter $ \epsilon _{11} = 1.0 $ and the distance parameter
@@ -660,7 +660,7 @@ \subsubsection{My first input}
660660$ \epsilon _{22} = 0.5 $ , and $ \sigma _{22} = 3.0 $ .
661661
662662\begin {note }
663- By default, LAMMPS calculates the {\color {blue} mixed atom type } force
663+ By default, LAMMPS calculates the {\color {blue} mixed} force
664664 field coefficients for different atom types using geometric averages:
665665 $ \epsilon _{ij} = \sqrt {\epsilon _{ii} \epsilon _{jj}}$ ,
666666 $ \sigma _{ij} = \sqrt {\sigma _{ii} \sigma _{jj}}$ . In the present case,
@@ -703,11 +703,11 @@ \subsubsection{My first input}
703703
704704\begin {note } {\color {blue}The `thermodynamic information' printed by
705705 LAMMPS using \lmpcmd {thermo\_ style custom} keywords refers to
706- instantaneous values of those thermodynamic properties at the
707- specified steps , not cumulative averages. But it is also possible
708- to reference a wide variety of custom data from compute styles, fix
709- styles, and variables which can be used for on-the-fly analysis
710- including cumulative and sliding window averages.}
706+ instantaneous values of the specified thermodynamic properties
707+ at each output step , not cumulative averages. However, LAMMPS also
708+ allows to reference a wide variety of custom data from compute styles, fix
709+ styles, and variables. These can be used for on-the-fly analysis,
710+ including cumulative and sliding- window averages.}
711711\end {note }
712712
713713You can now run LAMMPS {\color {blue}(see subsection \ref {running-lammps-label }
@@ -746,13 +746,9 @@ \subsubsection{My first input}
746746selected algorithm using the computed forces, aiming to reduce the
747747potential energy. By default, LAMMPS uses the conjugate gradient (CG)
748748algorithm~\cite {hestenes1952methods }. The simulation will stop as soon
749- as {\color {blue}one of the four minimizer criteria are met and LAMMPS
750- will output which stopping criterion ended the minimization and
751- and some system properties at the initial and the final step.}
752- % SG: I don't think that its true, its rather the algorithm
753- % will stop when one of the four criteria is met. Axel, what do you think?
754- % I propose to replace by "when specific convergence criteria are met"
755- % AK: yes this is correct.
749+ as {\color {blue}one of the four minimizer criteria is met. LAMMPS
750+ will then report which stopping criterion was satisfied, along with
751+ selected system properties at both the initial and final steps.}
756752Note that, except for trivial systems, minimization algorithms will find a
757753local minimum rather than the global minimum.
758754
@@ -1177,7 +1173,7 @@ \subsubsection{Improving the script}
11771173 name of the compute style: global data (no suffix), local data
11781174 (/local suffix), per-atom data (/atom suffix), per-chunk data
11791175 (/chunk suffix), per-gridpoint data (/grid suffix). In the example
1180- above the \lmpcmd {compute coord/atom} produces per-atom data, which
1176+ above, the \lmpcmd {compute coord/atom} produces per-atom data, which
11811177 is used as input for \lmpcmd {compute reduce} which returns global
11821178 data. For global data three kinds of data exists: scalars (single
11831179 values), vectors (one-dimensional arrays), or arrays
@@ -4108,6 +4104,7 @@ \subsubsection{Method 1: Free sampling}
41084104% into several variables
41094105% CA: I think we can remove the & (to avoid confusion) and warn the reader about the
41104106% continuation of the formula.
4107+ % SG: its tricky because in principle all commands from here can be copy-pasted...
41114108\ begin{lstlisting}
41124109variable U atom ${U0}*atan((x+${x0})/${dlt})&
41134110 -${U0}*atan((x-${x0})/${dlt})
@@ -4126,14 +4123,12 @@ \subsubsection{Method 1: Free sampling}
41264123in the NVT ensemble, maintaining a constant number of
41274124atoms $ N$ , constant volume $ V$ , and a temperature $ T$ that
41284125fluctuates around a target value.
4129- % SG: may be discuss the choice of constant "500" -> chosen for a relatiely weak coupling with thermostat (add a box?)
4130- % CA: I can propose something that can be put in a yellow box. Feel free to uncomment the lines below and modify it if you want to.
4131- % {\color{blue}You can find proposed in the LAMMPS documentation
4132- % values of 100x and 1000x the value of the timestep for the damping
4133- % constants of the thermostats and barostats. Here, a smaller
4134- % value is used in order to have the temperature of the system
4135- % relaxed to the target value more often.}
4136- % PS: I am not defining the damping constants because you had already done it in one of the firsts tutorials.
4126+
4127+ \begin {note }
4128+ {\color {blue}LAMMPS documentation suggests using damping constants for thermostats
4129+ that are approximately 100 times the timestep value. In this case, a value of 500
4130+ is used, resulting in a relatively weak coupling to the thermostat.}
4131+ \end {note }
41374132
41384133\begin {figure }
41394134\centering
0 commit comments