!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!! !!!!! !!!!! VERIFICATION CHECK: vc-dimer-continuity-c1 !!!!! !!!!! !!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Description: Determines whether a model has a continuous energy and first derivative, i.e. belongs to the C^1 continuity class, for all possible dimers. For a model supporting N species, there are N + N!/(2(N-2)! distinct dimers for all possible species combinations. (For example if N=3, there are 3+3!/2=6 dimers. If the species are A,B,C, the 6 dimers are AA, BB, CC, AB, AC, BC.) For each dimer, the equilibrium separation and cutoff are determined. The continuity across the cutoff is assessed. Then an analysis is performed to detect any discontinuities from half the equilibrium distance to the cutoff. Although the verification check only requires C^1 continuity to pass, continuity up to 3rd order is checked and reported. Author: Ellad Tadmor ----------------------------------------------------------------------------------------------------- Results for KIM Model : EAM_Dynamo_BonnyTerentyev_2014EAM2_W__MO_626183701337_000 Supported species : H He W ----------------------------------------------------------------------------------------------------- DIMER H--H ----------------------------------------------------------------------------------------------------- E n e r g y E x t r e m a distance energy #iter #fn-calls warnings Equilibrium separation: 3.35000000e+00 0.00000000e+00 33 92 none Cutoff separation: 2.99996536e+00 2.14164629e-16 29 31 none C u t o f f S m o o t h n e s s value continuous Energy 2.14164629e-16 yes 1st-Deriv -2.83989926e-29 yes 2nd-Deriv 1.81753553e-27 yes 3rd-Deriv -2.58493941e-25 yes C o n t i n u i t y Checking continuity for r = [1.67500,2.99997] at 133 points (Delta r = 0.01000) Local edge detection based on a normalized 5th-order local difference formula T^5 is used to determine the presence of discontinuities. The tolernace is |T^5|>1.00000. (For details see Anne Gelb and Eitan Tadmor, J. Sci. Comp., 28:279-306, 2006.) Checking Energy (energy vs. distance in file "dimer-energy-HH.dat") ... No discontinuities found. Checking 1st-Deriv (1st-deriv vs. distance in file "dimer-energy-deriv-1-HH.dat") ... No discontinuities found. Checking 2nd-Deriv (2nd-deriv vs. distance in file "dimer-energy-deriv-2-HH.dat") ... No discontinuities found. Checking 3rd-Deriv (3rd-deriv vs. distance in file "dimer-energy-deriv-3-HH.dat") ==> Suspected discontinuity encountered at r=1.99500000e+00 (|T^5| = 1.38950755e+00) ==> Suspected discontinuity encountered at r=2.19500000e+00 (|T^5| = 1.19997064e+00) ----------------------------------------------------------------------------------------------------- DIMER He--He ----------------------------------------------------------------------------------------------------- E n e r g y E x t r e m a distance energy #iter #fn-calls warnings Equilibrium separation: 1.60642981e+00 -1.85112552e-01 29 58 none Cutoff separation: 3.00214016e+00 2.22355657e-16 25 27 none C u t o f f S m o o t h n e s s value continuous Energy 2.22355657e-16 yes 1st-Deriv -4.40162820e-07 yes 2nd-Deriv -4.11337381e-04 ** NO ** 3rd-Deriv -1.92199360e-01 ** NO ** C o n t i n u i t y Checking continuity for r = [0.80321,3.00214] at 220 points (Delta r = 0.01000) Local edge detection based on a normalized 5th-order local difference formula T^5 is used to determine the presence of discontinuities. The tolernace is |T^5|>1.00000. (For details see Anne Gelb and Eitan Tadmor, J. Sci. Comp., 28:279-306, 2006.) Checking Energy (energy vs. distance in file "dimer-energy-HeHe.dat") ... No discontinuities found. Checking 1st-Deriv (1st-deriv vs. distance in file "dimer-energy-deriv-1-HeHe.dat") ... No discontinuities found. Checking 2nd-Deriv (2nd-deriv vs. distance in file "dimer-energy-deriv-2-HeHe.dat") ... No discontinuities found. Checking 3rd-Deriv (3rd-deriv vs. distance in file "dimer-energy-deriv-3-HeHe.dat") ==> Suspected discontinuity encountered at r=1.77321490e+00 (|T^5| = 1.99464416e+01) ==> Suspected discontinuity encountered at r=1.78321490e+00 (|T^5| = 7.97272796e+01) ==> Suspected discontinuity encountered at r=1.79321490e+00 (|T^5| = 1.19453905e+02) ==> Suspected discontinuity encountered at r=1.80321490e+00 (|T^5| = 7.94540055e+01) ==> Suspected discontinuity encountered at r=1.81321490e+00 (|T^5| = 1.97270394e+01) ==> Suspected discontinuity encountered at r=1.97321490e+00 (|T^5| = 1.04903483e+00) ==> Suspected discontinuity encountered at r=1.98321490e+00 (|T^5| = 4.19400247e+00) ==> Suspected discontinuity encountered at r=1.99321490e+00 (|T^5| = 6.28537677e+00) ==> Suspected discontinuity encountered at r=2.00321490e+00 (|T^5| = 4.18266777e+00) ==> Suspected discontinuity encountered at r=2.01321490e+00 (|T^5| = 1.03995831e+00) ==> Suspected discontinuity encountered at r=2.07321490e+00 (|T^5| = 3.99573319e+00) ==> Suspected discontinuity encountered at r=2.08321490e+00 (|T^5| = 1.59782668e+01) ==> Suspected discontinuity encountered at r=2.09321490e+00 (|T^5| = 2.39565240e+01) ==> Suspected discontinuity encountered at r=2.10321490e+00 (|T^5| = 1.59566182e+01) ==> Suspected discontinuity encountered at r=2.11321490e+00 (|T^5| = 3.97839345e+00) ----------------------------------------------------------------------------------------------------- DIMER W--W ----------------------------------------------------------------------------------------------------- E n e r g y E x t r e m a distance energy #iter #fn-calls warnings Equilibrium separation: 2.34822783e+00 -7.65269190e+00 32 64 none Cutoff separation: 5.45950146e+00 -2.88760060e-14 28 30 none C u t o f f S m o o t h n e s s value continuous Energy -2.88760060e-14 yes 1st-Deriv -4.03331746e-06 ** NO ** 2nd-Deriv 7.64089071e-03 ** NO ** 3rd-Deriv -8.47524499e+00 ** NO ** C o n t i n u i t y Checking continuity for r = [1.17411,5.45950] at 429 points (Delta r = 0.01000) Local edge detection based on a normalized 5th-order local difference formula T^5 is used to determine the presence of discontinuities. The tolernace is |T^5|>1.00000. (For details see Anne Gelb and Eitan Tadmor, J. Sci. Comp., 28:279-306, 2006.) Checking Energy (energy vs. distance in file "dimer-energy-WW.dat") ... No discontinuities found. Checking 1st-Deriv (1st-deriv vs. distance in file "dimer-energy-deriv-1-WW.dat") ... No discontinuities found. Checking 2nd-Deriv (2nd-deriv vs. distance in file "dimer-energy-deriv-2-WW.dat") ==> Suspected discontinuity encountered at r=1.97411392e+00 (|T^5| = 7.52662431e+00) ==> Suspected discontinuity encountered at r=1.98411392e+00 (|T^5| = 3.02643719e+01) ==> Suspected discontinuity encountered at r=1.99411392e+00 (|T^5| = 4.55936668e+01) ==> Suspected discontinuity encountered at r=2.00411392e+00 (|T^5| = 3.04584339e+01) ==> Suspected discontinuity encountered at r=2.01411392e+00 (|T^5| = 7.55844741e+00) ==> Suspected discontinuity encountered at r=2.55411392e+00 (|T^5| = 2.58041255e+00) ==> Suspected discontinuity encountered at r=2.56411392e+00 (|T^5| = 2.42927888e+00) ==> Suspected discontinuity encountered at r=2.61411392e+00 (|T^5| = 1.34288827e+00) ==> Suspected discontinuity encountered at r=2.63411392e+00 (|T^5| = 2.34663982e+00) ==> Suspected discontinuity encountered at r=2.64411392e+00 (|T^5| = 1.04823916e+00) ==> Suspected discontinuity encountered at r=2.68411392e+00 (|T^5| = 2.55441528e+00) ==> Suspected discontinuity encountered at r=2.69411392e+00 (|T^5| = 2.44591739e+00) ==> Suspected discontinuity encountered at r=3.81411392e+00 (|T^5| = 1.14407238e+00) ==> Suspected discontinuity encountered at r=3.88411392e+00 (|T^5| = 2.07567621e+00) ==> Suspected discontinuity encountered at r=3.89411392e+00 (|T^5| = 1.85020077e+00) ==> Suspected discontinuity encountered at r=3.94411392e+00 (|T^5| = 3.48254360e+00) ==> Suspected discontinuity encountered at r=3.95411392e+00 (|T^5| = 6.84183577e+00) ==> Suspected discontinuity encountered at r=3.96411392e+00 (|T^5| = 5.73075391e+00) ==> Suspected discontinuity encountered at r=3.97411392e+00 (|T^5| = 1.75038619e+00) ==> Suspected discontinuity encountered at r=4.02411392e+00 (|T^5| = 1.10080946e+00) ==> Suspected discontinuity encountered at r=4.03411392e+00 (|T^5| = 3.61128805e+00) ==> Suspected discontinuity encountered at r=4.04411392e+00 (|T^5| = 3.92952486e+00) ==> Suspected discontinuity encountered at r=4.05411392e+00 (|T^5| = 1.49583643e+00) ==> Suspected discontinuity encountered at r=4.12411392e+00 (|T^5| = 1.14605237e+00) ==> Suspected discontinuity encountered at r=4.13411392e+00 (|T^5| = 4.60025153e+00) ==> Suspected discontinuity encountered at r=4.14411392e+00 (|T^5| = 6.81184764e+00) ==> Suspected discontinuity encountered at r=4.15411392e+00 (|T^5| = 4.27387242e+00) ==> Suspected discontinuity encountered at r=4.29411392e+00 (|T^5| = 4.24317027e+00) ==> Suspected discontinuity encountered at r=4.30411392e+00 (|T^5| = 7.06032892e+00) ==> Suspected discontinuity encountered at r=4.31411392e+00 (|T^5| = 5.22493552e+00) ==> Suspected discontinuity encountered at r=4.32411392e+00 (|T^5| = 1.45180453e+00) Checking 3rd-Deriv (3rd-deriv vs. distance in file "dimer-energy-deriv-3-WW.dat") ==> Suspected discontinuity encountered at r=1.97411392e+00 (|T^5| = 1.47346632e+03) ==> Suspected discontinuity encountered at r=1.98411392e+00 (|T^5| = 7.38527154e+03) ==> Suspected discontinuity encountered at r=1.99411392e+00 (|T^5| = 1.48091274e+04) ==> Suspected discontinuity encountered at r=2.00411392e+00 (|T^5| = 1.48489226e+04) ==> Suspected discontinuity encountered at r=2.01411392e+00 (|T^5| = 7.44479551e+03) ==> Suspected discontinuity encountered at r=2.02411392e+00 (|T^5| = 1.49308194e+03) ==> Suspected discontinuity encountered at r=2.47411392e+00 (|T^5| = 3.21973708e+01) ==> Suspected discontinuity encountered at r=2.48411392e+00 (|T^5| = 1.27570792e+02) ==> Suspected discontinuity encountered at r=2.49411392e+00 (|T^5| = 1.89582040e+02) ==> Suspected discontinuity encountered at r=2.50411392e+00 (|T^5| = 1.25352401e+02) ==> Suspected discontinuity encountered at r=2.51411392e+00 (|T^5| = 3.10994332e+01) ==> Suspected discontinuity encountered at r=2.53411392e+00 (|T^5| = 1.64470054e+01) ==> Suspected discontinuity encountered at r=2.54411392e+00 (|T^5| = 1.20798628e+01) ==> Suspected discontinuity encountered at r=2.55411392e+00 (|T^5| = 2.15908065e+02) ==> Suspected discontinuity encountered at r=2.56411392e+00 (|T^5| = 4.07872939e+02) ==> Suspected discontinuity encountered at r=2.57411392e+00 (|T^5| = 2.99722003e+02) ==> Suspected discontinuity encountered at r=2.58411392e+00 (|T^5| = 7.87693755e+01) ==> Suspected discontinuity encountered at r=2.60411392e+00 (|T^5| = 1.84904817e+02) ==> Suspected discontinuity encountered at r=2.61411392e+00 (|T^5| = 7.38653541e+02) ==> Suspected discontinuity encountered at r=2.62411392e+00 (|T^5| = 1.10838418e+03) ==> Suspected discontinuity encountered at r=2.63411392e+00 (|T^5| = 7.39738254e+02) ==> Suspected discontinuity encountered at r=2.64411392e+00 (|T^5| = 1.85422513e+02) ==> Suspected discontinuity encountered at r=2.66411392e+00 (|T^5| = 2.40736275e+01) ==> Suspected discontinuity encountered at r=2.67411392e+00 (|T^5| = 2.42892223e+01) ==> Suspected discontinuity encountered at r=2.68411392e+00 (|T^5| = 1.38996053e+02) ==> Suspected discontinuity encountered at r=2.69411392e+00 (|T^5| = 3.26139066e+02) ==> Suspected discontinuity encountered at r=2.70411392e+00 (|T^5| = 2.56769522e+02) ==> Suspected discontinuity encountered at r=2.71411392e+00 (|T^5| = 7.02496085e+01) ==> Suspected discontinuity encountered at r=2.84411392e+00 (|T^5| = 8.22320368e+00) ==> Suspected discontinuity encountered at r=2.85411392e+00 (|T^5| = 3.29274090e+01) ==> Suspected discontinuity encountered at r=2.86411392e+00 (|T^5| = 4.91625505e+01) ==> Suspected discontinuity encountered at r=2.87411392e+00 (|T^5| = 3.31043236e+01) ==> Suspected discontinuity encountered at r=2.88411392e+00 (|T^5| = 8.98152981e+00) ==> Suspected discontinuity encountered at r=2.94411392e+00 (|T^5| = 7.12269203e+00) ==> Suspected discontinuity encountered at r=2.95411392e+00 (|T^5| = 2.64928089e+01) ==> Suspected discontinuity encountered at r=2.96411392e+00 (|T^5| = 3.66355396e+01) ==> Suspected discontinuity encountered at r=2.97411392e+00 (|T^5| = 2.08023657e+01) ==> Suspected discontinuity encountered at r=2.98411392e+00 (|T^5| = 2.12191391e+00) ==> Suspected discontinuity encountered at r=2.99411392e+00 (|T^5| = 2.04273802e+00) ==> Suspected discontinuity encountered at r=3.00411392e+00 (|T^5| = 1.06292123e+00) ==> Suspected discontinuity encountered at r=3.04411392e+00 (|T^5| = 1.06141976e+00) ==> Suspected discontinuity encountered at r=3.05411392e+00 (|T^5| = 2.75954381e+00) ==> Suspected discontinuity encountered at r=3.06411392e+00 (|T^5| = 7.38421973e+00) ==> Suspected discontinuity encountered at r=3.07411392e+00 (|T^5| = 1.73641444e+01) ==> Suspected discontinuity encountered at r=3.08411392e+00 (|T^5| = 3.48613645e+01) ==> Suspected discontinuity encountered at r=3.09411392e+00 (|T^5| = 4.36837055e+01) ==> Suspected discontinuity encountered at r=3.10411392e+00 (|T^5| = 2.83299550e+01) ==> Suspected discontinuity encountered at r=3.11411392e+00 (|T^5| = 8.04527223e+00) ==> Suspected discontinuity encountered at r=3.12411392e+00 (|T^5| = 1.38026111e+00) ==> Suspected discontinuity encountered at r=3.13411392e+00 (|T^5| = 1.40817901e+00) ==> Suspected discontinuity encountered at r=3.14411392e+00 (|T^5| = 1.47500020e+00) ==> Suspected discontinuity encountered at r=3.15411392e+00 (|T^5| = 1.65233221e+00) ==> Suspected discontinuity encountered at r=3.16411392e+00 (|T^5| = 1.70608953e+00) ==> Suspected discontinuity encountered at r=3.17411392e+00 (|T^5| = 1.13543666e+00) ==> Suspected discontinuity encountered at r=3.20411392e+00 (|T^5| = 1.33389723e+00) ==> Suspected discontinuity encountered at r=3.21411392e+00 (|T^5| = 1.60830190e+00) ==> Suspected discontinuity encountered at r=3.22411392e+00 (|T^5| = 1.81339637e+00) ==> Suspected discontinuity encountered at r=3.23411392e+00 (|T^5| = 1.76286546e+00) ==> Suspected discontinuity encountered at r=3.26411392e+00 (|T^5| = 1.59197475e+00) ==> Suspected discontinuity encountered at r=3.27411392e+00 (|T^5| = 2.20106296e+00) ==> Suspected discontinuity encountered at r=3.28411392e+00 (|T^5| = 2.39981447e+00) ==> Suspected discontinuity encountered at r=3.29411392e+00 (|T^5| = 1.45139679e+00) ==> Suspected discontinuity encountered at r=3.31411392e+00 (|T^5| = 1.46182296e+00) ==> Suspected discontinuity encountered at r=3.34411392e+00 (|T^5| = 1.60544834e+00) ==> Suspected discontinuity encountered at r=3.37411392e+00 (|T^5| = 1.66853072e+00) ==> Suspected discontinuity encountered at r=3.40411392e+00 (|T^5| = 1.81963236e+00) ==> Suspected discontinuity encountered at r=3.42411392e+00 (|T^5| = 1.05363538e+00) ==> Suspected discontinuity encountered at r=3.43411392e+00 (|T^5| = 2.12069530e+00) ==> Suspected discontinuity encountered at r=3.44411392e+00 (|T^5| = 1.11561045e+00) ==> Suspected discontinuity encountered at r=3.46411392e+00 (|T^5| = 2.09203164e+00) ==> Suspected discontinuity encountered at r=3.47411392e+00 (|T^5| = 2.53183864e+00) ==> Suspected discontinuity encountered at r=3.48411392e+00 (|T^5| = 5.14051931e+00) ==> Suspected discontinuity encountered at r=3.49411392e+00 (|T^5| = 6.95468165e+00) ==> Suspected discontinuity encountered at r=3.50411392e+00 (|T^5| = 1.96272845e+00) ==> Suspected discontinuity encountered at r=3.51411392e+00 (|T^5| = 5.35756681e+00) ==> Suspected discontinuity encountered at r=3.52411392e+00 (|T^5| = 7.98615955e+00) ==> Suspected discontinuity encountered at r=3.53411392e+00 (|T^5| = 7.44970877e+00) ==> Suspected discontinuity encountered at r=3.54411392e+00 (|T^5| = 5.45790727e+00) ==> Suspected discontinuity encountered at r=3.55411392e+00 (|T^5| = 2.74816039e+00) ==> Suspected discontinuity encountered at r=3.56411392e+00 (|T^5| = 2.03702860e+00) ==> Suspected discontinuity encountered at r=3.57411392e+00 (|T^5| = 2.53463717e+00) ==> Suspected discontinuity encountered at r=3.58411392e+00 (|T^5| = 1.74272373e+00) ==> Suspected discontinuity encountered at r=3.60411392e+00 (|T^5| = 1.18574862e+00) ==> Suspected discontinuity encountered at r=3.62411392e+00 (|T^5| = 3.42310275e+00) ==> Suspected discontinuity encountered at r=3.63411392e+00 (|T^5| = 3.26493166e+00) ==> Suspected discontinuity encountered at r=3.64411392e+00 (|T^5| = 1.03991920e+00) ==> Suspected discontinuity encountered at r=3.65411392e+00 (|T^5| = 8.26422077e+00) ==> Suspected discontinuity encountered at r=3.66411392e+00 (|T^5| = 1.26908574e+01) ==> Suspected discontinuity encountered at r=3.67411392e+00 (|T^5| = 8.66172360e+00) ==> Suspected discontinuity encountered at r=3.68411392e+00 (|T^5| = 5.08748053e+00) ==> Suspected discontinuity encountered at r=3.69411392e+00 (|T^5| = 2.01898354e+01) ==> Suspected discontinuity encountered at r=3.70411392e+00 (|T^5| = 5.26125070e+01) ==> Suspected discontinuity encountered at r=3.71411392e+00 (|T^5| = 6.54590615e+01) ==> Suspected discontinuity encountered at r=3.72411392e+00 (|T^5| = 3.94873665e+01) ==> Suspected discontinuity encountered at r=3.73411392e+00 (|T^5| = 9.35790969e+00) ==> Suspected discontinuity encountered at r=3.74411392e+00 (|T^5| = 4.58401950e+00) ==> Suspected discontinuity encountered at r=3.75411392e+00 (|T^5| = 1.84948212e+01) ==> Suspected discontinuity encountered at r=3.76411392e+00 (|T^5| = 2.83408210e+01) ==> Suspected discontinuity encountered at r=3.77411392e+00 (|T^5| = 1.93582294e+01) ==> Suspected discontinuity encountered at r=3.78411392e+00 (|T^5| = 4.96178536e+00) ==> Suspected discontinuity encountered at r=3.79411392e+00 (|T^5| = 3.68021859e+01) ==> Suspected discontinuity encountered at r=3.80411392e+00 (|T^5| = 1.89065804e+02) ==> Suspected discontinuity encountered at r=3.81411392e+00 (|T^5| = 3.88721582e+02) ==> Suspected discontinuity encountered at r=3.82411392e+00 (|T^5| = 3.98205590e+02) ==> Suspected discontinuity encountered at r=3.83411392e+00 (|T^5| = 1.99221711e+02) ==> Suspected discontinuity encountered at r=3.84411392e+00 (|T^5| = 3.25747325e+01) ==> Suspected discontinuity encountered at r=3.85411392e+00 (|T^5| = 1.02104927e+01) ==> Suspected discontinuity encountered at r=3.86411392e+00 (|T^5| = 6.25418012e+01) ==> Suspected discontinuity encountered at r=3.87411392e+00 (|T^5| = 2.11726629e+02) ==> Suspected discontinuity encountered at r=3.88411392e+00 (|T^5| = 3.06796301e+02) ==> Suspected discontinuity encountered at r=3.89411392e+00 (|T^5| = 2.03681888e+02) ==> Suspected discontinuity encountered at r=3.90411392e+00 (|T^5| = 5.13799468e+01) ==> Suspected discontinuity encountered at r=3.93411392e+00 (|T^5| = 1.06121972e+02) ==> Suspected discontinuity encountered at r=3.94411392e+00 (|T^5| = 4.37243460e+02) ==> Suspected discontinuity encountered at r=3.95411392e+00 (|T^5| = 6.75913671e+02) ==> Suspected discontinuity encountered at r=3.96411392e+00 (|T^5| = 4.64610540e+02) ==> Suspected discontinuity encountered at r=3.97411392e+00 (|T^5| = 1.19823099e+02) ==> Suspected discontinuity encountered at r=4.01411392e+00 (|T^5| = 7.25045362e+01) ==> Suspected discontinuity encountered at r=4.02411392e+00 (|T^5| = 2.60082592e+02) ==> Suspected discontinuity encountered at r=4.03411392e+00 (|T^5| = 3.00502377e+02) ==> Suspected discontinuity encountered at r=4.04411392e+00 (|T^5| = 6.48813899e+01) ==> Suspected discontinuity encountered at r=4.05411392e+00 (|T^5| = 9.39959864e+01) ==> Suspected discontinuity encountered at r=4.06411392e+00 (|T^5| = 4.59598208e+01) ==> Suspected discontinuity encountered at r=4.12411392e+00 (|T^5| = 3.50463738e+02) ==> Suspected discontinuity encountered at r=4.13411392e+00 (|T^5| = 1.76565591e+03) ==> Suspected discontinuity encountered at r=4.14411392e+00 (|T^5| = 3.55907323e+03) ==> Suspected discontinuity encountered at r=4.15411392e+00 (|T^5| = 3.58651536e+03) ==> Suspected discontinuity encountered at r=4.16411392e+00 (|T^5| = 1.80329483e+03) ==> Suspected discontinuity encountered at r=4.17411392e+00 (|T^5| = 3.56390664e+02) ==> Suspected discontinuity encountered at r=4.18411392e+00 (|T^5| = 5.68357358e+00) ==> Suspected discontinuity encountered at r=4.19411392e+00 (|T^5| = 1.41376165e+00) ==> Suspected discontinuity encountered at r=4.28411392e+00 (|T^5| = 4.54758312e+01) ==> Suspected discontinuity encountered at r=4.29411392e+00 (|T^5| = 1.89298180e+02) ==> Suspected discontinuity encountered at r=4.30411392e+00 (|T^5| = 2.95855654e+02) ==> Suspected discontinuity encountered at r=4.31411392e+00 (|T^5| = 2.05758827e+02) ==> Suspected discontinuity encountered at r=4.32411392e+00 (|T^5| = 5.37305262e+01) ==> Suspected discontinuity encountered at r=4.68411392e+00 (|T^5| = 3.28279977e+00) ==> Suspected discontinuity encountered at r=4.69411392e+00 (|T^5| = 4.92358451e+00) ==> Suspected discontinuity encountered at r=4.70411392e+00 (|T^5| = 3.28302261e+00) ==> Suspected discontinuity encountered at r=4.87411392e+00 (|T^5| = 4.04657487e+00) ==> Suspected discontinuity encountered at r=4.88411392e+00 (|T^5| = 1.32515513e+01) ==> Suspected discontinuity encountered at r=4.89411392e+00 (|T^5| = 1.84094338e+01) ==> Suspected discontinuity encountered at r=4.90411392e+00 (|T^5| = 1.17838000e+01) ==> Suspected discontinuity encountered at r=4.91411392e+00 (|T^5| = 2.87289038e+00) ==> Suspected discontinuity encountered at r=5.07411392e+00 (|T^5| = 2.78609572e+00) ==> Suspected discontinuity encountered at r=5.08411392e+00 (|T^5| = 4.17938825e+00) ==> Suspected discontinuity encountered at r=5.09411392e+00 (|T^5| = 2.78641778e+00) ==> Suspected discontinuity encountered at r=5.37411392e+00 (|T^5| = 1.61785046e+00) ==> Suspected discontinuity encountered at r=5.38411392e+00 (|T^5| = 6.45808614e+00) ==> Suspected discontinuity encountered at r=5.39411392e+00 (|T^5| = 9.65356860e+00) ==> Suspected discontinuity encountered at r=5.40411392e+00 (|T^5| = 6.39093954e+00) ==> Suspected discontinuity encountered at r=5.41411392e+00 (|T^5| = 1.56415765e+00) ----------------------------------------------------------------------------------------------------- DIMER H--He ----------------------------------------------------------------------------------------------------- E n e r g y E x t r e m a distance energy #iter #fn-calls warnings Equilibrium separation: 1.66973929e+00 -2.17312658e-01 30 60 none Cutoff separation: 3.00214146e+00 -2.22325255e-16 26 28 none C u t o f f S m o o t h n e s s value continuous Energy -2.22325255e-16 yes 1st-Deriv 6.17648873e-07 yes 2nd-Deriv 5.76847244e-04 ** NO ** 3rd-Deriv 2.69370637e-01 ** NO ** C o n t i n u i t y Checking continuity for r = [0.83487,3.00214] at 217 points (Delta r = 0.01000) Local edge detection based on a normalized 5th-order local difference formula T^5 is used to determine the presence of discontinuities. The tolernace is |T^5|>1.00000. (For details see Anne Gelb and Eitan Tadmor, J. Sci. Comp., 28:279-306, 2006.) Checking Energy (energy vs. distance in file "dimer-energy-HHe.dat") ... No discontinuities found. Checking 1st-Deriv (1st-deriv vs. distance in file "dimer-energy-deriv-1-HHe.dat") ... No discontinuities found. Checking 2nd-Deriv (2nd-deriv vs. distance in file "dimer-energy-deriv-2-HHe.dat") ... No discontinuities found. Checking 3rd-Deriv (3rd-deriv vs. distance in file "dimer-energy-deriv-3-HHe.dat") ==> Suspected discontinuity encountered at r=1.97486965e+00 (|T^5| = 3.25519775e+00) ==> Suspected discontinuity encountered at r=1.98486965e+00 (|T^5| = 1.30242282e+01) ==> Suspected discontinuity encountered at r=1.99486965e+00 (|T^5| = 1.95374265e+01) ==> Suspected discontinuity encountered at r=2.00486965e+00 (|T^5| = 1.30252404e+01) ==> Suspected discontinuity encountered at r=2.01486965e+00 (|T^5| = 3.25639501e+00) ==> Suspected discontinuity encountered at r=2.48486965e+00 (|T^5| = 1.59990891e+00) ==> Suspected discontinuity encountered at r=2.49486965e+00 (|T^5| = 2.39993239e+00) ==> Suspected discontinuity encountered at r=2.50486965e+00 (|T^5| = 1.59996154e+00) ----------------------------------------------------------------------------------------------------- DIMER H--W ----------------------------------------------------------------------------------------------------- E n e r g y E x t r e m a distance energy #iter #fn-calls warnings Equilibrium separation: 2.03397572e+00 -1.16403662e+01 31 62 none Cutoff separation: 4.90329759e+00 -2.16324078e-16 27 29 none C u t o f f S m o o t h n e s s value continuous Energy -2.16324078e-16 yes 1st-Deriv 1.35632241e-04 ** NO ** 2nd-Deriv 8.39962345e-02 ** NO ** 3rd-Deriv 2.75630523e+01 ** NO ** C o n t i n u i t y Checking continuity for r = [1.01699,4.90330] at 389 points (Delta r = 0.01000) Local edge detection based on a normalized 5th-order local difference formula T^5 is used to determine the presence of discontinuities. The tolernace is |T^5|>1.00000. (For details see Anne Gelb and Eitan Tadmor, J. Sci. Comp., 28:279-306, 2006.) Checking Energy (energy vs. distance in file "dimer-energy-HW.dat") ... No discontinuities found. Checking 1st-Deriv (1st-deriv vs. distance in file "dimer-energy-deriv-1-HW.dat") ... No discontinuities found. Checking 2nd-Deriv (2nd-deriv vs. distance in file "dimer-energy-deriv-2-HW.dat") ==> Suspected discontinuity encountered at r=1.97698786e+00 (|T^5| = 7.71012550e+00) ==> Suspected discontinuity encountered at r=1.98698786e+00 (|T^5| = 3.05673551e+01) ==> Suspected discontinuity encountered at r=1.99698786e+00 (|T^5| = 4.54678530e+01) ==> Suspected discontinuity encountered at r=2.00698786e+00 (|T^5| = 3.00609759e+01) ==> Suspected discontinuity encountered at r=2.01698786e+00 (|T^5| = 7.45003002e+00) ==> Suspected discontinuity encountered at r=3.81698786e+00 (|T^5| = 1.13603411e+00) ==> Suspected discontinuity encountered at r=3.86698786e+00 (|T^5| = 1.79911382e+00) ==> Suspected discontinuity encountered at r=3.87698786e+00 (|T^5| = 3.16052975e+00) ==> Suspected discontinuity encountered at r=3.88698786e+00 (|T^5| = 2.27379498e+00) ==> Suspected discontinuity encountered at r=3.94698786e+00 (|T^5| = 2.33103261e+00) ==> Suspected discontinuity encountered at r=3.95698786e+00 (|T^5| = 5.11516768e+00) ==> Suspected discontinuity encountered at r=3.96698786e+00 (|T^5| = 4.58149821e+00) ==> Suspected discontinuity encountered at r=3.97698786e+00 (|T^5| = 1.46395894e+00) ==> Suspected discontinuity encountered at r=4.04698786e+00 (|T^5| = 1.91567671e+00) ==> Suspected discontinuity encountered at r=4.12698786e+00 (|T^5| = 1.27338300e+00) ==> Suspected discontinuity encountered at r=4.13698786e+00 (|T^5| = 5.23469695e+00) ==> Suspected discontinuity encountered at r=4.14698786e+00 (|T^5| = 8.07440472e+00) ==> Suspected discontinuity encountered at r=4.15698786e+00 (|T^5| = 5.53914357e+00) ==> Suspected discontinuity encountered at r=4.16698786e+00 (|T^5| = 1.42601160e+00) ==> Suspected discontinuity encountered at r=4.28698786e+00 (|T^5| = 1.44124874e+00) ==> Suspected discontinuity encountered at r=4.30698786e+00 (|T^5| = 1.85858448e+00) ==> Suspected discontinuity encountered at r=4.31698786e+00 (|T^5| = 2.51821684e+00) Checking 3rd-Deriv (3rd-deriv vs. distance in file "dimer-energy-deriv-3-HW.dat") ==> Suspected discontinuity encountered at r=1.96698786e+00 (|T^5| = 2.80327882e+00) ==> Suspected discontinuity encountered at r=1.97698786e+00 (|T^5| = 1.73206062e+01) ==> Suspected discontinuity encountered at r=1.98698786e+00 (|T^5| = 1.03389235e+02) ==> Suspected discontinuity encountered at r=1.99698786e+00 (|T^5| = 1.85060613e+02) ==> Suspected discontinuity encountered at r=2.00698786e+00 (|T^5| = 1.50087829e+02) ==> Suspected discontinuity encountered at r=2.01698786e+00 (|T^5| = 5.51779813e+01) ==> Suspected discontinuity encountered at r=2.02698786e+00 (|T^5| = 6.90256587e+00) ==> Suspected discontinuity encountered at r=2.47698786e+00 (|T^5| = 3.37676620e+01) ==> Suspected discontinuity encountered at r=2.48698786e+00 (|T^5| = 1.32304571e+02) ==> Suspected discontinuity encountered at r=2.49698786e+00 (|T^5| = 1.95780637e+02) ==> Suspected discontinuity encountered at r=2.50698786e+00 (|T^5| = 1.28881768e+02) ==> Suspected discontinuity encountered at r=2.51698786e+00 (|T^5| = 3.18447899e+01) ==> Suspected discontinuity encountered at r=2.57698786e+00 (|T^5| = 1.10462618e+00) ==> Suspected discontinuity encountered at r=2.58698786e+00 (|T^5| = 1.09517740e+00) ==> Suspected discontinuity encountered at r=2.59698786e+00 (|T^5| = 1.03310055e+00) ==> Suspected discontinuity encountered at r=2.60698786e+00 (|T^5| = 1.03740250e+00) ==> Suspected discontinuity encountered at r=2.61698786e+00 (|T^5| = 2.49555385e+00) ==> Suspected discontinuity encountered at r=2.62698786e+00 (|T^5| = 2.96034097e+00) ==> Suspected discontinuity encountered at r=2.63698786e+00 (|T^5| = 5.60513186e+00) ==> Suspected discontinuity encountered at r=2.64698786e+00 (|T^5| = 1.64054969e+01) ==> Suspected discontinuity encountered at r=2.65698786e+00 (|T^5| = 1.37996577e+01) ==> Suspected discontinuity encountered at r=2.66698786e+00 (|T^5| = 4.10471115e+00) ==> Suspected discontinuity encountered at r=2.98698786e+00 (|T^5| = 1.01554209e+00) ==> Suspected discontinuity encountered at r=3.02698786e+00 (|T^5| = 1.17152281e+00) ==> Suspected discontinuity encountered at r=3.03698786e+00 (|T^5| = 1.25134987e+00) ==> Suspected discontinuity encountered at r=3.04698786e+00 (|T^5| = 1.27040150e+00) ==> Suspected discontinuity encountered at r=3.05698786e+00 (|T^5| = 1.28929691e+00) ==> Suspected discontinuity encountered at r=3.06698786e+00 (|T^5| = 1.28470586e+00) ==> Suspected discontinuity encountered at r=3.07698786e+00 (|T^5| = 8.47685790e+00) ==> Suspected discontinuity encountered at r=3.08698786e+00 (|T^5| = 2.98952539e+01) ==> Suspected discontinuity encountered at r=3.09698786e+00 (|T^5| = 4.38835550e+01) ==> Suspected discontinuity encountered at r=3.10698786e+00 (|T^5| = 2.94848808e+01) ==> Suspected discontinuity encountered at r=3.11698786e+00 (|T^5| = 8.35107396e+00) ==> Suspected discontinuity encountered at r=3.12698786e+00 (|T^5| = 1.40928032e+00) ==> Suspected discontinuity encountered at r=3.13698786e+00 (|T^5| = 1.44309741e+00) ==> Suspected discontinuity encountered at r=3.14698786e+00 (|T^5| = 1.49353302e+00) ==> Suspected discontinuity encountered at r=3.15698786e+00 (|T^5| = 1.55266539e+00) ==> Suspected discontinuity encountered at r=3.16698786e+00 (|T^5| = 1.61314671e+00) ==> Suspected discontinuity encountered at r=3.17698786e+00 (|T^5| = 1.73232148e+00) ==> Suspected discontinuity encountered at r=3.18698786e+00 (|T^5| = 1.70939334e+00) ==> Suspected discontinuity encountered at r=3.19698786e+00 (|T^5| = 1.02840112e+00) ==> Suspected discontinuity encountered at r=3.21698786e+00 (|T^5| = 1.15014084e+00) ==> Suspected discontinuity encountered at r=3.22698786e+00 (|T^5| = 1.58759833e+00) ==> Suspected discontinuity encountered at r=3.23698786e+00 (|T^5| = 1.76330908e+00) ==> Suspected discontinuity encountered at r=3.24698786e+00 (|T^5| = 1.85309588e+00) ==> Suspected discontinuity encountered at r=3.25698786e+00 (|T^5| = 1.67881025e+00) ==> Suspected discontinuity encountered at r=3.27698786e+00 (|T^5| = 6.53468497e+00) ==> Suspected discontinuity encountered at r=3.28698786e+00 (|T^5| = 1.01796209e+01) ==> Suspected discontinuity encountered at r=3.29698786e+00 (|T^5| = 6.51393934e+00) ==> Suspected discontinuity encountered at r=3.31698786e+00 (|T^5| = 1.87286780e+00) ==> Suspected discontinuity encountered at r=3.32698786e+00 (|T^5| = 2.28763503e+00) ==> Suspected discontinuity encountered at r=3.33698786e+00 (|T^5| = 2.17994367e+00) ==> Suspected discontinuity encountered at r=3.35698786e+00 (|T^5| = 1.53319130e+00) ==> Suspected discontinuity encountered at r=3.36698786e+00 (|T^5| = 3.26399076e+00) ==> Suspected discontinuity encountered at r=3.37698786e+00 (|T^5| = 2.17233124e+00) ==> Suspected discontinuity encountered at r=3.38698786e+00 (|T^5| = 1.25516175e+00) ==> Suspected discontinuity encountered at r=3.39698786e+00 (|T^5| = 3.81057795e+00) ==> Suspected discontinuity encountered at r=3.40698786e+00 (|T^5| = 3.06872855e+00) ==> Suspected discontinuity encountered at r=3.45698786e+00 (|T^5| = 1.19867658e+00) ==> Suspected discontinuity encountered at r=3.46698786e+00 (|T^5| = 2.08268438e+00) ==> Suspected discontinuity encountered at r=3.47698786e+00 (|T^5| = 1.79796480e+00) ==> Suspected discontinuity encountered at r=3.48698786e+00 (|T^5| = 1.44745364e+00) ==> Suspected discontinuity encountered at r=3.51698786e+00 (|T^5| = 1.28613572e+00) ==> Suspected discontinuity encountered at r=3.52698786e+00 (|T^5| = 1.82381481e+00) ==> Suspected discontinuity encountered at r=3.53698786e+00 (|T^5| = 2.34538576e+00) ==> Suspected discontinuity encountered at r=3.54698786e+00 (|T^5| = 1.60090032e+00) ==> Suspected discontinuity encountered at r=3.56698786e+00 (|T^5| = 1.63339529e+00) ==> Suspected discontinuity encountered at r=3.57698786e+00 (|T^5| = 2.50722891e+00) ==> Suspected discontinuity encountered at r=3.58698786e+00 (|T^5| = 1.70357722e+00) ==> Suspected discontinuity encountered at r=3.59698786e+00 (|T^5| = 1.92221657e+00) ==> Suspected discontinuity encountered at r=3.60698786e+00 (|T^5| = 6.12704277e+00) ==> Suspected discontinuity encountered at r=3.61698786e+00 (|T^5| = 9.42027034e+00) ==> Suspected discontinuity encountered at r=3.62698786e+00 (|T^5| = 6.43080403e+00) ==> Suspected discontinuity encountered at r=3.63698786e+00 (|T^5| = 8.73844777e+00) ==> Suspected discontinuity encountered at r=3.64698786e+00 (|T^5| = 3.74943059e+01) ==> Suspected discontinuity encountered at r=3.65698786e+00 (|T^5| = 7.92546874e+01) ==> Suspected discontinuity encountered at r=3.66698786e+00 (|T^5| = 8.37321103e+01) ==> Suspected discontinuity encountered at r=3.67698786e+00 (|T^5| = 4.42226239e+01) ==> Suspected discontinuity encountered at r=3.68698786e+00 (|T^5| = 2.87178776e+00) ==> Suspected discontinuity encountered at r=3.69698786e+00 (|T^5| = 2.65474513e+01) ==> Suspected discontinuity encountered at r=3.70698786e+00 (|T^5| = 4.08355258e+01) ==> Suspected discontinuity encountered at r=3.71698786e+00 (|T^5| = 2.79263337e+01) ==> Suspected discontinuity encountered at r=3.72698786e+00 (|T^5| = 7.17587598e+00) ==> Suspected discontinuity encountered at r=3.73698786e+00 (|T^5| = 3.51949445e+00) ==> Suspected discontinuity encountered at r=3.74698786e+00 (|T^5| = 1.30742722e+01) ==> Suspected discontinuity encountered at r=3.75698786e+00 (|T^5| = 1.67266064e+01) ==> Suspected discontinuity encountered at r=3.76698786e+00 (|T^5| = 6.80079510e+00) ==> Suspected discontinuity encountered at r=3.77698786e+00 (|T^5| = 1.82687907e+00) ==> Suspected discontinuity encountered at r=3.78698786e+00 (|T^5| = 1.45426002e+00) ==> Suspected discontinuity encountered at r=3.79698786e+00 (|T^5| = 5.64925824e+00) ==> Suspected discontinuity encountered at r=3.80698786e+00 (|T^5| = 2.30771466e+01) ==> Suspected discontinuity encountered at r=3.81698786e+00 (|T^5| = 3.54448425e+01) ==> Suspected discontinuity encountered at r=3.82698786e+00 (|T^5| = 2.42060332e+01) ==> Suspected discontinuity encountered at r=3.83698786e+00 (|T^5| = 6.19383115e+00) ==> Suspected discontinuity encountered at r=3.85698786e+00 (|T^5| = 1.33343883e+02) ==> Suspected discontinuity encountered at r=3.86698786e+00 (|T^5| = 7.09365612e+02) ==> Suspected discontinuity encountered at r=3.87698786e+00 (|T^5| = 1.50892243e+03) ==> Suspected discontinuity encountered at r=3.88698786e+00 (|T^5| = 1.60435857e+03) ==> Suspected discontinuity encountered at r=3.89698786e+00 (|T^5| = 8.52696049e+02) ==> Suspected discontinuity encountered at r=3.90698786e+00 (|T^5| = 1.81239497e+02) ==> Suspected discontinuity encountered at r=3.92698786e+00 (|T^5| = 3.85054408e+01) ==> Suspected discontinuity encountered at r=3.93698786e+00 (|T^5| = 2.98704874e+02) ==> Suspected discontinuity encountered at r=3.94698786e+00 (|T^5| = 8.22659072e+02) ==> Suspected discontinuity encountered at r=3.95698786e+00 (|T^5| = 1.06178272e+03) ==> Suspected discontinuity encountered at r=3.96698786e+00 (|T^5| = 6.57895200e+02) ==> Suspected discontinuity encountered at r=3.97698786e+00 (|T^5| = 1.58578311e+02) ==> Suspected discontinuity encountered at r=4.01698786e+00 (|T^5| = 1.30599085e+02) ==> Suspected discontinuity encountered at r=4.02698786e+00 (|T^5| = 5.51333867e+02) ==> Suspected discontinuity encountered at r=4.03698786e+00 (|T^5| = 8.84328436e+02) ==> Suspected discontinuity encountered at r=4.04698786e+00 (|T^5| = 6.49815509e+02) ==> Suspected discontinuity encountered at r=4.05698786e+00 (|T^5| = 1.98907480e+02) ==> Suspected discontinuity encountered at r=4.06698786e+00 (|T^5| = 1.26875343e+01) ==> Suspected discontinuity encountered at r=4.12698786e+00 (|T^5| = 1.46579923e+01) ==> Suspected discontinuity encountered at r=4.13698786e+00 (|T^5| = 5.97826453e+01) ==> Suspected discontinuity encountered at r=4.14698786e+00 (|T^5| = 9.14339724e+01) ==> Suspected discontinuity encountered at r=4.15698786e+00 (|T^5| = 6.21824686e+01) ==> Suspected discontinuity encountered at r=4.16698786e+00 (|T^5| = 1.58733514e+01) ==> Suspected discontinuity encountered at r=4.27698786e+00 (|T^5| = 4.32455248e+02) ==> Suspected discontinuity encountered at r=4.28698786e+00 (|T^5| = 2.20714304e+03) ==> Suspected discontinuity encountered at r=4.29698786e+00 (|T^5| = 4.51143992e+03) ==> Suspected discontinuity encountered at r=4.30698786e+00 (|T^5| = 4.61684148e+03) ==> Suspected discontinuity encountered at r=4.31698786e+00 (|T^5| = 2.36569968e+03) ==> Suspected discontinuity encountered at r=4.32698786e+00 (|T^5| = 4.85615143e+02) ----------------------------------------------------------------------------------------------------- DIMER He--W ----------------------------------------------------------------------------------------------------- E n e r g y E x t r e m a distance energy #iter #fn-calls warnings Equilibrium separation: 4.15000000e+00 0.00000000e+00 34 95 none Cutoff separation: 3.49964981e+00 8.09588232e-16 29 31 none C u t o f f S m o o t h n e s s value continuous Energy 8.09588232e-16 yes 1st-Deriv -3.47098798e-29 yes 2nd-Deriv 8.07793567e-28 yes 3rd-Deriv -5.16987883e-26 yes C o n t i n u i t y Checking continuity for r = [2.07500,3.49965] at 143 points (Delta r = 0.01000) Local edge detection based on a normalized 5th-order local difference formula T^5 is used to determine the presence of discontinuities. The tolernace is |T^5|>1.00000. (For details see Anne Gelb and Eitan Tadmor, J. Sci. Comp., 28:279-306, 2006.) Checking Energy (energy vs. distance in file "dimer-energy-HeW.dat") ... No discontinuities found. Checking 1st-Deriv (1st-deriv vs. distance in file "dimer-energy-deriv-1-HeW.dat") ... No discontinuities found. Checking 2nd-Deriv (2nd-deriv vs. distance in file "dimer-energy-deriv-2-HeW.dat") ... No discontinuities found. Checking 3rd-Deriv (3rd-deriv vs. distance in file "dimer-energy-deriv-3-HeW.dat") ==> Suspected discontinuity encountered at r=2.17500000e+00 (|T^5| = 3.71006707e+00) ==> Suspected discontinuity encountered at r=2.18500000e+00 (|T^5| = 1.48499055e+01) ==> Suspected discontinuity encountered at r=2.19500000e+00 (|T^5| = 2.22743800e+01) ==> Suspected discontinuity encountered at r=2.20500000e+00 (|T^5| = 1.48469039e+01) ==> Suspected discontinuity encountered at r=2.21500000e+00 (|T^5| = 3.70750873e+00) ===================================================================================================== SUMMARY of Model Continuity Results Across All Dimers: continuous ----------------------------- Energy yes 1st-Deriv ** NO ** 2nd-Deriv ** NO ** 3rd-Deriv ** NO ** ===================================================================================================== Continuity must be C^1 or higher to pass this verification check. Grade: F Comment: The model is C^0 continuous. This means that the model has continuous energy, but a discontinuous first derivative.