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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

weiyongkai

木蟲 (職業(yè)作家)

應(yīng)助: 13 (小學(xué)生)
金幣: 3192
散金: 209
紅花: 3
帖子: 3186
在線: 228.1小時
蟲號: 1104094
注冊: 2010-09-20
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

[求助] QE計算在某一壓強下的晶體參數(shù)要比相同條件下的VASP小的很多，是怎么回事呢？

輸入文件及勢函數(shù)（官方網(wǎng)站下的）如下：
#! /bin/sh
#
# SCF at dense k-mesh, good enough for electronic DOS
#
cat > al.scf.fit.in <<!
&control
calculation='vc-relax',
restart_mode='from_scratch',
prefix='al',
pseudo_dir = './ps',
outdir='./tmp',
wf_collect =.true.,
etot_conv_thr = 1.0d-6 ,
forc_conv_thr = 1.0d-4 ,
tstress = .true. ,
tprnfor = .true. ,
nstep = 200 ,
/
&system
   ibrav=14,
   a=3.169,b=3.169,c=3.169,
   cosab=0,cosac=0,cosbc=0,
   nat= 8, ntyp= 2,
ecutwfc =100.0,
ecutrho = 800,
occupations='smearing', smearing='methfessel-paxton', degauss=0.05,
la2F=.true.,
/
&electrons
diagonalization='david',
conv_thr =  1.0d-14,
mixing_beta = 0.7,
      mixing_ndim = 10,
diagonalization='david',
/
&IONS
ion_dynamics='bfgs',
ion_positions = 'from_input',
pot_extrapolation = "second_order",
wfc_extrapolation = "second_order",
/
&cell
            cell_dynamics = 'bfgs' ,
            cell_dofree = 'all' ,
            press = 250,
            cell_factor = 5.0,
/
ATOMIC_SPECIES
  H  1.008 H.pz-vbc.UPF
  Al 26.982  Al.pz-vbc.UPF
ATOMIC_POSITIONS {crystal}
  H 0.25 0.00 0.50
  H  -0.25  -0.00 0.50
  H 0.50 0.25  -0.00
  H 0.50  -0.25 0.00
  H 0.00 0.50 0.25
  H 0.00 0.50  -0.25
Al 0.00 0.00 0.00
Al 0.50 0.50 0.50
K_POINTS {automatic}
16 16 16  0 0 0
!
mpirun -np 12 pw.x<al.scf.fit.in > 25.out

B為QE計算的P-V結(jié)果，D為VASP計算的P-V結(jié)果
Graph1.png

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : ps.rar

2013-05-10 11:28:21, 16.66 K

凌云

1樓 2013-05-10 11:28:22

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

DFPT

金蟲 (小有名氣)

應(yīng)助: 37 (小學(xué)生)
金幣: 834.1
紅花: 6
帖子: 153
在線: 92.6小時
蟲號: 1758999
注冊: 2012-04-16
專業(yè): 玻璃材料

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
weiyongkai: 金幣+10, ★有幫助 2013-05-13 09:15:05

有沒有換算VASP與QE的壓強單位的不同?

贊一下

回復(fù)此樓

010-8393-9164

2樓2013-05-10 17:15:09

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

huazhorg

鐵桿木蟲 (正式寫手)

應(yīng)助: 21 (小學(xué)生)
金幣: 4953.8
散金: 5
紅花: 2
帖子: 381
在線: 363.2小時
蟲號: 629899
注冊: 2008-10-18
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
weiyongkai: 金幣+10, ★有幫助 2013-05-13 09:16:37

LDA贗勢會低估晶格常數(shù)。

贊一下(1人)

回復(fù)此樓

3樓2013-05-12 10:36:53

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

weiyongkai

木蟲 (職業(yè)作家)

應(yīng)助: 13 (小學(xué)生)
金幣: 3192
散金: 209
紅花: 3
帖子: 3186
在線: 228.1小時
蟲號: 1104094
注冊: 2010-09-20
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

引用回帖:

2樓: Originally posted by DFPT at 2013-05-10 17:15:09
有沒有換算VASP與QE的壓強單位的不同?

換算了啊，QE中的壓強單位為Kbar，設(shè)置的時候是按1GPa=10Kbar設(shè)置的啊，這也小的太多了�。�

