计算化学公社

标题: 在Gaussian输出文件中记录AtFile的内容导致Multiwfn无法绘制UV谱图 [打印本页]
作者
Author: liyuanhe211    时间: 2023-4-12 20:40
标题: 在Gaussian输出文件中记录AtFile的内容导致Multiwfn无法绘制UV谱图
本帖最后由 liyuanhe211 于 2023-4-12 20:50 编辑


当在Gaussian输入文件中需要以@符号引入其他文件内容时，在Multiwfn载入相应输出文件绘制其UV谱图时，Multiwfn会将其识别为Plain text file，而不是Gaussian output file，会导致程序崩溃。
如果在@文件时加入/N标示使Gaussian不在输出文件开头打印一遍at文件的内容则运行正常。


如附件crashed.out和crashed.gjf所示，如果输入文件中使用@/public1/home/sc50780/g09/ma-TZVP.basis引入ma-TZVP基组，载入Multiwfn中，运行11/3后程序崩溃
  1. 3
  2. Recognized as a plain text file
  3. forrtl: severe (59): list-directed I/O syntax error, unit 10, file D:\Gaussian\LHY_Spectrum\Multiwfn_plot_UV_Problem\crashed.out
  4. Image              PC                Routine            Line        Source
  5. Multiwfn.exe       00007FF7E48C5CAF  Unknown               Unknown  Unknown
  6. Multiwfn.exe       00007FF7E487820C  Unknown               Unknown  Unknown
  7. Multiwfn.exe       00007FF7E4876823  Unknown               Unknown  Unknown
  8. Multiwfn.exe       00007FF7E40940C1  Unknown               Unknown  Unknown
  9. Multiwfn.exe       00007FF7E408F607  Unknown               Unknown  Unknown
  10. Multiwfn.exe       00007FF7E4282199  Unknown               Unknown  Unknown
  11. Multiwfn.exe       00007FF7E524A24E  Unknown               Unknown  Unknown
  12. Multiwfn.exe       00007FF7E55173E4  Unknown               Unknown  Unknown
  13. KERNEL32.DLL       00007FFBCE4A7614  Unknown               Unknown  Unknown
  14. ntdll.dll          00007FFBCE6626A1  Unknown               Unknown  Unknown
复制代码
若使用@/public1/home/sc50780/g09/ma-TZVP.basis/N引入文件则可正常绘制（如附件successful.out和successful.gjf所示）

Gaussian 16 A.03 on CentOS 8 Stream
Multiwfn 3.8 dev 2023-Mar-8 on Windows 10
作者
Author: snljty2    时间: 2023-4-12 21:14
本帖最后由 snljty2 于 2023-4-12 21:28 编辑

之前我问过卢老师这个问题，记得是Multiwfn会查看输入文件的前几百行来找Gaussian字样，如果找到就认为是Gaussian的输出文件，否则就当纯文本从而报错。而因为被at的部分展开后太长，导致没找到Gaussian字样从而就报错。大概是这个意思，具体细节我记不清了。刚翻了一下源代码，是前500行找到"Gaussian, Inc"或者前200行找到“Entering Gaussian System”，就当Gaussian输出文件，找不到就继续。在util.f90。
作者
Author: sobereva    时间: 2023-4-12 21:20
就是2L说的那样

Multiwfn为了避免从头到尾搜索整个文件在某些情况下造成过高耗时，只从开头搜索一定行数的内容。一般也都建议用/N，要不然自己看输出文件也不便。
作者
Author: liyuanhe211    时间: 2023-4-13 00:58
snljty2 发表于 2023-4-12 21:14
之前我问过卢老师这个问题，记得是Multiwfn会查看输入文件的前几百行来找Gaussian字样，如果找到就认为是Ga ...

我的文件（crashed.out）第一行就有"Entering Gaussian System"，为什么没有命中呢？
  1. Entering Gaussian System, Link 0=/public1/home/sc50780/g16/g16
  2. Input=/public1/home/sc50780/Gaussian/LHY_Spectrum/A/Cis_Sub[SP_TD_wB97xD_maTZVP].gjf
  3. Output=/public1/home/sc50780/Gaussian/LHY_Spectrum/A/Cis_Sub[SP_TD_wB97xD_maTZVP].out
  4. AtFile(1): /public1/home/sc50780/g09/ma-TZVP.basis
  5.             -H     0
  6.             S   3   1.00
  7.                  34.0613410              0.60251978E-02
  8.                   5.1235746              0.45021094E-01
  9.                   1.1646626              0.20189726
复制代码
作者
Author: sobereva    时间: 2023-4-13 01:10
我今天刚更新的Multiwfn已经能兼容你的文件了
作者
Author: liyuanhe211    时间: 2023-4-13 03:09
sobereva 发表于 2023-4-13 01:10
我今天刚更新的Multiwfn已经能兼容你的文件了

测试可以了，谢谢



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3