Hi
The force field is not properly defined in your input, e.g. a SHELL_CHARGE is set but no CORE_CHARGE. It is difficult to check the other parameters, because you did not specify the force field parameters for CeO2 which you want to use. The &CHARGE sections for rigid ion potentials should not appear when you specify &SHELL sections for a core-shell model. There is also no need to specify initially in addition to the &COORD section a &SHELL_COORD and a &CORE_COORD section. The &KIND sections can be dropped, because these are not considered by the Fist module. I am also not sure, if the lattice vectors and the atomic coordinates are matching and thus will represent the fluorite structure of CeO2 with PBC properly.
Best
Matthias
--
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to
cp2k+uns...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/cp2k/747656f5-ceef-45f3-8762-8bf8b40265den%40googlegroups.com.
CORE_CHARGE and SHELL_CHARGE must sum up to the desired formal charge which is usually +4 or +3 for Ce. Just take one of the UO2 test inputs for the CP2K/Fist module and fill in the units and values from the force field parameter table in the literature. Both CeO2 and UO2 have the same fluorite structure and you can substitute U by Ce. I would also cross-check the force field parameters with other papers, because such tables contain often typos.
HTH
To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/46dbe775-65f8-4ed7-aab2-a45cfa7416afn%40googlegroups.com.