!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!! !!!!! !!!!! VERIFICATION CHECK: vc-forces-numerical-derivative !!!!! !!!!! !!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Description: Comparison of the analytical forces obtained from the model with forces computed by numerical differentiation using Richardson extrapolation for a randomly distorted non-periodic structure based off the equilibrium face-centered cubic (fcc) structure. Separate configurations are tested for each species supported by the model, as well as one containing a random distribution of all species. Configurations used for testing are provided as auxiliary files. Author: Ellad Tadmor ----------------------------------------------------------------------------------------------------- Results for KIM Model : Sim_LAMMPS_ReaxFF_ChenowethVanDuinPersson_2008_CHOV__SM_429148913211_001 Supported species : C H O V random seed = 13 perturbation amplitude = 0.250 number unit cells per side = 2 ----------------------------------------------------------------------------------------------------- Comparison of analytical forces obtained from the model, the force computed as a numerical derivative of the energy, the difference between them, and the uncertainty in the numerical estimate of the force. The computed equilibrium lattice constant for the crystal (a0) is given in the heading. MONOATOMIC STRUCTURE -- Species = C | a0 = 2.893721580505371 | (Configuration in file "config-C.xyz") ----------------------------------------------------------------------------------------------------- Part Spec Dir Force_model Force_numer |Force diff| uncertainty ----------------------------------------------------------------------------------------------------- 1 C 1 2.343056362658277e+00 2.343056340923895e+00 2.17344e-08 6.11945e-07 2 2.286267798857737e+00 2.286267813896410e+00 1.50387e-08 4.17381e-07 3 4.177354173275275e+00 4.177354381686920e+00 2.08412e-07 1.41267e-06 ----------------------------------------------------------------------------------------------------- 2 C 1 3.424308003071628e+00 3.424307983314120e+00 1.97575e-08 2.24797e-07 2 -4.139595757238918e+00 -4.139595800198744e+00 4.29598e-08 3.40858e-04 * 3 -7.219626138063542e+00 -7.219626098728815e+00 3.93347e-08 1.95140e-06 ----------------------------------------------------------------------------------------------------- 3 C 1 -1.415220043285152e+00 -1.415219607770922e+00 4.35514e-07 7.77602e-06 2 -1.403922612631774e+00 -1.403922719612918e+00 1.06981e-07 8.63924e-07 3 -1.388728058809186e+00 -1.388727629390297e+00 4.29419e-07 3.80689e-06 ----------------------------------------------------------------------------------------------------- 4 C 1 -2.453647595766995e+00 -2.453647673804066e+00 7.80371e-08 3.48200e-04 * 2 -2.675607557378011e+00 -2.675607572094681e+00 1.47167e-08 3.40181e-04 * 3 -5.086319600529071e+00 -5.086319483647073e+00 1.16882e-07 6.70704e-07 ----------------------------------------------------------------------------------------------------- 5 C 1 -1.401956018365375e+00 -1.401956253112506e+00 2.34747e-07 5.31167e-07 2 4.040372923214254e+00 4.040372792753584e+00 1.30461e-07 8.28039e-07 3 3.942728961198716e+00 3.942728938180918e+00 2.30178e-08 5.29426e-08 ----------------------------------------------------------------------------------------------------- 6 C 1 -2.982542765660998e+00 -2.982543204806993e+00 4.39146e-07 1.52426e-06 2 -1.083691928446810e+00 -1.083692081478499e+00 1.53032e-07 2.26174e-04 * 3 1.130075716412159e+00 1.130075971565026e+00 2.55153e-07 3.34923e-07 ----------------------------------------------------------------------------------------------------- 7 C 1 -4.093824667843245e+00 -4.093824684915258e+00 1.70720e-08 3.32028e-07 2 1.015628433081185e+01 1.015628424045346e+01 9.03584e-08 1.86208e-07 3 5.334086073742755e+00 5.334086293566635e+00 2.19824e-07 3.49722e-07 ----------------------------------------------------------------------------------------------------- 8 C 1 1.146023828595634e+01 1.146023793366760e+01 3.52289e-07 8.37673e-07 2 6.098853077854146e+00 6.098853128739350e+00 5.08852e-08 2.39060e-06 3 1.146275209379772e+00 1.146275337895564e+00 1.28516e-07 7.83060e-07 ----------------------------------------------------------------------------------------------------- 9 C 1 -6.531327898700781e+00 -6.531327703554995e+00 1.95146e-07 2.02394e-06 2 2.916706394471148e+00 2.916706344402521e+00 5.00686e-08 5.30687e-07 3 2.601442708777334e-01 2.601442901960860e-01 1.93184e-08 1.29299e-06 ----------------------------------------------------------------------------------------------------- 10 C 1 -8.549910640720025e-01 -8.549911822494430e-01 1.18177e-07 4.00406e-06 2 5.250389043237627e+00 5.250388989419569e+00 5.38181e-08 1.00635e-03 * 3 -1.557916596359359e+00 -1.557916511219614e+00 8.51397e-08 1.49903e-06 ----------------------------------------------------------------------------------------------------- 11 C 1 2.350834719423057e+01 2.350834810928716e+01 9.15057e-07 8.84096e-06 2 -9.808735991257856e-01 -9.808734836884428e-01 1.15437e-07 1.27937e-05 * 3 3.565672316270461e+01 3.565672339109943e+01 2.28395e-07 7.44282e-07 ----------------------------------------------------------------------------------------------------- 12 C 1 -8.294770763051660e-01 -8.294771374848284e-01 6.11797e-08 1.14690e-06 2 -2.128836890344073e+00 -2.128837046319176e+00 1.55975e-07 2.59277e-06 3 -6.096823541527592e+00 -6.096823599203689e+00 5.76761e-08 2.58502e-07 ----------------------------------------------------------------------------------------------------- 13 C 1 8.911081094801888e+00 8.911080639784148e+00 4.55018e-07 1.02675e-03 * 2 -8.187523146538240e+00 -8.187523065973975e+00 8.05643e-08 3.82639e-07 3 2.113852137151988e+00 2.113852013775960e+00 1.23376e-07 2.00747e-06 ----------------------------------------------------------------------------------------------------- 14 C 1 1.940999182284102e+00 1.940999088213708e+00 9.40704e-08 2.85281e-06 2 -4.088087054398648e+00 -4.088087054952855e+00 5.54206e-10 1.48047e-06 3 1.179053593205071e+00 1.179053571213057e+00 2.19920e-08 2.56755e-06 ----------------------------------------------------------------------------------------------------- 15 C 1 -8.230372570174252e+00 -8.230372330606254e+00 2.39568e-07 3.96442e-06 2 9.789304149133452e+00 9.789303919601890e+00 2.29532e-07 4.35113e-06 3 -2.742445726297773e+00 -2.742445785301507e+00 5.90037e-08 2.59370e-07 ----------------------------------------------------------------------------------------------------- 16 C 1 1.048852264266950e+01 1.048852211341241e+01 5.29257e-07 8.45275e-07 2 8.622015367939253e+00 8.622014812830448e+00 5.55109e-07 1.74299e-05 * 3 7.739888755039378e+00 7.739889186283918e+00 4.31245e-07 1.64816e-05 * ----------------------------------------------------------------------------------------------------- 17 C 1 3.161538769755230e+00 3.161538828030785e+00 5.82756e-08 1.04102e-06 2 1.824007842160408e+00 1.824007786797891e+00 5.53625e-08 7.83917e-07 3 1.568035913375297e+00 1.568035936322678e+00 2.29474e-08 2.60949e-06 ----------------------------------------------------------------------------------------------------- 18 C 1 -6.878085313115361e+00 -6.878085024525411e+00 2.88590e-07 4.00115e-05 * 2 -7.851120369189442e+00 -7.851120295300741e+00 7.38887e-08 3.76866e-06 3 8.331863471902537e+00 8.331863423790731e+00 4.81118e-08 2.28997e-06 ----------------------------------------------------------------------------------------------------- 19 C 1 -2.428296840938830e-01 -2.428297244961485e-01 4.04023e-08 5.81751e-07 2 4.399702891865987e+00 4.399702882610052e+00 9.25593e-09 1.66229e-07 3 -1.609144329869022e-01 -1.609144985580249e-01 6.55711e-08 4.40499e-07 ----------------------------------------------------------------------------------------------------- 20 C 1 3.875188839812279e+00 3.875188808668564e+00 3.11437e-08 4.35329e-07 2 -5.961293513387139e-02 -5.961300179572986e-02 6.66619e-08 2.68889e-06 3 -4.239476934577707e+00 -4.239476922730572e+00 1.18471e-08 6.13037e-07 ----------------------------------------------------------------------------------------------------- 21 C 1 4.586101019964175e+00 4.586100659222978e+00 3.60741e-07 6.96944e-07 2 -2.913727396771698e+00 -2.913727289423905e+00 1.07348e-07 1.64581e-06 3 4.038218058503071e-01 4.038218719812547e-01 6.61309e-08 2.46325e-06 ----------------------------------------------------------------------------------------------------- 22 C 1 7.384444507929069e-01 7.384444368215031e-01 1.39714e-08 3.29818e-07 2 -9.084587997005928e+00 -9.084588242264386e+00 2.45258e-07 3.02855e-06 3 -3.397562579401858e+00 -3.397562697356292e+00 1.17954e-07 1.08776e-07 ----------------------------------------------------------------------------------------------------- 23 C 1 -4.368634711463904e+00 -4.368634715610754e+00 4.14685e-09 6.86065e-08 2 1.563311498900333e+00 1.563311477794863e+00 2.11055e-08 1.24173e-07 3 -4.021875800507699e+00 -4.021875784745811e+00 1.57619e-08 1.25915e-07 ----------------------------------------------------------------------------------------------------- 24 C 1 -1.582688370938756e+00 -1.582688340994033e+00 2.99447e-08 1.63534e-06 2 -3.480503809610526e+00 -3.480503756884910e+00 5.27256e-08 5.48545e-07 3 -4.153790677232213e-01 -4.153790174591768e-01 5.02640e-08 7.40140e-07 ----------------------------------------------------------------------------------------------------- 25 C 1 2.493452333910392e+00 2.493452502495883e+00 1.68585e-07 1.29698e-06 2 1.859296184522082e+00 1.859296296452488e+00 1.11930e-07 2.81075e-06 3 2.470220723979357e+00 2.470220748141731e+00 2.41624e-08 2.69783e-06 ----------------------------------------------------------------------------------------------------- 26 C 1 3.835039970126414e+00 3.835040073754163e+00 1.03628e-07 6.18146e-07 2 -3.716029830777518e+00 -3.716030074008533e+00 2.43231e-07 1.01673e-03 * 3 7.971450351586915e-01 7.971449245410885e-01 1.10618e-07 1.37158e-06 ----------------------------------------------------------------------------------------------------- 27 C 1 -9.347076951822853e+00 -9.347076930433808e+00 2.13890e-08 8.99123e-08 2 9.687092994893661e+00 9.687092932584248e+00 6.23094e-08 2.98348e-06 3 4.836781435858258e-01 4.836781983521014e-01 5.47663e-08 1.64274e-07 ----------------------------------------------------------------------------------------------------- 28 C 1 8.914422219073394e+00 8.914422313457376e+00 9.43840e-08 7.45011e-07 2 -1.026505906724694e+01 -1.026505907069153e+01 3.44459e-09 1.44001e-06 3 1.277855967777989e+00 1.277855949769933e+00 1.80081e-08 1.58414e-07 ----------------------------------------------------------------------------------------------------- 29 C 1 -1.978625010281958e+01 -1.978625029547417e+01 1.92655e-07 5.29775e-06 2 1.841846684624594e+00 1.841846629370799e+00 5.52538e-08 3.56120e-06 3 -2.897409146079908e+01 -2.897409247177337e+01 1.01097e-06 9.37003e-06 ----------------------------------------------------------------------------------------------------- 30 C 1 -1.130220682490965e+01 -1.130220597690423e+01 8.48005e-07 3.55413e-06 2 -6.655751767148701e+00 -6.655751523465921e+00 2.43683e-07 1.88289e-06 3 -8.767191633332619e+00 -8.767191625978977e+00 7.35364e-09 1.65128e-05 * ----------------------------------------------------------------------------------------------------- 31 C 1 -6.229868289643939e+00 -6.229868515695773e+00 2.26052e-07 4.37893e-07 2 3.110573282525807e+00 3.110572828280836e+00 4.54245e-07 2.98019e-04 * 3 -3.036832063932976e+00 -3.036832071191335e+00 7.25836e-09 9.45349e-07 ----------------------------------------------------------------------------------------------------- 32 C 1 -1.149740420125205e+00 -1.149740375290846e+00 4.48344e-08 4.87993e-07 2 -4.731492746025453e+00 -4.731492770012283e+00 2.39868e-08 6.48712e-07 3 -9.076194797688624e-01 -9.076196315984179e-01 1.51830e-07 1.71305e-07 ----------------------------------------------------------------------------------------------------- * Starred lines are suspected outliers and are not included when determining the error. A calculation is considered an outlier if it has an uncertainty that lies at an abnormal distance from the other uncertainties in this set of calculations. Outliers are determined using the box plot construction with fences. An outlier could indicate a problem with the the numerical differentiation or problems with the potential energy, such as discontinuities. Maximum error obtained for particle = 20, direction = 2: |F_model - F_numer| error = ----------------------- = 1.11824e-06 max{|F_numer|, eps} Comparison of analytical forces obtained from the model, the force computed as a numerical derivative of the energy, the difference between them, and the uncertainty in the numerical estimate of the force. The computed equilibrium lattice constant for the crystal (a0) is given in the heading. MONOATOMIC STRUCTURE -- Species = H | a0 = 2.5463924407958984 | (Configuration in file "config-H.xyz") ----------------------------------------------------------------------------------------------------- Part Spec Dir Force_model Force_numer |Force diff| uncertainty ----------------------------------------------------------------------------------------------------- 1 H 1 -1.034253901745309e+00 -1.034253945090348e+00 4.33450e-08 1.73918e-07 2 -8.450471017690646e-01 -8.450470936736192e-01 8.09545e-09 1.34502e-07 3 -1.526888175492049e+00 -1.526888156366024e+00 1.91260e-08 2.24157e-07 ----------------------------------------------------------------------------------------------------- 2 H 1 1.770155238894070e+00 1.770156655531625e+00 1.41664e-06 9.50138e-06 2 1.116376206112292e-01 1.116370468108438e-01 5.73800e-07 2.80365e-05 * 3 -4.514185305471729e+00 -4.514189953617431e+00 4.64815e-06 7.83176e-06 ----------------------------------------------------------------------------------------------------- 3 H 1 -1.232327618701551e+00 -1.232327793243021e+00 1.74541e-07 3.45018e-06 2 -1.126050007240144e+00 -1.126049594938876e+00 4.12301e-07 5.29933e-06 3 1.939141250191755e+00 1.939141114293556e+00 1.35898e-07 1.76782e-06 ----------------------------------------------------------------------------------------------------- 4 H 1 -3.547686161926613e+00 -3.547685991136189e+00 1.70790e-07 2.43415e-06 2 -2.679702559395873e+00 -2.679702398451766e+00 1.60944e-07 1.87762e-06 3 -2.336909486187068e+00 -2.336909657608339e+00 1.71421e-07 1.45241e-06 ----------------------------------------------------------------------------------------------------- 5 H 1 3.143078968480343e-01 3.143079563571612e-01 5.95091e-08 3.27741e-07 2 5.671524535049494e-01 5.671524864219690e-01 3.29170e-08 9.35516e-07 3 -6.861646269671706e-01 -6.861646390975029e-01 1.21303e-08 1.53012e-06 ----------------------------------------------------------------------------------------------------- 6 H 1 3.489319586884462e-01 3.489311415937470e-01 8.17095e-07 5.83340e-06 2 -1.473545577040822e+00 -1.473545584809332e+00 7.76851e-09 2.82539e-05 * 3 5.333382672564390e-01 5.333358301567996e-01 2.43710e-06 7.01012e-07 ----------------------------------------------------------------------------------------------------- 7 H 1 -1.050716077382164e-01 -1.050715860466605e-01 2.16916e-08 2.13612e-07 2 3.791158331513158e-01 3.791160100503443e-01 1.76899e-07 7.89814e-07 3 -3.303899680655377e-01 -3.303900117637403e-01 4.36982e-08 2.34160e-07 ----------------------------------------------------------------------------------------------------- 8 H 1 3.120303440212434e+00 3.120326516403775e+00 2.30762e-05 4.00899e-04 * 2 -3.561577269816097e+00 -3.561575296877534e+00 1.97294e-06 1.06643e-05 3 5.821398296161483e-01 5.821654802403439e-01 2.56506e-05 1.44851e-04 * ----------------------------------------------------------------------------------------------------- 9 H 1 1.641148275053436e-01 1.641148369721486e-01 9.46681e-09 2.00423e-06 2 8.972960658253936e-01 8.972958540441294e-01 2.11781e-07 1.13873e-05 3 -3.056878567139375e-01 -3.056876558755547e-01 2.00838e-07 2.54185e-06 ----------------------------------------------------------------------------------------------------- 10 H 1 -3.252554980812577e+00 -3.252554956299023e+00 2.45136e-08 5.05076e-07 2 2.341220736543555e+00 2.341220700090457e+00 3.64531e-08 2.14285e-06 3 -1.333032226276360e+00 -1.333032076217048e+00 1.50059e-07 4.16322e-06 ----------------------------------------------------------------------------------------------------- 11 H 1 -2.195001384271058e+00 -2.195002439068425e+00 1.05480e-06 6.06775e-06 2 2.600374901366738e+00 2.600375641169416e+00 7.39803e-07 8.93336e-06 3 -1.819012527246159e-01 -1.819013601399066e-01 1.07415e-07 2.45102e-05 * ----------------------------------------------------------------------------------------------------- 12 H 1 1.505114612895090e+00 1.505114690512902e+00 7.76178e-08 1.57096e-06 2 -1.909928760625413e-01 -1.909926371388605e-01 2.38924e-07 7.34506e-07 3 -1.311797919935884e+00 -1.311797950088553e+00 3.01527e-08 1.43035e-06 ----------------------------------------------------------------------------------------------------- 13 H 1 4.368079954770432e-01 4.368078585987308e-01 1.36878e-07 3.55555e-05 * 2 3.328883054457170e+00 3.328881807164072e+00 1.24729e-06 6.17624e-06 3 -4.507064614814560e+00 -4.507064736333561e+00 1.21519e-07 1.38762e-05 ----------------------------------------------------------------------------------------------------- 14 H 1 5.849353688140488e-01 5.849352909506821e-01 7.78634e-08 3.56085e-07 2 9.674849053685403e-01 9.674848416039465e-01 6.37646e-08 8.46660e-07 3 -2.944571452495110e-03 -2.944512184660291e-03 5.92678e-08 2.00087e-06 ----------------------------------------------------------------------------------------------------- 15 H 1 3.798197397699600e+00 3.798201494751415e+00 4.09705e-06 4.04495e-06 2 -3.025116980343934e-01 -3.025105293441172e-01 1.16869e-06 1.15582e-04 * 3 1.785053429538393e+00 1.785053026121232e+00 4.03417e-07 1.61725e-05 ----------------------------------------------------------------------------------------------------- 16 H 1 -2.764310511571089e-01 -2.764311847865502e-01 1.33629e-07 2.82429e-06 2 5.237696753366040e-01 5.237692148202567e-01 4.60516e-07 2.72751e-06 3 -1.955586162105423e+00 -1.955586280614755e+00 1.18509e-07 2.77720e-05 * ----------------------------------------------------------------------------------------------------- 17 H 1 -4.566018211641752e-01 -4.566018311520448e-01 9.98787e-09 1.30448e-06 2 -2.250342959162033e+00 -2.250343013167352e+00 5.40053e-08 2.06335e-07 3 3.734949940164467e+00 3.734949874316263e+00 6.58482e-08 5.31666e-07 ----------------------------------------------------------------------------------------------------- 18 H 1 -4.600910420015627e+00 -4.600910559691645e+00 1.39676e-07 4.36547e-06 2 1.940929053706797e+00 1.940929343214842e+00 2.89508e-07 5.56010e-07 3 7.315529798182834e+00 7.315530048302237e+00 2.50119e-07 2.07184e-06 ----------------------------------------------------------------------------------------------------- 19 H 1 -2.126032920947007e+00 -2.126032946611870e+00 2.56649e-08 1.47821e-07 2 1.478435018591602e+00 1.478435014718944e+00 3.87266e-09 4.17105e-07 3 -4.365804898941314e-01 -4.365805237534230e-01 3.38593e-08 6.83034e-07 ----------------------------------------------------------------------------------------------------- 20 H 1 7.780983447084128e-01 7.780983251163803e-01 1.95920e-08 9.54026e-08 2 -1.852252024386150e+00 -1.852252041083915e+00 1.66978e-08 1.50338e-07 3 -1.954436217441871e-01 -1.954436199333953e-01 1.81079e-09 1.83535e-06 ----------------------------------------------------------------------------------------------------- 21 H 1 2.218026923161158e+00 2.218027029918638e+00 1.06757e-07 3.78148e-06 2 1.491174723631000e+00 1.491176231915755e+00 1.50828e-06 1.38427e-05 3 1.116835851866376e+00 1.116836500306022e+00 6.48440e-07 1.30648e-06 ----------------------------------------------------------------------------------------------------- 22 H 1 -7.043713343820611e-01 -7.043831837602862e-01 1.18494e-05 1.10120e-04 * 2 -1.587114196571534e+00 -1.587114326965973e+00 1.30394e-07 6.63870e-06 3 1.111097543550056e+00 1.111121108708915e+00 2.35652e-05 1.39220e-04 * ----------------------------------------------------------------------------------------------------- 23 H 1 -6.373213783575203e-01 -6.373214064290406e-01 2.80715e-08 3.74357e-08 2 1.950650525537065e-01 1.950651004496843e-01 4.78960e-08 2.44599e-07 3 -7.419719814516245e-01 -7.419720864045900e-01 1.04953e-07 2.08532e-07 ----------------------------------------------------------------------------------------------------- 24 H 1 -1.216794163641106e+00 -1.216794144882213e+00 1.87589e-08 4.51778e-07 2 9.286662585914353e-01 9.286663283519845e-01 6.97605e-08 8.91621e-08 3 -3.736436809662014e-01 -3.736436974005177e-01 1.64343e-08 6.96541e-07 ----------------------------------------------------------------------------------------------------- 25 H 1 3.764064146144053e+00 3.764064032407291e+00 1.13737e-07 1.47538e-06 2 -3.074091096773133e+00 -3.074091077482008e+00 1.92911e-08 1.34808e-06 3 -4.621789741594903e-01 -4.621789206566386e-01 5.35029e-08 1.68300e-07 ----------------------------------------------------------------------------------------------------- 26 H 1 -8.182228187731169e-02 -8.182262936435986e-02 3.47487e-07 5.03609e-07 2 1.355695738003844e+00 1.355696552339275e+00 8.14335e-07 3.39539e-05 * 3 6.070629205591360e-01 6.070625719886429e-01 3.48570e-07 9.61491e-07 ----------------------------------------------------------------------------------------------------- 27 H 1 4.617274176512736e-01 4.617275433234783e-01 1.25672e-07 1.68018e-06 2 1.207125852665440e-01 1.207125286586685e-01 5.66079e-08 6.30857e-07 3 1.103719404031142e+00 1.103719396500354e+00 7.53079e-09 7.81645e-07 ----------------------------------------------------------------------------------------------------- 28 H 1 6.545991370130764e-01 6.545991149545486e-01 2.20585e-08 1.29939e-07 2 -8.770520024520188e-01 -8.770519827464673e-01 1.97056e-08 6.02717e-07 3 8.507172027182694e-02 8.507170199892586e-02 1.82729e-08 2.41633e-07 ----------------------------------------------------------------------------------------------------- 29 H 1 4.069174064347518e-01 4.069182780875001e-01 8.71653e-07 1.46283e-06 2 1.753933317395117e+00 1.753933668051238e+00 3.50656e-07 1.33029e-05 3 1.853535317955899e+00 1.853534132039756e+00 1.18592e-06 4.94310e-06 ----------------------------------------------------------------------------------------------------- 30 H 1 7.889152992101446e-01 7.889152194352476e-01 7.97749e-08 3.09446e-07 2 -7.881884292491116e-02 -7.881766712302915e-02 1.17580e-06 1.85980e-06 3 9.847344220370222e-01 9.847344059711258e-01 1.60659e-08 1.20459e-05 ----------------------------------------------------------------------------------------------------- 31 H 1 -1.315175376205751e-01 -1.315176156650299e-01 7.80445e-08 1.12117e-07 2 -2.962443352334482e-01 -2.962422213762681e-01 2.11386e-06 2.54000e-06 3 -6.562697731524283e-01 -6.562688590309962e-01 9.14121e-07 8.90131e-06 ----------------------------------------------------------------------------------------------------- 32 H 1 4.834811530008391e-01 4.834811513043712e-01 1.69647e-09 2.35504e-06 2 -7.862044470433830e-01 -7.862044238830829e-01 2.31603e-08 6.62439e-07 3 -8.935690076466027e-01 -8.935691589640850e-01 1.51317e-07 1.63468e-06 ----------------------------------------------------------------------------------------------------- * Starred lines are suspected outliers and are not included when determining the error. A calculation is considered an outlier if it has an uncertainty that lies at an abnormal distance from the other uncertainties in this set of calculations. Outliers are determined using the box plot construction with fences. An outlier could indicate a problem with the the numerical differentiation or problems with the potential energy, such as discontinuities. Maximum error obtained for particle = 14, direction = 3: |F_model - F_numer| error = ----------------------- = 2.01282e-05 max{|F_numer|, eps} Comparison of analytical forces obtained from the model, the force computed as a numerical derivative of the energy, the difference between them, and the uncertainty in the numerical estimate of the force. The computed equilibrium lattice constant for the crystal (a0) is given in the heading. MONOATOMIC STRUCTURE -- Species = O | a0 = 2.6753597259521484 | (Configuration in file "config-O.xyz") ----------------------------------------------------------------------------------------------------- Part Spec Dir Force_model Force_numer |Force diff| uncertainty ----------------------------------------------------------------------------------------------------- 1 O 1 -1.599150830435909e+00 -1.599149997647265e+00 8.32789e-07 2.40698e-06 2 -1.141326689132194e+00 -1.141326923918049e+00 2.34786e-07 4.61310e-06 3 -2.830382822403254e+00 -2.830383855142014e+00 1.03274e-06 7.03337e-03 * ----------------------------------------------------------------------------------------------------- 2 O 1 -5.231629791829911e+00 -5.231629459558560e+00 3.32271e-07 4.17296e-06 2 -1.688581792120657e-01 -1.688579118578855e-01 2.67354e-07 6.56438e-06 3 -1.441451131220392e+00 -1.441451336858740e+00 2.05638e-07 3.14249e-06 ----------------------------------------------------------------------------------------------------- 3 O 1 6.582017555710391e+00 6.582017192930076e+00 3.62780e-07 6.16087e-06 2 -1.021327817132815e+01 -1.021327969814532e+01 1.52682e-06 8.42866e-06 3 1.268305510881899e+01 1.268305515388577e+01 4.50668e-08 4.04193e-06 ----------------------------------------------------------------------------------------------------- 4 O 1 -9.241454115299106e+00 -9.241453809869153e+00 3.05430e-07 8.55688e-06 2 -3.125117664159212e+00 -3.125117804417795e+00 1.40259e-07 5.97260e-06 3 -6.189060478473070e+00 -6.189060116372878e+00 3.62100e-07 5.25321e-02 * ----------------------------------------------------------------------------------------------------- 5 O 1 9.239870802615709e+00 9.239870320465222e+00 4.82150e-07 1.49764e-06 2 -4.182454124841539e+00 -4.182454273175100e+00 1.48334e-07 1.56479e-06 3 -9.248834969555508e+00 -9.248835241263585e+00 2.71708e-07 1.67138e-06 ----------------------------------------------------------------------------------------------------- 6 O 1 -1.172122850768205e-01 -1.172127214245619e-01 4.36348e-07 3.70307e-05 2 -1.167142834865408e+01 -1.167142843432617e+01 8.56721e-08 3.63578e-06 3 4.830521383735708e+00 4.830521370724824e+00 1.30109e-08 1.34127e-06 ----------------------------------------------------------------------------------------------------- 7 O 1 -2.019190153819203e-01 -2.019203461043708e-01 1.33072e-06 1.83884e-06 2 5.845330057133115e+00 5.845332932622860e+00 2.87549e-06 2.82828e-05 3 -6.897991485451882e+00 -6.897989973869295e+00 1.51158e-06 2.04897e-05 ----------------------------------------------------------------------------------------------------- 8 O 1 1.186832056682876e+00 1.186831644319268e+00 4.12364e-07 6.82237e-06 2 5.686229101862548e+00 5.686229077819179e+00 2.40434e-08 2.03780e-06 3 1.271464912312802e+01 1.271465000509102e+01 8.81963e-07 3.57015e-06 ----------------------------------------------------------------------------------------------------- 9 O 1 -3.656196242361507e+00 -3.656196305505063e+00 6.31436e-08 2.28887e-06 2 1.056618031711906e+00 1.056618213077419e+00 1.81366e-07 3.44227e-05 3 -5.737780685814641e+00 -5.737781892750063e+00 1.20694e-06 6.35104e-06 ----------------------------------------------------------------------------------------------------- 10 O 1 -3.845140809419255e+00 -3.845143147781596e+00 2.33836e-06 2.26094e-05 2 2.145115389567908e+00 2.145114397844569e+00 9.91723e-07 2.99644e-06 3 -5.782510614691246e+00 -5.782511556868807e+00 9.42178e-07 1.54325e-06 ----------------------------------------------------------------------------------------------------- 11 O 1 -9.743932481199838e+00 -9.743933201957635e+00 7.20758e-07 5.60711e-06 2 1.245293152852752e+01 1.245293279732891e+01 1.26880e-06 4.62113e-06 3 -3.134650139855637e+00 -3.134650288874913e+00 1.49019e-07 4.35083e-06 ----------------------------------------------------------------------------------------------------- 12 O 1 1.048167453099821e+01 1.048167426011493e+01 2.70883e-07 1.68321e-06 2 -1.125179648414459e+01 -1.125179672015828e+01 2.36014e-07 6.66284e-07 3 -3.879292697464040e+00 -3.879293249824549e+00 5.52361e-07 2.16392e-06 ----------------------------------------------------------------------------------------------------- 13 O 1 1.087791458382596e+01 1.087791557865425e+01 9.94828e-07 3.98215e-06 2 1.188997158848901e+01 1.188997266486198e+01 1.07637e-06 6.59482e-06 3 -3.971084167162898e+00 -3.971084845492786e+00 6.78330e-07 4.73247e-06 ----------------------------------------------------------------------------------------------------- 14 O 1 -4.401438143646928e+00 -4.401434418460581e+00 3.72519e-06 2.56927e-06 2 -1.211189917663944e+01 -1.211189922607059e+01 4.94312e-08 2.76376e-06 3 -2.069045914094102e+00 -2.069046186462504e+00 2.72368e-07 1.01666e-06 ----------------------------------------------------------------------------------------------------- 15 O 1 2.993407104949979e+00 2.993405833347875e+00 1.27160e-06 5.91641e-06 2 8.828661742361410e+00 8.828663054609544e+00 1.31225e-06 2.14605e-05 3 9.745894389128630e+00 9.745892950624317e+00 1.43850e-06 1.19958e-05 ----------------------------------------------------------------------------------------------------- 16 O 1 5.197890809659977e-01 5.197891477080587e-01 6.67421e-08 7.36638e-06 2 8.768168136859071e+00 8.768168005215472e+00 1.31644e-07 9.01901e-06 3 3.432988918412376e+00 3.432992056110535e+00 3.13770e-06 3.74874e-05 ----------------------------------------------------------------------------------------------------- 17 O 1 1.162482346052676e+00 1.162482231979991e+00 1.14073e-07 8.61012e-06 2 -8.631632936193439e+00 -8.631629422212146e+00 3.51398e-06 1.75563e-05 3 -1.152940525483565e+00 -1.152937677364797e+00 2.84812e-06 6.12919e-05 ----------------------------------------------------------------------------------------------------- 18 O 1 6.362740903314532e+00 6.362726428825769e+00 1.44745e-05 6.17303e-05 2 -7.399495909826581e-01 -7.399520824827196e-01 2.49150e-06 1.30022e-04 * 3 -3.247109155753505e+00 -3.247103736050641e+00 5.41970e-06 1.11067e-04 * ----------------------------------------------------------------------------------------------------- 19 O 1 4.874005072249278e+00 4.874003914651614e+00 1.15760e-06 2.31990e-06 2 7.549539878354281e+00 7.549534471351884e+00 5.40700e-06 6.93741e-05 * 3 -4.786215282267796e-01 -4.786272315260361e-01 5.70330e-06 2.46529e-05 ----------------------------------------------------------------------------------------------------- 20 O 1 4.659750910135095e+00 4.659753597457851e+00 2.68732e-06 5.64560e-05 2 3.113603135263215e+00 3.113603631382004e+00 4.96119e-07 1.82775e-05 3 1.577958368714440e+00 1.577942397167627e+00 1.59715e-05 7.89282e-05 * ----------------------------------------------------------------------------------------------------- 21 O 1 -9.727068075249445e+00 -9.727068010541586e+00 6.47079e-08 4.18956e-06 2 -9.203438735395794e-01 -9.203430148444290e-01 8.58695e-07 9.24578e-06 3 -9.335654691596689e+00 -9.335652886378673e+00 1.80522e-06 1.75277e-05 ----------------------------------------------------------------------------------------------------- 22 O 1 -3.155127137117335e+00 -3.155126486752353e+00 6.50365e-07 7.76951e-06 2 -1.910367011449624e+00 -1.910367663351491e+00 6.51902e-07 6.45905e-06 3 -6.445345008217124e+00 -6.445342978143029e+00 2.03007e-06 2.20005e-05 ----------------------------------------------------------------------------------------------------- 23 O 1 4.106903721188134e+00 4.106905656370450e+00 1.93518e-06 8.12116e-06 2 -2.886236239028158e+00 -2.886237996005744e+00 1.75698e-06 3.14517e-06 3 2.503132255628252e+00 2.503130425945193e+00 1.82968e-06 1.74419e-06 ----------------------------------------------------------------------------------------------------- 24 O 1 -6.912374290674564e+00 -6.912375183969645e+00 8.93295e-07 3.57934e-05 2 -3.354824674628261e+00 -3.354824247312764e+00 4.27315e-07 2.75936e-06 3 4.666203425058960e+00 4.666205232369955e+00 1.80731e-06 4.39792e-05 ----------------------------------------------------------------------------------------------------- 25 O 1 -9.025724397610290e+00 -9.025726495776844e+00 2.09817e-06 6.24254e-06 2 9.178560932793328e+00 9.178560488808712e+00 4.43985e-07 3.57587e-06 3 5.508864962974972e+00 5.508864687404965e+00 2.75570e-07 6.47578e-06 ----------------------------------------------------------------------------------------------------- 26 O 1 7.158390990303826e+00 7.158391748792839e+00 7.58489e-07 2.78328e-05 2 8.889105291720407e+00 8.889106215559854e+00 9.23839e-07 1.37003e-05 3 1.319357620208830e+01 1.319357583772922e+01 3.64359e-07 3.52857e-06 ----------------------------------------------------------------------------------------------------- 27 O 1 -1.946044994150778e-01 -1.946025648964190e-01 1.93452e-06 4.02861e-05 2 -7.495798700926037e+00 -7.495794867912569e+00 3.83301e-06 1.45060e-05 3 5.980668888238776e+00 5.980667382548266e+00 1.50569e-06 1.63060e-05 ----------------------------------------------------------------------------------------------------- 28 O 1 -2.601587725364086e+00 -2.601588597767368e+00 8.72403e-07 1.59081e-06 2 3.951095280100160e+00 3.951095063547307e+00 2.16553e-07 1.58114e-06 3 4.711488188402538e+00 4.711487595002447e+00 5.93400e-07 4.34527e-06 ----------------------------------------------------------------------------------------------------- 29 O 1 -3.854643416607398e+00 -3.854632884700349e+00 1.05319e-05 1.68180e-05 2 -4.762849350031165e+00 -4.762848363654089e+00 9.86377e-07 1.70896e-05 3 1.076257002162837e+00 1.076257498534497e+00 4.96372e-07 7.08375e-05 * ----------------------------------------------------------------------------------------------------- 30 O 1 1.131133460128261e+01 1.131133294841930e+01 1.65286e-06 1.77286e-05 2 -7.637759068196237e+00 -7.637759225961469e+00 1.57765e-07 5.91551e-06 3 -8.386540366601279e+00 -8.386540511859600e+00 1.45258e-07 2.39299e-06 ----------------------------------------------------------------------------------------------------- 31 O 1 1.835242986744497e+00 1.835240948395027e+00 2.03835e-06 1.26219e-06 2 -1.584094372314275e-01 -1.584064276839977e-01 3.00955e-06 6.53617e-06 3 1.864536433443248e+00 1.864540638899296e+00 4.20546e-06 1.60357e-05 ----------------------------------------------------------------------------------------------------- 32 O 1 -9.843153990330389e+00 -9.843154654050304e+00 6.63720e-07 3.27540e-06 2 3.009399625573983e+00 3.009398433893365e+00 1.19168e-06 3.21075e-05 3 -4.261498267870452e+00 -4.261498005309165e+00 2.62561e-07 1.27375e-05 ----------------------------------------------------------------------------------------------------- * Starred lines are suspected outliers and are not included when determining the error. A calculation is considered an outlier if it has an uncertainty that lies at an abnormal distance from the other uncertainties in this set of calculations. Outliers are determined using the box plot construction with fences. An outlier could indicate a problem with the the numerical differentiation or problems with the potential energy, such as discontinuities. Maximum error obtained for particle = 31, direction = 2: |F_model - F_numer| error = ----------------------- = 1.89989e-05 max{|F_numer|, eps} Comparison of analytical forces obtained from the model, the force computed as a numerical derivative of the energy, the difference between them, and the uncertainty in the numerical estimate of the force. The computed equilibrium lattice constant for the crystal (a0) is given in the heading. MONOATOMIC STRUCTURE -- Species = V | a0 = 3.8062543869018555 | (Configuration in file "config-V.xyz") ----------------------------------------------------------------------------------------------------- Part Spec Dir Force_model Force_numer |Force diff| uncertainty ----------------------------------------------------------------------------------------------------- 1 V 1 1.236818400066429e+00 1.236818397449306e+00 2.61712e-09 2.92558e-08 2 1.064483877061555e+00 1.064483875629727e+00 1.43183e-09 1.58331e-08 3 1.331186136633036e+00 1.331186134173490e+00 2.45955e-09 6.26181e-07 * ----------------------------------------------------------------------------------------------------- 2 V 1 -8.179829151383077e-01 -8.179829149945016e-01 1.43806e-10 1.70657e-08 2 5.556328921364101e-01 5.556328922924990e-01 1.56089e-10 1.73542e-08 3 8.201287504735014e-01 8.201287463780907e-01 4.09541e-09 1.75033e-08 ----------------------------------------------------------------------------------------------------- 3 V 1 -3.344272277568807e+00 -3.344272278063349e+00 4.94542e-10 9.83709e-09 2 2.186333328668939e+00 2.186333328516468e+00 1.52471e-10 5.79315e-09 3 9.747311591097568e-03 9.747306413208289e-03 5.17789e-09 5.24308e-09 ----------------------------------------------------------------------------------------------------- 4 V 1 -1.272298598259148e-01 -1.272298602877867e-01 4.61872e-10 1.38946e-08 2 -1.787072771689490e+00 -1.787072771670470e+00 1.90195e-11 1.31911e-08 3 7.287618446634054e-01 7.287618427771344e-01 1.88627e-09 1.70964e-08 ----------------------------------------------------------------------------------------------------- 5 V 1 1.072803234072963e+00 1.072803236746732e+00 2.67377e-09 1.23373e-08 2 8.471851137116805e-01 8.471851099807141e-01 3.73097e-09 3.42236e-08 3 2.761719394308597e-01 2.761719381583890e-01 1.27247e-09 1.63793e-08 ----------------------------------------------------------------------------------------------------- 6 V 1 -1.558104002795334e-02 -1.558104149383511e-02 1.46588e-09 4.30616e-08 2 -1.992941472287147e+00 -1.992941472518901e+00 2.31753e-10 1.64250e-08 3 -2.290450312665964e+00 -2.290450313408299e+00 7.42336e-10 3.69039e-08 ----------------------------------------------------------------------------------------------------- 7 V 1 8.613484117107570e-02 8.613484218920978e-02 1.01813e-09 2.01762e-08 2 1.026291811694606e+00 1.026291806403227e+00 5.29138e-09 1.42343e-08 3 -3.765528957084033e-01 -3.765528936667998e-01 2.04160e-09 2.77689e-08 ----------------------------------------------------------------------------------------------------- 8 V 1 -3.680187098855655e+00 -3.680187098955503e+00 9.98477e-11 1.67127e-09 2 -2.345569344494591e+00 -2.345569345322760e+00 8.28170e-10 5.92922e-09 3 2.143585072205023e+00 2.143585072248472e+00 4.34497e-11 3.75509e-09 ----------------------------------------------------------------------------------------------------- 9 V 1 4.300794206901682e-01 4.300794214126780e-01 7.22510e-10 4.07195e-08 2 -1.118536636292700e+00 -1.118536632391353e+00 3.90135e-09 7.63999e-09 3 2.043019282394235e-01 2.043019258093907e-01 2.43003e-09 9.23654e-09 ----------------------------------------------------------------------------------------------------- 10 V 1 1.490072216375582e+00 1.490072219066288e+00 2.69071e-09 1.14393e-08 2 2.420372777663660e+00 2.420372776244569e+00 1.41909e-09 1.06232e-08 3 -1.917743571406886e+00 -1.917743572002984e+00 5.96098e-10 3.74219e-07 * ----------------------------------------------------------------------------------------------------- 11 V 1 -2.366169199541178e+00 -2.366169199870478e+00 3.29300e-10 4.73601e-09 2 -2.395765647215116e+00 -2.395765646049682e+00 1.16543e-09 3.52670e-09 3 3.926855773037946e+00 3.926855772536209e+00 5.01737e-10 2.11420e-09 ----------------------------------------------------------------------------------------------------- 12 V 1 -7.961261760239516e-01 -7.961261805565635e-01 4.53261e-09 3.17780e-08 2 -1.539014902062953e-02 -1.539015032105198e-02 1.30042e-09 4.10407e-08 3 1.164397331297920e+00 1.164397330235429e+00 1.06249e-09 6.12553e-09 ----------------------------------------------------------------------------------------------------- 13 V 1 2.949990744466303e+00 2.949990746681013e+00 2.21471e-09 7.09920e-09 2 5.040263996946133e-01 5.040263993313374e-01 3.63276e-10 4.25136e-09 3 -2.732758300984427e+00 -2.732758298965312e+00 2.01912e-09 1.94752e-08 ----------------------------------------------------------------------------------------------------- 14 V 1 -1.257776717930858e+00 -1.257776716831080e+00 1.09978e-09 3.64553e-08 2 -1.267351650859141e+00 -1.267351653653094e+00 2.79395e-09 6.32832e-09 3 1.051054730025635e+00 1.051054724166802e+00 5.85883e-09 5.19537e-09 ----------------------------------------------------------------------------------------------------- 15 V 1 4.848014677727019e-01 4.848014666863310e-01 1.08637e-09 1.68209e-08 2 4.268426333653006e+00 4.268426327865334e+00 5.78767e-09 4.91070e-08 3 2.420947534627875e+00 2.420947537516115e+00 2.88824e-09 1.08244e-08 ----------------------------------------------------------------------------------------------------- 16 V 1 1.410980680969932e+00 1.410980679169253e+00 1.80068e-09 3.49250e-08 2 8.126228728963680e-01 8.126228704581261e-01 2.43824e-09 1.63359e-08 3 -2.102741998099563e+00 -2.102742001395484e+00 3.29592e-09 1.82881e-08 ----------------------------------------------------------------------------------------------------- 17 V 1 2.293375343117487e+00 2.293375340939489e+00 2.17800e-09 9.04042e-09 2 1.589287476137578e+00 1.589287473868798e+00 2.26878e-09 7.02419e-09 3 -1.095881189974912e+00 -1.095881189108177e+00 8.66735e-10 1.64510e-06 * ----------------------------------------------------------------------------------------------------- 18 V 1 9.362461526334571e-01 9.362461524869515e-01 1.46506e-10 1.12329e-09 2 -2.487151307147803e+00 -2.487151308460669e+00 1.31287e-09 7.14041e-09 3 1.745984708997371e+00 1.745984708554263e+00 4.43108e-10 4.32628e-09 ----------------------------------------------------------------------------------------------------- 19 V 1 -8.820411443296574e-01 -8.820411459158310e-01 1.58617e-09 2.60837e-08 2 1.311680068832415e+00 1.311680067474232e+00 1.35818e-09 3.13951e-08 3 -1.608531776367125e+00 -1.608531775402344e+00 9.64781e-10 2.38694e-08 ----------------------------------------------------------------------------------------------------- 20 V 1 1.653653588294270e+00 1.653653586784715e+00 1.50956e-09 2.87330e-08 2 4.775335022011706e-01 4.775335015816161e-01 6.19555e-10 6.95718e-09 3 -1.475118911366752e+00 -1.475118918340231e+00 6.97348e-09 2.13290e-09 ----------------------------------------------------------------------------------------------------- 21 V 1 -2.144740802772881e-01 -2.144740803088347e-01 3.15466e-11 5.04772e-09 2 1.716015091742733e+00 1.716015091485488e+00 2.57245e-10 8.09799e-09 3 2.133852267165524e+00 2.133852267466943e+00 3.01420e-10 8.84444e-09 ----------------------------------------------------------------------------------------------------- 22 V 1 4.997513775989209e-01 4.997513728324873e-01 4.76643e-09 1.27415e-08 2 2.259553112943092e-02 2.259553162166328e-02 4.92232e-10 2.23071e-08 3 6.066919674394597e-01 6.066919670585110e-01 3.80949e-10 1.21171e-08 ----------------------------------------------------------------------------------------------------- 23 V 1 1.815659150710353e-01 1.815659153814204e-01 3.10385e-10 1.57830e-08 2 -6.266643004057089e-01 -6.266643015331841e-01 1.12748e-09 2.57785e-08 3 -1.148786717138462e+00 -1.148786711737005e+00 5.40146e-09 1.40715e-08 ----------------------------------------------------------------------------------------------------- 24 V 1 -1.740394734438004e+00 -1.740394735347566e+00 9.09562e-10 5.82671e-09 2 1.461985132315852e+00 1.461985131683337e+00 6.32515e-10 9.02193e-09 3 -1.691867747159899e+00 -1.691867749550325e+00 2.39043e-09 1.63719e-08 ----------------------------------------------------------------------------------------------------- 25 V 1 1.525841746066325e+00 1.525841745707835e+00 3.58490e-10 3.00610e-07 * 2 -3.453935421311355e+00 -3.453935421945689e+00 6.34334e-10 3.59467e-09 3 -2.226993187524564e+00 -2.226993188987321e+00 1.46276e-09 1.93884e-08 ----------------------------------------------------------------------------------------------------- 26 V 1 1.562949170854197e+00 1.562949172429076e+00 1.57488e-09 9.77946e-09 2 -1.098857300421076e+00 -1.098857301965003e+00 1.54393e-09 2.72308e-09 3 -2.384651277664441e-01 -2.384651273770322e-01 3.89412e-10 4.17236e-08 ----------------------------------------------------------------------------------------------------- 27 V 1 -4.756639009549485e-01 -4.756639018736352e-01 9.18687e-10 2.82860e-08 2 3.624654268720210e-01 3.624654271678968e-01 2.95876e-10 1.44316e-08 3 -7.649677586206916e-01 -7.649677592499894e-01 6.29298e-10 1.45583e-08 ----------------------------------------------------------------------------------------------------- 28 V 1 -9.086265993529010e-01 -9.086266045756566e-01 5.22276e-09 1.86449e-08 2 1.440818937043707e+00 1.440818937004928e+00 3.87781e-11 8.25697e-09 3 3.347487936698132e+00 3.347487934261629e+00 2.43650e-09 7.22144e-09 ----------------------------------------------------------------------------------------------------- 29 V 1 -2.190275786016737e+00 -2.190275778266277e+00 7.75046e-09 4.58176e-08 2 -4.298110299312266e+00 -4.298110298798672e+00 5.13594e-10 2.13912e-08 3 -3.411919533068785e+00 -3.411919534376832e+00 1.30805e-09 6.41118e-09 ----------------------------------------------------------------------------------------------------- 30 V 1 2.443335990157602e+00 2.443335989369604e+00 7.87998e-10 4.12224e-09 2 1.868971850951559e+00 1.868971852410836e+00 1.45928e-09 2.13963e-08 3 -3.333621418088839e+00 -3.333621419844787e+00 1.75595e-09 1.20281e-08 ----------------------------------------------------------------------------------------------------- 31 V 1 3.181539560568025e+00 3.181539559250096e+00 1.31793e-09 1.19164e-08 2 -5.161998199610519e+00 -5.161998198668999e+00 9.41520e-10 2.71631e-08 3 2.660146550501251e+00 2.660146556706000e+00 6.20475e-09 1.02085e-08 ----------------------------------------------------------------------------------------------------- 32 V 1 -4.623138319664317e+00 -4.623138323599081e+00 3.93476e-09 1.14073e-08 2 4.112616075660237e+00 4.112616073850308e+00 1.80993e-09 2.19258e-08 3 1.845098662914252e+00 1.845098661728960e+00 1.18529e-09 2.29825e-08 ----------------------------------------------------------------------------------------------------- * Starred lines are suspected outliers and are not included when determining the error. A calculation is considered an outlier if it has an uncertainty that lies at an abnormal distance from the other uncertainties in this set of calculations. Outliers are determined using the box plot construction with fences. An outlier could indicate a problem with the the numerical differentiation or problems with the potential energy, such as discontinuities. Maximum error obtained for particle = 3, direction = 3: |F_model - F_numer| error = ----------------------- = 5.31212e-07 max{|F_numer|, eps} Comparison of analytical forces obtained from the model, the force computed as a numerical derivative of the energy, the difference between them, and the uncertainty in the numerical estimate of the force. The computed equilibrium lattice constant for the crystal (a0) is given in the heading. MIXED STRUCTURE -- Species = C H O V | a0 = 3.15543736694126 | (Configuration in file "config-CHOV.xyz") ----------------------------------------------------------------------------------------------------- Part Spec Dir Force_model Force_numer |Force diff| uncertainty ----------------------------------------------------------------------------------------------------- 1 C 1 2.085384257588141e-01 2.085390578644525e-01 6.32106e-07 3.33078e-05 2 4.236546926816691e+00 4.236549412251835e+00 2.48544e-06 2.87726e-05 3 3.010969397103258e+00 3.010988093632846e+00 1.86965e-05 3.63628e-05 ----------------------------------------------------------------------------------------------------- 2 H 1 1.223921227188162e-01 1.223552301619095e-01 3.68926e-05 5.24076e-03 2 6.993747462667040e-01 6.990686067698111e-01 3.06139e-04 8.67663e-04 3 1.533178668312799e+00 1.533525740164541e+00 3.47072e-04 3.51204e-04 ----------------------------------------------------------------------------------------------------- 3 H 1 6.295029422181494e-01 6.297206439607290e-01 2.17702e-04 1.29869e-04 2 -9.606042796989132e-01 -9.615214490448138e-01 9.17169e-04 2.61435e-05 3 1.066825602731800e+00 1.066861379284251e+00 3.57766e-05 9.93699e-05 ----------------------------------------------------------------------------------------------------- 4 O 1 -2.661864387656782e-01 -2.663494972485123e-01 1.63058e-04 2.41040e-03 2 -4.838445868657064e+00 -4.838407015616524e+00 3.88530e-05 6.55779e-05 3 -1.215784490869722e+00 -1.215189278239664e+00 5.95213e-04 1.49948e-04 ----------------------------------------------------------------------------------------------------- 5 H 1 1.831624667138100e+00 1.831631278784424e+00 6.61165e-06 5.40490e-03 2 1.633775656978551e+00 1.633447921789126e+00 3.27735e-04 1.44807e-03 3 7.400767251628310e-01 7.409110383832778e-01 8.34313e-04 4.35062e-05 ----------------------------------------------------------------------------------------------------- 6 C 1 -1.410741367518903e+00 -1.410585791230420e+00 1.55576e-04 4.86217e-04 2 -2.882947772097685e+00 -2.883192835419439e+00 2.45063e-04 1.70724e-04 3 6.560859389038144e-01 6.562074574904826e-01 1.21519e-04 1.34570e-04 ----------------------------------------------------------------------------------------------------- 7 V 1 4.904758625025494e+00 4.904882154316834e+00 1.23529e-04 1.76200e-03 2 -1.579717876426940e+00 -1.576652513099170e+00 3.06536e-03 1.32536e-03 3 -5.619259555474749e+00 -5.620525141878618e+00 1.26559e-03 9.63238e-03 ----------------------------------------------------------------------------------------------------- 8 C 1 4.383763436115045e-01 4.382056725907710e-01 1.70671e-04 3.30391e-04 2 3.257940586220940e+00 3.256713822938044e+00 1.22676e-03 1.45172e-03 3 -2.916513497947412e+00 -2.916319819276621e+00 1.93679e-04 2.60515e-04 ----------------------------------------------------------------------------------------------------- 9 C 1 6.515100814455113e+00 6.515112312317733e+00 1.14979e-05 3.89207e-05 2 3.138191611915904e+00 3.138190730347910e+00 8.81568e-07 2.85290e-05 3 2.301140776376645e+00 2.301182159031751e+00 4.13827e-05 3.83582e-05 ----------------------------------------------------------------------------------------------------- 10 C 1 -3.107801862276215e+00 -3.107776174602037e+00 2.56877e-05 7.27099e-05 2 -4.957202932740434e+00 -4.957229591282727e+00 2.66585e-05 2.29005e-04 3 8.199365411088535e-01 8.198602607590262e-01 7.62803e-05 1.53153e-05 ----------------------------------------------------------------------------------------------------- 11 H 1 -2.973294391774382e+00 -2.973264820362359e+00 2.95714e-05 4.61063e-03 2 -9.811971488794732e-02 -9.812690648160569e-02 7.19159e-06 6.85135e-03 3 -2.074682568686470e+00 -2.074667455709199e+00 1.51130e-05 3.76130e-04 ----------------------------------------------------------------------------------------------------- 12 H 1 6.655589986170944e-02 6.656233997470566e-02 6.44011e-06 1.87422e-04 2 -2.706420687196502e+00 -2.706661131690504e+00 2.40444e-04 4.33361e-05 3 2.137610456571217e+00 2.137768727597763e+00 1.58271e-04 1.64603e-05 ----------------------------------------------------------------------------------------------------- 13 H 1 -7.632374754251974e-01 -7.631587390784259e-01 7.87363e-05 2.14135e-04 2 1.616936675956933e+00 1.616907032135277e+00 2.96438e-05 8.95476e-06 3 5.573341146380578e-01 5.575462086269176e-01 2.12094e-04 6.76907e-04 ----------------------------------------------------------------------------------------------------- 14 V 1 -2.100817027652166e+00 -2.101366955152202e+00 5.49928e-04 1.43349e-05 2 -2.484402644733700e+00 -2.485164663462045e+00 7.62019e-04 4.79345e-04 3 3.465314004457126e+00 3.465466402777770e+00 1.52398e-04 2.55567e-03 ----------------------------------------------------------------------------------------------------- 15 O 1 8.744862935947304e-01 8.749896460696637e-01 5.03352e-04 2.17580e-04 2 2.878401588026537e+00 2.878708817346688e+00 3.07229e-04 6.97830e-04 3 -2.551716587611326e+00 -2.552393030465939e+00 6.76443e-04 2.13372e-04 ----------------------------------------------------------------------------------------------------- 16 O 1 2.660096230949422e+00 2.660325822500903e+00 2.29592e-04 6.17648e-05 2 4.857792200874462e-01 4.861299892831943e-01 3.50769e-04 9.66918e-05 3 6.049403275777621e-01 6.048038202323739e-01 1.36507e-04 8.90123e-04 ----------------------------------------------------------------------------------------------------- 17 V 1 4.756642732108314e+00 4.757676840516479e+00 1.03411e-03 4.45501e-03 2 -1.475860524715663e+01 -1.475806657983217e+01 5.38667e-04 2.70233e-03 3 -1.579426252582964e+01 -1.579440164076023e+01 1.39115e-04 3.48696e-05 ----------------------------------------------------------------------------------------------------- 18 O 1 -3.344587165785422e+00 -3.345157648617048e+00 5.70483e-04 2.75686e-03 2 1.964884590874370e+00 1.966067071195420e+00 1.18248e-03 1.00433e-02 * 3 2.688108077209897e+00 2.686364659394912e+00 1.74342e-03 3.43830e-02 * ----------------------------------------------------------------------------------------------------- 19 O 1 -2.075217419225757e+00 -2.075077018721714e+00 1.40401e-04 2.48875e-02 * 2 2.511398571544162e+00 2.515856845699281e+00 4.45827e-03 1.21899e-02 * 3 -2.109670828410017e+00 -2.109714562104577e+00 4.37337e-05 6.77004e-03 ----------------------------------------------------------------------------------------------------- 20 V 1 -2.851974466263338e+00 -2.852492558638695e+00 5.18092e-04 2.12823e-03 2 1.000276917060409e+01 1.000274204352586e+01 2.71271e-05 6.79701e-04 3 1.665931473580327e+01 1.665955861263753e+01 2.43877e-04 2.36804e-04 ----------------------------------------------------------------------------------------------------- 21 V 1 -1.796837424044015e+00 -1.796943020559687e+00 1.05597e-04 3.67756e-03 2 -3.257102669078769e+00 -3.255233783457654e+00 1.86889e-03 6.20853e-04 3 2.019080359965046e+00 2.018826092351514e+00 2.54268e-04 4.48698e-03 ----------------------------------------------------------------------------------------------------- 22 C 1 -1.772824904643407e+00 -1.772682301342269e+00 1.42603e-04 2.13958e-04 2 -3.563358303512176e+00 -3.563996628056310e+00 6.38325e-04 3.68405e-03 3 -4.692074076534205e+00 -4.692176600363746e+00 1.02524e-04 1.07613e-02 * ----------------------------------------------------------------------------------------------------- 23 C 1 -2.673912161326963e+00 -2.675466444490915e+00 1.55428e-03 4.20719e-03 2 -2.985030131221589e-01 -2.984388975523494e-01 6.41156e-05 1.05142e-03 3 -4.596897592418940e-01 -4.585126940056222e-01 1.17707e-03 2.56772e-03 ----------------------------------------------------------------------------------------------------- 24 V 1 6.298217442548293e+00 6.296630160808993e+00 1.58728e-03 9.66025e-03 2 1.635401874051733e+00 1.635719843740103e+00 3.17970e-04 8.98685e-03 3 5.187078515566296e+00 5.183811261865527e+00 3.26725e-03 1.98315e-04 ----------------------------------------------------------------------------------------------------- 25 V 1 -4.889373668216349e+00 -4.889295060732355e+00 7.86075e-05 3.91179e-04 2 9.060704798036742e+00 9.061488066537045e+00 7.83269e-04 2.17068e-03 3 -1.188510923093198e+01 -1.188450988651433e+01 5.99344e-04 1.37832e-03 ----------------------------------------------------------------------------------------------------- 26 H 1 -1.602170184073531e+00 -1.602226569853802e+00 5.63858e-05 1.24555e-04 2 -3.555363529877694e+00 -3.555458570846172e+00 9.50410e-05 1.50436e-03 3 1.921444887953409e+00 1.921456471726988e+00 1.15838e-05 1.51807e-04 ----------------------------------------------------------------------------------------------------- 27 V 1 8.250820346601307e+00 8.251269924876627e+00 4.49578e-04 1.48384e-04 2 7.986932153583928e+00 7.986641386184407e+00 2.90767e-04 6.47285e-04 3 3.418433735059011e+00 3.417958677986083e+00 4.75057e-04 1.09134e-03 ----------------------------------------------------------------------------------------------------- 28 O 1 1.538873656879421e+00 1.539116518752055e+00 2.42862e-04 5.35392e-03 2 -4.143264235756743e+00 -4.145622233379319e+00 2.35800e-03 1.82100e-05 3 -1.564695207281315e+00 -1.564820146908914e+00 1.24940e-04 2.59413e-04 ----------------------------------------------------------------------------------------------------- 29 O 1 4.851070235793314e-01 4.845385327697729e-01 5.68491e-04 7.42580e-04 2 2.001147846819050e+00 2.001140175883962e+00 7.67094e-06 2.26330e-03 3 1.886640851249670e+00 1.887582412229925e+00 9.41561e-04 6.66649e-04 ----------------------------------------------------------------------------------------------------- 30 C 1 -5.718884954542509e+00 -5.719240819563407e+00 3.55865e-04 2.60263e-04 2 -2.774311004285134e+00 -2.775133492912990e+00 8.22489e-04 1.33010e-04 3 -6.591940686509519e-01 -6.592280646224511e-01 3.39960e-05 4.29780e-04 ----------------------------------------------------------------------------------------------------- 31 O 1 -1.730493384488409e+00 -1.731441737060331e+00 9.48353e-04 1.78284e-03 2 1.185235694102472e+00 1.185017509619210e+00 2.18184e-04 7.53159e-03 3 2.163732147573009e+00 2.164006603406106e+00 2.74456e-04 6.59877e-03 ----------------------------------------------------------------------------------------------------- 32 H 1 -5.027392710262809e-01 -5.028150435916070e-01 7.57726e-05 7.02884e-05 2 -1.437051932657774e+00 -1.437063255312673e+00 1.13227e-05 1.31803e-05 3 -1.294593465854089e+00 -1.294501553796540e+00 9.19121e-05 1.95713e-04 ----------------------------------------------------------------------------------------------------- * Starred lines are suspected outliers and are not included when determining the error. A calculation is considered an outlier if it has an uncertainty that lies at an abnormal distance from the other uncertainties in this set of calculations. Outliers are determined using the box plot construction with fences. An outlier could indicate a problem with the the numerical differentiation or problems with the potential energy, such as discontinuities. Maximum error obtained for particle = 23, direction = 3: |F_model - F_numer| error = ----------------------- = 2.56714e-03 max{|F_numer|, eps} ===================================================================================================== Grade is based on maximum error across all systems (for a model supporting multiple species): Maximum error = 2.56714e-03 Machine precision = 2.22045e-16 Note: Lower grades can reflect errors in the analytical forces, a lack of smoothness in the potential, or loss of precision in the numerical derivative computation. Grade: D Comment: The letter grade D was assigned because the normalized error in the computation was 2.56714e-03 compared with a machine precision of 2.22045e-16. The letter grade was based on 'score=log10(error/eps)', with ranges A=[0, 7.5], B=(7.5, 10.0], C=(10.0, 12.5], D=(12.5, 15.0), F>15.0. 'A' is the best grade, and 'F' indicates failure.