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

返回列表

1984mjx

應(yīng)助: (幼兒園)
在線(xiàn):
蟲(chóng)號(hào): 629105

[交流] 【求助】PBE0 計(jì)算Si、Ge能帶結(jié)構(gòu) 已有6人參與

請(qǐng)教高手，用vasp5.2的pbe0計(jì)算硅鍺能帶結(jié)構(gòu)時(shí)，所得的能帶結(jié)構(gòu)圖上總會(huì)出現(xiàn)尖銳的峰值，顯然計(jì)算數(shù)據(jù)有問(wèn)題，用4.6版計(jì)算的能帶結(jié)構(gòu)和教科書(shū)上比較溫和，為何同樣的設(shè)置做PBE0時(shí)就出錯(cuò)了。請(qǐng)做過(guò)PBE0計(jì)算的大蝦指點(diǎn)一下

以下是我的個(gè)文件：
以ge為例,做優(yōu)化時(shí)的設(shè)置
INCAR:
SYSTEM =  Azoc MD

   PREC = Medium
   EDIFF = 1e-5
   EDIFFG  = -2e-2
   IALGO  =  48

   NPAR    = 2
   LPLANE = TRUE
   LCHARG  = .TRUE.

   NELMIN  = 3
   LREAL = F
   ROPT = 1E-3
   ISYM = 1
   NWRITE  = 1

  SPIN POLARIZATION SPECIFICATION
   ISPIN = 1

  IONIC RELAXATION
   NSW    = 200
   NBLOCK  = 1
   IBRION  = 2
   POTIM = 0.5
   LWAVE = .F.
   LCHG  = .T.
   ISIF = 2
   LCHARG  = .TRUE.
  ELECTRON OCCUPATION
   ISMEAR = 0 occupation forced by INCAR

  DOS Stuff
   SIGMA = 0.1    check T*S to make sure
   LHFCALC = .TRUE.
   AEXX = 0.2

KPOINTS:
0
Monkhorst-Pack
4 4 4
0.  0.  0.
POSCAR:
  Ge
5.6575
0.5 0.5 0.0
0.0 0.5 0.5
0.5 0.0 0.5
2
Cartesian
   0.0000000000000    0.00000000000    0.0000000000000
   0.2500000000000    0.25000000000    0.2500000000000

算能帶時(shí)設(shè)置
INCAR:
SYSTEM =  Azoc MD

   PREC = Medium
   EDIFF = 1e-5
   EDIFFG  = -2e-2
   IALGO  =  48

   NPAR    = 2
   LPLANE = TRUE
   LCHARG  = .TRUE.

   NELMIN  = 3
   LREAL = F
   ROPT = 1E-3
   ISYM = 1
   NWRITE  = 1
ICHARG = 11

  SPIN POLARIZATION SPECIFICATION
   ISPIN = 1

  IONIC RELAXATION
   NSW    = 1
   NBLOCK  = 1
   IBRION  = 1
   POTIM = 0.5
   LWAVE = .F.
   LCHG  = .T.
   ISIF = 2
   LCHARG  = .TRUE.
  ELECTRON OCCUPATION
   ISMEAR = 0 occupation forced by INCAR

  DOS Stuff
   SIGMA = 0.1    check T*S to make sure
   LHFCALC = .TRUE.
   AEXX = 0.2

KPOINTS:
Automatic mesh
19
rec
0.50000000  0.50000000  0.50000000    0.050
0.44444444  0.44444444  0.44444444    0.050
0.38888889  0.38888889  0.38888889    0.050
0.33333333  0.33333333  0.33333333    0.050
0.27777778  0.27777778  0.27777778    0.050
0.22222222  0.22222222  0.22222222    0.050
0.16666667  0.16666667  0.16666667    0.050
0.11111111  0.11111111  0.11111111    0.050
0.05555556  0.05555556  0.05555556    0.050
0.00000000  0.00000000  0.00000000    0.050
0.05555556  0.05555556  0.00000000    0.050
0.11111111  0.11111111  0.00000000    0.050
0.16666667  0.16666667  0.00000000    0.050
0.22222222  0.22222222  0.00000000    0.050
0.27777778  0.27777778  0.00000000    0.050
0.33333333  0.33333333  0.00000000    0.050
0.38888889  0.38888889  0.00000000    0.050
0.44444444  0.44444444  0.00000000    0.050
0.50000000  0.50000000  0.00000000    0.050

