计算化学公社

 找回密码 Forget password
 注册 Register

看文献,说multiwfn可以计算CT态中,泛函的选择。不知道怎么回事,谢谢指点

查看数: 16570 | 评论数: 22 | 收藏 Add to favorites 5
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2015-3-22 10:30

正文摘要:

图片是文献中的说法,我用高斯计算了phcz分子想做测试,但是不知道他们是怎么计算出来的?高斯输入文件有log和fchk。谢谢指点

回复 Reply

石头 发表于 Post on 2016-7-12 08:02:04
sobereva 发表于 2016-7-12 06:54
加不加隐式溶剂模型对于几何结构优化的影响通常是微不足道的

明白了,感谢Sob老师。
sobereva 发表于 Post on 2016-7-12 06:54:50
石头 发表于 2016-7-11 21:36
您好,请问您知道为何基态没有加溶剂优化吗?这一点我不太明白。


加不加隐式溶剂模型对于几何结构优化的影响通常是微不足道的
石头 发表于 Post on 2016-7-11 21:36:06
小范范1989 发表于 2016-4-13 07:28
同时,你应该看文献知道,基态是bl3lyp gas下优化。 利用OHF计算其他的激发能或者什么的时候,都是在tolu ...

您好,请问您知道为何基态没有加溶剂优化吗?这一点我不太明白。
小范范1989 发表于 Post on 2016-4-15 08:04:48
yezhonghua 发表于 2016-4-14 22:44
那你的PHCZ S0 opt后,计算后q值能够为0.097么?我怎么是0.5996

sob老师给你回复了,就是他那样子
yezhonghua 发表于 Post on 2016-4-14 22:44:42
小范范1989 发表于 2016-4-14 22:25
确实是,不过就是找那种贡献比较明显的。我一般都是按照默认的阈值找,然后哪个上面有分布,就用哪个原子

那你的PHCZ S0 opt后,计算后q值能够为0.097么?我怎么是0.5996
yezhonghua 发表于 Post on 2016-4-14 22:35:17
小范范1989 发表于 2016-4-13 07:27
1:我是采用每个原子对HOMO/LUMO的贡献度,进行差值计算,并取绝对值。得到的数处以2就是q了。这样不对, ...


你好,我下载你的PHCZ分子,你只做了opt分析,没有进行freq计算。
我是通过SCPC方法找到HOMO轨道成分的原子贡献值ai;LUMO轨道的成分的原子贡献值bi。两者相减后取绝对值,但是得到的数值为0.5996和文献的0.097相差太大。不知道你是怎么求的q值。跪求一下方法。
小范范1989 发表于 Post on 2016-4-14 22:25:31
yezhonghua 发表于 2016-4-14 21:54
事实上基本上所有的原子都会对HOMO/LUMO有贡献,只是有的原子贡献很小占0.005%左右,这类原子也应该算在 ...

确实是,不过就是找那种贡献比较明显的。我一般都是按照默认的阈值找,然后哪个上面有分布,就用哪个原子
yezhonghua 发表于 Post on 2016-4-14 22:09:34
本帖最后由 yezhonghua 于 2016-4-14 22:17 编辑
yezhonghua 发表于 2016-4-13 20:46
望兄台指点一下,目前q值的计算卡住了。
      我还是不太明白你的意思。假如HOMO分布在1-20,21-2 ...
yezhonghua 发表于 Post on 2016-4-14 21:54:52
小范范1989 发表于 2016-4-13 07:27
1:我是采用每个原子对HOMO/LUMO的贡献度,进行差值计算,并取绝对值。得到的数处以2就是q了。这样不对, ...

事实上基本上所有的原子都会对HOMO/LUMO有贡献,只是有的原子贡献很小占0.005%左右,这类原子也应该算在里面吧。那这样的话片段怎么去定义呢?
谢谢!
yezhonghua 发表于 Post on 2016-4-14 20:12:19
本帖最后由 yezhonghua 于 2016-4-14 20:22 编辑
小范范1989 发表于 2016-4-13 22:21
,选择离21%相近的泛函—B3LYP进行激发态的opt,输出文件中计算得到q值。
这句话我不太懂。前面你是对的 ...

