24小时热门版块排行榜    

查看: 1424  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

lightme36

银虫 (初入文坛)

[交流] 【求助】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)去计算原来小基组优化好的结构就可以计算,但稍有变化。我想问为什么基组大了反而结果不那么精确?要怎么解决这个问题?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangmt

至尊木虫 (著名写手)

我叫MT

lightme36(金币+2):谢谢指导,不过我没找到您说的那个例子,不知道您能不能告诉我一下那个例子的题目 2010-06-29 21:13:43
引用回帖:
Originally posted by lightme36 at 2010-06-29 10:48:34:
下面就是出错时的log文件,我用B3LYP/6-31G(d,p)//B3LYP/6-31G(d,p)算加溶The electronic state of the initial guess is 1-A.
Requested convergence on RMS density matrix=1.00D-08 within 128 cycles.

溶液中不容易收敛是比较正常的,上面那一句就可以看出来,已经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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

coolrainbow

木虫 (著名写手)

未来国家冻凉


小木虫(金币+0.5):给个红包,谢谢回帖交流
lightme36(金币+4):谢谢,我金币也不多,待会我把出错的文件发上来 2010-06-29 10:37:42
这表明:小基组时有误差抵消现象,所以结果反而好,这是量化计算中的常见现象;也可能是结果误差正好和实验接近而已。具体结果你可能需要用更加大的方法,或者干脆就用3-21g

至于溶剂模型的问题,你最好给出具体的出错信息
技术博客:http://hi.baidu.com/coolrainbow/blog
2楼2010-06-28 20:26:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kekexiliwolf

金虫 (正式写手)

小屁虫


小木虫(金币+0.5):给个红包,谢谢回帖交流
lightme36(金币+2):谢谢,我试试 2010-06-29 10:38:14
保留geom=connectivity产生的冗余坐标,但是计算方法里面把geom=connectivity删除试试,我这样算的时候原来的报错没了
我们在为了文章而搞科研
3楼2010-06-29 10:12:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lightme36

银虫 (初入文坛)

引用回帖:
Originally posted by coolrainbow at 2010-06-28 20:26:11:
这表明:小基组时有误差抵消现象,所以结果反而好,这是量化计算中的常见现象;也可能是结果误差正好和实验接近而已。具体结果你可能需要用更加大的方法,或者干脆就用3-21g

至于溶剂模型的问题,你最好给出具 ...

下面就是出错时的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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 321求调剂 +5 大米饭! 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 283求调剂 +10 小楼。 2026-03-12 14/700 2026-03-16 16:08 by 13811244083
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 085600材料与化工 求调剂 +13 enenenhui 2026-03-13 14/700 2026-03-16 15:19 by 了了了了。。
[基金申请] NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗 20+5 NSFC2026我来了 2026-03-10 14/700 2026-03-15 23:53 by 不负韶华的虎
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 学硕285求调剂 +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[考研] 341求调剂 +4 番茄头--- 2026-03-10 4/200 2026-03-13 23:12 by JourneyLucky
[考研] 材料371求调剂 +9 鳄鱼? 2026-03-11 11/550 2026-03-13 22:53 by JourneyLucky
[考研] 求调剂(材料与化工327) +4 爱吃香菜啦 2026-03-11 4/200 2026-03-13 22:11 by JourneyLucky
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] 290求调剂 +9 ADT 2026-03-11 9/450 2026-03-13 21:55 by JourneyLucky
[考研] 315求调剂 +9 小羊小羊_ 2026-03-11 10/500 2026-03-13 21:13 by SXNU李老师
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[考研] 土木第一志愿276求调剂,科研和技能十分丰富,求新兴方向的导师收留 +3 土木小天才 2026-03-12 3/150 2026-03-13 15:01 by JourneyLucky
[考研] 277求调剂 +4 anchor17 2026-03-12 4/200 2026-03-13 11:15 by 白夜悠长
[考研] 求调剂 资源与环境 285 +3 未名考生 2026-03-10 3/150 2026-03-13 10:31 by houyaoxu
[考研] 333求调剂 +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
信息提示
请填处理意见