计算化学公社

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

[新手求助] 高斯09加电场优化一直不收敛求解,b3lyp/def2tzvp 是文献有的

[复制链接 Copy URL]

4

帖子

0

威望

31

eV
积分
35

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
%chk=1.chk
#p opt=z-matrix b3lyp/def2tzvp nosymm field=xxyz-1000

Title Card Required

0 1
C              
C                  1    1.38308908
C                  2    1.40412040    1  120.98909973
C                  3    1.40349536    2  118.00076569    1    0.00000000    0
C                  4    1.38405086    3  121.05487153    2    0.00000000    0
C                  5    1.39662425    4  120.45003453    3    0.00042957    0
H                  1    1.08294028    2  119.54866122    3 -179.99943285    0
H                  2    1.08186247    1  119.75881234    6  180.00000000    0
H                  4    1.08194847    3  119.24621901    2 -180.00000000    0
H                  5    1.08257960    4  119.32621117    3 -179.99945641    0
C                  3    4.91775735    2  106.80921190    1 -179.99837937    0
C                 11    1.38124100    3  165.79711575    2 -179.99382599    0
C                 12    1.40534404   11  120.85979268    3    0.00000000    0
C                 13    1.40534197   12  118.28380604   11    0.00000000    0
C                 14    1.38124100   13  120.85640128   12    0.00000000    0
C                 11    1.40534197    3   44.94071448    2 -179.99394090    0
H                 11    1.08183735    3   74.21745602    2    0.00604346    0
H                 12    1.08183935   11  119.98140749    3  180.00000000    0
H                 14    1.08183735   13  119.15817050   12  180.00000000    0
H                 15    1.08183935   14  119.98140749   13  180.00000000    0
C                 13    4.91870632   12  106.69823578   11  179.99966120    0
C                 21    1.38405088   13  165.90667313   12  179.99881800    0
C                 22    1.39662380   21  120.45007242   13    0.00000000    0
C                 23    1.39787406   22  118.98918117   21    0.00000000    0
C                 24    1.38308908   23  120.51600982   22    0.00000000    0
C                 21    1.40349569   13   44.85183689   12  179.99881152    0
H                 21    1.08194810   13   74.39441491   12   -0.00116857    0
H                 22    1.08258005   21  119.32620927   13 -180.00000000    0
H                 24    1.08293987   23  119.93532733   22  180.00000000    0
H                 25    1.08186284   24  119.75881028   23  180.00000000    0
C                  3    1.41875300    2  121.00969652    1  180.00000000    0
C                 31    1.20980500    3  179.99444513    2    0.49249587    0
C                 13    1.41809501   12  120.86175731   11 -180.00000000    0
C                 33    1.20980500   13  179.98984427   12    0.26448280    0
S                  6    1.76925971    5  123.15401783    4 -179.99945325    0
S                 23    1.76925997   22  123.15403342   21 -180.00000000    0
H                 36    1.34217615   23   97.04579107   22   -0.00054292    0
H                 35    1.34217636    6   97.04583340    5   -0.00248852    0

错误提示
RMSDP=2.50D-04 MaxDP=4.14D-02 DE= 6.78D+00 OVMax= 5.89D-03

>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(RB3LYP) =  -9025.83345612     A.U. after  129 cycles
            NFock=128  Conv=0.25D-03     -V/T= 7.8113
KE= 1.325122569705D+03 PE=-1.426698716502D+04 EE= 2.220825291968D+03
Convergence failure -- run terminated.
Error termination via Lnk1e in /home-ys/Soft/Gaussian.D01.Linda/g09/l502.exe at Fri Dec 25 22:28:45 2020.
Job cpu time:       1 days  9 hours 48 minutes 51.0 seconds.
File lengths (MBytes):  RWF=    356 Int=      0 D2E=      0 Chk=     12 Scr=      1

655

帖子

1

威望

5322

eV
积分
5997

Level 6 (一方通行)

2#
发表于 Post on 2020-12-31 10:49:34 | 只看该作者 Only view this author
或许可以试试把电场从0逐步加强,每次读之前的收敛的波函数作为初猜?
现代化学以狄拉克的一句“一切化学问题业已解决”为嚆矢。滥觞于经验主义传统的期望正失去它们的借鉴意义。但面对看似不可达的通往天堂之阶梯,我想循伍德沃德“最好的模型是你底物的对映异构体”的信仰好过过早地振翮。
我们怀揣热忱的灵魂天然被赋予对第一性的追求,不屑于单一坐标的约束,钟情于势能面彼端的芬芳。但

1187

帖子

5

威望

2876

