版塊導(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í)熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

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的回帖

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)文件用自動(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的回帖

查看全部 20 個(gè)回答

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的回帖

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)文件用自動(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的回帖

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í)候，是會(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的回帖

查看全部 20 個(gè)回答

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

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

[考研] 材料考研調(diào)劑 +5	Gs大王 2026-04-02	5/250	2026-04-02 07:02 by yjolah
[考研] 318求調(diào)劑，計(jì)算材料方向 +7	吸喵有害笙命 2026-04-01	8/400	2026-04-01 23:54 by 義文wang
[考研] 279求調(diào)劑 +3	學(xué)而思兮知 2026-04-01	3/150	2026-04-01 23:28 by 欣喜777
[考研] 304求調(diào)劑 +12	素年祭語 2026-03-31	15/750	2026-04-01 22:41 by peike
[考研] 270求調(diào)劑 +7	小杰pp 2026-03-31	8/400	2026-04-01 21:36 by 七度不信任
[考研] 085601材料工程找調(diào)劑 +20	oatmealR 2026-03-29	21/1050	2026-04-01 21:00 by lijunpoly
[考研] 0856初試324分求調(diào)劑 +5	想上學(xué)求調(diào) 2026-04-01	5/250	2026-04-01 15:28 by 考研上岸呀一定?/a>
[考研] 0805求調(diào)劑 +6	是水分 2026-03-31	6/300	2026-04-01 09:05 by oooqiao
[考研] 299求調(diào)劑 +8	嗯嗯嗯嗯2 2026-03-27	8/400	2026-03-31 18:23 by lizhi8172
[考研] 340求調(diào)劑 +4	希望如此i 2026-03-31	4/200	2026-03-31 16:40 by 690616278
[基金申請] 面上5B能上會(huì)嗎？ +8	redcom 2026-03-29	8/400	2026-03-31 15:53 by niuailing
[考研] 262求調(diào)劑 +7	ZZ..000 2026-03-30	8/400	2026-03-31 10:05 by cal0306
[考研] 一志愿中海洋320化學(xué)工程與技術(shù)學(xué)碩求調(diào)劑 +8	披星河 2026-03-30	8/400	2026-03-31 08:53 by lbsjt
[考研] 085600 286分材料求調(diào)劑 +11	麻辣魷魚 2026-03-27	12/600	2026-03-30 19:33 by Wang200018
[考研] 328求調(diào)劑 +8	嗯滴的基本都 2026-03-27	8/400	2026-03-30 17:20 by Wang200018
[考研] 303求調(diào)劑 +7	DLkz1314. 2026-03-30	7/350	2026-03-30 16:05 by shuang5186
[考研] 284求調(diào)劑 +14	junqihahaha 2026-03-26	15/750	2026-03-30 14:12 by 探123
[考研] 356求調(diào)劑 +4	gysy?s?a 2026-03-28	4/200	2026-03-29 10:32 by 唐沐兒
[考研] 308求調(diào)劑 +7	墨墨漠 2026-03-27	7/350	2026-03-28 07:43 by 熱情沙漠
[考研] 315分求調(diào)劑 +7	26考研上岸版26 2026-03-26	7/350	2026-03-28 04:05 by fmesaito

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

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

goodluck1982

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

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

» 猜你喜歡

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

goodluck1982

goodluck1982

bjwang

highpeter

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

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