[CP2K-user] [CP2K:18170] williams potential

[Krack Matthias]

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<https://github.com/cp2k/cp2k/blob/master/tests/Fist/regtest-7-1/UO2-4x4x4-shell-forces.inp> 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

Matthias

From: cp2k at googlegroups.com <cp2k at googlegroups.com> on behalf of elephant gogogo! <qumingzizhemenan at gmail.com>
Date: Friday, 9 December 2022 at 22:02
To: cp2k <cp2k at googlegroups.com>
Subject: Re: [CP2K:18169] williams potential
Thank you matthias, I just get the potential from the literature, the literature only use the shell charge, so I don't know how to fix this mistake. But for the charge section, the potential has also columbic potential, so I should define the charge?

Best Regards

在2022年12月9日星期五 UTC+1 17:40:46<Matthias Krack> 写道：
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

From: cp... at googlegroups.com <cp... at googlegroups.com> on behalf of elephant gogogo! <qumingzi... at gmail.com>
Date: Friday, 9 December 2022 at 13:07
To: cp2k <cp... at googlegroups.com>
Subject: [CP2K:18166] williams potential
Dear CP2K users,

I have found 5 different potentials from the literatures, but when I use it in the cp2k, it will always have the mistakes, e.g.,geometry wrong. I assume, maybe the potential is not so accurate or not suitble for cp2k. So how could I deal with this problem? Please find one of my inputs in the attachment.

Thank you for your time!
--
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... at googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/747656f5-ceef-45f3-8762-8bf8b40265den%40googlegroups.com<https://groups.google.com/d/msgid/cp2k/747656f5-ceef-45f3-8762-8bf8b40265den%40googlegroups.com?utm_medium=email&utm_source=footer>.
--
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+unsubscribe at googlegroups.com<mailto:cp2k+unsubscribe at googlegroups.com>.
To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/46dbe775-65f8-4ed7-aab2-a45cfa7416afn%40googlegroups.com<https://groups.google.com/d/msgid/cp2k/46dbe775-65f8-4ed7-aab2-a45cfa7416afn%40googlegroups.com?utm_medium=email&utm_source=footer>.

-- 
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+unsubscribe at googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/ZRAP278MB0827957DD26B0581A9B5D58AF41F9%40ZRAP278MB0827.CHEP278.PROD.OUTLOOK.COM.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20221210/203a6205/attachment-0001.htm>