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

返回列表

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

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

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

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

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

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

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

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

[考研] 一志愿南航材料專(zhuān)317分求調(diào)劑 +3	炸呀炸呀炸薯?xiàng)l 2026-03-23	3/150	2026-03-23 20:47 by pswait
[考研] 化學(xué)308分求調(diào)劑 +3	你好明天你好 2026-03-23	3/150	2026-03-23 20:11 by macy2011
[考研] 336化工調(diào)劑 +4	王大坦1 2026-03-23	5/250	2026-03-23 18:32 by allen-yin
[考研] 一志愿上海交大生物與醫(yī)藥專(zhuān)碩324分，求調(diào)劑 +5	jiajunX 2026-03-22	5/250	2026-03-23 18:07 by YMU施老師
[考研] 上海電力大學(xué)材料防護(hù)與新材料重點(diǎn)實(shí)驗(yàn)室招收調(diào)劑研究生（材料、化學(xué)、電化學(xué)，環(huán)境） +3	我愛(ài)學(xué)電池 2026-03-23	3/150	2026-03-23 17:16 by AZMK
[考研] 招08考數(shù)學(xué) +6	laoshidan 2026-03-20	14/700	2026-03-23 14:37 by 15614359529
[考研] 08工科 320總分求調(diào)劑 +11	梨花珞晚風(fēng) 2026-03-17	11/550	2026-03-22 17:42 by luoyongfeng
[考研] 275求調(diào)劑 +6	shansx 2026-03-22	8/400	2026-03-22 15:27 by barlinike
[考研] 求調(diào)劑 +5	Zhangbod 2026-03-21	7/350	2026-03-22 13:13 by Zhangbod
[考研] 285求調(diào)劑 +6	ytter 2026-03-22	6/300	2026-03-22 12:09 by 星空星月
[考研] 280求調(diào)劑 +11	咕嚕曉曉 2026-03-18	12/600	2026-03-21 22:40 by ACS Nano——
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +3	邱gl 2026-03-20	7/350	2026-03-21 19:05 by 15709483992
[考研] 265求調(diào)劑 +3	Jack?k?y 2026-03-17	3/150	2026-03-21 03:17 by JourneyLucky
[考研] 301求調(diào)劑 +10	yy要上岸呀 2026-03-17	10/500	2026-03-21 03:14 by JourneyLucky
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3	晨昏線(xiàn)與星海 2026-03-18	3/150	2026-03-21 00:46 by JourneyLucky
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3	hisfailed 2026-03-19	6/300	2026-03-20 23:43 by hisfailed
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +9	Ncdx123456 2026-03-19	9/450	2026-03-20 23:41 by lovewei0727
[考研] 考研調(diào)劑求學(xué)校推薦 +3	伯樂(lè)29 2026-03-18	5/250	2026-03-20 22:59 by JourneyLucky
[考研] 0856調(diào)劑，是學(xué)校就去 +8	sllhht 2026-03-19	9/450	2026-03-20 14:25 by 無(wú)懈可擊111
[考研] 材料學(xué)碩318求調(diào)劑 +5	February_Feb 2026-03-19	5/250	2026-03-19 23:51 by 23Postgrad