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: Pb fcc LJ_ElliottAkerson_2015_Universal__MO_959249795837_003 [4.015757858753205] CELL_SIZE_MIN: 4 CELL_SIZE_MAX: 6 Smallest System Size: 256 Largest System Size: 864 [Calculation] Supercell Size: 4 Unrelaxed Cell: [[16.06303144 0. 0. ] [ 0. 16.06303144 0. ] [ 0. 0. 16.06303144]] Unrelaxed Cell Vector: [16.06303143501282, 0.0, 16.06303143501282, 0.0, 0.0, 16.06303143501282] Unrelaxed Cell Energy: -1830.655422713052 Energy of Unrelaxed Cell With Vacancy: -1830.655422713052 ========== Loop: 0 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 13:40:19 -1816.353427* 0.7831 FIRE: 1 13:40:19 -1816.414172* 0.3938 FIRE: 2 13:40:19 -1816.433800* 0.1808 FIRE: 3 13:40:19 -1816.437260* 0.1555 FIRE: 4 13:40:20 -1816.442047* 0.1090 FIRE: 5 13:40:20 -1816.445363* 0.0494 FIRE: 6 13:40:20 -1816.445841* 0.0822 FIRE: 7 13:40:20 -1816.445964* 0.0794 FIRE: 8 13:40:20 -1816.446192* 0.0737 FIRE: 9 13:40:20 -1816.446492* 0.0656 FIRE: 10 13:40:20 -1816.446821* 0.0554 FIRE: 11 13:40:20 -1816.447138* 0.0435 FIRE: 12 13:40:20 -1816.447408* 0.0307 FIRE: 13 13:40:20 -1816.447613* 0.0206 FIRE: 14 13:40:20 -1816.447770* 0.0253 FIRE: 15 13:40:20 -1816.447892* 0.0284 FIRE: 16 13:40:20 -1816.448012* 0.0297 FIRE: 17 13:40:20 -1816.448159* 0.0270 FIRE: 18 13:40:20 -1816.448328* 0.0219 FIRE: 19 13:40:20 -1816.448457* 0.0160 FIRE: 20 13:40:20 -1816.448448* 0.0097 FIRE: 21 13:40:20 -1816.448452* 0.0095 FIRE: 22 13:40:20 -1816.448459* 0.0090 FIRE: 23 13:40:20 -1816.448469* 0.0084 FIRE: 24 13:40:20 -1816.448481* 0.0075 FIRE: 25 13:40:20 -1816.448494* 0.0065 FIRE: 26 13:40:20 -1816.448506* 0.0054 FIRE: 27 13:40:20 -1816.448518* 0.0047 FIRE: 28 13:40:20 -1816.448529* 0.0043 FIRE: 29 13:40:20 -1816.448540* 0.0039 FIRE: 30 13:40:20 -1816.448548* 0.0035 FIRE: 31 13:40:20 -1816.448554* 0.0027 FIRE: 32 13:40:20 -1816.448557* 0.0018 FIRE: 33 13:40:20 -1816.448557* 0.0018 FIRE: 34 13:40:20 -1816.448557* 0.0017 FIRE: 35 13:40:20 -1816.448558* 0.0016 FIRE: 36 13:40:20 -1816.448558* 0.0015 FIRE: 37 13:40:20 -1816.448558* 0.0013 FIRE: 38 13:40:20 -1816.448559* 0.0012 FIRE: 39 13:40:20 -1816.448559* 0.0010 FIRE: 40 13:40:20 -1816.448559* 0.0008 Relaxation Completed. Steps: 40 Cell Size Relaxation... Optimization terminated successfully. Current function value: 7.055704 Iterations: 356 Function evaluations: 630 Current VFE: 7.055703950006318 Energy of Supercell: -1830.655422713052 Unrelaxed Cell Volume: 4144.599094677714 Current Relaxed Cell Volume: 4143.986164669759 Current Relaxation Volume: 0.6129300079546738 Current Cell: [[ 1.60622391e+01 0.00000000e+00 0.00000000e+00] [-4.75031844e-08 1.60622398e+01 0.00000000e+00] [ 2.44333291e-07 -2.28023416e-07 1.60622398e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 13:40:51 -1816.448721* 0.0017 FIRE: 1 13:40:51 -1816.448722* 0.0010 Relaxation Completed. Steps: 1 Cell Size Relaxation... Optimization terminated successfully. Current function value: 7.055703 Iterations: 136 Function evaluations: 303 Current VFE: 7.055703432484961 Energy of Supercell: -1830.655422713052 Unrelaxed Cell Volume: 4144.599094677714 Current Relaxed Cell Volume: 4143.985650463128 Current Relaxation Volume: 0.6134442145857975 Current Cell: [[ 1.60622390e+01 0.00000000e+00 0.00000000e+00] [-4.88869016e-08 1.60622384e+01 0.00000000e+00] [ 2.52291274e-07 -2.19676699e-07 1.60622392e+01]] ========== Loop: 2 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 13:41:06 -1816.448722* 0.0010 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 7.055703 Iterations: 112 Function evaluations: 268 Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 13:41:18 -1816.448722* 0.0010 FIRE: 1 13:41:18 -1816.448722* 0.0006 FIRE: 2 13:41:18 -1816.448722* 0.0002 FIRE: 3 13:41:18 -1816.448722* 0.0004 FIRE: 4 13:41:18 -1816.448722* 0.0004 FIRE: 5 13:41:18 -1816.448722* 0.0003 FIRE: 6 13:41:18 -1816.448722* 0.0002 FIRE: 7 13:41:18 -1816.448722* 0.0001 FIRE: 8 13:41:18 -1816.448722* 0.0001 FIRE: 9 13:41:18 -1816.448722* 0.0001 FIRE: 10 13:41:18 -1816.448722* 0.0001 FIRE: 11 13:41:18 -1816.448722* 0.0001 FIRE: 12 13:41:18 -1816.448722* 0.0001 FIRE: 13 13:41:18 -1816.448722* 0.0000 FIRE: 14 13:41:18 -1816.448722* 0.0000 FIRE: 15 13:41:18 -1816.448722* 0.0000 FIRE: 16 13:41:18 -1816.448722* 0.0000 FIRE: 17 13:41:18 -1816.448722* 0.0000 FIRE: 18 13:41:18 -1816.448722* 0.0000 FIRE: 19 13:41:18 -1816.448722* 0.0000 FIRE: 20 13:41:18 -1816.448722* 0.0000 Optimization terminated successfully. Current function value: 7.055703 Iterations: 182 Function evaluations: 429 --------------- Calculation Completed. Number Of Atoms in Supercell: 256 Vacancy Formation Energy (relaxed): 7.055702931178303 Vacancy Formation Energy (unrelaxed): 7.150997744980032 Unrelaxed Cell Volume: 4144.599094677714 Relaxed Cell Volume: 4143.985650463128 Relaxation Volume: 0.6134442145857975 Relaxed Cell Vector: [16.062237735941054, -4.9055742520122404e-08, 16.062237803150445, 2.6046426142608266e-07, -2.1395143197009538e-07, 16.06223787382732] Unrelaxed Cell Vector: [16.06303143501282, 0.0, 16.06303143501282, 0.0, 0.0, 16.06303143501282] Relaxed Cell: [[ 1.60622377e+01 0.00000000e+00 0.00000000e+00] [-4.90557425e-08 1.60622378e+01 0.00000000e+00] [ 2.60464261e-07 -2.13951432e-07 1.60622379e+01]] Unrelaxed Cell: [[16.06303144 0. 0. ] [ 0. 16.06303144 0. ] [ 0. 0. 16.06303144]] Supercell Size: 5 Unrelaxed Cell: [[20.07878929 0. 0. ] [ 0. 20.07878929 0. ] [ 0. 0. 20.07878929]] Unrelaxed Cell Vector: [20.078789293766025, 0.0, 20.078789293766025, 0.0, 0.0, 20.078789293766025] Unrelaxed Cell Energy: -3575.498872484064 Energy of Unrelaxed Cell With Vacancy: -3575.498872484064 ========== Loop: 0 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 13:41:40 -3561.196877* 0.7831 FIRE: 1 13:41:40 -3561.257630* 0.3939 FIRE: 2 13:41:40 -3561.277287* 0.1809 FIRE: 3 13:41:40 -3561.280752* 0.1555 FIRE: 4 13:41:40 -3561.285550* 0.1090 FIRE: 5 13:41:40 -3561.288879* 0.0493 FIRE: 6 13:41:40 -3561.289355* 0.0822 FIRE: 7 13:41:40 -3561.289476* 0.0794 FIRE: 8 13:41:40 -3561.289700* 0.0738 FIRE: 9 13:41:40 -3561.289993* 0.0656 FIRE: 10 13:41:40 -3561.290317* 0.0554 FIRE: 11 13:41:40 -3561.290628* 0.0435 FIRE: 12 13:41:40 -3561.290896* 0.0306 FIRE: 13 13:41:40 -3561.291105* 0.0201 FIRE: 14 13:41:40 -3561.291274* 0.0247 FIRE: 15 13:41:40 -3561.291417* 0.0279 FIRE: 16 13:41:41 -3561.291568* 0.0294 FIRE: 17 13:41:41 -3561.291756* 0.0269 FIRE: 18 13:41:41 -3561.291972* 0.0224 FIRE: 19 13:41:41 -3561.292152* 0.0163 FIRE: 20 13:41:41 -3561.292199* 0.0077 FIRE: 21 13:41:41 -3561.292203* 0.0074 FIRE: 22 13:41:41 -3561.292211* 0.0070 FIRE: 23 13:41:41 -3561.292222* 0.0063 FIRE: 24 13:41:41 -3561.292235* 0.0055 FIRE: 25 13:41:41 -3561.292248* 0.0046 FIRE: 26 13:41:41 -3561.292261* 0.0044 FIRE: 27 13:41:41 -3561.292273* 0.0041 FIRE: 28 13:41:41 -3561.292284* 0.0036 FIRE: 29 13:41:41 -3561.292294* 0.0034 FIRE: 30 13:41:41 -3561.292304* 0.0031 FIRE: 31 13:41:41 -3561.292313* 0.0025 FIRE: 32 13:41:41 -3561.292320* 0.0020 FIRE: 33 13:41:41 -3561.292323* 0.0030 FIRE: 34 13:41:41 -3561.292322* 0.0041 FIRE: 35 13:41:41 -3561.292323* 0.0040 FIRE: 36 13:41:41 -3561.292323* 0.0039 FIRE: 37 13:41:41 -3561.292324* 0.0036 FIRE: 38 13:41:41 -3561.292326* 0.0032 FIRE: 39 13:41:41 -3561.292327* 0.0028 FIRE: 40 13:41:42 -3561.292328* 0.0024 FIRE: 41 13:41:42 -3561.292330* 0.0019 FIRE: 42 13:41:42 -3561.292331* 0.0017 FIRE: 43 13:41:42 -3561.292332* 0.0016 FIRE: 44 13:41:42 -3561.292333* 0.0013 FIRE: 45 13:41:42 -3561.292334* 0.0009 Relaxation Completed. Steps: 45 Cell Size Relaxation... Optimization terminated successfully. Current function value: 7.055457 Iterations: 264 Function evaluations: 508 Current VFE: 7.055457154669057 Energy of Supercell: -3575.498872484064 Unrelaxed Cell Volume: 8094.920106792415 Current Relaxed Cell Volume: 8094.307024886691 Current Relaxation Volume: 0.6130819057243571 Current Cell: [[ 2.00782822e+01 0.00000000e+00 0.00000000e+00] [ 5.46472301e-07 2.00782821e+01 0.00000000e+00] [ 5.00864488e-07 -2.92472326e-07 2.00782828e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 13:42:26 -3561.292418* 0.0009 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 7.055457 Iterations: 117 Function evaluations: 282 Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 13:42:47 -3561.292418* 0.0009 FIRE: 1 13:42:47 -3561.292418* 0.0005 FIRE: 2 13:42:47 -3561.292419* 0.0005 FIRE: 3 13:42:47 -3561.292419* 0.0004 FIRE: 4 13:42:47 -3561.292420* 0.0002 FIRE: 5 13:42:47 -3561.292420* 0.0004 FIRE: 6 13:42:47 -3561.292420* 0.0004 FIRE: 7 13:42:47 -3561.292420* 0.0003 FIRE: 8 13:42:47 -3561.292420* 0.0001 FIRE: 9 13:42:48 -3561.292420* 0.0002 FIRE: 10 13:42:48 -3561.292420* 0.0002 FIRE: 11 13:42:48 -3561.292420* 0.0002 FIRE: 12 13:42:48 -3561.292420* 0.0002 FIRE: 13 13:42:48 -3561.292420* 0.0001 FIRE: 14 13:42:48 -3561.292420* 0.0001 FIRE: 15 13:42:48 -3561.292420* 0.0001 FIRE: 16 13:42:48 -3561.292420* 0.0001 FIRE: 17 13:42:48 -3561.292420* 0.0001 FIRE: 18 13:42:48 -3561.292420* 0.0001 FIRE: 19 13:42:48 -3561.292420* 0.0001 FIRE: 20 13:42:48 -3561.292420* 0.0001 Optimization terminated successfully. Current function value: 7.055455 Iterations: 170 Function evaluations: 414 --------------- Calculation Completed. Number Of Atoms in Supercell: 500 Vacancy Formation Energy (relaxed): 7.055455025647916 Vacancy Formation Energy (unrelaxed): 7.1509977449677535 Unrelaxed Cell Volume: 8094.920106792415 Relaxed Cell Volume: 8094.307024886691 Relaxation Volume: 0.6130819057243571 Relaxed Cell Vector: [20.07828152654301, 5.580058402809814e-07, 20.07828095225328, 5.045309579600651e-07, -2.9021010013447346e-07, 20.07828187635653] Unrelaxed Cell Vector: [20.078789293766025, 0.0, 20.078789293766025, 0.0, 0.0, 20.078789293766025] Relaxed Cell: [[ 2.00782815e+01 0.00000000e+00 0.00000000e+00] [ 5.58005840e-07 2.00782810e+01 0.00000000e+00] [ 5.04530958e-07 -2.90210100e-07 2.00782819e+01]] Unrelaxed Cell: [[20.07878929 0. 0. ] [ 0. 20.07878929 0. ] [ 0. 0. 20.07878929]] Supercell Size: 6 Unrelaxed Cell: [[24.09454715 0. 0. ] [ 0. 24.09454715 0. ] [ 0. 0. 24.09454715]] Unrelaxed Cell Vector: [24.094547152519233, 0.0, 24.094547152519233, 0.0, 0.0, 24.094547152519233] Unrelaxed Cell Energy: -6178.462051653514 Energy of Unrelaxed Cell With Vacancy: -6178.462051653514 ========== Loop: 0 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 13:43:20 -6164.160056* 0.7831 FIRE: 1 13:43:20 -6164.220809* 0.3939 FIRE: 2 13:43:20 -6164.240466* 0.1809 FIRE: 3 13:43:20 -6164.243931* 0.1555 FIRE: 4 13:43:20 -6164.248729* 0.1091 FIRE: 5 13:43:20 -6164.252059* 0.0493 FIRE: 6 13:43:20 -6164.252538* 0.0822 FIRE: 7 13:43:20 -6164.252659* 0.0794 FIRE: 8 13:43:20 -6164.252883* 0.0738 FIRE: 9 13:43:20 -6164.253178* 0.0656 FIRE: 10 13:43:20 -6164.253502* 0.0554 FIRE: 11 13:43:20 -6164.253814* 0.0435 FIRE: 12 13:43:21 -6164.254082* 0.0306 FIRE: 13 13:43:21 -6164.254291* 0.0201 FIRE: 14 13:43:21 -6164.254458* 0.0247 FIRE: 15 13:43:21 -6164.254599* 0.0279 FIRE: 16 13:43:21 -6164.254746* 0.0294 FIRE: 17 13:43:21 -6164.254930* 0.0268 FIRE: 18 13:43:21 -6164.255144* 0.0224 FIRE: 19 13:43:21 -6164.255327* 0.0164 FIRE: 20 13:43:21 -6164.255388* 0.0077 FIRE: 21 13:43:21 -6164.255393* 0.0075 FIRE: 22 13:43:21 -6164.255402* 0.0071 FIRE: 23 13:43:21 -6164.255415* 0.0064 FIRE: 24 13:43:21 -6164.255430* 0.0056 FIRE: 25 13:43:21 -6164.255446* 0.0048 FIRE: 26 13:43:21 -6164.255462* 0.0046 FIRE: 27 13:43:21 -6164.255478* 0.0042 FIRE: 28 13:43:21 -6164.255494* 0.0037 FIRE: 29 13:43:21 -6164.255510* 0.0037 FIRE: 30 13:43:21 -6164.255526* 0.0035 FIRE: 31 13:43:21 -6164.255541* 0.0030 FIRE: 32 13:43:21 -6164.255555* 0.0021 FIRE: 33 13:43:21 -6164.255565* 0.0030 FIRE: 34 13:43:21 -6164.255569* 0.0042 FIRE: 35 13:43:21 -6164.255569* 0.0041 FIRE: 36 13:43:21 -6164.255570* 0.0039 FIRE: 37 13:43:21 -6164.255571* 0.0036 FIRE: 38 13:43:21 -6164.255573* 0.0033 FIRE: 39 13:43:21 -6164.255575* 0.0030 FIRE: 40 13:43:22 -6164.255576* 0.0025 FIRE: 41 13:43:22 -6164.255578* 0.0021 FIRE: 42 13:43:22 -6164.255580* 0.0020 FIRE: 43 13:43:22 -6164.255582* 0.0018 FIRE: 44 13:43:22 -6164.255585* 0.0015 FIRE: 45 13:43:22 -6164.255587* 0.0011 FIRE: 46 13:43:22 -6164.255589* 0.0009 Relaxation Completed. Steps: 46 Cell Size Relaxation... Optimization terminated successfully. Current function value: 7.055417 Iterations: 223 Function evaluations: 447 Current VFE: 7.055416864982362 Energy of Supercell: -6178.462051653514 Unrelaxed Cell Volume: 13988.021944537306 Current Relaxed Cell Volume: 13987.409283487234 Current Relaxation Volume: 0.6126610500723473 Current Cell: [[ 2.40941951e+01 0.00000000e+00 0.00000000e+00] [ 9.39850785e-07 2.40941948e+01 0.00000000e+00] [-1.09916581e-07 -1.11783508e-06 2.40941962e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 13:44:21 -6164.255637* 0.0009 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 7.055417 Iterations: 111 Function evaluations: 266 Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 13:44:55 -6164.255637* 0.0009 FIRE: 1 13:44:55 -6164.255638* 0.0007 FIRE: 2 13:44:55 -6164.255640* 0.0006 FIRE: 3 13:44:55 -6164.255642* 0.0005 FIRE: 4 13:44:55 -6164.255644* 0.0003 FIRE: 5 13:44:55 -6164.255645* 0.0004 FIRE: 6 13:44:55 -6164.255646* 0.0002 FIRE: 7 13:44:55 -6164.255645* 0.0003 FIRE: 8 13:44:55 -6164.255645* 0.0003 FIRE: 9 13:44:55 -6164.255645* 0.0002 FIRE: 10 13:44:55 -6164.255645* 0.0001 FIRE: 11 13:44:55 -6164.255645* 0.0002 FIRE: 12 13:44:55 -6164.255645* 0.0002 FIRE: 13 13:44:56 -6164.255645* 0.0002 FIRE: 14 13:44:56 -6164.255645* 0.0001 FIRE: 15 13:44:56 -6164.255646* 0.0001 FIRE: 16 13:44:56 -6164.255646* 0.0001 FIRE: 17 13:44:56 -6164.255646* 0.0001 FIRE: 18 13:44:56 -6164.255646* 0.0001 FIRE: 19 13:44:56 -6164.255646* 0.0001 FIRE: 20 13:44:56 -6164.255646* 0.0001 Optimization terminated successfully. Current function value: 7.055408 Iterations: 181 Function evaluations: 429 --------------- Calculation Completed. Number Of Atoms in Supercell: 864 Vacancy Formation Energy (relaxed): 7.0554083140350485 Vacancy Formation Energy (unrelaxed): 7.15099774495684 Unrelaxed Cell Volume: 13988.021944537306 Relaxed Cell Volume: 13987.409283487234 Relaxation Volume: 0.6126610500723473 Relaxed Cell Vector: [24.09419461908191, 9.382736534129236e-07, 24.094193520028476, -1.1328569328550029e-07, -1.1153025357909038e-06, 24.09419453686455] Unrelaxed Cell Vector: [24.094547152519233, 0.0, 24.094547152519233, 0.0, 0.0, 24.094547152519233] Relaxed Cell: [[ 2.40941946e+01 0.00000000e+00 0.00000000e+00] [ 9.38273653e-07 2.40941935e+01 0.00000000e+00] [-1.13285693e-07 -1.11530254e-06 2.40941945e+01]] Unrelaxed Cell: [[24.09454715 0. 0. ] [ 0. 24.09454715 0. ] [ 0. 0. 24.09454715]] [Calculation Results Summary] Sizes: [4, 5, 6] Unrelaxed Formation Energy By Size: [7.150997744980032, 7.1509977449677535, 7.15099774495684] Formation Energy By Size: [7.055702931178303, 7.055455025647916, 7.0554083140350485] Relaxation Volume By Size: [0.6134442145857975, 0.6130819057243571, 0.6126610500723473] [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: [7.15099774 7.15099774] Fitting Results: (array([7.15099774e+00, 1.61008088e-09]), 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: [7.05570293 7.05545503] Fitting Results: (array([7.05519493, 0.0325122 ]), 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.61344421 0.61308191] Fitting Results: (array([0.61270178, 0.04751592]), 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: [7.15099774 7.15099774] Fitting Results: (array([7.15099774e+00, 3.23860038e-09]), array([], dtype=float64), 2, array([1.41424176, 0.00238316])) Fit with Size: [5 6] Orders: [0, 3] Matrix A (Ax = y): [[1. 0.008 ] [1. 0.00462963]] Data for Fitting: [7.05545503 7.05540831] Fitting Results: (array([7.05534415, 0.01385949]), 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.61308191 0.61266105] Fitting Results: (array([0.61208295, 0.12486926]), 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: [7.15099774 7.15099774 7.15099774] Fitting Results: (array([7.15099774e+00, 2.02445784e-09]), array([9.18956509e-24]), 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: [7.05570293 7.05545503 7.05540831] Fitting Results: (array([7.05526056, 0.02776885]), array([1.2068577e-09]), 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.61344421 0.61308191 0.61266105] Fitting Results: (array([0.61242961, 0.06718671]), array([2.07553619e-08]), 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: [7.15099774 7.15099774 7.15099774] Fitting Results: (array([ 7.15099774e+00, 1.02974223e-08, -2.87208262e-08]), 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: [7.05570293 7.05545503 7.05540831] Fitting Results: (array([ 7.05546471, -0.06703843, 0.32913759]), 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.61344421 0.61308191 0.61266105] Fitting Results: (array([ 0.61158297, 0.46035532, -1.36494328]), 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: [7.15099774 7.15099774 7.15099774] Fitting Results: (array([ 7.15099774e+00, 6.39097327e-09, -5.55175348e-08]), 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: [7.05570293 7.05545503 7.05540831] Fitting Results: (array([ 7.0554296 , -0.02227094, 0.636225 ]), 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.61344421 0.61308191 0.61266105] Fitting Results: (array([ 0.61172858, 0.27470323, -2.63844381]), 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: [7.15099774 7.15099774 7.15099774] Fitting Results: (array([ 7.15099774e+00, 5.10777273e-09, -1.48030617e-07]), 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: [7.05570293 7.05545503 7.05540831] Fitting Results: (array([ 7.05540698, -0.0075656 , 1.696415 ]), 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.61344421 0.61308191 0.61266105] Fitting Results: (array([ 0.61182239, 0.21371975, -7.035083 ]), 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: [[7.150997744954869, 7.150997744941847], [7.150997744949142], [7.150997744931332], [7.150997744934395], [7.150997744936364]] Formation Energy Fits By Size: [[7.055194928042261, 7.055344149731661], [7.055260557384936], [7.055464712950909], [7.055429601196209], [7.055406979916667]] Relaxation Volume Fits By Size: [[0.612701778394321, 0.6120829516492571], [0.6124296115754684], [0.6115829724247683], [0.6117285819066175], [0.611822393019756]] [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" 7.150997744941847 "source-unit" "eV" "source-std-uncert-value" 8.5509473137855e-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" [ "fcc" ] } "host-a" { "source-value" 4.015757858753205 "source-unit" "angstrom" } "host-b" { "source-value" 4.015757858753205 "source-unit" "angstrom" } "host-c" { "source-value" 4.015757858753205 "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" [ "Pb" ] } "reservoir-cohesive-potential-energy" { "source-value" 7.150997744971724 "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.015757858753205 "source-unit" "angstrom" } "reservoir-b" { "source-value" 4.015757858753205 "source-unit" "angstrom" } "reservoir-c" { "source-value" 4.015757858753205 "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" [ "Pb" ] } } { "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" 7.055344149731661 "source-unit" "eV" "source-std-uncert-value" 0.00012086607685951736 } "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.015757858753205 "source-unit" "angstrom" } "host-b" { "source-value" 4.015757858753205 "source-unit" "angstrom" } "host-c" { "source-value" 4.015757858753205 "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" [ "Pb" ] } "reservoir-cohesive-potential-energy" { "source-value" 7.150997744971724 "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.015757858753205 "source-unit" "angstrom" } "reservoir-b" { "source-value" 4.015757858753205 "source-unit" "angstrom" } "reservoir-c" { "source-value" 4.015757858753205 "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" [ "Pb" ] } } { "property-id" "tag:staff@noreply.openkim.org,2015-07-28:property/monovacancy-neutral-relaxation-volume-crystal-npt" "instance-id" 3 "relaxation-volume" { "source-value" 0.6120829516492571 "source-unit" "angstrom^3" "source-std-uncert-value" 0.0020641016892911277 } "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.015757858753205 "source-unit" "angstrom" } "host-b" { "source-value" 4.015757858753205 "source-unit" "angstrom" } "host-c" { "source-value" 4.015757858753205 "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" [ "Pb" ] } } ]