@@ -349,6 +349,14 @@ <h2>Imposed forcing</h2>
349
349
a group name 'ethanol_pull' (still to be defined) and 'water'. The spring constant
350
350
of the harmonic potential is 1000 kJ/mol/nm2, and the requested distance between the center-of-mass
351
351
of the two groups is 2 nm along the x dimension.
352
+ < br > < br >
353
+ Copy as well the previously created topol.top file. Modify the
354
+ first lines to adapt the path to the force field folder:
355
+ < pre > < code > #include "../preparation/ff/gromos54a7_atb.ff/forcefield.itp"
356
+ #include "../preparation/ff/BIPQ_GROMACS_G54A7FF_allatom.itp"
357
+ #include "../preparation/ff/h2o.itp"</ code > </ pre >
358
+ </ p >
359
+ < p >
352
360
< br > < br >
353
361
Let us create an index file :
354
362
< pre > < code > gmx make_ndx -f ../preparation/nvt.gro -o index.ndx</ code > </ pre >
@@ -408,9 +416,7 @@ <h2>Imposed forcing</h2>
408
416
</ p >
409
417
< p >
410
418
Then, during the production run, the average distance between the 2 groups
411
- is measured over time. Let us plot the histogram of this distance,
412
- toghether with the distribution for longer production run (if you have a
413
- good computer, feel free to run longer production run than 50 ps):
419
+ is measured over time. Let us plot the histogram of this distance:
414
420
</ p >
415
421
< p >
416
422
< img src ="../figures/03-EthanolAdsorption/probability_position.png "
@@ -419,6 +425,12 @@ <h2>Imposed forcing</h2>
419
425
distance between the two center-of-mass. Short (50 ps) and
420
426
long (1.5 ns) runs are compared.
421
427
</ p >
428
+ < p >
429
+ < b > Note:</ b > A longer run has been added for comparison.
430
+ If you have a good computer, feel free to run longer production run than 50 ps.
431
+ </ p >
432
+
433
+
422
434
< p >
423
435
The distribution is not centered around x = 2 nm because of the interaction between the
424
436
CO2 and the rest of the fluid.
@@ -440,11 +452,12 @@ <h2>Free energy profile calculation</h2>
440
452
< a href ="https://raw.githubusercontent.com/gromacstutorials/gromacstutorials.github.io/main/inputs/03-EthanolAdsorption/adsorption/inputs/pro.mdp " download target ="_blank "> pro.mdp</ a > ,
441
453
in it.
442
454
< br > < br >
443
- The only difference with the previous input scripts, is the command:
455
+ The only difference with the previous input scripts is the command:
444
456
< pre > < code > pull-coord1-init = to_be_replaced</ code > </ pre >
445
457
</ p >
446
458
< p >
447
- where the keyword 'to_be_replaced' is to be systematicaly replaced within a bash script. Create a bash script, call it 'run.sh', and
459
+ where the keyword 'to_be_replaced' is to be systematically replaced using a bash script.
460
+ Create a bash script, call it 'run.sh', and
448
461
copy the following lines in it:
449
462
</ p >
450
463
< pre > < code > #!/bin/bash
0 commit comments