计算化学公社

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

[Gaussian/gview] 在G16中使用MN15-D3(BJ)、MN15L-D3(0)的方法

[复制链接 Copy URL]

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

跳转到指定楼层 Go to specific reply
楼主
从G16开始支持比较新的明尼苏达系列MN15和MN15L,但G16目前版本(B.01)还尚未纳入相应的D3参数,正好今天有人问怎么对MN15计算加D3校正,我遂对之前写的《DFT-D色散校正的使用》http://sobereva.com/210 进行了扩展,其中4.4节介绍了如何实现。
北京科音自然科学研究中心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

131

帖子

0

威望

927

eV
积分
1058

Level 4 (黑子)

2#
发表于 Post on 2020-6-26 15:32:20 | 只看该作者 Only view this author
sob老师好!
我在做oniom计算时高层为dimer,想在高层使用MN15-D3(BJ)泛函。计算的关键词设置如下:
#p opt freq oniom(mn15/6-31g(d) em=gd3bj:uff=qeq)=embedcharge geom=connectivity
并且在~/.bashrc中自定义了如下参数:
export GAUSS_DFTD3_S6=1000000
export GAUSS_DFTD3_S8=786200
export GAUSS_DFTD3_ABJ1=2097100
export GAUSS_DFTD3_ABJ2=7592300
但是计算时还是报错:
Two-electron integral symmetry is turned on.
R6DS8: Unable to choose the S8 parameter, IExCor= 7252 IXCFnc= 73 ScaHFX=  0.440000 IDFTD=4
请问这是什么原因呢?谢谢老师。

164

帖子

0

威望

1991

eV
积分
2155

Level 5 (御坂)

3#
发表于 Post on 2020-6-27 00:43:28 | 只看该作者 Only view this author
暖空 发表于 2020-6-26 15:32
sob老师好!
我在做oniom计算时高层为dimer,想在高层使用MN15-D3(BJ)泛函。计算的关键词设置如下:
#p o ...

source ~/.bashrc了吗

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

4#
 楼主 Author| 发表于 Post on 2020-6-27 01:05:35 | 只看该作者 Only view this author
暖空 发表于 2020-6-26 15:32
sob老师好!
我在做oniom计算时高层为dimer,想在高层使用MN15-D3(BJ)泛函。计算的关键词设置如下:
#p o ...

运行一下bash命令,或者重新进入终端,才能令.bashrc里新加入的设置生效
北京科音自然科学研究中心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

131

帖子

0

威望

927

eV
积分
1058

Level 4 (黑子)

5#
发表于 Post on 2020-6-27 10:55:55 | 只看该作者 Only view this author
thanhtam 发表于 2020-6-27 00:43
source ~/.bashrc了吗

source了哈

131

帖子

0

威望

927

eV
积分
1058

Level 4 (黑子)

6#
发表于 Post on 2020-6-27 11:01:43 | 只看该作者 Only view this author
sobereva 发表于 2020-6-27 01:05
运行一下bash命令,或者重新进入终端,才能令.bashrc里新加入的设置生效

sob老师好~
我添加完环境变量以后source了,也重进了终端,并且echo $GAUSS_DFTD3_S6的时候有显示数值1000000,但是提交的时候还是相同的报错:
Two-electron integral symmetry is turned on.
R6DS8: Unable to choose the S8 parameter, IExCor= 7252 IXCFnc= 73 ScaHFX=  0.440000 IDFTD=4
Error termination via Lnk1e in /share/apps/g16-avx2/l301.exe at Sat Jun 27 10:54:22 2020.
请问还可能是哪里的设置出问题了呢?谢谢老师。

164

帖子

0

威望

1991

eV
积分
2155

Level 5 (御坂)

7#
发表于 Post on 2020-6-27 12:21:45 | 只看该作者 Only view this author
暖空 发表于 2020-6-27 11:01
sob老师好~
我添加完环境变量以后source了,也重进了终端,并且echo $GAUSS_DFTD3_S6的时候有显示数值10 ...

报错的是S8,你echo S6干啥

131

帖子

0

威望

927

eV
积分
1058

Level 4 (黑子)

8#
发表于 Post on 2020-6-27 15:58:56 | 只看该作者 Only view this author
thanhtam 发表于 2020-6-27 12:21
报错的是S8,你echo S6干啥

这只是测试环境变量有没有生效的一种方法,sob老师之前的帖子中有说过

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

9#
 楼主 Author| 发表于 Post on 2020-6-27 18:09:34 | 只看该作者 Only view this author
远程说不清楚
我不知道你要算什么,MN15能派上用场的地方很少,还不如直接改成别的泛函,比如M06-2X、wB97XD,或者自带了D3参数的泛函
北京科音自然科学研究中心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

