!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!! !!!!! !!!!! 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 : MEAM_LAMMPS_JeongLee_2020_PdC__MO_068985622065_002 Supported species : C Pd 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 = 3.379122734069824 | (Configuration in file "config-C.xyz") ----------------------------------------------------------------------------------------------------- Part Spec Dir Force_model Force_numer |Force diff| uncertainty ----------------------------------------------------------------------------------------------------- 1 C 1 1.092309830761578e+00 1.092309830772646e+00 1.10676e-11 1.10298e-10 2 1.319831826400909e+00 1.319831826394858e+00 6.05138e-12 2.10776e-10 3 1.328900853730179e+00 1.328900853719571e+00 1.06084e-11 1.02235e-10 ----------------------------------------------------------------------------------------------------- 2 C 1 1.492183345020174e-01 1.492183345100956e-01 8.07826e-12 1.67503e-10 2 1.026998445884790e-01 1.026998446215088e-01 3.30298e-11 8.04283e-11 3 -2.288644397149875e+00 -2.288644397166734e+00 1.68585e-11 1.57720e-10 ----------------------------------------------------------------------------------------------------- 3 C 1 -1.448196704790096e+00 -1.448196704773935e+00 1.61606e-11 9.36122e-11 2 -8.832314050943950e-01 -8.832314050510576e-01 4.33374e-11 2.80575e-10 3 1.658458342149162e-01 1.658458342282808e-01 1.33647e-11 1.64759e-10 ----------------------------------------------------------------------------------------------------- 4 C 1 -3.628311870237289e+00 -3.628311870186168e+00 5.11218e-11 1.41839e-10 2 -4.591922352379775e-01 -4.591922352594069e-01 2.14294e-11 3.03924e-10 3 -1.207589798044527e+00 -1.207589798060598e+00 1.60709e-11 3.75179e-10 ----------------------------------------------------------------------------------------------------- 5 C 1 6.858700793924168e-01 6.858700794102298e-01 1.78130e-11 5.88968e-10 2 -5.153106640680629e-01 -5.153106641331037e-01 6.50408e-11 5.62348e-10 3 -2.258894168046437e+00 -2.258894168016278e+00 3.01594e-11 2.56421e-10 ----------------------------------------------------------------------------------------------------- 6 C 1 2.110833937237838e-01 2.110833936347595e-01 8.90242e-11 4.30190e-10 2 -1.117054757952925e+00 -1.117054758023488e+00 7.05633e-11 2.96431e-10 3 -4.880986686415933e-01 -4.880986686293267e-01 1.22666e-11 3.30557e-10 ----------------------------------------------------------------------------------------------------- 7 C 1 3.431911963438815e-01 3.431911963110958e-01 3.27857e-11 1.02238e-10 2 -2.569399848199392e+00 -2.569399848166438e+00 3.29541e-11 1.53027e-09 3 -7.952248793086802e-01 -7.952248792880610e-01 2.06193e-11 5.75070e-10 ----------------------------------------------------------------------------------------------------- 8 C 1 5.473813321422672e-02 5.473813327558651e-02 6.13598e-11 7.50021e-09 * 2 1.563464211448300e+00 1.563464211549311e+00 1.01010e-10 3.08455e-10 3 7.934802631890240e-01 7.934802630930697e-01 9.59542e-11 4.30942e-10 ----------------------------------------------------------------------------------------------------- 9 C 1 -3.464296227710182e+00 -3.464296227698080e+00 1.21019e-11 2.41703e-10 2 -3.424179763805477e+00 -3.424179763811199e+00 5.72165e-12 1.19032e-10 3 5.112569275362000e-02 5.112569271763272e-02 3.59873e-11 3.66258e-10 ----------------------------------------------------------------------------------------------------- 10 C 1 -1.408664868571794e+00 -1.408664868549291e+00 2.25036e-11 7.12364e-11 2 7.681547630367758e+00 7.681547630379659e+00 1.19007e-11 1.52430e-10 3 -4.728738609288344e+00 -4.728738609310578e+00 2.22338e-11 1.50411e-09 ----------------------------------------------------------------------------------------------------- 11 C 1 -3.218321639478736e+00 -3.218321639504406e+00 2.56701e-11 2.67615e-10 2 -4.704251417317991e+00 -4.704251417299543e+00 1.84484e-11 4.89378e-10 3 -1.919712300068316e+00 -1.919712300006247e+00 6.20684e-11 5.26688e-10 ----------------------------------------------------------------------------------------------------- 12 C 1 -2.993055346980078e+00 -2.993055346976663e+00 3.41460e-12 3.56994e-11 2 1.325975401226325e+00 1.325975401249105e+00 2.27796e-11 6.88170e-11 3 9.093391701767757e-01 9.093391701663582e-01 1.04174e-11 4.49003e-10 ----------------------------------------------------------------------------------------------------- 13 C 1 4.590136086008497e+00 4.590136085976856e+00 3.16414e-11 5.44447e-10 2 -1.663281241978062e+00 -1.663281242013247e+00 3.51847e-11 5.20265e-09 * 3 -7.883992620333652e+00 -7.883992620280934e+00 5.27178e-11 4.94020e-10 ----------------------------------------------------------------------------------------------------- 14 C 1 1.546639148260378e+00 1.546639148284782e+00 2.44043e-11 6.16638e-11 2 1.382080238482530e-01 1.382080238557420e-01 7.48893e-12 1.17310e-10 3 -1.453227625699374e+00 -1.453227625704399e+00 5.02420e-12 2.83340e-11 ----------------------------------------------------------------------------------------------------- 15 C 1 8.683784953615403e+00 8.683784953656183e+00 4.07798e-11 4.97240e-10 2 -5.423129967584334e+00 -5.423129967568961e+00 1.53726e-11 4.54243e-10 3 -3.414020365718696e-01 -3.414020366355234e-01 6.36537e-11 4.01655e-10 ----------------------------------------------------------------------------------------------------- 16 C 1 -5.498095065224441e+00 -5.498095065214004e+00 1.04370e-11 1.79694e-10 2 1.051576874221785e+01 1.051576874230447e+01 8.66187e-11 8.10935e-10 3 6.363586992346232e+00 6.363586992311068e+00 3.51639e-11 1.09054e-10 ----------------------------------------------------------------------------------------------------- 17 C 1 -2.888011844567568e+00 -2.888011844566790e+00 7.78044e-13 3.13349e-11 2 -3.584056113560441e+00 -3.584056113557613e+00 2.82796e-12 9.98179e-10 3 -1.094685389439817e+00 -1.094685389385390e+00 5.44269e-11 1.61106e-08 * ----------------------------------------------------------------------------------------------------- 18 C 1 1.281709217754617e+00 1.281709217688285e+00 6.63321e-11 3.76281e-10 2 2.089589244584773e+00 2.089589244501123e+00 8.36495e-11 4.74874e-10 3 -1.760938972461792e+00 -1.760938972509566e+00 4.77740e-11 3.39455e-10 ----------------------------------------------------------------------------------------------------- 19 C 1 7.049561701727833e-01 7.049561700077323e-01 1.65051e-10 1.26717e-10 2 -2.849817364308621e+00 -2.849817364263389e+00 4.52327e-11 4.84583e-10 3 2.486429246392488e+00 2.486429246387603e+00 4.88543e-12 1.53013e-10 ----------------------------------------------------------------------------------------------------- 20 C 1 -2.087024485161280e+00 -2.087024485156936e+00 4.34408e-12 2.96097e-10 2 1.905507099197575e+00 1.905507099318834e+00 1.21259e-10 5.06043e-10 3 2.882481472440060e+00 2.882481472366223e+00 7.38374e-11 2.13777e-10 ----------------------------------------------------------------------------------------------------- 21 C 1 8.950419265577767e-01 8.950419265744980e-01 1.67213e-11 1.80105e-10 2 -4.380691683427155e+00 -4.380691683407482e+00 1.96723e-11 5.17534e-10 3 2.298247779319476e+00 2.298247779363215e+00 4.37392e-11 2.31993e-10 ----------------------------------------------------------------------------------------------------- 22 C 1 6.172679963637309e+00 6.172679963608412e+00 2.88969e-11 1.52664e-10 2 2.886466150278788e+00 2.886466150275025e+00 3.76277e-12 6.87411e-10 3 3.606962580263052e+00 3.606962580247949e+00 1.51035e-11 2.20431e-10 ----------------------------------------------------------------------------------------------------- 23 C 1 -1.080359095471988e+00 -1.080359095529962e+00 5.79741e-11 3.85982e-10 2 4.007065522212708e-01 4.007065522097684e-01 1.15024e-11 1.07530e-10 3 -1.587128474382808e+00 -1.587128474291059e+00 9.17488e-11 5.32274e-11 ----------------------------------------------------------------------------------------------------- 24 C 1 2.741909983231208e-01 2.741909982937535e-01 2.93673e-11 1.31372e-10 2 5.072927910943383e-01 5.072927911782833e-01 8.39450e-11 3.13842e-10 3 6.379951459816882e-01 6.379951459935250e-01 1.18369e-11 2.18723e-10 ----------------------------------------------------------------------------------------------------- 25 C 1 1.450094253463424e-01 1.450094253456683e-01 6.74100e-13 6.60135e-11 2 -1.320626123330087e+00 -1.320626123329089e+00 9.98313e-13 2.09909e-10 3 -1.732980176939414e+00 -1.732980176938621e+00 7.92921e-13 5.46104e-10 ----------------------------------------------------------------------------------------------------- 26 C 1 1.970944217322042e+00 1.970944217281858e+00 4.01847e-11 1.21971e-09 2 2.387488411099212e+00 2.387488411103014e+00 3.80185e-12 2.64175e-11 3 -2.776032599080335e+00 -2.776032599064033e+00 1.63021e-11 2.45415e-10 ----------------------------------------------------------------------------------------------------- 27 C 1 -4.445179751047403e+00 -4.445179751067007e+00 1.96039e-11 1.56406e-09 * 2 7.159566356237588e-01 7.159566356464281e-01 2.26693e-11 9.51292e-10 3 2.367241190755745e+00 2.367241190760277e+00 4.53149e-12 3.32429e-10 ----------------------------------------------------------------------------------------------------- 28 C 1 -2.436705572942337e+00 -2.436705573030499e+00 8.81628e-11 1.26396e-10 2 2.496988806146537e-01 2.496988806094667e-01 5.18702e-12 2.10273e-10 3 4.099964366343430e+00 4.099964366290188e+00 5.32419e-11 3.70690e-10 ----------------------------------------------------------------------------------------------------- 29 C 1 9.263003963243019e-02 9.263003968184508e-02 4.94149e-11 1.17569e-09 2 -3.175489440176560e+00 -3.175489440171847e+00 4.71356e-12 2.00833e-09 * 3 3.241923769180135e-01 3.241923769098066e-01 8.20688e-12 3.45514e-10 ----------------------------------------------------------------------------------------------------- 30 C 1 3.282266697808694e+00 3.282266697752148e+00 5.65468e-11 1.12901e-10 2 -5.810755344810068e-01 -5.810755344923734e-01 1.13666e-11 5.22737e-11 3 -2.973366236360493e-01 -2.973366236352757e-01 7.73548e-13 1.82284e-10 ----------------------------------------------------------------------------------------------------- 31 C 1 2.882510236076805e-01 2.882510236022719e-01 5.40856e-12 6.55556e-10 2 9.345825319623818e-01 9.345825319613036e-01 1.07825e-12 4.48782e-09 * 3 -5.600674464499613e-01 -5.600674464603420e-01 1.03807e-11 3.90945e-10 ----------------------------------------------------------------------------------------------------- 32 C 1 2.131571636198211e+00 2.131571636210841e+00 1.26308e-11 3.88594e-10 2 1.926003583747859e+00 1.926003583730776e+00 1.70832e-11 9.28631e-11 3 4.858901820718145e+00 4.858901820720079e+00 1.93356e-12 4.82886e-10 ----------------------------------------------------------------------------------------------------- * 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 = 9, direction = 3: |F_model - F_numer| error = ----------------------- = 7.03898e-10 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 = Pd | a0 = 3.889087677001953 | (Configuration in file "config-Pd.xyz") ----------------------------------------------------------------------------------------------------- Part Spec Dir Force_model Force_numer |Force diff| uncertainty ----------------------------------------------------------------------------------------------------- 1 Pd 1 1.020375888998226e+00 1.020375888991899e+00 6.32694e-12 1.08361e-10 2 7.112934708439705e-01 7.112934708455544e-01 1.58384e-12 5.37954e-11 3 9.536810698091642e-01 9.536810698181069e-01 8.94262e-12 7.83809e-11 ----------------------------------------------------------------------------------------------------- 2 Pd 1 -5.967062843059638e-01 -5.967062842806415e-01 2.53223e-11 4.95274e-11 2 -3.488388354551437e-01 -3.488388354688122e-01 1.36685e-11 1.46655e-11 3 8.224918960537896e-01 8.224918960251331e-01 2.86565e-11 1.43128e-10 ----------------------------------------------------------------------------------------------------- 3 Pd 1 -1.205742646838988e-01 -1.205742646659828e-01 1.79160e-11 1.59709e-10 2 9.852512107081168e-01 9.852512106815842e-01 2.65327e-11 2.48266e-10 3 7.443142490103483e-02 7.443142490219215e-02 1.15732e-12 9.85952e-11 ----------------------------------------------------------------------------------------------------- 4 Pd 1 7.630213788353342e-01 7.630213788247686e-01 1.05655e-11 8.67815e-11 2 -7.249470567748914e-01 -7.249470568330806e-01 5.81892e-11 1.84287e-10 3 8.900213743050112e-03 8.900213786946972e-03 4.38969e-11 1.60317e-10 ----------------------------------------------------------------------------------------------------- 5 Pd 1 -4.986705323313650e+00 -4.986705323326200e+00 1.25500e-11 1.25296e-10 2 2.788432954946788e-01 2.788432955151217e-01 2.04429e-11 3.63153e-11 3 -3.968656465444207e+00 -3.968656465448229e+00 4.02167e-12 1.30797e-10 ----------------------------------------------------------------------------------------------------- 6 Pd 1 2.085534726059986e+00 2.085534726081772e+00 2.17857e-11 3.38313e-10 2 -3.300088560581102e+00 -3.300088560624029e+00 4.29270e-11 2.05171e-10 3 2.351933619211747e-01 2.351933619181411e-01 3.03360e-12 1.32219e-09 * ----------------------------------------------------------------------------------------------------- 7 Pd 1 4.551985321063414e+00 4.551985321035839e+00 2.75753e-11 5.10453e-09 * 2 1.093070092615015e+00 1.093070092599573e+00 1.54416e-11 5.79292e-11 3 5.139192883320702e+00 5.139192883287097e+00 3.36051e-11 8.22303e-09 * ----------------------------------------------------------------------------------------------------- 8 Pd 1 1.576359332709979e+00 1.576359332709232e+00 7.47180e-13 4.67295e-11 2 -3.156533513842416e+00 -3.156533513833472e+00 8.94351e-12 7.01399e-10 3 -3.810153271313307e+00 -3.810153271337054e+00 2.37477e-11 1.27059e-10 ----------------------------------------------------------------------------------------------------- 9 Pd 1 1.227891356680826e+00 1.227891356687320e+00 6.49392e-12 1.15126e-10 2 3.099199124091758e-01 3.099199124046010e-01 4.57478e-12 1.66451e-11 3 8.279607753745756e-01 8.279607753952238e-01 2.06483e-11 5.71371e-11 ----------------------------------------------------------------------------------------------------- 10 Pd 1 3.117320995277617e-03 3.117320996038510e-03 7.60893e-13 6.31709e-10 2 -1.601727605995171e+00 -1.601727606001273e+00 6.10179e-12 1.09205e-10 3 8.104665674079830e-01 8.104665674090384e-01 1.05549e-12 8.58492e-11 ----------------------------------------------------------------------------------------------------- 11 Pd 1 -2.852101169312235e+00 -2.852101169233038e+00 7.91962e-11 4.55244e-10 2 3.719052197847839e-03 3.719052181860556e-03 1.59873e-11 5.31242e-11 3 -7.657681639836671e-01 -7.657681639604261e-01 2.32410e-11 8.63875e-10 ----------------------------------------------------------------------------------------------------- 12 Pd 1 5.854186361200180e-01 5.854186361289265e-01 8.90854e-12 1.50658e-10 2 -7.765792660932744e-01 -7.765792660882823e-01 4.99212e-12 1.45698e-10 3 -2.693482580794651e-01 -2.693482580460266e-01 3.34385e-11 7.11575e-11 ----------------------------------------------------------------------------------------------------- 13 Pd 1 -3.576517063605802e+00 -3.576517063646968e+00 4.11653e-11 9.78950e-11 2 2.968238693104468e+00 2.968238693144087e+00 3.96199e-11 2.43385e-10 3 -5.729614795120550e-01 -5.729614795215583e-01 9.50329e-12 1.47431e-08 * ----------------------------------------------------------------------------------------------------- 14 Pd 1 -1.292596914632911e+00 -1.292596914653400e+00 2.04887e-11 1.43556e-10 2 -3.653973171125955e-01 -3.653973170856095e-01 2.69861e-11 6.61070e-11 3 1.226524280207688e+00 1.226524280217862e+00 1.01741e-11 8.02258e-11 ----------------------------------------------------------------------------------------------------- 15 Pd 1 4.047688359958697e-01 4.047688360187323e-01 2.28626e-11 1.05545e-10 2 2.000624542206621e-01 2.000624542021237e-01 1.85384e-11 7.34318e-11 3 -1.137052142877192e+00 -1.137052142951906e+00 7.47142e-11 9.41090e-10 ----------------------------------------------------------------------------------------------------- 16 Pd 1 2.020698149688233e+00 2.020698149725879e+00 3.76459e-11 1.77518e-10 2 1.340036918528615e+00 1.340036918522234e+00 6.38067e-12 1.36408e-10 3 8.906140858712349e-01 8.906140858816683e-01 1.04334e-11 2.29635e-10 ----------------------------------------------------------------------------------------------------- 17 Pd 1 8.393182430972187e-01 8.393182430695353e-01 2.76834e-11 5.93728e-10 2 8.554380449356911e-01 8.554380449187134e-01 1.69776e-11 1.25901e-10 3 -1.297953729110115e-01 -1.297953729832956e-01 7.22842e-11 2.20367e-10 ----------------------------------------------------------------------------------------------------- 18 Pd 1 -3.036910649786084e+00 -3.036910649796325e+00 1.02407e-11 1.73951e-10 2 1.512388518414359e+00 1.512388518411281e+00 3.07820e-12 1.24886e-10 3 -1.271986798097760e+00 -1.271986798102942e+00 5.18208e-12 2.34913e-10 ----------------------------------------------------------------------------------------------------- 19 Pd 1 -2.975109405225704e+00 -2.975109405252352e+00 2.66476e-11 1.93424e-10 2 -3.176007761117679e+00 -3.176007761124520e+00 6.84075e-12 9.08136e-11 3 7.230009060618823e-01 7.230009060490398e-01 1.28425e-11 8.33615e-11 ----------------------------------------------------------------------------------------------------- 20 Pd 1 8.044262372945704e-01 8.044262372770169e-01 1.75535e-11 7.84279e-11 2 1.227558673947698e-01 1.227558673744776e-01 2.02921e-11 1.89401e-09 * 3 -4.801710727102744e-01 -4.801710727290657e-01 1.87912e-11 3.09955e-10 ----------------------------------------------------------------------------------------------------- 21 Pd 1 -1.799067170031781e+00 -1.799067170054890e+00 2.31088e-11 1.23379e-10 2 -3.297478764886163e+00 -3.297478764866103e+00 2.00604e-11 1.41734e-10 3 -1.871451644672190e+00 -1.871451644669791e+00 2.39875e-12 4.82431e-09 * ----------------------------------------------------------------------------------------------------- 22 Pd 1 4.398199673796460e+00 4.398199673844708e+00 4.82476e-11 4.88225e-10 2 2.331330809156130e+00 2.331330809148562e+00 7.56772e-12 2.74634e-10 3 -4.079592746247881e+00 -4.079592746231833e+00 1.60485e-11 2.20616e-10 ----------------------------------------------------------------------------------------------------- 23 Pd 1 -5.924567314750206e-01 -5.924567314845753e-01 9.55469e-12 1.69238e-10 2 -5.886046814148593e-01 -5.886046814063938e-01 8.46556e-12 1.05941e-10 3 1.076394893212113e+00 1.076394893222762e+00 1.06493e-11 5.31772e-10 ----------------------------------------------------------------------------------------------------- 24 Pd 1 1.580810885422407e+00 1.580810885425338e+00 2.93121e-12 1.99275e-10 2 4.607183160817904e+00 4.607183160815954e+00 1.94955e-12 2.49296e-10 3 5.116558059448624e+00 5.116558059461082e+00 1.24576e-11 3.16938e-09 * ----------------------------------------------------------------------------------------------------- 25 Pd 1 -2.671328477500001e-01 -2.671328477459799e-01 4.02023e-12 6.96151e-10 2 -6.389564301336171e-01 -6.389564300709674e-01 6.26498e-11 3.59532e-10 3 -1.536514020618650e+00 -1.536514020667336e+00 4.86866e-11 1.05012e-10 ----------------------------------------------------------------------------------------------------- 26 Pd 1 1.960300761346650e-01 1.960300761298220e-01 4.84299e-12 4.60023e-10 2 -1.871133766474777e-02 -1.871133764386649e-02 2.08813e-11 4.96552e-09 * 3 -5.456546601067552e-01 -5.456546600601336e-01 4.66216e-11 2.10764e-09 * ----------------------------------------------------------------------------------------------------- 27 Pd 1 -8.859478118988120e-01 -8.859478118689221e-01 2.98899e-11 5.84116e-10 2 -2.300757489732926e+00 -2.300757489737484e+00 4.55769e-12 1.20435e-08 * 3 3.959044757276017e-01 3.959044757111650e-01 1.64366e-11 1.21102e-09 * ----------------------------------------------------------------------------------------------------- 28 Pd 1 2.600790808636931e-01 2.600790808704067e-01 6.71363e-12 1.65075e-10 2 2.107409862023993e-01 2.107409862074433e-01 5.04399e-12 1.20844e-10 3 2.153066520328286e-01 2.153066520174278e-01 1.54008e-11 3.77781e-10 ----------------------------------------------------------------------------------------------------- 29 Pd 1 1.258936762733916e-01 1.258936762249190e-01 4.84726e-11 2.38285e-10 2 1.981434849532913e+00 1.981434849505330e+00 2.75830e-11 4.45626e-10 3 2.463135962604272e+00 2.463135962700908e+00 9.66356e-11 2.39426e-10 ----------------------------------------------------------------------------------------------------- 30 Pd 1 -6.652132146670524e-01 -6.652132146712907e-01 4.23828e-12 4.64414e-11 2 -2.210643215027425e-01 -2.210643215094554e-01 6.71288e-12 1.94780e-10 3 4.084401445975934e-01 4.084401446124191e-01 1.48257e-11 2.30885e-10 ----------------------------------------------------------------------------------------------------- 31 Pd 1 -7.828214404739937e-01 -7.828214404926952e-01 1.87015e-11 6.39876e-11 2 6.145848406469978e-03 6.145848399460733e-03 7.00924e-12 1.49870e-10 3 -4.062259965274502e-01 -4.062259965467359e-01 1.92856e-11 1.54127e-10 ----------------------------------------------------------------------------------------------------- 32 Pd 1 1.985931471133340e+00 1.985931471136229e+00 2.88924e-12 9.44200e-11 2 9.978397573241471e-01 9.978397572855346e-01 3.86124e-11 1.19062e-10 3 -5.428655591934440e-01 -5.428655592065843e-01 1.31403e-11 6.61092e-11 ----------------------------------------------------------------------------------------------------- * 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 = 4, direction = 3: |F_model - F_numer| error = ----------------------- = 4.93211e-09 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 Pd | a0 = 3.669341514229472 | (Configuration in file "config-CPd.xyz") ----------------------------------------------------------------------------------------------------- Part Spec Dir Force_model Force_numer |Force diff| uncertainty ----------------------------------------------------------------------------------------------------- 1 C 1 9.373187428342671e-01 9.373187428487821e-01 1.45151e-11 1.37167e-11 2 7.678307499248533e-01 7.678307499246219e-01 2.31482e-13 2.64781e-10 3 8.785146910647738e-01 8.785146910594370e-01 5.33673e-12 4.24623e-11 ----------------------------------------------------------------------------------------------------- 2 Pd 1 -6.375085343155209e-01 -6.375085343384476e-01 2.29267e-11 6.57794e-11 2 -2.895603105122999e+00 -2.895603105078491e+00 4.45080e-11 4.35096e-11 3 -2.896306280373283e-01 -2.896306280282462e-01 9.08201e-12 3.03438e-10 ----------------------------------------------------------------------------------------------------- 3 Pd 1 -3.580467719559798e+00 -3.580467719613282e+00 5.34843e-11 1.75616e-10 2 4.575463112956828e-01 4.575463112815157e-01 1.41671e-11 4.70481e-10 3 -6.224329093750716e-01 -6.224329093132770e-01 6.17946e-11 1.41589e-10 ----------------------------------------------------------------------------------------------------- 4 C 1 1.749030410910865e-01 1.749030411098378e-01 1.87513e-11 8.86178e-11 2 5.391824545475310e-01 5.391824545430810e-01 4.45000e-12 1.51588e-10 3 -2.625806212491734e+00 -2.625806212471280e+00 2.04547e-11 6.19140e-11 ----------------------------------------------------------------------------------------------------- 5 C 1 -2.320648499034568e+00 -2.320648499021470e+00 1.30980e-11 9.99911e-11 2 8.734846393300478e-01 8.734846393521097e-01 2.20619e-11 4.70570e-11 3 -7.474205486023479e-01 -7.474205485992564e-01 3.09142e-12 4.73072e-11 ----------------------------------------------------------------------------------------------------- 6 C 1 -2.654748052279718e-01 -2.654748052464408e-01 1.84691e-11 8.03161e-11 2 9.880340205752943e-01 9.880340205789563e-01 3.66196e-12 1.56192e-10 3 -4.787733617944536e-01 -4.787733617899939e-01 4.45965e-12 4.21553e-10 ----------------------------------------------------------------------------------------------------- 7 C 1 -1.037225566505761e+00 -1.037225566508327e+00 2.56595e-12 5.89730e-10 2 -8.431848378685197e-02 -8.431848378497643e-02 1.87554e-12 1.44465e-10 3 2.534194055215570e+00 2.534194055183595e+00 3.19744e-11 1.20353e-10 ----------------------------------------------------------------------------------------------------- 8 Pd 1 -7.232883252383856e+00 -7.232883252438785e+00 5.49285e-11 1.04263e-10 2 -1.012073556716679e+00 -1.012073556753644e+00 3.69649e-11 7.61171e-10 3 -7.523259712574864e+00 -7.523259712506521e+00 6.83436e-11 8.66258e-10 ----------------------------------------------------------------------------------------------------- 9 Pd 1 -3.435926402559185e+00 -3.435926402570413e+00 1.12284e-11 1.00999e-10 2 -7.882671555930212e+00 -7.882671555977965e+00 4.77529e-11 1.83072e-10 3 -4.364972286672235e+00 -4.364972286687319e+00 1.50839e-11 3.73963e-10 ----------------------------------------------------------------------------------------------------- 10 C 1 3.294353923618555e+00 3.294353923568953e+00 4.96017e-11 1.90163e-10 2 1.231116824492396e+00 1.231116824452368e+00 4.00282e-11 2.96168e-10 3 -6.681065486960731e-01 -6.681065486683532e-01 2.77199e-11 2.42302e-10 ----------------------------------------------------------------------------------------------------- 11 Pd 1 1.639269495944487e+00 1.639269495980463e+00 3.59752e-11 1.87720e-09 * 2 -2.550394377028975e+00 -2.550394376906711e+00 1.22264e-10 1.47574e-09 * 3 7.384938029567458e-01 7.384938030399318e-01 8.31860e-11 9.55001e-10 ----------------------------------------------------------------------------------------------------- 12 Pd 1 -3.309088295114469e+00 -3.309088295152222e+00 3.77525e-11 2.89051e-10 2 5.020968803484729e+00 5.020968803405004e+00 7.97256e-11 4.10093e-10 3 2.569239178026231e+00 2.569239178063305e+00 3.70739e-11 3.24034e-10 ----------------------------------------------------------------------------------------------------- 13 Pd 1 8.554311857465475e-01 8.554311857896544e-01 4.31069e-11 4.53860e-10 2 3.410764293914401e+00 3.410764293912546e+00 1.85452e-12 2.93696e-10 3 -2.173584662417472e+00 -2.173584662436209e+00 1.87370e-11 1.77240e-09 * ----------------------------------------------------------------------------------------------------- 14 Pd 1 -2.539564590955766e+00 -2.539564590860127e+00 9.56391e-11 1.44770e-10 2 -1.534723340914010e+00 -1.534723340847306e+00 6.67044e-11 3.13544e-10 3 -4.863336377365126e-01 -4.863336377070023e-01 2.95103e-11 6.32914e-11 ----------------------------------------------------------------------------------------------------- 15 C 1 3.428163337839647e+00 3.428163337893614e+00 5.39671e-11 1.26050e-10 2 -1.611850490740814e+00 -1.611850490705287e+00 3.55265e-11 1.83362e-10 3 4.589269303054151e-01 4.589269302895589e-01 1.58561e-11 1.34082e-10 ----------------------------------------------------------------------------------------------------- 16 C 1 -9.629302413539728e-01 -9.629302413498732e-01 4.09961e-12 6.34709e-10 2 1.604833774905992e+00 1.604833774903243e+00 2.74913e-12 1.15639e-10 3 1.927772754270588e+00 1.927772754330095e+00 5.95068e-11 1.23213e-09 ----------------------------------------------------------------------------------------------------- 17 C 1 1.608202329278602e+00 1.608202329310749e+00 3.21472e-11 1.99235e-10 2 -8.351622749978774e-01 -8.351622749638279e-01 3.40495e-11 8.36180e-11 3 1.525513476187226e+00 1.525513476155621e+00 3.16047e-11 3.16431e-11 ----------------------------------------------------------------------------------------------------- 18 C 1 5.087555228112421e-01 5.087555226645283e-01 1.46714e-10 5.18606e-10 2 7.211434322425329e-01 7.211434322376549e-01 4.87799e-12 9.10216e-11 3 -2.670392337615584e+00 -2.670392337520657e+00 9.49267e-11 8.14354e-10 ----------------------------------------------------------------------------------------------------- 19 Pd 1 -9.271080996701976e-01 -9.271080996619376e-01 8.25995e-12 1.21231e-10 2 1.688476800691088e-01 1.688476800226844e-01 4.64244e-11 3.47558e-10 3 -6.836061580708508e-01 -6.836061580380982e-01 3.27526e-11 9.44801e-10 ----------------------------------------------------------------------------------------------------- 20 Pd 1 9.982213259693088e-01 9.982213259738726e-01 4.56379e-12 3.99820e-10 2 -3.283489336150041e-01 -3.283489336391338e-01 2.41297e-11 2.32257e-10 3 -7.617325606642976e-01 -7.617325606467388e-01 1.75588e-11 2.86018e-10 ----------------------------------------------------------------------------------------------------- 21 Pd 1 3.859284303287591e+00 3.859284303237168e+00 5.04228e-11 4.01321e-10 2 3.758574149019930e-01 3.758574148903924e-01 1.16006e-11 1.26784e-10 3 9.604259518445475e-01 9.604259518288528e-01 1.56947e-11 2.49549e-09 * ----------------------------------------------------------------------------------------------------- 22 Pd 1 1.111277267125739e+01 1.111277267130869e+01 5.12994e-11 2.97461e-09 * 2 -3.379864704005255e+00 -3.379864704077494e+00 7.22395e-11 2.68351e-10 3 9.932666187314096e+00 9.932666187467461e+00 1.53365e-10 1.99916e-09 * ----------------------------------------------------------------------------------------------------- 23 C 1 -9.384786908139160e-01 -9.384786908073298e-01 6.58618e-12 6.65318e-11 2 3.669198591723847e-01 3.669198591729235e-01 5.38736e-13 3.03918e-10 3 -1.202003989106739e+00 -1.202003989105140e+00 1.59939e-12 2.65443e-11 ----------------------------------------------------------------------------------------------------- 24 C 1 1.422201775815447e-01 1.422201776221793e-01 4.06346e-11 3.17751e-11 2 -1.149653602157974e+00 -1.149653602285114e+00 1.27141e-10 3.04418e-10 3 -5.539009762271248e-01 -5.539009764336342e-01 2.06509e-10 7.77520e-10 ----------------------------------------------------------------------------------------------------- 25 C 1 -1.393029328080309e-01 -1.393029327235247e-01 8.45062e-11 3.43628e-10 2 1.514243218384046e+00 1.514243218350978e+00 3.30684e-11 9.45327e-10 3 -3.594362718328666e-01 -3.594362718017060e-01 3.11606e-11 2.93463e-10 ----------------------------------------------------------------------------------------------------- 26 Pd 1 -3.826309068205853e-01 -3.826309068046831e-01 1.59022e-11 5.68064e-11 2 1.689536617635533e+00 1.689536617712351e+00 7.68179e-11 7.48890e-11 3 4.670974431623135e+00 4.670974431781034e+00 1.57899e-10 1.28162e-09 ----------------------------------------------------------------------------------------------------- 27 C 1 -1.450604129993883e+00 -1.450604130073619e+00 7.97353e-11 2.94972e-10 2 1.476281872121258e+00 1.476281871993485e+00 1.27773e-10 3.67909e-10 3 7.572911123584021e-01 7.572911124356495e-01 7.72474e-11 3.79009e-10 ----------------------------------------------------------------------------------------------------- 28 Pd 1 2.047777649955543e+00 2.047777649970731e+00 1.51883e-11 2.16478e-10 2 -9.445638879952154e-01 -9.445638879971274e-01 1.91203e-12 1.74373e-10 3 -1.809877848267891e+00 -1.809877848266878e+00 1.01297e-12 1.30977e-10 ----------------------------------------------------------------------------------------------------- 29 Pd 1 -3.359001721596019e+00 -3.359001721619998e+00 2.39786e-11 1.67843e-09 * 2 -2.625101544800884e-01 -2.625101544860775e-01 5.98910e-12 2.14616e-10 3 -2.119597580845679e+00 -2.119597580932890e+00 8.72116e-11 1.09386e-10 ----------------------------------------------------------------------------------------------------- 30 Pd 1 -7.548615059604187e-02 -7.548615059637605e-02 3.34177e-13 2.05204e-09 * 2 4.595759823121642e-01 4.595759823234594e-01 1.12952e-11 3.04628e-10 3 2.544696670606569e-01 2.544696670997311e-01 3.90742e-11 3.67956e-10 ----------------------------------------------------------------------------------------------------- 31 C 1 3.896119327991126e+00 3.896119327965209e+00 2.59166e-11 6.27018e-11 2 -1.524023451171871e+00 -1.524023451143248e+00 2.86222e-11 5.22069e-11 3 2.741377942269195e-01 2.741377942273069e-01 3.87357e-13 5.56321e-11 ----------------------------------------------------------------------------------------------------- 32 C 1 -1.908462495897393e+00 -1.908462495910646e+00 1.32525e-11 5.74360e-10 2 4.329593969353871e+00 4.329593969370734e+00 1.68630e-11 1.02630e-10 3 2.658248198574818e+00 2.658248198578329e+00 3.51097e-12 3.64779e-10 ----------------------------------------------------------------------------------------------------- * 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 = 25, direction = 1: |F_model - F_numer| error = ----------------------- = 6.06636e-10 max{|F_numer|, eps} ===================================================================================================== Grade is based on maximum error across all systems (for a model supporting multiple species): Maximum error = 4.93211e-09 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: A Comment: The letter grade A was assigned because the normalized error in the computation was 4.93211e-09 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.