eV
积分
4163

Level 6 (一方通行)

3#
发表于 Post on 2020-12-31 12:01:42 | 只看该作者 Only view this author
首先不要用z-matrix,不然怎么知道电场到底是哪个方向。
然后标题有歧义,这个报错本身是自洽场不收敛导致的,和几何优化没有直接关系。
然后这么强的电场,甚至有必要考虑电子是不是会直接电离出去了。

4

帖子

0

威望

31

eV
积分
35

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2020-12-31 14:35:24 | 只看该作者 Only view this author
snljty 发表于 2020-12-31 12:01
首先不要用z-matrix,不然怎么知道电场到底是哪个方向。
然后标题有歧义,这个报错本身是自洽场不收敛导致 ...

好的,顺便我想问一直高斯如何中加偏置电压,谢谢

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

5#
发表于 Post on 2020-12-31 17:32:12 | 只看该作者 Only view this author
刘自余 发表于 2020-12-31 14:35
好的,顺便我想问一直高斯如何中加偏置电压,谢谢

只能加电场,不可能加电压
北京科音自然科学研究中心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

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

6#
发表于 Post on 2020-12-31 17:35:27 | 只看该作者 Only view this author
强电场下SCF不收敛是司空见惯的事。此时电子结构变得偏离平衡状态,gap减小,甚至有的电子会被极化得非常显著乃至近乎脱离体系。仔细看此文的讨论对你有重要帮助:
一篇文章深入揭示外电场对18碳环的超强调控作用
http://sobereva.com/570http://bbs.keinsci.com/thread-19364-1-1.html
强电场下的计算强烈建议带弥散函数,否则没法描述好被严重极化的电子,这也是为什么上文的计算都刻意带了弥散函数。

当前情况写=z-matrix、nosymm毫无意义
如果你真的对电场的方向敏感,必须用笛卡尔坐标写输入文件,并且用nosymm

搞清楚field关键词的使用,如果你是要加特定方向的均匀电场,根本就是这么写的,你加的是十六极场

记得仔细看此文
解决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

199

帖子

2

威望

3559

eV
积分
3798

Level 5 (御坂)

7#
发表于 Post on 2022-4-8 22:57:16 | 只看该作者 Only view this author
sobereva 发表于 2020-12-31 17:35
强电场下SCF不收敛是司空见惯的事。此时电子结构变得偏离平衡状态,gap减小,甚至有的电子会被极化得非常显 ...

请问社长,加外电场优化后也是需要正常算频率确保没有虚频吗?还是不需要

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

8#
发表于 Post on 2022-4-9 00:20:43 | 只看该作者 Only view this author
量化小菜鸡 发表于 2022-4-8 22:57
请问社长,加外电场优化后也是需要正常算频率确保没有虚频吗?还是不需要

需要
诸如下文的研究里的结构都是没虚频的
一篇文章深入揭示外电场对18碳环的超强调控作用
http://sobereva.com/570http://bbs.keinsci.com/thread-19364-1-1.html

电场影响的是势能面,优化、振动分析照常处理和考虑
北京科音自然科学研究中心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

199

帖子

2

威望

3559

eV
积分
3798

Level 5 (御坂)

9#
发表于 Post on 2022-4-10 09:57:03 | 只看该作者 Only view this author
sobereva 发表于 2022-4-9 00:20
需要
诸如下文的研究里的结构都是没虚频的
一篇文章深入揭示外电场对18碳环的超强调控作用

好的,谢谢社长。

49

帖子

0

威望

994

eV
积分
1043

Level 4 (黑子)

10#
发表于 Post on 2022-11-18 16:45:15 | 只看该作者 Only view this author
sobereva 发表于 2020-12-31 17:32
只能加电场,不可能加电压

老师问下,加电场和电压的区别是什么呢
还有如果要计算锂离子在电池放电状态下与溶剂分子的键能,能通过加电场来模拟?如果不能的话有什么其他方法吗

1万

帖子

0

威望

9888

eV
积分
22133

Level 6 (一方通行)

11#
发表于 Post on 2022-11-18 17:30:30 | 只看该作者 Only view this author
Smes 发表于 2022-11-18 09:45
老师问下,加电场和电压的区别是什么呢
还有如果要计算锂离子在电池放电状态下与溶剂分子的键能,能通过 ...

