计算化学公社

标题: 求助dushin运行报错 forrtl: severe (64): input conversion error, unit 15 [打印本页]

作者
Author: 狄拉克自带buff 时间: 2020-10-13 12:05
标题: 求助dushin运行报错 forrtl: severe (64): input conversion error, unit 15
老师，您好，我根据帖子使用Dushin分解重组能和计算Huang-Rhys因子http://bbs.keinsci.com/forum.php ... 7&fromuid=12600，在计算时出现了报错信息：forrtl: severe (64): input conversion error, unit 15, file D:\study\Dushin_win\BN.log
请问这是什么原因？谢谢老师的解答。

作者
Author: wzkchem5 时间: 2020-10-13 17:22
没有遇到过这个问题，但是从fortran角度分析，应该是BN.log里面有一个数太大，溢出了。
仔细检查一下BN.log里面有没有哪个本该是数字的地方显示为一串星号

作者
Author: sobereva 时间: 2020-10-13 17:48
原因就是Fortran从输入文件里读取数据失败。
既可能是2L所说的情况，也可能是当前情况输出文件内容异常，或者与dushin存在不兼容性。
确保输入文件合理性的前提下，自行调试代码，或者发邮件问作者

作者
Author: 狄拉克自带buff 时间: 2020-10-13 18:02

wzkchem5 发表于 2020-10-13 17:22
没有遇到过这个问题，但是从fortran角度分析，应该是BN.log里面有一个数太大，溢出了。
仔细检查一下BN.lo ...

好的，万分感谢您的回复！

作者
Author: 狄拉克自带buff 时间: 2020-10-13 18:02

sobereva 发表于 2020-10-13 17:48
原因就是Fortran从输入文件里读取数据失败。
既可能是2L所说的情况，也可能是当前情况输出文件内容异常， ...

万分感谢sob老师的解答！

作者
Author: 狄拉克自带buff 时间: 2020-10-13 22:26

sobereva 发表于 2020-10-13 17:48
原因就是Fortran从输入文件里读取数据失败。
既可能是2L所说的情况，也可能是当前情况输出文件内容异常， ...

老师，您好。我在对比了正常结束的文件和出现报错的文件详细的计算过程后，我认为主要原因是输出的.log文件中缺失了input orientation信息。请问老师，这是不是主要原因？如果是的话，怎样找回缺失的input orientation信息？
谢谢老师的耐心解答。

作者
Author: Daniel_Arndt 时间: 2020-10-14 07:55

狄拉克自带buff 发表于 2020-10-13 22:26
老师，您好。我在对比了正常结束的文件和出现报错的文件详细的计算过程后，我认为主要原因是输出的.log文 ...

如果你觉得原因是没有Input Orientation的话，首先看一下你的体系的原子数。如果体系的原子数超过50的话，你可以尝试在高斯的输入文件中加“geom=printinputorient”关键词。

高斯的默认设置应该是体系原子数超过50时不在输出文件中写Input Orientation。如果你的体系的原子数不超过50却又没有Input Orientation的话，你估计只能联系客服了。

作者
Author: 破晓时的曙光 时间: 2020-10-14 10:04
你用的是g16还是g09，有的版本的Dushin不支持g16，需要将输出文件开头的g16改成g09；另外Dushin也可能不能识别一些泛函，需要把输出文件的泛函名称改成b3lyp

作者
Author: 狄拉克自带buff 时间: 2020-10-14 10:54

Daniel_Arndt 发表于 2020-10-14 07:55
如果你觉得原因是没有Input Orientation的话，首先看一下你的体系的原子数。如果体系的原子数超过50的话 ...

好的，谢谢您！我的体系分子确实都超过了50个，用的是G09

作者
Author: 破晓时的曙光 时间: 2020-10-14 23:40

狄拉克自带buff 发表于 2020-10-14 10:54
好的，谢谢您！我的体系分子确实都超过了50个，用的是G09

另外，不知道你优化结构算频率的时候有没有加nosymm，有时候优化过程会改变分子取向，导致最后有问题

作者
Author: 狄拉克自带buff 时间: 2020-10-15 10:59

破晓时的曙光发表于 2020-10-14 23:40
另外，不知道你优化结构算频率的时候有没有加nosymm，有时候优化过程会改变分子取向，导致最后有问题

十分感谢您的回复！
我是之前就优化好的结构，只算了他的频率，所以没有加nosymm。
经过我的核实，原因确实是没有Input Orientation信息，在添加gemo=printinputorient关键词之后已经解决了。可以用dushin正常计算了！
再次感谢您的回复

作者
Author: tiantian_dong 时间: 2024-6-20 15:30
非常感谢！！

作者
Author: 陈AG 时间: 2025-4-28 09:27
本帖最后由陈AG 于 2025-4-28 09:28 编辑

狄拉克自带buff 发表于 2020-10-15 10:59
十分感谢您的回复！
我是之前就优化好的结构，只算了他的频率，所以没有加nosymm。
经过我的核实，原因 ...

老师，请问您说的Input Orientation没有，是freq计算输出文件搜索就搜不到吗，因为我现在也遇到dushin.out文件不对的问题（http://bbs.keinsci.com/thread-53162-1-1），我不知道是不是您说的这个原因（现在在加gemo=printinputorient，再算一遍），还是G09和G16输出的Input Orientation不同，sob老师示例示例的输出文件是
   Input orientation:
---------------------------------------------------------------------
Center    Atomic    Atomic          Coordinates (Angstroms)
Number    Number    Type          X          Y          Z
---------------------------------------------------------------------
   1       7          0       0.000000 0.000000 1.416074
   2       6          0       0.000000 1.143023 0.725415
   3       6          0       0.000000 1.200182 -0.672711
   4       6          0       0.000000 0.000000 -1.386843
   5       6          0       0.000000 -1.200182 -0.672711
   6       6          0       0.000000 -1.143023 0.725415
   7       1          0       0.000000 2.068993 1.313486
   8       1          0       0.000000 2.165189 -1.185337
   9       1          0       0.000000 0.000000 -2.480201
   10       1          0       0.000000 -2.165189 -1.185337
   11       1          0       0.000000 -2.068993 1.313486
---------------------------------------------------------------------
我的输出文件是
Electric dipole moment (input orientation):
(Debye = 10**-18 statcoulomb cm , SI units = C m)
               (au)          (Debye)       (10**-30 SI)
Tot       0.588754D+00    0.149646D+01    0.499166D+01
x       0.560947D+00    0.142579D+01    0.475591D+01
y       -0.176690D+00    -0.449100D+00    -0.149804D+01
z       -0.273936D-01    -0.696276D-01    -0.232253D+00

Dipole polarizability, Alpha (input orientation).
(esu units = cm**3 , SI units = C**2 m**2 J**-1)
Alpha(0;0):
            (au)          (10**-24 esu)    (10**-40 SI)
iso       0.902001D+03    0.133663D+03    0.148720D+03
aniso    0.671716D+03    0.995380D+02    0.110751D+03
xx       0.127720D+04    0.189261D+03    0.210581D+03
yx       -0.202472D+02    -0.300032D+01    -0.333831D+01
yy       0.924508D+03    0.136998D+03    0.152431D+03
zx       0.601086D+01    0.890718D+00    0.991057D+00
zy       -0.152489D+02    -0.225965D+01    -0.251420D+01
zz       0.504300D+03    0.747296D+02    0.831479D+02

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