ICODE-ADC/lammps
4
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@7497 f3b2605a-c512-4ea7-a41b-209d697bcdaa

Browse Source
This commit is contained in:
sjplimp
2012-01-10 01:54:14 +00:00
parent 36ae86996e
commit 34c1930201
104 changed files with 30506 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

161
examples/ellipse/log.ellipse.gayberne.9Jan12.linux.1 Normal file
View File

@ -0,0 +1,161 @@
LAMMPS (10 Jan 2012)
# GayBerne ellipsoids in LJ background fluid
units lj
atom_style ellipsoid
dimension 2
lattice sq 0.02
Lattice spacing in x,y,z = 7.07107 7.07107 7.07107
region box block 0 20 0 20 -0.5 0.5
create_box 2 box
Created orthogonal box = (0 0 -3.53553) to (141.421 141.421 3.53553)
1 by 1 by 1 MPI processor grid
create_atoms 1 box
Created 400 atoms
set group all type/fraction 2 0.1 95392
45 settings made for type/fraction
set type 1 mass 1.0
355 settings made for mass
set type 2 mass 1.5
45 settings made for mass
set type 1 shape 1 1 1
355 settings made for shape
set type 2 shape 3 1 1
45 settings made for shape
set group all quat/random 18238
400 settings made for quat/random
compute rot all temp/asphere
group spheroid type 1
355 atoms in group spheroid
variable dof equal count(spheroid)+2
compute_modify rot extra ${dof}
compute_modify rot extra 357
velocity all create 2.4 87287 loop geom
pair_style gayberne 1.0 3.0 1.0 4.0
pair_coeff 1 1 3.0 1.0 1 1 1 1 1 1 2.5
pair_coeff 1 2 3.0 1.0 1 1 1 0 0 0
pair_coeff 2 2 1.0 1.0 1 1 0.2 0 0 0
#pair_style resquared 4.0
#pair_coeff 1 1 3.0 1 1 1 1 1 1 1 2.5
#pair_coeff 1 2 3.0 1 1 1 1 0 0 0
#pair_coeff 2 2 1.0 1 1 1 0.2 0 0 0
neighbor 0.8 bin
thermo_style custom step c_rot epair etotal press vol
thermo 100
timestep 0.002
compute q all property/atom quatw quati quatj quatk
dump 1 all custom 100 dump.ellipse.gayberne id type x y z c_q[1] c_q[2] c_q[3] c_q[4]
#dump 1 all custom 100 dump.ellipse.resquared # id type x y z c_q[1] c_q[2] c_q[3] c_q[4]
fix 1 all npt/asphere temp 2.0 2.0 0.1 iso 0.0 1.0 1.0 mtk no pchain 0 tchain 1
fix 2 all enforce2d
compute_modify 1_temp extra ${dof}
compute_modify 1_temp extra 357
# equilibrate to shrink box around dilute system
run 2000
Memory usage per processor = 3.71598 Mbytes
Step rot E_pair TotEng Press Volume
0 2.2718861 0 2.394 0.04788 20000
100 1.7443002 0 1.8380563 0.03576216 20558.672
200 2.2770454 0 2.3994366 0.046545139 20620.298
300 1.8573283 0 1.9571597 0.04240689 18460.771
400 2.1709963 -0.00050890768 2.2871782 0.066789189 13689.461
500 1.987484 -0.043660052 2.0342913 0.099199049 8475.4405
600 2.088535 -0.24318266 1.8270941 0.22391583 4674.8592
700 2.2604089 -0.73325761 1.481207 0.57184109 2507.095
800 1.8731398 -1.6161576 0.22080738 1.3049285 1420.4278
900 2.1404423 -3.0146759 -0.84039269 3.3566833 922.07169
1000 1.5927835 -4.5754144 -2.9697785 2.5493591 743.19289
1100 1.8380043 -4.8820159 -3.0604394 2.2950846 708.014
1200 2.05516 -4.5701279 -2.4718822 1.466469 753.52878
1300 2.0000364 -4.0755503 -2.031995 0.35675413 842.14053
1400 2.1055032 -3.7342781 -1.6632699 0.58193577 925.51606
1500 1.9241806 -3.582658 -1.639377 0.24973032 1006.7339
1600 2.125676 -3.4733166 -1.3572836 0.72069214 1061.1355
1700 1.8811354 -3.2933821 -1.4183258 0.59272021 1119.2646
1800 2.1014782 -3.1254259 -1.0248797 0.95254771 1162.1817
1900 1.8977855 -3.1905908 -1.2727444 0.73157369 1191.4876
2000 2.0794708 -3.253147 -1.2230906 0.4949131 1187.4615
Loop time of 0.744775 on 1 procs for 2000 steps with 400 atoms
Pair time (%) = 0.5056 (67.8863)
Neigh time (%) = 0.0121155 (1.62673)
Comm time (%) = 0.0133355 (1.79054)
Outpt time (%) = 0.0179069 (2.40434)
Other time (%) = 0.195817 (26.2921)
Nlocal: 400 ave 400 max 400 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost: 232 ave 232 max 232 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs: 2868 ave 2868 max 2868 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Total # of neighbors = 2868
Ave neighs/atom = 7.17
Neighbor list builds = 138
Dangerous builds = 80
# run dynamics on dense system
unfix 1
fix 1 all nve/asphere
run 2000
Memory usage per processor = 3.72613 Mbytes
Step rot E_pair TotEng Press Volume
2000 2.0794708 -3.253147 -1.2230906 0.4949131 1187.4615
2100 2.0350673 -3.2084761 -1.1823685 0.33981976 1187.4615
2200 1.9674878 -3.1385619 -1.1467435 1.1720135 1187.4615
2300 2.0160349 -3.191283 -1.1652127 0.98222385 1187.4615
2400 2.1323063 -3.3152815 -1.1669962 0.76099867 1187.4615
2500 2.0866912 -3.269205 -1.176606 0.67971578 1187.4615
2600 2.0982825 -3.2774432 -1.1869918 0.90454676 1187.4615
2700 2.0482779 -3.2229456 -1.1928884 0.69425859 1187.4615
2800 1.9676351 -3.1400429 -1.1685864 0.7577946 1187.4615
2900 2.0421413 -3.2182319 -1.1866955 0.79375204 1187.4615
3000 2.1074018 -3.2869882 -1.2013002 0.18882145 1187.4615
3100 2.0467315 -3.2228162 -1.1821143 0.70508986 1187.4615
3200 2.1165018 -3.2983333 -1.1590455 0.5195734 1187.4615
3300 2.1523232 -3.3384994 -1.1803176 0.45229381 1187.4615
3400 2.0966131 -3.2797219 -1.1608863 0.75030508 1187.4615
3500 2.0623035 -3.2433952 -1.1409405 0.78443505 1187.4615
3600 2.0110509 -3.1874756 -1.1629634 0.54103399 1187.4615
3700 2.0105927 -3.1835726 -1.1894401 0.27445727 1187.4615
3800 1.9528229 -3.1230718 -1.1908005 0.53926039 1187.4615
3900 2.0657038 -3.2420623 -1.1608259 0.56113492 1187.4615
4000 2.0354627 -3.2069944 -1.2073941 0.5668484 1187.4615
Loop time of 0.780395 on 1 procs for 2000 steps with 400 atoms
Pair time (%) = 0.622527 (79.7707)
Neigh time (%) = 0.00560331 (0.71801)
Comm time (%) = 0.0139475 (1.78723)
Outpt time (%) = 0.0179484 (2.29991)
Other time (%) = 0.120369 (15.4241)
Nlocal: 400 ave 400 max 400 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost: 222 ave 222 max 222 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs: 2914 ave 2914 max 2914 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Total # of neighbors = 2914
Ave neighs/atom = 7.285
Neighbor list builds = 47
Dangerous builds = 0

