 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 552 | 回復(fù): 2 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
qinform木蟲 (著名寫手)
|
[求助]
求助 晶格優(yōu)化腳本 一個命令格式的修改
|
||
|
在vasp入門指南里 晶格優(yōu)化部分，有一個腳本 #!/bin/sh rm WAVECAR for i in 150 200 250 300 350 400 do cat > INCAR <<! SYSTEM = Si-Diamond ENCUT = $i ISTART = 0 ; ICHARG = 2 ISMEAR = -5 PREC = Accurate ! echo "ENCUT = $i eV" ; time vasp E=‘grep "TOTEN" OUTCAR | tail -1 | awk ’{printf "%12.6f \n", $5 }’‘ echo $i $E >>comment done 我想修改這個腳本，我用的程序輸出文件名是report, 。。。。。。。。。。。。。。。。。。。。。。。。 --------------------------------------------------- Ewald = -.24937674101229E+02 Alpha = 0.44417481059890E+01 E_extV = 0.00000000000000E+00 0.0000E+00 E_NSC = 0.55244485925725E+01 -.1520E-04 E[-rho*V_Hxc]= -.32491868458856E+01 0.1217E-04 E_Hxc = -.47988207613946E+01 0.3027E-05 -TS = -.20308418661012E-29 0.0000E+00 E_tot = -.23019485009947E+02 -.2602E-08 --------------------------------------------------- Zero temp. E_tot =-.23019485009947E+02 Using formula: E_tot(T)+TS/(N+2), N= 0 --------------------------------------------------- E_Hart,E_xc,E_ion =0.2374022815E+01 -.7172843576E+01 -.7078483387E+01 E_rhoVext,E_IVext =0.0000000000E+00 0.0000000000E+00 E_psiV,E_dDrho =-.3829296541E+01 0.0000000000E+00 ave(vtot):v0 =-.4925372855E+00 ave(V_ion_s(or p,d))=ave(V_Hatree)=0; ave(Vtot)=ave(V_xc)=v0 --------------------------------------------------- 弛豫之后輸出總能的一行 Zero temp. E_tot =-.23019485009947E+02 Using formula: E_tot(T)+TS/(N+2), N= 0 。。。。。。。。。。。。。。。。。。 E=‘grep "Zero temp. E_tot“ report | tail -1 | awk ’{printf "%12.6f \n", $5 }’‘ 不能正確輸出總能的數(shù)值，輸出為0。 試了grep "Zero temp. E_tot“ report | tail -1 能夠正確找出 Zero temp. E_tot =-.23019485009947E+02 Using formula: E_tot(T)+TS/(N+2), N= 0 但是， | awk ’{printf "%12.6f \n", $5}‘ 不知道怎么修改。 請幫忙。 |
回復(fù)此樓» 猜你喜歡
288資源與環(huán)境專碩求調(diào)劑，不限專業(yè)，有學(xué)上就行
已經(jīng)有20人回復(fù)
-
311（085601）求調(diào)劑
已經(jīng)有11人回復(fù)
-
材料調(diào)劑
已經(jīng)有12人回復(fù)
-
一志愿 南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩
已經(jīng)有3人回復(fù)
-
0856材料與化工調(diào)劑，339
已經(jīng)有6人回復(fù)
-
08工科275分求調(diào)劑
已經(jīng)有6人回復(fù)
-
334分 一志愿武理-080500 材料求調(diào)劑
已經(jīng)有10人回復(fù)
-
324分 085600材料與化工
已經(jīng)有7人回復(fù)
-
各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩
已經(jīng)有11人回復(fù)
-
309求調(diào)劑
已經(jīng)有6人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 在comsol里面怎么將一個模塊的計算結(jié)果用到另一個應(yīng)力模塊作為加載源
已經(jīng)有9人回復(fù)
- 如何在MS的模型中的晶胞一個棱上（只有一個原子）摻雜替代這個原子
已經(jīng)有13人回復(fù)
- BASH腳本的一個小問題
已經(jīng)有3人回復(fù)
- 關(guān)于優(yōu)化晶格常數(shù)
已經(jīng)有4人回復(fù)
- ENCUT測試腳本
已經(jīng)有17人回復(fù)
- 【求助】lammps最適合安裝在哪一個版本的linux系統(tǒng)中？
已經(jīng)有7人回復(fù)
- 【求助】vasp優(yōu)化時選擇了isif=3 怎么限制晶格的一個維度呢？
已經(jīng)有16人回復(fù)
- 【求助】VASP怎么做晶格參數(shù)的結(jié)構(gòu)優(yōu)化
已經(jīng)有17人回復(fù)
- 【求助】如何根據(jù)cell文件寫poscar
已經(jīng)有11人回復(fù)
- 【求助】vasp優(yōu)化晶格參數(shù)a
已經(jīng)有13人回復(fù)
- 【討論】lammps里的boundary的位置能夠隨意更改？
已經(jīng)有14人回復(fù)
- 【求助】如何從Murnaghan方程的結(jié)果中得到優(yōu)化的晶格參數(shù)？
已經(jīng)有11人回復(fù)
- 【求助成功】晶格優(yōu)化中，murn.f的使用方法
已經(jīng)有7人回復(fù)
- 【求助】侯老師的vasp入門的晶格優(yōu)化部分沒看懂，問一下
已經(jīng)有25人回復(fù)
- 【求助】求證晶格常數(shù)的優(yōu)化
已經(jīng)有15人回復(fù)
- 【求助】關(guān)于晶胞中取代一個原子的優(yōu)化問題
已經(jīng)有14人回復(fù)
- 【求助】幫忙編個用于VASP小程序
已經(jīng)有22人回復(fù)
- 【求助】如何固定某幾個原子不動？
已經(jīng)有20人回復(fù)
jerkwin
專家顧問 (正式寫手)
-
專家經(jīng)驗: +14 - 應(yīng)助: 454 (碩士)
- 金幣: 20699.1
- 散金: 148
- 紅花: 81
- 帖子: 813
- 在線: 2648.3小時
- 蟲號: 1023452
- 注冊: 2010-05-19
- 專業(yè): 理論和計算化學(xué)
- 管轄: 分子模擬
3樓2013-09-22 00:56:16
jerkwin
專家顧問 (正式寫手)
-
專家經(jīng)驗: +14 - 應(yīng)助: 454 (碩士)
- 金幣: 20699.1
- 散金: 148
- 紅花: 81
- 帖子: 813
- 在線: 2648.3小時
- 蟲號: 1023452
- 注冊: 2010-05-19
- 專業(yè): 理論和計算化學(xué)
- 管轄: 分子模擬
; printf "%12.6f \n", $5}' 2樓2013-09-21 22:10:29
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿 南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +3 | @taotao 2026-03-31 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 070300化學(xué)279求調(diào)劑 +9 | 哈哈哈^_^ 2026-03-31 | 11/550 |
|
|
|
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +11 | 柳溪邊 2026-03-26 | 12/600 |
|
|
|
[考研] 省雙一流重點一本大學(xué)招收調(diào)劑 +3 | wwwwffffff 2026-03-31 | 5/250 |
|
|
|
[考研] 0710生物學(xué)求調(diào)劑！ +5 | 敘述文 2026-03-31 | 5/250 |
|
|
[考研] 343求調(diào)劑 +8 | 愛羈絆 2026-03-28 | 8/400 |
|
|
|
[考研] 化學(xué)工程085602 305分求調(diào)劑 +28 | RichLi_ 2026-03-25 | 36/1800 |
|
|
|
[考研] 335求調(diào)劑 +3 | 321* 2026-03-31 | 3/150 |
|
|
|
[考研] 英一數(shù)一408，總分284，二戰(zhàn)真誠求調(diào)劑 +3 | 12.27 2026-03-30 | 5/250 |
|
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +11 | 多久上課 2026-03-27 | 12/600 |
|
|
|
[考研] 材料化工340求調(diào)劑 +3 | jhx777 2026-03-30 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +10 | 張zz111 2026-03-27 | 11/550 |
|
|
|
[考研] 085600，材料與化工321分求調(diào)劑 +10 | 大饞小子 2026-03-28 | 10/500 |
|
|
|
[考研] 279求調(diào)劑 +4 | 蝶舞輕繞 2026-03-29 | 4/200 |
|
|
|
[考研] 356求調(diào)劑 +3 | gysy?s?a 2026-03-28 | 3/150 |
|
|
|
[考研] 394求調(diào)劑 +3 | 好事多磨靜候佳?/a> 2026-03-26 | 5/250 |
|
|
|
[考研] 藥學(xué)105500求調(diào)劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
|
[考研] 一志愿 南京郵電大學(xué) 288分 材料考研 求調(diào)劑 +3 | jl0720 2026-03-26 | 3/150 |
|
|
|
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
|
[考研] 07化學(xué)303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|