计算化学公社

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

[Gaussian/gview] gaussian的ADMP模拟中报错L121

[复制链接 Copy URL]

86

帖子

0

威望

523

eV
积分
609

Level 4 (黑子)

本帖最后由 halouhapily 于 2022-6-23 17:25 编辑

各位老师好
我想要计算的内容为稀土离子碰撞有机分子的反应路径,采用的方法是高斯ADMP模块,首先我为了验证ADMP模块输入的正确与否,测试了单个有机分子的动力学,结果输出很顺利(结果见附件hccamd.gjf  hccamd.log),接着我在有机物分子的基础上添加了稀土离子,计算出现了L121的报错(结果见hccaeumd.gjf hccaeumd.log),我没有找到关于这类错误的解决办法,还请老师们不吝赐教,是不是我输入哪里有问题?

admp.rar

630.18 KB, 下载次数 Times of downloads: 12

6万

帖子

99

威望

6万

eV
积分
125153

管理员

公社社长

2#
发表于 Post on 2022-6-23 18:01:40 | 只看该作者 Only view this author
你的格式不对,看我培训班ppt里的例子


北京科音自然科学研究中心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

86

帖子

0

威望

523

eV
积分
609

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2022-6-23 21:39:03 | 只看该作者 Only view this author
sobereva 发表于 2022-6-23 18:01
你的格式不对,看我培训班ppt里的例子

老师,我按您的格式修改之后,还是报相同的错误
gjf文件:
# wb97xd/genecp ADMP(maxpoint=1000,StepSize=2000) IOp(1/82=300)

Title Card Required

3 1
C                  3.69583200   -0.30527000    0.00017000
C                  2.63115800   -1.19213300    0.00002700
C                  1.33765700   -0.67859700   -0.00012300
C                  1.09447700    0.69993100   -0.00021100
C                  2.19198400    1.57689700   -0.00007200
C                  3.48130300    1.07970800    0.00014600
H                  4.70963500   -0.69270500    0.00028300
H                  2.77680600   -2.26637500    0.00005300
C                 -0.27409600    1.12921100   -0.00023100
H                  2.00866000    2.64735800   -0.00008200
H                  4.32706700    1.75872500    0.00033000
C                 -1.28748400    0.23373300   -0.00004600
C                 -0.99618800   -1.19892200    0.00001400
H                 -0.52034300    2.18796400   -0.00037300
O                  0.31176400   -1.57458300   -0.00012800
O                 -1.83143400   -2.08345600    0.00011200
C                 -2.70713600    0.74663000   -0.00003900
O                 -3.67386600   -0.16392400   -0.00001900
O                 -2.93496100    1.93228000    0.00024800
H                 -3.27887900   -1.06463600    0.00027300
Eu                -7.46479772   -1.91266835    0.13079829

C H O 0
6-31G*
****
Eu
MWB52
****

Eu
MWB52


log文件:
TRJ-TRJ-TRJ-TRJ-TRJ-TRJ-TRJ-TRJ-TRJ-TRJ-TRJ-TRJ-TRJ-TRJ-TRJ-TRJ-TRJ-TRJ
Wanted a floating point number as input.
Found a string as input.
C H O 0                                                                        
    ?
Error termination via Lnk1e in /usr/local/G16-A03/g16/l121.exe at Thu Jun 23 20:00:54 2022.
Job cpu time:       0 days  0 hours  0 minutes  0.8 seconds.
Elapsed time:       0 days  0 hours  0 minutes  0.2 seconds.
File lengths (MBytes):  RWF=      6 Int=      0 D2E=      0 Chk=      1 Scr=      1

6万

帖子

99

威望

6万

eV
积分
125153

管理员

公社社长

4#
发表于 Post on 2022-6-23 22:04:20 | 只看该作者 Only view this author
halouhapily 发表于 2022-6-23 21:39
老师,我按您的格式修改之后,还是报相同的错误
gjf文件:
# wb97xd/genecp ADMP(maxpoint=1000,StepSi ...

根本就没按我的ppt写,空行都注明得那么清楚了
北京科音自然科学研究中心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

86

帖子

0

威望

523

eV
积分
609

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2022-6-24 10:23:30 | 只看该作者 Only view this author
sobereva 发表于 2022-6-23 22:04
根本就没按我的ppt写,空行都注明得那么清楚了

可以了可以了,谢谢老师,
是我的问题,我还以为第二个空行是解释说明呢。
万分感谢老师

86

帖子

0

威望

523

eV
积分
609

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2022-6-24 19:47:58 | 只看该作者 Only view this author
本帖最后由 halouhapily 于 2022-6-24 19:54 编辑
sobereva 发表于 2022-6-23 22:04
根本就没按我的ppt写,空行都注明得那么清楚了