161
examples/ellipse/log.ellipse.gayberne.9Jan12.linux.4 Normal file
View File

@ -0,0 +1,161 @@
LAMMPS (10 Jan 2012)
# GayBerne ellipsoids in LJ background fluid
units lj
atom_style ellipsoid
dimension 2
lattice sq 0.02
Lattice spacing in x,y,z = 7.07107 7.07107 7.07107
region box block 0 20 0 20 -0.5 0.5
create_box 2 box
Created orthogonal box = (0 0 -3.53553) to (141.421 141.421 3.53553)
2 by 2 by 1 MPI processor grid
create_atoms 1 box
Created 400 atoms
set group all type/fraction 2 0.1 95392
45 settings made for type/fraction
set type 1 mass 1.0
355 settings made for mass
set type 2 mass 1.5
45 settings made for mass
set type 1 shape 1 1 1
355 settings made for shape
set type 2 shape 3 1 1
45 settings made for shape
set group all quat/random 18238
400 settings made for quat/random
compute rot all temp/asphere
group spheroid type 1
355 atoms in group spheroid
variable dof equal count(spheroid)+2
compute_modify rot extra ${dof}
compute_modify rot extra 357
velocity all create 2.4 87287 loop geom
pair_style gayberne 1.0 3.0 1.0 4.0
pair_coeff 1 1 3.0 1.0 1 1 1 1 1 1 2.5
pair_coeff 1 2 3.0 1.0 1 1 1 0 0 0
pair_coeff 2 2 1.0 1.0 1 1 0.2 0 0 0
#pair_style resquared 4.0
#pair_coeff 1 1 3.0 1 1 1 1 1 1 1 2.5
#pair_coeff 1 2 3.0 1 1 1 1 0 0 0
#pair_coeff 2 2 1.0 1 1 1 0.2 0 0 0
neighbor 0.8 bin
thermo_style custom step c_rot epair etotal press vol
thermo 100
timestep 0.002
compute q all property/atom quatw quati quatj quatk
dump 1 all custom 100 dump.ellipse.gayberne id type x y z c_q[1] c_q[2] c_q[3] c_q[4]
#dump 1 all custom 100 dump.ellipse.resquared # id type x y z c_q[1] c_q[2] c_q[3] c_q[4]
fix 1 all npt/asphere temp 2.0 2.0 0.1 iso 0.0 1.0 1.0 mtk no pchain 0 tchain 1
fix 2 all enforce2d
compute_modify 1_temp extra ${dof}
compute_modify 1_temp extra 357
# equilibrate to shrink box around dilute system
run 2000
Memory usage per processor = 3.70503 Mbytes
Step rot E_pair TotEng Press Volume
0 2.2718861 0 2.394 0.04788 20000
100 1.7443002 0 1.8380563 0.03576216 20558.672
200 2.2770454 0 2.3994366 0.046545139 20620.298
300 1.8573283 0 1.9571597 0.04240689 18460.771
400 2.1709963 -0.00050890768 2.2871782 0.066789189 13689.461
500 1.987484 -0.043660052 2.0342913 0.099199049 8475.4405
600 2.088535 -0.24318266 1.8270941 0.22391583 4674.8592
700 2.2604089 -0.73325761 1.481207 0.57184109 2507.095
800 1.8731398 -1.6161576 0.22080738 1.3049285 1420.4278
900 2.1404423 -3.0146759 -0.84039269 3.3566833 922.07169
1000 1.5927835 -4.5754144 -2.9697785 2.5493591 743.19289
1100 1.8380043 -4.8820159 -3.0604394 2.2950846 708.014
1200 2.05516 -4.5701279 -2.4718822 1.466469 753.52878
1300 2.0000364 -4.0755503 -2.031995 0.35675413 842.14053
1400 2.1055032 -3.7342781 -1.6632699 0.58193577 925.51606
1500 1.9241806 -3.582658 -1.639377 0.24973032 1006.7339
1600 2.125676 -3.4733166 -1.3572836 0.72069214 1061.1355
1700 1.8811354 -3.2933821 -1.4183258 0.59272021 1119.2646
1800 2.1014782 -3.1254259 -1.0248797 0.95254771 1162.1817
1900 1.8977855 -3.1905908 -1.2727444 0.73157369 1191.4876
2000 2.0794708 -3.253147 -1.2230906 0.49491311 1187.4615
Loop time of 0.276916 on 4 procs for 2000 steps with 400 atoms
Pair time (%) = 0.13188 (47.6245)
Neigh time (%) = 0.00327903 (1.18413)
Comm time (%) = 0.0477025 (17.2263)
Outpt time (%) = 0.0184636 (6.66759)
Other time (%) = 0.0755909 (27.2974)
Nlocal: 100 ave 113 max 79 min
Histogram: 1 0 0 0 0 1 0 0 1 1
Nghost: 139.75 ave 151 max 132 min
Histogram: 2 0 0 0 0 1 0 0 0 1
Neighs: 717 ave 880 max 487 min
Histogram: 1 0 0 0 0 1 0 1 0 1
Total # of neighbors = 2868
Ave neighs/atom = 7.17
Neighbor list builds = 138
Dangerous builds = 80
# run dynamics on dense system
unfix 1
fix 1 all nve/asphere
run 2000
Memory usage per processor = 3.70587 Mbytes
Step rot E_pair TotEng Press Volume
2000 2.0794708 -3.253147 -1.2230906 0.49491311 1187.4615
2100 2.0350673 -3.2084761 -1.1823685 0.33981977 1187.4615
2200 1.9674878 -3.1385619 -1.1467435 1.1720134 1187.4615
2300 2.0160348 -3.1912829 -1.1652127 0.98222397 1187.4615
2400 2.1323061 -3.3152813 -1.1669962 0.76099924 1187.4615
2500 2.0866909 -3.2692047 -1.1766059 0.67971739 1187.4615
2600 2.0982828 -3.2774435 -1.1869916 0.90454525 1187.4615
2700 2.0482779 -3.2229456 -1.1928882 0.69426441 1187.4615
2800 1.9676376 -3.1400474 -1.1685904 0.7577817 1187.4615
2900 2.0421364 -3.2182286 -1.1866999 0.79377489 1187.4615
3000 2.1074314 -3.2870216 -1.2013314 0.18876892 1187.4615
3100 2.0467964 -3.2228886 -1.1821044 0.70505639 1187.4615
3200 2.1167631 -3.298607 -1.1590547 0.51900634 1187.4615
3300 2.1528874 -3.339087 -1.1804887 0.45092936 1187.4615
3400 2.0979325 -3.2807362 -1.1607408 0.75035729 1187.4615
3500 2.0604373 -3.2416646 -1.1425421 0.7913285 1187.4615
3600 2.013392 -3.1897322 -1.1643385 0.56365647 1187.4615
3700 1.9922643 -3.1651575 -1.1883771 0.3074551 1187.4615
3800 1.9456148 -3.1160905 -1.1794166 0.73639236 1187.4615
3900 2.0237054 -3.1999536 -1.1869795 0.77500395 1187.4615
4000 2.0613642 -3.2364067 -1.2375068 0.53583637 1187.4615
Loop time of 0.294889 on 4 procs for 2000 steps with 400 atoms
Pair time (%) = 0.158921 (53.8917)
Neigh time (%) = 0.00147468 (0.500079)
Comm time (%) = 0.0758675 (25.7275)
Outpt time (%) = 0.0184799 (6.26672)
Other time (%) = 0.0401463 (13.614)
Nlocal: 100 ave 120 max 68 min
Histogram: 1 0 0 0 0 0 1 0 1 1
Nghost: 136 ave 151 max 126 min
Histogram: 2 0 0 0 0 0 1 0 0 1
Neighs: 732.5 ave 972 max 402 min
Histogram: 1 0 0 0 0 1 0 1 0 1
Total # of neighbors = 2930
Ave neighs/atom = 7.325
Neighbor list builds = 47
Dangerous builds = 0

