司馬星空

應(yīng)助: 0 (幼兒園)
金幣: 319
散金: 120
紅花: 3
帖子: 36
在線(xiàn): 323.5小時(shí)
蟲(chóng)號(hào): 2497523
注冊(cè): 2013-06-05
專(zhuān)業(yè): 心理學(xué)史

[求助] Ag體內(nèi)做vc-relax出錯(cuò)，不知怎樣改參數(shù)？

自己做了Ag原子的贗勢(shì)，想測(cè)試下VC-relax的晶格常數(shù)變化大不大�？墒怯行╁e(cuò)誤，不知道怎么樣修改參數(shù)。
&control
calculation='vc-relax'
restart_mode='from_scratch',
pseudo_dir = '/',
outdir='/',
prefix='Ag'
/
&system
ibrav= 2, celldm(1) =7.7290
nat=1, ntyp=1,
ecutwfc =45
occupations='smearing', smearing='methfessel-paxton', degauss=0.04
/
&electrons
  mixing_beta = 0.4
mixing_ndim =12
conv_thr =  1.0d-8
diagonalization ='cg'
/
&ions
upscale=10
/
&CELL
cell_dynamics = 'bfgs' ,
press = 0.00 ,
cell_dofree='xyz',
/
ATOMIC_SPECIES
Ag  107.9 Ag.rel-pbe.newUPF
ATOMIC_POSITIONS
Ag    0.00000000  0.00000000 0.00000000
K_POINTS (automatic)
8 8 8 0 0 0

輸出結(jié)果：
CELL_PARAMETERS (alat=  7.72900000)
  -0.507305951 0.003652975 0.496347025
  -0.003652975 0.507305951 0.496347025
  -0.503652975 0.503652975  -0.007305951

ATOMIC_POSITIONS (alat)
Ag    0.000000000 0.000000000 0.000000000

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   from checkallsym : error #       2
   not orthogonal operation
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

回復(fù)此樓

» 猜你喜歡

277跪求調(diào)劑已經(jīng)有11人回復(fù)
一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑已經(jīng)有18人回復(fù)
求調(diào)劑已經(jīng)有6人回復(fù)
0703化學(xué) 已經(jīng)有14人回復(fù)
面上5B能上會(huì)嗎？已經(jīng)有7人回復(fù)
0703本科鄭州大學(xué)求調(diào)劑已經(jīng)有7人回復(fù)
286分調(diào)劑已經(jīng)有3人回復(fù)
0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可已經(jīng)有15人回復(fù)
317分一志愿南理工材料工程本科湖工大求調(diào)劑已經(jīng)有11人回復(fù)
332求調(diào)劑已經(jīng)有12人回復(fù)

1樓 2013-07-01 20:47:28

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

chuanghua304

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

★
sunyang1988: 金幣+1, 謝謝交流 2013-07-04 10:06:37

本帖內(nèi)容被屏蔽

2樓2013-07-02 14:34:58

已閱回復(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)

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

★ ★
感謝參與，應(yīng)助指數(shù) +1
sunyang1988: 金幣+2, 謝謝幫助 2013-07-04 10:06:46
sunyang1988: 應(yīng)助指數(shù)+1 2013-07-04 10:07:05

樓主的基矢（CELLDM和IBRAV)和原子坐標(biāo)設(shè)置都是合適的，因?yàn)橛?jì)算的是fcc的Ag原胞。
出現(xiàn)這樣的錯(cuò)誤，是vc-relax計(jì)算時(shí)最常見(jiàn)的，由于對(duì)稱(chēng)性檢查的問(wèn)題。

你可試試改用 ibrav=0,另外手動(dòng)設(shè)置初始的fcc的晶格基矢在輸入文件中：
CELL_PARAMETERS {alat}
  -0.507305951 0.003652975 0.496347025
  -0.003652975 0.507305951 0.496347025
  -0.503652975 0.503652975  -0.007305951

你要計(jì)算晶格常數(shù)的話(huà)，還是最好按計(jì)算一組“晶格常數(shù)--總能”的數(shù)據(jù)之后進(jìn)行擬合得到最優(yōu)的晶格常數(shù)，這樣的方式來(lái)做。

贊一下

