Enter the name of the KIM Model you wish to perform calculations for: Enter the name of the species you wish to simulate: Enter the lattice type of the crystal ('bcc', 'diamond', 'fcc', 'hcp', or 'sc'): Enter the lattice constant 'a' in meters: If the lattice type is hcp, enter the lattice constant 'c' in meters (if you are not using hcp as the lattice, simply put any value here, as it will be ignored): Parameters Input: W bcc Pair_Morse_Shifted_GirifalcoWeizer_LowCutoff_W__MO_489351836217_001 [3.18056891859] CELL_SIZE_MIN: 5 CELL_SIZE_MAX: 7 Smallest System Size: 250 Largest System Size: 686 [Calculation] Supercell Size: 5 Unrelaxed Cell: [[15.90284459 0. 0. ] [ 0. 15.90284459 0. ] [ 0. 0. 15.90284459]] Unrelaxed Cell Vector: [15.90284459295, 0.0, 15.90284459295, 0.0, 0.0, 15.90284459295] Unrelaxed Cell Energy: -2069.49025342 Energy of Unrelaxed Cell With Vacancy: -2069.49025342 ========== Loop: 0 Position Relaxation... Step Time Energy fmax FIRE: 0 22:03:44 -2052.934331 1.9847 FIRE: 1 22:03:44 -2053.242775 1.6176 FIRE: 2 22:03:44 -2053.666357 1.0343 FIRE: 3 22:03:44 -2054.053624 0.7323 FIRE: 4 22:03:44 -2054.376640 0.6515 FIRE: 5 22:03:44 -2054.630500 0.4620 FIRE: 6 22:03:44 -2054.782139 0.3714 FIRE: 7 22:03:44 -2054.844900 0.5461 FIRE: 8 22:03:44 -2054.864652 0.5073 FIRE: 9 22:03:44 -2054.900453 0.4334 FIRE: 10 22:03:44 -2054.945849 0.3322 FIRE: 11 22:03:44 -2054.993251 0.2702 FIRE: 12 22:03:44 -2055.035745 0.2156 FIRE: 13 22:03:44 -2055.068540 0.1446 FIRE: 14 22:03:44 -2055.089666 0.1399 FIRE: 15 22:03:44 -2055.100194 0.1414 FIRE: 16 22:03:44 -2055.099839 0.1724 FIRE: 17 22:03:44 -2055.100559 0.1708 FIRE: 18 22:03:44 -2055.101952 0.1677 FIRE: 19 22:03:44 -2055.103928 0.1631 FIRE: 20 22:03:44 -2055.106363 0.1568 FIRE: 21 22:03:45 -2055.109107 0.1490 FIRE: 22 22:03:45 -2055.111997 0.1395 FIRE: 23 22:03:45 -2055.114871 0.1284 FIRE: 24 22:03:45 -2055.117854 0.1141 FIRE: 25 22:03:45 -2055.120725 0.0961 FIRE: 26 22:03:45 -2055.123225 0.0736 FIRE: 27 22:03:45 -2055.125125 0.0466 FIRE: 28 22:03:45 -2055.126350 0.0349 FIRE: 29 22:03:45 -2055.127060 0.0428 FIRE: 30 22:03:45 -2055.127598 0.0504 FIRE: 31 22:03:45 -2055.128336 0.0628 FIRE: 32 22:03:45 -2055.129469 0.0717 FIRE: 33 22:03:45 -2055.130796 0.0655 FIRE: 34 22:03:45 -2055.131631 0.0429 FIRE: 35 22:03:45 -2055.131700 0.0416 FIRE: 36 22:03:45 -2055.131831 0.0390 FIRE: 37 22:03:45 -2055.132006 0.0353 FIRE: 38 22:03:45 -2055.132203 0.0306 FIRE: 39 22:03:45 -2055.132401 0.0252 FIRE: 40 22:03:45 -2055.132579 0.0193 FIRE: 41 22:03:45 -2055.132722 0.0135 FIRE: 42 22:03:45 -2055.132828 0.0107 FIRE: 43 22:03:45 -2055.132885 0.0083 FIRE: 44 22:03:45 -2055.132892 0.0111 FIRE: 45 22:03:45 -2055.132895 0.0110 FIRE: 46 22:03:45 -2055.132901 0.0108 FIRE: 47 22:03:45 -2055.132909 0.0104 FIRE: 48 22:03:45 -2055.132919 0.0099 FIRE: 49 22:03:45 -2055.132931 0.0094 WARNING: Max number of steps exceeded. Structure may be unstable. Relaxation Completed. Steps: 50 Cell Size Relaxation... Optimization terminated successfully. Current function value: 6.068910 Iterations: 363 Function evaluations: 673 Current VFE: 6.06891019743 Energy of Supercell: -2069.49025342 Unrelaxed Cell Volume: 4021.83681063 Current Relaxed Cell Volume: 4015.37946494 Current Relaxation Volume: 6.45734569256 Current Cell: [[ 1.58943293e+01 0.00000000e+00 0.00000000e+00] [ 1.00252780e-05 1.58943285e+01 0.00000000e+00] [ 8.48637362e-06 -1.01893244e-07 1.58943291e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax FIRE: 0 22:04:04 -2055.143382 0.0086 FIRE: 1 22:04:04 -2055.143405 0.0077 FIRE: 2 22:04:04 -2055.143436 0.0058 FIRE: 3 22:04:04 -2055.143462 0.0034 FIRE: 4 22:04:04 -2055.143473 0.0047 FIRE: 5 22:04:04 -2055.143475 0.0055 FIRE: 6 22:04:04 -2055.143477 0.0053 FIRE: 7 22:04:04 -2055.143479 0.0049 FIRE: 8 22:04:04 -2055.143482 0.0043 FIRE: 9 22:04:04 -2055.143486 0.0035 FIRE: 10 22:04:04 -2055.143490 0.0026 FIRE: 11 22:04:04 -2055.143492 0.0019 FIRE: 12 22:04:04 -2055.143494 0.0015 FIRE: 13 22:04:04 -2055.143495 0.0013 FIRE: 14 22:04:04 -2055.143496 0.0016 FIRE: 15 22:04:04 -2055.143496 0.0016 FIRE: 16 22:04:04 -2055.143496 0.0016 FIRE: 17 22:04:04 -2055.143496 0.0015 FIRE: 18 22:04:04 -2055.143496 0.0014 FIRE: 19 22:04:04 -2055.143496 0.0014 FIRE: 20 22:04:04 -2055.143497 0.0013 FIRE: 21 22:04:04 -2055.143497 0.0011 FIRE: 22 22:04:04 -2055.143497 0.0010 FIRE: 23 22:04:04 -2055.143498 0.0008 Relaxation Completed. Steps: 23 Cell Size Relaxation... Optimization terminated successfully. Current function value: 6.068795 Iterations: 258 Function evaluations: 499 Current VFE: 6.0687946752 Energy of Supercell: -2069.49025342 Unrelaxed Cell Volume: 4021.83681063 Current Relaxed Cell Volume: 4015.37778654 Current Relaxation Volume: 6.45902409067 Current Cell: [[ 1.58943249e+01 0.00000000e+00 0.00000000e+00] [-2.24645304e-07 1.58943266e+01 0.00000000e+00] [ 1.12478166e-06 -1.89547878e-07 1.58943288e+01]] ========== Loop: 2 Position Relaxation... Step Time Energy fmax FIRE: 0 22:04:17 -2055.143498 0.0008 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 6.068795 Iterations: 130 Function evaluations: 315 Step Time Energy fmax FIRE: 0 22:04:25 -2055.143498 0.0008 FIRE: 1 22:04:25 -2055.143498 0.0008 FIRE: 2 22:04:26 -2055.143498 0.0006 FIRE: 3 22:04:26 -2055.143499 0.0005 FIRE: 4 22:04:26 -2055.143499 0.0005 FIRE: 5 22:04:26 -2055.143499 0.0004 FIRE: 6 22:04:26 -2055.143499 0.0004 FIRE: 7 22:04:26 -2055.143500 0.0002 FIRE: 8 22:04:26 -2055.143500 0.0001 FIRE: 9 22:04:26 -2055.143500 0.0001 FIRE: 10 22:04:26 -2055.143500 0.0001 FIRE: 11 22:04:26 -2055.143500 0.0001 FIRE: 12 22:04:26 -2055.143500 0.0001 FIRE: 13 22:04:26 -2055.143500 0.0001 FIRE: 14 22:04:26 -2055.143500 0.0001 FIRE: 15 22:04:26 -2055.143500 0.0001 FIRE: 16 22:04:26 -2055.143500 0.0001 FIRE: 17 22:04:26 -2055.143500 0.0001 FIRE: 18 22:04:26 -2055.143500 0.0001 FIRE: 19 22:04:26 -2055.143500 0.0000 Optimization terminated successfully. Current function value: 6.068793 Iterations: 199 Function evaluations: 455 --------------- Calculation Completed. Number Of Atoms in Supercell: 250 Vacancy Formation Energy (relaxed): 6.0687927918 Vacancy Formation Energy (unrelaxed): 8.27796101368 Unrelaxed Cell Volume: 4021.83681063 Relaxed Cell Volume: 4015.37778654 Relaxation Volume: 6.45902409067 Relaxed Cell Vector: [15.894328068334316, -2.3048486112133357e-07, 15.894328449178339, 1.1681726561570755e-06, -1.9314891416666205e-07, 15.894326904992381] Unrelaxed Cell Vector: [15.90284459295, 0.0, 15.90284459295, 0.0, 0.0, 15.90284459295] Relaxed Cell: [[ 1.58943281e+01 0.00000000e+00 0.00000000e+00] [-2.30484861e-07 1.58943284e+01 0.00000000e+00] [ 1.16817266e-06 -1.93148914e-07 1.58943269e+01]] Unrelaxed Cell: [[15.90284459 0. 0. ] [ 0. 15.90284459 0. ] [ 0. 0. 15.90284459]] Supercell Size: 6 Unrelaxed Cell: [[19.08341351 0. 0. ] [ 0. 19.08341351 0. ] [ 0. 0. 19.08341351]] Unrelaxed Cell Vector: [19.08341351154, 0.0, 19.08341351154, 0.0, 0.0, 19.08341351154] Unrelaxed Cell Energy: -3576.07915791 Energy of Unrelaxed Cell With Vacancy: -3576.07915791 ========== Loop: 0 Position Relaxation... Step Time Energy fmax FIRE: 0 22:04:38 -3559.523236 1.9847 FIRE: 1 22:04:38 -3559.831678 1.6176 FIRE: 2 22:04:38 -3560.255253 1.0343 FIRE: 3 22:04:38 -3560.642515 0.7322 FIRE: 4 22:04:38 -3560.965559 0.6516 FIRE: 5 22:04:38 -3561.219425 0.4622 FIRE: 6 22:04:38 -3561.371050 0.3736 FIRE: 7 22:04:38 -3561.434128 0.5464 FIRE: 8 22:04:38 -3561.453937 0.5077 FIRE: 9 22:04:38 -3561.489912 0.4340 FIRE: 10 22:04:38 -3561.535751 0.3330 FIRE: 11 22:04:38 -3561.584165 0.2700 FIRE: 12 22:04:38 -3561.628741 0.2156 FIRE: 13 22:04:38 -3561.665381 0.1456 FIRE: 14 22:04:38 -3561.692807 0.1405 FIRE: 15 22:04:38 -3561.713418 0.1443 FIRE: 16 22:04:38 -3561.726829 0.1607 FIRE: 17 22:04:38 -3561.734143 0.2650 FIRE: 18 22:04:38 -3561.735401 0.2605 FIRE: 19 22:04:38 -3561.737829 0.2515 FIRE: 20 22:04:38 -3561.741255 0.2383 FIRE: 21 22:04:38 -3561.745447 0.2213 FIRE: 22 22:04:38 -3561.750132 0.2008 FIRE: 23 22:04:38 -3561.755022 0.1773 FIRE: 24 22:04:38 -3561.759844 0.1517 FIRE: 25 22:04:38 -3561.764820 0.1216 FIRE: 26 22:04:39 -3561.769623 0.0877 FIRE: 27 22:04:39 -3561.773923 0.0732 FIRE: 28 22:04:39 -3561.777521 0.0660 FIRE: 29 22:04:39 -3561.780509 0.0643 FIRE: 30 22:04:39 -3561.783293 0.0605 FIRE: 31 22:04:39 -3561.786401 0.0597 FIRE: 32 22:04:39 -3561.790105 0.0634 FIRE: 33 22:04:39 -3561.794007 0.0534 FIRE: 34 22:04:39 -3561.796935 0.0306 FIRE: 35 22:04:39 -3561.797607 0.0229 FIRE: 36 22:04:39 -3561.797727 0.0222 FIRE: 37 22:04:39 -3561.797945 0.0209 FIRE: 38 22:04:39 -3561.798222 0.0190 FIRE: 39 22:04:39 -3561.798512 0.0167 FIRE: 40 22:04:39 -3561.798770 0.0138 FIRE: 41 22:04:39 -3561.798965 0.0114 FIRE: 42 22:04:39 -3561.799084 0.0132 FIRE: 43 22:04:39 -3561.799136 0.0149 FIRE: 44 22:04:39 -3561.799129 0.0150 FIRE: 45 22:04:39 -3561.799134 0.0147 FIRE: 46 22:04:39 -3561.799144 0.0142 FIRE: 47 22:04:39 -3561.799158 0.0134 FIRE: 48 22:04:39 -3561.799174 0.0124 FIRE: 49 22:04:39 -3561.799193 0.0112 WARNING: Max number of steps exceeded. Structure may be unstable. Relaxation Completed. Steps: 50 Cell Size Relaxation... Optimization terminated successfully. Current function value: 5.995705 Iterations: 321 Function evaluations: 613 Current VFE: 5.99570457567 Energy of Supercell: -3576.07915791 Unrelaxed Cell Volume: 6949.73400877 Current Relaxed Cell Volume: 6943.1599879 Current Relaxation Volume: 6.57402086384 Current Cell: [[1.90773941e+01 0.00000000e+00 0.00000000e+00] [3.26106121e-05 1.90773934e+01 0.00000000e+00] [4.01790427e-05 9.21237172e-05 1.90773956e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax FIRE: 0 22:05:02 -3561.805492 0.0073 FIRE: 1 22:05:02 -3561.805512 0.0060 FIRE: 2 22:05:02 -3561.805536 0.0036 FIRE: 3 22:05:02 -3561.805548 0.0026 FIRE: 4 22:05:02 -3561.805549 0.0036 FIRE: 5 22:05:02 -3561.805550 0.0034 FIRE: 6 22:05:02 -3561.805553 0.0031 FIRE: 7 22:05:02 -3561.805557 0.0028 FIRE: 8 22:05:02 -3561.805561 0.0023 FIRE: 9 22:05:02 -3561.805564 0.0018 FIRE: 10 22:05:02 -3561.805567 0.0016 FIRE: 11 22:05:02 -3561.805568 0.0016 FIRE: 12 22:05:02 -3561.805569 0.0021 FIRE: 13 22:05:02 -3561.805571 0.0025 FIRE: 14 22:05:02 -3561.805573 0.0025 FIRE: 15 22:05:02 -3561.805577 0.0022 FIRE: 16 22:05:02 -3561.805580 0.0014 FIRE: 17 22:05:03 -3561.805582 0.0005 Relaxation Completed. Steps: 17 Cell Size Relaxation... Optimization terminated successfully. Current function value: 5.995615 Iterations: 322 Function evaluations: 577 Current VFE: 5.99561518668 Energy of Supercell: -3576.07915791 Unrelaxed Cell Volume: 6949.73400877 Current Relaxed Cell Volume: 6943.16510436 Current Relaxation Volume: 6.56890440171 Current Cell: [[ 1.90773981e+01 0.00000000e+00 0.00000000e+00] [ 2.61092597e-07 1.90773989e+01 0.00000000e+00] [-5.72645544e-07 3.20006775e-07 1.90774002e+01]] ========== Loop: 2 Position Relaxation... Step Time Energy fmax FIRE: 0 22:05:23 -3561.805582 0.0005 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 5.995615 Iterations: 145 Function evaluations: 316 Step Time Energy fmax FIRE: 0 22:05:33 -3561.805582 0.0005 FIRE: 1 22:05:33 -3561.805582 0.0004 FIRE: 2 22:05:33 -3561.805582 0.0003 FIRE: 3 22:05:33 -3561.805582 0.0003 FIRE: 4 22:05:33 -3561.805582 0.0003 FIRE: 5 22:05:33 -3561.805582 0.0003 FIRE: 6 22:05:33 -3561.805582 0.0003 FIRE: 7 22:05:33 -3561.805582 0.0002 FIRE: 8 22:05:33 -3561.805582 0.0002 FIRE: 9 22:05:33 -3561.805582 0.0002 FIRE: 10 22:05:33 -3561.805582 0.0002 FIRE: 11 22:05:34 -3561.805582 0.0002 FIRE: 12 22:05:34 -3561.805582 0.0002 FIRE: 13 22:05:34 -3561.805582 0.0002 FIRE: 14 22:05:34 -3561.805582 0.0001 FIRE: 15 22:05:34 -3561.805582 0.0001 FIRE: 16 22:05:34 -3561.805582 0.0001 FIRE: 17 22:05:34 -3561.805582 0.0001 FIRE: 18 22:05:34 -3561.805582 0.0001 FIRE: 19 22:05:34 -3561.805582 0.0001 Optimization terminated successfully. Current function value: 5.995614 Iterations: 186 Function evaluations: 436 --------------- Calculation Completed. Number Of Atoms in Supercell: 432 Vacancy Formation Energy (relaxed): 5.9956144686 Vacancy Formation Energy (unrelaxed): 8.27796101368 Unrelaxed Cell Volume: 6949.73400877 Relaxed Cell Volume: 6943.16510436 Relaxation Volume: 6.56890440171 Relaxed Cell Vector: [19.077397022340925, 2.5674115508306425e-07, 19.07740111596917, -5.85370917573937e-07, 3.271179177030677e-07, 19.077399944195356] Unrelaxed Cell Vector: [19.08341351154, 0.0, 19.08341351154, 0.0, 0.0, 19.08341351154] Relaxed Cell: [[ 1.90773970e+01 0.00000000e+00 0.00000000e+00] [ 2.56741155e-07 1.90774011e+01 0.00000000e+00] [-5.85370918e-07 3.27117918e-07 1.90773999e+01]] Unrelaxed Cell: [[19.08341351 0. 0. ] [ 0. 19.08341351 0. ] [ 0. 0. 19.08341351]] Supercell Size: 7 Unrelaxed Cell: [[22.26398243 0. 0. ] [ 0. 22.26398243 0. ] [ 0. 0. 22.26398243]] Unrelaxed Cell Vector: [22.26398243013, 0.0, 22.26398243013, 0.0, 0.0, 22.26398243013] Unrelaxed Cell Energy: -5678.6812554 Energy of Unrelaxed Cell With Vacancy: -5678.6812554 ========== Loop: 0 Position Relaxation... Step Time Energy fmax FIRE: 0 22:05:52 -5662.125333 1.9847 FIRE: 1 22:05:52 -5662.433776 1.6176 FIRE: 2 22:05:52 -5662.857350 1.0343 FIRE: 3 22:05:52 -5663.244610 0.7322 FIRE: 4 22:05:52 -5663.567644 0.6516 FIRE: 5 22:05:52 -5663.821495 0.4622 FIRE: 6 22:05:52 -5663.973108 0.3736 FIRE: 7 22:05:52 -5664.036162 0.5465 FIRE: 8 22:05:52 -5664.055972 0.5077 FIRE: 9 22:05:52 -5664.091950 0.4340 FIRE: 10 22:05:52 -5664.137796 0.3331 FIRE: 11 22:05:52 -5664.186222 0.2700 FIRE: 12 22:05:52 -5664.230815 0.2156 FIRE: 13 22:05:52 -5664.267504 0.1455 FIRE: 14 22:05:52 -5664.295101 0.1407 FIRE: 15 22:05:52 -5664.316317 0.1447 FIRE: 16 22:05:52 -5664.331469 0.1601 FIRE: 17 22:05:52 -5664.342758 0.2646 FIRE: 18 22:05:53 -5664.354554 0.3317 FIRE: 19 22:05:53 -5664.371549 0.3357 FIRE: 20 22:05:53 -5664.393896 0.2609 FIRE: 21 22:05:53 -5664.413999 0.1158 FIRE: 22 22:05:53 -5664.422717 0.1065 FIRE: 23 22:05:53 -5664.424440 0.0980 FIRE: 24 22:05:53 -5664.427466 0.0820 FIRE: 25 22:05:53 -5664.431111 0.0604 FIRE: 26 22:05:53 -5664.434662 0.0359 FIRE: 27 22:05:53 -5664.437617 0.0478 FIRE: 28 22:05:53 -5664.439795 0.0636 FIRE: 29 22:05:53 -5664.441250 0.0687 FIRE: 30 22:05:53 -5664.442144 0.0608 FIRE: 31 22:05:53 -5664.442360 0.0391 FIRE: 32 22:05:53 -5664.442424 0.0381 FIRE: 33 22:05:53 -5664.442547 0.0361 FIRE: 34 22:05:53 -5664.442717 0.0332 FIRE: 35 22:05:53 -5664.442918 0.0295 FIRE: 36 22:05:53 -5664.443132 0.0250 FIRE: 37 22:05:53 -5664.443340 0.0199 FIRE: 38 22:05:53 -5664.443523 0.0143 FIRE: 39 22:05:53 -5664.443678 0.0082 FIRE: 40 22:05:53 -5664.443779 0.0066 FIRE: 41 22:05:53 -5664.443803 0.0065 FIRE: 42 22:05:53 -5664.443804 0.0065 FIRE: 43 22:05:53 -5664.443806 0.0063 FIRE: 44 22:05:53 -5664.443809 0.0061 FIRE: 45 22:05:54 -5664.443813 0.0058 FIRE: 46 22:05:54 -5664.443818 0.0054 FIRE: 47 22:05:54 -5664.443823 0.0050 FIRE: 48 22:05:54 -5664.443829 0.0045 FIRE: 49 22:05:54 -5664.443836 0.0043 WARNING: Max number of steps exceeded. Structure may be unstable. Relaxation Completed. Steps: 50 Cell Size Relaxation... Optimization terminated successfully. Current function value: 5.955434 Iterations: 286 Function evaluations: 560 Current VFE: 5.95543410269 Energy of Supercell: -5678.6812554 Unrelaxed Cell Volume: 11035.9202084 Current Relaxed Cell Volume: 11029.3002807 Current Relaxation Volume: 6.61992769906 Current Cell: [[2.22595329e+01 0.00000000e+00 0.00000000e+00] [4.58022090e-05 2.22595268e+01 0.00000000e+00] [4.65361366e-05 7.36014833e-05 2.22595297e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax FIRE: 0 22:06:31 -5664.447860 0.0040 FIRE: 1 22:06:31 -5664.447870 0.0037 FIRE: 2 22:06:31 -5664.447886 0.0032 FIRE: 3 22:06:31 -5664.447904 0.0025 FIRE: 4 22:06:32 -5664.447921 0.0019 FIRE: 5 22:06:32 -5664.447936 0.0017 FIRE: 6 22:06:32 -5664.447949 0.0011 FIRE: 7 22:06:32 -5664.447957 0.0006 Relaxation Completed. Steps: 7 Cell Size Relaxation... Optimization terminated successfully. Current function value: 5.955338 Iterations: 245 Function evaluations: 480 Current VFE: 5.9553375401 Energy of Supercell: -5678.6812554 Unrelaxed Cell Volume: 11035.9202084 Current Relaxed Cell Volume: 11029.2990946 Current Relaxation Volume: 6.62111373307 Current Cell: [[ 2.22595288e+01 0.00000000e+00 0.00000000e+00] [ 1.15418639e-06 2.22595298e+01 0.00000000e+00] [ 1.82994707e-07 -7.94588305e-07 2.22595285e+01]] ========== Loop: 2 Position Relaxation... Step Time Energy fmax FIRE: 0 22:07:00 -5664.447957 0.0006 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 5.955338 Iterations: 143 Function evaluations: 322 Step Time Energy fmax FIRE: 0 22:07:17 -5664.447957 0.0006 FIRE: 1 22:07:17 -5664.447958 0.0006 FIRE: 2 22:07:17 -5664.447959 0.0005 FIRE: 3 22:07:17 -5664.447960 0.0005 FIRE: 4 22:07:17 -5664.447961 0.0005 FIRE: 5 22:07:17 -5664.447962 0.0005 FIRE: 6 22:07:17 -5664.447963 0.0004 FIRE: 7 22:07:17 -5664.447963 0.0005 FIRE: 8 22:07:18 -5664.447964 0.0004 FIRE: 9 22:07:18 -5664.447964 0.0003 FIRE: 10 22:07:18 -5664.447964 0.0003 FIRE: 11 22:07:18 -5664.447964 0.0003 FIRE: 12 22:07:18 -5664.447964 0.0002 FIRE: 13 22:07:18 -5664.447964 0.0002 FIRE: 14 22:07:18 -5664.447964 0.0001 FIRE: 15 22:07:18 -5664.447964 0.0001 FIRE: 16 22:07:18 -5664.447964 0.0001 FIRE: 17 22:07:18 -5664.447964 0.0001 FIRE: 18 22:07:18 -5664.447964 0.0001 FIRE: 19 22:07:18 -5664.447964 0.0001 Optimization terminated successfully. Current function value: 5.955330 Iterations: 188 Function evaluations: 434 --------------- Calculation Completed. Number Of Atoms in Supercell: 686 Vacancy Formation Energy (relaxed): 5.95533002751 Vacancy Formation Energy (unrelaxed): 8.27796101372 Unrelaxed Cell Volume: 11035.9202084 Relaxed Cell Volume: 11029.2990946 Relaxation Volume: 6.62111373307 Relaxed Cell Vector: [22.25952873442111, 1.117636276148319e-06, 22.25952929971172, 1.8838292884441215e-07, -8.157587361344482e-07, 22.25952725508955] Unrelaxed Cell Vector: [22.26398243013, 0.0, 22.26398243013, 0.0, 0.0, 22.26398243013] Relaxed Cell: [[ 2.22595287e+01 0.00000000e+00 0.00000000e+00] [ 1.11763628e-06 2.22595293e+01 0.00000000e+00] [ 1.88382929e-07 -8.15758736e-07 2.22595273e+01]] Unrelaxed Cell: [[22.26398243 0. 0. ] [ 0. 22.26398243 0. ] [ 0. 0. 22.26398243]] [Calculation Results Summary] Sizes: [5, 6, 7] Unrelaxed Formation Energy By Size: [8.27796101368358, 8.27796101368449, 8.277961013724962] Formation Energy By Size: [6.0687927918047535, 5.995614468597978, 5.9553300275056245] Relaxation Volume By Size: [6.459024090670027, 6.5689044017126434, 6.621113733068341] [Extrapolation] Fitting with 2 points, including orders [0, 3] Fit with data beginning 0 Fit with Size: [5 6] Orders: [0, 3] Matrix A (Ax = y): [[1. 0.008 ] [1. 0.00462963]] Data for Fitting: [8.27796101 8.27796101] Fitting Results: (array([ 8.27796101e+00, -2.69821242e-10]), array([], dtype=float64), 2, array([1.41424176, 0.00238316])) Fit with Size: [5 6] Orders: [0, 3] Matrix A (Ax = y): [[1. 0.008 ] [1. 0.00462963]] Data for Fitting: [6.06879279 5.99561447] Fitting Results: (array([ 5.89509479, 21.71224974]), array([], dtype=float64), 2, array([1.41424176, 0.00238316])) Fit with Size: [5 6] Orders: [0, 3] Matrix A (Ax = y): [[1. 0.008 ] [1. 0.00462963]] Data for Fitting: [6.45902409 6.5689044 ] Fitting Results: (array([ 6.71983889, -32.60185053]), array([], dtype=float64), 2, array([1.41424176, 0.00238316])) Fit with data beginning 1 Fit with Size: [6 7] Orders: [0, 3] Matrix A (Ax = y): [[1. 0.00462963] [1. 0.00291545]] Data for Fitting: [8.27796101 8.27796101] Fitting Results: (array([ 8.27796101e+00, -2.36108628e-08]), array([], dtype=float64), 2, array([1.41422363e+00, 1.21209807e-03])) Fit with Size: [6 7] Orders: [0, 3] Matrix A (Ax = y): [[1. 0.00462963] [1. 0.00291545]] Data for Fitting: [5.99561447 5.95533003] Fitting Results: (array([ 5.88681476, 23.50073757]), array([], dtype=float64), 2, array([1.41422363e+00, 1.21209807e-03])) Fit with Size: [6 7] Orders: [0, 3] Matrix A (Ax = y): [[1. 0.00462963] [1. 0.00291545]] Data for Fitting: [6.5689044 6.62111373] Fitting Results: (array([ 6.70991071, -30.45736174]), array([], dtype=float64), 2, array([1.41422363e+00, 1.21209807e-03])) Fitting with 3 points, including orders [0, 3] Fit with data beginning 0 Fit with Size: [5 6 7] Orders: [0, 3] Matrix A (Ax = y): [[1. 0.008 ] [1. 0.00462963] [1. 0.00291545]] Data for Fitting: [8.27796101 8.27796101 8.27796101] Fitting Results: (array([ 8.27796101e+00, -7.04428309e-09]), array([4.52926081e-22]), 2, array([1.73207406, 0.00365829])) Fit with Size: [5 6 7] Orders: [0, 3] Matrix A (Ax = y): [[1. 0.008 ] [1. 0.00462963] [1. 0.00291545]] Data for Fitting: [6.06879279 5.99561447 5.95533003] Fitting Results: (array([ 5.89138287, 22.23138336]), array([2.65918603e-06]), 2, array([1.73207406, 0.00365829])) Fit with Size: [5 6 7] Orders: [0, 3] Matrix A (Ax = y): [[1. 0.008 ] [1. 0.00462963] [1. 0.00291545]] Data for Fitting: [6.45902409 6.5689044 6.62111373] Fitting Results: (array([ 6.71538811, -31.97938265]), array([3.82317609e-06]), 2, array([1.73207406, 0.00365829])) Fitting with 3 points, including orders [0, 3, 4] Fit with data beginning 0 Fit with Size: [5 6 7] Orders: [0, 3, 4] Matrix A (Ax = y): [[1.00000000e+00 8.00000000e-03 1.60000000e-03] [1.00000000e+00 4.62962963e-03 7.71604938e-04] [1.00000000e+00 2.91545190e-03 4.16493128e-04]] Data for Fitting: [8.27796101 8.27796101 8.27796101] Fitting Results: (array([ 8.27796101e+00, -1.48796374e-07, 6.04291132e-07]), array([], dtype=float64), 3, array([1.73207481e+00, 3.75760407e-03, 3.42874264e-05])) Fit with Size: [5 6 7] Orders: [0, 3, 4] Matrix A (Ax = y): [[1.00000000e+00 8.00000000e-03 1.60000000e-03] [1.00000000e+00 4.62962963e-03 7.71604938e-04] [1.00000000e+00 2.91545190e-03 4.16493128e-04]] Data for Fitting: [6.06879279 5.99561447 5.95533003] Fitting Results: (array([ 5.87813407, 33.09289247, -46.30276399]), array([], dtype=float64), 3, array([1.73207481e+00, 3.75760407e-03, 3.42874264e-05])) Fit with Size: [5 6 7] Orders: [0, 3, 4] Matrix A (Ax = y): [[1.00000000e+00 8.00000000e-03 1.60000000e-03] [1.00000000e+00 4.62962963e-03 7.71604938e-04] [1.00000000e+00 2.91545190e-03 4.16493128e-04]] Data for Fitting: [6.45902409 6.5689044 6.62111373] Fitting Results: (array([ 6.69950212, -18.95587526, -55.51939268]), array([], dtype=float64), 3, array([1.73207481e+00, 3.75760407e-03, 3.42874264e-05])) Fitting with 3 points, including orders [0, 3, 5] Fit with data beginning 0 Fit with Size: [5 6 7] Orders: [0, 3, 5] Matrix A (Ax = y): [[1.00000000e+00 8.00000000e-03 3.20000000e-04] [1.00000000e+00 4.62962963e-03 1.28600823e-04] [1.00000000e+00 2.91545190e-03 5.94990183e-05]] Data for Fitting: [8.27796101 8.27796101 8.27796101] Fitting Results: (array([ 8.27796101e+00, -8.07168653e-08, 1.41661335e-06]), array([], dtype=float64), 3, array([1.73207409e+00, 3.66323826e-03, 1.50029206e-05])) Fit with Size: [5 6 7] Orders: [0, 3, 5] Matrix A (Ax = y): [[1.00000000e+00 8.00000000e-03 3.20000000e-04] [1.00000000e+00 4.62962963e-03 1.28600823e-04] [1.00000000e+00 2.91545190e-03 5.94990183e-05]] Data for Fitting: [6.06879279 5.99561447 5.95533003] Fitting Results: (array([ 5.88051603, 27.87641769, -108.54555044]), array([], dtype=float64), 3, array([1.73207409e+00, 3.66323826e-03, 1.50029206e-05])) Fit with Size: [5 6 7] Orders: [0, 3, 5] Matrix A (Ax = y): [[1.00000000e+00 8.00000000e-03 3.20000000e-04] [1.00000000e+00 4.62962963e-03 1.28600823e-04] [1.00000000e+00 2.91545190e-03 5.94990183e-05]] Data for Fitting: [6.45902409 6.5689044 6.62111373] Fitting Results: (array([ 6.7023582 , -25.21069649, -130.15169115]), array([], dtype=float64), 3, array([1.73207409e+00, 3.66323826e-03, 1.50029206e-05])) Fitting with 3 points, including orders [0, 3, 6] Fit with data beginning 0 Fit with Size: [5 6 7] Orders: [0, 3, 6] Matrix A (Ax = y): [[1.00000000e+00 8.00000000e-03 6.40000000e-05] [1.00000000e+00 4.62962963e-03 2.14334705e-05] [1.00000000e+00 2.91545190e-03 8.49985975e-06]] Data for Fitting: [8.27796101 8.27796101 8.27796101] Fitting Results: (array([ 8.27796101e+00, -5.82471826e-08, 4.59063679e-06]), array([], dtype=float64), 3, array([1.73207406e+00, 3.65852192e-03, 4.63568288e-06])) Fit with Size: [5 6 7] Orders: [0, 3, 6] Matrix A (Ax = y): [[1.00000000e+00 8.00000000e-03 6.40000000e-05] [1.00000000e+00 4.62962963e-03 2.14334705e-05] [1.00000000e+00 2.91545190e-03 8.49985975e-06]] Data for Fitting: [6.06879279 5.99561447 5.95533003] Fitting Results: (array([ 5.88206703, 26.15471708, -351.74961323]), array([], dtype=float64), 3, array([1.73207406e+00, 3.65852192e-03, 4.63568288e-06])) Fit with Size: [5 6 7] Orders: [0, 3, 6] Matrix A (Ax = y): [[1.00000000e+00 8.00000000e-03 6.40000000e-05] [1.00000000e+00 4.62962963e-03 2.14334705e-05] [1.00000000e+00 2.91545190e-03 8.49985975e-06]] Data for Fitting: [6.45902409 6.5689044 6.62111373] Fitting Results: (array([ 6.70421794, -27.27510397, -421.76585626]), array([], dtype=float64), 3, array([1.73207406e+00, 3.65852192e-03, 4.63568288e-06])) [Fitting Results Summary] Sizes: [5, 6, 7] Data Points Used: [2, 3, 3, 3, 3] Orders Included: [list([0, 3]) list([0, 3]) list([0, 3, 4]) list([0, 3, 5]) list([0, 3, 6])] Unrelaxed Formation Energy Fits By Size: [list([8.27796101368574, 8.277961013793798]) list([8.277961013734183]) list([8.277961013907086]) list([8.277961013876]) list([8.277961013855762])] Formation Energy Fits By Size: [list([5.895094793863398, 5.886814757616268]) list([5.891382873719454]) list([5.878134074453954]) list([5.880516026405586]) list([5.882067030410347])] Relaxation Volume Fits By Size: [list([6.7198388949030505, 6.709910706082758]) list([6.715388111962616]) list([6.69950212102922]) list([6.702358203784277]) list([6.704217937263782])] [Final Results] [ { "property-id" "tag:staff@noreply.openkim.org,2015-07-28:property/monovacancy-neutral-unrelaxed-formation-potential-energy-crystal-npt" "instance-id" 1 "unrelaxed-formation-potential-energy" { "source-value" 8.277961013793798 "source-unit" "eV" "source-std-uncert-value" 7.5125944940184065e-06 } "host-cauchy-stress" { "source-value" [ 0 0 0 0 0 0 ] "source-unit" "GPa" } "host-removed-atom" { "source-value" 1 } "host-short-name" { "source-value" [ "bcc" ] } "host-a" { "source-value" 3.18056891859 "source-unit" "angstrom" } "host-b" { "source-value" 3.18056891859 "source-unit" "angstrom" } "host-c" { "source-value" 3.18056891859 "source-unit" "angstrom" } "host-alpha" { "source-value" 90.0 "source-unit" "degree" } "host-beta" { "source-value" 90.0 "source-unit" "degree" } "host-gamma" { "source-value" 90.0 "source-unit" "degree" } "host-space-group" { "source-value" "Im-3m" } "host-wyckoff-multiplicity-and-letter" { "source-value" [ "2a" ] } "host-wyckoff-coordinates" { "source-value" [ [ 0.0 0.0 0.0 ] ] } "host-wyckoff-species" { "source-value" [ "W" ] } "reservoir-cohesive-potential-energy" { "source-value" 8.277961013694282 "source-unit" "eV" } "reservoir-short-name" { "source-value" [ "bcc" ] } "reservoir-cauchy-stress" { "source-value" [ 0.0 0.0 0.0 0.0 0.0 0.0 ] "source-unit" "GPa" } "reservoir-a" { "source-value" 3.18056891859 "source-unit" "angstrom" } "reservoir-b" { "source-value" 3.18056891859 "source-unit" "angstrom" } "reservoir-c" { "source-value" 3.18056891859 "source-unit" "angstrom" } "reservoir-alpha" { "source-value" 90.0 "source-unit" "degree" } "reservoir-beta" { "source-value" 90.0 "source-unit" "degree" } "reservoir-gamma" { "source-value" 90.0 "source-unit" "degree" } "reservoir-space-group" { "source-value" "Im-3m" } "reservoir-wyckoff-multiplicity-and-letter" { "source-value" [ "2a" ] } "reservoir-wyckoff-coordinates" { "source-value" [ [ 0.0 0.0 0.0 ] ] } "reservoir-wyckoff-species" { "source-value" [ "W" ] } } { "property-id" "tag:staff@noreply.openkim.org,2015-07-28:property/monovacancy-neutral-relaxed-formation-potential-energy-crystal-npt" "instance-id" 2 "relaxed-formation-potential-energy" { "source-value" 5.886814757616268 "source-unit" "eV" "source-std-uncert-value" 0.008680686413155789 } "host-cauchy-stress" { "source-value" [ 0 0 0 0 0 0 ] "source-unit" "GPa" } "host-removed-atom" { "source-value" 1 } "host-short-name" { "source-value" [ "bcc" ] } "host-a" { "source-value" 3.18056891859 "source-unit" "angstrom" } "host-b" { "source-value" 3.18056891859 "source-unit" "angstrom" } "host-c" { "source-value" 3.18056891859 "source-unit" "angstrom" } "host-alpha" { "source-value" 90.0 "source-unit" "degree" } "host-beta" { "source-value" 90.0 "source-unit" "degree" } "host-gamma" { "source-value" 90.0 "source-unit" "degree" } "host-space-group" { "source-value" "Im-3m" } "host-wyckoff-multiplicity-and-letter" { "source-value" [ "2a" ] } "host-wyckoff-coordinates" { "source-value" [ [ 0.0 0.0 0.0 ] ] } "host-wyckoff-species" { "source-value" [ "W" ] } "reservoir-cohesive-potential-energy" { "source-value" 8.277961013694282 "source-unit" "eV" } "reservoir-short-name" { "source-value" [ "bcc" ] } "reservoir-cauchy-stress" { "source-value" [ 0.0 0.0 0.0 0.0 0.0 0.0 ] "source-unit" "GPa" } "reservoir-a" { "source-value" 3.18056891859 "source-unit" "angstrom" } "reservoir-b" { "source-value" 3.18056891859 "source-unit" "angstrom" } "reservoir-c" { "source-value" 3.18056891859 "source-unit" "angstrom" } "reservoir-alpha" { "source-value" 90.0 "source-unit" "degree" } "reservoir-beta" { "source-value" 90.0 "source-unit" "degree" } "reservoir-gamma" { "source-value" 90.0 "source-unit" "degree" } "reservoir-space-group" { "source-value" "Im-3m" } "reservoir-wyckoff-multiplicity-and-letter" { "source-value" [ "2a" ] } "reservoir-wyckoff-coordinates" { "source-value" [ [ 0.0 0.0 0.0 ] ] } "reservoir-wyckoff-species" { "source-value" [ "W" ] } } { "property-id" "tag:staff@noreply.openkim.org,2015-07-28:property/monovacancy-neutral-relaxation-volume-crystal-npt" "instance-id" 3 "relaxation-volume" { "source-value" 6.709910706082758 "source-unit" "angstrom^3" "source-std-uncert-value" 0.010447744217773667 } "host-cauchy-stress" { "source-value" [ 0 0 0 0 0 0 ] "source-unit" "GPa" } "host-removed-atom" { "source-value" 1 } "host-short-name" { "source-value" [ "bcc" ] } "host-a" { "source-value" 3.18056891859 "source-unit" "angstrom" } "host-b" { "source-value" 3.18056891859 "source-unit" "angstrom" } "host-c" { "source-value" 3.18056891859 "source-unit" "angstrom" } "host-alpha" { "source-value" 90.0 "source-unit" "degree" } "host-beta" { "source-value" 90.0 "source-unit" "degree" } "host-gamma" { "source-value" 90.0 "source-unit" "degree" } "host-space-group" { "source-value" "Im-3m" } "host-wyckoff-multiplicity-and-letter" { "source-value" [ "2a" ] } "host-wyckoff-coordinates" { "source-value" [ [ 0.0 0.0 0.0 ] ] } "host-wyckoff-species" { "source-value" [ "W" ] } } ]