| 查看: 1962 | 回复: 10 | ||||
哈哈宝木虫 (小有名气)
|
[求助]
用td方法做了个三态电子激发的模拟，结果出现了负值的光谱，求助！ 已有2人参与
|
|
输入命令： #p genecp td=(triplets,nstates=5,root=1) scrf=(PCM,solvent=dichloromethane) guess=check geom=check 输出文件部分： Excited State 1: Triplet-A -2.7255 eV -454.90 nm f=-0.0000 <S**2>=2.000 160 ->172 0.14488 160 ->175 -0.13306 160 ->188 -0.10636 162 ->175 -0.11445 163 ->172 -0.14890 163 ->175 -0.16744 164 ->175 -0.10638 165 ->175 -0.13703 166 ->172 0.16643 166 ->175 0.16901 168 ->172 -0.13219 168 ->175 0.12184 169 ->172 0.31913 169 ->175 -0.29000 This state for optimization and/or second-order correction. Total Energy, E(TD-HF/TD-KS) = -2389.28684805 Copying the excited state density for this state as the 1-particle RhoCI density. Excited State 2: Triplet-A -1.8956 eV -654.05 nm f=-0.0000 <S**2>=2.000 161 ->183 -0.17066 162 ->183 -0.12002 163 ->173 -0.10733 163 ->176 -0.16327 163 ->177 -0.11776 166 ->173 -0.11830 166 ->176 -0.12736 167 ->173 -0.11415 167 ->176 0.10424 167 ->180 0.12808 167 ->181 -0.22879 168 ->181 -0.10721 170 ->173 0.40676 170 ->174 0.12346 170 ->176 -0.15587 170 ->177 -0.15584 170 ->181 -0.11602 161 <-183 -0.10694 167 <-181 -0.13329 170 <-173 0.18909 Excited State 3: Triplet-A -1.8451 eV -671.98 nm f=-0.0000 <S**2>=2.000 161 ->184 -0.11372 162 ->184 0.14627 162 ->185 -0.10421 164 ->174 -0.15460 164 ->176 0.13147 164 ->177 -0.14912 165 ->174 0.15307 168 ->174 -0.11961 168 ->182 -0.24302 169 ->174 -0.12514 171 ->173 -0.12809 171 ->174 0.37880 171 ->176 0.14148 171 ->177 -0.16627 171 ->182 -0.13710 168 <-182 -0.14448 171 <-174 0.17861 Excited State 4: Triplet-A 0.9456 eV 1311.19 nm f=0.0000 <S**2>=2.000 161 ->173 0.29650 161 ->176 -0.11632 161 ->177 -0.11074 162 ->173 0.24071 163 ->176 0.27866 163 ->177 0.21614 163 ->179 -0.12910 166 ->173 0.18136 166 ->176 0.17070 166 ->177 0.11937 167 ->180 0.19101 167 ->181 -0.36194 167 ->182 -0.10393 167 ->183 0.17433 168 ->181 -0.16695 170 ->176 0.12976 170 ->181 -0.13366 170 ->183 -0.33513 170 ->185 -0.17527 161 <-173 0.22769 162 <-173 0.18429 163 <-176 0.21414 163 <-177 0.16679 163 <-179 -0.10067 166 <-173 0.13643 166 <-176 0.12972 167 <-180 0.14530 167 <-181 -0.27557 167 <-183 0.13582 168 <-181 -0.12683 170 <-181 -0.10108 170 <-183 -0.25980 170 <-185 -0.13645 Excited State 5: Triplet-A 1.0726 eV 1155.88 nm f=0.0000 <S**2>=2.000 149 ->224 -0.11200 156 ->174 -0.10270 161 ->174 0.20902 161 ->177 -0.10223 162 ->174 -0.20383 162 ->176 -0.16692 162 ->177 0.19561 164 ->174 -0.17557 164 ->176 0.18163 164 ->177 -0.21583 165 ->174 0.20906 165 ->182 0.15438 167 ->182 -0.17451 168 ->182 0.36812 168 ->184 0.11847 169 ->182 0.12214 169 ->184 0.10383 171 ->176 0.10308 171 ->177 -0.10203 171 ->182 0.20442 171 ->183 -0.17469 171 ->184 -0.25734 171 ->185 0.17982 161 <-174 0.15533 162 <-174 -0.15094 162 <-176 -0.12356 162 <-177 0.14568 164 <-174 -0.12847 164 <-176 0.13334 164 <-177 -0.15918 165 <-174 0.15195 165 <-182 0.11704 167 <-182 -0.12848 168 <-182 0.27063 171 <-182 0.14912 171 <-183 -0.13054 171 <-184 -0.19307 171 <-185 0.13507 请问各位木虫有没有遇见过这种情况，我看文献说的是在t1结构的基础上，用TD方法做激发态电子性质的模拟，难道是理解错误，求高手指教！ |
回复此楼» 收录本帖的淘帖专辑推荐
 经验 | 微观分子理论和模拟 |
