计算化学公社

标题: 请问一下大家有siesta LDA和PBE赝势库的网址或者压缩包吗 [打印本页]

作者
Author: dior_w 时间: 2024-9-12 16:04
标题: 请问一下大家有siesta LDA和PBE赝势库的网址或者压缩包吗
我在官网上只找到了GGA的赝势库，但是计算感觉总是不尽如人意，想问问大家有没有LDA和PBE的赝势库？谢谢！

作者
Author: 卡开发发 时间: 2024-9-13 00:15
GGA的那个就是PBE的，SIESTA应该没有专门针对PW91做过赝势，LDA的应该只要把链接地址从GGA改成LDA就行了，但也不是很推荐使用。

SIESTA的稳定性其实还可以，但赝势和基组质量实在太糟糕了，也许你可以使用ONCV-Dojo的赝势或是专门优化的基组，但精度还是不尽如人意。所以尽管SIESTA我用了很多年但很少和周围的人推。

作者
Author: dior_w 时间: 2024-9-13 10:21

卡开发发发表于 2024-9-13 00:15
GGA的那个就是PBE的，SIESTA应该没有专门针对PW91做过赝势，LDA的应该只要把链接地址从GGA改成LDA就行了， ...

好的特别感谢前辈，这个ONCV-Dojo的赝势请问您有网站吗？还有最近在计算siesta的时候遇到了Block Chemical_species_label does not exist.这个问题，请问您知道怎么解决吗？谢谢！

作者
Author: 卡开发发 时间: 2024-9-13 11:11

dior_w 发表于 2024-9-13 10:21
好的特别感谢前辈，这个ONCV-Dojo的赝势请问您有网站吗？还有最近在计算siesta的时候遇到了Block Chemica ...

ONCV-Dojo赝势你可以到这里下载，不过我一直没下载下来过，而且实际上基组也不一定确保能适配Dojo的赝势。需要注意的是Dojo赝势使用需要用特定版本或是新版，比较折腾。
Block Chemical_species_label does not exist.这个问题我不确定，如果实在不行你可以打包你的输入文件，暂时看不出问题来。

作者
Author: 薛定谔的老猫 时间: 2024-10-12 15:44

卡开发发发表于 2024-9-13 00:15
GGA的那个就是PBE的，SIESTA应该没有专门针对PW91做过赝势，LDA的应该只要把链接地址从GGA改成LDA就行了， ...

我使用ONCV-Dojo的赝势库做了计算，和VASP进行了比较，还是存在差异，且定性存在错误，感觉赝势精度很难保证，准备用Openmx了，不知道会不会好一些？

作者
Author: 卡开发发 时间: 2024-10-12 16:25

薛定谔的老猫发表于 2024-10-12 15:44
我使用ONCV-Dojo的赝势库做了计算，和VASP进行了比较，还是存在差异，且定性存在错误，感觉赝势精度很难 ...

从delta测试的表现看openmx会比siesta现在的情况好不少。

作者
Author: 薛定谔的老猫 时间: 2024-10-12 16:41

卡开发发发表于 2024-10-12 16:25
从delta测试的表现看openmx会比siesta现在的情况好不少。

好的呢，刚刚编译完openmx，准备测试看下，Siesta赝势太坑了。

作者
Author: 卡开发发 时间: 2024-10-13 08:51

薛定谔的老猫发表于 2024-10-12 16:41
好的呢，刚刚编译完openmx，准备测试看下，Siesta赝势太坑了。

OpenMX的赝势基组或许你可以参考官网2019版本推荐赝势基组的来。我个人测试下来除了设置磁矩（等原子信息）比SIESTA麻烦一些，另外迭代的算法稳定性不如SIESTA，在用法上和SIESTA差的不太多。关于关于delta测试的方式你可以参考这里，虽然说delta测试表现好的程序或是赝势不一定实际用起来没问题，但数据难看的实际用起来多半会有坑。

作者
Author: 薛定谔的老猫 时间: 2024-10-17 15:19

卡开发发发表于 2024-10-13 08:51
OpenMX的赝势基组或许你可以参考官网2019版本推荐赝势基组的来。我个人测试下来除了设置磁矩（等原子信息 ...

这几天也在测试OpenMX，不知道您有没有测试过里面的OrderN计算，在我的机子上面，用它提供的例子，原子数小于1000时，跑起来还蛮快，但是原子数较大时，例如它的work/large3_example/Si12880*的例子，用48核跑，SCF第一步都算不了。不知道是不是内存不够了，也没其他报错，最后直接退出了。

作者
Author: 卡开发发 时间: 2024-10-17 17:50

薛定谔的老猫发表于 2024-10-17 15:19
这几天也在测试OpenMX，不知道您有没有测试过里面的OrderN计算，在我的机子上面，用它提供的例子，原子数 ...

我没有测例子。我看了下那个体系是1280原子，我个人感觉48核心够呛，但如果使用O(N)方法中Krylov子空间的方法1000多原子应该问题不大。

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