计算化学公社

标题: BDE计算结果和文献相差一倍求助 [打印本页]

作者
Author: betsy 时间: 2022-5-9 14:49
标题: BDE计算结果和文献相差一倍求助
各位老师大家好！现根据一篇文献计算BDE，除了没有进行potential energy surface simulation, 所有的步骤以及基组都是都与文献一致，但是与文献中的值相差一倍，不知道是什么原因，下面是我的input file 以及log文件，麻烦各位老师帮忙看看，谢谢！ps: 输出文件后，按照Shermo教程，进行了单点能校正。文献link：Exceptionally stable blue phosphorescent organic light-emitting diodes | Nature Photonics

(, 下载次数 Times of downloads: 12)

作者
Author: wzkchem5 时间: 2022-5-9 15:15
不要用guess=mix，只有计算具有不止一个单电子，且各个单电子的自旋取向不一致（有的向上，有的向下）的时候才应该用guess=mix，其他情况下用guess=mix反倒是错误的。仔细阅读http://sobereva.com/82了解guess=mix的原理，只有了解了原理才能确保不用错。
另外检查一下电荷和自旋多重度，我不知道你每个输入文件对应的是哪个结构，但是我觉得你大概率设错了。对于每个体系（断键前、断键后），都要数一下有多少个alpha不成对电子，多少个beta不成对电子，然后用(alpha不成对电子数-beta不成对电子数+1）算自旋多重度。
以后贴图问问题，如果既贴了不止一个结构又贴了不止一个输入文件，一定要说清楚哪个输入文件对应哪个结构。

作者
Author: betsy 时间: 2022-5-9 16:32
本帖最后由 betsy 于 2022-5-9 16:59 编辑

wzkchem5 发表于 2022-5-9 15:15
不要用guess=mix，只有计算具有不止一个单电子，且各个单电子的自旋取向不一致（有的向上，有的向下）的时 ...

谢谢老师，图片的问题下次注意。第一个图片是断键前的input file，第二个和第三个是断键后的两个自由基（将C-N键断掉），第三个是断键前的分子构型（Si:1, C: 48, N:2,H:34 ）。我计算的未断键前的基态自旋量子数为1，断键后为2.

作者
Author: wzkchem5 时间: 2022-5-9 17:00

betsy 发表于 2022-5-9 09:32
谢谢老师，图片的问题下次注意。第一个图片是断键前的input file，第二个和第三个是断键后的两个自由基 ...

因为那两个数里第一个是电荷，第二个是自旋多重度。你这里涉及的分子都是中性的，所以电荷必须是0。
此外注意自旋多重度和自旋量子数是两个不同的概念，自旋多重度等于自旋量子数的2倍加1，高斯要求设的是自旋多重度

作者
Author: betsy 时间: 2022-5-10 14:46

wzkchem5 发表于 2022-5-9 17:00
因为那两个数里第一个是电荷，第二个是自旋多重度。你这里涉及的分子都是中性的，所以电荷必须是0。
此 ...

嗯嗯谢谢老师！

作者
Author: betsy 时间: 2022-5-11 11:17
本帖最后由 betsy 于 2022-5-11 11:22 编辑

各位老师，重新设置了自旋多重度并进行了计算，但是得到的结果分别与文献相差0.2与0.7ev。不知道出现这么大差别的原因在哪？文献中显示对于发光材料的计算基于基态闭壳层计算，采用了GD3的校正。对于host材料全部的模拟采用了非限制性计算，对于这一块儿我理解的应该是非限制性开壳层，且对于host材料没有进行校正。我计算的步骤（1）host材料基态优化；（2）材料加电荷优化，（3）因为知道断键的位置，直接对两个片段分别进行优化。ps在对片段优化的过程中分别尝试了片段1（正电荷）+片段2（无电荷）；片段1（无电荷）+片段2（正电荷）这两种组合方式，采用能量低的进行BDE计算。麻烦各位老师帮忙看下，谢谢！下面的截图分别未文献中的计算细节以及输入文件~ (, 下载次数 Times of downloads: 12) (, 下载次数 Times of downloads: 18)

作者
Author: wzkchem5 时间: 2022-5-11 13:59

betsy 发表于 2022-5-11 04:17
各位老师，重新设置了自旋多重度并进行了计算，但是得到的结果分别与文献相差0.2与0.7ev。不知道出现这么大 ...

先和文献的理论级别严格统一了再比。如果理论级别严格统一以后能对上，再研究理论级别有差别会导致什么情况

作者
Author: betsy 时间: 2022-5-11 15:15

wzkchem5 发表于 2022-5-11 13:59
先和文献的理论级别严格统一了再比。如果理论级别严格统一以后能对上，再研究理论级别有差别会导致什么情 ...

谢谢老师。关于host 材料的计算我采用的理论级别与文献是一致的，但是文献对于BD-1/2（发光材料）采用了GD3进行校正。可能对于host材料也采用了这个校正，目前正在尝试。关于这个理论级别严格统一这个问题，现在不知道怎么更正了，麻烦老师帮忙更正，谢谢

作者
Author: wzkchem5 时间: 2022-5-11 15:51

betsy 发表于 2022-5-11 08:15
谢谢老师。关于host 材料的计算我采用的理论级别与文献是一致的，但是文献对于BD-1/2（发光材料）采用了 ...

文献用的是B3LYP/6-311G**，你用的是PBE0/6-31G(d,p)。把泛函和基组都改成和文献一样，其他方面（溶剂模型等等）也和文献严格统一，再算。

作者
Author: betsy 时间: 2022-5-12 09:36

wzkchem5 发表于 2022-5-11 15:51
文献用的是B3LYP/6-311G**，你用的是PBE0/6-31G(d,p)。把泛函和基组都改成和文献一样，其他方面（溶剂模 ...

收到谢谢老师~

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