| 查看: 1439 | 回复: 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)去计算原来小基组优化好的结构就可以计算，但稍有变化。我想问为什么基组大了反而结果不那么精确？要怎么解决这个问题？ |
回复此楼» 猜你喜欢
348 求调剂
已经有4人回复
-
招收 材料、新能源、冶金、化工、环境 等工科类考生
已经有2人回复
-
物理化学论文润色/翻译怎么收费?
已经有178人回复
-
招收 新能源、冶金、材料、化工、环境 等工科类考生
已经有0人回复
-
393求调剂推荐
已经有1人回复
-
能源环境新材料研究团队 招收 新能源、冶金、材料、化工、环境 等工科类考生
已经有0人回复
-
武汉纺织大学全国重点实验室杨应奎教授团队拟招17-19位硕博生！
已经有0人回复
-
材料与化工专硕接受跨专业调剂（重庆市属高校）
已经有0人回复
-
武汉纺织大学全国重点实验室杨应奎教授团队拟招17-19位硕博生！
已经有0人回复
-
能源环境新材料研究团队 招收 新能源、冶金、材料、化工、环境 等工科类研究生
已经有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
10