计算化学公社

标题: dushin分解重组能forrtl：severe（24）求助 [打印本页]

作者
Author: 起名真是难 时间: 2020-5-27 17:38
标题: dushin分解重组能forrtl：severe（24）求助
各位老师好，最近我想用dushin分解一个分子S0、S1与S2的重组能，计算正常终止并检查无虚频后按照sob老师的博文进行操作并使用博文中编译好的linux版本dushin程序。其中对于S0-S1与S0-S2的重组能分解，dushin程序正常运行并获得结果，但是进行S1-S2的重组能分解时遇到forrtl：severe（24），报错情况如图。

(, 下载次数 Times of downloads: 40)

之后有另外尝试读取S1-S1与S2-S2均能正常读取运行。想请教一下各位老师，为什么会产生这种情况，有没有什么比较好的解决办法。谢谢。

作者
Author: 狄拉克自带buff 时间: 2020-10-16 11:17
同问，我查了一下所有相关报错问题的帖子，我怀疑原因是由于计算体系的分子量有点大(＞50)但是我没有证据。
请问您问题解决了吗?

作者
Author: abin 时间: 2020-10-16 11:20
你把你的所有需要用dushin处理的文件，
通过奶牛快传分享一下，
我试试我这里是否可以处理？

作者
Author: 起名真是难 时间: 2020-10-24 17:34
不好意思太久没关注这个帖子了。后来我换了论坛里比较新的版本之后设置第一行2 0 'curvilinear and atom reordering flags'解决了，但是后面有时候用的时候还是会有各种问题。现在还是不知道为什么。今天又遇到nv is incorrect的报错。。。

作者
Author: scfslyzkf 时间: 2020-10-29 21:12

起名真是难发表于 2020-10-24 17:34
不好意思太久没关注这个帖子了。后来我换了论坛里比较新的版本之后设置第一行2 0 'curvilinear and atom ...

感谢我最近也遇到了这个问题
第一行设置2 0 和第一行设置 1 2 有什么区别吗？

作者
Author: 起名真是难 时间: 2020-10-30 09:58
README里面是这么讲的
the curvilinear flag can be 0 for standard rectilinear analysis, 1 for normal curvilinear,
         or 2 for a more agressive curviilinear analysis of the Duschinsky matrix
      the atom ordering flag
   if = 0 if every set of coordinates has atoms provided in the same order and orientation
      if = 1 then the orientation and order is determined for the second dtata set, thence fixed
      if = 2 then the orientation and order is always determined
      if = 3 then the order is read from input
   a failure of 0, 1, or 2 leads to a read from input

作者
Author: zhangyl 时间: 2020-12-19 10:17

起名真是难发表于 2020-10-30 09:58
README里面是这么讲的
the curvilinear flag can be 0 for standard rectilinear analysis, 1 for normal ...

您好，我使用相同的两个结构(S0，T1结构）计算。发现当第一行设置为1 2时计算所得重组能较小（659 cm-1 467 cm-1）,当第一行设置为2 0时计算的重组能非常大（174373 cm-1 181929 cm-1）。请问此时应如何考虑呢？设置为2 0的计算结果还是否可用呢？

作者
Author: 王寓于 时间: 2020-12-23 18:21

zhangyl 发表于 2020-12-19 10:17
您好，我使用相同的两个结构(S0，T1结构）计算。发现当第一行设置为1 2时计算所得重组能较小（659 cm-1 4 ...

从README里看不同的设置是有差别的，不过究竟有什么差别我也不清楚。

作者
Author: wugaiyan 时间: 2021-6-30 16:52

起名真是难发表于 2020-10-24 17:34
不好意思太久没关注这个帖子了。后来我换了论坛里比较新的版本之后设置第一行2 0 'curvilinear and atom ...

你好，对于“nv is incorrect的报错”，解决了吗? 请问你是怎么解决的？

作者
Author: ridley 时间: 2022-5-5 14:59

起名真是难发表于 2020-10-24 17:34
不好意思太久没关注这个帖子了。后来我换了论坛里比较新的版本之后设置第一行2 0 'curvilinear and atom ...

您好，可否分享一下新版本的dushin，谢谢！

作者
Author: snljty 时间: 2022-5-5 15:47

wugaiyan 发表于 2021-6-30 16:52
你好，对于“nv is incorrect的报错”，解决了吗? 请问你是怎么解决的？

看的利用dushin计算空穴重组能时报错：nv is incorrect4楼。顺便可以参考基于Gaussian FCHT计算分解重组能和计算Huang-Rhys因子的Python小脚本直接使用Gaussian分解重组能。

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