计算化学公社

标题: ORCA_ESD在计算Displacing mode时候要算多少个？ [打印本页]

作者
Author: ABetaCarw 时间: 2023-5-17 10:19
标题: ORCA_ESD在计算Displacing mode时候要算多少个？
本帖最后由 ABetaCarw 于 2023-5-17 11:24 编辑

请教有ORCA_ESD使用经验的前辈。我之前没有使用ESD模块的经验，按照手册的例子，测试了BIA.xyz文件，研究磷光发射。每个子计算中需要跑30个Displacing mode
然后再官方论坛看到一位前辈算某体系的时候https://orcaforum.kofo.mpg.de/viewtopic.php?f=11&t=8736&p=37542 看到的是算了102个
现在我按照手册BIA的例子计算自己的体系，大约100个原子，现在到了算Displacing mode的阶段。但是这个算的很慢很慢，8个小时算了20个Displacing mode 。我不清楚我的这种体系该算几个Displacing mode，我现在是不是该继续等下去，还是果断换方法换服务器。

简单来说，问题其实就是：需要计算的Displacing mode数量和体系有什么关系？3N-6 吗？

我的体系是，某个金属有机化学物抽取的金属核心带外部配体的结构，总原子数100原子，含有小于5个过渡金属。我想要研究该结构的磷光发射谱图的主峰右边的小峰。

补充问题：
1 我在输出文件中看到了
Error (ORCA_ESD): The K*K value is too large (>30)!
然后我看了GS和T1的结构，存在某个类似于叔丁基旋转的区别。ORCA论坛说这时候谐振近似不适用，不清楚这个对我计算有什么影响，我需要怎么办比较好？
2 如果我尝试了删掉DOHT，然后就不算Displacing mode。这种是不是就是只考虑FC而不考虑HT。能否对我的体系定性正确研究？

作者
Author: wzkchem5 时间: 2023-5-17 17:03
大概3N左右或者6N左右，具体数字我记不清了，找个小体系算一下便知
K*K太大，可能（但不一定）导致高估振动分辨光谱的短波长区域的强度。一般的解决方法是用tcutfreq去掉频率较低的振动模式，因为一般只有低频振动模式的非谐振效应才有明显影响。tcutfreq一般惯例设为100cm-1左右，具体设置值没有客观判断标准，在让计算出来谱图看起来合理的情况下，越小越好。这个会给计算结果引入一些经验性，但是没有更好的办法，这是基于谐振近似的费米黄金规则计算方法的通病，不是ORCA的问题。
你研究的不是主峰的话，我觉得还是考虑H-T效应保险一些。H-T效应的主要影响在于让本来对称性禁阻的磷光发射能够有非零的强度，同时对振动精细结构的形状也有一定影响。

作者
Author: xuhj199508 时间: 2023-5-17 17:17
Displacing mode数量一般是原子数的3倍左右。

作者
Author: kevin00 时间: 2023-11-8 15:45
请问orca esd可以开启多核计算吗？

作者
Author: wzkchem5 时间: 2024-3-3 05:53

kevin00 发表于 2023-11-8 08:45
请问orca esd可以开启多核计算吗？

不能，原因参见http://bbs.keinsci.com/forum.php ... 9&fromuid=21811

作者
Author: kevin00 时间: 2024-7-4 22:06

wzkchem5 发表于 2024-3-3 05:53
不能，原因参见http://bbs.keinsci.com/forum.php?mod=redirect&goto=findpost&ptid=41113&pid=270019&fr ...

好的，感谢老师回复

作者
Author: GYJ2580 时间: 2025-7-25 09:05
你好，我在计算磷光速率时出现了Error (ORCA_ESD): The K*K value is too large (>30)! 随后尝试使用了tcutfreq 100发现还是会出现Error (ORCA_ESD): The K*K value is too large (>30)! 请问我是继续改变tcutfreq的值或者有什么更好的解决方法吗

作者
Author: wzkchem5 时间: 2025-7-25 19:56

GYJ2580 发表于 2025-7-25 09:05
你好，我在计算磷光速率时出现了Error (ORCA_ESD): The K*K value is too large (>30)! 随后尝试使用了t ...

考虑改用AHAS或VG

作者
Author: GYJ2580 时间: 2025-7-31 08:47

wzkchem5 发表于 2025-7-25 19:56
考虑改用AHAS或VG

（1）按ESD荧光的写法写一个输入文件，但在%tddft里写irootmult triplet（也就是让tddft算三重态，但让ESD以为算的是单重态）。这个任务会产生一个文件xxx.ES.hess，内含用VG方法估算的激发态结构和Hessian
（2）把xxx.ES.hess重命名为xxx.TS.hess，另写一个AH的输入文件，把这个Hessian读进来，让ESD以为这个Hessian是激发态平衡结构的激发态Hessian
你好老师，我就是在尝试上述方法，在第二步将xxx.TS.hess读取进来时还是出现了Error (ORCA_ESD): The K*K value is too large (>30)! 请问可能是什么造成的

作者
Author: wal 时间: 2025-7-31 10:37

GYJ2580 发表于 2025-7-31 08:47
（1）按ESD荧光的写法写一个输入文件，但在%tddft里写irootmult triplet（也就是让tddft算三重态，但让ES ...

ESD还能这么玩的？为什么不直接VG呢

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