加电场相当于在分子两侧稍远、但不是特别远的地方放两个极板,一个正极一个负极,但不允许电子在极板和分子之间转移。需要加电场相关的关键词来计算,不需要且不能改变电荷。
加电压相当于在分子旁边放一个工作电极,在离分子很远的地方放一个参考电极,且允许工作电极和分子之间发生电子转移。需要通过改变体系的电荷来计算,不需要且不能加电场相关的关键词。
所以如果只加电场的话,电子云会向一个方向偏,可能带动分子结构也发生一些变化,但电子不会被吸到正极上,负极的电子也不会跑到分子上;如果只加电压的话,分子要么不变,要么得到整数个电子,要么失去整数个电子。
因为锂离子电池的电压只有几伏而已,电极间距离起码也是毫米量级,所以电场只有10^3 V/m量级,对分子的影响可以忽略。所以只有电压的影响,也就是说可以认为,锂离子只要不得到电子,与溶剂分子的键能就是不随电池电压而改变的,而如果得到电子,它就不叫锂离子了。如果实验发现锂电池放电过程中锂离子与溶剂结合能有变化,建议优先考虑锂离子浓度变化导致活度系数变化的可能性。

评分 Rate

参与人数
Participants 1
eV +3 收起 理由
Reason
hebrewsnabla + 3

查看全部评分 View all ratings

Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

49

帖子

0

威望

994

eV
积分
1043

Level 4 (黑子)

12#
发表于 Post on 2022-11-22 12:03:32 | 只看该作者 Only view this author
wzkchem5 发表于 2022-11-18 17:30
加电场相当于在分子两侧稍远、但不是特别远的地方放两个极板,一个正极一个负极,但不允许电子在极板和分 ...

谢谢老师,明白了,但是我又有几了个新的问题
1.锂电池循环过程中电解液溶剂分子在正负极界面发生氧化或还原成气体,这个时候可以理解为有电场影响吗?
如果要计算过渡态的话可以通过设置加电场情况下,溶剂分子的缺电子或得电子状态下(开壳层)与电极材料界面做过渡态优化。还是说这个时候仍然可以忽略电场的影响
2.老师我在计算溶剂氧化还原势的过程中,我发现部分溶剂的在氧化态或还原态结构会发生解离(分子部分结构与结构之间距离明显变长,我没做具体分析不敢说直接断键了),如果是真的断键了,那我是否需要直接优化氧化或还原产物,再计算最后的反应总的吉布斯自由能变最后得到氧化还原势呢?还是说仍然算一个整体产物进行继续的计算。还有个问题就是分解的有部分集团是不稳定状态的明显不可能是最终氧化还原反应的产物。但还是把不稳定的中间态当产物进行计算呢?
以上,感谢老师的指导,回答了我好几个困惑

1万

帖子

0

威望

9888

eV
积分
22133

Level 6 (一方通行)

13#
发表于 Post on 2022-11-22 15:06:48 | 只看该作者 Only view this author
Smes 发表于 2022-11-22 05:03
谢谢老师,明白了,但是我又有几了个新的问题
1.锂电池循环过程中电解液溶剂分子在正负极界面发生氧化或 ...

1. 锂电池里的电场对离子的宏观迁移或许有影响,但是对分子尺度发生的任何事情的影响都可以忽略。这个结论和是否生成气体无关。
2. 需要以断键后的状态为基准做计算。氧化还原势是和反应平衡常数挂钩的,而不是和固定分子结构,只转移电子的过程的平衡常数挂钩的,所以如果反应导致分子结构变化,必然要把这个结构变化考虑进去,哪怕结构变化大到成、断键也是如此。
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

49

帖子

0

威望

994

eV
积分
1043

Level 4 (黑子)

14#
发表于 Post on 2022-11-22 18:30:37 | 只看该作者 Only view this author
wzkchem5 发表于 2022-11-22 15:06
1. 锂电池里的电场对离子的宏观迁移或许有影响,但是对分子尺度发生的任何事情的影响都可以忽略。这个结 ...

老师关于2.如果断键,断键的片段要分别优化,分别计算自由能吗

1万

帖子

0

威望

9888

eV
积分
22133

Level 6 (一方通行)

15#
发表于 Post on 2022-11-22 20:56:18 | 只看该作者 Only view this author
Smes 发表于 2022-11-22 11:30
老师关于2.如果断键,断键的片段要分别优化,分别计算自由能吗

断键以后应该是先生成范德华络合物,然后范德华络合物可以解离,也可以不解离。所以问题归结于范德华络合物是解离的多还是不解离的多,这个用Gibbs自由能判断就行了。不过注意考虑浓度校正,因为解离前、解离后的产物的浓度并不处于标准态。
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

本版积分规则 Credits rule

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

GMT+8, 2026-2-24 14:36 , Processed in 0.700253 second(s), 27 queries , Gzip On.

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