| 5 | 1/1 | 返回列表 |
| 查看: 2117 | 回復(fù): 8 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||||
[求助]
Ag體內(nèi)做vc-relax出錯(cuò),不知怎樣改參數(shù)?
|
||||
|
自己做了Ag原子的贗勢,想測試下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 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% |
至尊木蟲 (職業(yè)作家)
|
樓主的基矢(CELLDM和IBRAV)和原子坐標(biāo)設(shè)置都是合適的,因?yàn)橛?jì)算的是fcc的Ag原胞。 出現(xiàn)這樣的錯(cuò)誤,是vc-relax計(jì)算時(shí)最常見的,由于對稱性檢查的問題。 你可試試改用 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ù),這樣的方式來做。 |
禁蟲 (職業(yè)作家)
|
本帖內(nèi)容被屏蔽 |
專家顧問 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +174 |
|
因?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的角。 |

專家顧問 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +174 |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京化工大學(xué)材料與化工(085600)296求調(diào)劑 +19 | 稻妻小編 2026-03-26 | 19/950 |
|
|---|---|---|---|---|
|
[考研] 346求調(diào)劑 一志愿070303有機(jī)化學(xué) +5 | 蘿卜燉青菜 2026-03-28 | 6/300 |
|
|
[考研] 071010 323 分求調(diào)劑 +3 | Baekzhy 2026-03-27 | 3/150 |
|
|
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷 +8 | cccchenso 2026-03-23 | 8/400 |
|
|
[考研] 考研調(diào)劑 +5 | Sanmu-124 2026-03-26 | 5/250 |
|
|
[考研] 367求調(diào)劑 +3 | 芋泥啵啵… 2026-03-28 | 3/150 |
|
|
[考研] 化學(xué)308分調(diào)劑 +7 | 你好明天你好 2026-03-30 | 8/400 |
|
|
[考研] 求調(diào)劑 +10 | 張zz111 2026-03-27 | 11/550 |
|
|
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +8 | 柳溪邊 2026-03-26 | 8/400 |
|
|
[考研] 298求調(diào)劑 +4 | 種圣賜 2026-03-28 | 4/200 |
|
|
[考研] 求調(diào)劑 +7 | 爭取九點(diǎn)睡 2026-03-28 | 8/400 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
[考研] 材料求調(diào)劑 一志愿哈工大總分298分,前三科223分 +5 | dongfang59 2026-03-27 | 5/250 |
|
|
[考研] 070300化學(xué)求調(diào)劑 +4 | 起個(gè)名咋這么難 2026-03-27 | 4/200 |
|
|
[考研] 336材料求調(diào)劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
[考研] 321求調(diào)劑 +6 | Ymlll 2026-03-24 | 6/300 |
|
|
[考研] 290分調(diào)劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
[考研] 302求調(diào)劑 +4 | 錦衣衛(wèi)藤椒 2026-03-25 | 4/200 |
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|
|
[考研] 0854人工智能方向招收調(diào)劑 +4 | 章小魚567 2026-03-24 | 4/200 |
|