@@ -1682,10 +1682,11 @@ \subsubsection{Breakable bonds}
16821682
16831683\begin {note }
16841684 Bonds cannot be displayed by the \lmpcmdnote {dump image} when using
1685- the \lmpcmdnote {atom\_ style atomic}, as it contains no bonds. A tip
1686- for displaying bonds with the present system using LAMMPS is provided
1687- at the end of the tutorial. You can also use external tools like VMD
1688- or OVITO (see \hyperref [tip-dynamic-bonds]{tip for tutorial 3}).
1685+ the \lmpcmdnote {atom\_ style atomic}, as it contains no bonds. A
1686+ \hyperref [tip-dynamic-bonds]{tip for displaying bonds} with the
1687+ present system using LAMMPS is provided at the end of the tutorial.
1688+ You can also use external tools like VMD or OVITO (see the
1689+ \hyperref [tip-external-viz]{tip for tutorial 3}).
16891690\end {note }
16901691
16911692\paragraph {Launch the deformation }
@@ -1767,6 +1768,22 @@ \subsubsection{Breakable bonds}
17671768\subsection {Tutorial 3: Polymer in water }
17681769\label {all-atom-label }
17691770
1771+ The goal of this tutorial is to use LAMMPS to solvate a small
1772+ hydrophilic polymer (PEG - polyethylene glycol) in a reservoir of water
1773+ (Fig.~\ref {fig:PEG }). Once the water reservoir is properly equilibrated
1774+ at the desired temperature and pressure, the polymer molecule is added
1775+ and a constant stretching force is applied to both ends of the polymer.
1776+ The evolution of the polymer length is measured as a function of time.
1777+ The GROMOS 54A7 force field~\cite {schmid2011definition } is used for the
1778+ PEG, the SPC/Fw model~\cite {wu2006flexible } is used for the water, and
1779+ the long-range Coulomb interactions are solved using the PPPM
1780+ solver~\cite {luty1996calculating }. This tutorial was inspired by a
1781+ publication by Liese and coworkers, in which molecular dynamics
1782+ simulations are compared with force spectroscopy experiments, see
1783+ Ref.\, ~\citenum {liese2017hydration}.
1784+
1785+ \subsubsection {Preparing the water reservoir }
1786+
17701787\begin {figure }
17711788\centering
17721789\includegraphics [width=0.55\linewidth ]{PEG}
@@ -1776,19 +1793,6 @@ \subsection{Tutorial 3: Polymer in water}
17761793\label {fig:PEG }
17771794\end {figure }
17781795
1779- \noindent The goal of this tutorial is to use LAMMPS to solvate a small hydrophilic
1780- polymer (PEG - polyethylene glycol) in a reservoir of water (Fig.~\ref {fig:PEG }).
1781- Once the water reservoir is properly equilibrated at the desired temperature and
1782- pressure, the polymer molecule is added and a constant stretching force is applied
1783- to both ends of the polymer. The evolution of the polymer length is measured as
1784- a function of time. The GROMOS 54A7 force field~\cite {schmid2011definition } is used
1785- for the PEG, the SPC/Fw model~\cite {wu2006flexible } is used for the water, and the
1786- long-range Coulomb interactions are solved using the PPPM solver~\cite {luty1996calculating }.
1787- This tutorial was inspired by a publication by Liese and coworkers, in which molecular
1788- dynamics simulations are compared with force spectroscopy experiments, see Ref.\, ~\citenum {liese2017hydration}.
1789-
1790- \subsubsection {Preparing the water reservoir }
1791-
17921796In this tutorial, the water reservoir is first prepared in the absence of the polymer.
17931797A rectangular box of water is created and equilibrated at ambient temperature and
17941798pressure. The SPC/Fw water model is used~\cite {wu2006flexible }, which is
@@ -1844,12 +1848,12 @@ \subsubsection{Preparing the water reservoir}
18441848\end {lstlisting }
18451849
18461850\begin {note }
1847- This tutorial uses type labels~\cite {typelabel_paper } to map each of the
1848- numeric atom types with a string (see the \flecmd {parameters.inc} file):
1849- \lmpcmdnote {labelmap atom 1 OE 2 C 3 HC 4 H 5 CPos 6 OAlc 7 OW 8 HW}
1850- Therefore, the oxygen and hydrogen atoms of water (respectively types 7 and 8)
1851- can be referred to as `OW' and `HW', respectively. Similar maps are used for
1852- the bond types, angle types, and dihedral types.
1851+ This tutorial uses type labels~\cite {typelabel_paper } to map each
1852+ numeric atom type to a string (see the \flecmd {parameters.inc} file):
1853+ \lmpcmdnote {labelmap atom 1 OE 2 C 3 HC 4 H 5 CPos 6 OAlc 7 OW 8 HW}
1854+ Therefore, the oxygen and hydrogen atoms of water (respectively types
1855+ 7 and 8) can be referred to as `OW' and `HW', respectively. Similar
1856+ maps are used for the bond types, angle types, and dihedral types.
18531857\end {note }
18541858
18551859Let us create water molecules. To do so, let us import a molecule template called
@@ -2194,13 +2198,14 @@ \subsubsection{Stretching the PEG molecule}
21942198\end {figure }
21952199
21962200\paragraph {Tip: using external visualization tools }
2201+ \label {tip-external-viz }
21972202
2198- Trajectories can be visualized using external tools such as VMD
2199- or OVITO~\cite {humphrey1996vmd , ovito_paper }. To do so, the IDs and
2203+ Trajectories can be visualized using external tools such as VMD or
2204+ OVITO~\cite {humphrey1996vmd , ovito_paper }. To do so, the IDs and
22002205positions of the atoms must be regularly written to a file during the
2201- simulation. This can be accomplished by adding a \lmpcmd {dump}
2202- command to the input file. For instance, create a duplicate of \flecmd {pull.lmp}
2203- and name it
2206+ simulation. This can be accomplished by adding a \lmpcmd {dump} command
2207+ to the input file. For instance, create a duplicate of
2208+ \flecmd {pull.lmp} and name it
22042209% do not wrap this line
22052210\href {\filepath tutorial3/solution/pull-with-tip.lmp}{\dwlcmd {pull-with-tip.lmp}}.
22062211% do not wrap this line
@@ -2212,13 +2217,14 @@ \subsubsection{Stretching the PEG molecule}
22122217named \flecmd {pull.lammpstrj}, which can be opened in OVITO or VMD.
22132218
22142219\begin {note }
2215- Since the trajectory dump file does not contain information about topology and atom
2216- types, it is usually preferred to first write out a data file and import it directly
2217- (in the case of OVITO) or convert it to a PSF file (for VMD). This allows the topology
2218- to be loaded before \emph {adding } the
2219- trajectory file to it. When using \lammpsgui {}, this process can be automated
2220- through the \guicmd {View in OVITO} or \guicmd {View in VMD} options in the \guicmd {Run} menu.
2221- Afterward, only the trajectory dump needs to be added.
2220+ Since the trajectory dump file does not contain information about
2221+ topology and elements, it is usually preferred to first write out a
2222+ data file and import it directly (in the case of OVITO) or convert it
2223+ to a PSF file (for VMD). This allows the topology to be loaded before
2224+ \emph {adding } the trajectory file to it. When using \lammpsgui {},
2225+ this process can be automated through the \guicmd {View in OVITO} or
2226+ \guicmd {View in VMD} options in the \guicmd {Run} menu. Afterwards
2227+ only the trajectory dump needs to be added.
22222228\end {note }
22232229
22242230\subsection {Tutorial 4: Nanosheared electrolyte }
@@ -2235,7 +2241,7 @@ \subsection{Tutorial 4: Nanosheared electrolyte}
22352241\label {fig:NANOSHEAR }
22362242\end {figure }
22372243
2238- \noindent The objective of this tutorial is to simulate an electrolyte
2244+ The objective of this tutorial is to simulate an electrolyte
22392245nanoconfined and sheared between two walls (Fig.~\ref {fig:NANOSHEAR }). The density
22402246and velocity profiles of the fluid in the direction normal to the walls are
22412247extracted to highlight the effect of confining a fluid on its local properties.
@@ -3214,7 +3220,7 @@ \subsection{Tutorial 6: Water adsorption in silica}
32143220\label {fig:GCMC }
32153221\end {figure }
32163222
3217- \noindent The objective of this tutorial is to combine molecular dynamics and
3223+ The objective of this tutorial is to combine molecular dynamics and
32183224grand canonical Monte Carlo simulations to compute the adsorption of water
32193225molecules in cracked silica material (Fig.~\ref {fig:GCMC }). This tutorial
32203226illustrates the use of the grand canonical ensemble in molecular simulation, an
@@ -3224,7 +3230,7 @@ \subsection{Tutorial 6: Water adsorption in silica}
32243230
32253231\subsubsection {Generation of the silica block }
32263232
3227- \noindent To begin this tutorial, select \guicmd {Start Tutorial 6} from the
3233+ To begin this tutorial, select \guicmd {Start Tutorial 6} from the
32283234\guicmd {Tutorials} menu of \lammpsgui {} and follow the instructions.
32293235The editor should display the following content corresponding to \flecmd {generate.lmp}:
32303236\ begin{lstlisting}
@@ -3394,21 +3400,21 @@ \subsubsection{Cracking the silica}
33943400
33953401\subsubsection {Adding water }
33963402
3397- \noindent To add the water molecules to the silica, we will employ
3398- the Monte Carlo method in the grand canonical ensemble (GCMC). In short, the
3399- system is placed into contact with a virtual reservoir of a given chemical potential
3400- $ \mu $ , and multiple attempts to insert water molecules at random positions are
3401- made. Each attempt is either accepted or rejected based on energy considerations.
3402- For further details, please refer to classical textbooks like Ref.~\citenum {frenkel2023understanding}.
3403+ To add the water molecules to the silica, we will employ the Monte Carlo
3404+ method in the grand canonical ensemble (GCMC). In short, the system is
3405+ placed into contact with a virtual reservoir of a given chemical
3406+ potential $ \mu $ , and multiple attempts to insert water molecules at
3407+ random positions are made. Each attempt is either accepted or rejected
3408+ based on energy considerations. For further details, please refer to
3409+ classical textbooks like Ref.~\citenum {frenkel2023understanding}.
34033410
34043411\paragraph {Using hydrid potentials }
34053412
3406- \noindent The first particularly of our system is that it
3407- combine water and silica, which necessitates the use
3408- of two force fields: Vashishta (for $ \text {SiO}_2 $ ), and
3409- TIP4P (for water). Here, the TIP4P/2005 model is employed for the
3410- water~\cite {abascal2005general }. Open the \flecmd {gcmc.lmp} file, which
3411- should contain the following lines:
3413+ The first particularly of our system is that it combine water and
3414+ silica, which necessitates the use of two force fields: Vashishta (for
3415+ $ \text {SiO}_2 $ ), and TIP4P (for water). Here, the TIP4P/2005 model is
3416+ employed for the water~\cite {abascal2005general }. Open the
3417+ \flecmd {gcmc.lmp} file, which should contain the following lines:
34123418\ begin{lstlisting}
34133419units metal
34143420boundary p p p
@@ -3639,21 +3645,24 @@ \subsection{Tutorial 7: Free energy calculation}
36393645\label {fig:US }
36403646\end {figure }
36413647
3642- \noindent The objective of this tutorial is to measure the free energy profile
3643- of particles through a barrier potential using two methods: free sampling and
3644- umbrella sampling~\cite {kastner2011umbrella , allen2017computer , frenkel2023understanding } (Fig.~\ref {fig:US }).
3645- To simplify the process and minimize computation time, the barrier potential will be
3646- imposed on the atoms using an additional force, mimicking the presence of a repulsive
3647- area in the middle of the simulation box without needing to simulate additional atoms.
3648- The procedure is valid for more complex systems and can be adapted to many other
3649- situations, such as measuring adsorption barriers near an interface or calculating
3650- translocation barriers through a membrane~\cite {wilson1997adsorption , makarov2009computer ,
3651- gravelle2021adsorption , loche2022molecular , hayatifar2024probing }.
3648+ The objective of this tutorial is to measure the free energy profile of
3649+ particles through a barrier potential using two methods: free sampling
3650+ and umbrella sampling~\cite {kastner2011umbrella , allen2017computer ,
3651+ frenkel2023understanding } (Fig.~\ref {fig:US }). To simplify the
3652+ process and minimize computation time, the barrier potential will be
3653+ imposed on the atoms using an additional force, mimicking the presence
3654+ of a repulsive area in the middle of the simulation box without needing
3655+ to simulate additional atoms. The procedure is valid for more complex
3656+ systems and can be adapted to many other situations, such as measuring
3657+ adsorption barriers near an interface or calculating translocation
3658+ barriers through a membrane~\cite {wilson1997adsorption ,
3659+ makarov2009computer , gravelle2021adsorption , loche2022molecular ,
3660+ hayatifar2024probing }.
36523661
36533662\subsubsection {Method 1: Free sampling }
3654- The most direct way to calculate a free energy profile is to extract the partition
3655- function from a classical (i.e.~unbiased) molecular dynamics simulation, and
3656- then estimate the Gibbs free energy by using
3663+ The most direct way to calculate a free energy profile is to extract the
3664+ partition function from a classical (i.e.~unbiased) molecular dynamics
3665+ simulation, and then estimate the Gibbs free energy by using
36573666\begin {equation }
36583667\Delta G = -RT \ln (p/p_0),
36593668\label {eq:G }
0 commit comments