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: Mo bcc Three_Body_Stillinger_Weber_MoS__MO_201919462778_000 [2.65344320983] CELL_SIZE_MIN: 5 CELL_SIZE_MAX: 7 Smallest System Size: 250 Largest System Size: 686 [Calculation] Supercell Size: 5 Unrelaxed Cell: [[13.26721605 0. 0. ] [ 0. 13.26721605 0. ] [ 0. 0. 13.26721605]] Unrelaxed Cell Vector: [13.267216049149999, 0.0, 13.267216049149999, 0.0, 0.0, 13.267216049149999] Unrelaxed Cell Energy: -2096.50045436 Energy of Unrelaxed Cell With Vacancy: -2096.50045436 ========== Loop: 0 Position Relaxation... Step Time Energy fmax FIRE: 0 21:59:14 -2079.728451 4.8561 FIRE: 1 21:59:14 -2081.391495 2.8983 FIRE: 2 21:59:14 -2082.189383 2.0197 FIRE: 3 21:59:14 -2082.304604 1.8085 FIRE: 4 21:59:14 -2082.488622 1.4185 FIRE: 5 21:59:14 -2082.678545 0.9140 FIRE: 6 21:59:14 -2082.836006 1.0248 FIRE: 7 21:59:14 -2082.966403 1.0475 FIRE: 8 21:59:14 -2083.091699 0.9032 FIRE: 9 21:59:14 -2083.210089 0.5956 FIRE: 10 21:59:14 -2083.295539 0.2423 FIRE: 11 21:59:14 -2083.294063 0.4601 FIRE: 12 21:59:14 -2083.297880 0.4424 FIRE: 13 21:59:14 -2083.304796 0.4078 FIRE: 14 21:59:14 -2083.313559 0.3584 FIRE: 15 21:59:14 -2083.322696 0.2966 FIRE: 16 21:59:14 -2083.330903 0.2263 FIRE: 17 21:59:14 -2083.337391 0.1574 FIRE: 18 21:59:14 -2083.342071 0.1696 FIRE: 19 21:59:14 -2083.345673 0.1631 FIRE: 20 21:59:14 -2083.348428 0.1318 FIRE: 21 21:59:14 -2083.350200 0.0941 FIRE: 22 21:59:14 -2083.350114 0.1154 FIRE: 23 21:59:14 -2083.350228 0.1132 FIRE: 24 21:59:14 -2083.350445 0.1088 FIRE: 25 21:59:14 -2083.350749 0.1023 FIRE: 26 21:59:14 -2083.351113 0.0939 FIRE: 27 21:59:14 -2083.351509 0.0838 FIRE: 28 21:59:14 -2083.351907 0.0724 FIRE: 29 21:59:14 -2083.352281 0.0599 FIRE: 30 21:59:14 -2083.352643 0.0454 FIRE: 31 21:59:14 -2083.352958 0.0293 FIRE: 32 21:59:14 -2083.353198 0.0268 FIRE: 33 21:59:14 -2083.353341 0.0221 FIRE: 34 21:59:14 -2083.353382 0.0186 FIRE: 35 21:59:14 -2083.353388 0.0183 FIRE: 36 21:59:14 -2083.353398 0.0176 FIRE: 37 21:59:14 -2083.353413 0.0166 FIRE: 38 21:59:14 -2083.353430 0.0153 FIRE: 39 21:59:14 -2083.353450 0.0137 FIRE: 40 21:59:14 -2083.353469 0.0119 FIRE: 41 21:59:14 -2083.353487 0.0099 FIRE: 42 21:59:14 -2083.353504 0.0075 FIRE: 43 21:59:14 -2083.353518 0.0048 FIRE: 44 21:59:14 -2083.353526 0.0044 FIRE: 45 21:59:14 -2083.353529 0.0049 FIRE: 46 21:59:14 -2083.353529 0.0048 FIRE: 47 21:59:14 -2083.353529 0.0048 FIRE: 48 21:59:14 -2083.353530 0.0046 FIRE: 49 21:59:14 -2083.353531 0.0045 WARNING: Max number of steps exceeded. Structure may be unstable. Relaxation Completed. Steps: 50 Cell Size Relaxation... Optimization terminated successfully. Current function value: 4.759717 Iterations: 328 Function evaluations: 604 Current VFE: 4.75971721315 Energy of Supercell: -2096.50045436 Unrelaxed Cell Volume: 2335.28238959 Current Relaxed Cell Volume: 2334.15692535 Current Relaxation Volume: 1.1254642341 Current Cell: [[1.32650844e+01 0.00000000e+00 0.00000000e+00] [2.20993755e-05 1.32650847e+01 0.00000000e+00] [2.90433988e-05 2.73948651e-05 1.32650840e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax FIRE: 0 21:59:39 -2083.354735 0.0044 FIRE: 1 21:59:39 -2083.354742 0.0028 FIRE: 2 21:59:39 -2083.354746 0.0032 FIRE: 3 21:59:39 -2083.354747 0.0029 FIRE: 4 21:59:39 -2083.354749 0.0025 FIRE: 5 21:59:39 -2083.354750 0.0019 FIRE: 6 21:59:40 -2083.354751 0.0015 FIRE: 7 21:59:40 -2083.354752 0.0016 FIRE: 8 21:59:40 -2083.354753 0.0015 FIRE: 9 21:59:40 -2083.354754 0.0011 FIRE: 10 21:59:40 -2083.354755 0.0005 Relaxation Completed. Steps: 10 Cell Size Relaxation... Optimization terminated successfully. Current function value: 4.759698 Iterations: 258 Function evaluations: 504 Current VFE: 4.75969790273 Energy of Supercell: -2096.50045436 Unrelaxed Cell Volume: 2335.28238959 Current Relaxed Cell Volume: 2334.15795753 Current Relaxation Volume: 1.12443205776 Current Cell: [[1.32650863e+01 0.00000000e+00 0.00000000e+00] [1.12840499e-07 1.32650863e+01 0.00000000e+00] [5.29538809e-08 8.64381499e-08 1.32650864e+01]] ========== Loop: 2 Position Relaxation... Step Time Energy fmax FIRE: 0 22:00:03 -2083.354755 0.0005 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 4.759698 Iterations: 117 Function evaluations: 283 Step Time Energy fmax FIRE: 0 22:00:16 -2083.354755 0.0005 FIRE: 1 22:00:16 -2083.354755 0.0004 FIRE: 2 22:00:16 -2083.354755 0.0004 FIRE: 3 22:00:16 -2083.354755 0.0004 FIRE: 4 22:00:16 -2083.354755 0.0002 FIRE: 5 22:00:16 -2083.354755 0.0002 FIRE: 6 22:00:16 -2083.354755 0.0002 FIRE: 7 22:00:16 -2083.354755 0.0001 FIRE: 8 22:00:16 -2083.354755 0.0001 FIRE: 9 22:00:16 -2083.354755 0.0001 FIRE: 10 22:00:16 -2083.354755 0.0001 FIRE: 11 22:00:16 -2083.354755 0.0001 FIRE: 12 22:00:16 -2083.354755 0.0001 FIRE: 13 22:00:16 -2083.354755 0.0001 FIRE: 14 22:00:16 -2083.354755 0.0001 FIRE: 15 22:00:16 -2083.354755 0.0000 FIRE: 16 22:00:16 -2083.354755 0.0000 FIRE: 17 22:00:16 -2083.354755 0.0000 FIRE: 18 22:00:16 -2083.354755 0.0000 FIRE: 19 22:00:16 -2083.354755 0.0000 Optimization terminated successfully. Current function value: 4.759698 Iterations: 168 Function evaluations: 402 --------------- Calculation Completed. Number Of Atoms in Supercell: 250 Vacancy Formation Energy (relaxed): 4.7596975582 Vacancy Formation Energy (unrelaxed): 8.38600181744 Unrelaxed Cell Volume: 2335.28238959 Relaxed Cell Volume: 2334.15795753 Relaxation Volume: 1.12443205776 Relaxed Cell Vector: [13.26508648557122, 1.1337230436791986e-07, 13.2650867166783, 5.059176614728805e-08, 8.953788725399615e-08, 13.26508687834066] Unrelaxed Cell Vector: [13.267216049149999, 0.0, 13.267216049149999, 0.0, 0.0, 13.267216049149999] Relaxed Cell: [[1.32650865e+01 0.00000000e+00 0.00000000e+00] [1.13372304e-07 1.32650867e+01 0.00000000e+00] [5.05917661e-08 8.95378873e-08 1.32650869e+01]] Unrelaxed Cell: [[13.26721605 0. 0. ] [ 0. 13.26721605 0. ] [ 0. 0. 13.26721605]] Supercell Size: 6 Unrelaxed Cell: [[15.92065926 0. 0. ] [ 0. 15.92065926 0. ] [ 0. 0. 15.92065926]] Unrelaxed Cell Vector: [15.920659258979999, 0.0, 15.920659258979999, 0.0, 0.0, 15.920659258979999] Unrelaxed Cell Energy: -3622.75278513 Energy of Unrelaxed Cell With Vacancy: -3622.75278513 ========== Loop: 0 Position Relaxation... Step Time Energy fmax FIRE: 0 22:00:34 -3605.980781 4.8561 FIRE: 1 22:00:34 -3607.643798 2.8983 FIRE: 2 22:00:34 -3608.442190 2.0193 FIRE: 3 22:00:34 -3608.557510 1.8082 FIRE: 4 22:00:34 -3608.741562 1.4182 FIRE: 5 22:00:34 -3608.931300 0.9142 FIRE: 6 22:00:34 -3609.088410 1.0245 FIRE: 7 22:00:34 -3609.218528 1.0470 FIRE: 8 22:00:34 -3609.343771 0.9026 FIRE: 9 22:00:34 -3609.462204 0.5951 FIRE: 10 22:00:34 -3609.547405 0.2395 FIRE: 11 22:00:34 -3609.545053 0.4609 FIRE: 12 22:00:34 -3609.548933 0.4431 FIRE: 13 22:00:34 -3609.555986 0.4085 FIRE: 14 22:00:35 -3609.564978 0.3589 FIRE: 15 22:00:35 -3609.574465 0.2971 FIRE: 16 22:00:35 -3609.583160 0.2267 FIRE: 17 22:00:35 -3609.590286 0.1600 FIRE: 18 22:00:35 -3609.595734 0.1729 FIRE: 19 22:00:35 -3609.600281 0.1672 FIRE: 20 22:00:35 -3609.604124 0.1361 FIRE: 21 22:00:35 -3609.607077 0.0930 FIRE: 22 22:00:35 -3609.608181 0.1139 FIRE: 23 22:00:35 -3609.608303 0.1116 FIRE: 24 22:00:35 -3609.608536 0.1073 FIRE: 25 22:00:35 -3609.608862 0.1008 FIRE: 26 22:00:35 -3609.609254 0.0925 FIRE: 27 22:00:35 -3609.609683 0.0826 FIRE: 28 22:00:35 -3609.610117 0.0713 FIRE: 29 22:00:35 -3609.610528 0.0590 FIRE: 30 22:00:35 -3609.610931 0.0448 FIRE: 31 22:00:36 -3609.611290 0.0299 FIRE: 32 22:00:36 -3609.611574 0.0303 FIRE: 33 22:00:36 -3609.611767 0.0286 FIRE: 34 22:00:36 -3609.611869 0.0242 FIRE: 35 22:00:36 -3609.611890 0.0239 FIRE: 36 22:00:36 -3609.611900 0.0235 FIRE: 37 22:00:36 -3609.611919 0.0226 FIRE: 38 22:00:36 -3609.611947 0.0214 FIRE: 39 22:00:36 -3609.611981 0.0198 FIRE: 40 22:00:36 -3609.612019 0.0179 FIRE: 41 22:00:36 -3609.612060 0.0157 FIRE: 42 22:00:36 -3609.612100 0.0133 FIRE: 43 22:00:36 -3609.612143 0.0105 FIRE: 44 22:00:36 -3609.612184 0.0102 FIRE: 45 22:00:36 -3609.612221 0.0097 FIRE: 46 22:00:36 -3609.612251 0.0089 FIRE: 47 22:00:36 -3609.612273 0.0076 FIRE: 48 22:00:36 -3609.612288 0.0096 FIRE: 49 22:00:36 -3609.612298 0.0108 WARNING: Max number of steps exceeded. Structure may be unstable. Relaxation Completed. Steps: 50 Cell Size Relaxation... Optimization terminated successfully. Current function value: 4.753787 Iterations: 388 Function evaluations: 692 Current VFE: 4.75378658856 Energy of Supercell: -3622.75278513 Unrelaxed Cell Volume: 4035.3679692 Current Relaxed Cell Volume: 4034.24993923 Current Relaxation Volume: 1.11802997803 Current Cell: [[ 1.59191880e+01 0.00000000e+00 0.00000000e+00] [-1.57866682e-07 1.59191889e+01 0.00000000e+00] [-3.41103508e-07 -1.70572312e-07 1.59191896e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax FIRE: 0 22:01:23 -3609.612997 0.0100 FIRE: 1 22:01:23 -3609.613025 0.0066 FIRE: 2 22:01:23 -3609.613047 0.0036 FIRE: 3 22:01:23 -3609.613045 0.0046 FIRE: 4 22:01:23 -3609.613047 0.0041 FIRE: 5 22:01:23 -3609.613051 0.0032 FIRE: 6 22:01:23 -3609.613054 0.0019 FIRE: 7 22:01:23 -3609.613056 0.0016 FIRE: 8 22:01:23 -3609.613057 0.0015 FIRE: 9 22:01:23 -3609.613057 0.0015 FIRE: 10 22:01:23 -3609.613057 0.0014 FIRE: 11 22:01:23 -3609.613057 0.0014 FIRE: 12 22:01:23 -3609.613058 0.0013 FIRE: 13 22:01:23 -3609.613058 0.0012 FIRE: 14 22:01:24 -3609.613058 0.0010 FIRE: 15 22:01:24 -3609.613058 0.0009 Relaxation Completed. Steps: 15 Cell Size Relaxation... Optimization terminated successfully. Current function value: 4.753725 Iterations: 117 Function evaluations: 289 Current VFE: 4.75372508497 Energy of Supercell: -3622.75278513 Unrelaxed Cell Volume: 4035.3679692 Current Relaxed Cell Volume: 4034.25109222 Current Relaxation Volume: 1.11687697987 Current Cell: [[ 1.59191906e+01 0.00000000e+00 0.00000000e+00] [-1.63952261e-07 1.59191901e+01 0.00000000e+00] [-3.35079570e-07 -1.76113484e-07 1.59191903e+01]] ========== Loop: 2 Position Relaxation... Step Time Energy fmax FIRE: 0 22:01:44 -3609.613058 0.0009 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 4.753725 Iterations: 111 Function evaluations: 270 Step Time Energy fmax FIRE: 0 22:02:02 -3609.613058 0.0009 FIRE: 1 22:02:02 -3609.613058 0.0008 FIRE: 2 22:02:02 -3609.613059 0.0007 FIRE: 3 22:02:02 -3609.613059 0.0005 FIRE: 4 22:02:02 -3609.613060 0.0003 FIRE: 5 22:02:02 -3609.613060 0.0001 FIRE: 6 22:02:02 -3609.613060 0.0001 FIRE: 7 22:02:02 -3609.613060 0.0001 FIRE: 8 22:02:02 -3609.613060 0.0000 FIRE: 9 22:02:02 -3609.613060 0.0000 FIRE: 10 22:02:02 -3609.613060 0.0000 FIRE: 11 22:02:02 -3609.613060 0.0000 FIRE: 12 22:02:02 -3609.613060 0.0000 FIRE: 13 22:02:02 -3609.613060 0.0000 FIRE: 14 22:02:02 -3609.613060 0.0000 FIRE: 15 22:02:02 -3609.613060 0.0000 FIRE: 16 22:02:03 -3609.613060 0.0000 FIRE: 17 22:02:03 -3609.613060 0.0000 FIRE: 18 22:02:03 -3609.613060 0.0000 FIRE: 19 22:02:03 -3609.613060 0.0000 Optimization terminated successfully. Current function value: 4.753724 Iterations: 159 Function evaluations: 390 --------------- Calculation Completed. Number Of Atoms in Supercell: 432 Vacancy Formation Energy (relaxed): 4.75372355526 Vacancy Formation Energy (unrelaxed): 8.38600181744 Unrelaxed Cell Volume: 4035.3679692 Relaxed Cell Volume: 4034.25109222 Relaxation Volume: 1.11687697987 Relaxed Cell Vector: [15.919190704263594, -1.6446461630591166e-07, 15.91919002213476, -3.403151789141921e-07, -1.7666383951807902e-07, 15.919190146186917] Unrelaxed Cell Vector: [15.920659258979999, 0.0, 15.920659258979999, 0.0, 0.0, 15.920659258979999] Relaxed Cell: [[ 1.59191907e+01 0.00000000e+00 0.00000000e+00] [-1.64464616e-07 1.59191900e+01 0.00000000e+00] [-3.40315179e-07 -1.76663840e-07 1.59191901e+01]] Unrelaxed Cell: [[15.92065926 0. 0. ] [ 0. 15.92065926 0. ] [ 0. 0. 15.92065926]] Supercell Size: 7 Unrelaxed Cell: [[18.57410247 0. 0. ] [ 0. 18.57410247 0. ] [ 0. 0. 18.57410247]] Unrelaxed Cell Vector: [18.574102468809997, 0.0, 18.574102468809997, 0.0, 0.0, 18.574102468809997] Unrelaxed Cell Energy: -5752.79724676 Energy of Unrelaxed Cell With Vacancy: -5752.79724676 ========== Loop: 0 Position Relaxation... Step Time Energy fmax FIRE: 0 22:02:29 -5736.025243 4.8561 FIRE: 1 22:02:29 -5737.688260 2.8983 FIRE: 2 22:02:29 -5738.486649 2.0193 FIRE: 3 22:02:29 -5738.601981 1.8082 FIRE: 4 22:02:29 -5738.786062 1.4182 FIRE: 5 22:02:29 -5738.975847 0.9142 FIRE: 6 22:02:29 -5739.133008 1.0245 FIRE: 7 22:02:30 -5739.263156 1.0470 FIRE: 8 22:02:30 -5739.388376 0.9026 FIRE: 9 22:02:30 -5739.506711 0.5950 FIRE: 10 22:02:30 -5739.591747 0.2388 FIRE: 11 22:02:30 -5739.589343 0.4610 FIRE: 12 22:02:30 -5739.593240 0.4433 FIRE: 13 22:02:30 -5739.600324 0.4086 FIRE: 14 22:02:30 -5739.609352 0.3591 FIRE: 15 22:02:30 -5739.618868 0.2973 FIRE: 16 22:02:30 -5739.627581 0.2270 FIRE: 17 22:02:30 -5739.634710 0.1596 FIRE: 18 22:02:30 -5739.640159 0.1727 FIRE: 19 22:02:30 -5739.644739 0.1673 FIRE: 20 22:02:31 -5739.648699 0.1368 FIRE: 21 22:02:31 -5739.651930 0.0931 FIRE: 22 22:02:31 -5739.653528 0.1146 FIRE: 23 22:02:31 -5739.653663 0.1124 FIRE: 24 22:02:31 -5739.653923 0.1079 FIRE: 25 22:02:31 -5739.654287 0.1014 FIRE: 26 22:02:31 -5739.654725 0.0930 FIRE: 27 22:02:31 -5739.655207 0.0830 FIRE: 28 22:02:31 -5739.655699 0.0716 FIRE: 29 22:02:31 -5739.656172 0.0592 FIRE: 30 22:02:31 -5739.656646 0.0450 FIRE: 31 22:02:31 -5739.657090 0.0303 FIRE: 32 22:02:32 -5739.657477 0.0309 FIRE: 33 22:02:32 -5739.657796 0.0297 FIRE: 34 22:02:32 -5739.658045 0.0258 FIRE: 35 22:02:32 -5739.658221 0.0245 FIRE: 36 22:02:32 -5739.658321 0.0267 FIRE: 37 22:02:32 -5739.658334 0.0262 FIRE: 38 22:02:32 -5739.658361 0.0252 FIRE: 39 22:02:32 -5739.658400 0.0238 FIRE: 40 22:02:32 -5739.658448 0.0219 FIRE: 41 22:02:32 -5739.658503 0.0196 FIRE: 42 22:02:32 -5739.658562 0.0170 FIRE: 43 22:02:32 -5739.658622 0.0153 FIRE: 44 22:02:32 -5739.658686 0.0143 FIRE: 45 22:02:32 -5739.658752 0.0131 FIRE: 46 22:02:33 -5739.658813 0.0114 FIRE: 47 22:02:33 -5739.658864 0.0093 FIRE: 48 22:02:33 -5739.658901 0.0068 FIRE: 49 22:02:33 -5739.658925 0.0101 WARNING: Max number of steps exceeded. Structure may be unstable. Relaxation Completed. Steps: 50 Cell Size Relaxation... Optimization terminated successfully. Current function value: 4.751872 Iterations: 362 Function evaluations: 653 Current VFE: 4.75187224506 Energy of Supercell: -5752.79724676 Unrelaxed Cell Volume: 6408.01487702 Current Relaxed Cell Volume: 6406.90118701 Current Relaxation Volume: 1.11369000894 Current Cell: [[ 1.85730267e+01 0.00000000e+00 0.00000000e+00] [ 1.97725749e-06 1.85730262e+01 0.00000000e+00] [-8.49032553e-07 8.64222673e-07 1.85730262e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax FIRE: 0 22:03:44 -5739.659373 0.0119 FIRE: 1 22:03:44 -5739.659409 0.0085 FIRE: 2 22:03:44 -5739.659445 0.0040 FIRE: 3 22:03:44 -5739.659456 0.0043 FIRE: 4 22:03:44 -5739.659458 0.0039 FIRE: 5 22:03:45 -5739.659462 0.0031 FIRE: 6 22:03:45 -5739.659467 0.0028 FIRE: 7 22:03:45 -5739.659471 0.0025 FIRE: 8 22:03:45 -5739.659474 0.0021 FIRE: 9 22:03:45 -5739.659476 0.0017 FIRE: 10 22:03:45 -5739.659478 0.0021 FIRE: 11 22:03:45 -5739.659480 0.0022 FIRE: 12 22:03:45 -5739.659482 0.0019 FIRE: 13 22:03:45 -5739.659484 0.0013 FIRE: 14 22:03:45 -5739.659484 0.0009 Relaxation Completed. Steps: 14 Cell Size Relaxation... Optimization terminated successfully. Current function value: 4.751761 Iterations: 143 Function evaluations: 323 Current VFE: 4.75176067539 Energy of Supercell: -5752.79724676 Unrelaxed Cell Volume: 6408.01487702 Current Relaxed Cell Volume: 6406.90208248 Current Relaxation Volume: 1.11279454508 Current Cell: [[ 1.85730272e+01 0.00000000e+00 0.00000000e+00] [ 2.01103138e-06 1.85730275e+01 0.00000000e+00] [-8.47754607e-07 8.59698346e-07 1.85730271e+01]] ========== Loop: 2 Position Relaxation... Step Time Energy fmax FIRE: 0 22:04:19 -5739.659484 0.0009 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 4.751761 Iterations: 129 Function evaluations: 296 Step Time Energy fmax FIRE: 0 22:04:52 -5739.659484 0.0009 FIRE: 1 22:04:52 -5739.659485 0.0006 FIRE: 2 22:04:52 -5739.659485 0.0003 FIRE: 3 22:04:52 -5739.659485 0.0005 FIRE: 4 22:04:52 -5739.659485 0.0004 FIRE: 5 22:04:52 -5739.659485 0.0003 FIRE: 6 22:04:52 -5739.659485 0.0002 FIRE: 7 22:04:52 -5739.659485 0.0002 FIRE: 8 22:04:52 -5739.659485 0.0003 FIRE: 9 22:04:52 -5739.659485 0.0003 FIRE: 10 22:04:52 -5739.659485 0.0002 FIRE: 11 22:04:52 -5739.659485 0.0002 FIRE: 12 22:04:53 -5739.659485 0.0002 FIRE: 13 22:04:53 -5739.659485 0.0002 FIRE: 14 22:04:53 -5739.659485 0.0001 FIRE: 15 22:04:53 -5739.659485 0.0001 FIRE: 16 22:04:53 -5739.659485 0.0001 FIRE: 17 22:04:53 -5739.659485 0.0001 FIRE: 18 22:04:53 -5739.659485 0.0001 FIRE: 19 22:04:53 -5739.659485 0.0001 Optimization terminated successfully. Current function value: 4.751760 Iterations: 170 Function evaluations: 413 --------------- Calculation Completed. Number Of Atoms in Supercell: 686 Vacancy Formation Energy (relaxed): 4.7517600015 Vacancy Formation Energy (unrelaxed): 8.38600181745 Unrelaxed Cell Volume: 6408.01487702 Relaxed Cell Volume: 6406.90208248 Relaxation Volume: 1.11279454508 Relaxed Cell Vector: [18.57302698093978, 1.9995534620320833e-06, 18.573026742154877, -8.790491480610926e-07, 8.8596794473868e-07, 18.573027047974037] Unrelaxed Cell Vector: [18.574102468809997, 0.0, 18.574102468809997, 0.0, 0.0, 18.574102468809997] Relaxed Cell: [[ 1.85730270e+01 0.00000000e+00 0.00000000e+00] [ 1.99955346e-06 1.85730267e+01 0.00000000e+00] [-8.79049148e-07 8.85967945e-07 1.85730270e+01]] Unrelaxed Cell: [[18.57410247 0. 0. ] [ 0. 18.57410247 0. ] [ 0. 0. 18.57410247]] [Calculation Results Summary] Sizes: [5, 6, 7] Unrelaxed Formation Energy By Size: [8.386001817437773, 8.386001817435954, 8.386001817448232] Formation Energy By Size: [4.759697558201879, 4.753723555255419, 4.751760001495313] Relaxation Volume By Size: [1.1244320577598046, 1.1168769798732683, 1.1127945450789412] [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: [8.38600182 8.38600182] Fitting Results: (array([8.38600182e+00, 5.39642505e-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: [4.75969756 4.75372356] Fitting Results: (array([4.74551751, 1.77250637]), 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.12443206 1.11687698] Fitting Results: (array([1.10649913, 2.24161652]), 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: [8.38600182 8.38600182] Fitting Results: (array([ 8.38600182e+00, -7.16342234e-09]), array([], dtype=float64), 2, array([1.41422363e+00, 1.21209807e-03])) Fit with Size: [6 7] Orders: [0, 3] Matrix A (Ax = y): [[1. 0.00462963] [1. 0.00291545]] Data for Fitting: [4.75372356 4.75176 ] Fitting Results: (array([4.74842041, 1.14547851]), 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.11687698 1.11279455] Fitting Results: (array([1.10585119, 2.38157031]), 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: [8.38600182 8.38600182 8.38600182] Fitting Results: (array([ 8.38600182e+00, -1.69596058e-09]), array([4.93222974e-23]), 2, array([1.73207406, 0.00365829])) Fit with Size: [5 6 7] Orders: [0, 3] Matrix A (Ax = y): [[1. 0.008 ] [1. 0.00462963] [1. 0.00291545]] Data for Fitting: [4.75969756 4.75372356 4.75176 ] Fitting Results: (array([4.74681887, 1.59050276]), array([3.26851448e-07]), 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.12443206 1.11687698 1.11279455] Fitting Results: (array([1.10620866, 2.28224006]), array([1.62834371e-08]), 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: [8.38600182 8.38600182 8.38600182] Fitting Results: (array([ 8.38600182e+00, -4.84735092e-08, 1.99413339e-07]), array([], dtype=float64), 3, array([1.73207481e+00, 3.75760407e-03, 3.42874264e-05])) Fit with Size: [5 6 7] Orders: [0, 3, 4] Matrix A (Ax = y): [[1.00000000e+00 8.00000000e-03 1.60000000e-03] [1.00000000e+00 4.62962963e-03 7.71604938e-04] [1.00000000e+00 2.91545190e-03 4.16493128e-04]] Data for Fitting: [4.75969756 4.75372356 4.75176 ] Fitting Results: (array([ 4.75146378, -2.21744531, 16.23333545]), 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.12443206 1.11687698 1.11279455] Fitting Results: (array([ 1.10517191, 3.1321812 , -3.62331086]), 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: [8.38600182 8.38600182 8.38600182] Fitting Results: (array([ 8.38600182e+00, -2.60075793e-08, 4.67475994e-07]), array([], dtype=float64), 3, array([1.73207409e+00, 3.66323826e-03, 1.50029206e-05])) Fit with Size: [5 6 7] Orders: [0, 3, 5] Matrix A (Ax = y): [[1.00000000e+00 8.00000000e-03 3.20000000e-04] [1.00000000e+00 4.62962963e-03 1.28600823e-04] [1.00000000e+00 2.91545190e-03 5.94990183e-05]] Data for Fitting: [4.75969756 4.75372356 4.75176 ] Fitting Results: (array([ 4.75062869, -0.38859585, 38.05510039]), 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.12443206 1.11687698 1.11279455] Fitting Results: (array([ 1.1053583 , 2.72397858, -8.49396965]), 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: [8.38600182 8.38600182 8.38600182] Fitting Results: (array([ 8.38600182e+00, -1.85926857e-08, 1.51488937e-06]), array([], dtype=float64), 3, array([1.73207406e+00, 3.65852192e-03, 4.63568288e-06])) Fit with Size: [5 6 7] Orders: [0, 3, 6] Matrix A (Ax = y): [[1.00000000e+00 8.00000000e-03 6.40000000e-05] [1.00000000e+00 4.62962963e-03 2.14334705e-05] [1.00000000e+00 2.91545190e-03 8.49985975e-06]] Data for Fitting: [4.75969756 4.75372356 4.75176 ] Fitting Results: (array([ 4.75008492, 0.215017 , 123.32027235]), 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.12443206 1.11687698 1.11279455] Fitting Results: (array([ 1.10547967, 2.58925106, -27.52531567]), 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: [list([0, 3]) list([0, 3]) list([0, 3, 4]) list([0, 3, 5]) list([0, 3, 6])] Unrelaxed Formation Energy Fits By Size: [list([8.386001817433458, 8.386001817469117]) list([8.386001817449444]) list([8.386001817506498]) list([8.386001817496242]) list([8.386001817489564])] Formation Energy Fits By Size: [list([4.745517507252043, 4.7484204139978115]) list([4.746818873297312]) list([4.751463783973209]) list([4.750628692896055]) list([4.750084924746529])] Relaxation Volume Fits By Size: [list([1.1064991256335208, 1.1058511914129991]) list([1.10620865829578]) list([1.1051719055366451]) list([1.1053582994285363]) list([1.1054796694975624])] [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" 8.386001817469117 "source-unit" "eV" "source-std-uncert-value" 6.738901003154583e-07 } "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.65344320983 "source-unit" "angstrom" } "host-b" { "source-value" 2.65344320983 "source-unit" "angstrom" } "host-c" { "source-value" 2.65344320983 "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" [ "Mo" ] } "reservoir-cohesive-potential-energy" { "source-value" 8.386001817430609 "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.65344320983 "source-unit" "angstrom" } "reservoir-b" { "source-value" 2.65344320983 "source-unit" "angstrom" } "reservoir-c" { "source-value" 2.65344320983 "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" [ "Mo" ] } } { "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" 4.7484204139978115 "source-unit" "eV" "source-std-uncert-value" 0.003043370050007155 } "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.65344320983 "source-unit" "angstrom" } "host-b" { "source-value" 2.65344320983 "source-unit" "angstrom" } "host-c" { "source-value" 2.65344320983 "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" [ "Mo" ] } "reservoir-cohesive-potential-energy" { "source-value" 8.386001817430609 "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.65344320983 "source-unit" "angstrom" } "reservoir-b" { "source-value" 2.65344320983 "source-unit" "angstrom" } "reservoir-c" { "source-value" 2.65344320983 "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" [ "Mo" ] } } { "property-id" "tag:staff@noreply.openkim.org,2015-07-28:property/monovacancy-neutral-relaxation-volume-crystal-npt" "instance-id" 3 "relaxation-volume" { "source-value" 1.1058511914129991 "source-unit" "angstrom^3" "source-std-uncert-value" 0.0007514953874930485 } "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.65344320983 "source-unit" "angstrom" } "host-b" { "source-value" 2.65344320983 "source-unit" "angstrom" } "host-c" { "source-value" 2.65344320983 "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" [ "Mo" ] } } ]