!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!! !!!!!
!!!!! 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_Mishin_NiAl__MO_101214310689_004
Supported species : Al Ni
-----------------------------------------------------------------------------------------------------
DIMER Al--Al
-----------------------------------------------------------------------------------------------------
E n e r g y E x t r e m a
distance energy #iter #fn-calls warnings
Equilibrium separation: 2.78970674e+00 -1.25550079e+00 32 64 none
Cutoff separation: 7.02439366e+00 -1.97380031e-02 28 30 none
C u t o f f S m o o t h n e s s
value continuous
Energy -1.97380031e-02 ** NO **
1st-Deriv 0.00000000e+00 yes
2nd-Deriv 0.00000000e+00 yes
3rd-Deriv 0.00000000e+00 yes
C o n t i n u i t y
Checking continuity for r = [1.39485,7.02439] at 563 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-AlAl.dat")
... No discontinuities found.
Checking 1st-Deriv (1st-deriv vs. distance in file "dimer-energy-deriv-1-AlAl.dat")
... No discontinuities found.
Checking 2nd-Deriv (2nd-deriv vs. distance in file "dimer-energy-deriv-2-AlAl.dat")
... No discontinuities found.
Checking 3rd-Deriv (3rd-deriv vs. distance in file "dimer-energy-deriv-3-AlAl.dat")
==> Suspected discontinuity encountered at r=1.56485337e+00 (|T^5| = 1.02446625e+00)
==> Suspected discontinuity encountered at r=1.59485337e+00 (|T^5| = 1.10443063e+00)
==> Suspected discontinuity encountered at r=1.62485337e+00 (|T^5| = 1.35452336e+00)
==> Suspected discontinuity encountered at r=1.63485337e+00 (|T^5| = 1.12271901e+00)
==> Suspected discontinuity encountered at r=1.68485337e+00 (|T^5| = 1.26396660e+00)
==> Suspected discontinuity encountered at r=1.69485337e+00 (|T^5| = 1.55243019e+00)
==> Suspected discontinuity encountered at r=1.70485337e+00 (|T^5| = 1.18391657e+00)
==> Suspected discontinuity encountered at r=1.74485337e+00 (|T^5| = 1.01388349e+00)
==> Suspected discontinuity encountered at r=1.75485337e+00 (|T^5| = 1.14356681e+00)
==> Suspected discontinuity encountered at r=1.76485337e+00 (|T^5| = 1.15974569e+00)
==> Suspected discontinuity encountered at r=1.77485337e+00 (|T^5| = 1.24092819e+00)
==> Suspected discontinuity encountered at r=1.89485337e+00 (|T^5| = 1.31724058e+00)
==> Suspected discontinuity encountered at r=1.90485337e+00 (|T^5| = 1.37406280e+00)
==> Suspected discontinuity encountered at r=4.90485337e+00 (|T^5| = 1.07596346e+00)
==> Suspected discontinuity encountered at r=5.69485337e+00 (|T^5| = 1.39281762e+00)
==> Suspected discontinuity encountered at r=5.70485337e+00 (|T^5| = 7.07778856e+00)
==> Suspected discontinuity encountered at r=5.71485337e+00 (|T^5| = 1.43827173e+01)
==> Suspected discontinuity encountered at r=5.72485337e+00 (|T^5| = 1.46185928e+01)
==> Suspected discontinuity encountered at r=5.73485337e+00 (|T^5| = 7.42865964e+00)
==> Suspected discontinuity encountered at r=5.74485337e+00 (|T^5| = 1.50110825e+00)
-----------------------------------------------------------------------------------------------------
DIMER Ni--Ni
-----------------------------------------------------------------------------------------------------
E n e r g y E x t r e m a
distance energy #iter #fn-calls warnings
Equilibrium separation: 2.28919415e+00 -1.87292575e+00 31 62 none
Cutoff separation: 7.81205003e+00 -4.50993003e-02 28 30 none
C u t o f f S m o o t h n e s s
value continuous
Energy -4.50993003e-02 ** NO **
1st-Deriv 0.00000000e+00 yes
2nd-Deriv 0.00000000e+00 yes
3rd-Deriv 0.00000000e+00 yes
C o n t i n u i t y
Checking continuity for r = [1.14460,7.81205] at 667 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-NiNi.dat")
... No discontinuities found.
Checking 1st-Deriv (1st-deriv vs. distance in file "dimer-energy-deriv-1-NiNi.dat")
... No discontinuities found.
Checking 2nd-Deriv (2nd-deriv vs. distance in file "dimer-energy-deriv-2-NiNi.dat")
... No discontinuities found.
Checking 3rd-Deriv (3rd-deriv vs. distance in file "dimer-energy-deriv-3-NiNi.dat")
==> Suspected discontinuity encountered at r=1.49459707e+00 (|T^5| = 1.68114854e+00)
==> Suspected discontinuity encountered at r=1.50459707e+00 (|T^5| = 2.70816245e+00)
==> Suspected discontinuity encountered at r=1.51459707e+00 (|T^5| = 4.16017593e+00)
==> Suspected discontinuity encountered at r=1.52459707e+00 (|T^5| = 4.44824080e+00)
==> Suspected discontinuity encountered at r=1.53459707e+00 (|T^5| = 2.12302503e+00)
==> Suspected discontinuity encountered at r=1.55459707e+00 (|T^5| = 1.09443824e+00)
==> Suspected discontinuity encountered at r=1.57459707e+00 (|T^5| = 2.37787214e+00)
==> Suspected discontinuity encountered at r=1.58459707e+00 (|T^5| = 3.69363470e+00)
==> Suspected discontinuity encountered at r=1.59459707e+00 (|T^5| = 4.58741438e+00)
==> Suspected discontinuity encountered at r=1.60459707e+00 (|T^5| = 3.54488211e+00)
==> Suspected discontinuity encountered at r=1.68459707e+00 (|T^5| = 1.03894078e+00)
==> Suspected discontinuity encountered at r=1.70459707e+00 (|T^5| = 3.51776646e+00)
==> Suspected discontinuity encountered at r=1.71459707e+00 (|T^5| = 2.65558187e+00)
==> Suspected discontinuity encountered at r=1.73459707e+00 (|T^5| = 2.30098414e+00)
==> Suspected discontinuity encountered at r=1.79459707e+00 (|T^5| = 2.10252708e+00)
==> Suspected discontinuity encountered at r=1.80459707e+00 (|T^5| = 1.96222557e+00)
==> Suspected discontinuity encountered at r=1.81459707e+00 (|T^5| = 1.13618355e+00)
==> Suspected discontinuity encountered at r=1.82459707e+00 (|T^5| = 3.75524874e+00)
==> Suspected discontinuity encountered at r=1.83459707e+00 (|T^5| = 3.23166978e+00)
==> Suspected discontinuity encountered at r=1.92459707e+00 (|T^5| = 1.50808782e+00)
==> Suspected discontinuity encountered at r=1.93459707e+00 (|T^5| = 3.70511795e+00)
==> Suspected discontinuity encountered at r=1.94459707e+00 (|T^5| = 4.43013887e+00)
==> Suspected discontinuity encountered at r=1.95459707e+00 (|T^5| = 3.57493665e+00)
==> Suspected discontinuity encountered at r=1.96459707e+00 (|T^5| = 1.75571690e+00)
==> Suspected discontinuity encountered at r=2.01459707e+00 (|T^5| = 1.10127125e+00)
==> Suspected discontinuity encountered at r=2.02459707e+00 (|T^5| = 1.32482056e+00)
==> Suspected discontinuity encountered at r=2.06459707e+00 (|T^5| = 1.40000333e+00)
==> Suspected discontinuity encountered at r=2.07459707e+00 (|T^5| = 2.79991345e+00)
==> Suspected discontinuity encountered at r=2.08459707e+00 (|T^5| = 2.88002662e+00)
==> Suspected discontinuity encountered at r=2.09459707e+00 (|T^5| = 1.30227701e+00)
==> Suspected discontinuity encountered at r=3.55459707e+00 (|T^5| = 1.02031287e+00)
==> Suspected discontinuity encountered at r=3.69459707e+00 (|T^5| = 1.62358326e+00)
==> Suspected discontinuity encountered at r=3.70459707e+00 (|T^5| = 3.18520231e+00)
==> Suspected discontinuity encountered at r=3.71459707e+00 (|T^5| = 2.98708005e+00)
==> Suspected discontinuity encountered at r=3.88459707e+00 (|T^5| = 1.07450897e+00)
==> Suspected discontinuity encountered at r=3.90459707e+00 (|T^5| = 1.08698567e+00)
==> Suspected discontinuity encountered at r=3.91459707e+00 (|T^5| = 3.72142699e+00)
==> Suspected discontinuity encountered at r=3.92459707e+00 (|T^5| = 5.77825772e+00)
==> Suspected discontinuity encountered at r=3.93459707e+00 (|T^5| = 3.99820680e+00)
==> Suspected discontinuity encountered at r=3.94459707e+00 (|T^5| = 1.05341067e+00)
==> Suspected discontinuity encountered at r=4.01459707e+00 (|T^5| = 1.11735548e+00)
==> Suspected discontinuity encountered at r=4.02459707e+00 (|T^5| = 2.03512207e+00)
==> Suspected discontinuity encountered at r=4.03459707e+00 (|T^5| = 1.82519326e+00)
==> Suspected discontinuity encountered at r=4.25459707e+00 (|T^5| = 1.28908699e+00)
==> Suspected discontinuity encountered at r=4.26459707e+00 (|T^5| = 2.71041380e+00)
==> Suspected discontinuity encountered at r=4.27459707e+00 (|T^5| = 2.83498177e+00)
==> Suspected discontinuity encountered at r=4.28459707e+00 (|T^5| = 1.47051158e+00)
==> Suspected discontinuity encountered at r=4.40459707e+00 (|T^5| = 2.25810224e+00)
==> Suspected discontinuity encountered at r=4.41459707e+00 (|T^5| = 3.51996646e+00)
==> Suspected discontinuity encountered at r=4.42459707e+00 (|T^5| = 2.44936544e+00)
-----------------------------------------------------------------------------------------------------
DIMER Al--Ni
-----------------------------------------------------------------------------------------------------
E n e r g y E x t r e m a
distance energy #iter #fn-calls warnings
Equilibrium separation: 2.46024622e+00 -1.66277523e+00 32 64 none
Cutoff separation: 7.14625833e+00 -3.24186517e-02 29 31 none
C u t o f f S m o o t h n e s s
value continuous
Energy -3.24186517e-02 ** NO **
1st-Deriv 0.00000000e+00 yes
2nd-Deriv 0.00000000e+00 yes
3rd-Deriv 0.00000000e+00 yes
C o n t i n u i t y
Checking continuity for r = [1.23012,7.14626] at 592 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-AlNi.dat")
... No discontinuities found.
Checking 1st-Deriv (1st-deriv vs. distance in file "dimer-energy-deriv-1-AlNi.dat")
... No discontinuities found.
Checking 2nd-Deriv (2nd-deriv vs. distance in file "dimer-energy-deriv-2-AlNi.dat")
... No discontinuities found.
Checking 3rd-Deriv (3rd-deriv vs. distance in file "dimer-energy-deriv-3-AlNi.dat")
==> Suspected discontinuity encountered at r=1.48012311e+00 (|T^5| = 1.00308857e+00)
==> Suspected discontinuity encountered at r=1.49012311e+00 (|T^5| = 2.82029976e+00)
==> Suspected discontinuity encountered at r=1.50012311e+00 (|T^5| = 3.05289048e+00)
==> Suspected discontinuity encountered at r=1.51012311e+00 (|T^5| = 2.53018246e+00)
==> Suspected discontinuity encountered at r=1.52012311e+00 (|T^5| = 2.43226790e+00)
==> Suspected discontinuity encountered at r=1.53012311e+00 (|T^5| = 1.07998808e+00)
==> Suspected discontinuity encountered at r=1.54012311e+00 (|T^5| = 1.12804303e+00)
==> Suspected discontinuity encountered at r=1.55012311e+00 (|T^5| = 2.03190559e+00)
==> Suspected discontinuity encountered at r=1.57012311e+00 (|T^5| = 1.49071124e+00)
==> Suspected discontinuity encountered at r=1.58012311e+00 (|T^5| = 2.73420084e+00)
==> Suspected discontinuity encountered at r=1.59012311e+00 (|T^5| = 3.03032663e+00)
==> Suspected discontinuity encountered at r=1.60012311e+00 (|T^5| = 2.48261979e+00)
==> Suspected discontinuity encountered at r=1.61012311e+00 (|T^5| = 1.89012059e+00)
==> Suspected discontinuity encountered at r=1.62012311e+00 (|T^5| = 1.85851710e+00)
==> Suspected discontinuity encountered at r=1.63012311e+00 (|T^5| = 1.38069602e+00)
==> Suspected discontinuity encountered at r=1.65012311e+00 (|T^5| = 1.15818031e+00)
==> Suspected discontinuity encountered at r=1.66012311e+00 (|T^5| = 3.06549499e+00)
==> Suspected discontinuity encountered at r=1.67012311e+00 (|T^5| = 3.81037103e+00)
==> Suspected discontinuity encountered at r=1.68012311e+00 (|T^5| = 2.17972555e+00)
==> Suspected discontinuity encountered at r=1.70012311e+00 (|T^5| = 2.43335225e+00)
==> Suspected discontinuity encountered at r=1.71012311e+00 (|T^5| = 2.60512406e+00)
==> Suspected discontinuity encountered at r=1.74012311e+00 (|T^5| = 1.94829187e+00)
==> Suspected discontinuity encountered at r=1.75012311e+00 (|T^5| = 3.15293059e+00)
==> Suspected discontinuity encountered at r=1.76012311e+00 (|T^5| = 3.94514325e+00)
==> Suspected discontinuity encountered at r=1.77012311e+00 (|T^5| = 3.63713249e+00)
==> Suspected discontinuity encountered at r=1.78012311e+00 (|T^5| = 2.94588577e+00)
==> Suspected discontinuity encountered at r=1.79012311e+00 (|T^5| = 1.50577333e+00)
==> Suspected discontinuity encountered at r=1.80012311e+00 (|T^5| = 1.00696444e+00)
==> Suspected discontinuity encountered at r=1.81012311e+00 (|T^5| = 2.49484995e+00)
==> Suspected discontinuity encountered at r=1.82012311e+00 (|T^5| = 1.40764147e+00)
==> Suspected discontinuity encountered at r=1.84012311e+00 (|T^5| = 1.99512924e+00)
==> Suspected discontinuity encountered at r=1.85012311e+00 (|T^5| = 1.38881364e+00)
==> Suspected discontinuity encountered at r=1.87012311e+00 (|T^5| = 1.39648946e+00)
==> Suspected discontinuity encountered at r=1.88012311e+00 (|T^5| = 2.40727441e+00)
==> Suspected discontinuity encountered at r=1.89012311e+00 (|T^5| = 2.39607445e+00)
==> Suspected discontinuity encountered at r=1.90012311e+00 (|T^5| = 1.66924800e+00)
==> Suspected discontinuity encountered at r=1.93012311e+00 (|T^5| = 1.86162526e+00)
==> Suspected discontinuity encountered at r=1.94012311e+00 (|T^5| = 3.00005550e+00)
==> Suspected discontinuity encountered at r=1.95012311e+00 (|T^5| = 2.58907773e+00)
==> Suspected discontinuity encountered at r=1.97012311e+00 (|T^5| = 1.12696994e+00)
==> Suspected discontinuity encountered at r=1.98012311e+00 (|T^5| = 2.17076959e+00)
==> Suspected discontinuity encountered at r=1.99012311e+00 (|T^5| = 1.68446920e+00)
==> Suspected discontinuity encountered at r=2.00012311e+00 (|T^5| = 1.01750141e+00)
==> Suspected discontinuity encountered at r=2.03012311e+00 (|T^5| = 2.02864335e+00)
==> Suspected discontinuity encountered at r=2.04012311e+00 (|T^5| = 1.22200291e+00)
==> Suspected discontinuity encountered at r=2.09012311e+00 (|T^5| = 1.40157725e+00)
==> Suspected discontinuity encountered at r=2.10012311e+00 (|T^5| = 2.38864273e+00)
==> Suspected discontinuity encountered at r=2.11012311e+00 (|T^5| = 2.93187431e+00)
==> Suspected discontinuity encountered at r=2.12012311e+00 (|T^5| = 2.81234270e+00)
==> Suspected discontinuity encountered at r=2.13012311e+00 (|T^5| = 2.19108104e+00)
==> Suspected discontinuity encountered at r=2.14012311e+00 (|T^5| = 1.08844183e+00)
==> Suspected discontinuity encountered at r=2.45012311e+00 (|T^5| = 1.04737844e+00)
==> Suspected discontinuity encountered at r=2.46012311e+00 (|T^5| = 1.02908145e+00)
==> Suspected discontinuity encountered at r=2.47012311e+00 (|T^5| = 1.05468467e+00)
==> Suspected discontinuity encountered at r=2.48012311e+00 (|T^5| = 1.05753604e+00)
==> Suspected discontinuity encountered at r=4.40012311e+00 (|T^5| = 3.75567239e+00)
==> Suspected discontinuity encountered at r=4.41012311e+00 (|T^5| = 7.59238787e+00)
==> Suspected discontinuity encountered at r=4.42012311e+00 (|T^5| = 7.73170474e+00)
==> Suspected discontinuity encountered at r=4.43012311e+00 (|T^5| = 3.87518122e+00)
==> Suspected discontinuity encountered at r=4.95012311e+00 (|T^5| = 1.80511746e+00)
==> Suspected discontinuity encountered at r=4.96012311e+00 (|T^5| = 1.86378276e+00)
==> Suspected discontinuity encountered at r=5.07012311e+00 (|T^5| = 1.13164787e+00)
==> Suspected discontinuity encountered at r=5.13012311e+00 (|T^5| = 1.95996987e+00)
==> Suspected discontinuity encountered at r=5.14012311e+00 (|T^5| = 1.95113315e+00)
=====================================================================================================
SUMMARY of Model Continuity Results Across All Dimers:
continuous
-----------------------------
Energy ** NO **
1st-Deriv yes
2nd-Deriv yes
3rd-Deriv ** NO **
=====================================================================================================
Continuity must be C^1 or higher to pass this verification check.
Grade: F
Comment: The model is C^-1 continuous. This means that the model has discontinuous energy.