版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

goodluck1982

鐵蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 821.3
紅花: 1
帖子: 42
在線: 16.6小時(shí)
蟲號(hào): 444269
注冊(cè): 2007-10-27
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[交流] vasp 計(jì)算 partial charge 時(shí)k點(diǎn)選取的問題已有9人參與

比如 graphene， SCF 計(jì)算時(shí)使用 PAW PBE勢(shì)，然后算能帶，然后根據(jù)能帶的波函數(shù)提取出M點(diǎn)價(jià)帶的 partial charge。然而使用兩種 k-line 取法時(shí)得到的結(jié)果卻很不一樣：
(1) band1: 第一種就是直接算 Gamma-K-M-Gamma 一個(gè)比較完整的能帶，每段41個(gè)點(diǎn)，那么 M 點(diǎn)是第 82和83 個(gè)k點(diǎn)。
(2) band2: 第二種就是直接算 M 到 -M，就取兩個(gè)點(diǎn)，那么 M 點(diǎn)就是第一個(gè) k 點(diǎn)。

然后 (1) 得到的 PARCHG.0004.0082 和（2）得到的 PARCHG.0004.0001 差別很大。（可以用 diff 或 vimdiff 比較）
這是怎么回事呢？該如何選取K點(diǎn)好呢？
而且即使在（1）中， PARCHG.0004.0082 和 PARCHG.0004.0083 也是有些許差別的，但這個(gè)可以看出是數(shù)值誤差。

附： scf 的 INCAR 如下

CODE:

system = graphene

ISTART=1

ICHARG=1

PREC=High

ENCUT=400

EDIFF=1e-6

ISMEAR=-5

KPOINTS 取 40x40x1
POSCAR 為

CODE:

graphene

2.46

1         0            0

0.5       0.8660254    0

0         0            4

C

2

Cartesian

0         0          0         

0.5       0.2886751  0

band1 INCAR:

CODE:

ISTART=1

ICHARG=11

PREC=High

ENCUT=400

EDIFF=1e-6

ISMEAR=1

LWAVE  = T

之后計(jì)算 partial charge 的 INCAR 為

CODE:

ISTART=1

ICHARG=1

PREC=High

ENCUT=400

EDIFF=1e-6

LPARD=T

IBAND=4

KPUSE=82 83 

LSEPB=T

LSEPK=T

KPOINTS為

CODE:

for band: Gamma-K-M-Gamma

41

line mode

Reciprocal

0           0           0   ! Gamma

0.6666667   0.3333333   0   ! K

0.6666667   0.3333333   0   ! K

0.5         0           0   ! M

0.5         0           0   ! M

0           0           0   ! Gamma

band2 的 INCAR 只改

CODE:

KPUSE=1   2

KPOINTS 相應(yīng)為

CODE:

for band

2

line mode

Reciprocal

0.5         0           0   ! M

-0.5         0           0   !- M

[ Last edited by goodluck1982 on 2013-1-31 at 18:00 ]

計(jì)算 partial charge 的 INCAR 中的 ICHARG 應(yīng)為 11

[ Last edited by goodluck1982 on 2015-3-19 at 17:44 ]

回復(fù)此樓

1樓 2013-01-31 17:54:15

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

回帖支持 ( 顯示支持度最高的前 50 名 )

bjwang

鐵桿木蟲 (正式寫手)

應(yīng)助: 4 (幼兒園)
金幣: 5860.3
散金: 94
紅花: 6
帖子: 792
在線: 556.3小時(shí)
蟲號(hào): 250077
注冊(cè): 2006-05-08
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

進(jìn)行Partial Charge分析必須要保留有自洽計(jì)算的WAVECAR才可以。http://www.gaoyang168.com/bbs/viewthread.php?tid=1404880
至于k點(diǎn)文件用自動(dòng)的不行嗎？

我一般都是計(jì)算某個(gè)能量區(qū)間的，采用和自洽相同的k點(diǎn)文件。沒計(jì)算過指定能帶的，指定計(jì)算能帶的k點(diǎn)文件要和計(jì)算能帶的k點(diǎn)文件一樣？

