计算化学公社

 找回密码 Forget password
 注册 Register
Views: 18376|回复 Reply: 25
打印 Print 上一主题 Last thread 下一主题 Next thread

[量化理论] 关于OHF方法对CT结果的影响

[复制链接 Copy URL]

193

帖子

0

威望

4229

eV
积分
4422

Level 6 (一方通行)

各位老师好,今天学习了文献Huang, S.; Zhang, Q.; Shiota, Y.; Nakagawa, T.; Kuwabara,K.; Yoshizawa, K.; Adachi, C., J. Chem. TheoryComput. 2013, 9, 3872-3877.有一些关于OHF方法的问题不懂,想向各位老师特别是Sob老师请教:

1.没有看明白文中是怎么操作,我想弄清楚对于一个D-A复合体,找到最优HF%的,找到最优HF%后是否怎么在Gaussian程序中应用。还是说不需要找到optimal HF%, 直接利用文中公式计算E0-0(S1)和 E0-0(T1)?


2.文中最后对OHF方法和长程矫正密度泛函结果进行了比较,长程矫正密度泛函不可以,或者不需要优化HF%吗?如果长程矫正密度泛函不优化HF%,是不是对于CT体系普遍表现不如HF%优化后的普通杂化泛函好?






eureka

1294

帖子

0

威望

6236

eV
积分
7530

Level 6 (一方通行)

2#
发表于 Post on 2016-11-22 17:27:29 | 只看该作者 Only view this author
现在做延时荧光方向的这么多了?嘿嘿
之前论坛和群里有过讨论,你搜索一下。
PS:如果真要明白,这个文献看一两遍是不行的。我感觉至少要精读好几遍吧
https://www.x-mol.com/groups/fan_jianzhong

strive for greatness

193

帖子

0

威望

4229

eV
积分
4422

Level 6 (一方通行)

3#
 楼主 Author| 发表于 Post on 2016-11-22 17:38:06 | 只看该作者 Only view this author
小范范1989 发表于 2016-11-22 17:27
现在做延时荧光方向的这么多了?嘿嘿
之前论坛和群里有过讨论,你搜索一下。
PS:如果真要明白,这个文献 ...

谢谢,我查下。我主要还是做溶剂化理论的,只是研究体系是CT体系。
eureka

1102

帖子

18

威望

6643

eV
积分
8105

Level 6 (一方通行)

計算化学の社畜

4#
发表于 Post on 2016-11-22 17:45:35 | 只看该作者 Only view this author

Stand on the shoulders of giants

5万

帖子

99

威望

5万

eV
积分
112492

管理员

公社社长

5#
发表于 Post on 2016-11-22 18:02:47 | 只看该作者 Only view this author
1 4L已经回了
2 长程校正泛函没法说HF成份,要说w参数,看
优化长程校正泛函w参数的简便工具optDFTw
http://sobereva.com/346
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

193

帖子

0

威望

4229

eV
积分
4422

Level 6 (一方通行)

6#
 楼主 Author| 发表于 Post on 2016-11-22 21:06:11 | 只看该作者 Only view this author
sobereva 发表于 2016-11-22 18:02
1 4L已经回了
2 长程校正泛函没法说HF成份,要说w参数,看
优化长程校正泛函w参数的简便工具optDFTw

谢谢老师。还有两个问题向您请教:
1.计算q只有您幻灯片中的那个公式吗?如果对于分子间的电荷转移,能不能通过计算受体分子基态和激发态的NPA电荷总和的改变量来确定?
2.在找到B3LYP的OHF值后,怎么来改变B3LYP中默认的HF%,应该使用什么关键词?
eureka

202

帖子

0

威望

2663

eV
积分
2865

Level 5 (御坂)

7#
发表于 Post on 2016-11-22 21:55:39 | 只看该作者 Only view this author
本帖最后由 赵云跳槽 于 2016-11-22 22:00 编辑
Frank 发表于 2016-11-22 21:06
谢谢老师。还有两个问题向您请教:
1.计算q只有您幻灯片中的那个公式吗?如果对于分子间的电荷转移,能 ...


1.
q的计算的前提是主要是 HOMO → LUMO跃迁,计算过程如下
D:\Desktop\NPh3.fchk
8  // Orbital composition analysis
-1 // Define fragment 1
a 7 // Add all basis functions in atom 7(N atom) into fragment1 假设中间N为7号原子
q  // Save fragment and return to upper menu
3  // Orbital composition analysis with Ros-Schuit (SCPA) partition
65 // the HOMO orbital of NPh3
==>Composition of the fragment:   33.068321%
66 // the LUMO orbital of NPh3
==>Composition of the fragment:    0.744701%
则q=33.07%-0.74%=32.33% (NPh3)
当然,也可以定义除7外的所有原子为一个片段,计算后的q值也是一样的

至于后一问,是否存在关系不清楚

2.
优化得到OHF后,不是修改B3LYP泛函,而是选择最接近OHF的HF%成分的DFT泛函来计算