贊一下

回復(fù)此樓

凌云

4樓2013-05-13 09:14:54

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

weiyongkai

木蟲 (職業(yè)作家)

應(yīng)助: 13 (小學(xué)生)
金幣: 3192
散金: 209
紅花: 3
帖子: 3186
在線: 228.1小時
蟲號: 1104094
注冊: 2010-09-20
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

引用回帖:

3樓: Originally posted by huazhorg at 2013-05-12 10:36:53
LDA贗勢會低估晶格常數(shù)。

是的，可是從結(jié)果看來低估的有些大，需要轉(zhuǎn)換勢函數(shù)么？可以用abinit勢函數(shù)轉(zhuǎn)化過來用么？

贊一下

回復(fù)此樓

凌云

5樓2013-05-13 09:16:29

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

goldenfisher

金蟲 (著名寫手)

1ST強帖: 3
應(yīng)助: 55 (初中生)
貴賓: 0.153
金幣: 5779.1
散金: 2908
紅花: 21
帖子: 1552
在線: 351.5小時
蟲號: 117764
注冊: 2005-11-26
性別: GG
專業(yè): 原子和分子物理

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
weiyongkai: 金幣+10, ★有幫助 2013-05-13 11:31:03
weiyongkai: 金幣+10, ★★★很有幫助, 謝謝提示，正在一個一個嘗試：） 2013-05-14 08:54:05

80%是贗勢的問題。
QE里面Al和H的贗勢那么多，選擇也非常多、

贊一下

回復(fù)此樓

6樓2013-05-13 10:18:42

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

weiyongkai

木蟲 (職業(yè)作家)

應(yīng)助: 13 (小學(xué)生)
金幣: 3192
散金: 209
紅花: 3
帖子: 3186
在線: 228.1小時
蟲號: 1104094
注冊: 2010-09-20
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

引用回帖:

6樓: Originally posted by goldenfisher at 2013-05-13 10:18:42
80%是贗勢的問題。
QE里面Al和H的贗勢那么多，選擇也非常多、

我試了其他的贗勢，但老提示錯誤，比如用H.blyp-van_ak.UPF和Al.blyp-n-van_ak.UPF提示：
task# 6
from read-pseudo-nl: error #6
Reading pseudofile (qfunc)
難道從官方網(wǎng)站下下來的勢函數(shù)不能直接用么？

贊一下

回復(fù)此樓

凌云

7樓2013-05-13 11:30:54

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

huazhorg

鐵桿木蟲 (正式寫手)

應(yīng)助: 21 (小學(xué)生)
金幣: 4953.8
散金: 5
紅花: 2
帖子: 381
在線: 363.2小時
蟲號: 629899
注冊: 2008-10-18
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
weiyongkai: 金幣+10, ★★★很有幫助, 謝謝幫助！ 2013-05-14 08:52:51

引用回帖:

7樓: Originally posted by weiyongkai at 2013-05-13 11:30:54
我試了其他的贗勢，但老提示錯誤，比如用H.blyp-van_ak.UPF和Al.blyp-n-van_ak.UPF提示：
task# 6
from read-pseudo-nl: error #6
Reading pseudofile (qfunc)
難道從官方網(wǎng)站下下來的勢函數(shù)不能直接用么？...

你試試PBE贗勢，比如Al.pbe-n-van.UPF和H.pbe-van_ak.UPF

贊一下

回復(fù)此樓

8樓2013-05-14 08:36:48

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

goldenfisher

金蟲 (著名寫手)

1ST強帖: 3
應(yīng)助: 55 (初中生)
貴賓: 0.153
金幣: 5779.1
散金: 2908
紅花: 21
帖子: 1552
在線: 351.5小時
蟲號: 117764
注冊: 2005-11-26
性別: GG
專業(yè): 原子和分子物理

引用回帖:

肯定可以直接用的。你的設(shè)置可能有些問題。另外，帶ak字樣的是從CPMD中轉(zhuǎn)換過來的，有些小問題。

贊一下

回復(fù)此樓

9樓2013-05-14 16:07:54

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

weiyongkai

木蟲 (職業(yè)作家)

