LAMMPS (7 Jan 2022) using 1 OpenMP thread(s) per MPI task log 2.1_Initialize.out #------------------------------------------------------------------------------- # Stage 2.1: Initialize LAMMPS run for 3-d periodic #------------------------------------------------------------------------------- units real boundary p p p atom_style full pair_style lj/class2/coul/cut 9.5 pair_modify mix sixthpower pair_modify tail yes bond_style class2 angle_style class2 dihedral_style class2 improper_style class2 special_bonds lj 0.0 0.0 1.0 coul 0.0 0.0 1.0 box tilt large read_data structure.dat Reading data file ... orthogonal box = (0 0 0) to (31.042985 31.042985 31.042985) 2 by 4 by 4 MPI processor grid reading atoms ... 3000 atoms scanning bonds ... 1 = max bonds/atom scanning angles ... 1 = max angles/atom reading bonds ... 2000 bonds reading angles ... 1000 angles Finding 1-2 1-3 1-4 neighbors ... special bond factors lj: 0 0 1 special bond factors coul: 0 0 1 2 = max # of 1-2 neighbors 1 = max # of 1-3 neighbors 2 = max # of special neighbors special bonds CPU = 0.001 seconds read_data CPU = 0.020 seconds include parameters.dat bond_coeff 1 0.97 563.28 -1428.2 1902.1 # h* o* angle_coeff 1 103.7 49.84 -11.6 -8 # h* o* h* dihedral_coeff 1 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 # xxx xxx xxx xxx improper_coeff 1 0.0 0.0 # xxx xxx xxx xxx angle_coeff 1 bb -9.5 0.97 0.97 # h* o* h* angle_coeff 1 ba 22.35 22.35 0.97 0.97 # h* o* h* dihedral_coeff 1 mbt 0.0000 0.0000 0.0000 0.0 # xxx xxx xxx xxx dihedral_coeff 1 ebt 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0 0.0 # xxx xxx xxx xxx dihedral_coeff 1 at 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0 0.0 # xxx xxx xxx xxx dihedral_coeff 1 aat 0.0000 0.0 0.0 # xxx xxx xxx xxx dihedral_coeff 1 bb13 0.0 0.0 0.0 # xxx xxx xxx xxx improper_coeff 1 aa 0.0 0.0 0.0 120.0 120.0 120.0 # xxx xxx xxx xxx pair_coeff 1 1 0 1.098 # hw - pair_coeff 2 2 0.148 3.7824 # o* - neighbor 2.0 bin neigh_modify delay 0 every 1 check yes variable R equal 0.00198722 variable sysvol equal vol variable sysmass equal mass(all)/6.0221367e+23 variable sysdensity equal v_sysmass/v_sysvol/1.0e-24 variable coulomb equal ecoul+elong variable etotal equal etotal variable pe equal pe variable ke equal ke variable evdwl equal evdwl variable epair equal epair variable ebond equal ebond variable eangle equal eangle variable edihed equal edihed variable eimp equal eimp variable lx equal lx variable ly equal ly variable lz equal lz variable Nthermo equal 0 variable cella equal lx variable cellb equal sqrt(ly*ly+xy*xy) variable cellc equal sqrt(lz*lz+xz*xz+yz*yz) variable cellalpha equal acos((xy*xz+ly*yz)/(v_cellb*v_cellc)) variable cellbeta equal acos(xz/v_cellc) variable cellgamma equal acos(xy/v_cellb) variable p equal press variable pxx equal pxx variable pyy equal pyy variable pzz equal pzz variable pyz equal pyz variable pxz equal pxz variable pxy equal pxy variable sxx equal -pxx variable syy equal -pyy variable szz equal -pzz variable syz equal -pyz variable sxz equal -pxz variable sxy equal -pxy variable fmax equal fmax variable fnorm equal fnorm variable time equal step*dt+0.000001 variable surfacetension equal 0.5*v_lz*(0.5*(v_sxx+v_syy)-v_szz) thermo_style custom step v_time press vol v_sysdensity temp ebond eangle edihed eimp evdwl ecoul etail elong pe ke thermo_modify flush yes # # Set up the fixed and movable groups # group movable union all 3000 atoms in group movable group fixed subtract all movable 0 atoms in group fixed compute sum_f1 movable reduce sum fx fy fz variable sum_fx equal c_sum_f1[1] variable sum_fy equal c_sum_f1[2] variable sum_fz equal c_sum_f1[3] # # Subsets # group subset_O(H2O) id 1 4 7 10 13 16 19 22 25 28 10 atoms in group subset_O(H2O) group subset_O(H2O) id 31 34 37 40 43 46 49 52 55 58 20 atoms in group subset_O(H2O) group subset_O(H2O) id 61 64 67 70 73 76 79 82 85 88 30 atoms in group subset_O(H2O) group subset_O(H2O) id 91 94 97 100 103 106 109 112 115 118 40 atoms in group subset_O(H2O) group subset_O(H2O) id 121 124 127 130 133 136 139 142 145 148 50 atoms in group subset_O(H2O) group subset_O(H2O) id 151 154 157 160 163 166 169 172 175 178 60 atoms in group subset_O(H2O) group subset_O(H2O) id 181 184 187 190 193 196 199 202 205 208 70 atoms in group subset_O(H2O) group subset_O(H2O) id 211 214 217 220 223 226 229 232 235 238 80 atoms in group subset_O(H2O) group subset_O(H2O) id 241 244 247 250 253 256 259 262 265 268 90 atoms in group subset_O(H2O) group subset_O(H2O) id 271 274 277 280 283 286 289 292 295 298 100 atoms in group subset_O(H2O) group subset_O(H2O) id 301 304 307 310 313 316 319 322 325 328 110 atoms in group subset_O(H2O) group subset_O(H2O) id 331 334 337 340 343 346 349 352 355 358 120 atoms in group subset_O(H2O) group subset_O(H2O) id 361 364 367 370 373 376 379 382 385 388 130 atoms in group subset_O(H2O) group subset_O(H2O) id 391 394 397 400 403 406 409 412 415 418 140 atoms in group subset_O(H2O) group subset_O(H2O) id 421 424 427 430 433 436 439 442 445 448 150 atoms in group subset_O(H2O) group subset_O(H2O) id 451 454 457 460 463 466 469 472 475 478 160 atoms in group subset_O(H2O) group subset_O(H2O) id 481 484 487 490 493 496 499 502 505 508 170 atoms in group subset_O(H2O) group subset_O(H2O) id 511 514 517 520 523 526 529 532 535 538 180 atoms in group subset_O(H2O) group subset_O(H2O) id 541 544 547 550 553 556 559 562 565 568 190 atoms in group subset_O(H2O) group subset_O(H2O) id 571 574 577 580 583 586 589 592 595 598 200 atoms in group subset_O(H2O) group subset_O(H2O) id 601 604 607 610 613 616 619 622 625 628 210 atoms in group subset_O(H2O) group subset_O(H2O) id 631 634 637 640 643 646 649 652 655 658 220 atoms in group subset_O(H2O) group subset_O(H2O) id 661 664 667 670 673 676 679 682 685 688 230 atoms in group subset_O(H2O) group subset_O(H2O) id 691 694 697 700 703 706 709 712 715 718 240 atoms in group subset_O(H2O) group subset_O(H2O) id 721 724 727 730 733 736 739 742 745 748 250 atoms in group subset_O(H2O) group subset_O(H2O) id 751 754 757 760 763 766 769 772 775 778 260 atoms in group subset_O(H2O) group subset_O(H2O) id 781 784 787 790 793 796 799 802 805 808 270 atoms in group subset_O(H2O) group subset_O(H2O) id 811 814 817 820 823 826 829 832 835 838 280 atoms in group subset_O(H2O) group subset_O(H2O) id 841 844 847 850 853 856 859 862 865 868 290 atoms in group subset_O(H2O) group subset_O(H2O) id 871 874 877 880 883 886 889 892 895 898 300 atoms in group subset_O(H2O) group subset_O(H2O) id 901 904 907 910 913 916 919 922 925 928 310 atoms in group subset_O(H2O) group subset_O(H2O) id 931 934 937 940 943 946 949 952 955 958 320 atoms in group subset_O(H2O) group subset_O(H2O) id 961 964 967 970 973 976 979 982 985 988 330 atoms in group subset_O(H2O) group subset_O(H2O) id 991 994 997 1000 1003 1006 1009 1012 1015 1018 340 atoms in group subset_O(H2O) group subset_O(H2O) id 1021 1024 1027 1030 1033 1036 1039 1042 1045 1048 350 atoms in group subset_O(H2O) group subset_O(H2O) id 1051 1054 1057 1060 1063 1066 1069 1072 1075 1078 360 atoms in group subset_O(H2O) group subset_O(H2O) id 1081 1084 1087 1090 1093 1096 1099 1102 1105 1108 370 atoms in group subset_O(H2O) group subset_O(H2O) id 1111 1114 1117 1120 1123 1126 1129 1132 1135 1138 380 atoms in group subset_O(H2O) group subset_O(H2O) id 1141 1144 1147 1150 1153 1156 1159 1162 1165 1168 390 atoms in group subset_O(H2O) group subset_O(H2O) id 1171 1174 1177 1180 1183 1186 1189 1192 1195 1198 400 atoms in group subset_O(H2O) group subset_O(H2O) id 1201 1204 1207 1210 1213 1216 1219 1222 1225 1228 410 atoms in group subset_O(H2O) group subset_O(H2O) id 1231 1234 1237 1240 1243 1246 1249 1252 1255 1258 420 atoms in group subset_O(H2O) group subset_O(H2O) id 1261 1264 1267 1270 1273 1276 1279 1282 1285 1288 430 atoms in group subset_O(H2O) group subset_O(H2O) id 1291 1294 1297 1300 1303 1306 1309 1312 1315 1318 440 atoms in group subset_O(H2O) group subset_O(H2O) id 1321 1324 1327 1330 1333 1336 1339 1342 1345 1348 450 atoms in group subset_O(H2O) group subset_O(H2O) id 1351 1354 1357 1360 1363 1366 1369 1372 1375 1378 460 atoms in group subset_O(H2O) group subset_O(H2O) id 1381 1384 1387 1390 1393 1396 1399 1402 1405 1408 470 atoms in group subset_O(H2O) group subset_O(H2O) id 1411 1414 1417 1420 1423 1426 1429 1432 1435 1438 480 atoms in group subset_O(H2O) group subset_O(H2O) id 1441 1444 1447 1450 1453 1456 1459 1462 1465 1468 490 atoms in group subset_O(H2O) group subset_O(H2O) id 1471 1474 1477 1480 1483 1486 1489 1492 1495 1498 500 atoms in group subset_O(H2O) group subset_O(H2O) id 1501 1504 1507 1510 1513 1516 1519 1522 1525 1528 510 atoms in group subset_O(H2O) group subset_O(H2O) id 1531 1534 1537 1540 1543 1546 1549 1552 1555 1558 520 atoms in group subset_O(H2O) group subset_O(H2O) id 1561 1564 1567 1570 1573 1576 1579 1582 1585 1588 530 atoms in group subset_O(H2O) group subset_O(H2O) id 1591 1594 1597 1600 1603 1606 1609 1612 1615 1618 540 atoms in group subset_O(H2O) group subset_O(H2O) id 1621 1624 1627 1630 1633 1636 1639 1642 1645 1648 550 atoms in group subset_O(H2O) group subset_O(H2O) id 1651 1654 1657 1660 1663 1666 1669 1672 1675 1678 560 atoms in group subset_O(H2O) group subset_O(H2O) id 1681 1684 1687 1690 1693 1696 1699 1702 1705 1708 570 atoms in group subset_O(H2O) group subset_O(H2O) id 1711 1714 1717 1720 1723 1726 1729 1732 1735 1738 580 atoms in group subset_O(H2O) group subset_O(H2O) id 1741 1744 1747 1750 1753 1756 1759 1762 1765 1768 590 atoms in group subset_O(H2O) group subset_O(H2O) id 1771 1774 1777 1780 1783 1786 1789 1792 1795 1798 600 atoms in group subset_O(H2O) group subset_O(H2O) id 1801 1804 1807 1810 1813 1816 1819 1822 1825 1828 610 atoms in group subset_O(H2O) group subset_O(H2O) id 1831 1834 1837 1840 1843 1846 1849 1852 1855 1858 620 atoms in group subset_O(H2O) group subset_O(H2O) id 1861 1864 1867 1870 1873 1876 1879 1882 1885 1888 630 atoms in group subset_O(H2O) group subset_O(H2O) id 1891 1894 1897 1900 1903 1906 1909 1912 1915 1918 640 atoms in group subset_O(H2O) group subset_O(H2O) id 1921 1924 1927 1930 1933 1936 1939 1942 1945 1948 650 atoms in group subset_O(H2O) group subset_O(H2O) id 1951 1954 1957 1960 1963 1966 1969 1972 1975 1978 660 atoms in group subset_O(H2O) group subset_O(H2O) id 1981 1984 1987 1990 1993 1996 1999 2002 2005 2008 670 atoms in group subset_O(H2O) group subset_O(H2O) id 2011 2014 2017 2020 2023 2026 2029 2032 2035 2038 680 atoms in group subset_O(H2O) group subset_O(H2O) id 2041 2044 2047 2050 2053 2056 2059 2062 2065 2068 690 atoms in group subset_O(H2O) group subset_O(H2O) id 2071 2074 2077 2080 2083 2086 2089 2092 2095 2098 700 atoms in group subset_O(H2O) group subset_O(H2O) id 2101 2104 2107 2110 2113 2116 2119 2122 2125 2128 710 atoms in group subset_O(H2O) group subset_O(H2O) id 2131 2134 2137 2140 2143 2146 2149 2152 2155 2158 720 atoms in group subset_O(H2O) group subset_O(H2O) id 2161 2164 2167 2170 2173 2176 2179 2182 2185 2188 730 atoms in group subset_O(H2O) group subset_O(H2O) id 2191 2194 2197 2200 2203 2206 2209 2212 2215 2218 740 atoms in group subset_O(H2O) group subset_O(H2O) id 2221 2224 2227 2230 2233 2236 2239 2242 2245 2248 750 atoms in group subset_O(H2O) group subset_O(H2O) id 2251 2254 2257 2260 2263 2266 2269 2272 2275 2278 760 atoms in group subset_O(H2O) group subset_O(H2O) id 2281 2284 2287 2290 2293 2296 2299 2302 2305 2308 770 atoms in group subset_O(H2O) group subset_O(H2O) id 2311 2314 2317 2320 2323 2326 2329 2332 2335 2338 780 atoms in group subset_O(H2O) group subset_O(H2O) id 2341 2344 2347 2350 2353 2356 2359 2362 2365 2368 790 atoms in group subset_O(H2O) group subset_O(H2O) id 2371 2374 2377 2380 2383 2386 2389 2392 2395 2398 800 atoms in group subset_O(H2O) group subset_O(H2O) id 2401 2404 2407 2410 2413 2416 2419 2422 2425 2428 810 atoms in group subset_O(H2O) group subset_O(H2O) id 2431 2434 2437 2440 2443 2446 2449 2452 2455 2458 820 atoms in group subset_O(H2O) group subset_O(H2O) id 2461 2464 2467 2470 2473 2476 2479 2482 2485 2488 830 atoms in group subset_O(H2O) group subset_O(H2O) id 2491 2494 2497 2500 2503 2506 2509 2512 2515 2518 840 atoms in group subset_O(H2O) group subset_O(H2O) id 2521 2524 2527 2530 2533 2536 2539 2542 2545 2548 850 atoms in group subset_O(H2O) group subset_O(H2O) id 2551 2554 2557 2560 2563 2566 2569 2572 2575 2578 860 atoms in group subset_O(H2O) group subset_O(H2O) id 2581 2584 2587 2590 2593 2596 2599 2602 2605 2608 870 atoms in group subset_O(H2O) group subset_O(H2O) id 2611 2614 2617 2620 2623 2626 2629 2632 2635 2638 880 atoms in group subset_O(H2O) group subset_O(H2O) id 2641 2644 2647 2650 2653 2656 2659 2662 2665 2668 890 atoms in group subset_O(H2O) group subset_O(H2O) id 2671 2674 2677 2680 2683 2686 2689 2692 2695 2698 900 atoms in group subset_O(H2O) group subset_O(H2O) id 2701 2704 2707 2710 2713 2716 2719 2722 2725 2728 910 atoms in group subset_O(H2O) group subset_O(H2O) id 2731 2734 2737 2740 2743 2746 2749 2752 2755 2758 920 atoms in group subset_O(H2O) group subset_O(H2O) id 2761 2764 2767 2770 2773 2776 2779 2782 2785 2788 930 atoms in group subset_O(H2O) group subset_O(H2O) id 2791 2794 2797 2800 2803 2806 2809 2812 2815 2818 940 atoms in group subset_O(H2O) group subset_O(H2O) id 2821 2824 2827 2830 2833 2836 2839 2842 2845 2848 950 atoms in group subset_O(H2O) group subset_O(H2O) id 2851 2854 2857 2860 2863 2866 2869 2872 2875 2878 960 atoms in group subset_O(H2O) group subset_O(H2O) id 2881 2884 2887 2890 2893 2896 2899 2902 2905 2908 970 atoms in group subset_O(H2O) group subset_O(H2O) id 2911 2914 2917 2920 2923 2926 2929 2932 2935 2938 980 atoms in group subset_O(H2O) group subset_O(H2O) id 2941 2944 2947 2950 2953 2956 2959 2962 2965 2968 990 atoms in group subset_O(H2O) group subset_O(H2O) id 2971 2974 2977 2980 2983 2986 2989 2992 2995 2998 1000 atoms in group subset_O(H2O) group subset_H(H2O) id 2 3 5 6 8 9 11 12 14 15 10 atoms in group subset_H(H2O) group subset_H(H2O) id 17 18 20 21 23 24 26 27 29 30 20 atoms in group subset_H(H2O) group subset_H(H2O) id 32 33 35 36 38 39 41 42 44 45 30 atoms in group subset_H(H2O) group subset_H(H2O) id 47 48 50 51 53 54 56 57 59 60 40 atoms in group subset_H(H2O) group subset_H(H2O) id 62 63 65 66 68 69 71 72 74 75 50 atoms in group subset_H(H2O) group subset_H(H2O) id 77 78 80 81 83 84 86 87 89 90 60 atoms in group subset_H(H2O) group subset_H(H2O) id 92 93 95 96 98 99 101 102 104 105 70 atoms in group subset_H(H2O) group subset_H(H2O) id 107 108 110 111 113 114 116 117 119 120 80 atoms in group subset_H(H2O) group subset_H(H2O) id 122 123 125 126 128 129 131 132 134 135 90 atoms in group subset_H(H2O) group subset_H(H2O) id 137 138 140 141 143 144 146 147 149 150 100 atoms in group subset_H(H2O) group subset_H(H2O) id 152 153 155 156 158 159 161 162 164 165 110 atoms in group subset_H(H2O) group subset_H(H2O) id 167 168 170 171 173 174 176 177 179 180 120 atoms in group subset_H(H2O) group subset_H(H2O) id 182 183 185 186 188 189 191 192 194 195 130 atoms in group subset_H(H2O) group subset_H(H2O) id 197 198 200 201 203 204 206 207 209 210 140 atoms in group subset_H(H2O) group subset_H(H2O) id 212 213 215 216 218 219 221 222 224 225 150 atoms in group subset_H(H2O) group subset_H(H2O) id 227 228 230 231 233 234 236 237 239 240 160 atoms in group subset_H(H2O) group subset_H(H2O) id 242 243 245 246 248 249 251 252 254 255 170 atoms in group subset_H(H2O) group subset_H(H2O) id 257 258 260 261 263 264 266 267 269 270 180 atoms in group subset_H(H2O) group subset_H(H2O) id 272 273 275 276 278 279 281 282 284 285 190 atoms in group subset_H(H2O) group subset_H(H2O) id 287 288 290 291 293 294 296 297 299 300 200 atoms in group subset_H(H2O) group subset_H(H2O) id 302 303 305 306 308 309 311 312 314 315 210 atoms in group subset_H(H2O) group subset_H(H2O) id 317 318 320 321 323 324 326 327 329 330 220 atoms in group subset_H(H2O) group subset_H(H2O) id 332 333 335 336 338 339 341 342 344 345 230 atoms in group subset_H(H2O) group subset_H(H2O) id 347 348 350 351 353 354 356 357 359 360 240 atoms in group subset_H(H2O) group subset_H(H2O) id 362 363 365 366 368 369 371 372 374 375 250 atoms in group subset_H(H2O) group subset_H(H2O) id 377 378 380 381 383 384 386 387 389 390 260 atoms in group subset_H(H2O) group subset_H(H2O) id 392 393 395 396 398 399 401 402 404 405 270 atoms in group subset_H(H2O) group subset_H(H2O) id 407 408 410 411 413 414 416 417 419 420 280 atoms in group subset_H(H2O) group subset_H(H2O) id 422 423 425 426 428 429 431 432 434 435 290 atoms in group subset_H(H2O) group subset_H(H2O) id 437 438 440 441 443 444 446 447 449 450 300 atoms in group subset_H(H2O) group subset_H(H2O) id 452 453 455 456 458 459 461 462 464 465 310 atoms in group subset_H(H2O) group subset_H(H2O) id 467 468 470 471 473 474 476 477 479 480 320 atoms in group subset_H(H2O) group subset_H(H2O) id 482 483 485 486 488 489 491 492 494 495 330 atoms in group subset_H(H2O) group subset_H(H2O) id 497 498 500 501 503 504 506 507 509 510 340 atoms in group subset_H(H2O) group subset_H(H2O) id 512 513 515 516 518 519 521 522 524 525 350 atoms in group subset_H(H2O) group subset_H(H2O) id 527 528 530 531 533 534 536 537 539 540 360 atoms in group subset_H(H2O) group subset_H(H2O) id 542 543 545 546 548 549 551 552 554 555 370 atoms in group subset_H(H2O) group subset_H(H2O) id 557 558 560 561 563 564 566 567 569 570 380 atoms in group subset_H(H2O) group subset_H(H2O) id 572 573 575 576 578 579 581 582 584 585 390 atoms in group subset_H(H2O) group subset_H(H2O) id 587 588 590 591 593 594 596 597 599 600 400 atoms in group subset_H(H2O) group subset_H(H2O) id 602 603 605 606 608 609 611 612 614 615 410 atoms in group subset_H(H2O) group subset_H(H2O) id 617 618 620 621 623 624 626 627 629 630 420 atoms in group subset_H(H2O) group subset_H(H2O) id 632 633 635 636 638 639 641 642 644 645 430 atoms in group subset_H(H2O) group subset_H(H2O) id 647 648 650 651 653 654 656 657 659 660 440 atoms in group subset_H(H2O) group subset_H(H2O) id 662 663 665 666 668 669 671 672 674 675 450 atoms in group subset_H(H2O) group subset_H(H2O) id 677 678 680 681 683 684 686 687 689 690 460 atoms in group subset_H(H2O) group subset_H(H2O) id 692 693 695 696 698 699 701 702 704 705 470 atoms in group subset_H(H2O) group subset_H(H2O) id 707 708 710 711 713 714 716 717 719 720 480 atoms in group subset_H(H2O) group subset_H(H2O) id 722 723 725 726 728 729 731 732 734 735 490 atoms in group subset_H(H2O) group subset_H(H2O) id 737 738 740 741 743 744 746 747 749 750 500 atoms in group subset_H(H2O) group subset_H(H2O) id 752 753 755 756 758 759 761 762 764 765 510 atoms in group subset_H(H2O) group subset_H(H2O) id 767 768 770 771 773 774 776 777 779 780 520 atoms in group subset_H(H2O) group subset_H(H2O) id 782 783 785 786 788 789 791 792 794 795 530 atoms in group subset_H(H2O) group subset_H(H2O) id 797 798 800 801 803 804 806 807 809 810 540 atoms in group subset_H(H2O) group subset_H(H2O) id 812 813 815 816 818 819 821 822 824 825 550 atoms in group subset_H(H2O) group subset_H(H2O) id 827 828 830 831 833 834 836 837 839 840 560 atoms in group subset_H(H2O) group subset_H(H2O) id 842 843 845 846 848 849 851 852 854 855 570 atoms in group subset_H(H2O) group subset_H(H2O) id 857 858 860 861 863 864 866 867 869 870 580 atoms in group subset_H(H2O) group subset_H(H2O) id 872 873 875 876 878 879 881 882 884 885 590 atoms in group subset_H(H2O) group subset_H(H2O) id 887 888 890 891 893 894 896 897 899 900 600 atoms in group subset_H(H2O) group subset_H(H2O) id 902 903 905 906 908 909 911 912 914 915 610 atoms in group subset_H(H2O) group subset_H(H2O) id 917 918 920 921 923 924 926 927 929 930 620 atoms in group subset_H(H2O) group subset_H(H2O) id 932 933 935 936 938 939 941 942 944 945 630 atoms in group subset_H(H2O) group subset_H(H2O) id 947 948 950 951 953 954 956 957 959 960 640 atoms in group subset_H(H2O) group subset_H(H2O) id 962 963 965 966 968 969 971 972 974 975 650 atoms in group subset_H(H2O) group subset_H(H2O) id 977 978 980 981 983 984 986 987 989 990 660 atoms in group subset_H(H2O) group subset_H(H2O) id 992 993 995 996 998 999 1001 1002 1004 1005 670 atoms in group subset_H(H2O) group subset_H(H2O) id 1007 1008 1010 1011 1013 1014 1016 1017 1019 1020 680 atoms in group subset_H(H2O) group subset_H(H2O) id 1022 1023 1025 1026 1028 1029 1031 1032 1034 1035 690 atoms in group subset_H(H2O) group subset_H(H2O) id 1037 1038 1040 1041 1043 1044 1046 1047 1049 1050 700 atoms in group subset_H(H2O) group subset_H(H2O) id 1052 1053 1055 1056 1058 1059 1061 1062 1064 1065 710 atoms in group subset_H(H2O) group subset_H(H2O) id 1067 1068 1070 1071 1073 1074 1076 1077 1079 1080 720 atoms in group subset_H(H2O) group subset_H(H2O) id 1082 1083 1085 1086 1088 1089 1091 1092 1094 1095 730 atoms in group subset_H(H2O) group subset_H(H2O) id 1097 1098 1100 1101 1103 1104 1106 1107 1109 1110 740 atoms in group subset_H(H2O) group subset_H(H2O) id 1112 1113 1115 1116 1118 1119 1121 1122 1124 1125 750 atoms in group subset_H(H2O) group subset_H(H2O) id 1127 1128 1130 1131 1133 1134 1136 1137 1139 1140 760 atoms in group subset_H(H2O) group subset_H(H2O) id 1142 1143 1145 1146 1148 1149 1151 1152 1154 1155 770 atoms in group subset_H(H2O) group subset_H(H2O) id 1157 1158 1160 1161 1163 1164 1166 1167 1169 1170 780 atoms in group subset_H(H2O) group subset_H(H2O) id 1172 1173 1175 1176 1178 1179 1181 1182 1184 1185 790 atoms in group subset_H(H2O) group subset_H(H2O) id 1187 1188 1190 1191 1193 1194 1196 1197 1199 1200 800 atoms in group subset_H(H2O) group subset_H(H2O) id 1202 1203 1205 1206 1208 1209 1211 1212 1214 1215 810 atoms in group subset_H(H2O) group subset_H(H2O) id 1217 1218 1220 1221 1223 1224 1226 1227 1229 1230 820 atoms in group subset_H(H2O) group subset_H(H2O) id 1232 1233 1235 1236 1238 1239 1241 1242 1244 1245 830 atoms in group subset_H(H2O) group subset_H(H2O) id 1247 1248 1250 1251 1253 1254 1256 1257 1259 1260 840 atoms in group subset_H(H2O) group subset_H(H2O) id 1262 1263 1265 1266 1268 1269 1271 1272 1274 1275 850 atoms in group subset_H(H2O) group subset_H(H2O) id 1277 1278 1280 1281 1283 1284 1286 1287 1289 1290 860 atoms in group subset_H(H2O) group subset_H(H2O) id 1292 1293 1295 1296 1298 1299 1301 1302 1304 1305 870 atoms in group subset_H(H2O) group subset_H(H2O) id 1307 1308 1310 1311 1313 1314 1316 1317 1319 1320 880 atoms in group subset_H(H2O) group subset_H(H2O) id 1322 1323 1325 1326 1328 1329 1331 1332 1334 1335 890 atoms in group subset_H(H2O) group subset_H(H2O) id 1337 1338 1340 1341 1343 1344 1346 1347 1349 1350 900 atoms in group subset_H(H2O) group subset_H(H2O) id 1352 1353 1355 1356 1358 1359 1361 1362 1364 1365 910 atoms in group subset_H(H2O) group subset_H(H2O) id 1367 1368 1370 1371 1373 1374 1376 1377 1379 1380 920 atoms in group subset_H(H2O) group subset_H(H2O) id 1382 1383 1385 1386 1388 1389 1391 1392 1394 1395 930 atoms in group subset_H(H2O) group subset_H(H2O) id 1397 1398 1400 1401 1403 1404 1406 1407 1409 1410 940 atoms in group subset_H(H2O) group subset_H(H2O) id 1412 1413 1415 1416 1418 1419 1421 1422 1424 1425 950 atoms in group subset_H(H2O) group subset_H(H2O) id 1427 1428 1430 1431 1433 1434 1436 1437 1439 1440 960 atoms in group subset_H(H2O) group subset_H(H2O) id 1442 1443 1445 1446 1448 1449 1451 1452 1454 1455 970 atoms in group subset_H(H2O) group subset_H(H2O) id 1457 1458 1460 1461 1463 1464 1466 1467 1469 1470 980 atoms in group subset_H(H2O) group subset_H(H2O) id 1472 1473 1475 1476 1478 1479 1481 1482 1484 1485 990 atoms in group subset_H(H2O) group subset_H(H2O) id 1487 1488 1490 1491 1493 1494 1496 1497 1499 1500 1000 atoms in group subset_H(H2O) group subset_H(H2O) id 1502 1503 1505 1506 1508 1509 1511 1512 1514 1515 1010 atoms in group subset_H(H2O) group subset_H(H2O) id 1517 1518 1520 1521 1523 1524 1526 1527 1529 1530 1020 atoms in group subset_H(H2O) group subset_H(H2O) id 1532 1533 1535 1536 1538 1539 1541 1542 1544 1545 1030 atoms in group subset_H(H2O) group subset_H(H2O) id 1547 1548 1550 1551 1553 1554 1556 1557 1559 1560 1040 atoms in group subset_H(H2O) group subset_H(H2O) id 1562 1563 1565 1566 1568 1569 1571 1572 1574 1575 1050 atoms in group subset_H(H2O) group subset_H(H2O) id 1577 1578 1580 1581 1583 1584 1586 1587 1589 1590 1060 atoms in group subset_H(H2O) group subset_H(H2O) id 1592 1593 1595 1596 1598 1599 1601 1602 1604 1605 1070 atoms in group subset_H(H2O) group subset_H(H2O) id 1607 1608 1610 1611 1613 1614 1616 1617 1619 1620 1080 atoms in group subset_H(H2O) group subset_H(H2O) id 1622 1623 1625 1626 1628 1629 1631 1632 1634 1635 1090 atoms in group subset_H(H2O) group subset_H(H2O) id 1637 1638 1640 1641 1643 1644 1646 1647 1649 1650 1100 atoms in group subset_H(H2O) group subset_H(H2O) id 1652 1653 1655 1656 1658 1659 1661 1662 1664 1665 1110 atoms in group subset_H(H2O) group subset_H(H2O) id 1667 1668 1670 1671 1673 1674 1676 1677 1679 1680 1120 atoms in group subset_H(H2O) group subset_H(H2O) id 1682 1683 1685 1686 1688 1689 1691 1692 1694 1695 1130 atoms in group subset_H(H2O) group subset_H(H2O) id 1697 1698 1700 1701 1703 1704 1706 1707 1709 1710 1140 atoms in group subset_H(H2O) group subset_H(H2O) id 1712 1713 1715 1716 1718 1719 1721 1722 1724 1725 1150 atoms in group subset_H(H2O) group subset_H(H2O) id 1727 1728 1730 1731 1733 1734 1736 1737 1739 1740 1160 atoms in group subset_H(H2O) group subset_H(H2O) id 1742 1743 1745 1746 1748 1749 1751 1752 1754 1755 1170 atoms in group subset_H(H2O) group subset_H(H2O) id 1757 1758 1760 1761 1763 1764 1766 1767 1769 1770 1180 atoms in group subset_H(H2O) group subset_H(H2O) id 1772 1773 1775 1776 1778 1779 1781 1782 1784 1785 1190 atoms in group subset_H(H2O) group subset_H(H2O) id 1787 1788 1790 1791 1793 1794 1796 1797 1799 1800 1200 atoms in group subset_H(H2O) group subset_H(H2O) id 1802 1803 1805 1806 1808 1809 1811 1812 1814 1815 1210 atoms in group subset_H(H2O) group subset_H(H2O) id 1817 1818 1820 1821 1823 1824 1826 1827 1829 1830 1220 atoms in group subset_H(H2O) group subset_H(H2O) id 1832 1833 1835 1836 1838 1839 1841 1842 1844 1845 1230 atoms in group subset_H(H2O) group subset_H(H2O) id 1847 1848 1850 1851 1853 1854 1856 1857 1859 1860 1240 atoms in group subset_H(H2O) group subset_H(H2O) id 1862 1863 1865 1866 1868 1869 1871 1872 1874 1875 1250 atoms in group subset_H(H2O) group subset_H(H2O) id 1877 1878 1880 1881 1883 1884 1886 1887 1889 1890 1260 atoms in group subset_H(H2O) group subset_H(H2O) id 1892 1893 1895 1896 1898 1899 1901 1902 1904 1905 1270 atoms in group subset_H(H2O) group subset_H(H2O) id 1907 1908 1910 1911 1913 1914 1916 1917 1919 1920 1280 atoms in group subset_H(H2O) group subset_H(H2O) id 1922 1923 1925 1926 1928 1929 1931 1932 1934 1935 1290 atoms in group subset_H(H2O) group subset_H(H2O) id 1937 1938 1940 1941 1943 1944 1946 1947 1949 1950 1300 atoms in group subset_H(H2O) group subset_H(H2O) id 1952 1953 1955 1956 1958 1959 1961 1962 1964 1965 1310 atoms in group subset_H(H2O) group subset_H(H2O) id 1967 1968 1970 1971 1973 1974 1976 1977 1979 1980 1320 atoms in group subset_H(H2O) group subset_H(H2O) id 1982 1983 1985 1986 1988 1989 1991 1992 1994 1995 1330 atoms in group subset_H(H2O) group subset_H(H2O) id 1997 1998 2000 2001 2003 2004 2006 2007 2009 2010 1340 atoms in group subset_H(H2O) group subset_H(H2O) id 2012 2013 2015 2016 2018 2019 2021 2022 2024 2025 1350 atoms in group subset_H(H2O) group subset_H(H2O) id 2027 2028 2030 2031 2033 2034 2036 2037 2039 2040 1360 atoms in group subset_H(H2O) group subset_H(H2O) id 2042 2043 2045 2046 2048 2049 2051 2052 2054 2055 1370 atoms in group subset_H(H2O) group subset_H(H2O) id 2057 2058 2060 2061 2063 2064 2066 2067 2069 2070 1380 atoms in group subset_H(H2O) group subset_H(H2O) id 2072 2073 2075 2076 2078 2079 2081 2082 2084 2085 1390 atoms in group subset_H(H2O) group subset_H(H2O) id 2087 2088 2090 2091 2093 2094 2096 2097 2099 2100 1400 atoms in group subset_H(H2O) group subset_H(H2O) id 2102 2103 2105 2106 2108 2109 2111 2112 2114 2115 1410 atoms in group subset_H(H2O) group subset_H(H2O) id 2117 2118 2120 2121 2123 2124 2126 2127 2129 2130 1420 atoms in group subset_H(H2O) group subset_H(H2O) id 2132 2133 2135 2136 2138 2139 2141 2142 2144 2145 1430 atoms in group subset_H(H2O) group subset_H(H2O) id 2147 2148 2150 2151 2153 2154 2156 2157 2159 2160 1440 atoms in group subset_H(H2O) group subset_H(H2O) id 2162 2163 2165 2166 2168 2169 2171 2172 2174 2175 1450 atoms in group subset_H(H2O) group subset_H(H2O) id 2177 2178 2180 2181 2183 2184 2186 2187 2189 2190 1460 atoms in group subset_H(H2O) group subset_H(H2O) id 2192 2193 2195 2196 2198 2199 2201 2202 2204 2205 1470 atoms in group subset_H(H2O) group subset_H(H2O) id 2207 2208 2210 2211 2213 2214 2216 2217 2219 2220 1480 atoms in group subset_H(H2O) group subset_H(H2O) id 2222 2223 2225 2226 2228 2229 2231 2232 2234 2235 1490 atoms in group subset_H(H2O) group subset_H(H2O) id 2237 2238 2240 2241 2243 2244 2246 2247 2249 2250 1500 atoms in group subset_H(H2O) group subset_H(H2O) id 2252 2253 2255 2256 2258 2259 2261 2262 2264 2265 1510 atoms in group subset_H(H2O) group subset_H(H2O) id 2267 2268 2270 2271 2273 2274 2276 2277 2279 2280 1520 atoms in group subset_H(H2O) group subset_H(H2O) id 2282 2283 2285 2286 2288 2289 2291 2292 2294 2295 1530 atoms in group subset_H(H2O) group subset_H(H2O) id 2297 2298 2300 2301 2303 2304 2306 2307 2309 2310 1540 atoms in group subset_H(H2O) group subset_H(H2O) id 2312 2313 2315 2316 2318 2319 2321 2322 2324 2325 1550 atoms in group subset_H(H2O) group subset_H(H2O) id 2327 2328 2330 2331 2333 2334 2336 2337 2339 2340 1560 atoms in group subset_H(H2O) group subset_H(H2O) id 2342 2343 2345 2346 2348 2349 2351 2352 2354 2355 1570 atoms in group subset_H(H2O) group subset_H(H2O) id 2357 2358 2360 2361 2363 2364 2366 2367 2369 2370 1580 atoms in group subset_H(H2O) group subset_H(H2O) id 2372 2373 2375 2376 2378 2379 2381 2382 2384 2385 1590 atoms in group subset_H(H2O) group subset_H(H2O) id 2387 2388 2390 2391 2393 2394 2396 2397 2399 2400 1600 atoms in group subset_H(H2O) group subset_H(H2O) id 2402 2403 2405 2406 2408 2409 2411 2412 2414 2415 1610 atoms in group subset_H(H2O) group subset_H(H2O) id 2417 2418 2420 2421 2423 2424 2426 2427 2429 2430 1620 atoms in group subset_H(H2O) group subset_H(H2O) id 2432 2433 2435 2436 2438 2439 2441 2442 2444 2445 1630 atoms in group subset_H(H2O) group subset_H(H2O) id 2447 2448 2450 2451 2453 2454 2456 2457 2459 2460 1640 atoms in group subset_H(H2O) group subset_H(H2O) id 2462 2463 2465 2466 2468 2469 2471 2472 2474 2475 1650 atoms in group subset_H(H2O) group subset_H(H2O) id 2477 2478 2480 2481 2483 2484 2486 2487 2489 2490 1660 atoms in group subset_H(H2O) group subset_H(H2O) id 2492 2493 2495 2496 2498 2499 2501 2502 2504 2505 1670 atoms in group subset_H(H2O) group subset_H(H2O) id 2507 2508 2510 2511 2513 2514 2516 2517 2519 2520 1680 atoms in group subset_H(H2O) group subset_H(H2O) id 2522 2523 2525 2526 2528 2529 2531 2532 2534 2535 1690 atoms in group subset_H(H2O) group subset_H(H2O) id 2537 2538 2540 2541 2543 2544 2546 2547 2549 2550 1700 atoms in group subset_H(H2O) group subset_H(H2O) id 2552 2553 2555 2556 2558 2559 2561 2562 2564 2565 1710 atoms in group subset_H(H2O) group subset_H(H2O) id 2567 2568 2570 2571 2573 2574 2576 2577 2579 2580 1720 atoms in group subset_H(H2O) group subset_H(H2O) id 2582 2583 2585 2586 2588 2589 2591 2592 2594 2595 1730 atoms in group subset_H(H2O) group subset_H(H2O) id 2597 2598 2600 2601 2603 2604 2606 2607 2609 2610 1740 atoms in group subset_H(H2O) group subset_H(H2O) id 2612 2613 2615 2616 2618 2619 2621 2622 2624 2625 1750 atoms in group subset_H(H2O) group subset_H(H2O) id 2627 2628 2630 2631 2633 2634 2636 2637 2639 2640 1760 atoms in group subset_H(H2O) group subset_H(H2O) id 2642 2643 2645 2646 2648 2649 2651 2652 2654 2655 1770 atoms in group subset_H(H2O) group subset_H(H2O) id 2657 2658 2660 2661 2663 2664 2666 2667 2669 2670 1780 atoms in group subset_H(H2O) group subset_H(H2O) id 2672 2673 2675 2676 2678 2679 2681 2682 2684 2685 1790 atoms in group subset_H(H2O) group subset_H(H2O) id 2687 2688 2690 2691 2693 2694 2696 2697 2699 2700 1800 atoms in group subset_H(H2O) group subset_H(H2O) id 2702 2703 2705 2706 2708 2709 2711 2712 2714 2715 1810 atoms in group subset_H(H2O) group subset_H(H2O) id 2717 2718 2720 2721 2723 2724 2726 2727 2729 2730 1820 atoms in group subset_H(H2O) group subset_H(H2O) id 2732 2733 2735 2736 2738 2739 2741 2742 2744 2745 1830 atoms in group subset_H(H2O) group subset_H(H2O) id 2747 2748 2750 2751 2753 2754 2756 2757 2759 2760 1840 atoms in group subset_H(H2O) group subset_H(H2O) id 2762 2763 2765 2766 2768 2769 2771 2772 2774 2775 1850 atoms in group subset_H(H2O) group subset_H(H2O) id 2777 2778 2780 2781 2783 2784 2786 2787 2789 2790 1860 atoms in group subset_H(H2O) group subset_H(H2O) id 2792 2793 2795 2796 2798 2799 2801 2802 2804 2805 1870 atoms in group subset_H(H2O) group subset_H(H2O) id 2807 2808 2810 2811 2813 2814 2816 2817 2819 2820 1880 atoms in group subset_H(H2O) group subset_H(H2O) id 2822 2823 2825 2826 2828 2829 2831 2832 2834 2835 1890 atoms in group subset_H(H2O) group subset_H(H2O) id 2837 2838 2840 2841 2843 2844 2846 2847 2849 2850 1900 atoms in group subset_H(H2O) group subset_H(H2O) id 2852 2853 2855 2856 2858 2859 2861 2862 2864 2865 1910 atoms in group subset_H(H2O) group subset_H(H2O) id 2867 2868 2870 2871 2873 2874 2876 2877 2879 2880 1920 atoms in group subset_H(H2O) group subset_H(H2O) id 2882 2883 2885 2886 2888 2889 2891 2892 2894 2895 1930 atoms in group subset_H(H2O) group subset_H(H2O) id 2897 2898 2900 2901 2903 2904 2906 2907 2909 2910 1940 atoms in group subset_H(H2O) group subset_H(H2O) id 2912 2913 2915 2916 2918 2919 2921 2922 2924 2925 1950 atoms in group subset_H(H2O) group subset_H(H2O) id 2927 2928 2930 2931 2933 2934 2936 2937 2939 2940 1960 atoms in group subset_H(H2O) group subset_H(H2O) id 2942 2943 2945 2946 2948 2949 2951 2952 2954 2955 1970 atoms in group subset_H(H2O) group subset_H(H2O) id 2957 2958 2960 2961 2963 2964 2966 2967 2969 2970 1980 atoms in group subset_H(H2O) group subset_H(H2O) id 2972 2973 2975 2976 2978 2979 2981 2982 2984 2985 1990 atoms in group subset_H(H2O) group subset_H(H2O) id 2987 2988 2990 2991 2993 2994 2996 2997 2999 3000 2000 atoms in group subset_H(H2O) log 2.2_Minimization.out #------------------------------------------------------------------------------- # Stage 2.2: Minimization #------------------------------------------------------------------------------- min_style cg min_modify dmax 0.05 line fast reset_timestep 0 thermo_style custom step fmax fnorm press vol v_sysdensity v_sxx v_syy v_szz v_syz v_sxz v_sxy pe v_cella v_cellb v_cellc v_cellalpha v_cellbeta v_cellgamma v_sum_fx v_sum_fy v_sum_fz WARNING: New thermo_style command, previous thermo_modify settings will be lost (../output.cpp:888) dump sci all custom 1000 2.2.xyz id mol type q xs ys zs thermo 100 minimize 0.0 1.0 1000 10000 generated 0 of 1 mixed pair_coeff terms from sixthpower mixing rule Neighbor list info ... update every 1 steps, delay 0 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 11.5 ghost atom cutoff = 11.5 binsize = 5.75, bins = 6 6 6 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair lj/class2/coul/cut, perpetual attributes: half, newton on pair build: half/bin/newton stencil: half/bin/3d bin: standard Setting up cg style minimization ... Unit style : real Current step : 0 Per MPI rank memory allocation (min/avg/max) = 9.037 | 9.042 | 9.048 Mbytes Step Fmax Fnorm Press Volume v_sysdensity v_sxx v_syy v_szz v_syz v_sxz v_sxy PotEng v_cella v_cellb v_cellc v_cellalpha v_cellbeta v_cellgamma v_sum_fx v_sum_fy v_sum_fz 0 1022.784 5447.0541 -20801.497 29915.098 1.0000033 20085.95 20900.914 21417.625 5014.9525 2027.8678 2306.8614 15356.302 31.042985 31.042985 31.042985 1.5707963 1.5707963 1.5707963 2.7284841e-12 -1.4779289e-12 -1.4779289e-12 100 11.058773 172.25016 -2958.6469 29915.098 1.0000033 3084.9697 3012.3335 2778.6377 -1518.8443 -869.18988 -254.99649 -7932.8753 31.042985 31.042985 31.042985 1.5707963 1.5707963 1.5707963 2.7711167e-13 1.3216095e-12 -9.8054898e-13 200 8.2137423 110.28754 -4926.5349 29915.098 1.0000033 5739.1421 3804.3808 5236.0817 -157.46257 800.50481 -1561.1082 -9865.8991 31.042985 31.042985 31.042985 1.5707963 1.5707963 1.5707963 -4.9737992e-14 -3.126388e-13 -7.8159701e-14 300 6.6453262 85.853767 -6238.3552 29915.098 1.0000033 7123.716 5048.8505 6542.4992 -136.00562 750.73195 -1436.6619 -11588.844 31.042985 31.042985 31.042985 1.5707963 1.5707963 1.5707963 4.1300297e-13 9.343637e-13 -4.3698378e-13 400 5.4986371 79.882299 -5550.9936 29915.098 1.0000033 5653.3442 4311.7633 6687.8733 -106.69487 1035.5172 -1805.6706 -10661.78 31.042985 31.042985 31.042985 1.5707963 1.5707963 1.5707963 -1.0036416e-13 -9.8676622e-13 3.3395509e-13 500 6.9493562 72.743741 -5311.2012 29915.098 1.0000033 5886.4729 4619.8394 5427.2912 -163.61603 770.58431 -1219.6041 -10345.829 31.042985 31.042985 31.042985 1.5707963 1.5707963 1.5707963 -5.4356519e-13 -5.3290705e-13 9.9475983e-14 600 6.1778537 74.270177 -4761.594 29915.098 1.0000033 5566.946 3172.7348 5545.1012 128.55252 826.70497 -1818.7999 -9627.1479 31.042985 31.042985 31.042985 1.5707963 1.5707963 1.5707963 -6.750156e-14 -2.8421709e-13 3.2152059e-13 700 6.4735221 71.258904 -5040.2467 29915.098 1.0000033 5530.8044 3312.3734 6277.5622 16.855569 603.73943 -1671.3659 -10014.937 31.042985 31.042985 31.042985 1.5707963 1.5707963 1.5707963 5.1514348e-14 1.1919354e-12 1.7852386e-13 800 7.4157532 72.717356 -5398.7814 29915.098 1.0000033 6530.7179 3006.9672 6658.6591 -362.17538 769.29866 -1694.9167 -10490.57 31.042985 31.042985 31.042985 1.5707963 1.5707963 1.5707963 -1.9095836e-13 1.5685231e-12 -2.4158453e-13 900 5.6250538 73.051649 -4662.0987 29915.098 1.0000033 5152.642 2347.9954 6485.6588 -346.37391 851.92202 -1934.9193 -9532.5987 31.042985 31.042985 31.042985 1.5707963 1.5707963 1.5707963 -2.1316282e-14 6.3238303e-13 -2.1138646e-13 1000 6.6453723 77.262896 -4619.5939 29915.098 1.0000033 5477.9185 1987.0636 6393.7996 -140.08004 850.39955 -1943.8382 -9456.2776 31.042985 31.042985 31.042985 1.5707963 1.5707963 1.5707963 1.2114754e-12 6.7856831e-13 8.8462571e-13 Loop time of 1.35411 on 32 procs for 1000 steps with 3000 atoms 99.2% CPU use with 32 MPI tasks x 1 OpenMP threads Minimization stats: Stopping criterion = max iterations Energy initial, next-to-last, final = 15356.3019669914 -8707.07024271778 -9456.27761502613 Force two-norm initial, final = 5447.0541 77.262896 Force max component initial, final = 1022.784 6.6453723 Final line search alpha, max atom move = 3.8290963e-05 0.0002544577 Iterations, force evaluations = 1000 1140 MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 0.74789 | 0.92513 | 1.098 | 8.9 | 68.32 Bond | 0.0051257 | 0.0058276 | 0.0074834 | 0.6 | 0.43 Neigh | 0.0077114 | 0.0079988 | 0.0082629 | 0.1 | 0.59 Comm | 0.19329 | 0.36 | 0.54238 | 14.1 | 26.59 Output | 0.0019852 | 0.0019901 | 0.0020607 | 0.0 | 0.15 Modify | 0 | 0 | 0 | 0.0 | 0.00 Other | | 0.05317 | | | 3.93 Nlocal: 93.75 ave 107 max 81 min Histogram: 3 3 0 8 5 3 3 0 3 4 Nghost: 3570.53 ave 3624 max 3513 min Histogram: 2 8 5 1 0 0 0 2 11 3 Neighs: 29820.3 ave 36426 max 23783 min Histogram: 3 1 1 8 6 3 5 2 0 3 Total # of neighbors = 954249 Ave neighs/atom = 318.083 Ave special neighs/atom = 2 Neighbor list builds = 7 Dangerous builds = 0 undump sci log 2.3_Velocities.out #------------------------------------------------------------------------------- # Stage 2.3: Set the initial velocities for $T #------------------------------------------------------------------------------- velocity all create 298.2 72489 dist gaussian mom yes rot no log 2.4_NPT.out #------------------------------------------------------------------------------- # Stage 2.4: NPT integration for 2 ns with a timestep of 1 fs # Temperature 298.2 K # Pressure 1 atm #------------------------------------------------------------------------------- # Fix shake is needed for water in pcff+ fix shaken all shake .0001 20 50000 a 1 b 1 Finding SHAKE clusters ... 0 = # of size 2 clusters 0 = # of size 3 clusters 0 = # of size 4 clusters 1000 = # of frozen angles find clusters CPU = 0.001 seconds reset_timestep 0 thermo_style custom step v_time press vol v_sysdensity temp ebond eangle edihed eimp evdwl ecoul etail elong pe ke thermo ${Nthermo} thermo 0 fix 1 movable npt temp 298.2 298.2 100 iso 1 1 100 drag 0 mtk yes nreset 40000 fix 2 movable ave/time 1 199999 200000 v_time c_thermo_temp c_thermo_press v_sysvol v_sysdensity v_etotal v_cella v_cellb v_cellc v_cellalpha v_cellbeta v_cellgamma v_pe v_ke v_evdwl v_coulomb v_sxx v_syy v_szz v_syz v_sxz v_sxy file 2.4_averages.txt off 1 fix 3 movable ave/time 200 1 200 v_time c_thermo_temp c_thermo_press v_sysvol v_sysdensity v_etotal v_cella v_cellb v_cellc v_cellalpha v_cellbeta v_cellgamma v_pe v_ke v_evdwl v_coulomb v_sxx v_syy v_szz v_syz v_sxz v_sxy file 2.4_instantaneous.txt restart 2000000 2.4.restart dump sci all custom 200000 2.4.xyz id mol type q xs ys zs timestep 1 run 2000000 generated 0 of 1 mixed pair_coeff terms from sixthpower mixing rule Setting up Verlet run ... Unit style : real Current step : 0 Time step : 1 SHAKE stats (type/ave/delta/count) on step 0 1 0.993009 0.058347 3000 1 98.5447 12.3003 1000 Per MPI rank memory allocation (min/avg/max) = 8.787 | 8.793 | 8.798 Mbytes Step v_time Press Volume v_sysdensity Temp E_bond E_angle E_dihed E_impro E_vdwl E_coul E_tail E_long PotEng KinEng 0 1e-06 -755.06453 29915.098 1.0000033 447.37459 0 0 0 0 2778.1061 -13273.508 -103.93236 0 -10495.402 2665.7438 SHAKE stats (type/ave/delta/count) on step 50000 1 0.969948 1.2518e-06 3000 1 103.7 9.29645e-05 1000 SHAKE stats (type/ave/delta/count) on step 100000 1 0.970113 1.26604e-06 3000 1 103.7 9.74779e-05 1000 SHAKE stats (type/ave/delta/count) on step 150000 1 0.969991 1.43983e-06 3000 1 103.7 0.000113031 1000 SHAKE stats (type/ave/delta/count) on step 200000 1 0.969885 1.522e-06 3000 1 103.7 0.000112396 1000 SHAKE stats (type/ave/delta/count) on step 250000 1 0.970039 1.6945e-06 3000 1 103.7 0.000105922 1000 SHAKE stats (type/ave/delta/count) on step 300000 1 0.969932 1.64484e-06 3000 1 103.7 0.000117213 1000 SHAKE stats (type/ave/delta/count) on step 350000 1 0.970065 1.18608e-06 3000 1 103.7 9.0553e-05 1000 SHAKE stats (type/ave/delta/count) on step 400000 1 0.970023 1.30631e-06 3000 1 103.7 9.7881e-05 1000 SHAKE stats (type/ave/delta/count) on step 450000 1 0.970234 2.59679e-06 3000 1 103.7 0.000167042 1000 SHAKE stats (type/ave/delta/count) on step 500000 1 0.970073 1.87086e-06 3000 1 103.7 0.000145164 1000 SHAKE stats (type/ave/delta/count) on step 550000 1 0.969999 1.14664e-06 3000 1 103.7 9.77849e-05 1000 SHAKE stats (type/ave/delta/count) on step 600000 1 0.969985 1.39233e-06 3000 1 103.7 0.000109542 1000 SHAKE stats (type/ave/delta/count) on step 650000 1 0.969971 1.12489e-06 3000 1 103.7 0.00010164 1000 SHAKE stats (type/ave/delta/count) on step 700000 1 0.96998 1.64048e-06 3000 1 103.7 0.000114895 1000 SHAKE stats (type/ave/delta/count) on step 750000 1 0.970111 1.54206e-06 3000 1 103.7 0.000120681 1000 SHAKE stats (type/ave/delta/count) on step 800000 1 0.970097 1.3116e-06 3000 1 103.7 0.000119828 1000 SHAKE stats (type/ave/delta/count) on step 850000 1 0.970031 1.54274e-06 3000 1 103.7 9.97522e-05 1000 SHAKE stats (type/ave/delta/count) on step 900000 1 0.969996 1.29209e-06 3000 1 103.7 9.41064e-05 1000 SHAKE stats (type/ave/delta/count) on step 950000 1 0.970054 1.51822e-06 3000 1 103.7 0.000155945 1000 SHAKE stats (type/ave/delta/count) on step 1000000 1 0.969957 1.47421e-06 3000 1 103.7 8.94468e-05 1000 SHAKE stats (type/ave/delta/count) on step 1050000 1 0.970113 1.45368e-06 3000 1 103.7 9.26296e-05 1000 SHAKE stats (type/ave/delta/count) on step 1100000 1 0.969825 1.72164e-06 3000 1 103.7 0.00014218 1000 SHAKE stats (type/ave/delta/count) on step 1150000 1 0.97009 2.26841e-06 3000 1 103.7 0.00013978 1000 SHAKE stats (type/ave/delta/count) on step 1200000 1 0.969967 1.23114e-06 3000 1 103.7 9.10834e-05 1000 SHAKE stats (type/ave/delta/count) on step 1250000 1 0.969878 1.28014e-06 3000 1 103.7 0.000106748 1000 SHAKE stats (type/ave/delta/count) on step 1300000 1 0.969951 1.20286e-06 3000 1 103.7 9.67506e-05 1000 SHAKE stats (type/ave/delta/count) on step 1350000 1 0.96999 1.1856e-06 3000 1 103.7 9.67865e-05 1000 SHAKE stats (type/ave/delta/count) on step 1400000 1 0.969908 1.47647e-06 3000 1 103.7 0.000138004 1000 SHAKE stats (type/ave/delta/count) on step 1450000 1 0.970069 1.86682e-06 3000 1 103.7 0.000112364 1000 SHAKE stats (type/ave/delta/count) on step 1500000 1 0.970054 1.90047e-06 3000 1 103.7 0.000133838 1000 SHAKE stats (type/ave/delta/count) on step 1550000 1 0.97003 1.54525e-06 3000 1 103.7 9.70636e-05 1000 SHAKE stats (type/ave/delta/count) on step 1600000 1 0.969943 1.84091e-06 3000 1 103.7 0.000129051 1000 SHAKE stats (type/ave/delta/count) on step 1650000 1 0.969972 1.56219e-06 3000 1 103.7 0.000141151 1000 SHAKE stats (type/ave/delta/count) on step 1700000 1 0.969974 1.29313e-06 3000 1 103.7 0.000102137 1000 SHAKE stats (type/ave/delta/count) on step 1750000 1 0.97 1.31759e-06 3000 1 103.7 9.67278e-05 1000 SHAKE stats (type/ave/delta/count) on step 1800000 1 0.970115 1.13215e-06 3000 1 103.7 9.58676e-05 1000 SHAKE stats (type/ave/delta/count) on step 1850000 1 0.969946 1.20226e-06 3000 1 103.7 9.94339e-05 1000 SHAKE stats (type/ave/delta/count) on step 1900000 1 0.970005 1.13906e-06 3000 1 103.7 9.31866e-05 1000 SHAKE stats (type/ave/delta/count) on step 1950000 1 0.969969 1.25437e-06 3000 1 103.7 9.48926e-05 1000 SHAKE stats (type/ave/delta/count) on step 2000000 1 0.970063 3.25047e-06 3000 1 103.7 0.000190044 1000 2000000 2000000 -4208.6691 26763.922 1.1177434 297.90029 0 0 0 0 2567.014 -16542.17 -116.16932 0 -13975.156 1775.0804 Loop time of 3059.55 on 32 procs for 2000000 steps with 3000 atoms Performance: 56.479 ns/day, 0.425 hours/ns, 653.692 timesteps/s 99.4% CPU use with 32 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 1662.7 | 1812.4 | 1984.5 | 180.5 | 59.24 Bond | 1.255 | 1.4843 | 2.2465 | 14.7 | 0.05 Neigh | 102.8 | 103.03 | 103.16 | 1.1 | 3.37 Comm | 431.53 | 588.54 | 733.14 | 293.8 | 19.24 Output | 0.0048943 | 0.006133 | 0.0073401 | 1.0 | 0.00 Modify | 517.31 | 532.8 | 551.5 | 49.9 | 17.41 Other | | 21.32 | | | 0.70 Nlocal: 93.75 ave 104 max 85 min Histogram: 1 0 5 9 4 5 4 2 1 1 Nghost: 3858.06 ave 3886 max 3830 min Histogram: 4 4 1 4 1 4 4 4 5 1 Neighs: 33255.7 ave 37161 max 28553 min Histogram: 1 1 2 7 1 7 3 5 4 1 Total # of neighbors = 1064182 Ave neighs/atom = 354.72733 Ave special neighs/atom = 2 Neighbor list builds = 89067 Dangerous builds = 0 undump sci restart 0 dump sci all custom 2000000 2.4.xyz id mol type q xs ys zs run 0 generated 0 of 1 mixed pair_coeff terms from sixthpower mixing rule Setting up Verlet run ... Unit style : real Current step : 2000000 Time step : 1 SHAKE stats (type/ave/delta/count) on step 2000000 1 0.970063 3.25047e-06 3000 1 103.7 0.000190044 1000 Per MPI rank memory allocation (min/avg/max) = 8.79 | 8.795 | 8.801 Mbytes Step v_time Press Volume v_sysdensity Temp E_bond E_angle E_dihed E_impro E_vdwl E_coul E_tail E_long PotEng KinEng 2000000 2000000 -3722.4463 26763.922 1.1177434 297.90029 0 0 0 0 2567.014 -16542.17 -116.16932 0 -13975.156 1775.0804 Loop time of 6.42872e-06 on 32 procs for 0 steps with 3000 atoms 109.4% CPU use with 32 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 0 | 0 | 0 | 0.0 | 0.00 Bond | 0 | 0 | 0 | 0.0 | 0.00 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0 | 0 | 0 | 0.0 | 0.00 Output | 0 | 0 | 0 | 0.0 | 0.00 Modify | 0 | 0 | 0 | 0.0 | 0.00 Other | | 6.429e-06 | | |100.00 Nlocal: 93.75 ave 104 max 85 min Histogram: 1 0 5 9 4 5 4 2 1 1 Nghost: 3858.06 ave 3886 max 3830 min Histogram: 4 4 1 4 1 4 4 4 5 1 Neighs: 33255.7 ave 37161 max 28553 min Histogram: 1 1 2 7 1 7 3 5 4 1 Total # of neighbors = 1064182 Ave neighs/atom = 354.72733 Ave special neighs/atom = 2 Neighbor list builds = 0 Dangerous builds = 0 undump sci unfix 1 unfix 2 unfix 3 unfix shaken log 2.5_NVT.out #------------------------------------------------------------------------------- # Stage 2.5: NVT integration for 2 ns with a timestep of 1 fs # Temperature 298.2 K #------------------------------------------------------------------------------- # Fix shake is needed for water in pcff+ fix shaken all shake .0001 20 50000 a 1 b 1 Finding SHAKE clusters ... 0 = # of size 2 clusters 0 = # of size 3 clusters 0 = # of size 4 clusters 1000 = # of frozen angles find clusters CPU = 0.001 seconds reset_timestep 0 thermo_style custom step v_time press vol v_sysdensity temp ebond eangle edihed eimp evdwl ecoul etail elong pe ke thermo ${Nthermo} thermo 0 fix 1 movable nvt temp 298.2 298.2 100 drag 0.0 fix 2 movable ave/time 1 200000 200000 v_time c_thermo_temp c_thermo_press v_sysvol v_sysdensity v_etotal v_pe v_ke v_evdwl v_coulomb v_sxx v_syy v_szz v_syz v_sxz v_sxy file 2.5_averages.txt off 1 fix 3 movable ave/time 200 1 200 v_time c_thermo_temp c_thermo_press v_sysvol v_sysdensity v_etotal v_pe v_ke v_evdwl v_coulomb v_sxx v_syy v_szz v_syz v_sxz v_sxy file 2.5_instantaneous.txt restart 2000000 2.5.restart dump sci all custom 200000 2.5.xyz id mol type q xs ys zs timestep 1 run 2000000 generated 0 of 1 mixed pair_coeff terms from sixthpower mixing rule Setting up Verlet run ... Unit style : real Current step : 0 Time step : 1 SHAKE stats (type/ave/delta/count) on step 0 1 0.97 2.10398e-10 3000 1 103.7 1.22977e-08 1000 Per MPI rank memory allocation (min/avg/max) = 8.79 | 8.795 | 8.801 Mbytes Step v_time Press Volume v_sysdensity Temp E_bond E_angle E_dihed E_impro E_vdwl E_coul E_tail E_long PotEng KinEng 0 1e-06 -3727.6426 26763.922 1.1177434 297.90029 0 0 0 0 2567.0096 -16552.361 -116.16932 0 -13985.351 1775.0804 SHAKE stats (type/ave/delta/count) on step 50000 1 0.97 1.33976e-06 3000 1 103.7 0.000109865 1000 SHAKE stats (type/ave/delta/count) on step 100000 1 0.97 1.58224e-06 3000 1 103.7 0.000111456 1000 SHAKE stats (type/ave/delta/count) on step 150000 1 0.97 1.18808e-06 3000 1 103.7 9.29505e-05 1000 SHAKE stats (type/ave/delta/count) on step 200000 1 0.97 1.13089e-06 3000 1 103.7 0.000100915 1000 SHAKE stats (type/ave/delta/count) on step 250000 1 0.97 1.3653e-06 3000 1 103.7 0.00015753 1000 SHAKE stats (type/ave/delta/count) on step 300000 1 0.97 1.74248e-06 3000 1 103.7 0.000113139 1000 SHAKE stats (type/ave/delta/count) on step 350000 1 0.97 1.336e-06 3000 1 103.7 0.000120271 1000 SHAKE stats (type/ave/delta/count) on step 400000 1 0.97 1.0765e-06 3000 1 103.7 9.69623e-05 1000 SHAKE stats (type/ave/delta/count) on step 450000 1 0.97 1.22561e-06 3000 1 103.7 8.87765e-05 1000 SHAKE stats (type/ave/delta/count) on step 500000 1 0.97 1.87196e-06 3000 1 103.7 0.00018829 1000 SHAKE stats (type/ave/delta/count) on step 550000 1 0.97 1.9046e-06 3000 1 103.7 0.00013532 1000 SHAKE stats (type/ave/delta/count) on step 600000 1 0.97 2.01355e-06 3000 1 103.7 0.000147222 1000 SHAKE stats (type/ave/delta/count) on step 650000 1 0.97 1.20712e-06 3000 1 103.7 0.000109142 1000 SHAKE stats (type/ave/delta/count) on step 700000 1 0.97 1.6872e-06 3000 1 103.7 0.000124612 1000 SHAKE stats (type/ave/delta/count) on step 750000 1 0.97 1.21637e-06 3000 1 103.7 8.99016e-05 1000 SHAKE stats (type/ave/delta/count) on step 800000 1 0.97 1.22755e-06 3000 1 103.7 0.000109335 1000 SHAKE stats (type/ave/delta/count) on step 850000 1 0.97 1.76037e-06 3000 1 103.7 0.0001126 1000 SHAKE stats (type/ave/delta/count) on step 900000 1 0.97 1.31264e-06 3000 1 103.7 0.000113133 1000 SHAKE stats (type/ave/delta/count) on step 950000 1 0.97 1.26568e-06 3000 1 103.7 0.00010345 1000 SHAKE stats (type/ave/delta/count) on step 1000000 1 0.97 1.17285e-06 3000 1 103.7 9.4665e-05 1000 SHAKE stats (type/ave/delta/count) on step 1050000 1 0.97 1.37583e-06 3000 1 103.7 0.000122114 1000 SHAKE stats (type/ave/delta/count) on step 1100000 1 0.97 1.3577e-06 3000 1 103.7 0.000126336 1000 SHAKE stats (type/ave/delta/count) on step 1150000 1 0.97 1.34461e-06 3000 1 103.7 0.000146663 1000 SHAKE stats (type/ave/delta/count) on step 1200000 1 0.97 1.59032e-06 3000 1 103.7 0.000128458 1000 SHAKE stats (type/ave/delta/count) on step 1250000 1 0.97 1.95638e-06 3000 1 103.7 0.000153848 1000 SHAKE stats (type/ave/delta/count) on step 1300000 1 0.97 1.01264e-06 3000 1 103.7 0.000107888 1000 SHAKE stats (type/ave/delta/count) on step 1350000 1 0.97 1.57361e-06 3000 1 103.7 0.000117202 1000 SHAKE stats (type/ave/delta/count) on step 1400000 1 0.97 1.12605e-06 3000 1 103.7 9.12038e-05 1000 SHAKE stats (type/ave/delta/count) on step 1450000 1 0.97 1.37627e-06 3000 1 103.7 0.000137695 1000 SHAKE stats (type/ave/delta/count) on step 1500000 1 0.97 1.45407e-06 3000 1 103.7 9.8738e-05 1000 SHAKE stats (type/ave/delta/count) on step 1550000 1 0.97 1.21104e-06 3000 1 103.7 0.000102444 1000 SHAKE stats (type/ave/delta/count) on step 1600000 1 0.97 1.59569e-06 3000 1 103.7 0.000128561 1000 SHAKE stats (type/ave/delta/count) on step 1650000 1 0.97 1.43513e-06 3000 1 103.7 0.000126453 1000 SHAKE stats (type/ave/delta/count) on step 1700000 1 0.97 1.63991e-06 3000 1 103.7 0.000122222 1000 SHAKE stats (type/ave/delta/count) on step 1750000 1 0.97 2.16811e-06 3000 1 103.7 0.00012987 1000 SHAKE stats (type/ave/delta/count) on step 1800000 1 0.97 1.47299e-06 3000 1 103.7 0.000108467 1000 SHAKE stats (type/ave/delta/count) on step 1850000 1 0.97 2.14574e-06 3000 1 103.7 0.000143598 1000 SHAKE stats (type/ave/delta/count) on step 1900000 1 0.97 1.263e-06 3000 1 103.7 9.80381e-05 1000 SHAKE stats (type/ave/delta/count) on step 1950000 1 0.97 1.23414e-06 3000 1 103.7 0.000101154 1000 SHAKE stats (type/ave/delta/count) on step 2000000 1 0.97 1.75612e-06 3000 1 103.7 0.000149725 1000 2000000 2000000 -445.09542 26763.922 1.1177434 293.40142 0 0 0 0 2712.8962 -13041.416 -116.16932 0 -10328.52 1748.2732 Loop time of 2999.26 on 32 procs for 2000000 steps with 3000 atoms Performance: 57.614 ns/day, 0.417 hours/ns, 666.832 timesteps/s 99.3% CPU use with 32 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 1626.4 | 1833.5 | 2021 | 232.1 | 61.13 Bond | 1.1703 | 1.3481 | 1.9964 | 12.9 | 0.04 Neigh | 101.94 | 102.16 | 102.3 | 1.0 | 3.41 Comm | 461.51 | 628.86 | 845.24 | 371.7 | 20.97 Output | 0.0038934 | 0.0051068 | 0.0063518 | 1.0 | 0.00 Modify | 392.58 | 413.67 | 444.9 | 92.1 | 13.79 Other | | 19.73 | | | 0.66 Nlocal: 93.75 ave 103 max 83 min Histogram: 1 2 3 4 5 4 1 3 4 5 Nghost: 3852.38 ave 3885 max 3823 min Histogram: 3 3 7 2 1 1 7 7 0 1 Neighs: 33271.1 ave 38972 max 29409 min Histogram: 5 4 4 3 4 5 2 1 2 2 Total # of neighbors = 1064675 Ave neighs/atom = 354.89167 Ave special neighs/atom = 2 Neighbor list builds = 87197 Dangerous builds = 0 undump sci restart 0 dump sci all custom 2000000 2.5.xyz id mol type q xs ys zs run 0 generated 0 of 1 mixed pair_coeff terms from sixthpower mixing rule Setting up Verlet run ... Unit style : real Current step : 2000000 Time step : 1 SHAKE stats (type/ave/delta/count) on step 2000000 1 0.97 1.75612e-06 3000 1 103.7 0.000149725 1000 Per MPI rank memory allocation (min/avg/max) = 8.79 | 8.796 | 8.804 Mbytes Step v_time Press Volume v_sysdensity Temp E_bond E_angle E_dihed E_impro E_vdwl E_coul E_tail E_long PotEng KinEng 2000000 2000000 -443.5811 26763.922 1.1177434 293.40142 0 0 0 0 2712.8962 -13041.416 -116.16932 0 -10328.52 1748.2732 Loop time of 6.88906e-06 on 32 procs for 0 steps with 3000 atoms 109.3% CPU use with 32 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 0 | 0 | 0 | 0.0 | 0.00 Bond | 0 | 0 | 0 | 0.0 | 0.00 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0 | 0 | 0 | 0.0 | 0.00 Output | 0 | 0 | 0 | 0.0 | 0.00 Modify | 0 | 0 | 0 | 0.0 | 0.00 Other | | 6.889e-06 | | |100.00 Nlocal: 93.75 ave 103 max 83 min Histogram: 1 2 3 4 5 4 1 3 4 5 Nghost: 3852.38 ave 3885 max 3823 min Histogram: 3 3 7 2 1 1 7 7 0 1 Neighs: 33271.1 ave 38972 max 29409 min Histogram: 5 4 4 3 4 5 2 1 2 2 Total # of neighbors = 1064675 Ave neighs/atom = 354.89167 Ave special neighs/atom = 2 Neighbor list builds = 0 Dangerous builds = 0 undump sci unfix 1 unfix 2 unfix 3 unfix shaken log 2.6_NVT.out #------------------------------------------------------------------------------- # Stage 2.6: NVT integration for 2 ns with a timestep of 1 fs # Temperature 298.2 K #------------------------------------------------------------------------------- # Fix shake is needed for water in pcff+ fix shaken all shake .0001 20 50000 a 1 b 1 Finding SHAKE clusters ... 0 = # of size 2 clusters 0 = # of size 3 clusters 0 = # of size 4 clusters 1000 = # of frozen angles find clusters CPU = 0.001 seconds reset_timestep 0 thermo_style custom step v_time press vol v_sysdensity temp ebond eangle edihed eimp evdwl ecoul etail elong pe ke thermo ${Nthermo} thermo 0 fix 1 movable nvt temp 298.2 298.2 100 drag 0.0 fix 2 movable ave/time 1 200000 200000 v_time c_thermo_temp c_thermo_press v_sysvol v_sysdensity v_etotal v_pe v_ke v_evdwl v_coulomb v_sxx v_syy v_szz v_syz v_sxz v_sxy v_surfacetension file 2.6_averages.txt off 1 fix 3 movable ave/time 200 1 200 v_time c_thermo_temp c_thermo_press v_sysvol v_sysdensity v_etotal v_pe v_ke v_evdwl v_coulomb v_sxx v_syy v_szz v_syz v_sxz v_sxy v_surfacetension file 2.6_instantaneous.txt restart 2000000 2.6.restart dump sci all custom 200000 2.6.xyz id mol type q xs ys zs timestep 1 run 2000000 generated 0 of 1 mixed pair_coeff terms from sixthpower mixing rule Setting up Verlet run ... Unit style : real Current step : 0 Time step : 1 SHAKE stats (type/ave/delta/count) on step 0 1 0.97 1.67122e-12 3000 1 103.7 1.44667e-10 1000 Per MPI rank memory allocation (min/avg/max) = 8.79 | 8.796 | 8.804 Mbytes Step v_time Press Volume v_sysdensity Temp E_bond E_angle E_dihed E_impro E_vdwl E_coul E_tail E_long PotEng KinEng 0 1e-06 -443.5505 26763.922 1.1177434 293.40142 0 0 0 0 2712.8961 -13041.401 -116.16932 0 -10328.505 1748.2732 SHAKE stats (type/ave/delta/count) on step 50000 1 0.97 1.07878e-06 3000 1 103.7 9.77554e-05 1000 SHAKE stats (type/ave/delta/count) on step 100000 1 0.97 1.71273e-06 3000 1 103.7 0.000133217 1000 SHAKE stats (type/ave/delta/count) on step 150000 1 0.97 1.42979e-06 3000 1 103.7 8.93603e-05 1000 SHAKE stats (type/ave/delta/count) on step 200000 1 0.97 1.2003e-06 3000 1 103.7 0.00010498 1000 SHAKE stats (type/ave/delta/count) on step 250000 1 0.97 1.45345e-06 3000 1 103.7 0.000102301 1000 SHAKE stats (type/ave/delta/count) on step 300000 1 0.97 1.97075e-06 3000 1 103.7 0.000111243 1000 SHAKE stats (type/ave/delta/count) on step 350000 1 0.97 1.3982e-06 3000 1 103.7 0.000113106 1000 SHAKE stats (type/ave/delta/count) on step 400000 1 0.97 1.48269e-06 3000 1 103.7 0.000119219 1000 SHAKE stats (type/ave/delta/count) on step 450000 1 0.97 1.53176e-06 3000 1 103.7 0.000103606 1000 SHAKE stats (type/ave/delta/count) on step 500000 1 0.97 1.53747e-06 3000 1 103.7 0.000102103 1000 SHAKE stats (type/ave/delta/count) on step 550000 1 0.97 1.46013e-06 3000 1 103.7 9.84499e-05 1000 SHAKE stats (type/ave/delta/count) on step 600000 1 0.97 1.21963e-06 3000 1 103.7 9.29744e-05 1000 SHAKE stats (type/ave/delta/count) on step 650000 1 0.97 1.32278e-06 3000 1 103.7 0.000113193 1000 SHAKE stats (type/ave/delta/count) on step 700000 1 0.97 1.44559e-06 3000 1 103.7 9.93277e-05 1000 SHAKE stats (type/ave/delta/count) on step 750000 1 0.97 1.1727e-06 3000 1 103.7 9.5025e-05 1000 SHAKE stats (type/ave/delta/count) on step 800000 1 0.97 1.41562e-06 3000 1 103.7 9.43937e-05 1000 SHAKE stats (type/ave/delta/count) on step 850000 1 0.97 1.24836e-06 3000 1 103.7 9.55229e-05 1000 SHAKE stats (type/ave/delta/count) on step 900000 1 0.97 1.01904e-06 3000 1 103.7 9.26796e-05 1000 SHAKE stats (type/ave/delta/count) on step 950000 1 0.97 1.39031e-06 3000 1 103.7 9.94295e-05 1000 SHAKE stats (type/ave/delta/count) on step 1000000 1 0.97 1.35109e-06 3000 1 103.7 9.86247e-05 1000 SHAKE stats (type/ave/delta/count) on step 1050000 1 0.97 2.73745e-06 3000 1 103.7 0.000153106 1000 SHAKE stats (type/ave/delta/count) on step 1100000 1 0.97 1.19604e-06 3000 1 103.7 0.0001034 1000 SHAKE stats (type/ave/delta/count) on step 1150000 1 0.97 1.56969e-06 3000 1 103.7 0.000109792 1000 SHAKE stats (type/ave/delta/count) on step 1200000 1 0.97 1.22839e-06 3000 1 103.7 0.000125893 1000 SHAKE stats (type/ave/delta/count) on step 1250000 1 0.97 1.21518e-06 3000 1 103.7 8.64665e-05 1000 SHAKE stats (type/ave/delta/count) on step 1300000 1 0.97 1.42037e-06 3000 1 103.7 0.000109974 1000 SHAKE stats (type/ave/delta/count) on step 1350000 1 0.97 1.155e-06 3000 1 103.7 9.07315e-05 1000 SHAKE stats (type/ave/delta/count) on step 1400000 1 0.97 1.30888e-06 3000 1 103.7 0.000103115 1000 SHAKE stats (type/ave/delta/count) on step 1450000 1 0.97 1.6989e-06 3000 1 103.7 0.000111774 1000 SHAKE stats (type/ave/delta/count) on step 1500000 1 0.97 1.09e-06 3000 1 103.7 0.000102603 1000 SHAKE stats (type/ave/delta/count) on step 1550000 1 0.97 1.41178e-06 3000 1 103.7 9.784e-05 1000 SHAKE stats (type/ave/delta/count) on step 1600000 1 0.97 1.31548e-06 3000 1 103.7 0.000106518 1000 SHAKE stats (type/ave/delta/count) on step 1650000 1 0.97 1.3859e-06 3000 1 103.7 0.000137791 1000 SHAKE stats (type/ave/delta/count) on step 1700000 1 0.97 1.49634e-06 3000 1 103.7 9.28148e-05 1000 SHAKE stats (type/ave/delta/count) on step 1750000 1 0.97 1.33172e-06 3000 1 103.7 0.000104593 1000 SHAKE stats (type/ave/delta/count) on step 1800000 1 0.97 1.64952e-06 3000 1 103.7 0.000116364 1000 SHAKE stats (type/ave/delta/count) on step 1850000 1 0.97 1.36319e-06 3000 1 103.7 9.11465e-05 1000 SHAKE stats (type/ave/delta/count) on step 1900000 1 0.97 1.32547e-06 3000 1 103.7 9.99401e-05 1000 SHAKE stats (type/ave/delta/count) on step 1950000 1 0.97 1.33513e-06 3000 1 103.7 0.000107758 1000 SHAKE stats (type/ave/delta/count) on step 2000000 1 0.97 1.16305e-06 3000 1 103.7 9.58328e-05 1000 2000000 2000000 -3165.0128 26763.922 1.1177434 298.16469 0 0 0 0 2815.19 -16726.646 -116.16932 0 -13911.456 1776.6558 Loop time of 3021.37 on 32 procs for 2000000 steps with 3000 atoms Performance: 57.193 ns/day, 0.420 hours/ns, 661.952 timesteps/s 99.3% CPU use with 32 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 1619.9 | 1833.3 | 2069.9 | 233.8 | 60.68 Bond | 1.1384 | 1.3646 | 1.9795 | 13.2 | 0.05 Neigh | 101.21 | 101.34 | 101.48 | 0.6 | 3.35 Comm | 407.68 | 624.65 | 844.3 | 386.1 | 20.67 Output | 0.0043531 | 0.0055861 | 0.0068386 | 1.0 | 0.00 Modify | 416.46 | 438.29 | 467.64 | 85.1 | 14.51 Other | | 22.45 | | | 0.74 Nlocal: 93.75 ave 106 max 83 min Histogram: 2 1 5 7 4 3 4 3 2 1 Nghost: 3864.62 ave 3882 max 3852 min Histogram: 5 2 2 5 6 4 4 1 1 2 Neighs: 33261.1 ave 38778 max 28036 min Histogram: 2 2 1 8 6 0 5 6 1 1 Total # of neighbors = 1064355 Ave neighs/atom = 354.785 Ave special neighs/atom = 2 Neighbor list builds = 86904 Dangerous builds = 0 undump sci restart 0 dump sci all custom 2000000 2.6.xyz id mol type q xs ys zs run 0 generated 0 of 1 mixed pair_coeff terms from sixthpower mixing rule Setting up Verlet run ... Unit style : real Current step : 2000000 Time step : 1 SHAKE stats (type/ave/delta/count) on step 2000000 1 0.97 1.16305e-06 3000 1 103.7 9.58328e-05 1000 Per MPI rank memory allocation (min/avg/max) = 8.79 | 8.796 | 8.804 Mbytes Step v_time Press Volume v_sysdensity Temp E_bond E_angle E_dihed E_impro E_vdwl E_coul E_tail E_long PotEng KinEng 2000000 2000000 -3163.9738 26763.922 1.1177434 298.16469 0 0 0 0 2815.19 -16726.646 -116.16932 0 -13911.456 1776.6558 Loop time of 1.34676e-05 on 32 procs for 0 steps with 3000 atoms 99.1% CPU use with 32 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 0 | 0 | 0 | 0.0 | 0.00 Bond | 0 | 0 | 0 | 0.0 | 0.00 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0 | 0 | 0 | 0.0 | 0.00 Output | 0 | 0 | 0 | 0.0 | 0.00 Modify | 0 | 0 | 0 | 0.0 | 0.00 Other | | 1.347e-05 | | |100.00 Nlocal: 93.75 ave 106 max 83 min Histogram: 2 1 5 7 4 3 4 3 2 1 Nghost: 3864.62 ave 3882 max 3852 min Histogram: 5 2 2 5 6 4 4 1 1 2 Neighs: 33261.1 ave 38778 max 28036 min Histogram: 2 2 1 8 6 0 5 6 1 1 Total # of neighbors = 1064355 Ave neighs/atom = 354.785 Ave special neighs/atom = 2 Neighbor list builds = 0 Dangerous builds = 0 undump sci unfix 1 unfix 2 unfix 3 unfix shaken Total wall time: 2:31:21