| 5 | 1/1 | 返回列表 |
| 查看: 4050 | 回復(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í)的指定輸出沒用,因此寫了個(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è)貌似沒辦法解決。下面是我測(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á)到最小后,就開始增加了。不知道大家有沒有遇到過(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)格沒包括G點(diǎn),但差異也太大了,還不知道是什么原因,是不是輸入文件里哪些參數(shù)沒設(shè)好? [ Last edited by zhang668 on 2010-10-24 at 18:56 ] |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿西安交大材料學(xué)碩(英一數(shù)二)347,求調(diào)劑到高分子/材料相關(guān)專業(yè) +4 | zju51 2026-03-31 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 309求調(diào)劑 +5 | gajsj 2026-03-25 | 6/300 |
|
|
[考研] 085602 307分 求調(diào)劑 +10 | 不知道叫什么! 2026-03-26 | 10/500 |
|
|
[考研] 309求調(diào)劑 +17 | 誰(shuí)不是少年 2026-03-29 | 17/850 |
|
|
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分,英語(yǔ)已過(guò)四六級(jí),六級(jí)521,求調(diào)劑 +10 | 美味可樂雞翅 2026-03-26 | 11/550 |
|
|
[考研] 254材料與化工求調(diào)劑 +3 | 翰冬林楠 2026-03-30 | 4/200 |
|
|
[考研] 求調(diào)劑 生物學(xué) 377分 +6 | zzll03 2026-03-31 | 6/300 |
|
|
[考研] 理學(xué)07化學(xué) 303求調(diào)劑 +8 | 睿08 2026-03-27 | 8/400 |
|
|
[考研] 材料與化工272求調(diào)劑 +25 | 阿斯蒂芬2004 2026-03-28 | 25/1250 |
|
|
[考研] 263求調(diào)劑 +3 | DDDDuu 2026-03-27 | 3/150 |
|
|
[考研] 322求調(diào)劑 +10 | 舊吢 2026-03-24 | 10/500 |
|
|
[考研] 一志愿大連理工大學(xué),機(jī)械工程學(xué)碩,341 +3 | 西瓜田的守望者 2026-03-30 | 3/150 |
|
|
[考研] 266分,求材料相關(guān)專業(yè)調(diào)劑 +10 | 哇呼哼呼哼 2026-03-30 | 12/600 |
|
|
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +6 | 星沉uu 2026-03-26 | 7/350 |
|
|
[考研] 083000環(huán)境科學(xué)與工程調(diào)劑,總分281 +4 | 橙子(勝意) 2026-03-30 | 4/200 |
|
|
[考研] 281求調(diào)劑 +5 | 亞克西good 2026-03-26 | 7/350 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +6 | 丹青奶蓋 2026-03-26 | 8/400 |
|
|
[考研] 085701求調(diào)劑初試286分 +5 | secret0328 2026-03-28 | 5/250 |
|
|
[考研] 本科新能源科學(xué)與工程,一志愿華理能動(dòng)285求調(diào)劑 +3 | AZMK 2026-03-27 | 5/250 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-27 | 8/400 |
|