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

登錄注冊

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

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

返回列表

goodluck1982

鐵蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 821.3
紅花: 1
帖子: 42
在線: 16.6小時(shí)
蟲號: 444269
注冊: 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勢，然后算能帶，然后根據(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í)
蟲號: 250077
注冊: 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)文件用自動的不行嗎？

我一般都是計(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í)
蟲號: 1599165
注冊: 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)。自動產(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。
按理說程序會判斷，我的KPOINTS中，有哪些個(gè)k點(diǎn)對應(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)對應(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í)
蟲號: 1599165
注冊: 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線段之后，我們再開始計(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í)
蟲號: 1599165
注冊: 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ì)算中依然自動用了自洽中的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í)
蟲號: 444269
注冊: 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í)
蟲號: 444269
注冊: 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)文件用自動的不行嗎？

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

要算指定的能帶和k點(diǎn)。自動產(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í)
蟲號: 898205
注冊: 2009-11-09
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

引用回帖:

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

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

贊一下

回復(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í)
蟲號: 898205
注冊: 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í)
蟲號: 898205
注冊: 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 的主題更新

返回列表

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

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

[考研] 0856材料專碩353求調(diào)劑 +3	NIFFFfff 2026-03-20	3/150	2026-03-21 10:23 by luoyongfeng
[考研] 材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué) +4	順順順mr 2026-03-18	5/250	2026-03-21 10:22 by luoyongfeng
[考研] 310求調(diào)劑 +3	baibai1314 2026-03-16	3/150	2026-03-21 03:56 by JourneyLucky
[考研] 265求調(diào)劑 +3	Jack?k?y 2026-03-17	3/150	2026-03-21 03:17 by JourneyLucky
[考研] 初始318分求調(diào)劑（有工作經(jīng)驗(yàn)） +3	1911236844 2026-03-17	3/150	2026-03-21 02:33 by JourneyLucky
[考研] 化學(xué)求調(diào)劑 +4	臨澤境llllll 2026-03-17	5/250	2026-03-21 02:23 by JourneyLucky
[考研] 265求調(diào)劑 +9	梁梁校校 2026-03-17	9/450	2026-03-21 02:17 by JourneyLucky
[考研] 材料 336 求調(diào)劑 +3	An@. 2026-03-18	4/200	2026-03-21 01:39 by JourneyLucky
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4	葵梓衛(wèi)隊(duì) 2026-03-18	6/300	2026-03-20 23:02 by JourneyLucky
[考研] 一志愿中南化學(xué)（0703）總分337求調(diào)劑 +8	niko- 2026-03-19	9/450	2026-03-20 21:57 by luoyongfeng
[考研] 求調(diào)劑 +3	eation27 2026-03-20	3/150	2026-03-20 19:32 by JourneyLucky
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +9	手機(jī)用戶 2026-03-17	9/450	2026-03-20 14:24 by 無懈可擊111
[考研] 一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5	1孫悟空 2026-03-17	6/300	2026-03-19 23:46 by zcl123
[考研] 081700化工學(xué)碩調(diào)劑 +3	【1】 2026-03-16	3/150	2026-03-19 23:40 by edmund7
[考研] 梁成偉老師課題組歡迎你的加入 +9	一鴨鴨喲 2026-03-14	11/550	2026-03-19 17:22 by ！本暗一次！
[考研] 一志愿福大288有機(jī)化學(xué)，求調(diào)劑 +3	小木蟲200408204 2026-03-18	3/150	2026-03-19 13:31 by houyaoxu
[考研] 311求調(diào)劑 +11	冬十三 2026-03-15	12/600	2026-03-18 14:36 by 星空星月
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3	麥田耕者 2026-03-14	3/150	2026-03-14 20:04 by 外星文明

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

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

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

» 猜你喜歡

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

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