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: Si diamond EDIP_LAMMPS_JiangMorganSzlufarska_2012_SiC__MO_667792548433_000 [5.425452411174774] CELL_SIZE_MIN: 3 CELL_SIZE_MAX: 5 Smallest System Size: 216 Largest System Size: 1000 [Calculation] Supercell Size: 3 Unrelaxed Cell: [[16.27635723 0. 0. ] [ 0. 16.27635723 0. ] [ 0. 0. 16.27635723]] Unrelaxed Cell Vector: [16.276357233524323, 0.0, 16.276357233524323, 0.0, 0.0, 16.276357233524323] Unrelaxed Cell Energy: -1003.8002409966629 Energy of Unrelaxed Cell With Vacancy: -1003.8002409966629 ========== Loop: 0 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:11:11 -995.963120* 0.4080 FIRE: 1 16:11:11 -995.980481* 0.4057 FIRE: 2 16:11:11 -996.007787* 0.3980 FIRE: 3 16:11:11 -996.036380* 0.3797 FIRE: 4 16:11:11 -996.062753* 0.3446 FIRE: 5 16:11:11 -996.088057* 0.2899 FIRE: 6 16:11:11 -996.112200* 0.2175 FIRE: 7 16:11:11 -996.131885* 0.1338 FIRE: 8 16:11:11 -996.145401* 0.0998 FIRE: 9 16:11:11 -996.151849* 0.0962 FIRE: 10 16:11:11 -996.152742* 0.0916 FIRE: 11 16:11:11 -996.154342* 0.0828 FIRE: 12 16:11:11 -996.156337* 0.0707 FIRE: 13 16:11:11 -996.158387* 0.0564 FIRE: 14 16:11:11 -996.160229* 0.0426 FIRE: 15 16:11:11 -996.161757* 0.0477 FIRE: 16 16:11:11 -996.163016* 0.0505 FIRE: 17 16:11:11 -996.164210* 0.0473 FIRE: 18 16:11:11 -996.165326* 0.0355 FIRE: 19 16:11:11 -996.166198* 0.0381 FIRE: 20 16:11:11 -996.166562* 0.0412 FIRE: 21 16:11:11 -996.166600* 0.0408 FIRE: 22 16:11:11 -996.166672* 0.0400 FIRE: 23 16:11:11 -996.166773* 0.0388 FIRE: 24 16:11:11 -996.166895* 0.0372 FIRE: 25 16:11:11 -996.167029* 0.0353 FIRE: 26 16:11:11 -996.167167* 0.0331 FIRE: 27 16:11:11 -996.167305* 0.0306 FIRE: 28 16:11:11 -996.167450* 0.0275 FIRE: 29 16:11:11 -996.167600* 0.0239 FIRE: 30 16:11:11 -996.167748* 0.0198 FIRE: 31 16:11:11 -996.167889* 0.0153 FIRE: 32 16:11:11 -996.168011* 0.0107 FIRE: 33 16:11:11 -996.168106* 0.0089 FIRE: 34 16:11:11 -996.168179* 0.0109 FIRE: 35 16:11:11 -996.168243* 0.0107 FIRE: 36 16:11:11 -996.168308* 0.0080 FIRE: 37 16:11:11 -996.168360* 0.0070 FIRE: 38 16:11:11 -996.168378* 0.0063 FIRE: 39 16:11:11 -996.168383* 0.0062 FIRE: 40 16:11:11 -996.168392* 0.0061 FIRE: 41 16:11:11 -996.168403* 0.0059 FIRE: 42 16:11:11 -996.168414* 0.0057 FIRE: 43 16:11:11 -996.168425* 0.0053 FIRE: 44 16:11:11 -996.168434* 0.0049 FIRE: 45 16:11:11 -996.168442* 0.0043 FIRE: 46 16:11:11 -996.168450* 0.0035 FIRE: 47 16:11:11 -996.168455* 0.0024 FIRE: 48 16:11:11 -996.168457* 0.0012 FIRE: 49 16:11:11 -996.168458* 0.0012 FIRE: 50 16:11:11 -996.168458* 0.0012 WARNING: Max number of steps exceeded. Structure may be unstable. Relaxation Completed. Steps: 50 Cell Size Relaxation... Optimization terminated successfully. Current function value: 2.977991 Iterations: 557 Function evaluations: 944 Current VFE: 2.9779908289739296 Energy of Supercell: -1003.8002409966629 Unrelaxed Cell Volume: 4311.929381051275 Current Relaxed Cell Volume: 4321.423593310269 Current Relaxation Volume: -9.494212258993684 Current Cell: [[ 1.62882950e+01 0.00000000e+00 0.00000000e+00] [ 6.36162160e-05 1.62882946e+01 0.00000000e+00] [-3.16367108e-05 -2.82171307e-05 1.62882939e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:11:13 -996.175027* 0.0125 FIRE: 1 16:11:13 -996.175034* 0.0120 FIRE: 2 16:11:13 -996.175047* 0.0109 FIRE: 3 16:11:13 -996.175063* 0.0094 FIRE: 4 16:11:13 -996.175079* 0.0077 FIRE: 5 16:11:13 -996.175094* 0.0058 FIRE: 6 16:11:13 -996.175108* 0.0040 FIRE: 7 16:11:13 -996.175118* 0.0023 FIRE: 8 16:11:13 -996.175124* 0.0016 FIRE: 9 16:11:13 -996.175127* 0.0014 FIRE: 10 16:11:13 -996.175127* 0.0024 FIRE: 11 16:11:13 -996.175128* 0.0024 FIRE: 12 16:11:13 -996.175128* 0.0023 FIRE: 13 16:11:13 -996.175129* 0.0022 FIRE: 14 16:11:13 -996.175130* 0.0020 FIRE: 15 16:11:13 -996.175131* 0.0018 FIRE: 16 16:11:13 -996.175132* 0.0015 FIRE: 17 16:11:13 -996.175133* 0.0012 FIRE: 18 16:11:13 -996.175134* 0.0008 Relaxation Completed. Steps: 18 Cell Size Relaxation... Optimization terminated successfully. Current function value: 2.977881 Iterations: 174 Function evaluations: 377 Current VFE: 2.9778814744388455 Energy of Supercell: -1003.8002409966629 Unrelaxed Cell Volume: 4311.929381051275 Current Relaxed Cell Volume: 4321.605523410978 Current Relaxation Volume: -9.67614235970268 Current Cell: [[ 1.62885229e+01 0.00000000e+00 0.00000000e+00] [ 6.28983242e-05 1.62885232e+01 0.00000000e+00] [-3.19853837e-05 -2.89221841e-05 1.62885231e+01]] ========== Loop: 2 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:11:14 -996.175136* 0.0006 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 2.977881 Iterations: 229 Function evaluations: 446 Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:11:15 -996.175136* 0.0006 FIRE: 1 16:11:15 -996.175136* 0.0005 FIRE: 2 16:11:15 -996.175136* 0.0005 FIRE: 3 16:11:15 -996.175137* 0.0004 FIRE: 4 16:11:15 -996.175137* 0.0003 FIRE: 5 16:11:15 -996.175137* 0.0002 FIRE: 6 16:11:15 -996.175137* 0.0002 FIRE: 7 16:11:15 -996.175137* 0.0001 FIRE: 8 16:11:15 -996.175137* 0.0001 FIRE: 9 16:11:15 -996.175137* 0.0001 FIRE: 10 16:11:15 -996.175137* 0.0001 FIRE: 11 16:11:15 -996.175137* 0.0001 FIRE: 12 16:11:15 -996.175137* 0.0001 FIRE: 13 16:11:15 -996.175137* 0.0001 FIRE: 14 16:11:15 -996.175137* 0.0001 FIRE: 15 16:11:15 -996.175137* 0.0001 FIRE: 16 16:11:15 -996.175137* 0.0001 FIRE: 17 16:11:15 -996.175137* 0.0001 FIRE: 18 16:11:15 -996.175137* 0.0000 FIRE: 19 16:11:15 -996.175137* 0.0000 FIRE: 20 16:11:15 -996.175137* 0.0000 Optimization terminated successfully. Current function value: 2.977881 Iterations: 342 Function evaluations: 651 --------------- Calculation Completed. Number Of Atoms in Supercell: 216 Vacancy Formation Energy (relaxed): 2.977880629504625 Vacancy Formation Energy (unrelaxed): 3.189897295977289 Unrelaxed Cell Volume: 4311.929381051275 Relaxed Cell Volume: 4321.605523410978 Relaxation Volume: -9.67614235970268 Relaxed Cell Vector: [16.28851844126815, 1.4887350618018785e-05, 16.2885185882761, -7.294543688062829e-06, -6.761834329490817e-06, 16.2885183666297] Unrelaxed Cell Vector: [16.276357233524323, 0.0, 16.276357233524323, 0.0, 0.0, 16.276357233524323] Relaxed Cell: [[ 1.62885184e+01 0.00000000e+00 0.00000000e+00] [ 1.48873506e-05 1.62885186e+01 0.00000000e+00] [-7.29454369e-06 -6.76183433e-06 1.62885184e+01]] Unrelaxed Cell: [[16.27635723 0. 0. ] [ 0. 16.27635723 0. ] [ 0. 0. 16.27635723]] Supercell Size: 4 Unrelaxed Cell: [[21.70180964 0. 0. ] [ 0. 21.70180964 0. ] [ 0. 0. 21.70180964]] Unrelaxed Cell Vector: [21.701809644699097, 0.0, 21.701809644699097, 0.0, 0.0, 21.701809644699097] Unrelaxed Cell Energy: -2379.378349029316 Energy of Unrelaxed Cell With Vacancy: -2379.378349029316 ========== Loop: 0 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:11:17 -2371.541228* 0.4080 FIRE: 1 16:11:17 -2371.558589* 0.4057 FIRE: 2 16:11:17 -2371.585895* 0.3980 FIRE: 3 16:11:17 -2371.614488* 0.3797 FIRE: 4 16:11:17 -2371.640861* 0.3446 FIRE: 5 16:11:17 -2371.666165* 0.2899 FIRE: 6 16:11:17 -2371.690308* 0.2175 FIRE: 7 16:11:17 -2371.709994* 0.1338 FIRE: 8 16:11:17 -2371.723515* 0.0998 FIRE: 9 16:11:17 -2371.729977* 0.0962 FIRE: 10 16:11:17 -2371.730871* 0.0916 FIRE: 11 16:11:17 -2371.732475* 0.0828 FIRE: 12 16:11:17 -2371.734478* 0.0707 FIRE: 13 16:11:17 -2371.736541* 0.0564 FIRE: 14 16:11:17 -2371.738406* 0.0426 FIRE: 15 16:11:17 -2371.739971* 0.0477 FIRE: 16 16:11:17 -2371.741288* 0.0505 FIRE: 17 16:11:17 -2371.742577* 0.0474 FIRE: 18 16:11:17 -2371.743846* 0.0359 FIRE: 19 16:11:17 -2371.744948* 0.0381 FIRE: 20 16:11:17 -2371.745632* 0.0414 FIRE: 21 16:11:17 -2371.745866* 0.0408 FIRE: 22 16:11:17 -2371.745942* 0.0404 FIRE: 23 16:11:17 -2371.746084* 0.0397 FIRE: 24 16:11:17 -2371.746274* 0.0385 FIRE: 25 16:11:17 -2371.746494* 0.0370 FIRE: 26 16:11:17 -2371.746721* 0.0352 FIRE: 27 16:11:17 -2371.746941* 0.0329 FIRE: 28 16:11:17 -2371.747146* 0.0302 FIRE: 29 16:11:17 -2371.747353* 0.0268 FIRE: 30 16:11:17 -2371.747561* 0.0226 FIRE: 31 16:11:17 -2371.747763* 0.0177 FIRE: 32 16:11:17 -2371.747945* 0.0122 FIRE: 33 16:11:17 -2371.748083* 0.0086 FIRE: 34 16:11:17 -2371.748172* 0.0132 FIRE: 35 16:11:17 -2371.748251* 0.0151 FIRE: 36 16:11:17 -2371.748364* 0.0134 FIRE: 37 16:11:17 -2371.748508* 0.0084 FIRE: 38 16:11:17 -2371.748633* 0.0078 FIRE: 39 16:11:17 -2371.748698* 0.0065 FIRE: 40 16:11:17 -2371.748707* 0.0059 FIRE: 41 16:11:17 -2371.748722* 0.0058 FIRE: 42 16:11:17 -2371.748739* 0.0057 FIRE: 43 16:11:17 -2371.748754* 0.0054 FIRE: 44 16:11:17 -2371.748765* 0.0051 FIRE: 45 16:11:17 -2371.748773* 0.0045 FIRE: 46 16:11:17 -2371.748781* 0.0039 FIRE: 47 16:11:17 -2371.748788* 0.0030 FIRE: 48 16:11:17 -2371.748793* 0.0019 FIRE: 49 16:11:17 -2371.748792* 0.0012 FIRE: 50 16:11:17 -2371.748792* 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: 2.979529 Iterations: 319 Function evaluations: 600 Current VFE: 2.979529379541418 Energy of Supercell: -2379.378349029316 Unrelaxed Cell Volume: 10220.86964397338 Current Relaxed Cell Volume: 10230.393568172862 Current Relaxation Volume: -9.52392419948228 Current Cell: [[2.17085473e+01 0.00000000e+00 0.00000000e+00] [3.51568555e-05 2.17085487e+01 0.00000000e+00] [7.36281950e-05 2.98916843e-05 2.17085486e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:11:19 -2371.751596* 0.0052 FIRE: 1 16:11:19 -2371.751599* 0.0049 FIRE: 2 16:11:19 -2371.751603* 0.0045 FIRE: 3 16:11:19 -2371.751607* 0.0038 FIRE: 4 16:11:19 -2371.751610* 0.0031 FIRE: 5 16:11:19 -2371.751613* 0.0023 FIRE: 6 16:11:19 -2371.751616* 0.0015 FIRE: 7 16:11:19 -2371.751617* 0.0008 Relaxation Completed. Steps: 7 Cell Size Relaxation... Optimization terminated successfully. Current function value: 2.979509 Iterations: 196 Function evaluations: 410 Current VFE: 2.979508992823412 Energy of Supercell: -2379.378349029316 Unrelaxed Cell Volume: 10220.86964397338 Current Relaxed Cell Volume: 10230.441700873313 Current Relaxation Volume: -9.572056899933159 Current Cell: [[2.17085823e+01 0.00000000e+00 0.00000000e+00] [3.74197980e-05 2.17085828e+01 0.00000000e+00] [2.18832960e-05 3.50966803e-05 2.17085818e+01]] ========== Loop: 2 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:11:20 -2371.751617* 0.0008 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 2.979509 Iterations: 228 Function evaluations: 457 Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:11:21 -2371.751617* 0.0008 FIRE: 1 16:11:21 -2371.751617* 0.0008 FIRE: 2 16:11:21 -2371.751618* 0.0008 FIRE: 3 16:11:21 -2371.751618* 0.0008 FIRE: 4 16:11:21 -2371.751619* 0.0008 FIRE: 5 16:11:21 -2371.751619* 0.0008 FIRE: 6 16:11:21 -2371.751619* 0.0007 FIRE: 7 16:11:21 -2371.751619* 0.0006 FIRE: 8 16:11:21 -2371.751620* 0.0005 FIRE: 9 16:11:21 -2371.751620* 0.0003 FIRE: 10 16:11:21 -2371.751620* 0.0002 FIRE: 11 16:11:21 -2371.751620* 0.0003 FIRE: 12 16:11:21 -2371.751620* 0.0003 FIRE: 13 16:11:21 -2371.751620* 0.0002 FIRE: 14 16:11:21 -2371.751620* 0.0001 FIRE: 15 16:11:21 -2371.751620* 0.0001 FIRE: 16 16:11:21 -2371.751620* 0.0001 FIRE: 17 16:11:21 -2371.751620* 0.0001 FIRE: 18 16:11:21 -2371.751620* 0.0001 FIRE: 19 16:11:21 -2371.751620* 0.0001 FIRE: 20 16:11:21 -2371.751620* 0.0001 Optimization terminated successfully. Current function value: 2.979505 Iterations: 411 Function evaluations: 784 --------------- Calculation Completed. Number Of Atoms in Supercell: 512 Vacancy Formation Energy (relaxed): 2.9795052122576635 Vacancy Formation Energy (unrelaxed): 3.18989729597979 Unrelaxed Cell Volume: 10220.86964397338 Relaxed Cell Volume: 10230.441700873313 Relaxation Volume: -9.572056899933159 Relaxed Cell Vector: [21.708598996914716, 3.5377568985977177e-06, 21.708598080224608, 6.698412655449215e-06, 2.2239662458130988e-06, 21.70859911110891] Unrelaxed Cell Vector: [21.701809644699097, 0.0, 21.701809644699097, 0.0, 0.0, 21.701809644699097] Relaxed Cell: [[2.17085990e+01 0.00000000e+00 0.00000000e+00] [3.53775690e-06 2.17085981e+01 0.00000000e+00] [6.69841266e-06 2.22396625e-06 2.17085991e+01]] Unrelaxed Cell: [[21.70180964 0. 0. ] [ 0. 21.70180964 0. ] [ 0. 0. 21.70180964]] Supercell Size: 5 Unrelaxed Cell: [[27.12726206 0. 0. ] [ 0. 27.12726206 0. ] [ 0. 0. 27.12726206]] Unrelaxed Cell Vector: [27.12726205587387, 0.0, 27.12726205587387, 0.0, 0.0, 27.12726205587387] Unrelaxed Cell Energy: -4647.223337947803 Energy of Unrelaxed Cell With Vacancy: -4647.223337947803 ========== Loop: 0 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:11:24 -4639.386217* 0.4080 FIRE: 1 16:11:24 -4639.403578* 0.4057 FIRE: 2 16:11:24 -4639.430884* 0.3980 FIRE: 3 16:11:24 -4639.459477* 0.3797 FIRE: 4 16:11:24 -4639.485850* 0.3446 FIRE: 5 16:11:25 -4639.511154* 0.2899 FIRE: 6 16:11:25 -4639.535297* 0.2175 FIRE: 7 16:11:25 -4639.554983* 0.1338 FIRE: 8 16:11:25 -4639.568504* 0.0998 FIRE: 9 16:11:25 -4639.574966* 0.0962 FIRE: 10 16:11:25 -4639.575860* 0.0916 FIRE: 11 16:11:25 -4639.577464* 0.0828 FIRE: 12 16:11:25 -4639.579467* 0.0707 FIRE: 13 16:11:25 -4639.581531* 0.0564 FIRE: 14 16:11:25 -4639.583395* 0.0426 FIRE: 15 16:11:25 -4639.584961* 0.0477 FIRE: 16 16:11:25 -4639.586279* 0.0505 FIRE: 17 16:11:25 -4639.587570* 0.0474 FIRE: 18 16:11:25 -4639.588844* 0.0359 FIRE: 19 16:11:25 -4639.589954* 0.0381 FIRE: 20 16:11:25 -4639.590656* 0.0414 FIRE: 21 16:11:25 -4639.590922* 0.0408 FIRE: 22 16:11:25 -4639.591000* 0.0404 FIRE: 23 16:11:25 -4639.591145* 0.0397 FIRE: 24 16:11:25 -4639.591341* 0.0386 FIRE: 25 16:11:25 -4639.591567* 0.0371 FIRE: 26 16:11:25 -4639.591803* 0.0352 FIRE: 27 16:11:25 -4639.592034* 0.0329 FIRE: 28 16:11:25 -4639.592251* 0.0302 FIRE: 29 16:11:25 -4639.592475* 0.0268 FIRE: 30 16:11:25 -4639.592704* 0.0227 FIRE: 31 16:11:25 -4639.592935* 0.0178 FIRE: 32 16:11:25 -4639.593151* 0.0123 FIRE: 33 16:11:25 -4639.593332* 0.0085 FIRE: 34 16:11:25 -4639.593471* 0.0133 FIRE: 35 16:11:25 -4639.593606* 0.0153 FIRE: 36 16:11:25 -4639.593783* 0.0136 FIRE: 37 16:11:25 -4639.593999* 0.0087 FIRE: 38 16:11:25 -4639.594203* 0.0073 FIRE: 39 16:11:25 -4639.594352* 0.0063 FIRE: 40 16:11:25 -4639.594407* 0.0042 FIRE: 41 16:11:25 -4639.594416* 0.0038 FIRE: 42 16:11:25 -4639.594431* 0.0032 FIRE: 43 16:11:25 -4639.594447* 0.0033 FIRE: 44 16:11:25 -4639.594462* 0.0035 FIRE: 45 16:11:25 -4639.594474* 0.0036 FIRE: 46 16:11:25 -4639.594484* 0.0037 FIRE: 47 16:11:25 -4639.594492* 0.0036 FIRE: 48 16:11:25 -4639.594500* 0.0034 FIRE: 49 16:11:25 -4639.594505* 0.0029 FIRE: 50 16:11:25 -4639.594509* 0.0021 WARNING: Max number of steps exceeded. Structure may be unstable. Relaxation Completed. Steps: 50 Cell Size Relaxation... Optimization terminated successfully. Current function value: 2.980161 Iterations: 384 Function evaluations: 685 Current VFE: 2.980160515465286 Energy of Supercell: -4647.223337947803 Unrelaxed Cell Volume: 19962.636023385523 Current Relaxed Cell Volume: 19972.179031767508 Current Relaxation Volume: -9.543008381984691 Current Cell: [[2.71315851e+01 0.00000000e+00 0.00000000e+00] [2.75251661e-05 2.71315842e+01 0.00000000e+00] [2.38326134e-05 3.03383790e-05 2.71315828e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:11:29 -4639.595954* 0.0015 FIRE: 1 16:11:29 -4639.595957* 0.0012 FIRE: 2 16:11:29 -4639.595961* 0.0010 FIRE: 3 16:11:29 -4639.595964* 0.0009 Relaxation Completed. Steps: 3 Cell Size Relaxation... Optimization terminated successfully. Current function value: 2.980150 Iterations: 210 Function evaluations: 426 Current VFE: 2.980150158469769 Energy of Supercell: -4647.223337947803 Unrelaxed Cell Volume: 19962.636023385523 Current Relaxed Cell Volume: 19972.184203788067 Current Relaxation Volume: -9.548180402543949 Current Cell: [[2.71315857e+01 0.00000000e+00 0.00000000e+00] [1.86747440e-05 2.71315874e+01 0.00000000e+00] [1.72917111e-05 2.17462999e-05 2.71315861e+01]] ========== Loop: 2 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:11:32 -4639.595964* 0.0009 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 2.980150 Iterations: 132 Function evaluations: 303 Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:11:33 -4639.595964* 0.0009 FIRE: 1 16:11:33 -4639.595965* 0.0008 FIRE: 2 16:11:33 -4639.595966* 0.0008 FIRE: 3 16:11:33 -4639.595968* 0.0007 FIRE: 4 16:11:33 -4639.595969* 0.0007 FIRE: 5 16:11:33 -4639.595971* 0.0006 FIRE: 6 16:11:33 -4639.595972* 0.0005 FIRE: 7 16:11:33 -4639.595972* 0.0004 FIRE: 8 16:11:33 -4639.595973* 0.0004 FIRE: 9 16:11:33 -4639.595974* 0.0003 FIRE: 10 16:11:33 -4639.595975* 0.0003 FIRE: 11 16:11:33 -4639.595975* 0.0005 FIRE: 12 16:11:33 -4639.595976* 0.0006 FIRE: 13 16:11:33 -4639.595976* 0.0006 FIRE: 14 16:11:33 -4639.595976* 0.0005 FIRE: 15 16:11:33 -4639.595976* 0.0005 FIRE: 16 16:11:33 -4639.595976* 0.0004 FIRE: 17 16:11:33 -4639.595976* 0.0003 FIRE: 18 16:11:33 -4639.595976* 0.0002 FIRE: 19 16:11:33 -4639.595976* 0.0001 FIRE: 20 16:11:33 -4639.595976* 0.0001 Optimization terminated successfully. Current function value: 2.980138 Iterations: 178 Function evaluations: 429 --------------- Calculation Completed. Number Of Atoms in Supercell: 1000 Vacancy Formation Energy (relaxed): 2.980138189950594 Vacancy Formation Energy (unrelaxed): 3.1898972959734238 Unrelaxed Cell Volume: 19962.636023385523 Relaxed Cell Volume: 19972.184203788067 Relaxation Volume: -9.548180402543949 Relaxed Cell Vector: [27.13159352557581, 1.91527645263705e-05, 27.131593322310366, 1.7671878768551214e-05, 2.0987019993105398e-05, 27.13159410119294] Unrelaxed Cell Vector: [27.12726205587387, 0.0, 27.12726205587387, 0.0, 0.0, 27.12726205587387] Relaxed Cell: [[2.71315935e+01 0.00000000e+00 0.00000000e+00] [1.91527645e-05 2.71315933e+01 0.00000000e+00] [1.76718788e-05 2.09870200e-05 2.71315941e+01]] Unrelaxed Cell: [[27.12726206 0. 0. ] [ 0. 27.12726206 0. ] [ 0. 0. 27.12726206]] [Calculation Results Summary] Sizes: [3, 4, 5] Unrelaxed Formation Energy By Size: [3.189897295977289, 3.18989729597979, 3.1898972959734238] Formation Energy By Size: [2.977880629504625, 2.9795052122576635, 2.980138189950594] Relaxation Volume By Size: [-9.67614235970268, -9.572056899933159, -9.548180402543949] [Extrapolation] Fitting with 2 points, including orders [0, 3] Fit with data beginning 0 Fit with Size: [3 4] Orders: [0, 3] Matrix A (Ax = y): [[1. 0.03703704] [1. 0.015625 ]] Data for Fitting: [3.1898973 3.1898973] Fitting Results: (array([ 3.18989730e+00, -1.16796195e-10]), array([], dtype=float64), 2, array([1.41470379, 0.01513535])) Fit with Size: [3 4] Orders: [0, 3] Matrix A (Ax = y): [[1. 0.03703704] [1. 0.015625 ]] Data for Fitting: [2.97788063 2.97950521] Fitting Results: (array([ 2.98069072, -0.07587241]), array([], dtype=float64), 2, array([1.41470379, 0.01513535])) Fit with Size: [3 4] Orders: [0, 3] Matrix A (Ax = y): [[1. 0.03703704] [1. 0.015625 ]] Data for Fitting: [-9.67614236 -9.5720569 ] Fitting Results: (array([-9.49610265, -4.86107228]), array([], dtype=float64), 2, array([1.41470379, 0.01513535])) Fit with data beginning 1 Fit with Size: [4 5] Orders: [0, 3] Matrix A (Ax = y): [[1. 0.015625] [1. 0.008 ]] Data for Fitting: [3.1898973 3.1898973] Fitting Results: (array([3.18989730e+00, 8.34939115e-10]), array([], dtype=float64), 2, array([1.41431223, 0.00539131])) Fit with Size: [4 5] Orders: [0, 3] Matrix A (Ax = y): [[1. 0.015625] [1. 0.008 ]] Data for Fitting: [2.97950521 2.98013819] Fitting Results: (array([ 2.9808023 , -0.08301347]), array([], dtype=float64), 2, array([1.41431223, 0.00539131])) Fit with Size: [4 5] Orders: [0, 3] Matrix A (Ax = y): [[1. 0.015625] [1. 0.008 ]] Data for Fitting: [-9.5720569 -9.5481804] Fitting Results: (array([-9.52312965, -3.13134392]), array([], dtype=float64), 2, array([1.41431223, 0.00539131])) Fitting with 3 points, including orders [0, 3] Fit with data beginning 0 Fit with Size: [3 4 5] Orders: [0, 3] Matrix A (Ax = y): [[1. 0.03703704] [1. 0.015625 ] [1. 0.008 ]] Data for Fitting: [3.1898973 3.1898973 3.1898973] Fitting Results: (array([3.18989730e+00, 7.88637324e-11]), array([1.77572363e-23]), 2, array([1.73240492, 0.02128544])) Fit with Size: [3 4 5] Orders: [0, 3] Matrix A (Ax = y): [[1. 0.03703704] [1. 0.015625 ] [1. 0.008 ]] Data for Fitting: [2.97788063 2.97950521 2.98013819] Fitting Results: (array([ 2.98073855, -0.0773405 ]), array([9.99671561e-10]), 2, array([1.73240492, 0.02128544])) Fit with Size: [3 4 5] Orders: [0, 3] Matrix A (Ax = y): [[1. 0.03703704] [1. 0.015625 ] [1. 0.008 ]] Data for Fitting: [-9.67614236 -9.5720569 -9.5481804 ] Fitting Results: (array([-9.50768965, -4.50546525]), array([5.86526274e-05]), 2, array([1.73240492, 0.02128544])) Fitting with 3 points, including orders [0, 3, 4] Fit with data beginning 0 Fit with Size: [3 4 5] Orders: [0, 3, 4] Matrix A (Ax = y): [[1. 0.03703704 0.01234568] [1. 0.015625 0.00390625] [1. 0.008 0.0016 ]] Data for Fitting: [3.1898973 3.1898973 3.1898973] Fitting Results: (array([ 3.18989730e+00, 3.97465572e-09, -1.03806121e-08]), array([], dtype=float64), 3, array([1.73243559e+00, 2.27351353e-02, 3.80051182e-04])) Fit with Size: [3 4 5] Orders: [0, 3, 4] Matrix A (Ax = y): [[1. 0.03703704 0.01234568] [1. 0.015625 0.00390625] [1. 0.008 0.0016 ]] Data for Fitting: [2.97788063 2.97950521 2.98013819] Fitting Results: (array([ 2.98086614, -0.10657104, 0.07788682]), array([], dtype=float64), 3, array([1.73243559e+00, 2.27351353e-02, 3.80051182e-04])) Fit with Size: [3 4 5] Orders: [0, 3, 4] Matrix A (Ax = y): [[1. 0.03703704 0.01234568] [1. 0.015625 0.00390625] [1. 0.008 0.0016 ]] Data for Fitting: [-9.67614236 -9.5720569 -9.5481804 ] Fitting Results: (array([ -9.53859356, 2.57483714, -18.86596448]), array([], dtype=float64), 3, array([1.73243559e+00, 2.27351353e-02, 3.80051182e-04])) Fitting with 3 points, including orders [0, 3, 5] Fit with data beginning 0 Fit with Size: [3 4 5] Orders: [0, 3, 5] Matrix A (Ax = y): [[1.00000000e+00 3.70370370e-02 4.11522634e-03] [1.00000000e+00 1.56250000e-02 9.76562500e-04] [1.00000000e+00 8.00000000e-03 3.20000000e-04]] Data for Fitting: [3.1898973 3.1898973 3.1898973] Fitting Results: (array([ 3.18989730e+00, 2.19002321e-09, -1.57372135e-08]), array([], dtype=float64), 3, array([1.73240774e+00, 2.14761184e-02, 2.65390853e-04])) Fit with Size: [3 4 5] Orders: [0, 3, 5] Matrix A (Ax = y): [[1.00000000e+00 3.70370370e-02 4.11522634e-03] [1.00000000e+00 1.56250000e-02 9.76562500e-04] [1.00000000e+00 8.00000000e-03 3.20000000e-04]] Data for Fitting: [2.97788063 2.97950521 2.98013819] Fitting Results: (array([ 2.98084585, -0.09318075, 0.11807795]), array([], dtype=float64), 3, array([1.73240774e+00, 2.14761184e-02, 2.65390853e-04])) Fit with Size: [3 4 5] Orders: [0, 3, 5] Matrix A (Ax = y): [[1.00000000e+00 3.70370370e-02 4.11522634e-03] [1.00000000e+00 1.56250000e-02 9.76562500e-04] [1.00000000e+00 8.00000000e-03 3.20000000e-04]] Data for Fitting: [-9.67614236 -9.5720569 -9.5481804 ] Fitting Results: (array([ -9.53367927, -0.66859515, -28.60117564]), array([], dtype=float64), 3, array([1.73240774e+00, 2.14761184e-02, 2.65390853e-04])) Fitting with 3 points, including orders [0, 3, 6] Fit with data beginning 0 Fit with Size: [3 4 5] Orders: [0, 3, 6] Matrix A (Ax = y): [[1.00000000e+00 3.70370370e-02 1.37174211e-03] [1.00000000e+00 1.56250000e-02 2.44140625e-04] [1.00000000e+00 8.00000000e-03 6.40000000e-05]] Data for Fitting: [3.1898973 3.1898973 3.1898973] Fitting Results: (array([ 3.18989730e+00, 1.60930322e-09, -3.27770270e-08]), array([], dtype=float64), 3, array([1.73240519e+00, 2.13086278e-02, 1.28423698e-04])) Fit with Size: [3 4 5] Orders: [0, 3, 6] Matrix A (Ax = y): [[1.00000000e+00 3.70370370e-02 1.37174211e-03] [1.00000000e+00 1.56250000e-02 2.44140625e-04] [1.00000000e+00 8.00000000e-03 6.40000000e-05]] Data for Fitting: [2.97788063 2.97950521 2.98013819] Fitting Results: (array([ 2.98083304, -0.08882355, 0.24592945]), array([], dtype=float64), 3, array([1.73240519e+00, 2.13086278e-02, 1.28423698e-04])) Fit with Size: [3 4 5] Orders: [0, 3, 6] Matrix A (Ax = y): [[1.00000000e+00 3.70370370e-02 1.37174211e-03] [1.00000000e+00 1.56250000e-02 2.44140625e-04] [1.00000000e+00 8.00000000e-03 6.40000000e-05]] Data for Fitting: [-9.67614236 -9.5720569 -9.5481804 ] Fitting Results: (array([ -9.53057587, -1.72400912, -59.5697268 ]), array([], dtype=float64), 3, array([1.73240519e+00, 2.13086278e-02, 1.28423698e-04])) [Fitting Results Summary] Sizes: [3, 4, 5] 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: [[3.189897295981613, 3.189897295966742], [3.1898972959752414], [3.1898972959582363], [3.1898972959609386], [3.1898972959626466]] Formation Energy Fits By Size: [[2.9806907185909597, 2.9808022976939954], [2.9807385547453573], [2.9808661393457334], [2.980845851037548], [2.9808330388754443]] Relaxation Volume Fits By Size: [[-9.496102645506745, -9.523129651184773], [-9.507689654016787], [-9.538593556495401], [-9.533679265148967], [-9.530575867034845]] [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" 3.189897295966742 "source-unit" "eV" "source-std-uncert-value" 1.1968519174846117e-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" [ "diamond" ] } "host-a" { "source-value" 5.425452411174774 "source-unit" "angstrom" } "host-b" { "source-value" 5.425452411174774 "source-unit" "angstrom" } "host-c" { "source-value" 5.425452411174774 "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" "Fd-3m" } "host-wyckoff-multiplicity-and-letter" { "source-value" [ "8a" ] } "host-wyckoff-coordinates" { "source-value" [ [ 0.0 0.0 0.0 ] ] } "host-wyckoff-species" { "source-value" [ "Si" ] } "reservoir-cohesive-potential-energy" { "source-value" 4.647223337947556 "source-unit" "eV" } "reservoir-short-name" { "source-value" [ "diamond" ] } "reservoir-cauchy-stress" { "source-value" [ 0.0 0.0 0.0 0.0 0.0 0.0 ] "source-unit" "GPa" } "reservoir-a" { "source-value" 5.425452411174774 "source-unit" "angstrom" } "reservoir-b" { "source-value" 5.425452411174774 "source-unit" "angstrom" } "reservoir-c" { "source-value" 5.425452411174774 "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" "Fd-3m" } "reservoir-wyckoff-multiplicity-and-letter" { "source-value" [ "8a" ] } "reservoir-wyckoff-coordinates" { "source-value" [ [ 0.0 0.0 0.0 ] ] } "reservoir-wyckoff-species" { "source-value" [ "Si" ] } } { "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.9808022976939954 "source-unit" "eV" "source-std-uncert-value" 6.49538447505258e-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" [ "diamond" ] } "host-a" { "source-value" 5.425452411174774 "source-unit" "angstrom" } "host-b" { "source-value" 5.425452411174774 "source-unit" "angstrom" } "host-c" { "source-value" 5.425452411174774 "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" "Fd-3m" } "host-wyckoff-multiplicity-and-letter" { "source-value" [ "8a" ] } "host-wyckoff-coordinates" { "source-value" [ [ 0.0 0.0 0.0 ] ] } "host-wyckoff-species" { "source-value" [ "Si" ] } "reservoir-cohesive-potential-energy" { "source-value" 4.647223337947556 "source-unit" "eV" } "reservoir-short-name" { "source-value" [ "diamond" ] } "reservoir-cauchy-stress" { "source-value" [ 0.0 0.0 0.0 0.0 0.0 0.0 ] "source-unit" "GPa" } "reservoir-a" { "source-value" 5.425452411174774 "source-unit" "angstrom" } "reservoir-b" { "source-value" 5.425452411174774 "source-unit" "angstrom" } "reservoir-c" { "source-value" 5.425452411174774 "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" "Fd-3m" } "reservoir-wyckoff-multiplicity-and-letter" { "source-value" [ "8a" ] } "reservoir-wyckoff-coordinates" { "source-value" [ [ 0.0 0.0 0.0 ] ] } "reservoir-wyckoff-species" { "source-value" [ "Si" ] } } { "property-id" "tag:staff@noreply.openkim.org,2015-07-28:property/monovacancy-neutral-relaxation-volume-crystal-npt" "instance-id" 3 "relaxation-volume" { "source-value" -9.523129651184773 "source-unit" "angstrom^3" "source-std-uncert-value" 0.022291051840446295 } "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" [ "diamond" ] } "host-a" { "source-value" 5.425452411174774 "source-unit" "angstrom" } "host-b" { "source-value" 5.425452411174774 "source-unit" "angstrom" } "host-c" { "source-value" 5.425452411174774 "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" "Fd-3m" } "host-wyckoff-multiplicity-and-letter" { "source-value" [ "8a" ] } "host-wyckoff-coordinates" { "source-value" [ [ 0.0 0.0 0.0 ] ] } "host-wyckoff-species" { "source-value" [ "Si" ] } } ]