forked from lijiext/lammps
210 lines
12 KiB
Groff
210 lines
12 KiB
Groff
LAMMPS (5 Oct 2016)
|
|
# Testsystem for core-shell model compared to Mitchel and Finchham
|
|
# Hendrik Heenen, June 2014
|
|
|
|
# ------------------------ INITIALIZATION ----------------------------
|
|
|
|
units metal
|
|
dimension 3
|
|
boundary p p p
|
|
atom_style full
|
|
|
|
# ----------------------- ATOM DEFINITION ----------------------------
|
|
|
|
fix csinfo all property/atom i_CSID
|
|
read_data data.coreshell fix csinfo NULL CS-Info
|
|
orthogonal box = (0 0 0) to (24.096 24.096 24.096)
|
|
1 by 2 by 2 MPI processor grid
|
|
reading atoms ...
|
|
432 atoms
|
|
scanning bonds ...
|
|
1 = max bonds/atom
|
|
reading bonds ...
|
|
216 bonds
|
|
1 = max # of 1-2 neighbors
|
|
0 = max # of 1-3 neighbors
|
|
0 = max # of 1-4 neighbors
|
|
1 = max # of special neighbors
|
|
|
|
group cores type 1 2
|
|
216 atoms in group cores
|
|
group shells type 3 4
|
|
216 atoms in group shells
|
|
|
|
neighbor 2.0 bin
|
|
comm_modify vel yes
|
|
|
|
# ------------------------ FORCE FIELDS ------------------------------
|
|
|
|
kspace_style ewald 1.0e-6
|
|
pair_style born/coul/long/cs 20.0 20.0 # A, rho, sigma=0, C, D
|
|
pair_coeff * * 0.0 1.000 0.00 0.00 0.00
|
|
pair_coeff 3 3 487.0 0.23768 0.00 1.05 0.50 #Na-Na
|
|
pair_coeff 3 4 145134.0 0.23768 0.00 6.99 8.70 #Na-Cl
|
|
pair_coeff 4 4 405774.0 0.23768 0.00 72.40 145.40 #Cl-Cl
|
|
|
|
bond_style harmonic
|
|
bond_coeff 1 63.014 0.0
|
|
bond_coeff 2 25.724 0.0
|
|
|
|
# ------------------------ Equilibration Run -------------------------------
|
|
|
|
reset_timestep 0
|
|
|
|
thermo 50
|
|
thermo_style custom step etotal pe ke temp press epair evdwl ecoul elong ebond fnorm fmax vol
|
|
|
|
compute CSequ all temp/cs cores shells
|
|
|
|
# output via chunk method
|
|
|
|
#compute prop all property/atom i_CSID
|
|
#compute cs_chunk all chunk/atom c_prop
|
|
#compute cstherm all temp/chunk cs_chunk temp internal com yes cdof 3.0
|
|
#fix ave_chunk all ave/time 100 1 100 c_cstherm file chunk.dump mode vector
|
|
|
|
thermo_modify temp CSequ
|
|
|
|
# velocity bias option
|
|
|
|
velocity all create 1427 134 dist gaussian mom yes rot no bias yes temp CSequ
|
|
Ewald initialization ...
|
|
WARNING: Using 12-bit tables for long-range coulomb (../kspace.cpp:316)
|
|
G vector (1/distance) = 0.175257
|
|
estimated absolute RMS force accuracy = 1.49563e-05
|
|
estimated relative force accuracy = 1.03866e-06
|
|
KSpace vectors: actual max1d max3d = 257 5 665
|
|
kxmax kymax kzmax = 5 5 5
|
|
Neighbor list info ...
|
|
1 neighbor list requests
|
|
update every 1 steps, delay 10 steps, check yes
|
|
max neighbors/atom: 2000, page size: 100000
|
|
master list distance cutoff = 22
|
|
ghost atom cutoff = 22
|
|
binsize = 11 -> bins = 3 3 3
|
|
velocity all scale 1427 temp CSequ
|
|
|
|
fix thermoberendsen all temp/berendsen 1427 1427 0.4
|
|
fix nve all nve
|
|
fix_modify thermoberendsen temp CSequ
|
|
|
|
# 2 fmsec timestep
|
|
|
|
timestep 0.002
|
|
run 500
|
|
Ewald initialization ...
|
|
WARNING: Using 12-bit tables for long-range coulomb (../kspace.cpp:316)
|
|
G vector (1/distance) = 0.175257
|
|
estimated absolute RMS force accuracy = 1.49563e-05
|
|
estimated relative force accuracy = 1.03866e-06
|
|
KSpace vectors: actual max1d max3d = 257 5 665
|
|
kxmax kymax kzmax = 5 5 5
|
|
Memory usage per processor = 15.2111 Mbytes
|
|
Step TotEng PotEng KinEng Temp Press E_pair E_vdwl E_coul E_long E_bond Fnorm Fmax Volume
|
|
0 -635.44099 -675.09865 39.657659 1427 -20613.612 -675.09865 1.6320365 1018.8211 -1695.5518 0 2.9851021e-14 4.2188475e-15 13990.5
|
|
50 -634.17349 -666.17817 32.004686 1151.6234 -4030.346 -667.9562 37.809075 989.23283 -1694.9981 1.7780241 9.7133824 2.1917971 13990.5
|
|
100 -631.93518 -661.91229 29.977104 1078.665 -3172.6554 -664.06722 39.333371 991.36935 -1694.7699 2.1549316 11.09144 2.1704351 13990.5
|
|
150 -630.16821 -663.12421 32.956001 1185.8545 142.23746 -665.48039 46.253735 982.97455 -1694.7087 2.3561823 11.467862 2.1651502 13990.5
|
|
200 -628.62229 -663.64928 35.026994 1260.375 -1623.495 -665.96883 41.661825 986.81676 -1694.4474 2.3195535 10.597709 2.9538127 13990.5
|
|
250 -627.31146 -662.01843 34.706969 1248.8595 -1510.677 -664.50329 42.10159 987.90193 -1694.5068 2.4848607 10.392642 1.6933333 13990.5
|
|
300 -626.57662 -660.97479 34.398174 1237.7482 -1638.7013 -663.48072 41.653268 989.53239 -1694.6664 2.5059279 11.204558 2.2814875 13990.5
|
|
350 -625.91972 -665.3583 39.438581 1419.1169 -1352.428 -667.61665 41.5825 985.53386 -1694.733 2.2583465 10.831573 1.8963653 13990.5
|
|
400 -625.24854 -660.77581 35.527271 1278.3764 -424.57407 -663.12561 43.89225 987.75056 -1694.7684 2.3497967 11.523588 2.1775908 13990.5
|
|
450 -624.67739 -662.88749 38.210098 1374.9125 1107.8255 -665.43273 46.964166 982.35231 -1694.7492 2.5452345 11.00224 2.0018738 13990.5
|
|
500 -624.30227 -661.98003 37.677766 1355.7576 -999.27326 -664.60512 42.373091 987.58943 -1694.5676 2.6250858 11.151835 2.8270448 13990.5
|
|
Loop time of 2.40619 on 4 procs for 500 steps with 432 atoms
|
|
|
|
Performance: 35.907 ns/day, 0.668 hours/ns, 207.797 timesteps/s
|
|
99.7% CPU use with 4 MPI tasks x no OpenMP threads
|
|
|
|
MPI task timing breakdown:
|
|
Section | min time | avg time | max time |%varavg| %total
|
|
---------------------------------------------------------------
|
|
Pair | 1.7876 | 1.9016 | 2.0542 | 7.0 | 79.03
|
|
Bond | 0.00064754 | 0.0006901 | 0.00074625 | 0.1 | 0.03
|
|
Kspace | 0.16854 | 0.31596 | 0.42594 | 16.5 | 13.13
|
|
Neigh | 0.051607 | 0.051768 | 0.051931 | 0.1 | 2.15
|
|
Comm | 0.11933 | 0.12412 | 0.12769 | 0.9 | 5.16
|
|
Output | 0.00039911 | 0.00042087 | 0.00047565 | 0.2 | 0.02
|
|
Modify | 0.0059447 | 0.0063215 | 0.0067348 | 0.4 | 0.26
|
|
Other | | 0.005341 | | | 0.22
|
|
|
|
Nlocal: 108 ave 119 max 102 min
|
|
Histogram: 1 1 1 0 0 0 0 0 0 1
|
|
Nghost: 6516 ave 6545 max 6491 min
|
|
Histogram: 1 0 0 1 1 0 0 0 0 1
|
|
Neighs: 74337.8 ave 81248 max 69989 min
|
|
Histogram: 1 0 2 0 0 0 0 0 0 1
|
|
|
|
Total # of neighbors = 297351
|
|
Ave neighs/atom = 688.312
|
|
Ave special neighs/atom = 1
|
|
Neighbor list builds = 22
|
|
Dangerous builds = 0
|
|
|
|
unfix thermoberendsen
|
|
|
|
# ------------------------ Dynamic Run -------------------------------
|
|
|
|
run 1000
|
|
Ewald initialization ...
|
|
WARNING: Using 12-bit tables for long-range coulomb (../kspace.cpp:316)
|
|
G vector (1/distance) = 0.175257
|
|
estimated absolute RMS force accuracy = 1.49563e-05
|
|
estimated relative force accuracy = 1.03866e-06
|
|
KSpace vectors: actual max1d max3d = 257 5 665
|
|
kxmax kymax kzmax = 5 5 5
|
|
Memory usage per processor = 15.213 Mbytes
|
|
Step TotEng PotEng KinEng Temp Press E_pair E_vdwl E_coul E_long E_bond Fnorm Fmax Volume
|
|
500 -624.30215 -661.98003 37.677882 1355.7618 -999.26443 -664.60512 42.373091 987.58943 -1694.5676 2.6250858 11.151835 2.8270448 13990.5
|
|
550 -624.30702 -662.13315 37.826126 1361.096 330.71041 -664.76759 45.316349 984.7201 -1694.804 2.6344423 10.876902 1.9726348 13990.5
|
|
600 -624.31001 -661.82602 37.516015 1349.9373 -978.11503 -664.34272 42.382797 987.88159 -1694.6071 2.5166949 10.21416 1.787289 13990.5
|
|
650 -624.30516 -659.19843 34.893268 1255.5631 548.43865 -661.95754 46.556222 986.06829 -1694.582 2.7591048 13.16005 2.1511562 13990.5
|
|
700 -624.30926 -661.91108 37.601819 1353.0248 -175.1339 -664.4528 44.284398 985.834 -1694.5712 2.5417263 10.972065 1.9847042 13990.5
|
|
750 -624.30955 -660.4096 36.100048 1298.9866 1222.5932 -663.20741 47.511377 984.09471 -1694.8135 2.7978144 11.761903 2.1481489 13990.5
|
|
800 -624.30447 -660.30683 36.002357 1295.4714 270.6639 -663.25589 45.210122 986.25515 -1694.7212 2.9490566 11.50988 4.0037749 13990.5
|
|
850 -624.30774 -663.65149 39.343757 1415.7049 339.356 -666.13397 45.12874 983.38616 -1694.6489 2.4824719 11.01305 2.1691298 13990.5
|
|
900 -624.30311 -658.78647 34.483353 1240.8132 309.88983 -661.8322 45.836842 987.17541 -1694.8445 3.045732 12.635589 1.9594437 13990.5
|
|
950 -624.303 -663.59545 39.292454 1413.8588 -649.78872 -665.94289 43.141981 985.63193 -1694.7168 2.347436 11.012602 2.0685819 13990.5
|
|
1000 -624.3032 -660.59776 36.294555 1305.9856 -753.87453 -663.22012 42.930848 988.41708 -1694.568 2.6223608 10.526403 1.8887197 13990.5
|
|
1050 -624.3023 -660.7017 36.399405 1309.7584 790.88591 -663.6376 46.453299 984.4785 -1694.5694 2.9359001 11.210429 2.1822375 13990.5
|
|
1100 -624.30545 -662.76291 38.457457 1383.8132 -853.19698 -665.29416 42.533854 986.84955 -1694.6776 2.5312537 10.517674 2.6445601 13990.5
|
|
1150 -624.31539 -663.03811 38.722722 1393.3582 -380.13868 -665.58847 43.476852 985.81618 -1694.8815 2.5503626 10.818389 1.9559813 13990.5
|
|
1200 -624.3093 -659.85835 35.549055 1279.1603 44.140107 -662.93041 44.773497 986.67501 -1694.3789 3.0720506 11.310516 2.2601607 13990.5
|
|
1250 -624.30277 -663.46462 39.161852 1409.1594 -352.57297 -666.0144 43.477124 984.93408 -1694.4256 2.5497829 10.720231 1.8458373 13990.5
|
|
1300 -624.30329 -660.15937 35.85608 1290.2079 -348.83001 -662.73568 43.884677 988.07413 -1694.6945 2.5763161 11.369766 2.2580977 13990.5
|
|
1350 -624.30469 -662.72631 38.421619 1382.5236 1338.9038 -665.20399 47.35379 982.28307 -1694.8408 2.4776776 10.587886 1.7677665 13990.5
|
|
1400 -624.30499 -662.12513 37.820139 1360.8806 132.00801 -664.73589 44.55252 985.34275 -1694.6312 2.6107622 11.193448 2.1613441 13990.5
|
|
1450 -624.30207 -661.42856 37.126485 1335.9209 2069.8075 -664.07373 49.544418 981.20445 -1694.8226 2.6451763 13.011189 3.3987593 13990.5
|
|
1500 -624.29727 -656.99924 32.701969 1176.7137 504.56835 -660.08591 46.605885 987.7602 -1694.452 3.0866635 13.271242 3.025272 13990.5
|
|
Loop time of 4.96339 on 4 procs for 1000 steps with 432 atoms
|
|
|
|
Performance: 34.815 ns/day, 0.689 hours/ns, 201.475 timesteps/s
|
|
99.7% CPU use with 4 MPI tasks x no OpenMP threads
|
|
|
|
MPI task timing breakdown:
|
|
Section | min time | avg time | max time |%varavg| %total
|
|
---------------------------------------------------------------
|
|
Pair | 3.6491 | 3.8484 | 4.2642 | 12.4 | 77.54
|
|
Bond | 0.0013638 | 0.0014506 | 0.0016019 | 0.2 | 0.03
|
|
Kspace | 0.3347 | 0.73521 | 0.92903 | 27.3 | 14.81
|
|
Neigh | 0.10979 | 0.11008 | 0.11042 | 0.1 | 2.22
|
|
Comm | 0.23785 | 0.25289 | 0.2601 | 1.8 | 5.10
|
|
Output | 0.0008378 | 0.00088245 | 0.00098872 | 0.2 | 0.02
|
|
Modify | 0.0022187 | 0.0024077 | 0.0026662 | 0.3 | 0.05
|
|
Other | | 0.01206 | | | 0.24
|
|
|
|
Nlocal: 108 ave 112 max 102 min
|
|
Histogram: 1 0 0 0 1 0 0 0 0 2
|
|
Nghost: 6500.75 ave 6612 max 6392 min
|
|
Histogram: 1 0 1 0 0 0 0 1 0 1
|
|
Neighs: 74275 ave 77543 max 69706 min
|
|
Histogram: 1 0 0 0 0 1 0 1 0 1
|
|
|
|
Total # of neighbors = 297100
|
|
Ave neighs/atom = 687.731
|
|
Ave special neighs/atom = 1
|
|
Neighbor list builds = 45
|
|
Dangerous builds = 0
|
|
Total wall time: 0:00:07
|