计算化学公社

 找回密码 Forget password
 注册 Register

使用optDFTw优化w参数,设置不同上限得到不同的结果

查看数: 1046 | 评论数: 18 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2024-3-21 20:12

正文摘要:

本帖最后由 2421886719 于 2024-3-21 20:13 编辑 各位老师好,        我通过optDFTw优化wB97XD泛函的w参数,使用的分子为Y6,结构见下图,        于是我通过在 ...

回复 Reply

wzkchem5 发表于 Post on 2024-4-13 16:23:34
qian. 发表于 2024-4-12 07:00
老师,我算了一下发现两个激发能差的还挺大,w=0.01的激发能0.5ev,w=0的激发能0.9ev,但是审稿意见提了 ...

可以考虑在开着cLR的情况下调w值
qian. 发表于 Post on 2024-4-12 14:00:23
wzkchem5 发表于 2024-4-2 17:58
算TDDFT单点,看一下激发能就行了

老师,我算了一下发现两个激发能差的还挺大,w=0.01的激发能0.5ev,w=0的激发能0.9ev,但是审稿意见提了这个方法,我用optDFTw算的值很小,为0.01008,已经到了老师你说的超出合理范围(而且J^2也没到我设的收敛限0.001,但是已经到了我设的下限0.01,w=0.01008时J^2为0.03795),然后我用这个值算了一下紫外,与实验对不上,现在想请教一下老师接下来该怎么弄,我不知道该怎么弄了,麻烦老师了
qian. 发表于 Post on 2024-4-2 18:21:22
wzkchem5 发表于 2024-4-2 17:58
算TDDFT单点,看一下激发能就行了

好的,谢谢老师的答复,麻烦啦
qian. 发表于 Post on 2024-4-2 18:17:39
wzkchem5 发表于 2024-4-2 17:58
算TDDFT单点,看一下激发能就行了

好的,谢谢老师的解答,麻烦啦
wzkchem5 发表于 Post on 2024-4-2 17:58:17
qian. 发表于 2024-4-2 01:51
老师,我还想问一下,这个算是算什么呢,可以直接算单点吗,因为计算的是团簇,原子数比较多。
计算的关 ...

算TDDFT单点,看一下激发能就行了
qian. 发表于 Post on 2024-4-2 08:51:13
wzkchem5 发表于 2024-4-1 16:44
引用第一次提出range separate hybrid(印象中是CAM-B3LYP?)的原文,说w极其小相当于让泛函变成普通杂 ...

老师,我还想问一下,这个算是算什么呢,可以直接算单点吗,因为计算的是团簇,原子数比较多。
计算的关键词是#p lc-blyp/def2svp IOp(3/107=0010000000,3/108=0010000000)以及#p lc-blyp/def2svp IOp(3/107=0000000000,3/108=0000000000)吗?
wzkchem5 发表于 Post on 2024-4-1 16:44:51
qian. 发表于 2024-4-1 08:31
老师,你好,我想知道w低于0.01,不需要再调的这个说明有没有什么文献根据呢,想了解一下,麻烦了

引用第一次提出range separate hybrid(印象中是CAM-B3LYP?)的原文,说w极其小相当于让泛函变成普通杂化泛函了,就行了。
至于如何说明小于0.01(而不是比如说0.001)才算是足够接近普通杂化泛函,可以通过把w设成0做一个计算来说明,如果w=0的结果和w=0.01的结果几乎一样,就说明w=0.01是足够接近普通杂化泛函了。
qian. 发表于 Post on 2024-4-1 15:31:42
wzkchem5 发表于 2024-3-22 01:08
w=0.01即使客观上是J^2的极小值点,我觉得也超出合理范围了,因为这样调出来的泛函几乎就是一个普通(非范 ...

老师,你好,我想知道w低于0.01,不需要再调的这个说明有没有什么文献根据呢,想了解一下,麻烦了
2421886719 发表于 Post on 2024-3-22 21:44:53
KSeGaSn 发表于 2024-3-22 21:08
当然可以。首先,调节w的范围分离泛函本来就不是必须使用的方法,在溶剂条件下则更是要小心。
不过还是 ...

谢谢老师的提醒,之后会关注一下benchmark型的文献,我就是看了很多不是benchmark类型但是期刊档次不低的文章,他们用的都是调控w的wB97XD泛函计算重组能,但是没给出w值的具体参数,我用自己调控的w又不能复现文章结果
KSeGaSn 发表于 Post on 2024-3-22 21:08:13
2421886719 发表于 2024-3-22 20:22
谢谢老师的耐心答复,想问一下如果我使用B3LYP计算出的激发能和实验测试的紫外吸收光谱数据比较吻合,可 ...

当然可以。首先,调节w的范围分离泛函本来就不是必须使用的方法,在溶剂条件下则更是要小心。
不过还是看看你这类分子的既往工作、benchmark型的文献中用的较多的理论方法是什么。
2421886719 发表于 Post on 2024-3-22 20:24:49
wzkchem5 发表于 2024-3-22 01:08
w=0.01即使客观上是J^2的极小值点,我觉得也超出合理范围了,因为这样调出来的泛函几乎就是一个普通(非范 ...

感谢老师的耐心答复
2421886719 发表于 Post on 2024-3-22 20:22:31
KSeGaSn 发表于 2024-3-21 21:41
10.1016/j.surfin.2021.101385这篇文献里的结果吧,看起来他没明确说是在什么条件下优化的w参数,很有可能 ...

谢谢老师的耐心答复,想问一下如果我使用B3LYP计算出的激发能和实验测试的紫外吸收光谱数据比较吻合,可不可以使用B3LYP计算该分子的重组能
2421886719 发表于 Post on 2024-3-22 20:19:59
让你变成回忆 发表于 2024-3-21 20:33
文献中这么小的w大概率是在PCM下优化的。

谢谢老师的耐心答复,但是我通过在template文件中添加关键词"scrf(PCM, solvent=chloroform)"和"scrf(SMD, solvent=chloroform)"分别进行了优化,然后在终端以"./optDFTw 0.005 0.015 0.01 0.0001"命令进行优化,两次优化结果都是非常靠近下限0.005,远远达不到文献里说的0.01
KSeGaSn 发表于 Post on 2024-3-22 13:40:21
wzkchem5 发表于 2024-3-22 01:08
w=0.01即使客观上是J^2的极小值点,我觉得也超出合理范围了,因为这样调出来的泛函几乎就是一个普通(非范 ...

是这样的,这篇文献里用B3LYP或者PBE0计算出来的激发能已经足够好了,甚至用范围分离泛函调控的结果都称不上是正向优化……

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

GMT+8, 2024-11-25 03:02 , Processed in 0.169590 second(s), 25 queries , Gzip On.

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