计算化学公社

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

[Gaussian/gview] 求助。金属cluster几何优化不收敛

[复制链接 Copy URL]

19

帖子

0

威望

179

eV
积分
198

Level 3 能力者

本帖最后由 lucl13 于 2021-5-6 09:21 编辑

如题。金属cluster的几何优化在气相中可以顺利完成,以气相优化的结构作为初始结构在scrf=smd下优化就无法收敛。常用的方法都试过了(calcfc甚至calcall, int=ultrafine),感觉没有用,该怎么办。

用的理论方法是PBE0。基组为def2svp。打算在计算单点能时使用def2tzvp

4290

帖子

4

威望

9558

eV
积分
13928

Level 6 (一方通行)

MOKIT开发者

2#
发表于 Post on 2021-5-5 17:32:58 | 只看该作者 Only view this author
建议改用IEFPCM做结构优化,而非SMD。至于算大基组的单点,IEFPCM或SMD随你挑。
自动做多参考态计算的程序MOKIT

19

帖子

0

威望

179

eV
积分
198

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2021-5-5 17:36:47 | 只看该作者 Only view this author
还想请教下,我的理论方法和基组的选用合理吗?是不是优化基组应该再好一点

444

帖子

0

威望

2580

eV
积分
3024

Level 5 (御坂)

娃娃儿鱼

4#
发表于 Post on 2021-5-5 17:41:14 | 只看该作者 Only view this author
本帖最后由 hdhxx123 于 2021-5-5 17:43 编辑
lucl13 发表于 2021-5-5 17:36
还想请教下,我的理论方法和基组的选用合理吗?是不是优化基组应该再好一点

不知道你是啥金属,可以用高斯中pbe1pbe(即pbe0泛函),基组没啥问题,想再高点用defTZVP(高斯关键词是TZVP)或者更好的def2tzvp。泛函看这篇:
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272http://bbs.keinsci.com/thread-536-1-1.html
真·探

19

帖子

0

威望

179

eV
积分
198

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2021-5-5 17:50:14 | 只看该作者 Only view this author
hdhxx123 发表于 2021-5-5 17:41
不知道你是啥金属,可以用高斯中pbe1pbe(即pbe0泛函),基组没啥问题,想再高点用defTZVP(高斯关键词是 ...

我用的就是pbe1pbe。打算单点用def2tzvp,优化的时候用的def2svp 是不是太低了

19

帖子

0

威望

179

eV
积分
198

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2021-5-5 17:54:08 | 只看该作者 Only view this author
hdhxx123 发表于 2021-5-5 17:41
不知道你是啥金属,可以用高斯中pbe1pbe(即pbe0泛函),基组没啥问题,想再高点用defTZVP(高斯关键词是 ...

体系为Au团簇加有机配体

19

帖子

0

威望

179

eV
积分
198

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2021-5-5 17:55:31 | 只看该作者 Only view this author
zjxitcc 发表于 2021-5-5 17:32
建议改用IEFPCM做结构优化,而非SMD。至于算大基组的单点,IEFPCM或SMD随你挑。

好的。谢谢

6万

帖子

99

威望

6万

eV
积分
125204

管理员

公社社长

8#
发表于 Post on 2021-5-6 02:26:50 | 只看该作者 Only view this author
lucl13 发表于 2021-5-5 17:50
我用的就是pbe1pbe。打算单点用def2tzvp,优化的时候用的def2svp 是不是太低了

发帖时候注意严谨。明明1L你写的是PBE

级别问题不大。不说多少个原子,没法在基组上进行准确的推荐

这种体系优化没必要带溶剂模型
北京科音自然科学研究中心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

19

帖子

0

威望

179

eV
积分
198

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2021-5-6 09:17:33 | 只看该作者 Only view this author
本帖最后由 lucl13 于 2021-5-6 13:11 编辑
sobereva 发表于 2021-5-6 02:26
发帖时候注意严谨。明明1L你写的是PBE

级别问题不大。不说多少个原子,没法在基组上进行准确的推荐

下次一定注意。谢谢sob老师。

我的体系是蛋白质和Au配位形成团簇(QM计算中抠出了结合位点区域,包括几个氨基酸残基和Au团簇),在实验中是以溶液的形式存在的。在这种情况下也不需要溶剂模型吗?还是说算单点的时候带溶剂模型就可以了?

关于体系大小的话。总原子数目为100个原子,其中Au原子数目为10个左右。希望sob老师给推荐一下 ☺️

1665

帖子

5

威望

4788

eV
积分
6553

Level 6 (一方通行)

喵星人

10#
发表于 Post on 2021-5-6 11:29:30 | 只看该作者 Only view this author
蛋白质和Au配位你是怎么处理的。。。

19

帖子

0

威望

179

eV
积分
198

Level 3 能力者

11#
 楼主 Author| 发表于 Post on 2021-5-6 15:00:43 | 只看该作者 Only view this author
喵星大佬 发表于 2021-5-6 11:29
蛋白质和Au配位你是怎么处理的。。。

初始结构从晶体结构中得到的

1665

帖子

5

威望

4788

eV
积分
6553

Level 6 (一方通行)

喵星人

12#
发表于 Post on 2021-5-6 15:21:05 | 只看该作者 Only view this author
lucl13 发表于 2021-5-6 15:00
初始结构从晶体结构中得到的

我的意思是蛋白质你总不能全都放进去优化吧。。。。是扣了个团簇还是怎么处理的

19

帖子

0

威望

179

eV
积分
198

Level 3 能力者

13#
 楼主 Author| 发表于 Post on 2021-5-6 19:31:44 | 只看该作者 Only view this author
喵星大佬 发表于 2021-5-6 15:21
我的意思是蛋白质你总不能全都放进去优化吧。。。。是扣了个团簇还是怎么处理的

对的 就是抠出了团簇和结合氨基酸残基

6万

帖子

99

威望

6万

eV
积分
125204

管理员

公社社长

14#
发表于 Post on 2021-5-7 03:00:31 | 只看该作者 Only view this author
lucl13 发表于 2021-5-6 09:17
下次一定注意。谢谢sob老师。

我的体系是蛋白质和Au配位形成团簇(QM计算中抠出了结合位点区域,包括 ...

PBE0-D3(BJ)。优化时候Au用SDD,其它用6-311G*。单点时候可以都用def2-TZVP。次要区域可以把基组降低一档以节约时间。
北京科音自然科学研究中心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

19

帖子

0

威望

179

eV
积分
198

Level 3 能力者

15#
 楼主 Author| 发表于 Post on 2021-5-7 15:47:03 | 只看该作者 Only view this author
sobereva 发表于 2021-5-7 03:00
PBE0-D3(BJ)。优化时候Au用SDD,其它用6-311G*。单点时候可以都用def2-TZVP。次要区域可以把基组降低一档 ...

谢谢sob老师

本版积分规则 Credits rule

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

GMT+8, 2026-2-27 01:34 , Processed in 0.168400 second(s), 20 queries , Gzip On.

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