| 查看: 5338 | 回復(fù): 22 | ||||
[求助]
VASP計(jì)算時(shí)需要考慮范德華力,我的體系有鉈元素,對(duì)應(yīng)的VDW_C6= 和VDW_R0=值是多少急求 已有2人參與
|
||||
|
VASP計(jì)算時(shí)需要考慮范德華力,我的體系有鉈元素,對(duì)應(yīng)的LVDW=.TRUE. VDW_RADIUS=30.0 VDW_SCALING=0.75 VDW_D=20.0 VDW_C6= 和 VDW_R0=值是多少?急求,謝謝!!!!!! |
VASP | 科研 |
金蟲 (小有名氣)
|
用我前面說的方法做就可以了吧。。。。 其實(shí)我發(fā)現(xiàn)新版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/) |
金蟲 (小有名氣)
|
http://www.thch.uni-bonn.de/tc/i ... mp;amp;lang=english 網(wǎng)頁中點(diǎn)擊: Tarball of the current versions of the code and coefficient file (V3.1 Rev 0 ) 可以下載到DFT-D3 的 source code 解壓縮以后,編譯可以得到dftd3的binary 然后把你要計(jì)算的POSCAR拷過來,運(yùn)行 dftd3 POSCAR -pbc -func pbe -old (相關(guān)參數(shù)設(shè)置見http://www.thch.uni-bonn.de/tc/downloads/DFT-D3/data/man.pdf 最后幾頁有寫,-pbc是指使用POSCAR的結(jié)構(gòu)文件,-func pbe指你使用PBE的functional,-old指切換到DFT-D2的模式) 然后可以得到C6和R0,但是貌似單位有些許問題,你可以對(duì)照http://cms.mpi.univie.ac.at/vasp/vasp/DFT_D2_method_Grimme.html 里提供的表格矯正以下單位,例如: 程序給出H的C6=2.43,然而網(wǎng)頁上是C6_H=0.14,程序給出的Tl是993.09,網(wǎng)頁沒有給出相關(guān)數(shù)值,于是C6_Tl=(0.14/2.43)*993.09=57.22 R0用類似的方法計(jì)算即可。注意一下,C6最好用上面Z=xx那里的數(shù)值,因?yàn)橄旅娴谋砀駜?nèi)的只有一位小數(shù),可能計(jì)算出來有誤差 單位是怎么變的,我也沒搞清楚,但用這種辦法進(jìn)行校正的話是可以給出正確數(shù)值的,我測(cè)試過H,C,N 都給出了網(wǎng)頁上的正確數(shù)值 |
|
我在網(wǎng)上查到, 發(fā)現(xiàn)DFT-D3的單位是atomic unit, 所以長(zhǎng)度是Bohr, 能量是Hartree. DFT-D3給的C6單位是Hartree*Bohr, 而VASP網(wǎng)站上的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的單位轉(zhuǎn)換為 993.09*0.0576526021 = 57.2542226 (J/mol)*(nm)^6 |
金蟲 (小有名氣)
木蟲 (小有名氣)
木蟲 (小有名氣)
|
不需要設(shè)置那些數(shù)值,只需要INCAR里面IVDW= 11 http://cms.mpi.univie.ac.at/vasp/vasp/DFT_D3_method.html |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工306分找調(diào)劑 +6 | 滄海輕舟e 2026-04-02 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 085600,320分求調(diào)劑 +5 | 大饞小子 2026-04-02 | 5/250 |
|
|
[考研] 初試301,代碼085701環(huán)境工程,本碩一致,四六級(jí)已過,有二區(qū)一作,共發(fā)表5篇論文 +6 | axibli 2026-04-01 | 6/300 |
|
|
[考研] 337求調(diào)劑 +11 | 《樹》 2026-03-29 | 11/550 |
|
|
[考研] 085600,材料與化工321分調(diào)劑 +8 | 大饞小子 2026-03-27 | 10/500 |
|
|
[考研] 070300化學(xué)求調(diào)劑 +14 | 小黃鴨寶 2026-03-30 | 14/700 |
|
|
[考研] 各位老師好,我的一志愿為北京科技大學(xué)085601材料專碩 +13 | Koxui 2026-03-28 | 13/650 |
|
|
[考研] 求調(diào)劑0703 +5 | 周嘉堯 2026-03-31 | 8/400 |
|
|
[考研] 085600 一志愿9 總分351 求調(diào)劑學(xué)校 +7 | czhcz 2026-03-31 | 9/450 |
|
|
[考研] 086000生物與醫(yī)藥298調(diào)劑求助 +4 | 元元青青 2026-03-31 | 6/300 |
|
|
[考研] 326求調(diào)劑 +4 | 崽崽仔 2026-03-31 | 4/200 |
|
|
[考研] 311(085601)求調(diào)劑 +12 | liziyeyeye 2026-03-28 | 13/650 |
|
|
[考研] 335求調(diào)劑 +3 | 321* 2026-03-31 | 4/200 |
|
|
[考研] 材料工程085601數(shù)二英一335求調(diào)劑 +5 | 雙馬尾痞老板2 2026-03-31 | 5/250 |
|
|
[考研] 353求調(diào)劑 +3 | 江上楓_26 2026-03-28 | 3/150 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-31 | 4/200 |
|
|
[考研] 南京大學(xué)化學(xué)調(diào)劑 +11 | 景隨風(fēng) 2026-03-29 | 16/800 |
|
|
[考研] 322求調(diào)劑 +10 | 宋明欣 2026-03-27 | 10/500 |
|
|
[考研] 298求調(diào)劑 +4 | 種圣賜 2026-03-28 | 4/200 |
|
|
[考研] 394求調(diào)劑 +3 | 好事多磨靜候佳?/a> 2026-03-26 | 5/250 |
|