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_HighCutoff_W__MO_646516726498_001 [3.1535765305199996] CELL_SIZE_MIN: 5 CELL_SIZE_MAX: 7 Smallest System Size: 250 Largest System Size: 686 [Calculation] Supercell Size: 5 Unrelaxed Cell: [[15.76788265 0. 0. ] [ 0. 15.76788265 0. ] [ 0. 0. 15.76788265]] Unrelaxed Cell Vector: [15.767882652599997, 0.0, 15.767882652599997, 0.0, 0.0, 15.767882652599997] Unrelaxed Cell Energy: -2205.95212843 Energy of Unrelaxed Cell With Vacancy: -2205.95212843 ========== Loop: 0 Position Relaxation... Step Time Energy fmax FIRE: 0 22:03:41 -2188.304511 2.2635 FIRE: 1 22:03:41 -2188.693471 1.8056 FIRE: 2 22:03:41 -2189.208643 1.0925 FIRE: 3 22:03:41 -2189.666814 0.8695 FIRE: 4 22:03:41 -2190.047399 0.7117 FIRE: 5 22:03:41 -2190.310226 0.4451 FIRE: 6 22:03:41 -2190.431557 0.5137 FIRE: 7 22:03:41 -2190.458235 0.4634 FIRE: 8 22:03:41 -2190.505710 0.3997 FIRE: 9 22:03:41 -2190.564143 0.3374 FIRE: 10 22:03:41 -2190.622861 0.2545 FIRE: 11 22:03:41 -2190.673338 0.1851 FIRE: 12 22:03:41 -2190.711105 0.1712 FIRE: 13 22:03:41 -2190.735712 0.1673 FIRE: 14 22:03:41 -2190.750203 0.2226 FIRE: 15 22:03:41 -2190.755934 0.3349 FIRE: 16 22:03:41 -2190.757387 0.3295 FIRE: 17 22:03:41 -2190.760190 0.3187 FIRE: 18 22:03:41 -2190.764149 0.3029 FIRE: 19 22:03:41 -2190.768996 0.2822 FIRE: 20 22:03:41 -2190.774413 0.2573 FIRE: 21 22:03:41 -2190.780061 0.2287 FIRE: 22 22:03:41 -2190.785617 0.1970 FIRE: 23 22:03:41 -2190.791312 0.1595 FIRE: 24 22:03:41 -2190.796737 0.1166 FIRE: 25 22:03:41 -2190.801460 0.0958 FIRE: 26 22:03:41 -2190.805203 0.0851 FIRE: 27 22:03:41 -2190.808033 0.0710 FIRE: 28 22:03:41 -2190.810383 0.0618 FIRE: 29 22:03:41 -2190.812842 0.0775 FIRE: 30 22:03:42 -2190.815779 0.0832 FIRE: 31 22:03:42 -2190.818964 0.0716 FIRE: 32 22:03:42 -2190.821408 0.0488 FIRE: 33 22:03:42 -2190.821723 0.0407 FIRE: 34 22:03:42 -2190.821850 0.0395 FIRE: 35 22:03:42 -2190.822082 0.0372 FIRE: 36 22:03:42 -2190.822379 0.0339 FIRE: 37 22:03:42 -2190.822693 0.0297 FIRE: 38 22:03:42 -2190.822978 0.0248 FIRE: 39 22:03:42 -2190.823205 0.0210 FIRE: 40 22:03:42 -2190.823362 0.0222 FIRE: 41 22:03:42 -2190.823472 0.0220 FIRE: 42 22:03:42 -2190.823548 0.0195 FIRE: 43 22:03:42 -2190.823615 0.0144 FIRE: 44 22:03:42 -2190.823693 0.0132 FIRE: 45 22:03:42 -2190.823781 0.0137 FIRE: 46 22:03:42 -2190.823858 0.0130 FIRE: 47 22:03:42 -2190.823901 0.0202 FIRE: 48 22:03:42 -2190.823913 0.0194 FIRE: 49 22:03:42 -2190.823935 0.0178 WARNING: Max number of steps exceeded. Structure may be unstable. Relaxation Completed. Steps: 50 Cell Size Relaxation... Optimization terminated successfully. Current function value: 6.295093 Iterations: 448 Function evaluations: 792 Current VFE: 6.29509326036 Energy of Supercell: -2205.95212843 Unrelaxed Cell Volume: 3920.30753729 Current Relaxed Cell Volume: 3914.38067032 Current Relaxation Volume: 5.92686697217 Current Cell: [[1.57599318e+01 0.00000000e+00 0.00000000e+00] [1.35925871e-07 1.57599332e+01 0.00000000e+00] [2.46988617e-07 4.92437968e-07 1.57599325e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax FIRE: 0 22:04:17 -2190.833227 0.0200 FIRE: 1 22:04:18 -2190.833273 0.0155 FIRE: 2 22:04:18 -2190.833329 0.0079 FIRE: 3 22:04:18 -2190.833355 0.0040 FIRE: 4 22:04:18 -2190.833350 0.0065 FIRE: 5 22:04:18 -2190.833353 0.0060 FIRE: 6 22:04:18 -2190.833357 0.0050 FIRE: 7 22:04:18 -2190.833362 0.0036 FIRE: 8 22:04:18 -2190.833366 0.0020 FIRE: 9 22:04:18 -2190.833368 0.0011 FIRE: 10 22:04:18 -2190.833368 0.0014 FIRE: 11 22:04:18 -2190.833368 0.0014 FIRE: 12 22:04:18 -2190.833368 0.0013 FIRE: 13 22:04:18 -2190.833369 0.0012 FIRE: 14 22:04:18 -2190.833369 0.0011 FIRE: 15 22:04:18 -2190.833369 0.0010 Relaxation Completed. Steps: 15 Cell Size Relaxation... Optimization terminated successfully. Current function value: 6.294951 Iterations: 155 Function evaluations: 343 Current VFE: 6.29495088363 Energy of Supercell: -2205.95212843 Unrelaxed Cell Volume: 3920.30753729 Current Relaxed Cell Volume: 3914.36735511 Current Relaxation Volume: 5.94018218567 Current Cell: [[1.57599144e+01 0.00000000e+00 0.00000000e+00] [1.38876607e-07 1.57599154e+01 0.00000000e+00] [2.44815303e-07 5.02465659e-07 1.57599140e+01]] ========== Loop: 2 Position Relaxation... Step Time Energy fmax FIRE: 0 22:04:34 -2190.833369 0.0010 FIRE: 1 22:04:34 -2190.833369 0.0007 Relaxation Completed. Steps: 1 Cell Size Relaxation... Optimization terminated successfully. Current function value: 6.294950 Iterations: 117 Function evaluations: 297 Current VFE: 6.29495044277 Energy of Supercell: -2205.95212843 Unrelaxed Cell Volume: 3920.30753729 Current Relaxed Cell Volume: 3914.36717264 Current Relaxation Volume: 5.94036465435 Current Cell: [[1.57599149e+01 0.00000000e+00 0.00000000e+00] [1.40612704e-07 1.57599143e+01 0.00000000e+00] [2.47875106e-07 5.08746700e-07 1.57599139e+01]] ========== Loop: 3 Position Relaxation... Step Time Energy fmax FIRE: 0 22:04:48 -2190.833369 0.0007 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 6.294950 Iterations: 115 Function evaluations: 279 Step Time Energy fmax FIRE: 0 22:05:00 -2190.833369 0.0007 FIRE: 1 22:05:00 -2190.833370 0.0005 FIRE: 2 22:05:00 -2190.833370 0.0004 FIRE: 3 22:05:00 -2190.833370 0.0005 FIRE: 4 22:05:00 -2190.833370 0.0005 FIRE: 5 22:05:00 -2190.833371 0.0004 FIRE: 6 22:05:00 -2190.833371 0.0002 FIRE: 7 22:05:00 -2190.833371 0.0004 FIRE: 8 22:05:00 -2190.833371 0.0003 FIRE: 9 22:05:00 -2190.833371 0.0003 FIRE: 10 22:05:00 -2190.833371 0.0002 FIRE: 11 22:05:00 -2190.833371 0.0001 FIRE: 12 22:05:00 -2190.833371 0.0001 FIRE: 13 22:05:00 -2190.833371 0.0001 FIRE: 14 22:05:00 -2190.833371 0.0001 FIRE: 15 22:05:00 -2190.833371 0.0001 FIRE: 16 22:05:00 -2190.833371 0.0001 FIRE: 17 22:05:00 -2190.833371 0.0001 FIRE: 18 22:05:00 -2190.833371 0.0001 FIRE: 19 22:05:00 -2190.833371 0.0001 Optimization terminated successfully. Current function value: 6.294949 Iterations: 189 Function evaluations: 436 --------------- Calculation Completed. Number Of Atoms in Supercell: 250 Vacancy Formation Energy (relaxed): 6.29494898153 Vacancy Formation Energy (unrelaxed): 8.82380851372 Unrelaxed Cell Volume: 3920.30753729 Relaxed Cell Volume: 3914.36717264 Relaxation Volume: 5.94036465435 Relaxed Cell Vector: [15.759915600642827, 1.396486549863576e-07, 15.759911204338824, 2.5518100420196005e-07, 5.137841023615399e-07, 15.759913413225663] Unrelaxed Cell Vector: [15.767882652599997, 0.0, 15.767882652599997, 0.0, 0.0, 15.767882652599997] Relaxed Cell: [[1.57599156e+01 0.00000000e+00 0.00000000e+00] [1.39648655e-07 1.57599112e+01 0.00000000e+00] [2.55181004e-07 5.13784102e-07 1.57599134e+01]] Unrelaxed Cell: [[15.76788265 0. 0. ] [ 0. 15.76788265 0. ] [ 0. 0. 15.76788265]] Supercell Size: 6 Unrelaxed Cell: [[18.92145918 0. 0. ] [ 0. 18.92145918 0. ] [ 0. 0. 18.92145918]] Unrelaxed Cell Vector: [18.92145918312, 0.0, 18.92145918312, 0.0, 0.0, 18.92145918312] Unrelaxed Cell Energy: -3811.88527793 Energy of Unrelaxed Cell With Vacancy: -3811.88527793 ========== Loop: 0 Position Relaxation... Step Time Energy fmax FIRE: 0 22:05:21 -3794.237661 2.2635 FIRE: 1 22:05:21 -3794.626623 1.8056 FIRE: 2 22:05:21 -3795.141804 1.0925 FIRE: 3 22:05:21 -3795.600041 0.8693 FIRE: 4 22:05:21 -3795.980765 0.7119 FIRE: 5 22:05:21 -3796.243737 0.4452 FIRE: 6 22:05:21 -3796.365356 0.5132 FIRE: 7 22:05:21 -3796.392034 0.4630 FIRE: 8 22:05:21 -3796.439495 0.3994 FIRE: 9 22:05:21 -3796.497909 0.3368 FIRE: 10 22:05:21 -3796.556706 0.2539 FIRE: 11 22:05:22 -3796.607766 0.1848 FIRE: 12 22:05:22 -3796.647385 0.1711 FIRE: 13 22:05:22 -3796.676175 0.1679 FIRE: 14 22:05:22 -3796.699053 0.2196 FIRE: 15 22:05:22 -3796.718478 0.3327 FIRE: 16 22:05:22 -3796.738709 0.4012 FIRE: 17 22:05:22 -3796.763978 0.3979 FIRE: 18 22:05:22 -3796.793352 0.3049 FIRE: 19 22:05:22 -3796.817167 0.1298 FIRE: 20 22:05:22 -3796.825591 0.1254 FIRE: 21 22:05:22 -3796.827755 0.1157 FIRE: 22 22:05:22 -3796.831524 0.0972 FIRE: 23 22:05:22 -3796.835978 0.0719 FIRE: 24 22:05:22 -3796.840148 0.0425 FIRE: 25 22:05:22 -3796.843338 0.0462 FIRE: 26 22:05:22 -3796.845317 0.0679 FIRE: 27 22:05:22 -3796.846264 0.0778 FIRE: 28 22:05:22 -3796.846504 0.0736 FIRE: 29 22:05:22 -3796.846588 0.0719 FIRE: 30 22:05:22 -3796.846748 0.0685 FIRE: 31 22:05:23 -3796.846971 0.0636 FIRE: 32 22:05:23 -3796.847236 0.0573 FIRE: 33 22:05:23 -3796.847520 0.0497 FIRE: 34 22:05:23 -3796.847800 0.0411 FIRE: 35 22:05:23 -3796.848053 0.0317 FIRE: 36 22:05:23 -3796.848281 0.0208 FIRE: 37 22:05:23 -3796.848450 0.0121 FIRE: 38 22:05:23 -3796.848532 0.0097 FIRE: 39 22:05:23 -3796.848519 0.0145 FIRE: 40 22:05:23 -3796.848523 0.0143 FIRE: 41 22:05:23 -3796.848529 0.0139 FIRE: 42 22:05:23 -3796.848538 0.0132 FIRE: 43 22:05:23 -3796.848550 0.0123 FIRE: 44 22:05:23 -3796.848563 0.0113 FIRE: 45 22:05:23 -3796.848577 0.0101 FIRE: 46 22:05:23 -3796.848592 0.0087 FIRE: 47 22:05:23 -3796.848608 0.0071 FIRE: 48 22:05:23 -3796.848625 0.0052 FIRE: 49 22:05:23 -3796.848642 0.0042 WARNING: Max number of steps exceeded. Structure may be unstable. Relaxation Completed. Steps: 50 Cell Size Relaxation... Optimization terminated successfully. Current function value: 6.207190 Iterations: 313 Function evaluations: 592 Current VFE: 6.20718976778 Energy of Supercell: -3811.88527793 Unrelaxed Cell Volume: 6774.29142444 Current Relaxed Cell Volume: 6768.22963991 Current Relaxation Volume: 6.06178453575 Current Cell: [[ 1.89158126e+01 0.00000000e+00 0.00000000e+00] [-7.32925640e-07 1.89158136e+01 0.00000000e+00] [ 2.97373910e-07 5.72974740e-07 1.89158150e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax FIRE: 0 22:06:11 -3796.854280 0.0038 FIRE: 1 22:06:12 -3796.854294 0.0031 FIRE: 2 22:06:12 -3796.854315 0.0027 FIRE: 3 22:06:12 -3796.854337 0.0026 FIRE: 4 22:06:12 -3796.854357 0.0022 FIRE: 5 22:06:12 -3796.854371 0.0016 FIRE: 6 22:06:12 -3796.854376 0.0022 FIRE: 7 22:06:12 -3796.854376 0.0020 FIRE: 8 22:06:12 -3796.854377 0.0015 FIRE: 9 22:06:12 -3796.854378 0.0009 Relaxation Completed. Steps: 9 Cell Size Relaxation... Optimization terminated successfully. Current function value: 6.207092 Iterations: 145 Function evaluations: 330 Current VFE: 6.20709175782 Energy of Supercell: -3811.88527793 Unrelaxed Cell Volume: 6774.29142444 Current Relaxed Cell Volume: 6768.22640042 Current Relaxation Volume: 6.06502402367 Current Cell: [[ 1.89158123e+01 0.00000000e+00 0.00000000e+00] [-7.39242072e-07 1.89158100e+01 0.00000000e+00] [ 3.04078282e-07 5.70805226e-07 1.89158098e+01]] ========== Loop: 2 Position Relaxation... Step Time Energy fmax FIRE: 0 22:06:36 -3796.854378 0.0009 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 6.207092 Iterations: 112 Function evaluations: 276 Step Time Energy fmax FIRE: 0 22:06:58 -3796.854378 0.0009 FIRE: 1 22:06:58 -3796.854378 0.0008 FIRE: 2 22:06:58 -3796.854379 0.0008 FIRE: 3 22:06:58 -3796.854380 0.0007 FIRE: 4 22:06:58 -3796.854381 0.0003 FIRE: 5 22:06:58 -3796.854381 0.0003 FIRE: 6 22:06:58 -3796.854381 0.0006 FIRE: 7 22:06:58 -3796.854381 0.0006 FIRE: 8 22:06:58 -3796.854381 0.0005 FIRE: 9 22:06:58 -3796.854381 0.0005 FIRE: 10 22:06:59 -3796.854382 0.0004 FIRE: 11 22:06:59 -3796.854382 0.0003 FIRE: 12 22:06:59 -3796.854382 0.0002 FIRE: 13 22:06:59 -3796.854382 0.0002 FIRE: 14 22:06:59 -3796.854382 0.0001 FIRE: 15 22:06:59 -3796.854382 0.0001 FIRE: 16 22:06:59 -3796.854382 0.0001 FIRE: 17 22:06:59 -3796.854382 0.0001 FIRE: 18 22:06:59 -3796.854382 0.0001 FIRE: 19 22:06:59 -3796.854382 0.0001 Optimization terminated successfully. Current function value: 6.207088 Iterations: 163 Function evaluations: 415 --------------- Calculation Completed. Number Of Atoms in Supercell: 432 Vacancy Formation Energy (relaxed): 6.20708759238 Vacancy Formation Energy (unrelaxed): 8.82380851373 Unrelaxed Cell Volume: 6774.29142444 Relaxed Cell Volume: 6768.22640042 Relaxation Volume: 6.06502402367 Relaxed Cell Vector: [18.91581231272172, -7.328615569442108e-07, 18.915810607761813, 3.020357438100712e-07, 5.915082552392e-07, 18.915812093602984] Unrelaxed Cell Vector: [18.92145918312, 0.0, 18.92145918312, 0.0, 0.0, 18.92145918312] Relaxed Cell: [[ 1.89158123e+01 0.00000000e+00 0.00000000e+00] [-7.32861557e-07 1.89158106e+01 0.00000000e+00] [ 3.02035744e-07 5.91508255e-07 1.89158121e+01]] Unrelaxed Cell: [[18.92145918 0. 0. ] [ 0. 18.92145918 0. ] [ 0. 0. 18.92145918]] Supercell Size: 7 Unrelaxed Cell: [[22.07503571 0. 0. ] [ 0. 22.07503571 0. ] [ 0. 0. 22.07503571]] Unrelaxed Cell Vector: [22.07503571364, 0.0, 22.07503571364, 0.0, 0.0, 22.07503571364] Unrelaxed Cell Energy: -6053.13264042 Energy of Unrelaxed Cell With Vacancy: -6053.13264042 ========== Loop: 0 Position Relaxation... Step Time Energy fmax FIRE: 0 22:07:33 -6035.485023 2.2635 FIRE: 1 22:07:33 -6035.873985 1.8056 FIRE: 2 22:07:33 -6036.389165 1.0925 FIRE: 3 22:07:33 -6036.847395 0.8693 FIRE: 4 22:07:33 -6037.228110 0.7119 FIRE: 5 22:07:33 -6037.491092 0.4453 FIRE: 6 22:07:33 -6037.612738 0.5133 FIRE: 7 22:07:33 -6037.639420 0.4631 FIRE: 8 22:07:33 -6037.686889 0.3994 FIRE: 9 22:07:33 -6037.745319 0.3368 FIRE: 10 22:07:33 -6037.804138 0.2539 FIRE: 11 22:07:33 -6037.855216 0.1851 FIRE: 12 22:07:33 -6037.894828 0.1713 FIRE: 13 22:07:33 -6037.923591 0.1679 FIRE: 14 22:07:33 -6037.946531 0.2198 FIRE: 15 22:07:33 -6037.966623 0.3329 FIRE: 16 22:07:33 -6037.989256 0.4021 FIRE: 17 22:07:34 -6038.020379 0.4006 FIRE: 18 22:07:34 -6038.060631 0.3108 FIRE: 19 22:07:34 -6038.099789 0.1384 FIRE: 20 22:07:34 -6038.123094 0.1233 FIRE: 21 22:07:34 -6038.126948 0.1877 FIRE: 22 22:07:34 -6038.129450 0.1700 FIRE: 23 22:07:34 -6038.133700 0.1367 FIRE: 24 22:07:34 -6038.138473 0.0912 FIRE: 25 22:07:34 -6038.142521 0.0390 FIRE: 26 22:07:34 -6038.145024 0.0305 FIRE: 27 22:07:34 -6038.145888 0.0545 FIRE: 28 22:07:34 -6038.145951 0.0531 FIRE: 29 22:07:34 -6038.146071 0.0502 FIRE: 30 22:07:34 -6038.146236 0.0461 FIRE: 31 22:07:34 -6038.146432 0.0408 FIRE: 32 22:07:34 -6038.146643 0.0346 FIRE: 33 22:07:34 -6038.146851 0.0276 FIRE: 34 22:07:34 -6038.147047 0.0201 FIRE: 35 22:07:35 -6038.147239 0.0136 FIRE: 36 22:07:35 -6038.147423 0.0112 FIRE: 37 22:07:35 -6038.147602 0.0146 FIRE: 38 22:07:35 -6038.147789 0.0159 FIRE: 39 22:07:35 -6038.147999 0.0168 FIRE: 40 22:07:35 -6038.148234 0.0182 FIRE: 41 22:07:35 -6038.148469 0.0158 FIRE: 42 22:07:35 -6038.148660 0.0098 FIRE: 43 22:07:35 -6038.148781 0.0160 FIRE: 44 22:07:35 -6038.148850 0.0167 FIRE: 45 22:07:35 -6038.148878 0.0101 FIRE: 46 22:07:35 -6038.148886 0.0095 FIRE: 47 22:07:35 -6038.148901 0.0085 FIRE: 48 22:07:35 -6038.148920 0.0070 FIRE: 49 22:07:35 -6038.148940 0.0052 WARNING: Max number of steps exceeded. Structure may be unstable. Relaxation Completed. Steps: 50 Cell Size Relaxation... Optimization terminated successfully. Current function value: 6.156257 Iterations: 384 Function evaluations: 701 Current VFE: 6.15625671514 Energy of Supercell: -6053.13264042 Unrelaxed Cell Volume: 10757.3238823 Current Relaxed Cell Volume: 10751.2039953 Current Relaxation Volume: 6.11988703962 Current Cell: [[ 2.20708507e+01 0.00000000e+00 0.00000000e+00] [-8.46101964e-06 2.20708519e+01 0.00000000e+00] [ 8.91872736e-06 -4.34256521e-06 2.20708435e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax FIRE: 0 22:09:16 -6038.152575 0.0039 FIRE: 1 22:09:16 -6038.152591 0.0030 FIRE: 2 22:09:16 -6038.152615 0.0024 FIRE: 3 22:09:16 -6038.152639 0.0023 FIRE: 4 22:09:16 -6038.152656 0.0021 FIRE: 5 22:09:16 -6038.152665 0.0018 FIRE: 6 22:09:16 -6038.152666 0.0021 FIRE: 7 22:09:16 -6038.152667 0.0020 FIRE: 8 22:09:16 -6038.152668 0.0019 FIRE: 9 22:09:16 -6038.152669 0.0017 FIRE: 10 22:09:16 -6038.152671 0.0014 FIRE: 11 22:09:16 -6038.152673 0.0013 FIRE: 12 22:09:16 -6038.152675 0.0013 FIRE: 13 22:09:16 -6038.152677 0.0012 FIRE: 14 22:09:16 -6038.152679 0.0011 FIRE: 15 22:09:16 -6038.152680 0.0008 Relaxation Completed. Steps: 15 Cell Size Relaxation... Optimization terminated successfully. Current function value: 6.156152 Iterations: 140 Function evaluations: 324 Current VFE: 6.15615183437 Energy of Supercell: -6053.13264042 Unrelaxed Cell Volume: 10757.3238823 Current Relaxed Cell Volume: 10751.2048172 Current Relaxation Volume: 6.11906514535 Current Cell: [[ 2.20708487e+01 0.00000000e+00 0.00000000e+00] [-8.18892900e-06 2.20708473e+01 0.00000000e+00] [ 9.08954995e-06 -4.48371928e-06 2.20708519e+01]] ========== Loop: 2 Position Relaxation... Step Time Energy fmax FIRE: 0 22:10:02 -6038.152680 0.0008 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 6.156152 Iterations: 117 Function evaluations: 283 Step Time Energy fmax FIRE: 0 22:10:44 -6038.152680 0.0008 FIRE: 1 22:10:44 -6038.152680 0.0007 FIRE: 2 22:10:44 -6038.152681 0.0006 FIRE: 3 22:10:44 -6038.152681 0.0006 FIRE: 4 22:10:44 -6038.152682 0.0006 FIRE: 5 22:10:44 -6038.152683 0.0004 FIRE: 6 22:10:44 -6038.152684 0.0003 FIRE: 7 22:10:44 -6038.152684 0.0003 FIRE: 8 22:10:44 -6038.152685 0.0002 FIRE: 9 22:10:44 -6038.152685 0.0001 FIRE: 10 22:10:45 -6038.152685 0.0001 FIRE: 11 22:10:45 -6038.152685 0.0001 FIRE: 12 22:10:45 -6038.152685 0.0001 FIRE: 13 22:10:45 -6038.152685 0.0001 FIRE: 14 22:10:45 -6038.152685 0.0001 FIRE: 15 22:10:45 -6038.152685 0.0001 FIRE: 16 22:10:45 -6038.152685 0.0001 FIRE: 17 22:10:45 -6038.152685 0.0001 FIRE: 18 22:10:45 -6038.152685 0.0001 FIRE: 19 22:10:45 -6038.152685 0.0001 Optimization terminated successfully. Current function value: 6.156147 Iterations: 190 Function evaluations: 449 --------------- Calculation Completed. Number Of Atoms in Supercell: 686 Vacancy Formation Energy (relaxed): 6.15614677248 Vacancy Formation Energy (unrelaxed): 8.8238085137 Unrelaxed Cell Volume: 10757.3238823 Relaxed Cell Volume: 10751.2048172 Relaxation Volume: 6.11906514535 Relaxed Cell Vector: [22.070847725070212, -8.32110599635245e-06, 22.070846910746923, 8.82898914963536e-06, -4.6553039415469825e-06, 22.070849842309656] Unrelaxed Cell Vector: [22.07503571364, 0.0, 22.07503571364, 0.0, 0.0, 22.07503571364] Relaxed Cell: [[ 2.20708477e+01 0.00000000e+00 0.00000000e+00] [-8.32110600e-06 2.20708469e+01 0.00000000e+00] [ 8.82898915e-06 -4.65530394e-06 2.20708498e+01]] Unrelaxed Cell: [[22.07503571 0. 0. ] [ 0. 22.07503571 0. ] [ 0. 0. 22.07503571]] [Calculation Results Summary] Sizes: [5, 6, 7] Unrelaxed Formation Energy By Size: [8.823808513723634, 8.823808513729091, 8.823808513704535] Formation Energy By Size: [6.294948981527341, 6.207087592379139, 6.156146772476859] Relaxation Volume By Size: [5.940364654345558, 6.06502402366641, 6.119065145352579] [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.82380851 8.82380851] Fitting Results: (array([ 8.82380851e+00, -1.61892750e-09]), 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.29494898 6.20708759] Fitting Results: (array([ 6.08639887, 26.06876381]), 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: [5.94036465 6.06502402] Fitting Results: (array([ 6.23625942, -36.98684584]), 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.82380851 8.82380851] Fitting Results: (array([8.82380851e+00, 1.43253792e-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: [6.20708759 6.15614677] Fitting Results: (array([ 6.06950727, 29.71735012]), 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.06502402 6.11906515] Fitting Results: (array([ 6.2109776 , -31.52597341]), 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.82380851 8.82380851 8.82380851] Fitting Results: (array([8.82380851e+00, 3.00965832e-09]), array([2.11325781e-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.29494898 6.20708759 6.15614677] Fitting Results: (array([ 6.07882641, 27.12781704]), array([1.10669001e-05]), 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: [5.94036465 6.06502402 6.11906515] Fitting Results: (array([ 6.22492564, -35.40175121]), array([2.47913858e-05]), 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.82380851 8.82380851 8.82380851] Fitting Results: (array([ 8.82380851e+00, 9.98356717e-08, -4.12770638e-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.29494898 6.20708759 6.15614677] Fitting Results: (array([ 6.05179833, 49.2857276 , -94.4594801 ]), 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: [5.94036465 6.06502402 6.11906515] Fitting Results: (array([ 6.18447247, -2.23780419, -141.37836618]), 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.82380851 8.82380851 8.82380851] Fitting Results: (array([ 8.82380851e+00, 5.33328836e-08, -9.67640207e-07]), 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.29494898 6.20708759 6.15614677] Fitting Results: (array([ 6.05665761, 38.64391159, -221.43724 ]), 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: [5.94036465 6.06502402 6.11906515] Fitting Results: (array([ 6.1917454 , -18.16550736, -331.42713859]), 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.82380851 8.82380851 8.82380851] Fitting Results: (array([ 8.82380851e+00, 3.79846107e-08, -3.13570724e-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.29494898 6.20708759 6.15614677] Fitting Results: (array([ 6.05982172, 35.13157387, -717.5832009 ]), 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: [5.94036465 6.06502402 6.11906515] Fitting Results: (array([ 6.19648115, -23.4224553 , -1074.01332753]), 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.823808513736587, 8.823808513662767]) list([8.823808513703495]) list([8.823808513585384]) list([8.823808513606618]) list([8.823808513620447])] Formation Energy Fits By Size: [list([6.086398871021719, 6.069507267761171]) list([6.078826406296333]) list([6.051798328858981]) list([6.056657605576108]) list([6.059821715432873])] Relaxation Volume Fits By Size: [list([6.236259421085164, 6.210977604283387]) list([6.224925644002421]) list([6.184472473765389]) list([6.191745397587755]) list([6.1964811496952095])] [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.823808513662767 "source-unit" "eV" "source-std-uncert-value" 5.061892807725925e-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.1535765305199996 "source-unit" "angstrom" } "host-b" { "source-value" 3.1535765305199996 "source-unit" "angstrom" } "host-c" { "source-value" 3.1535765305199996 "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.823808513730926 "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.1535765305199996 "source-unit" "angstrom" } "reservoir-b" { "source-value" 3.1535765305199996 "source-unit" "angstrom" } "reservoir-c" { "source-value" 3.1535765305199996 "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" 6.069507267761171 "source-unit" "eV" "source-std-uncert-value" 0.017708939625631873 } "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.1535765305199996 "source-unit" "angstrom" } "host-b" { "source-value" 3.1535765305199996 "source-unit" "angstrom" } "host-c" { "source-value" 3.1535765305199996 "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.823808513730926 "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.1535765305199996 "source-unit" "angstrom" } "reservoir-b" { "source-value" 3.1535765305199996 "source-unit" "angstrom" } "reservoir-c" { "source-value" 3.1535765305199996 "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.210977604283387 "source-unit" "angstrom^3" "source-std-uncert-value" 0.026556419133921243 } "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.1535765305199996 "source-unit" "angstrom" } "host-b" { "source-value" 3.1535765305199996 "source-unit" "angstrom" } "host-c" { "source-value" 3.1535765305199996 "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" ] } } ]