计算化学公社

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

[辅助/分析程序] 求助优化长程校正泛函w参数的optDFTw在集群上安装并使用的相关问题

[复制链接 Copy URL]

143

帖子

0

威望

1548

eV
积分
1691

Level 5 (御坂)

老师您好,我看http://bbs.keinsci.com/thread-4100-1-3.html 这个网址提供了优化长程校正泛函w参数的简便工具,因为g09是在集群上,我直接下网页上的安装包然后上传到集群上,下一步怎么让他变为可以被调用的呢?我试了chmod +x optDFTw 和optDFTw.f90,这两个变成了绿色,然后我如图编写了pbs文件,但是一提交就停了,也没有生成log文件,只有一个.o文件,如图2所示,请问应该怎么做才能使用这个小程序呢?

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

2#
发表于 Post on 2019-12-25 11:42:57 | 只看该作者 Only view this author
如置顶的新社员必读贴和论坛首页的公告栏所示,求助帖必须在帖子标题明确体现出“提问”、“求助”要素及具体内容(仔细看http://bbs.keinsci.com/thread-9348-1-1.html),我已把你的严重不恰当标题“优化长程校正泛函w参数的简便工具optDFTw”改了,以后务必注意。你的标题明显有非常严重的问题,会让别人以为你的帖子是分享一个新的叫做optDFTw的程序。在本论坛给帖子设置标题时一定要非常谨慎、重视。
北京科音自然科学研究中心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!

143

帖子

0

威望

1548

eV
积分
1691

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2019-12-25 15:00:14 | 只看该作者 Only view this author
sobereva 发表于 2019-12-25 11:42
如置顶的新社员必读贴和论坛首页的公告栏所示,求助帖必须在帖子标题明确体现出“提问”、“求助”要素及具 ...

谢谢老师,我已把标题修改了,请问应该操作才能在集群上使用这个小程序呢

1102

帖子

18

威望

6643

eV
积分
8105

Level 6 (一方通行)

計算化学の社畜

4#
发表于 Post on 2019-12-26 13:38:38 | 只看该作者 Only view this author
小西瓜123 发表于 2019-12-25 15:00
谢谢老师,我已把标题修改了,请问应该操作才能在集群上使用这个小程序呢

编译好的optDFTw 文件丢入g09程序所在目录,并给予其执行权限,然后作业提交脚本可以参考如下脚本进行修改:
  1. #!/bin/sh
  2. #PBS -N Gaussian
  3. #PBS -o job.log
  4. #PBS -e err.log
  5. #PBS -q yuan
  6. #PBS -l nodes=1:ppn=6:yuan
  7. #PBS -l mem=6000mb
  8. #PBS -l walltime=1440:00:00
  9. #PBS -j oe

  10. mkdir -p /tmp/$PBS_JOBID
  11. cp -r $PBS_O_WORKDIR/* /tmp/$PBS_JOBID
  12. cd /tmp/$PBS_JOBID

  13. optDFTw 0.002 0.3 0.151828


  14. cp -r /tmp/$PBS_JOBID/*  $PBS_O_WORKDIR
  15. rm -rf /tmp/$PBS_JOBID
复制代码
Stand on the shoulders of giants

143

帖子

0

威望

1548

eV
积分
1691

Level 5 (御坂)

5#
 楼主 Author| 发表于 Post on 2019-12-29 19:46:40 | 只看该作者 Only view this author
本帖最后由 小西瓜123 于 2019-12-30 14:34 编辑
冰释之川 发表于 2019-12-26 13:38
编译好的optDFTw 文件丢入g09程序所在目录,并给予其执行权限,然后作业提交脚本可以参考如下脚本进行修 ...

好的,谢谢老师,请问应该怎么编译呢?关于Linux系统了解的很少,今天老师给了一个编译后的版本,运行的时候,用默认值能运行,但是会报错, 然后如果是指定的值,有的结果只能显示一步,有的结果直接报错,报错如图 请问这是怎么回事呢?

490

帖子

2

威望

4877

eV
积分
5407

Level 6 (一方通行)

6#
发表于 Post on 2019-12-31 08:54:39 | 只看该作者 Only view this author
小西瓜123 发表于 2019-12-29 19:46
好的,谢谢老师,请问应该怎么编译呢?关于Linux系统了解的很少,今天老师给了一个编译后的版本,运行的 ...

在你使用默认参数的时候前面已经正常计算了N步了,说明optDFTw此时是能够正常运行的,报错的原因也已经提示了“error in finding Alpha virt. eigenvalues”,这不是optDFTw的报错,而是你的高斯任务没有正常结束,你需要打开N.out 或者N+1.out或者N-1.out看看究竟是什么原因导致的报错(极有可能是SCF没有收敛);
后面是指定w的范围,报错原因和上面类似,应该是去找高斯的输出文件。
PS:编译的话,可以使用,ifort optDFTw.f90 -o optDFTw

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
冰释之川 + 5 我很赞同

查看全部评分 View all ratings

143

帖子

0

威望

1548

eV
积分
1691

Level 5 (御坂)

7#
 楼主 Author| 发表于 Post on 2019-12-31 14:15:29 | 只看该作者 Only view this author
本帖最后由 小西瓜123 于 2019-12-31 14:17 编辑
让你变成回忆 发表于 2019-12-31 08:54
在你使用默认参数的时候前面已经正常计算了N步了,说明optDFTw此时是能够正常运行的,报错的原因也已经提 ...

谢谢老师,我刚刚看了N.out,N-1.out是normal的, N+1.out502不收敛了,这种情况应该怎么解决呢?

490

帖子

2

威望

4877

eV
积分
5407

Level 6 (一方通行)

8#
发表于 Post on 2019-12-31 14:42:09 | 只看该作者 Only view this author
小西瓜123 发表于 2019-12-31 14:15
谢谢老师,我刚刚看了N.out,N-1.out是normal的, N+1.out502不收敛了,这种情况应该怎么解决呢?

去论坛里或者http://sobereva.com/网站下找关于如何解决SCF不收敛的帖子。

143

帖子

0

威望

1548

eV
积分
1691

Level 5 (御坂)

9#
 楼主 Author| 发表于 Post on 2019-12-31 15:43:23 | 只看该作者 Only view this author
本帖最后由 小西瓜123 于 2019-12-31 15:46 编辑
让你变成回忆 发表于 2019-12-31 14:42
去论坛里或者http://sobereva.com/网站下找关于如何解决SCF不收敛的帖子。

谢谢老师,请问相应解决SCF不收敛的关键词是加在N+1.gjf ,N-1.gjf, N.gjf这几个文件吗?还是直接加在template.gjf 呢?那我可以分别优化正负离子然后把坐标放到对应的输入文件,然后再用optDFTw程序得到w吗?

490

帖子

2

威望

4877

eV
积分
5407

Level 6 (一方通行)

10#
发表于 Post on 2019-12-31 17:12:27 | 只看该作者 Only view this author
小西瓜123 发表于 2019-12-31 15:43
谢谢老师,请问相应解决SCF不收敛的关键词是加在N+1.gjf ,N-1.gjf, N.gjf这几个文件吗?还是直接加在temp ...

直接在template.gjf加就行了;
你优化正负离子干啥?整个optDFTw过程中,分子的坐标都是没变的,仅仅算了正负离子的单点而已。

143

帖子

0

威望

1548

eV
积分
1691

Level 5 (御坂)

11#
 楼主 Author| 发表于 Post on 2019-12-31 20:23:56 | 只看该作者 Only view this author
让你变成回忆 发表于 2019-12-31 17:12
直接在template.gjf加就行了;
你优化正负离子干啥?整个optDFTw过程中,分子的坐标都是没变的,仅仅算 ...

我现在在template.gjf 里加入了关键词,三个文件都Normal了,在等w的迭代计算结果,之前我们是手动计算优化w的,就是三种状态的分子结构是单独优化的,离子态的坐标应该和中性坐标有微小区别吧,然后从0到0.6分别算单点能,进而得到j^2最小值对应的w,

490

帖子

2

威望

4877

eV
积分
5407

Level 6 (一方通行)

12#
发表于 Post on 2020-1-2 08:32:24 | 只看该作者 Only view this author
小西瓜123 发表于 2019-12-31 20:23
我现在在template.gjf 里加入了关键词,三个文件都Normal了,在等w的迭代计算结果,之前我们是手动计算优 ...

还是那句话,优化w和结构优化没有关系,在某个结构下只对应一个最优的w值。
当然,如果你指的是优化w和结构优化到自洽另当别论。

评分 Rate

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

查看全部评分 View all ratings

34

帖子

0

威望

486

eV
积分
520

Level 4 (黑子)

13#
发表于 Post on 2020-9-23 11:04:31 | 只看该作者 Only view this author
本帖最后由 小潘潘 于 2020-9-23 11:08 编辑

老师,这是我提交作业的脚本文件,但输出的结果中还是没有修改w的上限和下限,麻烦老师看一下到底哪里出了错误

ZLZD`XSD3JF2BZ(5(JDJJ{2.png (12.32 KB, 下载次数 Times of downloads: 47)

ZLZD`XSD3JF2BZ(5(JDJJ{2.png

ILEA@J~`{[U34$Q2$_M]}DJ.png (16.64 KB, 下载次数 Times of downloads: 45)

ILEA@J~`{[U34$Q2$_M]}DJ.png

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

14#
发表于 Post on 2020-9-28 13:18:16 | 只看该作者 Only view this author
小潘潘 发表于 2020-9-23 11:04
老师,这是我提交作业的脚本文件,但输出的结果中还是没有修改w的上限和下限,麻烦老师看一下到底哪里出了 ...

后来optDFTw更新过一次,解决了一个bug。你重新下载最新版本再试试。如果本机能正常计算,PBS提交方式不行,那最好问管理员,PBS的事我不清楚
北京科音自然科学研究中心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!

34

帖子

0

威望

486

eV
积分
520

Level 4 (黑子)

15#
发表于 Post on 2020-9-29 10:10:55 | 只看该作者 Only view this author
sobereva 发表于 2020-9-28 13:18
后来optDFTw更新过一次,解决了一个bug。你重新下载最新版本再试试。如果本机能正常计算,PBS提交方式不 ...

好的,谢谢老师

本版积分规则 Credits rule

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

GMT+8, 2024-11-24 18:28 , Processed in 0.220661 second(s), 25 queries , Gzip On.

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