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: Al fcc MEAM_LAMMPS_ShimKoKim_2013_AlVH__MO_344724145339_002 [4.044650822877884] CELL_SIZE_MIN: 4 CELL_SIZE_MAX: 6 Smallest System Size: 256 Largest System Size: 864 [Calculation] Supercell Size: 4 Unrelaxed Cell: [[16.17860329 0. 0. ] [ 0. 16.17860329 0. ] [ 0. 0. 16.17860329]] Unrelaxed Cell Vector: [16.178603291511536, 0.0, 16.178603291511536, 0.0, 0.0, 16.178603291511536] Unrelaxed Cell Energy: -860.1599999999875 Energy of Unrelaxed Cell With Vacancy: -860.1599999999875 ========== Loop: 0 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:51:40 -856.072814* 0.1023 FIRE: 1 16:51:40 -856.076382* 0.0963 FIRE: 2 16:51:40 -856.082872* 0.0848 FIRE: 3 16:51:40 -856.091137* 0.0686 FIRE: 4 16:51:40 -856.099765* 0.0532 FIRE: 5 16:51:40 -856.107384* 0.0396 FIRE: 6 16:51:40 -856.112938* 0.0277 FIRE: 7 16:51:40 -856.115918* 0.0167 FIRE: 8 16:51:40 -856.116470* 0.0311 FIRE: 9 16:51:40 -856.116538* 0.0306 FIRE: 10 16:51:40 -856.116669* 0.0295 FIRE: 11 16:51:40 -856.116855* 0.0281 FIRE: 12 16:51:40 -856.117084* 0.0261 FIRE: 13 16:51:40 -856.117340* 0.0238 FIRE: 14 16:51:40 -856.117608* 0.0211 FIRE: 15 16:51:40 -856.117873* 0.0181 FIRE: 16 16:51:40 -856.118143* 0.0146 FIRE: 17 16:51:40 -856.118398* 0.0105 FIRE: 18 16:51:40 -856.118610* 0.0060 FIRE: 19 16:51:40 -856.118758* 0.0050 FIRE: 20 16:51:40 -856.118836* 0.0067 FIRE: 21 16:51:40 -856.118862* 0.0079 FIRE: 22 16:51:40 -856.118868* 0.0078 FIRE: 23 16:51:40 -856.118881* 0.0077 FIRE: 24 16:51:40 -856.118900* 0.0074 FIRE: 25 16:51:40 -856.118923* 0.0070 FIRE: 26 16:51:40 -856.118949* 0.0065 FIRE: 27 16:51:40 -856.118977* 0.0060 FIRE: 28 16:51:41 -856.119006* 0.0054 FIRE: 29 16:51:41 -856.119036* 0.0046 FIRE: 30 16:51:41 -856.119066* 0.0037 FIRE: 31 16:51:41 -856.119093* 0.0027 FIRE: 32 16:51:41 -856.119114* 0.0016 FIRE: 33 16:51:41 -856.119128* 0.0021 FIRE: 34 16:51:41 -856.119135* 0.0030 FIRE: 35 16:51:41 -856.119138* 0.0035 FIRE: 36 16:51:41 -856.119139* 0.0035 FIRE: 37 16:51:41 -856.119141* 0.0034 FIRE: 38 16:51:41 -856.119144* 0.0032 FIRE: 39 16:51:41 -856.119147* 0.0030 FIRE: 40 16:51:41 -856.119151* 0.0028 FIRE: 41 16:51:41 -856.119155* 0.0025 FIRE: 42 16:51:41 -856.119159* 0.0022 FIRE: 43 16:51:41 -856.119163* 0.0019 FIRE: 44 16:51:41 -856.119167* 0.0015 FIRE: 45 16:51:41 -856.119171* 0.0010 Relaxation Completed. Steps: 45 Cell Size Relaxation... Optimization terminated successfully. Current function value: 0.679613 Iterations: 365 Function evaluations: 657 Current VFE: 0.6796130952128578 Energy of Supercell: -860.1599999999875 Unrelaxed Cell Volume: 4234.704183686962 Current Relaxed Cell Volume: 4230.132977827147 Current Relaxation Volume: 4.571205859814654 Current Cell: [[ 1.61727800e+01 0.00000000e+00 0.00000000e+00] [-2.70915515e-05 1.61727796e+01 0.00000000e+00] [ 2.49343285e-05 7.78258885e-06 1.61727798e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:51:52 -856.120387* 0.0012 FIRE: 1 16:51:52 -856.120388* 0.0012 FIRE: 2 16:51:52 -856.120389* 0.0011 FIRE: 3 16:51:52 -856.120391* 0.0010 FIRE: 4 16:51:52 -856.120393* 0.0009 Relaxation Completed. Steps: 4 Cell Size Relaxation... Optimization terminated successfully. Current function value: 0.679607 Iterations: 249 Function evaluations: 485 Current VFE: 0.6796067802476955 Energy of Supercell: -860.1599999999875 Unrelaxed Cell Volume: 4234.704183686962 Current Relaxed Cell Volume: 4230.128700727983 Current Relaxation Volume: 4.575482958978682 Current Cell: [[ 1.61727743e+01 0.00000000e+00 0.00000000e+00] [-2.23090747e-06 1.61727743e+01 0.00000000e+00] [ 3.48207218e-06 1.51308249e-05 1.61727745e+01]] ========== Loop: 2 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:51:57 -856.120393* 0.0009 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 0.679607 Iterations: 161 Function evaluations: 350 Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:52:01 -856.120393* 0.0009 FIRE: 1 16:52:01 -856.120394* 0.0009 FIRE: 2 16:52:01 -856.120394* 0.0009 FIRE: 3 16:52:01 -856.120396* 0.0008 FIRE: 4 16:52:01 -856.120397* 0.0007 FIRE: 5 16:52:01 -856.120398* 0.0006 FIRE: 6 16:52:01 -856.120400* 0.0005 FIRE: 7 16:52:01 -856.120401* 0.0003 FIRE: 8 16:52:01 -856.120403* 0.0003 FIRE: 9 16:52:01 -856.120404* 0.0002 FIRE: 10 16:52:01 -856.120405* 0.0002 FIRE: 11 16:52:01 -856.120405* 0.0002 FIRE: 12 16:52:01 -856.120405* 0.0002 FIRE: 13 16:52:01 -856.120405* 0.0002 FIRE: 14 16:52:01 -856.120405* 0.0002 FIRE: 15 16:52:01 -856.120405* 0.0002 FIRE: 16 16:52:01 -856.120405* 0.0001 FIRE: 17 16:52:01 -856.120405* 0.0001 FIRE: 18 16:52:01 -856.120405* 0.0001 FIRE: 19 16:52:01 -856.120405* 0.0001 FIRE: 20 16:52:01 -856.120405* 0.0000 Optimization terminated successfully. Current function value: 0.679595 Iterations: 194 Function evaluations: 458 --------------- Calculation Completed. Number Of Atoms in Supercell: 256 Vacancy Formation Energy (relaxed): 0.6795951215776768 Vacancy Formation Energy (unrelaxed): 0.7271861367681822 Unrelaxed Cell Volume: 4234.704183686962 Relaxed Cell Volume: 4230.128700727983 Relaxation Volume: 4.575482958978682 Relaxed Cell Vector: [16.172759455236168, -2.286380286856423e-06, 16.172759153145662, 3.5538719682870086e-06, 1.4802901963675667e-05, 16.17275911345388] Unrelaxed Cell Vector: [16.178603291511536, 0.0, 16.178603291511536, 0.0, 0.0, 16.178603291511536] Relaxed Cell: [[ 1.61727595e+01 0.00000000e+00 0.00000000e+00] [-2.28638029e-06 1.61727592e+01 0.00000000e+00] [ 3.55387197e-06 1.48029020e-05 1.61727591e+01]] Unrelaxed Cell: [[16.17860329 0. 0. ] [ 0. 16.17860329 0. ] [ 0. 0. 16.17860329]] Supercell Size: 5 Unrelaxed Cell: [[20.22325411 0. 0. ] [ 0. 20.22325411 0. ] [ 0. 0. 20.22325411]] Unrelaxed Cell Vector: [20.22325411438942, 0.0, 20.22325411438942, 0.0, 0.0, 20.22325411438942] Unrelaxed Cell Energy: -1679.9999999998963 Energy of Unrelaxed Cell With Vacancy: -1679.9999999998963 ========== Loop: 0 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:52:12 -1675.912814* 0.1023 FIRE: 1 16:52:12 -1675.916382* 0.0963 FIRE: 2 16:52:12 -1675.922872* 0.0848 FIRE: 3 16:52:12 -1675.931136* 0.0686 FIRE: 4 16:52:12 -1675.939765* 0.0532 FIRE: 5 16:52:12 -1675.947386* 0.0396 FIRE: 6 16:52:12 -1675.952945* 0.0277 FIRE: 7 16:52:12 -1675.955939* 0.0167 FIRE: 8 16:52:12 -1675.956522* 0.0310 FIRE: 9 16:52:12 -1675.956591* 0.0305 FIRE: 10 16:52:12 -1675.956725* 0.0295 FIRE: 11 16:52:12 -1675.956913* 0.0280 FIRE: 12 16:52:12 -1675.957146* 0.0260 FIRE: 13 16:52:13 -1675.957406* 0.0237 FIRE: 14 16:52:13 -1675.957679* 0.0210 FIRE: 15 16:52:13 -1675.957948* 0.0180 FIRE: 16 16:52:13 -1675.958225* 0.0144 FIRE: 17 16:52:13 -1675.958487* 0.0103 FIRE: 18 16:52:13 -1675.958711* 0.0058 FIRE: 19 16:52:13 -1675.958879* 0.0048 FIRE: 20 16:52:13 -1675.958989* 0.0066 FIRE: 21 16:52:13 -1675.959061* 0.0080 FIRE: 22 16:52:13 -1675.959124* 0.0086 FIRE: 23 16:52:13 -1675.959197* 0.0093 FIRE: 24 16:52:13 -1675.959280* 0.0084 FIRE: 25 16:52:13 -1675.959348* 0.0061 FIRE: 26 16:52:13 -1675.959367* 0.0049 FIRE: 27 16:52:13 -1675.959378* 0.0048 FIRE: 28 16:52:13 -1675.959398* 0.0045 FIRE: 29 16:52:13 -1675.959426* 0.0042 FIRE: 30 16:52:13 -1675.959457* 0.0038 FIRE: 31 16:52:13 -1675.959490* 0.0033 FIRE: 32 16:52:13 -1675.959519* 0.0029 FIRE: 33 16:52:13 -1675.959544* 0.0026 FIRE: 34 16:52:13 -1675.959564* 0.0023 FIRE: 35 16:52:13 -1675.959578* 0.0016 FIRE: 36 16:52:13 -1675.959584* 0.0011 FIRE: 37 16:52:13 -1675.959585* 0.0014 FIRE: 38 16:52:13 -1675.959585* 0.0014 FIRE: 39 16:52:13 -1675.959586* 0.0013 FIRE: 40 16:52:13 -1675.959588* 0.0013 FIRE: 41 16:52:13 -1675.959590* 0.0012 FIRE: 42 16:52:14 -1675.959592* 0.0011 FIRE: 43 16:52:14 -1675.959594* 0.0009 Relaxation Completed. Steps: 43 Cell Size Relaxation... Optimization terminated successfully. Current function value: 0.679780 Iterations: 252 Function evaluations: 500 Current VFE: 0.6797802697026327 Energy of Supercell: -1679.9999999998963 Unrelaxed Cell Volume: 8270.90660876359 Current Relaxed Cell Volume: 8266.331032298396 Current Relaxation Volume: 4.57557646519308 Current Cell: [[ 2.02195239e+01 0.00000000e+00 0.00000000e+00] [ 3.04842625e-05 2.02195247e+01 0.00000000e+00] [-4.82818084e-06 2.99368308e-05 2.02195239e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:52:34 -1675.960220* 0.0009 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 0.679780 Iterations: 207 Function evaluations: 419 Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:52:50 -1675.960220* 0.0009 FIRE: 1 16:52:50 -1675.960220* 0.0009 FIRE: 2 16:52:50 -1675.960222* 0.0008 FIRE: 3 16:52:50 -1675.960223* 0.0007 FIRE: 4 16:52:50 -1675.960225* 0.0005 FIRE: 5 16:52:50 -1675.960226* 0.0004 FIRE: 6 16:52:50 -1675.960227* 0.0003 FIRE: 7 16:52:50 -1675.960227* 0.0002 FIRE: 8 16:52:50 -1675.960227* 0.0002 FIRE: 9 16:52:50 -1675.960227* 0.0002 FIRE: 10 16:52:50 -1675.960227* 0.0002 FIRE: 11 16:52:50 -1675.960228* 0.0002 FIRE: 12 16:52:50 -1675.960228* 0.0002 FIRE: 13 16:52:50 -1675.960228* 0.0001 FIRE: 14 16:52:50 -1675.960228* 0.0001 FIRE: 15 16:52:50 -1675.960228* 0.0001 FIRE: 16 16:52:50 -1675.960228* 0.0001 FIRE: 17 16:52:50 -1675.960228* 0.0001 FIRE: 18 16:52:50 -1675.960228* 0.0001 FIRE: 19 16:52:50 -1675.960228* 0.0001 FIRE: 20 16:52:50 -1675.960228* 0.0001 Optimization terminated successfully. Current function value: 0.679772 Iterations: 246 Function evaluations: 529 --------------- Calculation Completed. Number Of Atoms in Supercell: 500 Vacancy Formation Energy (relaxed): 0.679772201541482 Vacancy Formation Energy (unrelaxed): 0.7271861367673864 Unrelaxed Cell Volume: 8270.90660876359 Relaxed Cell Volume: 8266.331032298396 Relaxation Volume: 4.57557646519308 Relaxed Cell Vector: [20.219523141237868, 9.757514684947805e-07, 20.21952332174621, -9.217162335974998e-06, 1.7857789023418858e-06, 20.21952348376532] Unrelaxed Cell Vector: [20.22325411438942, 0.0, 20.22325411438942, 0.0, 0.0, 20.22325411438942] Relaxed Cell: [[ 2.02195231e+01 0.00000000e+00 0.00000000e+00] [ 9.75751468e-07 2.02195233e+01 0.00000000e+00] [-9.21716234e-06 1.78577890e-06 2.02195235e+01]] Unrelaxed Cell: [[20.22325411 0. 0. ] [ 0. 20.22325411 0. ] [ 0. 0. 20.22325411]] Supercell Size: 6 Unrelaxed Cell: [[24.26790494 0. 0. ] [ 0. 24.26790494 0. ] [ 0. 0. 24.26790494]] Unrelaxed Cell Vector: [24.267904937267303, 0.0, 24.267904937267303, 0.0, 0.0, 24.267904937267303] Unrelaxed Cell Energy: -2903.0399999996234 Energy of Unrelaxed Cell With Vacancy: -2903.0399999996234 ========== Loop: 0 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:53:00 -2898.952814* 0.1023 FIRE: 1 16:53:00 -2898.956382* 0.0963 FIRE: 2 16:53:00 -2898.962872* 0.0848 FIRE: 3 16:53:00 -2898.971136* 0.0686 FIRE: 4 16:53:00 -2898.979765* 0.0532 FIRE: 5 16:53:00 -2898.987386* 0.0396 FIRE: 6 16:53:00 -2898.992945* 0.0277 FIRE: 7 16:53:00 -2898.995939* 0.0167 FIRE: 8 16:53:00 -2898.996521* 0.0310 FIRE: 9 16:53:00 -2898.996590* 0.0305 FIRE: 10 16:53:00 -2898.996723* 0.0295 FIRE: 11 16:53:00 -2898.996911* 0.0280 FIRE: 12 16:53:00 -2898.997143* 0.0260 FIRE: 13 16:53:00 -2898.997403* 0.0237 FIRE: 14 16:53:00 -2898.997675* 0.0210 FIRE: 15 16:53:00 -2898.997943* 0.0180 FIRE: 16 16:53:00 -2898.998220* 0.0144 FIRE: 17 16:53:00 -2898.998482* 0.0103 FIRE: 18 16:53:00 -2898.998706* 0.0059 FIRE: 19 16:53:00 -2898.998875* 0.0048 FIRE: 20 16:53:01 -2898.998988* 0.0066 FIRE: 21 16:53:01 -2898.999066* 0.0080 FIRE: 22 16:53:01 -2898.999138* 0.0086 FIRE: 23 16:53:01 -2898.999228* 0.0093 FIRE: 24 16:53:01 -2898.999337* 0.0084 FIRE: 25 16:53:01 -2898.999444* 0.0062 FIRE: 26 16:53:01 -2898.999512* 0.0036 FIRE: 27 16:53:01 -2898.999523* 0.0035 FIRE: 28 16:53:01 -2898.999543* 0.0032 FIRE: 29 16:53:01 -2898.999570* 0.0030 FIRE: 30 16:53:01 -2898.999600* 0.0030 FIRE: 31 16:53:01 -2898.999631* 0.0029 FIRE: 32 16:53:01 -2898.999659* 0.0028 FIRE: 33 16:53:01 -2898.999681* 0.0025 FIRE: 34 16:53:01 -2898.999700* 0.0021 FIRE: 35 16:53:01 -2898.999712* 0.0014 FIRE: 36 16:53:01 -2898.999720* 0.0012 FIRE: 37 16:53:01 -2898.999724* 0.0012 FIRE: 38 16:53:01 -2898.999725* 0.0012 FIRE: 39 16:53:01 -2898.999726* 0.0011 FIRE: 40 16:53:01 -2898.999728* 0.0010 FIRE: 41 16:53:01 -2898.999730* 0.0009 Relaxation Completed. Steps: 41 Cell Size Relaxation... Optimization terminated successfully. Current function value: 0.679908 Iterations: 425 Function evaluations: 754 Current VFE: 0.679908057603825 Energy of Supercell: -2903.0399999996234 Unrelaxed Cell Volume: 14292.126619943492 Current Relaxed Cell Volume: 14287.555919130045 Current Relaxation Volume: 4.57070081344682 Current Cell: [[ 2.42653174e+01 0.00000000e+00 0.00000000e+00] [ 9.85025100e-09 2.42653177e+01 0.00000000e+00] [ 5.74756833e-07 -7.52994787e-07 2.42653179e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:53:39 -2899.000092* 0.0009 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 0.679908 Iterations: 92 Function evaluations: 250 Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:53:48 -2899.000092* 0.0009 FIRE: 1 16:53:48 -2899.000093* 0.0008 FIRE: 2 16:53:48 -2899.000096* 0.0007 FIRE: 3 16:53:48 -2899.000099* 0.0007 FIRE: 4 16:53:48 -2899.000103* 0.0007 FIRE: 5 16:53:48 -2899.000106* 0.0006 FIRE: 6 16:53:48 -2899.000108* 0.0005 FIRE: 7 16:53:48 -2899.000109* 0.0004 FIRE: 8 16:53:48 -2899.000110* 0.0004 FIRE: 9 16:53:48 -2899.000111* 0.0004 FIRE: 10 16:53:48 -2899.000111* 0.0004 FIRE: 11 16:53:49 -2899.000112* 0.0007 FIRE: 12 16:53:49 -2899.000113* 0.0008 FIRE: 13 16:53:49 -2899.000114* 0.0007 FIRE: 14 16:53:49 -2899.000115* 0.0004 FIRE: 15 16:53:49 -2899.000115* 0.0004 FIRE: 16 16:53:49 -2899.000115* 0.0003 FIRE: 17 16:53:49 -2899.000115* 0.0002 FIRE: 18 16:53:49 -2899.000116* 0.0001 FIRE: 19 16:53:49 -2899.000116* 0.0001 FIRE: 20 16:53:49 -2899.000116* 0.0001 Optimization terminated successfully. Current function value: 0.679884 Iterations: 164 Function evaluations: 403 --------------- Calculation Completed. Number Of Atoms in Supercell: 864 Vacancy Formation Energy (relaxed): 0.6798842030816559 Vacancy Formation Energy (unrelaxed): 0.7271861367667043 Unrelaxed Cell Volume: 14292.126619943492 Relaxed Cell Volume: 14287.555919130045 Relaxation Volume: 4.57070081344682 Relaxed Cell Vector: [24.26531675579922, 9.581119521155392e-09, 24.265317286713184, 5.987527358044935e-07, -7.676795786339326e-07, 24.26531638176378] Unrelaxed Cell Vector: [24.267904937267303, 0.0, 24.267904937267303, 0.0, 0.0, 24.267904937267303] Relaxed Cell: [[ 2.42653168e+01 0.00000000e+00 0.00000000e+00] [ 9.58111952e-09 2.42653173e+01 0.00000000e+00] [ 5.98752736e-07 -7.67679579e-07 2.42653164e+01]] Unrelaxed Cell: [[24.26790494 0. 0. ] [ 0. 24.26790494 0. ] [ 0. 0. 24.26790494]] [Calculation Results Summary] Sizes: [4, 5, 6] Unrelaxed Formation Energy By Size: [0.7271861367681822, 0.7271861367673864, 0.7271861367667043] Formation Energy By Size: [0.6795951215776768, 0.679772201541482, 0.6798842030816559] Relaxation Volume By Size: [4.575482958978682, 4.57557646519308, 4.57070081344682] [Extrapolation] Fitting with 2 points, including orders [0, 3] Fit with data beginning 0 Fit with Size: [4 5] Orders: [0, 3] Matrix A (Ax = y): [[1. 0.015625] [1. 0.008 ]] Data for Fitting: [0.72718614 0.72718614] Fitting Results: (array([7.27186137e-01, 1.04357094e-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: [0.67959512 0.6797722 ] Fitting Results: (array([ 0.67995799, -0.0232236 ]), 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: [4.57548296 4.57557647] Fitting Results: (array([ 4.57567457, -0.01226311]), array([], dtype=float64), 2, array([1.41431223, 0.00539131])) Fit with data beginning 1 Fit with Size: [5 6] Orders: [0, 3] Matrix A (Ax = y): [[1. 0.008 ] [1. 0.00462963]] Data for Fitting: [0.72718614 0.72718614] Fitting Results: (array([7.27186137e-01, 2.02365939e-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: [0.6797722 0.6798842] Fitting Results: (array([ 0.68003805, -0.03323123]), 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.57557647 4.57070081] Fitting Results: (array([4.56400349, 1.44662195]), array([], dtype=float64), 2, array([1.41424176, 0.00238316])) Fitting with 3 points, including orders [0, 3] Fit with data beginning 0 Fit with Size: [4 5 6] Orders: [0, 3] Matrix A (Ax = y): [[1. 0.015625 ] [1. 0.008 ] [1. 0.00462963]] Data for Fitting: [0.72718614 0.72718614 0.72718614] Fitting Results: (array([7.27186137e-01, 1.29311903e-10]), array([3.33109376e-26]), 2, array([1.73212763, 0.00796621])) Fit with Size: [4 5 6] Orders: [0, 3] Matrix A (Ax = y): [[1. 0.015625 ] [1. 0.008 ] [1. 0.00462963]] Data for Fitting: [0.67959512 0.6797722 0.6798842 ] Fitting Results: (array([ 0.6799932 , -0.02576852]), array([3.47403451e-10]), 2, array([1.73212763, 0.00796621])) Fit with Size: [4 5 6] Orders: [0, 3] Matrix A (Ax = y): [[1. 0.015625 ] [1. 0.008 ] [1. 0.00462963]] Data for Fitting: [4.57548296 4.57557647 4.57070081] Fitting Results: (array([4.5705415 , 0.35872835]), array([7.38268411e-06]), 2, array([1.73212763, 0.00796621])) Fitting with 3 points, including orders [0, 3, 4] Fit with data beginning 0 Fit with Size: [4 5 6] Orders: [0, 3, 4] Matrix A (Ax = y): [[1.00000000e+00 1.56250000e-02 3.90625000e-03] [1.00000000e+00 8.00000000e-03 1.60000000e-03] [1.00000000e+00 4.62962963e-03 7.71604938e-04]] Data for Fitting: [0.72718614 0.72718614 0.72718614] Fitting Results: (array([ 7.27186137e-01, 6.27401083e-10, -1.72919065e-09]), array([], dtype=float64), 3, array([1.73213142e+00, 8.29013951e-03, 1.01423660e-04])) Fit with Size: [4 5 6] Orders: [0, 3, 4] Matrix A (Ax = y): [[1.00000000e+00 1.56250000e-02 3.90625000e-03] [1.00000000e+00 8.00000000e-03 1.60000000e-03] [1.00000000e+00 4.62962963e-03 7.71604938e-04]] Data for Fitting: [0.67959512 0.6797722 0.6798842 ] Fitting Results: (array([ 0.68010274, -0.07663489, 0.17659017]), array([], dtype=float64), 3, array([1.73213142e+00, 8.29013951e-03, 1.01423660e-04])) Fit with Size: [4 5 6] Orders: [0, 3, 4] Matrix A (Ax = y): [[1.00000000e+00 1.56250000e-02 3.90625000e-03] [1.00000000e+00 8.00000000e-03 1.60000000e-03] [1.00000000e+00 4.62962963e-03 7.71604938e-04]] Data for Fitting: [4.57548296 4.57557647 4.57070081] Fitting Results: (array([ 4.55457387, 7.77389339, -25.74284805]), array([], dtype=float64), 3, array([1.73213142e+00, 8.29013951e-03, 1.01423660e-04])) Fitting with 3 points, including orders [0, 3, 5] Fit with data beginning 0 Fit with Size: [4 5 6] Orders: [0, 3, 5] Matrix A (Ax = y): [[1.00000000e+00 1.56250000e-02 9.76562500e-04] [1.00000000e+00 8.00000000e-03 3.20000000e-04] [1.00000000e+00 4.62962963e-03 1.28600823e-04]] Data for Fitting: [0.72718614 0.72718614 0.72718614] Fitting Results: (array([ 7.27186137e-01, 3.92206072e-10, -3.34253623e-09]), array([], dtype=float64), 3, array([1.73212782e+00, 7.99080777e-03, 5.44349833e-05])) Fit with Size: [4 5 6] Orders: [0, 3, 5] Matrix A (Ax = y): [[1.00000000e+00 1.56250000e-02 9.76562500e-04] [1.00000000e+00 8.00000000e-03 3.20000000e-04] [1.00000000e+00 4.62962963e-03 1.28600823e-04]] Data for Fitting: [0.67959512 0.6797722 0.6798842 ] Fitting Results: (array([ 0.6800839 , -0.05261607, 0.34134988]), array([], dtype=float64), 3, array([1.73212782e+00, 7.99080777e-03, 5.44349833e-05])) Fit with Size: [4 5 6] Orders: [0, 3, 5] Matrix A (Ax = y): [[1.00000000e+00 1.56250000e-02 9.76562500e-04] [1.00000000e+00 8.00000000e-03 3.20000000e-04] [1.00000000e+00 4.62962963e-03 1.28600823e-04]] Data for Fitting: [4.57548296 4.57557647 4.57070081] Fitting Results: (array([ 4.55732007, 4.27249257, -49.76108457]), array([], dtype=float64), 3, array([1.73212782e+00, 7.99080777e-03, 5.44349833e-05])) Fitting with 3 points, including orders [0, 3, 6] Fit with data beginning 0 Fit with Size: [4 5 6] Orders: [0, 3, 6] Matrix A (Ax = y): [[1.00000000e+00 1.56250000e-02 2.44140625e-04] [1.00000000e+00 8.00000000e-03 6.40000000e-05] [1.00000000e+00 4.62962963e-03 2.14334705e-05]] Data for Fitting: [0.72718614 0.72718614 0.72718614] Fitting Results: (array([ 7.27186137e-01, 3.14948598e-10, -8.91245805e-09]), array([], dtype=float64), 3, array([1.73212764e+00, 7.96793362e-03, 2.04739548e-05])) Fit with Size: [4 5 6] Orders: [0, 3, 6] Matrix A (Ax = y): [[1.00000000e+00 1.56250000e-02 2.44140625e-04] [1.00000000e+00 8.00000000e-03 6.40000000e-05] [1.00000000e+00 4.62962963e-03 2.14334705e-05]] Data for Fitting: [0.67959512 0.6797722 0.6798842 ] Fitting Results: (array([ 0.68007176, -0.0447263 , 0.9101671 ]), array([], dtype=float64), 3, array([1.73212764e+00, 7.96793362e-03, 2.04739548e-05])) Fit with Size: [4 5 6] Orders: [0, 3, 6] Matrix A (Ax = y): [[1.00000000e+00 1.56250000e-02 2.44140625e-04] [1.00000000e+00 8.00000000e-03 6.40000000e-05] [1.00000000e+00 4.62962963e-03 2.14334705e-05]] Data for Fitting: [4.57548296 4.57557647 4.57070081] Fitting Results: (array([ 4.55908935, 3.12234339, -132.68175674]), array([], dtype=float64), 3, array([1.73212764e+00, 7.96793362e-03, 2.04739548e-05])) [Fitting Results Summary] Sizes: [4, 5, 6] 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: [[0.7271861367665513, 0.7271861367657674], [0.7271861367662064], [0.7271861367651342], [0.7271861367653187], [0.727186136765437]] Formation Energy Fits By Size: [[0.6799579903559658, 0.6800380513511255], [0.6799932020635087], [0.6801027363935722], [0.6800838980997939], [0.6800717612438238]] Relaxation Volume Fits By Size: [[4.575674570073756, 4.56400348961954], [4.5705415003331025], [4.554573874947916], [4.557320071667943], [4.5590893504809165]] [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" 0.7271861367657674 "source-unit" "eV" "source-std-uncert-value" 2.3854522169131097e-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" [ "fcc" ] } "host-a" { "source-value" 4.044650822877884 "source-unit" "angstrom" } "host-b" { "source-value" 4.044650822877884 "source-unit" "angstrom" } "host-c" { "source-value" 4.044650822877884 "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" "Fm-3m" } "host-wyckoff-multiplicity-and-letter" { "source-value" [ "4a" ] } "host-wyckoff-coordinates" { "source-value" [ [ 0.0 0.0 0.0 ] ] } "host-wyckoff-species" { "source-value" [ "Al" ] } "reservoir-cohesive-potential-energy" { "source-value" 3.3599999999999453 "source-unit" "eV" } "reservoir-short-name" { "source-value" [ "fcc" ] } "reservoir-cauchy-stress" { "source-value" [ 0.0 0.0 0.0 0.0 0.0 0.0 ] "source-unit" "GPa" } "reservoir-a" { "source-value" 4.044650822877884 "source-unit" "angstrom" } "reservoir-b" { "source-value" 4.044650822877884 "source-unit" "angstrom" } "reservoir-c" { "source-value" 4.044650822877884 "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" "Fm-3m" } "reservoir-wyckoff-multiplicity-and-letter" { "source-value" [ "4a" ] } "reservoir-wyckoff-coordinates" { "source-value" [ [ 0.0 0.0 0.0 ] ] } "reservoir-wyckoff-species" { "source-value" [ "Al" ] } } { "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" 0.6800380513511255 "source-unit" "eV" "source-std-uncert-value" 6.894340392127255e-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" [ "fcc" ] } "host-a" { "source-value" 4.044650822877884 "source-unit" "angstrom" } "host-b" { "source-value" 4.044650822877884 "source-unit" "angstrom" } "host-c" { "source-value" 4.044650822877884 "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" "Fm-3m" } "host-wyckoff-multiplicity-and-letter" { "source-value" [ "4a" ] } "host-wyckoff-coordinates" { "source-value" [ [ 0.0 0.0 0.0 ] ] } "host-wyckoff-species" { "source-value" [ "Al" ] } "reservoir-cohesive-potential-energy" { "source-value" 3.3599999999999453 "source-unit" "eV" } "reservoir-short-name" { "source-value" [ "fcc" ] } "reservoir-cauchy-stress" { "source-value" [ 0.0 0.0 0.0 0.0 0.0 0.0 ] "source-unit" "GPa" } "reservoir-a" { "source-value" 4.044650822877884 "source-unit" "angstrom" } "reservoir-b" { "source-value" 4.044650822877884 "source-unit" "angstrom" } "reservoir-c" { "source-value" 4.044650822877884 "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" "Fm-3m" } "reservoir-wyckoff-multiplicity-and-letter" { "source-value" [ "4a" ] } "reservoir-wyckoff-coordinates" { "source-value" [ [ 0.0 0.0 0.0 ] ] } "reservoir-wyckoff-species" { "source-value" [ "Al" ] } } { "property-id" "tag:staff@noreply.openkim.org,2015-07-28:property/monovacancy-neutral-relaxation-volume-crystal-npt" "instance-id" 3 "relaxation-volume" { "source-value" 4.56400348961954 "source-unit" "angstrom^3" "source-std-uncert-value" 0.009548437661256952 } "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" [ "fcc" ] } "host-a" { "source-value" 4.044650822877884 "source-unit" "angstrom" } "host-b" { "source-value" 4.044650822877884 "source-unit" "angstrom" } "host-c" { "source-value" 4.044650822877884 "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" "Fm-3m" } "host-wyckoff-multiplicity-and-letter" { "source-value" [ "4a" ] } "host-wyckoff-coordinates" { "source-value" [ [ 0.0 0.0 0.0 ] ] } "host-wyckoff-species" { "source-value" [ "Al" ] } } ]