| 4 | 1/1 | 返回列表 |
| 查看: 903 | 回復(fù): 3 | ||||
gemucai木蟲 (正式寫手)
|
[交流]
交流一下vc-relax過(guò)程中出現(xiàn)的電荷數(shù)錯(cuò)誤的問(wèn)題 已有1人參與
|
|
最近做vc-relax出錯(cuò),報(bào)錯(cuò)如下: from electrons : error # 1 charge is wrong WARNING: integrated charge= 25.58637298, expected= 52.00000000 電荷量只有預(yù)期的一半。我查了PW論壇以前出現(xiàn)的相似問(wèn)題,他們的電荷跟預(yù)期相差無(wú)幾,我這個(gè)就差的太多了。 給出輸入文件: &control calculation='vc-relax' restart_mode='from_scratch', prefix='ret', wf_collect=.true. tstress=.true. tprnfor=.true. pseudo_dir = '...', outdir='...' / &system ibrav= 0, a=1, nat= 6, ntyp= 2, ecutwfc =50.0, nspin = 2, starting_magnetization=0.5, occupations='smearing', smearing='gaussian', degauss=0.02 / &electrons mixing_beta = 0.3 / &IONS bfgs_ndim = 3, pot_extrapolation = "second_order", wfc_extrapolation = "second_order", / &CELL / ATOMIC_SPECIES Mo 95.94 Mo.NCPP_LDA.UPF S 32.0655 S.pz-bhs.UPF ATOMIC_POSITIONS crystal Mo 0.6666666666666666 0.6666666666666666 0.0000000000000000 Mo 0.3333333333333333 0.3333333333333333 0.5000000000000000 S 0.3333333333333333 0.3333333333333333 0.1279036059999967 S 0.3333333333333333 0.3333333333333333 0.8720963940000033 S 0.6666666666666666 0.6666666666666666 0.6279036059999967 S 0.6666666666666666 0.6666666666666666 0.3720963940000033 K_POINTS {automatic} 9 9 5 0 0 0 CELL_PARAMETERS { cubic } 2.7366402800000000 -1.5800000000000001 0.0000000000000000 2.7366402800000000 1.5800000000000001 0.0000000000000000 0.0000000000000000 0.0000000000000000 12.2899999999999991 我贗勢(shì)的階段半徑是不是太小了,不過(guò)我做了測(cè)試,在50的時(shí)候就與45時(shí)的能量相差不到0.1Ry [ Last edited by gemucai on 2011-11-17 at 14:24 ] |
pwscf相關(guān) |
金蟲 (小有名氣)
木蟲 (正式寫手)
金蟲 (小有名氣)
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 308求調(diào)劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂(lè)29 2026-03-18 | 5/250 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 5/250 |
|
|
[考研] 290求調(diào)劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
[考研] 22408 344分 求調(diào)劑 一志愿 華電計(jì)算機(jī)技術(shù) +3 | solanXXX 2026-03-20 | 3/150 |
|
|
[考研] 085410人工智能專碩317求調(diào)劑(0854都可以) +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
[考研] 288求調(diào)劑,一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯! 2026-03-15 | 18/900 |
|
|
[考研] 0854可跨調(diào)劑,一作一項(xiàng)核心論文五項(xiàng)專利,省、國(guó)級(jí)證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考研] 26考研一志愿中國(guó)石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|