24小时热门版块排行榜    

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

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

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的回帖

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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 深圳大学硕士招生(2026秋,传感器方向,仅录取第一志愿) +4 xujiaoszu 2026-03-11 9/450 2026-03-17 10:29 by xujiaoszu
[考研] 278求调剂 +3 Yy7400 2026-03-13 3/150 2026-03-17 08:24 by laoshidan
[考研] 085601求调剂 +3 Du.11 2026-03-16 3/150 2026-03-16 20:42 by 无际的草原
[考研] 化学调剂0703 +8 啊我我的 2026-03-11 8/400 2026-03-16 17:23 by 我的船我的海
[考研] 0703化学调剂 ,六级已过,有科研经历 +7 曦熙兮 2026-03-15 7/350 2026-03-16 16:34 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
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[基金申请] 现在如何回避去年的某一个专家,不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 311求调剂 +5 牛乳糖的卡卡 2026-03-10 5/250 2026-03-14 00:05 by JourneyLucky
[考研] 285 求调剂 资源与环境 一志愿北京化工大学 +3 未名考生 2026-03-10 3/150 2026-03-13 23:04 by JourneyLucky
[考研] 求调剂(材料与化工327) +4 爱吃香菜啦 2026-03-11 4/200 2026-03-13 22:11 by JourneyLucky
[考研] 336求调剂 +6 Iuruoh 2026-03-11 6/300 2026-03-13 22:06 by JourneyLucky
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] 工科,求调剂 +3 我887 2026-03-11 3/150 2026-03-13 21:39 by JourneyLucky
[考研] 求调剂 +5 一定有学上- 2026-03-12 5/250 2026-03-13 18:31 by ms629
[考研] 【0856】化学工程(085602)313 分,本科学科评估A类院校化学工程与工艺,诚求调剂 +7 小刘快快上岸 2026-03-11 7/350 2026-03-13 16:06 by ruiyingmiao
[考研] 283求调剂,材料、化工皆可 +8 苏打水7777 2026-03-11 10/500 2026-03-13 09:06 by Linda Hu
信息提示
请填处理意见