计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: chanyeolhgh
打印 Print 上一主题 Last thread 下一主题 Next thread

[新手求助] 长链脂肪酸优化不收敛

[复制链接 Copy URL]

32

帖子

0

威望

83

eV
积分
115

Level 2 能力者

16#
 楼主 Author| 发表于 Post on 2020-9-14 15:06:06 | 只看该作者 Only view this author
本帖最后由 chanyeolhgh 于 2020-9-14 15:09 编辑
ldatea 发表于 2020-9-14 14:17
你第二个输出文件没正确贴上来

第二个输出文件 12.5.7z (193.26 KB, 下载次数 Times of downloads: 1)
不好意思老师 第一次上传附件 带来不便恳请原谅 您看我这次贴出来的对吗

12.5.7z

193.26 KB, 下载次数 Times of downloads: 11

115

帖子

1

威望

1122

eV
积分
1257

Level 4 (黑子)

17#
发表于 Post on 2020-9-14 23:29:46 | 只看该作者 Only view this author
本帖最后由 ldatea 于 2020-9-15 05:34 编辑

后面这个不收敛其实很容易就知道发生了什么,就是13到14结构变化太大,导致算14的时候读取13收敛的波函数作为初猜质量太差,出现scf在规定圈数不收敛
如果结构改变不大,上一个点的波函数作为下一个点的初猜大部分时候会有利于收敛,你前面没有报错就是因为收敛难度低了,但是13-14这一步,结构改变大,可能和你用gaussian默认的Harris泛函初猜质量差不多,收敛难度就回到了你最开始的128圈也不收敛的那种状态,出现128圈不收敛是很正常的。

用gaussviw打开.log,左下角勾选“read intermediate....”然后就可以看到13-14结构发生了巨大的改变。

你可以尝试提取13的坐标(guassview切到13的几何结构,直接保存为.gjf即可),然后opt=(notrust,maxstep=N)    (N一般用3,4,5)避免每次结构改变过大。
用14坐标继续优化也是可以的,但是我感觉14偏离平衡位置更远。


32

帖子

0

威望

83

eV
积分
115

Level 2 能力者

18#
 楼主 Author| 发表于 Post on 2020-9-15 09:49:18 | 只看该作者 Only view this author
ldatea 发表于 2020-9-14 23:29
后面这个不收敛其实很容易就知道发生了什么,就是13到14结构变化太大,导致算14的时候读取13收敛的波函数作 ...

好的老师 感谢您的耐心指导 我这就去试一试

99

帖子

0

威望

1556

eV
积分
1655

Level 5 (御坂)

19#
发表于 Post on 2020-9-15 10:32:12 | 只看该作者 Only view this author
您好,请问可以解释下为什么要带一个电荷吗?

6万

帖子

99

威望

6万

eV
积分
125156

管理员

公社社长

20#
发表于 Post on 2020-9-16 09:32:14 | 只看该作者 Only view this author
chanyeolhgh 发表于 2020-9-14 13:02
感谢Sob老师 给答疑的老师们带来不便深感抱歉 下次一定注意 但老师如果不写tight的话 那关键词int=ultraf ...

设置明显是错的
你当前这个体系明明是中性单重态,你却当阴离子自由基计算,显然电子结构古怪、SCF难收敛
北京科音自然科学研究中心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

32

帖子

0

威望

83

eV
积分
115

Level 2 能力者

21#
 楼主 Author| 发表于 Post on 2020-9-16 15:04:34 | 只看该作者 Only view this author
sobereva 发表于 2020-9-16 09:32
设置明显是错的
你当前这个体系明明是中性单重态,你却当阴离子自由基计算,显然电子结构古怪、SCF难收 ...

Sob老师 我想要计算脂肪酸的电子亲和能 想要分别计算中性和自由基阴离子的能量 自由基阴离子的结构优化都是在中性结构基础上只更改了电荷和自旋多重度 其他的都没改变 请问老师我应该怎样更改 才能合理的优化自由基阴离子呢

32

帖子

0

威望

83

eV
积分
115

Level 2 能力者

22#
 楼主 Author| 发表于 Post on 2020-9-16 15:07:14 | 只看该作者 Only view this author
naonao5205 发表于 2020-9-15 10:32
您好,请问可以解释下为什么要带一个电荷吗?

老师您好 我是想计算电子亲和能 想要分别计算中性和自由基阴离子能量 所以在优化过程中 添加了一个电荷 请问老师这样优化和计算自由基阴离子是不合理的吗 那我应该怎样更改呢 感谢老师

99

帖子

0

威望

1556

eV
积分
1655

Level 5 (御坂)

23#
发表于 Post on 2020-9-16 15:56:35 | 只看该作者 Only view this author
chanyeolhgh 发表于 2020-9-16 15:07
老师您好 我是想计算电子亲和能 想要分别计算中性和自由基阴离子能量 所以在优化过程中 添加了一个电荷  ...

这个目的的话可以不需要优化。

32

帖子

0

威望

83

eV
积分
115

Level 2 能力者

24#
 楼主 Author| 发表于 Post on 2020-9-16 16:04:07 | 只看该作者 Only view this author
naonao5205 发表于 2020-9-16 15:56
这个目的的话可以不需要优化。

感谢老师 但计算绝热电子亲和能不是在中性和阴离子自由基状态下分别优化后的能量差吗 如果不优化 直接计算单点能做差吗 还想请问老师 我的阴离子自由基体系的设置在中性体系下改变电荷和自旋多重度是合理的吗

6万

帖子

99

威望

6万

eV
积分
125156

管理员

公社社长

25#
发表于 Post on 2020-9-17 13:22:02 | 只看该作者 Only view this author
尝试换个泛函,比如wB97XD优化,或许能好点
单纯对于SCF不收敛这个问题,直接的解决方法就是下文说的,没有可补充的
解决SCF不收敛问题的方法
http://sobereva.com/61

另外也可以先尝试把烃链截短,看看能优化出什么样
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2026-2-24 23:40 , Processed in 0.233666 second(s), 22 queries , Gzip On.

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