[ Last edited by bjwang on 2013-1-31 at 22:07 ]

贊一下(1人)

回復(fù)此樓

It is a shame for the soul to be first to give way in this life,when thy body does not give way.

3樓2013-01-31 21:31:23

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

北漂2012

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 94.9
紅花: 1
帖子: 29
在線: 95.8小時(shí)
蟲號(hào): 1599165
注冊(cè): 2012-02-04
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

引用回帖:

4樓: Originally posted by goodluck1982 at 2013-02-02 17:08:37
要算指定的能帶和k點(diǎn)。自動(dòng)產(chǎn)生的k點(diǎn)里不一定有我想要的k點(diǎn)，那只能通過能帶方式指定了。...

我先自洽，然后畫出能帶。根據(jù)能帶上，我要計(jì)算fermi面附近一定能量區(qū)間的電荷密度。
然后做parital charge的計(jì)算，當(dāng)然用的是自洽的WAVRCAR，選擇費(fèi)米面附近-0.1eV-0.1eV的能量區(qū)間
其中INCAR的設(shè)置：
NBMOD=-3
EINT=-0.1 0.1
LPARD=.TRUE.
還有這個(gè)時(shí)候我用的KPOINTS，是畫能帶的KPOINTS。
按理說程序會(huì)判斷，我的KPOINTS中，有哪些個(gè)k點(diǎn)對(duì)應(yīng)的能量落在我指定的能量區(qū)間內(nèi)。
(1)我發(fā)現(xiàn)，程序判斷的k點(diǎn)和我從能帶圖上直接觀察到屬于指定能量區(qū)間的k點(diǎn)，不一致。比如從能帶上看，我看到的第5條能帶的第10個(gè)k點(diǎn)落在我選擇的能量區(qū)間；而程序自己找的卻是是第5條能帶的第2個(gè)k點(diǎn)。這樣的話，得到的patial電荷密度結(jié)果就應(yīng)該不可靠了。
(2)前面計(jì)算畫能帶得到的本征值文件EIGENVAL，和做partial charge charge得到的EIGENVAL文件，雖然他們的KPOINTS是一樣的,但是里面同一個(gè)k點(diǎn)對(duì)應(yīng)的能量值，卻是不同的。
我覺得是兩個(gè)過程的EIGENVAL的不同，導(dǎo)致問題(1)

贊一下(1人)

回復(fù)此樓

ibelieveicanfly,ibelieveeicantouchthesky.Ificanseeit,thenicandoit,ifijustbelieveit,thethereisnothingtoit.

7樓2014-02-23 12:47:59

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

北漂2012

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 94.9
紅花: 1
帖子: 29
在線: 95.8小時(shí)
蟲號(hào): 1599165
注冊(cè): 2012-02-04
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

引用回帖:

6樓: Originally posted by highpeter at 2013-10-07 21:38:41
弱弱的問一下，為什么82.83是同一個(gè)k點(diǎn)��？...

K-points
82
Line-mode
0 0 0    !G
0 0 0.5 !Z

0 0 0.5  !Z
0.5 0 0  ! X
比如上面是我畫能帶的KPOINTS文件。
把G-Z點(diǎn)之間線段，分成81(82-1)份，包括線段兩端，共有82個(gè)點(diǎn)。G是第一個(gè)點(diǎn)，Z是第82個(gè)點(diǎn)；
同理
把Z-X點(diǎn)之間線段，也分成81(82-1)份；包括線段兩端，共有82個(gè)點(diǎn)。接著G-Z線段之后，我們?cè)匍_始計(jì)數(shù)，這里Z是第83(82+1)點(diǎn)，X是第164(82+82)個(gè)點(diǎn)；
這樣，第82和83個(gè)點(diǎn)都是Z點(diǎn)。不信你可以看一下畫出來的能帶數(shù)據(jù)。

贊一下(1人)

回復(fù)此樓

ibelieveicanfly,ibelieveeicantouchthesky.Ificanseeit,thenicandoit,ifijustbelieveit,thethereisnothingtoit.

