计算化学公社

标题: 赝氢体系电子总数偶数，但提示奇数，任务中断 [打印本页]

作者
Author: zhoumeimei 时间: 2021-10-18 12:15
标题: 赝氢体系电子总数偶数，但提示奇数，任务中断
赝氢体系 1.5电子的赝氢21个（电子总数31.5），0.5电子赝氢35个（电子总数17.5），氢电子1个。所有电子总数31.5+17.5+1=50
但是orca stddft直接报错，提示电子数是奇数中断
请问，是必须每种赝氢的电子总数都要是偶数，才能判断为偶数正常运行吗？

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

作者
Author: zhoumeimei 时间: 2021-10-20 08:50

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

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

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

作者
Author: wzkchem5 时间: 2021-10-20 15:11

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的赝氢组成的氢分子。

作者
Author: mjj 时间: 2021-10-22 09:34

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

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

作者
Author: wzkchem5 时间: 2021-10-22 14:32

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

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

作者
Author: meltconstraint 时间: 2024-3-27 20:52

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

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

作者
Author: wzkchem5 时间: 2024-3-27 22:17

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

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

作者
Author: meltconstraint 时间: 2024-3-28 02:14

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

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

作者
Author: wzkchem5 时间: 2024-3-28 16:43

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)