作者indene (茚)
看板Chemistry
标题[实验] 理论计算之极化函数添加
时间Mon Nov 13 10:05:16 2017
各位版友安安
小弟最近再用gaussian 09跑有机发光小分子的模拟
常使用的方法是unrestriced DFT/B3LYP
基底函数是6-31++G(d,p)
我一般是跑含C H N O的分子(到第二周期)
现在有一个case 结构含有S 会到第三周期
我看很多paper是用6-31+G(d,p) 不会加更大的极化函数 或者加2个扩散函数
但找理论计算的文章 有人说要加到(3df,3pd)
Gaussian操作手册也这样写
所以我不懂到底用哪种才是对的
我也担心算太久 也不知道要用几个cpu算比较划算
总原子数大约50~60(氢也算进去) 有资源可以用御风者算
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 140.117.103.133
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Chemistry/M.1510538718.A.66E.html
1F：→ sexyrickysky: 48颗CPU，算50个transitions了不起一个礼拜就搞定了 11/13 11:50
那有硫原子是要加到哪里(d,p)够吗
我不太懂这含意以(d,p)来说 H要加到p轨域 而这个P轨域是指1p吗?
因为在gaussian view还可以加2p轨域 但1p轨域是什麽意思..
※ 编辑: indene (140.117.103.133), 11/13/2017 13:08:52
2F：推 sunev: basis set不一定要是原子的解，所以可以随便凑， 11/13 23:30
3F：→ sunev: 像是 R_{11}Y_{10}也是 okay的 11/13 23:31
4F：推 powernba: 我自己研究硫原子反应的时候也只用6-31+G(d) 11/14 16:52
5F：→ powernba: 不过我们老板有跟我们说如果反应过程中有牵扯到H，那使 11/14 16:53
6F：→ powernba: 用6-31+G(d,p)是比较适合的。 11/14 16:53
我只是要算萤光光谱跟吸收光谱 还有S1 T1能阶 所以不用用到(3df,3pd)吗?
※ 编辑: indene (218.173.156.89), 11/14/2017 19:36:09
7F：推 powernba: 我是觉得Gaussian算光谱本来就不太准，Basis set拉这麽 11/15 11:00
8F：→ powernba: 大，好像也没多大的意义...... 11/15 11:00
我不觉得 因为我有算一篇paper的萤光光谱 超级准
另外最近也开始有相关的paper算萤光光谱(以前没看过)不过他是用大基组
※ 编辑: indene (140.117.103.133), 11/15/2017 14:05:25
9F：推 sexyrickysky: 做任何研究之前，最基本的就是要做benchmark test 11/15 20:30
10F：→ sexyrickysky: 你现在要做的是，先找找近三年，多数人使用的方法， 11/15 20:32
11F：→ sexyrickysky: 并找寻是否有跟你系统有关的研究，而且有实验数据 11/15 20:32
12F：→ sexyrickysky: 用先前人家用过的方法算一遍，看看是否接近，如果有 11/15 20:33
13F：→ sexyrickysky: 所偏差，在试着加大基组。 11/15 20:34
14F：→ sexyrickysky: 基组要怎麽加，可以看一下几个LUMO+n (n=1~5)的轨域 11/15 20:36
15F：→ sexyrickysky: 分布。如果你看到相比於HOMO，LUMO的轨欲变得十分扩 11/15 20:37
16F：→ sexyrickysky: 散，或着LUMO轨欲有看到有涉及H原子的参与，或着不 11/15 20:37
17F：→ sexyrickysky: 硫原子中心的分布占大多数，表示要对硫原子与H原子 11/15 20:39
18F：→ sexyrickysky: 添加额外基组。至於要加diffuse or polarzied，如果 11/15 20:40
19F：→ sexyrickysky: 轨欲形状变化很多＝＝>添加极化函数。只是分布变广 11/15 20:41
20F：→ sexyrickysky: ==>加扩散函数。通常这两种都是要加的啦，只是谁 11/15 20:42
21F：→ sexyrickysky: 占的比重比高而已。 11/15 20:42
22F：→ sexyrickysky: 与其在这边纠结老半天，6-31+G(d,p)跟6-311G++(3df, 11/15 20:43
23F：→ sexyrickysky: 3dp)都算算看不就知道了。自己比较看看差异有没有 11/15 20:44
24F：→ sexyrickysky: 很多。再决定要用那个当作主要研究方法。 11/15 20:45
25F：→ sexyrickysky: 一开始的benchmark test本来就会耗掉很多时间，但是 11/15 20:46
26F：→ sexyrickysky: 可以省掉你後面很多功夫。不要在意消耗的SU跟时间， 11/15 20:47
27F：→ sexyrickysky: 你不自己算一次，永远不知道结果 11/15 20:48
非常谢谢您的回覆 因为我是第一次碰这个 实验室以前也没碰过
这方面的paper我看很多都是类似的方法(DFT/B3LYP 6-31G有些会加极化跟扩散)
最近有人用UDFT(但我也不知道为何发光材料要用这个)
虽然也有其他方法(PBEPBE WB97XD..等)
我是研究TADF材料(热延迟萤光) 发光过程是S0→...→S1→T1→S1→S0
我目前的LUMO非常集中在acceptor
但LUMO+1会往其他原子跑 +2又跑回来一点点 +3又跑出去
但都不会离分子很远 也没有跑到H
我想TADF主要涉及S1跟T1的交换 所以不太懂算LUMO+n的含意
不过我都会算^^
感恩 :)
※ 编辑: indene (218.173.145.91), 11/16/2017 12:15:46
28F：→ animadpig: 可以多参考Gaussian网页，有很多基本概念的讲解 11/16 20:10
谢谢你的建议，我目前会翻简体版的操作手册
29F：推 sexyrickysky: S1 or T1 excited state 不是单纯的HOMO到LUMO的tra 11/16 20:12
30F：→ sexyrickysky: nsition。而是可能涉及了HOMO-n到LUMO+n的transitio 11/16 20:12
31F：→ sexyrickysky: nal wave functions的linear combinations。所以你 11/16 20:12
32F：→ sexyrickysky: 用Gaussian 算吸收光谱，Gaussian 会给你一组数值， 11/16 20:12
33F：→ sexyrickysky: 告诉你S1 state,HOMO to LUMO到贡献占几%，HOMO-n 11/16 20:12
34F：→ sexyrickysky: 到LUMO+n占几%，这就是为什麽有时後我们要看一下LUM 11/16 20:12
35F：→ sexyrickysky: O+n轨域形状的原因。 11/16 20:12
原来如此 因为我看到大部分paper都是放HOMO LUMO 只有少部分会放+1-1的
您说的数值应该是Oscillator strength吧
另外不晓得您知不知道为何要用UDFT来算三重态 又该如何解释电子云分布的现象
(我也在小木虫找答案中@@)
谢谢^^
※ 编辑: indene (218.173.145.91), 11/16/2017 21:51:48
36F：推 sexyrickysky: 因为triplet state是open shell 组态，alpha跟beta 11/16 23:37
37F：→ sexyrickysky: 电子的数目不一样，意味着alpha电子跟beta电子受到 11/16 23:37
38F：→ sexyrickysky: 的interaction （eg. exchange) 会有所差异。所以再 11/16 23:37
39F：→ sexyrickysky: 处理上，会用两组MO各自处理alpha跟beta电子。而sin 11/16 23:37
40F：→ sexyrickysky: glet state因为alpha beta电子成对，所以用RDFT 11/16 23:37
感谢 大致理解
41F：推 sexyrickysky: Oscillator strength反应在光谱上的结果是吸收强度( 11/16 23:39
42F：→ sexyrickysky: intensity)。跟component transitions 不一样...... 11/16 23:39
抱歉 因为我只学过基本的的量子力学><(不是化学系)
Excited State 1: Singlet-A 3.2638 eV 379.88 nm f=0.0023
<S**2>=0.000
87 ->103 -0.20152 ← 您说的应该是这个
87 ->115 0.23170
88 ->102 -0.54675
88 ->108 -0.19811
88 ->112 0.11465
※ 编辑: indene (140.117.103.133), 11/17/2017 14:08:55