计算化学公社

标题: 求助，关于cfour输入文件读取不正常的问题 [打印本页]

作者
Author: commanderwilson 时间: 2025-8-19 19:27
标题: 求助，关于cfour输入文件读取不正常的问题
我对于一个六个原子的体系打算进行CCSD(T)水平的几何优化，输入文件：
(, 下载次数 Times of downloads: 1) (, 下载次数 Times of downloads: 4)
但是我运行并行计算的时候，它出现了严重的报错：
(, 下载次数 Times of downloads: 3)
我检查了一下，我明明计算的是六个原子的体系，但是似乎它只读取了我的第一个原子的坐标，导致它实际上在对单个原子做计算。
大佬们能否讲讲，我这个问题出在哪了？

作者
Author: zjxitcc 时间: 2025-8-19 20:09
本帖最后由 zjxitcc 于 2025-8-19 20:15 编辑

我看到你使用了MOKIT的fch2cfour小程序，我们在MOKIT手册上有写fch2cfour仅支持CFOUR >= v2.1。而您用的是V0.1。不同CFOUR大版本之间差距巨大，有的甚至无法读取轨道。自从我入行量子化学以来只见过CFOUR v1.0, v2.00beta和v2.1，这是我第一次见V0.1，不知道它年龄几何。。

另外，即使您未来使用CFOUR >= v2.1，也建议先从HF或CCSD(T)单点计算开始测试，如果一切没问题再考虑结构优化。以防CFOUR程序安装有问题，或者mpi并行有问题，又或者CFOUR有新的大改动版本release。

由于现在ORCA也有CCSD(T)解析导数，一般更推荐使用ORCA进行CCSD(T)结构优化。MOKIT也提供fch2mkl小程序用于Gaussian -> ORCA 传轨道。

作者
Author: sobereva 时间: 2025-8-20 02:44
除非你还要用CFOUR做振动分析计算，为了有解析Hessian而用CFOUR，否则非常不推荐用CFOUR干这事。更友好的ORCA和PSI4都有CCSD(T)解析梯度

作者
Author: Bacchus 时间: 2025-9-8 20:51
本帖最后由 Bacchus 于 2025-9-8 20:52 编辑

除了这两位老师说的情况外，还有一个可能，就是输入格式的问题，cfour的输入文件对于格式要求比较严格，在内坐标书写部分，要求必须是字符和字符之间的间隔为1个空格，我看你的空格很多，这就会导致程序再读取的时候读完第一个原子后，空一格读取第二个时什么都没有，所以不停的再优化第一个原子。解决的方法就是删去多余的空格。

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