版块导航
正在加载中...
客户端APP下载
登录
注册
帖子
帖子
用户
本版
应《网络安全法》要求,自2017年10月1日起,未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用,请尽快对帐号进行手机号验证,感谢您的理解与支持!
24小时热门版块排行榜
>
论坛更新日志
(434)
>
基金申请
(23)
>
导师招生
(14)
>
考博
(11)
>
职场人生
(9)
>
硕博家园
(9)
>
招聘信息布告栏
(8)
>
论文投稿
(5)
>
材料综合
(4)
>
虫友互识
(4)
>
休闲灌水
(3)
>
化工设备
(2)
>
教师之家
(2)
>
找工作
(2)
>
博后之家
(1)
>
绿色求助(高悬赏)
(1)
小木虫论坛-学术科研互动平台
»
计算模拟区
»
量子化学
»
Gaussian
»
请教一个老问题:结构优化时的对称性问题
18
1/1
返回列表
查看: 2450 | 回复: 17
只看楼主
@他人
存档
新回复提醒
(忽略)
收藏
在APP中查看
绍敏郡主
木虫
(著名写手)
应助: 289
(大学生)
金币: 1251.1
帖子: 1575
在线: 1866.7小时
虫号: 2611526
[交流]
请教一个老问题:结构优化时的对称性问题
一个分子应该是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
]
回复此楼
» 收录本帖的淘帖专辑推荐
量子化学和计算化学
» 猜你喜欢
钯炭、铂炭、钌炭等加氢催化剂
已经有27人回复
AgI的FT-IR谱图
已经有1人回复
物理化学论文润色/翻译怎么收费?
已经有87人回复
元素硫的化合价和氧化态
已经有1人回复
【急招】 北京工业大学 第二批博士申请-杰青课题组1-2名额,25日截止!!
已经有12人回复
1,6-己二胺甲基化合成四甲基己二胺
已经有2人回复
江苏科技大学2026年招聘启事
已经有0人回复
求助电子版书籍《半导体材料计算》
已经有0人回复
» 本主题相关商家推荐:
(我也要在这里推广)
高级回复
» 本主题相关价值贴推荐,对您同样有帮助:
怎样给出过渡态初始构型
已经有5人回复
gaussian保持对称性优化
已经有8人回复
如何让高斯在优化结构时保持其中某个单键不变
已经有7人回复
MS中对称性问题
已经有3人回复
结构优化 siesta 和 vasp
已经有7人回复
用castep进行结构优化的遇到的一些问题
已经有17人回复
请教一下那款免费软件的结构优化功能比较好用?
已经有7人回复
几何优化的对称性约束问题
已经有6人回复
晶体结构优化和画能带图的方法是否保持一致???
已经有12人回复
新手请教用NEB方法计算energy barrier.
已经有15人回复
分子对称性的选择问题
已经有4人回复
菜鸟问题:配合物分子式的书写 与 不对称结构单元 怎么看啊
已经有15人回复
结构优化中问题
已经有6人回复
关于结构优化中的对称性问题
已经有9人回复
QMolSymmetry:分子对称性识别,调整,优化(C++源码)
已经有6人回复
【求助】ISIF=3全局优化和对称性改变的问题
已经有7人回复
【求助】怎样解三方,六方,立方等对称性高的晶体结构?
已经有6人回复
【求助】用CASTEP计算晶体时结构问题
已经有10人回复
【求助】关于判断正反尖晶石结构的问题
已经有13人回复
【求助】Domol的优化问题
已经有10人回复
» 抢金币啦!回帖就可以得到:
查看全部散金贴
华南师范大学海洋环境科学课题组招聘青年英才和博士后
+
2
/120
【专业EPR】:自由基(超氧、羟基、单线态氧),空位缺陷;定量计算,VX: 761711562。
+
1
/89
Polymers期刊特刊Polymeric Materials for Advanced Drug Delivery组稿
+
1
/83
广东工业大学自动化学院国家特聘专家苏春翌教授招收2026年博士后及硕博研究生(推免)
+
1
/80
大连工业大学杰青/长江团队-生物质材料-储能电池方向招收2026级博士生
+
1
/34
海南大学国家高层次人才团队硕、博士招生
+
1
/32
湖南师范大学(211)—招收材料专业“申请-考核”制博士(有光学经验者优先)
+
1
/30
浙江理工大学国家级高层次人才团队招收能源光电催化材料方向博士生(截止5月31号)
+
2
/24
博士申请
+
1
/23
电力全国重点实验室双一流A类长江学者团队补招26年全日制学术博士1名(6月2日截止)
+
2
/22
上海大学微电子学院杨军教授团队招聘带编专任教师
+
1
/18
上海大学微电子学院杨军教授团队招聘带编专任教师
+
1
/17
军事医学研究院生物信息学方向胡明达老师团队招研究生 / 进修生 / 联培生
+
1
/10
上海大学微电子学院杨军教授团队招聘带编专任教师
+
1
/8
上海交通大学化学化工学院张智涛课题组诚聘博士后
+
1
/6
中科院博士后/特别研究助理招聘(光学工程、仪器科学、机械、电子、控制)
+
1
/4
电子科技大学材料学院SFT创新中心招收准备考硕和读博的科研助理 理工医交叉方向
+
1
/4
【通知】北京信息科技大学仪器科学与光电工程学院招聘优秀师资!
+
1
/4
上海交通大学海优建组招开山弟子
+
1
/1
上海理工大学,招收有机化学或者药物化学博士
+
1
/1
1楼
2013-11-21 20:16:21
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
绍敏郡主
木虫
(著名写手)
应助: 289
(大学生)
金币: 1251.1
帖子: 1575
在线: 1866.7小时
虫号: 2611526
又查了查文献,好像明白了。
据说要考虑相对论效应,D2h结构才是极小点,否则D2h结构是鞍点。
赞
一下
回复此楼
2楼
2013-11-22 02:53:04
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
userhung
禁虫
(文学泰斗)
应助: 1505
(讲师)
贵宾: 1.347
金币: 107941
帖子: 122982
在线: 4265.2小时
虫号: 119626
★
小木虫: 金币+0.5, 给个红包,谢谢回帖
blessing, blessing, blessing ! ~~~
赞
一下
回复此楼
3楼
2013-11-22 06:17:46
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
jinwei331
至尊木虫
(文坛精英)
应助: 3742
(副教授)
金币: 28683.7
帖子: 26959
在线: 1170.4小时
虫号: 2085922
★
小木虫: 金币+0.5, 给个红包,谢谢回帖
……………………
赞
一下
回复此楼
4楼
2013-11-22 06:33:01
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
sculhf
荣誉版主
(知名作家)
专家经验: +1
应助: 243
(大学生)
贵宾: 4.018
金币: 43335.3
帖子: 7682
在线: 2008小时
虫号: 497659
★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
绍敏郡主: 金币+5, 谢谢
2013-11-22 18:01:44
你可以用文献中的泛函与基组验证一下
另外,不同的泛函会对结果产生一定的影响
赞
一下
回复此楼
5楼
2013-11-22 17:53:19
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
zindo
金虫
(小有名气)
应助: 10
(幼儿园)
金币: 1700.4
帖子: 227
在线: 120.2小时
虫号: 47238
★ ★ ★ ★ ★ ★
小木虫: 金币+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的回帖
卡开发发
专家顾问
(著名写手)
专家经验: +224
应助: 47
(小学生)
金币: 3546.8
帖子: 2477
在线: 1254小时
虫号: 1369423
★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
绍敏郡主: 金币+5, 谢谢
2013-11-23 20:41:03
有的dft有的时候结构确实不好,我自己见到的一个案例就是金属离子的配合物几种DFT给出不同的最优构型,有些结果从数据上判断都是错的。
楼上两层说的有道理,不知道杂化泛函情况如何。
赞
一下
回复此楼
7楼
2013-11-23 18:37:32
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
绍敏郡主
木虫
(著名写手)
应助: 289
(大学生)
金币: 1251.1
帖子: 1575
在线: 1866.7小时
虫号: 2611526
引用回帖:
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
木虫
(著名写手)
QC强帖: 2
应助: 263
(大学生)
金币: 7161.3
帖子: 1605
在线: 317.6小时
虫号: 2069080
★ ★ ★ ★ ★ ★
小木虫: 金币+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的回帖
绍敏郡主
木虫
(著名写手)
应助: 289
(大学生)
金币: 1251.1
帖子: 1575
在线: 1866.7小时
虫号: 2611526
引用回帖:
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
木虫
(著名写手)
QC强帖: 2
应助: 263
(大学生)
金币: 7161.3
帖子: 1605
在线: 317.6小时
虫号: 2069080
★
小木虫: 金币+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的回帖
绍敏郡主
木虫
(著名写手)
应助: 289
(大学生)
金币: 1251.1
帖子: 1575
在线: 1866.7小时
虫号: 2611526
引用回帖:
13楼
:
Originally posted by
virtualzx
at 2013-11-24 00:33:05
基组有点太小了吧。...
我一般先用小基组优化一下,再用大基组优化一次,直接用大基组的话有点慢
赞
一下
回复此楼
14楼
2013-11-24 00:46:58
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
virtualzx
木虫
(著名写手)
QC强帖: 2
应助: 263
(大学生)
金币: 7161.3
帖子: 1605
在线: 317.6小时
虫号: 2069080
★
小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
14楼
:
Originally posted by
绍敏郡主
at 2013-11-23 11:46:58
我一般先用小基组优化一下,再用大基组优化一次,直接用大基组的话有点慢...
这样做没问题,但是就没法知道错误的结构到底是基组的问题还是方法的问题了。
赞
一下
回复此楼
15楼
2013-11-24 01:03:22
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
绍敏郡主
木虫
(著名写手)
应助: 289
(大学生)
金币: 1251.1
帖子: 1575
在线: 1866.7小时
虫号: 2611526
引用回帖:
15楼
:
Originally posted by
virtualzx
at 2013-11-24 01:03:22
这样做没问题,但是就没法知道错误的结构到底是基组的问题还是方法的问题了。...
嗯,有道理
回复此楼
16楼
2013-11-24 01:06:16
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
卡开发发
专家顾问
(著名写手)
专家经验: +224
应助: 47
(小学生)
金币: 3546.8
帖子: 2477
在线: 1254小时
虫号: 1369423
★
小木虫: 金币+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的回帖
绍敏郡主
木虫
(著名写手)
应助: 289
(大学生)
金币: 1251.1
帖子: 1575
在线: 1866.7小时
虫号: 2611526
引用回帖:
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
回复
花开半夏春暖
9楼
2013-11-23 19:17
回复
相关版块跳转
第一性原理
量子化学
计算模拟
分子模拟
仿真模拟
程序语言
我要订阅楼主
绍敏郡主
的主题更新
18
1/1
返回列表
如果回帖内容含有宣传信息,请如实选中。否则帐号将被全论坛禁言
普通表情
龙
兔
虎
猫
高级回复
(可上传附件)
百度网盘
|
360云盘
|
千易网盘
|
华为网盘
在新窗口页面中打开自己喜欢的网盘网站,将文件上传后,然后将下载链接复制到帖子内容中就可以了。
信息提示
关闭
请填处理意见
关闭
确定