计算化学公社

标题: VESTA导出POSCAR文件晶格常数为什么会变化 [打印本页]

作者
Author:
往事随风0809    时间: 2022-3-26 17:29
标题: VESTA导出POSCAR文件晶格常数为什么会变化
(, 下载次数 Times of downloads: 49) 图一:VESTA里Mg单胞晶格常数
(, 下载次数 Times of downloads: 24) 图二:VESTA导出POSCAR的晶格常数
求教为什么会变化,对计算有影响么

作者
Author:
ML-catalyze    时间: 2022-3-26 19:07
建议把INCAR贴出来看下
作者
Author:
往事随风0809    时间: 2022-3-26 19:35
ML-catalyze 发表于 2022-3-26 19:07
建议把INCAR贴出来看下

你好,这个是还没有算,只是从MS导出的Mg.cif用vesta打开是正常的晶格常数,再用VESTA导出成POSCAR的话晶格常数就变了,成了图2
作者
Author:
乐平    时间: 2022-3-26 19:46
本帖最后由 乐平 于 2022-3-26 13:49 编辑

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

第一个截图里给出的是晶胞参数的一种定义方式,给出了 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°  


作者
Author:
往事随风0809    时间: 2022-3-26 19:48
乐平 发表于 2022-3-26 19:46
你需要补充一下晶体结构的知识。

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

非常感谢
作者
Author:
wangyuxiao    时间: 2023-3-22 10:38
往事随风0809 发表于 2022-3-26 19:48
非常感谢

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




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3