|
1 | 1 | # 1) Initialization |
2 | | -units lj |
3 | | -dimension 3 |
4 | | -atom_style atomic |
5 | | -boundary p p p |
| 2 | +units lj |
| 3 | +dimension 3 |
| 4 | +atom_style atomic |
| 5 | +boundary p p p |
6 | 6 | # 2) System definition |
7 | | -pair_style lj/cut 4.0 |
8 | | -read_data min_coords.data |
9 | | -region cyl_in cylinder z 0 0 10 INF INF side in |
10 | | -region cyl_out cylinder z 0 0 10 INF INF side out |
11 | | -group grp_t1 type 1 |
12 | | -group grp_t2 type 2 |
13 | | -group grp_in region cyl_in |
14 | | -group grp_out region cyl_out |
15 | | -group grp_t1_in intersect grp_t1 grp_in |
16 | | -group grp_t2_out intersect grp_t2 grp_out |
17 | | -delete_atoms group grp_t1_in |
18 | | -delete_atoms group grp_t2_out |
| 7 | +pair_style lj/cut 4.0 |
| 8 | +read_data min_coords.data |
| 9 | +region cyl_in cylinder z 0 0 10 INF INF side in |
| 10 | +region cyl_out cylinder z 0 0 10 INF INF side out |
| 11 | +group grp_t1 type 1 |
| 12 | +group grp_t2 type 2 |
| 13 | +group grp_in region cyl_in |
| 14 | +group grp_out region cyl_out |
| 15 | +group grp_t1_in intersect grp_t1 grp_in |
| 16 | +group grp_t2_out intersect grp_t2 grp_out |
| 17 | +delete_atoms group grp_t1_in |
| 18 | +delete_atoms group grp_t2_out |
19 | 19 | # delete unnecessary groups |
20 | | -group grp_in delete |
21 | | -group grp_out delete |
22 | | -group grp_t1_in delete |
23 | | -group grp_t2_out delete |
24 | | -variable n1_in equal count(grp_t1,cyl_in) |
25 | | -variable n2_in equal count(grp_t2,cyl_in) |
26 | | -compute coor12 grp_t1 coord/atom cutoff 2.0 group grp_t2 |
27 | | -compute sumcoor12 grp_t1 reduce ave c_coor12 |
| 20 | +group grp_in delete |
| 21 | +group grp_out delete |
| 22 | +group grp_t1_in delete |
| 23 | +group grp_t2_out delete |
| 24 | +variable n1_in equal count(grp_t1,cyl_in) |
| 25 | +variable n2_in equal count(grp_t2,cyl_in) |
| 26 | +compute coor12 grp_t1 coord/atom cutoff 2.0 group grp_t2 |
| 27 | +compute sumcoor12 grp_t1 reduce ave c_coor12 |
28 | 28 | # 3) Settings |
29 | 29 | # 4) Visualization |
30 | | -thermo 1000 |
31 | | -thermo_style custom step temp pe ke etotal press v_n1_in v_n2_in c_sumcoor12 |
| 30 | +thermo 1000 |
| 31 | +thermo_style custom step temp pe ke etotal press v_n1_in v_n2_in c_sumcoor12 |
32 | 32 |
|
33 | | -dump mydmp all image 100 dump.md.*.ppm type type shiny 0.1 box no 0.01 view 0 0 zoom 1.8 fsaa yes |
34 | | -dump_modify mydmp adiam 1 1 adiam 2 3 acolor 1 turquoise acolor 2 royalblue backcolor white |
| 33 | +dump viz all image 100 myimage.md.*.ppm type type shiny 0.1 box no 0.01 view 0 0 zoom 1.8 fsaa yes |
| 34 | +dump_modify viz adiam 1 1 adiam 2 3 acolor 1 turquoise acolor 2 royalblue backcolor white |
35 | 35 |
|
36 | 36 | # 5) Run |
37 | | -velocity all create 1.0 49284 mom yes dist gaussian |
38 | | -fix mynve all nve |
39 | | -fix mylgv all langevin 1.0 1.0 0.1 1530917 zero yes |
40 | | -timestep 0.005 |
41 | | -run 300000 |
| 37 | +velocity all create 1.0 49284 mom yes dist gaussian |
| 38 | +fix mynve all nve |
| 39 | +fix mylgv all langevin 1.0 1.0 0.1 1530917 zero yes |
| 40 | +timestep 0.005 |
| 41 | +run 300000 |
0 commit comments