回復(fù)此樓

3樓2013-07-02 16:01:47

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

souledge

專(zhuān)家顧問(wèn) (著名寫(xiě)手)

專(zhuān)家經(jīng)驗(yàn): +174
1ST強(qiáng)帖: 19
應(yīng)助: 303 (大學(xué)生)
金幣: 11306.4
散金: 1123
紅花: 108
帖子: 2814
在線(xiàn): 916.1小時(shí)
蟲(chóng)號(hào): 941486
注冊(cè): 2010-01-12
專(zhuān)業(yè): 結(jié)構(gòu)陶瓷
管轄: 第一性原理

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

★
感謝參與，應(yīng)助指數(shù) +1
sunyang1988: 金幣+1, 專(zhuān)家考核, 謝謝交流 2013-07-04 10:07:36

因?yàn)镕CC的Primitive Cell的a, b, c三軸的夾角為60度，而不是90度，所以當(dāng)LZ使用&cell下的cell_dofree = 'xyz'時(shí)，就是固定了三個(gè)角度，就會(huì)出現(xiàn)角度不正交的毛病。
換句話(huà)說(shuō)，想要固定角度，除非是三個(gè)90度角，否則在QE中就會(huì)出現(xiàn)正交錯(cuò)誤，并在第一個(gè)SCF后結(jié)束vc-relax的計(jì)算。當(dāng)然，是否有方法克服這個(gè)直接設(shè)置cell_dofree的問(wèn)題，個(gè)人暫時(shí)不清楚。所以這并不表示QE就一定不能固定至少有一個(gè)非90的角。

贊一下

回復(fù)此樓

思想重于技巧，內(nèi)涵重于表象

4樓2013-07-02 18:34:30

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

souledge

專(zhuān)家顧問(wèn) (著名寫(xiě)手)

專(zhuān)家經(jīng)驗(yàn): +174
1ST強(qiáng)帖: 19
應(yīng)助: 303 (大學(xué)生)
金幣: 11306.4
散金: 1123
紅花: 108
帖子: 2814
在線(xiàn): 916.1小時(shí)
蟲(chóng)號(hào): 941486
注冊(cè): 2010-01-12
專(zhuān)業(yè): 結(jié)構(gòu)陶瓷
管轄: 第一性原理

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

順便補(bǔ)一句官方對(duì)于cell_dofree參數(shù)的解釋?zhuān)?fieldset>引用回帖:

BEWARE: if axis are not orthogonal, some of these options do not work ( symmetry is broken). If you are not happy with them, edit subroutine init_dofree in file Module/cell_base.f90

贊一下

回復(fù)此樓

思想重于技巧，內(nèi)涵重于表象

5樓2013-07-02 18:38:39

已閱回復(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)

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

★
sunyang1988: 金幣+1, 專(zhuān)家考核, 謝謝交流 2013-07-04 10:07:53

引用回帖:

4樓: Originally posted by souledge at 2013-07-02 19:34:30
因?yàn)镕CC的Primitive Cell的a, b, c三軸的夾角為60度，而不是90度，所以當(dāng)LZ使用&cell下的cell_dofree = 'xyz'時(shí)，就是固定了三個(gè)角度，就會(huì)出現(xiàn)角度不正交的毛病。
換句話(huà)說(shuō)，想要固定角度，除非是三個(gè)90度角 ...

你說(shuō)的cell_dofree='xyz'是一個(gè)問(wèn)題，不適用于FCC。改用cell_dofree='all'試試。

贊一下

回復(fù)此樓

6樓2013-07-02 20:14:12

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

chuanghua304

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

★
sunyang1988: 金幣+1, 謝謝交流 2013-07-04 10:08:06

本帖內(nèi)容被屏蔽

7樓2013-07-02 22:34:19

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

司馬星空

銀蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 319
散金: 120
紅花: 3
帖子: 36
在線(xiàn): 323.5小時(shí)
蟲(chóng)號(hào): 2497523
注冊(cè): 2013-06-05
專(zhuān)業(yè): 心理學(xué)史

引用回帖:

6樓: Originally posted by valenhou001 at 2013-07-02 20:14:12
你說(shuō)的cell_dofree='xyz'是一個(gè)問(wèn)題，不適用于FCC。改用cell_dofree='all'試試。...

