计算化学公社

标题: 用CDA分析过程中基组使用的疑问 [打印本页]

作者
Author: ch_pzq1102 时间: 2024-11-30 15:05
标题: 用CDA分析过程中基组使用的疑问
本人用multiwfn软件计算得到的CDA，对于HOMO轨道而言，片段1出现了超过100%的现象，不知道是否可取？望前辈们指点迷津。

体系为主族体系，小分子。
方法采用：b3lyp。基组分别用6-31g**,6-311g,6-311g**。以下是得到的结果图。

作者
Author: sobereva 时间: 2024-12-2 04:10
虚线是非占据轨道，箭头对应的是LUMO
Multiwfn的CDA分析用的是Mulliken方式做的轨道成分分析，超过100%很正常，尤其对于空轨道来说。当前情况当成100%就完了

作者
Author: ch_pzq1102 时间: 2024-12-2 10:34

sobereva 发表于 2024-12-2 04:10
虚线是非占据轨道，箭头对应的是LUMO
Multiwfn的CDA分析用的是Mulliken方式做的轨道成分分析，超过100%很 ...

好的，谢谢老师！

作者
Author: ZHUOHAO 时间: 2024-12-9 22:22

ch_pzq1102 发表于 2024-12-2 10:34
好的，谢谢老师！

你好，我最近也在做CDA分解。优化的时候也加了nosymm关键词得到相应的fch文件，但是在用multiwfn分析的时候，还是说坐标偏离。包括用了multiwfn安装软件下example\CDA\Pt(NH3)2Cl2去去算也有这个提示。请问你有什么建议解决吗？提示：Warning: The coordinate of the fragments deviates from the complex distinctly! The result may be fully meaningless. Please check input files of your quantum chemistry code to make the coordinate of the fragments fully consistent with the complex. If you are a Gaussian user, please do not forget to add "nosymm" to avoid automatically translating and rotating the system to standard orientation

作者
Author: lemon_electron 时间: 2024-12-10 01:34

ZHUOHAO 发表于 2024-12-9 22:22
你好，我最近也在做CDA分解。优化的时候也加了nosymm关键词得到相应的fch文件，但是在用multiwfn分析的时 ...

是把拆解后的片段又重新优化了吗？

作者
Author: 超限制抱怨 时间: 2024-12-10 09:41

sobereva 发表于 2024-12-2 04:10
虚线是非占据轨道，箭头对应的是LUMO
Multiwfn的CDA分析用的是Mulliken方式做的轨道成分分析，超过100%很 ...

想学习一下这个，sob老师来个论坛内链接呗，方便收藏

作者
Author: ZHUOHAO 时间: 2024-12-10 10:56

lemon_electron 发表于 2024-12-10 01:34
是把拆解后的片段又重新优化了吗？

对，拆解了重新算了一下但是没有opt

作者
Author: sobereva 时间: 2024-12-11 02:18

ZHUOHAO 发表于 2024-12-10 10:56
对，拆解了重新算了一下但是没有opt

算单个片段时必须带上nosymm，仔细看下文
使用Multiwfn做电荷分解分析(CDA)、绘制轨道相互作用图
http://sobereva.com/166
谈谈Gaussian中的对称性与nosymm关键词的使用
http://sobereva.com/297（http://bbs.keinsci.com/thread-1355-1-1.html）

作者
Author: ch_pzq1102 时间: 2024-12-11 15:03

ZHUOHAO 发表于 2024-12-9 22:22
你好，我最近也在做CDA分解。优化的时候也加了nosymm关键词得到相应的fch文件，但是在用multiwfn分析的时 ...

①片段是不能优化的，同时注意片段和整体的对应（参考sob老师的推文，通常来说直接从第一步优化后的结果提取就行）
②然后分别对Frag.1，Frag.2，complex只进行单点能计算就可以了（注意需要在同方法下进行）
本人还在初学，理解不一定对，仅供参考

。ps:主要为了抓紧发出去，拿折扣去参加培训课，这次的赶不上了可惜。。。。

作者
Author: ZHUOHAO 时间: 2024-12-11 17:38

sobereva 发表于 2024-12-11 02:18
算单个片段时必须带上nosymm，仔细看下文
使用Multiwfn做电荷分解分析(CDA)、绘制轨道相互作用图
http: ...

谢谢老师提醒，找到原因了。是我在把片段拖进Multiwfn时顺序错了，调整过来就解决了

作者
Author: ch_pzq1102 时间: 2024-12-11 19:27

ZHUOHAO 发表于 2024-12-9 22:22
你好，我最近也在做CDA分解。优化的时候也加了nosymm关键词得到相应的fch文件，但是在用multiwfn分析的时 ...

第一步优化后的复合物提取片段(具体参考sob老师的推文），注意片段编号对应
第二步分别对片段和复合物进行单点能计算（注意只需要计算单点能，且计算方法与第一步相同），获得chk/fch。然后就开始CDA分析吧。
然后现在给的贡献度相当乱，估计还要整理整理（目前在一条条对轨道。。）
本人初学，仅个人理解，多多交流，指点迷津。

ps:本来想早点投了，好参加sob老师的培训课，下次一定！

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