24小时热门版块排行榜    

查看: 1420  |  回复: 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的回帖

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的回帖
查看全部 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的回帖

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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 275求调剂 +3 太阳花天天开心 2026-03-16 3/150 2026-03-17 08:01 by 热情沙漠
[考研] [复试调剂]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
[考研] 一志愿华中师范071000,325求调剂 +6 RuitingC 2026-03-12 6/300 2026-03-16 14:50 by 可淡不可忘
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +7 Ncdx123456 2026-03-13 8/400 2026-03-16 12:15 by karry wen
[考研] 311求调剂 +6 冬十三 2026-03-15 6/300 2026-03-16 08:00 by wang_dand
[考博] 欢迎申博同学联系 +3 天道酬勤2026686 2026-03-10 7/350 2026-03-15 19:03 by 天道酬勤2026686
[考研] 294求调剂 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[考研] 290求调剂 +4 @将就将就看 2026-03-10 8/400 2026-03-14 14:23 by 千千运气
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 0856材料与化工309分求调剂 +6 ZyZy…… 2026-03-10 6/300 2026-03-14 00:38 by JourneyLucky
[考研] 311求调剂 +8 zchqwer 2026-03-10 8/400 2026-03-14 00:01 by JourneyLucky
[考研] [0860]321分求调剂,ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] 一志愿西南交大,材料专硕317求调剂 +5 lx8568 2026-03-11 5/250 2026-03-13 21:43 by peike
[考研] 281求调剂 +9 Koxui 2026-03-12 11/550 2026-03-13 20:50 by Koxui
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
[考研] 【0856】化学工程(085602)313 分,本科学科评估A类院校化学工程与工艺,诚求调剂 +7 小刘快快上岸 2026-03-11 7/350 2026-03-13 16:06 by ruiyingmiao
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂! +3 不想理你 2026-03-12 3/150 2026-03-13 14:18 by JourneyLucky
信息提示
请填处理意见