懇請(qǐng)高手指點(diǎn)，非常急，拜托了！

[ Last edited by 1984mjx on 2010-6-4 at 16:29 ]

回復(fù)此樓

1樓 2010-06-04 16:26:42

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

valenhou001

至尊木蟲(chóng) (職業(yè)作家)

1ST強(qiáng)帖: 13
應(yīng)助: 241 (大學(xué)生)
金幣: 25701.7
散金: 602
紅花: 166
帖子: 3782
在線(xiàn): 873.8小時(shí)
蟲(chóng)號(hào): 1007127
注冊(cè): 2010-04-27
專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★ ★ ★
1984mjx(金幣+1):謝謝參與
gavinliu7390(金幣+2):謝謝交流！ 2010-06-04 22:17:23

請(qǐng)參考vasp官方論壇提到的如何進(jìn)行hybrid functional來(lái)計(jì)算band structure的步驟：
下面的是從上面給拷貝下來(lái)的：
band structures using hybrid functionals have to be calculated the following way:
-) First perform a selfconsistent Hartree-Fock/HSE calculation using a conventional KPOINTS file.
-) Copy the IBZKPT file to KPOINTS, and explicitely add all desired k-points along high-symmetry lines of the BZ that are needed for the bandstructure plot. Add the points at the end of the KPOINTS file, but set the weights of these added k-points to 0.
-) Do not forget to set the number of k-points in KPOINTS correctly (to the number of the k-points used in the standard mesh PLUS the number of the k-points along the lines), such that all k-points are used for the calculations.
-) NKRED can not be used.
-) Perform a second VASP run:
It is recommended to use the Davidson algorithm, since it converges that eigen energies at the new k-points fastest. Since VASP terminates when the total energy is converged to a certain threshold, it is important to force VASP to do a minimum number of steps, so that the orbitals at the new k-points are fully converged (note: since their weight is zero, they do not contribute to the total energy).
This can be done using e.g.

ALGO = N ; NELMIN = 5 ! Davidson, minimum 5 scf-steps
IBRION = 1 ! Use simple charge mixer, since
Pulay might blow up

-) The KS-eigenvalues of the states along the high-symmetry lines are written in OUTCAR, EIGENVAL, vasprun.xml; please cut the k-points required for the bandstructure from one of these files and proceed as usual (using p4vasp or any other graphics tool you usually use to produce bandstructure-plots)

Note: A Hartee-Fock calculation can NOT be continued from an existing CHGCAR file, since the non-local exchange is not determined by the charge density but by the density matrix and/or the KS-orbitals.

贊一下(2人)

回復(fù)此樓

2樓2010-06-04 19:13:04

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

vasp001

金蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 69 (初中生)
金幣: 1513.7
散金: 134
紅花: 24
帖子: 702
在線(xiàn): 326小時(shí)
蟲(chóng)號(hào): 727265
注冊(cè): 2009-03-20
性別: GG
專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★ ★
1984mjx(金幣+1):謝謝參與
gavinliu7390(金幣+1):謝謝交流！ 2010-06-04 22:17:31
1984mjx(金幣+1): 2010-06-06 17:15:28

引用回帖:

Originally posted by 1984mjx at 2010-06-04 16:26:42:
請(qǐng)教高手，用vasp5.2的pbe0計(jì)算硅鍺能帶結(jié)構(gòu)時(shí)，所得的能帶結(jié)構(gòu)圖上總會(huì)出現(xiàn)尖銳的峰值，顯然計(jì)算數(shù)據(jù)有問(wèn)題，用4.6版計(jì)算的能帶結(jié)構(gòu)和教科書(shū)上比較溫和，為何同樣的設(shè)置做PBE0時(shí)就出錯(cuò)了。請(qǐng)做過(guò)PBE0計(jì)算的大蝦指 ...

看到了樓上的帖子了，實(shí)際上官方論壇的VASP001 和現(xiàn)在的賬號(hào)是一個(gè)人，都是我提的。對(duì)于Si這個(gè)體系而言，你用HSE06 得出的帶隙和實(shí)驗(yàn)值僅差0.1ev，而且，計(jì)算時(shí)間比PBE0 少很多，

具體原因，自己看看，這個(gè)庫(kù)侖截?cái)噙€是起的作用很大。此外，不一定用AlGO=Normal ,實(shí)際上ALGO=Damped 或者ALL 更快，我測(cè)試的結(jié)果是DAMPED最快。

