| 查看: 1418 | 回复: 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)去计算原来小基组优化好的结构就可以计算，但稍有变化。我想问为什么基组大了反而结果不那么精确？要怎么解决这个问题？ |
回复此楼» 猜你喜欢
汕头大学化学招调剂，相关内容如下。
已经有4人回复
-
国内青年破解费曼难题，从分子层面阐释液体蒸发（气液相变）微观过程
已经有5人回复
-
物理化学论文润色/翻译怎么收费?
已经有71人回复
-
研究发现一个新的临界系数，独立于临界压缩因子
已经有7人回复
-
200年来人们一直无法破解的难题： 气液临界点状态方程
已经有7人回复
-
关于引力延迟下双星系统角动量不守恒的疑惑---求高手解答！
已经有10人回复
-
七嗪类物质合成求助
已经有0人回复
-
半透性容器进行高湿或者低湿研究
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 高斯计算的问题？？
已经有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
20