!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!! !!!!! !!!!! 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 : Sim_LAMMPS_SMTBQ_SallesPolitanoAmzallag_2016_TiO__SM_349577644423_000 Supported species : O Ti random seed = 13 lattice constant (orig) = 5.000 perturbation amplitude = 0.500 number unit cells per side = 2 ------------------------------------------------------------------------------------------------------------------------- MONOATOMIC STRUCTURE -- Species = O (Configuration in file "config-O.xyz") ------------------------------------------------------------------------------------------------------------------------- Permutation list = [26, 59, 29, 61, 22, 44, 57, 45, 28, 12, 36, 33, 9, 53, 51, 63, 37, 35, 42, 50, 23, 31, 4, 20, 30, 40, 0, 41, 8, 2, 24, 21, 62, 11, 39, 17, 10, 16, 60, 34, 25, 27, 18, 47, 5, 7, 48, 43, 46, 58, 19, 14, 54, 13, 52, 56, 55, 6, 49, 1, 38, 3, 32, 15] 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)) = 154.495622064 V(r_1,...,r_N) = 154.495622064 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 -4.17008924e+00 -2.07872279e+00 -8.67375840e-01 | -4.17008924e+00 -2.07872279e+00 -8.67375840e-01 1 -1.20332366e-01 5.86587016e+00 -5.98811341e+00 | -1.20332366e-01 5.86587016e+00 -5.98811341e+00 2 1.70284425e+00 2.43866505e+00 1.25757946e+00 | 1.70284425e+00 2.43866505e+00 1.25757946e+00 3 -4.09993517e+00 -3.72398659e+00 -7.54481982e+00 | -4.09993517e+00 -3.72398659e+00 -7.54481982e+00 4 -2.45773085e-01 -7.74997248e+00 4.72782142e+00 | -2.45773085e-01 -7.74997248e+00 4.72782142e+00 5 3.55033907e+00 4.95309544e-01 2.19911422e+00 | 3.55033907e+00 4.95309544e-01 2.19911422e+00 6 -5.98678178e+00 -5.75750661e+00 -7.72816103e+00 | -5.98678178e+00 -5.75750661e+00 -7.72816103e+00 7 -3.86130796e+00 2.67849147e+01 -1.44118310e+01 | -3.86130796e+00 2.67849147e+01 -1.44118310e+01 8 -1.96208104e+00 -2.18092865e+00 -2.08515753e+00 | -1.96208104e+00 -2.18092865e+00 -2.08515753e+00 9 -2.94972229e+00 -2.22042312e+00 6.28930937e+00 | -2.94972229e+00 -2.22042312e+00 6.28930937e+00 10 -1.36563610e+00 1.80570070e-01 -1.14058964e+01 | -1.36563610e+00 1.80570070e-01 -1.14058964e+01 11 6.68732703e+00 3.76965910e+00 2.14389226e+00 | 6.68732703e+00 3.76965910e+00 2.14389226e+00 12 -4.95695214e+00 -8.89371072e+00 -1.33527875e+00 | -4.95695214e+00 -8.89371072e+00 -1.33527875e+00 13 8.37991226e+00 5.23807040e+00 9.35712377e-01 | 8.37991226e+00 5.23807040e+00 9.35712377e-01 14 1.40643947e+00 -1.18289155e+00 -1.56444331e+00 | 1.40643947e+00 -1.18289155e+00 -1.56444331e+00 15 3.44524188e+00 -1.17788390e+00 2.47216360e+00 | 3.44524188e+00 -1.17788390e+00 2.47216360e+00 16 -7.61439096e+00 -1.21726422e+00 -5.84092271e+00 | -7.61439096e+00 -1.21726422e+00 -5.84092271e+00 17 4.11222425e+00 4.20117889e-02 5.56296339e+00 | 4.11222425e+00 4.20117889e-02 5.56296339e+00 18 9.00525038e-01 1.85394011e+00 -2.10503207e+00 | 9.00525038e-01 1.85394011e+00 -2.10503207e+00 19 -1.74108295e+01 1.12231467e+01 -1.18497396e+01 | -1.74108295e+01 1.12231467e+01 -1.18497396e+01 20 -1.34521823e+00 2.90599432e+00 -1.14607723e+00 | -1.34521823e+00 2.90599432e+00 -1.14607723e+00 21 3.09669927e+00 5.08953635e-01 3.46375938e-01 | 3.09669927e+00 5.08953635e-01 3.46375938e-01 22 -1.43520419e+00 -2.14492366e+00 -4.46412235e+00 | -1.43520419e+00 -2.14492366e+00 -4.46412235e+00 23 4.68948196e-01 3.26076247e+00 4.09429843e+00 | 4.68948196e-01 3.26076247e+00 4.09429843e+00 24 -6.08513205e+00 -3.35456785e+00 -6.93459034e+00 | -6.08513205e+00 -3.35456785e+00 -6.93459034e+00 25 1.86210795e+00 6.97624244e+00 1.51935075e+00 | 1.86210795e+00 6.97624244e+00 1.51935075e+00 26 -1.74938301e+00 2.60939247e+00 -3.82910048e+00 | -1.74938301e+00 2.60939247e+00 -3.82910048e+00 27 3.08155910e+00 1.07965988e+00 1.71496980e+00 | 3.08155910e+00 1.07965988e+00 1.71496980e+00 28 1.53575947e+00 -6.95661954e+00 4.75771811e+00 | 1.53575947e+00 -6.95661954e+00 4.75771811e+00 29 2.68538446e+00 1.33045201e+00 2.15060618e+00 | 2.68538446e+00 1.33045201e+00 2.15060618e+00 30 1.21555807e+01 -8.94206407e+00 5.89139709e+00 | 1.21555807e+01 -8.94206407e+00 5.89139709e+00 31 3.83562511e+00 2.46179435e+00 4.06668439e+00 | 3.83562511e+00 2.46179435e+00 4.06668439e+00 32 -4.33019947e+00 -3.67226355e+00 -3.48596059e+00 | -4.33019947e+00 -3.67226355e+00 -3.48596059e+00 33 2.20137893e+00 1.79338191e+00 3.30362249e+00 | 2.20137893e+00 1.79338191e+00 3.30362249e+00 34 1.23779569e+01 -2.61097964e+01 1.40277220e+01 | 1.23779569e+01 -2.61097964e+01 1.40277220e+01 35 -1.36271897e+00 -3.17145820e+00 2.17276538e+00 | -1.36271897e+00 -3.17145820e+00 2.17276538e+00 36 -2.40198293e+00 -2.36407652e+00 1.65590071e-01 | -2.40198293e+00 -2.36407652e+00 1.65590071e-01 37 -1.32553000e+00 -2.04288423e+00 4.76671741e+00 | -1.32553000e+00 -2.04288423e+00 4.76671741e+00 38 -5.99234842e+00 -4.38662287e+00 -3.96313645e+00 | -5.99234842e+00 -4.38662287e+00 -3.96313645e+00 39 7.10574670e-01 -8.71739691e-01 8.20573814e+00 | 7.10574670e-01 -8.71739691e-01 8.20573814e+00 40 -3.32675029e+00 -2.90239223e+00 -2.49880373e+00 | -3.32675029e+00 -2.90239223e+00 -2.49880373e+00 41 1.00428073e+01 -7.33433467e+00 -3.07309878e+00 | 1.00428073e+01 -7.33433467e+00 -3.07309878e+00 42 -4.80394189e+00 -4.55028858e+00 -2.56563273e+00 | -4.80394189e+00 -4.55028858e+00 -2.56563273e+00 43 6.07098816e+00 5.35704859e+00 4.26461874e+00 | 6.07098816e+00 5.35704859e+00 4.26461874e+00 44 -3.21185934e+00 2.13003323e+00 -2.59503934e+00 | -3.21185934e+00 2.13003323e+00 -2.59503934e+00 45 1.13361358e+00 6.77196397e-01 1.21090684e+00 | 1.13361358e+00 6.77196397e-01 1.21090684e+00 46 -3.87975886e-01 5.36315522e+00 4.43171649e+00 | -3.87975886e-01 5.36315522e+00 4.43171649e+00 47 2.37697149e+00 -3.24446144e-01 1.77620952e+00 | 2.37697149e+00 -3.24446144e-01 1.77620952e+00 48 -6.96545187e+00 4.53945322e+00 7.57675629e+00 | -6.96545187e+00 4.53945322e+00 7.57675629e+00 49 -5.13399547e+00 -1.27508216e+00 2.03827255e+00 | -5.13399547e+00 -1.27508216e+00 2.03827255e+00 50 4.00965644e+00 4.86296459e+00 -2.44720101e-01 | 4.00965644e+00 4.86296459e+00 -2.44720101e-01 51 5.78065107e-01 -1.57111410e+00 -2.79845579e+00 | 5.78065107e-01 -1.57111410e+00 -2.79845579e+00 52 -1.41416260e+00 4.81981054e+00 -4.84885246e+00 | -1.41416260e+00 4.81981054e+00 -4.84885246e+00 53 2.23819966e+00 3.17306478e+00 5.38164355e+00 | 2.23819966e+00 3.17306478e+00 5.38164355e+00 54 -1.18777048e+00 1.68457490e-01 3.19916523e-01 | -1.18777048e+00 1.68457490e-01 3.19916523e-01 55 5.84154585e-02 7.98614652e-01 1.40300912e+00 | 5.84154585e-02 7.98614652e-01 1.40300912e+00 56 -6.81178655e-01 -3.29383508e+00 -1.57286390e+00 | -6.81178655e-01 -3.29383508e+00 -1.57286390e+00 57 9.84432079e+00 -6.58131138e+00 -1.14841834e+01 | 9.84432079e+00 -6.58131138e+00 -1.14841834e+01 58 6.97228946e-02 -8.79917414e-03 -2.23453171e+00 | 6.97228946e-02 -8.79917414e-03 -2.23453171e+00 59 2.11601054e+00 1.05645456e+00 2.90884254e+00 | 2.11601054e+00 1.05645456e+00 2.90884254e+00 60 6.60494387e-01 2.33690386e-01 3.25420580e-01 | 6.60494387e-01 2.33690386e-01 3.25420580e-01 61 1.22349105e+00 8.98323466e-01 1.44648297e+00 | 1.22349105e+00 8.98323466e-01 1.44648297e+00 62 -6.60095585e+00 1.26431966e+01 1.39943834e+01 | -6.60095585e+00 1.26431966e+01 1.39943834e+01 63 -1.33593648e-01 5.01655928e-01 6.15649490e-01 | -1.33593648e-01 5.01655928e-01 6.15649490e-01 ------------------------------------------------------------------------------------------------------------------------- PASS: Energies and forces are the same to within a relative error of 1e-08 ------------------------------------------------------------------------------------------------------------------------- ERROR: Unable to perform verification check. Message = ERROR: Cannot find a working configuration within a reasonable lattice constant range.