计算化学公社

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

[综合交流] 使用sobMECP程序结合Gaussian程序搜索极小能量交叉点

  [复制链接 Copy URL]

300

帖子

6

威望

2711

eV
积分
3131

Level 5 (御坂)

91#
发表于 Post on 2022-1-24 14:51:00 | 只看该作者 Only view this author
sobereva 发表于 2022-1-20 22:37
不是非得有guess(read),这只不过是让sobMECP运行像普通的Gaussian优化一样,每一步读取上一步收敛的波函 ...

社长您好。可以正常运行了,以下是运行中遇到的问题。
优化到某个点就开始振荡。这个点的S0/S1 gap已经很小了,我就以它为起点,再开一个优化任务,步长调成0.000500。但是优化完全是朝着反方向进行的。
  1. (base) [yzhangnn@login-0 sobMECP]$ grep State ReportFile
  2. Energy of First State:     -913.3140316190
  3. Energy of Second State:    -913.3140222100
  4. Energy of First State:     -913.3141242030
  5. Energy of Second State:    -913.3139641900
  6. Energy of First State:     -913.3142225020
  7. Energy of Second State:    -913.3138938610
  8. Energy of First State:     -913.3142930820
  9. Energy of Second State:    -913.3138421970
  10. Energy of First State:     -913.3143875760
  11. Energy of Second State:    -913.3137692700
  12. Energy of First State:     -913.3144722690
  13. Energy of Second State:    -913.3137019970
  14. Energy of First State:     -913.3144550870
  15. Energy of Second State:    -913.3137119090
  16. Energy of First State:     -913.3145141190
  17. Energy of Second State:    -913.3136749680
  18. Energy of First State:     -913.3145143560
  19. Energy of Second State:    -913.3136808920
  20. Energy of First State:     -913.3145192510
  21. Energy of Second State:    -913.3136830340
  22. Energy of First State:     -913.3145278840
  23. Energy of Second State:    -913.3136832340
复制代码

请问有什么参数可以用来帮助它收敛吗?至少不要让它越优化,gap越大。

5万

帖子

99

威望

5万

eV
积分
112381

管理员

公社社长

92#
 楼主 Author| 发表于 Post on 2022-1-24 18:33:15 | 只看该作者 Only view this author
Freeman 发表于 2022-1-24 14:51
社长您好。可以正常运行了,以下是运行中遇到的问题。
优化到某个点就开始振荡。这个点的S0/S1 gap已经 ...

可能是初始Hessian不好所致
可以从更早一些的结构作为起始结构试试
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

59

帖子

0

威望

962

eV
积分
1021

Level 4 (黑子)

93#
发表于 Post on 2022-3-8 19:26:45 | 只看该作者 Only view this author
sob老师,我刚刚安装测试了sobMECP,我的编译器是ifort,编译之后,用软件自带的例子做测试,用了Pt_coord 和 FeO的例子, prepare.sh 和  runfirst.sh 都可以正常运行,但是运行runMECP.sh的时候,报错信息为@: Expression Syntax.  此外没有其他的报错信息。检查了runMECP的文件,没有编程基础,不知道哪里有语法错误
还请sob老师帮忙看下

runMECP.sh

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

5万

帖子

99

威望

5万

eV
积分
112381

管理员

公社社长

94#
 楼主 Author| 发表于 Post on 2022-3-9 04:52:05 | 只看该作者 Only view this author
lianghuashou 发表于 2022-3-8 19:26
sob老师,我刚刚安装测试了sobMECP,我的编译器是ifort,编译之后,用软件自带的例子做测试,用了Pt_coord  ...

我不知道你怎么运行的
这是csh脚本,当前运行环境必须有csh
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

101

帖子

0

威望

882

eV
积分
983

Level 4 (黑子)

95#
发表于 Post on 2022-3-25 12:22:48 | 只看该作者 Only view this author
老师我之前运行sobMECP时都能正常运行,但我没注意sobMECP在后台运行时又运行了一个,导致之后运行时就出现了如下报错,自己尝试解决未果,希望老师能帮我排查下原因,教授我解决方案 Job0_A.log (62.58 KB, 下载次数 Times of downloads: 1) Job0_B.log (103.4 KB, 下载次数 Times of downloads: 1)

