计算化学公社

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

[综合交流] 简单量子化学问题答疑专帖

  [复制链接 Copy URL]

3

帖子

0

威望

68

eV
积分
71

Level 2 能力者

3961#
发表于 Post on 2024-11-1 01:18:01 | 只看该作者 Only view this author
各位老师好,我使用Gaussian16对一个过渡金属有机螯合物进行结构优化(使用TPSSh泛函和Def2TZVP基组,64核CPU,128GB内存),在第22步时达到几何优化收敛,并输出了“Normal termination of Gaussian”,但程序没能正常自动结束,进行执行了下一步,并卡在了这一步的cycle 1上不再进行,且无明显报错。保留该任务一个多小时后仍没有变化,在此期间CPU始终被“l1110.exe”所占用。附件为任务的输入和输出文件。

想向各位老师请教:
1. 在这种情况下手动终止(Ctrl + C)该任务是否合理?如果不合理,应当如何操作?
2. 依照各位老师的经验,这种情况可能是由于什么原因造成的?如果我想避免这种情况再次发生,有哪些方面可以进一步排查?

由于发现自己没有权限在本版块发帖,因此在答疑帖中发表回复,恳请各位老师赐教。

Gaussian_file.zip

363.47 KB, 下载次数 Times of downloads: 1

5万

帖子

99

威望

5万

eV
积分
112349

管理员

公社社长

3962#
 楼主 Author| 发表于 Post on 2024-11-1 01:24:56 | 只看该作者 Only view this author
