24小时热门版块排行榜    

查看: 1349  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 lightme36 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见