ok,thank you a lot! I will check it again!<div class="gmail_quote"><div dir="auto" class="gmail_attr">在2022年12月10日星期六 UTC+1 17:51:03<Matthias Krack> 写道： </div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"> <div lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"> <div> 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 <a href="https://github.com/cp2k/cp2k/blob/master/tests/Fist/regtest-7-1/UO2-4x4x4-shell-forces.inp" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=zh-CN&q=https://github.com/cp2k/cp2k/blob/master/tests/Fist/regtest-7-1/UO2-4x4x4-shell-forces.inp&source=gmail&ust=1670919140211000&usg=AOvVaw34aq-h2p99__Kcyq_deCfl"> UO2 test inputs</a> 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</div></div><div lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"><div> Matthias <div style="border:none;border-top:solid #b5c4df 1.0pt;padding:3.0pt 0cm 0cm 0cm"> From: <a href data-email-masked rel="nofollow">cp...@googlegroups.com</a> <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>> on behalf of elephant gogogo! <<a href data-email-masked rel="nofollow">qumingzi...@gmail.com</a>> Date: Friday, 9 December 2022 at 22:02 To: cp2k <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>> Subject: Re: [CP2K:18169] williams potential </div> 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? <div> </div> <div> Best Regards </div> <div> </div> <div> <div> 在2022年12月9日星期五 UTC+1 17:40:46<Matthias Krack> 写道： </div> <blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm"> <div> <div> 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 </div> </div> <div> <div> <div style="border:none;border-top:solid #b5c4df 1.0pt;padding:3.0pt 0cm 0cm 0cm"> From: <a href data-email-masked rel="nofollow">cp...@googlegroups.com</a> <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>> on behalf of elephant gogogo! <<a href data-email-masked rel="nofollow">qumingzi...@gmail.com</a>> Date: Friday, 9 December 2022 at 13:07 To: cp2k <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>> Subject: [CP2K:18166] williams potential </div> Dear CP2K users, <div> </div> <div> 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. </div> <div> </div> <div> Thank you for your time! </div> </div> </div> <div> <div> -- 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 <a href data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a>. To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/747656f5-ceef-45f3-8762-8bf8b40265den%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=zh-CN&q=https://groups.google.com/d/msgid/cp2k/747656f5-ceef-45f3-8762-8bf8b40265den%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1670919140211000&usg=AOvVaw2di2CKq3lxBtrlIHcDYjmu"> https://groups.google.com/d/msgid/cp2k/747656f5-ceef-45f3-8762-8bf8b40265den%40googlegroups.com</a>. </div> </div> </blockquote> </div> </div></div><div lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"><div>-- 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 <a href data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a>. </div></div><div lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"><div> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/46dbe775-65f8-4ed7-aab2-a45cfa7416afn%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=zh-CN&q=https://groups.google.com/d/msgid/cp2k/46dbe775-65f8-4ed7-aab2-a45cfa7416afn%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1670919140212000&usg=AOvVaw3U32y5gL-R-qjAfsFyH35m"> https://groups.google.com/d/msgid/cp2k/46dbe775-65f8-4ed7-aab2-a45cfa7416afn%40googlegroups.com</a>. </div> </div> </blockquote></div> -- 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 <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>. To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/6d2efa4c-2c77-4c5d-a6da-a877adf5c302n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/6d2efa4c-2c77-4c5d-a6da-a877adf5c302n%40googlegroups.com</a>.