dushin分解重组能forrtl：severe（24）求助

起名真是难

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

   

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

狄拉克自带buff

同问，我查了一下所有相关报错问题的帖子，我怀疑原因是由于计算体系的分子量有点大(＞50)但是我没有证据。
请问您问题解决了吗?

abin

你把你的所有需要用dushin处理的文件，
通过奶牛快传分享一下，
我试试我这里是否可以处理？

起名真是难

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

scfslyzkf

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

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

起名真是难

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

zhangyl

起名真是难 发表于 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的计算结果还是否可用呢？

王寓于

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

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

wugaiyan

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

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

ridley

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

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

snljty

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

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