应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP计算时需要考虑范德华力，我的体系有铊元素,对应的VDW_C6= 和VDW_R0=值是多少急求
231/3返回列表123下一页
查看: 5671  |  回复: 22
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

巨蟹水蜜桃

铜虫 (初入文坛)

[求助] VASP计算时需要考虑范德华力，我的体系有铊元素,对应的VDW_C6= 和VDW_R0=值是多少急求 已有2人参与

VASP计算时需要考虑范德华力，我的体系有铊元素，对应的LVDW=.TRUE.
VDW_RADIUS=30.0
VDW_SCALING=0.75
VDW_D=20.0
VDW_C6= 和
VDW_R0=值是多少？急求，谢谢!!!!!!
回复此楼

» 收录本帖的淘帖专辑推荐

VASP 科研

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-05-20 15:31:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

awdawdawdd

金虫 (小有名气)
引用回帖:
14楼: Originally posted by 微6凉 at 2015-07-28 11:34:50
请问pt的这些值是什么呢？希望能给与帮助...

用我前面说的方法做就可以了吧。。。。
其实我发现新版VASP里面集成了后面那些元素的值了，你可以看看源代码里面有(vdwforcefield.F)：

1234 FORMAT(/,"  DFT-D2 method for vdW energy calculation",&
              /,"  -------------------------------------------------------------------"
  1235 FORMAT(/,"  Parameters of vdW forcefield:",&
              /,"         C6(Jnm^6/mol)     R0(A)",&
              /,"   -----------------------------"
  1236 FORMAT("   ",A2,5X,F7.3,8X,F7.3)
  2250 FORMAT(/,"  IVDW = ",I3)
  1237 FORMAT("  VDW_RADIUS = ",F9.3," A"
  1233 FORMAT("  VDW_S6 = ",F9.3)
  1241 FORMAT("  VDW_SR = ",F9.3)
  1240 FORMAT("  VDW_D = ",F9.3)
  2246 FORMAT("  LVDW_EWALD = ",L1)
        

          if (historycounter==1) then
            ALLOCATE(cc_(T_INFO%NTYP),r0_(T_INFO%NTYP))
            tags=(/'H ','He','Li','Be','B ','C ','N ','O ','F ','Ne',&
                   'Na','Mg','Al','Si','P ','S ','Cl','Ar','K ','Ca',&
                   'Sc','Ti','V ','Cr','Mn','Fe','Co','Ni','Cu','Zn',&
                   'Ga','Ge','As','Se','Br','Kr','Rb','Sr',&
                   'Y ','Zr','Nb','Mo','Tc','Ru','Rh','Pd','Ag','Cd',&
                   'In','Sn','Sb','Te','I ','Xe','X ', &
                   'Cs','Ba','La','Ce','Pr','Nd','Pm',&
                   'Sm','Eu','Gd',&
                   'Tb','Dy','Ho','Er','Tm','Yb','Lu',&
                   'Hf','Ta','W ','Re','Os','Ir','Pt','Au','Hg',&
                   'Tl','Pb','Bi','Po','At','Rn'/)

            cc=(/0.14,0.08,1.61,1.61,3.13,1.75,1.23,0.70,0.75,0.63,&
                 5.71,5.71,10.79,9.23,7.84,5.57,5.07,4.61,10.8,10.8,&
                 10.8,10.8,10.8,10.8,10.8,10.8,10.8,10.8,10.8,10.8,&
                 16.99,17.10,16.37,12.64,12.47,12.01,24.67,24.67,&
                 24.67,24.67,24.67,24.67,24.67,24.67,24.67,24.67,24.67,24.67,&
                 37.32,38.71,38.44,31.74,31.50,29.99,29.99, &
                 315.275,226.994,176.252,140.68,140.68,140.68,140.68,&
                 140.68,140.68,140.68, &
                 140.68,140.68,140.68,140.68,140.68,140.68,140.68, &
                 105.112,81.24,81.24,81.24,81.24,81.24,81.24,81.24,57.364,&
                 57.254,63.162,63.540,55.283,57.171,56.64/)
            r0=(/1.001,1.012,0.825,1.408,1.485,1.452,1.397,1.342,1.287,1.243,&
                 1.144,1.364,1.639,1.716,1.705,1.683,1.639,1.595,1.485,1.474,&
                 1.562,1.562,1.562,1.562,1.562,1.562,1.562,1.562,1.562,1.562,&
                 1.650,1.727,1.760,1.771,1.749,1.727,1.628,1.606,&
                 1.639,1.639,1.639,1.639,1.639,1.639,1.639,1.639,1.639,1.639,&
                 1.672,1.804,1.881,1.892,1.892,1.881,1.881,&              
                 1.802,1.762,1.720,1.753,1.753,1.753,1.753,&
                 1.753,1.753,1.753,&
                 1.753,1.753,1.753,1.753,1.753,1.753,1.753,&
                 1.788,1.772,1.772,1.772,1.772,1.772,1.772,1.772,1.758,&
                 1.989,1.944,1.898,2.005,1.991,1.924/)
赞一下(3人) 回复此楼
15楼2015-07-28 12:32:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qlihytc

新虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
VASP manual 里面应该有个表格，表格里面没有就到原始文献（手册上就有文献信心）看看有没有，原始文献也没有的话就没有了。但是5.3以上貌似增加了D3
赞一下(1人) 回复此楼
2楼2015-05-20 15:37:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

awdawdawdd

金虫 (小有名气)
引用回帖:
10楼: Originally posted by 巨蟹水蜜桃 at 2015-05-21 16:39:43
@awdawdawdd应助我的Tl的范德华力参数C6=57.22 R0=1.986值跟你所说的vasp5.3.5里边的值一样吗...

http://www.thch.uni-bonn.de/tc/i ... mp;amp;lang=english
网页中点击：
Tarball of the current versions of the code and coefficient file (V3.1 Rev 0 )
可以下载到DFT-D3 的 source code
解压缩以后，编译可以得到dftd3的binary
然后把你要计算的POSCAR拷过来，运行 dftd3 POSCAR -pbc -func pbe -old （相关参数设置见http://www.thch.uni-bonn.de/tc/downloads/DFT-D3/data/man.pdf 最后几页有写,-pbc是指使用POSCAR的结构文件，-func pbe指你使用PBE的functional，-old指切换到DFT-D2的模式）

然后可以得到C6和R0，但是貌似单位有些许问题，你可以对照http://cms.mpi.univie.ac.at/vasp/vasp/DFT_D2_method_Grimme.html 里提供的表格矫正以下单位，例如：
程序给出H的C6=2.43，然而网页上是C6_H=0.14，程序给出的Tl是993.09，网页没有给出相关数值，于是C6_Tl=(0.14/2.43)*993.09=57.22
R0用类似的方法计算即可。注意一下，C6最好用上面Z=xx那里的数值，因为下面的表格内的只有一位小数，可能计算出来有误差

单位是怎么变的，我也没搞清楚，但用这种办法进行校正的话是可以给出正确数值的，我测试过H,C,N 都给出了网页上的正确数值
赞一下(1人) 回复此楼
11楼2015-05-21 21:57:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

limurphy

铁虫 (初入文坛)
引用回帖:
11楼: Originally posted by awdawdawdd at 2015-05-21 21:57:47
http://www.thch.uni-bonn.de/tc/index.php?section=downloads&subsection=getd3&lang=english
网页中点击：
Tarball of the current versions of the code and coefficient file (V3.1 Rev 0 )
可以下 ...

我在網上查到, 發現DFT-D3的單位是atomic unit, 所以長度是Bohr, 能量是Hartree.

DFT-D3給的C6單位是Hartree*Bohr, 而VASP網站上的C6單位是(J/mol)*(nm)^6
1 Hartree = 2625499.39949012 J/mol
1 Bohr = 0.052917725 nm
所以 1 Hartree*Bohr = 0.0576526021 (J/mol)*(nm)^6

而以 Tl為例, C6的單位轉換為 993.09*0.0576526021 = 57.2542226 (J/mol)*(nm)^6
赞一下(1人) 回复此楼
18楼2015-12-03 07:08:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

巨蟹水蜜桃

铜虫 (初入文坛)
引用回帖:
2楼: Originally posted by qlihytc at 2015-05-20 15:37:20
VASP manual 里面应该有个表格，表格里面没有就到原始文献（手册上就有文献信心）看看有没有，原始文献也没有的话就没有了。但是5.3以上貌似增加了D3

我没有找到啊你可以帮我找一下吗？谢谢
赞一下 回复此楼
3楼2015-05-20 21:02:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

awdawdawdd

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
巨蟹水蜜桃: 金币+15, ★★★很有帮助 2015-05-21 09:15:54
Tl
C6=57.22
R0=1.986

ref:
http://www.thch.uni-bonn.de/tc/i ... mp;amp;lang=english
http://www.thch.uni-bonn.de/tc/downloads/DFT-D3/data/man.pdf
http://cms.mpi.univie.ac.at/vasp/vasp/DFT_D2_method_Grimme.html
赞一下 回复此楼
4楼2015-05-20 23:29:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

巨蟹水蜜桃

铜虫 (初入文坛)
引用回帖:
4楼: Originally posted by awdawdawdd at 2015-05-20 23:29:08
Tl
C6=57.22
R0=1.986

ref:
http://www.thch.uni-bonn.de/tc/index.php?section=downloads&subsection=getd3&lang=english
http://www.thch.uni-bonn.de/tc/downloads/DFT-D3/data/man.pdf
http:// ...

我打开看了链接，也没有找到Tl的范德华力参数
赞一下 回复此楼
5楼2015-05-21 09:22:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuraining

木虫 (小有名气)
用vasp5.3.5，自带vdw-dft-d3，里面所有的元素vdw参数默认都有。
赞一下 回复此楼
6楼2015-05-21 11:50:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

巨蟹水蜜桃

铜虫 (初入文坛)
引用回帖:
6楼: Originally posted by yuraining at 2015-05-21 11:50:25
用vasp5.3.5，自带vdw-dft-d3，里面所有的元素vdw参数默认都有。

没有用过，你可以帮我截个图发过来吗？谢谢O(∩_∩)O谢谢
赞一下 回复此楼
7楼2015-05-21 13:47:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuraining

木虫 (小有名气)
引用回帖:
7楼: Originally posted by 巨蟹水蜜桃 at 2015-05-21 13:47:32
没有用过，你可以帮我截个图发过来吗？谢谢O(∩_∩)O谢谢...

不需要设置那些数值，只需要INCAR里面IVDW= 11
http://cms.mpi.univie.ac.at/vasp/vasp/DFT_D3_method.html
赞一下 回复此楼
8楼2015-05-21 15:12:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

巨蟹水蜜桃

铜虫 (初入文坛)
引用回帖:
8楼: Originally posted by yuraining at 2015-05-21 15:12:07
不需要设置那些数值，只需要INCAR里面IVDW= 11
http://cms.mpi.univie.ac.at/vasp/vasp/DFT_D3_method.html...

我们使用的是vasp5.2，你给我发的链接我找不到啊。
赞一下 回复此楼
9楼2015-05-21 16:33:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

巨蟹水蜜桃

铜虫 (初入文坛)
引用回帖:
9楼: Originally posted by 巨蟹水蜜桃 at 2015-05-21 16:33:47
我们使用的是vasp5.2，你给我发的链接我找不到啊。...

@awdawdawdd应助我的Tl的范德华力参数C6=57.22 R0=1.986值跟你所说的vasp5.3.5里边的值一样吗
赞一下 回复此楼
10楼2015-05-21 16:39:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 巨蟹水蜜桃 的主题更新
231/3返回列表123下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 可靠消息，8月21号(明天)放榜 +3 医学老男孩 2026-08-20 3/150 2026-08-20 09:50 by 突突小宝贝
[基金申请] 只有每年这种时候来逛逛小木虫 +20 yaoyewhu2008 2026-08-20 20/1000 2026-08-20 09:41 by 2026年国自然基�
[基金申请] 今天基金会出结果吗？20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[基金申请] 时间戳变了，能看出什么问题？ +12 基诺咪客 2026-08-17 16/800 2026-08-19 22:04 by kk 的小木虫
[基金申请] filecode，4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低？主要原因有哪些？ +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗？ +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 重要消息，中午系统在维护 +10 yuleib84 2026-08-18 11/550 2026-08-19 07:30 by 羽毛枫f
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护，明天很可能放榜！ +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
关闭
请填处理意见
关闭