往事随风0809 发表于 2022-3-26 19:48
非常感谢

您好，这个截图是VESTA导出来的POSCAR，优化后的晶格常数是CONTCAR里面的吗

乐平 发表于 2022-3-26 19:46
你需要补充一下晶体结构的知识。

第一个截图里给出的是晶胞参数的一种定义方式，给出了 a, b, c, 以及它 ...

非常感谢

你需要补充一下晶体结构的知识。

第一个截图里给出的是晶胞参数的一种定义方式，给出了 a, b, c, 以及它们之间的夹角 alpha = 90°,  beta = 90°, gamma = 120°
请注意，alpha 是 b 和 c 之间的夹角， beta 是 a 和 c 之间的夹角，gamma 是 a 和 b 之间的夹角。


第二个截图给出的是另一种定义方式，是矩阵表达或者说是向量表达。

你可以把
3.2093999386      0.0000000000    0.0000000000   看成是 a
-1.6046999693     2.7794218777     0.0000000000   看成是 b
0.0000000000     0.0000000000     5.2104997635   看成是 c
加上坐标系
       x                         y                         z

可以看出， a 在 x 轴的正方向上（因为 y, z 分量都为 0），c 在 z 轴的正方向（因为，x, y 分量都为 0）
b 在 x, y 平面内（因为 z 分量为 0）

那么你可以算出  |b| = √(-1.6046999693^2 + 2.7794218777^2) = 3.209399938572445 = |a|
也就是说矩阵表达的 a 和 b 长度是相等的，也就是和第一个截图里是一致的。

在计算机里浮点数运算会有误差，不过这个误差可以忽略。

你还可以验证 a 和 b 之间的夹角是 120°  

ML-catalyze 发表于 2022-3-26 19:07
建议把INCAR贴出来看下

你好，这个是还没有算，只是从MS导出的Mg.cif用vesta打开是正常的晶格常数，再用VESTA导出成POSCAR的话晶格常数就变了，成了图2

建议把INCAR贴出来看下