计算化学公社

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

[辅助/分析程序] 真香:优化w值的optDFTw上传到linux系统可以直接用吗?

[复制链接 Copy URL]

14

帖子

0

威望

49

eV
积分
63

Level 2 能力者

我想要优化WB97XD里的W值,网上下载到sob老师的很厉害的optDFTw。http://sobereva.com/346 解压之后上传到平时做计算的平台

对于帖子中写到的“在template.gjf准备好后直接启动optDFTw就可以进行对w的优化。”不是很理解,怎么启动?用什么命令?
多谢指教!

14

帖子

0

威望

49

eV
积分
63

Level 2 能力者

2#
 楼主 Author| 发表于 Post on 2019-12-4 18:34:22 | 只看该作者 Only view this author
直接输入 optDFTw *.gjf, 得到的结果是:./optDFTw: Permission denied

1187

帖子

5

威望

2876

eV
积分
4163

Level 6 (一方通行)

3#
发表于 Post on 2019-12-4 18:40:08 | 只看该作者 Only view this author
Babily777 发表于 2019-12-4 18:34
直接输入 optDFTw *.gjf, 得到的结果是:./optDFTw: Permission denied

chmod +x optDFTw

14

帖子

0

威望

49

eV
积分
63

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2019-12-4 19:09:47 | 只看该作者 Only view this author

非常感谢,可以运行了!然后出错了。。。
forrtl: severe (59): list-directed I/O syntax error, unit -5, file Internal List-Directed Read
Image              PC                Routine            Line        Source            
optDFTw            000000000048061A  Unknown               Unknown  Unknown
optDFTw            000000000047F116  Unknown               Unknown  Unknown
optDFTw            00000000004417B0  Unknown               Unknown  Unknown
optDFTw            000000000040B75E  Unknown               Unknown  Unknown
optDFTw            000000000040AC9F  Unknown               Unknown  Unknown
optDFTw            000000000042035B  Unknown               Unknown  Unknown
optDFTw            000000000041EB2E  Unknown               Unknown  Unknown
optDFTw            0000000000402D39  Unknown               Unknown  Unknown
optDFTw            0000000000402A9C  Unknown               Unknown  Unknown
libc.so.6          00002B9A4E4F0B35  Unknown               Unknown  Unknown
optDFTw            0000000000402999  Unknown               Unknown  Unknown

1187

帖子

5

威望

2876

eV
积分
4163

Level 6 (一方通行)

5#
发表于 Post on 2019-12-4 19:14:07 | 只看该作者 Only view this author
Babily777 发表于 2019-12-4 19:09
非常感谢,可以运行了!然后出错了。。。
forrtl: severe (59): list-directed I/O syntax error, unit  ...

上传输入文件,格式应该有问题。

14

帖子

0

威望

49

eV
积分
63

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2019-12-4 19:54:21 | 只看该作者 Only view this author
snljty 发表于 2019-12-4 19:14
上传输入文件,格式应该有问题。

每个小时只能发三个帖子,我被限制了。。。
我的输入文件是
%mem=24GB
%nprocshared=12
#  wb97xd/6-31g(d,p)

我查阅了论坛里的其他帖子,说是要把optDFTw放在gaussian安装的目录下,这是必须的吗?如果是,我想请教一下我怎么看我们计算平台上的gaussian安装在哪里呢?非常感谢!

1187

帖子

5

威望

2876

eV
积分
4163

Level 6 (一方通行)

7#
发表于 Post on 2019-12-4 20:00:24 | 只看该作者 Only view this author
Babily777 发表于 2019-12-4 19:54
每个小时只能发三个帖子,我被限制了。。。
我的输入文件是
%mem=24GB

不需要,但是你的g09或者g16得在$PATH里。就是终端输入g09或者g16得找得着命令。$GAUSS_EXEDIR这个路径也要正确设置。

1102

帖子

18

威望

6703

eV
积分
8165

Level 6 (一方通行)

計算化学の社畜

8#
发表于 Post on 2019-12-4 20:44:59 | 只看该作者 Only view this author
本帖最后由 冰释之川 于 2019-12-4 20:47 编辑

在template.gjf准备好后直接启动optDFTw就可以进行对w的优化。Brent优化算法需要给定初始的w范围以及初猜,给得越合适越可能用较少步数收敛,范围一定要能够把实际的w值囊括在此范围中。默认的w下限是0.05(不能写0,否则Gaussian没法运行),默认上限是0.6(一般足够大了),默认的w收敛限是0.0001,默认的初猜值是上下限的中间值。如果要自己设这些参数,需要以命令行方式运行,即:
                optDFTw [下限] [上限] [初猜] [收敛限]
没设的参数会自动用默认值。迭代次数上限是100,如果要改的话需要改源代码里的maxit参数。