謝謝各位。cell_dofree='all' 測(cè)試過(guò)是可行的。我后來(lái)想了一下，把Ag體內(nèi)的原胞取成了簡(jiǎn)單立方ibrav=1 ，
ATOMIC_POSITIONS
Ag  0.00000000  0.00000000 0.00000000
Ag  0.50000000  0.50000000 0.00000000
Ag  0.00000000  0.50000000 0.50000000
Ag  0.50000000  0.00000000 0.50000000
能正常運(yùn)行，輸出結(jié)果：

CELL_PARAMETERS (alat=  7.72900000)
1.024584086 0.000000000 0.000000000
0.000000000 1.024584086 0.000000000
0.000000000 0.000000000 1.024584086
基矢變了，說(shuō)明晶格常數(shù)也變了，變化后的應(yīng)該是7.729*1.024584086，不知道這樣是否有問(wèn)題？

贊一下

回復(fù)此樓

8樓2013-07-05 10:14:13

已閱回復(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)

是可以的。

回復(fù)此樓

9樓2013-07-05 10:32:54

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

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

[考研] 277跪求調(diào)劑 +7	1915668 2026-03-27	11/550	2026-03-30 12:49 by fangnagu
[考研] 327求調(diào)劑 +4	小卡不卡. 2026-03-29	4/200	2026-03-30 06:15 by wxiongid
[考研] 311求調(diào)劑 +6	冬十三 2026-03-24	6/300	2026-03-29 20:45 by 無(wú)際的草原
[考研] 總分293求調(diào)劑 +8	加一一九 2026-03-25	11/550	2026-03-29 19:53 by 無(wú)際的草原
[考研] 329求調(diào)劑，一志愿西北工業(yè)大學(xué)，材料工程（085601） +4	小小機(jī)靈蟲(chóng) 2026-03-29	10/500	2026-03-29 18:20 by 無(wú)際的草原
[考研] 一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑 +8	mikasa的圍巾 2026-03-28	8/400	2026-03-29 12:48 by 無(wú)際的草原
[考研] 340求調(diào)劑 +6	Amber00 2026-03-26	6/300	2026-03-29 12:06 by 無(wú)際的草原
[考研] 298求調(diào)劑 +4	種圣賜 2026-03-28	4/200	2026-03-29 08:42 by q1092522407
[考研] 083000學(xué)碩274求調(diào)劑 +8	Li李魚(yú) 2026-03-26	8/400	2026-03-28 20:33 by 加油向未來(lái)啊
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +8	邱gl 2026-03-27	8/400	2026-03-28 12:42 by 唐沐兒
[考研] 085600材料與化工調(diào)劑 +10	A-哆啦Z夢(mèng) 2026-03-23	16/800	2026-03-27 15:13 by caszguilin
[考研] 085601 材料工程 313分求調(diào)劑 +5	Ong3 2026-03-27	5/250	2026-03-27 12:24 by goldfish51
[考研] 316求調(diào)劑 +5	Pigcasso 2026-03-24	5/250	2026-03-27 12:10 by zhshch
[考研] 考研調(diào)劑 +9	小蠟新筆 2026-03-26	9/450	2026-03-27 11:10 by 不吃魚(yú)的貓
[論文投稿] Journal of Mechanical Science and Technology +3	Russ_ss 2026-03-25	5/250	2026-03-27 10:49 by 陸小果畫(huà)大餅
[考研] 343求調(diào)劑 +4	贈(zèng)我一本書(shū) 2026-03-23	4/200	2026-03-27 00:40 by wxiongid
[考研] 材料考研求調(diào)劑 +3	Dendel 2026-03-23	6/300	2026-03-26 17:51 by fmesaito
[考研] 303求調(diào)劑 +6	藍(lán)山月 2026-03-25	6/300	2026-03-25 22:47 by 418490947
[考研] 0854人工智能方向招收調(diào)劑 +4	章小魚(yú)567 2026-03-24	4/200	2026-03-25 13:29 by 2177681040
[考研] 求調(diào)劑 +6	研研，接電話(huà) 2026-03-24	7/350	2026-03-24 17:01 by barlinike