老师我还是之前的体系,ADMP模块计算稀土离子碰撞有机分子,但是结果不收敛是什么原因呢,我尝试将弥散基组去掉也不行,加上FullSCF关键词也不能解决
SCF Done:  E(RwB97XD) =  -719.342687741     A.U. after  129 cycles
            NFock=128  Conv=0.22D-06     -V/T= 2.0372
Convergence failure -- run terminated.
Error termination via Lnk1e in /usr/local/G16-A03/g16/l502.exe at Fri Jun 24 19:41:20 2022.
Job cpu time:       0 days  0 hours 15 minutes 46.0 seconds.
Elapsed time:       0 days  0 hours  4 minutes 19.6 seconds.
File lengths (MBytes):  RWF=     29 Int=      0 D2E=      0 Chk=      2 Scr=      1

admp.rar

7.2 KB, 下载次数 Times of downloads: 2

1万

帖子

0

威望

9886

eV
积分
22130

Level 6 (一方通行)

7#
发表于 Post on 2022-6-24 20:43:51 | 只看该作者 Only view this author
halouhapily 发表于 2022-6-24 12:47
老师我还是之前的体系,ADMP模块计算稀土离子碰撞有机分子,但是结果不收敛是什么原因呢,我尝试将弥散基 ...

你研究的是真空里的反应吗?
如果不是,必须加上溶剂,而且必须用显式/隐式混合溶剂模型。如果该加溶剂的时候不加溶剂,不仅scf难收敛,而且即使收敛了结果也不能用。
此外下次跑的时候输入文件写上#p,打出scf迭代具体过程,这样万一还是不收敛,比较好排查原因。
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
本团队长期招收研究生,有意者可私信联系

86

帖子

0

威望

523

eV
积分
609

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2022-6-24 21:45:36 | 只看该作者 Only view this author
wzkchem5 发表于 2022-6-24 20:43
你研究的是真空里的反应吗?
如果不是,必须加上溶剂,而且必须用显式/隐式混合溶剂模型。如果该加溶剂 ...

非常感谢老师您的帮助,我回去测试以下,我还有以下几个问题希望老师能不吝赐教:
1、目前我想做的工作是模拟实验反应路径,实验中的反应物为硝酸铕六水合物与香豆素配合物,反应溶剂为纯乙醇溶剂。在计算模拟中我为了简化计算模型将硝酸铕六水合物用三价铕离子代替了,这个简化合理吗?
2、关于基组和泛函的选择:我在文献中找到的基组和泛函是计算体系激发态(荧光光谱)使用的,这个基组和泛函的组合使用于动力学模拟吗?
3、如果想要节省机时,您建议我使用哪些基组和泛函呢?
4、对于自旋多重度的选择我拿捏不准,对于含有三价稀土离子的体系的自旋多重度应该怎么设置呢

1万

帖子

0

威望

9886

eV
积分
22130

Level 6 (一方通行)

9#
发表于 Post on 2022-6-25 00:34:00 | 只看该作者 Only view this author
halouhapily 发表于 2022-6-24 14:45
非常感谢老师您的帮助,我回去测试以下,我还有以下几个问题希望老师能不吝赐教:
1、目前我想做的工作 ...

1. 极其不合理,任何金属离子必须带着所有的配体,不管配体是不是溶剂
2. 应该去查这种体系结构优化用的泛函和基组。结构优化和MD对泛函和基组的要求是相似的,但是和激发态的要求有可能差别很大
3. 参见http://sobereva.com/272http://sobereva.com/336
4. 由这个离子的价层电子排布,通过Hund规则自己计算自旋多重度
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
本团队长期招收研究生,有意者可私信联系

86

帖子

0

威望

523

eV
积分
609

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2022-6-25 16:15:02 | 只看该作者 Only view this author
wzkchem5 发表于 2022-6-25 00:34
1. 极其不合理,任何金属离子必须带着所有的配体,不管配体是不是溶剂
2. 应该去查这种体系结构优化用的 ...

老师,我按照您说的加上了溶剂,ADMP模块运行了几步之后报了错误,
Error in input density:  TNEA=************ NAE=    53 TNEB=************ NBE=    53.
这个错误是什么原因导致的呢?
我尝试使用SMD和高斯自带的溶剂模型都出现同样的错误。

admp.zip

48.26 KB, 下载次数 Times of downloads: 5

6万

帖子

99

威望

6万

eV
积分
125153

管理员

公社社长

11#
发表于 Post on 2022-6-25 16:39:06 | 只看该作者 Only view this author
halouhapily 发表于 2022-6-25 16:15
老师,我按照您说的加上了溶剂,ADMP模块运行了几步之后报了错误,
Error in input density:  TNEA=**** ...

