@@ -391,18 +391,19 @@ \subsection{About \lammpsgui{}}
391391
392392\lammpsgui {} is a graphical text editor, enhanced for editing LAMMPS
393393input files and linked to the LAMMPS library, allowing it to run LAMMPS
394- directly. The text editor functions similarly to other graphical
394+ directly. The text editor { \color {blue}is similar} to other graphical
395395editors, such as Notepad or Gedit, but offers the following enhancements
396- specifically for LAMMPS:
396+ for { \color {blue} running} LAMMPS:
397397\begin {itemize }
398398 \item Wizard dialogs to set up these tutorials
399399 \item Auto-completion of LAMMPS commands and options
400400 \item Context-sensitive online help
401401 \item Syntax highlighting for LAMMPS input files
402402 \item Syntax-aware line indentation
403+ \item {\color {blue}Editor switches working directory to that of input file}
403404 \item Visualization using LAMMPS' built-in renderer
404405 \item Start and stop simulations via mouse or keyboard
405- \item Monitoring of simulation progress
406+ \item Monitoring of simulation progress { \color {blue} and parallelization}
406407 \item Dynamic capture of LAMMPS output in a text window
407408 \item Automatic plotting of thermodynamic data during runs
408409 \item Capture of `` dump image'' outputs for animations
@@ -1695,7 +1696,7 @@ \subsubsection{Unbreakable bonds}
16951696
16961697\begin {figure }
16971698\centering
1698- \includegraphics [width=\linewidth ]{CNT-unbreakable-length-energy}
1699+ \includegraphics [width=\linewidth ]{CNT-unbreakable-length-energy}\\ [-2ex]
16991700\caption {a) Evolution of the length $ L_\text {cnt}$ of the CNT with time,
17001701as simulated during \hyperref [carbon-nanotube-label]{Tutorial 2}.
17011702The CNT starts deforming at $ t = 5 \, \text {ps}$ , and $ L_\text {cnt-0}$ is the
@@ -1738,7 +1739,7 @@ \subsubsection{Unbreakable bonds}
17381739
17391740\begin {figure }
17401741\centering
1741- \includegraphics [width=0.55\linewidth ]{CNT-unbreakable-stress-strain}
1742+ \includegraphics [width=0.55\linewidth ]{CNT-unbreakable-stress-strain}\\ [-2ex]
17421743\caption {Stress applied on the CNT during deformation, $ F_\text {cnt}/A_\text {cnt}$ ,
17431744where $ F_\text {cnt}$ is the force and $ A_\text {cnt}$ the CNT surface area,
17441745as a function of the strain, $ \Delta L_\text {cnt} = (L_\text {cnt}-L_\text {cnt-0}/L_\text {cnt-0})$ , where
@@ -1891,7 +1892,7 @@ \subsubsection{Breakable bonds}
18911892
18921893\begin {figure }
18931894\centering
1894- \includegraphics [width=\linewidth ]{CNT-breakable-stress-energy}
1895+ \includegraphics [width=\linewidth ]{CNT-breakable-stress-energy}\\ [-2ex]
18951896\caption {a) Evolution of the total energy $ E$ of the CNT with time $ t$ .
18961897b) Stress applied on the CNT during deformation, $ F_\text {cnt}/A_\text {cnt}$ ,
18971898where $ F_\text {cnt}$ is the force and $ A_\text {cnt}$ the CNT surface area,
@@ -2079,7 +2080,7 @@ \subsubsection{Preparing the water reservoir}
20792080
20802081\begin {figure }
20812082\centering
2082- \includegraphics [width=\linewidth ]{PEG-density}
2083+ \includegraphics [width=\linewidth ]{PEG-density}\\ [-2ex]
20832084\caption {a) Temperature, $ T$ , of the water reservoir from \hyperref [all-atom-label]{Tutorial 3}
20842085as a function of the time, $ t$ . The horizontal dashed line is the target temperature of 300\, K.
20852086b) Evolution of the system density, $ \rho $ , with $ t$ .}
@@ -2656,6 +2657,16 @@ \subsubsection{System preparation}
26562657delete_atoms random fraction 0.15 yes H2O NULL 482793 mol yes
26572658\end {lstlisting }
26582659
2660+ \begin {figure }
2661+ \centering
2662+ \includegraphics [width=\linewidth ]{NANOSHEAR-system}
2663+ \caption {Side view of the system. Periodic images are represented in darker colors.
2664+ Water molecules are in red and white, $ \text {Na}^+$ ions in purple, $ \text {Cl}^-$
2665+ ions in lime, and wall atoms in gray. Note the absence of atomic defect at the
2666+ cell boundaries.}
2667+ \label {fig:NANOSHEAR-system }
2668+ \end {figure }
2669+
26592670To create an image of the system, add the following \lmpcmd {dump} image
26602671into \flecmd {create.lmp} (see also Fig.~\ref {fig:NANOSHEAR-system }):
26612672\ begin{lstlisting}
@@ -2675,24 +2686,15 @@ \subsubsection{System preparation}
26752686
26762687write_data create.data nocoeff
26772688\end {lstlisting }
2678- The \lmpcmd {run 0} command runs the simulation for 0 steps, which {\color {blue} serves} for
2679- creating the system and saving its state. The \lmpcmd {write\_ data} command
2689+ The \lmpcmd {run 0} command {\color {blue}initializes} the simulation {\color {blue}but
2690+ does not advance positions or velocities which is required for cleanly saving the
2691+ state}. The \lmpcmd {write\_ data} command
26802692generates a file called \lmpcmd {system.data} containing the information required
26812693to restart the simulation from the final configuration produced by this input
26822694file. With the \lmpcmd {nocoeff} option, the parameters from the force field are
26832695not included in the \flecmd {.data} file. Run the \flecmd {create.lmp} file using LAMMPS,
26842696and a file named \flecmd {create.data} will be created alongside \flecmd {create.lmp}.
26852697
2686- \begin {figure }
2687- \centering
2688- \includegraphics [width=\linewidth ]{NANOSHEAR-system}
2689- \caption {Side view of the system. Periodic images are represented in darker colors.
2690- Water molecules are in red and white, $ \text {Na}^+$ ions in purple, $ \text {Cl}^-$
2691- ions in lime, and wall atoms in gray. Note the absence of atomic defect at the
2692- cell boundaries.}
2693- \label {fig:NANOSHEAR-system }
2694- \end {figure }
2695-
26962698\paragraph {Energy minimization }
26972699
26982700Let us move the atoms and place them in more energetically favorable positions
@@ -2703,10 +2705,9 @@ \subsubsection{System preparation}
27032705% like that performed in the first tutorial; \hyperref[lennard-jones-label]{Lennard-Jones fluid}.
27042706% Instead, we will conduct a molecular dynamics simulation, employing certain techniques
27052707% to prevent the system from exploding due to overlapping atoms.
2706-
27072708Open the \flecmd {equilibrate.lmp} file that was downloaded alongside
2708- \flecmd {create.lmp} during the tutorial setup. {\color {blue} Simiarly
2709- to before, it} contains the following lines:
2709+ \flecmd {create.lmp} during the tutorial setup. {\color {blue}Same as
2710+ before, it} contains the following lines:
27102711\ begin{lstlisting}
27112712boundary p p f
27122713units real
@@ -2854,6 +2855,17 @@ \subsubsection{System preparation}
28542855The \lmpcmd {undump} command is used to cancel the previous \lmpcmd {dump} command.
28552856Then, a new \lmpcmd {dump} command with a larger dumping period is used.
28562857
2858+ \begin {figure }
2859+ \centering
2860+ \includegraphics [width=\linewidth ]{NANOSHEAR-equilibration}\\ [-2ex]
2861+ \caption {a)~Pressure, $ p$ , of the nanosheared electrolyte system
2862+ simulated in \hyperref [sheared-confined-label]{Tutorial 4} as a function of the
2863+ time, $ t$ . b)~Distance between the walls, $ \Delta z$ , as a function of $ t$ .
2864+ \textcolor {blue}{The orange line
2865+ shows the raw data, and the blue line represents a time-averaged curve.}}
2866+ \label {fig:NANOSHEAR-equilibration }
2867+ \end {figure }
2868+
28572869\begin {note }
28582870{\color {blue}
28592871Just like the \lmpcmd {undump} command can cancel an active \lmpcmd {dump}, other
@@ -2877,17 +2889,6 @@ \subsubsection{System preparation}
28772889variables \lmpcmd {walltopz} and \lmpcmd {wallbotz}, i.e.~the distance between the
28782890two centers of mass of the walls.
28792891
2880- \begin {figure }
2881- \centering
2882- \includegraphics [width=\linewidth ]{NANOSHEAR-equilibration}\\ [-2ex]
2883- \caption {a)~Pressure, $ p$ , of the nanosheared electrolyte system
2884- simulated in \hyperref [sheared-confined-label]{Tutorial 4} as a function of the
2885- time, $ t$ . b)~Distance between the walls, $ \Delta z$ , as a function of $ t$ .
2886- \textcolor {blue}{The orange line
2887- shows the raw data, and the blue line represents a time-averaged curve.}}
2888- \label {fig:NANOSHEAR-equilibration }
2889- \end {figure }
2890-
28912892Finally, let us run the simulation for 30~ps by adding a \lmpcmd {run} command
28922893to \flecmd {equilibrate.lmp}:
28932894\ begin{lstlisting}
@@ -3632,15 +3633,6 @@ \subsubsection{Cracking the silica}
36323633thermo_style custom step temp etotal vol density
36333634\end {lstlisting }
36343635
3635- \begin {figure }
3636- \centering
3637- \includegraphics [width=\linewidth ]{GCMC-cracked}
3638- \caption {Block of silica from \hyperref [gcmc-silica-label]{Tutorial 6}
3639- after deformation. Silicon atoms are represented in yellow,
3640- and oxygen atoms in red. The crack was induced by the
3641- imposed deformation of the box along the $ x$ -axis (i.e.,~the horizontal axis).}
3642- \label {fig:GCMC-cracked }
3643- \end {figure }
36443636Let us progressively increase the size of the box in the $ x$ direction,
36453637forcing the silica to deform and eventually crack. To achive this,
36463638the \lmpcmd {fix deform} command is used, with a rate
@@ -3666,6 +3658,15 @@ \subsubsection{Cracking the silica}
36663658NVT ensemble, using the \lmpcmd {fix nvt} command does not guarantee that
36673659a simulation actually samples the NVT ensemble.}
36683660\end {note }
3661+ \begin {figure }
3662+ \centering
3663+ \includegraphics [width=\linewidth ]{GCMC-cracked}
3664+ \caption {Block of silica from \hyperref [gcmc-silica-label]{Tutorial 6}
3665+ after deformation. Silicon atoms are represented in yellow,
3666+ and oxygen atoms in red. The crack was induced by the
3667+ imposed deformation of the box along the $ x$ -axis (i.e.,~the horizontal axis).}
3668+ \label {fig:GCMC-cracked }
3669+ \end {figure }
36693670
36703671\subsubsection {Adding water }
36713672
@@ -3917,6 +3918,16 @@ \subsubsection{Adding water}
39173918{\color {blue}The \lmpcmd {f\_ } keywords extract the Monte Carlo move statistics
39183919output by the \lmpcmd {fix gcmc} command.}
39193920
3921+ \begin {figure }
3922+ \centering
3923+ \includegraphics [width=\linewidth ]{GCMC-solvated}
3924+ \caption {Snapshot of the silica system after the adsorption of water molecules
3925+ during \hyperref [gcmc-silica-label]{Tutorial 6}.
3926+ The oxygen atoms of the water molecules are represented in cyan, the silicon
3927+ atoms in yellow, and the oxygen atoms of the solid in red.}
3928+ \label {fig:GCMC-solvated }
3929+ \end {figure }
3930+
39203931\begin {note }
39213932When using the pressure argument, LAMMPS ignores the value of the
39223933chemical potential (here $ \mu = -0.5 \, \text {eV}$ , which corresponds roughly to
@@ -3935,16 +3946,6 @@ \subsubsection{Adding water}
39353946of successfully inserting a molecule. Here, the short simulation duration was
39363947made possible by the use of a high pressure.
39373948
3938- \begin {figure }
3939- \centering
3940- \includegraphics [width=\linewidth ]{GCMC-solvated}
3941- \caption {Snapshot of the silica system after the adsorption of water molecules
3942- during \hyperref [gcmc-silica-label]{Tutorial 6}.
3943- The oxygen atoms of the water molecules are represented in cyan, the silicon
3944- atoms in yellow, and the oxygen atoms of the solid in red.}
3945- \label {fig:GCMC-solvated }
3946- \end {figure }
3947-
39483949\subsection {Tutorial 7: Free energy calculation }
39493950\label {umbrella-sampling-label }
39503951
0 commit comments