计算化学公社

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

[VASP] 请教几个概念:偶极修正,零点能和范德华力。相应地在什么情况下需要考虑这几个条件?

[复制链接 Copy URL]

1

帖子

0

威望

25

eV
积分
26

Level 2 能力者

本帖最后由 pcczhouwei 于 2018-3-7 19:27 编辑

1、偶极修正:非对称的平板体系存在偶极矩。周期性的偶极矩形成一个假的场,这个场可能影响表面的一些性质,如功函数。解决这个问题的办法是在晶胞中引入一个偶极场来消除非零偶极矩,从而最大限度的减少这种偶极相互作用。(是不是可以这么理解,一个表面结构4层一个周期,我构建5层,上面表面是对称的是不是就不需要考虑这个偶极修正?还有什么情况下可以不用考虑这个?)
2、零点能:物质在绝对温度为零度下在真空中产生的能量。(MS中可以通过计算声子谱得到,vasp就不清楚了,一般什么情况下需要考虑零点能?)
3、范德华力:又称分子作用力,产生于 分子或原子之间的静电相互作用,相当于是一种弱相互作用力。(通常可以通过DFT-D色散修正进行计算,一般氢键好像就要考虑这个修正,不知道还有其它什么情况也需要用到这个修正?例如:O2、CO这些分子吸附的时候是否需要考虑?)
希望大伙儿能给予指点,上述所说纯属个人观点,有错误之处,希望大家指出,谢谢!

3621

帖子

3

威望

1万

eV
积分
18425

Level 6 (一方通行)

第一原理惨品小作坊

2#
发表于 Post on 2018-3-13 08:07:54 | 只看该作者 Only view this author
1、除了slab模型,对于有偶极矩的分子也应当三个方向进行偶极矩修正。对称或反对称的表面本身没有slab-slab的偶极相互作用因此可以不做修正,但并非所有结构都能这样,如果强行做成对称或反对称的就可能破坏化学计量比,比如你说的4层不是对称的但能够保证化学计量比,如果强行加到5层就可能会破坏化学计量比的情况,这种只能使用4层并考虑修正。

2、VASP也能计算声子频率来得到零点能。我并不认为光计算零点能有什么实际用处,借助VASP+Phonopy声子频率计算完成后能够同时得到声子谱(或声子态密度)、零点能、焓、自由能(原子核贡献),声子谱(或声子态密度)可以用来判断虚频情况,而焓和自由能与SCF能量加和可以得到体系在谐振近似下的焓和自由能。注意三点:(1)Phonopy本身为固体计算设计,因此算分子的时候热力学量不包含平动部分,需要自行修正。(2)VASP的计算是周期边界条件的,这个周期性边界条件本身破坏体系对称性。(3)谐振近似在高温情况下可能不是好的近似。

3、具体概念你可以参考http://sobereva.com/83,原则上,DFT-D对氢键和VDW作用的体系都有改善。

评分 Rate

参与人数
Participants 3
eV +6 收起 理由
Reason
ycl0426 + 2
tjuptz + 2
sobereva + 2

查看全部评分 View all ratings

日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。
本周忙

200

帖子

0

威望

1501

eV
积分
1701

Level 5 (御坂)

3#
发表于 Post on 2018-12-18 19:19:06 | 只看该作者 Only view this author
卡开发发 发表于 2018-3-13 08:07
1、除了slab模型,对于有偶极矩的分子也应当三个方向进行偶极矩修正。对称或反对称的表面本身没有slab-slab ...

请问单层2D材料掺杂或者缺陷体系需要考虑偶极校正吗?体系主要是用在光电催化研究方面

3621

帖子

3

威望

1万

eV
积分
18425

Level 6 (一方通行)

第一原理惨品小作坊

4#
发表于 Post on 2018-12-18 20:39:54 | 只看该作者 Only view this author
灰飞的旋律 发表于 2018-12-18 19:19
请问单层2D材料掺杂或者缺陷体系需要考虑偶极校正吗?体系主要是用在光电催化研究方面

原则上应该需要。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。
本周忙

200

帖子

0

威望

1501

eV
积分
1701

Level 5 (御坂)

5#
发表于 Post on 2018-12-24 11:17:16 | 只看该作者 Only view this author

那么对于自发极化的体系研究它的内建电场大小和载流子迁移率也同样需要??

3621

帖子

3

威望

1万

eV
积分
18425

Level 6 (一方通行)

第一原理惨品小作坊

6#
发表于 Post on 2018-12-24 16:29:56 | 只看该作者 Only view this author
灰飞的旋律 发表于 2018-12-24 11:17
那么对于自发极化的体系研究它的内建电场大小和载流子迁移率也同样需要??

