| 查看: 1419 | 回复: 4 | |||
[交流]
【求助】TDDFT计算中基组的问题 已有4人参与
|
| 我由于计算一个比较大的天然化合物ECD谱，加上氢原子大概有80多个原子。我先用B3LYP/3-21G//B3LYP/3-21G计算，得到的结果和实验有点类似。于是我想用B3LYP/6-31G(d,p)//B3LYP/6-31G(d,p)计算得更加精确些，但是结果反而变得很差，和实验值差得很远，并且如果加入cosmor溶剂模型的话，在计算ECD谱的时候就出错，而我用B3LYP/6-31G(d,p)去计算原来小基组优化好的结构就可以计算，但稍有变化。我想问为什么基组大了反而结果不那么精确？要怎么解决这个问题？ |
回复此楼» 猜你喜欢
cif求助下载
已经有1人回复
-
高效微纳米储能材料与产业化推进实验室 苏州 招收硕士研究生（含调剂）
已经有5人回复
-
物理化学论文润色/翻译怎么收费?
已经有158人回复
-
高效微纳米储能材料与产业化推进实验室 苏州 招收硕士研究生（含调剂）
已经有0人回复
-
高效微纳米储能材料与产业化推进实验室 苏州 招收硕士研究生（含调剂）
已经有2人回复
-
化学化工学院 招收化工、化学、材料等相关方向研究生（学硕、专硕都有调剂名额）
已经有1人回复
-
高效微纳米储能材料与产业化推进实验室 苏州 招收硕士研究生（含调剂）
已经有0人回复
-
温州大学新能源材料微观溯源课题组招收2026级博士研究生（材料、化学方向）
已经有6人回复
-
材料与化工 323 英一+数二+物化，一志愿：哈工大 本人本科双一流
已经有5人回复
-
26申博
已经有2人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 高斯计算的问题？？
已经有5人回复
- Co的计算问题！新手求助
已经有10人回复
- Fluent计算时的迭代问题
已经有14人回复
- 求助环糊精包合小分子的构象（从环糊精的哪个面进去问题）， DFT计算判断
已经有3人回复
- 求教Gaussian用TDDFT计算单重态多重态的问题
已经有6人回复
- 求助一个体积计算问题
已经有10人回复
- 高斯计算中分子轨道问题
已经有6人回复
- SPSS19计算LD50问题
已经有9人回复
- 抗寒胁迫下，做转录组测序的取样部位问题。
已经有6人回复
- 【求助】求助关于TDDFT计算旋光度的问题
已经有7人回复
- 【求助】问如何选择基组精确计算有机小分子结构和能量的问题
已经有5人回复
- 【求助】咨询单点能计算采用优化计算时得出的结构和基组不能重复结果的问题
已经有3人回复
- 【求助】关于TDDFT计算激发态的问题
已经有7人回复
coolrainbow
木虫 (著名写手)
未来国家冻凉
- QC强帖: 11
- 应助: 1 (幼儿园)
- 贵宾: 0.1
- 金币: 3979.5
- 散金: 10
- 红花: 50
- 帖子: 1137
- 在线: 272.3小时
- 虫号: 239560
- 注册: 2006-04-07
- 性别: GG
- 专业: 理论和计算化学
2楼2010-06-28 20:26:11
kekexiliwolf
金虫 (正式写手)
小屁虫
- 应助: 14 (小学生)
- 金币: 1152.3
- 散金: 16
- 红花: 1
- 沙发: 1
- 帖子: 533
- 在线: 160.6小时
- 虫号: 928533
- 注册: 2009-12-15
- 性别: GG
- 专业: 药物设计与药物信息
3楼2010-06-29 10:12:10
|
下面就是出错时的log文件，我用B3LYP/6-31G(d,p)//B3LYP/6-31G(d,p)算加溶剂模型时ecd的错误，空气中也可以计算，结果很差。同时B3LYP/6-31G(d,p)//B3LYP/3-21G加溶剂化，这样也可以计算，结果尚可。 The electronic state of the initial guess is 1-A. Requested convergence on RMS density matrix=1.00D-08 within 128 cycles. Requested convergence on MAX density matrix=1.00D-06. Requested convergence on energy=1.00D-06. No special actions if energy rises. QIter: Convergence failure for iQ= 1! Last iteration charges are: 1 1-0.156939D-01 2-0.162575D-01 3-0.969792D-01 4-0.977627D-01 5-0.675217D-01 6-0.108387D+00 .............................. 4867-0.130186D+00 4868-0.439563D-01 4869-0.408914D-01 4870-0.719997D-01 Error termination via Lnk1e in C:\G03W\l502.exe at Fri Jun 11 14:16:53 2010. Job cpu time: 0 days 0 hours 15 minutes 48.0 seconds. File lengths (MBytes): RWF= 231 Int= 0 D2E= 0 Chk= 1 Scr= 1 |
4楼2010-06-29 10:48:34
zhangmt
至尊木虫 (著名写手)
我叫MT
- QC强帖: 5
- 应助: 99 (初中生)
- 金币: 6961.8
- 散金: 10406
- 红花: 49
- 帖子: 1761
- 在线: 763.2小时
- 虫号: 880392
- 注册: 2009-10-22
- 性别: GG
- 专业: 理论和计算化学
lightme36(金币+2):谢谢指导，不过我没找到您说的那个例子，不知道您能不能告诉我一下那个例子的题目 2010-06-29 21:13:43
|
溶液中不容易收敛是比较正常的，上面那一句就可以看出来，已经128cycles了，对于64位的gaussian来说，默认的502如果128cycles不收敛，就报错终止。 你应该用molden监控一下SCF的收敛过程，看看是振荡了，有收敛趋势，还是没有收敛的意思，然后分别采取不同的措施。本版有我以前写的一个例子，带图示，有解释。 试试这个地址 http://muchong.com/bbs/viewthread.php?tid=1668181 [ Last edited by zhangmt on 2010-6-29 at 22:44 ] |
5楼2010-06-29 12:40:50
