| 查看: 2438 | 回复: 5 | ||
| 当前主题已经存档。 | ||
cheerian铜虫 (正式写手)
|
[交流]
【讨论】vasp的输出文件解读--OUTCAR
|
|
|
vasp计算后出来一堆的输出文件，各有其用，随着输入的调整，输出也变化很大，但是不管做什么样的计算，其中最重要、最全面的输出文件莫过于OUTCAR了，但是正是因为全面，所以有太多的疑问，详细探讨、解读OUTCAR对于虫子们尤其是我们这些新虫大有裨益，希望大虫们不吝你们的口舌，说出你们的见解，谢谢！ [ Last edited by wuchenwf on 2009-6-13 at 23:04 ] |
回复此楼
cheerian
铜虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 5758.2
- 红花: 1
- 帖子: 427
- 在线: 83.5小时
- 虫号: 476233
- 注册: 2007-12-10
- 专业: 凝聚态物性I:结构、力学和
★ ★ ★ ★
qasd(金币+4,VIP+0):xiexie 3-7 17:49
qasd(金币+4,VIP+0):xiexie 3-7 17:49
|
转贴一些先（已实践）： 查看所计算体系的体积，使用下面的命令 grep ‘volume’ OUTCAR 得到的结果如下 volume/ion in A,a.u. = 32.92 222.17 volume of cell : 65.84 第一行给出体系的体积分别以Å3/atom, a.u.3/atom为单位给出的。 第二行给出体系的体积是以Å3/unit cell为单位给出的。 查看所计算体系的总能，使用下面的命令 当ISMEAR = -5时，Free energy TOTEN是与energy without entropy是相等，则用 grep ‘TOTEN’ OUTCAR 得到结果如下 free energy TOTEN = -7.910804 eV 当ISMEAR等于其他的值时，Free energy TOTEN是与energy without entropy是不相等，则用 grep ‘entropy=’ OUTCAR 得到结果如下 energy without entropy= -7.910804 energy(sigma->0) = -7.910804 在计算体系的结合能时，体系的总能取为energy without entropy后面的值。 查看所计算体系的费米能级，使用下面的命令 grep 'Fermi' OUTCAR | tail -1 得到的结果为 BZINTS: Fermi energy: 6.171330; 20.000000 electrons 上一行中第一个数就是体系的费米能级，第二个数就是体系的总价电子数。 注释：对半导体的体系，VASP取价带顶作为费米能级。对呈现金属性的体系，费米能级就是该体系的真实（具有物理意义的）费米能级。 查看所计算体系的倒格子基矢 在采用vi对OUTCAR编辑时，用下面的命令来查找 g/reciprocal lattice vectors 或 g/recip 查看所计算体系中原子的受力情况 在采用vi对OUTCAR编辑时，用下面的命令来查找 g/TOTAL-FORCE 原子所受的力的单位是eV/angstrom。 |
2楼2009-03-07 16:56:56
cheerian
铜虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 5758.2
- 红花: 1
- 帖子: 427
- 在线: 83.5小时
- 虫号: 476233
- 注册: 2007-12-10
- 专业: 凝聚态物性I:结构、力学和
疑问一-----期待解答
★ ★
wuchenwf(金币+2,VIP+0):这块用文字说可得挺多，建议楼主看一下谢希德的固体能带理论上面有非常详细的说明 3-7 19:18
wuchenwf(金币+2,VIP+0):这块用文字说可得挺多，建议楼主看一下谢希德的固体能带理论上面有非常详细的说明 3-7 19:18
|
FREE ENERGIE OF THE ION-ELECTRON SYSTEM (eV) --------------------------------------------------- free energy TOTEN = -0.314277 eV energy without entropy= -0.005753 energy(sigma->0) = -0.160015 FORCE on cell =-STRESS in cart. coord. units (eV/reduce length): Direction X Y Z XY YZ ZX -------------------------------------------------------------------------------------- Alpha Z 0.27 0.27 0.27 Ewald -30.64 -30.64 -30.64 0.00 0.00 0.00 Hartree 93.89 93.89 93.89 0.00 0.00 0.00 E(xc) -27.94 -27.94 -27.94 0.00 0.00 0.00 Local -147.85 -147.85 -147.85 0.00 0.00 0.00 n-local -20.54 -20.54 -20.54 0.00 0.00 0.00 augment 5.55 5.55 5.55 0.00 0.00 0.00 Kinetic 126.49 126.49 126.49 0.00 0.00 0.00 ------------------------------------------------------------------------------------- Total -0.77 -0.77 -0.77 0.00 0.00 0.00 in kB -2.42 -2.42 -2.42 0.00 0.00 0.00 哪位能给我们解读一下这一段的每行文字，怎么来的？  物理上有什么意义？主要用途在哪些方面？如果能够说出哪些参数会对这当中某个量有什么影响，那就五体投地感激了 |
3楼2009-03-07 17:02:28
cheerian
铜虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 5758.2
- 红花: 1
- 帖子: 427
- 在线: 83.5小时
- 虫号: 476233
- 注册: 2007-12-10
- 专业: 凝聚态物性I:结构、力学和
4楼2009-03-07 19:08:01
freer007
木虫 (正式写手)
- 1ST强帖: 1
- 应助: 1 (幼儿园)
- 金币: 2421.7
- 红花: 1
- 帖子: 676
- 在线: 146小时
- 虫号: 707811
- 注册: 2009-02-24
- 性别: GG
- 专业: 半导体材料
5楼2009-09-03 19:10:06
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
小木虫(金币+0.5):给个红包，谢谢回帖交流
|
支持！ |
6楼2009-09-03 22:43:45