Peng14 发表于 2024-11-1 01:18
各位老师好,我使用Gaussian16对一个过渡金属有机螯合物进行结构优化(使用TPSSh泛函和Def2TZVP基组,64核C ...

正在做freq任务,对大体系本来耗时就高,极为正常
北京科音自然科学研究中心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!

14

帖子

0

威望

124

eV
积分
138

Level 2 能力者

3963#
发表于 Post on 2024-11-6 16:39:11 | 只看该作者 Only view this author
本帖最后由 Solitude198 于 2024-11-7 10:58 编辑

没事了,过渡态找错了.jpg



energy.py

2.88 KB, 下载次数 Times of downloads: 0

338

帖子

0

威望

1105

eV
积分
1443

Level 4 (黑子)

3964#
发表于 Post on 2024-11-7 17:05:26 | 只看该作者 Only view this author
老师们好,我在升级Windows版ORCA时出现了如下问题:

之前在使用5.0.3版本,今天打算升级到6.0.1,在503存在的情况下,我安装了601(另一个路径),在运行601时出现了并行报错“Calling Command: mpiexec -np 6  C:\ORCA_6.0.1\orca_startup_mpi.exe D-gCP.int.tmp D-gCP”,但如果使用503版本来运行则可以顺利结束。后续我把503版本删除后(包括环境变量),依然存在并行报错的情况。

请问在这种升级ORCA版本的情境下,即使删除了旧版本的ORCA,也要对msmpi做一些处理才可以使得新版本ORCA正常运行吗?

1万

帖子

0

威望

7390

eV
积分
18143

Level 6 (一方通行)

3965#
发表于 Post on 2024-11-7 20:07:31 | 只看该作者 Only view this author
cokie 发表于 2024-11-7 10:05
老师们好,我在升级Windows版ORCA时出现了如下问题:

之前在使用5.0.3版本,今天打算升级到6.0.1,在503 ...

检查MPI的版本是否满足ORCA 6.0.1的要求
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

338

帖子

0

威望

1105

eV
积分
1443

Level 4 (黑子)

3966#
发表于 Post on 2024-11-7 20:25:10 | 只看该作者 Only view this author
wzkchem5 发表于 2024-11-7 20:07
检查MPI的版本是否满足ORCA 6.0.1的要求

下载的是最新的版本,10.1.3 (10.1.12498.52)
在问完问题之后,我还尝试了如下手段:
1. 卸载msmpi后,再次安装ORCA_6.0.1,之后再安装msmpi最新版本(即先安装ORCA,再安装msmpi)
2. 卸载ORCA6.0.1后重新安装(即先安装msmpi,再安装ORCA_6.0.1)
均在并行时报错,我的处理器是7950X,调用6核应该不存调用核心超过物理核心数的问题。

并且我在上述操作结束后,又重新解压了ORCA_5.0.3版本尝试并行,依旧没有问题,只有6.0.1版本并行报错。请问我还可以做哪些尝试呢?

(还有个支线问题,我在安装了ORCA_6.0.1之后,在用户变量中自动生成一个名为ORCADIR的变量,路径是C:\ORCA_6.0.1,请问老师了解这个变量的作用是什么么?)

1万

帖子

0

威望

7390

eV
积分
18143

Level 6 (一方通行)

3967#
发表于 Post on 2024-11-7 21:16:47 | 只看该作者 Only view this author
cokie 发表于 2024-11-7 13:25
下载的是最新的版本,10.1.3 (10.1.12498.52)
在问完问题之后,我还尝试了如下手段:
1. 卸载msmpi后, ...

如果有条件的话可以试一下openmpi,和ORCA的兼容性比msmpi好一些。如果还不行,可以在ORCA论坛的Installation版面问一下
ORCADIR这个变量我不知道是否有用,但是即使没用,留着应该也不碍事
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

338

帖子

0

威望

1105

eV
积分
1443

Level 4 (黑子)

3968#
发表于 Post on 2024-11-7 22:49:28 | 只看该作者 Only view this author
wzkchem5 发表于 2024-11-7 21:16
如果有条件的话可以试一下openmpi,和ORCA的兼容性比msmpi好一些。如果还不行,可以在ORCA论坛的Installa ...

感谢老师的建议,我去了ORCA论坛installation版块进行查找。对6.0.0及之后的windows版本,并行似乎是一个通病(已经有不少人反馈),目前大致知道的方向是:我的安装没有问题,会如此报错是因为任务判断我所给的内存不够。

对于ORCA_6_0_1
我之前所用的关键字是wb97m-v def2-tzvp def2/j def2-tzvp/c,370原子,6核/9000MB,并行报错。
之后我将关键词改为   B3LYP  Def2-svp def2/j  370原子,8核/7500MB,程序顺利进行。

对于ORCA_5.0.3
前些日子做过计算,所用的关键字是wb97m-v def2-tzvp def2/j def2-tzvp/c,370原子,6核/9000MB,任务成功运行并结束。

总结来说,对于win版6.0.0以上的ORCA版本,似乎变成了一个更可怕的内存怪兽:
1. 同级别任务对内存的需求比5系版本要求更大;
2. 如果分配内存不够,不同于原来直到某一步调用内存失败才报错,现在似乎多出一步预估步骤,如果预估认为内存分配不足,则直接从一开始就shut down了。

ORCA什么时候可以给内存不太大的机器一条活路

160

帖子

0

威望

619

eV
积分
779

Level 4 (黑子)

迷茫的本科生一枚

3969#
发表于 Post on 2024-11-7 23:19:36 | 只看该作者 Only view this author
请问老师,柔性扫描如何同时扫描键长和键角,即二者同步变化,我只知道同时扫描两个键长的操作,但是网上搜了一圈没得到如何同时扫描键长和键角的方法,
放下

14

帖子

0

威望

124

eV
积分
138

Level 2 能力者

3970#
发表于 Post on 2024-11-8 10:51:56 | 只看该作者 Only view this author
2877321934 发表于 2024-11-7 23:19
请问老师,柔性扫描如何同时扫描键长和键角,即二者同步变化,我只知道同时扫描两个键长的操作,但是网上搜 ...

http://sobereva.com/474

5 用Gaussian做柔性扫描的方法和实例
5.1 基本用法
……柔性扫描也是可以同时扫多个,只要写多行这种设定即可,对于二维柔性扫描也是可以用GaussView显示出曲面图。……

160

帖子

0

威望

619

eV
积分
779

Level 4 (黑子)

迷茫的本科生一枚

3971#
发表于 Post on 2024-11-8 10:56:50 | 只看该作者 Only view this author
Solitude198 发表于 2024-11-8 10:51
http://sobereva.com/474

5 用Gaussian做柔性扫描的方法和实例

谢谢回答 但我的意思是两者同步变化 而不是绘制出二维势能面 比如说
RCH(NSteps=20,StepSize=-0.10)=R(34,35)
ROR=R(36,20)
F(Frozen)=2.0*RCH-ROR
上面的写法就是同步变化CH和OR键,最终的扫描曲线还是一维,但是同步变化键长和键角我不知道怎么写
放下

14

帖子

0

威望

124

eV
积分
138

Level 2 能力者

3972#
发表于 Post on 2024-11-8 11:27:04 | 只看该作者 Only view this author
本帖最后由 Solitude198 于 2024-11-8 11:42 编辑
2877321934 发表于 2024-11-8 10:56
谢谢回答 但我的意思是两者同步变化 而不是绘制出二维势能面 比如说
RCH(NSteps=20,StepSize=-0.10)=R(3 ...

那要这样利用GIC https://gaussian.com/gic/
对于你的键角

scan(StepSize=-10.0,NSteps=13)=D(6,4,1)

同时开始键角和结束键角必然是已知的,比如 180 和 50
那么假设你的键长 R(2,5) 想要同时从 3.0 到 1.7 就可以固定

rcons(freeze)=R(2,5)-D(6,4,1)/100

这样他们的差会被固定在1.2
当然如果这两个没法相互计算的话就没办法了

1万

帖子

0

威望

7390

eV
积分
18143

Level 6 (一方通行)

3973#
发表于 Post on 2024-11-8 14:25:40 | 只看该作者 Only view this author
cokie 发表于 2024-11-7 15:49
感谢老师的建议,我去了ORCA论坛installation版块进行查找。对6.0.0及之后的windows版本,并行似乎是一个 ...

可以考虑降核数、增加maxcore。ORCA 6相比5有加速,即使因为内存问题不得不降核数,也未必会比ORCA 5慢
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

105

帖子

0

威望

498

eV
积分
603

Level 4 (黑子)

3974#
发表于 Post on 2024-11-8 16:48:08 | 只看该作者 Only view this author

相互作用能计算

本帖最后由 davi 于 2024-11-8 08:01 编辑

像图片这种相互作用能(结合能)是如何计算的呢,翻看帖子没找到合适的,有没有什么思路呢

屏幕截图 2024-11-08 165007.jpg (94.57 KB, 下载次数 Times of downloads: 2)

屏幕截图 2024-11-08 165007.jpg

160

帖子

0

威望

619

eV
积分
779

Level 4 (黑子)

迷茫的本科生一枚

3975#
发表于 Post on 2024-11-8 23:13:24 | 只看该作者 Only view this author
Solitude198 发表于 2024-11-8 11:27
那要这样利用GIC https://gaussian.com/gic/
对于你的键角

ok 谢谢解惑
放下

本版积分规则 Credits rule

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

GMT+8, 2024-11-23 05:08 , Processed in 0.200451 second(s), 24 queries , Gzip On.

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