161
examples/ellipse/log.ellipse.resquared.9Jan12.linux.1 Normal file
View File

@ -0,0 +1,161 @@
LAMMPS (10 Jan 2012)
# RESquared ellipsoids in LJ background fluid
units lj
atom_style ellipsoid
dimension 2
lattice sq 0.02
Lattice spacing in x,y,z = 7.07107 7.07107 7.07107
region box block 0 20 0 20 -0.5 0.5
create_box 2 box
Created orthogonal box = (0 0 -3.53553) to (141.421 141.421 3.53553)
1 by 1 by 1 MPI processor grid
create_atoms 1 box
Created 400 atoms
set group all type/fraction 2 0.1 95392
45 settings made for type/fraction
set type 1 mass 1.0
355 settings made for mass
set type 2 mass 1.5
45 settings made for mass
set type 1 shape 1 1 1
355 settings made for shape
set type 2 shape 3 1 1
45 settings made for shape
set group all quat/random 18238
400 settings made for quat/random
compute rot all temp/asphere
group spheroid type 1
355 atoms in group spheroid
variable dof equal count(spheroid)+2
compute_modify rot extra ${dof}
compute_modify rot extra 357
velocity all create 2.4 87287 loop geom
#pair_style gayberne 1.0 3.0 1.0 4.0
#pair_coeff 1 1 3.0 1.0 1 1 1 1 1 1 2.5
#pair_coeff 1 2 3.0 1.0 1 1 1 0 0 0
#pair_coeff 2 2 1.0 1.0 1 1 0.2 0 0 0
pair_style resquared 4.0
pair_coeff 1 1 3.0 1 1 1 1 1 1 1 2.5
pair_coeff 1 2 3.0 1 1 1 1 0 0 0
pair_coeff 2 2 1.0 1 1 1 0.2 0 0 0
neighbor 0.8 bin
thermo_style custom step c_rot epair etotal press vol
thermo 100
timestep 0.002
compute q all property/atom quatw quati quatj quatk
#dump 1 all custom 100 dump.ellipse.gayberne # id type x y z c_q[1] c_q[2] c_q[3] c_q[4]
dump 1 all custom 100 dump.ellipse.resquared id type x y z c_q[1] c_q[2] c_q[3] c_q[4]
fix 1 all npt/asphere temp 2.0 2.0 0.1 iso 0.0 1.0 1.0 mtk no pchain 0 tchain 1
fix 2 all enforce2d
compute_modify 1_temp extra ${dof}
compute_modify 1_temp extra 357
# equilibrate to shrink box around dilute system
run 2000
Memory usage per processor = 3.71598 Mbytes
Step rot E_pair TotEng Press Volume
0 2.2718861 0 2.394 0.04788 20000
100 1.7443002 0 1.8380563 0.03576216 20558.672
200 2.2770454 0 2.3994366 0.046545139 20620.298
300 1.8573283 0 1.9571597 0.04240689 18460.771
400 2.1708604 -2.7650654e-05 2.2875165 0.066837495 13689.498
500 1.997512 0.016920625 2.0564634 0.11517748 8482.9706
600 2.0622304 0.078679437 2.0821362 0.28827312 4776.2481
700 2.1106581 0.19687916 2.326279 0.85203631 2756.4602
800 1.745752 0.29267549 1.978786 1.6290595 1832.812
900 1.7925965 0.24064027 2.0273447 1.7295568 1516.7688
1000 1.8280597 0.16642381 1.9844257 1.4314533 1517.3656
1100 2.157879 0.10561715 2.2607067 0.98795205 1744.9669
1200 1.7966771 0.057760855 1.8309261 0.54369976 2182.8259
1300 2.2982472 0.036382605 2.3348878 0.45715755 2731.1166
1400 1.7666501 0.038109497 1.7892555 0.32977405 3230.78
1500 2.3048922 0.081153096 2.376275 0.45835895 3431.9991
1600 1.8588345 0.067622658 1.9190878 0.39494686 3276.3398
1700 2.1798404 0.1400767 2.3395072 0.71093539 2829.7905
1800 1.9307392 0.14286145 2.0646696 0.84046137 2320.9547
1900 2.0667728 0.16459518 2.1859494 1.1672332 1900.3276
2000 1.9311181 0.17205165 2.1209583 1.3796723 1657.9242
Loop time of 1.45024 on 1 procs for 2000 steps with 400 atoms
Pair time (%) = 1.21183 (83.5607)
Neigh time (%) = 0.0117629 (0.811095)
Comm time (%) = 0.00985646 (0.679641)
Outpt time (%) = 0.0178938 (1.23385)
Other time (%) = 0.198897 (13.7147)
Nlocal: 400 ave 400 max 400 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost: 207 ave 207 max 207 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs: 1886 ave 1886 max 1886 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Total # of neighbors = 1886
Ave neighs/atom = 4.715
Neighbor list builds = 179
Dangerous builds = 148
# run dynamics on dense system
unfix 1
fix 1 all nve/asphere
run 2000
Memory usage per processor = 3.718 Mbytes
Step rot E_pair TotEng Press Volume
2000 1.9311181 0.17205165 2.1209583 1.3796723 1657.9242
2100 1.9516972 0.15050034 2.1113621 1.1750685 1657.9242
2200 1.9410444 0.1617683 2.122019 1.2634696 1657.9242
2300 1.918524 0.185457 2.1121008 1.3486177 1657.9242
2400 1.9148973 0.18929207 2.1281553 1.3450222 1657.9242
2500 1.8892245 0.21635013 2.1016576 1.5081781 1657.9242
2600 1.9450941 0.15744769 2.0838764 1.2106283 1657.9242
2700 1.9229729 0.18066636 2.1170012 1.3384025 1657.9242
2800 1.9512215 0.15099708 2.0988078 1.2134332 1657.9242
2900 1.9211262 0.18279164 2.0994516 1.3620722 1657.9242
3000 1.9067873 0.19775142 2.0931298 1.4157494 1657.9242
3100 1.9114596 0.19296349 2.1126573 1.411347 1657.9242
3200 1.9427389 0.15998518 2.0865326 1.3075167 1657.9242
3300 1.9111615 0.19317456 2.1216596 1.4309787 1657.9242
3400 1.9538371 0.14833029 2.1210138 1.2276245 1657.9242
3500 1.9244333 0.17929251 2.082515 1.3186851 1657.9242
3600 1.9511825 0.15107526 2.1197982 1.2271872 1657.9242
3700 1.9369717 0.16601219 2.0970211 1.2812847 1657.9242
3800 1.9418085 0.16093887 2.1231359 1.2354724 1657.9242
3900 1.9450402 0.15760907 2.139479 1.2392343 1657.9242
4000 1.9350094 0.16819308 2.1142355 1.2676502 1657.9242
Loop time of 2.46236 on 1 procs for 2000 steps with 400 atoms
Pair time (%) = 2.30114 (93.4523)
Neigh time (%) = 0.00465512 (0.189051)
Comm time (%) = 0.0134063 (0.544447)
Outpt time (%) = 0.0178988 (0.726895)
Other time (%) = 0.125267 (5.08728)
Nlocal: 400 ave 400 max 400 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost: 209 ave 209 max 209 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs: 1903 ave 1903 max 1903 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Total # of neighbors = 1903
Ave neighs/atom = 4.7575
Neighbor list builds = 49
Dangerous builds = 0