评分 Rate

参与人数
Participants 1
eV +3 收起 理由
Reason
sobereva + 3

查看全部评分 View all ratings

1102

帖子

18

威望

6643

eV
积分
8105

Level 6 (一方通行)

計算化学の社畜

8#
发表于 Post on 2016-11-23 09:22:14 | 只看该作者 Only view this author
赵云跳槽 发表于 2016-11-22 21:55
1.
q的计算的前提是主要是 HOMO → LUMO跃迁,计算过程如下
D:\Desktop\NPh3.fchk

好像是修改B3LYP的 HF%成分再计算的,参见4L的PPT
Stand on the shoulders of giants

193

帖子

0

威望

4229

eV
积分
4422

Level 6 (一方通行)

9#
 楼主 Author| 发表于 Post on 2016-11-23 09:49:47 | 只看该作者 Only view this author
赵云跳槽 发表于 2016-11-22 21:55
1.
q的计算的前提是主要是 HOMO → LUMO跃迁,计算过程如下
D:\Desktop\NPh3.fchk

谢谢老师。
eureka

202

帖子

0

威望

2663

eV
积分
2865

Level 5 (御坂)

10#
发表于 Post on 2016-11-23 09:51:37 | 只看该作者 Only view this author
冰释之川 发表于 2016-11-23 09:22
好像是修改B3LYP的 HF%成分再计算的,参见4L的PPT

不是,是根据最接近OHF的HF%成分的DFT泛函来计算的

可以很明显查看到SI文件中S1构型优化时是使用的是最接近OHF的HF%成分的DFT泛函

193

帖子

0

威望

4229

eV
积分
4422

Level 6 (一方通行)

11#
 楼主 Author| 发表于 Post on 2016-11-23 09:54:20 | 只看该作者 Only view this author
sobereva 发表于 2016-11-22 18:02
1 4L已经回了
2 长程校正泛函没法说HF成份,要说w参数,看
优化长程校正泛函w参数的简便工具optDFTw

老师,7L老师的回答已看懂。现在还是有些问题不清楚,想请教您:
1.对于分子间的电荷转移,能不能通过计算受体分子基态和激发态的NPA电荷总和的改变量来确定?
2.对于分子间的电荷转移,用OHF方法和优化range-separated 泛函w参数的方法哪个效果更好?还是说,对于每个研究体系不一样,没有统一的结论?
eureka

1102

帖子

18

威望

6643

eV
积分
8105

Level 6 (一方通行)

計算化学の社畜

12#
发表于 Post on 2016-11-23 09:58:10 | 只看该作者 Only view this author
Frank 发表于 2016-11-23 09:54
老师,7L老师的回答已看懂。现在还是有些问题不清楚,想请教您:
1.对于分子间的电荷转移,能不能通过计 ...

1. 这个你要实践一下才知道,可能一样,也可能会有系统误差
2. 这个也不好说,有的体系就算用最优化调控RSE泛函,算出来的结果也是不太理想(不光光是本文中提到的一些分子)
Stand on the shoulders of giants

193

帖子

0

威望

4229

eV
积分
4422

Level 6 (一方通行)

13#
 楼主 Author| 发表于 Post on 2016-11-23 10:07:22 | 只看该作者 Only view this author
冰释之川 发表于 2016-11-23 09:58
1. 这个你要实践一下才知道,可能一样,也可能会有系统误差
2. 这个也不好说,有的体系就算用最优化调控 ...

谢谢老师,请问如果主导CT跃迁的轨道不是HOMO和LUMO,那怎么计算q呢?
eureka

1102

帖子

18

威望

6643

eV
积分
8105

Level 6 (一方通行)

計算化学の社畜

14#
发表于 Post on 2016-11-23 10:14:25 | 只看该作者 Only view this author
本帖最后由 冰释之川 于 2016-11-23 10:31 编辑
赵云跳槽 发表于 2016-11-22 21:55
1.
q的计算的前提是主要是 HOMO → LUMO跃迁,计算过程如下
D:\Desktop\NPh3.fchk


这样的话,是不是sob老师的PPT描述有误了呢? 我们只要根据最优的HF%选择最接近的泛函来计算吸收能即可了,不需要刻意调节HF%至OHF值。
我的理解是,激发态的构型优化选择接近OHF值的泛函;然后在这个构型基础上用OHF值下的B3LYP算激发能
Stand on the shoulders of giants

202

帖子

0

威望

2663

eV
积分
2865

Level 5 (御坂)

15#
发表于 Post on 2016-11-23 10:30:46 | 只看该作者 Only view this author
冰释之川 发表于 2016-11-23 10:14
这样的话,是不是sob老师的PPT描述有误了呢? 我们只要根据最优的HF%选择最接近的泛函来计算吸收能即可 ...

恩,文献的SI中的数据中都是使用的是接近OHF的泛函,不是修改的B3LYP中HF成份

本版积分规则 Credits rule

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

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

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