 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 551 | 回復(fù): 2 | |||
qinform木蟲 (著名寫手)
|
[求助]
求助 晶格優(yōu)化腳本 一個(gè)命令格式的修改
|
|
在vasp入門指南里 晶格優(yōu)化部分，有一個(gè)腳本 #!/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 我想修改這個(gè)腳本，我用的程序輸出文件名是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ù)此樓» 猜你喜歡
材料科學(xué)與工程求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿大連理工大學(xué)材料求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑
已經(jīng)有9人回復(fù)
-
材料求調(diào)劑 一志愿哈工大總分298分，前三科223分
已經(jīng)有8人回復(fù)
-
274求調(diào)劑
已經(jīng)有3人回復(fù)
-
合肥區(qū)域性重點(diǎn)一本招收調(diào)劑
已經(jīng)有4人回復(fù)
-
085600材料與化工調(diào)劑
已經(jīng)有10人回復(fù)
-
22408 359分調(diào)劑
已經(jīng)有7人回復(fù)
-
抱歉
已經(jīng)有4人回復(fù)
-
085602化工求調(diào)劑（331分）
已經(jīng)有8人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 在comsol里面怎么將一個(gè)模塊的計(jì)算結(jié)果用到另一個(gè)應(yīng)力模塊作為加載源
已經(jīng)有9人回復(fù)
- 如何在MS的模型中的晶胞一個(gè)棱上（只有一個(gè)原子）摻雜替代這個(gè)原子
已經(jīng)有13人回復(fù)
- BASH腳本的一個(gè)小問題
已經(jīng)有3人回復(fù)
- 關(guān)于優(yōu)化晶格常數(shù)
已經(jīng)有4人回復(fù)
- ENCUT測試腳本
已經(jīng)有17人回復(fù)
- 【求助】lammps最適合安裝在哪一個(gè)版本的linux系統(tǒng)中？
已經(jīng)有7人回復(fù)
- 【求助】vasp優(yōu)化時(shí)選擇了isif=3 怎么限制晶格的一個(gè)維度呢？
已經(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)于晶胞中取代一個(gè)原子的優(yōu)化問題
已經(jīng)有14人回復(fù)
- 【求助】幫忙編個(gè)用于VASP小程序
已經(jīng)有22人回復(fù)
- 【求助】如何固定某幾個(gè)原子不動(dòng)？
已經(jīng)有20人回復(fù)
jerkwin
專家顧問 (正式寫手)
-
專家經(jīng)驗(yàn): +14 - 應(yīng)助: 454 (碩士)
- 金幣: 20699.1
- 散金: 148
- 紅花: 81
- 帖子: 813
- 在線: 2648.3小時(shí)
- 蟲號(hào): 1023452
- 注冊: 2010-05-19
- 專業(yè): 理論和計(jì)算化學(xué)
- 管轄: 分子模擬
; printf "%12.6f \n", $5}' 
2樓2013-09-21 22:10:29
jerkwin
專家顧問 (正式寫手)
-
專家經(jīng)驗(yàn): +14 - 應(yīng)助: 454 (碩士)
- 金幣: 20699.1
- 散金: 148
- 紅花: 81
- 帖子: 813
- 在線: 2648.3小時(shí)
- 蟲號(hào): 1023452
- 注冊: 2010-05-19
- 專業(yè): 理論和計(jì)算化學(xué)
- 管轄: 分子模擬
3樓2013-09-22 00:56:16
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 qinform 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿華東師范大學(xué)有機(jī)化學(xué)專業(yè)，初試351分，復(fù)試被刷求調(diào)劑! +7 | 真名有冰 2026-03-29 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 303求調(diào)劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|
|
[考研] 一志愿：西北大學(xué)，英一數(shù)一408-284分求調(diào)劑 +5 | 12.27 2026-03-27 | 5/250 |
|
|
|
[考研] 311求調(diào)劑 +7 | 藍(lán)月亮亮 2026-03-30 | 7/350 |
|
|
|
[考研] 317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑 +12 | 芋泥小鈴鐺 2026-03-28 | 12/600 |
|
|
|
[考研] 298求調(diào)劑 +3 | 什么是胖頭魚 2026-03-30 | 5/250 |
|
|
|
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +13 | 丹青奶蓋 2026-03-26 | 15/750 |
|
|
|
[考研] 一志愿中南大學(xué)化學(xué)0703總分337求調(diào)劑 +6 | niko- 2026-03-27 | 6/300 |
|
|
|
[考研] 275求調(diào)劑 +15 | Micky11223 2026-03-25 | 20/1000 |
|
|
|
[考研] 調(diào)劑310 +12 | 溫柔的晚安 2026-03-25 | 13/650 |
|
|
|
[考研] 086000生物與醫(yī)藥調(diào)劑 +5 | Feisty。 2026-03-28 | 9/450 |
|
|
|
[考研] 298求調(diào)劑 +4 | 種圣賜 2026-03-28 | 4/200 |
|
|
|
[考研] 356求調(diào)劑 +3 | gysy?s?a 2026-03-28 | 3/150 |
|
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | 劉柯@ 2026-03-24 | 4/200 |
|
|
|
[考研] 314求調(diào)劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
|
[考研] 286求調(diào)劑 +4 | lim0922 2026-03-26 | 4/200 |
|
|
|
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4 | 小魚愛有機(jī) 2026-03-25 | 4/200 |
|