@@ -1664,7 +1664,7 @@ \subsubsection{Breakable bonds}
16641664
16651665After equilibration, let us set the velocity of the edges equal to
16661666$ 75 ~\text {m/s}$ (or $ 0.75 ~\text {\AA {}/ps}$ ) and run for a longer duration than
1667- previously. Add the following lines into \flecmd {breakable.lmp}:
1667+ previously. Add the following lines into \flecmd {breakable.lmp}:
16681668\ begin{lstlisting}
16691669velocity cnt_top set 0.75 0 0
16701670velocity cnt_bot set -0.75 0 0
@@ -1942,7 +1942,7 @@ \subsubsection{Solvating the PEG in water}
19421942\end {figure }
19431943
19441944Open the file named \flecmd {merge.lmp} that was downloaded
1945- alongside \flecmd {water.lmp} during the tutorial setup. It only contain one line:
1945+ alongside \flecmd {water.lmp} during the tutorial setup. It only contain one line:
19461946\ begin{lstlisting}
19471947read_restart water.restart
19481948\end {lstlisting }
@@ -3492,7 +3492,7 @@ \subsubsection{Adding water}
34923492with time. The \lmpcmd {compute\_ modify} command with the \lmpcmd {dynamic yes}
34933493option for water is used to specify that the number of molecules will not be constant.
34943494
3495- Finally, let us use the \textit {fix gcmc } and perform the grand canonical Monte
3495+ Finally, let us use the \lmpcmd {fix gcmc} and perform the grand canonical Monte
34963496Carlo steps. Add the following lines into \flecmd {gcmc.lmp}:
34973497\ begin{lstlisting}
34983498variable tfac equal 5.0/3.0
@@ -3505,12 +3505,10 @@ \subsubsection{Adding water}
35053505freedom. Here, 100 insertion and deletion attemps are made every 100 steps.
35063506
35073507\begin {note }
3508- At a pressure of $ p = 100 \ \text {bar}$ , the chemical potential of water
3509- vapor at $ T = 300 \ \text {K}$ can be calculated using as
3510- $ \mu = \mu ^\circ + RT \ln (\frac {p}{p_0}),$
3511- where $ \mu _0 $ is the standard chemical potential
3512- at $ p^\circ = 1 \, \text {bar}$ , \( R = 8.314\ \text {J/mol·K}\) is
3513- the gas constant, \( T = 300\ \text {K}\) is the temperature.
3508+ At a pressure of $ p = 100 \ \text {bar}$ , the chemical potential of water vapor at $ T = 300 \ \text {K}$
3509+ can be calculated using as $ \mu = \mu _0 + RT \ln (\frac {p}{p_0}),$ where $ \mu _0 $ is the standard
3510+ chemical potential (typically taken at a pressure $ p_0 = 1 \, \text {bar}$ ), \( R = 8.314\ \text {J/mol·K}\)
3511+ is the gas constant, \( T = 300\ \text {K}\) is the temperature.
35143512\end {note }
35153513
35163514Finally, let us print some information and run for 25\, ps:
@@ -3617,7 +3615,7 @@ \subsubsection{Method 1: Free sampling}
36173615to create a Weeks-Chandler-Andersen (WCA) potential, which is a truncated and
36183616purely repulsive LJ potential~\cite {weeks1971role }. It was calculated
36193617as $ 2 ^{1/6} \sigma $ . The potential is also shifted to be
3620- equal to 0 at the cut-off using the \lmpcmd {pair\_ modify}. The system of unit
3618+ equal to 0 at the cut-off using the \lmpcmd {pair\_ modify} command . The system of unit
36213619\lmpcmd {real}, in which energy is in kcal/mol, distance in Ångstrom, or time in
36223620femtosecond, has been chosen for practical reasons: the WHAM algorithm used in
36233621the second part of the tutorial automatically assumes the energy to be in kcal/mol.
@@ -4122,21 +4120,21 @@ \subsubsection{Reaction templates}
41224120The first reaction uses the prefix `M-M' for the pre-reaction template,
41234121post-reaction template, and reaction map file:
41244122\begin {itemize }
4125- \item \href {\filepath tutorial8/M-M_pre.mol}{\textit {M-M$ \_ $ pre.mol }},
4126- \item \href {\filepath tutorial8/M-M_post.mol}{\textit {M-M$ \_ $ post.mol }},
4127- \item \href {\filepath tutorial8/M-M.rxnmap}{\textit {M-M.rxnmap }}.
4123+ \item \href {\filepath tutorial8/M-M_pre.mol}{\dwlcmd {M-M$ \_ $ pre.mol}},
4124+ \item \href {\filepath tutorial8/M-M_post.mol}{\dwlcmd {M-M$ \_ $ post.mol}},
4125+ \item \href {\filepath tutorial8/M-M.rxnmap}{\dwlcmd {M-M.rxnmap}}.
41284126\end {itemize }
41294127The second reaction uses the prefix `M-P',
41304128\begin {itemize }
4131- \item \href {\filepath tutorial8/M-P_pre.lmpmol}{\textit {M-P$ \_ $ pre.mol }},
4132- \item \href {\filepath tutorial8/M-P_post.lmpmol}{\textit {M-P$ \_ $ post.mol }},
4133- \item \href {\filepath tutorial8/M-P.rxnmap}{\textit {M-P.rxnmap }}.
4129+ \item \href {\filepath tutorial8/M-P_pre.lmpmol}{\dwlcmd {M-P$ \_ $ pre.mol}},
4130+ \item \href {\filepath tutorial8/M-P_post.lmpmol}{\dwlcmd {M-P$ \_ $ post.mol}},
4131+ \item \href {\filepath tutorial8/M-P.rxnmap}{\dwlcmd {M-P.rxnmap}}.
41344132\end {itemize }
41354133The third reaction uses the prefix `P-P',
41364134\begin {itemize }
4137- \item \href {\filepath tutorial8/P-P_pre.lmpmol}{\textit {P-P$ \_ $ pre.mol }},
4138- \item \href {\filepath tutorial8/P-P_post.lmpmol}{\textit {P-P$ \_ $ post.mol }},
4139- \item \href {\filepath tutorial8/P-P.rxnmap}{\textit {P-P.rxnmap }}.
4135+ \item \href {\filepath tutorial8/P-P_pre.lmpmol}{\dwlcmd {P-P$ \_ $ pre.mol}},
4136+ \item \href {\filepath tutorial8/P-P_post.lmpmol}{\dwlcmd {P-P$ \_ $ post.mol}},
4137+ \item \href {\filepath tutorial8/P-P.rxnmap}{\dwlcmd {P-P.rxnmap}}.
41404138\end {itemize }
41414139Here, the file names for each reaction use the abbreviation `M' for monomer and `P'
41424140for polymer.
0 commit comments