8樓2014-02-23 12:59:47

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

北漂2012

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 94.9
紅花: 1
帖子: 29
在線: 95.8小時(shí)
蟲號(hào): 1599165
注冊(cè): 2012-02-04
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

引用回帖:

同學(xué)，在嗎。
VASP 做partial charge 計(jì)算時(shí)候。
可以通過畫能帶的KPOINTS，去指定某條能帶，某個(gè)k點(diǎn)嗎?
為什么我覺得我用KPOINTS去指定能帶，好像不起作用。
它實(shí)際計(jì)算中依然自動(dòng)用了自洽中的K點(diǎn)呢？
謝謝啦！

贊一下(1人)

回復(fù)此樓

ibelieveicanfly,ibelieveeicantouchthesky.Ificanseeit,thenicandoit,ifijustbelieveit,thethereisnothingtoit.

9樓2014-02-24 22:45:18

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

普通回帖

goodluck1982

鐵蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 821.3
紅花: 1
帖子: 42
在線: 16.6小時(shí)
蟲號(hào): 444269
注冊(cè): 2007-10-27
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

可能是算能帶時(shí) CHGCAR 文件用亂了。
后來重復(fù)了一遍，算能帶時(shí)使用相同的 scf 收斂后的CHGCAR，最后后 (1) 得到的 PARCHG.0004.0082 和（2）得到的 PARCHG.0004.0001 差別不大了，不過還是有數(shù)值誤差，在小數(shù)點(diǎn)后四五位左右，比 PARCHG.0004.0082 和 PARCHG.0004.0083 間的差別稍大。

不過，82 和 83 明明是同一個(gè)k點(diǎn)，算出來卻不完全一樣，讓人覺得不爽。

贊一下

回復(fù)此樓

2樓2013-01-31 18:11:12

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

goodluck1982

鐵蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 821.3
紅花: 1
帖子: 42
在線: 16.6小時(shí)
蟲號(hào): 444269
注冊(cè): 2007-10-27
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

3樓: Originally posted by bjwang at 2013-01-31 21:31:23
進(jìn)行Partial Charge分析必須要保留有自洽計(jì)算的WAVECAR才可以。http://www.gaoyang168.com/bbs/viewthread.php?tid=1404880
至于k點(diǎn)文件用自動(dòng)的不行嗎？

我一般都是計(jì)算某個(gè)能量區(qū)間的，采用和自洽相同的k點(diǎn)文件。沒 ...

要算指定的能帶和k點(diǎn)。自動(dòng)產(chǎn)生的k點(diǎn)里不一定有我想要的k點(diǎn)，那只能通過能帶方式指定了。

贊一下

回復(fù)此樓

4樓2013-02-02 17:08:37

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

highpeter

金蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 952.1
帖子: 43
在線: 113小時(shí)
蟲號(hào): 898205
注冊(cè): 2009-11-09
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

引用回帖:

我一般都是計(jì)算某個(gè)能量區(qū)間的，采用和自洽相同的k點(diǎn)文件。沒 ...

計(jì)算能量區(qū)間的時(shí)候，是會(huì)自動(dòng)幫我們選定能帶的吧。但是用這種方法算費(fèi)米能附近partial charge很麻煩，有些時(shí)候系統(tǒng)自動(dòng)給我們選的費(fèi)米能會(huì)非常夸張的大。
如果在INCAR的EINT設(shè)置費(fèi)米能，又有計(jì)算partial charge和自洽計(jì)算時(shí)費(fèi)米能會(huì)有一定浮動(dòng)的問題。

贊一下

回復(fù)此樓

5樓2013-10-07 21:27:00

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

highpeter

金蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 952.1
帖子: 43
在線: 113小時(shí)
蟲號(hào): 898205
注冊(cè): 2009-11-09
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

引用回帖:

