请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 1775|回复: 8

[VASP] 如何验证vaspsol是否编译成功

[复制链接]

407

帖子

0

威望

813

eV
积分
1220

Level 4 (黑子)

发表于 2018-11-10 18:55:17 | 显示全部楼层 |阅读模式
本帖最后由 小书童 于 2018-11-10 19:33 编辑

首先,如何验证vaspsol是否编译成功呢?在OUTCAR中会有什么关键词提示吗?
第二,对于vasosol的使用上,USAGE.md上说首先要 Do a vacuum calculation for your system然后存储波函数文件,这具体应该怎么操作呢?
第三,cavitation energy 这是指什么呢?
还望懂得的前辈们指点下,多谢多谢。

1773

帖子

3

威望

6626

eV
积分
8459

Level 6 (一方通行)

Ab Initio Amateur

发表于 2018-11-10 23:13:52 | 显示全部楼层
本帖最后由 卡开发发 于 2018-11-10 23:16 编辑

1、我没用过这个东西,但看了下例子https://github.com/henniggroup/VASPsol/blob/master/examples/CO/Solvation/INCAR
例子里面开了:
LSOL = .TRUE.
EB_K = 80
之后OSZICAR(https://github.com/henniggroup/V ... lvation/OSZICAR.ref)的显示除了DAV还有SOL。
2、真空计算的时候LWAVE = .TRUE.(这已经是默认的了,所以不用写),然后做溶剂化计算的时候开启ISTART=1
3、http://sobereva.com/42

评分

参与人数 2eV +7 收起 理由
小书童 + 5 谢谢
sobereva + 2

查看全部评分

满招损,谦受益。热衷于理论和方法研究水平不高但欢迎讨论。

47

帖子

0

威望

189

eV
积分
236

Level 3 能力者

发表于 2019-3-23 23:25:49 | 显示全部楼层
您好,我刚使用VASP,想请教一下VASPsol在服务器上如何编译,服务器是天津超算,谢谢

79

帖子

1

威望

759

eV
积分
858

Level 4 (黑子)

发表于 2019-3-24 13:36:55 | 显示全部楼层
F凌兰 发表于 2019-3-23 23:25
您好,我刚使用VASP,想请教一下VASPsol在服务器上如何编译,服务器是天津超算,谢谢

首先确认可以正常编译原版的vasp,下载:
https://github.com/henniggroup/VASPsol/blob/master/src/solvation.F
再重新编译就行了。方法如下:
微信截图_20190324102255.png

评分

参与人数 1eV +5 收起 理由
小书童 + 5 谢谢

查看全部评分

79

帖子

1

威望

759

eV
积分
858

Level 4 (黑子)

发表于 2019-3-24 13:37:39 | 显示全部楼层
本帖最后由 啦啦黑还黑 于 2019-3-24 13:40 编辑

成功编译后INCAR只需要加LSOL = T就是水溶液的溶剂化。OSZICAR里能看到每一步电子步都有一步溶剂化的计算:
  1.        N       E                     dE             d eps       ncg     rms          rms(c)
  2. SOL:   1  -0.27575E+00    0.32189E-01   -0.24356E+00   109
  3. RMM:   1    -0.145388318781E+02   -0.14539E+02   -0.46978E-01   826   0.253E+00    0.609E-01
  4. SOL:   2  -0.30548E+00    0.32137E-01   -0.27334E+00    98
  5. RMM:   2    -0.145312080255E+02    0.76239E-02   -0.36526E-02   891   0.104E+00    0.287E-01
  6. SOL:   3  -0.34279E+00    0.31882E-01   -0.31091E+00   100
  7. RMM:   3    -0.145346598726E+02   -0.34518E-02   -0.19833E-02   956   0.790E-01    0.177E-01
  8. SOL:   4  -0.38999E+00    0.31644E-01   -0.35834E+00   102
  9. RMM:   4    -0.145356933928E+02   -0.10335E-02   -0.12221E-02   871   0.630E-01    0.131E-01
  10. SOL:   5  -0.42000E+00    0.31572E-01   -0.38843E+00    99
  11. RMM:   5    -0.145346100972E+02    0.10833E-02   -0.63980E-03   984   0.444E-01    0.545E-02
  12. SOL:   6  -0.43956E+00    0.31622E-01   -0.40794E+00    98
  13. RMM:   6    -0.145349081511E+02   -0.29805E-03   -0.31742E-03   742   0.344E-01    0.477E-02
复制代码
换其他溶剂需要查表换EB_k, SIGMA_K, NC_K, TAU这些参数,
主要的就是介电常数:EB_k,vaspsol默认的是水,用的78.4

47

帖子

0

威望

189

eV
积分
236

Level 3 能力者

发表于 2019-3-26 11:27:23 | 显示全部楼层
啦啦黑还黑 发表于 2019-3-24 13:36
首先确认可以正常编译原版的vasp,下载:
https://github.com/henniggroup/VASPsol/blob/master/src/sol ...

好的,非常感谢

46

帖子

0

威望

803

eV
积分
849

Level 4 (黑子)

发表于 2019-9-6 10:27:37 | 显示全部楼层
您好,请问一下我调用了vasp.5.4.1.sol计算得到的OSZICAR里还是只有DAV行,没有新的SOL行,您知道这是什么原因吗?

62

帖子

0

威望

215

eV
积分
277

Level 3 能力者

发表于 2020-1-6 21:49:49 | 显示全部楼层
啦啦黑还黑 发表于 2019-3-24 13:36
首先确认可以正常编译原版的vasp,下载:
https://github.com/henniggroup/VASPsol/blob/master/src/sol ...

您给的连接我没找到下载的地方,麻烦您能发我一份吗?

34

帖子

0

威望

222

eV
积分
256

Level 3 能力者

发表于 2020-1-13 16:40:37 | 显示全部楼层
下午茶 发表于 2020-1-6 21:49
您给的连接我没找到下载的地方,麻烦您能发我一份吗?

https://github.com/henniggroup/VASPsol

其实5楼那个链接返回就是他了
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2020-2-23 16:10 , Processed in 0.166241 second(s), 27 queries .

快速回复 返回顶部 返回列表