VASP计算时需要考虑范德华力，我的体系有铊元素,对应的VDW_C6= 和VDW_R0=值是多少急求 - 第一性原理 - Vasp&MedeA

11楼2015-05-21 21:57:47

微6凉

铁虫 (小有名气)

应助: 1 (幼儿园)
金币: 202.8
帖子: 89
在线: 35.6小时
虫号: 2724445
注册: 2013-10-14
专业: 金属非晶态、准晶和纳米晶

引用回帖:

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:// ...

想问一下，您是怎么找到这个参数的，我怎么都没找到

12楼2015-07-28 10:04:01

微6凉

铁虫 (小有名气)

应助: 1 (幼儿园)
金币: 202.8
帖子: 89
在线: 35.6小时
虫号: 2724445
注册: 2013-10-14
专业: 金属非晶态、准晶和纳米晶

引用回帖:

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...

想问一下，你用这个方法算过结果么？怎么样

13楼2015-07-28 10:05:13

微6凉

铁虫 (小有名气)

应助: 1 (幼儿园)
金币: 202.8
帖子: 89
在线: 35.6小时
虫号: 2724445
注册: 2013-10-14
专业: 金属非晶态、准晶和纳米晶

引用回帖:

请问pt的这些值是什么呢？希望能给与帮助

14楼2015-07-28 11:34:50

awdawdawdd

金虫 (小有名气)

应助: 24 (小学生)
金币: 679
散金: 200
红花: 6
帖子: 111
在线: 67小时
虫号: 2237233
注册: 2013-01-12
专业: 凝聚态物性 II ：电子结构

引用回帖:

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

微6凉

铁虫 (小有名气)

应助: 1 (幼儿园)
金币: 202.8
帖子: 89
在线: 35.6小时
虫号: 2724445
注册: 2013-10-14
专业: 金属非晶态、准晶和纳米晶

引用回帖:

15楼: Originally posted by awdawdawdd at 2015-07-28 12:32:19
用我前面说的方法做就可以了吧。。。。
其实我发现新版VASP里面集成了后面那些元素的值了，你可以看看源代码里面有(vdwforcefield.F)：

1234 FORMAT(/," DFT-D2 method for vdW energy calculation&quo ...

因为是初学，对这部分还不熟悉，所以看不懂这些，不过你的这些数值对我来说有很大的帮助，非常感谢

16楼2015-07-28 18:35:13

sulisha_agrh

新虫 (小有名气)

应助: 2 (幼儿园)
金币: 30.4
散金: 324
红花: 1
帖子: 88
在线: 273.8小时
虫号: 3638590
注册: 2015-01-10
性别: MM
专业: 凝聚态物性I:结构、力学和

引用回帖:

为何不用vdw－DF

17楼2015-08-02 10:02:43

limurphy

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 31
红花: 1
帖子: 19
在线: 18.1小时
虫号: 2152990
注册: 2012-11-28
专业: 凝聚态物性 II ：电子结构

引用回帖:

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