计算化学公社

标题: 使用fchred的过程中遇到fortran报错，如何处理？ [打印本页]

作者
Author: yaonaite 时间: 2026-5-1 14:18
标题: 使用fchred的过程中遇到fortran报错，如何处理？
本帖最后由 yaonaite 于 2026-5-1 15:33 编辑

在gaussian任务正常结束后，fchk文件大小为72.4 GB，希望使用fchred进行转换。但遇到了fortran的报错，有大佬遇到过这样的问题吗？如何处理？

原计算任务：
#p td=(nstates=100) b3lyp/6-31g(d,p) em=gd3bj g09default

报错信息为
forrtl: severe (64): input conversion error, unit -5, file Internal List-Directed Read

使用该文件的意图是：作部分重要的轨道图像、读轨道能量、作电子激发分析

作者
Author: sobereva 时间: 2026-5-1 15:32
检查文件末尾是否完整
如果文件被破坏、数据有异常（被记录为星号）等都可能导致fchred读取失败
处理大体系之前先用同样的关键词搞小体系进行测试

作者
Author: yaonaite 时间: 2026-5-1 15:50

sobereva 发表于 2026-5-1 15:32
检查文件末尾是否完整
如果文件被破坏、数据有异常（被记录为星号）等都可能导致fchred读取失败
处理大体 ...

感谢社长的回答！
较小的体系（没有很小，fchk文件也在数十GB量级）进行的同关键词计算没有出现该问题，能够正常转换。
该文件似乎超过了geany能打开的极限，看不了数据有没有问题，是否有别的软件可看呢？

另外，如果使用同一博文中介绍的使用multiwfn转换的方法，不会报错，转换出的fch格式文件会比fchred报错后转换出的文件更小一些（359MB vs 717MB），请问老师这个文件是不是也足够用呢？

作者
Author: sobereva 时间: 2026-5-1 17:15

yaonaite 发表于 2026-5-1 15:50
感谢社长的回答！
较小的体系（没有很小，fchk文件也在数十GB量级）进行的同关键词计算没有出现该问题， ...

Multiwfn转出来的fch文件不会包含组态系数信息

用tail命令可以输出最后部分
ultraedit可以快速打开巨大文件

作者
Author: yaonaite 时间: 2026-5-1 23:20

sobereva 发表于 2026-5-1 17:15
Multiwfn转出来的fch文件不会包含组态系数信息

用tail命令可以输出最后部分

原fchk文件的末尾是：
  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00
Gaussian Version                         C N=          2
ES64L-G16RevC.01
看上去和能正常转换的文件末尾是一样的。

我对比了转换出的新的fchk文件，发现报错的那些文件最终写到Excited state spins字段的末尾就会截止，就是说会缺少Mulliken Charges及之后的字段。随后我使用tail命令输出了最后5000行（此时已经能完整看到原fchk文件的Mulliken Charges及之后的字段），并未发现有数据变成*的问题。

请问老师我有办法进一步定位问题吗？是不是把原fchk文件的这个字段直接粘贴进新的（转换过程中报错的）fchk文件末尾也能得到预期的新fchk文件？

作者
Author: sobereva 时间: 2026-5-2 07:25

yaonaite 发表于 2026-5-1 23:20
原fchk文件的末尾是：
0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000 ...

可以试试
程序原理也很简单，就是去掉两个字段，自己结合当前文件的情况写个类似的程序也可以

作者
Author: lyj714 时间: 2026-5-2 18:04
给个chk文件呗，我试一下。之前写过一个程序直接转成fchk只保留Multiwfn所需的信息。

作者
Author: yaonaite 时间: 2026-5-5 22:25

lyj714 发表于 2026-5-2 18:04
给个chk文件呗，我试一下。之前写过一个程序直接转成fchk只保留Multiwfn所需的信息。

哎呀，这个已经删掉了 : (

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