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 Sim_LAMMPS_TersoffZBL_HenrikssonBjorkasNordlund_2013_FeC__SM_473463498269_000 [2.8886456936597833] CELL_SIZE_MIN: 5 CELL_SIZE_MAX: 7 Smallest System Size: 250 Largest System Size: 686 [Calculation] Supercell Size: 5 Unrelaxed Cell: [[14.44322847 0. 0. ] [ 0. 14.44322847 0. ] [ 0. 0. 14.44322847]] Unrelaxed Cell Vector: [14.443228468298916, 0.0, 14.443228468298916, 0.0, 0.0, 14.443228468298916] Unrelaxed Cell Energy: -1044.6645492346443 Energy of Unrelaxed Cell With Vacancy: -1044.6645492346443 ========== Loop: 0 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:06:42 -1038.733055* 0.3406 FIRE: 1 16:06:42 -1038.748838* 0.3153 FIRE: 2 16:06:42 -1038.776475* 0.2857 FIRE: 3 16:06:42 -1038.809625* 0.2438 FIRE: 4 16:06:42 -1038.841926* 0.1920 FIRE: 5 16:06:42 -1038.869057* 0.1334 FIRE: 6 16:06:42 -1038.889290* 0.0918 FIRE: 7 16:06:42 -1038.902275* 0.0744 FIRE: 8 16:06:42 -1038.908047* 0.0641 FIRE: 9 16:06:43 -1038.908453* 0.0622 FIRE: 10 16:06:43 -1038.909227* 0.0583 FIRE: 11 16:06:43 -1038.910296* 0.0527 FIRE: 12 16:06:43 -1038.911561* 0.0456 FIRE: 13 16:06:43 -1038.912910* 0.0372 FIRE: 14 16:06:43 -1038.914228* 0.0291 FIRE: 15 16:06:43 -1038.915413* 0.0241 FIRE: 16 16:06:43 -1038.916475* 0.0182 FIRE: 17 16:06:43 -1038.917281* 0.0130 FIRE: 18 16:06:43 -1038.917732* 0.0146 FIRE: 19 16:06:43 -1038.917823* 0.0195 FIRE: 20 16:06:43 -1038.917843* 0.0192 FIRE: 21 16:06:43 -1038.917882* 0.0188 FIRE: 22 16:06:43 -1038.917937* 0.0182 FIRE: 23 16:06:43 -1038.918005* 0.0175 FIRE: 24 16:06:43 -1038.918082* 0.0165 FIRE: 25 16:06:43 -1038.918164* 0.0154 FIRE: 26 16:06:43 -1038.918246* 0.0142 FIRE: 27 16:06:43 -1038.918332* 0.0126 FIRE: 28 16:06:43 -1038.918414* 0.0108 FIRE: 29 16:06:43 -1038.918486* 0.0085 FIRE: 30 16:06:43 -1038.918539* 0.0060 FIRE: 31 16:06:44 -1038.918571* 0.0062 FIRE: 32 16:06:44 -1038.918584* 0.0065 FIRE: 33 16:06:44 -1038.918586* 0.0064 FIRE: 34 16:06:44 -1038.918589* 0.0063 FIRE: 35 16:06:44 -1038.918595* 0.0061 FIRE: 36 16:06:44 -1038.918602* 0.0059 FIRE: 37 16:06:44 -1038.918610* 0.0056 FIRE: 38 16:06:44 -1038.918619* 0.0053 FIRE: 39 16:06:44 -1038.918628* 0.0049 FIRE: 40 16:06:44 -1038.918638* 0.0045 FIRE: 41 16:06:44 -1038.918648* 0.0039 FIRE: 42 16:06:44 -1038.918658* 0.0033 FIRE: 43 16:06:44 -1038.918667* 0.0026 FIRE: 44 16:06:44 -1038.918674* 0.0020 FIRE: 45 16:06:44 -1038.918679* 0.0022 FIRE: 46 16:06:44 -1038.918684* 0.0027 FIRE: 47 16:06:44 -1038.918688* 0.0029 FIRE: 48 16:06:44 -1038.918692* 0.0026 FIRE: 49 16:06:44 -1038.918697* 0.0019 FIRE: 50 16:06:44 -1038.918700* 0.0011 WARNING: Max number of steps exceeded. Structure may be unstable. Relaxation Completed. Steps: 50 Cell Size Relaxation... Optimization terminated successfully. Current function value: 1.563970 Iterations: 408 Function evaluations: 726 Current VFE: 1.5639697391968639 Energy of Supercell: -1044.6645492346443 Unrelaxed Cell Volume: 3012.956374200818 Current Relaxed Cell Volume: 3008.6376617505525 Current Relaxation Volume: 4.318712450265593 Current Cell: [[1.44363241e+01 0.00000000e+00 0.00000000e+00] [6.64692972e-06 1.44363238e+01 0.00000000e+00] [2.47334851e-05 1.54459160e-05 1.44363250e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:07:07 -1038.921921* 0.0043 FIRE: 1 16:07:07 -1038.921923* 0.0039 FIRE: 2 16:07:08 -1038.921927* 0.0032 FIRE: 3 16:07:08 -1038.921930* 0.0022 FIRE: 4 16:07:08 -1038.921932* 0.0010 FIRE: 5 16:07:08 -1038.921933* 0.0011 FIRE: 6 16:07:08 -1038.921933* 0.0011 FIRE: 7 16:07:08 -1038.921933* 0.0010 Relaxation Completed. Steps: 7 Cell Size Relaxation... Optimization terminated successfully. Current function value: 1.563958 Iterations: 215 Function evaluations: 433 Current VFE: 1.5639577439524146 Energy of Supercell: -1044.6645492346443 Unrelaxed Cell Volume: 3012.956374200818 Current Relaxed Cell Volume: 3008.6280712098355 Current Relaxation Volume: 4.328302990982593 Current Cell: [[1.44363093e+01 0.00000000e+00 0.00000000e+00] [7.79890966e-06 1.44363088e+01 0.00000000e+00] [6.27767725e-06 2.09915871e-05 1.44363087e+01]] ========== Loop: 2 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:07:26 -1038.921933* 0.0010 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 1.563958 Iterations: 256 Function evaluations: 480 Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:07:44 -1038.921933* 0.0010 FIRE: 1 16:07:44 -1038.921933* 0.0009 FIRE: 2 16:07:44 -1038.921934* 0.0006 FIRE: 3 16:07:44 -1038.921934* 0.0004 FIRE: 4 16:07:44 -1038.921934* 0.0004 FIRE: 5 16:07:44 -1038.921934* 0.0005 FIRE: 6 16:07:45 -1038.921934* 0.0005 FIRE: 7 16:07:45 -1038.921934* 0.0004 FIRE: 8 16:07:45 -1038.921934* 0.0002 FIRE: 9 16:07:45 -1038.921934* 0.0001 FIRE: 10 16:07:45 -1038.921934* 0.0001 FIRE: 11 16:07:45 -1038.921934* 0.0001 FIRE: 12 16:07:45 -1038.921934* 0.0001 FIRE: 13 16:07:45 -1038.921934* 0.0001 FIRE: 14 16:07:45 -1038.921934* 0.0001 FIRE: 15 16:07:45 -1038.921934* 0.0001 FIRE: 16 16:07:45 -1038.921934* 0.0001 FIRE: 17 16:07:45 -1038.921934* 0.0001 FIRE: 18 16:07:45 -1038.921934* 0.0001 FIRE: 19 16:07:45 -1038.921934* 0.0000 FIRE: 20 16:07:45 -1038.921934* 0.0000 Optimization terminated successfully. Current function value: 1.563957 Iterations: 311 Function evaluations: 616 --------------- Calculation Completed. Number Of Atoms in Supercell: 250 Vacancy Formation Energy (relaxed): 1.5639565902017694 Vacancy Formation Energy (unrelaxed): 1.7528356744508073 Unrelaxed Cell Volume: 3012.956374200818 Relaxed Cell Volume: 3008.6280712098355 Relaxation Volume: 4.328302990982593 Relaxed Cell Vector: [14.436305443761889, -5.923019995741986e-08, 14.436305928226822, 7.720407928981705e-06, 2.76530000128849e-07, 14.436307228347387] Unrelaxed Cell Vector: [14.443228468298916, 0.0, 14.443228468298916, 0.0, 0.0, 14.443228468298916] Relaxed Cell: [[ 1.44363054e+01 0.00000000e+00 0.00000000e+00] [-5.92302000e-08 1.44363059e+01 0.00000000e+00] [ 7.72040793e-06 2.76530000e-07 1.44363072e+01]] Unrelaxed Cell: [[14.44322847 0. 0. ] [ 0. 14.44322847 0. ] [ 0. 0. 14.44322847]] Supercell Size: 6 Unrelaxed Cell: [[17.33187416 0. 0. ] [ 0. 17.33187416 0. ] [ 0. 0. 17.33187416]] Unrelaxed Cell Vector: [17.3318741619587, 0.0, 17.3318741619587, 0.0, 0.0, 17.3318741619587] Unrelaxed Cell Energy: -1805.1803410774642 Energy of Unrelaxed Cell With Vacancy: -1805.1803410774642 ========== Loop: 0 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:08:18 -1799.248847* 0.3406 FIRE: 1 16:08:19 -1799.264630* 0.3153 FIRE: 2 16:08:19 -1799.292267* 0.2857 FIRE: 3 16:08:19 -1799.325420* 0.2438 FIRE: 4 16:08:19 -1799.357722* 0.1920 FIRE: 5 16:08:19 -1799.384855* 0.1334 FIRE: 6 16:08:19 -1799.405090* 0.0918 FIRE: 7 16:08:19 -1799.418087* 0.0745 FIRE: 8 16:08:19 -1799.423885* 0.0642 FIRE: 9 16:08:19 -1799.424292* 0.0622 FIRE: 10 16:08:19 -1799.425070* 0.0584 FIRE: 11 16:08:19 -1799.426144* 0.0527 FIRE: 12 16:08:19 -1799.427419* 0.0456 FIRE: 13 16:08:19 -1799.428781* 0.0372 FIRE: 14 16:08:19 -1799.430120* 0.0293 FIRE: 15 16:08:19 -1799.431335* 0.0244 FIRE: 16 16:08:20 -1799.432446* 0.0188 FIRE: 17 16:08:20 -1799.433334* 0.0131 FIRE: 18 16:08:20 -1799.433920* 0.0138 FIRE: 19 16:08:20 -1799.434217* 0.0196 FIRE: 20 16:08:20 -1799.434322* 0.0259 FIRE: 21 16:08:20 -1799.434359* 0.0256 FIRE: 22 16:08:20 -1799.434430* 0.0249 FIRE: 23 16:08:20 -1799.434531* 0.0239 FIRE: 24 16:08:20 -1799.434655* 0.0225 FIRE: 25 16:08:20 -1799.434794* 0.0209 FIRE: 26 16:08:20 -1799.434939* 0.0190 FIRE: 27 16:08:20 -1799.435082* 0.0169 FIRE: 28 16:08:20 -1799.435226* 0.0143 FIRE: 29 16:08:21 -1799.435358* 0.0113 FIRE: 30 16:08:21 -1799.435464* 0.0078 FIRE: 31 16:08:21 -1799.435531* 0.0061 FIRE: 32 16:08:21 -1799.435559* 0.0071 FIRE: 33 16:08:21 -1799.435561* 0.0071 FIRE: 34 16:08:21 -1799.435566* 0.0069 FIRE: 35 16:08:21 -1799.435573* 0.0068 FIRE: 36 16:08:21 -1799.435582* 0.0066 FIRE: 37 16:08:21 -1799.435593* 0.0063 FIRE: 38 16:08:21 -1799.435605* 0.0060 FIRE: 39 16:08:22 -1799.435617* 0.0056 FIRE: 40 16:08:22 -1799.435632* 0.0051 FIRE: 41 16:08:22 -1799.435648* 0.0046 FIRE: 42 16:08:22 -1799.435665* 0.0040 FIRE: 43 16:08:22 -1799.435681* 0.0032 FIRE: 44 16:08:22 -1799.435698* 0.0025 FIRE: 45 16:08:22 -1799.435713* 0.0026 FIRE: 46 16:08:22 -1799.435727* 0.0028 FIRE: 47 16:08:22 -1799.435741* 0.0029 FIRE: 48 16:08:22 -1799.435756* 0.0027 FIRE: 49 16:08:23 -1799.435771* 0.0022 FIRE: 50 16:08:23 -1799.435783* 0.0013 WARNING: Max number of steps exceeded. Structure may be unstable. Relaxation Completed. Steps: 50 Cell Size Relaxation... Optimization terminated successfully. Current function value: 1.564031 Iterations: 569 Function evaluations: 952 Current VFE: 1.5640310621104163 Energy of Supercell: -1805.1803410774642 Unrelaxed Cell Volume: 5206.388614619008 Current Relaxed Cell Volume: 5202.069729448395 Current Relaxation Volume: 4.3188851706127025 Current Cell: [[ 1.73270804e+01 0.00000000e+00 0.00000000e+00] [ 4.98721768e-07 1.73270805e+01 0.00000000e+00] [-7.26180958e-07 5.18780758e-07 1.73270803e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:09:17 -1799.437652* 0.0033 FIRE: 1 16:09:17 -1799.437654* 0.0031 FIRE: 2 16:09:17 -1799.437656* 0.0025 FIRE: 3 16:09:17 -1799.437659* 0.0018 FIRE: 4 16:09:17 -1799.437661* 0.0009 Relaxation Completed. Steps: 4 Cell Size Relaxation... Optimization terminated successfully. Current function value: 1.564021 Iterations: 136 Function evaluations: 316 Current VFE: 1.5640214535180803 Energy of Supercell: -1805.1803410774642 Unrelaxed Cell Volume: 5206.388614619008 Current Relaxed Cell Volume: 5202.0644862601 Current Relaxation Volume: 4.324128358907728 Current Cell: [[ 1.73270744e+01 0.00000000e+00 0.00000000e+00] [ 4.97494474e-07 1.73270746e+01 0.00000000e+00] [-7.23661934e-07 5.34238798e-07 1.73270746e+01]] ========== Loop: 2 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:09:34 -1799.437661* 0.0009 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 1.564021 Iterations: 117 Function evaluations: 287 Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:09:52 -1799.437661* 0.0009 FIRE: 1 16:09:52 -1799.437662* 0.0009 FIRE: 2 16:09:52 -1799.437662* 0.0008 FIRE: 3 16:09:52 -1799.437663* 0.0007 FIRE: 4 16:09:52 -1799.437663* 0.0006 FIRE: 5 16:09:52 -1799.437664* 0.0004 FIRE: 6 16:09:52 -1799.437664* 0.0003 FIRE: 7 16:09:52 -1799.437664* 0.0003 FIRE: 8 16:09:52 -1799.437665* 0.0002 FIRE: 9 16:09:52 -1799.437665* 0.0002 FIRE: 10 16:09:52 -1799.437665* 0.0003 FIRE: 11 16:09:52 -1799.437665* 0.0002 FIRE: 12 16:09:52 -1799.437665* 0.0002 FIRE: 13 16:09:52 -1799.437665* 0.0002 FIRE: 14 16:09:53 -1799.437665* 0.0002 FIRE: 15 16:09:53 -1799.437665* 0.0001 FIRE: 16 16:09:53 -1799.437665* 0.0001 FIRE: 17 16:09:53 -1799.437665* 0.0001 FIRE: 18 16:09:53 -1799.437665* 0.0001 FIRE: 19 16:09:53 -1799.437665* 0.0001 FIRE: 20 16:09:53 -1799.437665* 0.0000 Optimization terminated successfully. Current function value: 1.564018 Iterations: 178 Function evaluations: 426 --------------- Calculation Completed. Number Of Atoms in Supercell: 432 Vacancy Formation Energy (relaxed): 1.5640176354916093 Vacancy Formation Energy (unrelaxed): 1.7528356744464872 Unrelaxed Cell Volume: 5206.388614619008 Relaxed Cell Volume: 5202.0644862601 Relaxation Volume: 4.324128358907728 Relaxed Cell Vector: [17.327069114088644, 5.121404549295904e-07, 17.327069335252574, -7.461159352202269e-07, 5.303456570222888e-07, 17.32706904468246] Unrelaxed Cell Vector: [17.3318741619587, 0.0, 17.3318741619587, 0.0, 0.0, 17.3318741619587] Relaxed Cell: [[ 1.73270691e+01 0.00000000e+00 0.00000000e+00] [ 5.12140455e-07 1.73270693e+01 0.00000000e+00] [-7.46115935e-07 5.30345657e-07 1.73270690e+01]] Unrelaxed Cell: [[17.33187416 0. 0. ] [ 0. 17.33187416 0. ] [ 0. 0. 17.33187416]] Supercell Size: 7 Unrelaxed Cell: [[20.22051986 0. 0. ] [ 0. 20.22051986 0. ] [ 0. 0. 20.22051986]] Unrelaxed Cell Vector: [20.220519855618484, 0.0, 20.220519855618484, 0.0, 0.0, 20.220519855618484] Unrelaxed Cell Energy: -2866.559523099658 Energy of Unrelaxed Cell With Vacancy: -2866.559523099658 ========== Loop: 0 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:10:20 -2860.628029* 0.3406 FIRE: 1 16:10:20 -2860.643812* 0.3153 FIRE: 2 16:10:20 -2860.671449* 0.2857 FIRE: 3 16:10:20 -2860.704602* 0.2438 FIRE: 4 16:10:20 -2860.736904* 0.1920 FIRE: 5 16:10:20 -2860.764037* 0.1334 FIRE: 6 16:10:20 -2860.784274* 0.0918 FIRE: 7 16:10:21 -2860.797272* 0.0745 FIRE: 8 16:10:21 -2860.803072* 0.0642 FIRE: 9 16:10:21 -2860.803481* 0.0623 FIRE: 10 16:10:21 -2860.804258* 0.0584 FIRE: 11 16:10:21 -2860.805334* 0.0528 FIRE: 12 16:10:21 -2860.806609* 0.0456 FIRE: 13 16:10:21 -2860.807973* 0.0372 FIRE: 14 16:10:21 -2860.809314* 0.0293 FIRE: 15 16:10:21 -2860.810532* 0.0244 FIRE: 16 16:10:21 -2860.811649* 0.0188 FIRE: 17 16:10:21 -2860.812545* 0.0131 FIRE: 18 16:10:22 -2860.813147* 0.0137 FIRE: 19 16:10:22 -2860.813473* 0.0196 FIRE: 20 16:10:22 -2860.813629* 0.0260 FIRE: 21 16:10:22 -2860.813668* 0.0256 FIRE: 22 16:10:22 -2860.813745* 0.0249 FIRE: 23 16:10:22 -2860.813853* 0.0239 FIRE: 24 16:10:22 -2860.813986* 0.0226 FIRE: 25 16:10:22 -2860.814137* 0.0210 FIRE: 26 16:10:22 -2860.814296* 0.0191 FIRE: 27 16:10:22 -2860.814454* 0.0170 FIRE: 28 16:10:22 -2860.814618* 0.0144 FIRE: 29 16:10:22 -2860.814776* 0.0114 FIRE: 30 16:10:23 -2860.814912* 0.0079 FIRE: 31 16:10:23 -2860.815016* 0.0062 FIRE: 32 16:10:23 -2860.815086* 0.0073 FIRE: 33 16:10:23 -2860.815133* 0.0078 FIRE: 34 16:10:23 -2860.815179* 0.0082 FIRE: 35 16:10:23 -2860.815242* 0.0079 FIRE: 36 16:10:23 -2860.815325* 0.0081 FIRE: 37 16:10:23 -2860.815412* 0.0071 FIRE: 38 16:10:23 -2860.815467* 0.0047 FIRE: 39 16:10:23 -2860.815457* 0.0042 FIRE: 40 16:10:23 -2860.815462* 0.0040 FIRE: 41 16:10:24 -2860.815471* 0.0037 FIRE: 42 16:10:24 -2860.815482* 0.0033 FIRE: 43 16:10:24 -2860.815495* 0.0027 FIRE: 44 16:10:24 -2860.815506* 0.0021 FIRE: 45 16:10:24 -2860.815515* 0.0018 FIRE: 46 16:10:24 -2860.815521* 0.0018 FIRE: 47 16:10:24 -2860.815525* 0.0019 FIRE: 48 16:10:24 -2860.815527* 0.0019 FIRE: 49 16:10:24 -2860.815530* 0.0016 FIRE: 50 16:10:24 -2860.815534* 0.0015 WARNING: Max number of steps exceeded. Structure may be unstable. Relaxation Completed. Steps: 50 Cell Size Relaxation... Optimization terminated successfully. Current function value: 1.564149 Iterations: 237 Function evaluations: 488 Current VFE: 1.5641490511993652 Energy of Supercell: -2866.559523099658 Unrelaxed Cell Volume: 8267.552290807052 Current Relaxed Cell Volume: 8263.226487862936 Current Relaxation Volume: 4.325802944116731 Current Cell: [[2.02169919e+01 0.00000000e+00 0.00000000e+00] [4.26402300e-05 2.02169926e+01 0.00000000e+00] [2.01722713e-05 7.47989438e-05 2.02169933e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:11:02 -2860.816716* 0.0016 FIRE: 1 16:11:02 -2860.816718* 0.0013 FIRE: 2 16:11:02 -2860.816721* 0.0011 FIRE: 3 16:11:02 -2860.816724* 0.0008 Relaxation Completed. Steps: 3 Cell Size Relaxation... Optimization terminated successfully. Current function value: 1.564141 Iterations: 335 Function evaluations: 604 Current VFE: 1.5641411021329077 Energy of Supercell: -2866.559523099658 Unrelaxed Cell Volume: 8267.552290807052 Current Relaxed Cell Volume: 8263.225269184184 Current Relaxation Volume: 4.327021622868415 Current Cell: [[ 2.02169918e+01 0.00000000e+00 0.00000000e+00] [ 2.27866553e-07 2.02169917e+01 0.00000000e+00] [-1.17495009e-07 2.72457637e-08 2.02169913e+01]] ========== Loop: 2 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:11:39 -2860.816724* 0.0008 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 1.564141 Iterations: 119 Function evaluations: 284 Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:11:59 -2860.816724* 0.0008 FIRE: 1 16:11:59 -2860.816724* 0.0008 FIRE: 2 16:11:59 -2860.816725* 0.0006 FIRE: 3 16:11:59 -2860.816726* 0.0005 FIRE: 4 16:11:59 -2860.816727* 0.0003 FIRE: 5 16:12:00 -2860.816727* 0.0003 FIRE: 6 16:12:00 -2860.816728* 0.0004 FIRE: 7 16:12:00 -2860.816728* 0.0004 FIRE: 8 16:12:00 -2860.816728* 0.0003 FIRE: 9 16:12:00 -2860.816728* 0.0003 FIRE: 10 16:12:00 -2860.816728* 0.0003 FIRE: 11 16:12:00 -2860.816728* 0.0003 FIRE: 12 16:12:00 -2860.816728* 0.0003 FIRE: 13 16:12:00 -2860.816728* 0.0002 FIRE: 14 16:12:00 -2860.816728* 0.0002 FIRE: 15 16:12:00 -2860.816728* 0.0001 FIRE: 16 16:12:01 -2860.816728* 0.0001 FIRE: 17 16:12:01 -2860.816728* 0.0001 FIRE: 18 16:12:01 -2860.816728* 0.0001 FIRE: 19 16:12:01 -2860.816728* 0.0002 FIRE: 20 16:12:01 -2860.816728* 0.0002 Optimization terminated successfully. Current function value: 1.564137 Iterations: 183 Function evaluations: 433 --------------- Calculation Completed. Number Of Atoms in Supercell: 686 Vacancy Formation Energy (relaxed): 1.5641365194765058 Vacancy Formation Energy (unrelaxed): 1.752835674446942 Unrelaxed Cell Volume: 8267.552290807052 Relaxed Cell Volume: 8263.225269184184 Relaxation Volume: 4.327021622868415 Relaxed Cell Vector: [20.21699058949076, 2.3313864063498304e-07, 20.216990213282404, -1.1787488255916275e-07, 2.7857823466735695e-08, 20.21699127676532] Unrelaxed Cell Vector: [20.220519855618484, 0.0, 20.220519855618484, 0.0, 0.0, 20.220519855618484] Relaxed Cell: [[ 2.02169906e+01 0.00000000e+00 0.00000000e+00] [ 2.33138641e-07 2.02169902e+01 0.00000000e+00] [-1.17874883e-07 2.78578235e-08 2.02169913e+01]] Unrelaxed Cell: [[20.22051986 0. 0. ] [ 0. 20.22051986 0. ] [ 0. 0. 20.22051986]] [Calculation Results Summary] Sizes: [5, 6, 7] Unrelaxed Formation Energy By Size: [1.7528356744508073, 1.7528356744464872, 1.752835674446942] Formation Energy By Size: [1.5639565902017694, 1.5640176354916093, 1.5641365194765058] Relaxation Volume By Size: [4.328302990982593, 4.324128358907728, 4.327021622868415] [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.75283567 1.75283567] Fitting Results: (array([1.75283567e+00, 1.28183729e-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: [1.56395659 1.56401764] Fitting Results: (array([ 1.56410149, -0.01811234]), 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.32830299 4.32412836] Fitting Results: (array([4.31839397, 1.2386271 ]), 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.75283567 1.75283567] Fitting Results: (array([ 1.75283567e+00, -2.65440847e-10]), 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.56401764 1.56413652] Fitting Results: (array([ 1.56433872, -0.06935336]), 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.32412836 4.32702162] Fitting Results: (array([ 4.33194245, -1.68784362]), 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.75283567 1.75283567 1.75283567] Fitting Results: (array([1.75283567e+00, 8.32669703e-10]), array([1.98954011e-24]), 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.56395659 1.56401764 1.56413652] Fitting Results: (array([ 1.56420784, -0.03298576]), array([2.18279168e-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.32830299 4.32412836 4.32702162] Fitting Results: (array([4.32446772, 0.38917802]), array([7.119756e-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: [1.75283567 1.75283567 1.75283567] Fitting Results: (array([ 1.75283567e+00, -8.56223027e-09, 4.00505889e-08]), 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.56395659 1.56401764 1.56413652] Fitting Results: (array([ 1.56458742, -0.34417313, 1.32659605]), 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.32830299 4.32412836 4.32702162] Fitting Results: (array([ 4.3461465 , -17.38331553, 75.76438655]), 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.75283567 1.75283567 1.75283567] Fitting Results: (array([ 1.75283567e+00, -4.05012627e-09, 9.38888490e-08]), 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.56395659 1.56401764 1.56413652] Fitting Results: (array([ 1.56451918, -0.19471866, 3.10988128]), 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.32830299 4.32412836 4.32702162] Fitting Results: (array([ 4.34224895, -8.84769093, 177.61114743]), 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.75283567 1.75283567 1.75283567] Fitting Results: (array([ 1.75283567e+00, -2.56090364e-09, 3.04253525e-07]), 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.56395659 1.56401764 1.56413652] Fitting Results: (array([ 1.56447474, -0.14539113, 10.07779253]), 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.32830299 4.32412836 4.32702162] Fitting Results: (array([ 4.33971107, -6.03050294, 575.56161596]), 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.752835674440553, 1.7528356744477158], [1.7528356744437643], [1.752835674455224], [1.7528356744531635], [1.7528356744518225]] Formation Energy Fits By Size: [[1.5641014889117193, 1.5643387158602666], [1.5642078371767176], [1.5645874215321307], [1.5645191774838916], [1.5644747404867279]] Relaxation Volume Fits By Size: [[4.318393974189509, 4.331942449762181], [4.324467722890231], [4.346146496733532], [4.342248951255618], [4.33971107107679]] [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.7528356744477158 "source-unit" "eV" "source-std-uncert-value" 4.582656401914176e-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" 2.8886456936597833 "source-unit" "angstrom" } "host-b" { "source-value" 2.8886456936597833 "source-unit" "angstrom" } "host-c" { "source-value" 2.8886456936597833 "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.1786581969384695 "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.8886456936597833 "source-unit" "angstrom" } "reservoir-b" { "source-value" 2.8886456936597833 "source-unit" "angstrom" } "reservoir-c" { "source-value" 2.8886456936597833 "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" 1.5643387158602666 "source-unit" "eV" "source-std-uncert-value" 0.0002487478883469394 } "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.8886456936597833 "source-unit" "angstrom" } "host-b" { "source-value" 2.8886456936597833 "source-unit" "angstrom" } "host-c" { "source-value" 2.8886456936597833 "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.1786581969384695 "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.8886456936597833 "source-unit" "angstrom" } "reservoir-b" { "source-value" 2.8886456936597833 "source-unit" "angstrom" } "reservoir-c" { "source-value" 2.8886456936597833 "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" 4.331942449762181 "source-unit" "angstrom^3" "source-std-uncert-value" 0.014248348645109602 } "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.8886456936597833 "source-unit" "angstrom" } "host-b" { "source-value" 2.8886456936597833 "source-unit" "angstrom" } "host-c" { "source-value" 2.8886456936597833 "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" ] } } ]