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 Sim_LAMMPS_TersoffZBL_DevanathanDiazdelaRubiaWeber_1998_SiC__SM_578912636995_000 [5.432005241513252] CELL_SIZE_MIN: 3 CELL_SIZE_MAX: 5 Smallest System Size: 216 Largest System Size: 1000 [Calculation] Supercell Size: 3 Unrelaxed Cell: [[16.29601572 0. 0. ] [ 0. 16.29601572 0. ] [ 0. 0. 16.29601572]] Unrelaxed Cell Vector: [16.296015724539757, 0.0, 16.296015724539757, 0.0, 0.0, 16.296015724539757] Unrelaxed Cell Energy: -999.9925184402206 Energy of Unrelaxed Cell With Vacancy: -999.9925184402206 ========== Loop: 0 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:19:41 -991.259582* 0.8731 FIRE: 1 16:19:41 -991.298224* 0.8149 FIRE: 2 16:19:41 -991.361002* 0.7058 FIRE: 3 16:19:41 -991.428220* 0.5600 FIRE: 4 16:19:41 -991.486574* 0.3987 FIRE: 5 16:19:41 -991.533525* 0.2465 FIRE: 6 16:19:41 -991.570914* 0.2255 FIRE: 7 16:19:41 -991.597688* 0.1621 FIRE: 8 16:19:41 -991.613027* 0.1845 FIRE: 9 16:19:41 -991.618304* 0.1987 FIRE: 10 16:19:41 -991.620484* 0.1871 FIRE: 11 16:19:41 -991.624396* 0.1651 FIRE: 12 16:19:41 -991.629283* 0.1348 FIRE: 13 16:19:41 -991.634319* 0.0992 FIRE: 14 16:19:41 -991.638863* 0.0623 FIRE: 15 16:19:41 -991.642631* 0.0733 FIRE: 16 16:19:41 -991.645695* 0.0783 FIRE: 17 16:19:41 -991.648501* 0.0730 FIRE: 18 16:19:41 -991.650968* 0.0550 FIRE: 19 16:19:41 -991.652639* 0.0339 FIRE: 20 16:19:41 -991.652840* 0.0466 FIRE: 21 16:19:41 -991.652918* 0.0462 FIRE: 22 16:19:41 -991.653069* 0.0454 FIRE: 23 16:19:41 -991.653279* 0.0442 FIRE: 24 16:19:41 -991.653535* 0.0426 FIRE: 25 16:19:41 -991.653817* 0.0407 FIRE: 26 16:19:41 -991.654110* 0.0383 FIRE: 27 16:19:41 -991.654396* 0.0357 FIRE: 28 16:19:41 -991.654691* 0.0323 FIRE: 29 16:19:41 -991.654979* 0.0283 FIRE: 30 16:19:41 -991.655243* 0.0234 FIRE: 31 16:19:41 -991.655463* 0.0178 FIRE: 32 16:19:41 -991.655625* 0.0116 FIRE: 33 16:19:41 -991.655722* 0.0126 FIRE: 34 16:19:41 -991.655771* 0.0169 FIRE: 35 16:19:42 -991.655781* 0.0165 FIRE: 36 16:19:42 -991.655801* 0.0158 FIRE: 37 16:19:42 -991.655829* 0.0148 FIRE: 38 16:19:42 -991.655863* 0.0134 FIRE: 39 16:19:42 -991.655901* 0.0119 FIRE: 40 16:19:42 -991.655941* 0.0102 FIRE: 41 16:19:42 -991.655981* 0.0084 FIRE: 42 16:19:42 -991.656023* 0.0077 FIRE: 43 16:19:42 -991.656066* 0.0067 FIRE: 44 16:19:42 -991.656104* 0.0054 FIRE: 45 16:19:42 -991.656137* 0.0057 FIRE: 46 16:19:42 -991.656162* 0.0068 FIRE: 47 16:19:42 -991.656178* 0.0077 FIRE: 48 16:19:42 -991.656190* 0.0082 FIRE: 49 16:19:42 -991.656203* 0.0080 FIRE: 50 16:19:42 -991.656220* 0.0070 WARNING: Max number of steps exceeded. Structure may be unstable. Relaxation Completed. Steps: 50 Cell Size Relaxation... Optimization terminated successfully. Current function value: 3.701403 Iterations: 427 Function evaluations: 760 Current VFE: 3.701403261905625 Energy of Supercell: -999.9925184402206 Unrelaxed Cell Volume: 4327.572029756306 Current Relaxed Cell Volume: 4336.282918302424 Current Relaxation Volume: -8.71088854611844 Current Cell: [[ 1.63069422e+01 0.00000000e+00 0.00000000e+00] [ 3.07692915e-05 1.63069428e+01 0.00000000e+00] [-1.27036975e-05 1.57424447e-05 1.63069422e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:19:44 -991.661520* 0.0189 FIRE: 1 16:19:44 -991.661539* 0.0178 FIRE: 2 16:19:44 -991.661571* 0.0159 FIRE: 3 16:19:44 -991.661610* 0.0133 FIRE: 4 16:19:44 -991.661646* 0.0102 FIRE: 5 16:19:44 -991.661676* 0.0070 FIRE: 6 16:19:44 -991.661697* 0.0041 FIRE: 7 16:19:44 -991.661710* 0.0034 FIRE: 8 16:19:44 -991.661720* 0.0039 FIRE: 9 16:19:44 -991.661728* 0.0037 FIRE: 10 16:19:44 -991.661734* 0.0030 FIRE: 11 16:19:44 -991.661737* 0.0038 FIRE: 12 16:19:44 -991.661738* 0.0037 FIRE: 13 16:19:44 -991.661739* 0.0036 FIRE: 14 16:19:44 -991.661742* 0.0034 FIRE: 15 16:19:44 -991.661745* 0.0032 FIRE: 16 16:19:44 -991.661747* 0.0029 FIRE: 17 16:19:44 -991.661749* 0.0025 FIRE: 18 16:19:44 -991.661751* 0.0021 FIRE: 19 16:19:44 -991.661752* 0.0015 FIRE: 20 16:19:44 -991.661753* 0.0008 Relaxation Completed. Steps: 20 Cell Size Relaxation... Optimization terminated successfully. Current function value: 3.701167 Iterations: 176 Function evaluations: 378 Current VFE: 3.701166878076151 Energy of Supercell: -999.9925184402206 Unrelaxed Cell Volume: 4327.572029756306 Current Relaxed Cell Volume: 4336.522923774192 Current Relaxation Volume: -8.95089401788664 Current Cell: [[ 1.63072431e+01 0.00000000e+00 0.00000000e+00] [ 3.06727214e-05 1.63072433e+01 0.00000000e+00] [-1.29865288e-05 1.58457368e-05 1.63072433e+01]] ========== Loop: 2 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:19:45 -991.661757* 0.0004 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 3.701167 Iterations: 218 Function evaluations: 429 Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:19:46 -991.661757* 0.0004 FIRE: 1 16:19:46 -991.661757* 0.0004 FIRE: 2 16:19:46 -991.661757* 0.0004 FIRE: 3 16:19:46 -991.661757* 0.0003 FIRE: 4 16:19:46 -991.661757* 0.0002 FIRE: 5 16:19:46 -991.661757* 0.0003 FIRE: 6 16:19:46 -991.661757* 0.0003 FIRE: 7 16:19:46 -991.661757* 0.0002 FIRE: 8 16:19:46 -991.661757* 0.0002 FIRE: 9 16:19:46 -991.661757* 0.0002 FIRE: 10 16:19:46 -991.661757* 0.0002 FIRE: 11 16:19:46 -991.661757* 0.0002 FIRE: 12 16:19:46 -991.661757* 0.0002 FIRE: 13 16:19:46 -991.661757* 0.0002 FIRE: 14 16:19:46 -991.661757* 0.0002 FIRE: 15 16:19:46 -991.661757* 0.0002 FIRE: 16 16:19:46 -991.661757* 0.0002 FIRE: 17 16:19:46 -991.661757* 0.0002 FIRE: 18 16:19:46 -991.661757* 0.0002 FIRE: 19 16:19:46 -991.661757* 0.0001 FIRE: 20 16:19:46 -991.661757* 0.0001 Optimization terminated successfully. Current function value: 3.701166 Iterations: 289 Function evaluations: 590 --------------- Calculation Completed. Number Of Atoms in Supercell: 216 Vacancy Formation Energy (relaxed): 3.701165978346239 Vacancy Formation Energy (unrelaxed): 4.103341059085551 Unrelaxed Cell Volume: 4327.572029756306 Relaxed Cell Volume: 4336.522923774192 Relaxation Volume: -8.95089401788664 Relaxed Cell Vector: [16.307235322298887, 1.1178303226676424e-05, 16.307233975026, -1.983577777842123e-05, 1.0228449657685518e-05, 16.30723499226708] Unrelaxed Cell Vector: [16.296015724539757, 0.0, 16.296015724539757, 0.0, 0.0, 16.296015724539757] Relaxed Cell: [[ 1.63072353e+01 0.00000000e+00 0.00000000e+00] [ 1.11783032e-05 1.63072340e+01 0.00000000e+00] [-1.98357778e-05 1.02284497e-05 1.63072350e+01]] Unrelaxed Cell: [[16.29601572 0. 0. ] [ 0. 16.29601572 0. ] [ 0. 0. 16.29601572]] Supercell Size: 4 Unrelaxed Cell: [[21.72802097 0. 0. ] [ 0. 21.72802097 0. ] [ 0. 0. 21.72802097]] Unrelaxed Cell Vector: [21.72802096605301, 0.0, 21.72802096605301, 0.0, 0.0, 21.72802096605301] Unrelaxed Cell Energy: -2370.3526363025494 Energy of Unrelaxed Cell With Vacancy: -2370.3526363025494 ========== Loop: 0 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:19:47 -2361.619700* 0.8731 FIRE: 1 16:19:47 -2361.658342* 0.8149 FIRE: 2 16:19:47 -2361.721120* 0.7058 FIRE: 3 16:19:47 -2361.788338* 0.5600 FIRE: 4 16:19:47 -2361.846692* 0.3987 FIRE: 5 16:19:47 -2361.893643* 0.2465 FIRE: 6 16:19:47 -2361.931032* 0.2255 FIRE: 7 16:19:47 -2361.957808* 0.1621 FIRE: 8 16:19:47 -2361.973151* 0.1844 FIRE: 9 16:19:47 -2361.978434* 0.1986 FIRE: 10 16:19:47 -2361.980612* 0.1871 FIRE: 11 16:19:47 -2361.984518* 0.1651 FIRE: 12 16:19:47 -2361.989394* 0.1347 FIRE: 13 16:19:47 -2361.994413* 0.0991 FIRE: 14 16:19:47 -2361.998930* 0.0621 FIRE: 15 16:19:47 -2362.002664* 0.0732 FIRE: 16 16:19:47 -2362.005688* 0.0782 FIRE: 17 16:19:47 -2362.008455* 0.0728 FIRE: 18 16:19:47 -2362.010906* 0.0549 FIRE: 19 16:19:47 -2362.012631* 0.0337 FIRE: 20 16:19:47 -2362.013034* 0.0463 FIRE: 21 16:19:47 -2362.013125* 0.0459 FIRE: 22 16:19:47 -2362.013299* 0.0451 FIRE: 23 16:19:47 -2362.013545* 0.0439 FIRE: 24 16:19:47 -2362.013846* 0.0423 FIRE: 25 16:19:47 -2362.014181* 0.0404 FIRE: 26 16:19:47 -2362.014533* 0.0381 FIRE: 27 16:19:47 -2362.014883* 0.0354 FIRE: 28 16:19:47 -2362.015252* 0.0322 FIRE: 29 16:19:47 -2362.015621* 0.0282 FIRE: 30 16:19:47 -2362.015973* 0.0234 FIRE: 31 16:19:47 -2362.016286* 0.0179 FIRE: 32 16:19:47 -2362.016546* 0.0118 FIRE: 33 16:19:47 -2362.016754* 0.0137 FIRE: 34 16:19:47 -2362.016934* 0.0184 FIRE: 35 16:19:47 -2362.017137* 0.0199 FIRE: 36 16:19:47 -2362.017398* 0.0172 FIRE: 37 16:19:47 -2362.017701* 0.0104 FIRE: 38 16:19:47 -2362.017957* 0.0070 FIRE: 39 16:19:47 -2362.018062* 0.0049 FIRE: 40 16:19:47 -2362.018070* 0.0045 FIRE: 41 16:19:47 -2362.018084* 0.0036 FIRE: 42 16:19:47 -2362.018101* 0.0025 FIRE: 43 16:19:47 -2362.018117* 0.0020 FIRE: 44 16:19:47 -2362.018131* 0.0020 FIRE: 45 16:19:47 -2362.018141* 0.0019 FIRE: 46 16:19:47 -2362.018148* 0.0020 FIRE: 47 16:19:47 -2362.018152* 0.0020 FIRE: 48 16:19:47 -2362.018154* 0.0016 FIRE: 49 16:19:47 -2362.018154* 0.0016 FIRE: 50 16:19:47 -2362.018155* 0.0015 WARNING: Max number of steps exceeded. Structure may be unstable. Relaxation Completed. Steps: 50 Cell Size Relaxation... Optimization terminated successfully. Current function value: 3.702578 Iterations: 319 Function evaluations: 596 Current VFE: 3.702577510900028 Energy of Supercell: -2370.3526363025494 Unrelaxed Cell Volume: 10257.948514977914 Current Relaxed Cell Volume: 10266.776627313895 Current Relaxation Volume: -8.828112335981132 Current Cell: [[ 2.17342525e+01 0.00000000e+00 0.00000000e+00] [-1.13346420e-06 2.17342524e+01 0.00000000e+00] [-4.44067876e-06 -2.18925445e-06 2.17342520e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:19:49 -2362.020464* 0.0035 FIRE: 1 16:19:49 -2362.020466* 0.0033 FIRE: 2 16:19:49 -2362.020469* 0.0030 FIRE: 3 16:19:49 -2362.020473* 0.0026 FIRE: 4 16:19:49 -2362.020475* 0.0021 FIRE: 5 16:19:49 -2362.020477* 0.0015 FIRE: 6 16:19:49 -2362.020478* 0.0011 FIRE: 7 16:19:49 -2362.020479* 0.0008 Relaxation Completed. Steps: 7 Cell Size Relaxation... Optimization terminated successfully. Current function value: 3.702563 Iterations: 164 Function evaluations: 362 Current VFE: 3.7025626224758525 Energy of Supercell: -2370.3526363025494 Unrelaxed Cell Volume: 10257.948514977914 Current Relaxed Cell Volume: 10266.808930858417 Current Relaxation Volume: -8.860415880502842 Current Cell: [[ 2.17342752e+01 0.00000000e+00 0.00000000e+00] [-1.15228572e-06 2.17342750e+01 0.00000000e+00] [-4.37393688e-06 -2.23664412e-06 2.17342751e+01]] ========== Loop: 2 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:19:50 -2362.020479* 0.0008 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 3.702563 Iterations: 119 Function evaluations: 291 Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:19:52 -2362.020479* 0.0008 FIRE: 1 16:19:52 -2362.020479* 0.0007 FIRE: 2 16:19:52 -2362.020480* 0.0006 FIRE: 3 16:19:52 -2362.020480* 0.0006 FIRE: 4 16:19:52 -2362.020481* 0.0006 FIRE: 5 16:19:52 -2362.020481* 0.0005 FIRE: 6 16:19:52 -2362.020481* 0.0005 FIRE: 7 16:19:52 -2362.020481* 0.0005 FIRE: 8 16:19:52 -2362.020481* 0.0005 FIRE: 9 16:19:52 -2362.020481* 0.0005 FIRE: 10 16:19:52 -2362.020481* 0.0005 FIRE: 11 16:19:52 -2362.020481* 0.0004 FIRE: 12 16:19:52 -2362.020481* 0.0004 FIRE: 13 16:19:52 -2362.020481* 0.0004 FIRE: 14 16:19:52 -2362.020481* 0.0003 FIRE: 15 16:19:52 -2362.020481* 0.0003 FIRE: 16 16:19:52 -2362.020481* 0.0003 FIRE: 17 16:19:52 -2362.020481* 0.0002 FIRE: 18 16:19:52 -2362.020482* 0.0002 FIRE: 19 16:19:52 -2362.020482* 0.0002 FIRE: 20 16:19:52 -2362.020482* 0.0001 Optimization terminated successfully. Current function value: 3.702560 Iterations: 177 Function evaluations: 438 --------------- Calculation Completed. Number Of Atoms in Supercell: 512 Vacancy Formation Energy (relaxed): 3.7025595951222385 Vacancy Formation Energy (unrelaxed): 4.10334105908305 Unrelaxed Cell Volume: 10257.948514977914 Relaxed Cell Volume: 10266.808930858417 Relaxation Volume: -8.860415880502842 Relaxed Cell Vector: [21.73428290953624, -1.1802914810465814e-06, 21.734282551134022, -4.392616111233137e-06, -2.2369107437039363e-06, 21.734282982539792] Unrelaxed Cell Vector: [21.72802096605301, 0.0, 21.72802096605301, 0.0, 0.0, 21.72802096605301] Relaxed Cell: [[ 2.17342829e+01 0.00000000e+00 0.00000000e+00] [-1.18029148e-06 2.17342826e+01 0.00000000e+00] [-4.39261611e-06 -2.23691074e-06 2.17342830e+01]] Unrelaxed Cell: [[21.72802097 0. 0. ] [ 0. 21.72802097 0. ] [ 0. 0. 21.72802097]] Supercell Size: 5 Unrelaxed Cell: [[27.16002621 0. 0. ] [ 0. 27.16002621 0. ] [ 0. 0. 27.16002621]] Unrelaxed Cell Vector: [27.16002620756626, 0.0, 27.16002620756626, 0.0, 0.0, 27.16002620756626] Unrelaxed Cell Energy: -4629.594992778961 Energy of Unrelaxed Cell With Vacancy: -4629.594992778961 ========== Loop: 0 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:19:53 -4620.862057* 0.8731 FIRE: 1 16:19:53 -4620.900699* 0.8149 FIRE: 2 16:19:53 -4620.963476* 0.7058 FIRE: 3 16:19:53 -4621.030694* 0.5600 FIRE: 4 16:19:53 -4621.089048* 0.3987 FIRE: 5 16:19:53 -4621.135999* 0.2465 FIRE: 6 16:19:53 -4621.173389* 0.2255 FIRE: 7 16:19:53 -4621.200164* 0.1621 FIRE: 8 16:19:53 -4621.215508* 0.1844 FIRE: 9 16:19:53 -4621.220791* 0.1986 FIRE: 10 16:19:53 -4621.222969* 0.1871 FIRE: 11 16:19:53 -4621.226875* 0.1651 FIRE: 12 16:19:53 -4621.231751* 0.1347 FIRE: 13 16:19:53 -4621.236770* 0.0991 FIRE: 14 16:19:53 -4621.241287* 0.0621 FIRE: 15 16:19:53 -4621.245022* 0.0732 FIRE: 16 16:19:53 -4621.248047* 0.0782 FIRE: 17 16:19:53 -4621.250815* 0.0728 FIRE: 18 16:19:53 -4621.253268* 0.0549 FIRE: 19 16:19:53 -4621.254993* 0.0337 FIRE: 20 16:19:53 -4621.255392* 0.0463 FIRE: 21 16:19:53 -4621.255482* 0.0459 FIRE: 22 16:19:53 -4621.255655* 0.0451 FIRE: 23 16:19:53 -4621.255900* 0.0439 FIRE: 24 16:19:53 -4621.256198* 0.0423 FIRE: 25 16:19:53 -4621.256532* 0.0404 FIRE: 26 16:19:53 -4621.256881* 0.0381 FIRE: 27 16:19:53 -4621.257229* 0.0354 FIRE: 28 16:19:53 -4621.257596* 0.0322 FIRE: 29 16:19:53 -4621.257965* 0.0282 FIRE: 30 16:19:53 -4621.258318* 0.0234 FIRE: 31 16:19:53 -4621.258639* 0.0179 FIRE: 32 16:19:53 -4621.258915* 0.0118 FIRE: 33 16:19:53 -4621.259151* 0.0136 FIRE: 34 16:19:53 -4621.259374* 0.0183 FIRE: 35 16:19:53 -4621.259630* 0.0199 FIRE: 36 16:19:53 -4621.259956* 0.0173 FIRE: 37 16:19:53 -4621.260326* 0.0105 FIRE: 38 16:19:53 -4621.260652* 0.0070 FIRE: 39 16:19:53 -4621.260836* 0.0053 FIRE: 40 16:19:53 -4621.260832* 0.0090 FIRE: 41 16:19:54 -4621.260848* 0.0084 FIRE: 42 16:19:54 -4621.260877* 0.0073 FIRE: 43 16:19:54 -4621.260913* 0.0057 FIRE: 44 16:19:54 -4621.260949* 0.0040 FIRE: 45 16:19:54 -4621.260981* 0.0032 FIRE: 46 16:19:54 -4621.261004* 0.0030 FIRE: 47 16:19:54 -4621.261019* 0.0039 FIRE: 48 16:19:54 -4621.261028* 0.0046 FIRE: 49 16:19:54 -4621.261035* 0.0047 FIRE: 50 16:19:54 -4621.261044* 0.0043 WARNING: Max number of steps exceeded. Structure may be unstable. Relaxation Completed. Steps: 50 Cell Size Relaxation... Optimization terminated successfully. Current function value: 3.703164 Iterations: 424 Function evaluations: 753 Current VFE: 3.703163660473365 Energy of Supercell: -4629.594992778961 Unrelaxed Cell Volume: 20035.05569331622 Current Relaxed Cell Volume: 20043.904282952564 Current Relaxation Volume: -8.848589636345423 Current Cell: [[ 2.71640248e+01 0.00000000e+00 0.00000000e+00] [-8.18513389e-07 2.71640238e+01 0.00000000e+00] [ 2.83369939e-07 -3.60303419e-07 2.71640236e+01]] ========== Loop: 1 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:19:59 -4621.262234* 0.0026 FIRE: 1 16:19:59 -4621.262239* 0.0023 FIRE: 2 16:19:59 -4621.262248* 0.0017 FIRE: 3 16:19:59 -4621.262257* 0.0012 FIRE: 4 16:19:59 -4621.262264* 0.0006 Relaxation Completed. Steps: 4 Cell Size Relaxation... Optimization terminated successfully. Current function value: 3.703133 Iterations: 116 Function evaluations: 291 Current VFE: 3.703133360590982 Energy of Supercell: -4629.594992778961 Unrelaxed Cell Volume: 20035.05569331622 Current Relaxed Cell Volume: 20043.899536642966 Current Relaxation Volume: -8.843843326747447 Current Cell: [[ 2.71640217e+01 0.00000000e+00 0.00000000e+00] [-8.41016014e-07 2.71640223e+01 0.00000000e+00] [ 2.89779746e-07 -3.55592974e-07 2.71640218e+01]] ========== Loop: 2 Position Relaxation... Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:20:03 -4621.262264* 0.0006 Relaxation Completed. Steps: 0 Cell Size Relaxation... Optimization terminated successfully. Current function value: 3.703133 Iterations: 111 Function evaluations: 281 Step Time Energy fmax *Force-consistent energies used in optimization. FIRE: 0 16:20:05 -4621.262264* 0.0006 FIRE: 1 16:20:05 -4621.262265* 0.0006 FIRE: 2 16:20:05 -4621.262266* 0.0005 FIRE: 3 16:20:05 -4621.262268* 0.0005 FIRE: 4 16:20:05 -4621.262270* 0.0004 FIRE: 5 16:20:05 -4621.262272* 0.0004 FIRE: 6 16:20:05 -4621.262273* 0.0003 FIRE: 7 16:20:05 -4621.262275* 0.0003 FIRE: 8 16:20:05 -4621.262277* 0.0003 FIRE: 9 16:20:05 -4621.262279* 0.0003 FIRE: 10 16:20:05 -4621.262281* 0.0002 FIRE: 11 16:20:05 -4621.262282* 0.0002 FIRE: 12 16:20:05 -4621.262282* 0.0003 FIRE: 13 16:20:05 -4621.262282* 0.0003 FIRE: 14 16:20:05 -4621.262282* 0.0003 FIRE: 15 16:20:05 -4621.262283* 0.0002 FIRE: 16 16:20:05 -4621.262283* 0.0002 FIRE: 17 16:20:05 -4621.262283* 0.0002 FIRE: 18 16:20:05 -4621.262283* 0.0002 FIRE: 19 16:20:05 -4621.262283* 0.0001 FIRE: 20 16:20:05 -4621.262283* 0.0001 Optimization terminated successfully. Current function value: 3.703115 Iterations: 184 Function evaluations: 433 --------------- Calculation Completed. Number Of Atoms in Supercell: 1000 Vacancy Formation Energy (relaxed): 3.703114660170286 Vacancy Formation Energy (unrelaxed): 4.103341059093509 Unrelaxed Cell Volume: 20035.05569331622 Relaxed Cell Volume: 20043.899536642966 Relaxation Volume: -8.843843326747447 Relaxed Cell Vector: [27.164021730932852, -8.412008472860487e-07, 27.164022177287553, 2.8941275128897484e-07, -3.7353466679935497e-07, 27.164021820805157] Unrelaxed Cell Vector: [27.16002620756626, 0.0, 27.16002620756626, 0.0, 0.0, 27.16002620756626] Relaxed Cell: [[ 2.71640217e+01 0.00000000e+00 0.00000000e+00] [-8.41200847e-07 2.71640222e+01 0.00000000e+00] [ 2.89412751e-07 -3.73534667e-07 2.71640218e+01]] Unrelaxed Cell: [[27.16002621 0. 0. ] [ 0. 27.16002621 0. ] [ 0. 0. 27.16002621]] [Calculation Results Summary] Sizes: [3, 4, 5] Unrelaxed Formation Energy By Size: [4.103341059085551, 4.10334105908305, 4.103341059093509] Formation Energy By Size: [3.701165978346239, 3.7025595951222385, 3.703114660170286] Relaxation Volume By Size: [-8.95089401788664, -8.860415880502842, -8.843843326747447] [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: [4.10334106 4.10334106] Fitting Results: (array([4.10334106e+00, 1.16796076e-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: [3.70116598 3.7025596 ] Fitting Results: (array([ 3.70357656, -0.06508567]), 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: [-8.95089402 -8.86041588] Fitting Results: (array([-8.79439129, -4.22557355]), 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: [4.10334106 4.10334106] Fitting Results: (array([ 4.10334106e+00, -1.37171516e-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: [3.7025596 3.70311466] Fitting Results: (array([ 3.70369702, -0.07279542]), 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: [-8.86041588 -8.84384333] Fitting Results: (array([-8.82645573, -2.17344967]), 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: [4.10334106 4.10334106 4.10334106] Fitting Results: (array([ 4.10334106e+00, -1.89228165e-10]), array([4.34326812e-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: [3.70116598 3.7025596 3.70311466] Fitting Results: (array([ 3.7036282 , -0.06667068]), array([1.16523043e-09]), 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: [-8.95089402 -8.86041588 -8.84384333] Fitting Results: (array([-8.80813795, -3.80368669]), array([8.25541301e-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: [4.10334106 4.10334106 4.10334106] Fitting Results: (array([ 4.10334106e+00, -6.28202068e-09, 1.62346749e-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: [3.70116598 3.7025596 3.70311466] Fitting Results: (array([ 3.70376595, -0.09822901, 0.08408939]), 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: [-8.95089402 -8.86041588 -8.84384333] Fitting Results: (array([ -8.84480188, 4.59627825, -22.38229829]), 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: [4.10334106 4.10334106 4.10334106] Fitting Results: (array([ 4.10334106e+00, -3.49095917e-09, 2.46120886e-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: [3.70116598 3.7025596 3.70311466] Fitting Results: (array([ 3.70374405, -0.08377238, 0.12748118]), 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: [-8.95089402 -8.86041588 -8.84384333] Fitting Results: (array([ -8.83897163, 0.74831866, -33.93200731]), 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: [4.10334106 4.10334106 4.10334106] Fitting Results: (array([ 4.10334106e+00, -2.58274675e-09, 5.12613681e-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: [3.70116598 3.7025596 3.70311466] Fitting Results: (array([ 3.70373021, -0.07906819, 0.26551422]), 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: [-8.95089402 -8.86041588 -8.84384333] Fitting Results: (array([ -8.83528981, -0.5038087 , -70.67263357]), 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: [[4.103341059081223, 4.103341059104481], [4.103341059091198], [4.10334105911779], [4.10334105911356], [4.103341059110891]] Formation Energy Fits By Size: [[3.703576558715533, 3.7036970234993816], [3.7036282043359305], [3.703765949226226], [3.7037440452442647], [3.7037302127765668]] Relaxation Volume Fits By Size: [[-8.794391293763306, -8.82645572936473], [-8.808137947484752], [-8.844801875507818], [-8.838971633701469], [-8.835289808561043]] [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" 4.103341059104481 "source-unit" "eV" "source-std-uncert-value" 1.87004206964095e-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.432005241513252 "source-unit" "angstrom" } "host-b" { "source-value" 5.432005241513252 "source-unit" "angstrom" } "host-c" { "source-value" 5.432005241513252 "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.629594992778737 "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.432005241513252 "source-unit" "angstrom" } "reservoir-b" { "source-value" 5.432005241513252 "source-unit" "angstrom" } "reservoir-c" { "source-value" 5.432005241513252 "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" 3.7036970234993816 "source-unit" "eV" "source-std-uncert-value" 7.141751574530256e-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.432005241513252 "source-unit" "angstrom" } "host-b" { "source-value" 5.432005241513252 "source-unit" "angstrom" } "host-c" { "source-value" 5.432005241513252 "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.629594992778737 "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.432005241513252 "source-unit" "angstrom" } "reservoir-b" { "source-value" 5.432005241513252 "source-unit" "angstrom" } "reservoir-c" { "source-value" 5.432005241513252 "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" -8.82645572936473 "source-unit" "angstrom^3" "source-std-uncert-value" 0.018346212479868854 } "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.432005241513252 "source-unit" "angstrom" } "host-b" { "source-value" 5.432005241513252 "source-unit" "angstrom" } "host-c" { "source-value" 5.432005241513252 "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" ] } } ]