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