计算化学公社

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

[综合交流] wb97xd中,关于w的一些问题?

[复制链接 Copy URL]

97

帖子

0

威望

337

eV
积分
434

Level 3 能力者

老师,你好

我现在在尝试使用wb97xd来优化激发态分子,在gview里面,当选择泛函wb97xd后,会有fitting set,老板告诉我这个一般不用设置,我也就没有动过。但是后来发现在关于wb97xd的paper (Jeng-Da Chai and Martin Head-Gordon,2008)中有提到w,但不是非常理解w的具体含义是什么?

另外,文章中有关于w的优化设置,我不知道需不需要设置w的值,或者fitting set来优化计算?如果需要,应该怎么设置呢?

3097

帖子

29

威望

1万

eV
积分
17221

Level 6 (一方通行)

2#
发表于 Post on 2015-11-3 12:41:37 | 只看该作者 Only view this author
本帖最后由 liyuanhe211 于 2015-11-3 16:00 编辑





239

帖子

0

威望

2356

eV
积分
2595

Level 5 (御坂)

3#
发表于 Post on 2015-11-3 14:17:27 | 只看该作者 Only view this author
本帖最后由 brothers 于 2015-11-3 14:28 编辑

没用过,但是看过方程,估计你看看方程也能猜出来。
下面是胡言乱语,因为不是做这个的,说错了还请见谅~

w是定义近程与远程范围的变量,w越大,远程所覆盖的范围也就越大,也就是采用高HF成分的范围越宽。
fitting set,估计是指的原作者在泛函优化过程中所选取的实验数据吧?不同set所获得的最终泛函的参数是有些不一样的。不可能需要自己重新优化一遍的吧~~

楼上给出的介绍是关于D(dispersion)部分的参数,跟长程校正的w是无关的。至于S6,SR6参数的含义,我没看文献,也不是特别清楚。但是damping function,应该就是对色散作用(-6次方)和排斥作用(-12次方)的参数,看看L-J势的数学形式即能作此推测。

坐等专家的详细解读~


评分 Rate

参与人数
Participants 3
eV +13 收起 理由
Reason
qianbao + 5 赞!现在自己也想学习使用这个泛函计算,非.
sobereva + 3
卡开发发 + 5 我很赞同

查看全部评分 View all ratings

十室之邑,必有忠信如丘焉,不如丘之好学也。

6万

帖子

99

威望

6万

eV
积分
125152

管理员

公社社长

4#
发表于 Post on 2015-11-3 17:10:16 | 只看该作者 Only view this author
fitting set和w完全没关系。

fitting set是指用什么密度拟合基组,用密度拟合基组可以省一些时间,但会带来精度的牺牲。高斯中密度拟合基组只能对纯泛函用,带来的加速效果不明显。虽然你用杂化泛函时也可以设密度拟合基组,但是实际计算中是忽略的,设了也没用。所以你不用管fitting set的设置。

w就是3L所解释的。
北京科音自然科学研究中心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

239

帖子

0

威望

2356

eV
积分
2595

Level 5 (御坂)

5#
发表于 Post on 2015-11-3 18:44:00 | 只看该作者 Only view this author
sobereva 发表于 2015-11-3 17:10
fitting set和w完全没关系。

fitting set是指用什么密度拟合基组,用密度拟合基组可以省一些时间,但会 ...

学习了,抽空看看这方面的介绍和文献~
十室之邑,必有忠信如丘焉,不如丘之好学也。

97

帖子

0

威望

337

eV
积分
434

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2015-11-4 14:10:52 | 只看该作者 Only view this author
brothers 发表于 2015-11-3 14:17
没用过,但是看过方程,估计你看看方程也能猜出来。
下面是胡言乱语,因为不是做这个的,说错了还请见谅~
...

谢谢,paper中也提到w是控制interelectronic distance的,以此来区分远近程中HF的成分比例

97

帖子

0

威望

337

eV
积分
434

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2015-11-4 14:16:03 | 只看该作者 Only view this author
sobereva 发表于 2015-11-3 17:10
fitting set和w完全没关系。

fitting set是指用什么密度拟合基组,用密度拟合基组可以省一些时间,但会 ...

谢谢,我看原paper中提到,w的优化是self-consistent的过程,但是我在outputfile中没有找到w参数converge的过程,不知道在哪里能看到?或者w的初始值需要自己用keyword设置什么吗?

6万

帖子

99

威望

6万

eV
积分
125152

管理员

公社社长

8#
发表于 Post on 2015-11-4 16:08:33 | 只看该作者 Only view this author
ccccn345 发表于 2015-11-4 14:16
谢谢,我看原paper中提到,w的优化是self-consistent的过程,但是我在outputfile中没有找到w参数converge ...

自己计算的时候不用管w,这是提出泛函过程中开发者事先拟合好的,找到的最理想的一个值。你直接用wB97XD就完了。
北京科音自然科学研究中心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

97

帖子

0

威望

337

eV
积分
434

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2015-11-5 12:58:41 | 只看该作者 Only view this author
sobereva 发表于 2015-11-4 16:08
自己计算的时候不用管w,这是提出泛函过程中开发者事先拟合好的,找到的最理想的一个值。你直接用wB97XD ...

谢谢

167

帖子

0

威望

1051

eV
积分
1218

Level 4 (黑子)