——以上内容均摘自《优化长程校正泛函w参数的简便工具optDFTw》http://sobereva.com/346
Stand on the shoulders of giants

1102

帖子

18

威望

6703

eV
积分
8165

Level 6 (一方通行)

計算化学の社畜

9#
发表于 Post on 2019-12-4 20:48:40 | 只看该作者 Only view this author
Babily777 发表于 2019-12-4 19:54
每个小时只能发三个帖子,我被限制了。。。
我的输入文件是
%mem=24GB

把optDFTw放在工作目录下(与template.gjf模板文件处于同一目录即可)
Stand on the shoulders of giants

14

帖子

0

威望

49

eV
积分
63

Level 2 能力者

10#
 楼主 Author| 发表于 Post on 2019-12-4 21:14:43 | 只看该作者 Only view this author
冰释之川 发表于 2019-12-4 20:48
把optDFTw放在工作目录下(与template.gjf模板文件处于同一目录即可)

谢谢!我写了一个sub文件成功把我的任务算上啦!
我想请问一下,我可以在template.gjf 中加入溶剂在有溶剂的情况下优化我的W吗?
比如scrf=(solvent=generic,read)
谢谢!

1102

帖子

18

威望

6703

eV
积分
8165

Level 6 (一方通行)

計算化学の社畜

11#
发表于 Post on 2019-12-4 22:02:42 | 只看该作者 Only view this author
Babily777 发表于 2019-12-4 21:14
谢谢!我写了一个sub文件成功把我的任务算上啦!
我想请问一下,我可以在template.gjf 中加入溶剂在有溶 ...

可以
Stand on the shoulders of giants

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

12#
发表于 Post on 2019-12-5 04:43:19 | 只看该作者 Only view this author
Babily777 发表于 2019-12-4 21:14
谢谢!我写了一个sub文件成功把我的任务算上啦!
我想请问一下,我可以在template.gjf 中加入溶剂在有溶 ...

一般不建议在溶剂模型下优化w,而应当真空下优化w,哪怕最终是要用于溶剂环境中的计算
北京科音自然科学研究中心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

14

帖子

0

威望

49

eV
积分
63

Level 2 能力者

13#
 楼主 Author| 发表于 Post on 2019-12-5 16:05:11 | 只看该作者 Only view this author
sobereva 发表于 2019-12-5 04:43
一般不建议在溶剂模型下优化w,而应当真空下优化w,哪怕最终是要用于溶剂环境中的计算

如您所说,加了溶剂之后读取输入文件时就报错了:
Using the following non-standard input for PCM:
  1 2 1.0 3 2.0 27 1.0
------------------------------------------------------------------------------
The following wrong keyword(s) have been found:






1.0
------------------------------------------------------------------------------
Error termination via Lnk1e in /home/scicons/Gaussian09.D01/l301.exe at Thu Dec  5 14:22:06 2019.

您看还能拯救一下吗?
我用没加溶剂优化出来的W(12.09)计算了后续激发态的能量,结果和文献报道的差距还挺大的。
scrf=(solvent=generic,read)
eps=3

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

14#
发表于 Post on 2019-12-5 17:17:07 | 只看该作者 Only view this author
Babily777 发表于 2019-12-5 16:05
如您所说,加了溶剂之后读取输入文件时就报错了:
Using the following non-standard input for PCM:
...

别写solvent=generic
北京科音自然科学研究中心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

4

帖子

0

威望

155

eV
积分
159

Level 3 能力者

15#
发表于 Post on 2023-9-29 15:50:10 | 只看该作者 Only view this author
本帖最后由 zhangwenyu 于 2023-9-29 16:03 编辑
snljty 发表于 2019-12-4 20:00
不需要,但是你的g09或者g16得在$PATH里。就是终端输入g09或者g16得找得着命令。$GAUSS_EXEDIR这个路径也 ...

老师您好,我也看了关于sob老师的推送http://sobereva.com/346,但是对于其中“Windows下运行之前需要在系统中添加GAUSS_EXEDIR环境变量使之指向g09.exe或g16.exe所在目录,并且在PATH环境变量里也添加这个目录,使得能通过命令行方式顺利调用Gaussian。”这项操作困扰了很久,我也出现了相同的问题,如图。也发现您也说了这个问题,但是我还是不太懂这个怎么设置,麻烦您能不能详细说一下。真的是刚刚入手计算化学,一个人成长太难了

输入文件.png (33.6 KB, 下载次数 Times of downloads: 40)

输入文件.png

报错.png (37.04 KB, 下载次数 Times of downloads: 36)

报错.png

本版积分规则 Credits rule

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

GMT+8, 2026-2-22 06:21 , Processed in 0.560074 second(s), 29 queries , Gzip On.

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