贊一下(2人)

回復(fù)此樓

一定要精通MaterialsStudios,VASP.

3樓2010-06-04 22:07:21

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

1984mjx

應(yīng)助: (幼兒園)
在線(xiàn):
蟲(chóng)號(hào): 629105

PBE0 計(jì)算Si、Ge能帶結(jié)構(gòu)

首先非常感謝樓上的同仁的幫助，我也是自己在摸索pbe0的計(jì)算，不太懂里面的具體設(shè)置，花了不少時(shí)間也沒(méi)解決好這個(gè)問(wèn)題，用PAE0做能帶計(jì)算的步驟和gga的不一樣嗎，我都是先做結(jié)構(gòu)優(yōu)化，然后將產(chǎn)生的CHGCAR/CHG文件考到新文件夾下，將ICHARG = 11設(shè)定后計(jì)算，最后在產(chǎn)生的文件EIGENVAL讀取能帶信息，請(qǐng)樓上的幫個(gè)忙，能把你做pBE0計(jì)算能帶的INCAR 、kpoints、poscar、貼出來(lái)共享一下嗎，不勝感激額，先行謝過(guò)了！

贊一下

回復(fù)此樓

4樓2010-06-06 17:14:41

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

vasp001

金蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 69 (初中生)
金幣: 1513.7
散金: 134
紅花: 24
帖子: 702
在線(xiàn): 326小時(shí)
蟲(chóng)號(hào): 727265
注冊(cè): 2009-03-20
性別: GG
專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★ ★
zzy870720z(金幣+2):感謝分享 2010-06-07 08:04:45
1984mjx(金幣+1):非常感謝！ 2010-06-08 09:20:46

引用回帖:

Originally posted by 1984mjx at 2010-06-06 17:14:41:
首先非常感謝樓上的同仁的幫助，我也是自己在摸索pbe0的計(jì)算，不太懂里面的具體設(shè)置，花了不少時(shí)間也沒(méi)解決好這個(gè)問(wèn)題，用PAE0做能帶計(jì)算的步驟和gga的不一樣嗎，我都是先做結(jié)構(gòu)優(yōu)化，然后將產(chǎn)生的CHGCAR/CHG文件 ...

