 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 4038 | 回復: 18 | ||||
| 【獎勵】 本帖被評價4次，作者zhang668增加金幣 3.2 個 | ||||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||||
[資源]
【原創(chuàng)】寫了個集群下應用測試K點或截斷能的腳本 已有4人參與
|
||||
|
由于我用的是學校服務器，scf生成的總能等信息都寫在了任務提交系統(tǒng)指定的結(jié)果文件output..+（變化ID） 文件中，提交任務時的指定輸出沒用，因此寫了個在指定任務提交系統(tǒng)下的腳本測試截斷能和K點。 #!/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 |
回復此樓» 收錄本帖的淘帖專輯推薦
 材料計算模擬實用技巧 | VASP |
» 猜你喜歡
299求調(diào)劑
已經(jīng)有8人回復
-
一志愿北京理工大學本科211材料工程294求調(diào)劑
已經(jīng)有6人回復
-
300求調(diào)劑，材料科學英一數(shù)二
已經(jīng)有8人回復
-
招收生物學/細胞生物學調(diào)劑
已經(jīng)有5人回復
-
070305高分子化學與物理 304分求調(diào)劑
已經(jīng)有7人回復
-
289求調(diào)劑
已經(jīng)有13人回復
-
一志愿哈爾濱工業(yè)大學材料與化工方向336分
已經(jīng)有9人回復
-
081200-11408-276學碩求調(diào)劑
已經(jīng)有6人回復
-
調(diào)劑求院校招收
已經(jīng)有5人回復
-
調(diào)劑310
已經(jīng)有8人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- BET測試中回滯環(huán)變寬,什么原因？
已經(jīng)有10人回復
- 中科院上海應用物理研究所核能項目招聘熔鹽質(zhì)量控制人員（熔鹽理化組）
已經(jīng)有7人回復
- 請教Windows系統(tǒng)下perl腳本的ARGV應用問題，謝謝
已經(jīng)有5人回復
- 氨基樹脂應用方面的求助
已經(jīng)有23人回復
- 求助怎樣寫論文綜述
已經(jīng)有7人回復
- 寫國基過程中，這樣的問題你們有遇到嗎？我前幾年沒遇到。
已經(jīng)有3人回復
- 原子熒光測試遇到問題，怎么辦？
已經(jīng)有5人回復
- 【活動】Show 腳本，贏獎勵
已經(jīng)有23人回復
- UL 黃卡申請 阻燃、RTI、紫外老化的測試方法及標準
已經(jīng)有131人回復
- 請教ORCA的并行及集群上的腳本文件
已經(jīng)有14人回復
- 【求助成功】vasp測試K點收斂性的腳本問題
已經(jīng)有20人回復
9樓2010-11-08 16:27:08
2樓2010-10-24 18:32:43
|
得到的Evi.dat文件中，同一組的i，E數(shù)據(jù)重復出現(xiàn)次數(shù)和并行節(jié)點數(shù)一樣多，這個貌似沒辦法解決。下面是我測試的結(jié)果：（第一列是截斷能，單位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)一個奇怪的現(xiàn)象，能量變化幅度不是隨之截斷能的增加而逐漸減小的，而是在120Ry達到最小后，就開始增加了。不知道大家有沒有遇到過這種情況？是什么原因呢？ [ Last edited by zhang668 on 2010-10-24 at 18:38 ] |
3樓2010-10-24 18:36:54
4樓2010-10-24 18:43:31
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 070305高分子化學與物理 304分求調(diào)劑 +7 | c297914 2026-03-28 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 081200-11408-276學碩求調(diào)劑 +6 | 崔wj 2026-03-26 | 6/300 |
|
|
|
[考研] 332求調(diào)劑 +4 | @MZB382400 2026-03-28 | 4/200 |
|
|
|
[考研] 一志愿華理，數(shù)一英一285求A區(qū)調(diào)劑 +8 | AZMK 2026-03-25 | 12/600 |
|
|
|
[考研] 0856求調(diào)劑 +11 | zhn03 2026-03-25 | 12/600 |
|
|
|
[考研] 藥學105500求調(diào)劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
|
[考研] 085405 考的11408求各位老師帶走 +3 | Qiu學ing 2026-03-28 | 3/150 |
|
|
|
[考研] 292求調(diào)劑 +14 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 15/750 |
|
|
|
[考研] 張芳銘-中國農(nóng)業(yè)大學-環(huán)境工程專碩-298 +4 | 手機用戶 2026-03-26 | 4/200 |
|
|
|
[考研] 0703化學求調(diào)劑，各位老師看看我�。�！ +5 | 祁祺祺 2026-03-25 | 5/250 |
|
|
|
[有機交流]
高溫高壓反應求助
 10+4
|
chibby 2026-03-25 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +4 | 零八# 2026-03-27 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | 劉柯@ 2026-03-24 | 4/200 |
|
|
|
[考研] 0703化學338求調(diào)劑！ +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
|
[考研] 276求調(diào)劑。有半年電池和半年高分子實習經(jīng)歷 +10 | 材料學257求調(diào)劑 2026-03-23 | 11/550 |
|
|
|
[考研]
材料調(diào)劑
5+4
|
想要一壺桃花水 2026-03-25 | 10/500 |
|
|
|
[考研] 一志愿中南大學化學學碩0703總分337求調(diào)劑 +7 | niko- 2026-03-22 | 7/350 |
|
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|
|
|
[考研] 材料專碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|