必须把和Eu配位的溶剂全都直接画出来才行,用隐式溶剂根本没法表现配位溶剂和Eu的显著的配位作用。越是电子结构诡异的情况SCF往往越难收敛。
另外,这体系用wB97XD没什么好处,建议用PBE0-D3(BJ),参考
使用Gaussian做镧系金属配合物的量子化学计算
http://sobereva.com/581http://bbs.keinsci.com/thread-21274-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

86

帖子

0

威望

523

eV
积分
609

Level 4 (黑子)

12#
 楼主 Author| 发表于 Post on 2022-7-5 20:40:41 | 只看该作者 Only view this author
sobereva 发表于 2022-6-25 16:39
必须把和Eu配位的溶剂全都直接画出来才行,用隐式溶剂根本没法表现配位溶剂和Eu的显著的配位作用。越是电 ...

感谢sob老师,动力学部分已经解决了,
老师,还有问题是
1、通过gview打开的potential energy(见图)含不含零点能呢?
2、为什么我的势能曲线这么不稳定,是不是因为我的时间设置不够长,目前跑了5000fs
3、如何将势能提取出来,我想通过origin作图
4、对于过渡态的搜索,可以直接用opt=ts搜索,然后通过判断是否有唯一虚频判断吗?还有什么过渡态搜索的方法吗?
C:\users\liyan\desktop\1.png

1.png (66.94 KB, 下载次数 Times of downloads: 10)

1.png

6万

帖子

99

威望

6万

eV
积分
125153

管理员

公社社长

13#
发表于 Post on 2022-7-6 07:33:58 | 只看该作者 Only view this author
halouhapily 发表于 2022-7-5 20:40
感谢sob老师,动力学部分已经解决了,
老师,还有问题是
1、通过gview打开的potential energy(见图) ...

1 没零点能的事。零点能来自于量子力学考虑核振动,跟ADMP没关系
2 称不上有多不稳定。小体系结合控温,能量没明显波动才怪
3 点右键,自己找
4 简谈Gaussian里找过渡态的关键词opt=TS和QST2、3
http://sobereva.com/460http://bbs.keinsci.com/thread-12093-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

86

帖子

0

威望

523

eV
积分
609

Level 4 (黑子)

14#
 楼主 Author| 发表于 Post on 2022-8-2 19:51:10 | 只看该作者 Only view this author
本帖最后由 halouhapily 于 2022-8-2 19:53 编辑
wzkchem5 发表于 2022-6-25 00:34
1. 极其不合理,任何金属离子必须带着所有的配体,不管配体是不是溶剂
2. 应该去查这种体系结构优化用的 ...

老师,对于稀土离子周围加配体的问题我还想请教一下,
因为反应物是硝酸铕六水合物,我尝试模拟的时候将全部的配体都加上,结果在跑动力学时,两个反应物弹开了,所以我在稀土离子周围只加了一个硝基和一个水分子,这样简化合理吗?还有显式溶剂的问题:我不确定在反应物周围加几个乙醇分子才算合理,也不清楚应该把乙醇分子加到什么位置?

1016

帖子

0

威望

2126

eV
积分
3142

Level 5 (御坂)

天上地下三山亖海五湖六合八荒九州之泛粤大典編輯委員會編外委員

15#
发表于 Post on 2022-8-3 14:27:44 | 只看该作者 Only view this author
halouhapily 发表于 2022-8-2 19:51
老师,对于稀土离子周围加配体的问题我还想请教一下,
因为反应物是硝酸铕六水合物,我尝试模拟的时候将 ...

王老师放大假去玩啦(画外音:别偷吃我的酸菜和香肠!)
个人观点:不要简化,简化不合理。应当将全部配体加上,配体弹开了可能是配合物结构不合理。如果跑得动,在配合物周围加多几个乙醇分子试试,初始结构可以用团簇模型试试。

但是我觉得这种体系最好用第一性原理软件跑,可以加更多溶剂分子,问题是你的f区元素要找个合适的赝势和赝势基组,我没做过我不知道,等方家指正。
ORCA大法好!CP2K大法好!
嶺南粤音 泛粤典 https://jyutjam.org/

山河自落蒼梧月,風雨猶驅草木兵。
san huɔɞ ki lɔk tʰɔŋ ŋ ᵑgut,fʊŋ ʝi ʝiu kʰui tʰɔu ᵐbʊk pɪŋ

本版积分规则 Credits rule

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

GMT+8, 2026-2-24 00:26 , Processed in 0.181589 second(s), 23 queries , Gzip On.

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