24小时热门版块排行榜    

查看: 1426  |  回复: 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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 环境工程调剂 +8 大可digkids 2026-03-16 8/400 2026-03-18 09:36 by zhukairuo
[考研] 298-一志愿中国农业大学-求调剂 +4 手机用户 2026-03-17 4/200 2026-03-18 09:30 by b不易
[考研] 0703化学336分求调剂 +5 zbzihdhd 2026-03-15 6/300 2026-03-18 04:34 by anny19840123
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 能源材料化学课题组招收硕士研究生8-10名 +3 脱颖而出 2026-03-16 6/300 2026-03-17 21:19 by z1z2z3879
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
[考研] 341求调剂 +5 捣蛋猪猪 2026-03-11 7/350 2026-03-17 19:09 by 捣蛋猪猪
[考研] 考研化学学硕调剂,一志愿985 +4 张vvvv 2026-03-15 6/300 2026-03-17 17:15 by ruiyingmiao
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 材料专硕326求调剂 +6 墨煜姒莘 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[考研] 085600材料与化工 +4 安全上岸! 2026-03-16 4/200 2026-03-17 14:02 by 勇敢太监王公公
[考博] 26申博 +4 八旬速览 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 材料专硕274一志愿陕西师范大学求调剂 +5 薛云鹏 2026-03-13 5/250 2026-03-17 10:15 by Sammy2
[考研] 267一志愿南京工业大学0817化工求调剂 +6 SUICHILD 2026-03-12 6/300 2026-03-17 09:24 by 雾散后相遇lc
[考研] 311求调剂 +5 26研0 2026-03-15 5/250 2026-03-16 16:21 by a不易
[考研] 294求调剂 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
[考研] 255求调剂 +3 李嘉慧, 2026-03-12 4/200 2026-03-14 16:58 by 有只狸奴
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生,2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
信息提示
请填处理意见