分享一篇审稿的意见，不知该说什么好三个审稿人，两个给中，我给拒了，最后是接受了，还是个很好的期刊AEM

xexlalalan 发表于 2021-7-10 18:56
我一直也感觉VASP的算法计算的不太靠谱，也没有能检验它算的波函数的稳定性的方法。

周期性体系的波函数稳定性检测至今都还没有成熟可用的代码

xexlalalan 发表于 2021-7-10 18:56
我一直也感觉VASP的算法计算的不太靠谱，也没有能检验它算的波函数的稳定性的方法。

磁性这个事情，初猜比较重要，各种优化也就那样了

乘风万里 发表于 2021-7-7 16:30
VASP自带的磁性优化算法不一定能优化出基态，单独体系一般都是手动设MAGMOM去算铁磁反铁磁的。而一般的高 ...

我一直也感觉VASP的算法计算的不太靠谱，也没有能检验它算的波函数的稳定性的方法。

xexlalalan 发表于 2021-7-1 12:16
QE我没有用过，不太清楚，我只用过VASP和CP2K。个人认为VASP最大的优势是磁性体系的计算。VASP不需要指定 ...

VASP自带的磁性优化算法不一定能优化出基态，单独体系一般都是手动设MAGMOM去算铁磁反铁磁的。而一般的高通量算晶体流程（比如Materials project），算的磁结构都不准的，优化出的绝大部分都是铁磁相。

xexlalalan 发表于 2021-7-1 05:16
QE我没有用过，不太清楚，我只用过VASP和CP2K。个人认为VASP最大的优势是磁性体系的计算。VASP不需要指定 ...

这个功能，用任何其他的软件搭配一个二三十行的脚本都可以做，无非就是不断改变自旋多重度直至能量最低

北大-陶豫 发表于 2021-6-12 16:11
以前做量子化学计算的，不太懂第一性原理，想请大佬简单介绍一下，VASP、CP2K、Quantum Expresso 三种程序 ...

QE我没有用过，不太清楚，我只用过VASP和CP2K。个人认为VASP最大的优势是磁性体系的计算。VASP不需要指定体系的自旋多重度，可以自动优化磁矩和自旋多重度（不知道他用的是什么算法）。因此VASP在计算物理，还有具有磁性性质的固体计算中具有很大的优势。
不太清楚QE是不是也能这么做

喵星大佬 发表于 2021-6-14 19:42
能利用局域性加速的前提是不能是导体，导体能带结构是连续的，OT方法不考虑smearing会有各种包括不限于不 ...

突然想起来sob老师做过dmol3、orca和某不能说软件的benchmark：
http://sobereva.com/508
根据这个结果，三四十个原子的体系，dmol3是和orca类似或略慢于orca的，如果两者有crossover point，可能在50个原子以上

卡开发发 发表于 2021-6-16 00:09
过渡金属体系还是最好16。转换不麻烦，推荐ase，直接xsd或car生成gjf，然后再改改就行了。

我都是导出mol然后用gv改的

喵星大佬 发表于 2021-6-15 20:19
数值没具体对比过，只能说只展开到四级距出来的结构基本合理。
不过在默认的8极下找过一些模型反应的过 ...

过渡金属体系还是最好16。转换不麻烦，推荐ase，直接xsd或car生成gjf，然后再改改就行了。

wzkchem5 发表于 2021-6-15 19:49
按你经验，16极矩和4极矩的结果（绝对能量、结构之类的）大概能差多少？

数值没具体对比过，只能说只展开到四级距出来的结构基本合理。
不过在默认的8极下找过一些模型反应的过渡态，极小点和过渡态结构跟XX在差不多级别(同泛函，DNP4.4/6-31G**)下优化出来的几乎是重合的。

所以我说即使是预优化的情况下其实也能得到相当不错的初猜结构，后面优化用XX收敛速度非常快，一般不超过20步(切成了杂化泛函)，唯一的麻烦事就是这两个中间导结构。。。。

wzkchem5 发表于 2021-6-15 17:00
RI对于大体系确实不如DMol那种multipolar expansion，因为RI涉及Coulomb metric的Cholesky分解，而Coulom ...

这只是DMol3最高展开到16，虽然如此，如果非要判断误差情况你可以很简单比较8极和16极的收敛性测试，我想这不难实现。这样的数值计算中，做好收敛性测试工作尤为重要，并不应该对着默认值随手一选。当然，可能很多程序给出的默认值都是经过一些常见体系测试的，但是也许也有一些程序的参数可能就是随手提供的。基组不同固然不好比较积分，但是不表示这些基组不能比较误差，能量、本征值等结果仍然可以比较，这也是基组赝势构造应该要保证的问题。

喵星大佬 发表于 2021-6-15 10:15
最高也就到16极距，一般以预优化为目的用的时候只设到四级距

按你经验，16极矩和4极矩的结果（绝对能量、结构之类的）大概能差多少？

wzkchem5 发表于 2021-6-15 17:00
RI对于大体系确实不如DMol那种multipolar expansion，因为RI涉及Coulomb metric的Cholesky分解，而Coulom ...

最高也就到16极距，一般以预优化为目的用的时候只设到四级距

卡开发发 发表于 2021-6-15 03:01
是，cp2k收敛方面问题现在我也不知道怎么样，以前算个有磁性的体系大概率是难收敛的，网格和混合参数全都 ...

RI对于大体系确实不如DMol那种multipolar expansion，因为RI涉及Coulomb metric的Cholesky分解，而Coulomb metric是一个O(N)*O(N)的dense matrix，所以这一步是O(N^3)的，且不能通过prescreening降低scaling。相比之下multipolar expansion做prescreening可以做到O(N^2)。至于是否用数值基，应该也会有一些影响，不过只影响prefactor。
再者我觉得DMol可能还有一个优势，就是数值基没办法算解析积分，导致别人不能拿解析积分的结果来比，所以DMol的库伦项可以做得很粗糙，比如多极展开好像只展开到8极矩。像BDF这种在高斯基上面做多极展开的程序，就要展开到256极矩，把总能量误差做到比如说0.1mHartree的量级，不然别人分分钟做一个解析积分的计算，告诉你说你的误差竟然有1mHartree，我怎么敢用？相比之下，DMol哪怕跟解析积分的结果差100个mHartree，也没人知道，谁也不知道是基组误差还是多极展开的误差。