2樓: Originally posted by goodluck1982 at 2013-01-31 18:11:12
可能是算能帶時(shí) CHGCAR 文件用亂了。
后來重復(fù)了一遍，算能帶時(shí)使用相同的 scf 收斂后的CHGCAR，最后后 (1) 得到的 PARCHG.0004.0082 和（2）得到的 PARCHG.0004.0001 差別不大了，不過還是有數(shù)值誤差，在小數(shù)點(diǎn) ...

弱弱的問一下，為什么82.83是同一個(gè)k點(diǎn)��？

贊一下

回復(fù)此樓

6樓2013-10-07 21:38:41

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

highpeter

金蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 952.1
帖子: 43
在線: 113小時(shí)
蟲號(hào): 898205
注冊(cè): 2009-11-09
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

引用回帖:

7樓: Originally posted by 北漂2012 at 2014-02-23 12:47:59
我先自洽，然后畫出能帶。根據(jù)能帶上，我要計(jì)算fermi面附近一定能量區(qū)間的電荷密度。
然后做parital charge的計(jì)算，當(dāng)然用的是自洽的WAVRCAR，選擇費(fèi)米面附近-0.1eV-0.1eV的能量區(qū)間
其中INCAR的設(shè)置：
NBMOD=- ...

你選EINT=-0.1 0.1的話，它應(yīng)該是能量-0.1到0.1的范圍，而不是費(fèi)米能附近的。

贊一下

回復(fù)此樓

10樓2014-02-26 12:27:00

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 goodluck1982 的主題更新

返回列表

普通表情龍兔虎貓高級(jí)回復(fù) (可上傳附件)

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

[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +12	yangfz 2026-03-17	12/600	2026-03-21 03:30 by JourneyLucky
[考研] 303求調(diào)劑 +5	睿08 2026-03-17	7/350	2026-03-21 03:11 by JourneyLucky
[考研] 295求調(diào)劑 +4	一志愿京區(qū)211 2026-03-18	6/300	2026-03-20 23:41 by JourneyLucky
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +9	Ncdx123456 2026-03-19	9/450	2026-03-20 23:41 by lovewei0727
[考研] 304求調(diào)劑 +7	司空. 2026-03-18	7/350	2026-03-20 23:08 by JourneyLucky
[考研] 考研調(diào)劑求學(xué)校推薦 +3	伯樂29 2026-03-18	5/250	2026-03-20 22:59 by JourneyLucky
[考研] 085600材料與化工 +8	安全上岸！ 2026-03-16	8/400	2026-03-20 22:13 by luoyongfeng
[考研] 290求調(diào)劑 +7	^O^乜 2026-03-19	7/350	2026-03-20 21:43 by JourneyLucky
[考研] 261求B區(qū)調(diào)劑，科研經(jīng)歷豐富 +3	牛奶很忙 2026-03-20	4/200	2026-03-20 19:34 by JourneyLucky
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +9	手機(jī)用戶 2026-03-17	9/450	2026-03-20 14:24 by 無懈可擊111
[考博] 招收博士1-2人 +3	QGZDSYS 2026-03-18	3/150	2026-03-20 11:58 by 呱呱呱呱叫
[考博] 申博26年 +3	八6八68 2026-03-19	3/150	2026-03-19 19:43 by nxgogo
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5	ioodiiij 2026-03-17	5/250	2026-03-19 18:22 by zcl123
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518
[考研] 化學(xué)工程321分求調(diào)劑 +15	大米飯！ 2026-03-15	18/900	2026-03-18 14:52 by haxia
[考研] 0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國級(jí)證書40+數(shù)一英一287 +8	小李0854 2026-03-16	8/400	2026-03-18 14:35 by 搏擊518
[考研] 304求調(diào)劑 +12	小熊joy 2026-03-14	13/650	2026-03-18 12:34 by Linda Hu
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[考研] 290求調(diào)劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm
[考研] 283求調(diào)劑 +3	聽風(fēng)就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠

24小時(shí)熱門版塊排行榜

[交流] vasp 計(jì)算 partial charge 時(shí)k點(diǎn)選取的問題 已有9人參與

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

» 猜你喜歡

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

[交流] vasp 計(jì)算 partial charge 時(shí)k點(diǎn)選取的問題已有9人參與

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