赝氢体系电子总数偶数，但提示奇数，任务中断

zhoumeimei

赝氢体系 1.5电子的赝氢21个（电子总数31.5），0.5电子赝氢35个（电子总数17.5），氢电子1个。所有电子总数31.5+17.5+1=50
但是orca stddft直接报错，提示电子数是奇数中断
请问，是必须每种赝氢的电子总数都要是偶数，才能判断为偶数正常运行吗？

wzkchem5

总电荷是按Z=1算的，而不是按你设的Z算的。这是因为orca代码要求总电荷必须是整数，如果按你设的Z算，万一你的所有Z加起来不是整数，总电荷就没法设了。
不过手册确实好像没把这一点写清楚，我更新一下orca手册

zhoumeimei

wzkchem5 发表于 2021-10-18 15:27
总电荷是按Z=1算的，而不是按你设的Z算的。这是因为orca代码要求总电荷必须是整数，如果按你设的Z算，万一 ...

那意思是orca不能算赝氢体系吗？

但是我把每种赝氢的个数都设置成偶数以后（即 H1.5*16个=24e  H0.5*36=18e），orca是可以正常跑了，只是现在还没有收敛。

wzkchem5

zhoumeimei 发表于 2021-10-20 01:50
那意思是orca不能算赝氢体系吗？

但是我把每种赝氢的个数都设置成偶数以后（即 H1.5*16个=24e  H0.5*3 ...

不是这样的，我的意思是，一个设成Z=1.5的赝氢，在算总电荷时，应当看成Z=1。不是说程序把它当成Z=1来用。
所以比如说，两个Z=1.5的赝氢组成的氢分子，总电荷应该写为0而不是1。但是计算结果确实对应于两个Z=1.5的赝氢组成的氢分子。

mjj

wzkchem5 发表于 2021-10-18 15:27
总电荷是按Z=1算的，而不是按你设的Z算的。这是因为orca代码要求总电荷必须是整数，如果按你设的Z算，万一 ...

老师，这里“总电荷是按Z=1算的”，总电荷是1的话不就是体系带电的意思吗？跟楼主说的赝氢电荷设置是一回事吗？

wzkchem5

mjj 发表于 2021-10-22 02:34
老师，这里“总电荷是按Z=1算的”，总电荷是1的话不就是体系带电的意思吗？跟楼主说的赝氢电荷设置是一回 ...

不是总电荷是1，而是赝氢的电荷看作1。楼主不是用Z=xxx设置赝氢电荷的嘛，我指的是这个Z

meltconstraint

wzkchem5 发表于 2021-10-22 14:32
不是总电荷是1，而是赝氢的电荷看作1。楼主不是用Z=xxx设置赝氢电荷的嘛，我指的是这个Z

老师，那这样该如何让程序运行呢，一直显示总电荷是奇数，程序运行不了呀

wzkchem5

meltconstraint 发表于 2024-3-27 13:52
老师，那这样该如何让程序运行呢，一直显示总电荷是奇数，程序运行不了呀

如果把所有的赝氢都改成普通氢，能正常运行吗

meltconstraint

wzkchem5 发表于 2024-3-27 22:17
如果把所有的赝氢都改成普通氢，能正常运行吗

不行，改成普通的氢，一是想要的钝化效果没有了，二是总电子数是奇数，还是带电体系了，不管是设置赝氢电荷量，还是用点电荷Q，还是用embeddeding potential，都是一样的。这个时候把自旋多重度硬改成2跑一下，结构就直接散架了。

wzkchem5

meltconstraint 发表于 2024-3-27 19:14
不行，改成普通的氢，一是想要的钝化效果没有了，二是总电子数是奇数，还是带电体系了，不管是设置赝氢电 ...

所以说明整个体系的带电量算错了。
先把赝氢换成普通氢，重新仔细计算这个团簇应该带多少电荷，然后在不改变电荷数的情况下把普通氢换回赝氢（而不管赝氢的Z是多少），再算