Prettyswan

木蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 2919
散金: 149
紅花: 4
帖子: 145
在線: 59.7小時(shí)
蟲號(hào): 958555
注冊(cè): 2010-03-01
性別: MM
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[交流] 【轉(zhuǎn)帖】POTCAR文件詳解已有8人參與

近期，我對(duì)POTCAR的格式有點(diǎn)茫然，在小木蟲論壇上也沒有找到對(duì)于POTCAR的解讀。今天我在http://blog.chinaunix.net/u2/65114/showart_574445.html找到了一個(gè)關(guān)于POTCAR詳解的博文，轉(zhuǎn)來給大家共同學(xué)習(xí)，下面是博文全文：

引用回帖:

VASP贗勢(shì)文件POTCAR詳解
  POTCAR 贗勢(shì)文件可以理解為分子力學(xué)模擬中的力場(chǎng)文件，但包括的信息更多。VASP4.6將各元素優(yōu)化的INCAR里的參數(shù)也包括在這里了，作為支持PREC的缺省選擇。通常各元素的POTCAR已經(jīng)包括在軟件包里了
  我們只需要按照POSCAR里的順序，將各元素的POTCAR按順序連接起來就可以了，如以下命令：
  cat file1 file2 file3 > POTCAR

  軟件包自帶的絕大多數(shù)贗勢(shì)是超軟贗勢(shì)（US-PP）了，但不少元素有兩個(gè)版本，如何選取呢？
  一個(gè)簡(jiǎn)單的辦法是看后綴：標(biāo)準(zhǔn)的沒有后綴 _h 硬一點(diǎn) _s 軟一點(diǎn)。_pv,_sv,_d 就是說semi-core的p,s或者d也當(dāng)做價(jià)態(tài)處理了。如果是數(shù)字的話，表示的可能是不同的半徑截距。

  也可以參考各版本同目錄下的V_RHFIN file、PSCTR file，這兩個(gè)文件告知該版本的贗勢(shì)是如何生成的。比如：
  V_RHFIN file
  Sc: 6p d2 s1
8 21. .002000 44.95590 125. .25E-05 .300 200FCA 12.00000
.7 1.0 0
1.0 .0 .5 -320.8847 2.0000
2.0 .0 .5 -34.4217 2.0000
2.0 1.0 1.5 -28.2366 6.0000
3.0 .0 .5 -3.7944 2.0000
3.0 1.0 1.5 -2.2591 6.0000
3.0 2.0 2.5 -.1113 2.0000
4.0 .0 .5 -.2699 1.0000
4.0 3.0 2.5 -.1000 .0000
  第一行是注釋行，給出基本的信息
  第二行是最重要的控制行
8 21. .00 2000 44.95590 125. .25E-05 .300 200 F CA 12.00000
J Z XION N AM H DELRVR PHI NC1 | CH QCOR | GREEN

  J - 軌道數(shù)
  Z - 原子序數(shù)
  XION - 離子化程度一般設(shè)為0
  N - 格點(diǎn)數(shù)
  AM - 原子質(zhì)量
  H - 決定格點(diǎn)間距
  DELRVR - 自洽收斂標(biāo)準(zhǔn)
  PHI - 線性擬合參數(shù)
  NC1 - 最大自洽循環(huán)次數(shù)
  GREEN - 是否存在初始的勢(shì)
  CH - 交換相關(guān)能(XC)類型
Slater-XC
HL Hedin Lundquist (1971)
CA Ceperly and Alder parameterized by J.Perdew and Zunger
WI Wigner interpolation
PB Perdew -Becke
PW Perdew -Wang 86
LM Langreth-Mehl-Hu
91 Perdew -Wang 91
  QCOR - 非價(jià)鍵電子數(shù)（core electrons）
  第三行開始是每個(gè)軌道的具體參數(shù)，依次為n l j(=l±1/2)、原子軌道能、占有率

  PSCTR file of LDA/H1.25
TITEL = US H
LULTRA = T use ultrasoft PP ?
RWIGS = 0.57 nn distance ! Wigner-Seitz radius

RCLOC = .65
NE = 100
LCOR = .TRUE.
QCUT = -1
RMAX = 3.0 ! core radius for proj-oper