161
examples/ellipse/log.ellipse.resquared.9Jan12.linux.4 Normal file
View File

@ -0,0 +1,161 @@
LAMMPS (10 Jan 2012)
# RESquared ellipsoids in LJ background fluid
units lj
atom_style ellipsoid
dimension 2
lattice sq 0.02
Lattice spacing in x,y,z = 7.07107 7.07107 7.07107
region box block 0 20 0 20 -0.5 0.5
create_box 2 box
Created orthogonal box = (0 0 -3.53553) to (141.421 141.421 3.53553)
2 by 2 by 1 MPI processor grid
create_atoms 1 box
Created 400 atoms
set group all type/fraction 2 0.1 95392
45 settings made for type/fraction
set type 1 mass 1.0
355 settings made for mass
set type 2 mass 1.5
45 settings made for mass
set type 1 shape 1 1 1
355 settings made for shape
set type 2 shape 3 1 1
45 settings made for shape
set group all quat/random 18238
400 settings made for quat/random
compute rot all temp/asphere
group spheroid type 1
355 atoms in group spheroid
variable dof equal count(spheroid)+2
compute_modify rot extra ${dof}
compute_modify rot extra 357
velocity all create 2.4 87287 loop geom
#pair_style gayberne 1.0 3.0 1.0 4.0
#pair_coeff 1 1 3.0 1.0 1 1 1 1 1 1 2.5
#pair_coeff 1 2 3.0 1.0 1 1 1 0 0 0
#pair_coeff 2 2 1.0 1.0 1 1 0.2 0 0 0
pair_style resquared 4.0
pair_coeff 1 1 3.0 1 1 1 1 1 1 1 2.5
pair_coeff 1 2 3.0 1 1 1 1 0 0 0
pair_coeff 2 2 1.0 1 1 1 0.2 0 0 0
neighbor 0.8 bin
thermo_style custom step c_rot epair etotal press vol
thermo 100
timestep 0.002
compute q all property/atom quatw quati quatj quatk
#dump 1 all custom 100 dump.ellipse.gayberne # id type x y z c_q[1] c_q[2] c_q[3] c_q[4]
dump 1 all custom 100 dump.ellipse.resquared id type x y z c_q[1] c_q[2] c_q[3] c_q[4]
fix 1 all npt/asphere temp 2.0 2.0 0.1 iso 0.0 1.0 1.0 mtk no pchain 0 tchain 1
fix 2 all enforce2d
compute_modify 1_temp extra ${dof}
compute_modify 1_temp extra 357
# equilibrate to shrink box around dilute system
run 2000
Memory usage per processor = 3.70503 Mbytes
Step rot E_pair TotEng Press Volume
0 2.2718861 0 2.394 0.04788 20000
100 1.7443002 0 1.8380563 0.03576216 20558.672
200 2.2770454 0 2.3994366 0.046545139 20620.298
300 1.8573283 0 1.9571597 0.04240689 18460.771
400 2.1708604 -2.7650654e-05 2.2875165 0.066837495 13689.498
500 1.997512 0.016920625 2.0564634 0.11517748 8482.9706
600 2.0622304 0.078679437 2.0821362 0.28827312 4776.2481
700 2.1106581 0.19687916 2.326279 0.85203631 2756.4602
800 1.745752 0.29267549 1.978786 1.6290595 1832.812
900 1.7925965 0.24064027 2.0273447 1.7295568 1516.7688
1000 1.8280597 0.16642381 1.9844257 1.4314533 1517.3656
1100 2.157879 0.10561715 2.2607067 0.98795205 1744.9669
1200 1.7966771 0.057760855 1.8309261 0.54369976 2182.8259
1300 2.2982472 0.036382605 2.3348878 0.45715755 2731.1166
1400 1.7666501 0.038109497 1.7892555 0.32977405 3230.78
1500 2.3048922 0.081153096 2.376275 0.45835895 3431.9991
1600 1.8588345 0.067622658 1.9190878 0.39494686 3276.3398
1700 2.1798404 0.1400767 2.3395072 0.71093539 2829.7905
1800 1.9307392 0.14286145 2.0646696 0.84046137 2320.9547
1900 2.0667728 0.16459518 2.1859494 1.1672332 1900.3276
2000 1.9311181 0.17205165 2.1209583 1.3796723 1657.9242
Loop time of 0.700645 on 4 procs for 2000 steps with 400 atoms
Pair time (%) = 0.315776 (45.0694)
Neigh time (%) = 0.00318927 (0.45519)
Comm time (%) = 0.202914 (28.961)
Outpt time (%) = 0.0198641 (2.83512)
Other time (%) = 0.158901 (22.6793)
Nlocal: 100 ave 109 max 96 min
Histogram: 2 0 1 0 0 0 0 0 0 1
Nghost: 116 ave 120 max 109 min
Histogram: 1 0 0 0 0 0 1 0 0 2
Neighs: 471.5 ave 499 max 455 min
Histogram: 2 0 0 0 1 0 0 0 0 1
Total # of neighbors = 1886
Ave neighs/atom = 4.715
Neighbor list builds = 179
Dangerous builds = 148
# run dynamics on dense system
unfix 1
fix 1 all nve/asphere
run 2000
Memory usage per processor = 3.70587 Mbytes
Step rot E_pair TotEng Press Volume
2000 1.9311181 0.17205165 2.1209583 1.3796723 1657.9242
2100 1.9516972 0.15050034 2.1113621 1.1750685 1657.9242
2200 1.9410444 0.1617683 2.122019 1.2634696 1657.9242
2300 1.918524 0.185457 2.1121008 1.3486177 1657.9242
2400 1.9148973 0.18929207 2.1281553 1.3450222 1657.9242
2500 1.8892245 0.21635014 2.1016576 1.5081781 1657.9242
2600 1.9450941 0.15744769 2.0838764 1.2106283 1657.9242
2700 1.9229729 0.18066636 2.1170012 1.3384026 1657.9242
2800 1.9512215 0.15099708 2.0988078 1.2134332 1657.9242
2900 1.9211262 0.18279165 2.0994516 1.3620722 1657.9242
3000 1.9067873 0.19775135 2.0931298 1.4157491 1657.9242
3100 1.9114597 0.19296344 2.1126573 1.4113468 1657.9242
3200 1.942739 0.15998509 2.0865326 1.3075162 1657.9242
3300 1.9111617 0.19317434 2.1216595 1.4309773 1657.9242
3400 1.9538374 0.14833002 2.121014 1.2276232 1657.9242
3500 1.9244329 0.17929295 2.0825158 1.3186873 1657.9242
3600 1.9511828 0.1510749 2.1197985 1.227186 1657.9242
3700 1.9369672 0.16601695 2.0970237 1.2813052 1657.9242
3800 1.9418088 0.16093858 2.1231366 1.2354749 1657.9242
3900 1.9450382 0.15761114 2.1394772 1.2392412 1657.9242
4000 1.9350119 0.16819038 2.1142388 1.2676423 1657.9242
Loop time of 0.718566 on 4 procs for 2000 steps with 400 atoms
Pair time (%) = 0.588593 (81.9121)
Neigh time (%) = 0.00123084 (0.171291)
Comm time (%) = 0.0711493 (9.90156)
Outpt time (%) = 0.0183669 (2.55604)
Other time (%) = 0.0392267 (5.45902)
Nlocal: 100 ave 107 max 97 min
Histogram: 1 2 0 0 0 0 0 0 0 1
Nghost: 117.5 ave 121 max 115 min
Histogram: 1 1 0 0 0 1 0 0 0 1
Neighs: 475.75 ave 511 max 446 min
Histogram: 1 0 0 1 1 0 0 0 0 1
Total # of neighbors = 1903
Ave neighs/atom = 4.7575
Neighbor list builds = 49
Dangerous builds = 0