使用fchred的过程中遇到fortran报错，如何处理？

yaonaite

在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


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

sobereva

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

yaonaite

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

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

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

sobereva

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

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

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

yaonaite

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文件？