| 查看: 622 | 回复: 2 | |||
qinform木虫 (著名写手)
|
[求助]
求助 晶格优化脚本 一个命令格式的修改
|
|
在vasp入门指南里 晶格优化部分，有一个脚本 #!/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 }’‘ 不能正确输出总能的数值，输出为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}‘ 不知道怎么修改。 请帮忙。 |
回复此楼» 猜你喜欢
时间戳又变了8-15
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有50人回复
-
filecode
已经有11人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 在comsol里面怎么将一个模块的计算结果用到另一个应力模块作为加载源
已经有9人回复
- 如何在MS的模型中的晶胞一个棱上（只有一个原子）掺杂替代这个原子
已经有13人回复
- BASH脚本的一个小问题
已经有3人回复
- 关于优化晶格常数
已经有4人回复
- ENCUT测试脚本
已经有17人回复
- 【求助】lammps最适合安装在哪一个版本的linux系统中？
已经有7人回复
- 【求助】vasp优化时选择了isif=3 怎么限制晶格的一个维度呢？
已经有16人回复
- 【求助】VASP怎么做晶格参数的结构优化
已经有17人回复
- 【求助】如何根据cell文件写poscar
已经有11人回复
- 【求助】vasp优化晶格参数a
已经有13人回复
- 【讨论】lammps里的boundary的位置能够随意更改？
已经有14人回复
- 【求助】如何从Murnaghan方程的结果中得到优化的晶格参数？
已经有11人回复
- 【求助成功】晶格优化中，murn.f的使用方法
已经有7人回复
- 【求助】侯老师的vasp入门的晶格优化部分没看懂，问一下
已经有25人回复
- 【求助】求证晶格常数的优化
已经有15人回复
- 【求助】关于晶胞中取代一个原子的优化问题
已经有14人回复
- 【求助】帮忙编个用于VASP小程序
已经有22人回复
- 【求助】如何固定某几个原子不动？
已经有20人回复
jerkwin
专家顾问 (正式写手)
-
专家经验: +14 - 应助: 454 (硕士)
- 金币: 20659.1
- 散金: 148
- 红花: 81
- 帖子: 813
- 在线: 2648.7小时
- 虫号: 1023452
- 注册: 2010-05-19
- 专业: 理论和计算化学
- 管辖: 分子模拟
; printf "%12.6f \n", $5}' 
2楼2013-09-21 22:10:29
jerkwin
专家顾问 (正式写手)
-
专家经验: +14 - 应助: 454 (硕士)
- 金币: 20659.1
- 散金: 148
- 红花: 81
- 帖子: 813
- 在线: 2648.7小时
- 虫号: 1023452
- 注册: 2010-05-19
- 专业: 理论和计算化学
- 管辖: 分子模拟
3楼2013-09-22 00:56:16