» 猜你喜欢
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有109人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有1人回复
-
实际气体内压强的微观机制或将改写--------从一个基本的力学佯谬说起
已经有1人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有5人回复
-
招收0703开头的学硕调剂生 线上复试
已经有5人回复
yongleli
木虫 (正式写手)
- 应助: 93 (初中生)
- 金币: 3769.1
- 散金: 369
- 红花: 24
- 帖子: 764
- 在线: 101.3小时
- 虫号: 303595
- 注册: 2006-12-03
- 专业: 原子和分子物理
【答案】应助回帖
感谢参与，应助指数 +1
|
很久没用TD-DFT了，尝试说几个问题： 1. 态数过少，计算不准确。 一般TD-DFT计算要包括进40个以上的态，然后增加，比如40->80->100->200直到再继续增加态结果不变为止。这是因为TDDFT是个微扰法，理论上说来要把所有激发态全用上才能得到准确结果。（前提是激发只是对体系的一个微扰） 2.基态是否为能量低态？也就是说你这个体系有单态、三态，用你的泛函／基组组合，能得到跟实验／文献一致的能级顺序吗？（比如三态能量比单态低）否则你的TD-DFT计算相当于站在高能态往低能态跃迁，激发能肯定会是负值。 3.楼下继续。。。。。。 |
2楼2014-09-02 09:32:05
哈哈宝
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2253.3
- 红花: 2
- 帖子: 108
- 在线: 56.8小时
- 虫号: 1237673
- 注册: 2011-03-18
- 性别: MM
- 专业: 无机材料化学
3楼2014-09-02 09:36:03
枪下游魂
木虫 (著名写手)
- 应助: 718 (博后)
- 金币: 8336.9
- 散金: 175
- 红花: 80
- 帖子: 1498
- 在线: 542.6小时
- 虫号: 2620258
- 注册: 2013-08-28
- 专业: 理论和计算化学
4楼2014-09-02 09:48:58
yongleli
木虫 (正式写手)
- 应助: 93 (初中生)
- 金币: 3769.1
- 散金: 369
- 红花: 24
- 帖子: 764
- 在线: 101.3小时
- 虫号: 303595
- 注册: 2006-12-03
- 专业: 原子和分子物理
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
哈哈宝: 金币+10, ★★★很有帮助 2014-09-02 10:09:21
哈哈宝: 金币+10, ★★★很有帮助 2014-09-02 10:09:21
|
根据楼下 枪下游魂 说的， 你的体系首先要优化到合理结构。这个过程要做benchmark。 也就是说，选择流行的泛函／基组组合，做基态计算，然后再做激发态计算。 貌似你在这部分工作上也要努力一下。推荐你一个土耳其科学家（Ulrike Salzner）的文章，他在TD-DFT上下了很多功夫。建议你认真学习，多做计算。 http://pubs.acs.org/action/doSea ... qsSearchArea=author 尤其是这一篇： http://pubs.acs.org/doi/abs/10.1 ... p;searchHistoryKey= |
5楼2014-09-02 09:58:58
哈哈宝
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2253.3
- 红花: 2
- 帖子: 108
- 在线: 56.8小时
- 虫号: 1237673
- 注册: 2011-03-18
- 性别: MM
- 专业: 无机材料化学
6楼2014-09-02 10:02:38
yongleli
木虫 (正式写手)
- 应助: 93 (初中生)
- 金币: 3769.1
- 散金: 369
- 红花: 24
- 帖子: 764
- 在线: 101.3小时
- 虫号: 303595
- 注册: 2006-12-03
- 专业: 原子和分子物理
7楼2014-09-02 10:49:06
哈哈宝
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2253.3
- 红花: 2
- 帖子: 108
- 在线: 56.8小时
- 虫号: 1237673
- 注册: 2011-03-18
- 性别: MM
- 专业: 无机材料化学
8楼2014-09-02 11:05:03
aioroslove
金虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 1033
- 散金: 65
- 帖子: 215
- 在线: 145.1小时
- 虫号: 1199768
- 注册: 2011-02-07
- 性别: GG
- 专业: 理论和计算化学
9楼2014-09-02 19:11:35
小范范1989
木虫 (著名写手)
- 应助: 123 (高中生)
- 金币: 3821.3
- 散金: 968
- 红花: 243
- 帖子: 1979
- 在线: 452.8小时
- 虫号: 2013403
- 注册: 2012-09-19
- 性别: GG
- 专业: 原子和分子物理
10楼2014-09-02 20:59:18
10