| 查看: 1448 | 回复: 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)去计算原来小基组优化好的结构就可以计算，但稍有变化。我想问为什么基组大了反而结果不那么精确？要怎么解决这个问题？ |
回复此楼» 猜你喜欢
求ENDNOTE破解软件安装包
已经有4人回复
-
研究生招生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有166人回复
-
河海大学 （211.双一流高校）钙钛矿光电实验室丁勇团队欢迎有志青年！
已经有3人回复
-
黑龙江省原子能研究院化学课题组招收2026级化学专业或化学相关专业硕士生
已经有0人回复
-
山西大学资源与环境工程研究所化工分离/环境功能材料小组招收2026年博士生
已经有3人回复
-
26博士申请
已经有0人回复
-
一志愿南师大0703化学 275求调剂
已经有6人回复
-
A 区过线就可申！重庆双非公办高校 0856调剂，接受跨专业！纳米材料 / 润滑热门方向
已经有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
