计算化学公社

标题: cp2k优化分子筛体系报错 [打印本页]

作者
Author: yaoyuan0711 时间: 2023-3-15 08:57
标题: cp2k优化分子筛体系报错
本帖最后由 yaoyuan0711 于 2023-3-15 10:36 编辑

各位老师好，我用cp2k和GFN1-xTB变胞优化一个1200个原子的分子筛体系，输入文件为multiwfn产生，报错提示Index to radix array not found，这个该怎么办？附件有完整输入输出文件
另外这种体系用DFT方法能算动吗？优化用半经验方法，单点用DFT方法，用于发文章，精度能否接受呢？
(, 下载次数 Times of downloads: 13)

作者
Author: sobereva 时间: 2023-3-15 11:40
报错可能是程序bug。尝试2023试试，也不行的话可以问开发者。或者不做变胞优化再试（如果本身就是合理的晶体结构），也可以改用PM6试试

单点用CP2K在DFT下做毫无问题，用Multiwfn默认设置，在比较好的双路服务器下SCF一步绝对超不过一分钟

作者
Author: biogon 时间: 2023-4-19 10:06
官方给的解决方法
by default CP2K uses the FFT roots available in the internal
FFT library. This is currently restricted to a maximum
of 1024. Your input requires a length of about 1200.
If you compiled cp2k with an external FFT library (probably FFTW3)
you can switch on extended lengths with the keyword
EXTENDED_FFT_LENGTHS in the GLOBAL section.

作者
Author: yaoyuan0711 时间: 2023-4-19 12:25

biogon 发表于 2023-4-19 10:06
官方给的解决方法
by default CP2K uses the FFT roots available in the internal
FFT library. This is ...

好的，谢谢您

作者
Author: 楚瑶 时间: 2023-9-16 22:30

sobereva 发表于 2023-3-15 11:40
报错可能是程序bug。尝试2023试试，也不行的话可以问开发者。或者不做变胞优化再试（如果本身就是合理的晶 ...

sob老师，请问这个SCF一步是指一轮迭代还是一小步的用时呢，我1300原子，GFN1-xTB，64核心，一小步SCF迭代需要45s，这个速度感觉有点慢，请问一下是正常速度吗。

作者
Author: sobereva 时间: 2023-9-18 09:15

楚瑶发表于 2023-9-16 22:30
sob老师，请问这个SCF一步是指一轮迭代还是一小步的用时呢，我1300原子，GFN1-xTB，64核心，一小步SCF迭 ...

我不知道你怎么定义的所谓的“一小步”。没有前提的情况下不存在这种说法

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