狄拉克自带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

非常感谢！！

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

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

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

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

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

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

你用的是g16还是g09，有的版本的Dushin不支持g16，需要将输出文件开头的g16改成g09；另外Dushin也可能不能识别一些泛函，需要把输出文件的泛函名称改成b3lyp

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

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

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

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

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

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

万分感谢sob老师的解答！

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

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

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

没有遇到过这个问题，但是从fortran角度分析，应该是BN.log里面有一个数太大，溢出了。
仔细检查一下BN.log里面有没有哪个本该是数字的地方显示为一串星号