司馬星空

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 319
散金: 120
紅花: 3
帖子: 36
在線: 323.5小時(shí)
蟲號(hào): 2497523
注冊(cè): 2013-06-05
專業(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ù)此樓

» 猜你喜歡

一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑已經(jīng)有23人回復(fù)
08工科，295，接受跨專業(yè)調(diào)劑已經(jīng)有3人回復(fù)
材料學(xué)碩333求調(diào)劑已經(jīng)有12人回復(fù)
334分一志愿武理材料求調(diào)劑已經(jīng)有16人回復(fù)
0703化學(xué)321分求調(diào)劑已經(jīng)有9人回復(fù)
材料專碩 085600求調(diào)劑已經(jīng)有7人回復(fù)
297求調(diào)劑已經(jīng)有15人回復(fù)
11408總分309，一志愿東南大學(xué)求調(diào)劑，不挑專業(yè) 已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有9人回復(fù)
化工專碩求調(diào)劑已經(jīng)有4人回復(fù)

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

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

chuanghua304

禁蟲 (職業(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

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

1ST強(qiáng)帖: 13
應(yīng)助: 241 (大學(xué)生)
金幣: 25701.7
散金: 602
紅花: 166
帖子: 3782
在線: 873.8小時(shí)
蟲號(hào): 1007127
注冊(cè): 2010-04-27
專業(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í)最常見的，由于對(duì)稱性檢查的問題。

你可試試改用 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ù)的話，還是最好按計(jì)算一組“晶格常數(shù)--總能”的數(shù)據(jù)之后進(jìn)行擬合得到最優(yōu)的晶格常數(shù)，這樣的方式來做。

贊一下

回復(fù)此樓

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

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

souledge

專家顧問 (著名寫手)

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

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

★
感謝參與，應(yīng)助指數(shù) +1
sunyang1988: 金幣+1, 專家考核, 謝謝交流 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)角度不正交的毛病。
換句話說，想要固定角度，除非是三個(gè)90度角，否則在QE中就會(huì)出現(xiàn)正交錯(cuò)誤，并在第一個(gè)SCF后結(jié)束vc-relax的計(jì)算。當(dāng)然，是否有方法克服這個(gè)直接設(shè)置cell_dofree的問題，個(gè)人暫時(shí)不清楚。所以這并不表示QE就一定不能固定至少有一個(gè)非90的角。

贊一下

回復(fù)此樓

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

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

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

souledge

專家顧問 (著名寫手)

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

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

順便補(bǔ)一句官方對(duì)于cell_dofree參數(shù)的解釋：

引用回帖:

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

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

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

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

★
sunyang1988: 金幣+1, 專家考核, 謝謝交流 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)角度不正交的毛病。
換句話說，想要固定角度，除非是三個(gè)90度角 ...

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

贊一下

回復(fù)此樓

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

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

chuanghua304

禁蟲 (職業(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的回帖

司馬星空

銀蟲 (初入文壇)

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

引用回帖:

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

謝謝各位。cell_dofree='all' 測(cè)試過是可行的。我后來想了一下，把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
基矢變了，說明晶格常數(shù)也變了，變化后的應(yīng)該是7.729*1.024584086，不知道這樣是否有問題？

贊一下

回復(fù)此樓

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

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

valenhou001

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

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

是可以的。

回復(fù)此樓

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

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

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

[考研] 286求調(diào)劑 +5	Faune 2026-03-30	5/250	2026-03-30 17:15 by 無際的草原
[考研] 292求調(diào)劑 +12	是妍子也是研子 2026-03-30	12/600	2026-03-30 17:07 by 晴空210210
[考研] 化學(xué)0703 調(diào)劑 306分一志愿211 +8	26要上岸 2026-03-28	8/400	2026-03-30 16:39 by herarysara
[考研] 332求調(diào)劑 +13	032500 2026-03-25	13/650	2026-03-30 16:25 by wangjy2002
[考研] 322求調(diào)劑 +9	宋明欣 2026-03-27	9/450	2026-03-30 16:15 by wangjy2002
[考研] 求調(diào)劑 +7	研研，接電話 2026-03-24	8/400	2026-03-30 14:14 by 白云朵朵飛
[考研] 375求調(diào)劑 +6	雨夏整夜 2026-03-29	6/300	2026-03-30 10:21 by herarysara
[考研] 290求調(diào)劑 +3	dfffsar 2026-03-29	3/150	2026-03-29 22:38 by 毛毛毛阿莫2
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +6	崔wj 2026-03-26	6/300	2026-03-29 01:11 by hanserlol
[考研] 壓國家一區(qū)線，求導(dǎo)師收留，有恩必謝！ +7	迷人的哈哈 2026-03-28	7/350	2026-03-28 16:47 by 催化大白
[考研] 復(fù)試調(diào)劑 +3	raojunqi0129 2026-03-28	3/150	2026-03-28 15:27 by 落睿可思
[考研] 308求調(diào)劑 +7	墨墨漠 2026-03-25	7/350	2026-03-27 14:47 by 狂炫麥當(dāng)當(dāng)
[考研] 315調(diào)劑 +4	0860求調(diào)劑 2026-03-26	5/250	2026-03-27 11:23 by wangjy2002
[考研] 調(diào)劑 +3	李嘉圖·S·路 2026-03-27	3/150	2026-03-27 11:19 by wangjy2002
[考研] 考研調(diào)劑 +9	小蠟新筆 2026-03-26	9/450	2026-03-27 11:10 by 不吃魚的貓
[考研] 281求調(diào)劑 +6	Koxui 2026-03-24	7/350	2026-03-26 15:37 by 無際的草原
[考研] 一志愿南京郵電大學(xué) 288分材料考研求調(diào)劑 +3	jl0720 2026-03-26	3/150	2026-03-26 13:39 by zzll406
[考研] 求b區(qū)院校調(diào)劑 +4	周56 2026-03-24	5/250	2026-03-25 17:12 by yishunmin
[考研] 340求調(diào)劑 +5	話梅糖111 2026-03-24	5/250	2026-03-25 06:53 by ilovexiaobin
[考研] 080500求調(diào)劑 +3	zzzzfan 2026-03-24	3/150	2026-03-24 16:38 by barlinike