計(jì)算能帶的方法見(jiàn)VASP 官方論壇，有別于常規(guī)的DFT 能帶計(jì)算。。你有問(wèn)題，再提問(wèn)吧。
Band like this :
> # Startparameter for this run:
> ISPIN = 2
> ISTART = 0
> ICHARG = 2
> NWRITE = 2
> PREC = Normal
>
> # Electronic Relaxation
> ENCUT = 400.0
> LREAL = AUTO
> EDIFF = 1E-04
> NELM = 80
> NELMIN = 5
>
> # Ionic relaxation
> EDIFFG = -1E-02
> NSW = 0
> POTIM = 0.50
> IBRION = 1
> ISIF = 2
> ISYM = 2
>
>
> # DOS related values:
> ISMEAR = 0; SIGMA = 0.20
>
>
> # Electronic relaxation 2 (details)
> ALGO = Normal
>
> # Write flags
> LWAVE =.TRUE.
> LCHARG =.FALSE.
> LORBIT = 11
> NBANDS=48
> # Accelorate computation
> NSIM = 4
> NPAR = 64
> LPLANE =.TRUE.
> LSCALU =.FALSE.
> MAXMIX = 80
>
> # Exchange correlation treatment:
> LHFCALC = .TRUE.
> HFSCREEN = 0.2
> AEXX = 0.25
> AGGAX = 0.75
> ALDAC = 1.0
> AGGAc = 1.0
>
>
> ENCUTFOCK= 0
> LMAXFOCK= 4
> # NKRED = 2
> TIME = 0.4
>
> KPOINTS is like this :
>
> Automatically generated mesh
> 56
> Reciprocal lattice
> 0.05000000000000 0.05000000000000 0.12500000000000 8 0.15000000000000
> 0.05000000000000 0.12500000000000 16 0.25000000000000 0.05000000000000
> 0.12500000000000 16 0.35000000000000 0.05000000000000 0.12500000000000
> 16 0.45000000000000 0.05000000000000 0.12500000000000 16
> 0.15000000000000 0.15000000000000 0.12500000000000 8 0.25000000000000
> 0.15000000000000 0.12500000000000 16 0.35000000000000 0.15000000000000
> 0.12500000000000 16 0.45000000000000 0.15000000000000 0.12500000000000
> 16 0.25000000000000 0.25000000000000 0.12500000000000 8
> 0.35000000000000 0.25000000000000 0.12500000000000 16 0.45000000000000
> 0.25000000000000 0.12500000000000 16 0.35000000000000 0.35000000000000
> 0.12500000000000 8 0.45000000000000 0.35000000000000 0.12500000000000
> 16 0.45000000000000 0.45000000000000 0.12500000000000 8
> 0.05000000000000 0.05000000000000 0.37500000000000 8 0.15000000000000
> 0.05000000000000 0.37500000000000 16 0.25000000000000 0.05000000000000
> 0.37500000000000 16 0.35000000000000 0.05000000000000 0.37500000000000
> 16 0.45000000000000 0.05000000000000 0.37500000000000 16
> 0.15000000000000 0.15000000000000 0.37500000000000 8 0.25000000000000
> 0.15000000000000 0.37500000000000 16 0.35000000000000 0.15000000000000
> 0.37500000000000 16 0.45000000000000 0.15000000000000 0.37500000000000
> 16 0.25000000000000 0.25000000000000 0.37500000000000 8
> 0.35000000000000 0.25000000000000 0.37500000000000 16 0.45000000000000
> 0.25000000000000 0.37500000000000 16 0.35000000000000 0.35000000000000
> 0.37500000000000 8 0.45000000000000 0.35000000000000 0.37500000000000
> 16 0.45000000000000 0.45000000000000 0.37500000000000 8 0.0 0.0 0.0 0
> 0.1 0.1 0.0 0
> 0.2 0.2 0.0 0
> 0.3 0.3 0.0 0
> 0.4 0.4 0.0 0
> 0.5 0.5 0.0 0
> 0.5 0.5 0.1 0
> 0.5 0.5 0.2 0
> 0.5 0.5 0.3 0
> 0.5 0.5 0.4 0
> 0.5 0.5 0.5 0
> 0.5 0.4 0.5 0
> 0.5 0.3 0.5 0
> 0.5 0.2 0.5 0
> 0.5 0.1 0.5 0
> 0.5 0.0 0.5 0
> 0.4 0.0 0.5 0
> 0.3 0.0 0.5 0
> 0.2 0.0 0.5 0
> 0.1 0.0 0.5 0
> 0.0 0.0 0.5 0
> 0.0 0.0 0.4 0
> 0.0 0.0 0.3 0
> 0.0 0.0 0.2 0
> 0.0 0.0 0.1 0
> 0.0 0.0 0.0 0

這是HSE06 的，PBE0 ，要去掉那個(gè)HFSCREEN 截?cái)唷?br />
[ Last edited by vasp001 on 2010-6-9 at 14:28 ]

贊一下(1人)

回復(fù)此樓

一定要精通MaterialsStudios,VASP.

5樓2010-06-06 23:18:33

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

wo8165827

