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: Fe bcc Morse_Shifted_GirifalcoWeizer_1959MedCutoff_Fe__MO_984358344196_004 [2.868830233812332] CELL_SIZE_MIN: 5 CELL_SIZE_MAX: 7 Smallest System Size: 250 Largest System Size: 686 [Calculation] Supercell Size: 5 Unrelaxed Cell: [[14.34415117 0. 0. ] [ 0. 14.34415117 0. ] [ 0. 0. 14.34415117]] Unrelaxed Cell Vector: [14.34415116906166, 0.0, 14.34415116906166, 0.0, 0.0, 14.34415116906166] Unrelaxed Cell Energy: -1032.6051866407813 Energy of Unrelaxed Cell With Vacancy: -1032.6051866407813 ========== Loop: 0 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:06:42 -1024.344345* 1.2423 FIRE: 1 16:06:42 -1024.465542* 1.1061 FIRE: 2 16:06:42 -1024.655582* 0.8642 FIRE: 3 16:06:42 -1024.851080* 0.5686 FIRE: 4 16:06:42 -1025.020836* 0.4962 FIRE: 5 16:06:42 -1025.167572* 0.4804 FIRE: 6 16:06:42 -1025.293926* 0.3344 FIRE: 7 16:06:42 -1025.388174* 0.2166 FIRE: 8 16:06:42 -1025.446614* 0.2113 FIRE: 9 16:06:42 -1025.469823* 0.3462 FIRE: 10 16:06:42 -1025.479231* 0.3236 FIRE: 11 16:06:42 -1025.496395* 0.2802 FIRE: 12 16:06:42 -1025.518432* 0.2203 FIRE: 13 16:06:42 -1025.541884* 0.1609 FIRE: 14 16:06:43 -1025.563570* 0.1398 FIRE: 15 16:06:43 -1025.581231* 0.1110 FIRE: 16 16:06:43 -1025.593861* 0.0732 FIRE: 17 16:06:43 -1025.602075* 0.0813 FIRE: 18 16:06:43 -1025.605163* 0.0881 FIRE: 19 16:06:43 -1025.605455* 0.0857 FIRE: 20 16:06:43 -1025.606019* 0.0812 FIRE: 21 16:06:43 -1025.606815* 0.0745 FIRE: 22 16:06:43 -1025.607789* 0.0660 FIRE: 23 16:06:43 -1025.608880* 0.0608 FIRE: 24 16:06:43 -1025.610021* 0.0561 FIRE: 25 16:06:43 -1025.611153* 0.0510 FIRE: 26 16:06:43 -1025.612336* 0.0456 FIRE: 27 16:06:43 -1025.613509* 0.0444 FIRE: 28 16:06:43 -1025.614619* 0.0415 FIRE: 29 16:06:43 -1025.615642* 0.0360 FIRE: 30 16:06:43 -1025.616598* 0.0276 FIRE: 31 16:06:43 -1025.617528* 0.0287 FIRE: 32 16:06:43 -1025.618455* 0.0234 FIRE: 33 16:06:43 -1025.619351* 0.0180 FIRE: 34 16:06:43 -1025.620153* 0.0231 FIRE: 35 16:06:43 -1025.620810* 0.0261 FIRE: 36 16:06:43 -1025.621262* 0.0210 FIRE: 37 16:06:43 -1025.621347* 0.0129 FIRE: 38 16:06:43 -1025.621382* 0.0122 FIRE: 39 16:06:43 -1025.621443* 0.0109 FIRE: 40 16:06:43 -1025.621520* 0.0090 FIRE: 41 16:06:43 -1025.621599* 0.0066 FIRE: 42 16:06:43 -1025.621666* 0.0064 FIRE: 43 16:06:43 -1025.621715* 0.0065 FIRE: 44 16:06:43 -1025.621744* 0.0064 FIRE: 45 16:06:43 -1025.621762* 0.0059 FIRE: 46 16:06:43 -1025.621777* 0.0053 FIRE: 47 16:06:43 -1025.621796* 0.0059 FIRE: 48 16:06:43 -1025.621821* 0.0051 FIRE: 49 16:06:43 -1025.621841* 0.0044 FIRE: 50 16:06:43 -1025.621847* 0.0046 WARNING: Max number of steps exceeded. Structure may be unstable. Relaxation Completed. Steps: 50 Cell Size Relaxation... Optimization terminated successfully. Current function value: 2.848605 Iterations: 460 Function evaluations: 795 Current VFE: 2.8486050193878327 Energy of Supercell: -1032.6051866407813 Unrelaxed Cell Volume: 2951.376129823061 Current Relaxed Cell Volume: 2946.4839115499435 Current Relaxation Volume: 4.892218273117578 Current Cell: [[ 1.43362228e+01 0.00000000e+00 0.00000000e+00] [-2.99841414e-07 1.43362226e+01 0.00000000e+00] [-1.26143721e-06 2.38107893e-06 1.43362180e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:07:03 -1025.626161* 0.0046 FIRE: 1 16:07:03 -1025.626169* 0.0041 FIRE: 2 16:07:03 -1025.626183* 0.0033 FIRE: 3 16:07:03 -1025.626196* 0.0026 FIRE: 4 16:07:03 -1025.626205* 0.0017 FIRE: 5 16:07:03 -1025.626209* 0.0011 FIRE: 6 16:07:03 -1025.626208* 0.0015 FIRE: 7 16:07:03 -1025.626208* 0.0014 FIRE: 8 16:07:03 -1025.626208* 0.0014 FIRE: 9 16:07:03 -1025.626209* 0.0012 FIRE: 10 16:07:03 -1025.626209* 0.0011 FIRE: 11 16:07:03 -1025.626210* 0.0009 Relaxation Completed. Steps: 11 Cell Size Relaxation... Optimization terminated successfully. Current function value: 2.848556 Iterations: 148 Function evaluations: 336 Current VFE: 2.8485560127467124 Energy of Supercell: -1032.6051866407813 Unrelaxed Cell Volume: 2951.376129823061 Current Relaxed Cell Volume: 2946.481535633809 Current Relaxation Volume: 4.894594189252075 Current Cell: [[ 1.43362178e+01 0.00000000e+00 0.00000000e+00] [-2.94072691e-07 1.43362172e+01 0.00000000e+00] [-1.28557485e-06 2.42810708e-06 1.43362168e+01]] ========== Loop: 2 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:07:07 -1025.626210* 0.0009 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 2.848556 Iterations: 138 Function evaluations: 318 Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:07:10 -1025.626210* 0.0009 FIRE: 1 16:07:10 -1025.626210* 0.0008 FIRE: 2 16:07:10 -1025.626211* 0.0006 FIRE: 3 16:07:10 -1025.626211* 0.0005 FIRE: 4 16:07:10 -1025.626211* 0.0004 FIRE: 5 16:07:10 -1025.626212* 0.0002 FIRE: 6 16:07:10 -1025.626212* 0.0003 FIRE: 7 16:07:10 -1025.626212* 0.0003 FIRE: 8 16:07:10 -1025.626212* 0.0003 FIRE: 9 16:07:10 -1025.626212* 0.0003 FIRE: 10 16:07:10 -1025.626212* 0.0003 FIRE: 11 16:07:10 -1025.626212* 0.0002 FIRE: 12 16:07:10 -1025.626212* 0.0002 FIRE: 13 16:07:10 -1025.626212* 0.0002 FIRE: 14 16:07:10 -1025.626212* 0.0001 FIRE: 15 16:07:10 -1025.626212* 0.0001 FIRE: 16 16:07:10 -1025.626212* 0.0001 FIRE: 17 16:07:10 -1025.626212* 0.0001 FIRE: 18 16:07:10 -1025.626212* 0.0001 FIRE: 19 16:07:10 -1025.626212* 0.0001 FIRE: 20 16:07:10 -1025.626212* 0.0001 Optimization terminated successfully. Current function value: 2.848554 Iterations: 181 Function evaluations: 429 --------------- Calculation Completed. Number Of Atoms in Supercell: 250 Vacancy Formation Energy (relaxed): 2.848554117687854 Vacancy Formation Energy (unrelaxed): 4.13042074655209 Unrelaxed Cell Volume: 2951.376129823061 Relaxed Cell Volume: 2946.481535633809 Relaxation Volume: 4.894594189252075 Relaxed Cell Vector: [14.336219404090834, -3.006386256033498e-07, 14.336217582014694, -1.3060536770738275e-06, 2.4700365737535295e-06, 14.33622031134729] Unrelaxed Cell Vector: [14.34415116906166, 0.0, 14.34415116906166, 0.0, 0.0, 14.34415116906166] Relaxed Cell: [[ 1.43362194e+01 0.00000000e+00 0.00000000e+00] [-3.00638626e-07 1.43362176e+01 0.00000000e+00] [-1.30605368e-06 2.47003657e-06 1.43362203e+01]] Unrelaxed Cell: [[14.34415117 0. 0. ] [ 0. 14.34415117 0. ] [ 0. 0. 14.34415117]] Supercell Size: 6 Unrelaxed Cell: [[17.2129814 0. 0. ] [ 0. 17.2129814 0. ] [ 0. 0. 17.2129814]] Unrelaxed Cell Vector: [17.212981402873993, 0.0, 17.212981402873993, 0.0, 0.0, 17.212981402873993] Unrelaxed Cell Energy: -1784.3417625149036 Energy of Unrelaxed Cell With Vacancy: -1784.3417625149036 ========== Loop: 0 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:07:21 -1776.080921* 1.2423 FIRE: 1 16:07:21 -1776.202119* 1.1061 FIRE: 2 16:07:21 -1776.392163* 0.8642 FIRE: 3 16:07:22 -1776.587682* 0.5686 FIRE: 4 16:07:22 -1776.757499* 0.4961 FIRE: 5 16:07:22 -1776.904344* 0.4803 FIRE: 6 16:07:22 -1777.030824* 0.3346 FIRE: 7 16:07:22 -1777.125171* 0.2165 FIRE: 8 16:07:22 -1777.183729* 0.2110 FIRE: 9 16:07:22 -1777.207005* 0.3455 FIRE: 10 16:07:22 -1777.216354* 0.3229 FIRE: 11 16:07:22 -1777.233415* 0.2798 FIRE: 12 16:07:22 -1777.255351* 0.2200 FIRE: 13 16:07:22 -1777.278823* 0.1609 FIRE: 14 16:07:22 -1777.300809* 0.1399 FIRE: 15 16:07:22 -1777.319299* 0.1114 FIRE: 16 16:07:22 -1777.333575* 0.0744 FIRE: 17 16:07:22 -1777.344817* 0.0823 FIRE: 18 16:07:22 -1777.352691* 0.0926 FIRE: 19 16:07:22 -1777.357205* 0.0993 FIRE: 20 16:07:22 -1777.359173* 0.1541 FIRE: 21 16:07:22 -1777.359903* 0.1513 FIRE: 22 16:07:22 -1777.361306* 0.1456 FIRE: 23 16:07:22 -1777.363272* 0.1372 FIRE: 24 16:07:22 -1777.365655* 0.1264 FIRE: 25 16:07:22 -1777.368285* 0.1135 FIRE: 26 16:07:22 -1777.370988* 0.0988 FIRE: 27 16:07:22 -1777.373604* 0.0828 FIRE: 28 16:07:22 -1777.376238* 0.0643 FIRE: 29 16:07:22 -1777.378693* 0.0437 FIRE: 30 16:07:22 -1777.380770* 0.0379 FIRE: 31 16:07:22 -1777.382342* 0.0339 FIRE: 32 16:07:22 -1777.383445* 0.0350 FIRE: 33 16:07:23 -1777.384322* 0.0334 FIRE: 34 16:07:23 -1777.385339* 0.0390 FIRE: 35 16:07:23 -1777.386768* 0.0405 FIRE: 36 16:07:23 -1777.388520* 0.0349 FIRE: 37 16:07:23 -1777.390039* 0.0225 FIRE: 38 16:07:23 -1777.390662* 0.0187 FIRE: 39 16:07:23 -1777.390746* 0.0179 FIRE: 40 16:07:23 -1777.390898* 0.0162 FIRE: 41 16:07:23 -1777.391089* 0.0138 FIRE: 42 16:07:23 -1777.391285* 0.0108 FIRE: 43 16:07:23 -1777.391458* 0.0083 FIRE: 44 16:07:23 -1777.391586* 0.0094 FIRE: 45 16:07:23 -1777.391666* 0.0102 FIRE: 46 16:07:23 -1777.391702* 0.0086 FIRE: 47 16:07:23 -1777.391695* 0.0068 FIRE: 48 16:07:23 -1777.391698* 0.0066 FIRE: 49 16:07:23 -1777.391704* 0.0064 FIRE: 50 16:07:23 -1777.391713* 0.0061 WARNING: Max number of steps exceeded. Structure may be unstable. Relaxation Completed. Steps: 50 Cell Size Relaxation... Optimization terminated successfully. Current function value: 2.817045 Iterations: 236 Function evaluations: 487 Current VFE: 2.8170450755130787 Energy of Supercell: -1784.3417625149036 Unrelaxed Cell Volume: 5099.977952334259 Current Relaxed Cell Volume: 5095.007756954457 Current Relaxation Volume: 4.9701953798012255 Current Cell: [[1.72073837e+01 0.00000000e+00 0.00000000e+00] [4.03685876e-05 1.72073873e+01 0.00000000e+00] [3.55044623e-05 7.88165999e-05 1.72073928e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:07:43 -1777.394297* 0.0054 FIRE: 1 16:07:43 -1777.394312* 0.0047 FIRE: 2 16:07:43 -1777.394336* 0.0036 FIRE: 3 16:07:43 -1777.394355* 0.0020 FIRE: 4 16:07:43 -1777.394362* 0.0011 FIRE: 5 16:07:43 -1777.394362* 0.0010 FIRE: 6 16:07:43 -1777.394363* 0.0010 Relaxation Completed. Steps: 6 Cell Size Relaxation... Optimization terminated successfully. Current function value: 2.816979 Iterations: 272 Function evaluations: 513 Current VFE: 2.816978846826032 Energy of Supercell: -1784.3417625149036 Unrelaxed Cell Volume: 5099.977952334259 Current Relaxed Cell Volume: 5095.010659736326 Current Relaxation Volume: 4.96729259793301 Current Cell: [[1.72073914e+01 0.00000000e+00 0.00000000e+00] [5.69169549e-08 1.72073912e+01 0.00000000e+00] [1.63189893e-07 7.72020532e-07 1.72073910e+01]] ========== Loop: 2 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:08:14 -1777.394363* 0.0010 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 2.816979 Iterations: 115 Function evaluations: 286 Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:08:32 -1777.394363* 0.0010 FIRE: 1 16:08:32 -1777.394364* 0.0009 FIRE: 2 16:08:32 -1777.394365* 0.0009 FIRE: 3 16:08:32 -1777.394366* 0.0008 FIRE: 4 16:08:32 -1777.394368* 0.0008 FIRE: 5 16:08:32 -1777.394369* 0.0009 FIRE: 6 16:08:32 -1777.394371* 0.0008 FIRE: 7 16:08:32 -1777.394372* 0.0006 FIRE: 8 16:08:32 -1777.394374* 0.0006 FIRE: 9 16:08:32 -1777.394375* 0.0004 FIRE: 10 16:08:32 -1777.394376* 0.0002 FIRE: 11 16:08:32 -1777.394375* 0.0004 FIRE: 12 16:08:32 -1777.394375* 0.0003 FIRE: 13 16:08:32 -1777.394375* 0.0003 FIRE: 14 16:08:32 -1777.394376* 0.0002 FIRE: 15 16:08:32 -1777.394376* 0.0002 FIRE: 16 16:08:32 -1777.394376* 0.0002 FIRE: 17 16:08:32 -1777.394376* 0.0002 FIRE: 18 16:08:32 -1777.394376* 0.0001 FIRE: 19 16:08:32 -1777.394376* 0.0001 FIRE: 20 16:08:32 -1777.394376* 0.0001 Optimization terminated successfully. Current function value: 2.816966 Iterations: 171 Function evaluations: 416 --------------- Calculation Completed. Number Of Atoms in Supercell: 432 Vacancy Formation Energy (relaxed): 2.8169661184920187 Vacancy Formation Energy (unrelaxed): 4.1304207465607305 Unrelaxed Cell Volume: 5099.977952334259 Relaxed Cell Volume: 5095.010659736326 Relaxation Volume: 4.96729259793301 Relaxed Cell Vector: [17.207394180071333, 5.6515813885511204e-08, 17.207391651845697, 1.6568255583791113e-07, 7.835216643486153e-07, 17.207391202961976] Unrelaxed Cell Vector: [17.212981402873993, 0.0, 17.212981402873993, 0.0, 0.0, 17.212981402873993] Relaxed Cell: [[1.72073942e+01 0.00000000e+00 0.00000000e+00] [5.65158139e-08 1.72073917e+01 0.00000000e+00] [1.65682556e-07 7.83521664e-07 1.72073912e+01]] Unrelaxed Cell: [[17.2129814 0. 0. ] [ 0. 17.2129814 0. ] [ 0. 0. 17.2129814]] Supercell Size: 7 Unrelaxed Cell: [[20.08181164 0. 0. ] [ 0. 20.08181164 0. ] [ 0. 0. 20.08181164]] Unrelaxed Cell Vector: [20.081811636686325, 0.0, 20.081811636686325, 0.0, 0.0, 20.081811636686325] Unrelaxed Cell Energy: -2833.4686321412732 Energy of Unrelaxed Cell With Vacancy: -2833.4686321412732 ========== Loop: 0 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:08:50 -2825.207791* 1.2423 FIRE: 1 16:08:50 -2825.328988* 1.1061 FIRE: 2 16:08:50 -2825.519032* 0.8642 FIRE: 3 16:08:50 -2825.714549* 0.5686 FIRE: 4 16:08:50 -2825.884363* 0.4961 FIRE: 5 16:08:50 -2826.031206* 0.4803 FIRE: 6 16:08:50 -2826.157697* 0.3345 FIRE: 7 16:08:50 -2826.252079* 0.2166 FIRE: 8 16:08:50 -2826.310689* 0.2110 FIRE: 9 16:08:50 -2826.334002* 0.3456 FIRE: 10 16:08:50 -2826.343353* 0.3230 FIRE: 11 16:08:50 -2826.360418* 0.2798 FIRE: 12 16:08:50 -2826.382358* 0.2200 FIRE: 13 16:08:50 -2826.405831* 0.1609 FIRE: 14 16:08:50 -2826.427805* 0.1399 FIRE: 15 16:08:50 -2826.446259* 0.1113 FIRE: 16 16:08:50 -2826.460471* 0.0742 FIRE: 17 16:08:50 -2826.471652* 0.0823 FIRE: 18 16:08:50 -2826.479621* 0.0928 FIRE: 19 16:08:50 -2826.484740* 0.0993 FIRE: 20 16:08:50 -2826.488340* 0.1543 FIRE: 21 16:08:50 -2826.492969* 0.1862 FIRE: 22 16:08:50 -2826.501032* 0.1823 FIRE: 23 16:08:50 -2826.512025* 0.1353 FIRE: 24 16:08:50 -2826.521561* 0.0517 FIRE: 25 16:08:50 -2826.525409* 0.0550 FIRE: 26 16:08:50 -2826.526351* 0.0506 FIRE: 27 16:08:50 -2826.527974* 0.0423 FIRE: 28 16:08:50 -2826.529858* 0.0311 FIRE: 29 16:08:50 -2826.531593* 0.0183 FIRE: 30 16:08:50 -2826.532928* 0.0224 FIRE: 31 16:08:50 -2826.533830* 0.0331 FIRE: 32 16:08:50 -2826.534398* 0.0378 FIRE: 33 16:08:50 -2826.534736* 0.0353 FIRE: 34 16:08:50 -2826.534791* 0.0250 FIRE: 35 16:08:50 -2826.534831* 0.0244 FIRE: 36 16:08:50 -2826.534908* 0.0231 FIRE: 37 16:08:50 -2826.535014* 0.0212 FIRE: 38 16:08:50 -2826.535140* 0.0187 FIRE: 39 16:08:50 -2826.535275* 0.0158 FIRE: 40 16:08:50 -2826.535406* 0.0125 FIRE: 41 16:08:51 -2826.535523* 0.0089 FIRE: 42 16:08:51 -2826.535625* 0.0063 FIRE: 43 16:08:51 -2826.535695* 0.0041 FIRE: 44 16:08:51 -2826.535717* 0.0042 FIRE: 45 16:08:51 -2826.535718* 0.0042 FIRE: 46 16:08:51 -2826.535719* 0.0040 FIRE: 47 16:08:51 -2826.535721* 0.0039 FIRE: 48 16:08:51 -2826.535724* 0.0037 FIRE: 49 16:08:51 -2826.535727* 0.0034 FIRE: 50 16:08:51 -2826.535730* 0.0031 WARNING: Max number of steps exceeded. Structure may be unstable. Relaxation Completed. Steps: 50 Cell Size Relaxation... Optimization terminated successfully. Current function value: 2.800832 Iterations: 375 Function evaluations: 666 Current VFE: 2.8008324623401677 Energy of Supercell: -2833.4686321412732 Unrelaxed Cell Volume: 8098.576100234481 Current Relaxed Cell Volume: 8093.576659751811 Current Relaxation Volume: 4.999440482669343 Current Cell: [[ 2.00776782e+01 0.00000000e+00 0.00000000e+00] [-4.80296201e-07 2.00776769e+01 0.00000000e+00] [-1.07131915e-06 1.09564506e-08 2.00776803e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:09:16 -2826.537379* 0.0024 FIRE: 1 16:09:16 -2826.537384* 0.0021 FIRE: 2 16:09:16 -2826.537393* 0.0019 FIRE: 3 16:09:16 -2826.537404* 0.0018 FIRE: 4 16:09:16 -2826.537415* 0.0016 FIRE: 5 16:09:16 -2826.537425* 0.0013 FIRE: 6 16:09:16 -2826.537435* 0.0014 FIRE: 7 16:09:16 -2826.537444* 0.0013 FIRE: 8 16:09:16 -2826.537454* 0.0008 Relaxation Completed. Steps: 8 Cell Size Relaxation... Optimization terminated successfully. Current function value: 2.800758 Iterations: 125 Function evaluations: 307 Current VFE: 2.8007578710071357 Energy of Supercell: -2833.4686321412732 Unrelaxed Cell Volume: 8098.576100234481 Current Relaxed Cell Volume: 8093.576443731238 Current Relaxation Volume: 4.999656503242477 Current Cell: [[ 2.00776778e+01 0.00000000e+00 0.00000000e+00] [-4.83960154e-07 2.00776778e+01 0.00000000e+00] [-1.10756129e-06 1.10620552e-08 2.00776793e+01]] ========== Loop: 2 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:09:40 -2826.537454* 0.0008 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 2.800758 Iterations: 104 Function evaluations: 269 Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:10:02 -2826.537454* 0.0008 FIRE: 1 16:10:02 -2826.537454* 0.0008 FIRE: 2 16:10:02 -2826.537456* 0.0006 FIRE: 3 16:10:02 -2826.537457* 0.0005 FIRE: 4 16:10:02 -2826.537459* 0.0005 FIRE: 5 16:10:02 -2826.537460* 0.0005 FIRE: 6 16:10:02 -2826.537462* 0.0005 FIRE: 7 16:10:02 -2826.537463* 0.0005 FIRE: 8 16:10:03 -2826.537465* 0.0004 FIRE: 9 16:10:03 -2826.537466* 0.0003 FIRE: 10 16:10:03 -2826.537467* 0.0002 FIRE: 11 16:10:03 -2826.537468* 0.0004 FIRE: 12 16:10:03 -2826.537468* 0.0004 FIRE: 13 16:10:03 -2826.537468* 0.0003 FIRE: 14 16:10:03 -2826.537468* 0.0003 FIRE: 15 16:10:03 -2826.537468* 0.0002 FIRE: 16 16:10:03 -2826.537468* 0.0001 FIRE: 17 16:10:03 -2826.537468* 0.0001 FIRE: 18 16:10:03 -2826.537468* 0.0001 FIRE: 19 16:10:03 -2826.537468* 0.0001 FIRE: 20 16:10:03 -2826.537468* 0.0001 Optimization terminated successfully. Current function value: 2.800743 Iterations: 171 Function evaluations: 421 --------------- Calculation Completed. Number Of Atoms in Supercell: 686 Vacancy Formation Energy (relaxed): 2.8007433558491357 Vacancy Formation Energy (unrelaxed): 4.130420746553227 Unrelaxed Cell Volume: 8098.576100234481 Relaxed Cell Volume: 8093.576443731238 Relaxation Volume: 4.999656503242477 Relaxed Cell Vector: [20.077676438592206, -4.940424556285391e-07, 20.077677531181724, -1.0937177253849528e-06, 1.1292512685985038e-08, 20.077676414945394] Unrelaxed Cell Vector: [20.081811636686325, 0.0, 20.081811636686325, 0.0, 0.0, 20.081811636686325] Relaxed Cell: [[ 2.00776764e+01 0.00000000e+00 0.00000000e+00] [-4.94042456e-07 2.00776775e+01 0.00000000e+00] [-1.09371773e-06 1.12925127e-08 2.00776764e+01]] Unrelaxed Cell: [[20.08181164 0. 0. ] [ 0. 20.08181164 0. ] [ 0. 0. 20.08181164]] [Calculation Results Summary] Sizes: [5, 6, 7] Unrelaxed Formation Energy By Size: [4.13042074655209, 4.1304207465607305, 4.130420746553227] Formation Energy By Size: [2.848554117687854, 2.8169661184920187, 2.8007433558491357] Relaxation Volume By Size: [4.894594189252075, 4.96729259793301, 4.999656503242477] [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: [4.13042075 4.13042075] Fitting Results: (array([ 4.13042075e+00, -2.56367457e-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: [2.84855412 2.81696612] Fitting Results: (array([2.77357601, 9.3722635 ]), 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: [4.89459419 4.9672926 ] Fitting Results: (array([ 5.06715305, -21.56985752]), 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: [4.13042075 4.13042075] Fitting Results: (array([4.13042075e+00, 4.37748416e-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: [2.81696612 2.80074336] Fitting Results: (array([2.77315189, 9.46387432]), 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: [4.9672926 4.9996565] Fitting Results: (array([ 5.05470063, -18.88013399]), 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: [4.13042075 4.13042075 4.13042075] Fitting Results: (array([ 4.13042075e+00, -5.48998502e-10]), array([4.00507105e-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: [2.84855412 2.81696612 2.80074336] Fitting Results: (array([2.77338588, 9.39885482]), array([6.97702548e-09]), 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: [4.89459419 4.9672926 4.9996565 ] Fitting Results: (array([ 5.06157066, -20.78912763]), array([6.01439611e-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: [4.13042075 4.13042075 4.13042075] Fitting Results: (array([ 4.13042075e+00, 4.16032672e-08, -1.79695693e-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: [2.84855412 2.81696612 2.80074336] Fitting Results: (array([ 2.77270724, 9.9552085 , -2.37174345]), 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: [4.89459419 4.9672926 4.9996565 ] Fitting Results: (array([ 5.04164567, -4.45440275, -69.63515873]), 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: [4.13042075 4.13042075 4.13042075] Fitting Results: (array([ 4.13042075e+00, 2.13587295e-08, -4.21252778e-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: [2.84855412 2.81696612 2.80074336] Fitting Results: (array([ 2.77282925, 9.68800761, -5.55997475]), 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: [4.89459419 4.9672926 4.9996565 ] Fitting Results: (array([ 5.04522791, -12.29950783, -163.24266596]), 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: [4.13042075 4.13042075 4.13042075] Fitting Results: (array([ 4.13042075e+00, 1.46770077e-08, -1.36509973e-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: [2.84855412 2.81696612 2.80074336] Fitting Results: (array([ 2.7729087 , 9.59981779, -18.01749552]), 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: [4.89459419 4.9672926 4.9996565 ] Fitting Results: (array([ 5.04756047, -14.88878945, -528.99952493]), 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: [[4.1304207465726, 4.130420746540466], [4.130420746558196], [4.130420746506776], [4.130420746516021], [4.130420746522042]] Formation Energy Fits By Size: [[2.773576009706531, 2.773151885527383], [2.7733858758668806], [2.7727072392105327], [2.7728292487508424], [2.7729086950576876]] Relaxation Volume Fits By Size: [[5.067153049417813, 5.054700625658576], [5.061570658083615], [5.041645665265833], [5.045227905027905], [5.047560474420488]] [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" 4.130420746540466 "source-unit" "eV" "source-std-uncert-value" 1.4515158000033795e-05 } "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" 2.868830233812332 "source-unit" "angstrom" } "host-b" { "source-value" 2.868830233812332 "source-unit" "angstrom" } "host-c" { "source-value" 2.868830233812332 "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" [ "Fe" ] } "reservoir-cohesive-potential-energy" { "source-value" 4.130420746562632 "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" 2.868830233812332 "source-unit" "angstrom" } "reservoir-b" { "source-value" 2.868830233812332 "source-unit" "angstrom" } "reservoir-c" { "source-value" 2.868830233812332 "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" [ "Fe" ] } } { "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" 2.773151885527383 "source-unit" "eV" "source-std-uncert-value" 0.00044488317219262065 } "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" 2.868830233812332 "source-unit" "angstrom" } "host-b" { "source-value" 2.868830233812332 "source-unit" "angstrom" } "host-c" { "source-value" 2.868830233812332 "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" [ "Fe" ] } "reservoir-cohesive-potential-energy" { "source-value" 4.130420746562632 "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" 2.868830233812332 "source-unit" "angstrom" } "reservoir-b" { "source-value" 2.868830233812332 "source-unit" "angstrom" } "reservoir-c" { "source-value" 2.868830233812332 "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" [ "Fe" ] } } { "property-id" "tag:staff@noreply.openkim.org,2015-07-28:property/monovacancy-neutral-relaxation-volume-crystal-npt" "instance-id" 3 "relaxation-volume" { "source-value" 5.054700625658576 "source-unit" "angstrom^3" "source-std-uncert-value" 0.013178480356131552 } "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" 2.868830233812332 "source-unit" "angstrom" } "host-b" { "source-value" 2.868830233812332 "source-unit" "angstrom" } "host-c" { "source-value" 2.868830233812332 "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" [ "Fe" ] } } ]