计算化学公社

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

[VASP] 求助VASP可以通过修改POTCAR来使原子带电吗

[复制链接 Copy URL]

30

帖子

0

威望

277

eV
积分
307

Level 3 能力者

各位论坛里的大佬,我想请问一下可以通过修改POTCAR赝势文件来让特定的原子带电子吗?我看过一篇文献里有说可以通过构建正负偶极板来模拟电场的作用,如下图所示,那么应该怎样在VASP里面做到让原子带电呢?麻烦各位给小弟指点一下迷津,谢谢大家

图1.png (41.99 KB, 下载次数 Times of downloads: 63)

图1.png

33

帖子

0

威望

124

eV
积分
157

Level 3 能力者

2#
发表于 Post on 2020-3-22 19:02:57 | 只看该作者 Only view this author
VASP下修改POTCAR和修改NELEC一样都是background charge. 具体到原子所带电荷只能是算出来啥样就是啥样。https://cms.mpi.univie.ac.at/vasp/vasp/NELECT.html

30

帖子

0

威望

277

eV
积分
307

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2020-3-22 20:37:21 | 只看该作者 Only view this author
谢谢大佬解惑,那我问下cp2k正常来说可以做到在原子上增减电子来使其带电吗?

3754

帖子

3

威望

1万

eV
积分
19675

Level 6 (一方通行)

围观吃瓜群众

4#
发表于 Post on 2020-3-22 21:14:36 | 只看该作者 Only view this author
MADAOLay 发表于 2020-3-22 20:37
谢谢大佬解惑,那我问下cp2k正常来说可以做到在原子上增减电子来使其带电吗?

一般这种模型可能都是通过点电荷来实现的,如果有程序支持识别点电荷,原则上可以做到。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

471

帖子

0

威望

1795

eV
积分
2266

Level 5 (御坂)

5#
发表于 Post on 2020-3-22 21:34:51 | 只看该作者 Only view this author
就我所知,拿现成POTCAR修改,没法实现你的想法,而且这样做实际意义存疑.....
顺便一问,图片有意思,能说下来源么?

30

帖子

0

威望

277

eV
积分
307

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2020-3-23 13:46:14 | 只看该作者 Only view this author
hakuna 发表于 2020-3-22 21:34
就我所知,拿现成POTCAR修改,没法实现你的想法,而且这样做实际意义存疑.....
顺便一问,图片有意思,能 ...

来自论文 Activation of CO2 on Copper Surfaces: The Synergy between Electric Field, Surface Morphology and Excess Electrons

30

帖子

0

威望

277

eV
积分
307

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2020-3-23 13:52:32 | 只看该作者 Only view this author
各位老师们,我再问下vasp还有其他方法可以做到对原子带电量的控制吗?其他软件也可以,各位有接触过类似的情况或者有类似的方法可以给小弟指点一下吗

3754

帖子

3

威望

1万

eV
积分
19675

Level 6 (一方通行)

围观吃瓜群众

8#
发表于 Post on 2020-3-23 14:27:51 | 只看该作者 Only view this author
MADAOLay 发表于 2020-3-23 13:52
各位老师们,我再问下vasp还有其他方法可以做到对原子带电量的控制吗?其他软件也可以,各位有接触过类似的 ...

1、电荷应该不是直接在原子上的,而是以外加电荷考虑的。
2、你发来的原文是CP2K做的,可以研究一下相关资料。
3、我了解的情况,QE挂Environ插件(https://github.com/environ-developers/Environ)也能做,原则上应该是设置EXTERNAL_CHARGES下面的参数,其中dim=2,外电荷的形式应该是一个Gaussian型的平板,你可以具体看看手册。

评分 Rate

参与人数
Participants 1
eV +4 收起 理由
Reason
MADAOLay + 4 谢谢

查看全部评分 View all ratings

日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

30

帖子

0

威望

277

eV
积分
307

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2020-4-16 23:50:10 | 只看该作者 Only view this author
卡开发发 发表于 2020-3-22 21:14
一般这种模型可能都是通过点电荷来实现的,如果有程序支持识别点电荷,原则上可以做到。

老师,时隔这么久再来问一下,模型中点电荷的添加是通过建模软件比如materials studio之类的还是通过在仿真软件中直接设置的?

3754

帖子

3

威望

1万

eV
积分
19675

Level 6 (一方通行)

围观吃瓜群众

10#
发表于 Post on 2020-4-17 09:47:04 | 只看该作者 Only view this author
MADAOLay 发表于 2020-4-16 23:50
老师,时隔这么久再来问一下,模型中点电荷的添加是通过建模软件比如materials studio之类的还是通过在仿 ...

QE的话挂了Environ的插件是可以直接通过设置参数给定一个Gaussian型的电荷面,不需要其他设置。MS的DMol3虽然能以点电荷来做,但是我不确定有没有什么其他的问题。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

18

帖子

0

威望

313

eV
积分
331

Level 3 能力者

11#
发表于 Post on 2020-4-18 13:05:02 | 只看该作者 Only view this author
VASP 修改价电子数可以带电 NELEM

14

帖子

0

威望

459

eV
积分
473

Level 3 能力者

12#
发表于 Post on 2020-6-22 12:07:19 | 只看该作者 Only view this author
aaronzjw 发表于 2020-4-18 13:05
VASP 修改价电子数可以带电 NELEM

好像没有查到这个参数,怎么使用呢?

15

帖子

0

威望

135

eV
积分
150

Level 3 能力者

13#
发表于 Post on 2020-6-22 18:56:50 | 只看该作者 Only view this author
ykr 发表于 2020-6-22 12:07
好像没有查到这个参数,怎么使用呢?

是NELECT哦,vasp手册的P73

23

帖子

0

威望

227

eV
积分
250

Level 3 能力者

14#
发表于 Post on 2021-10-9 22:16:42 | 只看该作者 Only view this author
请问各位老师,CP2K如何让原子带电?有类似的修改价电子数的参数或者方法吗?

23

帖子

0

威望

227

eV
积分
250

Level 3 能力者

15#
发表于 Post on 2021-10-9 22:25:37 | 只看该作者 Only view this author
卡开发发 发表于 2020-3-22 21:14
一般这种模型可能都是通过点电荷来实现的,如果有程序支持识别点电荷,原则上可以做到。

卡老师,我在CP2K的论坛上看到有人解答如何建立电场,他说“you can define a ghost atom with a charge and no basis functions.” 请问这应该如何实现呢?是需要在基组文件和赝势文件中自己写一个带正电的氢吗?

还想请问,在CP2K中如果想让某些具体的原子带电,比如N原子,是否有办法?能否像VASP那样修改价电子数?

感谢!

本版积分规则 Credits rule

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

GMT+8, 2025-8-15 21:02 , Processed in 0.183764 second(s), 24 queries , Gzip On.

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