5万

帖子

99

威望

5万

eV
积分
112381

管理员

公社社长

96#
 楼主 Author| 发表于 Post on 2022-3-25 15:43:23 | 只看该作者 Only view this author
好兄弟 发表于 2022-3-25 12:22
老师我之前运行sobMECP时都能正常运行,但我没注意sobMECP在后台运行时又运行了一个,导致之后运行时就出现 ...

重算
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

101

帖子

0

威望

882

eV
积分
983

Level 4 (黑子)

97#
发表于 Post on 2022-3-25 15:57:34 | 只看该作者 Only view this author

我重算了依然不行,,,

5万

帖子

99

威望

5万

eV
积分
112381

管理员

公社社长

98#
 楼主 Author| 发表于 Post on 2022-3-25 16:02:28 | 只看该作者 Only view this author
好兄弟 发表于 2022-3-25 12:22
老师我之前运行sobMECP时都能正常运行,但我没注意sobMECP在后台运行时又运行了一个,导致之后运行时就出现 ...

log文件看不出问题
之前的文件彻底清理干净,先把帖子里的例子跑出来再说
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

101

帖子

0

威望

882

eV
积分
983

Level 4 (黑子)

99#
发表于 Post on 2022-4-10 10:51:03 | 只看该作者 Only view this author
老师您好,在输出五个yes得到mecp下结构后,是否还需要通过查看某些数值等方式去验证是否计算的是正确的mecp下的结构?

5万

帖子

99

威望

5万

eV
积分
112381

管理员

公社社长

100#
 楼主 Author| 发表于 Post on 2022-4-11 00:26:31 | 只看该作者 Only view this author
好兄弟 发表于 2022-4-10 10:51
老师您好,在输出五个yes得到mecp下结构后,是否还需要通过查看某些数值等方式去验证是否计算的是正确的mec ...

达到收敛标准,而且肉眼看上去这结构就是你想要的,就够了。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

10

帖子

0

威望

157

eV
积分
167

Level 3 能力者

101#
发表于 Post on 2023-4-4 16:51:32 | 只看该作者 Only view this author
风影月瞳 发表于 2018-5-2 10:21
谢谢老师,已经搞定,我的高斯版本里TD能量项的关键词是TD-HF/TD-DFT,改一下就正常了。

请问要改成什么呀?我也遇到了同样的问题。

10

帖子

0

威望

157

eV
积分
167

Level 3 能力者

102#
发表于 Post on 2023-4-4 16:58:28 | 只看该作者 Only view this author
sui 发表于 2021-3-16 14:46
请问Sob老师,我运行了测试文件没有问题,算自己的文件时提示
ProgFile Exists - OK
First Input OK

遇到了相同问题,请问您是怎么解决的?

10

帖子

0

威望

157

eV
积分
167

Level 3 能力者

103#
发表于 Post on 2023-4-4 17:12:12 | 只看该作者 Only view this author
好兄弟 发表于 2022-3-25 12:22
老师我之前运行sobMECP时都能正常运行,但我没注意sobMECP在后台运行时又运行了一个,导致之后运行时就出现 ...

请问这个问题是怎么解决的呀?

22

帖子

0

威望

1083

eV
积分
1105

Level 4 (黑子)

104#
发表于 Post on 2023-6-11 22:21:38 | 只看该作者 Only view this author
sob老师,请教一下:
如果你使用了sobMECP发表了文章,请在文中以这种格式引用:Tian Lu, sobMECP program, http://sobereva.com/286 (accessed month day, year)

这里的accessed month day, year指的是那个时间?是“First release: 2015-Mar-18    Last update: 2022-May-27”中的 2022-May-27?

5万

帖子

99

威望

5万

eV
积分
112381

管理员

公社社长

105#
 楼主 Author| 发表于 Post on 2023-6-12 10:42:16 | 只看该作者 Only view this author
yejlyejl 发表于 2023-6-11 22:21
sob老师,请教一下:
如果你使用了sobMECP发表了文章,请在文中以这种格式引用:Tian Lu, sobMECP program ...

最后一次你访问sobMECP页面的时间
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

本版积分规则 Credits rule

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

GMT+8, 2024-11-25 07:42 , Processed in 0.169203 second(s), 22 queries , Gzip On.

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