计算化学公社

 找回密码 Forget password
 注册 Register
Views: 266|回复 Reply: 11
打印 Print 上一主题 Last thread 下一主题 Next thread

[Gaussian/gview] 结构优化能量收敛且没有虚频但是优化结构与初始构型严重跑偏

[复制链接 Copy URL]

10

帖子

0

威望

63

eV
积分
73

Level 2 能力者

各位老师好,我做的是一个四核物质的几何优化,优化结果显示能量收敛以及没有虚频 但是结构严重不合理 想问问各位老师这种情况的结果是可靠的吗
所用的基组是b3lyp/genecp Cu用的是SDD 其余原子是6-31g*

下面是该结构的输入文件的构型是


输出文件的结果是


下面是输入跟输出文件
fcore.gjf (9.37 KB, 下载次数 Times of downloads: 4) fcore.log.bz2 (1.21 MB, 下载次数 Times of downloads: 0)

306

帖子

0

威望

688

eV
积分
994

Level 4 (黑子)

2#
发表于 Post on 2026-1-27 22:26:13 | 只看该作者 Only view this author
B3LYP泛函用在这里不合适吧,可以用TPSS或PBE0等,参考简谈量子化学计算中DFT泛函的选择

不过这里我也没有看出来“结构严重不合理”在哪

115

帖子

0

威望

471

eV
积分
586

Level 4 (黑子)

3#
发表于 Post on 2026-1-27 22:37:28 | 只看该作者 Only view this author
KazusaT 发表于 2026-1-27 22:26
B3LYP泛函用在这里不合适吧,可以用TPSS或PBE0等,参考简谈量子化学计算中DFT泛函的选择

不过这里我也没 ...

B3LYP-D3(BJ)的确不是最佳,但一般来说也不是不能接受;PBE0-D3(BJ)或TPSSh-D3(BJ)得到的结果肯定更合理,但多数情况也不至于总体几何结构上与B3LYP-D3(BJ)优化的结果相差太显著

115

帖子

0

威望

471

eV
积分
586

Level 4 (黑子)

4#
发表于 Post on 2026-1-27 22:40:58 | 只看该作者 Only view this author
本帖最后由 UW_0728. 于 2026-1-27 22:44 编辑

光从计算级别看没什么特别不合理的地方,优化成啥就是啥(好奇你口中“结构严重不合理”具体是指什么地方);不放心的话按2L建议的改用PBE0-D3(BJ)优化下看看会怎么样

还有就是你opt关键词里咋那么冗长,最不理解的是干嘛opt里面加个noraman?

10

帖子

0

威望

63

eV
积分
73

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2026-1-27 22:46:07 | 只看该作者 Only view this author
KazusaT 发表于 2026-1-27 22:26
B3LYP泛函用在这里不合适吧,可以用TPSS或PBE0等,参考简谈量子化学计算中DFT泛函的选择

不过这里我也没 ...

好的 谢谢

10

帖子

0

威望

63

eV
积分
73

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2026-1-27 22:48:53 | 只看该作者 Only view this author
UW_0728. 发表于 2026-1-27 22:40
光从计算级别看没什么特别不合理的地方,优化成啥就是啥(好奇你口中“结构严重不合理”具体是指什么地方) ...

你好 是这样的 因为初始结构可以看出是对称的 那么就是存在两类Cu原子 然后优化之后感觉是不对称了 后期要进行反应的话那么就是四类Cu 这里是不清楚是显示方式的问题还是说就是存在四类Cu了呢

115

帖子

0

威望

471

eV
积分
586

Level 4 (黑子)

7#
发表于 Post on 2026-1-27 22:56:05 | 只看该作者 Only view this author
TanJing 发表于 2026-1-27 22:48
你好 是这样的 因为初始结构可以看出是对称的 那么就是存在两类Cu原子 然后优化之后感觉是不对称了 后期 ...

改用PBE0-D3(BJ)重新优化一下看看吧。另外把没用的冗余关键词自觉去掉

10

帖子

0

威望

63

eV
积分
73

Level 2 能力者

8#
 楼主 Author| 发表于 Post on 2026-1-27 23:00:23 | 只看该作者 Only view this author
UW_0728. 发表于 2026-1-27 22:56
改用PBE0-D3(BJ)重新优化一下看看吧。另外把没用的冗余关键词自觉去掉

好的 谢谢老师

306

帖子

0

威望

688

eV
积分
994

Level 4 (黑子)

9#
发表于 Post on 2026-1-28 00:12:48 | 只看该作者 Only view this author
TanJing 发表于 2026-1-27 22:48
你好 是这样的 因为初始结构可以看出是对称的 那么就是存在两类Cu原子 然后优化之后感觉是不对称了 后期 ...

这里没有理由说低对称的结构就是不合理的,结构上没有明显不符合化学直觉的地方
如果你有实验证据(光谱、晶体)确信是两类Cu那可以考虑限制性优化,没有其他证据的话优化出来是什么样子就是什么样子
当然你也可以考虑做个简单的构象搜索

10

帖子

0

威望

63

eV
积分
73

Level 2 能力者

10#
 楼主 Author| 发表于 Post on 2026-1-28 14:46:01 | 只看该作者 Only view this author
KazusaT 发表于 2026-1-28 00:12
这里没有理由说低对称的结构就是不合理的,结构上没有明显不符合化学直觉的地方
如果你有实验证据(光谱 ...

好的 谢谢大佬 这里实验结果显示的是该复合物是一个中心对称的二聚体单元 之前那个结果我画的初始结构没有画好对称性所以输出文件即为C1  第二次在 GView 里手动把初始结构调成了 C₂对称,加了 symm=loose 关键词,优化后输出显示为 C₂点群 这里是否还需要加上opt=symmetry再做一次限制型优化呢

6万

帖子

99

威望

6万

eV
积分
125116

管理员

公社社长

11#
发表于 Post on 2026-1-28 15:23:30 | 只看该作者 Only view this author
注意把下文看了
实验测定分子结构的方法以及将实验结构用于量子化学计算需要注意的问题
http://sobereva.com/569
很多分子晶体中存在的构象都需要晶体环境稳定化,简单抠出来优化和实验初始结构没可比性
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

10

帖子

0

威望

63

eV
积分
73

Level 2 能力者

12#
 楼主 Author| 发表于 Post on 2026-1-29 15:15:43 | 只看该作者 Only view this author
sobereva 发表于 2026-1-28 15:23
注意把下文看了
实验测定分子结构的方法以及将实验结构用于量子化学计算需要注意的问题
http://sobereva. ...

好的 谢谢sob老师

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-16 20:16 , Processed in 0.199651 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list