| 查看: 1303 | 回复: 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)去计算原来小基组优化好的结构就可以计算，但稍有变化。我想问为什么基组大了反而结果不那么精确？要怎么解决这个问题？ |
回复此楼» 猜你喜欢
PbS量子点紫外吸收
已经有0人回复
-
北京交通大学王熙教授团队招收2026年入学催化方向博士研究生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有178人回复
-
祈福论文小修
已经有21人回复
-
香港理工大学计算材料与物理课题组招聘博后
已经有15人回复
-
香港理工大学计算材料与物理课题组招聘博后
已经有15人回复
-
西班牙Institut Català de Nanociència i Nanotecnologia, CSIC-BIST招收公派博士生
已经有0人回复
-
比利时鲁汶大学与国家留学基金委共同资助博士研究生CSC-KU Leuven PhD Scholarship
已经有0人回复
-
英国利物浦大学屠昕和雍雪课题组招收CSC博士生
已经有26人回复
-
北京交通大学王熙教授团队招收2026年入学催化方向博士研究生（还有1个名额）
已经有0人回复
-
华北电力大学（北京）（第一性原理计算）博士招生——学博，专博各1人
已经有44人回复
» 本主题相关价值贴推荐，对您同样有帮助:
高斯计算的问题？？
已经有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