Model Extended KIM ID = === Verification check vc-forces-numerical-derivative start (2019-07-11 05:34:07) === !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!! !!!!! !!!!! 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 : EAM_Dynamo_GrocholaRusso_2005_Au__MO_557267801129_000 Supported species : Au 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 = Au | a0 = 4.070110511779786 | (Configuration in file "config-Au.xyz") ----------------------------------------------------------------------------------------------------- Part Spec Dir Force_model Force_numer |Force diff| uncertainty ----------------------------------------------------------------------------------------------------- 1 Au 1 1.010035478066550e+00 1.010035478248741e+00 1.82191e-10 7.08920e-10 2 5.853361883983091e-01 5.853361882965370e-01 1.01772e-10 2.70803e-10 3 6.107981919235184e-01 6.107981919576725e-01 3.41541e-11 5.17394e-10 ----------------------------------------------------------------------------------------------------- 2 Au 1 -2.383956862420664e-01 -2.383956863374488e-01 9.53824e-11 6.27464e-10 2 9.094702133205792e-01 9.094702135436059e-01 2.23027e-10 3.39009e-10 3 1.716628995296191e+00 1.716628995262520e+00 3.36711e-11 1.52029e-09 ----------------------------------------------------------------------------------------------------- 3 Au 1 5.626057979404463e-01 5.626057980836231e-01 1.43177e-10 1.06406e-09 2 1.344738905900944e+00 1.344738905892509e+00 8.43547e-12 1.06025e-09 3 -8.092711085067503e-02 -8.092711095745178e-02 1.06777e-10 8.02609e-10 ----------------------------------------------------------------------------------------------------- 4 Au 1 1.625755602230025e+00 1.625755602222720e+00 7.30527e-12 2.69406e-10 2 5.094516718434593e-01 5.094516718524831e-01 9.02378e-12 6.29183e-10 3 6.288548083390690e-01 6.288548085025952e-01 1.63526e-10 2.15924e-10 ----------------------------------------------------------------------------------------------------- 5 Au 1 -5.076846020969046e-02 -5.076846033237224e-02 1.22682e-10 5.58647e-10 2 8.205760812744242e-01 8.205760812623130e-01 1.21112e-11 1.18123e-09 3 1.234515462158848e+00 1.234515462148410e+00 1.04381e-11 2.89965e-10 ----------------------------------------------------------------------------------------------------- 6 Au 1 -1.225499176368516e+00 -1.225499176485351e+00 1.16835e-10 7.15474e-10 2 5.723484676715150e-01 5.723484679490419e-01 2.77527e-10 4.93248e-10 3 1.327661047318559e+00 1.327661047110428e+00 2.08130e-10 1.49758e-09 ----------------------------------------------------------------------------------------------------- 7 Au 1 2.155447791681604e-01 2.155447787401994e-01 4.27961e-10 8.79883e-10 2 -4.952476399015704e-01 -4.952476398285936e-01 7.29768e-11 1.20109e-09 3 -2.044810616377883e+00 -2.044810616372842e+00 5.04130e-12 1.74625e-09 ----------------------------------------------------------------------------------------------------- 8 Au 1 -1.360997042240869e+00 -1.360997042248439e+00 7.56972e-12 2.65459e-09 2 -8.269358432550863e-01 -8.269358432273595e-01 2.77268e-11 1.76341e-10 3 -4.559719973731516e-01 -4.559719971808194e-01 1.92332e-10 5.27112e-10 ----------------------------------------------------------------------------------------------------- 9 Au 1 7.736070008933479e-01 7.736070010410833e-01 1.47735e-10 2.39326e-10 2 -9.212004147999007e-01 -9.212004149410974e-01 1.41197e-10 6.21978e-10 3 8.192442742210895e-01 8.192442742997562e-01 7.86667e-11 3.38864e-10 ----------------------------------------------------------------------------------------------------- 10 Au 1 -2.667473567586988e-01 -2.667473567408332e-01 1.78655e-11 5.65171e-11 2 -6.539405650819328e-01 -6.539405650848392e-01 2.90645e-12 1.80358e-09 3 2.385256836434579e-01 2.385256837405480e-01 9.70901e-11 1.68802e-10 ----------------------------------------------------------------------------------------------------- 11 Au 1 -2.211928815928783e-01 -2.211928818471826e-01 2.54304e-10 7.66101e-10 2 -9.220721421407614e-01 -9.220721421595092e-01 1.87478e-11 1.00982e-09 3 1.495946255649577e-01 1.495946254371385e-01 1.27819e-10 1.27232e-09 ----------------------------------------------------------------------------------------------------- 12 Au 1 1.135388350018985e+00 1.135388349999098e+00 1.98872e-11 1.17538e-09 2 -5.854236208282324e-01 -5.854236207271111e-01 1.01121e-10 2.06581e-09 3 4.743069050612508e-01 4.743069049577546e-01 1.03496e-10 2.74246e-10 ----------------------------------------------------------------------------------------------------- 13 Au 1 -1.216579834709021e+00 -1.216579834734909e+00 2.58884e-11 4.74127e-10 2 -1.468570271236509e+00 -1.468570271189326e+00 4.71827e-11 1.18789e-09 3 -1.350132891955259e+00 -1.350132891877289e+00 7.79701e-11 2.41628e-09 ----------------------------------------------------------------------------------------------------- 14 Au 1 -7.462804711486671e-01 -7.462804711948027e-01 4.61355e-11 5.18783e-10 2 -4.726866419386256e-01 -4.726866419893917e-01 5.07661e-11 1.19956e-10 3 7.926520298775220e-01 7.926520298943868e-01 1.68647e-11 5.99376e-10 ----------------------------------------------------------------------------------------------------- 15 Au 1 1.403046488578206e+00 1.403046488890122e+00 3.11916e-10 4.88249e-10 2 -9.459631952691474e-01 -9.459631952951477e-01 2.60003e-11 2.88915e-10 3 7.121348607615837e-01 7.121348608358827e-01 7.42990e-11 6.17613e-10 ----------------------------------------------------------------------------------------------------- 16 Au 1 5.377367642948258e-01 5.377367642203033e-01 7.45225e-11 4.17045e-10 2 8.817543014813919e-01 8.817543015193378e-01 3.79459e-11 5.75627e-10 3 1.162365586659678e+00 1.162365586735229e+00 7.55511e-11 7.88702e-10 ----------------------------------------------------------------------------------------------------- 17 Au 1 6.672251735390408e-01 6.672251735610630e-01 2.20223e-11 1.73670e-09 2 9.238164716991275e-01 9.238164717842245e-01 8.50970e-11 3.72080e-10 3 -5.975927191570196e-01 -5.975927191757228e-01 1.87032e-11 9.06993e-10 ----------------------------------------------------------------------------------------------------- 18 Au 1 2.529770496398832e-01 2.529770497323426e-01 9.24595e-11 4.82115e-10 2 8.397001727971976e-01 8.397001727970264e-01 1.71196e-13 1.67995e-10 3 -2.762817055654834e-01 -2.762817055025472e-01 6.29362e-11 2.40032e-10 ----------------------------------------------------------------------------------------------------- 19 Au 1 6.184402971049444e-01 6.184402970650520e-01 3.98924e-11 8.50860e-10 2 7.308316174711932e-01 7.308316174437678e-01 2.74254e-11 1.01819e-09 3 -7.685939271212660e-01 -7.685939272344184e-01 1.13152e-10 6.93155e-10 ----------------------------------------------------------------------------------------------------- 20 Au 1 1.168730746107295e+00 1.168730746144483e+00 3.71880e-11 6.28672e-10 2 7.536074746287948e-01 7.536074745786541e-01 5.01408e-11 4.32240e-10 3 -1.097667753709363e+00 -1.097667753614884e+00 9.44786e-11 4.08801e-10 ----------------------------------------------------------------------------------------------------- 21 Au 1 -1.757804399257012e+00 -1.757804399446472e+00 1.89460e-10 6.17770e-10 2 1.067654784255358e+00 1.067654784518768e+00 2.63410e-10 4.29048e-10 3 1.118596046100145e+00 1.118596045967923e+00 1.32222e-10 3.08258e-09 ----------------------------------------------------------------------------------------------------- 22 Au 1 -8.658581124160958e-01 -8.658581123302337e-01 8.58621e-11 1.06520e-09 2 1.984025298969458e+00 1.984025298869371e+00 1.00087e-10 3.47810e-10 3 1.133196896985482e+00 1.133196897078827e+00 9.33451e-11 3.65869e-10 ----------------------------------------------------------------------------------------------------- 23 Au 1 -9.502956571978515e-01 -9.502956572672546e-01 6.94031e-11 1.86653e-09 2 7.475412359517735e-01 7.475412359440046e-01 7.76890e-12 1.76112e-09 3 -9.035015051565125e-01 -9.035015049306814e-01 2.25831e-10 1.12566e-09 ----------------------------------------------------------------------------------------------------- 24 Au 1 -9.848324567165800e-02 -9.848324561580218e-02 5.58558e-11 2.73079e-10 2 4.329536426733211e-01 4.329536421225023e-01 5.50819e-10 5.27575e-10 3 -1.734494974684470e+00 -1.734494974673203e+00 1.12670e-11 6.68314e-11 ----------------------------------------------------------------------------------------------------- 25 Au 1 6.430549221181275e-01 6.430549222380004e-01 1.19873e-10 1.89894e-09 2 -8.615521012180714e-01 -8.615521012481495e-01 3.00782e-11 1.11579e-10 3 4.318961943707738e-01 4.318961945076730e-01 1.36899e-10 2.38921e-10 ----------------------------------------------------------------------------------------------------- 26 Au 1 -4.093165603812706e-01 -4.093165603789321e-01 2.33852e-12 2.60368e-10 2 -5.284638687589702e-01 -5.284638686614483e-01 9.75219e-11 1.59949e-09 3 -9.213915279917609e-01 -9.213915281993488e-01 2.07588e-10 5.54132e-10 ----------------------------------------------------------------------------------------------------- 27 Au 1 1.126209486583526e+00 1.126209486539931e+00 4.35945e-11 3.05102e-10 2 -2.351754513885423e+00 -2.351754513888300e+00 2.87637e-12 7.07763e-10 3 -9.764270468245770e-01 -9.764270471989360e-01 3.74359e-10 1.03522e-09 ----------------------------------------------------------------------------------------------------- 28 Au 1 -2.866752683529629e-01 -2.866752683591847e-01 6.22180e-12 1.26997e-09 2 2.764284101967420e-01 2.764284101703677e-01 2.63743e-11 1.39009e-09 3 -8.251490452189295e-01 -8.251490450822185e-01 1.36711e-10 2.33843e-10 ----------------------------------------------------------------------------------------------------- 29 Au 1 -6.829482978400141e-01 -6.829482979587291e-01 1.18715e-10 1.55472e-09 2 -6.889960903052490e-01 -6.889960902810636e-01 2.41854e-11 3.37146e-10 3 6.726329488318566e-01 6.726329489398188e-01 1.07962e-10 4.25634e-10 ----------------------------------------------------------------------------------------------------- 30 Au 1 -1.095252901987213e+00 -1.095252901947436e+00 3.97771e-11 8.11753e-10 2 -4.267843129823182e-01 -4.267843130043267e-01 2.20085e-11 1.69930e-09 3 2.732590227618633e-01 2.732590227472980e-01 1.45653e-11 8.76256e-10 ----------------------------------------------------------------------------------------------------- 31 Au 1 -4.345704702079245e-01 -4.345704701930962e-01 1.48282e-11 6.20054e-10 2 -8.900646941435937e-01 -8.900646941665912e-01 2.29975e-11 3.50591e-10 3 -7.439529887507017e-01 -7.439529889328775e-01 1.82176e-10 1.00728e-09 ----------------------------------------------------------------------------------------------------- 32 Au 1 1.673078862990472e-01 1.673078862959732e-01 3.07401e-12 8.35373e-10 2 -3.405790227881978e-01 -3.405790229141918e-01 1.25994e-10 6.57984e-10 3 -7.199677691387943e-01 -7.199677692267612e-01 8.79669e-11 1.86054e-10 ----------------------------------------------------------------------------------------------------- Maximum error obtained for particle = 5, direction = 1: |F_model - F_numer| error = ----------------------- = 2.41650e-09 max{|F_numer|, eps} ===================================================================================================== Grade is based on maximum error across all systems (for a model supporting multiple species): Maximum error = 2.41650e-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 2.41650e-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. === Verification check vc-forces-numerical-derivative end (2019-07-11 05:34:11) ===