ghifi37 发表于 2021-3-24 09:51
对于复杂的体系，量化软件里也应该要有寻找能量最低自旋多重度的算法

比较电子能量就能解决的问题，问题在于如何搞准这个东西，涉及到不黑箱的多参考实际上问题很多

对于复杂的体系，量化软件里也应该要有寻找能量最低自旋多重度的算法

jiangning198511 发表于 2021-3-23 09:59
以前想研究一下VASP是如何实现的，看了看代码没搞明白

那好啊，我抽空研究下，到时候欢迎一块讨论。

卡开发发 发表于 2021-3-23 09:16
也疑惑这个问题，不过找了诸多资料我不是很明确具体的算法，有时间拆解某个程序看看再一块讨论。不过从程 ...

以前想研究一下VASP是如何实现的，看了看代码没搞明白

jiangning198511 发表于 2021-3-23 09:07
卡卡大佬，在VASP中自旋极化计算中程序可以自动优化到一个自旋状态，当然也可以指定自旋，但可能不是基态 ...

也疑惑这个问题，不过找了诸多资料我不是很明确具体的算法，有时间拆解某个程序看看再一块讨论。不过从程序设计来说，这些程序可能是把密度和自旋密度分开放而不是alpha和beta电子密度，应该是有些说法。

卡开发发 发表于 2021-3-23 06:32
我指的是都尽可能手工制定磁矩。如果指的是磁矩初始猜测的值，不同程序可能会不同，但这种差异需要看手册 ...

卡卡大佬，在VASP中自旋极化计算中程序可以自动优化到一个自旋状态，当然也可以指定自旋，但可能不是基态。这种自旋优化的算法是如何实现的我一直没有理解，能否详细解答一下

Penson 发表于 2021-3-22 19:21
嗯嗯，cp2k和vasp应该都有各自的规则吧？

我指的是都尽可能手工制定磁矩。如果指的是磁矩初始猜测的值，不同程序可能会不同，但这种差异需要看手册，vasp全按照1猜，也有程序按0猜的。

喵星大佬 发表于 2021-3-22 19:02
GTO基处理PBC稳定性差是特指GTO还是包括STO，NAO在内的所有原子中心基

具体取决于积分的形式。在real space下的基函数进行积分需要对不同的周期R求和，但这个求和在数值上的处理是有限的，这种截断可能就会导致数值上一定程度的不稳定。数值基组本身有截断可以严格避免这个问题，当然也可以通过类似于Ewald的形式去做。

Penson 发表于 2021-3-22 19:21
嗯嗯，cp2k和vasp应该都有各自的规则吧？

对。同类程序都是相同规则。

卡开发发 发表于 2021-3-21 22:51
几乎没有同类软件能帮忙找合适的磁序。

嗯嗯，cp2k和vasp应该都有各自的规则吧？

卡开发发 发表于 2021-3-22 09:21
1、固体计算当中，涉及到有磁性的体系，对原子磁矩进行初始猜测是必须进行的，这应该作为基本常识。不同的 ...

GTO基处理PBC稳定性差是特指GTO还是包括STO，NAO在内的所有原子中心基

1、固体计算当中，涉及到有磁性的体系，对原子磁矩进行初始猜测是必须进行的，这应该作为基本常识。不同的程序对开启自旋极化的情况有一套默认的规则来确定原子初始磁矩，但这种默认的方式并不一定合理，有些体系可能多种磁序都是能算出来，但并尽然是我们需要的结果。当然，很复杂的体系甚至可能磁序的排列并不是随便猜几次就能得到，这种情况建议慎重一些，有些体系可能磁性对结构会比较敏感。

2、GTO的基组本身处理PBC稳定性上就要稍差一些（有专门讨论基组方面的文献），如果基组本身还要弥散一些稳定性就要更差一些，这个和长程上的截断的处理方式有一定的关系，所以同样等级的基组GTO会出现难收敛的情况也并不意外。

3、赝势和基组的问题不要小看，如果审稿人真的针对数据询问轻则补算单点重则全部重来，delta test的数据当然只是参考但总比没参考要来得强。

4、窄带的情况可能对k依赖会很大，k不充分的情况还做固定占据的处理可能会得到完全不正确的结果。

喵星大佬 发表于 2021-3-20 02:25
cp2k算磁性体系可能会死人吧

不复杂，&KIND里设个MAGNETIZATION指定初猜中此类原子alpha减beta电子数（或者在&BS里指定初猜的原子电子组态），MULTIPLICITY设成实际整体自旋多重度即可，并无困难，速度也很快。若真碰上很复杂的体系，可能需要写脚本将各相关原子的MAGNETIZATION设置做个排列组合。

Penson 发表于 2021-3-21 21:32
cp2k磁性得自己设定，软件没法帮忙帮忙找合适的吧？

几乎没有同类软件能帮忙找合适的磁序。