!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!! !!!!! !!!!! 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_000 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.343056362663442e+00 2.343056337190397e+00 2.54730e-08 7.95140e-07 2 2.286267798854501e+00 2.286267773195394e+00 2.56591e-08 1.08535e-07 3 4.177354173269533e+00 4.177353922866284e+00 2.50403e-07 5.54970e-07 ----------------------------------------------------------------------------------------------------- 2 C 1 3.424308003052464e+00 3.424307982443882e+00 2.06086e-08 2.28410e-07 2 -4.139595757221547e+00 -4.139595720925910e+00 3.62956e-08 3.40931e-04 * 3 -7.219626138056181e+00 -7.219626317813777e+00 1.79758e-07 7.78907e-07 ----------------------------------------------------------------------------------------------------- 3 C 1 -1.415220043305832e+00 -1.415219496262032e+00 5.47044e-07 9.09363e-06 * 2 -1.403922612627726e+00 -1.403922601784559e+00 1.08432e-08 9.67217e-07 3 -1.388728058794796e+00 -1.388728361125487e+00 3.02331e-07 7.32084e-06 ----------------------------------------------------------------------------------------------------- 4 C 1 -2.453647595739028e+00 -2.453647678074436e+00 8.23354e-08 3.98056e-07 2 -2.675607557352134e+00 -2.675607541326635e+00 1.60255e-08 4.67544e-07 3 -5.086319600557945e+00 -5.086319649157522e+00 4.85996e-08 4.68387e-07 ----------------------------------------------------------------------------------------------------- 5 C 1 -1.401956018394078e+00 -1.401956731845113e+00 7.13451e-07 1.85211e-06 2 4.040372923210039e+00 4.040373184622105e+00 2.61412e-07 2.30127e-04 * 3 3.942728961204490e+00 3.942728963288881e+00 2.08439e-09 1.65352e-07 ----------------------------------------------------------------------------------------------------- 6 C 1 -2.982542765641677e+00 -2.982542817163320e+00 5.15216e-08 4.93947e-06 2 -1.083691928446763e+00 -1.083691956908687e+00 2.84619e-08 3.69466e-07 3 1.130075716408048e+00 1.130075739902235e+00 2.34942e-08 6.36929e-07 ----------------------------------------------------------------------------------------------------- 7 C 1 -4.093824667826520e+00 -4.093824680759282e+00 1.29328e-08 6.92566e-08 2 1.015628433082176e+01 1.015628422051226e+01 1.10309e-07 1.92560e-06 3 5.334086073777645e+00 5.334086091401645e+00 1.76240e-08 7.31844e-07 ----------------------------------------------------------------------------------------------------- 8 C 1 1.146023828595281e+01 1.146023849724112e+01 2.11288e-07 1.49848e-06 2 6.098853077863545e+00 6.098852946478503e+00 1.31385e-07 3.40196e-04 * 3 1.146275209386697e+00 1.146275183508195e+00 2.58785e-08 8.90901e-07 ----------------------------------------------------------------------------------------------------- 9 C 1 -6.531327898699722e+00 -6.531327773049600e+00 1.25650e-07 1.01192e-06 2 2.916706394473638e+00 2.916706421114176e+00 2.66405e-08 1.89126e-07 3 2.601442708801232e-01 2.601441442973816e-01 1.26583e-07 2.67133e-06 ----------------------------------------------------------------------------------------------------- 10 C 1 -8.549910640932252e-01 -8.549911449540682e-01 8.08608e-08 6.05777e-07 2 5.250389043234924e+00 5.250389215843451e+00 1.72609e-07 3.12096e-06 3 -1.557916596350334e+00 -1.557916553710988e+00 4.26393e-08 4.23073e-06 ----------------------------------------------------------------------------------------------------- 11 C 1 2.350834719421043e+01 2.350834654705520e+01 6.47155e-07 1.78604e-05 * 2 -9.808735991293626e-01 -9.808736275496754e-01 2.84203e-08 5.81005e-07 3 3.565672316270442e+01 3.565672265424475e+01 5.08460e-07 3.26212e-06 ----------------------------------------------------------------------------------------------------- 12 C 1 -8.294770763016395e-01 -8.294774698330513e-01 3.93531e-07 3.45928e-04 * 2 -2.128836890347824e+00 -2.128837032325213e+00 1.41977e-07 2.77335e-06 3 -6.096823541543060e+00 -6.096823551983316e+00 1.04403e-08 1.22037e-06 ----------------------------------------------------------------------------------------------------- 13 C 1 8.911081094809019e+00 8.911080822122841e+00 2.72686e-07 4.42868e-06 2 -8.187523146504237e+00 -8.187523082804734e+00 6.36995e-08 1.92705e-06 3 2.113852137142314e+00 2.113852364752632e+00 2.27610e-07 3.37712e-04 * ----------------------------------------------------------------------------------------------------- 14 C 1 1.940999182282095e+00 1.940999163141523e+00 1.91406e-08 4.18580e-07 2 -4.088087054401336e+00 -4.088087043562681e+00 1.08387e-08 6.51485e-07 3 1.179053593198213e+00 1.179053311832355e+00 2.81366e-07 1.36793e-06 ----------------------------------------------------------------------------------------------------- 15 C 1 -8.230372570183270e+00 -8.230372700429795e+00 1.30247e-07 2.61876e-06 2 9.789304149101394e+00 9.789303991267760e+00 1.57834e-07 2.40233e-07 3 -2.742445726285029e+00 -2.742445809936467e+00 8.36514e-08 8.00427e-07 ----------------------------------------------------------------------------------------------------- 16 C 1 1.048852264264293e+01 1.048852299541390e+01 3.52771e-07 1.52096e-05 * 2 8.622015367930079e+00 8.622016391293144e+00 1.02336e-06 2.94886e-06 3 7.739888755041107e+00 7.739888629800053e+00 1.25241e-07 1.43189e-06 ----------------------------------------------------------------------------------------------------- 17 C 1 3.161538769770973e+00 3.161538715864191e+00 5.39068e-08 8.05167e-07 2 1.824007842149314e+00 1.824007910412342e+00 6.82630e-08 1.46375e-06 3 1.568035913365545e+00 1.568035858948273e+00 5.44173e-08 3.05413e-07 ----------------------------------------------------------------------------------------------------- 18 C 1 -6.878085313109154e+00 -6.878085235409199e+00 7.77000e-08 2.34571e-07 2 -7.851120369194744e+00 -7.851120213227263e+00 1.55967e-07 4.97297e-07 3 8.331863471910705e+00 8.331863593705812e+00 1.21795e-07 1.74705e-07 ----------------------------------------------------------------------------------------------------- 19 C 1 -2.428296840974531e-01 -2.428296567204108e-01 2.73770e-08 3.22058e-07 2 4.399702891864878e+00 4.399702859277501e+00 3.25874e-08 1.84970e-07 3 -1.609144329869874e-01 -1.609144649789743e-01 3.19920e-08 2.25144e-07 ----------------------------------------------------------------------------------------------------- 20 C 1 3.875188839816483e+00 3.875188862658297e+00 2.28418e-08 9.89407e-07 2 -5.961293511681837e-02 -5.961300530890759e-02 7.01921e-08 6.89740e-07 3 -4.239476934572235e+00 -4.239476898645231e+00 3.59270e-08 6.40442e-07 ----------------------------------------------------------------------------------------------------- 21 C 1 4.586101019964572e+00 4.586100914940447e+00 1.05024e-07 9.28303e-07 2 -2.913727396771900e+00 -2.913727433962393e+00 3.71905e-08 6.02060e-07 3 4.038218058508810e-01 4.038214585699023e-01 3.47281e-07 8.07894e-07 ----------------------------------------------------------------------------------------------------- 22 C 1 7.384444507913124e-01 7.384444858046876e-01 3.50134e-08 9.09679e-07 2 -9.084587997007981e+00 -9.084588088820372e+00 9.18124e-08 1.99039e-06 3 -3.397562579405941e+00 -3.397562594171784e+00 1.47658e-08 2.06306e-07 ----------------------------------------------------------------------------------------------------- 23 C 1 -4.368634711459892e+00 -4.368634723740997e+00 1.22811e-08 5.25427e-08 2 1.563311498903365e+00 1.563311507002038e+00 8.09867e-09 2.52728e-07 3 -4.021875800521383e+00 -4.021875756553165e+00 4.39682e-08 3.41811e-07 ----------------------------------------------------------------------------------------------------- 24 C 1 -1.582688370942182e+00 -1.582688328603989e+00 4.23382e-08 2.56452e-07 2 -3.480503809617153e+00 -3.480503829531432e+00 1.99143e-08 2.44145e-06 3 -4.153790677230790e-01 -4.153791063544072e-01 3.86313e-08 8.57169e-07 ----------------------------------------------------------------------------------------------------- 25 C 1 2.493452333899941e+00 2.493452457823155e+00 1.23923e-07 3.47058e-04 * 2 1.859296184484432e+00 1.859296157517106e+00 2.69673e-08 2.88029e-07 3 2.470220723959586e+00 2.470220391046498e+00 3.32913e-07 3.34809e-04 * ----------------------------------------------------------------------------------------------------- 26 C 1 3.835039970147344e+00 3.835040098421678e+00 1.28274e-07 2.20604e-06 2 -3.716029830779966e+00 -3.716029464052496e+00 3.66727e-07 1.19564e-06 3 7.971450351131651e-01 7.971450653499905e-01 3.02368e-08 3.36141e-04 * ----------------------------------------------------------------------------------------------------- 27 C 1 -9.347076951824118e+00 -9.347076989949672e+00 3.81256e-08 6.04620e-07 2 9.687092994909252e+00 9.687092965457847e+00 2.94514e-08 1.83145e-06 3 4.836781435917183e-01 4.836782049252014e-01 6.13335e-08 2.47551e-07 ----------------------------------------------------------------------------------------------------- 28 C 1 8.914422219072291e+00 8.914422244152640e+00 2.50803e-08 4.86663e-07 2 -1.026505906725354e+01 -1.026505905801764e+01 9.23590e-09 2.41529e-07 3 1.277855967795142e+00 1.277855952269456e+00 1.55257e-08 7.12758e-08 ----------------------------------------------------------------------------------------------------- 29 C 1 -1.978625010281759e+01 -1.978624950940793e+01 5.93410e-07 2.64439e-05 * 2 1.841846684623683e+00 1.841846478433255e+00 2.06190e-07 2.04875e-06 3 -2.897409146077961e+01 -2.897409107950966e+01 3.81270e-07 1.17203e-06 ----------------------------------------------------------------------------------------------------- 30 C 1 -1.130220682489521e+01 -1.130220789779415e+01 1.07290e-06 2.29275e-06 2 -6.655751767141250e+00 -6.655751759014384e+00 8.12687e-09 8.21791e-07 3 -8.767191633317886e+00 -8.767191626328023e+00 6.98986e-09 2.20238e-06 ----------------------------------------------------------------------------------------------------- 31 C 1 -6.229868289635354e+00 -6.229868618475920e+00 3.28841e-07 4.61754e-07 2 3.110573282522463e+00 3.110573287638042e+00 5.11558e-09 9.23800e-07 3 -3.036832063935527e+00 -3.036832080199926e+00 1.62644e-08 2.20529e-07 ----------------------------------------------------------------------------------------------------- 32 C 1 -1.149740420110154e+00 -1.149740433580821e+00 1.34707e-08 3.68546e-07 2 -4.731492746032979e+00 -4.731492766199908e+00 2.01669e-08 2.64136e-07 3 -9.076194797693131e-01 -9.076193911210765e-01 8.86482e-08 7.52499e-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.17746e-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.034253901744184e+00 -1.034253918764757e+00 1.70206e-08 2.05363e-07 2 -8.450471017701440e-01 -8.450470764011430e-01 2.53690e-08 2.37268e-07 3 -1.526888175491671e+00 -1.526888203900924e+00 2.84093e-08 2.57734e-08 ----------------------------------------------------------------------------------------------------- 2 H 1 1.770155238871000e+00 1.770153433410875e+00 1.80546e-06 9.19409e-07 2 1.116376206478238e-01 1.116370991775719e-01 5.21470e-07 7.38265e-06 3 -4.514185305463033e+00 -4.514184930410449e+00 3.75053e-07 8.07131e-06 ----------------------------------------------------------------------------------------------------- 3 H 1 -1.232327618716385e+00 -1.232328937325567e+00 1.31861e-06 4.84840e-06 2 -1.126050007219055e+00 -1.126049678516616e+00 3.28702e-07 5.77239e-06 3 1.939141250195860e+00 1.939141222905282e+00 2.72906e-08 3.19102e-06 ----------------------------------------------------------------------------------------------------- 4 H 1 -3.547686161919594e+00 -3.547686394901039e+00 2.32981e-07 5.28409e-07 2 -2.679702559411613e+00 -2.679702402732070e+00 1.56680e-07 3.65780e-05 * 3 -2.336909486185431e+00 -2.336909447964012e+00 3.82214e-08 1.94443e-06 ----------------------------------------------------------------------------------------------------- 5 H 1 3.143078968607829e-01 3.143078071954592e-01 8.96653e-08 4.58901e-07 2 5.671524535073481e-01 5.671524533504622e-01 1.56886e-10 2.70628e-07 3 -6.861646269701227e-01 -6.861647035969693e-01 7.66268e-08 6.34697e-07 ----------------------------------------------------------------------------------------------------- 6 H 1 3.489319586564747e-01 3.489305007303304e-01 1.45793e-06 2.87663e-06 2 -1.473545577035512e+00 -1.473545684972997e+00 1.07937e-07 2.82948e-06 3 5.333382672329842e-01 5.333368517158348e-01 1.41552e-06 1.95698e-05 ----------------------------------------------------------------------------------------------------- 7 H 1 -1.050716077406576e-01 -1.050715807532935e-01 2.69874e-08 7.26737e-07 2 3.791158331478283e-01 3.791158217117799e-01 1.14360e-08 3.65086e-08 3 -3.303899680648130e-01 -3.303899937125658e-01 2.56478e-08 1.03701e-06 ----------------------------------------------------------------------------------------------------- 8 H 1 3.120303440205182e+00 3.120273869430605e+00 2.95708e-05 7.49087e-05 * 2 -3.561577269829616e+00 -3.561574955556080e+00 2.31427e-06 3.25735e-05 * 3 5.821398296254291e-01 5.821453620810572e-01 5.53246e-06 5.87931e-05 * ----------------------------------------------------------------------------------------------------- 9 H 1 1.641148275119400e-01 1.641148788702565e-01 5.13583e-08 2.60492e-07 2 8.972960658513720e-01 8.972959360284406e-01 1.29823e-07 2.39780e-06 3 -3.056878567338279e-01 -3.056879429030951e-01 8.61693e-08 7.35477e-07 ----------------------------------------------------------------------------------------------------- 10 H 1 -3.252554980818667e+00 -3.252555134792765e+00 1.53974e-07 3.14190e-07 2 2.341220736527118e+00 2.341220707314392e+00 2.92127e-08 3.22002e-06 3 -1.333032226278943e+00 -1.333032284685027e+00 5.84061e-08 1.27181e-06 ----------------------------------------------------------------------------------------------------- 11 H 1 -2.195001384143433e+00 -2.195001385404183e+00 1.26075e-09 2.37769e-06 2 2.600374901475929e+00 2.600373756194096e+00 1.14528e-06 2.25387e-06 3 -1.819012526257891e-01 -1.819012546487197e-01 2.02293e-09 2.61046e-06 ----------------------------------------------------------------------------------------------------- 12 H 1 1.505114612878601e+00 1.505114424328176e+00 1.88550e-07 5.63458e-07 2 -1.909928761281701e-01 -1.909927683005272e-01 1.07828e-07 1.78293e-06 3 -1.311797919968250e+00 -1.311798066960263e+00 1.46992e-07 2.28070e-06 ----------------------------------------------------------------------------------------------------- 13 H 1 4.368079954720031e-01 4.368036928520110e-01 4.30262e-06 6.43765e-05 * 2 3.328883054382310e+00 3.328883997335461e+00 9.42953e-07 8.68259e-05 * 3 -4.507064614687441e+00 -4.507061583495875e+00 3.03119e-06 6.39986e-06 ----------------------------------------------------------------------------------------------------- 14 H 1 5.849353688356461e-01 5.849352791745497e-01 8.96611e-08 3.76813e-07 2 9.674849053674327e-01 9.674850070718763e-01 1.01704e-07 3.45104e-06 3 -2.944571432922189e-03 -2.944802949863120e-03 2.31517e-07 5.99866e-07 ----------------------------------------------------------------------------------------------------- 15 H 1 3.798197397610771e+00 3.798196789272909e+00 6.08338e-07 1.77467e-05 2 -3.025116980140427e-01 -3.025155966754021e-01 3.89866e-06 8.63105e-05 * 3 1.785053429492054e+00 1.785053144647065e+00 2.84845e-07 6.99873e-06 ----------------------------------------------------------------------------------------------------- 16 H 1 -2.764310511700536e-01 -2.764292098623882e-01 1.84131e-06 5.27131e-06 2 5.237696753348182e-01 5.237688928729020e-01 7.82462e-07 5.78458e-06 3 -1.955586162124481e+00 -1.955586309667648e+00 1.47543e-07 3.74000e-06 ----------------------------------------------------------------------------------------------------- 17 H 1 -4.566018211500398e-01 -4.566018633961320e-01 4.22461e-08 4.21615e-06 2 -2.250342959169774e+00 -2.250342918131792e+00 4.10380e-08 1.98789e-06 3 3.734949940161069e+00 3.734949889167452e+00 5.09936e-08 1.09386e-07 ----------------------------------------------------------------------------------------------------- 18 H 1 -4.600910420079130e+00 -4.600909296312942e+00 1.12377e-06 1.13060e-05 2 1.940929053742325e+00 1.940930070287140e+00 1.01654e-06 2.10378e-05 3 7.315529798148755e+00 7.315531304798479e+00 1.50665e-06 1.06026e-05 ----------------------------------------------------------------------------------------------------- 19 H 1 -2.126032920974399e+00 -2.126033391566448e+00 4.70592e-07 7.09406e-07 2 1.478435018606335e+00 1.478434961710533e+00 5.68958e-08 6.38390e-07 3 -4.365804898798168e-01 -4.365805217212587e-01 3.18414e-08 2.15503e-07 ----------------------------------------------------------------------------------------------------- 20 H 1 7.780983447306847e-01 7.780983601196393e-01 1.53890e-08 2.31537e-07 2 -1.852252024411228e+00 -1.852251926801702e+00 9.76095e-08 1.32056e-06 3 -1.954436217496991e-01 -1.954436460036874e-01 2.42540e-08 2.43394e-07 ----------------------------------------------------------------------------------------------------- 21 H 1 2.218026923169131e+00 2.218027370438914e+00 4.47270e-07 5.69595e-06 2 1.491174723631093e+00 1.491175609935057e+00 8.86304e-07 1.06183e-05 3 1.116835851862586e+00 1.116835237021135e+00 6.14841e-07 2.14344e-06 ----------------------------------------------------------------------------------------------------- 22 H 1 -7.043713343554905e-01 -7.043784426234657e-01 7.10827e-06 1.01979e-04 * 2 -1.587114196548905e+00 -1.587114134735699e+00 6.18132e-08 1.76678e-06 3 1.111097543537409e+00 1.111086945622207e+00 1.05979e-05 1.60959e-04 * ----------------------------------------------------------------------------------------------------- 23 H 1 -6.373213783629770e-01 -6.373213705775690e-01 7.78541e-09 1.87045e-07 2 1.950650525534144e-01 1.950649474673687e-01 1.05086e-07 1.59668e-06 3 -7.419719814546499e-01 -7.419720282782795e-01 4.68236e-08 1.83812e-06 ----------------------------------------------------------------------------------------------------- 24 H 1 -1.216794163641506e+00 -1.216794195223585e+00 3.15821e-08 4.13799e-07 2 9.286662585938654e-01 9.286663311988808e-01 7.26050e-08 2.65261e-07 3 -3.736436809624131e-01 -3.736437910449293e-01 1.10083e-07 8.11654e-07 ----------------------------------------------------------------------------------------------------- 25 H 1 3.764064146183354e+00 3.764064262415047e+00 1.16232e-07 3.07402e-06 2 -3.074091096748224e+00 -3.074091240442077e+00 1.43694e-07 9.21158e-07 3 -4.621789741685296e-01 -4.621786586220369e-01 3.15546e-07 1.83716e-06 ----------------------------------------------------------------------------------------------------- 26 H 1 -8.182228191548001e-02 -8.182241424215261e-02 1.32327e-07 1.54300e-06 2 1.355695737907616e+00 1.355695748055114e+00 1.01475e-08 1.04146e-06 3 6.070629205063400e-01 6.070629178139930e-01 2.69235e-09 1.82634e-06 ----------------------------------------------------------------------------------------------------- 27 H 1 4.617274176895738e-01 4.617274413427879e-01 2.36532e-08 1.77003e-06 2 1.207125852828348e-01 1.207126175814200e-01 3.22986e-08 9.21758e-08 3 1.103719404005444e+00 1.103719413871975e+00 9.86653e-09 8.62318e-07 ----------------------------------------------------------------------------------------------------- 28 H 1 6.545991370125519e-01 6.545991113230690e-01 2.56895e-08 1.05807e-07 2 -8.770520024509981e-01 -8.770520091696136e-01 6.71862e-09 7.82831e-07 3 8.507172027207302e-02 8.507174132968727e-02 2.10576e-08 7.56544e-07 ----------------------------------------------------------------------------------------------------- 29 H 1 4.069174064801466e-01 4.069193197186423e-01 1.91324e-06 2.04114e-05 2 1.753933317378733e+00 1.753934815137202e+00 1.49776e-06 1.11725e-05 3 1.853535317934892e+00 1.853534363848029e+00 9.54087e-07 5.76026e-06 ----------------------------------------------------------------------------------------------------- 30 H 1 7.889152991921473e-01 7.889152182908032e-01 8.09013e-08 1.42663e-07 2 -7.881884291236826e-02 -7.881599658094846e-02 2.84633e-06 3.80304e-06 3 9.847344220484038e-01 9.847337253974027e-01 6.96651e-07 9.19721e-06 ----------------------------------------------------------------------------------------------------- 31 H 1 -1.315175376514616e-01 -1.315178205034677e-01 2.82852e-07 8.90906e-07 2 -2.962443352422692e-01 -2.962427297230456e-01 1.60552e-06 4.07762e-06 3 -6.562697731494400e-01 -6.562691435392617e-01 6.29610e-07 5.06995e-06 ----------------------------------------------------------------------------------------------------- 32 H 1 4.834811530234638e-01 4.834809573822283e-01 1.95641e-07 1.54754e-05 2 -7.862044470462843e-01 -7.862045437483044e-01 9.67020e-08 1.02362e-06 3 -8.935690076320285e-01 -8.935686115769965e-01 3.96055e-07 1.88504e-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 = 14, direction = 3: |F_model - F_numer| error = ----------------------- = 7.86188e-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.6753616333007812 | (Configuration in file "config-O.xyz") ----------------------------------------------------------------------------------------------------- Part Spec Dir Force_model Force_numer |Force diff| uncertainty ----------------------------------------------------------------------------------------------------- 1 O 1 -1.599128119223265e+00 -1.599127286450307e+00 8.32773e-07 2.40685e-06 2 -1.141309538344101e+00 -1.141309772788627e+00 2.34445e-07 4.61238e-06 3 -2.830349108837620e+00 -2.830350141607622e+00 1.03277e-06 7.03319e-03 * ----------------------------------------------------------------------------------------------------- 2 O 1 -5.231858778364951e+00 -5.231859166562520e+00 3.88198e-07 7.16083e-06 2 -1.688684063568031e-01 -1.688681362351535e-01 2.70122e-07 6.45064e-06 3 -1.441588684717872e+00 -1.441588889518727e+00 2.04801e-07 3.14047e-06 ----------------------------------------------------------------------------------------------------- 3 O 1 6.582048303690513e+00 6.582047940817931e+00 3.62873e-07 6.16539e-06 2 -1.021328495916701e+01 -1.021328648558698e+01 1.52642e-06 8.42640e-06 3 1.268305387763344e+01 1.268305392256655e+01 4.49331e-08 4.04175e-06 ----------------------------------------------------------------------------------------------------- 4 O 1 -9.241380472613789e+00 -9.241380165110630e+00 3.07503e-07 8.54512e-06 2 -3.125149423937394e+00 -3.125149563001134e+00 1.39064e-07 5.93508e-06 3 -6.189036532429689e+00 -6.189036170313779e+00 3.62116e-07 5.25314e-02 * ----------------------------------------------------------------------------------------------------- 5 O 1 9.239745818738065e+00 9.239745336515586e+00 4.82222e-07 1.49807e-06 2 -4.182671866669200e+00 -4.182672015087545e+00 1.48418e-07 1.56520e-06 3 -9.249005019430822e+00 -9.249005291411343e+00 2.71981e-07 1.67097e-06 ----------------------------------------------------------------------------------------------------- 6 O 1 -1.171190016610826e-01 -1.171194373517336e-01 4.35691e-07 3.70201e-05 2 -1.167131634017070e+01 -1.167131642586959e+01 8.56989e-08 3.63608e-06 3 4.830423842914408e+00 4.830423829871322e+00 1.30431e-08 1.34091e-06 ----------------------------------------------------------------------------------------------------- 7 O 1 -2.017699558888091e-01 -2.017712857084621e-01 1.32982e-06 1.83795e-06 2 5.845235631378082e+00 5.845238508922037e+00 2.87754e-06 2.82688e-05 3 -6.897938376202765e+00 -6.897936864141017e+00 1.51206e-06 2.04914e-05 ----------------------------------------------------------------------------------------------------- 8 O 1 1.186792115088809e+00 1.186791702367037e+00 4.12722e-07 6.81669e-06 2 5.686193226834543e+00 5.686193202413133e+00 2.44214e-08 2.03627e-06 3 1.271443598997783e+01 1.271443687306680e+01 8.83089e-07 3.56700e-06 ----------------------------------------------------------------------------------------------------- 9 O 1 -3.656226177328317e+00 -3.656226240841203e+00 6.35129e-08 2.28898e-06 2 1.056590943622761e+00 1.056591125187851e+00 1.81565e-07 3.44087e-05 3 -5.737802200793453e+00 -5.737803407948878e+00 1.20716e-06 6.35100e-06 ----------------------------------------------------------------------------------------------------- 10 O 1 -3.845082925428450e+00 -3.845085266304038e+00 2.34088e-06 2.26183e-05 2 2.145158314259848e+00 2.145157322743978e+00 9.91516e-07 2.99791e-06 3 -5.782555771342253e+00 -5.782556713196161e+00 9.41854e-07 1.54281e-06 ----------------------------------------------------------------------------------------------------- 11 O 1 -9.744117139686770e+00 -9.744117858470462e+00 7.18784e-07 5.60562e-06 2 1.245277080027767e+01 1.245277206991510e+01 1.26964e-06 4.62507e-06 3 -3.134884108356737e+00 -3.134884257500755e+00 1.49144e-07 4.34781e-06 ----------------------------------------------------------------------------------------------------- 12 O 1 1.048148950615064e+01 1.048148923524358e+01 2.70907e-07 1.68339e-06 2 -1.125167715208239e+01 -1.125167738807509e+01 2.35993e-07 6.66652e-07 3 -3.879235823019266e+00 -3.879236375374870e+00 5.52356e-07 2.16219e-06 ----------------------------------------------------------------------------------------------------- 13 O 1 1.087791268907075e+01 1.087791401093528e+01 1.32186e-06 4.05448e-06 2 1.189002535791625e+01 1.189002645992039e+01 1.10200e-06 6.73626e-06 3 -3.971089255459870e+00 -3.971089934195550e+00 6.78736e-07 4.73371e-06 ----------------------------------------------------------------------------------------------------- 14 O 1 -4.401499740475585e+00 -4.401496017161842e+00 3.72331e-06 2.56705e-06 2 -1.211189795095222e+01 -1.211189800040748e+01 4.94553e-08 2.76290e-06 3 -2.069040476609762e+00 -2.069040748939000e+00 2.72329e-07 1.01723e-06 ----------------------------------------------------------------------------------------------------- 15 O 1 2.993570996592928e+00 2.993569724899614e+00 1.27169e-06 5.92085e-06 2 8.828696327838184e+00 8.828697639903398e+00 1.31207e-06 2.14644e-05 3 9.745804882451123e+00 9.745803471927589e+00 1.41052e-06 2.00034e-05 ----------------------------------------------------------------------------------------------------- 16 O 1 5.199616183554583e-01 5.199616845845381e-01 6.62291e-08 7.34769e-06 2 8.768473978640895e+00 8.768473850575830e+00 1.28065e-07 9.03923e-06 3 3.432571561875937e+00 3.432574522134572e+00 2.96026e-06 3.82139e-05 ----------------------------------------------------------------------------------------------------- 17 O 1 1.162425128850666e+00 1.162425014634994e+00 1.14216e-07 8.60699e-06 2 -8.631640203569839e+00 -8.631636690288294e+00 3.51328e-06 1.75585e-05 3 -1.152922554267487e+00 -1.152919704836732e+00 2.84943e-06 6.13032e-05 ----------------------------------------------------------------------------------------------------- 18 O 1 6.362607084008934e+00 6.362603765934488e+00 3.31807e-06 8.50407e-05 * 2 -7.399559559223181e-01 -7.399584436117427e-01 2.48769e-06 1.30032e-04 * 3 -3.247025180647324e+00 -3.247019760113421e+00 5.42053e-06 1.11123e-04 * ----------------------------------------------------------------------------------------------------- 19 O 1 4.873998861599495e+00 4.873997704295617e+00 1.15730e-06 2.31948e-06 2 7.549505859540271e+00 7.549500449301033e+00 5.41024e-06 6.94062e-05 * 3 -4.786535341622480e-01 -4.786592379875403e-01 5.70383e-06 2.46545e-05 ----------------------------------------------------------------------------------------------------- 20 O 1 4.659721853880193e+00 4.659724539801098e+00 2.68592e-06 5.64641e-05 2 3.113592334289757e+00 3.113592830670741e+00 4.96381e-07 1.82695e-05 3 1.577930496018311e+00 1.577914520378257e+00 1.59756e-05 7.89164e-05 * ----------------------------------------------------------------------------------------------------- 21 O 1 -9.727130635345141e+00 -9.727130571291765e+00 6.40534e-08 4.18837e-06 2 -9.203791658906092e-01 -9.203783070227090e-01 8.58868e-07 9.24771e-06 3 -9.335609957819697e+00 -9.335608152800512e+00 1.80502e-06 1.75273e-05 ----------------------------------------------------------------------------------------------------- 22 O 1 -3.154747412848879e+00 -3.154746759233141e+00 6.53616e-07 7.76510e-06 2 -1.910518639093415e+00 -1.910519331833368e+00 6.92740e-07 7.51842e-06 3 -6.445023210098299e+00 -6.445021292111542e+00 1.91799e-06 1.14210e-05 ----------------------------------------------------------------------------------------------------- 23 O 1 4.106854231529131e+00 4.106856162945345e+00 1.93142e-06 8.11856e-06 2 -2.886205269702124e+00 -2.886207026493751e+00 1.75679e-06 3.14222e-06 3 2.503098613501055e+00 2.503096783026670e+00 1.83047e-06 1.74523e-06 ----------------------------------------------------------------------------------------------------- 24 O 1 -6.912344372394558e+00 -6.912345264204729e+00 8.91810e-07 3.57812e-05 2 -3.354878845451402e+00 -3.354880150446486e+00 1.30500e-06 2.83601e-06 3 4.666212730068728e+00 4.666214537246636e+00 1.80718e-06 4.39806e-05 ----------------------------------------------------------------------------------------------------- 25 O 1 -9.025627551383257e+00 -9.025629649380578e+00 2.09800e-06 6.24217e-06 2 9.178564051576677e+00 9.178563605714624e+00 4.45862e-07 3.57870e-06 3 5.508874051285388e+00 5.508873718874176e+00 3.32411e-07 6.57861e-06 ----------------------------------------------------------------------------------------------------- 26 O 1 7.157983904350215e+00 7.157984662029970e+00 7.57680e-07 2.78426e-05 2 8.889453007992218e+00 8.889453930582015e+00 9.22590e-07 1.37022e-05 3 1.319359582459264e+01 1.319359546045655e+01 3.64136e-07 3.69123e-06 ----------------------------------------------------------------------------------------------------- 27 O 1 -1.946617868569626e-01 -1.946598516561165e-01 1.93520e-06 4.02782e-05 2 -7.495741568990852e+00 -7.495737733263022e+00 3.83573e-06 1.43279e-05 3 5.980612091670800e+00 5.980610302503486e+00 1.78917e-06 1.66975e-05 ----------------------------------------------------------------------------------------------------- 28 O 1 -2.601590748919671e+00 -2.601591621502930e+00 8.72583e-07 1.59108e-06 2 3.951038161916563e+00 3.951037945368546e+00 2.16548e-07 1.58170e-06 3 4.711456044587807e+00 4.711455451360336e+00 5.93227e-07 4.34541e-06 ----------------------------------------------------------------------------------------------------- 29 O 1 -3.854575996453311e+00 -3.854565453916138e+00 1.05425e-05 1.68233e-05 2 -4.763470332322268e+00 -4.763469349419937e+00 9.82902e-07 1.71155e-05 3 1.076853125366335e+00 1.076853622432977e+00 4.97067e-07 7.08339e-05 * ----------------------------------------------------------------------------------------------------- 30 O 1 1.131148838153574e+01 1.131148672993863e+01 1.65160e-06 1.77271e-05 2 -7.637592417143754e+00 -7.637592575003794e+00 1.57860e-07 5.91421e-06 3 -8.386515483412575e+00 -8.386515628406066e+00 1.44993e-07 2.39268e-06 ----------------------------------------------------------------------------------------------------- 31 O 1 1.835253440460930e+00 1.835251402271460e+00 2.03819e-06 1.26268e-06 2 -1.584014455522282e-01 -1.583984352347803e-01 3.01032e-06 6.53026e-06 3 1.864527476730309e+00 1.864531682280699e+00 4.20555e-06 1.60214e-05 ----------------------------------------------------------------------------------------------------- 32 O 1 -9.842993119029655e+00 -9.842993783120743e+00 6.64091e-07 3.27437e-06 2 3.009661485234898e+00 3.009660294605218e+00 1.19063e-06 3.21040e-05 3 -4.261175331066406e+00 -4.261175068473654e+00 2.62593e-07 1.27346e-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.90047e-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.236818400066450e+00 1.236818397450855e+00 2.61560e-09 2.92914e-08 2 1.064483877061583e+00 1.064483875627062e+00 1.43452e-09 1.58244e-08 3 1.331186136633049e+00 1.331186132701657e+00 3.93139e-09 6.27724e-07 * ----------------------------------------------------------------------------------------------------- 2 V 1 -8.179829151383133e-01 -8.179829149945678e-01 1.43746e-10 1.70812e-08 2 5.556328921363886e-01 5.556328922869147e-01 1.50526e-10 1.73624e-08 3 8.201287504734930e-01 8.201287463741030e-01 4.09939e-09 1.75112e-08 ----------------------------------------------------------------------------------------------------- 3 V 1 -3.344272277568798e+00 -3.344272278064010e+00 4.95211e-10 9.83624e-09 2 2.186333328668974e+00 2.186333328517116e+00 1.51858e-10 5.79621e-09 3 9.747311591108584e-03 9.747306420206804e-03 5.17090e-09 5.23251e-09 ----------------------------------------------------------------------------------------------------- 4 V 1 -1.272298598258877e-01 -1.272298602377128e-01 4.11825e-10 1.39640e-08 2 -1.787072771689490e+00 -1.787072771671914e+00 1.75755e-11 1.32088e-08 3 7.287618446634103e-01 7.287618427654062e-01 1.89800e-09 1.70914e-08 ----------------------------------------------------------------------------------------------------- 5 V 1 1.072803234072957e+00 1.072803236749339e+00 2.67638e-09 1.23378e-08 2 8.471851137116982e-01 8.471851099757503e-01 3.73595e-09 3.42236e-08 3 2.761719394308441e-01 2.761719381587093e-01 1.27213e-09 1.63806e-08 ----------------------------------------------------------------------------------------------------- 6 V 1 -1.558104002795581e-02 -1.558104149357688e-02 1.46562e-09 4.30576e-08 2 -1.992941472287133e+00 -1.992941472518707e+00 2.31573e-10 1.64194e-08 3 -2.290450312665977e+00 -2.290450313409036e+00 7.43059e-10 3.68935e-08 ----------------------------------------------------------------------------------------------------- 7 V 1 8.613484117105011e-02 8.613484218864662e-02 1.01760e-09 2.01674e-08 2 1.026291811694634e+00 1.026291806403838e+00 5.29080e-09 1.42075e-08 3 -3.765528957083471e-01 -3.765528936664333e-01 2.04191e-09 2.77516e-08 ----------------------------------------------------------------------------------------------------- 8 V 1 -3.680187098855653e+00 -3.680187098943001e+00 8.73484e-11 1.68827e-09 2 -2.345569344494637e+00 -2.345569345325393e+00 8.30755e-10 5.91195e-09 3 2.143585072205000e+00 2.143585072234341e+00 2.93419e-11 3.73370e-09 ----------------------------------------------------------------------------------------------------- 9 V 1 4.300794206901559e-01 4.300794214117947e-01 7.21639e-10 4.07005e-08 2 -1.118536636292701e+00 -1.118536632391404e+00 3.90130e-09 7.64140e-09 3 2.043019282394194e-01 2.043019257969299e-01 2.44249e-09 9.26092e-09 ----------------------------------------------------------------------------------------------------- 10 V 1 1.490072216375615e+00 1.490072219066756e+00 2.69114e-09 1.14450e-08 2 2.420372777663663e+00 2.420372776243478e+00 1.42019e-09 1.06346e-08 3 -1.917743571406900e+00 -1.917743571996253e+00 5.89353e-10 3.74182e-07 * ----------------------------------------------------------------------------------------------------- 11 V 1 -2.366169199541193e+00 -2.366169199872335e+00 3.31142e-10 4.75486e-09 2 -2.395765647215105e+00 -2.395765645840635e+00 1.37447e-09 3.73681e-09 3 3.926855773037938e+00 3.926855772535330e+00 5.02608e-10 2.09525e-09 ----------------------------------------------------------------------------------------------------- 12 V 1 -7.961261760239560e-01 -7.961261805524893e-01 4.52853e-09 3.18585e-08 2 -1.539014902063099e-02 -1.539015090272376e-02 1.88209e-09 6.49138e-09 3 1.164397331297934e+00 1.164397330236011e+00 1.06192e-09 6.12380e-09 ----------------------------------------------------------------------------------------------------- 13 V 1 2.949990744466291e+00 2.949990746682468e+00 2.21618e-09 7.11486e-09 2 5.040263996946233e-01 5.040263993316305e-01 3.62993e-10 4.25466e-09 3 -2.732758300984429e+00 -2.732758298966643e+00 2.01779e-09 1.94656e-08 ----------------------------------------------------------------------------------------------------- 14 V 1 -1.257776717930842e+00 -1.257776716832672e+00 1.09817e-09 3.64874e-08 2 -1.267351650859132e+00 -1.267351653688155e+00 2.82902e-09 6.30398e-09 3 1.051054730025605e+00 1.051054724165908e+00 5.85970e-09 5.24769e-09 ----------------------------------------------------------------------------------------------------- 15 V 1 4.848014677727216e-01 4.848014666890634e-01 1.08366e-09 1.68207e-08 2 4.268426333653016e+00 4.268426327866499e+00 5.78652e-09 4.91072e-08 3 2.420947534627874e+00 2.420947537517404e+00 2.88953e-09 1.08101e-08 ----------------------------------------------------------------------------------------------------- 16 V 1 1.410980680969917e+00 1.410980679173004e+00 1.79691e-09 3.49048e-08 2 8.126228728963369e-01 8.126228704640713e-01 2.43227e-09 1.63400e-08 3 -2.102741998099545e+00 -2.102742001399663e+00 3.30012e-09 1.82936e-08 ----------------------------------------------------------------------------------------------------- 17 V 1 2.293375343117481e+00 2.293375340925936e+00 2.19155e-09 9.05240e-09 2 1.589287476137597e+00 1.589287473764620e+00 2.37298e-09 6.95103e-09 3 -1.095881189974877e+00 -1.095881189218528e+00 7.56349e-10 1.64849e-06 * ----------------------------------------------------------------------------------------------------- 18 V 1 9.362461526334613e-01 9.362461524730620e-01 1.60399e-10 1.11318e-09 2 -2.487151307147811e+00 -2.487151308457550e+00 1.30974e-09 7.18848e-09 3 1.745984708997354e+00 1.745984708553966e+00 4.43387e-10 4.31747e-09 ----------------------------------------------------------------------------------------------------- 19 V 1 -8.820411443296646e-01 -8.820411459155808e-01 1.58592e-09 2.60843e-08 2 1.311680068832431e+00 1.311680067475500e+00 1.35693e-09 3.13988e-08 3 -1.608531776367124e+00 -1.608531775402359e+00 9.64764e-10 2.38694e-08 ----------------------------------------------------------------------------------------------------- 20 V 1 1.653653588294230e+00 1.653653586782445e+00 1.51178e-09 2.87533e-08 2 4.775335022011183e-01 4.775335015800508e-01 6.21067e-10 6.95344e-09 3 -1.475118911366724e+00 -1.475118918336118e+00 6.96939e-09 2.14107e-09 ----------------------------------------------------------------------------------------------------- 21 V 1 -2.144740802772822e-01 -2.144740803090199e-01 3.17377e-11 5.04535e-09 2 1.716015091742761e+00 1.716015091483358e+00 2.59403e-10 8.13887e-09 3 2.133852267165491e+00 2.133852267493739e+00 3.28249e-10 8.89469e-09 ----------------------------------------------------------------------------------------------------- 22 V 1 4.997513775989529e-01 4.997513728308620e-01 4.76809e-09 1.27097e-08 2 2.259553112941213e-02 2.259553162437688e-02 4.94965e-10 2.22207e-08 3 6.066919674394641e-01 6.066919670579800e-01 3.81484e-10 1.21384e-08 ----------------------------------------------------------------------------------------------------- 23 V 1 1.815659150710701e-01 1.815659153750709e-01 3.04001e-10 1.57938e-08 2 -6.266643004057507e-01 -6.266643015351057e-01 1.12935e-09 2.57715e-08 3 -1.148786717138411e+00 -1.148786711733322e+00 5.40509e-09 1.40233e-08 ----------------------------------------------------------------------------------------------------- 24 V 1 -1.740394734438005e+00 -1.740394735347362e+00 9.09357e-10 5.82165e-09 2 1.461985132315859e+00 1.461985131685528e+00 6.30331e-10 8.99240e-09 3 -1.691867747159898e+00 -1.691867749552979e+00 2.39308e-09 1.63850e-08 ----------------------------------------------------------------------------------------------------- 25 V 1 1.525841746066327e+00 1.525841745823830e+00 2.42498e-10 3.01181e-07 * 2 -3.453935421311343e+00 -3.453935421948079e+00 6.36736e-10 3.60721e-09 3 -2.226993187524576e+00 -2.226993188987115e+00 1.46254e-09 1.93613e-08 ----------------------------------------------------------------------------------------------------- 26 V 1 1.562949170854191e+00 1.562949172424160e+00 1.56997e-09 9.80257e-09 2 -1.098857300421052e+00 -1.098857301965768e+00 1.54472e-09 2.71527e-09 3 -2.384651277664308e-01 -2.384651273773629e-01 3.89068e-10 4.16895e-08 ----------------------------------------------------------------------------------------------------- 27 V 1 -4.756639009549244e-01 -4.756639018707142e-01 9.15790e-10 2.83130e-08 2 3.624654268720389e-01 3.624654271710896e-01 2.99051e-10 1.44395e-08 3 -7.649677586206954e-01 -7.649677592481297e-01 6.27434e-10 1.45491e-08 ----------------------------------------------------------------------------------------------------- 28 V 1 -9.086265993529117e-01 -9.086266045631369e-01 5.21023e-09 1.86524e-08 2 1.440818937043685e+00 1.440818936813880e+00 2.29805e-10 8.04937e-09 3 3.347487936698132e+00 3.347487934259372e+00 2.43876e-09 7.21917e-09 ----------------------------------------------------------------------------------------------------- 29 V 1 -2.190275786016730e+00 -2.190275778301580e+00 7.71515e-09 4.57828e-08 2 -4.298110299312251e+00 -4.298110298801234e+00 5.11017e-10 2.13881e-08 3 -3.411919533068779e+00 -3.411919534377377e+00 1.30860e-09 6.38955e-09 ----------------------------------------------------------------------------------------------------- 30 V 1 2.443335990157579e+00 2.443335989378056e+00 7.79523e-10 4.10996e-09 2 1.868971850951532e+00 1.868971852405360e+00 1.45383e-09 2.13968e-08 3 -3.333621418088859e+00 -3.333621419842316e+00 1.75346e-09 1.20025e-08 ----------------------------------------------------------------------------------------------------- 31 V 1 3.181539560567999e+00 3.181539559250067e+00 1.31793e-09 1.19149e-08 2 -5.161998199610539e+00 -5.161998198669647e+00 9.40893e-10 2.71585e-08 3 2.660146550501219e+00 2.660146556651890e+00 6.15067e-09 1.01577e-08 ----------------------------------------------------------------------------------------------------- 32 V 1 -4.623138319664333e+00 -4.623138323591897e+00 3.92756e-09 1.14447e-08 2 4.112616075660220e+00 4.112616073848137e+00 1.81208e-09 2.19102e-08 3 1.845098662914237e+00 1.845098661728255e+00 1.18598e-09 2.29695e-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.30495e-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.1554128589923494 | (Configuration in file "config-CHOV.xyz") ----------------------------------------------------------------------------------------------------- Part Spec Dir Force_model Force_numer |Force diff| uncertainty ----------------------------------------------------------------------------------------------------- 1 C 1 7.815688028447560e-01 7.459765825080081e-01 3.55922e-02 3.33347e-03 2 4.212076139228996e+00 4.149605747478576e+00 6.24704e-02 1.14596e-03 3 2.981828502911518e+00 2.963198583180875e+00 1.86299e-02 1.63717e-03 ----------------------------------------------------------------------------------------------------- 2 H 1 -4.571919226523720e-01 -4.772367473997269e-01 2.00448e-02 2.06855e-03 2 6.579346783406703e-01 6.516216481505375e-01 6.31303e-03 2.50908e-03 3 1.780910287105187e+00 1.796994958233179e+00 1.60847e-02 1.42631e-03 ----------------------------------------------------------------------------------------------------- 3 H 1 4.416894532624104e-01 4.103036243960856e-01 3.13858e-02 1.92191e-02 2 -1.227766675575923e+00 -1.138600636803293e+00 8.91660e-02 1.87480e-01 * 3 1.213373386695059e+00 1.188867548910618e+00 2.45058e-02 1.92584e-02 ----------------------------------------------------------------------------------------------------- 4 O 1 -7.402561140631941e-01 -6.230938851491110e-01 1.17162e-01 1.45975e-01 2 -4.525232854104099e+00 -4.499340411806078e+00 2.58924e-02 4.63170e-03 3 -2.263636533161463e+00 -2.210266808271022e+00 5.33697e-02 4.37340e-03 ----------------------------------------------------------------------------------------------------- 5 H 1 1.986298855715020e+00 1.882427364200545e+00 1.03871e-01 3.02522e-02 2 8.119895027938367e-01 8.054599346876340e-01 6.52957e-03 1.71481e-03 3 1.434749850232996e+00 1.429274092621149e+00 5.47576e-03 1.31379e-01 ----------------------------------------------------------------------------------------------------- 6 C 1 -9.332231975120060e-01 -8.782192946427529e-01 5.50039e-02 1.63793e-03 2 -1.876155202274089e+00 -1.932817030595843e+00 5.66618e-02 2.74920e-02 3 7.543730665390695e-01 7.617760517618695e-01 7.40299e-03 6.67368e-03 ----------------------------------------------------------------------------------------------------- 7 V 1 4.258463650859290e+00 4.279251127816324e+00 2.07875e-02 4.26153e-02 2 -1.965349987921357e+00 -1.803318719303767e+00 1.62031e-01 3.09164e-01 * 3 -6.348132661939038e+00 -6.309759616296569e+00 3.83730e-02 2.50125e-04 ----------------------------------------------------------------------------------------------------- 8 C 1 6.364007554110777e-01 6.898465082292955e-01 5.34458e-02 6.40397e-03 2 3.615897223432640e+00 3.362149709654889e+00 2.53748e-01 1.56473e-03 3 -3.034772554472556e+00 -2.998617919415446e+00 3.61546e-02 1.37554e-02 ----------------------------------------------------------------------------------------------------- 9 C 1 6.460008588334182e+00 6.451201419407960e+00 8.80717e-03 1.10148e-02 2 3.151787820862651e+00 3.157282109160611e+00 5.49429e-03 7.04325e-03 3 2.241276123637857e+00 2.225960171337057e+00 1.53160e-02 2.55047e-03 ----------------------------------------------------------------------------------------------------- 10 C 1 -3.123589823375978e+00 -3.123269976010855e+00 3.19847e-04 2.00761e-02 2 -4.961086469775970e+00 -4.959974033431735e+00 1.11244e-03 2.74366e-02 3 8.554773562860776e-01 8.381257494276161e-01 1.73516e-02 1.69652e-01 * ----------------------------------------------------------------------------------------------------- 11 H 1 -3.142842761272851e+00 -3.122800301105705e+00 2.00425e-02 1.17839e-03 2 -3.065948102647995e-01 -2.932670680560900e-01 1.33277e-02 2.30821e-02 3 -1.899056041193263e+00 -1.898813411054385e+00 2.42630e-04 9.81143e-04 ----------------------------------------------------------------------------------------------------- 12 H 1 2.284167605048476e-02 6.468174584572280e-02 4.18401e-02 1.45171e-02 2 -2.783750236902003e+00 -2.805822042612370e+00 2.20718e-02 3.13285e-02 3 2.159335654292953e+00 2.159416417382677e+00 8.07631e-05 1.40786e-03 ----------------------------------------------------------------------------------------------------- 13 H 1 -7.377923805980278e-01 -7.357766591021431e-01 2.01572e-03 3.74949e-04 2 1.634725979372122e+00 1.623721338239327e+00 1.10046e-02 4.64696e-02 3 4.460285985708065e-01 5.249519312396589e-01 7.89233e-02 1.19444e-01 ----------------------------------------------------------------------------------------------------- 14 V 1 -2.123702235531213e+00 -2.106862320366585e+00 1.68399e-02 7.09595e-03 2 -2.490668376237226e+00 -2.492264720177279e+00 1.59634e-03 8.77550e-03 3 3.449547862894084e+00 3.464731619031230e+00 1.51838e-02 8.82542e-03 ----------------------------------------------------------------------------------------------------- 15 O 1 9.379517004592277e-01 9.356413532101966e-01 2.31035e-03 3.23095e-02 2 2.845255762960155e+00 2.851576801351466e+00 6.32104e-03 4.08177e-02 3 -2.536601025392491e+00 -2.551650589099172e+00 1.50496e-02 8.49672e-03 ----------------------------------------------------------------------------------------------------- 16 O 1 2.688459828809469e+00 2.657416977680380e+00 3.10429e-02 4.96587e-03 2 5.073878592756044e-01 4.863094138128176e-01 2.10784e-02 8.17765e-03 3 6.241038141414224e-01 6.062238251764491e-01 1.78800e-02 1.81167e-03 ----------------------------------------------------------------------------------------------------- 17 V 1 5.087399430504562e+00 5.029262004879661e+00 5.81374e-02 1.10913e-02 2 -1.472500561307611e+01 -1.472434896067862e+01 6.56652e-04 5.63219e-03 3 -1.581957472692129e+01 -1.581020513467606e+01 9.36959e-03 1.28422e-02 ----------------------------------------------------------------------------------------------------- 18 O 1 -2.970232914299733e+00 -2.993377981718492e+00 2.31451e-02 6.86247e-03 2 2.069728318921646e+00 2.103401891314533e+00 3.36736e-02 9.41539e-02 3 3.263038747580807e+00 3.223487908055457e+00 3.95508e-02 7.34715e-02 ----------------------------------------------------------------------------------------------------- 19 O 1 -2.075166437899808e+00 -2.099644918898271e+00 2.44785e-02 2.42320e-02 2 2.454729856900848e+00 2.474902986101865e+00 2.01731e-02 1.30315e-01 3 -2.070078251983990e+00 -2.044145218359189e+00 2.59330e-02 6.07411e-02 ----------------------------------------------------------------------------------------------------- 20 V 1 -2.847509959277799e+00 -2.842813989072444e+00 4.69597e-03 2.05902e-03 2 1.001223522446058e+01 1.000864579485421e+01 3.58943e-03 3.62480e-03 3 1.669408263045769e+01 1.669203415452121e+01 2.04848e-03 2.12197e-03 ----------------------------------------------------------------------------------------------------- 21 V 1 -1.926187969228354e+00 -1.847458617382073e+00 7.87294e-02 4.61528e-02 2 -3.329070582366978e+00 -3.242837333294219e+00 8.62332e-02 2.33502e-04 3 1.944769637847743e+00 2.009653638592996e+00 6.48840e-02 1.04680e-02 ----------------------------------------------------------------------------------------------------- 22 C 1 -1.726848157953148e+00 -1.741259095584401e+00 1.44109e-02 8.24183e-03 2 -3.495611520272896e+00 -3.540967006548143e+00 4.53555e-02 3.40874e-02 3 -4.684590927906001e+00 -4.695659054980390e+00 1.10681e-02 3.32509e-02 ----------------------------------------------------------------------------------------------------- 23 C 1 -2.579939484852683e+00 -2.654321234269821e+00 7.43817e-02 1.36345e-01 2 -3.069827371680238e-01 -2.986546928657787e-01 8.32804e-03 8.47005e-03 3 -4.057333359376222e-01 -4.272246688975879e-01 2.14913e-02 3.33133e-01 * ----------------------------------------------------------------------------------------------------- 24 V 1 6.141543143905478e+00 6.246948220505538e+00 1.05405e-01 1.96932e-02 2 1.704479911710969e+00 1.650513160883238e+00 5.39668e-02 1.80076e-03 3 5.279885348665615e+00 5.207638069519687e+00 7.22473e-02 1.34177e-01 ----------------------------------------------------------------------------------------------------- 25 V 1 -4.798236184901178e+00 -4.874888930528881e+00 7.66527e-02 3.53821e-03 2 9.209360766363959e+00 9.123194436852947e+00 8.61663e-02 4.01069e-02 3 -1.194111844964698e+01 -1.184695565511239e+01 9.41628e-02 7.38745e-02 ----------------------------------------------------------------------------------------------------- 26 H 1 -1.635630869333867e+00 -1.611669691646208e+00 2.39612e-02 2.12273e-03 2 -3.606354935717907e+00 -3.565099174559064e+00 4.12558e-02 1.09492e-02 3 1.964263740169780e+00 1.941995656262719e+00 2.22681e-02 5.16263e-03 ----------------------------------------------------------------------------------------------------- 27 V 1 8.316127757718981e+00 8.260947310089280e+00 5.51804e-02 2.40117e-02 2 8.075307001736807e+00 7.986439788596154e+00 8.88672e-02 1.89015e-02 3 3.348062275143199e+00 3.393959448489979e+00 4.58972e-02 1.41479e-02 ----------------------------------------------------------------------------------------------------- 28 O 1 1.513930732166928e+00 1.558329234167254e+00 4.43985e-02 8.05642e-02 2 -4.160902798271293e+00 -4.167963380587325e+00 7.06058e-03 2.15389e-02 3 -1.509459016526275e+00 -1.551339941399415e+00 4.18809e-02 4.10677e-02 ----------------------------------------------------------------------------------------------------- 29 O 1 5.619593659276481e-01 4.844040145636443e-01 7.75554e-02 1.84352e-02 2 1.953327310950089e+00 1.959858147263585e+00 6.53084e-03 3.86465e-01 * 3 1.858574294004346e+00 1.877248842156988e+00 1.86745e-02 1.25228e-02 ----------------------------------------------------------------------------------------------------- 30 C 1 -5.795196311722608e+00 -5.762801214742944e+00 3.23951e-02 1.06708e-01 2 -2.829631481223758e+00 -2.847447508250043e+00 1.78160e-02 3.36645e-01 * 3 -6.462554410893828e-01 -6.640605300517615e-01 1.78051e-02 2.36839e-03 ----------------------------------------------------------------------------------------------------- 31 O 1 -1.674001791467561e+00 -1.681353534159127e+00 7.35174e-03 1.76258e-01 * 2 1.124726767354712e+00 1.165070646999589e+00 4.03439e-02 1.63488e-02 3 2.125340384810608e+00 2.166089539670497e+00 4.07492e-02 2.09774e-03 ----------------------------------------------------------------------------------------------------- 32 H 1 -5.470952260271446e-01 -5.511479838655531e-01 4.05276e-03 1.39344e-01 2 -1.450785843513844e+00 -1.436893581061133e+00 1.38923e-02 2.59200e-03 3 -1.260012595816461e+00 -1.285199752961506e+00 2.51872e-02 4.42114e-03 ----------------------------------------------------------------------------------------------------- * 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 = 12, direction = 1: |F_model - F_numer| error = ----------------------- = 6.46861e-01 max{|F_numer|, eps} ===================================================================================================== Grade is based on maximum error across all systems (for a model supporting multiple species): Maximum error = 6.46861e-01 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: F Comment: The letter grade F was assigned because the normalized error in the computation was 6.46861e-01 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.