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

返回列表

weiyongkai

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

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

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

輸入文件及勢(shì)函數(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計(jì)算的P-V結(jié)果，D為VASP計(jì)算的P-V結(jié)果
Graph1.png

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲(chóng)僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶(hù)自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶(hù)本人，如對(duì)版權(quán)有異議，請(qǐng)聯(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

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

應(yīng)助: 37 (小學(xué)生)
金幣: 834.1
紅花: 6
帖子: 153
在線(xiàn): 92.6小時(shí)
蟲(chóng)號(hào): 1758999
注冊(cè): 2012-04-16
專(zhuān)業(yè): 玻璃材料

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

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

有沒(méi)有換算VASP與QE的壓強(qiáng)單位的不同?

贊一下

回復(fù)此樓

010-8393-9164

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

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

huazhorg

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

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

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

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

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

贊一下(1人)

回復(fù)此樓

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

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

weiyongkai

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

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

引用回帖:

2樓: Originally posted by DFPT at 2013-05-10 17:15:09
有沒(méi)有換算VASP與QE的壓強(qiáng)單位的不同?

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

贊一下

回復(fù)此樓

凌云

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

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

weiyongkai

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

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

引用回帖:

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

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

贊一下

回復(fù)此樓

凌云

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

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

goldenfisher

金蟲(chóng) (著名寫(xiě)手)

1ST強(qiáng)帖: 3
應(yīng)助: 55 (初中生)
貴賓: 0.153
金幣: 5779.1
散金: 2908
紅花: 21
帖子: 1552
在線(xiàn): 351.5小時(shí)
蟲(chóng)號(hào): 117764
注冊(cè): 2005-11-26
性別: GG
專(zhuān)業(yè): 原子和分子物理

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

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

80%是贗勢(shì)的問(wèn)題。
QE里面Al和H的贗勢(shì)那么多，選擇也非常多、

贊一下

回復(fù)此樓

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

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

weiyongkai

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

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

引用回帖:

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

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

贊一下

回復(fù)此樓

凌云

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

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

huazhorg

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

應(yīng)助: 21 (小學(xué)生)
金幣: 4953.8
散金: 5
紅花: 2
帖子: 381
在線(xiàn): 363.2小時(shí)
蟲(chóng)號(hào): 629899
注冊(cè): 2008-10-18
性別: GG
專(zhuān)業(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
我試了其他的贗勢(shì)，但老提示錯(cuò)誤，比如用H.blyp-van_ak.UPF和Al.blyp-n-van_ak.UPF提示：
task# 6
from read-pseudo-nl: error #6
Reading pseudofile (qfunc)
難道從官方網(wǎng)站下下來(lái)的勢(shì)函數(shù)不能直接用么？...

你試試PBE贗勢(shì)，比如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

金蟲(chóng) (著名寫(xiě)手)

1ST強(qiáng)帖: 3
應(yīng)助: 55 (初中生)
貴賓: 0.153
金幣: 5779.1
散金: 2908
紅花: 21
帖子: 1552
在線(xiàn): 351.5小時(shí)
蟲(chóng)號(hào): 117764
注冊(cè): 2005-11-26
性別: GG
專(zhuān)業(yè): 原子和分子物理

引用回帖:

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

贊一下

回復(fù)此樓

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

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

weiyongkai

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

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

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

引用回帖:

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

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

贊一下

回復(fù)此樓

凌云

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

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

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

返回列表

普通表情龍兔虎貓

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

[考博] 申博26年 +4	八6八68 2026-03-19	4/200	2026-03-24 15:49 by 小Ben呵呵
[考研] 299求調(diào)劑 +7	某某某某位 2026-03-21	7/350	2026-03-24 15:24 by cuifj
[考研] 0703化學(xué)調(diào)劑，求導(dǎo)師收 +4	天天好運(yùn)來(lái)上岸?/a> 2026-03-24	4/200	2026-03-24 15:23 by cuifj
[考研] 300求調(diào)劑，材料科學(xué)英一數(shù)二 +3	leaflight 2026-03-24	3/150	2026-03-24 15:10 by 冠c哥
[考研] 081700 調(diào)劑 267分 +9	迷人的哈哈 2026-03-23	9/450	2026-03-24 11:58 by 544594351
[考研] 0856材料專(zhuān)碩353求調(diào)劑 +5	NIFFFfff 2026-03-20	5/250	2026-03-24 11:46 by 544594351
[考研] 一志愿武理材料工程348求調(diào)劑 +6	￣^￣゜汗 2026-03-19	9/450	2026-03-23 19:53 by pswait
[考研] 350求調(diào)劑 +6	weudhdk 2026-03-19	6/300	2026-03-23 15:47 by tangyuan0840221
[考研] 287求調(diào)劑 +8	晨昏線(xiàn)與星海 2026-03-19	9/450	2026-03-22 17:01 by i_cooler
[考研] 324求調(diào)劑 +6	lucky呀呀呀鴨 2026-03-20	6/300	2026-03-22 16:01 by ColorlessPI
[考研] 275求調(diào)劑 +6	shansx 2026-03-22	8/400	2026-03-22 15:27 by barlinike
[考博] 招收博士1-2人 +3	QGZDSYS 2026-03-18	4/200	2026-03-22 10:25 by QGZDSYS
[考研] 材料學(xué)碩333求調(diào)劑 +3	北道巷 2026-03-18	3/150	2026-03-21 18:17 by 學(xué)員8dgXkO
[考研] 307求調(diào)劑 +3	余意卿 2026-03-18	3/150	2026-03-21 17:31 by ColorlessPI
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3	研究生何瑤明 2026-03-18	3/150	2026-03-21 01:30 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é)碩 354求調(diào)劑211或者雙一流 +3	我想要讀研究生 2026-03-20	3/150	2026-03-20 20:13 by JourneyLucky
[考研] 求調(diào)劑 +3	@taotao 2026-03-20	3/150	2026-03-20 19:35 by JourneyLucky
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518
[考研] 0703化學(xué)調(diào)劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011

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

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

» 本帖附件資源列表

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

» 猜你喜歡

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

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

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

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

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

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

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