!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!! !!!!! !!!!! VERIFICATION CHECK: vc-permutation-symmetry !!!!! !!!!! !!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Description: Check whether a model is invariant with respect to atom permutations that preserve species, i.e. swapping any two atoms with the same species must not change the energy or forces. This must be true for all models. The check is performed for a randomly distorted non-periodic diamond cube base structure. Separate configurations are tested for each species supported by the model, as well as one containing a random distribution of all species. The energy and forces of each configuration are compared with one where each atom is randomly swapped with another atom of the same species. The energy and forces must remain unchanged. The verification check will pass if these conditions are satisfied for all configurations that the model is able to compute. Configurations used for testing are provided as auxiliary files. Author: Ellad Tadmor ------------------------------------------------------------------------------------------------------------------------- Results for KIM Model : Morse_Shifted_GirifalcoWeizer_1959MedCutoff_Al__MO_279544746097_002 Supported species : Al random seed = 13 lattice constant (orig) = 5.000 perturbation amplitude = 0.500 number unit cells per side = 2 ------------------------------------------------------------------------------------------------------------------------- MONOATOMIC STRUCTURE -- Species = Al (Configuration in file "config-Al.xyz") ------------------------------------------------------------------------------------------------------------------------- Permutation list = [3, 51, 32, 0, 62, 27, 31, 46, 48, 43, 53, 49, 41, 17, 52, 60, 37, 13, 45, 34, 55, 24, 30, 63, 21, 42, 54, 5, 33, 59, 22, 6, 2, 28, 19, 50, 56, 16, 47, 61, 58, 12, 25, 9, 57, 18, 7, 38, 8, 11, 35, 1, 14, 10, 26, 20, 36, 44, 40, 29, 15, 39, 4, 23] Energy requirement: V(r_(PI 1),...,r_(PI N)) = V(r_1,...,r_N), where r_i is the position of atom i, V is the potential energy, and PI is the permutation operator (which randomly swaps all atoms with other atoms of the same species). V(r_(PI 1),...,r_(PI N)) = 133.70945454063573 V(r_1,...,r_N) = 133.7094545406357 Forces requirement: f_(PI i)(r_(PI 1),...,r_(PI N)) = f_i(r_1,...,r_N), where r_i is the position of atom i, f_i is the force on atom i and PI is the permutation operator (which randomly swaps all atoms with other atoms of the same species). i f_(PI i)(r_(PI 1),...,r_(PI N)) f_i(r_1,...,r_N) ------------------------------------------------------------------------------------------------------------------------- 0 -5.83476919e+00 -2.71293664e+00 -7.56678604e-01 | -5.83476919e+00 -2.71293664e+00 -7.56678604e-01 1 -7.66925568e+00 -1.10776593e+01 4.15981740e+00 | -7.66925568e+00 -1.10776593e+01 4.15981740e+00 2 -2.40359090e+00 8.79295258e+00 -1.96684765e+01 | -2.40359090e+00 8.79295258e+00 -1.96684765e+01 3 1.05283577e+01 1.24866371e+01 5.24095710e+00 | 1.05283577e+01 1.24866371e+01 5.24095710e+00 4 -4.42588196e+00 -1.12198742e+01 5.96417765e+00 | -4.42588196e+00 -1.12198742e+01 5.96417765e+00 5 4.65643221e+00 8.50322421e-01 4.48737632e+00 | 4.65643221e+00 8.50322421e-01 4.48737632e+00 6 -1.31933711e+00 -8.41686477e-01 -1.74298960e+00 | -1.31933711e+00 -8.41686477e-01 -1.74298960e+00 7 -1.47525436e+01 3.44674257e+01 -2.10059993e+01 | -1.47525436e+01 3.44674257e+01 -2.10059993e+01 8 -5.51232052e+00 -4.38872979e+00 -7.74969833e+00 | -5.51232052e+00 -4.38872979e+00 -7.74969833e+00 9 6.04202838e+00 1.25268731e+00 3.30380870e+00 | 6.04202838e+00 1.25268731e+00 3.30380870e+00 10 -4.72987375e+00 -4.06256255e+00 -7.74040503e+00 | -4.72987375e+00 -4.06256255e+00 -7.74040503e+00 11 7.84743566e+00 5.66557444e+00 5.47966430e+00 | 7.84743566e+00 5.66557444e+00 5.47966430e+00 12 -3.73159056e+00 -1.25829872e+01 -1.11051339e+00 | -3.73159056e+00 -1.25829872e+01 -1.11051339e+00 13 1.18994156e+01 4.78942614e+00 6.00088223e-01 | 1.18994156e+01 4.78942614e+00 6.00088223e-01 14 -4.19465707e+00 1.02167947e+01 3.86344539e-01 | -4.19465707e+00 1.02167947e+01 3.86344539e-01 15 1.82839883e+00 -1.14275247e+00 5.58331957e-02 | 1.82839883e+00 -1.14275247e+00 5.58331957e-02 16 -5.31417741e+00 -4.17112061e+00 -4.13837610e+00 | -5.31417741e+00 -4.17112061e+00 -4.13837610e+00 17 3.84081600e+00 -1.12253139e+00 6.62656022e+00 | 3.84081600e+00 -1.12253139e+00 6.62656022e+00 18 -1.04145783e+00 1.70513283e+00 -5.16583305e+00 | -1.04145783e+00 1.70513283e+00 -5.16583305e+00 19 -2.55949526e+01 1.89705233e+01 -1.62612706e+01 | -2.55949526e+01 1.89705233e+01 -1.62612706e+01 20 -9.27810721e+00 -6.61081275e+00 -9.99442146e+00 | -9.27810721e+00 -6.61081275e+00 -9.99442146e+00 21 1.47608213e+01 9.48412987e+00 6.48417577e+00 | 1.47608213e+01 9.48412987e+00 6.48417577e+00 22 -5.01308726e+00 -6.34730650e+00 -1.11588512e+01 | -5.01308726e+00 -6.34730650e+00 -1.11588512e+01 23 3.91845231e+00 8.15293636e+00 1.13542306e+01 | 3.91845231e+00 8.15293636e+00 1.13542306e+01 24 -4.91028048e+00 -2.20838684e+00 -4.00527832e+00 | -4.91028048e+00 -2.20838684e+00 -4.00527832e+00 25 -6.67984665e+00 1.14627873e+01 -5.62727399e+00 | -6.67984665e+00 1.14627873e+01 -5.62727399e+00 26 -5.23661566e+00 1.43963229e+00 -2.74706865e+00 | -5.23661566e+00 1.43963229e+00 -2.74706865e+00 27 4.77132414e+00 -6.43809120e-01 2.57877319e+00 | 4.77132414e+00 -6.43809120e-01 2.57877319e+00 28 9.19206887e+00 -1.20949312e+01 1.00446166e+01 | 9.19206887e+00 -1.20949312e+01 1.00446166e+01 29 3.32294147e+00 -3.88417159e-01 1.27574580e+00 | 3.32294147e+00 -3.88417159e-01 1.27574580e+00 30 2.63236818e+01 -1.56647234e+01 1.66130978e+01 | 2.63236818e+01 -1.56647234e+01 1.66130978e+01 31 -1.45944946e+01 6.78285239e+00 -7.45864786e+00 | -1.45944946e+01 6.78285239e+00 -7.45864786e+00 32 -6.30409010e+00 -6.87847048e+00 -6.13044274e+00 | -6.30409010e+00 -6.87847048e+00 -6.13044274e+00 33 1.36724322e+00 8.12886282e+00 7.66054826e-02 | 1.36724322e+00 8.12886282e+00 7.66054826e-02 34 1.74667127e+01 -3.85006111e+01 2.05217501e+01 | 1.74667127e+01 -3.85006111e+01 2.05217501e+01 35 -1.02686909e+01 -1.63949295e+01 1.57603367e+01 | -1.02686909e+01 -1.63949295e+01 1.57603367e+01 36 2.84258970e+00 -3.23010125e+00 6.31680782e+00 | 2.84258970e+00 -3.23010125e+00 6.31680782e+00 37 -5.03429391e+00 -5.46663579e+00 4.76603469e+00 | -5.03429391e+00 -5.46663579e+00 4.76603469e+00 38 -7.96711112e+00 -6.19190112e+00 -7.74519134e+00 | -7.96711112e+00 -6.19190112e+00 -7.74519134e+00 39 -5.35143012e+00 -6.30808473e+00 1.45978716e+01 | -5.35143012e+00 -6.30808473e+00 1.45978716e+01 40 9.08681340e-01 -3.86400890e-01 2.52344687e+00 | 9.08681340e-01 -3.86400890e-01 2.52344687e+00 41 -1.22785704e+01 2.15260705e+01 -1.85640194e+01 | -1.22785704e+01 2.15260705e+01 -1.85640194e+01 42 -2.07206448e+01 -1.35346080e+01 -1.69108197e+01 | -2.07206448e+01 -1.35346080e+01 -1.69108197e+01 43 1.68354124e+01 1.86328943e+01 2.02073049e+01 | 1.68354124e+01 1.86328943e+01 2.02073049e+01 44 1.50724436e+01 -2.63225964e+01 1.19234511e+01 | 1.50724436e+01 -2.63225964e+01 1.19234511e+01 45 4.34021654e-01 1.34934013e+00 5.53645078e-01 | 4.34021654e-01 1.34934013e+00 5.53645078e-01 46 -2.57259656e+00 9.36654138e-01 -2.39281142e+00 | -2.57259656e+00 9.36654138e-01 -2.39281142e+00 47 6.61985192e-01 7.16969109e+00 3.64640756e+00 | 6.61985192e-01 7.16969109e+00 3.64640756e+00 48 -1.11708824e+01 -1.82169589e+00 -2.14087162e+00 | -1.11708824e+01 -1.82169589e+00 -2.14087162e+00 49 -2.25191163e-01 -7.71567596e-02 7.40976785e+00 | -2.25191163e-01 -7.71567596e-02 7.40976785e+00 50 -2.83641420e-01 -7.08602322e+00 -1.25474295e+01 | -2.83641420e-01 -7.08602322e+00 -1.25474295e+01 51 6.95976846e+00 1.19103113e+01 1.02824661e+01 | 6.95976846e+00 1.19103113e+01 1.02824661e+01 52 1.12020641e+01 1.30025909e+01 -8.39676727e+00 | 1.12020641e+01 1.30025909e+01 -8.39676727e+00 53 -2.60540249e+00 -2.22420081e+00 5.98946699e+00 | -2.60540249e+00 -2.22420081e+00 5.98946699e+00 54 -2.50497960e+00 -2.51707091e+00 -2.44558209e+00 | -2.50497960e+00 -2.51707091e+00 -2.44558209e+00 55 3.28239180e+00 1.71370291e+00 1.88488817e+00 | 3.28239180e+00 1.71370291e+00 1.88488817e+00 56 1.05507625e+01 1.35781513e+01 -1.18916884e+01 | 1.05507625e+01 1.35781513e+01 -1.18916884e+01 57 1.35127941e+01 -2.32860483e+01 -2.60338290e+01 | 1.35127941e+01 -2.32860483e+01 -2.60338290e+01 58 1.48364376e+01 -4.69390283e+00 9.78351482e+00 | 1.48364376e+01 -4.69390283e+00 9.78351482e+00 59 4.26413585e-01 -2.77495372e-01 -5.56986552e-01 | 4.26413585e-01 -2.77495372e-01 -5.56986552e-01 60 -6.19062769e+00 -4.88213086e+00 -6.80119481e+00 | -6.19062769e+00 -4.88213086e+00 -6.80119481e+00 61 8.38771729e+00 8.47151424e+00 6.29314047e+00 | 8.38771729e+00 8.47151424e+00 6.29314047e+00 62 -7.76522858e+00 2.41713227e+01 2.05327637e+01 | -7.76522858e+00 2.41713227e+01 2.05327637e+01 63 -1.99392199e-01 2.50370752e-01 1.16447860e+00 | -1.99392199e-01 2.50370752e-01 1.16447860e+00 ------------------------------------------------------------------------------------------------------------------------- PASS: Energies and forces are the same to within a relative error of 1e-08 ------------------------------------------------------------------------------------------------------------------------- ========================================================================================================================= To pass this verification check the model must be invariant with respect to perumutation symmetry for all configurations it was able to compute. Grade: P Comment: Model energy has permutation symmetry for all configurations the model was able to compute.