zjxitcc 发表于 2022-6-6 16:12 麻烦老师了,实在是太感谢了,我大概理解了这里的轨道选取方法 |
本帖最后由 zjxitcc 于 2022-6-6 16:14 编辑 baobao123 发表于 2022-6-6 16:03 你可以看我5L更新后的回复,我列出了最重要的5个UNO轨道,这样在你自己动手前就有大概印象。这个体系比较简单,UNO就可以搞定,Molpro程序自身就支持UNO。当然,你用MOKIT自动做也可以。 |
zjxitcc 发表于 2022-6-6 15:19 好的老师,我这就试试MOKIT软件 |
wzkchem5 发表于 2022-6-6 15:36 好的谢谢老师 |
baobao123 发表于 2022-6-6 07:51 如果可能的话优先用楼上说的方法。如果因为某种原因必须在成键的结构上选轨道,那么需要你利用化学常识来推断文献选的轨道在成键以后分别变成了什么轨道。比如碳自由基的单占轨道和氧气的其中一个单占轨道成键以后,显然就会变成C-O的sigma和sigma*轨道。你可以同理推断氧气的另一个单占轨道变成了什么轨道。 |
本帖最后由 zjxitcc 于 2022-6-6 16:11 编辑 baobao123 发表于 2022-6-6 14:51 可以借助开源程序MOKIT来全自动/半自动确定活性轨道(可与Molpro联用,不冲突),看我举的几个例子 http://bbs.keinsci.com/thread-322-6-1.html 其中80层 http://bbs.keinsci.com/thread-17905-1-1.html “按理来说应该从产物自由基那端扫描,使得C-O键长变长,然后变成O2和环烷烃自由基来处理这个过程”这不是最妥当的做法,事实上,从断键扫至成键才是更合理的做法。因为断键时重要的UNO轨道、GVB轨道一目了然(就是HONO, LUNO或其附近),而成键情况下体系多组态/多参考态特征较小,UNO轨道中的HONO, LUNO(或其附近)很可能不是目标轨道,底下的次要UNO是离域的 就更难找了,此时需要用到 局域UNO轨道;GVB轨道一定包含你要的轨道,但其按占据数排序时 你得点好几次鼠标才能轮到、看到你要的轨道,多耗了十几秒。 我随手搭了一个结构,算了个UNO,发现这个体系比较简单,SOMO附近几个轨道就恰好满足文中所说(7e,5o),轨道顺序不需做任何调换,见下图
上面5个轨道即表示5o。其中32为单占据轨道,29~31为occupied UNO轨道,33为unoccupied UNO轨道,程序会默认识别为(7e,5o)。附件压缩包中是UNO轨道、GVB轨道和CASSCF(3,3)自然轨道。值得一提的是,MOKIT自动算出的结果为CASSCF(3e,3o),暗示并不需要(7e,5o)。当然了,想算多大活性空间、想取什么结果最终仍是用户决定的。
cyclohexane_O2.7z
(791.06 KB, 下载次数 Times of downloads: 40)
|
| 参与人数Participants 1 | eV +5 | 收起 理由Reason |
|---|---|---|
|
| + 5 | 牛! |
wzkchem5 发表于 2022-6-5 22:58 老师我画着看了一下还是比较懵,这里选取轨道的话有什么技巧吗 |
wzkchem5 发表于 2022-6-5 22:58 好的好的谢谢老师 |
| 先画轨道图,确认文献选取的活性轨道在R自由基和O2反应后分别变成了哪些轨道,再选择那些轨道。 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-18 18:11 , Processed in 0.247416 second(s), 26 queries , Gzip On.