131

帖子

0

威望

927

eV
积分
1058

Level 4 (黑子)

10#
发表于 Post on 2020-6-28 08:56:10 | 只看该作者 Only view this author
sobereva 发表于 2020-6-27 18:09
远程说不清楚
我不知道你要算什么,MN15能派上用场的地方很少,还不如直接改成别的泛函,比如M06-2X、wB97 ...

感谢回复。
我是做发光的,常见的几种不同HF成分的泛函我都测试了,发现只有mn15的数据和实验值最符合(吸收和发射),所以就选择了mn15。现在想测试一下高层是dimer时发光性能有什么变化,所以需要在mn15的基础上加上色散校正。

131

帖子

0

威望

927

eV
积分
1058

Level 4 (黑子)

11#
发表于 Post on 2020-6-28 09:25:13 | 只看该作者 Only view this author
sobereva 发表于 2020-6-27 18:09
远程说不清楚
我不知道你要算什么,MN15能派上用场的地方很少,还不如直接改成别的泛函,比如M06-2X、wB97 ...

sob老师好。
我还有一个问题,就是已经通过环境变量定义了MN15的DFT-D3参数,那我用其他加色散校正的泛函计算时是直接用其内置的参数还是我自定义的DFT-D3参数呢?

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

12#
 楼主 Author| 发表于 Post on 2020-6-29 01:03:16 | 只看该作者 Only view this author
暖空 发表于 2020-6-28 09:25
sob老师好。
我还有一个问题,就是已经通过环境变量定义了MN15的DFT-D3参数,那我用其他加色散校正的泛 ...

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

131

帖子

0

威望

927

eV
积分
1058

Level 4 (黑子)

13#
发表于 Post on 2020-7-1 11:17:25 | 只看该作者 Only view this author

sob老师好。
我给高斯的技术支持发邮件咨询了,是因为G16中改为了用IOP定义色散,我测试了一个MN15-D3(BJ)的计算任务,可以正常结束,以下是高斯技术的回复邮件:
欢迎来信咨询。环境变量定义色散的方法是用于G09的,而在G16中变为支持直接使用IOp定义色散,IOp(3/174=1000000,3/175=786200,3/177=2097100,3/178=7592300) 。
注意IOp是内部参数调试选项,因此多步计算中IOP只能适用于第一步,即例如"opt freq"中IOp只能用于Opt一步,对Freq步骤无效(见gaussian.com/iop/)。因此需要分成两个任务,各自使用同样的IOp。

希望可以帮助到有需要的小伙伴们。

评分 Rate

参与人数
Participants 2
eV +5 收起 理由
Reason
xiaobin + 2 谢谢
thanhtam + 3

查看全部评分 View all ratings

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

14#
 楼主 Author| 发表于 Post on 2020-7-1 20:23:24 | 只看该作者 Only view this author
G16 <=B.01是支持环境变量方式设置D3参数的,<=B.01版手册里搜GAUSS_DFTD3_S8都能搜到,客服的说法有误导性。应当是现在的C.01才把这个给撤了而改为必须用IOp(我没测试)。

已在《DFT-D色散校正的使用》(http://sobereva.com/210)中增加以下补充
-----
2020-Jul-1补充:下文说的通过环境变量自定义DFT-D3参数的做法只适合Gaussian 16 B.01及之前,从G16 C.01开始需要通过IOp直接定义,具体来说,把以下IOp设为NNNNNNNN代表把相应参数设为NNNNNNNN/1000000
IOp(3/174):S6
IOp(3/175):S8
IOp(3/176):SR6
IOp(3/177):A1
IOp(3/178):A2
-----
北京科音自然科学研究中心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

1187

帖子

5

威望

2876

eV
积分
4163

Level 6 (一方通行)

15#
发表于 Post on 2020-12-11 23:32:32 | 只看该作者 Only view this author
给不爱看手册的小伙伴们一个提醒,用IOp的时候,如果某个参数是0,选项应该写成=-1而不是=0。
顺便,我的A.03 Linux avx2版也只能用IOp,环境变量读不到

评分 Rate

参与人数
Participants 7
eV +22 收起 理由
Reason
Xx_DerEnTh + 5 谢谢
wangzuwei + 2 谢谢
xiaobin + 3 谢谢
wang7344412 + 5 谢谢
jitou11 + 5 同样的情况
mt13 + 1 谢谢
达斯维达尼亚 + 1 谢谢

查看全部评分 View all ratings

本版积分规则 Credits rule

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

GMT+8, 2026-2-23 01:33 , Processed in 0.180815 second(s), 21 queries , Gzip On.

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