计算化学公社

 找回密码 Forget password
 注册 Register

求助:高斯分别用M062X和B3LYP-GD3优化苯基对位氢迁移的产物结构差异较大,怎么解...

查看数: 1906 | 评论数: 18 | 收藏 Add to favorites 2
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2023-6-2 10:38

正文摘要:

本帖最后由 czoa 于 2023-6-2 17:59 编辑 请教各位老师,我想算图一的苯基对位氢迁移反应的能垒,最开始我使用的级别是CCSD(T)/aug-cc-pVTZ//M06-2X/6-31g(d,p),算出来的能垒和文献相符,但产物的相对能量比文献 ...

回复 Reply

czoa 发表于 Post on 2023-6-14 17:06:24
wjc404 发表于 2023-6-14 16:43
1楼给的产物结构,当自旋多重度为2时,有两种可能的电子结构,如下图所示:

可能M06-2X的优化任务中SCF ...

非常感谢您,请问您能分享一下从哪里看到的吗?
wjc404 发表于 Post on 2023-6-14 16:43:41
本帖最后由 wjc404 于 2023-6-23 22:48 编辑

1楼给的产物结构,当自旋多重度为2时,有两种可能的电子结构,如下图所示:


可能M06-2X的优化任务中SCF收敛到了上图右边的结构,而ωB97XD和B3LYP的优化任务中波函数收敛到了左边的结构。

06.23 续:
尝试去计算这两种结构,结果只找到了左边的那种,其FIC-NEVPT2(7,7)/def2-TZVP优化出的几何结构如下:
  1.   C        1.160        -0.785        -0.200
  2.   H        1.933        -1.238        -0.808
  3.   C        -0.055        -1.307        -0.007
  4.   C        -1.236        -0.708        0.202
  5.   H        -2.022        -1.098        0.836
  6.   C        -1.261        0.713        -0.233
  7.   C        -0.068        1.289        -0.088
  8.   H        -2.170        1.251        -0.467
  9.   C        1.266        0.670        0.224
  10.   H        1.517        0.852        1.274
  11.   H        2.049        1.140        -0.374
复制代码
以上结构的(7e,7o)活性空间自然轨道见下图:



感谢ORCA5和Multiwfn软件。


czoa 发表于 Post on 2023-6-14 08:52:13
wzkchem5 发表于 2023-6-13 21:54
随便找一本有机化学教材,都会讲怎么画共振式。这样你可以知道有多少个单电子,但是不知道几个单电子自旋 ...

好的,感谢王老师的回复。
wzkchem5 发表于 Post on 2023-6-13 21:54:51
czoa 发表于 2023-6-13 14:30
感谢老师的回复,学生对于自旋多重度的判断来自于网上的零散学习,目前知道是alpha电子数-beta电子数+1, ...

随便找一本有机化学教材,都会讲怎么画共振式。这样你可以知道有多少个单电子,但是不知道几个单电子自旋向上,几个自旋向下。然后再按Hund规则判断自旋向上和向下各有几个。Hund规则极个别时候有例外,但这个可以暂时不用管。
假如要算金属配合物,尤其是过渡金属配合物,那么还需要学配体场理论
wzkchem5 发表于 Post on 2023-6-13 21:14:24
czoa 发表于 2023-6-13 13:28
感谢老师,也就是说对我这个体系自旋多重度是偶数吗?

对。画出最主要的共振式,数一下有几个单电子,其中几个自旋向上,几个自旋向下,两者相减再加1就是自旋多重度。这个规则有一些例外,但是做普通的DFT计算一般不用管
czoa 发表于 Post on 2023-6-13 20:28:23
wzkchem5 发表于 2023-6-13 20:09
自己学一下怎么判断一个体系的自旋多重度,就知道怎么回事了。这是最基本的常识

感谢老师,也就是说对我这个体系自旋多重度是偶数吗?
wzkchem5 发表于 Post on 2023-6-13 20:09:35
czoa 发表于 2023-6-13 12:50
感谢sob老师,我对产物结构设置0,1 使用6-311+g(d,p) um062x guess=mix nosymm stable关键词测试波函数 ...

自己学一下怎么判断一个体系的自旋多重度,就知道怎么回事了。这是最基本的常识
czoa 发表于 Post on 2023-6-13 19:50:33
sobereva 发表于 2023-6-6 09:38
stable关键词测试一下便知

感谢sob老师,我对产物结构设置0,1 使用6-311+g(d,p) um062x guess=mix nosymm stable关键词测试波函数稳定性时,程序直接报对我这个体系(C6H5):The combination of multiplicity 1 and    41 electrons is impossible.就停了
sobereva 发表于 Post on 2023-6-6 09:38:26
czoa 发表于 2023-6-5 09:01
sob老师不好意思,再打扰一下,感觉我这个反应从电子数来讲,似乎不太有单重态的可能。

stable关键词测试一下便知
sobereva 发表于 Post on 2023-6-4 02:48:43
czoa 发表于 2023-6-3 15:56
感谢sob老师,我按照您的建议,将wb97xd优化好的结构和波函数作为M06-2X的初猜,最后算出的产物能量与文献 ...

这种体系大概率是双自由基体系,当闭壳层计算本来就不对
仔细看
谈谈片段组合波函数与自旋极化单重态
http://sobereva.com/82

每个级别计算的时候都要确保波函数的稳定性
czoa 发表于 Post on 2023-6-3 16:03:14
sobereva 发表于 2023-6-3 11:22
调控w是算激发态、超极化率才用得着的,跟你的情况没关系

sob老师不好意思,我本来是想回复的时候再发一个附件问您的,不小心弄成发帖了,请老师帮我看看我的附件。再次感谢。

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-25 08:22 , Processed in 0.174052 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list