VASP-ICONST文件用法

ICONST文件在液体的NPT模拟中有一定的作用，他可以保持盒子的形状，避免液体盒子被压扁。

ICONST基本结构

前面M行称为primitive parameters，最后一行称为complex parameters（但事实上不是这样的，不止最后一行是complex parameters）。
其中的一些设置参数意义：


需要注意的是：所有的complex parameters必须放在primitive parameters后面。
最常用的是保持盒子是正交的：比如我想保持我的盒子是正交的，且a1:a2:a3=1:1:2。为了固定盒子的边长比例以及角度，我们定义这样一个表达式：

c _i * a _i + c _j * a _j = 0 首先定义夹角不变： LA 1 2 0 x方向和y方向的夹角不变 LA 1 3 0 x方向和z方向的夹角不变 LA 2 3 0 y方向和z方向的夹角不变 再是定义盒子的边长：

LR 1 0
LR 2 0
LR 3 0

这前面几个定义的都是primitive parameters。接下来定义complex parameter。
前面有6行primitive，那么对应complex就应该有6列。

S 1 0 0 0 0 0 0
S 0 1 0 0 0 0 0
S 0 0 1 0 0 0 0
S 0 0 0 1 -1 0 0
S 0 0 0 2 0 -1 0
S 0 0 0 0 2 -1 0

前面三行对应的是固定角度。后面三行定义的是边长。我们为了固定a1:a2:a3=1:1:2，根据flag S 的定义，这行的系数要和primitive进行线性组合。primitive指的是a1, a2, a3。这个c_i就是该行的系数。
从第四行开始，c_x*a_x+c_y*a_y=0。又ax:ay = 1:1, 则cx : cy = 1 : -1。
第五行：c_x*a_x+c_z*a_z=0。又ax:az = 1:2, 则cx : cz = 2: -1。
第六行：c_y*a_y+c_z*a_z=0。又ay:az = 1:2, 则cy : cz = 2: -1。