 |
 |
 |
| 24小時(shí)熱門(mén)版塊排行榜 |
- >論壇更新日志 (701)
- >考研 (83)
- >蟲(chóng)友互識(shí) (43)
- >導(dǎo)師招生 (28)
- >考博 (13)
- >碩博家園 (12)
- >論文投稿 (12)
- >基金申請(qǐng) (11)
- >教師之家 (9)
- >無(wú)機(jī)/物化 (5)
- >文獻(xiàn)求助 (5)
- >公派出國(guó) (4)
- >休閑灌水 (3)
- >招聘信息布告欄 (2)
- >論文道賀祈福 (2)
- >文獻(xiàn)求助完結(jié)子版 (2)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 550 | 回復(fù): 2 | |||
qinform木蟲(chóng) (著名寫(xiě)手)
|
[求助]
求助 晶格優(yōu)化腳本 一個(gè)命令格式的修改
|
|
在vasp入門(mén)指南里 晶格優(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}‘ 不知道怎么修改。 請(qǐng)幫忙。 |
回復(fù)此樓» 猜你喜歡
309求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷
已經(jīng)有8人回復(fù)
-
求調(diào)劑
已經(jīng)有5人回復(fù)
-
求化學(xué)調(diào)劑
已經(jīng)有5人回復(fù)
-
085600，材料與化工321分求調(diào)劑
已經(jīng)有10人回復(fù)
-
0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分
已經(jīng)有5人回復(fù)
-
一志愿南開(kāi)大學(xué)0710生物學(xué)359求調(diào)劑
已經(jīng)有3人回復(fù)
-
085600，專(zhuān)業(yè)課化工原理，320分求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料與化工272求調(diào)劑
已經(jīng)有16人回復(fù)
-
290求調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 在comsol里面怎么將一個(gè)模塊的計(jì)算結(jié)果用到另一個(gè)應(yīng)力模塊作為加載源
已經(jīng)有9人回復(fù)
- 如何在MS的模型中的晶胞一個(gè)棱上（只有一個(gè)原子）摻雜替代這個(gè)原子
已經(jīng)有13人回復(fù)
- BASH腳本的一個(gè)小問(wèn)題
已經(jīng)有3人回復(fù)
- 關(guān)于優(yōu)化晶格常數(shù)
已經(jīng)有4人回復(fù)
- ENCUT測(cè)試腳本
已經(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文件寫(xiě)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入門(mén)的晶格優(yōu)化部分沒(méi)看懂，問(wèn)一下
已經(jīng)有25人回復(fù)
- 【求助】求證晶格常數(shù)的優(yōu)化
已經(jīng)有15人回復(fù)
- 【求助】關(guān)于晶胞中取代一個(gè)原子的優(yōu)化問(wèn)題
已經(jīng)有14人回復(fù)
- 【求助】幫忙編個(gè)用于VASP小程序
已經(jīng)有22人回復(fù)
- 【求助】如何固定某幾個(gè)原子不動(dòng)？
已經(jīng)有20人回復(fù)
jerkwin
專(zhuān)家顧問(wèn) (正式寫(xiě)手)
-
專(zhuān)家經(jīng)驗(yàn): +14 - 應(yīng)助: 454 (碩士)
- 金幣: 20699.1
- 散金: 148
- 紅花: 81
- 帖子: 813
- 在線(xiàn): 2648.3小時(shí)
- 蟲(chóng)號(hào): 1023452
- 注冊(cè): 2010-05-19
- 專(zhuān)業(yè): 理論和計(jì)算化學(xué)
- 管轄: 分子模擬
; printf "%12.6f \n", $5}' 
2樓2013-09-21 22:10:29
jerkwin
專(zhuān)家顧問(wèn) (正式寫(xiě)手)
-
專(zhuān)家經(jīng)驗(yàn): +14 - 應(yīng)助: 454 (碩士)
- 金幣: 20699.1
- 散金: 148
- 紅花: 81
- 帖子: 813
- 在線(xiàn): 2648.3小時(shí)
- 蟲(chóng)號(hào): 1023452
- 注冊(cè): 2010-05-19
- 專(zhuān)業(yè): 理論和計(jì)算化學(xué)
- 管轄: 分子模擬
3樓2013-09-22 00:56:16
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 qinform 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 各位老師好，我的一志愿為北京科技大學(xué)085601材料專(zhuān)碩 +9 | Koxui 2026-03-28 | 9/450 |
|
|---|---|---|---|---|
|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +14 | Mr. Z 2026-03-25 | 14/700 |
|
|
[考研] 298求調(diào)劑 +3 | 種圣賜 2026-03-29 | 3/150 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +6 | 崔wj 2026-03-26 | 6/300 |
|
|
|
[考研] 356求調(diào)劑 +3 | gysy?s?a 2026-03-28 | 3/150 |
|
|
|
[考研] 學(xué)碩274求調(diào)劑 +9 | Li李魚(yú) 2026-03-26 | 9/450 |
|
|
|
[考研] 316求調(diào)劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
|
[考研] 壓國(guó)家一區(qū)線(xiàn)，求導(dǎo)師收留，有恩必謝！ +7 | 迷人的哈哈 2026-03-28 | 7/350 |
|
|
|
[考研] 本科新能源科學(xué)與工程，一志愿華理能動(dòng)285求調(diào)劑 +3 | AZMK 2026-03-27 | 5/250 |
|
|
|
[考研] 070300求調(diào)劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|
|
[考研] 085701環(huán)境工程，267求調(diào)劑 +16 | minht 2026-03-26 | 16/800 |
|
|
|
[考研] 286求調(diào)劑 +12 | PolarBear11 2026-03-26 | 12/600 |
|
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +7 | hanamiko 2026-03-27 | 7/350 |
|
|
|
[考研] 085405 考的11408求各位老師帶走 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
|
[考研] 279 分 求調(diào)劑 +4 | 睡個(gè)好覺(jué)_16 2026-03-24 | 4/200 |
|
|
|
[考研] 343求調(diào)劑 +4 | 贈(zèng)我一本書(shū) 2026-03-23 | 4/200 |
|
|
|
[考研] 327求調(diào)劑 +7 | prayer13 2026-03-23 | 7/350 |
|
|
|
[考研] 一志愿河工大 081700 276求調(diào)劑 +4 | 地球繞著太陽(yáng)轉(zhuǎn) 2026-03-23 | 4/200 |
|
|
|
[考研] 一志愿 南京郵電大學(xué) 288分 材料考研 求調(diào)劑 +3 | jl0720 2026-03-26 | 3/150 |
|
|
|
[考研] 309求調(diào)劑 +4 | gajsj 2026-03-25 | 5/250 |
|