10#
发表于 Post on 2018-5-30 15:08:08 | 只看该作者 Only view this author
楼主您好,我想请教一下,想在高斯用wB97XD泛函优化,是可以在命令行直接输入就ok吗,不用先对wB87XD求值计算吗?

490

帖子

2

威望

5306

eV
积分
5836

Level 6 (一方通行)

11#
发表于 Post on 2018-5-30 16:47:02 | 只看该作者 Only view this author
qianbao 发表于 2018-5-30 15:08
楼主您好,我想请教一下,想在高斯用wB97XD泛函优化,是可以在命令行直接输入就ok吗,不用先对wB87XD求值计 ...

这种范围分离泛函,都是有默认的w的数值的。直接用就可以,你可以在输出文件找到默认的w值是多少。
估计你指的求值,是只针对某一体系,如何得到最优的w的值。这里所谓的最优,根据判断标准的不同,也有所不同。这就是经常听到的优化w的方法。具体实现见:
优化长程校正泛函w参数的简便工具optDFTw
http://sobereva.com/346

评分 Rate

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

查看全部评分 View all ratings

167

帖子

0

威望

1051

eV
积分
1218

Level 4 (黑子)

12#
发表于 Post on 2018-5-30 19:52:32 | 只看该作者 Only view this author
让你变成回忆 发表于 2018-5-30 16:47
这种范围分离泛函,都是有默认的w的数值的。直接用就可以,你可以在输出文件找到默认的w值是多少。
估计 ...

谢谢您!首先,我是想用此泛函像这篇文章这样优化,看是否可以得到同样的数据,所以我是可以在服务器中直接使用这个泛函和基组吗?(DOI: 10.1021/acs.chemmater.6b02930
其次,我还想问一下对于sob老师的这个博文中http://sobereva.com/346,这一处,我如果想计算的话,是需要这个界面体系编辑这个命令的,这样理解对吗?
使用optDFTw程序前首先要编辑一个Gaussian的长程校正泛函的单点任务文件作为模板,存到当前目录下template.gjf中,比如
%mem=60GB
%nproc=16
# LC-wPBE/6-311+G**

test

0 1
C                  0.00000000    0.00000000   -0.52710800
H                  0.00000000    0.93885600   -1.11413900
H                  0.00000000   -0.93885600   -1.11413900
O                  0.00000000    0.00000000    0.67386600

最后,非常感谢您!

微信截图_20180530194746.png (89.51 KB, 下载次数 Times of downloads: 76)

微信截图_20180530194746.png

微信截图_20180530194805.png (131.45 KB, 下载次数 Times of downloads: 86)

微信截图_20180530194805.png

490

帖子

2

威望

5306

eV
积分
5836

Level 6 (一方通行)

13#
发表于 Post on 2018-5-30 20:25:45 | 只看该作者 Only view this author
本帖最后由 让你变成回忆 于 2018-5-30 20:28 编辑
qianbao 发表于 2018-5-30 19:52
谢谢您!首先,我是想用此泛函像这篇文章这样优化,看是否可以得到同样的数据,所以我是可以在服务器中直 ...

(1)Bredas他们这篇CM的文章应该是直接用的wB97XD泛函,没有进行w的调控(至于行不行,有你问的时间,你都可以试出来究竟行不行了)。他们的很多文章都有调控过程,比如你在他们的文章中看到了sob的博文中提到的公式,通常就是在调控,根据J或者J^2,通常称为IP调控或者gap调控;
(2)我一般是在Linux下进行计算,所以不清楚你指的这个界面体系是啥意思。不过原理应该是一样的,即你需要准备一个template.gjf文件,这个文件中的分子结构部分就是你要进行计算的体系(通常是在一定的级别下优化好了的)。关键词就是一个单点任务的关键词(当然,如果你要在PCM下优化w,相应地加上PCM的关键词即可),因为所谓的调控过程就是不停地计算N、N+1、N-1电子体系单点能。

评分 Rate

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

查看全部评分 View all ratings

167

帖子

0

威望

1051

eV
积分
1218

Level 4 (黑子)

14#
发表于 Post on 2018-5-31 09:07:27 | 只看该作者 Only view this author
让你变成回忆 发表于 2018-5-30 20:25
(1)Bredas他们这篇CM的文章应该是直接用的wB97XD泛函,没有进行w的调控(至于行不行,有你问的时间,你 ...

非常感谢您的耐心回答!
首先,我已经通入到Linux下的服务器内,只是由于最初输入ω不识别导致出错,现已更换成w
其次,我是想根据这篇文章学习计算相互作用能的,而此前想先学会优化w;
最后,我的专业术语太匮乏,我自己说的界面体系其实就是文章中的PBDTTPD-C60;再次感谢您!

6万

帖子

99

威望

6万

eV
积分
125152

管理员

公社社长

15#
发表于 Post on 2018-5-31 17:40:52 | 只看该作者 Only view this author
qianbao 发表于 2018-5-31 09:07
非常感谢您的耐心回答!
首先,我已经通入到Linux下的服务器内,只是由于最初输入ω不识别导致出错,现 ...


计算相互作用能根本用不着优化w,就用wB97XD默认的就完了,调控w根本不是为了更好计算相互作用能的目的。计算弱相互作用的常识仔细看
谈谈“计算时是否需要加DFT-D3色散校正?”
http://sobereva.com/413
乱谈DFT-D
http://sobereva.com/83
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2026-2-23 18:01 , Processed in 0.208696 second(s), 28 queries , Gzip On.

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