文献中PHCZ的计算在S0和S1态的q值是不一样的相差很大。
1、仔细看了文章,他应该是基于基态进行B3LYP计算得到的q(第一个q)值,TD计算EVA,通过不同的泛函得到EVA-OHF曲线,根据得到的q值在曲线上找到EVA。
2、计算EVE的时候是在甲苯中的,文献提到S1的激发态opt的泛函开始选取与基态q值相近的泛函,然后进行opt后,得到q值(第二个q)。最后,通过不同的泛函得到EVE-OHF曲线,根据q值找到EVE。
3、我选取了PXT-TRZ分子按照文献进行了计算,就是第一个q值计算为0.950(文献0.905)不太符合,所以想请教你如何计算q值的。
4、EVE计算得到的q(第二个q)也不是很确定,请教一下兄台。

这篇文章的计算量很大,我算了PXZ-TRZ分子,计算结果能对上,不管是EVA还是EVE,但是就是q的计算应该是有问题的。
小范范1989 发表于 Post on 2016-4-13 22:21:09
yezhonghua 发表于 2016-4-13 20:53
EVA在gas下计算的;EVE在甲苯中计算的。就是如何确定在甲苯中计算的q值。
我觉得EVE甲苯中的q值是不是这 ...

,选择离21%相近的泛函—B3LYP进行激发态的opt,输出文件中计算得到q值。
这句话我不太懂。前面你是对的。通过前面你不是确定出来了OHF了?然后用这个OHF来处理激发态,他没有用OHF来优化激发态,文献都是基于基态,来近似计算激发态。以为大的有机分子涉及到激发态的优化很费事,所以,他就基于基态,利用他里面的公式,求得S1和T1能量,找到了能差。
3didadi 发表于 Post on 2016-4-13 21:20:47
看你们讨论好热烈,搜了一下原文,还没有仔细看。但是这篇文章的Fig.4相关性那么高啊啊?看激发态计算的文章不多,但是也没见过这么好的结果呢。他定义的这个参数有什么特殊之处才达到这么好的效果?请各位老师讲解一下啊,非常感谢!! sob老师能否也大致介绍一下? @sobereva
yezhonghua 发表于 Post on 2016-4-13 20:53:53
小范范1989 发表于 2016-4-13 07:28
同时,你应该看文献知道,基态是bl3lyp gas下优化。 利用OHF计算其他的激发能或者什么的时候,都是在tolu ...

EVA在gas下计算的;EVE在甲苯中计算的。就是如何确定在甲苯中计算的q值。
我觉得EVE甲苯中的q值是不是这样计算的:选择和gas中 OHF相近的泛函,比如gas中计算的q=0.5,则泛函为0.5*42=21;那么EVE甲苯中的计算,选择离21%相近的泛函—B3LYP进行激发态的opt,输出文件中计算得到q值。
谢谢!
yezhonghua 发表于 Post on 2016-4-13 20:46:54
本帖最后由 yezhonghua 于 2016-4-14 20:17 编辑
小范范1989 发表于 2016-4-13 07:27
1:我是采用每个原子对HOMO/LUMO的贡献度,进行差值计算,并取绝对值。得到的数处以2就是q了。这样不对, ...

     望兄台指点一下,目前q值的计算卡住了。
      我还是不太明白你的意思。假如HOMO分布在1-20,21-24,28原子上,那么是否可以定义为1-20为一个片段;21-24为一个片度;28为一个原子。在选项-1里输入a 1-20原子定义为被考虑的片段,然后选3进入SCPA分析,再输入64,得到的Composition of the fragment就是a_i,再输入65,得到的就是b_i,a_i -b_i就是1-20原子的q。然后再计算21-24,28原子的q。同样,对于LUMO分布的原子也是执行上述操作。最后,所有的q相加就是最后的q了。
     谢谢!

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

GMT+8, 2025-8-14 08:06 , Processed in 0.372072 second(s), 25 queries , Gzip On.

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