应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】TDDFT计算中基组的问题
51/1返回列表
查看: 1428  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-06-28 19:33:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

coolrainbow

木虫 (著名写手)

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

至于溶剂模型的问题，你最好给出具体的出错信息
赞一下(1人) 回复此楼
技术博客：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删除试试，我这样算的时候原来的报错没了
赞一下(1人) 回复此楼
我们在为了文章而搞科研
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的回帖

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 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工 +5 安全上岸！ 2026-03-16 5/250 2026-03-18 15:33 by cmz0325
[考研] 266求调剂 +5 阳阳哇塞 2026-03-14 9/450 2026-03-18 15:05 by stone_128
[考研] 311求调剂 +6 26研0 2026-03-15 6/300 2026-03-18 14:43 by haxia
[考研] 化工学硕306求调剂 +10 42838695 2026-03-12 10/500 2026-03-18 14:42 by haxia
[考研] 0703化学调剂 ，六级已过，有科研经历 +10 曦熙兮 2026-03-15 10/500 2026-03-18 14:19 by 007_lilei
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[考研] 26考研求调剂 +6 丶宏Sir 2026-03-13 6/300 2026-03-17 16:13 by 醉在风里
[考研] 290求调剂 +6 孔志浩 2026-03-12 11/550 2026-03-17 14:41 by 周舟舟77
[考研] 一志愿南京大学，080500材料科学与工程，调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 085600调剂 +5 漾漾123sun 2026-03-12 6/300 2026-03-16 15:58 by 漾漾123sun
[考研] 285求调剂 +6 ytter 2026-03-12 6/300 2026-03-16 15:05 by njzyff
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 315求调剂 +9 小羊小羊_ 2026-03-11 10/500 2026-03-13 21:13 by SXNU李老师
[考研] 材料专硕350 求调剂 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 321求调剂（食品/专硕） +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
信息提示
关闭
请填处理意见
关闭