z方向引入artificial的周期性引发的问题原则上肯定要消除,但不见得所有研究这个问题都很重要罢了,如果涉及到电势问题,其实还是尽量修正为好。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。
本周忙

200

帖子

0

威望

1501

eV
积分
1701

Level 5 (御坂)

7#
发表于 Post on 2018-12-24 19:43:36 | 只看该作者 Only view this author
卡开发发 发表于 2018-12-24 16:29
z方向引入artificial的周期性引发的问题原则上肯定要消除,但不见得所有研究这个问题都很重要罢了,如果 ...

嗯,我对比看看

7

帖子

0

威望

33

eV
积分
40

Level 2 能力者

8#
发表于 Post on 2021-5-6 22:16:43 | 只看该作者 Only view this author
卡开发发 发表于 2018-12-24 16:29
z方向引入artificial的周期性引发的问题原则上肯定要消除,但不见得所有研究这个问题都很重要罢了,如果 ...

卡哥,有个疑问。计算TiAl3(001)极性表面时,有文献为了消除偶极矩效应,取得7层但是破坏了化学计量比。那是不是可以构建4层或者8层的表面模型符合化学计量比然后施加偶极修正?第二问题就是极性表面怎么计算表面能?

3621

帖子

3

威望

1万

eV
积分
18425

Level 6 (一方通行)

第一原理惨品小作坊

9#
发表于 Post on 2021-5-7 09:11:52 | 只看该作者 Only view this author
魂牵梦绕 发表于 2021-5-6 22:16
卡哥,有个疑问。计算TiAl3(001)极性表面时,有文献为了消除偶极矩效应,取得7层但是破坏了化学计量比 ...
那是不是可以构建4层或者8层的表面模型符合化学计量比然后施加偶极修正?


第二问题就是极性表面怎么计算表面能?

符合计量比的情况直接按照表面能定义做。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。
本周忙

7

帖子

0

威望

33

eV
积分
40

Level 2 能力者

10#
发表于 Post on 2021-5-7 17:28:36 | 只看该作者 Only view this author
卡开发发 发表于 2021-5-7 09:11


符合计量比的情况直接按照表面能定义做。

大佬,目的:计算极性表面其中一个面的表面能  做法:第一步,切四层原子层,符合化学计量比,用公式算出两个面表面能的和;第二步,将优化前的四层原子层其中的A面两层原子固定且在A面悬挂氢原子钝化,利用公式只算一个面,试图算出单独B面的表面能;第三步将第二步中的A面换成B面再算一次,A面表面能。如果单独算出来的A面和B面的表面能之和与第一步总的表面能相差无几的话,是不是说这个方法是成立的?

3621

帖子

3

威望

1万

eV
积分
18425

Level 6 (一方通行)

第一原理惨品小作坊

11#
发表于 Post on 2021-5-7 19:20:35 | 只看该作者 Only view this author
魂牵梦绕 发表于 2021-5-7 17:28
大佬,目的:计算极性表面其中一个面的表面能  做法:第一步,切四层原子层,符合化学计量比,用公式算出 ...

利用公式如何算单个的面?表面总是成对出现的,计算上定义的表面能是将一个块体(可能是超晶格)劈裂成两个表面在单位面积上能量的变化。氢原子钝化的过程已经引入了额外的能量贡献。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。
本周忙

91

帖子

0

威望

5244

eV
积分
5335

Level 6 (一方通行)

12#
发表于 Post on 2022-11-25 16:31:39 | 只看该作者 Only view this author
卡开发发 发表于 2021-5-7 19:20
利用公式如何算单个的面?表面总是成对出现的,计算上定义的表面能是将一个块体(可能是超晶格)劈裂成两 ...

卡哥,VASP中OUTCAR里面读取dipolmoment数值是补偿偶极矩么,与实际值正好互为相反数么

30

帖子

0

威望

281

eV
积分
311

Level 3 能力者

13#
发表于 Post on 2023-8-7 12:24:35 | 只看该作者 Only view this author
Jack 发表于 2022-11-25 16:31
卡哥,VASP中OUTCAR里面读取dipolmoment数值是补偿偶极矩么,与实际值正好互为相反数么

同问,请问这个问题解决了吗

7

帖子

0

威望

41

eV
积分
48

Level 2 能力者

14#
发表于 Post on 2024-9-14 09:26:15 | 只看该作者 Only view this author
框架里的金属钠离子,吸附水分子算个吸附能需要加偶极校正吗,加了试了几次发现不收敛,删了就能收敛

本版积分规则 Credits rule

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

GMT+8, 2024-11-23 09:44 , Processed in 0.185052 second(s), 28 queries , Gzip On.

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