| 5 | 1/1 | 返回列表 |
| 查看: 4048 | 回復(fù): 18 | ||||
| 【獎(jiǎng)勵(lì)】 本帖被評(píng)價(jià)4次,作者zhang668增加金幣 3.2 個(gè) | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||||
zhang668木蟲 (正式寫手)
|
[資源]
【原創(chuàng)】寫了個(gè)集群下應(yīng)用測(cè)試K點(diǎn)或截?cái)嗄艿哪_本 已有4人參與
|
|||
|
由于我用的是學(xué)校服務(wù)器,scf生成的總能等信息都寫在了任務(wù)提交系統(tǒng)指定的結(jié)果文件output..+(變化ID) 文件中,提交任務(wù)時(shí)的指定輸出沒(méi)用,因此寫了個(gè)在指定任務(wù)提交系統(tǒng)下的腳本測(cè)試截?cái)嗄芎蚄點(diǎn)。 #!/bin/sh for i in 60 70 80 90 100 110 120 130 140 150 160 170 180 190 200 do cat > scf.in < &control calculation= 'scf' restart_mode='from_scratch' pseudo_dir = './' outdir ='./' prefix = 'BiAlO3' tprnfor = .true. / &system ibrav = 1, celldm(1)=7.09, nat=5, ntyp = 3, ecutwfc = $i, ecutrho = 300, / &electrons diagonalization = 'david' conv_thr = 1.0e-8 mixing_beta = 0.7 / ATOMIC_SPECIES Bi 208.98 Bi.pbe-d-mt.UPF Al 26.98 Al.pbe-sp-van.UPF O 15.9994 O.pbe-van_ak.UPF ATOMIC_POSITIONS Bi 0 0 0 Al 0.5 0.5 0.5 O 0.5 0.5 0 O 0.5 0 0.5 O 0 0.5 0.5 K_POINTS AUTOMATIC 5 5 5 0 0 0 ! pw.x -inp scf.in E=`grep 'Harris-Foulkes estimate ' output.* |tail -1 |awk '{printf "%7.6f \n", $4}'` echo $i $E >>Evi.dat done |
材料計(jì)算模擬實(shí)用技巧 | VASP |
木蟲 (正式寫手)
|
得到的Evi.dat文件中,同一組的i,E數(shù)據(jù)重復(fù)出現(xiàn)次數(shù)和并行節(jié)點(diǎn)數(shù)一樣多,這個(gè)貌似沒(méi)辦法解決。下面是我測(cè)試的結(jié)果:(第一列是截?cái)嗄,單位Ry) 60 -397.677201 70 -397.680237 80 -397.685099 90 -397.689479 100 -397.691504 110 -397.692147 120 -397.692890 130 -397.694475 140 -397.696578 150 -397.698359 160 -397.699449 170 -397.700197 180 -397.701290 190 -397.703283 200 -397.706246 發(fā)現(xiàn)一個(gè)奇怪的現(xiàn)象,能量變化幅度不是隨之截?cái)嗄艿脑黾佣饾u減小的,而是在120Ry達(dá)到最小后,就開始增加了。不知道大家有沒(méi)有遇到過(guò)這種情況?是什么原因呢? [ Last edited by zhang668 on 2010-10-24 at 18:38 ] |
木蟲 (正式寫手)
木蟲 (正式寫手)
|
K點(diǎn)測(cè)試結(jié)果,截?cái)嗄苡玫?20Ry,其中用的N X N X N 0 0 0的K網(wǎng)格,(第一列為N) 3 -387.643933 4 -387.652958 5 -387.616159 6 -387.646384 7 -387.615051 8 -387.646694 9 -387.659891 偶數(shù)K網(wǎng)格沒(méi)包括G點(diǎn),但差異也太大了,還不知道是什么原因,是不是輸入文件里哪些參數(shù)沒(méi)設(shè)好? [ Last edited by zhang668 on 2010-10-24 at 18:56 ] |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 323分 食品與營(yíng)養(yǎng)調(diào)劑 +3 | 嘿ooo 2026-03-31 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 085404 22408 315分 +3 | zhuangyan123 2026-03-31 | 3/150 |
|
|
[考研] 一志愿:西北大學(xué),英一數(shù)一408-284分求調(diào)劑 +5 | 12.27 2026-03-27 | 5/250 |
|
|
[考研] 309求調(diào)劑 +15 | 誰(shuí)不是少年 2026-03-29 | 15/750 |
|
|
[考研] 求調(diào)劑 +10 | 家佳佳佳佳佳 2026-03-29 | 10/500 |
|
|
[考研] 329求調(diào)劑,一志愿西北工業(yè)大學(xué),材料工程(085601) +5 | 小小機(jī)靈蟲 2026-03-29 | 11/550 |
|
|
[考研] 求調(diào)劑 +10 | 張zz111 2026-03-27 | 11/550 |
|
|
[考研] 279求調(diào)劑 +4 | 蝶舞輕繞 2026-03-29 | 4/200 |
|
|
[考研] 0856求調(diào)劑 +13 | zhn03 2026-03-25 | 14/700 |
|
|
[碩博家園] 招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑 +4 | IceGuo 2026-03-26 | 5/250 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +6 | 崔wj 2026-03-26 | 6/300 |
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
[考研] 材料277求調(diào)劑 +7 | min3 2026-03-24 | 7/350 |
|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +7 | 閆旭東 2026-03-28 | 9/450 |
|
|
[考研] 復(fù)試調(diào)劑,一志愿南農(nóng)083200食品科學(xué)與工程 +5 | XQTJZ 2026-03-26 | 5/250 |
|
|
[考研] 調(diào)劑推薦 +5 | 清酒714 2026-03-26 | 6/300 |
|
|
[考研] 0703化學(xué)338求調(diào)劑! +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
[考研]
材料調(diào)劑
5+4
|
想要一壺桃花水 2026-03-25 | 10/500 |
|