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: Ba bcc Morse_Shifted_GirifalcoWeizer_1959MedCutoff_Ba__MO_229241184339_004 [5.034461110830307] CELL_SIZE_MIN: 5 CELL_SIZE_MAX: 7 Smallest System Size: 250 Largest System Size: 686 [Calculation] Supercell Size: 5 Unrelaxed Cell: [[25.17230555 0. 0. ] [ 0. 25.17230555 0. ] [ 0. 0. 25.17230555]] Unrelaxed Cell Vector: [25.172305554151535, 0.0, 25.172305554151535, 0.0, 0.0, 25.172305554151535] Unrelaxed Cell Energy: -435.9917542026768 Energy of Unrelaxed Cell With Vacancy: -435.9917542026768 ========== Loop: 0 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 18:32:42 -432.503820* 0.3423 FIRE: 1 18:32:42 -432.513669* 0.3374 FIRE: 2 18:32:42 -432.532726* 0.3279 FIRE: 3 18:32:42 -432.559774* 0.3138 FIRE: 4 18:32:42 -432.593145* 0.2956 FIRE: 5 18:32:42 -432.630893* 0.2738 FIRE: 6 18:32:42 -432.670996* 0.2489 FIRE: 7 18:32:42 -432.711573* 0.2215 FIRE: 8 18:32:42 -432.755129* 0.1894 FIRE: 9 18:32:42 -432.799995* 0.1530 FIRE: 10 18:32:42 -432.844624* 0.1133 FIRE: 11 18:32:42 -432.888228* 0.1213 FIRE: 12 18:32:42 -432.931048* 0.1314 FIRE: 13 18:32:42 -432.973693* 0.1272 FIRE: 14 18:32:42 -433.015238* 0.1047 FIRE: 15 18:32:42 -433.052023* 0.0631 FIRE: 16 18:32:42 -433.078036* 0.0542 FIRE: 17 18:32:43 -433.093208* 0.0517 FIRE: 18 18:32:43 -433.101001* 0.0728 FIRE: 19 18:32:43 -433.103659* 0.0696 FIRE: 20 18:32:43 -433.108652* 0.0633 FIRE: 21 18:32:43 -433.115392* 0.0544 FIRE: 22 18:32:43 -433.123119* 0.0435 FIRE: 23 18:32:43 -433.131028* 0.0350 FIRE: 24 18:32:43 -433.138390* 0.0308 FIRE: 25 18:32:43 -433.144666* 0.0254 FIRE: 26 18:32:43 -433.149980* 0.0181 FIRE: 27 18:32:43 -433.153787* 0.0140 FIRE: 28 18:32:43 -433.155771* 0.0183 FIRE: 29 18:32:43 -433.155919* 0.0175 FIRE: 30 18:32:43 -433.156049* 0.0172 FIRE: 31 18:32:43 -433.156301* 0.0164 FIRE: 32 18:32:43 -433.156660* 0.0159 FIRE: 33 18:32:43 -433.157104* 0.0155 FIRE: 34 18:32:43 -433.157608* 0.0150 FIRE: 35 18:32:43 -433.158145* 0.0144 FIRE: 36 18:32:43 -433.158687* 0.0137 FIRE: 37 18:32:43 -433.159262* 0.0127 FIRE: 38 18:32:43 -433.159835* 0.0113 FIRE: 39 18:32:43 -433.160367* 0.0095 FIRE: 40 18:32:43 -433.160821* 0.0072 FIRE: 41 18:32:43 -433.161187* 0.0048 FIRE: 42 18:32:43 -433.161477* 0.0048 FIRE: 43 18:32:43 -433.161726* 0.0049 FIRE: 44 18:32:43 -433.161972* 0.0055 FIRE: 45 18:32:43 -433.162234* 0.0077 FIRE: 46 18:32:43 -433.162498* 0.0082 FIRE: 47 18:32:43 -433.162709* 0.0067 FIRE: 48 18:32:43 -433.162783* 0.0039 FIRE: 49 18:32:43 -433.162803* 0.0038 FIRE: 50 18:32:43 -433.162838* 0.0035 WARNING: Max number of steps exceeded. Structure may be unstable. Relaxation Completed. Steps: 50 Cell Size Relaxation... Optimization terminated successfully. Current function value: 1.083208 Iterations: 243 Function evaluations: 507 Current VFE: 1.0832080848375085 Energy of Supercell: -435.9917542026768 Unrelaxed Cell Volume: 15950.304719948243 Current Relaxed Cell Volume: 15919.885546579602 Current Relaxation Volume: 30.41917336864026 Current Cell: [[2.51562903e+01 0.00000000e+00 0.00000000e+00] [4.14702302e-05 2.51562937e+01 0.00000000e+00] [6.57692221e-05 6.77903976e-05 2.51562955e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 18:32:51 -433.164579* 0.0035 FIRE: 1 18:32:51 -433.164585* 0.0034 FIRE: 2 18:32:51 -433.164596* 0.0033 FIRE: 3 18:32:51 -433.164611* 0.0032 FIRE: 4 18:32:51 -433.164630* 0.0030 FIRE: 5 18:32:51 -433.164652* 0.0028 FIRE: 6 18:32:51 -433.164676* 0.0026 FIRE: 7 18:32:51 -433.164700* 0.0023 FIRE: 8 18:32:51 -433.164725* 0.0020 FIRE: 9 18:32:51 -433.164751* 0.0016 FIRE: 10 18:32:51 -433.164773* 0.0014 FIRE: 11 18:32:51 -433.164791* 0.0013 FIRE: 12 18:32:51 -433.164802* 0.0010 FIRE: 13 18:32:51 -433.164807* 0.0008 Relaxation Completed. Steps: 13 Cell Size Relaxation... Optimization terminated successfully. Current function value: 1.082980 Iterations: 214 Function evaluations: 438 Current VFE: 1.0829804653510564 Energy of Supercell: -435.9917542026768 Unrelaxed Cell Volume: 15950.304719948243 Current Relaxed Cell Volume: 15919.932891170507 Current Relaxation Volume: 30.371828777735573 Current Cell: [[2.51563215e+01 0.00000000e+00 0.00000000e+00] [6.27811829e-05 2.51563154e+01 0.00000000e+00] [2.76159139e-05 2.96972287e-06 2.51563174e+01]] ========== Loop: 2 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 18:32:58 -433.164807* 0.0008 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 1.082980 Iterations: 198 Function evaluations: 405 Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 18:33:05 -433.164807* 0.0008 FIRE: 1 18:33:05 -433.164807* 0.0008 FIRE: 2 18:33:05 -433.164808* 0.0007 FIRE: 3 18:33:05 -433.164810* 0.0007 FIRE: 4 18:33:05 -433.164812* 0.0007 FIRE: 5 18:33:05 -433.164814* 0.0007 FIRE: 6 18:33:05 -433.164816* 0.0006 FIRE: 7 18:33:05 -433.164818* 0.0006 FIRE: 8 18:33:05 -433.164821* 0.0006 FIRE: 9 18:33:05 -433.164823* 0.0005 FIRE: 10 18:33:05 -433.164825* 0.0004 FIRE: 11 18:33:05 -433.164827* 0.0004 FIRE: 12 18:33:05 -433.164829* 0.0003 FIRE: 13 18:33:05 -433.164831* 0.0004 FIRE: 14 18:33:05 -433.164833* 0.0004 FIRE: 15 18:33:05 -433.164835* 0.0003 FIRE: 16 18:33:05 -433.164837* 0.0003 FIRE: 17 18:33:05 -433.164839* 0.0002 FIRE: 18 18:33:05 -433.164840* 0.0003 FIRE: 19 18:33:05 -433.164840* 0.0003 FIRE: 20 18:33:05 -433.164840* 0.0003 Optimization terminated successfully. Current function value: 1.082947 Iterations: 302 Function evaluations: 615 --------------- Calculation Completed. Number Of Atoms in Supercell: 250 Vacancy Formation Energy (relaxed): 1.0829471148403513 Vacancy Formation Energy (unrelaxed): 1.7439670168107 Unrelaxed Cell Volume: 15950.304719948243 Relaxed Cell Volume: 15919.932891170507 Relaxation Volume: 30.371828777735573 Relaxed Cell Vector: [25.15632256647084, 3.3826441376424043e-07, 25.156328612796717, 1.473455687348461e-06, 8.071954361051523e-06, 25.156330183278104] Unrelaxed Cell Vector: [25.172305554151535, 0.0, 25.172305554151535, 0.0, 0.0, 25.172305554151535] Relaxed Cell: [[2.51563226e+01 0.00000000e+00 0.00000000e+00] [3.38264414e-07 2.51563286e+01 0.00000000e+00] [1.47345569e-06 8.07195436e-06 2.51563302e+01]] Unrelaxed Cell: [[25.17230555 0. 0. ] [ 0. 25.17230555 0. ] [ 0. 0. 25.17230555]] Supercell Size: 6 Unrelaxed Cell: [[30.20676666 0. 0. ] [ 0. 30.20676666 0. ] [ 0. 0. 30.20676666]] Unrelaxed Cell Vector: [30.206766664981842, 0.0, 30.206766664981842, 0.0, 0.0, 30.206766664981842] Unrelaxed Cell Energy: -753.3937512619738 Energy of Unrelaxed Cell With Vacancy: -753.3937512619738 ========== Loop: 0 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 18:33:33 -749.905817* 0.3423 FIRE: 1 18:33:33 -749.915667* 0.3374 FIRE: 2 18:33:33 -749.934726* 0.3279 FIRE: 3 18:33:33 -749.961778* 0.3138 FIRE: 4 18:33:33 -749.995154* 0.2956 FIRE: 5 18:33:33 -750.032910* 0.2738 FIRE: 6 18:33:33 -750.073025* 0.2489 FIRE: 7 18:33:33 -750.113618* 0.2215 FIRE: 8 18:33:33 -750.157198* 0.1895 FIRE: 9 18:33:33 -750.202097* 0.1530 FIRE: 10 18:33:33 -750.246769* 0.1133 FIRE: 11 18:33:33 -750.290423* 0.1213 FIRE: 12 18:33:33 -750.333292* 0.1315 FIRE: 13 18:33:33 -750.375966* 0.1274 FIRE: 14 18:33:33 -750.417506* 0.1049 FIRE: 15 18:33:34 -750.454252* 0.0634 FIRE: 16 18:33:34 -750.480208* 0.0542 FIRE: 17 18:33:34 -750.495280* 0.0517 FIRE: 18 18:33:34 -750.502824* 0.0725 FIRE: 19 18:33:34 -750.505450* 0.0693 FIRE: 20 18:33:34 -750.510388* 0.0631 FIRE: 21 18:33:34 -750.517063* 0.0543 FIRE: 22 18:33:34 -750.524741* 0.0433 FIRE: 23 18:33:34 -750.532647* 0.0351 FIRE: 24 18:33:34 -750.540100* 0.0309 FIRE: 25 18:33:34 -750.546606* 0.0256 FIRE: 26 18:33:34 -750.552406* 0.0184 FIRE: 27 18:33:34 -750.557088* 0.0139 FIRE: 28 18:33:34 -750.560449* 0.0187 FIRE: 29 18:33:34 -750.562475* 0.0186 FIRE: 30 18:33:34 -750.563393* 0.0270 FIRE: 31 18:33:35 -750.563612* 0.0267 FIRE: 32 18:33:35 -750.564039* 0.0260 FIRE: 33 18:33:35 -750.564647* 0.0249 FIRE: 34 18:33:35 -750.565400* 0.0236 FIRE: 35 18:33:35 -750.566256* 0.0219 FIRE: 36 18:33:35 -750.567169* 0.0199 FIRE: 37 18:33:35 -750.568091* 0.0177 FIRE: 38 18:33:35 -750.569071* 0.0150 FIRE: 39 18:33:35 -750.570050* 0.0118 FIRE: 40 18:33:35 -750.570967* 0.0088 FIRE: 41 18:33:35 -750.571754* 0.0077 FIRE: 42 18:33:35 -750.572394* 0.0085 FIRE: 43 18:33:35 -750.572943* 0.0089 FIRE: 44 18:33:35 -750.573501* 0.0085 FIRE: 45 18:33:35 -750.574173* 0.0084 FIRE: 46 18:33:35 -750.574995* 0.0087 FIRE: 47 18:33:35 -750.575850* 0.0068 FIRE: 48 18:33:35 -750.576457* 0.0035 FIRE: 49 18:33:35 -750.576539* 0.0042 FIRE: 50 18:33:35 -750.576563* 0.0040 WARNING: Max number of steps exceeded. Structure may be unstable. Relaxation Completed. Steps: 50 Cell Size Relaxation... Optimization terminated successfully. Current function value: 1.072194 Iterations: 185 Function evaluations: 411 Current VFE: 1.0721935676968997 Energy of Supercell: -753.3937512619738 Unrelaxed Cell Volume: 27562.126556070587 Current Relaxed Cell Volume: 27531.448459815558 Current Relaxation Volume: 30.67809625502923 Current Cell: [[3.01955564e+01 0.00000000e+00 0.00000000e+00] [5.62919415e-05 3.01955572e+01 0.00000000e+00] [2.25181009e-05 8.43742982e-05 3.01955522e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 18:34:16 -750.577591* 0.0040 FIRE: 1 18:34:16 -750.577598* 0.0040 FIRE: 2 18:34:16 -750.577611* 0.0039 FIRE: 3 18:34:16 -750.577631* 0.0037 FIRE: 4 18:34:17 -750.577655* 0.0035 FIRE: 5 18:34:17 -750.577683* 0.0033 FIRE: 6 18:34:17 -750.577712* 0.0030 FIRE: 7 18:34:17 -750.577742* 0.0027 FIRE: 8 18:34:17 -750.577774* 0.0023 FIRE: 9 18:34:17 -750.577805* 0.0019 FIRE: 10 18:34:17 -750.577831* 0.0014 FIRE: 11 18:34:17 -750.577850* 0.0009 Relaxation Completed. Steps: 11 Cell Size Relaxation... Optimization terminated successfully. Current function value: 1.071934 Iterations: 252 Function evaluations: 487 Current VFE: 1.0719340756022575 Energy of Supercell: -753.3937512619738 Unrelaxed Cell Volume: 27562.126556070587 Current Relaxed Cell Volume: 27531.425297573835 Current Relaxation Volume: 30.701258496752416 Current Cell: [[ 3.01955428e+01 0.00000000e+00 0.00000000e+00] [-2.28013656e-06 3.01955492e+01 0.00000000e+00] [ 2.58285940e-05 1.49679109e-06 3.01955484e+01]] ========== Loop: 2 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 18:34:48 -750.577850* 0.0009 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 1.071934 Iterations: 175 Function evaluations: 368 Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 18:35:18 -750.577850* 0.0009 FIRE: 1 18:35:19 -750.577851* 0.0009 FIRE: 2 18:35:19 -750.577851* 0.0009 FIRE: 3 18:35:19 -750.577852* 0.0008 FIRE: 4 18:35:19 -750.577854* 0.0008 FIRE: 5 18:35:19 -750.577855* 0.0008 FIRE: 6 18:35:19 -750.577857* 0.0007 FIRE: 7 18:35:19 -750.577859* 0.0007 FIRE: 8 18:35:19 -750.577861* 0.0007 FIRE: 9 18:35:19 -750.577863* 0.0006 FIRE: 10 18:35:19 -750.577866* 0.0006 FIRE: 11 18:35:19 -750.577869* 0.0006 FIRE: 12 18:35:19 -750.577872* 0.0005 FIRE: 13 18:35:19 -750.577874* 0.0004 FIRE: 14 18:35:19 -750.577877* 0.0003 FIRE: 15 18:35:19 -750.577878* 0.0002 FIRE: 16 18:35:19 -750.577880* 0.0003 FIRE: 17 18:35:19 -750.577881* 0.0003 FIRE: 18 18:35:19 -750.577881* 0.0003 FIRE: 19 18:35:19 -750.577882* 0.0003 FIRE: 20 18:35:19 -750.577882* 0.0003 Optimization terminated successfully. Current function value: 1.071902 Iterations: 173 Function evaluations: 435 --------------- Calculation Completed. Number Of Atoms in Supercell: 432 Vacancy Formation Energy (relaxed): 1.0719020129325827 Vacancy Formation Energy (unrelaxed): 1.7439670168085968 Unrelaxed Cell Volume: 27562.126556070587 Relaxed Cell Volume: 27531.425297573835 Relaxation Volume: 30.701258496752416 Relaxed Cell Vector: [30.195560640464485, -2.351240213015531e-06, 30.195558241338674, 2.4127638117785264e-05, 1.5825571330647567e-06, 30.19555870299167] Unrelaxed Cell Vector: [30.206766664981842, 0.0, 30.206766664981842, 0.0, 0.0, 30.206766664981842] Relaxed Cell: [[ 3.01955606e+01 0.00000000e+00 0.00000000e+00] [-2.35124021e-06 3.01955582e+01 0.00000000e+00] [ 2.41276381e-05 1.58255713e-06 3.01955587e+01]] Unrelaxed Cell: [[30.20676666 0. 0. ] [ 0. 30.20676666 0. ] [ 0. 0. 30.20676666]] Supercell Size: 7 Unrelaxed Cell: [[35.24122778 0. 0. ] [ 0. 35.24122778 0. ] [ 0. 0. 35.24122778]] Unrelaxed Cell Vector: [35.24122777581215, 0.0, 35.24122777581215, 0.0, 0.0, 35.24122777581215] Unrelaxed Cell Energy: -1196.3613735324586 Energy of Unrelaxed Cell With Vacancy: -1196.3613735324586 ========== Loop: 0 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 18:35:41 -1192.873439* 0.3423 FIRE: 1 18:35:41 -1192.883289* 0.3374 FIRE: 2 18:35:41 -1192.902348* 0.3279 FIRE: 3 18:35:41 -1192.929400* 0.3138 FIRE: 4 18:35:41 -1192.962777* 0.2956 FIRE: 5 18:35:41 -1193.000532* 0.2738 FIRE: 6 18:35:41 -1193.040647* 0.2489 FIRE: 7 18:35:41 -1193.081241* 0.2215 FIRE: 8 18:35:41 -1193.124821* 0.1895 FIRE: 9 18:35:41 -1193.169722* 0.1530 FIRE: 10 18:35:41 -1193.214398* 0.1133 FIRE: 11 18:35:41 -1193.258061* 0.1213 FIRE: 12 18:35:41 -1193.300945* 0.1315 FIRE: 13 18:35:41 -1193.343647* 0.1273 FIRE: 14 18:35:41 -1193.385225* 0.1049 FIRE: 15 18:35:41 -1193.422013* 0.0634 FIRE: 16 18:35:41 -1193.447992* 0.0542 FIRE: 17 18:35:41 -1193.463066* 0.0517 FIRE: 18 18:35:41 -1193.470580* 0.0724 FIRE: 19 18:35:41 -1193.473205* 0.0692 FIRE: 20 18:35:41 -1193.478140* 0.0630 FIRE: 21 18:35:41 -1193.484811* 0.0541 FIRE: 22 18:35:41 -1193.492482* 0.0432 FIRE: 23 18:35:41 -1193.500378* 0.0350 FIRE: 24 18:35:41 -1193.507811* 0.0309 FIRE: 25 18:35:41 -1193.514287* 0.0256 FIRE: 26 18:35:41 -1193.520037* 0.0184 FIRE: 27 18:35:41 -1193.524660* 0.0140 FIRE: 28 18:35:41 -1193.527979* 0.0188 FIRE: 29 18:35:41 -1193.530071* 0.0187 FIRE: 30 18:35:41 -1193.531305* 0.0270 FIRE: 31 18:35:41 -1193.532455* 0.0360 FIRE: 32 18:35:41 -1193.534500* 0.0395 FIRE: 33 18:35:41 -1193.537963* 0.0354 FIRE: 34 18:35:41 -1193.542139* 0.0234 FIRE: 35 18:35:41 -1193.545310* 0.0099 FIRE: 36 18:35:41 -1193.546266* 0.0126 FIRE: 37 18:35:41 -1193.546629* 0.0113 FIRE: 38 18:35:41 -1193.547257* 0.0087 FIRE: 39 18:35:41 -1193.547987* 0.0056 FIRE: 40 18:35:41 -1193.548653* 0.0041 FIRE: 41 18:35:41 -1193.549146* 0.0030 FIRE: 42 18:35:41 -1193.549446* 0.0057 FIRE: 43 18:35:41 -1193.549596* 0.0075 FIRE: 44 18:35:41 -1193.549652* 0.0077 FIRE: 45 18:35:42 -1193.549667* 0.0076 FIRE: 46 18:35:42 -1193.549695* 0.0072 FIRE: 47 18:35:42 -1193.549734* 0.0066 FIRE: 48 18:35:42 -1193.549781* 0.0059 FIRE: 49 18:35:42 -1193.549831* 0.0051 FIRE: 50 18:35:42 -1193.549879* 0.0041 WARNING: Max number of steps exceeded. Structure may be unstable. Relaxation Completed. Steps: 50 Cell Size Relaxation... Optimization terminated successfully. Current function value: 1.066877 Iterations: 255 Function evaluations: 511 Current VFE: 1.0668765581253865 Energy of Supercell: -1196.3613735324586 Unrelaxed Cell Volume: 43767.63615153794 Current Relaxed Cell Volume: 43736.89436620522 Current Relaxation Volume: 30.74178533272061 Current Cell: [[3.52329827e+01 0.00000000e+00 0.00000000e+00] [9.50751515e-06 3.52329719e+01 0.00000000e+00] [3.63748780e-05 5.15975903e-05 3.52329700e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 18:36:32 -1193.550530* 0.0043 FIRE: 1 18:36:32 -1193.550534* 0.0042 FIRE: 2 18:36:32 -1193.550542* 0.0041 FIRE: 3 18:36:32 -1193.550554* 0.0038 FIRE: 4 18:36:32 -1193.550568* 0.0036 FIRE: 5 18:36:32 -1193.550585* 0.0032 FIRE: 6 18:36:32 -1193.550602* 0.0028 FIRE: 7 18:36:32 -1193.550621* 0.0024 FIRE: 8 18:36:32 -1193.550640* 0.0019 FIRE: 9 18:36:32 -1193.550660* 0.0013 FIRE: 10 18:36:32 -1193.550680* 0.0006 Relaxation Completed. Steps: 10 Cell Size Relaxation... Optimization terminated successfully. Current function value: 1.066727 Iterations: 208 Function evaluations: 439 Current VFE: 1.066726764151099 Energy of Supercell: -1196.3613735324586 Unrelaxed Cell Volume: 43767.63615153794 Current Relaxed Cell Volume: 43736.83558360654 Current Relaxation Volume: 30.80056793140102 Current Cell: [[ 3.52329581e+01 0.00000000e+00 0.00000000e+00] [ 1.62833457e-05 3.52329598e+01 0.00000000e+00] [ 2.99521994e-05 -6.13806230e-07 3.52329594e+01]] ========== Loop: 2 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 18:37:06 -1193.550680* 0.0006 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 1.066727 Iterations: 195 Function evaluations: 407 Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 18:37:30 -1193.550680* 0.0006 FIRE: 1 18:37:30 -1193.550681* 0.0006 FIRE: 2 18:37:30 -1193.550682* 0.0006 FIRE: 3 18:37:30 -1193.550684* 0.0006 FIRE: 4 18:37:30 -1193.550687* 0.0006 FIRE: 5 18:37:30 -1193.550691* 0.0006 FIRE: 6 18:37:30 -1193.550695* 0.0005 FIRE: 7 18:37:30 -1193.550699* 0.0005 FIRE: 8 18:37:30 -1193.550704* 0.0005 FIRE: 9 18:37:30 -1193.550710* 0.0005 FIRE: 10 18:37:30 -1193.550718* 0.0004 FIRE: 11 18:37:30 -1193.550726* 0.0004 FIRE: 12 18:37:30 -1193.550734* 0.0004 FIRE: 13 18:37:30 -1193.550744* 0.0003 FIRE: 14 18:37:30 -1193.550753* 0.0003 FIRE: 15 18:37:30 -1193.550763* 0.0003 FIRE: 16 18:37:30 -1193.550771* 0.0004 FIRE: 17 18:37:30 -1193.550778* 0.0004 FIRE: 18 18:37:30 -1193.550783* 0.0003 FIRE: 19 18:37:30 -1193.550784* 0.0002 FIRE: 20 18:37:30 -1193.550784* 0.0002 Optimization terminated successfully. Current function value: 1.066622 Iterations: 204 Function evaluations: 468 --------------- Calculation Completed. Number Of Atoms in Supercell: 686 Vacancy Formation Energy (relaxed): 1.0666221180756565 Vacancy Formation Energy (unrelaxed): 1.7439670168221255 Unrelaxed Cell Volume: 43767.63615153794 Relaxed Cell Volume: 43736.83558360654 Relaxation Volume: 30.80056793140102 Relaxed Cell Vector: [35.23295292275895, 1.6654827330863574e-05, 35.232943392151114, 2.876895345392426e-05, -6.352542271557061e-07, 35.232954237399156] Unrelaxed Cell Vector: [35.24122777581215, 0.0, 35.24122777581215, 0.0, 0.0, 35.24122777581215] Relaxed Cell: [[ 3.52329529e+01 0.00000000e+00 0.00000000e+00] [ 1.66548273e-05 3.52329434e+01 0.00000000e+00] [ 2.87689535e-05 -6.35254227e-07 3.52329542e+01]] Unrelaxed Cell: [[35.24122778 0. 0. ] [ 0. 35.24122778 0. ] [ 0. 0. 35.24122778]] [Calculation Results Summary] Sizes: [5, 6, 7] Unrelaxed Formation Energy By Size: [1.7439670168107, 1.7439670168085968, 1.7439670168221255] Formation Energy By Size: [1.0829471148403513, 1.0719020129325827, 1.0666221180756565] Relaxation Volume By Size: [30.371828777735573, 30.701258496752416, 30.80056793140102] [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: [1.74396702 1.74396702] Fitting Results: (array([1.74396702e+00, 6.24054815e-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: [1.08294711 1.07190201] Fitting Results: (array([1.05673017, 3.27711815]), 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: [30.37182878 30.7012585 ] Fitting Results: (array([ 31.15377185, -97.74288366]), 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: [1.74396702 1.74396702] Fitting Results: (array([ 1.74396702e+00, -7.89251454e-09]), 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: [1.07190201 1.06662212] Fitting Results: (array([1.05764214, 3.08013268]), 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: [30.7012585 30.80056793] Fitting Results: (array([ 30.96947217, -57.93415271]), 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: [1.74396702 1.74396702 1.74396702] Fitting Results: (array([ 1.74396702e+00, -1.84796440e-09]), array([6.02940186e-23]), 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: [1.08294711 1.07190201 1.06662212] Fitting Results: (array([1.057139 , 3.21994036]), array([3.22585713e-08]), 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: [30.37182878 30.7012585 30.80056793] Fitting Results: (array([ 31.07115075, -86.18784248]), array([0.00131745]), 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: [1.74396702 1.74396702 1.74396702] Fitting Results: (array([ 1.74396702e+00, -5.35672986e-08, 2.20480239e-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: [1.08294711 1.07190201 1.06662212] Fitting Results: (array([1.05859824, 2.02364512, 5.0998232 ]), 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: [30.37182878 30.7012585 30.80056793] Fitting Results: (array([ 30.77625476, 155.57108181, -1030.62164792]), 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: [1.74396702 1.74396702 1.74396702] Fitting Results: (array([ 1.74396702e+00, -2.87279692e-08, 5.16862210e-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: [1.08294711 1.07190201 1.06662212] Fitting Results: (array([ 1.05833589, 2.59819179, 11.95529314]), 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: [30.37182878 30.7012585 30.80056793] Fitting Results: (array([ 30.829273 , 39.46112714, -2416.04138566]), 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: [1.74396702 1.74396702 1.74396702] Fitting Results: (array([ 1.74396702e+00, -2.05297337e-08, 1.67492893e-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: [1.08294711 1.07190201 1.06662212] Fitting Results: (array([ 1.05816506, 2.78782127, 38.74198178]), 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: [30.37182878 30.7012585 30.80056793] Fitting Results: (array([ 3.08637957e+01, 1.13896731e+00, -7.82935477e+03]), 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: [[0, 3], [0, 3], [0, 3, 4], [0, 3, 5], [0, 3, 6]] Unrelaxed Formation Energy Fits By Size: [[1.7439670168057078, 1.7439670168451353], [1.7439670168233834], [1.7439670168864694], [1.7439670168751273], [1.7439670168677426]] Formation Energy Fits By Size: [[1.0567301696526812, 1.0576421394213573], [1.057139003485527], [1.0585982367594844], [1.058335886679667], [1.058165057866699]] Relaxation Volume Fits By Size: [[31.153771847050283, 30.969472166708883], [31.071150747871073], [30.776254759899363], [30.82927300404371], [30.863795744494798]] [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" 1.7439670168451353 "source-unit" "eV" "source-std-uncert-value" 0.00010464607544236876 } "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" 5.034461110830307 "source-unit" "angstrom" } "host-b" { "source-value" 5.034461110830307 "source-unit" "angstrom" } "host-c" { "source-value" 5.034461110830307 "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" [ "Ba" ] } "reservoir-cohesive-potential-energy" { "source-value" 1.7439670168109402 "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" 5.034461110830307 "source-unit" "angstrom" } "reservoir-b" { "source-value" 5.034461110830307 "source-unit" "angstrom" } "reservoir-c" { "source-value" 5.034461110830307 "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" [ "Ba" ] } } { "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" 1.0576421394213573 "source-unit" "eV" "source-std-uncert-value" 0.0009618071122003408 } "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" 5.034461110830307 "source-unit" "angstrom" } "host-b" { "source-value" 5.034461110830307 "source-unit" "angstrom" } "host-c" { "source-value" 5.034461110830307 "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" [ "Ba" ] } "reservoir-cohesive-potential-energy" { "source-value" 1.7439670168109402 "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" 5.034461110830307 "source-unit" "angstrom" } "reservoir-b" { "source-value" 5.034461110830307 "source-unit" "angstrom" } "reservoir-c" { "source-value" 5.034461110830307 "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" [ "Ba" ] } } { "property-id" "tag:staff@noreply.openkim.org,2015-07-28:property/monovacancy-neutral-relaxation-volume-crystal-npt" "instance-id" 3 "relaxation-volume" { "source-value" 30.969472166708883 "source-unit" "angstrom^3" "source-std-uncert-value" 0.196033436607354 } "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" 5.034461110830307 "source-unit" "angstrom" } "host-b" { "source-value" 5.034461110830307 "source-unit" "angstrom" } "host-c" { "source-value" 5.034461110830307 "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" [ "Ba" ] } } ]