金蟲(chóng) (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1446.6
散金: 63
紅花: 4
帖子: 297
在線(xiàn): 90.6小時(shí)
蟲(chóng)號(hào): 612039
注冊(cè): 2008-09-25
性別: GG
專(zhuān)業(yè): 金屬材料的合金相、相變及

★ ★
1984mjx(金幣+1):謝謝參與
xiaohunhun(金幣+1):xiexie 2010-06-08 23:29:15
1984mjx(金幣+1): 2010-06-10 21:30:37

態(tài)密度計(jì)算過(guò)程中，你可以計(jì)算出來(lái)的效果不是很平滑，那就需要調(diào)節(jié)參數(shù)來(lái)計(jì)算，從結(jié)構(gòu)優(yōu)化開(kāi)始到DOS計(jì)算都是有可以調(diào)的參數(shù)。具體可以參考侯老師的說(shuō)明書(shū)，

贊一下(2人)

回復(fù)此樓

計(jì)算的魅力

6樓2010-06-08 17:55:34

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

1984mjx

應(yīng)助: (幼兒園)
在線(xiàn):
蟲(chóng)號(hào): 629105

引用回帖:

Originally posted by wo8165827 at 2010-06-08 17:55:34:
態(tài)密度計(jì)算過(guò)程中，你可以計(jì)算出來(lái)的效果不是很平滑，那就需要調(diào)節(jié)參數(shù)來(lái)計(jì)算，從結(jié)構(gòu)優(yōu)化開(kāi)始到DOS計(jì)算都是有可以調(diào)的參數(shù)。具體可以參考侯老師的說(shuō)明書(shū)，

非常感謝您的指點(diǎn)，再麻煩您一下，能給我發(fā)份侯老師的說(shuō)明書(shū)嗎，再次感謝郵箱：1984mjx@163.com

贊一下

回復(fù)此樓

7樓2010-06-10 21:33:32

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

wo8165827

金蟲(chóng) (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1446.6
散金: 63
紅花: 4
帖子: 297
在線(xiàn): 90.6小時(shí)
蟲(chóng)號(hào): 612039
注冊(cè): 2008-09-25
性別: GG
專(zhuān)業(yè): 金屬材料的合金相、相變及

好的。不過(guò)壇子里搜下就有了。。。我發(fā)給你一份吧。。

贊一下

回復(fù)此樓

計(jì)算的魅力

8樓2010-06-11 08:37:04

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

lbing9002

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 21 (小學(xué)生)
金幣: 2197.6
散金: 42
紅花: 1
帖子: 471
在線(xiàn): 138小時(shí)
蟲(chóng)號(hào): 1038678
注冊(cè): 2010-06-09
性別: GG
專(zhuān)業(yè): 計(jì)算機(jī)科學(xué)

強(qiáng)帖

頂��！

回復(fù)此樓

窮則獨(dú)善其身，達(dá)則兼濟(jì)天下。

9樓2012-07-15 18:58:29

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

sg18408926

10樓2012-07-25 16:33:47

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

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

返回列表

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

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

[考研] 302求調(diào)劑 +6	呼呼呼。。。。 2026-03-17	6/300	2026-03-17 23:56 by 星空星月
[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[考研] 341求調(diào)劑 +5	搗蛋豬豬 2026-03-11	7/350	2026-03-17 19:09 by 搗蛋豬豬
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-16	8/400	2026-03-17 19:03 by Wangjingyue
[考研] 332求調(diào)劑 +6	Zz版 2026-03-13	6/300	2026-03-17 17:03 by ruiyingmiao
[考研] 求調(diào)劑，總分315，考的生物醫(yī)藥，一志愿湖南師范大學(xué)。調(diào)劑到任何專(zhuān)業(yè)都可以 +4	小丁想進(jìn)步 2026-03-11	5/250	2026-03-17 16:05 by 外星文明
[考研] 290求調(diào)劑 +6	孔志浩 2026-03-12	11/550	2026-03-17 14:41 by 周舟舟77
[考研] 材料工程專(zhuān)碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學(xué)員h26Tkc
[基金申請(qǐng)] 國(guó)自科面上基金字體 +6	iwuli 2026-03-12	7/350	2026-03-16 21:18 by sculhf
[考研] 機(jī)械專(zhuān)碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 環(huán)境工程調(diào)劑 +6	大可digkids 2026-03-16	6/300	2026-03-16 17:16 by barlinike
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運(yùn)氣yunqi
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 255求調(diào)劑 +3	李嘉慧， 2026-03-12	4/200	2026-03-14 16:58 by 有只貍奴
[考研] 材料080500調(diào)劑求收留 +3	一顆meteor 2026-03-13	3/150	2026-03-14 10:54 by peike
[考研] 招收0805（材料）調(diào)劑 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、
[考研] 281求調(diào)劑 +9	Koxui 2026-03-12	11/550	2026-03-13 20:50 by Koxui
[考研] 考研調(diào)劑 +4	芬達(dá)46 2026-03-12	4/200	2026-03-13 16:04 by ruiyingmiao
[考研] 工科278分求調(diào)劑 +5	周慢熱啊 2026-03-12	7/350	2026-03-13 15:49 by JourneyLucky
[考博] 福州大學(xué)楊黃浩課題組招收2026年專(zhuān)業(yè)學(xué)位博士研究生，2026.03.20截止 +3	Xiangyu_ou 2026-03-12	3/150	2026-03-13 09:36 by duanwu655

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

[交流] 【求助】PBE0 計(jì)算Si、Ge能帶結(jié)構(gòu) 已有6人參與

» 收錄本帖的淘帖專(zhuān)輯推薦

» 猜你喜歡

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

PBE0 計(jì)算Si、Ge能帶結(jié)構(gòu)

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