24小时热门版块排行榜    

查看: 2796  |  回复: 17

绍敏郡主

木虫 (著名写手)


[交流] 请教一个老问题:结构优化时的对称性问题

一个分子应该是D2h对称性的,我画好结构以后去计算,第一轮很快就能收敛,但是第二轮就不行了,而且能量差的极其夸张(从第一轮的-3000变为-4000),RMSDP也只在10(-2)数量级波动。
当我把初始结构稍微调一下(某两个原子只移动0.01埃),对称性变为Ci,用这个结构计算就正常了。
而文献里,不止一个人都说是用D2h对称性计算的。看来最后的稳定结构确实应该是D2h的。
请问我的哪里出问题了?
谢谢。

是因为初始结构离稳定结构差得太远了吗?

以上用的是DFT方法。后来发现用HF方法可以按D2h对称性正常优化结束。用这个优化的结构再用DFT方法计算,还是不能收敛。按说HF优化的结构已经比较接近稳定结构了吧,至少比手搭的好吧。。

[ Last edited by 绍敏郡主 on 2013-11-22 at 01:58 ]
回复此楼

» 收录本帖的淘帖专辑推荐

量子化学和计算化学

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

绍敏郡主

木虫 (著名写手)


又查了查文献,好像明白了。
据说要考虑相对论效应,D2h结构才是极小点,否则D2h结构是鞍点。
2楼2013-11-22 02:53:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

userhung

禁虫 (文学泰斗)



小木虫: 金币+0.5, 给个红包,谢谢回帖
blessing, blessing, blessing ! ~~~
3楼2013-11-22 06:17:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jinwei331

至尊木虫 (文坛精英)



小木虫: 金币+0.5, 给个红包,谢谢回帖
……………………
4楼2013-11-22 06:33:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
绍敏郡主: 金币+5, 谢谢 2013-11-22 18:01:44
你可以用文献中的泛函与基组验证一下
另外,不同的泛函会对结果产生一定的影响
5楼2013-11-22 17:53:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zindo

金虫 (小有名气)


★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
绍敏郡主: 金币+5, 谢谢 2013-11-23 10:56:09
同意楼上的,另外也可以考虑反其道而行之,先不限制优化(nosymm opt),看是否收敛,取其优化构型用Gaussview等工具限制成D2h,再用popt优化。
祝好运!
6楼2013-11-23 09:58:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
绍敏郡主: 金币+5, 谢谢 2013-11-23 20:41:03
有的dft有的时候结构确实不好,我自己见到的一个案例就是金属离子的配合物几种DFT给出不同的最优构型,有些结果从数据上判断都是错的。

楼上两层说的有道理,不知道杂化泛函情况如何。
7楼2013-11-23 18:37:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

绍敏郡主

木虫 (著名写手)


引用回帖:
7楼: Originally posted by 卡开发发 at 2013-11-23 18:37:32
有的dft有的时候结构确实不好,我自己见到的一个案例就是金属离子的配合物几种DFT给出不同的最优构型,有些结果从数据上判断都是错的。

楼上两层说的有道理,不知道杂化泛函情况如何。

谢谢。
用B3LYP和PBEPBE都不行,用HF可以。
10楼2013-11-23 20:43:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

virtualzx

木虫 (著名写手)


★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
绍敏郡主: 金币+5 2013-11-24 00:48:21
引用回帖:
10楼: Originally posted by 绍敏郡主 at 2013-11-23 07:43:24
谢谢。
用B3LYP和PBEPBE都不行,用HF可以。...

用的什么基组呢?
11楼2013-11-23 22:17:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

绍敏郡主

木虫 (著名写手)


引用回帖:
11楼: Originally posted by virtualzx at 2013-11-23 22:17:19
用的什么基组呢?...

C原子3-21g*,Y原子LanL2DZ或者CEP-4G
12楼2013-11-23 22:28:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

virtualzx

木虫 (著名写手)



小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
12楼: Originally posted by 绍敏郡主 at 2013-11-23 09:28:38
C原子3-21g*,Y原子LanL2DZ或者CEP-4G...

基组有点太小了吧。
13楼2013-11-24 00:33:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

绍敏郡主

木虫 (著名写手)


引用回帖:
13楼: Originally posted by virtualzx at 2013-11-24 00:33:05
基组有点太小了吧。...

我一般先用小基组优化一下,再用大基组优化一次,直接用大基组的话有点慢
14楼2013-11-24 00:46:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

virtualzx

木虫 (著名写手)



小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
14楼: Originally posted by 绍敏郡主 at 2013-11-23 11:46:58
我一般先用小基组优化一下,再用大基组优化一次,直接用大基组的话有点慢...

这样做没问题,但是就没法知道错误的结构到底是基组的问题还是方法的问题了。
15楼2013-11-24 01:03:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

绍敏郡主

木虫 (著名写手)


引用回帖:
15楼: Originally posted by virtualzx at 2013-11-24 01:03:22
这样做没问题,但是就没法知道错误的结构到底是基组的问题还是方法的问题了。...

嗯,有道理
16楼2013-11-24 01:06:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
10楼: Originally posted by 绍敏郡主 at 2013-11-23 20:43:24
谢谢。
用B3LYP和PBEPBE都不行,用HF可以。...

现在是scf不收敛了?还是几何优化振荡?

如果是scf不收敛建议试试采用b3lyp以及scf=verysleazy先优化,再恢复到正常精度,也许能成功。
17楼2013-11-24 01:15:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

绍敏郡主

木虫 (著名写手)


引用回帖:
17楼: Originally posted by 卡开发发 at 2013-11-24 01:15:00
现在是scf不收敛了?还是几何优化振荡?

如果是scf不收敛建议试试采用b3lyp以及scf=verysleazy先优化,再恢复到正常精度,也许能成功。...

是不收敛。
谢谢建议。
18楼2013-11-24 01:17:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
郑克宁8楼
2013-11-23 19:17   回复  
2013-11-23 19:17   回复  
相关版块跳转 我要订阅楼主 绍敏郡主 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 10/500 2026-08-29 16:49 by zICmwzsBXjbN
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 14:23 by jCd0dEvKHShX
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 14:13 by jCd0dEvKHShX
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +5 gdd2018 2026-08-28 10/500 2026-08-29 14:10 by Jacob678
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 11:49 by jCd0dEvKHShX
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +14 晴天加油 2026-08-26 15/750 2026-08-29 11:06 by biily
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见