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 MEAM_LAMMPS_JelinekGrohHorstemeyer_2012_AlSiMgCuFe__MO_262519520678_002 [2.850998893380165] CELL_SIZE_MIN: 5 CELL_SIZE_MAX: 7 Smallest System Size: 250 Largest System Size: 686 [Calculation] Supercell Size: 5 Unrelaxed Cell: [[14.25499447 0. 0. ] [ 0. 14.25499447 0. ] [ 0. 0. 14.25499447]] Unrelaxed Cell Vector: [14.254994466900826, 0.0, 14.254994466900826, 0.0, 0.0, 14.254994466900826] Unrelaxed Cell Energy: -1069.9995330904171 Energy of Unrelaxed Cell With Vacancy: -1069.9995330904171 ========== Loop: 0 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 15:58:12 -1063.319628* 0.7301 FIRE: 1 15:58:12 -1063.370282* 0.6869 FIRE: 2 15:58:12 -1063.459168* 0.6026 FIRE: 3 15:58:12 -1063.565868* 0.4860 FIRE: 4 15:58:12 -1063.670198* 0.3515 FIRE: 5 15:58:12 -1063.759784* 0.2270 FIRE: 6 15:58:12 -1063.831818* 0.2526 FIRE: 7 15:58:12 -1063.888747* 0.2391 FIRE: 8 15:58:12 -1063.938316* 0.1961 FIRE: 9 15:58:12 -1063.979845* 0.1240 FIRE: 10 15:58:12 -1064.001538* 0.1250 FIRE: 11 15:58:12 -1064.004267* 0.1216 FIRE: 12 15:58:12 -1064.009402* 0.1148 FIRE: 13 15:58:12 -1064.016344* 0.1048 FIRE: 14 15:58:12 -1064.024304* 0.0918 FIRE: 15 15:58:12 -1064.032418* 0.0761 FIRE: 16 15:58:12 -1064.039867* 0.0581 FIRE: 17 15:58:12 -1064.045990* 0.0381 FIRE: 18 15:58:12 -1064.050722* 0.0209 FIRE: 19 15:58:12 -1064.053327* 0.0223 FIRE: 20 15:58:12 -1064.053533* 0.0379 FIRE: 21 15:58:12 -1064.053619* 0.0376 FIRE: 22 15:58:12 -1064.053786* 0.0370 FIRE: 23 15:58:12 -1064.054025* 0.0360 FIRE: 24 15:58:12 -1064.054321* 0.0347 FIRE: 25 15:58:12 -1064.054658* 0.0331 FIRE: 26 15:58:12 -1064.055018* 0.0311 FIRE: 27 15:58:12 -1064.055382* 0.0289 FIRE: 28 15:58:12 -1064.055770* 0.0261 FIRE: 29 15:58:12 -1064.056158* 0.0226 FIRE: 30 15:58:12 -1064.056517* 0.0182 FIRE: 31 15:58:12 -1064.056819* 0.0131 FIRE: 32 15:58:12 -1064.057048* 0.0080 FIRE: 33 15:58:12 -1064.057208* 0.0089 FIRE: 34 15:58:12 -1064.057321* 0.0091 FIRE: 35 15:58:12 -1064.057413* 0.0106 FIRE: 36 15:58:12 -1064.057501* 0.0137 FIRE: 37 15:58:12 -1064.057584* 0.0146 FIRE: 38 15:58:12 -1064.057635* 0.0125 FIRE: 39 15:58:12 -1064.057650* 0.0121 FIRE: 40 15:58:12 -1064.057677* 0.0114 FIRE: 41 15:58:12 -1064.057715* 0.0103 FIRE: 42 15:58:12 -1064.057757* 0.0090 FIRE: 43 15:58:12 -1064.057801* 0.0074 FIRE: 44 15:58:12 -1064.057840* 0.0057 FIRE: 45 15:58:12 -1064.057873* 0.0040 FIRE: 46 15:58:12 -1064.057900* 0.0030 FIRE: 47 15:58:12 -1064.057918* 0.0024 FIRE: 48 15:58:12 -1064.057924* 0.0024 FIRE: 49 15:58:12 -1064.057925* 0.0024 FIRE: 50 15:58:12 -1064.057926* 0.0023 WARNING: Max number of steps exceeded. Structure may be unstable. Relaxation Completed. Steps: 50 Cell Size Relaxation... Optimization terminated successfully. Current function value: 1.654908 Iterations: 457 Function evaluations: 815 Current VFE: 1.6549081965019923 Energy of Supercell: -1069.9995330904171 Unrelaxed Cell Volume: 2896.6842583156404 Current Relaxed Cell Volume: 2890.599126841865 Current Relaxation Volume: 6.085131473775164 Current Cell: [[ 1.42450055e+01 0.00000000e+00 0.00000000e+00] [-7.87669124e-06 1.42450055e+01 0.00000000e+00] [-3.80835332e-06 -1.24288271e-05 1.42450057e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 15:58:18 -1064.064627* 0.0049 FIRE: 1 15:58:18 -1064.064632* 0.0044 FIRE: 2 15:58:18 -1064.064640* 0.0035 FIRE: 3 15:58:18 -1064.064648* 0.0023 FIRE: 4 15:58:18 -1064.064653* 0.0010 FIRE: 5 15:58:18 -1064.064655* 0.0013 FIRE: 6 15:58:18 -1064.064656* 0.0021 FIRE: 7 15:58:18 -1064.064656* 0.0020 FIRE: 8 15:58:18 -1064.064657* 0.0019 FIRE: 9 15:58:18 -1064.064658* 0.0018 FIRE: 10 15:58:18 -1064.064658* 0.0016 FIRE: 11 15:58:18 -1064.064659* 0.0014 FIRE: 12 15:58:18 -1064.064660* 0.0012 FIRE: 13 15:58:18 -1064.064661* 0.0009 Relaxation Completed. Steps: 13 Cell Size Relaxation... Optimization terminated successfully. Current function value: 1.654874 Iterations: 235 Function evaluations: 456 Current VFE: 1.6548739776828825 Energy of Supercell: -1069.9995330904171 Unrelaxed Cell Volume: 2896.6842583156404 Current Relaxed Cell Volume: 2890.587290572709 Current Relaxation Volume: 6.096967742931156 Current Cell: [[ 1.42449862e+01 0.00000000e+00 0.00000000e+00] [-8.42028991e-06 1.42449860e+01 0.00000000e+00] [-5.80644075e-06 -1.66616689e-06 1.42449861e+01]] ========== Loop: 2 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 15:58:22 -1064.064661* 0.0009 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 1.654874 Iterations: 176 Function evaluations: 363 Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 15:58:24 -1064.064661* 0.0009 FIRE: 1 15:58:24 -1064.064661* 0.0008 FIRE: 2 15:58:24 -1064.064662* 0.0007 FIRE: 3 15:58:24 -1064.064662* 0.0006 FIRE: 4 15:58:24 -1064.064663* 0.0004 FIRE: 5 15:58:24 -1064.064663* 0.0004 FIRE: 6 15:58:24 -1064.064664* 0.0005 FIRE: 7 15:58:24 -1064.064664* 0.0005 FIRE: 8 15:58:24 -1064.064665* 0.0004 FIRE: 9 15:58:24 -1064.064665* 0.0002 FIRE: 10 15:58:24 -1064.064665* 0.0001 FIRE: 11 15:58:24 -1064.064665* 0.0001 FIRE: 12 15:58:24 -1064.064665* 0.0001 FIRE: 13 15:58:24 -1064.064665* 0.0000 FIRE: 14 15:58:24 -1064.064665* 0.0000 FIRE: 15 15:58:24 -1064.064665* 0.0000 FIRE: 16 15:58:24 -1064.064665* 0.0000 FIRE: 17 15:58:24 -1064.064665* 0.0000 FIRE: 18 15:58:24 -1064.064665* 0.0000 FIRE: 19 15:58:24 -1064.064665* 0.0001 FIRE: 20 15:58:24 -1064.064665* 0.0001 Optimization terminated successfully. Current function value: 1.654870 Iterations: 215 Function evaluations: 484 --------------- Calculation Completed. Number Of Atoms in Supercell: 250 Vacancy Formation Energy (relaxed): 1.6548698887042974 Vacancy Formation Energy (unrelaxed): 2.3999067319114147 Unrelaxed Cell Volume: 2896.6842583156404 Relaxed Cell Volume: 2890.587290572709 Relaxation Volume: 6.096967742931156 Relaxed Cell Vector: [14.244986147506866, -1.1819279694028673e-07, 14.244986343407572, -7.197414181430022e-07, -3.451959653600986e-06, 14.244986266696122] Unrelaxed Cell Vector: [14.254994466900826, 0.0, 14.254994466900826, 0.0, 0.0, 14.254994466900826] Relaxed Cell: [[ 1.42449861e+01 0.00000000e+00 0.00000000e+00] [-1.18192797e-07 1.42449863e+01 0.00000000e+00] [-7.19741418e-07 -3.45195965e-06 1.42449863e+01]] Unrelaxed Cell: [[14.25499447 0. 0. ] [ 0. 14.25499447 0. ] [ 0. 0. 14.25499447]] Supercell Size: 6 Unrelaxed Cell: [[17.10599336 0. 0. ] [ 0. 17.10599336 0. ] [ 0. 0. 17.10599336]] Unrelaxed Cell Vector: [17.10599336028099, 0.0, 17.10599336028099, 0.0, 0.0, 17.10599336028099] Unrelaxed Cell Energy: -1848.9591931802488 Energy of Unrelaxed Cell With Vacancy: -1848.9591931802488 ========== Loop: 0 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 15:58:28 -1842.279288* 0.7301 FIRE: 1 15:58:28 -1842.329948* 0.6869 FIRE: 2 15:58:28 -1842.418876* 0.6026 FIRE: 3 15:58:28 -1842.525674* 0.4860 FIRE: 4 15:58:28 -1842.630109* 0.3514 FIRE: 5 15:58:28 -1842.719677* 0.2264 FIRE: 6 15:58:28 -1842.791465* 0.2523 FIRE: 7 15:58:28 -1842.847983* 0.2393 FIRE: 8 15:58:28 -1842.897166* 0.1969 FIRE: 9 15:58:28 -1842.938740* 0.1244 FIRE: 10 15:58:28 -1842.961421* 0.1256 FIRE: 11 15:58:28 -1842.964264* 0.1221 FIRE: 12 15:58:28 -1842.969635* 0.1152 FIRE: 13 15:58:28 -1842.976944* 0.1050 FIRE: 14 15:58:28 -1842.985415* 0.0919 FIRE: 15 15:58:28 -1842.994193* 0.0761 FIRE: 16 15:58:28 -1843.002463* 0.0579 FIRE: 17 15:58:28 -1843.009560* 0.0380 FIRE: 18 15:58:28 -1843.015535* 0.0228 FIRE: 19 15:58:28 -1843.019671* 0.0207 FIRE: 20 15:58:28 -1843.021646* 0.0374 FIRE: 21 15:58:28 -1843.021774* 0.0599 FIRE: 22 15:58:28 -1843.021953* 0.0591 FIRE: 23 15:58:28 -1843.022300* 0.0576 FIRE: 24 15:58:28 -1843.022795* 0.0555 FIRE: 25 15:58:28 -1843.023409* 0.0526 FIRE: 26 15:58:28 -1843.024106* 0.0490 FIRE: 27 15:58:28 -1843.024847* 0.0449 FIRE: 28 15:58:28 -1843.025592* 0.0402 FIRE: 29 15:58:28 -1843.026376* 0.0344 FIRE: 30 15:58:28 -1843.027141* 0.0275 FIRE: 31 15:58:28 -1843.027813* 0.0193 FIRE: 32 15:58:28 -1843.028322* 0.0121 FIRE: 33 15:58:28 -1843.028630* 0.0111 FIRE: 34 15:58:28 -1843.028773* 0.0128 FIRE: 35 15:58:28 -1843.028848* 0.0160 FIRE: 36 15:58:28 -1843.028877* 0.0157 FIRE: 37 15:58:28 -1843.028933* 0.0151 FIRE: 38 15:58:28 -1843.029013* 0.0142 FIRE: 39 15:58:28 -1843.029112* 0.0130 FIRE: 40 15:58:28 -1843.029224* 0.0116 FIRE: 41 15:58:28 -1843.029344* 0.0104 FIRE: 42 15:58:28 -1843.029465* 0.0091 FIRE: 43 15:58:28 -1843.029591* 0.0075 FIRE: 44 15:58:28 -1843.029713* 0.0062 FIRE: 45 15:58:29 -1843.029818* 0.0047 FIRE: 46 15:58:29 -1843.029892* 0.0035 FIRE: 47 15:58:29 -1843.029928* 0.0059 FIRE: 48 15:58:29 -1843.029935* 0.0084 FIRE: 49 15:58:29 -1843.029938* 0.0083 FIRE: 50 15:58:29 -1843.029943* 0.0081 WARNING: Max number of steps exceeded. Structure may be unstable. Relaxation Completed. Steps: 50 Cell Size Relaxation... Optimization terminated successfully. Current function value: 1.645350 Iterations: 310 Function evaluations: 574 Current VFE: 1.6453500752227228 Energy of Supercell: -1848.9591931802488 Unrelaxed Cell Volume: 5005.470398369431 Current Relaxed Cell Volume: 4999.356526818578 Current Relaxation Volume: 6.113871550853219 Current Cell: [[1.70990260e+01 0.00000000e+00 0.00000000e+00] [4.90056792e-05 1.70990257e+01 0.00000000e+00] [1.71428910e-05 6.14370285e-05 1.70990259e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 15:58:36 -1843.033845* 0.0096 FIRE: 1 15:58:36 -1843.033861* 0.0090 FIRE: 2 15:58:36 -1843.033890* 0.0078 FIRE: 3 15:58:36 -1843.033924* 0.0062 FIRE: 4 15:58:36 -1843.033956* 0.0042 FIRE: 5 15:58:36 -1843.033981* 0.0027 FIRE: 6 15:58:36 -1843.033997* 0.0021 FIRE: 7 15:58:36 -1843.034005* 0.0022 FIRE: 8 15:58:36 -1843.034012* 0.0027 FIRE: 9 15:58:36 -1843.034022* 0.0032 FIRE: 10 15:58:36 -1843.034035* 0.0029 FIRE: 11 15:58:36 -1843.034046* 0.0018 FIRE: 12 15:58:36 -1843.034047* 0.0012 FIRE: 13 15:58:36 -1843.034047* 0.0011 FIRE: 14 15:58:36 -1843.034048* 0.0010 FIRE: 15 15:58:36 -1843.034048* 0.0009 Relaxation Completed. Steps: 15 Cell Size Relaxation... Optimization terminated successfully. Current function value: 1.645147 Iterations: 336 Function evaluations: 602 Current VFE: 1.6451467004549158 Energy of Supercell: -1848.9591931802488 Unrelaxed Cell Volume: 5005.470398369431 Current Relaxed Cell Volume: 4999.34893489065 Current Relaxation Volume: 6.121463478781152 Current Cell: [[1.70990173e+01 0.00000000e+00 0.00000000e+00] [2.35298391e-07 1.70990171e+01 0.00000000e+00] [3.79347303e-07 1.47971082e-06 1.70990174e+01]] ========== Loop: 2 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 15:58:43 -1843.034048* 0.0009 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 1.645147 Iterations: 114 Function evaluations: 279 Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 15:58:46 -1843.034048* 0.0009 FIRE: 1 15:58:46 -1843.034049* 0.0009 FIRE: 2 15:58:46 -1843.034049* 0.0008 FIRE: 3 15:58:46 -1843.034050* 0.0007 FIRE: 4 15:58:46 -1843.034050* 0.0006 FIRE: 5 15:58:46 -1843.034051* 0.0005 FIRE: 6 15:58:46 -1843.034052* 0.0005 FIRE: 7 15:58:46 -1843.034052* 0.0004 FIRE: 8 15:58:46 -1843.034053* 0.0003 FIRE: 9 15:58:46 -1843.034053* 0.0002 FIRE: 10 15:58:46 -1843.034054* 0.0001 FIRE: 11 15:58:46 -1843.034054* 0.0002 FIRE: 12 15:58:46 -1843.034054* 0.0002 FIRE: 13 15:58:46 -1843.034054* 0.0002 FIRE: 14 15:58:46 -1843.034054* 0.0002 FIRE: 15 15:58:46 -1843.034054* 0.0002 FIRE: 16 15:58:46 -1843.034054* 0.0001 FIRE: 17 15:58:46 -1843.034054* 0.0001 FIRE: 18 15:58:46 -1843.034054* 0.0000 FIRE: 19 15:58:46 -1843.034054* 0.0000 FIRE: 20 15:58:46 -1843.034054* 0.0000 Optimization terminated successfully. Current function value: 1.645141 Iterations: 180 Function evaluations: 424 --------------- Calculation Completed. Number Of Atoms in Supercell: 432 Vacancy Formation Energy (relaxed): 1.6451412518181314 Vacancy Formation Energy (unrelaxed): 2.399906731912097 Unrelaxed Cell Volume: 5005.470398369431 Relaxed Cell Volume: 4999.34893489065 Relaxation Volume: 6.121463478781152 Relaxed Cell Vector: [17.099016429390673, 2.402209564331744e-07, 17.099016834642153, 3.9193586408684764e-07, 1.4963379882993822e-06, 17.099016764745834] Unrelaxed Cell Vector: [17.10599336028099, 0.0, 17.10599336028099, 0.0, 0.0, 17.10599336028099] Relaxed Cell: [[1.70990164e+01 0.00000000e+00 0.00000000e+00] [2.40220956e-07 1.70990168e+01 0.00000000e+00] [3.91935864e-07 1.49633799e-06 1.70990168e+01]] Unrelaxed Cell: [[17.10599336 0. 0. ] [ 0. 17.10599336 0. ] [ 0. 0. 17.10599336]] Supercell Size: 7 Unrelaxed Cell: [[19.95699225 0. 0. ] [ 0. 19.95699225 0. ] [ 0. 0. 19.95699225]] Unrelaxed Cell Vector: [19.956992253661156, 0.0, 19.956992253661156, 0.0, 0.0, 19.956992253661156] Unrelaxed Cell Energy: -2936.0787188001086 Energy of Unrelaxed Cell With Vacancy: -2936.0787188001086 ========== Loop: 0 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 15:58:51 -2929.398814* 0.7301 FIRE: 1 15:58:51 -2929.449474* 0.6869 FIRE: 2 15:58:51 -2929.538405* 0.6026 FIRE: 3 15:58:51 -2929.645212* 0.4860 FIRE: 4 15:58:51 -2929.749657* 0.3514 FIRE: 5 15:58:51 -2929.839231* 0.2264 FIRE: 6 15:58:51 -2929.911011* 0.2523 FIRE: 7 15:58:51 -2929.967505* 0.2392 FIRE: 8 15:58:51 -2930.016660* 0.1969 FIRE: 9 15:58:51 -2930.058237* 0.1242 FIRE: 10 15:58:51 -2930.081030* 0.1256 FIRE: 11 15:58:51 -2930.083888* 0.1221 FIRE: 12 15:58:51 -2930.089287* 0.1152 FIRE: 13 15:58:51 -2930.096640* 0.1050 FIRE: 14 15:58:52 -2930.105172* 0.0919 FIRE: 15 15:58:52 -2930.114030* 0.0760 FIRE: 16 15:58:52 -2930.122407* 0.0579 FIRE: 17 15:58:52 -2930.129650* 0.0379 FIRE: 18 15:58:52 -2930.135853* 0.0231 FIRE: 19 15:58:52 -2930.140359* 0.0210 FIRE: 20 15:58:52 -2930.142917* 0.0372 FIRE: 21 15:58:52 -2930.143888* 0.0597 FIRE: 22 15:58:52 -2930.144089* 0.0590 FIRE: 23 15:58:52 -2930.144480* 0.0575 FIRE: 24 15:58:52 -2930.145038* 0.0554 FIRE: 25 15:58:52 -2930.145734* 0.0525 FIRE: 26 15:58:52 -2930.146531* 0.0490 FIRE: 27 15:58:52 -2930.147387* 0.0449 FIRE: 28 15:58:52 -2930.148260* 0.0403 FIRE: 29 15:58:52 -2930.149198* 0.0346 FIRE: 30 15:58:52 -2930.150147* 0.0278 FIRE: 31 15:58:52 -2930.151039* 0.0198 FIRE: 32 15:58:52 -2930.151806* 0.0133 FIRE: 33 15:58:52 -2930.152413* 0.0136 FIRE: 34 15:58:52 -2930.152884* 0.0149 FIRE: 35 15:58:52 -2930.153305* 0.0153 FIRE: 36 15:58:52 -2930.153770* 0.0194 FIRE: 37 15:58:52 -2930.154324* 0.0205 FIRE: 38 15:58:52 -2930.154890* 0.0172 FIRE: 39 15:58:52 -2930.155245* 0.0097 FIRE: 40 15:58:52 -2930.155122* 0.0073 FIRE: 41 15:58:52 -2930.155140* 0.0071 FIRE: 42 15:58:52 -2930.155175* 0.0065 FIRE: 43 15:58:52 -2930.155221* 0.0057 FIRE: 44 15:58:52 -2930.155271* 0.0048 FIRE: 45 15:58:52 -2930.155321* 0.0038 FIRE: 46 15:58:52 -2930.155366* 0.0030 FIRE: 47 15:58:52 -2930.155401* 0.0034 FIRE: 48 15:58:52 -2930.155431* 0.0038 FIRE: 49 15:58:52 -2930.155456* 0.0039 FIRE: 50 15:58:52 -2930.155479* 0.0036 WARNING: Max number of steps exceeded. Structure may be unstable. Relaxation Completed. Steps: 50 Cell Size Relaxation... Optimization terminated successfully. Current function value: 1.640778 Iterations: 371 Function evaluations: 672 Current VFE: 1.640777736187374 Energy of Supercell: -2936.0787188001086 Unrelaxed Cell Volume: 7948.501604818125 Current Relaxed Cell Volume: 7942.37320694915 Current Relaxation Volume: 6.1283978689752985 Current Cell: [[1.99518623e+01 0.00000000e+00 0.00000000e+00] [1.88404761e-05 1.99518613e+01 0.00000000e+00] [5.59120354e-05 4.58632297e-06 1.99518621e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 15:59:03 -2930.157943* 0.0027 FIRE: 1 15:59:03 -2930.157949* 0.0025 FIRE: 2 15:59:03 -2930.157960* 0.0022 FIRE: 3 15:59:03 -2930.157973* 0.0018 FIRE: 4 15:59:03 -2930.157985* 0.0015 FIRE: 5 15:59:03 -2930.157995* 0.0012 FIRE: 6 15:59:03 -2930.158002* 0.0014 FIRE: 7 15:59:03 -2930.158008* 0.0015 FIRE: 8 15:59:03 -2930.158013* 0.0018 FIRE: 9 15:59:03 -2930.158020* 0.0017 FIRE: 10 15:59:03 -2930.158028* 0.0012 FIRE: 11 15:59:03 -2930.158033* 0.0004 Relaxation Completed. Steps: 11 Cell Size Relaxation... Optimization terminated successfully. Current function value: 1.640687 Iterations: 225 Function evaluations: 445 Current VFE: 1.6406872184838903 Energy of Supercell: -2936.0787188001086 Unrelaxed Cell Volume: 7948.501604818125 Current Relaxed Cell Volume: 7942.369757376984 Current Relaxation Volume: 6.1318474411409625 Current Cell: [[1.99518591e+01 0.00000000e+00 0.00000000e+00] [7.47323479e-07 1.99518591e+01 0.00000000e+00] [1.27786167e-06 8.26518071e-06 1.99518589e+01]] ========== Loop: 2 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 15:59:11 -2930.158033* 0.0004 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 1.640687 Iterations: 140 Function evaluations: 327 Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 15:59:16 -2930.158033* 0.0004 FIRE: 1 15:59:16 -2930.158034* 0.0004 FIRE: 2 15:59:16 -2930.158034* 0.0003 FIRE: 3 15:59:16 -2930.158034* 0.0003 FIRE: 4 15:59:16 -2930.158035* 0.0002 FIRE: 5 15:59:16 -2930.158035* 0.0002 FIRE: 6 15:59:16 -2930.158035* 0.0002 FIRE: 7 15:59:16 -2930.158035* 0.0002 FIRE: 8 15:59:16 -2930.158035* 0.0002 FIRE: 9 15:59:16 -2930.158035* 0.0002 FIRE: 10 15:59:16 -2930.158035* 0.0002 FIRE: 11 15:59:16 -2930.158035* 0.0002 FIRE: 12 15:59:16 -2930.158035* 0.0001 FIRE: 13 15:59:16 -2930.158035* 0.0001 FIRE: 14 15:59:16 -2930.158035* 0.0001 FIRE: 15 15:59:16 -2930.158035* 0.0001 FIRE: 16 15:59:16 -2930.158035* 0.0001 FIRE: 17 15:59:16 -2930.158035* 0.0001 FIRE: 18 15:59:16 -2930.158035* 0.0001 FIRE: 19 15:59:16 -2930.158035* 0.0001 FIRE: 20 15:59:16 -2930.158035* 0.0001 Optimization terminated successfully. Current function value: 1.640685 Iterations: 200 Function evaluations: 451 --------------- Calculation Completed. Number Of Atoms in Supercell: 686 Vacancy Formation Energy (relaxed): 1.640685299370034 Vacancy Formation Energy (unrelaxed): 2.3999067319114147 Unrelaxed Cell Volume: 7948.501604818125 Relaxed Cell Volume: 7942.369757376984 Relaxation Volume: 6.1318474411409625 Relaxed Cell Vector: [19.951859501839486, 1.1248959233497227e-06, 19.951858624947608, 1.2906500002633481e-06, 2.1267302932944075e-07, 19.95185964922343] Unrelaxed Cell Vector: [19.956992253661156, 0.0, 19.956992253661156, 0.0, 0.0, 19.956992253661156] Relaxed Cell: [[1.99518595e+01 0.00000000e+00 0.00000000e+00] [1.12489592e-06 1.99518586e+01 0.00000000e+00] [1.29065000e-06 2.12673029e-07 1.99518596e+01]] Unrelaxed Cell: [[19.95699225 0. 0. ] [ 0. 19.95699225 0. ] [ 0. 0. 19.95699225]] [Calculation Results Summary] Sizes: [5, 6, 7] Unrelaxed Formation Energy By Size: [2.3999067319114147, 2.399906731912097, 2.3999067319114147] Formation Energy By Size: [1.6548698887042974, 1.6451412518181314, 1.640685299370034] Relaxation Volume By Size: [6.096967742931156, 6.121463478781152, 6.1318474411409625] [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: [2.39990673 2.39990673] Fitting Results: (array([ 2.39990673e+00, -2.02365937e-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.65486989 1.64514125] Fitting Results: (array([1.63177774, 2.88651864]), array([], dtype=float64), 2, array([1.41424176, 0.00238316])) Fit with Size: [5 6] Orders: [0, 3] Matrix A (Ax = y): [[1. 0.008 ] [1. 0.00462963]] Data for Fitting: [6.09696774 6.12146348] Fitting Results: (array([ 6.15511147, -7.26796558]), 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: [2.39990673 2.39990673] Fitting Results: (array([2.39990673e+00, 3.97886493e-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.64514125 1.6406853 ] Fitting Results: (array([1.63310667, 2.59946933]), array([], dtype=float64), 2, array([1.41422363e+00, 1.21209807e-03])) Fit with Size: [6 7] Orders: [0, 3] Matrix A (Ax = y): [[1. 0.00462963] [1. 0.00291545]] Data for Fitting: [6.12146348 6.13184744] Fitting Results: (array([ 6.14950835, -6.05769294]), 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: [2.39990673 2.39990673 2.39990673] Fitting Results: (array([ 2.39990673e+00, -2.80716761e-11]), array([2.99504695e-25]), 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.65486989 1.64514125 1.6406853 ] Fitting Results: (array([1.6323735 , 2.80319856]), array([6.84998696e-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: [6.09696774 6.12146348 6.13184744] Fitting Results: (array([ 6.15259961, -6.91666701]), array([1.21770804e-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: [2.39990673 2.39990673 2.39990673] Fitting Results: (array([ 2.39990673e+00, 3.61709698e-09, -1.55394046e-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.65486989 1.64514125 1.6406853 ] Fitting Results: (array([1.63449991, 1.05994451, 7.4315162 ]), array([], dtype=float64), 3, array([1.73207481e+00, 3.75760407e-03, 3.42874264e-05])) Fit with Size: [5 6 7] Orders: [0, 3, 4] Matrix A (Ax = y): [[1.00000000e+00 8.00000000e-03 1.60000000e-03] [1.00000000e+00 4.62962963e-03 7.71604938e-04] [1.00000000e+00 2.91545190e-03 4.16493128e-04]] Data for Fitting: [6.09696774 6.12146348 6.13184744] Fitting Results: (array([ 6.14363412, 0.43333398, -31.3331562 ]), 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: [2.39990673 2.39990673 2.39990673] Fitting Results: (array([ 2.39990673e+00, 1.86642585e-09, -3.64283485e-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.65486989 1.64514125 1.6406853 ] Fitting Results: (array([ 1.63411761, 1.89717999, 17.42137934]), array([], dtype=float64), 3, array([1.73207409e+00, 3.66323826e-03, 1.50029206e-05])) Fit with Size: [5 6 7] Orders: [0, 3, 5] Matrix A (Ax = y): [[1.00000000e+00 8.00000000e-03 3.20000000e-04] [1.00000000e+00 4.62962963e-03 1.28600823e-04] [1.00000000e+00 2.91545190e-03 5.94990183e-05]] Data for Fitting: [6.09696774 6.12146348 6.13184744] Fitting Results: (array([ 6.14524599, -3.09666304, -73.4529517 ]), 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: [2.39990673 2.39990673 2.39990673] Fitting Results: (array([ 2.39990673e+00, 1.28861582e-09, -1.18048667e-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.65486989 1.64514125 1.6406853 ] Fitting Results: (array([ 1.63386867, 2.17351006, 56.45522476]), array([], dtype=float64), 3, array([1.73207406e+00, 3.65852192e-03, 4.63568288e-06])) Fit with Size: [5 6 7] Orders: [0, 3, 6] Matrix A (Ax = y): [[1.00000000e+00 8.00000000e-03 6.40000000e-05] [1.00000000e+00 4.62962963e-03 2.14334705e-05] [1.00000000e+00 2.91545190e-03 8.49985975e-06]] Data for Fitting: [6.09696774 6.12146348 6.13184744] Fitting Results: (array([ 6.14629556, -4.26174066, -238.02953919]), 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: [[2.3999067319130343, 2.3999067319102543], [2.399906731911788], [2.3999067319073406], [2.39990673190814], [2.399906731908662]] Formation Energy Fits By Size: [[1.63177773961186, 1.6331066715842932], [1.63237349658818], [1.6344999067166843], [1.6341176074113613], [1.633868673862066]] Relaxation Volume Fits By Size: [[6.155111467586092, 6.149508353500957], [6.152599605147249], [6.143634120992344], [6.145245991777926], [6.14629555872737]] [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" 2.3999067319102543 "source-unit" "eV" "source-std-uncert-value" 1.919113856276846e-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.850998893380165 "source-unit" "angstrom" } "host-b" { "source-value" 2.850998893380165 "source-unit" "angstrom" } "host-c" { "source-value" 2.850998893380165 "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.279998132361639 "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.850998893380165 "source-unit" "angstrom" } "reservoir-b" { "source-value" 2.850998893380165 "source-unit" "angstrom" } "reservoir-c" { "source-value" 2.850998893380165 "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.6331066715842932 "source-unit" "eV" "source-std-uncert-value" 0.0013932364541336167 } "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.850998893380165 "source-unit" "angstrom" } "host-b" { "source-value" 2.850998893380165 "source-unit" "angstrom" } "host-c" { "source-value" 2.850998893380165 "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.279998132361639 "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.850998893380165 "source-unit" "angstrom" } "reservoir-b" { "source-value" 2.850998893380165 "source-unit" "angstrom" } "reservoir-c" { "source-value" 2.850998893380165 "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" 6.149508353500957 "source-unit" "angstrom^3" "source-std-uncert-value" 0.005881551077161624 } "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.850998893380165 "source-unit" "angstrom" } "host-b" { "source-value" 2.850998893380165 "source-unit" "angstrom" } "host-c" { "source-value" 2.850998893380165 "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" ] } } ]