Description
l E TYP RCUT TYP RCUT(cutoff radius)
0 0 15 0.80 23 1.25
0 0.5 15 0.80 23 1.25
1 -0.2 15 0.80 23 1.25
  最重要的地方上面已經(jīng)用顏色標(biāo)出來啦：）
  說明一下，TYP是指贗勢(shì)的類型,RCUT是半徑截距，TYP可取的值如下：
  正則
  1 BHS
  2 TM
  3 VAN
  6 XNC
  7 RRKJ wave function possibly with node
  15 RRKJ wave function strictly no node
  非正則 +8

  最后一個(gè)問題是LDA or GGA。貌似沒有定論目前。
  這個(gè)最好是兩個(gè)一起做做看啦�；蛘呖次墨I(xiàn)別人驗(yàn)證過哪個(gè)數(shù)據(jù)好。
  其實(shí)據(jù)說目前最好的是PAW(P.E.Blochl,Phys.Rev.B 50,17953(1994).,Phys.Rev.B 59,1758(1999))。

原文地址 http://new.quantumchemistry.net/ ... dPAW/200512/99.html

希望大家能在自己的理解基礎(chǔ)上，說說自己的理解，好好交流交流。

[ Last edited by Prettyswan on 2010-11-17 at 18:22 ]

回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

» 猜你喜歡

調(diào)劑已經(jīng)有6人回復(fù)
280求調(diào)劑已經(jīng)有6人回復(fù)
302求調(diào)劑一志愿華中師范大學(xué) 已經(jīng)有5人回復(fù)
346求調(diào)劑已經(jīng)有4人回復(fù)
材料專碩322分已經(jīng)有9人回復(fù)
372分材料與化工（085600）一志愿湖南大學(xué)求調(diào)劑已經(jīng)有5人回復(fù)
318求調(diào)劑，計(jì)算材料方向已經(jīng)有10人回復(fù)
310求調(diào)劑已經(jīng)有14人回復(fù)
315求調(diào)劑已經(jīng)有3人回復(fù)
265求調(diào)劑已經(jīng)有13人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

twcplot 導(dǎo)入cas，dat文件顯示“out of memory” 已經(jīng)有10人回復(fù)
大家來幫我看看這兩個(gè)INCAR 文件有什么區(qū)別哦。。為啥計(jì)算結(jié)果不一樣呢已經(jīng)有4人回復(fù)
VASP中如何獲得POTCAR文件？已經(jīng)有8人回復(fù)
怎樣把CONTCAR拷貝成下一次計(jì)算的POTCAR文件再進(jìn)行計(jì)算已經(jīng)有6人回復(fù)
【求助】POTCAR庫中有5個(gè)函數(shù)，分別代表什么意思？已經(jīng)有3人回復(fù)
【求助成功】如何合并POTCAR文件已經(jīng)有8人回復(fù)
【求助】怎么建立POTCAR，vasp中找不到贗勢(shì)文件已經(jīng)有4人回復(fù)
【原創(chuàng)】CHGCAR轉(zhuǎn)化成MS可讀文件修正已經(jīng)有18人回復(fù)
【求助】vasp能夠提供同位素的POTCAR嗎？已經(jīng)有8人回復(fù)
【求助】關(guān)于POTCAR中EATOM 已經(jīng)有5人回復(fù)

1樓 2010-11-17 18:12:21

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wuli8

榮譽(yù)版主 (知名作家)

…………

1ST強(qiáng)帖: 2
應(yīng)助: 35 (小學(xué)生)
貴賓: 12.924
金幣: 20190.4
散金: 15888
紅花: 88
沙發(fā): 4
帖子: 7840
在線: 1114.6小時(shí)
蟲號(hào): 465889
注冊(cè): 2007-11-23
專業(yè): 物理學(xué)I
管轄: 計(jì)算模擬

很好

回復(fù)此樓

…………

2樓2010-11-17 20:15:34

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

sunyang1988

至尊木蟲 (著名寫手)

1ST強(qiáng)帖: 1
應(yīng)助: 24 (小學(xué)生)
貴賓: 0.984
金幣: 9397.7
散金: 1593
紅花: 23
沙發(fā): 1
帖子: 1744
在線: 1255.8小時(shí)
蟲號(hào): 965316
注冊(cè): 2010-03-09
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

奶昔啊奶昔！

回復(fù)此樓

心之所向無懼無悔

3樓2010-11-17 21:16:29

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zzy870720z

榮譽(yù)版主 (文壇精英)

1ST強(qiáng)帖: 3
應(yīng)助: 47 (小學(xué)生)
貴賓: 9.05
金幣: 30914.3
散金: 5613
紅花: 68
沙發(fā): 99
帖子: 12592
在線: 23567.6小時(shí)
蟲號(hào): 745488
注冊(cè): 2009-04-10
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
管轄: 分子模擬