應(yīng)助: 13 (小學(xué)生)
金幣: 3192
散金: 209
紅花: 3
帖子: 3186
在線: 228.1小時
蟲號: 1104094
注冊: 2010-09-20
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

★
ben_ladeng: 金幣+1, 歡迎分享您的經(jīng)驗 2013-05-15 07:37:29

引用回帖:

9樓: Originally posted by goldenfisher at 2013-05-14 16:07:54
肯定可以直接用的。你的設(shè)置可能有些問題。另外，帶ak字樣的是從CPMD中轉(zhuǎn)換過來的，有些小問題。...

我也懷疑是我的設(shè)置的問題，另外我覺得勢函數(shù)的選擇應(yīng)該要根據(jù)所要計算的要求進行選擇，謝謝

贊一下

回復(fù)此樓

凌云

10樓2013-05-14 16:54:22

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

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

返回列表

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

[考研] 一志愿華東理工大學(xué)081700，初試分?jǐn)?shù)271 +4	kotoko_ik 2026-03-23	5/250	2026-03-23 23:23 by 呆呆師姐
[考研] 291求調(diào)劑 +8	hhhhxn.. 2026-03-23	8/400	2026-03-23 23:15 by peike
[考研] 384求調(diào)劑 +3	子系博 2026-03-22	6/300	2026-03-23 21:45 by 子系博
[考研] 265求調(diào)劑 +10	梁梁校校 2026-03-17	10/500	2026-03-23 21:17 by 一切OK
[考研] 材料專碩英一數(shù)二306 +8	z1z2z3879 2026-03-18	8/400	2026-03-23 20:49 by baobaoye
[考研] 308求調(diào)劑 +3	墨墨漠 2026-03-21	3/150	2026-03-22 16:54 by i_cooler
[考研] 324求調(diào)劑 +6	lucky呀呀呀鴨 2026-03-20	6/300	2026-03-22 16:01 by ColorlessPI
[基金申請] 山東省面上項目限額評審 +4	石瑞0426 2026-03-19	4/200	2026-03-22 08:50 by Wei_ren
[考研] 求調(diào)劑 +3	13341 2026-03-20	3/150	2026-03-21 18:28 by 學(xué)員8dgXkO
[考研] 材料學(xué)碩333求調(diào)劑 +3	北道巷 2026-03-18	3/150	2026-03-21 18:17 by 學(xué)員8dgXkO
[考研] 302求調(diào)劑 +12	呼呼呼。。。。 2026-03-17	12/600	2026-03-21 17:29 by ColorlessPI
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +9	Ncdx123456 2026-03-19	9/450	2026-03-20 23:41 by lovewei0727
[考研] 317求調(diào)劑 +5	申子申申 2026-03-19	9/450	2026-03-20 22:26 by JourneyLucky
[考研] 329求調(diào)劑 +9	想上學(xué)吖吖 2026-03-19	9/450	2026-03-20 22:01 by luoyongfeng
[考研] 295材料求調(diào)劑，一志愿武漢理工085601專碩 +5	Charlieyq 2026-03-19	5/250	2026-03-20 20:35 by JourneyLucky
[考研] 261求B區(qū)調(diào)劑，科研經(jīng)歷豐富 +3	牛奶很忙 2026-03-20	4/200	2026-03-20 19:34 by JourneyLucky
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生（接收計算機專業(yè)調(diào)劑） +3	123阿標(biāo) 2026-03-17	3/150	2026-03-20 15:58 by 飛行琦
[考研] 0856調(diào)劑，是學(xué)校就去 +8	sllhht 2026-03-19	9/450	2026-03-20 14:25 by 無懈可擊111
[考研] 生物學(xué)調(diào)劑招人�。�！ +3	山海天嵐 2026-03-17	4/200	2026-03-19 21:34 by 怎么釋懷
[考博] 申博26年 +3	八6八68 2026-03-19	3/150	2026-03-19 19:43 by nxgogo

24小時熱門版塊排行榜

[求助] QE計算在某一壓強下的晶體參數(shù)要比相同條件下的VASP小的很多，是怎么回事呢？

» 本帖附件資源列表

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

» 猜你喜歡

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

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

[求助] QE計算在某一壓強下的晶體參數(shù)要比相同條件下的VASP小的很多，是怎么回事呢？