@@ -35,7 +35,7 @@ Basic LAMMPS parameters
35
35
To begin this tutorial, if you are using LAMMPS--GUI, select ``Start Tutorial 7 ``
36
36
from the ``Tutorials `` menu and follow the instructions. Alternatively, if you are
37
37
not using LAMMPS--GUI, create a new folder and add a file named
38
- **free-energy .lmp **. Open the file in a text editor and paste in the following
38
+ **free-sampling .lmp **. Open the file in a text editor and paste in the following
39
39
content:
40
40
41
41
.. code-block :: lammps
@@ -66,7 +66,7 @@ System creation and settings
66
66
----------------------------
67
67
68
68
Let us define the simulation box and randomly add atoms by addying the
69
- following lines to **free-energy .lmp **:
69
+ following lines to **free-sampling .lmp **:
70
70
71
71
.. code-block :: lammps
72
72
@@ -119,7 +119,7 @@ frequently overcome the energy barrier due to thermal agitation.
119
119
120
120
We impose the force :math: `F(x)` to the atoms in the simulation
121
121
using the ``fix addforce `` command. Add the following
122
- lines to **free-energy .lmp **:
122
+ lines to **free-sampling .lmp **:
123
123
124
124
.. code-block :: lammps
125
125
@@ -161,7 +161,7 @@ VMD-compatible trajectory to a file:
161
161
dump_modify viz1 backcolor white acolor 1 cyan adiam 1 3 boxcolor black
162
162
163
163
# Option 2
164
- dump viz2 all atom 50000 free-energy .lammpstrj
164
+ dump viz2 all atom 50000 free-sampling .lammpstrj
165
165
166
166
Finally, let us perform an equilibration of 50000 steps,
167
167
using a timestep of :math: `2 \,\text {fs}`, corresponding to a total duration of :math: `100 \,\text {ps}`:
@@ -194,10 +194,11 @@ Run and data acquisition
194
194
195
195
Once the system is equilibrated, we will record the density profile of
196
196
the atoms along the :math: `x`-axis using the ``ave/chunk `` command.
197
- Add the following line to **free-energy .lmp **:
197
+ Add the following line to **free-sampling .lmp **:
198
198
199
199
.. code-block :: lammps
200
200
201
+ # undump viz2 # Uncomment this line if you're using Option 2
201
202
reset_timestep 0
202
203
203
204
thermo 200000
@@ -287,7 +288,7 @@ unbiased free energy profile.
287
288
LAMMPS input script
288
289
-------------------
289
290
290
- If you are using LAMMPS--GUI, open the file named **free-energy .lmp **.
291
+ If you are using LAMMPS--GUI, open the file named **umbrella-sampling .lmp **.
291
292
Alternatively, if you are not using LAMMPS--GUI, create a new input file
292
293
and paste in the following content:
293
294
@@ -355,7 +356,7 @@ umbrella sampling run:
355
356
dump_modify viz1 backcolor white acolor 1 cyan acolor 2 red adiam 1 3 adiam 2 3 boxcolor black
356
357
357
358
# Option 2
358
- dump viz2 all atom 50000 free-energy .lammpstrj
359
+ dump viz2 all atom 50000 umbrella-sampling .lammpstrj
359
360
360
361
timestep 2.0
361
362
run 50000
0 commit comments