支持一下，占樓學(xué)習(xí)

回復(fù)此樓

博學(xué)、審問、慎思、明辨、篤學(xué)

4樓2010-11-17 21:28:54

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

liushaking

鐵桿木蟲 (正式寫手)

應(yīng)助: 20 (小學(xué)生)
金幣: 5283.4
散金: 222
紅花: 15
帖子: 588
在線: 1160.8小時(shí)
蟲號(hào): 839616
注冊(cè): 2009-09-03
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

不錯(cuò)

回復(fù)此樓

勤快點(diǎn)就好了

5樓2011-03-05 10:22:22

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

bingmou

金蟲 (著名寫手)

1ST強(qiáng)帖: 5
應(yīng)助: 2 (幼兒園)
貴賓: 0.446
金幣: 431.9
散金: 171
紅花: 17
帖子: 1207
在線: 114.8小時(shí)
蟲號(hào): 487243
注冊(cè): 2007-12-28
專業(yè): 材料科學(xué)/半導(dǎo)體
管轄: 第一性原理

★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流

其實(shí)這篇文章里面說的在手冊(cè)里面都有

贊一下(1人)

回復(fù)此樓

6樓2011-03-05 17:20:28

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 071000生物學(xué)調(diào)劑 +6	知昭蔓 2026-04-02	6/300	2026-04-02 15:28 by sophie2180
[考研] 化學(xué)308分調(diào)劑 +14	你好明天你好 2026-03-30	15/750	2026-04-02 13:39 by 嘻嘻笑22
[考研] 275求調(diào)劑 +13	jjjjjjjjjjl 2026-03-27	13/650	2026-04-02 13:07 by yulian1987
[考研] 材料專業(yè)求調(diào)劑 +10	月月鳥木 2026-04-01	10/500	2026-04-02 12:57 by wxiongid
[考研] 找調(diào)劑 +6	楚喬喬 2026-04-01	6/300	2026-04-02 12:06 by xiaoranmu
[考研] 生醫(yī)工0831調(diào)劑求推薦 +5	小熊睿睿_s 2026-03-27	7/350	2026-04-02 11:30 by 紡大楊老師
[考研] 生物學(xué)308分求調(diào)劑（一志愿華東師大） +5	相信必會(huì)光芒萬?/a> 2026-03-31	6/300	2026-04-02 10:16 by guoweigw
[考研] 各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩 +13	Koxui 2026-03-28	13/650	2026-04-02 09:35 by 噠噠噠呱呱呱
[考研] 土木304求調(diào)劑 +6	兔突突突， 2026-03-31	7/350	2026-04-02 09:06 by coolminer
[考研] 安全工程 285 求調(diào)劑 +3	Xinyu56 2026-04-01	4/200	2026-04-01 21:50 by 靜靜靜靜靜靜靜?/a>
[考研] 085600 一志愿9 總分351 求調(diào)劑學(xué)校 +7	czhcz 2026-03-31	9/450	2026-04-01 19:24 by 無際的草原
[考研] 材料科學(xué)與工程339求調(diào)劑 +11	hyz0119 2026-03-31	12/600	2026-04-01 18:40 by 偉大河北
[考研] 0856調(diào)劑 +7	曲聽筠 2026-03-30	7/350	2026-04-01 08:51 by JourneyLucky
[考研] 考研調(diào)劑求助 +7	13287130938 2026-03-31	7/350	2026-03-31 16:39 by 690616278
[考研] 英一數(shù)一總分334求調(diào)劑 +4	陳陽坤 2026-03-31	4/200	2026-03-31 14:22 by 記事本2026
[考研] 085600，專業(yè)課化工原理，321分求調(diào)劑 +5	大饞小子 2026-03-28	5/250	2026-03-29 08:56 by qingfeng258
[考研] 藥學(xué)105500求調(diào)劑 +3	Ssun。。 2026-03-28	3/150	2026-03-28 11:24 by lxf170613
[考研] 340求調(diào)劑 +5	jhx777 2026-03-27	5/250	2026-03-28 04:18 by fmesaito
[考研] 復(fù)試調(diào)劑，一志愿南農(nóng)083200食品科學(xué)與工程 +5	XQTJZ 2026-03-26	5/250	2026-03-27 14:49 by 狂炫麥當(dāng)當(dāng)
[考研] 一志愿吉大071010，316分求調(diào)劑 +3	xgbiknn 2026-03-27	3/150	2026-03-27 10:36 by guoweigw